A simple model for predicting soil temperature in snow-covered and seasonally frozen soil: model description and testing

Microbial processes in soil are moisture, nutrient and temperature dependent and, consequently, accurate calculation of soil temperature is important for modelling nitrogen processes. Microbial activity in soil occurs even at sub-zero temperatures so that, in northern latitudes, a method to calculate soil temperature under snow cover and in frozen soils is required. This paper describes a new and simple model to calculate daily values for soil temperature at various depths in both frozen and unfrozen soils. The model requires four parameters average soil thermal conductivity, specific beat capacity of soil, specific heat capacity due to freezing and thawing and an empirical snow parameter. Precipitation, air temperature and snow depth (measured or calculated) are needed as input variables. The proposed model was applied to five sites in different parts of Finland representing different climates and soil types. Observed soil temperatures at depths of 20 and 50 cm (September 1981-August 1990) were used for model calibration. The calibrated model was then tested using observed soil temperatures from September 1990 to August 2001. R-2-values of the calibration period varied between 0.87 and 0.96 at a depth of 20 cm and between 0.78 and 0.97 at 50 cm. R-2 -values of the testing period were between 0.87 and 0.94 at a depth of 20cm. and between 0.80 and 0.98 at 50cm. Thus, despite the simplifications made, the model was able to simulate soil temperature at these study sites. This simple model simulates soil temperature well in the uppermost soil layers where most of the nitrogen processes occur. The small number of parameters required means, that the model is suitable for addition to catchment scale models.

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.

The Dubawnt Lake palaeo-ice stream: evidence for dynamic ics sheet behaviour on the Canadian Shield and insights regarding the controls on ice-stream location and vigour

We report evidence for a major ice stream that operated over the northwestern Canadian Shield in the Keewatin Sector of the Laurentide Ice Sheet during the last deglaciation 9000-8200 (uncalibrated) yr BP. It is reconstructed at 450 km in length, 140 km in width, and had an estimated catchment area of 190000 km. Mapping from satellite imagery reveals a suite of bedforms ('flow-set') characterized by a highly convergent onset zone, abrupt lateral margins, and where flow was presumed to have been fastest, a remarkably coherent pattern of mega-scale glacial lineations with lengths approaching 13 km and elongation ratios in excess of 40:1. Spatial variations in bedform elongation within the flow-set match the expected velocity field of a terrestrial ice stream. The flow pattern does not appear to be steered by topography and its location on the hard bedrock of the Canadian Shield is surprising. A soft sedimentary basin may have influenced ice-stream activity by lubricating the bed over the downstream crystalline bedrock, but it is unlikely that it operated over a pervasively deforming till layer. The location of the ice stream challenges the view that they only arise in deep bedrock troughs or over thick deposits of 'soft' fine-grained sediments. We speculate that fast ice flow may have been triggered when a steep ice sheet surface gradient with high driving stresses contacted a proglacial lake. An increase in velocity through calving could have propagated fast ice flow upstream (in the vicinity of the Keewatin Ice Divide) through a series of thermomechanical feedback mechanisms. It exerted a considerable impact on the Laurentide Ice Sheet, forcing the demise of one of the last major ice centres.

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.

Eutrophication control in river-systems: an application of INCA-P to the River Lugg

The Integrated Catchments model of Phosphorus dynamics (INCA-P) was applied to the River Lugg to determine the key factors controlling delivery of phosphorus to the main channel and to quantify the relative contribution of diffuse and point sources to the in-stream phosphorus (P) load under varying hydrological conditions. The model is able to simulate the seasonal variations and inter-annual variations in the in-stream total-phosphorus concentrations. However, difficulties in simulating diffuse inputs arise due to equifinality in the model structure and parameters. The River Lugg is split into upper and lower reaches. The upper reaches are dominated by grassland and woodland, so the patterns in the stream-water total-phosphorus concentrations are typical of diffuse source inputs; application of the model leads to estimates of the relative contribution to the in-stream P load from diffuse and point sources as 9:1. In the lower reaches, which are more intensively cultivated and urbanised, the stream-water total-phosphorus concentration dynamics are influenced more by point-sources; the simulated relative diffuse/point contribution to the in-stream P load is 1: 1. The model set-up and simulations are used to identify the key source-areas of P in the catchment, the P contribution of the Lugg to the River Wye during years with contrasting precipitation inputs, and the uptake and release of P from within-reach sediment. In addition, model scenarios are run to identify the impacts of likely P reductions at sewage treatment works on the in-stream soluble-reactive P concentrations and the suitability of this as a management option is assessed for reducing eutrophication.

