USGS Science⌕ Search

SEARCH · USGS Science

Results for “Ore and Energy Resource Geology”

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,664 records · Page 24Linked to original sources
USGS1995

Evaluation of ground-water flow and hydrologic budget for Lake Five-O, a seepage lake in northwestern Florida

Temporal and spatial distributions of ground-water inflow to, and leakage from Lake Five-O, a softwater, seepage lake in northwestern Florida, were evaluated using hydrologic data and simulation models of the shallow ground-water system adjacent to the lake. The simulation models indicate that ground-water inflow to the lake and leakage from the lake to the ground-water system are the dominant components in the total inflow (precipitation plus ground-water inflow) and total outflow (evaporation plus leakage) budgets of Lake Five-O. Simlulated ground-water inflow and leakage were approximately 4 and 5 times larger than precipitation inputs and evaporative losses, respectively, during calendar years 1989-90. Exchanges of water between Lake Five-O and the ground-water system were consistently larger than atmospheric-lake exchanges. A consistent pattern of shallow ground-water inflow and deep leakage was also evident throughout the study period. The mean time of travel from ground-water that discharges at Lake Five-O (time from recharge at the water table to discharge at the lake) was estimated to be within a range of 3 to 6 years. Flow-path evaluations indicated that the intermediate confining unit probably has a negligible influence on the geochemistry of ground-water inflow to Lake Five-O. The hydrologic budgets and flow-path evaluations provide critical information for developing geochemical budgets for Lake Five-O and for improving the understanding of the relative importance of various processes that regulate the acid-neutralizing capacity of softwater seepage lakes in Florida.

Florida↗
USGS1996

Pesticides in surface water from three agricultural basins in south-central Georgia, 1993-97

Twenty-two of 43 pesticides analyzed were detected in 128 water samples collected from the Tucsawhatchee Creek, the Little River, and the Withlacoochee River. These streams drain agricultural basins in south-central Georgia and were sampled from March 1993 through June 1995. Herbicides were detected more frequently than insecticides. The most frequently detected herbicides were atrazine and metolachlor and the most frequently detected insecticide was carbaryl. Pesticide concentrations in the three streams were low and did not exceed U.S. Environmental Protection Agency drinking water standards. The maximum pesticide concentration was 2.6 mg/L (micrograms per liter) for propargite, a miticide detected in only one sample. The maximum concentrations of the remaining 21 pesticides were less than 0.25 mg/L. The median concentrations were equal to the method detection limit for all pesticides except atrazine (0.008 mg/L) and metolachlor (0.012 mg/L). The ratio of herbicide detections to nondetections was largest in the planting season, smaller in the harvest season and smallest in the fallow season for the three basins. The same pattern existed for the insecticide ratios in the Little River and the Withlacoochee River. Pairwise correlations between concentrations of atrazine and metolachlor, and four parameters (discharge, and concentrations of dissolved organic carbon, suspended organic carbon, and suspended sediment) were evaluated for each stream. The strongest correlations existed between metolachlor and mean daily discharge, and metolachlor and sediment in the Withlacoochee River. The only significant correlation for the Little River was between atrazine and suspended sediment.

Georgia↗
USGS2026

A targeted approach for mapping groundwater discharge to surface water and fish thermal refuge in four Lake Ontario tributaries

The duration, magnitude, and frequency of heatwaves are predicted to increase in the coming decades, a combination that can reduce the survival of many fish species. Across the world, there is broad interest in identifying thermal refuge for heat-intolerant fish species and exploring opportunities to enhance or protect these areas. Because deeper groundwater maintains a relatively constant temperature, groundwater-influenced areas along streams can provide cool-water refuge for fish during periods of extreme heat. A targeted approach was developed for identifying existing cold-water zones and areas of substantial groundwater discharge in four high priority Lake Ontario tributaries. Our approach included: (1) predicting where groundwater discharge is most likely with a simple geospatial model and (2) using model predictions to select field sites for intensive high-resolution study, including ground-based mapping of groundwater features (springs, seeps, tributaries) as well as drone-based optical and thermal infrared surveys. Results from field sites were used to both verify model performance and map different types and aerial extents of thermal anomalies. Geospatial modelling successfully predicted regions of widespread groundwater upwelling, later verified and mapped by field and drone surveys. Comparison of model and field survey results further highlighted specific geospatial layers, such as soil/bedrock types and topographic wetness index, as being particularly useful for predicting groundwater influence on streams in the study area. In addition, a comparison of geospatial model results with a model of fish abundances along the studied streams showed significant positive correlations for many heat-intolerant fish species over a wide geographic area. The approach developed in this study can be applied to other watersheds to highlight areas of probable groundwater discharge and could be used by fishery and water resource managers to support cold-water fish habitat management decision-making and resource conservation.

