Riverine Export Research Articles

Synthesis of nutrient and sediment export patterns in the Chesapeake Bay watershed: Complex and non-stationary concentration-discharge relationships

Derived from river monitoring data, concentration-discharge (C-Q) relationships are useful indicators of riverine export dynamics. A top-down synthesis of C-Q patterns was conducted for suspended sediment (SS), total phosphorus (TP), and total nitrogen (TN) for nine major tributaries (15 monitoring sites) to Chesapeake Bay, which represent diverse characteristics in terms of land use, physiography, and hydrological settings. Model coefficients from the recently-developed Weighted Regressions on Time, Discharge, and Season (WRTDS) method were used to make informative interpretation of C-Q relationships. Unlike many previous C-Q studies that focused on stormflow conditions, this approach allows simultaneous examination of various discharge conditions within an uncertainty framework. This synthesis on WRTDS coefficients (i.e., the sensitivity of concentration to discharge) has offered new insights on the complexity of watershed function. Results show that watershed export has been dominated by mobilization patterns for SS and TP (particulate-dominated species) and chemostasis patterns for TN (dissolved-dominated species) under many river discharge conditions. Among nine possible modalities of low-flow vs. high-flow patterns, the three most frequent modalities are mobilization vs. mobilization (17 cases), chemostasis vs. mobilization (13 cases), and chemostasis vs. chemostasis (7 cases), representing 82% of all 45 watershed-constituent pairs. The general lack of dilution patterns may suggest that none of these constituents has been supply-limited in these watersheds. For many watershed-constituent combinations, results show clear temporal non-stationarity in C-Q relationships under selected time-invariant discharges, reflecting major changes in dominant watershed sources due to anthropogenic actions. These results highlight the potential pitfalls of assuming fixed C-Q relationships in the record. Overall, this work demonstrates the utility of WRTDS model coefficients for interpretation of river water-quality data and for generation of sensible hypotheses on dominant processes in different watersheds. The approach is readily adaptable to other river systems, where long-term discretely-sampled data are available, to decipher complex interactions between hydrological and biogeochemical processes.

Do Regional Aerosols Contribute to the Riverine Export of Dissolved Black Carbon?

AbstractThe fate of black carbon (BC), a stable form of thermally altered organic carbon produced during biomass and fuel combustion, remains an area of uncertainty in the global carbon cycle. The transfer of photosynthetically derived BC into extremely long‐term oceanic storage is of particular significance and rivers are the key linkage between terrestrial sources and oceanic stores. Significant fluvial fluxes of dissolved BC to oceans result from the slow release of BC from degrading charcoal stocks; however, these fluvial fluxes may also include undetermined contributions of aerosol BC, produced by biomass and fossil fuel combustion, which are deposited in river catchments following atmospheric transport. By investigation of the Paraíba do Sul River catchment in Southeast Brazil we show that aerosol deposits can be substantial contributors to fluvial fluxes of BC. We derived spatial distributions of BC stocks within the catchment associated with soil charcoal and with aerosol from both open biomass burning and fuel combustion. We then modeled the fluvial concentrations of dissolved BC (DBC) in scenarios with varying rates of export from each stock. We analyzed the ability of each scenario to reproduce the variability in DBC concentrations measured in four data sets of river water samples collected between 2010 and 2014 and found that the best performing scenarios included a 5–18% (135–486 Mg DBC year−1) aerosol contribution. Our results suggest that aerosol deposits of BC in river catchments have a shorter residence time in catchments than charcoal BC and, therefore, contribute disproportionately (with respect to stock magnitude) toward fluvial fluxes of BC.

Antropogenic input of nitrogen and riverine export from a Mediterranean catchment. The Celone, a temporary river case study

