首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
We address impact cratering on Io and Europa, with the emphasis on the origin of small craters on Europa as secondary to the primary impacts of comets on Io, Europa, and Ganymede. In passing we also address the origin of secondary craters generated by Zunil, a well-studied impact crater on Mars that is a plausible analog to impact craters on Io. At nominal impact rates, and taking volcanic resurfacing into account, we find that there should be 1.3 impact craters on Io, equally likely to be of any diameter between 100 m and 20 km. The corresponding model age of Europa's surface is between 60 and 100 Ma. This range of ages does not include a factor three uncertainty stemming from the uncertain sizes and numbers of comets. The mass of basaltic impact ejecta from Io to reach Europa is found to meet or exceed the micrometeoroid flux as a source of rock-forming elements to Europa's ice crust. To describe impact ejecta in more detail we adapt models for impact-generated spalls and Grady-Kipp fragments originally developed by Melosh. Our model successfully reproduces the observed size-number distributions of small craters on both Mars and Europa. However, the model predicts that a significant fraction of the 200-500 m diameter craters on Europa are not traditional secondary craters but are instead sesquinary craters caused by impact ejecta from Io that had gone into orbit about Jupiter. This prediction is not supported by observation, which implies that high speed spalls usually break up into smaller fragments that make smaller sesquinary craters. Iogenic basalts are also interesting because they provide stratigraphic horizons on Europa that in principle could be used to track historic motions of the ice, and they provide materials suitable to radiometric dating of Europa's surface.  相似文献   

2.
We present a total of 289 new astrometric observations of the inner jovian satellites, Amalthea and Thebe, obtained using the Cassini ISS narrow angle camera. Observations were made using image sequences from 2000 December 11-12 (inbound) and 2001 January 15-16 (outbound), at phase angles of approximately 2° and 122°, respectively. Target distances were of order 284 RJ, giving a maximum resolution of approximately 100 km/pixel. Centroided line and sample values for 239 observations of Amalthea and 50 of Thebe are provided, together with estimated camera pointing information for each image. Orbit fitting using a uniformly precessing Keplerian ellipse model, taking into account the oblateness of Jupiter up to terms in J6, gave RMS fit residuals of 0.364 and 0.443 pixel for Amalthea and Thebe, respectively (equivalent to 0.450 and 0.547 arcsec). RMS residuals relative to the JPL JUP230 ephemeris were 0.306 and 0.604 pixel (equivalent to 0.378 and 0.746 arcsec), for Amalthea and Thebe. The fitted orbital parameters confirm the relatively high inclinations of these satellites (0.374°±0.002° and 1.076°±0.003°, respectively), equivalent to maximum vertical displacements above Jupiter's equatorial plane of 1188±6 and 4240±12 km, respectively, consistent with current estimates of the half-thicknesses of the Amalthea and Thebe gossamer rings [Ockert-Bell, M.E., Burns, J.A., Dauber, I.J., Thomas, P.C., Veverka, J., Belton, M.J.S., Klaasen, K.P., 1999. Icarus 138, 188-213].  相似文献   

3.
We examine the effects of Io ejecta on the surface and environment of Europa. We find that the observed sulfur on the trailing side of Europa, when interpreted as a deposit in equilibrium between implanation of, and sputtering by, corotating Io ejecta, implies a slow loss of material from Europa by sputtering. From this we infer that the spectrum of particles sputtered from water ice is soft. The quantity of observed sulfur and its confinement to the trailing hemisphere appear to exclude significant implantation and sputtering by energetic heavy ions. We also conclude that the contribution from Europa to the magnetospheric plasma (even at Europa itself) is negligible compared to the matter ejected from Io.  相似文献   

4.
《Planetary and Space Science》2006,54(9-10):1014-1023
Faint rings of micrometre-sized dust particles embrace many planets in the Solar system. As a rule, they are replenished by ejecta from embedded atmosphereless moons. On a number of occasions, the ejecta are generated by hypervelocity meteoroid impacts into the moons. Small ejecta fragments are then swiftly shifted into rings by an array of non-gravitational forces, e.g. radiation pressure or plasma drag. A significant fraction of ejecta mass, however, is contained in relatively big, multi-micrometre fragments which are subject to gravity only. Having escaped from the satellite, they stay close to its orbit and form a belt around planet. This belt is itself a source of ring dust through collisional disruption of its particles. Here the contributions of belts to the respective rings are estimated for selected satellites of Jupiter and Saturn. The belts under review could supply substantially more dust to rings than the direct ejecta from satellites and should be taken into account when estimating ring dust budgets. The belts are very difficult to observe, however, and some of them remain a theoretical proposition. We find an appealing evidence for the belts due to Amalthea and Thebe around Jupiter, and for the belt due to Enceladus around Saturn.  相似文献   

