USGS ScienceSearch

SEARCH · USGS Science

Results for “Earth Surface Processes and Landforms”

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,665 records · Page 79Linked to original sources

Geochemical disequilibrium at the brittle-ductile transition

We investigate the microtextural, microchemical, and isotopic effects of late-stage ductile deformation in quartzite mylonites and kyanite–muscovite–quartz veins from the Raft River shear zone (Utah). Quartz microstructures record pervasive disequilibrium, expressed by unannealed features including undulatory extinction, deformation lamellae, and poorly defined fabrics, typical of waning deformation in shear zones. Microchemical disequilibrium is best preserved in kyanite quartzites, where CL-zoned kyanite records repeated fractures, overgrowth, and mineral precipitation, and in muscovite from muscovite-poor quartzite mylonites that shows minor-element zoning consistent with syn-deformational overgrowth on detrital cores. In contrast, muscovite from muscovite-rich kyanite quartzites exhibits minimal chemical zoning. These microchemical variations correlate with 40 Ar/ 39 Ar age systematics. Chemically zoned muscovite preserves variable single-step ages, including ∼150 Ma ages, reflecting retention of detrital cores. In contrast, syndeformational muscovite consistently yields Miocene ages, indicating recrystallization and new growth below argon closure temperatures that reset inherited isotopic signatures. Similar trends are observed in quartzite mylonites, where increasing quartz recrystallization and stronger crystallographic preferred orientations occur toward deeper structural levels. Together, these observations indicate increasing retrograde deformation and recrystallization with depth in the Raft River shear zone and demonstrate that strain-driven recrystallization exerts a first-order control on muscovite ages. We suggest that apparent thermochronologic gradients in retrograde shear zones may reflect recrystallization gradients rather than temperature gradients. Where cooling limits the thermal driving force for recrystallization, isotopic relics are preserved, and local deformation and fluid availability control re-equilibration. Consequently, isotopic disequilibrium—particularly in the 40 Ar/ 39 Ar system—may be the rule rather than the exception of retrograde tectonic environments.

Idaho, Utah

A new water temperature modeling approach to predict thermal habitat suitability for nonnative cichlids in Florida rivers

As global temperatures increase, the spatiotemporal arrangement of thermal habitats in Florida rivers may shift, creating the potential for greater dispersal and establishment of nonnative tropical freshwater fishes. To understand how water temperature changes may affect the spatial distribution of these nonnative species, more effective water temperature prediction models are necessary. Currently, most models employ either a generalized air–water temperature relationship or require expensive and complicated tools to measure hydrometeorological factors (e.g. groundwater input). Thus, we developed a novel modeling approach that is accurate, accessible, and cost-effective in allowing fisheries managers to project water temperatures in rivers across Central and North Florida. To characterize the potential for nonnative fishes to spread northward, we evaluated two hardy and abundant species currently found primarily in South Florida: Mayan Cichlid ( Mayaheros urophthalmus ) and Oscar ( Astronotus ocellatus ). Our results show an increase in thermally suitable winter days for both species in 10 of 11 rivers studied, consistent with predicted water temperature warming under 16 climate-change scenarios spanning different levels of air temperature warming (+1 °C, +2 °C, +3 °C, +4 °C) and precipitation/groundwater thermal sensitivity (0, 0.33, 0.66, 1). Considering resource limitations, fisheries managers can use our water temperature modeling approach to predict effects of climate change on Mayan Cichlid and Oscar survival, growth, and dispersal and take actions to manage potential northward movement of these species.

Florida

Statewide surficial geologic map of Nebraska underscores Quaternary landscape evolution from the High Plains to the Central Lowland

Surficial geologic mapping in Nebraska has been conducted primarily at the 1:24,000 scale since the mid-1990s, although there have also been limited efforts to map generalized Quaternary and surficial geology within the state and the region. We compiled and evaluated disparate maps (1:24,000–1:1,000,000) and datasets—including geologic, soil and soil parent material, and geomorphic maps as well as LiDAR derivatives—to produce a single 1:500,000 scale surficial geologic map of Nebraska that is the first of its kind. This new map provides a coarse-scale surficial geologic map that will be incorporated into a nationwide U.S. Geological Survey Quaternary geologic map. It also reflects the variation and uniqueness of physical landscapes in the state, where the Great Plains and Central Lowland physiographic provinces meet, further developing a richer interdisciplinary understanding of regional geomorphology in the heart of North America.

