SEARCH · Geology Reports

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 433 records · Page 24Linked to original sources

USGS1997

Characterization of stormwater runoff from the Naval Air Station and Naval Wepons Industrial Reserve Plant, Dallas, Texas, 1994-96

The characterization of stormwater runoff from the Naval Air Station (NAS) and the Naval Weapons Industrial Reserve Plant (NWIRP), Dallas, Texas, is necessary to determine if runoff from the facilities is contributing to off-site contamination of surface waters, A network of five fixed sites and four grab sites was established to collect stormwater-runoff samples from a substantial part of the drainage area of each facility. Fixed sites were instrumented to measure and store precipitation, stage, discharge, and runoff-volume data and to collect flow-weighted composite samples during a storm. Grab and composite samples were collected for six storms at each of the five fixed sites from October 1994 to March 1996. The grab samples were analyzed for about 100 properties and constituents including specific conductance, pH, water temperature, bacteria, trace elements, oil and grease, total phenols, and volatile organic compounds. The composite samples were analyzed for about 220 properties and constituents including specific conductance, pH, chemical oxygen demand, biochemical oxygen demand, major ions, suspended and dissolved solids, nutrients, trace elements, total organic carbon, volatile organic compounds, semivolatile organic compounds, and organochlorine and organophosphorus pesticides. Grab samples were collected for two storms (September 18,1995, and October 2,1995) at each of the four grab sites. The grab samples were analyzed for about 80 constituents including specific conductance, pH, water temperature, trace elements, and volatile organic compounds. Composite samples were collected for two of the six storms sampled at the fixed sites and analyzed for aquatic toxicity. Fathead minnow growth and survival toxicity tests and water flea reproduction and survival toxicity tests were done. Median event-mean concentrations computed for 12 selected constituents in samples from NAS and NWIRP fixed sites were compared to median event-mean concentrations for residential, commercial, industrial, and highway land uses within the Dallas-Fort Worth area computed from data collected for the National Pollutant Discharge Elimination System program. NAS and NWIRP median event-mean concentrations also were compared to those for residential and commercial land uses from the Nationwide Urban Runoff Program.

Open-File Report↗

USGS2005

Volcano surveillance by ACR silver fox

Recent growth in the business of unmanned air vehicles (UAVs) both in the US and abroad has improved their overall capability, resulting in a reduction in cost, greater reliability and adoption into areas where they had previously not been considered. Uses in coastal and border patrol, forestry and agriculture have recently been evaluated in an effort to expand the observed area and reduce surveillance and reconnaissance costs for information gathering. The scientific community has both contributed and benefited greatly in this development. A larger suite of light-weight miniaturized sensors now exists for a range of applications which in turn has led to an increase in the gathering of information from these autonomous vehicles. In October 2004 the first eruption of Mount St Helens since 1986 caused tremendous interest amoUg people worldwide. Volcanologists at the U.S. Geological Survey rapidly ramped up the level of monitoring using a variety of ground-based sensors deployed in the crater and on the flanks of the volcano using manned helicopters. In order to develop additional unmanned sensing methods that can be used in potentially hazardous and low visibility conditions, a UAV experiment was conducted during the ongoing eruption early in November. The Silver Fox UAV was flown over and inside the crater to perform routine observation and data gathering, thereby demonstrating a technology that could reduce physical risk to scientists and other field operatives. It was demonstrated that UAVs can be flown autonomously at an active volcano and can deliver real time data to a remote location. Although still relatively limited in extent, these initial flights provided information on volcanic activity and thermal conditions within the crater and at the new (2004) lava dome. The flights demonstrated that readily available visual and infrared video sensors mounted in a small and relatively low-cost aerial platform can provide useful data on volcanic phenomena. This was made possible by utilizing GPS and computer-controlled flight direction and stabilization to acquire and track target areas within the Mount St. Helens crater. It was also determined that additional light-weight sensor development will be needed to enable autonomous measurements of volcanic gasses and imaging in poor-weather conditions. Copyright ?? 2005 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Conference Paper↗