In a catchmnent, nitrogen (N) export from terrestrial ecosystems to rivers is controlled by the hydrological processes and N balance. In the present paper, the soil system N budget and riverine export were quantified in a Mediterranean watershed, the Celone (South-East, Italy). The study area (72km2) consists of agricultural land (70%) and deciduous forests with three small residential areas. Major N inputs derived from fertilizers and animal manure, corresponding for the whole watershed area to 68 and 12kgNha−1yr−1, respectively. N input from point sources was 1% of total input and atmospheric depositions measured in a gauging station near the study area was quantified in ∼6kgNha−1yr−1. Crop N uptake was the main N output from agricultural land; it was estimated in ∼37kgNha−1yr−1 by using data on crop yields provided by local farmers. Total flux of N in surface water was quantified for a year at the outlet of the study area by using continuous measures of flow and discrete measures of N concentrations carried out with a different frequency during the normal and low flow and when floods occurred. The hydrological regime of the stream, which is a temporary river, plays an important role in N transport. Water quality was found to vary considerably through the year in terms of both nutrient concentrations and loads. Riverine N export was quantified in 41% of total N input, and it was mainly transported during flood events, ∼60% of the annual load was delivered during floods occurred in 38 days. Organic nitrogen and nitrate were the main N forms in surface water, and the contribution per hectare was about 24 and 14kgN, respectively. On a yearly basis, the difference between N inputs and outputs including riverine export was estimated in about 4kgNha−1yr−1 for the whole watershed area. This amount partly accumulates in soils in different N forms and the remaining part, mainly in form of nitrate, percolates through unsaturated soil towards groundwater. This study reports an important analysis of N pollution in a Mediterranean watershed with a temporary river system and limited data availability. Data acquisition and handling have proved to be an important challenge to overcome in N balance quantification. The results and the methodology of the present work can be useful for understanding nitrogen loss dynamics and for functional water management and land use planning.

Snowmelt controls on concentration‐discharge relationships and the balance of oxidative and acid‐base weathering fluxes in an alpine catchment, East River, Colorado

AbstractAlthough important for riverine solute and nutrient fluxes, the connections between biogeochemical processes and subsurface hydrology remain poorly characterized. We investigate these couplings in the East River, CO, a high‐elevation shale‐dominated catchment in the Rocky Mountains, using concentration‐discharge (C‐Q) relationships for major cations, anions, and organic carbon. Dissolved organic carbon (DOC) displays a positive C‐Q relationship with clockwise hysteresis, indicating mobilization and depletion of DOC in the upper soil horizons and emphasizing the importance of shallow flow paths during snowmelt. Cation and anion concentrations demonstrate that carbonate weathering, which dominates solute fluxes, is promoted by both sulfuric acid derived from pyrite oxidation in the shale bedrock and carbonic acid derived from subsurface respiration. Sulfuric acid weathering dominates during base flow conditions when waters infiltrate below the inferred pyrite oxidation front, whereas carbonic acid weathering plays a dominant role during snowmelt as a result of shallow flow paths. Differential C‐Q relationships between solutes suggest that infiltrating waters approach calcite saturation before reaching the pyrite oxidation front, after which sulfuric acid reduces carbonate alkalinity. This reduction in alkalinity results in CO2 outgassing when waters equilibrate to surface conditions, and reduces the riverine export of carbon and alkalinity by roughly 33% annually. Future changes in snowmelt dynamics that control the balance of carbonic and sulfuric acid weathering may substantially alter carbon cycling in the East River. Ultimately, we demonstrate that differential C‐Q relationships between major solutes can provide unique insights into the complex subsurface flow and biogeochemical dynamics that operate at catchment scales.

Nitrate concentrations and source identification in a Mediterranean river system

As most of the rural areas in Mediterranean Region, the Celone river basin (Apulia, Italy) is characterized by diffuse pollution from agriculture. Nutrient export from terrestrial ecosystem to the river is controlled by hydrological processes and nutrient input from point and non-point sources. In the present paper, anthropogenic nitrogen input and the riverine export of nitrate-N (N-NO3) were quantified through several survey campaigns, continuous measures of streamflow and discrete determinations of concentrations. A different frequency for water sampling was used, respectively, during floods and during normal- or low-flow state. Nitrate-N concentrations in surface waters were found to vary considerably across the seasons. The riverine N-NO3 export was estimated to be about 14 kg N ha−1 year−1. Major N inputs derived from fertilizers and animal manure correspond to 68.2 and 24 kg N ha−1 year−1, respectively. Nitrate-N fluxes in stream during flood events (21 days) accounted for about 63% of the annual loading. This study clearly demonstrates the importance of flood event contributions to the annual nutrient loading delivered to the reservoir and that more efficient fertilizers and management practices are needed to reduce N input.

