Modelling Subarctic watershed dissolved organic carbon response to hydroclimatic regime.

Shifts in hydroclimatic regimes associated with global climate change may impact freshwater availability and quality. In high latitudes of the northern hemisphere, where vast quantities of carbon are stored terrestrially, explaining landscape-scale carbon (C) budgets and associated pollutant transfer is necessary for understanding the impact of changing hydroclimatic regimes. We used a dynamic modelling approach to simulate streamflow, DOC concentration, and DOC export in a northern Canadian catchment that has undergone notable climate warming, and will continue to for the remainder of this century. The Integrated Catchment model for Carbon (INCA-C) was successfully calibrated to a multi-year period (2012-2016) that represents a range in hydrologic conditions. The model was subsequently run over 30-year periods representing baseline and two future climate scenarios. Average discharge is predicted to decrease under an elevated temperature scenario (22-27 % of baseline) but increase (116-175 % of baseline) under an elevated temperature and precipitation scenario. In the latter scenario the nival hydroclimatic regime is expected to shift to a combined nival and pluvial regime. Average DOC flux over 30 years is predicted to decrease (24-27 % of baseline) under the elevated temperature scenario, as higher DOC concentrations are offset by lower runoff. Under the elevated temperature and precipitation scenario, results suggest an increase in carbon export of 64-81 % above baseline. These increases are attributed to greater connectivity of the catchment. The largest increase in DOC export is expected to occur in early winter. These predicted changes in DOC export, particularly under a climate that is warmer and wetter could be part of larger ecosystem change and warrant additional monitoring efforts in the region.

Small artificial waterbodies are widespread and persistent emitters of methane and carbon dioxide.

Inland waters play an active role in the global carbon cycle and emit large volumes of the greenhouse gases (GHGs), methane (CH4 ) and carbon dioxide (CO2 ). A considerable body of research has improved emissions estimates from lakes, reservoirs and rivers but recent attention has been drawn to the importance of small, artificial waterbodies as poorly quantified but potentially important emission hotspots. Of particular interest are emissions from drainage ditches and constructed ponds. These waterbody types are prevalent in many landscapes and their cumulative surface areas can be substantial. Furthermore, GHG emissions from constructed waterbodies are anthropogenic in origin and form part of national emissions reporting, whereas emissions from natural waterbodies do not (according to Intergovernmental Panel on Climate Change guidelines). Here, we present GHG data from two complementary studies covering a range of land uses. In the first, we measured emissions from nine ponds and seven ditches over a full year. Annual emissions varied considerably: 0.1-44.3 g CH4  m-2  year-1 and -36-4421 g CO2  m-2  year-1 . In the second, we measured GHG concentrations in 96 ponds and 64 ditches across seven countries, covering subtropical, temperate and sub-arctic biomes. When CH4 emissions were converted to CO2  equivalents, 93% of waterbodies were GHG sources. In both studies, GHGs were positively related to nutrient status (C, N, P), and pond GHG concentrations were highest in smallest waterbodies. Ditch and pond emissions were larger per unit area when compared to equivalent natural systems (streams, natural ponds). We show that GHG emissions from natural systems should not be used as proxies for those from artificial waterbodies, and that artificial waterbodies have the potential to make a substantial but largely unquantified contribution to emissions from the Agriculture, Forestry and Other Land Use sector, and the global carbon cycle.

A 25-year retrospective analysis of factors influencing success of aluminum treatment for lake restoration.

For more than 50 years, aluminum (Al)-salts have been used with varying degrees of success to inactivate excess mobile phosphorus (P) in lake sediments and restore lake water quality. Here, we analyzed the factors influencing effectiveness and longevity of Al-treatments performed in six Swedish lakes over the past 25 years. Trends in post-treatment measurements of total phosphorus (TP), Chlorophyll a (Chl_a), Secchi disk depth (SD) and internal P loading rates (Li) were analyzed and compared to pre-treatment conditions. All measured water quality parameters improved significantly during at least the first 4 years post-treatment and determination of direct effects of Al-treatment on sediment P release (Li) was possible for three lakes. Improvements in TP (-29 to -80%), Chl_a (-50 to -78%), SD (7 to 121%) and Li (-68 to -94%) were observed. Treatment longevity, determined via decreases in surface water TP after treatment, varied from 7 to >47 years. Lake type, Al dose, and relative watershed area were related to longevity. In addition, greater binding efficiency between Al and P was positively related to treatment longevity, which has not previously been shown. Our findings also demonstrate that adequate, long-term monitoring programs, including proper determination of external loads, are crucial to document the effect of Al-treatment on sediment P release and lake water quality.

