首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The presence and movement of plasma density fluctuations in the F-region of the ionosphere are studied by monitoring phase and amplitude of radio waves propagating through the region. In this paper, we have used weak scattering theory and assumed the plasma density fluctuations to behave like phase changing diffraction screen. Appropriate relations for scintillation index S4, and phase variance δ? are derived and computed for different parameters of the plasma density irregularities of the ionosphere. SROSS-C2 satellite in situ measurements of plasma density fluctuations, which provide direct information about the structure and morphology of irregularities that are responsible for scintillation of radio waves, were used first time to develop a scintillation model for low latitude. It is observed that the scintillation index S4 and phase variance δ? depends on the strength of the plasma turbulence. Finally, the results obtained from modeling are compared and discussed with the available recent results.  相似文献   

2.
The ionospheric scintillation and TEC (Total Electron Content) variations are studied using GPS (Global Positioning System) measurements at an Indian low latitude station Surat (21.16°N, 72.78°E; Geomagnetic: 12.90°N, 147.35°E), situated near the northern crest of the equatorial anomaly region. The results are presented for data collected during the initial phase of current rising solar activity (low to moderate solar activity) period between January 2009 and December 2011. The results show that within a total number of 656 night-time scintillation events, 340 events are observed with TEC depletions, Rate of change of TEC (ROT) fluctuations and enhancement of Rate of change of TEC Index (ROTI). A comparison of night-time scintillation events from the considered period reveal strong correlation amongst the duration of scintillation activity in S4 index, TEC depletion, ROT fluctuations and ROTI enhancement in the year 2011, followed by the year 2010 and least in 2009. The statistical analyses of scintillation activity with enhancement of ROTI also show that about 70–96% scintillation activity took place in equinox and winter months. Moreover, from a nocturnal variation in occurrence of scintillation with (S4 ? 0.2) and enhancement of ROTI with (ROTI ? 0.5), a general trend of higher occurrence in pre-midnight hours of equinox and winter seasons is observed in both indices during the year 2011 and 2010, while no significant trend is observed in the year 2009. The results suggest the presence of F-region ionospheric irregularities with scale sizes of few kilometers and few hundred meters over Surat and are found to be influenced by solar and magnetic activity.  相似文献   

3.
This study characterizes equatorial scintillations at L-band frequency over Lagos, Nigeria during the minimum and ascending phases of solar cycle 24. Three years (2009–2011) of amplitude scintillation data were used for the investigation. The data were grouped on daily, monthly, seasonal, and yearly scales at three levels of scintillation (weak (0.3 ? S4 < 0.4), moderate (0.4 ? S4 < 0.7), and intense (S4 ? 0.7)). To ensure reliable statistical inferences, three data cut-off criteria were adopted. Scintillations were observed to have a daily trend of occurrence during the hours of 1900–0200 LT, and higher levels of scintillations were localized within the hours of 2000–2300 LT. On monthly basis, September and October recorded the highest occurrences of scintillation, while January recorded the least. Scintillations were recorded during all the months of 2011, except January. Surprisingly, pockets of scintillation events (weak levels) were also observed during the summer months (May, June, and July). Seasonally, equinoxes recorded the highest occurrences of scintillation, while June solstice recorded the least occurrences. Scintillation activity also increases with solar and geomagnetic activity. On a scintillation active day, the number of satellites available to the receiver’s view reduces as the duration of observation reduces. These results may support the development of future models that could provide real-time predictability of African equatorial scintillations, with a view to supporting the implementation of GNSS-based navigation for aviation applications in Africa.  相似文献   

