SEARCH · USGS Science

Results for “Canadian Journal of Fisheries and Aquatic Sciences”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,688 records · Page 62Linked to original sources

USGS2025

Water supply in the conterminous United States, Alaska, Hawaii, and Puerto Rico, water years 2010–20

We present an assessment of water supply across the conterminous United States (CONUS), Alaska, Hawaii, and Puerto Rico covering water years 2010–20. Our analysis drew on two national hydrologic models, the National Hydrologic Model Precipitation-Runoff Modeling System and the Weather Research and Forecasting model hydrologic modeling system. Both models produced estimates of streamflow, evapotranspiration, soil moisture, snow water equivalent, and other hydrologic states and fluxes. The models were driven by the bias-adjusted 4-kilometer-resolution, long-term regional hydroclimate simulation over the conterminous United States dataset (CONUS404). We assessed spatial and temporal error distributions by comparing monthly simulations at the 12-digit hydrologic unit code and regional scale from both models against external benchmarking datasets. Results showed that average annual rainfall across the CONUS was 857 millimeters per year for the period of analysis, with water year 2012 the driest year (729 millimeters) and water year 2019 the wettest year (995 millimeters). Key interannual variability results included the following: (1) the California–Nevada hydrologic region had the highest variability in precipitation and snow accumulation, and (2) the Texas hydrologic region was among hydrologic regions with the highest variability in precipitation. We related interannual variability in precipitation to storage volumes in soil moisture, snow water equivalent, and lakes and reservoirs to highlight areas with little storage and large year-to-year variability in precipitation. These areas included the Southern High Plains, Central High Plains, Texas, Souris–Red–Rainy, Mississippi Embayment, and Midwest regions. Our analysis of groundwater-level data showed that several of these areas overlap aquifers where groundwater levels were considerably lower than historical averages, including the Colorado Plateaus aquifers, the Rio Grande aquifer system, and the Central and Southern regions of the High Plains aquifer. Many of these lowered groundwater levels are continuations of decades-long declines from overpumping that started well before the assessment period. The resulting water budgets and their analyses provide a high-resolution foundational assessment of the mean state and variability of the terrestrial hydrologic cycle across the CONUS and Alaska, Hawaii, and Puerto Rico to support a wide range of water resource management applications.

Professional Paper↗

USGS2024

Insights on gas hydrate formation and growth within an interbedded sand reservoir from well logging at the Qiongdongnan Basin, South China Sea

Although variable well log resolution and its control on saturation estimation has been studied, it has not been directly applied to a specific location to explore the nature of gas hydrate within a sand reservoir. We applied in-situ measurements of resistivities, neutron porosity, and gamma ray at two sites in the Qiongdongnan Basin, South China Sea (QDN-W05–2021 and QDN-W08–2021) to investigate the reservoir parameters of a hydrate-bearing sand reservoir. Our results show that gas hydrate is distributed in 5 zones with a total thickness of 10.7 m and an average saturation of 69% at the QDN-W05–2021 site, while they are distributed in 2 zones with a total thickness of 4.3 m and an average saturation of 49% at the QDN-W08–2021 site. We found that variances in saturations estimated from lateral-extra deep button (RX), phase shift (P40H-P40L), and attenuation (A40H-A40L) resistivities within the laterally mapped continuous sand body were affected by the nature of gas hydrate occurrences. Results indicate gas hydrate forms and accumulates at the center of the sand layer and tends to be less or not present toward the top and base. Integrated with seismic data, the in-situ measurements provide insights in the evolution of a mushroom-shaped, hydrate-gas reservoir system. In the system, free gas is likely horizontally transported from the top-center of the gas chimney to the surrounding areas in the early stage dominated by a warm-gas environment, whereas hydrate forms in the opposite pathway starting from the surrounding areas in the following stage with temperature reducing. Our study suggests that high-resolution in-situ measurements not only are a tool to identify the physical properties, but also can be used to help explain the physical process of hydrate growth and accumulation.

Marine Geology↗

USGS1976

Effect of deicing chemicals on the hydrologic environment in Massachusetts; evaluation of surface-water data collection

