Magma chambers and localization of deformation during thrusting

Terra Nova, 22, 390–395, 2010 Abstract We present the results of coupled analogue and numerical models that provide new insights into the relationships between volcanoes and thrusts. The effects of both upper‐crustal magma chambers and the load of volcanoes on the geometry of thrust systems were investigated. Analogue modelling points to a strong influence exerted by a magma chamber on thrust geometry, which, as suggested by the numerical models used to rationalize these results, is related to the stress redistribution around the weak heterogeneity. The low‐viscosity body below a volcanic edifice localizes compressional deformation and causes a curvature of the thrusts towards the magma chamber, opposite to the direction of tectonic transport. In these conditions, the volcanic load has a negligible effect on the structural geometry. These results are in contrast with those of previous studies, where intrusions or the load of major volcanoes generated a curvature of the thrusts away from volcanic edifices in the direction of tectonic transport.     

Integrating computer simulations and experiments to study explosive volcanic eruptions

Explosive eruptions cannot be directly measured because they are inaccessible, which is the reason why computer simulations are developing more and more in volcanology. Calculations of volcanic hazard take substantial advantage of the simulations, particularly when they are validated experimentally. Here, an integrated approach is proposed, which combines 3D simulations and large-scale experiments, in order to provide an effective tool that helps detailing actual explosive eruptions.     

Shallow landslides in deeply weathered slates of the Sierra de Collcerola (Catalonian Coastal Range, Spain)

Stability conditions in an area located NW of Barcelona (Spain) are discussed. Here, several mass movements were observed, mainly affecting weathered Paleozoic slates. Many of these failures involved slopes cut along recent infrastructures: debris flows, wedge and plane failures, generally surficial, occurred more frequently. After a detailed geological and geomorphologic survey, geomechanic characterization was carried out, according to RMR and SMR classifications. This rating gave a prediction of slope behaviour, in fairly good agreement with the real observed one.

Stability numerical analysis was carried out for the main cut slopes, based upon the Limit Equilibrium Method. First of all, the deterministic factor of safety was computed using the mean values of parameters. After that, a simulation technique based upon the Monte Carlo Method was applied in order to obtain factor of safety distributions. The probability of failure was estimated as P(F<1).

Finally, results from deterministic and probabilistic approaches were compared. The effectiveness of different possible remedial measures was highlighted by means of a sensitivity analysis, which showed that the more important parameters in the study area are the geometrical ones (height, slope and failure plane angles). The final technical solutions adopted are briefly outlined.     

Liquefaction of soil in the Emilia-Romagna region after the 2012 Northern Italy earthquake sequence

At the end of May 2012, the Po plain region in northern Italy was shaken by a long sequence of seismic events. The 2012 Northern Italy earthquake sequence counted two mainshocks, about 1,600 aftershocks and lasted for several weeks. Although the mainshocks, which occurred on May 20 and May 29, 2012, registered a moment magnitude of 5.9 and 5.8, respectively, these two events caused widespread soil liquefaction and substantial damages to the built environment. This paper reports lessons learnt from a field investigation conducted in the areas affected by the earthquake sequence. Based on the field observations, it was concluded that despite the relatively low magnitudes of the shocks, most of the damages occurred as a consequence of liquefaction phenomena and/or absence of retrofitting of historical structures. The latter comprise churches, tower bells, towers, castles and fortresses. It was found that the occurrence of liquefaction was mainly associated with the presence of saturated alluvial soil deposits which were characterised by high liquefaction susceptibility. It was noted that these highly liquefiable soils were mainly located in proximity of ancient river courses that were artificially diverted in the eighteenth century to mitigate flooding and other hydrological risks.     

Exploiting historical rainfall and landslide data in a spatial database for the derivation of critical rainfall thresholds

