SEARCH · USGS Science

Results for “Environmental Modeling and Software”

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.

At least 145 records · Page 8Linked to original sources

USGS2006

Hierarchical models of animal abundance and occurrence

Much of animal ecology is devoted to studies of abundance and occurrence of species, based on surveys of spatially referenced sample units. These surveys frequently yield sparse counts that are contaminated by imperfect detection, making direct inference about abundance or occurrence based on observational data infeasible. This article describes a flexible hierarchical modeling framework for estimation and inference about animal abundance and occurrence from survey data that are subject to imperfect detection. Within this framework, we specify models of abundance and detectability of animals at the level of the local populations defined by the sample units. Information at the level of the local population is aggregated by specifying models that describe variation in abundance and detection among sites. We describe likelihood-based and Bayesian methods for estimation and inference under the resulting hierarchical model. We provide two examples of the application of hierarchical models to animal survey data, the first based on removal counts of stream fish and the second based on avian quadrat counts. For both examples, we provide a Bayesian analysis of the models using the software WinBUGS.

Journal of Agricultural, Biological, and Environme↗

USGS2016

Using inferential sensors for quality control of Everglades Depth Estimation Network water-level data

The Everglades Depth Estimation Network (EDEN), with over 240 real-time gaging stations, provides hydrologic data for freshwater and tidal areas of the Everglades. These data are used to generate daily water-level and water-depth maps of the Everglades that are used to assess biotic responses to hydrologic change resulting from the U.S. Army Corps of Engineers Comprehensive Everglades Restoration Plan. The generation of EDEN daily water-level and water-depth maps is dependent on high quality real-time data from water-level stations. Real-time data are automatically checked for outliers by assigning minimum and maximum thresholds for each station. Small errors in the real-time data, such as gradual drift of malfunctioning pressure transducers, are more difficult to immediately identify with visual inspection of time-series plots and may only be identified during on-site inspections of the stations. Correcting these small errors in the data often is time consuming and water-level data may not be finalized for several months. To provide daily water-level and water-depth maps on a near real-time basis, EDEN needed an automated process to identify errors in water-level data and to provide estimates for missing or erroneous water-level data. The Automated Data Assurance and Management (ADAM) software uses inferential sensor technology often used in industrial applications. Rather than installing a redundant sensor to measure a process, such as an additional water-level station, inferential sensors, or virtual sensors, were developed for each station that make accurate estimates of the process measured by the hard sensor (water-level gaging station). The inferential sensors in the ADAM software are empirical models that use inputs from one or more proximal stations. The advantage of ADAM is that it provides a redundant signal to the sensor in the field without the environmental threats associated with field conditions at stations (flood or hurricane, for example). In the event that a station does malfunction, ADAM provides an accurate estimate for the period of missing data. The ADAM software also is used in the quality assurance and quality control of the data. The virtual signals are compared to the real-time data, and if the difference between the two signals exceeds a certain tolerance, corrective action to the data and (or) the gaging station can be taken. The ADAM software is automated so that, each morning, the real-time EDEN data are compared to the inferential sensor signals and digital reports highlighting potential erroneous real-time data are generated for appropriate support personnel. The development and application of inferential sensors is easily transferable to other real-time hydrologic monitoring networks.

Florida↗

USGS2007

Development of the Hydroecological Integrity Assessment Process for Determining Environmental Flows for New Jersey Streams

