Spatial and temporal variation of natural recharge in the semi-arid valley of Aguascalientes,Mexico

The investigation involves a temporal and spatial analysis to characterize recharge in the Aguascalientes valley in central Mexico. The results are based on a two-part methodological strategy using a geospatial and numerical analysis. Results from the geospatial analysis are based on the analytical hierarchical process (AHP) method, which involves developing a zoning map that qualifies the conditions for groundwater recharge in the valley. In addition, one-dimensional numerical modeling based on the van Genuchten equation was applied for various soil column configurations to estimate the transit time of recharge through the vadose zone. The analysis was performed over a 20-year period from 1995 to 2015. The results of the geospatial analysis indicate that the optimal area for groundwater recharge is the lower zone of the valley, particularly the portion influenced by surface discontinuities. The shortest estimated period of recharge for water to pass through the vadose zone is approximately 4.25 years.     

Identification of arid phases during the last 50 cal. ka BP from the Fuentillejo maar‐lacustrine record (Campo de Calatrava Volcanic Field,Spain)

Geochemical (element analysis, molecular analysis of organic compounds), physical, palynological, mineralogical and sedimentary facies analysis were performed to characterise the sedimentary record in Fuentillejo maar‐lake in the Central Spanish Volcanic Field of Campo de Calatrava, in order to reconstruct the palaeoenvironmental and palaeoclimatic processes which controlled vegetation patterns and deposition of different sedimentary facies. The upper 20 m of core FUENT‐1 show variations in clastic input, water chemistry, vegetation and organic fraction sources in the lake throughout the Late Pleistocene and Holocene. The temporal framework provided by ¹⁴C accelerator mass spectrometry dating allows assigning the sequence to the last 50 cal. ka BP. Arid phases identified in the FUENT‐1 sequence are correlated to Heinrich events (HE) and to stadials of the Dansgaard/Oeschger (D/O) cycles. Siliciclastic facies with high magnetic susceptibility values, high Juniperus pollen content, a low Paq index (aquatic macrophysics proxy index), a decrease in the relative percentage of the n‐C₂₇ and an increase in the n‐C₃₁ alkanes are indicative of arid and colder climatic events related to HE 2, HE 1 and the Younger Dryas (YD). Similar short cold and arid phases during the Holocene were identified at 9.2–8.6, 7.5–7 and 5.5–5 cal. ka BP. In dolomite–mud facies, the pollen data show an increase in the herbs component, mainly – Chenopodiaceae, Artemisia and Ephedra – steppe taxa; a low Paq index, a decrease in the relative percentage of the n‐C₂₇ alkane and an increase in the n‐C₃₁ alkane are also observed. This facies was probably the result of lower lake levels and more saline–alkaline conditions, which can be interpreted as linked to arid–warm periods. These warm and arid phases were more frequent during Marine Isotope Stage (MIS) 3 and the interstadials of MIS 2. HE 4, HE 2, HE 1 and the YD in core FUENT‐1 were immediately followed by increases of warm steppe pollen assemblages that document rapid warming similar to the D/O cycles but do not imply increasing humidity in the area. Fuentillejo hydrology is controlled by changes in the atmospheric and oceanic systems that operated on the North Atlantic region at millennial scale during the last 50 cal. ka BP. Copyright © 2009 John Wiley & Sons, Ltd.     

Analyzing state-dependent model–data comparison in multi-regime systems

An approach to analyze regime change in spatial time series datasets is followed and extended to jointly analyze a dynamical model depicting regime shift and observational data informing the same process. We analyze changes in the joint model–data regime and covariability within each regime. The method is applied to two observational datasets of equatorial sea surface temperature (TAO/TRITON array and satellite) and compared with the predicted data by the ECCO-JPL modeling system.     

Tracing salinization processes in coastal aquifers using an isotopic and geochemical approach: comparative studies in western Morocco and southwest Portugal

