# Advection surface-flux balance controls the seasonal steric sea level amplitude

Antoine Hochet, William Llovel, Thierry Huck, Florian Sévellec

► To cite this version:

Antoine Hochet, William Llovel, Thierry Huck, Florian Sévellec. Advection surface-flux balance controls the seasonal steric sea level amplitude. Scientific Reports, 2024, 14 (1), pp.10644. 10.1038/s41598-024-61447-y. hal-04586366

HAL Id: hal-04586366

https://hal.univ-brest.fr/hal-04586366

Submitted on 24 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Advection surface-flux balance controls the seasonal steric sea level amplitude

Antoine Hochet, William Llovel, Thierry Huck, Florian Sévellec

► To cite this version:

Antoine Hochet, William Llovel, Thierry Huck, Florian Sévellec. Advection surface-flux balance controls the seasonal steric sea level amplitude. Scientific Reports, 2024, 14 (1), pp.10644. 10.1038/s41598-024-61447-y. hal-04586366

HAL Id: hal-04586366

https://hal.univ-brest.fr/hal-04586366v1

Submitted on 24 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Advection surface-flux balance controls the seasonal steric sea level amplitude

Antoine Hochet<fnref n="1" />✉, William Llovel<fnref n="1" />, Thierry Huck<fnref n="1" /> & Florian Sévellec<fnref n="1" />,²

Along with the mean sea level rise due to climate change, the sea level exhibits natural variations at a large number of different time scales. One of the most important is the one linked with the seasonal cycle. In the Northern Hemisphere winter, the sea level is as much as 20 cm below its summer values in some locations. It is customary to associate these variations with the seasonal cycle of the sea surface net heat flux which drives an upper-ocean thermal expansion creating a positive steric sea level anomaly. Here, using a novel framework based on steric sea level variance budget applied to observations and to the Estimating the Circulation and Climate of the Ocean state estimate, we demonstrate that the steric sea level seasonal cycle amplitude results from a balance between the seasonal sea surface net heat flux and the oceanic advective processes. Moreover, for up to 50% of the ocean surface, surface heat fluxes act to damp the seasonal steric sea level cycle amplitude, which is instead forced by oceanic advection processes. We also show that eddies play an important role in damping the steric sea level seasonal cycle. Our study contributes to a better understanding of the steric sea level mechanisms which is crucial to ensure accurate and reliable climate projections.

**Keywords** Seasonal steric sea level, Physical mechanisms, Surface heat fluxes, Ocean advective processes

Global mean sea level (GMSL) rise is one of the most emblematic and objective consequences of current global warming. The rate of GMSL rise reached 3 mm per year over the 1993–2015 period<fnref n="1" />. However, this trend is modulated by a seasonal signal which reaches a magnitude of 4 mm at the end of boreal summer², based on satellite altimetry data. The main drivers of GMSL rise are global ocean mass increase (i.e., barystatic sea level³) and global ocean density change (global mean steric sea level). The latter is primarily controlled by global ocean warming. Although global mean sea level is highly relevant as a climate index, coastal population are more affected by regional sea level variations. These regional variations are usually much larger than the yearly increase of GMSL. Locally, the sea level seasonal cycle can indeed have amplitudes up to 20 cm in some regions⁴, and are thus comparable to the magnitude of the GMSL rise over the last century⁵,⁶. It is therefore of primary importance to understand their mechanisms, in order to ensure that they are correctly represented in climate models and that future estimates are accurate.

Once corrected for the atmospheric pressure effect, these regional sea level variations are mostly associated with steric sea level (SSL) variations: they owe their existence to the seasonal variations of the oceanic density field⁸. Previous studies based on the scale analysis of the SSL variations equation⁸⁻¹² have found that these variations are mainly driven by the seasonal variations of net heat flux to the ocean. During the Northern and Southern Hemisphere summer months, this anomalous net heat flux (i.e., “warming”) leads to a thermal expansion of the water column inducing a positive anomaly of SSL. Changes in the steric component of the sea level η are assumed to be given at first order by¹⁰:

$$ \partial_t \eta \approx \frac{\alpha Q}{\rho_0 C_p}, \qquad (1) $$

where Q is the anomalous sea surface net heat flux, α the thermal expansion coefficient, ρ₀ the reference density, and Cₚ the specific heat of seawater. The reconstruction of the SSL from the time integration of Eq. (1) is generally found to be well correlated with the true SSL everywhere except at low latitudes¹⁰, where the effect of wind and baroclinic Rossby wave propagation needs to be included¹³⁻¹⁹. These positive correlations are often interpreted as the signature of the leading role of the sea surface seasonal heat flux anomalies in driving the seasonal SSL variations⁴,⁸,¹⁰,¹²,²⁰⁻²⁴. The underlying assumption behind Eq. (1) is that the SSL seasonal cycle is a passive local

<footnote n="1">¹Univ Brest, CNRS, IFREMER, IRD, Laboratoire d'Océanographie Physique et Spatiale (LOPS, UMR 6523), IUEM, Brest, France. ²INRIA, CNRS, ODYSSEY Team-Project, Brest, France. ³email: antoine.hochet@univ-brest.fr</footnote>
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Advection surface-flux balance controls the seasonal steric sea level amplitude

Antoine Hochet<fnref n="1" />✉, William Llovel<fnref n="1" />, Thierry Huck<fnref n="1" /> & Florian Sévellec<fnref n="1" />,²

Along with the mean sea level rise due to climate change, the sea level exhibits natural variations at a large number of different time scales. One of the most important is the one linked with the seasonal cycle. In the Northern Hemisphere winter, the sea level is as much as 20 cm below its summer values in some locations. It is customary to associate these variations with the seasonal cycle of the sea surface net heat flux which drives an upper-ocean thermal expansion creating a positive steric sea level anomaly. Here, using a novel framework based on steric sea level variance budget applied to observations and to the Estimating the Circulation and Climate of the Ocean state estimate, we demonstrate that the steric sea level seasonal cycle amplitude results from a balance between the seasonal sea surface net heat flux and the oceanic advective processes. Moreover, for up to 50% of the ocean surface, surface heat fluxes act to damp the seasonal steric sea level cycle amplitude, which is instead forced by oceanic advection processes. We also show that eddies play an important role in damping the steric sea level seasonal cycle. Our study contributes to a better understanding of the steric sea level mechanisms which is crucial to ensure accurate and reliable climate projections.

**Keywords** Seasonal steric sea level, Physical mechanisms, Surface heat fluxes, Ocean advective processes

Global mean sea level (GMSL) rise is one of the most emblematic and objective consequences of current global warming. The rate of GMSL rise reached 3 mm per year over the 1993–2015 period<fnref n="1" />. However, this trend is modulated by a seasonal signal which reaches a magnitude of 4 mm at the end of boreal summer², based on satellite altimetry data. The main drivers of GMSL rise are global ocean mass increase (i.e., barystatic sea level³) and global ocean density change (global mean steric sea level). The latter is primarily controlled by global ocean warming. Although global mean sea level is highly relevant as a climate index, coastal population are more affected by regional sea level variations. These regional variations are usually much larger than the yearly increase of GMSL. Locally, the sea level seasonal cycle can indeed have amplitudes up to 20 cm in some regions⁴, and are thus comparable to the magnitude of the GMSL rise over the last century⁵,⁶. It is therefore of primary importance to understand their mechanisms, in order to ensure that they are correctly represented in climate models and that future estimates are accurate.

Once corrected for the atmospheric pressure effect, these regional sea level variations are mostly associated with steric sea level (SSL) variations: they owe their existence to the seasonal variations of the oceanic density field⁸. Previous studies based on the scale analysis of the SSL variations equation⁸⁻¹² have found that these variations are mainly driven by the seasonal variations of net heat flux to the ocean. During the Northern and Southern Hemisphere summer months, this anomalous net heat flux (i.e., “warming”) leads to a thermal expansion of the water column inducing a positive anomaly of SSL. Changes in the steric component of the sea level η are assumed to be given at first order by¹⁰:

$$ \partial_t \eta \approx \frac{\alpha Q}{\rho_0 C_p}, \qquad (1) $$

where Q is the anomalous sea surface net heat flux, α the thermal expansion coefficient, ρ₀ the reference density, and Cₚ the specific heat of seawater. The reconstruction of the SSL from the time integration of Eq. (1) is generally found to be well correlated with the true SSL everywhere except at low latitudes¹⁰, where the effect of wind and baroclinic Rossby wave propagation needs to be included¹³⁻¹⁹. These positive correlations are often interpreted as the signature of the leading role of the sea surface seasonal heat flux anomalies in driving the seasonal SSL variations⁴,⁸,¹⁰,¹²,²⁰⁻²⁴. The underlying assumption behind Eq. (1) is that the SSL seasonal cycle is a passive local

<footnote n="1">¹Univ Brest, CNRS, IFREMER, IRD, Laboratoire d'Océanographie Physique et Spatiale (LOPS, UMR 6523), IUEM, Brest, France. ²INRIA, CNRS, ODYSSEY Team-Project, Brest, France. ³email: antoine.hochet@univ-brest.fr</footnote>
<figure><img src="image_3.png" /></figure>

**Figure 1.** Observations reveal that the sea level seasonal cycle cannot be solely explained by surface net heat flux variations. Time average (1993–2014) of the product between the two terms of Eq. (1) and seasonal variations in SSL ($\eta$): (a) left hand side (i.e. $\eta \frac{\partial \eta}{\partial t}$) (in cm$^2$ year$^{-1}$) where $\eta$ is approximated by satellite altimetry$^7$ and (b) right hand side (i.e. $\eta \frac{\alpha Q}{\rho_0 C_p}$) (in cm$^2$ year$^{-1}$) where Q is obtained from the ERA5 reanalysis. The rhs term (a) is different from the rhs term (b), thus demonstrating that important terms are missing in Eq. (1).

<figure><img src="image_4.png" /></figure>

**Figure 2.** Seasonal SSL cycle in ECCOv4r3 1993–2014. (left) Seasonal SSL anomaly (in cm) averaged for January–February–March (a) and July–August–September (b). (right) Log10 of the variance (cm$^2$) of the seasonal SSL anomalies in ECCOv4r3 over the 1993–2014 period (c) and percentage of the total variance explained by the seasonal cycle (d).

response of the ocean to the seasonal cycle in surface net heat flux. Multiplying the SSL time variations equation by $\eta$ yields an equation for the time tendencies of $\eta$ squared i.e. for the amplitude of the seasonal cycle. This equation allows us to determine which of the terms are acting to increase or decrease the amplitude of the cycle. Considering Eq. (1), if $\eta \frac{\alpha Q}{\rho_0 C_p} > 0$, the amplitude increases and if $\eta \frac{\alpha Q}{\rho_0 C_p} < 0$, the amplitude decreases (where the over bar represents the time mean). As Eq. (1) has only two terms, they should remain equal when multiplied by $\eta$. However, we show in Fig. 1, using a combination of satellite altimetry (AVISO$^7$) and surface net heat flux from the ECMWF ERA5 reanalysis$^{25}$ that, when multiplied by $\eta$ and time averaged over the period 1993–2014, the left hand side term of Eq. (1) is not equal to its right hand side term as it should be if Eq. (1) was complete (see “Methods” for more details about this calculation). Moreover, the magnitude of the right hand side term is around 30 cm$^2$ year$^{-1}$ (Fig. 1b) everywhere while the SSL variance magnitude is around 100 cm$^2$ (see Fig. 2c). It follows that without any other term in Eq. (1) a time scale of 3 to 4 years (100 cm$^2$/30 cm$^2$ year$^{-1}$) would be sufficient to completely modify the present SSL seasonal cycle. Additionally, Fig. 1b shows that large regions of the ocean have negative values of $\eta \frac{\alpha Q}{\rho_0 C_p}$ suggesting that the amplitude of the SSL seasonal cycle is, in large parts of the ocean and across all latitudes, damped by the seasonal heat flux, implying that several important terms
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Observations reveal that the sea level seasonal cycle cannot be solely explained by surface net heat flux variations. Time average (1993–2014) of the product between the two terms of Eq. (1) and seasonal variations in SSL ($\eta$): (a) left hand side (i.e. $\eta \frac{\partial \eta}{\partial t}$) (in cm$^2$ year$^{-1}$) where $\eta$ is approximated by satellite altimetry$^7$ and (b) right hand side (i.e. $\eta \frac{\alpha Q}{\rho_0 C_p}$) (in cm$^2$ year$^{-1}$) where Q is obtained from the ERA5 reanalysis. The rhs term (a) is different from the rhs term (b), thus demonstrating that important terms are missing in Eq. (1).</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Seasonal SSL cycle in ECCOv4r3 1993–2014. (left) Seasonal SSL anomaly (in cm) averaged for January–February–March (a) and July–August–September (b). (right) Log10 of the variance (cm$^2$) of the seasonal SSL anomalies in ECCOv4r3 over the 1993–2014 period (c) and percentage of the total variance explained by the seasonal cycle (d).</figcaption></figure>

