Co-seismic deformation for the 2015 M_(W)7.8 Gorkha earthquake(Nepal)using near-field GPS data

Seasonal variations and common mode errors affect the precision of the Global Positioning System(GPS)time series.In this paper,we explore to improve the precision of coordinate time series,thereby providing a better detection of weak or transient deformation signals,particularly co-seismic signals.Based on 97 GPS stations,including the campaign and continuous GPS stations in Nepal and southern Tibet,we first consider seasonal variations and common errors,then obtain co-seismic deformation of the 2015 Gorkha earthquake in Nepal and southern Tibet.Our co-seismic rupture model is characterized by a shallow ramp and a deeper detachment fault,in agreement with the relocated aftershock sequence.Our results indicate that the earthquake rupture is mainly distributed in the upper-crustal fault,and the maximum slip is up to 8.0 m at~15.0 km depth located in the approximate-80 km east of the epicenter.The average slip is more than 5 m,and the total modelled magnitude is M_(W)7.84,consistent with the observed seismic moment.Our rupture model for the 2015 Gorkha earthquake suggests that the rupture zone is not only in the upper crustal Main Himalayan Thrust(MHT),but also spreads to the northern segment of the MHT.

Analysis of the Wenchuan Ms8.0 Earthquake's Co-seismic Stress and Displacement Change by Using the Finite Element Method

The variation of in situ stress before and after earthquakes is an issue studied by geologists. In this paper, on the basis of the fault slip dislocation model of Wenchuan Ms8.0 earthquake, the changes of co-seismic displacement and the distribution functions of stress tensor around the Longmen Shan fault zone are calculated. The results show that the co-seismic maximum surface displacement is 4.9 m in the horizontal direction and 6.5 m in the vertical direction, which is almost consistent with the on-site survey and GPS observations. The co-seismic maximum horizontal stress in the hanging wall and footwall decreased sharply as the distance from the Longmen Shan fault zone increased. However, the vertical stress and minimum horizontal stress increased in the footwall and in some areas of the hanging wall. The study of the co-seismic displacement and stress was mainly focused on the long and narrow region along the Longmen Shan fault zone, which coincides with the distribution of the earthquake aftershocks. Therefore, the co-seismic stress only affects the aftershocks, and does not affect distant faults and seismic activities. The results are almost consistent with in situ stress measurements at the two sites before and after Wenchuan Ms8.0 earthquake. Along the fault plane, the co-seismic shear stress in the dip direction is larger than that in the strike direction, which indicates that the faulting mechanism of the Longmen Shan fault zone is a dominant thrust with minor strike-slipping. The results can be used as a reference value for future studies of earthquake mechanisms.

Surface Rupture and Co-seismic Displacement Produced by the Ms 8.0 Wenchuan Earthquake of May ^(12)th,2008,Sichuan,China:Eastwards Growth of the Qinghai-Tibet Plateau

An earthquake of Ms 8 struck Wenchuan County, western Sichuan, China, on May 12^th, 2008 and resulted in long surface ruptures （〉300 km）. The first-hand observations about the surface ruptures produced by the earthquake in the worst-hit areas of Yingxiu, Beichuan and Qingchuan, ascertained that the causative structure of the earthquake was in the central fault zones of the Longmenshan tectonic belt. Average co-seismic vertical displacements along the individual fault of the Yingxiu-Beiehuan rupture zone reach 2.514 m and the cumulative vertical displacements across the central and frontal Longmenshan fault belt is about 5-6 m. The surface rupture strength was reduced from north of Beichuan to Qingchuan County and shows 2-3 m dextral strike-slip component. The Wenchuan thrust-faulting earthquake is a manifestation of eastward growth of the Tibetan Plateau under the action of continuous convergence of the Indian and Eurasian continents.

Co-seismic slip distribution of the 2011 Tohoku(MW 9.0)earthquake inverted from GPS and space-borne gravimetric data

