Long-Term RST Analysis of Anomalous TIR Sequences in Relation with Earthquakes Occurred in Greece in the Period 2004–2013

Real-time integration of multi-parametric observations is expected to accelerate the process toward improved, and operationally more effective, systems for time-Dependent Assessment of Seismic Hazard (t-DASH) and earthquake short-term (from days to weeks) forecast. However, a very preliminary step in this direction is the identification of those parameters (chemical, physical, biological, etc.) whose anomalous variations can be, to some extent, associated with the complex process of preparation for major earthquakes. In this paper one of these parameters (the Earth’s emitted radiation in the Thermal InfraRed spectral region) is considered for its possible correlation with M ≥ 4 earthquakes occurred in Greece in between 2004 and 2013. The Robust Satellite Technique (RST) data analysis approach and Robust Estimator of TIR Anomalies (RETIRA) index were used to preliminarily define, and then to identify, significant sequences of TIR anomalies (SSTAs) in 10 years (2004–2013) of daily TIR images acquired by the Spinning Enhanced Visible and Infrared Imager on board the Meteosat Second Generation satellite. Taking into account the physical models proposed for justifying the existence of a correlation among TIR anomalies and earthquake occurrences, specific validation rules (in line with the ones used by the Collaboratory for the Study of Earthquake Predictability—CSEP—Project) have been defined to drive a retrospective correlation analysis process. The analysis shows that more than 93 % of all identified SSTAs occur in the prefixed space–time window around (M ≥ 4) earthquake's time and location of occurrence with a false positive rate smaller than 7 %. Molchan error diagram analysis shows that such a correlation is far to be achievable by chance notwithstanding the huge amount of missed events due to frequent space/time data gaps produced by the presence of clouds over the scene. Achieved results, and particularly the very low rate of false positives registered on a so long testing period, seems already sufficient (at least) to qualify TIR anomalies (identified by RST approach and RETIRA index) among the parameters to be considered in the framework of a multi-parametric approach to t-DASH.     

Vegetation turnover in a braided river: frequency and effectiveness of floods of different magnitude

This work addresses the temporal dynamics of riparian vegetation in large braided rivers, exploring the relationship between vegetation erosion and flood magnitude. In particular, it investigates the existence of a threshold discharge, or a range of discharges, above which erosion of vegetated patches within the channel occurs. The research was conducted on a 14 km long reach of the Tagliamento River, a braided river in north‐eastern Italy. Ten sets of aerial photographs were used to investigate vegetation dynamics in the period 1954–2011. By using different geographic information system (GIS) procedures, three aspects of geomorphic‐vegetation dynamics and interactions were addressed: (i) long‐term (1954–2011) channel evolution and vegetation dynamics; (ii) the relationship between vegetation erosion/establishment and flow regime; (iii) vegetation turnover, in the period 1986–2011. Results show that vegetation turnover is remarkably rapid in the study reach with 50% of in‐channel vegetation persisting for less than 5–6 years and only 10% of vegetation persisting for more than 18–19 years. The analysis shows that significant vegetation erosion is determined by relatively frequent floods, i.e. floods with a recurrence interval of c. 1–2.5 years, although some differences exist between sub‐reaches with different densities of vegetation cover. These findings suggest that the erosion of riparian vegetation in braided rivers may not be controlled solely by very large floods, as is the case for lower energy gravel‐bed rivers. Besides flow regime, other factors seem to play a significant role for in‐channel vegetation cover over long time spans. In particular, erosion of marginal vegetation, which supplies large wood elements to the channel, increased notably over the study period and was an important factor for in‐channel vegetation trends. Copyright © 2014 John Wiley & Sons, Ltd.     

Carbon dioxide degassing and thermal energy release in the Monte Amiata volcanic-geothermal area (Italy)

