Thermal anomalies associated with shallow gas hydrates in the K-2 mud volcano, Lake Baikal

Thermal measurements and hydrate mapping in the vicinity of the K-2 mud volcano in Lake Baikal have revealed a particular type of association of thermal anomalies (29–121?mW?m^–2) near hydrate-forming layers. Detailed coring within K-2 showed that hydrates are restricted to two distinct zones at sub-bottom depths exceeding 70–300?cm. Temperature data from stations with hydrate recovery and degassing features all display low thermal gradients. Otherwise, the thermal gradients within the mud volcano are generally increased. These findings imply a more complicated thermal regime than often assumed for mud volcanoes, with important roles for both fluids and hydrates. The coexistence of neighbouring low and high thermal anomalies is interpreted to result from discharging and recharging fluid activity, rather than hydrate thermodynamics. It is suggested that hydrates play a key role in controlling the fluid circulation pattern at an early stage. At a later stage, the inflow of undersaturated lake water would favour the dissolution of structure I hydrates and the formation of structure II hydrates, the latter having been observed on top of structure I hydrates in the K-2 mud volcano.     

Quaternary alluvial-fan development, climate and morphologic dating of fault scarps in Laguna Salada, Baja California, Mexico

Late Quaternary slip across the Cañada David detachment has produced an extensive array of Quaternary scarps cutting alluvial-fans along nearly the entire length (~ 60 km) of the range-bounding detachment. Eight regional alluvial-fan surfaces (Q₁ [youngest] to Q₈ [oldest]) are defined and mapped along the entire Sierra el Mayor range-front. Terrestrial cosmogenic nuclide ¹⁰Be concentrations from individual boulders on alluvial-fan surfaces Q₄ and Q₇ yield surface exposure ages of 15.5 ± 2.2 ka and 204 ± 11 ka, respectively. Formation of the fans is probably tectonic, but their evolution is strongly moderated by climate, with surfaces developing as the hydrological conditions have changed in response to climate change on Milankovitch timescales. Systematic mapping reveals that the fault scarp array along active range-bounding faults in Sierras Cucapa and El Mayor can be divided into individual rupture zones, based on cross-cutting relationships with alluvial-fans. Quantitative morphological ages of the Laguna Salada fault-scarps, derived from linear diffusive degradation modeling, are consistent with the age of the scarps based on cross-cutting relationships. The weighted means of the maximum mass diffusivity constant for all scarps with offsets < 4 m is 0.051 and 0.066 m²/ka for the infinite and finite-slope solutions of the diffusion equation, respectively. This estimate is approximately an order of magnitude smaller than the lowest diffusivity constants documented in other regions and it probably reflects the extreme aridity and other microclimatic conditions that characterize the eastern margin of Laguna Salada.     

Seismic stratigraphy of Lago Puyehue (Chilean Lake District): new views on its deglacial and Holocene evolution

Prior to the collection of a series of sediment cores, a high- and very-high-resolution reflection seismic survey was carried out on Lago Puyehue, Lake District, South-Central Chile. The data reveal a complex bathymetry and basin structure, with three sub-basins separated by bathymetric ridges, bedrock islands and interconnected channels. The sedimentary infill reaches a thickness of >200 m. It can be sub-divided into five seismic-stratigraphic units, which are interpreted as: moraine, ice-contact or outwash deposits (Unit I), glacio-lacustrine sediments rapidly deposited in a proglacial or subglacial lake at the onset of deglaciation (Unit II), lacustrine fan deposits fed by sediment-laden meltwater streams in a proglacial lake (Unit III), distal deposits of fluvially derived sediment in an open, post-glacial lake (Unit IV) and authigenic lacustrine sediments, predominantly of biogenic origin, that accumulated in an open, post-glacial lake (Unit V). This facies succession is very similar to that observed in other glacial lakes, and minor differences are attributed to an overall higher depositional energy and higher terrigenous input caused by the strong seismic and volcanic activity in the region combined with heavy precipitation. A long sediment core (PU-II core) penetrates part of Unit V and its base is dated as 17,915 cal. yr. BP. Extrapolation of average sedimentation rates yields an age of ca. 24,750 cal. yr. BP for the base of Unit V, and of ca. 28,000 cal. yr. BP for the base of Unit IV or for the onset of open-water conditions. This is in contrast with previous glacial-history reconstructions based on terrestrial records, which date the complete deglaciation of the basin as ca. 14,600 cal. yr. BP. This discrepancy cannot be easily explained and highlights the need for more lacustrine records from this region. This is the second in a series of eight papers published in this special issue dedicated to the 17,900 year multi-proxy lacustrine record of Lago Puyehue, Chilean Lake District. The papers in this special issue were collected by M. De Batist, N. Fagel, M.-F. Loutre and E. Chapron.     

