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 451 records · Page 25Linked to original sources
USGS2020

Groundwater withdrawals and regional flow paths at and near Willow Grove and Warminster, Pennsylvania—Data compilation and preliminary simulations for conditions in 1999, 2010, 2013, 2016, and 2017

In 2014, groundwater samples from residential and public supply wells in the vicinity of two former U.S. Navy bases at Willow Grove and Warminster, and an active Air National Guard Station at Horsham, Bucks and Montgomery Counties, Pennsylvania, were found to have concentrations of perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS), which are per- and polyfluoroalkyl substances (PFAS), above U.S. Environmental Protection Agency (EPA) provisional health advisory (HA) levels for drinking water. Five supply wells near the bases were shut down because of PFAS contamination. In 2016, after EPA established a Lifetime HA for PFAS in drinking water that is lower than the provisional HA in place in 2014, at least 13 additional supply wells near the bases were shut down because of PFAS contamination. At the request of the U.S. Navy, and in consultation with other Federal and State agencies and local stakeholders, the U.S. Geological Survey used historical and recent data on well withdrawals, recharge rates, aquifer properties, groundwater levels, and stream base flow to evaluate regional groundwater-flow paths from identified areas of PFAS groundwater contamination or potential PFAS sources at the bases. Groundwater withdrawals near the bases from public supply and other large wells decreased substantially from the 1990s to 2017, increasing the proportion of groundwater recharge that discharged to local streams. A preliminary groundwater-flow model, calibrated using 1,009 groundwater levels and 17 stream base flow estimates, simulated regional flow paths from the bases and showed that recharge at the bases discharged to withdrawal wells and local streams, generally within a mile or two of the bases. Supply and remediation wells at the bases captured some of the recharge on base areas of possible PFAS contamination, whereas other base recharge was simulated to flow to nearby public supply wells and streams, depending on water use and aquifer recharge conditions between 1999 and 2017. The locations of many residential wells near the bases that were identified by the Navy and Air National Guard as having elevated PFAS concentrations were generally consistent with the simulated flow paths from possible sources at the bases. However, there are some areas of observed PFAS contamination where no flow paths from base sources were simulated. Additionally, no data were available on PFAS concentrations in groundwater in some areas of simulated flow paths from base sources. Data and models used for this study are provided in this report and in digital data releases to support further investigations and model revisions.

Pennsylvania↗
USGS2016

Potential postwildfire debris-flow hazards—A prewildfire evaluation for the Jemez Mountains, north-central New Mexico

Wildfire can substantially increase the probability of debris flows, a potentially hazardous and destructive form of mass wasting, in landscapes that have otherwise been stable throughout recent history. Although the exact location, extent, and severity of wildfire or subsequent rainfall intensity and duration cannot be known, probabilities of fire and debris‑flow occurrence for given locations can be estimated with geospatial analysis and modeling. The purpose of this report is to provide information on which watersheds might constitute the most serious potential debris - flow hazards in the event of a large - scale wildfire and subsequent rainfall in the Jemez Mountains. Potential probabilities and estimated volumes of postwildfire debris flows in both the unburned and previously burned areas of the Jemez Mountains and surrounding areas were estimated using empirical debris - flow models developed by the U.S. Geological Survey in combination with fire behavior and burn probability models developed by the U.S. Forest Service. Of the 4,998 subbasins modeled for this study, computed debris-flow probabilities in 671 subbasins were greater than 80 percent in response to the 100 - year recurrence interval, 30-minute duration rainfall event. These subbasins ranged in size from 0.01 to 6.57 square kilometers (km 2 ), with an average area of 0.29 km 2 , and were mostly steep, upstream tributaries to larger channels in the area. Modeled debris-flow volumes in 465 subbasins were greater than 10,000 cubic meters (m 3 ), and 14 of those subbasins had modeled debris‑flow volumes greater than 100,000 m 3 . The rankings of integrated relative debris - flow hazard indexes for each subbasin were generated by multiplying the individual subbasin values for debris - flow volume, debris‑flow probability, and average burn probability. The subbasins with integrated hazard index values in the top 2 percent typically are large, upland tributaries to canyons and channels primarily in the Upper Rio Grande and Rio Grande-Santa Fe watershed areas. No subbasins in this group have basin areas less than 1.0 km 2 . Many of these areas already had significant mass‑wasting episodes following the Las Conchas Fire in 2011. Other subbasins with integrated hazard index values in the top 2 percent are scattered throughout the Jemez River watershed area, including some subbasins in the interior of the Valles Caldera. Only a few subbasins in the top integrated hazard index group are in the Rio Chama watershed area. This prewildfire assessment approach is valuable to resource managers because the analysis of the debris-flow threat is made before a wildfire occurs, which facilitates prewildfire management, planning, and mitigation. In north‑central New Mexico, widespread watershed restoration efforts are being done to safeguard vital watersheds against the threat of catastrophic wildfire. This study was designed to help select ideal locations for the restoration efforts that could have the best return on investment.