The quaternary volcanic complex of Mount Amiata is located in southern Tuscany (Italy) and represents the most recent manifestation of the Tuscan Magmatic Province. The region is characterised by a large thermal anomaly and by the presence of numerous CO₂-rich gas emissions and geothermal features, mainly located at the periphery of the volcanic complex. Two geothermal systems are located, at increasing depths, in the carbonate and metamorphic formations beneath the volcanic complex. The shallow volcanic aquifer is separated from the deep geothermal systems by a low permeability unit (Ligurian Unit). A measured CO₂ discharge through soils of 1.8 × 10⁹ mol a⁻¹ shows that large amounts of CO₂ move from the deep reservoir to the surface. A large range in δ¹³C_TDIC (−21.07 to +3.65) characterises the waters circulating in the aquifers of the region and the mass and isotopic balance of TDIC allows distinguishing a discharge of 0.3 × 10⁹ mol a⁻¹ of deeply sourced CO₂ in spring waters. The total natural CO₂ discharge (2.1 × 10⁹ mol a⁻¹) is slightly less than minimum CO₂ output estimated by an indirect method (2.8 × 10⁹ mol a⁻¹), but present-day release of 5.8 × 10⁹ mol a⁻¹ CO₂ from deep geothermal wells may have reduced natural CO₂ discharge. The heat transported by groundwater, computed considering the increase in temperature from the infiltration area to the discharge from springs, is of the same order of magnitude, or higher, than the regional conductive heat flow (>200 mW m⁻²) and reaches extremely high values (up to 2700 mW m⁻²) in the north-eastern part of the study area. Heat transfer occurs mainly by conductive heating in the volcanic aquifer and by uprising gas and vapor along fault zones and in those areas where low permeability cover is lacking. The comparison of CO₂ flux, heat flow and geological setting shows that near surface geology and hydrogeological setting play a central role in determining CO₂ degassing and heat transfer patterns.     

Assessing the results of scenarios of climate and land use changes on the hydrology of an Italian catchment: modelling study

Hydrological models are recognized as valid scientific tools to study water quantity and quality and provide support for the integrated management and planning of water resources at different scales. In common with many catchments in the Mediterranean, the study catchment has many problems such as the increasing gap between water demand and supply, water quality deterioration, scarcity of available data, lack of measurements and specific information. The application of hydrological models to investigate hydrological processes in this type of catchments is of particular relevance for water planning strategies to address the possible impact of climate and land use changes on water resources. The distributed catchment scale model (DiCaSM) was selected to study the impact of climate and land use changes on the hydrological cycle and the water balance components in the Apulia region, southern Italy, specifically in the Candelaro catchment (1780 km²). The results obtained from this investigation proved the ability of DiCaSM to quantify the different components of the catchment water balance and to successfully simulate the stream flows. In addition, the model was run with the climate change scenarios for southern Italy, i.e. reduced winter rainfall by 5–10%, reduced summer rainfall by 15–20%, winter temperature rise by 1·25–1·5 °C and summer temperature rise by 1·5–1·75 °C. The results indicated that by 2050, groundwater recharge in the Candelaro catchment would decrease by 21–31% and stream flows by 16–23%. The model results also showed that the projected durum wheat yield up to 2050 is likely to decrease between 2·2% and 10·4% due to the future reduction in rainfall and increase in temperature. In the current study, the reliability of the DiCaSM was assessed when applied to the Candelaro catchment; those parameters that may cause uncertainty in model output were investigated using a generalized likelihood uncertainty estimation (GLUE) methodology. The results showed that DiCaSM provided a small level of uncertainty and subsequently, a higher confidence level. Copyright © 2010 John Wiley & Sons, Ltd.     

On characterizing the temporal dominance patterns of model parameters and processes

Diagnostic analyses of hydrological models intend to improve the understanding of how processes and their dynamics are represented in models. Temporal patterns of parameter dominance could be precisely characterized with a temporally resolved parameter sensitivity analysis. In this way, the discharge conditions are characterized, that lead to a parameter dominance in the model. To achieve this, the analysis of temporal dynamics in parameter sensitivity is enhanced by including additional information in a three‐tiered framework on different aggregation levels. Firstly, temporal dynamics of parameter sensitivity provide daily time series of their sensitivities to detect variations in the dominance of model parameters. Secondly, the daily sensitivities are related to the flow duration curve (FDC) to emphasize high sensitivities of model parameters in relation to specific discharge magnitudes. Thirdly, parameter sensitivities are monthly averaged separately for five segments of the FDC to detect typical patterns of parameter dominances for different discharge magnitudes. The three methodical steps are applied on two contrasting catchments (upland and lowland catchment) to demonstrate how the temporal patterns of parameter dynamics represent different hydrological regimes. The discharge dynamic in the lowland catchment is controlled by groundwater parameters for all discharge magnitudes. In contrast, different processes are relevant in the upland catchment, because the dominances of parameters from fast and slow runoff components in the upland catchment are changing over the year for the different discharge magnitudes. The joined interpretation of these three diagnostic steps provides deeper insights of how model parameters represent hydrological dynamics in models for different discharge magnitudes. Thus, this diagnostic framework leads to a better characterization of model parameters and their temporal dynamics and helps to understand the process behaviour in hydrological models. Copyright © 2015 John Wiley & Sons, Ltd.     

Water‐bottom multiple attenuation by Kirchhoff extrapolation