4.
This study presents a statistical analysis of GPS L-band scintillation with data observed from July 2008 to March 2012 at the northern crest of equatorial anomaly stations in Guangzhou and Shenzhen of South China. The variations of the scintillation with local time, season, solar activity and duration of scintillation patches were investigated. The relationship between the scintillation and TEC depletion was also reported. Our results revealed that GPS scintillation occurred from 19:30 LT (pre-midnight) to 03:00 LT (post-midnight). During quiet solar activity years, the scintillation was only observed in pre-midnight hours of equinox months and patches durations were mostly less than 60 min. During high solar activity years, more scintillation occurred in the pre-midnight hours of equinox and winter months; and GPS scintillation started to occur in the post-midnight hours of summer and winter. The duration of scintillation patches extended to 180 min in high solar activity years. Solar activity had a larger effect to strong scintillations (S4 > 0.6) than to weak scintillations (0.6 ? S4 > 0.2). Strong scintillations were accompanied by TEC depletion especially in equinox months. We also discussed the relationship between TEC depletion and plasma bubble.  相似文献   

5.
The occurrence of ionospheric scintillations with S4 ? 0.2 was studied using GPS measurements at Guilin, China (25.29°N, 110.33°E; geomagnetic: 15.04°N, 181.98°E), a station located near the northern crest of the equatorial anomaly. The results are presented for data collected from January 2009 to March 2010. The results show that nighttime amplitude scintillations only took place in February and March of the considered years, while daytime amplitude scintillations occurred in August and December of 2009. Nighttime amplitude scintillations, observed in the south of Guilin, always occurred with phase scintillations, TEC (Total Electron Content) depletions, and ROT (Rate Of change of TEC) fluctuations. However, TEC depletions and ROT fluctuations were weak during daytime amplitude scintillations, and daytime amplitude scintillations always took place simultaneously for most of the GPS satellites which appeared over Guilin in different azimuth directions. Ground-based GPS scintillation/TEC observations recorded at Guilin and signal-to-noise-ratio (SNR) measurements obtained from GPS-COSMIC radio occultation indicate that nighttime and daytime scintillations are very likely caused by ionospheric F region irregularities and sporadic E, respectively. Moreover, strong daytime amplitude scintillations may be associated with the plasma density enhancements in ionospheric E region caused by the Perseid and Geminid meteor shower activities.  相似文献   

6.
The occurrence of ionospheric irregularities at high latitudes, with dimensions of several kms down to decameter scale size shows strong correlation with geomagnetic disturbance, season and solar activity. Transionospheric radio waves propagating through these irregularities experience rapid random fluctuations in phase and/or amplitude of the signal at the receiver, termed scintillation, which can degrade GNSS services. Thus, investigation and prediction of this scintillation effect is very important. To investigate such scintillation effects, a GISTM (GPS Ionospheric Scintillation and TEC Monitoring) NovAtel dual frequency (L1/L2) GPS receiver has been installed at Trondheim, Norway (63.41°63.41° N, 10.4°10.4° E), capable of collecting scintillation indices at a 1 min rate as well as the raw data (phase and intensity) of the satellite signals at a 50 Hz sampling rate and TEC (Total Electron Content) at a 1 Hz rate. Many researchers have reported that both phase and amplitude scintillation is closely associated with TEC fluctuations or associated with a significant developing enhancement or depletion in the TEC. In this study, a novel analogous phase index is developed which provides samples at a 1 min rate. Generally the scintillation indices can help in estimating the irregularity scintillation effect at a one minute rate, but such procedures are time consuming if DFTs of the phase and/or amplitude at a 50 Hz data are required. In this study, instead, this analogous phase index is estimated from 1 Hz rate TEC values obtained from the raw signals and is then compared for weak, moderate and strong scintillation at Trondheim for one year of data collected from the installed GPS receiver. The spectral index of the irregularities (that is the inverse power law of their spatial spectrum) is determined from the resultant phase scintillation psd. The correlations of the scintillation indices and spectral indices with the analogous phase index have been investigated under different geomagnetic conditions (represented by the Kp index) and an approximate linear correlation of phase scintillation with the analogous phase index was found. Then a principal advantage of this index is that it achieves this correlation without requiring a high sampling data rate and the need for DFTs. Thus, the index seems a good candidate for developing a simple means of ionospheric scintillation prediction which could also be utilized in the development of alerts using regional mappings.  相似文献   

