Search USGSSearch

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 199 records · Page 11Linked to original sources

Chapter J: Issues and challenges in the application of geostatistics and spatial-data analysis to the characterization of sand-and-gravel resources

Sand-and-gravel (aggregate) resources are a critical component of the Nation's infrastructure, yet aggregate-mining technologies lag far behind those of metalliferous mining and other sectors. Deposit-evaluation and site-characterization methodologies are antiquated, and few serious studies of the potential applications of spatial-data analysis and geostatistics have been published. However, because of commodity usage and the necessary proximity of a mine to end use, aggregate-resource exploration and evaluation differ fundamentally from comparable activities for metalliferous ores. Acceptable practices, therefore, can reflect this cruder scale. The increasing use of computer technologies is colliding with the need for sand-and-gravel mines to modernize and improve their overall efficiency of exploration, mine planning, scheduling, automation, and other operations. The emergence of megaquarries in the 21st century will also be a contributing factor. Preliminary research into the practical applications of exploratory-data analysis (EDA) have been promising. For example, EDA was used to develop a linear-regression equation to forecast freeze-thaw durability from absorption values for Lower Paleozoic carbonate rocks mined for crushed aggregate from quarries in Oklahoma. Applications of EDA within a spatial context, a method of spatial-data analysis, have also been promising, as with the investigation of undeveloped sand-and-gravel resources in the sedimentary deposits of Pleistocene Lake Bonneville, Utah. Formal geostatistical investigations of sand-and-gravel deposits are quite rare, and the primary focus of those studies that have been completed is on the spatial characterization of deposit thickness and its subsequent effect on ore reserves. A thorough investigation of a gravel deposit in an active aggregate-mining area in central Essex, U.K., emphasized the problems inherent in the geostatistical characterization of particle-size-analysis data. Beyond such factors as common drilling methods jeopardizing the accuracy of the size-distribution curve, the application of formal geostatistical principles has other limitations. Many of the variables used in evaluating gravel deposits, including such sedimentologic parameters as sorting and such United Soil Classification System parameters as gradation coefficient, are nonadditive. Also, uniform sampling methods, such as drilling, are relatively uncommon, and sampling is generally accomplished by a combination of boreholes, water-well logs, test pits, trenches, stratigraphic columns from exposures, and, possibly, some geophysical cross sections. When evaluated in consideration of the fact that uniform mining blocks are also uncommon in practice, subsequent complexities in establishment of the volume/variance relation are inevitable. Several approaches exist to confront the limitations of geostatistical methods in evaluating sand-and-gravel deposits. Initially, we must acknowledge the practical requirements of the aggregate industry, as well as the limitations of the data collected by that industry, as a function of what the industry requires at the practical level, and consider that broader acceptance of formal geostatistics may require modifications of typical exploration and sampling protocols. Future investigations should utilize data from the full spectrum of sand-and-gravel deposits (flood plain, glacial, catastrophic flood, and marine), integrate such other disci plines as sedimentology and geophysics into the research, develop commodity-specific approaches to nonadditive variables, and include the results of comparative drilling.

Bulletin

Suspended-sediment concentrations, bedload, particle sizes, surrogate measurements, and annual sediment loads for selected sites in the lower Minnesota River Basin, water years 2011 through 2016

Accurate measurements of fluvial sediment are important for assessing stream ecological health, calculating flood levels, computing sediment budgets, and managing and protecting water resources. Sediment-enriched rivers in Minnesota are a concern among Federal, State, and local governments because turbidity and sediment-laden waters are the leading impairments and affect more than 6,000 miles of rivers in Minnesota. The suspended sediment in the lower Minnesota River is deleterious, contributing about 75 to 90 percent of the suspended sediment being deposited into Lake Pepin. The Saint Paul District of the U.S. Army Corps of Engineers and the Lower Minnesota River Watershed District collaborate to maintain a navigation channel on the lower 14.7 miles of the Minnesota River through scheduled dredging operations. The Minnesota Pollution Control Agency has adopted a sediment-reduction strategy to reduce sediment in the Minnesota River by 90 percent by 2040. The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, the Minnesota Pollution Control Agency, and the Lower Minnesota River Watershed District, collected suspended-sediment, bedload, and particle-size samples at five sites in the lower Minnesota River Basin during water years 2011 through 2014 and surrogate measurements of acoustic backscatter at one of these sites on the lower Minnesota River during water years 2012 through 2016 to quantify sediment loads and improve understanding of sediment-transport relations. Annual sediment loads were computed for calendar years 2011 through 2014. Data collected from water years 2011 through 2014 indicated that two tributaries, Le Sueur River and High Island Creek, had the highest sediment yield and concentrations of suspended sediment. These tributaries also had greater stream gradients than the sites on the Minnesota River. Suspended fines were greater than suspended sand at all sites in the study area. The range of median particle sizes matched the range for stream gradients from greatest to smallest. Bedload ranged from 3 to 20 percent of the total load at the Le Sueur River, Minnesota River at Mankato, and High Island Creek and was less than 1 percent of the total load at the Minnesota River near Jordan and at Fort Snelling State Park. The reach of the Minnesota River between Mankato and Jordan is a major source of sediment, with the sediment yield at Jordan being two and a half times greater than at Mankato. Between Jordan and Fort Snelling, the sediment yield decreases substantially, which indicates that the Minnesota River in this reach is a sink for sediment. Surrogate measurements (acoustic backscatter) collected with suspended-sediment concentration data from water years 2012 through 2016 from the Minnesota River at Fort Snelling State Park indicated strong relations between the acoustic backscatter and suspended-sediment concentrations. These results point to the dynamic nature of sediment aggradation, degradation, and transport in the Minnesota River Basin. The analyses described in this report will improve the understanding of sediment-transport relations and sediment budgets in the Minnesota River Basin.

