Study of the 2013 Lushan M7.0 earthquake coseismic ionospheric disturbances

On April 20, 2013, an earthquake of M7.0 occurred in Lushan, Sichuan province, China. This paper investigates the coseismic ionospheric anomalies using GPS (Global Positioning System) data from 23 reference stations in Sichuan province that are a part of the Crustal Movement Observation Network of China (CMONOC). The recorded results show that a clear ionospheric anomaly occurred within 15 min after the earthquake near the epicenter, and the occurrence time of the anomalies recorded by various stations is related to the distance from the epicenter. The maximum anomaly is 0.25 TECu, with a 2 min duration and the distance of the recording station to the epicenter is 83 km. Acoustic waves generated by the crustal vertical movement during the earthquake propagate up to the height of the ionosphere lead to the ionospheric anomaly, and the propagation speed of the acoustic wave is calculated as 0.72 ± 0.04 km/s based on the propagation time and propagation distance, consistent with the average speed of sound waves within a 0–450 km atmospheric height.     

Interseismic slip rate of the Garze–Yushu fault belt in the Tibetan Plateau from C-band InSAR observations between 2003 and 2010

InSAR time series techniques can provide high-spatial resolution deformation fields across an active fault belt, even for zones with heavy vegetation coverage. An interseismic deformation map across the Garze–Yushu fault belt in the Tibetan Plateau, ∼300 km by ∼100 km, is derived from C-band Envisat/ASAR imagery collected between 2003 and 2010. Unlike previous research, we obtain a lookup figure which relates the slip rate with the fault locking depth, the dip angle and the rake angle. The estimated slip rate changes significantly with the locking depth and the rake angle, but relatively little with the dip angle. When considering the focal mechanism solutions of historical earthquake along the Garze–Yushu fault, the interseismic slip rate of the Garze–Yushu fault is close to a value of 6.4 mm/yr, which is between the highest (18.2 mm/yr) and the lowest (3.1 mm/yr) slip rate from GPS estimations, but slightly less than the minimum value (∼ 7 mm/yr) from the geological estimations. The earthquake recurrence interval on the Yushu part of Garze–Yushu fault equals 272 yr, and the April 14, 2010 Mw 6.9 earthquake has not completely released the accumulated strain energy between 1738 and 2010.     

GPS total electron content and surface latent heat flux variations before the 11 March 2011 M9.0 Sendai earthquake

The present paper describes the variations of the GPS total electron content (TEC) from the International GNSS service network and surface latent heat flux (SLHF) from the Scientific Computing Division of the National Center for Atmospheric Research (NCAR) before the 11 March 2011 M9.0 Sendai earthquake, respectively. The analysis shows pronounced enhancements in the GPS TEC and SLHF a few days prior to the earthquake event. The maximum increase in the GPS TEC was about 30 TECu with an extended spatial distribution on a geomagnetically quiet day (Dst ? −20 nT, between two moderate geomagnetic storms), 8 March, 3 days prior to the earthquake. This giant positive disturbance was possibly associated with the impending disastrous earthquake and contributed from the enhanced solar radiation. Moreover, there were several anomalous regions of SLHF on the global map, but an area of enhanced SLHF very close to the epicenter. The purpose of this paper is to report the existence of the changes in surface and ionosphere, and show the potential application of multi-source data to identify seismic precursors.     

Optimal combination of InSAR and GPS for measuring interseismic crustal deformation

