Lithological control of land subsidence induced by groundwater withdrawal in new urban AREAS (Granada Basin,SE Spain). Multiband DInSAR monitoring

Ground subsidence in the southeastern border of the Granada Basin (SE Spain) has been studied using remote sensing techniques. Over the last decades, the region has experienced a huge urban expansion, which has caused a substantial increase in water supply requirements. Water needs are exclusively met by groundwater by means of numerous pumping wells, which exploit a confined detrital aquifer of alluvial fan deposits with a heterogeneous facies distribution. A general piezometric level decline (up to 50 m) has been recorded in the aquifer during the past 30 years that has induced the generation of a subsiding area with oval shape oriented WNW‐ESE just where the new urban areas and pumping wells are located. Subsidence has been monitored by exploiting synthetic aperture radar (SAR) images from ENVISAT (2003–2009) and Cosmo‐SkyMed (2011–2014). A new approach, which combines A‐DInSAR and small‐area persistent scatterer interferometry (PSI) analysis, has been applied obtaining a good accuracy regarding temporal and spatial dimension of the subsidence. ENVISAT data (2003–2009) reveal subsidence rates up to 10–15 mm/year, and Cosmo‐SkyMed (2011–2014) values slightly lower; up to 10 mm/year. Temporal variations in the subsidence velocity are in accordance with the rainfall pattern and piezometric fluctuations in the aquifer. The sector with highest rates of subsidence does not correspond to the area with more intense groundwater exploitation but to the area with greater presence of clays in the confining layer of the aquifer. There is a clear lithological control in the spatial distribution of the ground subsidence. This work integrates detailed geological and hydrogeological data with differential SAR interferometry monitoring with the aim to better understand subsidence processes in detrital aquifers with small‐scale heterogeneity. Copyright © 2016 John Wiley & Sons, Ltd.     

Modelling of coupled fluid‐mechanical problems in fractured geological media using enriched finite elements

Geological environments, such as petroleum reservoirs, normally exhibit physical discontinuities, for example, fractures and faults. Because of the reduced thickness of these discontinuities, finite element formulations with strong discontinuity have been applied to the numerical modelling of geological environments. Until now, two relevant characteristics of petroleum reservoirs have not been addressed by these formulations. The first is the pore pressure jump in the direction normal to a discontinuity in a fluid‐mechanical coupling condition, which is present primarily in sealing faults owing to the contrast of permeability with the porous medium. The absence of this jump can affect the prediction of the deformability of a physical discontinuity. Furthermore, reservoir models frequently use coarse meshes. Thus, the method used to evaluate the pore pressure in the discontinuity may exhibit a strong dependence relative to the mesh refinement. Based on these characteristics, in this study, a formulation of an enriched finite element for application to coupled fluid‐mechanical problems with pre‐existing physical discontinuities saturated by a single fluid is presented. The formulation employs discontinuous interpolation functions and enables the reproduction of jumps of displacement and pore pressure associated with a discontinuity inside the element without the need to discretise it. An approximation to estimate the pore pressure in the discontinuity was developed, one which seeks to minimise the influence of refinement. The element's response is verified by comparison with a one‐dimensional analytical solution and simple examples that are simulated using commercial software. Copyright © 2015 John Wiley & Sons, Ltd.     

Volume change response and fabric evolution of granular MX80 bentonite along different hydro-mechanical stress paths

Despite the increasing understanding of bentonite behaviour, there is still missing evidence on how different hydro-mechanical loadings, including sequences of hydration and compression, affect the fabric and the volume change behaviour of the material. It is generally assumed that the interplay between the behaviour of clay assemblages and the overall fabric of the material is the reason of having final states that are dependent on the stress path followed. Here the results of an experimental campaign aiming to study these factors are reported and discussed. Free swelling and swelling pressure tests were performed, both followed by compression to a relatively high stress. The experimental program involved various samples that were dismantled at intermediate states in order to perform microstructural observations by means of mercury intrusion porosimetry and electronic scanning microscopy. It was observed that while the void ratio at a given stress level depends on the stress path, subsequent compression led to a unique virgin compression line. The data obtained at the microscale gave further insight for an interpretation of the volume change behaviour observed at the macroscale, showing that at high stress the material tends to recover the same fabric regardless of the path to saturation.

     

Short-term effects of the prestige oil spill on the peregrine falcon (Falco peregrinus)

We have monitored the distribution, population status, breeding success, turnover rate and diet of a Peregrine Falcon population in Bizkaia (North of Spain) since 1997. On the 13th November 2002, the tanker Prestige sunk off La Coruña (NW Spain) causing an oil spill that affected the whole of the Cantabrian Coast and the Southwest of France. The total number of birds affected by the Prestige oil spill was expected to be between 115,000 and 230,000, some of them raptors. The loss of clutches during the incubation period increased significantly and was correlated with the loss of females. Moreover, the turnover rate of the population increased from 21% to 30%. The polycyclic aromatic hydrocarbon concentrations in the eggs, collected from five nests after they were deserted, ranged from 21.20 ng/g to 461.08 ng/g, values which are high enough to cause the death of the embryos and poisoning of adult birds. The effects of pollution reached inland since some inland-breeding falcons prey on shorebirds that use rivers during their migratory flights. As the Prestige oil spill has clearly resulted in increased rates of adult mortality and reduced fertility, we suggest that the environmental authorities urgently undertake measures aimed at protecting the Peregrine Falcon in Bizkaia.     

Geo-environmental mapping using physiographic analysis: constraints on the evaluation of land instability and groundwater pollution hazards in the Metropolitan District of Campinas,Brazil

