首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 437 毫秒
1.
In this study we examine a set of lunar mare domes located in the Hortensius/Milichius/T. Mayer region and in northern Mare Tranquillitatis with respect to their formation along crustal fractures, their rheologic properties, the dimensions of their feeder dikes, and the importance of magma evolution processes during dome formation. Many of these domes display elongated summit vents oriented radially with respect to major impact basins, and several dome locations are also aligned in these preferential directions. Analysis of Clementine UV/VIS and Lunar Prospector gamma ray spectrometer data reveals that the examined mare domes formed from low-Si basaltic lavas of high FeO and low to moderate TiO2 content. Based on their morphometric properties (diameter, height, volume) obtained by photoclinometric and shape from shading analysis of telescopic CCD images, we derive rheologic quantities (lava viscosity during eruption, effusion rate, duration of the effusion process, magma rise speed) and the dimensions of the feeder dikes. We establish three rheologic groups characterised by specific combinations of rheologic properties and dike dimensions, where the most relevant discriminative parameter is the lava viscosity η. The first group is characterised by and contains the domes with elongated vents in the Milichius/T. Mayer region and two similar domes in northern Mare Tranquillitatis. The second group with comprises the very low aligned domes in northern Mare Tranquillitatis, and the third group with the relatively steep domes near Hortensius and in the T. Mayer region. The inferred dike dimensions in comparison to lunar crustal thickness data indicate that the source regions of the feeder dikes are situated within the upper crust for six of the domes in northern Mare Tranquillitatis, while they are likely to be located in the lower crust and in the upper mantle for the other examined domes. By comparing the time scale of magma ascent with the time scale on which heat is conducted from the magma into the host rock, we find evidence that the importance of magma evolution processes during ascent such as cooling and crystallisation increases with lava viscosity. We conclude that different degrees of evolution of initially fluid basaltic magma are able to explain the broad range of lava viscosities inferred for the examined mare domes. The spectral data reveal that differences in TiO2 content may additionally account for the systematic difference in lava viscosity between the two examined lunar regions. We show that the described mechanisms are likely to be valid also for other lunar mare domes situated near Cauchy and Arago, regarded for comparison. On the other hand, we find for the Gruithuisen and Mairan highland domes that despite their inferred high lava viscosities of , no significant magma cooling in the dike occurred during ascent, supporting previous findings that the highland domes were formed during a specific phase of non-mare volcanism by highly silicic viscous lavas.  相似文献   

2.
In this study we examine the spectral and morphometric properties of the four important lunar mare dome fields near Cauchy, Arago, Hortensius, and Milichius. We utilize Clementine UV-vis multispectral data to examine the soil composition of the mare domes while employing telescopic CCD imagery to compute digital elevation maps in order to determine their morphometric properties, especially flank slope, height, and edifice volume. After reviewing previous attempts to determine topographic data for lunar domes, we propose an image-based 3D reconstruction approach which is based on a combination of photoclinometry and shape from shading. Accordingly, we devise a classification scheme for lunar mare domes which is based on a principal component analysis of the determined spectral and morphometric features. For the effusive mare domes of the examined fields we establish four classes, two of which are further divided into two subclasses, respectively, where each class represents distinct combinations of spectral and morphometric dome properties. As a general trend, shallow and steep domes formed out of low-TiO2 basalts are observed in the Hortensius and Milichius dome fields, while the domes near Cauchy and Arago that consist of high-TiO2 basalts are all very shallow. The intrusive domes of our data set cover a wide continuous range of spectral and morphometric quantities, generally characterized by larger diameters and shallower flank slopes than effusive domes. A comparison to effusive and intrusive mare domes in other lunar regions, highland domes, and lunar cones has shown that the examined four mare dome fields display such a richness in spectral properties and 3D dome shape that the established representation remains valid in a more global context. Furthermore, we estimate the physical parameters of dome formation for the examined domes based on a rheologic model. Each class of effusive domes defined in terms of spectral and morphometric properties is characterized by its specific range of values for lava viscosity, effusion rate, and duration of the effusion process. For our data set we report lava viscosities between about 102 and , effusion rates between 25 and , and durations of the effusion process between three weeks and 18 years. Lava viscosity decreases with increasing R415/R750 spectral ratio and thus TiO2 content; however, the correlation is not strong, implying an important influence of further parameters like effusion temperature on lava viscosity.  相似文献   

