SEARCH · USGS Science

Results for “Scientific Monthly”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

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

At least 757 records · Page 42Linked to original sources

USGS2008

Drier forest composition associated with hydrologic change in the Apalachicola River floodplain, Florida

Forests of the Apalachicola River floodplain had shorter flood durations, were drier in composition, and had 17 percent fewer trees in 2004 than in 1976. The change to drier forest composition is expected to continue for at least 80 more years. Floodplain drying was caused by large declines in river levels resulting from erosion of the river channel after 1954 and from decreased flows in spring and summer months since the 1970s. Water-level declines have been greatest at low and medium flows, which are the most common flows (occurring about 80 percent of the time). Water levels have remained relatively unchanged during large floods which continue to occur about three times per decade. A study conducted by the U.S. Geological Survey compared temporal changes in hydrologic conditions, forest composition, forest characteristics, and individual species of trees, as well as estimated the potential for change in composition of floodplain forests in the nontidal reach of the Apalachicola River. The study was conducted with the cooperation of the Florida Department of Environmental Protection and the Northwest Florida Water Management District. Forest composition and field observations from studies conducted in 1976-1984 (termed “1976 data”) were used as baseline data for comparison with data from plots sampled in 2004-2006 (“2004 data”). Flood durations were shorter in all periods subsequent to 1923-1976. The periods of record used to calculate flood durations for forest data were subsets of the complete record available (1923-2004). At sampled plots in all forest types and reaches combined, flood durations changed an average of more than 70 percent toward the baseline flood duration of the next drier forest type. For all forest types, changes in flood durations toward the next drier type were greatest in the upper reach (95.9 percent) and least in the lower reach (42.0 percent). All forests are expected to be 38.2 percent drier in species composition by 2085, the year when the median age of surviving 2004 subcanopy trees will reach the median age (99 years) of the 2004 large canopy trees. The change will be greatest for forests in the upper reach (45.0 percent). Forest composition changes from pre-1954 to 2085 were calculated using Floodplain Indices from 1976 and 2004 tree-size classes and replicate plots. Species composition in high bottomland hardwood forests is expected to continue to change, and some low bottomland hardwood forests are expected to become high bottomland hardwood forests. Organisms associated with floodplain forests will be affected by the changes in tree species, which will alter the timing of leaf-out, fruiting, and leaf-drop, the types of fruit and debris produced, and soil chemistry. Swamps will contain more bottomland hardwood species, but will also have an overall loss of tree density. The density of trees in swamps significantly decreased by 37 percent from 1976 to 2004. Of the estimated 4.3 million (17 percent) fewer trees that existed in the nontidal floodplain in 2004 than in 1976, 3.3 million trees belonged to four swamp species: popash, Ogeechee tupelo, water tupelo, and bald cypress. Water tupelo, the most important tree in the nontidal floodplain in terms of basal area and density, has declined in number of trees by nearly 20 percent since 1976. Ogeechee tupelo, the species valuable to the tupelo honey industry, has declined in number of trees by at least 44 percent. Greater hydrologic variability in recent years may be the reason swamps have had a large decrease in tree density. Drier conditions are detrimental for the growth of swamp species, and periodic large floods kill invading bottomland hardwood trees. The loss of canopy density in swamps may result in the swamp floor being exposed to more light with an increase in the amount of ground cover present, which in turn, would reduce tree replacement. The microclimate of the swamp floor would become warmer due to the decrease in shade and inundation. Soils would become dehydrated more quickly in dry periods and debris would decompose more quickly. A loss of tree density in swamps would lead to a decrease in tree and leaf litter biomass, which would have additional effects on swamp organisms. The loss of litter would result in a loss of substrate for benthic organisms in the floodplain and, ultimately, in the downstream waters of the river and estuary.

Florida↗

USGS2015

Updated numerical model with uncertainty assessment of 1950-56 drought conditions on brackish-water movement within the Edwards aquifer, San Antonio, Texas

In 2010, the U.S. Geological Survey, in cooperation with the San Antonio Water System, began a study to assess the brackish-water movement within the Edwards aquifer (more specifically the potential for brackish-water encroachment into wells near the interface between the freshwater and brackish-water transition zones, referred to in this report as the transition-zone interface) and effects on spring discharge at Comal and San Marcos Springs under drought conditions using a numerical model. The quantitative targets of this study are to predict the effects of higher-than-average groundwater withdrawals from wells and drought-of-record rainfall conditions of 1950–56 on (1) dissolved-solids concentration changes at production wells near the transition-zone interface, (2) total spring discharge at Comal and San Marcos Springs, and (3) the groundwater head (head) at Bexar County index well J-17. The predictions of interest, and the parameters implemented into the model, were evaluated to quantify their uncertainty so the results of the predictions could be presented in terms of a 95-percent credible interval. The model area covers the San Antonio and Barton Springs segments of the Edwards aquifer; the history-matching effort was focused on the San Antonio segment. A previously developed diffuse-flow model of the Edwards aquifer, which forms the basis for the model in this assessment, is primarily based on a conceptualization in which flow in the aquifer is predominately through a network of numerous small fractures and openings. Primary updates to this model include an extension of the active area downdip, a conversion to an 8-layer SEAWAT variable-density flow and transport model to simulate dissolved-solids concentration effects on water density, history matching to 1999–2009 conditions, and parameter estimation in a highly parameterized context using automated methods in PEST (a model-independent Parameter ESTimation code). In addition to the best-fit parameter values derived from history matching, the uncertainty of model parameters was also estimated by using linear uncertainty analysis. Comparison of “prior” (before history matching) and “posterior” (after history matching) variances of parameters indicate that the information within the observation dataset used for history matching informs many parameters. The concentration threshold parameters were well-informed by the observation dataset as their posterior distributions were much narrower than their prior distributions. The transition-zone scaling parameters of hydraulic conductivity, effective porosity, and specific storage were all informed by the observation dataset, as evidenced by the difference between the prior and posterior variances. Saline-zone scaling parameters, alternatively, were not informed by the observation dataset for effective porosity and specific storage. Resulting posterior drier-month, wetter-month, and annual recharge multiplier parameter variances are important to understanding how well recharge is estimated and implemented within the model. The shifts of the posterior distributions left and right indicate that there were zones where less or more water was needed in the model. The widths of the distributions were not decreased substantially, indicating that many of the best-fit recharge parameters are not statistically different from the initial values specified in the history-matching effort. Recharge from rainfall is the driving force behind groundwater flow and heads in the aquifer; therefore, an increase in understanding of this process would benefit model development by potentially decreasing the uncertainty of this parameter. The history-matching effort was most helpful in informing the parameters in the model that control discharge at springs, namely, the spring orifice (drain) altitude and drain conductance parameters for each spring. The uncertainty assessment of the predictive model (a hypothetical recurrence of 1950–56 drought conditions and higher-than-average groundwater withdrawals from wells) provided insights into the potential effects of these conditions on dissolved-solids concentration changes at production wells near the transition-zone interface, discharges at Comal and San Marcos Springs, and heads at Bexar County index well J-17. Results at the 25 production wells near the transition-zone interface indicate that the uncertainty of model input parameters based on expert knowledge yielded an upper bound of the 95-percent credible interval of dissolved-solids concentrations that exceeds the secondary drinking water standards of 1,000 milligrams per liter (mg/L) of the Texas Commission on Environmental Quality (TCEQ) for many wells. However, the history-matching process provided key information to inform prediction-sensitive model parameters and therefore, contributed to a substantial decrease of the upper bound of the 95-percent credible interval to below the secondary drinking water standards. Reductions in dissolved-solids concentration changes were on the order of 400 mg/L to 1,300 mg/L. The reduction in uncertainty in regards to this prediction implies that this prediction of dissolved-solids concentration change can be made with some certainty using this current model and that those parameters that control this prediction are informed by the observation dataset. Even though predictive uncertainty was reduced for this prediction, dissolved-solids concentration changes were still greater than zero, indicating a minimal increase in concentration at these 25 production wells during the 7-year simulation period is likely. However, this minimal concentration increase indicates a small potential for movement of the brackish-water transition zone near these wells during the 7-year simulation period of drought-ofrecord (1950–56) rainfall conditions with higher-than-average groundwater withdrawals by wells. Predictive results of total spring discharge during the 7-year period, as well as head predictions at Bexar County index well J-17, were much different than the dissolved-solids concentration change results at the production wells. These upper bounds are an order of magnitude larger than the actual prediction which implies that (1) the predictions of total spring discharge at Comal and San Marcos Springs and head at Bexar County index well J-17 made with this model are not reliable, and (2) parameters that control these predictions are not informed well by the observation dataset during historymatching, even though the history-matching process yielded parameters to reproduce spring discharges and heads at these locations during the history-matching period. Furthermore, because spring discharges at these two springs and heads at Bexar County index well J-17 represent more of a cumulative effect of upstream conditions over a larger distance (and longer time), many more parameters (with their own uncertainties) are potentially controlling these predictions than the prediction of dissolved-solids concentration change at the prediction wells, and therefore contributing to a large posterior uncertainty.

Texas↗

USGS2015

