Search USGS⌕ Search

SEARCH · Search USGS

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 325 records · Page 18Linked to original sources
USGS1997

Analysis of bottom sediment to estimate nonpoint-source phosphorus loads for 1981-96 in Hillsdale Lake, northeast Kansas

Bottom sediment in Hillsdale Lake, northeast Kansas, was analyzed as a means of estimating the annual load of total phosphorus deposited in the lake from nonpoint sources. Topographic, bathymetric, and sediment-core data were used to estimate the total mass of phosphorus in the lake-bottom sediment. Available streamflow and water-quality data were used to compute the mean annual mass of phosphorus (dissolved plus suspended) exiting the lake as well as the mean annual load of phosphorus added to the lake from point sources. A simple mass balance then was used to compute the mean annual load of phosphorus from nonpoint sources. Mean annual sediment deposition from 1981 through 1996 was estimated to be 265 million pounds (120 million kilograms). The total mass of phosphorus in the lake-bottom sediment was estimated to be 924,000 kilograms, with a mean annual load of 62,000 kilograms. The mean annual mass of phosphorus exiting in the lake out-flow was estimated to be about 8,000 kilograms. The mean annual loads of phosphorus added to the lake from point and nonpoint sources were estimated to be 5,000 and 65,000 kilograms, respectively. Thus, the contribution to the total mean annual phosphorus load in Hillsdale Lake from point sources is about 7 percent and from nonpoint sources, about 93 percent.

Kansas↗
USGS2003

Surface-water hydrology of the Gulf Intracoastal Waterway in South-Central Louisiana, 1996-99

The flow of freshwater and suspended sediment from the Lower Atchafalaya River (LAR) and Wax Lake Outlet (WLO) into and along the Gulf Intracoastal Waterway (GIWW) and selected adjacent surface-water bodies between Cypremort and Larose in south-central Louisiana, from October 1996 to December 1999, was characterized using instantaneous and computed continuous discharge measurements and measurements of suspended- sediment concentrations. The GIWW parallels the entire Louisiana coast near the wetland/ upland interface. Following natural hydraulic gradients, the GIWW captures water and sediment from the southward flowing LAR and the WLO where it crosses those waterways, and distributes this freshwater and sediment to points east and west. East of Morgan City, La., an average of 12,200 ft3/s (cubic feet per second) of water flowed from the LAR into the Avoca Island Cutoff Channel. The LAR was the primary source of water to the GIWW east of Morgan City. Drainage from the Verret Subbasin through Bayou Boeuf contributed an average of 1,000 ft3/s to the eastward flow in the GIWW. Eastward flow in the GIWW near Bay Wallace east of Morgan City and to the west of the Houma Navigation Canal (HNC) at Houma, La., averaged about 5,700 ft3/s. Average flow in the GIWW east of the HNC at Houma was 2,610 ft3/s to the east, and 2,200 ft3/s east of Bayou Lafourche at Larose, also to the east. Measured discharge in the GIWW was always to the west between the LAR and WLO. Water entered this stretch of the GIWW from the LAR. The WLO was the primary source of water to the GIWW west of WLO. Discharge in the GIWW averaged 9,460 ft3/s west of WLO south of Calumet and 8,230 ft3/s east of Jaws Bay west of Franklin. Average discharge in the GIWW west of Jaws Bay near Cypremort was 3,310 ft3/s and at Cypremort was 1,350 ft3/s. Average discharge was to the west at all four locations, but discharge as high as 2,830 ft3/s was measured flowing eastward toward Jaws Bay in the GIWW at Cypremort. In bayous and canals in most of coastal Louisiana, including the GIWW, stage narrowly fluctuates around the Gulf of Mexico level. Where the GIWW crosses the LAR and WLO, stage can reach 3 ft (feet) or more above the North American Vertical Datum of 1988 (NAVD88). Flow in the GIWW results from these differences in stage. When the LAR at Morgan City reached 3 to 4 ft above NAVD88, flow in the GIWW became more predictable. Discharge at most sites between the HNC and Jaws Bay increased in varying amounts as stage of the LAR at Morgan City increased beyond 3 ft above NAVD88. At sites in the GIWW east of HNC, discharge did not increase predictably. For all measurements made when the LAR at Morgan City was 3 ft or more above NAVD88, average discharge was about 3,100 ft3/s in the GIWW east of HNC at Houma and 2,880 ft3/s east of Bayou Lafourche at Larose. The LAR at Morgan City is 3 ft or more above NAVD88 for about 7 months in a normal year. When the LAR at Morgan City was less than 3 ft above NAVD88, water in the GIWW flowed along the prevailing water-level gradients, to the east between Bay Wallace and the HNC and to the west between WLO and Jaws Bay. However, local runoff and drainage from areas adjacent to the GIWW became more significant in maintaining flow at low LAR stage. Discharge was consistently higher in the GIWW west of the HNC at Houma than farther west near Bay Wallace east of Morgan City, when the LAR at Morgan City was less than 3 ft above NAVD88. Westward flow in the GIWW between Houma and Morgan City was observed near Bay Wallace east of Morgan City but was never observed west of the HNC at Houma. Discharge in the GIWW east of Jaws Bay west of Franklin, La., frequently was higher than in the GIWW west of WLO south of Calumet, La., at low LAR stage. East of the LAR, suspended-sediment concentrations averaged about 162 mg/L (milligrams per liter) at the two sites closest to the LAR, Avoca Island Cutoff Channel and Bayou Penchant south of Morgan Ci