3.
The non-thermal escape of neutral O atoms from Mars at the current epoch is largely due to dissociative recombination of :
  相似文献   

4.
The giant planets of our solar system possess envelopes consisting mainly of hydrogen and helium but are also significantly enriched in heavier elements relatively to our Sun. In order to better constrain how these heavy elements have been delivered, we quantify the amount accreted during the so-called “late heavy bombardment”, at a time when planets were fully formed and planetesimals could not sink deep into the planets. On the basis of the “Nice model”, we obtain accreted masses (in terrestrial units) equal to for Jupiter, and for Saturn. For the two other giant planets, the results are found to depend mostly on whether they switched position during the instability phase. For Uranus, the accreted mass is with an inversion and without an inversion. Neptune accretes in models in which it is initially closer to the Sun than Uranus, and otherwise. With well-mixed envelopes, this corresponds to an increase in the enrichment over the solar value of 0.033±0.001 and 0.074±0.007 for Jupiter and Saturn, respectively. For the two other planets, we find the enrichments to be 2.1±1.4 (w/ inversion) or 1.2±0.7 (w/o inversion) for Uranus, and 2.0±1.2 (w/ inversion) or 2.7±1.6 (w/o inversion) for Neptune. This is clearly insufficient to explain the inferred enrichments of ∼4 for Jupiter, ∼7 for Saturn and ∼45 for Uranus and Neptune.  相似文献   

5.
We present CCD photometric observations of 23 main-belt asteroids, of which 8 have never been observed before; thus, the data of these objects are the first in the literature. The majority showed well-detectable light variations, exceeding 0m1. We have determined synodic periods for 756 Lilliana (936), 1270 Datura (34), 1400 Tirela (1336), 1503 Kuopio (998), 3682 Welther (359), 7505 Furushu (414) and 11436 1969 QR (123), while uncertain period estimates were possible for 469 Argentina (123), 546 Herodias (104) and 1026 Ingrid (53). The shape of the lightcurves of 3682 Welther changed on a short time-scale and showed dimmings that might be attributed to eclipses in a binary system. For the remaining objects, only lower limits of the periods and amplitudes were concluded.  相似文献   

6.
We present mid-infrared observations of the binary L5-Trojan system (617) Patroclus-Menoetius before, during, and after two shadowing events, using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. For the first time, we effectively observe changes in asteroid surface temperature in real time, allowing the thermal inertia to be determined very directly. A new detailed binary thermophysical model is presented which accounts for the system’s known mutual orbit, arbitrary component shapes, and thermal conduction in the presence of eclipses.We obtain two local thermal-inertia values, representative of the respective shadowed areas: and . The average thermal inertia is estimated to be , potentially with significant surface heterogeneity. This first thermal-inertia measurement for a Trojan asteroid indicates a surface covered in fine regolith. Independently, we establish the presence of fine-grained (<a few μm) silicates on the surface, based on emissivity features near 10 and similar to those previously found on other Trojans.We also report V-band observations and report a lightcurve with complete rotational coverage. The lightcurve has a low amplitude of peak-to-peak, implying a roughly spherical shape for both components, and is single-periodic with a period equal to the period of the mutual orbit, indicating that the system is fully synchronized.The diameters of Patroclus and Menoetius are 106±11 and , respectively, in agreement with previous findings. Taken together with the system’s known total mass, this implies a bulk mass density of , significantly below the mass density of L4-Trojan asteroid (624) Hektor and suggesting a bulk composition dominated by water ice.All known physical properties of Patroclus, arguably the best studied Trojan asteroid, are consistent with those expected in icy objects with devolatilized surface (extinct comets), consistent with what might be implied by recent dynamical modeling in the framework of the Nice Model.  相似文献   