Minnesota

Source contributions to suspended sediment and particulate selenium export from the Loutsenhizer Arroyo and Sunflower Drain watersheds in Colorado

Selenium in aquatic ecosystems of the lower Gunnison River Basin in Colorado is affecting the recovery of populations of endangered, native fish species. Dietary exposure is the primary pathway for bioaccumulation of selenium in fish, and particulate selenium can be consumed directly by fish or by the invertebrates on which fish feed. Although selenium can be incorporated into particulate matter via biogeochemical processes, particulate selenium can also enter aquatic ecosystems of the lower Gunnison River Basin from sediments derived from the selenium-rich Mancos Shale. The U.S. Geological Survey, in cooperation with the Colorado Water Conservation Board, conducted this study during 2018–19 to identify sources of selenium-rich suspended sediments from two watersheds underlain by Mancos Shale: Loutsenhizer Arroyo and Sunflower Drain, which is a locally known agricultural drainage near the municipality of Delta, Colorado. A multipronged approach (fieldwork, laboratory work, and computer modeling) referred to as “sediment fingerprinting” was used to evaluate sources of suspended sediments in the streams flowing out of the two studied watersheds. Four potential source types for suspended sediments were identified and sampled (using soil plugs) within the watersheds: rangelands, agricultural fields, arroyo walls, and streambanks. The sediment fingerprinting approach used elemental concentrations and naturally occurring fallout radionuclides as tracers to apportion percent contributions from the four source types of suspended sediments found in streamflow from both watersheds. To determine the dominant sources of suspended sediment in streamflow from both watersheds, a mathematical “unmixing” model was used. Unmixing models apportion source percentages to samples of material in which those sources are mixed. These models used elemental and isotopic data in the suspended sediments to unmix them into proportional contributions from source types. The results indicated that arroyo walls and streambanks generally dominated as sources of the suspended sediment. Arroyo walls and streambanks were channel-adjacent sources, with sediments mobilized by water flowing within the stream channel. These sources accounted for greater than 50 percent of suspended sediment in all but one sample and accounted for 100 percent of suspended sediment in 5 of the 11 samples collected. Rangeland and agricultural field sources were located in uplands outside of stream channels and were detected more often during the non-irrigation season. Rangeland and agricultural field sources each were found in 5 of the 11 samples collected. Concentrations of selenium in sediment-source samples were comparatively greater in streambanks and lower in rangelands, with agricultural fields and arroyo walls being intermediate. As a result, source apportionments for particulate selenium skewed towards sources adjacent to stream channels more than for suspended sediments. Water imports for irrigation have changed the hydrology of the watersheds, and a notable fraction of imported water passes through the watersheds rapidly. The rapid flowthrough water during the irrigation season likely contributes heavily to sediment erosion and transport in Loutsenhizer Arroyo and Sunflower Drain, particularly from channel-adjacent sources of sediment. Decreases in irrigation season streamflow, at least in Loutsenhizer Arroyo, may have decreased sediment erosion and transport during the 2018–20 irrigation seasons compared to the 2015–17 seasons.

Colorado

Geochemical evolution of a high arsenic, alkaline pit-lake in the Mother Lode Gold District, California

