Modeling of Rayleigh waves dispersion in the Sannio region (Southern Italy) from seismic active refraction data

Short period surface waves, recorded during a seismic refractionsurvey in the Sannio region (Southern Italy), have been modeled to infera shallow velocity model for the area. Based on the decrease of resolutionwith depth, due to the bias on group velocity estimates arising frominterference of the Rayleigh waves with higher modes, we carried out aprocedure of fitting, with synthetic seismograms, of selected filtered traceswith a gaussian filter, having a width at half height equal to 1 Hz and acentral frequency lying in the range [1,4] Hz. We estimated the likelihoodbetween synthetic and observed seismograms by measuring their semblance.In this way we were able to infer a more refined local velocity modelcharacterized by a high Vp and Vs vertical gradient in the sedimentarycover. Two ad hoc resolution studies, based on group velocity andamplitude data respectively, indicate that the local velocity model is a goodvelocity model also for the entire studied area. The increase in the numberof available data when using amplitude information allows us to make amore selective choice in the model parameter space (Vp and Vs of eachlayer) and to solve for the Vp/Vs ratio. The inferred Vp velocity in thehalf-space is equal to 2.8 km/s. This value is in excellent agreement withthat inferred by other authors (3 km/s) by modeling P-wave travel timevs. distance. The best-fit model furnish low Vp/Vs for the sedimentarycover so indicating a high degree of the sediment's compaction in thestudied area. The inferred shallow high-velocity gradient indicates thatthe shallow sedimentary layer in the area could trap and focus the energytraveling into it.     

Migrating sand waves

A simple mathematical model is described, which reproduces the major features of sand waves' appearance and growth and in particular predicts their migration speed. The model is based on the linear stability analysis of the flat configuration of the sea bottom subject to tidal currents. Attention is focused on the prediction of the complex growth rate that bottom perturbations undergo because of both oscillatory fluid motions and residual currents. While the real part _r of controls the amplification or decay of the amplitude of the bedforms, the imaginary part _i is related to their migration speed. Previous works on the migration of the sand waves (Németh etal. 2002) consider a forcing tide made up by the M2 constituent (oscillatory period equal to 12 h) plus the residual current Z0 and predict always a downcurrent migration of the bedforms. However, field cases exist of upcurrent-migrating sand waves (downcurrent/upcurrent-migrating sand waves mean bedforms moving in the direction of the steady residual tidal current or in the opposite direction, respectively). The inclusion of a tide constituent characterized by a period of 6 h (M4) is the main novelty of the present work, which allows for the prediction of the migration of sand waves against the residual current Z0. Indeed, the M4 tide constituent, as does also the residual current Z0, breaks the symmetry of the problem forced only by the M2 tide constituent, and induces sand-wave migration. The model proposed by Besio etal. (2003a) forms the basis for the present analysis. Previous works on the subject (Gerkema 2000; Hulscher 1996a,b; Komarova and Hulscher 2000) are thus improved by using a new solution procedure (Besio etal. 2003a) which allows for a more accurate evaluation of the growth rate for arbitrary values of the parameter r, which is the ratio between the horizontal tidal excursion and the perturbation wavelength. Responsible Editor: Jens Kappenberg     

Seismic fragility curves of as-built single-span masonry arch bridges

Masonry arch bridges are crucial elements in the railway transportation network throughout Europe. Although significant advances in seismic risk assessment of various bridge types have been made by developing fragility curves of generalized classes of structures, there are no comparable tools for masonry arch structures. In this context, this paper presents the construction of seismic fragility curves of single-span masonry bridges according to the limit analysis method. An iterative procedure is implemented to define the capacity curve of the equivalent single degree of freedom system through non-linear kinematic analysis. The process involves determination of the collapse mechanism, calculation of the limit load multiplier, and definition of the thrust line. The intrinsic variability of the seismic action is incorporated with the use of different sets of elastic spectra compatible with EC 8 Type-1 spectrum for various types of soil, with peak ground acceleration varying over the range 0.05–1.5 g. The fragility curves of the generalized classes of single-span masonry bridges are finally obtained from the effective ranges of the main geometric and material parameters affecting arch bridge capacity.     

Risk assessment strategies for the reactivation of earth flows in the Northern Apennines (Italy)

The aim of this paper is to analyze the reactivation mechanism of ancient earth flows, with a view to gleaning information that can subsequently be utilized to formulate a risk-reduction strategy. All considerations made herein are the result of direct experience and observation of actual events which have occurred over the past few decades in the Northern Apennines. Particular attention has been paid to the analysis of the evolution of landslides during actual reactivation, acknowledging a typical, recurring succession of events that precede the failure of the slope. The hazard assessment of these large landslide bodies, which are of slope scale, constitutes a thorny problem, especially in view of the inapplicability of traditional deterministic models such as limit equilibrium stability analysis. Nevertheless, a site-specific assessment of probability of reactivation of these large and ancient earth flows is fundamental to effective land-use planning.     

The calcium carbonate saturation state in cyanobacterial mats throughout Earth’s history