Morphology,structure and kinematics of a rainfall controlled slow‐moving Andean landslide,Peru

The large slow‐moving landslide of Maca is located in the upper Colca valley (southern Peru), a region characterized by a well pronounced rainy period, and intense and recurrent sustained seismicity. The landslide, developed in deep lacustrine deposits, has recently accelerated, threatening the Maca village. This work aims at understanding the rupture mechanism and the causes of the recent landslide reactivation/acceleration. We present a multidisciplinary characterization of the Maca landslide that includes: (i) geological and morphological mapping in the field; (ii) remote sensing analysis using an historical aerial photograph of 1955 and the Pléiades satellite images (2013); (iii) global positioning system (GPS) including time‐series of surveys over 13 years, and continuous measurements over 14 months; (iv) a geophysical campaign with deep electrical resistivity tomography profiles acquired across the landslide mass. Our study shows that this 60 Mm³ landslide, which can be classified as a clay/silt compound landslide, moved by 15 m between 2001 and 2014 with a large inter‐annual velocity variation (up to a factor of 500) depending on the rainfall intensity. We suggest that these dramatic changes in velocity are the result of the combination of a threshold mechanism and the short intense rainy season in Peru. This study reveals three main driving factors acting at different timescales: (i) over several decades, the river course has significantly changed, causing the Maca landslide reactivation in the 1980s due to the erosion of its toe; (ii) at the year scale, a minimum amount of rainfall is required to trigger the motion and this amount controls the landslide velocity; (iii) transient changes in slide velocity may occur anytime due to earthquakes. This study particularly highlights the non‐linear behaviour of the motion with rainfall. Copyright © 2016 John Wiley & Sons, Ltd.     

Review of self-potential methods in hydrogeophysics

The self-potential (SP) method is a passive geophysical method based on the natural occurrence of electrical fields on the Earth's surface. Combined with other geophysical methods, SP surveys are especially useful for localizing and quantifying groundwater flows and pollutant plume spreading, and estimating pertinent hydraulic properties of aquifers (water table, hydraulic conductivity). Laboratory experiments have shown that the involved coupling coefficients mainly depend on the fluid chemistry, conductivity and pH, and on the soil or rock properties. The interpretation of SP observations can be done qualitatively, for instance, by correlation of SP gradients with water fluxes (through electrokinetics) or salt fluxes (through electro-diffusion). In recent years, the interpretation has been improved with the help of modelling or/and inversion of the Poisson equation and endeavours to estimate hydraulic parameters by means of the intensity of electric current sources caused by underground flows.     

Influence of shallow infiltration on time-lapse ERT: Experience of advanced interpretation

Previous time-lapse Electrical Resistivity Tomography (ERT) studies have experienced difficulties in reconstructing reliable calculated resistivity changes in the subsurface. Increases or decreases of resistivity appear in the calculated ERT image where no changes were noted in the subsurface, leading to erroneous hydrological interpretations of the geophysical results. In this article, we investigate how a variation of actual resistivity with time and at shallow depth can influence time-lapse ERT results and produce resistivity artefacts at depth. We use 1 and 2-D numerical modelling to simulate infiltration scenarios. Using a standard time-lapse inversion, we demonstrate the resistivity artefact production according to the electrode spacing parameter. We used an advanced inversion methodology with a decoupling line at shallow depth to attenuate or remove resistivity artefacts. We also applied this methodology to a field data set obtained in a semi-arid environment in Burkina Faso, West Africa. Here, time-lapse ERT shows several resistivity artefacts of calculated resistivity if a standard inversion is used. We demonstrate the importance of a dense sampling of shallow resistivity variations at shallow depth. Advanced interpretation allows us to significantly attenuate or remove the resistivity artefact production at intermediate depth and produce reliable interpretation of hydrological processes.     

Contribution of geophysical surveys to groundwater modelling of a porous aquifer in semiarid Niger: An overview

Subsurface geophysical surveys were carried out using a large range of methods in an unconfined sandstone aquifer in semiarid south-western Niger for improving both the conceptual model of water flow through the unsaturated zone and the parameterization of numerical a groundwater model of the aquifer. Methods included: electromagnetic mapping, electrical resistivity tomography (ERT), resistivity logging, time domain electromagnetic sounding (TDEM), and magnetic resonance sounding (MRS). Analyses of electrical conductivities, complemented by geochemical measurements, allowed us to identify preferential pathways for infiltration and drainage beneath gullies and alluvial fans. The mean water content estimated by MRS (13%) was used for computing the regional groundwater recharge from long-term change in the water table. The ranges in permeability and water content obtained with MRS allowed a reduction of the degree of freedom of aquifer parameters used in groundwater modelling.