USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Ecological Engineering Design”

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,703 records · Page 60Linked to original sources
USGS2025

Characterizing the scale of regional landslide triggering from storm hydrometeorology

Rainfall strongly affects landslide triggering; however, understanding how storm characteristics relate to the severity of landslides at the regional scale has thus far remained unclear, despite the societal benefits that would result from defining this relationship. As mapped landslide inventories typically cover a small region relative to a storm system, here we develop a dimensionless index for landslide-inducing rainfall, A * , based on extremes of modeled soil water relative to its local climatology. We calibrate A * using four landslide inventories, comprising over 11 000 individual landslides over four unique storm events, and find that a common threshold can be applied to estimate regional shallow-landslide-triggering potential across diverse climatic regimes in California (USA). We then use the spatial distribution of A * , along with topography, to calculate the landslide potential area (LPA) for nine landslide-inducing storm events over the past 20 years, and we test whether atmospheric metrics describing the strength of landfalling storms, such as integrated water vapor transport, correlate with the magnitude of hazardous landslide-inducing rainfall. We find that although the events with the largest LPA do occur during exceptional atmospheric river (AR) storms, the strength of landfalling atmospheric rivers does not scale neatly with landslide potential area, and even exceptionally strong ARs may yield minimal landslide impacts. Other factors, such as antecedent soil moisture driven by storm frequency and mesoscale precipitation features within storms, are instead more likely to dictate the patterns of landslide-generating rainfall throughout the state.

California↗
USGS2024

Status and understanding of groundwater quality in the Mojave Basin Domestic-Supply Aquifer study unit, 2018—California GAMA Priority Basin Project

