Hydrogeology and sustainable future groundwater abstraction from the Agua Verde aquifer in the Atacama Desert,northern Chile

The hyper-arid conditions prevailing in Agua Verde aquifer in northern Chile make this system the most important water source for nearby towns and mining industries. Due to the growing demand for water in this region, recharge is investigated along with the impact of intense pumping activity in this aquifer. A conceptual model of the hydrogeological system is developed and implemented into a two-dimensional groundwater-flow numerical model. To assess the impact of climate change and groundwater extraction, several scenarios are simulated considering variations in both aquifer recharge and withdrawals. The estimated average groundwater lateral recharge from Precordillera (pre-mountain range) is about 4,482 m³/day. The scenarios that consider an increase of water withdrawal show a non-sustainable groundwater consumption leading to an over-exploitation of the resource, because the outflows surpasses inflows, causing storage depletion. The greater the depletion, the larger the impact of recharge reduction caused by the considered future climate change. This result indicates that the combined effects of such factors may have a severe impact on groundwater availability as found in other groundwater-dependent regions located in arid environments. Furthermore, the scenarios that consider a reduction of the extraction flow rate show that it may be possible to partially alleviate the damage already caused to the aquifer by the continuous extractions since 1974, and it can partially counteract climate change impacts on future groundwater availability caused by a decrease in precipitation (and so in recharge), if the desalination plant in Taltal increases its capacity.     

U–Pb evidence for a polyorogenic evolution of the HP–HT units of the NW Iberian Massif

. A isotope dilution thermal ionisation mass spectrometry U-Pb geochronological study was carried out on the high-pressure and high-temperature units (HP-HT units) overlying the oceanic suture in the Allochthonous Complexes of the NW Iberian Variscan Belt. The rocks investigated are seven granulite- to eclogite-facies paragneisses and one leucosome within mafic high-pressure granulites in the Ordenes and Cabo Ortegal Complexes of NW Spain. U-Pb dating of zircon, monazite, titanite and rutile reveal the presence of a pervasive Early Ordovician metamorphic event at ca. 500-480 Ma and a later Early Devonian event at ca. 400-380 Ma. The U-Pb ages, in conjunction with petrological and structural data, indicate that the high-pressure event recorded by these rocks is Early Ordovician in age. Monazite ages in the paragneisses suggest that peak metamorphic conditions were reached at ca. 500-485 Ma. Subsequently, the rock ensemble underwent exhumation accompanied by partial melting and zircon growth at ca. 485-470 Ma. Melting of mafic granulites was coeval with this latter episode as indicated by zircon crystallisation age in the leucosomes dated at ca. 486 Ma. Based on these data and on the general features of magmatism and metamorphic evolution, it is proposed that this process took place at a convergent plate boundary within a peri-Gondwanan oceanic domain. Monazite, titanite and rutile data in some of the samples studied show evidence of a second metamorphic episode that took place between ca. 400 and 380 Ma (with a peak at ca. 390-385 Ma). This Early Devonian event, at variance with the previous one, was not pervasive, but, rather, was localised in areas of intense Variscan tectonothermal reworking. It is claimed that this later metamorphic event was recorded by the U-Pb system in areas where monazite and titanite growth was enhanced by fluid circulation in highly strained rocks (Variscan shear zones). According to previous structural studies and Ar-Ar dating of fabrics, this Early Devonian episode took place as the HP-HT units were deformed and thrusted upon the ophiolitic units in the early stages of the Variscan collision.     

From brightening to dimming in sunshine duration over the eastern and central Tibetan Plateau (1961–2005)

The Tibetan Plateau (TP) with an average elevation of over 4,000 m asl is the highest and most extensive highland in the world. We used monthly mean sunshine duration from the Chinese Meteorological Administration to examine the spatial and temporal variability of sunshine duration at 71 stations with elevations above 2,000 m asl in the eastern and central TP during the 1961–2005 period. The temporal evolution of the mean annual sunshine duration series shows a significant increase from 1961 to 1982 at a rate of 49.8 h/decade, followed by a decrease from 1983 to 2005 at a rate of ?65.1 h/decade, with an overall significant decrease at a rate of ?20.6 h/decade during the whole 1961–2005 period, which is mainly due to the summer and spring seasons. This confirms the evidence that sunshine duration in the TP ranges from brightening to dimming in accordance with sunshine duration trends in the rest of China. The surface solar radiation downwards from ERA-40 reanalysis data in the same region confirms the brightening/dimming phenomenon shown by the sunshine duration before/after the 1980s. Otherwise, additional climatic variables such as low cloud amount, total cloud amount, precipitation, relative humidity and water vapor pressure, in most cases, exhibit significant negative correlation with sunshine duration in the TP on an annual and seasonal basis before and after 1982, respectively. The trends of these variables suggest that changes in some of them might be related to the brightening and dimming detected with the use of sunshine duration measurements over the TP. We also hypothesize that the impact of anthropogenic aerosols upon the climatic variables analyzed cannot be rejected, especially in the significant increase in low cloud cover since approximately 1980.     