7.
This study presents results on the investigation of the diurnal, monthly and seasonal variability of Total Electron Content (TEC), phase (σΦσΦ) and amplitude (S4) scintillation indices over Ugandan (Low latitude) region. Scintillation Network Decision Aid (SCINDA) data was obtained from Makerere (0.34°N, 32.57°E) station, Uganda for two years (2011 and 2012). Data from two dual frequency GPS receivers at Mbarara (0.60°S, 30.74°E) and Entebbe (0.04°N, 32.44°E) was used to study TEC climatology during the same period of scintillation study. The results show that peak TEC values were recorded during the months of October–November, and the lowest values during the months of July–August. The diurnal peak of TEC occurs between 10:00 and 14:00 UT hours. Seasonally, the ascending and descending phases of TEC were observed during the equinoxes (March and September) and solstice (June and December), respectively. The scintillations observed during the study were classified as weak (0.1≤S4,σΦσΦ0.3) and strong (0.3<<S4,σΦσΦ1.0). The diurnal scintillation pattern showed peaks between 17:00 and 22:00 UT hour, while the seasonal pattern follows the TEC pattern mentioned above. Amplitude scintillation was more dominant than phase scintillation during the two years of the study. Scintillation peaks occur during the months of March–April and September–October, while the least scintillations occur during the months of June–July. Therefore, the contribution of this study is filling the gap in the current documentation of amplitude scintillation without phase scintillation over the Ugandan region. The scintillations observed have been attributed to wave-like structures which have periods of about 2–3 h, in the range of that of large scale travelling ionospheric disturbances (LSTIDs).  相似文献   

8.
Rain drop size distribution (DSD) was measured at four places in Southern India {Thiruvananthapuram, Kochi, Munnar and Sriharikota (SHAR)} using a Joss–Waldvogel (JW) impact type disdrometer. The data for each minute were corrected for dead time errors and rain rate was computed from the corrected data. The data for a whole month were then sorted according to rain rate (R) into several classes ranging from 0.1 to >100 mm/h. The average DSD in each class was computed, and the lognormal distribution function was fitted to the average. In all the cases, the function fitted the data very well. The fit parameters were found to have dependence on rain rate. The total number of drops (NT), the geometric mean diameter (Dg) and the standard geometric deviation (σ) were also computed from the fit parameters. The standard geometric deviation (σ) was found to be more or less constant with rain rate at all the sites and in all months. The other two parameters (NT and Dg) were found to vary exponentially with rain rate except in Munnar, a high altitude station. At Thiruvananthapuram, in most of the months, NT increased exponentially with rain rate up to some value of R, which was different in different months, and then remained more or less constant or decrease slightly. In all cases, the variation of NT and Dg was such that NTDg3 increased linearly with rain rate.  相似文献   

9.
We have analysed a sample of 328 time-integrated GRB prompt emission spectra taken via the Konus instrument on board the US GGS-Wind spacecraft between 2002 and 2004 using a couple of two-components models, Cut-off Power Law (CPL) + Power Law (PL) and blackbody (BB) + PL. The spectra show clear deviation from the Band function. The PL term is interpreted as the low energy tail of a nonthermal emission mechanism. The distributions of corresponding index β give values β < −2/3 consistent with synchrotron and synchrotron self-Compton mechanisms. The distribution of low energy index α associated with the CPL term shows clear discordance with synchrotron models for 31.4% of the analysed GRBs with values exceeding that for the line of death, α = −2/3. Then, a set of nonthermal radiation mechanisms producing harder slopes, i.e., α > −2/3, are presented and discussed. For the remaining majority (68.6%) of GRBs with CPL index α < −2/3, we show that optically thin synchrotron produced by a power law electron distribution of type, N(γ) ∼ γp, γ1 < γ < γ2, for finite energy range (γ2 ≠ ∞) is a likely emission mechanism with α ∼−(p + 1)/2 in the frequency range ν1 ? ν ? ν2 (where ν2 = η2ν1 with η = γ2/γ1), such that for p > 1/3, one gets α < −2/3. We also show that corresponding spectra in terms of Fν and νFν functions are peaked around frequency ν2 instead of ν1, respectively for p < 1 and p < 3. Besides, thermal emission is examined taking a single Planck function for fitting the low energy range. It can be interpreted as an early emission from the GRB fireball photosphere with observed mean temperature, kT′ ∼ 16.8 keV. Furthermore, we have performed a statistical comparison between the CPL + PL and BB + PL models finding comparable χ2-values for an important fraction of GRBs, which makes it difficult to distinguish which model and specific radiation mechanism (possible thermal or nonthermal γ-ray emissions) are best suitable for describing the reported data. Therefore, additional information for those bursts, such as γ-ray polarization, would be highly desirable in future determinations of GRBs observational data.  相似文献   