The Harvard orebody at the Jamestown gold mine, located along the Melones fault zone in the southern Mother Lode gold district, California, was mined in an open-pit operation from 1987 to 1994. Dewatering during mining produced a hydrologic cone of depression; recovery toward the premining ground-water configuration produced a monomictic pit lake with alkaline Ca-Mg-HCO 3 -SO 4 –type pit water, concentrations of As up to 1,200 μ g/L, and total dissolved solids (TDS) up to 2,000 mg/L. In this study, pit-wall rocks were mapped and chemically analyzed to provide a context for evaluating observed variability in the composition of the pit-lake waters in relationship to seasonal weather patterns. An integrated hydrogeochemical model of pit-lake evolution based on observations of pit-lake volume, water composition (samples collected between 1998–2000, 2004), and processes occurring on pit walls was developed in three stages using the computer code PHREEQC. Stage 1 takes account of seasonally variable water fluxes from precipitation, evaporation, springs, and ground water, as well as lake stratification and mixing processes. Stage 2 adds CO 2 fluxes and wall-rock interactions, and stage 3 assesses the predictive capability of the model. Two major geologic units in fault contact comprise the pit walls. The hanging wall is composed of interlayered slate, metavolcanic and metavolcaniclastic rocks, and schists; the footwall rocks are chlorite-actinolite and talc-tremolite schists generated by metasomatism of greenschist-facies mafic and ultramafic igneous rocks. Alteration in the ore zone provides evidence for mineralizing fluids that introduced CO 2 , S, and K 2 O, and redistributed SiO 2 . Arsenian pyrite associated with the alteration weathers to produce goethite and jarosite on pit walls and in joints, as well as copiapite and hexahydrite efflorescences that accumulate on wall-rock faces during dry California summers. All of these pyrite weathering products incorporate arsenic at concentrations from <100 up to 1,200 ppm. In the pit lake, pH and TDS reach seasonal highs in the summer epilimnion; pH is lowest in the summer hypolimnion. Arsenic and bicarbonate covary in the hypolimnion, rising as stratification proceeds and declining during winter rains. The computational model suggests that water fluxes alone do not account for this seasonal variability. Loss of CO 2 to the atmosphere, interaction with pit walls including washoff of efflorescent salts during the first flush and seasonal rainfall, and arsenic sorption appear to contribute to the observed pit-lake characteristics.

California

Lake Michigan Diversion Accounting land cover change estimation by use of the National Land Cover Dataset and raingage network partitioning analysis

The U.S. Army Corps of Engineers (USACE), Chicago District, is responsible for monitoring and computation of the quantity of Lake Michigan water diverted by the State of Illinois. As part of this effort, the USACE uses the Hydrological Simulation Program&ndash;FORTRAN (HSPF) with measured meteorological data inputs to estimate runoff from the Lake Michigan diversion special contributing areas (SCAs), the North Branch Chicago River above Niles and the Little Calumet River above South Holland gaged basins, and the Lower Des Plaines and the Calumet ungaged that historically drained to Lake Michigan. These simulated runoffs are used for estimating the total runoff component from the diverted Lake Michigan watershed, which is accountable to the total diversion by the State of Illinois. The runoff is simulated from three interpreted land cover types in the HSPF models: impervious, grass, and forest. The three land cover data types currently in use were derived from aerial photographs acquired in the early 1990s. This study used the National Land Cover Dataset (NLCD) and developed an automated process for determining the area of the three land cover types, thereby allowing faster updating of future models, and for evaluating land cover changes by use of historical NLCD datasets. The study also carried out a raingage partitioning analysis so that the segmentation of land cover and rainfall in each modeled unit is directly applicable to the HSPF modeling. Historical and existing impervious, grass, and forest land acreages partitioned by percentages covered by two sets of raingages for the Lake Michigan diversion SCAs, gaged basins, and ungaged basins are presented.

Illinois

Water resources data for California, water year 1976; Volume 1: Colorado River basin, southern Great Basin from Mexican border to Mono Lake basin, and Pacific Slope basins from Tijuana River to Santa Maria River

Water-resources data for the 1976 water year for California consist of records of stage, discharge, and water quality of streams; stage, contents, and water quality of lakes and reservoirs; records of water levels in selected observation wells; and selected chemical analyses of ground water. Records for a few pertinent streamflow and water-quality stations in bordering States are also included. The records were collected and computed by the Water Resources Division of the U.S. Geological Survey under the direction of Lee R. Peterson, district chief; Winchell Smith, assistant district chief for hydrologic data; and Leonard N. Jorgensen, chief of the basic-data section. These data, a contribution to the National Water Data System, were collected by the Geological Survey and cooperating local, State, and Federal agencies in California.

