Uncertainty assessment of a process-based integrated catchment model of phosphorus

Despite the many models developed for phosphorus concentration prediction at differing spatial and temporal scales, there has been little effort to quantify uncertainty in their predictions. Model prediction uncertainty quantification is desirable, for informed decision-making in river-systems management. An uncertainty analysis of the process-based model, integrated catchment model of phosphorus (INCA-P), within the generalised likelihood uncertainty estimation (GLUE) framework is presented. The framework is applied to the Lugg catchment (1,077 km2), a River Wye tributary, on the England–Wales border. Daily discharge and monthly phosphorus (total reactive and total), for a limited number of reaches, are used to initially assess uncertainty and sensitivity of 44 model parameters, identified as being most important for discharge and phosphorus predictions. This study demonstrates that parameter homogeneity assumptions (spatial heterogeneity is treated as land use type fractional areas) can achieve higher model fits, than a previous expertly calibrated parameter set. The model is capable of reproducing the hydrology, but a threshold Nash-Sutcliffe co-efficient of determination (E or R 2) of 0.3 is not achieved when simulating observed total phosphorus (TP) data in the upland reaches or total reactive phosphorus (TRP) in any reach. Despite this, the model reproduces the general dynamics of TP and TRP, in point source dominated lower reaches. This paper discusses why this application of INCA-P fails to find any parameter sets, which simultaneously describe all observed data acceptably. The discussion focuses on uncertainty of readily available input data, and whether such process-based models should be used when there isn’t sufficient data to support the many parameters.

Simulating metals and mine discharges in river basins using a new integrated catchment model for metals: pollution impacts and restoration strategies in the Aries-Mures river system in Transylvania, Romania

The INtegrated CAtchment (INCA) model has been developed to simulate the impact of mine discharges on river systems. The model accounts for the key kinetic chemical processes operating as well as the dilution, mixing and redistribution of pollutants in rivers downstream of mine discharges or acid rock drainage sites. The model is dynamic and simulates the day-to-day behaviour of hydrology and eight metals (cadmium, mercury, copper, zinc, lead, arsenic, manganese and chromium) as well as cyanide and ammonia. The model is semi-distributed and can simulate catchments, sub-catchment and in-stream river behaviour. The model has been applied to the Roia Montan Mine in Transylvania, Romania, and used to assess the impacts of old mine adits on the local catchments as well as on the downstream Aries and Mures river system. The question of mine restoration is investigated and a set of clean-up scenarios investigated. It is shown that the planned restoration will generate a much improved water quality from the mine and also alleviate the metal pollution of the river system.

Over-parameterised, uncertain 'mathematical marionettes' - How can we best use catchment water quality models? An example of an 80-year catchment-scale nutrient balance

The Integrated Catchment Model of Nitrogen (INCA-N) was applied to the River Lambourn, a Chalk river-system in southern England. The model's abilities to simulate the long-term trend and seasonal patterns in observed stream water nitrate concentrations from 1920 to 2003 were tested. This is the first time a semi-distributed, daily time-step model has been applied to simulate such a long time period and then used to calculate detailed catchment nutrient budgets which span the conversion of pasture to arable during the late 1930s and 1940s. Thus, this work goes beyond source apportionment and looks to demonstrate how such simulations can be used to assess the state of the catchment and develop an understanding of system behaviour. The mass-balance results from 1921, 1922, 1991, 2001 and 2002 are presented and those for 1991 are compared to other modelled and literature values of loads associated with nitrogen soil processes and export. The variations highlighted the problem of comparing modelled fluxes with point measurements but proved useful for identifying the most poorly understood inputs and processes thereby providing an assessment of input data and model structural uncertainty. The modelled terrestrial and instream mass-balances also highlight the importance of the hydrological conditions in pollutant transport. Between 1922 and 2002, increased inputs of nitrogen from fertiliser, livestock and deposition have altered the nitrogen balance with a shift from possible reduction in soil fertility but little environmental impact in 1922, to a situation of nitrogen accumulation in the soil, groundwater and instream biota in 2002. In 1922 and 2002 it was estimated that approximately 2 and 18 kg N ha(-1) yr(-1) respectively were exported from the land to the stream. The utility of the approach and further considerations for the best use of models are discussed. (C) 2008 Elsevier B.V. All rights reserved.

Quantifying sediment transport processes in periglacial mountain environments at a catchment scale using geomorphic process units