High spatial resolution measurements of interseismic deformation along major faults are critical for understanding the earthquake cycle and for assessing earthquake hazard. We propose a new remove/filter/restore technique to optimally combine GPS and InSAR data to measure interseismic crustal deformation, considering the spacing of GPS stations in California and the characteristics of interseismic signal and noise using InSAR. To constrain the longer wavelengths (>40 km) we use GPS measurements, combined with a dislocation model, and for the shorter wavelength information we rely on InSAR measurements. Expanding the standard techniques, which use a planar ramp to remove long wavelength error, we use a Gaussian filter technique. Our method has the advantage of increasing the signal-to-noise ratio, controlling the variance of atmosphere error, and being isotropic. Our theoretical analysis indicates this technique can improve the signal-to-noise ratio by up to 20%. We test this method along three segments of the San Andreas Fault (Southern section near Salton Sea, Creeping section near Parkfield and Mojave/Big Bend section near Los Angeles), and find improvements of 26%, 11% and 8% in these areas, respectively. Our data shows a zone of uplift to the west of the Creeping section of the San Andreas Fault and an area of subsidence near the city of Lancaster. This work suggests that after only 5 years of data collection, ALOS interferograms will provide a major improvement in measuring details of interseismic deformation.     

首都圈GPS地形变监测网的布设与观测

系统地介绍了首都圈ＧＰＳ地形变监测网布设的科学思路和原则,该网的基本情况,以及为提高观测精度而采取的重大技术措施。这个ＧＰＳ监测网是全国规模最大、精度最高的监测网,现已完成布设和观测,预期精度将优于５×１０－８将在地震监测预报中发挥重要的作用。     

Ground-based and satellite DC-ULF electric field anomalies around Wenchuan M8.0 earthquake

The ground-based and satellite DC-ULF electric field data were analyzed around Wenchuan M8.0 earthquake on May 12, 2008 in China. The results show that ground electric field anomalies occurred at 3 stations located to the north and south of the epicenter with the amplitude of 3–100 mV/km. The change shapes and their amplitude of ground electric field anomalies are different largely due to their individual underground layer conductivity, water level and so on. The analysis of long time series illustrates that the abnormal geoelectric field started since March 2008. Onboard the DEMETER satellite, the ULF waveforms of electric field were collected and processed by wavelet transform method. The disturbances in the ionosphere were about 3–5 mV/m at a frequency band lower than 0.5 Hz. When the ground and space electric field anomalies were compared, their occurrence time and spatial distribution points are consistent with each other, including the long time anomalies from March 2008 and the short term ones 1–2 days before the Wenchuan earthquake. Finally, the coupling mechanism was discussed.     

Study of seasonal and long-term vertical deformation in Nepal based on GPS and GRACE observations

Lithospheric deformation signal can be detected by combining data from continuous global positioning system (CGPS) and satellite observations from the Gravity Recovery and Climate Experiment (GRACE). In this paper, we use 2.5- to 19-year-long time series from 35 CGPS stations to estimate vertical deformation rates in Nepal, which is located in the southern side of the Himalaya. GPS results were compared with GRACE observations. Principal component analysis was conducted to decompose the time series into three-dimensional principal components (PCs) and spatial eigenvectors. The top three high-order PCs were calculated to correct common mode errors. Both GPS and GRACE observations showed significant seasonal variations. The observed seasonal GPS vertical variations are in good agreement with those from the GRACE-derived results, particularly for changes in surface pressure, non-tidal oceanic mass loading, and hydrologic loading. The GPS-observed rates of vertical deformation obtained for the region suggest both tectonic impact and mass decrease. The rates of vertical crustal deformation were estimated by removing the GRACE-derived hydrological vertical rates from the GPS measurements. Most of the sites located in the southern part of the Main Himalayan Thrust subsided, whereas the northern part mostly showed an uplift. These results may contribute to the understanding of secular vertical crustal deformation in Nepal.     

Seismoionospheric fountain-effect as analogue of active space experiment

Space investigation includes modelling and initiation of natural processes in ionospheric-magnetospheric plasma by active experimental methods with artificial influence as a specific part. For the first time we present here the possibility of triggering the unique equatorial ionosphere phenomenon known as fountain-effect by lithospheric processes before an earthquake. Anomalous influence of earthquake preparation processes with epicenter close to magnetic equator on latitudinal distribution of f_oF - F-region penetration frequencies and spread-F is discussed and considered by us as an analogue of an active experiment. Using Alouette-1 data it is shown that approximately one day before the earthquake f_oF2 dependence on magnetic latitude looks like a double-crest curve with minimum close to epicenter. Simul taneously oblique and spread traces appeared at topside ionograms. The effect is evident during evening hours. Such kind of f_oF2-dependence is unusual and is not observed during the preceding days, nor some days later. So, a possible reason for the observations might be connected with generation of an anomalous electric field near the epicenter during earthquake preparation stage, which initiates a phenomenon similar to the natural fountain-effect.     

