Search USGS⌕ Search

SEARCH · Search USGS

Results for “Environmental Modelling and Software”

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 181 records · Page 10Linked to original sources
USGS2018

VS2DRTI: Simulating heat and reactive solute transport in variably saturated porous media

Variably saturated groundwater flow, heat transport, and solute transport are important processes in environmental phenomena, such as the natural evolution of water chemistry of aquifers and streams, the storage of radioactive waste in a geologic repository, the contamination of water resources from acid-rock drainage, and the geologic sequestration of carbon dioxide. Up to now, our ability to simulate these processes simultaneously with fully coupled reactive transport models has been limited to complex and often difficult-to-use models. To address the need for a simple and easy-to-use model, the VS2DRTI software package has been developed for simulating water flow, heat transport, and reactive solute transport through variably saturated porous media. The underlying numerical model, VS2DRT, was created by coupling the flow and transport capabilities of the VS2DT and VS2DH models with the equilibrium and kinetic reaction capabilities of PhreeqcRM. Flow capabilities include two-dimensional, constant-density, variably saturated flow; transport capabilities include both heat and multicomponent solute transport; and the reaction capabilities are a complete implementation of geochemical reactions of PHREEQC. The graphical user interface includes a preprocessor for building simulations and a postprocessor for visual display of simulation results. To demonstrate the simulation of multiple processes, the model is applied to a hypothetical example of injection of heated waste water to an aquifer with temperature-dependent cation exchange. VS2DRTI is freely available public domain software.

Groundwater↗
USGS2016

Groundwater contaminant plume maps and volumes, 100-K and 100-N Areas, Hanford Site, Washington

This study provides an independent estimate of the areal and volumetric extent of groundwater contaminant plumes which are affected by waste disposal in the 100-K and 100-N Areas (study area) along the Columbia River Corridor of the Hanford Site. The Hanford Natural Resource Trustee Council requested that the U.S. Geological Survey perform this interpolation to assess the accuracy of delineations previously conducted by the U.S. Department of Energy and its contractors, in order to assure that the Natural Resource Damage Assessment could rely on these analyses. This study is based on previously existing chemical (or radionuclide) sampling and analysis data downloaded from publicly available Hanford Site Internet sources, geostatistically selected and interpreted as representative of current (from 2009 through part of 2012) but average conditions for groundwater contamination in the study area. The study is limited in scope to five contaminants—hexavalent chromium, tritium, nitrate, strontium-90, and carbon-14, all detected at concentrations greater than regulatory limits in the past. All recent analytical concentrations (or activities) for each contaminant, adjusted for radioactive decay, non-detections, and co-located wells, were converted to log-normal distributions and these transformed values were averaged for each well location. The log-normally linearized well averages were spatially interpolated on a 50 × 50-meter (m) grid extending across the combined 100-N and 100-K Areas study area but limited to avoid unrepresentative extrapolation, using the minimum curvature geostatistical interpolation method provided by SURFER ® data analysis software. Plume extents were interpreted by interpolating the log-normally transformed data, again using SURFER ® , along lines of equal contaminant concentration at an appropriate established regulatory concentration . Total areas for each plume were calculated as an indicator of relative environmental damage. These plume extents are shown graphically and in tabular form for comparison to previous estimates. Plume data also were interpolated to a finer grid (10 × 10 m) for some processing, particularly to estimate volumes of contaminated groundwater. However, hydrogeologic transport modeling was not considered for the interpolation. The compilation of plume extents for each contaminant also allowed estimates of overlap of the plumes or areas with more than one contaminant above regulatory standards. A mapping of saturated aquifer thickness also was derived across the 100-K and 100–N study area, based on the vertical difference between the groundwater level (water table) at the top and the altitude of the top of the Ringold Upper Mud geologic unit, considered the bottom of the uppermost unconfined aquifer. Saturated thickness was calculated for each cell in the finer (10 × 10 m) grid. The summation of the cells’ saturated thickness values within each polygon of plume regulatory exceedance provided an estimate of the total volume of contaminated aquifer, and the results also were checked using a SURFER ® volumetric integration procedure. The total volume of contaminated groundwater in each plume was derived by multiplying the aquifer saturated thickness volume by a locally representative value of porosity (0.3). Estimates of the uncertainty of the plume delineation also are presented. “Upper limit” plume delineations were calculated for each contaminant using the same procedure as the “average” plume extent except with values at each well that are set at a 95-percent upper confidence limit around the log-normally transformed mean concentrations, based on the standard error for the distribution of the mean value in that well; “lower limit” plumes are calculated at a 5-percent confidence limit around the geometric mean. These upper- and lower-limit estimates are considered unrealistic because the statistics were increased or decreased at each well simultaneously and were not adjusted for correlation among the well distributions (i.e., it is not realistic that all wells would be high simultaneously). Sources of the variability in the distributions used in the upper- and lower-extent maps include time varying concentrations and analytical errors. The plume delineations developed in this study are similar to the previous plume descriptions developed by U.S. Department of Energy and its contractors. The differences are primarily due to data selection and interpolation methodology. The differences in delineated plumes are not sufficient to result in the Hanford Natural Resource Trustee Council adjusting its understandings of contaminant impact or remediation.