7.
David A. Minton  Renu Malhotra 《Icarus》2010,207(2):744-7225
The cumulative effects of weak resonant and secular perturbations by the major planets produce chaotic behavior of asteroids on long timescales. Dynamical chaos is the dominant loss mechanism for asteroids with diameters in the current asteroid belt. In a numerical analysis of the long-term evolution of test particles in the main asteroid belt region, we find that the dynamical loss history of test particles from this region is well described with a logarithmic decay law. In our simulations the loss rate function that is established at persists with little deviation to at least . Our study indicates that the asteroid belt region has experienced a significant amount of depletion due to this dynamical erosion—having lost as much as ∼50% of the large asteroids—since 1 Myr after the establishment of the current dynamical structure of the asteroid belt. Because the dynamical depletion of asteroids from the main belt is approximately logarithmic, an equal amount of depletion occurred in the time interval 10-200 Myr as in 0.2-4 Gyr, roughly ∼30% of the current number of large asteroids in the main belt over each interval. We find that asteroids escaping from the main belt due to dynamical chaos have an Earth-impact probability of ∼0.3%. Our model suggests that the rate of impacts from large asteroids has declined by a factor of 3 over the last 3 Gyr, and that the present-day impact flux of objects on the terrestrial planets is roughly an order of magnitude less than estimates currently in use in crater chronologies and impact hazard risk assessments.  相似文献   

8.
Using the S-band radar at Arecibo Observatory, we observed six new M-class main-belt asteroids (MBAs), and re-observed one, bringing the total number of Tholen M-class asteroids observed with radar to 19. The mean radar albedo for all our targets is , significantly higher than the mean radar albedo of every other class (Magri, C., Nolan, M.C., Ostro, S.J., Giorgini, J.D. [2007]. Icarus 186, 126-151). Seven of these objects (Asteroids 16 Psyche, 129 Antigone, 216 Kleopatra, 347 Pariana, 758 Mancunia, 779 Nina, 785 Zwetana) have radar albedos indicative of a very high metal content , and consistent with a remnant iron/nickel core interpretation (irons) or exotic high metal meteorite types such as CB. We propose designating these high radar albedo objects as Mm. Two asteroids, 110 Lydia and 678 Fredegundis, have more moderate radar albedos , but exhibit high values at some rotation phases suggesting a significant metal content. The remaining 10 objects have moderate radar albedos at all rotation phases. Most of our targets have visible/near-infrared spectra (Hardersen, P.S., Gaffey, M.J., Abell, P.A. [2005]. Icarus 175, 141-158; Fornasier, S., Clark, B.E., Dotto, E., Migliorini, A., Ockert-Bell, M., Barucci, M.A. [2009]. Icarus, submitted for publication) that indicate the presence of at least some silicate phases. All of the non-Mm asteroids show a positive correlation between visual and radar albedo but the reasons for this are not clear. All of the higher radar albedo targets (the 7 Mm asteroids, Lydia, and Fredegundis) show moderate to large variations in radar albedo with rotation phase. We suggest that their high radar reflectivity exaggerates irregularities in the asteroid shape to cause this behavior. One-third of our targets show evidence for asteroid-scale concavities or bifurcation. Based on all the evidence available, we suggest that most Tholen M-class asteroids are not remnant iron cores or enstatite chondrites, but rather collisional composites of silicates and irons with compositions more analogous to stony-iron meteorites and high-iron carbonaceous chondrites.  相似文献   

9.
We investigate the effects of atmospheric gravity waves on the vertical and horizontal structure of the ionosphere of Jupiter. The presented non-linear, two-dimensional model of the jovian ionosphere allows for spatially and temporally varying neutral wind and temperature fields and tracks the time evolution of six ionospheric species, , and . An analytical approach is used to validate the model results for linear, small-amplitude waves and to elucidate the mechanisms that leads to perturbations in the density of the main ion species, H+ and . We demonstrate that the long-lived H+ ions are perturbed directly by wave dynamics whereas short-lived ions such as are perturbed by chemical interactions with other perturbed ion species. The model is then applied using larger gravity wave amplitudes consistent with observations. Atmospheric gravity waves propagating at high altitudes create layers of enhanced electron density similar to the system of layers observed during the J0-ingress radio occultation of the Galileo spacecraft. Our best fit to the J0-ingress observation is achieved using an 82 min period forcing wave with horizontal and vertical wavelengths of 500 km and 60 km respectively, and peaks at 510 km above the 1 bar pressure level. We further investigate the effects of the wave-induced ion flux on the background ionospheric structure and demonstrate that in the presence of a gravity wave the background density profiles of the H+ and ions are significantly modified. We also find that the column density of has variations that can exceed 10% as the wave propagates.  相似文献   