Professional Paper↗
USGS2020

Identifying areas of degrading and improving groundwater-quality conditions in the State of California, USA, 1974-2014

Areas of improving and degrading groundwater-quality conditions in the State of California were assessed using spatial weighting of a new metric for scoring wells based on constituent concentrations and the direction and magnitude of a trend slope (Sen). Individual well scores were aggregated across 2135 equal-area grid cells covering the entire groundwater resource used for public supply in the state. Spatial weighting allows results to be aggregated locally (well or grid cell), regionally (groundwater basin), provincially, or statewide. Results differentiate degrading (increasing concentration trends) areas with low to moderate concentrations (unimpaired) from degrading areas with moderate to high concentrations (impaired). Results also differentiate improving areas (decreasing concentration trends) in the same manner. Multi-year to decadal groundwater-quality trends were computed from periodic, inorganic water-quality data for 38 constituents collected between 1974 and 2014 for compliance monitoring of nearly 13,000 public-supply wells (PSWs) in the State of California. Mann-Kendall (MK) rank correlations and Sen’s slope estimator were used to detect statistically significant trends for the entire period of recorded data (long-term trend), for the period since 2000 (recent trend), for different pumping seasons (seasonal trend), and for reversals of trends. Statewide, the most frequently detected trends since 2000 were for nitrate (36%), gross alpha/uranium (10%), arsenic (14%), total dissolved solids (TDS) (23%), and the major ions that contribute to TDS (19–28%). The Transverse and Selected Peninsular Ranges (TSPR) and the San Joaquin Valley (SJV) hydrogeologic provinces had the largest percentage of areas with moderate to high nitrate concentrations and groundwater quality trends. Improving nitrate concentrations in parts of the TSPR is associated with long-term managed aquifer recharge that has replaced historical, agriculturally affected groundwater with low-nitrate recharge in parts of the TSPR. This example suggests that application of dilute, excess surface water to agricultural fields during the winter could improve groundwater-quality in the SJV over the long term.

California↗
USGS2017

Perspectives on chemical oceanography in the 21st century: Participants of the COME ABOARD Meeting examine aspects of the field in the context of 40 years of DISCO

