Multiscale multiphysics data-informed modeling for three-dimensional ocean acoustic simulation and prediction

Author Posting. © Acoustical Society of America, 2019. This article is posted here by permission of Acoustical Society of America for personal use, not for redistribution. The definitive version was published in Journal of the Acoustical Society of America 146(3), (2019): 1996-2015, doi:10.1121/1.5126012.Three-dimensional (3D) underwater sound field computations have been used for a few decades to understand sound propagation effects above sloped seabeds and in areas with strong 3D temperature and salinity variations. For an approximate simulation of effects in nature, the necessary 3D sound-speed field can be made from snapshots of temperature and salinity from an operational data-driven regional ocean model. However, these models invariably have resolution constraints and physics approximations that exclude features that can have strong effects on acoustics, example features being strong submesoscale fronts and nonhydrostatic nonlinear internal waves (NNIWs). Here, work to predict NNIW fields to improve 3D acoustic forecasts using an NNIW model nested in a tide-inclusive data-assimilating regional model is reported. The work was initiated under the Integrated Ocean Dynamics and Acoustics project. The project investigated ocean dynamical processes that affect important details of sound-propagation, with a focus on those with strong intermittency (high kurtosis) that are challenging to predict deterministically. Strong internal tides and NNIW are two such phenomena, with the former being precursors to NNIW, often feeding energy to them. Successful aspects of the modeling are reported along with weaknesses and unresolved issues identified in the course of the work.This work was supported by Department of Defense Multidisciplinary University Initiative (MURI) Grant No. N00014-11-1-0701, managed by the Office of Naval Research Ocean Acoustics Program, and National Science Foundation Grant No. OCE-1060430. Final manuscript preparation was supported by ONR Ocean Acoustics Grant Nos. N00014-17-1-2624 and N00014-17-1-2692. P.F.J.L. also thanks ONR and NSF for research support under Grant Nos. N00014-13-1-0518 (Multi-DA) and OCE-1061160 (ShelfIT) to MIT, respectively. The MSEAS-based series of simulations for the New Jersey shelf region examined here was accelerated toward completion by the interest in realistic 3D acoustic fields expressed by Dr. Ivars Kirsteins at the Naval Undersea Warfare Center.2020-03-3

A three-dimensional parabolic equation model of sound propagation using higher-order operator splitting and Padé approximants

Author Posting. © Acoustical Society of America, 2012. This article is posted here by permission of Acoustical Society of America for personal use, not for redistribution. The definitive version was published in Journal of the Acoustical Society of America 132 (2012): EL364-EL370, doi:10.1121/1.4754421.An alternating direction implicit (ADI) three-dimensional fluid parabolic equation solution method with enhanced accuracy is presented. The method uses a square-root Helmholtz operator splitting algorithm that retains cross-multiplied operator terms that have been previously neglected. With these higher-order cross terms, the valid angular range of the parabolic equation solution is improved. The method is tested for accuracy against an image solution in an idealized wedge problem. Computational efficiency improvements resulting from the ADI discretization are also discussed.This work was sponsored by the Office of Naval Research under Grant Nos. N00014-10-1-0040 and N00014-11-1-0701

Two-Spinon and Orbital Excitations of the Spin-Peierls System TiOCl

We combine high-resolution resonant inelastic x-ray scattering with cluster calculations utilizing a recently derived effective magnetic scattering operator to analyze the polarization, excitation energy, and momentum dependent excitation spectrum of the low-dimensional quantum magnet TiOCl in the range expected for orbital and magnetic excitations (0 - 2.5 eV). Ti 3d orbital excitations yield complete information on the temperature-dependent crystal-field splitting. In the spin-Peierls phase we observe a dispersive two-spinon excitation and estimate the inter- and intra-dimer magnetic exchange coupling from a comparison to cluster calculations

Phylogenetic relationships of cone snails endemic to Cabo Verde based on mitochondrial genomes