Geo-environmental terrain assessments and territorial zoning are useful tools for the formulation and implementation of environmental management instruments (including policy-making, planning, and enforcement of statutory regulations). They usually involve a set of procedures and techniques for delimitation, characterisation and classification of terrain units. However, terrain assessments and zoning exercises are often costly and time-consuming, particularly when encompassing large areas, which in many cases prevent local agencies in developing countries from properly benefiting from such assessments. In the present paper, a low-cost technique based on the analysis of texture of satellite imagery was used for delimitation of terrain units. The delimited units were further analysed in two test areas situated in Southeast Brazil to provide estimates of land instability and the vulnerability of groundwater to pollution hazards. The implementation incorporated procedures for inferring the influences and potential implications of tectonic fractures and other discontinuities on ground behaviour and local groundwater flow. Terrain attributes such as degree of fracturing, bedrock lithology and weathered materials were explored as indicators of ground properties. The paper also discusses constraints on- and limitations of- the approaches taken.     

An intercomparison of observed and simulated extreme rainfall and temperature events during the last half of the twentieth century: part 2: historical trends

We analyze historical simulations of variability in temperature and rainfall extremes in the twentieth century, as derived from various global models run informing the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC-AR4). On the basis of three indices of climate extremes, we compare observed and modeled trends in time and space, including the direction and significance of the changes at the scale of South America south of 10° S. The climate extremes described warm nights, heavy rainfall amounts and dry spells. The reliability of the GCM simulations is suggested by similarity between observations and simulations in the case of warm nights and extreme rainfall in some regions. For any specific extreme temperature index, minor differences appear in the spatial distribution of the changes across models in some regions, while substantial differences appear in regions in the interior of tropical and subtropical South America. The differences are in the relative magnitude of the trends. Consensus and significance are less strong when regional patterns are considered, with the exception of the La Plata Basin, where observed and simulated trends in warm nights and extreme rainfall are evident.     

Sensitivity of Spatial Analysis Neural Network Training and Interpolation to Structural Parameters

A Spatial Analysis Neural Network (SANN) algorithm was applied for the analysis of geospatial data, on the basis of nonparametric statistical analysis and the concepts of traditional Artificial Neural Networks. SANN consists of a number of layers in which the neurons or nodes between layers are interconnected successively in a feed-forward direction. The Gaussian Kernel Function layer has several nodes, and each node has a transfer or an activation function that only responds (or activates) when the input pattern falls within its receptive field, which is defined by its smoothing parameter or width. The activation widths are functions of the model structural parameters, including the number of the nearest neighbor points P and a control factor F. The estimation method is based on two operational modes, namely, a training-validation mode in which the model structure is constructed and validated, and an interpolation mode. In this paper we discuss the effect of varying F and P upon the accuracy of the estimation in a two-dimensional domain for different input field sizes, using spatial data of wheat crop yield from Eastern Colorado. Crop yield is estimated as a function of the two-dimensional Cartesian coordinates (easting and northing). The results of the research led to the conclusion that optimal values of F and P depend on the sample size, i.e., for small data sets F=1.5 and P=7 while for large data sets F=2.5 and P=9. In addition, the accuracy of the interpolated field varies with the sample size. As expected for small sample sizes, the interpolated field and its variability may be significantly underestimated.     

Estimates of Reynolds stress in a highly energetic shelf sea

In recent years the use of Acoustic Doppler Current Profilers (ADCPs) to estimate Reynolds stresses, using the so-called variance method, has become popular; and although there was great effort in studying the uncertainties on this technique, there were no reports in the main literature of its validity using independent measurements. This work reports on the comparison of ADCP and Acoustic Current Velocimeter (ADV) estimates of Reynolds stresses. The comparison of the ADCP and ADV is encouraging during periods when no strong waves were present with both the explained variance of 0.8 and the slope of the regression being 0.97. Nevertheless, when strong waves are present the method breaks down and the comparison between ADCP and ADV is very poor with R² =0.04.     

U–Pb SHRIMP zircon dating of Grenvillian metamorphism in Western Sierras Pampeanas (Argentina): Correlation with the Arequipa-Antofalla craton and constraints on the extent of the Precordillera Terrane

Metamorphism of Grenvillian age (ca. 1.2 Ga; U–Pb zircon dating) is recognized for the first time in the Western Sierras Pampeanas (Sierra de Maz). Conditions reached granulite facies (ca. 780 °C and ca. 780 MPa). Comparing geochronological and petrological characteristics with other outcrops of Mesoproterozoic basement, particularly in the northern and central Arequipa-Antofalla craton, we suggest that these regions were part of a single continental crustal block from Mesoproterozoic times, and thus autochthonous or parautochthonous to Gondwana.     

Stability, solubility and maximum metal binding capacity in metal–humic complexes involving humic substances extracted from peat and organic compost

The aim of this work is to investigate the influence of pH and the metal:humic substances (HS) ratio on HS complexing capacity and the stability and solubility of metal–HS complexes in solution. We selected four HS with different physicochemical properties and studied their interaction with Cu(II), Zn(II) and Fe(II) at different pH and metal:HS ratios. The selected HS were a humic acid and a whole humic extract (containing the humic and fulvic acids) extracted from black peat, and a fulvic acid and a whole humic extract extracted from a compost of grape solid wastes.Our results showed that HS complexing capacity significantly varied as a function of pH, thus indicating the influence of both functional group ionisation and molecular conformation on this property. As was expected, total acidity affected the complexing capacity of the selected HS.The results related to stability and complexing capacity indicated the possible presence of two binding patterns, one at acid-neutral pH probably involving carboxylates, and another at alkaline pH probably involving carboxylates and phenolic groups. The relationship between these binding patterns and the strength of the binding process varied according to the complexed metal.Complex solubility was greatly affected by the ratio between the concentration of free ionised functional groups and the molecular weight in the HS studied.