Critical rainfall thresholds for landslides are powerful tools for preventing landslide hazard. The thresholds are commonly estimated empirically starting from rainfall events that triggered landslides in the past. The creation of the appropriate rainfall–landslide database is one of the main efforts in this approach. In fact, an accurate agreement between the landslide and rainfall information, in terms of location and timing, is essential in order to correctly estimate the rainfall–landslide relationships. A further issue is taking into account the average moisture conditions prior the triggering event, which reasonably may be crucial in determining the sufficient amount of precipitation. In this context, the aim of this paper is exploiting historical landslide and rainfall data in a spatial database for the derivation of critical rainfall thresholds for landslide occurrence in Sicily, southern Italy. The hourly rainfall events that caused landslides occurred in the twentieth century were specifically identified and reconstructed. A procedure was proposed to automatically convert rain guages charts recorded on paper tape into digital format and then to provide the cumulative rainfall hyetograph in digital format. This procedure is based on a segmentation followed by signal recognition techniques which allow to digitalize and to recognize the hyetograph automatically. The role of rainfall prior to the landslide events was taken into account by including in the analysis the rainfall occurred 5, 15 and 30 days before each landslide. Finally, cumulated rainfall duration thresholds for different exceedance probability levels were determined. The obtained thresholds resulted in agreement with the regional curves proposed by other authors for the same area; antecedent rainfall turned out to be particularly important in triggering landslides.     

Topographical controls on small‐volume pyroclastic flows

Emplacement of small‐volume (<0·1 km³) pyroclastic flows is significantly influenced by topography. The Arico ignimbrite on Tenerife (Canary Islands) is a characteristic small‐volume pyroclastic flow deposit emplaced on high relief topography. The pyroclastic flow flowed down pre‐existing valleys on the southern slopes of the island. In proximal areas deep (up to 100 m) valleys acted as efficient conduits for the pyroclastic flow, which was mostly channelled; in this particular area the ignimbrite corresponds to a homogeneous, moderately welded deposit, consisting of flattened pumices in an abundant ashy matrix with a relatively low lithic fragment content. In intermediate zones significant changes occur in the steepness of the slope and, although still channelled, here the pyroclastic flow was influenced by hydraulic jumps. In this area, two different units can be clearly distinguished in the ignimbrite: the lower unit is composed of a lithic‐rich ground‐layer deposit that formed at the turbulent, highly concentrated head of the flow; the upper unit consists of a well welded pumice‐rich deposit that occasionally reveals a basal layer formed by shearing with the lower part. This division into two units is maintained as far as distal areas near the present‐day coastline, where the slope is very gentle or null and the ignimbrite is not channelled. The ground layer is not found in distal areas. The ignimbrite here only consists of the upper unit, which is occasionally repeated due to a surging process provoked by the lower flow speed, as the pyroclastic flow spread out of the channelled zone. A theoretical model on how topography controlled the deposition of the Arico ignimbrite is derived by interpreting the observed lithological and sedimentological variations in terms of changes in topography and bedrock morphology. This new model is of general applicability and will help to explain other deposits of similar characteristics.     

The performance of churches in the 2012 Emilia earthquakes

In this study the damage suffered by churches during the 2012 Emilia seismic sequence in Italy is analysed, based on surveys and inspections carried out in the area. Similarly to what was observed after other Italian earthquakes, the damage to churches was severe. However, the Emilia churches present some characteristic features such as the use of unreinforced clay brick masonry. In order to appropriately address the performance of this class of buildings, typical architectural layouts and construction techniques are described. Such techniques are interpreted also in the light of the local seismic catalogue. Fifty churches are then selected and their damage is studied, with reference to typical local-collapse mechanisms of different macro-elements. The study highlights that the damage is often concentrated at the top section of the façade, in the clerestory walls, in the vaults and in the bell towers. Structural analyses are performed to explain some of the observations. The overturning of the top section of the façade is analytically addressed, modelling the friction interlocking. With reference to the case study of San Francesco in Mirandola, non-linear static and dynamic analyses allow us to correlate the directionality of damage to the higher seismic demand along the NS direction, to point out the negligible role of the large vertical component of ground motion and to emphasise the relevance of the buttresses for the seismic response of the façade.     

Effects of air temperature on the cyclic behavior of elastomeric seismic isolators