Europe-wide estuarine export and surface water concentrations of PFOS and PFOA

The STREAM-EU model was used to predict the water concentrations, estuarine export and retention of perfluorooctane sulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) in the eleven most populated European river catchments to provide a European-wide perspective on the contamination by these substances.Emissions of PFOS and PFOA to those catchments were calculated based on population, wealth and wastewater treatment plant (WWTP) coverage and efficiency using a previously published method and used as model input. Our estimated emissions showed the lowest values for the Thames catchment (PFOS: 0.4 ton/y; PFOA: 0.2 ton/y) and the highest values for the Rhine for PFOS (1.6 ton/y) and for the Dnieper for PFOA (1.7 ton/y).The model predicted concentrations agreed reasonable well with the existing range of measurements, apart from for PFOA in the River Po, where there is a known historical industrial contamination, and PFOS in the Rhone River, where results were much higher than the few measurements available. It was concerning that the model predicted that the surface water EQS for PFOS (0.65 ng/L) was exceeded by a wide margin in all the eleven studied European river catchments.The total calculated riverine export to the seas from the eleven catchments was 4.5 ton/y of PFOS and 3.7 ton/y of PFOA with highest exported quantities from the Rhine (PFOS: 1.0 ton/y; PFOA: 1.0 ton/y) and Danube estuaries (PFOS: 0.9 ton/y; PFOA: 0.7 ton/y). For the seas where the rivers discharge, riverine discharge of PFOS was estimated to be 2.5–30 times more important as an input than atmospheric deposition, whereas for PFOA the opposite was true (atmospheric deposition was 2–10 times more important) except for very small seas.

Changes in anthropogenic nitrogen and phosphorus inputs to the St. Lawrence sub‐basin over 110 years and impacts on riverine export

AbstractHuman activities have increased the flow of nitrogen (N) and phosphorus (P) over much of the Earth, leading to increased agricultural production, but also the degradation of air, soil, and water quality. Here we quantify the sources of anthropogenic N and P inputs to 76 watersheds of the St. Lawrence Basin (SLB) throughout the 20th century using NANI/NAPI (net anthropogenic N/P input to watersheds), a mass balance modeling approach, and estimate the fraction of these inputs exported to adjacent rivers. Our results show that since 1901, NANI and NAPI increased 4.5‐fold and 3.8‐fold, respectively, with a peak in 1991 mainly due to high atmospheric N deposition and P fertilizer application. However, the relative increase over the course of the last century was much higher in certain watersheds, particularly those where there was greater urbanization. Ranges in NANI and NAPI vary greatly among watersheds (110 to 9351 kg N km−2 yr−1 and 0.16 to 1938 kg P km−2 yr−1, respectively in 2011) and are strongly related to riverine fluxes (R2 = 0.87 and 0.71 for N and P, respectively). Our results suggest that 22% of NANI (ranging from 11% to 68% across watersheds) and 17% of NAPI (ranging from 3% to 173%) are exported to rivers. Predominant sources of inputs vary spatially and through time largely due to changes in farming practices. By tracking the main sources of inputs to specific watersheds and through time, our work provides insights for N and P management. Reduction strategies will likely need to be watershed specific, although through time, our results clearly show the large‐scale impact of targeted legislation.

Earthquake-triggered increase in biospheric carbon export from a mountain belt

