Source time functions of the Gonghe，China earthquake retrieved from long－period digital waveform data using empirical Green＇s function technique

An earthquake of Ms= 6, 9 occurred at the Gonghe, Qinghai Province, China on April 26, 1990. Three larger aftershocks took place at the same region, Ms= 5. 0 on May 7, 1990, Ms= 6. 0 on Jan. 3, 1994 and Ms= 5. 7on Feb. 16, 1994. The long-period recordings of the main shock from China Digital Seismograph Network (CDSN) are deconvolved for the source time functions by the correspondent0 recordings of the three aftershocks asempirical Green's functions (EGFs). No matter which aftershock is taken as EGF, the relative source time functions (RSTFs) Obtained are nearly identical. The RSTFs suggest the Ms= 6. 9 event consists of at least two subevents with approximately equal size whose occurrence times are about 30 s apart, the first one has a duration of 12 s and a rise time of about 5 s, and the second one has a duration of 17 s and a rise time of about & s. COmParing the RSTFs obtained from P- and SH-phases respectively, we notice that those from SH-phases are a slightly more complex than those from p-phases, implying other finer subevents exist during the process of the main shock. It is interesting that the results from the EGF deconvolution of long-Period way form data are in good agreement with the results from the moment tensor inversion and from the EGF deconvolution of broadband waveform data. Additionally, the two larger aftershocks are deconvolved for their RSTFs. The deconvolution results show that the processes of the Ms= 6. 0 event on Jan. 3, 1994 and the Ms= 5. 7 event on Feb. 16,1994 are quite simple, both RSTFs are single impulses.The RSTFs of the Ms= 6. 9 main shock obtained from different stations are noticed to be azimuthally dependent, whose shapes are a slightly different with different stations. However, the RSTFs of the two smaller aftershocks are not azimuthally dependent. The integrations of RSTFs over the processes are quite close to each other, i. e., the scalar seismic moments estimated from different stations are in good agreement. Finally the scalar seismic moments of the three aftershocks are compared. The relative scalar seismic moment Of the three aftershocks deduced from the relative scalar seismic moments of the Ms=6. 9 main shock are very close to those inverted directly from the EGF deconvolution. The relative scalar seismic moment of the Ms =6. 9 main shock calculated using the three aftershocks as EGF are 22 (the Ms= 6. 0 aftershock being EGF), 26 (the Ms= 5. 7 aftershock being EGF) and 66 (the Ms= 5. 5 aftershock being EGF), respectively. Deducingfrom those results, the relative scalar sesimic moments of the Ms= 6. 0 to the Ms= 5. 7 events, the Ms= 6. 0 tothe Ms= 5. 5 events and the Ms= 5. 7 to the Ms= 5. 5 events are 1. 18, 3. 00 and 2. 54, respectively. The correspondent relative scalar seismic moments calculated directly from the waveform recordings are 1. 15, 3. 43, and 3. 05.

Source parameters of the Gonghe,Qinghai Province,China,earthquake from inversion of digital broadband waveform data

An earthquake of M S=6.9 occurred in Gonghe County, Qinghai Province, China on April 26, 1990.This earthquake was followed by three larger aftershocks of M S=5.5 on May 7, 1990, M S=6.0 on Jan.3, 1994, and M S=5.7 on Feb.16, 1994, consecutively. The moment tensors of these earthquakes as function of time were obtained by the technique of moment tensor inversion in frequency domain . The results inverted indicate that these earthquakes had a very similar focal mechanism of predominantly reverse faulting on a plane striking NWW, dipping to SSW.The scalar seismic moments of these earthquakes are M 0=9.4×10 18 Nm for the M S=6.9 event, 8.0×10 16 Nm for the M S=5.5 event, 4.9×10 17 Nm for the M S =6.0 event and 2.9×10 17 Nm for the M S=5.7 event, respectively. The results inverted also show that the source processes of these events were significantly different. The main shock had a very complex process, consisting of two distinct sub events with comparable sizes. The first sub event occurred in the first 12s, having a seismic moment of 4.7×10 18 Nm, and the second one continued from 31s to 41s, having a seismic moment of 2.5×10 18 Nm. In addition, a much smaller sub event, having a seismic moment of about 2.1×10 18 Nm, may exist in the interval of 12 s and 31 s, In contrast, the source processes of the three aftershocks are quite simple. The source time function of each of aftershocks is a single impulse, suggestting that each of aftershocks consists of a mainly uninterrupted rupture. The rise times and total rupture durations are 4 s and 11 s for the M S=5.5 event, 6 s and 16 s for the M S= 6.0 event and 6 s and 13 s for the M S=5.7 event, respectively.

