Shuangxi Zhang

Recovering period of postseismic fluid pressure in fault valve

The present study aims to reveal the recovering period of the postseismic fluid pressure in fault zone, offering an insight into earthquake recurrence. Numerical modeling is performed based on a 2D simple layered fault-valve model to simulate the fluid activities within the earthquake fault. In order to demonstrate the features of postseismic fluid pressure in natural state, the interference of tectonic movements is not considered. The recovering period of postseismic fluid pressure includes a suddenchanging period and a much longer fluctuating period. Modeling results show that fault permeability and porosity are sensitive parameters and reversely proportional to the recovering period of the fluid pressure in earthquake fault zone. When the permeability reduces from 10-15 to 10-18 m2, the recovering period increases from 400 to 2 000 yrs, correspondently. The upper and lower fluid pressures are separated by the valve seal, causing their fluctuations in opposite tendencies.

Fault slip model of 2013 Lushan Earthquake retrieved based on GPS coseismic displacements

Lushan Earthquake (~Mw 6.6) occurred in Sichuan Province of China on 20 April 2013, was the largest earthquake in Longmenshan fault belt since 2008 Wenchuan Earthquake. To better understand its rupture pattern, we focused on the influences of fault parameters on fault slips and performed fault slip inversion using Akaike’s Bayesian Information Criterion (ABIC) method. Based on GPS coseismic data, our inverted results showed that the fault slip was mainly confined at depths. The maximum slip amplitude is about 0.7 m, and the scalar seismic moment is about 9.47×1018 N·m. Slip pattern reveals that the earthquake occurred on the thrust fault with large dip-slip and small strike-slip, such a simple fault slip represents no second sub-event occurred. The Coulomb stress changes (ΔCFF) matched the most aftershocks with negative anomalies. The inverted results demonstrated that the source parameters have significant impacts on fault slip distribution, especially on the slip direction and maximum displacement.

Distribution of Intra-Crustal Low Velocity Zones beneath Yunnan from Seismic Ambient Noise Tomography

Previous studies have reached consensus that low velocity zones are widespread in the crust beneath Yunnan region. However, the relationships between the low velocity zones and large faults, earthquake distribution are less investigated by available studies. By analyzing the seismic ambient noise recorded by Yunnan Seismic Networks and Tengchong volcano array, we construct a 3D crustal shear wave velocity model for the Yunnan region, which provides more details of the distribution of intra-crustal low velocity zones in all of Yunnan. The distribution of low velocity zones shows different features at different depths. At shallow depths, the results are well correlated with near surface geological features. With increasing depth, the low velocity zones are gradually concentrated to the northern part of our study area, most likely reflecting variations in crustal thickness beneath the Yunnan region. The low velocity zones are truncated at depth by several large faults in Yunnan. It is interesting that most strong earthquakes (Ms≥5.0) occurred in Yunnan are distributed in low velocity zones or the transition zones between low and high velocity anomalies within the upper-to-middle crust. The crustal structure is composed of a brittle, seismically active upper-to-middle crust and a warm, aseismic lower crust.

Crustal Structure of Yunnan Province of China from Teleseismic Receiver Functions: Implications for Regional Crust Evolution

Yunnan Province is located on the southeastern margin of Tibet and represents an important marker in understanding the tectonic evolution of Tibetan Plateau. In this study, we calculated teleseismic P-wave receiver functions at 49 permanent broadband seismic stations in Yunnan Province and estimated crustal thickness and the bulk crust ratios of P-wave to S-wave velocities using the H-κ method together with more detailed crustal structural profiles from the common conversion point stacking method. There is a significant transition of Moho interface and lower crustal composition along latitude 26°N in northwestern Yunnan. Decrease of crustal thickness with a concomitant increase of Poisson’s ratio occurs at station CUX. An interesting phenomenon is that a step-like Moho fashion is observed at several stations, which might correspond to local thermal activities, such as partial melt/lower crust delamination. Our results show changes in crustal properties appear to be associated with varieties in upper mantle structure and compositions, combined with other previous studies. We propose the controlling factor of the dynamic processes below 26°N is the result of eastern forward subduction of the Indian Plate; the northern part is controlled by the redirected material flow from the SE Tibet.

Short-Impending Earthquake Anomaly Index Extraction of GNSS Continuous Observation Data in Yunnan, Southwestern China

This paper presents a comprehensive area expansion prediction index method to apply GNSS for short-impending prediction of earthquakes. Based on continuous GNSS observation data from Yunnan Province, a displacement field was detected after data cycle-slip repair using precision data processing software and geophysical field effect model correction. The Yunnan area was divided into 56 grid cells for displacement field interpolation to obtain a more uniform displacement field and a strain field variation time series. The pre-earthquake response of each grid-cell expansion time series was evaluated and synthesized to extract a short-impending earthquake anomaly identification index. The results show that this index indicated occurrence times and hypocenter for earthquakes of magnitude M≥5. Fourteen earthquakes were predicted accurately, and there were five false reports. This index can therefore be used for the short-impending prediction of earthquakes.

Finite-difference modeling of surface waves in poroelastic media and stress mirror conditions

During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surfacewave phenomenon. To accurately determine the effects of this interface on surface-wave propagation, the boundary conditions must be accurately modeled. In this paper, we present a numerical approach based on the dynamic poroelasticity for a space–time-domain staggeredgrid finite-difference simulation in porous media that contain a free-surface boundary. We propose a generalized stess mirror formulation of the free-surface boundary for solids and fluids in porous media for the grid mesh on which lays the free-surface plane. Its analog is that used for elastic media, which is suitable for precise and stable Rayleigh-type surface-wave modeling. The results of our analysis of first kind of Rayleigh (R1) waves obtained by this model demonstrate that the discretization of the mesh in a similar way to that for elastic media can realize stable numerical solutions with acceptable precision. We present numerical examples demonstrating the efficiency and accuracy of our proposed method.