An attenuation study using earthquakes from the 1997 Umbria-Marche sequence

We studied spatial and temporal characteristics of seismic attenuation inCentral Italy using S- and coda- waves recorded by the MarchesanSeismograph Network from earthquakes located in the epicentral area ofthe 1997 Umbria-Marche sequence. The amplitude decay of the S waveswith distance was defined calculating empirical attenuation functions at 15frequencies between 1 and 25 Hz. We analyzed separately foreshocks andaftershocks and we found the same attenuation functions, suggesting thatthe possible temporal variations could be confined in a small area. Thefrequency dependence of Q _S was approximated by the equation Q _S=18 · f ^2.0between 1 and 10 Hz. At higher frequencies (10–25 Hz), the frequencydependence of Q _s weakens, having an average value of Q _S=990. We also estimated Q from coda waves (Q _C) using the single-scattering models of Aki andChouet (1975) and Sato (1977). We found that Q _C=77 · f ^0.6, (between 2 and 20Hz) at the western side of the mountain chain, using either foreshocks oraftershocks. This relation is consistent with previous estimates of Q _Creported for the Central Apennines. For a volume sampling the Colfioritobasin, the Apennines and the Marche region we found that Q _C=55 · f ^0.8,indicating highattenuation below the mountain belt. To detect small temporal changes ofQ, we calculated spectral ratios of 5 temporal doublets located in theepicentral area and recorded at the closest station. We found temporalchanges of Q that vary from 27% to 56%, depending on the locationof the doublets. This variability suggests that the temporal change ofattenuation may depend on the spatial variation of Q and perhaps on thespatial distribution of tectonic stress in the epicentral area.     

The ecohydrological imprint of deforestation in the semiarid Chaco: insights from the last forest remnants of a highly cultivated landscape

The semiarid Chaco plains present one of the highest rates of forest clearing and agricultural expansion of the world. In other semiarid plains, such massive vegetation replacements initiated a groundwater recharge and salt mobilization process that, after decades, raised regional water tables and salts to the surface, degrading agricultural and natural ecosystems. Indirect evidence suggests that this process (known as dryland salinity) began in the Chaco plains. Multiple approaches (deep soil profiles, geoelectric surveys and monitoring of groundwater salinity, level and isotopic composition) were combined to assess the dryland salinity status in one of the oldest and most active agricultural hotspots of the region, where isolated forest remnants occupy an extremely flat cultivated matrix. Full vadose moisture and chloride profiles from paired agriculture‐forest stands (17 profiles, six sites) revealed the following: a generalized onset of deep drainage with cultivation (32 to >87 mm year^?1), full leaching of native chloride pools (13.7 ± 2.5 kg m^?2) down to the water table after >40 years following clearing and differential groundwater table rises (0.7 to 2 m shallower water tables under agriculture than under neighbouring forests). Continuous level monitoring showed abrupt water table rises under annual crops (up to 2.6 m in 15 days) not seen under forests or pastures. Varying deep drainage rates and groundwater isotopic composition under agricultural plots suggest that these pulses are strongly modulated by crop choices and sequences. In contrast to other dryland salinity‐affected areas of the world, forest remnants in the study area (10–20% of the area) are not only surviving the observed hydrological shifts but also sustaining active salty groundwater transpirative discharge, as evidenced by continuous water table records. The overall impact of these forest remnants on lowering neighbouring water tables would be limited by the low hydraulic conductivity of the sediments. As highly cultivated areas of the Chaco evolve to new hydrological conditions of shallower saline water tables, innovative crop rotations that minimize recharge, enhance transpirative discharge and tolerate salinity will be needed. Copyright © 2016 John Wiley & Sons, Ltd.     

Selective drought-induced decline of pine species in southeastern Spain

