skip to main content


Title: An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data
SUMMARY

Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

 
more » « less
NSF-PAR ID:
10478364
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Volume:
236
Issue:
2
ISSN:
0956-540X
Format(s):
Medium: X Size: p. 778-797
Size(s):
["p. 778-797"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Our understanding of plate boundary deformation has been enhanced by transient signals observed against the backdrop of time‐independent secular motions. We make use of a new analysis of displacement time series from about 1,000 continuous Global Positioning System (GPS) stations in California from 1999 to 2018 to distinguish tectonic and nontectonic transients from secular motion. A primary objective is to define a high‐resolution three‐dimensional reference frame (datum) for California that can be rapidly maintained with geodetic data to accommodate both secular and time‐dependent motions. To this end, we compare the displacements to those predicted by a horizontal secular fault slip model for the region and construct displacement and strain rate fields. Over the past 19 years, California has experienced 19 geodetically detectable earthquakes and widespread postseismic deformation. We observe postseismic strain rate variations as large as 1,000 nstrain/year with moment releases equivalent up to an Mw6.8 earthquake. We find significant secular differences up to 10 mm/year with the fault slip model, from the Mendocino Triple Junction to the southern Cascadia subduction zone, the northern Basin and Range, and the Santa Barbara channel. Secular vertical uplift is observed across the Transverse Ranges, Coastal Ranges, Sierra Nevada, as well as large‐scale postseismic uplift after the 1999 Mw7.1 Hector Mine and 2010 Mw7.2 El Mayor‐Cucapah earthquakes. We also identify areas of vertical land motions due to anthropogenic, natural, and magmatic processes. Finally, we demonstrate the utility of the kinematic datum by improving the accuracy of high‐spatial‐resolution 12‐day repeat‐cycle Sentinel‐1 Interferometric Synthetic Aperture Radar displacement and velocity maps.

     
    more » « less
  2. null (Ed.)
    Slow slip events (SSEs) at the northern Hikurangi subduction margin, New Zealand, are among the best-documented shallow SSEs on Earth. International Ocean Discovery Program Expeditions 372 and 375 were undertaken to investigate the processes and in situ conditions that underlie subduction zone SSEs at the northern Hikurangi Trough. We accomplished this goal by (1) coring and geophysical logging at four sites, including penetration of an active thrust fault (the Pāpaku fault) near the deformation front, the upper plate above the SSE source region, and the incoming sedimentary succession in the Hikurangi Trough and atop the Tūranganui Knoll seamount; and (2) installing borehole observatories in the Pāpaku fault and in the upper plate overlying the slow slip source region. Logging-while-drilling (LWD) data for this project were acquired as part of Expedition 372, and coring, wireline logging, and observatory installations were conducted during Expedition 375. Northern Hikurangi subduction margin SSEs recur every 1–2 y and thus provide an ideal opportunity to monitor deformation and associated changes in chemical and physical properties throughout the slow slip cycle. In situ measurements and sampling of material from the sedimentary section and oceanic basement of the subducting plate reveal the rock properties, composition, lithology, and structural character of material that is transported downdip into the SSE source region. A recent seafloor geodetic experiment raises the possibility that SSEs at northern Hikurangi may propagate to the trench, indicating that the shallow thrust fault (the Pāpaku fault) targeted during Expeditions 372 and 375 may also lie in the SSE rupture area and host a portion of the slip in these events. Hence, sampling and logging at this location provides insights into the composition, physical properties, and architecture of a shallow fault that may host slow slip. Expeditions 372 and 375 were designed to address three fundamental scientific objectives: 1. Characterize the state and composition of the incoming plate and shallow fault near the trench, which comprise the protolith and initial conditions for fault zone rock at greater depth and which may itself host shallow slow slip; 2. Characterize material properties, thermal regime, and stress conditions in the upper plate directly above the SSE source region; and 3. Install observatories in the Pāpaku fault near the deformation front and in the upper plate above the SSE source to measure temporal variations in deformation, temperature, and fluid flow. The observatories will monitor volumetric strain (via pore pressure as a proxy) and the evolution of physical, hydrological, and chemical properties throughout the SSE cycle. Together, the coring, logging, and observatory data will test a suite of hypotheses about the fundamental mechanics and behavior of SSEs and their relationship to great earthquakes along the subduction interface. 
    more » « less
  3. Abstract

    Large earthquakes and tsunamis in subduction zone forearcs occur via slip on the shallow plate boundary and upper plate faults, but the locations, geometries, and slip histories of these faults can be difficult to constrain in regions with minimal subsurface geophysical and stratigraphic data. Here, we test a new approach to quantify the submarine seafloor geomorphic response to forearc deformation in order to identify structures that contribute to active deformation, to interpret their geometry and kinematics, and to evaluate their relative rates, magnitudes, and timing of deformation. We develop a workflow that uses filtered bathymetric digital elevation models, where long wavelength topography has been removed, to isolate the slope, relief, curvature, ridgelines, and trough lines associated with faults, fault‐related folds, and slope failures. We apply these methods to the Kumano region of the Nankai accretionary prism, southeastern Japan, where existing constraints on fault geometry, kinematics, and deformation history allow us to both evaluate the efficacy of our approach and to identify the lateral continuity of deformation processes. Our bathymetric analyses yield a high‐resolution tectono‐geomorphic map of active structures and reveal along strike variations in strain accumulation and out‐of‐sequence deformation. These metrics also demonstrate the importance of a strike‐slip fault system at the seaward edge of the Kumano Basin as a primary structure that accommodates deformation and partitions strain in the Nankai forearc. These results show the utility of using a submarine tectono‐geomorphic approach to evaluate active deformation in forearcs, particularly in regions with limited geophysical and core data.

     
    more » « less
  4. null (Ed.)
    Slow slip events (SSEs) at the northern Hikurangi subduction margin, New Zealand, are among the best-documented shallow SSEs on Earth. International Ocean Discovery Program Expedition 375 was undertaken to investigate the processes and in situ conditions that underlie subduction zone SSEs at the northern Hikurangi Trough by (1) coring at four sites, including an active fault near the deformation front, the upper plate above the high-slip SSE source region, and the incoming sedimentary succession in the Hikurangi Trough and atop the Tūranganui Knoll Seamount, and (2) installing borehole observatories in an active thrust near the deformation front and in the upper plate overlying the slow slip source region. Logging-while-drilling (LWD) data for this project were acquired as part of Expedition 372 (26 November 2017–4 January 2018; see the Expedition 372 Preliminary Report for further details on the LWD acquisition program). Northern Hikurangi subduction margin SSEs recur every 1–2 years and thus provide an ideal opportunity to monitor deformation and associated changes in chemical and physical properties throughout the slow slip cycle. Sampling of material from the sedimentary section and oceanic basement of the subducting plate reveals the rock properties, composition, lithology, and structural character of material that is transported downdip into the SSE source region. A recent seafloor geodetic experiment raises the possibility that SSEs at northern Hikurangi may propagate all the way to the trench, indicating that the shallow thrust fault zone targeted during Expedition 375 may also lie in the SSE rupture area. Hence, sampling at this location provides insights into the composition, physical properties, and architecture of a shallow fault that may host slow slip. Expedition 375 (together with the Hikurangi subduction LWD component of Expedition 372) was designed to address three fundamental scientific objectives: (1) characterize the state and composition of the incoming plate and shallow plate boundary fault near the trench, which comprise the protolith and initial conditions for fault zone rock at greater depth and which may itself host shallow slow slip; (2) characterize material properties, thermal regime, and stress conditions in the upper plate above the core of the SSE source region; and (3) install observatories at an active thrust near the deformation front and in the upper plate above the SSE source to measure temporal variations in deformation, temperature, and fluid flow. The observatories will monitor volumetric strain (via pore pressure as a proxy) and the evolution of physical, hydrological, and chemical properties throughout the SSE cycle. Together, the coring, logging, and observatory data will test a suite of hypotheses about the fundamental mechanics and behavior of SSEs and their relationship to great earthquakes along the subduction interface. 
    more » « less
  5. null (Ed.)
    Slow slip events (SSEs) at the northern Hikurangi subduction margin, New Zealand, are among the best-documented shallow SSEs on Earth. International Ocean Discovery Program Expedition 375 aims to investigate the processes and in situ conditions that underlie subduction zone SSEs at northern Hikurangi through coring of the frontal thrust, upper plate, and incoming sedimentary succession and through installation of borehole observatories in the frontal thrust and upper plate above the slow slip source area. Logging-while-drilling (LWD) data for this project will be acquired as part of Expedition 372 (beginning in November 2017; see the Expedition 372 Scientific Prospectus for further details on the LWD acquisition program). Northern Hikurangi subduction margin SSEs recur every 2 years and thus provide an excellent setting to monitor deformation and associated chemical and physical properties surrounding the SSE source area throughout the slow slip cycle. Sampling material from the sedimentary section and oceanic basement of the subducting plate and from the primary active thrust in the outer wedge near the trench will reveal the rock properties, composition, and lithologic and structural character of the material transported downdip to the known SSE source region. A recent seafloor geodetic experiment shows the possibility that SSEs at northern Hikurangi may propagate all the way to the trench, indicating that the shallow fault zone target for Expedition 375 may lie within the SSE rupture area. Four primary sites are planned for coring, and observatories will be installed at two of these sites. Expedition 375 (together with the Hikurangi subduction component of Expedition 372) is designed to address three fundamental scientific objectives: (1) characterize the state and composition of the incoming plate and shallow plate boundary fault near the trench, which comprise the protolith and initial conditions for fault zone rock at greater depth; (2) characterize material properties, thermal regime, and stress conditions in the upper plate above the SSE source region; and (3) install observatories at the frontal thrust and in the upper plate above the SSE source to measure temporal variations in deformation, fluid flow, and seismicity. The observatories will monitor deformation and the evolution of physical, hydrological, and chemical properties throughout the SSE cycle. Together, the coring, logging, and observatory data will test a suite of hypotheses about the fundamental mechanics and behavior of slow slip events and their relationship to great earthquakes along the subduction interface. 
    more » « less