USGS2005

Development of a comprehensive watershed model applied to study stream yield under drought conditions

We developed a model code to simulate a watershed's hydrology and the hydraulic response of an interconnected stream-aquifer system, and applied the model code to the Lower Republican River Basin in Kansas. The model code links two well-known computer programs: MODFLOW (modular 3-D flow model), which simulates ground water flow and stream-aquifer interaction; and SWAT (soil water assessment tool), a soil water budget simulator for an agricultural watershed. SWAT represents a basin as a collection of subbasins in terms of soil, land use, and weather data, and simulates each subbasin on a daily basis to determine runoff, percolation, evaporation, irrigation, pond seepage, and crop growth. Because SWAT applies a lumped hydrologic model to each sub-basin, spatial heterogeneities with respect to factors such as soil type and land use are not resolved geographically, but can instead be represented statistically. For the Republican River Basin model, each combination of six soil types and three land uses, referred to as a hydrologic response unit (HRU), was simulated with a separate execution of SWAT. A spatially weighted average was then taken over these results for each hydrologic flux and time step by a separate program, SWBAVG. We wrote a package for MODFLOW to associate each subbasin with a subset of aquifer grid cells and stream reaches, and to distribute the hydrologic fluxes given for each subbasin by SWAT and SWBAVG over MODFLOW's stream-aquifer grid to represent tributary flow, surface and ground water diversions, ground water recharge, and evapotranspiration from ground water. The Lower Republican River Basin model was calibrated with respect to measured ground water levels, streamflow, and reported irrigation water use. The model was used to examine the relative contributions of stream yield components and the impact on stream yield and base flow of administrative measures to restrict irrigation water use during droughts. Model results indicate that tributary flow is the dominant component of stream yield and that reduction of irrigation water use produces a corresponding increase in base flow and stream yield. However, the increase in stream yield resulting from reduced water use does not appear to be of sufficient magnitude to restore minimum desirable streamflows.

Groundwater↗

USGS2009

Comparison of groundwater flow in Southern California coastal aquifers

Development of the coastal aquifer systems of Southern California has resulted in overdraft, changes in streamflow, seawater intrusion, land subsidence, increased vertical flow between aquifers, and a redirection of regional flow toward pumping centers. These water-management challenges can be more effectively addressed by incorporating new understanding of the geologic, hydrologic, and geochemical setting of these aquifers. Groundwater and surface-water flow are controlled, in part, by the geologic setting. The physiographic province and related tectonic fabric control the relation between the direction of geomorphic features and the flow of water. Geologic structures such as faults and folding also control the direction of flow and connectivity of groundwater flow. The layering of sediments and their structural association can also influence pathways of groundwater flow and seawater intrusion. Submarine canyons control the shortest potential flow paths that can result in seawater intrusion. The location and extent of offshore outcrops can also affect the flow of groundwater and the potential for seawater intrusion and land subsidence in coastal aquifer systems. As coastal aquifer systems are developed, the source and movement of ground-water and surface-water resources change. In particular, groundwater flow is affected by the relative contributions of different types of inflows and outflows, such as pump-age from multi-aquifer wells within basal or upper coarse-grained units, streamflow infiltration, and artificial recharge. These natural and anthropogenic inflows and outflows represent the supply and demand components of the water budgets of ground-water within coastal watersheds. They are all significantly controlled by climate variability related to major climate cycles, such as the El Niño–Southern Oscillation and the Pacific Decadal Oscillation. The combination of natural forcings and anthropogenic stresses redirects the flow of groundwater and either mitigates or exacerbates the potential adverse effects of resource development, such as declining water levels, sea-water intrusion, land subsidence, and mixing of different waters. Streamflow also has been affected by development of coastal aquifer systems and related conjunctive use. Saline water is the largest water-quality problem in Southern California coastal aquifer systems. Seawater intrusion is a significant source of saline water, but saline water is also known to come from other sources and processes. Seawater intrusion is typically restricted to the coarse-grained units at the base of fining-upward sequences of terrestrial deposits, and at the top of coarsening upward sequences of marine deposits. This results in layered and narrow intrusion fronts. Maintaining the sustainability of Southern California coastal aquifers requires joint management of surface water and groundwater (conjunctive use). This requires new data collection and analyses (including research drilling, modern geohydrologic investigations, and development of detailed computer groundwater models that simulate the supply and demand components separately), implementation of new facilities (including spreading and injection facilities for artificial recharge), and establishment of new institutions and policies that help to sustain the water resources and better manage regional development.