California

Water resources data for California, water year 1976; Volume 2: Pacific Slope basins from Arroyo Grande to Oregon State line except Central Valley

Water-resources data for the 1976 water year for California consist of records of stage, discharge, and water quality of streams; stage, contents, and water quality of lakes and reservoirs; records of water levels in selected observation wells; and selected chemical analyses of ground water. Records for a few pertinent streamflow and water-quality stations in bordering States are also included. The records were collected and computed by the Water Resources Division of the U.S. Geological Survey under the direction of Lee R. Peterson, district chief; Winchell Smith, assistant district chief for hydrologic data; and Leonard N. Jorgensen, chief of the basic-data section. These data, a contribution to the National Water Data System, were collected by the Geological Survey and cooperating local, State, and Federal agencies in California.

California

Water resources data for California, water year 1976; Volume 3, Southern Central Valley basins and the Great Basin from Walker River to Truckee River

Water-resources data for the 1976 water year for California consist of records of stage, discharge, and water quality of streams; stage, contents, and water quality of lakes and reservoirs; records of water levels in selected observation wells; and selected chemical analyses of ground water. Records for a few pertinent streamflow and water-quality stations in bordering States are also included. The records were collected and computed by the Water Resources Division of the U.S. Geological Survey under the direction of Lee R. Peterson, district chief; Winchell Smith, assistant district chief for hydrologic data; and Leonard N. Jorgensen, chief of the basic-data section. These data, a contribution to the National Water Data System, were collected by the Geological Survey and cooperating local, State, and Federal agencies in California.

California

Simulation of the shallow hydrologic system in the vicinity of Middle Genesee Lake, Wisconsin, using analytic elements and parameter estimation

Middle Genesee Lake is a ground-water flow-through lake located in a developing area in southeastern Wisconsin. Because the lake is in good connection with the shallow ground-water system, hydrologic stresses to the shallow ground-water system could adversely affect the lake system. In order to assess the effects of potential stresses on the lake, a study was completed by the U.S. Geological Survey, in cooperation with the Middle Genesee Lake Management District. The objective of the study was to identify areas that contribute ground water to the lake and estimate the hydrologic budget of the lake and hydraulic parameters affecting ground-water flow. A two-dimensional, steady-state analytic element model of the lake and surrounding area was developed using the computer code GFLOW. A parameter estimation model, UCODE, was used to optimize the calibration to measured water levels and streamflow. The calibrated model was used to evaluate the effect of three hypothetical stress scenarios on the stage of Middle Genesee Lake; the simulations were linked to UCODE, which formally incorporated parameter uncertainty into 95-percent confidence intervals around the simulated value. The scenarios included: (1) pumping from upgradient irrigation wells, (2) pumping from Lower Genesee Lake to lower lake levels, and (3) reduction in recharge resulting from development. The results of the simulations demonstrated that lake levels could be affected by hydrologic stresses in the shallow hydrologic system, with effects ranging from a 2.7 feet decline in lake stage resulting from pumping in Lower Genesee Lake to a 0.1 feet decline in lake stage from development in part of the upgradient recharge area. The range of lake stage decline increased when parameter uncertainty was included, from a decline of 3.1 feet for pumping from Lower Genesee Lake to no reduction in lake stage for the development in the recharge area. Whereas these simulated effects are within the natural variation in lake stage, they represent a systematic reduction of ground-water flow to the lake. Therefore, these hypothetical stresses are expected to establish a new, lower, baseline lake stage over which the natural variation due to climatic effects are added and subtracted.

Wisconsin

Occurrence of phosphorus, other nutrients, and triazine herbicides in water from the Hillsdale Lake basin, Northeast Kansas, May 1994 through May 1995