New York↗
USGS2025

Toward a new framework to evaluate process-based model configurations and quantify data worth prior to calibration

Model criticism, discrimination, and selection methods often rely on calibrated model outputs. Because calibration can be computationally expensive, model criticism can first be undertaken by assessing model outputs obtained from limited prior parameter ensembles. However, such prior-based methods are often heuristic and do not formalize the notion of balancing model consistency with data and model complexity (i.e., model adequacy). We present a new framework to discriminate among candidate models prior to calibration that formalizes prior-to-calibration model adequacy into a metric to implicitly balance prior model output data coverage with model complexity represented by prior output (co)variance. The prior model adequacy metric “Mahalanobis distance deviation” quantifies the deviation of (a) the set of squared Mahalanobis distances of data from a prior model output distribution from (b) the set of squared Mahalanobis distances of data from their own distribution. A new data worth metric “discernment value” is also presented which quantifies the value of data for screening less-adequate models prior to calibration. Discernment value is calculated from the change in variance of a weighted average of prior model outputs from all candidate models due to less-adequate model outputs receiving lower weight. The framework is demonstrated using a one-dimensional groundwater flow model with eight possible configurations. A synthetic data network is used to test the framework. Results show the framework identifies the candidate models most similar to the true model used to create the synthetic data. Discernment values show variation in the value of different data types and locations for screening less-adequate models.

Water Resources Research↗
USGS1995

Methods for estimating the magnitude and frequency of peak discharges of rural, unregulated streams in Virginia

Methods are presented for estimating the peak discharges of rural, unregulated streams in Virginia. A Pearson Type III distribution is fitted to the logarithms of the unregulated annual peak-discharge records from 363 stream-gaging stations in Virginia to estimate the peak discharge at these stations for recurrence intervals of 2 to 500 years. Peak-discharge characteristics for 284 unregulated stations are divided into eight regions based on physiographic province, and regressed on basin characteristics, including drainage area, main channel length, main channel slope, mean basin elevation, percentage of forest cover, mean annual precipitation, and maximum rainfall intensity. Regression equations for each region are computed by use of the generalized least-squares method, which accounts for spatial and temporal correlation between nearby gaging stations. This regression technique weights the significance of each station to the regional equation based on the length of records collected at each cation, the correlation between annual peak discharges among the stations, and the standard deviation of the annual peak discharge for each station. Drainage area proved to be the only significant explanatory variable in four regions, while other regions have as many as three significant variables. Standard errors of the regression equations range from 30 to 80 percent. Alternate equations using drainage area only are provided for the five regions with more than one significant explanatory variable. Methods and sample computations are provided to estimate peak discharges at gaged and engaged sites in Virginia for recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years, and to adjust the regression estimates for sites on gaged streams where nearby gaging-station records are available.

Virginia↗
USGS2025

Quantifying groundwater response and uncertainty in beaver-influenced mountainous floodplains using machine learning-based model calibration

Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model-data integration workflow to address uncertainty in groundwater response to beaver-induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning-based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down-valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down-valley underflow, from 2% during dry pond periods to 20% during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down-valley underflow, minimizes the effect of beaver-induced inundation on water quality. We emphasize the need for field-scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response.

Colorado↗
USGS2024

Beyond the wedge: Impact of tidal streams on salinization of groundwater in a coastal aquifer stressed by pumping and sea-level rise