New Mexico↗
USGS2013

Improving regression-model-based streamwater constituent load estimates derived from serially correlated data

A regression-model based approach is a commonly used, efficient method for estimating streamwater constituent load when there is a relationship between streamwater constituent concentration and continuous variables such as streamwater discharge, season and time. A subsetting experiment using a 30-year dataset of daily suspended sediment observations from the Mississippi River at Thebes, Illinois, was performed to determine optimal sampling frequency, model calibration period length, and regression model methodology, as well as to determine the effect of serial correlation of model residuals on load estimate precision. Two regression-based methods were used to estimate streamwater loads, the Adjusted Maximum Likelihood Estimator (AMLE), and the composite method, a hybrid load estimation approach. While both methods accurately and precisely estimated loads at the model’s calibration period time scale, precisions were progressively worse at shorter reporting periods, from annually to monthly. Serial correlation in model residuals resulted in observed AMLE precision to be significantly worse than the model calculated standard errors of prediction. The composite method effectively improved upon AMLE loads for shorter reporting periods, but required a sampling interval of at least 15-days or shorter, when the serial correlations in the observed load residuals were greater than 0.15. AMLE precision was better at shorter sampling intervals and when using the shortest model calibration periods, such that the regression models better fit the temporal changes in the concentration–discharge relationship. The models with the largest errors typically had poor high flow sampling coverage resulting in unrepresentative models. Increasing sampling frequency and/or targeted high flow sampling are more efficient approaches to ensure sufficient sampling and to avoid poorly performing models, than increasing calibration period length.

Illinois↗
USGS1979

A digital model for simulation of ground-water hydrology in the Houston area, Texas

This report documents the construction and calibration of a digital model for the simulation of hydrologic conditions in the Chicot and Evangeline aquifers in the Houston area of southeastern Texas. The model is a five-layer finite-difference model, with a grid pattern of 63 x 67 nodes representing an area of 27,000 square miles, for simulation of three- dimensional ground-water flow. The hydrologic properties and processes modeled were ground-water withdrawals, transmissivities, storage coefficients of the aquifers and clays, quantity of water derived from storage in the clays, and vertical hydraulic conductivity and vertical leakage. The model, which simulates water-level declines, changes in storage in the clay layers, and land-surface subsidence, was calibrated by use of historical records from 1890 to 1975. It is very sensitive to variations in transmissivities and to variations in water-table and artesian storage. It is less sensitive to variations in clay storage.

Texas↗
USGS2000

Application of nonlinear least-squares regression to ground-water flow modeling, west-central Florida

A nonlinear least-squares regression technique for estimation of ground-water flow model parameters was applied to an existing model of the regional aquifer system underlying west-central Florida. The regression technique minimizes the differences between measured and simulated water levels. Regression statistics, including parameter sensitivities and correlations, were calculated for reported parameter values in the existing model. Optimal parameter values for selected hydrologic variables of interest are estimated by nonlinear regression. Optimal estimates of parameter values are about 140 times greater than and about 0.01 times less than reported values. Independently estimating all parameters by nonlinear regression was impossible, given the existing zonation structure and number of observations, because of parameter insensitivity and correlation. Although the model yields parameter values similar to those estimated by other methods and reproduces the measured water levels reasonably accurately, a simpler parameter structure should be considered. Some possible ways of improving model calibration are to: (1) modify the defined parameter-zonation structure by omitting and/or combining parameters to be estimated; (2) carefully eliminate observation data based on evidence that they are likely to be biased; (3) collect additional water-level data; (4) assign values to insensitive parameters, and (5) estimate the most sensitive parameters first, then, using the optimized values for these parameters, estimate the entire data set.