Towards an improved understanding of the nitrate dynamics in lowland, permeable river-systems: Applications of INCA-N

The Integrated Catchment Model of Nitrogen (INCA-N) was applied to the Lambourn and Pang river-systems to integrate current process-knowledge and available-data to test two hypotheses and thereby determine the key factors and processes controlling the movement of nitrate at the catchment-scale in lowland, permeable river-systems: (i) that the in-stream nitrate concentrations were controlled by two end-members only: groundwater and soil-water, and (ii) that the groundwater was the key store of nitrate in these river-systems. Neither hypothesis was proved true or false. Due to equifinality in the model structure and parameters at least two alternative models provided viable explanations for the observed in-stream nitrate concentrations. One model demonstrated that the seasonal-pattern in the stream-water nitrate concentrations was controlled mainly by the mixing of ground- and soil-water inputs. An alternative model demonstrated that in-stream processes were important. It is hoped further measurements of nitrate concentrations made in the catchment soil- and ground-water and in-stream may constrain the model and help determine the correct structure, though other recent studies suggest that these data may serve only to highlight the heterogeneity of the system. Thus when making model-based assessments and forecasts it is recommend that all possible models are used, and the range of forecasts compared. In this study both models suggest that cereal production contributed approximately 50% the simulated in-stream nitrate toad in the two catchments, and the point-source contribution to the in-stream load was minimal. (c) 2006 Elsevier B.V. All rights reserved.

Nutrient monitoring, simulation and management within a major lowland UK river system: the Kennet

This paper describes an assessment of the nitrogen and phosphorus dynamics of the River Kennet in the south east of England. The Kennet catchment (1200 km(2)) is a predominantly groundwater fed river impacted by agricultural and sewage sources of nutrient (nitrogen and phosphorus) pollution. The results from a suite of simulation models are integrated to assess the key spatial and temporal variations in the nitrogen (N) and phosphorus (P) chemistry, and the influence of changes in phosphorous inputs from a Sewage Treatment Works on the macrophyte and epiphyte growth patterns. The models used are the Export Co-efficient model, the Integrated Nitrogen in Catchments model, and a new model of in-stream phosphorus and macrophyte dynamics: the 'Kennet' model. The paper concludes with a discussion on the present state of knowledge regarding the water quality functioning, future research needs regarding environmental modelling and the use of models as management tools for large, nutrient impacted riverine systems. (C) 2003 IMACS. Published by Elsevier B.V. All rights reserved.

Impacts of forestry on nitrogen in upland and lowland catchments: a comparison of the River Severn at Plynlimon in mid-Wales and the Bedford Ouse in south-east England using the INCA Model

The impacts of afforestation at Plynlimon in the Severn catchment, mid-Wales. and in the Bedford Ouse catchment in south-east England are evaluated using the INCA model to simulate Nitrogen (N) fluxes and concentrations. The INCA model represents the key hydrological and N processes operating in catchments and simulates the daily dynamic behaviour as well as the annual fluxes. INCA has been applied to five years of data front the Hafren and Hore headwater sub-catchments (6.8 km(2) area in total) of the River Severn at Plytilimon and the model was calibrated and validated against field data. Simulation of afforestation is achieved by altering the uptake rate parameters in the model. INCA simulates the daily N behaviour in the catchments with good accuracy as well as reconstructing the annual budgets for N release following clearfelling a four-fold increase in N fluxes was followed by a slow recovery after re-afforestation. For comparison, INCA has been applied to the large (8380 km(2)) Bedford Ouse catchment to investigate the impact of replacing 20% arable land with forestry. The reduction in fertiliser inputs from arable farming and the N uptake by the forest are predicted to reduce the N flux reaching the main river system, leading to a 33% reduction in N-Nitrate concentrations in the river water.