Saltwater intrusion (SWI) is a well-studied phenomenon that threatens the freshwater supplies of coastal communities around the world. The development and advancement of numerical models has led to improved assessment of the risk of salinization. However, these studies often fail to include the impact of surface waters as potential sources of aquifer salinity and how they may impact SWI. Based on field-collected data, we developed a regional, variable-density groundwater model using SEAWAT for east Dover, Delaware. In this location, major users of groundwater from the surficial aquifer are the City of Dover and irrigation for agriculture. Our model includes salinized marshland and tidal streams, along with irrigation and municipal pumping wells. Model scenarios were run for 100 years and included changes in pumping rates and sea-level rise (SLR). We examined how these drivers of SWI affect the extent and location of salinization in the surficial aquifer by evaluating differences in chloride concentration near surface waters and the subsurface freshwater-saltwater interface. We found the presence of the marsh inverts the typical freshwater-saltwater wedge interface and that the edge of the interface did not migrate farther inland. Additionally, we found that tidal streams are the dominant pathways of SWI at our site with salinization from streams being exacerbated by SLR. Our results also show that spatial distribution of pumping affects both the magnitude and extent of salinization, with an increase in concentrated pumping leading to more intensive salinization than a more widely distributed increase of the same total pumping volume.

Delaware↗
USGS1994

Hydrogeologic atlas of aquifers in Indiana

Aquifers in 12 water-management basins of Indiana are identified in a series of 104 hydrogeologic sections and 12 maps that show the thickness and configuration of aquifers. The vertical distribution of water-bearing units and a generalized potentiometric profile are shown along 3,500 miles of section lines that were constructed from drillers' logs of more than 4,200 wells. The horizontal scale of the sections is 1:125,000. Maps of aquifers showing the areal distribution of each aquifer type were drawn at a scale of 1:500,000. Unconsolidated aquifers are the most widely used aquifers in Indiana and include surficial, buried, and discontinuous layers of sand and gravel. Most of the surficial sand and gravel is in large outwash plains in northern Indiana and along the major rivers. Buried sand and gravel aquifers are interbedded with till deposits in much of the northern two-thirds of Indiana. Discontinuous sand and gravel deposits are present as isolated lenses, primarily in glaciated areas. The bedrock aquifers generally have lower yields than most of the sand and gravel aquifers; however, bedrock aquifers are areally widespread and are an important source of water. Bedrock aquifer types consist of carbonates; sandstones; complexly interbedded sandstones, siltstones, shales, limestones, and coals; and an upper weathered zone in low permeability rock. Carbonate aquifers underlie about one-half of Indiana and are the most productive of the bedrock aquifers. The other principal bedrock aquifer type, sandstone, underlies large areas in the southwestern one-fifth of Indiana. No aquifer is known to be present in the southeastern corner of Indiana.

Water-Resources Investigations Report↗
USGS2026

A roadmap for identifying and interpreting physical processes and national water model prediction bias associated with baseflow index regimes across the contiguous United States

Understanding how groundwater–surface water interactions shape streamflow variability is critical for diagnosing low flow behavior and prediction bias in continental scale hydrologic models. We present a process informed framework that links observed baseflow (BF) dynamics, watershed attributes, and National Water Model (NWM) performance across the contiguous United States. Using daily observed streamflow from 797 reference quality streamgages, we developed monthly baseflow index (BFI) signatures using a streamgage specific, calibrated digital filter. Hierarchical clustering of these signatures identified seven distinct BFI regimes capturing regional and seasonal variability. We evaluated NWM v3.0 retrospective streamflow performance within each regime using multiple hydrograph and flow duration curve-based metrics. Model skill varied systematically across regimes: mixed flow systems were simulated most accurately, while predominantly BF dominated and quickflow dominated regimes exhibited substantially poorer performance. Across nearly all regimes, the NWM underestimated observed BFI magnitude and frequently failed to reproduce seasonal BF patterns, indicating systematic biases in simulated low flow contributions. To relate these regimes to potential process controls, we trained a Random Forest classifier using static watershed attributes and applied Shapley Additive Explanations to identify features most strongly associated with each regime. Results highlight regionally varying influences, including the dominant role of snow fraction and seasonal runoff timing in snow dominated basins and the importance of evapotranspiration and aridity in quickflow dominated systems. Collectively, these findings demonstrate how hydrologic signatures combined with interpretable machine learning can diagnose regime specific model biases and generate process-based hypotheses about limitations in large scale hydrologic prediction systems.