Effective relaxation time and viscosity of the earth inferred from the postseismic vertical deformations of the 1990 MS=7.0 Gonghe earthquake

The postseismic vertical deformation rates of the 1990 Gonghe M S=7.0 earthquake appears to have decreased exponentially. Based on Okada′s coseismic surface displacement solution caused by a uniform fault slip in an elastic homogeneous half space, we derived its postseismic surface displacement by using a single layer standard linear solid model, and further derived a simplified formula for determining the effective relaxation time and viscosity of the earth, which is independent of the dislocation parameters of the causative fault. From the postseismic vertical deformation of the 1990 Gonghe earthquake, we inferred that the effective relaxation time defined by τ = η/μ is 2.6 years, and the effective viscosity η is about 10 18 Pa·s.

Near-field surface deformation during the April 20,2013,Ms7.0 Lushan earthquake measured by 1-Hz GNSS

The April 20, 2013, Ms7.0 Lushan earthquake was successfully recorded by closely spaced Continuous Global Positioning System （CGPS） stations owned by the Crustal Movement Observation Network of Chi- na （CMONC）. The 1-Hz GNSS data from eight CGPS stations, which are located between 30 km and 200 km from the hypocenter, were processed within quasi-real-time. The near-field surface deformation indicated the following characteristics ： the near-field movements were limited to several centimeters ; the peak of the deformation wave was significantly larger than the static permanent offset; at the beginning of the event, the north wall of the fault moved to the southeast as the south wall moved to the southwest ; station SCTQ, which was the closest station to the hypocenter at 30 km, had the largest static permanent displacement of 2 cm; the peaks of the deformation waves were 1.5 cm, 5 cm and 3 cm, to the east, the south and vertically upward, respectively ; and the peaks of velocity and acceleration, derived from the deformation, were 3.4 cm/s and 5.3 cm/s^2,respectively.

Cround vertical deformation in earthquake of Gonghe,Qinghai Province

ＣｒｏｕｎｄｖｅｒｔｉｃａｌｄｅｆｏｒｍａｔｉｏｎｉｎｅａｒｔｈｑｕａｋｅｏｆＧｏｎｇｈｅ，ＱｉｎｇｈａｉＰｒｏｖｉｎｃｅＳｈｏｕ－ＷｅｎＧＯＮＧ（巩守文）ａｎｄＦｅｎｇ－ＹｉｎｇＧＵＯ（郭逢英）（ＳｅｃｏｎｄＣｒｕｓｔａｌＤｅｆｏｒｍａｔｉｏｎＭｏ...

Numerical simulation of precursors of Gonghe M7.0 earthquake in Qinghai Province

Based on a simplified media model of Gonghe area, the precursor characteristics of Gonghe M 7.0 earthquake in 1990 are simulated in this paper by using the constitutive relationship of binary medium (solid and water). The results show that the simulated state distribution and extension variation of media are identical with the spatial and temporal distribution of reliable anomalies before Gonghe earthquake. The study also suggests that the development of the Gonghe earthquake has experienced a series of processes such as elastic deformation of large scale, early nonelastic dilatation, strain softening and elastic recovery in the neighbour region of Xining, nonelasscale, early nonelastic dilatation, strain softening and elastic recovery in the neighbour region of nonelastic dilatation of high dense block near seismic source and earthquake occurrence. Therefore, it can be concluded that this earthquake development is neither a simple process for fissures developing and linking up with one another, nor a process of nonelastic volume expanion and water flowing into medium around seismic source, conrary, it is a complicated medium saate changing process, and the style of such changing depends on geological structure environment of seismic source and its neighbour region. It is considered that different earthquakes are accompanied by different geological conditions, the spatial and temporal behavior of their precusors are certainly different.

Seismic response of cracking features in Jubao Mountain during the aftershocks of Jiuzhaigou Ms7.0 earthquake

