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).

An application of the GLUE methodology for estimating the parameters of the INCA-N model

The conceptual and parameter uncertainty of the semi-distributed INCA-N (Integrated Nutrients in Catchments-Nitrogen) model was studied using the GLUE (Generalized Likelihood Uncertainty Estimation) methodology combined with quantitative experimental knowledge, the concept known as 'soft data'. Cumulative inorganic N leaching, annual plant N uptake and annual mineralization proved to be useful soft data to constrain the parameter space. The INCA-N model was able to simulate the seasonal and inter-annual variations in the stream-water nitrate concentrations, although the lowest concentrations during the growing season were not reproduced. This suggested that there were some retention processes or losses either in peatland/wetland areas or in the river which were not included in the INCA-N model. The results of the study suggested that soft data was a way to reduce parameter equifinality, and that the calibration and testing of distributed hydrological and nutrient leaching models should be based both on runoff and/or nutrient concentration data and the qualitative knowledge of experimentalist. (c) 2006 Elsevier B.V. All rights reserved.

Interactions between mass balance, atmospheric circulation, and recent climate change on the Djankuat Glacier, Caucasus Mountains, Russia

This paper reports recent changes in the mass balance record from the Djankuat Glacier, central greater Caucasus, Russia, and investigates possible relationships between the components of mass balance, local climate, and distant atmospheric forcing. The results clearly show that a strong warming signal has emerged in the central greater Caucasus, particularly since the 1993/1994 mass balance year, and this has led to a significant increase in the summer ablation of Djankuat. At the same time, there has been no compensating consistent increase in winter precipitation and accumulation leading to the strong net loss of mass and increase in glacier runoff. Interannual variability in ablation and accumulation is partly associated with certain major patterns of Northern Hemisphere climatic variability. The positive phase of the North Pacific (NP) teleconnection pattern forces negative geopotential height and temperature anomalies over the Caucasus in summer and results in reduced summer melt, such as in the early 1990s, when positive NP extremes resulted in a temporary decline in ablation rates. The positive phase of the NP is related to El Nino-Southern Oscillation, and it is possible that a teleconnection between the tropical Pacific sea surface temperatures and summer air temperatures in the Caucasus is bridged through the NP pattern. More recently, the NP pattern was predominantly negative, and this distant moderating forcing on summer ablation in the Caucasus was absent. Statistically significant correlations are observed between accumulation and the Scandinavian (SCA) teleconnection pattern. The frequent occurrence of the positive SCA phase at the beginning of accumulation season results in lower than average snowfall and reduced accumulation. The relationship between the North Atlantic Oscillation (NAO), Arctic Oscillation, and accumulation is weak, although positive precipitation anomalies in the winter months are associated with the negative phase of the NAO. A stronger positive correlation is observed between accumulation on Djankuat and geopotential height over the Bay of Biscay unrelated to the established modes of the Northern Hemisphere climatic variability. These results imply that the mass balance of Djankuat is sensitive to the natural variability in the climate system. Distant forcing, however, explains only 16% of the variance in the ablation record and cannot fully explain the recent increase in ablation and negative mass balance.

Evaluation of two hybrid metric-conceptual models, for simulating phosphorus transfer from agricultural land in the river enborne, a lowland UK catchment

Across Europe, elevated phosphorus (P) concentrations in lowland rivers have made them particularly susceptible to eutrophication. This is compounded in southern and central UK by increasing pressures on water resources, which may be further enhanced by the potential effects of climate change. The EU Water Framework Directive requires an integrated approach to water resources management at the catchment scale and highlights the need for modelling tools that can distinguish relative contributions from multiple nutrient sources and are consistent with the information content of the available data. Two such models are introduced and evaluated within a stochastic framework using daily flow and total phosphorus concentrations recorded in a clay catchment typical of many areas of the lowland UK. Both models disaggregate empirical annual load estimates, derived from land use data, as a function of surface/near surface runoff, generated using a simple conceptual rainfall-runoff model. Estimates of the daily load from agricultural land, together with those from baseflow and point sources, feed into an in-stream routing algorithm. The first model assumes constant concentrations in runoff via surface/near surface pathways and incorporates an additional P store in the river-bed sediments, depleted above a critical discharge, to explicitly simulate resuspension. The second model, which is simpler, simulates P concentrations as a function of surface/near surface runoff, thus emphasising the influence of non-point source loads during flow peaks and mixing of baseflow and point sources during low flows. The temporal consistency of parameter estimates and thus the suitability of each approach is assessed dynamically following a new approach based on Monte-Carlo analysis. (c) 2004 Elsevier B.V. All rights reserved.