A comprehensive study of the 2016 Mw 6.0 Italy earthquake based on high-rate (10 Hz) GPS data

The 2016 M_w 6.0 Italy earthquake is successfully recorded by the near-field 10?Hz GPS and 200?Hz Strong Motion (SM) stations, providing valuable data for this study. A comprehensive study of this earthquake is carried out based on GPS data, which contains coseismic deformations analysis, noise analysis, seismic wave picking, and magnitude determination. The noise of most GPS-derived displacement waveforms can be described as a combination of white noise, flicker noise, and random walk noise after the earthquake occurrence, and the spectral indices vary significantly for most stations, implying that the seismic signals have affected the noise characteristic of GPS-derived displacement waveforms. S-transform is employed to assess the GPS capability to detect the seismic arrival time. The SM station AMT and the GPS station AMAT are in good agreement in seismic wave picking, and the difference is only 1.2?s in the north component, suggesting that the outcome of seismic wave picking using GPS data is reliable. Then, a classic empirical formula is employed to determine the moment magnitude. A robust moment magnitude (Mw 5.90) can be estimated by the nine GPS stations with about 23.9?s. If four GPS stations near the epicenter is chosen to determine the magnitude, it only take 13.0?s to retrieve a reliable preliminary (Mw 5.82) magnitude, which is 5.4?s ahead of nine stations. In addition, Cross Wavelet Transform (XWT) is adopted to measuring the correlation and phase relationship between GPS and SM records. The result of XWT analysis indicates 10?Hz GPS is capable of capturing reliable and accurate coseismic dynamic deformations, as evidenced by the XWT-based semblance being close to 1 between GPS and SM records. The above results confirm the capability of 10?Hz GPS to capture coseismic dynamic deformations, detect seismic arrival time, and determine earthquake magnitude. Moreover, rapid magnitude determination based on 10?Hz GPS data can be regarded as an important supplement to Earthquake Early Warning (EEW).     

Change in refractivity of the atmosphere and large variation in TEC associated with some earthquakes,observed from GPS receiver

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.     

玉树地震前的电离层异常现象分析

分析了玉树地震前地基电离层探测临界频率、GPS TEC和卫星探测原位等离子体参量等多个参数的扰动变化信息, 研究了不同高度异常变化的时空关联性. 分析发现, 在地震前一天的4月13日, 多个电离层参量出现同步扰动异常, 电离层临界频率f₀F₂异常相对滑动中值增大40%, 异常空间上存在从震中东南向西南漂移的特性; GPS TEC异常增强15TECU (1TECU=1016m-2)左右, 分布于震中南部经度15°范围内, 且有明显的磁共轭效应; DEMETER观测的原位氧离子密度N_i(O+) 4月13日为1-4月中最强的一天, 异常分布偏向赤道区, 但仅局限在30°-50°左右的经度范围内. 综合三个参量的异常特征发现, 无论是空间的局地性还是时间上的密切关联均反映这次电离层扰动可能与玉树地震孕育有关. 结合其他观测信息, 进一步探讨了这次地震孕育过程的地震电离层耦合机理.      

Present-day crustal movement of the Chinese mainland based on Global Navigation Satellite System data from 1998 to 2018