5.
I. Kulyk  K. Jockers 《Icarus》2004,170(1):24-34
We present the results of photometric measurements of the inner jovian satellites Thebe, Amalthea and Metis based on extensive optical observations taken from October 1999 to January 2002. The observations were made in the phase angle range from 8.1° to 0.3°. The Two-Channel Focal Reducer of the Max-Planck Institute for Aeronomy attached to the 2-m RCC telescope at Terskol Observatory (Pik Terskol, Northern Caucasus) was used in coronagraph mode. The observations were performed at a wavelength of 0.887 μm. Mean observational uncertainties corresponding to 1σ rms errors were 3% for the leading and trailing sides of Amalthea, 7 and 9% for the leading and trailing sides of Thebe and 9% for the leading side of Metis after taking into account the longitude brightness variations. Photometric data calibrated on an absolute scale were used to evaluate the near-opposition behavior of satellite brightness. All three satellites exhibit significant opposition brightening, but the strength of this effect, measured as the ratios of intensities at α1=1.6° and α2=6.7° does not vary significantly among these satellites. In order to measure the opposition surge parameters the empirical law proposed by Karkoschka and Hapke's model were used. The parameters of the satellite opposition effects are presented and discussed. The values of geometric albedos calculated with best-fit Hapke parameters are 0.096, 0.157, and 0.24 for Thebe, Amalthea, and Metis respectively. We found that the average leading/trailing ratios of surface reflectance at the measured phase angles are 1.53±0.05, 1.25±0.04, 1.04±0.08 for Amalthea, Thebe, and Metis.  相似文献   

6.
We use the specific scintillations of jovian decametric radio sources (modulation lanes), which are produced by plasma inhomogeneities in the vicinity of that planet, to probe the inner magnetosphere of Jupiter. The positions and frequency drift of 1762 lanes have been measured on the DAM spectra from archives. A special 3D algorithm is used for space localization of field-aligned magnetospheric inhomogeneities by the frequency drift of modulation lanes. As a result, the main regions of the lane formation are found: the Io plasma torus; the magnetic shell of the Gossamer Ring at Thebe and Amalthea orbits; and the region above the magnetic anomaly in the northern magnetosphere. It is shown that modulation lanes reveal the depleted magnetic tubes in practically unvisited, innermost regions of the jovian magnetosphere. The local and probably temporal plasma enhancement is found at the magnetic shell of Thebe satellite. Hence, the modulation lanes are a valuable instrument for remote sensing of those parts of jovian magnetosphere, which are not studied yet in situ.  相似文献   

7.
Tenuous dust clouds of Jupiter's Galilean moons Io, Europa, Ganymede and Callisto have been detected with the in-situ dust detector on board the Galileo spacecraft. The majority of the dust particles have been sensed at altitudes below five radii of these lunar-sized satellites. We identify the particles in the duut clouds surrounding the moons by their impact direction, impact velocity, and mass distribution. Average particle sizes are between 0.5 and 1 μm, just above the detector threshold, indicating a size distribution with decreasing numbers towards bigger particles. Our results imply that the particles have been kicked up by hypervelocity impacts of micrometeoroids onto the satellites' surfaces. The measured radial dust density profiles are consistent with predictions by dynamical modeling for satellite ejecta produced by interplanetary impactors (Krivov et al., 2003, Planet. Space Sci. 51, 251-269), assuming yield, mass and velocity distributions of the ejecta from laboratory measurements. A comparison of all four Galilean moons (data for Ganymede published earlier; Krüger et al., 2000, Planet. Space Sci. 48, 1457-1471) shows that the dust clouds of the three outer Galilean moons have very similar properties and are in good agreement with the model predictions for solid ice-silicate surfaces. The dust density in the vicinity of Io, however, is more than an order of magnitude lower than expected from theory. This may be due to a softer, fluffier surface of Io (volcanic deposits) as compared to the other moons. The log-log slope of the dust number density in the clouds vs. distance from the satellite center ranges between −1.6 and −2.8. Appreciable variations of number densities obtained from individual flybys with varying geometry, especially at Callisto, are found. These might be indicative of leading-trailing asymmetries of the clouds due to the motion of the moons with respect to the field of impactors.  相似文献   