Data obtained by GRACE(Gravity Recovery and Climate Experiment) have been used to invert for the seismic source parameters of megathrust earthquakes under the assumption of either uniform slip over an entire fault or a point-like seismic source.Herein, we further extend the inversion of GRACE long-wavelength gravity changes to heterogeneous slip distributions during the 2011 Tohoku earthquake using three fault models:(Ⅰ) a constant-strike and constant-dip fault,(Ⅱ) a variable dip fault, and(Ⅲ) a realistically varying strike fault. By removing the post-seismic signal from the time series, and taking the effect of ocean water redistribution into account, we invert for slip models I, II, and III using co-seismic gravity changes measured by GRACE, de-striped by DDK3 decorrelation filter. The total seismic moments of our slip models, with respective values of 4.9×10^(22) Nm, 5.1×10^(22) Nm, and 5.0×10^(22) Nm, are smaller than those obtained by other studies relying on GRACE data. The resulting centroids are also located at greater depths(20 km, 19.8 km,and 17.4 km, respectively). By combining onshore GPS, GPS-Acoustic, and GRACE data, we obtain a jointly inverted slip model with a seismic moment of 4.8×10^(22) Nm, which is larger than the seismic moment obtained using only the GPS displacements. We show that the slip inverted from low degree space-borne gravimetric data, which contains information at the ocean region, is affected by the strike of the arcuate trench. The space-borne gravimetric data help us constrain the source parameters of a megathrust earthquake within the frame of heterogeneous slip models.

Co-seismic Faults and Geological Hazards and Incidence of Active Fault of Wenchuan Ms 8.0 Earthquake,Sichuan,China

There are two co-seismic faults which developed when the Wenchuan earthquake happened. One occurred along the active fault zone in the central Longmen Mts. and the other in the front of Longmen Mts. The length of which is more than 270 kin and about 80 km respectively. The co-seismic fault shows a reverse flexure belt with strike of N45°-60°E in the ground, which caused uplift at its northwest side and subsidence at the southeast. The fault face dips to the northwest with a dip angle ranging from 50° to 60°. The vertical offset of the co-seismic fault ranges 2.5-3.0 m along the Yingxiu- Beichuan co-seismic fault, and 1.5-1.1 m along the Doujiangyan-Hanwang fault. Movement of the coseismic fault presents obvious segmented features along the active fault zone in central Longmen Mts. For instance, in the section from Yingxiu to Leigu town, thrust without evident slip occurred; while from Beichuan to Qingchuan, thrust and dextral strike-slip take place. Main movement along the front Longmen Mts. shows thrust without slip and segmented features. The area of earthquake intensity more than IX degree and the distribution of secondary geological hazards occurred along the hanging wall of co-seismic faults, and were consistent with the area of aftershock, and its width is less than 40km from co-seismic faults in the hanging wall. The secondary geological hazards, collapses, landslides, debris flows et al., concentrated in the hanging wall of co-seismic fault within 0-20 km from co-seismic fault.

Co-seismic displacements of 2011 Japan Mw9.0 earthquake recorded by far-field GPS stations

Co-seismic displacements of the 2011 Mw9.0 Japan earthquake recorded by GPS stations in China and surrounding areas showed a movement toward the epicenter. The horizontal displacements were up to 1 - 3 cm in northeastern China, 3 -8 mm in the North China, and 2 cm in the Korean peninsula. The vertical movements in China were small uplifts.

Deformation of the Most Recent Co-seismic Surface Ruptures Along the Garzê–Yushu Fault Zone(Dangjiang Segment)and Tectonic Implications For the Tibetan Plateau

The Garzê–Yushu strike-slip fault in central Tibet is the locus of strong earthquakes（M 〉 7）. The deformation and geometry of the co-seismic surface ruptures are reflected in the surface morphology of the fault and depend on the structure of the upper crust as well as the pre-existing tectonics. Therefore, the most recent co-seismic surface ruptures along the Garzê–Yushu fault zone（Dangjiang segment） reveal the surface deformation of the central Tibetan Plateau. Remote sensing images and field investigations suggest a 85 km long surface rupture zone（striking NW-NWW）, less than 50 m wide, defined by discontinuous fault scarps, right-stepping en echelon tensional cracks and left-stepping mole tracks that point to a left-lateral strike-slip fault. The gullies that cross fault scarps record systematic left-lateral offsets of 1.8 m to 5.0 m owing to the most recent earthquake, with moment magnitude of about M 7.5, in the Dangjiang segment. Geological and geomorphological features suggest that the spatial distribution of the 1738 co-seismic surface rupture zone was controlled by the pre-existing active Garzê–Yushu fault zone（Dangjiang segment）. We confirm that the Garzê–Yushu fault zone, a boundary between the Bayan Har Block to the north and the Qiangtang Block to the south, accommodates the eastward extrusion of the Tibetan Plateau and generates strong earthquakes that release the strain energy owing to the relative motion between the Bayan Har and Qiangtang Blocks.

