Search USGS⌕ Search

SEARCH · Search USGS

Results for “Physics of Fluids”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9Linked to original sources
USGS2016

Wave attenuation in the shallows of San Francisco Bay

Waves propagating over broad, gently-sloped shallows decrease in height due to frictional dissipation at the bed. We quantified wave-height evolution across 7 km of mudflat in San Pablo Bay (northern San Francisco Bay), an environment where tidal mixing prevents the formation of fluid mud. Wave height was measured along a cross shore transect (elevation range−2mto+0.45mMLLW) in winter 2011 and summer 2012. Wave height decreased more than 50% across the transect. The exponential decay coefficient λ was inversely related to depth squared (λ=6×10 −4 h −2 ). The physical roughness length scale k b , estimated from near-bed turbulence measurements, was 3.5×10 −3 m in winter and 1.1×10 −2 m in summer. Estimated wave friction factor f w determined from wave-height data suggests that bottom friction dominates dissipation at high Re w but not at low Re w . Predictions of near-shore wave height based on offshore wave height and a rough formulation for f w were quite accurate, with errors about half as great as those based on the smooth formulation for f w . Researchers often assume that the wave boundary layer is smooth for settings with fine-grained sediments. At this site, use of a smooth f w results in an underestimate of wave shear stress by a factor of 2 for typical waves and as much as 5 for more energetic waves. It also inadequately captures the effectiveness of the mudflats in protecting the shoreline through wave attenuation.

California↗
USGS2007

BAT3 Analyzer: Real-time data display and interpretation software for the multifunction bedrock-aquifer transportable testing tool (BAT3)

The BAT3 Analyzer provides real-time display and interpretation of fluid pressure responses and flow rates measured during geochemical sampling, hydraulic testing, or tracer testing conducted with the Multifunction Bedrock-Aquifer Transportable Testing Tool (BAT3) (Shapiro, 2007). Real-time display of the data collected with the Multifunction BAT3 allows the user to ensure that the downhole apparatus is operating properly, and that test procedures can be modified to correct for unanticipated hydraulic responses during testing. The BAT3 Analyzer can apply calibrations to the pressure transducer and flow meter data to display physically meaningful values. Plots of the time-varying data can be formatted for a specified time interval, and either saved to files, or printed. Libraries of calibrations for the pressure transducers and flow meters can be created, updated and reloaded to facilitate the rapid set up of the software to display data collected during testing with the Multifunction BAT3. The BAT3 Analyzer also has the functionality to estimate calibrations for pressure transducers and flow meters using data collected with the Multifunction BAT3 in conjunction with corroborating check measurements. During testing with the Multifunction BAT3, and also after testing has been completed, hydraulic properties of the test interval can be estimated by comparing fluid pressure responses with model results; a variety of hydrogeologic conceptual models of the formation are available for interpreting fluid-withdrawal, fluid-injection, and slug tests.

Open-File Report↗
USGS2011

Direction of unsaturated flow in a homogeneous and isotropic hillslope

The distribution of soil moisture in a homogeneous and isotropic hillslope is a transient, variably saturated physical process controlled by rainfall characteristics, hillslope geometry, and the hydrological properties of the hillslope materials. The major driving mechanisms for moisture movement are gravity and gradients in matric potential. The latter is solely controlled by gradients of moisture content. In a homogeneous and isotropic saturated hillslope, absent a gradient in moisture content and under the driving force of gravity with a constant pressure boundary at the slope surface, flow is always in the lateral downslope direction, under either transient or steady state conditions. However, under variably saturated conditions, both gravity and moisture content gradients drive fluid motion, leading to complex flow patterns. In general, the flow field near the ground surface is variably saturated and transient, and the direction of flow could be laterally downslope, laterally upslope, or vertically downward. Previous work has suggested that prevailing rainfall conditions are sufficient to completely control these flow regimes. This work, however, shows that under time-varying rainfall conditions, vertical, downslope, and upslope lateral flow can concurrently occur at different depths and locations within the hillslope. More importantly, we show that the state of wetting or drying in a hillslope defines the temporal and spatial regimes of flow and when and where laterally downslope and/or laterally upslope flow occurs.

Water Resources Research↗
USGS1971

A paleohydrologic model for mineralization of the White Pine copper deposit, northern Michigan

