the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Sea level reconstruction reveals improved separations of regional climate and trend patterns over the last seven decades
Abstract. Rapidly rising sea level is one of the major adverse consequences of anthropogenic climate change. Sea level rise poses an existential threat to coastal populations, particularly for urban settlements with accelerating growth rates. Contemporary empirical sea level reconstructions have been used to conflate short-term (~3 decades) satellite altimetry geocentric sea level data and long-term (50 years or longer) tide gauge records to better estimate reliable sea level rise towards multi-decadal to centennial time scales. However, adequate separations and quantifications of low-frequency climate patterns and sea level trends globally at regional scales remain elusive. Here, we propose a new sea level reconstruction framework that incorporates Empirical Orthogonal Function (EOF) into the contemporary Cyclostationary EOF with Reduced Space Optimal Interpolation (CSEOF-OI) algorithm to better reconstruct sea level fields. Using 225 selected long-term gap-filled tide gauge records with vertical land motion adjusted and satellite altimetry, our global reconstructed monthly sea level time series, January 1950– January 2022, exhibit distinct delineations between modeled climate patterns and sea level trends at regional scales. The separated sea level patterns include trends, modulated annual cycles, the El Niño Southern Oscillation (ENSO), and the Pacific Decadal Oscillation (PDO). The third principal component of the reconstructed sea level exhibits a Pearson correlation coefficient of 0.87 with the Niño 3.4 ENSO index, and the fourth principal component correlates at 0.75 with the PDO index, indicating excellent agreement. The global mean sea level trend, accounting for the predominant climate periodicities, is 1.9 ± 0.2 mm yr⁻¹ (95 % confidence), and the estimate during the satellite altimetry era (January 1993–December 2021) is 3.2 ± 0.3 mm yr⁻¹ (95 % confidence). Compared with previous studies, we conclude that our 72-year sea-level reconstruction allows us to better separate the ENSO and PDO climate patterns, as well as the sea level they induced. Finally, we show that the short-term (5-year) rates of ENSO and PDO patterns significantly affect sea level both on a global and regional scale, altering global mean sea level trends by up to 1.1 ± 0.5 mm yr¹ (January 2011–January 2016). Over the past seven decades, the climate pattens exerted a minor impact on sea level trends, but substantially modulated apparent regional sea level accelerations, particularly in the western Pacific (e.g., 0.09 ± 0.05 mm yr⁻² at the Kuroshio Current), and in the east and central equatorial Pacific Ocean (e.g., −0.04 ± 0.03 mm yr⁻² near Costa Rica). The reconstructed sea level and analysis results datasets are available at https://doi.org/10.5281/zenodo.15288817 (Wang, 2025).
- Preprint
(33773 KB) - Metadata XML
-
Supplement
(31767 KB) - BibTeX
- EndNote
Status: final response (author comments only)
- RC1: 'Comment on essd-2025-251', Anonymous Referee #1, 24 Jun 2025
-
RC2: 'Comment on essd-2025-251', Anonymous Referee #2, 07 Oct 2025
This manuscript introduces a modified spatial-temporal sea-level reconstruction method/product that combines two complementary decomposition techniques with reduced-space optimal interpolation, extending geocentric (absolute) sea-level anomalies back to 1950. The manuscript provides two main improvements. First, the reconstruction products support more cleanly separating climate-related sea level variability (e.g., ENSO, PDO) from the long-term trend within a spatiotemporal framework. Second, after quantifying climate-related variabilities, the authors assess their influence on regional short-term trend estimates as well as on regional long-term trends and accelerations. Other minor contributions include testing different gap-filling methods for tide gauge records to increase record availability. Comparison between different ways to correct vertical land motion at tide gauges to better match tide-gauge and satellite altimetry sea level. Overall, the datasets and their derived results are likely to benefit the sea-level and climate communities. The authors have provided both the products and core plotting/analysis codes at the data availability link, which facilitates a user-friendly experience. However, I still have several comments/concerns that, once addressed, should make the manuscript suitable for publication.
- The manuscript tests three gap-filling approaches (regularized EM, PPCA, AR) for tide gauges by using 48 nearly complete 70-year records and concludes that regularized EM performs best, as illustrated in Figure 2. Since the reconstruction process uses many more gap-filled tide gauges, please include a robustness test on a broader sample, e.g., considering setting the gauge record length shorter to include a larger set of near-complete gauge records, and report whether the ranking of methods is stable.
- In the section GNSS-derived vertical land motion at tide-gauge locations compared with the results from the altimetry minus tide-gauge (SA-TG) approach, this comparison should provide the GNSS data selection criterion used in Figure 3, e.g., collocation criteria, time spans, etc. More generally, in this section, please refocus SA-TG approach, which was used in the reconstruction, and present the other two correction approaches briefly. This will avoid letting side methods overshadow the main workflow and will make the procedure easier to follow.
- For the validation on the modified sea level reconstruction product in Figure 6d, the manuscript emphasizes applying long-term tide gauge records excluded from the reconstruction as independent validations, which is good. However, it should also report the comparison between the gauges used in the reconstruction process and their collocated reconstruction series to assess whether potential overfitting exists. For Figure 6g, since produced from this study in 1° × 1° resolution, please also compare the full-resolution 70-year trends with other gridded reconstruction over their common period if comparable-resolution product exists.
- In Figure 9, the 70-year combined ENSO–PDO induced sea-level acceleration is spatially prominent. To better serve users who are interested in multi-decadal sea level acceleration studies, the authors should also include separate spatial maps showing the ENSO-induced and PDO-induced acceleration coefficients, respectively.
- Add a flowchart of the modified reconstruction procedure to highlight the steps implemented in the reconstruction section more clearly.
Minor suggestions:
Line 25, exhibit -> exhibits
Line 68, spatiotemporal sea level changes -> spatiotemporal changes in sea level
Line 103, the estimate of -> estimating, difficult -> challenging
Line 112, mirrors -> was similar to
Line 134, PSMSL, first time showing in the main text, should provide the full name.
Line 162, mitigate -> address
Line 189, needed-> necessary
Line 218, long-> largely
Line 225, to be estimated -> estimation
Line 227 to Line 229, “An issue that might be regarded as a mere matter of language is examined here … entails circular logic or tautological reasoning.” Should simplify/reorganize these two sentences.
Line 234, not -> rather than
Line 324, corresponding the -> the corresponding
Line 335, In Figure 4, the author should add the trend values of the CSEOF trend mode estimated regional trends, and trend values from the further refined trend mode. Also, anomlies -> anomalies
Line 393, initial and terminal -> near beginning and end
Line 450, and respective reconstruction -> and their respective reconstructions
Line 464 to Line 466 “This validation was performed over 40 years (January 1982 to December 2021; 224 gauges) and 65 years (January 1967 to December 2021; 34 gauges), yielding a consistent median r0.85 (85.19% for 40 years and 84.76% for 65 years).” Since the reconstructed sea level is from January 1950 to January 2022, consider updating this comparison to also include January 2022.
Line 541, in line -> consistent
Line 560, before -> prior to the
Line 577, suggesting -> suggest, meet -> meeting
Line 578, partly -> have been partly
Line 586, ENSO, PDO‐induced sea level -> ENSO and PDO-induced sea level
Line 621, increased -> increasing
Line 627, evident by most of the strong rising areas concentrating on the low-latitude Pacific -> as evident by the concentration of the strongest rising areas in the low-latitude Pacific.
Line 641, contrasted by -> in contrast to
Line 642, western -> West
Line 652, the critical need to account for both ENSO and PDO influences for understanding short-term sea level changes. -> the importance of considering both ENSO and PDO influences in understanding short-term sea level changes.
Line 658, during -> over
Line 705, generally between 1 and 2 mm yr⁻¹ -> typically ranging from 1 to 2 mm yr⁻¹
Line 730, estimate -> estimating
Line 750, scales -> timescales
Line 762, multi-decades -> multiple decades
Citation: https://doi.org/10.5194/essd-2025-251-RC2
Data sets
Modified Sea Level Reconstruction Reveals Improved Separation of Climate and Trend Patterns Shengdao Wang https://doi.org/10.5281/zenodo.15288817
Model code and software
Modified Sea Level Reconstruction Reveals Improved Separation of Climate and Trend Patterns Shengdao Wang https://doi.org/10.5281/zenodo.15288817
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
1,353 | 72 | 23 | 1,448 | 38 | 26 | 45 |
- HTML: 1,353
- PDF: 72
- XML: 23
- Total: 1,448
- Supplement: 38
- BibTeX: 26
- EndNote: 45
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
The manuscript presents a global sea level reconstruction spanning 1950–2022 at a 1° × 1° resolution, using a combined dataset of tide gauge records (corrected for vertical land motion) and satellite altimetry, implemented through an enhanced CSEOF-OI framework incorporating EOF decomposition. The authors emphasize improved separation of climate modes (e.g., ENSO, PDO) and long-term trends.
While the study contributes to an important topic and provides a potentially valuable dataset to the sea level and climate research community, several key aspects of the methodology and validation require further clarification. In particular, the validation appears limited to coarser 5° × 5° spatial scales, and insufficient attention is given to the accuracy of reconstructed data in the pre-altimetry era (especially before 1993) and in open-ocean regions lacking tide gauge constraints. These issues, along with others detailed below, should be addressed to ensure the robustness and usability of the dataset.