The questions that chemical oceanographers prioritize over the coming decades, and the methods we use to address these questions, will define our field's contribution to 21st century science. In recognition of this, the U.S. National Science Foundation and National Oceanic and Atmospheric Administration galvanized a community effort (the Chemical Oceanography MEeting: A BOttom-up Approach to Research Directions, or COME ABOARD) to synthesize bottom-up perspectives on selected areas of research in Chemical Oceanography. Representing only a small subset of the community, COME ABOARD participants did not attempt to identify targeted research directions for the field. Instead, we focused on how best to foster diverse research in Chemical Oceanography, placing emphasis on the following themes: strengthening our core chemical skillset; expanding our tools through collaboration with chemists, engineers, and computer scientists; considering new roles for large programs; enhancing interface research through interdisciplinary collaboration; and expanding ocean literacy by engaging with the public. For each theme, COME ABOARD participants reflected on the present state of Chemical Oceanography, where the community hopes to go and why, and actionable pathways to get there. A unifying concept among the discussions was that dissimilar funding structures and metrics of success may be required to accommodate the various levels of readiness and stages of knowledge development found throughout our community. In addition to the science, participants of the concurrent Dissertations Symposium in Chemical Oceanography (DISCO) XXV, a meeting of recent and forthcoming Ph.D. graduates in Chemical Oceanography, provided perspectives on how our field could show leadership in addressing long-standing diversity and early-career challenges that are pervasive throughout science. Here we summarize the COME ABOARD Meeting discussions, providing a synthesis of reflections and perspectives on the field.

Marine Chemistry↗
USGS2023

Reproducibility starts at the source: R, Python, and Julia Packages for retrieving USGS hydrologic data

Much of modern science takes place in a computational environment, and, increasingly, that environment is programmed using R, Python, or Julia. Furthermore, most scientific data now live on the cloud, so the first step in many workflows is to query a cloud database and load the response into a computational environment for further analysis. Thus, tools that facilitate programmatic data retrieval represent a critical component in reproducible scientific workflows. Earth science is no different in this regard. To fulfill that basic need, we developed R, Python, and Julia packages providing programmatic access to the U.S. Geological Survey’s National Water Information System database and the multi-agency Water Quality Portal. Together, these packages create a common interface for retrieving hydrologic data in the Jupyter ecosystem, which is widely used in water research, operations, and teaching. Source code, documentation, and tutorials for the packages are available on GitHub. Users can go there to learn, raise issues, or contribute improvements within a single platform, which helps foster better engagement and collaboration between data providers and their users.

Water↗
USGS2017

Climate change-driven cliff and beach evolution at decadal to centennial time scales

Here we develop a computationally efficient method that evolves cross-shore profiles of sand beaches with or without cliffs along natural and urban coastal environments and across expansive geographic areas at decadal to centennial time-scales driven by 21st century climate change projections. The model requires projected sea level rise rates, extrema of nearshore wave conditions, bluff recession and shoreline change rates, and cross-shore profiles representing present-day conditions. The model is applied to the ~470-km long coast of the Southern California Bight, USA, using recently available projected nearshore waves and bluff recession and shoreline change rates. The results indicate that eroded cliff material, from unarmored cliffs, contribute 11% to 26% to the total sediment budget. Historical beach nourishment rates will need to increase by more than 30% for a 0.25 m sea level rise (~2044) and by at least 75% by the year 2100 for a 1 m sea level rise, if evolution of the shoreline is to keep pace with rising sea levels.

Conference Paper↗
USGS2008

Hydrological response to timber harvest in northern Idaho: Implications for channel scour and persistence of salmonids