The objective of this study is to develop predictive relationships that can be used to describe the impact of highway deicing salts on the hydrologic environment. A method, presented in an earlier report, for estimating yearly mean chloride concentrations from estimated or actual runoff and salt-application data was applied to the period 1971 through1974 and the results compared with yearly mean chloride concentrations computed from records of specific conductance. Estimates were found to be from minus 63 percent to plus 56 percent in error when the error was expressed as the difference between estimated and computed values as a percentage of the computed values. Chloride concentrations used in this study are computed from specific conductance/chloride concentration relationships and records of specific conductance. Preliminary results from graphic analyses and least-squares regression analysis show that relationships between measured values of specific conductance and chloride concentration have correlation coefficients ranging from 0.30 to 0.97. Analysis of streamflow for all major dissolved constituents is recommended with the purpose of attempting to describe the variations in the specific conductance/chloride concentration relationships.

Massachusetts↗

USGS1983

Technique for estimating depth of floods in Tennessee

Estimates of flood depths are needed for design of roadways across flood plains and for other types of construction along streams. Equations for estimating flood depths in Tennessee were derived using data for 150 gaging stations. The equations are based on drainage basin size and can be used to estimate depths of the 10-year and 100-year floods for four hydrologic areas. A method also was developed for estimating depth of floods having recurrence intervals between 10 and 100 years. Standard errors range from 22 to 30 percent for the 10-year depth equations and from 23 to 30 percent for the 100-year depth equations.

Tennessee↗

USGS2026

Effects of human activity on Khumbu Glacier: Towards a sustainable Everest Base Camp

The Everest Base Camp (EBC), a critical staging area for mountaineers attempting to summit Mount Everest, has seen a significant increase in human activity over the past few decades. The increasing number of climbers, trekkers, and support staff, particularly at the campsite area, has intensified environmental impacts on the Khumbu Glacier. This study investigates the effects of human activity at EBC, focusing on three contributors to glacier melt: global climate change, local fossil fuel consumption, and human urinary discharge. Through field data collected during the spring 2023 climbing season, this study examines the energy released through the use of liquefied petroleum gas (LPG), kerosene, and petrol for cooking, heating, and electricity, as well as the heat generated by human urinary discharge. The total heat energy released from these activities at EBC in the spring of 2023 was 849,174 ± 179,774 MJ, which would be sufficient to melt 2492 ± 528 tons of glacier ice and snow. Trends in land surface temperature at EBC and adjacent environments during the 1991–2023 period were calculated from Landsat 5, 7, and 9 satellite data. The 32-year Landsat record reveals that EBC surface temperature increased by 0.28 °C year −1 , which was roughly twice the rate of warming at the surface of Khumbu Glacier adjacent to EBC and about 15% higher than the temperature increase of the debris-covered terrain immediately to the north of EBC. The findings suggest that addressing these anthropogenic influences could help to preserve the Khumbu Glacier and support the sustainability of mountaineering in the region. This paper also takes a transformative approach and explores stable and safer locations in case there is a necessity to relocate the current EBC. Two sites have been identified southwest of the existing EBC, which, unlike the current supraglacial site, are situated on stable ground.

Everest Base Camp, Khumbu Glacier↗

USGS2025

Earthquake magnitude and source parameter estimation with a distributed acoustic sensing dataset in the Gorda subduction zone

Distributed acoustic sensing (DAS) systems offer a cost‐effective way to create large‐scale strainmeter arrays for seismological applications using fiber‐optic cables. DAS‐based strain measurements are known to be influenced by various factors, bringing into question their general reliability for accurate earthquake characterization. A 15‐km‐long DAS deployment in northern California was operational within 3 days of the 2022 M w 6.4 Ferndale earthquake and ran continuously throughout the aftershock sequence. We utilize these aftershock data to validate DAS‐based strain measurements in two ways. We first test the accuracy of DAS‐based magnitude estimates from peak dynamic strains by comparing them with magnitude and attenuation scaling relations derived independently from traditional borehole strainmeter (BSM) data. We demonstrate that DAS‐based magnitudes are comparable to BSM‐based magnitudes when corrections for variations in site response along the fiber‐optic cable are properly made. Magnitude errors are spatially correlated, potentially because of factors such as finite‐fault effects (e.g., stress drop) or more complex, unmodeled path attenuation or because of wave propagation effects in heterogeneous media. We then apply more advanced source characterization methodology to the DAS data using a time‐domain empirical Green’s function (EGF) deconvolution approach to measure details of the moment rate history. The EGF approach using DAS data depends on careful treatment of distorting factors such as anthropogenic sources of noise and optical phase wrapping but successfully isolates source spectra for moderate‐magnitude earthquakes: source spectral ratios obtained from DAS data, broadband seismometer data, and BSM data in the same region show consistent results, revealing differences in directivity and spectral shape among earthquakes. Although further research is needed to refine source‐time‐function estimation techniques for DAS data, particularly for larger magnitude events, these case studies demonstrate the clear potential of DAS for earthquake source characterization.

