首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
A simple technique has been used to investigate gravity waves in the Antarctic and Arctic mesosphere. Data from two meteor radars are used, one at Rothera in the Antarctic (68°S) and one at Esrange in the Arctic (68°N). Observations have been made from 2005 to 2008 in the Antarctic and from 2000 to 2008 in the Arctic. The results reveal a seasonal behaviour with a semi-annual cycle. Over both sites wave activity maximises at the solstices and minimises at the equinoxes. Monthly-mean gravity-wave activity increases with height in all seasons except in summer when gravity-wave variances show little or no increase with height at heights below about 90 km. Above about 90 km in summer there is a rapid increase in wave activity. We suggest that this summer-time behaviour is a consequence of rapid wave growth as gravity waves ascend into the more stable air of the lower thermosphere. Differences between the Arctic and Antarctic gravity-wave field are also evident. In particular, a higher level of gravity-wave activity is observed in the Antarctic spring compared to the Arctic spring. We suggest that these inter-hemispheric differences may be the result of differences in the generation of gravity waves in the polar troposphere and stratosphere.  相似文献   

2.
We analyzed meteor decay times measured by a VHF radar at King Sejong Station by classifying strong and weak meteors according to their estimated electron line densities. The height profiles of monthly averaged decay times show a peak whose altitude varies with season at altitudes of 80?85 km. The higher peak during summer is consistent with colder temperatures that cause faster chemical reactions of electron removal. By adopting temperature dependent empirical recombination rates from rocket experiments and meteor electron densities of 2×105?2×106 cm?3 in a decay time model, we are able to account for decreasing decay times below the peak for all seasons without invoking meteor electron removal by hypothetical icy particles.  相似文献   

3.
We show examples of common volume observations of three metals by lidar focusing on the altitude of the topside of the meteoric metal layer as described by Höffner and Friedman (H&F) [The mesospheric metal layer topside: a possible connection to meteoroids, Atmos. Chem. Phys. 4 (2004) 801–808]. In contrast to H&F, we will focus on time scales of a few hours and less whereas the previous study examined the seasonally averaged climatological state on time scales of several days or weeks, and we examine the entire topside, whereas H&F focused on data at 113 km. The examples, taken under different observation conditions in 1997 and 1998 at Kühlungsborn, Germany (54°N, 15°E), show that the metal layers can often be observed at altitudes as high as 130 km if the signal is integrated over a period of several hours. Under such conditions it is possible to derive reasonably good metal abundance ratios from nocturnally averaged data, which, in turn, allow the discussion of metal abundance ratios to broaden from a single altitude as discussed in H&F to an altitude range extending as high as 130 km. The examples herein show, for the first time, that it is possible to track the transition in the metal abundance ratios from the main layer to an altitude region that has not been studied in the past by lidar. On shorter time scales, small structures are detectable and observable, sometimes above 120 km, resulting in, on average, a broad but weak topside layer above 105 km. In particular, the example of 26–27 October 1997, obtained during enhanced meteor activity, is an indication that this broad layer may result from meteor ablation occurring in this altitude range during the observation. Ratios of metal densities for Ca, Fe, K, and Na are remarkably consistent above about 110 km and in close agreement with the results of H&F. They are less consistent with ratios measured in individual meteor trails and appear to have little relation to the ratios measured in CI meteorites. Finally, it is the temporal smoothing of descending sporadic metal atom layers on top of an undisturbed background metal layer that is the basis of the summer topside extension as described by H&F.  相似文献   

4.
A study on variability of the equatorial ionosphere was carried out at fixed heights below the F2 peak for two different levels of solar activity. The study covered height range of 100 km up to the peak of F2 layer using a real height step increase of 10 km. The variability index used is the percentage ratio of standard deviation over the average value for the month. Daytime minimum variability of between 3% and 10% was observed at height range of about 150–210 km during low solar activity and between 2% and 7% at height range of 160–220 km during high solar activity. The nighttime maximum of between 70% and 187% was observed at height range of about 210–250 km during low solar activity and between 42% and 127% at height range of 210–250 km during high solar activity. The height range at which daytime minimum was observed falls within the F1 height of the ionosphere. The result obtained is consistent with previous works carried out in the low latitude locations for American sector.  相似文献   

