Quantifying soil accumulation of atmospheric mercury using fallout radionuclide chronometry.

Soils are a principal global reservoir of mercury (Hg), a neurotoxic pollutant that is accumulating through anthropogenic emissions to the atmosphere and subsequent deposition to terrestrial ecosystems. The fate of Hg in global soils remains uncertain, however, particularly to what degree Hg is re-emitted back to the atmosphere as gaseous elemental mercury (GEM). Here we use fallout radionuclide (FRN) chronometry to directly measure Hg accumulation rates in soils. By comparing these rates with measured atmospheric fluxes in a mass balance approach, we show that representative Arctic, boreal, temperate, and tropical soils are quantitatively efficient at retaining anthropogenic Hg. Potential for significant GEM re-emission appears limited to a minority of coniferous soils, calling into question global models that assume strong re-emission of legacy Hg from soils. FRN chronometry poses a powerful tool to reconstruct terrestrial Hg accumulation across larger spatial scales than previously possible, while offering insights into the susceptibility of Hg mobilization from different soil environments.

Forest growth responds more to air pollution than soil acidification.

The forests of central Europe have undergone remarkable transitions in the past 40 years as air quality has improved dramatically. Retrospective analysis of Norway spruce (Picea abies) tree rings in the Czech Republic shows that air pollution (e.g. SO2 concentrations, high acidic deposition to the forest canopy) plays a dominant role in driving forest health. Extensive soil acidification occurred in the highly polluted "Black Triangle" in Central Europe, and upper mineral soils are still acidified. In contrast, acidic atmospheric deposition declined by 80% and atmospheric SO2 concentration by 90% between the late 1980s and 2010s. In this study we oserved that annual tree ring width (TRW) declined in the 1970s and subsequently recovered in the 1990s, tracking SO2 concentrations closely. Furthermore, recovery of TRW was similar in unlimed and limed stands. Despite large increases in soil base saturation, as well as soil pH, as a result of repeated liming starting in 1981, TRW growth was similar in limed and unlimed plots. TRW recovery was interrupted in 1996 when highly acidic rime (originating from more pronounced decline of alkaline dust than SO2 from local power plants) injured the spruce canopy, but recovered soon to the pre-episode growth. Across the long-term site history, changes in soil chemistry (pH, base saturation, Bc/Al soil solution ratio) cannot explain observed changes in TRW at the two study sites where we tracked soil chemistry. Instead, statistically significant recovery in TRW is linked to the trajectory of annual SO2 concentrations or sulfur deposition at all three stands.

Indirect nitrous oxide emission factors of fluvial networks can be predicted by dissolved organic carbon and nitrate from local to global scales.

Streams and rivers are important sources of nitrous oxide (N2 O), a powerful greenhouse gas. Estimating global riverine N2 O emissions is critical for the assessment of anthropogenic N2 O emission inventories. The indirect N2 O emission factor (EF5r ) model, one of the bottom-up approaches, adopts a fixed EF5r value to estimate riverine N2 O emissions based on IPCC methodology. However, the estimates have considerable uncertainty due to the large spatiotemporal variations in EF5r values. Factors regulating EF5r are poorly understood at the global scale. Here, we combine 4-year in situ observations across rivers of different land use types in China, with a global meta-analysis over six continents, to explore the spatiotemporal variations and controls on EF5r values. Our results show that the EF5r values in China and other regions with high N loads are lower than those for regions with lower N loads. Although the global mean EF5r value is comparable to the IPCC default value, the global EF5r values are highly skewed with large variations, indicating that adopting region-specific EF5r values rather than revising the fixed default value is more appropriate for the estimation of regional and global riverine N2 O emissions. The ratio of dissolved organic carbon to nitrate (DOC/NO3 - ) and NO3 - concentration are identified as the dominant predictors of region-specific EF5r values at both regional and global scales because stoichiometry and nutrients strictly regulate denitrification and N2 O production efficiency in rivers. A multiple linear regression model using DOC/NO3 - and NO3 - is proposed to predict region-specific EF5r values. The good fit of the model associated with easily obtained water quality variables allows its widespread application. This study fills a key knowledge gap in predicting region-specific EF5r values at the global scale and provides a pathway to estimate global riverine N2 O emissions more accurately based on IPCC methodology.

Linking soil phosphorus with forest litterfall resistance and resilience to cyclone disturbance: A pantropical meta-analysis.