Water quality and ecology of the River Lee: mass balance and a review of temporal and spatial data

A regional overview of the water quality and ecology of the River Lee catchment is presented. Specifically, data describing the chemical, microbiological and macrobiological water quality and fisheries communities have been analysed, based on a division into river, sewage treatment works, fish-farm, lake and industrial samples. Nutrient enrichment and the highest concentrations of metals and micro-organics were found in the urbanised, lower reaches of the Lee and in the Lee Navigation. Average annual concentrations of metals were generally within environmental quality standards although, oil many occasions, concentrations of cadmium, copper, lead, mercury and zinc were in excess of the standards. Various organic substances (used as herbicides, fungicides, insecticides, chlorination by-products and industrial solvents) were widely detected in the Lee system. Concentrations of ten micro-organic substances were observed in excess of their environmental quality standards, though not in terms of annual averages. Sewage treatment works were the principal point source input of nutrients. metals and micro-organic determinands to the catchment. Diffuse nitrogen sources contributed approximately 60% and 27% of the in-stream load in the upper and lower Lee respectively, whereas approximately 60% and 20% of the in-stream phosphorus load was derived from diffuse sources in the upper and lower Lee. For metals, the most significant source was the urban runoff from North London. In reaches less affected by effluent discharges, diffuse runoff from urban and agricultural areas dominated trends. Flig-h microbiological content, observed in the River Lee particularly in urbanised reaches, was far in excess of the EC Bathing Water Directive standards. Water quality issues and degraded habitat in the lower reaches of the Lee have led to impoverished aquatic fauna but, within the mid-catchment reaches and upper agricultural tributaries, less nutrient enrichment and channel alteration has permitted more diverse aquatic fauna.

Late-20th-century changes in glacier extent in the Caucasus Mountains, Russia/Georgia

Glaciers occupy an area of similar to 1600 km(2) in the Caucasus Mountains. There is widespread evidence of retreat since the Little Ice Age, but an up-to-date regional assessment of glacier change is lacking. In this paper, satellite imagery (Landsat Thematic Mapper and Enhanced Thematic Mapper Plus) is used to obtain the terminus position of 113 glaciers in the central Caucasus in 1985 and 2000, using a manual delineation process based on a false-colour composite (bands 5, 4, 3). Measurements reveal that 94% of the glaciers have retreated, 4% exhibited no overall change and 2% advanced. The mean retreat rate equates to similar to 8 m a(-1), and maximum retreat rates approach similar to 38 m a(-1). The largest (>10 km(2)) glaciers retreated twice as much (similar to 12 m a(-1)) as the smallest (<1 km(2)) glaciers (similar to 6 m a(-1)), and glaciers at lower elevations generally retreated greater distances. Supraglacial debris cover has increased in association with glacier retreat, and the surface area of bare ice has reduced by similar to 10% between 1985 and 2000. Results are compared to declassified Corona imagery from the 1960s and 1970s and detailed field measurements and mass-balance data for Djankuat glacier, central Caucasus. It is concluded that the decrease in glacier area appears to be primarily driven by increasing temperatures since the 1970s and especially since the mid-1990s. Continued retreat could lead to considerable changes in glacier runoff, with implications for regional water resources.

Modeling the Stream Water Nitrate Dynamics in 60,000-km(2) European Catchment, the Garonne, Southwest France

