Geology ReportsSearch

SEARCH · Geology Reports

Results for “Natural Hazards and Earth Systems Sciences (NHESS)”

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.

1,668 records · Page 40Linked to original sources

Long‐period ground motions from dynamic rupture simulations of large earthquakes on the creeping Hayward–Calaveras–Rodgers Creek fault system

he Hayward, Calaveras, and Rodgers Creek faults in the San Francisco Bay region of California have a high probability of producing a large earthquake in the next decades. Although these faults creep, the creep is insufficient to keep up with their relatively rapid slip rates on their deepest sections, so they have been storing tectonic strain since their last large earthquakes, with the Hayward’s and Rodgers Creek’s more than 150 yr ago. We do not know what the next large Hayward–Calaveras–Rodgers Creek earthquakes will look like or how strongly they will shake the San Francisco Bay region. Harris et al. (2021) used the 3D dynamic (spontaneous) rupture method to simulate large earthquakes on these creeping faults. In this article, we examine the resulting simulated long‐period ( T > 1 s) ground shaking from 0 to 50 km distance, for earthquakes nucleating on the Hayward fault and earthquakes nucleating on the Rodgers Creek fault. We compare these simulated long‐period ground motions with the Boore et al. (2014) well‐established empirically based ground‐motion model suitable for the slowest material velocity in our 3D velocity structure. We find that the simulated long‐period ground motions from the creeping‐fault earthquake scenarios produce a reasonable agreement with the empirical expectations if frictional cohesion is included only where it is appropriate.

California

Are field observations of surface rupture useful? An example from the 2023 Mw 7.8 Pazarcık, Turkey (Türkiye), earthquake

Field investigations have long been an important component of the scientific response to surface‐faulting earthquakes. However, in light of advances in remote data and models, the question arises whether field‐based observations of surface rupture remain useful for understanding rupture processes and seismic hazards. We approach this question using a field‐based study of the central 2023 M w 7.8 Pazarcık, Turkey (Türkiye), earthquake rupture, at the intersection of the east Anatolian fault (EAF) and Narlı fault. Our field observations include the surface rupture expression and extent of the central EAF and northernmost Narlı fault in generally forested and steep terrain and 68 measurements of left‐lateral surface displacement. These data improve our understanding of the Pazarcık rupture complexity, resolve the surface geometry of the Narlı fault–EAF connection, and exhibit a clear (>2 m) change in surface displacement across this intersection zone that confirms remote‐based coseismic slip models. Our study shows that focusing field efforts in areas of obscured or low‐resolution remote data can yield essential data for refining rupture extent, documenting perishable on‐fault displacement, and improving postearthquake situational awareness. A comparison of similarly large‐magnitude continental surface‐rupturing earthquakes indicates that displacement uncertainties relate to a complex set of factors, including measurement methods, rupture complexity, and displacement magnitude. Our study validates the need for postearthquake field observations, which, when driven by clear motivating questions and knowledge of methodological strengths and limitations, provide high‐resolution rupture data that complement remote‐based models.

Pazarcık

Assessing legacy nitrogen in groundwater using numerical models of the Long Island aquifer system, New York

Nitrogen transported along groundwater flow paths in coastal aquifers can contribute substantially to nitrogen loading into surface water receptors, particularly in hydrologic systems dominated by groundwater discharge. Nitrogen entrained in the aquifer is a function of land use and associated nitrogen sources at the time of groundwater recharge, which may differ considerably from present-day sources. Legacy nitrogen can result in substantial discrepancies between observed present-day nitrogen loading to surface water receptors and loading estimated from present-day sources. Additionally, legacy nitrogen can continue to discharge into surface waters after nitrogen mitigation actions have been undertaken. Here, we use a numerical modeling framework to compare three methods of estimating time-varying historical nitrogen loads to four water bodies (receptors) on eastern Long Island, New York. The methods span a range of data requirements and process complexity, from instantaneous receptor loads calculated from steady-state groundwater contributing areas, to transient loads estimated by explicitly simulating legacy groundwater nitrogen transport over a century with large changes in nitrogen sources and hydrologic conditions. The effects of legacy nitrogen on estimated receptor loads varied temporally and spatially within the study area. Depending on antecedent nitrogen inputs and hydrologic conditions, historical annual nitrogen loads estimated from transient simulations accounting for legacy nitrogen can be quite similar (<10% difference) or substantially different (±100%) from those estimated from simpler instantaneous methods. Continued input of present-day nitrogen sources using methods that account for legacy nitrogen results in asymptotic increases in receptor nitrogen loads over time, indicating that simulated present-day receptor nitrogen loads are not in equilibrium with present-day inputs. For these receptors in disequilibrium, models simulating transient groundwater nitrogen transport could be used to account for legacy nitrogen lag times to help resource managers evaluate the potential effectiveness of proposed nitrogen mitigation actions.