response of the ocean to the seasonal cycle in surface net heat flux. Multiplying the SSL time variations equation by $\eta$ yields an equation for the time tendencies of $\eta$ squared i.e. for the amplitude of the seasonal cycle. This equation allows us to determine which of the terms are acting to increase or decrease the amplitude of the cycle. Considering Eq. (1), if $\eta \frac{\alpha Q}{\rho_0 C_p} > 0$, the amplitude increases and if $\eta \frac{\alpha Q}{\rho_0 C_p} < 0$, the amplitude decreases (where the over bar represents the time mean). As Eq. (1) has only two terms, they should remain equal when multiplied by $\eta$. However, we show in Fig. 1, using a combination of satellite altimetry (AVISO$^7$) and surface net heat flux from the ECMWF ERA5 reanalysis$^{25}$ that, when multiplied by $\eta$ and time averaged over the period 1993–2014, the left hand side term of Eq. (1) is not equal to its right hand side term as it should be if Eq. (1) was complete (see “Methods” for more details about this calculation). Moreover, the magnitude of the right hand side term is around 30 cm$^2$ year$^{-1}$ (Fig. 1b) everywhere while the SSL variance magnitude is around 100 cm$^2$ (see Fig. 2c). It follows that without any other term in Eq. (1) a time scale of 3 to 4 years (100 cm$^2$/30 cm$^2$ year$^{-1}$) would be sufficient to completely modify the present SSL seasonal cycle. Additionally, Fig. 1b shows that large regions of the ocean have negative values of $\eta \frac{\alpha Q}{\rho_0 C_p}$ suggesting that the amplitude of the SSL seasonal cycle is, in large parts of the ocean and across all latitudes, damped by the seasonal heat flux, implying that several important terms
are missing in Eq. (1). This strongly contrasts with the Sea Surface Temperature (SST) seasonal cycle amplitude which is virtually everywhere forced by the seasonal net heat flux i.e. $\frac{SSTQ}{\rho_0 C_p} > 0$ (see Fig. S1 of the supporting information) as already noticed by several authors$^{17,26}$. Understanding what are the missing term(s) in the equation controlling the amplitude of the SSL seasonal cycle and explaining why the atmospheric net heat flux can dampen the SSL seasonal cycle in some regions are the two main objectives of this work.

While the interannual variability and long-term trend of SSL have been often investigated in the previous years$^{27-30}$, the seasonal cycle of SSL has received much less attention. Hence, here, we apply for the first time a powerful diagnostic to characterize the drivers of the SSL seasonal cycle and describe the dynamics of its sources and sinks. This diagnostic is based on steric sea level variance budget. It has recently been developed and applied to interannual steric sea level variations to understand their mechanisms$^{30}$. Variance budgets are a common tool in physics and have been widely applied to various variables (density, temperature, salinity) in the oceanographic literature in the past$^{31-37}$. This diagnostic is constructed in a similar way to the kinetic energy budget, and is a rigorous way to assess the mechanisms controlling the steric sea level variability.

## Budget of seasonal SSL variance

To investigate the balance controlling the amplitude of the seasonal cycle of regional SSL, we compute the budget of variance over 1993–2014 based on the ECCOv4r3 state estimate$^{38}$. The seasonal cycle of each variable is obtained from the time mean of the 22 years monthly time series (see “Methods”). At mid and high latitudes, in the Northern (Southern) Hemisphere, the seasonal SSL anomaly is negative (positive) (~ 10 cm) in boreal winter while it is positive (negative) in summer (Fig. 2a,b). Close to the equator in each hemisphere, the anomaly can be positive or negative over the two periods of time, depending on the location. Western boundary currents, such as the Kuroshio and the Gulf Stream, also have a strong signature in SSL seasonal cycle. The largest values of the seasonal SSL variance are of the order of $10^2$ cm$^2$ (Fig. 2c) and are found at low latitude in the Eastern Pacific, in the North Indian Ocean and the South Equatorial Current regime of the Indian Ocean and in the northern hemisphere western boundary currents. The seasonal SSL variance represents a large proportion (i.e., above 50%) of the total variance in numerous regions (Fig. 2d). This is particularly true in the Atlantic Ocean across all latitudes and in western boundary currents of the Pacific Ocean as well as in the Arabian Sea and western Bay of Bengal.

The variance budget for the seasonal cycle of the SSL ($\eta$) is composed of four different terms (see section “SSL variance budget” in “Methods”):

$$ \frac{1}{2} \frac{\partial \eta^2}{\partial t} = \text{VAR}_{\text{adv}} + \text{VAR}_{\text{dif}} + \text{VAR}_{\text{flu}}. \quad (2) $$

The overline represents the time average over the ECCOv4r3 period, the left hand side term is the averaged time evolution variance of the seasonal cycle of the SSL. On the right hand side, VAR$_{adv}$ is the effect of the oceanic advective terms, VAR$_{dif}$ the effect of parameterized diffusion (including isoneutral and dianeutral mixing), convective adjustment and background vertical mixing, and VAR$_{flu}$ is the effect of the seasonal cycle of the ocean surface buoyancy fluxes. For instance, a positive VAR$_{flu}$ indicates a local source of SSL variance because it acts to increase the SSL variance (i.e. the SSL amplitude). On the contrary, if it is negative, it acts to decrease the variance: it is a sink of SSL variance. The most striking feature in the budget (Fig. 3) is the strong local compensation between the two dominant terms, the buoyancy forcing term (VAR$_{flu}$) and the advective term (VAR$_{adv}$), both reaching values up to $\pm 30$ cm$^2$ year$^{-1}$. It is shown in the supplementary information file that VAR$_{flu}$ is largely dominated by the net heat flux, in agreement with the literature (Fig. S2). On the other hand, the diffusive term (VAR$_{dif}$) is approximately one order of magnitude smaller (bounded by $\pm$ cm$^2$ year$^{-1}$), and can then be ignored in the overall main balance of the cycle. The term associated with the time tendency of the seasonal cycle SSL variance is also an order of magnitude smaller than the other terms, but is not exactly zero. This is explained by the methodology used to extract the seasonal cycle as explained in section “Methods”. This term remains much smaller than the dominant terms of the budget and can be neglected. Figure 3 demonstrates that the main balance in the seasonal SSL variance budget is between the advective and buoyancy forcing terms. It contrasts strongly with Eq. (1), from which a constant amplitude could only be attained if $\overline{\eta Q} = 0$, which is clearly not the case in Fig. 1b as well as in Fig. 3b. The fact that $\overline{\eta Q} \neq 0$ implies that $\frac{\partial \eta}{\partial t}$ is not exactly in phase with Q. The time lag between these two terms is estimated to be between 2 and 3 weeks depending on the location (supplementary file Fig. S3). Although this time lag may seem small, it represents between 15% and 25% of the time lag required for these two terms to be in quadrature, which is about 13 weeks (i.e. a quarter of an annual period). This relatively small time lag is sufficient to induce the large values of VAR$_{flu}$ found in Figs. 1b and 3b. Moreover, the seasonal cycle of the ocean surface buoyancy fluxes acts as a sink of seasonal SSL variance over 49% of the ocean surface. In these regions, the seasonal cycle of the buoyancy flux damps the seasonal cycle of SSL instead of sustaining it in regions where it is positive. The locations where VAR$_{flu}$ acts as a sink are relatively symmetric with respect to the equator: they are found between 30° N and 60° N and between 30° S and 60° S in every oceans, in the eastern half of the Pacific, Atlantic, and Indian Oceans at low latitudes. Very similar patterns are found using ERA5 and satellite altimetry (Fig. 1b) in the eddy-rich regions, such as the Southern Ocean, which are not resolved by the laminar resolution of the ECCO v4. We find that the variance budget is very similar when all frequencies resolved by ECCO v4 are considered (i.e. when the mean seasonal cycle is not extracted) (see Fig. S4 of the supporting information). This first shows that the SSL variance fluxes linked with the mean seasonal cycle dominate the SSL variance budget at all frequencies and secondly demonstrates that our results do not depend on the methodology used to extract the mean seasonal cycle. The horizontal average of each term, given in the titles of Fig. 3, shows that globally, the surface buoyancy flux is a source (0.9 cm$^2$ year$^{-1}$) that is partially balanced by advection
are missing in Eq. (1). This strongly contrasts with the Sea Surface Temperature (SST) seasonal cycle amplitude which is virtually everywhere forced by the seasonal net heat flux i.e. $\frac{SSTQ}{\rho_0 C_p} > 0$ (see Fig. S1 of the supporting information) as already noticed by several authors$^{17,26}$. Understanding what are the missing term(s) in the equation controlling the amplitude of the SSL seasonal cycle and explaining why the atmospheric net heat flux can dampen the SSL seasonal cycle in some regions are the two main objectives of this work.

While the interannual variability and long-term trend of SSL have been often investigated in the previous years$^{27-30}$, the seasonal cycle of SSL has received much less attention. Hence, here, we apply for the first time a powerful diagnostic to characterize the drivers of the SSL seasonal cycle and describe the dynamics of its sources and sinks. This diagnostic is based on steric sea level variance budget. It has recently been developed and applied to interannual steric sea level variations to understand their mechanisms$^{30}$. Variance budgets are a common tool in physics and have been widely applied to various variables (density, temperature, salinity) in the oceanographic literature in the past$^{31-37}$. This diagnostic is constructed in a similar way to the kinetic energy budget, and is a rigorous way to assess the mechanisms controlling the steric sea level variability.

## Budget of seasonal SSL variance

To investigate the balance controlling the amplitude of the seasonal cycle of regional SSL, we compute the budget of variance over 1993–2014 based on the ECCOv4r3 state estimate$^{38}$. The seasonal cycle of each variable is obtained from the time mean of the 22 years monthly time series (see “Methods”). At mid and high latitudes, in the Northern (Southern) Hemisphere, the seasonal SSL anomaly is negative (positive) (~ 10 cm) in boreal winter while it is positive (negative) in summer (Fig. 2a,b). Close to the equator in each hemisphere, the anomaly can be positive or negative over the two periods of time, depending on the location. Western boundary currents, such as the Kuroshio and the Gulf Stream, also have a strong signature in SSL seasonal cycle. The largest values of the seasonal SSL variance are of the order of $10^2$ cm$^2$ (Fig. 2c) and are found at low latitude in the Eastern Pacific, in the North Indian Ocean and the South Equatorial Current regime of the Indian Ocean and in the northern hemisphere western boundary currents. The seasonal SSL variance represents a large proportion (i.e., above 50%) of the total variance in numerous regions (Fig. 2d). This is particularly true in the Atlantic Ocean across all latitudes and in western boundary currents of the Pacific Ocean as well as in the Arabian Sea and western Bay of Bengal.

The variance budget for the seasonal cycle of the SSL ($\eta$) is composed of four different terms (see section “SSL variance budget” in “Methods”):

$$ \frac{1}{2} \frac{\partial \eta^2}{\partial t} = \text{VAR}_{\text{adv}} + \text{VAR}_{\text{dif}} + \text{VAR}_{\text{flu}}. \qquad (2) $$

The overline represents the time average over the ECCOv4r3 period, the left hand side term is the averaged time evolution variance of the seasonal cycle of the SSL. On the right hand side, VAR$_{adv}$ is the effect of the oceanic advective terms, VAR$_{dif}$ the effect of parameterized diffusion (including isoneutral and dianeutral mixing), convective adjustment and background vertical mixing, and VAR$_{flu}$ is the effect of the seasonal cycle of the ocean surface buoyancy fluxes. For instance, a positive VAR$_{flu}$ indicates a local source of SSL variance because it acts to increase the SSL variance (i.e. the SSL amplitude). On the contrary, if it is negative, it acts to decrease the variance: it is a sink of SSL variance. The most striking feature in the budget (Fig. 3) is the strong local compensation between the two dominant terms, the buoyancy forcing term (VAR$_{flu}$) and the advective term (VAR$_{adv}$), both reaching values up to $\pm 30$ cm$^2$ year$^{-1}$. It is shown in the supplementary information file that VAR$_{flu}$ is largely dominated by the net heat flux, in agreement with the literature (Fig. S2). On the other hand, the diffusive term (VAR$_{dif}$) is approximately one order of magnitude smaller (bounded by $\pm$ cm$^2$ year$^{-1}$), and can then be ignored in the overall main balance of the cycle. The term associated with the time tendency of the seasonal cycle SSL variance is also an order of magnitude smaller than the other terms, but is not exactly zero. This is explained by the methodology used to extract the seasonal cycle as explained in section “Methods”. This term remains much smaller than the dominant terms of the budget and can be neglected. Figure 3 demonstrates that the main balance in the seasonal SSL variance budget is between the advective and buoyancy forcing terms. It contrasts strongly with Eq. (1), from which a constant amplitude could only be attained if $\overline{\eta Q} = 0$, which is clearly not the case in Fig. 1b as well as in Fig. 3b. The fact that $\overline{\eta Q} \neq 0$ implies that $\frac{\partial \eta}{\partial t}$ is not exactly in phase with Q. The time lag between these two terms is estimated to be between 2 and 3 weeks depending on the location (supplementary file Fig. S3). Although this time lag may seem small, it represents between 15% and 25% of the time lag required for these two terms to be in quadrature, which is about 13 weeks (i.e. a quarter of an annual period). This relatively small time lag is sufficient to induce the large values of VAR$_{flu}$ found in Figs. 1b and 3b. Moreover, the seasonal cycle of the ocean surface buoyancy fluxes acts as a sink of seasonal SSL variance over 49% of the ocean surface. In these regions, the seasonal cycle of the buoyancy flux damps the seasonal cycle of SSL instead of sustaining it in regions where it is positive. The locations where VAR$_{flu}$ acts as a sink are relatively symmetric with respect to the equator: they are found between 30° N and 60° N and between 30° S and 60° S in every oceans, in the eastern half of the Pacific, Atlantic, and Indian Oceans at low latitudes. Very similar patterns are found using ERA5 and satellite altimetry (Fig. 1b) in the eddy-rich regions, such as the Southern Ocean, which are not resolved by the laminar resolution of the ECCO v4. We find that the variance budget is very similar when all frequencies resolved by ECCO v4 are considered (i.e. when the mean seasonal cycle is not extracted) (see Fig. S4 of the supporting information). This first shows that the SSL variance fluxes linked with the mean seasonal cycle dominate the SSL variance budget at all frequencies and secondly demonstrates that our results do not depend on the methodology used to extract the mean seasonal cycle. The horizontal average of each term, given in the titles of Fig. 3, shows that globally, the surface buoyancy flux is a source ($0.9$ cm$^2$ year$^{-1}$) that is partially balanced by advection
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Seasonal SSL variance budget. (a) Time mean of the time tendencies of the SSL seasonal variance (in cm² year⁻¹). (b) Contribution of the ocean surface buoyancy fluxes. (c) Contribution of the diffusive terms on the budget. (d) Contribution of the advective term. Positive values (red) indicates that the corresponding term acts to increase the seasonal SSL variance, negative values (blue) that it acts to decrease the variance. The advective term (VAR<sub>adv</sub>) is further decomposed into the advection by (e) the resolved “laminar” velocities (VAR<sub>adv</sub><sup>laminar</sup>) and (f) the eddy-induced velocities from the GM parameterization (VAR<sub>adv</sub><sup>eddy</sup>). The global average of each term is given in the corresponding titles.</figcaption></figure>