Groundwater quality in the western part of the Mojave Desert in San Bernardino County, California, was investigated in 2018 as part of the California State Water Resources Control Board Groundwater Ambient Monitoring and Assessment Program Priority Basin Project. The Mojave Basin Domestic-Supply Aquifer study unit (MOBS) region was divided into two study areas—floodplain and regional—to assess differences between the two major aquifers used for drinking water supply in the area. This assessment characterized the quality of ambient groundwater and not the quality of treated drinking water. The study included three components: (1) a status assessment, which characterized the quality of groundwater resources used for domestic drinking-water supply in the floodplain and regional study areas; (2) a brief understanding assessment, which evaluated factors that could potentially affect the quality of groundwater used by domestic wells in the region; and (3) a comparative assessment between the groundwater resources used by domestic wells and public-supply wells in the two study areas. The domestic-well assessment was based on data collected by the U.S. Geological Survey from 48 domestic wells in January–May 2018. The public-supply assessment was based on data for samples from 322 public-supply wells in 2008–18, either collected by the U.S. Geological Survey or compiled from the California State Water Resources Control Boards Division of Drinking Water publicly available database. Concentrations of water-quality constituents in ambient groundwater were compared to regulatory and non-regulatory benchmarks typically used by the State of California and Federal agencies as health-based or aesthetic standards for public drinking water. Relative concentrations, defined as the measured concentration divided by the benchmark concentration, were classified as high (greater than 1.0), moderate (greater than 0.5 for inorganic constituents or 0.1 for organic and special-interest constituents, and not high), or low (concentrations lower than moderate). The floodplain and regional study areas were divided into 15 and 35 grid cells, respectively, and grid-based methods were used to compute the areal proportions of the two study areas with high, moderate, or low relative concentrations of individual constituents and classes of constituents. For the domestic-supply assessment, one or more inorganic constituents with health-based benchmarks were detected at high relative concentrations in 58 percent of the regional study area and 13 percent of the floodplain study area. The inorganic constituents with health-based benchmarks detected at high relative concentrations in the regional study area were arsenic, chromium and hexavalent chromium, fluoride, adjusted gross alpha particle activity, uranium, molybdenum, strontium, and nitrate; only arsenic was detected at high relative concentrations in the floodplain study area. One or more inorganic constituents with secondary maximum contaminant level benchmarks were detected at high concentrations in 15 and 6.7 percent of the regional and floodplain study areas, respectively. The constituents detected at high relative concentrations in the regional study area were total dissolved solids, chloride, sulfate, and iron; only total dissolved solids and sulfate were detected at high relative concentrations in the floodplain study area. Organic constituents were not detected at moderate or high relative concentrations in either the regional or floodplain study areas. Volatile organic compounds were detected at low relative concentrations in 21 and 27 percent of the regional and floodplain study areas, respectively, and pesticides were detected at low relative concentrations in 9.1 and 20 percent of the regional and floodplain study areas, respectively. The only individual organic constituent detected in more than 10 percent of either study area was the trihalomethane trichloromethane. Total coliform bacteria were detected in 15 and 27 percent of the grid wells in the regional and floodplain study areas, respectively. The greater prevalence of high relative concentrations of many inorganic constituents in the regional study area compared to the floodplain area likely indicates the greater diversity of geologic material at depth in aquifer material and generally finer-grained alluvium compared to the floodplain study area combined with generally older groundwater that has had more contact time with aquifer materials. In general, trace element concentrations (1) increased with increasing groundwater age, (2) increased with distance from recharge sources in the mountains, and (3) increased with closer proximity to some types of geological units. In general, groundwater from domestic wells in the floodplain study area is young, with most samples containing a component of modern groundwater based on tritium and unadjusted carbon-14 activities, whereas groundwater from domestic wells in the regional study area generally is old, with most samples having unadjusted carbon-14 ages of 5,000–40,000 years. Public-supply wells in MOBS generally were deeper than domestic wells and presumably are in contact with older, more weathered alluvium that may have more mobile trace elements, such as arsenic or uranium. However, only 26 percent of the public-supply regional study area had high relative concentrations of inorganic constituents, compared to 58 percent for the domestic regional study area. The percentages of the public-supply and domestic floodplain study areas with high relative concentrations of inorganic constituents were 11 and 13 percent, respectively. The ages of groundwater used by public-supply and domestic wells in each study area were similar, which was not expected given the greater depth of the public-supply wells. Three potential factors may contribute to these results: (1) greater spatial footprint of domestic well network, which may result in domestic wells pumping groundwater from fractured bedrock or mineralized areas not used by public-supply wells; (2) greater pumping rates in public-supply wells, resulting in more water being withdrawn from coarse-grained, heterogeneous alluvium than finer-grained layers, which may have higher concentrations of (or more mobile) inorganic constituents; and (3) a greater degree of well management with public-supply wells, which may include pausing use of or decommissioning wells if treating or blending water is not feasible to lower constituent concentrations.

California↗
USGS2024

Rainfall intensification amplifies exposure of American Southwest to conditions that trigger postfire debris flows

Short-duration, high-intensity rainfall can initiate deadly and destructive debris flows after wildfire. Methods to estimate the conditions that can trigger debris flows exist and guidance to determine how often those thresholds will be exceeded under the present climate are available. However, the limited spatiotemporal resolution of climate models has hampered efforts to characterize how rainfall intensification driven by global warming may affect debris-flow hazards. We use novel, dynamically downscaled (3.75-km), convection-permitting simulations of short-duration (15-min) rainfall to evaluate threshold exceedance for late 21st-century climate scenarios in the American Southwest. We observe significant increases in the frequency and magnitude of exceedances for regions dominated by cool- and warm-season rainfall. We also observe an increased frequency of exceedance in regions where postfire debris flows have not been documented, and communities are unaccustomed to the hazard. Our findings can inform planning efforts to increase resiliency to debris flows under a changing climate.

Arizona, California, Nevada, New Mexico, Utah↗
USGS2026

Apatite (U-Th)/He thermochronology from Marsh Creek anticline reconciles Cenozoic and Holocene strain patterns and elucidates the Cenozoic canning displacement zone of Arctic Alaska