EarthArXiv

Borehole geophysical time-series logging to monitor passive ISCO treatment of residual chlorinated-ethenes in a confining bed, NAS Pensacola, Florida

In-situ chemical oxidation (ISCO) is a common method to remediate chlorinated ethene contaminants in groundwater. Monitoring the effectiveness of ISCO can be hindered because of insufficient observations to assess oxidant delivery. Advantageously, potassium permanganate, one type of oxidant, provides the opportunity to use its strong electrical signal as a surrogate to track oxidant delivery using time-series borehole geophysical methods, like electromagnetic (EM) induction logging. Here we report a passive ISCO (P-ISCO) experiment, using potassium permanganate cylinders emplaced in boreholes, at a chlorinated ethene contamination site, Naval Air Station Pensacola, Florida. The contaminants are found primarily at the base of a shallow sandy aquifer in contact with an underlying silty-clay confining bed. We used results of the time-series borehole logging collected between 2017 and 2022 in 4 monitoring wells to track oxidant delivery. The EM-induction logs from the monitoring wells showed an increase in EM response primarily along the contact, likely from pooling of the oxidant, during P-ISCO treatment in 2021. Interestingly, concurrent natural gamma-ray (NGR) logging showed a decrease in NGR response at 3 of the 4 wells possibly from the formation of manganese precipitates coating sediments. The coupling of time-series logging and well-chemistry data allowed for an improved assessment of passive ISCO treatment effectiveness.

Florida

Using the D-Claw software package to model lahars in the Middle Fork Nooksack River drainage and beyond, Mount Baker, Washington

Lahars, or volcanic mudflows, are the most hazardous eruption-related phenomena that will affect communities living along rivers that originate on Mount Baker. In the past 15,000 years, the largest lahars from Mount Baker have affected the Middle Fork Nooksack River drainage and beyond. Here we use the physics-based D-Claw software package to model nine lahar scenarios that are initiated as water-saturated landslides between Sherman Crater and the Roman Wall on the Mount Baker edifice and flow down the Middle Fork Nooksack River. The scenarios range in volume from 1 to 260 million cubic meters and have an initial hydraulic permeability from 10 −12 to 10 −10 meters squared. Model output includes data such as flow depth, velocity, runout distance, area inundated, arrival time, and sediment concentration as well as information that allows scientists to calculate other important hydrologic characteristics such as lahar discharge. These data are important to officials who have the responsibility to plan for, or take mitigation measures against, future Mount Baker lahars. To check the validity of the D-Claw results, we compare the scenarios to known geologic information. We also compare D-Claw results with empirical models that have been used in the past to determine potential inundation areas, runout distances, and arrival times. These comparisons highlight similarities and differences between empirical and physics-based models. We also present D-Claw scenario-based animations to help scientists, officials, and lay people alike to visualize how future lahars could affect communities.

Washington

Using public participatory geographic information systems (PPGIS) to explore uses and values for Mojave Trails National Monument, California