($-0.5 \text{ cm}^2 \text{ year}^{-1}$) over the period 1993–2014 (global mean diffusion is $-0.2 \text{ cm}^2 \text{ year}^{-1}$, being a net sink, and the variance time tendency is $0.2 \text{ cm}^2 \text{ year}^{-1}$).

To investigate the large-scale mechanisms of seasonal SSL variations, we also compute the variance equation for the hemispheric average of the SSL as follows (see “Methods”):

$$ \frac{1}{2} \frac{\partial \langle \eta \rangle^2}{\partial t} = \mathrm{VAR}_{\mathrm{adv}}^{\mathrm{H}} + \mathrm{VAR}_{\mathrm{dif}}^{\mathrm{H}} + \mathrm{VAR}_{\mathrm{flu}}^{\mathrm{H}}, \quad (3) $$

where $\langle \eta \rangle$ is the hemispheric average of the seasonal SSL, VAR<sub>adv</sub><sup>H</sup>, VAR<sub>dif</sub><sup>H</sup>, VAR<sub>flu</sub><sup>H</sup> respectively the contribution of oceanic advection, diffusion and buoyancy flux. The hemispheric averages are computed over several regions bounded by different latitudes. For the northern hemisphere the averaging region is defined by all locations with latitudes north of a given latitude, which varies from 0° N to 70° N. Similarly, for the southern hemisphere, the averaging region is defined by all locations with latitudes south of a given latitude, which varies from 70° S to 0° S. The results show that the large-scale seasonal SSL variations obey a similar balance as for the local variance budget, with the oceanic advective terms and the buoyancy flux term being the two main contributors to the budget in both hemispheres (Fig. 4a,b). The balance in both hemispheres is nearly symmetric around 8° N, although the amplitudes of the averaged SSL variance budget terms are about twice as large in the Northern Hemisphere (NH) as in the Southern Hemisphere (SH). When the region encompasses the entire hemisphere (corresponding to latitude 0° N for the northern hemisphere and 0° S for the southern hemisphere in Fig. 4), the main source is attributed to the buoyancy flux term, which is mainly due to net heat flux ($3 \text{ cm}^2 \text{ year}^{-1}$ in the
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Seasonal SSL variance budget. (a) Time mean of the time tendencies of the SSL seasonal variance (in cm² year⁻¹). (b) Contribution of the ocean surface buoyancy fluxes. (c) Contribution of the diffusive terms on the budget. (d) Contribution of the advective term. Positive values (red) indicates that the corresponding term acts to increase the seasonal SSL variance, negative values (blue) that it acts to decrease the variance. The advective term (VAR<sub>adv</sub>) is further decomposed into the advection by (e) the resolved “laminar” velocities (VAR<sub>adv</sub><sup>laminar</sup>) and (f) the eddy-induced velocities from the GM parameterization (VAR<sub>adv</sub><sup>eddy</sup>). The global average of each term is given in the corresponding titles.</figcaption></figure>

($-0.5 \text{ cm}^2 \text{ year}^{-1}$) over the period 1993–2014 (global mean diffusion is $-0.2 \text{ cm}^2 \text{ year}^{-1}$, being a net sink, and the variance time tendency is $0.2 \text{ cm}^2 \text{ year}^{-1}$).

To investigate the large-scale mechanisms of seasonal SSL variations, we also compute the variance equation for the hemispheric average of the SSL as follows (see “Methods”):

$$ \frac{1}{2} \frac{\partial \langle \eta \rangle^2}{\partial t} = \text{VAR}_{\text{adv}}^{\text{H}} + \text{VAR}_{\text{dif}}^{\text{H}} + \text{VAR}_{\text{flu}}^{\text{H}}, \quad (3) $$

where $\langle \eta \rangle$ is the hemispheric average of the seasonal SSL, VAR<sub>adv</sub><sup>H</sup>, VAR<sub>dif</sub><sup>H</sup>, VAR<sub>flu</sub><sup>H</sup> respectively the contribution of oceanic advection, diffusion and buoyancy flux. The hemispheric averages are computed over several regions bounded by different latitudes. For the northern hemisphere the averaging region is defined by all locations with latitudes north of a given latitude, which varies from 0° N to 70° N. Similarly, for the southern hemisphere, the averaging region is defined by all locations with latitudes south of a given latitude, which varies from 70° S to 0° S. The results show that the large-scale seasonal SSL variations obey a similar balance as for the local variance budget, with the oceanic advective terms and the buoyancy flux term being the two main contributors to the budget in both hemispheres (Fig. 4a,b). The balance in both hemispheres is nearly symmetric around 8° N, although the amplitudes of the averaged SSL variance budget terms are about twice as large in the Northern Hemisphere (NH) as in the Southern Hemisphere (SH). When the region encompasses the entire hemisphere (corresponding to latitude 0° N for the northern hemisphere and 0° S for the southern hemisphere in Fig. 4), the main source is attributed to the buoyancy flux term, which is mainly due to net heat flux ($3 \text{ cm}^2 \text{ year}^{-1}$ in the
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Variance budget for the spatially averaged Seasonal SSL. (a) Terms of the variance budget (in cm² year⁻¹) for the spatially averaged SSL over Southern Hemisphere regions with latitudes south of the latitude given by the abscissa value. (b) same but for Northern Hemisphere regions with latitudes north of the abscissa value. The Contributions of the advective term (blue line), time mean of the time tendencies of the SSL seasonal variance (orange line), diffusive terms (green line) and of the ocean surface buoyancy fluxes (red line) are shown. (c) Variance of the regionally averaged SSL seasonal cycle for latitudes south of the abscissa value (in cm²), (d) Same but for the northern hemisphere (i.e. latitudes north of the abscissa value)..</figcaption></figure>

SH, 5 cm² year⁻¹ in the NH), while the main sink is associated with the advective term (−3 cm² year⁻¹ in the SH, −5 cm² year⁻¹ in the NH). Although the amplitude of the seasonal variations of the hemispheric averages of the advective term (⟨adv⟩, see Eq. (14) in “Methods”), is almost 3 times weaker than that of the hemispheric averages of the surface fluxes (⟨flu⟩), ⟨adv⟩ is almost in phase with ⟨η⟩ and can balance the effect of ⟨flu⟩, which is almost in quadrature with ⟨η⟩ (see Fig. S5 of the supplementary file). For a region bounded to the North by 30° S for the Southern hemisphere and to the South by 30° N for the Northern hemisphere, the advection term is the main source of SSL seasonal variance (1 cm² year⁻¹ in the SH, 2.5 cm² year⁻¹ in the NH) and is balanced by the surface buoyancy flux term. The variance of the regionally averaged SSL is between 3 and 8 cm² in the NH and between 1 and 4 cm² in the SH (Fig. 4c,d) As for the local variance budget, calculating the ratio of the regionally averaged SSL variance to the budget term values gives the time scale that would be sufficient to completely change the characteristics of the current regional SSL seasonal cycle. This time scale is less than one year for the entire NH and SH regions, and about 3 years with a southern/northern boundary located at 30° S and 30° N. This demonstrates the importance of the balance between the advective terms and the surface buoyancy fluxes in determining the amplitude of the regional SSL seasonal cycle.

## The vertical structure of the density anomaly controls the sign of VAR<sub>flu</sub>

As shown in “Methods”, VAR<sub>flu</sub> is proportional to the product of the vertically-averaged density seasonal anomaly (ρ<sub>BT</sub>) and the vertically-averaged buoyancy anomaly which is mostly due to the net heat flux (Q<sub>BT</sub>, see Fig. S2):

$$ VAR_{\text{flu}} \propto \overline{\rho_{\text{BT}} Q_{\text{BT}}}, \qquad (4) $$

this formula demonstrates that the sign of VAR<sub>flu</sub> depends on the correlation between these two fields. We show (Fig. 5) the density anomaly in the first 300 m as a function of time and depth, averaged North of 20° N in regions where VAR<sub>flu</sub> > 0 (left column) and in regions where VAR<sub>flu</sub> < 0 (right column). We average North of 20° N to ensure that the density anomalies have the same phase everywhere, but the same results hold for the Southern Hemisphere. The propagation of the density anomaly to deeper depths is not instantaneous during the seasonal cycle (top panels Fig. 5) which implies that the vertically-averaged density anomalies are not necessarily in phase with the surface density anomaly. On the contrary, the atmospheric buoyancy flux is confined to the first 50 m and has a coherent phase with depth. Thus the sign of VAR<sub>flu</sub> is controlled by the vertical structure of the density anomalies. The surface density is almost everywhere positively correlated with Q<sub>BT</sub> (supplementary Fig. S6). However, regions where the vertically-averaged density anomaly is influenced by shifted deeper density anomalies are associated with negative values of VAR<sub>flu</sub> (right column Fig. 5). The comparison of the time evolution of the vertically-averaged density (orange line, bottom panels Fig. 5) with the vertical integral of the seasonal buoyancy flux (blue line, Fig. 5) reveals a different behavior between regions where VAR<sub>flu</sub> > 0 and regions where VAR<sub>flu</sub> < 0. In regions where VAR<sub>flu</sub> is a source (positive), the vertically-averaged density anomaly follows the vertical integral of the seasonal cycle of buoyancy flux (blue line) with a delay of approximately two
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Variance budget for the spatially averaged Seasonal SSL. (a) Terms of the variance budget (in cm² year⁻¹) for the spatially averaged SSL over Southern Hemisphere regions with latitudes south of the latitude given by the abscissa value. (b) same but for Northern Hemisphere regions with latitudes north of the abscissa value. The Contributions of the advective term (blue line), time mean of the time tendencies of the SSL seasonal variance (orange line), diffusive terms (green line) and of the ocean surface buoyancy fluxes (red line) are shown. (c) Variance of the regionally averaged SSL seasonal cycle for latitudes south of the abscissa value (in cm²), (d) Same but for the northern hemisphere (i.e. latitudes north of the abscissa value)..</figcaption></figure>

SH, 5 cm² year⁻¹ in the NH), while the main sink is associated with the advective term (−3 cm² year⁻¹ in the SH, −5 cm² year⁻¹ in the NH). Although the amplitude of the seasonal variations of the hemispheric averages of the advective term (⟨adv⟩, see Eq. (14) in “Methods”), is almost 3 times weaker than that of the hemispheric averages of the surface fluxes (⟨flu⟩), ⟨adv⟩ is almost in phase with ⟨η⟩ and can balance the effect of ⟨flu⟩, which is almost in quadrature with ⟨η⟩ (see Fig. S5 of the supplementary file). For a region bounded to the North by 30° S for the Southern hemisphere and to the South by 30° N for the Northern hemisphere, the advection term is the main source of SSL seasonal variance (1 cm² year⁻¹ in the SH, 2.5 cm² year⁻¹ in the NH) and is balanced by the surface buoyancy flux term. The variance of the regionally averaged SSL is between 3 and 8 cm² in the NH and between 1 and 4 cm² in the SH (Fig. 4c,d) As for the local variance budget, calculating the ratio of the regionally averaged SSL variance to the budget term values gives the time scale that would be sufficient to completely change the characteristics of the current regional SSL seasonal cycle. This time scale is less than one year for the entire NH and SH regions, and about 3 years with a southern/northern boundary located at 30° S and 30° N. This demonstrates the importance of the balance between the advective terms and the surface buoyancy fluxes in determining the amplitude of the regional SSL seasonal cycle.

## The vertical structure of the density anomaly controls the sign of VARflu

As shown in “Methods”, VARflu is proportional to the product of the vertically-averaged density seasonal anomaly (ρBT) and the vertically-averaged buoyancy anomaly which is mostly due to the net heat flux (QBT, see Fig. S2):

