Supraglacial weathering crust dynamics inferred from cryoconite hole hydrology

Water levels in cryoconite holes were monitored at high resolution over a 3‐week period on Austre Brøggerbreen (Svalbard). These data were combined with melt and energy balance modelling, providing insights into the evolution of the glacier's near‐surface hydrology and confirming that the hydrology of the near‐surface, porous ice known as the ‘weathering crust’ is dynamic and analogous to a shallow‐perched aquifer. A positive correlation between radiative forcing of melt and drainage efficiency was found within the weathering crust. This likely resulted from diurnal contraction and dilation of interstitial pore spaces driven by variations in radiative and turbulent fluxes in the surface energy balance, occasionally causing ‘sudden drainage events’. A linear decrease in water levels in cryoconite holes was also observed and attributed to cumulative increases in near‐surface ice porosity over the measurement period. The transport of particulate matter and microbes between cryoconite holes through the porous weathering crust is shown to be dependent upon weathering crust hydraulics and particle size. Cryoconite holes therefore yield an indication of the hydrological dynamics of the weathering crust and provide long‐term storage loci for cryoconite at the glacier surface. This study highlights the importance of the weathering crust as a crucial component of the hydrology, ecology and biogeochemistry of the glacier ecosystem and glacierized regions and demonstrates the utility of cryoconite holes as natural piezometers on glacier surfaces. Copyright © 2015 John Wiley & Sons, Ltd.     

Seasonal sediment fluxes forcing supraglacial melting on the Wright Lower Glacier,McMurdo Dry Valleys,Antarctica

Cold‐based glaciers exist in low temperature and low humidity environments in which shortwave radiation is the largest source of energy to the glacier surface and the energy budget is very sensitive to the surface albedo. Consequently, the presence of relatively low volumes of debris on glacier surfaces has a significant impact on the timing, magnitude and rate of ablation at the surface. The aim of this study is to understand how the presence of sediment on the glacier surface at the start of the melt season can affect meltwater generation and delivery on a cold‐based glacier. A combination of field measurements, energy balance modelling and chemical mixing modelling were used on the Wright Lower Glacier, McMurdo Dry Valleys, Antarctica, between October 2005 and January 2006 to address this aim. In this system, sediment was transported onto the glacier surface during the winter months (March–October) by foehn winds, which reduced surface albedo at the start of the summer melt season. The areas of the glacier on which sediment accumulated began to melt earlier than other parts of the glacier and experienced a longer melt season. Over the study period, the total ablation on the dirty surfaces was nine times greater than for clean ice. Ablation on the dirty surfaces is dominated by melting, whereas sublimation dominates the clean ice. As the sediment was unevenly distributed over the glacier surface, the variation in melt amount and timing drove the development of a cryoconite hole system. Copyright © 2012 John Wiley & Sons, Ltd.     

Using synchronous fluorescence spectroscopy and principal components analysis to monitor dissolved organic matter dynamics in a glacier system

The molecular characteristics of dissolved organic matter (DOM) reflect both its source material and its biogeochemical history. In glacial systems, DOM characteristics might be expected to change over the course of a melt season as changes in the glacier drainage system cause the mobilization of DOM from different OM pools. To test this hypothesis we used Principal Components Analysis (PCA) of synchronous fluorescence spectra to detect and describe changes in the DOM in meltwater from a glacier system in the Coast Mountains of northern British Columbia, Canada. For most of the melt season, the dominant component of subglacially routed meltwater DOM is characterized by a tyrosine‐like fluorophore. This DOM component is most likely derived from supraglacial snowmelt. During periods of high discharge, a second component of DOM is present which is humic in character and similar to DOM sampled from a nearby non‐glacial stream. This DOM component is inferred to be derived from a moss‐covered soil environment that has been glacially overrun. It is probably entrained into glacial melt waters when the supraglacial meltwater flux exceeds the capacity of the principal subglacial drainage channels and water floods areas of the glacier bed that are normally isolated from the subglacial drainage system. Another source of DOM also appears to be mobilized during periods of high air temperatures. It is characterized by both humic and proteinaceous fluorophores and may be derived from the drainage of supraglacial cryoconite holes. Copyright © 2009 John Wiley & Sons, Ltd.     