The potential for forest harvest to increase snowmelt rates in maritime snow climates is well recognized. However, questions still exist about the magnitude of peak flow increases in basins larger than 10 km2 and the geomorphic and biological consequences of these changes. In this study, we used observations from two nearly adjacent small basins (13 and 30 km2) in the Coeur d'Alene River basin, one with recent, relatively extensive, timber harvest, and the other with little disturbance in the last 50 years to explore changes in peak flows due to timber harvest and their potential effects on fish. Peak discharge was computed for a specitic rain-on-snow event using a series of physical models that linked predicted values of snowmelt input to a runoff-routing model. Predictions indicate that timber harvest caused a 25% increase in the peak flow of the modelled event and increased the frequency of events of this magnitude from a 9-year recurrence interval to a 3-6-year event. These changes in hydrologic regime, with larger discharges at shorter recurrence intervals, are predicted to increase the depth and frequency of streambed scour, causing up to 15% added mortality of bull trout (Salvelinus confluentus) embryos. Mortality from increased scour, although not catastrophic, may have contributed to the extirpation of this species from the Coeur d'Alene basin, given the widespread timber harvest that occurred in this region. Copyright ?? 2008 John Wiley & Sons, Ltd.

Hydrological Processes↗
USGS2005

Use of the Biotic Ligand Model to predict metal toxicity to aquatic biota in areas of differing geology

This work evaluates the use of the biotic ligand model (BLM), an aquatic toxicity model, to predict toxic effects of metals on aquatic biota in areas underlain by different rock types. The chemical composition of water, soil, and sediment is largely derived from the composition of the underlying rock. Geologic source materials control key attributes of water chemistry that affect metal toxicity to aquatic biota, including: 1) potentially toxic elements, 2) alkalinity, 3) total dissolved solids, and 4) soluble major elements, such as Ca and Mg, which contribute to water hardness. Miller (2002) compiled chemical data for water samples collected in watersheds underlain by ten different rock types, and in a mineralized area in western Colorado. He found that each rock type has a unique range of water chemistry. In this study, the ten rock types were grouped into two general categories, igneous and sedimentary. Water collected in watersheds underlain by sedimentary rock has higher mean pH, alkalinity, and calcium concentrations than water collected in watersheds underlain by igneous rock. Water collected in the mineralized area had elevated concentrations of calcium and sulfate in addition to other chemical constituents. Miller's water-chemistry data were used in the BLM (computer program) to determine copper and zinc toxicity to Daphnia magna. Modeling results show that waters from watersheds underlain by different rock types have characteristic ranges of predicted LC 50 values (a measurement of aquatic toxicity) for copper and zinc, with watersheds underlain by igneous rock having lower predicted LC 50 values than watersheds underlain by sedimentary rock. Lower predicted LC 50 values suggest that aquatic biota in watersheds underlain by igneous rock may be more vulnerable to copper and zinc inputs than aquatic biota in watersheds underlain by sedimentary rock. For both copper and zinc, there is a trend of increasing predicted LC 50 values with increasing dissolved organic carbon (DOC) concentrations. Predicted copper LC 50 values are extremely sensitive to DOC concentrations, whereas alkalinity appears to have an influence on zinc toxicity at alkalinities in excess of about 100 mg/L CaCO 3 . These findings show promise for coupling the BLM (computer program) with measured water-chemistry data to predict metal toxicity to aquatic biota in different geologic settings and under different scenarios. This approach may ultimately be a useful tool for mine-site planning, mitigation and remediation strategies, and ecological risk assessment.

Conference Paper↗
USGS2012

Influence of welded boundaries in anelastic media on energy flow, and characteristics of P, S-I, and S-II waves: Observational evidence for inhomogeneous body waves in low-loss solids