$$ VAR_{flu} \propto \overline{\rho_{BT} Q_{BT}}, \qquad (4) $$

this formula demonstrates that the sign of VARflu depends on the correlation between these two fields. We show (Fig. 5) the density anomaly in the first 300 m as a function of time and depth, averaged North of 20° N in regions where VARflu > 0 (left column) and in regions where VARflu < 0 (right column). We average North of 20° N to ensure that the density anomalies have the same phase everywhere, but the same results hold for the Southern Hemisphere. The propagation of the density anomaly to deeper depths is not instantaneous during the seasonal cycle (top panels Fig. 5) which implies that the vertically-averaged density anomalies are not necessarily in phase with the surface density anomaly. On the contrary, the atmospheric buoyancy flux is confined to the first 50 m and has a coherent phase with depth. Thus the sign of VARflu is controlled by the vertical structure of the density anomalies. The surface density is almost everywhere positively correlated with QBT (supplementary Fig. S6). However, regions where the vertically-averaged density anomaly is influenced by shifted deeper density anomalies are associated with negative values of VARflu (right column Fig. 5). The comparison of the time evolution of the vertically-averaged density (orange line, bottom panels Fig. 5) with the vertical integral of the seasonal buoyancy flux (blue line, Fig. 5) reveals a different behavior between regions where VARflu > 0 and regions where VARflu < 0. In regions where VARflu is a source (positive), the vertically-averaged density anomaly follows the vertical integral of the seasonal cycle of buoyancy flux (blue line) with a delay of approximately two
<figure><img src="image_7.png" /><figcaption>Figure 5. The vertical structure of density anomalies controls the sign of the surface buoyancy flux contribution to SSL variance (VAR<sub>flu</sub>). Top panels: density anomaly (in kg m<sup>-3</sup>) as a function of depth and time (month). Bottom panels: time variation of the vertically-averaged density (orange line), vertically-averaged net heat flux (blue line) and product of the two terms, i.e., ρ<sub>BT</sub> Q<sub>BT</sub> (black line, α VAR<sub>flu</sub>). The surface density (orange dashed) is also shown. Each time series is normalized, red shading indicates where VAR<sub>flu</sub> is positive and blue shading where it is negative. The left and right columns show respectively the Northern Hemisphere values averaged at locations where VAR<sub>flu</sub> > 0 and where VAR<sub>flu</sub> < 0.</figcaption></figure>

months. The time series of the product between ρ<sub>BT</sub> and Q<sub>BT</sub> is then strongly positive in boreal winter and summer (months 11 to 3 and months 6 to 9) and weakly negative for the remaining periods, resulting in an overall positive time average (VAR<sub>flu</sub> > 0). Moreover, the vertically-averaged density is almost indistinguishable from the surface density (orange dashed line). In regions where VAR<sub>flu</sub> is a sink (negative), the vertically-averaged density anomaly is shifted away from the vertical integral of the seasonal buoyancy flux because of the influence of deeper density anomalies. The delay between the vertical integral of the seasonal cycle of buoyancy flux and the vertically-averaged density is then larger (~4 months) than for the previous case and induces large negative values of the product between ρ<sub>BT</sub> and Q<sub>BT</sub> in boreal spring and autumn (months 3 to 6 and months 9 to 12). Unlike the previous case, the vertically-averaged density is also shifted with respect to the surface density. It is therefore the influence of deeper density anomalies that explains why VAR<sub>flu</sub> is a sink (negative) over large areas of the ocean.

## Eddies dampen the seasonal SSL cycle

To understand the role of eddies in VAR<sub>adv</sub> (Fig. 3d), this term is decomposed into two terms (see section “Methods”):

$$VAR_{adv} = VAR_{adv}^{\text{laminar}} + VAR_{adv}^{\text{eddy}} \qquad (5)$$

VAR<sub>adv</sub><sup>laminar</sup> is linked to the advection of density by the resolved “laminar” currents, whereas VAR<sub>adv</sub><sup>eddy</sup> is linked to the advection by the eddy-induced velocities as parameterized in Gent and McWilliams<sup>39</sup>. VAR<sub>adv</sub><sup>eddy</sup> is important almost everywhere (Fig. 3e) and can be a source or a sink of seasonal SSL variability depending on the region. Its pattern is similar to that of VAR<sub>adv</sub> (Fig. 3d), except around the equator and in western intensified boundary currents. When horizontally averaged, it is the largest source for the seasonal cycle (5.2 cm<sup>2</sup> year<sup>-1</sup>). At the global scale, this source is almost exactly compensated by the sink made by the term linked with eddies VAR<sub>adv</sub><sup>eddy</sup> (~5.7 cm<sup>2</sup> year<sup>-1</sup>). The Gent & McWilliams parameterization<sup>39</sup> mimics the fact that isopycnals slopes are baroclinically unstable and eddy kinetic energy is created through the release of potential energy. This eddy parameterization leads to a damping of seasonal SSL variance. Locally, VAR<sub>adv</sub><sup>eddy</sup> (Fig. 3f) is strongly negative (~-30 cm<sup>2</sup> year<sup>-1</sup>) in the equatorial regions of the Pacific, Atlantic, and Indian Oceans. There, it significantly compensates the laminar advection term. In line with our results, but for a different timescale, it has recently
<figure><img src="image_7.png" /><figcaption>Figure 5. The vertical structure of density anomalies controls the sign of the surface buoyancy flux contribution to SSL variance (VAR<sub>flu</sub>). Top panels: density anomaly (in kg m<sup>-3</sup>) as a function of depth and time (month). Bottom panels: time variation of the vertically-averaged density (orange line), vertically-averaged net heat flux (blue line) and product of the two terms, i.e., ρ<sub>BT</sub> Q<sub>BT</sub> (black line, α VAR<sub>flu</sub>). The surface density (orange dashed) is also shown. Each time series is normalized, red shading indicates where VAR<sub>flu</sub> is positive and blue shading where it is negative. The left and right columns show respectively the Northern Hemisphere values averaged at locations where VAR<sub>flu</sub> > 0 and where VAR<sub>flu</sub> < 0.</figcaption></figure>

months. The time series of the product between ρ<sub>BT</sub> and Q<sub>BT</sub> is then strongly positive in boreal winter and summer (months 11 to 3 and months 6 to 9) and weakly negative for the remaining periods, resulting in an overall positive time average (VAR<sub>flu</sub> > 0). Moreover, the vertically-averaged density is almost indistinguishable from the surface density (orange dashed line). In regions where VAR<sub>flu</sub> is a sink (negative), the vertically-averaged density anomaly is shifted away from the vertical integral of the seasonal buoyancy flux because of the influence of deeper density anomalies. The delay between the vertical integral of the seasonal cycle of buoyancy flux and the vertically-averaged density is then larger (~4 months) than for the previous case and induces large negative values of the product between ρ<sub>BT</sub> and Q<sub>BT</sub> in boreal spring and autumn (months 3 to 6 and months 9 to 12). Unlike the previous case, the vertically-averaged density is also shifted with respect to the surface density. It is therefore the influence of deeper density anomalies that explains why VAR<sub>flu</sub> is a sink (negative) over large areas of the ocean.

## Eddies dampen the seasonal SSL cycle

To understand the role of eddies in VAR<sub>adv</sub> (Fig. 3d), this term is decomposed into two terms (see section “Methods”):

$$VAR_{adv} = VAR_{adv}^{\text{laminar}} + VAR_{adv}^{\text{eddy}} \qquad (5)$$

VAR<sub>adv</sub><sup>laminar</sup> is linked to the advection of density by the resolved “laminar” currents, whereas VAR<sub>adv</sub><sup>eddy</sup> is linked to the advection by the eddy-induced velocities as parameterized in Gent and McWilliams<sup>39</sup>. VAR<sub>adv</sub><sup>eddy</sup> is important almost everywhere (Fig. 3e) and can be a source or a sink of seasonal SSL variability depending on the region. Its pattern is similar to that of VAR<sub>adv</sub> (Fig. 3d), except around the equator and in western intensified boundary currents. When horizontally averaged, it is the largest source for the seasonal cycle (5.2 cm<sup>2</sup> year<sup>-1</sup>). At the global scale, this source is almost exactly compensated by the sink made by the term linked with eddies VAR<sub>adv</sub><sup>eddy</sup> (~5.7 cm<sup>2</sup> year<sup>-1</sup>). The Gent & McWilliams parameterization<sup>39</sup> mimics the fact that isopycnals slopes are baroclinically unstable and eddy kinetic energy is created through the release of potential energy. This eddy parameterization leads to a damping of seasonal SSL variance. Locally, VAR<sub>adv</sub><sup>eddy</sup> (Fig. 3f) is strongly negative (~-30 cm<sup>2</sup> year<sup>-1</sup>) in the equatorial regions of the Pacific, Atlantic, and Indian Oceans. There, it significantly compensates the laminar advection term. In line with our results, but for a different timescale, it has recently
been shown that, in the Equatorial Pacific, ENSO variability is inhibited by mesoscale eddies⁴⁰. It also has large values in the Gulf Stream and Kuroshio regions where it can either be a source or a sink of SSL variance.

## Discussion

Based on a new variance budget framework, we have demonstrated in this article that the amplitude of the SSL seasonal cycle is controlled both at local scale and at the hemispheric scale by a balance between the oceanic advection and surface buoyancy forcing terms. The diffusive and variance tendency terms are an order of magnitude smaller. At mid latitudes and in the eastern parts of low latitudes regions, the seasonal cycle of the buoyancy fluxes even acts to damp the amplitude of the seasonal cycle of SSL instead of sustaining it. In these regions, the main source of SSL seasonal variability is associated with oceanic advective terms. Whether the surface buoyancy fluxes are a source or a sink of SSL cycle variance depends on the vertical structure of the density anomalies. Buoyancy fluxes are a source of variance in regions where the vertically-averaged density is controlled by surface values of density. On the contrary, buoyancy fluxes become a sink of variance in regions where density anomalies are controlled by sub-surface density anomalies. The horizontal average of the local variance budget shows that the seasonal ocean surface buoyancy flux is a source for the seasonal cycle of SSL (+ cm² year⁻¹), while advective terms (and diffusive terms) are a sink (respectively -0.5 cm² year⁻¹ and -0.2 cm² year⁻¹). We also show for the first time that eddies are a major sink of the seasonal cycle of SSL close to the equator. The horizontal resolution of the ECCOv4 reanalysis (1° in average) does not, however, explicitly resolve eddies and their effect is only parameterized. Further studies are thus needed to better understand their exact contribution to the seasonal cycle of SSL, in particular for western boundary currents and the Southern ocean where the turbulent field is known to imprint itself strongly on the SSL⁴¹. In higher resolution models, global or regional, a similar methodology could be used to study the seasonal cycle mechanisms at different space scales and determine whether different mechanisms are at play. In this work we have decomposed the oceanic advective terms into two parts: resolved and eddy induced advection. A previous study of the mechanisms of interannual variations in steric sea level has shown that the advective term can be decomposed into several terms each associated with different physical mechanisms³⁰. A perspective is thus to apply the same framework to the study of the seasonal variations in SSL. Although we have shown that ECCOv4 is able to reproduce most of the patterns associated with the buoyancy flux term in the variance equation (Figs. 1a and 3b) giving us confidence in our results, a limitation of this work is the use of a single model which may be associated with substantial error model. Future work should therefore focus on reproducing our results in different numerical models. Another limitation of our study is the focus on the steric sea level. Steric sea level is the major component of sea level variations at low and mid latitudes but seasonal variations in manometric sea level are important at high latitudes and in shallow water regions. Therefore, future studies should also investigate the mechanisms of this component. One important implication of our results is that the amplitude of the sea level seasonal cycle is potentially much more sensitive to anthropogenic climate change than previously assumed. It will indeed be sensitive to changes in the advective and net heat flux seasonal cycle (whether as a source or a sink), as well as to changes in the oceanic mean circulation and eddy field. Numerous studies have demonstrated that important changes in oceanic mean circulation are underway⁴². Studying how the SSL seasonal cycle is modified by anthropogenic climate change is thus a subject of primary importance and we hope that the methodology developed in this work will help to shed light on this matter.

## Methods

### Assessment of the SSL variance budget from Eq. (1) in observations

In this subsection, we describe the methodology used to evaluate the SSL variance equation stemming from Eq. (1) in observations and to obtain Fig. 1. We use the sea level anomaly (SLA) observed by satellite altimetry (AVISO⁷) to approximate the seasonal steric sea level. It has been shown⁸ that seasonal regional variations of sea level are mostly due to the steric sea level and we also show, using the ECCO v4r3 state estimate (more details about this state estimate are given in the following subsection), that this approximation holds almost everywhere at low and mid latitudes except in some semi-enclosed seas (see Fig. S7 of the supplementary file). The AVISO product used has a 1/4° horizontal resolution and we select the period 1993–2014 to be consistent with the analysis performed in the remainder of this article. Following a common practice in oceanography e.g.²⁷⁻²⁹, the mean SLA seasonal cycle is derived by first removing the 1993–2014 trend and then by computing the time mean for each individual month resulting in 12 monthly values. The time tendencies of $\eta$ are obtained by differentiating the daily values of $\eta$ at the start and end of each month divided by the number of days. The mean seasonal cycle of $\eta$ time tendencies is then computed from its monthly time series. The time average of the product between $\eta$ and $\partial_t\eta$ gives Fig. 1a. To obtain Fig. 1b, the net heat flux is first computed from the sum of the net short wave, net long wave, latent and sensible heat fluxes, given by ERA5²⁵ on a 1/4° grid. The seasonal cycle of Q is extracted following the same procedure as for $\eta$. To compute the term $\eta \frac{\alpha Q}{\rho_0 C_p}$ shown in Fig. 1, we use the following coefficients: $\rho_0 = 1028 \text{ kg m}^{-3}$ and $C_p = 4000 \text{ J kg}^{-1} \text{ K}^{-1}$. The thermal expansion coefficient $\alpha$ is computed from the ECCOv4r3 surface temperature and surface salinity, then time averaged over 1993–2014 as well as zonally averaged so that the resulting coefficient is only a function of latitude. We have also checked (Fig. S8 of the supplementary file) that a different choice of net heat flux database (OAFlux⁴³) leads to very similar results.