The northeastern Brooks Range is a vast fold-thrust belt that records the northernmost expression of Cenozoic deformation in Alaska, likely with complex dynamic linkages to distant tectonic elements. The principal focus of this manuscript is the Marsh Creek anticline (MCA), which is a large and recently active structure within the broader northeastern Brooks Range. The MCA is a complex, polyphase structure mostly buried beneath the Arctic coastal plain. We present apatite (U-Th)/He (AHe) dates from six outcrop samples, all near the structural crest of the anticline, as well as supporting detrital zircon U/Pb (ZUPb) age spectra from five outcrop samples. The AHe sample transect spans ~3 km of structural relief. The four shallowest AHe samples are not reset. However, the two deepest samples exhibit low intrasample variability among replicates and strong date versus grain radius correlations, suggesting post-depositional burial heating to within the AHe partial retention zone. Inverse thermal history modeling of these data resolves a protracted phase of exhumational cooling that began around 33–28 Ma. Detrital ZUPb age spectra from Paleogene strata support an interpretation of diverse inherited thermal histories recorded by the AHe data, and they are also consistent with extant reconstructions of Paleogene sediment dispersal. Placed into regional context, the results constrain a phase of increasing structural relief in the northeastern Brooks Range fold-thrust belt that began by the Eocene. The relief developed in a zone characterized by diminishing structural relief to the west called the Canning displacement zone. Although this zone was active since the Eocene, several structures including the studied part of the MCA were rapidly exhumed in the Oligocene. Lastly, we argue for kinematic compatibility between Cenozoic contraction and active strike-slip in the northeastern Brooks Range, with the temporal change in strain style likely reflecting the increase in structural relief and topography.

Alaska↗
USGS2000

Ground-water quality atlas of Oakland County, Michigan

The U.S. Geological Survey (USGS), in cooperation with Oakland County Health Division (OCHD), collected 140 water samples from 38 wells in Oakland County during 1998 to better understand ground-water quality. OCHD had observed temporal variations in concentrations of various constituents, so two additional sets of samples were collected to evaluate potential short-term variability related to sample collection procedures and long-term seasonal variability. Replicate samples from 28 wells were analyzed in the Michigan Department of Environmental Quality (MDEQ) Drinking Water Laboratory to compare MDEQ’s analytical results to those obtained from the USGS National Water Quality Laboratory. Several additional databases describing population, land use, water supply, soils, geology, and flows of ground water and surface water are presented in the first part of the report to assist in interpreting the water-quality data. Maps created from these databases are provided in the first portion of the report as an extension of the study-area description. The U.S. Environmental Protection Agency (USEPA) has established Maximum Contaminant Levels (MCL) and Secondary Maximum Contaminant Levels (SMCL) for which samples were analyzed in this study. Water from the 38 wells sampled by the USGS did not exceed the SMCL or MCL for sulfate, fluoride, or nitrite. However, water from 26 wells exceeded the SMCL for iron, water from 12 wells exceeded the SMCL for manganese, and water from 12 wells exceeded the SMCL for dissolved solids. Water from two wells exceeded the MCL for nitrate, although nitrate concentrations in water from most wells was below the detection limit. Water from seven wells exceeded the SMCL for chloride, and water from all wells contained detectable concentrations of chloride. Water from five wells exceeded the MCL for arsenic, and most of the wells sampled contained detectable concentrations of arsenic. These five wells were identified from previous MDEQ analyses to have elevated arsenic concentrations, and were sampled to obtain additional chemistry information. Replicate samples were collected from 26 of the 38 wells for analysis at the MDEQ Drinking-Water Laboratory to compare the results with the USGS National Water Quality Laboratory. The results of the replicate analyses indicate close agreement between the laboratories, with mean differences for nitrate, chloride, and arsenic of 0.10 milligrams per liter (mg/L) as nitrogen, 6.8 mg/L, and 0.0008 mg/L, respectively between the USGS and MDEQ analyses. Potential health effects associated with ingesting nitrate, chloride, and arsenic are provided with the water-quality data, along with references for further information.

Michigan↗
USGS2025

False positives in the identification of dynamic earthquake triggering

Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with global M ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial-temporal windows in Southern California tested in DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 , thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023, https://doi.org/10.1029/2023jb026487 ), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor-Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.

JGR Solid Earth↗
USGS2024

Layered intrusions in the Precambrian: Observations and perspectives