The mechanical properties of elastomers can change significantly due to air temperature variations. In particular, prolonged exposure to subzero temperatures can result in rubber crystallization, with a considerable increase in the shear stiffness of the material. As a result, the seismic response of structures with elastomeric isolators can be strongly influenced by air temperature. Current seismic codes, indeed, require an upper and lower bound analysis, using suitable modification factors, to account for the changes in the cyclic behavior of elastomeric isolators due to air temperature variations. In this study, the sensitivity of the cyclic behavior of elastomeric isolators to air temperature variations is investigated based on the experimental results of an extensive test program on six different elastomeric compounds for seismic isolators, characterized by a shear modulus ranging from 0.5 to 1.2 MPa at 100% shear strain and 20°C. The cyclic tests have been performed on small-size specimens, subjected to shear strain amplitudes and frequency of loading typical for elastomeric seismic isolators, at seven different air temperatures, ranging from 40 to −20°C. The effects of rubber crystallization due to prolonged exposure to low-temperatures have been also investigated. A finite element model for the evaluation of the temperature contour map inside a full-size elastomeric isolator exposed to low air temperatures has been also developed. In the paper, the experimental outcomes are compared with the modification factors provided by the current seismic codes to account for the temperature effects on the mechanical properties of elastomeric isolators.     

First study of fault trench stratigraphy at Mt. Etna volcano, Southern Italy: understanding Holocene surface faulting along the Moscarello fault

Paleoseismology, the study of past earthquakes based on their geological record in the stratigraphy and landscape, is a successful newly developing field of research. The application of fault trench studies in volcanic environments is one of the youngest branches of paleoseismology. In this paper, we present the results of the first exploratory trenches excavated at Mt. Etna in Sicily, the largest European volcano. Modern surface faulting at Etna is a very well known feature, which poses significant hazard to the local community, both in terms of ground displacement of essential lifelines and ground shacking from frequent damaging earthquakes. However, while the geomorphology and the seismicity of the active fault in the Etna region consistently show very high rates of tectonic activity, the Holocene cumulative throw and slip-rates, along with the nature (coseismic vs. creeping fault slip), dimension and timing of the displacement events, are still poorly constrained. For this purpose, we selected as a sample area the Moscarello fault, one of the most outstanding segment of the Timpe system of active normal faults in the volcano’s lower eastern flank. Displaced landforms and volcanic units at the Fondo Macchia basin, in the central sector of this fault, indicate some hundreds of meters of vertical offset in the last ca. 80 kyr, with a long-term slip-rate substantially higher than 1.5–2.0 mm/yr. According to the historical sources and instrumental observations, the Moscarello fault ruptured four times in the last 150 years during shallow (H < 5 km) and moderate magnitude (M < 4.8) earthquakes. These events were associated with severe damage in a narrow epicentral area (macroseismic intensities up to the IX–X grade of the MSK scale) and extensive surface faulting (end-to-end rupture length up to 6 km, vertical offsets up to 90 cm). This clearly indicates very high modern rates of deformation along this fault. We conducted trench investigations at the Fondo Macchia site, in a point where eyewitnesses observed ca. 20 cm of coseismic vertical displacement after the April 21, 1971, M_s=3.7, earthquake. The excavated sections provided direct stratigraphic evidence for a vertical slip-rate of 1.4–2.7 mm/yr in the last ca. 6 kyr. This should be regarded as a minimum slip-rate for the central section of the fault. We explored a single scarp at a single site, while we know from recent historical observations that several parallel scarps may rupture coseismically at Fondo Macchia. Thus, the relevant deformation rate documented for the modern period might be likely extended back in the past to a time-span of some thousands of years at least. As expected, for such a volcanic environment, the activity rates of the Moscarello fault are also significantly higher than for the Apennines normal faults, typically showing slip-rates lower than 1 mm/yr. The agriculturally reworked trench hangingwall stratigraphy did not allow to recognize individual displacement events. Nevertheless, the sedimentary structures observed in the trench footwall strongly suggest that, as for the last 150–200 years of detailed historical record, fault behavior at Fondo Macchia is governed by coseismic surface displacement rather than fault creep. This research confirms that paleoseismology techniques can be effectively applied also in active volcanic environments, typically characterized by rheology and, consequently, seismicity and fault dynamics very different from those of other tectonic environments in which paleoseismology has been firstly developed and is today extensively applied.