Many people ascribe a variety of values to public lands and waters, but some values are more difficult to assess and quantify than others. Public participatory geographic information systems (PPGIS) are tools that have been used to help quantify and map the public’s diverse values for a landscape. This work describes the first known Office of Management and Budget–approved use of PPGIS by a Department of the Interior bureau. The U.S. Geological Survey developed an internet-based application to aid in gathering PPGIS data, called Values Mapping for Planning in Regional Ecosystems (VaMPIRE). Further, this work describes the first pilot of the VaMPIRE application in coordination with the Bureau of Land Management to collect spatial data and other survey data regarding the public’s uses of and values for locations within Mojave Trails National Monument. We emailed the link to the VaMPIRE application to an interested party email list in 2024 with 207 valid emails and received 74 responses; we also received 47 responses from members of an off-roading social media group. Of the list of 16 value options, recreation was the most popular value for the monument, followed by wilderness and inspirational. Over 1,000 points were placed throughout the monument, indicating locations people use or value, with the locations spread throughout the entire monument. Additionally, most survey respondents stated their ability to receive benefits in locations they mapped would not change in response to a hypothetical scenario related to recreational facility development. This report describes exploratory results from the first use of the VaMPIRE tool in Mojave Trails National Monument and includes reflections on how the process went and considerations for future use of VaMPIRE.

California

High‐resolution surface deformation and slip distribution observations for the 2023 Kahramanmaraş, Türkiye, earthquake sequence help constrain the rupture process

Splay, or branch, faults are a common geometric feature of earthquake surface ruptures and may provide constraints on the rupture behavior of an earthquake. The 2023 M w 7.8 Pazarcık and M w 7.5 Elbistan, Türkiye, earthquakes are examples of ruptures with multiple small splays, and the Pazarcık earthquake nucleated on a splay fault, the Narlı fault, before rupturing bilaterally on the East Anatolian fault (EAF). Here, we present 3‐m‐resolution surface displacement from subpixel correlation of Planet Dove optical images for the entirety of both ruptures with corresponding surface slip distributions. For a 30‐km‐long study region spanning the Narlı‐EAF intersection, we compare surface slip derived from five data sets with different resolutions (on‐the‐ground, WorldView, Planet Dove, Sentinel‐2, and Sentinel‐1) to elucidate complementary information. In addition, we integrate information from the surface expression of faulting with published dynamic rupture simulations and rupture process studies to constrain a rupture evolution for the Pazarcık earthquake that is consistent across data sets. This work highlights the complementary nature of disparate surface slip data sets and the role that high‐resolution surface displacement information, including from fault splays, can play in constraining nonunique rupture models and refining understanding of the earthquake rupture process.

Bulletin of the Seismological Society of America

Archean to Mesozoic–Cenozoic seismic crustal structure: Implications for geological and biological evolution

We use >4500 measurements of crustal structure to investigate the seismic structure of continental crust, Archean to Mesozoic–Cenozoic. The mean crustal thickness of continents, including their margins, is 36.5 km. We find that Archean, Paleoproterozoic, and Mesoproterozoic crust have similar mean crustal thickness (41 km), seismic velocities, Vp/Vs ratio, and density. Crusts of these ages span 3 Ga, from 4.0 to 1.0 Ga, and their similarity of physical properties (thickness, density, and seismic velocities) suggests that the process of crustal formation may have been similar during this time period, which covers 66% of Earth's history. Neoproterozoic crustal properties differ significantly from Mesoproterozoic and older crusts. The mean thickness of Neoproterozoic crust is 32 km, some 9 km thinner than the mean thickness of Archean, Paleoproterozoic, and Mesoproterozoic crust. A lithospheric root with a thickness of 150–200 km underlies Archean, Paleoproterozoic, and Mesoproterozoic crusts, and this root resists lithospheric rifting and crustal extension and thinning. Based on these observations, particularly the lithospheric thickness, we conclude that Archean, Paleoproterozoic, and Mesoproterozoic lithosphere are unique and together form the stable cratonic nuclei, defined as the thick (150–250 km), long-lived cores of continents. Higher mantle temperatures during the Archean, Paleoproterozoic, and Mesoproterozoic may have played a key role in the formation of the thick lithospheric roots. The second key finding is that the similarities of mean seismic properties indicate that the process of crustal formation operating in the Archean continued during the Paleoproterozoic and Mesoproterozoic. The thin (mean value < 135 km) lithospheric root beneath Neoproterozoic and younger crust may be related to the steady decrease in mantle temperature through time. Neoproterozoic and Paleozoic crust have similar physical properties, and these eras are characterized by pronounced biodiversification, including the renowned Garden of Ediacara, the Cambrian Explosion, and the Great Ordovician Biodiversification. Mesozoic–Cenozoic crust is the most diverse and reflects the current tectonic and magmatic processes of crustal formation.

Book chapter