Note that the methodology used in this article and in many different studies to extract the mean seasonal cycle implies that the product of $\eta$ and $\partial_t\eta$ is not exactly zero. The mean seasonal cycle of the time tendencies of $\eta$ is indeed not exactly equal to the time tendencies of the mean seasonal cycle of $\eta$. If the amplitude of the seasonal cycle increases or decreases over the 22 years of the ECCOv4 period, then the product $\eta\partial_t\eta$ will reflect the mean rate of this change. Additionally, the various approximations used by AVISO to construct daily SLA on a regular 2D grid from an ensemble of inter-calibrated altimeter missions and the use of SLA to approximate SSL seasonal variations also contribute to the fact that $\eta \frac{\partial\eta}{\partial t}$ is not exactly zero. However, Fig. 1a shows that this
been shown that, in the Equatorial Pacific, ENSO variability is inhibited by mesoscale eddies⁴⁰. It also has large values in the Gulf Stream and Kuroshio regions where it can either be a source or a sink of SSL variance.

## Discussion

Based on a new variance budget framework, we have demonstrated in this article that the amplitude of the SSL seasonal cycle is controlled both at local scale and at the hemispheric scale by a balance between the oceanic advection and surface buoyancy forcing terms. The diffusive and variance tendency terms are an order of magnitude smaller. At mid latitudes and in the eastern parts of low latitudes regions, the seasonal cycle of the buoyancy fluxes even acts to damp the amplitude of the seasonal cycle of SSL instead of sustaining it. In these regions, the main source of SSL seasonal variability is associated with oceanic advective terms. Whether the surface buoyancy fluxes are a source or a sink of SSL cycle variance depends on the vertical structure of the density anomalies. Buoyancy fluxes are a source of variance in regions where the vertically-averaged density is controlled by surface values of density. On the contrary, buoyancy fluxes become a sink of variance in regions where density anomalies are controlled by sub-surface density anomalies. The horizontal average of the local variance budget shows that the seasonal ocean surface buoyancy flux is a source for the seasonal cycle of SSL (+ cm² year⁻¹), while advective terms (and diffusive terms) are a sink (respectively -0.5 cm² year⁻¹ and -0.2 cm² year⁻¹). We also show for the first time that eddies are a major sink of the seasonal cycle of SSL close to the equator. The horizontal resolution of the ECCOv4 reanalysis (1° in average) does not, however, explicitly resolve eddies and their effect is only parameterized. Further studies are thus needed to better understand their exact contribution to the seasonal cycle of SSL, in particular for western boundary currents and the Southern ocean where the turbulent field is known to imprint itself strongly on the SSL⁴¹. In higher resolution models, global or regional, a similar methodology could be used to study the seasonal cycle mechanisms at different space scales and determine whether different mechanisms are at play. In this work we have decomposed the oceanic advective terms into two parts: resolved and eddy induced advection. A previous study of the mechanisms of interannual variations in steric sea level has shown that the advective term can be decomposed into several terms each associated with different physical mechanisms³⁰. A perspective is thus to apply the same framework to the study of the seasonal variations in SSL. Although we have shown that ECCOv4 is able to reproduce most of the patterns associated with the buoyancy flux term in the variance equation (Figs. 1a and 3b) giving us confidence in our results, a limitation of this work is the use of a single model which may be associated with substantial error model. Future work should therefore focus on reproducing our results in different numerical models. Another limitation of our study is the focus on the steric sea level. Steric sea level is the major component of sea level variations at low and mid latitudes but seasonal variations in manometric sea level are important at high latitudes and in shallow water regions. Therefore, future studies should also investigate the mechanisms of this component. One important implication of our results is that the amplitude of the sea level seasonal cycle is potentially much more sensitive to anthropogenic climate change than previously assumed. It will indeed be sensitive to changes in the advective and net heat flux seasonal cycle (whether as a source or a sink), as well as to changes in the oceanic mean circulation and eddy field. Numerous studies have demonstrated that important changes in oceanic mean circulation are underway⁴². Studying how the SSL seasonal cycle is modified by anthropogenic climate change is thus a subject of primary importance and we hope that the methodology developed in this work will help to shed light on this matter.

## Methods

### Assessment of the SSL variance budget from Eq. (1) in observations

In this subsection, we describe the methodology used to evaluate the SSL variance equation stemming from Eq. (1) in observations and to obtain Fig. 1. We use the sea level anomaly (SLA) observed by satellite altimetry (AVISO⁷) to approximate the seasonal steric sea level. It has been shown⁸ that seasonal regional variations of sea level are mostly due to the steric sea level and we also show, using the ECCO v4r3 state estimate (more details about this state estimate are given in the following subsection), that this approximation holds almost everywhere at low and mid latitudes except in some semi-enclosed seas (see Fig. S7 of the supplementary file). The AVISO product used has a 1/4° horizontal resolution and we select the period 1993–2014 to be consistent with the analysis performed in the remainder of this article. Following a common practice in oceanography e.g.²⁷⁻²⁹, the mean SLA seasonal cycle is derived by first removing the 1993–2014 trend and then by computing the time mean for each individual month resulting in 12 monthly values. The time tendencies of $\eta$ are obtained by differentiating the daily values of $\eta$ at the start and end of each month divided by the number of days. The mean seasonal cycle of $\eta$ time tendencies is then computed from its monthly time series. The time average of the product between $\eta$ and $\partial_t\eta$ gives Fig. 1a. To obtain Fig. 1b, the net heat flux is first computed from the sum of the net short wave, net long wave, latent and sensible heat fluxes, given by ERA5²⁵ on a 1/4° grid. The seasonal cycle of Q is extracted following the same procedure as for $\eta$. To compute the term $\eta \frac{\alpha Q}{\rho_0 C_p}$ shown in Fig. 1, we use the following coefficients: $\rho_0 = 1028 \text{ kg m}^{-3}$ and $C_p = 4000 \text{ J kg}^{-1} \text{ K}^{-1}$. The thermal expansion coefficient $\alpha$ is computed from the ECCOv4r3 surface temperature and surface salinity, then time averaged over 1993–2014 as well as zonally averaged so that the resulting coefficient is only a function of latitude. We have also checked (Fig. S8 of the supplementary file) that a different choice of net heat flux database (OAFlux⁴³) leads to very similar results.

Note that the methodology used in this article and in many different studies to extract the mean seasonal cycle implies that the product of $\eta$ and $\partial_t\eta$ is not exactly zero. The mean seasonal cycle of the time tendencies of $\eta$ is indeed not exactly equal to the time tendencies of the mean seasonal cycle of $\eta$. If the amplitude of the seasonal cycle increases or decreases over the 22 years of the ECCOv4 period, then the product $\eta\partial_t\eta$ will reflect the mean rate of this change. Additionally, the various approximations used by AVISO to construct daily SLA on a regular 2D grid from an ensemble of inter-calibrated altimeter missions and the use of SLA to approximate SSL seasonal variations also contribute to the fact that $\eta \frac{\partial\eta}{\partial t}$ is not exactly zero. However, Fig. 1a shows that this
term has amplitudes and patterns that are very different from the term linked with the net heat flux (Fig. 1a) and it therefore demonstrates that Eq. (1) is not complete.

## The ECCOv4r3 dataset

Seasonal SSL variance budgets are computed using the ECCOv4r3 state estimate which covers the period 1992–2017. This state estimate is the output of the Massachusetts Institute of Technology general circulation model (MITgcm) assimilating available observations for the period 1992 to 2017³⁸. The advantage of ECCOv4 is that it satisfies the equation of motion and conservation laws hence making it possible to compute tracer budgets. The solution used in this article is computed on the LLC90 grid which has an average horizontal resolution of 1° and 50 vertical levels. Outputs of the model consist of one month average and the closed budget can be obtained for the 1993–2014 period. Thus we compute the SSL budget over this 1993–2014 period. Snapshots at the start and end of each month are also provided by ECCO in order to compute the tracer time tendencies required to close the budgets. ECCO has already been used in several past studies to compute budget of steric sea level variations²⁷⁻³⁰.

For each variable from the monthly model outputs, at each gridpoint, the seasonal cycle is obtained by first removing the linear trend over the 1993–2014 period, and then by computing the 12 monthly values as average of the respective monthly anomalies over the 22 years. The SSL seasonal cycle studied here is thus an average of the seasonal cycle over the 22 years of the 1993–2014 period.

## Seasonal SSL variance budget

The seasonal SSL anomaly $\eta$ is expressed as the vertical integral of the seasonal density anomaly $\rho$ as follows:

$$ \eta = -\frac{1}{\rho_0} \int_{-H}^{0} \rho \, dz, \qquad (6) $$

where $\rho_0 = 1029 \text{ kg m}^{-3}$. The evolution equation for $\eta$ is then simply obtained as the time derivative of (6):

$$ \frac{\partial \eta}{\partial t} = -\frac{1}{\rho_0} \int_{-H}^{0} \frac{\partial \rho}{\partial t} dz. \qquad (7) $$

Following³⁰, the time evolution of $\rho$ can be decomposed into the time evolution of the potential temperature $\theta$ and salinity $S$ as:

$$ \frac{\partial \rho}{\partial t} = -\rho_0 \alpha \frac{\partial \theta}{\partial t} + \rho_0 \beta \frac{\partial S}{\partial t}, \qquad (8) $$

where $\alpha = -\frac{1}{\rho_0} \frac{\partial \rho}{\partial \theta}$ is the thermal expansion coefficient and $\beta = \frac{1}{\rho_0} \frac{\partial \rho}{\partial S}$ is the haline contraction coefficient and vary in space and time according to the temperature, salinity and pressure fields. Then, the ECCO V4r3 state estimate gives all the necessary terms to decompose the potential temperature and salinity evolution equation into advection, diffusion and surface fluxes terms (net heat flux for the potential temperature and freshwater flux for the salinity):

$$ \left( \frac{\partial \theta}{\partial t}, \frac{\partial S}{\partial t} \right) = (\text{adv}_{\theta}, \text{adv}_{S}) + (\text{dif}_{\theta}, \text{dif}_{S}) + (f_{\theta}, f_{S}) \qquad (9) $$

where $adv_{\theta}$, $adv_{S}$ are respectively the advective terms for temperature and salinity, $dif_{\theta}$, $dif_{S}$, the diffusive terms for temperature and salinity and $f_{\theta}$, $f_{S}$, the atmospheric forcing terms for temperature and salinity.

Using this decomposition (i.e. Eq. (9)) of the temperature and salinity evolution equation and Eq. (8), the time evolution of $\rho$ is decomposed itself into advection $adv_{\rho} = -\rho_0 \alpha adv_{\theta} + \rho_0 \beta adv_{S}$, parametrized diffusion $dif_{\rho} = -\rho_0 \alpha dif_{\theta} + \rho_0 \beta dif_{S}$, and buoyancy fluxes from the atmosphere $flu_{\rho} = -\rho_0 \alpha f_{\theta} + \rho_0 \beta f_{S}$:

$$ \frac{\partial \rho}{\partial t} = \text{adv}_{\rho} + \text{dif}_{\rho} + \text{flu}_{\rho}. \qquad (10) $$

Then, inserting Eq. (10) in Eq. (7) gives:

$$ \frac{\partial \eta}{\partial t} = \text{adv} + \text{dif} + \text{flu}. \qquad (11) $$

where $X = \text{adv, dif, or flu}$ is related to $X_{\rho}$ through the following formula:

$$ X = -\frac{1}{\rho_0} \int_{-H}^{0} X_{\rho} dz. \qquad (12) $$

Following previous work on density variance³¹,³³,³⁴,⁴⁴, temperature variance⁴⁵,⁴⁶ and also steric sea level variance³⁰, multiplying equation (11) by $\eta$ and computing the time average (over all months of the mean seasonal cycle) gives the budget of the seasonal cycle steric sea level variance:
term has amplitudes and patterns that are very different from the term linked with the net heat flux (Fig. 1a) and it therefore demonstrates that Eq. (1) is not complete.

## The ECCOv4r3 dataset

Seasonal SSL variance budgets are computed using the ECCOv4r3 state estimate which covers the period 1992–2017. This state estimate is the output of the Massachusetts Institute of Technology general circulation model (MITgcm) assimilating available observations for the period 1992 to 2017³⁸. The advantage of ECCOv4 is that it satisfies the equation of motion and conservation laws hence making it possible to compute tracer budgets. The solution used in this article is computed on the LLC90 grid which has an average horizontal resolution of 1° and 50 vertical levels. Outputs of the model consist of one month average and the closed budget can be obtained for the 1993–2014 period. Thus we compute the SSL budget over this 1993–2014 period. Snapshots at the start and end of each month are also provided by ECCO in order to compute the tracer time tendencies required to close the budgets. ECCO has already been used in several past studies to compute budget of steric sea level variations²⁷⁻³⁰.