While tropical cyclone regimes are shifting with climate change, the mechanisms underpinning the resistance (ability to withstand disturbance-induced change) and resilience (capacity to return to pre-disturbance reference) of tropical forest litterfall to cyclones remain largely unexplored pantropically. Single-site studies in Australia and Hawaii suggest that litterfall on low-phosphorus (P) soils is more resistant and less resilient to cyclones. We conducted a meta-analysis to investigate the pantropical importance of total soil P in mediating forest litterfall resistance and resilience to 22 tropical cyclones. We evaluated cyclone-induced and post-cyclone litterfall mass (g/m2 /day), and P and nitrogen (N) fluxes (mg/m2 /day) and concentrations (mg/g), all indicators of ecosystem function and essential for nutrient cycling. Across 73 case studies in Australia, Guadeloupe, Hawaii, Mexico, Puerto Rico, and Taiwan, total litterfall mass flux increased from ~2.5 ± 0.3 to 22.5 ± 3 g/m2 /day due to cyclones, with large variation among studies. Litterfall P and N fluxes post-cyclone represented ~5% and 10% of the average annual fluxes, respectively. Post-cyclone leaf litterfall N and P concentrations were 21.6 ± 1.2% and 58.6 ± 2.3% higher than pre-cyclone means. Mixed-effects models determined that soil P negatively moderated the pantropical litterfall resistance to cyclones, with a 100 mg P/kg increase in soil P corresponding to a 32% to 38% decrease in resistance. Based on 33% of the resistance case studies, total litterfall mass flux reached pre-disturbance levels within one-year post-disturbance. A GAMM indicated that soil P, gale wind duration and time post-cyclone jointly moderate the short-term resilience of total litterfall, with the nature of the relationship between resilience and soil P contingent on time and wind duration. Across pantropical forests observed to date, our results indicate that litterfall resistance and resilience in the face of intensifying cyclones will be partially determined by total soil P.

A general pattern of trade-offs between ecosystem resistance and resilience to tropical cyclones.

Tropical cyclones drive coastal ecosystem dynamics, and their frequency, intensity, and spatial distribution are predicted to shift with climate change. Patterns of resistance and resilience were synthesized for 4138 ecosystem time series from n = 26 storms occurring between 1985 and 2018 in the Northern Hemisphere to predict how coastal ecosystems will respond to future disturbance regimes. Data were grouped by ecosystems (fresh water, salt water, terrestrial, and wetland) and response categories (biogeochemistry, hydrography, mobile biota, sedentary fauna, and vascular plants). We observed a repeated pattern of trade-offs between resistance and resilience across analyses. These patterns are likely the outcomes of evolutionary adaptation, they conform to disturbance theories, and they indicate that consistent rules may govern ecosystem susceptibility to tropical cyclones.

Light and flow regimes regulate the metabolism of rivers.

Mean annual temperature and mean annual precipitation drive much of the variation in productivity across Earth's terrestrial ecosystems but do not explain variation in gross primary productivity (GPP) or ecosystem respiration (ER) in flowing waters. We document substantial variation in the magnitude and seasonality of GPP and ER across 222 US rivers. In contrast to their terrestrial counterparts, most river ecosystems respire far more carbon than they fix and have less pronounced and consistent seasonality in their metabolic rates. We find that variation in annual solar energy inputs and stability of flows are the primary drivers of GPP and ER across rivers. A classification schema based on these drivers advances river science and informs management.

Land Use Overrides Stream Order and Season in Driving Dissolved Organic Matter Dynamics Throughout the Year in a River Network.

Anthropogenic land use has increased nutrient concentrations and altered dissolved organic matter (DOM) character and its bioavailability. Despite widespread recognition that DOM character and its reactivity can vary temporally, the relative influence of land use and stream order on DOM characteristics is poorly understood across seasons and the entire flow regime. We examined DOM character and 28-day bioavailable dissolved organic carbon (BDOC) across a river network to determine the relative roles of land use and stream order in driving variability in DOM character and bioavailability throughout the year. DOM in 1st-order streams was distinct from higher stream orders with lower DOC concentrations, less aromatic (specific ultraviolet absorbance at 254 nm (SUVA254)), more autochthonous (fluorescence index), and more recently produced (ß/α) DOM. Across all months, variability in DOM character was primarily explained by land use, rather than stream order or season. Land use and stream order explained the most DOM variation in transitional and winter months and the least during dry months. BDOC was greater in watersheds with less aromatic (SUVA254) and more recent allochthonous DOM (ß/α) and more development and impervious surface. With continued development, the bioavailability of DOM in the smallest and most impacted watersheds is expected to increase.

