Multiridge Euler deconvolution

Potential field interpretation can be carried out using multiscale methods. This class of methods analyses a multiscale data set, which is built by upward continuation of the original data to a number of altitudes conveniently chosen. Euler deconvolution can be cast into this multiscale environment by analysing data along ridges of potential fields, e.g., at those points along lines across scales where the field or its horizontal or vertical derivative respectively is zero. Previous work has shown that Euler equations are notably simplified along any of these ridges. Since a given anomaly may generate one or more ridges we describe in this paper how Euler deconvolution may be used to jointly invert data along all of them, so performing a multiridge Euler deconvolution. The method enjoys the stable and high‐resolution properties of multiscale methods, due to the composite upward continuation/vertical differentiation filter used. Such a physically‐based field transformation can have a positive effect on reducing both high‐wavenumber noise and interference or regional field effects. Multiridge Euler deconvolution can also be applied to the modulus of an analytic signal, gravity/magnetic gradient tensor components or Hilbert transform components. The advantages of using multiridge Euler deconvolution compared to single ridge Euler deconvolution include improved solution clustering, increased number of solutions, improvement of accuracy of the results obtainable from some types of ridges and greater ease in the selection of ridges to invert. The multiscale approach is particularly well suited to deal with non‐ideal sources. In these cases, our strategy is to find the optimal combination of upward continuation altitude range and data differentiation order, such that the field could be sensed as approximately homogeneous and then characterized by a structural index close to an integer value. This allows us to estimate depths related to the top or the centre of the structure.     

Spectral-element simulations of elastic wave propagation in exploration and geotechnical applications

We apply the spectral-element method (SEM), a high-order finite-element method (FEM) to simulate seismic wave propagation in complex media for exploration and geotechnical problems. The SEM accurately treats geometrical complexities through its flexible FEM mesh and accurately interpolates wavefields through high-order Lagrange polynomials. It has been a numerical solver used extensively in earthquake seismology. We demonstrate the applicability of SEM for selected 2D exploration and geotechnical velocity models with an open-source SEM software package SPECFEM2D. The first scenario involves a marine survey for a salt dome with the presence of major internal discontinuities, and the second example simulates seismic wave propagation for an open-pit mine with complex surface topography. Wavefield snapshots, synthetic seismograms, and peak particle velocity maps are presented to illustrate the promising use of SEM for industrial problems.     

Carbon dioxide degassing and thermal energy release in the Monte Amiata volcanic-geothermal area (Italy)

The quaternary volcanic complex of Mount Amiata is located in southern Tuscany (Italy) and represents the most recent manifestation of the Tuscan Magmatic Province. The region is characterised by a large thermal anomaly and by the presence of numerous CO₂-rich gas emissions and geothermal features, mainly located at the periphery of the volcanic complex. Two geothermal systems are located, at increasing depths, in the carbonate and metamorphic formations beneath the volcanic complex. The shallow volcanic aquifer is separated from the deep geothermal systems by a low permeability unit (Ligurian Unit). A measured CO₂ discharge through soils of 1.8 × 10⁹ mol a⁻¹ shows that large amounts of CO₂ move from the deep reservoir to the surface. A large range in δ¹³C_TDIC (−21.07 to +3.65) characterises the waters circulating in the aquifers of the region and the mass and isotopic balance of TDIC allows distinguishing a discharge of 0.3 × 10⁹ mol a⁻¹ of deeply sourced CO₂ in spring waters. The total natural CO₂ discharge (2.1 × 10⁹ mol a⁻¹) is slightly less than minimum CO₂ output estimated by an indirect method (2.8 × 10⁹ mol a⁻¹), but present-day release of 5.8 × 10⁹ mol a⁻¹ CO₂ from deep geothermal wells may have reduced natural CO₂ discharge. The heat transported by groundwater, computed considering the increase in temperature from the infiltration area to the discharge from springs, is of the same order of magnitude, or higher, than the regional conductive heat flow (>200 mW m⁻²) and reaches extremely high values (up to 2700 mW m⁻²) in the north-eastern part of the study area. Heat transfer occurs mainly by conductive heating in the volcanic aquifer and by uprising gas and vapor along fault zones and in those areas where low permeability cover is lacking. The comparison of CO₂ flux, heat flow and geological setting shows that near surface geology and hydrogeological setting play a central role in determining CO₂ degassing and heat transfer patterns.     

