2019 CGU Report Discussion on the Inverted Tectonic Characteristics and Seismic Control Mechanisms of the Maodong Fault Zone in Eastern China

Topic 27: Regional Seismic Tectonics and Their Seismic Hazard and Seismological Geological Disasters
Location:Beijing, Beijing International Conference Center, Conference Room 10 (3rd Floor, Room 301)
Time:October 28, 2019 14:55-15:10
Title:Inverted Tectonic Characteristics and Seismic Control Mechanisms of the Maodong Fault Zone in Eastern China: Insights from Discrete Element Numerical Simulation
Speaker:Wu Zhenyun, Yin Hongwei, Li Changsheng Donghua University of Science and Technology, Nanjing University

Abstract:By analyzing the geological structure and historical seismic characteristics of the Maodong Fault Zone, combined with discrete element simulation experiments, the tectonic inversion and fault activity characteristics of the faulted basin under compressive stress were analyzed. The study explored the inverted tectonic characteristics of the Maodong Fault Zone, as a faulted basin, under the current regional stress field and the mechanism of normal fault inversion for seismic control. The results indicate that:(1)Under compressive stress, tectonic inversion occurs in the faulted basin, which reactivates pre-existing normal faults and leads to the deformation of shallow strata, but the pre-existing faults in the basin generally maintain a normal fault shape;(2)During the process of tectonic inversion, stress-strain convergence is more likely to occur preferentially along the basal fault plane, gradually propagating into the sedimentary strata, ultimately triggering structural deformation at the shallow surface;(3)The Maodong Fault Zone is an active fault group, with the Maodong Fault and other boundary faults (basal faults) controlling basin activity. Along such faults, stress-strain convergence is evident, making them capable of generating medium to large earthquakes. Additionally, other concealed secondary pre-existing faults within the faulted basin may also have the potential for activity, with stress-strain convergence occurring in the strata between fault blocks under compressive-torsional stress fields, making them capable of generating shallow-source small earthquakes. Keywords: Faulted Basin; Discrete Element Simulation; Maodong Fault Zone; Tectonic Inversion; Fault Activity

The eastern part of China is influenced by the combined effects of the subduction of the Pacific Plate to the west and the subduction of the Philippine Plate to the northwest towards the Eurasian continent. The main characteristics of the modern tectonic stress field are expressed as compressive stress in the northeast-east to southwest-west direction (Figure 1a). Affected by the dextral strike-slip movement along the Tan-Lu Fault Zone and the dextral strike-slip movement along the southeastern coastal fault zones, the eastern part of the Lower Yangtze Block is mainly located in the compressive-torsional transition zone between these forces (Figure 1a). The modern tectonic stress field in the Lower Yangtze Block shows distinct compressive-torsional characteristics, leading to the development of both normal and reverse fault structures within the faulted basins of the region (Figure 1b).

Figure 1: Regional Location Map of the Study Area and Distribution of Historical Earthquakes
Figure 1: Regional Location Map of the Study Area and Distribution of Historical Earthquakes

Maodong Faulted Basin

Within the faulted basin (e.g., the Zhiqi Bridge Depression), high-angle faults are distributed, and these faults generally exhibit normal fault characteristics (Figure 2a). On the western side of the depression, the faults dip to the southeast (SEE); on the eastern side of the depression, the faults dip to the northwest (NWW) (Figure 2a). Previous studies along the Maodong Fault have conducted detailed geological surveys, resulting in representative surface geological profiles of the Maodong Faulted Basin (Figures 2b to 2d). These geological profiles indicate that along the Maodong Fault, the strata have undergone varying degrees of folding deformation. In addition, along the Maodong Fault, some Quaternary sedimentary strata in the faulted basin have been disrupted by faulting (Figure 2b) or have experienced co-tectonic deformation (Figure 2d), further demonstrating that the reverse faulting activity along the Maodong Fault has persisted into the Quaternary.

Figure 2: Exploration Profile of the Maodong Faulted Basin
Figure 2: Exploration Profile of the Maodong Faulted Basin

Model Design