Multi-scale Decomposition of Co-seismic Deformation from High Resolution DEMs:a Case Study of the 2004 Mid-Niigata Earthquake

Decomposing co-seismic deformation is an immediate need for researchers who are interested in earthquake inversion analysis and geo-hazard mapping. However, conventional InSAR or digital elevation models （DEMs） imagery analyses only provide the displacement in the Line-of-Sight （LOS） direction or elevation changes. The 2004 Mid-Niigata earthquake in Japan provides lessons on how to decompose co-seismic deformation from two sets of DEMs. If three adjacent points undergo a rigid-body-translation movement, their co-seismic deformation can be decomposed by solving simultaneous equations. Although this method has been successfully used to discuss tectonic deformations, the algorithm needed improvement and a more rigorous algorithm, including a new definition of nominal plane, DEMs comparability improvement and matrix condition check is provided. Even with these procedures, the obtained decomposed displacement often showed remarkable scatter prompting the use of the moving average method, which was used to determine both tectonic and localized displacement characteristics. A cut-off window and a pair of band-pass windows were selected according to the regional geology and construction activities to ease the tectonic and localized displacement calculations, respectively. The displacement field of the tectonic scale shows two major clusters of large lateral components, and coincidently major visible landslides were found mostly within them. The localized displacement helps to reveal hidden landslides in the target area. As far as the Kizawa hamlet is concerned, the obtained vectors show down-slope movements, which are consistent with the observed traces of dislocations that were found in the Kizawa tunnel and irrigation wells. The method proposed has great potential to be applied to understanding post-earthquake rehabilitation in other areas.

Pre-seismic and Co-seismic Crustal Movements of the M 7.3 Kyushu Earthquake on April 16,2016 in Japan

The time series of coordinates of a large number of GPS stations in the world,processed by Prof. Geoffrey Blewitt with GIPSY software are available at http://geodesy. unr. edu.Based on the time series of coordinates in the global reference frame of IGS08 at more than250 stations of continuous GPS observations,downloaded from the website,the co-seismic displacements of the M7. 3 Kyushu earthquake on April 16,2016 in Japan and the preseismic strain accumulations and displacements in the regional reference frame were obtained. The station of continuous GPS observation at BJFS near Beijing has been quite stable in displacement in the eastern part of China for more than 17 years since the beginning of its operation,and this station is used as the core station in the regional reference frame for the pre-seismic displacement of the Kyushu earthquake of M7. 3. The main feature of the pre-seismic displacements of the Kyushu earthquake is characterized by locking in the crust at and near the epicenter. The anomalous pre-seismic strain accumulation developed in an area of anomalous accumulation of the shear strain component of γ1 on the northeast side of the epicenter,with increasing size of the area and increasing magnitude in γ1. The largest area covered by the anomalous γ1 is about 2000 km2. The change in the E component at BJFS since November 26,2015 was caused by the replacement of the receiver and the antenna at the station. In order to study the shortterm change in displacements at stations at and near the epicenter,the time series at 3 stations with continuous GPS observations,2 at SUWN and DAEJ in south Korea and 1 at BJSH near Beijing were analyzed. The analysis shows that the displacements at the 3 stations have been quite stable in the same manner in east Asia. Thus,BJSH is used as the core station in the regional reference frame of displacement and the displacement time series show that there were no significant short term anomalies before the earthquake.

Numerical simulation of co-seismic deformation of 2011 Japan Mw9.0 earthquake