Layered intrusions are plutonic bodies of cumulates that form by the crystallization of mantle-derived melts. These intrusions are characterized by igneous layering distinguishable by shifts in mineralogy, texture, or composition. Layered intrusions have been fundamental to our understanding of igneous petrology; however, it is their status as important repositories of critical metals – such as platinum-group elements, chromium, and vanadium – that has predominantly driven associated research in recent decades. Many layered intrusions were emplaced during the Precambrian, predominantly at the margins of ancient cratons during intervals of supercontinent accretion and destruction. It appears that large, layered intrusions require rigid crust to ensure their preservation, and their geometry and layering is primarily controlled by the nature of melt emplacement. Layered intrusions are best investigated by integrating observations from various length-scales. At the macroscale, intrusion geometries can be discerned, and their presence understood in the context of the regional geology. At the mesoscale, the layering of an intrusion may be characterized, intrusion-host rock contact relationships studied, and the nature of stratiform mineral occurrences described. At the microscale, the mineralogy and texture of cumulate rocks and any mineralization are elucidated, particularly when novel microtextural and mineral chemical datasets are integrated. For example, here we demonstrate how mesoscale observations and microscale datasets can be combined to understand the petrogenesis of the perplexing snowball oiks outcrop located in the Upper Banded Series of the Stillwater Complex. Our data suggest that the orthopyroxene oikocrysts did not form in their present location, but rather formed in a dynamic magma chamber where crystals were transported either by convective currents or within crystal-rich slurries. Critical metals may be transported to the level of a nascent intrusion as dissolved components in the melt. Alternatively, ore minerals are entrained from elsewhere in a plumbing system, potentially facilitated by volatile-rich phases. There are many ore-forming processes propounded by researchers to occur at the level of emplacement; however, each must address the arrival of the ore mineral, its concentration of metals, and its accumulation into orebodies. In this contribution, several of these processes are described as well as our perspectives on the future of layered intrusion research.

Precambrian Research↗
USGS2025

Implications of physics-based M9 ground motions on liquefaction-induced damage in the Cascadia Subduction Zone: Looking forward and backward

Given the likelihood of future M9 Cascadia Subduction Zone (CSZ) earthquakes, various estimates of the resulting, regional ground motions have been made, including a suite of 30 physics-based simulations that reflect key modeling uncertainties. However, because the last CSZ interface rupture occurred in 1700 CE, the shaking expected in such an event is especially uncertain, as are the impacts to the built and living environments. Like other coseismic impacts, soil liquefaction poses a significant threat and must be considered by any scenario study used to inform planning and response, or to focus mitigation resources. Liquefaction is also notable for its potential to “ground truth” ground-motion estimates, given that its presence or absence in the geologic record can provide constraint on the intensities of shaking in past events. It is thus an important phenomenon looking both forward and backward. Accordingly, using recent physics-based simulations, this study (1) predicts liquefaction in M9 CSZ ruptures at 400 locations in Oregon, Washington, and British Columbia (BC) using an array of cone-penetration-test based models and (2) uses paleoliquefaction evidence at ten sites spanning from Southern Oregon to Vancouver, BC to constrain possible ground-motion intensities experienced in the 1700 CE earthquake. The forward predictions indicate that liquefaction in M9 events could be pervasive in the region and affect numerous population hubs, with the potential for damage across hundreds of square kilometers. The backward analyses suggest that 1700 CE ground-motion intensities may have been less than expected from M9 simulations in some northern portions of the CSZ (e.g. Seattle), given the paucity of 1700 CE liquefaction evidence in these areas. Ultimately, further discovery and analysis of CSZ paleoliquefaction, or lack thereof, will confirm or modify this possibility and the conclusions drawn herein.

Oregon, Washington↗
USGS2026

Lithium brine potential of North Dakota: Results from the 2025 sampling program

