Long-term deforestation in NW Spain: linking the Holocene fire history to vegetation change and human activities

The Holocene fire regime is thought to have had a key role in deforestation and shrubland expansion in Galicia (NW Spain) but the contribution of past societies to vegetation burning remains poorly understood. This may be, in part, due to the fact that detailed fire records from areas in close proximity to archaeological sites are scarce. To fill this gap, we performed charcoal analysis in five colluvial soils from an archaeological area (Campo Lameiro) and compared the results to earlier studies from this area and palaeo-ecological literature from NW Spain. This analysis allowed for the reconstruction of the vegetation and fire dynamics in the area during the last ca 11 000 yrs. In the Early Holocene, Fabaceae and Betula sp. were dominant in the charcoal record. Quercus sp. started to replace these species around 10 000 cal BP, forming a deciduous forest that prevailed during the Holocene Thermal Maximum until ～5500 cal BP. Following that, several cycles of potentially fire-induced forest regression with subsequent incomplete recovery eventually led to the formation of an open landscape dominated by shrubs (Erica sp. and Fabaceae). Major episodes of forest regression were (1) ～5500–5000 cal BP, which marks the mid-Holocene cooling after the Holocene Thermal Maximum, but also the period during which agropastoral activities in NW Spain became widespread, and (2) ～2000–1500 cal BP, which corresponds roughly to the end of the Roman Warm Period and the transition from the Roman to the Germanic period. The low degree of chronological precision, which is inherent in fire history reconstructions from colluvial soils, made it impossible to distinguish climatic from human-induced fires. Nonetheless, the abundance of synanthropic pollen indicators (e.g. Plantago lanceolata and Urtica dioica) since at least ～6000 cal BP strongly suggests that humans used fire to generate and maintain pasture.     

Low-order stabilized finite element for the full Biot formulation in soil mechanics at finite strain

This article presents a novel finite element formulation for the Biot equation using low-order elements. Additionally, an extra degree of freedom is introduced to treat the volumetric locking steaming from the effective response of the medium; its balance equation is also stabilized. The accuracy of the proposed formulation is demonstrated by means of numerical analyses.     

Spatiotemporal estimation of snow depth using point data from snow stakes,digital terrain models,and satellite data

Snow availability in Alpine catchments plays an important role in water resources management. In this paper, we propose a method for an optimal estimation of snow depth (areal extension and thickness) in Alpine systems from point data and satellite observations by using significant explanatory variables deduced from a digital terrain model. It is intended to be a parsimonious approach that may complement physical‐based methodologies. Different techniques (multiple regression, multicriteria analysis, and kriging) are integrated to address the following issues: We identify the explanatory variables that could be helpful on the basis of a critical review of the scientific literature. We study the relationship between ground observations and explanatory variables using a systematic procedure for a complete multiple regression analysis. Multiple regression models are calibrated combining all suggested model structures and explanatory variables. We also propose an evaluation of the models (using indices to analyze the goodness of fit) and select the best approaches (models and variables) on the basis of multicriteria analysis. Estimation of the snow depth is performed with the selected regression models. The residual estimation is improved by applying kriging in cases with spatial correlation. The final estimate is obtained by combining regression and kriging results, and constraining the snow domain in accordance with satellite data. The method is illustrated using the case study of the Sierra Nevada mountain range (Southern Spain). A cross‐validation experiment has confirmed the efficiency of the proposed procedure. Finally, although it is not the scope of this work, the snow depth is used to asses a first estimation of snow water equivalent resources.     

Pre‐Alpine discordant granitic dikes in the metamorphic core of the Betic Cordillera: tectonic implications

Pegmatite dikes bearing andalusite crosscut foliation S₂ in Alpujarride gneisses and schists. Post‐S₂ andalusite is transposed by a foliation S₃, defined by fibrolite, which affects the dikes. The dikes represent highly differentiated granitic magmas with low REE and Zr contents and a positive Eu anomaly. U‐Pb SHRIMP dating of magmatic zircons provided Pan‐African ages (cores) and late Variscan ages (rims). However, U‐rich rims also provided metamorphic Alpine ages, supporting a polyorogenic tectonometamorphic history for pre‐Mesozoic Alpujarride rocks.     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part 1: Saturated state