For each variable from the monthly model outputs, at each gridpoint, the seasonal cycle is obtained by first removing the linear trend over the 1993–2014 period, and then by computing the 12 monthly values as average of the respective monthly anomalies over the 22 years. The SSL seasonal cycle studied here is thus an average of the seasonal cycle over the 22 years of the 1993–2014 period.

## Seasonal SSL variance budget

The seasonal SSL anomaly $\eta$ is expressed as the vertical integral of the seasonal density anomaly $\rho$ as follows:

$$ \eta = -\frac{1}{\rho_0} \int_{-H}^{0} \rho \, dz, \qquad (6) $$

where $\rho_0 = 1029 \text{ kg m}^{-3}$. The evolution equation for $\eta$ is then simply obtained as the time derivative of (6):

$$ \frac{\partial \eta}{\partial t} = -\frac{1}{\rho_0} \int_{-H}^{0} \frac{\partial \rho}{\partial t} dz. \qquad (7) $$

Following³⁰, the time evolution of $\rho$ can be decomposed into the time evolution of the potential temperature $\theta$ and salinity $S$ as:

$$ \frac{\partial \rho}{\partial t} = -\rho_0 \alpha \frac{\partial \theta}{\partial t} + \rho_0 \beta \frac{\partial S}{\partial t}, \qquad (8) $$

where $\alpha = -\frac{1}{\rho_0} \frac{\partial \rho}{\partial \theta}$ is the thermal expansion coefficient and $\beta = \frac{1}{\rho_0} \frac{\partial \rho}{\partial S}$ is the haline contraction coefficient and vary in space and time according to the temperature, salinity and pressure fields. Then, the ECCO V4r3 state estimate gives all the necessary terms to decompose the potential temperature and salinity evolution equation into advection, diffusion and surface fluxes terms (net heat flux for the potential temperature and freshwater flux for the salinity):

$$ \left( \frac{\partial \theta}{\partial t}, \frac{\partial S}{\partial t} \right) = (\text{adv}_{\theta}, \text{adv}_{S}) + (\text{dif}_{\theta}, \text{dif}_{S}) + (f_{\theta}, f_{S}) \qquad (9) $$

where $adv_{\theta}$, $adv_{S}$ are respectively the advective terms for temperature and salinity, $dif_{\theta}$, $dif_{S}$, the diffusive terms for temperature and salinity and $f_{\theta}$, $f_{S}$, the atmospheric forcing terms for temperature and salinity.

Using this decomposition (i.e. Eq. (9)) of the temperature and salinity evolution equation and Eq. (8), the time evolution of $\rho$ is decomposed itself into advection $adv_{\rho} = -\rho_0 \alpha adv_{\theta} + \rho_0 \beta adv_{S}$, parametrized diffusion $dif_{\rho} = -\rho_0 \alpha dif_{\theta} + \rho_0 \beta dif_{S}$, and buoyancy fluxes from the atmosphere $flu_{\rho} = -\rho_0 \alpha f_{\theta} + \rho_0 \beta f_{S}$:

$$ \frac{\partial \rho}{\partial t} = \text{adv}_{\rho} + \text{dif}_{\rho} + \text{flu}_{\rho}. \qquad (10) $$

Then, inserting Eq. (10) in Eq. (7) gives:

$$ \frac{\partial \eta}{\partial t} = \text{adv} + \text{dif} + \text{flu}. \qquad (11) $$

where $X = \text{adv, dif, or flu}$ is related to $X_{\rho}$ through the following formula:

$$ X = -\frac{1}{\rho_0} \int_{-H}^{0} X_{\rho} dz. \qquad (12) $$

Following previous work on density variance³¹,³³,³⁴,⁴⁴, temperature variance⁴⁵,⁴⁶ and also steric sea level variance³⁰, multiplying equation (11) by $\eta$ and computing the time average (over all months of the mean seasonal cycle) gives the budget of the seasonal cycle steric sea level variance:
$$ \begin{aligned} \overline{\eta \frac{\partial \eta}{\partial t}} &= \overline{\frac{1}{2} \frac{\partial \eta^2}{\partial t}} = \overline{\eta \text{adv}} + \overline{\eta \text{dif}} + \overline{\eta \text{flu}}, \\ &= \text{VAR}_{\text{adv}} + \text{VAR}_{\text{dif}} + \text{VAR}_{\text{flu}}, \end{aligned} \quad (13) $$

where the overline represents the average over the seasonal cycle and $\overline{\eta X}$ is a 2D field, it is positive (negative) when $\overline{\eta X}$ is a source (sink) of seasonal SSL variance. By examining the sign and relative intensity of the seasonal SSL variance budget terms, it is then possible to determine which term is locally driving or damping the seasonal variation of the SSL. $\overline{\eta X}$ can be a source (sink) in two cases: 1) if $\eta > 0$ and $-\int_{-H}^{0} X_{\rho} dz > 0$ ($-\int_{-H}^{0} X_{\rho} dz < 0$) because the last term acts to increase (decrease) the positive anomaly of $\eta$ or 2) if $\eta < 0$ and $-\int_{-H}^{0} X_{\rho} dz < 0$ ($-\int_{-H}^{0} X_{\rho} dz > 0$) because the last term acts to increase (decrease) the magnitude of the negative anomaly of $\eta$. Note that, similar to the analysis performed on SLA from AVISO, in practice the term $\frac{\partial \eta^2}{\partial t}$ is not exactly zero because the seasonal cycle of the SSL time tendencies is not exactly equal to the time tendencies of the SSL seasonal cycle. However, Figure 3a shows that this term remains one order of magnitude smaller than the other terms.

## Variance budget for the spatially-averaged seasonal steric sea level

In this subsection, we compute the variance budget for the spatially-averaged seasonal steric sea level. To this end, the seasonal steric sea level anomaly $\eta$ is first spatially averaged $\langle \eta \rangle$, where $\langle \cdot \rangle$ represents the spatial average over the Northern hemisphere for latitudes North of a given limit value or over the Southern hemisphere for latitudes South of a given boundary. Then the equation for the evolution of $\eta$ is also spatially averaged over the same region as follows:

$$ \frac{\partial \langle \eta \rangle}{\partial t} = \langle \text{adv} \rangle + \langle \text{dif} \rangle + \langle \text{flu} \rangle \quad (14) $$

Multiplying this equation by $\langle \eta \rangle$ gives the budget for the variance of $\langle \eta \rangle$:

$$ \begin{aligned} \frac{1}{2} \frac{\partial \langle \eta \rangle^2}{\partial t} &= \overline{\langle \eta \rangle \langle \text{adv} \rangle} + \overline{\langle \eta \rangle \langle \text{dif} \rangle} + \overline{\langle \eta \rangle \langle \text{flu} \rangle}, \\ &= \text{VAR}_{\text{adv}}^{\text{H}} + \text{VAR}_{\text{dif}}^{\text{H}} + \text{VAR}_{\text{flu}}^{\text{H}}, \end{aligned} \quad (15) $$

$\overline{\frac{1}{2} \frac{\partial \langle \eta \rangle^2}{\partial t}}$ is the variance tendency term, $\text{VAR}_{\text{adv}}^{\text{H}} = \overline{\langle \eta \rangle \langle \text{adv} \rangle}$, $\text{VAR}_{\text{dif}}^{\text{H}} = \overline{\langle \eta \rangle \langle \text{dif} \rangle}$ and $\text{VAR}_{\text{flu}}^{\text{H}} = \overline{\langle \eta \rangle \langle \text{flu} \rangle}$ are respectively the terms associated with advection, diffusion and surface buoyancy fluxes for the variance equation of the SSL hemispheric average.

## Vertically-averaged density and VARflu

Inserting Eqs. (6) in (12), and recognizing that the atmospheric density flux is mostly due to the net atmospheric heat flux $q$ (in W m$^{-3}$) (see supporting information Fig. S2), VARflu can be written as:

$$ \text{VAR}_{\text{flu}} = \overline{\eta \text{flu}} = \frac{1}{\rho_0^2} \int_{-H}^{0} \rho \, dz \int_{-H}^{0} -\frac{\alpha}{c_p} q \, dz. \quad (16) $$

Then, defining the vertically-averaged density anomaly as $\rho_{\text{BT}} = \frac{1}{H} \int_{-H}^{0} \rho dz$ and the vertically-averaged net heat flux as $Q_{\text{BT}} = \frac{1}{H} \int_{-H}^{0} -\frac{\alpha}{c_p} q dz$ ($q$ is non-zero below the surface because of the penetrating nature of the shortwave term) in this equation one recovers Eq. (4).

## Advective term decomposition

To decompose the advective term in the seasonal SSL budget, we first decompose the advective term in the density evolution equation (10) into the contributions of laminar (resolved) velocities and eddy-induced velocities linked to the Gent and McWilliams parameterization:

$$ \text{adv}_{\rho} = \text{adv}_{\rho}^{\text{laminar}} + \text{adv}_{\rho}^{\text{eddy}}, \quad (17) $$

where $\text{adv}_{\rho}^{\text{laminar}}$ is the ECCO v4r3 resolved part of the seasonal advective term and $\text{adv}_{\rho}^{\text{eddy}}$ the eddy induced part of the seasonal advective term. Using Eq. (17), $\text{VAR}_{\text{adv}} = \eta \text{adv}$ can thus be decomposed into two terms $\text{VAR}_{\text{adv}}^{\text{laminar}} = \eta \text{adv}_{\text{laminar}}$ and $\text{VAR}_{\text{adv}}^{\text{eddy}} = \eta \text{adv}_{\text{eddy}}$, as shown by formula (5) where the notation $\text{VAR}_{\text{adv}}^{\text{X}} = \eta \text{adv}_{\text{X}}$ is given by:

$$ \text{VAR}_{\text{adv}}^{\text{X}} = \overline{\eta \text{adv}_{\text{X}}} = -\frac{\eta}{\rho_0} \int_{-H}^{0} \text{adv}_{\rho}^{\text{X}} dz, \quad (18) $$

where X stands for laminar or eddy.

Received: 28 February 2024; Accepted: 6 May 2024

Published online: 09 May 2024
$$ \begin{aligned} \overline{\eta \frac{\partial \eta}{\partial t}} &= \overline{\frac{1}{2} \frac{\partial \eta^2}{\partial t}} = \overline{\eta \text{adv}} + \overline{\eta \text{dif}} + \overline{\eta \text{flu}}, \\ &= \text{VAR}_{\text{adv}} + \text{VAR}_{\text{dif}} + \text{VAR}_{\text{flu}}, \end{aligned} \quad (13) $$

where the overline represents the average over the seasonal cycle and $\overline{\eta X}$ is a 2D field, it is positive (negative) when $\overline{\eta X}$ is a source (sink) of seasonal SSL variance. By examining the sign and relative intensity of the seasonal SSL variance budget terms, it is then possible to determine which term is locally driving or damping the seasonal variation of the SSL. $\overline{\eta X}$ can be a source (sink) in two cases: 1) if $\eta > 0$ and $-\int_{-H}^{0} X_{\rho} dz > 0$ ($-\int_{-H}^{0} X_{\rho} dz < 0$) because the last term acts to increase (decrease) the positive anomaly of $\eta$ or 2) if $\eta < 0$ and $-\int_{-H}^{0} X_{\rho} dz < 0$ ($-\int_{-H}^{0} X_{\rho} dz > 0$) because the last term acts to increase (decrease) the magnitude of the negative anomaly of $\eta$. Note that, similar to the analysis performed on SLA from AVISO, in practice the term $\frac{\partial \eta^2}{\partial t}$ is not exactly zero because the seasonal cycle of the SSL time tendencies is not exactly equal to the time tendencies of the SSL seasonal cycle. However, Figure 3a shows that this term remains one order of magnitude smaller than the other terms.

## Variance budget for the spatially-averaged seasonal steric sea level

In this subsection, we compute the variance budget for the spatially-averaged seasonal steric sea level. To this end, the seasonal steric sea level anomaly $\eta$ is first spatially averaged $\langle \eta \rangle$, where $\langle \cdot \rangle$ represents the spatial average over the Northern hemisphere for latitudes North of a given limit value or over the Southern hemisphere for latitudes South of a given boundary. Then the equation for the evolution of $\eta$ is also spatially averaged over the same region as follows:

$$ \frac{\partial \langle \eta \rangle}{\partial t} = \langle \text{adv} \rangle + \langle \text{dif} \rangle + \langle \text{flu} \rangle \quad (14) $$

Multiplying this equation by $\langle \eta \rangle$ gives the budget for the variance of $\langle \eta \rangle$:

$$ \begin{aligned} \frac{1}{2} \frac{\partial \langle \eta \rangle^2}{\partial t} &= \overline{\langle \eta \rangle \langle \text{adv} \rangle} + \overline{\langle \eta \rangle \langle \text{dif} \rangle} + \overline{\langle \eta \rangle \langle \text{flu} \rangle}, \\ &= \text{VAR}_{\text{adv}}^{\text{H}} + \text{VAR}_{\text{dif}}^{\text{H}} + \text{VAR}_{\text{flu}}^{\text{H}}, \end{aligned} \quad (15) $$