California State Waters Map Series: Offshore of Refugio Beach, California

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within California’s State Waters. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data, acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. The Offshore of Refugio Beach map lies within the western Santa Barbara Channel region of the Southern California Bight. This geologically complex region forms a major biogeographic transition zone, separating the cold-temperate Oregonian province north of Point Conception from the warm-temperate California province to the south. The map area is in the southern part of the Western Transverse Ranges geologic province, which is north of the California Continental Borderland. Significant clockwise rotation—at least 90°—since the early Miocene has been proposed for the Western Transverse Ranges province, and geodetic studies indicate that the region is presently undergoing north-south shortening. Uplift rates (as much as 0.5 mm/yr) that are based on studies of onland marine terraces provide further evidence of significant shortening. The coastal zone of this map area lies at the steep flank of the Santa Ynez Mountains. The crest of the range, which lies about 8 km from the shoreline (north of the map area), has a maximum elevation of about 780 m. This area is largely open space, partly used for livestock grazing, with no significant towns or population centers. Highway 101 crosses the map area, adjacent to and within a few hundred meters of the shoreline. The most significant developments are the recreational state beaches at El Capitan Beach and Refugio Beach. The beaches are subject to erosion each winter during storm-wave attack, and then they undergo gradual recovery or accretion during the more gentle wave climate of the late spring, summer, and fall months. The Offshore of Refugio Beach map area lies in the west-central part of the Santa Barbara littoral cell, which is characterized by west-to-east transport of sediment from Point Arguello on the northwest to Hueneme and Mugu Canyons on the southeast. Longshore drift rates have been reported to range from about 160,000 to 800,000 tons/yr, averaging 400,000 tons/yr. Sediment supply to the western and central part of the littoral cell, including the map area, is mainly from relatively small coastal watersheds. Within the map area, these coastal watersheds include (from east to west) Cañada del Capitan, Tajiquas Creek, Arroyo Hondo, Cañada del Molino, and several unnamed canyons and creeks. The Santa Ynez and Santa Maria Rivers, the mouths of which are 80 to 120 km northwest of the map area, are not considered significant sediment sources because Point Conception and Point Arguello provide obstacles to downcoast sediment transport and also because much of their sediment load is trapped in dams. The Ventura and Santa Clara Rivers, the mouths of which are about 70 km to the southeast of the map area, are not sediment sources for the map area. The offshore part of the map area consists of relatively flat and shallow continental shelf, which dips gently seaward (about 0.8° to 1.0°) so that water depths at the shelf break, roughly coincident with the California’s State Waters limit, are about 80 to 100 m. This part of the Santa Barbara Channel is relatively well protected from large Pacific swells from the north and northwest by Point Conception and from the south and southwest by offshore islands and banks. The shelf is underlain by variable amounts of upper Quaternary marine and fluvial sediments deposited as sea level fluctuated in the late Pleistocene. In the map area, the shelf break is at depths of about 90 m and lies about 5.6 to 6.4 km offshore. Beyond the shelf break, the slope is steep (as much as about 7°) and unstable. Several submarine landslides, including the large (130 km 2 ) Goleta landslide complex, have been documented offshore of Goleta, a few kilometers east of the map area. This compound slump complex may have been initiated more than 200,000 year ago, but it also includes three recent failures that may have been generated 8,000 to 10,000 years ago. A local, 10-m-high tsunami may have been generated from these failure events. Small folds related to local faulting are superimposed on the homocline that makes up the south flank of the Santa Ynez Mountains. One of these superimposed anticlines hosts the Molino gas field, which was discovered in 1962 and subsequently developed through directional drilling from onshore wells. The Oligocene Sespe and Vaqueros Formations are the reservoirs in the Molino gas field, and the map area includes numerous seafloor hydrocarbon seeps. Seafloor habitats in the broad Santa Barbara Channel region consist of significant amounts of soft, unconsolidated sediment interspersed with isolated areas of rocky habitat that support kelp-forest communities nearshore and rocky-reef communities in deep water. The potential marine benthic habitat types mapped in the Offshore of Refugio Beach map area are directly related to its Quaternary geologic history, geomorphology, and active sedimentary processes. These potential habitats, which lie primarily within the Shelf (continental shelf) but also partly within the Flank (basin flank or continental slope) megahabitats, primarily are composed of soft sediment interrupted by a few carbonate mounds. This homogeneous seafloor of sediment and low-relief bedrock provides promising habitat for groundfish, crabs, shrimp, and other marine benthic organisms.

California↗

USGS2017

Water quality and quantity and simulated surface-water and groundwater flow in the Laurel Hill Creek Basin, southwestern Pennsylvania, 1991–2007