Background: Due to their great species and ecological diversity as well as their capacity to produce hundreds of different toxins, cone snails are of interest to evolutionary biologists, pharmacologists and amateur naturalists alike. Taxonomic identification of cone snails still relies mostly on the shape, color, and banding patterns of the shell. However, these phenotypic traits are prone to homoplasy. Therefore, the consistent use of genetic data for species delimitation and phylogenetic inference in this apparently hyperdiverse group is largely wanting. Here, we reconstruct the phylogeny of the cones endemic to Cabo Verde archipelago, a well-known radiation of the group, using mitochondrial (mt) genomes. Results: The reconstructed phylogeny grouped the analyzed species into two main clades, one including Kalloconus from West Africa sister to Trovaoconus from Cabo Verde and the other with a paraphyletic Lautoconus due to the sister group relationship of Africonus from Cabo Verde and Lautoconus ventricosus from Mediterranean Sea and neighboring Atlantic Ocean to the exclusion of Lautoconus endemic to Senegal (plus Lautoconus guanche from Mauritania, Morocco, and Canary Islands). Within Trovaoconus, up to three main lineages could be distinguished. The clade of Africonus included four main lineages (named I to IV), each further subdivided into two monophyletic groups. The reconstructed phylogeny allowed inferring the evolution of the radula in the studied lineages as well as biogeographic patterns. The number of cone species endemic to Cabo Verde was revised under the light of sequence divergence data and the inferred phylogenetic relationships. Conclusions: The sequence divergence between continental members of the genus Kalloconus and island endemics ascribed to the genus Trovaoconus is low, prompting for synonymization of the latter. The genus Lautoconus is paraphyletic. Lautoconus ventricosus is the closest living sister group of genus Africonus. Diversification of Africonus was in allopatry due to the direct development nature of their larvae and mainly triggered by eustatic sea level changes during the Miocene-Pliocene. Our study confirms the diversity of cone endemic to Cabo Verde but significantly reduces the number of valid species. Applying a sequence divergence threshold, the number of valid species within the sampled Africonus is reduced to half.Spanish Ministry of Science and Innovation [CGL2013-45211-C2-2-P, CGL2016-75255-C2-1-P, BES-2011-051469, BES-2014-069575, Doctorado Nacional-567]info:eu-repo/semantics/publishedVersio

An adaptive Metropolis-Hastings scheme: sampling and optimization

We propose an adaptive Metropolis-Hastings algorithm in which sampled data are used to update the proposal distribution. We use the samples found by the algorithm at a particular step to form the information-theoretically optimal mean-field approximation to the target distribution, and update the proposal distribution to be that approximatio. We employ our algorithm to sample the energy distribution for several spin-glasses and we demonstrate the superiority of our algorithm to the conventional MH algorithm in sampling and in annealing optimization.Comment: To appear in Europhysics Letter

Time-evolving acoustic propagation modeling in a complex ocean environment

During naval operations, sonar performance estimates often need to be computed in-situ with limited environmental information. This calls for the use of fast acoustic propagation models. Many naval operations are carried out in challenging and dynamic environments. This makes acoustic propagation and sonar performance behavior particularly complex and variable, and complicates prediction. Using data from a field experiment, we have investigated the accuracy with which acoustic propagation loss (PL) can be predicted, using only limited modeling capabilities. Environmental input parameters came from various sources that may be available in a typical naval operation. The outer continental shelf shallow-water experimental area featured internal tides, packets of nonlinear internal waves, and a meandering water mass front. For a moored source/receiver pair separated by 19.6 km, the acoustic propagation loss for 800 Hz pulses was computed using the peak amplitude. The variations in sound speed translated into considerable PL variability of order 15 dB. Acoustic loss modeling was carried out using a data-driven regional ocean model as well as measured sound speed profile data for comparison. The acoustic model used a two-dimensional parabolic approximation (vertical and radial outward wavenumbers only). The variance of modeled propagation loss was less than that measured. The effect of the internal tides and sub-tidal features was reasonably well modeled; these made use of measured sound speed data. The effects of nonlinear waves were not well modeled, consistent with their known three-dimensional effects but also with the lack of measurements to initialize and constrain them.Netherlands. Ministry of DefenceUnited States. Office of Naval Research (Grant N00014-12-1-0944 (ONR6.2))United States. Office of Naval Research (Grant N00014-08-1-1097 (ONR6.1))United States. Office of Naval Research (Grant N00014-08-1-0680 (PLUS-SEAS)

Dark Matter Direct Detection Signals inferred from a Cosmological N-body Simulation with Baryons

We extract at redshift z=0 a Milky Way sized object including gas, stars and dark matter (DM) from a recent, high-resolution cosmological N-body simulation with baryons. Its resolution is sufficient to witness the formation of a rotating disk and bulge at the center of the halo potential. The phase-space structure of the central galactic halo reveals the presence of a dark disk component, that is co-rotating with the stellar disk. At the Earth's location, it contributes to around 25% of the total DM local density, whose value is rho_DM ~ 0.37 GeV/cm^3. The velocity distributions also show strong deviations from pure Gaussian and Maxwellian distributions, with a sharper drop of the high velocity tail. We give a detailed study of the impact of these features on the predictions for DM signals in direct detection experiments. In particular, the question of whether the modulation signal observed by DAMA is or is not excluded by limits set by other experiments (CDMS, XENON and CRESST...) is re-analyzed and compared to the case of a standard Maxwellian halo, in both the elastic and the inelastic scattering scenarios. We find that the compatibility between DAMA and the other experiments is improved. In the elastic scenario, the DAMA modulation signal is slightly enhanced in the so-called channeling region, as a result of several effects. For the inelastic scenario, the improvement of the fit is mainly attributable to the departure from a Maxwellian distribution at high velocity.Comment: 39 page