Geology ReportsSearch

USGS · 70157048

A case study demonstrating analysis of stormflows, concentrations, and loads of nutrients in highway runoff and swale discharge with the Stochastic Empirical Loading and Dilution Model (SELDM)

Abstract

Decisionmakers need information about the quality and quantity of stormwater runoff, the risk for adverse effects of runoff on receiving waters, and the potential effectiveness of mitigation measures to reduce these risks. The Stochastic Empirical Loading and Dilution Model (SELDM) uses Monte Carlo methods to generate stormflows, concentrations, and loads from a highway site and an upstream basin to provide needed risk-based information. SELDM was designed to help inform water-management decisions for streams and lakes receiving runoff from a highway or other land-use site. The purpose of this paper is to provide a brief description of SELDM and a hypothetical case study demonstrating the type of risk-based information that SELDM can provide. Total nitrogen (TN) and total phosphorus (TP) were selected as example constituents because nutrients are a common concern throughout the Nation and data for receiving waters, highway runoff, and the performance of best management practices (BMPs) are readily available for these constituents. The case study is hypothetical, but was formulated by using actual data from selected monitoring sites in New England. Data representing streamflow and water-quality were collected at U.S. Geological Survey (USGS) streamgage 01208950 Sasco Brook near Southport, CT, which has a drainage area of 7.38 square miles. In this hypothetical case study a 4-lane highway would replace the current 2-lane road and would have a contributing area of 2.2 acres between the topographic basin divides. Concentrations of TN and TP in highway runoff were simulated with data from USGS highway-runoff monitoring station 423027071291301 along State Route 2 in Littleton Massachusetts. Results of a highway-runoff analysis are shown in relation to three hypothetical discharge criteria for TN and two hypothetical discharge criteria for TP. The risks for exceeding TN discharge criteria of 3, 5, and 8 mg/L for highway runoff are 7.4, 0.83, and 0.13 percent of 1,721 runoff events that may occur during a stochastic 30-year simulation. If a grassy swale is used to treat the runoff, the risks for TN exceedances are reduced to 3.2, 0.33 and 0.03 percent, respectively. The risks for exceeding TP discharge criteria of 0.1 and 0.5 mg/L for highway runoff are 49 and 1.2 percent, respectively. If a grassy swale is used to treat the runoff, the risks for TP exceedances are 57 and 0.8 percent, respectively. The risks for the 0.1 mg/L criterion increase because swales can be a source of TP if pavement concentrations are low. The risks for the 0.5 mg/L criterion decrease because the swale is effective for reducing high TP concentrations. Although the results are mixed for storm-event concentrations, the grassy swale effectively reduces annual loads. Annual loads from the swale are, on average, about 49 percent of highway loads for TN and 62 percent of highway loads of TP because the swale reduces high runoff concentrations and stormflow volumes. Analysis of upstream and downstream concentrations indicates that runoff from the site of interest does not have a substantial effect on instream stormflow concentrations in this example simulation.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Gregory E. Granato, Susan C. Jones. 2015. A case study demonstrating analysis of stormflows, concentrations, and loads of nutrients in highway runoff and swale discharge with the Stochastic Empirical Loading and Dilution Model (SELDM). https://pubs.usgs.gov/publication/70157048

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Time-to-depth conversion of seismic-reflection data from eastern Lake Superior and implications for the eastern arm of the Midcontinent Rift

Seismic-reflection data were acquired in the mid 1980s along several lines across eastern Lake Superior by industry and the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE) (Fig. 1). The lines form part of a larger network of crossing lines over the entire lake, which can be used to develop three-dimensional geologic models of the Mesoproterozoic Midcontinent Rift that lies below. To better interpret these lines, we developed velocity models to convert seismic reflections versus two-way travel time (TWTT) to reflections versus depth. In addition, the velocity models themselves provide insights into the structure of the Midcontinent Rift by recognizing common velocity ranges for certain rock types (Grauch, 2023).

