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.     

Assessing the Effect of Data Imports on the Completeness of OpenStreetMap – A United States Case Study

The assessment of OpenStreetMap (OSM) data quality has become an interdisciplinary research area over the recent years. The question of whether the OSM road network should be updated through periodic data imports from public domain data, or whether the currency of OSM data should rather rely on more traditional data collection efforts by active contributors, has led to perpetual debates within the OSM community. A US Census TIGER/Line 2005 import into OSM was accomplished in early 2008, which generated a road network foundation for the active community members in the US. In this study we perform a longitudinal analysis of road data for the US by comparing the development of OSM and TIGER/Line data since the initial TIGER/Line import. The analysis is performed for the 50 US states and the District of Columbia, and 70 Urbanized Areas. In almost all tested states and Urbanized Areas, OSM misses roads for motorized traffic when compared with TIGER/Line street data, while significant contributions could be observed in pedestrian related network data in OSM compared with corresponding TIGER/Line data. We conclude that the quality of OSM road data could be improved through new OSM editor tools allowing contributors to trace current TIGER/Line data.     

Tsunami hazard and exposure on the global scale

In the aftermath of the 2004 Indian Ocean tsunami, a large increase in the activity of tsunami hazard and risk mapping is observed. Most of these are site-specific studies with detailed modelling of the run-up locally. However, fewer studies exist on the regional and global scale. Therefore, tsunamis have been omitted in previous global studies comparing different natural hazards. Here, we present a first global tsunami hazard and population exposure study. A key topic is the development of a simple and robust method for obtaining reasonable estimates of the maximum water level during tsunami inundation. This method is mainly based on plane wave linear hydrostatic transect simulations, and validation against results from a standard run-up model is given. The global hazard study is scenario based, focusing on tsunamis caused by megathrust earthquakes only, as the largest events will often contribute more to the risk than the smaller events. Tsunamis caused by non-seismic sources are omitted. Hazard maps are implemented by conducting a number of tsunami scenario simulations supplemented with findings from literature. The maps are further used to quantify the number of people exposed to tsunamis using the Landscan population data set. Because of the large geographical extents, quantifying the tsunami hazard assessment is focusing on overall trends.     

Southern Hemisphere extra-tropical forcing: a new paradigm for El Ni?o-Southern Oscillation

The main goal of this paper is to shed additional light on the reciprocal dynamical linkages between mid-latitude Southern Hemisphere climate and the El Ni?o-Southern Oscillation (ENSO) signal. While our analysis confirms that ENSO is a dominant source of interannual variability in the Southern Hemisphere, it is also suggested here that subtropical dipole variability in both the Southern Indian and Atlantic Oceans triggered by Southern Hemisphere mid-latitude variability may also provide a controlling influence on ENSO in the equatorial Pacific. This subtropical forcing operates through various coupled air?Csea feedbacks involving the propagation of subtropical sea surface temperature (SST) anomalies into the deep tropics of the Atlantic and Indian Oceans from boreal winter to boreal spring and a subsequent dynamical atmospheric response to these SST anomalies linking the three tropical basins at the beginning of the boreal spring. This atmospheric response is characterized by a significant weakening of the equatorial Atlantic and Indian Inter-Tropical Convergence Zone (ITCZ). This weakened ITCZ forces an equatorial ??cold Kelvin wave?? response in the middle to upper troposphere that extends eastward from the heat sink regions into the western Pacific. By modulating the vertical temperature gradient and the stability of the atmosphere over the equatorial western Pacific Ocean, this Kelvin wave response promotes persistent zonal wind and convective anomalies over the western equatorial Pacific, which may trigger El Ni?o onset at the end of the boreal winter. These different processes explain why South Atlantic and Indian subtropical dipole time series indices are highly significant precursors of the Ni?o34 SST index several months in advance before the El Ni?o onset in the equatorial Pacific. This study illustrates that the atmospheric internal variability in the mid-latitudes of the Southern Hemisphere may significantly influence ENSO variability. However, this surprising relationship is observed only during recent decades, after the so-called 1976/1977 climate regime shift, suggesting a possible linkage with global warming or decadal fluctuations of the climate system.     

