USGS Science⌕ Search

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.

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

At least 397 records · Page 22Linked to original sources
USGS2019

Development of a flood-inundation map library and precipitation-runoff modeling for the Clear Fork Mohican River in and near Bellville, Ohio

The U.S. Geological Survey (USGS), in cooperation with the Muskingum Watershed Conservancy District, led hydrologic and hydraulic analyses within the Clear Fork Mohican River Basin in and near Bellville, Ohio. The analyses included the development of digital flood-inundation maps for an approximately 2.5-mile reach of the Clear Fork Mohican River and the development of a precipitation-runoff model for a portion of the Clear Fork Mohican River Basin. Data collection for the study involved the installation and operation of 2 streamgages (Clear Fork Mohican River at Bellville, Ohio, and Cedar Fork above Bellville, Ohio); 1 lake-level gage (Clear Fork Reservoir near Lexington, Ohio); 2 precipitation gages (Clear Fork Reservoir near Lexington, Ohio, and Rain Gage at Cedar Fork above Bellville, Ohio); and 12 submersible pressure transducers on Clear Fork Mohican River and 4 of its tributaries. Data collection also included field surveys of hydraulic structures and channel cross sections. Flood profiles were computed for the 2.5-mile reach of the Clear Fork Mohican River by means of a one-dimensional step-backwater model. The model was calibrated to 16 measured events and to a portion (stages 9 to 11 feet) of the current stage-streamflow relation at the USGS streamgage Clear Fork Mohican River at Bellville, Ohio, and to stage recorded at a submersible pressure transducer site near the downstream study limit. After calibration the step-backwater model was used to compute nine flood profiles for stages ranging from 9 to 17 feet. The flood profiles were then used in combination with a digital elevation model to delineate the area that would be inundated at each stage. A precipitation-runoff model was developed and calibrated using data from the streamgage, precipitation gage, and 11 submersible pressure transducers. The modeling included data during 10 runoff events that were used for model calibration and validation, with focus on 6 events. The Nash-Sutcliffe model efficiency coefficients for six peak streamflow events ranged from 0.459 to 0.851. The models produced by this study can be used to assess possible flood mitigation options and define flood hazard areas that could contribute to the protection of life and property. The availability of flood-inundation maps, internet information from USGS streamgages, and forecasted stages from the National Weather Service could provide emergency management personnel and residents with information on forecasting floods, appropriate flood response activities, and post-flood recovery efforts.

Ohio↗
USGS1997

A comparison of drainage basin nutrient inputs with instream nutrient loads for seven rivers in Georgia and Florida, 1986-90

Instream nutrient loads of the Altamaha, Suwannee, St. Johns, Satilla, Ogeechee, Withlacoochee, and Ochlockonee River Basins were computed and compared with nutrient inputs for each basin for the period 1986-90. Nutrient constituents that were considered included nitrate, ammonia, organic nitrogen, and total phosphorus. Sources of nutrients considered for this analysis included atmospheric deposition, fertilizer, animal waste, wastewater-treatment plant discharge, and septic discharge. The mean nitrogen input ranged from 2,400 kilograms per year per square kilometer (kg/yr)km2 in the Withlacoochee River Basin to 5,470 (kg/yr)km2 in the Altamaha River Basin. The Satilla and Ochlockonee River Basins also had large amounts of nitrogen input per unit area, totaling 5,430 and 4,920 (kg/yr)km2, respectively. Fertilizer or animal waste, as sources of nitrogen, predominated in all basins. Atmospheric deposition contributed less than one-fourth of the mean total nitrogen input to all basins and was consistently the third largest input in all but the Ogeechee River Basin, where it was the second largest. The mean total phosphorus input ranged from 331 (kg/yr)km2 in the Withlacoochee River Basin to 1,380 (kg/yr)km2 in both the Altamaha and Satilla River Basins. The Ochlockonee River Basin had a phosphorus input of 1,140 (kg/yr)km2. Per unit area, the Suwannee River discharged the highest instream mean total nitrogen and phosphorus loads and also discharged higher instream nitrate loads per unit area than the other six rivers. Phosphorus loads in stream discharge were highest in the Suwannee and Ochlockonee Rivers. The ratio of nutrient outputs to inputs for the seven studied rivers ranged from 4.2 to 14.9 percent, with the St. Johns (14.9 percent) and Suwannee (12.1 percent) Rivers having significantly higher percentages than those from the other basins. The output/input percentages for mean total phosphorus ranged from 1.0 to 7.0 percent, with the St. Johns (6.2 percent) and Suwannee (7.0 percent) Rivers exporting the highest percentage of phosphorus. Although instream nutrient loads constitute only one of the various pathways nutrients may take in leaving a river basin, only a relatively small part of nutrient input to the basin leaves the basin in stream discharge for the major coastal rivers examined in this study. The actual amount of nutrient transported in a river basin depends on the ways in which nutrients are physically handled, geographically distributed, and chemically assimilated within a river basin.

Florida, Georgia↗
USGS1997