8.
In a previous paper (The Rotation of Europa, Henrard, Celest. Mech. Dyn. Astr., 91, 131–149, 2005) we have developed a semi-analytical theory of Europa, one of the Galilean satellites of Jupiter. It is based on a synthetic theory of the orbit of Europa and is developed in the framework of Hamiltonian formalism. It was assumed that Europa is a rigid body and Jupiter a point mass. Several additional effects should be investigated in order to complete the theory. The present contribution considers the effect of the shape of Jupiter and of the gravitational pull of Io. The sensitivity of the main theory to a change in the values of the moments of inertia of Europa is also considered.  相似文献   

9.
The distribution of ejecta from impact craters significantly affects the surface characters of satellites and asteroids. In order to understand better the distinctive features seen on Phobos, Deimos, and Amalthea, we study the dynamics of nearby debris but include several factors — planetary tides plus satellite rotation and nonspherical shape-that complicate the problem. We have taken several different approaches to investigate the behavior of ejecta from satellites near planets. For example, we have calculated numerically the usual pseudoenergy (Jacobi) integral. This is done in the framework of a restricted three-body problem where we model the satellites as triaxial ellipsoids rather than point masses as in past work. Iso-contours of this integral show that Deimos and Amalthea are entirely enclosed by their Roche lobes, and the surfaces of their model ellipsoids lie nearly along equipotentials. Presumably this was once also the case for Phobos, before tidal evolution brought it so close to Mars. Presently the surface of Phobos overflows its Roche lobe, except for the regions within a few kilometers of the sub- and anti-Mars points. Thus most surface material on Phobos is not energetically bound: nevertheless it is retained by the satellite because local gravity has an inward component everywhere. Similar situations probably prevail for the newly discovered satellite of Jupiter (J14) and for the several objects found just outside Saturn's rings. We have also examined the fate of crater ejecta from the satellites of Mars by numerical integration of trajectories for particles leaving their surfaces in the equatorial plane. The ejecta behavior depends dramatically on the longitude of the primary impact, as well as on the speed and direction of ejection. Material thrown farther than a few degrees of longitude remains in flight for an appreciable time. Over intervals of an hour or more, the satellites travel through substantial arcs of their orbits, so that the Coriolis effect then becomes important. For this reason the limit of debris deposition is elongated toward the west while debris thrown to the east escapes at lower ejection velocities. We display some typical trajectories, which include many interesting special effects, such as loops, cusps, “folded” ejecta blankets, and even a temporary satellite of Deimos. Besides being important for understanding the formation of surface features on satellites, our work is perhaps pertinent to regolith development on small satellites and asteroids, and also to the budgets of dust belts around planets.  相似文献   

10.
W.-H. Ip 《Icarus》1996,120(2):317-325
The newly detected oxygen atmosphere of Europa is modeled by invoking charged particle sputtering with H2O and O2molecules as the main ejecta. The magnetospheric corotating ions could provide the required source strength (∼3 × 1026sec−1) of O2molecules if a fraction (∼20%) of the exospheric ions were recycled to Europa's surface where they produce additional sputtering product. Two exospheric components are expected to form: an extended corona with a size of a few satellite radii which is composed of sputtered molecules in ballistic motion, and a thermal population with a surface density of 108–109cm−3and a scale height of about 20 km. The electron impact ionization of this exosphere would lead to an Io-like interaction with the jovian magnetosphere with a field-aligned Birkeland current of about 5 × 105A.  相似文献   