Mapping of erosion rates in marly badlands based on a coupling of anatomical changes in exposed roots with slope maps derived from LiDAR data

Black marls form very extensive outcrops in the Alps and constitute some of the most eroded terrains, thus causing major problems of sedimentation in artificial storage systems (e.g. reservoirs) and river systems. In the experimental catchments near Draix (France), soil erosion rates have been measured in the past at the plot scale through a detailed monitoring of surface elevation changes and at the catchment scale through continuous monitoring of sediment yield in traps at basin outlets. More recently, erosion rates have been determined by means of dendrogeomorphic techniques in three monitored catchments of the Draix basin. A total of 48 exposed roots of Scots pine have been sampled and anatomical variations in annual growth rings resulting from denudation analysed. At the plot scale, average medium‐term soil erosion rates derived from exposed roots vary between 1·8 and 13·8 mm yr^?1 (average: 5·9 mm yr^?1) and values are significantly correlated with slope angle. The dendrogeomorphic record of point‐scale soil erosion rates matches very well with soil erosion rates measured in the Draix basins. Based on the point‐scale measurements and dendrogeomorphic results obtained at the point scale, a linear regression model involving slope angle was derived and coupled to high‐resolution slope maps obtained from a LiDAR‐generated digital elevation model so as to generate high‐resolution soil erosion maps. The resulting regression model is statistically significant and average soil erosion rates obtained from the areal erosion map (5·8, 5·2 and 6·2 mm yr^?1 for the Roubine, Moulin and Laval catchments, respectively) prove to be well in concert with average annual erosion rates measured in traps at the outlet of these catchments since 1985 (6·3, 4·1 and 6·4 mm yr^?1). This contribution demonstrates that dendrogeomorphic analyses of roots clearly have significant potential and that they are a powerful tool for the quantification and mapping of soil erosion rates in areas where measurements of past erosion is lacking. Copyright © 2011 John Wiley & Sons, Ltd.     

Coupled discrete element and finite volume solution of two classical soil mechanics problems

One dimensional solutions for the classic critical upward seepage gradient/quick condition and the time rate of consolidation problems are obtained using coupled routines for the finite volume method (FVM) and discrete element method (DEM), and the results compared with the analytical solutions. The two phase flow in a system composed of fluid and solid is simulated with the fluid phase modeled by solving the averaged Navier–Stokes equation using the FVM and the solid phase is modeled using the DEM. A framework is described for the coupling of two open source computer codes: YADE-OpenDEM for the discrete element method and OpenFOAM for the computational fluid dynamics. The particle–fluid interaction is quantified using a semi-empirical relationship proposed by Ergun [12]. The two classical verification problems are used to explore issues encountered when using coupled flow DEM codes, namely, the appropriate time step size for both the fluid and mechanical solution processes, the choice of the viscous damping coefficient, and the number of solid particles per finite fluid volume.     

Extension tardi-orogénique et formation des bassins intracontinentaux: le bassin stéphanien des Cévennes

Résumé

Le bassin carbonifère des Cévennes, voisin du décrochement sénestre de Villefort, a été étudié en intégrant les données sédimentologiques, les données structurales du socle et du remplissage, ainsi que les données de la pétrologie métamorphique du socle. Le remplissage du bassin est contrôlé par des failles subméridiennes induisant une subsidence localisée en début d’ouverture. Cette subsidence devient plus régionale en fin d’histoire. Un grenat mangane- sifère est associé aux minéraux du faciès schistes verts, des blastes d’andalousite recoupent l’ensemble. Cette association minéralogique indique un gradient de température de 50°C/km. Ce gradient est semblable au gradient estimé dans le bassin sur la matière organique. La schistosité porte trois familles de linéation : 1- une famille à biotite-chlorite-quartz, orientée N45 ± 20° et associée à la mise en place des nappes cévenoles; 2- une famille à minéraux phylliteux dont l’orientation est située autour de N90 et associée à l’ouverture du bassin; 3- une famille à minéraux phylliteux orientée NO horizontale sur des plans de schistosité redressés près de la faille de Villefort. La schistosité régionale est affectée par des plis asymétriques, des bandes de cisaillement et des zones cata- clasées, le tout ayant une cinématique en faille normale vers l’Est. Les relations entre ces différentes structures suggèrent un continuum de déformation depuis des niveaux ductiles jusqu’à des niveaux fragiles. Le remplissage du bassin est affecté par des décrochements parallèles à la faille de Ville- fort. Le bassin des Cévennes s’est donc ouvert dans un contexte extensif est-ouest sur une croûte épaissie. La région a ensuite été reprise par des décrochements nord-sud dont celui de Villefort.     

