The effects of the 2004 tsunami on a coastal aquifer in Sri Lanka

On December 26, 2004, the earthquake off the southern coast of Sumatra in the Indian Ocean generated far-reaching tsunami waves, resulting in severe disruption of the coastal aquifers in many countries of the region. The objective of this study was to examine the impact of the tsunami on groundwater in coastal areas. Field investigations on the east coast of Sri Lanka were carried out along a transect located perpendicular to the coastline on a 2.4 km wide sand stretch bounded by the sea and a lagoon. Measurements of groundwater table elevation and electrical conductivity (EC) of the groundwater were carried out monthly from October 2005 to August 2007. The aquifer system and tsunami saltwater intrusion were modeled using the variable-density flow and solute transport code HST3D to understand the tsunami plume behavior and estimate the aquifer recovery time. EC values reduced as a result of the monsoonal rainfall following the tsunami with a decline in reduction rate during the dry season. The upper part of the saturated zone (down to 2.5 m) returned to freshwater conditions (EC < 1000 μS/cm) 1 to 1.5 years after the tsunami, according to field observations. On the basis of model simulations, it may take more than 15 years for the entire aquifer (down to 28 m) to recover completely, although the top 6 m of the aquifer may become fresh in about 5 years.     

Simulating Martian regolith in the laboratory

Regolith and dust cover the surfaces of the Solar Systems solid bodies, and thus constitute the visible surface of these objects. The topmost layers also interact with space or the atmosphere in the case of Mars, Venus and Titan. Surface probes have been proposed, studied and flown to some of these worlds. Landers and some of the mechanisms they carry, e.g. sampling devices, drills and subsurface probes (“moles”) will interact with the porous surface layer. The absence of true extraterrestrial test materials in ample quantities restricts experiments to the use of soil or regolith analogue materials. Several standardized soil simulants have been developed and produced and are commonly used for a variety of laboratory experiments. In this paper we intend to give an overview of some of the most important soil simulants, and describe experiments (penetrometry, thermal conductivity, aeolian transport, goniometry, spectroscopy and exobiology) made in various European laboratory facilities.     

Ion microprobe studies of water in silicate melts. Concentration-dependent water diffusion in obsidian

The ion microprobe at Johnson Space Center has been calibrated for in situ water determinations on a 10-μm scale over the range 0.2 wt.% H₂O to 1.8, 6.8, and 3.7 wt.%, for basaltic, albitic, and rhyolitic glasses, respectively. The basalt glass calibration curve differs substantially from those of albite and rhyolite glasses, indicating a need to carefully match composition and/or melt structure between H₂O standards and unknowns.A value for the diffusivity of water as a function of concentration and time has been calculated from water diffusion profiles measured in rhyolite glasses prepared at 850°C and 700 barsP_t(H₂O) [1]. Transient diffusion into a semi-infinite medium is described by the equation:?(φ/2)?¸/?φ=?(D_w?¸/?φ)/?φ #x003B8;=1, φ=0, θ→ 0, θ→∞, wherex =distance from the cylinder edge,t =time,C₀ =initial concentration,C_s =concentration at the edge,C =concentration at x,θ = C ? C₀/C_s ? C₀,φ = x/t^1/2, andD_w =diffusivity of water. An iterative technique has been used to calculate solutions to the diffusion equation as a function ofD_w [2]. Comparison of these solutions with the ion probe data indicate that, for0.2wt.% ≤ C ≤ 3.7wt.%H₂O,D_w can be described by an exponential function of θ, of the formD_w = D₀exp(bθ), withD₀ (i.e.,D_w at 0.2%) = (0.8?2.2) × 10^?8 cm²/s and2 ≤ b ≤ 4.     

Ion microprobe studies of water in silicate melts: Temperature-dependent water diffusion in obsidian