A total of more than 260 continuous stations and 2000 campaign stations from the Crustal Movement Observation Network of China (CMONOC) project, covering the Chinese mainland and its surrounding areas during the period of 1998–2018, are processed using the Bernese Global Navigation Satellite System (GNSS) software via a state-of-the-art method. We obtain the coordinate time series of all the stations given in the reference frame ITRF2014, estimate the coseismic deformation, and remove outliers. Lastly, we present the latest, most complete, and most accurate contemporary horizontal velocity field with respect to the stable Eurasian plate, irrespective of the postseismic deformations. This study shows that the signal of tectonic movement in Western China is stronger than that in Eastern China particularly in the Tibetan Plateau, with a rate of 18–32?mm/a. Moreover, the signal decays sharply from south to north. However, North China and South China move coherently to the ESE direction mostly at a rate of 4–10?mm/a and have not experienced any abrupt velocity gradients in their interiors. Meanwhile, Northeast China has the lowest velocity of only 2–4?mm/a in addition to the coastal areas that have slightly larger velocities. The densified and continuous observation of GNSS stations are of great significance to the study of the present-day crustal movement and tectonic deformation characteristics of the Chinese mainland. This would help to provide better constraints on the kinematics and dynamics of the region.     

Noise analysis of continuous GPS coordinate time series for CMONOC

The Crustal Movement Observation Network of China (CMONOC) is one of the major scientific infrastructures, mainly using Global Positioning System (GPS) measurements, to monitor crustal deformation in the Chinese mainland. In this paper, decade-long coordinate time series of 26 continuous GPS sites of CMONOC are analyzed for their noise content using maximum likelihood estimation (MLE). We study the noise properties of continuous GPS time series of CMONOC for the unfiltered, filtered solutions and also the common mode signals in terms of power law plus white noise model. In the spatial filtering, we remove for every time series a common mode error that was estimated from a modified stacking of position residuals from other sites within ∼1000 km of the selected site. We find that the common mode signal in our network has a combination of spatially correlated flicker noise and a common white noise with large spatial extent. We demonstrate that for the unfiltered solutions of CMONOC continuous GPS sites the main colored noise is a flicker process, with a mean spectral index of ∼1. For the filtered solutions, the main colored noise is a general power law process, indicating that a major number of the filtered regional solutions have a combination of noise sources or local effects. The velocity uncertainties from CMONOC continuous GPS coordinate time series may be underestimated by factors of 8–16 if a pure white noise model is assumed. In addition, using a white plus flicker noise model, the median values of velocity errors for the unfiltered solutions are 0.16 (north), 0.17 (east) and 0.58 (vertical) mm/yr, and the median values for the filtered solutions are 0.09 (north), 0.10 (east) and 0.40 (vertical) mm/yr.     

Pre-earthquake ionospheric anomalies before three major earthquakes by GPS-TEC and GIM-TEC data during 2015–2017

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.     

Determination of crustal deformations based on GPS observing-session duration in Marmara region,Turkey

This paper reports on the estimation of changes of a GPS-based crustal deformation field as function of the observing session duration. This investigation was carried out on the Marmara Continuous GPS Network (Marmara region) by using GPS data collected on a set of 10 stations during the period 2002.50–2005.47 The GPS observations were processed in the ITRF 2005 reference frame using Bernese v5.0 software.     

The Western Crete geodetic infrastructure: Long-range power-law correlations in GPS time series using Detrended Fluctuation Analysis

The Global Positioning System makes it possible, nowadays, to measure crustal displacements with unprecedented accuracy. These measurements can improve our understanding of the spatio-temporal evolution of motion along tectonic plate boundaries, as well as deepen our comprehension for the nature of earthquake fault behavior and earthquake cycle. The Hellenic subduction zone and the earthquake-prone Crete, with large crustal motions of 3–4 cm/yr, may provide such a natural laboratory for this type of investigation. In this work, we examine the statistical structure and behavior of time series, as produced by permanent GPS sites, established on the islands of Crete and Gavdos.     

Anomalous enhancement of electric field derived from ionosonde data before the great Wenchuan earthquake