Washington↗
USGS1995

Near real-time monitoring of seismic events and status of portable digital recorders using satellite telemetry

Near real-time monitoring of seismic events and status of portable 16-bit digital recorders has been established for arrays near Parkfield, Mammoth Lakes, and San Francisco, California. This monitoring system provides near real-time seismic event identification (rough location and magnitude) and a cost-effective means to maintain arrays at near 100% operational level. Principal objectives in the design of this system have been portability and low-cost telemetry. The system has been developed to use portable digital seismic recorders ( GEOS —General Earthquake Observation System) and portable data collection platforms ( DCP 's) for the Geostationary Operational Environmental Satellite (GEOS) telemetry system. Data are transferred asynchronously from the GEOS seismic system through a microprocessor-controlled interface every 10 min. The interface stores, determines priority, converts, and synchronously transfers these data to a Sutron Corp. model 8004 DCP for transmission through the GEOS satellite telemetry system. Event parameters include trigger time, peak amplitude, time of peak amplitude, and event duration. Instrument configuration parameters, transmitted at system start-up time and every 24 hr, include recording parameters, trigger parameters, GEOS software version, clock reference, and location parameter. Instrument status includes battery voltage, number of events, and percentage of tape usage. These data are transmitted as appropriate to the U.S. Geological Survey satellite downlink and computers located in Menlo Park, California, where they are processed and displayed.

Bulletin of the Seismological Society of America↗
USGS2008

User's Manual for the Object User Interface (OUI): An Environmental Resource Modeling Framework

The Object User Interface is a computer application that provides a framework for coupling environmental-resource models and for managing associated temporal and spatial data. The Object User Interface is designed to be easily extensible to incorporate models and data interfaces defined by the user. Additionally, the Object User Interface is highly configurable through the use of a user-modifiable, text-based control file that is written in the eXtensible Markup Language. The Object User Interface user's manual provides (1) installation instructions, (2) an overview of the graphical user interface, (3) a description of the software tools, (4) a project example, and (5) specifications for user configuration and extension.

Open-File Report↗
USGS2024

Simulated mean monthly groundwater-transported nitrogen loads in watersheds on the north shore of Long Island Sound, 1993–2022

Elevated nitrogen loads are pervasive in the Long Island Sound, an estuary that receives freshwater and nutrients from both surface-water and groundwater discharge. Surface-water nitrogen loads to the Long Island Sound are relatively well characterized, but less is known about groundwater-transported nitrogen loads. Prior work on the northern shore of Long Island Sound (Connecticut and areas of New York and Rhode Island) suggested that groundwater travel times are relatively short (median less than 2 years) and that decade-long nutrient legacies are not widespread. Because the travel times are short, groundwater flow and nutrient loads likely vary substantially between months. In the current study, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency’s Long Island Sound Study and the Connecticut Department of Energy and Environmental Protection, developed a set of models to better characterize spatial and temporal patterns of groundwater-transported nitrogen loading from atmospheric deposition, septic systems, and fertilizers within the study area. The models provide an estimate, with uncertainty, of groundwater-transported nitrogen loads in the study area, filling a key gap in the nitrogen budget for Long Island Sound. The models also highlight the spatial and temporal variation in nitrogen loading throughout the study area. The modeling workflow involved four models. (1) A soil-water-balance model was developed by using the Soil-Water-Balance software to simulate groundwater recharge across the study area for water years 2005 through 2022. The simulated mean monthly recharge from the soil-water-balance model was used as input into a groundwater-flow model. (2) The groundwater-flow model was developed by using the MODFLOW 6 software and data for water years 1993 through 2022 and simulates average monthly hydrologic conditions. The groundwater-flow model was calibrated by using the Iterative Ensemble Smoother method within the PEST++ software. The Iterative Ensemble Smoother method generates an ensemble of sets of parameter values, with each set producing reasonable simulated hydrologic parameter values. (3) An ensemble of MODPATH particle-tracking simulations were run to generate particle flow paths and travel times, with each simulation using a different set of the flow model parameters. (4) A nitrogen load model uses the MODPATH simulation outputs to track nitrogen from the land surface through multiple attenuation zones until it discharges into fresh or saline surface water. As with the groundwater-flow model, the nitrogen model simulated average monthly groundwater-transported nitrogen loads for water years 1993 through 2022. One novel aspect of the nitrogen load model is that the nitrogen attenuation parameters were calibrated to observed nitrogen loads. Across the ensemble of simulated nitrogen loads, the median study-area-wide monthly simulated nitrogen loads from the aquifer to Long Island Sound throughout the year ranged from 900 to 18,600 kilograms of nitrogen per day, with a median load of 5,100 kilograms of nitrogen per day. The simulated loads were based on average monthly conditions for water years 1993 through 2022. Loads were highest during the winter and early spring and lowest during the late summer. However, simulated travel times for groundwater and nitrogen loads discharged to Long Island Sound during summer were longer than travel times for groundwater and loads discharged during the winter, indicating that, on average, groundwater discharged during summer traveled along different, and longer, flow paths, than groundwater discharged during winter. This indicates that summer loads would respond more slowly to changes in nitrogen inputs at the water table than winter loads. Over the entire study area, approximately 15 percent of the simulated load is from atmospheric deposition sources, 30 to 40 percent is from fertilizer, and 50 to 60 percent is from septic systems. The final analysis of the study involved simulating the change in groundwater-transported nitrogen load in response to upgrading septic systems or reducing fertilizing inputs to areas of turf grass. Both management interventions reduced the groundwater-transported nitrogen load, and reductions were greater in areas with greater loads from septic systems or turf-grass fertilizers. The delay between management actions and substantial reductions in groundwater-transported nitrogen loads varied seasonally; loads during the late summer months remained elevated longer than the winter loads.

Connecticut, Rhode Island↗
USGS2012

Approaches in highly parameterized inversion-PESTCommander, a graphical user interface for file and run management across networks

Models of environmental systems have become increasingly complex, incorporating increasingly large numbers of parameters in an effort to represent physical processes on a scale approaching that at which they occur in nature. Consequently, the inverse problem of parameter estimation (specifically, model calibration) and subsequent uncertainty analysis have become increasingly computation-intensive endeavors. Fortunately, advances in computing have made computational power equivalent to that of dozens to hundreds of desktop computers accessible through a variety of alternate means: modelers have various possibilities, ranging from traditional Local Area Networks (LANs) to cloud computing. Commonly used parameter estimation software is well suited to take advantage of the availability of such increased computing power. Unfortunately, logistical issues become increasingly important as an increasing number and variety of computers are brought to bear on the inverse problem. To facilitate efficient access to disparate computer resources, the PESTCommander program documented herein has been developed to provide a Graphical User Interface (GUI) that facilitates the management of model files ("file management") and remote launching and termination of "slave" computers across a distributed network of computers ("run management"). In version 1.0 described here, PESTCommander can access and ascertain resources across traditional Windows LANs: however, the architecture of PESTCommander has been developed with the intent that future releases will be able to access computing resources (1) via trusted domains established in Wide Area Networks (WANs) in multiple remote locations and (2) via heterogeneous networks of Windows- and Unix-based operating systems. The design of PESTCommander also makes it suitable for extension to other computational resources, such as those that are available via cloud computing. Version 1.0 of PESTCommander was developed primarily to work with the parameter estimation software PEST; the discussion presented in this report focuses on the use of the PESTCommander together with Parallel PEST. However, PESTCommander can be used with a wide variety of programs and models that require management, distribution, and cleanup of files before or after model execution. In addition to its use with the Parallel PEST program suite, discussion is also included in this report regarding the use of PESTCommander with the Global Run Manager GENIE, which was developed simultaneously with PESTCommander.

Techniques and Methods↗
USGS2016

Using spatial capture–recapture to elucidate population processes and space-use in herpetological studies

The cryptic behavior and ecology of herpetofauna make estimating the impacts of environmental change on demography difficult; yet, the ability to measure demographic relationships is essential for elucidating mechanisms leading to the population declines reported for herpetofauna worldwide. Recently developed spatial capture–recapture (SCR) methods are well suited to standard herpetofauna monitoring approaches. Individually identifying animals and their locations allows accurate estimates of population densities and survival. Spatial capture–recapture methods also allow estimation of parameters describing space-use and movement, which generally are expensive or difficult to obtain using other methods. In this paper, we discuss the basic components of SCR models, the available software for conducting analyses, and the experimental designs based on common herpetological survey methods. We then apply SCR models to Red-backed Salamander ( Plethodon cinereus ), to determine differences in density, survival, dispersal, and space-use between adult male and female salamanders. By highlighting the capabilities of SCR, and its advantages compared to traditional methods, we hope to give herpetologists the resource they need to apply SCR in their own systems.

Journal of Herpetology↗
USGS2014

Hydraulic modeling development and application in water resources engineering

The use of modeling has become widespread in water resources engineering and science to study rivers, lakes, estuaries, and coastal regions. For example, computer models are commonly used to forecast anthropogenic effects on the environment, and to help provide advanced mitigation measures against catastrophic events such as natural and dam-break floods. Linking hydraulic models to vegetation and habitat models has expanded their use in multidisciplinary applications to the riparian corridor. Implementation of these models in software packages on personal desktop computers has made them accessible to the general engineering community, and their use has been popularized by the need of minimal training due to intuitive graphical user interface front ends. Models are, however, complex and nontrivial, to the extent that even common terminology is sometimes ambiguous and often applied incorrectly. In fact, many efforts are currently under way in order to standardize terminology and offer guidelines for good practice, but none has yet reached unanimous acceptance. This chapter provides a view of the elements involved in modeling surface flows for the application in environmental water resources engineering. It presents the concepts and steps necessary for rational model development and use by starting with the exploration of the ideas involved in defining a model. Tangible form of those ideas is provided by the development of a mathematical and corresponding numerical hydraulic model, which is given with a substantial amount of detail. The issues of model deployment in a practical and productive work environment are also addressed. The chapter ends by presenting a few model applications highlighting the need for good quality control in model validation.

Book chapter↗
USGS2014

An application of Social Values for Ecosystem Services (SolVES) to three national forests in Colorado and Wyoming

Despite widespread recognition that social-value information is needed to inform stakeholders and decision makers regarding trade-offs in environmental management, it too often remains absent from ecosystem service assessments. Although quantitative indicators of social values need to be explicitly accounted for in the decision-making process, they need not be monetary. Ongoing efforts to map such values demonstrate how they can also be made spatially explicit and relatable to underlying ecological information. We originally developed Social Values for Ecosystem Services (SolVES) as a tool to assess, map, and quantify nonmarket values perceived by various groups of ecosystem stakeholders. With SolVES 2.0 we have extended the functionality by integrating SolVES with Maxent maximum entropy modeling software to generate more complete social-value maps from available value and preference survey data and to produce more robust models describing the relationship between social values and ecosystems. The current study has two objectives: (1) evaluate how effectively the value index, a quantitative, nonmonetary social-value indicator calculated by SolVES, reproduces results from more common statistical methods of social-survey data analysis and (2) examine how the spatial results produced by SolVES provide additional information that could be used by managers and stakeholders to better understand more complex relationships among stakeholder values, attitudes, and preferences. To achieve these objectives, we applied SolVES to value and preference survey data collected for three national forests, the Pike and San Isabel in Colorado and the Bridger–Teton and the Shoshone in Wyoming. Value index results were generally consistent with results found through more common statistical analyses of the survey data such as frequency, discriminant function, and correlation analyses. In addition, spatial analysis of the social-value maps produced by SolVES provided information that was useful for explaining relationships between stakeholder values and forest uses. Our results suggest that SolVES can effectively reproduce information derived from traditional statistical analyses while adding spatially explicit, social-value information that can contribute to integrated resource assessment, planning, and management of forests and other ecosystems.

Colorado;Wyoming↗
USGS1991

System requirements specification for the U.S. Geological Survey's National Water Information System II

The U.S. Geological Survey's Water Resources Division is designing and developing a new computer software system for processing and storing hydrologic data. This system, the National Water Information System (NWIS-II) will replace the current water-data and information systems: the National Water Data Storage and Retrieval System (WATSTORE), the National Water Data Exchange (NAWDEX), and the National Water-Use Information System, as well as the current National Water Information System (NWIS-I). The new system will be utilized across the Nation on the Distributed Information System (DIS-II), a network of 32-bit microcomputers on a wide-area network of desktop workstations and related hardware. The purpose of the NWIS-II System Requirements Specification is to integrate requirements described by eight User Groups representing the disciplines for Surface Water, Ground Water, Water Quality, Water Sediment, Water Use, Biology, Geographic Information, and National Water Data Exchange. The requirements specified in this document will serve as the framework for subsequent design and development activities throughout the life cycle of NWIS-II to ensure that the development products reflect the needs of the users. Specification details and associated data-base plans are contained in the description of the integrated functional requirements, the logical data model and the data dictionary, the plan for transferring data from the existing system to NWIS-II, and the descriptions of performance requirements and design constraints. The NWIS-II will be consistent across hydrologic disciplines and the other major Division software systems. The color-graphics desktop environment of NWIS-II will utilize multi-window and multi-tasking capabilities with multilevel of help for the users. Data will be distributed among nodes and be automatically accessible during input, editing, verification and analysis of data. Hydrologic data will be cataloged independent of the Environmental Protection Agency's storage and retrieval (STORET) codes. Character, digital, and graphical data will be entered into NWIS-II from digital and analog recorders, external files, or keyboard through standard and user-defined input forms. Data will be verified manually or automatically during input, retrieval, and editing of data. Numerous reference lists will be maintained so that entry of codes on forms will not be required. User-defined series of computer actions for standard programs will be supported for processing basic data and other frequently repeated processes. Data will be verified and exchanged through links to other DIS-II software, including a Geographic Information System, and selected external software. Quality assurance of data will include aging of data through change on status type and tracking access violations. Output of data will be to a variety of media in standard and user-defined formats. The data base will be preserved by operations such as backup, recovery, audit trails, archiving, and data histories. In addition, an index of water data from USGS and other agencies will be maintained.

Open-File Report↗
USGS2004

Phosphorus and suspended sediment load estimates for the Lower Boise River, Idaho, 1994-2002

The U.S. Geological Survey used LOADEST, newly developed load estimation software, to develop regression equations and estimate loads of total phosphorus (TP), dissolved orthophosphorus (OP), and suspended sediment (SS) from January 1994 through September 2002 at four sites on the lower Boise River: Boise River below Diversion Dam near Boise, Boise River at Glenwood Bridge at Boise, Boise River near Middleton, and Boise River near Parma. The objective was to help the Idaho Department of Environmental Quality develop and implement total maximum daily loads (TMDLs) by providing spatial and temporal resolution for phosphorus and sediment loads and enabling load estimates made by mass balance calculations to be refined and validated. Regression models for TP and OP generally were well fit on the basis of regression coefficients of determination (R2), but results varied in quality from site to site. The TP and OP results for Glenwood probably were affected by the upstream wastewater-treatment plant outlet, which provides a variable phosphorus input that is unrelated to river discharge. Regression models for SS generally were statistically well fit. Regression models for Middleton for all constituents, although statistically acceptable, were of limited usefulness because sparse and intermittent discharge data at that site caused many gaps in the resulting estimates. Although the models successfully simulated measured loads under predominant flow conditions, errors in TP and SS estimates at Middleton and in TP estimates at Parma were larger during high- and low-flow conditions. This shortcoming might be improved if additional concentration data for a wider range of flow conditions were available for calibrating the model. The average estimated daily TP load ranged from less than 250 pounds per day (lb/d) at Diversion to nearly 2,200 lb/d at Parma. Estimated TP loads at all four sites displayed cyclical variations coinciding with seasonal fluctuations in discharge. Estimated annual loads of TP ranged from less than 8 tons at Diversion to 570 tons at Parma. Annual loads of dissolved OP peaked in 1997 at all sites and were consistently higher at Parma than at the other sites. The ratio of OP to TP varied considerably throughout the year at all sites. Peaks in the OP:TP ratio occurred primarily when flows were at their lowest annual stages; estimated seasonal OP:TP ratios were highest in autumn at all sites. Conversely, when flows were high, the ratio was low, reflecting increased TP associated with particulate matter during high flows. Parma exhibited the highest OP:TP ratio during all seasons, at least 0.60 in spring and nearly 0.90 in autumn. Similar OP:TP ratios were estimated at Glenwood. Whereas the OP:TP ratio for Parma and Glenwood peaked in November or December, decreased from January through May, and increased again after June, estimates for Diversion showed nearly the opposite pattern ? ratios were highest in July and lowest in January and February. This difference might reflect complex biological and geochemical processes involving nutrient cycling in Lucky Peak Lake, but further data are needed to substantiate this hypothesis. Estimated monthly average SS loads were highest at Diversion, about 400 tons per day (ton/d). Average annual loads from 1994 through 2002 were 144,000 tons at Diversion, 33,000 tons at Glenwood, and 88,000 tons at Parma. Estimated SS loads peaked in the spring at all sites, coinciding with high flows. Increases in TP in the reach from Diversion to Glenwood ranged from 200 to 350 lb/d. Decreases in TP were small in this reach only during high flows in January and February 1997. Decreases in SS, were large during high-flow conditions indicating sediment deposition in the reach. Intermittent data at Middleton indicated that increases and decreases in TP in the reach from Glenwood to Middleton were during low- and high-flow conditions, respectively. All constituents increased in the r

Scientific Investigations Report↗
USGS2021

Assessment of water quality and discharge in the Herring River, Wellfleet, Massachusetts, November 2015 to September 2017

The U.S. Geological Survey, Cape Cod National Seashore of the National Park Service, and Friends of Herring River cooperated from 2015 to 2017 to assess nutrient concentrations and fluxes across the ocean-estuary boundary at a dike on the Herring River in Wellfleet, Massachusetts. The purpose of this assessment was to characterize environmental conditions prior to a future removal of the dike, which has restricted saltwater inputs into the Herring River watershed for more than 100 years. Water temperature, dissolved oxygen, pH, and specific conductance were monitored continuously, and flow-weighted composite samples were collected approximately twice per month at the ocean-estuary boundary. Bidirectional discharge was computed for the U.S. Geological Survey Herring River at Chequessett Neck Road at Wellfleet, Massachusetts, streamgage (011058798) by using a stage-area rating and index-velocity ratings developed with acoustic Doppler current profile measurements made upstream and downstream from the dike. LOADEST regression modeling software was used to estimate nutrient fluxes (loads) from composite, paired nutrient concentration and discharge data in conjunction with continuous discharge data. Temperature, dissolved oxygen, pH, and specific conductance were also monitored continuously on two tributaries to the Herring River, Pole Dike Creek and Bound Brook, from late-May 2016 to mid-June 2017. Composite or discrete water samples were collected from the tributaries approximately twice per month in most months from late-May 2016 to mid-June 2017 and analyzed for total nitrogen, total phosphorus, and dissolved organic carbon. Flow-weighted concentrations of ammonium, nitrate, and total nitrogen on the Herring River at the dike on the ebb tide generally varied between 0.01 and 0.1, 0.003 and 0.03, and 0.3 and 0.7 milligram per liter as nitrogen, respectively. Flow-weighted concentrations of orthophosphate, total dissolved phosphorus, and total phosphorus generally varied between 0.002 and 0.02, 0.003 and 0.06, and 0.03 and 0.1 milligram per liter as phosphorus, respectively, on the ebb tide. Flow-weighted concentrations of silicate and dissolved organic carbon on the ebb tide generally varied between 0.08 and 3.0 milligrams per liter of silica (silicon dioxide), and 1.7 and 5.6 milligrams per liter of carbon, respectively. Ebb tide concentrations of nitrate were highest in winter and lowest in summer. By contrast, ebb tide concentrations of phosphorus species were highest in late summer and early fall and lowest in winter. Silica and dissolved organic carbon did not exhibit systematic variation in seasonal concentrations. There was uncertainty in estimates of nutrient fluxes, but the LOADEST-estimated fluxes indicated that annual (and in almost all cases seasonal) exports (ebb tides) exceeded inputs (flood tides). Ebb tide concentrations of ammonium, nitrate, total nitrogen, and silica were positively correlated with antecedent cumulative 7-day precipitation.

Massachusetts↗
USGS2014

Basin-scale simulation of current and potential climate changed hydrologic conditions in the Lake Michigan Basin, United States

The Great Lakes Restoration Initiative (GLRI) is the largest public investment in the Great Lakes in two decades. A task force of 11 Federal agencies developed an action plan to implement the initiative. The U.S. Department of the Interior was one of the 11 agencies that entered into an interagency agreement with the U.S. Environmental Protection Agency as part of the GLRI to complete scientific projects throughout the Great Lakes basin. The U.S. Geological Survey, a bureau within the Department of the Interior, is involved in the GLRI to provide scientific support to management decisions as well as measure progress of the Great Lakes basin restoration efforts. This report presents basin-scale simulated current and forecast climatic and hydrologic conditions in the Lake Michigan Basin. The forecasts were obtained by constructing and calibrating a Precipitation-Runoff Modeling System (PRMS) model of the Lake Michigan Basin; the PRMS model was calibrated using the parameter estimation and uncertainty analysis (PEST) software suite. The calibrated model was used to evaluate potential responses to climate change by using four simulated carbon emission scenarios from eight general circulation models released by the World Climate Research Programme’s Coupled Model Intercomparison Project phase 3. Statistically downscaled datasets of these scenarios were used to project hydrologic response for the Lake Michigan Basin. In general, most of the observation sites in the Lake Michigan Basin indicated slight increases in annual streamflow in response to future climate change scenarios. Monthly streamflows indicated a general shift from the current (2014) winter-storage/snowmelt-pulse system to a system with a more equally distributed hydrograph throughout the year. Simulated soil moisture within the basin illustrates that conditions within the basin are also expected to change on a monthly timescale. One effect of increasing air temperature as a result of the changing climate was the appreciable increase in the length of the growing season in the Lake Michigan Basin. The increase in growing season will cause an increase in evapotranspiration across the Lake Michigan Basin, which will directly affect soil moisture and late growing season streamflows. Output from the Lake Michigan Basin PRMS model is available through an online dynamic web mapping service available at (http://pubs.usgs.gov/sir/2014/5175/). The map service includes layers for the each of the 8 global climate models and 4 carbon emission scenarios combinations for 12 hydrologic model state variables. The layers are pre-rendered maps of annual hydrologic response from 1977 through 2099 that provide an easily accessible online method to examine climate change effects across the Lake Michigan Basin.

Indiana, Illinois, Michigan, Wisconsin↗
USGS2024

Evaluation of an impulse-response emulator for groundwater contaminant transport modeling

There is a significant need to develop decision support tools capable of delivering accurate representations of environmental conditions, such as ground and surface water solute concentrations, in a timely and computationally efficient manner. Such tools can be leveraged to assess a large number of potential management strategies for mitigating non-point source pollutants. Here, we assess the effectiveness of the impulse-response emulation approach to approximate process-based groundwater model estimates of solute transport from MODFLOW and MT3D over a wide range of model inputs and parameters, with the goal of assessing where in parameter space the assumptions underlying this emulation approach are valid. The impulse-response emulator was developed using the sensitivity analysis utilities in the PEST++ software suite and is capable of approximating MODFLOW/MT3D estimates of solute transport over a large portion of the parameter space tested, except in cases where the Courant number is above 0.5. Across all runs tested, the highest percent errors were at the plume fronts. These results suggest that the impulse-response approach may be suitable for emulation of solute transport models for a wide range of cases, except when high-resolution outputs are needed, or when very low concentrations at plume edges are of particular interest.

Groundwater↗
USGS2014

Users' guide to system dynamics model describing Coho salmon survival in Olema Creek, Point Reyes National Seashore, Marin County, California

The system dynamics model described in this report is the result of a collaboration between U.S. Geological Survey (USGS) scientists and National Park Service (NPS) San Francisco Bay Area Network (SFAN) staff, whose goal was to develop a methodology to integrate inventory and monitoring data to better understand ecosystem dynamics and trends using salmon in Olema Creek, Marin County, California, as an example case. The SFAN began monitoring multiple life stages of coho salmon (Oncorhynchus kisutch) in Olema Creek during 2003 (Carlisle and others, 2013), building on previous monitoring of spawning fish and redds. They initiated water-quality and habitat monitoring, and had access to flow and weather data from other sources. This system dynamics model of the freshwater portion of the coho salmon life cycle in Olema Creek integrated 8 years of existing monitoring data, literature values, and expert opinion to investigate potential factors limiting survival and production, identify data gaps, and improve monitoring and restoration prescriptions. A system dynamics model is particularly effective when (1) data are insufficient in time series length and/or measured parameters for a statistical or mechanistic model, and (2) the model must be easily accessible by users who are not modelers. These characteristics helped us meet the following overarching goals for this model: Summarize and synthesize NPS monitoring data with data and information from other sources to describe factors and processes affecting freshwater survival of coho salmon in Olema Creek. Provide a model that can be easily manipulated to experiment with alternative values of model parameters and novel scenarios of environmental drivers. Although the model describes the ecological dynamics of Olema Creek, these dynamics are structurally similar to numerous other coastal streams along the California coast that also contain anadromous fish populations. The model developed for Olema can be used, at least as a starting point, for other watersheds. This report describes each of the model elements with sufficient detail to guide the primary target audience, the NPS resource specialist, to run the model, interpret the results, change the input data to explore hypotheses, and ultimately modify and improve the model. Running the model and interpreting the results does not require modeling expertise on the part of the user. Additional companion publications will highlight other aspects of the model, such as its development, the rationale behind the methodological approach, scenario testing, and discussions of its use. System dynamics models consist of three basic elements: stocks , flows , and converters . Stocks are measurable quantities that can change over time, such as animal populations. Flows are any processes or conditions that change the quantity in a stock over time (Ford, 1999), are expressed in the model as a rate of change, and are diagrammed as arrows to or from stocks. Converters are processes or conditions that change the rate of flows. A converter is connected to a flow with an arrow indicating that it alters the rate of change. Anything that influences the rate of change (such as different environmental conditions, other external factors, or feedbacks from other stocks or flows) is modeled as a converter. For example, the number of fish in a population is appropriately modeled as a stock. Mortality is modeled as a flow because it is a rate of change over time used to determine the number of fish in the population. The density-dependent effect on mortality is modeled as a converter because it influences the rate of morality. Together, the flow and converter change the number, or stock, of juvenile coho. The instructions embedded in the stocks, flows, converters, and the sequence in which they are linked are processed by the simulation software with each completed sequence composing a model run. At each modeled time step within the model run, the stock counts will go up, down, or stay the same based on the modeled flows and the influence of converters on those flows. The model includes a user-friendly interface to change model parameters, which allows park staff and others to conduct sensitivity analyses, incorporate future knowledge, and implement scenarios for various future conditions. The model structure incorporates place holders for relationships that we hypothesize are significant but data are currently lacking. Future climate scenarios project stream temperatures higher than any that have ever been recorded at Olema Creek. Exploring climate change impacts on coho survival is a high priority for park staff, therefore the model provides the user with the option to experiment with hypothesized effects and to incorporate effects based on future observations.

California↗
USGS2017

Purpose, processes, partnerships, and products: four Ps to advance participatory socio-environmental modeling

Including stakeholders in environmental model building and analysis is an increasingly popular approach to understanding ecological change. This is because stakeholders often hold valuable knowledge about socio-environmental dynamics and collaborative forms of modeling produce important boundary objects used to collectively reason about environmental problems. Although the number of participatory modeling (PM) case studies and the number of researchers adopting these approaches has grown in recent years, the lack of standardized reporting and limited reproducibility have prevented PM's establishment and advancement as a cohesive field of study. We suggest a four-dimensional framework (4P) that includes reporting on dimensions of (1) the Purpose for selecting a PM approach (the why ); (2) the Process by which the public was involved in model building or evaluation (the how ); (3) the Partnerships formed (the who ); and (4) the Products that resulted from these efforts (the what ). We highlight four case studies that use common PM software-based approaches (fuzzy cognitive mapping, agent-based modeling, system dynamics, and participatory geospatial modeling) to understand human–environment interactions and the consequences of ecological changes, including bushmeat hunting in Tanzania and Cameroon, agricultural production and deforestation in Zambia, and groundwater management in India. We demonstrate how standardizing communication about PM case studies can lead to innovation and new insights about model-based reasoning in support of ecological policy development. We suggest that our 4P framework and reporting approach provides a way for new hypotheses to be identified and tested in the growing field of PM.

Ecological Applications↗
USGS2000

Monitoring beach morphology changes using small-format aerial photography and digital softcopy photogrammetry

Current methods of monitoring beach morphology changes commonly involve the establishment of Global Positioning System profiles that are surveyed on a regular basis. Although this method produces precise measurements of coastal topography, it is costly in time and effort and may result in large data gaps between profiles. Much of our understanding of coastal dynamics is thus limited by profile spacing and survey frequency. Softcopy photogrammetry is increasingly being used as an alternative to assess shoreline change and study beach morphodynamics. This method of producing three-dimensional topographic models and orthophotographs from digitized aerial photography can aid in filling in the gaps between established profiles. A limiting factor to this technology is the cost of obtaining high-resolution aerial photography. We have developed an aerial mapping system designed to collect data in an efficient and cost-effective way. We use a small-format aerial photography system that can be mounted on a variety of small aircraft on short notice. After a flight, the photographs are scanned, and softcopy photogrammetry software is used to create both DTMs (Digital Terrain Models) and orthophotographs. The DTMs can be compared with existing profiles for accuracy, and volumetric changes can be computed. The orthophotos are used to make precise measurements of the position and morphology of shoreline features. This aerial mapping system is advantageous over previous methods of beach morphology change monitoring because it allows for rapid response to storm events and provides a cost-effective method of establishing a continual monitoring program in erosion hazard areas.

Environmental Geosciences↗
USGS2026

3D semantic mapping of surface geological features

Semantic mapping in 3D is fundamental to a wide range of geoscientific studies and applications, including geomorphology, hazard assessment, and environmental monitoring. However, automatically segmenting geological features from large-scale photogrammetric datasets remains a significant challenge. We present a methodology to address this gap. Using overlapping images collected over environments of interest, Structure-from-Motion (SfM) produces georeferenced point clouds and estimates camera poses. Existing large vision models, such as Segment Anything Model, segment objects in the images, generating pixel-segmentation associations. To produce pixel-point associations, we project the points back onto the camera image planes. As objects are independently segmented across multiple images with different perspectives, we develop a segmentation mosaicking algorithm to build probabilistic point-segmentation associations that combines the pixel-segmentation associations and pixel-point associations. Our methodology is validated using both synthetic data generated by Kubric and real-world UAV-SfM data. The implementation is designed to be compatible with existing SfM software, including Agisoft and OpenDroneMap, for photogrammetry mapping in geoscience studies. As a case study, we apply our method to the semantic mapping of precariously balanced rocks (PBRs), which provide upper-bound constraints on historical ground motion shaking intensity. To support object-level identification of PBRs, we additionally integrated Grounding DINO, enabling text-prompted segmentation of features of interest within UAV imagery. This case study demonstrates the effectiveness of our method in generating a 3D semantic map of PBRs, enabling spatial distribution of PBR fragility for earthquake hazard analysis.

Computers & Geosciences↗