Search USGSSearch

SEARCH · Search USGS

Results for “Trends”

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

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

At least 361 records · Page 20Linked to original sources

Assessment of surface water chloride and conductivity trends in areas of unconventional oil and gas development — Why existing national data sets cannot tell us what we would like to know

Heightened concern regarding the potential effects of unconventional oil and gas development on regional water quality has emerged, but the few studies on this topic are limited in geographic scope. Here we evaluate the potential utility of national and publicly available water-quality data sets for addressing questions regarding unconventional oil and gas development. We used existing U.S. Geological Survey and U.S. Environmental Protection Agency data sets to increase understanding of the spatial distribution of unconventional oil and gas development in the U.S. and broadly assess surface water quality trends in these areas. Based on sample size limitations, we were able to estimate trends in specific conductance (SC) and chloride (Cl − ) from 1970 to 2010 in 16% ( n = 155) of the watersheds with unconventional oil and gas resources. We assessed these trends relative to spatiotemporal distributions of hydraulically fractured wells. Results from this limited analysis suggest no consistent and widespread trends in surface water quality for SC and Cl − in areas with increasing unconventional oil and gas development and highlight limitations of existing national databases for addressing questions regarding unconventional oil and gas development and water quality.

Water Resources Research

The recent warming trend in North Greenland

The Arctic is among the fastest warming regions on Earth, but it is also one with limited spatial coverage of multidecadal instrumental surface air temperature measurements. Consequently, atmospheric reanalyses are relatively unconstrained in this region, resulting in a large spread of estimated 30 year recent warming trends, which limits their use to investigate the mechanisms responsible for this trend. Here we present a surface temperature reconstruction over 1982–2011 at NEEM (North Greenland Eemian Ice Drilling Project, 51°W, 77°N), in North Greenland, based on the inversion of borehole temperature and inert gas isotope data. We find that NEEM has warmed by 2.7 ± 0.33°C over the past 30 years, from the long-term 1900–1970 average of −28.55 ± 0.29°C. The warming trend is principally caused by an increase in downward longwave heat flux. Atmospheric reanalyses underestimate this trend by 17%, underlining the need for more in situ observations to validate reanalyses.

Geophysical Research Letters

Factors driving nutrient trends in streams of the Chesapeake Bay watershed

Despite decades of effort toward reducing nitrogen and phosphorus flux to Chesapeake Bay, water-quality and ecological responses in surface waters have been mixed. Recent research, however, provides useful insight into multiple factors complicating the understanding of nutrient trends in bay tributaries, which we review in this paper, as we approach a 2025 total maximum daily load (TMDL) management deadline. Improvements in water quality in many streams are attributable to management actions that reduced point sources and atmospheric nitrogen deposition and to changes in climate. Nutrient reductions expected from management actions, however, have not been fully realized in watershed streams. Nitrogen from urban nonpoint sources has declined, although water-quality responses to urbanization in individual streams vary depending on predevelopment land use. Evolving agriculture, the largest watershed source of nutrients, has likely contributed to local nutrient trends but has not affected substantial changes in flux to the bay. Changing average nitrogen yields from farmland underlain by carbonate rocks, however, may suggest future trends in other areas under similar management, climatic, or other influences, although drivers of these changes remain unclear. Regardless of upstream trends, phosphorus flux to the bay from its largest tributary has increased due to sediment infill in the Conowingo Reservoir. In general, recent research emphasizes the utility of input reductions over attempts to manage nutrient fate and transport at limiting nutrients in surface waters. Ongoing research opportunities include evaluating effects of climate change and conservation practices over time and space and developing tools to disentangle and evaluate multiple influences on regional water quality.

Chesapeake Bay watershed

Strong variation in Brook Trout trends across geology, elevation, and stream size in Shenandoah National Park

