首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Wavenumber spectra of the martian atmosphere covering zonal wavenumbers s=1-6 were obtained as a function of latitude and season for the first time from the temperatures measured by the Thermal Emission Spectrometer onboard the Mars Global Surveyor. The stationary component tends to peak at s=2, where the martian topography has large amplitude, and drops rapidly at higher wavenumbers. The transient component in the middle and high latitudes tends to peak at s=1, which is lower than the most unstable wavenumber based on linear theories, and exhibits spectral slopes much flatter than the stationary component. In the equatorial region, the spectra of the transient component are almost flat, indicating that the organization of large-scale structures is less efficient in this region. The spectral shapes are similar between the 0.5 and 2.2 hPa surfaces, except that the slopes are slightly steeper at 0.5 than at 2.2 hPa, probably due to selective vertical transmission at low wavenumbers. The seasonal variation is relatively large in the middle and high latitudes, where the maximum power occurs in winter and the minimum occurs in summer, with an exception that the transient component is maximum in spring in the southern hemisphere. Intensification of s=1 transient waves is observed around the period of the initiation of global dust storms.  相似文献   

2.
We have used the complete set of Mars Global Surveyor (MGS) Mars Daily Global Maps (MDGMs) to study martian weather in the southern hemisphere, focusing on curvilinear features, including frontal events and streaks. “Frontal events” refer to visible events that are morphologically analogous to terrestrial baroclinic storms. MDGMs show that visible frontal events were mainly concentrated in the 210-300°E (60-150°W) sector and the 0-60°E sector around the southern polar cap during Ls = 140-250° and Ls = 340-60°. The non-uniform spatial and temporal distributions of activity were also shown by MGS Thermal Emission Spectrometer transient temperature variations near the surface. “Streaks” refer to long curvilinear features in the polar hood or over the polar cap. They are an indicator of the shape of the polar vortex. Streaks in late winter usually show wavy segments between the 180° meridian and Argyre. Model results suggest that the zonal wave number m = 3 eastward traveling waves are important for their formation.  相似文献   

3.
The Viking Orbiters imaged early morning, long, linear wave clouds along the flanks of the Tharsis volcanoes during late northern spring and early summer. These clouds are believed to be a product of either an atmospheric bore wave or a hydraulic jump generated by nightly katabatic winds. The Mars Regional Atmospheric Modeling System was used to study the interaction of the katabatic flows with the surrounding atmosphere to determine what mechanism is responsible for the clouds. Simulations at Ls=90°, 100°, 142°, 180°, 270°, and 358° were conducted focusing on the eastern flank of Olympus Mons. Model results compare well with Viking observations and closely approximate theoretical treatments of atmospheric bores. Strong downslope flows are simulated during the night, with a bore wave forming on and behind a well-defined katabatic front. The observed seasonality of the clouds was reproduced in the simulations; the bore was deeper and faster during northern summer and weakest during the winter. When the bore was strong, it was undular in form, and generated vertically propagating gravity waves in the atmosphere above. During the winter, the atmospheric structure was such that any gravity waves generated damped with height. Less atmospheric water vapor abundance during northern winter, as compared to the summer, is also a factor in the seasonality of the wave clouds. This study concludes that bore waves are the most likely mechanism for the generation of the observed linear wave clouds.  相似文献   

4.
The interval from Ls = 330° in Mars Year (MY) 26 until Ls = 84° in MY 27 has been used to compare and validate measurements from the Mars Global Surveyor Thermal Emission Spectrometer (TES) and the Mars Express Planetary Fourier Spectrometer (PFS). We studied differences between atmospheric temperatures observed by the two instruments. The best agreement between atmospheric temperatures was found at 50 Pa between 40°S and 40°N latitude, where differences were within ±5 K. For other atmospheric levels, differences as large as ∼25 K were observed between the two instruments at some locations. The largest temperature differences occurred mainly over the Hellas Planitia, Argyre Planitia, Tharsis and Valles Marineris regions.On this basis we report on the variability of the martian atmosphere during the 5.5 martian years of Mars climatology obtained by combining the two data sets from TES and PFS. Atmospheric temperatures at 50 Pa responded to the global-scale dust storms of MY 25 and in MY 28 raising temperatures from ∼220 K to ∼250 K during the daytime. An atmospheric temperature of ∼140 K at 50 Pa was observed poleward of 70°N during northern winter and poleward of 60°S during southern winter each year in both the PFS and TES results. Water vapor observed by the two spectrometers showed consistent seasonal and latitudinal variations.  相似文献   

