Reconstructing historic Glacial Lake Outburst Floods through numerical modelling and geomorphological assessment: Extreme events in the Himalaya

Recession of high‐mountain glaciers in response to climatic change frequently results in the development of moraine‐dammed glacial lakes. Moraine dam failure is often accompanied by the release of large volumes of water and sediment, termed a Glacial Lake Outburst Flood (GLOF). Chukhung Glacier is a small (~3 km²) receding valley glacier in Mt. Everest (Sagarmatha) National Park, Nepal. Unlike many Himalayan glaciers, which possess a thick mantle of supraglacial debris, its surface is relatively clean. The glacier terminus has receded 1.3 km from its maximum Holocene position, and in doing so provided the space for an ice‐contact moraine‐dammed lake to develop. The lake had a maximum volume of 5.5 × 10⁵ m³ and drained as a result of breaching of the terminal moraine. An estimated 1.3 × 10⁵ m³ of material was removed from the terminal moraine during breach development. Numerical dam‐breach modelling, implemented within a Generalised Likelihood Uncertainty Estimation (GLUE) framework, was used to investigate a range of moraine‐dam failure scenarios. Reconstructed outflow peak discharges, including failure via overtopping and piping mechanisms, are in the range 146–2200 m³ s^‐1. Results from two‐dimensional hydrodynamic GLOF modelling indicate that maximum local flow depths may have exceeded 9 m, with maximum flow velocities exceeding 20 m s^‐1 within 700 m of the breach. The floodwaters mobilised a significant amount of material, sourced mostly from the expanding breach, forming a 300 m long and 100 m wide debris fan originating at the breach exit. moraine‐dam. These results also suggest that inundation of the entire floodplain may have been achieved within ten minutes of initial breach development, suggesting that debris fan development was rapid. We discuss the key glaciological and geomorphological factors that have determined the evolution of a hazardous moraine‐dammed lake complex and the subsequent generation of a GLOF and its geomorphological impact. © 2014 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Concentration–discharge relationships in the coal mined region of the New River basin and Indian Fork sub‐basin,Tennessee, USA

For many basins, identifying changes to water quality over time and understanding current hydrologic processes are hindered by fragmented and discontinuous water‐quality and hydrology data. In the coal mined region of the New River basin and Indian Fork sub‐basin, muted and pronounced changes, respectively, to concentration–discharge (C–Q) relationships were identified using linear regression on log‐transformed historical (1970s–1980s) and recent (2000s) water‐quality and streamflow data. Changes to C–Q relationships were related to coal mining histories and shifts in land use. Hysteresis plots of individual storms from 2007 (New River) and the fall of 2009 (Indian Fork) were used to understand current hydrologic processes in the basins. In the New River, storm magnitude was found to be closely related to the reversal of loop rotation in hysteresis plots; a peak‐flow threshold of 25 cubic meters per second (m³/s) segregates hysteresis patterns into clockwise and counterclockwise rotational groups. Small storms with peak flow less than 25 m³/s often resulted in dilution of constituent concentrations in headwater tributaries like Indian Fork and concentration of constituents downstream in the mainstem of the New River. Conceptual two or three component mixing models for the basins were used to infer the influence of water derived from spoil material on water quality. Copyright © 2012 John Wiley & Sons, Ltd.     

Geochemical surveillance of active volcanoes: data on the fumaroles of Vulcano (Aeolian Islands,Italy)

Temperature and chemical composition of condensates of higher temperature fumaroles, located on the rim of the crater of Vulcano, were studied over a period of one year. Maximum temperature changed from 216°C (June 1977) to 291°C (May 1978). HCl is the main constituent of condensates which, in comparison with data from fumaroles of other volcanoes, contain also noteworthy quantities of boron and bromine. Since these chemical characteristics suggested the possible inflow of water of marine origin into the volcanic conduit, experimental data on seawater-rock reactions at high temperatures and pressures were taken into account. A model is drafted for the volcanic system, in which the occurrence of an aquifer interposed between the magma chamber and the surface is considered. The vaporization of this aquifer would feed the fumaroles giving rise, through different mixing processes with surface waters, to the observed differences in temperature and chemical composition. However, data available at present are not sufficient to support the proposed model, which must be mainly considered as a working hypothesis.     

ON THE CHARACTERISTICS OF NON-LINEAR SOIL RESPONSE AND DYNAMIC SOIL PROPERTIES USING VERTICAL ARRAY DATA IN JAPAN

Non-linear response of the soil is investigated by comparing the spectral ratios (uphole/downhole) using weak and strong motions. Data from seven vertical arrays in Japan are analysed in this study. The frequency-dependent transfer function of soil is calculated as a ratio of the spectrum at uphole to the spectrum at downhole, considering the horizontal component of shear wave. In spectral ratio analysis auto- and cross-spectra are employed. The reduction in the predominant frequency of the transfer function with increases in excitation level reflects the non-linear response of the soil. Results of analysis demonstrate a significant non-linear ground response at six sites with surface PGA exceeding 90 gal. However, the results of one site show the linear response up to 130 gal surface PGA. Furthermore, the in situ strain-dependent soil behaviour is examined through the shear modulus – shear strain relationship. When compared, the actual and laboratory results of the shear strain – shear modulus relationship are in agreement. Additionally, a good consistency between the tendency of reduction in shear modulus ratio with shear strain increases, and reduction of predominant frequency with ground motion increases, confirms the significance of non-linearity in site effects study.     