Objective Landscape context structures fish abundance and dynamics, and understanding trends in fish abundance across the landscape is often prerequisite for effective conservation. In this study, we evaluated the status and trends of Brook Trout Salvelinus fontinalis in Shenandoah National Park to understand how these are structured across bedrock geology, elevation, and stream size. Methods We used long-term monitoring data from 94 sites in Shenandoah National Park to evaluate trends in Brook Trout abundance over a 27-year period (1996–2022) and assess the importance of local environmental covariates using a hierarchical Bayesian N-mixture model based on depletion sampling. Focal covariates were chosen for their demonstrated importance in structuring fish populations in Shenandoah National Park and elsewhere. Bedrock geology controls sensitivity to acid deposition, watershed area is related to stream habitat features such as complexity and flow variability, and elevation creates gradients in temperature. Result Models revealed significant decreases in adult Brook Trout abundance over time (95% credible intervals < 0) for 31 of 94 sites (33%), and at least three sites exhibited apparent extirpations over the study period. Estimated Brook Trout abundance declined by 50% or more in approximately 70% of streams across the park over the study period. Sites with the warmest water temperatures exhibited the fastest declines in abundance. However, large watersheds on poorly buffered bedrock exhibited significant gains in abundance over time, suggesting some recovery from acid deposition due to improvements in air quality. Conclusion Our analysis revealed large and divergent changes in Brook Trout abundance over recent decades and suggests the importance of local water temperature and acid sensitivity as probable causal mechanisms. These results highlight the importance of considering local factors when evaluating long-term trends in stream fish populations. Results of this study can assist the development of targeted conservation actions within Shenandoah National Park and elsewhere.

Virginia

Historical cover trends in a sagebrush steppe ecosystem from 1985 to 2013: Links with climate, disturbance, and management

Understanding the causes and consequences of component change in sagebrush steppe is crucial for evaluating ecosystem sustainability. The sagebrush ( Artemisia spp.) steppe ecosystem of the northwest USA has been impacted by the invasion of exotic grasses, increasing fire return intervals, changing land management practices, and fragmentation, often lowering the overall resilience to change. We utilized contemporary and historical Landsat imagery, field data, and regression tree models to produce fractional cover maps of rangeland components (shrub, sagebrush, herbaceous, bare ground, and litter) through the last 30 years. Our main goals were to (1) investigate rangeland component trends over 30 years, (2) evaluate the magnitude and direction of trends in components and climate drivers and their relationship, and (3) assess component trends influenced by climate. Results indicated that over the study period, shrub, sage, herbaceous, and litter cover decreased, while bare ground cover increased. Measured rates of change ranged from − 0.14% decade −1 for shrub cover to 0.05% decade −1 for bare ground, whereas herbaceous and litter cover trends were negligible. Net landscape cover changes were consistent with expectations of climate change and disturbance producing a loss of biotic cover, and converting a portion of shrub and sagebrush to herbaceous cover. Overall, fire and related successional recovery was the greatest change agent for all components in terms of area and cover change, while increasing minimum temperature, at a rate of 0.66°C decade −1 , was found to be the most significant climate driver.

California, Nevada, Oregon

Extreme events, trends, and variability in Northern Hemisphere lake-ice phenology (1855-2005)

Often extreme events, more than changes in mean conditions, have the greatest impact on the environment and human well-being. Here we examine changes in the occurrence of extremes in the timing of the annual formation and disappearance of lake ice in the Northern Hemisphere. Both changes in the mean condition and in variability around the mean condition can alter the probability of extreme events. Using long-term ice phenology data covering two periods 1855–6 to 2004–5 and 1905–6 to 2004–5 for a total of 75 lakes, we examined patterns in long-term trends and variability in the context of understanding the occurrence of extreme events. We also examined patterns in trends for a 30-year subset (1975–6 to 2004–5) of the 100-year data set. Trends for ice variables in the recent 30-year period were steeper than those in the 100- and 150-year periods, and trends in the 150-year period were steeper than in the 100-year period. Ranges of rates of change (days per decade) among time periods based on linear regression were 0.3−1.6 later for freeze, 0.5−1.9 earlier for breakup, and 0.7−4.3 shorter for duration. Mostly, standard deviation did not change, or it decreased in the 150-year and 100-year periods. During the recent 50-year period, standard deviation calculated in 10-year windows increased for all ice measures. For the 150-year and 100-year periods changes in the mean ice dates rather than changes in variability most strongly influenced the significant increases in the frequency of extreme lake ice events associated with warmer conditions and decreases in the frequency of extreme events associated with cooler conditions.