Jiuzhaigou is a world-heritage site located in the plateau area of Northwest Sichuan Province,China.Serious slope failures in the epicentral area were triggered by the Ms7.0 Jiuzhaigou earthquake occurred on August 8,2017.The source areas of the hazards are usually concentrated near ridge crests,revealingthe possible occurrence of ground motion amplification phenomena.To explore the role of the amplification of ground motions in the formation of earthquake-triggered slope failures,two seismometers were installed,on the next day after the main shock,at the bottom of the slopeof Jubao Mountain near the seismogenic fault.The two monitoring sites are located at elevations of 1414 m(J1)and 1551 m(J2,the top of the mountain).Five aftershocks were recorded by the monitoring instruments.We compared the mean levels of the peak ground acceleration(PGA)observed at different locations,and investigated the directional variations inthe shaking energy by analyzing the polar diagrams of the Arias intensity(Ia).Then,in order to identify the directional resonance phenomenonandtheir frequencies and amplification coefficients,we examined the horizontal-to-vertical spectral ratio(HVSR)and the standard spectral ratio(SSR).Polar diagrams of theArias intensity(Ia)indicated that the site response of Jubao Mountain showed a pronounced directivity(in theEW direction)with shaking maxima near the hill top oriented orthogonally to the elongation of the relief.We observed anobvious resonance phenomenonat site J2 at relatively low frequencies(2.5-9 Hz)and very weak spectral amplifications at site J1 at high frequencies(5-15 Hz),which suggested that the predominant frequency of monitoring site J2 was obviously attenuated and that the difference in the spectra was related to the influences of the local-scale site conditions of the whole mountain.The results of spectral ratio analyses(HVSR and SSR)showed that the direction of resonance was concentrated around an EW orientation,and the amplification factors near the hill top were larger than 2.It suggests that geologic factors also play a significant role in the anisotropic amplifications affecting the tops of slopes besides the topographic effects.

Emergency Response to the MS7.0 Jiuzhaigou,Sichuan Earthquake,and Characteristics of Seismic Disasters in the Stricken Area in Gansu Province

This paper introduces the response process of the Gansu Earthquake Agency during the Jiuzhaigou M_S7.0 earthquake in Sichuan Province,including earthquake emergency disposal procedures,information reports,disaster investigation and intensity assessment,seismic monitoring and trend determination,and emergency dissemination. This paper reveals the characteristics of earthquake damage in the quake-hit areas of Gansu Province,draws some corresponding conclusions and summarizes the countermeasures for recovery and reconstruction in the quake-hit areas of Gansu Province.

Characteristics and Analysis of the Seismic Damage from the MS7.0 Jiuzhaigou Earthquake in the Area between Songpan and Jiuzhaigou

This paper expounds the features of the buildings and analyzes the seismic disaster characteristics of the Jiuzhaigou M_S7. 0 earthquake in the area between Songpan and Jiuzhaigou. New buildings (especially the frame structure) had good anti-seismic performance,but brick-wood structures and brick-concrete structures sustained large amounts of damage in the earthquake. By computing the seismic damage index,we found that the seismic damage index of the frame structure was far less than that of civil structures and brick-wood structures. The seismic damage index of frame structures were all zero in the Ⅵ area,and increased rapidly with the increase of intensity,but the increasing range was reduced. We also discussed how to evaluate the intensity in areas where there was a lack of buildings or there was only one structure type,which can be referenced in future field work.

A rapid stability assessment of China's IGS sites after the Ms7.0 Lushan earthquake

A rapid and accurate assessment of the stability of surveying and mapping reference points is important for post - disaster rescue, disaster relief and reconstruction activities. Using Precise Point Positioning （PPP） technology, a rapid assessment of the stability of the IGS sites in China was performed after the Ms7.0 Lushan earthquake using rapid precise ephemeris and rapid precise satellite clock products. The results show that the earthquake had a very small impact and did not cause significant permanent deformation at the IGS sites. Most of the sites were unaffected and remained stable after the earthquake.

Dislocation parameters of Gonghe earthquake jointly inferred by using genetic algorithms and least squares method