Optimization of aluminum treatment efficiency to control internal phosphorus loading in eutrophic lakes.

Historical accumulation of phosphorus (P) in lake sediment often contributes to and sustains eutrophic conditions in lakes, even when external sources of P are reduced. The most cost-effective and commonly used method to restore the balance between P and P-binding metals in the sediment is aluminum (Al) treatment. The binding efficiency of Al, however, has varied greatly among treatments conducted over the past five decades, resulting in substantial differences in the amount of P bound per unit Al. We analyzed sediment from seven previously Al treated Swedish lakes to investigate factors controlling binding efficiency. In contrast to earlier work, lake morphology was negatively correlated to binding efficiency, meaning that binding efficiency was higher in lakes with steeply sloping bathymetry than in lakes with more gradually sloping bottoms. This was likely due to Al generally being added directly into the sediment, and not to the water column. Higher binding efficiencies were detected when Al was applied directly into the sediment, whereas the lowest binding efficiency was detected where Al was instead added to the water column. Al dose, mobile sediment P and lake morphology together explained 87% of the variation in binding efficiency among lakes where Al was added directly into the sediment. This led to the development of a model able to predict the optimal Al dose to maximize binding efficiency based on mobile sediment P mass and lake morphology. The predictive model can be used to evaluate cost-effectiveness and potential outcomes when planning Al-treatment using direct sediment application to restore water quality in eutrophic lakes.

Optimizing land management strategies for maximum improvements in lake dissolved oxygen concentrations.

Eutrophication and anoxia are unresolved issues in many large waterbodies. Globally, management success has been inconsistent, highlighting the need to identify approaches which reliably improve water quality. We used a process-based model chain to quantify effectiveness of terrestrial nutrient control measures on in-lake nitrogen, phosphorus, chlorophyll and dissolved oxygen (DO) concentrations in Lake Simcoe, Canada. Across a baseline period of 2010-2016 hydrochemical outputs from catchment models INCA-N and INCA-P were used to drive the lake model PROTECH, which simulated water quality in the three main basins of the lake. Five terrestrial nutrient control strategies were evaluated. Effectiveness differed between catchments, and water quality responses to nutrient load reductions varied between deep and shallow lake basins. Nutrient load reductions were a significant driver of increased DO concentrations, however strategies which reduced tributary inflow had a greater impact on lake restoration, associated with changes in water temperature and chemistry. Importantly, when multiple strategies were implemented simultaneously, resultant large flow reductions induced warming throughout the water column. Negative impacts of lake warming on DO overwhelmed the positive effects of nutrient reduction, and limited the effectiveness of lake restoration strategies. This study indicates that rates of lake recovery may be accelerated through a coordinated management approach, which considers strategy interactions, and the potential for temperature change-induced physical and biological feedbacks. Identified impacts of flow and temperature on rates of lake recovery have implications for management sustainability under a changing climate.

Modelling metaldehyde in catchments: a River Thames case-study.

The application of metaldehyde to agricultural catchment areas to control slugs and snails has caused severe problems for drinking water supply in recent years. In the River Thames catchment, metaldehyde has been detected at levels well above the EU and UK drinking water standards of 0.1 µg l-1 at many sites across the catchment between 2008 and 2015. Metaldehyde is applied in autumn and winter, leading to its increased concentrations in surface waters. It is shown that a process-based hydro-biogeochemical transport model (INCA-contaminants) can be used to simulate metaldehyde transport in catchments from areas of application to the aquatic environment. Simulations indicate that high concentrations in the river system are a direct consequence of excessive application rates. A simple application control strategy for metaldehyde in the Thames catchment based on model results is presented.

Fate and transport of polychlorinated biphenyls (PCBs) in the River Thames catchment - Insights from a coupled multimedia fate and hydrobiogeochemical transport model.

