Merging Active and Passive Data Sets in Travel-Time Tomography: The Case Study of Campi Flegrei Caldera (Southern Italy)

We propose a strategy for merging both active and passive data sets in linearized tomographic inversion. We illustrate this in the reconstruction of 3D images of a complex volcanic structure, the Campi Flegrei caldera, located in the vicinity of the city of Naples, southern Italy. The caldera is occasionally the site of significant unrests characterized by large ground uplifts and seismicity. The P and S velocity models of the caldera structure are obtained by a tomographic inversion based on travel times recorded during two distinct experiments. The first data set is composed of 606 earthquakes recorded in 1984 and the second set is composed of recordings for 1528 shots produced during the SERAPIS experiment in 2001. The tomographic inversion is performed using an improved method based on an accurate finite-difference traveltime computation and a simultaneous inversion of both velocity models and earthquake locations. In order to determine the adequate inversion parameters and relative data weighting factors, we perform massive synthetic simulations allowing one to merge the two types of data optimally. The proper merging provides high resolution velocity models, which allow one to reliably retrieve velocity anomalies over a large part of the tomography area. The obtained images confirm the presence of a high P velocity ring in the southern part of the bay of Pozzuoli and extends its trace inland as compared to previous results. This annular anomaly represents the buried trace of the rim of the Campi Flegrei caldera. Its shape at 1.5 km depth is in good agreement with the location of hydrothermalized lava inferred by gravimetric data modelling. The Vp/Vs model confirms the presence of two characteristic features. At about 1 km depth a very high Vp/Vs anomaly is observed below the town of Pozzuoli and is interpreted as due to the presence of rocks that contain fluids in the liquid phase. A low Vp/Vs body extending at about 3–4 km depth below a large part of the caldera is interpreted as the top of formations that are enriched in gas under supercritical conditions

Converted phase identification and retrieval of Vp/Vs ratios from move-out reflection analysis: application to the Campi Flegrei caldera

Here, we propose a method for the determination of Vp/Vs ratios in a horizontally layered propagation media using maximization of a coherency function along theoretical travel-times of PS reflected phases. The theoretical travel-times are computed using the information about the propagation media that is extracted by velocity analysis or by topographic analysis performed on the first arrivals. The method is also a valid tool for the identification of the PS phases associated with a fixed seismic reflector, and it is particularly suitable for data that is stored in common mid-point and common conversion point bin- ning; for this kind of data the hypothesis of horizontally and layered media can usually be verified. We applied the method to both simulated and real datasets. The use of the real data that was acquired in the Campi Flegrei caldera (southern Italy) allowed us to estimate a relatively high Vp/Vs ratio (3.5 ± 0.6) for a very shallow layer (maximum depth, 600 m). This hypothesis has been tested by theoretical rock physical modeling of the Vp/Vs ratios as a function of porosity, suggesting that the shallow layer appears to be formed of unconsolidated, water-saturated, volcanic and marine sediments that filled Pozzuoli Bay during the post-caldera activity

On the use of simulated annealing method and cross-validation theory for deconvolution of seismograms Bull

Abstract In order to retrieve the apparent source time function (ASTF) seen at a given station, one must take into account propagation effects, site, and instrumental influences. Removing these effects can be performed by a deconvolution of the mainshock seismogram by a seismogram of a smaller event recorded at the same station. This smaller event must occur nearby the mainshock, and the associated seismogram is considered as an empirical Green's function. We propose a deconvolution based on the simulated annealing method, and we compare it with the often-used spectral division technique. We show on both synthetic and real signals that the simulated annealing deconvolution (SAD) provides stable and positive ASTF, whereas results from spectral division are very sensitive to an ad hoc parameter called water level. Finally, the application of cross-validation analysis between the three components of the seismogram in addition to the SAD allows us to estimate errors on the ASTF

Challenges in shallow target reconstruction by 3D elastic full-waveform inversion - Which initial model?