The natural flow regime paradigm and parallel stream ecological concepts and theories have established the benefits of maintaining or restoring the full range of natural hydrologic variation for physiochemical processes, biodiversity, and the evolutionary potential of aquatic and riparian communities. A synthesis of recent advances in hydroecological research coupled with stream classification has resulted in a new process to determine environmental flows and assess hydrologic alteration. This process has national and international applicability. It allows classification of streams into hydrologic stream classes and identification of a set of non-redundant and ecologically relevant hydrologic indices for 10 critical sub-components of flow. Three computer programs have been developed for implementing the Hydroecological Integrity Assessment Process (HIP): (1) the Hydrologic Indices Tool (HIT), which calculates 171 ecologically relevant hydrologic indices on the basis of daily-flow and peak-flow stream-gage data; (2) the New Jersey Hydrologic Assessment Tool (NJHAT), which can be used to establish a hydrologic baseline period, provide options for setting baseline environmental-flow standards, and compare past and proposed streamflow alterations; and (3) the New Jersey Stream Classification Tool (NJSCT), designed for placing unclassified streams into pre-defined stream classes. Biological and multivariate response models including principal-component, cluster, and discriminant-function analyses aided in the development of software and implementation of the HIP for New Jersey. A pilot effort is currently underway by the New Jersey Department of Environmental Protection in which the HIP is being used to evaluate the effects of past and proposed surface-water use, ground-water extraction, and land-use changes on stream ecosystems while determining the most effective way to integrate the process into ongoing regulatory programs. Ultimately, this scientifically defensible process will help to quantify the effects of anthropogenic changes and development on hydrologic variability and help planners and resource managers balance current and future water requirements with ecological needs.

Scientific Investigations Report↗

USGS2007

Prediction of lethal/effective concentration/dose in the presence of multiple auxiliary covariates and components of variance

Predictors of the percentile lethal/effective concentration/dose are commonly used measures of efficacy and toxicity. Typically such quantal-response predictors (e.g., the exposure required to kill 50% of some population) are estimated from simple bioassays wherein organisms are exposed to a gradient of several concentrations of a single agent. The toxicity of an agent may be influenced by auxiliary covariates, however, and more complicated experimental designs may introduce multiple variance components. Prediction methods lag examples of those cases. A conventional two-stage approach consists of multiple bivariate predictions of, say, medial lethal concentration followed by regression of those predictions on the auxiliary covariates. We propose a more effective and parsimonious class of generalized nonlinear mixed-effects models for prediction of lethal/effective dose/concentration from auxiliary covariates. We demonstrate examples using data from a study regarding the effects of pH and additions of variable quantities 2???,5???-dichloro-4???- nitrosalicylanilide (niclosamide) on the toxicity of 3-trifluoromethyl-4- nitrophenol to larval sea lamprey (Petromyzon marinus). The new models yielded unbiased predictions and root-mean-squared errors (RMSEs) of prediction for the exposure required to kill 50 and 99.9% of some population that were 29 to 82% smaller, respectively, than those from the conventional two-stage procedure. The model class is flexible and easily implemented using commonly available software. ?? 2007 SETAC.

Environmental Toxicology and Chemistry↗

USGS2019

Evaluating associations between environmental variables and Escherichia coli levels for predictive modeling at Pawtuckaway Beach in Nottingham, New Hampshire, from 2015 to 2017

From 2015 through 2017, the U.S. Geological Survey in cooperation with the New Hampshire Department of Health and Human Services and the New Hampshire Department of Environmental Services studied occurrences of high levels of Escherichia coli ( E. coli ) bacteria at the Pawtuckaway State Park Beach in Nottingham, New Hampshire. Historic data collected by the New Hampshire Department of Environmental Services indicated that E. coli concentrations in the water typically increased through the beach season to levels considered potentially harmful to beachgoers. During the three beach seasons that were studied, E. coli samples were collected three to four times per week, and water-quality and meteorological data were collected continuously. The Virtual Beach software was used to generate a predictive model for each year of the study (2015–2017), and the model for each of these years was tested with data from the other two. Additionally, data from all study years were combined to generate a comprehensive model to help identify independent variables that might characterize environmental conditions relative to E. coli levels during multiple seasons. The accuracy of the models in predicting the occurrence of high E. coli levels was marginal, but the models did provide insights into the likely mechanisms for increased E. coli levels during the seasons. Variables most important in explaining high E. coli levels were the presence of geese at the beach, the progression of the season, the number of visitors at the beach, and wind vectors relative to beach orientation.

