Mountain glaciation and paleoclimate reconstruction in the Picos de Europa (Iberian Peninsula,SW Europe)

Geomorphic mapping and stratigraphic analysis of a lake core document the late Quaternary glacial history of the Central and Eastern Massifs of the Picos de Europa, northwestern Spain. The distribution of glacial deposits indicates that at their most advanced positions glaciers occupied 9.1 km², extended as far as 7 km down-valley and had an estimated equilibrium-line altitude (ELA) ranging between 1666 and 1722 m. Radiocarbon dating of sediment deposited in a lake dammed by moraines of this advance show that the maximum glacial extent was prior to 35,280 ± 440 cal yr BP. This advance was followed by two subsequent but less extensive late Pleistocene advances, recorded by multiple moraines flanking both massifs and sedimentary characteristics in the lake deposits. The last recognized glacial episode is the 19th-century maximum extent of small Little Ice Age glaciers in the highest cirques above 2200 m.     

10,000 years of climate control over carbon accumulation in an Iberian bog(southwestern Europe)

The northwest region of the Iberian Peninsula is home to a unique ecosystem of bogs,which are particularly sensitive to projected climate cha nge.In this context,the rate of carbon(C)accumulation in Chao de Veiga Mol.an intact raised bog,was analysed.Changes in the accumulation rate over the past10 millennia were determined in a peat core of 847 cm in depth,with a high mean rate of peat growth(11 yr cm^-1,0.09 cm yr^-1).An age-depth model was generated from 22 14C dates and fallout radionuclides.Chronological,stratigraphical and physico-chemical data confirmed the existence of a single cycle of peat formation throughout the Holocene and the formation of ombrotrophic peat 9500 years ago.The total mean C content was 50.2%,and over 10 millennia 583 kg C m^-2 accumulated at a mean rate of 35.3 g C m^-2 yr^-1,with a long-term(apparent)rate of carbon accumulation in the catotelm of59.9 g C m^-2 yr^-1.These values are much higher than reported for other Iberian peatlands and are amongst the highest documented for peatlands in the northern hemisphere.The dynamics of C accumulation and other measured parameters reveals important variations throughout the Holocene.They could be associated with the main climatic events described in the northern hemisphere and are highly consistent with models established for northern latitudes.The Chao de Veiga Mol raised bog is unique and of great potential value for carrying out high resolution palaeoenvironmental studies,especially in relation to regional and Iocal modulations in southern Europe.     

The Weichselian Late-glacial in southwestern Europe (Iberian Peninsula,Pyrenees, Massif Central,northern Apennines)

This paper presents a summary of Late-glacial environmental changes in southwestern Europe (lberian Peninsula, Pyrenees, Massif Central and the northern Apennines). The emphasis is on palaeoclimatic interpretations inferred from key sites in the region from which the most detailed records are available and which have been radiocarbon dated. The earliest evidence for climatic improvement following the end of the last glacial stage is dated to ca. 15 ka BP and is found at a few sites only. By 13 ka BP, a more widespread and marked climatic improvement is evident, although it is difficult to be precise in the timing and magnitude of the event. There are significant variations in detail between the Late-glacial records, but evidence for a significant cooling correlated with the Younger Dryas event is widespread throughout the region. Just two sites in the region provide evidence for an earlier, less emphatic phase of climatic cooling, which is tentatively equated with the ‘Older Dryas’ of continental northern Europe. Dry conditions appear to have predominated throughout the region in the later part of the Younger Dryas and the early Holocene.     

A procedure to estimate the seismic hazard in an urban area: an application to Acireale (Eastern Sicily)