Elastic full-waveform inversion (FWI) is a powerful tool for high-resolution subsurface multiparameter characterization. However, 3D FWI applied to land data for near-surface applications is particularly challenging because the seismograms are dominated by highly energetic, dispersive, and complex-scattered surface waves (SWs). In these conditions, a successful deterministic FWI scheme requires an accurate initial model. Our study, primarily focused on field data analysis for 3D applications, aims at enhancing the resolution in the imaging of complex shallow targets, by integrating devoted SW analysis techniques with a 3D spectral-element-based elastic FWI. From dispersion curves, extracted from seismic data recorded over a sharp-interface shallow target, we build different initial S-wave (VS) and P-wave (VP) velocity models (laterally homogeneous and laterally variable), using a specific data transform. Starting from these models, we carry out 3D FWI tests on synthetic and field data, using a relatively straightforward inversion scheme. The field data processing before FWI consists of band-pass filtering and muting of noisy traces. During FWI, a weighting function is applied to the far-offset traces. We test 2D and 3D acquisition layouts, with different positions of the sources and variable offsets. The 3D FWI workflow enriches the overall content of the initial models, allowing a reliable reconstruction of the shallow target, especially when using laterally variable initial models. Moreover, a 3D acquisition layout guarantees a better reconstruction of the target's shape and lateral extension. In addition, the integration of model-oriented (preliminary monoparametric FWI) and data-oriented (time windowing) strategies into the main optimization scheme has produced further improvement of the FWI results

A field expansions method for scattering by periodic multilayered media

The interaction of acoustic and electromagnetic waves with periodic structures plays an important role in a wide range of problems of scientific and technological interest. This contribution focuses upon the robust and high-order numerical simulation of a model for the interaction of pressure waves generated within the earth incident upon layers of sediment near the surface. Herein described is a boundary perturbation method for the numerical simulation of scattering returns from irregularly shaped periodic layered media. The method requires only the discretization of the layer interfaces (so that the number of unknowns is an order of magnitude smaller than finite difference and finite element simulations), while it avoids not only the need for specialized quadrature rules but also the dense linear systems characteristic of boundary integral/element methods. The approach is a generalization to multiple layers of Bruno and Reitich’s “Method of Field Expansions” for dielectric structures with two layers. By simply considering the entire structure simultaneously, rather than solving in individual layers separately, the full field can be recovered in time proportional to the number of interfaces. As with the original field expansions method, this approach is extremely efficient and spectrally accurate

LSST: from Science Drivers to Reference Design and Anticipated Data Products

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.Comment: 57 pages, 32 color figures, version with high-resolution figures available from https://www.lsst.org/overvie

Study of large hemispherical photomultiplier tubes for the ANTARES neutrino telescope

The ANTARES neutrino telescope, to be immersed depth in the Mediterranean Sea, will consist of a 3 dimensional matrix of 900 large area photomultiplier tubes housed in pressure resistant glass spheres. The selection of the optimal photomultiplier was a critical step for the project and required an intensive phase of tests and developments carried out in close collaboration with the main manufacturers worldwide. This paper provides an overview of the tests performed by the collaboration and describes in detail the features of the PMT chosen for ANTARES

Microseismic Full Waveform Modeling in Anisotropic Media with Moment Tensor Implementation

Seismic anisotropy which is common in shale and fractured rocks will cause travel-time and amplitude discrepancy in different propagation directions. For microseismic monitoring which is often implemented in shale or fractured rocks, seismic anisotropy needs to be carefully accounted for in source location and mechanism determination. We have developed an efficient finite-difference full waveform modeling tool with an arbitrary moment tensor source. The modeling tool is suitable for simulating wave propagation in anisotropic media for microseismic monitoring. As both dislocation and non-double-couple source are often observed in microseismic monitoring, an arbitrary moment tensor source is implemented in our forward modeling tool. The increments of shear stress are equally distributed on the staggered grid to implement an accurate and symmetric moment tensor source. Our modeling tool provides an efficient way to obtain the Green’s function in anisotropic media, which is the key of anisotropic moment tensor inversion and source mechanism characterization in microseismic monitoring. In our research, wavefields in anisotropic media have been carefully simulated and analyzed in both surface array and downhole array. The variation characteristics of travel-time and amplitude of direct P- and S-wave in vertical transverse isotropic media and horizontal transverse isotropic media are distinct, thus providing a feasible way to distinguish and identify the anisotropic type of the subsurface. Analyzing the travel-times and amplitudes of the microseismic data is a feasible way to estimate the orientation and density of the induced cracks in hydraulic fracturing. Our anisotropic modeling tool can be used to generate and analyze microseismic full wavefield with full moment tensor source in anisotropic media, which can help promote the anisotropic interpretation and inversion of field data