California↗

USGS1997

Comparison of two methods for estimating discharge and nutrient loads from tidally affected reaches of the Myakka and Peace Rivers, west-central Florida

The Myakka and Peace River Basins constitute more than 60 percent of the total inflow area and contribute more than half the total tributary inflow to the Charlotte Harbor estuarine system. Water discharge and nutrient enrichment have been identified as significant concerns in the estuary, and consequently, it is important to accurately estimate the magnitude of discharges and nutrient loads transported by inflows from both rivers. Two methods for estimating discharge and nutrient loads from tidally affected reaches of the Myakka and Peace Rivers were compared. The first method was a tidal-estimation method, in which discharge and nutrient loads were estimated based on stage, water-velocity, discharge, and water-quality data collected near the mouths of the rivers. The second method was a traditional basin-ratio method in which discharge and nutrient loads at the mouths were estimated from discharge and loads measured at upstream stations. Stage and water-velocity data were collected near the river mouths by submersible instruments, deployed in situ, and discharge measurements were made with an acoustic Doppler current profiler. The data collected near the mouths of the Myakka River and Peace River were filtered, using a low-pass filter, to remove daily mixed-tide effects with periods less than about 2 days. The filtered data from near the river mouths were used to calculate daily mean discharge and nutrient loads. These tidal-estimation-method values were then compared to the basin-ratio-method values. Four separate 30-day periods of differing streamflow conditions were chosen for monitoring and comparison. Discharge and nutrient load estimates computed from the tidal-estimation and basin-ratio methods were most similar during high-flow periods. However, during high flow, the values computed from the tidal-estimation method for the Myakka and Peace Rivers were consistently lower than the values computed from the basin-ratio method. There were substantial differences between discharges and nutrient loads computed from the tidal-estimation and basin-ratio methods during low-flow periods. Furthermore, the differences between the methods were not consistent. Discharges and nutrient loads computed from the tidal-estimation method for the Myakka River were higher than those computed from the basin-ratio method, whereas discharges and nutrients loads computed by the tidal-estimation method for the Peace River were not only lower than those computed from the basin-ratio method, but they actually reflected a negative, or upstream, net movement. Short-term tidal measurement results should be used with caution, because antecedent conditions can influence the discharge and nutrient loads. Continuous tidal data collected over a 1- or 2-year period would be necessary to more accurately estimate the tidally affected discharge and nutrient loads for the Myakka and Peace River Basins.

Florida↗

USGS2021

Analysis of Escherichia coli, total recoverable iron, and dissolved selenium concentrations, loading, and identifying data gaps for selected 303(d) listed streams, Grand Valley, western Colorado, 1980–2018

