CoolFace
Datasetpublic

LEAP/ChaosBench

ChaosBench: A Multi-Channel, Physics-Based Benchmark for Subseasonal-to-Seasonal Climate Prediction NeurIPS 2024 Oral ChaosBench is a benchmark project to improve and extend the predictability range of deep weather emulators to the subseasonal-to-seasonal (S2S) range. Predictability at this scale is more challenging due to its: (1) double sensitivities to intial condition (in weather-scale) and boundary condition (in climate-scale), (2) butterfly effect, and… See the full description on the dataset page: https://huggingface.co/datasets/LEAP/ChaosBench.

sourceHugging Facegpl-3.0updated 1y agoView on Hugging Face
7likes5.8kdownloads
Dataset Card

ChaosBench: A Multi-Channel, Physics-Based Benchmark for Subseasonal-to-Seasonal Climate Prediction

<div align="center"> <h3>NeurIPS 2024 Oral</h3> </div>

<div align="center" style="display: flex; justify-content: center; gap: 10px;"> <a href="https://leap-stc.github.io/ChaosBench"><img src="https://img.shields.io/badge/View-Documentation-blue?style=for-the-badge" alt="Homepage"/></a> <a href="https://arxiv.org/abs/2402.00712"><img src="https://img.shields.io/badge/ArXiV-2402.00712-b31b1b.svg?style=for-the-badge" alt="arXiv"/></a> <a href="https://huggingface.co/datasets/LEAP/ChaosBench"><img src="https://img.shields.io/badge/Dataset-HuggingFace-ffd21e?style=for-the-badge" alt="Huggingface Dataset"/></a> <a href="https://github.com/leap-stc/ChaosBench/blob/main/LICENSE"><img src="https://img.shields.io/badge/License-GNU%20GPL-green?style=for-the-badge" alt="License Badge"/></a> </div>

ChaosBench is a benchmark project to improve and extend the predictability range of deep weather emulators to the subseasonal-to-seasonal (S2S) range. Predictability at this scale is more challenging due to its: (1) _double sensitivities to intial condition (in weather-scale) and boundary condition (in climate-scale), (2) butterfly effect, and our (3) inherent lack of understanding of physical processes operating at this scale. Thus, given the high socioeconomic stakes_ for accurate, reliable, and stable S2S forecasts (e.g., for disaster/extremes preparedness), this benchmark is timely for DL-accelerated solutions.

🪜 Leaderboard

View our interactive leaderboard dashboard here

You can access the full score and checkpoints in logs/<MODEL_NAME> within the following subdirectory:

  • Scores: eval/<METRIC>.csv
  • Model checkpoints: lightning_logs/

✨ Features

1️⃣ _Diverse Observations_. Spanning over 45 years (1979-), we include ERA5/LRA5/ORAS5 reanalysis for a fully-coupled Earth system emulation (atmosphere-terrestrial-sea-ice)

2️⃣ _Diverse Baselines_. Wide selection of physics-based forecasts from leading national agencies in Europe, the UK, America, and Asia

3️⃣ _Differentiable Physics Metrics_. Introduces two differentiable physics-based metrics to minimize the decay of power spectra at long forecasting horizon (blurriness)

4️⃣ _Large-Scale Benchmarking_. Systematic evaluation (deterministic, probabilistic, physics-based) for state-of-the-art ML-based weather emulators like ViT/ClimaX, PanguWeather, GraphCast, and FourcastNetV2

🏁 Getting Started

_NOTE:_ Only need the dataset? Jump directly to Step 2. If you find any problems, feel free to contact us or raise a GitHub issue.

Step 0: Clone the ChaosBench Github repository

Step 1: Install package dependencies

$ cd ChaosBench
$ pip install -r requirements.txt

Step 2: Initialize the data space by running

$ cd data/
$ wget https://huggingface.co/datasets/LEAP/ChaosBench/resolve/main/process.sh
$ chmod +x process.sh

Step 3: Download the data

# Required for inputs and climatology (e.g., for normalization; 1979-)
$ ./process.sh era5
$ ./process.sh lra5
$ ./process.sh oras5
$ ./process.sh climatology

# Optional: control (deterministic) forecasts (2018-)
$ ./process.sh ukmo
$ ./process.sh ncep
$ ./process.sh cma
$ ./process.sh ecmwf

# Optional: perturbed (ensemble) forecasts (2022-)
$ ./process.sh ukmo_ensemble
$ ./process.sh ncep_ensemble
$ ./process.sh cma_ensemble
$ ./process.sh ecmwf_ensemble

# Optional: state-of-the-art (deterministic) forecasts (2022-)
$ ./process.sh panguweather
$ ./process.sh fourcastnetv2
$ ./process.sh graphcast

🔍 Dataset Overview

All data has daily and 1.5-degree resolution.

  1. 1._ERA5 Reanalysis for Surface-Atmosphere (1979-2023). The following table indicates the 48 variables (channels) that are available for Physics-based models. Note that the Input ERA5 observations contains ALL_ fields, including the unchecked boxes:
Parameters/Levels (hPa)10009258507005003002001005010
Geopotential height, z ($gpm$)&check;&check;&check;&check;&check;&check;&check;&check;&check;&check;
Specific humidity, q ($kg kg^{-1}$)&check;&check;&check;&check;&check;&check;&check;&nbsp;&nbsp;&nbsp;
Temperature, t ($K$)&check;&check;&check;&check;&check;&check;&check;&check;&check;&check;&check;
U component of wind, u ($ms^{-1}$)&check;&check;&check;&check;&check;&check;&check;&check;&check;&check;
V component of wind, v ($ms^{-1}$)&check;&check;&check;&check;&check;&check;&check;&check;&check;&check;
Vertical velocity, w ($Pas^{-1}$)&nbsp;&nbsp;&nbsp;&nbsp;&check;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
  1. 1._LRA5 Reanalysis_ for Terrestrial (1979-2023)
