Trans-dimensional inversion of modal dispersion data on the New England Mud Patch

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Bonnel, J., Dosso, S. E., Eleftherakis, D., & Chapman, N. R. Trans-dimensional inversion of modal dispersion data on the New England Mud Patch. IEEE Journal of Oceanic Engineering, 45(1), (2020): 116-130, doi:10.1109/JOE.2019.2896389.This paper presents single receiver geoacoustic inversion of two independent data sets recorded during the 2017 seabed characterization experiment on the New England Mud Patch. In the experimental area, the water depth is around 70 m, and the seabed is characterized by an upper layer of fine grained sediments with clay (i.e., mud). The first data set considered in this paper is a combustive sound source signal, and the second is a chirp emitted by a J15 source. These two data sets provide differing information on the geoacoustic properties of the seabed, as a result of their differing frequency content, and the dispersion properties of the environment. For both data sets, source/receiver range is about 7 km, and modal time-frequency dispersion curves are estimated using warping. Estimated dispersion curves are then used as input data for a Bayesian trans-dimensional inversion algorithm. Subbottom layering and geoacoustic parameters (sound speed and density) are thus inferred from the data. This paper highlights important properties of the mud, consistent with independent in situ measurements. It also demonstrates how information content differs for two data sets collected on reciprocal tracks, but with different acoustic sources and modal content.10.13039/100000006-Office of Naval Research 10.13039/100007297-Office of Naval Research Globa

Basis Pursuit Receiver Function

Receiver functions (RFs) are derived by deconvolution of the horizontal (radial or transverse) component of ground motion from the vertical component, which segregates the PS phases. Many methods have been proposed to employ deconvolution in frequency as well as in time domain. These methods vary in their approaches to impose regularization that addresses the stability problem. Here, we present application of a new time-domain deconvolution technique called basis pursuit deconvolution (BPD) that has recently been applied to seismic exploration data. Unlike conventional deconvolution methods, the BPD uses an L1 norm constraint on model reflectivity to impose sparsity. In addition, it uses an overcomplete wedge dictionary based on a dipole reflectivity series to define model constraints, which can achieve higher resolution than that obtained by the traditional methods. We demonstrate successful application of BPD based RF estimation from synthetic data for a crustal model with a near-surface thin layer of thickness 5, 7, 10, and 15 km. The BPD can resolve these thin layers better with much improved signal-to-noise ratio than the conventional methods. Finally, we demonstrate application of the BPD receiver function (BPRF) method to a field dataset from Kutch, India, where near-surface sedimentary layers are known to be present. The BPRFs are able to resolve reflections from these layers very well.Jackson Chair funds at the Jackson School of Geosciences, University of Texas, AustinCouncil of Scientific and Industrial Research twelfth five year plan project at the Council of Scientific and Industrial Research National Geophysical Research Institute (CSIR-NGRI), HyderabadInstitute for Geophysic

Highly efficient Bayesian joint inversion for receiver-based data and its application to lithospheric structure beneath the southern Korean Peninsula

With the deployment of extensive seismic arrays, systematic and efficient parameter and uncertainty estimation is of increasing importance and can provide reliable, regional models for crustal and upper-mantle structure.We present an efficient Bayesian method for the joint inversion of surface-wave dispersion and receiver-function data that combines trans-dimensional (trans-D) model selection in an optimization phase with subsequent rigorous parameter uncertainty estimation. Parameter and uncertainty estimation depend strongly on the chosen parametrization such that meaningful regional comparison requires quantitative model selection that can be carried out efficiently at several sites. While significant progress has been made for model selection (e.g. trans-D inference) at individual sites, the lack of efficiency can prohibit application to large data volumes or cause questionable results due to lack of convergence. Studies that address large numbers of data sets have mostly ignored model selection in favour of more efficient/simple estimation techniques (i.e. focusing on uncertainty estimation but employing ad-hoc model choices). Our approach consists of a two-phase inversion that combines trans-D optimization to select the most probable parametrization with subsequent Bayesian sampling for uncertainty estimation given that parametrization. The trans-D optimization is implemented here by replacing the likelihood function with the Bayesian information criterion (BIC). The BIC provides constraints on model complexity that facilitate the search for an optimal parametrization. Parallel tempering (PT) is applied as an optimization algorithm. After optimization, the optimal model choice is identified by the minimum BIC value from all PT chains. Uncertainty estimation is then carried out in fixed dimension. Data errors are estimated as part of the inference problem by a combination of empirical and hierarchical estimation. Data covariance matrices are estimated from data residuals (the difference between prediction and observation) and periodically updated. In addition, a scaling factor for the covariance matrix magnitude is estimated as part of the inversion. The inversion is applied to both simulated and observed data that consist of phase- and group-velocity dispersion curves (Rayleigh wave), and receiver functions. The simulation results show that model complexity and important features are well estimated by the fixed dimensional posterior probability density. Observed data for stations in different tectonic regions of the southern Korean Peninsula are considered. The results are consistent with published results, but important features are better constrained than in previous regularized inversions and are more consistent across the stations. For example, resolution of crustal and Moho interfaces, and absolute values and gradients of velocities in lower crust and upper mantle are better constrained