eastern Lake Superior

Revisiting the utility of regional-scale, high-quality geophysical data in mineral exploration - A case study featuring the Mammoth Magnetic Anomaly, Pinal County, Arizona

Regional aeromagnetic surveys passively measure the total magnetic intensity (TMI) and are a foundational tool used in mineral exploration (Airo, 2015). With the increased global demand and the number of critical mineral resources required for manufacturing high-tech devices, developing high-quality, regional-scale geophysical surveys could aid critical mineral exploration efforts and geologic mapping. In 2019, the U. S. Geological Survey launched the Earth Mapping Resources Initiative (Earth MRI) to modernize the geologic and geophysical mapping of regions that have the potential to contain critical mineral resources within the United States. In support of planning Earth MRI geophysical surveys, Drenth and Grauch (2019) defined five aeromagnetic data quality rankings (rank 1 through rank 5) applying them to the airborne geophysical survey inventory of the United States (Johnson et al., 2021). Rank 1 aeromagnetic surveys are of the highest quality, meeting modern standards and allowing best practices for qualitative and quantitative interpretation; whereas rank 5 aeromagnetic surveys are of the lowest quality, being useful only for qualitative interpretation of broad features. Through the Earth MRI effort, 48 high-quality, regional-scale rank 1 and 2 airborne magnetic and radiometric geophysical surveys have been planned, collected, or publicly release through May 2025 (U. S. Geological Survey, 2025). Here, a portion of a rank 1 Earth MRI aeromagnetic survey in southeast Arizona is presented and compared to a legacy rank 5 aeromagnetic survey over the Mammoth Magnetic Anomaly (MMA), demonstrating how modern, high-quality aeromagnetic data improves our view of crustal geology, aiding mineral exploration.

Arizona

Early Miocene volcanic rocks and associated tectonics, Lava Hills and southern Bristol Mountains, California

Volcanic rocks of latest Oligocene to early Miocene age form an east-west belt across part of the central eastern Mojave Desert from the Whipple Mountains on the east to the Rosamond Hills on the west. We term this the central belt because it is separated from northern and southern belts by swaths with no volcanic rocks. Limited geochronologic data indicate that much of the belt is latest Oligocene and early Miocene in age, about 24 to 19 Ma, a finding that is consistent with these rocks being overlain by the 18.8 Ma Peach Spring Tuff in many places. We describe Miocene geology in a central area of the belt, in the Lava Hills, southern Bristol Mountains, and southern Old Dad Mountains. Sedimentary basins formed coeval with early andesite to rhyolite volcanism, progressing from fluvial and lacustrine tuffaceous sandstone to volcanic lavas, tuffs, and breccias, indicating that early basins formed proximal to volcanic edifices. Higher strata are fluvial and lacustrine with lavas punctuating the sequence. Although basins may partly have been formed within topographic lows bounded by volcanic domes, plateaus, and stratovolcanoes, consistent stratigraphic sections over wide areas indicate that tectonic basin evolution affected broad areas. The volcanic section is capped by local basalt flows and the regional Peach Spring Tuff. Limited data on normal faults support interpretations of early extensional basin development caused by northeast-southwest oriented stretching. Later extension caused stratal rotations, tilting early deposits down to the southwest. This tilted and subsequently beveled basin architecture was overlain by the youngest volcanic deposits, primarily rhyolite and basalt. The Peach Spring Tuff, 18.8 Ma, lies within this upper unit. Similar stratigraphic and structural relations are exposed in the nearby Marble Mountains and Van Winkle Mountain sections, reinforcing that a broad area underwent similar volcanism and tectonism. In our study area the upper unit is only very gently tilted except near dextral strike-slip faults of the eastern California shear zone. These late Miocene to Recent faults are represented as four main fault zones spaced about 5 km apart, representing distributed shear north of the Bristol Lake basin.

California