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.     

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.     

An external–internal mode coupling for a 3D hydrodynamical model for applications at regional scale (MARS)

This paper presents a 3D model in sigma coordinates. Although the principles it is based on have been established for some time, some original aspects for this type of 3D mode splitting model are presented here. The model was designed to simulate flows in various coastal areas from the regional scale down to the inshore scale of small bays or estuaries where circulation is generally driven by a mix of processes. The processes to be modeled enable simplifications of the Navier–Stokes equations on the classic Boussinesq and hydrostatic hypotheses. These equations are transformed within a sigma framework to make free surface processing easier. The main point of our demonstration focuses on the original aspect of the coupling between barotropic and baroclinic modes especially designed for ADI. It explains how full consistency of the transport calculated within the 2D and 3D equation sets was obtained. Lastly, we describe the physical processes simulated on a realistic configuration at a regional scale in the Bay of Biscay.     

Interpretation of Ni2p XPS spectra of Ni conductors and Ni insulators

Ni2p_3/2 X-ray photoelectron spectral peak binding energies of Ni metal, NiS, and NiAs (all conductors) span a range of about 0.5 eV and are, consequently, insensitive to formal Ni oxidation state and to the nature of the ligand to which Ni is bonded, relative to other metals (e.g., Fe). Ni2p_3/2 peak structures and binding energies reflect two energetic contributions. The major contribution is that associated with the electrostatic field produced by ejection of the Ni(2p) photoelectron, the minor contribution is the relaxation energy associated with filling unoccupied, conduction band 3d⁹ and 4s Ni metal orbitals. These conduction band orbitals become localized on the Ni photoion (and sometimes filled) in response to the field created by the photoemission event. Because only the core Ni2p electron and nonbonding orbitals of predominantly metallic character are affected, the main peak of all three conductors are affected similarly, leading to similar Ni2p_3/2 main peak binding energies. NiO, Ni(OH)₂, and NiSO₄ are insulators in which Ni is divalent and is bonded to oxygen. Although Ni is bonded to oxide in these phases, Ni2p binding energies differ substantially, and reflect primarily the nature of the ligand (O²⁻, OH⁻, SO₄ ²⁻) to which Ni is bonded. The influence of the ligand is the result of charge (electron) transfer from valence band bonding orbitals of dominantly ligand character, to unoccupied conduction band orbitals localized on Ni photoions. Relaxation energy resulting from charge transfer is acquired by the emitted photoelectron, thus Ni2p_3/2 photopeak binding energies of these insulators reflect the nature of the ligand to which Ni is bonded. The Ni2p main peak binding energy of these conductors and insulators is a poor guide to Ni oxidation states. The Ni2p_3/2 binding energies of insulators reflect, however, the nature of the ligand in the first coordination sphere of Ni. The intensity of the Doniach–Sunjic contribution to Ni2p XPS spectra of NiS and NiAs is dependent on the nature of the ligand. The Doniach–Sunjic contribution to ligand XPS core-level photopeaks (e.g., S2p of NiS and As3d of NiAs) has not been explained and is poorly understood. Received 24 May 1999 / Revised, accepted: 30 June 1999     

LA‐ICP‐MS Analysis of Small Samples: Carbonate Reference Materials and Larval Fish Otoliths

This article proposes a methodology to analyse the composition of very small carbonate samples such as larval fish otoliths. The chemical composition of otoliths, which are carbonate structures in the inner ear, is often used to explore population dynamics in fishes. Recent advances in laser ablation‐inductively coupled plasma‐mass spectrometry have suggested its potential application to this field. In this study, analyses were performed using a 193 nm ArF Resonetics LA system, coupled to an Agilent 7700X‐ICP‐MS, with the following ablation parameters: a beam diameter of 5 μm, energy of 3 mJ, 2.7 J cm^?2, laser repetition rate of 10 Hz and translation speed of 2.5 μm s^?1. NIST SRM 610 glass was used as the primary calibration material. Performing this protocol, characterisation of a USGS GP‐4 reference material was achieved with suitable precision and accuracy, but the USGS MACS‐3 reference material appeared more heterogeneous under the ablation conditions tested. Calibration was performed using two different beam diameters (5 and 11 μm). Capelin (Mallotus villosus) otoliths measuring between 10 and 20 μm in diameter were tested. Even though a smaller beam diameter and lower energy were used compared with those normally employed to analyse larger otoliths, the method was successful.