Geology topics

Kathryn Zerbe Materna

Publications and source records attributed to Kathryn Zerbe Materna.

6 recordsLinked to original sources

USGS2024

Relatively stable pressure effects and time-increasing thermal contraction control Heber geothermal field deformation

Due to geological complexities and observational gaps, it is challenging to identify the governing physical processes of geothermal field deformation including ground subsidence and earthquakes. In the west and east regions of the Heber Geothermal Field (HGF), decade-long subsidence was occurring despite injection of heat-depleted brines, along with transient reversals between uplift and subsidence. These observed phenomena contradict current knowledge that injection leads to surface uplift. Here we show that high-yield production wells at the HGF center siphon fluid from surrounding regions, which can cause subsidence at low-rate injection locations. Moreover, the thermal contraction effect by cooling increases with time and eventually overwhelms the pressure effects of pressure fluctuation and poroelastic responses, which keep relatively stable during geothermal operations. The observed subsidence anomalies result from the siphoning effect and thermal contraction. We further demonstrate that thermal contraction dominates long-term trends of surface displacement and seismicity growth, while pressure effects drive near-instantaneous changes.

California↗

USGS2024

Subduction intraslab-interface fault interactions in the 2022 Mw 6.4 Ferndale, California earthquake sequence

The Mendocino triple junction, the intersection of the Pacific, North American, and Gorda plates, activates a collection of disparate faults that reconcile Cascadia subduction with San Andreas transform motion. The December 20, 2022, Mw 6.4 Ferndale, California earthquake occurred within this complex zone as strike-slip faulting within the subducting Gorda slab. Here, we analyze the seismic and geodetic signatures of the mainshock and aftershock sequence to illuminate its role within complex tectonic surroundings. We find aftershocks on varied fault structures within the uppermost Gorda slab, yet seismicity on the subduction interface itself was notably absent. Nevertheless, we identify small but coherent postseismic deformation that is well modeled by aseismic slip on this interface, likely triggered by stresses generated at the updip limit of coseismic rupture. This sequence demonstrates the potential for interactions between intra-slab earthquakes and slip on the subduction megathrust, highlighting the need to consider this and other subduction zones as coupled systems of interacting faults.

California↗

USGS2023

Slip deficit rates on southern Cascadia faults resolved with viscoelastic earthquake cycle modeling of geodetic deformation

The fore‐arc of the southern Cascadia subduction zone (CSZ), north of the Mendocino triple junction (MTJ), is home to a network of Quaternary‐active crustal faults that accumulate strain due to the interaction of the North American, Juan de Fuca (Gorda), and Pacific plates. These faults, including the Little Salmon and Mad River fault (LSF and MRF) zones, are located near the most populated parts of California’s north coast and show paleoseismic evidence for three slip events of several‐meter scale in the past 1700 yr. However, the geodetic slip rates of these faults are poorly constrained. In this work, we analyze a new compilation of interseismic geodetic velocities from Global Navigation Satellite Systems, leveling, and tide gauge data near the MTJ to constrain present‐day slip deficit rates on upper‐plate faults and coupling on the megathrust. We construct Green’s functions for interseismic slip deficit for discrete faults embedded in an elastic plate overlying a viscoelastic mantle. We then use a constrained least‐squares inversion to determine best‐fitting slip rates on the major faults and investigate slip rate trade‐offs between faults. Results indicate that the LSF and MRF systems together accumulate 4–5 mm/yr of reverse‐slip deficit, although their separate slip rates cannot be determined independently. Modeling of the horizontal and vertical velocities suggests that the southernmost CSZ is coupled interseismically to deeper than 25 km depth. We also find that 6–17 mm/yr of right‐lateral slip deficit extends north of the MTJ and into the southern Cascadia fore‐arc. These results reinforce the notion that both the southernmost Cascadia megathrust and the smaller fore‐arc faults above it contribute to regional seismic hazard.

California, Oregon, Washington↗

USGS2023

Quantification of geodetic strain rate uncertainties and implications for seismic hazard estimates

