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.     

Analysis of fibrous zeolites in the volcanic deposits of the Viterbo Province,Italy

Since Etrurian age, the Viterbo Province (Central Italy) is famous for its ancient towns carved out of local ignimbrite deposits which geologically represent the sedimentation of pumice-rich, volcano-related pyroclastic flows. Almost the entire study area is geologically characterized by a thick succession of ignimbrites, tephra fallouts and lava flow deposits locally subjected to zeolitization. Among zeolites, fibrous erionite represents a well-known health hazard and so this work aims at locating and quantifying the presence of fibrous zeolites in volcanic deposits of the Viterbo area, suggesting at the same time a standard operational procedure useful in other areas showing the same possible hazard. 41 samples collected in the Viterbo area were analyzed by X-ray powder diffraction, optical and electron microscopy. Fibrous zeolites were detected in five samples and quantified in amounts ranging from 0.35 to 1.66% vol. They generally occur in tetragonal or orthorhombic prismatic habit and showed chemical composition mainly consistent with K-phillipsite. Fibrous zeolites occurred with aspect ratios (mean = 6.3), comparable or lower than those of erionite fibers reported in previous studies and mean diameters >3 μm.     

LMDZ5B: the atmospheric component of the IPSL climate model with revisited parameterizations for clouds and convection

Based on a decade of research on cloud processes, a new version of the LMDZ atmospheric general circulation model has been developed that corresponds to a complete recasting of the parameterization of turbulence, convection and clouds. This LMDZ5B version includes a mass-flux representation of the thermal plumes or rolls of the convective boundary layer, coupled to a bi-Gaussian statistical cloud scheme, as well as a parameterization of the cold pools generated below cumulonimbus by re-evaporation of convective precipitation. The triggering and closure of deep convection are now controlled by lifting processes in the sub-cloud layer. An available lifting energy and lifting power are provided both by the thermal plumes and by the spread of cold pools. The individual parameterizations were carefully validated against the results of explicit high resolution simulations. Here we present the work done to go from those new concepts and developments to a full 3D atmospheric model, used in particular for climate change projections with the IPSL-CM5B coupled model. Based on a series of sensitivity experiments, we document the differences with the previous LMDZ5A version distinguishing the role of parameterization changes from that of model tuning. Improvements found previously in single-column simulations of case studies are confirmed in the 3D model: (1) the convective boundary layer and cumulus clouds are better represented and (2) the diurnal cycle of convective rainfall over continents is delayed by several hours, solving a longstanding problem in climate modeling. The variability of tropical rainfall is also larger in LMDZ5B at intraseasonal time-scales. Significant biases of the LMDZ5A model however remain, or are even sometimes amplified. The paper emphasizes the importance of parameterization improvements and model tuning in the frame of climate change studies as well as the new paradigm that represents the improvement of 3D climate models under the control of single-column case studies simulations.     

Control of deep convection by sub-cloud lifting processes: the ALP closure in the LMDZ5B general circulation model

Recently, a new conceptual framework for deep convection scheme triggering and closure has been developed and implemented in the LMDZ5B general circulation model, based on the idea that deep convection is controlled by sub-cloud lifting processes. Such processes include boundary-layer thermals and evaporatively-driven cold pools (wakes), which provide an available lifting energy that is compared to the convective inhibition to trigger deep convection, and an available lifting power (ALP) at cloud base, which is used to compute the convective mass flux assuming the updraft vertical velocity at the level of free convection. While the ALP closure was shown to delay the local hour of maximum precipitation over land in better agreement with observations, it results in an underestimation of the convection intensity over the tropical ocean both in the 1D and 3D configurations of the model. The specification of the updraft vertical velocity at the level of free convection appears to be a key aspect of the closure formulation, as it is weaker over tropical ocean than over land and weaker in moist mid-latitudes than semi-arid regions. We propose a formulation making this velocity increase with the level of free convection, so that the ALP closure is adapted to various environments. Cloud-resolving model simulations of observed oceanic and continental case studies are used to evaluate the representation of lifting processes and test the assumptions at the basis of the ALP closure formulation. Results favor closures based on the lifting power of sub-grid sub-cloud processes rather than those involving quasi-equilibrium with the large-scale environment. The new version of the model including boundary-layer thermals and cold pools coupled together with the deep convection scheme via the ALP closure significantly improves the representation of various observed case studies in 1D mode. It also substantially modifies precipitation patterns in the full 3D version of the model, including seasonal means, diurnal cycle and intraseasonal variability.     

Effectiveness of geophysical surveys for water wells relocation in overexploited aquifers (the example of Maggiore and Traversola Valleys,Northwestern Italy)