Pre-eruptive characteristics of “suspect” silicic magmas in Carlin-type Au-forming systems

World-class Carlin-type Au deposits hosted in sedimentary rock were formed when profuse Eocene silicic magmatism swept across northern Nevada in response to arc migration. Carlin-type Au deposits formed along with porphyry/skarn Cu-Mo-W-Au deposits, epithermal Ag-Au deposits, and distal disseminated Ag-Au deposits. But unlike these other Au-bearing deposits that have clear associations with igneous intrusions, Carlin-type ore deposits appear to have formed distant from concealed plutons, and their origin remains controversial. Despite decades of abundant geophysical, geochronological, and geochemical studies suggesting the involvement of magmas, concrete evidence for magmatic involvement is still lacking. Consequently, the involvement of contemporaneous igneous systems remains inferred based on age, proximity, and variable isotopic, geochemical, and geophysical clues. A recent synthesis of deposit models postulates that Carlin-type Au deposits are intrusion-related, but that the causative magmas reside deeper (∼6–12 km) than in typical porphyry and peripheral systems (∼3–5 km), meaning that Carlin-type deposits are perhaps more distal expressions of igneous intrusions. We investigate a collection of “suspect” magmatic systems over a ∼7 m.y. timespan (∼41–34 Ma) that are contemporaneous with and near known Carlin-type ore deposits. We report results of a multifaceted array of in situ geochemical analyses (FTIR, EMP, SHRIMP-RG, LA-ICP-MS) of quartz-hosted melt inclusions, biotite, and quartz to better characterize the pre-eruptive characteristics of these magmas. We also report results of thermobarometry and thermodynamic phase equilibria modeling to help place constraints on magmatic reservoir depths and processes. Rather than a single “flavor” of silicic magma, we observe a surprisingly broad compositional spectrum of rhyolites, with one end of the spectrum exhibiting more arc-like (I-type) characteristics and the other end displaying more post-subduction, thick-crust extensional (A-type) characteristics. This broad compositional spectrum suggests a more complex picture of silicic crustal magmatism operating over a narrow span of time during slab rollback. Despite this spectrum, magmatic systems in this study are consistently ferroan and generally peraluminous, which we interpret as an expression of the relatively elevated geotherm at the time and incorporation of variable amounts of highly peraluminous metasedimentary crustal components. The silicic magma spectrum encompasses a range of mineralization associations, including subduction-related Cu-Mo-W-Au-Ag and post-subduction, thick-crust extensional rare-metal Mo-Sn-W-F-Be-Ag-Au, consistent with the prolific and diverse array of ore deposits that formed during this time. Carlin-type Au deposition appears to be associated with nearly the entire magmatic spectrum. This apparent indifference to silicic magma “flavor” would seem to imply that if magmas are involved in Carlin-type Au deposit genesis, they perhaps do not need to be compositionally specialized and/or possibly are only relevant as heat sources driving circulation to remobilize and redistribute metals.

Nevada

Simulation of the effects of development of the ground-water flow system of Long Island, New York

Extensive development on Long Island since the late 19th century and projections of increased urbanization and ground-water use makes effective water-resource management essential for preservation of the island's hydrologic environment and maintenance of a reliable source of water supply. This report presents results of a ground-water flow simulation analysis of the effects of development on the Long Island ground-water system. It describes ground-water levels, stream-flow, and the ground-water budget for the predevelopment period (pre-1900), the 1960's drought, and a more recent (1968-83) period with significant hydrologic stress. The report also presents estimated effects of a proposed water-supply strategy for the year 2020. Long Island has three major aquifers-the upper glacial (water-table), the Magothy, and the Lloyd aquifers-that are separated to varying degrees by confining units. Before development, recharge from precipitation entered the ground-water system at a rate of more than 1.1 billion gallons per day. An equal amount discharged to streams (41 percent), the shore (52 percent), and subsea boundaries (7 percent) . Urbanization and withdrawal of more than 400 Mgal/d (million gallons per day) from wells have resulted in local effects that include declines in ground-water levels, drying up and burial of streams and wetlands, reduction of ground-water recharge by increased overland flow to the ocean, a general decrease in ground-water discharge, and salt water intrusion. In some areas, the reduction in recharge is mitigated by leakage from water-supply and wastewater disposal lines, and infiltration of storm water through recharge basins. During 1968-83, a net loss of 240 Mgal/d from the ground-water system caused a decrease in ground-water discharge to streams (135 Mgal/d), to the shore (82 Mgal/d), and to subsea boundaries (23Mgal/d).The greatest adverse effects have been in western Long Island, where the most severe development has occurred. This analysis shows stream base flow to be highly sensitive to water-table fluctuations, and long streams to be more sensitive than short ones. A water-supply scenario for the year 2020 was simulated that employs redistribution of pumping centers to mitigate extreme local effects . Although the net stress on the ground-water system was projected to increase 57 Mgal/d (24 percent) above that of 1968-83, redistribution of ground-water withdrawals across the island would allow recovery of cones of depression in western Long Island, thereby reducing the threat of salt water intrusion and increasing base flow of some streams . The increased stress would cause a net decrease in base flow island wide of 44 Mgal/d; total base flow would be 281 Mgal/d - 39 percent below predevelopment levels or 14 percent below 1968-83 levels. The most severe effects would be in Nassau and western Suffolk Counties.

