首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 390 毫秒
1.
We investigate the ionospheric total electron content (TEC) anomalies occurred in the Qinghai-Tibet region before three large earthquakes (M > 7.0). The temporal and spatial TEC variations were used to detect the ionospheric possible precursors of these earthquakes. We identified two TEC enhancements in the afternoon local time 9 days and 2–3 days before each earthquake, between which a TEC decrement occurred 3–6 days before earthquakes. These anomalies happened in the area of about 30° in latitude and the maximum is localized equatorward from the epicenters. These TEC anomalies can be found in all three earthquakes regardless the geomagnetic conditions. The features of these anomalies have something in common and may have differences from those caused by geomagnetic storms. Our results suggest that these ionospheric TEC perturbations may be precursors of the large earthquakes.  相似文献   

2.
This paper investigates the capacity of the latest version of the International Reference Ionosphere (IRI-2012) model in predicting the vertical Total Electron Content (vTEC) over Ethiopian regions during solar minimum (2009) and solar maximum (2013) phases. This has been carried out by comparing the IRI-2012 modeled and experimental vTEC inferred from eight ground based dual frequency GPS (Global Positioning System) receivers installed recently at different regions of the country. In this work, the diurnal, monthly and seasonal variation in the measured vTEC have been analyzed and compared with the IRI-2012 modeled vTEC. During the solar minimum phase, the lowest and highest diurnal peak of the experimental vTEC are observed in July and October, respectively. In general, the diurnal variability of vTEC has shown minimum values around 0300 UT (0600 LT) and maximum values between around 1000 and 1300 UT (1300 and 1600 LT) during both solar activity phases. Moreover, the maximum and minimum monthly and seasonal mean hourly vTEC values are observed in October and July and in the March equinox and June solstice, respectively. It is also shown that the IRI-2012-model better predicts the diurnal vTEC in the time interval of about 0000–0300 UT (0300–0600 LT) during the solar minimum phase. However, the model generally overestimates the diurnal vTEC except in the time interval of about 0900–1500 UT (1200–1800 LT) during the solar maximum phase. The overall result of this work shows that the diurnal vTEC prediction performance of the model is generally better during the solar minimum phase than during solar maximum phase. Regarding the monthly and seasonal prediction capacity of the model, there is a good agreement between the modeled and measured monthly and seasonal mean hourly vTEC values in January and December solstice, respectively. Another result of the work depicts that unlike the GPS–TEC the IRI-2012 TEC does not respond to the effect resulted from geomagnetic storms.  相似文献   

3.
Observations of ionospheric vertical total electron content (vTEC) from European ground-based Global Navigation Satellite Systems (GNSS) receivers during the period January 2008–January 2010 are used to investigate, for the first time, vTEC sensitivity to weak geomagnetic disturbances under extreme solar minimum conditions. This study shows a significant number of events for the period in question, all of which exhibited some form of exceptionally large values of vTEC during small-magnitude geomagnetic disturbances. To illustrate our point on the importance of vTEC enhancements during the extreme solar minimum and its relevance for the current GNSS and future Galileo applications, we present in this paper the results associated with two significant events that both occurred in equinoctial months. The 10–12 October 2009 event of anomalous TEC enhancement at two distant mid-latitude locations HERS (0.3 E; 50.9 N) and NICO (33.4 E, 35.1 N) is discussed in the context of strong vTEC variations during the well established ionospheric storm on 11 October 2008. We conclude with a short summary of the new findings and their consequences on ionospheric monitoring and modelling for operational communication and navigation systems.  相似文献   

4.
Employing a dual frequency GPS-receiver, ionospheric total electron (TEC) measurements have been in progress at Agra station in India (Geograph. lat. 27.2°N, long. 78°E) since 1 April 2006. In this paper, the TEC data have been analyzed for a period of one month from 1 April-1 May 2013 to examine the effect of multiple earthquakes, some of which occurred on the same day of 16 April 2013, and others occurred in the same month of April, 2013 in India and neighboring countries. We process the data using quartile and epoch analysis based statistical techniques and show that out of all the earthquakes, the one of the largest magnitude (M = 7.8) that occurred on Pakistan-Iran border caused anomalous enhancements and depletions in TEC 1–9 days before the occurrence of main sock. The E × B drift mechanism is suggested for the anomalies to occur in which the seismogenic electric field E is generated in a process suggested by Pulinets (2004).  相似文献   

5.
A precise determination of ionospheric total electron content (TEC) anomaly variations that are likely associated with large earthquakes as observed by global positioning system (GPS) requires the elimination of the ionospheric effect from irregular solar electromagnetic radiation. In particular, revealing the seismo-ionospheric anomalies when earthquakes occurred during periods of high solar activity is of utmost importance. To overcome this constraint, a multiresolution time series processing technique based on wavelet transform applicable to global ionosphere map (GIM) TEC data was used to remove the nonlinear effect from solar radiation for the earthquake that struck Tohoku, Japan, on 11 March, 2011. As a result, it was found that the extracted TEC have a good correlation with the measured solar extreme ultraviolet flux in 26–34 nm (EUV26–34) and the 10.7 cm solar radio flux (F10.7). After removing the influence of solar radiation origin in GIM TEC, the analysis results show that the TEC around the forthcoming epicenter and its conjugate were significantly enhanced in the afternoon period of 8 March 2011, 3 days before the earthquake. The spatial distributions of the TEC anomalous and extreme enhancements indicate that the earthquake preparation process had brought with a TEC anomaly area of size approximately 1650 and 5700 km in the latitudinal and longitudinal directions, respectively.  相似文献   

