首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 359 毫秒
1.
The plasmaspheric electron content (PEC) was estimated by comparison of GPS TEC observations and FORMOSAT-3/COSMIC radio occultation measurements at the extended solar minimum of cycle 23/24. Results are retrieved for different seasons (equinoxes and solstices) of the year 2009. COSMIC-derived electron density profiles were integrated up to the height of 700 km in order to retrieve estimates of ionospheric electron content (IEC). Global maps of monthly median values of COSMIC IEC were constructed by use of spherical harmonics expansion. The comparison between two independent measurements was performed by analysis of the global distribution of electron content estimates, as well as by selection specific points corresponded to mid-latitudes of Northern America, Europe, Asia and the Southern Hemisphere. The analysis found that both kinds of observations show rather similar diurnal behavior during all seasons, certainly with GPS TEC estimates larger than corresponded COSMIC IEC values. It was shown that during daytime both GPS TEC and COSMIC IEC values were generally lower at winter than in summer solstice practically over all specific points. The estimates of PEC (h > 700 km) were obtained as a difference between GPS TEC and COSMIC IEC values. Results of comparative study revealed that for mid-latitudinal points PEC estimates varied weakly with the time of a day and reached the value of several TECU for the condition of solar minimum. Percentage contribution of PEC to GPS TEC indicated the clear dependence from the time with maximal values (more than 50–60%) during night-time and lesser values (25–45%) during day-time.  相似文献   

2.
We examined performance of two empirical profile-based ionospheric models, namely IRI-2016 and NeQuick-2, in electron content (EC) and total electron content (TEC) representation for different seasons and levels of solar activity. We derived and analyzed EC estimates in several representative altitudinal intervals for the ionosphere and the plasmasphere from the COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations. It allows us to estimate a quantitative impact of the ionospheric electron density profiles formulation in several altitudinal intervals and to examine the source of the model-data discrepancies of the EC specification from the bottom-side ionosphere towards the GPS orbit altitudes. The most pronounced model-data differences were found at the low latitude region as related to the equatorial ionization anomaly appearance. Both the IRI-2016 and NeQuick-2 models tend to overestimate the daytime ionospheric EC and TEC at low latitudes during all seasons of low solar activity. On the contrary, during high solar activity the model results underestimated the EC/TEC observations at low latitudes. We found that both models underestimated the EC for the topside ionosphere and plasmasphere regions for all levels of solar activity. For low solar activity, the underestimated EC from the topside ionosphere and plasmasphere can compensate the overestimation of the ionospheric EC and, consequently, can slightly decrease the resulted model overestimation of the ground-based TEC. For high solar activity, the underestimated EC from the topside ionosphere and plasmasphere leads to a strengthening of the model underestimation of the ground-based TEC values. We demonstrated that the major source of the model-data discrepancies in the EC/TEC domain comes from the topside ionosphere/plasmasphere system.  相似文献   

3.
The ionosphere/plasmasphere electron content (PEC) variations during strong geomagnetic storms in November 2004 were estimated by combining of mid-latitude Kharkov incoherent scatter radar observations and GPS TEC data derived from global TEC maps. The comparison between two independent measurements was performed by analysis of the height-temporal distribution for specific location corresponding to the mid-latitudes of Europe. The percentage contribution of PEC to GPS TEC indicated the clear dependence from the time with maximal values (more than 70%) during night-time. During day-time the lesser values (30–45%) were observed for quiet geomagnetic conditions and rather high values of the PEC contribution to GPS TEC (up to 90%) were observed during strong negative storm. These changes can be explained by the competing effects of electric fields and winds, which tend to raise the layer to the region with lower loss rate and movement of the ionospheric plasma to the plasmasphere.  相似文献   