Water-Resources Investigations Report↗
USGS2000

MODFLOW-2000, the U.S. Geological Survey Modular Ground-Water Model -Documentation of the Hydrogeologic-Unit Flow (HUF) Package

This report documents the Hydrogeologic-Unit Flow (HUF) Package for the groundwater modeling computer program MODFLOW-2000. The HUF Package is an alternative internal flow package that allows the vertical geometry of the system hydrogeology to be defined explicitly within the model using hydrogeologic units that can be different than the definition of the model layers. The HUF Package works with all the processes of MODFLOW-2000. For the Ground-Water Flow Process, the HUF Package calculates effective hydraulic properties for the model layers based on the hydraulic properties of the hydrogeologic units, which are defined by the user using parameters. The hydraulic properties are used to calculate the conductance coefficients and other terms needed to solve the ground-water flow equation. The sensitivity of the model to the parameters defined within the HUF Package input file can be calculated using the Sensitivity Process, using observations defined with the Observation Process. Optimal values of the parameters can be estimated by using the Parameter-Estimation Process. The HUF Package is nearly identical to the Layer-Property Flow (LPF) Package, the major difference being the definition of the vertical geometry of the system hydrogeology. Use of the HUF Package is illustrated in two test cases, which also serve to verify the performance of the package by showing that the Parameter-Estimation Process produces the true parameter values when exact observations are used.

Open-File Report↗
USGS2026

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗
USGS2012

Dissolved oxygen analysis, TMDL model comparison, and particulate matter shunting—Preliminary results from three model scenarios for the Klamath River upstream of Keno Dam, Oregon