11.
Observations of an occultation of Europa by Io are fitted by a model light curve. The model has five free parameters, namely the radius of Europa, the impact parameter, the brightness ratio of the satellites, the time of midevent and the mean relative velocity. The model assumes a fixed value for the radius of Io and for the solar phase angle α, and that Europa has a uniform surface brightness. The OC residuals of the best fitting light curve are very small (~0.002 mag) and of a purely random nature; there is no evidence of albedo features. Taking α = 0 does not affect significantly the quality of the fit. Six mutual eclipses were also observed, and their times of minima agree well with the predictions of Aksnes Icarus21 (1974). For two events these predictions differ by about 20 min from those of Brinkmann and Millis Sky & Telescope45 (1973).  相似文献   

12.
13.
We measured the velocity distributions of impact ejecta with velocities higher than ∼100 m s−1 (high-velocity ejecta) for impacts at variable impact angle α into unconsolidated targets of small soda-lime glass spheres. Polycarbonate projectiles with mass of 0.49 g were accelerated to ∼250 m s−1 by a single-stage light-gas gun. The impact ejecta are detected by thin aluminum foils placed around the targets. We analyzed the holes on the aluminum foils to derive the total number and volume of ejecta that penetrated the aluminum foils. Using the minimum velocity of the ejecta for penetration, determined experimentally, the velocity distributions of the high-velocity ejecta were obtained at α=15°, 30°, 45°, 60°, and 90°. The velocity distribution of the high-velocity ejecta is shown to depend on impact angle. The quantity of the high-velocity ejecta for vertical impact (α=90°) is considerably lower than derived from a power-law relation for the velocity distribution on the low-velocity ejecta (less than 10 m s−1). On the other hand, in oblique impacts, the quantity of the high-velocity ejecta increases with decreasing impact angle, and becomes comparable to those derived from the power-law relation. We attempt to scale the high-velocity ejecta for oblique impacts to a new scaling law, in which the velocity distribution is scaled by the cube of projectile radius (scaled volume) and a horizontal component of impactor velocity (scaled ejection velocity), respectively. The high-velocity ejecta data shows a good correlation between the scaled volume and the scaled ejection velocity.  相似文献   

14.
Impact craters are formed by the displacement and ejection of target material. Ejection angles and speeds during the excavation process depend on specific target properties. In order to quantify the influence of the constitutive properties of the target and impact velocity on ejection trajectories, we present the results of a systematic numerical parameter study. We have carried out a suite of numerical simulations of impact scenarios with different coefficients of friction (0.0–1.0), porosities (0–42%), and cohesions (0–150 MPa). Furthermore, simulations with varying pairs of impact velocity (1–20 km s−1) and projectile mass yielding craters of approximately equal volume are examined. We record ejection speed, ejection angle, and the mass of ejected material to determine parameters in scaling relationships, and to calculate the thickness of deposited ejecta by assuming analytical parabolic trajectories under Earth gravity. For the resulting deposits, we parameterize the thickness as a function of radial distance by a power law. We find that strength—that is, the coefficient of friction and target cohesion—has the strongest effect on the distribution of ejecta. In contrast, ejecta thickness as a function of distance is very similar for different target porosities and for varying impact velocities larger than ~6 km s−1. We compare the derived ejecta deposits with observations from natural craters and experiments.  相似文献   

15.
We present results of polarimetric observations of the Galilean satellites Io, Europa, Ganymede, and Callisto at phase angles ranging from 0.19° to 2.22°. The observations in the UBVR filters were performed using a one-channel photoelectric polarimeter attached to 70-cm telescope of the Chuguev Observation Station (Ukraine) on November 19-December 7, 2000. We have observed the polarization opposition effect for Io, Europa, and Ganymede to be a sharp secondary spike of negative polarization with an amplitude of about −0.4% centered at phase angles of 0.2°-0.7° and superimposed on the regular negative polarization branch. Although these minima for Io, Europa, and Ganymede show many similarities, they also exhibit a number of distinctions. The polarization opposition effect appears to be wavelength-dependent, at least for Europa and Ganymede. No polarization opposition effect was found for Callisto. The results obtained are discussed within the framework of different mechanisms of light scattering.  相似文献   