A general computer code, developed to calculate anelastic reflection-refraction coefficients, energy flow, and the physical characteristics for general P , S -I, and S -II waves, quantitatively describes physical characteristics for wave fields in anelastic media that do not exist in elastic media. Consideration of wave fields incident on boundaries between anelastic media shows that scattered wave fields experience reductions in phase and energy speeds, increases in maximum attenuation and Q −1 , and directions of maximum energy flow distinct from phase propagation. Each of these changes in physical characteristics are shown to vary with angle of incidence. Finite relaxation times for anelastic media result in energy flow due to interaction of superimposed radiation fields and contribute to energy flow across anelastic boundaries for all angles of incidence. Agreement of theoretical and numerical results with laboratory measurements argues for the validity of the theoretical and numerical formulations incorporating inhomogeneous wave fields. The agreement attests to the applicability of the model and helps confirm the existence of inhomogeneous body waves and their associated set of distinct physical characteristics in the earth. The existence of such body waves in layered, low-loss anelastic solids implies the need to reformulate some seismological models of the earth. The exact anelastic formulation for a liquid-solid interface with no low-loss approximations predicts the existence of a range of angles of incidence or an anelastic Rayleigh window, through which significant amounts of energy are transmitted across the boundary. The window accounts for the discrepancy apparent between measured reflection data presented in early textbooks and predictions based on classical elasticity theory. Characteristics of the anelastic Rayleigh window are expected to be evident in certain sets of wide-angle, ocean-bottom reflection data and to be useful in estimating Q −1 for some ocean bottom reflectors.

Journal of Geophysical Research B: Solid Earth↗
USGS1975

The Black Hills-Rapid City flood of June 9-10, 1972: A description of the storm and flood

On June 9, 1972, an almost stationary group of thunderstorms formed over the eastern Black Hills of South Dakota near Rapid City and produced record amounts of rainfall and flood discharges. Nearly 15 inches of rain fell in about 6 hours near Nemo, S. Dak., and more than 10 inches of rain fell over a 60-square-mile area. The resulting floods were the highest ever recorded in South Dakota. At least 18 of the 27 streams where peak flows were computed experienced flows that exceeded the 50-year flood. At least 237 people died in the Black Hills flood and 8 people were still listed as missing 6 months after the flood. Another 3,057 people were injured, and total damage is estimated to have exceeded $160 million. As documented by radar and' satellite images and by radiosonde observations, as well as by much other weather data, the unusual and excessive rainfall was partly the result of a strong, low-level, easterly airflow that forced moist air upslope over the Black Hills. This sustained orographic effect helped the air to rise, cool, become very unstable, and release its moisture in repeating -thunderstorms. Another important contributing factor was the unusually light winds at higher levels. Those light winds did not disperse the moist air or move the thunderstorms along to prevent the extreme concentration of rainfall. Instead, these conditions allowed heavy rainfall to remain almost stationary along the eastern slopes of the Black Hills. Rainfall data from 24 observing stations and from over 200 miscellaneous sites show that all reported rains greater than 4 inches fell on the eastern slopes of the Black Hills. At least 6 inches of rain fell over a 300-square-mile area that lies primarily between 4,000 and 5,000 feet elevation. Almost all the flood peaks occurred between 2230 MDT on June 9 and 0100 MDT on June 10, 1972, in a flood belt about 40 miles long and 20 miles wide along the eastern slopes of the Black Hills. This belt extended from Sturgis, S. Dak., on the north to Hermosa, S. Dak., on the south, with Rapid City near the center. To document the flood, peak discharge determinations were made at 49 sites. Records show that about 13,000 acre-feet of water flowed through Rapid City during the 2 days of flooding. At one point during the night of June 9, the floodwaters rose about 3.5 feet in 15 minutes. Coming off the slopes of the Black Hills, the flood peak traveled the 22 miles between Deer Creek and Rapid City in about 3.5 hours.

South Dakota↗
USGS1991

Simulation of the effects of proposed tide gates on circulation, flushing, and water quality in residential canals, Cape Coral, Florida