5.
To study the occurrence characteristics of equatorial spread-F irregularities and their latitudinal extent, simultaneous digital ionosonde data (January–December 2001) from Trivandrum (8.2°N), Waltair (17.7°N) and Delhi (28.6°N) and 4 GHz scintillation data from Sikandarabad (26.8°N) and Chenglepet (10.4°N), and 250 MHz scintillation data from Bhopal (23.2°N) for equinoxes period are analysed. It is noted that except summer months, occurrence of spread F is always maximum at Trivandrum, minimum at Delhi and moderate at Waltair. During equinoxes and winter months. Their occurrences at higher latitude station are always conditional to their prior occurrences at lower latitudes indicating their association with the generation of equatorial plasma bubble and associated irregularities. Scintillation occurrences also follow the similar pattern. During the summer months, the spread-F occurrences are highest at equatorial location Trivandrum, moderate at Delhi and minimum at Waltair and seem to be caused by irregularities generated locally especially over Delhi.To gain forecasting capability, night-to-night occurrences of spread-F/scintillation at these locations are examined in relation to post sunset rise of h’F and upward ExB drift velocity over the magnetic equator using Trivandrum ionosonde data. It is noted that except the summer months, the spread-F at Trivandrum, Waltair and Delhi are observed only when equatorial ExB (h’F) is more than about 15 m/s (325 km), 20 m/s (350 km) and 25 m/s (375 km), respectively. With these threshold values their corresponding success rate of predictions are more than 90%, 50% and 15% at the respective locations. Whereas in the case of GHz scintillations near equator are observed only when ExB (h’F) is more than 15 m/s (325 km), whereas for low latitude, the same should be 30 m/s (400 km) and their success rate of prediction is about 90% and 30%, respectively. The intensity of 4 GHz scintillation at low latitude is also found to be positively correlated with equatorial upward ExB drift velocity values, whereas correlation is poor with that of equatorial scintillations. In conclusions, near magnetic equator threshold values of ExB or h’F can be successfully used for the night-to-night prediction of spread-F/scintillations occurrences, whereas these are necessary but not sufficient for their prediction at higher latitudes. For that some other controlling parameters like background electron density, neutral winds, gravity waves, etc. should also be examined.  相似文献   

6.
We report observations of a noctilucent cloud (NLC) over central Alaska by a ground-based lidar and camera on the night of 9–10 August 2005. The lidar at Poker Flat Research Range (PFRR), Chatanika (65°N, 147°W) measured a maximum integrated backscatter coefficient of 2.4×10?6 sr?1 with a peak backscatter coefficient of 2.6×10?9 m?1 sr?1 corresponding to an aerosol backscatter ratio of 120 at an altitude of 82.1 km. The camera at Donnelly Dome, 168 km southeast of PFRR, recorded an extensive NLC display across the sky with distinct filamentary features corresponding to wave structures measured by the lidar. The occurrence of the maximum integrated backscatter coefficient corresponded to the passage of a bright cloud band to the southwest over PFRR. The camera observations indicate that the cloud band had a horizontal width of 50 km and a length of 150 km. The horizontal scale of the cloud band was confirmed by medium-frequency radar wind measurements that reported mesopause region winds of 30 m/s to the southwest during the period when the cloud band passed over PFRR. Comparison of these measurements with current NLC microphysical models suggests a lower bound on the water vapor mixing ratio at 83 km of 7–9 ppmv and a cloud ice mass of 1.5–1.8×103 kg. Satellite measurements show that this NLC display occurred during a burst of cloud activity that began on 5 August and lasted for 10 days. This cloud appeared 10 days after a launch of the space shuttle. We discuss the appearance of NLCs in August over several years at this lower polar latitude site in terms of planetary wave activity and space shuttle launches.  相似文献   