Factors controlling the temporal variability in dissolved oxygen regime of salmon spawning gravels

Dissolved oxygen (DO) plays a critical role in the development of the juvenile stages of benthic spawning fish and salmonids in particular. Factors influencing the DO regime within spawning gravels include the accumulation of fine sediment, penetration of groundwater or surface water into the gravels, thermal regime, and the consumption of oxygen by sediment and its associated organic fractions. In this field study, we quantify the DO regime within an artificial salmon redd at high temporal resolution. The environment within the redd is shown to be complex, with large variations in DO. Application of a numerical model (SIDO‐UK) enables for the first time, the quantification of the relative contributions to DO consumption from thermal regime, sediment accumulation and sediment oxygen demand. Sediment accumulation is shown to have a minor impact on DO in the redd whereas upwelling groundwater is identified as the most likely cause of the major changes in DO. Bed mobility has a minor impact on DO regime of the redd. The effects of fine sediment and oxygen supply on salmon embryo survival are estimated. Implications for river catchment management and prospects for future research are discussed. Copyright © 2012 John Wiley & Sons, Ltd.     

Range imaging: a new method for high‐resolution topographic measurements in small‐ and medium‐scale field sites

Topographic measurements are essential for the study of earth surface processes. Three‐dimensional data have been conventionally obtained through terrestrial laser scanning or photogrammetric methods. However, particularly in steep and rough terrain, high‐resolution field measurements remain challenging and often require new creative approaches. In this paper, range imaging is evaluated as an alternative method for obtaining surface data in such complex environments. Range imaging is an emerging time‐of‐flight technology, using phase shift measurements on a multi‐pixel sensor to generate a distance image of a surface. Its suitability for field measurements has yet not been tested. We found ambient light and surface reflectivity to be the main factors affecting error in distance measurements. Low‐reflectivity surfaces and strong illumination contrasts under direct exposure to sunlight lead to noisy distance measurements. However, regardless of lighting conditions, the accuracy of range imaging was markedly improved by averaging multiple images of the same scene. For medium ambient lighting (shade) and a light‐coloured surface the measurement uncertainty was approximately 9 mm. To further test the suitability of range imaging for field applications we measured a reach of a steep mountain stream with a horizontal resolution of approximately 1 cm (in the focal plane of the camera), allowing for the interpolation of a digital elevation model on a 2 cm grid. Comparison with an elevation model obtained from terrestrial laser scanning for the same site revealed that both models show similar degrees of topographic detail. Despite limitations in measurement range and accuracy, particularly at bright ambient lighting, range imaging offers three‐dimensional data in real time and video mode without the need of post‐processing. Therefore, range imaging is a useful complement or alternative to existing methods for high‐resolution measurements in small‐ to medium‐scale field sites. Copyright © 2012 John Wiley & Sons, Ltd.     

Factors affecting rill erosion of unpaved loess roads in China

This study explores factors that affect road surface erosion in a small watershed on the Loess Plateau. Global positioning system (GPS)‐assisted field surveys and geographical information system methods were applied. The results show that road surface rills in the watershed are more easily formed on main roads, which are disturbed by intensive human activities. Secondary unpaved road networks occupied the largest road surface area and contributed 49% of the total road surface rill volumes. Spatial analysis reveals that roads near residential areas or leading to other human‐disturbed land‐use types are at high risk of soil loss. In each road segment, slope gradient, road segment length and drainage area have impacts on surface rill formation and development. Among these factors, slope gradients have been verified as a controlling factor of rill erosion intensification. Both road segment length (R = 0.83, N = 82) and drainage area (R = 0.72 for road segment and 0.76 for upslope drainage areas, N = 82) significantly influence total road surface rill volumes. The interaction variable of road segment length multiplied by slope is more closely correlated with road segment soil loss than that of the independent variables alone. Linear equations composed of slope gradient, road segment length and upslope drainage area are proposed. The new equation performs much better at predicting surface soil loss from secondary road segments compared with the previous models, which have not considered upslope drainage areas. The relationships and equations from this study will be helpful for road erosion evaluation in a small watershed of the study area. Copyright © 2014 John Wiley & Sons, Ltd.     