contiguous United States↗
USGS2000

Hydrology and water and sediment quality at James Campbell National Wildlife Refuge near Kahuku, Island of Oahu, Hawaii

The James Campbell National Wildlife Refuge occupies two lowland marsh and pond complexes on the northern coastal plain of Oahu: the mostly natural ponds and wetlands of the Punamano Unit and the constructed ponds of the Kii Unit. The U.S. Fish and Wildlife Service manages the Refuge primarily to protect and enhance habitat for four endangered species of Hawaiian waterbirds. Kii Unit is fed by artesian wells and rainfall, whereas Punamano Unit is fed naturally by rainfall, runoff, and ground-water seepage. Streams drain from the uplands into lowland ditches that pass through Kii Unit on their way to the ocean. A high-capacity pump transfers water from the inner ditch terminus at Kii to the ocean outlet channel. Stormwaters also exit the inner ditch system over flood-relief swales near the outlet pump and through a culvert with a one-way valve. A hydrologic investigation was done from November 1996 through February 1998 to identify and quantify principal inflows and outflows of water to and from the Refuge, identify hydraulic factors affecting flooding, document ground-water/surface-water interactions, determine the adequacy of the current freshwater supply, and determine water and sediment quality. These goals were accomplished by installing and operating a network of stream-gaging stations, meteorology stations, and shallow ground-water piezometers, by computing water budgets for the two Refuge units, and by sampling and analyzing water and pond-bottom sediments for major ions, trace metals, and organic compounds. Streamflow during the study was dominated by winter stormflows, followed by a gradual recession of flow into summer 1997, as water that had been stored in alluvial fans drained to lowland ditches. Outflow at the ditch terminus in 1997 was 125 million gallons greater than measured inflow to the coastal plain, mainly reflecting gains from ground water along the ditches between outlying gages and the ditch terminus. Of the measured 1997 outflow, 98 percent was through the Kii outlet pump, with the outlet culvert valve only opening for brief periods during storms. Large volumes of stormflow overflowed the flood-relief swales unmeasured. The largest storm of the study, in November 1996, was estimated to have a flood frequency of about 3 to 4 years. Streamflow exceeded culvert capacity and overtopped Kamehameha Highway at Kalaeokahipa Stream and Hospital ditch. Slight overbank flooding in Kii ditch resulted strictly from high discharge. Minor overbank flooding farther out on the coastal plain probably was caused mainly by the small hydraulic gradients available to convey stormflows along the lowland ditches. Stormwaters flooded Kii ponds and flowed back upstream along Punamano ditch into Punamano marsh, introducing suspended sediment and possibly other contaminants to the Refuge. Two smaller storms in January 1997 resulted in smaller flows and no overbank flooding. The Kii outlet pump ran continuously for 7 days during the November 1996 storm and for 1 to 2 days during the January 1997 storms. During all three storms, the outlet culvert valve opened and the inner ditches overtopped the flood-relief swales, allowing free outflow of water from the inner ditch. Backwater effects hindered drainage during the January 1997 storms at Hospital ditch at Kamehameha Highway, and at Punamano ditch at Nudist Camp Road (where the backflow into Punamano marsh in November 1996 constituted an extreme backwater effect). A probable marine backwater effect was imposed at the ocean outlet ditch during the November 1996 storm through a combination of high spring tides and wave setup from large surf. Whether this backwater effect propagated upstream in the ditches to affect inland sites could not be determined conclusively. A sand plug may have built up in the ocean outlet channel before the November 1996 storm, but if so, it probably washed out prior to, or early in the storm, and was not present at the time of peak stage at inland sites. A season-long buildup of the sand plug in late 1997 was inferred from rising water levels in the outlet ditch. Seawater flows up the outlet channel or over the sand berm and into the outer ditch system on most high tides, and particularly during spring high tides. Ponds and ditches of the Refuge and surrounding lowlands have mud- and clay-lined bottoms that form an effective confining unit and inhibit interaction with an underlying shallow limestone aquifer. At Kii Unit, pond levels are higher than adjacent ditch levels and underlying ground-water levels, establishing lateral and downward head gradients that could foster seepage losses from the ponds. Regional ground-water discharge from the Koolau aquifer to the coastal-plain sediments is mostly diffuse, but is concentrated where ridges of Koolau Basalt plunge beneath coastal-plain sediments near Punamano Unit and at the head of Hospital ditch. Kii ditch gains brackish ground water downstream of Kamehameha Highway. Wastewater disposal from the sewage treatment plant adjoining Kii Unit poses little or no threat to Refuge habitat. Disposal is at six injection wells located 0.45 mi away at Kahuku, and the wells inject into confined limestone aquifers that do not extend to Kii Unit. The natural freshwater supply to Punamano Unit is adequate for maintaining the wildlife habitat, judging from stable pond levels and low salinities there. A monthly water budget for Punamano showed an apparent annual deficit in measurable flows in 1997, requiring unmeasured ground-water gains equalling 51 inches of water. The freshwater supply to Kii Unit is inadequate according to Refuge managers, because there is not enough water to manipulate levels adequately in the ponds during most of the year, and particularly during the driest months. This is confirmed by monthly deficits in the water budget for the Kii ponds during summer months. However, the Kii budget showed an annual surplus in measurable flows for 1997 equalling 24 inches of water. Unmeasured losses are required to explain the apparent annual surplus, such as discharge to the ditches through pond water-control structures and downward and/or lateral ground-water seepage. The apparent surplus at Kii is strictly hydrologic and is not a surplus in a management sense; it cannot be stored or used to supply the Refuge, but instead reflects losses from the system that render this amount of water unavailable for use. The budget year, 1997, was drier than normal (24 percent below long-term mean rainfall) and so the measured potential evaporation for 1997 was probably higher than the long-term mean. Few metals or organic compounds of potential concern were detected in pond and ditch waters and in pond-bottom sediments. Detected pesticides were at trace levels or just above minimum reporting limits. Exceptions that exceeded quality guidelines for freshwater sediment were copper and zinc in sediment from Kii ponds C and D, and copper in sediment from Punamano north pond. Therefore, urban and agricultural runoff probably have contributed little in the way of harmful metals or organic compounds to the Refuge, although the potential for such contribution remains from periodic flooding of the ponds by ditch stormflows. Salinity was low throughout most Refuge waters, qualifying as fresh to slightly brackish and suitable for the environmental needs of Refuge fauna. Higher salinities have been observed in ditches during past periods of sugarcane cultivation and saltwater aquaculture, however. Resumption of saltwater aquaculture could raise ditch salinities if saltwater effluents are disposed directly into the ditches, as they were in the past.

