A joint analysis of river runoff and meteorological forcing in the Karakoram,upper Indus Basin

Satellite‐geodetic altimetry investigations in the Karakoram have indicated slight mass gain or loss of the glaciers during the early part of 21st century. Equivalent discharge in the upper Indus Basin due to these mass changes has been estimated at 5 to 10% of mean annual flow. However, satellite altimetry and geodetic glacier mass estimates in the extreme topography of the Karakoram have not yet been counter‐validated by hydrological analysis. Therefore, we present a first cross validation of three to five decades of river flow data from the three major watersheds in the Karakoram, with matching series of monthly precipitation, temperature, and evaporation provided by six atmospheric reanalysis products for 1979–2014. The analyses suggest that in most cases river flows have been increasing steadily from the end of the 1960s and 1970s to the middle of the 1990s and have stabilized or are in decline since then. Hunza watershed in Karakoram West shows consistently declining flows over the first half of the analysis period and stable flows during the second half for most of the summer melting season, suggesting mass accumulation. Rising river flows in the Shyok and Shigar watersheds, followed by stabilizing or slightly declining flows from 1995 onward, can be explained by consistently increasing precipitation during the first half of the analysis period, and successive stabilization or minor decline thereof. Flow data do not necessarily suggest considerable loss or gain of glacial mass in the Karakoram during the late 90s and early 2000s as suggested by satellite‐based altimetry studies.     

Glacial geomorphology and ice ages in Tibet and the surrounding mountains

Abstract   Since 1973 new data have been obtained on the maximum extent of glaciation in High Asia. Evidence for an ice sheet covering Tibet during the last glacial period means a radical rethinking about glaciation in the northern hemisphere. The ice sheet's subtropical latitude, vast size (2.4 million km²) and high elevation (6000 m a.s.l.) are supposed to have resulted in a substantial, albedo-induced cooling of the Earth's atmosphere and the disruption of summer monsoon circulation. Moraines were found to reach down to 460 m a.s.l. on the southern flank of the Himalayas and to 2300 m a.s.l. on the northern slope of the Tibetan Plateau, in the Qilian Shan region. On the northern slopes of the Karakoram, Aghil and Kuen-Lun Mountains, moraines occur as far down as 1900 m a.s.l. In southern Tibet, radiographic analyses of erratics suggest a former ice thickness of at least 1200 m. Glacial polish and roches moutonnées in the Himalayas and Karakoram suggest former glaciers as thick as 1200–2700 m. On the basis of this evidence, a 1100–1600-m lower equilibrium line has been reconstructed, resulting in an ice sheet of 2.4 million km², covering almost all of Tibet. Radiometric ages, obtained by different methods, classify this glaciation as isotope stage 3–2 in age (Würmian, the last glacial period, ca  60 000–18 000 years ago).     

Geomorphology under ice streams: Moving from form to process

Ice streams are integral components of an ice sheet's mass balance and directly impact on sea level. Their flow is governed by processes at the ice‐bed interface which create landforms that, in turn, modulate ice stream dynamics through their influence on bed topography and basal shear stresses. Thus, ice stream geomorphology is critical to understanding and modelling ice streams and ice sheet dynamics. This paper reviews developments in our understanding of ice stream geomorphology from a historical perspective, with a focus on the extent to which studies of modern and palaeo‐ice streams have converged to take us from a position of near‐complete ignorance to a detailed understanding of their bed morphology. During the 1970s and 1980s, our knowledge was limited and largely gleaned from geophysical investigations of modern ice stream beds in Antarctica. Very few palaeo‐ice streams had been identified with any confidence. During the 1990s, however, glacial geomorphologists began to recognise their distinctive geomorphology, which included distinct patterns of highly elongated mega‐scale glacial lineations, ice stream shear margin moraines, and major sedimentary depocentres. However, studying relict features could say little about the time‐scales over which this geomorphology evolved and under what glaciological conditions. This began to be addressed in the early 2000s, through continued efforts to scrutinise modern ice stream beds at higher resolution, but our current understanding of how landforms relate to processes remains subject to large uncertainties, particularly in relation to the mechanisms and time‐scales of sediment erosion, transport and deposition, and how these lead to the growth and decay of subglacial bedforms. This represents the next key challenge and will require even closer cooperation between glaciology, glacial geomorphology, sedimentology, and numerical modelling, together with more sophisticated methods to quantify and analyse the anticipated growth of geomorphological data from beneath active ice streams. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