AcronymsLong NameUnits
asnsnow albedo(0 - 1)
d2m2-meter dewpoint temperatureK
etotal evaporationm of water equivalent
essnow evaporationm of water equivalent
evabsevaporation from bare soilm of water equivalent
evaowevaporation from open waterm of water equivalent
evatcevaporation from top of canopym of water equivalent
evavtevaporation from vegetation transpirationm of water equivalent
falforecaste albedo(0 - 1)
lai\_hvleaf area index, high vegetation$m^2 m^{-2}$
lai\_lvleaf area index, low vegetation$m^2 m^{-2}$
pevpotential evaporationm
rorunoffm
rsnsnow density$kg m^{-3}$
sdsnow depthm of water equivalent
sdesnow depth water equivalentm
sfsnowfallm of water equivalent
sktskin temperatureK
slhfsurface latent heat flux$J m^{-2}$
smltsnowmeltm of water equivalent
snowcsnowcover\%
spsurface pressurePa
srcskin reservoir contentm of water equivalent
srosurface runoffm
sshfsurface sensible heat flux$J m^{-2}$
ssrnet solar radiation$J m^{-2}$
ssrddownload solar radiation$J m^{-2}$
ssrosub-surface runoffm
stl1soil temperature level 1K
stl2soil temperature level 2K
stl3soil temperature level 3K
stl4soil temperature level 4K
strnet thermal radiation$J m^{-2}$
strddownward thermal radiation$J m^{-2}$
swvl1volumetric soil water layer 1$m^3 m^{-3}$
swvl2volumetric soil water layer 2$m^3 m^{-3}$
swvl3volumetric soil water layer 3$m^3 m^{-3}$
swvl4volumetric soil water layer 4$m^3 m^{-3}$
t2m2-meter temperatureK
tptotal precipitationm
tsntemperature of snow layerK
u1010-meter u-wind$ms^{-1}$
v1010-meter v-wind$ms^{-1}$
  1. 1._ORAS Reanalysis_ for Sea-Ice (1979-2023)
AcronymsLong NameUnits
iicethicsea ice thicknessm
iicevelusea ice zonal velocity$ms^{-1}$
iicevelvsea ice meridional velocity$ms^{-1}$
ileadfrasea ice concentration(0-1)
so14chgtdepth of 14$^\circ$ isothermm
so17chgtdepth of 17$^\circ$ isothermm
so20chgtdepth of 20$^\circ$ isothermm
so26chgtdepth of 26$^\circ$ isothermm
so28chgtdepth of 28$^\circ$ isothermm
sohefldonet downward heat flux$W m^{-2}$
sohtc300heat content at upper 300m$J m^{-2}$
sohtc700heat content at upper 700m$J m^{-2}$
sohtcbtmheat content for total water column$J m^{-2}$
sometauymeridonial wind stress$N m^{-2}$
somxl010mixed layer depth 0.01m
somxl030mixed layer depth 0.03m
sosalinesalinityPractical Salinity Units (PSU)
sossheigsea surface heightm
sosstsstsea surface temperature$^\circ C$
sowaflupnet upward water flux$kg/m^2/s$
sozotauxzonal wind stress$N m^{-2}$

💡 Baseline Models

In addition to climatology and persistence, we evaluate the following:

  1. 1._Physics-based models (including control/perturbed forecasts)_:
  2. 2.[x] UKMO: UK Meteorological Office
  3. 3.[x] NCEP: National Centers for Environmental Prediction
  4. 4.[x] CMA: China Meteorological Administration
  5. 5.[x] ECMWF: European Centre for Medium-Range Weather Forecasts
  6. 6._Data-driven models_:
  7. 7.[x] Lagged-Autoencoder
  8. 8.[x] Fourier Neural Operator (FNO)
  9. 9.[x] ResNet
  10. 10.[x] UNet
  11. 11.[x] ViT/ClimaX
  12. 12.[x] PanguWeather
  13. 13.[x] GraphCast
  14. 14.[x] Fourcastnetv2

🏅 Evaluation Metrics

We divide our metrics into 3 classes: (1) Deterministic-based, which cover evaluation used in conventional deterministic forecasting tasks, (2) Physics-based, which are aimed to construct a more physically-faithful and explainable data-driven forecast, and (3) Probabilistic-based, which account for the skillfulness of ensemble forecasts.

  1. 1._Deterministic-based:_
  2. 2.[x] RMSE
  3. 3.[x] Bias
  4. 4.[x] Anomaly Correlation Coefficient (ACC)
  5. 5.[x] Multiscale Structural Similarity Index (MS-SSIM)
  6. 6._Physics-based:_
  7. 7.[x] Spectral Divergence (SpecDiv)
  8. 8.[x] Spectral Residual (SpecRes)
  1. 1._Probabilistic-based:_
  2. 2.[x] RMSE Ensemble
  3. 3.[x] Bias Ensemble
  4. 4.[x] ACC Ensemble
  5. 5.[x] MS-SSIM Ensemble
  6. 6.[x] SpecDiv Ensemble
  7. 7.[x] SpecRes Ensemble
  8. 8.[x] Continuous Ranked Probability Score (CRPS)
  9. 9.[x] Continuous Ranked Probability Skill Score (CRPSS)
  10. 10.[x] Spread
  11. 11.[x] Spread/Skill Ratio