Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Modeling Report”

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.

At least 1,243 records · Page 69Linked to original sources
USGS2023

Panel review of Ground Motion Characterization Model in 2023 NSHM

The 2023 National Seismic Hazard Model (NSHM; Petersen et al., 2023) has two major components – a seismic source characterization (SSC) model and a ground motion characterization (GMC) model. The US Geological Survey (USGS) established separate panels to review and provide input on these two models. Both panels are advisory, meaning that they provide input on technical issues for consideration by the USGS NSHM team, but they do not have decision making authority. Here, we report on the activities and recommendations of the Ground Motion Characterization Panel, made up of the authors of this review. Final modeling decisions are presented in separate USGS documents, including Petersen et al., (2023) and Moschetti et al., (2023). Where modeling decisions depart from our recommendations, the rationale is explained in those publications.

Report↗
USGS2020

Characterization of peak streamflow and stages at selected streamgages in eastern and northeastern Oklahoma from the May to June 2019 flood event—With an emphasis on flood peaks downstream from dams and on tributaries to the Arkansas River

As much as 22 inches of rain fell in Oklahoma in May 2019, resulting in historic flooding along the Arkansas River and its tributaries in eastern and northeastern Oklahoma. The flooding along the Arkansas River and its tributaries that began in May continued into June 2019. Peaks of record were measured at nine U.S. Geological Survey (USGS) and U.S. Army Corps of Engineers (USACE) streamgages on various streams in eastern and northeastern Oklahoma. This report documents the peak streamflows and stages for 38 selected streamgages in eastern and northeastern Oklahoma and is a followup to a previous report by the USGS that documented flood peaks associated with the May 2019 flood event. Most of the flood peaks occurred from May 26 to June 4, 2019. This report includes data from streamgages on tributaries to the Arkansas River and uses modeling methods to extend the period of record for Arkansas River streamgages. The historic flooding caused homes to fall into the river as a result of bank erosion, forced some towns to be evacuated, and resulted in the highest flood depths in Tulsa, Oklahoma, since 1986. Several USGS and USACE streamgages along the Arkansas River and its tributaries recorded new peaks of record.

Oklahoma↗
USGS2020

Estimating rupture dimensions of three major earthquakes in Sichuan, China, for early warning and rapid loss estimates

Large earthquakes like in Wenchuan in 2008, MW 7.9, Sichuan, China, provide opportunity for earthquake early warning (EEW) as many heavily shaken areas are far (~50 km) from the epicenter and warning time could be long enough (≥ 5 s) to take effective preventative action. On the other hand, earthquakes with magnitudes larger than ~M 6.5 are challenging for EEW since source dimensions need to be defined in order to adequately estimate shaking. The Finite-Fault Rupture Detector (FinDer) is an approach to identify fault rupture extents from real-time strong motion and/or broadband records. In this study, we playback local and regional on-scale strong motion waveforms recorded during the 2008 MW 7.9 Wenchuan, 2013 MW 6.6 Lushan, and 2017 MW 6.5 Jiuzhaigou earthquakes to study the performance of FinDer for the current layout of the China Strong Motion Network. Overall, the FinDer line-source models agree well with the observed spatial distribution of aftershocks and fault models determined from waveform inversion. However, since FinDer models are constructed to characterize seismic ground motions (as needed for EEW) instead of source parameters, the rupture length can be overestimated for events radiating high levels of high-frequency motions, as is the case in the Lushan earthquake. If the set of strong motion data used had been available in real-time, 50% to 80% of sites experiencing shaking of intensity MMI IV-VII (light to very strong) and 30% experiencing VIII-IX (severe to violent) could have been issued a warning with 10 s and 5 s, respectively, before the arrival of the destructive S-wave. We also show that loss estimates after devastating earthquakes based on the FinDer line-source are more accurate compared to a point-source model. For the Wenchuan earthquake, for example, they predict a four to six times larger number of fatalities and injured, which is consistent with official reports. At the same time, these losses could be provided 1/2~3 hours faster than if based on more complex inversion rupture models.