帕隆藏布流域冰湖变化及危险性评估

冰湖作为区域气候变化的灵敏指示器和主要冰川灾害的启动器,认识其空间分布及变化特征对探讨冰湖对气候变化的响应规律及冰湖溃决危险性评估具有重要意义.基于1968-1980年地形图数据和1994-2016年Landsat TM/OLI遥感影像资料,综合利用RS、GIS技术和数理统计方法分析帕隆藏布流域面积≥ 0.01 km²冰湖时空分布及其动态变化,并对潜在危险性冰湖进行判别和评估.结果表明：2016年帕隆藏布流域共有冰湖351个,面积50.48 km²,且面积和数量分别以面积>1 km²和面积<0.1 km²的冰湖为主,这些冰湖主要分布于海拔2800~5400 m之间.近50年来帕隆藏布流域冰湖总体呈数量增多、面积增加态势;海拔<3000 m的冰湖相对稳定,而海拔>4500 m的冰湖数量和面积增加则相对迅速.近50年间帕隆藏布流域冰川面积减少591.34 km²,气候变暖导致的冰川末端退缩和冰川融水增加为冰湖形成和扩张提供了发育空间和物质来源.切毛措、光谢错等9个冰湖为潜在危险性冰湖,预计未来一段时间内帕隆藏布流域冰湖溃决可能处于活跃阶段,其形成和暴发也将更加频繁.     

Seafloor evidence for ice shelf flow across the Alaska–Beaufort margin of the Arctic Ocean

New imagery of ～14 100 km² of seafloor along a 640 km stretch of the Alaska and Beaufort margins (ABM) in water depths from 250 to 2800 m depicts a repetitive association of glaciogenic bedforms (lineations and iceberg scours), broad erosional bathymetric features and adjacent downslope turbidite gullies. These bedforms have styles, depths and orientations similar to features discovered earlier on the Chukchi Borderland, up to 800 km northwest of the ABM. Lineations occur across the surface of a flattened bathymetric bench interpreted to have formed by an ice shelf sliding along the continental slope and scraping the seafloor at temporary grounding locations. The glacial geology of surrounding areas suggests that an ice shelf probably flowed from the mouths of overdeepened glacial troughs in the Canadian Arctic Archipelago westward along the ABM and across the Chukchi Borderland. This curved pathway indicates an obstruction to ice flow in the central Canada Basin, possibly caused by either a basin‐wide ice shelf or by a pile‐up of mega‐bergs originating from the Eurasian side of the Arctic Ocean. The ice shelf that affected the ABM may have formed between Oxygen Isotopic Stage 4 to 5b, possibly correlating to an inferred intra‐Stage 5 widespread Beringian glaciation. Evidence for glaciogenic features on the ABM corroborates suggestions that large ice volumes and extents existed in the Arctic during Pleistocene glacial periods. These findings have far‐reaching implications for Arctic climate studies, ocean circulation, sediment stratigraphy and the stability of circum‐Arctic continental ice masses. Copyright © 2007 John Wiley & Sons, Ltd.     

Holocene calcrete crust deposits on the moraine of Batura Glacier, northern Pakistan

Abstract   Calcretes can be observed on the surface of old moraines around Batura Glacier in the upper Hunza Valley, Karakoram Mountains, Pakistan. They develop as a calcareous crust cementing small gravels under boulders. In order to understand the genesis of the calcrete crust, a variety of methods were employed: (i) study of mineralogy and geochemistry of a calcrete crust precipitated on the lateral moraine using X-ray diffractometer and electron probe microanalysis; (ii) analysis of solute chemistry of surface water and ice bodies around the Batura Glacier; and (iii) accelerator mass spectrometry ¹⁴C dating of the crust itself. The results indicate that the calcrete crust has definite laminated layers composed of a fine-grain and compact calcite layer, and a mineral fragment layer. The chemical composition of the calcite layer is approximately 60% CaO and 1% MgO. The mineral fragment layer consists of rounded grain materials up to 0.2 mm in diameter. It shows a graded bedding structure with fine grains of quartz, albite and muscovite. Meanwhile, as the Paleozoic Pasu limestone is distributed around the terminal of Batura Glacier, Ca cations dissolve in the melt water of the glacier. Accordingly, the calcrete crust is precipitated by decreases in CO₂ partial pressure from glacier ice and evaporation of the melt water, including high concentration of Ca²⁺ at ephemeral streams and small ponds stagnating between the moraine and glacial ice. On the basis of the AMS ¹⁴C age, the calcrete is considered to have formed approximately 8200 calibrated years bp under the Batura glacial stage.     