Monitoring and modeling dissolved oxygen dynamics through continuous longitudinal sampling: a case study in Wen‐Rui Tang River,Wenzhou, China

Synoptic water sampling at a fixed site monitoring station provides only limited ‘snap‐shots’ of the complex water quality dynamics within a surface water system. However, water quality often changes rapidly in both spatial and temporal dimensions, especially in highly polluted urban rivers. In this study, we designed and applied a continuous longitudinal sampling technique to monitor the fine‐scale spatial changes of water quality conditions, assess water pollutant sources, and determine the assimilative capacity for biochemical oxygen demand (BOD) in an urban segment of the hypoxic Wen‐Rui Tang River in eastern China. The continuous longitudinal sampling was capable of collecting dissolved oxygen (DO) data every 5 s yielding a ~11 m sampling interval with a precision of ±0.1 mg L^?1. The Streeter and Phelps BOD‐DO model was used to calculate: (1) the oxygen consumption coefficient (K₁) required for calibration of water quality models, (2) BOD assimilative capacity, and (3) BOD source and load identification. In the 2014 m river segment sampled, the oxygen consumption coefficient (K₁) was 0.428 d^?1 (20°C), the total BOD discharge was 916 kg d^?1, and the BOD assimilative capacity was 382 kg d^?1 when the minimum DO level was set to 2 mg L^?1. In addition, the longitudinal analysis identified eight major drainage outlets (BOD point sources), which were verified by field observations. This new approach provides a simple, cost‐effective method of evaluating BOD‐DO dynamics over large spatial areas with rapidly changing water quality conditions, such as urban environments. It represents a major breakthrough in the development and application of water quality sampling techniques to obtain spatially distributed DO and BOD in real time. Copyright © 2012 John Wiley & Sons, Ltd.     

On the potential of MetOp ASCAT‐derived soil wetness indices as a new aperture for hydrological monitoring and prediction: a field evaluation over Luxembourg

In situ soil moisture data from the Bibeschbach experimental catchment in Luxembourg are used to evaluate relative surface soil moisture observed with the MetOp‐A Advanced Scatterometer (ASCAT). Filtered and bias‐corrected surface soil wetness indices (SWIs) derived from coarse‐resolution (25 km) C‐band scatterometer observations are shown to be highly correlated (r = 0.86) with catchment‐averaged soil moisture measured in the field. The combination of ASCAT and ENVISAT Advanced Synthetic Aperture Radar (ASAR) data sets yields high‐resolution (1 km) relative surface soil moisture that is equally well correlated with in situ measurements. It is concluded that for soil moisture monitoring applications at a catchment scale, the two soil moisture products are equivalent. The best correlation between the SWI derived from ASCAT and ASCAT‐ASAR with in situ soil moisture observations at ca. 5 cm depth is obtained with a characteristic time length parameter T equal to 288 h. These results suggest that satellite‐derived surface soil wetness may serve as proxy for soil storage that enables the monitoring of abrupt switches in river system dynamics to appear when an effective field capacity is exceeded and rapid subsurface stormflow is initiated. In catchments where soil moisture is the main controlling factor of rapid subsurface flow, MetOp ASCAT–derived SWI has the potential to monitor how a river system approaches a critical threshold. Copyright © 2011 John Wiley & Sons, Ltd.     

Two‐dimensional joint inversions of cross‐hole resistivity data and resolution analysis of combined arrays

In this study, a new two‐dimensional inversion algorithm was developed for the inversion of cross‐hole direct current resistivity measurements. In the last decades, various array optimisation methods were suggested for resistivity tomography. However, researchers have still collected data by using classical electrode arrays in most cross‐hole applications. Therefore, we investigated the accuracy of both the individual and the joint inversion of the classical cross‐hole arrays by using both synthetic and field data with the developed algorithm. We showed that the joint inversion of bipole–bipole, pole–bipole, bipole–pole, and pole–tripole electrode arrays gives inverse solutions that are closer to the real model than the individual inversions of the electrode array datasets for the synthetic data inversion. The model resolution matrix of the suggested arrays was used to analyse the inversion results. This model resolution analysis also showed the advantage of the joint inversion of bipole–bipole, pole–bipole, bipole–pole, and pole–tripole arrays. We also used sensitivity sections from each of the arrays and their superpositions to explain why joint inversion gives better resolution than the any individual inversion result.     