Nebraska

Effects of stochastically-simulated near-fault ground motions on soil liquefaction

The scarcity of historically recorded near-fault ground motions poses a challenge to systematically understanding the influence of near-fault effects on various types of seismic demands for engineering purposes. In particular, the current state of knowledge of the influence of ground-shaking intensity on soil liquefaction and its consequences does not specifically account for the effects of near-fault ground motion characteristics. In this study, the influence of near-fault ground motions on liquefaction triggering and lateral spreading are investigated using non-linear modeling of a hypothetical liquefiable soil column in the finite-element computational platform OpenSees subjected to simulated ground motion time series that represent strong earthquake shaking in the near field. The simulated ground motion time series and resulting datasets are based on a parametric stochastic model and are developed for a range of source and path parameters to represent a realistic variability of ground motion characteristics. Dependencies between ground motion intensity measures (IMs) and liquefaction demand parameters are investigated for near-fault pulse and nonpulse-like ground motion sets. Evolutionary IMs, such as cumulative absolute velocity (CAV) and the time-varying magnitude-adjusted peak ground acceleration (PGAM), are considered in developing liquefaction triggering probability density functions. Post-liquefaction triggering responses such as lateral spreading displacements are examined in relation to PGA M and CAV. The ground motion simulations are validated by comparing their liquefaction-capacity PGA M fragilities and post-triggering CAV vulnerability relationships to historical records from the 1994 Northridge earthquake in California, USA. Finally, a path forward for future studies that includes finding systematic differences in the IM-liquefaction demand relationships between near-fault and far-field stochastic ground motion sets is outlined.

Conference Paper

Forecast, monitor, adapt: A multi-agency strategy to protect people from postfire debris flows

In 2020, a wildfire burned across Glenwood Canyon in Colorado, USA. A history of postfire debris flows in the region and a hazard assessment for the burn area indicated that potentially life-threatening debris flows could be triggered by rainfall within months of a wildfire. As a result, four government agencies evaluated strategies to help mitigate hazards, including the loss of human life, that may be associated with debris-flow events. After the fire, 26 large debris flows occurred in the summer of 2021 and three sediment-laden flows occurred in the summer of 2023, but there were no major injuries or fatalities reported. We found that integrating hazard assessment/ forecasting, monitoring, and adaptation scenarios was a successful strategy for reducing postfire debris-flow risks to human life (including injuries and fatalities). Weather forecasts and estimates of debris-flow triggering rainfall thresholds, likelihood, and volume were used to anticipate the timing, location, and magnitude of debris-flow events. Rainfall monitoring and detailed recordkeeping of storms that triggered debris flows were used to validate and update debris-flow warning thresholds that varied with time following the wildfire. Although the governmental agencies working in this burn area had distinct and differing agency mandates, they were able to integrate information to reduce the risk of debris-flow events to human life.

Colorado

Seismic response of rock towers at the Trona Pinnacles (U.S.A.) to the 2019 Ridgecrest earthquake sequence: Theory, observations, and models

We analyze the seismic response of a class of fragile geologic features (FGFs), referred to as rock towers (RTs) at the Trona Pinnacles, a group of RTs in southern California that suffered strong shaking during the 2019 Ridgecrest earthquake sequence. FGFs, including RTs, may provide maximum constraints on past earthquake shaking intensity, and thereby support probabilistic seismic hazard assessments (PSHAs). In a rare case study, we explore the hypothesis that RT structural integrity is time dependent, as damage accumulates progressively. We develop finite‐element method (FEM) models of the RTs using photogrammetric shape models. We validate the models by comparing numerical simulations of their response to broadband ground shaking with low‐intensity seismic recordings obtained at the Pinnacles. Results of our simulations are in good agreement with the seismic recordings of actual earthquake aftershocks. We next use the results of the FEM models to analyze the response and evolution of RTs. Our analyses elucidate the influence of geometry over their seismic response, providing a rationale that may explain the rarity of slender RTs at Trona: high‐aspect‐ratio structures that respond in bending develop detrimental tensile stresses that crack the rock, whereas low‐aspect‐ratio ones’ response also includes shearing, which does not compromise material integrity as much as tension. Field measurements with a rebound hammer support this finding, suggesting that the material around the base of slender rocks has been weakened relative to other parts of the RT. We also study how to define simplified mechanical models (“archetypes”) to predict the natural frequencies of RTs. Results from our work illuminate the fundamental mechanisms of seismic response and progressive failure of RTs, and open new avenues of research to potentially incorporate these geologic features as long‐return period constraints on PSHA, in ways analogous to those of the widely used precariously balanced rocks.