Co-seismic displacements associated with the Mw9.0 earthquake on March 11, 2011 in Japan are numerically simulated on the basis of a finite-fault dislocation model with PSGRN/PSCMP software. Compared with the inland GPS observation, 90% of the computed eastward, northward and vertical displacements have residuals less than 0.10 m, suggesting that the simulated results can be, to certain extent, used to demon- strate the co-seismic deformation in the near field. In this model, the maximum eastward displacement increa- ses from 6 m along the coast to 30 m near the epicenter, where the maximum southward displacement is 13 m. The three-dimensional display shows that the vertical displacement reaches a maximum uplift of 14.3 m, which is comparable to the tsunami height in the near-trench region. The maximum subsidence is 5.3 m.

Decomposing InSAR LOS displacement into co-seismic dislocation with a linear in-terpolation model: A case study of the Kunlun Mountain M_s=8.1 earthquake

It has always been a difficult problem to extract horizontal and vertical displacement components from the InSAR LOS （Line of Sight） displacement since the advent of monitoring ground surface deformation with InSAR technique. Having tried to fit the firsthand field investigation data with a least squares model and obtained a preliminary result, this paper, based on the previous field data and the InSAR data, presents a linear cubic interpolation model which well fits the feature of earthquake fracture zone. This model inherits the precision of investigation data; moreover make use of some advantages of the InSAR technique, such as quasi-real time observation, continuous recording and all-weather measurement. Accordingly, by means of the model this paper presents a method to decompose the InSAR slant range co-seismic displacement （i.e. LOS change） into horizontal and vertical displacement components. Approaching the real motion step by step, finally a serial of curves representing the co-seismic horizontal and vertical displacement component along the main earthquake fracture zone are approximately obtained.

Hazard Assessment of Co-seismic Landslides Based on Information Value Method:A Case in 2018 MW6.6 Hokkaido Earthquake,Japan

An MW6.6 earthquake occurred in eastern Hokkaido,Japan on September 6th,2018.Based on the pre-earthquake image from Google Earth and the post-earthquake image from high resolution(3 m)planet satellite,we manually interpret 9293 coseismic landslides and select 7 influencing factors of seismic landslide,such as elevation,slope,slope direction,road distance,flow distance,peak ground acceleration(PGA)and lithology.Then,9293 landslide points are randomly divided into training samples and validation samples with a proportion of 7:3.In detail,the training sample has 6505 landslide points and the validation sample has 2788 landslide points.The hazard risk assessment of seismic landslide is conducted by using the information value method and the study area is further divided into five risk grades,including very low risk area,low risk area,moderate risk area high risk area and very high risk area.The results show that there are 7576 landslides in high risk area and very high risk area,accounting for81.52%of the total landslide number,and the landslide area is 22.93 km^2,accounting for 74.35%of the total area.The hazard zoning is in high accordance with the actual situation.The evaluation results are tested by using the curve of cumulative percentage of hazardous area and cumulative percentage of landslides number.The results show that the success rate of the information value method is 78.50%and the prediction rate is 78.43%.The evaluation results are satisfactory,indicating that the hazard risk assessment results based on information value method may provide scientific reference for landslide hazard risk assessment as well as the disaster prevention and mitigation in the study area.

Co-seismic deformation derived from GPS observations during April 20th,2013 Lushan Earthquake,Sichuan,China

We process the standard 30 s, static GPS data and the 1 s, high-rate GPS （HRGPS） data provided by the Crustal Movement Observation Network of China with GAMIT/GLOBK software package, and obtain the co- seismic displacements of near field and far field, and the epoch-by-epoch time series of HRGPS during Lushan earthquake. GPS data from about 20 sites in Sichuan province, which located between 40 and 450 km from the epicenter, are analyzed so as to study the characteristics of the static displacements and the dynamic crustal defor- mations, with periods ranging from several minutes to over a month. The result shows that： the static displacements caused by Lushan earthquake are limited to several centi- meters; the nearest station SCTQ at 43 km from the epi- center has the largest static displacement of about 2 cm, while the other stations generally have insignificant dis- placements of less than 5 mm. the stations in the east ofSichuan-Yunnan region shifts 5-10 mm toward the southwest, and the stations in the middle-west of Sichuan Basin moves indistinctively 1-2 mm toward the northwest; station SCTQ has the largest kinematic displacement of about 4 and 3 cm peak-to-peak on the north and east component, respectively, and is much greater than the static permanent displacement; for the stations located at a distance greater than 150 km from the epicenter, the kinematic motions are generally insignificant; exception- ally, station SCNC and station SCSN in central Sichuan Basin have significant kinematic motions although they are more than 200 km away from the epicenter.