The behavior of granular materials is known to depend on its loose or dense nature, which in turns depends both on density and confining pressure. Many models developed in the past require the use of different sets of constitutive parameters for the same material under different confining pressures. The purpose of this paper is to extend a basic generalized plasticity model for sands proposed by Pastor, Zienkiewicz and Chan by modifying the main ingredients of the model flow—rule, loading–unloading discriminating direction and plastic modulus—to include a dependency on the state parameter. The proposed model is tested against the available experimental data on three different sands, using for each of them a single set of material parameters, finding a reasonably good agreement between experiments and predictions. Copyright © 2010 John Wiley & Sons, Ltd.     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part II: Unsaturated soil modeling

The aim of this paper is to extend the generalized plasticity state parameter‐based model presented in part 1 to reproduce the hydro‐mechanical behavior of unsaturated soils. The proposed model is based on two pairs of stress–strain variables and a suitable hardening law taking into account the bonding—debonding effect of suction and degree of saturation. A generalized state parameter for unsaturated state is proposed to reproduce soil behavior using a single set of material parameters. Generalized plasticity gives a suitable framework to reproduce not only monotonic stress path but also cyclic behavior. The hydraulic hysteresis during a drying—wetting cycle and the void ratio effect on the hydraulic behavior is introduced. Comparison between model simulations and a series of experimental data available, both cohesive and granular, are given to illustrate the accuracy of the enhanced generalized plasticity equation. Copyright © 2010 John Wiley & Sons, Ltd.     

Throughfall patterns in sugarcane and riparian forest: understanding the effect of sugarcane age and land use conversion

Sugarcane is an annual crop with a dynamic canopy that changes over time mainly because of genetic adaptation. There is uncertainty about the temporal trends of throughfall (TF) in this important commercial crop. In the present paper, we used troughs to measure TF in a third and fourth ratoon and subsequently in a fourth and fifth ratoon. Additional measurements were carried out in an adjacent riparian forest. There were no significant differences between cycles of sugarcane, growth phases and riparian forest. The TF results for ratoon crop and riparian forest in 2011/2012 were 76% and 79.5% of gross rainfall, respectively, while in 2012/2013, they were 79% and 78%, respectively. However, TF was remarkably lower in the riparian forest relative to ratoon from the second half of the culm formation and elongation phase (280 days after harvest) until harvest. Copyright © 2016 John Wiley & Sons, Ltd.     

P–T–t–d evolution of orogenic middle crust of the Roc de Frausa Massif (Eastern Pyrenees): a result of horizontal crustal flow and Carboniferous doming?

Structural, petrological and textural studies are combined with phase equilibria modelling of metapelites from different structural levels of the Roc de Frausa Massif in the Eastern Pyrenees. The pre‐Variscan lithological succession is divided into the Upper, Intermediate and Lower series by two orthogneiss sheets and intruded by Variscan igneous rocks. Structural analysis reveals two phases of Variscan deformation. D1 is marked by tight to isoclinal small‐scale folds and an associated flat‐lying foliation (S1) that affects the whole crustal section. D2 structures are characterized by tight upright folds facing to the NW with steep NE–SW axial planes. D2 heterogeneously reworks the D1 fabrics, leading to an almost complete transposition into a sub‐vertical foliation (S2) in the high‐grade metamorphic domain. All structures are affected by late open to tight, steeply inclined south‐verging NW–SE folds (F3) compatible with steep greenschist facies dextral shear zones of probable Alpine age. In the micaschists of the Upper series, andalusite and sillimanite grew during the formation of the S1 foliation indicating heating from 580 to 640 °C associated with an increase in pressure. Subsequent static growth of cordierite points to post‐D1 decompression. In the Intermediate series, a sillimanite–biotite–muscovite‐bearing assemblage that is parallel to the S1 fabric is statically overgrown by cordierite and K‐feldspar. This sequence points to ~1 kbar of post‐D1 decompression at 630–650 °C. The Intermediate series is intruded by a gabbro–diorite stock that has an aureole marked by widespread migmatization. In the aureole, the migmatitic S1 foliation is defined by the assemblage biotite–sillimanite–K‐feldspar–garnet. The microstructural relationships and garnet zoning are compatible with the D1 pressure peak at ~7.5 kbar and ~750 °C. Late‐ to post‐S2 cordierite growth implies that F2 folds and the associated S2 axial planar leucosomes developed during nearly isothermal decompression to <5 kbar. The Lower series migmatites form a composite S1–S2 fabric; the garnet‐bearing assemblage suggests peak P–T conditions of >5 kbar at suprasolidus conditions. Almost complete consumption of garnet and late cordierite growth points to post‐D2 equilibration at <4 kbar and <750 °C. The early metamorphic history associated with the S1 fabric is interpreted as a result of horizontal middle crustal flow associated with progressive heating and possible burial. The upright F2 folding and S2 foliation are associated with a pressure decrease coeval with the intrusion of mafic magma at mid‐crustal levels. The D2 tectono‐metamorphic evolution may be explained by a crustal‐scale doming associated with emplacement of mafic magmas into the core of the dome.     