The fate of persistent organic pollutants (POPs) in riverine environments is strongly influenced by hydrology (including flooding) and fluxes of sediments and organic carbon. Coupling multimedia fate models (MMFMs) and hydrobiogeochemical transport models offers unique opportunities for understanding the environmental behaviour of POPs. While MMFMs are widely used for simulating the fate and transport of legacy and emerging pollutants, they use greatly simplified representations of climate, hydrology and biogeochemical processes. Using additional information about weather, river flows and water chemistry in hydrobiogeochemical transport models can lead to new insights about POP behaviour in rivers. As most riverine POPs are associated with suspended sediments (SS) or dissolved organic carbon (DOC), coupled models simulating SS and DOC can provide additional insights about POPs behaviour. Coupled simulations of river flow, DOC, SS and POP dynamics offer the possibility of improved predictions of contaminant fate and fluxes by leveraging the additional information in routine water quality time series. Here, we present an application of a daily time step dynamic coupled multimedia fate and hydrobiogeochemical transport model (The Integrated Catchment (INCA) Contaminants model) to simulate the behaviour of selected PCB congeners in the River Thames (UK). This is a follow-up to an earlier study where a Level III fugacity model was used to simulate PCB behaviour in the Thames. While coupled models are more complex to apply, we show that they can lead to much better representation of POPs dynamics. The present study shows the importance of accurate sediment and organic carbon simulations to successfully predict riverine PCB transport. Furthermore, it demonstrates the important impact of short-term weather variation on PCB movement through the environment. Specifically, it shows the consequences of the severe flooding, which occurred in early 2014 on sediment PCB concentrations in the River Thames.

An INCA model for pathogens in rivers and catchments: Model structure, sensitivity analysis and application to the River Thames catchment, UK.

Pathogens are an ongoing issue for catchment water management and quantifying their transport, loss and potential impacts at key locations, such as water abstractions for public supply and bathing sites, is an important aspect of catchment and coastal management. The Integrated Catchment Model (INCA) has been adapted to model the sources and sinks of pathogens and to capture the dominant dynamics and processes controlling pathogens in catchments. The model simulates the stores of pathogens in soils, sediments, rivers and groundwaters and can account for diffuse inputs of pathogens from agriculture, urban areas or atmospheric deposition. The model also allows for point source discharges from intensive livestock units or from sewage treatment works or any industrial input to river systems. Model equations are presented and the new pathogens model has been applied to the River Thames in order to assess total coliform (TC) responses under current and projected future land use. A Monte Carlo sensitivity analysis indicates that the input coliform estimates from agricultural sources and decay rates are the crucial parameters controlling pathogen behaviour. Whilst there are a number of uncertainties associated with the model that should be accounted for, INCA-Pathogens potentially provides a useful tool to inform policy decisions and manage pathogen loading in river systems.

Rainfall runoff modelling of the Upper Ganga and Brahmaputra basins using PERSiST.

There are ongoing discussions about the appropriate level of complexity and sources of uncertainty in rainfall runoff models. Simulations for operational hydrology, flood forecasting or nutrient transport all warrant different levels of complexity in the modelling approach. More complex model structures are appropriate for simulations of land-cover dependent nutrient transport while more parsimonious model structures may be adequate for runoff simulation. The appropriate level of complexity is also dependent on data availability. Here, we use PERSiST; a simple, semi-distributed dynamic rainfall-runoff modelling toolkit to simulate flows in the Upper Ganges and Brahmaputra rivers. We present two sets of simulations driven by single time series of daily precipitation and temperature using simple (A) and complex (B) model structures based on uniform and hydrochemically relevant land covers respectively. Models were compared based on ensembles of Bayesian Information Criterion (BIC) statistics. Equifinality was observed for parameters but not for model structures. Model performance was better for the more complex (B) structural representations than for parsimonious model structures. The results show that structural uncertainty is more important than parameter uncertainty. The ensembles of BIC statistics suggested that neither structural representation was preferable in a statistical sense. Simulations presented here confirm that relatively simple models with limited data requirements can be used to credibly simulate flows and water balance components needed for nutrient flux modelling in large, data-poor basins.

Assessing the impacts of climate change and socio-economic changes on flow and phosphorus flux in the Ganga river system.

Anthropogenic climate change has impacted and will continue to impact the natural environment and people around the world. Increasing temperatures and altered rainfall patterns combined with socio-economic factors such as population changes, land use changes and water transfers will affect flows and nutrient fluxes in river systems. The Ganga river, one of the largest river systems in the world, supports approximately 10% global population and more than 700 cities. Changes in the Ganga river system are likely to have a significant impact on water availability, water quality, aquatic habitats and people. In order to investigate these potential changes on the flow and water quality of the Ganga river, a multi-branch version of INCA Phosphorus (INCA-P) model has been applied to the entire river system. The model is used to quantify the impacts from a changing climate, population growth, additional agricultural land, pollution control and water transfers for 2041-2060 and 2080-2099. The results provide valuable information about potential effects of different management strategies on catchment water quality.