An investigation of the occurrence of phosporus, other nutrients, and triazine herbicides in water samples from the Hillsdale Lake Basin in northeast Kansas was conducted from May 1994 through May 1995. Point-source and nonpoint-source contributions of these water-quality constituents were estimated by conducting synoptic sampling at 48 sites in the basin during five periods of low- flow conditions. Samples were collected for the determination of nutrients, including total phosphorus as phosphorus, dissolved orthophosphate as phosphorus, total nitrite plus nitrate as nitrogen, and total ammonia plus organic nitrogen as nitrogen, and for selected triazine herbicides. On the basis of criteria developed by the Kansas Department of Health and Environment, the Hillsdale Water-Quality Protection Project established a goal to maintain water quality in the tributaries of the Hillsdale Lake Basin at a mean annual low-flow total phosphorus concentration of 0.05 mg/L (milligrams per liter). The mean low- flow total phosphorus concentration of water samples collected in the Big Bull Creek (which includes drainage from Martin Creek), Rock Creek, Little Bull Creek, Wade Branch, and Smith Branch subbasins during low-flow conditions ranged from 0.05 to 4.9 mg/L during this study. Of the 44 sites sampled during low flow, 95 percent had low-flow total phosphorus concentrations larger than the 0.05-mg/L criterion. Discharges from wastewater- treatment plants located in Big Bull Creek and Martin Creek subbasins and the Little Bull Creek subbasin affected nutrient concentrations. Nutrient concentrations in water samples collected from the subbasins not affected by point-source discharges generally were smaller than those in the Big Bull Creek and Little Bull Creek subbasins. Estimated annual low-flow phosphorus loads computed at sampling sites located at the outlet of the subbasins show that the Big Bull Creeksubbasin, which includes drainage from the Martin Creek subbasin, had the largest estimate annual low-flow load, 2,740 kg/yr (kilograms per year).Rock Creek, Little Bull Creek, Wade Branch, and Smith Branch subbasins contributed less annual low-flow phosphorus load, 175, 161, 234, and 22kg/yr, respectively. With the exception of the Smith Branch subbasin, the largest triazine herbicide concentrations occurred in water samples collectedduring May 1994 and May 1995. During May 1994, 10 of 17 sampling sites in the Big Bull Creek and Martin Creek subbasins, 5 of 6 sites in theRock Creek subbasin, and 4 of 10 sites in the Little Bull Creek subbasin had triazine herbicide concentrations in water larger than the U.S.Environmental Protection Agency's Maximum Contaminant Level (MCL), which is an annual mean 3.0 ug/L (micrograms per liter) for atrazine indrinking water. During May 1995, 7 of 19 sites in the Big Bull Creek and Martin Creek subbasins, 5 of 6 sites in the Rock Creek subbasin, 1 of 12 sites in the Little Bull Creek subbasin, and 2 of 4 sites in the Wade Branch subbasin had samples with trazine herbicide concentrations larger than the MCL.Water samples collected in the Rock Creek subbasins had the largest mean triazine herbicide concentrations during May 1994 and May 1995, 6.4 and 4.5 ug/L, respectively.

Water-Resources Investigations Report

Suspended sediment fluxes in a tidal wetland: Measurement, controlling factors, and error analysis

Suspended sediment fluxes to and from tidal wetlands are of increasing concern because of habitat restoration efforts, wetland sustainability as sea level rises, and potential contaminant accumulation. We measured water and sediment fluxes through two channels on Browns Island, at the landward end of San Francisco Bay, United States, to determine the factors that control sediment fluxes on and off the island. In situ instrumentation was deployed between October 10 and November 13, 2003. Acoustic Doppler current profilers and the index velocity method were employed to calculate water fluxes. Suspended sediment concentrations (SSC) were determined with optical sensors and cross-sectional water sampling. All procedures were analyzed for their contribution to total error in the flux measurement. The inability to close the water balance and determination of constituent concentration were identified as the main sources of error; total error was 27% for net sediment flux. The water budget for the island was computed, with an unaccounted input of 0.20 m 3 s −1 (22% of mean inflow), after considering channel flow, change in water storage, evapotranspiration, and precipitation. The net imbalance may be a combination of groundwater seepage, overland flow, and flow through minor channels. Change of island water storage, caused by local variations in water surface elevation, dominated the tidally averaged water flux. These variations were mainly caused by wind and barometric pressure change, which alter regional water levels throughout the Sacramento-San Joaquin River Delta. Peak instantaneous ebb flow was 35% greater than peak flood flow, indicating an ebbdominant system, though dominance varied with the spring-neap cycle. SSC were controlled by wind-wave resuspension adjacent to the island and local tidal currents that mobilized sediment from the channel bed. During neap tides sediment was imported onto the island but during spring tides sediment was exported because the main channel became ebb dominant. Over the 34-d monitoring period 14,000 kg of suspended sediment were imported through the two channels. The water imbalance may affect the sediment balance if the unmeasured water transport pathways are capable of transporting large amounts of sediment. We estimate a maximum of 2,800 kg of sediment may have been exported through unmeasured pathways, giving a minimum ent import of 11,200 kg. Sediment flux measurements provide insight on tidal to fortnightly marsh sedimentation processes, especially in complex systems where sedimentation is spatially and temporally variable.

California

Estimating the “thickness” of the Boulder Batholith, Montana, from heat-flow and heat-productivity data

