Geology topics

David A. Lockner

Publications and source records attributed to David A. Lockner.

At least 19 recordsLinked to original sources

USGS2026

Velocity-independent dry friction on mica: A realization of ideal Amontons-Coulomb friction

The Amontons-Coulomb friction law assumes that the frictional force between materials is independent of sliding velocity. However, as Coulomb noted, this is a rough approximation, and a second-order dependence of friction on the logarithm of sliding velocity is incorporated in a commonly used ‘rate- and state-dependent’ friction representation. Here we conduct shear experiments on mica, a layer-structured mineral, at temperatures ranging from 25 to 200ºC and under normal stress of 100 MPa. The friction coefficient clearly depends on the logarithmic sliding velocity at 25ºC, but rate sensitivity decreases with increasing temperature until at 200ºC, the friction coefficient is independent of sliding velocity. Our findings could initiate the development of velocity-independent frictional materials, realizing the ideal Amontons-Coulomb friction.

Physical Review Letters↗

USGS2026

Earthquake rupture arrest from depth‐dependent frictional stability on the Pütürge segment of the East Anatolian Fault Zone, Türkiye (Turkey)

Determining why earthquake ruptures stop is a central challenge in earthquake science and seismic hazard assessment. The Pütürge segment of the East Anatolian Fault Zone, Türkiye, exhibits shallow creep (<∼3 km depth) yet hosts M > 6.5 earthquake ruptures at greater depth. Here, we evaluate whether variations in frictional stability along this segment aided arrest of the 2020 M 6.7 Elaziğ and 2023 M 7.8 Pazarcık earthquake ruptures. Analysis of Sentinel‐1 Synthetic Aperture Radar imagery indicates the 2023 M 7.8 rupture propagated laterally into a metamorphic massif within the Pütürge segment, where slip rapidly decayed below detection limits. Creepmeters along this segment recorded no significant surface afterslip (<3 mm) in the following year. To investigate this fault‐slip behavior, we conducted triaxial friction experiments on Pütürge fault gouge sampled from an outcrop exposure. The gouge, composed primarily of muscovite, quartz, and calcite, is velocity strengthening at conditions approximating 0–2.5 km depth and velocity weakening at 4–5 km depth. This transition to velocity‐weakening friction is associated with enhanced comminution and shear localization observed microstructurally. Our results suggest that depth‐dependent frictional stability of the Pütürge fault segment facilitates rupture nucleation and propagation at depth while maintaining rupture‐arresting behavior near Earth’s surface.

Pütürge segment of the East Anatolian Fault Zone↗

USGS2025

Comparative properties of saponitic fault gouge and serpentinite muds cored from mud volcanoes of the Mariana subduction zone

We obtained 12 core samples for physical and chemical characterization from three serpentinite mud volcanoes (Yinazao, Asùt Tesoru, and Fantangisña) located on the forearc of the Mariana subduction system, that were drilled during International Ocean Discovery Program Expedition 366. Two samples from the Fantangisña mud volcano are interpreted to be clay-rich fault gouges derived from the subduction channel. Their bulk compositions are intermediate between the serpentinites and oceanic basalts. The oceanic crustal materials in the gouges have been thoroughly metasomatized and the serpentinites extensively altered to the trioctahedral, Mg-rich smectite clays saponite and corrensite. The only relict phases in clasts of crustal rock are accessory Ti- and P-bearing minerals. The two fault gouge samples have lower frictional strengths ( μ < 0.2) than the serpentinites ( μ = 0.2–0.4), and their measured permeabilities are also somewhat lower. Their physical and compositional properties correspond to saponitic gouges from other faults that juxtapose serpentinite against crustal rocks, in particular gouges from the two creeping traces of the San Andreas Fault recovered in the core from the San Andreas Fault Observatory at Depth. The décollement beneath Fantangisña mud volcano is thus expected to be very weak and likely characterized by stable slip.

Geophysics, Geochemistry, Geosystems↗

USGS2025

Evolution of permeability and strength recovery of shear fracture under hydrothermal conditions