Trends in the frequency of synoptic types in central-southern Chile in the period 1961–2012 using the Jenkinson and Collison synoptic classification

Atmospheric circulation patterns in southern Chile (42° 30′ S) were studied in order to determine and analyse the most characteristic synoptic types and their recent trends, as well as to gain an understanding of how they are associated with low-frequency variability patterns. According to the Jenkinson and Collison (J&C) classification method, a 16-point grid of sea-level pressure data was employed. The findings reveal that some synoptic types show statistically significant trends with a 95% confidence level, positively for anticyclonic westerly hybrids (AW) and advective types for third and fourth quadrant wind flows (W, NW, and N) and negatively for SW and cyclonic hybrids (CS and CSW). A model has been constructed of the linear regression of some weather types with teleconnections that most affect Chile: the undetermined types (U), AW were associated with El Niño or the warm phase of the Pacific Decadal Oscillation (PDO), whereas the cyclonic northerly and cyclonic northeasterly types (CN and CNE) were associated with La Niña or cool phase of the PDO. The weather types associated with Antarctic Oscillation (AAO) in its positive phase are anticyclonic northerly and northeasterly and northerly advection types, while in its negative phase are cyclonic southwesterly and advection types.     

Long-lasting floods buffer the thermal regime of the Pampas

The presence of large water masses influences the thermal regime of nearby land shaping the local climate of coastal areas by the ocean or large continental lakes. Large surface water bodies have an ephemeral nature in the vast sedimentary plains of the Pampas (Argentina) where non-flooded periods alternate with flooding cycles covering up to one third of the landscape for several months. Based on temperature records from 17 sites located 1 to 700 km away from the Atlantic coast and MODIS land surface temperature data, we explore the effects of floods on diurnal and seasonal thermal ranges as well as temperature extremes. In non-flooded periods, there is a linear increase of mean diurnal thermal range (DTR) from the coast towards the interior of the region (DTR increasing from 10 to 16 K, 0.79 K/100 km, r ² = 0.81). This relationship weakens during flood episodes when the DTR of flood-prone inland locations shows a decline of 2 to 4 K, depending on surface water coverage in the surrounding area. DTR even approaches typical coastal values 500 km away from the ocean in the most flooded location that we studied during the three flooding cycles recorded in the study period. Frosts-free periods, a key driver of the phenology of both natural and cultivated ecosystems, are extended by up to 55 days during floods, most likely as a result of enhanced ground heat storage across the landscape (~2.7 fold change in day-night heat transfer) combined with other effects on the surface energy balance such as greater night evaporation rates. The reduced thermal range and longer frost-free periods affect plant growth development and may offer an opportunity for longer crop growing periods, which may not only contribute to partially compensating for regional production losses caused by floods, but also open avenues for flood mitigation through higher plant evapotranspirative water losses.     

Coupled thermo-hydro-bio-geochemical reactive transport model of CERBERUS heating and radiation experiment in Boom clay