The spatial and temporal dynamics in the stream water NO3-N concentrations in a major European river-system, the Garonne (62,700 km(2)), are described and related to variations in climate, land management, and effluent point-sources using multivariate statistics. Building on this, the Hydrologiska Byrans Vattenbalansavdelning (HBV) rainfall-runoff model and the Integrated Catchment Model of Nitrogen (INCA-N) are applied to simulate the observed flow and N dynamics. This is done to help us to understand which factors and processes control the flow and N dynamics in different climate zones and to assess the relative inputs from diffuse and point sources across the catchment. This is the first application of the linked HBV and INCA-N models to a major European river system commensurate with the largest basins to be managed tinder the Water Framework Directive. The simulations suggest that in the lowlands, seasonal patterns in the stream water NO3-N concentrations emerge and are dominated by diffuse agricultural inputs, with an estimated 75% of the river load in the lowlands derived from arable farming. The results confirm earlier European catchment studies. Namely, current semi-distrubuted catchment-scale dynamic models, which integrate variations in land cover, climate, and a simple representation of the terrestrial and in-stream N cycle, are able to simulate seasonal NO3-N patterns at large spatial (> 300 km(2)) and temporal (>= monthly) scales using available national datasets.

The sediment trapping efficiency of the macro-tidal Daly Estuary, tropical Australia

Field studies were carried out on the water and sediment dynamics in the tropical, macro-tidal, Daly Estuary. The estuary is shallow, very-turbid, about 100 km long, and the entrance is funnel-shape. In the wet, high flow season, normal tidal ranges can be suppressed in the estuary, depending on inflow rates, and freshwater becomes dominant up to the mouth. At that time a fraction of the fine sediment load is exported offshore as a bottom-tagging nepheloid layer after the sediment falls out of suspension of the thin, near-surface, river plume. The remaining fraction and the riverine coarse sediment form a large sediment bar 10 km long, up to 6 m in height and extending across the whole width of the channel near the mouth. This bar, as well as shoals in the estuary, partially pond the mid- to upper-estuary. This bar builds up from the deposition of riverine sediment during a wet season with high runoff and can raise mean water level by up to 2 m in the upper estuary in the low flow season. This ponding effect takes about three successive dry years to disappear by the sediment forming the bar being redistributed all over the estuary by tidal pumping of fine and coarse sediment in the dry season, which is the low flow season. The swift reversal of the tidal currents from ebb to flood results in macro-turbulence that lasts about 20 min. Bed load transport is preferentially landward and occurs only for water currents greater than 0.6 m s(-1). This high value of the threshold velocity suggests that the sand may be cemented by the mud. The Daly Estuary thus is a leaky sediment trap with an efficiency varying both seasonally and inter-annually. (c) 2006 Elsevier Ltd. All rights reserved.

The impact of global freshwater forcing on the thermohaline circulation: adjustment of North Atlantic convection sites in a CGCM

On the time scale of a century, the Atlantic thermohaline circulation (THC) is sensitive to the global surface salinity distribution. The advection of salinity toward the deep convection sites of the North Atlantic is one of the driving mechanisms for the THC. There is both a northward and a southward contributions. The northward salinity advection (Nsa) is related to the evaporation in the subtropics, and contributes to increased salinity in the convection sites. The southward salinity advection (Ssa) is related to the Arctic freshwater forcing and tends on the contrary to diminish salinity in the convection sites. The THC changes results from a delicate balance between these opposing mechanisms. In this study we evaluate these two effects using the IPSL-CM4 ocean-atmosphere-sea-ice coupled model (used for IPCC AR4). Perturbation experiments have been integrated for 100 years under modern insolation and trace gases. River runoff and evaporation minus precipitation are successively set to zero for the ocean during the coupling procedure. This allows the effect of processes Nsa and Ssa to be estimated with their specific time scales. It is shown that the convection sites in the North Atlantic exhibit various sensitivities to these processes. The Labrador Sea exhibits a dominant sensitivity to local forcing and Ssa with a typical time scale of 10 years, whereas the Irminger Sea is mostly sensitive to Nsa with a 15 year time scale. The GIN Seas respond to both effects with a time scale of 10 years for Ssa and 20 years for Nsa. It is concluded that, in the IPSL-CM4, the global freshwater forcing damps the THC on centennial time scales.

In-street wind direction variability in the vicinity of a busy intersection in central London