A “standard procedure” to characterize the seismic hazard of a given area was proposed. It is based on a multidisciplinary approach implying: (1) the knowledge of the seismic history of the area; (2) detailed geological surveys; (3)seismic noise measurements; (4) simulations of earthquake scenarios. The downtown of Acireale, a typical baroque town located on Eastern Sicily, was chosen as the “test area”. A catalog of the local seismogenic faults (able to generate earthquakes in historical times) has been compiled, as well as a seismic catalog for the effects of both local and regional earthquakes. The analysis of both catalogs allowed us to make the following conclusions: (1) the most important seismogenic faults affecting the Acireale municipality do not affect the downtown, while the related local earthquakes attenuate their energy (and intensity) in short (few km) distances; (2) the highest seismic intensity (degree X) experienced in Acireale downtown was caused by the 1693 regional earthquake; (3) over the last 140 years, the downtown has experienced the highest intensity value of VII only once, while six times the intensity was VI. On the whole, this implies a moderate seismic hazard. The estimation of the seismic hazard has been also approached by the experimental method of recording seismic noise. Measurements have been performed at seven different sites, where drills gave detailed information on the shallow subsurface geology to obtain HV (horizontal/vertical) spectral ratios. On the whole, the highest site amplification factor was moderate (about 7). A further investigation based on synthetic seismograms (and spectra) produced by simulating two given earthquake scenarios was also performed. The two scenarios are, respectively, representative of the largest expected earthquake in the area (the 1693 shock) and of a moderate (magnitude ca. 5.5) local earthquake (as the 1818 one). Moderate to strong locally expected accelerations were evidenced.     

Geomorphic classification of rivers: a new methodology applied in an Atlantic Region (Galicia,NW Iberian Peninsula)

The classification of rivers based on geomorphological criteria played, in the past, a secondary role for management decision making, although in the last years they have experienced a shift from the scientific field to that of the technical. Currently, managers require the most simplified form of classifications in order to use them in plans and planning projects, management and restoration. On the one hand, this means that classifications should be directed toward simplifying the diversity of fluvial environments in a number of manageable types, and on the other hand, to apply to each geomorphic type a management model. In this study, we have developed a method of “Geomorphic Classification of Rivers” according to specific stream power variables and median grain size. The new method is dynamic (ability to readjust and gain robustness with the incorporation of new data) and predictive. We obtained six types directly from the method, although we added a seventh type (bedrock rivers) for its special singularity. Each geomorphic type presents concrete values of specific stream power and median grain size, which additionally involves a hydraulic geometry adjusted to terms of balance. The change from a geomorphic type to another would reflect a geomorphologic imbalance in the form of greater power and sediment size. Our results have been compared with six commonly used classifications (e.g., Rosgen Classification System or River Styles Framework). “Geomorphic Classification of Rivers” stands as a flexible tool that allows the development of a “personalized” geomorphic classification for rivers of the same geomorphological province. Through the temporary revision of various sites that act as control points, we can learn, should the case arise, the intensity and geomorphic change of the site. “Geomorphic Classification of Rivers” acts as an alert system for any geomorphologic disturbance. Its simple application and interpretation facilitate the implementation in the administrative environment, or its attachment to other commonly used classifications.     

Biogeochemical characterization of a Mediterranean shallow lake using stable isotopes: Laguna del Cristo (NW Iberian Peninsula)

The present multi-isotopic study (δ¹⁸O–δD_water, δ³⁴S–δ¹⁸O_{dissolved-sulphate}, δ¹³C_{dissolved-inorganic-carbon}, δ¹³C–δ¹⁸O_{shells-modern-gastropods}, δ¹³C_plants, and δ¹³C_{sedimentary-organic-matter}) is aimed at assessing the hydrogeochemical changes and biogeochemical dynamics in a Mediterranean shallow lake fed by a Quaternary–Tertiary aquifer, the “Laguna del Cristo” (NW Iberian Peninsula), a system sensitive to climate fluctuations, between 2010 and 2011. Lake water is of the bicarbonate type, and there are no major pollutants. δ¹⁸O-δD_water values plot on a local evaporation line (δD=5.29δ¹⁸O–12.29) indicating that evaporative enrichment had a significant impact on lake water isotopic features. Periods of high water levels are characterized by lower δ³⁴S–δ¹⁸O_{dissolved-sulphate} and δ¹³C_{dissolved-inorganic-carbon} values and suggest sulphate derived from weathering of sulphides in the catchment area, delivered to the lake by surface run-off or via groundwater, and in situ decay of organic matter. During lower water levels, sulphate reduction and enhanced primary productivity lead to higher δ³⁴S_{dissolved-sulphate} and δ¹³C_{dissolved-inorganic-carbon} values. Evaporation induced enrichment in ¹⁸O_sulphate, ¹³C_{dissolved-inorganic-carbon} and ¹³C–¹⁸O_{shells-Galba-Gyraulus}. δ¹³C_plant confirms the C3 photosynthetic pathway. Enrichment in ¹³C submerged aquatic plants indicates that HCO₃^? is the main carbon source, except for ¹³C-depleted Potamogeton. The TOC, δ¹³C_org values, and TC/TN ratios in sediments all confirm the autochthonous character of organic matter contribution. This study provides a baseline for isotopic research into shallow, flow-through lakes fed by siliciclastic aquifers, and stresses the importance of evaporation and refilling (direct precipitation and groundwater discharge) in controlling the solute chemistry and stable isotopic composition in temperate regions with contrasting seasonal climates. The results also provide a snapshot of modern lake isotope variability that can be applied to paleoenvironmental reconstructions.     