Laurel Hill Creek is considered one of the most pristine waterways in southwestern Pennsylvania and has high recreational value as a high-quality cold-water fishery; however, the upper parts of the basin have documented water-quality impairments. Groundwater and surface water are withdrawn for public water supply and the basin has been identified as a Critical Water Planning Area (CWPA) under the State Water Plan. The U.S. Geological Survey, in cooperation with the Somerset County Conservation District, collected data and developed modeling tools to support the assessment of water-quality and water-quantity issues for a basin designated as a CWPA. Streams, springs, and groundwater wells were sampled for water quality in 2007. Streamflows were measured concurrent with water-quality sampling at main-stem sites on Laurel Hill Creek and tributaries in 2007. Stream temperatures were monitored continuously at five main-stem sites from 2007 to 2010. Water usage in the basin was summarized for 2003 and 2009 and a Water-Analysis Screening Tool (WAST) developed for the Pennsylvania State Water Plan was implemented to determine whether the water use in the basin exceeded the “safe yield” or “ the amount of water that can be withdrawn from a water resource over a period of time without impairing the long-term utility of a water resource .” A groundwater and surface-water flow (GSFLOW) model was developed for Laurel Hill Creek and calibrated to the measured daily streamflow from 1991 to 2007 for the streamflow-gaging station near the outlet of the basin at Ursina, Pa. The CWPA designation requires an assessment of current and future water use. The calibrated GSFLOW model can be used to assess the hydrologic effects of future changes in water use and land use in the basin. Analyses of samples collected for surface-water quality during base-flow conditions indicate that the highest nutrient concentrations in the main stem of Laurel Hill Creek were at sites in the northeastern part of the basin where agricultural activity is prominent. All of the total nitrogen (N) and a majority of the total phosphorus (P) concentrations in the main stem exceeded regional nutrient criteria levels of 0.31 and 0.01 milligrams per liter (mg/L), respectively. The highest total N and total P concentrations in the main stem were 1.42 and 0.06 mg/L, respectively. Tributary sites with the highest nutrient concentrations are in subbasins where treated wastewater is discharged, such as Kooser Run and Lost Creek. The highest total N and total P concentrations in subbasins were 3.45 and 0.11 mg/L, respectively. Dissolved chloride and sodium concentrations were highest in the upper part of the basin downstream from Interstate 76 because of road deicing salts. The mean base-flow concentrations of dissolved chloride and sodium were 117 and 77 mg/L, respectively, in samples from the main stem just below Interstate 76, and the mean concentrations in Clear Run were 210 and 118 mg/L, compared to concentrations less than 15 mg/L in tributaries that were not affected by highway runoff. Water quality in forested tributary subbasins underlain by the Allegheny and Pottsville Formations was influenced by acidic precipitation and, to a lesser extent, the underlying geology as indicated by pH values less than 5.0 and corresponding specific conductance ranging from 26 to 288 microsiemens per centimeter at 25 degrees Celsius for some samples; in contrast, pH values for main stem sites ranged from 6.6 to 8.5. Manganese (Mn) was the only dissolved constituent in the surface-water samples that exceeded the secondary maximum contaminant level (SMCL). More than one-half the samples from the main stem had Mn concentrations exceeding the SMCL level of 50 micrograms per liter (μg/L), whereas only 19 percent of samples from tributaries exceeded the SMCL for Mn. Stream temperatures along the main stem of Laurel Hill Creek became higher moving downstream. During the summer months of June through August, the daily mean temperatures at the five sites exceeded the limit of 18.9 degrees Celsius (°C) for a cold-water fishery. The maximum instantaneous values for each site ranged from 27.2 to 32.8 °C. Water-quality samples collected at groundwater sites (wells and springs) indicate that wells developed within the Mauch Chunk Formation had the best water quality, whereas wells developed within the Allegheny and Pottsville Formations yielded the poorest water quality. Waters from the Mauch Chunk Formation had the highest median pH (7.6) and alkalinity (80 mg/L calcium carbonate) values. The lowest pH and alkalinity median values were in waters from the Allegheny and Pottsville Formations. Groundwater samples collected from wells in the Allegheny and Pottsville Formations also had the highest concentrations of dissolved iron (Fe) and dissolved Mn. Seventy-eight percent of the groundwater samples collected from the Allegheny Formation exceeded the SMCL of 300 μg/L for Fe and 50 μg/L for Mn. Forty-three and 62 percent of the groundwater samples collected from the Pottsville Formation exceeded the SMCL for iron and Mn, respectively. The highest Fe and Mn concentrations for surface waters were measured for tributaries draining the Pottsville Formation. The highest median Fe concentration for tributaries was in samples from streams draining the Allegheny Formation. During base-flow conditions, the streamflow per unit area along the main stem of Laurel Hill Creek was lowest in the upper parts of the basin [farthest upstream site 0.07 cubic foot per second per square mile (ft 3 /s/mi 2 )] and highest (two sites averaging about 0.20 (ft 3 /s/mi 2 ) immediately downstream from Laurel Hill Lake in the center of the basin. Tributaries with the highest streamflow per unit area were those subbasins that drain the western ridge of the Laurel Hill Creek Basin. The mean streamflow per unit area for tributaries draining areas that extend into the western ridge and draining eastern or central sections was 0.24 and 0.05 ft 3 /s/mi 2 , respectively. In general, as the drainage area increased for tributary basins, the streamflow per unit area increased. Criteria established by the Pennsylvania Department of Environmental Protection indicate that the safe yield of water withdrawals from the Laurel Hill Creek Basin is 1.43 million gallons per day (Mgal/d). Water-use data for 2009 indicate that net (water withdrawals subtracted by water discharges) water withdrawals from groundwater and surface-water sources in the basin were approximately 1.93 Mgal/d. Water withdrawals were concentrated in the upper part of the basin with approximately 80 percent of the withdrawals occurring in the upper 36 mi 2 of the basin. Three subbasins—Allen Creek, Kooser Run, and Shafer Run— in the upper part were affected the most by water withdrawals such that safe yields were exceeded by more than 1,000 percent in the first two and more than 500 percent in the other. In the subbasin of Shafer Run, intermittent streamflow characterizes sections that historically have been perennial. The GSFLOW model of the Laurel Hill Creek Basin is a simple one-layer representation of the groundwater flow system. The GSFLOW model was primarily calibrated to reduce the error term associated with base-flow periods. The total amount of observed streamflow at the Laurel Hill Creek at Ursina, Pa. streamflow-gaging station and the simulated streamflow were within 0.1 percent over the entire modeled period; however, annual differences between simulated and observed streamflow showed a range of -27 to 24 percent from 1992 to 2007 with nine of the years having less than a 10-percent difference. The primary source of simulated streamflow in the GSFLOW model was the subsurface (interflow; 62 percent), followed by groundwater (25 percent) and surface runoff (13 percent). Most of the simulated subsurface flow that reached the stream was in the form of slow flow as opposed to preferential (fast) interflow.

Pennsylvania↗

USGS2019

Geologic map of the central-southeast flank of Mauna Loa Volcano, Island of Hawaii, Hawaii

Mauna Loa, the largest volcano on Earth, has erupted 33 times since written descriptions became available in 1832. Some eruptions began with only brief seismic unrest, while others followed several months to a year of increased seismicity. Once underway, its eruptions can produce lava flows that may reach the sea in less than 24 hours, severing roads and utilities. For example, lava flows erupted from the Southwest Rift Zone (SWRZ) in 1950 advanced at an average rate of 9.3 km per hour, and all three lobes reached the ocean within approximately 24 hours (Finch and Macdonald, 1953). Near the eruptive vents, the flows must have traveled even faster. In terms of eruption frequency, pre-eruption warning, and rapid flow emplacement, Mauna Loa poses an enormous volcanic-hazard threat to the Island of Hawai‘i. Volcanic hazards on Mauna Loa may be anticipated, and risk substantially mitigated, by documenting the past activity to refine our knowledge of the hazards and by alerting the public and local government officials of our findings and their implications for hazards assessments and risk. From the geologic record, we may deduce several generalized facts about the geologic history of the Northeast Rift Zone (NERZ). The middle to uppermost segments of the rift zone were more active in the past 4,000 years than the lower portion of the rift zone. This may be due to buttressing of the lower east rift zone by Mauna Kea and Kīlauea volcanoes. The historical flows that erupted on the north side of the rift zone advanced toward Hilo. This flank of the volcano may be more vulnerable to inundation. Lockwood (1990) noted that the vents of historical activity are migrating to the south. The volcano appears to have a self-regulating mechanism that evenly distributes long-term activity across its flanks. The geologic record also supports this notion; the time prior to the historical period (Age Group 1, pre-A.D. 1832 to 1,000 yrs B.P.; orange units) is dominated by activity on the south side of the NERZ. Although most Mauna Loa eruptions begin in the summit area at the 12,000-ft elevation (Lockwood and Lipman, 1987), the central-southeast flank has not been the source of any activity. All flows originated from the summit or the upper reaches of the Northeast Rift Zone (NERZ) or the Southwest Rift Zone (SWRZ). The NERZ was the source of eight flank eruptions since 1843. The NERZ extends from the 13,680-ft-high summit towards Hilo (population ~60,000; second-largest city in State of Hawaii). The northern portion of the map area is built entirely on flows erupted from the NERZ. The SWRZ extends from the summit towards Kalae (South Point) at sea level. The southern portion of the map area is built entirely on flows erupted from the SWRZ. The map area extends from the 10,350 -ft elevation on Mauna Loa’s east flank toward the Hawaii Volcanoes National Park and the town of Volcano (population approx. 2,000) in the northeast. At the south boundary of the map area is the town of Pāhala (population approx. 900). This map includes areas adjacent to and downslope of the NERZ and regions east of and directly downslope of Moku‘āweoweo, Mauna Loa’s summit caldera. The map encompasses 506 km 2 of the southeast flank (fig. 1) of Mauna Loa from 10,350 -ft elevation to sea level. The map of the central-southeast flank of Mauna Loa shows the distribution and relations of volcanic and surficial sedimentary deposits separated into 15 age groups ranging from a period greater than 50,000 yr B.P. to A.D. 1984. It incorporates previously reported work published in generalized small-scale maps (Lockwood and Lipman, 1987; Lockwood, 1995; Wolfe and Morris, 1996). This map is the second in a series of five maps that will cover Mauna Loa volcano. See SIM 2932-A at https://doi.org/10.3133/sim2932A . NOTE: Map sheet 1 contains lines and type with overprint. This feature may be turned on or off in the Adobe Acrobat page display preferences.

Hawaii↗

USGS2004

Hydrology and cycling of nitrogen and phosphorus in Little Bean Marsh: A remnant riparian wetland along the Missouri River in Platte County, Missouri, 1996–97

The lack of concurrent water-quality and hydrologic data on riparian wetlands in the Midwestern United States has resulted in a lack of knowledge about the water-quality functions that these wetlands provide. Therefore, Little Bean Marsh, a remnant riparian wetland along the Missouri River, was investigated in 1996 and 1997 primarily to determine the magnitude and character of selected water-quality benefits that can be produced in such a wetland and to identify critical processes that can be managed in remnant or restored riparian wetlands for amelioration of water quality. Little Bean Marsh averages 69 hectares in size, has a maximum depth of about 1 meter, and the majority of the marsh is covered by macrophytes. In 1997, 41 percent of the water received by Little Bean Marsh was from direct precipitation, 14 percent was from ground-water seepage, 30 percent from watershed runoff, and 15 percent was backflow from Bean Lake. Although, Little Bean Marsh was both a ground-water recharge and discharge area, discharge to the marsh was three times the recharge to ground water. Ground-water levels closely tracked marsh water levels indicating a strong hydraulic connection between ground water and the marsh. Reduced surface runoff and ground-water availability are stabilizing influences on marsh hydrology and probably contribute to the persistence of emergent vegetation. The rapid hydraulic connection between Little Bean Marsh and ground water indicates that the hydrologic regime of most wetlands along the lower Missouri River is largely a function of the altitude of the marsh bottom relative to the altitude of the water table. More water was lost from the marsh through evapotranspiration (59 percent) than all other pathways combined. This is partially because the transpiration process of abundant macrophytes can greatly contribute to the evapotranspiration above that lost from open water surfaces. Surface outflow accounted for 36 percent and ground-water seepage accounted for only 5 percent of the losses. Large residence times allows the marsh to greatly affect water quality before water escapes as ground-water recharge or surface outflow. The shallowness of Little Bean Marsh and ion exclusion during ice formation caused the highest specific conductances of 1,100 to 1,300 microsiemens per centimeter at 25 degrees Celsius to occur during the winter. This concentration of dissolved solutes under ice can make wetlands more vulnerable to toxic contaminants than deeper surface-water bodies. Dissolved oxygen was less than 5 mg/L (milligrams per liter) for 3 to 4 months and near 0 mg/L for about 1 month in summer. Despite depths of less than 1 meter, temperature stratification persisted more than 3 months during the summers of 1996 and 1997, preventing mixing and contributing to periods of anoxia. Shallow depths and extended periods of anoxia in the marsh limit the ability of some organisms to escape high-temperature stress. Turbidity in Little Bean Marsh usually was low for several reasons: sediment loadings from the largely flood-plain drainage were low, emergent vegetation shade out algae and shield the water from wind, and high concentrations of bivalent cations increase flocculation rates of inorganic suspended material. The high concentrations of bivalent cations was largely because of a substantial amount of ground-water seepage into the marsh. Dissolved organic nitrogen was the dominant nitrogen species in Little Bean Marsh. Denitrification and biotic uptake kept more than 62 percent of nitrate (NO3) and 43 percent of ammonium (NH4) concentrations in marsh samples less than a detection limit of 0.005 mg/L. This contrasts with the Missouri River where inorganic NO3 dominates. Consequently, artificial flood-plain drainage that bypasses riparian wetlands likely deliver substantially more biotically available inorganic nitrogen to receiving waters than surface water that has been routed through a remnant wetland. Average total nitrogen concentrations in Little Bean Marsh were substantially less than those at other Missouri River wetlands, roughly one-half the mean concentrations in the Missouri River, but roughly twice the average nitrogen values in reservoirs of the glaciated plains of Missouri. The largest concentrations of nearly all species of nitrogen and phosphorus and the most intense period of hypereutrophy coincided with a phytoplankton bloom and senescence of River Bulrush (Scirpus fluviatilis) and common cattail (Typha latifolia) in September 1997. The rapid leaching of nitrogen that occurs soon after macrophyte senescence combined with a recent destratification of the marsh probably provided nitrogen to the nitrogen-limited open-water areas and triggered a phytoplankton bloom. Despite the rarity of runoff events, surface runoff from the watershed, combined with atmospheric deposition, contributed more than seven times the 530 kg (kilograms) of nitrogen that escaped Little Bean Marsh in surface outflow during 1997. Atmospheric deposition alone was more than 530 kg. Seepage to ground water contained less than 1.5 percent of the nitrogen leaving the marsh in surface outflow. The slow decay rate of Scirpus fluviatilis and reducing conditions in bottom sediments make burial of organic nitrogen a substantial sink of nitrogen. Denitrification experiments indicate that denitrification rates were limited by NO3 in the water column. Consequently, decomposition and nitrification of NH4 and organic nitrogen are the rate limiting steps of nitrogen removal in Little Bean Marsh. The NO3-limited rates of denitrification also indicate that Little Bean Marsh has a large unused capacity for nitrogen removal. These data indicate that the vast extent of riparian marshes along the Missouri and Mississippi Rivers may have had a substantial role in limiting NO3 loads to the Gulf of Mexico before agricultural development of flood plains. Drainage and removal of riparian marshes may be a major cause of the increased NO3 loads to the Gulf of Mexico. Periods of anoxia had much larger effects on phosphorus release than the other variables. The largest concentrations of phosphorus occurred in late summer and corresponded with senescing macrophytes, periods of anoxia, and a large algal bloom in Little Bean Marsh. Low water levels prevented the escape of phosphorus in surface outflow during these periods of highest phosphorus concentrations. Dry weather in late summer is typical and probably makes the correspondence of low water levels, anoxia, and consequent low phosphorus release a common occurrence in marshes along the Missouri River. Little Bean Marsh retained more than 95 percent of the phosphorus it received. The amount of phosphorus in surface inflows to the marsh were more than one order of magnitude greater than that escaping in surface outflows. The long hydraulic residence time of the marsh and large contributions of iron from ground water (that provide many sorption sites for phosphorus) make the marsh an effective sediment and phosphorus trap.

Missouri↗

USGS2025

Assessment of channel morphology, hydraulics, and bedload transport along the Siletz River, western Oregon

Significant Findings Chinook salmon ( Oncorhynchus tshawytscha ) and Pacific lamprey ( Entosphenus tridentatus ) are native, anadromous fish species in the Siletz River Basin, western Oregon, that face many threats to their survival in freshwater and the ocean. The Confederated Tribes of Siletz Indians of Oregon seek to mitigate freshwater threats to Chinook salmon and Pacific lamprey, where possible, with habitat conservation and restoration efforts. This study was conducted to assist the Confederated Tribes of Siletz Indians of Oregon in documenting and understanding the hydrogeomorphic processes shaping present-day habitat conditions and assessing future habitat implications for Chinook salmon and Pacific lamprey along the main-stem Siletz River. As such, this study focused on understanding geomorphic processes and patterns of channel change, including lateral and vertical adjustments in channel position and changes in bed-material sediment (sands, gravels, and cobbles that mantle the channel bed), which collectively determine overall patterns of channel morphology and fluvial habitats. Objective One was to evaluate lateral changes in channel position, vertical changes in bed elevation, and longitudinal patterns in bed-material particle size along the Siletz River using detailed channel maps developed from aerial photographs collected from 1939 to 2016, long-term records of stage and discharge collected by the U.S. Geological Survey (USGS) near the City of Siletz, and sediment particle size data. Objective Two was to assess hydraulic conditions using one- and two-dimensional hydraulic models and transport capacity of bed-material sediment using bedload transport models and sediment particle size data for a range of discharge conditions. Objective Three was to identify potential burrowing habitat for lamprey larvae (PBH) along the Siletz River network and provide insights in local factors influencing PBH along the main-stem Siletz River. The overall findings are synthesized to describe habitat implications for Chinook salmon and Pacific lamprey under present-day and future conditions. Results of Objective One, an evaluation of changes in channel position and bed elevations and longitudinal patterns in bed-material particle size along the Siletz River, include the following From 1939 to 2016, channel planform and the mapped area of gravel bars did not change considerably along the 97.2-kilometers (km) of Siletz River between Elk Creek and Millport Slough, except for in short sections generally less than 1 km long. This inherent lateral channel stability results from the resistant bedrock and terraces that bound most of the channel and limit lateral changes in channel position. Intermittent sections along the study reach where the active channel widened at channel bends displayed noticeable planform changes and increases in mapped bar area. From water year (WY; a 12-month period from October 1 through September 30 and named for the year in which it ends) 1906 to 2021, changes in the stage-discharge relation interpreted as rising and lowering channel bed levels were observed at the USGS streamgage 14305500 (Siletz River at Siletz, OR) in response to floods (such as high flows in February 1996 and January 2002 that exceeded 0.1 and 0.667 annual exceedance probabilities [AEP] events, respectively). However, the rating curve representing the stage-discharge relation did not change in response to high-magnitude floods between 2007 and 2021. Along the approximately 54-km of the Siletz River between Moonshine Park and the Bulls Bag area, surficial particle distributions varied considerably between sampling sites in response to changes in channel width and gradient, sediment inputs from tributaries, and basin geology. Despite this variability, median particle sizes tended to decrease in size in the downstream direction over the 54 km. Results of Objective Two, an evaluation of hydraulic and bedload transport conditions along the Siletz River, include the following The most substantial increases in maximum and mean water velocity and bedload sediment transport capacity occurred at events between the 0.995 and 0.50 AEPs. Events of these magnitudes occur approximately every 1–2 years. Smaller events (0.995 AEP) are generally contained by the banks of the main channel, whereas larger events (0.50 AEP) generally spill over the banks and inundate high-elevation bar and low-elevation floodplain surfaces. Multiple smaller floods within a WY that exceed the 0.995 and 0.50 AEPs (such as in WY 1996) can transport as much or more bed-material sediment than a single, higher magnitude event (such as the maximum event recorded in WY 2000 with an estimated AEP of around 0.002). Bedload transport capacity generally exceeds sediment supply (greater than 2 millimeters [mm]) for most of the study area from Wildcat Creek to the City of Siletz as evidenced by substantial in-channel bedrock. Despite overall conditions of limited sediment supply, transport capacity still varies considerably within and between years with discharge magnitude and spatially in relation to local hydraulics imposed by bedrock, channel morphology, and human infrastructure (such as bridges). Results of Objective Three, an analysis of PBH for lamprey larvae, include the following About 28 percent of the Siletz River network meets the mean annual suspended sediment loads and channel slope criteria for PBH. Along the main-stem Siletz River, in-channel bedrock outcrops and high transport capacity are expected to further constrain PBH. Together, these results suggest that most of the Siletz River between Wildcat Creek and the City of Siletz has had only modest vertical and lateral change between the 1930s and 2010s because of the bedrock in and along the main channel and the river’s relatively high transport capacity relative to bed-material sediment supply. However, localized sections of the Siletz River where the active channel widens, particularly at channel bends, exhibited some change in channel planform and the locations and area of gravel bars. In the future, moderate increases in autumn-winter discharge may not result in substantial changes in coarse gravel bars along the Siletz River but may result in selective transport of finer bed-material sediment (gravel, sands, and silts) that provide spawning habitats for Chinook salmon and Pacific lamprey and burrowing habitats for lamprey larvae. Assuming no substantial changes in bed-material sediment supply, increased bedload transport capacity may cause frequent entrainment of lamprey larvae that are burrowed in coarse sand deposits, suspension and downstream transport of salmon eggs incubating in gravels, and reductions in the areas of spawning gravels for Chinook salmon and Pacific lamprey. Exact implications of current and future discharge conditions for these species along the Siletz River depends on many factors, including sediment supply, local hydraulics, and the timing of flood events relative to fish life stages.

Oregon↗

USGS2020

Spatial and temporal patterns in streamflow, water chemistry, and aquatic macroinvertebrates of selected streams in Fairfax County, Virginia, 2007–18

Urbanization substantially alters the landscape in ways that can impact stream hydrology, water chemistry, and the health of aquatic communities. Stormwater best management practices (BMPs) are the primary tools used to mitigate the effects of urban stressors such as increased runoff, decreased baseflow, and increased nutrient and sediment transport. To date, Fairfax County Virginia’s stormwater management program has made substantial investments into the implementation of both structural and nonstructural BMPs aimed at restoring and protecting watersheds. The U.S. Geological Survey (USGS), in cooperation with Fairfax County, Virginia, established a long-term water-resources monitoring program to evaluate the watershed-scale effects of these investments. Monitoring began at 14 stations in 2007 and was expanded to 20 stations in 2013. This report utilized the first 10 years of data collection to (1) assess water quantity and quality, as well as ecological condition; (2) compute annual nutrient and sediment loads; and (3) evaluate trends in streamflow, water quality, and ecological condition. Efforts are underway to link the biotic and abiotic patterns described herein to watershed management practices as well as factors such as land use change, public works infrastructure, and climate. Hydrologic, chemical, and benthic macroinvertebrate community conditions in the streams monitored were similar to those observed in other studies of urban streams. Multidecadal trends in baseflow indices and runoff ratios at long-term Chesapeake Bay Non-tidal Network streamgages (CB-NTN) indicate a decrease in groundwater recharge and increase in storm runoff as a result of urbanization. Streamflow yields varied spatially with land cover, geology, and soil characteristics, whereas flashiness was positively related to impervious area. Dissolved oxygen typically was lowest in the Coastal Plain and across all Triassic Lowlands streams, and highest in the Piedmont. Dissolved oxygen concentrations generally were above Virginia’s minimum criterion of 4.0 milligrams per liter (mg/L), most violations occurred at Paul Spring Branch in the Coastal Plain during the warmest months of the year owing to increased chemical and biological oxygen demand. Typical pH values of the monitored streams centered on neutrality (pH = 7); however, diurnal fluctuations were most prevalent in the continuous pH data at Flatlick Branch (FLAT; a Triassic Lowlands station), as a result of increased photosynthesis catalyzed by phosphorus-rich geology. Specific conductance (SC) varied spatially owing to geology (highest at Triassic Lowlands stations) and anthropogenic disturbance (watersheds with high impervious land cover). Specific conductance typically was inversely related to streamflow except in winter months following deicing road salt applications, when values increased by several orders of magnitude. A significant increase in SC of about 2 percent per year was observed from the combined trend result of all monitoring stations over the 10-year period. Significant SC increases occurred at nearly all monitoring stations. Increasing trends were observed during winter and nonwinter months, which suggests that salts applied to deice roadways and other impervious surfaces are stored in the environment and released year-round. Suspended-sediment (SS) concentrations in monthly samples did not vary significantly between most stations, but typically were highest in the spring and lowest in the fall as a result of seasonal differences in streamflow and climate. Suspended-sediment yields ranged from 62 to 1,428 tons per square mile (ton/mi 2 ), with a median of 302 ton/mi 2 . Annual loads were greatest during the wettest water years (October 1-September 30; 2008, 2011, and 2014), with the greatest interannual variability occurring at Difficult Run above Fox Lake (DIFF) and South Fork Little Difficult Run (SFLIL). Suspended sediment was primarily composed of silts and clays; however, the proportion of sand in suspended sediment was related positively to streamflow. Cross-correlation analyses suggested the dominant sources of SS were streambank erosion and resuspension of in-channel material at DIFF and FLAT; whereas, upland sources and erosion of upper streambanks were more common at Dead Run (DEAD), Long Branch (LONG), and SFLIL. Median total phosphorus (TP) concentrations ranged from 0.016 to 0.077 mg/L, with a networkwide median of 0.022 mg/L, were highest in the warm season (April-September), and were composed primarily of dissolved phosphorous. Although TP concentrations were relatively low across the network, the highest concentrations were consistently at stations located in the Triassic Lowlands, owing to phosphorous-rich geology, and in the Coastal Plain, owing to the low-phosphorous sorptive capacity of those soils. A significant increase in TP concentration occurred in a few stations, but the combined trend results from all stations demonstrated a significant increase of about 4 percent per year. Networkwide increases were also observed in total dissolved phosphorus, orthophosphate, and total particulate phosphorus. The composition of TP shifted from dissolved to particulate as streamflow increased and for this reason loads primarily were composed of particulate phosphorous. Median annual TP loads were highest at FLAT and DEAD and ranged from 247 to 642 pounds per square mile (lbs/mi 2 ) networkwide. Interannual variability in phosphorous yields was apparent at most stations; the highest loading years were also the wettest years during the study period and coincident with the highest peak annual flows. Total nitrogen (TN) concentrations typically were low throughout the network with exceptions occurring at stations located in watersheds with a high density of septic infrastructure. Elevated TN concentrations also were observed in some watersheds without a high density of septic systems and may be attributable to geologic and soil properties that limit denitrification as well as other unknown anthropogenic inputs. Total nitrogen typically was dominated by nitrate during baseflows; however, the proportion of particulate nitrogen increased during stormflows. Total nitrogen yields were similar across stations, with medians ranging from about 3,600 to 6,300 lbs/mi 2 and were related to annual streamflow volume. Total nitrogen concentrations and flow-normalized concentrations decreased over the 10-year period at 7 stations, with median reductions of about 2.5 percent. Increasing trends were observed at the two stations with the highest median TN concentration (Captain Hickory Run and SFLIL, 3–5 mg/L), both watersheds contain a high density of septic infrastructure. The combined trend results from all stations revealed no trend in TN and a declining trend in nitrate of about 2 percent per year. Overall, benthic community metrics indicated that streams throughout Fairfax County were initially of poor health; however, many metrics show an improving trend (from poor to fair based on the Fairfax County Index of Biological Integrity [IBI]). Significant increasing trends in IBI occurred at the network-scale and at 4 individual stations; additionally, scores improved by at least 1 qualitative category (for example, poor to fair, fair to good) at 11 of the 14 stations between 2009 (the first year all 14 stations were sampled) and 2017. Changes in all metrics suggest that the biodiversity, function, and condition of streams in Fairfax County are improving, but some of these improvements are driven by increased diversity and percent composition of organisms that are tolerant of the urban environment.

Virginia↗

USGS2016

Simulated effects of groundwater withdrawals from aquifers in Ocean County and vicinity, New Jersey

Rapid population growth since the 1930s in Ocean County and vicinity, New Jersey, has placed increasing demands upon the area’s freshwater resources. To examine effects of groundwater withdrawals, a three-dimensional groundwater-flow model was developed to simulate the groundwater-flow systems of five area aquifers: the unconfined Kirkwood-Cohansey aquifer system and Vincentown aquifer, and three confined aquifers— the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, and the Piney Point aquifer. The influence of withdrawals is evaluated by using transient groundwater-flow model simulations that incorporate three withdrawal schemes. These are (1) no-withdrawal conditions; (2) 2000–03 withdrawal conditions, using reported monthly withdrawals at all production wells from January 2000 through December 2003; and (3) maximum-allocation withdrawal conditions using the maximum withdrawal allowed by New Jersey Department of Environmental Protection permits at each well. Particle tracking analysis, using results from model simulations, delineated particle flow paths from production wells to the point of recharge, and estimated particle travel times. Compared with no-withdrawal conditions, 2000–03 withdrawal conditions reduced the amount of groundwater flow out of the Kirkwood-Cohansey aquifer system into streams, increased the net flow of water into other layers, reduced net flow into or out of storage, and reduced flow from the Kirkwood-Cohansey aquifer system to constant head cells. Freshwater discharging to the Barnegat Bay-Little Egg Harbor estuary from streams and groundwater is essential to maintaining the ecology of the bay. Examination of selected stress periods indicates that simulated base flow in streams flowing into the Barnegat Bay-Little Egg Harbor estuary is reduced by as much as 49 cubic feet per second for 2000 to 2003 withdrawal conditions when compared with no-withdrawal conditions. In the three confined aquifers, water levels during periods of low recharge and high withdrawals, and high recharge and low withdrawals, were examined to determine seasonal effects on the confined flow systems. The simulated potentiometric surface of the Rio Grande water-bearing zone and the Atlantic City 800-foot sand during selected stress periods indicates substantial declines from no-withdrawal conditions to 2000–03 conditions as a result of groundwater withdrawals. Cones of depression in Toms River Township, Seaside Heights and Seaside Park Boroughs, and Barnegat Light Borough developed in the potentiometric surface of the Piney Point aquifer in response to withdrawals. Maximum-allocation withdrawals decreased flow out of the Kirkwood-Cohansey aquifer system to constant head cells, increased flow out of the aquifer system to adjacent and lower layers, and reduced groundwater discharge to streams when compared with 2000–03 withdrawal conditions. Increases in withdrawals from the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, and the Piney Point aquifer result in an increase in simulated net groundwater flow into these aquifers. Base-flow reduction from 2000–03 conditions to maximum-allocation conditions of 25 to 29 cubic feet per second in all streams draining to the Barnegat Bay-Little Egg Harbor also is indicated. Potentiometric surfaces of the Rio Grande water-bearing zone, Atlantic City 800-foot sand, and the Piney Point aquifer during two stress periods of simulated maximum-allocation withdrawal conditions indicated the expansion of several cones of depression developed during 2000–03 withdrawals. Simulation of average 2000–03 withdrawal conditions indicated the extent to which the groundwater-flow system is susceptible to potential saltwater intrusion into near-shore wells. Travel time from recharge to discharge location ranged from 11 to approximately 50,700 years in near-shore Kirkwood-Cohansey aquifer system wells. Those in Seaside Heights Borough, in Island Beach State Park (Berkeley Township), and in Ship Bottom Borough have particle travel times from 140 to 12,000 years and flow paths that originated under Barnegat Bay or the Atlantic Ocean from the simulation of average maximum-allocation withdrawal conditions. Travel time along flow paths to wells screened in the Rio Grande water-bearing zone and the Atlantic City 800-foot sand from recharge to discharge point ranged from nearly 530 years to greater than 3.73 million years from the simulation of average 2000–03 withdrawal conditions. Particle tracking indicated that most wells screened in these aquifers derived a large part of their recharge from the Oswego River Basin, with a small portion of flow originating either beneath Barnegat Bay or to the east beneath the Atlantic Ocean. Travel time along flow paths that start beneath either Barnegat Bay or the Atlantic Ocean ranged from 2,300 to approximately 134,000 years from the simulation of average maximum-allocation withdrawal conditions."

New Jersey↗

USGS2012

Computing daily mean streamflow at ungaged locations in Iowa by using the Flow Anywhere and Flow Duration Curve Transfer statistical methods

The U.S. Geological Survey (USGS) maintains approximately 148 real-time streamgages in Iowa for which daily mean streamflow information is available, but daily mean streamflow data commonly are needed at locations where no streamgages are present. Therefore, the USGS conducted a study as part of a larger project in cooperation with the Iowa Department of Natural Resources to develop methods to estimate daily mean streamflow at locations in ungaged watersheds in Iowa by using two regression-based statistical methods. The regression equations for the statistical methods were developed from historical daily mean streamflow and basin characteristics from streamgages within the study area, which includes the entire State of Iowa and adjacent areas within a 50-mile buffer of Iowa in neighboring states. Results of this study can be used with other techniques to determine the best method for application in Iowa and can be used to produce a Web-based geographic information system tool to compute streamflow estimates automatically. The Flow Anywhere statistical method is a variation of the drainage-area-ratio method, which transfers same-day streamflow information from a reference streamgage to another location by using the daily mean streamflow at the reference streamgage and the drainage-area ratio of the two locations. The Flow Anywhere method modifies the drainage-area-ratio method in order to regionalize the equations for Iowa and determine the best reference streamgage from which to transfer same-day streamflow information to an ungaged location. Data used for the Flow Anywhere method were retrieved for 123 continuous-record streamgages located in Iowa and within a 50-mile buffer of Iowa. The final regression equations were computed by using either left-censored regression techniques with a low limit threshold set at 0.1 cubic feet per second (ft3/s) and the daily mean streamflow for the 15th day of every other month, or by using an ordinary-least-squares multiple linear regression method and the daily mean streamflow for the 15th day of every other month. The Flow Duration Curve Transfer method was used to estimate unregulated daily mean streamflow from the physical and climatic characteristics of gaged basins. For the Flow Duration Curve Transfer method, daily mean streamflow quantiles at the ungaged site were estimated with the parameter-based regression model, which results in a continuous daily flow-duration curve (the relation between exceedance probability and streamflow for each day of observed streamflow) at the ungaged site. By the use of a reference streamgage, the Flow Duration Curve Transfer is converted to a time series. Data used in the Flow Duration Curve Transfer method were retrieved for 113 continuous-record streamgages in Iowa and within a 50-mile buffer of Iowa. The final statewide regression equations for Iowa were computed by using a weighted-least-squares multiple linear regression method and were computed for the 0.01-, 0.05-, 0.10-, 0.15-, 0.20-, 0.30-, 0.40-, 0.50-, 0.60-, 0.70-, 0.80-, 0.85-, 0.90-, and 0.95-exceedance probability statistics determined from the daily mean streamflow with a reporting limit set at 0.1 ft 3 /s. The final statewide regression equation for Iowa computed by using left-censored regression techniques was computed for the 0.99-exceedance probability statistic determined from the daily mean streamflow with a low limit threshold and a reporting limit set at 0.1 ft 3 /s. For the Flow Anywhere method, results of the validation study conducted by using six streamgages show that differences between the root-mean-square error and the mean absolute error ranged from 1,016 to 138 ft 3 /s, with the larger value signifying a greater occurrence of outliers between observed and estimated streamflows. Root-mean-square-error values ranged from 1,690 to 237 ft 3 /s. Values of the percent root-mean-square error ranged from 115 percent to 26.2 percent. The logarithm (base 10) streamflow percent root-mean-square error ranged from 13.0 to 5.3 percent. Root-mean-square-error observations standard-deviation-ratio values ranged from 0.80 to 0.40. Percent-bias values ranged from 25.4 to 4.0 percent. Untransformed streamflow Nash-Sutcliffe efficiency values ranged from 0.84 to 0.35. The logarithm (base 10) streamflow Nash-Sutcliffe efficiency values ranged from 0.86 to 0.56. For the streamgage with the best agreement between observed and estimated streamflow, higher streamflows appear to be underestimated. For the streamgage with the worst agreement between observed and estimated streamflow, low flows appear to be overestimated whereas higher flows seem to be underestimated. Estimated cumulative streamflows for the period October 1, 2004, to September 30, 2009, are underestimated by -25.8 and -7.4 percent for the closest and poorest comparisons, respectively. For the Flow Duration Curve Transfer method, results of the validation study conducted by using the same six streamgages show that differences between the root-mean-square error and the mean absolute error ranged from 437 to 93.9 ft 3 /s, with the larger value signifying a greater occurrence of outliers between observed and estimated streamflows. Root-mean-square-error values ranged from 906 to 169 ft 3 /s. Values of the percent root-mean-square-error ranged from 67.0 to 25.6 percent. The logarithm (base 10) streamflow percent root-mean-square error ranged from 12.5 to 4.4 percent. Root-mean-square-error observations standard-deviation-ratio values ranged from 0.79 to 0.40. Percent-bias values ranged from 22.7 to 0.94 percent. Untransformed streamflow Nash-Sutcliffe efficiency values ranged from 0.84 to 0.38. The logarithm (base 10) streamflow Nash-Sutcliffe efficiency values ranged from 0.89 to 0.48. For the streamgage with the closest agreement between observed and estimated streamflow, there is relatively good agreement between observed and estimated streamflows. For the streamgage with the poorest agreement between observed and estimated streamflow, streamflows appear to be substantially underestimated for much of the time period. Estimated cumulative streamflow for the period October 1, 2004, to September 30, 2009, are underestimated by -9.3 and -22.7 percent for the closest and poorest comparisons, respectively.

Illinois;Iowa;Minnesota;Missouri;Nebraska;Wisconsi↗

USGS2021

Assessment of runoff volume reduction associated with soil amendments added to portions of highway median-strip catchments in Ohio, 2018–20

The U.S. Geological Survey installed 10 rain gages and 12 calibrated H-flumes to measure rainfall and runoff volumes at 10 locations in Ohio Department of Transportation highway median-strip catchments. Data were collected to facilitate comparisons of rainfall and runoff volumes at study sites before and after stormwater best management practices (BMPs) were installed and between sites with different BMPs. The BMP treatments comprised removing the top layer of the existing soil, rototilling the remaining soil to a 6-inch depth, mixing the soils with one of two soil amendments (compost with sand or shale) at one of two thicknesses (4 inches or 6 inches), topping with a compost blanket, seeding, and installing erosion control matting. The overall treatment used at a given study site is referred to as “BMP.” At two locations where soil amendments were installed, a second “control” site was installed to measure runoff from an adjacent catchment in the same median strip where no soil amendment was installed. This no-treatment option (no soil amendment) was considered its own class of BMP. Rainfall and runoff data were collected during periods when air temperatures were above freezing (including all months except January, February, and parts of December and March) from 2018 to 2020. The data collection period for each study site was divided into “pre-BMP” and “post-BMP” periods. Equipment to measure rainfall and runoff was installed and data were collected from April to December 2018 before installation of soil amendments (the pre-BMP period). The post-BMP period started between April and May of 2019 at the first measured rainfall after soil amendments were installed. Rainfall and runoff monitoring continued through September 2020. For control sites, the post-BMP periods were assigned to start with the first measured rainfall in the 2019 data collection season. A rainfall-runoff “event” was defined as beginning at the time of the first measured rainfall and ending when rainfall and runoff (if any) ceased and remained ceased for at least 3 hours. A value referred to as “event runoff percentage,” defined as the total volume of runoff during an event expressed as a percentage of the total volume of rainfall falling over the catchment, was computed for each event. The distribution of rainfall totals associated with events was similar between the pre-BMP and post-BMP periods; however, there were appreciable between-site differences in the distribution of event runoff percentages during the pre-BMP and post-BMP periods. Empirical distribution function (EDF) tests were performed with and without data from events that resulted in no runoff to determine whether the distribution of event runoff percentages changed from the pre-BMP period to the post-BMP period. The null hypothesis that the EDFs of event runoff percentages were equal in the pre-BMP and post-BMP periods was rejected ( α =0.05) in at least one of the two tests for four sites (one site with a shale amendment and three sites with sand amendments). Mean event runoff percentages at each of those four sites decreased from the pre-BMP period to the post-BMP period. The null hypothesis that the EDFs of event runoff percentages were equal was not rejected for the other six sites’ draining catchments with soil amendments or the two control sites. EDF tests performed on event rainfall totals indicated no statistically significant changes between the pre-BMP and post-BMP period distributions for any of the sites. Double-mass analyses of cumulative runoff were performed for two pairs of closely spaced sites (each pair located in a common median strip): one site in each pair drained a catchment where soil amendments were installed, and the other (a control) drained a catchment without soil amendments. Those double-mass analyses indicated a small reduction in runoff from the pre-BMP to post-BMP period at the site whose catchment received the sand and compost amendment, but no perceptible reduction in runoff at the site whose catchment received the shale and compost amendment. Regression analyses indicated that (a) three rainfall factors (event rainfall totals, total rainfall for the previous 7 days, and a cross product of the factors) and the intercept term were the four most important factors explaining event runoff percentages, (b) the effect of amendment type on event runoff percentage was small in comparison to the rainfall and intercept terms, (c) event runoff percentages tended to be lower for sites with shale amendments than sites with sand amendments; however, event runoff percentages tended to be lower for control sites than for sites with shale or sand amendments, and (d) event runoff percentages increased with increasing amendment thickness. The counterintuitive results that event runoff percentages increased with increasing amendment thickness and that control sites tended to have lower event runoff percentages than sites draining soil-amended catchments likely reflects unmeasured factors that existed at the sites before BMPs were installed rather than the effect of the BMP treatments. Although not definitive, some support for the conclusion that the sand amendment was generally more effective at reducing runoff than the shale amendment was provided by results from the EDF tests, double-mass analyses, and runoff statistics.

Ohio↗

USGS2014

California State Waters Map Series: Offshore of Coal Oil Point, California

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within the 3-nautical-mile limit of California’s State Waters. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data, acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. The Offshore of Coal Oil Point map area lies within the central Santa Barbara Channel region of the Southern California Bight. This geologically complex region forms a major biogeographic transition zone, separating the cold-temperate Oregonian province north of Point Conception from the warm-temperate California province to the south. The map area is in the southern part of the Western Transverse Ranges geologic province, which is north of the California Continental Borderland. Significant clockwise rotation—at least 90°—since the early Miocene has been proposed for the Western Transverse Ranges province, and geodetic studies indicate that the region is presently undergoing north-south shortening. Uplift rates (as much as 2.0 mm/yr) that are based on studies of onland marine terraces provide further evidence of significant shortening. The cities of Goleta and Isla Vista, the main population centers in the map area, are in the western part of a contiguous urban area that extends eastward through Santa Barbara to Carpinteria. This urban area is on the south flank of the east-west-trending Santa Ynez Mountains, on coalescing alluvial fans and uplifted marine terraces underlain by folded and faulted Miocene bedrock. In the map area, the relatively low-relief, elevated coastal bajada narrows from about 2.5 km wide in the east to less than 500 m wide in the west. Several beaches line the actively utilized coastal zone, including Isla Vista County Park beach, Coal Oil Point Reserve, and Goleta Beach County Park. The beaches are subject to erosion each winter during storm-wave attack, and then they undergo gradual recovery or accretion during the more gentle wave climate of the late spring, summer, and fall months. The Offshore of Coal Oil Point map area lies in the central part of the Santa Barbara littoral cell, which is characterized by littoral drift to the east-southeast. Longshore drift rates have been reported to range from about 160,000 to 800,000 tons/yr, averaging 400,000 tons/yr. Sediment supply to the western and central parts of the littoral cell, including the map area, is largely from relatively small transverse coastal watersheds. Within the map area, these coastal watersheds include (from east to west) Las Llagas Canyon, Gato Canyon, Las Varas Canyon, Dos Pueblos Canyon, Eagle Canyon, Tecolote Canyon, Winchester Canyon, Ellwood Canyon, Glen Annie Canyon, and San Jose Creek. The Santa Ynez and Santa Maria Rivers, the mouths of which are about 100 to 140 km northwest of the map area, are not significant sediment sources because Point Conception and Point Arguello provide obstacles to downcoast sediment transport and also because much of their sediment load is trapped in dams. The Ventura and Santa Clara Rivers, the mouths of which are about 45 to 55 km southeast of the map area, are much larger sediment sources. Still farther east, eastward-moving sediment in the littoral cell is trapped by Hueneme and Mugu Canyons and then transported to the deep-water Santa Monica Basin. The offshore part of the map area consists of a relatively flat and shallow continental shelf, which dips gently seaward (about 0.8° to 1.0°) so that water depths at the shelf break, roughly coincident with the California’s State Waters limit, are about 90 m. This part of the Santa Barbara Channel is relatively well protected from large Pacific swells from the north and northwest by Point Conception and from the south and southwest by offshore islands and banks. The shelf is underlain by variable amounts of upper Quaternary marine and fluvial sediments deposited as sea level fluctuated in the late Pleistocene. The large (130 km2) Goleta landslide complex lies along the shelf break in the southern part of the map area. This compound slump complex may have been initiated more than 200,000 years ago, but it also includes three recent failures that may have been generated between 8,000 to 10,000 years ago. A local, 5- to 10-m-high tsunami may have been generated from these failure events. The map area has had a long history of hydrocarbon development, which began in 1928 with discovery of the Ellwood oil field. Subsequent discoveries in the offshore include South Ellwood offshore oil field, Coal Oil Point oil field, and Naples oil and gas field. Development of South Ellwood offshore field began in 1966 from platform “Holly,” the last platform to be installed in California’s State Waters. The area also is known for “the world’s most spectacular marine hydrocarbon seeps,” and large tar seeps are exposed on beaches east of the mouth of Goleta Slough. Offshore seeps adjacent to South Ellwood oil field release about 40 tons per day of methane and about 19 tons per day of ethane, propane, butane, and higher hydrocarbons. Seafloor habitats in the broad Santa Barbara Channel region consist of significant amounts of soft sediment and isolated areas of rocky habitat that support kelp-forest communities nearshore and rocky-reef communities in deep water. The potential marine benthic habitat types mapped in the Offshore of Coal Oil Point map area are directly related to its Quaternary geologic history, geomorphology, and active sedimentary processes. These potential habitats, which lie primarily within the Shelf (continental shelf) but also partly within the Flank (basin flank or continental slope) megahabitats, range from soft, unconsolidated sediment to hard sedimentary bedrock. This heterogeneous seafloor provides promising habitat for rockfish, groundfish, crabs, shrimp, and other marine benthic organisms.

California↗

USGS2008

Occurrence and distribution of algal biomass and Its relation to nutrients and selected basin characteristics in Indiana streams, 2001-2005

Algal biomass and nutrient data were gathered at 322 randomly selected sites on 261 streams in the West Fork White River, Whitewater River, East Fork White River, Upper Wabash River, Kankakee River, Lower Wabash River, Tributaries to the Great Lakes, and Tributaries to the Ohio River Basins in Indiana from May through October for years 2001 through 2005. Basin characteristics (land use and drainage area), substrate, turbidity, and nutrient concentrations were determined for the basin and sampling sites. The relations of the seasonal algal biomass parameters periphyton chlorophyll a (CHLa), ash-free dry mass (AFDM), seston CHLa, and particulate organic carbon (POC) to concentrations of the seasonal nutrients nitrate, total Kjeldahl nitrogen (TKN), total nitrogen (TN), and total phosphorus (TP) were determined using Spearman’s rho. The effects of streamflow were determined using data collected at U.S. Geological Survey (USGS) streamflow-gaging stations spatially located throughout the study basins. Throughout the 5-year study, the magnitude and frequency of stream discharge varied monthly and annually and greatly influenced algal biomass concentrations through scour and algal drift. Algal biomass median concentrations in Indiana streams consisted of periphyton CHLa, 41.2 milligrams per square meter (mg/m 2 ); AFDM, 52.1 grams per square meter (g/m 2 ); seston CHLa, 2.44 micrograms per liter (µg/L); and POC, 0.75 milligrams per liter (mg/L). Approximately 32 percent of the periphyton CHLa and 6 percent of the seston CHLa samples would be considered eutrophic (nutrient enriched). To ascertain seasonal variability, samples were collected in the spring (May), summer (June through August), and fall (September through October). The highest median concentration of periphyton CHLa was in the spring, 63.2 mg/m 2 , while the highest median concentrations of AFDM, seston CHLa, and POC were in the summer 55.4 g/m 2 , 2.96 µg/L, and 0.81 mg/L respectively. There were no significant differences among seasons for periphyton CHLa and AFDM; there were significant differences among seasons for seston CHLa and POC. There were no significant relations with nutrients and periphyton or seston CHLa parameters. The only significant positive relations were observed between summer POC and summer TP as well as summer POC and summer TKN. Positive relations also related spring POC and spring TP. These significant relations with TP are most likely related to phosphorus associated within seston algal cells and attached to sediment. Drainage area and land use were analyzed to understand the effect of site location on algal growth. Study basins varied in size (headwater streams, 0–51 km 2 ; wadable streams, 52–2,590 km 2 ; and boatable streams, 2,591–38,900 km 2 ) and were dominated by agricultural land use. Basin characteristics (land use, drainage area) as well as substrate type, turbidity, and nutrients, affected the concentration of algal biomass parameters. Of the eight basins in which samples were collected during the 5-year study, the Whitewater River Basin (2002) had the highest median concentration of periphyton CHLa (63.1 mg/m 2 ), the Tributaries to the Great Lakes (2005) exhibited the highest median concentration for AFDM (160 g/m 2 ), the East Fork White River Basin (2002) had nearly twice the median concentration (4.01 µg/L) of seston CHLa as the other basins, and the West Fork White River Basin (2001) exhibited the highest median concentration of POC (1.10 mg/L). Of the eight major basins sampled, 15–45 percent of the periphyton CHLa and up to 20 percent of the seston CHLa samples were eutrophic. Samples collected at headwater and wadable streams were the most eutrophic for periphyton CHLa (31–36 percent) and 28 percent of samples collected at boatable streams were eutrophic for seston CHLa. As basin size increased, seston CHLa and POC concentrations increased while periphyton CHLa and AFDM concentrations decreased. The median turbidity values ranged from 6.95 NTU for headwater streams to 8.27 NTU for wadable streams, and 17.0 NTU for boatable streams. In addition, the types and availability of periphytic substrates (epilithic, epipsammic, or epidendric) were an important factor when comparing periphyton CHLa and AFDM concentrations in the study due to the periphytic substrates individual susceptibility to bed movement and scouring. Periphyton CHLa median concentrations ranged from 53.8 mg/m 2 for epilithic substrates, to 41.8 mg/m 2 for epipsammic substrates, and 17.2 mg/m 2 for epidendric substrates. Higher AFDM concentrations were collected from epipsammic substrates during years of low stream discharge velocity, which enhanced the settling of organic matter on epipsammic substrates. AFDM median concentrations ranged from 141 g/m 2 for epipsammic to 28.8 g/m 2 for epilithic and 22.9 g/m 2 for epidendric substrates. The seasonal values for nutrients (nitrate, TKN, TN, and TP) and algal biomass (periphyton CHLa, AFDM, seston CHLa, and POC) were compared to published U. S. Environmental Protection Agency (USEPA) values for their respective ecoregions. Algal biomass values either were greater than the 25th percentile published USEPA values or extended the range of data in Aggregate Nutrient Ecoregions VI, VII, IX and USEPA Level III Ecoregions 54, 55, 56, 71, and 72. If the values for the 25th percentile proposed by the USEPA were adopted as nutrient water-quality criteria, then about 71 percent of the nutrient samples and 57 percent of the CHLa samples within the eight study basins would be considered nutrient enriched.

Scientific Investigations Report↗

USGS2018

Groundwater contributions to excessive algal growth in the East Fork Carson River, Carson Valley, west-central Nevada, 2010 and 2012

Excessive algal growth and low dissolved oxygen concentrations were observed during low streamflow conditions during summer months along a 5,800-foot reach of the East Fork Carson River in Carson Valley, west-central Nevada. Algal growth from nutrient enrichment of a stream reduces aquatic diversity, threatens fish ecology and stream health, and can be a recreational nuisance. In response to concerns that groundwater discharging to the 5,800-foot reach of the East Fork Carson River may be a source of nutrients to the stream, the U.S. Geological Survey, in cooperation with the Carson Water Subconservancy District and the Nevada Division of Environmental Protection, conducted studies during the summers of 2010 and 2012 to gain an improved understanding of the contributions of nutrients to the stream from groundwater, characterize algal conditions and algal effects on water quality, assess potential sources of nitrate in groundwater discharging to the stream, and evaluate nitrate reduction in groundwater from denitrification. A reconnaissance study in the summer of 2010 along the 5,800-foot study reach located a subreach with clear evidence of nutrient-rich groundwater discharging to the stream. At the subreach, nitrate plus nitrite (referred to hereafter as nitrate) concentrations in groundwater discharging to the stream were high (average 2.75 milligrams per liter as nitrogen) along the right bank. The stream at this location had the highest stream nitrate concentrations (average 0.056 milligrams per liter as nitrogen) compared to other locations upstream and downstream of the subreach. As a result, the 2012 study focused on a 405-foot subreach of the East Fork Carson River centered where results from the 2010 study found the highest stream and groundwater concentrations of nitrate, as well as the greatest observed contributions of groundwater discharge to the stream. Groundwater nutrient concentrations were much higher than stream nutrient concentrations during the summer of 2012 during low streamflow conditions at the 405-foot subreach of the East Fork Carson River. Average groundwater nitrate and orthophosphate concentrations along the right bank of the 405‑foot subreach were 9 and 12 times higher, respectively, than in the stream at this subreach. Groundwater discharge rates to the study reach based on different methods varied from 0.09 to 1.2 cubic feet per second per mile. Estimated groundwater discharge rates to the right bank of the study subreach were used to calculate groundwater nutrient load estimates to the subreach right bank, which were found to be low when compared to stream nutrient loads. Elevated algal biomass levels above nuisance thresholds were observed during the summers of 2010 and 2012. The study reach was characterized as mesotrophic-eutrophic during the 2010 study and eutrophic during the 2012 study. The presence of algae caused daily dissolved oxygen and pH fluctuations in the stream, resulting in exceedances of the State of Nevada water-quality standards owing to low dissolved oxygen concentrations and high pH concentrations, although the standards might not have been applicable during 2012 because of extremely low streamflow. The addition of nutrients to the stream from the constant supply in groundwater discharge sustains the growth of algae during low streamflow conditions. In the summer when streamflow is low or very low, nutrient-rich groundwater discharge enters the stream through the sediment-water interface at the streambed. Because the attached algae is thick and stream velocity is low, the nutrient-rich water pools at the sediment-water interface. Higher nutrient concentrations at the streambed create a favorable microenvironment for algae attached to the substrate to consume available nutrients from the groundwater before the groundwater mixes with overlying stream water. The source of nitrate in groundwater in this subreach is anthropogenic because nitrate concentrations are greater than background groundwater nitrate concentrations in Douglas County, high groundwater nitrate concentrations are only found at the right bank of the stream near a housing development, and organic wastewater compounds indicative of human-derived sources were also detected in groundwater wells on the right bank of the stream. Nitrogen and oxygen isotope concentrations of nitrate in shallow groundwater were used to determine the specific source of the nitrate, but the isotopic values indicated denitrification was occurring. Further investigation is needed to determine the specific anthropogenic source of the nitrate in the groundwater because the denitrification present in all samples obscures the original source of nitrogen.

Nevada↗

USGS2008

Evaluation of four structural best management practices for highway runoff in Beaufort and Colleton Counties, South Carolina, 2005–2006

From 2005 to 2006, the U.S. Geological Survey worked cooperatively with the South Carolina Department of Transportation in Beaufort and Colleton Counties, South Carolina, to assess the performance of four different structural devices that served as best management practices (BMPs). These structural devices were installed to mitigate the effects of stormwater runoff on waterways near State roads. The South Carolina Department of Transportation is required to address the quality of stormwater runoff from State-maintained roadways as part of the National Pollutant Discharge Elimination System stormwater program mandated in the Clean Water Act. The performance assessment of the four structural best management practices was based on stormflow measurements and chemical analyses of stormwater-quality samples collected during a 20-month period from March 2005 through October 2006, which represented a range of seasons and rainfall intensities. A total of 49 sample sets that included stormwater from the inlet and outlet of each of the four structural devices were collected as flow-weighted composites to provide event-mean concentrations of suspended sediment, nutrients, and trace metals. In addition, each set included grab samples that were collected to provide the first flush concentrations of oil and grease and fecal-indicator bacteria. A tiered statistical approach was used in the data analysis. Performances of the four structural BMPs were assessed individually based on how well the BMPs were able to reduce the selected constituents. Descriptive statistics and nonparametric Wilcoxon signed rank tests were applied to event-mean concentrations and loads in the water entering the inlet and in the water leaving the outlet of each BMP for each constituent to identify if significant reductions occurred. If significant reductions occurred, the BMP was considered efficient at reducing the constituent. To quantify efficiency, a simplistic approach was applied to compute mean and geometric mean efficiency ratios for the significantly reduced constituents in each BMP. Each BMP performance was ranked based on its computed efficiency ratios, however, the computed efficiency ratios were not sufficient to determine if statistical differences occurred among the performances of the four BMPs. Consequently, a more complex approach was used to apply statistical comparison tests to reduction percentages computed for individual storms (a modified removal efficiency of individual storm-load approach) to determine if differences in event-mean concentrations, loads, and reduction percentages for significantly reduced constituents occurred among the four structural BMPs. Overall, the four BMPs were efficient in reducing suspended-sediment event-mean concentrations and loads in the stormwater entering the inlets of the BMPs to significantly lower event-mean concentrations before discharging the stormwater from the outlets. The cumulative suspended-sediment event-mean load in stormwater entering the BMPs from the storms sampled during the data-collection period was 1,026 kilograms (1.13 tons). The BMPs removed a cumulative suspended-sediment load of 558 kilograms (0.62 ton). The BMPs tended to preferentially trap the sand-size fraction of the sediment, thereby releasing a greater percentage of fine-grained (silt and clay) sediment in the water discharging from the outlet. The preferential trapping of fine-grained sediment by the BMPs could explain, at least in part, why the BMPs were not successful at significantly reducing these constituents. In general, the four BMPs were not successful at significantly reducing fecal bacteria, nutrients, and total organic carbon (including associated properties of biochemical oxygen demand and chemical oxygen demand). Three of the four BMPs significantly lowered oil and grease concentrations before the stormwater discharged from the outlet. Additionally, only one BMP was effective at reducing all total and particulate trace-metal event-mean concentrations and particulate trace-metal event-mean loads in stormwater entering the inlet. With respect to trace-metal event-mean concentrations, however, minimal or no improvement in outlet water quality was observed for the four BMPs, and the majority of the outlet concentrations were above the established acute and chronic aquatic-life criteria by the South Carolina Department of Health and Environmental Control. No statistical differences among the removal-efficiency of the four BMPs were determined for suspended-sediment event-mean concentrations, total suspended solids event-mean concentrations, or oil and grease concentrations. These statistical findings indicated that differences among the mean efficiency ratios were not significant among the BMPs for these properties. Additionally, one BMP generally had statistically greater removal efficiency for total and particulate cadmium, copper, lead, and zinc than one or more of the other three BMPs. Statistical correlation tests were unable to identify a single major factor that would explain the high variability in inlet and outlet water concentrations and in removal efficiencies estimated by reduction percentage. Highly variable inlet and outlet concentrations for each BMP that produced highly variable reduction percentages were probably the result of multiple interacting factors, particularly rainfall intensity, the amount of rainfall between sampling events, traffic density, and the period of time since the last maintenance (clean out) of the BMP.

South Carolina↗

USGS2018

Simulation of groundwater flow and analysis of projected water use for the Rush Springs aquifer, western Oklahoma

The U.S. Geological Survey, in cooperation with the Bureau of Reclamation and the Oklahoma Water Resources Board, (1) quantified the groundwater resources of the Rush Springs aquifer in western Oklahoma by developing a numerical groundwater-flow model, (2) evaluated the effects of estimated equal-proportionate-share (EPS) pumping rates on aquifer storage and streamflow for time periods of 20, 40, and 50 years into the future, (3) assessed the uncertainty in the EPS scenario results, and (4) evaluated the effects of (a) projected groundwater-use rates extended 50 years into the future and (b) sustained hypothetical drought conditions over a 10-year period on stream base flow and groundwater in storage. The Rush Springs aquifer is an important source of water for municipal and irrigation use by many communities and agricultural users in the study area. The study area is composed of about 4,970 square miles (3,181,003 acres) of Rush Springs aquifer bedrock deposits located in 14 counties. The study area also includes the alluvium and terrace deposits of the Canadian and Washita Rivers, as well as alluvium along the Little Washita River, Deer Creek, and a number of smaller tributaries of the Washita River that overlie the bedrock. A numerical groundwater-flow model of the Rush Springs aquifer was constructed by using MODFLOW with the Newton solver. Groundwater flow was simulated for January 1979–December 2015 by using monthly stress periods, and an initial steady-state stress period was configured to represent mean annual inflows and outflows. The model was calibrated to groundwater-level observations at selected wells, monthly base flow at nine streamgages, stream seepage as estimated for the conceptual water budget, and Fort Cobb Reservoir stage. The EPS scenarios for the Rush Springs aquifer were run for periods of 20, 40, and 50 years. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were 0.82, 0.49, and 0.43 acre-foot per acre per year, respectively. Given the 2,954,545-acre aquifer area used for the EPS scenarios, the 20-year rate corresponds to an annual yield of about 2,422,727 acre-feet per year. Groundwater storage at the end of the 20-year EPS scenario was about 13,321,000 acre-feet, or about 31,516,437 acre-feet (70 percent) less than the starting EPS scenario storage. This decrease in storage was equivalent to a mean groundwater-level decline of about 152 feet. Water availability under the EPS pumping rate was primarily from the western area of the model. Saturation was sustained though the entire EPS scenario where the aquifer was sufficiently thick or a shallow hydraulic gradient was present. Fort Cobb Reservoir stage was below the dead-pool stage after about 5 years of 20-year EPS pumping. An uncertainty analysis was conducted to assess the uncertainty in the EPS scenario results. An ensemble of 400 random sets of possible parameter values was performed for the uncertainty analysis by using a multivariate normal distribution centered on the calibrated parameter values. The parameter bounds for the uncertainty analysis were determined by using the posterior covariance matrix, which allows for the incorporation of knowledge gained during the calibration process as well as observation uncertainty and the correlation between estimated parameters. The uncertainty results indicate a 95-percent confidence interval for the 20-year EPS pumping rate between 0.73 and 0.95 acre-foot per acre per year. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Rush Springs aquifer. The effects of well withdrawals were evaluated by comparing changes in groundwater storage between four 50-year scenarios using (1) no groundwater use, (2) mean groundwater use for the study period (1979–2015), (3) increasing groundwater use, and (4) groundwater use at the 2015 rate. The increasing-use scenario assumed a 38-percent increase in pumping over 50 years on the basis of 2010–60 demand projections for western Oklahoma. Simulated groundwater storage changes ranged between an increase of 6.3 percent for the scenario with no groundwater use, and 0.9 percent for the scenario with 2015 groundwater-use rates. For the Fort Cobb Reservoir surface watershed, simulated groundwater storage changes ranged between an increase of 23.6 percent for the scenario with no groundwater use and a decrease of 4.0 percent for the increasing groundwater-use scenario. Groundwater-level changes were generally greater in areas with a large concentration of groundwater wells and groundwater use such as the Fort Cobb Reservoir surface watershed. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on the Rush Springs aquifer groundwater storage and Fort Cobb Reservoir stage and storage. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model. To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–1992), and upstream inflows to the Canadian and Washita Rivers and associated tributaries were reduced by 37 percent. Groundwater storage at the end of the hypothetical drought period in December 1992 was about 42,983,000 acre-feet, or about 3,525,000 acre-feet (7.6 percent) less than the groundwater storage of the calibrated numerical model. This change in groundwater storage is equivalent to a mean groundwater-level decline of 15.8 feet. Simulated mean base-flow declines at the Canadian and Washita River streamgages were between 39 and 59 percent during the drought period. The minimum stage in Fort Cobb Reservoir at the end of the hypothetical drought period was 1,311 feet, indicating a storage capacity of only 10 percent of active conservation pool storage. The Fort Cobb Reservoir storage declines mostly resulted from reduced base flows in Cobb, Lake, and Willow Creeks upstream from the reservoir.