Contemporary glacigenic inputs to the dust cycle

The importance of glacigenic dust in the Earth's system during glacial periods is widely acknowledged. Under contemporary conditions, the world's largest dust sources are in low‐lying, hot, arid regions and this is where most aeolian research is focused. However the processes of dust production and emissions are still operating in cold climate regions, particularly in proglacial areas. This paper assesses current understanding of the relationship between glacierised landscapes and dust emissions and inputs to the global dust cycle. It focuses on how elements in the glacial and aeolian geomorphic sub‐systems interact to determine the magnitude, frequency and timing of aeolian dust emissions, and on feedback mechanisms between the systems. Where they have been measured, dust emission intensity and deposition rates in glacierised catchments are very high, in some cases far exceeding those in lower latitudes, however, few studies span long time scales. The impact of future glacier retreat on the balance between sediment supply, availability and aeolian transport capacity and implications for glacigenic dust emissions is also considered. This balance depends on relative spatial and temporal changes in meltwater suspended sediment concentration and wind strengths, which promote dust emissions, and patterns and rates of soil development and vegetation succession on recently‐deglaciated terrain which protect sediments from deflation. Retreat of the Antarctic ice sheet could mean that in future glacigenic contributions to the dust cycle exceed those of non‐glacigenic sources in the southern hemisphere. Copyright © 2012 John Wiley & Sons, Ltd.     

An automated algorithm for river ice monitoring over the Susquehanna River using the MODIS data

Reliable and prompt information on river ice condition and extent is needed to make accurate hydrological forecasts to predict ice jams breakups and issue timely flood warnings. This study presents a technique to detect and monitor river ice using observations from the MODIS instrument onboard the Terra satellite. The technique incorporates a threshold‐based decision tree image classification algorithm to process MODIS data and to determine the extent of ice. To differentiate between ice‐covered and ice‐free pixels within the riverbed, the algorithm combines observations in the visible and near‐infrared spectral bands. The developed technique presents the core of the MODIS‐based river ice mapping system, which has been developed to support National Oceanic and Atmospheric Administration NWS's operations. The system has been tested over the Susquehanna River in northeastern USA, where ice jam events leading to spring floods are a frequent occurrence. The automated algorithm generates three products: daily ice maps, weekly composite ice maps and running cloud‐free composite ice maps. The performance of the system was evaluated over nine winter seasons. The analysis of the derived products has revealed their good agreement with the aerial photography and with in situ observations‐based ice charts. The probability of ice detection determined from the comparison of the product with the high‐resolution Landsat imagery was equal to 91%. A consistent inverse relationship was found between the river discharge and the ice extent. The correlation between the discharge and the ice extent as determined from the weekly composite product reached 0.75. The developed CREST River Ice Observation System has been implemented at National Oceanic and Atmospheric Administration–Cooperative Remote Sensing Science and Technology Center as an operational Web tool allowing end users and forecasters to assess ice conditions on the river. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydro‐climatic impacts on the ice cover of the lower Peace River

Since the late 1960s, a paucity of ice‐jam flooding in the lower Peace River has resulted in prolonged dry periods and considerable reduction in the area covered by lakes and ponds that provide habitat for aquatic life in the Peace–Athabasca Delta (PAD) region. Though major ice jams occur at breakup, antecedent conditions play a significant role in their frequency and severity. These conditions are partly defined by the mode of freezeup and the maximum thickness that is attained during the winter, shortly before the onset of spring and development of positive net heat fluxes to the ice cover. Data from hydrometric gauge records and from field surveys are utilized herein to study these conditions. It is shown that freezeup flows are considerably larger at the present time than before regulation, and may be responsible for more frequent formation of porous accumulation covers. Despite a concomitant rise in winter temperatures, solid‐ice thickness has increased since the 1960s. Using a simple ice growth model, specifically developed for the study area, it is shown that porous accumulation covers enhance winter ice growth via accelerated freezing into the porous accumulation. Coupled with a reduction in winter snowfall, this effect can not only negate, but reverse, the effect of warmer winters on ice thickness, thus explaining present conditions. The present model is also shown to be a useful prediction tool, especially for extrapolating incomplete data to the end of the winter. Copyright © 2007 Crown in the right of Canada. Published by John Wiley & Sons, Ltd.     

