Magnetic anisotropy produced by magma flow: theoretical model and experimental data from Ferrar dolerite sills (Antarctica)

Volcanic rocks forming sills, dykes or lava flows may display a magnetic anisotropy derived from the viscous flow during their emplacement. We model a sill as a steady-state flow of a Bingham fluid, driven by a pressure gradient in a horizontal conduit. The magma velocity as a function of depth is calculated from the motion and constitutive equations. Vorticity and strain rate are determined for a reference system moving with the fluid. The angular velocity and the orientation of an ellipsoidal magnetic grain immersed in the fluid are calculated as functions of time or strain. Magnetic susceptibility is then calculated for a large number of grains with a uniform distribution of initial orientations. It is shown that the magnetic lineation oscillates in the vertical plane through the magma flow direction, and that the magnetic foliation plane changes periodically from horizontal to vertical. The results are compared with the magnetic fabric of Ferrar dolerite sills (Victoria Land, East Antarctica) derived from low-field susceptibility measurements.     

The effect of crystallization on the rheology and dynamics of lava flows

The dynamics of a lava flow is studied by a two-dimensional model describing a viscous fluid with Bingham rheology, flowing down a slope. The temperature in the flow is calculated assuming that heat is transferred through the plug by conduction and is lost by radiation to the atmosphere at the top of the flow. Taken into account is that the increasing crystallization takes place in the flow as a consequence of cooling. The lava viscosity and yield stress are expressed as a function of crystallization degree as well as of temperature: in particular it is assumed that yield stress reaches a maximum value above the solidus temperature, according to experimental data. Dynamical variables, such as velocity and thickness of the flow, are calculated for different values of the maximum crystallization degree and the flow rate. The model shows how the lava flow dynamics is affected by cooling and crystallization. The cooling of the flow is controlled by the increase of yield stress, which produces a thicker plug and makes the heat loss slower. The increasing crystallization has two opposing effects on viscosity: it produces an increase of viscosity, but at the same time produces an increase of yield stress and hence reduces the heat loss and keeps the internal temperature high. As a consequence, lava flows are significantly affected by the dependence of yield stress on temperature and scarcely by the maximum crystallization degree.     

The historic and artistic heritage facing the earthquake risk: The Italian case

The problem of protection against earthquakes in Italy is made difficult by the presence of a huge historic and artistic heritage. Such a heritage is mainly made up of ancient buildings and monuments situated in the urban centres, which are densely distributed throughout Italy. Therefore, very complex problems are met in the determination of the value parameter concurring with the determination of seismic risk, in addition to hazard and vulnerability. An indication of the monetary value of a building is not sufficient as far as the cultural heritage is concerned: different criteria are necessary in order to distinguish which are the strategic buildings. If we consider that there are more than 2000 museums in Italy, most of which are placed inside historic buildings, it appears that museums should receive the highest priority in future initiatives for seismic rehabilitation.     

A viscoelastic shear zone model of compressional and extensional plate boundaries

A model is proposed describing the mechanical evolution of a shear zone along compressional and extensional plate boundaries, subject to constant strain rate. The shear zones are assumed as viscoelastic with Maxwell rheology and with elastic and rheological parameters depending on temperature and petrology. Stress and strain are computed as functions of time and depth. For both kinds of boundaries the model reproduces the existence of a shallow seismogenic zone, characterized by a stress concentration. The thickness of the seismogenic layer is evaluated considering the variations of shear stress and frictional strength on faults embedded in the shear zone. Assuming that a fault dislocation takes place, the brittle-ductile transition is assumed to occur at the depth at which the time derivative of total shear stress changes from positive to negative values. The effects of different strain rates and geothermal gradients on the depth of the brittle-ductile transition are studied. The model predictions are consistent with values inferred from seismicity data of different boundary zones.     

Shallow earthquakes in a viscoelastic shear zone with depth-dependent friction and rheology

