Subglacial drainage system structure and morphology of Brewster Glacier,New Zealand

A global positioning system and ground penetrating radar surveys is used to produce digital elevation models of the surface and bed of Brewster Glacier. These are used to derive maps of subglacial hydraulic potential and drainage system structure using three different assumptions about the subglacial water pressure (P_w): (i) P_w = ice overburden; (ii) P_w = half ice overburden; (iii) P_w = atmospheric. Additionally, 16 dye‐tracing experiments at 12 locations were performed through a summer melt season. Dye return curve shape, together with calculations of transit velocity, dispersivity and storage, are used to infer the likely morphology of the subglacial drainage system. Taken together, the data indicate that the glacier is underlain by a channelised but hydraulically inefficient drainage system in the early summer in which water pressures are close to ice overburden. By mid‐summer, water pressures are closer to half‐ice overburden and the channelised drainage system is more hydraulically efficient. Surface streams that enter the glacier close to the location of major subglacial drainage pathways are routed quickly to the channels and then to the glacier snout. Streams that enter the glacier further away from the drainage pathways are routed slowly to the channels and then to the snout because they first flow through a distributed drainage system. Copyright © 2008 John Wiley & Sons, Ltd.     

Benchmarking observational uncertainties for hydrology: rainfall,river discharge and water quality

This review and commentary sets out the need for authoritative and concise information on the expected error distributions and magnitudes in observational data. We discuss the necessary components of a benchmark of dominant data uncertainties and the recent developments in hydrology which increase the need for such guidance. We initiate the creation of a catalogue of accessible information on characteristics of data uncertainty for the key hydrological variables of rainfall, river discharge and water quality (suspended solids, phosphorus and nitrogen). This includes demonstration of how uncertainties can be quantified, summarizing current knowledge and the standard quantitative results available. In particular, synthesis of results from multiple studies allows conclusions to be drawn on factors which control the magnitude of data uncertainty and hence improves provision of prior guidance on those uncertainties. Rainfall uncertainties were found to be driven by spatial scale, whereas river discharge uncertainty was dominated by flow condition and gauging method. Water quality variables presented a more complex picture with many component errors. For all variables, it was easy to find examples where relative error magnitudes exceeded 40%. We consider how data uncertainties impact on the interpretation of catchment dynamics, model regionalization and model evaluation. In closing the review, we make recommendations for future research priorities in quantifying data uncertainty and highlight the need for an improved ‘culture of engagement’ with observational uncertainties. Copyright © 2012 John Wiley & Sons, Ltd.     

Graphical tools based on Turc-Budyko plots to detect changes in catchment behaviour

Abstract

This paper aims at revisiting the use of non-dimensional representations for catchment change and model robustness analysis. As well as being helpful for catchment classification according to hydroclimatic conditions, the Turc-Budyko representation enables visualization of temporal changes in these conditions and/or modification of the hydrological behaviour of the catchment. This brings a new perspective to hydrological analysis, different from the classical time series plots, which helps when interpreting catchment functioning with respect to hydroclimatic constraints. These tools do not require statistical analyses of observed series and are therefore very simple to implement. Four case studies are considered here to illustrate the use of the Q/P = f(P/E_p) graph. When combined with the inter-annual Turc-Mezentsev formula, this visual framework enables anticipation of potential difficulties for models based on the same hypotheses, solely using the analysis of observed data.     

Growing season evapotranspiration in boreal fens in the Athabasca Oil Sands Region: Variability and environmental controls

Current efforts to assess changes to the wetland hydrology caused by growing anthropogenic pressures in the Athabasca Oil Sands Region (AOSR) require well-founded spatial and temporal estimates of actual evapotranspiration (ET), which is the dominant component of the water budget in this region. This study assessed growing season (May–September) and peak growing season (July) ET variability at a treed moderate-rich fen and treed poor fen (in 2013–2018), open poor fen (in 2011–2014), and saline fen (in 2015–2018) using eddy covariance technique and a set of complementary environmental data. Seasonal fluctuations in ET were positively related to net radiation, air temperature and vapour pressure deficit and followed trends typical for the Boreal Plains (BP) and AOSR with highest rates in June–July. However, no strong effect of water table position on ET was found. Strong surface control on ET is evident from lower ET values than potential evapotranspiration (PET); the lowest ET/PET was observed at saline fen, followed by open fen, moderately treed fen, and heavily treed fen, suggesting a strong influence of vegetation on water loss. In most years PET exceeded precipitation (P), and positive relations between P/PET and ET were observed with the highest July ET rates occurring under P/PET ~1. However, during months with P/PET > 1, increased P/PET was associated with decreased July ET. With respect to 30-year mean values of air temperature and P in the area, both dry and wet, cool and warm growing seasons (GS) were observed. No clear trends between ET values and GS wetness/coldness were found, but all wet GS were characterized by peak growing seasons with high daily ET variability.     