Estimates of minimum thickness of the Boulder batholith, computed using the linear relation between heat flow and heat productivity and assuming constant heat productivity with depth, are highly nonspecific. They can vary between about 3 and 20 km, depending on values of surface-rock heat productivity and values of assumed contribution of nonbatholith heat sources (such as lower crustal and upper mantle) to the measured surface heat flow used in the calculations. Models involving radiogenic heat sources decreasing with depth in the batholith lead to significantly greater estimates of thickness by as much as a factor of two or more. A reappraisal of data and arguments related to earlier conflicting estimates, based only on one heat-flow determination and within the context of several newly published additional heat-flow measurements, indicates that the previous differences of opinion are negligible and acceptable, in view of the enormous uncertainties inherent in the method of estimation.

Geology

Present and potential sediment yields in the Yampa River Basin, Colorado and Wyoming

Average annual suspended- and total-sediment loads in streamflow were determined by the flow-duration sediment-transport-curve method at 18 sites in the Yampa River basin, Colorado and Wyoming. These computations indicate that about 2.0 million tons of sediment are carried by the Yampa River at Deerlodge Park during an average year. Significant areal differences in the sediment yield from various parts of the basin also were determined. The lower Little Snake River subbasin contributes about 60 percent of the total basin sediment yield, although it represents less than 35 percent of the area and supplies less than 3 percent of the streamflow. In contrast, the upland (eastern) one-third of the basin contributes only about 14 percent of the sediment yield but 76 percent of the streamflow. Projected economic development of the basin, especially surface mining of coal, will impact the physical environment. Depending upon the amount and location of land disturbed, an estimated 10 ,000 to 30,000 tons per year of additional sediment will be contributed to the main-stem Yampa River. (Woodard-USGS)

Water-Resources Investigations Report

Consistency of patterns in concentration‐discharge plots

Concentration‐discharge (c‐Q) plots have been used to infer how flow components such as event water, soil water, and groundwater mix to produce the observed episodic hydrochemical response of small catchments. Because c‐Q plots are based only on observed streamflow and solute concentration, their interpretation requires assumptions about the relative volume, hydrograph timing, and solute concentration of the streamflow end‐members. Evans and Davies [1998] present a taxonomy of c‐Q loops resulting from three‐component conservative mixing. Their analysis, based on a fixed template of end‐member hydrograph volume, timing, and concentration, suggests a unique relationship between c‐Q loop form and the rank order of end‐member concentrations. Many catchments exhibit variability in component contributions to storm flow in response to antecedent conditions or rainfall characteristics, but the effects of such variation on c‐Q relationships have not been studied systematically. Starting with a “baseline” condition similar to that assumed by Evans and Davies [1998] , we use a simple computer model to characterize the variability in c‐Q plot patterns resulting from variation in end‐member volume, timing, and solute concentration. Variability in these three factors can result in more than one c‐Q loop shape for a given rank order of end‐member solute concentrations. The number of resulting hysteresis patterns and their relative frequency depends on the rank order of solute concentrations and on their separation in absolute value. In ambiguous cases the c‐Q loop shape is determined by the relative “prominence” of the event water versus soil water components. This “prominence” is broadly defined as a capacity to influence the total streamflow concentration and may result from a combination of end‐member volume, timing, or concentration. The modeling results indicate that plausible hydrological variability in field situations can confound the interpretation of c‐Q plots, even when fundamental end‐member mixing assumptions are satisfied.

Water Resources Research

Methods for estimating the magnitude and frequency of peak streamflows at ungaged sites in and near the Oklahoma Panhandle

