Projected stream water fluxes of NO3 and total organic carbon from the Storgama headwater catchment, Norway, under climate change and reduced acid deposition

Fluctuations in the 20-year record of nitrate (NO3) and total organic carbon (TOC) concentrations and fluxes in runoff at the small headwater catchment Storgama, southern Norway, were related to climate and acid deposition. The long-term decline in NO3 related to reduced NO3 deposition and increased winter discharge, whereas the long-term increase in TOC related to reduced sulfur deposition. Multiple regression models describing long-term trends and seasonal variability in these records were used to project future concentrations given scenarios of climate change and acid deposition. All scenarios indicated reduced NO3 fluxes and increased TOC fluxes; the largest projected changes for the period 2071-2100 were -86% and +24%, respectively. Uncertainties are that the predicted future temperatures are considerably higher than the historical record. Also, nonlinear responses of ecosystem processes (nitrogen [N] mineralization) to temperature, N-enrichment of soils, and step-changes in environmental conditions may affect future leaching of carbon and N.     

Use of modeling to protect,plan, and manage water resources in catchment areas

The degradation of water resources by diffuse pollution, mainly due to nitrate and pesticides, is an important matter for public health. Restoration of the quality of natural water catchments by focusing on their catchment areas is therefore a national priority in France. To consider catchment areas as homogeneous and to expend an equal effort on the entire area inevitably leads to a waste of time and money, and restorative actions may not be as efficient as intended. The variability of the pedological and geological properties of the area is actually an opportunity to invest effort on smaller areas, simply because every action is not equally efficient on every kind of pedological or geological surface. Using this approach, it is possible to invest in a few selected zones that will be efficient in terms of environmental results. The contributive hydraulic areas (CHA) concept is different from that of the catchment area. Because the transport of most of the mobile and persistent pollutants is primarily driven by water circulation, the concept of the CHA is based on the water pathway from the surface of the soil in the catchment area to the well. The method uses a three-dimensional hydrogeological model of surface and groundwater integrated with a geographic information system called Watermodel. The model calculates the contribution (m³/h or %) of each point of the soil to the total flow pumped in a well. Application of this model, partially funded by the Seine Normandy Basin Agency, to the catchment of the Dormelles Well in the Cretaceous chalk aquifer in the Orvanne valley, France (catchment area of 23,000 ha at Dormelles, county 77), shows that 95 % of the water pumped at the Dormelles Well comes from only 26 % of the total surface area of the catchment. Consequently, an action plan to protect the water resource will be targeted at the 93 farmers operating in this source area rather than the total number of farmers (250) across the entire 23,000 ha. Another model, developed from Epiclès© software, permits the calculation of the under-root nitrate concentrations for each field based on soil type, climate, and farming practices. When the Watermodel and Epiclès© are coupled, nitrate transfers from the soil to the catchment and the river can be modeled. In this study, the initial pollution due to the actual farming practices was simulated and we were also able to estimate the efficiency of the agronomic action plan by testing several scenarios and calculating the time needed to reach the target nitrate concentration in the well.     

Integrated catchment modeling for nutrient reduction: scenarios showing impacts, potential, and cost of measures

A hydrological-based model (HBV-NP) was applied to a catchment (1900 km2) in the southern part of Sweden. Careful characterization of the present load situation and the potential for improved treatment or reduced soil leaching were analyzed. Several scenarios were modeled to find strategies to reach the Swedish environmental goals of reducing anthropogenic nitrogen load by 30% and phosphorus load by 20%. It was stated that the goals could be reached by different approaches that would affect different polluters and social sectors. However, no single measure was enough by itself. Instead, a combination of measures was necessary to achieve the goals. The nitrogen goal was the most difficult to attain. In order to be cost-effective, these measures should be applied to areas contributing the most to the net loading of the sea. This strategy could reduce the costs by 70%-80% when compared with implementing the measures in the entire catchment. Integrated catchment models may thus be helpful tools for reducing costs in environmental control programs.     