The research record on the quantification of sediment transport processes in periglacial mountain environments in Scandimvia dates back to the 1950s. A wide range of measurements is. available, especially from the Karkevagge region of northern Sweden. Within this paper satellite image analysis and tools provided by geographic information systems (GIS) are exploited in order to extend and improve this research and to complement geophysical methods. The processes of interest include mass movements such as solifluction, slope wash, dirty avalanches and rock-and boulder falls. Geomorphic process units have been derived in order to allow quantification via GIS techniques at a catchment scale. Mass movement rates based on existing Field measurements are employed in the budget calculation. In the Karkevagge catch ment. 80% of the area can be identified either as a source area for sediments or as a zone where sediments are deposited. The overall budget for the slopes beneath the rockwalls in the Karkevagge is approximately 680 t a(-1) whilst about 150 : a-1 are transported into the fluvial System.

Floodplain environmental change during the Younger Dryas and Holocene in Northwest Europe: Insights from the lower Kennet Valley, south central England

Many lowland rivers across northwest Europe exhibit broadly similar behavioural responses to glacial-interglacial transitions and landscape development. Difficulties exist in assessing these, largely because the evidence from many rivers remains limited and fragmentary. Here we address this issue in the context of the river Kennet, a tributary of the Thames, since c. 13,000 cal BP. Some similarities with other rivers are present, suggesting that regional climatic shifts are important controls. The Kennet differs from the regional pattern in a number of ways. The rate of response to sudden climatic change, particularly at the start of the Holocene and also mid-Holocene forest clearance, appears very high. This may reflect abrupt shifts between two catchment scale hydrological states arising from contemporary climates, land use change and geology. Stadial hydrology is dominated by nival regimes, with limited winter infiltration and high spring and summer runoff. Under an interglacial climate, infiltration is more significant. The probable absence of permafrost in the catchment means that a lag between the two states due to its gradual decay is unlikely. Palaeoecology, supported by radiocarbon dates, suggests that, at the very start of the Holocene, a dramatic episode of fine sediment deposition across most of the valley floor occurred, lasting 500-1000 years. A phase of peat accumulation followed as mineral sediment supply declined. A further shift led to tufa deposition, initially in small pools, then across the whole floodplain area, with the river flowing through channels cut in tufa and experiencing repeated avulsion. Major floods, leaving large gravel bars that still form positive relief features on the floodplain, followed mid-Holocene floodplain stability. Prehistoric deforestation is likely to be the cause of this flooding, inducing a major environmental shift with significantly increased surface runoff. Since the Bronze Age, predominantly fine sediments were deposited along the valley with apparently stable channels and vertical floodplain accretion associated with soil erosion and less catastrophic flooding. The Kennet demonstrates that, while a general pattern of river behaviour over time, within a region, may be identifiable, individual rivers are likely to diverge from this. Consequently, it is essential to understand catchment controls, particularly the relative significance of surface and subsurface hydrology. (c) 2005 Elsevier B.V. All rights reserved.

Modelling stream and soil water nitrate dynamics during experimentally increased nitrogen deposition in a coniferous forest catchment at Gardsjon, Sweden

Increased atmospheric deposition of inorganic nitrogen (N) may lead to increased leaching of nitrate (NO3-) to surface waters. The mechanisms responsible for, and controls on, this leaching are matters of debate. An experimental N addition has been conducted at Gardsjon, Sweden to determine the magnitude and identify the mechanisms of N leaching from forested catchments within the EU funded project NITREX. The ability of INCA-N, a simple process-based model of catchment N dynamics, to simulate catchment-scale inorganic N dynamics in soil and stream water during the course of the experimental addition is evaluated. Simulations were performed for 1990-2002. Experimental N addition began in 1991. INCA-N was able to successfully reproduce stream and soil water dynamics before and during the experiment. While INCA-N did not correctly simulate the lag between the start of N addition and NO 2 3 breakthrough, the model was able to simulate the state change resulting from increased N deposition. Sensitivity analysis showed that model behaviour was controlled primarily by parameters related to hydrology and vegetation dynamics and secondarily by in-soil processes.

Catchment-scale modelling of flow and nutrient transport in the Chalk unsaturated zone

The unsaturated zone exerts a major control on the delivery of nutrients to Chalk streams, yet flow and transport processes in this complex, dual-porosity medium have remained controversial. A major challenge arises in characterising these processes, both at the detailed mechanistic level and at an appropriate level for inclusion within catchment-scale models for nutrient management. The lowland catchment research (LOCAR) programme in the UK has provided a unique set of comprehensively instrumented groundwater-dominated catchments. Of these, the Pang and Lambourn, tributaries of the Thames near Reading, have been a particular focus for research into subsurface processes and surface water-groundwater interactions. Data from LOCAR and other sources, along with a new dual permeability numerical model of the Chalk, have been used to explore the relative roles of matrix and fracture flow within the unsaturated zone and resolve conflicting hypotheses of response. From the improved understanding gained through these explorations, a parsimonious conceptualisation of the general response of flow and transport within the Chalk unsaturated zone was formulated. This paper summarises the modelling and data findings of these explorations, and describes the integration of the new simplified unsaturated zone representation with a catchment-scale model of nutrients (INCA), resulting in a new model for catchment-scale flow and transport within Chalk systems: INCA-Chalk. This model is applied to the Lambourn, and results, including hindcast and forecast simulations, are presented. These clearly illustrate the decadal time-scales that need to be considered in the context of nutrient management and the EU Water Framework Directive. (C) 2007 Elsevier B.V. All rights reserved.