4.
This paper presents small scale (duration ?1 h, ΔTEC ? 1TECU) night-time total electron content (TEC) enhancements observed at the equatorial anomaly region in China, for the first time. The data is from a GPS receiver chain established in 2005 by Institute of Center for Space Science and Applied Research, Chinese Academy of Sciences and a GPS receiver of International GPS Service (IGS), located between Fuzhou (26.1°N, 119.3°E) and Nanning (22.8°N, 108.3°E). Two other GPS observations of IGS taken at higher latitude were also used to investigate the localization of such phenomenon. The characteristics of the night-time TEC enhancement are examined with two case studies. The TEC increases about 1–3TECU, intermittently. While the night-time TEC enhancement mainly occurs at the equatorial anomaly region, it can be observed at middle latitude as well. The spatial size of the enhancement region is less than 5° in longitude. The primary statistical study shows that the TEC enhancement is more often observed in spring and autumn, but rarely in summer. It has no dependence on geomagnetic activity. The enhancement can occur both before and after midnight.  相似文献   

5.
By using data from GPS receivers we detected huge-amplitude solitary large-scale traveling acoustic-gravity waves (LS AGW) which manifested themselves as perturbations of total electron content (TEC) of duration of about 40 min. Originated in the auroral area after significant alterations of geomagnetic field intensity during geomagnetic storms on 29–30 October 2003, LS disturbances propagated with a velocity about 1000–1200 m/s and caused generation of secondary small-scale (SS) waves with time period of 2–10 min. Such SS structure followed the solitary intensive AGW at a distance more than 4000 km. However, we observed such phenomenon only within the territory with high values of “vertical” TEC and steep gradients of TEC. Apparently, these conditions are necessary for generation of SS waves due to propagation of LS AGW.  相似文献   

6.
The Constellation Observing System for Meteorology Ionosphere and Climate (COSMIC) is a six satellite radio occultation mission that was launched in April 2006. The close proximity of these satellites during some months after launch provides a unique opportunity to evaluate the precision of Global Positioning System (GPS) radio occultation (RO) retrievals of ionospheric electron density from nearly collocated and simultaneous observations. RO data from 30 consecutive days during July and August 2006 are divided into ten groups in terms of daytime or nighttime and latitude. In all cases, the best precision values (about 1%) are found at the F peak height and they slightly degrade upwards. For all daytime groups, it is seen that electron density profiles above about 120 km height exhibit a substantial improvement in precision. Nighttime groups are rather diverse: in particular, the precision becomes better than 10% above different levels between 120 and 200 km height. Our overall results show that up to 100–200 km (depending on each group), the uncertainty associated with the precision is in the order of the measured electron density values. Even worse, the retrieved values tend sometimes to be negative. Although we cannot rely directly on electron density values at these altitudes, the shape of the profiles could be indicative of some ionospheric features (e.g. waves and sporadic E layers). Above 200 km, the profiles of precision are qualitatively quite independent from daytime or latitude. From all the nearly collocated pairs studied, only 49 exhibited a difference between line of sight angles of both RO at the F peak height larger than 10°. After analyzing them we find no clear indications of a significant representativeness error in electron density profiles due to the spherical assumption above 120 km height. Differences in precision between setting and rising GPS RO may be attributed to the modification of the processing algorithms applied to rising cases during the initial period of the COSMIC mission.  相似文献   

7.
Analysis of a long-time series of hourly median characteristics of the ionospheric plasma at two mid-latitude locations in the Northern and Southern hemisphere, Juliusruh (54.6N; 13.4E) and Hobart (42.9S; 147.3E), reveals patterns of their synchronous and independent variability. We studied timelines of GPS vTEC, ionogram-derived F2-layer peak electron density NmF2, ionospheric equivalent slab thickness τ, and their ratios at two locations during the complete 23rd solar cycle and its following period of the extremely low solar activity in 2008–2009. This study has also involved the comparative analysis of the observed data versus the model predictions by IRI-2012. During the high solar activity in 2000–2002, seasonal variations show a complicated cross-hemisphere behavior influenced by the winter and semi-annual anomalies, with the largest noon-time values of TEC and NmF2 observed around equinoxes. Strength of the winter anomaly in NmF2 was significantly greater at Juliusruh in comparison with Hobart. The winter anomaly in GPS vTEC values was much weaker than in NmF2 for the Northern hemisphere mid-latitudes and was entirely absent at the Southern hemisphere. Cross-hemisphere analysis of the equivalent slab thickness shows its clear seasonal dependence for all levels of solar activity: the day-time maximum τmax is observed during local summer, whereas the day-time minimum τmin is observed during local winter. The night-time values of τ were higher compared to the day-time values during the winter and equinox seasons. Comparative model-data study shows rather good IRI performance of the day-time NmF2 for mid-latitudes of both hemispheres and rather noticeable overestimations for the mid-night NmF2 values during high solar activity. Analysis of IRI vTEC demonstrates the model limitations, related with the absence of the plasmaspheric part, and actual demand in a reliable and standard ionosphere–plasmasphere model for analysis of GPS vTEC.  相似文献   