Snow interception modelling: Isolated observations have led to many land surface models lacking appropriate temperature sensitivities

When formulating a hydrologic model, scientists rely on parameterizations of multiple processes based on field data, but literature review suggests that more frequently people select parameterizations that were included in pre-existing models rather than re-evaluating the underlying field experiments. Problems arise when limited field data exist, when “trusted” approaches do not get reevaluated, and when sensitivities fundamentally change in different environments. The physics and dynamics of snow interception by conifers is just such a case, and it is critical to simulation of the water budget and surface albedo. The most commonly used interception parameterization is based on data from four trees from one site, but results from this field study are not directly transferable to locations with relatively warmer winters, where the dominant processes differ dramatically. Here, we combine a literature review with model experiments to demonstrate needed improvements. Our results show that the choice of model form and parameters can vary the fraction of snow lost through interception by as much as 30%. In most simulations, the warming of mean winter temperatures from −7 to 0°C reduces the modelled fraction of snow under the canopy compared to the open, but the magnitude of simulated decrease varies from about 10% to 40%. The range of results is even larger when considering models that neglect the melting of in-canopy snow in higher-humidity environments where canopy sublimation plays less of a role. Thus, we recommend that all models represent canopy snowmelt and include representation of increased loading due to increased adhesion and cohesion when temperatures rise from −3 to 0°C. In addition to model improvements, field experiments across climates and forest types are needed to investigate how to best model the combination of dynamically changing forest cover and snow cover to better understand and predict changes to albedo and water supplies.     

Using computational intelligence techniques to model subglacial water systems

Measurements of water pressure beneath Trapridge Glacier, Yukon Territory, Canada show that the basal water system is highly heterogeneous. Three types of behaviour were recorded: pressure records which are strongly correlated, records which are strongly anticorrelated, and records which alternate between strong correlation and strong anticorrelation. We take the pressure in bore-holes that are connected to the evacuation route for basal water as the forcing, and the other pressures as the response to this forcing. Previous work (Murray and Clarke 1995) has shown that these relationships can be modelled using low-order nonlinear differential equations optimized by inversion. However, despite optimizing the model parameters we cannot be sure that the final model forms are themselves optimal. Computational intelligence techniques provide alternative methods for fitting models and are robust to missing or noisy data, applicable to non-smooth models, and attempt to derive optimal model forms as well as optimal model parameters. Four computational intelligence techniques have been used and the results compared with the more conventional mathematical model. These methods were genetic programming, artificial neural networks, fuzzy logic and self-organizing maps. We compare each technique and offer an evaluation of their suitability for modelling the pressure data. The evaluation criteria are threefold: (1) goodness of fit and an ability to predict subsequent data under different surface weather conditions; (2) interpretability, and the extent and significance of any new insights offered into the physics of the glacier; (3) computation time. The results suggest that the suitability of the computational intelligence techniques to model these data increases with the complexity of the system to be modelled.     

A near‐annual palaeohydrological study based on testate amoebae from a sub‐alpine mire: surface wetness and the role of climate during the instrumental period

We present the first testate amoeba‐based palaeohydrological reconstruction from the Swiss Alps, and the first depth to the water table (DWT) calibration dataset for this region. Compared to existing models, our new calibration dataset performs well (RMSEP = 4.88), despite the length of the water table gradient covered (53 cm). The present‐day topography and vegetation of the study mire Mauntschas suggest that it is partly ombrotrophic (large Sphagnum fuscum hummocks, one of which was the coring site) but mostly under the minerotrophic influence of springs in the mire and runoff from the surrounding area. Ombrotrophic Sphagnum fuscum hummocks developed at the sampling site only during the last 50 years, when testate amoebae indicate a shift towards dry and/or acid conditions. Prior to AD 1950 the water table was much higher, suggesting that the influence of the mineral‐rich water prevented the development of ombrotrophic hummocks. The reconstructed DWT correlated with Pinus cembra pollen accumulation rates, suggesting that testate amoebae living on the mire and P. cembra growing outside of it partly respond to the same factor(s). Finally, temperature trends from the nearby meteorological station paralleled trends in reconstructed DWT. However, contrary to other studies made on raised bogs of northwestern Europe, the highest correlation was observed for winter temperature, despite the fact that testate amoebae would more logically respond to moisture conditions during the growing season. The observed correlation with winter temperature might reflect a control of winter severity on surface moisture during at least the first part of the growing season, through snow melt and soil frost phenomena influencing run‐off. More ecohydrological work on sub‐alpine mires is needed to understand the relationships between climate, testate amoebae and peatland development. Copyright © 2009 John Wiley & Sons, Ltd.     