Recent trends on glacier area retreat over the group of Nevados Caullaraju-Pastoruri (Cordillera Blanca,Peru) using Landsat imagery

The Cordillera Blanca, located in the central zone of the Andes Mountains in Peru, has shown a retreat in its glaciers. This paper presents a trend analysis of the glacier area over the groups of Nevados Caullaraju-Pastoruri from 1975 to 2010 using Landsat-5 Thematic Mapper (TM) imagery. In the case of the Nevados Pastoruri/Tuco, the study period was extended back to 1957 by using an aerial photograph taken that year. The extent of clean glacier ice was estimated using Normalized Difference Snow Index (NDSI) thresholds. Moreover, the estimation of debris-covered glacier ice was retrieved by means of a decision tree classification method using NDSI, Normalized Difference Vegetation Index (NDVI) and Land Surface Temperature (LST). Area estimations derived from Landsat imagery were compared to the glacier ground-truth data in 1975 and 2010. Results show a statistically significant (p < 0.05) decreasing trend over the whole study area. Total glacier area decreased at a rate of 4.5 km² per decade from 1975 to 2010, with a total loss of 22.5 km² (58%). Lower decreasing rates were found for the period 1987–2010: 3.5 km² per decade with a total loss of 7.7 km² (32.5%). In the case of the Nevados Pastoruri/Tuco, decreasing rates of clean ice extent were constant for the periods 1957–2010, 1975–2010 and 1987–2010, with values close to 1.4 km² per decade and a total loss between 1957 and 2010 estimated at about 5 km² (54%). This work shows an evident area decrease in the Caullaraju-Pastoruri tropical glaciers, which needs to be included in a future hydrological scenario of local adaptability and water management.     

Development of a parameterization for simulating the urban temperature hazard using satellite observations in climate model

Urban surface temperature is hazardously higher than surrounding regions (so-called urban heat island effect UHI). Accurately simulating urbanization-induced temperature hazard is critical for realistically representing urban regions in the land surface-atmosphere climate system. However, inclusion of urban landscapes in regional or global climate models has been overlooked due to the coarse spatial resolution of these models as well as the lack of observations for urban physical properties. Recently, National Aeronautics and Space Administration (NASA) Earth Observing System (EOS) Moderate Resolution Imaging Spectroradiometer (MODIS) observations illustrate important urban physical properties, including skin temperature, surface albedo, surface emissivity, and leaf area index, It is possible to identify the unique urban features globally and thus simulate global urban processes. An urban scheme is designed to represent the urban-modified physical parameters (albedo, emissivity, land cover, roughness length, thermal and hydraulic properties) and to include new, unique physical processes that exist in urban regions. The urban scheme is coupled with National Center for Atmospheric Research (NCAR) Community Land Model Version 2 (CLM2) and single column coupled NCAR Community Atmosphere Model CAM2/CLM2 to assess the mechanisms responsible for UHI. There are two-steps in our model development. First, satellite observations of albedo, emissivity, LAI, and in situ observed thermal properties are updated in CLM2 to represent the first-order urban effects. Second, new terms representing the urban anthropogenic heat flux, storage heat flux, and roughness length are calculated in the model. Model simulations suggest that human activity-induced surface temperature hazard results in overlying atmosphere instability and convective rainfall, which may enhance the possibility of urban flood hazard.

Study of vertical wind profiles in an urban area with complex terrain (Tehran)