Sichuan↗
USGS2025

Estimating agricultural irrigation water consumption for the High Plains aquifer region with integrated energy- and water-balance evapotranspiration modeling approaches

Estimation of irrigation water use provides essential information for the management and conservation of agricultural water resources. Conventionally, water use data are created based on reports and surveys from water users, whereas manual records may not be complete due to lacking flow meters, measurement gaps, inconsistent methods across regions, and time- and cost-consuming data processing. Alternatively, spatially explicit estimation of irrigation water use can be conducted efficiently using remote sensing evapotranspiration (ET) modeling approaches. In this study, we created a gridded blue water evapotranspiration (BWET) dataset to estimate historical irrigation water consumption (1986 – 2020) in the croplands across the United States High Plains aquifer region. The BWET data were generated by integrating an energy-balance ET model [Operational Simplified Surface Energy Balance model (SSEBop)] and a water-balance ET model [Vegetation ET model (VegET)]. BWET in croplands indicates crop consumptive use of irrigation water extracted from surface water and groundwater resources. The BWET estimates were compared with reported irrigation water use data for all counties within the aquifer region. The results revealed high agreement between growing season (May – September) BWET and annual water withdrawal at county level. Specifically, correlation coefficients of volumetric BWET and water withdrawal were 0.90 and 0.96, respectively, for the entire aquifer region and western Kansas. The timeseries of BWET and water withdrawal showed similar temporal trends and high covariations. The BWET estimates were systematically lower than the water withdrawal measurements, which was primarily attributed to blue water losses in the irrigation system. The irrigation efficiency, calculated as the ratio of BWET to water withdrawal depth, was 0.57 and 0.74 for the entire aquifer region and western Kansas, respectively. This study demonstrates the capability of using satellite-based ET models (e.g., SSEBop and VegET) to efficiently estimate crop water consumption and evaluate irrigation efficiency at landscape, county, and regional scales.

Colorado, Kansas, Nebraska, New Mexico, Oklahoma, ↗
USGS2014

Simulation of groundwater flow and interaction of groundwater and surface water on the Lac du Flambeau Reservation, Wisconsin

