Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of the Nevada Water Resources Association”

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

Tapwater exposures, residential risk, and mitigation in a PFAS-impacted-groundwater community

Tapwater (TW) safety and sustainability are priorities in the United States. Per/polyfluoroalkyl substance(s) (PFAS) contamination is a growing public-health concern due to prolific use, widespread TW exposures, and mounting human-health concerns. Historically-rural, actively-urbanizing communities that rely on surficial-aquifer private wells incur elevated risks of unrecognized TW exposures, including PFAS, due to limited private-well monitoring and contaminant-source proliferation in urbanizing landscapes. Here, a broad-analytical-scope TW-assessment was conducted in a hydrologically-vulnerable, Mississippi River alluvial-island community, where PFAS contamination of the shallow-alluvial drinking-water aquifer has been documented, but more comprehensive contaminant characterization to inform decision-making is currently lacking. In 2021, we analyzed 510 organics, 34 inorganics, and 3 microbial groups in 11 residential and community locations to assess (1) TW risks beyond recognized PFAS issues, (2) day-to-day and year-to-year risk variability, and (3) suitability of the underlying sandstone aquifer as an alternative source to mitigate TW-PFAS exposures. Seventy-six organics and 25 inorganics were detected. Potential human-health risks of detected TW exposures were explored based on cumulative benchmark-based toxicity quotients ( ∑ TQ ). Elevated risks ( ∑ TQ ≥ 1) from organic and inorganic contaminants were observed in all alluvial-aquifer-sourced synoptic samples but not in sandstone-aquifer-sourced samples. Repeated sampling at 3 sites over 52–55 h indicated limited variability in risk over the short-term. Comparable PFAS-specific ∑ TQ for spatial-synoptic, short-term (3 days) temporal, and long-term (3 years quarterly) temporal samples indicated that synoptic results provided useful insight into the risks of TW-PFAS exposures at French Island over the long-term. No PFAS detections in sandstone-aquifer-sourced samples over a 3 year period indicated no PFAS-associated risk and supported the sandstone aquifer as an alternative drinking-water source to mitigate community TW-PFAS exposures. This study illustrated the importance of expanded contaminant monitoring of private-well TW, beyond known concerns (in this case, PFAS), to reduce the risks of a range of unrecognized contaminant exposures.

Wisconsin↗
USGS2026

Non-destructive sampling of larval amphibians shows tail tissues reflect whole-body methylmercury concentrations and trophic-related differences in bioaccumulation

Mercury (Hg), and specifically methylmercury (MeHg), is a contaminant of global concern to humans and wildlife, but there is still limited understanding of its effects on many taxa. Amphibians are closely associated with water, where Hg is converted to MeHg, and high concentrations of MeHg can reduce survival of amphibians in the wild. However, limited non-lethal proxies exist for estimating MeHg bioaccumulation in amphibian larvae, which could be used to reduce destructive sampling and allow sampling of imperiled species. We evaluated whether tail clips were indicative of whole-body MeHg concentrations of larval amphibians, whether dragonflies were an effective bioindicator of Hg concentrations of larval amphibians, and compared MeHg bioaccumulation in predatory caudate (salamanders, newts) and non-predatory anuran (frogs, toads) larvae. Tail-clip and whole-body MeHg concentrations were strongly correlated, especially for caudates. There was high variability among sites and species, but caudates had higher MeHg than anurans, and MeHg concentrations of caudates and anurans were strongly correlated. Dragonflies were a weak indicator of MeHg for anurans but strong indicators for larval caudates, likely because caudates and dragonflies are carnivorous and occupy similar trophic positions. Our study revealed that sampling tails can be an effective non-destructive index of whole-body MeHg for larval amphibians and demonstrates differences in MeHg bioaccumulation related to trophic position in larval amphibians and dragonflies. All post-embryo life stages for a broad suite of amphibian species can now be sampled non-lethally for MeHg, allowing for a broader examination of population-level effects and focus on imperiled species that cannot be sampled destructively.

conterminous United States↗
USGS2026

Riverscape heterogeneity shapes population diversity for a migratory fish