Impacts of climate change and socio-economic scenarios on flow and water quality of the Ganges, Brahmaputra and Meghna (GBM) river systems: low flow and flood statistics.

The potential impacts of climate change and socio-economic change on flow and water quality in rivers worldwide is a key area of interest. The Ganges-Brahmaputra-Meghna (GBM) is one of the largest river basins in the world serving a population of over 650 million, and is of vital concern to India and Bangladesh as it provides fresh water for people, agriculture, industry, conservation and for the delta system downstream. This paper seeks to assess future changes in flow and water quality utilising a modelling approach as a means of assessment in a very complex system. The INCA-N model has been applied to the Ganges, Brahmaputra and Meghna river systems to simulate flow and water quality along the rivers under a range of future climate conditions. Three model realisations of the Met Office Hadley Centre global and regional climate models were selected from 17 perturbed model runs to evaluate a range of potential futures in climate. In addition, the models have also been evaluated using socio-economic scenarios, comprising (1) a business as usual future, (2) a more sustainable future, and (3) a less sustainable future. Model results for the 2050s and the 2090s indicate a significant increase in monsoon flows under the future climates, with enhanced flood potential. Low flows are predicted to fall with extended drought periods, which could have impacts on water and sediment supply, irrigated agriculture and saline intrusion. In contrast, the socio-economic changes had relatively little impact on flows, except under the low flow regimes where increased irrigation could further reduce water availability. However, should large scale water transfers upstream of Bangladesh be constructed, these have the potential to reduce flows and divert water away from the delta region depending on the volume and timing of the transfers. This could have significant implications for the delta in terms of saline intrusion, water supply, agriculture and maintaining crucial ecosystems such as the mangrove forests, with serious implications for people's livelihoods in the area. The socio-economic scenarios have a significant impact on water quality, altering nutrient fluxes being transported into the delta region.

Dynamic modeling of the Ganga river system: impacts of future climate and socio-economic change on flows and nitrogen fluxes in India and Bangladesh.

This study investigates the potential impacts of future climate and socio-economic change on the flow and nitrogen fluxes of the Ganga river system. This is the first basin scale water quality study for the Ganga considering climate change at 25 km resolution together with socio-economic scenarios. The revised dynamic, process-based INCA model was used to simulate hydrology and water quality within the complex multi-branched river basins. All climate realizations utilized in the study predict increases in temperature and rainfall by the 2050s with significant increase by the 2090s. These changes generate associated increases in monsoon flows and increased availability of water for groundwater recharge and irrigation, but also more frequent flooding. Decreased concentrations of nitrate and ammonia are expected due to increased dilution. Different future socio-economic scenarios were found to have a significant impact on water quality at the downstream end of the Ganga. A less sustainable future resulted in a deterioration of water quality due to the pressures from higher population growth, land use change, increased sewage treatment discharges, enhanced atmospheric nitrogen deposition, and water abstraction. However, water quality was found to improve under a more sustainable strategy as envisaged in the Ganga clean-up plan.

Uncertainty assessments and hydrological implications of climate change in two adjacent agricultural catchments of a rapidly urbanizing watershed.

Lake Simcoe is the most important inland lake in Southern Ontario. The watershed is predominantly agricultural and under increasing pressure from urbanization, leading to changing runoff patterns in rivers draining to the lake. Uncertainties in rainfall-runoff modeling in tributary catchments of the Lake Simcoe Watershed (LSW) can be an order of magnitude larger than pristine watersheds, hampering water quality predictions and export calculations. Here we conduct a robust assessment to constrain the uncertainty in hydrological simulations and projections in the LSW using two representative adjacent agricultural catchments. Downscaled CGCM 3 projections using A1B and A2 emission scenarios projected increases of 4°C in air temperature and a 26% longer growing season. The fraction of precipitation falling as snow will decrease. Spring runoff is an important event in LSW but individual HBV best calibrated parameter sets under-predicted peak flows by up to 32%. Using an ensemble of behavioral parameter sets achieved credible representations of present day hydrology and constrained uncertainties in future projections. Parameter uncertainty analysis showed that the catchments differ in terms of their snow accumulation/melt and groundwater dynamics. Human activities exacerbate the differences in hydrological response. Model parameterization in one catchment could not generate credible hydrological simulations in the other. We cautioned against extrapolating results from monitored to ungauged catchments in managed watersheds like the LSW.