California

California State Waters Map Series—Benthic habitat characterization in the region offshore Humboldt Bay, California

Coastal and Marine Ecological Classification Standard (CMECS) geoform, substrate, and biotic component geographic information system (GIS) products were developed for the California State Waters of northern California in the region offshore of Humboldt Bay. The study was motivated by interest in development of offshore wind-energy capacity and infrastructure in Federal waters offshore. This project, carried out by the U.S. Geological Survey (USGS), resulted in four data releases for individual map blocks that are part of the “California State Waters Map Series”: (1) Offshore of Arcata, (2) Offshore of Eureka, (3) Offshore of the Eel River, and (4) Offshore of Cape Mendocino. The study area consists of 436 square kilometers of multibeam echo sounder (MBES) data acquired by Fugro Pelagos, Inc., in 2007. Towed camera-sled video was acquired in 2009 and 2010 to supervise the classification of the MBES data into habitats, and single channel sparker data were collected to calculate sediment thickness above the transgressive unconformity. Using video observations of habitat as ground truth, derivatives of the MBES data were classified into 3 seafloor character types (hard-rugose, hard-flat, and soft-flat), 26 induration-slope-depth groups, and 15 geoforms. The study area substrate is predominantly soft-flat sediment (mud and fine sand) covering 73.6 percent of the area. Hard-flat substrate areas, predominantly coarse sediment in scour depressions, cover 5.4 percent of the study area. The hard-rugose substrate areas are primarily outcrops of layered sedimentary bedrock and constitute 20.9 percent of the study area. Fifteen geoforms were identified in the analysis. The predominant geoforms mirror the seafloor character results, shelf geoforms, rock outcrop geoforms, and scour depression geoforms. Rock and scour areas are restricted to the southern portion of the study area off Cape Mendocino where uplift has exposed bedrock. On the flat shelf area post-transgressive sediment varies in thickness from 1.7 meters (m) nearshore to 28.1 m offshore.

California

Bedrock geologic map of the Woodstock quadrangle, Grafton County, New Hampshire

The bedrock geology of the Woodstock 7.5-minute quadrangle consists of highly deformed metasedimentary rocks of the Central Maine trough, including the Silurian Rangeley and Perry Mountain Formations and the Devonian Littleton Formation. The central, northern, and eastern parts of the quadrangle are underlain by the oldest rocks in the area, the Rangeley Formation. In the southwest and south-central part of the quadrangle, metaturbidites of the Perry Mountain Formation and subsequent Littleton Formation overly the Rangeley Formation in a deformed F 1 synform, herein informally called the Bagley Brook basin. The metasedimentary rocks were intruded by widespread syn- to post-tectonic granitoids of the Devonian New Hampshire Plutonic Suite and minor post-metamorphic Jurassic-Cretaceous mafic dikes of the White Mountain Plutonic-Volcanic Suite. The metasedimentary rocks were affected by at least two episodes of deformation in the Devonian Acadian orogeny. The dominant regional foliation is second-generation (S 2 /D 2 ) and formed during the development of sillimanite-muscovite mineral assemblages. Large bodies of the Early Devonian Kinsman Granodiorite intruded the metasedimentary rocks semi-concordantly during D 2 deformation. Dikes of the Late Devonian Concord Granite cut the Kinsman Granodiorite and the metasedimentary rocks and were emplaced either syn- or post-D 2 . The map pattern in the Rangeley Formation is dominated by northeast to northwest trending, moderately to steeply north-dipping F 2 and F 3 folds. Map-scale F 1 folds are defined by the Bagley Brook basin. Previous division of Rangeley Formation stratigraphy in this region into “upper” and “lower” parts was not corroborated by 1:24,000-scale mapping of lithodemic units, and rocks previously mapped as part of the Smalls Falls and Madrid Formations are here reassigned to the Rangeley Formation. Some rocks previously mapped as the lower part of the Littleton Formation are now assigned to the Perry Mountain Formation. The Littleton Formation on this map is approximately equivalent to rocks previously mapped as the upper part of the same formation. Steeply dipping fractures in the quadrangle show a preferred northeast orientation, consistent with subsurface fracture orientations in the well fields near Mirror Lake. Jurassic-Cretaceous mafic dikes and normal faults show preferred northeast orientations, similar to the fractures, suggesting that the extensional stress field that controlled dike orientation during the Mesozoic also produced the dominant brittle fabrics in the area.