On geological time scales, the erosion of carbon from the terrestrial biosphere and its burial in sediments can counter CO2 emissions from the solid Earth. Earthquakes may increase the erosion of this biospheric carbon and supply it to mountain rivers by triggering landslides, which rapidly strip hillslopes of vegetation and soil. Over the long term, elevated river sediment loads may promote more efficient carbon burial. However, riverine export of earthquake-mobilized carbon has remained poorly constrained. Here we quantify biospheric carbon discharge by the Zagunao River following a large earthquake with a unique set of samples collected before and after the A.D. 2008 M w 7.9 Wenchuan (China) earthquake. Radioactive and stable carbon isotopes are used to isolate the biospheric carbon, accounting for rock-derived organic carbon inputs. River discharge of biospheric carbon doubled in the downstream reaches, characterized by moderate landslide impact, following the earthquake. The rapid export of carbon from earthquake-triggered landslides appears to outpace its degradation on hillslopes while sediment loads are elevated. This means that enhanced river discharge of biospheric carbon following large earthquakes can link active tectonics to CO2 drawdown.

Multi-decadal increases in dissolved organic carbon and alkalinity flux from the Mackenzie drainage basin to the Arctic Ocean

Riverine exports of organic and inorganic carbon (OC, IC) to oceans are intricately linked to processes occurring on land. Across high latitudes, thawing permafrost, alteration of hydrologic flow paths, and changes in vegetation may all affect this flux, with subsequent implications for regional and global carbon (C) budgets. Using a unique, multi-decadal dataset of continuous discharge coupled with water chemistry measurements for the Mackenzie River, we show major increases in dissolved OC (DOC) and IC (as alkalinity) fluxes since the early 1970s, for a watershed that covers 1.8 M km2 of northwestern Canada, and provides substantial inputs of freshwater and biogeochemical constituents to the Arctic Ocean. Over a 39-year period of record, DOC flux at the Mackenzie mouth increased by 39.3% (44.5 ± 22.6 Gmol), while alkalinity flux increased by 12.5% (61.5 ± 60.1 Gmol). Isotopic analyses and substantial increases in sulfate flux indicate that increases in alkalinity are driven by accelerating sulfide oxidation, a process that liberates IC from rock and soils in the absence of CO2 consumption. Seasonal and sub-catchment trends suggest that permafrost thaw plays an important role in the observed increases in DOC and alkalinity: sub-catchment increases for all constituents are confined to northern, permafrost-affected regions, while observed increases in autumn to winter are consistent with documented landscape-scale changes that have resulted from changing thaw dynamics. This increase in DOC and sulfide-derived alkalinity represents a substantial intensification of land-to-ocean C mobilization, at a level that is significant within the regional C budget. The change we observe, for example, is similar to current and projected future rates of CO2 consumption by weathering in the Mackenzie basin.

Crop rotation and Raccoon River nitrate

Previous research predicted that the biofuel-driven expansion of corn (Zea mays L.) area would increase riverine export of nitrate-nitrogen (NO3-N) in the Mississippi River Basin. Accurate information about water quality trends in agricultural watersheds is needed to better inform agricultural policy and to help quantify the effectiveness of field and landscape management practices. This study was designed to (1) characterize possible NO3-N trends within the Raccoon River watershed (RRWS) of central Iowa, and (2) explore links between the relative areas planted to corn and soybean (Glycine max) and water quality. We examined NO3-N concentration and loading data from more than 60 main stem river and tributary locations in the RRWS for the period 1999 to 2014. In addition, we assessed the role of climate, crop rotations, and simplified annual N budgets on NO3-N concentrations and loads to show that expansion of corn area has not increased Raccoon River NO3-N levels. Nitrate-N concentrations have not increased as corn area increased 19% and fertilizer N inputs increased 24% since 1999. We conclude that expansion of corn area at the expense of soybean may be affecting water quality. Better management of soybean in a corn–soybean rotation should reduce NO3-N export from the watershed, and reducing throughput of water in this artificially drained system will improve water quality.

Modeling forest/agricultural and residential nitrogen budgets and riverine export dynamics in catchments with contrasting anthropogenic impacts in eastern China between 1980–2010