Hawaii↗
USGS2025

Flood-inundation maps for 14.8 miles of Little and Big Papillion Creeks in Omaha, Nebraska, 2023

Digital flood-inundation map libraries for two reaches that constitute 14.8 miles of Little and Big Papillion Creeks in Omaha, Nebraska, were created by the U.S. Geological Survey (USGS) in cooperation with the Papio-Missouri River Natural Resource District. The flood-inundation maps, which can be accessed through the USGS Flood Inundation Mapping Program website at https://www.usgs.gov/mission-areas/water-resources/science/flood-inundation-mapping-fim-program , depict estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) at Little Papillion Creek at Irvington, Nebr. (USGS station 06610750), Little Papillion Creek at Ak-Sar-Ben at Omaha, Nebr. (USGS station 06610765), and Big Papillion Creek at Q Street at Omaha, Nebr. (USGS station 06610770) streamgages. Near-real-time stages at these streamgages may be obtained from the USGS National Water Information System database at https://doi.org/10.5066/F7P55KJN or from the National Weather Service Advanced Hydrologic Prediction Service at https://water.weather.gov/ahps/ . Flood profiles were computed for two different reaches that constitute 14.8 miles of stream length in the study area by using hydraulic models. The models were calibrated by adjusting roughness coefficients to best represent the current (2022) stage-streamflow relation at the streamgages within the study reach. The hydraulic models were then used to compute water-surface profiles at 1-foot stage intervals for selected stage ranges to represent various flooding scenarios at the streamgages in each reach. The simulated water-surface profiles then were combined with a digital elevation model using a geographic information system, which had a 10-foot grid spacing to delineate the flooding extents and water depths for each stage. The availability of these flood-inundation maps, along with information regarding current stage from the USGS streamgages, can provide emergency management personnel and residents with information that is critical for flood response activities and post flood recovery efforts.