The Lac du Flambeau Band of Lake Superior Chippewa and Indian Health Service are interested in improving the understanding of groundwater flow and groundwater/surface-water interaction on the Lac du Flambeau Reservation (Reservation) in southwest Vilas County and southeast Iron County, Wisconsin, with particular interest in an understanding of the potential for contamination of groundwater supply wells and the fate of wastewater that is infiltrated from treatment lagoons on the Reservation. This report describes the construction, calibration, and application of a regional groundwater flow model used to simulate the shallow groundwater flow system of the Reservation and water-quality results for groundwater and surface-water samples collected near a system of waste-water-treatment lagoons. Groundwater flows through a permeable glacial aquifer that ranges in thickness from 60 to more than 200 feet (ft). Seepage and drainage lakes are common in the area and influence groundwater flow patterns on the Reservation. A two-dimensional, steady-state analytic element groundwater flow model was constructed using the program GFLOW. The model was calibrated by matching target water levels and stream base flows through the use of the parameter-estimation program, PEST. Simulated results illustrate that groundwater flow within most of the Reservation is toward the Bear River and the chain of lakes that feed the Bear River. Results of analyses of groundwater and surface-water samples collected downgradient from the wastewater infiltration lagoons show elevated levels of ammonia and dissolved phosphorus. In addition, wastewater indicator chemicals detected in three downgradient wells and a small downgradient stream indicate that infiltrated wastewater is moving southwest of the lagoons toward Moss Lake. Potential effects of extended wet and dry periods (within historical ranges) were evaluated by adjusting precipitation and groundwater recharge in the model and comparing the resulting simulated lake stage and water budgets to stages and water budgets from the calibrated model. Simulated lake water budgets and water level changes illustrate the importance of understanding the position of a lake within the hydrologic system (headwater or downstream), the type of lake (surface-water drainage or seepage lake), and the role of groundwater in dampening the effects of large-scale changes in weather patterns on lake levels. Areas contributing recharge to drinking-water supply wells on the Reservation were delineated using forward particle tracking from the water table to the well. Monte Carlo uncertainty analyses were used to produce maps showing the probability of groundwater capture for areas around each well nest. At the Main Pumphouse site near the Village of Lac du Flambeau, most of the area contributing recharge to the wells occurs downgradient from a large wetland between the wells and the wastewater infiltration lagoons. Nonetheless, a small potential for the wells to capture infiltrated wastewater is apparent when considering uncertainty in the model parameter values. At the West Pumphouse wells south of Flambeau Lake, most of the area contributing recharge is between the wells and Tippecanoe Lake. The extent of infiltrated wastewater from two infiltration lagoons was tracked using the groundwater flow model and Monte Carlo uncertainty analyses. Wastewater infiltrated from the lagoons flows predominantly south toward Moss Lake as it integrates with the regional groundwater flow system. The wastewater-plume-extent simulations support the area-contributing-recharge simulations, indicating that there is a possibility, albeit at low probability, that some wastewater could be captured by water-supply wells. Comparison of simulated water-table contours indicate that the lagoons may mound the water table approximately 4 ft, with diminishing levels of mounding outward from the lagoons. Four scenarios, representing potential alternatives for wastewater management, were simulated (at current discharge rates) to evaluate the potential extent of wastewater in the aquifer and discharge to surface-water bodies associated with each management scenario. Wastewater simulated to infiltrate through a hypothetical diffuser below a wetland south of the current lagoons appears to discharge to the overlying wetland and would likely discharge to Moss Lake as overland flow. Wastewater simulated to discharge to a small lake (Mindy Lake) between Moss and Fence Lakes appears to spread radically over a large area between the lakes. Wastewater simulated to discharge to lagoons south and northeast of the current lagoons also appears to spread radially, but the areas of the aquifer with the highest probability of encountering waste-water contamination would likely be between the lagoons and the nearest lake, where the wastewater would eventually discharge. Probability results for the wastewater-plume-extent scenarios are sensitive to the number of mathematical water particles used to represent infiltrating wastewater and the level of detail in the synthetic grid used for the probability analysis. Thus, probability results from wastewater-plume-extent simulations are qualitative only; however, it is expected that illustrations of relatively high or low probability will be useful as a general guide for decision making. Management problems requiring quantitative estimates of probability are best re-cast into problems evaluating the area that contributes recharge to the location of interest, which is not dependent upon the number of simulated particles or the resolution of a synthetic grid.

Wisconsin↗
USGS2019

Regionalization of surface-water statistics using multiple linear regression

This report serves as a reference document in support of the regionalization of surface-water statistics using multiple linear regression. Streamflow statistics are quantitative characterizations of hydrology and are often derived from observed streamflow records. In the absence of observed streamflow records, as at unmonitored or ungaged locations, other techniques are required. Multiple linear regression is one tool that is widely used to regionalize or transfer information from gaged to ungaged locations. This report provides the background to support regression-based regionalization of streamflow statistics. This background includes tools for data assembly, exploratory data analysis, model estimation in a least-squares framework, and model evaluation.

Techniques and Methods↗
USGS2006

Scaling local species-habitat relations to the larger landscape with a hierarchical spatial count model

Much of what is known about avian species-habitat relations has been derived from studies of birds at local scales. It is entirely unclear whether the relations observed at these scales translate to the larger landscape in a predictable linear fashion. We derived habitat models and mapped predicted abundances for three forest bird species of eastern North America using bird counts, environmental variables, and hierarchical models applied at three spatial scales. Our purpose was to understand habitat associations at multiple spatial scales and create predictive abundance maps for purposes of conservation planning at a landscape scale given the constraint that the variables used in this exercise were derived from local-level studies. Our models indicated a substantial influence of landscape context for all species, many of which were counter to reported associations at finer spatial extents. We found land cover composition provided the greatest contribution to the relative explained variance in counts for all three species; spatial structure was second in importance. No single spatial scale dominated any model, indicating that these species are responding to factors at multiple spatial scales. For purposes of conservation planning, areas of predicted high abundance should be investigated to evaluate the conservation potential of the landscape in their general vicinity. In addition, the models and spatial patterns of abundance among species suggest locations where conservation actions may benefit more than one species. ?? 2006 Springer Science+Business Media B.V.