$\overline{\frac{1}{2} \frac{\partial \langle \eta \rangle^2}{\partial t}}$ is the variance tendency term, $\text{VAR}_{\text{adv}}^{\text{H}} = \overline{\langle \eta \rangle \langle \text{adv} \rangle}$, $\text{VAR}_{\text{dif}}^{\text{H}} = \overline{\langle \eta \rangle \langle \text{dif} \rangle}$ and $\text{VAR}_{\text{flu}}^{\text{H}} = \overline{\langle \eta \rangle \langle \text{flu} \rangle}$ are respectively the terms associated with advection, diffusion and surface buoyancy fluxes for the variance equation of the SSL hemispheric average.

## Vertically-averaged density and VARflu

Inserting Eqs. (6) in (12), and recognizing that the atmospheric density flux is mostly due to the net atmospheric heat flux $q$ (in W m$^{-3}$) (see supporting information Fig. S2), VARflu can be written as:

$$ \text{VAR}_{\text{flu}} = \overline{\eta \text{flu}} = \frac{1}{\rho_0^2} \int_{-H}^{0} \rho \, dz \int_{-H}^{0} -\frac{\alpha}{c_p} q \, dz. \quad (16) $$

Then, defining the vertically-averaged density anomaly as $\rho_{\text{BT}} = \frac{1}{H} \int_{-H}^{0} \rho dz$ and the vertically-averaged net heat flux as $Q_{\text{BT}} = \frac{1}{H} \int_{-H}^{0} -\frac{\alpha}{c_p} q dz$ ($q$ is non-zero below the surface because of the penetrating nature of the shortwave term) in this equation one recovers Eq. (4).

## Advective term decomposition

To decompose the advective term in the seasonal SSL budget, we first decompose the advective term in the density evolution equation (10) into the contributions of laminar (resolved) velocities and eddy-induced velocities linked to the Gent and McWilliams parameterization:

$$ \text{adv}_{\rho} = \text{adv}_{\rho}^{\text{laminar}} + \text{adv}_{\rho}^{\text{eddy}}, \quad (17) $$

where $\text{adv}_{\rho}^{\text{laminar}}$ is the ECCO v4r3 resolved part of the seasonal advective term and $\text{adv}_{\rho}^{\text{eddy}}$ the eddy induced part of the seasonal advective term. Using Eq. (17), $\text{VAR}_{\text{adv}} = \eta \text{adv}$ can thus be decomposed into two terms $\text{VAR}_{\text{adv}}^{\text{laminar}} = \eta \text{adv}_{\text{laminar}}$ and $\text{VAR}_{\text{adv}}^{\text{eddy}} = \eta \text{adv}_{\text{eddy}}$, as shown by formula (5) where the notation $\text{VAR}_{\text{adv}}^{\text{X}} = \eta \text{adv}_{\text{X}}$ is given by:

$$ \text{VAR}_{\text{adv}}^{\text{X}} = \overline{\eta \text{adv}_{\text{X}}} = -\frac{\eta}{\rho_0} \int_{-H}^{0} \text{adv}_{\rho}^{\text{X}} dz, \quad (18) $$

where X stands for laminar or eddy.

Received: 28 February 2024; Accepted: 6 May 2024

Published online: 09 May 2024
References

1. Oppenheimer, M. et al. Sea Level Rise and Implications for Low Lying Islands, Coasts and Communities (The Intergovernmental Panel on Climate Change, 2019).

2. Leuliette, E. W. & Willis, J. K. Balancing the sea level budget. *Oceanography* **24**(2), 122–129 (2011).

3. Gregory, J. M., Griffies, S. M. & Hughes, C. W. Concepts and terminology for sea level: Mean, variability and change, both local and global. *Surv. Geophys.* **40**, 1251–1289. https://doi.org/10.1007/s10712-019-09525-z (2019).

4. Storto, A. et al. Steric sea level variability (1993–2010) in an ensemble of ocean reanalyses and objective analyses. *Clim. Dyn.* **9(3)**, 709–729 (2017).

5. Gulev, S. K. et al. Changing state of the climate system. In *Climate Change 2021: The Physical Science Basis* 287–422 (Cambridge University Press, 2021). https://doi.org/10.1017/9781009157896.004.

6. Lombard, A. et al. Estimation of steric sea level variations from combined GRACE and Jason-1 data. *Earth Planet. Sci. Lett.* **254**(1–2), 194–202 (2007).

7. Copernicus Climate Change Service, Climate Data Store, 2018: Sea level gridded data from satellite observations for the global ocean from 1993 to present. Copernicus Climate Change Service (C3S) Climate Data Store (CDS). https://doi.org/10.24381/cds.4c328c78. Accessed 25-07-2023 (2018).

8. Gill, A. & Niller, P. The theory of the seasonal variability in the ocean. *Deep Sea Res. Oceanogr. Abstr.* **20**(2), 141–177 (1973).

9. Cheney, R., Miller, L., Agreen, R., Doyle, N. & Lillibridge, J. Topex/poseidon: The 2-cm solution. *J. Geophys. Res. Oceans* **99**(C12), 24555–24563. https://doi.org/10.1029/94JC02050 (1994).

10. Stammer, D. Steric and wind-induced changes in TOPEX/POSEIDON large-scale sea surface topography observations. *J. Geophys. Res. Oceans* **102**(C9), 20987–21009. https://doi.org/10.1029/97J01475 (1997).

11. Vinogradov, S. V., Ponte, R. M., Heimbach, P. & Wunsch, C. The mean seasonal cycle in sea level estimated from a data-constrained general circulation model. *J. Geophys. Res. Oceans* **113**(C3), 03032. https://doi.org/10.1029/2007JC004496 (2008).

12. Piecuch, C. G. & Ponte, R. M. Mechanisms of global-mean steric sea level change. *J. Clim.* **27**(2), 824–834. https://doi.org/10.1175/JCLI-D-13-00373.1 (2014).

13. Périgaud, C. & Delecluse, P. Annual sea level variations in the southern tropical Indian Ocean from Geosat and shallow-water simulations. *J. Geophys. Res. Oceans* **97**(C12), 20169–20178. https://doi.org/10.1029/92JC01961 (1992).

14. Vivier, F., Kelly, K. A. & Thompson, L. Contributions of wind forcing, waves, and surface heating to sea surface height observations in the Pacific Ocean. *J. Geophys. Res. Oceans* **104**(C9), 20767–20788. https://doi.org/10.1029/1999JC900096 (1999).

15. Bowen, M. M., Sutton, P. J. & Roemmich, D. Wind-driven and steric fluctuations of sea surface height in the southwest Pacific. *Geophys. Res. Lett.* **33**(14), 14617 (2006).

16. Cabanes, C., Huck, T. & Colin de Verdière, A. Contributions of wind forcing and surface heating to interannual sea level variations in the Atlantic Ocean. *J. Phys. Oceanogr.* **36**(9), 1739–1750. https://doi.org/10.1175/JPO2935.1 (2006).

17. Piecuch, C. G. & Ponte, R. M. Importance of circulation changes to Atlantic heat storage rates on seasonal and interannual time scales. *J. Clim.* **25**(1), 350–362. https://doi.org/10.1175/JCLI-D-11-00123.1 (2012).

18. Giglio, D., Roemmich, D. & Cornuelle, B. Understanding the annual cycle in global steric height. *Geophys. Res. Lett.* **40**(16), 4349–4354. https://doi.org/10.1002/grl.50774 (2013).

19. Roberts, C. D. et al. Surface flux and ocean heat transport convergence contributions to seasonal and interannual variations of ocean heat content. *J. Geophys. Res. Oceans* **122**(1), 726–744. https://doi.org/10.1002/2016JC012278 (2017).

20. Ferry, N., Reverdin, G. & Oschlies, A. Seasonal sea surface height variability in the North Atlantic Ocean. *J. Geophys. Res. Oceans* **105**(C3), 6307–6326. https://doi.org/10.1029/1999JC900296 (2000).

21. García-Lafuente, J., Del Río, J., Alvarez Fanjul, E., Gomis, D. & Delgado, J. Some aspects of the seasonal sea level variations around Spain. *J. Geophys. Res. Oceans* **109**(C9), 09008. https://doi.org/10.1029/2003JC002070 (2004).

22. Mork, K. A. & Skagseth, O. Annual sea surface height variability in the nordic seas. *Geophys. Monogr. Am. Geophys. Union* **158**, 51 (2005).

23. Roemmich, D. & Gilson, J. The 2004–2008 mean and annual cycle of temperature, salinity, and steric height in the global ocean from the Argo Program. *Prog. Oceanogr.* **82**(2), 81–100 (2009).

24. Torres, R. R. & Tsimplis, M. N. Seasonal sea level cycle in the Caribbean Sea. *J. Geophys. Res. Oceans* **117**(C7), 07011. https://doi.org/10.1029/2012JC008159 (2012).

25. Hersbach, H. et al. The ERA5 global reanalysis. *Q. J. R. Meteorol. Soc.* **146**(730), 1999–2049 (2020).

26. Roberts, C. D. et al. Surface flux and ocean heat transport convergence contributions to seasonal and interannual variations of ocean heat content. *J. Geophys. Res. Oceans* **122**(1), 726–744. https://doi.org/10.1002/2016JC012278 (2017).

27. Piecuch, C. & Ponte, R. Mechanisms of interannual steric sea level variability. *Geophys. Res. Lett.* **38**(15), 15605. https://doi.org/10.1029/2011GL048440 (2011).

28. Forget, G. & Ponte, R. M. The partition of regional sea level variability. *Prog. Oceanogr.* **137**, 173–195 (2015).

29. Meyssignac, B., Piecuch, C., Merchant, C., Racault, M.-F., Palanisamy, H., MacIntosh, C., Sathyendranath, S. & Brewin, R. Causes of the regional variability in observed sea level, sea surface temperature and ocean colour over the period 1993–2011. Integrative Study of the Mean Sea Level and Its Components, 191–219 (2017).

30. Hochet, A., Llovel, W., Sévellec, F. & Huck, T. Sources and sinks of interannual steric sea level variability. *J. Geophys. Res. Oceans* **128**(4), 2022–019335. https://doi.org/10.1029/2022JC019335 (2023).

31. Colin de Verdière, A. & Huck, T. Baroclinic instability: An oceanic wavemaker for interdecadal variability. *J. Phys. Oceanogr.* **29**(5), 893–910 (1999).

32. Arzel, O., Huck, T. & Colin de Verdière, A. The different nature of the interdecadal variability of the thermohaline circulation under mixed and flux boundary conditions. *J. Phys. Oceanogr.* **36**(9), 1703–1718 (2006).

33. Buckley, M. W., Ferreira, D., Campin, J.-M., Marshall, J. & Tulloch, R. On the relationship between decadal buoyancy anomalies and variability of the Atlantic meridional overturning circulation. *J. Clim.* **25**(23), 8009–8030 (2012).

34. Arzel, O., Huck, T. & Colin de Verdière, A. The internal generation of the Atlantic Ocean interdecadal variability. *J. Clim.* **31**(16), 6411–6432 (2018).

35. Osborn, T. R. & Cox, C. S. Oceanic fine structure. *Geophys. Fluid Dyn.* **3**(4), 321–345. https://doi.org/10.1080/0309192708236085 (1972).

36. MacCready, P., Geyer, W. R. & Burchard, H. Estuarine exchange flow is related to mixing through the salinity variance budget. *J. Phys. Oceanogr.* **48**(6), 1375–1384 (2018).

37. Gregg, M. C., D'Asaro, E. A., Riley, J. J. & Kunze, E. Mixing efficiency in the ocean. *Annu. Rev. Mar. Sci.* **10**(1), 443–473. https://doi.org/10.1146/annurev-marine-121916-063643 (2018).

38. Forget, G. et al. Ecco version 4: An integrated framework for non-linear inverse modeling and global ocean state estimation. *Geosci. Model Dev.* **8**(10), 3071–3104 (2015).

39. Gent, P. R. & McWilliams, J. C. Isopycnal mixing in ocean circulation models. *J. Phys. Oceanogr.* **20**(1), 150–155 (1990).

40. Wang, S. et al. El Niño/Southern Oscillation inhibited by submesoscale ocean eddies. *Nat. Geosci.* **15**(2), 112–117. https://doi.org/10.1038/s41561-021-00890-2 (2022).

41. Carret, A., Llovel, W., Penduff, T. & Molines, J.-M. Atmospherically forced and chaotic interannual variability of regional sea level and its components over 1993–2015. *J. Geophys. Res. Oceans* **126**(4), 2020–017123. https://doi.org/10.1029/2022JC019335 (2021).
References

1. Oppenheimer, M. et al. Sea Level Rise and Implications for Low Lying Islands, Coasts and Communities (The Intergovernmental Panel on Climate Change, 2019).

2. Leuliette, E. W. & Willis, J. K. Balancing the sea level budget. *Oceanography* **24**(2), 122–129 (2011).

3. Gregory, J. M., Griffies, S. M. & Hughes, C. W. Concepts and terminology for sea level: Mean, variability and change, both local and global. *Surv. Geophys.* **40**, 1251–1289. https://doi.org/10.1007/s10712-019-09525-z (2019).