New York

Structural controls on splay fault rupture dynamics during Cascadia megathrust earthquakes

Great subduction earthquakes ( M w ≥ 8.0) can generate devastating tsunamis by rapidly displacing the seafloor and overlying water column. These potentially tsunamigenic seafloor offsets result from coseismic fault slip and deformation beneath or within the accretionary wedge. The mechanics of these shallow rupture phenomena and their dependence on subduction zone properties remain unresolved, partly due to the sparsity of offshore observations of shallow megathrust earthquake deformation. Here, we analyze how offshore structure influences shallow rupture mechanics and slip partitioning using 3D dynamic earthquake simulations of the Cascadia subduction zone (CSZ) megathrust with and without variably dipping seaward- or landward-vergent splay faults in the wedge that sole into the megathrust. Resulting tradeoffs between splay and megathrust slip reveal structural controls on rupture partitioning, with greater splay slip leading to less shallow megathrust slip updip. Gently dipping and seaward-vergent splays host more slip than those with steeper, landward-vergent splays. To isolate the underlying mechanisms, we compare models with Andersonian and plunging principal stresses. Results suggest distinct static and dynamic processes control the dip- and vergence-dependence of splay rupture: static (mis)alignment relative to far-field tectonic loading favors slip on more optimally oriented, shallowly dipping splay faults. In contrast, dynamic stress interactions of an updip-propagating megathrust rupture front with the free surface and potential branch faults favor forward branching onto seaward-vergent splays and inhibit backward branching onto landward-vergent splays. Resulting seafloor displacements suggest splay fault structure may influence coseismic tsunami source processes, highlighting the importance of dynamically viable rupture scenarios in subduction hazard assessments.

Cascadia subduction zone

Comparative assessment of STIC sensors, streamflow and rain gauges for quantifying river connectivity in intermittent systems

In intermittent stream systems, including those occurring in Texas, USA, the severity of low-flow conditions, duration of seasonal disconnection, and frequency of no-flow events have been amplified by drought. Documentation of these no-flow events is necessary to evaluate ecosystem health. However, many intermittent reaches remain un-gauged given that perennial river sec-tions are often prioritized for gauge placement. Our objectives were to 1) document stream flow using Stream Temperature, Intermittency, and Conductivity (STIC) loggers to determine the frequency and duration of no-flow events in intermittent tributaries of the Colorado River, Texas and 2) compare logger data to publicly available data from streamflow discharge and precipitation gauge networks to understand differences among these data types for drying event characterization. We use these comparisons to summarize benefits and limitations of the application of in-stream data loggers. STIC loggers were deployed at 19 sites, one in each pool and riffle habitat of a stream reach. STIC loggers recorded a measurement of relative conductance every six hours from June 2022 to March 2024, which was used to determine the presence or absence of flow connectivity in a reach. No-flow duration among intermittent reaches varied between 37 and 270 days across tributaries during an ongoing drought in the study area. Overall, logger data was more precise than discharge data for characterizing no-flow events or precipitation data when documenting presence of water in the stream channel due to runoff. Lack of discharge gauges in intermittent tributaries left large sections of stream reaches undocumented and resulted in mischaracterization of flow patterns. Drought severity across the tributaries did not follow longitudinal patterns that would be expected by the climatic precipitation gradient of the study area. More research is needed to determine if factors such as population size affect severity. Likewise, precipitation data did not correlate well with logger water presence data, lacking consideration for groundwater recharge, soil hydrophobicity, and surface compaction. This study shows that to monitor no-flow events, detailed spatial datasets are necessary and that STIC loggers are useful tools that provide data to fill spatial information gaps and facilitate more accurate flow characterization and water presence data in intermittent systems.

