SEARCH · USGS Science

Results for “Computer contribution”

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

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25Linked to original sources

USGS2008

Simulation of the Quantity, Variability, and Timing of Streamflow in the Dennys River Basin, Maine, by Use of a Precipitation-Runoff Watershed Model

The U.S. Geological Survey (USGS), in cooperation with the Maine Department of Marine Resources Bureau of Sea Run Fisheries and Habitat, began a study in 2004 to characterize the quantity, variability, and timing of streamflow in the Dennys River. The study included a synoptic summary of historical streamflow data at a long-term streamflow gage, collecting data from an additional four short-term streamflow gages, and the development and evaluation of a distributed-parameter watershed model for the Dennys River Basin. The watershed model used in this investigation was the USGS Precipitation-Runoff Modeling System (PRMS). The Geographic Information System (GIS) Weasel was used to delineate the Dennys River Basin and subbasins and derive parameters for their physical geographic features. Calibration of the models used in this investigation involved a four-step procedure in which model output was evaluated against four calibration data sets using computed objective functions for solar radiation, potential evapotranspiration, annual and seasonal water budgets, and daily streamflows. The calibration procedure involved thousands of model runs and was carried out using the USGS software application Luca (Let us calibrate). Luca uses the Shuffled Complex Evolution (SCE) global search algorithm to calibrate the model parameters. The SCE method reliably produces satisfactory solutions for large, complex optimization problems. The primary calibration effort went into the Dennys main stem watershed model. Calibrated parameter values obtained for the Dennys main stem model were transferred to the Cathance Stream model, and a similar four-step SCE calibration procedure was performed; this effort was undertaken to determine the potential to transfer modeling information to a nearby basin in the same region. The calibrated Dennys main stem watershed model performed with Nash-Sutcliffe efficiency (NSE) statistic values for the calibration period and evaluation period of 0.79 and 0.76, respectively. The Cathance Stream model had an NSE value of 0.68. The Dennys River Basin models make use of limited streamflow-gaging station data and provide information to characterize subbasin hydrology. The calibrated PRMS watershed models of the Dennys River Basin provide simulated daily streamflow time series from October 1, 1985, through September 30, 2006, for nearly any location within the basin. These models enable natural-resources managers to characterize the timing and quantity of water moving through the basin to support many endeavors including geochemical calculations, water-use assessment, Atlantic salmon population dynamics and migration modeling, habitat modeling and assessment, and other resource-management scenario evaluations. Characterizing streamflow contributions from subbasins in the basin and the relative amounts of surface- and ground-water contributions to streamflow throughout the basin will lead to a better understanding of water quantity and quality in the basin. Improved water-resources information will support Atlantic salmon protection efforts.

Scientific Investigations Report↗

USGS2024

Methods to evaluate and improve the modeling of rupture directivity in assessment of seismic hazard

In recent years, there have been several advancements related to the modelling of near-source effects of earthquake rupture on strong ground shaking, leading to an improved characterization of ground motions and resulting seismic hazard. Some of these modifications have stemmed from physics-based numerical modelling of the earthquake rupture process, using physics-based dynamic rupture simulations. These contributions have led to a better understanding of how fault rupture characteristics, geometry, and the style of faulting can interact with the hypocenter-dependence on the path from source to site that may ultimately guide the development of seismic directivity models. Moving forward, the application of modern techniques can be used to incorporate these source characteristics and near-fault ground motion behavior that contribute to the azimuthally varying effects that result in rupture directivity. One example is the application of machine learning methods to support more automated integration of new predictor variables in model development and open more evaluation opportunities to access residuals. Here, we utilize several techniques to take advantage of the plethora of synthetic data and its ability to supplement preexisting trends observed in data. We showcase two examples of how models can be either developed, expanded upon, or constrained using artificial neural network model (ANNs). We evaluate the performance of the ANN with existing methods, comparing misfit, potential limitations, and ability to continue to improve upon these methods in the future. One approach uses a set of simulations with corresponding synthetic ground motions from the Southern California Earthquake Center (SCEC) CyberShake study to develop a ground motion model adapted to incorporate seismic directivity information using an ANN. This large database (TBs) enables us to train the model to capture magnitude, period, and distance variations and how these parameters relate to amplification from hypocenters located along finite-faults. In some cases, there is reduced misfit from better representing source features that aren’t included in base ground motion models that neglect hypocenter location (e.g. azimuthal variation, source-to-site terms). Another ANN method uses a shallow-layered neural network model to better fit a hypocenter-independent model. This method adjusts the median and aleatory variability to account for the averaged impact of various hypocenter distributions to fit the underlying directivity adjustment model. This method serves as a template to apply to other directivity models, improving computational efficiency and more readily enabling integration in hazard codes.

California↗

USGS2009

Hydrodynamic simulations of physical aquatic habitat availability for Pallid Sturgeon in the Lower Missouri River, at Yankton, South Dakota, Kenslers Bend, Nebraska, Little Sioux, Iowa, and Miami, Missouri, 2006-07

The objective of this study was to assess the sensitivity of habitat availability in the Lower Missouri River to discharge variation, with emphasis on habitats that might support spawning of the endangered pallid sturgeon. We constructed computational hydrodynamic models for four reaches that were selected because of evidence that sturgeon have spawned in them. The reaches are located at Miami, Missouri (river mile 259.6–263.5), Little Sioux, Iowa (river mile 669.6–673.5), Kenslers Bend, Nebraska (river mile 743.9–748.1), and Yankton, South Dakota reach (river mile 804.8–808.4). The models were calibrated for a range of measured flow conditions, and run for a range of discharges that might be affected by flow modifications from Gavins Point Dam. Model performance was assessed by comparing modeled and measured water velocities. A selection of derived habitat units was assessed for sensitivity to hydraulic input parameters (drag coefficient and lateral eddy viscosity). Overall, model results were minimally sensitive to varying eddy viscosity; varying lateral eddy viscosity by 20 percent resulted in maximum change in habitat units of 5.4 percent. Shallow-water habitat units were most sensitive to variation in drag coefficient with 42 percent change in unit area resulting from 20 percent change in the parameter value; however, no habitat unit value changed more than 10 percent for a 10 percent variation in drag coefficient. Sensitivity analysis provides guidance for selecting habitat metrics that maximize information content while minimizing model uncertainties. To assess model sensitivities arising from topographic variation from sediment transport on an annual time scale, we constructed separate models from two complete independent surveys in 2006 and 2007. The net topographic change was minimal at each site; the ratio of net topographic change to water volume in the reaches at 95 percent exceedance flow was less than 5 percent, indicating that on a reach-average basis, annual topographic change contributed little to habitat area variation. Net erosion occurred at Yankton (the upstream reach) and because erosion was distributed uniformly, there was little affect on many habitat metrics. Topographic change was spatially nonuniform at Little Sioux and Kenslers Bend reaches. Shallow water habitat units and some reach-scale patch statistics (edge density, patch density, and Simpson’s Diversity Index) were affected by these changes. Erosion dominated at the downstream reach but habitat metrics did not vary substantially from 2006 to 2007. Among habitat metrics that were explored, zones of convergent flow were identified as areas that most closely correspond to spawning habitats of other sturgeon species, as identified in the scientific literature, and that are consistent with sparse data on pallid sturgeon spawning locations in the Lower Missouri River. Areas of convergent zone habitat varied little with discharges that would be associated with spring pulsed flows, and relations with discharge changed negligibly between 2006 and 2007. Other habitat measures show how physical habitat varies with discharge and among the four reaches. Wake habitats defined by velocity gradients seem to correspond with migration pathways of adult pallid sturgeon. Habitats with low Froude-number correspond to low energy areas that may accumulate passively transporting particles, organic matter, and larval fish. Among the modeled reaches, Yankton had substantially longer water residence time for equivalent flow exceedances than the other three modeled reaches. Longer residence times result from greater flow resistance in the relatively wide, shallow channel and may be associated with longer residence times of passively transported particulate materials.