The negative impacts of severe drought on the growth and vigor of tree species and their relationship with forest decline have not been properly evaluated taking into account the differential responses to such stress of trees, sites and species. We evaluated these responses by quantifying the changes in radial growth of plantations of four pine species (Pinus sylvestris, Pinus nigra, Pinus pinaster, Pinus halepensis) which showed distinct decline and defoliation levels in southeastern Spain. We used dendrochronological methods, defoliation records, linear mixed models of basal area increment and dynamic factor analysis to quantify the responses of trees at the species and individual scales to site conditions and drought stress. In the region a temperature rise and a decrease in spring precipitation have led to drier conditions during the late twentieth century characterized by severe droughts in the 1990s and 2000s. As expected, the defoliation levels and the reductions in basal area increment were higher in those species more vulnerable to drought-induced xylem embolism (P. sylvestris) than in those more resistant (P. halepensis). Species adapted to xeric conditions but with high growth rates, such as P. pinaster, were also vulnerable to drought-induced decline. The reduction in basal area increment and the defoliation events occurred after consecutive severe droughts. A decrease in spring precipitation, which is the main driver of radial growth, is the most plausible cause of recent forest decline. The sharp growth reduction and widespread defoliation of the most affected pine plantations of Scots pine make their future persistence in drought-prone sites unlikely under the forecasted warmer and drier conditions.     

The Early Palaeozoic high-grade metamorphism at the active continental margin of West Gondwana in the Andes (NW Argentina/N Chile)

The evolution of the Early Palaeozoic orogen of West Gondwana in the Cambrian to Ordovician basement of the Andes between ~18° and 32° S is investigated for pressure and temperature conditions and age of metamorphism. It is characterized by mid-crust temperatures commonly above the wet granite solidus (~650°C). Widespread felsic migmatite and rare granulite formed at pressures of ca 0.5?C0.7?GPa, locally 1.0?GPa. These rocks represent the deepest exhumed sections of the Early Palaeozoic crust. High pressure?Clow temperature rocks are absent. The crystallization ages, compiled from the literature in combination with new data, for near peak metamorphic conditions of these high-grade metamorphic rocks in NW Argentina and N Chile are ~530?C500?Ma and ~470?C420?Ma. Both age groups are spatially overlapping. Radiogenic isotope signatures (Sr, Nd, Pb) are used to characterize the Early Palaeozoic basement. The Pb and Sr isotope compositions of the Early Palaeozoic basement indicate mixing arrays between pre-Palaeozoic unradiogenic and radiogenic crust. Crustal residence ages (Sm?CNd T_DM) indicate a prominent event of crust formation around ~2?Ga, which is known continent-wide. This material was recycled during Midproterozoic and Early Palaeozoic orogenies without prominent additions of new crust present in the isotope record, i.e. accretion of compositional exotic material is absent.     

Analysis of a GPS Network Based on Functional Data Analysis

This paper demonstrates the usefulness of approaching the dynamic study of the precise positioning of a network of permanent global positioning system (GPS) stations through functional data analysis. The displacement data for each GPS station, obtained from observations of the global navigation satellite system, are a discrete sample of the positioning curve. The aim of this paper is to reconstruct the original functions in order to use them as functional data. In the method presented in this paper, the geodetic series are obtained first by processing the GPS data with respect to a reference station. Second, for each station, a cleaning process is applied to eliminate the values considered as outliers, and the missing values are imputed by using a Kalman filter. Finally, the original functions are reconstructed by using smoothing techniques and by evaluating several bases of functions. Moreover, these functions are treated with statistical techniques for functional data. This procedure is applied to the permanent stations of the south of the Iberian peninsula and the north of Africa (SPINA) network. The topocentric series: east, north and up are analysed. In the analysis of the positioning curves, there is observed a synchronized behaviour of the functions in those periods of time with important seismic activity. This behaviour also appears in the analysis of the second principal component of the East and up dimensions. Furthermore, the first two principal components of the East coordinate enable us to make a classification of the stations in the SPINA network. The classification made is consistent with the previous knowledge of the tectonic plates in the studied area.     

Middle Ordovician subduction of continental crust in the Scandinavian Caledonides: an example from Tjeliken,Seve Nappe Complex,Sweden