Pertinent physical properties of the upper Keweenawan rocks can be measured or inferred within a sufficiently narrow range to make the quantitative evaluation of various paleohydrologic models for the origin of the White Pine copper deposit feasible. The approach is illustrated here by calculations for models that involve lateral migration of fluids through the subjacent Copper Harbor Conglomerate to the site of the deposit and stripping of copper from these solutions where they percolated upward through the Nonesuch Shale. The calculations reveal limitations to theories of origin that would not be evident from purely qualitative consideration; some of these limitations could be useful to exploration. For example, if the water was yielded by compaction of the Copper Harbor Conglomerate and contained 50 ppm Cu, there must have been significant convergence of solution paths toward White Pine. Surface water entering the Copper Harbor Conglomerate on the north limb of the Lake Superior syncline is an adequate source if it could be shown that the point of entry was once significantly higher in altitude than the water table at White Pine; this model implies a major copper deposit at great depth north of the axis of the syncline. Ground water entering the Nonesuch Shale up dip from White Pine is not a possible source of mineralizing solutions.

Michigan↗
USGS2014

Resolving terrestrial ecosystem processes along a subgrid topographic gradient for an earth-system model

Soil moisture is a crucial control on surface water and energy fluxes, vegetation, and soil carbon cycling. Earth-system models (ESMs) generally represent an areal-average soil-moisture state in gridcells at scales of 50–200 km and as a result are not able to capture the nonlinear effects of topographically-controlled subgrid heterogeneity in soil moisture, in particular where wetlands are present. We addressed this deficiency by building a subgrid representation of hillslope-scale topographic gradients, TiHy (Tiled-hillslope Hydrology), into the Geophysical Fluid Dynamics Laboratory (GFDL) land model (LM3). LM3-TiHy models one or more representative hillslope geometries for each gridcell by discretizing them into land model tiles hydrologically coupled along an upland-to-lowland gradient. Each tile has its own surface fluxes, vegetation, and vertically-resolved state variables for soil physics and biogeochemistry. LM3-TiHy simulates a gradient in soil moisture and water-table depth between uplands and lowlands in each gridcell. Three hillslope hydrological regimes appear in non-permafrost regions in the model: wet and poorly-drained, wet and well-drained, and dry; with large, small, and zero wetland area predicted, respectively. Compared to the untiled LM3 in stand-alone experiments, LM3-TiHy simulates similar surface energy and water fluxes in the gridcell-mean. However, in marginally wet regions around the globe, LM3-TiHy simulates shallow groundwater in lowlands, leading to higher evapotranspiration, lower surface temperature, and higher leaf area compared to uplands in the same gridcells. Moreover, more than four-fold larger soil carbon concentrations are simulated globally in lowlands as compared with uplands. We compared water-table depths to those simulated by a recent global model-observational synthesis, and we compared wetland and inundated areas diagnosed from the model to observational datasets. The comparisons demonstrate that LM3-TiHy has the capability to represent some of the controls of these hydrological variables, but also that improvement in parameterization and input datasets are needed for more realistic simulations. We found large sensitivity in model-diagnosed wetland and inundated area to the depth of conductive soil and the parameterization of macroporosity. With improved parameterization and inclusion of peatland biogeochemical processes, the model could provide a new approach to investigating the vulnerability of Boreal peatland carbon to climate change in ESMs.

Hydrology and Earth System Sciences↗
USGS1990

Textural development of clayey and quartzofeldspathic fault gouges relative to their sliding behavior

Many of the secondary fault structures developed during triaxial friction experiments have been generally correlated with the structures of natural fault zones. Therefore, any physical differences that can be found between laboratory samples that slide stably and those that show stick-slip motion may help to identify the cause of earthquakes. We have examined petrographically the run products of many triaxial friction experiments using clayey and quartzofeldspathic gouges, which comprise the principal types of natural fault gouge material. The examined samples were tested under a wide range of temperature, confining and fluid pressure, and velocity conditions. The clayey and quartzofeldspathic gouges show some textural differences, owing to their different mineral contents and grain sizes and shapes. In the clayey gouges, for example, a clay mineral fabric and kink band sets are commonly developed, whereas in the quartzofeldspathic gouges fracturing and crushing of the predominately quartz and feldspar grains are important processes. For both types of gouge, however, and whatever the pressure-temperature-velocity conditions of the experiments, the transition from stable sliding to stick-slip motion is correlated with: (i) a change from pervasive deformation of the gouge layer to localized slip in subsidiary shears; and (ii) an increase in the angle betweem the shears that crosscut the gouge layer (Riedel shears) and ones that form along the gouge-rock cylinder boundaries (boundary shears). This suggests that the localization of shear within a fault zone combined with relatively high Riedel-shear angles are somehow connected with earthquakes. Secondary fracture sets similar to Riedel shears have been identified at various scales in major strike-slip faults such as the San Andreas of the western United States (Wallace, 1973) and the Luhuo and Fuyun earthquake faults of China (Deng and Zhang, 1984; Deng et al. , 1986). The San Andreas also contains locked and creeping sections that correspond to the stick-slip and stably sliding experimental samples, respectively. We plan to study the physical structure of the San Andreas fault, to see if the experimentally observed differences related to sliding behavior can also be distinguished in the field.

