// NATURE NEWS — SPAZIO & SCIENZA
Addendum: Sea level much higher than assumed in most coastal hazard assessments
Nature
(2026) Cite this article
The Original Article was published on 04 March 2026
Addendum to: Nature https://doi.org/10.1038/s41586-026-10196-1 Published online 4 March 2026
This Addendum provides the configuration details for the geoid models used in the methodology underlying the coastal exposure meta-analysis in the original article1. It expands upon the description in the original publication to ensure full reproducibility of our results. The specific configuration details are outlined below and in Extended Data Tables 1–5.
Geoid model information was obtained from the openly accessible calculation service of the International Centre for Global Earth Models from GFZ Helmholtz Centre for Geosciences (ICGEM)2. Point data on geoid height anomaly to the WGS84 ellipsoid (in m) was obtained for the entire globe at a resolution of 0.085° for the EGM96 (ref. 3), EGM2008 (ref. 4) and GOCO06s (ref. 5) geoid models, respectively6 (data downloaded on 22/08/2023 for EGM96 and EGM2008 and on 10/01/2024 for GOCO06s) (Extended Data Tables 1–3). Specifically, the output grid from ICGEM documents the maximum latitudes of 90.000 °N and 89.945 °S (2118 latitude parallels), and the maximum longitudes of 180 °W and 179.975 °E (4236 longitude parallels). In total, 8,971,848 points were downloaded for each geoid model. We included the zero-degree term as given by ICGEM into the geoid model configuration7 and for each geoid model refer to a tide-free system configuration6. The reference ellipsoid WGS84 underlying these geoid models is defined by the geocentric gravitational constant GM = 3.98600441800 × 1014 m3/s2, a radius of a = 6378137.00 m, a flattening of f = 1/298.25722356300, an angular velocity of ω = 7.292115 × 10–5 s–1 and a normal potential of U0 = 6.263685171456948 × 107 m2/s2. We did not apply any filters in the configuration6. Consequently, the geoid model data we obtained from the ICGEM calculation service include the EGM96 (with a maximum used degree of 360) characterized by a weighted mean of –7.1416531 × 10–3 m, a maximum of 8.5891824 × 101 m, a minimum of –1.0651816 × 102 m and a weighted root mean square of 3.0566520 × 101 m (Extended Data Table 1); EGM2008 (with a maximum used degree of 2190) characterized by a weighted mean of –8.2435599 × 10–3 m, a maximum of 8.6064395 × 101 m, a minimum of –1.0648962 × 102 m and a weighted root mean square of 3.0570532 × 101 m (Extended Data Table 2); and GOCO06s (with a maximum used degree of 300) characterized by a weighted mean of –7.1516965 × 10–3 m, a maximum of 8.5419839 × 101 m, a minimum of –1.0651832 × 102 m and a weighted root mean square of 3.0569639 × 101 m (Extended Data Table 3).
For each geoid model, the height anomaly points were interpolated into a global raster by using multiquadric radial basis functions, which gave the most accurate interpolation results and is also used by gravitational field modelling studies8,9. The interpolations were conducted using the geostatistical wizard in ArcGIS Pro (version 3.1.3). For each interpolation, the kernel parameter was optimized. The standard neighbourhood type was automatically set with a minimum of 10 and a maximum of 15 neighbours. A sector type of one sector was used.
The processing was similar for our additional computations, involving the GOCO2025s (ref. 10) and EGM-DIR R4 (refs. 11,12) investigations, which were added in the later process of our article (Extended Data Tables 4 and 5). For both the GOCO2025s and the DIR R4 geoid models, we downloaded height anomaly data to the WGS84 ellipsoid as specified above6 (date of download for both geoid models: 03/10/2025). We used a grid step of 1.000° and a global coverage, extending from 90.000 °N to 90.000 °S and from 180.000 °W to 180.000 °E, thus obtaining 65,341 grid points. As for the previously described geoid models, we included the zero-degree term as given by ICGEM into the geoid model configuration7 and for each geoid model refer to