Texas

Simulated ground-motion records for the seismic assessment of monumental masonry structures

Earthquakes are natural disasters that can cause widespread devastation and loss of life. Simulated ground-motion records can be useful in regions with limited seismic stations or a history of damaging but infrequent earthquakes. This is especially true in areas with a high concentration of heritage masonry structures, which are especially susceptible to damage, as simulated records can be crucial in predicting their seismic response and protecting these buildings from seismic damage. Despite the importance of simulated earthquakes, few studies have investigated how effective they are compared to real earthquakes when assessing the structural response of heritage buildings. To address this knowledge gap, we employ two different simulation methods of the Mw 6.2 Faial earthquake, which occurred on July 9th, 1998, in the Azores, to replicate the recorded time-series at four available stations within an epicentral distance of 150 km. The study has two objectives: first, to validate the simulated records of the 1998 Faial earthquake using alternative stochastic ground-motion simulation approaches, and second, to determine how these approaches impact the seismic assessment of historic masonry structures. To accomplish these objectives, this study uses real and simulated ground-motion datasets to conduct non-linear response history analyses of the São Francisco Church, a monumental structure in Horta that sustained damage during the Faial earthquake. The results show that both simulation approaches yield structural responses similar to the observed records.

Conference Paper

Transmissivity and water quality of water-producing zones in the intermediate aquifer system, Sarasota County, Florida

The intermediate aquifer system is an important water source in Sarasota County, Florida, because the quality of water in it is usually better than that in the underlying Upper Floridan aquifer. The intermediate aquifer system consists of a group of up to three water-producing zones separated by less-permeable units that restrict the vertical movement of ground water between zones. The diverse lithology, that makes up the intermediate aquifer system, reflects the variety of depositional environments that occurred during the late Oligocene and Miocene epochs. Slight changes in the depositional environment resulted in aquifer heterogeneity, creating both localized connection between water-producing zones and abrupt culmination of water-producing zones that are not well documented. Aquifer heterogeneity results in vertical and areal variability in hydraulic and water-quality properties. The uppermost water-producing zone is designated producing zone 1 but is not extensively used because of its limited production capability and limited areal extent. The second water-producing zone is designated producing zone 2, and most of the domestic- and irrigation-supply wells in the area are open to this zone. Additionally, producing zone 2 is utilized for public supply in southern coastal areas of Sarasota County. Producing zone 3 is the lowermost and most productive water-producing zone in the intermediate aquifer system. Public-supply well fields serving the cities of Sarasota and Venice, as well as the Plantation and Mabry Carlton Reserve well fields, utilize producing zone 3. Heads within the intermediate aquifer system generally increase with aquifer depth. However, localized head-gradient reversals occur in the study area, coinciding with sites of intense ground-water withdrawals. Heads in producing zones 1, 2, and 3 range from 1 to 23, 0.2 to 34, and 7 to 42 feet above sea level, respectively. Generally, an upward head gradient exists between producing zones 3 and 2. However, an upward head gradient between producing zones 2 and 1 does not consistently occur throughout Sarasota County, probably the result of greater ground-water withdrawals from producing zone 2 than from producing zone 1. The transmissivity of the intermediate aquifer system is spatially variable. Specific-capacity data from selected wells penetrating producing zones 2 and 3, were used to estimate transmissivity. Estimated transmissivity values for producing zones 2 and 3 range from about 100 to 26,000 feet squared per day and from about 1,300 to 6,200 feet squared per day, respectively. Because the capacity of specific water-producing zones is highly variable from site to site, estimating the performance of a specific water-producing zone as a water resource is difficult. Water samples collected during the study were analyzed for major-ion concentrations. Generally, bicarbonate type water from rock interaction occurs in northern Sarasota County; enriched calcium-magnesium-sulfate type water from deeper aquifers occurs in central Sarasota County; and sodium-chloride type water from saltwater mixing occurs in southern Sarasota County. In some areas of northern Sarasota County, the major-ion concentrations in water are lower in producing zone 2 than in producing zone 1. Major-ion concentrations in water are higher in producing zone 3 throughout the study area. A major objective of the study was to evaluate hydraulic and water-quality data to determine distinctions that could be used to characterize a particular producing zone. However, data indicate that both hydraulic and water-quality properties are highly variable within and between zones, and are more related to the degree of connection between and areal extent of water-producing zones than to aquifer depth and distance from the coast.