Final disposal of high-level radioactive waste in deep repositories in clay formations is being considered by several countries. Repository safety assessment requires the use of numerical models of groundwater flow, solute transport and chemical processes. These models are being developed from data and knowledge gained from in situ experiments such as the CERBERUS experiment carried out at the HADES facility excavated in the Boom clay formation at Mol (Belgium). This long-term experiment is aimed at evaluating the effect of heating and radiation in Boom clay. The test was performed in a cased well drilled at 223 m depth and lasted from 1989 to 1994. A ⁶⁰Co source of 400 TBq and two heaters were emplaced inside the well. Dose rate, temperature, porewater pressure and pH/Eh were measured in situ during the experiment and gas and porewater samples were taken for chemical analyses. Here a coupled thermo-hydro-geochemical (THC) model of the CERBERUS experiment is presented which accounts simultaneously for heating, radiation, solute diffusion and a suite of geochemical reactions including: aqueous complexation, acid–base, redox, mineral dissolution/precipitation, cation exchange and gas dissolution/ex-solution. Computed results indicate that heating and radiation causes a slight oxidation, a decrease in pH, slight changes in porewater chemistry and pyrite dissolution near the well. THC model results follow the general evolution of chemical data, but cannot fit SO₄ data. Model discrepancies are partly overcome when microbially-mediated Fe and SO₄ reduction are taken into account in a coupled thermo-hydro-bio-geochemical (THBC) model. This THBC model captures the trends of geochemical data, improves the fit to dissolved SO₄ and predicts pyrite precipitation, a process observed near the CERBERUS well. The ability of the THBC numerical model to reproduce the overall trends of geochemical data of the CERBERUS experiment provides confidence in such a model as a suitable tool for the long-term prediction of geochemistry in the near field of a HLW repository in clay. However, the small number of available chemical data throughout the experiment and the lack of DOC and microbial data allow only a partial validation of the THBC model.     

Using CORINE land cover and the point survey LUCAS for area estimation

CORINE land cover 2000 (CLC2000) is a European land cover map produced by photo-interpretation of Landsat ETM+ images. Its direct use for area estimation can be strongly biased and does not generally report single crops. CLC areas need to be calibrated to give acceptable statistical results.LUCAS (land use/cover area frame survey) is a point survey carried out in 2001 and 2003 in the European Union (EU15) on a systematic sample of clusters of points. LUCAS is especially useful for area estimation in geographic units that do not coincide with administrative regions, such as set of coastal areas defined with a 10 km buffer. Some variance estimation issues with systematic sampling of clusters are analysed.The contingency table obtained overlaying CLC and LUCAS gives the fine scale composition of CLC classes. Using CLC for post-stratification of LUCAS is equivalent to the direct calibration estimator when the sampling units are points. Stratification is easier to adapt to a scheme in which the sampling units are the clusters of points used in LUCAS 2001/2003.     

Testing spatial heterogeneity in geographically weighted principal components analysis

We propose a method to evaluate the existence of spatial variability in the covariance structure in a geographically weighted principal components analysis (GWPCA). The method, that is extensive to locally weighted principal components analysis, is based on performing a statistical hypothesis test using the eigenvectors of the PCA scores covariance matrix. The application of the method to simulated data shows that it has a greater statistical power than the current statistical test that uses the eigenvalues of the raw data covariance matrix. Finally, the method was applied to a real problem whose objective is to find spatial distribution patterns in a set of soil pollutants. The results show the utility of GWPCA versus PCA.     

Evidence of episodic long-lived eruptions in the Yuma,Ginsburg, Jesús Baraza and Tasyo mud volcanoes,Gulf of Cádiz

High-resolution single channel and multichannel seismic reflection profiles and multibeam bathymetric and backscatter data collected during several cruises over the period 1999 to 2007 have enabled characterising not only the seabed morphology but also the subsurface structural elements of the Yuma, Ginsburg, Jesús Baraza and Tasyo mud volcanoes (MVs) in the Gulf of Cádiz at 1,050–1,250 m water depth. These MVs vary strongly in morphology and size. The data reveal elongated cone-shaped edifices, rimmed depressions, and scarps interpreted as flank failures developed by collapse, faulting, compaction and gravitational processes. MV architecture is characterised by both extrusive and intrusive complexes, comprising stacked edifices (including seabed cones and up to four buried bicones) underlain by chaotic vertical zones and downward-tapering cones suggesting feeder systems. These intrusive structures represent the upper layer of the feeder system linking the fluid mud sources with the constructional edifices. The overall architecture is interpreted to be the result of successive events of mud extrusion and outbuilding alternating with periods of dormancy. Each mud extrusion phase is connected with the development of an edifice, represented by a seabed cone or a buried bicone. In all four MVs, the stacked edifices and the intrusive complexes penetrate Late Miocene–Quaternary units and are rooted in the Gulf of Cádiz wedge emplaced during the late Tortonian. Major phases of mud extrusion and outbuilding took place since the Late Pliocene, even though in the Yuma and Jesús Baraza MVs mud volcanism started in the Late Miocene shortly after the emplacement of the Gulf of Cádiz wedge. This study shows that fluid venting in the eastern sector of the Gulf of Cádiz promoted the outbuilding of large long-lived mud volcanoes active since the Late Miocene, and which have been reactivated repeatedly until recent times.