We tested whether a larger seasonal mean sea-level contrast also occurs in monthly highest water at five Mid-Atlantic gauges. We used fixed NOAA records and paired September and February in each year. Between 1980-1999 and 2000-2020, the equal-gauge mean contrast rose by 10.36 cm for monthly mean level and 9.24 cm for monthly highest water. The change in highest water minus monthly mean was -1.12 cm, with a 95% three-year block interval from -7.78 to 5.67 cm. The data do not show a clear added seasonal change above the monthly mean. The gap result remains uncertain with other block lengths, common years, and an extension through 2025. The one-year interval for monthly highest change includes zero. The test does not identify storm change, flood frequency, or a physical cause. Missing observations, sampling changes, short periods, and slow ocean variation limit inference. Fixed paired evidence, exact raw-source hashes, code, bootstrap draws, and arithmetic checks support reproduction.
Filters: arxiv: physics.ao-ph (Atmospheric and Oceanic Physics) ×
We tested whether positive autumn Indian Ocean Dipole Mode Index (DMI) improves the probability forecast of a cold next-year ENSO winter. We used fixed HadISST1.1 DMI and NOAA ERSSTv5 Nino 3.4 data. Each logistic fit used earlier complete targets only. The Pacific baseline allowed unequal responses to warm and cold conditions. The test had 44 winters from 1981 to 2024 and 19 cold outcomes at a fixed -0.5 degrees C cutoff. Adding positive DMI increased the Brier score from 0.16335 to 0.16431. The score reduction was -0.59%, with a conditional 95% three-year block interval of -23.25% to +19.59%. All planned check intervals included zero. This gives no clear added skill for the tested model. It does not show equal skill or rule out a physical Indian Ocean effect. A short record, revised data, model choice and conditional uncertainty limit the result. The package contains code, aggregate evidence, source hashes, two solver checks and a figure. External data and record-level tables are regenerated outside the publication.
Question: did the apparent temperature sensitivity (days per °C) of spring phenology at Japan Meteorological Agency (JMA) stations change between 1953–1986 and 1987–2025, do heat-sum models with the same temperatures predict the change, and which model part must change to reproduce it? This is an independent replication of [citation 1], [citation 3] and [citation 2], with a new decomposition. Data: JMA records (regular species) for 8 spring events at 33–82 stations; daily station temperature from GHCN-Daily, with 6,415 gap station-months filled from JMA daily pages and each month anchored to the official JMA monthly mean. Method: per-station OLS slope of event day on the 60-day pre-event mean temperature, for observed days and for days from forcing-only (thermal-time) and two chill-forcing models fitted by grid search with one requirement per station. Uncertainty: year resampling shared by all stations and a permutation test over random year partitions. Results: for cherry first flowering, cherry full bloom, ginkgo budburst and wisteria flowering, the observed sensitivity became 0.73–1.04 days/°C stronger, while the forcing-only model predicted −0.08 to +0.14. Observed minus model: −0.86 (95% year CI −1.50 to −0.27), −0.72 (−1.32 to 0.02), −1.06 (−1.73 to −0.11) and −1.06 (−1.61 to −0.39); permutation p 0.002, 0.017, 0.013, 0.005 (Holm over 8 events: 0.016, 0.085, 0.078, 0.035). Mulberry, azalea, dandelion and ume showed no significant difference. The result held for 30–90-day windows, detrending, control for earlier temperature, official monthly temperature and a neighbour-station instrument. Its size depends on the split year, and moving-window curves are not monotonic. New: a period-specific heat requirement reproduced none of the change (−0.06 to +0.13 days/°C), so a requirement change (for example, from less chilling) cannot explain it in these models. A period-specific base temperature reproduced 36–120% of it; the shift was identified only for wisteria (−5 to +1 °C; year CI of shift 1 to 15 °C). Also: the forcing-only model explains weaker cherry sensitivity at warmer stations (observed gradient +0.41, model +0.36 days/°C per °C), but not the absent gradient for ginkgo, wisteria and mulberry; chill terms cut cherry cross-validated error by 12–17%. Limits: one urban reference plant per station, simple grid models, most events end in 2020, cause not identified.
Question: did the apparent temperature sensitivity (days per °C) of spring phenology at Japan Meteorological Agency (JMA) stations change between 1953–1986 and 1987–2025, do heat-sum models with the same temperatures predict the change, and which model part must change to reproduce it? This is an independent replication of [citation 1], [citation 3] and [citation 2], with a new decomposition. Data: JMA records (regular species) for 8 spring events at 33–82 stations; daily station temperature from GHCN-Daily, with 6,415 gap station-months filled from JMA daily pages and each month anchored to the official JMA monthly mean. Method: per-station OLS slope of event day on the 60-day pre-event mean temperature, for observed days and for days from forcing-only (thermal-time) and two chill-forcing models fitted by grid search with one requirement per station. Uncertainty: year resampling shared by all stations and a permutation test over random year partitions. Results: for cherry first flowering, cherry full bloom, ginkgo budburst and wisteria flowering, the observed sensitivity became 0.73–1.04 days/°C stronger, while the forcing-only model predicted −0.08 to +0.14. Observed minus model: −0.86 (95% year CI −1.50 to −0.27), −0.72 (−1.32 to 0.02), −1.06 (−1.73 to −0.11) and −1.06 (−1.61 to −0.39); permutation p 0.002, 0.017, 0.013, 0.005 (Holm over 8 events: 0.016, 0.085, 0.078, 0.035). Mulberry, azalea, dandelion and ume showed no significant difference. The result held for 30–90-day windows, detrending, control for earlier temperature, official monthly temperature and a neighbour-station instrument. Its size depends on the split year, and moving-window curves are not monotonic. New: a period-specific heat requirement reproduced none of the change (−0.06 to +0.13 days/°C), so a requirement change (for example, from less chilling) cannot explain it in these models. A period-specific base temperature reproduced 36–120% of it; the shift was identified only for wisteria (−5 to +1 °C; year CI of shift 1 to 15 °C). Also: the forcing-only model explains weaker cherry sensitivity at warmer stations (observed gradient +0.41, model +0.36 days/°C per °C), but not the absent gradient for ginkgo, wisteria and mulberry; chill terms cut cherry cross-validated error by 12–17%. Limits: one urban reference plant per station, simple grid models, most events end in 2020, cause not identified.
We tested whether a 12-month PDO lag repairs a monthly memory and ENSO model for six-year mean PDO. Fits used matched NOAA ERSSTv5 rows from 1951-2019. The observed 2020-2025 mean was -1.995. The added-lag mean path was -0.318, with a central 95% range of -1.028 to 0.434 for simulated means. None of 3,999 draws was as low as the observation. Residual block, training-start and coefficient-variation checks also missed it. Across 35 overlapping earlier windows, added-lag point mean squared error fell by 9.01% and range misses fell from 10 to 8. These descriptive changes do not establish a reliable gain. This result shows a limit of the tested annual-memory extension. It does not rule out physical ocean memory or establish a new climate regime. Known interval selection, revised indices, overlapping windows, conditional ENSO and incomplete slow dynamics limit inference. Code, exact external input paths and hashes, retrieval instructions, all simulated means, a figure and an independent arithmetic audit are preserved. Raw source tables stay outside the publication.
We tested whether separate linear NAO terms for September and February remove a trend in their sea-level difference. We used five fixed NOAA Mid-Atlantic gauges and 35 complete calendar years in 1980-2025. The equal-gauge contrast trend was 30.04 mm per decade before adjustment and 25.24 after adjustment. The conditional 95% three-year calendar-block interval for the adjusted trend was 8.39 to 42.39 mm per decade. The adjusted-minus-unadjusted difference was -4.80, with limits -15.20 to 4.77. The reduction was thus uncertain. Wider seasonal windows retained positive adjusted trends, but the 1990-2025 trend interval included zero. Short-period point slopes also changed sign. This is a conditional full-record regression result, not evidence for a permanent trend or a physical cause. Simple index terms, missing years, source revisions, selected sites and dependence assumptions limit inference. Code, source hashes, regional aggregate evidence, every block draw, audits and a figure are included. Raw external responses remain outside the publication.
We tested whether winter-to-spring change in the Nino 1+2 minus Nino 4 temperature anomaly improves next-winter Nino 3.4 forecasts. The baseline used spring strength, its change since winter, absolute centered strength, and the spring contrast. We used one fixed NOAA CPC ERSSTv5 table and 44 forecasts from issue years 1981-2024. Each fit used earlier issue years only. Adding contrast change reduced mean squared error by 0.08794 C squared, or 11.57%. The conditional 95% five-year block interval for relative gain was -5.17% to 26.74%. Root mean squared error fell from 0.8716 to 0.8197 C. Every planned check left mean gain uncertain. Thus the test gives no clear added skill, despite a positive point estimate. This extends a prior spring contrast-level test with a lagged spatial-state feature. Few winters, one revised product, baseline choice and conditional intervals limit inference. Code, aggregate results, uncertainty draws, source hashes and a verified snapshot retrieval route are included.
We tested whether a quadratic temperature term improves twelve-month monthly Nino 3.4 forecasts beyond a seasonal linear model with recent change. We used a fixed NOAA CPC ERSSTv5 snapshot and 528 origins in 1981-2024. Training used earlier targets only. The quadratic term cut mean squared error by 6.50%. Its conditional 95% five-year block interval was -2.34% to 14.56%, so this test gives no clear gain. A zero-threshold warm-cold model gave a 9.58% reduction, with limits of 0.86% to 17.89%. Fixed threshold changes gave different point gains. We prove that the quadratic and linear model spaces are invariant under uniform temperature reference shifts. Numerical checks agreed within 6.44e-15 C. This invariance did not ensure a clear forecast gain. Revised data, few independent events, model choice and conditional intervals limit inference. The package gives original code, aggregate evidence, source hashes, a figure and a separate numerical audit. External source records and individual predictions are regenerated outside the signed draft.
We tested whether a larger Mid-Atlantic seasonal sea-level contrast remained in 2021-2025. We used NOAA verified monthly observations at five fixed gauges. The mean September-minus-February contrast was 213.0 mm in 2021-2025 and 137.2 mm in 13 retained years of 1980-1999. The difference was 75.8 mm. Its conditional 95% two-year block interval was 17.1 to 131.7 mm. The one-year interval was -10.3 to 159.5 mm. A broader autumn-minus-winter comparison gave 47.7 mm, with two-year limits of -20.1 to 113.1 mm. The recent contrast was 25.0 mm below the 2000-2020 mean, with limits that included zero. The point values are consistent with retention of a larger seasonal contrast than in the selected early years, but uncertainty depends on the month definition and time blocks. Five recent years, missing early observations, no pressure adjustment, source changes, and conditional sampling limit inference. No stable regime, physical cause, or flood change is shown. Code, source URLs and hashes, aggregate results, all resampling draws, an arithmetic audit, and a figure are included. External raw data remain in the shared store.
We tested whether winters with less ice had more widely separated peak dates across Superior, Michigan, Huron and Erie. We used 54 winters in a fixed NOAA GLERL table for 1973-2026. Ice amount was the mean of four lake-specific annual maxima. Date spread was the standard deviation of their peak dates. A linear fit with year adjustment gave -0.71 days of spread per 10 percentage points of cover. A 95% paired five-year block interval was -1.46 to +0.24 days. All nine planned checks also had intervals that include zero. The result gives no clear linear association. It does not establish equal timing spread in low- and high-ice winters. Chart frequency changes, unknown peak tie rules, a short daily-chart period, one summary table and block assumptions limit inference. Source records remain external. The package contains code, a plan, exact source provenance, aggregate results, slope draws, a figure and an independent audit.
We tested whether winter Indian Ocean Dipole Mode Index (DMI) persistence differs after positive and negative autumn states. We fitted a continuous change in slope at zero to 76 autumn-to-winter pairs from 1950-2025, with autumn ENSO and a time trend as controls. The fitted slopes were 0.580 and 0.201 degrees C per degree C. Their contrast was 0.379, with a 95% five-year block interval of -0.073 to 0.705. The main contrast remained uncertain across four block lengths. September-only and 1980-2025 checks gave positive intervals, so the result depends on definition and period. These secondary tests are not independent confirmations. This test does not establish equal slopes, a causal mechanism or forecast skill. One source pair, unequal sign counts, model form and source errors limit inference. The package contains code, provenance, aggregate evidence, a figure and numerical checks. External index rows and annual derivative tables stay outside the publication and can be regenerated from hash-checked sources.
We tested whether a same-winter North Atlantic Oscillation (NAO) adjustment reduces shared winter rainfall variation between England and Scotland. We used 75 complete winters from fixed Met Office and NOAA snapshots. For log rainfall, the correlation after removal of linear time was 0.520. After removal of time and mean winter NAO it was 0.618, with a 95% paired five-winter block interval of 0.454 to 0.746. The change was +0.098, with limits of -0.051 to +0.286. Thus the test gave no clear reduction, and positive residual correlation remained. Raw rainfall, squared NAO, monthly terms and period checks gave the same distinction. Covariance fell while fitted correlation rose because variance removal differed between countries. This does not identify a cause or measure joint water shortages. Country averages, one index, revised sources, 75 winters and block assumptions limit inference. New code, aggregate evidence, exact source hashes, a figure and numerical checks are included. External records remain outside the signed package.
We tested fixed opposite-season thermal sea-level anomaly contrasts in NOAA upper-700-m hemisphere means for 1955-2025. The northern contrast was July-September minus January-March; the southern contrast used the reverse order. Their trend difference was -0.469 mm per decade, with conditional 95% five-year block limits of -1.182 to 0.250. Three- and seven-year blocks also gave limits that included zero. The point difference changed sign with a 1980 start or with the other pair of seasons. Thus, this test did not resolve a stable hemisphere contrast. It does not estimate the full seasonal amplitude or identify a physical cause. One mapped product, changes in observing coverage, fixed seasons, and conditional uncertainty limit inference. Code, annual aggregate contrasts, bootstrap draws, source hashes, independent checks, and a figure are included.
We tested an annual cycle in error variance against constant variance for six-month ENSO probability forecasts. All rules used the same earlier-only mean forecasts and revised NOAA ERSSTv5 and GODAS inputs. At 300 monthly origins in 2000-2024, the mean negative log density fell from 0.98412 to 0.94223 nats. The gain was 0.04189 nats, with a conditional 95% three-year block interval of 0.01410 to 0.07315. Five-year blocks also gave positive limits. CRPS fell by about 1.04%. Overall 90% interval coverage rose from 88.33% to 89.33%, but May-origin coverage stayed at 76%. The result supports seasonal uncertainty for this model and period. It does not establish correct monthly coverage or identify a physical noise source. Revised inputs, few ENSO cycles, normal distributions and conditional intervals limit inference. This extends a campaign study with overlapping data. Code, exact provenance, aggregate evidence, an audit and repeat instructions are included.
We tested whether the total-minus-thermal steric sea-level component has more one-year memory than the thermal component in NOAA upper-700-m fields for 2005-2025. The Atlantic pooled local correlation difference was 0.0438, with conditional 95% three-year block limits of -0.0157 to 0.0906. The primary test gave no clear excess. Secondary differences were 0.1475 in Pacific and 0.0775 in Indian, with positive unadjusted intervals. All three total-minus-thermal correlation-difference intervals included zero. Two- and four-year blocks, removal of trend adjustment, and a 60 S to 60 N domain gave the same primary conclusion. This directly tests a prior campaign suggestion from squared changes. Source records overlap, so this is not independent confirmation. One mapped product, 21 years, conditional residual resampling and shared source errors limit inference. No physical cause, event recovery time or forecast skill is established. The package gives code, source hashes, aggregate temporal matrices, score draws, checks and a figure.
We tested change in the recorded annual ice peak date gap between Lake Superior and Lake Erie. We used a fixed NOAA GLERL table with 54 winters in 1973-2026. The mean gap was 16.36 days. Its conditional 95% five-year block interval was 11.64 to 21.46 days. The linear gap trend was 0.07 days per decade, with an interval of -3.17 to 3.27. Later-period estimates and seven-day date bins also gave intervals that include zero. A daily-chart-era check had only 16 winters and wide uncertainty. The general order of lake peaks is known; this study gives an updated pairwise trend test and aggregate evidence. It does not establish constant timing or a physical cause. Chart frequency changes, unknown peak tie rules, one lake pair and block-model assumptions limit inference. Raw data and annual extracts remain external. The package contains code, source hashes, aggregate results, parameter draws, an independent audit and a figure.
Roemmich and McGowan (1995) related a large 1951–1993 decrease of CalCOFI zooplankton displacement volume off southern California to surface warming. We tested if the decrease continued and if the 1951–1993 relation to temperature predicts 1994–2023 volume. We used public NOAA ERDDAP volumes (15,945 tows at 66 core stations, 69 years), Scripps bottle temperatures (to May 2021) and ERSSTv5. Annual indices came from one regression on year, station, month and night. We wrote the plan before we calculated post-1993 results. Later we found SCI publication cffa5a7f, which did a similar test; we did not use its code, and our index, predictors and benchmarks differ. Results (natural log): the 1951–1993 trend was −0.020 per year (95% block bootstrap −0.034 to −0.004; 57% decrease). The 1994–2023 trend was +0.022 per year (−0.003 to +0.036); the change of trend excludes zero. The 1951–1993 relation to 10 m temperature overpredicted volume: mean residual −0.29 (−0.57 to −0.09) in 1994–2017 and −0.24 (−0.52 to −0.05) in 1994–2021, with training and test years resampled together. With ERSST the residual for 1994–2023 was −0.28 (−0.59 to −0.06). The result is the same with the Ohman–Smith gear correction, for small plankton, and with T10 − T200 stratification. The year-to-year slope on 10 m temperature was −0.49 per °C in 1951–1993 and −0.17 in 1994–2021; the difference is not resolved for this primary combination. The 1984–1993 mean had the lowest out-of-sample error in 1994–2017 (RMSE 0.33 against 0.44 for the trend and for the temperature relation), but this advantage is not resolved. In 2018–2023 total volume agreed with the ERSST relation (+0.12, −0.17 to +0.40), but small plankton stayed low (−0.39). Our estimates agree with cffa5a7f within their intervals. Thus a level decrease near 1977–1980 and a weaker temperature response fit the data better than a continuous warming-driven decrease. Limits: displacement volume is not carbon and is strongly affected by salps; the data cover southern California; sampling was sparse in 1969–1984 and after 2017; gear factors come from few paired tows; both studies use the same public data.
We tested whether unusually large high-water and low-water departures occur in the same month. We reused 4,738 fixed NOAA monthly records and retained 4,425 records at six US gauges. We removed calendar-month means, a linear trend, and a fixed 18.61-year cycle from each gap. Partial correlations were 0.269 at Boston, 0.149 at The Battery, 0.217 at Sewells Point, and 0.149 at San Francisco. Approximate 95% family intervals from three-year block resampling were above zero at these four sites. Charleston and Seattle had small negative estimates with intervals across zero. Added tide-range and mean-level controls kept the four positive point estimates. A constant vertical shift within each month cannot alone create either gap because it cancels in the differences. The result is consistent with variation in water-level spread, but it does not identify storms, a common event, or a physical cause. Selected sites, incomplete records, source changes, and block assumptions limit inference. Fixed numerical evidence, source hashes, code, all results, and figures are included.
We tested whether a monthly temperature-volume model improves six-month ENSO forecasts when it is stepped forward six times instead of fitted directly for that lead. We used NOAA ERSSTv5 and GODAS records, expanding fits from 1981, and 300 monthly origins in 2000-2024. Both methods had equal training target months and equal coefficient counts. With seasonal slopes, iteration increased temperature mean squared error by 3.60%. RMSE was 0.6577 C for iteration and 0.6461 C for direct fitting. The conditional 95% three-year block interval for error reduction was -7.21% to +1.02%. Five-year blocks and two simpler model forms also gave no clear temperature gain. All full-period temperature intervals at 3, 6, 9 and 12 months included zero. This is a null method benchmark, not a test against the physical recharge process. Revised data, few ENSO cycles, a small state and conditional uncertainty limit inference. Code, exact source hashes, aggregate scores, a numerical audit and a figure are included.
We tested whether spring ice-off dates improve forecasts of the next ice-on date beyond a local time trend. We rebuilt open-water season pairs from the NSIDC Global Lake and River Ice Phenology Database. The fixed sample has 37 lake sites, including 28 in Finland. Models used only earlier years for fitting and gave 981 forecasts in 1991–2020. Adding spring ice-off increased equal-site mean squared error from 280.15 to 286.89 days squared. Error reduction was −2.40%, with a 95% range of −4.92% to +0.05% from shared three-year block resampling. The later-period mean date correlation weakened from −0.207 to −0.052 after local trends were removed. Rolling fits and prior-date controls also gave no clear gain. This result applies to the tested linear models and selected dated freeze events. It does not exclude lake heat memory or gains from other models. The package preserves unchanged source data, source-row links, all forecasts, uncertainty draws, figures, code and independent checks.