Efforts are underway to identify actions that would improve water quality in the Link River to Keno Dam reach of the Upper Klamath River in south-central Oregon. To provide further insight into water-quality improvement options, three scenarios were developed, run, and analyzed using previously calibrated CE-QUAL-W2 hydrodynamic and water-quality models. Additional scenarios are under development as part of this ongoing study. Most of these scenarios evaluate changes relative to a "current conditions" model, but in some cases a "natural conditions" model was used that simulated the reach without the effect of point and nonpoint sources and set Upper Klamath Lake at its Total Maximum Daily Load (TMDL) targets. These scenarios were simulated using a model developed by the U.S. Geological Survey (USGS) and Watercourse Engineering, Inc. for the years 2006–09, referred to here as the "USGS model." Another model of the reach was developed by Tetra Tech, Inc. for years 2000 and 2002 to support the Klamath River TMDL process; that model is referred to here as the "TMDL model." The three scenarios described in this report included (1) an analysis of whether this reach of the Upper Klamath River would be in compliance with dissolved oxygen standards if sources met TMDL allocations, (2) an application of more recent datasets to the TMDL model with comparison to results from the USGS model, and (3) an examination of the effect on dissolved oxygen in the Klamath River if particulate material were stopped from entering Klamath Project diversion canals. Updates and modifications to the USGS model are in progress, so in the future these scenarios will be reanalyzed with the updated model and the interim results presented here will be superseded. Significant findings from this phase of the investigation include: * The TMDL analysis used depth-averaged dissolved oxygen concentrations from model output for comparison with dissolved oxygen standards. The Oregon dissolved oxygen standards do not specify whether the numeric criteria are based on depth-averaged dissolved oxygen concentration; this was an interpretation of the standards rule by the Oregon Department of Environmental Quality (ODEQ). In this study, both depth-averaged and volume-averaged dissolved oxygen concentrations were calculated from model output. Results showed that modeled depth-averaged concentrations typically were lower than volume-averaged dissolved oxygen concentrations because depth-averaging gives a higher weight to small volume areas near the channel bottom that often have lower dissolved oxygen concentrations. Results from model scenarios in this study are reported using volume-averaged dissolved oxygen concentrations. * Under all scenarios analyzed, violations of the dissolved oxygen standard occurred most often in summer. Of the three dissolved oxygen criteria that must be met, the 30-day standard was violated most frequently. Under the base case (current conditions), fewer violations occurred in the upstream part of the reach. More violations occurred in the down-stream direction, due in part to oxygen demand from the decay of algae and organic matter from Link River and other inflows. * A condition in which Upper Klamath Lake and its Link River outflow achieved Upper Klamath Lake TMDL water-quality targets was most effective in reducing the number of violations of the dissolved oxygen standard in the Link River to Keno Dam reach of the Klamath River. The condition in which point and nonpoint sources within the Link River to Keno Dam reach met Klamath River TMDL allocations had no effect on dissolved oxygen compliance in some locations and a small effect in others under current conditions. On the other hand, meeting TMDL allocations for nonpoint and point sources was predicted to be important in meeting dissolved oxygen criteria when Upper Klamath Lake and Link River also met Upper Klamath TMDL water-quality targets. * The location of greatest dissolved oxygen improvement from nutrient and organic matter reductions was downstream from point and nonpoint source inflows because time and distance are required for decay to occur and for oxygen demand to be exerted. * After assessing compliance with dissolved oxygen standards at all 102 model segments in the Link River to Keno Dam reach, it was determined that the seven locations used by ODEQ appear to be a representative subset of the reach for dissolved oxygen analysis. * The USGS and TMDL models were qualitatively compared by running both models for the 2006–09 period but preserving the essential characteristics of each, such as organic matter partitioning, bathymetric representation, and parameter rates. The analysis revealed that some constituents were not greatly affected by the differing algorithms, rates, and assumptions in the two models. Conversely, other constituents, especially organic matter, were simulated differently by the two models. Organic matter in this river system is best represented by a mixture of relatively labile particulate material and a substantial concentration of refractory dissolved material. In addition, the use of a first-order sediment oxygen demand, as in the USGS model, helps to capture the seasonal and dynamic effect of settled organic and algal material. * Simulation of shunting (diverting) particulate material away from the intake of four Klamath Project diversion canals, so that the material stayed in the river and out of the Project area, caused higher concentrations of particulate material to occur in the river. In all cases modeled, the increase in in-river particulate material also produced decreased dissolved oxygen concentrations and an increase in the number of days when dissolved oxygen standards were violated. * If particulate material were shunted back into the river at the Klamath Project diversion canals, less organic matter and nutrients would be taken into the Klamath Project area and the Lost River basin, resulting in return flows to the Klamath River via Lost River Diversion Channel that may have reduced nutrient concentrations. Model scenarios bracketing potential end-member nutrient concentrations showed that the composition of the return flows had little to no effect on dissolved oxygen compliance under simulated conditions.

Oregon↗
USGS2009

Geologic Criteria for the Assessment of Sedimentary Exhalative (Sedex) Zn-Pb-Ag Deposits

Sedex deposits account for more than 50 percent of the world's zinc and lead reserves and furnish more than 25 percent of the world's production of these two metals. This report draws on previous syntheses as well as on topical studies of deposits in sedex basins to determine the characteristics and processes that produced sedex deposits. This analysis also uses studies of the tectonic, sedimentary, and fluid evolution of modern and ancient sedimentary basins and mass balance constraints to identify the hydrothermal processes that are required to produce sedex deposits. This report demonstrates how a genetic model can be translated into geologic criteria that can be used in the U.S. Geological Survey National Assessments for sedex zinc-lead-silver deposits to define permissive tracts, assess the relative prospectivity of permissive tracts, and map favorability within permissive tracts.

Open-File Report↗
USGS2011

Development of a precipitation-runoff model to simulate unregulated streamflow in the South Fork Flathead River Basin, Montana