Geothermal energy is a clean and renewable resource that depends on the ability to move water through hot rock. In many locations, the ability to move water through rock requires the presence of extensive natural or human-made systems of fractures. However, these fracture systems are influenced by a variety of complex processes that occur at the temperature and pressure conditions found in geothermal reservoirs. To make geothermal systems more efficient, it is important to understand how these fractures evolve over time in response to these processes. This study explored how fractures in rock change under high-temperature and pressure conditions, like those found in geothermal reservoirs, through laboratory experiments and numerical simulations. In general fractures, tend to close due to pressure from the surrounding rock preventing the flow of water. One set of laboratory experiments revealed that fractures in granitic rock can weaken over time at elevated temperatures. These weak fractures can then slip slowly potentially keeping the fracture open for a longer period. This study also showed that at the highest temperature examined (250 °C) large differences between the chemical composition of the rock and the fluid moving through it could improve water flow temporarily likely due to the water dissolving parts of the rock. Flow rate predicted by models shows overall good agreement with the flow rate measured in our experiments but with some long time-scale variations that may be due to mechanical wear of the fracture surface. A second set of experiments found that, over the temperature range examined, the closure of fractures was caused by brittle failure driven by high stresses at points where fracture surfaces were in contact. This result suggests that fracture closure in geothermal reservoirs could be counteracted by maintaining high pore pressures. Additionally, we also tested a new method using electrical resistance to measure fracture closure in controlled systems with high precision.

Final Technical Report↗

USGS2023

Time-dependent weakening of granite at hydrothermal conditions

The evolution of a fault's frictional strength during the interseismic period is a critical component of the earthquake cycle, yet there have been relatively few studies that examine the time-dependent evolution of strength at conditions representative of seismogenic depths. Using a simulated fault in Westerly granite, we examined how frictional strength evolves under hydrothermal conditions up to 250°C during slide-hold-slide experiments. At temperatures ≤100°C, frictional strength generally increases with hold duration but, at 200 and 250°C, an initial increase in strength transitions to rapid time-dependent weakening for holds longer than 14 hr. Forward modeling of long hold periods at 250°C using the rate and state friction constitutive equations requires a second, strongly negative, state variable with a long evolution distance. This implies that significant hydrothermal alteration is occurring at 250°C, consistent with microstructural observations of dissolution and secondary mineral precipitation.

Geophysical Research Letters↗

USGS2023

Interlaboratory comparison of testing hydraulic, elastic, and failure properties in compression: Lessons learned

Many geoscientific problems require us to exploit synergies of experimental and numerical approaches, which in turn lead to questions regarding the significance of experimental details for validation of numerical codes. We report results of an interlaboratory comparison regarding experimental determination of mechanical and hydraulic properties of samples from five rock types, three sandstone varieties with porosities ranging from 5% to 20%, a marble, and a granite. The objective of this study was to build confidence in the participating laboratories’ testing approaches and to establish tractable standards for several physical properties of rocks. We addressed the issue of sample-to-sample variability by investigating the variability of basic physical properties of samples of a particular rock type and by performing repeat tests. Compressive strength of the different rock types spans an order of magnitude and shows close agreement between the laboratories. However, differences among stress–strain relations indicate that the external measurement of axial displacement and the determination of system stiffness require special attention, apparently more so than the external load measurement. Furthermore, post-failure behavior seems to exhibit some machine-dependence. The different methods used for the determination of hydraulic permeability, covering six orders of magnitude for the sample suite, yield differences in absolute values and pressure dependence for some rocks but not for others. The origin of the differences in permeability, in no case exceeding an order of magnitude, correlate with the compressive strength and potentially reflect a convolution of end plug–sample interaction, sample-to-sample variability, heterogeneity on sample scale, and/or anisotropy, the last two aspects are notably not accounted for by the applied evaluation procedures. Our study provides an extensive data set apt for “benchmarking” considerations, be it regarding new laboratory equipment or numerical modeling approaches.

Environmental Earth Sciences↗

USGS2023

Strength recovery in quartzite is controlled by changes in friction in experiments at hydrothermal conditions up to 200°C