Geomorphic evolution of small river–lake‐systems in northeast Germany during the Late Quaternary

This study investigates the post‐glacial development of four small river–lake systems in the Weichselian belt of northern central Europe. The valleys investigated are part of an immature drainage system characterized by frequent and abrupt changes in flow direction and the presence of numerous stagnant‐ice depressions in the valley course. The depressions contain thick sedimentary sequences which provide excellent archives for the reconstruction of the post‐glacial valley development. Study results indicate that the valleys reuse segments of former subglacial meltwater channels. During the Late Pleniglacial these channels carried meltwater streams. Stagnant‐ice melting occurred in stages from the Oldest Dryas to the early Holocene and was often followed by the formation of lakes in the valley course. Flow reversals occurred during the Late‐glacial–Holocene transition and were in response to general base‐level lowering caused by stagnant‐ice melting, headwater erosion and lake overspills. Lacustrine deposition typically started during the early Late‐glacial comprising mainly silicate gyttjas, whereas organic gyttjas and peats accumulated during the Allerød. The Younger Dryas is associated with a marked increase in fluvial and aeolian sedimentation, and lake‐level high stands. This was followed by early Holocene lake‐level low stands and a subsequent stabilization phase with decreasing silicate input and increasing organic lacustrine deposition. In general, dramatic changes in Late Pleniglacial to early Holocene sedimentation suggest that small‐scale catastrophic events played a more important role in triggering geomorphic changes then previously recognized. Infilling continued until peat accumulation and terrestrialization of lake basins became widespread during the mid‐ to late Holocene. Beginning in the late Holocene anthropogenic influences become important mainly involving an increase in sediment supply due to forest clearing and land use, followed by mill stowage, river course correction and anthropogenic lake‐level manipulations. Copyright © 2007 John Wiley & Sons, Ltd.     

Geothermal control on flow patterns in the Last Glacial Maximum ice sheet of Iceland

Because it is located both on the Mid‐Atlantic Ridge and on a mantle plume, Iceland is a region of intense tectonics and volcanism. During the last glaciation, the island was covered by an ice sheet approximately 1000 m thick. A reconstruction of the ice flow lines, based on glacial directional features, shows that the ice sheet was partly drained through fast‐flowing streams. Fast flow of the ice streams has been recorded in megascale lineations and flutes visible on the currently deglaciated bedrock, and is confirmed by simple mass balance considerations. Locations of the major drainage routes correlate with locations of geothermal anomalies, suggesting that ice stream activity was favoured by lubrication of the bed by meltwater produced in regions of high geothermal heat flux. Similar control of ice flow by geothermal activity is expected in ice sheets currently covering tectonically and volcanically active area such as the West Antarctic ice sheet. Copyright © 2000 John Wiley & Sons, Ltd.     

Size,shape and spatial arrangement of mega‐scale glacial lineations from a large and diverse dataset

