OS - Regionalizing the sea-level budget with machine learning techniques

OS - Regionalizing the sea-level budget with machine learning techniques
Articles
Volume 19, issue 1
Article
Assets
Peer review
Metrics
Related articles
Articles
Volume 19, issue 1
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.
Articles
Volume 19, issue 1
Article
Assets
Peer review
Metrics
Related articles
Research article
Highlight paper
16 Jan 2023
Research article |
Highlight paper
16 Jan 2023
Regionalizing the sea-level budget with machine learning techniques
Regionalizing the sea-level budget with machine learning techniques
Regionalizing the sea-level budget with machine learning techniques
Carolina M. L. Camargo et al.
Carolina M. L. Camargo
Riccardo E. M. Riva
Tim H. J. Hermans
Eike M. Schütt
Marta Marcos
Ismael Hernandez-Carrasco
and
Aimée B. A. Slangen
Carolina M. L. Camargo
CORRESPONDING AUTHOR
carolina.camargo@nioz.nl
Department of Estuarine and Delta Systems, NIOZ Royal Netherlands Institute for Sea Research, Yerseke, The Netherlands
Department of Geoscience and Remote Sensing, Delft University of Technology, Delft, The Netherlands
Riccardo E. M. Riva
Department of Geoscience and Remote Sensing, Delft University of Technology, Delft, The Netherlands
Tim H. J. Hermans
Department of Estuarine and Delta Systems, NIOZ Royal Netherlands Institute for Sea Research, Yerseke, The Netherlands
Department of Geoscience and Remote Sensing, Delft University of Technology, Delft, The Netherlands
Eike M. Schütt
Department of Geography, Kiel University, Kiel, Germany
Marta Marcos
Mediterranean Institute for Advanced
Studies (IMEDEA), Spanish National Research Council – University of the Balearic Islands (CSIC-UIB), Esporles, Spain
Ismael Hernandez-Carrasco
Mediterranean Institute for Advanced
Studies (IMEDEA), Spanish National Research Council – University of the Balearic Islands (CSIC-UIB), Esporles, Spain
Aimée B. A. Slangen
Department of Estuarine and Delta Systems, NIOZ Royal Netherlands Institute for Sea Research, Yerseke, The Netherlands
Abstract
Attribution of sea-level change to its different drivers is typically done using a sea-level budget approach. While the global mean sea-level budget is considered closed, closing the budget on a finer spatial scale is more complicated due to, for instance, limitations in our observational system and the spatial processes contributing to regional sea-level change.
Consequently, the regional budget has been mainly analysed on a basin-wide scale.
Here we investigate the sea-level budget at sub-basin scales, using two machine learning techniques to extract domains of coherent sea-level variability: a neural network approach (self-organizing map, SOM) and a network detection approach (
-MAPS).
The extracted domains provide more spatial detail within the ocean basins and indicate how sea-level variability is connected among different regions.
Using these domains we can close, within 1
uncertainty, the sub-basin regional sea-level budget from 1993–2016 in 100 % and 76 % of the SOM and
-MAPS regions, respectively.
Steric variations dominate the temporal sea-level variability and determine a significant part of the total regional change.
Sea-level change due to mass exchange between ocean and land has a relatively homogeneous contribution to all regions. In highly dynamic regions (e.g. the Gulf Stream region) the dynamic mass redistribution is significant.
Regions where the budget cannot be closed highlight processes that are affecting sea level but are not well captured by the observations, such as the influence of western boundary currents.
The use of the budget approach in combination with machine learning techniques leads to new insights into regional sea-level variability and its drivers.
Editorial statement
Closing the sea level budget (between observed rise and the sum of its causes) has been a challenge, is an ongoing effort and has primarily concerned the global mean. Here, the authors use machine learning to identify sub-areas with similar trends to close the sea level budget on a regional level, with much reduced errors compared with 1-degree grid points.
Closing the sea level budget (between observed rise and the sum of its causes) has been a...
Download & links
Article (PDF, 13021 KB)
Download & links
Article
(13021 KB)
Full-text XML
BibTeX
EndNote
Share
How to cite.
Camargo, C. M. L., Riva, R. E. M., Hermans, T. H. J., Schütt, E. M., Marcos, M., Hernandez-Carrasco, I., and Slangen, A. B. A.: Regionalizing the sea-level budget with machine learning techniques, Ocean Sci., 19, 17–41, https://doi.org/10.5194/os-19-17-2023, 2023.
Received: 02 Sep 2022
Discussion started: 13 Sep 2022
Revised: 01 Dec 2022
Accepted: 08 Dec 2022
Published: 16 Jan 2023
The sea-level budget
Sea-level change will be one of the major challenges of the coming centuries for coastal communities worldwide
Fox-Kemper et al.
2021
. Global mean sea-level change has been rising at a rate of 1.6 mm yr
−1
since 1900 and 3.3 mm yr
−1
since 1993
Frederikse et al.
2020
However, sea level does not change uniformly: it displays strong spatial and temporal variations
Hamlington et al.
2020
Ocean dynamics, land ice mass changes and associated gravitational effects, vertical land movement, and the inverse barometer effect are some of the processes responsible for these regional differences
(e.g.
Stammer et al.
2013
Slangen et al.
2017
Understanding the regional variability of the processes driving sea-level change is critical for improving our understanding of its causes, constraining sea-level projections, and better preparing for the impacts of climate change.
The attribution of sea-level change to its different drivers is typically done using a sea-level budget approach
Chambers et al.
2017
WCRP Global Sea Level Budget Group
2018
. For 1993–2018, about one-third of the observed rate of global mean change can be attributed to thermal expansion of the oceans, while the rest is due to the effect of water and ice mass exchanges between land and ocean
Frederikse et al.
2020
. Since the observed rate of sea-level change matches, within uncertainties, with the sum of the contributions of the various sources, the global mean sea-level budget for the period 1993–2018 is considered to be closed
WCRP Global Sea Level Budget Group
2018
Frederikse et al.
2020
Chen et al.
2020
Barnoud et al.
2021
However, locally attributing the drivers of sea-level change for this same period still leads to large differences between the total measured change and the sum of the contributions
(e.g.
Slangen et al.
2014
Royston et al.
2020
This is partly due to the spatial resolution of the current observational systems of the sea-level budget components and of the processes in question, which still limits the closure of the budget on a local spatial scale, for instance on a 1
resolution
Royston et al.
2020
Consequently, the regional sea-level budget has mainly been analysed on a basin-wide scale
(e.g.
Purkey et al.
2014
Frederikse et al.
2018
2020
Royston et al.
2020
and has not been closed on sub-basin scales consistently for the entire world.
The sea-level budget has also been analysed for individual coastline stretches characterized by coherent variability
Rietbroek et al.
2016
Frederikse et al.
2016
2017
Dangendorf et al.
2021
and at individual tide gauges
Wang et al.
2021
The basin-scale sea-level features extracted by
Thompson and Merrifield
2014
have been frequently used in regional sea-level budget studies
Purkey et al.
2014
Frederikse et al.
2018
2020
Royston et al.
2020
Although these publications have made significant advances in understanding the regional sea-level change, the basin scale is still too large to really understand the causes of local variations.
In this paper, we argue that understanding the spatial structure of contemporary sea-level change is a key point to move towards a budget with finer spatial resolution.
By identifying smaller physically coherent regions, some of the effects of small-scale variability can be removed, allowing for closure of the budget at a sub-basin scale.
Machine learning techniques, such as complex and neural networks, can be used to identify such spatial structures, determining the ideal resolution and regions of common sea-level variability and change.
While machine learning methods have widely been used in oceanography
(e.g.
Richardson et al.
2003
Liu et al.
2006
Hernández-Carrasco and Orfila
2018
Sonnewald et al.
2019
Falasca et al.
2019
2020
Novi et al.
2021
, only a few examples analysing sea-surface height can be found
(e.g.
Liu et al.
2016
Ma et al.
2016
Sonnewald et al.
2018
Here, we apply two unsupervised machine learning techniques – self-organizing map (SOM) and
-MAPS – to extract coherent spatial features (domains) in sea-level change observations.
In this study we use the extracted domains to analyse the sea-level budget on a sub-basin scale during the satellite altimetry period (1993–2016) by using state-of-the-art estimates of sea-level change and its components.
We limit our analysis to 2016 because of the temporal span of the hydrological models used to obtain the land water storage contribution to sea-level change. Additionally, instrumental problems (e.g. in Argo salinity data and satellite drifts) have raised questions about the performance and closure of the global mean sea-level budget after 2016
Chen et al.
2020
Barnoud et al.
2021
Cazenave and Moreira
2022
We hypothesize that by investigating the budget in covariant and physically coherent regions, we can resolve the discrepancies (i.e. close the budget) that appear in an increased-resolution sea-level budget (e.g.
1×1
).
Data and methods
In this section we introduce the data sets used for each of the different components of the sea-level budget (Sect.
2.1
). We also describe the trend and budget analysis (Sect.
2.2
) and introduce the machine learning techniques used to extract coherent regions (domains) of sea-level variability and change (Sect.
2.3
).
2.1
The components of the regional sea-level budget
For the budget, we compare the total observed sea-level change
total
to the sum of the drivers of sea-level change
drivers
(1)
total
drivers
where
stands for the rate of sea-level change.
Total sea-level change (
total
) can be measured by tide gauges and satellite altimeters. Satellite altimeters measure geocentric or absolute change (
geo(sat)
), that is, the sea surface height in relation to the reference ellipsoid
Gregory et al.
2019
. On the other hand, tide gauges measure sea surface height in reference to a terrestrial landmark (
rel(TG)
), registering the relative sea-level change. The latter is affected by vertical land motion (VLM) due to, for instance, land subsidence and tectonics
Wöppelmann and Marcos
2015
, while geocentric sea level cannot differentiate if the change is from the solid Earth or the ocean.
The relationship between geocentric and relative sea-level change is
(2)
total
geo(sat)
rel(TG)
VLM
From here onward, when we use
total
, we are referring to the geocentric sea-level change derived from satellite altimetry (Fig.
a). We use multi-mission gridded level-4 data from four distribution centres: CMEMS
CMEMS
2022
, JPL MEaSUREs
Zlotnicki et al.
2019
, SLcci
SLcci
2022
, and CSIRO
CSIRO
2022
. All of these products use the same reference ellipsoid model (GRS80/WGS) and have a monthly temporal resolution, except for JPL MEaSUREs time series, which provides sea surface height data every 5 d and was averaged into monthly means. All data are regridded to a
1×1
map, selected within 66
S to 66
N latitude, and combined into an ensemble mean to avoid systematic errors.
We apply a glacial isostatic adjustment (GIA) correction to the altimetry data from ICE-6G VM5a
Argus et al.
2014
Peltier et al.
2015
by removing the rate of change of the geoid (i.e. Drad
Dsea) from the trends.
Sea-level change expresses changes in the volume of the ocean. These can be caused by changes in the ocean density, mass, or area. Density-driven changes, known as steric sea-level changes, are caused by variations in the ocean temperature and salinity
Gill and Niller
1973
MacIntosh et al.
2017
. All sea-level variations not driven by density changes are known as manometric sea-level change
Gregory et al.
2019
Thus, Eq. (
) can be rewritten to
(3)
total
drivers
SSL
MAN
where
SSL
and
MAN
refer to steric and manometric sea-level change, respectively.
For steric sea-level change (
SSL
, Fig.
c), we use the estimates of
Camargo et al.
2020
, which are based on 15 different ocean temperature and salinity data sets down to 2000 m depth, using Argo floats
Roemmich and Gilson
2009
Gaillard et al.
2016
Li et al.
2017
Lu et al.
2019
, multiple in-situ observations
Ishii and Kimoto
2009
Guinehut et al.
2012
Cabanes et al.
2013
Good et al.
2013
Gaillard et al.
2016
Ishii et al.
2017
Cheng et al.
2019
Szekely et al.
2019
, and ocean reanalyses
Blockley et al.
2014
Maclachlan et al.
2015
Storto and Masina
2016
Garric and Parent
2017
Carton et al.
2018
Zuo et al.
2019
. We complement these data with the deep-ocean steric estimate of
Purkey et al.
2019
, updated from
Purkey and Johnson
2010
Manometric sea-level change (
MAN
), also referred to as the bottom pressure term
Gregory et al.
2019
, can be further divided into (i)
GRD
, i.e. sea-level change due to the gravitational, rotational, and viscoelastic deformation (GRD) response of the Earth to water and ice mass exchanges between land and ocean, and (ii)
DSL
, i.e. the dynamic redistribution of ocean mass due to ocean circulation, atmosphere, and ocean bottom pressure changes as a result of the steric change of the oceans
Landerer et al.
2007
as follows:
(4)
MAN
GRD
DSL
The GRD component (
GRD
, Fig.
d) reflects how the mass loss of continental ice stored in glaciers and ice sheets and variations in land water storage affect sea level. The GRD effect can be split between responses due to contemporary changes and due to the response of the Earth to the last ice age, known as the post-glacial rebound or GIA. The integrated response of the GRD effect over the oceans, i.e. the global mean, is known as barystatic sea-level change (
BSL
Gregory et al.
2019
). For the GRD component, we use the estimates from
Camargo et al.
2022
, which includes the geocentric sea level response to changes in the Antarctic and Greenland ice sheets, glaciers, and terrestrial water storage. These are based on a suite of different estimates of land mass change and computed solving the sea-level equation following
Farrell and Clark
1976
and
Slangen et al.
2014
The dynamic component (
DSL
, Fig.
e) refers to mass changes driven by bottom pressure changes, that is, the redistribution of mass that was already in the oceans. Note that, by our definition, the dynamic sea-level change (
DSL
) is part of the ocean dynamic change
Gregory et al.
2019
the latter also including the effect of local steric anomalies (
SSL
). That is, the dynamic term here is the residual of the sterodynamic sea-level change with the steric contribution removed
Gregory et al.
2019
DSL
is computed from the sea-surface height of five ocean reanalyses (Table
) by first removing the time-varying global mean from the sea-surface height and then by removing the local steric anomaly. This procedure is done in each ocean reanalysis individually, and we then combine the five estimates into an ensemble. We acknowledge that this method introduces some circularity to the budget analysis: the reanalysis, used to obtain
DSL
, assimilates satellite sea-surface height, and in the budget analysis we compare this estimate with satellite sea-surface height (
total
). Compared with the
estimated from Gravity Recovery and Climate Experiment Satellite
(GRACE,
Tapley et al.
2004
DSL
sea-level trends from 2005 to 2015 agree on large-scale patterns and the magnitude of dynamic changes (Fig.
A1
in the Appendix). Note that our budget components do not incorporate GRACE mass changes over the oceans, hence it is an independent estimate for validation. More detail on the estimation and validation of
DSL
is given in Appendix
Finally, Eq. (
) can be rewritten as follows:
(5)
total
SSL
GRD
DSL
such that the total observed sea-level change (Fig.
a) can be compared with the sum of the components (Fig.
b). The ensemble mean of each term of Eq. (
), used throughout this paper for the sea-level budget analysis, is shown in Fig.
, where
total
is the geocentric sea-level change from satellite altimetry corrected for the inverted barometer and GIA (
GIA
) effects;
SSL
is the full-depth steric sea-level change;
GRD
is the contemporary ocean mass redistribution due to the land–ocean mass exchange already corrected for
GIA
effects; and
DSL
is the mass redistribution due to purely ocean dynamics. A summary of the budget components and data set sources is given in Table
. Note that all of the used data sets have been homogenized to a monthly temporal resolution and a
spatial resolution.
CMEMS
2022
Zlotnicki et al.
2019
SLcci
2022
CSIRO
2022
Camargo et al.
2020
Purkey and Johnson
2010
Camargo et al.
2022
Carton et al.
2018
Storto and Masina
2016
Garric and Parent
2017
Blockley et al.
2014
Maclachlan et al.
2015
Zuo et al.
2019
Table 1
Summary of the sea-level budget components and data sources used in this paper.
Download Print Version
Download XLSX
2.2
Computing trends and uncertainties
Our sea-level budget includes the comparison of sea-level time series, trends, and associated uncertainties.
We assume that sea-level trends are the sum of a deterministic model (including annual and semi-annual signals) and stochastic noise (temporal uncertainty).
We use the software Hector
Bos et al.
2013
to compute the trends and the associated 1
uncertainty for each of the budget components.
Following
Bos et al.
2014
Royston et al.
2018
, and
Camargo et al.
2020
2022
, we test eight different noise models to describe the auto-correlation between the residuals of the regression. Using the Akaike and Bayesian information criteria
Akaike
1974
Schwarz
1978
, we select the best-performing noise model at each grid cell. More information on the noise-model analysis can be found in
Camargo et al.
2020
2022
. For the GRD component, in addition to the temporal uncertainties, we also consider the spatial, structural, and intrinsic uncertainties
Camargo et al.
2022
. Note that, unlike for the identification of the domains (Sec.
2.3
), the time series used to estimate trends and uncertainties include seasonality and global mean trends.
We assume independence of the terms and sum the trends linearly and uncertainties in quadrature. For each sea-level domain we take the area-weighted spatial average of the time series, trend, and uncertainties. Performance of the sea-level budget is evaluated by (i) the magnitude of the residual, (ii) the Pearson's correlation coefficient (
) between the altimetry time series and the budget components, and (iii) the normalized root-mean-squared error (nRMSE). The nRMSE measures the distance between the true value, in this case altimetry, and the modelled value, in this case the sum of the budget components. Contrary to
, nRMSE closer to 0 indicates better performance.
Figure 1
Regional sea-level trends for 1993–2016 (mm yr
−1
) for
(a)
altimetry,
(b)
sum of sea-level components,
(c)
full-depth steric,
(d)
GRD effect, and
(e)
dynamic sea-level change. The black contour line indicates global mean sea-level change.
2.3
Clustering techniques
To answer our research questions, we must first identify regions with similar sea-level variability. To do so, we use two different machine learning pattern detection algorithms, one based on an neural network approach, self-organizing map (SOM), and one based on a deep network detection method,
-MAPS.
The methodological differences in these two techniques leads to different patterns of sea-level change in terms of geographical location, region size, and ocean coverage. Hence, by using both methods, we can (i) find prevailing sea-level modes, (ii) compare the patterns and sea-level budget for the different methods, and (iii) balance out the advantages and disadvantages of using a single method.
Both methods are used to reduce the dimensionality of the data, transforming high-dimensional input data into low-dimensional features
Liu et al.
2006
Falasca et al.
2020
For both clustering techniques we use
monthly satellite altimetry time series
CMEMS
2022
from 1993–2019 as input. Note we use a longer time series than the ones for the budget analysis, as longer time series can resolve the temporal variability better. However, additional tests (not shown) showed that the clustering is not strongly affected by the extra 3 years of data.
We pre-process the input data by removing the global mean trend and seasonality and by applying a spatial Gaussian filter of 300 km half-width to remove small-scale variability. Note that after the domain identification for the budget analysis, global mean trend, seasonality, and small-scale variability are included in the time series. Smaller seas, such as the Mediterranean, Baltic, Black, and Caspian seas have been removed from the data prior to the clustering.
2.3.1
Self-organizing map (SOM)
SOM
Kohonen
1982
is a feature extraction and classification method based on an unsupervised neural network
Liu et al.
2006
, which was demonstrated to be more powerful than conventional feature extraction methods
(e.g.
Liu and Weisberg
2005
. The ability of SOM to extract patterns of sea level variability from satellite altimetry data has been shown in previous works
(e.g.
Hardman-Mountford et al.
2003
Liu et al.
2008
2016
Iskandar
2009
Weisberg and Liu
2017
Nickerson et al.
2022
. To analyse sea level data, SOM can be applied either in the spatial domain, focusing on the characteristic spatial patterns, or in the time domain, focusing on the characteristic time series
Liu et al.
2016
. The latter results in regionalizing the sea-level variability and is pursued here to analyse global sea level data. We use the MATLAB SOM toolbox
Vesanto et al.
2000
and follow
Liu et al.
2006
and
Hernández-Carrasco and Orfila
2018
to choose the parameters. We apply the SOM algorithm in the time domain in order to extract the spatial patterns, herein referred to as domains, based on coherent temporal sea-level variability.
Before initializing the SOM, the 3D input data (time, lat, long) is concatenated to 2D
(time, lat times long;
Richardson et al.
2003
Liu et al.
2016
and normalized to have unit variance.
The network is initialized linearly, based on the first two principal components of the time series, and trained in a batch mode; that is, at each step of the training process, all input data vectors are simultaneously used to update the network. Training is performed over 10 iterations, which is necessary to stabilize and converge the network, while avoiding overfitting of the SOM
Liu et al.
2006
. We use the “Epanechikov function” as a neighbourhood function, which returns the most accurate SOM patterns, a hexagonal lattice, and a neighbourhood radius (determining the radius of cells that are updated during the training process) of two cells at the beginning, decreasing linearly to 1 during the training process.
We tested different SOM parameters and verified that this combination gave the smallest quantification errors by computing the averaged Eulerian distance between each data input vector and the best matching unit (BMU). SOM domains do not need to be geographically contiguous, that is, different non-connected regions can be assigned to a single domain. Initially, the strong sea-level variability of the equatorial Pacific Ocean dominated the clustering, hindering pattern identification in the Atlantic Ocean (Fig.
A8
). To overcome this issue we perform the clustering analysis on the Atlantic Ocean and Indo-Pacific Ocean basins separately. We select a map size of
3×3
neurons (i.e. neural network nodes) in each basin, leading to a total of 18 domains. Using different map sizes (e.g. Fig.
A8
) led to more “patchy” results, hence we used a map size of
3×3
neurons as a compromise between the amount of detail and the interpretability of the domains.
2.3.2
-MAPS
-MAPS
Fountalis et al.
2018
is a complex network methodology that reduces the spatiotemporal dimensionality of a field by identifying regions (domains) with similar dynamics and their connectivity
Bracco et al.
2018
Falasca et al.
2020
. Here we focus only on the domains identification (dimensionality reduction) function of the
-MAPS method.
-MAPS domains are spatially continuous (i.e. grid cells need to be physically connected to be clustered in the same domain) and are potentially overlapping regions that have a highly correlated temporal activity
Falasca et al.
2019
Formally, each input grid cell is associated with a time series, including the
nearest neighbours, based on the haversine distance (angular distance between two points on a sphere). The local homogeneity, defined as the average Pearson cross-correlation between a grid cell and its
neighbours, is computed and tested against a threshold value
. If the local homogeneity is greater than
, with a statistical significance level of 0.1, then the grid cell is considered a core, which then is expanded to identity domains
Fountalis et al.
2018
Falasca et al.
2019
Novi et al.
2021
. Each domain expands to adjacent cells as long as the local homogeneity continues to be higher than
To choose the optimal neighbourhood size
, we follow a heuristic approach, testing
values from 4 to 25 following
Falasca et al.
2019
. As in
-MAPS not every grid point needs to belong to a domain (in contrast to SOM), we then choose the
value taking into account the amount of unclustered cells (i.e. the one with most of the ocean belonging to domains). We also use the normalized mutual information (NMI) matrix
Falasca et al.
2019
to identify the
value with high NMI for it and its neighbouring
values, meaning that the results are less sensitive to the chosen
value. These parameters led to the use of
=5
Identifying domains of sea-level variability
Both clustering methods successfully reduce the dimensionality of the input data, despite the higher number of domains identified by
-MAPS (Fig.
).
SOM identified 18 coherent domains, with a domain area varying from 3.84 to 34.51 million square kilometres, and an average and total size of 17.61 and 316.90 million square kilometres, respectively.
-MAPS identified 92 coherent domains, with a domain area varying from 0.03 to 24.15 million square kilometres, with average and total size of 2.53 and 242.01 million square kilometres, respectively.
Despite the methodological differences, we find that prominent sea-level features are clustered in a similar way by SOM and
-MAPS (Fig.
). Some of the patterns identified can be linked with known oceanic patterns, as we will discuss below. However, we note that covariability does not imply a common forcing and that some patterns may be statistically separated or grouped without a clear physical reason.
It is also important to note that these clustering methods do not account for auto-correlation in time, meaning the time lag in the progression of a signal across the ocean basin. Since we use monthly data, signals that propagate faster than a month (typically barotropic) will be more clearly correlated in our clustering. On the other hand, slower propagating signals, such as the first baroclinic mode, will lose correlation in space and will not be represented in the identified domains.
The central Pacific domain, where the variability is dominated by El Niño–Southern Oscillation (ENSO) events, covers a similar region in both methods. The “ENSO tongue”, starting from the coast of Peru and Ecuador and spreading west until the central Pacific, is identified by both methods (SOM domain 12 (pink) and
-MAPS domain 45 (light green)).
The western tropical Pacific Ocean (WTPO), influenced by ENSO and the Pacific decadal oscillation (PDO), is also identified as a single domain by both methods (SOM domain 16 (light green) and
-MAPS domain 89 (light brown)). The WTPO domains matches with the region of significant spatial correlation between steric and coastal sea level found by
Dangendorf et al.
2021
for western Australia. In the SOM clustering, the WTPO domain incorporates the Leeuwin Current
(western Australia,
Pattiaratchi and Siji
2020
in the Indian Ocean, which is affected by waves travelling through the tropical Australasian seas
Feng et al.
2004
. While this connection is not captured by
-MAPS, the coherence along the western coast of Australia is featured in a single domain (
-MAPS domain 92, light pink).
The Kuroshio Extension region is also identified in both methods (SOM domain 10 (brown) and
-MAPS domain 88 (brown)), reflecting how strong boundary currents influence the sea-level variability. Another example is the North Atlantic, which has similar clustering in both methods, especially in the domain south of Greenland (SOM domain 9 (light purple) and
-MAPS domain 33 (purple)), which is marked by decadal-scale sea-level change reflecting the strength and shape of the wind-driven subpolar gyre and the Atlantic meridional overturning circulation
Chafik et al.
2019
. Within these domains, density anomalies are known to flow southward from the Labrador Sea into the subpolar gyre through coastally trapped waves
Dangendorf et al.
2021
. Another region identified in both methods is the northwestern European Shelf (SOM domain 8 (purple) and
-MAPS domain 66 (grey)), which is part of a domain that extends along the whole western European coast, continuing down to the Canary islands and well into the Atlantic. This connection could be related to the hypothesis that coastally trapped waves and longshore winds cause a coherent region of sea-level variability from around the latitude of the Canary Islands up to the Norwegian Sea
Calafat et al.
2012
Chafik et al.
2019
Hughes et al.
2019
Hermans et al.
2020
Dangendorf et al.
2021
. These features are in a separate
-MAPS domain (53, green) to the northwestern European Shelf. It is important to note that coherent features smaller than 300 km are not captured in the domains because of the spatial filtering applied before the clustering analysis.
As SOM domains do not need to be contiguous, possible pseudo-teleconnections between different ocean regions (within the Atlantic Ocean and Indo-Pacific Ocean basins) come out of the analysis.
For example, areas adjacent to the ENSO tongue domain (both to north and south) are clustered together in domain 18 (light blue) or in domain 15 (moss green), indicating how the ENSO signal is propagated through the Pacific, possibly through coastally trapped waves
Hughes et al.
2019
in the coastal domains (15) or via atmospheric teleconnections.
However, not every region classified into the same SOM domain results from a clear connection.
For example, SOM domain 17 (blue) groups the ocean adjacent to South Africa, the region below the Kuroshio Extension (offshore of Taiwan) and a region south of Australia and New Zealand. Another example is SOM domain 7 (salmon-pink), which implies a connection between the Atlantic Caribbean Sea and the western part of the South Atlantic Gyre (capturing parts of the Brazil Current).
These regions have been classified together because they have a similar behaviour in terms of sea-level variability but probably different forcing.
Further investigation, with ocean currents, ocean–atmospheric oscillations, and ocean waves are necessary to explore and quantify the physical connection behind these patterns.
Unlike SOM, every
-MAPS domain is assigned a unique number and not every pixel needs to be clustered (Fig.
a, white regions). Consequently, this method yields a larger number of domains with smaller size while avoiding pseudo-teleconnections.
The dominant sea-level modes are clear on
-MAPS clustering, reflecting the influence, for example, of ENSO and western boundary currents on sea level. For example, the entire Caribbean (domain 87, brown) and Gulf of Mexico (domain 82, red) is in a single domain, highlighting the similarity in that region. The same goes for the Equatorial Atlantic (domain 86, light purple), the ENSO region (domains 45, 89, and 62, shown as light green, light brown and light brown, respectively), and the Kuroshio current (domain 88, brown).
Figure 2
Domains of coherent sea-level variability:
(a)
-MAPS method (92 domains) and
(b)
self-organizing map (SOM) method (18 domains). Numbers indicate the domain code, and domain names are given in the main text and in Table
A1
in the Appendix.
-MAPS domains with codes in magenta indicate selected domains for Fig.
. For visibility, small domains have not been labelled. Given the large number of domains,
-MAPS has a repeating colour pallet, but since
-MAPS domains need to be continuous, repeated colours do not indicate the same domain. White regions in
-MAPS indicate incoherent regions, which were not incorporated in any domain.
As shown in
Royston et al.
2020
, the components of the sea-level budget have a similar spectral power to the total observed sea-surface height of altimetry between wavelengths of approximately 3000 and 10 000 km.
The clustering techniques applied here not only reduce the dimensionality of the data but also average out sea-level variability in regions of coherent variability, being ideal for a regional budget analysis (next section).
The regional sea-level budget on different spatial scales
4.1
Sea-level trend budget closure
We investigate the trends of the sea-level budget on different spatial scales, from a finer (
1×1
) to coarser scale (
-MAPS and SOM domains; Fig.
). The residuals (i.e. the difference between the total sea-level change and the sum of the components) decrease towards a coarser spatial scale: for 1
, they range from
−8.2
to 21.1 mm yr
−1
, while for
-MAPS they range from
−1.2
to 3.8 mm yr
−1
and for SOM from 0.1 to 0.7 mm yr
−1
. This shows an improvement in the budget closure (i.e. total and sum of components agree within uncertainties) by using the pattern detection algorithms. The budget closes in all 18 SOM domains (100 % of SOM ocean area), in 70 out of 92 of the
-MAPS domains (94 % of
-MAPS ocean area, 229.9 million square kilometres), and in 72 % of the grid cells in the 1
budget (75 % of the ocean area) (Fig.
). There is a clear relation between spatial scale of the region considered for the trend, and the residuals of the budget (see also Fig.
A6
). The good closure in the 1
budget is likely an artefact of the large uncertainties of the observations, which on a local scale can be up to 18.9 mm yr
−1
(see Fig.
A3
). This is in line with
Royston et al.
2020
, who found that local biases of steric estimates combined with the resolution limitation of GRACE observations over the oceans hinder the budget closure at 1
resolution.
When the regional domains based on SOM and
-MAPS are considered, the uncertainties show a 5-fold reduction compared to the 1
resolution, reaching up to 3.6 mm yr
−1
and an average value of 1.6 mm yr
−1
(Fig.
A3
), while the budget still closes.
Consequently, there is a better match between the total observed rate of sea-level change with the sum of the components for the clustered regions (scatter points in Fig.
, right column), with a reduction in the spread of the scatter points and a move closer to the
1:1
line (dashed black line) for the coarser resolutions.
The dashed pink lines in Fig.
indicate the half-width of the 95 % confidence interval of the uncertainty of the residuals, showing a slightly larger width for 1
and a smaller one for SOM and
-MAPS. Even when the component uncertainties (grey error bars) are considered, the scattered values are mostly within the width of the 95 % confidence interval for the SOM domains, confirming the improvement of the budget for this case.
There is a strong linear correlation between the total and the sum of the drivers, with Pearson's
varying from 0.81 for the 1
budget and
-MAPS to 0.98 for SOM. The RMSE also decreases for the coarser scales, from 1.01 mm yr
−1
for the 1
budget to 0.47 mm yr
−1
for the SOM domains.
The altimetry trends are generally larger than the sum of the sea-level change drivers, as indicated by the positive residuals and scatter points above the
1:1
line on Fig.
. This is true for more than half of the
-MAPS domains and for all SOM domains except one: SOM domain 9 (south of Greenland) is marked by a negative residual, that is, the sum of the drivers is larger than the observed altimetry trend. Several
-MAPS domains, such as southwest of Australia (domains 92 and 67), the southeastern Pacific (domains 37 and 54), the Gulf Current (domain 82), and the Brazil–Falklands confluence zone (domains 80 and 69), also have a negative residual.
This might indicate a larger temporal variability or regime shifts in this region or could be due to the ocean dynamics contribution, such as the effect of the subpolar gyre around the south of Greenland
Chafik et al.
2019
, as we will see in the next section (Sect.
4.2
).
Figure 3
Sea-level budget residuals (maps
a, c, e
) and comparison between total sea-level change (
-axis) and sum of components (
axis) (scatter plots
b, d, f
) for 1
(a, b)
-MAPS domains
(c, d)
, and SOM domains
(e, f)
. Grey lines indicate the uncertainties (1
) of the components. In the scatter plots every point indicates one region (grid in case of 1
), dashed pink lines indicate the half-width of the 95 % confidence interval of the residuals uncertainty, and grey error bars indicate the component uncertainty. The use of
indicates that coefficient is statistically significant (
value
<0.01
).
4.2
Explaining the sea-level budget contributions
In this section, we investigate which components dominate the trend and temporal variability in each of the different domains. For comparison and discussion purposes, we choose 18
-MAPS domains (magenta numbers, Fig.
b) located close to the 18 SOM domains. Trends for all
-MAPS domains are available online as an interactive map (see the caption of Fig.
).
As shown previously, we find a good match between total observed sea-level change and the sum of components (Fig.
a, b, green stars and purple triangles, respectively) for all SOM and
-MAPS domains.
The largest budget uncertainties, considering both altimetry and the sum of components, is seen in the WTPO domain (SOM 16,
-MAPS 89). These uncertainties may be related to (i) poor performance of standard altimetry products in these shallow regions, (ii) poor Argo float coverage in the region
Kleinherenbrink et al.
2017
, influencing both the steric and dynamic components, and (iii) large internal variability due to ENSO events in this region, which may contribute to large temporal uncertainties in the steric and altimetry components
Kleinherenbrink et al.
2017
Wagner and Böning
2021
. This region is also within the Indian-south Pacific basin
Thompson and Merrifield
2014
, which was the only basin in which the regional budget from 2005–2015 could not be closed
Royston et al.
2020
The GRD component (Fig.
, blue) has a relatively comparable contribution to all regions, contributing about 1.5 mm yr
−1
of sea-level rise. The dynamic and steric components, however, show a strong regionally varying contribution (Fig.
, red and yellow, respectively). For example, for SOM (
-MAPS) domains 10 (88), 13 (92), 14 (67), and 16 (89), more than 50 % of the total trend is due to steric variations. On the other hand, for SOM domains 1 and 18 and
-MAPS domains 39, 45, and 62, the steric trend explains less than 20 %. The dynamic component shows a small contribution for most of SOM domains, and in some domains even shows a negative trend (e.g. SOM domains 11, 12, and 14). An exception is the Gulf Stream domain (SOM 1,
-MAPS 82), where almost half of the total trend is explained by the dynamic component. This dominance of the dynamic component reflects the influence of the strong western boundary current on sea level in this region. The south of Greenland domain (SOM 9,
-MAPS 33) also includes a relatively large dynamic contribution, with a trend of
0.49±0.21
mm yr
−1
, reflecting the influence of the subpolar gyre in this region. The dynamic component also has a significant contribution to other
-MAPS domains, such as domains 24, 69, 39, 66, and 67. Domain 67, located southwest of Australia, shows a large negative dynamic trend, which can be related to the influence of the West Australian Current.
Regarding the temporal evolution (Figs.
c, d and
A5
), both the SOM (solid lines) and the
-MAPS (dashed lines) time series show a similar behaviour. The steric component dominates the temporal sea-level variability, with a good match to the altimetry. The time series of the ENSO tongue domain (SOM 12 and
-MAPS 45, Fig.
d) shows the clear response of sea level to ENSO, with peaks coinciding with strong ENSO events, such as the El Niño events of 1997 and 2015. The prominent contribution of the dynamic component to the total trend in the Gulf Stream domain (SOM 1 and
-MAPS 82) is not reflected in the time series (Fig.
c). Hence, while the dynamic component has a significant impact on the overall change, it does not contribute to the seasonal to interannual sea-level variability. This is true for all other domains (Fig.
A5
), except for SOM (
-MAPS) domain 2 (24) and 18 (85), where we find a better match between the dynamic and altimetry time series.
Figure 4
Sea-level budget trends (mm yr
−1
) for
(a)
SOM and
(b)
-MAPS domains and
(c, d)
time series for two example domains, where solid and dashed lines indicate SOM and
-MAPS time series, respectively. The location of each domain is shown in Fig.
(domain numbers in magenta for
-MAPS). For comparison,
-MAPS domains are matched to the SOM domains, for example SOM domain 12 to
-MAPS domain 45. A bar plot for all other
-MAPS domains can be found in Fig.
A4
. Error bars indicate the 1
uncertainty of the trend. Time series for all SOM domains and for the 18
-MAPS domains in
(b)
are shown in Fig.
A5
. An interactive budget map is available at
(last access: 9 January 2023) for both SOM and
-MAPS.
4.3
Sea-level budget performance
Here, we investigate the closure of the budget considering (i) the components included in the budget, (ii) the size of the domains and the clustering method, and (iii) the data sets used for each component.
To illustrate the performance of the budget considering the domains used and the components included in the budget, we show how the Pearson's correlation coefficient (
) and the normalized root-mean-squared error (nRMSE) change when these factors vary (Fig.
).
This figure firstly shows that the budget closure improves when more components are included in the budget.
While we get a poorer performance when only considering the dynamic or the GRD component, the budget with only steric components already performs relatively well. The improved correlation and lower RMSE with the steric component is not surprising as the seasonal cycle is predominantly steric.
The budget performance is enhanced by the addition of the dynamic and GRD components, shown by the narrowing of the box-and-whiskers plot.
The figure also shows an improvement in the budget closure for
-MAPS and SOM domains, in relation to the 1
resolution, regardless of the budget combination.
There are two possible reasons why a coarser spatial resolution leads to decreasing uncertainties and a better budget closure:
(i) the spatial scale of the process itself, as changes in long-term sea level typically occur on a coarser resolution than 1
and/or (ii) there is a mismatch in the exact location between the sum of the components and altimetry observations on a finer spatial scale, resulting from the limited resolution of the observations compared to a coarser scale when such mismatches are partially averaged out. Additionally, the averaging of more samples leads to a smaller standard error. However, the measurement errors between altimetry and the sum of components will only compensate each other if they are uncorrelated in the spatial scale being analysed. The relationship between the spatial scale of the domains and the performance of the budget is further confirmed in Fig.
A6
, which shows how the residual of the budget decreases when larger regions are considered.
Note, however, that simply upscaling the resolution of the observations – i.e. considering
2×2
or
5×5
blocks – does not have the same effect on budget performance as the domains derived by machine learning (Fig.
A7
): there is demonstrated added value of considering regions that are physically coherent, rather than artificial blocks, for the budget analysis. That is, spatially averaging over areas of similar variability reduces the unexplained variance of the observations.
Figure 5
The effect on budget closure for different component combinations and spatial resolutions.
(a)
Pearson's correlation coefficient (
) and
(b)
normalized root-mean-squared error (nRMSE) (in mm yr
−1
) between total sea-level (altimetry) and the different components included in the budget (
axis) for the different spatial resolutions (1
in green,
-MAPS in red, and SOM in blue). Boxes represent the quartiles of the distribution, extending from the lower to upper quartile values of the data, with a line at the median, while the whiskers (not error bars) show the full distribution.
When it comes to the data sets of the different sea-level change drivers, sea-level budget studies often use the ensemble mean of several data sets for each component
(e.g.
WCRP Global Sea Level Budget Group
2018
or compute a range of budget combinations by varying the data set for each component to find the combination that returns the best budget closure
(e.g.
Gregory et al.
2013
. The latter approach can result in a budget closure for the wrong reasons
Royston et al.
2020
. On the other hand, while the ensemble mean approach may reduce the systematic biases of using individual data sets
Storto et al.
2017
, it may also hinder the real variability of the process being analysed
Rougier
2016
. Alternatively, the budget can be analysed with the data sets closest to the ensemble means, according to the RMSE analysis, which retains the true variability of an individual data set
Rougier
2016
Royston et al.
2020
All the results presented so far were computed using the ensemble means for each component, considering 15 steric, 5 dynamic, 4 barystatic, and 4 altimetry data sets. Considering all single data sets plus the ensemble of each component we can obtain 2400 possible budget combinations (
16
). To illustrate the dependence of the budget closure on the data set used, we now also discuss the residuals of each SOM domain considering all 2400 possible data set combinations (Fig.
). The residual value shows a large spread for the different budget combinations, ranging from about
−2
to 2 mm yr
−1
, and 33 % of the combinations would result in non-closure of the budget (i.e. the sum of the components does not match with the altimetry values, indicated in red).
The residuals of the ensemble combinations (used throughout this study and indicated by the filled blue squares in Fig.
) are comparable with the residuals of the combinations using the data sets with the smallest RMSE to the ensemble mean (indicated with filled purple triangles in Fig.
).
With the exception of the domains 14 and 16, we see that the ensemble and the RMSE combination have a similar residual value. This indicates that the closure of the budget is not an artefact of the data set choice.
Figure 6
Budget residuals (mm yr
−1
) for all possible data set combinations for every SOM domain (separated by vertical dashed lines). The ensemble mean combination (used for the main analysis) is indicated with squares. The RMSE combination, that is, the budget combination that uses the individual data set with smallest RMSE in relation to the ensemble mean for each component, is indicated with triangles.
Discussion and conclusions
Sea-level budget assessments are important tools for understanding the processes driving sea-level change, for detecting temporal changes in sea-level and its components, for identifying missing contributions to the budget, and for validating and constraining climate models used in sea-level projections
Cazenave and Moreira
2022
. In particular, understanding the processes on a finer spatial scale is essential for local sea-level projections and coastal management planning. In this study, we investigated the regional sea-level budget for 1993–2016 on a global scale.
Regional sea-level budget closure tends to be difficult due to the complex physical processes acting on different spatial scales. To overcome this spatial resolution issue, we applied a neural network approach, SOM, and a deep-network detection method,
-MAPS, to identify domains of coherent sea-level variability (Fig.
). Note, however, that the coherent patterns will be different based on whether total sea surface height or the individual components (steric, dynamic, GRD) are considered. Hence, depending on the purpose of the study, it is important to first remove the unwanted components from total sea-surface height and then perform the clustering. The identified patterns reflect, among other factors, the influence of natural internal climate modes
Han et al.
2017
, such as the ENSO, PDO, and North Atlantic oscillation. This indicates the potential of using machine learning and pattern detection algorithms, such as SOM and
-MAPS, to isolate the effects of natural climate modes from anthropogenic forcing on sea-level change. The domains also suggest how sea-level variability may be transferred between ocean regions. For example, the northwestern European Shelf SOM domain extends south down to the Strait of Gibraltar, possibly reflecting how coastally trapped waves propagate sea-level variability into the North Sea
Calafat et al.
2013
Dangendorf et al.
2014
2021
Hughes et al.
2019
Hermans et al.
2020
. Additionally, highly energetic ocean regions, such as the Kuroshio Current, the Gulf Stream and the Malvinas Confluence Zone, are also extracted as single features, matching the spectrum of sea-level variability in those zones
Hughes and Williams
2010
Compared with the basin regions of
Thompson and Merrifield
2014
, we have identified more and smaller domains, especially in the Southern Hemisphere. This means our domains can provide an additional level of spatial detail compared to ocean basins, while remaining large enough to provide a consistently closing regional sea-level budget. Using the domains identified with SOM and
-MAPS, we presented a regional sea-level budget assessment on an average scale of about
5×10
km
, with the largest regions about
30×10
km
. The performance of the budget improves from finer (1
resolution) to coarser scale (SOM domains), with a residual spread of 0.6 mm yr
−1
for SOM compared to 29.2 mm yr
−1
for 1
resolution. We also showed that the budget closes better when all components (steric, dynamic and GRD) are included, highlighting the importance of including the deep steric and dynamic contributions to regional sea-level change. Despite the large uncertainties at a regional scale (compared to the global mean)
Royston et al.
2020
, we were able to identify dominant drivers in most domains. The
-MAPS regions where the budget cannot be closed highlight processes that are affecting sea level but are not well captured by the observations, such as the influence of western boundary currents and dynamic processes (e.g. the Malvinas Confluence zone). They may also be related to the quality of global data sets in continental shelves and close to the coast or to instrumental noise.
The GRD component has a relatively homogeneous contribution, independent of the domain, in agreement with
Frederikse et al.
2020
. The steric contribution dominates the seasonal and interannual variability and results in the prevailing sea-level trend in most domains, especially for domains in the Southern Hemisphere and equatorial regions. The dynamic component is important in some regions, particularly in the Gulf Stream domain. The domains where the dynamic component plays an important role coincide with the coastal polygons of
Rietbroek et al.
2016
where a large part of the budget could not be explained solely by the sum of steric and land–ocean mass exchange. Hence, our analysis sheds light on the unexplained variance of previous sea-level budget studies. Note that the sea-level analysis in coastal regions is more challenging
Dangendorf et al.
2021
, since some of the dominant coastal ocean dynamics are not properly represented in the global data sets
Liu and Weisberg
2007
Here we showed that pattern detection techniques based on machine learning, such as SOM and
-MAPS, are powerful approaches for identifying and understanding features of global sea-level change and variability. The domains identified in this research highlight that different ocean regions are interconnected, revealing how large-scale circulation controls regional sea level. These domains are not only a good starting point for a regional sea-level budget analysis but also have the potential to separate natural and anthropogenic forcings of sea-level change in a detection and attribution approach, building on previous work
(e.g.
Marcos and Amores
2014
Slangen et al.
2014
2016
. Future work may include multiple linear regressions with climate modes to explore this potential. Additionally, these domains can also be used for coastal sea-level reconstructions (e.g. as
Dangendorf et al.
2021
) and for pattern scaling in sea-level projections
Bilbao et al.
2015
Appendix A:
Dynamic sea-level change estimation and validation
The dynamic redistribution of mass due to ocean circulation and atmospheric redistribution effects is known as dynamic sea-level change (
DSL
Landerer et al.
2007
Gregory et al.
2019
).
DSL
refers to mass changes driven by bottom pressure changes, that is, the redistribution of mass that was already in the oceans and includes mass exchange at any point by mass redistribution by wind stress and by non-linear interaction due to density changes. Note that, by our definition, the dynamic sea-level change (
DSL
) is part of the ocean dynamic sea-level change
Gregory et al.
2019
the latter also including the effect of local steric anomalies (
SSL
).
When the ocean dynamic component (
) is considered together with the global mean steric sea-level change (
SSL
), then it is known as sterodynamic sea-level change (
SDSL
Gregory et al.
2019
Dangendorf et al.
2021
Wang et al.
2021
).
By decomposing the steric component in a global mean (denoted with the overline bar) and local anomaly component (denoted by the prime symbol), we can write the sterodynamic equation as follows:
(A1)
SDSL
SSL
DSL
SSL
SSL
To obtain
, we use the sea-surface height of five ocean reanalysis data sets (SODA
Carton et al.
2018
, C-GLORS
Storto and Masina
2016
, GLORYS
Garric and Parent
2017
, FOAM-GloSea
Blockley et al.
2014
Maclachlan et al.
2015
and ORAS
Zuo et al.
2019
). As ocean reanalyses conserve mass
Griffies and Greatbatch
2012
, the sea-surface height of a reanalysis does not include the GRD component, but it does include the steric effect. We acknowledge that this method introduces some circularity to the budget analysis: the reanalysis, used to obtain
DSL
, assimilates satellite sea-surface height, and in the budget analysis we compare this estimate with satellite sea-surface height (
total
). Following
Wang et al.
2021
, we compute ocean dynamic sea-level change by removing the time-varying global mean from the reanalysis' sea surface height:
(A2)
rea
rea
Since we are interested purely in the dynamic part of
, that is, the dynamic sea-level change (
DSL
), we must remove the steric local anomaly (
SSL
) as follows:
(A3)
DSL
SSL
SSL
SSL
where the steric estimate has been computed with the ocean temperature and salinity of the respective reanalysis. We then compute the ensemble mean of the five dynamic estimates.
Figure A1
Dynamic sea-level change (
DSL
) estimated from
(a)
GRACE (average of JPL and CSR mascons),
(b)
an ensemble of ocean reanalysis, and
(c)
the difference between GRACE and the reanalysis.
To validate our estimate of
DSL
, we compare it with
DSL
estimated from the Gravity Recovery and Climate Experiment Satellite
(GRACE,
Tapley et al.
2004
. GRACE measures total mass changes, which can be used to derive estimates of manometric sea-level change over the oceans, meaning the change in response to both the dynamic ocean mass redistribution (
DSL
) and to mass redistribution due to the land–ocean mass exchange (
GRD
Chambers et al.
2004
Royston et al.
2020
. We use GRACE mass concentrations (mascons) products over the oceans from two different processing centres: RL06 from the Center for Spatial Research
(CSR,
Save et al.
2016
2020
and RL06 v02 from the Jet Propulsion Laboratory
(JPL,
Watkins et al.
2015
Wiese et al.
2019
. In order to obtain the
DSL
, we then remove the GRD patterns obtained for the same data sets by
Camargo et al.
2022
. Note that we use GRACE dynamic sea-level change for validation purposes but not in our budget analysis, as this data set only starts in 2002.
Qualitatively,
DSL
obtained from GRACE (Fig.
A1
a) and from ocean reanalysis (Fig.
A1
b) agree on large-scale patterns and magnitude of dynamic changes despite local differences (Fig.
A1
c).
The main differences are in the region surrounding Indonesia and Japan, related to the signature of the 2004 Sumatra (Indonesia) and 2011 Tohoku (Japan) mega-thrust earthquakes
Chen et al.
2007
Ghobadi-Far et al.
2020
on GRACE observations. To a lesser extent, we also see the effect of the 2010 Maule (Chile) earthquake and tsunami
Ghobadi-Far et al.
2020
. Another strong divergence is seen in the South Atlantic, where the positive trends of GRACE are not represented in the reanalysis, possibly suggesting that a source of dynamic sea-level change is not well parameterized in the reanalysis. Alternatively, this divergence might also be an artefact of the GRACE spherical harmonic solutions and low-degree corrections.
Figure A2
Uncertainties of the regional sea-level trends (Fig. 1) for 1993–2016 (mm yr
−1
) for
(a)
altimetry,
(b)
sum of sea-level components,
(c)
full-depth steric,
(d)
GRD effect, and
(f)
dynamic sea-level change. The black contour line indicates global mean sea-level change.
Figure A3
Distribution histograms of the altimetry (
a, c, e
, green blocks), sum of the components (
a, c, e
, purple blocks), and residuals trend (
b, d, f
, grey) and uncertainty (right columns, pink) for the
1×1
budget
(a, b)
-MAPS domains
(c, d)
, and SOM domains
(e, f)
. The dashed pink lines indicate the 95 % confidence interval of the residuals uncertainty, with the interval width reported in the panel titles, and was used as a reference for the residuals scatters in Fig. 3.
Table A1
Names of SOM and
-MAPS domains.
Download Print Version
Download XLSX
Figure A4
Extension of Fig. 4, showing the trend contribution for each of the
-MAPS domains. An interactive budget map is available at
(last access: 9 January 2023).
Figure A5
Extension of Fig. 4, showing the time series for each of the SOM domains and corresponding
-MAPS domains.
Figure A6
Scatter plot of the budget residuals (i.e. altimetry minus sum of components) against the area of each domain for
-MAPS (red) and SOM (blue). Stars and circles indicate domains in which the sea-level budget is open and closed, respectively. As the domain area increases, the residuals converge towards 0. All the SOM residuals are within
±1
mm yr
−1
, as are 74.2 % of the
-MAPS domains.
Figure A7
Same as Fig. 5 but including sea-level budget considering blocks of
2×2
and
5×5
. There is no clear improvement from the
1×1
budget to the
2×2
and
5×5
, showing the added value of using the
-MAPS and SOM domains.
Figure A8
SOM domains for a neural map of
4×4
(a)
and of
(b)
9×9
, using the entire ocean as input for the clustering.
Note that even with a larger neural map the SOM patterns are still different from the
-MAPS domains (Fig. 2), proving that the differences between the extracted patterns are not just a function of the number of SOM neurons but are due to differences between the two methods. It does, however, leads to fewer regions that are geographically distant being clustered in the same domain.
Code and data availability
Sea-level trends and scripts used for the budget analysis are available at
Camargo et al.
2022
) and
carocamargo
2023
). Interactive maps of the sea-level budget are available at
(last access 9 January 2023,
Camargo
2023
).
Author contributions
CMLC performed the research and drafted the article. CMLC, REMR, and ABAS designed the study.
The self-organizing map method and analysis was applied by CMLC, IHC, and MM. CMLC and EMS performed
-MAPS analysis. All authors contributed to the interpretation of the results and the writing of the manuscript.
Competing interests
The contact author has declared that neither of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Financial support
This research has been supported by the Netherlands Space Office User Support programme (grant no. ALWGO.2017.002)
Review statement
This paper was edited by John M. Huthnance and reviewed by Samantha Royston and three anonymous referees.
References
Akaike, H.: A New Look at the Statistical Model Identification, IEEE
T. Automat. Contr., 19, 716–723,
, 1974.
Argus, D. F., Peltier, W. R., Drummond, R., and Moore, A. W.: The Antarctica
component of postglacial rebound model ICE-6G_C (VM5a) based on GPS
positioning, exposure age dating of ice thicknesses, and relative sea level
histories, Geophys. J. Int., 198, 537–563,
, 2014.
Barnoud, A., Pfeffer, J., Guérou, A., Frery, M. L., Siméon, M.,
Cazenave, A., Chen, J., Llovel, W., Thierry, V., Legeais, J. F., and Ablain,
M.: Contributions of Altimetry and Argo to Non-Closure of the Global Mean
Sea Level Budget Since 2016, Geophys. Res. Lett., 48, 1–10,
, 2021.
Bilbao, R. A. F., Gregory, J. M., and Bouttes, N.: Analysis of the regional
pattern of sea level change due to ocean dynamics and density change for
1993–2099 in observations and CMIP5 AOGCMs, Clim. Dynam., 45,
2647–2666,
, 2015.
Blockley, E. W., Martin, M. J., McLaren, A. J., Ryan, A. G., Waters, J., Lea, D. J., Mirouze, I., Peterson, K. A., Sellar, A., and Storkey, D.: Recent development of the Met Office operational ocean forecasting system: an overview and assessment of the new Global FOAM forecasts, Geosci. Model Dev., 7, 2613–2638,
, 2014.
Bos, M. S., Fernandes, R. M., Williams, S. D., and Bastos, L.: Fast error
analysis of continuous GNSS observations with missing data, J. Geodesy, 87, 351–360,
, 2013.
Bos, M. S., Williams, S. D., Araújo, I. B., and Bastos, L.: The effect
of temporal correlated noise on the sea level rate and acceleration
uncertainty, Geophys. J. Int., 196, 1423–1430,
, 2014.
Bracco, A., Falasca, F., Nenes, A., Fountalis, I., and Dovrolis, C.: Advancing
climate science with knowledge-discovery through data mining, npj Clim. Atmos. Sci., 1, 1–6,
, 2018.
Cabanes, C., Grouazel, A., von Schuckmann, K., Hamon, M., Turpin, V., Coatanoan, C., Paris, F., Guinehut, S., Boone, C., Ferry, N., de Boyer Montégut, C., Carval, T., Reverdin, G., Pouliquen, S., and Le Traon, P.-Y.: The CORA dataset: validation and diagnostics of in-situ ocean temperature and salinity measurements, Ocean Sci., 9, 1–18,
, 2013.
Calafat, F. M., Chambers, D. P., and Tsimplis, M. N.: Mechanisms of decadal
sea level variability in the eastern North Atlantic and the Mediterranean
Sea, J. Geophys. Res.-Oceans, 117, 1–14,
, 2012.
Calafat, F. M., Chambers, D. P., and Tsimplis, M. N.: Inter-annual to decadal
sea-level variability in the coastal zones of the Norwegian and Siberian
Seas: The role of atmospheric forcing, J. Geophys. Res.-Oceans,, 118, 1287–1301,
, 2013.
Camargo, C. M. L.: Interactive maps for the regional sea-level budget,
, last access 9 January 2023.
Camargo, C. M., Riva, R. E., Hermans, T. H., and Slangen, A. B.: Exploring
Sources of Uncertainty in Steric Sea-Level Change Estimates, J. Geophys. Res.-Oceans, 125, 1–18,
, 2020.
Camargo, C. M. L., Riva, R. E. M., Hermans, T. H. J., and Slangen, A. B. A.: Trends and uncertainties of mass-driven sea-level change in the satellite altimetry era, Earth Syst. Dynam., 13, 1351–1375,
, 2022a.
Camargo, C. M. L., Riva, R. E. M., Hermans, T. H. J., Schütt, E. M., Marcos, M., Hernandez-Carrasco, I., and Slangen, A. B. A.: Regional Sea-level Budget from 1993–2016 (1.1), Zenodo, [data set], https://doi.org/10.5281/zenodo.7385561 (last access: 9 January 2023), 2022b.
carocamargo: carocamargo/regionalSLB: code for regional sea-level budget manuscript (publish), Zenodo, [code], https://doi.org/10.5281/zenodo.7515755, last access 9 January 2023.
Carton, J. A., Chepurin, G. A., and Chen, L.: SODA3: A new ocean climate
reanalysis, J. Climate, 31, 6967–6983,
, 2018.
Cazenave, A. and Moreira, L.: Contemporary sea-level changes from global to
local scales: a review, Philos. T. Roy. Soc. A, 478, 1–30,
2022.
Chafik, L., Nilsen, J. E. Ø., Dangendorf, S., Reverdin, G., and Frederikse,
T.: North Atlantic Ocean Circulation and Decadal Sea Level Change During the
Altimetry Era, Sci. Rep., 9, 1–9,
, 2019.
Chambers, D. P., Wahr, J., and Nerem, R. S.: Preliminary observations of
global ocean mass variations with GRACE, Geophys. Res. Lett., 31,
1–4,
, 2004.
Chambers, D. P., Cazenave, A., Champollion, N., Dieng, H., Llovel, W.,
Forsberg, R., von Schuckmann, K., and Wada, Y.: Evaluation of the Global
Mean Sea Level Budget between 1993 and 2014, Surv. Geophys., 38,
309–327,
, 2017.
Chen, J., Tapley, B., Wilson, C., Cazenave, A., Seo, K. W., and Kim, J. S.:
Global Ocean Mass Change From GRACE and GRACE Follow-On and Altimeter and
Argo Measurements, Geophys. Res. Lett., 47, 1–9,
, 2020.
Chen, J. L., Wilson, C. R., Tapley, B. D., and Grand, S.: GRACE detects
coseismic and postseismic deformation from the Sumatra-Andaman earthquake,
Geophys. Res. Lett., 34, 1–5,
, 2007.
Cheng, L., Abraham, J., Hausfather, Z., and Trenberth, K. E.: How fast are the
oceans warming?, Sci. Clim. Change, 128, 1–8,
, 2019.
CMEMS: Global ocean gridded l4 sea surface heights and derived variables
reprocessed (1993-Ongoing) SEALEVEL_GLO_PHY_L4_MY_008_047,
, 2022.
CSIRO: Combined TOPEX Poseidon, Jason-1, Jason-2 OSTM, Jason-3 near-global
gridded monthly-average sea level product.,
(last access: 9 January 2023), 2022.
Dangendorf, S., Calafat, F. M., Arns, A., Wahl, T., Haigh, I. D., and Jensen,
J.: Mean sea level variability in the North Sea: processes and
implications, J. Geophys. Res.-Oceans, 119, 6820–6841,
, 2014.
Dangendorf, S., Frederikse, T., Chafik, L., Klinck, J. M., Ezer, T., and
Hamlington, B. D.: Data-driven reconstruction reveals large-scale ocean
circulation control on coastal sea level, Nat. Clim. Change, 11,
514–520,
, 2021.
Falasca, F., Bracco, A., Nenes, A., and Fountalis, I.: Dimensionality
Reduction and Network Inference for Climate Data Using
-MAPS:
Application to the CESM Large Ensemble Sea Surface Temperature, J. Adv. Model. Earth Sy., 11, 1479–1515,
, 2019.
Falasca, F., Crétat, J., Braconnot, P., and Bracco, A.: Spatiotemporal
complexity and time-dependent networks in sea surface temperature from mid-
to late Holocene, Europ. Phys. J. Plus, 135, 1–21,
, 2020.
Farrell, W. E. and Clark, J. A.: On Postglacial Sea Level, Geophys. J. Roy. Astro. Soc., 46, 647–667, 1976.
Feng, M., Li, Y., and Meyers, G.: Multidecadal variations of Fremantle sea
level: Footprint of climate variability in the tropical Pacific, Geophys. Res. Lett., 31, 2–5,
, 2004.
Fountalis, I., Dovrolis, C., Bracco, A., Dilkina, B., and Keilholz, S.:
-MAPS: from spatio-temporal data to a weighted and lagged network
between functional domains, Appl. Net. Sci., 3, 3–21,
, 2018.
Fox-Kemper, B., Hewitt, H. T., Xiao, C., Aðalgeirsdóttir, G., Drijfhout,
S. S., Edwards, T. L., Hemer, M., Kopp, R. E., Krinner, G., Mix, A., Notz,
D., Nowicki, S., Nurhati, I., Ruiz, L., Sallée, J.-B., Slangen, A.
B. A., and Yu, Y.: Ocean, Cryosphere and Sea Level Change, in: The Physical
Science Basis. Contribution of Working Group I to the Sixth Assessment Report
of the Intergovernmental Panel on Climate Change, edited by Masson-Delmotte,
V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S., Caud,
N., Chen, Y., Goldfard, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy,
E., Matthews, J. B. R., Maaycock, T. K., Waterfield, T., Yeleçi, O.,
Yu, R., and Zhou, B., September, chap. 9, Cambridge
University Press, Cambridge, 1211–1362,
, 2021.
Frederikse, T., Riva, R., Kleinherenbrink, M., Wada, Y., van den Broeke, M.,
and Marzeion, B.: Closing the sea level budget on a regional scale: Trends
and variability on the Northwestern European continental shelf, Geophys.
Res. Lett., 43, 10810–864872,
, 2016.
Frederikse, T., Simon, K., Katsman, C. A., and Riva, R.: The sea-level budget
along the Northwest Atlantic coast: GIA, mass changes, and large-scale ocean
dynamics., J. Geophys. Res.-Oceans, 122, 5486–5501,
, 2017.
Frederikse, T., Jevrejeva, S., Riva, R. E., and Dangendorf, S.: A consistent
sea-level reconstruction and its budget on basin and global scales over
1958-2014, J. Climate, 31, 1267–1280,
, 2018.
Frederikse, T., Landerer, F., Caron, L., Adhikari, S., Parkes, D., Humphrey,
V. W., Dangendorf, S., Hogarth, P., Zanna, L., and Cheng, L.: The causes of
sea-level rise since 1900, Nature, 584, 393–397,
, 2020.
Gaillard, F., Reynaud, T., Thierry, V., Kolodziejczyk, N., and Von
Schuckmann, K.: In situ-based reanalysis of the global ocean temperature
and salinity with ISAS: Variability of the heat content and steric height,
J. Climate, 29, 1305–1323,
, 2016.
Garric, G. and Parent, L.: Quality Information Document For products
GLOBAL-REANALYSIS-PHY-001-025, Tech. Rep., Copernicus Marine Environment
Monitoring Service,
, 2017.
Ghobadi-Far, K., Han, S. C., Allgeyer, S., Tregoning, P., Sauber, J.,
Behzadpour, S., Mayer-Gürr, T., Sneeuw, N., and Okal, E.: GRACE
gravitational measurements of tsunamis after the 2004, 2010, and 2011 great
earthquakes, J. Geodesy, 94, 1–9,
2020.
Gill, A. E. and Niller, P. P.: The theory of the seasonal variability in the
ocean, Deep Sea Res., 20, 141–177,
, 1973.
Good, S. A., Martin, M. J., and Rayner, N. A.: EN4: Quality controlled ocean
temperature and salinity profiles and monthly objective analyses with
uncertainty estimates, J. Geophys. Res.-Oceans, 118,
6704–6716,
, 2013.
Gregory, J. M., White, N. J., Church, J. A., Bierkens, M. F., Box, J. E., Van
Den Broeke, M. R., Cogley, J. G., Fettweis, X., Hanna, E., Huybrechts, P.,
Konikow, L. F., Leclercq, P. W., Marzeion, B., Oerlemans, J., Tamisiea,
M. E., Wada, Y., Wake, L. M., and Van De Wal, R. S.: Twentieth-century
global-mean sea level rise: Is the whole greater than the sum of the parts?,
J. Climate, 26, 4476–4499,
, 2013.
Gregory, J. M., Griffies, S. M., Hughes, C. W., Lowe, J. A., Church, J. A.,
Fukimori, I., Gomez, N., Kopp, R. E., Landerer, F., Cozannet, G. L., Ponte,
R. M., Stammer, D., Tamisiea, M. E., and van de Wal, R. S.: Concepts and
Terminology for Sea Level: Mean, Variability and Change, Both Local and
Global, Surv. Geophys., 40, 1251–1289,
, 2019.
Griffies, S. M. and Greatbatch, R. J.: Physical processes that impact the
evolution of global mean sea level in ocean climate models, Ocean Modell.,
51, 37–72,
, 2012.
Guinehut, S., Dhomps, A.-L., Larnicol, G., and Le Traon, P.-Y.: High resolution 3-D temperature and salinity fields derived from in situ and satellite observations, Ocean Sci., 8, 845–857,
, 2012.
Hamlington, B. D., Gardner, A. S., Ivins, E., Lenaerts, J. T., Reager, J. T.,
Trossman, D. S., Zaron, E. D., Adhikari, S., Arendt, A., Aschwanden, A.,
Beckley, B. D., Bekaert, D. P., Blewitt, G., Caron, L., Chambers, D. P.,
Chandanpurkar, H. A., Christianson, K., Csatho, B., Cullather, R. I.,
DeConto, R. M., Fasullo, J. T., Frederikse, T., Freymueller, J. T., Gilford,
D. M., Girotto, M., Hammond, W. C., Hock, R., Holschuh, N., Kopp, R. E.,
Landerer, F., Larour, E., Menemenlis, D., Merrifield, M., Mitrovica, J. X.,
Nerem, R. S., Nias, I. J., Nieves, V., Nowicki, S., Pangaluru, K., Piecuch,
C. G., Ray, R. D., Rounce, D. R., Schlegel, N. J., Seroussi, H., Shirzaei,
M., Sweet, W. V., Velicogna, I., Vinogradova, N., Wahl, T., Wiese, D. N., and
Willis, M. J.: Understanding of Contemporary Regional Sea-Level Change and
the Implications for the Future, Rev. Geophys., 58, 1–39,
, 2020.
Han, W., Meehl, G. A., Stammer, D., Hu, A., Hamlington, B., Kenigson, J.,
Palanisamy, H., and Thompson, P.: Spatial Patterns of Sea Level Variability
Associated with Natural Internal Climate Modes, Surv. Geophys., 38,
217–250,
, 2017.
Hardman-Mountford, N. J., Richardson, A. J., Boyer, D. C., Kreiner, A., and
Boyer, H. J.: Relating sardine recruitment in the Northern Benguela to
satellite-derived sea surface height using a neural network pattern
recognition approach, Prog. Oceanogr., 59, 241–255,
, 2003.
Hermans, T. H., Le Bars, D., Katsman, C. A., Camargo, C. M., Gerkema, T.,
Calafat, F. M., Tinker, J., and Slangen, A. B.: Drivers of Interannual Sea
Level Variability on the Northwestern European Shelf, J. Geophys. Res.-Oceans, 125, 1–18,
, 2020.
Hernández-Carrasco, I. and Orfila, A.: The Role of an Intense Front on
the Connectivity of the Western Mediterranean Sea: The Cartagena-Tenes
Front, J. Geophys. Res.-Oceans, 123, 4398–4422,
, 2018.
Hughes, C. W. and Williams, S. D.: The color of sea level: Importance of
spatial variations in spectral shape for assessing the significance of
trends, J. Geophys. Res.-Oceans, 115, 1–18,
, 2010.
Hughes, C. W., Fukumori, I., Griffies, S. M., Huthnance, J. M., Minobe, S.,
Spence, P., Thompson, K. R., and Wise, A.: Sea Level and the Role of Coastal
Trapped Waves in Mediating the Influence of the Open Ocean on the Coast,
Surv. Geophys., 40, 1467–1492,
2019.
Ishii, M. and Kimoto, M.: Reevaluation of historical ocean heat content
variations with time-varying XBT and MBT depth bias corrections, J.
Oceanogr., 65, 287–299,
, 2009.
Ishii, M., Yasui, S., Fukuda, Y., Sato, K., Hirahara, S., and Suzuki, T.:
Accuracy of Global Upper Ocean Heat Content Estimation Expected from Present
Observational Data Sets, Sola, 13, 163–167,
2017.
Iskandar, I.: Variability of Satellite-Observed Sea Surface Height in the
Tropical Indian Ocean: Comparison of Eof and Som Analysis, Makara J. Sci., 13, 173–179,
, 2009.
Kleinherenbrink, M., Riva, R., Frederikse, T., Merrifield, M., and Wada, Y.:
Trends and interannual variability of mass and steric sea level in the
Tropical Asian Seas, J. Geophys. Res.-Oceans, 122,
6254–6276,
, 2017.
Kohonen, T.: Self-organized formation of topologically correct feature maps,
Biol. Cybernet., 43, 59–69,
, 1982.
Landerer, F. W., Jungclaus, J. H., and Marotzke, J.: Ocean bottom pressure
changes lead to a decreasing length-of-day in a warming climate, Geophys.
Res. Lett., 34, 1–5,
2007.
Li, H., Xu, F., Zhou, W., Wang, D., Wright, J. S., Liu, Z., and Lin, Y.:
Development of a global gridded Argo data set with Barnes successive
corrections, J. Geophys. Res.-Oceans, 122, 866–889,
, 2017.
Liu, Y. and Weisberg, R. H.: Patterns of ocean current variability on the West
Florida Shelf using the self-organizing map, J. Geophys. Res.-Oceans, 110, 1–12,
, 2005.
Liu, Y. and Weisberg, R. H.: Ocean currents and sea surface heights estimated
across the west Florida shelf, J. Phys. Ocean., 37,
1697–1713,
, 2007.
Liu, Y., Weisberg, R. H., and Mooers, C. N.: Performance evaluation of the
self-organizing map for feature extraction, J. Geophys. Res.-Oceans, 111, 1–14,
, 2006.
Liu, Y., Weisberg, R. H., and Yuan, Y.: Patterns of upper layer circulation
variability in the South China Sea from satellite altimetry using the
self-organizing map, Acta Ocean. Sin., 27, 129–144, 2008.
Liu, Y., Weisberg, R. H., Vignudelli, S., and Mitchum, G. T.: Patterns of the
loop current system and regions of sea surface height variability in the
eastern Gulf of Mexico revealed by the self-organizing maps, J. Geophys. Res.-Oceans, 121, 2347–2366,
2016.
Lu, S., Liu, Z., Li, H., Li, Z., Wu, X., and Sun, C.: Manual of Global Ocean
Argo gridded data set (BOA_Argo) (Version 2019), Tech. Rep., Second
Institue of Oceanography SOA,
, 2019.
Ma, F., Diao, Y. N., and Luo, D. H.: Spatiotemporal characteristics of the sea
level anomaly in the Kuroshio Extension using a self-organizing map,
Atmos. Ocean. Sci. Lett., 9, 471–478,
, 2016.
MacIntosh, C. R., Merchant, C. J., and von Schuckmann, K.: Uncertainties in
Steric Sea Level Change Estimation During the Satellite Altimeter Era:
Concepts and Practices, Surv. Geophys., 38, 59–87,
, 2017.
Maclachlan, C., Arribas, A., Peterson, K. A., Maidens, A., Fereday, D., Scaife,
A. A., Gordon, M., Vellinga, M., Williams, A., Comer, R. E., Camp, J.,
Xavier, P., Madec, G., and National, F.: Global Seasonal forecast system
version 5 (GloSea5): a high-resolution seasonal forecast system, Q. J. Roy. Meteor. Soc., 5, 1072–1084,
, 2015.
Marcos, M. and Amores, A.: Quantifying anthropogenic and natural contributions
to thermosteric sea level rise, Geophys. Res. Lett., 41,
2502–2507,
, 2014.
Nickerson, A. K., Weisberg, R. H., and Liu, Y.: On the Evolution of the Gulf
of Mexico Loop Current Through Its Penetrative, Ring Shedding and Retracted
States, Adv. Space Res., 69, 4058–4077,
, 2022.
Novi, L., Bracco, A., and Falasca, F.: Uncovering marine connectivity through
sea surface temperature, Sci. Rep., 11, 1–9,
, 2021.
Pattiaratchi, C. and Siji, P.: Variability in ocean currents around
Australia, in: State and Trends of Australia's Ocean Report, edited by:
Richardson, A., Eriksen, R., Moltmann, T., Hodgson-Johnston, I., and Wallis,
J., January 2020, Chap. 1.4 Variab, Integrated Marine Observing
System (IMOS), 1–6,
, 2020.
Peltier, W. R., Argus, D. F., and Drummond, R.: Space geodesy constrains ice
age terminal deglaciation: The global ICE-6G_C (VM5a) model, J. Geophys. Res.-Sol. Ea., 120, 450–487,
, 2015.
Purkey, S. G. and Johnson, G. C.: Warming of Global Abyssal and Deep Southern
Ocean Waters between the 1990s and 2000s: Contributions to Global Heat and
Sea Level Rise Budgets, J. Climate, 23, 6336–6351,
, 2010.
Purkey, S. G., Johnson, G. C., and Chambers, D. P.: Sea contributions of ocean
mass and deep steric changes to sea level rise between 1993 and 2013,
J. Geophys. Res.-Oceans, 119, 7509–7522,
, 2014.
Purkey, S. G., Johnson, G. C., Talley, L. D., Sloyan, B. M., Wijffels, S. E.,
Smethie, W., Mecking, S., and Katsumata, K.: Unabated Bottom Water Warming
and Freshening in the South Pacific Ocean, J. Geophys. Res.-Oceans, 124, 1778–1794,
, 2019.
Richardson, A. J., Risi En, C., and Shillington, F. A.: Using
self-organizing maps to identify patterns in satellite imagery, Prog.
Oceanogr., 59, 223–239,
, 2003.
Rietbroek, R., Brunnabend, S.-E., Kusche, J., Schröter, J., and Dahle,
C.: Revisiting the contemporary sea-level budget on global and regional
scales, P. Natl. Acad. Sci. USA, 113, 1504–1509,
, 2016.
Roemmich, D. and 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, 81–100,
, 2009.
Rougier, J.: Ensemble averaging and mean squared error, J. Climate,
29, 8865–8870,
, 2016.
Royston, S., Watson, C. S., Legresy, B., King, M. A., Church, J. A., and Bos,
M. S.: Sea-Level Trend Uncertainty With Pacific Climatic Variability and
Temporally-Correlated Noise, J. Geophys. Res.-Oceans, 123,
1–16,
, 2018.
Royston, S., Vishwakarma, B. D., Westaway, R., Rougier, J., Sha, Z., and
Bamber, J.: Can We Resolve the Basin-Scale Sea Level Trend Budget From GRACE
Ocean Mass?, J. Geophys. Res.-Oceans, 125, 1–16,
, 2020.
Save, H.: CSR GRACE and GRACE-FO RL06 Mascon Solutions v02, CSR, 2,
, 2020.
Save, H., Bettadpur, S., and Tapley, B. D.: High-resolution CSR GRACE RL05
mascons, J. Geophys. Res.-Sol. Ea., 121, 7547–7569,
, 2016.
Schwarz, G.: Estimating the dimension of a model, Ann. Stat., 6,
461–464,
, 1978.
Slangen, A. B. A., van de Wal, R. S. W., Wada, Y., and Vermeersen, L. L. A.: Comparing tide gauge observations to regional patterns of sea-level change (1961–2003), Earth Syst. Dynam., 5, 243–255,
, 2014.
Slangen, A. B., Church, J. A., Agosta, C., Fettweis, X., Marzeion, B., and
Richter, K.: Anthropogenic forcing dominates global mean sea-level rise
since 1970, Nat. Clim. Change, 6, 701–705,
2016.
Slangen, A. B. A., Adloff, F., Jevrejeva, S., Leclercq, P. W., Marzeion, B.,
Wada, Y., and Winkelmann, R.: A Review of Recent Updates of Sea-Level
Projections at Global and Regional Scales, Surv. Geophys., 38,
385–406,
, 2017.
SLcci: Time series of gridded sea level anomalies,
2022.
Sonnewald, M., Wunsch, C., and Heimbach, P.: Linear predictability: A sea
surface height case study, J. Climate, 31, 2599–2611,
, 2018.
Sonnewald, M., Wunsch, C., and Heimbach, P.: Unsupervised Learning Reveals
Geography of Global Ocean Dynamical Regions, Earth Space Sci., 6,
784–794,
, 2019.
Stammer, D., Cazenave, A., Ponte, R. M., and Tamisiea, M. E.: Causes for
Contemporary Regional Sea Level Changes, Annu. Rev. Mar. Sci., 5, 21–46,
, 2013.
Storto, A. and Masina, S.: C-GLORSv5: an improved multipurpose global ocean eddy-permitting physical reanalysis, Earth Syst. Sci. Data, 8, 679–696,
, 2016.
Storto, A., Masina, S., Balmaseda, M., Guinehut, S., Xue, Y., Szekely, T.,
Fukumori, I., Forget, G., Chang, Y. S., Good, S. A., Köhl, A.,
Vernieres, G., Ferry, N., Peterson, K. A., Behringer, D., Ishii, M., Masuda,
S., Fujii, Y., Toyoda, T., Yin, Y., Valdivieso, M., Barnier, B., Boyer, T.,
Lee, T., Gourrion, J., Wang, O., Heimback, P., Rosati, A., Kovach, R.,
Hernandez, F., Martin, M. J., Kamachi, M., Kuragano, T., Mogensen, K., Alves,
O., Haines, K., and Wang, X.: Steric sea level variability (1993–2010) in
an ensemble of ocean reanalyses and objective analyses, Clim. Dynam.,
49, 709–729,
, 2017.
Szekely, T., Gourrion, J., Pouliquen, S., and Reverdin, G.: The CORA 5.2 dataset for global in situ temperature and salinity measurements: data description and validation, Ocean Sci., 15, 1601–1614,
, 2019.
Tapley, B. D., Bettadpur, S., Watkins, M., and Reigber, C.: The Gravity
Recovery and Climate Experiment: Mission Overview and Early Results,
Geophys. Res. Lett., 31, 1–4,
, 2004.
Thompson, P. R. and Merrifield, M. A.: A unique asymmetry in the pattern of
recent sea level change, Geophys. Res. Lett., 41, 7675–7683,
, 2014.
Vesanto, J., Himberg, J., Alhoniemi, E., and Parhankangaas, J.: SOM Toolbox
for Matlab 5, Tech. Rep., Helsinki University of Technology, ISBN 951-22-4951-0, 2000.
Wagner, P. and Böning, C. W.: Decadal sea-level variability in the Australasian Mediterranean Sea, Ocean Sci., 17, 1473–1487,
, 2021.
Wang, J., Church, J. A., Zhang, X., Gregory, J. M., Zanna, L., and Chen, X.:
Evaluation of the Local Sea-Level Budget at Tide Gauges Since 1958,
Geophys. Res. Lett., 48, 1–12,
, 2021.
Watkins, M. M., Wiese, D. N., Yuan, D.-N., Boening, C., and Landerer, F. W.:
Improved methods for observing Earth's time variable mass distribution with
GRACE using spherical cap mascons, J. Geophys. Res.-Sol. Ea., 120, 1648–2671,
, 2015.
WCRP Global Sea Level Budget Group: Global sea-level budget 1993–present, Earth Syst. Sci. Data, 10, 1551–1590,
, 2018.
Weisberg, R. H. and Liu, Y.: On the Loop Current Penetration into the Gulf of
Mexico, J. Geophys. Res.-Oceans, 122, 9679–9694,
, 2017.
Wiese, D. N., Yuan, D.-N., Boening, C., Landerer, F. W., and Watkins, M.:. JPL GRACE Mascon Ocean, Ice, and Hydrology Equivalent Water Height RL06 CRI Filtered Version 02. Ver. 02. PO.DAAC, CA, USA, [data set], https://doi.org/10.5067/TEMSC-3JC62, 2019.
Wöppelmann, G. and Marcos, M.: Vertical land motion as a key to
understanding sea level change and variability, Rev. Geophys., 54,
64–92,
, 2015.
Zlotnicki, V., Qu, Z., Willis, J. K., Ray, R., and Hausman, J.: JPL MEASURES
Gridded Sea Surface Height Anomalies Version JPL1812, Tech. Rep., PO.DAAC,
CA, USA,
, 2019.
Zuo, H., Balmaseda, M. A., Tietsche, S., Mogensen, K., and Mayer, M.: The ECMWF operational ensemble reanalysis–analysis system for ocean and sea ice: a description of the system and assessment, Ocean Sci., 15, 779–808,
, 2019.
Articles
Abstract
The sea-level budget
Data and methods
Identifying domains of sea-level variability
The regional sea-level budget on different spatial scales
Discussion and conclusions
Appendix A:
Dynamic sea-level change estimation and validation
Code and data availability
Author contributions
Competing interests
Disclaimer
Financial support
Review statement
References
Article
(13021 KB)
Full-text XML
BibTeX
EndNote
Editorial statement
Closing the sea level budget (between observed rise and the sum of its causes) has been a challenge, is an ongoing effort and has primarily concerned the global mean. Here, the authors use machine learning to identify sub-areas with similar trends to close the sea level budget on a regional level, with much reduced errors compared with 1-degree grid points.
Closing the sea level budget (between observed rise and the sum of its causes) has been a...
Short summary
Sea-level change is mainly caused by variations in the ocean’s temperature and salinity and land ice melting. Here, we quantify the contribution of the different drivers to the regional sea-level change. We apply machine learning techniques to identify regions that have similar sea-level variability. These regions reduce the observational uncertainty that has limited the regional sea-level budget so far and highlight how large-scale ocean circulation controls regional sea-level change.
Sea-level change is mainly caused by variations in the ocean’s temperature and salinity and...
Share
Sections
Abstract
The sea-level budget
Data and methods
Identifying domains of sea-level variability
The regional sea-level budget on different spatial scales
Discussion and conclusions
Appendix A:
Dynamic sea-level change estimation and validation
Code and data availability
Author contributions
Competing interests
Disclaimer
Financial support
Review statement
References