首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 969 毫秒
1.
The radiosonde data available from British Atmospheric Data Centre (BADC) for the latitudinal occupancy of 58° north through 45° south were analyzed to observe the variation of temperature and water vapor density. These two climatological parameters are largely assumed to be the influencing factors in determining the millimeter wave window frequencies over the chosen range of latitudes in between the two successive maxima occurring at 60 and 120 GHz. It is observed that between temperature and water vapor density, the later one is influencing mostly in determining the window frequency. It is also observed that the minima is occurring at 75 GHz through 94 GHz over the globe during the month January–February and 73 GHz through 85 GHz during the month July–August, depending on the latitudinal occupancy. It is observed that the large abundance of water vapor is mainly held responsible for shifting of minima towards the low value of frequencies. Hence, it is becoming most important to look at the climatological parameters in determining the window frequency at the place of choice.  相似文献   

2.
Increasing demand of bandwidth in communication satellites has forced satellite links to be designed in Ku bands and above. But at these frequencies, rain and other tropospheric elements result in large attenuation. To mitigate the tropospheric attenuation of microwave satellite signals above 10 GHz using any standard Fade Mitigation Technique (FMT), it is essential to have a priori knowledge about the level of attenuation. Hence, short-term rain attenuation prediction models play a key role in maintaining the link in which necessary compensation can be applied depending on the early information of attenuation. This paper presents a method of attenuation prediction using Adaptive Artificial Neural Network. Here In situ Learning Algorithm (ILA) has been used to enable the system to track the non-stationary nature of the attenuation. To validate this, Ku Band data, collected at three different sites in India have been used for the purpose of prediction. The performance of the algorithm is determined through the estimation of prediction accuracy by comparing the predicted values with the measured data. Results obtained using the mentioned technique shows considerably good accuracy even up to 20 s of prediction interval with acceptable ratio between the under and over predictions. The prediction performance is evaluated for different prediction intervals. Furthermore the present model is also compared with the persistence model and the relative performance is quantified.  相似文献   

3.
The microstructure of rain has been studied with observations using a vertical looking Micro Rain Radar (MRR) at Ahmedabad (23.06°N, 72.62°E), a tropical location in the Indian region. The rain height, derived from the bright band signature of melting layer of radar reflectivity profile, is found to be variable between the heights 4600 m and 5200 m. The change in the nature of rain, classified on the basis of radar reflectivity, is also observed through the MRR. It has been found that there are three types of rain, namely, convective, mixed and stratiform rain, prevailing with different vertical rain microstructures, such as, Drop Size Distribution (DSD), mean drop size, rain rate, liquid water content and average fall speed of the drops at different heights. It is observed that the vertical DSD profile is more inhomogeneous for mixed and stratiform type rain than for convective type rain. It is also found that the large number of drops of size <0.5 mm is present in convective rain whereas in stratiform rain, drops concentration is appreciable up to 1 mm. A comparison of measurements taken by ground based Disdrometer and that from the 200 m level obtained from MRR shows good agreement for rain rate and DSD at smaller rain rate values. The results may be useful for understanding rain structures over this region.  相似文献   

4.
A lidar has been operated in São José dos Campos, Brazil (23.2°S, 45.8°W) since 1972, mainly dedicated to the study of mesospheric sodium at the 589 nm resonant line. The molecular Rayleigh scattering can also be used provided we limit the height to ∼75 km where the sodium scattering begins. Nevertheless, the weak signal obtained only permits the determination of density and temperature profiles by accumulating a large number of shots giving only nocturnal average profiles. Temporal variations in density and temperature on the scale of hours can however, be obtained by performing a superposed epoch analysis for a given time interval and covering a period of several days. In this way we obtained hourly mean profiles grouped by months, seasons and overall, with data acquired from 1993 to 2004. The difference between the hourly temperatures and the nocturnal means shows for some months, with enough data coverage, downward propagating structures that apparently have tidal origin. The seasonal averages show a recurrent feature with high temperatures before and low temperatures after midnight above 50 km. Some similarity is found with the GSWM model, but the observed temperature amplitudes are twice of that for the model.  相似文献   