Retention of sediments and nutrients in Jackson Creek wetland near Delavan Lake, Wisconsin, 1993-95

Excessive inflow of sediment and nutrients has contributed to eutrophication of Delavan Lake in southeastern Wisconsin. A wetland restoration project at the mouth of the lake's principal inflowing stream system was implemented in 1992, primarily for the purpose of reducing sediment and nutrient loading to the lake. This function of the wetland system, which includes three retention ponds, was examined in a 32-month study (February 1993-September 1995) that included continuous monitoring of discharge of two inflowing streams and the wetland outflow, intermittent monitoring of discharge of a third inflowing stream, monthly and storm-runoff sampling to determine sediment and nutrient concentrations, and measurement of sediment accumulation in three constructed sediment retention ponds. Inflow and outflow loads of sediment and nutrients were computed from the discharge and concentration data. The wetland system consistently retained suspended sediments during the study period. Volumetric surveys of the ponds showed sediment accumulation each year, and during the 32 months, 46 percent of the inflow sediment load was retained in the wetland and the retention ponds. Retention of nutrients in the wetland system, however, was of lesser magnitude and much greater seasonal variability. Over time periods of one year or more, most nutrient forms were retained in the wetland system. Nutrient retention, relative to the 32-month inflow loads was: 19 percent of total phosphorus, 11 percent of dissolved orthophosphate, 8 percent of total ammonia plus organic nitrogen, and 0.4 percent of dissolved nitrite plus nitrate nitrogen. Of the nutrient forms measured, only dissolved ammonia, whose 32-month outflow from the wetland system was 22 percent greater than inflow, showed a long-term net release from the system. However, net releases over shorter time periods were commonly observed for all nutrients, and these occurred frequently during the growing season, suggesting probable availability to algal and macrophyte communities downstream. Awareness of the variability and complexity of the nutrient-trapping function of a wetland system can help maintain realistic expectations of the benefits of wetland restoration projects and can be valuable in developing more effective management practices.

Wisconsin↗
USGS2009

Description and evaluation of numerical groundwater flow models for the Edwards Aquifer, south-central Texas

A substantial number of public water system wells in south-central Texas withdraw groundwater from the karstic, highly productive Edwards aquifer. However, the use of numerical groundwater flow models to aid in the delineation of contributing areas for public water system wells in the Edwards aquifer is problematic because of the complex hydrogeologic framework and the presence of conduit-dominated flow paths in the aquifer. The U.S. Geological Survey, in cooperation with the Texas Commission on Environmental Quality, evaluated six published numerical groundwater flow models (all deterministic) that have been developed for the Edwards aquifer San Antonio segment or Barton Springs segment, or both. This report describes the models developed and evaluates each with respect to accessibility and ease of use, range of conditions simulated, accuracy of simulations, agreement with dye-tracer tests, and limitations of the models. These models are (1) GWSIM model of the San Antonio segment, a FORTRAN computer-model code that pre-dates the development of MODFLOW; (2) MODFLOW conduit-flow model of San Antonio and Barton Springs segments; (3) MODFLOW diffuse-flow model of San Antonio and Barton Springs segments; (4) MODFLOW Groundwater Availability Modeling [GAM] model of the Barton Springs segment; (5) MODFLOW recalibrated GAM model of the Barton Springs segment; and (6) MODFLOW–DCM (dual conductivity model) conduit model of the Barton Springs segment. The GWSIM model code is not commercially available, is limited in its application to the San Antonio segment of the Edwards aquifer, and lacks the ability of MODFLOW to easily incorporate newly developed processes and packages to better simulate hydrologic processes. MODFLOW is a widely used and tested code for numerical modeling of groundwater flow, is well documented, and is in the public domain. These attributes make MODFLOW a preferred code with regard to accessibility and ease of use. The MODFLOW conduit-flow model incorporates improvements over previous models by using (1) a user-friendly interface, (2) updated computer codes (MODFLOW–96 and MODFLOW–2000), (3) a finer grid resolution, (4) less-restrictive boundary conditions, (5) an improved discretization of hydraulic conductivity, (6) more accurate estimates of pumping stresses, (7) a long transient simulation period (54 years, 1947–2000), and (8) a refined representation of high-permeability zones or conduits. All of the models except the MODFLOW–DCM conduit model have limitations resulting from the use of Darcy’s law to simulate groundwater flow in a karst aquifer system where non-Darcian, turbulent flow might actually dominate. The MODFLOW–DCM conduit model is an improvement in the ability to simulate karst-like flow conditions in conjunction with porous-media-type matrix flow. However, the MODFLOW–DCM conduit model has had limited application and testing and currently (2008) lacks commercially available pre- and post-processors. The MODFLOW conduit-flow and diffuse-flow Edwards aquifer models are limited by the lack of calibration for the northern part of the Barton Springs segment (Travis County) and their reliance on the use of the calibrated hydraulic conductivity and storativity values from the calibrated Barton Springs segment GAM model. The major limitation of the Barton Springs segment GAM and recalibrated GAM models is that they were calibrated to match measured water levels and springflows for a restrictive range of hydrologic conditions, with each model having different hydraulic conductivity and storativity values appropriate to the hydrologic conditions that were simulated. The need for two different sets of hydraulic conductivity and storativity values increases the uncertainty associated with the accuracy of either set of values, illustrates the non-uniqueness of the model solution, and probably most importantly demonstrates the limitations of using a one-layer model to represent the heterogeneous hydrostratigraphic units composing the Edwards aquifer. In general, the best matches or agreement between groundwater flow directions inferred by numerical model simulation, and by dye-tracer tests, are observed where model outputs accurately reproduce the configuration of the potentiometric surface with regard to the positions of major and minor groundwater troughs and divides. None of the models, with the possible exception of the MODFLOW–DCM conduit model, has a documented capability to accurately simulate travel times for conduit-dominated velocities in the Edwards aquifer. Public water system assessments of wells in the Barton Springs segment of the Edwards aquifer, and elsewhere where conduit-flow conditions are thought to dominate aquifer hydraulic behavior, might be enhanced by use of either the MODFLOW–DCM model or the newly developed U.S. Geological Survey MODFLOW Conduit-Flow Process module for MODFLOW–2005, because each incorporates a type of dual or triple hydraulic conductivity approach and has the capability to explicitly simulate turbulent flow and conduit hydraulic characteristics.