A cost-effectiveness analysis of water security and water quality: impacts of climate and land-use change on the River Thames system.

The catchment of the River Thames, the principal river system in southern England, provides the main water supply for London but is highly vulnerable to changes in climate, land use and population. The river is eutrophic with significant algal blooms with phosphorus assumed to be the primary chemical indicator of ecosystem health. In the Thames Basin, phosphorus is available from point sources such as wastewater treatment plants and from diffuse sources such as agriculture. In order to predict vulnerability to future change, the integrated catchments model for phosphorus (INCA-P) has been applied to the river basin and used to assess the cost-effectiveness of a range of mitigation and adaptation strategies. It is shown that scenarios of future climate and land-use change will exacerbate the water quality problems, but a range of mitigation measures can improve the situation. A cost-effectiveness study has been undertaken to compare the economic benefits of each mitigation measure and to assess the phosphorus reductions achieved. The most effective strategy is to reduce fertilizer use by 20% together with the treatment of effluent to a high standard. Such measures will reduce the instream phosphorus concentrations to close to the EU Water Framework Directive target for the Thames.

The interactive responses of water quality and hydrology to changes in multiple stressors, and implications for the long-term effective management of phosphorus.

Soluble reactive phosphorus (SRP) plays a key role in eutrophication, a global problem decreasing habitat quality and in-stream biodiversity. Mitigation strategies are required to prevent SRP fluxes from exceeding critical levels, and must be robust in the face of potential changes in climate, land use and a myriad of other influences. To establish the longevity of these strategies it is therefore crucial to consider the sensitivity of catchments to multiple future stressors. This study evaluates how the water quality and hydrology of a major river system in the UK (the River Thames) respond to alterations in climate, land use and water resource allocations, and investigates how these changes impact the relative performance of management strategies over an 80-year period. In the River Thames, the relative contributions of SRP from diffuse and point sources vary seasonally. Diffuse sources of SRP from agriculture dominate during periods of high runoff, and point sources during low flow periods. SRP concentrations rose under any future scenario which either increased a) surface runoff or b) the area of cultivated land. Under these conditions, SRP was sourced from agriculture, and the most effective single mitigation measures were those which addressed diffuse SRP sources. Conversely, where future scenarios reduced flow e.g. during winters of reservoir construction, the significance of point source inputs increased, and mitigation measures addressing these issues became more effective. In catchments with multiple point and diffuse sources of SRP, an all-encompassing effective mitigation approach is difficult to achieve with a single strategy. In order to attain maximum efficiency, multiple strategies might therefore be employed at different times and locations, to target the variable nature of dominant SRP sources and pathways.

Problems with the reconciliation of good ecological status and public participation in the Water Framework Directive.

The Water Framework Directive (WFD) is an ambitious piece of legislation focused on achieving good ecological status as defined by deviations from reference conditions. Achieving good ecological status depends on collaboration between stakeholders, scientists and the public. However, public participation is restricted to consultations about implementing measures to achieve good ecological status, not in the goal setting. There are multiple, competing interpretations of good ecological status. This study addresses two of the pillars of the WFD, good ecological status and public participation. We argue that these two pillars are currently at odds when defining reference conditions for surface waters, and it is unclear how they can work together in practice. We also contend that there is an intention in the WFD to integrate these two pillars, but there is no legal support for their connection. In a case study of a small boreal lake in Sweden, we show that local people possess a great deal of historical knowledge, which they use to conceptualize reference conditions. Their conceptualizations are compared with fish and water chemistry monitoring by the regulatory authority as well as paleolimnological reconstructions of water quality dating back to the beginning of the 20th century. The knowledge that the local people have corresponds to the historical data available for the lake, particularly with water clarity. We highlight the subjective nature of the concept of 'undisturbed state' to show that it varies depending on values, knowledge and perceptions of lay-people, scientists and relevant authorities. The subjectivity of the concept of undisturbed state promises to be a way of linking the two pillars of the WFD.

Using the INCA-Hg model of mercury cycling to simulate total and methyl mercury concentrations in forest streams and catchments.