New Hampshire

The geology of Canadian potash: A critical mineral for feeding the world

Potash, potassium-bearing water-soluble salt, is the primary global economic source of potassium. Potash is recognized as a critical mineral in Canada as it is the largest source of potassium used in fertilizers. It is essential for global agricultural productivity and food security. Canada is the world’s largest potash exporter with vast deposits in the widely mined Prairie Evaporite of Saskatchewan, which formed in the epicontinental Elk Point Basin during the Middle Devonian. Potash is also found in the Windsor Group of Atlantic Canada where it formed in a series of tectonically active basins during the Mississippian that have undergone substantial post-depositional subsurface alteration and deformation. Potash deposits were mined in New Brunswick up until 2016. Both deposits are salt giants, recording times in the geologic record of extensive and long-lasting evaporite genesis under arid conditions in restricted seas. This paper reviews the geological and economic significance of Canadian potash, including (1) the genesis of each deposit, (2) diagenetic, erosional, and tectonic modification, and (3) exploration and mining in each basin. Underdeveloped regions, possible undiscovered resources, environmental considerations, and the importance of sustainable practices in light of climate change and socioeconomic risks are also addressed.

Facets

Origin of the Pd/Pt ratio of the J-M Reef, Stillwater Complex, Montana, USA

The J-M reef of the Stillwater Complex is characterized by a high Pd/Pt ratio (mean ~3.8 with a standard error of 0.03) with a homogeneous geospatial distribution at the deposit scale. In this contribution, we demonstrate that the Pd/Pt ratio of the reef is the product of equilibration of an immiscible sulfide liquid with a silicate melt rich in Pd relative to Pt. Despite the high tenors of the J-M reef sulfides (avg 2,700 ppm Pt and 770 ppm Pt), numerical modeling shows that the parental melts did not have extraordinary Pd and Pt concentrations. Instead, the initial composition of a plausible parental silicate melt can have Pd and Pt contents well within the expected range of a normal, mantle-derived partial melt (i.e., ~10–20 ppb for both Pd and Pt with Pd/Pt of ~1). The relative differences in the partitioning behavior of Pt and Pd between sulfide liquid and silicate melt are unlikely to produce a consistent Pd/Pt ratio across a wide range of silicate melt to sulfide liquid mass ratios (i.e., R factors). Instead, the pre-emplacement fractionation of Pt alloy from S-undersaturated silicate magma accounts for the homogeneous and high Pd/Pt ratio of the J-M reef. We show that batch equilibration of sulfide liquid with silicate melt can produce the high Pd/Pt ratios of the reef if the partition coefficients between sulfide liquid and silicate melt for Pd and Pt are extremely high (>10 6 ). In an alternative model, Pd enrichment could be achieved by sulfide upgrading in resident footwall mush even if the partition coefficients between sulfide liquid and silicate melt are relatively small (between 10 4 and 10 6 ) because the instantaneous mass ratio of silicate melt to sulfide liquid is small (R ≈ 100–700), so the partitioning behavior of Pt and Pd has little impact on the composition of sulfide liquid.

Montana

Subsurface vertical connectivity shapes solute transport to montane streams: Insights from watershed-scale geophysics-informed modeling

Headwater catchments provide essential water and nutrients to downstream ecosystems. The timing and evolution of their transport is shaped in part by the “invisible” subsurface structure but have been largely unexplored at the watershed scale. This raises a fundamental question: how and to what extent does subsurface structure influence streamflow and solute transport in montane watersheds? Here, we test controls of three-dimensionally resolved subsurface structure on solute transport, which reflects source waters and water residence time. We integrate the geophysically mapped architecture of Coal Creek, a headwater catchment in the Colorado River Basin, into numerical flow-transport models. Our results reveal greater vertical connectivity elongates flow paths and enhances deep groundwater contributions to streams, therefore substantially influencing the timing and magnitude of solute transport. In contrast, watershed-integrated outlet discharge can be predicted without high-resolution subsurface information. These findings demonstrate incorporating subsurface architecture is critical for predicting ecosystem health and contaminant transport.