Climatic Change

Temporal trends in algae, benthic invertebrate, and fish assemblages in streams and rivers draining basins of varying land use in the south-central United States, 1993-2007

Site-specific temporal trends in algae, benthic invertebrate, and fish assemblages were investigated in 15 streams and rivers draining basins of varying land use in the south-central United States from 1993–2007. A multivariate approach was used to identify sites with statistically significant trends in aquatic assemblages which were then tested for correlations with assemblage metrics and abiotic environmental variables (climate, water quality, streamflow, and physical habitat). Significant temporal trends in one or more of the aquatic assemblages were identified at more than half (eight of 15) of the streams in the study. Assemblage metrics and abiotic environmental variables found to be significantly correlated with aquatic assemblages differed between land use categories. For example, algal assemblages at undeveloped sites were associated with physical habitat, while algal assemblages at more anthropogenically altered sites (agricultural and urban) were associated with nutrient and streamflow metrics. In urban stream sites results indicate that streamflow metrics may act as important controls on water quality conditions, as represented by aquatic assemblage metrics. The site-specific identification of biotic trends and abiotic–biotic relations presented here will provide valuable information that can inform interpretation of continued monitoring data and the design of future studies. In addition, the subsets of abiotic variables identified as potentially important drivers of change in aquatic assemblages provide policy makers and resource managers with information that will assist in the design and implementation of monitoring programs aimed at the protection of aquatic resources.

Arkansas, Colorado, Kansas, Kentucky, Louisiana, M

Evaluation of procedures for prediction of unconventional gas in the presence of geologic trends

This study extends the application of local spatial nonparametric prediction models to the estimation of recoverable gas volumes in continuous-type gas plays to regimes where there is a single geologic trend. A transformation is presented, originally proposed by Tomczak, that offsets the distortions caused by the trend. This article reports on numerical experiments that compare predictive and classification performance of the local nonparametric prediction models based on the transformation with models based on Euclidean distance. The transformation offers improvement in average root mean square error when the trend is not severely misspecified. Because of the local nature of the models, even those based on Euclidean distance in the presence of trends are reasonably robust. The tests based on other model performance metrics such as prediction error associated with the high-grade tracts and the ability of the models to identify sites with the largest gas volumes also demonstrate the robustness of both local modeling approaches. ?? International Association for Mathematical Geology 2009.

Natural Resources Research

The long-term trends (1982-2006) in vegetation greenness of the alpine ecosystem in the Qinghai-Tibetan Plateau

The increased rate of annual temperature in the Qinghai-Tibetan Plateau exceeded all other areas of the same latitude in recent decades. The influence of the warming climate on the alpine ecosystem of the plateau was distinct. An analysis of alpine vegetation under changes in climatic conditions was conducted in this study. This was done through an examination of vegetation greenness and its relationship with climate variability using the Advanced Very High Resolution Radiometer satellite imagery and climate datasets. Vegetation in the plateau experienced a positive trend in greenness, with 18.0 % of the vegetated areas exhibiting significantly positive trends, which were primarily located in the eastern and southwestern parts of the plateau. In grasslands, 25.8 % of meadows and 14.1 % of steppes exhibited significant upward trends. In contrast, the broadleaf forests experienced a trend of degradation. Temperature, particularly summer temperature, was the primary factor promoting the vegetation growth in the plateau. The wetter and warmer climate in the east contributed to the favorable conditions for vegetation. The alpine meadow was mostly sensitive to temperature, while the steppes were sensitive to both temperature and precipitation. Although a warming climate was expected to be beneficial to vegetation growth in the alpine region, the rising temperature coupled with reduced precipitation in the south did not favor vegetation growth due to low humidity and poor soil moisture conditions.

Qinghai-tibetan Plateau

Strontium isotopic constraints on the origin of ore-forming fluids of the Viburnum Trend, southeast Missouri