Landscape Ecology↗
USGS2020

Carrying capacity of a population diffusing in a heterogeneous environment

The carrying capacity of the environment for a population is one of the key concepts in ecology and it is incorporated in the growth term of reaction-diffusion equations describing populations in space. Analysis of reaction-diffusion models of populations in heterogeneous space have shown that, when the maximum growth rate and carrying capacity in a logistic growth function vary in space, conditions exist for which the total population size at equilibrium (i) exceeds the total population that which would occur in the absence of diffusion and (ii) exceeds that which would occur if the system were homogeneous and the total carrying capacity, computed as the integral over the local carrying capacities, was the same in the heterogeneous and homogeneous cases. We review here work over the past few years that has explained these apparently counter-intuitive results in terms of the way input of energy or another limiting resource (e.g., a nutrient) varies across the system. We report on both mathematical analysis and laboratory experiments confirming that total population size in a heterogeneous system with diffusion can exceed that in the system without diffusion. We further report, however, that when the resource of the population in question is explicitly modeled as a coupled variable, as in a reaction-diffusion chemostat model rather than a model with logistic growth, the total population in the heterogeneous system with diffusion cannot exceed the total population size in the corresponding homogeneous system in which the total carrying capacities are the same.

Mathematics and Computers in Modern Science - Acou↗
USGS2007

A coupled remote sensing and simplified surface energy balance approach to estimate actual evapotranspiration from irrigated fields

Accurate crop performance monitoring and production estimation are critical for timely assessment of the food balance of several countries in the world. Since 2001, the Famine Early Warning Systems Network (FEWS NET) has been monitoring crop performance and relative production using satellite-derived data and simulation models in Africa, Central America, and Afghanistan where ground-based monitoring is limited because of a scarcity of weather stations. The commonly used crop monitoring models are based on a crop water-balance algorithm with inputs from satellite-derived rainfall estimates. These models are useful to monitor rainfed agriculture, but they are ineffective for irrigated areas. This study focused on Afghanistan, where over 80 percent of agricultural production comes from irrigated lands. We developed and implemented a Simplified Surface Energy Balance (SSEB) model to monitor and assess the performance of irrigated agriculture in Afghanistan using a combination of 1-km thermal data and 250m Normalized Difference Vegetation Index (NDVI) data, both from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor. We estimated seasonal actual evapotranspiration (ETa) over a period of six years (2000-2005) for two major irrigated river basins in Afghanistan, the Kabul and the Helmand, by analyzing up to 19 cloud-free thermal and NDVI images from each year. These seasonal ETa estimates were used as relative indicators of year-to-year production magnitude differences. The temporal water-use pattern of the two irrigated basins was indicative of the cropping patterns specific to each region. Our results were comparable to field reports and to estimates based on watershed-wide crop water-balance model results. For example, both methods found that the 2003 seasonal ETa was the highest of all six years. The method also captured water management scenarios where a unique year-to-year variability was identified in addition to water-use differences between upstream and downstream basins. A major advantage of the energy-balance approach is that it can be used to quantify spatial extent of irrigated fields and their water-use dynamics without reference to source of water as opposed to a water-balance model which requires knowledge of both the magnitude and temporal distribution of rainfall and irrigation applied to fields. ?? 2007 by MDPI.

Sensors↗
USGS2015

U.S. Geological Survey Noble Gas Laboratory’s standard operating procedures for the measurement of dissolved gas in water samples