This report presents the results of a cooperative study by the U.S. Geological Survey and the Oklahoma Department of Transportation to estimate the magnitude and frequency of peak streamflows from regional regression equations for ungaged stream sites in and near the Oklahoma Panhandle. These methods relate basin characteristics (physiographic and climatic attributes) to selected peak streamflow frequency statistics with the 50-, 20-, 10-, 4-, 2-, 1-, and 0.2-percent annual exceedance probabilities. These relations were developed based on data from 32 selected streamflow-gaging stations in the Oklahoma Panhandle and in neighboring parts of Colorado, Kansas, New Mexico, and Texas. The basin characteristics for the selected streamflow-gaging stations were determined by using a geographic information system and the Oklahoma StreamStats application. Peak-streamflow frequency statistics were computed from annual peak-streamflow records from the irrigated period of record from water year 1978 through water year 2014. Generalized-least-squares multiple-linear regression analysis was used to formulate regression relations between peak-streamflow frequency statistics and basin characteristics. Contributing drainage area was the only basin characteristic determined to be statistically significant for all percentage of annual exceedance probabilities and was the only basin characteristic used in regional regression equations for estimating peak-streamflow frequency statistics on unregulated streams in and near the Oklahoma Panhandle. The regression model pseudo-coefficient of determination, converted to percent, for the Oklahoma Panhandle regional regression equations ranged from about 38 to 63 percent. The standard errors of prediction and the standard model errors for the Oklahoma Panhandle regional regression equations ranged from about 84 to 148 percent and from about 76 to 138 percent, respectively. These errors were comparable to those reported for regional peak-streamflow frequency regression equations for the High Plains areas of Texas and Colorado. The root mean square errors for the Oklahoma Panhandle regional regression equations (ranging from 3,170 to 92,000 cubic feet per second) were less than the root mean square errors for the Oklahoma statewide regression equations (ranging from 18,900 to 412,000 cubic feet per second); therefore, the Oklahoma Panhandle regional regression equations produce more accurate peak-streamflow statistic estimates for the irrigated period of record in the Oklahoma Panhandle than do the Oklahoma statewide regression equations. The regression equations developed in this report are applicable to streams that are not substantially affected by regulation, impoundment, or surface-water withdrawals. These regression equations are intended for use for stream sites with contributing drainage areas less than or equal to about 2,060 square miles, the maximum value for the independent variable used in the regression analysis.

Colorado, Kansas, New Mexico, Oklahoma, Texas

Factors affecting nutrient trends in major rivers of the Chesapeake Bay Watershed

Trends in nutrient loads and flow-adjusted concentrations in the major rivers entering Chesapeake Bay were computed on the basis of water-quality data collected between 1985 and 1998 at 29 monitoring stations in the Susquehanna, Potomac, James, Rappahannock, York, Patuxent, and Choptank River Basins. Two computer models-the Chesapeake Bay Watershed Model (WSM) and the U.S. Geological Survey's 'Spatially Referenced Regressions on Watershed attributes' (SPARROW) Model-were used to help explain the major factors affecting the trends. Results from WSM simulations provided information on temporal changes in contributions from major nutrient sources, and results from SPARROW model simulations provided spatial detail on the distribution of nutrient yields in these basins. Additional data on nutrient sources, basin characteristics, implementation of management practices, and ground-water inputs to surface water were analyzed to help explain the trends. The major factors affecting the trends were changes in nutrient sources and natural variations in streamflow. The dominant source of nitrogen and phosphorus from 1985 to 1998 in six of the seven tributary basins to Chesapeake Bay was determined to be agriculture. Because of the predominance of agricultural inputs, changes in agricultural nutrient sources such as manure and fertilizer, combined with decreases in agricultural acreage and implementation of best management practices (BMPs), had the greatest impact on the trends in flow-adjusted nutrient concentrations. Urban acreage and population, however, were noted to be increasing throughout the Chesapeake Bay Watershed, and as a result, delivered loads of nutrients from urban areas increased during the study period. Overall, agricultural nutrient management, in combination with load decreases from point sources due to facility upgrades and the phosphate detergent ban, led to downward trends in flow-adjusted nutrient concentrations atmany of the monitoring stations in the watershed. The loads of nutrients, however, were not reduced significantly at most of the monitoring stations. This is due primarily to higher streamflow in the latter years of the monitoring period, which led to higher loading in those years. Results of this study indicate a need for more detailed information on BMP effectiveness under a full range of hydrologic conditions and in different areas of the watershed; an internally consistent fertilizer data set; greater consideration of the effects of watershed processes on nutrient transport; a refinement of current modeling efforts; and an expansion of the non-tidal monitoring network in the Chesapeake Bay Watershed.

Delaware, Maryland, New York, Pennsylvania, Virgin

Streamflow and nutrient data for the Yazoo River below Steele Bayou near Long Lake, Mississippi, 1996-2000