We have measured 87 Sr 86 Sr "> 87Sr86Sr and Rb and Sr concentrations in several minerals, primarily sulfides, spanning the paragenesis of hydrothermal mineralization in the Viburnum Trend in southeast Missouri. Separate measurements were made for fluid inclusions opened by crushing or thermal decrepitation and for the solids. For comparison, measurements were also made on samples of probable local aquifers, the Bonneterre Formation and the Lamotte Sandstone. For some of the samples, concentrations of K, Ca, Cl, Na, and Mg are also reported. In several cases 87 Sr 86 Sr ratios "> 87Sr86Srratios are different (higher) in the solids than in the fluid inclusions. We have investigated the possibility that either type of sample gives spurious results, e.g., that the fluid inclusions are secondary or contaminated by host dolomite, or that Sr in the solids reflects a detrital rather than an authigenic source. Consideration of mass balance, overall solute chemistry, and examination of non-sulfide dissolution residue, however, suggests that both types of sample reflect primary fluid Sr composition. We thus adopt the working hypothesis that Sr isotopic composition in fluids at the time of hydrothermal mineralization was highly variable. The observed results for the Viburnum Trend do not conform well to expected trends for Mississippi Valley-type (MVT) mineralization based primarily on analyses of gangue carbonates, barite, and fluorite. Fluid inclusion Sr in some of the Viburnum Trend samples is more radiogenic than in the host dolomite but only moderately so; in other samples, notably main-stage octahedral galena, fluid inclusion Sr composition is within the range observed for the host dolomite. In contrast, Sr in some of the sulfides is very radiogenic, much more so than previously reported for MVT minerals, and is very radiogenic early in the paragenesis (pyrite, chalcopyrite, sphalerite), less radiogenic during main-stage ore deposition, and again more radiogenic in later paragenetic stages.

Geochimica et Cosmochimica Acta

Secular trends in the geologic record and the supercontinent cycle

Geologic secular trends are used to refine the timetable of supercontinent assembly, tenure, and breakup. The analysis rests on what is meant by the term supercontinent, which here is defined broadly as a grouping of formerly dispersed continents. To avoid the artificial pitfall of an all-or-nothing definition, quantitative measures of “supercontinentality” are presented: the number of continents, and the area of the largest continent, which both can be gleaned from global paleogeographic maps for the Phanerozoic. For the secular trends approach to be viable in the deep past when the very existence of supercontinents is debatable and reconstructions are fraught with problems, it must first be calibrated in the Phanerozoic against the well-constrained Pangea supercontinent cycle. The most informative geologic variables covering both the Phanerozoic and Precambrian are the abundances of passive margins and of detrital zircons. Both fluctuated with size of the largest continent during the Pangea supercontinent cycle and can be quantified back to the Neoarchean. The tenure of Pangea was a time represented in the rock record by few zircons and few passive margins. Thus, previously documented minima in the abundance of detrital zircons (and orogenic granites) during the Precambrian (Condie et al., 2009a, Gondwana Research 15, 228–242) now can be more confidently interpreted as marking the tenures of supercontinents. The occurrences of carbonatites, granulites, eclogites, and greenstone-belt deformation events also appear to bear the imprint of Precambrian supercontinent cyclicity. Together, these secular records are consistent with the following scenario. The Neoarchean continental assemblies of Superia and Sclavia broke up at ca. 2300 and ca. 2090 Ma, respectively. Some of their fragments collided to form Nuna by about 1750 Ma; Nuna then grew by lateral accretion of juvenile arcs during the Mesoproterozoic, and was involved in a series of collisions at ca. 1000 Ma to form Rodinia. Rodinia broke up in stages from ca. 1000 to ca. 520 Ma. Before Rodinia had completely come apart, some of its pieces had already been reassembled in a new configuration, Gondwana, which was completed by 530 Ma. Gondwana later collided with Laurentia, Baltica, and Siberia to form Pangea by about 300 Ma. Breakup of Pangea began at about 180 Ma (Early Jurassic) and continues today. In the suggested scenario, no supercontinent cycle in Earth history corresponded to the ideal, in which all the continents were gathered together, then broke apart, then reassembled in a new configuration. Nuna and Gondwana ended their tenures not by breakup but by collision and name change; Rodinia's assembly overlapped in time with its disassembly; and Pangea spalled Tethyan microcontinents throughout much of its tenure. Many other secular trends show a weak or uneven imprint of the supercontinent cycle, no imprint at all. Instead, these secular trends together reveal aspects of the shifting background against which the supercontinents came and went, making each cycle unique. Global heat production declined; plate tectonics sped up through the Proterozoic and slowed down through the Phanerozoic; the atmosphere and oceans became oxidized; life emerged as a major geochemical agent; some rock types went extinct or nearly so (BIF, massif-type anorthosite, komatiite); and other rock types came into existence or became common (blueschists, bioclastic limestone, coal).