Mega‐scale glacial lineations (MSGLs) are a characteristic landform on ice stream beds. Solving the puzzle of their formation is key to understanding how ice interacts with its bed and how this, in turn, influences the dynamics of ice streams. However, a comprehensive and detailed characterization of this landform's size, shape and spatial arrangement, which might serve to test and refine formational theories, is largely lacking. This paper presents a detailed morphometric analysis and comparison of 4043 MSGLs from eight palaeo‐ice stream settings: three offshore (Norway and Antarctica), four onshore (Canada), and one from under a modern ice stream in West Antarctica. The length of MSGLs is lower than previously suggested (mode 1000–2000 m; median 2892 m), and they initiate and terminate at various locations on an ice stream bed. Their spatial arrangement reveals a pattern that is characterized by an exceptional parallel conformity (80% of all mapped MSGLs have an azimuth within 5° from the mean values), and a fairly constant lateral spacing (mode 200–300 m; median 330 m), which we interpret as an indication that MSGLs are a spatially self‐organized phenomenon. Results show that size, shape and spatial arrangement of MSGLs are consistent both within and also generally between different ice stream beds. We suggest this results from a common mechanism of formation, which is largely insensitive to local factors. Although the elongation of MSGLs (mode 6–8; median 12.2) is typically higher than features described as drumlins, these values and those of their width (mode 100–200 m; median 268 m) overlap, which suggests the two landforms are part of a morphological continuum and may share a similar origin. We compare their morphometry with explicit predictions made by the groove‐ploughing and rilling instability theories of MSGL formation. Although the latter was most compatible, neither is fully supported by observations. © 2014 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Amplification of hydrological model uncertainties in projected climate simulations of the Upper Indus Basin: Does it matter where the water is coming from?

There is a growing appreciation of the uncertainties in the estimation of snow-melt and glacier-melt as a result of climate change in high elevation catchments. Through a detailed examination of three hydrological models in two catchments, and interpretation of results from previous studies, we observed that many variations in estimated streamflow could be explained by the selection of a best parameter set from the possible good model parameters. The importance of understanding changing glacial dynamics is critically important for our study areas in the Upper Indus Basin where Pakistan's policymakers are planning infrastructure to meet the future energy and water needs of hundreds of millions of people downstream. Yet, the effect of climate on glacial runoff and climate on snowmelt runoff is poorly understood. With the HBV model, for example, we estimated glacial melt as between 56% and 89% for the Hunza catchment. When rainfall was a scaled parameter, the models estimated glacial melt as between 20% and 100% of streamflow. These parameter sets produced wildly different projections of future climate for RCP8.5 scenarios in 2046–2075 compared to 1976–2005. Assuming no glacial shrinkage, for one climate projection, we found that the choice among good parameter sets resulted in projected values of future streamflow across a range from +54% to +125%. Parameter selection was the most significant source of uncertainty in the glaciated catchment and amplified climate model uncertainty, whereas climate model choice was more important in the rainfall dominated catchment. Although the study focuses on Pakistan, the overall conclusions are instructive for other similar regions in the world. We suggest that modellers of glaciated catchments should present results from at least the book-ends: models with low sensitivity to ice-melt and models with high sensitivity to ice-melt. This would reduce confusion among decision makers when they are faced with similar contrasting results.     

Internal structure and sources of selected ice jams on the lower Vistula River

The lower stretch of the Vistula is the most ice‐jammed river section on the North European Plain. Since 1982, the structure of hanging dams has been studied by means of a mechanical non‐core sampler. In this article, a selected of field research results of the hanging dams' structure and the degree of filling of the cross section with ice obtained during surface ice‐jam events in the years 1995–2014 are presented, along with an explanation of their causes. Surface ice jams occurred during spring snowmelt surges and ice breakups and also during freeze‐up and ice‐covered periods. Their main cause was changes in the river flow and was also those affected by anthropogenic sources. A characteristic feature of the analysed cross sections was the considerable share of the underhanging ice dam's firm accumulation with ice floes, when the cross section would be filled with ice in excess of 70%. In most cases, due to low river discharge, there was no substantial flooding damage. Copyright © 2016 John Wiley & Sons, Ltd.     

Wavelet analysis of inter‐annual variability in the runoff regimes of glacial and nival stream catchments,Bow Lake,Alberta