6.
Earthquake prediction stimulates the searches for a correlation between seismic activity and ionospheric anomalies. Contrary to common focuses on strong earthquakes, we report the ionospheric disturbances, 2 days before a moderate Ms = 4.7 Chongqing earthquake (29.4°N, 105.5°E, depth = 7.0 km, occurred at 21:21 LT, 10 September, 2010) with the data of ground-based ionosondes and IGS receivers. The data covering the period under the quiet geomagnetic conditions and a geomagnetic storm was analyzed with upper and lower bounds. It is found that there were significant enhancements of foF2 and total electron content (TEC) on the afternoon of 8 September, 2010, with a limited area close to the epicentre, which was different from the feature of ionospheric perturbations triggered by the geomagnetic storm on 15 September. Taking into account the heliogeomagnetical condition, we conclude that the observed ionospheric enhancements were very likely associated with the forthcoming moderate Chongqing earthquake, which implies that the relationship between the amplitudes of ionospheric disturbances and earthquakes is very complicated.  相似文献   

7.
Total electron content (TEC) measured simultaneously using Global Positioning System (GPS) ionospheric monitors installed at some locations in Nigeria during the year 2011 (Rz = 55.7) was used to study the diurnal, seasonal, and annual TEC variations. The TEC exhibits daytime maximum, seasonal variation and semiannual variations. Measured TEC were compared with those predicted by the improved versions of the International Reference Ionosphere (IRI) and NeQuick models. The models followed the diurnal and seasonal variation patterns of the observed values of TEC. However, IRI model produced better estimates of TEC than NeQuick at all locations.  相似文献   

8.
A powerful earthquake of Mw = 7.7 struck the Saravan region (28.107° N, 62.053° E) in Iran on 16 April 2013. Up to now nomination of an automated anomaly detection method in a non linear time series of earthquake precursor has been an attractive and challenging task. Artificial Neural Network (ANN) and Particle Swarm Optimization (PSO) have revealed strong potentials in accurate time series prediction. This paper presents the first study of an integration of ANN and PSO method in the research of earthquake precursors to detect the unusual variations of the thermal and total electron content (TEC) seismo-ionospheric anomalies induced by the strong earthquake of Saravan. In this study, to overcome the stagnation in local minimum during the ANN training, PSO as an optimization method is used instead of traditional algorithms for training the ANN method. The proposed hybrid method detected a considerable number of anomalies 4 and 8 days preceding the earthquake. Since, in this case study, ionospheric TEC anomalies induced by seismic activity is confused with background fluctuations due to solar activity, a multi-resolution time series processing technique based on wavelet transform has been applied on TEC signal variations. In view of the fact that the accordance in the final results deduced from some robust methods is a convincing indication for the efficiency of the method, therefore the detected thermal and TEC anomalies using the ANN + PSO method were compared to the results with regard to the observed anomalies by implementing the mean, median, Wavelet, Kalman filter, Auto-Regressive Integrated Moving Average (ARIMA), Support Vector Machine (SVM) and Genetic Algorithm (GA) methods. The results indicate that the ANN + PSO method is quite promising and deserves serious attention as a new tool for thermal and TEC seismo anomalies detection.  相似文献   

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

10.
Earthquake (EQ) anomalies in the form of enhancement and depletion in ionospheric Total Electron Content (TEC) from Global Positioning System (GPS) may considerably alarm about short and long term precursors of the impending main shock. In this paper, TEC anomalies are investigated from permanent GPS ground-stations in Turkey associated to Mw ≥ 6.0 EQs occurred in 2011–2012. Temporal and spatial analyses of TEC at 2 h sampling have shown significant evidences about EQ induced ionospheric anomalies during 10–14 h of UT (Universal Time) within 5 days before Mw 6.0 Greece, and Mw 7.1, Turkish EQ. Spatial analyses have manifested arrival of TEC anomalies at UT = 10 h to epicenter of both EQs, which linger above epicenter during UT = 12–14 h and left seismogenic zone after UT = 14 h before every EQ during Kp < 3 and Dst = 0 nT. Meanwhile, a geomagnetic storm (Dst < -100 nT) induce perturbation two days after the Mw 7.1 Turkish EQ, showing no relation with epicenter during spatial analysis. It also shows that TEC can be useful to distinguish geomagnetic storm variations to successfully detect EQ precursors. These anomalies during quiet storm (Kp < 3; Dst = 0 nT) conditions may be effective to link the lithosphere and ionosphere in severe seismic zones to detect EQ precursors before future EQs. Interpretation of TEC anomalies and it enhancements over EQ epicenters during UT = 12–14 h for both EQs have shown that EQs anomalies only occurred in particular time. Whereas, geomagnetic storm effect occurred during whole abnormal day over the Earth.  相似文献   

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

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

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