16.
We present an analysis of the observations of the Deep Impact event performed by the OSIRIS narrow angle camera aboard the Rosetta spacecraft over two weeks, in an effort to characterize the cometary dust grains ejected from the nucleus of Comet 9P/Tempel 1. We adopt a Monte Carlo approach to generate calibrated synthetic images, and a linear combination of them is fitted to the calibrated images so as to determine the physical parameters of the dust cloud. Our model considers spherical olivine particles with a density of 3780 kg m−3. It incorporates constraints on the direction of the cone of emission coming from additional images obtained at Pic du Midi observatory, and constraints on the dust terminal velocities coming from the physics of the impact. We find that the slope of the differential dust size distribution of grains with radii <20 μm (β>0.008) is 3.1±0.3, a value typical of cometary dust tails. This shows that there is no evidence in our data for an enhancement in sub-micron particles in the ejecta compared to the typical dust distribution of active comets. We estimate the mass of particles with radii <1.4 μm (β>0.14) to be 1.5±0.2×105 kg. These particles represent more than 80% of the cross-section of the observed dust cloud. The mass carried by larger particles depends whether the gas significantly increases the kinetic energy of the grains in the inner coma; it lies in the range 1-14×106 kg for particles with radii <100 μm (β>0.002). We obtain the distribution of terminal velocities reached by the dust after the dust-gas interaction which is very well constrained between 10 and 600 m s−1. It is characterized by Gaussian with a maximum at about 190 m s−1 and a width at half maximum of 150 m s−1.  相似文献   

17.
Galileo was the first artificial satellite to orbit Jupiter. During its late orbital mission the spacecraft made two passages through the giant planet’s gossamer ring system. The impact-ionization dust detector on board successfully recorded dust impacts during both ring passages and provided the first in-situ measurements from a dusty planetary ring. During the first passage—on 5 November 2002 while Galileo was approaching Jupiter—dust measurements were collected until a spacecraft anomaly at 2.33RJ (Jupiter radii) just 16 min after a close flyby of Amalthea put the spacecraft into a safing mode. The second ring passage on 21 September 2003 provided ring dust measurements down to about 2.5RJ and the Galileo spacecraft was destroyed shortly thereafter in a planned impact with Jupiter. In all, a few thousand dust impacts were counted with the instrument accumulators during both ring passages, but only a total of 110 complete data sets of dust impacts were transmitted to Earth. Detected particle sizes range from about 0.2 to 5 μm, extending the known size distribution by an order of magnitude towards smaller particles than previously derived from optical imaging [Showalter, M.R., de Pater, I., Verbanac, G., Hamilton, D.P., Burns, J.A., 2008. Icarus 195, 361-377; de Pater, I., Showalter, M.R., Macintosh, B., 2008. Icarus 195, 348-360]. The grain size distribution increases towards smaller particles and shows an excess of these tiny motes in the Amalthea gossamer ring compared to the Thebe ring. The size distribution for the Amalthea ring derived from our in-situ measurements for the small grains agrees very well with the one obtained from images for large grains. Our analysis shows that particles contributing most to the optical cross-section are about 5 μm in radius, in agreement with imaging results. The measurements indicate a large drop in particle flux immediately interior to Thebe’s orbit and some detected particles seem to be on highly-tilted orbits with inclinations up to 20°. Finally, the faint Thebe ring extension was detected out to at least 5RJ, indicating that grains attain higher eccentricities than previously thought. The drop interior to Thebe, the excess of submicron grains at Amalthea, and the faint ring extension indicate that grain dynamics is strongly influenced by electromagnetic forces. These findings can all be explained by a shadow resonance as detailed by Hamilton and Krüger [Hamilton, D.P., Krüger, H., 2008. Nature 453, 72-75].  相似文献   

18.
We present new observations of Jupiter's ring system at a wavelength of 2.2 μm obtained with the 10-m W.M. Keck telescopes on three nights during a ring plane crossing: UT 19 December 2002, and 22 and 26 January 2003. We used conventional imaging, plus adaptive optics on the last night. Here we present detailed radial profiles of the main ring, halo and gossamer rings, and interpret the data together with information extracted from radio observations of Jupiter's synchrotron radiation. The main ring is confined to a 800-km-wide annulus between 128,200 and 129,000 km, with a ∼5000 km extension on the inside. The normal optical depth is 8×10−6, 15% of which is provided by bodies with radii a?5 cm. These bodies are as red as Metis. Half the optical depth, τ≈4×10−6, is attributed to micron-sized dust, and the remaining τ≈3×10−6 to grains tens to hundreds of μm in size. The inward extension consists of micron-sized (a?10 μm) dust, which probably migrates inward under Poynting-Robertson drag. The inner limit of this extension falls near the 3:2 Lorentz resonance (at orbital radius r=122,400 km), and coincides with the outer limit of the halo. The gossamer rings appear to be radially confined, rather than broad sheets of material. The Amalthea ring is triangularly shaped, with a steep outer dropoff over ∼5000 km, extending a few 1000 km beyond the orbit of Amalthea, and a more gradual inner dropoff over 15,000-20,000 km. The inner edge is near the location of the synchronous orbit. The optical depth in the Amalthea ring is ∼5×10−7, up to 20% of which is comprised of macroscopic material. The optical depth in the Thebe ring is a factor of 3 smaller.  相似文献   