This study quantified the long-term response of riverine total nitrogen (TN) export to changes in net anthropogenic nitrogen inputs to forest/agricultural (NANIFA) and residential (NANIR) systems across three catchments affected by low (LD), medium (MD), and high (HD) human impacts in eastern China. Annual NANIFA increased by 63–87% in 1980–1999, followed by 0% (LD), −23% (MD) and −40% (HD) changes of NANIFA in 2000–2010, resulting in a net increase of 56–78% in NANIFA in 1980–2010. Annual NANIR increased by 101–152% in the three catchments in 1980–2010. Land-use showed a 58–65% increase in developed land area (D%) and a 96–108% increase in agricultural lands with improved drainage systems (AD%) over the study period. In response to changes in NANIFA, NANIR and land-use, riverine TN flux continuously increased 3.0- to 6.1-fold in the three catchments over the past 31 years. For each catchment, an empirical model incorporating annual NANIFA, NANIR, water discharge, D%, and AD% was developed (R2=0.93−0.97) for predicting and quantifying sources of annual riverine TN fluxes. The model estimated that NANIFA, NANIR and other N sources (e.g., natural background, legacy, and industrial N sources) contributed 27–90%, 0–45%, and 10–28% of riverine TN fluxes, respectively. Model results were consistent with spatio-temporal changes of riverine chloride, ammonium, nitrate, dissolve oxygen and pH, as well as changes in available N levels in agricultural soils. In terms of N source management, reduction of NANIFA in catchment LD and NANIR in catchment HD would have the greatest impact on reducing riverine TN fluxes. Furthermore, changes in land use and climate as well as legacy N should be considered in developing N pollution control strategies.

Chromium fluxes and speciation in ultramafic catchments and global rivers

The role of ultramafic rock weathering in global chromium (Cr) budgets and how these may be linked to changes in Earth's climate over geological time was investigated by compiling dissolved Cr speciation and fluxes in (i) 73 of the worlds large rivers, representing ~35% of global river discharge to oceans, and in (ii) an ultramafic catchment (Putah Creek watershed) in the California Coast Range Mountains, USA. Weathering of ultramafic rocks creates ideal conditions for Cr release and redox cycling. Alkaline river water draining ultramafic catchments is naturally enriched in Cr (up to 582nmol/L). Chromium concentrations increase with the extent of rock–water interaction (as indicated by correlations with Mg and HCO3− concentrations and pH), whereas Cr cycling depends on the availability of electron donors and acceptors (e.g. Fe(II), organics, Mn(III/IV)-oxides). Thus, dissolved Cr is exported from ultramafic watersheds as both the soluble hexavalent Cr(VI) species (MgCrO4, CrO42−), and also as trivalent Cr(III) species (CrOH(CO3)22−), Cr(OH)3). The latter Cr(III) species were previously thought to have low solubility. Ultramafic catchments have higher area normalized Cr and major ion fluxes for a given runoff when compared to global rivers and may thus have a disproportional impact on global Cr-budgets, long-term carbon cycling and moderation of the Earth's climate over million-year timescales. Riverine export fluxes of Cr are linearly correlated with fluxes of Ca+Mg and alkalinity in both ultramafic catchments and global rivers. This suggests that silicate weathering is a key control on Cr fluxes and that major element weathering fluxes may be a useful proxy for estimating Cr and other trace element fluxes. Globally, present-day riverine dissolved Cr fluxes to oceans are spatially variable and estimated to be 1.7×109mol/yr, three times higher than previously reported. However, when geochemical reactions in estuaries are considered riverine Cr fluxes may be lower. Throughout Earth history, Cr weathering fluxes also may have varied with changes in the global distribution of ultramafic rocks, Cr concentrations in continental rocks, or climactic conditions (e.g. atmospheric CO2 levels or runoff). Such variation should be considered when interpreting both modern and past Cr seawater isotopic compositions and residence times.

Isotopic evidence on multiple sources of nitrogen in the northern Jiulong River, Southeast China

