Atmospheric forcing for all regional ocean simulations, in two forms: the raw CMIP6 archive (each model’s own published output, unmodified) and the bias-corrected version of the same fields (quantile delta mapping calibrated against ERA5 reanalysis over the historical period). The raw section below comes first – it’s the reference the bias correction is calibrated against; everything after the divider is the bias-corrected side.

Raw CMIP6 Data

Meteorological forcing fields exactly as published by the source model – no bias correction applied. This is the reference the bias-corrected fields further down this page are calibrated against, fetched directly from Pangeo/GCS or ESGF (see ocean-data’s DataLoader) rather than through the bias-correction pipeline.

GFDL-ESM4 decadal trends

GFDL-ESM4, domain-averaged decadal means. Solid black = historical (1990–2014); dashed = ssp126/ssp370 (2015–2099). The dotted vertical line at 2015 marks where the underlying CMIP6 experiment switches from the model’s own historical run to SSP-forced projections – a real change in data source, not a plotting artifact.

3-hourly source availability

Which core variables actually publish 3-hourly output at all (checked directly against Pangeo and ESGF), independent of whether they’ve been fetched into the raw archive yet:

VariableGFDL-ESM4MPI-ESM1-2-HRCNRM-ESM2-1 (hist)CNRM-ESM2-1 (ssp126)CNRM-ESM2-1 (ssp370)
tas✅✅✅✅✅
pr✅✅✅✅✅
uas✅✅✅❌❌
vas✅✅✅❌❌
huss✅✅✅✅✅
rsds✅❌✅✅❌
rlds✅❌✅✅❌
rsus✅❌✅❌❌
rlus✅❌✅❌❌
ps✅❌✅❌❌

GFDL-ESM4 publishes all 10 at 3-hourly for every period and is the only model fetched into the raw archive so far. MPI-ESM1-2-HR is capped at 5/10 (no radiation or surface-pressure fields at 3-hourly, in any period). CNRM-ESM2-1’s historical run has full coverage, but its ScenarioMIP submission is considerably sparser – ssp370 publishes only tas/pr/huss at 3-hourly.


Bias-Corrected CMIP6 Data

Diagnostics for the bias correction of CMIP6 atmospheric forcing fields applied to all regional ocean simulations. Results are scenario-specific and valid across all model domains.

Bias correction is applied using quantile delta mapping (QDM) calibrated against ERA5 reanalysis over the historical period. The plots below show calibration performance and projected trends under each SSP scenario.

SSP Scenario Reference

ScenarioNarrativeRadiative forcingApprox. warming by 2100
SSP1-1.9Sustainability — very low emissions1.9 W m⁻²1.0–1.8 °C
SSP1-2.6Sustainability — low emissions2.6 W m⁻²1.3–2.4 °C
SSP2-4.5Middle of the road4.5 W m⁻²2.1–3.5 °C
SSP3-7.0Regional rivalry — high emissions7.0 W m⁻²2.8–4.6 °C
SSP5-8.5Fossil-fuelled development8.5 W m⁻²3.3–5.7 °C

Warming ranges from IPCC AR6 (2021), relative to 1850–1900.

Variable reference

Short nameFull name
evspsblEvaporation (evspsbl)
hussNear-surface Specific Humidity (huss)
net_lwNet Longwave Radiation (net_lw)
net_swNet Shortwave Radiation (net_sw)
peNet Freshwater Flux, P-E (pe)
prPrecipitation (pr)
pslSea-level Pressure (psl)
rldsDownwelling Longwave Radiation (rlds)
rlusUpwelling Longwave Radiation (rlus)
rsdsDownwelling Shortwave Radiation (rsds)
rsusUpwelling Shortwave Radiation (rsus)
tasNear-surface Air Temperature (tas)
uasEastward Wind (uas)
vasNorthward Wind (vas)

Variable availability

+ = bias-corrected data available, − = not available. Each model’s three columns are historical / SSP1-2.6 / SSP3-7.0, in that order.

VariableCNRM-ESM2-1GFDL-ESM4MPI-ESM1-2-HR
evspsbl+ + ++ + ++ + +
huss+ + ++ + ++ + +
net_lw+ − [4] ++ + ++ + +
net_sw+ − − [3]+ + ++ + +
pe+ + +− − − [5]+ + +
pr+ + ++ + ++ + +
psl+ + ++ + ++ + +
rlds+ − [2] ++ + ++ + +
rlus+ + ++ + ++ + +
rsds+ + ++ + ++ + +
rsus+ − − [1]+ + ++ + +
tas+ + ++ + ++ + +
uas+ + ++ + ++ + +
vas+ + ++ + ++ + +