Through early lithification, cyanobacterial mats produced vast amounts of CaCO₃ on Precambrian carbonate platforms (before 540 Myr ago). The superposition of lithified cyanobacterial mats forms internally laminated, macroscopic structures known as stromatolites. Similar structures can be important constituents of Phanerozoic carbonate platforms (540 Myr to present). Early lithification in modern marine cyanobacterial mats is thought to be driven by a metabolically-induced increase of the CaCO₃ saturation state (Ω_CaCO3) in the mat. However, it is uncertain which microbial processes produce the Ω_CaCO3 increase and to which extent similar Ω_CaCO3 shifts were possible in Precambrian oceans whose chemistry differed from that of the modern ocean. I developed a numerical model that calculates Ω_CaCO3 in cyanobacterial mats and used it to tackle these questions. The model is first applied to simulate Ω_CaCO3 in modern calcifying cyanobacterial mats forming at Highborne Cay (Bahamas); it shows that while cyanobacterial photosynthesis increases Ω_CaCO3 considerably, sulphate reduction has a small and opposite effect on mat Ω_CaCO3 because it is coupled to H₂S oxidation with O₂ which produces acidity. Numerical experiments show that the magnitude of the Ω_CaCO3 increase is proportional to DIC in DIC-limited waters (DIC < 3-10 mM), is proportional to pH when ambient water DIC is not limiting and always proportional to the concentration of Ca²⁺ in ambient waters. With oceanic Ca²⁺ concentrations greater than a few millimolar, an appreciable increase in Ω_CaCO3 occurs in mats under a wide range of environmental conditions, including those supposed to exist in the oceans of the past 2.8 Gyr. The likely lithological expression is the formation of the microsparitic stromatolite microtexture—indicative of CaCO₃ precipitation within the mats under the control of microbial activity—which is found in carbonate rocks spanning from the Precambrian to recent. The model highlights the potential for an increase in the magnitude of the Ω_CaCO3 shift in cyanobacterial mats throughout Earth’s history produced by a decrease in salinity and temperature of the ocean, a decrease in atmospheric pCO₂ and an increase in solar irradiance. Such a trend would explain how the formation of the microsparitic stromatolite microtexture was possible as the Ω_CaCO3 of the ocean decreased from the Paleoproterozoic to the Phanerozoic.     

A numerical Round Robin on tunnels under seismic actions

Although the seismic behaviour of shallow circular tunnels in soft ground is generally safer than aboveground structures, some tunnels were recently damaged during earthquakes. In some cases, damage was associated with strong ground shaking and site amplification, which increased the stress level in the tunnel lining. Pseudo-static and simplified dynamic analyses enable to assess transient changes in internal forces during shaking. Nevertheless, experimental evidences of permanent changes in internal loads in the tunnel lining would suggest that a full dynamic analysis including plastic soil behaviour should be performed when modelling the dynamic interaction between the tunnel and the ground. While sophisticated numerical methods can be used to predict seismic internal forces on tunnel structures during earthquakes, the accuracy of their predictions should be validated against field measurements, but the latter are seldom available. A series of centrifuge tests were therefore carried out at the University of Cambridge (UK) on tunnel models in sand, in the framework of a research project funded by the Italian Civil Protection Department. A numerical Round Robin on Tunnel Tests was later promoted among some research groups to predict the observed behaviour by means of numerical modelling. In this paper, the main results of five selected numerical predictions are summarized and compared with the experimental results.     

Contrasting patterns of precipitation seasonality during the Holocene in the south‐ and north‐central Mediterranean

Pollen‐based quantitative estimates of seasonal precipitation from Lake Pergusa and lake‐level data from Lake Preola in Sicily (southern Italy) allow three successive periods to be distinguished within the Holocene: an early Holocene period before ca. 9800 cal a BP with rather dry climate conditions in winter and summer, a mid‐Holocene period between ca. 9800 and 4500 cal a BP with maximum winter and summer wetness, and a late Holocene period after 4500 cal a BP with declining winter and summer wetness. This evolution observed in the south‐central Mediterranean shows strong similarities to that recognized in the eastern Mediterranean. But, it contrasts with that reconstructed in north‐central Italy, where the mid‐Holocene appears to be characterized by a winter (summer) precipitation maximum (minimum), while the late Holocene coincided with a decrease (increase) in winter (summer) precipitation. Maximum precipitation at ca. 10 000–4500 cal a BP may have resulted from (i) increased local convection in response to a Holocene insolation maximum at 10 000 cal a BP and then (ii) the gradual weakening of the Hadley cell activity, which allowed the winter rainy westerlies to reach the Mediterranean area more frequently. After 4500 cal a BP, changes in precipitation seasonality may reflect non‐linear responses to orbitally driven insolation decrease in addition to seasonal and inter‐hemispheric changes of insolation. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical Modelling of the Anisotropic Mechanical Behaviour of Opalinus Clay at the Laboratory-Scale Using FEM/DEM

The Opalinus Clay (OPA) is an argillaceous rock formation selected to host a deep geologic repository for high-level nuclear waste in Switzerland. It has been shown that the excavation damaged zone (EDZ) in this formation is heavily affected by the anisotropic mechanical response of the material related to the presence of bedding planes. In this context, the purpose of this study is twofold: (i) to illustrate the new developments that have been introduced into the combined finite-discrete element method (FEM/DEM) to model layered materials and (ii) to demonstrate the effectiveness of this new modelling approach in simulating the short-term mechanical response of OPA at the laboratory-scale. A transversely isotropic elastic constitutive law is implemented to account for the anisotropic elastic modulus, while a procedure to incorporate a distribution of preferentially oriented defects is devised to capture the anisotropic strength. Laboratory results of indirect tensile tests and uniaxial compression tests are used to calibrate the numerical model. Emergent strength and deformation properties, together with the simulated damage mechanisms, are shown to be in strong agreement with experimental observations. Subsequently, the calibrated model is validated by investigating the effect of confinement and the influence of the loading angle with respect to the specimen anisotropy. Simulated fracture patterns are discussed in the context of the theory of brittle rock failure and analyzed with reference to the EDZ formation mechanisms observed at the Mont Terri Underground Research Laboratory.