Continuous wavelet analyses of hourly time series of air temperature, stream discharge, and precipitation are used to compare the seasonal and inter‐annual variability in hydrological regimes of the two principal streams feeding Bow Lake, Banff National Park, Alberta: the glacial stream draining the Wapta Icefields, and the snowmelt‐fed Bow River. The goal is to understand how water sources and flow routing differ between the two catchments. Wavelet spectra and cross‐wavelet spectra were determined for air temperature and discharge from the two streams for summers (June–September) 1997–2000, and for rainfall and discharge for the summers of 1999 and 2000. The diurnal signal of the glacial runoff was orders of magnitude higher in 1998 than in other years, indicating that significant ice exposure and the development of channelized glacial drainage occurred as a result of the 1997–98 El Niño conditions. Early retreat of the snowpack in 1997 and 1998 led to a significant summer‐long input of melt runoff from a small area of ice cover in the Bow River catchment; but such inputs were not apparent in 1999 and 2000, when snow cover was more extensive. Rainfall had a stronger influence on runoff and followed quicker flow paths in the Bow River catchment than in the glacial catchment. Snowpack thickness and catchment size were the primary controls on the phase relationship between temperature and discharge at diurnal time scales. Wavelet analysis is a fast and effective means to characterize runoff, temperature, and precipitation regimes and their interrelationships and inter‐annual variability. The technique is effective at identifying inter‐annual and seasonal changes in the relative contributions of different water sources to runoff, and changes in the time required for routing of diurnal meltwater pulses through a catchment. However, it is less effective at identifying changes/differences in the type of the flow routing (e.g. overland flow versus through flow) between or within catchments. Copyright © 2003 John Wiley & Sons, Ltd.     

Hydrograph separation and precipitation source identification using stable water isotopes and conductivity: River Ganga at Himalayan foothills

The observed retreat of several Himalayan glaciers and snow packs is a cause of concern for the huge population in southern Asia that is dependent on the glacial‐fed rivers emanating from Himalayas. There is considerable uncertainty about how cryospheric recession in the Himalayan region will respond to climate change, and how the water resource availability will be affected. As a first step towards quantifying the contribution of glacier‐melt water, hydrograph separation of River Ganga at Rishikesh into its constituent components, namely (i) surface runoff, (ii) glacial ice‐melt and (iii) groundwater discharge has been done in this paper. A three‐component mixing model has been employed using the values of δ¹⁸O and electrical conductivity (EC) of the river water, and its constituents, to estimate the time‐varying relative fraction of each component. The relative fraction of the surface runoff peaks (70–90%) during winter, due to the near‐zero contribution of glacial ice‐melt, essentially represents the melting of surface snow from the catchment. The contribution of glacial ice‐melt to the stream discharge peaks during summer and monsoon reaches a maximum value of ～40% with an average of 32%. The fraction of groundwater discharge varies within a narrow range (15 ± 5%) throughout the year. On the basis of the variation in the d‐excess values of river water, it is also suggested that the snow‐melt and ice‐melt component has a significant fraction derived from winter precipitation with moisture source from mid‐latitude westerlies (also known as western disturbances). Copyright © 2010 John Wiley & Sons, Ltd.     

Tree ring analysis of hydro‐climatic thresholds that trigger ice jams on the Mistassini River,Quebec

This study documents the spatiotemporal variations in the frequency and magnitude of ice jams in the Mistassini River and applies that information to the identification of the hydro‐climatic threshold conditions associated with major events. Ice jams cause severe upheavals in water flow, which result in flooding upstream of the ice jam front, and therefore represent a significant geohazard to riparian populations. To analyze the spatiotemporal variations in the magnitude–frequency of ice jams, the Mistassini River was first divided into six different sites representing different geomorphological contexts. A 50‐year ice jam chronology was constructed from 85 damaged trees from all of the study sites. This chronology was then coupled with hydro‐climatic variables to construct classification trees, which helped identify the conditions and hydro‐climatic thresholds favourable to the triggering of ice jams in a predictive model. The results indicate complex interactions between the characteristics of flow, ice cover and river morphology that affect the frequency and magnitude of ice jam events on the Mistassini River. These factors affect the frequency and magnitude of ice jam events. The triggering of extreme events seems particularly influenced by exceptional ice conditions and sites with high sinuosity and islands. Copyright © 2015 John Wiley & Sons, Ltd.     

A multi‐dimensional analysis of pro‐glacial landscape change at Sólheimajökull,southern Iceland