10.
This work is a continuation of the previous article and it focuses on low solar activity and modeling effort. NeQuick model uses Epstein layer formalism to model each part of the profile. We study the diurnal and seasonal variations of B2bot, ΔB2 (B2best − B2NeQuick2) and R (B2best/B2NeQuick 2) at Hainan station during low solar activity. The results show it is possible to improve the B2bot parameter of the NeQuick model at that region during low solar activity. Then, we use a function ?(t) with LT in different seasons to correct the B2bot formula of NeQuick 2. The correction shows that (1) By the correction formula, the B2bot of NeQuick is improved. The maximum standard deviation is improved for 9 km. (2) The correction formula is more effective in summer than in equinox and winter and performs better during early morning hours than during the rest of the day.  相似文献   

11.
Multi-channel Global Positioning System (GPS) carrier phase signals, received by the six low Earth orbiting (LEO) satellites from the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) program, were used to undertake active limb sounding of the Earth’s atmosphere and ionosphere via radio occultation. In the ionospheric radio occultation (IRO) data processing, the standard Shell inversion technique (SIT), transformed from the traditional Abel inversion technique (AIT), is widely used, and can retrieve good electron density profiles. In this paper, an alternative SIT method is proposed. The comparison between different inversion techniques will be discussed, taking advantage of the availability of COSMIC datasets. Moreover, the occultation results obtained from the SIT and alternative SIT at 500 km and 800 km, are compared with ionosonde measurements. The electron densities from the alternative SIT show excellent consistency to those from the SIT, with strong correlations over 0.996 and 0.999 at altitudes of 500 km and 800 km, respectively, and the peak electron densities (NmF2) from the alternative SIT are equivalent to the SIT, with 0.839 vs. 0.844, and 0.907 vs. 0.909 correlation coefficients when comparing to those by the ionosondes. These results show that: (1) the NmF2 and hmF2 retrieved from the SIT and alternative SIT are highly consistent, and in a good agreement with those measured by ionosondes, (2) no matter which inversion technique is used, the occultation results at the higher orbits (∼800 km) are better than those at the lower orbits (∼500 km).  相似文献   

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.
In this paper we present a new mechanism of the main energy conversion of the solar flare. Since a flare inducing prominence (flux tube) rises Vz ? 300 km s−1, the plasmas below it cannot continuously eject with Alfvén speeds of VA = 3000 km s−1 but probably with Vz ≈ ±100 km s−1. Plasma up and downflows with VA will within a short duration be blocked between the chromosphere where reconnected flux tubes are piling up, and the slowly rising flux rope. Hence the Petschek slow shock mechanism is difficult to be realized as a major energy converting mechanism.  相似文献   