New Hampshire↗

USGS2015

Social Values for Ecosystem Services, version 3.0 (SolVES 3.0): documentation and user manual

The geographic information system (GIS) tool, S ocial V alues for E cosystem S ervices (SolVES), was developed to incorporate quantified and spatially explicit measures of social values into ecosystem service assessments. SolVES 3.0 continues to extend the functionality of SolVES, which was designed to assess, map, and quantify the social values of ecosystem services. Social values—the perceived, nonmarket values the public ascribes to ecosystem services, particularly cultural services, such as aesthetics and recreation—can be evaluated for various stakeholder groups. These groups are distinguishable by their attitudes and preferences regarding public uses, such as motorized recreation and logging. As with previous versions, SolVES 3.0 derives a quantitative 10-point, social-values metric—the value index—from a combination of spatial and nonspatial responses to public value and preference surveys. The tool also calculates metrics characterizing the underlying environment, such as average distance to water and dominant landcover. SolVES 3.0 is integrated with Maxent maximum entropy modeling software to generate more complete social-value maps and offer robust statistical models describing the relationship between the value index and explanatory environmental variables. A model’s goodness of fit to a primary study area and its potential performance in transferring social values to similar areas using value-transfer methodology can be evaluated. SolVES 3.0 provides an improved public-domain tool for decision makers and researchers to evaluate the social values of ecosystem services and to facilitate discussions among diverse stakeholders regarding the tradeoffs among ecosystem services in a variety of physical and social contexts ranging from forest and rangeland to coastal and marine.

Open-File Report↗

USGS2020

Social Values for Ecosystem Services, version 4.0 (SolVES 4.0)—Documentation and user manual

The geographic information system tool, Social Values for Ecosystem Services (SolVES), was developed to incorporate quantified and spatially explicit measures of social values into ecosystem service assessments. SolVES 4.0 provides an open-source version of SolVES, which was designed to assess, map, and quantify the social values of ecosystem services. Social values—the perceived, nonmarket values the public ascribes to ecosystem services, particularly cultural services, such as aesthetics and recreation—can be evaluated for various stakeholder groups. These groups are distinguishable by factors such as their attitudes and preferences regarding public uses (for example, motorized recreation and logging). As with previous versions, SolVES 4.0 derives a quantitative 10-point, social-values metric—the value index—from a combination of spatial and nonspatial responses to public value and preference surveys. The tool also calculates metrics characterizing the underlying environment, such as average distance to water and dominant landcover. SolVES 4.0 has been developed with Python using a QGIS user interface and a PostgreSQL database for required data. SolVES is integrated with Maxent maximum entropy modeling software to generate more complete social-value maps and offer robust statistical models describing the relation between the value index and explanatory environmental variables. A model’s goodness of fit to a primary study area and its potential performance in transferring social values to similar areas using value-transfer methods can be evaluated. SolVES 4.0 provides an improved open-source, public-domain tool for decision makers and researchers to evaluate the social values of ecosystem services and to facilitate discussions among diverse stakeholders regarding the tradeoffs among ecosystem services in a variety of biophysical and social contexts including mountain, forest, coastal, riparian, agricultural, and urban environments around the globe.

Techniques and Methods↗

USGS2012

Social values for ecosystem services (SolVES): Documentation and user manual, version 2.0