5.
We have characterized the annual behavior of martian atmospheric traveling waves in the MGS TES data set from the first two martian years of mapping. There is a high degree of repeatability between the two years. They are dominated by strong low zonal wavenumber waves with high amplitudes near the polar jets, strongest in late northern fall and early northern winter. The m=1 waves have amplitudes up to about 20 K, are vertically extended, and occasionally extend even into the tropics. Periods for m=1 range from 2.5 to 30 sols. Much weaker waves were identified in the south, with amplitudes less than about 3.5 K. Traveling waves with m=2 and m=3 are also seen, but their amplitudes are typically limited to less than 4 K, and are generally more confined near the surface. In the north, they are more evident in fall and spring rather than winter solstice, which is clearly dominated by m=1 waves. Some evidence of storm tracks has been identified in the data, with accentuated weather-related temperature perturbations near longitudes 200° to 320° E for both the southern and northern hemispheres near latitude ±65° at the surface. Some evidence was also found for a sharpening of longitudinal gradients into what may be frontal systems. EP flux divergences show the waves extracting energy from the zonal mean winds. When the m=1 waves were strongest, decelerations of the zonal jet of order 30 m/(s sol) were measured. Above 1 scale height, the waves extract energy from the jet predominately through barotropic processes, but their character is overall mixed barotropic/baroclinic. Inertial instabilities may exist at altitude on the equatorward flanks of the polar jets, and marginal stability extends through to the tropics. This may explain the coordination of the tropical behavior of the waves with that centered along the polar jet, consistent with the ideas expressed in Wilson et al. (2002, Geophys. Res. Lett. 29, #1684) and similar to those in Barnes et al. (1993, J. Geophys. Res. 98, 3125-3148). Throughout the year, there exist large regions with the meridional gradient of PV less than zero, but they are strongest near winter solstice. Poleward of the winter jet, the regions of instability reach the surface, equatorward they do not. These regions, satisfying a necessary criterion for instability, likely explain the genesis of the waves, and perhaps also their bimodal character between surface (faster waves) and altitude (slow m=1 waves).  相似文献   

6.
The middle atmospheric dynamics on Venus are investigated using a middle atmosphere general circulation model. The magnitude of the superrotation is sensitive to the amplitude of the planetary-scale waves. In particular, the critical level absorptions of the forced planetary-scale waves might contribute to the maintenance of the superrotation near the cloud base. In the case of strong 5.5-day wave forcing, the superrotation with zonal wind speed higher than 100 m s?1 is maintained by the forced wave. Four-day and 5.5-day waves are found near the equatorial cloud top and base, respectively. The planetary-scale waves have a Y-shaped pattern maintained by the amplitude modulation in the presence of strong thermal tides.The polar hot dipole is unstable and its dynamical behavior is complex near the cloud top in this model. The dipole merges into a monopole or breaks up into a tripole when the divergent eddies with high zonal wavenumbers are predominant in the hot dipole region. A cold collar is partly enhanced by a cold phase of slowly propagating waves with zonal wavenumber 1. Although such a complex dipole behavior has not been observed yet, it is likely to occur under a dynamical condition similar to the present simulation. Thus, the dynamical approach using a general circulation model might be useful for analyzing Venus Express and ground-based observation data.  相似文献   

7.
David P. Hinson  Huiqun Wang 《Icarus》2010,206(1):290-1294
We have investigated the near-surface meteorology in the northern hemisphere of Mars through detailed analysis of data obtained with Mars Global Surveyor in January-August 2005. The season in the northern hemisphere ranged from midsummer through winter solstice of Mars Year (MY) 27. We examined composite, wide-angle images from the Mars Orbiter Camera and compiled a catalog of the dust storms that occurred in this interval. As in previous martian years, activity in the northern hemisphere was dominated by regional “flushing” dust storms that sweep southward through the major topographic basins, most frequently in Acidalia Planitia. We also used atmospheric profiles retrieved from radio occultation experiments to characterize eddy activity near the surface at high northern latitudes. There are strong correlations between the two sets of observations, which allowed us to identify three factors that influence the timing and location of the regional dust storms: (1) transitions among baroclinic wave modes, which strongly modulate the intensity of meridional winds near the surface, (2) storms zones, which impose strong zonal variations on the amplitude of some baroclinic eddies, and (3) stationary waves, which further modulate the wind field near the surface. The flushing dust storms ceased abruptly in midautumn, possibly in response to source depletion, CO2 condensation, a shift in the period of the baroclinic eddies, and changes in the tidal wind field near the surface. Our results extend the meteorological record of the northern hemisphere, substantiate the findings of previous investigations, and further illuminate the climatic impact of baroclinic eddies.  相似文献   

8.
Huiqun Wang 《Icarus》2007,189(2):325-343
Data from the third Mars Global Surveyor (MGS) mapping year (MY 26, 2003-2005) are used to investigate dust storms originating in the northern hemisphere. Flushing dust storms, which originate as frontal dust storms at the northern polar vortex edge and propagate southward through topographic channels, are observed immediately before and after a quiescent period that occurs around the northern winter solstice (240°<Ls<300°). Both the pre- and post-solstice active periods can be further divided into two sub-periods. The most vigorous of these flushing storms occurred during Ls 210-220° and Ls 310-320°. The lifted dust crossed the equator and accumulated in the southern hemisphere. These major dust storms enhanced the Hadley circulation and suppressed the lower-level baroclinic eddies in the northern mid and high latitudes. The 2-3 sol wave number m=3 traveling waves show the best correlation with flushing dust storms and can combine with other wave modes to produce storm tracks and fronts within individual sub-periods.  相似文献   

9.
Terry Z. Martin 《Icarus》1981,45(2):427-446
A Mars average data set (MADS) has been constructed from thermal and albedo measurements of the Viking Infrared Thermal Mapper; by merging information from all longitudes; and, ensuring reasonably complete longitudinal sampling, a representation of mean Mars behavior is obtained. Brightness temperatures at 7, 9, 11, 15, and 20 μm and albedo information in the band 0.3–3.0 μm have been binned using 2° latitude strips, 24 times of day, 3 emission angle intervals, and 23 nonoverlapping Ls periods covering 1.43 Mars years starting at Ls = 84°. The MADS is ideally suited to parametric study of latitudinal, diurnal, angular, and seasonal dependences. Data are presented for surface thermal and albedo behavior in clear and dusty atmospheric conditions; the thermal response of the atmospheric temperature to a major dust storm is found to be consistent with Mariner 9 data from the 1971 storm. Examples of use of the MADS, which is available through the Mars Consortium, indicate how averaged data reveal specific surface and atmospheric phenomena.  相似文献   

10.
The seasonal variation of neutron emissions from Mars in different spectral intervals measured by the HEND neutron detector for the entire Martian year are analyzed. Based on these data, the spatial variations of the neutron emissions from the planet are globally mapped as a function of season, and the dynamics of seasonal variation of neutron fluxes with different energies is analyzed in detail. No differences were found between seasonal regimes of neutron fluxes in different energy ranges in the southern hemisphere of Mars, while the regime of fast neutrons (with higher energies) during the northern winter strongly differs from that during the southern winter. In winter (L s = 270°–330°), the fast neutron fluxes are noticeably reduced in the northern hemisphere (along with the consecutive thickening of the seasonal cap of solid carbon dioxide). This provides evidence of a temporary increase in the water content in the effective layer of neutron generation. According to the obtained estimates, the observed reduction of the flux of fast neutrons in the effective layer corresponds to an increase in the water abundance of up to 5% in the seasonal polar cap (70°–90°N), about 3% at mid-latitudes, and from 1.5 to 2% at low latitudes. The freezing out of atmospheric water at the planetary surface (at middle and high latitudes) and the hydration of salt minerals composing the Martian soil are considered as the main processes responsible for the temporary increase in the water content in the soil and upper layer of the seasonal polar cap. The meridional atmospheric transport of water vapor from the summer southern to the winter northern hemisphere within the Hadley circulation cell is a basic process that delivers water to the subsurface soil layer and ensures the observed scale of the seasonal increase in water abundance. In the summer northern hemisphere, the similar Hadley circulation cell transports mainly dry air masses to the winter southern hemisphere. The point is that the water vapor becomes saturated at lower heights during aphelion, and the bulk of the atmospheric water mass is captured in the near-equatorial cloudy belt and, thus, is only weakly transferred to the southern hemisphere. This phenomenon, known as the Clancy effect, was suggested by Clancy et al. (1996) as a basic mechanism for the explanation of the interhemispheric asymmetry of water storage in permanent polar caps. The asymmetry of seasonal meridional circulation of the Martian atmosphere seems to be another factor determining the asymmetry of the seasonal water redistribution in the “atmosphere-regolith-seasonal polar caps” system, found in the peculiarities of the seasonal regime of the neutron emission of Mars.  相似文献   

11.
Michael D Smith 《Icarus》2004,167(1):148-165
We use infrared spectra returned by the Mars Global Surveyor Thermal Emission Spectrometer (TES) to retrieve atmospheric and surface temperature, dust and water ice aerosol optical depth, and water vapor column abundance. The data presented here span more than two martian years (Mars Year 24, Ls=104°, 1 March 1999 to Mars Year 26, Ls=180°, 4 May 2003). We present an overview of the seasonal (Ls), latitudinal, and longitudinal dependence of atmospheric quantities during this period, as well as an initial assessment of the interannual variability in the current martian climate. We find that the perihelion season (Ls=180°-360°) is relatively warm, dusty, free of water ice clouds, and shows a relatively high degree of interannual variability in dust optical depth and atmospheric temperature. On the other hand, the aphelion season (Ls=0°-180°) is relatively cool, cloudy, free of dust, and shows a low degree of interannual variability. Water vapor abundance shows a moderate amount of interannual variability at all seasons, but the most in the perihelion season. Much of the small amount of interannual variability that is observed in the aphelion season appears to be caused by perihelion-season planet-encircling dust storms. These dust storms increase albedo through deposition of bright dust on the surface causing cooler daytime surface and atmospheric temperatures well after dust optical depth returns to prestorm values.  相似文献   

12.
Barney J. Conrath 《Icarus》1981,48(2):246-255
Wave-like perturbations are found in the Mariner 9 IRIS atmospheric temperature data during late Northern Hemisphere winter in a latitude band between 45°N and 65°N. The nature of the data base prevents a unique separation of spatial and temporal behavior, but Fourier analysis of the data constrains the waves to discrete combinations of planetary wavenumber and period. One major spectral component possesses a meridional amplitude cross section with a maximum near the 1-mbar level at 60°N and is strongly correlated with the circumpolar jet observed in thermal winds calculated from the mean meridional temperature cross section. This feature is consistent with the low-wavenumber baroclinic waves observed in Viking Lander data, and the vertical structure reflects the behavior anticipated for a vertically penetrating quasi-geostrophic disturbance. Other possible origins for the wave cannotbe ruled out, however. Among these is a stationary wave forced by wavenumber-2 topographic relief.  相似文献   

13.
Dry convective instabilities in Mars’s middle atmosphere are detected and mapped using temperature retrievals from Mars Climate Sounder observations spanning 1.5 martian years. The instabilities are moderately frequent in the winter extratropics. The frequency and strength of middle atmospheric convective instability in the northern extratropics is significantly higher in MY 28 than in MY 29. This may have coupled with changes to the northern hemisphere mid-latitude and tropical middle atmospheric temperatures and contributed to the development of the 2007 global dust storm. We interpret these instabilities to be the result of gravity waves saturating within regions of low stability created by the thermal tides. Gravity wave saturation in the winter extratropics has been proposed to provide the momentum lacking in general circulation models to produce the strong dynamically-maintained temperature maximum at 1-2 Pa over the winter pole, so these observations could be a partial control on modeling experiments.  相似文献   

14.
In this paper, we have analyzed neutron spectroscopy data gathered by the High Energy Neutron Detector (HEND) instrument onboard Mars Odyssey for comparison of polar regions. It is known that observation of the neutron albedo of Mars provides important information about the distribution of water-ice in subsurface layers and about peculiarities of the CO2 seasonal cycle. It was found that there are large water-rich permafrost areas with contents of up to ∼50% water by mass fraction at both the north and south Mars polar regions. The water-ice layers at high northern latitudes are placed close to the surface, but in the south they are covered by a dry and relatively thick (10-20 cm) layer of soil. Analysis of temporal variations of neutron flux between summer and winter seasons allowed the estimation of the masses of the CO2 deposits which seasonally condense at the polar regions. The total mass of the southern seasonal deposition was estimated as 6.3×1015 kg, which is larger than the total mass of the seasonal deposition at the north by 40-50%. These results are in good agreement with predictions from the NASA Ames Research Center General Circulation Model (GCM). But, the dynamics of the condensation and sublimation processes are not quite as consistent with these models: the peak accumulation of the condensed mass of CO2 occurred 10-15 degrees of Ls later than is predicted by the GCM.  相似文献   

15.
Bruce A. Cantor 《Icarus》2007,186(1):60-96
From 15 September 1997 through 21 January 2006, only a single planet-encircling martian dust storm was observed by MGS-MOC. The onset of the storm occurred on 26 June 2001 (Ls=184.7°), earliest recorded to date. It was initiated in the southern mid-to-low latitudes by a series of local dust storm pulses that developed along the seasonal cap edge in Malea and in Hellas basin (Ls=176.2°-184.4°). The initial expansion of the storm, though asymmetric, was very rapid in all directions (3-32 m s−1). The main direction of propagation, however, was to the east, with the storm becoming planet encircling in the southern hemisphere on Ls=192.3°. Several distinct centers of active dust lifting were associated with the storm, with the longest persisting for 86 sols (Syria-Claritas). These regional storms helped generate and sustain a dust cloud (“haze”), which reached an altitude of about 60 km and a peak opacity of τdust∼5.0. By Ls=197.0°, the cloud had encircled the entire planet between 59.0° S and 60.0° N, obscuring all but the largest volcanoes. The decay phase began around Ls∼200.4° with atmospheric dust concentrations returning to nominal seasonal low-levels at Ls∼304.0°. Exponential decay time constants ranged from 30-117 sols. The storm caused substantial regional albedo changes (darkening and brightening) as a result of the redistribution (removal and deposition) of a thin veneer of surface dust at least 0.1-11.1 μm thick. It also caused changes in meteorological phenomena (i.e., dust storms, dust devils, clouds, recession of the polar caps, and possibly surface temperatures) that persisted for just a few weeks to more than a single Mars year. The redistribution of dust by large annual regional storms might help explain the long period (∼30 years) between the largest planet-encircling dust storms events.  相似文献   

16.
17.
We study the propagation of gravity waves in the martian atmosphere using a linearized one-dimensional full-wave model. Calculations are carried out for atmospheric parameters characteristic of Mars Orbiter Laser Altimeter (on Mars Global Surveyor MGS) observations of apparent gravity waves in high latitude clouds and MGS radio occultation measurements of temperature variations with height suggestive of gravity wave activity. Waves that reach the thermosphere produce fluctuations in density comparable in amplitude with the density variations detected in Mars Odyssey aerobraking data. Gravity waves of modest amplitude are found to deposit momentum and generate significant heating and cooling in the martian atmosphere. The largest heating and cooling effects occur in the thermosphere, at altitudes between about 130 and 150 km, with heating occurring at the lower altitudes and cooling taking place above.  相似文献   

18.
We investigate the interaction of three consecutive large-scale coronal waves with a polar coronal hole, simultaneously observed on-disk by the Solar TErrestrial Relations Observatory (STEREO)-A spacecraft and on the limb by the PRoject for On-Board Autonomy 2 (PROBA2) spacecraft on 27 January 2011. All three extreme ultraviolet (EUV) waves originate from the same active region, NOAA 11149, positioned at N30E15 in the STEREO-A field of view and on the limb in PROBA2. For the three primary EUV waves, we derive starting velocities in the range of ≈?310 km?s?1 for the weakest up to ≈?500 km?s?1 for the strongest event. Each large-scale wave is reflected at the border of the extended coronal hole at the southern polar region. The average velocities of the reflected waves are found to be smaller than the mean velocities of their associated direct waves. However, the kinematical study also reveals that in each case the ending velocity of the primary wave matches the initial velocity of the reflected wave. In all three events, the primary and reflected waves obey the Huygens–Fresnel principle, as the incident angle with ≈?10° to the normal is of the same magnitude as the angle of reflection. The correlation between the speed and the strength of the primary EUV waves, the homologous appearance of both the primary and the reflected waves, and in particular the EUV wave reflections themselves suggest that the observed EUV transients are indeed nonlinear large-amplitude MHD waves.  相似文献   

19.
We show that the peak velocity of Jupiter’s visible-cloud-level zonal winds near 24°N (planetographic) increased from 2000 to 2008. This increase was the only change in the zonal velocity from 2000 to 2008 for latitudes between ±70° that was statistically significant and not obviously associated with visible weather. We present the first automated retrieval of fast (∼130 m s−1) zonal velocities at 8°N planetographic latitude, and show that some previous retrievals incorrectly found slower zonal winds because the eastward drift of the dark projections (associated with 5-μm hot spots) “fooled” the retrieval algorithms.We determined the zonal velocity in 2000 from Cassini images from NASA’s Planetary Data System using a global method similar to previous longitude-shifting correlation methods used by others, and a new local method based on the longitudinal average of the two-dimensional velocity field. We obtained global velocities from images acquired in May 2008 with the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST). Longer-term variability of the zonal winds is based on comparisons with published velocities based on 1979 Voyager 2 and 1995-1998 HST images. Fluctuations in the zonal wind speeds on the order of 10 m s−1 on timescales ranging from weeks to months were found in the 1979 Voyager 2 and the 1995-1998 HST velocities. In data separated by 10 h, we find that the east-west velocity uncertainty due to longitudinal fluctuations are nearly 10 m s−1, so velocity fluctuations of 10 m s−1 may occur on timescales that are even smaller than 10 h. Fluctuations across such a wide range of timescales limit the accuracy of zonal wind measurements. The concept of an average zonal velocity may be ill-posed, and defining a “temporal mean” zonal velocity as the average of several zonal velocity fields spanning months or years may not be physically meaningful.At 8°N, we use our global method to find peak zonal velocities of ∼110 m s−1 in 2000 and ∼130 m s−1 in 2008. Zonal velocities from 2000 Cassini data produced by our local and global methods agree everywhere, except in the vicinity of 8°N. There, the local algorithm shows that the east-west velocity has large variations in longitude; vast regions exceed ∼140 m s−1. Our global algorithm, and all of the velocity-extraction algorithms used in previously-published studies, found the east-west drift velocities of the visible dark projections, rather than the true zonal velocity at the visible-cloud level. Therefore, the apparent increase in zonal winds between 2000 and 2008 at 8°N is not a true change in zonal velocity.At 7.3°N, the Galileo probe found zonal velocities of 170 m s−1 at the 3-bar level. If the true zonal velocity at the visible-cloud level at this latitude is ∼140 m s−1 rather than ∼105 m s−1, then the vertical zonal wind shear is much less than the currently accepted value.  相似文献   