This report addresses the standard operating procedures used by the U.S. Geological Survey’s Noble Gas Laboratory in Denver, Colorado, U.S.A., for the measurement of dissolved gases (methane, nitrogen, oxygen, and carbon dioxide) and noble gas isotopes (helium-3, helium-4, neon-20, neon-21, neon-22, argon-36, argon-38, argon-40, kryton-84, krypton-86, xenon-103, and xenon-132) dissolved in water. A synopsis of the instrumentation used, procedures followed, calibration practices, standards used, and a quality assurance and quality control program is presented. The report outlines the day-to-day operation of the Residual Gas Analyzer Model 200, Mass Analyzer Products Model 215–50, and ultralow vacuum extraction line along with the sample handling procedures, noble gas extraction and purification, instrument measurement procedures, instrumental data acquisition, and calculations for the conversion of raw data from the mass spectrometer into noble gas concentrations per unit mass of water analyzed. Techniques for the preparation of artificial dissolved gas standards are detailed and coupled to a quality assurance and quality control program to present the accuracy of the procedures used in the laboratory.

Techniques and Methods↗
USGS1998

Hydrogeologic framework of the Puget Sound aquifer system, Washington and British Columbia

This report presents the generalized hydrogeologic framework of the Puget Sound aquifer system in Washington and British Columbia. The framework includes a conceptual model of the division of the aquifer system into regional hydrogeologic units for describing on a regional basis the ground-water flow in the system. The conceptual model is based on an analysis of historical data and on results of cross-sectional numerical models of ground-water flow.

British Columbia, Washington↗
USGS1984

An assessment of cumulative impacts of coal mining on the hydrology in part of the Powder River structural basin, Wyoming: A progress report

The U.S. Geological Survey and the Wyoming Department of Environmental Quality are involved in a cooperative effort to assess the probable cumulative impacts of coal mining on the hydrology of a part of the Powder River Structural Basin in Wyoming. It was assumed that the principal impacts on the ground-water system due to mining will occur in the relatively shallow aquifers which can be grouped into three homogeneous aquifers, namely, the Wyodak coal, the overburden, and the under burden. Emphasis of this report is on the results of analysis of surface-water resources in the Caballo Creek drainage. A surface-water model of the Caballo Creek drainage was developed using the Hydrological Simulation Program-Fortran model to help assess the impacts of mining activities on streamflow. The Caballo Creek drainage was divided into 10 land segments and 6 stream reaches in the modeling process. Three simulation runs show little, if any, change in streamflow between pre- and post-mining conditions and very little change between pre-mining and during-mining conditions. The principal reason for the absence of change is the high infiltration rate used in the model for all three conditions.

Wyoming↗
USGS2018

Golden Eagle Monitoring Plan for the Desert Renewable Energy Conservation Plan

This report describes options for monitoring the status and population trends of the golden eagle (Aquila chrysaetos) within the Desert Renewable Energy Conservation Plan (DRECP) area of Southern California in maintaining stable or increasing population in the planning area. The report profiles the ecology of golden eagles in the region and provides a range of potential sampling options to address monitoring needs and objectives. This approach also focused on links between changes in human land-use, golden eagle nesting and foraging habitat conditions, and population dynamics. The report outlines how monitoring data from demographic, prey, and habitat studies were used to develop a predictive demographic model for golden eagles in the DRECP area. Results from the model simulations suggest increases in renewable energy development could have negative consequences for population trajectories. Results also suggest site-specific conservation actions could reduce the magnitude of negative impacts to the local population of eagles. A monitoring framework is proposed including: (1) annual assessments of site-occupancy and reproduction by territorial pairs of golden eagles (including rates at which sites become colonized or vacated over time); (2) estimates of survival, movements, and intensity of use of landscapes by breeding and non-breeding golden eagles; (3) periodic (conducted every two to four years) assessments of nesting and foraging habitats, prey populations, and associations with land-use and management activities; and (4) updating the predictive demographic model with new information obtained on eagles and associated population stressors. The results of this research were published in the Journal of Rapture Research, Wiens, David,Inman, Rich D., Esque, Todd C., Longshore, Kathleen M. and Nussear, Kenneth (2017). Spatial Demographic Models to Inform Conservation Planning of Golden Eagles in Renewable Energy Landscapes. 51(3):234-257.

Report↗
USGS2003

Use of boundary fluxes when simulating solute transport with the MODFLOW ground-water transport process

