Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength

In this study, we present a new method called <i>BathySent</i> to retrieve shallow bathymetry from space that is based on the joint measurement of ocean wave celerity (c) and wavelength (λ). We developed the method to work with Sentinel 2 data, exploiting the time lag between two Sentine...

Full description

Bibliographic Details
Main Authors: Marcello de Michele, Daniel Raucoules, Deborah Idier, Farid Smai, Michael Foumelis
Format: Article
Language:English
Published: MDPI AG 2021-05-01
Series:Remote Sensing
Subjects:
Online Access:https://www.mdpi.com/2072-4292/13/11/2149
_version_ 1827691044525834240
author Marcello de Michele
Daniel Raucoules
Deborah Idier
Farid Smai
Michael Foumelis
author_facet Marcello de Michele
Daniel Raucoules
Deborah Idier
Farid Smai
Michael Foumelis
author_sort Marcello de Michele
collection DOAJ
description In this study, we present a new method called <i>BathySent</i> to retrieve shallow bathymetry from space that is based on the joint measurement of ocean wave celerity (c) and wavelength (λ). We developed the method to work with Sentinel 2 data, exploiting the time lag between two Sentinel 2 spectral bands, acquired quasi-simultaneously, from a single satellite dataset. Our method was based on the linear dispersion law, which related water depth to wave celerity and wavelength: when the water depth was less than about half the dominant wavelength, the wave celerity and wavelength decreased due to decreasing water depth (h) as the waves propagated towards the coast. Instead of using a best weighted <inline-formula><math display="inline"><semantics><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">c</mi><mo>,</mo><mi mathvariant="sans-serif">λ</mi><mo stretchy="false">)</mo></mrow></semantics></math></inline-formula> fit with the linear dispersion relation to retrieve h, we proposed solving the linear dispersion relation for each <inline-formula><math display="inline"><semantics><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">c</mi><mo>,</mo><mi mathvariant="sans-serif">λ</mi><mo stretchy="false">)</mo></mrow></semantics></math></inline-formula> pair to find multiple h-values within the same resolution cell. Then, we calculated the weighted averaged h-value for each resolution cell. To improve the precision of the final bathymetric map, we stacked the bathymetry values from <i>N</i>-different datasets acquired from the same study area on different dates. We first tested the algorithm on a set of images representing simulated ocean waves, then we applied it to a real set of Sentinel 2 data obtained of our study area, Gâvres peninsula (France, 47°,67 lat.; −3°35 lon.). A comparison with in situ bathymetry yielded good results from the synthetic images (r<sup>2</sup> = 0.9) and promising results with the Sentinel 2 images (r<sup>2</sup> = 0.7) in the 0–16 m depth zone.
first_indexed 2024-03-10T10:52:31Z
format Article
id doaj.art-4561a1e092e04fa4a51aebb1bb5c440f
institution Directory Open Access Journal
issn 2072-4292
language English
last_indexed 2024-03-10T10:52:31Z
publishDate 2021-05-01
publisher MDPI AG
record_format Article
series Remote Sensing
spelling doaj.art-4561a1e092e04fa4a51aebb1bb5c440f2023-11-21T22:07:40ZengMDPI AGRemote Sensing2072-42922021-05-011311214910.3390/rs13112149Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and WavelengthMarcello de Michele0Daniel Raucoules1Deborah Idier2Farid Smai3Michael Foumelis4BRGM, French Geological Survey, Direction des Risques et de la Prevention des Risques—3 av. C. Guillemin, 45000 Orleans, FranceBRGM, French Geological Survey, Direction des Risques et de la Prevention des Risques—3 av. C. Guillemin, 45000 Orleans, FranceBRGM, French Geological Survey, Direction des Risques et de la Prevention des Risques—3 av. C. Guillemin, 45000 Orleans, FranceBRGM, French Geological Survey, Direction des Risques et de la Prevention des Risques—3 av. C. Guillemin, 45000 Orleans, FranceBRGM, French Geological Survey, Direction des Risques et de la Prevention des Risques—3 av. C. Guillemin, 45000 Orleans, FranceIn this study, we present a new method called <i>BathySent</i> to retrieve shallow bathymetry from space that is based on the joint measurement of ocean wave celerity (c) and wavelength (λ). We developed the method to work with Sentinel 2 data, exploiting the time lag between two Sentinel 2 spectral bands, acquired quasi-simultaneously, from a single satellite dataset. Our method was based on the linear dispersion law, which related water depth to wave celerity and wavelength: when the water depth was less than about half the dominant wavelength, the wave celerity and wavelength decreased due to decreasing water depth (h) as the waves propagated towards the coast. Instead of using a best weighted <inline-formula><math display="inline"><semantics><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">c</mi><mo>,</mo><mi mathvariant="sans-serif">λ</mi><mo stretchy="false">)</mo></mrow></semantics></math></inline-formula> fit with the linear dispersion relation to retrieve h, we proposed solving the linear dispersion relation for each <inline-formula><math display="inline"><semantics><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">c</mi><mo>,</mo><mi mathvariant="sans-serif">λ</mi><mo stretchy="false">)</mo></mrow></semantics></math></inline-formula> pair to find multiple h-values within the same resolution cell. Then, we calculated the weighted averaged h-value for each resolution cell. To improve the precision of the final bathymetric map, we stacked the bathymetry values from <i>N</i>-different datasets acquired from the same study area on different dates. We first tested the algorithm on a set of images representing simulated ocean waves, then we applied it to a real set of Sentinel 2 data obtained of our study area, Gâvres peninsula (France, 47°,67 lat.; −3°35 lon.). A comparison with in situ bathymetry yielded good results from the synthetic images (r<sup>2</sup> = 0.9) and promising results with the Sentinel 2 images (r<sup>2</sup> = 0.7) in the 0–16 m depth zone.https://www.mdpi.com/2072-4292/13/11/2149bathymetrySentinel 2cross correlation
spellingShingle Marcello de Michele
Daniel Raucoules
Deborah Idier
Farid Smai
Michael Foumelis
Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength
Remote Sensing
bathymetry
Sentinel 2
cross correlation
title Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength
title_full Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength
title_fullStr Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength
title_full_unstemmed Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength
title_short Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength
title_sort shallow bathymetry from multiple sentinel 2 images via the joint estimation of wave celerity and wavelength
topic bathymetry
Sentinel 2
cross correlation
url https://www.mdpi.com/2072-4292/13/11/2149
work_keys_str_mv AT marcellodemichele shallowbathymetryfrommultiplesentinel2imagesviathejointestimationofwavecelerityandwavelength
AT danielraucoules shallowbathymetryfrommultiplesentinel2imagesviathejointestimationofwavecelerityandwavelength
AT deborahidier shallowbathymetryfrommultiplesentinel2imagesviathejointestimationofwavecelerityandwavelength
AT faridsmai shallowbathymetryfrommultiplesentinel2imagesviathejointestimationofwavecelerityandwavelength
AT michaelfoumelis shallowbathymetryfrommultiplesentinel2imagesviathejointestimationofwavecelerityandwavelength