Pro‐glacial landscapes are some of the most active on Earth. Previous studies of pro‐glacial landscape change have often been restricted to considering either sedimentological, geomorphological or topographic parameters in isolation and are often mono‐dimensional. This study utilized field surveys and digital elevation model (DEM) analyses to quantify planform, elevation and volumetric pro‐glacial landscape change at Sólheimajökull in southern Iceland for multiple time periods spanning from 1960 to 2010. As expected, the most intense geomorphological changes persistently occurred in the ice‐proximal area. During 1960 to 1996 the pro‐glacial river was relatively stable. However, after 2001 braiding intensity was higher, channel slope shallower and there was a shift from overall incision to aggradation. Attributing these pro‐glacial river channel changes to the 1999 jökulhlaup is ambiguous because it coincided with a switch from a period of glacier advance to that of glacier retreat. Furthermore, glacier retreat (of ~40 m yr^?1) coincided with ice‐marginal lake development and these two factors have both altered the pro‐glacial river channel head elevation. From 2001 to 2010 progressive increase in channel braiding and progressive downstream incision occurred; these together probably reflecting stream power due to increased glacier ablation and reduced sediment supply due to trapping of sediment by the developing ice‐marginal lake. Overall, this study highlights rapid spatiotemporal pro‐glacial landscape reactions to changes in glacial meltwater runoff regimes, glacier terminus position, sediment supply and episodic events such as jökuhlaups. Recognizing the interplay of these controlling factors on pro‐glacial landscapes will be important for understanding the geological record and for landscape stability assessments. Copyright © 2014 John Wiley & Sons, Ltd.     

Ice‐jam flood risk assessment and mapping

In northern regions, river ice‐ jam flooding can be more severe than open‐water flooding causing property and infrastructure damages, loss of human life and adverse impacts on aquatic ecosystems. Very little has been performed to assess the risk induced by ice‐related floods because most risk assessments are limited to open‐water floods. The specific objective of this study is to incorporate ice‐jam numerical modelling tools (e.g. RIVICE, Monte‐Carlo simulation) into flood hazard and risk assessment along the Peace River at the Town of Peace River (TPR) in Alberta, Canada. Adequate historical data for different ice‐jam and open‐water flooding events were available for this study site and were useful in developing ice‐affected stage‐frequency curves. These curves were then applied to calibrate a numerical hydraulic model, which simulated different ice jams and flood scenarios along the Peace River at the TPR. A Monte‐Carlo analysis was then carried out to acquire an ensemble of water level profiles to determine the 1 : 100‐year and 1 : 200‐year annual exceedance probability flood stages for the TPR. These flood stages were then used to map flood hazard and vulnerability of the TPR. Finally, the flood risk for a 200‐year return period was calculated to be an average of $32/m²/a ($/m²/a corresponds to a unit of annual expected damages or risk). Copyright © 2016 John Wiley & Sons, Ltd.     

The source and fate of sediment and mercury in Hunza River basin,Northern Areas,Pakistan

This article presents results of mercury in surface waters from Hunza River basin, Northern Areas, Pakistan. Small‐scale gold mining activities along the Hunza and Gilgit rivers are long known to be discharging mercury in the amalgamation and roasting processes. Previous studies reported high mercury concentrations in soils close to mining operations as well as serious health problems for miners. However, none of the studies have focused on the level of contamination in aqueous environments. This is the first study on the investigation of source and fate of sediment and river‐borne mercury in the Hunza River. The samples collected near gold panning sites showed higher mercury concentrations than critical levels established by the U.S. Environmental Protection Agency. The observed dissolved mercury concentrations ranged from 5.10 to 25.25 ng/l, whereas particulate‐bound mercury ranged from 4.85 to 154.62 ng/l. Particulate‐phase mercury corresponded to more than 75% of the total observed mercury concentrations for all of the sampled rivers. Thus, suspended sediments represented the major pathway of the riverine mercury transport. A mass balance calculation suggested an annual mercury flux of 48.6 g/km² into the Hunza River basin. The samples collected from the most affected river, the Shimsal River, averaged to have 108 ng/l total mercury. This amount was close to the average soil mercury data of 151 ng/l as reported by the Pakistan Mineral Development Corporation in 2001. The dominant source of contamination was shown to be the leaching of large quantities of mercury from the mercury‐rich sediment and flood plain soil into the rivers, rather than the direct release from mining activities. Significant decrease in both dissolved and particulate‐bound mercury concentration downstream of Attabad Lake suggested that mercury is being accumulated or consumed in the lake. Although minimization or elimination of mercury loses from the mining process seems important for the well‐being of the miners, preventing the remobilization of accumulated mercury is equally important in mercury control in this region. Copyright © 2014 John Wiley & Sons, Ltd.