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 505 records · Page 28Linked to original sources

USGS2003

Evaluation of Streamflow, Water Quality, and Permitted and Nonpermitted Loads and Yields in the Raritan River Basin, New Jersey, Water Years 1991-98

Seventeen water-quality constituents were analyzed in samples collected from 21 surface-water sampling sites in the Raritan River Basin during water years 1991-97. Loads were computed for seven constituents. Thirteen constituents have associated instream water-quality standards that are used as reference levels when evaluating the data. Nine of the 13 constituents did not meet water-quality reference levels in all samples at all sites. The constituents that most commonly failed to meet the water-quality reference levels in the 801 samples analyzed were total phosphorus (greater than 0.1 mg/L (milligrams per liter) in 32 percent of samples), fecal coliform bacteria (greater than 400 counts/100 milliliters in 29 percent), hardness (less than 50 mg/L in 21 percent), pH (greater than 8.5 or less than 6.5 in 17 percent), and water temperature in designated trout waters (greater than 20 degrees Celsius in 12 percent of samples). Concentrations of chloride, total dissolved solids, nitrate plus nitrite, and sulfate did not exceed water-quality reference levels in any sample. Results from previous studies on pesticides and volatile organic compounds in streamwater during 1996-98, and organic compounds and trace elements in sediments during 1976-93, were summarized for this study. Concentrations of pesticides in some samples exceeded the relevant standards. Water-quality data varied significantly as season and streamflow changed. Concentrations or values of 12 constituents were significantly higher in the growing season than in the nongrowing season at 1 to 21 sites, and concentrations of 6 constituents were significantly higher in the nongrowing season at 1 to 21 sites. Concentrations or values of seven constituents decreased significantly with increased streamflow, indicating a more significant contribution from base flow or permitted sources than from runoff. Concentrations or values of four constituents increased with increased flow, indicating a more significant contribution from runoff than from base flow or permitted sources. Phosphorus concentrations increased with flow at two sites with no point sources and decreased with flow at five sites with four or more permitted point sources. Concentrations of five constituents did not vary significantly with changes in streamflow at any of the sites. Concentrations of constituents differed significantly between sites. The sites with the most desirable values for the most constituents were Mulhockaway Creek, Spruce Run, Millstone River at Manalapan, Manalapan Brook, and Lamington River at Pottersville. The sites with the least desirable values for the most constituents were Millstone River at Blackwells Mills, Matchaponix Brook, Raritan River at Bound Brook, Neshanic River, and Millstone River at Grovers Mill. The total instream loads of seven constituents - total ammonia plus organic nitrogen (TKN), biochemical oxygen demand (BOD), total dissolved solids (TDS), nitrate plus nitrite (NO3+NO2), total organic carbon (TOC), total phosphorus, and total suspended solids (TSS) - were analyzed at low, median, and high flows. The quantities of total instream load that originated from facilities with permits issued by the New Jersey Department of Environmental Protection to discharge effluent to streams (permitted sources) and from other sources (nonpermitted sources) were estimated for each sampling site. TOC and TSS loads primarily were contributed by nonpermitted sources at all flows. BOD and TDS loads primarily were contributed by nonpermitted sources at median and high flows. At low flow, permitted sources contributed more than one-third of the TDS load at 10 sites and more than one-third of the BOD load at 3 sites. Permitted sources contributed more than one-third of the total phosphorus load at 15 and 14 sites at low and median flows, respectively. Permitted sources accounted for more than one-third of total instream load of NO3+NO2 at low- and median-flow conditions at nearly

Water-Resources Investigations Report↗

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↗

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↗

USGS2002

Simulation of runoff and water quality for 1990 and 2008 land-use conditions in the Reedy Creek watershed, east-central Florida