Fine sediment delivery and transfer in lowland catchments: modelling suspended sediment concentrations in response to hydrological forcing

Fine sediment delivery to and storage in stream channel reaches can disrupt aquatic habitats, impact river hydromorphology, and transfer adsorbed nutrients and pollutants from catchment slopes to the fluvial system. This paper presents a modelling toot for simulating the time-dependent response of the fine sediment system in catchments, using an integrated approach that incorporates both land phase and in-stream processes of sediment generation, storage and transfer. The performance of the model is demonstrated by applying it to simulate in-stream suspended sediment concentrations in two lowland catchments in southern England, the Enborne and the Lambourn, which exhibit contrasting hydrological and sediment responses due to differences in substrate permeability. The sediment model performs well in the Enborne catchment, where direct runoff events are frequent and peak suspended sediment concentrations can exceed 600 mg l(-1). The general trends in the in-stream concentrations in the Lambourn catchment are also reproduced by the model, although the observed concentrations are low (rarely exceeding 50 mg l(-1)) and the background variability in the concentrations is not fully characterized by the model. Direct runoff events are rare in this highly permeable catchment, resulting in a weak coupling between the sediment delivery system and the catchment hydrology. The generic performance of the model is also assessed using a generalized sensitivity analysis based on the parameter bounds identified in the catchment applications. Results indicate that the hydrological parameters contributing to the sediment response include those controlling (1) the partitioning of runoff between surface and soil zone flows and (2) the fractional loss of direct runoff volume prior to channel delivery. The principal sediment processes controlling model behaviour in the simulations are the transport capacity of direct runoff and the in-stream generation, storage and release of the fine sediment fraction. The in-stream processes appear to be important in maintaining the suspended sediment concentrations during low flows in the River Enborne and throughout much of the year in the River Lambourn. Copyright (c) 2007 John Wiley & Sons, Ltd.

Uncertainties in annual riverine phosphorus load estimation: Impact of load estimation methodology, sampling frequency, baseflow index and catchment population density

Models developed to identify the rates and origins of nutrient export from land to stream require an accurate assessment of the nutrient load present in the water body in order to calibrate model parameters and structure. These data are rarely available at a representative scale and in an appropriate chemical form except in research catchments. Observational errors associated with nutrient load estimates based on these data lead to a high degree of uncertainty in modelling and nutrient budgeting studies. Here, daily paired instantaneous P and flow data for 17 UK research catchments covering a total of 39 water years (WY) have been used to explore the nature and extent of the observational error associated with nutrient flux estimates based on partial fractions and infrequent sampling. The daily records were artificially decimated to create 7 stratified sampling records, 7 weekly records, and 30 monthly records from each WY and catchment. These were used to evaluate the impact of sampling frequency on load estimate uncertainty. The analysis underlines the high uncertainty of load estimates based on monthly data and individual P fractions rather than total P. Catchments with a high baseflow index and/or low population density were found to return a lower RMSE on load estimates when sampled infrequently than those with a tow baseflow index and high population density. Catchment size was not shown to be important, though a limitation of this study is that daily records may fail to capture the full range of P export behaviour in smaller catchments with flashy hydrographs, leading to an underestimate of uncertainty in Load estimates for such catchments. Further analysis of sub-daily records is needed to investigate this fully. Here, recommendations are given on load estimation methodologies for different catchment types sampled at different frequencies, and the ways in which this analysis can be used to identify observational error and uncertainty for model calibration and nutrient budgeting studies. (c) 2006 Elsevier B.V. All rights reserved.

Phosphorus dynamics and export in streams draining micro-catchments: Development of empirical models