10.
Jane L. Fox  Aleksander Ha? 《Icarus》2010,208(1):176-191
The production of energetic and escaping neutral O atoms at the current epoch in the martian thermosphere is thought to be dominated by the dissociative recombination process:
  相似文献   

11.
《Icarus》2009,199(2):458-476
On September 8, 2001 around 2 h UT, the largest uranian moon, Titania, occulted Hipparcos star 106829 (alias SAO 164538, a V=7.2, K0 III star). This was the first-ever observed occultation by this satellite, a rare event as Titania subtends only 0.11 arcsec on the sky. The star's unusual brightness allowed many observers, both amateurs or professionals, to monitor this unique event, providing fifty-seven occultations chords over three continents, all reported here. Selecting the best 27 occultation chords, and assuming a circular limb, we derive Titania's radius: (1-σ error bar). This implies a density of using the value derived by Taylor [Taylor, D.B., 1998. Astron. Astrophys. 330, 362-374]. We do not detect any significant difference between equatorial and polar radii, in the limit , in agreement with Voyager limb image retrieval during the 1986 flyby. Titania's offset with respect to the DE405 + URA027 (based on GUST86 theory) ephemeris is derived: ΔαTcos(δT)=−108±13 mas and ΔδT=−62±7 mas (ICRF J2000.0 system). Most of this offset is attributable to a Uranus' barycentric offset with respect to DE405, that we estimate to be: and ΔδU=−85±25 mas at the moment of occultation. This offset is confirmed by another Titania stellar occultation observed on August 1st, 2003, which provides an offset of ΔαTcos(δT)=−127±20 mas and ΔδT=−97±13 mas for the satellite. The combined ingress and egress data do not show any significant hint for atmospheric refraction, allowing us to set surface pressure limits at the level of 10-20 nbar. More specifically, we find an upper limit of 13 nbar (1-σ level) at 70 K and 17 nbar at 80 K, for a putative isothermal CO2 atmosphere. We also provide an upper limit of 8 nbar for a possible CH4 atmosphere, and 22 nbar for pure N2, again at the 1-σ level. We finally constrain the stellar size using the time-resolved star disappearance and reappearance at ingress and egress. We find an angular diameter of 0.54±0.03 mas (corresponding to projected at Titania). With a distance of 170±25 parsecs, this corresponds to a radius of 9.8±0.2 solar radii for HIP 106829, typical of a K0 III giant.  相似文献   