The rate of fault zone restrengthening between earthquakes can be influenced by both frictional and cohesive healing processes. Friction is dependent on effective normal stress while cohesion is independent of normal stress, potentially explaining—in part—the lack of depth dependence of earthquake stress drops. Although amenable to laboratory testing, few studies have systematically addressed the normal stress dependence of restrengthening rate. This is partially due to difficulty in separating relative contributions of friction and cohesion in recovery of fault strength. We present results from a series of slide-hold-slide tests on thin layers (≤10 𝜇m) of ultrafine quartz gouge that develop during shearing of initially bare-surface quartzite. Tests were conducted at 10 MPa constant pore pressure, 20–200 MPa constant effective normal stress, and temperatures of 22°–200°C. Restrengthening, defined as the difference between peak shear stress measured after resumption of sliding and steady-state sliding shear stress, increases with the log of hold duration. The 200°C healing rate, 0.014 per e-fold increase in time, is comparable to that determined from seismological observations along the Calaveras Fault, California. Construction of Mohr-Coulomb failure envelopes shows that changes in cohesion are small (<1 MPa) and independent of hold durations to 10 5 s, indicating that the increased strength is due to changes in the friction coefficient. These experimental results are inconsistent with the hypothesis that cohesive healing explains the depth independence of earthquake stress drop, but higher temperatures, longer time-scales, and more complex mineralogy could facilitate cohesive healing in natural fault systems.

Journal of Geophysical Research - Solid Earth↗

USGS2023

Effect of thermal and mechanical processes on hydraulic transmissivity evolution

Fracture healing is a critical component of enhanced geothermal systems, the earthquake cycle, and induced seismicity. Accordingly, there is significant interest in understanding the process of healing and its effects on fluid transport. The creation, reactivation, and sustainability of fracture networks depend on complex coupling among thermal, hydraulic, mechanical, and chemical processes. We use laboratory slide-hold-slide experiments, at temperatures from 22 to 200 ˚C, to examine effects of fracture reactivation and quasi-static loading on the evolution of fluid transport properties of simulated fractures in Westerly granite. At all temperatures, the in-plane hydraulic transmissivity consistently decays during hold periods resulting in an overall reduction in transmissivity. During the first three to fifteen hours of an experiment, transmissivity decreases rapidly due to the generation of wear products, development of a sliding surface, and compaction of the resulting gouge. Once the sliding surface has developed, the long-term transmissivity decay rate at 22 and 100 ˚C is significantly lower than the transmissivity decay rate during the initial 3 to 15 hours of the experiment. However, at 200 °C, the decay of hydraulic transmissivity remains high throughout the experiment. The long-term decay of hydraulic transmissivity can be fitted with a power law model with more rapid reduction of hydraulic transmissivity at higher temperature. Periods of sliding on the fracture surface result in transient increases in the transmissivity, due to shear dilation, as is expected for Coulomb materials. These transients are superimposed on the long-term decay. When sliding ceases and a new hold period commences, there is a rapid reduction in transmissivity and return to the long-term rate of transmissivity decay. The rate of decay of the transmissivity transients is inversely proportional to temperature, in contrast to the long-term decay and the expected behavior for processes like subcritical crack growth and indentation creep. The higher decay rates that are observed during the initial 3-15 hours of the tests and following sliding, are associated with times that the porosity of the gouge is expected to be high. The difference in decay rates suggests that when the gouge is driven far from equilibrium by active shearing, densification may be dominated by a different mechanism from long-term compaction.

Conference Paper↗

USGS2022

Friction in clay-bearing faults increases with the ionic radius of interlayer cations

Smectite can dramatically reduce the strength of crustal faults and may cause creep on natural faults without great earthquakes; however, the frictional mechanism remains unexplained. Here, our shear experiments reveal systematic increase in shear strength with the increase of the ionic radius of interlayer cations among lithium-, sodium-, potassium-, rubidium-, and cesium-montmorillonites, a smectite commonly found in faults. Using density-functional-theory calculations, we find that relatively small sodium ions fit in the ditrigonal cavities on the montmorillonite surfaces, resulting in weakening of interlayer repulsion during sliding. On the other hand, relatively large potassium ions do not fit in the ditrigonal cavities, resulting in a larger resistance to sliding due to electrostatic repulsion between potassium ions. Calculated shear strength is consistent with our shear experiments by considering the partial dehydration of the frictional contact area. These results provide the basis for developing a quantitative model of smectite-bearing fault rheology.

Communications Earth & Environment↗

USGS2022

Electrical properties and anisotropy of schists and fault rocks from New Zealand’s Southern Alps under confining pressure