Increased nutrient loading to the Gulf of Mexico from off-continent flux has been identified as contributing to the increase in the areal extent of the low dissolved-oxygen zone that develops annually off the Louisiana and Texas coast. The proximity of the Yazoo River Basin in northwestern Mississippi to the Gulf of Mexico, and the intensive agricultural activities in the basin have led to speculation that the Yazoo River Basin contributes a disproportionate amount of nitrogen and phosphorus to the Mississippi River and ultimately to the Gulf of Mexico. An empirical measurement of the flux of nitrogen and phosphorus from the Yazoo Basin has not been possible due to the hydrology of the lower Yazoo River Basin. Streamflow for the Yazoo River below Steele Bayou is affected by backwater from the Mississippi River. Flow at the gage is non-uniform and varying, with bi-directional and reverse flows possible. Streamflow was computed by using remote sensing and acoustic and conventional discharge and velocity measurement techniques. Streamflow from the Yazoo River for the 1996-2000 period accounted for 2.8 percent of the flow of the Mississippi River for the same period. Water samples from the Yazoo River were collected from February 1996 through December 2000 and were analyzed for total nitrogen, nitrate, total phosphorus, and orthophosphorus as part of the U.S. Geological Survey National Water-Quality Assessment Program. These data were used to compute annual loads of nitrogen and phosphorus discharged from the Yazoo River for the period 1996-2000. Annual loads of nitrogen and phosphorus were calculated by two methods. The first method used multivariate regression and the second method multiplied the mean annual concentration by the total annual flow. Load estimates based on the product of the mean annual concentration and the total annual flow were within the 95 percent confidence interval for the load calculated by multivariate regression in 10 of 20 cases. The Yazoo River loads, compared to average annual loads in the Mississippi River, indicated that the Yazoo River was contributing 1.4 percent of the total nitrogen load, 0.7 percent of the nitrate load, 3.4 percent of the total phosphorus load, and 1.6 percent of the orthophosphorus load during 1996 - 2000. The total nitrogen, nitrate, and orthophosphorus loads in the Yazoo River Basin were less than expected, whereas the total phosphorus load was slightly higher than expected based on discharge.

Water-Resources Investigations Report

Skill assessment of a total water level and coastal change forecast during the landfall of a hurricane

The Total Water Level and Coastal Change Forecast (TWL&CC Forecast) provides coastal communities with 6-day notice of potential elevated water levels and coastal change (i.e., dune erosion, overwash, or inundation) on sandy beaches that threatens safety, infrastructure, or resources. This continuously operating model provides hourly information for select regions along U.S. Gulf of Mexico and Atlantic Ocean coastlines. The objective of this work is to assess the skill of forecasts during a period of elevated water levels along the coasts of North Carolina (NC) and South Carolina, USA caused by Hurricane Isaias in August 2020, using a combination of observations and model hindcasts. Water levels and waves were observed throughout the storm at three locations near Wrightsville Beach, NC, which provided information to assess forecast skill; a wave buoy offshore, a tide gage at a local pier, and a pressure sensor deployed at the pier. In addition to observations, the non-hydrostatic phase-resolving model SWASH (Simulating WAves till SHore) was forced with hourly wave energy spectra derived from a coupled Delft3D-SWAN simulation during the peak of Isaias, to complement observations by computing nearshore wave height and wave-induced setup and runup at the shoreline. During the storm peak, SWASH-simulated water levels at the sensor position were comparable to those at the maximum landward extent (bias = −0.05 m; gain = 0.26; r 2 = 0.99), suggesting that observations at the USGS sensor location were a useful proxy for total water level (TWL; sum of tide, surge and wave runup) at the shoreline that are predicted by the TWL&CC Forecast. The TWL forecast at Wrightsville Beach was consistent with observations from the USGS sensor (bias = −0.38 m and −0.74 m, scatter index = 0.22 and 0.28 for the two forecast model grids considered, respectively; weighted regression considering model uncertainty explained 95 percent of variability in observed TWL). Observed TWL was within the confidence interval of the TWL&CC Forecast for the 5 h at the storm peak. Forecast mean water levels (MWL; sum of tide, surge and wave setup) and tide gage observations were also consistent (bias = 0.07 m and 0.02 m for the forecast model grids; scatter index = 0.46; r 2 = 0.80). Forecast MWL at the storm peak was within 0.06 m of the observed MWL from the tide gage for both sites. In the region where Isaias made landfall, eight additional pressure sensors were compared to the peak TWL forecast (bias = 0.14 m; scatter index = 0.18). Forecast TWL explained 90 percent of observed variability in TWL when considering uncertainty of the forecast with a weighted regression. The results demonstrate that wave-driven water levels contributed a significant portion of the forecast TWL during Isaias (52 percent during the three peak hours of the storm), and that TWL were represented using the forecast model. Mean absolute error of the coastal change forecast and observed overwash is 0.4 and 0.14 for the two forecast model grids considered. The skill demonstrated by this computationally efficient method indicates that the forecasting system can provide fast and reliable predictions of TWL across hundreds of km of coastline at sub-km resolution, days to hours in advance of when storms threaten coastal regions.

Coastal Engineering