Developing channel and floodplain dimensions with limited data: a case study in the Amazon Basin

This research builds on the concept of hydraulic geometry and presents a methodology for estimating bankfull discharge and the hydraulic geometry coefficients and exponents for a station using limited data; only stage‐discharge and Landsat imagery. The approach is implemented using 82 streamflow gauging locations in the Amazon Basin. Using the estimated values for the hydraulic geometry relations, bankfull discharge, discharge data above bankfull and upstream drainage area at each site, relationships for estimating channel and floodplain characteristics as a function of drainage area are developed. Specifically, this research provides relationships for estimating bankfull discharge, bankfull depth, bankfull width, and floodplain width as a function of upstream drainage area in the Amazon Basin intended for providing reasonable cross‐section estimates for large scale hydraulic routing models. The derived relationships are also combined with a high resolution drainage network to develop relationships for estimating cumulative upstream channel lengths and surface areas as a function of the specified minimum channel width ranging from 2 m to 1 km (i.e. threshold drainage areas ranging from 1 to 431,000 km²). At the finest resolution (i.e. all channels greater than 2 m or a threshold area of 1 km²), the Amazon Basin contains approximately 4.4 million kilometers of channels with a combined surface area of 59,700 km². The intended use of these relationships is for partitioning total floodable area (channels versus lakes and floodplain lakes) obtained from remote sensing for biogeochemical applications (e.g. quantifying CO₂ evasion in the Amazon Basin). Copyright © 2011 John Wiley & Sons, Ltd.     

A network‐index‐based version of TOPMODEL for use with high‐resolution digital topographic data