This report documents the development of a precipitation-runoff model for the South Fork Flathead River Basin, Mont. The Precipitation-Runoff Modeling System model, developed in cooperation with the Bureau of Reclamation, can be used to simulate daily mean unregulated streamflow upstream and downstream from Hungry Horse Reservoir for water-resources planning. Two input files are required to run the model. The time-series data file contains daily precipitation data and daily minimum and maximum air-temperature data from climate stations in and near the South Fork Flathead River Basin. The parameter file contains values of parameters that describe the basin topography, the flow network, the distribution of the precipitation and temperature data, and the hydrologic characteristics of the basin soils and vegetation. A primary-parameter file was created for simulating streamflow during the study period (water years 1967-2005). The model was calibrated for water years 1991-2005 using the primary-parameter file. This calibration was further refined using snow-covered area data for water years 2001-05. The model then was tested for water years 1967-90. Calibration targets included mean monthly and daily mean unregulated streamflow upstream from Hungry Horse Reservoir, mean monthly unregulated streamflow downstream from Hungry Horse Reservoir, basin mean monthly solar radiation and potential evapotranspiration, and daily snapshots of basin snow-covered area. Simulated streamflow generally was in better agreement with observed streamflow at the upstream gage than at the downstream gage. Upstream from the reservoir, simulated mean annual streamflow was within 0.0 percent of observed mean annual streamflow for the calibration period and was about 2 percent higher than observed mean annual streamflow for the test period. Simulated mean April-July streamflow upstream from the reservoir was about 1 percent lower than observed streamflow for the calibration period and about 4 percent higher than observed for the test period. Downstream from the reservoir, simulated mean annual streamflow was 17 percent lower than observed streamflow for the calibration period and 12 percent lower than observed streamflow for the test period. Simulated mean April-July streamflow downstream from the reservoir was 13 percent lower than observed streamflow for the calibration period and 6 percent lower than observed streamflow for the test period. Calibrating to solar radiation, potential evapotranspiration, and snow-covered area improved the model representation of evapotranspiration, snow accumulation, and snowmelt processes. Simulated basin mean monthly solar radiation values for both the calibration and test periods were within 9 percent of observed values except during the month of December (28 percent different). Simulated basin potential evapotranspiration values for both the calibration and test periods were within 10 percent of observed values except during the months of January (100 percent different) and February (13 percent different). The larger percent errors in simulated potential evaporation occurred in the winter months when observed potential evapotranspiration values were very small; in January the observed value was 0.000 inches and in February the observed value was 0.009 inches. Simulated start of melting of the snowpack occurred at about the same time as observed start of melting. The simulated snowpack accumulated to 90-100 percent snow-covered area 1 to 3 months earlier than observed snowpack. This overestimated snowpack during the winter corresponded to underestimated streamflow during the same period. In addition to the primary-parameter file, four other parameter files were created: for a "recent" period (1991-2005), a historical period (1967-90), a "wet" period (1989-97), and a "dry" period (1998-2005). For each data file of projected precipitation and air temperature, a single parameter file can be used to simulate a s

Montana;Idaho↗
USGS1983

Results of an adaptive environmental assessment modeling workshop concerning potential impacts of drilling muds and cuttings on the marine environment