8.
The GRAS radio occultation instrument is flying on Metop-A and belongs to the EPS (EUMETSAT Polar System). GRAS observes GPS satellites in occultation. Within this work, validation of GRAS closed-loop bending angle data against co-located ECMWF profiles extracted from model fields and occultations from the COSMIC constellation of radio occultation instruments is shown. Results confirm the high data quality and robustness, where GRAS shows lower bending angle noise against ECMWF than COSMIC and in terms of occultations per day, one GRAS ≈ two COSMIC satellites. This is partly due to the operational setup of EPS. For the investigation we focus on two observation periods where updates in the ECMWF (March 2009) and COSMIC processing (October 2009) have improved the statistics further. Bending angles biases agree to within 0.5% against ECMWF and to within 0.1% against COSMIC after the updates for altitudes between 8 and 40 km. In addition, we also analyze the impact of the Metop orbit processing on the derived GRAS bending angle data, where different GPS and Metop orbit solutions are analyzed. Results show that a batch based orbit processing would improve in particular the bending angle bias behavior at higher altitudes. Requirements for the operational processing of GRAS data are briefly outlined, options to ease the use of other positioning system satellites in the near future are discussed. A simplified analysis on the observation of several of these systems, e.g. GPS and Galileo, from one platform shows that about 16% of occultations are found within 300 km, ±3 h, thus providing similar information. A constellation of 2 GRAS like instruments would have only about 10% close-by.  相似文献   

9.
The present study reports the analysis of GPS based TEC for our station Surat (21.16°N, 72.78°E) located at the northern crest of equatorial anomaly region in India at times close to some earthquake events (M ? 5) during the year 2009 in India and its neighbouring regions. The TEC data used in the study are obtained from GPS Ionospheric Scintillation and TEC Monitoring (GISTM) system. The TEC data has been analysed corresponding to 11 earthquakes in low solar activity period and quiet geomagnetic condition. We found that, out of 11 cases of earthquakes (M > 5) there were seven cases in which enhancement in TEC occurred on earthquake day and in other four cases there was depletion in TEC on earthquake day. The variation in refractivity prior to earthquake was significant for the cases in which the epicentre lied within a distance of 600 km from the receiving station. By looking into the features on temporal enhancement and depletion of TEC a prediction was made 3–2 days prior to an earthquake (on 28 October 2009 in Bhuj – India). The paper includes a brief discussion on the method of potentially identifying an impending earthquake from ionospheric data.  相似文献   

10.
We present the spatial maps of the ionosphere–plasmasphere slab thickness τ (ratio of the vertical total electron content, TEC, to the F-region peak electron density, NmF2) during the intense ionospheric storms of October–November 2003. The model-assisted technology for estimate of the upper boundary of the ionosphere, hup, from the slab thickness components in the bottomside and topside ionosphere – eliminating the plasmasphere contribution of τ – is applied at latitudes 35° to 70°N and longitudes −10° to 40°E, from the data of 20 observatories of GPS-TEC and ionosonde networks, for selected days and hours of October and November 2003. The daily–hourly values of NmF2, hmF2 and TECgps are used as the constrained parameters for the International Reference Ionosphere extended to the plasmasphere, IRI-Plas, during the ionospheric quiet days, positive and negative storm phases for estimate of τ and hup. Good correlation has been found between the slab thickness and the upper boundary of the ionosphere for the intense ionospheric storms at October–November 2003. During the negative phase of the ionospheric storm, when the ionospheric plasma density is exhausted, the nighttime upper boundary of the ionosphere is greatly uplifted towards the magnetosphere tail, while the daytime upper boundary of the ionosphere is reduced below 500 km over the Earth.  相似文献   