In response to the need for incorporating quantified and spatially explicit measures of social values into ecosystem services assessments, the Rocky Mountain Geographic Science Center (RMGSC), in collaboration with Colorado State University, developed a geographic information system (GIS) application, Social Values for Ecosystem Services (SolVES). With version 2.0 (SolVES 2.0), RMGSC has improved and extended the functionality of SolVES, which was designed to assess, map, and quantify the perceived social values of ecosystem services. Social values such as aesthetics, biodiversity, and recreation can be evaluated for various stakeholder groups as distinguished by their attitudes and preferences regarding public uses, such as motorized recreation and logging. As with the previous version, SolVES 2.0 derives a quantitative, 10-point, social-values metric, the Value Index, from a combination of spatial and nonspatial responses to public attitude and preference surveys and calculates metrics characterizing the underlying environment, such as average distance to water and dominant landcover. Additionally, SolVES 2.0 integrates Maxent maximum entropy modeling software to generate more complete social value maps and to produce robust statistical models describing the relationship between the social values maps and explanatory environmental variables. The performance of these models can be evaluated for a primary study area, as well as for similar areas where primary survey data are not available but where social value mapping could potentially be completed using value-transfer methodology. SolVES 2.0 also introduces the flexibility for users to define their own social values and public uses, model any number and type of environmental variable, and modify the spatial resolution of analysis. With these enhancements, SolVES 2.0 provides an improved public domain tool for decisionmakers and researchers to evaluate the social values of ecosystem services and to facilitate discussions among diverse stakeholders regarding the tradeoffs among different ecosystem services in a variety of physical and social contexts ranging from forest and rangeland to coastal and marine.

Open-File Report↗

USGS2014

Virtual Beach 3: user's guide

Virtual Beach version 3 (VB 3 ) is a decision support tool that constructs site-specific statistical models to predict fecal indicator bacteria (FIB) concentrations at recreational beaches. VB 3 is primarily designed for beach managers responsible for making decisions regarding beach closures or the issuance of swimming advisories due to pathogen contamination. However, researchers, scientists, engineers, and students interested in studying relationships between water quality indicators and ambient environmental conditions will find VB 3 useful. VB 3 reads input data from a text file or Excel document, assists the user in preparing the data for analysis, enables automated model selection using a wide array of possible model evaluation criteria, and provides predictions using a chosen model parameterized with new data. With an integrated mapping component to determine the geographic orientation of the beach, the software can automatically decompose wind/current/wave speed and magnitude information into along-shore and onshore/offshore components for use in subsequent analyses. Data can be examined using simple scatter plots to evaluate relationships between the response and independent variables (IVs). VB 3 can produce interaction terms between the primary IVs, and it can also test an array of transformations to maximize the linearity of the relationship The software includes search routines for finding the "best" models from an array of possible choices. Automated censoring of statistical models with highly correlated IVs occurs during the selection process. Models can be constructed either using previously collected data or forecasted environmental information. VB 3 has residual diagnostics for regression models, including automated outlier identification and removal using DFFITs or Cook's Distances.

Report↗

USGS2012

Exploring similarities among many species distributions

Collecting species presence data and then building models to predict species distribution has been long practiced in the field of ecology for the purpose of improving our understanding of species relationships with each other and with the environment. Due to limitations of computing power as well as limited means of using modeling software on HPC facilities, past species distribution studies have been unable to fully explore diverse data sets. We build a system that can, for the first time to our knowledge, leverage HPC to support effective exploration of species similarities in distribution as well as their dependencies on common environmental conditions. Our system can also compute and reveal uncertainties in the modeling results enabling domain experts to make informed judgments about the data. Our work was motivated by and centered around data collection efforts within the Great Smoky Mountains National Park that date back to the 1940s. Our findings present new research opportunities in ecology and produce actionable field-work items for biodiversity management personnel to include in their planning of daily management activities.

Conference Paper↗

USGS2019

Let’s agree to disagree: Comparing auto-acoustic identification programs for northeastern bats