Riverine export accounts for a large portion of estuarine and coastal nutrients and could lead to severe eutrophication. However, nitrogen (N) sources at the catchment scale remain unclear because of spatial and temporal variations. The stable isotope of 15N, which has been proven to be effective in deducing sources and testing biogeochemical behaviours, is applied in this study to explore multiple sources of riverine N and their nutrient concentrations in the northern section of Jiulong River catchment seasonally. Results show that drastic seasonal variation in external nitrate sources occurs in the river; manure and sewage dominate during base flows, whereas organic N in soil and atmospheric deposition dominate during storm flows. The external sources change throughout the year depending on the environmental conditions. Furthermore, riverine nitrate import in the northern Jiulong River is dominated by an external process because of the increasing NO3−-δN15 value accompanied by an increasing NO3−-N concentration (p < 0.05). Manure and sewage show an increasing proportion of contribution to the downstream section and proximity to the estuary with high NO3−-δN15 values. The findings in this study highlight the fact that the isotopic composition of riverine NO3−-N provides effective evidence for identifying nitrate sources in catchments. The study could serve as a critical reference for nutrient management and eutrophication reversal.

Global carbon export from the terrestrial biosphere controlled by erosion

Riverine export of particulate organic carbon (POC) to the ocean affects the atmospheric carbon inventory over a broad range of timescales. On geological timescales, the balance between sequestration of POC from the terrestrial biosphere and oxidation of rock-derived (petrogenic) organic carbon sets the magnitude of the atmospheric carbon and oxygen reservoirs. Over shorter timescales, variations in the rate of exchange between carbon reservoirs, such as soils and marine sediments, also modulate atmospheric carbon dioxide levels. The respective fluxes of biospheric and petrogenic organic carbon are poorly constrained, however, and mechanisms controlling POC export have remained elusive, limiting our ability to predict POC fluxes quantitatively as a result of climatic or tectonic changes. Here we estimate biospheric and petrogenic POC fluxes for a suite of river systems representative of the natural variability in catchment properties. We show that export yields of both biospheric and petrogenic POC are positively related to the yield of suspended sediment, revealing that POC export is mostly controlled by physical erosion. Using a global compilation of gauged suspended sediment flux, we derive separate estimates of global biospheric and petrogenic POC fluxes of 157(+74)(-50) and 43(+61)(-25) megatonnes of carbon per year, respectively. We find that biospheric POC export is primarily controlled by the capacity of rivers to mobilize and transport POC, and is largely insensitive to the magnitude of terrestrial primary production. Globally, physical erosion rates affect the rate of biospheric POC burial in marine sediments more strongly than carbon sequestration through silicate weathering. We conclude that burial of biospheric POC in marine sediments becomes the dominant long-term atmospheric carbon dioxide sink under enhanced physical erosion.

Neogene cratonic erosion fluxes and landform evolution processes from regional regolith mapping (Burkina Faso, West Africa)

The regionally correlated and dated regolith–paleolandform sequence of Sub-Saharan West Africa offers a unique opportunity to constrain continental-scale regolith dynamics as the key part of the sediment routing system. In this study, a regolith mapping protocol is developed and applied at the scale of Southwestern Burkina Faso. Mapping combines field survey and remote sensing data to reconstruct the topography of the last pediplain that formed over West Africa in the Early and Mid-Miocene (24–11Ma). The nature and preservation pattern of the pediplain are controlled by the spatial variation of bedrock lithology and are partitioned among large drainage basins. Quantification of pediplain dissection and drainage growth allows definition of a cratonic background denudation rate of 2m/My and a minimum characteristic timescale of 20Ma for shield resurfacing. These results may be used to simulate minimum export fluxes of drainage basins of constrained size over geological timescales. Background cratonic denudation results in a clastic export flux of ~4t/km2/year, which is limited by low denudation efficiency of slope processes and correlatively high regolith storage capacity of tropical shields. These salient characteristics of shields' surface dynamics would tend to smooth the riverine export fluxes of shields through geological time.

Effects of a tectonically-triggered wave of incision on riverine exports and soil mineralogy in the Luquillo Mountains of Puerto Rico