Drilling fluids or "muds" are essential components of modern drilling operations. They provide integrity for the well bore, a medium for removal of formation cuttings, and lubrication and cooling of the drill bit and pipe. The modeling workshop described in this report was conducted September 14-18, 1981 in Gulf Breeze, Florida to consider potential impacts of discharged drilling muds and cuttings on the marine environment. The broad goals of the workshop were synthesis of information on fate and effects, identification of general relationships between drilling fluids and the marine environment, and identification of site-specific variables likely to determine impacts of drilling muds and cuttings in various marine sites. The workshop was structured around construction of a model simulating fate and effects of discharges from a single rig into open water areas of the Gulf of Mexico, and discussion of factors that might produce different fate and effects in enclosed areas such as bays and estuaries. The simulation model was composed of four connected submodels. A Discharge/Fate submodel dealt with the discharge characteristics of the rig and the subsequent fate of discharged material. Three effects submodels then calculated biological responses at distances away from the rig for the water column, soft bottom benthos (assuming the rig was located over a soft bottom environment), and hard bottom benthos (assuming the rig was located over a hard bottom environment). The model focused on direct linkages between the discharge and various organisms rather than on how the marine ecosystem itself is interconnected. Behavior of the simulation model indicated relatively localized effects of drilling muds and cuttings discharged from a single platform into open water areas. Water column fate and effects were dominated by rapid dilution. Effects from deposition of spent mud and cuttings were spatially limited with relatively rapid recovery, especially in soft bottom benthic communities which were conceptualized as being adapted to frequent storms. This behavior was generated by the set of assumptions about linkages and functional relationships used to construct the model. Areas of uncertainty included methods for extrapolating 96-hr LC50 so results to exposures of varying lengths and concentrations; recovery rates of benthic communities; responses to various depths and rates of burial; fate and effects of the plume in relationship to stratification layers; and long-term and sub-lethal effects of slightly elevated concentrations of discharged materials. Evaluation of the assumptions of the Soft Bottom Submodel suggest that the assumptions used may have been relatively liberal estimates of resiliency of these communities. Discussion of "closed" water bodies such as bays and estuaries indicated several reasons to expect different and more complex fate and effects behavior in these areas. These factors included different species and communities (such as aquatic macrophytes and oyster beds), more complex circulation and stratification patterns, and potentially more active resuspension processes. Much of the possible difference in behavior in these areas centers around the extent to which they are “closed” or in the relative residence times of water and sediments in these areas as they determine the long-term dispersion of discharged material. Despite the complexity and variability of these areas, a large body of knowledge (such as that concerning fate and physical effects of dredge spoil) that could be effectively employed in analysis of potential fate and physical effects in enclosed areas was identified.

Report↗
USGS2017

P- and S-wave velocity models incorporating the Cascadia subduction zone for 3D earthquake ground motion simulations, Version 1.6—Update for Open-File Report 2007–1348

In support of earthquake hazard studies and ground motion simulations in the Pacific Northwest, three-dimensional P- and S-wave velocity ( V P and V S , respectively) models incorporating the Cascadia subduction zone were previously developed for the region encompassed from about 40.2°N. to 50°N. latitude, and from about 122°W. to 129°W. longitude. This report describes updates to the Cascadia velocity property volumes of model version 1.3 (V1.3), herein called model version 1.6 (V1.6). As in model V1.3, the updated V1.6 model volume includes depths from 0 kilometers (mean sea level) to 60 kilometers, and it is intended to be a reference for researchers who have used, or are planning to use, this model in their Earth science investigations. To this end, it is intended that the V P and V S property volumes of model V1.6 will be considered a template for a community velocity model of the Cascadia region as additional results become available. With the recent and ongoing development of the National Crustal Model, we envision any future versions of this model will be directly integrated with that effort.

Oregon, Vancouver, Washington↗
USGS2011

Regional groundwater-flow model of the Redwall-Muav, Coconino, and alluvial basin aquifer systems of northern and central Arizona

A numerical flow model (MODFLOW) of the groundwater flow system in the primary aquifers in northern Arizona was developed to simulate interactions between the aquifers, perennial streams, and springs for predevelopment and transient conditions during 1910 through 2005. Simulated aquifers include the Redwall-Muav, Coconino, and basin-fill aquifers. Perennial stream reaches and springs that derive base flow from the aquifers were simulated, including the Colorado River, Little Colorado River, Salt River, Verde River, and perennial reaches of tributary streams. Simulated major springs include Blue Spring, Del Rio Springs, Havasu Springs, Verde River headwater springs, several springs that discharge adjacent to major Verde River tributaries, and many springs that discharge to the Colorado River. Estimates of aquifer hydraulic properties and groundwater budgets were developed from published reports and groundwater-flow models. Spatial extents of aquifers and confining units were developed from geologic data, geophysical models, a groundwater-flow model for the Prescott Active Management Area, drill logs, geologic logs, and geophysical logs. Spatial and temporal distributions of natural recharge were developed by using a water-balance model that estimates recharge from direct infiltration. Additional natural recharge from ephemeral channel infiltration was simulated in alluvial basins. Recharge at wastewater treatment facilities and incidental recharge at agricultural fields and golf courses were also simulated. Estimates of predevelopment rates of groundwater discharge to streams, springs, and evapotranspiration by phreatophytes were derived from previous reports and on the basis of streamflow records at gages. Annual estimates of groundwater withdrawals for agriculture, municipal, industrial, and domestic uses were developed from several sources, including reported withdrawals for nonexempt wells, estimated crop requirements for agricultural wells, and estimated per capita water use for exempt wells. Accuracy of the simulated groundwater-flow system was evaluated by using observational control from water levels in wells, estimates of base flow from streamflow records, and estimates of spring discharge. Major results from the simulations include the importance of variations in recharge rates throughout the study area and recharge along ephemeral and losing stream reaches in alluvial basins. Insights about the groundwater-flow systems in individual basins include the hydrologic influence of geologic structures in some areas and that stream-aquifer interactions along the lower part of the Little Colorado River are an effective control on water level distributions throughout the Little Colorado River Plateau basin. Better information on several aspects of the groundwater flow system are needed to reduce uncertainty of the simulated system. Many areas lack documentation of the response of the groundwater system to changes in withdrawals and recharge. Data needed to define groundwater flow between vertically adjacent water-bearing units is lacking in many areas. Distributions of recharge along losing stream reaches are poorly defined. Extents of aquifers and alluvial lithologies are poorly defined in parts of the Big Chino and Verde Valley sub-basins. Aquifer storage properties are poorly defined throughout most of the study area. Little data exist to define the hydrologic importance of geologic structures such as faults and fractures. Discharge of regional groundwater flow to the Verde River is difficult to identify in the Verde Valley sub-basin because of unknown contributions from deep percolation of excess surface water irrigation.