7.
In situ measurements of electron density were made over Trivandrum (8.5°N, 76.9°E) during nighttime to study E-region plasma density irregularities. Irregularities, with vertical scale sizes from a few km to 15 cm, were detected during rocket ascent and descent. Electron density profiles during ascent and descent of an earlier nighttime rocket flight from Trivandrum are also presented. Some of the important results are as follows: (i) horizontal gradients in electron density exist in 110–120 km region with horizontal scale size of at least 40 km, (ii) based on the presence/absence of electron density structures during ascent and descent of both flights, the horizontal distance over which the gradient drift instability operates is found to be at least 80 km and 90 km, for both the flights, (iii) observed irregularities in regions of negative density gradient are suggested to be produced through the gradient drift instability (GDI) driven by vertical polarization electric field as well as by electric field produced through wind shears and those in positive gradient regions by wind driven GDI, (iv) largest irregularity amplitude (≈30%) was associated with steepest gradients and so was the presence of smallest vertical scale sizes (12 m to 15 cm), which were absent at other altitudes, (v) the spectral index of irregularities was in the range of ?2.2±0.2 for large scales (few kilometers>λ>50 m), ?3.25±0.25 for medium scales (50 m>λ>10 m) and ?2.6±0.1 for smaller scales (10 m>λ>1 m) and (vi) irregularities in large and medium scales are expected to be produced directly through GDI and the small and sum-meter scales through non-linear GDI.  相似文献   

8.
The seasonal/annual characteristics of the high-altitude sporadic metal atom layers are presented on the basis of extensive Na and Fe lidar measurements at 30°N during the past several years. It is found that the extremely high sporadic Na (Nas) and Fe (Fes) layers above 105 km occurred mostly during summer. They had long durations (a few hours) and broad layer widths (much larger than 2 km). Their absolute peak densities could be comparable to or even larger than those of the corresponding main layers on a few nights. By using all the raw data profiles including sporadic layers, we have constructed the contour plots of Na and Fe densities versus month and altitude at 30°N. The Na and Fe layers both exhibit evidence for summer topside extension, which is consistent with the earlier observations for K and Ca at different latitudes. The summer topside extension of mean metal atom layers might represent a universal phenomenon that is alike for different atom species, different geographic locations and different measurement years. The extremely high sporadic metal atom layers above 105 km occurring during summer give rise to the phenomenon.  相似文献   

9.
The Michelson Interferometer for Passive Atmospheric Sounding (MIPAS), a limb sounding Fourier transform interferometer in the mid-infrared band, on board the polar-orbiting ENVIronmental SATellite (ENVISAT) was launched in March 2002 by the European Space Agency (ESA). For the MIPAS data validation, three balloon flights with MIPAS-B, a balloon version of MIPAS, were carried out on the 24th and 25th of September 2002 from Aire sur l’Adour (France 44°N, 0°E), on the 20th and 21th of March, and on the 2nd and 3rd of July 2003 from Esrange, Kiruna (Sweden 68°N, 21°E). The MIPAS operational data version 4.61 for the temperature vertical profiles were compared with the correlative MIPAS-B measurements via the coincident comparison and the trajectory comparison approaches, respectively. The precision of the MIPAS temperature was estimated to be 1.04–2.48 K in the region 123–8 hPa (15–33 km) and beyond the expectation. The systematic difference shows that the absolute accuracy of the MIPAS temperature is within the total combined errors between 356 and 3 hPa (8–39 km). The agreements in the middle stratosphere are better than in the lower stratosphere and upper troposphere. The maximum difference of 1.9 K (0.9%) was found at 228 hPa (11 km). The biases of MIPAS measurements are ?0.39 K (?0.19%) and +0.25 K (+0.12%) (“+” positive bias, “?” negative bias) corresponding to the altitude regions 143–31 hPa (14–24 km) and 31–5 hPa (24–37 km), respectively.  相似文献   

10.
Rayleigh lidar observations at Gadanki (13.5°N, 79.2°E) show an enhancement of the nightly mean temperature by 10–15 K at altitudes 70–80 km and of gravity wave potential energy at 60–70 km during the 2009 major stratospheric warming event. An enhanced quasi-16-day wave activity is observed at 50–70 km in the wavelet spectrum of TIMED–SABER temperatures, possibly due to the absence of a critical level in the low-latitude stratosphere because of less westward winds caused by this warming event. The observed low-latitude mesospheric warming could be due to wave breaking, as waves are damped at 80 km.  相似文献   