Lithium (Li), classified as a critical mineral and key in energy storage applications because of its high energy density, faces unprecedented demand. This has driven interest in unconventional sources of Li, especially brines from sedimentary basins. Accordingly, the spatial distribution of basinal brinehosted Li resources across North America is becoming better defined, including the Williston Basin underlying Saskatchewan, Canada. This has led to growing investment and exploration by industry throughout the 2020s; however, despite this interest, the Li potential of basinal brines of the Williston Basin in North Dakota remain underexplored. This work presents the results from the 2025 joint brine sampling campaign between the North Dakota Geological Survey, the United States Geological Survey (USGS), and The University of Texas at Austin. Twenty-seven fluid samples were collected from producing oil and water source wells from eight stratigraphic intervals across western North Dakota and analyzed for major cation and anion concentrations. Similar to southeastern Saskatchewan, the highest Li concentrations were identified in the Devonian Duperow Formation, with concentrations up to 179 mg/L being observed in a well commingled between the Duperow and Red River Formations. Some samples from the overlying Birdbear Formation exceeded 75 mg/L, whereas Li concentrations were lower in the Madison Group (33-44 mg/L) and the Three Forks (50 mg/L), Dawson Bay (43 mg/L), Winnipegosis (23 mg/L), and Red River (40.0 mg/L) Formations. This work indicates Li resource potential exists across the Williston Basin and advances our understanding of the Li distribution in basinal brines which could provide a new domestic source of this critical mineral.

North Dakota↗
USGS2024

Agricultural return flow dynamics on a reach of the East River, Colorado, as assessed by mass balance

The U.S. Geological Survey, in cooperation with the Upper Gunnison River Water Conservancy District, studied historical streamflow in a reach of the East River, Colorado, to gain a preliminary understanding of return flow dynamics. Return flow is agricultural irrigation water that is not consumed by evapotranspiration and instead reaches streams by surface and subsurface flow paths. The study reach had a contributing area of 50 square miles and contained 5.23 square miles of pastures irrigated with water diverted from the East River and its tributaries. By comparing upstream inflows to downstream outflows, the net water balance of the study reach from 1994 to 2023 was assessed. Two general hydrologic conditions for the study reach were identified. One hydrologic condition was characterized by a net loss or consumption of water, termed here as general deficit. This general deficit condition extended about 16 years, from 1997 to 2012. During general deficit years, there was usually a notable net loss of streamflow from April through July, and a small net gain, possibly related to return flows, occurred in August about 75 days after the minimums for losses. The second hydrologic condition was characterized by a net gain of water, termed here as general surplus. This second condition extended about 10 years, from 2014 to 2023. During general surplus years, two separate transitions from net loss to net gain commonly occurred during June through August. Losses during general surplus years were smaller than losses during general deficit years, the respective gains were larger, and times between losses and gains were about 18 and 22 days. Differences between the two hydrologic conditions could reflect interactions among irrigation water, available capacity to store additional shallow groundwater, and streamflow. However, deciphering the causes for the shifts between the two general hydrologic conditions was beyond the scope of this report.

Colorado↗
USGS2026

Distributed volcanic fields—Patterns, products, processes, and plumbing