20.
We present a seismological method for probing the solar atmosphere above sunspot umbrae with three-minute oscillations. Our technique allows us to estimate both the vertical distance between atmospheric layers and the wave-propagation speed, without specifying any additional parameters, in particular, the phase speed of the wave or the emission formation heights. Our method uses the projected wave paths of slow magnetohydrodynamic waves that propagate through the atmospheric layers of different heights and are guided by the magnetic field. The length of the projected wave path depends upon the distance between the layers and the inclination angle of the magnetic field with respect to the line of sight, allowing us to estimate the distance between the layers from measured projected wave paths and the local magnetic-field vector. In turn, the wave-propagation delay registered at different heights allows for the calculation of the phase speed. We estimated the vertical distance between the emission layers at the temperature minimum (1600 Å) and transition region (304 Å), as well as the average phase speed above the sunspot umbrae, for three active regions. We found that the distance between the 1600 Å emission layer and the transition region above the sunspot umbrae lies in the range of 500?–?800 km. The average phase speed between these layers was found to be about 30 km?s?1, giving a sound speed of 6 km?s?1. The temperature between the layers has been roughly estimated as 3000 K and corresponds to the region of the temperature minimum. The results obtained are consistent with the semiempirical model of the sunspot-umbrae atmosphere by Fontenla et al. (Astrophys. J.707, 482, 2009).  相似文献   

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

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