Colorado

Simulation of groundwater flow in the Silurian aquifer, eastern Iowa, 2020–45

The Silurian aquifer is an important water source for municipalities, industry, and rural households and communities in eastern Iowa, including Johnson County. Increasing demand for groundwater from the aquifer indicated that analytical tools may be beneficial to quantify groundwater resources and inform water-management decisions for the aquifer. The U.S. Geological Survey, in cooperation with the Johnson County Board of Supervisors, developed conceptual and numerical groundwater models to simulate water levels in the Silurian aquifer, determine groundwater budgets, and forecast changes in groundwater levels through 2045. A MODFLOW numerical model was constructed to match water levels for 2020 through 2022. Model performance was assessed by comparing simulated and observed water-level hydrographs at several pumping and monitoring wells. Hydrograph comparison indicates that the model simulates water levels accurately at monitoring wells that are not near areas of pumping and generally overpredicts water levels near pumping wells, while accurately simulating the drawdown trend over time at pumping wells. This transient model framework was then used to predict water levels through 2045 in a scenario of moderate drought and increased groundwater withdrawals. Results from the predictive model indicate as much as 13 meters of additional drawdown in 2045 from 2020 water levels in the area of pumping wells.

Iowa

Coastal Science Navigator companion guide—Discover the U.S. Geological Survey coastal science products you need

The Coastal Science Navigator is an online gateway to a wide variety of U.S. Geological Survey (USGS) coastal change hazards-related information, data, and tools relevant to stakeholders’ scientific and decision-making needs. The products within the Coastal Science Navigator provide data related to past, present, and future threats to our coastlines. The filter search allows users to see all available products and identify relevant options by time scale, geographic scope, coastal hazard theme, and other filters. The guided search suggests products based on users’ answers to a short series of questions. A comprehensive summary is available for each product. The idea for the Coastal Science Navigator arose in 2020 in response to stakeholder feedback identifying the need for a central source for USGS coastal science information. It was published in July 2023 and initially included 55 products. Regular updates are planned to integrate other existing and new products. This guide introduces some of the many coastal change hazards-related products available through the USGS. In it, we showcase the products included in the Coastal Science Navigator’s initial publication in July 2023. While it is not representative of all the information, tools, and data available, we hope it serves as a compelling snapshot of what the USGS has to offer and encourages you to explore the Coastal Science Navigator to discover more of the products you need. To navigate this guide, the products have been organized by the time scale they are best suited for—past, present, or future—although many products cover multiple time scales. An additional section features software, one of the many product types available as filters within the Coastal Science Navigator. Other products include downloadable data, websites, and geonarratives (web pages that combine text, images, and interactive maps into narratives you can scroll through). Featured geographic scopes are also highlighted within this guide, detailing some of the many regions in which the USGS conducts research and illustrating another way to filter products within the Coastal Science Navigator.

Circular

Black rails in Baja California

Summarizing records from 1905 to 2025, we document the presence and distribution of the California black rail ( Laterallus jamaicensis coturniculus ) in northwestern México, particularly Baja California and the Colorado River Delta. This subspecies has a restricted distribution in the western U.S. and northwestern México. We found few records over the past century: sparse detections in the early 1900s and an increase in reported records after 2000. Over 120 years, we compiled 35 documented records of California black rail occurrence in Baja California with a minimum reported count of 56 rails. Our call-broadcast survey methods reconfirmed previously occupied sites and documented the species at a previously unreported site. Surveys in the Baja California portion of the Colorado River Delta yielded 20 detections across nine survey events between 2000 and 2022. In April 2003, standardized surveys detected 12 California black rails: 7 in Bahía de San Quintín, 3 in Arroyo San Telmo, and 2 in Arroyo El Rosario. Our recent survey work (2022–2025) recorded maximum monthly counts of 31 detections at Bahía de San Quintín and 15 detections at Estero de Punta Banda, the latter representing a previously undocumented location for the subspecies. These surveys also yielded the first photographic and acoustic documentation of California black rail in México. The contribution of spatial and temporal data provided by citizen science platforms increased our understanding of the subspecies’ distribution. Enforceable protection of this subspecies’ habitat is needed to safeguard its persistence.