11.
The temporal and seasonal variations of Total Electron Content (TEC) are studied at Agra (Geographic Lat. 27.17°N, Long. 78.89°E, Dip: 41.4°), India, which is in the equatorial anomaly region, for a period of 12 months from 01 January to 31 December, 2007 using a Global Positioning System (GPS) receiver. The mean TEC values show a minimum at 0500 h LT (LT = UT + 5.5 h) and a peak value at about 1400 h LT. The lowest TEC values are observed in winter whereas largest values are observed in equinox and summer. Anomalous variations are found during the period of magnetic disturbances. These results are compared with the TEC derived from IRI-2007 using three different options of topside electron density, NeQuick, IRI01-corr, and IRI-2001. A good agreement is found between the TEC obtained at Agra and those derived from IRI models.  相似文献   

12.
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.  相似文献   

13.
Ionosphere response to severe geomagnetic storms that occurred in 2001–2003 was analyzed using data of global ionosphere maps (GIM), altimeter data from the Jason-1 and TOPEX satellites, and data of GPS receivers on-board CHAMP and SAC-C satellites. This allowed us to study in detail ionosphere redistribution due to geomagnetic storms, dayside ionospheric uplift and overall dayside TEC increase. It is shown that after the interplanetary magnetic field turns southward and intensifies, the crests of the equatorial ionization anomaly (EIA) travel poleward and the TEC value within the EIA area increases significantly (up to ∼50%). GPS data from the SAC-C satellite show that during the main phase of geomagnetic storms TEC values above the altitude of 715 km are 2–3 times higher than during undisturbed conditions. These effects of dayside ionospheric uplift occur owing to the “super-fountain effect” and last few hours while the enhanced interplanetary electric field impinged on the magnetopause.  相似文献   

14.
The International Reference Ionosphere model extended to the plasmasphere, IRI-Plas, presents global electron density profiles and total electron content, TECiri, up to the altitude of the GPS satellites (20,000 km). The model code is modified by input of GPS-derived total electron content, TECgps, so that the topside scale height, Hsc, is obtained minimizing in one step the difference between TECiri and TECgps observation. The topside basis scale height, Hsc, presents the distance in km above the peak height at which the peak plasma density, NmF2, decays by a factor of e (∼2.718). The ionosonde derived F2 layer peak density and height and GPS-derived TECgps data are used with IRI-Plas code during the main phase of more than 100 space weather storms for a period of 1999–2006. Data of seven stations are used for the analysis, and data from five other stations served as testing database. It is found that the topside basis scale height is growing (depressing) when the peak electron density (critical frequency foF2) and electron content are decreased (increased) compared to the median value, and vice versa. Relative variability of the scale height, rHsc, and the instantaneous Hsc are inferred analytically in a function of the instantaneous foF2, median fmF2 and median Hmsc avoiding a reference to geomagnetic indices. Results of validation suggest reliability of proposed algorithm for implementation in an operational mode.  相似文献   