The evolution of slump-induced destruction of Kukuy Griva slope (Lake Baikal) revealed on the base of the data of seismic and acoustic surveys

As a result of careful interpretation of the data of recent seismo-acoustic surveys, two major seismic complexes were identified in the structure of the upper sedimentary section of the north-western slope of Kukuy Griva (ridge). They are composed of several inner seismo-facies, which are very characteristic and separated with well-expressed regional reflector. Seismic and acoustic data revealed evidences of numerous landslide processes which were different in age and peculiarities of manifestation. For the first time, detailed schemes of seismo-facies distribution are compiled for major seismic complexes at the area. The boundary between seismic complexes is dated as 150 thousands years. It is shown that sliding was more intensive during Late Pleistocene and Holocene, reflecting activation of tectonic movements in this part of Baikal rift system.     

Geochemical processes controlling the spatial distribution of selenium in soil and water,west central South Dakota,USA

Selenium (Se) is essential in the human diet, but has a low threshold for toxic concentration. It is recommended that nutrients such as Se should be consumed through foods as part of a normal diet. Se concentrations in crops and meat depend on the amount of labile Se in the soil where crops and forage are grown. Therefore, managing agriculture for optimal Se in grain crops and forage requires an understanding of the distribution and mobility of Se. Elevated concentrations of Se occur in waters, soils, and forage 120 km west of Pierre in west central South Dakota, USA. The research site lies in an elevated, dissected plain where soils developed on gently dipping Pierre Shale. Soils were sampled along catena transects and waters collected from soil, ponds, and shallow borings in areas of known elevated forage and crop Se. Soil extracts from saturated-paste extraction and acid (aqua regia and hydrofluoric acid) extraction were analyzed. Selenate was the dominant Se species in both acid and saturated-paste extracts; selenite and organic Se were below detection (<0.2 ppb) in the same soil extracts. On average, 98% of soil Se was not water-soluble. The distribution of total Se shows much less spatial variation than water-soluble Se in the areas sampled. Se shows correlation with organic carbon in soils and waters, suggesting its association with organic carbon. Ca shows some correlation to Se in acid extracts, but not in saturated-paste extracts or in waters. Total Se shows no significant correlation to Na, Mg, and total S in the soils. Se in saturated-paste extracts and water samples shows good correlations with Na, Mg, and SO₄, suggesting that evaporitic Na–Mg–sulfate minerals may temporarily concentrate water-soluble Se in shallow soils. The dissolution and precipitation of these Na–Mg–sulfate salts together with pH and oxidation–reduction conditions apparently control water-soluble Se distribution and mobilization in shallow subsurface environments.     

Investigation on strength and stability of jointed rock mass using three‐dimensional numerical manifold method

This paper proposes a numerical model for jointed rock masses within the 3‐D numerical manifold method (NMM) framework equipped with a customized contact algorithm. The strength of rock sample containing a few sets of discontinuities is first investigated. The results of models with simple geometries are compared with the available analytical solutions to verify the developed computer code, whereas models with complex geometries are simulated to better understand the fundamental behavior and failure mechanism of jointed rock mass. Furthermore, the stability of jointed rock mass in an underground excavation is studied, where rock failure process is determined by the 3‐D NMM simulation. The simulation results provide valuable guidance on excavation process design and stabilization design in rock engineering practice. Copyright © 2012 John Wiley & Sons, Ltd.     

Characterising the distribution of near-shore submarine groundwater discharge along a coastline using 222Rn and electrical conductivity

The spatial and temporal distribution of near-shore fresh submarine groundwater discharge (SGD) was characterised from the coastal aquifers of the Willunga Basin, South Australia, an extensive aquifer system that supports an important viticultural region. Measurements of electrical conductivity (EC) and ²²²Rn (radon) activity were collected at 19 sites along the coastline during the Southern Hemisphere spring (2011) and summer (2013). At each site, samples were collected from the surf zone as well asporewater from beach sediment in the intertidal zone. Surf-zone radon activity ranged from <5 to 70mBq L^–1, and intertidal porewater radon ranged over two orders of magnitude (220–36 940 mBq L^–1) along the Willunga Basin coastline during both surveys. Overall, surf-zone and porewater EC was lower in the spring 2011 survey than in the summer 2013 survey. Porewater EC was similar to that of coastal water at most sites along the coastline, except at three sites where porewater EC was found to be lower than coastal water during both surveys, and three sites where evaporated seawater was observed in the summer survey. Based on the patterns in radon and EC along the coastline, two sites of localised fresh SGD were identified, in addition to a groundwater spring that is known to discharge to the coast. The results indicate that near-shore fresh SGD occurs as localised seeps rather than diffuse seepage along the entire coastline. The apparent absence of groundwater discharge at most locations is also consistent with current evidence suggesting that extensive groundwater pumping within the basin has resulted in seawater intrusion across much of the coastline. These observations also suggest that previous studies are likely to have over-estimated SGD rates from the Willunga Basin because they assumed that SGD occurred along the entire coastline.