Magnetotelluric models spanning the Pacific–Australian Plate boundary in New Zealand’s South Island indicate a localized zone of low electrical resistivity that is spatially coincident with theductile mid-crustal part of the Alpine Fault Zone (AFZ). We explored the source of this anomaly bymeasuring the electrical properties of samples collected from surface outcrops approaching the AFZthat have accommodated a gradient of systematic strain and deformation conditions. We investigatedthe effects of tectonite fabric, fluid saturated pore/fracture networks and surface conductivity on the bulk electrical response and the anisotropy of resistivity measured under increasing confining pressures up to 200 MPa. We find that porosity and resistivity increase while porosity and the change in anisotropy of resistivity with confining pressure (δ (ρk/ρ⊥)/δ (peff)) decreases approaching theAFZ, indicating the electrical response is controlled by pore fluid conductivity and modified during progressive metamorphism. Conversely, Alpine mylonites exhibit relatively low resistivities at low porosities, and lower δ (ρk/ρ⊥)/δ (peff) than the schists. These findings indicate a transition in both the porosity distribution and electrical charge transport processes in rocks that have experienced progressive grain size reduction and mixing of phases during development of mylonitic fabrics due to creep shear strain within the AFZ.

Southern Alps↗

USGS2022

Impact of fluid-rock interaction on strength and hydraulic transmissivity rvolution in shear fractures under hydrothermal conditions

Reactivated shear fractures contribute to the creation of pervasive fracture networks in geothermal systems. The creation, reactivation, and sustainability of fracture networks depend on complex coupling among thermal, hydraulic, mechanical, and chemical (THMC) processes. However, most laboratory experiments focus either solely on how fluid transport properties evolve in stationary fractures at elevated temperatures or on how strength evolves in shear faults at room temperature. Studies that examine the combined evolution of these properties, especially under hydrothermal conditions, are limited. Laboratory restrengthening tests, referred to as slide-hold-slide (SHS), measure the difference between the steady sliding strength and the peak shear strength immediately following a hold period at constant stress or displacement. We present the results of five SHS experiments on Westerly granite bare surfaces. These experiments are combined with in-plane flow tests and conducted at an effective confining pressure of 20 MPa and temperatures of 23 and 200˚C. Consistent with previous studies, we observe a greater reduction in fluid transmissivity at 200˚C than at 23˚C and, generally, strength increasing with log(hold duration). However, there are deviations from the expected frictional healing behavior that likely relate to mineralogy and pore fluid chemical equilibrium. For example, under no-flow conditions where fluid is allowed to equilibrate with the rock, restrengthening rate of coefficient of friction is initially 0.02 per decade change in hold time (decade -1 ) for holds ≤ 5x10 4 s. Contrary to expected behavior, time-dependent reduction in frictional healing occurs at a rate of -0.03 decade -1 for holds longer than 5x10 4 s. For flow-through tests in which undersaturated fluid is continuously introduced to the fault surface, this reduction begins earlier, starting at 5x10 3 s, but at the same rate of -0.03 decade -1 . The earlier onset of reduction in frictional healing under the continuous flow conditions suggests that the process underlying the behavior is related to chemical equilibrium. The mechanisms behind the complex healing behavior are not fully understood but it is apparent that fault zone restrengthening under hydrothermal conditions is the result of multiple interacting processes.

Conference Paper↗

USGS2021

Weakening of peridotite sheared at hydrothermal conditions

We conducted triaxial friction tests at hydrothermal conditions (25°C–350°C) on gouges of peridotite and its principal mineral constituents olivine and orthopyroxene. Pore-fluid chemistry was varied by the use of peridotite, granite, or quartzite driving blocks (representing wall rock) housing the gouge layer. Samples sheared at slow rates initially strengthen to a peak value, and then weaken toward a residual strength. The transition is accompanied by a change from velocity-weakening to velocity-strengthening behavior marked by a series of small stress drops. The extent of weakening varies with the ultramafic mineralogy and with the chemical environment established by the driving block lithology. The strengths of olivine and olivine-rich peridotite gouges decrease substantially (to μ ∼ 0.25–0.30), and that of orthopyroxene to a lesser extent, at temperatures ≥200°C when sheared between crustal driving blocks. Less weakening is observed in the peridotite-block experiments; the minimum strength of the peridotite gouges ( μ ∼ 0.5) occurs at 250°C, the temperature at which olivine hydration rates are near their maximum in ultramafic rocks. The strength reductions in all experiments are attributed to solution-transfer (pressure solution) processes that come to predominate over cataclastic mechanisms during shear. The lower pH of fluids in contact with silica-saturated crustal rocks enhances the weakening of olivine-rich gouges. In these short-duration experiments, secondary phyllosilicate mineral growth was of a limited extent and varied with gouge and wall-rock mineralogy and with temperature. Over geologic time spans, however, the alteration assemblages will assume an increasingly important role in fault-zone behavior.