The Second Crustal Deformation Monitoring Center, China Seismological Bureau, has detected a marked uplift associated with the Gonghe Ms=7.0 earthquake on April 26, 1990, Qinghai Province. From the observed vertical deformations and using a rectangular uniform slip model in a homogeneous elastic half space, we first employ genetic algorithms (GA) to infer the approximate global optimal solution, and further use least squares method to get more accurate global optimal solution by taking the approximate solution of GA as the initial parameters of least squares. The inversion results show that the causative fault of Gonghe Ms=7.0 earthquake is a right-lateral reverse fault with strike NW60°, dip SW and dip angle 37°, the coseismic fracture length, width and slip are 37 km, 6 km and 2.7 m respectively. Combination of GA and least squares algorithms is an effective joint inversion method, which could not only escape from local optimum of least squares, but also solve the slow convergence problem of GA after reaching adjacency of global optimal solution.

Comparison of Computational Methods for Instrumental Intensities Using the MS7.0 Jiuzhaigou Earthquake Data

After the occurrence of an earthquake,strong motion observation networks can record ground motion at distributed observation stations. Based on the ground motion parameters from these records,the spatial distribution of seismic intensity can be quickly determined,and the degree of damage in different areas can be estimated. This information provides the technical basis for the emergency response,so as to ensure that rescue teams can reach extreme earthquake areas and carry out the search and rescue operation in an accurate and timely manner to reduce casualties and property loss. In this paper,we introduced 7 intensity algorithms and compared the results with the records of the Jiuzhaigou M_S7. 0 earthquake. We found that the differences between the instrumental intensities calculated by each method and the macro intensities were within a 1-degree range,which suggested good practicality of these different methods. The results calculated by the industrial standards-based calculation method and the integrated test showed good consistency.

Research on the Dynamic Evolution of Deformation Field and the Acquisition of Seismogenic Information Before Lijiang Earthquake With Ms7.0

By using the integral dynamic adjustment of the "segmentation rate," this paper processes the leveling repetition measurement data of the Dianxi(western Yunnan)area since the 1980s,makes the fitting of the velocity surface for vertical deformation with multi-surface function,and splits information of the vertical deformation field.The dynamic evolution of the deformation field and vertical deformation information norm-shear deformation,extension deformation,volumetric deformation,and quasi acceleration,which are the distribution characteristics of variation,are studied by using the information split before the Lijiang earthquake with MS7.0(February 3,1996)in the north of the Dianxi(western Yunnan)leveling network.The results indicate that in Dianxi(western Yunnan)area,the acute vertical differential movement of the large-scale deformation velocity high gradient belt and high information norm distribution in the Lijiang,Yongsheng,and its adjacent areas in the middle stage of the 1980s is the reflection

Geometry and tectonic deformation of the seismogenic structure for the 8 August 2017 M_S 7.0 Jiuzhaigou earthquake sequence,northern Sichuan, China

To reveal the geometry of the seismogenic structure of the Aug. 8, 2017 M_S 7.0 Jiuzhaigou earthquake in northern Sichuan,data from the regional seismic network from the time of the main event to Oct. 31, 2017 were used to relocate the earthquake sequence by the tomoDD program, and the focal mechanism solutions and centroid depths of the M_L ≥ 3.5 events in the sequence were determined using the CAP waveform inversion method. Further, the segmental tectonic deformation characteristics of the seismogenic faults were analyzed preliminarily by using strain rosettes and areal strains(As). The results indicate:(1) The relocated M_S 7.0 Jiuzhaigou earthquake sequence displays a narrow ～ 38 km long NNW-SSE-trending zone between the NW-striking Tazang Fault and the nearly NSstriking Minjiang Fault, two branches of the East Kunlun Fault Zone. The spatial distribution of the sequence is narrow and deep for the southern segment, and relatively wide and shallow for the northern segment. The initial rupture depth of the mainshock is 12.5 km, the dominant depth range of the aftershock sequence is between 0 and 10 km with an average depth of 6.7 km. The mainshock epicenter is located in the middle of the aftershock region, showing a bilateral rupture behavior. The centroid depths of 32 M_L ≥ 3.5 events range from 3 to 12 km with a mean of about 7.3 km, consistent with the predominant focal depth of the whole sequence.(2) The geometric structure of the seismogenic fault on the southern section of the aftershock area(south of the mainshock) is relatively simple, with overall strike of ～150° and dip angle ～75°, but the dip angle and dip-orientation exhibit some variation along the segment. The seismogenic structure on the northern segment is more complicated; several faults, including the Minjiang Fault, may be responsible for the aftershock activities. The overall strike of this section is ～159° and dip angle is ～59°, illustrating a certain clockwise rotation and a smaller dip angle than the southern segment. The differences between the two segments demonstrate variation of the geometric structure along the seismogenic faults.(3) The focal mechanism solutions of 32 M_L ≥ 3.5 events in the earthquake sequence have obvious segmental characteristics. Strike-slip earthquakes are dominant on the southern segment, while 50% of events on the northern segment are thrusting and oblique thrusting earthquakes, revealing significant differences in the kinematic features of the seismogenic faults between the two segments.(4) The strain rosettes for the mainshock and the entire sequence of 31 M_L ≥ 3.5 aftershocks correspond to strike-slip type with NWW-SEE compressional white lobes and NNE-SSW extensional black lobes of nearly similar size. The strain rosette and As value of the entire sequence of 22 M_L ≥ 3.5 events on the southern segment are the same as those of the M_S 7.0 mainshock,indicating that the tectonic deformation here is strike-slip. However, the strain rosette of the entire sequence of 10 M_L ≥ 3.5 events on the northern segment show prominent white compressional lobes and small black extensional lobes, and the related As value is up to 0.52,indicating that the tectonic deformation of this segment is oblique thrusting with a certain strike-slip component. Differences between the two segments all reveal distinctly obvious segmental characteristics of the tectonic deformation of the seismogenic faults for the Jiuzhaigou earthquake sequence.