Scientific Investigations Report↗

USGS2019

A comparison of hydrocarbon-related landscape disturbance patterns along the New York-Pennsylvania border, 2004–2013

Executive Summary The New York-Pennsylvania area has a long history of hydrocarbon extraction, and the addition of shale gas extraction methods contributes to landscape disturbance borne by previously developed oil and non-shale gas resources. The main unconventional extraction method used to extract shale gas from the Marcellus Shale located in New York and Pennsylvania is hydraulic fracturing, or “fracking,” although other conventional methods are used extensively. All forms of hydrocarbon extraction disturb the surrounding landscape to some extent, primarily in the form of land clearance and degradation, road construction, and pipeline development, although the effects of these disturbances are not fully understood. In this study, landscape-change metrics and indicators are used to analyze change in a 10-county region along the New York-Pennsylvania border—the New York counties of Allegany, Steuben, Chemung, Tioga, and Broome, and the Pennsylvania counties of McKean, Potter, Tioga, Bradford, and Susquehanna. This 10-county region was selected due to the differences in policies between the States of New York and Pennsylvania. While fracking occurred extensively in Pennsylvania over the past 10 years or more, the State of New York issued a temporary moratorium against hydraulic fracturing in 2010—citing repercussions that might affect air quality, water quality, and public health—and officially banned hydraulic fracturing in June 2015. The quantification of landscape disturbance due to hydrocarbon extraction activities is presented in this report as land-use and land-cover (LULC) change between 2004 and 2013 and defined using specific disturbance categories (including well sites, roads, and pipelines) to compare the disturbances and changes, by county, on both sides of the New York-Pennsylvania border. The quantification was accomplished by gathering the signatures of disturbance from high-resolution aerial images, comparing the derived totals of disturbance, and then computing landscape metrics in a geographic information system (GIS) environment. The collected data represent a summation of landscape disturbance from oil and gas development, as some of the data represented were established decades earlier. The Analytical Tools Interface for Landscape Assessments (ATtILA) software was used to calculate land-cover area and landscape metrics for each shale gas, non-shale gas, oil, and other infrastructure types associated with hydrocarbons across each county and both five-county regions in the study area. The three primary metrics used to describe changes in forest structure were (1)forest area, (2) interior forest area, and (3) forest edge area. The changes in metrics were subsequently evaluated using the Pearson correlation coefficient. Overall, the disturbed-area footprint in the Pennsylvania region is considerably larger than the disturbed-area footprint in the New York region (13,687.9 hectares [ha] in Pennsylvania; 3,840.5 ha in New York). Disturbance per site is similar, with 1.2 disturbed ha per site in New York and 1.6 disturbed ha per site in Pennsylvania. In the New York-Pennsylvania 10-county region, hydrocarbon-development and extraction disturbance strongly correlate with a reduction in the percentage of forest for the entire region. This observation also appears to be true in the New York five-county region for forest area. This form of disturbance in the New York five-county region shows significantly correlated changes in forest metrics (–0.4 percent total forest area), particularly in the percentage of interior forest (–1.2 percent total area) and forest edge (+0.7 percent total area). On the other hand, gas and hydrocarbon-development and extraction disturbance (1.0 percent total area) in the Pennsylvania five-county region strongly correlates with a total decline in forest area and agricultural land area (–0.8 percent combined total area) but not with either land-cover class separately.

New York, Pennsylvania↗

USGS2001

Total Phosphorus Loads for Selected Tributaries to Sebago Lake, Maine