Habitat patch dynamics can scale up to influence population demography and diversity with implications for resilience to environmental stochasticity. But how the spatial arrangement and size of habitat patches interact with other components of habitat heterogeneity to shape population diversity at larger spatial scales is not well understood. For riverine fishes, there is increasing evidence that tributary streams provide critical demographic support to main stem rivers. However, the extent to which main stem rivers rely on demographic contributions from tributaries, and the factors underlying this dependence, have not been assessed. Here, we used genetic stock identification to evaluate the effect of tributaries on population diversity of Yellowstone cutthroat trout ( Oncorhynchus virginalis bouvieri ) occupying the main stem Snake River, Wyoming, USA. We found that the main stem relied almost entirely on tributaries for demographic support, but main stem composition varied spatially among river sections. Distance between habitat patches, catchment area, and groundwater availability acted in concert to determine the contribution of specific tributaries to the main stem, but contributions were ultimately modulated by habitat connectivity. We also found evidence for multi-scale spatial structure in tributary contributions, providing insight into untested drivers of main stem river population diversity. Our results demonstrate how spatially discrete and distributed riverscape attributes influence population diversity at broader spatial scales, illustrating how ecosystem resilience emerges from the dynamic, two-way exchange of individuals and energy across habitat networks. Management plans for large rivers that address the ecological contributions of tributaries may be needed to achieve optimal outcomes. Similarly, conservation strategies that exclusively focus on headwater streams may fail to capture the broader habitat requirements necessary to maintain robust cold-water fish populations and associated recreational fisheries, particularly under global environmental change.

Wyoming↗
USGS2026

Bottom-up characterization of geologic methane emissions in the San Juan Basin in the southwestern USA

Methane is a potent greenhouse gas that plays an important role in atmospheric chemistry and global warming. The current global methane budget has large uncertainties, and a better understanding of the budget would help to guide strategies for reducing anthropogenic emissions to fight climate change. Natural geologic methane emissions are a particularly poorly constrained source, with top-down estimates from 14 C in ice cores suggesting much lower geologic emissions than bottom-up scaling of direct flux measurements. Our study aims to contribute to resolving this discrepancy through improved bottom-up characterization of geologic methane seepage in the San Juan Basin in southwestern Colorado and northwestern New Mexico, USA. We performed 983 new flux chamber measurements in this basin during summer 2022 and winter 2023 field campaigns. Our results, in combination with prior measurements, suggest that natural seepage in the San Juan Basin only occurs on or near the Fruitland coal outcrop. Specifically, our new measurements confirm previous measurements of seepage along the northwestern exposure of the Fruitland outcrop in Colorado (a known hydrodynamic overpressure region) and for the first time, identified seepage locations along the southernmost Fruitland outcrop exposure in New Mexico, in association with a coal cleat and a fault. Overall, seepage along the Fruitland coal outcrop is heterogeneously distributed, with both positive and negative (interpreted as microbial soil sink) methane fluxes. Features that are hypothesized to be predictive of seepage (e.g., faults) were not associated with positive methane fluxes in areas outside of the Fruitland outcrop. Our best estimate for total geologic methane seepage in the San Juan Basin from spatial interpolation and statistical upscaling is approximately 0.14 Tg CH 4 /yr, with a range from 0.029 to 0.48 Tg CH 4 /yr. This best-estimate value is lower than a previous bottom-up estimate from a gridded seepage inventory, but higher than a previous top-down estimate.

Colorado, New Mexico↗
USGS2026

Reassessing water availability in the Nile River Basin using satellite-based irrigation water use accounting

Existing estimates of irrigation areas and water use in the Nile River Basin vary widely and are based on inconsistent data and methodologies. This study leverages advances in remote sensing to provide an up-to-date and consistent basin-wide reassessment of water use in the Nile River Basin. Land cover and actual evapotranspiration data from 2013 to 2022 were used to quantify irrigated areas, irrigation water use, and the average naturalized water yield. Across the basin, 7.1 million hectares of irrigated land consume an average of 72.5 ± 2.6 billion cubic meters (BCM) of blue water annually. Egypt and Sudan together account for 92% of the irrigated area and 97% of the irrigation water use. The basin’s average naturalized yield, evaluated at Dongola station, is determined to be between 108 and 114.4 ± 3.2 BCM. Although conservative, this updated analysis indicates the basin’s average naturalized yield is 14%–20% higher than the previous naturalized yield estimate (95 BCM) and 29%–36% higher than the commonly cited observed flow of 84 BCM, which is frequently treated as the total basin-wide available water. Importantly, the study underscores that the 84 BCM represents flow at the Aswan gauging station after upstream consumptive uses, not the total water available in the basin under naturalized conditions. This distinction is critical for accurate water accounting, planning, and governance in this water-scarce basin. This study demonstrates the value of open-source remote sensing resources in data-scarce regions while emphasizing the need for region specific validation and bias correction for improved accuracy.

Nile River Basin↗
USGS2025

Data gap analysis for estimation of agricultural return flows in the Upper Gunnison River Basin, Colorado