Upper crustal velocity and seismogenic environment of the M7.0 Jiuzhaigou earthquake region in Sichuan, China

On August 8,2017,a magnitude 7.0 earthquake occurred in Jiuzhaigou County,Sichuan Province,China.The deep seismogenic environment and potential seismic risk in the eastern margin of Tibetan Plateau have once again attracted the close attention of seismologists and scholars at home and abroad.The post-earthquake scientific investigation could not identify noticeable surface rupture zones in the affected area;the complex tectonic background and the reason(s)for the frequent seismicity in the Jiuzhaigou earthquake region are unclear.In order to reveal the characteristics of the deep medium and the seismogenic environment of the M7.0 Jiuzhaigou earthquake region,and to interpret the tectonic background and genesis of the seismicity comprehensively,in this paper,we have reviewed all available observation data recorded by the regional digital seismic networks and large-scale,dense mobile seismic array(China Array)for the northern section of the North-South Seismic Belt around Jiuzhaigou earthquake region.Using double-difference seismic tomography method to invert the three-dimensional P-wave velocity structure characteristics of the upper crust around the Jiuzhaigou earthquake region,we have analyzed and discussed such scientific questions as the relationship between the velocity structure characteristics and seismicity in the Jiuzhaigou earthquake region,its deep tectonic environment,and the ongoing seismic risk in this region.We report that:the P-wave velocity structure of the upper crust around the Jiuzhaigoug earthquake region exhibits obvious lateral inhomogeneity;the distribution characteristics of the shallow P-wave velocity structure are closely related to surface geological structure and formation lithology;the M7.0 Jiuzhaigou earthquake sequence is closely related to the velocity structure of the upper crust;the mainshock of the M7.0 earthquake occurred in the upper crust;the inhomogeneous variation of the velocity structure of the Jiuzhaigou earthquake area and its surrounding medium appears to be the deep structural factor controlling the spatial distribution of the mainshock and its sequence.The 3D P-wave velocity structure also suggests that the crustal low-velocity layer of northeastern SGB(Songpan-GarzêBlock)stretches into MSM(Minshan Mountain),and migrates to the northeast,but the tendency to emerge as a shallow layer is impeded by the high-velocity zone of Nanping Nappe tectonics and the Bikou Block.Our results reveal an uneven distribution of high-and low-velocity structures around the Tazang segment of the East Kunlun fault zone.Given that the rupture caused by the Jiuzhaigou earthquake has enhanced the stress fields at both ends of the seismogenic fault,it is very important to stay vigilant to possible seismic hazards in the large seismic gap at the Maqu-Maqên segment of the East Kunlun fault zone.

Which velocity model is more suitable for the 2017 M_S7.0 Jiuzhaigou earthquake?