California↗

USGS2000

Regional evaluation of evapotranspiration in the Everglades

Nine sites in the Florida Everglades were selected and instrumented for collection of data necessary for evapotranspiration-determination using the Bowen-ratio energy-budget method. The sites were selected to represent the sawgrass or cattail marshes, wet prairie, and open-water areas that constitute most of the natural Everglades system. At each site, measurements necessary for evapotranspiration (ET) calculation and modeling were automatically made and stored on-site at 15- or 30-minute intervals. Data collected included air temperature and humidity at two heights, wind speed and direction, incoming solar radiation, net solar radiation, water level and temperature, soil moisture content, soil temperature, soil heat flux, and rainfall. Data summarized in this report were collected from January 1996 through December 1997, and the development of site-specific and regional models of ET for this period is described. Latent heat flux is the energy flux density equivalent of the ET rate. Modified Priestley-Taylor models of latent heat flux as a function of selected independent variables were developed at each site. These models were used to fill in periods of missing latent heat flux measurement, and to develop regional models of the entire Everglades region. The regional models may be used to estimate ET in wet prairie, sawgrass or cattail marsh, and open-water portions of the natural Everglades system. The models are not applicable to forested areas or to the brackish areas adjacent to Florida Bay. Two types of regional models were developed. One type of model uses measurements of available energy at a site, together with incoming solar energy and water depth, to estimate hourly ET. This available-energy model requires site data for net radiation, water heat storage, and soil heat flux, as well as data for incoming solar radiation and water depth. The other type of model requires only incoming solar energy, air temperature, and water depth data to provide estimates of hourly ET. The second model thus uses data that are more readily available than the data required for the available-energy model. Computed ET mean annual totals for all nine sites for the 1996-97 period ranged from 42.4 inches per year at a site where the water level is below land surface for several months each year to 57.4 inches per year at an open-water site with no emergent vegetation. Although the density of photosynthetically-active plant leaves has been shown to relate directly to ET in some studies, it does not appear to relate directly to ET in the Everglades, based on comparison of annual ET data with leaf-area index, defined as the Normalized Difference Vegetation Index (NDVI), data from satellite imagery. NDVI and ET appear to be inversely related in the Everglades. The greatest ET rates occurred at open-water sites where the NDVI data indicated the lowest leaf-area index. Among the remaining vegetated sites, there is no clear relation between ET and NDVI, though the highest ET rate corresponded to the lowest NDVI and one of the lowest ET rates corresponded to the highest NDVI value. The variation in ET follows a seasonal pattern, with lowest monthly ET totals occurring in December through February, and highest ET occurring in May through August. The greatest range in monthly ET among all nine sites for the 2-year period occurred at site 3: from 1.81 inches in December 1997 to 6.84 inches in July 1996. A study to compare the Bowen-ratio/energy balance method of ET measurement with the eddy-correlation method was done at one site from June 22, 1998, through September 28, 1998. This comparison indicated that both methods gave comparable values of the Bowen ratio, but there was a considerable difference in available energy measured by the two methods. The mean of all 30-minute measured turbulent heat fluxes from the eddy-correlation apparatus for June 22 through September 29, 1998, was 137.4 watts per square meter, and the mean of the corresponding measured energy was 163.6 watts per square meter, or about 20 percent greater. The disagreement in mean energy fluxes measured by the two methods is problematical and is not fully understood. Although the difference seems to be related to friction velocity, and is practically non-existent at values of friction velocity greater than 0.3 meter per second, the "correctness" of either method cannot be determined with the data available.

Florida↗

USGS2024

Geophysical constraints on continental rejuvenation in central China: Implications for outward growth of the Tibetan Plateau