We present results from fast-response wind measurements within and above a busy intersection between two street canyons (Marylebone Road and Gloucester Place) in Westminster, London taken as part of the DAPPLE (Dispersion of Air Pollution and Penetration into the Local Environment; www.dapple.org.uk) 2007 field campaign. The data reported here were collected using ultrasonic anemometers on the roof-top of a building adjacent to the intersection and at two heights on a pair of lamp-posts on opposite sides of the intersection. Site characteristics, data analysis and the variation of intersection flow with the above-roof wind direction (θref) are discussed. Evidence of both flow channelling and recirculation was identified within the canyon, only a few metres from the intersection for along-street and across-street roof-top winds respectively. Results also indicate that for oblique rooftop flows, the intersection flow is a complex combination of bifurcated channelled flows, recirculation and corner vortices. Asymmetries in local building geometry around the intersection and small changes in the background wind direction (changes in 15-min mean θref of 5–10 degrees) were also observed to have profound influences on the behaviour of intersection flow patterns. Consequently, short time-scale variability in the background flow direction can lead to highly scattered in-street mean flow angles masking the true multi-modal features of the flow and thus further complicating modelling challenges.

Global hydrology modelling and uncertainty: running multiple ensembles with a campus grid.

Uncertainties associated with the representation of various physical processes in global climate models (GCMs) mean that, when projections from GCMs are used in climate change impact studies, the uncertainty propagates through to the impact estimates. A complete treatment of this ‘climate model structural uncertainty’ is necessary so that decision-makers are presented with an uncertainty range around the impact estimates. This uncertainty is often underexplored owing to the human and computer processing time required to perform the numerous simulations. Here, we present a 189-member ensemble of global river runoff and water resource stress simulations that adequately address this uncertainty. Following several adaptations and modifications, the ensemble creation time has been reduced from 750 h on a typical single-processor personal computer to 9 h of high-throughput computing on the University of Reading Campus Grid. Here, we outline the changes that had to be made to the hydrological impacts model and to the Campus Grid, and present the main results. We show that, although there is considerable uncertainty in both the magnitude and the sign of regional runoff changes across different GCMs with climate change, there is much less uncertainty in runoff changes for regions that experience large runoff increases (e.g. the high northern latitudes and Central Asia) and large runoff decreases (e.g. the Mediterranean). Furthermore, there is consensus that the percentage of the global population at risk to water resource stress will increase with climate change.

Northern hemisphere snow: measurement, modelling and predictability