Decades of dredging and filling of Florida's low-lying coastal wetlands have produced thousands of miles of residential tidal canals and adjacent waterfront property. Typically, these canals are poorly flushed, and over time, accumulated organic-rich bottom materials, contribute to an increasingly severe degraded water quality. One-dimensional hydrodynamic and constituent-transport models were applied to two dead-end canal systems to determine the effects of canal system interconnection using tide gates on water circulation and constituent flushing. The model simulates existing and possible future circulation and flushing conditions in about 29 miles of the approximately 130 miles of tidally influenced canals in Cape Coral, located on the central west coast of peninsular Florida. Model results indicate that tidal water-level differences between the two canal systems can be converted to kinetic energy, in the form of increased water circulation, but the use of one-way tide gate interconnections. Computations show that construction of from one to four tide gates will cause replacement of a volume of water equivalent to the total volume of canals in both systems in 15 to 9 days, respectively. Because some canals flush faster than others, 47 and 21 percent of the original canal water will remain in both systems 50 days after start of operation of one and four tide gates, respectively. Some of the effects that such increased flushing are expected to have include reduced density stratification and associated dissolved-oxygen depletion in canal bottom waters, increased localized reaeration, and more efficient discharge of stormwater runoff entering the canals.

Florida↗
USGS2012

Statistical averaging of marine magnetic anomalies and the aging of oceanic crust

Visual comparison of Mesozoic and Cenozoic magnetic anomalies in the North Pacific suggests that older anomalies contain less short-wavelength information than younger anomalies in this area. To test this observation, magnetic profiles from the North Pacific are examined from crust of three ages: 0–2.1, 29.3–33.1, and 64.9–70.3 m.y, B.P. For each time period, at least nine profiles were analyzed by (1) calculating the power density spectrum of each profile, (2) averaging the spectra together, and (3) computing a ‘recording filter’ for each time period by assuming a hypothetical seafloor model. The model assumes that the top of the source is acoustic basement, the source thickness is 0.5 km, and the time scale of geomagnetic reversals is according to Ness et al. (1980). The calculated power density spectra of the three recording filters are complex in shape but show an increase of attenuation of short-wavelength information as the crust ages. These results are interpreted using a multilayer model for marine magnetic anomalies in which the upper layer, corresponding to pillow basalt of seismic layer 2A, acts as a source of noise to the magnetic anomalies. As the ocean crust ages, this noisy contribution by the pillow basalts becomes less significant to the anomalies. Consequently, magnetic sources below layer 2A must be faithful recorders of geomagnetic reversals.

Journal of Geophysical Research Solid Earth↗
USGS2005

Sediment-transport investigations of the upper Yellowstone River, Montana, 1999 through 2001: Data collection, analysis, and simulation of sediment transport