Geodetic velocity data provide first-order constraints on crustal surface strain rates, which in turn are linked to seismic hazard. Estimating the 2-D surface strain tensor everywhere requires knowledge of the surface velocity field everywhere, while geodetic data such as Global Navigation Satellite System (GNSS) only have spatially scattered measurements on the surface of the Earth. To use these data to estimate strain rates, some type of interpolation is required. In this study, we review methodologies for strain rate estimation and compare a suite of methods, including a new implementation based on the geostatistical method of kriging, to compare variation between methods with uncertainty based on one method. We estimate the velocity field and calculate strain rates in southern California using a GNSS velocity field and five different interpolation methods to understand the sources of variability in inferred strain rates. Uncertainty related to data noise and station spacing (aleatoric uncertainty) is minimal where station spacing is dense and maximum far from observations. Differences between methods, related to epistemic uncertainty, are usually highest in areas of high strain rate due to differences in how gradients in the velocity field are handled by different interpolation methods. Parameter choices, unsurprisingly, have a strong influence on strain rate field, and we propose the traditional L -curve approach as one method for quantifying the inherent trade-off between fit to the data and models that are reflective of tectonic strain rates. Doing so, we find total variability between five representative strain rate models to be roughly 40 per cent, a much lower value than roughly 100 per cent that was found in previous studies (Hearn et al .). Using multiple methods to tune parameters and calculate strain rates provides a better understanding of the range of acceptable models for a given velocity field. Finally, we present an open-source Python package (Materna et al .) for calculating strain rates, Strain_2D, which allows for the same data and model grid to be used in multiple strain rate methods, can be extended with other methods from the community, and provides an interface for comparing strain rate models, calculating statistics and estimating strain rate uncertainty for a given GNSS data set.

Geophysical Journal International↗

USGS2023

Dense geophysical observations reveal a triggered, concurrent multi-fault rupture at the Mendocino Triple Junction

A central question of earthquake science is how far ruptures can jump from one fault to another, because cascading ruptures can increase the shaking of a seismic event. Earthquake science relies on earthquake catalogs and therefore how complex ruptures get documented and cataloged has important implications. Recent investments in geophysical instrumentation allow us to resolve increasingly complex, multi-fault ruptures for even moderate-sized earthquakes. We combine dense seismic and geodetic measurements to reveal an enigmatic rupture in late 2021 at the Mendocino Triple Junction in northern California. We show that rupture was dynamically triggered, yet concurrent, on two distinct faults roughly 30 km apart. Thus, this rupture combines features of complex ruptures usually considered to be single earthquakes, and triggered ruptures considered as multiple earthquakes. This event illustrates that moderate-sized earthquakes can exhibit similar complexity to that more commonly documented for large earthquakes.

California↗

USGS2022

Detection of aseismic slip and poroelastic reservoir deformation at the North Brawley Geothermal Field from 2009 to 2019

The North Brawley Geothermal Field, located within the Brawley Seismic Zone of Southern California, presents a case study for understanding seismic hazards linked to fluid injection and geothermal energy extraction. An earthquake swarm near the geothermal field in 2012 included two earthquakes with magnitudes greater than 5 and was potentially preceded by a years-long aseismic slip transient. To better understand ground deformation around the geothermal field, including its evolution with time and its physical mechanisms, we analyze deformation before, during, and after the swarm using ground- and satellite-based geodetic techniques between 2009 and 2019. We integrate observations from GNSS, Sentinel-1, TerraSAR-X, UAVSAR, and leveling surveys into a single deformation history. Modeling of this new collection of observations at the North Brawley Geothermal Field provides evidence for 80% more pre-swarm aseismic slip than previously recognized from 2009 to 2012. During the 2012 Brawley swarm, our geodetic slip inversions closely match the results of seismic waveform inversions from the swarm events. After the 2012 swarm, surface deformation is dominated by poroelastic deformation of a shallow fluid reservoir at <1 km depth rather than fault slip. The deformation history and seismicity catalogs at North Brawley suggest a cessation of fault-related slip during the ∼7 years after the 2012 earthquake swarm.

California↗