Nebraska↗
USGS2026

Top-down targeted network analysis of critical mineral commodities applied to international geochemistry database

The global demand for critical mineral commodities is rapidly increasing, making domestic production an important factor in supporting the economy and national security. Large scale, publicly available geochemical databases allow for the application of data informatics methods to interrogate critical mineral commodities data for correlations in deposit formation and distribution, particularly for identifying enrichment of multiple critical mineral commodities at the same deposit. In this study, we applied network analysis to the Critical Minerals Mapping Initiative (CMMI) ore geochemistry (Critical Minerals in Ores, CMiO) database to identify the high concentration (defined as 10× bulk crustal abundance) co-occurrence of different critical mineral commodities across a mineral system hierarchy from deposit environments to individual deposits. Identifying patterns or unique outliers in enrichment in network communities will allow for the location of secondary critical mineral commodity resources from under-utilized deposits. We find trends in the enrichment of critical mineral commodities in network-communities between the elements praseodymium (Pr), neodymium (Nd), terbium (Tb), and dysprosium (Dy) across multiple CMiO database deposit environments and groups down to specific deposit types and sites. A separate trend in network community deposition is observed as well between iridium (Ir) and platinum (Pt) in deposit environments, groups, types, and sites. Network analysis focused on critical minerals in magmatic-hydrothermal deposits identified multiple deposit sites from different deposit types within the CMiO database with concentrations of Dy, Nd, Tb, Pr, Ir, and Pt that are at least ten times greater than the crustal average. This approach can be applied to any target element(s) or deposit(s) of interest, allowing broad investigation of co-enriched critical mineral commodities.

Journal of Geochemical Exploration↗
USGS2026

Spectral characteristics and mapping of lithium-rich playas in the Basin and Range Province, western United States

This report presents a comprehensive set of prospectivity maps for lithium in brines located within the playas of the Basin and Range Province of western United States. The prospectivity assessment is based on criteria established in the descriptive deposit model outlined in a previous U.S. Geological Survey report, with a focus on regional remote sensing methodologies. The playas within the study area have been systematically ranked according to several key factors, including size, the presence and abundance of source rocks, vegetation as an indicator of water availability, documented prospects, and data products derived from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) sensor. Notably, the deposits in the playa of Clayton Valley make it stand out as the sole lithium brine-producing playa in the United States. The ASTER sensor provides optimal regional coverage, characterized by a 60-kilometer swath width and unique multispectral capabilities that are not available with other remote sensing instruments. This report and the accompanying maps are intended to serve as resources for further multidisciplinary assessments of lithium or other elements or commodities where prospectivity can be inferred from the presence of proxies such as alteration minerals.

Arizona, Baja California, California, Chihuahua, I↗
USGS2026

Asynchronous landslide seasonality across the United States

Mid-range landslide outlooks can facilitate weather-related landslide preparedness and disaster response planning, but seasonal landslide activity remains poorly quantified at continental scales. Leveraging >55,000 reported landslides from across the United States (U.S.), we used circular statistics to quantify landslide seasonality in 67 National Weather Service County Warning Areas (CWAs). We found regional differences in landslide season timing and duration, with transitions between domains variably corresponding to climate class or river basin. We assessed differences in seasonality by movement type for slides, flows, and falls, detecting apparent, but uncertain, differences between slide and fall seasonalities in 27 of 35 (77%) of CWAs with both types reported. In the Pacific Northwest, where long records exist, we found a credible shift toward a later mean landslide season in western Washington from 1990 to 2020, but no trend in western Oregon. Our results can provide emergency planners a resource to assess seasonal landslide probability nationwide.

Geophysical Research Letters↗
USGS2024

Managing basin-scale carbon sequestration: A tragedy of the commons approach