This report describes modifications to a U.S. Geological Survey (USGS) threedimensional solute-transport model (MODFLOWGWT), which is incorporated into the USGS MODFLOW ground-water model as the Ground- Water Transport (GWT) Process. The modifications improve the capability of MODFLOW-GWT to accurately simulate solute transport in simulations that represent a nonzero flux across an aquifer boundary. In such situations, the new Boundary Flux Package (BFLX) will allow the user flexibility to assign the flux to specific cell faces, although that flexibility is limited for certain types of fluxes (such as recharge and evapotranspiration, which can only be assigned to the top face if either is to be represented as a boundary flux). The approach is consistent with that used in the MODPATH model. The application of the BFLX Package was illustrated using a test case in which the Lake Package was active. The results using the BFLX Package showed noticeably higher magnitudes of velocity in the cells adjacent to the lake than previous results without the BFLX Package. Consequently, solute was transported slightly faster through the lake-aquifer system when the BFLX Package is active. However, the overall solute distributions did not differ greatly from simulations made without using the BFLX Package.

Open-File Report↗
USGS2011

Rapid estimation of the economic consequences of global earthquakes

The U.S. Geological Survey's (USGS) Prompt Assessment of Global Earthquakes for Response (PAGER) system, operational since mid 2007, rapidly estimates the most affected locations and the population exposure at different levels of shaking intensities. The PAGER system has significantly improved the way aid agencies determine the scale of response needed in the aftermath of an earthquake. For example, the PAGER exposure estimates provided reasonably accurate assessments of the scale and spatial extent of the damage and losses following the 2008 Wenchuan earthquake (Mw 7.9) in China, the 2009 L'Aquila earthquake (Mw 6.3) in Italy, the 2010 Haiti earthquake (Mw 7.0), and the 2010 Chile earthquake (Mw 8.8). Nevertheless, some engineering and seismological expertise is often required to digest PAGER's exposure estimate and turn it into estimated fatalities and economic losses. This has been the focus of PAGER's most recent development. With the new loss-estimation component of the PAGER system it is now possible to produce rapid estimation of expected fatalities for global earthquakes (Jaiswal and others, 2009). While an estimate of earthquake fatalities is a fundamental indicator of potential human consequences in developing countries (for example, Iran, Pakistan, Haiti, Peru, and many others), economic consequences often drive the responses in much of the developed world (for example, New Zealand, the United States, and Chile), where the improved structural behavior of seismically resistant buildings significantly reduces earthquake casualties. Rapid availability of estimates of both fatalities and economic losses can be a valuable resource. The total time needed to determine the actual scope of an earthquake disaster and to respond effectively varies from country to country. It can take days or sometimes weeks before the damage and consequences of a disaster can be understood both socially and economically. The objective of the U.S. Geological Survey's PAGER system is to reduce this time gap to more rapidly and effectively mobilize response. We present here a procedure to rapidly and approximately ascertain the economic impact immediately following a large earthquake anywhere in the world. In principle, the approach presented is similar to the empirical fatality estimation methodology proposed and implemented by Jaiswal and others (2009). In order to estimate economic losses, we need an assessment of the economic exposure at various levels of shaking intensity. The economic value of all the physical assets exposed at different locations in a given area is generally not known and extremely difficult to compile at a global scale. In the absence of such a dataset, we first estimate the total Gross Domestic Product (GDP) exposed at each shaking intensity by multiplying the per-capita GDP of the country by the total population exposed at that shaking intensity level. We then scale the total GDP estimated at each intensity by an exposure correction factor, which is a multiplying factor to account for the disparity between wealth and/or economic assets to the annual GDP. The economic exposure obtained using this procedure is thus a proxy estimate for the economic value of the actual inventory that is exposed to the earthquake. The economic loss ratio, defined in terms of a country-specific lognormal cumulative distribution function of shaking intensity, is derived and calibrated against the losses from past earthquakes. This report describes the development of a country or region-specific economic loss ratio model using economic loss data available for global earthquakes from 1980 to 2007. The proposed model is a potential candidate for directly estimating economic losses within the currently-operating PAGER system. PAGER's other loss models use indirect methods that require substantially more data (such as building/asset inventories, vulnerabilities, and the asset values exposed at the time of earthquake) to implement on a global basis and will thus take more time to develop and implement within the PAGER system.