Texas↗
USGS1986

Field observations of bed shear stress and sediment resuspension on continental shelves, Alaska and California

Bed shear stress was estimated using wave and current measurements obtained with the GEOPROBE bottom-tripod system during resuspension events in Norton Sound, Alaska, and on the northern California shelf. The boundary-layer model of Grant and Madsen (1979, Journal of Geophysical Research, 84, 1797-1808) was used to compute the bed shear stress under combined wave-generated and quasi-steady currents. Resuspension events were identified by sudden, large increases in light scattering at 1.9 m above the sea floor. The shear-stress values were used to compute the Shields parameter (??). The results for Norton Sound are in excellent agreement with the Shields threshold criterion; the data for the California shelf plot somewhat above the Shields threshold curve, though generally within the scatter envelope. Although the surface sediments in each area contain substantial fine-grained fractions (mean diameters were 0.007 cm in Norton Sound and 0.002 cm on the California shelf), the results do not indicate significant cohesion, because the sediment was entrained at bed shear-stress values close to those predicted by the modified Shields curve for cohesionless fine-grained particles. We suspect that frequent wave stirring and observed plowing of the surface sediment by benthonic animals maintain a high water content and contribute to the ease with which these materials are resuspended.

Alaska, California↗
USGS1971

An evaluation of procedures used in computing chemical denudation rates