The Tragedy of the Commons is a well studied problem in the literature of ecology, economics, and environmental policy which illustrates the deleterious consequences of managing common pool resources when individual and social incentives are misaligned. In this work, we apply a simple model of carbon sequestration in a deep saline aquifer by two neighboring geologic CO 2 storage (GCS) operators to begin investigating if a Tragedy of the Commons framework applies to GCS. Specifically, we consider the pressure space as a “commons” because the injection by each firm at its own well increases the downhole injection pressure at both wells. We assume that a firm will decrease its injection rate if the downhole pressure at its well exceeds a predefined maximum (i.e., exceeds the “pressure limit”). With this assumption in place, we find that the same injection flowrates are optimal for both wells, regardless of whether they are owned by the same firm or competing firms. This suggests that GCS may not be best represented by a pure Tragedy of the Commons framework under our initial assumptions. However, there could be economic incentives or contractual obligations that may result in either or both GCS operators being unwilling to reduce their injection rates. Thus, we conclude the conference paper with a discussion of future extensions of our approach that may demonstrate closer alignment with the Tragedy of the Commons, including explicit definitions of pore-space rights, firm uncertainty regarding the parameters of the Theis equation, and the potential role of unitization.

Conference Paper↗
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↗
USGS2026

Where will the cat cross the road? Comparing camera and GPS-based models for identifying wildlife corridors

Designing effective wildlife corridors is a critical conservation challenge in fragmented landscapes. GPS-based step selection functions strongly predict dispersal corridors and connectivity, but GPS collaring can be expensive and invasive. Camera-based occupancy models are widely used for connectivity analyses but may involve trade-offs in data resolution. Despite widespread use of both approaches, few studies have directly compared them using concurrent datasets. We developed a stacked single-species, single-season occupancy model and a Circuitscape connectivity surface for mountain lions (Puma concolor) on Washington’s Olympic Peninsula, USA, and compared them with a connectivity surface from an existing integrated step selection function. Both models predicted mountain lion GPS locations well, with binned Spearman rank correlations of 1 for Circuitscape and 0.96 for the step selection function, though step selection better identified habitat use by dispersers. Connectivity predictions were moderately correlated across the landscape ( r = 0.26), but agreement was strongest in human-dominated areas most critical for corridor planning. We conclude that GPS-based approaches are advantageous when data collection is feasible and the focus is on dispersal or fine-scale movement. However, camera-based approaches may be preferable for multi-species monitoring, large spatial and temporal scales, noninvasive sampling, when resources are limited, or when fine-scale or dispersal-specific inference is not required.

Washington↗
USGS1997

An optimized network for phosphorus load monitoring for Lake Okeechobee, Florida

Phosphorus load data were evaluated for Lake Okeechobee, Florida, for water years 1982 through 1991. Standard errors for load estimates were computed from available phosphorus concentration and daily discharge data. Components of error were associated with uncertainty in concentration and discharge data and were calculated for existing conditions and for 6 alternative load-monitoring scenarios for each of 48 distinct inflows. Benefit-cost ratios were computed for each alternative monitoring scenario at each site by dividing estimated reductions in load uncertainty by the 5-year average costs of each scenario in 1992 dollars. Absolute and marginal benefit-cost ratios were compared in an iterative optimization scheme to determine the most cost-effective combination of discharge and concentration monitoring scenarios for the lake. If the current (1992) discharge-monitoring network around the lake is maintained, the water-quality sampling at each inflow site twice each year is continued, and the nature of loading remains the same, the standard error of computed mean-annual load is estimated at about 98 metric tons per year compared to an absolute loading rate (inflows and outflows) of 530 metric tons per year. This produces a relative uncertainty of nearly 20 percent. The standard error in load can be reduced to about 20 metric tons per year (4 percent) by adopting an optimized set of monitoring alternatives at a cost of an additional $200,000 per year. The final optimized network prescribes changes to improve both concentration and discharge monitoring. These changes include the addition of intensive sampling with automatic samplers at 11 sites, the initiation of event-based sampling by observers at another 5 sites, the continuation of periodic sampling 12 times per year at 1 site, the installation of acoustic velocity meters to improve discharge gaging at 9 sites, and the improvement of a discharge rating at 1 site.

Florida↗