Continental rejuvenation results from the tectonic reactivation of crustal structures and lithospheric reworking by mantle flow. Geochemical observations and field mapping have traditionally provided the primary evidence for the secular evolution of crustal composition and tectonic processes during continental rejuvenation. Nonetheless, the impact of continental rejuvenation on the observed present-day strain rate and orogenic-scale lithospheric structure has not been well constrained. The pre-existing E-W–trending Central China Orogenic Belt has been overprinted by the N-S–trending Central Longitudinal Seismic Belt and constitutes the intracontinental West Qinling Syntaxis in central China, where the tectonic setting changes eastward from contraction to extension. Combining updated global positioning system data and high-resolution crustal seismic tomography, we reveal a modern continental rejuvenation process within the West Qinling Syntaxis in central China. The northward extrusion of the Tibetan Plateau's weak lithospheric layer (middle-lower crust and lithospheric mantle) of southwestern China relative to the rigid Sichuan Basin/Ordos Block of the eastern West Qinling Syntaxis results in regional dextral shearing that shapes the Central Longitudinal Seismic Belt and defines the eastern Tibetan Plateau margin. The pre-existing E-W–trending Central China Orogenic Belt has been preserved above the brittle-ductile transition zone, and the northward movement of the deep lithospheric layer drives the deformation of the upper crust in the West Qinling Syntaxis. Our results, along with previous studies, suggest the presence of an intracontinental lithospheric interchange structure in central China. The continental rejuvenation of the West Qinling Syntaxis results from a combination of fault reactivation in the upper crust (Stage I, Eocene–Oligocene) and reworking of the deep lithosphere (Stage II, middle–late Miocene) related to the plateau-wide shift in stress accommodation ultimately driven by the redistribution of mass outward from the central Tibetan Plateau. At present, the transition zone between the high- and low-velocity anomalies along the Central Longitudinal Seismic Belt not only shapes the landscape boundary but controls the size and recurrence interval of earthquakes within the West Qinling Syntaxis in central China.

GSA Bulletin↗

USGS2026

Modeled groundwater and surface-water interactions surrounding Mobile Bay, Alabama, 2008–15

The U.S. Geological Survey, in cooperation with the Gulf Coast Ecosystem Restoration Council, has used MODFLOW 6 to develop a groundwater-flow model to simulate groundwater and surface-water interactions in the Mobile Bay, Alabama, area. The model results indicated that, on average, groundwater discharge near the coastline is equal to 2.5 percent of the surface water that flows into the bay. The model was also used to determine how changes in recharge, sea level, and groundwater pumping affect groundwater levels and discharge rates. The results indicate that more groundwater discharge occurred in the winter and spring when recharge was higher, sea level was lower, and groundwater pumping was lower than during the summer and fall. Additionally, the amount of emergent groundwater was closely related to sea level; when sea level was higher, there was more area with emergent groundwater. Furthermore, the depth of nonemergent groundwater was related to trends in recharge and pumping. During periods of increased recharge and reduced pumping, a significant portion of the model area exhibited a depth to the water table of less than 1 meter below the land surface. Conversely, during periods of decreased recharge and heightened pumping, much of the model area showed a depth to the water table ranging from 1 to 5 meters.

Alabama, Mississippi↗

USGS2025

Statistical approaches for modeling correlated grade and tonnage distributions and applications for mineral resource assessments

Correlations between grade and tonnage exist in mineral resource data compiled from published reports, but they are not always addressed during quantitative assessment of undiscovered mineral resources. Failure to account for correlated grade and tonnage distributions can result in geologically unrealistic assessment results. Current software tools simulate univariate ore tonnage and multivariate resource grades of undiscovered deposits independently. As a result, analysts are forced to rely on ad-hoc solutions to minimize the correlation issues by: 1) creating subsets of data with restricted criteria; 2) truncating grade and tonnage distributions; and 3) testing model robustness using exploratory data analysis. While these methods represent pragmatic solutions, the statistical solutions presented here provide additional options to address real correlations in grade and tonnage data used for mineral resource assessments. We present a modified version of the MapMark4 package in R that introduces two alternatives for modeling grade and tonnage distributions, consisting of a multivariate solution that accounts for correlations between ore tonnage and metal grades and an empirical solution that utilizes simple random sampling with replacement to reproduce coupled grades and tonnages from the input data. We present simulations for contained ore and metal for three case studies representing tungsten skarn, komatiite-hosted nickel, and sediment-hosted carbonate amagmatic zinc-lead (Mississippi Valley-type) deposits. Employing the methods presented here yields quantitative mineral resource assessment results that more closely reflect the empirical distributions of grades and tonnages observed in nature and expands the applicability of these tools for ongoing critical mineral resource assessments.