11.
We report observations of seasonal and local time variation of the averaged electron and iron concentrations, as well as simultaneous measurements of the two species, above the Arecibo Observatory (18.35°N, 66.75°N), Puerto Rico. The average Fe profile between 21:00 and 24:00 LT has a single peak at about 85 km with the exception of the summer when an additional peak exists at about 95 km. The higher Fe peak in the summer is correlated with higher electron concentrations in this season. The three nights of simultaneous measurements of electron and iron concentrations show that narrow layers of Fe and electrons are well correlated. Comparison of the climatological and simultaneous Fe and electron data suggests that recombination of Fe+ plays an important role in determining the Fe profile in the upper part of the Fe layer. Above 93 km, the Fe concentration appears to increase after sunset if the electron concentration exceeds about 4000 electrons cm−3. The average rate of Fe production is about 0.1 atom cm−3 s−1 for all seasons at 100 km in the early evening hours. A chemical model reveals that the concentration of Fe+ must be 50–80% of the total ionization over Arecibo for typical equinox conditions to explain the observed rate of Fe production. These high relative Fe+ concentrations are consistent with in situ observations that Fe+ is usually the dominant ion in sporadic E layers in the nighttime lower E region. This suggests that the source of Fe+ is provided by sporadic E layers descending over Arecibo after sunset. The Fe density between 80 and 85 km decreases during the night, for all seasons. This is attributed to the formation of stable molecular Fe species, such as FeOH, due to the increase in O3 and decrease in atomic O and H during the night at these altitudes.  相似文献   

12.
Explosion deep seismic sounding data sections of high quality had been obtained with RV Meteor in the Reykjanes Iceland Seismic Project (RRISP77 [Angenheister, G., Gebrande, H., Miller, H., Goldflam, P., Weigel, W., Jacoby, W.R., Pálmason, G., Björnsson, S., Einarsson, P., Pavlenkova, N.I., Zverev, S., Litvinenko, I.V., Loncarecic, B., Solomon, S., 1980. Reykjanes Ridge Iceland Seismic Experiment (RRISP 77). J. Geophys. 47, 228–238]) which close an information gap near 62°N. Preliminary results were presented by Weigel [Weigel, W., 1980. Aufbau des Reykjanes Rückens nach refraktionsseismischen Messungen. In: Weigel, W. (Ed.), Reykjanes Rücken, Island, Norwegischer Kontinentalrand. Abschlusskolloquium, Hamburg zur Meteor-Expedition, vol. 45. DFG, Bonn, pp. 53–61], and here we report on the data and results of interpretation. Clear refracted phases to 90 km distance permit crustal and uppermost mantle structure to be modelled by ray tracing. The apparent P-wave velocities are around 4.5, 6–6.5, 7–7.6 and 8.2–8.7 km/s, but no wide-angle reflections have been clearly seen. Accompanying sparker reflection data reveal thin sediment ponds in the axial zone and up to 400 m thick sediments at 10 Ma crustal age. Ray tracing reveals the following model below the sediments: (1) a distinct, 1–2 km thick upper crust (layer 2A) with Vp increasing with age (to 10 Ma) from <3.4 to 4.9 km/s and with a vertical gradient of 0.1–0.2 km/s/km, (2) a lower crust or layer 3 beginning at depths of 2 (axis) to 4 km (10 Ma age) below sea level with 6.1–6.8 km/s and similar vertical gradients as above, (3) the lower crust bottoms at 5.2–9.5 km depth below sea level (0–10 Ma) with a marked discontinuity, underneath which (4) Vp rises from about 7.5–7.8 km/s (0–10 Ma) with a positive vertical gradient of, again, 0.1–0.2 km/s/km such that 8 km/s would be reached at 12 km and deeper near the axis. Our preferred interpretation is that the mantle begins at the distinct discontinuity (“Moho”), but a deeper “Moho” of Vp ≈ 8 km/s cannot be excluded. From Iceland southward to 60°N several experiments show a decrease of crustal thickness from 14 to 8 km. Velocity trends with age across the ridge reflect cooling and filling of cracks, and thickness trends probably suggest volcanic productivity variations as previously suggested.Gravity inversion concentrates on a profile across the ridge with the above seismic a priori information; with 0.2–0.5 km depth uncertainty it leads to a good fit (±2.5 mGal where seismic data exist). Best fitting densities are (in kg/m3) for sediments, 2180; upper crust, 2450–2570; lower crust, 2850–2940; mantle lithosphere, 3215–3240 with a deficit for an asthenospheric wedge of no more than −100 kg/m3. The morphological ridges and troughs superimposed on the SE ridge flank are partly correlated, partly anti-correlated with the Bouguer anomaly and suggest that variable crustal density variations accompany the morphology variations.  相似文献   