We document the long-term response of a tropical mountain stream to tectonics and show how this response exerts a first-order spatial control on the delivery of weatherable minerals to forest soils in its catchment. These minerals, in turn, affect nutrient availability and biogeochemical cycles. The Luquillo Mountains is a rainforest-covered isolated massif on Puerto Rico Island. Instead of displaying typical concave-up long profiles of most mountain streams, the rivers draining the southern flank of these mountains display a systematic downstream increase in gradient. Using concentrations of in situ-produced cosmogenic 10Be in quartz from river sediments, we find that the downstream steepening is due to the propagation of an erosion wave along the rivers from the coastal plain upstream, toward the headwaters. Decreased 10Be concentration along the steepened reaches (knickpoint faces or knickzones) results from faster denudation of the valley sides around the steepened reaches. The upper portions of the steepened reaches (knickpoint lips) cluster in elevation around the altitude of a flat-lying, dissected surface traceable around the Luquillo Mountains. Inspection of similar platforms over other parts of Puerto Rico and in the surrounding seas suggests that the dissected surface is an uplifted shore platform formed in the Early Pliocene (∼4 My ago). Upstream of the knickpoints, rivers possess alluvial reaches still graded to this uplifted platform. They represent relict profiles of the rivers that once drained an island that was surrounded by the now-uplifted platform. We conclude that the knickpoints initiated when the platform started to rise above the Caribbean Sea. The knickpoints then started propagating upstream, defining the front of a slow-moving wave of erosion. 10Be-derived catchment-scale denudation rates measured above and below the knickpoints indicate a 30–210% increase in denudation associated to the passage of the knickpoints. Mineralogical analysis of the soils show that soils upstream of the knickpoints overlie deeply depleted saprolite and only contain recalcitrant minerals, whereas soils downstream of the knickpoints experience an input of fresh, weatherable minerals. This influx considerably increases nutrient availability (Porder et al., 2015) in a forest where atmospheric inputs are otherwise the main source of nutrients (McDowell et al., 1990; Zarin and Johnson, 1995; Pett-Ridge, 2009).

A lagged variable model for characterizing temporally dynamic export of legacy anthropogenic nitrogen from watersheds to rivers.

Legacy nitrogen (N) sources originating from anthropogenic N inputs (NANI) may be a major cause of increasing riverine N exports in many regions, despite a significant decline in NANI. However, little quantitative knowledge exists concerning the lag effect of NANI on riverine N export. As a result, the N leaching lag effect is not well represented in most current watershed models. This study developed a lagged variable model (LVM) to address temporally dynamic export of watershed NANI to rivers. Employing a Koyck transformation approach used in economic analyses, the LVM expresses the indefinite number of lag terms from previous years' NANI with a lag term that incorporates the previous year's riverine N flux, enabling us to inversely calibrate model parameters from measurable variables using Bayesian statistics. Applying the LVM to the upper Jiaojiang watershed in eastern China for 1980-2010 indicated that ~97% of riverine export of annual NANI occurred in the current year and succeeding 10 years (~11 years lag time) and ~72% of annual riverine N flux was derived from previous years' NANI. Existing NANI over the 1993-2010 period would have required a 22% reduction to attain the target TN level (1.0 mg N L(-1)), guiding watershed N source controls considering the lag effect. The LVM was developed with parsimony of model structure and parameters (only four parameters in this study); thus, it is easy to develop and apply in other watersheds. The LVM provides a simple and effective tool for quantifying the lag effect of anthropogenic N input on riverine export in support of efficient development and evaluation of watershed N control strategies.

In-stream sources and links between particulate and dissolved black carbon following a wildfire

The occurrence of wildfires is expected to increase with the progression of climate change. These natural burn events can drastically alter the geomorphology and hydrology of affected areas and are one of the primary sources of black carbon (BC) in the environment. BC can be mobilized from soils and charcoal in fire-affected watersheds, potentially impacting downstream water quality. In June of 2012, the High Park Fire burned a large portion of the Cache La Poudre River watershed located in the Colorado Rocky Mountains. Seasonal riverine export of BC in both the dissolved (DBC) and particulate (PBC) phase was compared between burned and unburned sections of the watershed during the year following the High Park Fire. There was little difference in overall DBC concentration between sites, however seasonal changes in DBC quality reflected a shift in hydrology and associated DBC source between peak and base flow conditions. PBC export was substantially larger in fire-affected areas of the watershed during periods of overland flow. Our findings suggest that export processes of BC in the particulate and dissolved phase are decoupled in burned watersheds and that, in addition to DBC, the export of PBC could be a significant contributor to the cycling of charcoal in freshwater ecosystems.