Geophysical investigation using resistivity and GPR methods: a case study of a lubricant oil waste disposal area in the city of Ribeirão Preto,São Paulo,Brazil

Geophysics has been shown to be effective in identifying areas contaminated by waste disposal, contributing to the greater efficiency of soundings programs and the installation of monitoring wells. In the study area, four trenches were constructed with a total volume of about 25,000 m³. They were almost totally filled with re-refined lubricating oil waste for approximately 25 years. No protection liners were used in the bottoms and laterals of the disposal trenches. The purpose of this work is to evaluate the potential of the resistivity and ground penetrating radar (GPR) methods in characterizing the contamination of this lubricant oil waste disposal area in Ribeirão Preto, SP, situated on the geological domain of the basalt spills of the Serra Geral Formation and the sandstones of the Botucatu Formation. Geophysical results were shown in 2D profiles. The geophysical methods used enabled the identification of geophysical anomalies, which characterized the contamination produced by the trenches filled with lubricant oil waste. Conductive anomalies (smaller than 185 Ωm) immediately below the trenches suggest the action of bacteria in the hydrocarbons, as has been observed in several sites contaminated by hydrocarbons in previously reported cases in the literature. It was also possible to define the geometry of the trenches, as evidenced by the GPR method. Direct sampling (chemical analysis of the soil and the water in the monitoring well) confirmed the contamination. In the soil analysis, low concentrations of several polycyclic aromatic hydrocarbons (PAHs) were found, mainly naphthalene and phenanthrene. In the water samples, an analysis verified contamination of the groundwater by lead (Pb). The geophysical methods used in the investigation provided an excellent tool for environmental characterization in this study of a lubricant oil waste disposal area, and could be applied in the study of similar areas.     

U?Pb and 40Ar?39Ar geochronology of Palaeozoic units in the northern Apennines: determining protolith age and alpine evolution using the Calamita Schist and Ortano Porphyroid

In the northern Apennines, the Palaeozoic basement involved in the Late Oligocene–Middle Miocene nappe stack contains metamorphic units for which hypothetical ages have been assigned on the basis of lithological correlations with the Palaeozoic formations of the Variscan chain in Sardinia. This uncertainty concerning the age poses limitations to reconstructing the Palaeozoic stratigraphy, defining the Alpine and pre‐Alpine histories and correlations with other domains of the Variscan chain. We present the U Pb age of detrital zircon and the ⁴⁰Ar ³⁹Ar age of metamorphic muscovite for the Calamita Schist and Ortano Porphyroid, two metamorphic units of undetermined Palaeozoic age cropping out in the eastern Elba Island. The radioisotopic data allows us to: (i) define the Early Carboniferous and Middle Ordovician ages for the Calamita Schist and Ortano Porphyroid, respectively, as well as their derivation (flysch deposit and magmatic rocks); (ii) pose some constraints concerning their alpine tectonic and metamorphic histories. These new data generate a more precise reconstruction of the Palaeozoic sequence in the northern Apennines, and they document that the Palaeozoic basement involved in the alpine deformation underwent internal stacking with an inversion of the original sequence. Copyright © 2010 John Wiley & Sons, Ltd.     

Comparing models of debris-flow susceptibility in the alpine environment

Debris-flows are widespread in Val di Fassa (Trento Province, Eastern Italian Alps) where they constitute one of the most dangerous gravity-induced surface processes. From a large set of environmental characteristics and a detailed inventory of debris flows, we developed five models to predict location of debris-flow source areas. The models differ in approach (statistical vs. physically-based) and type of terrain unit of reference (slope unit vs. grid cell). In the statistical models, a mix of several environmental factors classified areas with different debris-flow susceptibility; however, the factors that exert a strong discriminant power reduce to conditions of high slope-gradient, pasture or no vegetation cover, availability of detrital material, and active erosional processes. Since slope and land use are also used in the physically-based approach, all model results are largely controlled by the same leading variables.Overlaying susceptibility maps produced by the different methods (statistical vs. physically-based) for the same terrain unit of reference (grid cell) reveals a large difference, nearly 25% spatial mismatch. The spatial discrepancy exceeds 30% for susceptibility maps generated by the same method (discriminant analysis) but different terrain units (slope unit vs. grid cell). The size of the terrain unit also led to different susceptibility maps (almost 20% spatial mismatch). Maps based on different statistical tools (discriminant analysis vs. logistic regression) differed least (less than 10%). Hence, method and terrain unit proved to be equally important in mapping susceptibility.Model performance was evaluated from the percentages of terrain units that each model correctly classifies, the number of debris-flow falling within the area classified as unstable by each model, and through the metric of ROC curves. Although all techniques implemented yielded results essentially comparable; the discriminant model based on the partition of the study area into small slope units may constitute the most suitable approach to regional debris-flow assessment in the Alpine environment.     

