cs.LGAug 17, 2026

OceanDepths: A Global Dataset of Paired Subsurface and Surface Ocean Observations

Authors: Simon Donike, Ruben Cartuyvels, Antonino Ian Ferola, Elisa Carli, Diego Fernandez Prieto, Marie-Helene Rio

Organizations: Image Processing Laboratory, University of Valencia, Av. de Blasco Ibáñez 13, 46010, Spain · ESA-ESRIN, Via Galileo Galilei 1, 00044 Frascati RM, Italy

Read PDF
Ask AI
Cite

Copy citation

DOINot available for this paper

Abstract

Despite comprising over 70% of its surface, the world's oceans are critically underobserved compared to the land surface or the atmosphere. Understanding the global ocean requires jointly observing its surface and subsurface structure, yet no standardized, high-resolution dataset couples satellite surface fields to co-located in situ depth profiles in an AI-ready format. Existing resources either consist of model-reconstructed gridded products rather than observations, cover only a single variable or basin, or operate at resolutions too coarse for mesoscale dynamics. We introduce OceanDepths, the first open, global, regridded AI-ready dataset that pairs satellite-derived sea surface temperature (SST), sea surface salinity (SSS), and sea surface height (SSH) L4 products with co-located EN4 subsurface temperature and salinity profiles, complemented by matched GLORYS12 ocean reanalysis data to support comparisons or multi-stage learning. The dataset spans 2000-2024 at 0.1 degrees x 0.1 degrees spatial resolution and at weekly temporal resolution, covering the entire globe's sea surface and with over 9.5 million paired profiles interpolated to 50 standardized depth levels. We provide a configurable system to split the globe in equally sized spatial patches. The 4D multivariate structure, high resolution, long temporal extent, and extreme sparsity of subsurface observations (approximately 0.01% per depth level) make OceanDepths a challenging testbed for novel AI methods. We demonstrate subsurface state reconstruction as an example task with simple baseline models, but also envision OceanDepths to support the development of observation-based forecast methods and other related tasks. Available at: https://huggingface.co/datasets/ESA-philab/OceanDepths.

Figures & tables

Figure 1: Global distribution of the EN4 density per week and examples of surface observations and reanalysis data. OceanDepths aggregates all profiles gathered between 2000 and 2024 with sea surface temperature, salinity and height observational products accompanied by the GLORYS12V1 reanalysis product. The SST, ADT, SSS, GLORYS12, and profile T and S are shown for the white square around New York on the map in the middle.
DatasetSurfaceSubsurf.Depth (m)LevelsSRTRPeriod
Interpolation / Analysis
ARMOR3D [ 13 ]–T, S1500501/8 ∘D1993–…
ISAS [ 15 ]–T, S55001870.5 ∘M2002–2020
WOA [ 19 ]–T, S, B55001020.25–1 ∘M1971–2022
IAPv4 [ 6 ]–T60001191 ∘M1940–2023
RG-Clim [ 24 ]–T, S2000581 ∘M2004–2018
Table 1: Comparison of existing public global ocean subsurface–surface datasets. Variables: T = temperature, S = salinity, U/V = currents, H = sea surface height or ADT (absolute dynamic topography), SST = sea surface temperature, SSS = sea surface salinity, W = wind, B = biological variables, nutrients, oxygen, SR/TR = spatial/temporal resolution G = Global. † : Preserves native sensor resolutions (2 km–0.25 ∘ ) without resampling to a common grid.
PropertyValue
Spatial coverageGlobal, 3,600 ×\times 1,800 px at 0.1 \SIUnitSymbolDegree0.1\text{\,}\mathrm{\SIUnitSymbolDegree} resolution (4,338,138 ocean cells)
Temporal coverageJanuary 2000 – July 2024, weekly resolution (1283 target weeks)
Total Argo profiles9,485,977 (9.4 M valid temperature; 7.0 M valid salinity)
Surface variables per sampleSST (OSTIA), SSS & surface density (MULTIOBS), ADT (DUACS)
Depth levels50 GLORYS standard levels (0.49 m to 5,728 m)
Patch dataset358 K patch–date samples ( 1282128^{2} px, no overlap) – 343 K train / 15 K eval
Table 2: Dataset summary statistics.
Figure 2: EN4–GLORYS12 profile alignment example. Example from the South Atlantic in September 2018. The depth-interpolated profiles closely track the GLORYS12 reanalysis across the full water column for both modalities, with occasional discrepancies (which is expected due to limitations of the reanalysis method and/or uncertainties in the profile measurements).
Figure 3: Global overlapping patch grid coverage for the training dataset. Retained patches ( ≤\leq 30% land, transparent blue) and force-included patches for enclosed basins (in red boxes) cover all major ocean regions. Discarded patches ( >> 30% land) are excluded. Transparency indicates overlap density from the 75% overlap tiling.
Figure 4: Patch land-fraction filtering with land areas (black), EN4 locations and SST. a: Kept by default. b: Argentinian coast, rejected by default. c: Aegean Sea, included due to relaxed region-based rules.
 TemperatureSalinity
MethodRMSEMAER2R^{2}RMSEMAER2R^{2}
Climatology0.9740.5100.9640.5430.1700.717
IDW0.9790.5210.9640.5660.1740.678
LSTM2.4201.6450.7010.6230.2950.633
1D CNN3.0732.3770.3700.7970.5360.338
3D U-Net1.1010.6450.9370.5180.1770.738
Table 3: Subsurface reconstruction results (evaluation year 2018, week 25, ∣Pval∣=0.2∣P∣|\mathcal{P}_{\mathrm{val}}|=0.2|\mathcal{P}| ). Performance is reported for temperature and salinity against held-out EN4 profiles, averaged across depth levels no deeper than 2000 m. Best values highlighted in bold.
Figure 5: RMSE per depth against held-out EN4 samples. T and S reconstruction error by depth level, averaged over the global 2018 week-25 validation set.
Appendix figures & tables2 assets