Northern hemisphere snow water equivalent (SWE) distribution from remote sensing (SSM/I), the ERA40 reanalysis product and the HadCM3 general circulation model are compared. Large differences are seen in the February climatologies, particularly over Siberia. The SSM/I retrieval algorithm may be overestimating SWE in this region, while comparison with independent runoff estimates suggest that HadCM3 is underestimating SWE. Treatment of snow grain size and vegetation parameterizations are concerns with the remotely sensed data. For this reason, ERA40 is used as `truth' for the following experiments. Despite the climatology differences, HadCM3 is able to reproduce the distribution of ERA40 SWE anomalies when assimilating ERA40 anomaly fields of temperature, sea level pressure, atmospheric winds and ocean temperature and salinity. However when forecasts are released from these assimilated initial states, the SWE anomaly distribution diverges rapidly from that of ERA40. No predictability is seen from one season to another. Strong links between European SWE distribution and the North Atlantic Oscillation (NAO) are seen, but forecasts of this index by the assimilation scheme are poor. Longer term relationships between SWE and the NAO, and SWE and the El Ni\~no-Southern Oscillation (ENSO) are also investigated in a multi-century run of HadCM3. SWE is impacted by ENSO in the Himalayas and North America, while the NAO affects SWE in North America and Europe. While significant connections with the NAO index were only present in DJF (and to an extent SON), the link between ENSO and February SWE distribution was seen to exist from the previous JJA ENSO index onwards. This represents a long lead time for SWE prediction for hydrological applications such as flood and wildfire forecasting. Further work is required to develop reliable large scale observation-based SWE datasets with which to test these model-derived connections.

An assessment of the fine sediment dynamics in an upland river system: INCA-Sed modifications and implications for fisheries

There is a need for better links between hydrology and ecology, specifically between landscapes and riverscapes to understand how processes and factors controlling the transport and storage of environmental pollution have affected or will affect the freshwater biota. Here we show how the INCA modelling framework, specifically INCA-Sed (the Integrated Catchments model for Sediments) can be used to link sediment delivery from the landscape to sediment changes in-stream. INCA-Sed is a dynamic, process-based, daily time step model. The first complete description of the equations used in the INCA-Sed software (version 1.9.11) is presented. This is followed by an application of INCA-Sed made to the River Lugg (1077 km2) in Wales. Excess suspended sediment can negatively affect salmonid health. The Lugg has a large and potentially threatened population of both Atlantic salmon (Salmo salar) and Brown Trout (Salmo trutta). With the exception of the extreme sediment transport processes, the model satisfactorily simulated both the hydrology and the sediment dynamics in the catchment. Model results indicate that diffuse soil loss is the most important sediment generation process in the catchment. In the River Lugg, the mean annual Guideline Standard for suspended sediment concentration, proposed by UKTAG, of 25 mg l− 1 is only slightly exceeded during the simulation period (1995–2000), indicating only minimal effect on the Atlantic salmon population. However, the daily time step simulation of INCA-Sed also allows the investigation of the critical spawning period. It shows that the sediment may have a significant negative effect on the fish population in years with high sediment runoff. It is proposed that the fine settled particles probably do not affect the salmonid egg incubation process, though suspended particles may damage the gills of fish and make the area unfavourable for spawning if the conditions do not improve.

Mitigation options for sediment and phosphorus loss from winter-sown arable crops

Sediment and P inputs to freshwaters from agriculture are a major problem in the United Kingdom (UK). This study investigated mitigation options for diffuse pollution losses from arable land. Field trials were undertaken at the hillslope scale over three winters at three UK sites with silt (Oxyaquic Hapludalf), sand (Udic Haplustept), and clay (Typic Haplaquept) soils. None of the mitigation treatments was effective in every year trialled, but each showed overall average reductions in losses. Over five site years, breaking up the compaction in tramlines (tractor wheel tracks) using a tine reduced losses of sediment and P to losses similar to those observed from areas without tramlines, with an average reduction in P loss of 1.06 kg TP ha(-1) Compared to traditional plowing, TP losses under minimum tillage were reduced by 0.30 kg TT ha(-1) over five site years, TP losses under contour cultivation were reduced by 0.30 kg TP ha(-1) over two site years, and TP losses using in-field barriers were reduced by 0.24 kg TP ha(-1) over two site years. In one site year, reductions in losses due to crop residue incorporation were nor significant. Each of the mitigation options trialled. is associated with a small cost at the farm-scale of up to 5 pound ha(-1), or with cost savings. The results indicate that each of the treatments his the potential to be a cost-effective mitigation option, but that tramline management is the most promising treatment, because tramlines dominate sediment and P transfer in surface runoff from arable hillslopes.

Influence of vegetation on the local climate and hydrology in the tropics: sensitivity to soil parameters

Land use change with accompanying major modifications to the vegetation cover is widespread in the tropics, due to increasing demands for agricultural land, and may have significant impacts on the climate. This study investigates (1) the influence of vegetation on the local climate in the tropics; (2) how that influence varies from region to region; and (3) how the sensitivity of the local climate to vegetation, and hence land use change, depends on the hydraulic characteristics of the soil. A series of idealised experiments with the Hadley Centre atmospheric model, HadAM3, are described in which the influence of vegetation in the tropics is assessed by comparing the results of integrations with and without tropical vegetation. The sensitivity of the results to the soil characteristics is then explored by repeating the experiments with a differing, but equally valid, description of soil hydraulic parameters. The results have shown that vegetation has a significant moderating effect on the climate throughout the tropics by cooling the surface through enhanced latent heat fluxes. The influence of vegetation is, however, seasonally dependent, with much greater impacts during the dry season when the availability of surface moisture is limited. Furthermore, there are significant regional variations both in terms of the magnitude of the cooling and in the response of the precipitation. Not all regions show a feedback of vegetation on the local precipitation; this result has been related both to vegetation type and to the prevailing meteorological conditions. An important finding has been the sensitivity of the results to the specification of the soil hydraulic parameters. The introduction of more freely draining soils has changed the soil-moisture contents of the control, vegetated system and has reduced, significantly, the climate sensitivity to vegetation and by implication, land use change. Changes to the soil parameters have also had an impact on the soil hydrology and its interaction with vegetation, by altering the partitioning between fast and slow runoff processes. These results raise important questions about the representation of highly heterogeneous soil characteristics in climate models, as well as the potential influence of land use change on the soil characteristics themselves.