Groundwater geochemistry of the Chihuahua City region in the Rio Conchos Basin (northern Mexico) and implications for water resources management

The Chihuahua City region, located in the semiarid-arid northern highlands of Mexico, has experienced intensive groundwater abstraction during the last 40 years to meet water demands in the region. A geochemical survey was carried out to investigate the evolution from baseline to modern conditions of a 130-km flow path including the El Sauz–Chihuahua–Aldama–San Diego de Alcalá regions. The research approach included the use of major chemical elements, chlorofluorocarbons and environmental isotope (¹⁸O, ²H, ¹³C and ¹⁴C) tracers. Stable isotopes indicate that groundwater evolves from the evaporation of local rainfall and surface water. Groundwater located at the lower end of the flow section is up to 6000 years old and older groundwater in the order of 9000 years BP was found in a deep well located in the upper part of the flow system, implying contribution from a neighbour basin. The background groundwater chemistry upstream of Chihuahua City results from feldspar weathering. Beyond Chihuahua City the chemical conditions are strongly modified owing to disposal of sewage from public and industrial water supplies into the Rio Chuviscar, subsequent allocation of this water to agricultural irrigation areas and direct infiltration under the river bed. As a consequence, anions like chloride and sulphate are mainly related to surface sources. Nitrate is controlled in part by sewage from public supply and industry and in part by agricultural practices. Arsenic and fluoride are related to weathering of rock formations of local mineralized ranges and subsequent enrichment of the basin-fill by magmatic processes. The results of this study have implications for groundwater management in an arid region that depends entirely on groundwater for domestic, industrial and agricultural water consumption. Copyright © 2008 John Wiley & Sons, Ltd.     

Improvement of distributed snowmelt energy balance modeling with MODIS‐based NDSI‐derived fractional snow‐covered area data

Describing the spatial variability of heterogeneous snowpacks at a watershed or mountain‐front scale is important for improvements in large‐scale snowmelt modelling. Snowmelt depletion curves, which relate fractional decreases in snow‐covered area (SCA) against normalized decreases in snow water equivalent (SWE), are a common approach to scale‐up snowmelt models. Unfortunately, the kinds of ground‐based observations that are used to develop depletion curves are expensive to gather and impractical for large areas. We describe an approach incorporating remotely sensed fractional SCA (FSCA) data with coinciding daily snowmelt SWE outputs during ablation to quantify the shape of a depletion curve. We joined melt estimates from the Utah Energy Balance Snow Accumulation and Melt Model (UEB) with FSCA data calculated from a normalized difference snow index snow algorithm using NASA's moderate resolution imaging spectroradiometer (MODIS) visible (0·545–0·565 µm) and shortwave infrared (1·628–1·652 µm) reflectance data. We tested the approach at three 500 m² study sites, one in central Idaho and the other two on the North Slope in the Alaskan arctic. The UEB‐MODIS‐derived depletion curves were evaluated against depletion curves derived from ground‐based snow surveys. Comparisons showed strong agreement between the independent estimates. Copyright © 2010 John Wiley & Sons, Ltd.     

The role of lithology,catchment size and the alluvial zone on the hydrogeochemistry of two intermittent Mediterranean streams

Hydrology and solute concentrations of two intermittent Mediterranean streams draining two nested catchments were compared. The two catchments were mainly underlain by granitic rocks and different types of sericitic schists. Only the lowland catchment had an alluvial zone and a well‐developed riparian forest. The rainfall–runoff relationship and the correlation between daily flow concentrations showed that hydrological behaviour was similar at both sites during most of the year. However, reverse fluxes were detected during the wetting and drying up periods only in the stream with an alluvial zone. The intermittence in stream flow also had effects on absolute solute concentrations, temporal solute dynamics and streamwater stoichiometry. Streamwater chemistry was not affected by drainage area, except for cations produced mainly by bedrock dissolution (i.e. calcium and magnesium) that increased with increasing catchment size. Differences in the relationship among cations and anions were detected between the two streams, which could be attributed to biogeochemical processes occurring in the alluvial zone. The multivariate model used in this study showed that stoichiometry was more useful than absolute concentrations when analyzing the influence of different lithologies on streamwater chemistry. Such differences were amplified in autumn, likely due to a low hydrological connectivity between the two nested catchments. Copyright © 2007 John Wiley & Sons, Ltd.