Earth-Science Reviews

Long-term trends of surface-water mercury and methylmercury concentrations downstream of historic mining within the Carson River watershed

The Carson River is a vital water resource for local municipalities and migratory birds travelling the Pacific Flyway. Historic mining practices that used mercury (Hg) to extract gold from Comstock Lode ore has left much of the river system heavily contaminated with Hg, a practice that continues in many parts of the world today. Between 1998 and 2013, the United States Geological Survey (USGS) collected and analyzed Carson River water for Hg and methylmercury (MeHg) concentrations resulting in a sixteen year record of unfiltered total mercury (uf.THg), filtered (dissolved) Hg (f.THg), total methylmercury (uf.MeHg), filtered MeHg (f.MeHg), and particulate-bound THg (p.THg) and MeHg (p.MeHg) concentrations. This represents one of the longest continuous records of Hg speciation data for any riverine system, thereby providing a unique opportunity to evaluate long-term trends in concentrations and annual loads. During the period of analysis, uf.THg concentration and load trended downward at rates of −0.85% and −1.8% per year, respectively. Conversely, the f.THg concentration increased at a rate of 1.7% per year between 1998 and 2005, and 4.9% per year between 2005 and 2013. Trends in flow-normalized partition coefficients for both Hg and MeHg suggest a statistically significant shift from the particulate to the filtered phase. The upwardly accelerating f.THg concentration and observed shift from the solid phase to the aqueous phase among the pools of Hg and MeHg within the river water column signals an increased risk of deteriorating ecological conditions in the lower basin with respect to Hg contamination. More broadly, the 16-year trend analysis, completed 140 years after the commencement of major Hg releases to the Carson River, provides a poignant example of the ongoing legacy left behind by gold and silver mining techniques that relied on Hg amalgamation, and a cautionary tale for regions still pursuing the practice in other countries.

Nevada

A bootstrap method for estimating uncertainty of water quality trends

Estimation of the direction and magnitude of trends in surface water quality remains a problem of great scientific and practical interest. The Weighted Regressions on Time, Discharge, and Season (WRTDS) method was recently introduced as an exploratory data analysis tool to provide flexible and robust estimates of water quality trends. This paper enhances the WRTDS method through the introduction of the WRTDS Bootstrap Test (WBT), an extension of WRTDS that quantifies the uncertainty in WRTDS-estimates of water quality trends and offers various ways to visualize and communicate these uncertainties. Monte Carlo experiments are applied to estimate the Type I error probabilities for this method. WBT is compared to other water-quality trend-testing methods appropriate for data sets of one to three decades in length with sampling frequencies of 6–24 observations per year. The software to conduct the test is in the EGRETci R-package.

Environmental Modelling and Software

The inevitability of large shallow craters on Callisto and Ganymede: Implications for crater depth-diameter trends