4. Storto, A. et al. Steric sea level variability (1993–2010) in an ensemble of ocean reanalyses and objective analyses. *Clim. Dyn.* **9(3)**, 709–729 (2017).

5. Gulev, S. K. et al. Changing state of the climate system. In *Climate Change 2021: The Physical Science Basis* 287–422 (Cambridge University Press, 2021). https://doi.org/10.1017/9781009157896.004.

6. Lombard, A. et al. Estimation of steric sea level variations from combined GRACE and Jason-1 data. *Earth Planet. Sci. Lett.* **254**(1–2), 194–202 (2007).

7. Copernicus Climate Change Service, Climate Data Store, 2018: Sea level gridded data from satellite observations for the global ocean from 1993 to present. Copernicus Climate Change Service (C3S) Climate Data Store (CDS). https://doi.org/10.24381/cds.4c328c78. Accessed 25-07-2023 (2018).

8. Gill, A. & Niller, P. The theory of the seasonal variability in the ocean. *Deep Sea Res. Oceanogr. Abstr.* **20**(2), 141–177 (1973).

9. Cheney, R., Miller, L., Agreen, R., Doyle, N. & Lillibridge, J. Topex/poseidon: The 2-cm solution. *J. Geophys. Res. Oceans* **99**(C12), 24555–24563. https://doi.org/10.1029/94JC02050 (1994).

10. Stammer, D. Steric and wind-induced changes in TOPEX/POSEIDON large-scale sea surface topography observations. *J. Geophys. Res. Oceans* **102**(C9), 20987–21009. https://doi.org/10.1029/97J01475 (1997).

11. Vinogradov, S. V., Ponte, R. M., Heimbach, P. & Wunsch, C. The mean seasonal cycle in sea level estimated from a data-constrained general circulation model. *J. Geophys. Res. Oceans* **113**(C3), 03032. https://doi.org/10.1029/2007JC004496 (2008).

12. Piecuch, C. G. & Ponte, R. M. Mechanisms of global-mean steric sea level change. *J. Clim.* **27**(2), 824–834. https://doi.org/10.1175/JCLI-D-13-00373.1 (2014).

13. Périgaud, C. & Delecluse, P. Annual sea level variations in the southern tropical Indian Ocean from Geosat and shallow-water simulations. *J. Geophys. Res. Oceans* **97**(C12), 20169–20178. https://doi.org/10.1029/92JC01961 (1992).

14. Vivier, F., Kelly, K. A. & Thompson, L. Contributions of wind forcing, waves, and surface heating to sea surface height observations in the Pacific Ocean. *J. Geophys. Res. Oceans* **104**(C9), 20767–20788. https://doi.org/10.1029/1999JC900096 (1999).

15. Bowen, M. M., Sutton, P. J. & Roemmich, D. Wind-driven and steric fluctuations of sea surface height in the southwest Pacific. *Geophys. Res. Lett.* **33**(14), 14617 (2006).

16. Cabanes, C., Huck, T. & Colin de Verdière, A. Contributions of wind forcing and surface heating to interannual sea level variations in the Atlantic Ocean. *J. Phys. Oceanogr.* **36**(9), 1739–1750. https://doi.org/10.1175/JPO2935.1 (2006).

17. Piecuch, C. G. & Ponte, R. M. Importance of circulation changes to Atlantic heat storage rates on seasonal and interannual time scales. *J. Clim.* **25**(1), 350–362. https://doi.org/10.1175/JCLI-D-11-00123.1 (2012).

18. Giglio, D., Roemmich, D. & Cornuelle, B. Understanding the annual cycle in global steric height. *Geophys. Res. Lett.* **40**(16), 4349–4354. https://doi.org/10.1002/grl.50774 (2013).

19. Roberts, C. D. et al. Surface flux and ocean heat transport convergence contributions to seasonal and interannual variations of ocean heat content. *J. Geophys. Res. Oceans* **122**(1), 726–744. https://doi.org/10.1002/2016JC012278 (2017).

20. Ferry, N., Reverdin, G. & Oschlies, A. Seasonal sea surface height variability in the North Atlantic Ocean. *J. Geophys. Res. Oceans* **105**(C3), 6307–6326. https://doi.org/10.1029/1999JC900296 (2000).

21. García-Lafuente, J., Del Río, J., Alvarez Fanjul, E., Gomis, D. & Delgado, J. Some aspects of the seasonal sea level variations around Spain. *J. Geophys. Res. Oceans* **109**(C9), 09008. https://doi.org/10.1029/2003JC002070 (2004).

22. Mork, K. A. & Skagseth, O. Annual sea surface height variability in the nordic seas. *Geophys. Monogr. Am. Geophys. Union* **158**, 51 (2005).

23. Roemmich, D. & Gilson, J. The 2004–2008 mean and annual cycle of temperature, salinity, and steric height in the global ocean from the Argo Program. *Prog. Oceanogr.* **82**(2), 81–100 (2009).

24. Torres, R. R. & Tsimplis, M. N. Seasonal sea level cycle in the Caribbean Sea. *J. Geophys. Res. Oceans* **117**(C7), 07011. https://doi.org/10.1029/2012JC008159 (2012).

25. Hersbach, H. et al. The ERA5 global reanalysis. *Q. J. R. Meteorol. Soc.* **146**(730), 1999–2049 (2020).

26. Roberts, C. D. et al. Surface flux and ocean heat transport convergence contributions to seasonal and interannual variations of ocean heat content. *J. Geophys. Res. Oceans* **122**(1), 726–744. https://doi.org/10.1002/2016JC012278 (2017).

27. Piecuch, C. & Ponte, R. Mechanisms of interannual steric sea level variability. *Geophys. Res. Lett.* **38**(15), 15605. https://doi.org/10.1029/2011GL048440 (2011).

28. Forget, G. & Ponte, R. M. The partition of regional sea level variability. *Prog. Oceanogr.* **137**, 173–195 (2015).

29. Meyssignac, B., Piecuch, C., Merchant, C., Racault, M.-F., Palanisamy, H., MacIntosh, C., Sathyendranath, S. & Brewin, R. Causes of the regional variability in observed sea level, sea surface temperature and ocean colour over the period 1993–2011. Integrative Study of the Mean Sea Level and Its Components, 191–219 (2017).

30. Hochet, A., Llovel, W., Sévellec, F. & Huck, T. Sources and sinks of interannual steric sea level variability. *J. Geophys. Res. Oceans* **128**(4), 2022–019335. https://doi.org/10.1029/2022JC019335 (2023).

31. Colin de Verdière, A. & Huck, T. Baroclinic instability: An oceanic wavemaker for interdecadal variability. *J. Phys. Oceanogr.* **29**(5), 893–910 (1999).

32. Arzel, O., Huck, T. & Colin de Verdière, A. The different nature of the interdecadal variability of the thermohaline circulation under mixed and flux boundary conditions. *J. Phys. Oceanogr.* **36**(9), 1703–1718 (2006).

33. Buckley, M. W., Ferreira, D., Campin, J.-M., Marshall, J. & Tulloch, R. On the relationship between decadal buoyancy anomalies and variability of the Atlantic meridional overturning circulation. *J. Clim.* **25**(23), 8009–8030 (2012).

34. Arzel, O., Huck, T. & Colin de Verdière, A. The internal generation of the Atlantic Ocean interdecadal variability. *J. Clim.* **31**(16), 6411–6432 (2018).

35. Osborn, T. R. & Cox, C. S. Oceanic fine structure. *Geophys. Fluid Dyn.* **3**(4), 321–345. https://doi.org/10.1080/0309192708236085 (1972).

36. MacCready, P., Geyer, W. R. & Burchard, H. Estuarine exchange flow is related to mixing through the salinity variance budget. *J. Phys. Oceanogr.* **48**(6), 1375–1384 (2018).

37. Gregg, M. C., D'Asaro, E. A., Riley, J. J. & Kunze, E. Mixing efficiency in the ocean. *Annu. Rev. Mar. Sci.* **10**(1), 443–473. https://doi.org/10.1146/annurev-marine-121916-063643 (2018).

38. Forget, G. et al. Ecco version 4: An integrated framework for non-linear inverse modeling and global ocean state estimation. *Geosci. Model Dev.* **8**(10), 3071–3104 (2015).

39. Gent, P. R. & McWilliams, J. C. Isopycnal mixing in ocean circulation models. *J. Phys. Oceanogr.* **20**(1), 150–155 (1990).

40. Wang, S. et al. El Niño/Southern Oscillation inhibited by submesoscale ocean eddies. *Nat. Geosci.* **15**(2), 112–117. https://doi.org/10.1038/s41561-021-00890-2 (2022).

41. Carret, A., Llovel, W., Penduff, T. & Molines, J.-M. Atmospherically forced and chaotic interannual variability of regional sea level and its components over 1993–2015. *J. Geophys. Res. Oceans* **126**(4), 2020–017123. https://doi.org/10.1029/2022JC019335 (2021).
42. Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S.L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. et al. Climate change 2021: The physical science basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, vol. 2 (2021).

43. Yu, L. & Weller, R. A. Objectively analyzed air-sea heat fluxes for the global ice-free oceans (1981–2005). *Bull. Am. Meteorol. Soc.* **88**(4), 527–540 (2007).

44. Hochet, A., Huck, T. & Colin de Verdière, A. Large-scale baroclinic instability of the mean oceanic circulation: A local approach. *J. Phys. Oceanogr.* **45**(11), 2738–2754. https://doi.org/10.1175/JPO-D-15-0084.1 (2015).

45. Boucharel, J. et al. A surface layer variance heat budget for ENSO. *Geophys. Res. Lett.* **42**(9), 3529–3537. https://doi.org/10.1002/2015GL063843 (2015).

46. Hochet, A. et al. Direct temporal cascade of temperature variance in eddy-permitting simulations of multidecadal variability. *J. Clim.* **33**(21), 9409–9425. https://doi.org/10.1175/JCLI-D-19-0921.1 (2020).

47. Met Office: Cartopy: A Cartographic Python Library with a Matplotlib Interface. Exeter, Devon (2010–2015). https://scitools.org.uk/cartopy

## Acknowledgements

The datasets analysed during the current study are available in the ECCO repository, https://web.corral.tacc.utexas.edu/OceanProjects/ECCO/ECCOv4/Release3/. ECCO v4r3 salinity and temperature budget computations were greatly facilitated by the existence of the ecco_v4 py Python library, and by several tutorials which can be found at: https://ecco-v4-python-tutorial.readthedocs.io/. We warmly acknowledge the authors of this work. All maps were created using the Cartopy package (https://scitools.org.uk/cartopy)⁴⁷.

## Author contributions

AH wrote the main manuscript text, prepared all figures. All authors reviewed the manuscript.

## Funding

AH is supported by an ESA Living Planet fellowship (PACIFIC project). This work was supported by ISblue project, Interdisciplinary graduate school for the blue planet (ANR-17- EURE-0015) and co-funded by a grant from the French government under the program "Investissements d'Avenir" embedded in France 2030. This work was also supported by LEFE program (LEFE IMAGO/GMMC project VACARM).

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-024-61447-y.

**Correspondence** and requests for materials should be addressed to A.H.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024
42. Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S.L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. et al. Climate change 2021: The physical science basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, vol. 2 (2021).

43. Yu, L. & Weller, R. A. Objectively analyzed air-sea heat fluxes for the global ice-free oceans (1981–2005). *Bull. Am. Meteorol. Soc.* **88**(4), 527–540 (2007).

44. Hochet, A., Huck, T. & Colin de Verdière, A. Large-scale baroclinic instability of the mean oceanic circulation: A local approach. *J. Phys. Oceanogr.* **45**(11), 2738–2754. https://doi.org/10.1175/JPO-D-15-0084.1 (2015).

45. Boucharel, J. et al. A surface layer variance heat budget for ENSO. *Geophys. Res. Lett.* **42**(9), 3529–3537. https://doi.org/10.1002/2015GL063843 (2015).

46. Hochet, A. et al. Direct temporal cascade of temperature variance in eddy-permitting simulations of multidecadal variability. *J. Clim.* **33**(21), 9409–9425. https://doi.org/10.1175/JCLI-D-19-0921.1 (2020).

47. Met Office: Cartopy: A Cartographic Python Library with a Matplotlib Interface. Exeter, Devon (2010–2015). https://scitools.org.uk/cartopy

## Acknowledgements

The datasets analysed during the current study are available in the ECCO repository, https://web.corral.tacc.utexas.edu/OceanProjects/ECCO/ECCOv4/Release3/. ECCO v4r3 salinity and temperature budget computations were greatly facilitated by the existence of the ecco_v4 py Python library, and by several tutorials which can be found at: https://ecco-v4-python-tutorial.readthedocs.io/. We warmly acknowledge the authors of this work. All maps were created using the Cartopy package (https://scitools.org.uk/cartopy)⁴⁷.

## Author contributions

AH wrote the main manuscript text, prepared all figures. All authors reviewed the manuscript.

## Funding

AH is supported by an ESA Living Planet fellowship (PACIFIC project). This work was supported by ISblue project, Interdisciplinary graduate school for the blue planet (ANR-17- EURE-0015) and co-funded by a grant from the French government under the program "Investissements d'Avenir" embedded in France 2030. This work was also supported by LEFE program (LEFE IMAGO/GMMC project VACARM).

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-024-61447-y.

**Correspondence** and requests for materials should be addressed to A.H.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024