Study of Shallow Low-Enthalpy Geothermal Resources Using Integrated Geophysical Methods

The paper is focused on low enthalpy geothermal exploration performed in south Italy and provides an integrated presentation of geological, hydrogeological, and geophysical surveys carried out in the area of municipality of Lecce. Geological and hydrogeological models were performed using the stratigraphical data from 51 wells. A ground-water flow (direction and velocity) model was obtained. Using the same wells data, the ground-water annual temperature was modeled. Furthermore, the ground surface temperature records from ten meteorological stations were studied. This allowed us to obtain a model related to the variations of the temperature at different depths in the subsoil. Integrated geophysical surveys were carried out in order to explore the low-enthalpy geothermal fluids and to evaluate the results of the model. Electrical resistivity tomography (ERT) and self-potential (SP) methods were used. The results obtained upon integrating the geophysical data with the models show a low-enthalpy geothermal resource constituted by a shallow ground-water system.     

Improved 2-D attenuation analysis for Northern Italy using a merged dataset from selected regional seismic networks

A merged, high-quality waveform dataset from different seismic networks has been used to improve our understanding of lateral seismic attenuation for Northern Italy. In a previous study on the same region, Morasca et al. (Bull Seismol Soc Am 98:1936–1946, 2008) were able to resolve only a small area due to limited data coverage. For this reason, the interpretation of the attenuation anomalies was difficult given the complexity of the region and the poor resolution of the available data. In order to better understand the lateral changes in the crustal structure and thickness of this region, we selected 770 earthquakes recorded by 54 stations for a total of almost 16,000 waveforms derived from seismic networks operating totally or partially in Northern Italy. Direct S-wave and coda attenuation images were obtained using an amplitude ratio technique that eliminates source terms from the formulation. Both direct and early-coda amplitudes are used as input for the inversions, and the results are compared. Results were obtained for various frequency bands ranging between 0.3 and 25.0 Hz and in all cases show significant improvement with respect to the previous study since the resolved area has been extended and more crossing paths have been used to image smaller scale anomalies. Quality-factor estimates are consistent with the regional tectonic structure exhibiting a general trend of low attenuation under the Po Plain basin and higher values for the Western Alps and Northern Apennines. The interpretation of the results for the Eastern Alps is not simple, possibly because our resolution for this area is still not adequate to resolve small-scale structures.     

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.     

Numerical reconstruction of trajectory of small-size surface drifter in the Mediterranean sea

In this study, we addressed the effects of wind-induced drift on Lagrangian trajectories of surface sea objects using high-resolution ocean forecast and atmospheric data. Application of stochastic Leeway model for prediction of trajectories drift was considered for the numerical reconstruction of the Elba accident that occurred during the period 21.06.2009–22.06.2009: a person on an inflatable raft was lost in the vicinity of the Elba Island coast; from the initial position, the person on a raft drifted southwards in the open sea and later he was found on a partially deflated raft during rescue operation. For geophysical forcing, we used high-resolution currents from the Mediterranean Forecasting System and atmospheric wind from the European Centre for Medium-Range Weather Forecasts. To investigate the effect of wind on trajectory behavior, numerical simulations were performed using different categories of drifter-like particles similar to a person on an inflatable raft. An algorithm of spatial clustering was used to differentiate the most probable search areas with a high density of particles. Our results showed that the simulation scenarios using particles with characteristics of draft-limited sea drifters provided better prediction of an observed trajectory in terms of the probability density of particles.