5.
Laboratory characterization of dielectric properties of terrestrial analogues of lunar soil (JSC-1A) and comparison with lunar samples returned from various Apollo missions is made at different as well as normalized bulk density. Here measurements of dielectric constants and losses were made at four microwave frequencies such as 1.7 GHz, 2.5 GHz, 6.6 GHz and 31.6 GHz. Complex permittivity of lunar simulant was measured at temperature ranging from −190 °C to + 200 °C using Wave-Guide cell method. Comparison of permittivity of JSC-1A with Apollo sample also has been done at similar microwave frequencies. The investigations reveal that dielectric constant and loss factor of terrestrial analogues of lunar soil are temperature dependent. As temperature is gradually increased both these parameter (storage factor and loss factor) also gradually increases. These temperatures were chosen because the Moon undergoes at that extremes level of temperature. It is scorching heat at 110 °C during the day and freezing cold at −180 °C during night. The measured value of ε can be useful for designing passive as well as active sensors.  相似文献   

6.
Areas with dimensions of 1000–3000 km in which the total ozone content (TOC) decreases fast are called ozone mini-holes. They are generated mainly dynamically in two ways, either by poor-ozone air mass transport from the tropics to higher latitudes by planetary wave activity or, they are connected with strong adiabatic uplifting of the tropopause height. An ozone mini-hole, generated by the second mechanism, was observed over the Balkan Peninsula on 19/21 March 2005. In the middle of March, the polar vortex was strongly disturbed by Rossby waves, reaching up to the lower stratosphere. Warming episodes over a geographical area, covering the Barents Sea and the Polar Sea north from Central Siberia, displaced a polar vortex fragment extremely southwards. However, the vorticity was weak and the stratospheric temperatures did not reach low values, providing conditions for ozone chemical destruction via heterogenic reactions. At the same time, a Rossby wave ridge was located below the European polar fragment. In the period from 13 to 19 March, the thermal tropopause over Sofia was uplifted almost by 3 km. Ozone distributions observed by the SCIAMACHY instrument on 18–21 March show a fast TOC decrease westwards from Ireland, which was moving eastwards during the next days, increasing the area in which the ozone content decreased. On 20/21 March low ozone content was observed above the Stara Zagora (42°N, 25°E) ground-based station by means of the GASCOD instrument, using DOAS technique. The TOMS Earth probe instrument detected 237 DU over Sofia. This is a record low March value from the beginning of the TOMS instrument measurements in 1978. In March/April the ozone distribution was characterized by its mean annual maximum of 360 DU at 42°N.  相似文献   

7.
The effects of a major stratospheric sudden warming (SSW) at extratropical latitudes have been investigated with wind and temperature observations over a Brazilian station, Cachoeira Paulista (22.7°S, 45°W) during September–October 2002. In response to the warming at polar latitudes a corresponding cooling at tropical and extratropical latitudes is prominent in the stratosphere. A conspicuous signature of latitudinal propagation of a planetary wave of zonal wavenumbers 1 and 2 from polar to low latitude has been observed during the warming period. The polar vortex which split into two parts of different size is found to travel considerably low latitude. Significant air mass mixing between low and high latitudes is caused by planetary wave breaking. The meridional wind exhibits oscillations of period 2–4 days during the warming period in the stratosphere. No wave feature is evident in the mesosphere during the warming period, although a 12–14 day periodicity is observed after 2 weeks of the warming event, indicating close resemblance to the results of other simultaneous investigations carried out from high latitude Antarctic stations. Convective activity over the present extratropical station diminishes remarkably during the warming period. This behavior is possibly due to destabilization and shift of equatorial convective active regions towards the opposite hemisphere in response to changes in the mean meridional circulation in concert with the SSW.  相似文献   

8.
A sodium lidar, capable of measuring temperature in the 80–100 km region, has been in operation at São José dos Campos (23° S, 46 W) since March 2007. Good quality data have been obtained for late autumn, winter and spring, but weather conditions make it extremely difficult to make measurements from mid-November to mid- February. We find the temperature structure to be strongly modulated by tides and gravity waves, but average profiles typically show a primary mesopause height close to 100 km with temperatures around 180 K, and a tendency for a secondary minimum of about 185 K to occur close to 90 km. Vertical temperature gradients greater than 50 K/km are sometimes seen even on profiles averaged over several hours. The strongest gradients are always positive and are frequently associated with strong gradients in sodium concentration. On the other hand, we frequently see rapid changes in the temperature profile, suggesting that models and non-local temperature measurements, as made by satellite radiometers, for example, are of little use in applications such as the analysis of gravity wave propagation seen in airglow images.  相似文献   