Nitrogen dynamics of a mountain forest on dolomitic limestone--a scenario-based risk assessment

The dominant nitrogen (N) fluxes were simulated in a mountain forest ecosystem on dolomitic bedrock in the Austrian Alps. Based on an existing small-scale climate model the simulation encompassed the present situation and a 50-yr projection. The investigated scenarios were current climate, current N deposition (SC1) and future climate (+2.5 degrees C and +10% annual precipitation) with three levels of N deposition (SC2, 3, 4). The microbially mediated N transformation, including the emission of nitrogen oxides, was calculated with PnET-N-DNDC. Soil hydrology was calculated with HYDRUS and was used to estimate the leaching of nitrate. The expected change of the forest ecosystem due to changes of the climate and the N availability was simulated with PICUS. The incentive for the project was the fact that forests on dolomitic limestone stock on shallow Rendzic Leptosols that are rich in soil organic matter are considered highly sensitive to the expected environmental changes. The simulation results showed a strong effect due to increased temperatures and to elevated levels of N deposition. The outflux of N, both as nitrate (6-25kg Nha(-1)yr(-1)) and nitrogen oxides (1-2kg Nha(-1)yr(-1)), from the forest ecosystem are expected to increase. Temperature exerts a stronger effect on the N(2)O emission than the increased rate of N deposition. The main part of the N emission will occur as N(2) (15kg Nha(-1)yr(-1)). The total N loss is partially offset by increased rates of N uptake in the biomass due to an increase in forest productivity.     

Impact of Climate Change on Ecological Quality Indicators and Biogeochemical Fluxes in the Baltic Sea: A Multi-Model Ensemble Study

Multi-model ensemble simulations using three coupled physical-biogeochemical models were performed to calculate the combined impact of projected future climate change and plausible nutrient load changes on biogeochemical cycles in the Baltic Sea. Climate projections for 1961-2099 were combined with four nutrient load scenarios ranging from a pessimistic business-as-usual to a more optimistic case following the Helsinki Commission's (HELCOM) Baltic Sea Action Plan (BSAP). The model results suggest that in a future climate, water quality, characterized by ecological quality indicators like winter nutrient, summer bottom oxygen, and annual mean phytoplankton concentrations as well as annual mean Secchi depth (water transparency), will be deteriorated compared to present conditions. In case of nutrient load reductions required by the BSAP, water quality is only slightly improved. Based on the analysis of biogeochemical fluxes, we find that in warmer and more anoxic waters, internal feedbacks could be reinforced. Increased phosphorus fluxes out of the sediments, reduced denitrification efficiency and increased nitrogen fixation may partly counteract nutrient load abatement strategies.     

Climate change impacts on the leaching of a heavy metal contamination in a small lowland catchment

The Keersop catchment (43km(2)) in the south of The Netherlands has been contaminated by the emissions of four zinc ore smelters. The objective of this study was to assess the effects of future projected climate change on the hydrology and the leaching of heavy metals (i.e. Cd and Zn) in the catchment. The numerical, quasi-2D, unsaturated zone Soil Water Atmosphere Plant model was used with 100-year simulated daily time series of precipitation and potential evapotranspiration. The time series are representative of stationary climates for the periods 1961-1990 ("baseline") and 2071-2100 ("future"). The time series of future climate were obtained by downscaling the results of eight regional climate model (RCM) experiments, driven by the SRES A2 emissions scenario, using change factors for a series of climate statistics and applying them to stochastic weather generator models. The time series are characterized by increased precipitation in winter, less precipitation in summer, and higher air temperatures (between 2°C and 5°C) throughout the year. Future climate scenarios project higher evapotranspiration rates, more irrigation, less drainage, lower discharge rates and lower groundwater levels, due to increased evapotranspiration and a slowing down of the groundwater system. As a result, lower concentrations of Cd and Zn in surface water are projected. The reduced leaching of heavy metals, due to drying of the catchment, showed a positive impact on a limited aspect of surface water quality.     

Modeling the impact of potential wetlands on phosphorus retention in a Swedish catchment