Impacts of climate change scenarios on dissolved oxygen in the River Thames, UK

A water quality model is used to assess the impact of possible climate change on dissolved oxygen (DO) in the Thames. The Thames catchment is densely populated and, typically, many pressures are anthropogenic. However, that same population also relies on the river for potable water supply and as a disposal route for treated wastewater. Thus, future water quality will be highly dependent on future activity. Dynamic and stochastic modelling has been used to assess the likely impacts on DO dynamics along the river system and the probability distributions associated with future variability. The modelling predictions indicate that warmer river temperatures and drought act to reduce dissolved oxygen concentrations in lowland river systems

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.

A baseline appraisal of water-dependant ecosystem services, the roles they play within desakota livelihood systems and their potential sensitivity to climate change

This report forms part of a larger research programme on 'Reinterpreting the Urban-Rural Continuum', which conceptualises and investigates current knowledge and research gaps concerning 'the role that ecosystems services play in the livelihoods of the poor in regions undergoing rapid change'. The report aims to conduct a baseline appraisal of water-dependant ecosystem services, the roles they play within desakota livelihood systems and their potential sensitivity to climate change. The appraisal is conducted at three spatial scales: global, regional (four consortia areas), and meso scale (case studies within the four regions). At all three scales of analysis water resources form the interweaving theme because water provides a vital provisioning service for people, supports all other ecosystem processes and because water resources are forecast to be severely affected under climate change scenarios. This report, combined with an Endnote library of over 1100 scientific papers, provides an annotated bibliography of water-dependant ecosystem services, the roles they play within desakota livelihood systems and their potential sensitivity to climate change. After an introductory, section, Section 2 of the report defines water-related ecosystem services and how these are affected by human activities. Current knowledge and research gaps are then explored in relation to global scale climate and related hydrological changes (e.g. floods, droughts, flow regimes) (section 3). The report then discusses the impacts of climate changes on the ESPA regions, emphasising potential responses of biomes to the combined effects of climate change and human activities (particularly land use and management), and how these effects coupled with water store and flow regime manipulation by humans may affect the functioning of catchments and their ecosystem services (section 4). Finally, at the meso-scale, case studies are presented from within the ESPA regions to illustrate the close coupling of human activities and catchment performance in the context of environmental change (section 5). At the end of each section, research needs are identified and justified. These research needs are then amalgamated in section 6.

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.

Calibrating flood inundation models: methods for assessing model performance

Satellite observed data for flood events have been used to calibrate and validate flood inundation models, providing valuable information on the spatial extent of the flood. Improvements in the resolution of this satellite imagery have enabled indirect remote sensing of water levels by using an underlying LiDAR DEM to extract the water surface elevation at the flood margin. Further to comparison of the spatial extent, this now allows for direct comparison between modelled and observed water surface elevations. Using a 12.5m ERS-1 image of a flood event in 2006 on the River Dee, North Wales, UK, both of these data types are extracted and each assessed for their value in the calibration of flood inundation models. A LiDAR guided snake algorithm is used to extract an outline of the flood from the satellite image. From the extracted outline a binary grid of wet / dry cells is created at the same resolution as the model, using this the spatial extent of the modelled and observed flood can be compared using a measure of fit between the two binary patterns of flooding. Water heights are extracted using points at intervals of approximately 100m along the extracted outline, and the students T-test is used to compare modelled and observed water surface elevations. A LISFLOOD-FP model of the catchment is set up using LiDAR topographic data resampled to the 12.5m resolution of the satellite image, and calibration of the friction parameter in the model is undertaken using each of the two approaches. Comparison between the two approaches highlights the sensitivity of the spatial measure of fit to uncertainty in the observed data and the potential drawbacks of using the spatial extent when parts of the flood are contained by the topography.