With the declines in abundance and changing distribution of white-nose syndrome–affected bat species, increased reliance on acoustic monitoring is now the new “normal.” As such, the ability to accurately identify individual bat species with acoustic identification programs has become increasingly important. We assessed rates of disagreement between the three U.S. Fish and Wildlife Service–approved acoustic identification software programs (Kaleidoscope Pro 4.2.0, Echoclass 3.1, and Bat Call Identification 2.7d) and manual visual identification using acoustic data collected during summers from 2003 to 2017 at Fort Drum, New York. We assessed the percentage of agreement between programs through pairwise comparisons on a total nightly count level, individual file level (e.g., individual echolocation pass call file), and grouped maximum likelihood estimate level (e.g., probability values that a species is misclassified as present when in fact it is absent) using preplanned contrasts, Akaike Information Criterion, and annual confusion matrices. Interprogram agreement on an individual file level was low, as measured by Cohen's Kappa (0.2–0.6). However, site-night level pairwise comparative analysis indicated that program agreement was higher (40–90%) using single season occupancy metrics. In comparing analytical outcomes of our different datasets (i.e., how comparable programs and visual identification are regarding the relationship between environmental conditions and bat activity), we determined high levels of congruency in the relative rankings of the model as well as the relative level of support for each individual model. This indicated that among individual software packages, when analyzing bat calls, there was consistent ecological inference beyond the file-by-file level at the scales used by managers. Depending on objectives, we believe our results can help users choose automated software and maximum likelihood estimate thresholds more appropriate for their needs and allow for better cross-comparison of studies using different automated acoustic software.

Journal of Fish and Wildlife Management↗

USGS2012

Approaches in highly parameterized inversion - PEST++, a Parameter ESTimation code optimized for large environmental models

An object-oriented parameter estimation code was developed to incorporate benefits of object-oriented programming techniques for solving large parameter estimation modeling problems. The code is written in C++ and is a formulation and expansion of the algorithms included in PEST, a widely used parameter estimation code written in Fortran. The new code is called PEST++ and is designed to lower the barriers of entry for users and developers while providing efficient algorithms that can accommodate large, highly parameterized problems. This effort has focused on (1) implementing the most popular features of PEST in a fashion that is easy for novice or experienced modelers to use and (2) creating a software design that is easy to extend; that is, this effort provides a documented object-oriented framework designed from the ground up to be modular and extensible. In addition, all PEST++ source code and its associated libraries, as well as the general run manager source code, have been integrated in the Microsoft Visual Studio® 2010 integrated development environment. The PEST++ code is designed to provide a foundation for an open-source development environment capable of producing robust and efficient parameter estimation tools for the environmental modeling community into the future.

Techniques and Methods↗

USGS2009

Mercury Loads in the South River and Simulation of Mercury Total Maximum Daily Loads (TMDLs) for the South River, South Fork Shenandoah River, and Shenandoah River: Shenandoah Valley, Virginia

Due to elevated levels of methylmercury in fish, three streams in the Shenandoah Valley of Virginia have been placed on the State's 303d list of contaminated waters. These streams, the South River, the South Fork Shenandoah River, and parts of the Shenandoah River, are downstream from the city of Waynesboro, where mercury waste was discharged from 1929-1950 at an industrial site. To evaluate mercury contamination in fish, this total maximum daily load (TMDL) study was performed in a cooperative effort between the U.S. Geological Survey, the Virginia Department of Environmental Quality, and the U.S. Environmental Protection Agency. The investigation focused on the South River watershed, a headwater of the South Fork Shenandoah River, and extrapolated findings to the other affected downstream rivers. A numerical model of the watershed, based on Hydrological Simulation Program-FORTRAN (HSPF) software, was developed to simulate flows of water, sediment, and total mercury. Results from the investigation and numerical model indicate that contaminated flood-plain soils along the riverbank are the largest source of mercury to the river. Mercury associated with sediment accounts for 96 percent of the annual downstream mercury load (181 of 189 kilograms per year) at the mouth of the South River. Atmospherically deposited mercury contributes a smaller load (less than 1 percent) as do point sources, including current discharge from the historic industrial source area. In order to determine how reductions of mercury loading to the stream could reduce methylmercury concentrations in fish tissue below the U.S. Environmental Protection Agency criterion of 0.3 milligrams per kilogram, multiple scenarios were simulated. Bioaccumulation of mercury was expressed with a site-specific exponential relation between aqueous total mercury and methylmercury in smallmouth bass, the indicator fish species. Simulations indicate that if mercury loading were to decrease by 98.9 percent from 189 to 2 kilograms per year, fish tissue methylmercury concentrations would drop below 0.3 milligrams per kilogram. Based on the simulations, the estimated maximum load of total mercury that can enter the South River without causing fish tissue methylmercury concentrations to rise above 0.3 milligrams per kilogram is 2.03 kilograms per year for the South River, and 4.12 and 6.06 kilograms per year for the South Fork Shenandoah River and Shenandoah River, respectively.