The Gunnison River and many tributaries in the Upper Gunnison River Basin provide water to irrigate agricultural crops. The application of irrigation water can recharge some aquifers locally by water percolating below the root zone and eventually flowing back to the stream or river through the subsurface. Diverting surface water for irrigation reduces streamflow during the irrigation season but can provide temporary storage of water and supplement streamflow after the snowmelt runoff season. Understanding the timing and quantity of agricultural return flows could help resource managers make informed decisions and adapt to potential changes in water management and availability that could affect irrigation practices. In 2024, the U.S. Geological Survey, in cooperation with the Upper Gunnison River Water Conservancy District, began a study to characterize agricultural return flows in the Upper Gunnison River Basin by using endmember mixing analysis and developing a groundwater model. Both approaches require data from multiple sources, but data gaps exist in the East River study reach and other reaches of interest (Ohio Creek, Tomichi Creek, and Cochetopa Creek). The East River Basin, which is the initial focus of the study, has fewer data gaps than the other basins. Data gaps could be addressed by installing additional surface water and groundwater monitoring sites, making regular streamflow measurements on tributaries, and completing tests to characterize local aquifer properties.

Colorado↗
USGS2024

New U-Pb geochronology and geochemistry of Paleozoic metaigneous rocks from western Yukon and eastern Alaska, cross-border synthesis, and implications for tectonic models

The tectonic evolution of and relation between the Yukon-Tanana terrane and the Lake George assemblage, as well as other associated tectonic assemblages in western Yukon and eastern Alaska, have been debated for decades. The Yukon-Tanana terrane is widely considered to be an allochthonous rifted fragment derived from the Laurentian continental margin, whereas the Lake George assemblage and associated assemblages are currently interpreted to be part of the parautochthonous continental margin of western North America (Laurentia). To address these topics, we present 40 new U-Pb zircon ages and 20 new whole-rock geochemical analyses. We incorporate these data into a new compilation of available geological mapping for a large area that straddles the Alaska-Yukon border, together with 34 previously published U-Pb age determinations and an extensive geochemical database of metaigneous rocks from Late Devonian to Early Mississippian and middle to late Permian assemblages in this area. Magmatism in the Lake George assemblage and related assemblages occurred in two pulses from about 371 to 360 and from about 358 to 347 million years ago (Ma); geochemical discrimination diagrams indicate a large crustal component, possibly indicative of arc magmatism, for felsic metaigneous rocks and a range of tectonic environments for mafic rocks. Magmatism in the Fortymile River and related assemblages, and parts of the Nasina assemblage—all parts of the Yukon-Tanana terrane—are mainly Early Mississippian and span a crystallization age range from about 361 to 343 Ma; geochemical discrimination diagrams for these rocks indicate primarily arc geochemical signatures for both mafic and felsic rocks. Middle to late Permian crystallization ages (about 261–253 Ma) are indicated for felsic metaigneous rocks in the Klondike assemblage and some of the felsic metaigneous rocks in the Nasina assemblage. Based on our mapping, we propose the existence of a possible unconformity between the Mississippian and Permian felsic metavolcanic rocks within the Nasina assemblage that is marked by sporadic occurrences of stretched-pebble conglomerate. Our combined database supports the well-established model of a magmatic arc comprising the Fortymile River and Finlayson assemblages of the rifted Yukon-Tanana terrane continental fragment on which a middle to late Permian arc (Klondike assemblage) was later built. The assemblages of the Yukon-Tanana terrane were subsequently intruded by Late Triassic to Early Jurassic granitoids, presumably during reaccretion of the Yukon-Tanana terrane to the continental margin. Permian and Late Triassic to Early Jurassic intrusions have not been mapped in the now structurally lower plate Lake George assemblage; their absence is one of the lines of evidence that have been used to support the parautochthonous, rather than allochthonous, origin of the Lake George assemblage and related assemblages. Our new data, together with previously published ranges of igneous crystallization ages and geochemical tectonic signatures of the Late Devonian to Early Mississippian magmatic rocks in the Lake George assemblage and associated assemblages and in the Fortymile River, Nasina, and correlated assemblages of the Yukon-Tanana terrane, indicate that the currently accepted interpretation of the Lake George assemblage and associated rocks being part of parauthochthonous North America is not the only possible interpretation of this tectonic entity. Approximately half of the dated intrusive rocks in the Lake George assemblage are contemporaneous with the metaigneous rocks of the Yukon-Tanana terrane arc (<361 Ma). We speculate that our approximately 361 Ma U-Pb age for quartz syenite in part of the North American continental margin in south-central Yukon defines the beginning of rifting of the Laurentian margin. Although the currently favored model of prolonged middle Paleozoic subduction and extension in both the Yukon-Tanana terrane and parautochthonous North America allows for simultaneous middle Paleozoic magmatism on both sides of the Slide Mountain Ocean, we now propose an alternative hypothesis in which the Lake George assemblage represents a deeper part of the rifted Yukon-Tanana terrane arc. If this is the case, the absence of Permian and Late Triassic to Early Jurassic arc rocks in the Lake George assemblage could be explained either by the arcs of these ages not being wide enough to have affected the Lake George assemblage or by tectonic displacement of these arc rocks away from the Lake George assemblage. Our approximately 259 Ma U-Pb zircon age and geochemical analyses of metarhyolite in the Seventymile terrane in Alaska, which comprises remnants of the back-arc basin that separated the Yukon-Tanana terrane from the Laurentian continental margin, confirm the presence of a late middle Permian volcanic arc component to the terrane. Our approximately 319 Ma U-Pb zircon age from the Chicken assemblage (as redefined in this study) in eastern Alaska, combined with previously reported fossil ages and a U-Pb zircon age from this assemblage, indicate that it is a Late Mississippian to Early Pennsylvanian arc assemblage. We propose several other relatively young, locally developed arc assemblages outboard of the ancient continental margin of Laurentia that may correlate with the Chicken assemblage, but we consider its origin to remain an enigma.