Complex craters with diameters ( D ) ≥ 40 km on Callisto and Ganymede are shallower than would be expected from simply extrapolating the depth-diameter trend from smaller ( D ≤ 40 km) craters. This unusual depth-diameter ( d-D ) trend, and associated changes in crater morphology, have been hypothesized to result from rheological transitions, including the existence of an ocean, within the moons' ice shell. Simulations of impact crater formation can reproduce the observed shallow depths but require heat fluxes roughly twice the maximum radiogenic flux to do so. Here we demonstrate that the d-D trends on Callisto and Ganymede can instead be explained as a direct consequence of viscous relaxation under radiogenic heating. We use numerical simulations of viscous relaxation to show that if craters form at the depth expected from an extrapolation of the complex crater d-D trend, they will evolve to the observed depths over timescales of 200 Myrs to 1 Gyrs. Large craters (e.g., D ≥ 80 km) younger than 200 Myrs, which would retain greater depths, should be relatively rare. If we instead assume that the craters formed at their observed depths, as proposed by previous impact modeling, they quickly become much shallower than observed. We find excellent agreement between observed crater depths on Ganymede and our simulated crater depths by assuming a pure-water ice composition and a diurnally averaged surface temperature of 120 K, but require either larger-grained or “dirty” ice with a modestly higher viscosity to match observations at Callisto, where the surface temperature is warmer (130 K). We favor the latter explanation because it is consistent with the existence of a dusty lag on Callisto's surface and the absence of a similar lag on Ganymede. Our results predict that, for a given crater diameter, post-relaxation crater depth should increase with increasing latitude, a hypothesis best tested on Callisto, whose relatively quiescent geologic history best preserves the signature of viscous relaxation under radiogenic heating.

Icarus

Trends in streamflow in the Yukon River Basin from 1944 to 2005 and the influence of the Pacific Decadal Oscillation

Streamflow characteristics in the Yukon River Basin of Alaska and Canada have changed from 1944 to 2005, and some of the change can be attributed to the two most recent modes of the Pacific Decadal Oscillation (PDO). Seasonal, monthly, and annual stream discharge data from 21 stations in the Yukon River Basin were analyzed for trends over the entire period of record, generally spanning 4-6 decades, and examined for differences between the two most recent modes of the PDO: cold-PDO (1944-1975) and warm-PDO (1976-2005) subsets. Between 1944 and 2005, average winter and April flow increased at 15 sites. Observed winter flow increases during the cold-PDO phase were generally limited to sites in the Upper Yukon River Basin. Positive trends in winter flow during the warm-PDO phase broadened to include stations in the Middle and Lower Yukon River drainage basins. Increases in winter streamflow most likely result from groundwater input enhanced by permafrost thawing that promotes infiltration and deeper subsurface flow paths. Increased April flow may be attributed to a combination of greater baseflow (from groundwater increases), earlier spring snowmelt and runoff, and increased winter precipitation, depending on location. Calculated deviations from long-term mean monthly discharges indicate below-average flow in the winter months during the cold PDO and above-average flow in the winter months during the warm PDO. Although not as strong a signal, results also support the reverse response during the summer months: above-average flow during the cold PDO and below-average flow during the warm PDO. Changes in the summer flows are likely an indirect consequence of the PDO, resulting from earlier spring snowmelt runoff and also perhaps increased summer infiltration and storage in a deeper active layer. Annual discharge has remained relatively unchanged in the Yukon River Basin, but a few glacier-fed rivers demonstrate positive trends, which can be attributed to enhanced glacier melting. A positive trend in annual flow during the warm PDO near the mouth of the Yukon River suggests that small increases in flow throughout the Yukon River Basin have resulted in an additive effect manifested in the downstream-most streamflow station. Many of the identified changes in streamflow patterns in the Yukon River Basin show a correlation to the PDO regime shift. This work highlights the importance of considering proximate climate forcings as well as global climate change when assessing hydrologic changes in the Arctic.

Journal of Hydrology

Long-term trends in abundance and potential drivers for eight species of coastal birds in the U.S. South Atlantic

The U.S. South Atlantic coastal region is used by many marine birds for foraging, reproduction, and migration. We developed standardized indices of relative abundance from long–term (1980–2016), semi-structured monitoring data (eBird) for eight species: Brown Pelican ( Pelecanus occidentalis ), Double-Crested Cormorant ( Nannopterum auritum ), White Ibis ( Eudocimus albus ), Wood Stork ( Mycteria americana ), Piping Plover ( Charadrius melodus ), American Oystercatcher ( Haematopus palliatus ), Clapper Rail ( Rallus crepitans ), and Northern Gannet ( Morus bassanus ). Following a period of stable or declining abundance from the 1980s through the 1990s, most species have shown stable or slightly upward trends through the late 2000s; Brown Pelican and Piping Plover have shown some evidence of recent declines. Species–specific correlations between abundance indices developed from presence/absence data and those developed from count data were positive for all species and ranged from 0.53 to 0.86. Dynamic factor analysis identified common trends in abundance among several species, in particular, Brown Pelican, Double–Crested Cormorant, and White Ibis. Model performance was improved with inclusion of an indicator of sea level rise, but not forage fish abundance or temperature, indicating habitat availability mediated by changing water levels may explain some of the underlying abundance trends. Our results provide baseline information on long–term trends for several important coastal birds that can help inform research, monitoring and conservation efforts in the U.S. South Atlantic region.