In southern Sweden, wetlands are constructed to remove nitrogen (N) in agricultural catchments. The possible effects of such wetlands on riverine phosphorus (P) were also estimated using input-output data from three well-monitored wetlands. This was done to formulate a simple model for removal of P that is dependent on inflow characteristics. Next, the N- and P-reducing effects of wetlands were modeled on a catchment scale (1900 km2) using the HBV-NP model and various assumptions about the wetland area and location. All three wetlands functioned as sinks for total P (tot-P) and for total suspended solids (TSS) with a removal of 10% to 31% and 28% to 50%, respectively. Mean P-removal rates of 17-49 kg ha(-1) yr(-1) were well simulated with the model. Catchment scale simulations indicated that wetlands were more efficient (in percentage of load) as traps for P than for N and that this may motivate the construction of wetlands for P removal far upstream from the catchment outlet.     

Recovery from eutrophication: experiences of reduced phosphorus input to the four largest lakes of Sweden

In-lake concentration changes of phosphorus (P) and nitrogen (N) in lakes V?ttern, V?nern, M?laren, and Hj?lmaren in response to diminished input has been examined from the mid-1960s onwards. In the former two deep and oligotrophic lakes with slow water renewal, drastic reductions in P-input from the middle of 1970s caused just minor reductions in P-concentration over a very long time. At the same time accumulation occurred in the water mass of inorganic N and possible reasons are discussed. In the latter two mesotrophic to hypertrophic lakes, two shallow basins in L. Hj?lmaren showed slow recovery due to release of P from sediments. The same basins and two basins in L. M?laren have suffered from N-deficiency, particularly during the pre-phosphorus reduction years, and nitrogen fixation was indicated. In two L. M?laren basins recovery of in-lake P concentrations was better than expected in comparison to the so-called IMSA-model for lake recovery from nutrient pollution. In the other five lakes/basins chlorophyll concentrations after 20 years were similar compared to those modeled.     

Spatial patterns concentrations in upland Wales in relation to catchment forest cover and forest age

Data on nitrate nitrogen were collected weekly during 1984 from 136 sites on streams in upland Wales. Mean nitrate concentrations in summer (0.02 to 1.5 mg litre(-1)) were significantly lower (P<0.001) than in winter (0.02 to 1.26 mg litre(-1)), particularly at sites with mature conifers (>30 years old). Mean concentrations increased significantly with the average age of conifers on each catchment (P<0.001), and with increasing areal cover by trees over 30 years old (P<0.001). Nitrate concentrations increased significantly with stream total hardness (P<0.001), possibly reflecting nitrogen mineralisation in soils of higher base status. Concentrations also increased with stream chloride (P<0.001), which is predominantly atmospherically derived, implying that increased nitrate occurred where general atmospheric inputs of solutes were increased. After accounting for variation in hardness, residual nitrate concentrations still increased with the average age of the conifers (P<0.001), and with catchment cover by mature trees (P<0.001). We infer that some additional nitrate under older conifers is thus independent of catchment sources associated with increasing hardness. Two possibilities are increased inputs and decreased retention of nitrogen within the ecosystem of maturing conifer forest. Residual nitrate after accounting for variations in chloride also increased significantly with conifer age (P<0.01) and cover (P<0.01), a pattern implying that some sources of nitrate may also be independent of increased sea-salt deposition. We allude to the possibility that additional nitrogen deposition adds to nitrogen throughputs from maturing forests, and we discuss the potential ecological role of additional nitrogen in runoff.     

Parameter precision in the HBV-NP model and impacts on nitrogen scenario simulations in the Rönneå river, southern Sweden