Alaska↗
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↗
USGS2026

Coastal and Marine Hazards and Resources Program decadal science strategy 2020–2030

Overview The decadal science strategy of the Coastal and Marine Hazards and Resources Program (CMHRP) of the U.S. Geological Survey for 2020 to 2030 describes the CMHRP's vision and mission and the strategic framework needed to support key program goals: Conduct research and develop science-based tools that lead to safer, more productive coastal communities and improved stewardship of natural resources.

Circular↗
USGS2023

An inventory of three-dimensional geologic models—U.S. Geological Survey, 2004–22

A database of spatial footprints and characteristics of three-dimensional geological models that were constructed by the U.S. Geological Survey between 2004 and 2022 was compiled as part of ongoing development of subsurface geologic information by the USGS National Cooperative Geologic Mapping Program. This initial inventory resulted in the compilation of 38 three-dimensional geological models that vary widely in their spatial extent, the type and purpose of the model, the number of subsurface units characterized by the model, and the software platforms used to create the model. This Data Report provides the scientific rationale and explanation of the contents of a companion USGS digital data release of spatial data and attributes associated with each three-dimensional model.

Data Report↗
USGS2024

On connecting hydro-social parameters to vegetation greenness differences in an evolving groundwater-dependent ecosystem

Understanding groundwater-dependent ecosystems (i.e., areas with a relatively shallow water table that plays a major role in supporting vegetation health) is key to sustaining water resources in the western United States. Groundwater-dependent ecosystems (GDEs) in Colorado have non-pristine temporal and spatial patterns, compared to agro-ecosystems, which make it difficult to quantify how these ecosystems are impacted by changes in water availability. The goal of this study is to examine how key hydrosocial parameters perturb GDE water use in time and in space. The temporal approach tests for the additive impacts of precipitation, surface water discharge, surface water mass balance as a surrogate for surface–groundwater exchange, and groundwater depth on the monthly Landsat normalized difference vegetation index (NDVI). The spatial approach tests for the additive impacts of river confluences, canal augmentation, development, perennial tributary confluences, and farmland modification on temporally integrated NDVI. Model results show a temporal trend (monthly, 1984–2019) is identifiable along segments of the Arkansas River at resolutions finer than 10 km. The temporal impacts of river discharge correlate with riparian water use sooner in time compared to precipitation, but this result is spatially variable and dependent on the covariates tested. Spatially, areal segments of the Arkansas River that have confluences with perennial streams have increased cumulative vegetation density. Quantifying temporal and spatial dependencies between the sources and effects of GDEs could aid in preventing the loss of a vulnerable ecosystem to increased water demand, changing climate, and evolving irrigation methodologies.

Colorado↗
USGS2025

Ages and trace element fertility of porphyry-related mineralization in the Philipsburg polymetallic district, Montana, with a comparison to Butte

The Philipsburg mining district is a Mo-Cu porphyry system with associated Cordilleran polymetallic veins. Geochronology was employed to date the porphyry (~66 Ma, U/Pb in zircon) and molybdenite mineralization from the veins (~76 Ma, Re-Os). Age results suggest that the two-mineralization events model proposed by Lund et al. (2018) for the Butte district can be applied to the Philipsburg deposit. Furthermore, common fertility indicators from zircon trace elements were analyzed with variable success to characterize the existing Mo-Cu mineralization of the Philipsburg porphyry. This raises concerns about the benefits of this method and/or the importance of sampling, which could impact exploration for similar porphyry deposits.