Geochemistry, Geophysics, Geosystems↗

USGS2021

A geology and geodesy based model of dynamic earthquake rupture on the Rodgers Creek‐Hayward‐Calaveras Fault System, California

The Hayward fault in California's San Francisco Bay area produces large earthquakes, with the last occurring in 1868. We examine how physics‐based dynamic rupture modeling can be used to numerically simulate large earthquakes on not only the Hayward fault, but also its connected companions to the north and south, the Rodgers Creek and Calaveras faults. Equipped with a wealth of images of this fault system, including those of its 3D geology and 3D geometry, in addition to inferences about its interseismic creep‐rate pattern and rock‐friction behavior, we use a finite‐element computer code to perform 3D dynamic earthquake rupture simulations. We find that the rock properties affect the locations and amount of slip produced in our simulated large earthquakes. Crucial factors that control rupture behavior in our modeling are the earthquake nucleation locations, the fault geometry, and the data that reveal where the fault system is creeping or locked. Our findings suggest that large Rodgers Creek‐Hayward‐Calaveras‐Northern Calaveras (RC‐H‐C‐NC) fault‐system earthquakes may result from dynamic rupture that starts in a locked part of the fault system, but is then stopped by the creeping parts, leading to high‐magnitude‐6 earthquakes; or, from dynamic rupture that starts in a locked part of the fault system, then cascades through some of the creeping parts, leading to magnitude‐7 earthquakes.

California↗

USGS2021

Evolution of fluid transmissivity and strength recovery of shear fractures under hydrothermal conditions

Geothermal systems rely on the presence of long-lived and high-volume, permeable fracture systems. The creation, reactivation, and sustainability of these systems depend on complex coupling among thermal, hydraulic, mechanical, and chemical (THMC) processes occurring in geothermal reservoirs. In part due to a paucity of experimental data, the evolution of fractures at geothermal conditions in response to THMC processes is poorly understood, particularly during the process of shear. We present preliminary results of triaxial slide-hold-slide experiments, with hold periods ranging in duration from 103 s to 106 s, to constrain rates and mechanisms of healing and sealing. Experiments were conducted on simulated fault gouge composed of Westerly granite and on bare surfaces of Westerly granite. The tests were run at temperatures of 22˚ and 200˚C with confining and average pore pressures of 30 MPa and 10 MPa, respectively. We used an axial displacement rate of 0.1 μm/s during sliding periods. Deionized water was continuously flowed along the simulated fracture so we could determine in-plane transmissivity during the tests. In gouge and bare surface experiments conducted at 200˚C, we observe significant decreases in transmissivity over the course of the experiments. For the hydrothermal gouge experiment we measured an order of magnitude net reduction in transmissivity from 1.73x10-18 to 0.17x10-18 m3, over the course of 220 hours while in the room temperature gouge experiment transmissivity only decreased by 0.35x10-18 m3 over the same amount of time. In the experiments, we observe an up to 16% recovery in transmissivity during sliding periods. At room temperature the friction data showed limited fault re-strengthening with time; healing rates are on the order of 0.1 MPa/decade. A similar healing rate was observed at 200˚C in the gouge but we observe an increase in the healing rate, to 0.75 MPa/decade, for a bare surface experiment at 200˚C. The differences in the healing rate of the gouge and bare surface experiments suggest that the generation of fine particles by grinding down of asperities on the bare surface promote quartz dissolution and reprecipitation at elevated temperatures. Further work is needed to test this possibility and provide better constraints on factors influencing the evolution of fluid transport properties and strength of shear fractures at geothermal conditions.

Conference Paper↗

USGS2021

Strength recovery and sealing under hydrothermal conditions