Distributed volcanic fields—also referred to as “monogenetic fields” or simply “volcanic fields”—are found on every continent, on the seafloors, and in every tectonic setting on Earth. Yet, our understanding of what constitutes a distributed volcanic field, and our knowledge of the common products, spatiotemporal patterns in activity, and plumbing systems of these fields, are commonly poor and (or) debated. Here we provide an overview of the current understanding of distributed volcanic fields (DVFs) through review and compilation of published data, with an emphasis on understanding the characteristics of fields that are essential to prepare for future hazardous volcanic activity. Distributed fields can erupt restricted to variable magma compositions. Activity in these fields produces nearly all types of volcanoes, including scoria and spatter cones, maars, low shields, lava fields, stratovolcanoes, and calderas, but monogenetic volcanoes are the defining characteristic of DVFs. The number of vents and lifespans of DVFs vary greatly; vent numbers range from 2 to greater than 1,000, and fields can remain active for 150,000 to 10 million years. Some characteristics of DVFs appear field dependent. Although the volumes of single monogenetic eruptions within DVFs globally are variable (from less than 0.1 to a few cubic kilometers), each field may have a characteristic range in eruption volumes. Eruptions within distributed volcanic fields can have a random distribution and timing, but spatial and (or) temporal vent clustering and alignment are common, and eruptions may be time or volume predictable. Average eruption recurrence intervals in DVFs span over an order of magnitude, from less than 3 to greater than 100 thousand years. Exposures of the shallow portions of older systems reveal that many magma batches end their ascent before reaching the surface, with some dikes stopping within only meters of the surface. Petrologic insights into the deeper plumbing systems of DVF volcanoes indicate variable depths of magma crystallization and (or) stalling, and that magma transport and storage systems may be unique to a particular DVF and its underlying crustal structure. Both seismic and petrologic data for mafic DVF eruptions corroborate intrusions of magma into storage systems in the months to years before eruption, as well as variable ascent rates through the crust (days to weeks). Few distributed fields are closely monitored, as they erupt infrequently, but these observations imply that the amount of time available to install monitoring instrumentation after the first felt seismicity may be short, as precursory activity is frequently on the order of days to months for modern DVF eruptions. Satellite-based geodesy may provide useful warning in such cases. This chapter concludes with outstanding questions and areas for future research that will aid in our understanding and preparedness for future DVF eruptions.

Professional Paper↗
USGS2026

Fully contained laboratory earthquakes: The effect of asperity aspect ratio and free surfaces

Corner frequency ( f c ) and seismic moment ( M 0 ) are key parameters derived from seismic signals that are used to characterize earthquake stress drop, rupture area, and slip. These parameters are also affected by fault geometry and boundary conditions. However, the systematic study of these effects in laboratory settings has been challenging. This study presents laboratory earthquake experiments that examine how rupture dynamics are influenced by (a) the aspect ratio of rectangular PMMA velocity-weakening (VW) asperities surrounded by the Teflon velocity-strengthening (VS) patches, and (b) whether the sides of a VW asperity are confined with VS patches or are free surfaces. We found that increasing confinement by reducing free surfaces or increasing the VW asperity aspect ratio stabilizes fault slip, so that higher normal stress is required to transition from aseismic to seismic slip. Increased confinement and high aspect ratios also reduced M 0 and increased f c , both of which were determined from the radiated seismic waves. M 0 and f c were primarily controlled by the shorter dimension of the VW asperity. Analysis of high-frequency acoustic emission signals revealed that ruptures on high-aspect-ratio VW asperities propagated more unidirectionally, whereas ruptures on square VW asperities were more complex. Further, the high-aspect-ratio asperities were more likely to be eroded by surrounding VS regions while low-aspect-ratio asperities were more likely to rupture into the VS surroundings. These results demonstrate that both the confinement from surrounding stable areas and the geometry of the seismogenic patch can affect rupture nucleation, propagation, and seismic source characteristics.

JGR Solid Earth↗
USGS2026

Determination of baseline water-quality conditions using continuous specific conductance and discrete measurements of total dissolved solids and selected major ions, Arkansas River Basin, Colorado

From 2022 through 2024, the U.S. Geological Survey (USGS) in cooperation with the Regional Resources Planning Group (Southeastern Colorado Water Conservancy District, Colorado Springs Utilities, Pueblo Water, Lower Arkansas Valley Water Conservancy District, Aurora Water, and Upper Arkansas Water Conservancy District) established a basinwide water-quality monitoring network in the Arkansas River Basin in Colorado. The study leveraged ongoing continuous measurements of streamflow and specific conductance (SC) at 13 sites along the main stem Arkansas River and Fountain Creek from the USGS 07108600 Arkansas River at Granite, Colo., (Ark Granite) site to the USGS 07137500 Arkansas River near Coolidge, Kansas, (Ark Coolidge) site. Baseline water-quality conditions were established using data from October 1989 through September 2022 to help identify changes in water quality that could result from changes in land and water use. Total dissolved solids (TDS) concentration was the main indicator of general water quality in this study because it can be accurately estimated from SC and is sensitive to changes in streamflow, water operations, and source contributions. Median TDS concentrations in Arkansas River samples (October 1989 through December 2023) ranged from 75 milligrams per liter (mg/L) at Ark Granite near the headwaters to 3,660 mg/L at Ark Coolidge near the Colorado–Kansas State line. Distinct seasonal differences in median streamflow and TDS were apparent during the winter and summer months. Seasonal Kendall tests showed no significant trends in monthly mean streamflow (October 1995 through September 2022) at Arkansas River sites upstream from Pueblo Reservoir. Upward trends in TDS concentrations and loads occurred only at Ark Granite. Downward trends in streamflow, TDS concentrations, and TDS loads occurred at Arkansas River sites from USGS 07109500 Arkansas River near Avondale, Colo., to Ark Coolidge. Downward trends in streamflow and SC (surrogate for TDS) were more prevalent during winter base flow months and may indicate changes in groundwater quantity, and possibly groundwater quality, potentially affecting the river water-quality condition.