The upper Yellowstone River in Montana is an important State and national water resource, providing recreational, agricultural, and commercial benefits. Floods in 1996 and 1997, with recorded peak discharges having recurrence intervals close to 100 years, caused substantial streambank erosion and hill- slope mass wasting. Large quantities of sand-, gravel-, and cobble-sized material entrained by the flood flows became flood-bar deposits, creating a source of sediment available for transport during future floods. The flood damage and resulting sedimentation raised concerns about potential streambank-stabilization projects and how the river and riparian corridor might be managed in the future. The U.S. Geological Survey, in cooperation with the Park Conservation District, the Montana Department of Transportation, and the U.S. Army Corps of Engineers, investigated sediment transport in the upper Yellowstone River near Livingston from 1999 through 2001 as part of a cumulative effects study to provide a scientific basis for future river management decisions. The purpose of this report is to present the results of data collection, analysis, and simulation of sediment transport for the upper Yellowstone River. The study area included a 13.5-mile study reach of the upper Yellowstone River where substantial sediment transport occurred in 1996 and 1997. In this study area, the upper Yellowstone River is a high gradient, coarse-bed stream having a slope of about 0.0028 foot per foot or more than 14 feet per mile. The study area drains about 3,551 square miles, and runoff results primarily from snowmelt during the spring and summer months. As part of sediment-transport investigations, the U.S. Geological Survey surveyed river cross sections, characterized streambed-material particle size using particle counts and sieve analyses, and collected bedload- and suspended-sediment data during three runoff seasons (1999-2001). Data were collected for stream discharges that ranged from 2,220 cubic feet per second (typical of pre- and post-runoff discharge) to 25,100 cubic feet per second (about 125 percent of bankfull discharge). The distribution of streambed-material particle size was determined, and sediment-transport curves for bedload discharge, suspended-sediment discharge, and total-sediment discharge were developed. The threshold values of streamflow and average stream velocity needed for initiation of bedload transport for selected sediment-size classes showed that little to no bedload was transported for an average stream velocity below about 3 feet per second, and the only particle size transported as bedload at that velocity was sand. Over the range of stream discharges sampled and with silt- and finer-sized particles excluded, bedload discharge averaged about 18 percent of the total-sediment discharge, equal to bedload discharge plus suspended-sediment discharge. At the lowest and highest stream discharges sampled, bedload was, respectively, less than about 2 percent and about 30 percent of the total-sediment discharge. Over the range of stream discharges sampled, the sand-sized part of the total suspended-sediment discharge averaged about 48 percent, where the total suspended-sediment discharge included sand-, silt- and finer-sized particles. At the lowest and highest stream discharges sampled, the sand-sized part of the total suspended-sediment discharge was, respectively, less than about 16 percent and about 50 percent of the total suspended-sediment discharge. The sediment-transport curves were compared to curves for selected sites in the western United States having drainage areas ranging from 21 square miles to over 20,000 square miles. Daily sediment loads transported at bankfull discharge were calculated for each site and results were plotted in relation to drainage area. Results based on the 1999-2001 data-collection period indicate that the estimated daily bedload transported at bankfull discharge in the upper Yellowstone River exceeded the envelope line that bounds the upper end of the data for other selected sites in the Northern Rocky Mountains and is similar in magnitude to that for selected sites in Alaska having braided channels and glacial and snowmelt runoff. Similar comparisons for suspended sediment indicate that daily suspended-sediment load at bankfull discharge is relatively high in the upper Yellowstone River, plotting slightly above the envelope line that bounds the upper end of the data for other selected sites in the Northern Rocky Mountains. Sediment data were used to develop individual transport equations for seven size classes of sediment ranging from small cobbles to very fine sand. A step-wise regression procedure relating sediment discharge to important hydraulic variables showed that average stream velocity was the only significant variable at the 95-percent confidence level. Bedload and suspended-sediment data and equations indicate that more sand is transported for a given velocity than any other particle size, and very little sand-size sediment load is transported below an average stream velocity of about 2.5 feet per second. Transport of coarser-sized sediment (limited to bedload) becomes very little for an average velocity less than about 3.5 feet per second. Results for the 1999-2001 data-collection period indicate that sediment transport in the upper Yellowstone River tends to be limited more by the transport capacity of the stream (capacity or transport limited), than to the availability of sediment in the watershed (supply limited). Sediment data collected and analyzed were used to simulate sediment transport in the study reach using the BRIdge Stream Tube model for Alluvial River Simulation, or BRI-STARS computer model. The model was calibrated and verified using selected data from historical runoff periods. Simulated total-sediment loads, on a reach-averaged basis, were in good agreement with the total-sediment loads determined from the transport curve for the 2-year flood hydrograph but were considerably smaller for the total-sediment loads determined from the transport curve for the 50-, 100-, and 500-year flood hydrographs. The differences probably were largely due to the inability of the model to simulate streambank erosion, hillslope mass-wasting, and other channel-widening processes, which had supplied substantial quantities of sediment to the channel during the 1996 and 1997 floods, and probably continued to contribute to the sediment load in the subsequent years (1999-2001) when the data were collected. Furthermore, the transport curve was applied beyond the measured data for the highest discharges, and may thus be unreliable. Also, the transport curve derived from only limited data may not apply over the full duration of the hydrograph and sediment might be transported over only a portion of the hydrograph, especially for rivers like the upper Yellowstone where snowmelt runoff predominates. The true sediment discharge is, therefore, unknown and might be closer to the simulated values than to the values calculated from the transport curve.

Montana↗
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↗
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↗