While there is significant evidence for healing in natural faults, geothermal reservoirs, and lab experiments, the thermal, hydraulic, mechanical, and chemical interactions that influence healing are poorly understood. We present preliminary results of triaxial slide-hold-slide experiments to constrain rates and mechanisms of healing. Experiments were conducted on gouge composed of Westerly granite and on bare surfaces of Westerly granite and Eureka quartzite. Tests were run at 22, 100, and 200˚C. In some experiments, we also determined the in-plane fluid transmissivity. In bare surface experiments we observe that restrengthening depends on both time and temperature. At 200˚C the simulated fractures restrengthen at a rate of ∆ µ /∆log(t hold ) = 0.009/decade while at 22˚C the healing rate is ~ 0.002/decade. In the gouge experiments restrengthening appears to be independent of temperature. This may be related to the heterogenous mineral composition and thickness of the gouge layer which could allow shearing to be accommodated in unhealed zones. In the experiments, an overall reduction in fluid transmissivity is observed but sliding periods are often associated with increases in the fluid transmissivity. The transmissivity reduction tends to be greater at 200˚C relative to room temperature. Our preliminary results suggest that multiple healing mechanisms are operating under hydrothermal conditions.

Conference Paper↗

USGS2020

Direct evidence for fluid pressure, dilatancy, and compaction affecting slip in isolated faults

Earthquake instability occurs as a result of strength loss during sliding on a fault. It has been known for over 50 years that fault compaction or dilatancy may cause significant weakening or strengthening by dramatically changing the fluid pressure trapped in faults. Despite this fundamental importance, we have no real understanding of the exact conditions that lead to compaction or dilation during nucleation or rupture. To date, no direct measurements of pore pressure changes during slip in hydraulically isolated faults have been reported. We show direct examples of fluid pressure variations during nucleation and rupture using a miniature pressure transducer embedded in an experimental fault. We demonstrate that fluids are not only significant in controlling fault behavior, but can provide the dominant mechanism controlling fault stability. The effect of fluid pressure changes can exceed frictional variations predicted by rate- and state-dependent friction laws, exerting fundamental controls on earthquake rupture initiation.

Geophysical Research Letters↗

USGS2020

Near-fault velocity spectra from laboratory failures and their relation to natural ground motion

We compared near-fault velocity spectra recorded during laboratory experiments to that of natural earthquakes. We fractured crystalline rock samples at room temperature and intermediate confining pressure (50 MPa). Subsequent slip events were generated on the fracture surfaces under higher confinement (300 MPa). Velocity spectra from rock fracture resemble the inverse frequency (1/ f ) decay of natural earthquake velocity. This spectrum can be attributed to fault creation via seismic fracturing over a wide range of spatial scales. In contrast, subsequent slips on the rough fracture surfaces are depleted in high frequency energy and falloff approximately as 1/ f 2 . The 1/ f 2 spectrum is more consistent with a slider-block model obeying static-kinetic friction than a natural earthquake. The depleted high frequency content precludes the rough fault experiments from being directly analogous to natural sources. The suppression of high frequencies may have resulted from two possible factors: (1) the presence of a well-developed shear zone and coseismic damping of the fault motion by dissipation within it or, in our favored interpretation, (2) a smaller amount of energy dissipated by shearing relative to the total energy release at elevated confining pressure. In context of the latter explanation, a unifying concept that applies to these experiments, earthquakes, ground motion, and models of complex radiated motion is that high frequency radiated energy is relatively enhanced when total energy release is nearly balanced within the source region by dissipative processes. This near-critical energy release condition can be accessed at low normal stress in laboratory experiments.

Journal of Geophysical Research↗

USGS2019

The susceptibility of Oklahoma’s basement to seismic reactivation

Recent widespread seismicity in Oklahoma is attributed to the reactivation of pre-existing, critically stressed and seismically unstable faults due to decades of wastewater injection. However, the structure and properties of the reactivated faults remain concealed by the sedimentary cover. Here, we explore the major ingredients needed to induce earthquakes in Oklahoma by characterizing basement faults in the field, in seismic surveys and via rock-mechanics experiments. Outcrop and satellite mapping reveal widespread fault and fracture systems with trends that display a marked similarity to the trends of recent earthquake lineaments. Our three-dimensional seismic analyses show steeply dipping basement-rooted faults that penetrate the overlying sedimentary sequences, representing pathways for wastewater migration. Experimental stability analysis indicates that Oklahoma’s basement rocks become seismically unstable at conditions relevant to the dominant hypocentral depths of the recent earthquakes. These analyses demonstrate that the geometry, structure and mechanical stability of Oklahoma’s basement make it critically susceptible to seismic reactivation.

Oklahoma↗