Colorado↗
USGS2026

Groundwater source water assessment area delineation of public supply wells in the Treasure Valley, western Idaho

The Idaho Department of Environmental Quality (IDEQ) is responsible for assessing the source water for public water supply systems including the delineation of source water assessment areas for wells, which are defined as the subsurface areas contributing water to a well during a period of time. Part of the IDEQ assessment plan is to update source water assessments over time. The U.S. Geological Survey (USGS) developed a new groundwater-flow model for part of the western Snake River Plain aquifer in southwestern Idaho. This groundwater- flow model provided a new opportunity for IDEQ to update source water assessment for public water supply wells within the model boundary. The USGS, in cooperation with IDEQ, used the USGS Treasure Valley Groundwater Flow Model along with a MODPATH particle tracking post-processing program to reassess public-water supply source water assessment areas. The USGS created an automated script that requests a well location, depth, and pumpage rate, and then incorporates that information into the Treasure Valley Groundwater Flow Model combined with the MODPATH model to create 3-, 6-, and 10-year time of travel zones that are output to a geographic information system shapefile. The TVSWA_Delineator script provides a fast, transparent tool for IDEQ to reassess or create new public-water supply source water assessment area delineations based on new and updated information and is available as a Python script and a Windows executable with the compiled script and all necessary components. TVSWA_Delineator utilizes a quadtree approach to rediscretize the Treasure Valley Groundwater Flow Model grid around the pumping well location from 1-mile cell lengths telescoped down to 1/16 mile. This rediscretization increases the source water assessment area delineation resolution near the pumping well and minimizes the effect of weak sinks in the MODPATH model. The stress inputs for each month of the final 10 years (2006 to 2015) of the Treasure Valley Groundwater Flow model historical simulation were averaged for the particle tracking analysis. The MODPATH model used backwards tracking of particles released at each stress period to create the source water assessment area delineations. This report documents the creation of the TVSWA_Delineator Python script, including adaptations made to the Treasure Valley Groundwater Flow Model, creation of a companion MODPATH model, and post-processing results to provide source water assessment areas.

Idaho↗
USGS2026

The geologic history of the Chehalis Forearc Basin, Washington State, USA

The Chehalis basin is located between the Cascade arc and the Coast Range in southwest Washington State. It consists of sedimentary and volcanic rocks deposited throughout the Cenozoic and is underlain by the Siletzia terrane, a thick, large igneous province accreted in the Eocene. Here, we constrain evolution of the Chehalis basin depocenter using isochore maps derived from isostatic gravity anomalies, borehole data, and stratigraphy for several time periods: the Eocene (47.6–36.8 Ma), the Oligocene (36.8–20 Ma), and the Neogene (20–0 Ma). Our results suggest that local subsidence is driven by deformation on faults that bound and intersect the basin. We see northward depocenter migration, interpreted as a shift in deformation from the northwest striking Cedar Creek fault to the west striking Doty fault. We interpret these data in terms of the long-term north-south shortening and clockwise rotation of the Cascadia forearc. During the Eocene, the Cedar Creek fault was preferentially aligned with north-south shortening, but became less active when clockwise rotation brought the Doty fault into east-west alignment with the stress field. Sediment accumulation rates decreased from 196 m/Myr in the Paleogene (∼40–20 Ma) to 27 m/Myr in the Neogene (∼20–0 Ma) as Miocene uplift of the Coast Range limited accommodation space and shifted the depositional environment from marine to fluvial. Our results are consistent with the geologic evolution of the Portland and Tualatin basins to the south and reveal a depositional and structural history uniquely shaped by clockwise rotation of the Cascadia forearc.