Applied Computing and Geosciences↗

USGS2024

Groundwater and surface-water interactions in the Lower Duwamish Waterway, Seattle, Washington

The U.S. Geological Survey (USGS), in cooperation with the Washington State Department of Ecology (Ecology), conducted a study to describe the current understanding of the regional groundwater system of the lower Duwamish River valley and groundwater and surface-water interactions in the lower Duwamish Waterway. The lower Duwamish Waterway is the final 5-mile (mi) reach of the Duwamish River before it empties into Elliott Bay in Puget Sound near Seattle, Washington. A nearshore site (hereinafter referred to as “Nearshore Site” to distinguish the particular site from general discussions of nearshore areas) along the western shoreline of the Duwamish River, about 1.5 mi upstream from the river mouth, was selected for focused groundwater data collection by USGS. Data loggers were deployed in seven groundwater wells and one stilling well in the Duwamish River to measure specific conductance, temperature, and depth at 15-minute intervals for a period of about 2 years. At the Nearshore Site during 2020–22, water levels in the shallow wells were 3–8 feet (ft) higher than water levels in the deep wells, providing evidence for a low-permeability layer between the shallow and deep aquifers in this area. The shallow wells had a pronounced seasonal variability, with high water levels in winter and low water levels in summer. Data from the deep wells showed far less seasonal variability, with slight increases in winter and a near-constant water level from spring to autumn. The deep wells had a strong hydraulic connection to the Duwamish River, as evidenced by the synchronous water-level variability during the tidal cycle, whereas the shallow wells had minimal to no tidal response. The potentiometric maps developed for the Nearshore Site and surrounding areas indicate large differences in groundwater-flow directions for the shallow and deep aquifers at low and high tides. For the shallow aquifer, flow is toward the lower Duwamish Waterway near the Nearshore Site, regardless of the tidal condition. For the deep aquifer, a potentiometric trough forms parallel to the shoreline during high tide, indicating that groundwater flow converges from the uplands to the west and the Duwamish River to the east. The geometry of the potentiometric surfaces between the nearshore-most well and the shoreline is complex and is further confounded by intermittent shoreline armoring and other buried infrastructure, which could serve as either a barrier or a conduit to flow. Groundwater and surface-water interactions in the lower Duwamish Waterway are inherently complex as a result of three overarching factors. First, water levels in the lower reaches of the Duwamish River vary daily by 11–16 ft because of tides from Puget Sound, which create large swings in the hydraulic gradient in the nearshore groundwater system. Second, the density and chemical composition of water in the Duwamish River change daily with the tides and seasonally, which constrains how river water entering the nearshore sediments interacts with discharging groundwater. Third, the nearshore subsurface and shoreline conditions are heterogenous because of extensive shoreline armoring over the past century, which governs the flow of groundwater and infiltrating river water. These unique features of groundwater and surface-water interactions in the lower Duwamish Waterway thus govern the transport of terrestrial contaminants to the lower Duwamish Waterway. Furthermore, the heterogenous aquifer properties in the lower Duwamish Waterway contribute to spatially and temporally dynamic contaminant-transport processes.

Washington↗

USGS2024

Upper crustal seismic velocity structure of the Hayward fault zone, San Francisco Bay, California, USA: Results from the 2016 East Bay Seismic Experiment (EBSI-16)