Tributaries to the Colorado River in the Grand Valley in western Colorado (segment COLCLC13b) have been placed on the State of Colorado 303(d) list as impaired for Escherichia coli (E. coli) , total recoverable iron, and dissolved selenium. The Colorado Department of Public Health and Environment Water Quality Control Division is required to develop total maximum daily loads for these constituents in these tributaries. The U.S. Geological Survey, in cooperation with the Grand Valley Drainage District and Colorado Water Conservation Board, conducted a study to (1) characterize concentrations, loads, and load reductions for E. coli , total recoverable iron, and dissolved selenium using existing data and (2) identify water-quality data gaps to inform future monitoring strategies. This study analyzed water-quality and streamflow data for 3 main-stem sites (2 sites along the Colorado River and 1 site along the Gunnison River) and 29 selected sites on tributaries to the Colorado River. Sample data were available at five sites along Adobe Creek and at six sites along Leach Creek, the two tributaries in the study area that are impaired for E. coli . All geometric mean E. coli concentrations at sites along Adobe Creek and Leach Creek exceeded the State recreational use standard of 126 colony forming units per 100 milliliters (CFU/100 mL). In Adobe Creek, E. coli concentrations in samples ranged from 45.7 to more than 2,420 CFU/100 mL (method upper reporting limit for undiluted samples), and geometric mean concentrations at sites ranged from 301 to 1,180 CFU/100 mL. The E. coli concentrations generally increased in the downstream direction in Adobe Creek; however, increases were not seen between all sites. The largest downstream increase in E. coli concentration was measured between the two most upstream sites. In Leach Creek, concentrations of E. coli in samples ranged from 25.9 to more than 2,420 CFU/100 mL, and geometric mean concentrations at sites ranged from 160 to 259 CFU/100 mL. The E. coli concentrations showed no consistent downgradient increase in Leach Creek. In fact, some of the highest E. coli concentrations were measured at the most upstream site, Leach Creek at Summer Hill Drive. Total recoverable iron concentrations and loads were evaluated at 15 tributary sites for samples collected from August 1993 to February 2018. Median total recoverable iron concentrations ranged from 211 to 4,670 micrograms per liter (µg/L). The chronic aquatic-life water-quality standard (1,000 µg/L) was exceeded in most irrigation season (April through October) samples but was rarely exceeded in nonirrigation season (November through March) samples. Concentrations were often an order of magnitude higher in samples collected during irrigation season than in samples collected during nonirrigation season. None of the sites had enough concurrent total recoverable iron and streamflow data to compute annual loads. As with E. coli , the lack of concurrent total recoverable iron and streamflow information represents a data gap, which needs to be addressed to compute annual loads. Dissolved selenium concentrations and loads were evaluated at 20 tributary sites using discrete water-quality data collected 1991–2018. Dissolved selenium concentrations were higher during nonirrigation season than during irrigation season at tributary sites. However, irrigation season dissolved selenium loads were generally higher than nonirrigation selenium loads, because streamflows were higher during irrigation season. Regression analysis was used to estimate daily dissolved selenium concentrations and loads at three main-stem sites for water years (WYs) 1980–2018 (Gunnison River near Grand Junction and Colorado River near Colorado-Utah State Line) and WYs 2002–18 (Colorado River near Cameo). A trend analysis of dissolved selenium concentrations and loads was completed for these sites from the same respective starting dates but ending in 2017. A continuing downward trend in dissolved selenium concentration was observed at all sites and across all seasonal designations of the analysis. The dissolved selenium concentration decreased by 0.12 µg/L from WY 2002 to 2017 at Colorado River near Cameo, representing an 18-percent decrease during the time period. The dissolved selenium concentration at Gunnison River near Grand Junction decreased by 4.2 µg/L from WY 1980 to 2017, representing a 56-percent decrease overall. During the same time period, dissolved selenium concentration at Colorado River near Colorado-Utah State Line decreased by 3.8 µg/L, representing a 56-percent decrease overall. A downward trend in dissolved selenium load was also observed at all sites and across all seasonal designations of the analysis. The relative contribution of dissolved selenium from the Grand Valley near Grand Junction was estimated by comparing loads at main-stem sites bracketing the study area. The two upstream sites, Colorado River near Cameo and Gunnison River near Grand Junction, contributed 60,300 cumulative pounds and 251,000 cumulative pounds, respectively, during WYs 2002–18. At the furthest downstream site, Colorado River near Colorado-Utah State Line, 490,000 cumulative pounds were estimated during the same time period, indicating that the region between Whitewater and State line contributed approximately 179,000 cumulative pounds or a mean annual load of 10,500 lb/yr. Grand Valley dissolved selenium contributions appear to be stable during WYs 2002–18.