Hydrologic and water-quality data have been collected within the 177-square-mile Reedy Creek, Florida, watershed, beginning as early as 1939, but the data have not been used to evaluate relations among land use, hydrology, and water quality. A model of the Reedy Creek watershed was developed and applied to the period January 1990 to December 1995 to provide a computational foundation for evaluating the effects of future land-use changes on hydrology and water quality in the watershed. The Hydrological Simulation Program-Fortran (HSPF) model was used to simulate hydrology and water quality of runoff for pervious land areas, impervious land areas, and stream reaches. Six land-use types were used to characterize the hydrology and water quality of pervious and impervious land areas in the Reedy Creek watershed: agriculture, rangeland, forest, wetlands, rapid infiltration basins, and urban areas. Hydrologic routing and water-quality reactions were simulated to characterize hydrologic and water-quality processes and the movement of runoff and its constituents through the main stream channels and their tributaries. Because of the complexity of the stream system within the Reedy Creek Improvement District (RCID) (hydraulic structures, retention ponds) and the anticipated difficulty of modeling the system, an approach of calibrating the model parameters for a subset of the gaged watersheds and confirming the usefulness of the parameters by simulating the remainder of the gaged sites was selected for this study. Two sub-watersheds (Whittenhorse Creek and Davenport Creek) were selected for calibration because both have similar land use to watersheds within the RCID (with the exception of urban areas). Given the lack of available rainfall data, the hydrologic calibration of the Whittenhorse Creek and Davenport Creek sub-watersheds was considered acceptable (for monthly data, correlation coefficients, 0.86 and 0.88, and coefficients of model-fit efficiency, 0.72 and 0.74, respectively). The hydrologic model was tested by applying the parameter sets developed for Whittenhorse Creek and Davenport Creek to other land areas within the Reedy Creek watershed, and by comparing the simulated results to observed data sets for Reedy Creek near Vineland, Bonnet Creek near Vineland, and Reedy Creek near Loughman. The hydrologic model confirmation for Reedy Creek near Vineland (correlation coefficient, 0.91, and coefficient of model fit efficiency, 0.78, for monthly flows) was acceptable. Flows for Bonnet Creek near Vineland were substantially under simulated. Consideration of the ground-water contribution to Bonnet Creek could improve the water balance simulation for Bonnet Creek near Vineland. On longer time scales (monthly or over the 72-month simulation period), simulated discharges for Reedy Creek near Loughman agreed well with observed data (correlation coefficient, 0.88). For monthly flows the coefficient of model-fit efficiency was 0.77. On a shorter time scale (less than a month), however, storm volumes were greatly over simulated and low flows (less than 8 cubic feet per second) were greatly under simulated. A primary reason for the poor results at low flows is the diversion of an unknown amount of water from the RCID at the Bonnet Creek near Kissimmee site. Selection of water-quality constituents for simulation was based primarily on the availability of water-quality data. Dissolved oxygen, nitrogen, and phosphorus species were simulated. Representation of nutrient cycling in HSPF also required simulation of biochemical oxygen demand and phytoplankton populations. The correlation coefficient for simulated and observed daily mean dissolved oxygen concentration values at Reedy Creek near Vineland was 0.633. Simulated time series of total phosphorus, phosphate, ammonia nitrogen, and nitrate nitrogen generally agreed well with periodically observed values for the Whittenhorse Creek and Davenport Creek sites. Simulated water-quality constituents at the Bonnet Creek and Reedy Creek near Vineland sites varied as to how well the values agreed with periodically observed constituent concentrations. Simulated water-quality constituent concentrations for the Reedy Creek near Loughman site generally agreed well with observed constituent concentrations. Simulation of a future land-use scenario for the Reedy Creek watershed was based on the hydrologic and water-quality simulations, projected 2008 land use within the RCID, and assuming no change in existing land use for other areas within the Reedy Creek watershed but external to the RCID. The percentages of forest and urban-impervious land use showed the most change between existing and future land use; forest areas decreased by 50 percent and urban-impervious areas increased by 300 percent. Simulated values of mean total phosphorus, phosphate, ammonia nitrogen, and nitrate nitrogen concentrations for existing and future land-use simulations were within 0.01 milligrams per liter of each other. The simulated maximum daily load increased an average of 10 percent for all constituents. Maximum daily nitrate nitrogen load increased about 17 percent, the greatest increase of all daily constituent loads. Duration curves of daily total phosphorus, phosphate, ammonia nitrogen, and nitrate nitrogen load indicated an increase in the likelihood of exceeding a given load throughout the range of daily constituent loads at Reedy Creek near Loughman.

Florida↗