Scientific Investigations Report↗
USGS2016

Model description and evaluation of the mark-recapture survival model used to parameterize the 2012 status and threats analysis for the Florida manatee ( Trichechus manatus latirostris )

This report provides supporting details and evidence for the rationale, validity and efficacy of a new mark-recapture model, the Barker Robust Design, to estimate regional manatee survival rates used to parameterize several components of the 2012 version of the Manatee Core Biological Model (CBM) and Threats Analysis (TA). The CBM and TA provide scientific analyses on population viability of the Florida manatee subspecies ( Trichechus manatus latirostris ) for U.S. Fish and Wildlife Service’s 5-year reviews of the status of the species as listed under the Endangered Species Act. The model evaluation is presented in a standardized reporting framework, modified from the TRACE (TRAnsparent and Comprehensive model Evaluation) protocol first introduced for environmental threat analyses. We identify this new protocol as TRACE-MANATEE SURVIVAL and this model evaluation specifically as TRACE-MANATEE SURVIVAL, Barker RD version 1. The longer-term objectives of the manatee standard reporting format are to (1) communicate to resource managers consistent evaluation information over sequential modeling efforts; (2) build understanding and expertise on the structure and function of the models; (3) document changes in model structures and applications in response to evolving management objectives, new biological and ecological knowledge, and new statistical advances; and (4) provide greater transparency for management and research review.

Open-File Report↗
USGS2021

Joint species distribution models of Everglades wading birds to inform restoration planning

Restoration of the Florida Everglades, a substantial wetland ecosystem within the United States, is one of the largest ongoing restoration projects in the world. Decision-makers and managers within the Everglades ecosystem rely on ecological models forecasting indicator wildlife response to changes in the management of water flows within the system. One such indicator of ecosystem health, the presence of wading bird communities on the landscape, is currently assessed using three species distribution models that assume perfect detection and report output on different scales that are challenging to compare against one another. We sought to use current advancements in species distribution modeling to improve models of Everglades wading bird distribution. Using a joint species distribution model that accounted for imperfect detection, we modeled the presence of nine species of wading bird simultaneously in response to annual hydrologic conditions and landscape characteristics within the Everglades system. Our resulting model improved upon the previous model in three key ways: 1) the model predicts probability of occupancy for the nine species on a scale of 0–1, making the output more intuitive and easily comparable for managers and decision-makers that must consider the responses of several species simultaneously; 2) through joint species modeling, we were able to consider rarer species within the modeling that otherwise are detected in too few numbers to fit as individual models; and 3) the model explicitly allows detection probability of species to be less than 1 which can reduce bias in the site occupancy estimates. These improvements are essential as Everglades restoration continues and managers require models that consider the impacts of water management on key indicator wildlife such as the wading bird community.

