首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 659 毫秒
1.
The Earth's ionosphere and especially its equatorial part is a highly dynamical medium. Geostationary satellites are known to be a powerful tool for ionospheric studies. Recent developments in BDS-GEO satellites allow such studies on the new level due to the best noise pattern in TEC estimations, which corresponds to those of GPS/GLONASS systems. Here we used BDS-GEO satellites to demonstrate their capability for studying equatorial ionosphere variability on different time scales. Analyzing data from the equatorial SIN1 IGS station we present seasonal variations in geostationary slant TEC for the periods of high (October 2013 - October 2014) and low (January 2017 - January 2018) solar activity, which show semi-annual periodicity with amplitudes about 10 TECU during solar maximum and about 5 TECU during the solar minimum. The 27-day variations are also prominent in geostationary slant TEC variations, which correlates quite well with the variations in solar extreme UV radiation. We found semi-annual pattern in small scale ionospheric disturbances evaluated based on geostationary ROTI index: maximal values correspond to spring and fall equinoxes and minimum values correspond to summer and winter solstices. The seasonal asymmetry in ROTI values was observed: spring equinox values were almost twice as higher than fall equinox ones. We also present results on the 2017 May 28–29 G3 geomagnetic storm, when ~30 TECU positive anomaly was recorded, minor and final major sudden stratospheric warmings in February and March 2016, with positive daytime TEC anomalies up to 15–20 TECU, as well as the 2017 September 6 X9.3 solar flare with 2 TECU/min TEC rate. Our results show the large potential of geostationary TEC estimations with BDS-GEO signals for continuous monitoring of space weather effects in low-latitude and equatorial ionosphere.  相似文献   

2.
The responses of the ionospheric F region using GPS–TEC measurements during two moderate geomagnetic storms at equatorial, low-, and mid-latitude regions over the South American and African sectors in May 2010, during the ascending phase of solar cycle 24, are investigated. The first moderate geomagnetic storm studied reached a minimum Dst value of −64 nT at 1500 UT on 02 May 2010 and the second moderate geomagnetic storm reached a minimum Dst value of −85 nT at 1400 UT on 29 May 2010. In this paper, we present vertical total electron content (VTEC) and phase fluctuations (in TECU/min) from Global Positioning System (GPS) observations from the equatorial to mid-latitude regions in the South American and African sectors. Our results obtained during these two moderate geomagnetic storms from both sectors show significant positive ionospheric storms during daytime hours at the equatorial, low-, and mid-latitude regions during the main and recovery phases of the storms. The thermospheric wind circulation change towards the equator is a strong indicator that suggests an important mechanism is responsible for these positive phases at these regions. A pre-storm event that was observed in the African sector from low- to the mid-latitude regions on 01 May 2010 was absent in the South American sector. This study also showed that there was no generation or suppression of ionospheric irregularities by storm events. Therefore, knowledge about the suppression and generation of ionospheric irregularities during moderate geomagnetic storms is still unclear.  相似文献   

3.
This paper presents traveling ionospheric disturbances (TIDs) observations from GPS measurements over the South African region during the geomagnetically disturbed period of 29–31 October 2003. Two receiver arrays, which were along two distinct longitudinal sectors of about 18°-20° and 27°-28° were used in order to investigate the amplitude, periods and virtual propagation characteristics of the storm induced ionospheric disturbances. The study revealed a large sudden TEC increase on 28 October 2003, the day before the first of the two major storms studied here, that was recorded simultaneously by all the receivers used. This pre-storm enhancement was linked to an X-class solar flare, auroral/magnetospheric activities and vertical plasma drift, based on the behaviour of the geomagnetic storm and auroral indices as well as strong equatorial electrojet. Diurnal trends of the TEC and foF2 measurements revealed that the geomagnetic storm caused a negative ionospheric storm; these parameters were depleted between 29 and 31 October 2003. Large scale traveling ionospheric disturbances were observed on the days of the geomagnetic storms (29 and 31 October 2003), using line-of-sight vertical TEC (vTEC) measurements from individual satellites. Amplitude and dominant periods of these structures varied between 0.08–2.16 TECU, and 1.07–2.13 h respectively. The wave structures were observed to propagate towards the equator with velocities between 587.04 and 1635.09 m/s.  相似文献   

4.
Moderate geomagnetic storms occurred during January 22–25, 2012 period. The geomagnetic storms are characterized by different indices and parameters. The SYM-H value on January 22 increased abruptly to 67 nT at sudden storm commencement (SSC), followed by a sharp decrease to −87 nT. A second SSC on January 24 followed by a shock on January 25 was also observed. These SSCs before the main storms and the short recovery periods imply the geomagnetic storms are CME  -driven. The sudden jump of solar wind dynamic pressure and IMF BzBz are also consistent with occurrence of CMEs. This is also reflected in the change in total electron content (TEC) during the storm relative to quiet days globally. The response of the ionospheric to geomagnetic storms can also be detected from wave components that account for the majority of TEC variance during the period. The dominant coherent modes of TEC variability are diurnal and semidiurnal signals which account upto 83% and 30% of the total TEC variance over fairly exclusive ionospheric regions respectively. Comparison of TEC anomalies attributed to diurnal (DW1) and semidiurnal (SW2) tides, as well as stationary planetary waves (SPW1) at 12 UTC shows enhancement in the positive anomalies following the storm. Moreover, the impact of the geomagnetic storms are distinctly marked in the daily time series of amplitudes of DW1, SW2 and SPW1. The abrupt changes in amplitudes of DW1 (5 TECU) and SW2 (2 TECU) are observed within 20°S–20°N latitude band and along 20°N respectively while that of SPW1 is about 3 TECU. Coherent oscillation with a period of 2.4 days between interplanetary magnetic field and TEC was detected during the storm. This oscillation is also detected in the amplitudes of DW1 over EIA regions in both hemispheres. Eventhough upward coupling of quasi two day wave (QTDWs) of the same periodicity, known to have caused such oscillation, are detected in both ionosphere and upper stratosphere, this one can likely be attributed to the geomagnetic storm as it happens after the storm commencement. Moreover, further analysis has indicated that QTDWs in the ionosphere are strengthened as a result of coherent oscillation of interplanetary magnetic field with the same frequency as QTDWs. On the otherhand, occurrences of minor SSW and geomagnetic storms in quick succession complicated clear demarcation of attribution of the respective events to variability of QTDWs amplitudes over upper stratosphere.  相似文献   

5.
From September 7 to 8, 2017, a G4-level strong geomagnetic storm occurred, which seriously impacted on the Earth’s ionosphere. In this work, the global ionospheric maps released by Chinese Academy of Sciences are used to investigate the ionospheric responses over China and its adjacent regions during the strong storm. The prominent TEC enhancements, which mainly associated with the neutral wind and eastward prompt penetration electric field, are observed at equatorial ionization anomaly crests during the main phase of the storm on 8 September 2017. Compared with those on 8 September, the TEC enhancements move to lower-latitude regions during the recovery phase on 9 September. A moderate storm occurred well before the start of the strong storm causes similar middle-latitude TEC enhancements on 7 September. However, the weak TEC depletion is observed at middle and low latitude on 9–10 September, which could be associated with the prevailing westward disturbance electric field or storm-time neural composition changes. In addition, the storm-time RMS and STD values of the ionospheric TEC grids over China increase significantly due to the major geomagnetic storm. The maximum of the RMS reaches 12.0 TECU, while the maximum of the STD reaches 8.3 TECU at ~04UT on 8 September.  相似文献   

6.
利用武汉电离层观象台研制的GPS TEC的现报方法及现报系统,对东亚地区GPS台网的观测数据进行处理分析,特别对2000年7月14-18日和2003年10月28日至11月1日两次特大磁暴期间的数据进行了对比考察,文中分析了两次磁暴间的电离层响应,得到对应不同磁暴时段电离层TEC的不同变化情况,着重揭示了TEC赤道异常峰的压缩和移动以及赤道异常随时间的压缩—反弹—恢复的过程,并结合高纬电离层的部分响应机制进行了说明,结果显示,两次磁暴期的电离层响应表现出了各自不同的特点,从而反映出因季节变化引起的高纬电离层暴时能量注入的不同而造成的全球性电离层扰动的不同形态,由此看出,磁暴期间电离层TEC的变化直接与太阳扰动发生的时间及其对高纬电离层的耦合有关,若短时期内连续发生多次磁暴,则电离层反应更加复杂,不能简单地当做单一磁暴叠加处理。  相似文献   