Physics and Chemistry of the Earth↗
USGS2014

Simultaneous Gaussian and exponential inversion for improved analysis of shales by NMR relaxometry

Nuclear magnetic resonance (NMR) relaxometry is commonly used to provide lithology-independent porosity and pore-size estimates for petroleum resource evaluation based on fluid-phase signals. However in shales, substantial hydrogen content is associated with solid and fluid signals and both may be detected. Depending on the motional regime, the signal from the solids may be best described using either exponential or Gaussian decay functions. When the inverse Laplace transform, the standard method for analysis of NMR relaxometry results, is applied to data containing Gaussian decays, this can lead to physically unrealistic responses such as signal or porosity overcall and relaxation times that are too short to be determined using the applied instrument settings. We apply a new simultaneous Gaussian-Exponential (SGE) inversion method to simulated data and measured results obtained on a variety of oil shale samples. The SGE inversion produces more physically realistic results than the inverse Laplace transform and displays more consistent relaxation behavior at high magnetic field strengths. Residuals for the SGE inversion are consistently lower than for the inverse Laplace method and signal overcall at short T2 times is mitigated. Beyond geological samples, the method can also be applied in other fields where the sample relaxation consists of both Gaussian and exponential decays, for example in material, medical and food sciences.

Journal of Magnetic Resonance↗
USGS1997

Review of electric and magnetic fields accompanying seismic and volcanic activity