Supplementary material from the paper’s appendix.

Appendix

Figure 6: Random selection of surface-level training patches from the dataset, illustrating the diversity of ocean regions, SST gradients, and coastline configurations encountered during training.
Figure 7: EN4 depth index compared to depth in meters. Heatmap of corrected depth (m, log scale) vs. EN4 level index, with median (white) and P10/P90 (pink dotted lines) curves. The irregular, profile-dependent depth sampling motivates the interpolation onto the fixed 50-level GLORYS12 coordinate.

Explore similar work

CardsList
Apr 3, 2026

High-resolution probabilistic estimation of three-dimensional regional ocean dynamics from sparse surface observations

The ocean interior regulates Earth's climate but remains sparsely observed due to limited in situ measurements, while satellite observations are restricted to the surface. We present a depth-aware generative framework for reconstructing high-resolution three-dimensional ocean states from extremely sparse surface data. Our approach employs a conditional denoising diffusion probabilistic model (DDPM) trained on sea surface height and temperature observations with up to 99.9 percent sparsity, without reliance on a background dynamical model. By incorporating continuous depth embeddings, the model learns a unified vertical representation of the ocean states and generalizes to previously unseen depths. Applied to the Gulf of Mexico, the framework accurately reconstructs subsurface temperature, salinity, and velocity fields across multiple depths. Evaluations using statistical metrics, spectral analysis, and heat transport diagnostics demonstrate recovery of both large-scale circulation and multiscale variability. These results establish generative diffusion models as a scalable approach for probabilistic ocean reconstruction in data-limited regimes, with implications for climate monitoring and forecasting.

Niloofar Asefi, Tianning Wu, Ruoying He +1

Department of Electrical and Computer Engineering, University of California, Santa Cruz, CA 95064, USA · Department of Marine, Earth, and Atmospheric Sciences, North Carolina State University, Raleigh, NC 27695, USA · Department of Applied Mathematics, University of California, Santa Cruz, CA 95064, USA

Conditional Diffusion Models
AskPDF
Apr 21, 2026

An Adaptive Spatiotemporal Clustering Framework for 3D Ocean Subsurface Temperature Reconstruction

The reconstruction of ocean subsurface temperature (OST) using satellite remote sensing data holds significant scientific value for advancing the understanding of ocean dynamics and climate variability. However, the scarcity of subsurface observations, combined with the high degree of nonlinearity and spatiotemporal heterogeneity in subsurface processes, poses substantial challenges to the accuracy and generalization capability of traditional reconstruction methods. To address these limitations, this study proposes an adaptive framework that could capture both vertical structural dependencies and temporal variation patterns of OST via spatio-temporal clustering. By incorporating this framework with various deep learning models, e.g., dual-path convolutional neural networks (DP-CNN), Attention U-Net, and Vision Transformer (ViT), the OST field can be accurately reconstructed at a global scale only using surface observations, i.e., sea surface temperature (SST), sea surface salinity (SSS), sea surface height (SSH), and sea surface wind (SSW). Experimental results demonstrate that multiple deep learning methods using the proposed framework largely outperform their original counterparts, yielding improvements in RMSE ranging from 12.4% to 27.2%. This study provides a reliable solution for subsurface temperature reconstruction, offering important implications for meteorological modeling and climate change assessment.

Ming Shan Loo, Wengen Li, Xudong Jiang +4

School of Computer Science and Technology, Tongji University, Shanghai, China · The Key Laboratory of Road and Traffic Engineering, Ministry of Education, College of Transportation, Tongji University, Shanghai, China · Project Management Office of China National Scientific Seafloor Observatory, Tongji University, Shanghai, China

Clustering
AskPDF
May 11, 2026

AxiomOcean: Forecasting the Three-Dimensional Structure of the Upper Ocean

Short-term ocean forecast skill depends strongly on the three-dimensional ocean structure of the upper ocean, which governs stratification, subsurface heat storage, and the response of the ocean to atmospheric forcing. However, AI ocean forecasting models often fail to preserve this vertical structure, resulting in over-smoothed subsurface features and weak physical consistency under strong forcing. Here, we present AxiomOcean, a global AI ocean forecasting model that explicitly represents vertical hierarchy and cross-layer dependence within the water column. By combining a fully three-dimensional encoder-backbone-decoder architecture with surface atmospheric forcing, AxiomOcean jointly predicts upper-ocean temperature, salinity, and three-dimensional currents at global 1/12° resolution down to 643 m depth. In 10-day forecasts, AxiomOcean outperforms an advanced AI comparison model across variables and lead times, reducing day-1 RMSE by approximately 20 to 35% while maintaining higher anomaly correlation. The gain is not achieved through excessive smoothing: AxiomOcean better preserves eddy kinetic energy, temperature and salinity variance. Its advantage also extends through the water column and remains evident across the equatorial Pacific, Kuroshio Extension, and Southern Ocean, yielding a more realistic reconstruction of upper-ocean heat content. These results show that explicitly preserving upper-ocean three-dimensional structure can improve both forecast accuracy and physical fidelity in AI ocean prediction.

Sensen Wu, Yifan Chen, Guantao Pu +9

School of Earth Sciences & Zhejiang Key Laboratory of Geographic Information Science, Zhejiang University, Hangzhou, China · State Key Laboratory of Ocean Sensing, Zhejiang University, Hangzhou, China · The State Key Lab of Brain-Machine Intelligence, Zhejiang University, Hangzhou, China +1

Spatiotemporal Forecasting
AskPDF