This paper describes the preliminary development of a network‐index approach to modify and to extend the classic TOPMODEL. Application of the basic Beven and Kirkby form of TOPMODEL to high‐resolution (2·0 m) laser altimetric data (based upon the UK Environment Agency's light detection and ranging (LIDAR) system) to a 13·8 km² catchment in an upland environment identified many saturated areas that remained unconnected from the drainage network even during an extreme flood event. This is shown to be a particular problem with using high‐resolution topographic data, especially over large appreciable areas. To deal with the hydrological consequences of disconnected areas, we present a simple network index modification in which saturated areas are only considered to contribute when the topographic index indicates continuous saturation through the length of a flow path to the point where the path becomes a stream. This is combined with an enhanced method for dealing with the problem of pits and hollows, which is shown to become more acute with higher resolution topographic data. The paper concludes by noting the implications of the research as presented for both methodological and substantive research that is currently under way. Copyright © 2004 John Wiley & Sons, Ltd.     

High resolution mapping of supra‐glacial drainage pathways reveals link between micro‐channel drainage density,surface roughness and surface reflectance

This paper reports on the use of a small unmanned aerial vehicle (sUAV) carrying a standard compact camera, to construct a high resolution orthomosaic (OM) and digital elevation model (DEM) over the lower reaches of the glacier Midtre Lovénbreen, Svalbard. Structure from Motion (SfM) techniques were used to build the OM and DEM, and together these reveal insights into the nature of supra‐glacial drainage. Major meandering supra‐glacial drainage pathways show clear dynamism, via meander cutoffs and abandoned channels. In addition, the imagery reveals a very extensive network of smaller channels that may well carry substantial amounts of water. This network of channels is in part controlled by the structure of the glacier, but in turn, these channels have a significant impact on the ice surface. Roughness of the ice surface is higher where channels are most extensive. In addition, we find a relationship between channel density and surface reflectance, such that greater channel density is associated with lower reflectance values. Given the role of surface reflectance and roughness in the energy balance of glaciers, it is therefore apparent that extensive networks of small supra‐glacial channels across such glaciers have the potential to have an important impact on energy exchanges between the atmosphere and the ice surface. Copyright © 2015 John Wiley & Sons, Ltd.     

Hydrological influences on hyporheic water quality: implications for salmon egg survival

The spatial and temporal variability of groundwater–surface‐water (GW–SW) interactions was investigated in an intensively utilized salmon spawning riffle. Hydrochemical tracers, were used along with high‐resolution hydraulic head and temperature data to assess hyporheic dynamics. Surface and subsurface hydrochemistry were monitored at three locations where salmon spawning had been observed in previous years. Temperature and hydraulic head were monitored in three nests of three piezometers located to characterize the head, the run and the tail‐out of the riffle feature. Hydrochemical gradients between surface and subsurface water indicated increasing GW influence with depth into the hyporheic zone. Surface water was characterized by high dissolved oxygen (DO) concentrations, low alkalinity and conductivity. Hyporheic water was generally characterized by high levels of alkalinity and conductivity indicative of longer residence times, and low DO, indicative of reducing conditions. Hydrochemical and temperature gradients varied spatially over the riffle in response to changes in local GW–SW interactions at the depths investigated. Groundwater inputs dominated the head and tail of the riffle. The influence of SW increased in the area of accelerating flow and decreasing water depth through the run of the riffle. Temporal GW–SW interactions also varied in response to changing hydrological conditions. Gross changes in hyporheic hydrochemistry were observed at the weekly scale in response to changing flow conditions and surface water inputs to the hyporheic zone. During low flows, caused by freezing or dry weather, hyporheic hydrochemistry was dominated by GW inputs. During higher flows hyporheic hydrochemistry indicated that SW contributions increased. In addition, high‐resolution hydraulic head data indicated that rapid changes in GW–SW interactions occurred during hydrological events. The spatial, and possibly the temporal, variability of GW–SW interactions had a marked effect on the survival of salmon ova. It is concluded that hyporheic dynamics and their effect on stream ecology should be given increased consideration by fisheries and water resource managers. Copyright © 2004 John Wiley & Sons, Ltd.     

Multitemporal ALSM change detection,sediment delivery,and process mapping at an active earthflow

Remote mapping and measurement of surface processes at high spatial resolution is among the frontiers in Earth surface process research. Remote measurements that allow meter‐scale mapping of landforms and quantification of landscape change can revolutionize the study of landscape evolution on human timescales. At Mill Gulch in northern California, USA, an active earthflow was surveyed in 2003 and 2007 by airborne laser swath mapping (ALSM), enabling meter‐scale quantification of landscape change. We calculate four‐year volumetric flux from the earthflow and compare it to long‐term catchment average erosion rates from cosmogenic radionuclide inventories from adjacent watersheds. We also present detailed maps of changing features on the earthflow, from which we can derive velocity estimates and infer dominant process. These measurements rely on proper digital elevation model (DEM) generation and a simple surface‐matching technique to align the multitemporal data in a manner that eliminates systematic error in either dataset. The mean surface elevation of the earthflow and an opposite slope that was directly influenced by the earthflow decreased 14 ± 1 mm/yr from 2003 to 2007. By making the conservative assumption that these features were the dominant contributor of sediment flux from the entire Mill Gulch drainage basin during this time interval, we calculate a minimum catchment‐averaged erosion rate of 0·30 ± 0·02 mm/yr. Analysis of beryllium‐10 (¹⁰Be) concentrations in fluvial sand from nearby Russian Gulch and the South Fork Gualala River provide catchment averaged erosion rates of 0·21 ± 0·04 and 0·23 ± 0·03 mm/yr respectively. From translated landscape features, we can infer surface velocities ranging from 0·5 m/yr in the wide upper ‘source’ portion of the flow to 5 m/yr in the narrow middle ‘transport’ portion of the flow. This study re‐affirms the importance of mass wasting processes in the sediment budgets of uplifting weak lithologies. Copyright © 2011 John Wiley & Sons, Ltd.     

The estimation of shear-wave statics using in situ measurements in marine near-surface sediments

Shear‐wave statics in marine seismic exploration data are routinely too large to be estimated using conventional techniques. Near‐surface unconsolidated sediments are often characterized by low values of V_s and steep velocity gradients. Minor variations in sediment properties at these depths correspond to variations in the shear‐wave velocity and will produce significant static shifts. It is suggested that a significant proportion of the shear‐wave statics solution can be estimated by performing a separate high‐resolution survey to target near‐surface unconsolidated sediments. Love‐wave, shear‐wave refraction and geotechnical measurements were individually used to form high‐resolution near‐surface shear‐wave velocity models to estimate the shear‐wave statics for a designated survey line. Comparisons with predicted statics revealed that shear‐wave statics could not be estimated using a velocity model predicted by substituting geotechnical measurements into empirical relationships. Empirical relationships represent a vast simplification of the factors that control V_s and are therefore not sufficiently sensitive to estimate shear‐wave statics. Refraction measurements are potentially sensitive to short‐wavelength variations in sediment properties when combined with accurate navigational data. Statics estimated from Love‐wave data are less sensitive, and sometimes smoothed in appearance, since interpreted velocity values represent an average both laterally and vertically over the receiver array and the frequency–depth sensitivity range, respectively. For the survey site, statics estimated from near‐surface irregularities using shear‐wave refraction measurements represent almost half the total statics solution. More often, this proportion will be greater when bedrock relief is less.     

Bank erosion in agricultural drainage networks: new challenges from structure‐from‐motion photogrammetry for post‐event analysis

Drainage channels are an integral part of agricultural landscapes, and their impact on catchment hydrology is strongly recognized. In cultivated and urbanized floodplains, channels have always played a key role in flood protection, land reclamation, and irrigation. Bank erosion is a critical issue in channels. Neglecting this process, especially during flood events, can result in underestimation of the risk in flood‐prone areas. The main aim of this work is to consider a low‐cost methodology for the analysis of bank erosion in agricultural drainage networks, and in particular for the estimation of the volumes of eroded and deposited material. A case study located in the Veneto floodplain was selected. The research is based on high‐resolution topographic data obtained by an emerging low‐cost photogrammetric method (structure‐from‐motion or SfM), and results are compared to terrestrial laser scanning (TLS) data. For the SfM analysis, extensive photosets were obtained using two standalone reflex digital cameras and an iPhone5® built‐in camera. Three digital elevation models (DEMs) were extracted at the resolution of 0.1 m using SfM and were compared with the ones derived by TLS. Using the different DEMs, the eroded areas were then identified using a feature extraction technique based on the topographic parameter Roughness Index (RI). DEMs derived from SfM were effective for both detecting erosion areas and estimating quantitatively the deposition and erosion volumes. Our results underlined how smartphones with high‐resolution built‐in cameras can be competitive instruments for obtaining suitable data for topography analysis and Earth surface monitoring. This methodology could be potentially very useful for farmers and/or technicians for post‐event field surveys to support flood risk management. Copyright © 2015 John Wiley & Sons, Ltd.     

High‐frequency logging technologies reveal state‐dependent hyporheic process dynamics: implications for hydroecological studies

This paper presents a novel method for assessing hyporheic water quality dynamics using advances in sensor technology. High‐resolution (15 min) dissolved oxygen (DO) and hydraulic head data were combined to assess groundwater–surface water (GW–SW) interactions in the hyporheic zone. DO concentrations varied at fine temporal and spatial scales, depending on the relative contributions of GW and SW. The effect of sample frequency on observed patterns of variability was assessed with reference to studies of the ecology of salmon spawning habitat. Conventional approaches fail to capture the full range of temporal variability in hyporheic water quality and demonstrate the need to reassess the interpretations of previous studies of the hyporheic zone. © Crown Copyright 2006. Reproduced with the permission of Her Majesty's Stationery Office. Published by John Wiley & Sons, Ltd.     

Variability of unsaturated Bromide fluxes as measured through a layered volcanic vadose zone in New Zealand

The measured drainage fluxes through a layered volcanic vadose zone exhibited high spatial variability as a consequence of heterogeneous flow conditions. The drainage flux variability was quantified using automated equilibrium tension lysimeters, installed in close‐proximity and resulted in high variability in the Br masses recovered from a conservative tracer experiment. The primary cause of the heterogeneous flow was attributed to textural changes occurring at the interface between volcanic layers, resulting in development of funnel‐flow patterns, and further enhanced by the existence of hydrophobic conditions. The Br recoveries in individual automated equilibrium tension lysimeters were used to determine the corresponding variable sizes of the surface areas contributing drainage to the lysimeters. The tracer experiment confirmed the existence of unsaturated lateral transport occurring at the sloping interface of the coarse Taupo Ignimbrite material with the silty Palaeosol layer at approximately 4.2 m depth. This study demonstrates that measurements of both flux and solute concentrations at multiple locations are essential when heterogeneous flow is suspected to be present, to be able to determine reliable estimates of contaminant leaching through the vadose zone at the plot scale. Copyright © 2013 John Wiley & Sons, Ltd.     

Improvement of streams hydro‐geomorphological assessment using LiDAR DEMs

Hydro‐geomorphological assessments are an essential component for riverine management plans. They usually require costly and time‐consuming field surveys to characterize the spatial variability of key variables such as flow depth, width, discharge, water surface slope, grain size and unit stream power throughout the river corridor. The objective of this research is to develop automated tools for hydro‐geomorphological assessments using high‐resolution LiDAR digital elevation models (DEMs). More specifically, this paper aims at developing geographic information system (GIS) tools to extract channel slope, width and discharge from 1 m‐resolution LiDAR DEMs to estimate the spatial distribution of unit stream power in two contrasted watersheds in Quebec: a small agricultural stream (Des Fèves River) and a large gravel‐bed river (Matane River). For slope, the centreline extracted from the raw LiDAR DEM was resampled at a coarser resolution using the minimum elevation value. The channel width extraction algorithm progressively increased the centerline from the raw DEM until thresholds of elevation differences and slopes were reached. Based on the comparison with over 4000 differential global positioning system (GPS) measurements of the water surface collected in a 50 km reach of the Matane River, the longitudinal profile and slope estimates extracted from the raw and resampled LiDAR DEMs were in very good agreement with the field measurements (correlation coefficients ranging from 0 · 83 to 0 · 87) and can thus be used to compute stream power. The extracted width also corresponded very well to the channel as seen from ortho‐photos, although the presence of bars in the Matane River increased the level of error in width estimates. The estimated maximum unit stream power spatial patterns corresponded well with field evidence of bank erosion, indicating that LiDAR DEMs can be used with confidence for initial hydro‐geomorphological assessments. Copyright © 2013 John Wiley & Sons, Ltd.     

Structure,morphology and water flux of a subglacial drainage system,Midtdalsbreen, Norway

Digital elevation models of the surface and bed of Midtdalsbreen, Norway are used to calculate subglacial hydraulic potential and infer drainage system structure for a series of subglacial water pressure assumptions ranging from atmospheric to ice overburden. A distributed degree‐day model is used to calculate the spatial distribution of melt on the glacier surface throughout a typical summer, which is accumulated along the various drainage system structures to calculate water fluxes beneath the glacier and exiting the portals for the different water pressure assumptions. In addition, 78 dye‐tracing tests were performed from 33 injection sites and numerous measurements of water discharge were made on the main proglacial streams over several summer melt seasons. Comparison of the calculated drainage system structures and water fluxes with dye tracing results and measured proglacial stream discharges suggests that the temporally and spatially averaged steady‐state water pressures beneath the glacier are ~70% of ice overburden. Analysis of the dye return curves, together with the calculated subglacial water fluxes shows that the main drainage network on the eastern half of the glacier consists of a hydraulically efficient system of broad, low channels (average width/height ratio ≈ 75). The smaller drainage network on the west consists of a hydraulically inefficient distributed system, dominated by channels that are exceptionally broad and very low (average width/height ratio ≈ 350). The even smaller central drainage network also consists of a hydraulically inefficient distributed system, dominated by channels that are very broad and exceptionally low (average width/height ratio ≈ 450). The channels beneath the western and central glacier must be so broad and low that they can essentially be thought of as a linked cavity system. Copyright © 2011 John Wiley & Sons, Ltd.