New observations of magnetic, electric and electromagnetic field variations, possibly related to recent volcanic and seismic events, have been obtained on Mt. Unzen in Japan, Reunion Island in Indian Ocean, the Long Valley volcanic caldera in California, and for faults in China and Russia, California and several other locations. For volcanic events, contributions from different physical processes can be identified during the various eruption stages. Slow processes (weeks to months) include near-surface thermal demagnetization effects, piezomagnetic effects, and effects from rotation/displacement of magnetized material. Rapid processes (seconds to days) include piezomagnetic effects from instantaneous stress redistribution with explosive eruptions and electrokinetic effects from rupture of high pressure fluid compartments commonly encountered in volcanic regions. For seismic events, the observed coseismic offsets are instantaneous, provided care has been taken to ensure sensors are insensitive to seismic shaking and are in regions of low magnetic field gradient. Simple piezomagnetic dislocation models based on geodetically and seismically determined fault parameters generally match the observed signals in size and sign. Electrokinetic effects resulting from rupture of fluid filled compartments at hydrostatic to lithostatic pore pressures can generate transient signals in the frequency band 100 Hz to 0.01 Hz. However, large-scale fluid driven processes are not evident in near-field measurements in the epicentral region minutes to weeks before large earthquakes. The subset of ionospheric disturbances generated by trapped atmospheric pressure waves (also termed gravity waves and/or acoustic waves, traveling ionospheric disturbances or TID's) that are excited by earthquakes and volcanic eruptions are common and propagate to great distances. These are known and expected consequences of earthquakes, volcanic explosions (and other atmospheric disturbances), that must be identified and their effects removed from VLF/ULF electromagnetic field records before associating new observations of ionospheric disturbances with earthquake activity.

Surveys in Geophysics↗
USGS2007

Dynamic triggering

A growing body of evidence demonstrates that dynamic stresses propagating as seismic waves from large earthquakes are capable of triggering additional earthquakes ranging from aftershocks in the near-field (within one or two source dimensions of the mainshock epicenter) to remotely triggered earthquakes at distances exceeding 10 000 km. Most of the triggered earthquakes are small (generally M ≤ 3) except within the near field, where dynamic stresses may trigger slip on subadjacent fault segments leading to complex rupture comprising of several large earthquakes of comparable magnitude. Crustal surface waves with periods of 15–30 s and peak dynamic stresses greater than ∼0.01 MPa seem to be most efficient in triggering remote seismicity. Current models for dynamic triggering fall under two broad groups: one appealing to Coulomb failure with various friction laws, and the other appealing to the activation of crustal fluids either hydrous or magmatic. No single model appears capable of accounting for the wide variation observed in the nature of triggered activity. Spatial sampling of dynamic triggering on a global scale is still woefully inadequate because of the limited distribution of adequate seismic networks. From the limited data currently available, it appears that extensional stress regimes hosting geothermal and volcanic activity are more susceptible to remote dynamic triggering than compressional stress regimes, although remote triggering is not limited to extensional regimes. Instances of remote triggering in the few areas with continuous, high-resolution deformation instrumentation (all volcanic or geothermal areas) include distinctive deformation transients, suggesting that the locally triggered seismicity in these areas may be a secondary response to a more fundamental aseismic process that likely involves some form of fluid transport or phase change. Recent evidence for triggering by solid Earth tides and ocean loading in convergent plate margins provides a low-frequency, low-amplitude reference point for the spectrum of stresses capable of dynamic triggering. Remaining challenges include establishing better sampling of the distribution of triggered seismicity and better constraints on physical models for the triggering process.

Book chapter↗
USGS2014

Appendix C—Report on research in the field of ground water being conducted by oil companies

In view of the shortness of time since the appointment of the writer to the Committee on Ground Water this report is confined to the technology and problems in the Gulf Coast Oil Province. Of course, many of the methods and practices would apply to most parts of the country however, some would differ materially from one region to another. The writer wishes to acknowledge the suggestions and comments by F. H. LAHEE and PAUL WEAVER. Having been stationed in Houston, Texas, in the heart of the Gulf Coast Area for four and a half years, the writer has had an opportunity to view the great similarity of the problems confronting the petroleum geologist and engineer and the ground‐water hydrologist. Both groups deal with the accumulation, movement, and withdrawal of fluid from underground strata, yet each group is content to study its own literature and use its own terminology without much concern for the other. The petroleum and ground‐water engineer, independently of one another, have developed mathematical formulas for the determination of permeability from field‐data. These formulas use the same basic principles of physics and the initial papers on the subject by both groups were published within two years of one another. Because of the similarity in the technology and problems of the petroleum engineer and the ground‐water hydrologist there is a definite need for closer cooperation. Some of the problems are so closely related that their solution rests in cooperative studies.

Eos, Transactions, American Geophysical Union↗
USGS2016

Comment on “The reduction of friction in long-runout landslides as an emergent phenomenon” by Brandon C. Johnson et al.

Results from a highly idealized, 2-D computational model indicate that dynamic normal-stress rarefactions might cause friction reduction in long-runout landslides, but the physical relevance of the idealized dynamics has not been confirmed by experimental tests. More importantly, the model results provide no evidence that refutes alternative hypotheses about friction reduction mechanisms. One alternative hypothesis, which is strongly supported by field evidence, experimental data, and the predictions of a well-constrained computational model, involves development of high pore fluid pressures in deforming landslide material or overridden bed material. However, no scientific basis exists for concluding that a universal mechanism is responsible for friction reduction in all long-runout landslides.

Journal of Geophysical Research F: Earth Surface↗
USGS2019

Valid debris-flow models must avoid hot starts

Debris-flow experiments and models commonly use “hot-start” initial conditions in which downslope motion begins when a large force imbalance is abruptly imposed. By contrast, initiation of natural debris flows almost invariably results from small perturbations of static force balances that apply to debris masses poised in steep channels or on steep slopes. Models that neglect these static balances may violate physical law. Here we assess how the effects of hot starts are manifested in physical experiments, analytical dam-break models, and numerical models in which frictional resistance is too small to satisfy static force balances in debris-flow source areas. We then outline a numerical modeling framework that avoids use of hot starts. In this framework an initial static force balance is gradually perturbed by increasing pore-fluid pressure that may trigger the onset of debris motion. Subsequent increases in pore-fluid pressure, driven by debris motion, may then reduce the debris frictional strength, leading to high flow mobility.

Conference Paper↗
USGS1986

Shallow subsurface temperature surveys in the Basin and Range province, U.S.A.-I. Review and evaluation

Temperature surveys at depths of 1–2 m have had varying success in geothermal exploration in the Basin and Range province. The most successful surveys have identified patterns of near-surface thermal-fluid flow within areas of less than 2 km 2 . Results have been less consistent in larger areas where zones of hydrothermal upflow are less well known, nongeothermal perturbing factors are significant and lateral variations in shallow subsurface temperature are small. Nongeothermal perturbations can be minimized by use of mean annual temperatures instead of synoptic temperatures, by physically based simulation of ground temperatures or by statistical modeling.

Geothermics↗
USGS2021

Development of soil radiocarbon profiles in a reactive transport framework

Today, there is a greater appreciation for the importance of the physical protection of carbon (C) through interactions with mineral surfaces, isolation from microbes, and the important role of transport in shaping soil properties and controlling moisture limitations on decomposition. As our paradigm for soil organic carbon (SOC) preservation changes, so too should our representation of the underlying processes in soil models. Reactive transport models (RTMs) provide a framework capable of assessing the interactive influence of soil chemistry and transport processes on the accumulation and turnover of SOC. In this study, we present new developments in the isotopically enabled RTM “CrunchTope,” which is capable of explicitly tracking the three isotopes of carbon ( 12 C, 13 C, and 14 C) and their fractionation between multiple coexisting and interacting solid, liquid and gas phases. This modeling framework opens the door to new applications of depth-resolved RTMs models in application to SOC and deeper subsurface carbon reservoirs. Here, we demonstrate SOC accumulation and radiocarbon aging for long-timescale models of soil development in CrunchTope. Our goal is to assess advantages and limitations of such an approach and to identify the type and complexity of reaction networks that are required to adequately apply this model to SOC dynamics. We assess the behavior of this model relative to a high-resolution dataset of SOC content, stable isotope composition, and radiocarbon ages as well as physical and hydrologic data measured from a chronosequence of soils located near Santa Cruz, California. Starting from a previously published model using a simplified reaction network with a single class of carbon, we sequentially incorporate multiple C reservoirs subject to both reactivity and transport pathways. Our results indicate that multiple SOC pools with different mean ages of C do not inherently emerge as a result of including reactions which are conventionally expected to provide a diversity of transit times, i.e., sorption and complexation of SOC on mineral surfaces. Instead, transit times emerge as a result of the timescales of the reactions represented in the reaction network. For mineral associated C, the RTM framework imposes dynamic equilibrium with the fluid phase dissolved organic C, such that no distinction in radiocarbon ages is achieved between these pools. Aged C can be produced by including a solid-phase C reservoir, with a rate-limited solubilization coefficient. Aging of SOC in this way is more akin to selective preservation than to mineral protection and, while such a mechanism may be at play in many soils, mineral protection is thought to be at least as important. As such, our results indicate that additional parameterization is required to reproduce the heterogeneity of carbon transit times that result from organo-mineral interactions. These efforts show the promise of a modeling approach where the varied transit time of soil C emerges from the dynamic physical and hydrologic properties of the model rather than from the a priori assignment of operationally defined pools.

Geochimica et Cosmochimica Acta↗
USGS2020

Post-1978 tumescence at Long Valley Caldera, California: A geophysical perspective

Long Valley Caldera has been restless since at least 1978. Prominent symptoms of this unrest include earthquake swarms and tumescence (inflation) centered on the resurgent dome. Over the years, interpretations of physical processes underlying this unrest have varied considerably. Results from a collection of geophysical studies infer the presence and/or active intrusion of magma in the crust. Geologic evidence, however, does not support recent magmatic activity in the caldera, leading to an interpretation that the caldera volcanic system is moribund, and the current unrest is a result of second boiling (aqueous fluids released during crystallization of the rhyolitic magma that produced the Bishop Tuff). Here, we examine the collective constraints provided by geophysical studies over the past four decades. Although the current geophysical evidence does not conclusively discriminate between unrest driven by recent crustal magmatic intrusion versus second boiling, it does provide evidence that a large volume of partial melt persists within the mid- and lower crust. This implies that the Long Valley Caldera system as a whole is long-lived, and the magma reservoir remains at least partially molten. In the shallow crust, the possibility of small pockets of magma remains. In aggregate, the data suggest commonality with other large, long-lived silicic caldera systems, such as Yellowstone and Campi Flegrei. Although the possibility of eruption within Long Valley Caldera remains unlikely, the geophysical evidence argues that Long Valley has an active magmatic system at depth and we must retain the possibility of eruptive hazards.

California↗
USGS2013

Iron mineralogy and bioaccessibility of dust generated from soils as determined by reflectance spectroscopy and magnetic and chemical properties--Nellis Dunes recreational area, Nevada

Atmospheric mineral dust exerts many important effects on the Earth system, such as atmospheric temperatures, marine productivity, and melting of snow and ice. Mineral dust also can have detrimental effects on human health through respiration of very small particles and the leaching of metals in various organs. These effects can be better understood through characterization of the physical and chemical properties of dust, including certain iron oxide minerals, for their extraordinary radiative properties and possible effects on lung inflammation. Studies of dust from the Nellis Dunes recreation area near Las Vegas, Nevada, focus on characteristics of radiative properties (capacity of dust to absorb solar radiation), iron oxide mineral type and size, chemistry, and bioaccessibility of metals in fluids that simulate human gastric, lung, and phagolysosomal fluids. In samples of dust from the Nellis Dunes recreation area with median grain sizes of 2.4, 3.1, and 4.3 micrometers, the ferric oxide minerals goethite and hematite, at least some of it nanosized, were identified. In one sample, in vitro bioaccessibility experiments revealed high bioaccessibility of arsenic in all three biofluids and higher leachate concentration and bioaccessibility for copper, uranium, and vanadium in the simulated lung fluid than in the phagolysosomal fluid. The combination of methods used here to characterize mineral dust at the Nellis Dunes recreation area can be applied to global dust and broad issues of public health.

Nevada↗
USGS2006

Quantity-activity relationship of denitrifying bacteria and environmental scaling in streams of a forested watershed

The spatial variability of subreach denitrification rates in streams was evaluated with respect to controlling environmental conditions, molecular examination of denitrifying bacteria, and dimensional analysis. Denitrification activities ranged from 0 and 800 ng-N gsed-1 d-1 with large variations observed within short distances (<50 m) along stream reaches. A log-normal probability distribution described the range in denitrification activities and was used to define low (16% of the probability distributibn), medium (68%), and high (16%) denitrification potential groups. Denitrifying bacteria were quantified using a competitive polymerase chain reaction (cPCR) technique that amplified the nirK gene that encodes for nitrite reductase. Results showed a range of nirK quantities from 103 to 107 gene-copy-number gsed.-1 A nonparametric statistical test showed no significant difference in nirK quantifies among stream reaches, but revealed that samples with a high denitrification potential had significantly higher nirK quantities. Denitrification activity was positively correlated with nirK quantities with scatter in the data that can be attributed to varying environmental conditions along stream reaches. Dimensional analysis was used to evaluate denitrification activities according to environmental variables that describe fluid-flow properties, nitrate and organic material quantities, and dissolved oxygen flux. Buckingham's pi theorem was used to generate dimensionless groupings and field data were used to determine scaling parameters. The resulting expressions between dimensionless NO3- flux and dimensionless groupings of environmental variables showed consistent scaling, which indicates that the subreach variability in denitrification rates can be predicted by the controlling physical, chemical, and microbiological conditions. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2017

Poroelastic properties of the Arbuckle Group in Oklahoma derived from well fluid level response to the 3 September 2016 M w 5.8 Pawnee and 7 November 2016 M w 5.0 Cushing earthquakes

The Arbuckle Group (Arbuckle) is a basal sedimentary unit that is the primary target for saltwater disposal in Oklahoma. Thus, the reservoir characteristics of the Arbuckle, including how the poroelastic properties change laterally and over time are of significant interest. We report observations of fluid level changes in two monitoring wells in response to the 3 September 2016 M w 5.8 Pawnee and the 7 November 2016 M w 5.0 Cushing earthquakes. We investigate the relationship between static strain resulting from these events and the fluid level changes observed in the wells. We model the fluid level response by estimating static strains from a set of earthquake source parameters and spatiotemporal poroelastic properties of the Arbuckle in the neighborhood of the monitoring wells. Results suggest that both the direction of the observed fluid level step and the amplitude can be predicted from the computed volumetric strain change and a reasonable set of poroelastic parameters. Modeling results indicate that poroelastic parameters differ at the time of the Pawnee and Cushing earthquakes, with a moderately higher Skempton’s coefficient required to fit the response to the Cushing earthquake. This may indicate that dynamic shaking resulted in physical alteration of the Arbuckle at distances up to ∼50 km from the Pawnee earthquake.

Oklahoma↗