Despite being less general than 3D surface‐related multiple elimination (3D‐SRME), multiple prediction based on wavefield extrapolation can still be of interest, because it is less CPU and I/O demanding than 3D‐SRME and moreover it does not require any prior data regularization. Here we propose a fast implementation of water‐bottom multiple prediction that uses the Kirchhoff formulation of wavefield extrapolation. With wavefield extrapolation multiple prediction is usually obtained through the cascade of two extrapolation steps. Actually by applying the Fermat’s principle (i.e., minimum reflection traveltime) we show that the cascade of two operators can be replaced by a single approximated extrapolation step. The approximation holds as long as the water bottom is not too complex. Indeed the proposed approach has proved to work well on synthetic and field data when the water bottom is such that wavefront triplications are negligible, as happens in many practical situations.     

Chlorinated compounds and polybrominated diphenyl ethers (PBDEs) in mussels (Mytilus galloprovincialis) collected from Apulia Region coasts

This project was carried out to assess the levels and spatial distribution of organochlorine compounds in the coastal marine environment, using mussels as bioindicators to evaluate the coastal water quality. Levels of polychlorobiphenils (PCB), chlorinated pesticides (DDT isomers, HCH isomers, Aldrin, Dieldrin, alfa-Endosulfan, Hexachlorobenzene, Pentachlorobenzene) and polybrominated diphenyl ethers (PBDEs) were determined in tissues from mussels (Mytilus galloprovincialis) collected along the Apulia Region coasts (Mediterranean Sea). Results indicate that contamination by organochlorine compounds is higher in mussels sampled in the Ionian Sea than in those from the Adriatic Sea, with PCB levels up to seven times higher in mussels from Ionian than from the Adriatic Sea. Although PCB levels were above the maximum values indicated by both European Community (EC) and National regulation in several sample sites, the PCB concentrations were particularly high in some stations, suggesting that these locations require a much specific attention. Conversely, results on the mussel contamination by PBDEs highlight their ubiquitous environmental distribution, and underline the need to establish the maximum level for these compounds in foodstuff, according to European Regulations.     

Landslide HotSpot Mapping by means of Persistent Scatterer Interferometry

Landslide detection and mapping represent fundamental requirements for every hazard and risk evaluation and consequent improvement of the management strategies for such natural hazards. Optical and radar remote sensing can be used to observe landslide-induced ground deformation, ranging from regional to local scales. This work presents a methodology called Landslide HotSpot Mapping; this approach integrates cartographic, thematic and optical data with Persistent Scatterer Interferometry for the identification of extremely slow to very slow moving landslides, and for the evaluation of their state of activity and intensity. This methodology scans wide areas to detect hotspots, which are narrow unstable zones characterized by higher landslide hazard. To these hotspots, priority has to be given when planning field surveys and in situ validation campaigns, so that field work time and effort can be optimized and significantly reduced. The approach is tested in Central Calabria, over a 4,470?km² area located in southern Italy. ENVISAT ascending images acquired between 2003 and 2009 and processed with the Persistent Scatterer Pairs (PSP) technique are used to analyse deformation patterns. Combining conventional photo-interpretation with the analysis of PSP data, 64 new landslides are identified and the spatial (boundaries) and temporal (activity) information of 980 pre-mapped phenomena (23.6% of updated inventory) are updated. 1,012 active (continuous or reactivated) landslides are identified and 4 hotspot areas selected: San Fili, Rende, Lago, Catanzaro. Urgent field checks have to be organized for these hotspots to validate the satellite-based observations and to design appropriate mitigation measures to reduce impacts on the elements at risk.     

Middle and late Pleistocene glaciations in the Campo Felice Basin (central Apennines, Italy)

The present paper refers to research conducted in the tectonic-karst depression of Campo Felice in the central Apennines (Italy), in which glacial, alluvial and lacustrine sediments have been preserved. Stratigraphic interpretations of sediments underlying the Campo Felice Plain are based on evidence obtained from nine continuous-core boreholes. The boreholes reach a depth of 120 m and provide evidence of five sedimentation cycles separated by erosion surfaces. Each cycle is interpreted as an initial response to a mainly warm stage, characterized by low-energy alluvial and colluvial deposition, pedogenesis, and limited episodes of marsh formation. In turn, a mainly cold stage follows during which a lake formed, and glaciers developed and expanded, leading to deposition of glacial and fluvioglacial deposits. The chronological framework is established by eleven accelerator mass spectrometer ¹⁴C ages and three ³⁹Ar-⁴⁰Ar ages on leucites from interbedded tephra layers. These age determinations indicate five glacial advances that respectively occurred during marine oxygen isotope stages 2, 3-4, 6, 10 and 14.