Analysis of the co-seismic responses of the fluid well pattern system in Jiangsu Province to the Wenchuan and Tohoku earthquakes

This thesis discusses the earthquake reflecting ability of the observation well pattern system of Jiangsu Province, China, which has been digitally renovated, and probes into the cause of the major differences in the earthquake reflecting abilities of well holes at different measurement points. This is achieved through the analysis of the co-seismic responses to the Wenchuan (2008; Ms8.0, China) and Tohoku (2011; Ms9.0, Japan) earthquakes. We found that the co-seismic response of water level from regional well holes in Jiangsu Province was stronger than that of water temperature. The water-level co-seismic response follows a consistent law and is closely related to the earthquake magnitude. The co-seismic response of water temperature strongly varied among well points, and was more often manifested as a slow restorative change. The co-seismic responses also varied based on tectonic elements. The response in central and northern Jiangsu was weaker than that of southern Jiangsu, possibly due to the thicker loess cover layer in central Jiangsu which makes it less effective at capturing the micro-changes of stress-strain states relative to the hilly land in the south. The more complicated geological structure in southern Jiangsu makes it contribute to greater changes in the state of underground water after a minor disturbance.

Modeling co-seismic thermal infrared brightness anomalies in petroliferous basins surrounding the North and East of the Qinghai–Tibet Plateau

Co-seismic gas leakage usually occurs on the edge of seismic faults in petroliferous basins,and it may have an impact on the local environment,such as the greenhouse effect,which can cause thermal infrared brightness anomalies.Using wavelet transform and power spectrum estimation methods,we processed brightness temperature data from the Chinese geostationary meteorological satellite FY-C/E.We report similarities between the co-seismic thermal infrared brightness(CTIB)anomalies before,during and after earthquakes that occurred at the edges of the Sichuan,Tarim,Qaidam,and Junggar basins surrounding the North and East of the Qinghai–Tibet Plateau in western China.Additionally,in each petroliferous basin,the area of a single CTIB anomaly accounted for 50%to 100%of the basin area,and the spatial distribution similarities in the CTIB anomalies existed before,during and after these earthquakes.To better interpret the similarities,we developed a basin warming effect model based on geological structures and topography.The model suggests that in a petroliferous basin with a subsurface gas reservoir,gas leakage could strengthen with the increasing stress before,during,and even after an earthquake.The accumulation of these gases,such as the greenhouse gases CH4 and CO2,results in the CTIB anomalies.In addition,we conclude that the CTIB anomalies are strengthened by the high mountains(altitude^5000 m)around the basins and the basins’independent climatic conditions.This work provides a new perspective from which to understand the CTIB anomalies in petroliferous basins surrounding the North and East of the Qinghai–Tibet Plateau.

Co-seismic responses of water wells in Jiangsu province to 2008 Wenchuan and 2011 Japan earthquakes

Co-seismic water-level and temperature changes of the 2008 magnitude - 8.0 Wenchuan and the 2011 magnitude-9.0 Japan earthquakes recorded at 10 observation wells in Jiangsu province are presented and analyzed. The data show that water level responded more regularly with earthquake magnitude and dis- tance than water temperature. The response was different for wells located in different tectonic units, being weaker in central and northern plain, which has a relatively thick surface layer of loess, than southern Jiangsu, which is hilly.

Co-seismic Coulomb stress changes on the northern Tanlu fault zone caused by the Tohoku-Oki M_(W)9.0 earthquake

Based on the spherical earth dislocation theory and a fault slip model of the Tohoku-Oki M_(W)9.0 earthquake,the co-seismic Coulomb failure stress changes(ΔCFS)on the northern Tanlu fault zone at depths of 0–40 km are calculated.By comparing two sets of results from the spherical earth dislocation theory and the semi-infinite space one,the effect of earth curvature on the calculation results is analyzed quantitatively.First,we systematically summarize previous researches related to the northern Tanlu fault zone,divide the fault zone as detailed as possible,give the geometric parameters of each segment,and establish a segmented structural model of the northern Tanlu fault zone.Second,we calculate the Coulomb stress changes on the northern Tanlu fault zone by using the spherical earth dislocation theory.The result shows the Coulomb stress changes are no more than 0.003 MPa,which proves the great earthquake did not significantly change the stress state of the fault zone.Finally,we quantitatively analyze the disparities between the results of semi-infinite space dislocation theory and the spherical earth one.The average disparity between them is about 7.7%on the northern Tanlu fault zone and is 16.8%on the Fangzheng graben,the maximum disparity on this graben reaches up to 25.5%.It indicates that the effect of earth curvature can not be ignored.So it’s necessary to use the spherical earth dislocation theory instead of the semi-infinite space one to study the Coulomb stress change in the far field.