9.
This paper presents an investigation into the variability and predictability of the maximum height of the ionospheric F2 layer, hmF2 over the South African region. Data from three South African stations, namely Madimbo (22.4°S, 26.5°E, dip angle: −61.47°), Grahamstown (33.3°S, 26.5°E, dip angle: −64.08°) and Louisvale (28.5°S, 21.2°E, dip angle: −65.44°) were used in this study. The results indicate that hmF2 shows a larger variability around midnight than during the daytime for all seasons. Monthly median hmF2 values were used in all cases and were compared with predictions from the IRI-2007 model, using the URSI (Union Radio-Scientifique Internationale) coefficient option. The analysis covers the diurnal and seasonal hourly hmF2 values for the selected months and time sectors e.g. January, July, April and October for 2003 and 2005. The time ranges between (03h00–23h00 UT; LT = UT + 2h) representing the local sunrise, midday, sunset and midnight hours. The time covers sunrise, midday, sunrise, and midnight hours (03–06h00 UT, 07–11h00 UT, sunrise 16–18h00 UT and 22–23h00 UT; LT = UT + 2h). The dependence of the results on solar activity levels was also investigated. The IRI-2007 predictions follow fairly well the diurnal and seasonal variation patterns of the observed hmF2 values at all the stations. However, the IRI-2007 model overestimates and underestimates the hmF2 value during different months for all the solar activity periods.  相似文献   

10.
The occurrence of mid-latitude spread F (SF) over South Africa has not been extensively studied since the installation of the DPS-4 digisondes in 1996 and 2000 at Grahamstown (33.32 °S, 26.50 °E) and Madimbo (22.38 °S, 30.88 °E) respectively. This study is intended to quantify the probability of occurrence of F region disturbances associated with SF over South Africa. A study was conducted using data for 8 years (2001–2008) over Madimbo (with a time resolution of 30 min) and Grahamstown (with a variable time resolution of 15 and 30 min). In this study, SF has been classified into frequency SF (FSF), range SF (RSF) and mixed SF (MSF). The SF events were identified by manually identifying ionograms showing SF and tabulating them according to type for further statistical analysis. The results show that the diurnal pattern of SF peaks strongly between 01:00 and 02:00 local time, LT (LT = UT + 2 h), where UT is the universal time. This pattern is true for all seasons and types of SF at Madimbo and Grahamstown in 2001 and 2005, except for RSF which had peaks during autumn and spring in 2001 at Madimbo. The probability of both MSF and FSF tends to increase with decreasing solar activity, with a peak in 2005 (a moderate solar activity period). The seasonal peaks of MSF and FSF are more frequent during winter months at both Madimbo and Grahamstown. In this study, SF was evident in ∼0.03% and ∼0.06% of the available ionograms at Madimbo and Grahamstown respectively during the 8 years.  相似文献   

11.
The Indian MST radar facility at Gadanki (13.5°N, 79.2°E) has been utilised to study the propagation of gravity waves from the troposphere/lower stratosphere to the mesosphere and their interaction with the radar backscattered signal variations. The main objective is to correlate vertically propagating gravity waves derived from the tropospheric velocity fields with the dynamics of mesospheric scattering centres. The tropospheric wind velocities and signal strengths over the entire height range have been subjected to power spectral and wavelet analysis to determine the predominant wave periods/amplitudes and the coupling between the lower atmosphere and mesosphere. Results show that (a) the gravity waves are clearly detectable near tropopause heights, (b) while relatively higher period gravity waves (20–50 min) interact with mesospheric scattering centres, the lower period waves (<20 min) are absorbed in the troposphere itself, (c) the mesospheric scattering layers are affected by gravity waves of complementary periods.  相似文献   

