Origin and type of rainfall for recharge of a karstic aquifer in the western Mediterranean: a case study from the Sierra de Gador–Campo de Dalias (southeast Spain)

Isotope signatures in precipitation from the Global Network for Isotopes in Precipitation around the Mediterranean basin and literature data are compared with isotopic data from a large karstic aquifer in southeast Spain to explain the origin and type of the precipitation events dominating recharge. Analysis of the deuterium excess d at the scale of the Mediterranean basin and at the regional scale allows us to understand the isotopic context of the study area: Campo de Dalias and the Sierra de Gador (Almería province). The origin of precipitation can be determined from its d value. The d value changes as a function of the initial evaporation condition. It depends on the relative humidity and temperature during the evaporation producing the water vapour of the clouds. The water vapour, which dominates the study area, is generated in two areas: the Atlantic Ocean (d = 10‰) and the western Mediterranean basin (d = 15‰). With increasing precipitation volume, the western Mediterranean character dominates. These heavier storms contribute mainly to recharge, as illustrated by the d value of 13·6‰ in deep groundwater of the Campo de Dalias. Weighted d values increase with the volume of precipitation, giving a significant relationship for the southern and eastern coasts of the Iberian Peninsula. This selectivity of d to monthly precipitation was used to estimate the return period of precipitation leading to aquifer recharge at 0·9–4·9 years. Moderate rainfall, which occurs more frequently, still represents ～60–90% of the total precipitation. One of the challenges to meet ever‐growing water demands is to increase recharge from moderate events yielding intermediate quantities per event, but forming the bulk of the annual precipitation. Copyright © 2006 John Wiley & Sons, Ltd.     

InSAR datum connection using GNSS-augmented radar transponders

Deformation estimates from Interferometric Synthetic Aperture Radar (InSAR) are relative: they form a ‘free’ network referred to an arbitrary datum, e.g. by assuming a reference point in the image to be stable. However, some applications require ‘absolute’ InSAR estimates, i.e. expressed in a well-defined terrestrial reference frame, e.g. to compare InSAR results with those of other techniques. We propose a methodology based on collocated InSAR and Global Navigation Satellite System (GNSS) measurements, achieved by rigidly attaching phase-stable millimetre-precision compact active radar transponders to GNSS antennas. We demonstrate this concept through a simulated example and practical case studies in the Netherlands.     

Numerical modelling of multiphase immiscible flow in double‐porosity featured groundwater systems

A numerical model describing the flow of multiphase, immiscible fluids in a deformable, double‐porosity featured soil has been developed. The model is focused on the modelling of the secondary porosity features in soil, which is more relevant to groundwater contamination problems. The non‐linear saturation and relative permeabilities were expressed as functions of the capillary pressures. The governing partial differential equations in terms of soil displacement and fluid pressures were solved numerically. Galerkin's weighted‐residual finite element method was employed to obtain the spatial discretization whereas temporal discretization was achieved using a fully implicit scheme. The model was verified against established, peer‐reviewed works, and the assumption that the immiscible fluids (non‐aqueous phase liquids) will flow preferentially through the secondary porosity features in soil was validated. Copyright © 2011 John Wiley & Sons, Ltd.     

Density function evaluation from borehole gravity meter data – regularized spectral domain deconvolution approach

We present a new method of transforming borehole gravity meter data into vertical density logs. This new method is based on the regularized spectral domain deconvolution of density functions. It is a novel alternative to the “classical” approach, which is very sensitive to noise, especially for high‐definition surveys with relatively small sampling steps. The proposed approach responds well to vertical changes of density described by linear and polynomial functions. The model used is a vertical cylinder with large outer radius (flat circular plate) crossed by a synthetic vertical borehole profile. The task is formulated as a minimization problem, and the result is a low‐pass filter (controlled by a regularization parameter) in the spectral domain. This regularized approach is tested on synthetic datasets with noise and gives much more stable solutions than the classical approach based on the infinite Bouguer slab approximation. Next, the tests on real‐world datasets are presented. The properties and presented results make our proposed approach a viable alternative to the other processing methods of borehole gravity meter data based on horizontally layered formations.     

Geochemistry of a confirmed Precambrian impact ejecta deposit: The Grænsesø spherule layer,South Greenland

Samples from a single outcrop of the Graenseso spherule layer, Midternaes, South Greenland, consist of a spherule‐bearing dolomixtite with matrix‐supported intraclasts up to 1 m in size. In addition to field observations, we performed mineralogical and whole rock geochemical analysis, including electron microprobe, neutron activation analysis, X‐ray fluorescence, and mass spectrometry of the horizon and the overlying and underlying strata. We show that the spherules are petrographically similar to those in the Zaonega spherule layer, Karelia, Russia. Our petrographic and chemical results are consistent with the previous suggestion that the Grænsesø layer correlates with the Zaonega layer, and it is possible that both layers are related to the Vredefort impact event. The samples containing spherules, as well as the overlying rocks, have elevated REEs compared to the underlying pre‐impact layer, suggestive of a new continental source of sediment that may be coincident with the impact event. Zircons separated from the lower part of the Grænsesø spherule layer display complex age patterns suggesting that they have genetically different origins based on distinctly different Th/U ratios. Crystallization ages of all groups are ≥ 2.8 Ga, with ~2.8 Ga representing a time of major crustal growth globally. Therefore, we cannot conclusively determine in this study if the zircons are locally derived or if they are transported with the ejecta. The spherule layer was deposited by a high‐energy, subaqueous debris flow, an origin that is consistent with the hypothesis that the layer was deposited by impact‐induced waves and/or currents.