Summary. Most crustal earthquakes of the world are observed to occur within a seismogenic layer which extends from the Earth's surface to a depth of a few tens of kilometres at most. A model is proposed in which the shear zone along a transcurrent plate margin is represented as a viscoelastic medium with depth-dependent power-law rheology. A frictional resistance linearly increasing with depth is assumed on a vertical transcurrent fault within the shear zone. Such a model is able to reproduce a continuous transition from the brittle behaviour of the upper crust to the ductile behaviour at depth. Assuming that the shear zone is subjected to a constant strain rate from the opposite motions of the two adjacent plates, it is found that there exists a maximum depth H below which tectonic stress can never reach the frictional threshold: this may be identified as the maximum depth of earthquake nucleation. The value of H is consistent with observations for plausible values of the model parameters. The stress evolution in the shear zone is calculated in the linear approximation of the constitutive equation. A change in rigidity with depth, which is also introduced in the model, may reproduce the high vertical gradient of shear stress, which has been measured across the San Andreas fault, and the fact that most earthquakes are nucleated at some depth in the seismogenic layer. A crack which drops the ambient stress to the dynamic frictional level is then introduced in the model. To this aim, a crack solution is employed without a stress singularity at its edges, which is compatible with a frictional stress threshold criterion for fracture. A constraint on the vertical friction gradient is obtained if such cracks are assumed to be entirely confined within the seismogenic layer.     

Seismic-induced rockfalls and landslide dam following the October 30, 2016 earthquake in Central Italy

On October 30, 2016, a seismic event and its aftershocks produced diffuse landslides along the SP 209 road in the Nera River Gorge (Central Italy). Due to the steep slopes and the outcropping of highly fractured and bedded limestone, several rockfalls were triggered, of which the main event occurred on the slope of Mount Sasso Pizzuto. The seismic shock acted on a rock wedge that, after an initial slide, developed into a rockfall. The debris accumulation blocked the SP 209 road and dammed the Nera River, forming a small lake. The river discharge was around 3.6 m³/s; the water overtopped the dam and flooded the road. By a preliminary topographic survey, we estimated that the debris accumulation covers an area of about 16,500 m², while the volume is around 70,000 m³. The maximum volume occupied by the pre-existing talus mobilized by the rockfall is about 20% of the total volume. Besides blocking the road, the rockfall damaged a bridge severely, while, downstream of the dam, the water flow caused erosion of a road embankment. A rockfall protection gallery, a few hundred meters downstream of the dam, was damaged during the event. Other elastic nets and rigid barriers were not sufficient to protect the road from single-block rockfalls, with volumes around 1–2 m³. Considering the geological and geomorphological conditions, as well as the high seismicity and the socioeconomic importance of the area, a review of the entire rockfall protection systems is required to ensure protection of critical infrastructure and local communities.     

Effects of geological complexities on coseismic displacement: hints from 2D numerical modelling

By means of a 2D finite‐element procedure, we tested how heterogeneities at the scale of seismogenic fault affect the displacement. We defined one or more slip distributions for two typical normal‐faulting earthquakes in the Central Apennines, computed the displacement occurred within different structures including lateral heterogeneities, and compared the different displacement profiles to isolate the effect of the crustal structure. To understand at what magnitude the heterogeneities affect the observation significantly, we compared the predicted coseismic displacement with GPS and DInSAR data for the Colfiorito 1997 earthquake. We find that heterogeneities significantly affect the observable horizontal coseismic displacement for the larger magnitudes, whereas for smaller quakes, they affect horizontal displacement close to the fault trace only.     

Viscous Newtonian laminar flow in a rectangular channel: application to Etna lava flows

 We introduce a 3D model for near-vent channelized lava flows. We assume the lava to be an isothermal Newtonian liquid flowing in a rectangular channel down a constant slope. The flow velocity is calculated with an analytical steady-state solution of the Navier-Stokes equation. The surface velocity and the flow rate are calculated as functions of the flow thickness for different flow widths, and the results are compared with those of a 2D model. For typical Etna lava flow parameters, the influence of levees on the flow dynamics is significant when the flow width is less than 25 m. The model predicts the volume flow rate corresponding to the surface velocity, taking into account that both depend on flow thickness. The effusion rate is a critical parameter to evaluate lava flow hazard. We propose a model to calculate the effusion rate given the lava flow width, the topograhic slope, the lava density, the surface flow velocity, and either the lava viscosity or the flow thickness. Received: 20 January 1998 / Accepted: 8 January 1999