Oklahoma↗

USGS2016

Sediment loads and transport at constructed chutes along the Missouri River - Upper Hamburg Chute near Nebraska City, Nebraska, and Kansas Chute near Peru, Nebraska

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, monitored suspended sediment within constructed Missouri River chutes during March through October 2012. Chutes were constructed at selected river bends by the U.S. Army Corps of Engineers to help mitigate aquatic habitat lost through the creation and maintenance of the navigation channel on the Missouri River. The restoration and development of chutes is one method for creating shallow-water habitat within the Missouri River to meet requirements established by the amended 2000 Biological Opinion. Understanding geomorphic channel-evolution processes and sediment transport is important for the design of chutes, monitoring and maintenance of existing chutes, and characterizing the habitat that the chutes provide. This report describes the methods used to monitor suspended sediment at two Missouri River chutes and presents the results of the data analysis to help understand the suspended-sediment characteristics of each chute and the effect the chutes have on the Missouri River. Upper Hamburg chute, near Nebraska City, Nebraska, and Kansas chute, near Peru, Nebraska, were selected for monitoring. At each study site, monthly discrete samples were collected from April through October in the Missouri River main-channel transects upstream from the chute inlet, downstream from the chute outlet, at the outlet (downstream transect) of both chutes, and at the inlet (upstream transect) of Kansas chute. In addition, grab samples from all chute sampling locations were collected using autosamplers. Suspended-sediment concentration (SSC) and grain-size metrics were determined for all samples (discrete and grab). Continuous water-quality monitors recorded turbidity and water temperature at 15-minute intervals at the three chute sampling locations. Two acoustic Doppler velocimeters, one within each chute, measured water depth and current velocities continuously. The depth and velocity data were used to estimate streamflow within each chute. The sampling design was developed to understand the suspended-sediment differences within each chute and between the chute and the Missouri River main channel during discrete sampling. The sampling design also allowed for site-specific surrogate relations between SSC and turbidity to be developed, which could be used to compute real-time estimates of SSC and sediment loads within the chutes. Real-time estimates of SSC and sediment loads enable a better understanding of sediment transport within the chutes during times when physical samples are not collected, including periods of high flow. High flows during the summer of 2011 resulted in substantial alterations to both studied chutes; therefore, the U.S. Army Corps of Engineers repaired and modified both chutes during 2012. These unforeseen repairs and modifications within the chutes added uncertainty to the analysis because concentrations were altered by construction equipment and flow alteration. Daily suspended-sediment and suspended-silt loads were estimated based on surrogate relations with turbidity. A linear regression was used to estimate equal-width increment (EWI)-equivalent SSC from autosampler SSC before using the model-calibration dataset to determine the best-fit model for prediction of SSC from the turbidity and, in some cases, discharge. Correlation between suspended-sand concentration ( SSandC ) in EWI samples and concurrent samples collected by an autosampler was low; therefore, SSandC was excluded from development of surrogate relations because a large part of the calibration dataset was from autosamples. Instead, SSandC was estimated as SSC minus suspended-silt-clay concentration ( SSiltC ). At all sites, the best-fit models included the base-10 logarithm of concentration and turbidity, and at Kansas chute upstream, the base-10 logarithm of streamflow was also included in the best-fit models. These surrogate models were used to estimate continuous time series of SSC and SSiltC . Estimated concentrations of suspended sediment were used to estimate instantaneous and daily loads for total suspended sediment, suspended silt-clay, and suspended sand. Estimated daily suspended-sediment loads were not significantly different between upstream and downstream transects within the Kansas chute, and most individual daily loads within the chute were not significantly different between upstream and downstream transects when evaluated using overlap in daily 95-percent confidence intervals. The comparison of daily load values for upstream and downstream chute transects, as estimated from turbidity-based surrogate models for Kansas chute, documents the daily dynamic nature of sediment transport within the chute with a temporal resolution that is not practical with discrete suspended-sediment sampling alone. Comparisons of concentrations and loads from EWI samples collected from different transects within a study site resulted in few significant differences, but comparisons are limited by small sample sizes and large within-transect variability. When comparing the Missouri River upstream transect to the chute inlet transect, similar results were determined in 2012 as were determined in 2008—the chute inlet affected the amount of sediment entering the chute from the main channel. In addition, the Kansas chute is potentially affecting the sediment concentration within the Missouri River main channel, but small sample size and construction activities within the chute limit the ability to fully understand either the effect of the chute in 2012 or the effect of the chute on the main channel during a year without construction. Finally, some differences in SSC were detected between the Missouri River upstream transects and the chute downstream transects; however, the effect of the chutes on the Missouri River main-channel sediment transport was difficult to isolate because of construction activities and sampling variability.

