skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Lithology and Fault‐Related Stress Variations Along the TCDP Boreholes: The Stress State Before and After the 1999 Chi‐Chi Earthquake
Abstract Understanding the stress state before and after an earthquake is essential to study how stress on faults evolves during the seismic cycle. This study integrates wellbore failure analysis, laboratory experiments, and edge dislocation model to study the stress state before and after the Chi‐Chi earthquake. The post‐earthquake in‐situ stress state observed along boreholes of the Taiwan Chelungpu‐fault Drilling Project (TCDP) is heterogeneous due to lithological variations. Along the borehole, we observe that drilling‐induced tensile fractures are only present in sandstones, whereas breakouts are mostly present in silt‐rich rocks. Laboratory experiments on TCDP cores also show that tensile and compressive strength are weaker in sandstones than in silt‐rich rocks. These observations imply that both maximum and minimum horizontal principal stresses are higher in silt‐rich intervals. Extended leak‐off tests in the TCDP borehole also show lower minimum horizontal stress in sand‐rich intervals, consistent with the above observations. We combine these observations to estimate a profile of stress magnitudes along the well which explains the variability of stress states found in previous studies. The stress heterogeneity we observed underlines the importance of acknowledging the spatial scale that the stress data represent. We then use an edge dislocation model constrained by GPS surface displacements obtained during Chi‐Chi earthquake to calculate the coseismic stress changes. Our inferred pre‐earthquake stress magnitudes, obtained by subtracting the coseismic stress change from the post‐earthquake stress, suggest subcritical stress state before the earthquake despite the large displacements observed during the Chi‐Chi earthquake in the region where TCDP encountered the fault.  more » « less
Award ID(s):
1727661 1829597
PAR ID:
10367732
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
127
Issue:
2
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Following large earthquakes, viscoelastic stress relaxation may contribute to postseismic deformation observed at Earth's surface. Mechanical representations of viscoelastic deformation require a constitutive relationship for the lower crust/upper mantle material where stresses are diffused and, for non‐linear rheologies, knowledge of absolute stress level. Here, we describe a kinematic approach to representing geodetically observed postseismic motions that does not require an assumed viscoelastic rheology. The core idea is to use observed surface motions to constrain time‐dependent displacement boundary conditions applied at the base of the elastic upper crust by viscoelastic motions in the lower crust/upper mantle, approximating these displacements as slip on a set of dislocation elements. Using three‐dimensional forward models of viscoelastically modulated postseismic deformation in a thrust fault setting, we show how this approach can accurately represent surface motions and recover predicted displacements at the base of the elastic layer. Applied to the 1999 Chi‐Chi (Taiwan) earthquake, this kinematic approach can reproduce geodetically observed displacements and estimates of the partitioning between correlated postseismic deformation mechanisms. Specifically, we simultaneously estimate afterslip on the earthquake source fault that is similar to previous estimates, along with slip on dislocations at the base of the elastic layer that mimic predictions from viscous stress dissipation models in which viscosity is inferred to vary three‐dimensionally. A use case for the dislocation approach to modeling viscoelastic deformation is the estimation of spatiotemporally variable fault slip processes, including across sequential interseismic phases of the earthquake cycle, without assuming a lower crust/upper mantle rheology. 
    more » « less
  2. Abstract Quantifying off-fault deformation in the near field remains a challenge for earthquake monitoring using geodetic observations. We propose an automated change detection strategy using geometric primitives generated using a deep neural network, random sample consensus and least squares adjustment. Using mobile laser scanning point clouds of vineyards acquired after the magnitude 6.0 2014 South Napa earthquake, our results reveal centimeter-level horizontal ground deformation over three kilometers along a segment of the West Napa Fault. A fault trace is detected from rows of vineyards modeled as planar primitives from the accumulated coseismic response, and the postseismic surface displacement field is revealed by tracking displacements of vineyard posts modeled as cylindrical primitives. Interpreted from the detected changes, we summarized distributions of deformation versus off-fault distances and found evidence of off-fault deformation. The proposed framework using geometric primitives is shown to be accurate and practical for detection of near-field off-fault deformation. 
    more » « less
  3. null (Ed.)
    The fault damage zone is a well-known structure of localized deformation around faults. Its material properties evolve over earthquake cycles due to coseismic damage accumulation and interseismic healing. We will present fully dynamic earthquake cycle simulations to show how the styles of earthquake nucleation and rupture propagation change as fault zone material properties vary temporally. First, we will focus on the influence of fault zone structural maturity quantified by near-fault seismic wave velocities in simulations. The simulations show that immature fault zones promote small and moderate subsurface earthquakes with irregular recurrence intervals, whereas mature fault zones host pulse-like earthquake rupture that can propagate to the surface, extend throughout the seismogenic zone, and occur at regular intervals. The interseismic healing in immature fault zones plays a key role in allowing the development of aseismic slip episodes including slow-slip events and creep, which can propagate into the seismogenic zone, and thus limit the sizes of subsequent earthquakes by releasing fault stress. In the second part, we will discuss how the precursory changes of seismic wave velocities of fault damage zones may affect earthquake nucleation process. Both laboratory experiments and seismic observations show that the abrupt earthquake failure can be preceded by accelerated fault deformation and the accompanying velocity reduction of near-fault rocks. We will use earthquake cycle simulations to systematically test the effects of timing and amplitudes of such precursory velocity changes. Our simulations will provide new insights into the interplay between fault zone structure and earthquake nucleation process, which can be used to guide future real-time monitoring of major fault zones. 
    more » « less
  4. null (Ed.)
    The Nankai Trough Seismogenic Zone Experiment (NanTroSEIZE) program is a coordinated, multiexpedition drilling project designed to investigate fault mechanics and seismogenesis along subduction megathrusts through direct sampling, in situ measurements, and long-term monitoring in conjunction with allied laboratory and numerical modeling studies. The fundamental scientific objectives of the NanTroSEIZE project include characterizing the nature of fault slip and strain accumulation, fault and wall rock composition, fault architecture, and state variables throughout the active plate boundary system. Site C0002 is in the Kumano forearc basin above the seismogenic, and presumably locked, portion of the plate boundary thrust system. The Kumano Basin sedimentary sequence and uppermost part of the accretionary prism were drilled, logged, and sampled during Integrated Ocean Drilling Program Expeditions 314 (logging while drilling [LWD] to 1401.5 mbsf), 315 (coring to 1057 mbsf), 338 (LWD to 2005 mbsf and coring to 1120 mbsf), and 348 (LWD to 3058.5 mbsf, with limited coring from 2163-2218.5 mbsf). International Ocean Discovery Program (IODP) Expedition 358 aims to reach and sample the megasplay fault/plate boundary fault at Site C0002 by extending the riser borehole (Hole C0002P) established during previous Integrated Ocean Drilling Program NanTroSEIZE expeditions. Evidence suggests that the drilling target at the megasplay will reach a region where megathrust seismogenic processes are active. A new borehole will be “kicked off” from Hole C0002P and will be extended to ~5000 meters below seafloor (mbsf) and then across the high-amplitude seismic reflector identified as the main plate boundary fault. Continuous LWD with drilling mud gas analysis and limited coring at the anticipated plate boundary fault depth, will leave the cased borehole completed at ~5200 mbsf. A suite of analyses on cuttings, mud gases, and limited cores will address the four primary scientific objectives: (1) determine the composition, stratigraphy, and deformational history of the Miocene accretionary prism; (2) reconstruct its thermal, diagenetic, and metamorphic history; (3) determine horizontal stress orientations and magnitudes; and (4) investigate the mechanical and hydrological properties of the upper plate of the seismogenic plate boundary. The main scientific objective is to log and sample the hanging wall, the fault zone, and into the footwall. The main contingency plan is to leave a cased borehole in good condition for future installation of a long-term borehole monitoring system (LTBMS). These operations will extend drilling conducted during Integrated Ocean Drilling Program Expeditions 326, 338, and 348. The entire expedition will cover a period of 164 days, beginning on 7 October 2018 and ending on 21 March 2019. 
    more » « less
  5. Abstract We study the mechanical response of two‐dimensional vertical strike‐slip fault to coseismic damage evolution and interseismic healing of fault damage zones by simulating fully dynamic earthquake cycles. Our models show that fault zone structure evolution during the seismic cycle can have pronounced effects on mechanical behavior of locked and creeping fault segments. Immature fault damage zone models exhibit small and moderate subsurface earthquakes with irregular recurrence intervals and abundance of slow‐slip events during the interseismic period. In contrast, mature fault damage zone models host pulse‐like earthquake ruptures that can propagate to the surface and extend throughout the seismogenic zone, resulting in large stress drop, characteristic rupture extents, and regular recurrence intervals. Our results suggest that interseismic healing and coseismic damage accumulation in fault zones can explain the observed differences of earthquake behaviors between mature and immature fault zones and indicate a link between regional seismic hazard and fault structural maturity. 
    more » « less