13.
Long-term changes of the temperature of the middle atmosphere are investigated using a data bank obtained by Russian rocketsondes at Heiss Island (80.6°N, 58°E). The major interest of the data series is that it is one of the longest and uninterrupted records obtained at high latitudes in the northern hemisphere over 25 years, from 1969 to 1994. Previous estimates using this dataset has shown the largest trends. The revised analysis performed here took into account all possible discontinuities in the data series, such as a change in the time-of-measurement, T0, and in the type of sensor. For this purpose, some data were filtered out, and a statistical model based on multiple regression analyses included step functions to take into account such discontinuities. The temperature responses to different sources of variability (solar activity, volcanic aerosols) were retrieved for summer and winter periods. The response to the 11-year solar cycle in the winter period is found to be largely positive in the stratosphere (∼+4 K) and largely negative in the mesosphere (∼−8 K), with a smaller and opposite response in summer. This response depends on the phase of the QBO, as already shown by previous studies. The response to volcanic aerosols is found to be significantly positive in the upper mesosphere, in good agreement with numerical simulations and with observations above France. The long-term trend resulting from this reanalysis indicates a cooling of the middle atmosphere, increasing with altitude from −2 K/decade at 40 km to a maximum of −6 K/decade around 65 km. This result is slightly larger than the trend observed at mid-latitude but quite smaller than previous estimates.  相似文献   

14.
The existence of the F3-layer has been observed at Brazilian equatorial stations. This paper reports on a 1-year observations at Parit Raja, Malaysia (Lat. 1°52′N, Long. 103°48′E). The greatest number of appearances of the F3-layer is around local noon while the least is after local dawn. Its occurrence is more pronounced during winter and the equinoxes. The mean height of reflection for the layer is about 600 km reaching a maximum of about 900 km during the winter of 2004/2005. It is seen that the critical frequency of the F2-layer decreases with the appearance of the F3-layer.  相似文献   

15.
The temporal variability in currents, temperature, and particulate matter concentration were measured in the Mississippi Canyon axis where the thalweg was 300 m deep from May–July and August–November 1998 using current meters, thermographs, a light-scattering sensor, and sediment traps. Canyon sediments were sampled by coring and observed using an ROV video camera. Currents in the upper Mississippi Canyon generally oscillated up/down canyon with diurnal periodicity and were bottom-intensified. Mean current speed at 3.5 mab was approximately 8 cm s?1 during both deployments, reaching maximum speeds of over 50 cm s?1 under normal conditions. Based on current velocities, critical bed shear stress for resuspension of canyon-floor sediments was exceeded about 30% of the time during both deployments. In late September, Hurricane Georges passed 150 km NE of the study site, significantly intensifying current velocities, bed shear stress, resuspension, trap fluxes and temperature fluctuations. As the hurricane passed, maximum current speed reached 68 cm?s and temperature decreased ～7 °C in less than two hours. Critical bed shear stress for sediment resuspension was exceeded approximately 50% of the time during the five days of hurricane influence. Further evidence for sediment resuspension was the five-fold (and perhaps 70–130 fold) increase in trap fluxes and compositional similarities between canyon surface sediment and material collected by traps.  相似文献   

16.
The paper is focused on the global spatial structure, seasonal and interannual variability of the ～5-day Rossby (W1) and ～6-day Kelvin (E1) waves derived from the SABER/TIMED temperature measurements for 6 full years (January 2002–December 2007). The latitude structure of the ～5-day W1 wave is related to the gravest symmetric wave number 1 Rossby wave. The vertical structure of the ～5-day Rossby wave amplitude consists of double-peaked maxima centred at ～80–90 km and ～105–110 km. This wave has a vertically propagating phase structure from the stratosphere up to 120 km altitude with a mean vertical wavelength of ～50–60 km. The ～6-day E1 wave is an equatorially trapped wave symmetric about the equator and located between 20°N and 20°S. Its seasonal behaviour indicates some equinoctial and June solstice amplifications, while the vertical phase structure indicates that this is a vertically propagating wave between 20–100 km altitudes with a mean vertical wavelength of ～25 km.  相似文献   