Aquifer overexploitation is a common problem in drinking water management. This will become more and more important given the general reduction of water resources. To overcome overexploitation, one of the most adopted solutions is the relocation of extraction wells. To establish the positions of new wells, a precise knowledge of the hydrogeological setting is required. Specific surveys are therefore necessary to obtain information over wide investigation zones. Geophysical methods, particularly electromagnetic and electrical, can be useful with this aim. In the present paper a case history on the combined use of Electric Resistivity Tomography and Time Domain Electromagnetic soundings is reported. Surveys have been performed within the Maggiore and Traversola Valleys, to investigate the uppermost part of the Quaternary deposits, hosting the near surface aquifer. The electromagnetic data have been inverted with a spatially constrained approach by assuming a quasi 1D model of the subsurface. Geophysical surveys allowed for depicting the depth and lateral continuity of the supposed aquifer level in the surveyed area up to a depth of about 200 m and proposing potential positions for well relocation.     

Monitoring Seismo-volcanic and Infrasonic Signals at Volcanoes: Mt. Etna Case Study

Volcanoes generate a broad range of seismo-volcanic and infrasonic signals, whose features and variations are often closely related to volcanic activity. The study of these signals is hence very useful in the monitoring and investigation of volcano dynamics. The analysis of seismo-volcanic and infrasonic signals requires specifically developed techniques due to their unique characteristics, which are generally quite distinct compared with tectonic and volcano-tectonic earthquakes. In this work, we describe analysis methods used to detect and locate seismo-volcanic and infrasonic signals at Mt. Etna. Volcanic tremor sources are located using a method based on spatial seismic amplitude distribution, assuming propagation in a homogeneous medium. The tremor source is found by calculating the goodness of the linear regression fit (R ²) of the log-linearized equation of the seismic amplitude decay with distance. The location method for long-period events is based on the joint computation of semblance and R ² values, and the location method of very long-period events is based on the application of radial semblance. Infrasonic events and tremor are located by semblance–brightness- and semblance-based methods, respectively. The techniques described here can also be applied to other volcanoes and do not require particular network geometries (such as arrays) but rather simple sparse networks. Using the source locations of all the considered signals, we were able to reconstruct the shallow plumbing system (above sea level) during 2011.     

Aeromechanic analysis of pyroclastic density currents past a building

An aeromechanic analysis of pyroclastic density currents (PDCs) past a building is carried out on the results of a computer simulation. The analysis shows that PDCs strongly interact with buildings, resulting in turbulent boundary layer separation and recirculation. These results could be used to better assess the hazard of PDCs impacting urban areas and be of service to civil protection authorities and urban planners who work in active volcanic areas.     

Parameterization of historical earthquakes in Switzerland

Macroseismic earthquake parameters of historical events have been reassessed in the framework of the update of the Earthquake Catalogue of Switzerland ECOS-09. The Bakun and Wentworth method (Bakun and Wentworth 1997) has been used to assess location, magnitude, and, when possible, focal depth. We apply a two-step procedure. Intensity attenuation is assessed first by fitting a model with a logarithmic and a linear term, using a set of 111 earthquakes. The magnitude range is 3 and 5.8. Then, intensity to magnitude relation is developed. A subset of the 111 events, all having an instrumental moment magnitude, was used to perform this intensity to magnitude calibration. Five final calibration strategies were developed based on different intensity calibration datasets, regionalized or non-regionalized models, and fixed or variable source depth. The final assessment of the macroseismic earthquake parameters is based on an expert judgment procedure, using the results derived from all five strategies, and taking into consideration the historical knowledge available for the particular earthquake. A bootstrap procedure has been applied to assess the uncertainty of parameters. Indicative lower and upper bounds of uncertainty are derived from distributions of location and magnitude for a number of events, obtained through bootstrap sampling of the intensity field and of the single intensity values. The final uncertainties are given in terms of parameter uncertainty classes already used in previous versions of the earthquake catalogue of Switzerland.     

Alluvial terraces on the Ionian coast of northern Calabria, southern Italy: Implications for tectonic and sea level controls

In this paper we present the results of an integrated geomorphological, pedological and stratigraphical study carried out along the Ionian coast of northern Calabria (southern Italy). This area is characterised by the occurrence of five orders of alluvial terraces that are striking features of the landscape, where large and steep catchments debouch from the mountain front to the hilly coastal belt.Field investigations indicate that the deposits of all five terraces are suggestive of shallow gravel-bed braided streams.On the basis of the age of the Pleistocene substratum and morphostratigraphic correlation with marine terraces cropping out in the nearby areas, each order has been associated to specific marine oxygen isotope stages.Consequently, we focused on the interplay of allocyclic factors influencing stream aggradation/degradation. Soil features and other climatic proxies suggest that climate didn't play an important role with respect to tectonic and base-level changes in controlling fluvial dynamics.In particular, we recognised that during the middle Pleistocene the study area experienced a period of subaerial landscape modelling, as suggested by the thick and complex alluvial sequence of the highest terrace (T1). The onset of regional uplift marks a change in the geomorphic scenario, with tectonic and eustatically driven changes in base-level working together in causing switches in fluvial aggradational/erosional phases (T2–T5 terraces). Because of the uplift, river dissection occurred during phases of sea level fall, whereas aggradation phases occurred during periods of climate amelioration (sea level rise) just before highstands were attained.As a consequence, the stepped terraces in the study area reflect the interplay between tectonics (uplift) and sea level changes, in which terraces define episodes of relative sea level fall during the late Quaternary.