12.
The analysis of energetic particles and magnetic field measurements from the Ulysses spacecraft has shown that in a series of events, the energy density contained in the suprathermal tail particle distribution is comparable to or larger than that of the magnetic field, creating conditions of high-beta plasma. In this work we analyze periods of high-beta suprathermal plasma occurrences (βep > 1) in interplanetary space, using the ratio (βep) of the energetic particle (20 keV to ∼5 MeV) and magnetic field energy densities from measurements covering the entire Ulysses mission lifetime (1990–2009) in order to reveal new or to reconfirm some recently defined interesting characteristics. The main key-results of the work are summarized as follows: (i) we verify that high-beta events are detected within well identified regions corresponding mainly to the vicinity of shock surfaces and magnetic structures, and associated with energetic particle intensity enhancements due to (a) reacceleration at shock-fronts and (b) unusually large magnetic field depressions. (ii) We define three considerable features for the high-beta events, concentrated on the next points: (a) there is an appreciable solar-activity influence on the high-beta events, during the maximum and middle solar-cycle phase, (b) the annual peak magnitude and the number of occurrences of high events are well correlated with the sunspot number, (c) the high-beta suprathermal plasma events present a spatial distribution in heliographic latitudes (HL) up to ∼±80°, and a specific important concentration on the low (−25° ? HL < −6°, 6° < HL ? 25°) and median (−45° ? HL < −25°, 25° < HL ? 45°) latitudes. We also reconfirm by a statistical analysis the results of Marhavilas and Sarris (2011), that the high-beta suprathermal plasma (βep > 1) events are characterized by a very large parameter βep (up to 1732.5), a great total duration (406 days) and a large percentage of the Ulysses-mission lifetime (which is equal to 6.34% of the total duration with usable measurements, and 11.3% of the duration in presence of suprathermal particles events).  相似文献   

13.
The nighttime vertical E × B drifts velocities of the F2-region were inferred from the hourly hmF2 values obtained from ionosonde data over an African equatorial station, Ilorin (8.50oN, 4.68oE; dip lat. 2.95o) during period of low solar activity. For each season, the plasma drift Vz is characterized by an evening upward enhancement, then by a downward reversal at 1900 LT till around 0000 LT, except for June solstice. This was explained using the Rayleigh–Taylor (R-T) instability mechanism. The occasional drift differences in Vz obtained by inferred and direct measurement over Ilorin and Jicamarca, respectively are reflective of the importance of chemistry and divergent transport system due to both the E region electric and magnetic fields instead of simple motions. The pre-reversal enhancement (PRE) magnitude is higher during the equinoctial months than the solsticial months over Jicamarca, highest during December solstice and the equinoctial months over Ilorin, suggesting the dominance of higher E × B fountain during equinoxes at both stations. The lowest PRE magnitude was in June solstice. The appearance of post-noon peak in NmF2 around 1700 LT is highest during the equinoctial months and lowest during the solsticial period. A general sharp drop in NmF2 around 1800 LT is distinct immediately after sunset, lowest during June solstice and highest in March equinox. Our result suggests that between 0930 and 2100 LT, the general theory that vertical drifts obtained by digisonde measurements only match the E × B drift if the F layer is higher than 300 km is reliable, but does not hold for the nighttime period of 2200–0600 LT under condition of solar minima. Hence, the condition may not be sufficient for the representation of vertical plasma drift at nighttime during solar minima. This assertion may still be tentative, as more equatorial stations needed to be studied for better confirmation.  相似文献   

14.
To investigate the Mesosphere and Lower Thermosphere (MLT) region, several ground-based instruments called SATI (Spectral Airglow Temperature Imager) were designed and built to measure airglow emission and temperature in the upper mesosphere. One SATI instrument was installed at Resolute Bay (74.7°N, 94.9°W) and has monitored the polar MLT region since November, 2001. In October 2007 another SATI instrument was installed at Eureka (80.0°N, 86.3°W) at the Polar Environment Atmospheric Research Laboratory (PEARL) as part of the Canadian Network for the Detection of Atmospheric Change (CANDAC) project. SATI is a spatial scanning Fabry–Perot spectrometer measuring column emission rates for several rotational lines of OH and O2 airglow at 87 and 94 km height. The rotational temperatures are inferred from the ratios of these lines. The measurements are divided into 12 sectors with an annular field of view. The phase differences between the sectors yield information on the horizontal atmospheric wave direction and wavelength. Horizontal perturbations of 2–8 h period have correlatively been observed and investigated at both locations. Short-periodic oscillations identified as gravity waves with periods between 2 and 8 h propagate in southward and eastward directions, but in opposite directions in some cases. The wave propagation characteristics are often different at the two locations; the relationship with the lower mean wind is considered.  相似文献   