Shifting stoichiometry: Long-term trends in stream-dissolved organic matter reveal altered C:N ratios due to history of atmospheric acid deposition.

Dissolved organic carbon (DOC) and nitrogen (DON) are important energy and nutrient sources for aquatic ecosystems. In many northern temperate, freshwater systems DOC has increased in the past 50 years. Less is known about how changes in DOC may vary across latitudes, and whether changes in DON track those of DOC. Here, we present long-term DOC and DON data from 74 streams distributed across seven sites in biomes ranging from the tropics to northern boreal forests with varying histories of atmospheric acid deposition. For each stream, we examined the temporal trends of DOC and DON concentrations and DOC:DON molar ratios. While some sites displayed consistent positive or negative trends in stream DOC and DON concentrations, changes in direction or magnitude were inconsistent at regional or local scales. DON trends did not always track those of DOC, though DOC:DON ratios increased over time for ~30% of streams. Our results indicate that the dissolved organic matter (DOM) pool is experiencing fundamental changes due to the recovery from atmospheric acid deposition. Changes in DOC:DON stoichiometry point to a shifting energy-nutrient balance in many aquatic ecosystems. Sustained changes in the character of DOM can have major implications for stream metabolism, biogeochemical processes, food webs, and drinking water quality (including disinfection by-products). Understanding regional and global variation in DOC and DON concentrations is important for developing realistic models and watershed management protocols to effectively target mitigation efforts aimed at bringing DOM flux and nutrient enrichment under control.

Distinctive Patterns and Controls of Nitrous Oxide Concentrations and Fluxes from Urban Inland Waters.

Inland waters are significant sources of nitrous oxide (N2O), a powerful greenhouse gas. However, considerable uncertainty exists in the estimates of N2O efflux from global inland waters due to a lack of direct measurements in urban inland waters, which are generally characterized by high carbon and nitrogen concentrations and low carbon-to-nitrogen ratios. Herein, we present direct measurements of N2O concentrations and fluxes in lakes and rivers of Beijing, China, during 2018-2020. N2O concentrations and fluxes in the waters of Beijing exceeded previous estimates of global rivers due to the high carbon and nutrient concentrations and high aquatic productivity. In contrast, the N2O emission factor (N2O-N/DIN, median 0.0005) was lower than global medians and the N2O yield (ΔN2O/(ΔN2O + ΔN2), average 1.6%) was higher than those typically observed in rivers and streams. The positive relationship between N2O emissions and denitrifying bacteria as well as the Michaelis-Menten relationship between N2O emissions and NO3--N concentrations suggested that bacteria control the net production of N2O in waters of Beijing with N saturation, leading to a low N2O emission factor. However, low carbon-to-nitrogen ratios are beneficial for N2O accumulation during denitrification, resulting in high N2O yields. This study demonstrates the significant N2O emissions and their distinctive patterns and controls in urban inland waters and suggests that N2O emission estimates based on nitrogen loads and simple emission factor values are not appropriate for urban inland water systems.

Experimental nitrogen and phosphorus enrichment stimulates multiple trophic levels of algal and detrital-based food webs: a global meta-analysis from streams and rivers.

Anthropogenic increases in nitrogen (N) and phosphorus (P) concentrations can strongly influence the structure and function of ecosystems. Even though lotic ecosystems receive cumulative inputs of nutrients applied to and deposited on land, no comprehensive assessment has quantified nutrient-enrichment effects within streams and rivers. We conducted a meta-analysis of published studies that experimentally increased concentrations of N and/or P in streams and rivers to examine how enrichment alters ecosystem structure (state: primary producer and consumer biomass and abundance) and function (rate: primary production, leaf breakdown rates, metabolism) at multiple trophic levels (primary producer, microbial heterotroph, primary and secondary consumers, and integrated ecosystem). Our synthesis included 184 studies, 885 experiments, and 3497 biotic responses to nutrient enrichment. We documented widespread increases in organismal biomass and abundance (mean response = +48%) and rates of ecosystem processes (+54%) to enrichment across multiple trophic levels, with no large differences in responses among trophic levels or between autotrophic or heterotrophic food-web pathways. Responses to nutrient enrichment varied with the nutrient added (N, P, or both) depending on rate versus state variable and experiment type, and were greater in flume and whole-stream experiments than in experiments using nutrient-diffusing substrata. Generally, nutrient-enrichment effects also increased with water temperature and light, and decreased under elevated ambient concentrations of inorganic N and/or P. Overall, increased concentrations of N and/or P altered multiple food-web pathways and trophic levels in lotic ecosystems. Our results indicate that preservation or restoration of biodiversity and ecosystem functions of streams and rivers requires management of nutrient inputs and consideration of multiple trophic pathways.