Open-File Report↗
USGS2019

Flood-frequency estimates for Ohio streamgages based on data through water year 2015 and techniques for estimating flood-frequency characteristics of rural, unregulated Ohio streams

Estimates of the magnitudes of annual peak streamflows with annual exceedance probabilities of 0.5, 0.2, 0.1, 0.04, 0.02, 0.01, and 0.002 (equivalent to recurrence intervals of 2-, 5-, 10-, 25-, 50-, 100-, and 500-years, respectively) were computed for 391 streamgages in Ohio and adjacent states based on data collected through the 2015 water year. The flood-frequency estimates were computed following guidance outlined in Bulletin 17C, developed by the Advisory Committee on Water Information. The Bulletin 17C guidelines retain the basic statistical framework of the superseded Bulletin 17B guidelines; however, the Bulletin 17C guidelines add several enhancements including an improved method of moments approach for fitting the log-Pearson Type III (LPIII) distribution to the flood peaks (called the expected moments algorithm), a generalization of the Grubbs Beck low-outlier test (called the Multiple Grubbs Beck test) that permits identification of multiple potentially influential low floods, and new methods for estimating regional skew and uncertainty. Equations for estimating flood-frequency characteristics at ungaged sites on rural, unregulated streams in Ohio were developed with a two-step process involving ordinary least-squares and generalized least-squares regression techniques. Data from 333 streamgages with 10 or more years of unregulated record were screened for redundancy and a regression dataset was selected that was composed of flood-frequency and basin-characteristic data for 275 streamgages in Ohio and adjacent states. Two sets of equations were developed—one set, referred to as the “simple model,” uses regression region and drainage area as regressor variables, and a second set, referred to as the “full model,” uses regression region, drainage area, main-channel slope, and the percentage of the watershed covered by water and wetlands as regressor variables. The average standard errors of prediction ranged from about 40.5 to 46.5 percent for the simple-model equations and from about 37.2 to 40.3 percent for the full-model equations. For sites meeting the rural, unregulated criteria, flood-frequency estimates determined by means of LPIII analyses are reported along with weighted flood-frequency estimates, computed as a function of the LPIII estimates and the regression estimates. For sites with homogenous periods of regulation, flood-frequency estimates determined by means of LPIII analyses are reported. Ninety-five percent confidence limits are reported for all estimates. Values of regressor variables were determined from digital spatial datasets by means of a geographic information system (GIS). The GIS datasets and the new full-model equations have been incorporated into Ohio’s StreamStats application, a web-based, GIS-backed system designed to facilitate the estimation of streamflow statistics at ungaged locations on streams. Seasonal patterns in peak flows were assessed for 295 streamgages in Ohio. Annual peak flows occurred most frequently between January and April, with March having the highest frequency of occurrence. The month with the fewest number of annual peaks was October. Peak-of-record flows occurred most frequently in March, followed by January (months in which two of Ohio’s most severe widespread floods in recent history occurred). None of the peak-of-record flows occurred in October and only two occurred in November. Temporal trend in annual peak flows were assessed for 133 streamgages on unregulated streams in Ohio with 30 or more years of systematic record. Trends were assessed by computing the rank correlation (as measured with the two-sided Kendall’s tau statistic) between time and annual peak flows. Weak but statistically significant trends were indicated at 15 of the 133 streamgages. Of the 15 streamgages with significant trend in annual peak flows, 12 had an upward trend (positive tau) and 3 had a downward trend (negative tau). All 12 streamgages with positive tau values were at latitudes north of 40°33', and streamgages with negative tau values were at latitudes south of 40°33'.

Ohio↗
USGS1995

User's guide to PHREEQC, a computer program for speciation, reaction-path, advective-transport, and inverse geochemical calculations