Baja California

U.S.-Mexico Borderland & vegetation community map

People on both sides of the United States-Mexico border need a high-resolution, binational vegetation community map that spans the entire United States-Mexico borderlands. Traditionally, mapping efforts in this region were impeded by complex logistics related to the international border, differing national needs and plans, and resource allocations and priorities. To address this need, scientists from the U.S. Geological Survey (USGS) Southwest Biological Science Center partnered with the Sonoran Joint Venture, the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program, data engineers from the Department of Biosystems Engineering at the University of Arizona, and collaborators from the Wildlands Network, the Borderlands Program to produce the first prototype land cover map within the overlapping Mojave Desert, Sonoran Desert, and the North American Bird Conservation Initiative’s Bird Conservation Region 33 (BCR33) using Landsat satellite data . BCR33 is an area of high biodiversity, providing habitat for bird species of concern and other wildlife. The land cover map supports FWS recovery plan efforts related to conservation planning activities for many species, including Yellow-billed Cuckoo ( Coccyzus americanus ), Cactus Ferruginous Pygmy-Owl ( Glaucidium brasilianum cactorum ), Southwestern Willow Flycatcher ( Empidonax traillii extimus ), Yuma Ridgway’s Rail ( Rallus obsoletus yumanensis ), Bendire’s thrasher ( Toxostoma bendirei ), LeConte’s thrasher ( Toxostoma lecontei ), Masked Bobwhite ( Colinus virginianus ridgwayi ), jaguar ( Panthera onca ), and endangered plants such as Bartram’s stonecrop ( Graptopetalum bartramii ) and the Pima pineapple cactus ( Coryphantha robustispina ssp. robustispina ). In 2024, a Phase-II map for the full BCR33 region was completed, increasing the understanding of the binational nature of natural communities. The published map and associated paper can be found here .

Borderland

Irrigated agriculture influences selenium levels in an endangered marsh bird

Selenium bioaccumulation in aquatic food webs poses risks to wildlife, particularly in wetlands receiving irrigation runoff. The Salton Sea, California’s largest lake, is primarily sustained by agricultural drainage. This drainage creates wetland habitat along the lakeshore that many bird species depend on, including the federally endangered Yuma Ridgway’s rail ( Rallus obsoletus yumanensis ). However, these marshes may pose an ecological trap – attracting rails despite high selenium exposure. We captured rails during the 2020–2023 breeding seasons and compared rail selenium levels within three types of marshes (fed with irrigation runoff, Colorado River water, or groundwater). We collected blood, breast feathers, and head feathers of rails in all three water sources for selenium comparisons. We tagged adult rails with GPS transmitters to locate nests and foraging locations where we collected eggshells, unhatched eggs, and prey. We assessed selenium exposure by collecting multiple prey species commonly eaten by rails in all three water sources. Selenium concentrations varied among sampling locations. Selenium concentrations in most sample types were predominately influenced by water source and marsh inflow velocity (sometimes in combination with marsh size). Distance to inflow, however, did not influence selenium concentrations in any sample type. Selenium concentrations were highest in agricultural-fed marshes compared to river-fed and spring-fed marshes. Increased marsh inflow velocities resulted in lower selenium concentrations. Given the risk of an ecological trap, our results suggest that supplementing wetlands with Colorado River water could mitigate selenium bioaccumulation in Yuma Ridgway’s rails.

California

Broadband stochastic simulation of earthquake ground motions with multiple strong phases with an application to the 2023 Kahramanmaraş, Turkey (Türkiye), earthquake