12.
We study the Jupiter family comet (JFC) population assumed to come from the Scattered Disk and transferred to the Jupiter’s zone through gravitational interactions with the Jovian planets. We shall define as JFCs those with orbital periods and Tisserand parameters in the range 2<T?3.1, while those comets coming from the same source, but that do not fulfill the previous criteria (mainly because they have periods ) will be called ‘non-JFCs’. We performed a series of numerical simulations of fictitious comets with a purely dynamical model and also with a more complete dynamical-physical model that includes besides nongravitational forces, sublimation and splitting mechanisms. With the dynamical model, we obtain a poor match between the computed distributions of orbital elements and the observed ones. However with the inclusion of physical effects in the complete model we are able to obtain good fits to observations. The best fits are attained with four splitting models with a relative weak dependence on q, and a mass loss in every splitting event that is less when the frequency is high and vice versa. The mean lifetime of JFCs with radii and is found to be of about 150-200 revolutions (∼. The total population of JFCs with radii within Jupiter’s zone is found to be of 450±50. Yet, the population of non-JFCs with radii in Jupiter-crossing orbits may be ∼4 times greater, thus leading to a whole population of JFCs + non-JFCs of ∼2250±250. Most of these comets have perihelia close to Jupiter’s orbit. On the other hand, very few non-JFCs reach the Earth’s vicinity (perihelion distances ) which gives additional support to the idea that JFCs and Halley-type comets have different dynamical origins. Our model allows us to define the zones of the orbital element space in which we would expect to find a large number of JFCs. This is the first time, to our knowledge, that a physico-dynamical model is presented that includes sublimation and different splitting laws. Our work helps to understand the role played by these erosion effects in the distribution of the orbital elements and lifetimes of JFCs.  相似文献   

13.
Solar light gets scattered at cloud top level in Venus’ atmosphere, in the visible range, which corresponds to the altitude of 67 km. We present Doppler velocity measurements performed with the high resolution spectrometer MTR of the Solar telescope THEMIS (Teide Observatory, Canary Island) on the sodium D2 solar line . Observations lasted only 49 min because of cloudy weather. However, we could assess the instrumental velocity sensitivity, per pixel of 1 arcsec, and give a value of the amplitude of zonal wind at equator at .  相似文献   

14.
Koji Matsuo 《Icarus》2009,202(1):90-578
The martian atmosphere seasonally exchanges CO2 with the surface by repeating condensation and sublimation, causing seasonal growth and decay of the polar CO2 snowcaps. These processes leave two kinds of geodetic signatures, i.e. seasonal changes of the martian gravity field and of surface elevation of the snow-covered regions. Here we study gradual increase of the volume density of the martian snow due to compaction, by combining these two data sets during 1999-2001 covering three martian winters. We found that light fresh snow of slowly becomes denser reaching or more immediately before it thaws. The maximum snow density varies slightly from year to year, and between hemispheres. In the second southern winter, the density became as high as . This might have been caused by a dust storm activity, e.g. increased mixing of silicate particles and/or enhancement of sintering.  相似文献   

15.
We have obtained numerically integrated orbits for Saturn's coorbital satellites, Janus and Epimetheus, together with Saturn's F-ring shepherding satellites, Prometheus and Pandora. The orbits are fit to astrometric observations acquired with the Hubble Space Telescope and from Earth-based observatories and to imaging data acquired from the Voyager spacecraft. The observations cover the 38 year period from the 1966 Saturn ring plane crossing to the spring of 2004. In the process of determining the orbits we have found masses for all four satellites. The densities derived from the masses for Janus, Epimetheus, Prometheus, and Pandora in units of g cm−3 are , , , and , respectively.  相似文献   

16.
The VIMS instrument onboard Cassini observed the north polar region of Titan at 113° phase angle, 28 December 2006. On this spectral image, a vast polar cloud can be seen northward to 62°N, and elsewhere, the haze appears as the dominant source of scattering. Because the surface does not appear in the wavelength range between 0.3 and , this spectro-image is ideal to study airborn scatterers both in methane bands and windows. In this work, we study this image, along with another image taken at 13° phase angle. This image probe both the atmosphere and the surface from pole to pole. First, we characterise the spatial distribution of the haze layer above 100 km between 80°S and 70°N. We find a north south asymmetry with a haze opacity increasing by a factor 3 from the south pole to the equator, then a constant value up to about 30°N and a decrease of a factor 2 between 30°N and about 60°N. Beyond 60°N, we can see the influence of the north polar cloud, even in the band, but no polar haze accumulation. The fact that the north polar region is still in the polar night is a possible explanation. No haze accumulation is observed in the southern polar region. Secondly, we partly identify the origin of spectral features in the 2.8-μm methane window, which are found to be due to deuterated methane (CH3D). This allows the analyse of this window and to retrieve the opacity of scatterers layer below 80 km (presumably made of aerosols and condensate droplets) between 35°N and 60°N. Finally, we constrained the values and the spectral behaviour of the imaginary part of the aerosol refractive index in the range between 0.3 and . To do so, we studied the 2.8-μm window with the image taken at 113° phase angle. To complete the analysis, we studied the transmission through the haze layer in the 3.4-μm band observed in solar occultation mode with VIMS, and we analysed the single scattering albedo retrieved with DISR instrument between 0.4 and . The imaginary part of the refractive index that we find for Titan aerosols follows Khare et al. (Khare, B.N. et al. [1984]. Icarus 60, 127-137) optical constant up to and becomes constant beyond this wavelength at least up to . It also has a prominent peak at and a secondary peak at , which indicates material rich in C-H bonds, with much less N-H bonds than in Khare et al. (1984) tholins.  相似文献   

17.
18.
The Alice ultraviolet spectrograph onboard the New Horizons spacecraft observed two occultations of the bright star χ Ophiucus by Jupiter’s atmosphere on February 22 and 23, 2007 during the approach phase of the Jupiter flyby. The ingress occultation probed the atmosphere at 32°N latitude near the dawn terminator, while egress probed 18°N latitude near the dusk terminator. A detailed analysis of both the ingress and egress occultations, including the effects of molecular hydrogen, methane, acetylene, ethylene, and ethane absorptions in the far ultraviolet (FUV), constrains the eddy diffusion coefficient at the homopause level to be  cm2 s−1, consistent with Voyager measurements and other analyses (Festou, M.C., Atreya, S.K., Donahue, T.M., Sandel, B.R., Shemansky, D.E., Broadfoot, A.L. [1981]. J. Geophys. Res. 86, 5717-5725; Vervack Jr., R.J., Sandel, B.R., Gladstone, G.R., McConnell, J.C., Parkinson, C.D. [1995]. Icarus 114, 163-173; Yelle, R.V., Young, L.A., Vervack Jr., R.J., Young, R., Pfister, L., Sandel, B.R. [1996]. J. Geophys. Res. 101 (E1), 2149-2162). However, the actual derived pressure level of the methane homopause for both occultations differs from that derived by [Festou et al., 1981] and [Yelle et al., 1996] from the Voyager ultraviolet occultations, suggesting possible changes in the strength of atmospheric mixing with time. We find that at 32°N latitude, the methane concentration is  cm−3 at 70,397 km, the methane concentration is  cm−3 at 70,383 km, the acetylene concentration is  cm−3 at 70,364 km, and the ethane concentration is  cm−3 at 70,360 km. At 18°N latitude, the methane concentration is  cm−3 at 71,345 km, the methane concentration is  cm−3 at 71,332 km, the acetylene concentration is cm−3 at 71,318 km, and the ethane concentration is  cm−3 at 71,315 km. We also find that the H2 occultation light curve is best reproduced if the atmosphere remains cold in the microbar region such that the base of the thermosphere is located at a lower pressure level than that determined by in situ instruments aboard the Galileo probe (Seiff, A., Kirk, D.B., Knight, T.C.D., Young, R.E., Mihalov, J.D., Young, L.A., Milos, F.S., Schubert, G., Blanchard, R.C., Atkinson, D. [1998]. J. Geophys. Res. 103 (E10), 22857-22889) - the Sieff et al. temperature profile leads to too much absorption from H2 at high altitudes. However, this result is highly model dependent and non-unique. The observations and analysis help constrain photochemical models of Jupiter’s atmosphere.  相似文献   

19.
Numerical simulations of asteroid breakups, including both the fragmentation of the parent body and the gravitational interactions between the fragments, have allowed us to reproduce successfully the main properties of asteroid families formed in different regimes of impact energy, starting from a non-porous parent body. In this paper, using the same approach, we concentrate on a single regime of impact energy, the so-called catastrophic threshold usually designated by , which results in the escape of half of the target’s mass. Thanks to our recent implementation of a model of fragmentation of porous materials, we can characterize for both porous and non-porous targets with a wide range of diameters. We can then analyze the potential influence of porosity on the value of , and by computing the gravitational phase of the collision in the gravity regime, we can characterize the collisional outcome in terms of the fragment size and ejection speed distributions, which are the main outcome properties used by collisional models to study the evolutions of the different populations of small bodies. We also check the dependency of on the impact speed of the projectile.In the strength regime, which corresponds to target sizes below a few hundreds of meters, we find that porous targets are more difficult to disrupt than non-porous ones. In the gravity regime, the outcome is controlled purely by gravity and porosity in the case of porous targets. In the case of non-porous targets, the outcome also depends on strength. Indeed, decreasing the strength of non-porous targets make them easier to disrupt in this regime, while increasing the strength of porous targets has much less influence on the value of . Therefore, one cannot say that non-porous targets are systematically easier or more difficult to disrupt than porous ones, as the outcome highly depends on the assumed strength values. In the gravity regime, we also confirm that the process of gravitational reaccumulation is at the origin of the largest remnant’s mass in both cases. We then propose some power-law relationships between and both target’s size and impact speed that can be used in collisional evolution models. The resulting fragment size distributions can also be reasonably fitted by a power-law whose exponent ranges between −2.2 and −2.7 for all target diameters in both cases and independently on the impact velocity (at least in the small range investigated between 3 and 5 km/s). Then, although ejection velocities in the gravity regime tend to be higher from porous targets, they remain on the same order as the ones from non-porous targets.  相似文献   

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

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