On Aug.8, 2017, an M_S 7.0 earthquake struck Jiuzhaigou, a county of Sichuan province, China. A number of investigations and studies have been conducted, some of which involved local velocity models. However, the suitability of these models has not been properly addressed. Here we collect 11 already-existing models, including those used in studies of the 2017 M_S 7.0 Jiuzhaigou earthquake,choose 10 local stations surrounding the earthquake, and employ the same technique(TRIT) to relocate the hypocenter. And furthermore, we choose a more suitable model from the 11 already-existed models by analyzing the relocation process and the relocated results for reasonability. Finally, our conclusion is that the model Fang 2018 is more suitable and the hypocenter parameters, 103.801°E,33.192°N and 15.8 km for longitude, latitude and depth, respectively, and 2017-08-08 13:19:46.66 for its origin time, based on this model should be recommended for the 2017 M_S7.0 Jiuzhaigou earthquake.

Directional seismic response to the complex topography:A case study of 2013 Lushan Ms 7.0 earthquake

Azimuthal variations in site response can provide a good insight into the site amplification and seismic conditions of geohazard occurrences.In this study,multiple directional site response methods,including D-Arias(Directional-Arias),D-SER(Directional-Shaking energy ratio),D-HVSR(Directional-Horizontal to vertical spectral ratio)and D-SSR(Directional–Standard spectral ratio),are adopted to analyse seismic data of the 2013 Lushan Ms 7.0 earthquake captured by the self-established Lengzhuguan(LZG)station which consists of the complex topography of isolated ridge,large mountain and some typical micro-reliefs.The results show that the isolated ridge could cause stronger site responses than the large mountain,and whose pronounced response direction is roughly perpendicular to its ridgeline.With the growth of elevation,the siteresonant frequency decreases.The different microreliefs on the mountain cause different site responses,which present as protruding slope>linear slope.The site response mainly exists on the surficial layer of the mountain and shows that with the increase of the distance to mountain surface,the site response gets weaker,the site resonant frequency gets higher,and the pronounced response direction is perpendicular to its ridgeline.

Coseismic displacements and inospheric changes of the 2013 Ms7. 0 Lushan earthquake from GPS measurements

By inverting GPS data recorded at stations of the Crustal Movement Observation Network of China （CMONOC） near the 2013 Lushan Ms7.0 earthquake, we found a horizontal displacement of 22 mm at a site about 32 kin SW of the epicenter and vertical displacements of as much as 12.4 mm at several sites. The vertical displacements were generally uplift on the west side of the nearby Longmenshan fault zone and subsidence on the east side. We also found coseismic ionospheric disturbances about 0.5 to 0.9 TECU in amplitude that lasted for about one hour.

Anomaly Feature Extraction of the Gravity Field before the Jiuzhaigou Earthquake in 2017

The M_S7. 0 Jiuzhaigou earthquake occurred on August 8,2017. The earthquake occurred in the vicinity of the Tazang fault,the Minjiang fault and the Huya fault,where the focal mechanism is of the strike slip type. The static and dynamic anomalies of the gravity field can provide important physical field information for studying the structural properties of deep crust. Multi-scale decomposition techniques are used to separate Bouguer gravity at different depths and give some explanation to gravity variations at different time space scales. The results indicate that the wavelet multi-scale results of the EGM2008 model and the measured gravity data are consistent. Through comparative analysis,it is found that the Jiuzhaigou earthquake occurred in the stress enhanced region. The variation of gravity field at different time scales has a certain scientific significance for further understanding potential earthquake risk trend.

Characteristics of the Variation of Cross Fault Short-Baseline and Short-Range Leveling in Western Yunnan Earthquake Test Site Before Lijiang M_S7.0 Earthquake

Results of analysis of variation of cross fault short-baseline and short-range leveling in Western Yunnan Earthquake Test Site (WYETS), results show that among five observation stations of cross fault short-baseline and short-range leveling in WYETS before the Lijiang MS7.0 Earthquake occurred in February 1996 only Yongsheng observation station (epicentral distance 82 km) located at Chenghai fault shows great variation about one year before the earthquake. And the nearest observation station, Lijiang (epicentral distance 42 km); presents great coseismic variation, but does not show obvious anomalous variation before the earthquake. There are no significant variations related to the earthquake at the other three observation stations. Two methods are used in analysis of the observed data and some valuable results have been obtained.