In this paper monthly trends of vertical wind profiles within and above an urban area with complex topography (Tehran) have been investigated using data from a Sodar, a 100-m and two 28-m towers and some surface stations. It includes the occurrence, evolution, dissipation time, peak time and maximum height for katabatic–anabatic winds due to topography, return flow and urban circulation. Vertical distributions of wind and the heat and momentum fluxes up to 600 m were also considered.     

The transition from an active to a passive margin (SW end of the South Shetland Trench, Antarctic Peninsula)

The lateral ending of the South Shetland Trench is analysed on the basis of swath bathymetry and multichannel seismic profiles in order to establish the tectonic and stratigraphic features of the transition from an northeastward active to a southwestward passive margin style. This trench is associated with a lithospheric-scale thrust accommodating the internal deformation in the Antarctic Plate and its lateral end represents the tip-line of this thrust. The evolutionary model deduced from the structures and the stratigraphic record includes a first stage with a compressional deformation, predating the end of the subduction in the southwestern part of the study area that produced reverse faults in the oceanic crust during the Tortonian. The second stage occurred during the Messinian and includes distributed compressional deformation around the tip-line of the basal detachment, originating a high at the base of the slope and the collapse of the now inactive accretionary prism of the passive margin. The initial subduction of the high at the base of the slope induced the deformation of the accretionary prism and the formation of another high in the shelf—the Shelf Transition High. The third stage, from the Early Pliocene to the present-day, includes the active compressional deformation of the shelf and the base-of-slope around the tip-line of the basal detachment, while extensional deformations are active in the outer swell of the trench.     

The role of rainfall in the landslide hazard: the case of the Avigliano urban area (Southern Apennines, Italy)

Mass movements varying in type and size, some of which are periodically reactivated, affect the urban area of Avigliano. The disturbed and remoulded masses consist of sandy–silty or silty–clayey plastic material interbedded with stone fragments and conglomerate blocks. Five landslides that were markedly liable to rainfall-associated instability phenomena were selected.

The relationships between landslides and rainfall were investigated using a hydrological and statistical model based on long-term series of daily rainfall data. The model was used to determine the return period of cumulative daily rainfall over 1–180 days. The resulting hydrological and statistical findings are discussed with the aim of identifying the rainfall duration most critical to landslides.

The concept of a precipitation threshold was generalized by defining some probability classes of cumulative rainfall. These classes indicate the thresholds beyond which reactivation is likely to occur. The probability classes are defined according to the return period of the cumulative rainfall concomitant with landslide reactivation.     

New evidence for the presence of reindeer (Rangifer tarandus) on the Iberian Peninsula in the Pleistocene: an archaeopalaeontological and chronological reassessment

Reindeer (Rangifer tarandus) populations reached the Iberian Peninsula during the end of the Middle Pleistocene and there are numerous examples of this species from Late Pleistocene Mousterian and Upper Palaeolithic sites. In this paper, new evidence of reindeer in the east Cantabrian region is presented that further supports the timing of this species’ first appearance, and previous assessments are updated. To date, the presence of this species has been identified at 55 sites in the Iberian Peninsula, nearly as many as those of mammoth (Mammuthus primigenius) and woolly rhinoceros (Coelodonta antiquitatis) combined. Most of the sites with presence of reindeer (50) are located in the Cantabrian region with a clear increase in the density of sites and remains towards the Pyrenees. The remaining five sites with evidence of reindeer are located on the other side of the Pyrenees in the NW corner of Catalonia. In contrast, archaeological evidence of reindeer in the form of art (both parietal and portable) is more scarce and scattered. Evidence for the representation of these animals has been found outside the northern fringe of the Iberian Peninsula, which could reflect either long‐distance cultural communication or the movement of human groups.     

Tree failure prediction model (TFPM): machine learning techniques comparison in failure hazard assessment of Platanus orientalis in urban forestry

Natural Hazards - Trees are generally harmed by multitude factors consisting of ecological condition and anthropogenic pressures in the cities. This study compares the multilayer perceptron (MLP)...     

3-D electromagnetic imaging of a Palaeozoic plate-tectonic boundary segment in SW Iberian Variscides (S Alentejo, Portugal)