15.
The response of the ionospheric F-region in the equatorial and low latitude regions in the Brazilian sector during the super geomagnetic storm on 06–07 April 2000 has been studied in the present investigation. The geomagnetic storm reached a minimum Dst of −288 nT at 0100 UT on 07 April. In this paper, we present vertical total electron content (VTEC) and phase fluctuations (in TECU/min) from GPS observations obtained at Imperatriz (5.5°S, 47.5°W; IMPZ), Brasília (15.9°S, 47.9°W; BRAZ), Presidente Prudente (22.12°S, 51.4°W; UEPP), and Porto Alegre (30.1°S, 51.1°W; POAL) during the period 05–08 April. Also, several GPS-based TEC maps are presented from the global GPS network, showing widespread and drastic TEC changes during the different phases of the geomagnetic storm. In addition, ion density measurements on-board the satellite Defense Meteorological Satellite Program (DMSP) F15 orbiting at an altitude of 840 km and the first Republic of China satellite (ROCSAT-1) orbiting at an altitude of 600 km are presented. The observations indicate that one of the orbits of the DMSP satellite is fairly close to the 4 GPS stations and both the DMSP F15 ion-density plots and the phase fluctuations from GPS observations show no ionospheric irregularities in the Brazilian sector before 2358 UT on the night of 06–07 April 2000. During the fast decrease of Dst on 06 April, there is a prompt penetration of electric field of magnetospheric origin resulting in decrease of VTEC at IMPZ, an equatorial station and large increase in VTEC at POAL, a low latitude station. This resulted in strong phase fluctuations on the night of 06–07 April, up to POAL. During the daytime on 07 April during the recovery phase, the VTEC observations show positive ionospheric storm at all the GPS stations, from IMPZ to POAL, and the effect increasing from IMPZ to POAL. This is possibly linked to the equatorward directed meridional wind. During the daytime on 08 April (the recovery phase continues), the VTEC observations show very small negative ionospheric storm at IMPZ but the positive ionospheric storm effect is observed from BRAZ to POAL possibly linked to enhancement of the equatorial ionospheric anomaly.  相似文献   

16.
On 21 June 2010 the TerraSAR-X satellite was joined by the TanDEM-X satellite. A Global Positioning System (GPS) radio occultation (RO) experiment using the twin satellites has been carried out to estimate the precision of GPS atmospheric soundings. For the Day Of Year (DOY) 330–336, 2011, we analyze phase and amplitude data recorded by GPS receivers separated by a few hundred meters in a low earth orbit and derive collocated atmospheric refractivity profiles. In the altitude range 10–20 km the standard deviation between TerraSAR-X and TanDEM-X refractivity does not exceed 0.15%. The standard deviation is rapidly increasing for lower and higher altitudes; close to the surface and at an altitude of 30 km the standard deviation reaches 0.8% and 0.5%, respectively. Systematic deviations between TerraSAR-X and TanDEM-X refractivity in the considered altitude range (0–30 km) are negligible. The results confirm the anticipated high precision of the GPS RO technique. However, the difference in the retrieved refractivity in the lower troposphere for different Open Loop (OL) signal tracking parameters, altered onboard TanDEM-X for DOY 49–55, 2012, calls for an in depth analysis. At the moment we can not exclude that a potential bias in the OL Doppler model introduces a bias in our retrieved refractivity at altitudes <8<8 km.  相似文献   

17.
We present an analysis of the ionosphere and thermosphere response to Solar Proton Events (SPE) and magnetospheric proton precipitation in January 2005, which was carried out using the model of the entire atmosphere EAGLE. The ionization rates for the considered period were acquired from the AIMOS (Atmospheric Ionization Module Osnabrück) dataset. For numerical experiments, we applied only the proton-induced ionization rates of that period, while all the other model input parameters, including the electron precipitations, corresponded to the quiet conditions. In January 2005, two major solar proton events with different energy spectra and proton fluxes occurred on January 17 and January 20. Since two geomagnetic storms and several sub-storms took place during the considered period, not only solar protons but also less energetic magnetospheric protons contributed to the calculated ionization rates. Despite the relative transparency of the thermosphere for high-energy protons, an ionospheric response to the SPE and proton precipitation from the magnetotail was obtained in numerical experiments. In the ionospheric E layer, the maximum increase in the electron concentration is localized at high latitudes, and at heights of the ionospheric F2 layer, the positive perturbations were formed in the near-equatorial region. An analysis of the model-derived results showed that changes in the ionospheric F2 layer were caused by a change in the neutral composition of the thermosphere. We found that in the recovery phase after both solar proton events and the enhancement of magnetospheric proton precipitations associated with geomagnetic disturbances, the TEC and electron density in the F region and in topside ionosphere/plasmasphere increase at low- and mid-latitudes due to an enhancement of atomic oxygen concentration. Our results demonstrate an important role of magnetospheric protons in the formation of negative F-region ionospheric storms. According to our results, the topside ionosphere/plasmasphere and bottom-side ionosphere can react to solar and magnetospheric protons both with the same sign of disturbances or in different way. The same statement is true for TEC and foF2 disturbances. Different disturbances of foF2 and TEC at high and low latitudes can be explained by topside electron temperature disturbances.  相似文献   