A carbon balance of Norway: terrestrial and aquatic carbon fluxes

Northern landscapes accumulate carbon in vegetation and soils while rivers transport significant amounts of land-derived carbon to coastal areas. Here, we quantify carbon sources and sinks in main ecosystems (forests, peatlands, mountains, agricultural areas, lakes) in Norway for 1990–2008, and compare riverine carbon transport with terrestrial carbon accumulation in Norway’s four major discharge areas. Mean annual carbon accumulation (6.0 ± 0.9 Tg C; 19 g C m−2) in terrestrial ecosystems balanced 40 % of national greenhouse gas emissions. The area-normalized terrestrial sink strength declined in the following order (in g C m−2 year−1): tree biomass (40 ± 3) > peatlands (19 ± 15) > forest soils (9 ± 1) ≫ lakes (2 ± 1) > mountains (0.5 ± 0.3), while agricultural soils were sources of carbon (−36 ± 74). The most precise estimate in the carbon balance was for tree biomass, because of the underlying forest inventory data. Poor data on land management and soil type in agricultural soils, and on (former) drainage and peatland type resulted in high uncertainty in carbon loss and uptake estimates in agricultural soils and peatlands, which impacted the uncertainty in total terrestrial carbon accumulation. Also, carbon losses from disturbance in organic soil types were poorly constrained. Riverine coastal inputs of land-derived organic carbon (OC) were 1.0 ± 0.1 Tg C year−1 (3.0 g C m−2 year−1), with highest area-specific riverine export in western (4.5 g C m−2 year−1) and lowest (1.7 g C m−2 year−1) in subarctic Norway. In west and middle Norway, river OC export was approximately equal to carbon accumulation in soils and wetlands, while it was 50 % of soil and wetland carbon accumulation in southeast and subarctic Norway. Lateral aquatic transport of carbon is not explicitly accounted for in forest soil carbon accumulation estimates, although aquatic fluxes represent a climate-dependent carbon loss from soil carbon pools. The lack of methods that adequately account for lateral fluxes in carbon balances adds considerable uncertainty to soil carbon sink estimates. Climate warming and associated changes in precipitation may result in substantial alterations of terrestrial and aquatic carbon fluxes, with uncertain implications for the terrestrial carbon sink of northern landscapes.

Influence of legacy phosphorus, land use, and climate change on anthropogenic phosphorus inputs and riverine export dynamics

A quantitative understanding of riverine phosphorus (P) export in response to changes in anthropogenic P inputs (NAPI), land use and climate is critical for developing effective watershed P control measures. This study indicated that annual riverine TP export for the six catchments of the Yongan River watershed in eastern China increased 4.1–30.3-fold over the 1980–2010 period. Increased riverine TP export resulted from a 61–85 % increase in NAPI and a 2.6–14.6-fold increase in riverine export fraction of NAPI due to 36–43, 30–125, and 65–76 % increases in developed land area (D%), drained agricultural land area (DA%), and storm events, respectively. For the 31-year cumulative record, 1.6–14 % of NAPI was exported by rivers, 40–64 % was stored in the upper 20 cm of agricultural soils, and 30–55 % was retained in other landscape positions. An empirical model that incorporates annual NAPI, precipitation, D%, and DA% accounted for 94 % of the variation in annual riverine TP fluxes across the six catchments and 31 years. The model estimated that NAPI and legacy P contributed 42–92 % and 8–58 % of annual riverine TP flux, respectively. The model forecasts an 8–18 % increase in riverine TP flux by 2030 due to a 4 % increase in precipitation with no changes in NAPI and land use compared to the 2000–2010 baseline condition. Enhanced export of NAPI and legacy P by changes in land use and climate will delay the decrease in riverine P flux in response to NAPI reductions and should be considered in developing and assessing watershed P management strategies.