In SW Iberian Variscides, the boundary between the South Portuguese Zone (SPZ) and the Ossa Morena Zone (OMZ) corresponds to a major tectonic suture that includes the Beja Acebuches Ophiolite Complex (BAOC) and the Pulo do Lobo Antiform Terrane (PLAT). Three sub-parallel and approximately equidistant MT profiles were performed, covering a critical area of this Palaeozoic plate-tectonic boundary in Portugal; the profiles, running roughly along an NE–SW direction, are sub-perpendicular to the main Variscan tectonic features. Results of the three-dimensional (3-D) modelling of MT data allow to generate, for the first time, a 3-D electromagnetic imaging of the OMZ–SPZ boundary, which reveals different conductive and resistive domains that display morphological variations in depth and are intersected by two major sub-vertical corridors; these corridors coincide roughly with the NE–SW, Messejana strike–slip fault zone and with the WNW–ESE, Ferreira–Ficalho thrust fault zone. The distribution of the shallow resistive domains is consistent with the lithological and structural features observed and mapped, integrating the expected electrical features produced by igneous intrusions and metamorphic sequences of variable nature and age. The development in depth of these resistive domains suggests that: (1) a significant vertical displacement along an early tectonic structure, subsequently re-taken by the Messejana fault-zone in Late-Variscan times, has to be considered to explain differences in deepness of the base of the Precambrian–Cambrian metamorphic pile; (2) hidden, syn- to late-collision igneous bodies intrude the meta-sedimentary sequences of PLAT; (3) the roots of BAOC are inferred from 12 km depth onwards, forming a moderate resistive band located between two middle-crust conductive layers extended to the north (in OMZ) and to the south (in SPZ). These conductive layers overlap the Iberian Reflective Body (evidenced by the available seismic reflection data) and are interpreted as part of an important middle-crust décollement developed immediately above or coinciding with the top of a graphite-bearing granulitic basement.     

Analysis of seismic hazard in landslide-prone regions: criteria and example for an area of Daunia (southern Italy)

In relation to the assessment of earthquake-induced landslide hazard, this paper discusses general principles and describes implementation criteria for seismic hazard estimates in landslide-prone regions. These criteria were worked out during the preparation of a hazard map belonging to the official Italian geological cartography and they are proposed as guidelines for future compilation of similar maps. In the presented case study, we used a procedure for the assessment of seismic hazard impact on slope stability adopting Arias intensity Ia as seismic shaking parameter and critical acceleration a _c as parameter representing slope strength to failures induced by seismic shaking. According to this procedure, after a preliminary comparison of estimated historical maximum values of Ia with values proposed in literature as landslide-triggering thresholds, a probabilistic approach, based on the Newmark’s model, is adopted: it allows to estimate the minimum critical acceleration a _c required for a slope to keep under a prefixed value, the probability of failures induced by seismic shakings expected in a given time interval. In this way, one can prepare seismic hazard maps where seismic shaking is expressed in an indirect way through a parameter (the critical acceleration) representing the “strength” that seismic shakings mobilise in slope materials (strength demand) with a prefixed exceedance probability. This approach was applied to an area of Daunia (Apulia—southern Italy) affected by frequent landslide phenomena. The obtained results indicate that shakings with a significant slope destabilisation potential can be expected particularly in the north-western part of the area, which is exposed to the seismic activity of Apennine tectonic structures.     

Hydrothermal alteration of felsic volcanic rocks associated with massive sulphide deposition in the northern Iberian Pyrite Belt (SW Spain)

