Fault creep reduces seismic hazard and serves as a window into plate boundary processes; however, creep rates are typically constrained with sparse measurements. We use differential lidar topography (11–13 year time span) to measure a spatially dense surface deformation field along a 150 km section of the Central San Andreas and Calaveras faults. We use an optimized windowed‐iterative‐closest‐point approach to resolve independent creep rates every 400 m at 1–2 km apertures. Rates vary from <10 mm/year along the creeping fault ends to over 30 mm/year along much of the central 100 km of the fault. Creep rates are 3–8 mm/year higher than most rates from alignment arrays and creepmeters, likely due to the larger aperture of the topographic differencing. Creep is often focused along discrete fault traces, but strain is sometimes distributed in areas of complex fault geometry, such as Mustang Ridge. Our observations constrain shallow seismic moment accumulation and the location of the creeping fault trace.
Imaging tectonic creep along active faults is critical for measuring strain accumulation and ultimately understanding the physical processes that guide creep and the potential for seismicity. We image tectonic deformation along the central creeping section of the San Andreas Fault at the Dry Lake Valley paleoseismic site (36.468°N, 121.055°W) using three data sets with varying spatial and temporal scales: (1) an Interferometric Synthetic Aperture Radar (InSAR) velocity field with an ~100‐km footprint produced from Sentinel‐1 satellite imagery, (2) light detection and ranging (lidar) and structure‐from‐motion 3‐D topographic differencing that resolves a decade of deformation over a 1‐km aperture, and (3) surface fractures that formed over the 3‐ to 4‐m wide fault zone during a drought from late 2012 to 2014. The InSAR velocity map shows that shallow deformation is localized to the San Andreas Fault. We demonstrate a novel approach for differencing airborne lidar and structure‐from‐motion topography that facilitates resolving deformation along and adjacent to the San Andreas Fault. The 40‐m resolution topographic differencing resolves a 2.5 ± 0.2 cm/yr slip rate localized to the fault. The opening‐mode fractures accommodate
- Award ID(s):
- 2054160
- NSF-PAR ID:
- 10381333
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 125
- Issue:
- 10
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Observations of surface deformation within 1–2 km of a surface rupture contain invaluable information about the coseismic behavior of the shallow crust. We investigate the oblique strike‐slip 2016 M7 Kumamoto, Japan, earthquake, which ruptured the Futagawa‐Hinagu Fault. We solve for variable fault slip in an inversion of differential lidar topography, satellite optical image correlation, and Interferometric Synthetic Aperture Radar (InSAR)‐derived surface displacements. The near‐fault differential lidar pose several challenges. The model fault geometry must follow the surface trace at the sub‐kilometer scale. Integration of displacement datasets with different sensitivities to the 3D deformation field and varying spatial distribution permits additional complexity in the inferred slip but introduces ambiguity that requires careful selection of the regularization. We infer a Mw
earthquake. The maximum slip of 6.9 m occurred at 4.5‐km depth, suggesting an on‐fault slip deficit in the upper several kilometers of the crust that likely reflects distributed and inelastic deformation within the shallow fault zone. -
more » « less
Abstract Measuring crustal strain and seismic moment accumulation, is crucial for understanding the growth and distribution of seismic hazards along major fault systems. Here, we develop a methodology to integrate 4.5 years (2015–2019.5) of Sentinel‐1 Interferometric Synthetic Aperture Radar (InSAR) and continuous Global Navigation Satellite System (GNSS) time series to achieve 6 to 12‐day sampling of surface displacements at ∼500 m spatial resolution over the entire San Andreas fault system. Numerous interesting deformation signals are identified with this product (video link:
https://www.youtube.com/watch?v=SxNLQKmHWpY ). We decompose the line‐of‐sight InSAR displacements into three dimensions by combining the deformation azimuth from a GNSS‐derived interseismic fault model. We then construct strain rate maps using a smoothing interpolator with constraints from elasticity. The resulting deformation field reveals a wide array of crustal deformation processes including, on‐ and off‐fault secular and transient tectonic deformation, creep rates on all the major faults, and vertical signals associated with hydrological processes. The strain rate maps show significant off‐fault components that were not captured by GNSS‐only models. These results are important in assessing the seismic hazard in the region. -
more » « less
Abstract The Calaveras Fault (CF) branches from the San Andreas Fault (SAF) near San Benito, extending sub‐parallel to the SAF for about 50 km with only 2–6 km separation and diverging northeastward. Both the SAF and CF are partially coupled, exhibit spatially variable aseismic creep and have hosted moderate to large earthquakes in recent decades. Understanding how slip partitions among the main fault strands of the SAF system and establishing their degree of coupling is crucial for seismic hazard evaluation. We perform a timeseries analysis using more than 5 years of Sentinel‐1 data covering the Bay Area (May 2015–October 2020), specifically targeting the spatiotemporal variations of creep rates around the SAF‐CF junction. We derive the surface creep rates from cross‐fault InSAR timeseries differences along the SAF and CF including adjacent Sargent and Quien Sabe Faults. We show that the variable creep rates (0–20 mm/yr) at the SAF‐CF junction are to first order controlled by the angle between the fault strike and the background stress orientation. We further examine the spatiotemporal variation of creep rates along the SAF and CF and find a multi‐annual coupling increase during 2016–2018 the subparallel sections of both faults, with the CF coupling change lagging behind the SAF by 3–6 months. Similar temporal variations are also observed in both
b ‐values inferred from declustered seismicity and aseismic slip rates inferred from characteristic repeating earthquakes. -
more » « less
Abstract Observations of shallow fault creep reveal increasingly complex time‐dependent slip histories that include quasi‐steady creep and triggered as well as spontaneous accelerated slip events. Here we report a recent slow slip event on the southern San Andreas fault triggered by the 2017
M w 8.2 Chiapas (Mexico) earthquake that occurred 3,000 km away. Geodetic and geologic observations indicate that surface slip on the order of 10 mm occurred on a 40‐km‐long section of the southern San Andreas fault between the Mecca Hills and Bombay Beach, starting minutes after the Chiapas earthquake and continuing for more than a year. Both the magnitude and the depth extent of creep vary along strike. We derive a high‐resolution map of surface displacements by combining Sentinel‐1 Interferometric Synthetic Aperture Radar acquisitions from different lines of sight. Interferometric Synthetic Aperture Radar‐derived displacements are in good agreement with the creepmeter data and field mapping of surface offsets. Inversions of surface displacement data using dislocation models indicate that the highest amplitudes of surface slip are associated with shallow (<1 km) transient slip. We performed 2‐D simulations of shallow creep on a strike‐slip fault obeying rate‐and‐state friction to constrain frictional properties of the top few kilometers of the upper crust that can produce the observed behavior.