Florida

A novel drive-point multilevel system to investigate PFAS and other contaminants of global concern in the hyporheic zone of a wastewater effluent dominated stream

Contaminants found in treated wastewater discharged to streams, including pharmaceuticals and per- and polyfluoroalkyl substances (PFAS), are of global concern due to their deleterious effects on aquatic ecosystems and potential impacts to human health. Hyporheic zones have strong potential for contaminant attenuation. Assessing this potential requires collection of physical and biogeochemical data within the hyporheic zone. This study tested the applicability of a novel drive-point multilevel system (DP-MLS) for quantifying head profiles and characterizing contaminant concentrations in the hyporheic zone of a temperate region effluent dominated stream (EDS). DP-MLS, each with 4 ports, were installed in the stream bed at two sites, DS-1 and DS-2, 0.2 and 4.7 km downstream of the effluent outfall, respectively. Head profiles were measured and groundwater collected for analysis of pharmaceuticals and PFAS temporally over two years. The DP-MLS withstood rapid changes in stage, ice formation, and floating debris. Vertical hydraulic gradients (VHG) were generally upward but varied in magnitude indicating heterogeneity in hydraulic conductivity and variability in flow conditions. Upward VHG were also about 2X larger at DS-1 than at DS-2. Contaminant concentration profiles consistently showed penetration of pharmaceuticals and PFAS to 1 m below the bed at DS-2 while there was less penetration, lower groundwater concentrations, and more temporal variability in concentrations at DS-1. Integration of the physical and chemical data suggests weaker upwelling conditions at DS-2 are more easily reversed during periods of high stream stage, which could facilitate migration of wastewater contaminants into the bed. However, further studies incorporating other transport processes and reach scale dynamics are required to fully characterize these exchanges. Overall, this study demonstrates the efficacy of these novel DP-MLSs for characterization of the hyporheic zone and provides new insights into the occurrence, composition, and persistence of wastewater derived contaminants in the hyporheic zone of a well-studied EDS.

Iowa

Hidden system identification: Basin modeling as a tool for examining sedimentary geothermal resource potential

Three-dimensional (3D) geologic and temperature models have been developed for the onshore U.S. Gulf Coast. The results from these models identify areas of moderate- to high-temperature (90°-150°C and >150°C; respectively) geothermal resources at depths <6 km. This modeling study addresses the fundamental challenge of predicting where opportune temperature and lithology coincide. Unlike traditional geothermal systems with surface expressions of hydrothermal circulation (e.g., hot springs, fumaroles, sinter), sedimentary geothermal systems (SGS) are generally hidden. Historically, simplified efforts to predict subsurface temperatures in sedimentary basins have focused on linear temperature extrapolation that does not consider the variable thermal properties of different lithologies or lithologic changes with depth (e.g., compaction, lithification). Therefore, the need to understand basin architecture and predict temperatures in 3D within SGS is paramount to identifying geothermal resources and determining economic feasibility. Basin modeling software has long been used to characterize the subsurface conditions of sedimentary basins, including temperature, in the pursuit of finding hydrocarbons. This tool can also be adapted to evaluate the potential of geothermal resources in a sedimentary basin by predicting the confluence of desirable temperatures and reservoir lithologies. In this work, PetroMod basin modeling software was used to create a regional geologic model of the onshore U.S. Gulf Coast, covering over 500,000 km 2 calibrated to temperature data from wells. Inputs include structural surfaces from commercial databases, lithology information derived from published literature, and corrected bottom-hole temperatures (BHT) from over 6,000 wells. The resulting 3D geologic model can be used to predict temperatures throughout the basin. Maps were exported showing the depth, depositional unit, and reservoir lithology at which temperatures of 90°C and 150°C were reached, revealing over 400,000 km 2 of moderate- to high-temperature resources at depths <6 km. These maps function as a first-order screening tool to identify areas where low-, moderate-, or high-grade resource potential may exist, based on temperature and if optimal reservoir lithologies or depositional units of interest are present. Depending on the success criteria of a project, the same maps can be exported for any isotherm or incorporate other 1407 Gardner and Birdwell subsurface properties. The methodology employed in this work can be applied in any sedimentary basin with available subsurface data. Further calibration incorporating other data, including pressure and porosity, can expand the utility of basin modeling for geothermal evaluations. Basin modeling is a powerful but underutilized tool for identifying prospective geothermal resources in sedimentary basins.