17.
The preliminary results of observation of the lightning electrical characteristics during tropical summer thunderstorms locally known as “Nor-Wester”, at a hilly place in North-East India (23.50°N, 91.25°E) are reported here. Some distinct peaks are observed in the VLF range between 1.5 to 6 kHz. Average cloud conductivity is found to be 8.12×10?10 S m?1. Specific characteristics of integrated field intensity of sferics (IFIS) at 14 and 17 kHz are also studied. Average enhancement of electrical activity during thunderstorms is found to be 170 dB kHz?1 compared to the normal weather sferics intensity.  相似文献   

18.
We invert measurements of coseismic displacements from 139 continuously recorded GPS sites from the 2010, Jiashian, Taiwan earthquake to solve for fault geometry and slip distribution using an elastic uniform stress drop inversion. The earthquake occurred at a depth of ~ 23 km in an area between the Western Foothills fold-and-thrust belt and the crystalline high mountains of the Central Range, providing an opportunity to examine the deep fault structure under Taiwan. The inferred rupture plane is oblique to the prominent orientation of thrust faults and parallel to several previously recognized NW-striking transfer zones that appear to connect stepping thrusts. We find that a fault striking 318°–344° with dip of 26°–41° fits the observations well with oblique reverse-sinistral slip under a low stress drop of about 0.5 MPa. The derived geodetic moment of 2.92 × 1018 N-m is equivalent to a Mw = 6.24 earthquake. Coseismic slip is largely concentrated within a circular patch with a 10-km radius at the depth between 10 and 24 km and maximum slip of 190 mm. We suggest this earthquake ruptured the NW-striking Chishan transfer fault zone, which we interpret as a listric NE-dipping lateral ramp with oblique slip connecting stepping thrust faults (ramps). The inferred slip on the lateral ramp is considerably deeper than the 7–15 km deep detachment identified in previous studies of western Taiwan. We infer an active basal detachment under western Taiwan at a depth of at least ~ 20–23 km based on these inversion results. The earthquake may have nucleated at the base of the lateral ramp near the intersection with the basal detachment. Coulomb stress change calculations suggest that this earthquake moved several NE-striking active thrust faults in western Taiwan nearer to failure.  相似文献   

19.
We present the first results of gravity wave signatures on polar mesospheric clouds (PMCs) during the summer of 2007, in the northern hemisphere polar region. The Cloud Imaging and Particle Size (CIPS) experiment has one of the three instruments on board the NASA Aeronomy of Ice in the Mesosphere (AIM) spacecraft, which was launched into a sun-synchronous orbit on April 25, 2007. CIPS is a four-camera, wide-field (120°×80°) imager designed to measure PMC morphology and particle properties. One of the objectives of AIM is to investigate gravity wave effects on PMC formation and evolution. CIPS images show distinct wave patterns and structures in PMCs that are similar to ground-based photographs of noctilucent clouds (NLCs). The observed horizontal wavelengths of the waves were found to vary between 15 and 320 km, with smaller-wavelength structures of less than 50 km being the most common. In this paper we present examples of individual quasi-monochromatic wave events observed by CIPS and statistics on the wave patterns observed in the northern hemisphere during the summer months of 2007, together with a map showing the geographic locations of gravity wave events observed from CIPS.  相似文献   

20.
Continuous magnetotelluric (MT) measurements were conducted from May 2008 to July 2009 at Sakurajima, one of the most active volcanoes in Japan. Two observation sites were established at locations 3.3 km east and 3 km west–northwest of the summit crater. At both observation sites, the high-quality component of the impedance tensor (Zyx) showed variations in apparent resistivity of approximately ± 20% and phase change of ± 2°, which continued for 20–180 days in the frequency range between 320 and 4 Hz. The start of the period of changes in apparent resistivity approximately coincided with the start of uplift in the direction of the summit crater, as observed by a tiltmeter, which is one of the most reliable pieces of equipment with which to detect magma intrusion beneath a volcano. A 2D inversion of MT impedance suggests that the resistivity change occurred at a depth around sea level. One of the possible implications of the present finding is that the degassed volatiles migrated not only vertically through the conduit but also laterally through a fracture network, mixing with shallow groundwater beneath sea level and thereby causing the observed resistivity change.  相似文献   

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

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