Florida, Georgia, North Carolina, South Carolina

An approach for decomposing river water-quality trends into different flow classes

A number of statistical approaches have been developed to quantify the overall trend in river water quality, but most approaches are not intended for reporting separate trends for different flow conditions. We propose an approach called FN 2Q , which is an extension of the flow-normalization (FN) procedure of the well-established WRTDS (“Weighted Regressions on Time, Discharge, and Season”) method. The FN 2Q approach provides a daily time series of low-flow and high-flow FN flux estimates that represent the lower and upper half of daily riverflow observations that occurred on each calendar day across the period of record. These daily estimates can be summarized into any time period of interest (e.g., monthly, seasonal, or annual) for quantifying trends. The proposed approach is illustrated with an application to a record of total nitrogen concentration (632 samples) collected between 1985 and 2018 from the South Fork Shenandoah River at Front Royal, Virginia (USA). Results show that the overall FN flux of total nitrogen has declined in the period of 1985–2018, which is mainly attributable to FN flux decline in the low-flow class. Furthermore, the decline in the low-flow class was highly correlated with wastewater effluent loads, indicating that the upgrades of treatment technology at wastewater treatment facilities have likely led to water-quality improvement under low-flow conditions. The high-flow FN flux showed a spike around 2007, which was likely caused by increased delivery of particulate nitrogen associated with sediment transport. The case study demonstrates the utility of the FN 2Q approach toward not only characterizing the changes in river water quality but also guiding the direction of additional analysis for capturing the underlying drivers. The FN 2Q approach (and the published code) can easily be applied to widely available river monitoring records to quantify water-quality trends under different flow conditions to enhance understanding of river water-quality dynamics.

Virginia

The statistical power to detect regional temporal trends in riverine contaminants in the Chesapeake Bay Watershed, USA

Chemical contamination of riverine ecosystems is largely a result of urbanization, industrialization, and agricultural activities occurring on adjacent terrestrial landscapes. Land management activities (e.g., Best Management Practices) are an important tool used to reduce point and non-point sources of pollution. However, the ability to confidently make inferences about the efficacy of land management activities on reducing in-stream chemical concentrations is poorly understood. We estimated regional temporal trends and components of variation for commonly used herbicides (atrazine and metolachlor), total estrogenicity, and riverine sediment concentrations of total PCBs for rivers in the Chesapeake Bay Watershed, USA. We then used the estimated variance components to perform a power analysis and evaluated the statistical power to detect regional temporal trends under different monitoring scenarios. Scenarios included varying the magnitude of the annual contaminant decline, the number of sites sampled each year, the number of years sampled, and sampling frequency. Monitoring for short time periods (e.g., 5 years) was inadequate for detecting regional temporal trends, regardless of the number of sites sampled or the magnitude of the annual declines. Even when monitoring over a 20-year period, sampling a relatively large number of sites each year was required (e.g., > 50 sites) to achieve adequate statistical power for smaller trend magnitudes (declines of 5 – 7%/year). Annual sampling frequency had little impact on power for any monitoring scenario. All sampling scenarios were underpowered for sediment total PCBs. Power was greatest for total estrogenicity, suggesting that this aggregate measure of estrogenic activity may be a useful indicator. This study provides information that can be used to help (1) guide the development of monitoring programs aimed at detecting regional declines in riverine chemical contaminant concentrations in response to land management actions, and (2) set expectations for the ability to detect changes over time.

Chesapeake Bay watershed