Resolving a paradox-high mercury deposition, but low bioaccumulation in northeastern Puerto Rico.

At a "clean air" trade winds site in northeastern Puerto Rico, we found an apparent paradox: atmospheric total mercury (THg) deposition was highest of any site in the USA Mercury Deposition Network, but assimilation into the local food web was quite low. Avian blood THg concentrations (n = 31, from eight species in five foraging guilds) ranged widely from 0.2 to 32 ng g-1 (median of 4.3 ng g-1). Within this population, THg was significantly greater at a low-elevation site near a wetland compared to an upland montane site, even when the comparison was limited to a single species. Overall, however, THg concentrations were approximately an order of magnitude lower than comparable populations in the continental U.S. In surface soil and sediment, potential rates of demethylation were 3 to 9-fold greater than those for Hg(II)-methylation (based on six radiotracer amendment incubations), but rates of change of ambient MeHg pools showed a slight net positive Hg(II)-methylation. Thus, the resolution of the paradox is that MeHg degradation approximately keeps pace with MeHg production in this landscape. Further, any net production of MeHg is subject to frequent flushing by high rainfall on chronically wet soils. The interplay of these microbial processes and hydrology appears to shield the local food web from adverse effects of high atmospheric mercury loading. This scenario may play out in other humid tropical ecosystems as well, but it is difficult to evaluate because coordinated studies of Hg deposition, methylation, and trophic uptake have not been conducted at other tropical sites.

Trace metals in Northern New England streams: Evaluating the role of road salt across broad spatial scales with synoptic snapshots.

Mobilization of trace metals from soils to surface waters can impact both human and ecosystem health. This study resamples a water sample archive to explore the spatial pattern of streamwater total concentrations of arsenic, cadmium, copper, lead, and zinc and their associations with biogeochemical controls in northern New England. Road deicing appears to result in elevated trace metal concentrations, as trace metal concentrations are strongly related to sodium concentrations and are most elevated when the sodium: chloride ratio is near 1.0 (~halite). Our results are consistent with previous laboratory and field studies that indicate cation exchange as a metal mobilization mechanism when road salt is applied to soils containing metals. This study also documents associations among sodium, chloride, dissolved organic carbon, iron, and metal concentrations, suggesting cation exchange mechanisms related to road deicing are not the only mechanisms that increase trace metal concentrations in surface waters. In addition to cation exchange, this study considers dissolved organic carbon complexation and oxidation-reduction conditions affecting metal mobility from soils in a salt-rich environment. These observations demonstrate that road deicing has the potential to increase streamwater trace metal concentrations across broad spatial scales and increase risks to human and ecosystem health.

Nutrient export and elemental stoichiometry in an urban tropical river.

Nutrient inputs to surface waters are particularly varied in urban areas, due to multiple nutrient sources and complex hydrologic pathways. Because of their close proximity to coastal waters, nutrient delivery from many urban areas can have profound impacts on coastal ecology. Relatively little is known about the temporal and spatial variability in stoichiometry of inorganic nutrients such as dissolved silica, nitrogen, and phosphorus (Si, N, and P) and dissolved organic matter in tropical urban environments. We examined nutrient stoichiometry of both inorganic nutrients and organic matter in an urban watershed in Puerto Rico served by municipal sanitary sewers and compared it to two nearby forested catchments using samples collected weekly from each river for 6 yr. Urbanization caused large increases in the concentration and flux of nitrogen and phosphorus (2- to 50-fold), but surprisingly little change in N:P ratio. Concentrations of almost all major ions and dissolved silica were also significantly higher in the urban river than the wildland rivers. Yield of dissolved organic carbon (DOC) was not increased dramatically by urbanization, but the composition of dissolved organic matter shifted toward N-rich material, with a larger increase in dissolved organic nitrogen (DON) than DOC. The molar ratio of DOC:DON was about 40 in rivers draining forested catchments but was only 10 in the urban river. Inclusion of Si in the assessment of urbanization's impacts reveals a large shift in the stoichiometry (Si:N and Si:P) of nutrient inputs. Because both Si concentrations and watershed exports are high in streams and rivers from many humid tropical catchments with siliceous bedrock, even the large increases in N and P exported from urban catchments result in delivery of Si, N, and P to coastal waters in stoichiometric ratios that are well in excess of the Si requirements of marine diatoms. Our data suggest that dissolved Si, often neglected in watershed biogeochemistry, should be included in studies of urban as well as less developed watersheds due to its potential significance for marine and lacustrine productivity.