The seismo-ionospheric coupling is a hot topic of modern scientific research. One of the most reasonable mechanisms is the hypothesis of quasistatic electric field. With a preliminary analytical method, electric field in F2-region before Wenchuan earthquake was derived from the five lower latitude ground-based ionosondes: Lhasa (29.63°N, 91.17°E), Chongqing (29.50°N, 106.40°E), Kunming (25.00°N, 102.70°E), Guangzhou (23.15°N, 113.35°E) and Haikou (23.15°N, 110.35°E). The results show that there was an anomalous enhancement of electric field in F2-region close to the epicenter with maximal amplitude ∼2 mV/m (about 10 times of the background) at 15:00 LT on 9 May, which is in line with the amplitudes of anomalous electric field related to some previous earthquakes observed by satellites.     

Application of SAR data to monitoring earth surface changes and displacement

Surface changes and displacement caused by earthquake and volcanic activities can be detected by satellite radar interferometric technology using L-band JERS-1 SAR images. For the detection of the coseismic deformation of the Hyogoken-nanbu earthquake, interferometric analysis shows 0.6 meter surface displacement along the radar line-of-sight direction around Kobe City and 1.1 meter displacement in Awaji Island where epicenter is located. For the detection of volcanic body deformation in Iwo-jima, the analysis shows concentric circular subsiding displacement that agrees with the results of ground measurement.     

Seismic deformation of the Mw 8.0 Wenchuan earthquake from high-rate GPS observations

High-rate GPS positioning has been recognized as a powerful tool in estimating epoch-wise station displacement which is particularly useful for seismology. In this study, station displacements during the 12 May 2008 M_w 8.0 Wenchuan earthquake are derived from the 1-Hz GPS data collected at a set of stations in China. The impacts of integer ambiguity resolution and station environment-dependent effects are investigated in order to yield more accurate results. The position accuracy of horizontal components of better than 1 cm suggests that GPS can sense the rapid position oscillation of about 2 cm in amplitude. Temporal and spatial analysis is applied to the surface displacement at station XANY and the characteristics of the movements due to Rayleigh and Love waves are detected and discussed. The comparison of GPS-derived displacement with relevant synthetic data computed based on a recently published rapture model shows a reasonable agreement in waveform. The various differences in amplitude need further investigation and also imply that rapture inversion might be improved if GPS-derived displacement is assimilated.     

Influence of different data fusion methods on the accuracy of three-dimensional displacements field

Combining multi-source data can improve the accuracy and the spatial resolution of the three-dimensional (3-D) displacements field. How to effectively integrate multi-source data to obtain high-precision and high spatial resolution 3-D displacements field is worthy of further study. The stochastic model and fusion model of integrating multi-source data affect the accuracy of data fusion. In this paper, based on the least squares method, the effects of different stochastic models and data fusion models on the 3-D displacements field’s accuracy are studied. The optimal method for estimating large-scale 3-D displacements field from integrated InSAR, leveling and GPS measurements is obtained. Then we realize the integrating InSAR, leveling and GPS measurements to obtain the high-precision 3-D displacements velocity field in Tianjin (China) from 2016 to 2018. The results are validated with GPS measurements at 6 independent stations, with the root mean squares (RMS) residuals of the discrepancies being 2.39 mm/yr, 2.54 mm/yr and 2.83 mm/yr in eastern, northern and vertical directions, respectively. By comparing different stochastic models, the 3-D displacements field obtained from multi-source data is optimized by the variance component estimation-least squares method, which is better than weighted least squares (WLS) method. By comparing different data fusion models, the accuracy of the horizontal displacements velocity is better than that of interpolated GPS results. The horizontal displacements component has a great influence on the vertical displacements velocity accuracy in the process of acquiring the 3-D displacements velocity by integrating InSAR, GPS and leveling measurements. This study provides a reference method for integrating multi-source data to obtain 3-D displacements field. This method effectively utilizes the advantages of GPS, InSAR and leveling measurements, and extends the limitations of single technical in describing surface-time scale applications. The 3-D displacements information with a large spatial scale and high spatial resolution provide a reliable data basis for studying the crustal movement and its dynamic mechanism.