Environmental stable and radioactive isotopes (δ²H, δ¹³C, δ¹⁸O; ³H and ¹⁴C), together with physical and geochemical data, were used in the determination of the origins of groundwater salinization and geochemical evolution processes in coastal regions. Two case studies on the Atlantic Coast are discussed, one located in the Essaouira sedimentary basin, western Morocco, and the second, in the Lower Tagus–Sado sedimentary basin, southwest Portugal. In both regions, groundwater degradation occurs by salinization increase to different concentrations and in relation to different origins. The main quality issues for the groundwater resources are related to seawater intrusion, dissolution of diapiric structures intruding the aquifer layers, brine dissolution at depth, and/or evaporation of irrigation water. Anthropogenic pollution ascribed to agricultural activities is another source for groundwater degradation, affecting mainly the shallow aquifers. The apparent ¹⁴C age of the analysed samples ranges from 2.9?±?0.3 up to 45.6?±?0.6 pmC in the Miocene groundwater samples from the basin in Portugal; at the Essaouira basin in Morocco, the ¹⁴C content varies from 60 to 86 pmC. In most of the water samples, the ³H concentration is below the detection limit. In both basins, the isotopic results together with the geochemical data provided an effective label for tracing the mineralization origin and groundwater degradation processes. Further, the isotopic signatures were used in the identification of a paleoclimate (colder period), recorded in the stable isotopic composition and corroborated with the ¹⁴C data.     

Numerical modeling of surface runoff and erosion due to moving rainstorms at the drainage basin scale

A physically-based distributed erosion model (MEFIDIS) was applied to evaluate the consequences of storm movement on runoff and erosion from the Alenquer basin in Portugal. Controlled soil flume laboratory experiments were also used to test the model. Nine synthetic circular storms were used, combining three storm diameters (0.5, 1 and 2 times the Alenquer basin’s axial length) with three speeds of storm movement (0.5, 1 and 2 m/s); storm intensities were synthesized in order to maintain a constant rainfall depth of 50 mm. The model was applied to storms moving downstream as well as upstream along the basin’s axis. In all tests, downstream-moving storms caused significantly higher peak runoff (56.5%) and net erosion (9.1%) than did upstream-moving storms. The consequences for peak runoff were amplified as the storm intensity increased. The hydrograph shapes were also different: for downstream-moving storms, runoff started later and the rising limb was steeper, whereas for upstream moving storms, runoff started early and the rising limb was less steep. Both laboratory and model simulations on the Alenquer basin showed that the direction of storm movement, especially in case of extreme rainfall events, significantly affected runoff and soil loss.     

Carbon isotopes and geochemical processes in CO2-rich cold mineral water,N-Portugal

This paper summarizes a new outlook on the conceptual model of Melgaço–Messegães CO₂-rich cold (≈18 °C) mineral water systems, issuing in N of Portugal, based on their isotopic (²H, ³H, ¹³C, ¹⁴C and ¹⁸O) and geochemical features. Stable isotopes indicate the meteoric origin of these CO₂-rich mineral waters. Based on the isotopic fractionation with the altitude, a recharge altitude between 513 up to 740 m a.s.l. was estimated, corroborating the tritium results. The lowest ³H content (0 TU) is found in the groundwater samples with the highest mineralization. The mineral waters circulation are mainly related to a granitic and granodioritic environment inducing two different groundwater types (Ca/Na–HCO₃ and Na/Ca–HCO₃), indicating different underground flow paths. Calcium dissolution is controlled by hydrolysis of rock-matrix silicate minerals (e.g. Ca-plagioclases) and not associated to anthropogenic sources. The shallow dilute groundwaters exhibit signatures of anthropogenic origins (e.g. NO₃) and higher Na/Ca ratios. The stable isotopes together with the geochemistry provided no indication of mixing between the regional shallow cold dilute groundwater and mineral water systems. The heavy isotopic signatures identified in the δ¹³C data (δ¹³C = 4.7 ‰, performed on the total dissolved inorganic carbon (TDIC) of CO₂-rich mineral waters) could be derived from a deep-seated (upper mantle) source or associated to methanogenesis (CH₄ source). The negligible ¹⁴C content (≈2 pmC) determined in the TDIC of the mineral waters, corroborates the hypothesis of a mantle-derived carbon source to the mineral groundwater systems or dissolution of carbonate layers at depth.     

Geostatistical seismic inversion for non‐stationary patterns using direct sequential simulation and co‐simulation