The Seve Nappe Complex of the Scandinavian Caledonides is thought to be derived from the distal passive margin of Baltica which collided with Laurentia in the Scandian Phase of the Caledonian Orogeny at 430–400 Ma. Parts of the Seve Nappe Complex were affected by pre-Scandian high- and ultrahigh-pressure metamorphism, in a tectonic framework that is still unclear, partly due to uncertainties about the exact timing. Previous age determinations yielded between ~ 505 and ~ 446 Ma, with a general trend of older ages in the North (Norrbotten) than in the South (Jämtland). New age determinations were performed on eclogite and garnet–phengite gneiss at Tjeliken in northern Jämtland. Thermodynamic modelling yielded peak metamorphic conditions of 25–27 kbar/680–760 °C for the garnet–phengite gneiss, similar to published peak metamorphic conditions of the eclogite (25–26 kbar/650–700 °C). Metamorphic rims of zircons from the garnet–phengite gneiss were dated using secondary ion mass spectrometry and yielded a concordia age of 458.9 ± 2.5 Ma. Lu–Hf garnet-whole rock dating yielded 458 ± 1.0 Ma for the eclogite. Garnet in the eclogite shows prograde major-element zoning and concentration of Lu in the cores, indicating that this age is related to garnet growth during pressure increase, i.e. subduction. The identical ages from both rock types, coinciding with published Sm–Nd ages from the eclogite, confirm subduction of the Seve Nappe Complex in Northern Jämtland during the Middle Ordovician in a fast subduction–exhumation cycle.     

Calibration and validation of rockfall modelling at regional scale: application along a roadway in Mallorca (Spain) and organization of its management

The Tramuntana range, in the northwestern sector of the island of Mallorca (Balearic Islands, Spain), is frequently affected by rockfalls which have caused significant damage, mainly along the road network. In this work, we present the procedure we have applied to calibrate and validate rockfall modelling in this region, using 103 cases of the available detailed rockfall inventory (630 rockfall events collected since the eighteenth century). We have exploited STONE (Guzzetti et al. 2002), a GIS-based rockfall simulation software which computes 2D and 3D rockfall trajectories starting from a DTM and maps of the dynamic rolling friction coefficient and of the normal and tangential energy restitution coefficients. The appropriate identification of these parameters determines the accuracy of the simulation. To calibrate them, we have selected 40 rockfalls along the range which include a wide variety of outcropping lithologies. Coefficients values have been changed in numerous attempts in order to select those where the extent and shape of the simulation matched the field mapping. Best results were summarized with the average statistical values for each parameter and for each geotechnical unit, determining that mode values represent more precisely the data. Initially, for the validation stage, 10 well-known rockfalls exploited in the calibration phase have been selected. Confidence tests have been applied to their modelling results taking into account not only the success but also the mistakes. The best accuracy is obtained when the rockfall has a preferential trajectory and worse results when the rockfall follows two or more trajectories. Additionally, the greater the rockfall runout length, the less precise the simulation is. We have further validated the calibrated parameters along the Ma-road (111 km), the main transportation corridor in the range, using 63 rockfall events that occurred during the past 18 years along the road. Of the rockfalls where source areas were properly identified, 81.5 % are well represented by STONE modelling, as the travel paths and the depositional areas are successfully ascertained. Results of the analysis have been used by the Road Maintenance Service of Mallorca to assess hazard and risk posed by rockfall at regional scale to design the road management plan.     

Backtracking drifting objects using surface currents from high-frequency (HF) radar technology

In this work, the benefits of high-frequency (HF) radar ocean observation technology for backtracking drifting objects are analysed. The HF radar performance is evaluated by comparison of trajectories between drifter buoys versus numerical simulations using a Lagrangian trajectory model. High-resolution currents measured by a coastal HF radar network combined with atmospheric fields provided by numerical models are used to backtrack the trajectory of two dataset of surface-drifting buoys: group I (with drogue) and group II (without drogue). A methodology based on optimization methods is applied to estimate the uncertainty in the trajectory simulations and to optimize the search area of the backtracked positions. The results show that, to backtrack the trajectory of the buoys in group II, both currents and wind fields were required. However, wind fields could be practically discarded when simulating the trajectories of group I. In this case, the optimal backtracked trajectories were obtained using only HF radar currents as forcing. Based on the radar availability data, two periods ranging between 8 and 10?h were selected to backtrack the buoy trajectories. The root mean squared error (RMSE) was found to be 1.01?km for group I and 0.82?km for group II. Taking into account these values, a search area was calculated using circles of RMSE radii, obtaining 3.2 and 2.11?km² for groups I and II, respectively. These results show the positive contribution of HF radar currents for backtracking drifting objects and demonstrate that these data combined with atmospheric models are of value to perform backtracking analysis of drifting objects.