Nonlinear time-warping made simple: a step-by-step tutorial on underwater acoustic modal separation with a single hydrophone

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Bonnel, J., Thode, A., Wright, D., & Chapman, R. Nonlinear time-warping made simple: a step-by-step tutorial on underwater acoustic modal separation with a single hydrophone. The Journal of the Acoustical Society of America, 147(3), (2020): 1897, doi:10.1121/10.0000937.Classical ocean acoustic experiments involve the use of synchronized arrays of sensors. However, the need to cover large areas and/or the use of small robotic platforms has evoked interest in single-hydrophone processing methods for localizing a source or characterizing the propagation environment. One such processing method is “warping,” a non-linear, physics-based signal processing tool dedicated to decomposing multipath features of low-frequency transient signals (frequency f  1 km). Since its introduction to the underwater acoustics community in 2010, warping has been adopted in the ocean acoustics literature, mostly as a pre-processing method for single receiver geoacoustic inversion. Warping also has potential applications in other specialties, including bioacoustics; however, the technique can be daunting to many potential users unfamiliar with its intricacies. Consequently, this tutorial article covers basic warping theory, presents simulation examples, and provides practical experimental strategies. Accompanying supplementary material provides matlab code and simulated and experimental datasets for easy implementation of warping on both impulsive and frequency-modulated signals from both biotic and man-made sources. This combined material should provide interested readers with user-friendly resources for implementing warping methods into their own research.This work was supported by the Office of Naval Research (Task Force Ocean, project N00014-19-1-2627) and by the North Pacific Research Board (project 1810). Original warping developments were supported by the French Delegation Generale de l'Armement

A Fast Blind Impulse Detector for Bernoulli-Gaussian Noise in Underspread Channel

The Bernoulli-Gaussian (BG) model is practical to characterize impulsive noises that widely exist in various communication systems. To estimate the BG model parameters from noise measurements, a precise impulse detection is essential. In this paper, we propose a novel blind impulse detector, which is proven to be fast and accurate for BG noise in underspread communication channels.Comment: v2 to appear in IEEE ICC 2018, Kansas City, MO, USA, May 2018 Minor erratums added in v

An array-based receiver function deconvolution method: methodology and application

Receiver functions (RFs) estimated on dense arrays have been widely used for the study of Earth structures across multiple scales. However, due to the ill-posedness of deconvolution, RF estimation faces challenges such as non-uniqueness and data overfitting. In this paper, we present an array-based RF deconvolution method in the context of emerging dense arrays. We propose to exploit the wavefield coherency along a dense array by joint inversions of waveforms from multiple events and stations for RFs with a minimum number of phases required by data. The new method can effectively reduce the instability of deconvolution and help retrieve RFs with higher fidelity. We test the algorithm on synthetic waveforms and show that it produces RFs with higher interpretability than those by the conventional RF estimation practice. Then we apply the method to real data from the 2016 Incorporated Research Institutions for Seismology (IRIS) community wavefield experiment in Oklahoma and are able to generate high-resolution RF profiles with only three teleseismic earthquakes recorded by the temporary deployment. This new method should help enhance RF images derived from short-term high-density seismic profiles

An array-based receiver function deconvolution method: methodology and application