7.
This work presents, for the first time, the analysis of the occurrence of ionospheric irregularities during geomagnetic storms at Tucumán, Argentina, a low latitude station in the Southern American longitudinal sector (26.9°S, 294.6°E; magnetic latitude 15.5°S) near the southern crest of the equatorial ionization anomaly (EIA). Three geomagnetic storms occurred on May 27, 2017 (a month of low occurrence rates of spread-F), October 12, 2016 (a month of transition from low to high occurrence rates of spread-F) and November 7, 2017 (a month of high occurrence rates of spread-F) are analyzed using Global Positioning System (GPS) receivers and ionosondes. The rate of change of total electron content (TEC) Index (ROTI), GPS Ionospheric L-band scintillation, the virtual height of the F-layer bottom side (h'F) and the critical frequency of the F2 layer (foF2) are considered. Furthermore, each ionogram is manually examined for the presence of spread-F signatures.The results show that, for the three events studied, geomagnetic activity creates favorable conditions for the initiation of ionospheric irregularities, manifested by ionogram spread-F and TEC fluctuation. Post-midnight irregularities may have occurred due to the presence of eastward disturbance dynamo electric fields (DDEF). For the May storm, an eastward over-shielding prompt penetration electric field, (PPEF) is also acting. A possibility is that the PPEF is added to the DDEF and produces the uplifting of the F region that helps trigger the irregularities. Finally, during October and November, strong GPS L band scintillation is observed associated with strong range spread-F (SSF), that is, irregularities extending from the bottom-side to the topside of the F region.  相似文献   

8.
The study investigated the effects of intense geomagnetic storms of 2015 on the occurrences of large scale ionospheric irregularities over the African equatorial/low-latitude region. Four major/intense geomagnetic storms of 2015 were analyzed for this study. These storms occurred on 17th March 2015 (?229?nT), 22nd June 2015 (?204?nT), 7th October 2015 (?124?nT), and 20th December 2015 (?170?nT). Total Electron Content (TEC) data obtained from five African Global Navigation Satellite Systems (GNSS) stations, grouped into eastern and western sectors were used to derive the ionospheric irregularities proxy indices, e.g., rate of change of TEC (ROT), ROT index (ROTI) and ROTI daily average (ROTIAVE). These indices were characterized alongside with the disturbance storm time (Dst), the Y component of the Interplanetary Electric Field (IEFy), polar cap (PC) index and the H component of the Earth’s magnetic field from ground-based magnetometers. Irregularities manifested in the form of fluctuations in TEC. Prompt penetration of electric field (PPEF) and disturbance dynamo electric field (DDEF) modulated the behaviour of irregularities during the main and recovery phases of the geomagnetic storms. The effect of electric field over both sectors depends on the local time of southward turning of IMF Bz. Consequently, westward electric field inhibited irregularities during the main phase of March and October 2015 geomagnetic storms, while for the June 2015 storm, eastward electric field triggered weak irregularities over the eastern sector. The effect of electric field on irregularities during December 2015 storm was insignificant. During the recovery phase of the storms, westward DDEF suppressed irregularities.  相似文献   

9.
The dual-frequency satellite-based measurements from Global Positioning System (GPS) may provide feasible ways of studying and potentially detecting of earthquake (EQ) related anomalies in the ionosphere. In this paper, GPS based Total Electron Content (TEC) data are studied for three major M?>?7.0 EQs in Nepal and Iran-Iraq border during 2015–2017 by implementing statistical procedures on temporal and spatial scale. Previous studies presented different time interval of pre-seismic ionospheric anomalies, however, this study showed that EQs ionospheric precursors may occur within 10?days. Furthermore, the ionospheric anomalies on the suspected day occurred during UT?=?08:00–12:00?h before the main shock. The Global Ionospheric Map TEC (GIM-TEC) data retrieved over the epicenter of M7.8 (Nepal EQ) showed a significant increase of 6 TECU on April 24, 2015 (one day before the main shock), which is recorded by the ground GPS station data of Islamabad (station lies within the EQ preparation zone). Furthermore, the spatial GIM-TEC result imply significant anomalies over the epicenter during the time interval between UT?=?08:00–12:00?h (LT?=?13:00–17:00). For M7.3 (Nepal EQ), the TEC anomalies were detected on May 10, 2015 (2?days before the event) in the temporal data. The spatial TEC data imply the huge clouds over the epicenter at about UT?=?08:00–12:00?h on May 10, 2015, that may be associated with this EQ in the quiet geomagnetic storm conditions. Similarly, temporal and spatial TEC showed anomaly on November 3, 2017, during UT?=?08:00–12:00 (9?days before the Iran-Iraq border EQ) after implementing the statistical method on it. Conversely, there exists a short-term but low magnitude TEC anomaly synchronized with a geomagnetic storm on November 7–8, 2017 (4 to 5?days prior to M7.3 Iran-Iraq border EQ). The diurnal and hourly GIM-TEC and VTEC data also imply the execution of ionospheric anomalies within 10?days prior to all events. All these positive anomalies in TEC may be due to the existence of a huge energy from the epicenter during the EQ preparation period.  相似文献   

10.
Data from the archive of the International GNSS Services (IGS) were used to study the seasonal variations of Total Electron Content (TEC) over three stations located at different latitudes in the southern hemisphere during the geomagnetic storms of 11 January, 6 April, 8 June, and 13 October 2000, representing storms that occurred in summer, autumn equinox, winter and spring equinox, respectively. The percentage TEC deviation with respect to reference values differs substantially from season to season. A strong seasonal anomaly and clear equinoctial asymmetry in TEC response to the storms were observed. Weak and short-lived positive TEC deviations as well as strong and long-lasting negative trends were observed in summer storm during the main and recovery phases respectively over the high and low latitudes whereas in winter storm, the highest positive TEC deviations was recorded during the main phase over the entire latitudes. TEC enhancement dominated all the stations during the autumn (March) equinox storm while TEC depletion was majorly observed during the spring (September) equinox. All these variations find their explanations in the thermospheric composition change and circulation. Future work with direct or modeled measurement of atomic Oxygen to molecular Nitrogen ratio (O/N2), large number of storms and other possible factors such as variations in storm’s intensity and local time dependence of the storm onset is expected to validate the observations in this study.  相似文献   

11.
This study investigates the morphology of the GPS TEC responses in the African Equatorial Ionization Anomaly (EIA) region to intense geomagnetic storms during the ascending and maximum phases of solar cycle 24 (2012–2014). Specifically, eight intense geomagnetic storms with Dst ≤ ?100 nT were considered in this investigation using TEC data obtained from 13 GNSS receivers in the East African region within 36–42°E geographic longitude; 29°N–10°S geographic latitude; ± 20°N magnetic latitude. The storm-time behavior of TEC shows clear positive and negative phases relative to the non-storm (median) behavior, with amplitudes being dependent on the time of sudden commencement of the storm and location. When a storm starts in the morning period, total electron content increases for all stations while a decrease in total electron content is manifested for a storm that had its sudden commencement in the afternoon period. The TEC and the EIA crest during the main phase of the storm is significantly impacted by the geomagnetic storm, which experiences an increase in the intensity of TEC while the location and spread of the crest usually manifest a poleward expansion.  相似文献   

12.
利用行星际太阳风参数与太阳活动指数、地磁活动指数、电离层总电子含量格点化地图数据,首次基于一种能处理时间序列的深度学习递归神经网络(Recurrent Neural Network,RNN),建立提前24h的单站电离层TEC预报模型.对北京站(40°N,115°E)的预测结果显示,RNN对扰动电离层的预测误差低于反向传播神经网络(Back Propagation Neural Network,BPNN)0.49~1.46TECU,将太阳风参数加入预报因子模型后对电离层正暴预测准确率的提升可达16.8%.RNN对2001和2015年31个强电离层暴预报的均方根误差比BPNN低0.2TECU,将太阳风参数加入RNN模型可使31个事件的平均预报误差降低0.36~0.47TECU.研究结果表明深度递归神经网络比BPNN更适用于电离层TEC的短期预报,且在预报因子中加入太阳风数据对电离层正暴的预报效果有明显改善.   相似文献   

13.
Results pertaining to the response of the low latitude ionosphere to a major geomagnetic storm that occurred on 24 August 2005 are presented. The dual frequency GPS data have been analyzed to retrieve vertical total electron content at two Indian low latitude stations (IGS stations) Hyderabad (Geographic latitude 17°20′N, Geographic longitude 78°30′E, Geomagnetic latitude 8.65°N) and Bangalore (Geographic latitude 12°58′N, Geographic longitude 77°33′E, Geomagnetic latitude 4.58°N). These results show variation of GPS derived total electron content (TEC) due to geomagnetic storm effect, local low latitude electrodynamics response to penetration of high latitude convection electric field and effect of modified fountain effect on GPS–TEC in low latitude zone.  相似文献   

14.
Total electron data (TEC) from GPS nowadays can be used as a tool for understanding the space weather phenomena. The development of prediction model for TEC is quiet crucial and challenging due to the dynamic behavior of the ionosphere, since it depends on different factors such as seasonal, diurnal and spatial variations, solar geomagnetic conditions etc. In this paper, an attempt is made for predicting the GPS derived TEC values for different GNSS stations over India using a hybrid method based on Ensemble empirical mode decomposition (EEMD) and Long Short-Term Memory (LSTM) deep learning method. The daily TEC time series data from the IISc Bangalore (Latitude 13.021, Longitude 77.570), Lucknow (Latitude 26.912, Longitude 80.956) and Hyderabad (Latitude 17.417, Longitude 78.551) stations over India during the period 2008 to 2015 of solar cycle 23 and 24 is used for analysis. The assessment of model performance for testing predicted output compared with LSTM and EMD-LSTM models, and their comparison results show that the hybrid EEMD-LSTM model presents better than the other models.  相似文献   

15.
基于加拿大地区高纬度电离层观测网的电离层闪烁观测数据,分析了2018年8月26日地磁暴事件引发的北半球高纬度地区电离层总电子含量(TEC)异常变化、TEC变化率指数(ROTI)及电离层相位闪烁的变化特征.结果表明:加拿大地区最大异常值约6 TECU,磁暴引发全球电离层TEC异常峰值高达20 TECU;加拿大地区电离层相位闪烁发生率最大增至12.6%,而磁静日期间约为1%;强电离层闪烁期间,电离层相位闪烁指数与ROTI之间具有较强的一致性.对GPS双频精密单点定位(Precise Point Positioning,PPP)结果进行分析发现:无闪烁期间定位误差随测站纬度的增高呈现出增大趋势,但均方根误差小于0.4m;闪烁发生期间各测站的定位误差均显著增大,水平和垂直方向定位误差分别增至约0.9m及1.7m.   相似文献   

16.
Nitric Oxide is a very important trace species which plays a significant role acting as a natural thermostat in Earth’s thermosphere during strong geomagnetic activity. In this paper, we present various aspects related to the variation in the NO Infrared radiative flux (IRF) exiting the thermosphere by utilizing the TIMED/SABER (Thermosphere Ionosphere Mesosphere Energetics and Dynamics/ Sounding of the Atmosphere using Broadband Emission Radiometry) observational data during the Halloween storm which occurred in late October 2003. The Halloween storm comprised of three intense-geomagnetic storms. The variability of NO infrared flux during these storm events and its connection to the strength of the geomagnetic storms were found to be different in contrast to similar super storms. The connection between the quantum of energy outflux from the upper atmosphere into space in terms of NO IRF and the duration of storms is established. The NO radiative cooling, and the closely correlated depletion in O/N2 ratio are controlled by the Joule heating intensity (proxied by AE-index). The collisional excitation rate of NO, calculated using the modelled datasets of WACCM-X (Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension), correlates well with the observed pattern of radiative emission by NO. Observational datasets from TIMED/GUVI (Global Ultra-Violet Imager) and MIT Haystack observatory madrigal GNSS (Global navigation satellite system) total electron content (TEC) database shows that the TEC and O/N2 enhancement in low-mid northern hemispheric latitudes are mainly controlled by the z-component of Interplanetary magnetic field (IMF-Bz). The penetration of eastward electric field during the storm events is found to be responsible for the overall enhancement of TEC. The contribution of enhanced day-side TEC in observed variation of O/N2 ratio by GUVI is also reported. It is also seen that during substorms related events the night-time polar region experiences more cooling due to NO than the daytime polar region. The connections between the mid- and low-latitude enhancement in NO IRF with the propagation of LSTIDs (Large-scale traveling ionospheric disturbances) in combination with the O/N2 variability, and the altitudinal variation in NO flux with the progression of the storm is also investigated. This study presents the evidence on the role of diffusion processes in the large scale enhancement of NO in the mesospheric altitudes.  相似文献   

17.
磁暴期间电离层扰动的GPS台网观测分析   总被引:1,自引:3,他引:1  
给出了一种利用GPS台网观测获取TEC快速变化的计算方法,并将该方法用于东亚一澳大利亚扇区的GPS台网观测数据,分析了2000年7月14—18H大磁暴期间的电离层响应,揭示出电离层暴期间赤道异常峰的压缩和移动等特性.计算结果表明,在站点分布不均匀、原始观测数据不足且随时间跳变等多种不利因素的影响下,这种新的算法仍能保持很好的计算稳定性,并能快速地提取给定时空范围内的三维TEC短时变化的特征,适用于研究电离层暴等情况下引起的TEC扰动.  相似文献   

18.
Global Navigation Satellite System (GNSS) measurements of the Total Electron Content (TEC) from local (Dourbes, 50.1°N, 04.6°E) and European IGS (International GNSS Service) stations were used to obtain the TEC changes during the geomagnetic storms of the latest solar activity cycle. A common epoch analysis, with respect to geomagnetic storm intensity, season, and latitude, was performed on data representing nearly 300 storm events. In general, the storm-time behaviour of TEC shows clear positive and negative phases, relative to the non-storm (median) behaviour, with amplitudes that tend to increase during more intense storms. The most pronounced positive phase is observed during winter, while the strongest and yet shortest negative phase is detected during equinox. Average storm-time patterns in the TEC behaviour are deduced for potential use in ionosphere prediction services.  相似文献   

19.
The present study reveals the features of ionospheric parameters variations during the geomagnetic storm of September 7–8, 2017. In particular, parameters of vertical (foF2, foEs) and oblique ionospheric sounding (MOF, modes), absorption level, Total Electron Content (TEC) and particle fluxes at high altitudes were under analysis. The storm was characterized by two Dst-index mimima and can be considered as a sequence of two storms: first - with Dstmin?=??142?nT at 02 UT on September 8th and second - with Dstmin?=??122?nT and at 15 UT on September 8th. It was found that these two storms had different impacts on the ionosphere and HF propagation at mid- and high-latitudes of Northern Hemisphere. The signals of vertical and oblique ionospheric sounding were present in all ionograms before the first storm. Further, at the maximum of the first storm these signals were totally absorbed. Then, before the second storm and during its maximum the signals were detected again in the ionograms due to the low absorption. GOES satellite data showed the significant burst of electrons and protons only during the first storm and small particle fluxes - during the second storm. This feature was also confirmed with GPS data: TEC increased during the first storm and decreased during the second storm.  相似文献   

20.
The equatorial ionosphere has been known to become highly disturbed and thus rendering space-based navigation unreliable during space weather events, such as geomagnetic storms. Modern navigation systems, such as the Global Positioning System (GPS) use radio-wave signals that reflect from or propagate through the ionosphere as a means of determining range or distance. Such systems are vulnerable to effects caused by geomagnetic storms, and their performance can be severely degraded. This paper analyses total electron content (TEC) and the corresponding GPS scintillations using two GPS SCINDA receivers located at Makerere University, Uganda (Lat: 0.3o N; Lon: 32.5o E) and at the University of Nairobi, Kenya (Lat: 1.3o S; Lon: 36.8o E), both in East Africa. The analysis shows that the scintillations actually correspond to plasma bubbles. The occurrence of plasma bubbles at one station was correlated with those at the other station by using observations from the same satellite. It was noted that some bubbles develop at one station and presumably “die off” before reaching the other station. The paper also discusses the effects of the geomagnetic storm of the 24–25 October 2011 on the ionospheric TEC at the two East African stations. Reductions in the diurnal TEC at the two stations during the period of the storm were observed and the TEC depletions observed during that period showed much deeper depletions than on the non-storm days. The effects during the storm have been attributed to the uplift of the ionospheric plasma, which was then transported away from this region by diffusion along magnetic field lines.  相似文献   

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

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