Fluid and silicate glass inclusions in ultramafic and mafic xenoliths from Hierro,Canary Islands: implications for mantle metasomatism

Fluid and solid inclusions have been studied in selected samples from a series of spinel-bearing Crdiopside-and Al-augite-series ultramafic (harzburgites, lherzolites, and olivine-clinopyroxene-rich rocks), and gabbroic xenoliths from Hierro, Canary Islands. In these samples several generations of fluid inclusions and ultramafic-and mafic-glass inclusions may be texturally related to different stages of crystal growth. The fluid inclusions consist of pure, or almost pure, CO₂. The solid inclusions in the ultramafic xenoliths comprise early inclusions of devitrified ultramafic glass, sulphide inclusions, as well as polyphase inclusions (spinel+clinopyroxene±glass±other silicates) believed to have formed from trapped basaltic melts. Vitreous basaltic glass±CO₂±sulphide±silicates are common as secondary inclusions in the ultramafic xenoliths, and as primary inclusions in the gabbroic xenoliths. Microthermometry gives minimum trapping temperatures of 1110° C for the early ultramafic-and mafic-glass inclusions, and a maximum of 1260–1280° C for late inclusions of host basaltic glass. In most samples the CO₂ inclusions show a wide range in homogenization temperatures (-40 to +31° C) as a result of decrepitation during ascent. The lowest homogenization temperatures of about-40° C, recorded in some of the smallest CO₂ inclusions, indicate a minimum depth of origin of 35 km (12 kbar) for both the Cr-diopside-and Al-augite-series xenoliths. The gabbroic xenoliths originate from a former magma chamber at a depth of 6–12 km.Contribution no. 100 of the Norwegian programme of the International Lithosphere Project     

Global and local sensitivity analyses for coupled thermo–hydro–mechanical problems

This paper focuses on the sensitivity analysis for coupled thermo–hydro–mechanical problems employing both local and global sensitivity methods. A derivative‐based method is used in the local sensitivity approach, whereas the random balance designs method is used for the global sensitivity analysis. The main goal is to investigate the effect of uncertainties in the constitutive parameters on the results from nonlinear coupled thermo–hydro–mechanical analyses of unsaturated soil behavior whose modeling generally involves large sets of constitutive relations. Knowing the parameter sensitivity allows to qualitatively assess the validity of the results obtained by computational simulations of high‐risk situations, for example, emerging nuclear waste repositories. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

Stability of massive ground ice bodies in University Valley,McMurdo Dry Valleys of Antarctica: Using stable O–H isotope as tracers of sublimation in hyper-arid regions

To date, studies of the stability of subsurface ice in the McMurdo Dry Valleys of Antarctica have been mainly based on climate-based vapor diffusion models. In University Valley (1800 m), a small glacier is found at the base of the head of the valley, and adjacent to the glacier, a buried body of massive ice was uncovered beneath 20–40 cm of loose cryotic sediments and sandstone boulders. This study assesses the origin and stability of the buried body of massive ice by measuring the geochemistry and stable O–H isotope composition of the ice and applies a sublimation and molecular diffusion model that accounts for the observed trends. The results indicate that the buried massive ice body represents an extension of the adjacent glacier that was buried by a rock avalanche during a cold climate period. The contrasting δ¹⁸O profiles and regression slope values between the uppermost 6 cm of the buried massive ice (upward convex δ¹⁸O profile and S_D-18O = 5.1) and that below it (progressive increase in δ¹⁸O and S_D-18O = 6.4) suggest independent post-depositional processes affected the isotope composition of the ice. The upward convex δ¹⁸O profile in the uppermost 6 cm is consistent with the ice undergoing sublimation. Using a sublimation and molecular diffusion model, and assuming that diffusion occurred through solid ice, the sublimation rate needed to fit the measured δ¹⁸O profile is 0.2 ? 10^? 3 mm yr^? 1, a value that is more similar to net ice removal rates derived from ³He data from cobbles in Beacon Valley till (7.0 ? 10^? 3 mm yr^? 1) than sublimation rates computed based on current climate (0.1–0.2 mm yr^?1). We suggest that the climate-based sublimation rates are offset due to potential ice recharge mechanisms or to missing parameters, particularly the nature and thermo-physical properties of the overlying sediments (i.e., temperature, humidity, pore structure and ice content, grain size).     

Multi-gas Emissions Pathways to Meet Climate Targets

So far, climate change mitigation pathways focus mostly on CO₂ and a limited number of climate targets. Comprehensive studies of emission implications have been hindered by the absence of a flexible method to generate multi-gas emissions pathways, user-definable in shape and the climate target. The presented method ‘Equal Quantile Walk’ (EQW) is intended to fill this gap, building upon and complementing existing multi-gas emission scenarios. The EQW method generates new mitigation pathways by ‘walking along equal quantile paths’ of the emission distributions derived from existing multi-gas IPCC baseline and stabilization scenarios. Considered emissions include those of CO₂ and all other major radiative forcing agents (greenhouse gases, ozone precursors and sulphur aerosols). Sample EQW pathways are derived for stabilization at 350 ppm to 750 ppm CO₂ concentrations and compared to WRE profiles. Furthermore, the ability of the method to analyze emission implications in a probabilistic multi-gas framework is demonstrated. The probability of overshooting a 2 ^∘C climate target is derived by using different sets of EQW radiative forcing peaking pathways. If the probability shall not be increased above 30%, it seems necessary to peak CO₂ equivalence concentrations around 475 ppm and return to lower levels after peaking (below 400 ppm). EQW emissions pathways can be applied in studies relating to Article 2 of the UNFCCC, for the analysis of climate impacts, adaptation and emission control implications associated with certain climate targets. See for EQW-software and data.     