18.
The total electron content (TEC) estimation by the Global Positioning System (GPS) can be seriously affected by the differential code biases (DCB), referred to as inter-frequency biases (IFB), of the satellite and receiver so that an accuracy of GPS–TEC value is dependent on the error of DCBs estimation. In this paper, we proposed the singular value decomposition (SVD) method to estimate the DCB of GPS satellites and receivers using the Korean GPS network (KGN) in South Korea. The receiver DCBs of about 49 GPS reference stations in KGN were determined for the accurate estimation of the regional ionospheric TEC. They obtained from the daily solution have large biases ranging from +5 to +27 ns for geomagnetic quiet days. The receiver DCB of SUWN reference station was compared with the estimates of IGS and JPL global ionosphere map (GIM). The results have shown comparatively good agreement at the level within 0.2 ns. After correction of receiver DCBs and knowing the satellite DCBs, the comparison between the behavior of the estimated TEC and that of GIMs was performed for consecutive three days. We showed that there is a good agreement between KASI model and GIMs.  相似文献   

19.
GPS observations from EUREF permanent GPS network were used to observe the response of TEC (Total Electron Content) to the total solar eclipse on October 3, 2005, under quiet geomagnetic conditions of the daytime ionosphere. The effect of the eclipse was detected in diurnal variations and more distinctly in the variations of TEC along individual satellite passes. The trough-like variations with a gradual decrease and followed by an increase of TEC at the time of the eclipse were observed over a large region. The depression of TEC amounted to 3–4 TECU. The maximum depression was observed over all stations located at the maximum path of the solar eclipse. The delay of a minimum level of TEC with respect to the maximum phase of the eclipse was about 20–30 min.  相似文献   

20.
The ionosphere induces a time delay in transionospheric radio signals such as the Global Positioning System (GPS) signal. The Total Electron Content (TEC) is a key parameter in the mitigation of ionospheric effects on transionospheric signals. The delay in GPS signal induced by the ionosphere is proportional to TEC along the path from the GPS satellite to a receiver. The diurnal monthly and seasonal variations of ionospheric electron content were studied during the year 2010, a year of extreme solar minimum (F10.7 = 81 solar flux unit), with data from the GPS receiver and the Digisonde Portable Sounder (DPS) collocated at Ilorin (Geog. Lat. 8.50°N, Long. 4.50°E, dip −7.9°). The diurnal monthly variation shows steady increases in TEC and F2-layer critical frequency (foF2) from pre-dawn minimum to afternoon maximum and then decreases after sunset. TEC show significant seasonal variation during the daytime between 0900 and 1900 UT (LT = UT + 1 h) with a maximum during the March equinox (about 35 TECU) and minimum during the June solstice (about 24 TECU). The GPS-TEC and foF2 values reveal a weak seasonal anomaly and equinoctial asymmetry during the daytime. The variations observed find their explanations in the amount of solar radiation and neutral gas composition. The measured TEC and foF2 values were compared with last two versions of the International Reference Ionosphere (IRI-2007 and IRI-2012) model predictions using the NeQuick and CCIR (International Radio Consultative Committee) options respectively in the model. In general, the two models give foF2 close to the experimental values, whereas significant discrepancies are found in the predictions of TEC from the models especially during the daytime. The error in height dependent thickness parameter, daytime underestimation of equatorial drift and contributions of electrons from altitudes above 2000 km have been suggested as the possible causes.  相似文献   

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

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