14.
First comparison of in situ density fluctuations measured by the DEMETER satellite with ground based GPS receiver measurements at the equatorial anomaly station Bhopal (geographic coordinates (23.2°N, 77.6°E); geomagnetic coordinates (14.29°N, 151.12°E)) for the low solar activity year 2005, are presented in this paper. Calculation of the diurnal maximum of the strength of the equatorial electrojet, which can serve as precursor to ionospheric scintillations in the anomaly region is also done. The Langmuir Probe experiment and Plasma Analyzer onboard DEMETER measure the electron and ion densities respectively. Irregularities in electron density distribution cause scintillations on transionospheric links and there exists a close relationship between an irregularity and scintillation. In 40% of the cases, DEMETER detects the irregularity structures (dNe/Ne ? 5% and dNi/Ni (O+) ? 5%) and GPS L band scintillations (S4 ? 0.2) are also observed around the same time, for the low solar activity period. It is found that maximum irregularity intensity is obtained in the geomagnetic latitude range of 10–20° for both electron density and ion density. As the GPS signals pass through this irregularity structure, scintillations are recorded by the GPS receiver installed at the equatorial anomaly station, Bhopal it is interesting to note that in situ density fluctuations observed on magnetic flux tubes that pass over Bhopal can be used as indicator of ionospheric scintillations at that site. Many cases of density fluctuations and associated scintillations have been observed during the descending low solar activity period. The percentage occurrence of density irregularities and scintillations shows good correspondence with diurnal maximum of the strength of electrojet, however this varies with different seasons with maximum correspondence in summer (up to 66%) followed by equinox (up to 50%) and winter (up to 46%). Also, there is a threshold value of EEJ strength to produce density irregularities ((dNe/Ne)max ? 5%) and for moderate to strong scintillations (S4 ? 0.3) to occur. For winter this value is found to be ∼40 nT whereas for equinox and summer it is around 50 nT.  相似文献   

15.
We report work utilizing 15-min resolution ionospheric data obtained with DPS-4 digisonde in 2003–2011 to study the seasonal variations in amplitudes and phases of the most powerful spectral components of the F2 layer critical frequency (foF2) and peak height (hmF2) fluctuations over Irkutsk (52.5°N, 104.0°E). We show that fluctuations of both parameters contain quasi-harmonic components with periods of Tn = 24/n h (n = 1–7). The number of distinct spectral peaks varies from 3 in summer to 7 in winter. Amplitude and phase characteristics of the diurnal (n = 1) and semidiurnal (n = 2) components is studied using the data sets extracted from the original data sets with band-pass filter. It has been found that the amplitudes of diurnal/semidiurnal foF2 and diurnal hmF2 components are maximum in winter and minimum in summer. Amplitudes of the diurnal components vary gradually; those of the foF2 semidiurnal one, abruptly, thus forming a narrow winter maximum in November–January. The phase (local time of maximum) of the diurnal foF2 component increases gradually by 4–6 h from winter to summer. The phase of the semidiurnal foF2 component is nearly stable in winter/summer and sharply decreases (increases) by 2–3 h near the spring (autumn) equinox. The phase of the diurnal component of hmF2 (local time of minimum) varies slightly between 1130 and 1300 LT; that of the semidiurnal one decreases (increases) by 4–6 h from January to March (from September to November). The results obtained show that the main features of seasonal variations in the diurnal and semidiurnal components of the mid-latitude F2 layer parameters recur consistently during the solar activity growth and decline phases.  相似文献   

16.
Recently a new model for explaining the observations of preferential heating of heavy ions in the polar solar corona was proposed ( and ). In that model the ion energization mechanism is the ion reflection off supercritical quasi-perpendicular collisionless shocks in the corona and the subsequent acceleration by the motional electric field E = −V × B/c. The mechanism of heavy ion reflection is based on ion gyration in the magnetic overshoot of the shock. The acceleration due to the motional electric field is perpendicular to the magnetic field, giving rise to large temperature anisotropy with T ? T, in agreement with SoHO observations. Such a model is tested here by means of a one dimensional test particle simulation where ions are launched toward electric and magnetic profiles representing the shock transition. We study the dynamics of O5+, as representative of coronal heavy ions for Alfvénic Mach numbers of 2–4, as appropriate to solar corona. It is found that O5+ ions are easily reflected and gain more than mass proportional energy with respect to protons.  相似文献   