Retention behaviour of humic substances onto mineral surfaces and consequences upon thorium (IV) mobility: case of iron oxides

The influence of humic acids (HA) on the retention of Th(IV) onto the surface of Fe oxides has been investigated. Thorium (IV) is here considered as an analogue of tetravalent radionuclides (Zr, Hf, U, Np, Pu). Hematite (α–Fe₂O₃) and gœthite (α–FeOOH) have been chosen as model phases because of their widespread presence in the earth's crust. Retention experiments were performed by batch procedure at constant ionic strength and varying pH in the range 2–10; the HA concentration varied in the range 1–100 mg l⁻¹. Compared with the former data obtained on the retention of Th(IV) with the same hematite sample free of HA, the study of a ternary system (i.e. hematite–HA–Th), clearly shows the influence of HA on the retention mechanisms of Th(IV) even in neutral to slightly alkaline medium. In particular, it is illustrated that the retention of Th(IV) onto Fe oxides closely depends upon the ratio between humic and surface (Fe oxide) sites. This ratio has to be greater than a critical value in order to impede the sorption of Th(IV).     

Receiver function method in reflection seismology

The receiver function method was originally developed to analyse earthquake data recorded by multicomponent (3C) sensors and consists in deconvolving the horizontal component by the vertical component. The deconvolution process removes travel path effects from the source to the base of the target as well as the earthquake source signature. In addition, it provides the possibility of separating the emergent P and PS waves based on adaptive subtraction between recorded components if plane waves of constant ray parameters are considered. The resulting receiver function signal is the local PS wave's impulse response generated at impedance contrasts below the 3C receiver.We propose to adapt this technique to the wide‐angle multi‐component reflection acquisition geometry. We focus on the simplest case of land data reflection acquisition. Our adapted version of the receiver function approach consists in a multi‐step procedure that first removes the P wavefield recorded on the horizontal component and next removes the source signature. The separation step is performed in the τ?p domain while the source designature can be achieved in either the τ?p or the t?x domain. Our technique does not require any a priori knowledge of the subsurface. The resulting receiver function is a pure PS‐wave reflectivity response, which can be used for amplitude versus slowness or offset analysis. Stack of the receiver function leads to a high‐quality S wave image.     

Contrasting regional discharge evolutions in the Amazon basin (1974–2004)

Former hydrological studies in the Amazon Basin generally describe annual discharge variability on the main stem. However, the downstream Amazon River only represents the mean state of the Amazonian hydrological system. This study therefore uses a new data set including daily discharge in 18 sub-basins to analyze the variability of regional extremes in the Amazon basin, after recalling the diversity of the hydrological annual cycles within the Amazon basin. Several statistical tests are applied in order to detect trends and breaks in the time series. We show that during the 1974–2004 period, the stability of the mean discharge on the main stem in Óbidos is explained by opposite regional features that principally involve Andean rivers: a decrease in the low stage runoff, particularly important in the southern regions, and an increase in the high stage runoff in the northwestern region. Both features are observed from the beginning of the nineties. These features are also observed in smaller meridian sub-basins in Peru and Bolivia. Moreover we show that the changes in discharge extremes are related to the regional pluriannual rainfall variability and the associated atmospheric circulation as well as to tropical large-scale climatic indicators.