Florida↗
USGS2015

Development of conceptual ecological models linking management of the Missouri River to pallid sturgeon population dynamics

This report documents the process of developing and refining conceptual ecological models (CEMs) for linking river management to pallid sturgeon ( Scaphirhynchus albus ) population dynamics in the Missouri River. The refined CEMs are being used in the Missouri River Pallid Sturgeon Effects Analysis to organize, document, and formalize an understanding of pallid sturgeon population responses to past and future management alternatives. The general form of the CEMs, represented by a population-level model and component life-stage models, was determined in workshops held in the summer of 2013. Subsequently, the Missouri River Pallid Sturgeon Effects Analysis team designed a general hierarchical structure for the component models, refined the graphical structure, and reconciled variation among the components and between models developed for the upper river (Upper Missouri & Yellowstone Rivers) and the lower river (Missouri River downstream from Gavins Point Dam). Importance scores attributed to the relations between primary biotic characteristics and survival were used to define a candidate set of working dominant hypotheses about pallid sturgeon population dynamics. These CEMs are intended to guide research and adaptive-management actions to benefit pallid sturgeon populations in the Missouri River.

Missouri River, Yellowstone River↗
USGS2002

Distinguishing sediment waves from slope failure deposits: Field examples, including the 'humboldt slide', and modelling results

Migrating sediment waves have been reported in a variety of marine settings, including submarine levee-fan systems, floors of fjords, and other basin or continental slope environments. Examination of such wave fields reveals nine diagnostic characteristics. When these characteristics are applied to several features previously attributed to submarine landslide deformation, they suggest that the features should most likely be reinterpreted as migrating sediment-wave fields. Sites that have been reinterpreted include the 'Humboldt slide' on the Eel River margin in northern California, the continental slope in the Gulf of Cadiz, the continental shelf off the Malaspina Glacier in the Gulf of Alaska, and the Adriatic shelf. A reassessment of all four features strongly suggests that numerous turbidity currents, separated by intervals of ambient hemipelagic sedimentation, deposited the wave fields over thousands of years. A numerical model of hyperpycnal discharge from the Eel River, for example, shows that under certain alongshore-current conditions, such events can produce turbidity currents that flow across the 'Humboldt slide', serving as the mechanism for the development of migrating sediment waves. Numerical experiments also demonstrate that where a series of turbidity currents flows across a rough seafloor (i.e. numerical steps), sediment waves can form and migrate upslope. Hemipelagic sedimentation between turbidity current events further facilitates the upslope migration of the sediment waves. Physical modelling of turbidity currents also confirms the formation and migration of seafloor bedforms. The morphologies of sediment waves generated both numerically and physically in the laboratory bear a strong resemblance to those observed in the field, including those that were previously described as submarine landslides.

Marine Geology↗
USGS2017

Assessing models of arsenic occurrence in drinking water from bedrock aquifers in New Hampshire

Three existing multivariate logistic regression models were assessed using new data to evaluate the capacity of the models to correctly predict the probability of groundwater arsenic concentrations exceeding the threshold values of 1, 5, and 10 micrograms per liter (µg/L) in New Hampshire, USA. A recently released testing dataset includes arsenic concentrations from groundwater samples collected in 2004–2005 from a mix of 367 public-supply and private domestic wells. The use of this dataset to test three existing logistic regression models demonstrated enhanced overall predictive accuracy for the 5 and 10 μg/L models. Overall accuracies of 54.8, 76.3, and 86.4 percent were reported for the 1, 5, and 10 μg/L models, respectively. The state was divided by counties into northwest and southeast regions. Regional differences in accuracy were identified; models had an average accuracy of 83.1 percent for the counties in the northwest and 63.7 percent in the southeast. This is most likely due to high model specificity in the northwest and regional differences in arsenic occurrence. Though these models have limitations, they allow for arsenic hazard assessment across the region. The introduction of well-type (public or private), well depth, and casing length as explanatory variables may be appropriate measures to improve model performance. Our findings indicate that the original models generalize to the testing dataset, and should continue to serve as an important vehicle of preventative public health that may be applied to other groundwater contaminants in New Hampshire.

New Hampshire↗