17.
A laboratory experiment helps to understand the light scattering property of regolith like samples with known compositions and other physical parameters. The laboratory data so obtained can be compared with the existing in situ data on celestial objects like asteroids. Further, it may be analyzed with the help of various theoretical models to understand the light scattering processes from regolith more clearly. In this work we have performed laboratory based photometry of the light scattered from the surfaces of powdered alumina (Al2O3) at various tilt angles of the sample and at large phase angles, with the particles having diameter 0.3 μm. The wavelength of observation was 632.8 nm. These data have been fitted by a surface scattering model originally suggested by Hapke. Instead of using empirical Henyey–Greenstein phase function to fix the values of albedo and phase function to be used within Hapke formula, we have used Mie theory for the same. This approach helped us to determine the single particle properties such as particle diameter and complex refractive index from surface scattering phase curve alone. Mie theory depends only on the size parameter X(=2π(radius/wavelength)) and complex refractive index (nk) of the material. Since the absorption coefficient (k) for alumina is known to be very low but not exactly zero, the best fit to the experimental data was obtained by least square technique with k as a free parameter, as the other parameters are known. Finally, we compare our results with other published results and discuss the scope of application of the method we adopted.  相似文献   

18.
We present and discuss here the first version of a data base of extreme solar and heliospheric events. The data base contains now 87 extreme events mostly since 1940. An event is classified as extreme if one of the three critical parameters passed a lower limit. The critical parameters were the X-ray flux (parameter R), solar proton flux (parameter S) and geomagnetic disturbance level (parameter G). We find that the five strongest extreme events based on four variables (X-rays SEP, Dst, Ap) are completely separate except for the October 2003 event which is one the five most extreme events according to SEP, Dst and Ap. This underlines the special character of the October 2003 event, making it unique within 35 years. We also find that the events based on R and G are rather separate, indicating that the location of even extreme flares on the solar disk is important for geomagnetic effects. We also find that S = 3 events are not extreme in the same sense as R > 3 and G > 3 events, while S = 5 events are missing so far. This suggests that it might be useful to rescale the classification of SEP fluxes.  相似文献   

19.
The geometries, electron affinities and/or electron detachment energies for the CnS and CnO (n = 2–8) molecules and their anions were calculated by using the RCCSD (T) method. The CnS (even n = 4, 6, and 8) and CnO (even n = 6 and 8) anions are found to be substantially more stable than their corresponding neutral species. Several anions are potentially detectable as interstellar molecules.  相似文献   

20.
To investigate the effect of nitrogen sources on Closed Ecological Systems (CESs), three nitrogen sources (NaNO3, sodium nitrate; NH4Cl, ammonium chloride; and NH4NO3, ammonium nitrate) were each tested in freshwater CESs consisting of a chemically defined medium, three species of green algae (Ankistrodesmus, Scenedesmus, and Selenastrum), the grazer Daphnia magna, and associated microbes, under 12 h light/12 h dark cycles. It had been hypothesized that the development of high pH in earlier CESs was the result of nitrate utilization, and that ammonium might result in acid conditions, while ammonium nitrate might result in more moderate pH. The three nitrogen sources supported similar densities of algae (estimated by in vivo fluorescence) and similar Daphnia populations. The experiments showed that pH levels rapidly increased when grazers were absent or at low abundances irrespective of the nitrogen source. Consequently, it is hypothesized that carbon cycles, rather than nitrogen sources, are responsible for the pH dynamics. Oxygen diurnal (light:dark) cycles tended to come into balance more quickly than pH. It may be more feasible to convert O2 data to energy units (using “oxycalorific” values) than CO2 data since CO2 dynamics may include other chemical reactions than just photosynthesis and respiration. The feasibility of sustaining grazer populations for at least several weeks in small, simple CESs was demonstrated, along with the ability to monitor algae-grazer dynamics, and the recording of O2 and pH measurements.  相似文献   

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

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