Search

Following the linearized attenuation tomography from our previous study (Pejić et al., 2017, https://doi.org/10.1002/2016JB013692), we perform hierarchical transdimensional Bayesian tomography of the upper ≈400 km of the inner core, using 398 globally distributed t∗ estimates. The results are in good agreement with the ones obtained through linearized tomography: they show more complex attenuation pattern than the purely hemispherical one, and the noise estimated from the hierarchical inversion is in good agreement with estimates obtained from the Discrepancy Principle in the previous study. The attenuation pattern we observe gives more weight to the geodynamical models that couple the thermal anomalies of the lowermost mantle to the inner core boundary. Key Points: Transdimensional sampling of Voronoi cells on a spherical surface applied to a problem in global tomographyNo hemispherical pattern is observed in attenuation structure of the upper inner coreThe study provides a strong seismological hint for a strongly attenuating inner core [ABSTRACT FROM AUTHOR]

Most linear inverse problems require regularization to ensure that robust and meaningful solutions can be found. Typically, Tikhonov-style regularization is used, whereby a preference is expressed for models that are somehow 'small' and/or 'smooth'. The strength of such preferences is expressed through one or more damping parameters, which control the character of the solution, and which must be set by the user. However, identifying appropriate values is often regarded as a matter of art, guided by various heuristics. As a result, such choices have often been the source of controversy and concern. By treating these as hyperparameters within a hierarchical Bayesian framework, we are able to obtain solutions that encompass the range of permissible regularization parameters. Furthermore, we show that these solutions are often well-approximated by those obtained via standard analysis using certain regularization choices which are—in a certain sense—optimal. We obtain algorithms for determining these optimal values in various cases of common interest, and show that they generate solutions with a number of attractive properties. A reference implementation of these algorithms, written in Python, accompanies this paper. [ABSTRACT FROM AUTHOR]

In geophysical inversion, inferences of Earth's properties from sparse data involve a trade-off between model complexity and the spatial resolving power. A recent Markov chain Monte Carlo (McMC) technique formalized by Green, the so-called trans-dimensional samplers, allows us to sample between these trade-offs and to parsimoniously arbitrate between the varying complexity of candidate models. Here we present a novel framework using transdimensional sampling over tree structures. This new class of McMC sampler can be applied to 1-D, 2-D and 3-D Cartesian and spherical geometries. In addition, the basis functions used by the algorithm are flexible and can include more advanced parametrizations such as wavelets, both in Cartesian and Spherical geometries, to permit Bayesian multiscale analysis. This new framework offers greater flexibility, performance and efficiency for geophysical imaging problems than previous sampling algorithms. Thereby increasing the range of applications and in particular allowing extension to trans-dimensional imaging in 3-D. Examples are presented of its application to 2-D seismic and 3-D teleseismic tomography including estimation of uncertainty. [ABSTRACT FROM AUTHOR]

We conduct an experiment to investigate whether linearity in the observed velocity gradient in the volume of the inner core sampled by the PKP ray paths beneath Central America is a robust approximation. Instead of solving an optimization problem, we approach it within the Bayesian inference. This is an ensemble approach, where model specification is relaxed so that instead of only one solution, groups of reasonable models are acceptable. Furthermore, in transdimensional Bayesian inference used here, the number of basis functions needed to model observations is by itself an unknown. Our modeling reveals that in the ensemble of models, the most likely are those containing only 2 nodes (linear trend). Thus our result justifies the assumption used for the determination of inner core rotation with respect to the rest of the mantle that the observed gradient is constant in its nature (linear). Recent observations in seismology suggest that it is likely that the spatial variability in elastic parameters is a widespread phenomenon in the inner core. Future array observations will further constrain spatial extent and magnitude of velocity changes and show whether there is a significant difference between these observations in the two quasi-hemispheres of the inner core. [ABSTRACT FROM AUTHOR]

This paper applies parallel tempering within a Bayesian formulation for strongly nonlinear geoacoustic inverse problems. Bayesian geoacoustic inversion consists of sampling the posterior probability density (PPD) of seabed parameters to estimate integral properties, such as marginal probability distributions, based on ocean acoustic data and prior information. This sampling is usually carried out using the Markov-chain Monte Carlo method of Metropolis-Hastings sampling. However, standard sampling methods can be very inefficient for strongly nonlinear problems involving multi-modal PPDs with the potential to miss important regions of the parameter space and to significantly underestimate parameter uncertainties. Parallel tempering achieves efficient/effective sampling of challenging parameter spaces with the ability to transition freely between multiple PPD modes by running parallel Markov chains at a series of increasing sampling temperatures with probabilistic interchanges between chains. The approach is illustrated for inversion of (simulated) acoustic reverberation data for which the PPD is highly multi-modal. While Metropolis-Hastings sampling gives poor results even with very large sample sizes, parallel tempering provides efficient, convergent sampling of the PPD. Methods to enhance the efficiency of parallel tempering are also considered. [ABSTRACT FROM AUTHOR]

The article discusses various papers published in this issue including one by Biswas and Sharma on self-potential anomaly over a two-dimensional inclined structure using very fast simulated-annealing global optimization, one by Coles et al. on post-inverison model uncertainty and one by Young et al. on Bayesian tomographic inversion scheme.