Washington↗
USGS2024

Shallow faulting and folding beneath south‐central Seattle, Washington State, from land‐based high‐resolution seismic‐reflection imaging

The geologic framework of the Seattle fault zone (SFZ) has been extensively studied, but the structure and fault strand locations in the central portion of the fault zone through the city of Seattle have remained controversial. Much of what is known about the SFZ has come from light detection and ranging (lidar)‐topographic surveys and paleoseismic investigations of fault scarps primarily west of Puget Sound, regional gravity and aeromagnetic modeling, and multiscale marine seismic imaging in waters both west and east of Seattle. We analyze ∼24 km of land‐based P ‐wave seismic‐reflection data that fill in a critical gap in our understanding of the SFZ beneath the urban areas of West Seattle, south‐central Seattle, and Mercer Island. These data image deformed strata in the upper 1 km, including upwarped Tertiary rock and younger sediments. Collectively, these data provide evidence for multiple Quaternary‐active thrust faults, back thrusts, and sub‐basins within the SFZ beneath the city of Seattle. The results indicate that multiple and potentially active back thrusts in the upper ∼500 m extend across the length of the SFZ and the entire urban corridor that may be analogous to those on Bainbridge Island west of Puget Sound.

Washington↗
USGS1977

Exploratory water well, St. George Island, Florida

A 1,026-foot well was drilled into the Floridan aquifer on St. George Island in search of a supply of water for residential development in 1971. This is the only well on the island on record to have been drilled below 200 feet for water. The well was abandoned; water flowing from the well was saline. It had a specific conductance of 7,500 micromhos per centimeter at 25° Celsius. The specific conductances of samples taken at drilling depths 217-1,026 feet ranged from 3,140 to 39,600, with the maximum value at 470 ft. A composite geophysical and lithologic log of the well shows spontaneous potential, single-point resistance, gamma-ray, neutron, gamma-gamma, and caliper logs. Well-site notes, a written description of cuttings samples, a graph showing the relation of specific conductance of water samples to the depth of the well at time of sampling, and four chemical analyses of the water are included in the report.

Florida↗
USGS2026

Baseflow and snowmelt sustained streamflow in the Upper Colorado River Basin, 1986-2020

The Upper Colorado River Basin (UCRB) faces substantial water availability limitations. Although most streamflow originates as snowmelt, the partitioning of snowmelt between surface runoff and groundwater recharge and subsequent groundwater discharge to streams is highly uncertain. On average, over half of the streamflow in the UCRB is estimated to originate from groundwater discharge to streams, highlighting the importance of baseflow in sustaining surface water. However, the historical patterns of baseflow and streamflow, along with their variability over space and time and their specific sources, remain unknown at the basin scale. This study addresses those gaps by characterizing the sources and transport pathways of both baseflow and streamflow in the UCRB at a seasonal timestep from 1986 to 2020, including the lagged delivery of subsurface water to streams beyond the current season, using coupled models of baseflow and streamflow. Between 1986 and 2020, on average 63% of UCRB streamflow originated from baseflow. About half of this baseflow took longer than one season to reach streams, and outside the snowmelt season, baseflow was the dominant source of streamflow. Snowmelt was a key source of both baseflow and streamflow. Current season snowmelt contributed 33% of streamflow via runoff, and 22% of the 29% of streamflow that originated as current season baseflow via subsurface flow to streams. Over the study period, baseflow index (BFI) declined in headwaters and increased at mid-elevations. Springtime increases in BFI demonstrate the increasingly important role baseflow plays in water supply. Identifying the sources, locations, and timing of water that contributed to the UCRB outlet can inform management of water resources in the basin.

Arizona, Colorado, New Mexico, Utah, Wyoming↗