The HBV-NP model is a newly developed water quality model that describes the turnover and fluxes of both nitrogen and phosphorous. It is based on the conceptual precipitation/runoff HBV model. The HBV-NP model was applied for simulation of nitrogen for the R?nne? catchment in southern Sweden. The catchment was divided into 64 subcatchments in the model. Discharge measurements from six stations and nitrogen measurements from 12 stations were used in the calibration of parameters in the model. Eight automatic calibrations were performed with different combinations of time periods, objective functions, and levels of the nitrogen load in the model. A regionally extended interpretation of the Nash-Sutcliffe R2 criterion was used in the calibration. In the evaluation of the criterion, the errors were summed over both time steps and sampling points. Scenario simulations of combined measures for reduction of nitrogen load into the sea by 30% were thereafter performed with the eight sets of parameters established by calibration. The model parameters were not uniquely defined by the calibration. However, the simulated relative reduction of nitrogen load into the sea was relatively insensitive to the choice of parameter set, given the available input sources, variables, and data.     

Salting our landscape: an integrated catchment model using readily accessible data to assess emerging road salt contamination to streams

A new integrated catchment model for salinity has been developed to assess the transport of road salt from upland areas in watersheds to streams using readily accessible landscape, hydrologic, and meteorological data together with reported salt applications. We used Fishkill Creek (NY) as a representative watershed to test the model. Results showed good agreement between modeled and measured stream water chloride concentrations. These results suggest that a dominant mode of catchment simulation that does not entail complex deterministic modeling is an appropriate method to model salinization and to assess effects of future applications of road salt to streams. We heuristically increased and decreased salt applications by 100% and results showed that stream chloride concentrations increased by 13% and decreased by 7%, respectively. The model suggests that future management of salt application can reduce environmental concentrations, albeit over some time.     

Long-term nutrient enrichment of an estuarine system: Ythan, Scotland (1958-1993)

A comprehensive review is given of past and current annually averaged nutrient concentrations (total oxidised nitrogen [TON], phosphate and silicate) in the River Ythan and its estuary. TON concentrations in the River Ythan have increased from ca 100-150 microm in the late 1960s to ca 500-550 microm in the early 1990s, changes which are also shown in the estuary. The interpretation of the phosphate data is more complex but concentrations in the estuary also appear to have increased. Silicate concentrations have remained constant with time. The results of 11 surveys of the estuary completed during 1993 are given. TON and silicate behave conservatively in the estuary, with concentrations determined by those of the fresh water end-member and hence by processes within the catchment. Phosphate and ammonia concentrations at low salinity are affected by the discharge from the Ellon sewage treatment facility at the head of the estuary. Throughout the estuary there is evidence of an input of ammonia to the water column from the sediments, but it represents only ca 0.2-2% of the nitrogen entering from the River Ythan. Despite the increased nutrient loading to the estuary there is no evidence that this has resulted in higher chlorophyll concentrations in the water column, this observation is attributed to the relatively short flushing period of the estuary.     

Oryzalin fate and transport in runoff water in Mediterranean vineyards

An experimental study was conducted in a 91.4-ha Mediterranean vineyard catchment in southern France to characterize the fate and transport of oryzalin in runoff water and thus to assess the risk of contamination of surface waters. Oryzalin concentrations in soil were monitored on two fields, one no-till and one tilled from March 1998 to March 2000. Concentrations in solution and on solid phase of runoff water were measured at the outlets of both fields and the catchment. The droughts in the two summer periods reduced the dissipation of oryzalin and increased its field half-life up to 35 days. Consequently, oryzalin was detected throughout the year in runoff water, with maximum dissolved concentrations > 600 microg l(-1) at the field scale. Oryzalin transport essentially occurred in solution. At the no-till field, seasonal losses were 2.29% and 1.89% of the applied amount in 1998 and 1999, respectively. The corresponding values at the tilled field were 1.56% and 0.29%, since tillage reduced total losses by reducing surface runoff. At the catchment scale, oryzalin concentrations were smaller than those at the field scale, due to dilution effects and staggering of application. Large part of the overland flow from the fields reinfiltrated in the ditches before reaching the outlet of the catchment. As a result, seasonal oryzalin losses were <0.2% of the applied amount.     

Estimating catchment nutrient flow with the HBV-NP model: sensitivity to input data

