Geostatistical Seismic Inversion for Temperature and Salinity in the Madeira Abyssal Plain

A two-dimensional multichannel seismic reflection profile acquired in the Madeira Abyssal Plain during June 2016 was used in a modelling workflow comprising seismic oceanography processing, geostatistical inversion and Bayesian classification to predict the probability of occurrence of distinct water masses. The seismic section was processed to image in detail the fine scale structure of the water column using seismic oceanography. The processing sequence was developed to preserve, as much as possible, the relative seismic amplitudes of the data and enhance the shallow structure of the water column by effectively suppressing the direct arrival. The migrated seismic oceanography section shows an eddy at the expected Mediterranean Outflow Water depths, steeply dipping reflectors, which indicate the possible presence of frontal activity or secondary dipping eddy structures, and strong horizontal reflections between intermediate water masses suggestive of double diffuse processes. We then developed and applied an iterative geostatistical seismic oceanography inversion methodology to predict the spatial distribution of temperature and salinity. Due to the lack of contemporaneous direct measurements of temperature and salinity we used a global ocean model as spatial constraint during the inversion and nearby contemporaneous ARGO data to infer the expected statistical properties of both model parameters. After the inversion, Bayesian classification was applied to all temperature and salinity models inverted during the last iteration to predict the spatial distribution of three distinct water masses. A preliminary interpretation of these probabilistic models agrees with the expected ocean dynamics of the region.