Scientific Investigations Report↗

USGS2013

NCWin — A Component Object Model (COM) for processing and visualizing NetCDF data

NetCDF (Network Common Data Form) is a data sharing protocol and library that is commonly used in large-scale atmospheric and environmental data archiving and modeling. The NetCDF tool described here, named NCWin and coded with Borland C + + Builder, was built as a standard executable as well as a COM (component object model) for the Microsoft Windows environment. COM is a powerful technology that enhances the reuse of applications (as components). Environmental model developers from different modeling environments, such as Python, JAVA, VISUAL FORTRAN, VISUAL BASIC, VISUAL C + +, and DELPHI, can reuse NCWin in their models to read, write and visualize NetCDF data. Some Windows applications, such as ArcGIS and Microsoft PowerPoint, can also call NCWin within the application. NCWin has three major components: 1) The data conversion part is designed to convert binary raw data to and from NetCDF data. It can process six data types (unsigned char, signed char, short, int, float, double) and three spatial data formats (BIP, BIL, BSQ); 2) The visualization part is designed for displaying grid map series (playing forward or backward) with simple map legend, and displaying temporal trend curves for data on individual map pixels; and 3) The modeling interface is designed for environmental model development by which a set of integrated NetCDF functions is provided for processing NetCDF data. To demonstrate that the NCWin can easily extend the functions of some current GIS software and the Office applications, examples of calling NCWin within ArcGIS and MS PowerPoint for showing NetCDF map animations are given.

Cartographic Journal↗

USGS2010

Bayesian Inference: with ecological applications

This text provides a mathematically rigorous yet accessible and engaging introduction to Bayesian inference with relevant examples that will be of interest to biologists working in the fields of ecology, wildlife management and environmental studies as well as students in advanced undergraduate statistics.. This text opens the door to Bayesian inference, taking advantage of modern computational efficiencies and easily accessible software to evaluate complex hierarchical models.

Book↗

USGS2020

msocc: Fit and analyse computationally efficient multi‐scale occupancy models in R

Environmental DNA (eDNA) sampling is a promising tool for the detection of rare and cryptic taxa, such as aquatic pathogens, parasites and invasive species. Environmental DNA sampling workflows commonly rely on multi‐stage hierarchical sampling designs that induce complicated dependencies within the data. This complex dependence structure can be intuitively modelled with Bayesian multi‐scale occupancy models. However, current software for such models are computationally demanding, impeding their use. We present an r package, msocc , that implements a data augmentation strategy to fit fully Bayesian, computationally efficient multi‐scale occupancy models. The msocc package allows users to fit multi‐scale occupancy models, to estimate and visualize posterior summaries of site, sample and replicate‐level occupancy, and to compare different models using Bayesian information criterion. Additionally, we provide a supplemental web application that allows users to investigate study design for multi‐scale occupancy models and acts as a graphical user interface to the msocc package. The utility of the msocc package is illustrated on a published dataset and the functions in msocc are compared to the primary Bayesian toolkit for multi‐scale occupancy modelling, eDNAoccupancy , using various computational benchmarks. These benchmarks indicate that msocc is capable of fitting models 50 times faster than eDNAoccupancy . We hope that access to software that efficiently fits, analyses and conducts study design investigations for multi‐scale occupancy models facilitates their implementation by the research and wildlife management communities.