River–aquifer interactions in a semi‐arid environment stressed by groundwater abstraction

Rivers and aquifers are, in many cases, a connected resource and as such the interactions between them need to be understood and quantified for the resource to be managed appropriately. The objective of this paper is to advance the understanding of river–aquifer interactions processes in semi‐arid environments stressed by groundwater abstraction. This is performed using data from a specific catchment where records of precipitation, evapotranspiration, river flow, groundwater levels and groundwater abstraction are analysed using basic statistics, hydrograph analysis and a simple mathematical model to determine the processes causing the spatial and temporal changes in river–aquifer interactions. This combined approach provides a novel but simple methodology to analyse river–aquifer interactions, which can be applied to catchments worldwide. The analysis revealed that the groundwater levels have declined (~ 3 m) since the onset of groundwater abstraction. The decline is predominantly due to the abstraction rather than climatic changes (r = 0.84 for the relationship between groundwater abstraction and groundwater levels; r = 0.92 for the relationship between decline in groundwater levels and magnitude of seasonal drawdown). It is then demonstrated that, since the onset of abstraction, the river has changed from being gaining to losing during low‐flow periods, defined as periods with flow less than 0.5, 1.0 or 1.5 GL/day (1 GL/day = 1 × 10⁶ m³/day). If defined as < 1.0 GL/day, low‐flow periods constitute approximately 65% of the river flows; the periods where the river is losing at low‐flow conditions are thus significant. Importantly, there was a significant delay (> 10 years) between the onset of groundwater abstraction and the changeover from gaining to losing conditions. Finally, a relationship between the groundwater gradient towards the river and the river flow at low‐flow is demonstrated. The results have important implications for water management as well as water ecology and quality. Copyright © 2012 John Wiley & Sons, Ltd.     

CO2 sequestration and extreme Mg depletion in serpentinized peridotite clasts from the Devonian Solund basin, SW-Norway

The conglomerates of the Solund Devonian basin of SW-Norway contain numerous (locally up to 20 vol.%) peridotitic clasts with concentric mm- to 10-cm thick zones of varying red to black color. The peridotite clasts show a clear, alteration-related textural evolution. The least-altered rocks are partly serpentinized peridotites, showing a typical mesh texture with veins of serpentine, magnesite and Ni-rich magnetite surrounding olivine (Fo₉₁) relicts and its Mg-depleted, clay-like alteration product (deweylite assemblage). In the more advanced ophicarbonate stage, the mesh cells contain calcite, silica and are surrounded by talc. In the final stage, quartz, calcite, and hematite dominate the mineralogy and occur together with minor amounts of chromite, talc, Cr-chlorite, and Cr-hydroandradite. In tandem with this textural evolution is a decrease in MgO from 40 to 2.5 wt% and a CaO increase from 1 to 35 wt%. All peridotite clasts are characterized by high Cr and Ni concentrations. The chemistry and the textural evolution show that the clasts formed by an extreme Mg-mobilization from the peridotite, with development of secondary porosity and subsequent precipitation of calcite. MgO removed from the clasts after burial is in part consumed by replacement reactions in the sediment matrix around the clasts where Mg-free minerals (e.g., almandine) are replaced by Mg-bearing minerals (e.g., talc). Calculated apparent ⁸⁷Sr/⁸⁶Sr ratios of the clasts at 385 Ma (0.7124-0.7139), corresponding to the inferred age of sediment deposition and incipient clast alteration, indicate interaction with diagenetic basinal fluids. We explain the reaction history as a three stage process involving (a) partial serpentinization of olivine in an oceanic environment (b) breakdown of olivine relicts to the deweylite assemblage resulting in mobilization of MgO under (near-) surface conditions in a tropical Devonian climate and (c) further Mg-mobilization and replacement of the deweylite assemblage by calcite and quartz after diagenesis. Sedimentary basins with abundant weathered peridotite represent potential sites for a permanent CO₂ storage by formation of calcite in a low-temperature environment.     

The influence of regional circulation patterns on wet and dry mineral dust and sea salt deposition over Greenland

Annually resolved ice core records from different regions over the Greenland ice sheet (GrIS) are used to investigate the spatial and temporal variability of calcium (Ca²⁺, mainly from mineral dust) and sodium (Na⁺, mainly from sea salt) deposition. Cores of high common inter-annual variability are grouped with an EOF analysis, resulting in regionally representative Ca²⁺ and Na⁺ records for northeastern and central Greenland. Utilizing a regression and validation method with ERA-40 reanalysis data, these common records are associated with distinct regional atmospheric circulation patterns over the North American Arctic, Greenland, and Central to Northern Europe. These patterns are interpreted in terms of transport and deposition of the impurities. In the northeastern part of the GrIS sea salt records reflect the intrusion of marine air masses from southeasterly flow. A large fraction of the Ca²⁺ variability in this region is connected to a circulation pattern suggesting transport from the west and dry deposition. This pattern is consistent with the current understanding of a predominantly Asian source of the dust deposited over the GrIS. However, our results also indicate that a significant fraction of the inter-annual dust variability in NE and Central Greenland is determined by the frequency and intensity of wet deposition during the season of high atmospheric dust loading, rather than representing the variability of the Asian dust source and/or long-range transport to Greenland. The variances in the regional proxy records explained by the streamfunction patterns are high enough to permit reconstructions of the corresponding regional deposition regimes and the associated circulation patterns.