The temperature dependence of water diffusivity in rhyolite melts over the range 650–950°C and [P_T(H₂O] = 700 bars is evaluated from water concentration-distance profiles measured in glass with an ion microprobe. Diffusivities are exponentially dependent on concentration over this temperature range and vary from about 10^?8 cm²/s at 650°C to about 10^?7 cm²/s at 950°C at 2 wt.% water. Water solubility also varies with temperature at a rate of ?0.14 wt. per 100°C increase. The avtivation energy (E_a) appears to be constant at 19 ± 1kal/mole for 1, 2,and 3 wt.% H₂O. Comparison of these data with results for cation diffusion indicates that this value is a minimum E_a for diffusion of any species in a rhyolite melt.Compensation plots of log₁₀D₀ (the frequency factor) versus E_a indicate that hydrous rhyolite melts follow the same trend as anhydrous basalts. D₀ increases for H₂O and Ca²⁺ [1] as E_a decreases. This suggests that these molecules may diffuse by different mechanisms than do monovalent cations, and that hydration of the melt affects diffusion of Ca²⁺ and H₂O differently than it does monovalent cation diffusion. The results imply that dramatic increases in cation diffusivities by hydration [1] may occur with additions of less than 1 wt.% H₂O.     

Strain localisation in bimineralic rocks: Experimental deformation of synthetic calcite–anhydrite aggregates

Deformation of synthetic calcite–anhydrite aggregates to large shear strains (up to γ = 12.4 at 600 °C, 300 MPa confining pressure and a constant angular displacement rate corresponding to a shear strain rate of 10^− 3 s^− 1) resulted in the first experimental observation of strain localisation from initially homogeneous rocks. In contrast to experiments on pure calcite and anhydrite, which deformed homogeneously to large strains (γ ≥ 5), all experiments on calcite–anhydrite mixtures resulted in heterogeneous deformation at γ > 1 and the formation of narrow localised bands in the microstructures at γ > 4. In these bands, the amount of strain is at least twice as large as in the rest of the sample and individual grains of the same phase cluster and align, thereby forming microstructural layering similar to planar fabrics in natural mylonites. A switch in deformation mechanism in anhydrite from dislocation creep to diffusion creep and/or grain boundary sliding occurs simultaneously with strain localisation. It is concluded that deformation-induced heterogeneous phase distributions cause local strength differences initiating strain localisation in the calcite–anhydrite mixtures. The study suggests that the presence of two phases in combination with a change in deformation mechanism may be responsible for strain localisation in natural poly-mineralic mylonites.     

Crystal field transitions in Co2[Al4Si5]O18 cordierite and CoAl2O4 spinel determined by neutron spectroscopy

Inelastic magnetic neutron scattering has been used to determine the energy of the ⁴ A ₂→⁴ T ₂ transition in CoAl₂O₄ spinel and the δ₁ transition in Co₂[Al₄Si₅]O₁₈ cordierite. The observed crystal field splitting in Co-spinel is 485 meV (3900 cm⁻¹), which corresponds to a crystal field stabilization energy of 56.2 kJmol⁻¹. The transition energy of the δ₁ transition in Co-cordierite has been determined to be 21 meV (170 cm⁻¹). The present data demonstrate that magnetic neutron scattering can be used to measure crystal field transitions at energies of interest in the study of 3d-containing silicates. It may be used to measure transition energies when the use of optical spectroscopy is inappropriate. Received: 30 January 1997 / Accepted: 5 July 1997     

Global Warming and the Summertime Evapotranspiration Regime of the Alpine Region

Changes of the summer evapotranspiration regime under increased levels of atmospheric greenhouse gases are discussed for three Alpine river basins on the basis of a new set of simulations carried out with a high-resolution hydrological model. The climate change signal was inferred from the output of two simulations with a state-of-the-art global climate model (GCM), a reference run valid for 1961–1990 and a time-slice simulation valid for 2071–2100 under forcing from the A2 IPCC emission scenario. In this particular GCM experiment and with respect to the Alpine region summer temperature was found to increase by 3 to 4 ^∘C, whereas precipitation was found to decrease by 10 to 20%. Global radiation and water vapor pressure deficit were found to increase by about 5% and 2 hPa, respectively. On this background, an overall increase of potential evapotranspiration of about 20% relative to the baseline was predicted by the hydrological model, with important variations between but also within individual basins. The results of the hydrological simulations also revealed a reduction in the evapotranspiration efficiency that depends on altitude. Accordingly, actual evapotranspiration was found to increase at high altitudes and to the south of the Alps, but to decrease in low elevation areas of the northern forelands and in the inner-Alpine domain. Such a differentiation does not appear in the GCM scenario, which predicts an overall increase in evapotranspiration over the Alps. This underlines the importance of detailed simulations for the quantitative assessment of the regional impact of climate change on the hydrological cycle.     

The importance of alternative conceptual models for simulation of concentrations in a multi-aquifer system

Four different conceptual models based on alternative geological interpretations were formulated for a shallow 600 km² aquifer system in Denmark comprising Quaternary deposits. Each of the four models was calibrated against groundwater heads and discharge measurements through inverse modeling. Subsequently, the transport capabilities of the four models were compared to 32 concentration measurements of environmental tracers (tritium ³H, helium-3 ³He, chlorofluorocarbons CFC11, CFC12 and CFC113). The flow simulations showed only minor differences in spatial head distribution associated with alternative conceptualizations despite the complexity of the aquifer system and the significant differences in geological interpretations. The models, however, showed major differences in predictions of the age of the groundwater and environmental tracer concentrations, differences that are seen as an effect of model structure uncertainty, because no additional calibrations to these data were performed. A single conceptualization may be adequate in characterizing the natural behavior of a field system after calibration, because the calibration procedure is able to compensate for errors in the data or in the conceptual model through biased parameter values. However, once extrapolation beyond the calibration base is attempted, different conceptual model formulations result in significantly different results. Consequently, it is crucial to take model conceptual uncertainty into account when making predictions beyond the calibration base.