Colorado↗

USGS2006

Suspended-sediment yields and stream-channel processes on Judy's Branch watershed in the St. Louis Metro East region in Illinois

Judy's Branch watershed, a small basin (8.64 square miles) in the St. Louis Metro East region in Illinois, was selected as a pilot site to determine suspended-sediment yields and stream-channel processes in the bluffs and American Bottoms (expansive low-lying valley floor in the region). Suspended-sediment and stream-chan-nel data collected and analyzed for Judy's Branch watershed are presented in this report to establish a baseline of data for water-resource managers to evaluate future stream rehabilitation and manage-ment alternatives. The sediment yield analysis determines the amount of sediment being delivered from the watershed and two subwatersheds: an urban tributary and an undeveloped headwater (pri-marily agricultural). The analysis of the subwater-sheds is used to compare the effects of urbanization on sediment yield to the river. The stream-channel contribution to sediment yield was determined by evaluation of the stream-channel processes operat-ing on the streambed and banks of Judy's Branch watershed. Bank stability was related to hydrologic events, bank stratigraphy, and channel geometry through model development and simulation. The average suspended-sediment yield from two upland subwatersheds (drainage areas of 0.23 and 0.40 sq.mi. was 1,163 tons per square mile per year (tons/sq.mi.-year) between July 2000 and June 2004. The suspended-sediment yield at the Route 157 station was 2,523 tons/sq.mi.-year, near the outlet of Judy's Branch watershed (drainage area = 8.33 sq.mi.). This is approximately 1,360 tons/sq.mi.-year greater than the average at the upland stations for the same time period. This result is unexpected in that, generally, the suspended-sediment yield decreases as the watershed area increases because of sediment stored in the channel and flood plain. The difference indicates a possible increase in yield from a source, such as bank retreat, and supports the concept that land-use changes increase stream-flows that may in turn result in higher rates of bank retreat. Utilizing both bank-rod data and resurveyed cross-section data, it was determined that approxi-mately half of the suspended- sediment yield at Route 157 during July 2000-June 2004 came from bank retreat. Given that bank retreat can be a substantial portion of the sediment yield, understanding bank stability processes is important. Bank stability can be assessed mathematically by computing the factor of safety, which is defined by the ratio of the shear strength (resisting force) along the failure surface and the shear stress (driving gravitational force). Once the factor of safety falls below one, the bank theoretically becomes unstable. Bank-stability conditions were related to hydrologic events, bank type, and channel geometry through model develop-ment and simulation. The most common type of bank in the watershed consists of cohesive alluvial soil deposits overlying a stiff glacial till. A stabil-ity chart for different bank types was developed using a bank-stability analysis. Banks steeper than 70 degrees and higher than from 10 to 11.5 feet (depending on bank type) become at risk for mass failure in the watershed under conditions that pro-mote saturation of the bank and a sudden drop in the river level.

Scientific Investigations Report↗

USGS2017

Observations and 3D hydrodynamics-based modeling of decadal-scale shoreline change along the Outer Banks, North Carolina