Annual total phosphorus (TP) export data from 108 European micro-catchments were analyzed against descriptive catchment data on climate (runoff), soil types, catchment size, and land use. The best possible empirical model developed included runoff, proportion of agricultural land and catchment size as explanatory variables but with a low explanation of the variance in the dataset (R-2 = 0.37). Improved country specific empirical models could be developed in some cases. The best example was from Norway where an analysis of TP-export data from 12 predominantly agricultural micro-catchments revealed a relationship explaining 96% of the variance in TP-export. The explanatory variables were in this case soil-P status (P-AL), proportion of organic soil, and the export of suspended sediment. Another example is from Denmark where an empirical model was established for the basic annual average TP-export from 24 catchments with percentage sandy soils, percentage organic soils, runoff, and application of phosphorus in fertilizer and animal manure as explanatory variables (R-2 = 0.97).

Sensitivity analysis of a catchment-scale nitrogen model

There are now considerable expectations that semi-distributed models are useful tools for supporting catchment water quality management. However, insufficient attention has been given to evaluating the uncertainties inherent to this type of model, especially those associated with the spatial disaggregation of the catchment. The Integrated Nitrogen in Catchments model (INCA) is subjected to an extensive regionalised sensitivity analysis in application to the River Kennet, part of the groundwater-dominated upper Thames catchment, UK The main results are: (1) model output was generally insensitive to land-phase parameters, very sensitive to groundwater parameters, including initial conditions, and significantly sensitive to in-river parameters; (2) INCA was able to produce good fits simultaneously to the available flow, nitrate and ammonium in-river data sets; (3) representing parameters as heterogeneous over the catchment (206 calibrated parameters) rather than homogeneous (24 calibrated parameters) produced a significant improvement in fit to nitrate but no significant improvement to flow and caused a deterioration in ammonium performance; (4) the analysis indicated that calibrating the flow-related parameters first, then calibrating the remaining parameters (as opposed to calibrating all parameters together) was not a sensible strategy in this case; (5) even the parameters to which the model output was most sensitive suffered from high uncertainty due to spatial inconsistencies in the estimated optimum values, parameter equifinality and the sampling error associated with the calibration method; (6) soil and groundwater nutrient and flow data are needed to reduce. uncertainty in initial conditions, residence times and nitrogen transformation parameters, and long-term historic data are needed so that key responses to changes in land-use management can be assimilated. The results indicate the general, difficulty of reconciling the questions which catchment nutrient models are expected to answer with typically limited data sets and limited knowledge about suitable model structures. The results demonstrate the importance of analysing semi-distributed model uncertainties prior to model application, and illustrate the value and limitations of using Monte Carlo-based methods for doing so. (c) 2005 Elsevier B.V. All rights reserved.

Integrated Nitrogen CAtchment model (INCA) applied to a tropical catchment in the Atlantic Forest, Sao Paulo, Brazil

Stream-water flows and in-stream nitrate and ammonium concentrations in a small (36.7 ha) Atlantic Forest catchment were simulated using the Integrated Nitrogen in CAtchments (INCA) model version 1.9.4. The catchment, at Cunha, is in the Serra do Mar State Park, SE Brazil and is nearly pristine because the nearest major conurbations, Sao Paulo and Rio, are some 450 km distant. However, intensive farming may increase nitrogen (N) deposition and there are growing pressures for urbanisation. The mean-monthly discharges and NO3-N concentration dynamics were simulated adequately for the calibration and validation periods with (simulated) loss rates of 6.55 kg.ha(-1) yr(-1) for NO3-N and 3.85 kg.ha(-1) yr(-1) for NH4-N. To investigate the effects of elevated levels of N deposition in the future, various scenarios for atmospheric deposition were simulated; the highest value corresponded to that in a highly polluted area of Atlantic Forest in Sao Paulo City. It was found that doubling the atmospheric deposition generated a 25% increase in the N leaching rate, while at levels approaching the highly polluted Sao Paulo deposition rate, five times higher than the current rate, leaching increased by 240%, which would create highly eutrophic conditions, detrimental to downstream water quality. The results indicate that the INCA model can be useful for estimating N concentration and fluxes for different atmospheric deposition rates and hydrological conditions.

Hydraulic modeling of runoff over a rough surface under partial inundation

This study presents a numerical method to derive the Darcy- Weisbach friction coefficient for overland flow under partial inundation of surface roughness. To better account for the variable influence of roughness with varying levels of emergence, we model the flow over a network which evolves as the free surface rises. This network is constructed using a height numerical map, based on surface roughness data, and a discrete geometry skeletonization algorithm. By applying a hydraulic model to the flows through this network, local heads, velocities, and Froude and Reynolds numbers over the surface can be estimated. These quantities enable us to analyze the flow and ultimately to derive a bulk friction factor for flow over the entire surface which takes into account local variations in flow quantities. Results demonstrate that although the flow is laminar, head losses are chiefly inertial because of local flow disturbances. The results also emphasize that for conditions of partial inundation, flow resistance varies nonmonotonically but does generally increase with progressive roughness inundation.