Give and Take: A Watershed Acid Rain Mitigation Experiment Increases Baseflow Nitrogen Retention but Increases Stormflow Nitrogen Export.

In many temperate forested watersheds, hydrologic nitrogen export has declined substantially in recent decades, and many of these watersheds show enduring effects from historic acid deposition. A watershed acid remediation experiment in New Hampshire reversed many of these legacy effects of acid deposition and also increased watershed nitrogen export, suggesting that these two phenomena may be coupled. Here we examine stream nitrate dynamics in this watershed acid remediation experiment for indicators of nitrogen saturation in the terrestrial and aquatic ecosystems. Post-treatment, the (positive) slope of the relationship between nitrate concentration and discharge increased by a median of 82% ( p = 0.004). This resulted in greater flushing of nitrate during storm events, a key indicator of early stage nitrogen saturation. Hysteretic behavior of the concentration-discharge relationship indicated that the mobilization of soil nitrate pools was responsible for this increased flushing. In contrast to this evidence for nitrogen saturation in the terrestrial ecosystem, we found that nitrogen uptake increased, post-treatment, in the aquatic ecosystem, substantially attenuating growing-season nitrate flux by up to 71.1% ( p = 0.025). These results suggest that, as forests slowly recover from acid precipitation, terrestrial, and aquatic ecosystem nitrogen balance may be substantially altered.

Nitrogen removal rates in a frigid high-altitude river estimated by measuring dissolved N2 and N2O.

Rivers are important sites of both nitrogen removal and emission of nitrous oxide (N2O), a powerful greenhouse gas. Previous measurements have focused on nitrogen-rich temperate rivers, with cold, low-nitrogen river systems at high-altitude receiving less attention. Here, nitrogen removal rates were estimated by directly measuring dissolved N2 and N2O of the Yellow River in its source region of the Tibetan Plateau, a frigid high-altitude environment. We measured the dissolved N2 and N2O using N2:Ar ratio method and headspace equilibrium technique, respectively. Dissolved N2 in the river water ranged from 337 to 513 µmol N2 L-1, and dissolved N2O ranged from 10.4 to 15.4 nmol N2O L-1. Excess dissolved N2 (△N2) ranged from -8.6 to 10.5 µmol N2 L-1, while excess dissolved N2O (△N2O) ranged from 2.1 to 8.3 nmol N2O L-1; they were relatively low compared with most other rivers in the world. However, N2 removal fraction (△N2/DIN, average 21.6%) and EF5r values (N2O - N/NO3 - N, range 1.6 × 10-4-5.0 × 10-2) were comparable with many other rivers despite the high altitude for the Yellow River source region. Furthermore, the EF5r values increased with altitude. Estimated fluxes of N2 and N2O to the atmosphere from the river surface ranged from -67.5 to 93.5 mmol N m-2 d-1 and from 4.8 to 93.8 µmol N m-2 d-1, respectively, and the nitrogen removal from rivers was estimated to be 1.87 × 107 kg N yr-1 for the Yellow River source region. This is the first report of nitrogen removal for a frigid high-altitude river; the results suggest that N removal and N2O emission from cold high-altitude rivers should be considered in the global nitrogen budget.

The next generation of site-based long-term ecological monitoring: Linking essential biodiversity variables and ecosystem integrity.

Global change effects on biodiversity and human wellbeing call for improved long-term environmental data as a basis for science, policy and decision making, including increased interoperability, multifunctionality, and harmonization. Based on the example of two global initiatives, the International Long-Term Ecological Research (ILTER) network and the Group on Earth Observations Biodiversity Observation Network (GEO BON), we propose merging the frameworks behind these initiatives, namely ecosystem integrity and essential biodiversity variables, to serve as an improved guideline for future site-based long-term research and monitoring in terrestrial, freshwater and coastal ecosystems. We derive a list of specific recommendations of what and how to measure at a monitoring site and call for an integration of sites into co-located site networks across individual monitoring initiatives, and centered on ecosystems. This facilitates the generation of linked comprehensive ecosystem monitoring data, supports synergies in the use of costly infrastructures, fosters cross-initiative research and provides a template for collaboration beyond the ILTER and GEO BON communities.

