Advancing ocean subsurface thermal structure estimation in the Pacific Ocean: A multi-model ensemble machine learning approach

Estimation of the ocean subsurface thermal structure (OSTS) is important for understanding thermodynamic processes and climate variability. In the present study, a novel multi-model ensemble machine learning (Ensemble-ML) model is developed to retrieve subsurface thermal structure in the Pacific Ocean by integrating sea surface data with Argo observations. The Ensemble-ML model integrates four individual machine learning models to enhance estimation accuracy and reliability. Our results exhibit good agreement between the satellite sea surface temperature (SST) and sea surface salinity (SSS) data and Argo observations, providing validation for the utilization of these datasets in the Ensemble-ML model. The Ensemble-ML model exhibits better performance compared to individual machine learning models, with an average root mean square error (RMSE) of 0.3273 °C and an average coefficient of determination (R²) of 0.9905. Notably, incorporating geographical information as input variables enhance model performance, emphasizing the importance of considering spatial context in OSTS estimation. The Ensemble-ML model accurately captures the spatial distribution of OSTS across depths and seasons in the Pacific Ocean, effectively reproducing critical temperature features while maintaining strong agreement with Argo observations. Nevertheless, its performance shows relative weakness within the thermocline layer and the equatorial Pacific region (spanning from 10°S to 10°N latitude), which are characterized by complex circulation systems. Despite these challenges, the Ensemble-ML model effectively reproduces the spatial distribution of OSTS of the Pacific Ocean. This indicates the potential of machine learning models, particularly ensemble models, for enhancing OSTS estimation in the Pacific Ocean and other regions, offering valuable insights for future research and applications in physical oceanography.