15.
A dual dipole antenna has been installed at low latitude station Kolhapur (Geographic 16.8°N, 74.25°E), Maharashtra, India for the study of cosmic radio noise absorption using Solid State Riometer (which operates at 30 MHz) during pre phase of 24th solar maxima. The aim for this type of study over Kolhapur was to know the response of lower (D region) ionosphere over low latitude by cosmic radio noise absorption using riometer technique during quite period as well as sudden ionospheric disturbances (SID). The observations are being taken for 3 years. Two different sites (∼40 km away from each other) were used for the installation of riometer equipment assuming minimum local noise. It is found that solar noise to cosmic radio noise hence resulting in signal saturation. The night time signal is relatively free of interference but sometimes local noise is responsible for spike-like signatures. Hence it is concluded that Kolhapur (a low latitude station) is not suitable for the study of cosmic radio noise absorption on 30 MHz with riometer and dual dipole antenna. Proper choice for operating frequency of riometer and antenna gain is suggested for low latitude use of this technique for ionospheric deviative and nondeviative absorption studies.  相似文献   

16.
It is shown in this paper for the first time that the intensity of the daytime thermospheric O(1D) 630.0 nm airglow as measured by the ground-based dayglow photometer over Trivandrum (8.5°N; 77°E; dip lat. 0.5°N), a geomagnetic dip equatorial station, exhibit a direct correlation with the electron density at 180 km. This altitude is about ∼40 km lower than the believed centroid of the O(1D) 630.0 nm dayglow emission i.e. 220 km. This observation is contrary to the understanding of the behavior of O(1D) 630.0 nm dayglow over equatorial/low latitudes. Over these latitudes, the variations of the measured intensity of O(1D) 630.0 nm dayglow are known to be associated with the changes in the electron density at altitudes around 220 km, the centroid of this emission. In this context, the present results indicating the lowering of the peak altitude of O(1D) 630.0 nm emission from ∼220 to ∼180 km over the dip equator is new. Recent results on solar XUV flux indicate that this could be an important parameter that controls the O(1D) 630.0 nm dayglow excitation rates through modulations in the neutral and ionic composition in lower thermosphere-ionosphere region. However, the lowering of the centroid of O(1D) 630.0 nm emission, as shown in this study, has been ascribed primarily to the fountain effect associated with the equatorial ionization anomaly.  相似文献   

17.
The Incoherent Scatter Radar measurement over Jicamarca, together with the IRI model-2007 measurements were compared with ground-based digisonde inferred E × B drift over Ilorin in the African region during year of solar minima (F10.7 = 81). Seasonally, Ilorin pre-reversal enhancement (PRE) had peak drift velocities of 7.2, 3.7 and 7.9 m/s for March equinox, September equinox and December solstice respectively, while Jicamarca drifts indicated 13.0, 10.5 and 5.2 m/s; as well as the IRI model with 14.3, 8.4 and 0.7 m/s in similar order. PRE value was insignificant during June solstice. The PRE magnitude of the IRI-model during the equinoxes is twice the value obtained at Ilorin. The daytime E × B drift peaked over Ilorin 1–2 h earlier than both the modeled and Jicamarca observations. This could be due to the difference in sunset time at the conjugate points corresponding to the altitude of the observation. During the evening time PRE, the respective correlation coefficients (R) for Vz–F10.7 relation over Jicamarca, Ilorin and the modeled observations are −0.5559, 0.4796 and −0.4979. Similarly, the Vz–Ap relation exhibit excellent anti-correlation coefficient (R = −0.8637) for the IRI-model, −0.4827 over Jicamarca and 0.3479 for Ilorin. Annual mean drift velocities over Jicamarca, Ilorin and IRI model measurements respectively are 10, 5.6 and 10 m/s for the peak PRE observation; 15, 16 and 21 m/s for the daytime pre-sunrise peak values; and −21, −9 and −16 m/s for the nighttime downward reversals. The root-mean square (RMS) deviation between IRI-model and the Ilorin drift between 2000 and 0500 h is 4.37, 2.03, 3.71 and 2.42 m/s for March equinox, June solstice, September equinox and December solstice respectively. For Jicamarca–Ilorin drift relation, RMS deviation is 5.48, 2.30, 3.47 and 1.27 m/s in the same order respectively. Annual hmF2 inferred drift over Ilorin during daytime is higher by a factor of ≈2 and 3 at Jicamarca and IRI model measurements respectively; and by a factor of ≈5 for both during the night-time period. The limitations in using hmF2 to infer drifts are discussed.  相似文献   