We developed Vp, Vs, Vp / Vs ratio, and Poisson’s ratio models of the uppermost crust (<4 km depth) from the eastern San Francisco (SF) Bay (California, USA) to near the Calaveras fault along a 15-km-long, linear profile. Upper crustal velocities are highly variable beneath, west, and well east of the Hayward fault. We observe eight notable features, from west to east: (1) Near San Francisco Bay, there is an ~2-km-wide structure with high Vp / Vs ratios (up to 5) and Poisson’s ratios (up to 0.48) extending from the surface to the base of our model, which we suggest the structure is a near-vertical fault that lies along a straight-line projection between the Silver Creek fault to the south and the Point Richmond fault to the north. The structure may be part of an ~90-km-long fault along the eastern SF Bay. (2) The western East Bay Plain, the lower lying area between the bay and the hills, includes up to 800 m of low-velocity sediments ( Vp ~1600–3000 m/s, Vs ~500 m/s to ~1000 m/s), underlain by higher velocity basement rocks ( Vp ~3000–5800 m/s; Vs ~1000–1500 m/s). (3) Between ~1 km and 3 km east of the Bay shoreline, sediments thin in a series of steps (likely faults) toward the Hayward fault. (4) Between ~3 km west and ~1 km east of the Hayward fault (at the East Chabot fault) at depths greater than 1 km, basement Vp (up to 6000 m/s) and Vs (up to 2800 m/s) are high, and Vp / Vs ratios (<2) and Poisson’s ratios (<0.3) are low, suggesting crystalline rocks. Furthermore, a near-vertical zone of low Vp / Vs ratios and Poisson’s ratios is between near-surface traces of the Hayward and East Chabot faults, likely corresponding to the San Leandro Gabbro of Ponce et al. (2003) . (5) Eastward of the East Chabot fault in the upper 1.5 km, basement Vp (~3000 m/s to ~4200 m/s) and Vs (~1200–2000 m/s) are lower than those west of the fault. (6) In the eastern Hayward/Oakland Hills, there are zones of laterally varying, high- and low-velocity ( Vp ~2500–3000 m/s) Jurassic–Cretaceous and Tertiary sediments in the shallow subsurface that likely extend much deeper than imaged. (7) Seismic energy that propagates westward from sources east of the Hayward fault (HF) appear weaker than energy that propagates eastward from sources west of the HF, suggesting that the HF acts as a partial barrier to shallow seismic energy propagation into the more populated eastern SF Bay area. (8) Unlike many fault zones, it appears that the active trace of the Hayward fault (in our study area) is not cored by a prominent, low-velocity zone relative to rocks to the east and west of the active trace. However, the active trace does mark a prominent change from relatively higher velocities to the west and lower velocities to the east.

California↗

USGS2024

Iodine-129 in the eastern Snake River Plain aquifer at and near the Idaho National Laboratory, Idaho, 2021–22

Between the 1950s and 1980s, wastewater generated at the Idaho National Laboratory contained Iodine-129 ( 129 I); this wastewater was discharged directly into the eastern Snake River Plain (ESRP) aquifer through a deep disposal well, unlined infiltration ponds, or leaked from distribution systems below industrial facilities. During 2021–22, the U.S. Geological Survey, in cooperation with the U.S. Department of Energy and the Idaho Department of Environmental Quality Idaho National Laboratory Oversight Program, collected groundwater samples from 64 monitoring wells in the ESRP aquifer, 6 of which are part of a multilevel monitoring system, to determine the concentration of 129 I in the groundwater. These samples were analyzed by accelerator mass spectrometry as part of a long-term ongoing study to track trends and occurrences of this carcinogenic, long-lived radionuclide in the environment. Concentrations ranged from slightly above the locally determined background concentration of 5.4×10 −6 picocuries per liter, to just below the U.S. Environmental Protection Agency’s maximum contaminant level of 1 picocurie per liter. Discharge of wastewater containing 129 I has been discontinued to the aquifer, and long-term trends from a subset (n=15) of sampled wells show decreasing 129 I concentrations over the last three decades. Concentrations of 129 I in groundwater from monitoring wells near facilities at the Idaho National Laboratory are affected by episodic recharge from an ephemeral surface-water source and by the fracture-flow dominated hydrologic regime in the ESRP aquifer. The spatially focused sampling effort has also identified a low-level 129 I plume that affects long-term water quality near and downgradient from the Advanced Test Reactor Complex in the southwestern part of the facility that had not been clearly defined in previous sampling efforts, although the definition of the plume is somewhat limited by available data.

Idaho↗

USGS2024

Peak streamflow trends and their relation to changes in climate in Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin

Flood-frequency analysis, also called peak-flow frequency or flood-flow frequency analysis, is essential to water resources management applications including critical structure design and floodplain mapping. Federal guidelines for doing 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 is that for drainage basins without major hydrologic alterations, statistical properties of the distribution of annual peak streamflows (peak flows) are stationary; that is, the mean, variance, and skew are constant. The stationarity assumption has been widely accepted within the flood-frequency community; however, a better understanding of long-term climatic persistence and concerns about potential climate change and land-use change has caused a reexamination of the stationarity assumption. This work is part of that reexamination. The stationarity assumption is a concern because flood-frequency analyses that do not incorporate observed trends and abrupt changes may result in a poor representation of the true flood risk. Bulletin 17C does not offer guidance for incorporating nonstationarities when estimating floods, and it describes a need for studies that incorporate changing climate or basin characteristics. In response to this need and a history of concern regarding nonstationarity peak flows in the region, this study was done to assess potential nonstationarity in peak flows in the north-central United States. This report summarizes the methods used to detect hydroclimatic changes in peak-flow data in the study region. Four periods were selected for analysis of peak flow, daily streamflow, and climate data. The periods are (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. The starting point for these analyses was the initial data analysis of peak flow described in Bulletin 17C, which includes plotting the peak flow and checking for autocorrelation, monotonic trends, and changes points. Analyses were added to examine additional features in the data. Results are provided in a U.S. Geological Survey data release. The study limitations are documented for users of the results.

Scientific Investigations Report↗

USGS1999

Effects of snowmobile use on snowpack chemistry in Yellowstone National Park, 1998

Snowmobile use in Yellowstone National Park has increased substantially in the past three decades. In areas of greatest snowmobile use, elevated levels of by-products of gasoline combustion such as ammonium and benzene have been detected in snowpack samples. Annual snowpacks and snow-covered roadways trap deposition from local and regional atmospheric emissions. Snowpack samples representing most of the winter precipitation were collected at about the time of maximum annual snow accumulation at a variety of locations in the park to observe effects of a range of snowmobile traffic levels. Concentrations of organic and inorganic compounds in snow samples from pairs of sites located directly in and off snow-packed roadways used by snowmobiles were compared to concentrations in samples collected at nearby off-road sites. Concentrations of ammonium were 2 to 5 times higher for the in-road snow compared to off-road snow for each pair of sites. Thus, concentrations decreased rapidly with distance from roadways. In addition, concentrations of ammonium, nitrate, sulfate, benzene, and toluene in snow were positively correlated with snowmobile use.

Wyoming↗

USGS2025

Synergy between geology and geophysics in graphite mineral resource assessment

Graphite is designated as a critical mineral by the U.S. Government due to its essential role in modern technology and its vulnerability to supply chain disruption. To evaluate domestic graphite resources, the U.S. Geological Survey (USGS) conducted two airborne electromagnetic (AEM) surveys as part of the Earth Mapping Resources Initiative (Earth MRI) over prospective areas in the southern Appalachians of Alabama and the Seward Peninsula of Alaska. These surveys aimed to refine areas of interest by revealing subsurface electrical resistivity patterns indicative of graphite-bearing rocks. The results of the AEM surveys were integrated with existing geological knowledge and supplemented by field observations to assess the likelihood of undiscovered resources in these regions. This multidisciplinary approach combined geophysical techniques with traditional geological methods, enhancing the understanding of graphite potential in areas that may have been previously overlooked. The findings from this initiative are expected to contribute significantly to the ongoing efforts to secure a reliable domestic supply of graphite, thereby reducing import reliance and supporting national interests in mineral resource management.

Conference Paper↗

USGS2026

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

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

Iowa↗

USGS2025

Application of Hydrologic Simulation Program—FORTRAN (HSPF) as part of an integrated hydrologic model for the Salinas Valley, California

The U.S. Geological Survey (USGS), in cooperation with the Monterey County Water Resources Agency, conducted studies to help evaluate the surface-water and groundwater resources of the Salinas Valley study area, consisting of the entire Salinas River watershed and several smaller, adjacent coastal watersheds draining into Monterey Bay. The Salinas Valley study area is a highly productive agricultural region that depends on the coordinated use of surface water and groundwater to meet demand for irrigation and public water supply. To continue to meet these demands, a better understanding of the historical water balance and the effects of water-resource development on the long-term sustainability of water resources in the Salinas Valley study area is needed.

California↗