According to the physical and orbital characteristics in Carme group, Ananke group, and Pasiphae group of Jupiter’s moons, the distributions of physical and orbital properties in these three groups are investigated by using one-sample Kolmogorov–Smirnov nonparametric test. Eight key characteristics of the moons are found to mainly obey the Birnbaum–Saunders distribution, logistic distribution, Weibull distribution, and t location-scale distribution. Furthermore, for the moons’ physical and orbital properties, the probability density curves of data distributions are generated; the differences of three groups are also demonstrated. Based on the inferred results, one can predict some physical or orbital features of moons with missing data or even new possible moons within a reasonable range. In order to better explain the feasibility of the theory, a specific example is illustrated. Therefore, it is helpful to predict some of the properties of Jupiter’s moons that have not yet been discovered with the obtained theoretical distribution inference.
National Natural Science Foundation of China116722591130218711571301Ministry of Land and Resources Research of China in the Public Interest201411007Top-notch Academic Programs Project of Jiangsu Higher Education InstitutionsPPZY2015B1091. Introduction
There are 69 (the number has been refreshed to 79 by a team from Carnegie Institution for Science in July 2018. https://sites.google.com/carnegiescience.edu/sheppard/moons/jupitermoons) confirmed moons of Jupiter, around 65 of which have been well investigated [1, 2]. Considering the formation of Jupiter’s moons is influenced by diverse factors, which results in their physical characteristics differing greatly [3], Jupiter’s moons are divided into two basic categories: regular and irregular. The regular satellites are so named because they have prograde and near-circular orbits of low inclination, and they are in turn split into two groups: Inner satellites and Galilean [4]. The irregular satellites are actually the objects whose orbits are far more distant and eccentric. They form families that share similar orbits (semi-major axis, inclination, and eccentricity) and composition. These families, which are considered to be part of collisions, arise when the larger parent bodies were shattered by impacts from asteroids captured by Jupiter’s gravitational field. That is to say, at the early time of moons’ formation of the Jupiter, mass of the original moon’s ring was still sufficient to absorb the asteroid’s power and put it into orbit. So, part of the irregular moons might be created by the captured asteroids and then collided with other moons [5, 6], thus forming the various groups we see today. The identification of satellite families is tentative (please see [7, 8] for more details), and these families bear the names of their largest members. The most detailed modelling of the collisional origin of the families was reported in [9, 10].
According to this identification scheme [11], 60 moons were classified into 8 different groups, including Small Inner Regulars and Rings, Galileans, Themisto group, Carpo group, Himalia group, Carme group, and Ananke group as well as Pasiphae group, in addition to 9 satellites that do not belong to any of previous groups. The detailed information about all the groups of Jupiter’s moons can be found in Appendix A.
In recent years, many scientists have paid considerable attention to astronomical observation, physical research, and deep space exploration of small bodies, including asteroids, comets, and satellites, and so on. For example, planetary scientist Carry collected mass and volume estimates of 17 near-Earth asteroids, 230 main-belt and Trojan asteroids, 12 comets, and 28 trans-Neptunian objects from the known literature [12]. The accuracy and biases affecting the methods used to estimate these quantities were discussed and best-estimates were strictly selected. For the asteroids in retrograde orbit, there are at least 50 known moons of Jupiter’s that are retrograde, some of which are thought to be asteroids or comets that originally formed near the gas giant and were captured when they got too close. Kankiewicz and Włodarczyk selected the 25 asteroids with the best-determined orbital elements and then estimated their dynamical lifetimes by using the latest observational data, including astrometry and physical properties [13]. However, few researchers have tried to extrapolate the distribution of Jupiter’s satellites through statistical methods as yet.
In this paper, distributions of physical and orbital properties for the moons of Jupiter will be conducted by using one-sample Kolmogorov-Smirnov (K-S) test and maximum likelihood estimation [14–16]. Based on the analysis of satellites’ data, it is found surprisingly that the physical and orbital characteristics obey some distribution, such as Birnbaum-Saunders distribution [17], logistic distribution, Weibull distribution, and t location-scale distribution. Furthermore, the probability density curves of the data distribution are generated, and the differences of physical and orbital characteristics in the three groups are presented. In addition, the results of theoretical inference results are then proved to be feasible through one concrete example. Therefore, the results may be helpful to astronomers to discover new moons of Jupiter in the future.
2. Method of Distribution Inference
In statistics, the K-S test, one type of nonparametric test, is used to determine whether a sample comes from a population with a specific distribution. The null hypothesis of one-sample K-S test is that the Cumulative Distribution Function (CDF) of the data follows the adopted CDF. For one-sample case, null distribution of statistic can be obtained from the null hypothesis that the sample is extracted from a reference distribution. The two-sided test for “unequal” CDF tests the null hypothesis against the alternative that the CDF of the data is different from the adopted CDF. The test statistic is the maximum absolute difference between the empirical CDF calculated by x and the hypothetical CDF:(1)Dn=supxFnx-Fx,where F(x) is a given CDF and(2)Fnx=1n∑i=1nI-∞,xXiis the empirical distribution function of the observations Xi. Here I-∞,x(Xi) is the indicator function with the following form(3)I-∞,xXi=1,Xi≤x,0,otherwise.
According to Glivenko–Cantelli theorem [18], if the sample comes from distribution F(x), then Dn will almost surely converge to zero when n→∞. Therefore, we only focus on three satellite groups that have more than ten satellites in their groups, respectively. Of all these groups, Themisto group and Carpo group only contain one satellite, respectively, and only 4 moons were found separately in Small Inner Regulars and Rings, as well as in Galileans group. Moreover, there are 5 moons in Himalia group. For these 5 groups, there is no sufficient data for distribution inference, so we will focus on the Carme group, the Ananke group, and the Pasiphae group.
In the following sections, in order to use the one-sample K-S test, three sets of observed data from Jupiter’s moons will be tested against some commonly used distributions in statistics. The list of these distributions is shown in Table 1.
A list of common distributions.
Name of Selected Distributions
Parameters
Meaning
Beta
a
first shape parameter
b
second shape parameter
Birnbaum-Saunders
β
scale parameter
γ
shape parameter
Gamma
a
shape parameter
b
scale parameter
Logistic
μ
mean
σ
scale parameter
Nakagami
μ
shape parameter
ω
scale parameter
Normal
μ
mean
σ
standard deviation
Rice
s
noncentrality parameter
σ
scale parameter
t location-scale
μ
location parameter
σ
scale parameter
υ
shape parameter
Weibull
A
scale parameter
B
shape parameter
For the 9 continuous distributions in Table 1, the one-sample K-S test will be used to select the distribution with the highest confidence level. In order to characterize these distributions with well-defined parameter values, maximum likelihood estimation is also used. In addition, the parameter values of these distributions can be calculated from the observed data. However, when confidence level (typically set to 0.05) decreases, the rejection domain of the test becomes smaller, so the observed values that initially fall into the rejection domain may eventually fall into the acceptance domain. This situation will bring some trouble in practical application. To this end, we adopt p value, which represents the obtained confidence level by using the one-sample K-S test. In addition, the use of p value not only avoids determining the level of significance in advance, but also makes it easy to draw conclusions about the test by comparing the p value and significance level of the test. If the p value is greater than 0.05, we declare that the null hypothesis can be accepted. Furthermore, if the p values of several distributions are all greater than 0.05, the distribution with the largest p value should be selected, and the corresponding distribution will be the most appropriate one to fit the observed data. Our results can be found in the tables in Appendix B.
3. Distribution Inference of Satellite Groups
In this section, distributions of several diverse physical and orbital properties for the Carme group, Ananke group, and Pasiphae group are inferred sequentially.
3.1. Carme Group
There are 15 moons in the Carme group (please see Appendix A.). Due to the lack of enough data from S/2010 J1, the number of adopted moons is 14, which means the length of each data set is 14 from the mathematical perspective. In addition, considering all the mean densities have been calculated at being 2.60 g/cm3, the surface gravity of Carme and other moons in this group is 0.017 m/s2 and 0.001 m/s2, respectively [19]. Therefore, these data obviously do not obey the distributions.
Based on the previous method of statistical distribution inference and MATLAB 2016a (Intel Core i5-3230 M, CPU 2.60 GHz), distribution inferences in the Carme group can be found in Table 2, and the last column represents the p value. The smaller the p value, the greater the significance because it tells us that hypothesis under consideration may not be sufficient to explain the observations. The hypothesis will be rejected if any of these probabilities is less than or equal to a small, fixed but arbitrarily predefined threshold value. The null hypothesis here refers to data obeying a particular distribution, and the alternative hypothesis assumes that the data does not obey the distribution.
The distribution inference in each physical and orbital characteristic (Carme group).
Characteristics
Distribution inference
Parameter estimates
p-value
Semi-major axis (107 km)
logistic distribution
μ = 2.3326, σ = 0.00651346
0.9988
Mean orbit velocity (103 km/h)
logistic distribution
μ = 8.21989, σ = 0.0112805
0.6607
Orbit eccentricity (10−1)
Birnbaum-Saunders distribution
β = 2.54254, γ = 0.0330888
0.6245
Inclination of orbit (102°)
t location-scale distribution
μ= 1.6511, σ = 0.0017, υ = 0.8751
0.6662
Equatorial radius (km)
t location-scale distribution
μ= 1.65708, σ = 0.440683, υ = 1.14501
0.7119
Escape velocity (km/h)
t location-scale distribution
μ = 7.32997, σ = 1.58419, υ= 1.06821
0.6619
From Table 2, semi-major axis and mean orbit velocity obey the logistic distribution, of which parameter μ denotes the average orbital eccentricity and parameter σ plays a key role in representing the variance of the data set, through the variance formula σ2π2/3, so the variances of these logistic distributions are 1.3957E-04 and 4.1863E-04, respectively. As the parameter μ increases, it indicates that the average semi-major axis and mean orbit velocity of the Carme group increase. As the parameter σ increases, the data in these two characteristics gradually disperse, and the discrepancy between the data and the average value increases. The meaning of decreasing the parameter μ and σ shares the same principle as that of increasing the parameters. The orbit eccentricity follows the Birnbaum–Saunders distribution, which is unimodal with a median of β. The mean value and the variance of the distribution can be calculated by the following relationships:(4)μ=β1+γ22,σ2=βγ21+5γ24.
Therefore, the mean value of orbital eccentricity can be calculated to be 0.2543 and the variance is 7.08E-9.
The inclination of orbit, equatorial radius, and escape velocity are subject to the t location-scale distribution, which contains the scale parameter σ, the location parameter μ, and the shape parameter υ. Without loss of generality, we assume that the data vector x obeys the t location-scale distribution, and then we have (x-μ)/σ~t(υ), which obeys Student’s t-distribution; here υ represents the degrees of freedom. As can also be seen in Table 2, the inclination of orbit obeys the t location-scale distribution with parameters (1.6511, 0.0017, 0.8751) and the mean inclination is 165.11°. The equatorial radius follows the t location-scale distribution with parameters (1.65708, 0.440683, 1.14501). Therefore, the moons in Carme group have an average equatorial radius of 1.65708 km. In addition, the escape velocity characteristic obeys the t location-scale distribution with parameters (7.32997, 1.58419, 1.06821) and the average escape velocity is 7.32997 km/h. As the parameter μ changes, the average equatorial radius and escape velocity of the group also change accordingly. In the t location-scale distribution, the variance is σ2υ/(υ-2), and when the shape parameter υ is greater than two, the variance of the distribution is defined. Therefore, the variances of these two specific t location-scale distributions cannot be defined.
3.2. Ananke Group
There are 10 moons in the Ananke group, of which we do not have enough data about S/2010 J2. Thus, the length of each characteristic of the remaining moons in this group will be 9.
From our discussion on the t location-scale distribution in Carme group, it is easy to understand the distribution inference of semi-major axis, mean orbit velocity, equatorial radius, surface gravity, and escape velocity in the Ananke group (see Table 3). However, compared with Table 2, it is noted that the orbit eccentricity and mass properties in the Ananke group are subject to the Weibull distribution. Parameters A and B represent the scale and shape parameters of the distribution [20], respectively, which together determine the mean and variance of the distribution.
The distribution inference in each orbital or physical characteristic (Ananke group).
Characteristics
Distribution inference
Parameter estimates
p-value
Semi-major axis (107 km)
t location-scale distribution
μ = 2.11665, σ = 0.00915557, υ = 1.52402
0.9304
Mean orbit velocity (103 km/h)
Weibull distribution
A = 2.29923, B = 21.9834
0.9788
Orbit eccentricity (10−1)
t location-scale distribution
μ = 8.72481, σ = 0.0203158, υ = 1.36328
0.6017
Inclination of orbit (102°)
t location-scale distribution
μ= 1.4879, σ = 0.0127, υ = 1.4538
0.9075
Equatorial radius (km)
t location-scale distribution
μ = 1.80724, σ = 0.660297, υ = 1.19615
0.7231
Mass (1013 kg)
Weibull distribution
A = 38.3573, B = 0.368827
0.4304
Surface gravity ((m/s2)/103)
t location-scale distribution
μ = 0.167044, σ = 0.0636087, υ = 1.40815
0.4210
Escape velocity (km/h)
t location-scale distribution
μ = 8.24903, σ = 2.61085, υ = 1.15837
0.7262
3.3. Pasiphae Group
There are 19 moons in the Pasiphae group; due to the lack of data about three moons, S/2011 J2, S/2017 J1, and S/2016 J1, the remaining 16 will be studied in this subsection.
The distributions in Table 4 are inferred to be similar to the distribution in Tables 2 and 3. Therefore, we can easily understand the parameters of these distributions.
The distribution inference in each orbital or physical characteristic (Pasiphae group).
Characteristics
Distribution inference
Parameter estimates
p-value
Semi-major axis (107 km)
t location-scale distribution
μ = 2.38241, σ = 0.0273627, υ = 0.730397
0.3730
Mean orbit velocity (103 km/h)
logistic distribution
μ = 8.22803, σ = 0.226046
0.4015
Orbit eccentricity (10−1)
logistic distribution
μ = 2.95251, σ = 0.599373
0.9550
Inclination of orbit (102°)
logistic distribution
μ = 1.5136, σ = 0.0339
0.8987
Equatorial radius (km)
Weibull distribution
A = 3.69781, B = 0.781648
0.0859
Mass (1013 kg)
Weibull distribution
A = 72.5121, B = 0.266585
0.1096
Surface gravity ((m/s2)/103)
Weibull distribution
A = 3.09876, B = 0.843794
0.0594
Escape velocity (km/h)
Weibull distribution
A = 17.0302, B = 0.807236
0.0722
4. Comparison of Data Properties
Table 5 is given according to the previous distribution inference.
Distribution inference summary.
Characteristics
Carme Group
Ananke Group
Pasiphae Group
Semi-major axis
logistic distribution
t location-scale distribution
t location-scale distribution
Mean orbit velocity
logistic distribution
t location-scale distribution
logistic distribution
Orbit eccentricity
Birnbaum–Saunders distribution
Weibull distribution
logistic distribution
Inclination of orbit
t location-scale distribution
t location-scale distribution
logistic distribution
Equatorial radius
t location-scale distribution
t location-scale distribution
Weibull distribution
Mass
None
Weibull distribution
Weibull distribution
Surface gravity
None
t location-scale distribution
Weibull distribution
Escape velocity
t location-scale distribution
t location-scale distribution
Weibull distribution
Based on the previous distribution inference, the properties of the moons’ data can be compared more specifically and conveniently. According to the distribution of the specific parameters, we get the following probability density function (PDF) diagram.
As shown in Figure 1, the semi-major axis in the Ananke group is the smallest, followed by the Carme group and the Pasiphae group. In addition, the PDF of the Pasiphae group is relatively flat, indicating a large dispersion in semi-major axis around Jupiter, while the data in the Carme and Ananke groups differ slightly. As can be seen from Figure 2, the trend of mean orbital velocity distribution in Carme group is the closest to each other, with the mean value being the smallest among the three groups, followed with the Pasiphae and Ananke groups. Obviously, the curve of the Pasiphae group is the flattest, which shows the mean orbit velocity in the Pasiphae group differing greatly. Figure 3 shows that the mean value in the Ananke group is the smallest, followed by the Carme and Pasiphae groups. It is obvious that the PDF curve of the Pasiphae group is flatter compared to the others, which means the orbital eccentricity of the Pasiphae group has a relatively large dispersion. In Figure 4, although the inclinations of orbit in the Ananke group and Pasiphae group are relatively close in value compared with the Carme group, and the data looks more dispersed than those in Carme group, the inclination of orbit in the Ananke group has the same distribution as the Carme group.
The PDF curves of semi-major axis around Jupiter.
The PDF curves of mean orbit velocity.
The PDF curves of orbit eccentricity.
The PDF curves of inclination.
Figure 5 shows the PDF curves of equatorial radii of these three groups. Data attributes are similar; most of the moons’ radiuses are less than 4 km. Figure 6 indicates similarities between the Ananke and Pasiphae groups. As can be seen from these curves, most of the moons in these two groups are of relatively small mass.
The PDF curves of equatorial radius.
The PDF curves of mass.
Figure 7 illustrates the surface gravity PDF curves in the Ananke and Pasiphae groups. The difference in density between the two groups indicates that the surface gravity of the Ananke moons is higher than that of the Pasiphae group. The PDF plots corresponding to the escape velocity are shown in Figure 8, where it is clear that the Pasiphae group is flatter than the other two groups. The Carme group also has similar escape velocity density to the Ananke group, but the former average escape velocity is smaller.
The PDF curves of surface gravity.
The PDF curves of escape velocity.
5. Verification of Rationality of Theoretical Results
In this section, we take the semi-major axis and mean orbit velocity of the moons in the Carme group as concrete examples to illustrate the rationality of the statistical inferences in the preceding sections. As can be seen in Table 2, semi-major axis (sma) and mean orbit velocity (mov) obey the logistic distribution with parameters (2.3326, 0.00651346) and (8.21989, 0.0112805), respectively. So the corresponding predicted PDFs can be written as (5)fpre,smav;μ,σ=e-v-2.3326/0.006513460.006513461+e-v-2.3326/0.006513462and(6)fpre,movv;μ,σ=e-v-8.21989/0.001128050.001128051+e-v-8.21989/0.001128052.
Note that the orbital and physical properties are not independent of each other. For instance, the semi-major axis r is related to the mean orbit velocity v through the relation v=GM/r, where GM is the mass parameter. Then the PDF of mean orbit velocity can also be derived analytically as follows(7)fana,movv;μ,σ=2GMv3fpre,smaGMv2;2.3326,0.00651346=2GMv3e-GM-2.3326v2/0.00651346v20.006513461+e-GM-2.3326v2/0.00651346v22Although the PDF of mean orbit velocity obtained by different methods has different mathematical representations, from Figure 9, we can find that fana,movv;μ,σ (PDF curve represented by red circles) obtained by the analytical method is in good agreement with fpre,movv;μ,σ (PDF curve represented by blue circles) obtained by statistical inference.
Comparison of statistical prediction distribution and analytical derivation distribution.
However, what we need to pay attention to here is that some physical features and orbital elements are mixed and can be linked by some mathematical formulas similar to the above. Theoretically, the distribution of another variable in the formula can be solved by a known defined distribution; although this may be a complex process because the probability density function may contain some transcendental functions and gamma functions, it can be achieved. Yet the known distribution becomes uncertain now; that is, the distribution exists with a certain probability. So, there will be a certain risk when we calculate the distribution of the linked variable based on this uncertain distribution, especially when the possibility of the inferred distribution is not very high.
In addition, in order to further show that the results of the KS test agree well with the actual observed results, we also compared the best-fit CDFs (the CDFs of inferred statistically) and the observed CDFs. From Figures 10 and 11, it can be seen that the best-fit CDF and the observed CDF of semi-major axis agree better than the case of mean orbital velocity. This should be due to the fact that the p value of the former is 0.9988, which is obviously larger than that of the latter 0.6607.
Comparison of the observed CDF and the best-fit CDF of semi-major axis.
Comparison of the observed CDF and best-fit CDF of mean orbit velocity.
6. Application of the Distribution Inference
As mentioned previously, identification of irregular satellite families is tentative. However, after the distribution was inferred, it is found that these features of moons in the three different groups obey some selected distributions. Furthermore, the obtained population distribution can also be used to predict the characteristic data. Considering the rationality of proving this prediction method, suppose some characteristic data of a moon is unknown in a given irregular moon group. One specific feature of the moon can be predicted by using other moons’ data and one-sample K-S method. Here is an example:
Assume that some characteristic data of the moon Erinome (S/2000 J4) in Carme group is poorly known. Now, we try to predict the equatorial radius of the Erinome. First, we use one-sample K-S method to find the most appropriate continuous distribution for these data. The inference results are as shown in Table 6.
Distribution inference of Carme group without Erinome (S/2000 J4).
p-value
h
Parameter Values
Beta
4.02E-05
1
a = 0.202265
b = 0.173517
Birnbaum-Saunders
0.0307
1
β = 2.37288
γ = 0.943409
Gamma
0.0308
1
a = 1.05985
b = 3.11366
Logistic
0.0463
1
μ = 1.95282
σ = 1.72324
Nakagami
0.0039
1
μ = 0.282272
ω = 43.5177
Normal
0.0037
1
μ = 3.3
σ = 5.94531
Rice
5.73E-08
1
s = 0.138623
σ = 4.66385
t location-scale
0.6806
0
μ = 1.67257
σ = 0.487941
υ = 1.15623
Weibull
0.0698
0
A = 3.06061
B = 0.897697
Note: h = 0 and h = 1 indicate acceptance of the null hypothesis and rejection of the null hypothesis, respectively.
From the p values displayed in Table 6, it becomes clear that the corresponding p value 0.6806 is the largest and h=0, so the best distribution for the remaining characteristic data in the Carme group is t location-scale distribution with the PDF(8)fx∣μ,σ,υ=Γυ+1/2σπυΓυ/2υ+x-μ/σ2υ-υ+1/2,-∞<x<+∞,where Γ(·) is the Gamma function.
The distribution implies that average equatorial radius of the moons in Carme group is 1.67257 km with confidence interval [1.23204,2.1131] (we believe that the confidence interval will be smaller and shorter with continuous progress of observation technology) under the level of significance 0.05. For the true value of Erinome’s equatorial radius being 1.6 km, it is easy to find that we can use the estimated value to predict the equatorial radius or as a reference to study other relevant physical and orbital characteristics.
7. Conclusions
By using the one-sample K-S nonparametric test method of statistical inference, the distribution laws of the physical and orbital properties of Jupiter’s moons are investigated statistically in this paper. The physical and orbital characteristics of Jupiter’s moons are found to obey the Birnbaum–Saunders distribution, the logistic distribution, the Weibull distribution, and the t location-scale distribution.
In addition, the probability density curves of the data distributions are also generated, and the differences in the physical and orbital characteristics of the three groups are explained in more detail.
Through a specific example, we find that some moons’ missing data can be inferred by using the aforementioned distributional model and probability density function. More importantly, with the help of the distribution, it can be even helpful to predict the physical or orbital features of the undiscovered moon.
If future observations will allow for the expansion of the number of Jupiter’s moons, we believe that the distribution laws will be slightly modified as potential newly discovered distribution functions fit the increased sample better, and the distributions will probably tend to be more uniform; i.e., some of the different properties follow the same distribution obviously, but these will not change dramatically over a long period of time.
AppendixA.
See Table 7.
Information of Jupiter’s moons (source: the Jupiter satellite and moon page. Carnegie Institution. http://home.dtm.ciw.edu/users/sheppard/satellites/jupsatdata.html).
Number
Name
Designation
The Year of Discovery
Small Inner Regulars and Rings
XVI
Metis
1979
XV
Adrastea
1979
V
Amalthea
1892
XIV
Thebe
1979
Galileans
I
Io
1610
II
Europa
1610
III
Ganymede
1610
IV
Callisto
1610
Themisto Prograde Irregular Group
XVIII
Themisto
S/2000 J1
2000
Himalia Prograde Irregular Group
XIII
Leda
1974
VI
Himalia
1904
X
Lysithea
1938
VII
Elara
1905
LIII
Dia
S/2000 J11
2000
Carpo Prograde Irregular Group
XLVI
Carpo
S/2003 J20
2003
Retrograde Irregular Groups
---
Ananke Retrograde Irregular Group
LII
S/2010 J2
2010
XXXIII
Euanthe
S/2001 J7
2001
XXIX
Thyone
S/2001 J2
2001
XL
Mneme
S/2003 J21
2003
XXII
Harpalyke
S/2000 J5
2000
XXX
Hermippe
S/2001 J3
2001
XXVII
Praxidike
S/2000 J7
2000
XLII
Thelxinoe
S/2003 J22
2003
LX
S/2003 J3
2003
XXIV
Iocaste
S/2000 J3
2000
XII
Ananke
1951
Carme Retrograde Irregular Group
XLIII
Arche
S/2002 J1
2002
XXXVIII
Pasithee
S/2001 J6
2001
L
Herse
S/2003 J17
2003
XXI
Chaldene
S/2000 J10
2000
XXXVII
Kale
S/2001 J8
2001
XXVI
Isonoe
S/2000 J6
2000
XXXI
Aitne
S/2001 J11
2001
XXV
Erinome
S/2000 J4
2000
LI
S/2010 J1
2010
XX
Taygete
S/2000 J9
2000
XI
Carme
1938
XXIII
Kalyke
S/2000 J2
2000
XLVII
Eukelade
S/2003 J1
2003
LVII
S/2003 J5
2003
XLIV
Kallichore
S/2003 J11
2003
Pasiphae or No Strong Clustering Retrograde Irregular Group
XXXIV
Euporie
S/2001 J10
2001
LV
S/2003 J18
2003
LIV
S/2016 J1
2016
XXXV
Orthosie
S/2001 J9
2001
XLV
Helike
S/2003 J6
2003
LVIII
S/2003 J15
2003
XXXII
Eurydome
S/2001 J4
2001
XXVIII
Autonoe
S/2001 J1
2001
LVI
S/2011 J2
2011
XXXVI
Sponde
S/2001 J5
2001
LIX
S/2017 J1
2017
VIII
Pasiphae
1908
XIX
Megaclite
S/2000 J8
2000
IX
Sinope
1914
XXXIX
Hegemone
S/2003 J8
2003
XLI
Aoede
S/2003 J7
2003
XVII
Callirrhoe
S/1999 J1
1999
XLVIII
Cyllene
S/2003 J13
2003
XLIX
Kore
S/2003 J14
2003
The new Jupiter satellites discovered but yet to be numbered or named.
All numbered or named satellites are in the table above
S/2003 J2
2003
S/2003 J4
2003
S/2003 J9
2003
S/2003 J10
2003
S/2003 J12
2003
S/2003 J16
2003
S/2003 J19
2003
S/2003 J23
2003
S/2011 J1
2011
B. Distribution Inference Results
See Tables 8, 9, and 10
Carme group.
Beta
Birnbaum-Saunders
Gamma
Logistic
Nakagami
Normal
Rice
t location-scale
Weibull
Semi-major axis (107 km)
p-value
3.55E-04
0.9898
null
0.9987
0.0136
0.9891
0.9882
null
0.7911
h
1
0
null
0
1
0
0
null
0
parameters
a = 27.2522b = 0.272087
β = 2.33268γ = 0.00500346
Matrix is close to singular.
μ = 2.3326σ = 0.00651346
μ = 9985.19ω = 5.44166
μ = 2.33271σ = 0.0121133
s = 2.33268σ = 0.0116728
The calculated Hessian is not positive definite.
A = 2.33856B = 201.843
Mean orbit velocity (103 km/h)
p-value
8.80E-03
0.5488
null
0.6607
7.2972E-5
0.5598
0.5461
The calculated Hessian is not positive definite.
0.2983
h
1
0
null
0
0
0
0
null
0
parameters
a = 68.2281b = 0.245874
β = 8.22086γ = 0.00232688
Matrix is close to singular.
μ = 8.21989σ = 0.0112805
μ = 46161.3ω = 67.5832
μ = 8.22088σ = 0.0198551
s = 8.22086σ = 0.0191329
MLE (Maximum Likelihood Estimation) did not converge.
A = 8.23058B = 452.354
Orbit eccentricity (10−1)
p-value
4.80E-03
0.6245
4.0418E-7
0.5695
0.6128
0.6118
0.6002
null
0.3726
h
1
0
0
0
0
0
0
null
0
parameters
a = 4.79853b = 0.247005
β = 2.54254γ = 0.0330888
a = 909.093b = 0.00279831
μ = 2.53584σ = 0.0506215
μ = 226.315ω = 6.47876
μ = 2.54393σ = 0.0880083
s = 2.54251σ = 0.0848305
The calculated Hessian is not positive definite.
A = 2.58664B = 30.5827
Inclination of orbit (102°)
p-value
0.0994
0.0984
0.0991
0.2697
0.0998
0.1005
0.1005
0.6662
0.2271
h
0
0
0
0
0
0
0
0
0
parameters
a = 32931.5b = 166729
β = 1.6494γ = 0.005040
a = 39419.4189b = 0.0000418
μ = 1.6502σ = 0.004132
μ = 9866.0175ω = 2.7205
μ = 1.6494σ = 0.008298
s = 1.6494σ = 0.008298
μ = 1.6511σ = 0.0017υ = 0.8751
A = 1.6533B = 212.8357
Equatorial radius (km)
p-value
1.78E-05
0.0254
0.0249
0.0374
0.0024
0.0024
2.01E-08
0.7119
0.0574
h
1
1
1
1
1
1
1
0
0
parameters
a = 0.206976b = 0.185628
β = 2.32292γ = 0.913758
a = 1.10481b = 2.87704
μ = 1.90719σ = 1.606
μ = 0.286619ω = 40.5921
μ = 3.17857σ = 5.73011
s = 0.116687σ = 4.50454
μ = 1.65708σ = 0.440683υ = 1.14501
A = 2.98203B = 0.911426
Mass (1013 kg)
p-value
6.83E-05
8.58E-05
9.31E-04
0.0028
6.02E-04
3.10E-04
6.68E-12
null
0.0478
h
1
1
1
1
1
1
1
null
1
parameters
a = 0.0841809b = 0.128804
β = 153.808γ = 6.88668
a = 0.160966b = 5891.52
μ = 151.943σ = 941.547
μ = 0.0699297ω = 1.44943e+7
μ = 948.336σ = 3522.9
s = 113.882σ = 2491.35
MLE did not converge.
A = 31.522B = 0.309664
Surface gravity ((m/s2)/102)
p-value
1.64E-05
0.0171
0.0232
0.0162
0.0023
0.0024
3.21E-08
null
0.0346
h
1
1
1
1
1
1
1
null
1
parameters
a = 0.211325b = 0.186944
β = 0.180961γ = 0.87834
a = 1.16534b = 0.208399
μ = 0.149513σ = 0.120199
μ = 0.298776ω = 0.224286
μ = 0.242857σ = 0.421927
s = 0.001165σ = 0.334886
The calculated Hessian is not positive definite.
A = 0.232103B = 0.933741
Escape velocity (km/h)
p-value
1.40E-05
0.0126
0.0141
0.0287
0.0015
0.0018
1.37E-08
0.6619
0.0358
h
1
1
1
1
1
1
1
0
1
parameters
a = 0.210438b = 0.186748
β = 10.5296γ = 0.871682
a = 1.16971b = 12.0298
μ = 8.93572σ = 7.99281
μ = 0.296761ω = 769.929
μ = 14.0714σ = 24.8177
s = 0.703387σ = 19.6157
μ = 7.32997σ = 1.58419υ = 1.06821
A = 13.4252B = 0.932058
Ananke group.
Beta
Birnbaum-Saunders
Gamma
Logistic
Nakagami
Normal
Rice
tlocation-scale
Weibull
Semi-major axis (107 km)
p-value
0.0328
0.1347
0.1371
0.6785
0.1395
0.1364
0.1421
0.9304
0.684
h
1
0
0
0
0
0
0
0
0
parameters
a = 17.098b = 0.178068
β = 2.1074γ = 0.0142595
Matrix is close to singular.
μ = 2.1139σ = 0.0123376
μ = 2.10761ω = 0.0311751
μ = 2.10761σ = 0.0311751
s = 2.1074σ = 0.0295767
μ = 2.11665σ = 0.00915557υ = 1.52402
A = 2.11821B = 141.939
Mean orbit velocity (103 km/h)
p-value
4.30E-3
0.91
0.9211
0.9652
0.9312
0.9485
0.9408
0.7032
0.9788
h
1
0
0
0
0
0
0
0
0
parameters
a = 2.619b = 0.185662
β = 2.24183γ = 0.052298
a = 368.072b = 0.006.09908
μ = 2.24867σ = 0.0683413
μ = 92.6977ω = 5.05313
μ = 2.2449σ = 0.1227135
s = 2.24187σ = 0.116491
MLE did not converge.
A = 2.29923B = 21.9834
Orbit eccentricity (10−1)
p-value
8.99E-05
0.1634
0.1616
0.6557
0.1598
0.151
0.1581
0.6017
0.0821
h
1
0
0
0
0
0
0
0
0
parameters
a = 8.26353b = 0.198125
β = 8.75391γ = 0.0087315
Matrix is close to singular.
μ = 8.73706σ = 0.0313777
μ = 3235.86ω = 76.6427
μ = 8.75424σ = 0.0813835
s = 8.7539σ = 0.0772086
μ = 8.72481σ = 0.0203158υ = 1.36328
A = 8.79917B = 82.6918
Inclination of orbit (102°)
p-value
0.5027
0.4904
0.4991
0.7114
0.5077
0.5167
0.5167
0.9075
0.6297
h
0
0
0
0
0
0
0
0
0
parameters
a = 1951.25b = 11193.4
β = 1.4841γ = 0.0210
a = 2286.47b = 0.000649
μ = 1.4863σ = 0.0154
μ = 575.3654ω = 2.2045
μ = 1.4844σ = 0.0309
s = 1.4841σ = 0.0308
μ = 1.4879σ = 0.0127υ = 1.4538
A = 1.4986B = 52.7307
Equatorial radius (km)
p-value
7.43E-04
0.5045
0.4118
0.2689
0.158
0.1036
2.40E-03
0.7231
0.4889
h
1
0
0
0
0
0
1
0
0
parameters
a = 0.225054b = 0.144156
β = 2.32358γ = 0.827948
a = 1.45368b = 2.11188
μ = 2.18108σ = 1.39766
μ = 0.395632ω = 23.241
μ = 3.07σ = 3.91806
s = 0.126333σ = 3.40798
μ = 1.80724σ = 0.660297υ = 1.19615
A = 3.18574B = 1.07908
Mass (1013 kg)
p-value
0.0025
0.0341
0.0668
0.03
0.037
6.00E-03
2.74E-08
null
0.4304
h
1
1
0
1
1
1
1
null
0
parameters
a = 0.0964716b = 0.106068
β = 60.0159γ = 4.33548
a = 0.219599b = 1411.32
μ = 78.0193σ = 298.807
μ = 0.0927262ω = 898663
μ = 309.925σ = 944.346
s = 22.4209σ = 670.178
Maximum likelihood estimation did not converge.
A = 38.3573B = 0.368827
Surface gravity ((m/s2)/102)
p-value
7.52E-04
0.2802
0.2227
0.3042
0.0982
0.0933
3.50E-03
0.421
0.2732
h
1
0
0
0
0
0
1
0
0
parameters
a = 0.246544b = 0.14788
β = 0.198988γ = 0.732471
a = 1.82168b = 0.137236
μ = 0.190652σ = 0.101607
μ = 0.482908ω = 0.129
μ = 0.25σ = 0.271825
s = 0.0041522σ = 0.253955
μ = 0.167044σ = 0.0636087υ = 1.40815
A = 0.270122B = 1.21203
Escape velocity (km/h)
p-value
6.67E-04
0.3295
0.2788
0.2429
0.1138
0.0989
1.70E-03
0.7262
0.3567
h
1
0
0
0
0
0
1
0
0
parameters
a = 0.230441b = 0.145225
β = 10.6647γ = 0.778473
a = 1.58042b = 8.66859
μ = 9.82501σ = 5.97056
μ = 0.418467ω = 445.3
μ = 13.7σ = 16.9184
s = 0.531109σ = 14.9178
μ = 8.24903σ = 2.61085υ = 1.15837
A = 14.419B = 1.11631
Pasiphae group.
Beta
Birnbaum-Saunders
Gamma
Logistic
Nakagami
Normal
Rice
t location-scale
Weibull
Semi-major axis (107 km)
p-value
0.0012
0.1039
0.1069
0.2952
0.1092
0.1238
0.1125
0.373
0.1848
h
1
0
0
0
0
0
0
0
0
parameters
a = 4.22446b = 0.262587
β = 2.28804γ = 0.068061
a = 222.433b = 0.0103103
μ = 2.32σ = 0.0823946
μ = 57.3364ω = 5.28181
μ = 2.29334σ = 0.154515
s = 2.28844σ = 0.14977
μ = 2.38241σ = 0.0273627υ = 0.730397
A = 2.35522B = 23.6908
Mean orbit velocity (103 km/h)
p-value
9.38E-04
0.2085
0.2025
0.4015
0.1961
0.2048
0.1906
null
0.1794
h
1
0
0
0
0
0
0
null
0
parameters
a = 3.19406b = 0.286862
β = 8.27432γ = 0.0465993
a = 455.023b = 0.0182041
μ = 8.22803σ = 0.226046
μ = 112.484ω = 68.7679
μ = 8.28331σ = 0.406311
s = 8.27394σ = 0.393632
The calculated Hessian is not positive definite.
A = 8.48321B = 19.8563
Orbit eccentricity (10−1)
p-value
0.0132
0.3639
0.5771
0.955
0.7272
0.9138
0.8914
null
0.9072
h
1
0
0
0
0
0
0
null
0
parameters
a = 0.766602b = 0.237408
β = 2.66229γ = 0.426688
a = 6.69584b = 0.434085
μ = 2.95251σ = 0.599373
μ = 2.02049ω = 9.47173
μ = 2.90656σ = 1.04493
s = 2.68061σ = 1.06913
MLE did not converge.
A = 3.24837B = 3.30468
Inclination of orbit (102°)
p-value
0.5345
0.5633
0.5423
0.8987
0.5206
0.5002
0.5002
0.5004
0.2133
h
0
0
0
0
0
0
0
0
0
parameters
a = 554.177b = 3091.06
β = 1.5191γ = 0.0389
a = 655.7110b = 0.0023
μ = 1.5136σ = 0.0339
μ = 162.6092ω = 2.3148
μ = 1.5203σ = 0.0599
s = 1.5191σ = 0.0599
μ = 1.5203σ = 0.0599υ = 5657176.32
A = 1.5508B = 23.9875
Equatorial radius (km)
p-value
1.85E-04
0.0224
0.0377
0.0153
0.0139
0.0079
2.09E-08
null
0.0859
h
1
1
1
1
1
1
1
null
0
parameters
a = 0.203821b = 0.201972
β = 2.80607γ = 1.20376
a = 0.770973b = 5.83678
μ = 2.50052σ = 2.9666
μ = 0.242429ω = 81.8925
μ = 4.5σ = 8.10876
s = 0.192009σ = 6.39782
The calculated Hessian is not positive definite.
A = 3.69781B = 0.781648
Mass (1013 kg)
p-value
0.0012
0.0051
0.0111
0.0047
0.0091
8.16E-04
4.36E-11
null
0.1096
h
1
1
1
1
1
1
1
null
0
parameters
a = 0.0820299b = 0.133182
β = 185.784γ = 8.45355
a = 0.147015b = 17059.2
μ = 722.941σ = 2437.85
μ = 0.0664969ω = 6.36374E+7
μ = 2507.97σ = 7838.61
s = 215.387σ = 5639.23
The calculated Hessian is not positive definite.
A = 72.5121B = 0.266585
Surface gravity ((m/s2)/102)
p-value
9.34E-5
0.0186
0.0238
0.0139
0.0062
0.0055
2.15E-08
null
0.0594
h
1
1
1
1
1
1
1
null
0
parameters
a = 0.216072b = 0.206344
β = 0.232932γ = 1.08145
a = 0.896179b = 0.390547
μ = 0.205059σ = 0.214293
μ = 0.270001ω = 0.44625
μ = 0.35σ = 0.587651
s = 0.0105316σ = 0.472316
The calculated Hessian is not positive definite.
A = 0.309876B = 0.843794
Escape velocity (km/h)
p-value
1.55E-04
0.0157
0.0296
0.0142
0.011
0.0074
2.39E-08
null
0.0722
h
1
1
1
1
1
1
1
null
0
parameters
a = 0.208766b = 0.203804
β = 12.9337γ = 1.14416
a = 0.822159b = 24.4022
μ = 11.3591σ = 12.8644
μ = 0.253336ω = 1566.94
μ = 20.0625σ = 35.2429
s = 0.946989σ = 27.9843
The calculated Hessian is not positive definite.
A = 17.0302B = 0.807236
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors acknowledge the support of National Natural Science Foundation of China (NSFC) through grant Nos. 11672259, 11302187, and 11571301; the Ministry of Land and Resources Research of China in the Public Interest through grant No. 201411007; and Top-notch Academic Programs Project of Jiangsu Higher Education Institutions (TAPP) through grant No. PPZY2015B109.
SheppardS. S.2017Carnegie InstitutionJupiterD. B.moons2017https://solarsystem.nasa.gov/planets/jupiter/moonsJewittD.HaghighipourN.Irregular satellites of the planets: Products of capture in the early solar system2007452612952-s2.0-3564896519110.1146/annurev.astro.44.051905.092459SheppardS. S.Outer irregular satellites of the planets and their relationship with asteroids, comets and Kuiper Belt objects200512293193342-s2.0-6354914959810.1017/S1743921305006824CarrollB. W.OstlieD. A.20062ndPearson10.1063/1.881827Elkins-Tanton LTJupiter and Saturn2010JewittD.SheppardS.PorcoC.BagenalF.DowlingT. E.McKinnonW. B.Jupiter's outer satellites and Trojans20041Cambridge, UKCambridge planetary science263280WikipediansMoons of Jupiter2018PediaPress GmbH7894NesvornýD.AlvarellosJ. L. A.DonesL.LevisonH. F.Orbital and collisional evolution of the irregular satellites200312613984292-s2.0-003858484610.1086/375461NesvornýD.BeaugéC.DonesL.Collisional origin of families of irregular satellites20041273176817832-s2.0-214269953510.1086/382099SheppardS. S.JewittD. C.An abundant population of small irregular satellites around Jupiter200342369372612632-s2.0-034592745510.1038/nature0158412748634CarryB.Density of asteroids2012731981182-s2.0-8487061807810.1016/j.pss.2012.03.009KankiewiczP.WłodarczykI.Dynamical lifetimes of asteroids in retrograde orbits201746844143415010.1093/mnras/stx558HoggR. V.CraigA. T.2004Pearson Education Asia Limited and Higher Education PressMR0137186XieH.CuiX.WanB.ZhangJ.Statistical analysis of radio interference of 1000 kV UHV AC double-circuit transmission lines in foul weather201622475510.17775/CSEEJPES.2016.00021WilcoxR. R.Some practical reasons for reconsidering the Kolmogorov-Smirnov test19975019202-s2.0-004003558010.1111/j.2044-8317.1997.tb01098.xZbl0886.62054BirnbaumZ. W.SaundersS. C.A new family of life distributions19696231932710.2307/3212003MR0253493Zbl0209.49801AthreyaK. B.LahiriS. N.2006New York, NY, USASpringerSpringer Texts in StatisticsMR2247694List of moons in orbit around Jupiter, 2018, https://www.universeguide.com/planet/jupiter#themoonsJiangR.MurthyD. N. P.A study of Weibull shape parameter: Properties and significance20119612161916262-s2.0-8005317514510.1016/j.ress.2011.09.003