Description: Amplitudes and shapes of seismic patterns derived from tomographic images often are strongly biased with respect to real structures in the earth. In particular, tomography usually provides continuous velocity distributions, whereas major velocity changes in the earth often occur on first-order interfaces. We propose an approach that constructs a realistic structure of the earth that combines forward modeling and tomographic inversion (FM&TI). Using available a priori information, we first construct a synthetic model with realistic patterns. Then we compute synthetic times and invert them using the same tomographic code and the same parameters as in the case of observed data processing. We compare the reconstruction result with the tomographicimage of observed data inversion. If a discrepancy is observed, we correct the synthetic model and repeat the FM&TI process. After several trials, we obtain similar results of synthetic and observed data inversion. In this case, the derived synthetic model adequately represents the real structure of the earth. In a working scheme of this approach, we three authors used two different synthetic models with a realistic setup. One of us created models, but the other two performed the reconstruction with no knowledge of the models. We discovered that the synthetic models derived by FM&TI were closer to the true model than the tomographic inversion result. Our reconstruction results from modeling marine data acquired in the Musicians Seamount Province in the Pacific Ocean indicate the capacity and limitations of FM&TI.

Description: General solutions of inverse problems often can be obtained by introducing probability distributions to sample the model space. We have developed a simple approach to define an a priori space in a tomographic study and retrieve the velocity/depth posterior distribution by a Monte Carlo method. Utilizing a fitting routine designed for very low statistics to set up and analyze the obtained tomography results, we can statistically separate the velocity/depth model space derived from inverting seismic refraction data. A profile acquired in the Lesser Antilles subduction zone reveals the effectiveness of this approach. Resolution analysis of the structural heterogeneity includes a divergence analysis that can dissect long, wide-angle profiles for deep crust and upper mantle studies. The complete information of any parameterized physical system is contained in the a posteriori distribution. Methods for analyzing and displaying key properties of the a posteriori distributions of highly nonlinear inverse problems are therefore essential in the scope of any interpretation. It is possible to map velocity variations in their extent and structure by measuring the total as well as relative divergence of the velocity structure in the a posteriori space. We have applied the divergence analysis to a part of the transect where a backstop structure has been identified, and the method resolves shallow features and returns information concerning the confidence level of results. Assuming a relationship between forearc and backstop, we can obtain a structural image in accordance with previous interpretations. ©2010 Society of Exploration Geophysicists

Description: To constrain seismic anisotropy under and around the Alps in Europe, we study SKS shear wave splitting from the region densely covered by the AlpArray seismic network. We apply a technique based on measuring the splitting intensity, constraining well both the fast orientation and the splitting delay. Four years of teleseismic earthquake data were processed, from 723 temporary and permanent broad-band stations of the AlpArray deployment including ocean-bottom seismometers, providing a spatial coverage that is unprecedented. The technique is applied automatically (without human intervention), and it thus provides a reproducible image of anisotropic structure in and around the Alpine region. As in earlier studies, we observe a coherent rotation of fast axes in the western part of the Alpine chain, and a region of homogeneous fast orientation in the Central Alps. The spatial variation of splitting delay times is particularly interesting though. On one hand, there is a clear positive correlation with Alpine topography, suggesting that part of the seismic anisotropy (deformation) is caused by the Alpine orogeny. On the other hand, anisotropic strength around the mountain chain shows a distinct contrast between the Western and Eastern Alps. This difference is best explained by the more active mantle flow around the Western Alps. The new observational constraints, especially the splitting delay, provide new information on Alpine geodynamics.

Description: Probing seismic anisotropy of the lithosphere provides valuable clues on the fabric of rocks. We present a 3-D probabilistic model of shear wave velocity and radial anisotropy of the crust and uppermost mantle of Europe, focusing on the mountain belts of the Alps and Apennines. The model is built from Love and Rayleigh dispersion curves in the period range 5–149 s. Data are extracted from seismic ambient noise recorded at 1521 broad-band stations, including the AlpArray network. The dispersion curves are first combined in a linearized least squares inversion to obtain 2-D maps of group velocity at each period. Love and Rayleigh maps are then jointly inverted at depth for shear wave velocity and radial anisotropy using a Bayesian Monte Carlo scheme that accounts for the trade-off between radial anisotropy and horizontal layering. The isotropic part of our model is consistent with previous studies. However, our anisotropy maps differ from previous large scale studies that suggested the presence of significant radial anisotropy everywhere in the European crust and shallow upper mantle. We observe instead that radial anisotropy is mostly localized beneath the Apennines while most of the remaining European crust and shallow upper mantle is isotropic. We attribute this difference to trade-offs between radial anisotropy and thin (hectometric) layering in previous studies based on least-squares inversions and long period data (〉30 s). In contrast, our approach involves a massive data set of short period measurements and a Bayesian inversion that accounts for thin layering. The positive radial anisotropy (VSH 〉 VSV) observed in the lower crust of the Apennines cannot result from thin layering. We rather attribute it to ductile horizontal flow in response to the recent and present-day extension in the region.