Pre-registration: gauge trend, land motion, splice to AR6 (written 13 Nov 2026 before any trend number was computed)
I have looked at the raw annual values (they sit between 6.80 and 7.12 m on the PSMSL datum) but have fitted nothing.
Quantity
Relative sea level at Constanța harbour: the sea surface measured against the gauge's own benchmark. It includes any sinking or rising of the quay. AR6 total is also relative sea level, so no land-motion correction is applied to either; the two are compared as they are.
Primary trend
OLS of annual RLR values on year, over the 59 years present, 1933–1996. Reported in mm/year with a 95% interval from the OLS standard error inflated for AR(1) residuals (effective sample size n_eff = n (1−r1)/(1+r1), r1 = lag-1 residual autocorrelation measured on the present consecutive pairs, t distribution with n_eff−2 degrees of freedom). The rise is called distinguishable from zero only if that interval excludes 0.
Sensitivities (all reported, none can replace the primary)
- Drop flagged years 1985 and 1992.
- First half (1933–1964) vs second half (1965–1996) slopes (is it one straight line?).
- Leave-ten-years-out jackknife (all 54 contiguous 10-year blocks removed in turn): range of the slope.
- Quadratic term: reported, with its AR(1)-corrected interval; no acceleration is claimed unless it excludes 0.
Pre-registered test (can fail, and is shown either way)
Prediction: the Constanța relative rise 1933–1996, after removing the AR6 background land motion at location 379, matches the global-mean rise. Global comparison, quoted from AR6 WG1 Chapter 2 (checked in the search result of 13 Nov 2026): 1.3 [0.6–2.1] mm/year for 1901–1971 and 1.9 [0.8–2.9] mm/year for 1971–2006; the comparison interval for 1933–1996 is therefore 0.6–2.9 mm/year (central 1.3–1.9). Land motion is the AR6 verticallandmotion median at location 379 read as the 2020–2100 slope from the ssp245 file (sign convention: see Correction 1 below). Pass = (gauge trend minus sinking rate) lies within 0.6–2.9 mm/year. Fail = outside; shown on the page either way. AR6's background land motion at tide-gauge sites is partly estimated from tide-gauge records, so the check is not independent of this record; the page says so.
Splice to AR6 (gauge stops 1996, AR6 baseline is the 1995–2014 mean)
The page's zero is "the 1995–2014 mean sea level", as in AR6. The gauge cannot measure it (no data after 1996). Rule: zero is the primary OLS line evaluated at 2004.5, the middle of the baseline period. The uncertainty of that anchor (95% prediction interval of the line at 2004.5, AR(1)-corrected) is drawn as a band and added in quadrature to every "cm higher than in your birth year" statement as the stated margin. Years 1997–2019 are drawn as a gap (dashed, no data), never as a measured line. The gauge year-to-year scatter (river flow, wind) stays visible as dots.
Reader's number
- Born 1933–1996: N = level(2026) − line(birth year), line from the primary fit; level(2026) from AR6 SSP2-4.5 median linearly between 2020 and 2030 (scenarios differ by under 1 cm by 2030).
- Born 1997–2025: N = level(2026) − level(birth), level before 2020 interpolated linearly from 0 at 2004.5 to the AR6 2020 median. Marked "modelled, no gauge".
- Born before 1933: the gauge starts in 1933; N is given since 1933 and says so.
- By age 80: AR6 median and likely range (17th–83rd percentile) for SSP1-2.6, SSP2-4.5 and SSP5-8.5, plus the margin above. A range, never one number.
- Total vs total-minus-land-motion: shown side by side as a sensitivity.
Black Sea caveat
AR6 projections are built from global models; the Black Sea is enclosed (Bosphorus) and its regional dynamics (Danube inflow, basin-scale response) are poorly represented. The caveat sits beside the fan on the page.
Terrain
Copernicus GLO-30 is a surface model with a stated vertical accuracy of a few metres; heights are relative to the EGM2008 geoid, not to the Black Sea mean level. The shore view shows land lower than a stated height above today's sea, not a flood map: storm surge, lagoon and river levels, sea walls, and DEM error are not modelled, and the page says so.
Corrections and deviations after the first results (13 Nov 2026)
- Sign of land motion (error in this file, found on review). I wrote "negative = sinking". AR6's component files give each process as a contribution to *relative sea level change*, so a negative value means relative sea level is pushed down, i.e. the land is rising. Check: at location 379, ssp245, 2100, the total with land motion has median 433 mm and the total without it (
ar6-regional_novlm-confidence.zip) has median 452 mm: the difference, about −19 mm, matches the land-motion component (−17 mm at 2100). The land-motion slope at Constanța is −0.21 mm/year, so the land is rising by about 0.2 mm/year and the geocentric rise is 1.19 − (−0.21) = 1.40 mm/year. The test statistic and its pass were computed this way from the start; only the words (here, and in the first page draft, "sinking") were wrong. The page says "land rising". - Uncertainty of the zero (deviation). The splice rule above says "95% prediction interval of the line at 2004.5". After seeing it (first run ±17.5 cm) I judged it the wrong quantity: the baseline is a 20-year mean, not a single year. The page uses the 95% interval of the line itself plus the noise of a 20-year mean, ±6.7 cm in total (line ±5.5 cm, 20-year-mean noise ±3.8 cm). All three numbers are in
data/analysis.json (anchor.pi95_mm, ci95_mean_mm, ci95_baseline_mean_noise_mm, ci95_total_mm). - Review fixes to the statistics (13 Nov 2026). (a) The residual standard deviation printed on the page was inflated by the autocorrelation factor (8.2 cm); the page now prints the plain standard deviation of the residuals, 6.8 cm. (b) The lag-1 autocorrelation was computed as a sum over consecutive pairs divided by the sum over all years; it is now the pairwise correlation (0.217 instead of 0.197), which widens the slope margin from ±1.21 to ±1.24 mm/year. The conclusion (the 95% interval just includes zero) is unchanged.
- "Likely" range of the reader's result. The first draft added the full 95% slope margin to AR6's 17th–83rd percentile range. The page now combines the AR6 spread with one standard deviation of the slope margin (margin ÷ 2.02) in quadrature, so the result is a roughly two-thirds range like AR6's own.
- Data after 1996. Found on review: NIMRD (the Romanian marine research institute) describes a continuation of this gauge, monthly 1998–2018 and daily 2007–2018, with national levelling offsets from the 1933 datum, as
NIMRD_Constanta_SeaLevel.zip, 'restricted, download upon request' (EO4SIBS dataset user manual, p. 29, https://sdn2.rmri.ro/esa_app_external_users/user_manual.pdf). No public datum-tied download was found (EMODnet, CMEMS, PSMSL, GLOSS searched). The page uses only the public PSMSL series and says so. - Reader rule for 1997–2004 (clarification). Between 1996 and 2004.5 the page continues the gauge trend line (the zero is defined there); only from 2004.5 does it ramp linearly to the AR6 value for 2020. The reader result for birth years 1997 and later says it is an estimate that combines the two.
- Uncertainty of "N cm higher". The anchor's uncertainty cancels (both ends are measured against the same zero), so the stated margin is the slope's 95% margin times the years before 2004.5, added to the AR6 17th–83rd percentile range. This is a simple linear sum, deliberately wide.
- AR6 quantiles. The dataset's quantile 0.5 is the midpoint of the 0.49 and 0.51 values (384, 433, 481 mm at 2100, ssp245), because AR6 merged two probability ranges there. The page uses the 17th, 50th and 83rd percentiles only, and the 95th for one statement of a high end in the reader table's data, not the 5th.
Back