The initial model for the discrete element simulation experiment in this study is shown in Figure 5. The model includes two types of basic particles representing rock grains, with radii set to 40 m and 60 m, respectively. All particles are randomly distributed within a closed square box, and they are allowed to settle and compact under the influence of gravity. In the experimental design, pre-existing normal faults (weakened faults) within the basin are defined, and the hanging wall and footwall of the faults are allowed to deform freely, so as to observe the deformation patterns of the hanging and footwalls during the compressive inversion process. To facilitate observation, the particles are marked with different colors: red represents the pre-existing normal fault, gray represents the basement strata, and blue, purple, and yellow represent the sedimentary strata. The left and bottom boundaries of the model are fixed, while the right boundary is pushed leftward to generate a compressive effect. In this experiment, the experimental model parameters were set based on previous designs. The particle density of the basement rock layer is 2500 kg·m^-3, and the particle density of the cover rock layer is 2200 kg·m^-3. At the same time, the friction coefficient for each layer is set to 0.3, and the particles exhibit cohesive forces. The friction coefficient between the pre-existing weakened faults is set to 0, and there is no cohesive force between the faults and any particles. The Poisson’s ratio for all particles in the model is set to 0.2.

Figure 3: Discrete Element Numerical Simulation Experimental Model Design for Tectonic Inversion in Faulted Basins
Figure 3: Discrete Element Numerical Simulation Experimental Model Design for Tectonic Inversion in Faulted Basins

In the initial stage of the experiment, the model did not exhibit significant tectonic deformation (Figure 4a1). However, strain analysis of the model revealed the development of volumetric and shear strains, primarily distributed along the pre-existing basement fault surface. This indicates that, in this stage, compressive stress mainly concentrated along the pre-existing basement fault, while the stress concentration in the sedimentary layers (cover layers) was minimal, resulting in almost no strain distribution (Figures 4a2 to 4a3). As compressive stress continued to act (e.g., when the compressive displacement reached 4%), the internal structural deformation and strain distribution remained concentrated along the basement fault, with only slight volumetric strain starting to develop in the shallow layers of the sedimentary cover (Figure 4b). When the compressive displacement reached 6.5%, the pre-existing basement fault began to reverse and activate, resulting in displacement (the fault maintained a normal fault character overall). Strain began to expand gradually from the basement fault surface into the sedimentary cover layers, and slight surface topography changes occurred within the model (Figure 4c), indicating that the stress had propagated to the surface and caused minor surface deformation.

Experimental Results

As the compressive displacement continued to increase (8.5%–13%), the degree of subsidence within the model began to increase, while the reverse displacement of the pre-existing basement fault showed little change, and the overall structural morphology also remained nearly unchanged (Figures 4d to 4e). On the side near the compressive end, a triangular strain region formed in the upper sedimentary layers (cover layers) above the concealed pre-existing basement fault (Figures 4d2 to 4d3). In this region, both volumetric strain and shear strain gradually increased with the increase in compressive displacement, indicating that compressive stress along the pre-existing basement fault gradually concentrated in this area. Additionally, along the large fault at the left basement boundary, strain also gradually extended into the shallow layers, and the fault reverse displacement started to gradually increase. When the compressive displacement increased to 15.5%, the large fault along the left basement boundary exhibited significant displacement, and reverse fault-related folds began to develop (Figures 4g1 to 4i1). In the basement horst and graben basin, the pre-existing faults experienced some degree of rotation and distortion due to the compressive stress. At this stage, the pre-existing faults in both the basement and the sedimentary layers (cover layers) were activated simultaneously, causing the volumetric strain and shear strain in the sedimentary layers to continue increasing. When the compressive displacement reached 20%, significant drag folding structural deformation occurred in the sedimentary layers above the large basement fault, while the cover layers inside the graben basin exhibited an overall folding pattern, forming a “folded above and faulted below” and “convex above and concave below” structural morphology.

Figure 4: Evolution of Structural Normal-Reverse Inversion in a Graben Basin Based on Discrete Element Simulation
Figure 4: Evolution of Structural Normal-Reverse Inversion in a Graben Basin Based on Discrete Element Simulation

Implications of the Research Results for the Maodong Graben Zone

  • The Maodong Graben Basin exhibits characteristics of both normal and reverse structural inversion
  • The reverse fault inversion controls the occurrence of earthquakes in this region.

Acknowledgements The numerical simulation experiment in this paper was conducted using the RICEBAL discrete element computational program and the subsequent stress-strain analysis code provided by Professor Julia Morgan from Rice University, USA. Constructive feedback and suggestions were provided by two anonymous reviewers during the manuscript review process. We would like to express our sincere gratitude to them.

The same experiment can also be conducted using VBOX

References

  1. Wu, Z.Y., Yin, H.W., Li, C.S., et al. Experimental Simulation of Normal-Reverse Inversion in Graben Basins and Its Implications for the Maodong Fault Zone.Journal of Nanjing University (Natural Sciences),2019,5(55):869-878.

Translator: Chi Yu