Notes

  1. rsus (upwelling shortwave) is confirmed permanently absent from the CMIP6 archive for CNRM-ESM2-1 – checked across all 5 ESGF Solr mirrors, the Globus ESGF2-US-1.5 catalog, and Pangeo, at both 3-hourly and daily frequency, both scenarios: zero results everywhere. Not a network issue, not retried further.
  2. rlds (downwelling longwave) for CNRM-ESM2-1/SSP1-2.6 specifically is confirmed genuinely blocked – all 3 access paths exhausted (ORNL replica 403s, CNRM’s own node unreachable, Pangeo has no daily-frequency rlds for this model/experiment). SSP3-7.0 is unaffected and has real rlds data.
  3. net_sw = rsds − rsus; blocked because rsus is permanently unavailable for CNRM-ESM2-1 (note 1). The pipeline’s pseudo_tcc radiation path (deriving cloud fraction from rsds alone) is used for CNRM instead of net_sw/net_lw/components.
  4. net_lw = rlds − rlus; blocked for SSP1-2.6 specifically because rlds is blocked there (note 2) – SSP3-7.0 has both ingredients and a real net_lw.
  5. pe (bias-corrected P-E composite) is skipped entirely for GFDL-ESM4: its CMIP6 hfls (needed to derive evspsbl from latent heat flux) is only published at monthly (Amon) frequency, not daily, for this model – a structural limitation, not a download failure. river_projection.py falls back to separately-corrected pr/evspsbl for GFDL runs instead.

Summary

Calibration quality — RMSE and Correlation heatmap across models and variables Calibration quality: RMSE (left) and Correlation (right) for each model and variable against ERA5. Darker red = higher RMSE; darker green = higher correlation.

Future mean by model and scenario Future area-mean of the bias-corrected field for each model and SSP scenario. Black tick marks show the ERA5 calibration period mean.

Future period-mean trajectory Area-mean of the bias-corrected field for each ~decade window within the future period, per model and SSP scenario — shows how the projected change develops over time rather than just its end-state value.

Calibration statistics

RMSE

VarCNRM-ESM2-1GFDL-ESM4MPI-ESM1-2-HR
evspsbl0.00010.00010.0001
huss0.00210.00220.0022
net_lw16.538017.527017.7580
net_sw46.200042.214042.9430
pe0.0000—0.0000
pr0.00000.00000.0000
psl644.1400630.9100687.8000
rlds27.411032.435027.7440
rlus37.346042.491036.4010
rsds49.174045.920046.0560
rsus—16.783016.4710
tas6.90726.94257.1517
uas1.94492.18991.9771
vas1.25341.12871.0989

Bias

VarCNRM-ESM2-1GFDL-ESM4MPI-ESM1-2-HR
evspsbl-0.0001-0.0001-0.0001
huss+0.0002-0.0001-0.0001
net_lw-3.5944+5.4586+7.0397
net_sw+13.3820-8.3224-5.2298
pe-0.0000—-0.0000
pr+0.0000+0.0000+0.0000
psl-11.1920-66.8880-47.0640
rlds+2.2628-0.3260+4.7042
rlus+5.8572-5.7845-2.3355
rsds+15.0320-5.8594-6.8518
rsus—+2.4631-1.6219
tas+0.3457-1.6123-1.1889
uas-0.5071-0.5832-0.0971
vas-0.1669-0.0121-0.1706

Correlation

VarCNRM-ESM2-1GFDL-ESM4MPI-ESM1-2-HR
evspsbl0.30290.30550.2767
huss0.51900.51480.4797
net_lw0.20750.16570.1869
net_sw0.84790.83460.8237
pe0.0330—0.0323
pr0.40760.49090.4790
psl0.0360-0.0467-0.1794
rlds0.56710.50840.5468
rlus0.47300.41690.4384
rsds0.86000.83670.8359
rsus—0.40070.3147
tas0.50130.51010.4723
uas-0.2500-0.2875-0.1838
vas0.40650.38870.4043

Future mean by model and scenario

CNRM-ESM2-1

VarhistoricalSSP1-2.6SSP3-7.0
evspsbl8.216e-057.686e-058.532e-05
huss0.006010.0059540.006443
net_lw-52.59—-50.16
net_sw103.7——
pe—1.077e-051.159e-05
pr3.626e-053.567e-053.893e-05
psl1.024e+051.024e+051.024e+05
rlds——311.8
rlus—1.284e+061.313e+06
rsds—118.4116.2
rsus———
tas284.6284.6285.9
uas1.1021.3511.402
vas0.54860.47370.5904

GFDL-ESM4

VarhistoricalSSP1-2.6SSP3-7.0
evspsbl7.83e-058.995e-058.954e-05
huss0.0060580.0062450.006724
net_lw-52.4-53.17-51.92
net_sw104.1106.6104.9
pe———
pr3.505e-053.79e-053.986e-05
psl1.023e+051.024e+051.024e+05
rlds—307.9313.4
rlus—363.3371.4
rsds—121.5118.9
rsus—16.7614.78
tas284.7285.5286.8
uas1.1441.3241.35
vas0.6690.54550.6379

MPI-ESM1-2-HR

VarhistoricalSSP1-2.6SSP3-7.0
evspsbl8.215e-059.337e-059.026e-05
huss0.0061540.0063770.006704
net_lw-53.15-53.32-52.3
net_sw104.9106.2105
pe—1.038e-051.096e-05
pr3.531e-053.786e-053.91e-05
psl1.02e+051.02e+051.02e+05
rlds—309.1313.5
rlus—367.41.344e+06
rsds—120.4118.8
rsus—15.625.326e+04
tas9.052287.4288.2
uas0.8691.2231.284
vas0.56150.50960.5673

Model Pages

Bias correction diagnostics and river flow projections for each model and scenario: