Eustatic and tectonic controls on Quaternary Ras Leona marine terraces (Strait of Gibraltar, northern Morocco)

Well-preserved Quaternary staircased marine terraces appear on Ras Leona limestone relief. This is a peculiar sector of the Betic-Rif Cordillera, lying in the four-way junction between the Atlantic and the Mediterranean, and Europe and Africa. The age and altitude correlation of the Ras Leona terraces with travertine-covered lateral equivalent terraces fashioned in the neighbouring Beni Younech area, and comparison with those along the Moroccan Atlantic coasts, would suggest that the Ras Leona terraces were mainly formed by eustatic factors. The importance of the eustasy is supported by further comparisons with Spanish and Moroccan Mediterranean terraces and with different marine terraces developed on passive-margin coasts around the world. A tectonic event occurred mainly during the period between the formation of the Maarifian and the Ouljian terraces (i.e., between 370 and 150 ka). The moderate Quaternary tectonic uplift deduced from the marine terraces and its comparison with uplifted marine terraces developed in active subduction setting disagrees with the model of an active eastwards subduction below the Gibraltar tectonic arc.     

Rhizosphere and flooding regime as key factors for the mobilisation of arsenic and potentially harmful metals in basic,mining-polluted salt marsh soils

The occurrence of mining areas in the vicinities of salt marshes may affect their ecological functions and facilitate the transfer of pollutants into the food chain. The mobilisation of metals in salt marsh soils is controlled by abiotic (pH, redox potential) and biotic (influence of rhizosphere) factors. The effect of the rhizosphere of two plant species (Sarcocornia fruticosa and Phragmites australis) and different flooding regimes on potentially harmful metals and As mobilisation from salt marsh soil polluted by mining activities were investigated (total concentrations: 536 mg kg⁻¹ As, 37 mg kg⁻¹ Cd, 6746 mg kg⁻¹ Pb, 15,320 mg kg⁻¹ Zn). The results show that the changes in redox conditions (from 300 mV to −100 mV) and pH after flooding and rewetting periods may mobilise the contaminant elements into soil solution (e.g., 100 μg L⁻¹ Cd, 30 μg L⁻¹ Pb, 7 mg L⁻¹ Zn), where they are available for plants or may be leached from the soil. Drying periods generated peaks of concentrations in the soil solution (up to 120 μg L⁻¹ Cd and 50 μg L⁻¹ Pb). The risk assessment of As and metal-polluted salt marshes should take into account flood dynamics in order to prevent metal(loid) mobilisation.     

Projected climate change impacts on spatial distribution of bioclimatic zones and ecoregions within the Kailash Sacred Landscape of China,India, Nepal

Rapidly accelerating climate change in the Himalaya is projected to have major implications for montane species, ecosystems, and mountain farming and pastoral systems. A geospatial modeling approach based on a global environmental stratification is used to explore potential impacts of projected climate change on the spatial distribution of bioclimatic strata and ecoregions within the transboundary Kailash Sacred Landscape (KSL) of China, India and Nepal. Twenty-eight strata, comprising seven bioclimatic zones, were aggregated to develop an ecoregional classification of 12 ecoregions (generally defined by their potential dominant vegetation type), based upon vegetation and landcover characteristics. Projected climate change impacts were modeled by reconstructing the stratification based upon an ensemble of 19 Earth System Models (CIMP5) across four Representative Concentration Pathways (RCP) emission scenarios (i.e. 63 impact simulations), and identifying the change in spatial distribution of bioclimatic zones and ecoregions. Large and substantial shifts in bioclimatic conditions can be expected throughout the KSL area by the year 2050, within all bioclimatic zones and ecoregions. Over 76 % of the total area may shift to a different stratum, 55 % to a different bioclimatic zone, and 36.6 % to a different ecoregion. Potential impacts include upward shift in mean elevation of bioclimatic zones (357 m) and ecoregions (371 m), decreases in area of the highest elevation zones and ecoregions, large expansion of the lower tropical and sub-tropical zones and ecoregions, and the disappearance of several strata representing unique bioclimatic conditions within the KSL, with potentially high levels of biotic perturbance by 2050, and a high likelihood of major consequences for biodiversity, ecosystems, ecosystem services, conservation efforts and sustainable development policies in the region.     

Impact of anthropogenic CO 2 on the next glacial cycle