Receiver functions (RFs) estimated on dense arrays have been widely used for the study of Earth structures across multiple scales. However, due to the ill-posedness of deconvolution, RF estimation faces challenges such as non-uniqueness and data overfitting. In this paper, we present an array-based RF deconvolution method in the context of emerging dense arrays. We propose to exploit the wavefield coherency along a dense array by joint inversions of waveforms from multiple events and stations for RFs with a minimum number of phases required by data. The new method can effectively reduce the instability of deconvolution and help retrieve RFs with higher fidelity. We test the algorithm on synthetic waveforms and show that it produces RFs with higher interpretability than those by the conventional RF estimation practice. Then we apply the method to real data from the 2016 Incorporated Research Institutions for Seismology (IRIS) community wavefield experiment in Oklahoma and are able to generate high-resolution RF profiles with only three teleseismic earthquakes recorded by the temporary deployment. This new method should help enhance RF images derived from short-term high-density seismic profiles

Spherical deconvolution of multichannel diffusion MRI data with non-Gaussian noise models and spatial regularization

Spherical deconvolution (SD) methods are widely used to estimate the intra-voxel white-matter fiber orientations from diffusion MRI data. However, while some of these methods assume a zero-mean Gaussian distribution for the underlying noise, its real distribution is known to be non-Gaussian and to depend on the methodology used to combine multichannel signals. Indeed, the two prevailing methods for multichannel signal combination lead to Rician and noncentral Chi noise distributions. Here we develop a Robust and Unbiased Model-BAsed Spherical Deconvolution (RUMBA-SD) technique, intended to deal with realistic MRI noise, based on a Richardson-Lucy (RL) algorithm adapted to Rician and noncentral Chi likelihood models. To quantify the benefits of using proper noise models, RUMBA-SD was compared with dRL-SD, a well-established method based on the RL algorithm for Gaussian noise. Another aim of the study was to quantify the impact of including a total variation (TV) spatial regularization term in the estimation framework. To do this, we developed TV spatially-regularized versions of both RUMBA-SD and dRL-SD algorithms. The evaluation was performed by comparing various quality metrics on 132 three-dimensional synthetic phantoms involving different inter-fiber angles and volume fractions, which were contaminated with noise mimicking patterns generated by data processing in multichannel scanners. The results demonstrate that the inclusion of proper likelihood models leads to an increased ability to resolve fiber crossings with smaller inter-fiber angles and to better detect non-dominant fibers. The inclusion of TV regularization dramatically improved the resolution power of both techniques. The above findings were also verified in brain data

Post-WMAP Assessment of Infrared Cutoff in the Primordial Spectrum from Inflation

The recent Cosmic Microwave Background (CMB) measurements indicate that there is power deficiency of the CMB anisotropies at large scales compared with the $\Lambda$CDM model. We have investigated the possibility of explaining such effects by a class of primordial power spectra which have infrared cutoffs close to the horizon scale. The primordial power spectrum recovered by direct deconvolution of the observed CMB angular spectrum indicates that the data prefers a sharp infrared cutoff with a localized excess (bump) just above the cutoff. We have been motivated to assess plausible extensions of simplest inflationary scenarios which readily accommodate similar form of infrared cutoff. We carry out a complete Bayesian analysis of the parameter space using {\it Markov Chain Monte Carlo} technique with such a class of primordial power spectra. We show that primordial power spectrum that have features such as an infrared cutoff followed by a subsequent excess in power give better fit to the observed data compared to a nearly scale-invariant power law or power spectrum with just a monotonic infrared cutoff. However, there is substantial room for improvement in the match to data and calls for exploration of other mechanisms that may lead to infrared cutoff even closer to that recovered by direct deconvolution approach.Comment: Changes to match version accepted for publication in PR

Dynamic Occupancy Grid Prediction for Urban Autonomous Driving: A Deep Learning Approach with Fully Automatic Labeling

Long-term situation prediction plays a crucial role in the development of intelligent vehicles. A major challenge still to overcome is the prediction of complex downtown scenarios with multiple road users, e.g., pedestrians, bikes, and motor vehicles, interacting with each other. This contribution tackles this challenge by combining a Bayesian filtering technique for environment representation, and machine learning as long-term predictor. More specifically, a dynamic occupancy grid map is utilized as input to a deep convolutional neural network. This yields the advantage of using spatially distributed velocity estimates from a single time step for prediction, rather than a raw data sequence, alleviating common problems dealing with input time series of multiple sensors. Furthermore, convolutional neural networks have the inherent characteristic of using context information, enabling the implicit modeling of road user interaction. Pixel-wise balancing is applied in the loss function counteracting the extreme imbalance between static and dynamic cells. One of the major advantages is the unsupervised learning character due to fully automatic label generation. The presented algorithm is trained and evaluated on multiple hours of recorded sensor data and compared to Monte-Carlo simulation