Co-seismic changes of well water level and volume strain meter in capital area and its vicinity,due to the Nov.14,2001 Ms8.1 Kunlun Mountain earthquake,China

The Kunlunshan Mountain Ms8.1 earthquake, occurred in Nov.14, 2001, is the first event with magnitude more than 8 in the China earthquake monitoring history, specifically at the beginning of digital techniques in precursor monitoring networks. Any investigation of recorded data on this earthquake is very important for testing the operation of the digital monitoring networks and understanding the preparation, occurrence, and adjustment of stress/strain of strong continental earthquakes. In this paper we investigated the coseismic response changes of well water level of groundwater and volume strain meter of bore hole in digital earthquake monitoring network of Capital area and its vicinity, due to the Nov.14, 2001 Ms8.1 Kunlun Mountain earthquake. The responding time, shapes or manners, amplitudes, and lasting time of well water level and strain-meters to seismic wave are studied in comparison. Then we discussed the possibility that the response changes of groundwater to strong distant earthquakes can be understood as one kind of observing evidence of stress/strain changes induced by distant earthquake.

Co-seismic fault geometry and slip distribution of the 26 December 2004, giant Sumatra–Andaman earthquake constrained by GPS, coral reef, and remote sensing data

We analyze co-seismic displacement field of the 26 December 2004, giant Sumatra–Andaman earthquake derived from Global Position System observations,geological vertical measurement of coral head, and pivot line observed through remote sensing. Using the co-seismic displacement field and AK135 spherical layered Earth model, we invert co-seismic slip distribution along the seismic fault. We also search the best fault geometry model to fit the observed data. Assuming that the dip angle linearly increases in downward direction, the postfit residual variation of the inversed geometry model with dip angles linearly changing along fault strike are plotted. The geometry model with local minimum misfits is the one with dip angle linearly increasing along strike from 4.3oin top southernmost patch to 4.5oin top northernmost path and dip angle linearly increased. By using the fault shape and geodetic co-seismic data, we estimate the slip distribution on the curved fault. Our result shows that the earthquake ruptured ＊200-km width down to a depth of about 60 km.0.5–12.5 m of thrust slip is resolved with the largest slip centered around the central section of the rupture zone78N–108N in latitude. The estimated seismic moment is8.2 9 1022 N m, which is larger than estimation from the centroid moment magnitude（4.0 9 1022 N m）, and smaller than estimation from normal-mode oscillation data modeling（1.0 9 1023 N m）.

Determination of Smoothing Factor for the Inversion of Co-seismic Slip Distribution

For the determination of the smoothing factor (also known as the regularization parameter) in the co-seismic slip distribution inversion, the compromise curve between the model roughness and the data fitting residual is generally used to determine (in order to distinguish the method proposed in this paper, the method is called “L curve” according to its shape). Based on the L-curve, the Eclectic Intersection curve as a new method is proposed to determine the smoothing factor in this paper. The results of the simulated experiment show that the inversion accuracy of the parameters of the seismic slip distribution with the smoothing factor determined by the Eclectic Intersection curve method is better than that of the L curve method. Moreover, the Eclectic Intersection curve method and the L curve method are used to determine the smoothing factor of L’Aquila earthquake and the Taiwan Meinong earthquake slip distribution inversion respectively, and the inversion results are compared and analyzed. The analysis results show that the L’Aquila and the Taiwan Meinong actual earthquake slip distribution results are in the range of other scholars at home and abroad, and compared with the L curve method, the Eclectic Intersection curve method has advantages of high computation efficiency, no need to depend on data fitting degree and more appropriate of smoothing factor and so on.