Regularized tomographic inversion with geological constraints

Reflection tomography is the industry standard tool for velocity model building, but it is also an ill‐posed inverse problem as its solution is not unique. The usual way to obtain an acceptable result is to regularize tomography by feeding the inversion with some a priori information. The simplest regularization forces the solution to be smooth, implicitly assuming that seismic velocity exhibits some degree of spatial correlation. However, velocity is a rock property; thus, the geometry and structure of rock formations should drive correlation in velocity depth models. This observation calls for constraints driven by geological models. In this work, we present a set of structural constraints that feed reflection tomography with geometrical information. These constraints impose the desired characteristics (flatness, shape, position, etc.) on imaged reflectors but act on the velocity update. Failure to respect the constraints indicates either velocity inaccuracies or wrong assumptions concerning the constraints. Reflection tomography with structural constraints is a flexible framework that can be specialized in order to achieve different goals: among others, to flatten the base of salt bodies or detachment surfaces, to recover the horizontalness of oil–water contacts, or to impose the co‐location of the same imaged horizon between PP and PS images. The straightforward application of structural constraints is that of regularizing tomography through geological information, particularly at the latest stages of the depth imaging workflow, when the depth migration structural setting reached a consistent geological interpretation. Structural constraints are also useful in minimizing the well‐to‐seismic mis‐ties. Moreover, they can be used as a tool to check the consistency of interpreters' hypothesis with seismic data. Indeed, inversion with structural constraints will preserve image focusing only if the interpreters' insights are consistent with the data. Results from synthetic and real data demonstrate the effectiveness of reflection tomography with structural constraints.     

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.     

Rate of soil erosion processes on vineyards in central Piedmont (NW Italy)

The ‘Tertiary Basin of Piedmont’ is a hilly region, mainly composed of marine sediments, such as marls, silts and sands. The slopes, largely devoted to grape production, are usually kept bare of vegetation and are thus prone to soil erosion processes. For two years we have measured soil loss in relation to rainfall on experimental plots located in vineyards. In all the plots considered erosion started with low rainfall intensities (0·07 mm/min), and above 0·4 mm/min the amount of soil loss dramatically increased (to over 1800 g/m in one event). Most of the erosion occurs during summer rainstorms, but the behaviour of the soil under erosive rainfalls is very different from one site to another, depending only in part upon the various rainfall rates, soil and geometrical characteristics of the plots. It can be seen that soil loss varies from nearly negligible values (20 g/m/yr) to unacceptably high levels (to over 4·7 kg/m/yr), according to the different types of land cultivation. Deep ploughing and heavy herbicide treatment give rise to accelerated erosion processes leading to soil losses much higher than in other vineyard plots in Europe.     

Vertical variations of the chlorophyll maximum during a red tide in a shallow lagoon

The results of eight sets of repeated observations on the vertical variations of the chlorophyll maximum layer in a shallow lagoon during a red tide show that these were more frequently hydrologically induced, rather than due to active vertical migrations of the red tide-forming organism. These results are discussed and compared to those existing in the literature, with special regard to the role of light and nitrogen in conditioning vertical migrations in red tide-forming dinoflagellates.     

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.     

Evaluating Earthquake-Triggered Landslide Hazard at the Basin Scale Through Gis in the Upper Sele River Valley

To evaluate techniques for assessing earthquake-triggeredlandslide hazard in the Southern Apennines (Italy), a GIS-based analysis was used to modelseismically induced slope deformations. Geological, geotechnical, geomorphological and seismologicaldata were integrated into a standard earthquake slope stability model. The model assessed thelandslide potential that existed during the 1980 Irpinian earthquake in the Upper Sele river Valley.The standard Newmark displacement analysis, widely used for predicting the location of shallowunstable slopes, does not take into account errors and/or uncertainties in the input parameters.Therefore, a probabilistic Newmark displacement analysis technique has been used. Probabilistictechniques allow, e.g., an estimation of the probability that a slope will exceed a certain criticalvalue of Newmark displacement. In our probabilistic method, a Monte-Carlo based simulation modelis used in conjunction with a GIS. The random variability of geotechnical data is modelled by probabilitydensity functions (pdfs), while for the seismic input three different regression laws wereconsidered. Input probability distributions are sampled and the resulting values input into empiricalrelations for estimating Newmark displacement. The outcome is a map in which to each siteis related a spatial probability distribution for the expected displacement in response to seismic loading.Results of the experiments show a high grade of uncertainty in the application of the Newmarkanalysis both for the deterministic and probabilistic approach in a complex geological setting suchas the high Sele valley, quite common in the Southern Apennines. They show a strong dependence onthe reliability of the spatial data used in input, so that, when the model is used at basin scale,results are strongly influenced by local environmental condition (e.g., topography, lithology, groundwatercondition) and decrease the model performance.     

Considering large earthquake clustering in seismic risk analysis

Earthquakes cluster in space and time resulting in nonlinear damage effects. We compute earthquake interactions using the Coulomb stress transfer theory and dynamic vulnerability from the concept of ductility capacity reduction. We combine both processes in the generic multi-risk framework where risk scenarios are simulated using a variant of the Markov chain Monte Carlo method. We apply the proposed approach to the thrust fault system of northern Italy, considering earthquakes with characteristic magnitudes in the range ~[6, 6.5], different levels of tectonic loading \(\dot{\tau }\) = {10⁻⁴, 10⁻³, 10⁻²} bar/year and a generic stock of fictitious low-rise buildings with different ductility capacities μ _Δ = {2, 4, 6}. We describe the process’ stochasticity by non-stationary Poisson earthquake probabilities and by binomial damage state probabilities. We find that earthquake clustering yields a tail fattening of the seismic risk curve, the effect of which is amplified by damage-dependent fragility due to clustering. The impact of clustering alone is in average more important than dynamic vulnerability, the spatial extent of the former phenomenon being greater than of the latter one.

     

Setting-up local “mapping geoids” with the aid of GPS/LEV traverses. Application to the geoids of Sardinia and Calabria

The problem of fitting a pre-existing geoid on local GPS/LEV points is first taken into account, a possible way for its solution is proposed, based on techniques similar to the aerophotogrammetric absolute orientation.This approach is then used to compute themapping geoids of Sardinia and Calabria i.e. geoids with decimetric accuracy, needed for the heigthing of large scale maps. The pre-existing geoids EAGGI by Brennecke et al. (1983), and EGG2 by Denker et al (1993) are utilized, fitted on points of GPS/LEV traverses.The complete input-output data, the contour maps of both geoids, and a substantial checking on some Calabrian GPS/TRIG points of the entreprise IGM.95 are finally reported.