Spatio-temporal methods for estimating subsurface ocean thermal response to tropical cyclones

TL;DR

This paper introduces a comprehensive spatio-temporal modeling framework that accurately estimates subsurface ocean thermal responses to tropical cyclones using in situ data, addressing noise, seasonality, and correlation challenges.

Contribution

It advances previous methods by integrating seasonal mean estimation, Gaussian processes, and nonparametric regression into an ANOVA model with continuous time handling and rigorous uncertainty quantification.

Findings

Revealed detailed 3D ocean cooling patterns along TC tracks

Identified subsurface warming on the right side of storm tracks

Provided a global-scale characterization of ocean response to TCs

Abstract

Tropical cyclones (TCs), driven by heat exchange between the air and sea, pose a substantial risk to many communities around the world. Accurate characterization of the subsurface ocean thermal response to TC passage is crucial for accurate TC intensity forecasts and for an understanding of the role that TCs play in the global climate system. However, that characterization is complicated by the high-noise ocean environment, correlations inherent in spatio-temporal data, relative scarcity of in situ observations, and the entanglement of the TC-induced signal with seasonal signals. We present a general methodological framework that addresses these difficulties, integrating existing techniques in seasonal mean field estimation, Gaussian process modeling, and nonparametric regression into an ANOVA decomposition model. Importantly, we improve upon past work by properly handling seasonality, providing rigorous uncertainty quantification, and treating time as a continuous variable, rather than producing estimates that are binned in time. This ANOVA model is estimated using in situ subsurface temperature profiles from the Argo fleet of autonomous floats through a multi-step procedure, which (1) characterizes the upper ocean seasonal shift during the TC season; (2) models the variability in the temperature observations; (3) fits a thin plate spline using the variability estimates to account for heteroskedasticity and correlation between the observations. This spline fit reveals the ocean thermal response to TC passage. Through this framework, we obtain new scientific insights into the interaction between TCs and the ocean on a global scale, including a three-dimensional characterization of the near-surface and subsurface cooling along the TC storm track and the mixing-induced subsurface warming on the track's right side.