14.
We investigated the ionospheric anomalies observed before the Tohoku earthquake, which occurred near the northeast coast of Honshu, Japan on 11 March, 2011. Based on data from a ground-based Global Positioning System (GPS) network on the Korean Peninsula, ionospheric anomalies were detected in the total electron content (TEC) during the daytime a few days before earthquake. Ionospheric TEC anomalies appeared on 5, 8 and 11 March. In particular, the ionospheric disturbances on 8 March evidenced a remarkable increase in TEC. The GPS TEC variation associated with the Tohoku earthquake was an increase of approximately 20 total electron content units (TECU), observed simultaneously in local and global TEC measurements. To investigate these pre-earthquake ionospheric anomalies, space weather conditions such as the solar activity index (F10.7) and geomagnetic activity indices (the Kp and Dst indices) were examined. We also created two-dimensional TEC maps to visual the spatial variations in the ionospheric anomalies preceding the earthquake.  相似文献   

15.
This paper presents the vertical total electron content vTEC variations for three African stations, located at mid-low and equatorial latitudes, and operating since more than 10 years. The vTEC of the middle latitude GPS station in Alexandria, Egypt (31.2167°N; 29.9667°E, geographic) is compared to the vTEC of two others GPS stations: the first one in Rabat/Morocco (33.9981°N; 353.1457°E, geographic), and the second in Libreville/Gabon (0.3539°N; 9.6721°E, geographic). Our results discussed the diurnal, seasonal, and solar cycle dependences of vTEC at the local ionospheric conditions, during different phases of solar cycle in the light of the classification of Legrand and Simon. The vTEC over Alexandria exhibits the well-known equinoctial asymmetry which changes with the phases of the solar cycle; the spring vTEC is larger than that of autumn during the maximum, decreasing and minimum phases of solar cycle 23. During the increasing phase of solar cycle 24, it is the contrary. The diurnal variation of the vTEC presents multiple maxima during the equinox from 2005 to 2008 and during the summer solstice from 2006 to 2012. A nighttime vTEC enhancement and winter anomaly are also observed. During the deep solar minimum (2006–2009) the diurnal variation of the vTEC observed over Alexandria is similar to the diurnal variation observed during quiet magnetic period at equatorial latitudes. We observed also that the amplitude of vTEC at Libreville is larger than the amplitude of vTEC observed at Alexandria and Rabat, indeed Libreville is near the southern crest of the Equatorial Ionization anomaly. Finally, the correlation coefficient between vTEC and the sunspot number Rz is high and changes with solar cycle phases.  相似文献   

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

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

18.
2017年9月8日发生了一次强磁暴,Kp指数最大值达到8.利用区域电离层格网模型(Regional Ionosphere Map,RIM)和区域ROTI(Rate of TEC Index)地图,分析了磁暴期间中国及其周边地区电离层TEC扰动特征和低纬地区电离层不规则体的产生与发展情况,同时利用不同纬度IGS(International GNSS Service)测站BJFS(39.6°N,115.9°E),JFNG(30.5°N,114.5°E)和HKWS(22.4°N,114.3°E)的GPS双频观测值,获取各测站的ROTI和DROT(Standard Deviation of Differential ROT)指数变化趋势.结果表明:此次磁暴发生期间电离层扰动先以正相扰动为主,主要发生在中低纬区域,dTEC(differential TEC)最大值达到14.9TECU,随后电离层正相扰动逐渐衰减,在低纬区域发生电离层负相扰动,dTEC最小值达到-7.2TECU;在12:30UT-13:30UT时段,中国南部低纬地区发生明显的电离层不规则体事件;相比BJFS和JFNG两个测站,位于低纬的HKWS测站的ROTI和DROT指数变化更为剧烈,这表明电离层不规则体结构存在纬度差异.   相似文献   

19.
Variations of the ionospheric F2 region critical frequency (foF2) have been investigated statistically before earthquakes during 1980–2008 periods in Japan area. Ionosonde data was taken from Kokubunji station which is in the earthquake preparation zone for all earthquakes. Standard Deviations and Inter-Quartile Range methods are applied to the foF2 data. It is observed that there are anomalous variations in foF2 before earthquakes. These variations can be regarded as ionospheric precursors and may be used for earthquake prediction.  相似文献   

20.
The ionospheric total electron content (TEC), derived by analyzing dual frequency signals from the Global Positioning System (GPS) recorded near the Indian equatorial anomaly region, Varanasi (geomagnetic latitude 14°, 55′N, geomagnetic longitude 154°E) is studied. Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the solar minimum period from May 2007 to April 2008. It is found that the daily maximum TEC near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semiannual variation is seen with two maxima occurring in both spring and autumn. Statistical studies indicate that the variation of EIA crest in TEC is poorly correlated with Dst-index (r = −0.03) but correlated well with Kp-index (r = 0.82). The EIA crest in TEC is found to be more developed around 12:30 LT.  相似文献   

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

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