We present a new, catchment-scale, process-based dynamic model for simulating mercury (Hg) in soils and surface waters. The Integrated Catchments Model for Mercury (INCA-Hg) simulates transport of gaseous, dissolved and solid Hg and transformations between elemental (Hg(0)), ionic (Hg(II)) and methyl (MeHg) Hg in natural and semi-natural landscapes. The mathematical description represents the model as a series of linked, first-order differential equations describing chemical and hydrological processes in catchment soils and waters which we believe control surface water Hg dynamics. The model simulates daily time series between one and 100 years long and can be applied to catchments ranging in size from <1 to ~10,000 km(2). Here we present applications of the model to two boreal forest headwater catchments in central Canada where we were able to reproduce observed patterns of stream water total mercury (THg) and MeHg fluxes and concentrations. Model performance was assessed using Monte Carlo techniques. Simulated in-stream THg and MeHg concentrations were sensitive to hydrologic controls and terrestrial and aquatic process rates.

Modelling the long term impact of climate change on the carbon budget of Lake Simcoe, Ontario using INCA-C.

This study presents a process-based model of dissolved organic carbon concentration ([DOC]) in catchments draining into Lake Simcoe, Ontario. INCA-C, the Integrated Catchment model for Carbon, incorporates carbon biogeochemical processes in a terrestrial system with hydrologic flow paths to simulate watershed wide [DOC]. The model successfully simulates present-day inter-annual and seasonal [DOC] dynamics in tributaries draining catchments with mixed or contrasting land cover in the Lake Simcoe watershed (LSW). The sensitivity of INCA-C to soil moisture, hydrologic controls and land uses within a watershed demonstrates its significance as a tool to explore pertinent environmental issues specific to the LSW. Projections of climate change under A1B and A2 SRES scenarios suggest a continuous monotonic increase in [DOC] in surface waters draining into Lake Simcoe. Large variations in seasonal DOC dynamics are predicted to occur during summer with a possibility of displacement of summer [DOC] maxima towards winter and a prolongation of summer [DOC] levels into the autumn. INCA-C also predicts possible increases in dissolved inorganic carbon in some tributaries with rising temperature suggesting increased CO(2) emissions from rivers as climate changes.

Assessment of uncertainty in long-term mass balances for acidification assessments: a MAGIC model exercise.

Long-term (1860-2010) catchment mass balance calculations rely on models and assumptions which are sources of uncertainty in acidification assessments. In this article, we report on an application of MAGIC to model acidification at the four Swedish IM forested catchments that have been subject to differing degrees of acidification stress. Uncertainties in the modeled mass balances were mainly associated with the deposition scenario and assumptions about sulfate adsorption and soil mass. Estimated base cation (BC) release rates (weathering) varied in a relatively narrow range of 47-62 or 42-47 meq m(-2) year(-1), depending on assumptions made about soil cation exchange capacity and base saturation. By varying aluminum solubility or introducing a dynamic weathering feedback that allowed BC release to increase at more acidic pHs, a systematic effect on predicted changes in acid neutralizing capacity (ΔANC ca. 10-41 µeq l(-1)) and pH (ca. ΔpH = 0.1-0.6) at all sites was observed. More robust projections of future changes in pH and ANC are dependent on reducing uncertainties in BC release rates, the timing, and extent of natural acidification through BC uptake by plants, temporal changes in soil element pools, and fluxes of Al between compartments.

Simulating dissolved organic carbon dynamics at the swedish integrated monitoring sites with the integrated catchments model for carbon, INCA-C.

Surface water concentrations of dissolved organic carbon ([DOC]) are changing throughout the northern hemisphere due to changes in climate, land use and acid deposition. However, the relative importance of these drivers is unclear. Here, we use the Integrated Catchments model for Carbon (INCA-C) to simulate long-term (1996-2008) streamwater [DOC] at the four Swedish integrated monitoring (IM) sites. These are unmanaged headwater catchments with old-growth forests and no major changes in land use. Daily, seasonal and long-term variations in streamwater [DOC] driven by runoff, seasonal temperature and atmospheric sulfate (SO4(2-)) deposition were observed at all sites. Using INCA-C, it was possible to reproduce observed patterns of variability in streamwater [DOC] at the four IM sites. Runoff was found to be the main short-term control on [DOC]. Seasonal patterns in [DOC] were controlled primarily by soil temperature. Measured SO4(2-) deposition explained some of the long-term [DOC] variability at all sites.