Stochastic ground motion simulation models are often less accurate at lower frequencies than at higher frequencies when fitting recorded data unless supplemented by a deterministic forward directivity velocity pulse model. Moreover, time-modulated stochastic models, which adjust ground motion amplitudes over time, typically use functions that fail to capture multiple strong-motion phases. The February 2023 Turkey (Türkiye) earthquake exhibited diverse recordings, including near-fault and far-field motions with pulse-like and non-pulse-like characteristics, along with single and multiple strong-motion phases. To better represent such a diverse set of recordings, this study enhances a fully non-stationary site-based stochastic model without combining it with a deterministic model. Improvements include a new band-pass filter with upper- and lower-frequency limits, which refines the representation of the low-frequency content. Moreover, a time-modulating function that can represent energy arrival in multiple strong phases is introduced. The reference model’s parameters are identified by fitting to the energy content, zero-level crossings, and cumulative counts of positive-minima and negative-maxima of a target accelerogram. This fitting procedure is modified to address the increased number of parameters. These improvements broaden the reference model’s applicability while preserving its simplicity, a key aspect appealing to engineering practitioners. The improved model’s applicability is demonstrated by simulating a dataset from the February 2023 Türkiye earthquake, and the accuracy is tested using a pulse-like Next Generation Attenuation Relationships for Western United States dataset. Validations are performed based on total energy, zero-level crossings, Fourier amplitude spectrum, elastic response spectra, and peak ground motion parameters. Validations are performed schematically in the time and frequency domains and quantitatively using goodness-of-fit scores, various validation-metrics errors, and inter-period correlations. Overall, the improved stochastic model can effectively simulate a set of diverse ground motion recordings, including near-fault pulse-like records, records with multiple strong phases, and far-field motions across a broad frequency range.

Kahramanmaraş

Evaluation of stream capture related to groundwater pumping, middle Humboldt River Basin, Nevada

Historical, future, and potential stream capture from groundwater pumping in the middle Humboldt River Basin (MHRB), Nevada, is estimated using a calibrated numerical groundwater flow model. The model was developed to estimate (1) stream capture, which is the change in flux between the groundwater system and the Humboldt River and tributaries, and (2) change in streamflow, which is the change in streamflow estimated for the Imlay gage on the Humboldt River (U.S. Geological Survey streamgage 10333000). Historical stream capture for water years (WYs) 1961–2015 is estimated using recorded and estimated groundwater pumping during that period. Future (predictive) stream capture was based on historical stresses (WYs 1961–2015) using a scenario that simulated non-mine pumping from WY 2015 at a uniform rate for 100 years into the future. Potential stream capture throughout the middle Humboldt River Basin from groundwater pumping during varying durations of time are presented in a series of capture maps. Maps also are presented that show the potential to capture from groundwater evapotranspiration, as well as the storage changes for pumping duration of 100 years. Estimates of historical stream capture from the mainstem Humboldt River during the early 1960s are less than 400 acre-feet per year (acre-ft/yr) when groundwater withdrawals and pumping rates were relatively small compared to more recent times. In the late 1980s and early 1990s, groundwater withdrawals increased and estimated historical stream capture also increased from about 4,000 acre-ft/yr in the late 1980s and early 1990s to as much as 18,800 acre-feet (acre-ft) in WY 1998. In WY 2015, estimated historical stream capture declined to about 13,000 acre-ft because of decreasing groundwater withdrawals and lower streamflow during the drought of WYs 2012–15, resulting in less stream water available for capture. Stream capture was estimated for 100 years into the future based on WY 2015 non-mine pumping rates and mine-dewatering activity through WY 2015. Stream capture is forecast to increase to about 23,000 acre-ft/yr, and streamflow in the Humboldt River could decrease by as much as 19,000 acre-ft/yr. Pumping for mine-dewatering and the associated discharge of that water affects streamflow in the Humboldt River at Imlay, Nevada (U.S. Geological Survey streamgage 10333000). Historically, from WYs 1991 to 2015, streamflow was greater at Imlay gage during active mine-dewatering from mine-water discharge operations and increased by as much as 105,000 acre-ft in WY 1998. The increase was attributed mostly to the discharge of groundwater from mine-related dewatering operations directly into the mainstem Humboldt River or its tributaries, with some of this increase associated with return flows from discharge to rapid infiltration basins. Results indicate that streamflow at Imlay gage is expected to decrease by as much as 1,600 acre-ft/yr 30 years after mine-related pumping and discharge are discontinued. The streamflow reductions at the Imlay gage are expected to then decrease to around 500 acre-ft/yr, 100 years after mine-related pumping and discharge are discontinued. Potential capture maps were produced for pumping durations of 10, 25, 50, and 100 years. Capture map results indicate that areas of greater potential stream capture occur adjacent to the Humboldt River and for upstream tributaries areas north of the Humboldt River.

Nevada