Computations of chemical denudation rates (that is, the average rate of lowering of the earth's crust by chemical processes) should be based only upon those components of dissolved loads of streams that are derived from chemical weathering of rocks and soils, even though identification of those components is difficult and imprecise. Past computations of chemical denudation rates exaggerated by 1.4 to 2.4 times the significance of this process because they were calculated from total dissolved loads. Identification of separate denudation and nondenudation components in a water analysis can be done most readily for streams in crystalline silicate rock terrane where commonly 25 to 55 percent of the dissolved solids are HCO 3 −1 , SO 4 −2 , Cl −1 , and NO 3 −1 that cannot be derived from the underlying rock; instead these ions are derived either directly from the atmosphere or from atmosphere-biosphere interactions. Most of the remaining solutes are derived from chemical weathering or rainfall, or both. Existing data are often insufficient to separate these two components, but atmospheric contributions to dissolved loads of dilute natural waters are substantial. The commonly observed increased concentration of dissolved solids in chemically polluted rivers is produced by concentration of natural dissolved loads by evaporation, and by addition of solutes not present in natural waters. Concentration of natural solutes represents a modification of natural denudation processes that does not change the rate of denudation. Addition of new solutes represents man-induced, accelerated chemical denudation. Extrapolation of modern chemical denudation rates to geologic situations is extremely imprecise because of uncertainties in assessing man's influence on present dissolved loads, in reconstructing the number and magnitude of past climatic changes, and in evaluating the effects of evolving flora and fauna.

Bulletin of the Geological Society of America↗
USGS2011

Mechanics of flow and sediment transport in delta distributary channels

Predicting the planform and dimensions of a channel downstream from a confluence of two smaller channels with known sediment and water supplies is a fundamental, well-studied problem in geomorphology and engineering. An analogous but less well understood problem is found well downstream of such confluences, where large river channels split into two or more distributary channels on a river delta. In this case, both the flow and sediment supplies in the downstream distributaries are set by the dynamics near the bifurcation of the upstream channel and by downstream boundary conditions. Over time, the pattern of erosion and deposition in the distributary channels gives rise to variations in the amount of water and sediment routed into them. In the simplest case, this results in channel switching on deltas, but in a more general sense these dynamics produce a rich suite of interesting morphologic change contributing both to the evolution of the channel distributary network and the overall evolution of the delta. As part of a study to develop a better understanding of these processes, we conducted a field study measuring the detailed morphology of the Hong-Luoc junction on the Red River downstream of Hanoi, Vietnam. This junction was selected for such a study because it has a 1,000-year history, modern observations suggest that it is currently switching (changing the proportion of sediment and streamflow provided to each of the distributary channels), and hydrologic configuration of the junction allows for the study of two bifurcations and one confluence in a single junction complex. In this paper, our morphologic observations are used in computational flow models to show how flow and sediment transport changes as a function of total discharge upstream of the junction. This is a key component of understanding how the junction attains stability over a range of flows or how imbalances in the distribution of flow and sediment transport lead to destabilization of the channel bifurcation.

Conference Paper↗
USGS2001

Hydrologic conditions and water quality in an agricultural area in Kleberg and Nueces Counties, Texas, 1996-98

During 1996–98, rainfall and runoff were monitored on a 49,680-acre agricultural watershed in Kleberg and Nueces Counties in South Texas. Nineteen rainfall samples were analyzed for selected nutrients, and runoff samples from 29 storms were analyzed for major ions, nutrients, and pesticides. Loads of nutrients in rainfall and loads of nutrients and pesticides in runoff were computed. For a 40,540-acre part of the watershed (lower study area), constituent loads entering the watershed in rainfall, in runoff from the upper study area, and from agricultural chemical applications to the lower study area were compared with runoff loads exiting the lower study area. Total rainfall for 1996–98 averaged 25.86 inches per year, which is less than the long-term annual average rainfall of 29.80 inches for the area. Rainfall and runoff during 1996–98 were typical of historical patterns, with periods of below average rainfall and runoff interspersed with extreme events. Five individual storms accounted for about 38 percent of the total rainfall and 94 percent of the total runoff. During the 3-year study, the total nitrogen runoff yield from the lower study area was 1.3 pounds per acre per year, compared with 49 pounds per acre per year applied as fertilizer and 3.1 pounds per acre per year from rainfall. While almost all of the fertilizer and rainfall nitrogen was ammonia and nitrate, most of the nitrogen in runoff was particulate organic nitrogen, associated with crop residue. Total nitrogen exiting the lower study area in surface-water runoff was about 2.5 percent of the nitrogen inputs (fertilizer and rainfall nitrogen). Annual deposition of total nitrogen entering the lower study area in rainfall exceeded net yields of total nitrogen exiting the watershed in runoff because most of the rainfall does not contribute to runoff. During the study, the total phosphorus runoff yield from the lower study area was 0.48 pound per acre per year compared with 4.2 pounds per acre per year applied as fertilizer and 0.03 pound per acre per year from rainfall. Twenty-one pesticides were detected in runoff with varying degrees of frequency during the study. The herbicide atrazine was detected in all runoff samples. All of the most frequently detected pesticides (atrazine, trifluralin, simazine, pendimethalin, and diuron) exhibited higher concentrations during the pre-harvest period (March–May) than during the post-harvest period (August–October). During 1996–98, an average of 0.37 pound per acre per year of atrazine was applied to the lower study area. During the same period, 0.0027 pound per acre per year of atrazine and its breakdown product deethylatrazine exited the lower study area in runoff (about 0.7 percent of the total atrazine applied to the cropland). During 1997, when heavy rainfall occurred during the months of April and May, the atrazine plus deethylatrazine exiting the lower study area was 1.8 percent of the applied atrazine. The 1996–98 average sediment yield was 610 pounds per acre per year. Sediment loads from the study area are associated with large storm events. Of the 45,300 tons of sediment transported from the study area during 1996–98 about 87 percent was transported during the three largest runoff events (April 1997, October 1997, and October 1998). Runoff-weighted average concentrations were computed for selected nutrients and pesticides. The 1996–98 runoff-weighted concentrations for total nitrog

Texas↗
USGS2016

Map visualization of groundwater withdrawals at the sub-basin scale

A simple method is proposed to visualize the magnitude of groundwater withdrawals from wells relative to user-defined water-resource metrics. The map is solely an illustration of the withdrawal magnitudes, spatially centered on wells—it is not capture zones or source areas contributing recharge to wells. Common practice is to scale the size (area) of withdrawal well symbols proportional to pumping rate. Symbols are drawn large enough to be visible, but not so large that they overlap excessively. In contrast to such graphics-based symbol sizes, the proposed method uses a depth-rate index (length per time) to visualize the well withdrawal rates by volumetrically consistent areas, called “footprints”. The area of each individual well’s footprint is the withdrawal rate divided by the depth-rate index. For example, the groundwater recharge rate could be used as a depth-rate index to show how large withdrawals are relative to that recharge. To account for the interference of nearby wells, composite footprints are computed by iterative nearest-neighbor distribution of excess withdrawals on a computational and display grid having uniform square cells. The map shows circular footprints at individual isolated wells and merged footprint areas where wells’ individual footprints overlap. Examples are presented for depth-rate indexes corresponding to recharge, to spatially variable stream baseflow (normalized by basin area), and to the average rate of water-table decline (scaled by specific yield). These depth-rate indexes are water-resource metrics, and the footprints visualize the magnitude of withdrawals relative to these metrics.

Hydrogeology Journal↗
USGS2015

Analyzing high resolution topography for advancing the understanding of mass and energy transfer through landscapes: A review

The study of mass and energy transfer across landscapes has recently evolved to comprehensive considerations acknowledging the role of biota and humans as geomorphic agents, as well as the importance of small-scale landscape features. A contributing and supporting factor to this evolution is the emergence over the last two decades of technologies able to acquire high resolution topography (HRT) (meter and sub-meter resolution) data. Landscape features can now be captured at an appropriately fine spatial resolution at which surface processes operate; this has revolutionized the way we study Earth-surface processes. The wealth of information contained in HRT also presents considerable challenges. For example, selection of the most appropriate type of HRT data for a given application is not trivial. No definitive approach exists for identifying and filtering erroneous or unwanted data, yet inappropriate filtering can create artifacts or eliminate/distort critical features. Estimates of errors and uncertainty are often poorly defined and typically fail to represent the spatial heterogeneity of the dataset, which may introduce bias or error for many analyses. For ease of use, gridded products are typically preferred rather than the more information-rich point cloud representations. Thus many users take advantage of only a fraction of the available data, which has furthermore been subjected to a series of operations often not known or investigated by the user. Lastly, standard HRT analysis work-flows are yet to be established for many popular HRT operations, which has contributed to the limited use of point cloud data. In this review, we identify key research questions relevant to the Earth-surface processes community within the theme of mass and energy transfer across landscapes and offer guidance on how to identify the most appropriate topographic data type for the analysis of interest. We describe the operations commonly performed from raw data to raster products and we identify key considerations and suggest appropriate work-flows for each, pointing to useful resources and available tools. Future research directions should stimulate further development of tools that take advantage of the wealth of information contained in the HRT data and address the present and upcoming research needs such as the ability to filter out unwanted data, compute spatially variable estimates of uncertainty and perform multi-scale analyses. While we focus primarily on HRT applications for mass and energy transfer, we envision this review to be relevant beyond the Earth-surface processes community for a much broader range of applications involving the analysis of HRT.

Earth-Science Reviews↗
USGS2023

Mode and provenance of sediment deposition on a transgressive marsh

In this study, we use foraminifera as environmental indicators to aid in computing the historical volumetric inputs of estuarine sediments to adjacent marsh. These data can help assess the importance of estuarine sediment inputs to marsh accretion. The Grand Bay system (GBS), located on the southern coast of Alabama and Mississippi, has been described as a “self-cannibalizing bay-marsh complex” due to the disproportionately large amount of suspended sediment exported from the GBS relative to the amount of sediment imported into tidal channels and available to the marsh. Despite this sediment limitation in the marshes, mass sediment budgets along shoreline-proximal marsh sites suggest depositional fluxes on the marsh are nearly equivalent to erosional losses at decadal scales. Geochronologies and foraminiferal census data for two shore-normal transects show that estuarine sediment contributes significant portions (28-65%) of this depositional flux. Contrast in shoreline morphology (i.e., orientation, irregularity, and exposure) and adjacent environments (i.e., tidal creeks and mud flats) may influence the proportion of estuarine sediment delivered to marsh. Improvements of sediment provenance tracers in these environments will continually improve the understanding of sediment reworking and transport in transgressive marsh-estuarine systems such as GBS and the role they play in marsh resiliency to sea-level change.

Alabama, Mississippi↗
USGS2019

Regional Global Navigation Satellite System networks for crustal deformation monitoring

Regional networks of Global Navigation Satellite System (GNSS) stations cover seismically and volcanically active areas throughout the United States. Data from these networks have been used to produce high‐precision, three‐component velocity fields covering broad geographic regions as well as position time series that track time‐varying crustal deformation. This information has contributed to assessing interseismic strain accumulation and related seismic hazard, revealed previously unknown occurrences of aseismic fault slip, constrained coseismic slip estimates, and enabled monitoring of volcanic unrest and postseismic deformation. In addition, real‐time GNSS data are now widely available. Such observations proved invaluable for tracking the rapidly evolving eruption of Kīlauea in 2018. Real‐time earthquake source modeling using GNSS data is being incorporated into tsunami warning systems, and a vigorous research effort is focused on quantifying the contribution that real‐time GNSS can make to improve earthquake early warnings as part of the Advanced National Seismic System ShakeAlert system. Real‐time GNSS data can also aid in the tracking of ionospheric disturbances and precipitable water vapor for weather forecasting. Although regional GNSS and seismic networks generally have been established independently, their spatial footprints often overlap, and in some cases the same institution operates both types of networks. Further integration of GNSS and seismic networks would promote joint use of the two data types to better characterize earthquake sources and ground motion as well as offer opportunities for more efficient network operations. Looking ahead, upgrading network stations to leverage new GNSS technology could enable more precise positioning and robust real‐time operations. New computational approaches such as machine learning have the potential to enable full utilization of the large amounts of data generated by continuous GNSS networks. Development of seafloor Global Positioning System‐acoustic networks would provide unique information for fundamental and applied research on subduction zone seismic hazard and, potentially, monitoring.

Seismological Research Letters↗
USGS1964

Water quality of the Swatara Creek Basin, PA

The Swatara Creek of the Susquehanna River Basin is the farthest downstream sub-basin that drains acid water (pH of 4.5 or less) from anthracite coal mines. The Swatara Creek drainage area includes 567 square miles of parts of Schuylkill, Berks, Lebanon, and Dauphin Counties in Pennsylvania. To learn what environmental factors and dissolved constituents in water were influencing the quality of Swatara Creek, a reconnaissance of the basin was begun during the summer of 1958. Most of the surface streams and the wells adjacent to the principal tributaries of the Creek were sampled for chemical analysis. Effluents from aquifers underlying the basin were chemically analyzed because ground water is the basic source of supply to surface streams in the Swatara Creek basin. When there is little runoff during droughts, ground water has a dominating influence on the quality of surface water. Field tests showed that all ground water in the basin was non-acidic. However, several streams were acidic. Sources of acidity in these streams were traced to the overflow of impounded water in unworked coal mines. Acidic mine effluents and washings from coal breakers were detected downstream in Swatara Creek as far as Harper Tavern, although the pH at Harper Tavern infrequently went below 6.0. Suspended-sediment sampling at this location showed the mean daily concentration ranged from 2 to 500 ppm. The concentration of suspended sediment is influenced by runoff and land use, and at Harper Tavern it consisted of natural sediments and coal wastes. The average daily suspended-sediment discharge there during the period May 8 to September 30, 1959, was 109 tons per day, and the computed annual suspended-sediment load, 450 tons per square mile. Only moderate treatment would be required to restore the quality of Swatara Creek at Harper Tavern for many uses. Above Ravine, however, the quality of the Creek is generally acidic and, therefore, of limited usefulness to public supplies, industries and recreation. In general, the quality of Swatara Creek improves after it mixes with water from the Upper Little and Lower Little Swatara Creeks, which converge with the main stream near Pine Grove. Jonestown is the first downstream location where Swatara Creek contains bicarbonate ion most of the time, and for the remaining downstream length of the stream, the concentration of bicarbonate progressively increases. Before the stream enters the Susquehanna River, chemical and diluting processes contributed by tributaries change the acidic calcium sulfate water, which characterizes the upper Swatara, to a calcium bicarbonate water. A major tributary to Swatara Creek is Quittapahilla Creek, which drains a limestone region and has alkaline characteristics. Effluents from a sewage treatment plant are discharged into this stream west of Lebanon. Adjacent to the Creek are limestone quarries and during the recovery of limestone, ground water seeps into the mining areas. This water is pumped to upper levels and flows over the land surface into Quittapahilla Creek. As compared with the 1940's, the quality of Swatara Creek is better today, and the water is suitable for more uses. In large part, this improvement is due to curtailment of anthracite coal mining and because of the controls imposed on new mines, stripping mines, and the related coal mining operations, by the Pennsylvania Sanitary Water Board. Thus, today (1962) smaller amounts of coal mine wastes are more effectively flushed and scoured away with each successive runoff during storms that affect the drainage basin. Natural processes neutralizing acid water in the stream by infiltration of alkaline ground water through springs and through the streambed are also indicated.

Pennsylvania↗
USGS2021

Regression relations and long-term water-quality constituent concentrations, loads, yields, and trends in the North Fork Ninnescah River, south-central Kansas, 1999–2019

Cheney Reservoir, in south-central Kansas, is the primary water supply for the city of Wichita, Kansas. The North Fork Ninnescah River is the largest tributary to Cheney Reservoir and contributes about 70 percent of the inflow. The U.S. Geological Survey, in cooperation with the City of Wichita, has been continuously monitoring water quality (including water temperature, specific conductance, pH, dissolved oxygen, and turbidity) on the North Fork Ninnescah River upstream from Cheney Reservoir (U.S. Geological Survey site 07144780) since November 1998. Continued data collection would be beneficial to update and describe changing water-quality conditions in the drainage basin and in the reservoir over time. Regression models were developed to describe relations between discretely measured constituent concentrations and continuously measured physical properties. The models updated in this report include total suspended solids (TSS), suspended-sediment concentration (SSC), nitrate plus nitrite, nitrate, orthophosphate (OP), total phosphorus (TP), and total organic carbon (TOC). Daily computed concentrations for TSS, TP, and nitrate plus nitrite during 1999–2019 were compared with Cheney Reservoir Task Force (CRTF) goals for base-flow and runoff conditions. CRTF goals for base-flow concentrations were exceeded more frequently (70 to 99.9 percent of the time) than runoff goals (0 to 11 percent of the time). Except for 2012, annual mean TSS concentrations exceeded the base-flow goal every year. Nitrate plus nitrite and TP annual mean concentrations exceeded the base-flow goals every year. TSS and nitrate plus nitrite annual mean concentrations during runoff conditions never exceeded the CRTF runoff goal. TP annual mean concentrations during runoff conditions only exceeded the CRTF runoff goal during 2002. Sedimentation is progressively reducing the storage capacity of Cheney Reservoir. During 1999–2019, 55 percent of the computed suspended-sediment load was transported during the top 1 percent of loading days (76 days); 22 percent of the total load was transported in the top 10 loading days, indicating that substantial parts of suspended-sediment loads continue to be delivered during disproportionately small periods in Cheney Reservoir. Successful sediment management efforts necessitate reduction techniques that account for these large load events. Flow-normalized concentrations and fluxes were computed during 1999 through 2019 using Weighted Regressions on Time, Discharge, and Season (WRTDS) statistical models and WRTDS bootstrap tests. Flow-normalized concentrations of TSS, SSC, OP, TP, and TOC had upward trend probabilities; conversely, nitrate plus nitrite had a downward trend. Flow-normalized fluxes for OP, TP, and TOC had an upward trend. No discernible patterns were identified for flow-normalized flux of TSS or suspended sediment. Nitrate plus nitrite flow-normalized flux indicated a downward trend. Flow-normalized concentrations for TSS were less than the CRTF long-term goal of 100 milligrams per liter (mg/L), but the upward trend indicated the long-term goal may be exceeded if no changes are made. Flow-normalized TP concentrations exceeded the CRTF long-term goal (0.1 mg/L) and were assigned a very likely upward trend. Flow-normalized nitrate plus nitrite concentrations exceeded the CRTF long-term goal of 1.2 mg/L during the beginning of the study period, then were less than the CRTF goal for the remainder of the study; however, during 2010–19 flow-normalized concentrations increased by 6 percent. Linking water-quality changes to causal factors requires consistent monitoring before, during, and after changes; this presents challenges related to length and frequency of data collection and available concomitant land-use and conservation practice data. As such, attribution of water-quality trends to land-use changes or conservation practices was not possible for this study because of a lack of land-use and conservation practice data. Additionally, because precipitation frequency and intensity are projected to continue to increase in the Great Plains region, accounting for extreme episodic events may be an important consideration in future sediment and nutrient load reduction plans.

Kansas↗
USGS1961

Hydraulic and hydrologic aspects of flood-plain planning

The valid incentives compelling occupation of the flood plain, up to and eve n into the stream channel, undoubtedly have contributed greatly to the development of the country. But the result has been a heritage of flood disaster, suffering, and enormous costs. Flood destruction awakened a consciousness toward reduction and elimination of flood hazards, originally manifested in the protection of existing developments. More recently, increased knowledge of the problem has shown the impracticability of permitting development that requires costly flood protect/on. The idea of flood zoning, or flood-plain planning, has received greater impetus as a result of this realization. This study shows how hydraulic and hydrologic data concerning the flood regimen of a stream can be used in appraising its flood potential and the risk inherent in occupation of its flood plain. The approach involves the study of flood magnitudes as recorded or computed; flood frequencies based1 on experience shown by many years of gaging-station record; use of existing or computed stagedischarge relations and flood profiles; and, where required, the preparation of flood-zone maps to show the areas inundated by floods of several magnitudes and frequencies. The planner can delineate areas subject to inundation by floods o* specific recurrence intervals for three conditions: (a) for the immediate vicinity of a gaging station; (b) for a gaged stream at a considerable distance from a gaging station; and (c) for an ungaged stream. The average depth for a flood of specific frequency can be estimated on the basis of simple measurements of area of drainage basin, width of channel, and slope of streambed. This simplified approach should be useful in the initial stages of flood-plain planning. Brief discussions are included on various types of flood hazards, the effects of urbanization on flood runoff, and zoning considerations.

Pennsylvania↗
USGS1993

Geomorphic and hydraulic assessment of the Bear River in and near Evanston, Wyoming

Geomorphic and hydraulic characteristics of the Bear River in and near Evanston, Wyoming, were assessed to assist planners in stabilizing the river channel. Present-day channel instability is the result of both human-made and natural factors. The primary factor is channelization of the river in Evanston, where several meander loops were cut off artificially during early development of the city. Other contributing factors include channel-width constrictions, bank stabilization, isolated bend cutoffs upstream from the city, and flooding in 1983 and 1984. A geomorphic analysis of bankfull-channel pattern, based on four aerial photographs taken during 1946-86, quantified geomorphic properties (reach sinuosity, bend sinuosity, bend radius of curvature, and bed length) that are characteristic of the study reach. The reach sinuosity of reach 2 (the channelized reach in Evanston) was 1.18 in 1986 and remained about the same throughout the period (1946-86). The reach sinuosity of reach 2 prior to channelization was substantially larger, about 2.3 as determined from maps prepared before 1946. Hydraulic analysis of the present-day channel (surveyed 1981-87) using a one-dimensional water-surface-profile computer model identified a bankfull discharge for the study reach of 3,600 cu ft/sec. A comparison of bankfull hydraulic properties for reaches 1, 2, and 3 indicated that the effects in reach 2 of channelization and channel-width constriction--increased slope, faster velocities, and greater hydraulic radii. The present-day channel slope in reach 2 is 0.00518 ft/ft, whereas a more stable slope would be between 0.00431 ft/ft (present-day slope in reach 1) and 0.00486 ft/ft (present-day slope in reach 3).

Water-Resources Investigations Report↗
USGS1958

Hydraulic and hydrologic aspects of flood-plain planning

The valid incentives compelling occupation of the flood plain, up to and even into the stream channel, undoubtedly have contributed greatly to the development of the country. But the result has been a heritage of flood disaster, suffering, and enormous costs. Flood destruction awakened a consciousness toward reduction and elimination of flood hazards, originally manifested in the protection of existing developments. More recently, increased information and realization of the problem have shown the impracticability of permitting developments that require costly flood protection. The idea of flood zoning, or flood-plain planning, has received greater impetus as a result of such realization. This study shows how hydraulic and hydrologic data concerning the flood regimen of a stream can be used in appraising its flood potential and the risk inherent in occupation of its flood plain. The approach involves the study of flood magnitudes as recorded or computed; flood frequencies based on the experience shown by many years of gaging-station record; use of existing or computed stage-discharge relations and flood profiles; and, where required, the preparation of flood-zone maps to show the areas inundated by floods of several magnitudes and frequencies. Methods are outlined that will enable the planner to delineate areas subject to inundation by floods of specific recurrence intervals for three conditions; (1) for the immediate vicinity of a gaging station, (2) for a gaged stream at a considerable distance from a gaging station, and (3) for an ungaged stream. Curves have been developed from which the average depth for a flood of specific frequency can be estimated on the basis of the average depth for bankfull flow; this simplified approach should be very useful in the initial stages of flood-plain planning. Brief discussions are included on various types of flood hazards, the effects of urbanization on flood runoff, and zoning considerations.

Pennsylvania↗
USGS2022

Hydraulics of freshwater mussel habitat in select reaches of the Big River, Missouri

The Big River is a tributary to the Meramec River in south-central Missouri. It drains an area that has been historically one of the largest lead producers in the world, and associated mine wastes have contaminated sediments in much of the river corridor. This study investigated hydraulic conditions in four study reaches to evaluate the potential contribution of physical habitat dynamics to mechanical and physiological stress on native mussel populations. We quantified hydraulic conditions and relative bed stability in previously identified and delineated mussel habitats (MHs) and in the surrounding reaches to refine understanding of the reach-scale (about 1 kilometer) hydraulic characteristics that affect the distribution of mussel aggregations in the river. Two-dimensional hydrodynamic models were compiled for discharge scenarios from base flow (90-percent flow exceedance) to the approximate bankfull discharge (2-year mean return interval peak flow) for the reaches. Discharge, velocity, and water-surface elevation data were collected at all four study reaches at various discharges to calibrate the models across a range of discharges. Shields values to predict incipient motion of the substrate were computed for the MHs and surrounding reaches using bed-surface sediment data collected during this study and previous studies. The distributions of hydraulic values at the range of simulated discharge scenarios were significantly different among the MHs. Depth values in the MHs ranged from 0.03 to 5.7 meters, with parts remaining dry at some lower flow scenarios (for example, 90- and 50-percent flow exceedance). MH velocities and bed shear stresses (shear stresses) reached 3.1 meters per second and 31 newtons per square meter, respectively. Through the range of simulated discharges, velocity and shear stress within the MHs were limited by reach-scale hydraulic behavior. Our calculations predicted sand mobility within at least 50 percent of the wetted area of all four MHs for discharges from the 50-percent exceedance flow to the approximate bankfull discharge, whereas 50th-percentile (median) particle size fraction mobility was only predicted within a small area of one of the MHs at the 2-year peak discharge. These results indicate that finer size fractions are mobile within the four MHs, but the larger framework grains of the substrate are predominantly stable at the most frequent discharges. Our results indicate that suitable mussel habitat on the Big River cannot be identified within a narrow range of velocities, depths, and shear stresses. However, the consistent patterns of sediment mobility and the slow increase of hydraulic forces with increasing discharge within all the MHs indicate that flushing flows at low discharges and coarse sediment stability at higher discharges are important for habitat suitability in the Big River. These patterns of sediment mobility are comparable among the robust and depauperate MHs, indicating that the depauperate beds are likely not impaired by bed instability or siltation. Coarse sediment stability up to bankfull discharges further indicates that bed instability is not widespread in these modeled reaches and is likely not related to the spatial distribution of mussels in these locations.

Missouri↗