PHREEQC is a computer program written in the C programming language that is designed to perform a wide variety of aqueous geochemical calculations. PHREEQC is based on an ion-association aqueous model and has capabilities for (1) speciation and saturation-index calculations, (2) reaction-path and advective-transport calculations involving specified irreversible reactions, mixing of solutions, mineral and gas equilibria surface-complex-ation reactions, and ion-exchange reactions, and (3) inverse modeling, which finds sets of mineral and gas mole transfers that account for composition differences between waters, within specified compositional uncertainties. PHREEQC is derived from the Fortran program PHREEQE, but it has been completely rewritten in C with the addition many new capabilities. New features include the capabilities to use redox couples to distribute redox elements among their valence states in speciation calculations; to model ion-exchange and surface-complex-ation reactions; to model reactions with a fixed-pressure, multicomponent gas phase (that is, a gas bubble); to calculate the mass of water in the aqueous phase during reaction and transport calculations; to keep track of the moles of minerals present in the solid phases and determine automatically the thermodynamically stable phase assemblage; to simulate advective transport in combination with PHREEQC's reaction-modeling capability; and to make inverse modeling calculations that allow for uncertainties in the analytical data. The user interface is improved through the use of a simplified approach to redox reactions, which includes explicit mole-balance equations for hydrogen and oxygen; the use of a revised input that is modular and completely free format; and the use of mineral names and standard chemical symbolism rather than index numbers. The use of (2) eliminates nearly all limitations on army sizes, including numbers of elements, aqueous species, solutions, phases, and lengths of character strings. A new equation solver that optimizes a set of equalities subject to both equality and inequality constraints is used to determine the thermodynamically stable set of phases in equilibrium with a solution. A more complete Newton-Raphson formulation, master-species switching, and scaling of the algebraic equations reduce the number of failures of the numerical method in PHREEQC relative to PHREEQE. This report presents the equations that are the basis for chemical equilibrium and inverse-modeling calculations in PHREEQC, describes the input for the program, and presents twelve examples that demonstrate most of the program's capabilities.

Water-Resources Investigations Report↗
USGS2026

Working group on American Eel (WGAMEEL; outputs from 2024 meeting)

The Working Group on American Eel (WGAMEEL) met virtually three times in 2022-2024 to address the five Terms of Reference (ToRs) of its three-year term. The first two ToRs tasked WGAMEEL with listing and evaluating data on American eel landings, abundance indices, and spatial and habitat data and also to describe assessment methods used in Canada and the US. Canada-wide American eel trajectory was estimated using 12 fishery-independent datasets. Generally, longer datasets had more negative trends than shorter ones. Limiting data to the post-2000 years produced fewer negative trends that did not differ from zero, suggesting the observed declines occurred pre-2000. Spatial modelling for American eel requires knowledge and mapping that covers the breadth of habitat types occupied by the species, including freshwater, estuarine, and marine environments. In recent years there has been an expansion of online databases with data from the aquatic environment, particularly in freshwater, with estuarine and marine data less consistently documented. This report broadly compiles abiotic data series of relevance to American eel. A larger challenge for spatial modelling will be acquiring enough high quality, georeferenced biological data sets with suitable observations to assess occurrence, abundance, and trends over time in a spatial framework. The third ToR considered Indigenous Knowledge Systems for American eel. A survey reaching First Nations representatives from four Canadian provinces confirmed the cultural importance of eels in Indigenous communities, and that Indigenous knowledge possessed by the participant groups are place-based and contextual, especially regarding threats impacting eels. The final two ToRs focused on identifying stock assessment modelling approaches applicable to American and European eel, and assessing whether any of these approaches might be appropriate for American eel management moving forward. WGAMEEL evaluated the various approaches for assessing American eel or providing management advice. Two approaches that could be completed in the next few years because of their minimal data needs are index-based methods and catch-only method. A suite of approaches considered by WGAMEEL that would take more time and data were spatial or habitat models, management strategy evaluation, and spawner-per-recruit (SPR) models potentially paired with meta-population models.

ICES Scientific Reports↗