Massive sulphide deposits of the northern Iberian Pyrite Belt (IPB) are mainly hosted by felsic volcanic rocks of rhyolitic to dacitic composition. Beneath most of the massive ores of this area (e.g., Concepción, San Miguel, Aguas Teñidas Este or San Telmo deposits) there is usually a wide hydrothermal alteration halo associated with stockwork-type mineralization. Within these alteration envelopes there are two principal rock types: (1) chlorite-rich rocks, linked to the inner and more intensely altered zones and dominantly comprising chlorite+pyrite+quartz+sericite (+carbonate+rutile+zircon+chalcopyrite), and (2) sericite-rich rocks, more common in the peripheral zones and showing a dominant paragenesis of sericite+quartz+pyrite+chlorite (+carbonate+rutile+zircon+sphalerite). Mass-balance calculations comparing altered and least-altered felsic volcanic rocks suggest that sericitization was accompanied by moderate enrichment in Mg, Fe and H₂O, with depletion in Si, Na and K, and a slight net mass loss of about 3%. Chloritization shows an overall pattern which is similar to that of the sericitic alteration, but with large gains in Fe, Mg and H₂O (and minor enrichment in Si, S and Mn), and a significant loss of Na and K and a minor loss of Ca and Rb. However, chloritization has involved a much larger net mass change (mass gain of about 28%). Only a few elements such as Nb, Y, Zr, Ti, P and LREE appear to have remained inert during hydrothermal alteration, whilst Ti and Al have undergone very minor mobilization. The results point to the severity of the physico-chemical conditions that prevailed during the waxing stage of the ore-forming hydrothermal systems. Further, mineralogical and geochemical studies of the altered footwall rocks in the studied deposits indicate that hydrothermal ore-bearing fluids reacted with host rocks in a multi-stage process which produced a succession of mineralogical and chemical changes as the temperature increased.     

Geological and Palaeontological context of three new Barremian (Lower Cretaceous) vertebrate sites in the Iberian Peninsula (Cuenca Province,Central Spain)

Three new Lower Cretaceous vertebrate sites (Vadillos-1, Vadillos-2, El Tobar) have been recently discovered and studied in the Cuenca Province (Central Spain). They are located in deposits of “Wealden” facies belonging to the El Collado Sandstone and Clay Formation. In these outcrops, micro and macroremains corresponding to plants, invertebrates and vertebrates have been collected and subsequently assigned to macrophytes, charophytes (e.g., Atopochara trivolvis triquetra, Globator maillardii trochiliscoides, Clavator harrisii harrisii), ostracods (e.g., Cypridea gr. modesta, Cypridea cf. C. isasae, Cypridea sp. aff. C. moneta, Cypridea sp. 1, Cypridea sp. 2), molluscs (Unionoida, Viviparus sp.), fishes, amphibians, turtles (cf. Eucryptodira), crocodyliforms (Neosuchia) and dinosaurs (ankylosaurs, ornithopods, theropods). Among the vertebrate remains, scales, teeth, plates, osteoderms, phalanges, ribs, vertebrae and other incomplete bones, as well as eggshell fragments have been identified. This rich and diverse assemblage was deposited in an upper Barremian alluvial-palustrine muddy floodplain crossed by braided sandy channels.     

Electrical resistivity imaging over natural (in situ) geological samples using physical model studies

To understand the characteristic responses of natural geological samples, viz., black granite, green marble, and graphite sheet, and to have “an a priori” knowledge of their physical property through electrical resistivity imaging, the physical model laboratory setup has been established to conduct scale model studies over targets of finite dimensions and resistivities. The present experiment involves IRIS make SYSCAL Pro-96 measuring system using 48 electrodes with 2-cm interelectrode separation in the laboratory model tank. In the present communication, we have presented the 2D cross section images using Wenner, Wenner–Schlumberger, and dipole–dipole array configurations over the resistive (granite, marble) and conductive (graphite) sheets. In the case of resistive target (black granite sheet, green marble), the combined usage of dipole–dipole and Wenner–Schlumberger arrays provided more accurate measures on target parameters, i.e., the combined usage of both the arrays is preferable in searching high-resistive targets beneath the low-resistive ones over burden. The shape of the resistive target (green marble sheet) is inappropriate when the thickness of the target is greater than one half of the minimum array separation. As the thickness of the target increases, the signatures of the target become feeble, and hence, the shape of the resistive target is not properly reflected in the corresponding tomogram. The response over graphite sheet indicates that the true parameters of the target are not reflected in the cross section, and the existence of the low-resistive (high-conductive) target in the host medium (water) deviates the resistivity of the medium. The target parameters from the cross section using dipole–dipole array are somewhat correlated with true parameters in the case of thin targets at shallow depths. In the case of the sequence of layers of gravel–marble gravel–sand gravel simulated in a small model tank in the physical model laboratory, the thickness of the high-resistive marble layer beneath the low-resistive gravel layer is enhanced conspicuously because of the significant resistivity contrast between gravel and marble.