18.
Global Positioning System (GPS) receiver on the CHAllenging Mini-satellite Payload (CHAMP) and the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument, one of four on board the TIMED satellite, provide middle atmosphere temperature profiles by Radio Occultation (RO) and limb viewing infrared emission measurements, respectively. These temperature profiles retrieved by two different techniques in the stratosphere are compared with each other using more than 1300 correlative profiles in March, September and December 2005. The over-all mean differences averaged over 15 and 35 km are approximately −2 K and standard deviation is less than 3 K. Below 20 km of altitude, relatively small mean temperature differences ∼1 K are observed in wide latitudinal range except for June (during the SABER nighttime observation). In the middle to low latitudes, between 30°S and 30°N, the temperature difference increases with height from ∼0–1 K at 15 km, to ∼−4 K at 35 km of altitude. Large temperature differences about −4 to −6 K are observed between 60°S and 30°N and 31–35 km of altitude for all months and between 0° and 30°N below 16 km during June (nighttime).  相似文献   

19.
A new and original stereo imaging method is introduced to measure the altitude of the OH nightglow layer and provide a 3D perspective map of the altitude of the layer centroid. Near-IR photographs of the OH layer are taken at two sites separated by a 645 km distance. Each photograph is processed in order to provide a satellite view of the layer. When superposed, the two views present a common diamond-shaped area. Pairs of matched points that correspond to a physical emissive point in the common area are identified in calculating a normalized cross-correlation coefficient (NCC). This method is suitable for obtaining 3D representations in the case of low-contrast objects. An observational campaign was conducted in July 2006 in Peru. The images were taken simultaneously at Cerro Cosmos (12°09′08.2″ S, 75°33′49.3″ W, altitude 4630 m) close to Huancayo and Cerro Verde Tellolo (16°33′17.6″ S, 71°39′59.4″ W, altitude 2272 m) close to Arequipa. 3D maps of the layer surface were retrieved and compared with pseudo-relief intensity maps of the same region. The mean altitude of the emission barycenter is located at 86.3 km on July 26. Comparable relief wavy features appear in the 3D and intensity maps. It is shown that the vertical amplitude of the wave system varies as exp (Δz/2H) within the altitude range Δz = 83.5–88.0 km, H being the scale height. The oscillatory kinetic energy at the altitude of the OH layer is comprised between 3 × 10−4 and 5.4 × 10−4 J/m3, which is 2–3 times smaller than the values derived from partial radio wave at 52°N latitude.  相似文献   

20.
Using the GPS ionospheric scintillation data at Hainan station (19.5°N, 109.1°E) in the eastern Asia equatorial regions and relevant ionospheric and geomagnetic data from July 2003 to June 2005, we investigate the response of L-band ionospheric scintillation activity over this region to different strong magnetic storm conditions (Dst < −100 nT) during the descending phase of the solar cycle. These strong storms and corresponding scintillations mainly took place in winter and summer seasons. When the main phase developed rapidly and reached the maximum near 20–21 LT (LT = UT + 8) after sunset, scintillations might occur in the following recovery phase. When the main phase maximum occurred shortly after midnight near 01–02 LT, following the strong scintillations in the pre-midnight main phase, scintillations might also occur in the post-midnight recovery phase. When the main phase maximum took place after 03 LT to the early morning hours no any scintillation could be observed in the latter of the night. Moreover, when the main phase maximum occurred during the daytime hours, scintillations could also hardly be observed in the following nighttime recovery phase, which might last until the end of recovery phase. Occasionally, scintillations also took place in the initial phase of the storm. During those scintillations associated with the nighttime magnetic storms, the height of F layer base (h’F) was evidently increased. However, the increase of F layer base height does not always cause the occurrence of scintillations, which indicates the complex interaction of various disturbance processes in ionosphere and thermosphere systems during the storms.  相似文献   

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

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