8000 years of vegetation history in the northern Iberian Peninsula inferred from the palaeoenvironmental study of the Zalama ombrotrophic bog (Basque‐Cantabrian Mountains,Spain)

This paper focuses on pollen, spores, non‐pollen palynomorphs (NPPs) and certain geochemical elements from the ombrotrophic blanket bog of Zalama (Basque‐Cantabrian Mountains, northern Iberian Peninsula), with the support of a robust chronology based on 17 AMS ¹⁴C dates. The main results related to the last 8000 years show that, during the early middle Holocene, pines and deciduous forests were the most extensive tree formations. At the beginning of the succession, pines reach 44%, showing regional presence, whereas after 7600 cal. a BP, deciduous forests were particularly abundant. From c. 6500 cal. a BP the pollen diagram constructed from our samples shows the first anthropogenic evidence, linked with the new economic practices related to the Neolithic of the Basque‐Cantabrian Mountains. From 3300 cal. a BP the expansion of Fagus sylvatica is particularly clear, and has since then become one of the dominant forest species in this region. We also discuss the Holocene evolution of other noteworthy plant communities in southwestern Europe, such as Taxus baccata, Juglans and shrublands.     

Upper Cretaceous marine‐continental transition (Leonese Area,NW Spain) defined from integrated outcrop and seismic stratigraphy

The Upper Cretaceous succession of the Leonese Area (NW Spain) comprises mixed clastic and carbonate sediments. This succession is divided into two lithostratigraphic units, the Voznuevo Member and the Boñar Formation, which represent fluvial, shoreface, intertidal, subtidal and open‐shelf sedimentary environments. Regional seismic interpretation and sequence stratigraphic analysis have allowed the study of lateral and vertical changes in the sedimentary record and the definition of third‐order levels of stratigraphic cyclicity. On the basis of these data, the succession can be divided into two second‐order depositional sequences (DS‐1 and DS‐2), incorporating three system tracts in a lowstand to transgressive to highstand system tract succession (LST–TST–HST). These sequences are composed of fluvial systems at the base with palaeocurrents that flowed westward and south‐westward. The upper part of DS‐1 (Late Albian–Middle Turonian) shows evidence of intertidal to subtidal and offshore deposits. DS‐2 (Late Turonian–Campanian) comprises intertidal to subtidal, tidal flat, shallow marine and lacustrine deposits and interbedded fluvial deposits. Two regressive–transgressive cycles occurred in the area related to eustatic controls. The evolution of the basin can be explained by base‐level changes and associated shifts in depositional trends of successive retrogradational episodes. By using isobath and isopach maps, the main palaeogeographic features of DS‐1 and DS‐2 were constrained, namely coastline positions, the existence and orientation of corridors through which fluvial networks were channelled and the location of the main depocentres of the basin. Sedimentation on the Upper Cretaceous marine platform was mainly controlled by (i) oscillations of sea level and (ii) the orientation of Mesozoic faults, which induced sedimentation along depocentres. Copyright © 2013 John Wiley & Sons, Ltd.