Conference Paper

The systematics of stable hydrogen (δ2H) and oxygen (δ18O) isotopes and tritium (3H) in the hydrothermal system of the Yellowstone Plateau volcanic field, USA

To improve our understanding of hydrothermal activity on the Yellowstone Plateau volcanic field, we collected and analyzed a large data set of δ 2 H, δ 18 O, and the 3 H concentrations of circum-neutral and alkaline waters. We find that (a) hot springs are fed by recharge throughout the volcanic plateau, likely focused through fractured, permeable tuff units. Previous work had stressed the need for light δ 2 H water recharge restricted to the northern part of the plateau or recharge during past cold periods. However, new data from the Y-7 drill hole suggests that recharge is not restricted to a certain area or a cold period. (b) δ 18 O values of thermal waters in the geyser basins are shifted from the global meteoric water line by temperature-dependent water-rock reactions with higher subsurface temperatures resulting in a greater shift. (c) Large temporal variations in the isotopic composition of meteoric water recharge and small temporal variability in the isotopic composition of hot spring discharge implies that the volume of groundwater in, and around the Yellowstone caldera is substantially larger than the volume of annual water recharge. (d) Hot springs discharged through different rhyolitic units correlate with identifiable differences in δ 2 H and δ 18 O compositions, 3 H concentrations, and water chemistry that imply equilibration at different temperatures and travel along different flow paths. (e) Based on measured 3 H concentrations, we calculate that hot spring waters in the central part of the geyser basins mostly contain <2% post-1950 meteoric water, whereas waters discharged at the basin margins contain larger fractions of post-1950s meteoric water.

Wyoming

High frequency and region-scale simulations of large (Mw7+) earthquakes on the southern Whidbey Island fault, Washington, USA

We simulate ground shaking in western Washington State from hypothetical M w 7.0–7.5 earthquakes on the southern Whidbey Island fault (SWIF). Ground motions are modeled considering kinematic source distributions on a complex fault plane, a 3D seismic velocity model, and region‐specific soil velocity models. We run simulations with varying model resolutions, including regional‐scale simulations with a maximum‐modeled frequency of ∼1 Hz and local‐scale simulations with a maximum‐modeled frequency of ∼2.5 Hz. Additional local‐scale simulations are run considering high‐resolution surface topography. We explore how source parameters (i.e., magnitude, hypocenter location, and dip direction) and 3D velocity structure impact peak shaking intensity and its variability. In particular, we find that earthquakes on the SWIF would likely produce strong shaking throughout the populated Puget Lowland, including in the cities of Everett, Seattle, Bellevue, and Tacoma, Washington. Simulated short‐period ( T ≤ 2 s) spectral accelerations are strong throughout the Puget Lowland, and long‐period shaking ( T ≥ 5 s) is strong in the deep regional sedimentary basins, especially the Everett and Seattle basins. Source parameters strongly influence intra‐ and interevent variability in response, primarily through changes in source and site geometry, as well as rupture directivity. We also note a potential coupling between rupture directivity and basin effects, wherein directivity pulses are seemingly guided through the region’s deep, interconnected sedimentary basins. Overall, this work highlights the impacts of 3D source, path, and site effects on seismic hazard in the U.S. Pacific Northwest and substantially expands the catalog of simulated ground motions for Puget Sound area crustal faults.

Washington