Geostatistical seismic inversion methods are routinely used in reservoir characterisation studies because of their potential to infer the spatial distribution of the petro‐elastic properties of interest (e.g., density, elastic, and acoustic impedance) along with the associated spatial uncertainty. Within the geostatistical seismic inversion framework, the retrieved inverse elastic models are conditioned by a global probability distribution function and a global spatial continuity model as estimated from the available well‐log data for the entire inversion grid. However, the spatial distribution of the real subsurface elastic properties is complex, heterogeneous, and, in many cases, non‐stationary since they directly depend on the subsurface geology, i.e., the spatial distribution of the facies of interest. In these complex geological settings, the application of a single distribution function and a spatial continuity model is not enough to properly model the natural variability of the elastic properties of interest. In this study, we propose a three‐dimensional geostatistical inversion technique that is able to incorporate the reservoir's heterogeneities. This method uses a traditional geostatistical seismic inversion conditioned by local multi‐distribution functions and spatial continuity models under non‐stationary conditions. The procedure of the proposed methodology is based on a zonation criterion along the vertical direction of the reservoir grid. Each zone can be defined by conventional seismic interpretation, with the identification of the main seismic units and significant variations of seismic amplitudes. The proposed method was applied to a highly non‐stationary synthetic seismic dataset with different levels of noise. The results of this work clearly show the advantages of the proposed method against conventional geostatistical seismic inversion procedures. It is important to highlight the impact of this technique in terms of higher convergence between real and inverted reflection seismic data and the more realistic approximation towards the real subsurface geology comparing with traditional techniques.     

Out‐of‐plane behaviour of a full scale stone masonry façade. Part 1: specimen and ground motion selection

The out‐of‐plane response of walls in existing stone masonry buildings is one of the major causes of vulnerability commonly observed in post‐earthquake damage surveys. In this context, a shaking table (ST) test campaign was carried out on a full‐scale masonry façade mainly focusing on the characterization of its out‐of‐plane overturning behaviour. The structure tested on the ST is a partial reproduction of an existing building from Azores, damaged during the 9 July 1998 Faial earthquake. The definition of the tested specimen as well as the selection of the input ground motion is reported in this paper. A specific emphasis is given to the definition of the time‐history to be applied during the tests because it was felt as an essential and crucial part of the work to obtain the desired overturning behaviour. The accelerogram to be imposed was selected from a large set of accelerograms (74) by means of a step‐by‐step procedure on the basis of several numerical analyses resorting to the rocking response of rigid blocks. A companion paper (Part 2) focuses on the ST test results and detailed data interpretation. Copyright © 2013 John Wiley & Sons, Ltd.     

Estimating quality factor from surface seismic data: A comparison of current approaches

The performances of the spectral ratio (SR), frequency centroid shift (FCS), and frequency peak shift (FPS) methods to estimate the effective quality factor Q are compared. These methods do not demand true amplitude data and their implementations were done following an “as simple as possible” approach to highlight their intrinsic potentials and limitations. We use synthetic zero-offset seismic data generated with a simple layer-cake isotropic model. The methods can be ranked from simple to complex in terms of automation as: FPS, FCS and SR. This is a consequence of: (i) peak identification consists basically of a sorting procedure, (ii) centroid estimation involves basically the evaluation of two well-behaved integrals, and (iii) implementation of the SR method involves at least choosing a usable frequency bandwidth and fitting a gradient. The methods can be ranked from robust to sensitive in the presence of noise content in the sequence SR, FCS, and FPS. This is consequence of: (i) the gradient estimate associated to the SR method averages out the noise content in the entire usable frequency bandwidth, (ii) in the presence of moderate-to-high noise level, the centroid estimation is biassed towards overestimating Q due to noise contribution in the tail of the amplitude spectrum, and (iii) peak identification is unstable due to local noise fluctuation in the amplitude spectrum around the peak frequency. Regarding the stability of the estimates relative to the attenuation amount, SR and FCS methods show similar behaviours, whereas FPS method presents an inferior performance. This fact is an indirect consequence of the sensitivity of FPS method to the noise content because the higher is the attenuation the lower is the signal-to-noise ratio. Finally, regarding the robustness of the methods to the presence of dipping layers, only SR and FCS methods provide good estimates, at least to typical dips in non-faulted sedimentary layers, with the estimates obtained with SR method being more accurate that those obtained with FCS method. Except in relation to the automation complexity, which is less important than the performances of the methods, SR method was superior or showed similar performance to FCS method in all scenarios we tried.