Montana↗
USGS2024

Peak streamflow trends in South Dakota and their relation to changes in climate, water years 1921–2020

Peak-flow (flood) frequency analysis is essential to water-resources management applications, including the design of critical infrastructure such as bridges and culverts, and floodplain mapping. Federal guidelines for performing peak-flow flood frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C, which documents the guidelines for determining annual peak streamflow frequency, is that, for basins without major hydrologic alterations (for example, regulation, diversion, and urbanization), statistical properties of the distribution of annual peak streamflows are stationary; that is, the mean, variance, and skew are constant through time. Nonstationarity is a statistical property of a peak-flow series such that the long-term (on the order of decades) distributional properties change one or more times either gradually or abruptly through time. Individual nonstationarities may be attributed to one source such as flow regulation, land-use change, or climate but are often the result of a combination of sources, making detection and attribution of nonstationarities challenging. In response to a growing concern regarding nonstationarity in peak streamflows in the region, the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, assessed the potential nonstationarity in peak streamflows in the north-central United States. This chapter characterizes the effects of natural hydroclimatic shifts and potential climate change on annual peak streamflows in the State of South Dakota. Annual peak and daily streamflow as well as model-simulated gridded climatic data were examined for temporal monotonic trends, change points, and other statistical properties indicative of changing climatic and environmental conditions. Changes in annual peak and daily flows were evaluated among 13, 35, and 81 qualifying U.S. Geological Survey streamgages for the 75-, 50-, and 30-year trend periods through water year 2020 (the period from October 1, 2019, to September 30, 2020) in South Dakota, respectively. No qualifying streamgages were in the 100-year trend period in the State. Statistical tests for autocorrelation (independent and identically distributed assumption), monotonic trends, and change points in the median and scale are analyzed to evaluate potential stationarity violations (nonstationarity) for performing at-site peak-flow flood-frequency analysis. The trends are reported using a likelihood approach as an alternative to simply reporting significant trends with an arbitrary p -value cutoff point. A distinct east-west spatial pattern of likely upward and downward monotonic trends and change points, respectively, was detected in 75- and 50-year trend periods, but an inconsistent spatial pattern was detected in the 30-year trend period. Additionally, change points in the median annual peak streamflows were detected in the late 1970s and early 1980s in the western part of the State, but in the east, the change point was more commonly detected in 1992–93. A similar east-west spatial pattern of likely upward and downward trends was detected in the annual peak-flow timing, the day of the year of the annal peak streamflow. In the western part of the State, the annual peak streamflows are arriving earlier, but in the east, the annual peak streamflows are arriving later. A peaks-over-threshold (POT) analysis where, on average, there are two events per year (POT2) and four events per year (POT4) was also used to evaluate changes in the frequency (count) of daily streamflows exceeding the threshold. Similar to detected changes in the annual peak streamflow, an east-west likely upward or downward change corresponding to an increase or decrease, respectively, in the frequency of daily streamflow greater than a POT2 and POT4 threshold was detected. A monthly water-balance model was used to evaluate hydroclimatic variation in annual and seasonal precipitation, snowfall, potential evapotranspiration, and soil moisture storage for all qualifying streamgages in the 75-, 50-, and 30-year trend periods. Detected trends in the annual hydroclimatic metrics for the 75- and 50-year trend periods indicate a spatially consistent statewide increase in precipitation, decrease in snowfall, increase in potential evapotranspiration, and increase in soil moisture storage. Furthermore, detected trends in seasonal precipitation in the 75- and 50-year trend periods highlight a pronounced change in precipitation in winter and later into the summer season, especially in the 50-year trend period in the eastern part of the State. Statewide increases in seasonal soil moisture storage were also detected, highlighting year-round increasing flood magnitudes, particularly in the eastern part of the State. Based on the results of these stationarity tests for the qualifying streamgages in South Dakota among the 75-, 50-, and 30-year trend periods, consistent temporal and spatial patterns of nonstationarity were detected among the 75- and 50-year trend periods. Furthermore, when nonstationarity is detected in daily streamflow, increased streamflow and volume (increasing frequency in POT), as well as potentially bridge scour, may have implications on culvert and highway design in the eastern part of South Dakota. Thus, when performing at-site peak-flow flood-frequency analyses in South Dakota, potential nonstationarities and alternative approaches are important considerations.

South Dakota↗