Nebraska↗

USGS2013

Estimating nitrate concentrations in groundwater at selected wells and springs in the surficial aquifer system and Upper Floridan aquifer, Dougherty Plain and Marianna Lowlands, Georgia, Florida, and Alabama, 2002-50

Groundwater from the surficial aquifer system and Upper Floridan aquifer in the Dougherty Plain and Marianna Lowlands in southwestern Georgia, northwestern Florida, and southeastern Alabama is affected by elevated nitrate concentrations as a result of the vulnerability of the aquifer, irrigation water-supply development, and intensive agricultural land use. The region relies primarily on groundwater from the Upper Floridan aquifer for drinking-water and irrigation supply. Elevated nitrate concentrations in drinking water are a concern because infants under 6 months of age who drink water containing nitrate concentrations above the U.S. Environmental Protection Agency maximum contaminant level of 10 milligrams per liter as nitrogen can become seriously ill with blue baby syndrome. In response to concerns about water quality in domestic wells and in springs in the lower Apalachicola–Chattahoochee–Flint River Basin, the Florida Department of Environmental Protection funded a study in cooperation with the U.S. Geological Survey to examine water quality in groundwater and springs that provide base flow to the Chipola River. A three-dimensional, steady-state, regional-scale groundwater-flow model and two local-scale models were used in conjunction with particle tracking to identify travel times and areas contributing recharge to six groundwater sites—three long-term monitor wells (CP-18A, CP-21A, and RF-41) and three springs (Jackson Blue Spring, Baltzell Springs Group, and Sandbag Spring) in the lower Apalachicola–Chattahoochee–Flint River Basin. Estimated nitrate input to groundwater at land surface, based on previous studies of nitrogen fertilizer sales and atmospheric nitrate deposition data, were used in the advective transport models for the period 2002 to 2050. Nitrate concentrations in groundwater samples collected from the six sites during 1993 to 2007 and groundwater age tracer data were used to calibrate the transport aspect of the simulations. Measured nitrate concentrations (as nitrogen) in wells and springs sampled during the study ranged from 0.37 to 12.73 milligrams per liter. Average apparent ages of groundwater calculated from measurements of chlorofluorocarbon, sulfur hexafluoride, and tritium from wells CP-18A, CP-21A,and RF-41 were about 23, 29, and 32 years, respectively. Average apparent ages of groundwater from Baltzell Springs Group, Sandbag Spring, and Jackson Blue Spring were about 16, 18, and 19 years, respectively. Simulated travel times of particles from the six selected sites ranged from less than 1 day to 511 years; both the minimum and maximum particle travel times were estimated for water from Jackson Blue Spring. Median simulated travel times of particles were about 30, 38, and 62 years for Jackson Blue Spring, Sandbag Spring, and Baltzell Springs Group, respectively. Study results indicated that travel times for approximately 50 percent of the particles from all spring sites were less than 50 years. The median simulated travel times of particles arriving at receptor wells CP-18A, CP-21A, and RF-41 were about 50, 35, and 36 years, respectively. All particle travel times were within the same order of magnitude as the tracer-derived average apparent ages for water, although slightly older than the measured ages. Travel time estimates were substantially greater than the measured age for groundwater reaching well CP-18A, as confirmed by the average apparent age of water determined from tracers. Local-scale particle-tracking models were used to predict nitrate concentrations in the three monitor wells and three springs from 2002 to 2050 for three nitrogen management scenarios: (1) fixed input of nitrate at the 2001 level, (2) reduction of nitrate inputs of 4 percent per year (from the previous year) from 2002 to 2050, and (3) elimination of nitrate input after 2001. Simulated nitrate concentrations in well CP-21A peaked at 7.82 milligrams per liter in 2030, and concentrations in background well RF-41 peaked at 1.10 milligrams per liter in 2020. The simulated particle travel times were longer than indicated by age dating analysis for groundwater in well CP-18A; to account for the poor calibration fit at this well, nitrate concentrations were shifted 21 years. With the shift, simulated nitrate concentrations in groundwater at CP-18A peaked at 13.76 milligrams per liter in 2026. For groundwater in Baltzell Springs Group, Jackson Blue Spring, and Sandbag Spring, simulated nitrate concentrations peaked at 3.77 milligrams per liter in 2006, 3.51 milligrams per liter in 2011, and 0.81 milligram per liter in 2018, respectively, under the three management scenarios. In management scenario 3 (elimination of nitrate input after 2001), simulated nitrate concentrations in Baltzell Springs Group decreased to less than background concentrations (0.10 milligram per liter) by 2033, and in Sandbag Spring concentrations decreased to less than background by 2041. Simulations using nitrate management scenarios 1 (fixed input of nitrate at 2001 levels) and 2 (reduction of 4.0 percent per year) indicate that nitrate concentrations in groundwater may remain above background concentrations through 2050 at all sites.

Alabama, Florida, Georgia↗