The model of Paillard and Parrenin (Earth Planet Sci Lett 227(3–4):263–271, 2004) has been recently optimized for the last eight glacial cycles, leading to two different relaxation models with model-data correlations between 0.8 and 0.9 (García-Olivares and Herrero (Clim Dyn 1–25, 2012b)). These two models are here used to predict the effect of an anthropogenic CO ₂ pulse on the evolution of atmospheric CO ₂, global ice volume and Antarctic ice cover during the next 300 kyr. The initial atmospheric CO ₂ condition is obtained after a critical data analysis that sets 1300 Gt as the most realistic carbon Ultimate Recoverable Resources (URR), with the help of a global compartmental model to determine the carbon transfer function to the atmosphere. The next 20 kyr will have an abnormally high greenhouse effect which, according to the CO ₂ values, will lengthen the present interglacial by some 25 to 33 kyr. This is because the perturbation of the current interglacial will lead to a delay in the future advance of the ice sheet on the Antarctic shelf, causing that the relative maximum of boreal insolation found 65 kyr after present (AP) will not affect the developing glaciation. Instead, it will be the following insolation peak, about 110 kyr AP, which will find an appropriate climatic state to trigger the next deglaciation.     

Linear system techniques applied to the Fuenmayor Karst Spring, Huesca (Spain)

Black box identification techniques are very useful tools to study the relationship between rainfall and the discharge of karst aquifers. Autocorrelation, cross-correlation, cross-spectral, impulse response identification (Wiener–Hopf and error minimization methods) and parametric identification of transfer functions techniques have been used by different authors. A review of these techniques has been realized to assess the possibilities and limitations. The reviewed techniques are applied to a time rainfall–discharge series of Fuenmayor karst spring, (southern central Pyrenees, Spain) to establish the relative quality of the obtained models. The quality is evaluated with the Nash–Sutcliffe model efficiency coefficient. In addition, effective rainfall has been used to improve the performance of every model. The whitening technique has been used in the cross-correlation technique with good results. The studied spring shows a Fontestorbes type behavior. The linearity of the system has been analyzed by the coherence function. Fuenmayor has an acceptable linear response, except for high frequency events. The parametric identification provides the best efficiency value E = 0.8164 and provides a proper transfer function, with five parameters (one zero, two poles, the gain and the delay). The prediction power of this transfer function is validated with another time series, with an efficiency value of E = 0.7904.     

Spatial modelling and uncertainty assessment of pyroclastic cover thickness in the Sorrento Peninsula

Assessing spatial variability of soil thickness is a critical issue for understanding and predicting slope processes. The present work was aimed at estimating the spatial scales at which the variation of pyroclastic cover thickness occurs in a sample area in the Sorrento Peninsula (Italy). Stochastic simulation was used to understand the spatial variability of pyroclastic cover thickness on Mount Pendolo and to assess its spatial uncertainty. In the study area, covering about 0.7 km², thickness measurements were collected using electrical resistivity tomography profiles, continuous core drillings and steel rod penetrometric tests. Variographic analysis revealed the occurrence of an anisotropic behaviour along the N50 and N140 directions. In the latter anisotropic direction, a nested variogram was fitted including (1) a long-range component which could be related to large-scale factors, like the curvature of the slope and contributing area and (2) a shorter scale variation which is probably associated with the occurrence of denudation processes or to the articulate cover/bedrock interface. To assess the spatial variability and uncertainty of pyroclastic cover thickness, a stochastic simulation algorithm was used and 500 equally probable images of cover thickness were yielded. The results showed that a better thickness distribution map can be drawn by simulating the data collected on the slope and at the footslope separately. The approach also allowed delineating the areas characterized by greater uncertainty, suggesting supplementary measurements to further improve the cover thickness distribution model, thus reducing the uncertainty.     

Multi-Modal Monitoring of Slip Along Frictional Discontinuities

Seismic wave transmission and digital image correlation (DIC) are employed to study slip processes along frictional discontinuities. A series of biaxial compression experiments are performed on gypsum specimens with non-homogeneous contact surfaces. The specimens are composed of two blocks with perfectly mated contact surfaces with a smooth surface with low frictional strength on the upper half and a rough surface with high frictional strength on the lower half. Compressional, P, and shear, S, wave pulses were transmitted through the discontinuity while digital images of the specimen surface were acquired during the test. A distinct peak in the amplitude of transmitted wave occurs prior to the peak shear strength and is considered a “precursor” to the failure. Precursors indicate that slip initiates from the smooth surface and extends to the rough surface as the shear load is increased. From the DIC data, slip is identified as a jump in the displacement field along the fracture that initiates from the smooth surface and propagates to the rough surface. Precursors are associated with an increase in the rate of slip across the discontinuity and are a measure of the reduction in the fracture shear stiffness.