The dynamic catchment model HBV-N has been further developed by adding routines for phosphorus transport and is now called the HBV-NP model. The model was shown to satisfactorily simulate nutrient dynamics in the R?nne? catchment (1,900 km2). Its sensitivity to input data was tested, and results demonstrated the increased sensitivity to the selection of input data on a subcatchment scale when compared with the catchment scale. Selection of soil and land use databases was found to be critical in some subcatchments but did not have a significant impact on a catchment scale. Although acceptable on a catchment scale, using templates and generalization, with regards to emissions from point sources and rural households, significantly decreased model performance in certain subcatchments when compared with using more detailed local information. A division into 64 subcatchments resulted in similar model performance at the catchment outlet when compared with a lumped approach. Adjusting the imported matrixes of the regional leaching of nitrogen, from agricultural land, against mean subcatchment water percolation did not have a significant impact on the model performance.     

Surface runoff and nitrogen (N) loss in a bamboo (Phyllostachys pubescens) forest under different fertilization regimes

Nitrogen (N) losses from agricultural fields have been extensively studied. In contrast, surface runoff and N losses have rarely been considered for bamboo forests that are widespread in regions such as southern China. The thriving of bamboo industries has led to increasing fertilizer use in bamboo forests. In this study, we evaluated surface runoff and N losses in runoff following different fertilization treatments under field conditions in a bamboo (Phyllostachys pubescens) forest in the catchment of Lake Taihu in Jiangsu, China. Under three different fertilization regimes, i.e., control, site-specific nutrient management (SSNM), and farmer's fertilization practice (FFP), the water runoff rate amounted to 356, 361, and 342 m³?ha^?1 and accounted for 1.91, 1.98, and 1.85 % of the water input, respectively, from June 2009 to May 2010. The total N losses via surface runoff ranged from 1.2 to 1.8 kg?ha^?1. Compared with FFP, the SSNM treatment reduced total nitrogen (TN) and dissolved nitrogen (DN) losses by 31 and 34 %, respectively. The results also showed that variations in N losses depended mainly on runoff fluxes, not N concentrations. Runoff samples collected from all treatments throughout the year showed TN concentrations greater than 0.35 mg?L^?1, with the mean TN concentration in the runoff from the FFP treatment reaching 8.97 mg?L^?1. The loss of NO₃ ^?–N was greater than the loss of NH₄ ⁺–N. The total loss of dissolved organic nitrogen (DON) reached 23–41 % of the corresponding DN. Therefore, DON is likely the main N species in runoff from bamboo forests and should be emphasized in the assessment and management of N losses in bamboo forest.     

Manipulation of precipitation in small headwater catchments at Storgama, Norway: effects on leaching of organic carbon and nitrogen species