19.
Libyan Desert Glass contains meteoritic material and, therefore, its origin is most likely associated with an impact event. However, the impact crater has not been found. We performed numerical simulations of impacts of stony and cometary bodies in order to confirm the version that this glass was formed from silica heated by radiation from aerial bursts near the ground. Asteroids were treated as strengthless bodies from dunite with a density of 3.3 g cm?3, and comets as icy bodies with a density of 1 g cm?3. The simulations based on hydrodynamic equations included the equations of radiation transfer. Melting and vaporization of a silica target under action of radiation incident on a planar surface were modeled using a one‐dimensional hydrodynamic equation of energy and equations of radiation transfer in two‐flux approximation. We selected those variants of simulations in which a crater is not formed, a fireball touches the earth surface, and the area of a molten target corresponds to the area of the Libyan Desert Glass strewn field. Appropriate options include the impact of an asteroid with a diameter of 300 m, an entry speed of 35 km s?1, and an entry angle of 8°, and cometary bodies with diameters from 150 to 300 m, speeds of 50–70 km s?1, and entry angles from 15° to 45°. Impact options with crater formation are also discussed. The maximum depth of molten silica at ground zero reaches 10 cm with the cometary impacts and 3–4 cm with the asteroidal impact. Melting occurs during a period of time from 50 to 400 s.  相似文献   

20.
Two types of trial three-layer models have been constructed for the satellites Io and Europa. In the models of the first type (Io1 and E1), the cores are assumed to consist of eutectic Fe-FeS melt with the densities ρ 1 = 5.15 g cm?3 (Io1) and 5.2 g cm?3 (E1). In the models of the second type (Io3 and E3), the cores consist of FeS with an admixture of nickel and have the density ρ 1 = 4.6 g cm?3. The approach used here differs from that used previously both in chosen model chemical composition of these satellites and in boundary conditions imposed on the models. The most important question to be answered by modeling the internal structure of the Galilean satellites is that of the condensate composition at the formation epoch of Jupiter’s system. Jupiter’s core and the Galilean satellites were formed from the condensate. Ganymede and Callisto were formed fairly far from Jupiter in zones with temperatures below the water condensation temperature, water was entirely incorporated into their bodies, and their modeling showed the mass ratio of the icy (I) component to the rock (R) component in them to be I/R ~ 1. The R composition must be clarified by modeling Io and Europa. The models of the second type (Io3 and E3), in which the satellite cores consist of FeS, yield 25.2 (Io3) and 22.8 (E3) for the core masses (in weight %). In discussing the R composition, we note that, theoretically, the material of which the FeS+Ni core can consist in the R accounts for ~25.4% of the satellite mass. In this case, such an important parameter as the mantle silicate iron saturation is Fe# = 0.265. The Io3 and E3 models agree well with this theoretical prediction. The models of the first and second types differ markedly in core radius; thus, in principle, the R composition in the formation zone of Jupiter’s system can be clarified by geophysical studies. Another problem studied here is that of the error made in modeling Io and Europa using the Radau-Darvin formula when passing from the Love number k 2 to the nondimensional polar moment of inertia $\bar C$ . For Io, the Radau-Darvin formula underestimates the true value of $\bar C$ by one and a half units in the third decimal digit. For Europa, this effect is approximately a factor of 3 smaller, which roughly corresponds to a ratio of the small parameters for the satellites under consideration α Io/α Europa ~ 3.4. In modeling the internal structure of the satellites, the core radius depends strongly on both the mean moment of inertia I* and k 2. Therefore, the above discrepancy in $\bar C$ for Io is appreciable.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号