Long-term decadal-scale shoreline change is an important parameter for quantifying the stability of coastal systems. The decadal-scale coastal change is controlled by processes that occur on short time scales (such as storms) and long-term processes (such as prevailing waves). The ability to predict decadal-scale shoreline change is not well established and the fundamental physical processes controlling this change are not well understood. Here we investigate the processes that create large-scale long-term shoreline change along the Outer Banks of North Carolina, an uninterrupted 60 km stretch of coastline, using both observations and a numerical modeling approach. Shoreline positions for a 24-yr period were derived from aerial photographs of the Outer Banks. Analysis of the shoreline position data showed that, although variable, the shoreline eroded an average of 1.5 m/yr throughout this period. The modeling approach uses a three-dimensional hydrodynamics-based numerical model coupled to a spectral wave model and simulates the full 24-yr time period on a spatial grid running on a short (second scale) time-step to compute the sediment transport patterns. The observations and the model results show similar magnitudes (O(10 5 m 3 /yr)) and patterns of alongshore sediment fluxes. Both the observed and the modeled alongshore sediment transport rates have more rapid changes at the north of our section due to continuously curving coastline, and possible effects of alongshore variations in shelf bathymetry. The southern section with a relatively uniform orientation, on the other hand, has less rapid transport rate changes. Alongshore gradients of the modeled sediment fluxes are translated into shoreline change rates that have agreement in some locations but vary in others. Differences between observations and model results are potentially influenced by geologic framework processes not included in the model. Both the observations and the model results show higher rates of erosion (∼−1 m/yr) averaged over the northern half of the section as compared to the southern half where the observed and modeled averaged net shoreline changes are smaller (<0.1 m/yr). The model indicates accretion in some shallow embayments, whereas observations indicate erosion in these locations. Further analysis identifies that the magnitude of net alongshore sediment transport is strongly dominated by events associated with high wave energy. However, both big- and small- wave events cause shoreline change of the same order of magnitude because it is the gradients in transport, not the magnitude, that are controlling shoreline change. Results also indicate that alongshore momentum is not a simple balance between wave breaking and bottom stress, but also includes processes of horizontal vortex force, horizontal advection and pressure gradient that contribute to long-term alongshore sediment transport. As a comparison to a more simple approach, an empirical formulation for alongshore sediment transport is used. The empirical estimates capture the effect of the breaking term in the hydrodynamics-based model, however, other processes that are accounted for in the hydrodynamics-based model improve the agreement with the observed alongshore sediment transport.

North Carolina↗

USGS2025

Methods for estimating selected low-flow statistics at gaged and ungaged stream sites in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Conservation and Recreation, Office of Water Resources, computed selected at-site streamflow statistics at U.S. Geological Survey streamgages in and near Massachusetts and developed regional regression equations for estimating selected streamflows at ungaged stream sites in Massachusetts. Two sets of regional regression equations were developed: (1) the “mainland” equations, for mainland Massachusetts excluding the area covered by the second set, and (2) the “southeastern” equations, for the Plymouth-Carver-Kingston-Duxbury aquifer area in southeastern Massachusetts and for Cape Cod. The regression equations and at-site statistics may be used by Federal, State, and local water managers in addressing water-resources issues relevant in Massachusetts. Regional regression analyses for the mainland equations were developed to estimate the following 27 streamflow statistics: 99-, 98-, 95-, 90-, 85-, 80-, 75-, 70-, 60-, and 50-percent flow durations; monthly June, July, August, and September 90- and 50-percent flow durations; February, June, and August median of the monthly means; harmonic mean; and medians of the following annual low-flow frequency statistics: 7-day; 7-day, 2-year; 7-day, 10-year; 30-day, 2-year; and 30-day, 10-year. The analyses used 81 streamgages with minimal to no regulations in and near Massachusetts. The regression analyses determined that four basin characteristics—drainage area, combined hydrologic soils A and B, streamflow variability index, and annual mean temperature—were the only significant explanatory variables for the different mainland equations. Regional regression equations were developed for the Plymouth-Carver-Kingston-Duxbury aquifer area in southeastern Massachusetts and Cape Cod, because surface-water drainage areas and groundwater contributing areas do not always coincide in this area of the State. The regression analyses to estimate 10 flow durations from the 99th to 50th percentiles used 18 streamflow sites with some occasional minor regulations—because there are few unregulated streams in southeastern Massachusetts. The analyses determined that groundwater contributing area and storage (combined water bodies and wetlands) were the only significant explanatory variables in the southeastern equations.

Massachusetts↗

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↗