Methods in Ecology and Evolution↗

USGS2022

Coupling validation effort with in situ bioacoustic data improves estimating relative activity and occupancy for multiple species with cross-species misclassifications

The increasing complexity and pace of ecological change requires natural resource managers to consider entire species assemblages. Acoustic recording units (ARUs) require minimal cost and effort to deploy and inform relative activity, or encounter rates, for multiple species simultaneously. ARU-based surveys require post-processing of the recordings via software algorithms that assign a species label to each recording. The automated classification process can result in cross-species misidentifications that should be accounted for when employing statistical modelling for conservation decision-making. Using simulation and ARU-based detection counts from 17 bat species in British Columbia, Canada, we investigate three strategies for adjusting statistical inference for species misclassification: (a) ‘coupling’ ambiguous and unambiguous detections by validating a subset of survey events post-hoc, (b) using a calibration dataset on the software algorithm's (in)accuracy for species identification or (c) specifying informative Bayesian priors on classification probabilities. We explore the impact of different Bayesian prior specifications for the classification probabilities on posterior estimation. We then consider how the quantity of data validated post-hoc impacts model convergence and resulting inferences for bat species relative activity as related to nightly conditions and yearly site occupancy after accounting for site-level environmental variables. Coupled methods resulted in less bias and uncertainty when estimating relative activity and species classification probabilities relative to calibration approaches. We found that species that were difficult-to-detect and those that were often inaccurately identified by the software required more validation effort than more easily detected and/or identified species. Our results suggest that, when possible, acoustic surveys should rely on coupled validated detection information to account for false-positive detections, rather than uncoupled calibration datasets. However, if the assemblage of interest contains a large number of rarely detected or less prevalent species, an intractable amount of effort may be required, suggesting there are benefits to curating a calibration dataset that is representative of the observation process. Our findings provide insights into the practical challenges associated with statistical analyses of ARU data and possible analytical solutions to support reliable and cost-effective decision-making for wildlife conservation/management in the face of known sources of observation errors.

British Columbia↗

USGS2025

Assessing nonpoint-source uranium pollution in an irrigated stream-aquifer system

Uranium (U) in rocks and soils of arid and semi-arid environments can be mobilized by irrigation and fertilization, posing environmental and health risks. Elevated U, along with selenium (Se) and nitrate (NO 3 ) co-constituents, necessitates careful monitoring and management. We developed a distributed-parameter numerical model to assess U pollution in an irrigated stream-aquifer system, applying it to a 552 km 2 region in Colorado's Lower Arkansas River Valley (LARV) over 14 years. A MODFLOW model, describing groundwater and stream flow, was coupled with an RT3D-OTIS model to portray reactive U transport. Calibration using the PESTPP-iES iterative ensemble smoother (iES) software indicated good agreement with observed U concentrations. The model revealed substantial and variable U levels across the LARV, highlighting potential hotspots and possible contributing factors, such as geological composition of the bedrock and near-surface shale and aquifer sediments derived from them, irrigation practices, and riparian landscape. U levels exceed the chronic standard (85th percentile = 30 μg/L, set by the US Environmental Protection Agency), which is the permissible regulatory threshold, in groundwater across 44 % of the region and along the river by an average factor of 2.9. Simulated average U concentrations in the non-riparian aquifer and river are 124 μg/L and 60 μg/L, respectively, compared with 112 μg/L and 62 μg/L for measured values. The average 85th percentile U concentration is 222 μg/L in the aquifer and 82 μg/L in the river. Average simulated U mass loading to the river is 0.17 kg/day per km, compared to an estimated 0.23 kg/day per km. Findings provide a baseline for comparing future simulated outcomes of alternative best management practices (BMPs) for U pollution mitigation and offer a methodology applicable to other irrigated regions.

Colorado↗