LAGOS-NE: a multi-scaled geospatial and temporal database of lake ecological context and water quality for thousands of US lakes.

Understanding the factors that affect water quality and the ecological services provided by freshwater ecosystems is an urgent global environmental issue. Predicting how water quality will respond to global changes not only requires water quality data, but also information about the ecological context of individual water bodies across broad spatial extents. Because lake water quality is usually sampled in limited geographic regions, often for limited time periods, assessing the environmental controls of water quality requires compilation of many data sets across broad regions and across time into an integrated database. LAGOS-NE accomplishes this goal for lakes in the northeastern-most 17 US states.LAGOS-NE contains data for 51 101 lakes and reservoirs larger than 4 ha in 17 lake-rich US states. The database includes 3 data modules for: lake location and physical characteristics for all lakes; ecological context (i.e., the land use, geologic, climatic, and hydrologic setting of lakes) for all lakes; and in situ measurements of lake water quality for a subset of the lakes from the past 3 decades for approximately 2600-12 000 lakes depending on the variable. The database contains approximately 150 000 measures of total phosphorus, 200 000 measures of chlorophyll, and 900 000 measures of Secchi depth. The water quality data were compiled from 87 lake water quality data sets from federal, state, tribal, and non-profit agencies, university researchers, and citizen scientists. This database is one of the largest and most comprehensive databases of its type because it includes both in situ measurements and ecological context data. Because ecological context can be used to study a variety of other questions about lakes, streams, and wetlands, this database can also be used as the foundation for other studies of freshwaters at broad spatial and ecological scales.

Drivers of nitrogen transfer in stream food webs across continents.

Studies of trophic-level material and energy transfers are central to ecology. The use of isotopic tracers has now made it possible to measure trophic transfer efficiencies of important nutrients and to better understand how these materials move through food webs. We analyzed data from thirteen 15 N-ammonium tracer addition experiments to quantify N transfer from basal resources to animals in headwater streams with varying physical, chemical, and biological features. N transfer efficiencies from primary uptake compartments (PUCs; heterotrophic microorganisms and primary producers) to primary consumers was lower (mean 11.5%, range <1% to 43%) than N transfer efficiencies from primary consumers to predators (mean 80%, range 5% to >100%). Total N transferred (as a rate) was greater in streams with open compared to closed canopies and overall N transfer efficiency generally followed a similar pattern, although was not statistically significant. We used principal component analysis to condense a suite of site characteristics into two environmental components. Total N uptake rates among trophic levels were best predicted by the component that was correlated with latitude, DIN:SRP, GPP:ER, and percent canopy cover. N transfer efficiency did not respond consistently to environmental variables. Our results suggest that canopy cover influences N movement through stream food webs because light availability and primary production facilitate N transfer to higher trophic levels.

A case study characterizing animal fecal sources in surface water using a mitochondrial DNA marker.

Water quality impairment by fecal waste in coastal watersheds is a public health issue. The present study provided evidence for the use of a mitochondrial (mtDNA) marker to detect animal fecal sources in surface water. The accurate identification of fecal pollution is based on the notion that fecal microorganisms preferentially inhabit a host animal's gut environment. In contrast, mtDNA host-specific markers are inherent to eukaryotic host cells, which offers the advantage by detecting DNA from the host rather than its fecal bacteria. The present study focused on sampling water presumably from non-point sources (NPS), which can increase bacterial and nitrogen concentrations to receiving water bodies. Stream sampling sites located within the Piscataqua River Watershed (PRW), New Hampshire, USA, were sampled from a range of sites that experienced nitrogen inputs such as sewer and septic systems and suburban runoff. Three mitochondrial (mtDNA) gene marker assays (human, bovine, and canine) were tested from surface water. Nineteen sites were sampled during an 18-month period. Analyses of the combined single and multiplex assay results showed that the proportion of occurrence was highest for bovine (15.6%; n = 77) compared to canine (5.6%; n = 70) and human (5.7%; n = 107) mtDNA gene markers. For the human mtDNA marker, there was a statistically significant relationship between presence vs. absence and land use (Fisher's test p = 0.0031). This result was evident particularly for rural suburban septic, which showed the highest proportion of presence (19.2%) compared to the urban sewered (3.3%), suburban sewered (0%), and agricultural (0%) as well as forested septic (0%) sites. Although further testing across varied land use is needed, our study provides evidence for using the mtDNA marker in large watersheds.