Berta Biescas

Modelling seismic oceanography experiments by using first- and second-order Complex Frequency Shifted Perfectly Matched Layers

This work investigates the ability of modelling seismic oceanography experiments by using underwater acoustic propagation equations. Seismic oceanography tries to retrieve the fine structure of the ocean water masses by processing the acoustic waves reflected in the low-contrast interfaces of fronts, eddies, internal waves or thermohaline intrusions. Since the reflectivity of such interfaces is of order 10 ―3 ―10 ―4 , the absorption capability of the numerical boundaries becomes crucial. Complex Frequency Shifted offers a better alternative to classical Perfectly Matched Layer formulation, but has not yet been extended to acoustic equations. Here, first- and second-order Complex Frequency Shifted Perfectly Matched Layers equations are proposed which can provide reflection coefficients of order 10 ―5 . Therefore, a numerical Finite-Difference Time-Domain (FDTD) scheme combined with the proposed CFS-PML equations is able to model such experiments.

Characterization of the sub-mesoscale energy cascade in the Alboran Sea thermocline from spectral analysis of high-resolution MCS data

Part of the kinetic energy that maintains ocean circulation cascades down to small scales until it is dissipated through mixing. While most steps of this downward energy cascade are well understood, an observational gap exists at horizontal scales of 103-101 m that prevents characterizing a key step in the chain: the transition from anisotropic internal wave motions to isotropic turbulence. Here we show that this observational gap can be covered using high-resolution multichannel seismic (HR-MCS) data. Spectral analysis of acoustic reflectors imaged in the Alboran Sea thermocline shows that this transition is likely caused by shear instabilities. In particular, we show that the averaged horizontal wavenumber spectra of the reflectors vertical displacements display three subranges that reproduce theoretical spectral slopes of internal waves [λx > 100 m], Kelvin-Helmholtz-type shear instabilities [100 m > λx > 33 m], and turbulence [λx < 33 m], indicating that the whole chain of events is occurring continuously and simultaneously in the surveyed area.

Markov Chain Monte Carlo inversion of temperature and salinity structure of an internal solitary wave packet from marine seismic data.

Marine seismic reflection technique is used to observe the strong ocean dynamic process of nonlinear internal solitary waves (ISWs or solitons) in the near-surface water. Analysis of ISWs is problematical because of their transient nature and limitations of classical physical oceanography methods. This work explores a Markov Chain Monte Carlo (MCMC) approach to recover the temperature and salinity of ISW field using the seismic reflectivity data and in situ hydrographic data. The MCMC approach is designed to directly sample the posterior probability distributions of temperature and salinity which are the solutions of the system under investigation. The principle improvement is the capability of incorporating uncertainties in observations and prior models which then provide quantified uncertainties in the output model parameters. We tested the MCMC approach on two acoustic reflectivity data sets one synthesized from a CTD cast and the other derived from multichannel seismic reflections. This method finds the solutions faithfully within the significantly narrowed confidence intervals from the provided priors. Combined with a low frequency initial model interpreted from seismic horizons of ISWs, the MCMC method is used to compute the finescale temperature, salinity, acoustic velocity, and density of ISW field. The statistically derived results are equivalent to the conventional linearized inversion method. However, the former provides us the quantified uncertainties of the temperature and salinity along the whole section whilst the latter does not. These results are the first time ISWs have been mapped with sufficient detail for further analysis of their dynamic properties.

Fine-scale thermohaline ocean structure retrieved with 2-D prestack full-waveform inversion of multichannel seismic data: Application to the Gulf of Cadiz (SW Iberia)

This work demonstrates the feasibility of 2D time-domain, adjoint-state acoustic full-waveform inversion (FWI) to retrieve high-resolution models of ocean physical parameters such as sound speed, temperature and salinity. The proposed method is first described and then applied to pre-stack multi-channel seismic (MCS) data acquired in the Gulf of Cadiz (SW Iberia) in 2007 in the framework of the Geophysical Oceanography project. The inversion strategy flow includes specifically-designed data pre-conditioning for acoustic noise reduction, followed by the inversion of sound speed in the shotgather domain. We show that the final sound speed model has a horizontal resolution of ∼ 70m, which is two orders of magnitude better than that of the initial model constructed with coincident eXpendable Bathy Thermograph (XBT) data, and close to the theoretical resolution of O(λ). Temperature (T) and salinity (S) are retrieved with the same lateral resolution as sound speed by combining the inverted sound speed model with the thermodynamic equation of seawater and a local, depth-dependent T-S relation derived from regional conductivity-temperature-depth (CTD) measurements of the National Oceanic and Atmospheric Administration (NOAA) database. The comparison of the inverted T and S models with XBT and CTD casts deployed simultaneously to the MCS acquisition shows that the thermohaline contrasts are resolved with an accuracy of 0.18oC for temperature and 0.08 PSU for salinity. The combination of oceanographic and MCS data into a common, pseudo-automatic inversion scheme allows to quantitatively resolve submeso-scale features that ought to be incorporated into larger-scale ocean models of oceans structure and circulation. This article is protected by copyright. All rights reserved.

Synthetic Modeling for an Acoustic Exploration System for Physical Oceanography

AbstractMarine multichannel seismic (MCS) data, used to obtain structural reflection images of the earth’s subsurface, can also be used in physical oceanography exploration. This method provides vertical and lateral resolutions of O(10–100) m, covering the existing observational gap in oceanic exploration. All MCS data used so far in physical oceanography studies have been acquired using conventional seismic instrumentation originally designed for geological exploration. This work presents the proof of concept of an alternative MCS system that is better adapted to physical oceanography and has two goals: 1) to have an environmentally low-impact acoustic source to minimize any potential disturbance to marine life and 2) to be light and portable, thus being installed on midsize oceanographic vessels. The synthetic experiments simulate the main variables of the source, shooting, and streamer involved in the MCS technique. The proposed system utilizes a 5-s-long exponential chirp source of 208 dB relative to ...