The streamflow and water-quality datacollection networks of the Portland Water District (PWD) and the U.S. Geological Survey (USGS) as of February 2000 were analyzed in terms of their applicability for estimating total phosphorus loads for selected tributaries to Sebago Lake in southern Maine. The long-term unit-area mean annual flows for the Songo River and for small, ungaged tributaries are similar to the long-term unit-area mean annual flows for the Crooked River and other gaged tributaries to Sebago Lake, based on a regression equation that estimates mean annual streamflows in Maine. Unit-area peak streamflows of Sebago Lake tributaries can be quite different, based on a regression equation that estimates peak streamflows for Maine. Crooked River had a statistically significant positive relation (Kendall's Tau test, p=0.0004) between streamflow and total phosphorus concentration. Panther Run had a statistically significant negative relation (p=0.0015). Significant positive relations may indicate contributions from nonpoint sources or sediment resuspension, whereas significant negative relations may indicate dilution of point sources. Total phosphorus concentrations were significantly larger in the Crooked River than in the Songo River (Wilcoxon rank-sum test, p<0.0001). Evidence was insufficient, however, to indicate that phosphorus concentrations from medium-sized drainage basins, at a significance level of 0.05, were different from each other or that concentrations in small-sized drainage basins were different from each other (Kruskal-Wallis test, p= 0.0980, 0.1265). All large- and medium-sized drainage basins were sampled for total phosphorus approximately monthly. Although not all small drainage basins were sampled, they may be well represented by the small drainage basins that were sampled. If the tributaries gaged by PWD had adequate streamflow data, the current PWD tributary monitoring program would probably produce total phosphorus loading data that would represent all gaged and ungaged tributaries to Sebago Lake. Outside the PWD tributary-monitoring program, the largest ungaged tributary to Sebago Lake contains 1.5 percent of the area draining to the lake. In the absence of unique point or nonpoint sources of phosphorus, ungaged tributaries are unlikely to have total phosphorus concentrations that differ significantly from those in the small tributaries that have concentration data. The regression method, also known as the rating-curve method, was used to estimate the annual total phosphorus load for Crooked River, Northwest River, and Rich Mill Pond Outlet for water years 1996-98. The MOVE.1 method was used to estimate daily streamflows for the regression method at Northwest River and Rich Mill Pond Outlet, where streamflows were not continuously monitored. An averaging method also was used to compute annual loads at the three sites. The difference between the regression estimate and the averaging estimate for each of the three tributaries was consistent with what was expected from previous studies.

Water-Resources Investigations Report↗

USGS2022

Documentation for the Skeletal Storage, Compaction, and Subsidence (CSUB) Package of MODFLOW 6

This report describes the skeletal storage, compaction and subsidence (CSUB) package of MODFLOW 6. The CSUB package simulates the vertical compaction of compressible sediments and land subsidence. The package simulates groundwater storage changes and elastic compaction in coarse-grained aquifer sediments. The CSUB package also simulates groundwater storage changes and elastic and inelastic compaction in fne-grained, compressible interbeds, or in extensive confning units. The package can account for effective stress-dependent changes in storage properties. The CSUB package can also explicitly account for the contribution of water compressibility to groundwater storage changes. Compaction of compressible sediments is formulated using Terzaghi’s elastoplastic model and assumes the total compaction is a small fraction of the total initial thickness of compressible sediments. Compaction is controlled by head or pore-pressure changes and overburden stress changes associated with water-table changes, and thus by effective stress changes within coarse-and fne-grained compressible sediments. If the stress in a compressible unit is less than the preconsolidation stress, compaction is elastic (recoverable). If the stress in a compressible sediment is greater than the preconsolidation stress, compaction is inelastic (irrecoverable) and permanent land subsidence occurs. The propagation of head changes within fne-grained, compressible interbeds is represented numerically using a transient, one-dimensional (vertical) groundwater fow equation. This equation accounts for delayed release of water from storage or uptake of water into storage in the interbeds. Vertical hydraulic conductivity, elastic and inelastic skeletal specifc storage, and interbed thickness control the timing of interbed storage changes. Interbeds that are thin, have a relatively large vertical hydraulic conductivity, or relatively small specifc-storage values equilibrate quickly with heads/pore pressures in surrounding coarse-grained sediments and can be represented as no-delay interbeds that use the simulated groundwater head in a cell to calculate interbed compaction and do not need to be solved numerically using a vertically discretized interbed and the vertical groundwater fow equation. In addition to the applicability to confned groundwater fow systems, several features of the CSUB package make it applicable to shallow, unconfned groundwater fow systems. Geostatic stress can be treated as a function of water-table elevation, and compaction is a function of computed changes in effective stress. The porosity, void ratio, and thickness of shallow and deep coarse-grained aquifer sediments, fne-grained interbeds, and extensive confning units can vary in time based on calculated strain.

Techniques and Methods↗

USGS2008

Mercury in precipitation in Indiana, January 2004–December 2005

Mercury in precipitation was monitored during 2004–2005 at five locations in Indiana as part of the National Atmospheric Deposition Program–Mercury Deposition Network (NADP–MDN). Monitoring stations were operated at Roush Lake near Huntington, Clifty Falls State Park near Madison, Fort Harrison State Park near Indianapolis, Monroe County Regional Airport near Bloomington, and Indiana Dunes National Lakeshore near Porter. At these monitoring stations, precipitation amounts were measured continuously and weekly samples were collected for analysis of mercury by methods achieving detection limits as low as 0.05 ng/L (nanograms per liter). Wet deposition was computed as the product of mercury concentration and precipitation. The data were analyzed for seasonal patterns, temporal trends, and geographic differences. In the 2 years, 520 weekly samples were collected at the 5 monitoring stations and 448 of these samples had sufficient precipitation to compute mercury wet deposition. The 2-year mean mercury concentration at the five monitoring stations (normalized to the sample volume) was 10.6 ng/L. As a reference for comparison, the total mercury concentration in 41 percent of the samples analyzed was greater than the statewide Indiana water-quality standard for mercury (12 ng/L, protecting aquatic life) and 99 percent of the concentrations exceeded the most conservative Indiana water-quality criterion (1.3 ng/L, protecting wild mammals and birds). The normalized annual mercury concentration at Clifty Falls in 2004 was the fourth highest in the NADP–MDN in eastern North America that year. In 2005, the mercury concentrations at Clifty Falls and Indiana Dunes were the ninth highest in the NADP–MDN in eastern North America. At the five monitoring stations during the study period, the mean weekly total mercury deposition was 0.208 µg/m 2 (micrograms per square meter) and mean annual total mercury deposition was 10.8 µg/m 2 . The annual mercury deposition at Clifty Falls in 2004 and 2005 was in the top 25 percent of the NADP–MDN stations in eastern North America. Mercury concentrations and deposition varied at the five monitoring stations during 2004–2005. Mercury concentrations in wet-deposition samples ranged from 1.2 to 116.6 ng/L and weekly mercury deposition ranged from 0.002 to 1.74 µg/m 2 . Data from weekly samples exhibited seasonal patterns. During April through September, total mercury concentrations and deposition were higher than the median for all samples. Annual precipitation at four of the five monitoring stations was within 10 percent of normal both years, with the exception of Indiana Dunes, where precipitation was 23 percent below normal in 2005. Episodes of high mercury deposition, which were the top 10 percent of weekly mercury deposition at the five monitoring stations, contributed 39 percent of all mercury deposition during 2004–2005. Mercury deposition more than 1.04 µg/m 2 (5 times the mean weekly deposition) was recorded for 12 samples. These episodes of highest mercury deposition were recorded at all five monitoring stations, but the most (7 of 12) were at Clifty Falls and contributed 34.4 percent of the total deposition at that station during 2004–2005. Weekly samples with high mercury deposition may help to explain the differences in annual mercury deposition among the five monitoring stations in Indiana. A statistical evaluation of the monitoring data for 2001–2005 indicated several statistically significant temporal trends. A statewide (5-station) decrease (p = 0.007) in mercury deposition and a statewide decrease (p = 0.059) in mercury concentration were shown. Decreases in mercury deposition (p = 0.061 and p = 0.083) were observed at Roush Lake and Bloomington. A statistically significant trend was not observed for precipitation at the five monitoring stations during this 5-year period. A potential explanation for part of the statewide decrease in mercury concentration and mercury deposition was a 28 percent decrease in the total estimated annual mercury emissions in Indiana between 2002 and 2005. Mercury deposition statistically was correlated most closely to precipitation in the 448 samples, 2004–2005, and this relation was demonstrated by statewide maps of annual precipitation and annual mercury deposition based on precipitation data from 127 National Weather Service Cooperative Observer Program stations. However, one area in southeastern Indiana in the vicinity of Clifty Falls exhibited high mercury deposition that might be related more to mercury concentration than to precipitation. This is because areas with the same range of precipitation as southeastern Indiana were mapped with less mercury deposition. Other data demonstrate a geographic difference for mercury in precipitation in the vicinity of the Clifty Falls monitoring station. The weekly mercury concentrations at Clifty Falls were statistically higher than concentrations at Roush Lake, Fort Harrison, and Bloomington. Clifty Falls data ranked highest among the five monitoring stations for mercury concentration and mercury deposition, 2004–2005, and in the previous 3 years. Episodes of high mercury deposition were recorded most often at Clifty Falls in 2004–2005 and in the previous 3 years. Statistical trends in mercury concentration or mercury deposition were not observed for the Clifty Falls data. A potential explanation for this geographic difference is that annual mercury emissions from sources in the vicinity of Clifty Falls were higher than those at the other stations. Other factors may help explain the differences in total mercury concentrations, such as the types of mercury emissions, mercury transport from stationary sources outside Indiana, and meteorological conditions. Additional data are needed to assign a localized or regional boundary to the area affected by high deposition of mercury near Clifty Falls.

Indiana↗

USGS2016

Groundwater-flow model for the Wood River Valley aquifer system, south-central Idaho

A three-dimensional numerical model of groundwater flow was developed for the Wood River Valley (WRV) aquifer system, Idaho, to evaluate groundwater and surface-water availability at the regional scale. This mountain valley is located in Blaine County and has a drainage area of about 2,300 square kilometers (888 square miles). The model described in this report can serve as a tool for water-rights administration and water-resource management and planning. The model was completed with support from the Idaho Department of Water Resources, and is part of an ongoing U.S. Geological Survey effort to characterize the groundwater resources of the WRV. A highly reproducible approach was taken for constructing the WRV groundwater-flow model. The collection of datasets, source code, and processing instructions used to construct and analyze the model was distributed as an R statistical-computing and graphics package. Flow in the WRV aquifer was simulated using the MODFLOW-USG groundwater flow model. The transient flow model simulates groundwater flow between 1995 and 2010. The model uses a 100-meter (328-feet) uniform grid spacing with 54,922 active model cells distributed over three model layers. A confining unit in the south-central part of the Bellevue fan necessitated the use of a multi-layer model. Specified-flow boundaries were used to simulate the groundwater inflows from each of the major tributary basins (also known as tributary basin underflow) and the areal recharge of precipitation and applied irrigation. Head‑dependent flow boundaries were used to simulate the stream-aquifer flow exchange in river reaches and the groundwater discharge at the outlet boundaries of Stanton Crossing and Silver Creek. The model was calibrated by adjusting aquifer hydraulic properties to match simulated and measured water levels and stream-aquifer flow exchange, using the parameter-estimation program PEST. The model reasonably simulated the measured water-table elevation, orientation, and gradients. Stream-aquifer flow exchange along river reaches also was reasonably simulated by the model. Inflow into the WRV aquifer system originates from three sources (from largest to smallest): Streamflow loss from the Big Wood River and Silver Creek, Areal recharge of precipitation and applied irrigation, and Tributary basin underflow. Outflow from the WRV aquifer system originates from five sources (from largest to smallest): Aquifer discharge into the Big Wood River and Silver Creek, Production-well pumping, Discharge across the water table into the vadose zone, Subsurface outflow beneath Silver Creek near Picabo, and Subsurface outflow beneath the Big Wood River near Stanton Crossing. Temporal changes in aquifer storage are most affected by areal recharge and groundwater pumping, and also contribute to changes in streamflow gains.

Idaho↗

USGS2012

Assessing the vulnerability of public-supply wells to contamination: Rio Grande aquifer system in Albuquerque, New Mexico

This fact sheet highlights findings from the vulnerability study of a public-supply well in Albuquerque, New Mexico (hereafter referred to as “the study well”). The study well produces about 3,000 gallons of water per minute from the Rio Grande aquifer system. Water samples were collected at the study well, at two other nearby public-supply wells, and at monitoring wells installed in or near the simulated zone of contribution to the study well. Untreated water samples from the study well contained arsenic at concentrations exceeding the Maximum Contaminant Level (MCL) of 10 micrograms per liter (µg/L) established by the U.S. Environmental Protection Agency for drinking water. Volatile organic compounds (VOCs) and nitrate also were detected, although at concentrations at least an order of magnitude less than established drinking-water standards, where such standards exist. Overall, study findings point to four primary influences on the movement and (or) fate of contaminants and the vulnerability of the public-supply well in Albuquerque: (1) groundwater age (how long ago water entered, or recharged, the aquifer), (2) groundwater development (introduction of manmade recharge and discharge sources), (3) natural geochemical conditions of the aquifer, and (4) seasonal pumping stresses. Concentrations of the isotope carbon-14 indicate that groundwater from most sampled wells in the local study area is predominantly water that entered, or recharged, the aquifer more than 6,000 years ago. However, the additional presence of the age tracer tritium in several groundwater samples at concentrations above 0.3 tritium units indicates that young (post-1950) recharge is reaching the aquifer across broad areas beneath Albuquerque. This young recharge is mixing with the thousands-of-years-old water, is migrating to depths as great as 245 feet below the water table, and is traveling to some (but not all) of the public-supply wells sampled. Most groundwater samples containing a fraction of young water also contain manmade VOCs, including chloroform (a byproduct of drinking-water chlorination), which indicates that the source of young recharge is, at least in part, infiltration of chlorinated municipal-supply water from leaking waterlines and sewerlines or from turf watering. Other likely manmade, urban recharge sources are seepage from constructed ponds and unlined portions of a stormwater diversion channel. A regional-scale computer-model simulation of groundwater flow and transport to the public-supply well shows that manmade sources of recharge and discharge that were added after about 1930 have greatly altered directions of groundwater flow near Albuquerque and have caused water levels to decline by as much as 120 feet. Local-scale simulations show that seasonal changes in the pumping schedule of the study well affect the age and quality of water produced by the well. Increased pumping during the summer causes significant volumes of water to flow downward from the shallow to the intermediate zones of the aquifer, causing a higher fraction of young water to be produced by the well in the summer than in the winter months and a corresponding increase in VOC detections in the summer relative to the winter. During the winter when the study-well pump is idle for several hours each day, old, high-arsenic water from the deep zone of the aquifer travels up the wellbore and exits into the intermediate zone of the aquifer. When the pump is activated in the winter (for a relatively short time each day), some of the leaked, high-arsenic water is recaptured by the well. This results in a higher arsenic concentration (commonly more than 12 µg/L) in water produced in the winter than in the summer, and a smaller fraction of young water being produced by the well in the winter than in the summer (6 percent in the winter, compared to 11 percent in the summer). Knowledge of the vertical flow direction (both natural and pumping-enhanced) in the vicinity of a long-screened well, coupled with understanding of variations in contaminant concentrations with depth in the aquifer, can help water managers predict the positive or negative effect that wellbore flow will have on water quality and can lead to development of strategies to mitigate contamination (such as changes in pumping schedules or development of devices to inhibit wellbore flow when the pump is off).

New Mexico↗

USGS2013

Simulations of groundwater flow, transport, and age in Albuquerque, New Mexico, for a study of transport of anthropogenic and natural contaminants (TANC) to public-supply wells

Vulnerability to contamination from manmade and natural sources can be characterized by the groundwater-age distribution measured in a supply well and the associated implications for the source depths of the withdrawn water. Coupled groundwater flow and transport models were developed to simulate the transport of the geochemical age-tracers carbon-14, tritium, and three chlorofluorocarbon species to public-supply wells in Albuquerque, New Mexico. A separate, regional-scale simulation of transport of carbon-14 that used the flow-field computed by a previously documented regional groundwater flow model was calibrated and used to specify the initial concentrations of carbon-14 in the local-scale transport model. Observations of the concentrations of each of the five chemical species, in addition to water-level observations and measurements of intra-borehole flow within a public-supply well, were used to calibrate parameters of the local-scale groundwater flow and transport models. The calibrated groundwater flow model simulates the mixing of “young” groundwater, which entered the groundwater flow system after 1950 as recharge at the water table, with older resident groundwater that is more likely associated with natural contaminants. Complexity of the aquifer system in the zone of transport between the water table and public-supply well screens was simulated with a geostatistically generated stratigraphic realization based upon observed lithologic transitions at borehole control locations. Because effective porosity was simulated as spatially uniform, the simulated age tracers are more efficiently transported through the portions of the simulated aquifer with relatively higher simulated hydraulic conductivity. Non-pumping groundwater wells with long screens that connect aquifer intervals having different hydraulic heads can provide alternate pathways for contaminant transport that are faster than the advective transport through the aquifer material. Simulation of flow and transport through these wells requires time discretization that adequately represents periods of pumping and non-pumping. The effects of intra-borehole flow are not fully represented in the simulation because it employs seasonal stress periods, which are longer than periods of pumping and non-pumping. Further simulations utilizing daily pumpage data and model stress periods may help quantify the relative effects of intra-borehole versus advective aquifer flow on the transport of contaminants near the public-supply wells. The fraction of young water withdrawn from the studied supply well varies with simulated pumping rates due to changes in the relative contributions to flow from different aquifer intervals. The advective transport of dissolved solutes from a known contaminant source to the public-supply wells was simulated by using particle-tracking. Because of the transient groundwater flow field, scenarios with alternative contaminant release times result in different simulated-particle fates, most of which are withdrawn from the aquifer at wells that are between the source and the studied supply well. The relatively small effective porosity required to simulate advective transport from the simulated contaminant source to the studied supply well is representative of a preferential pathway and not the predominant aquifer effective porosity that was estimated by the calibration of the model to observed chemical-tracer concentrations.

New Mexico↗

USGS2009

Simulation of the Regional Ground-Water-Flow System and Ground-Water/Surface-Water Interaction in the Rock River Basin, Wisconsin

A regional, two-dimensional, areal ground-water-flow model was developed to simulate the ground-water-flow system and ground-water/surface-water interaction in the Rock River Basin. The model was developed by the U.S. Geological Survey (USGS), in cooperation with the Rock River Coalition. The objectives of the regional model were to improve understanding of the ground-water-flow system and to develop a tool suitable for evaluating the effects of potential regional water-management programs. The computer code GFLOW was used because of the ease with which the model can simulate ground-water/surface-water interactions, provide a framework for simulating regional ground-water-flow systems, and be refined in a stepwise fashion to incorporate new data and simulate ground-water-flow patterns at multiple scales. The ground-water-flow model described in this report simulates the major hydrogeologic features of the modeled area, including bedrock and surficial aquifers, ground-water/surface-water interactions, and ground-water withdrawals from high-capacity wells. The steady-state model treats the ground-water-flow system as a single layer with hydraulic conductivity and base elevation zones that reflect the distribution of lithologic groups above the Precambrian bedrock and a regionally significant confining unit, the Maquoketa Formation. In the eastern part of the Basin where the shale-rich Maquoketa Formation is present, deep ground-water flow in the sandstone aquifer below the Maquoketa Formation was not simulated directly, but flow into this aquifer was incorporated into the GFLOW model from previous work in southeastern Wisconsin. Recharge was constrained primarily by stream base-flow estimates and was applied uniformly within zones guided by regional infiltration estimates for soils. The model includes average ground-water withdrawals from 1997 to 2006 for municipal wells and from 1997 to 2005 for high-capacity irrigation, industrial, and commercial wells. In addition, the model routes tributary base flow through the river network to the Rock River. The parameter-estimation code PEST was linked to the GFLOW model to select the combination of parameter values best able to match more than 8,000 water-level measurements and base-flow estimates at 9 streamgages. Results from the calibrated GFLOW model show simulated (1) ground-water-flow directions, (2) ground-water/surface-water interactions, as depicted in a map of gaining and losing river and lake sections, (3) ground-water contributing areas for selected tributary rivers, and (4) areas of relatively local ground water captured by rivers. Ground-water flow patterns are controlled primarily by river geometries, with most river sections gaining water from the ground-water-flow system; losing sections are most common on the downgradient shore of lakes and reservoirs or near major pumping centers. Ground-water contributing areas to tributary rivers generally coincide with surface watersheds; however the locations of ground-water divides are controlled by the water table, whereas surface-water divides are controlled by surface topography. Finally, areas of relatively local ground water captured by rivers generally extend upgradient from rivers but are modified by the regional flow pattern, such that these areas tend to shift toward regional ground-water divides for relatively small rivers. It is important to recognize the limitations of this regional-scale model. Heterogeneities in subsurface properties and in recharge rates are considered only at a very broad scale (miles to tens of miles). No account is taken of vertical variations in properties or pumping rates, and no provision is made to account for stacked ground-water-flow systems that have different flow patterns at different depths. Small-scale flow systems (hundreds to thousands of feet) associated with minor water bodies are not considered; as a result, the model is not currently designed for simulating site-specifi

Scientific Investigations Report↗

USGS2012

Simulation of streamflow, evapotranspiration, and groundwater recharge in the middle Nueces River watershed, south Texas, 1961-2008

The U.S. Geological Survey—in cooperation with the U.S. Army Corps of Engineers, Fort Worth District; City of Corpus Christi; Guadalupe–Blanco River Authority; San Antonio River Authority; and San Antonio Water System— configured, calibrated, and tested a watershed model for a study area consisting of about 7,726 square miles of the middle Nueces River watershed in south Texas. The purpose of the model is to contribute to the understanding of watershed processes and hydrologic conditions in the middle Nueces River watershed. The model simulates streamflow, evapotranspiration, and groundwater recharge by using a numerical representation of physical characteristics of the landscape and meteorological and streamflow data. Model simulations of streamflow, evapotranspiration, and groundwater recharge were performed for various periods of record depending upon available gaged data for input and comparison, starting as early as 1961. Because of the large size of the study area, the middle Nueces River watershed was divided into eight subwatersheds, and separate Hydrological Simulation Program—FORTRAN models were developed for each subwatershed. Simulation of the overall study area involved running simulations in downstream order. Output from the model was summarized by subwatershed, point locations, stream and reservoir reaches, and the Carrizo– Wilcox aquifer outcrop area. Four long-term U.S. Geological Survey streamflow-gaging stations were used for streamflow model calibration and testing with data from 1990 to 2008. Monthly evaporation estimates from 2001 to 2008 and waterlevel data from 1961 to 2008 at Lake Corpus Christi also were used for model calibration. Additionally, evapotranspiration data for 2006–8 from a U.S. Geological Survey meteorological station in Medina County were used for calibration. Streamflow calibrations were considered poor to very good. The 2000–8 calibration results were characterized as good to very good for total flow volumes and for the volume of the highest 10 percent of daily flows. Calibration results for streamflow volumes of the lowest 50 percent of daily flows were considered poor. The daily streamflow calibration at U.S. Geological Survey streamflow-gaging station 08210000 Nueces River near Three Rivers, Tex., had the lowest (best) root mean square error, and U.S. Geological Survey streamflow-gaging station 08194500 Nueces River near Tilden, Tex., had the highest root mean square error expressed as a percentage of the mean flow rate. The mean daily reservoir volume during 1961–2008 was 182,000 acre-feet. Simulated mean daily reservoir volume was within 9 percent of this computed volume. Selected results of the model include streamflow yields for the subwatersheds and water-balance information for the Carrizo–Wilcox aquifer outcrop area. For the entire model domain, the area-weighted mean streamflow yield from 1961 to 2008 was 1.12 inches/year. The mean annual rainfall on the outcrop area during the 1961–2008 simulation period was 21.7 inches. Of this rainfall, an annual mean of 20.1 inches (about 93 percent) was simulated as evapotranspiration, 1.2 inches (about 6 percent) was simulated as groundwater recharge, and 0.5 inches (about 2 percent) was simulated as surface runoff.

Texas↗

USGS2015

Effect of land-applied biosolids on surface-water nutrient yields and groundwater quality in Orange County, North Carolina

Land application of municipal wastewater biosolids is the most common method of biosolids management used in North Carolina and the United States. Biosolids have characteristics that may be beneficial to soil and plants. Land application can take advantage of these beneficial qualities, whereas disposal in landfills or incineration poses no beneficial use of the waste. Some independent studies and laboratory analysis, however, have shown that land-applied biosolids can pose a threat to human health and surface-water and groundwater quality. The effect of municipal biosolids applied to agriculture fields is largely unknown in relation to the delivery of nutrients, bacteria, metals, and contaminants of emerging concern to surface-water and groundwater resources. Therefore, the North Carolina Department of Environment and Natural Resources (NCDENR) collaborated with the U.S. Geological Survey (USGS) through the 319 Nonpoint Source Program to better understand the transport of nutrients and bacteria from biosolids application fields to groundwater and surface water and to provide a scientific basis for evaluating the effectiveness of the current regulations. The USGS conducted a paired agricultural watershed study in the Collins Creek and Cane Creek Reservoir watersheds in Orange County, North Carolina. Field activities were conducted from March 2011 through May 2013 at two field study sites, including biosolids field application sites owned by Orange County Water and Sewer Authority (OWASA) in the Collins Creek watershed and a background study site in the Cane Creek watershed that has no fields receiving biosolids applications. Samples of biosolids source material and soil were collected from the land-application fields for laboratory analyses. Soil samples were also collected from a background agricultural field in the Cane Creek watershed that has never received land-applied municipal biosolids. Shallow groundwater samples were collected quarterly from new monitoring wells installed by NCDENR along the edge of the biosolids land-application fields and a background agricultural field for laboratory analyses. Two surface-water monitoring sites were established on Collins Creek to compute continuous streamflow and collect discrete baseflow and stormwater runoff water-quality data upstream and downstream from the biosolids land-application fields. Surface water-quality samples were also collected for baseflow and stormwater runoff conditions at an existing USGS streamgage on Cane Creek to monitor water-quality conditions in the background study watershed. The study primarily focused on nutrients and bacteria; however, data for field properties and water-quality constituents, including metals, major ions, and contaminants of emerging concern (household-, industrial-, and agricultural-use compounds, pharmaceutical compounds, hormones, and antibiotics) also were collected and used in the analyses. There were no exceedances of the 10 elements with designated U.S. Environmental Protection Agency (EPA) ceiling concentrations for land-applied biosolids in any of the biosolids samples. Treatment processes and storage techniques used by OWASA are effective in eliminating Escherichia coli and fecal coliform bacteria from biosolids. Copper, molybdenum, total Kjeldahl nitrogen, and total phosphorus were elevated in the soil from biosolids land-application fields relative to the background field. The relative richness of these constituents in the biosolids land-application fields is consistent with biosolids being the source of the elevated concentrations given the relatively high concentrations of these constituents in the biosolids samples that were collected. Shallow groundwater in the transitional zone wells, which were located adjacent to and topographically downgradient from all the biosolids land-application fields, were found to be statistically different and had higher nitrate concentrations (medians greater than 12 milligrams per liter) than all the other wells sampled as part of the study. Surface-water nutrient concentrations and yields, primarily nitrate, were higher at the monitoring site on Collins Creek downstream from the biosolids land-application fields than the other study sites that drained watersheds without biosolids land application. The largest differences in concentrations between sites were measured at baseflow conditions, which indicate that the main cause of these differences, particularly between Cane Creek and the Collins Creek site downstream from the OWASA application fields, is related to nitrate contribution from the shallow groundwater. Contaminants of emerging concern were detected in approximately 40 percent of the laboratory analyses of the biosolids samples and more frequently in soil samples from the biosolids land-application fields (approximately 40 percent of laboratory analyses) relative to the soil samples from the background field (approximately 12 percent of laboratory analyses). However, contaminants of emerging concern detected in the laboratory analysis for this study do not appear to be good indicators of human-waste contaminants derived from land-applied biosolids in groundwater or surface-water because the number of detections and concentrations at the background wells and surface-water monitoring sites are similar to or higher than those at wells and monitoring sites adjacent to or downstream from the biosolids land-application fields. The data, analysis, and conclusions associated with this study can be used by regulatory agencies, resource managers, and wastewater-treatment operators to (1) better understand the quantity and characteristics of nutrients, bacteria, metals, and contaminants of emerging concern that are transported away from biosolids land-application fields to surface water and groundwater under current regulations for the purposes of establishing effective total maximum daily loads (TMDLs) and restoring impaired water resources, (2) assess how well existing regulations protect waters of the State and potentially recommend effective changes to regulations or land-application procedures, and (3) establish a framework for developing guidance on effective techniques for monitoring and regulatory enforcement of permitted biosolids land-application fields.

North Carolina↗

USGS2022

Distribution of streamflow, sediment, and nutrients entering Galveston Bay from the Trinity River, Texas, 2016–19

The U.S. Geological Survey (USGS), in cooperation with the Texas Water Development Board, collected streamflow and water-quality data at USGS monitoring stations in the lower Trinity River Basin from January 2016 to December 2019 to characterize streamflow, nutrients, and suspended sediment entering Galveston Bay from the Trinity River. Results from previous studies indicate that water from the main channel of the Trinity River is diverted into surrounding wetlands and water bodies and is stored or discharged directly into Galveston Bay through distributary channels in the delta. This study provides an assessment of the distribution of streamflow in the various channels that form the delta of the Trinity River to evaluate the effects of streamflow diversions on the eventual supply of freshwater, nutrients, and suspended sediment to Galveston Bay. Instantaneous streamflow data and continuous streamflow records from USGS monitoring stations in the delta of the Trinity River were used to quantify freshwater inflow into Galveston Bay and assess the distribution of streamflow in the lowermost reaches of the Trinity River Basin. In this report, periods in which releases from Lake Livingston caused a rise in streamflow farther downstream at USGS station 08067000 Trinity River at Liberty, Tex. (hereinafter referred to as the “Liberty site”) that did not exceed 20,000 cubic feet per second (ft 3 /s) are referred to as “low-flow events,” and periods in which streamflow at the Liberty site exceeded 20,000 ft 3 /s are referred to as “high-flow events.” During this study, it was estimated that only about 55 percent of the total water volume released from Lake Livingston was accounted for at USGS station 08067252 Trinity River at Wallisville, Tex. (hereinafter referred to as the “Wallisville site”), which is approximately 8 river miles upstream from where the Trinity River enters Galveston Bay. The difference in water volumes between what is released from Lake Livingston and what is measured at the Wallisville site is consistent with findings from previous studies and indicates that a large part of the volume released from Lake Livingston does not reach Galveston Bay through the main channel of the Trinity River. To assess the distribution of streamflow and estimate the amount of water diverted from the main channel of the Trinity River into distributary channels, instantaneous streamflow measurements were made at USGS station 08067230 Old River Lake near Wallisville, Tex. (hereinafter referred to as the “Old River Lake site”) and the Wallisville site during a range of hydrologic conditions. Results indicate that a large portion of the freshwater inflow was likely delivered to Galveston Bay through pathways other than the main channel of the Trinity River, including Old River Lake. When streamflow at the Liberty site, located upstream from the Wallisville site, exceeded approximately 40,000 ft 3 /s, Old River Lake and its network of hydrologically connected channels likely became the primary pathway for freshwater inflow entering Galveston Bay. Water quality was characterized from discrete samples collected during a range of hydrologic conditions at the Old River Lake site and the Wallisville site in order to evaluate the effects of streamflow diversions on the supply of suspended sediment and nutrients into Galveston Bay. Suspended-sediment concentrations were typically higher at the Wallisville site than at the Old River Lake site, likely because of lower water velocities at the Old River Lake site than at the Wallisville site; low water velocities allow suspended sediment to settle, thus reducing concentrations. Suspended-sediment loads were also typically higher at the Wallisville site than at the Old River Lake site during high-flow events. However, when streamflows at the Liberty site exceeded approximately 60,000 ft 3 /s, suspended-sediment loads were higher at the Old River Lake, which likely became the primary pathway for suspended-sediment delivery into Galveston Bay. Suspended-sediment concentrations and loads were computed at the Wallisville and Liberty sites for the duration of 11 hydrologic events representing different streamflows by using the regression equations developed for each monitoring station. Overall, approximately 25 percent of the total sediment load measured during events at the Liberty site was measured at the Wallisville site, indicating that only a portion of the suspended-sediment load from the Liberty site reached Galveston Bay through the main channel of the Trinity River during the measured events. Based on data from discrete samples, some of this sediment load was diverted into Old River Lake and associated distributary channels. Results from analysis of nutrient samples indicate that streamflow conditions affect the nitrogen concentrations in the delta of the Trinity River. At the Old River Lake site, nitrate plus nitrite and total dissolved nitrogen concentrations were typically lower during low-flow conditions than during high-flow events; low-flow conditions represent low-flow events or tidal-flow conditions (during low-flow conditions the streamflow at the Liberty site was less than 20,000 ft 3 /s). Lower concentrations of nitrate plus nitrite and total dissolved nitrogen at the Old River Lake site may be associated with various physical and biogeochemical processes, including the transformation and biological uptake of nitrate, nitrite, and other species of nitrogen resulting from extended water residence times and relatively small inputs of nitrogen from the upstream reaches of the Trinity River Basin. During high-flow events, the proportions of nitrogen species were similar among sites, indicating that the travel path through wetlands and channels surrounding Old River Lake likely does not affect the relative concentrations of the various nitrogen species present in freshwater inflow to Galveston Bay. Results from analysis of nutrient samples also indicate that the pathways for nutrient delivery from the Trinity River into Galveston Bay are dependent on event magnitude. When streamflows at the Liberty site were low (approximately 20,000 ft 3 /s), the main channel of the Trinity River was the primary pathway for nitrogen and phosphorus entering Galveston Bay. Once streamflow at the Liberty site exceeded 20,000 ft 3 /s, however, the contribution of nutrient loading through Old River Lake to Galveston Bay increased proportionally to the nutrient loading in the main channel, and when streamflow at the Liberty site exceeded approximately 50,000 ft 3 /s, Old River Lake likely became the primary pathway for nutrient delivery into Galveston Bay.

Texas↗

USGS2014

Evaluation of a mass-balance approach to determine consumptive water use in northeastern Illinois

A principal component of evaluating and managing water use is consumptive use. This is the portion of water withdrawn for a particular use, such as residential, which is evaporated, transpired, incorporated into products or crops, consumed by humans or livestock, or otherwise removed from the immediate water environment. The amount of consumptive use may be estimated by a water (mass)-balance approach; however, because of the difficulty of obtaining necessary data, its application typically is restricted to the facility scale. The general governing mass-balance equation is: Consumptive use = Water supplied - Return flows . This study explored a mass-balance field-based computation of consumptive use in a residential setting at the scale of a sanitary sewer service area (sewershed). In addition, the feasibility (cost and difficulty) and relative uncertainties (accuracies) associated with applying the approach at this scale were evaluated. The study was conducted during 2011–13 within a 3.5-square mile (mi 2 ) sewershed confined to a predominantly residential area of Elk Grove Village, Illinois. Following background evaluation of the geohydrologic setting, sewershed infrastructure, and possible components of supplied and returned water, the identified primary components were 1. public water deliveries by the Elk Grove Village Department of Public Works, 2. self-served groundwater withdrawals in an included unincorporated neighborhood with public sanitary sewer service, 3. return flows to the sanitary sewer system, and 4. direct return of water discharged from swimming pools to Salt Creek. Water volumes principally were reported for deliveries, measured for sanitary sewer returns by using an acoustic Doppler current-velocity meter, and estimated for domestic withdrawals and swimming pool discharges to storm sewers. All water volumes required some degree of estimation. Observation wells were installed adjacent to sewer pipelines (lines) to determine the depth of the water table relative to that of the sewer lines and to collect water samples for detection of optical brighteners, as they are routinely discharged as clotheswashing waste to sanitary sewers. These data provided qualitative information on gains (inflow and infiltration) and losses (exfiltration) of sewer flow by pipe leakage, which might otherwise not be considered in the sewer flow return measurements. Hydrographs of sewer flow also were evaluated to identify and estimate storm-associated inputs to sewer flow. The volume of sanitary sewer return flow (778 million gallons per year [Mgal/yr]) was determined to substantially exceed the volume of supplied water (566 Mgal/yr), thus, for this study setting, voiding the utility of the applied mass-balance approach for estimating consumptive water use. Mass-balance components, including sanitary sewer flow and supplied-water use, were estimated within reasonable limits of uncertainty. Evidence of a water table that is typically shallower than the area’s sewer lines, yet is sometimes depressed near more deeply buried sewer lines, suggests groundwater infiltration into the sewers contributes to the excess volume of return flow. Technical obstacles and project resources precluded accurate quantification of infiltration volumes and other gains and losses to sanitary sewer flow. As estimated from various simplified methods, a minimum of 26 percent of return flow measured in the sanitary sewer represented groundwater infiltration and stormwater inflow; separately, about 2 percent of return flow was estimated as inflow. On the basis of the alternative winter base-rate method, consumptive use in the sewershed was estimated as 13 percent, which compares favorably with that used by the State of Illinois for Lake Michigan allocation accounting (10 percent) and other States and Canadian Provinces in the Great Lakes region (generally 10-15 percent). The study also provided other findings considered useful to studies of water use and to performance evaluation of sanitary sewer infrastructure. In urban residential settings, the comparatively small volumes of nonpublic sources of water (self-supplied) and direct (nonsanitary) return flow potentially can be ignored in the estimation of consumptive use. An acoustic Doppler current-velocity meter can be used in sanitary sewers to accurately measure discharge and reasonably estimate storm-associated inflows. Hourly to daily patterns of water use can be readily identified and quantified in the return flow record for the sanitary sewers. Relative volumes of infiltration gains (and exfiltration losses) can be substantial, even in sewer systems of communities making significant investments in system upgrades to limit sewer line leakage. Monitoring of optical brighteners in groundwater (and potentially in sanitary sewer flow) can provide a useful means of identifying probable leakage from (and to) sewer lines. Accurate quantification of gains and losses to sanitary sewer flow at the sewershed scale will require additional research effort and technical advances. Under ideal conditions, accurate quantification of consumptive use at the sewershed scale by the described mass-balance approach might be possible. Under most prevailing conditions, quantification likely would be more costly and time consuming than that of the present study, given the freely contributed technical support of the host community and relatively appropriate conditions of the study area. Essentials to quantification of consumptive use are a fully cooperative community, storm and sanitary sewers that are separate, and newer sewer infrastructure and (or) a robust program for limiting infiltration, exfiltration, and inflow.

Illinois↗

USGS2016

Generalized sediment budgets of the Lower Missouri River, 1968–2014

Sediment budgets of the Lower Missouri River were developed in a study led by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers. The scope of the study included the development of a long-term (post-impoundment, 1968–2014) average annual sediment budget and selected annual, monthly, and daily sediment budgets for a reach and period that adequate data were available. Included in the analyses were 31 main-stem and tributary stations of the Lower Missouri River and two Mississippi River stations—the Mississippi River below Grafton, Illinois, and the Mississippi River at St. Louis, Missouri. Long-term average annual suspended-sediment loads of Missouri River main-stem stations ranged from 0.33 million tons at the Missouri River at Yankton, South Dakota, station to 71.2 million tons at Missouri River at Hermann, Mo., station. Gaged tributary gains accounted for 9–36 percent of the local reach budgets and cumulative gaged tributary contributions accounted for 84 percent of the long-term average suspended-sediment load of the Missouri River at Hermann, Mo., station. Although the sediment budgets for seven defined main-stem reaches generally were incomplete—missing bedload, reach storage, and ungaged tributary contributions—the budget residuals (net result of sediment inputs and outputs) for six of the seven reaches ranged from -7.0 to 1.7 million tons, or from -9.2 to 4.0 percent of the reach output suspended-sediment load, and were within the 10 percent reported measurement error of annual suspended-sediment loads for large rivers. The remaining reach, downstream from Gavin’s Point Dam, extended from Yankton, S. Dak., to Sioux City, Iowa, and had a budget residual of -9.8 million tons, which was -88 percent of the suspended-sediment load at Sioux City. The Lower Missouri River reach from Omaha, Nebraska, to Nebraska City, Nebr., had periods of concurrent sediment data for each primary budget component with which to analyze and determine a suspended-sediment budget for selected annual, monthly, and daily time increments. The temporal changes in the cumulative annual budget residuals were poorly correlated with the comparatively steady 1968–2011 annual stage trends at the Missouri River at Nebraska City, Nebr., station. An accurate total sediment budget is developed by having concurrent data available for all primary suspended and bedload components for a reach of interest throughout a period. Such a complete budget, with concurrent record for suspended-sediment load and bedload components, is unavailable for any reach and period in the Lower Missouri River. The primary data gaps are in bedload data, and also in suspended-sediment gains and losses including ungaged tributary inputs and sediment storage. Bedload data gaps in the Missouri River Basin are much more prevalent than suspended-sediment data gaps, and the first step in the development of reach bedload budgets is the establishment of a standardized bedload monitoring program at main-stem stations. The temporal changes in flow-adjusted suspended-sediment concentrations analyzed at main-stem Missouri River stations indicated an overall downward change in concentrations between 1968 and 2014. Temporary declines in flow-adjusted suspended-sediment concentrations during and following large floods were evident but generally returned to near pre-flood values within about 6 months. Data uncertainties associated with the development of a sediment budget include uncertainties associated with the collection of suspended-sediment and bedload data and the computation of suspended-sediment loads. These uncertainties vary depending on the frequency of data collection, the variability of conditions being represented by the discrete samples, and the statistical approach to suspended-sediment load computations. The coefficients of variation of suspended-sediment loads of Missouri River tributary stations for 1968–2014 were greater, 75.0 percent, than the main-stem stations, 47.1 percent. The lower coefficient of variation at main-stem stations compared to tributaries, primarily is the result of the lower variability in streamflow and sediment discharge identified at main-stem stations. To obtain similar accuracy between suspended-sediment loads at main-stem and tributary stations, a longer period of record is required of the tributary stations. During 1968–2014, however, the Missouri River main-stem station record was much more complete (87 percent) than the tributary station record (28 percent).

Missouri River Basin↗

USGS2010

Approaches to highly parameterized inversion: A guide to using PEST for model-parameter and predictive-uncertainty analysis

Analysis of the uncertainty associated with parameters used by a numerical model, and with predictions that depend on those parameters, is fundamental to the use of modeling in support of decisionmaking. Unfortunately, predictive uncertainty analysis with regard to models can be very computationally demanding, due in part to complex constraints on parameters that arise from expert knowledge of system properties on the one hand (knowledge constraints) and from the necessity for the model parameters to assume values that allow the model to reproduce historical system behavior on the other hand (calibration constraints). Enforcement of knowledge and calibration constraints on parameters used by a model does not eliminate the uncertainty in those parameters. In fact, in many cases, enforcement of calibration constraints simply reduces the uncertainties associated with a number of broad-scale combinations of model parameters that collectively describe spatially averaged system properties. The uncertainties associated with other combinations of parameters, especially those that pertain to small-scale parameter heterogeneity, may not be reduced through the calibration process. To the extent that a prediction depends on system-property detail, its postcalibration variability may be reduced very little, if at all, by applying calibration constraints; knowledge constraints remain the only limits on the variability of predictions that depend on such detail. Regrettably, in many common modeling applications, these constraints are weak. Though the PEST software suite was initially developed as a tool for model calibration, recent developments have focused on the evaluation of model-parameter and predictive uncertainty. As a complement to functionality that it provides for highly parameterized inversion (calibration) by means of formal mathematical regularization techniques, the PEST suite provides utilities for linear and nonlinear error-variance and uncertainty analysis in these highly parameterized modeling contexts. Availability of these utilities is particularly important because, in many cases, a significant proportion of the uncertainty associated with model parameters-and the predictions that depend on them-arises from differences between the complex properties of the real world and the simplified representation of those properties that is expressed by the calibrated model. This report is intended to guide intermediate to advanced modelers in the use of capabilities available with the PEST suite of programs for evaluating model predictive error and uncertainty. A brief theoretical background is presented on sources of parameter and predictive uncertainty and on the means for evaluating this uncertainty. Applications of PEST tools are then discussed for overdetermined and underdetermined problems, both linear and nonlinear. PEST tools for calculating contributions to model predictive uncertainty, as well as optimization of data acquisition for reducing parameter and predictive uncertainty, are presented. The appendixes list the relevant PEST variables, files, and utilities required for the analyses described in the document.

Michigan↗

USGS2017

Analysis of the variability in ground-motion synthesis and inversion

In almost all past inversions of large-earthquake ground motions for rupture behavior, the goal of the inversion is to find the “best fitting” rupture model that predicts ground motions which optimize some function of the difference between predicted and observed ground motions. This type of inversion was pioneered in the linear-inverse sense by Olson and Apsel (1982), who minimized the square of the difference between observed and simulated motions (“least squares”) while simultaneously minimizing the rupture-model norm (by setting the null-space component of the rupture model to zero), and has been extended in many ways, one of which is the use of nonlinear inversion schemes such as simulated annealing algorithms that optimize some other misfit function. For example, the simulated annealing algorithm of Piatanesi and others (2007) finds the rupture model that minimizes a “cost” function which combines a least-squares and a waveform-correlation measure of misfit. All such inversions that look for a unique “best” model have at least three problems. (1) They have removed the null-space component of the rupture model—that is, an infinite family of rupture models that all fit the data equally well have been narrowed down to a single model. Some property of interest in the rupture model might have been discarded in this winnowing process. (2) Smoothing constraints are commonly used to yield a unique “best” model, in which case spatially rough rupture models will have been discarded, even if they provide a good fit to the data. (3) No estimate of confidence in the resulting rupture models can be given because the effects of unknown errors in the Green’s functions (“theory errors”) have not been assessed. In inversion for rupture behavior, these theory errors are generally larger than the data errors caused by ground noise and instrumental limitations, and so overfitting of the data is probably ubiquitous for such inversions. Recently, attention has turned to the inclusion of theory errors in the inversion process. Yagi and Fukahata (2011) made an important contribution by presenting a method to estimate the uncertainties in predicted large-earthquake ground motions due to uncertainties in the Green’s functions. Here we derive their result and compare it with the results of other recent studies that look at theory errors in a Bayesian inversion context particularly those by Bodin and others (2012), Duputel and others (2012), Dettmer and others (2014), and Minson and others (2014). Notably, in all these studies, the estimates of theory error were obtained from theoretical considerations alone; none of the investigators actually measured Green’s function errors. Large earthquakes typically have aftershocks, which, if their rupture surfaces are physically small enough, can be considered point evaluations of the real Green’s functions of the Earth. Here we simulate smallaftershock ground motions with (erroneous) theoretical Green’s functions. Taking differences between aftershock ground motions and simulated motions to be the “theory error,” we derive a statistical model of the sources of discrepancies between the theoretical and real Green’s functions. We use this model with an extended frequency-domain version of the time-domain theory of Yagi and Fukahata (2011) to determine the expected variance 2 τ caused by Green’s function error in ground motions from a larger (nonpoint) earthquake that we seek to model. We also differ from the above-mentioned Bayesian inversions in our handling of the nonuniqueness problem of seismic inversion. We follow the philosophy of Segall and Du (1993), who, instead of looking for a best-fitting model, looked for slip models that answered specific questions about the earthquakes they studied. In their Bayesian inversions, they inductively derived a posterior probability-density function (PDF) for every model parameter. We instead seek to find two extremal rupture models whose ground motions fit the data within the error bounds given by 2 τ , as quantified by using a chi-squared test described below. So, we can ask questions such as, “What are the rupture models with the highest and lowest average rupture speed consistent with the theory errors?” Having found those models, we can then say with confidence that the true rupture speed is somewhere between those values. Although the Bayesian approach gives a complete solution to the inverse problem, it is computationally demanding: Minson and others (2014) needed 1010 forward kinematic simulations to derive their posterior probability distribution. In our approach, only about107 simulations are needed. Moreover, in practical application, only a small set of rupture models may be needed to answer the relevant questions—for example, determining the maximum likelihood solution (achievable through standard inversion techniques) and the two rupture models bounding some property of interest. The specific property that we wish to investigate is the correlation between various rupturemodel parameters, such as peak slip velocity and rupture velocity, in models of real earthquakes. In some simulations of ground motions for hypothetical large earthquakes, such as those by Aagaard and others (2010) and the Southern California Earthquake Center Broadband Simulation Platform (Graves and Pitarka, 2015), rupture speed is assumed to correlate locally with peak slip, although there is evidence that rupture speed should correlate better with peak slip speed, owing to its dependence on local stress drop. We may be able to determine ways to modify Piatanesi and others’s (2007) inversion’s “cost” function to find rupture models with either high or low degrees of correlation between pairs of rupture parameters. We propose a cost function designed to find these two extremal models.

Open-File Report↗