Projected changes in climate in Southern Norway include increases in summer and autumn precipitation. This may affect leaching of dissolved organic matter (DOM) from soils. Effects of experimentally added extra precipitation (10 mm week) during the growing season of 3 years (2004-2006) to small headwater catchments at Storgama (59 degrees 0'N, 550-600 m a.s.l.) on leaching of total organic carbon (TOC) and total organic nitrogen (TON) were assessed. Extra precipitation did not have a significant effect on average TOC and TON concentrations in runoff. Thus, fluxes of TOC and TON increased nearly proportionally with water fluxes. This suggests that a store of adsorbed and potentially mobile TOC and TON in catchment soils buffers the concentration of DOM in runoff. The size and dynamics of the pool of TOC and TON depends on the balance between production and leaching rates. Infrequent short droughts had only small effects on TOC and TON fluxes in runoff from the reference catchments.     

Increased nitrogen in runoff and soil following 13 years of experimentally increased nitrogen deposition to a coniferous-forested catchment at Gårdsjön, Sweden

Beginning in 1991, we have added nitrogen (N) to the 0.5-ha, N-poor, coniferous-forested catchment G2 NITREX at G?rdsj?n, Sweden, to investigate the consequences of chronic elevated N deposition. We have added 40 kg N ha-1 yr-1 in fortnightly doses of NH4NO3 to the ambient 15 kg N ha-1 yr-1 by means of a sprinkling system. NO3 concentrations in runoff increased during 13 years from<1 to 70 microeq L-1, and in 2004 comprised about 10% of N input. Inhibition of NO3 immobilisation due to increased availability of NH4 might explain the increased leaching of NO3. C and N pools in the forest floor increased but C/N ratio has not changed. The increase in NO3 leaching thus occurred independently of change in C/N ratio. The results from G?rdsj?n demonstrate that increased leaching of inorganic N and decrease in C/N ratio respond to increased N deposition at greatly different time scales.     

Scale-dependence of land use effects on water quality of streams in agricultural catchments

The influence of land use on water quality in streams is scale-dependent and varies in time and space. In this study, land cover patterns and stocking rates were used as measures of agricultural development in two pasture and one native grassland catchment in New Zealand and were related to water quality in streams of various orders. The amount of pasture per subcatchment correlated well to total nitrogen and nitrate in one catchment and turbidity and total phosphorous in the other catchment. Stocking rates were only correlated to total phosphorous in one pasture catchment but showed stronger correlations to ammonium, total phosphorous and total nitrogen in the other pasture catchment. Winter and spring floods were significant sources of nutrients and faecal coliforms from one of the pasture catchments into a wetland complex. Nutrient and faecal coliform concentrations were better predicted by pastural land cover in fourth-order than in second-order streams. This suggests that upstream land use is more influential in larger streams, while local land use and other factors may be more important in smaller streams. These temporal and spatial scale effects indicate that water-monitoring schemes need to be scale-sensitive.     

Riparian Zone Influence on Stream Water Dissolved Organic Carbon Concentrations at the Swedish Integrated Monitoring Sites

Short-term variability in stream water dissolved organic carbon (DOC) concentrations is controlled by hydrology, climate and atmospheric deposition. Using the Riparian flow-concentration Integration Model (RIM), we evaluated factors controlling stream water DOC in the Swedish Integrated Monitoring (IM) catchments by separating out hydrological effects on stream DOC dynamics. Model residuals were correlated with climate and deposition-related drivers. DOC was most strongly correlated to water flow in the northern catchment (Gammtratten). The southern Aneboda and Kindla catchments had pronounced seasonal DOC signals, which correlated weakly to flow. DOC concentrations at Gårdsjön increased, potentially in response to declining acid deposition. Soil temperature correlated strongly with model residuals at all sites. Incorporating soil temperature in RIM improved model performance substantially (20–62% lower median absolute error). According to the simulations, the RIM conceptualization of riparian processes explains between 36% (Kindla) and 61% (Aneboda) of the DOC dynamics at the IM sites.     

A method to estimate the impact of clear-cutting on nutrient concentrations in boreal headwater streams

Large-scale forestry operations, like clear-cutting, may impair surface water quality if not done with environmental considerations in mind. Catchment and country level estimates of nutrient loads from forestry are generally based on specific export values, i.e., changes in annual exports due to the implemented forestry operations expressed in kg ha⁻¹. We introduce here a specific concentration approach as a method to estimate the impact of clear-cutting on nutrient concentrations and export in headwater streams. This new method is potentially a more dynamic and flexible tool to estimate nutrient loads caused by forestry, because variation in annual runoff can be taken into account in load assessments. We combined water quality data from eight boreal headwater catchment pairs located in Finland and Sweden, where the effect of clear-cutting on stream water quality has been studied experimentally. Statistically significant specific concentration values could be produced for total nitrogen, nitrate, ammonium, and phosphate. The significant increases in the concentrations of these nutrients occurred between 2 and 6 years after clear-cutting. Significant specific concentration values could not be produced for total phosphorus and total organic carbon with the whole dataset, although in some single studies significant increases in their concentrations after clear-cutting were observed. The presented method enables taking into account variation in runoff, temporal dynamics of effects, and the proportional size of the treated area in load calculations. The number of existing studies considering large site-specific variation in responses to clear-cutting is small, and therefore further empirical studies are needed to improve predictive capabilities of the specific concentration values.