Please see the excellent upgrade deepsphere-weather, a scalable and flexible framework to apply convolutions on spherical unstructured grids for weather/climate applications.

Geometric deep learning for medium-range weather prediction

Wentao Feng, Michaël Defferrard, Gionata Ghiggi, Icíar Lloréns Jover, Natalie Bolón Brun

The code in this repository provides a framework for a deep learning medium range weather prediction method based on graph spherical convolutions.

[June 2020]: The results obtained with this code are detailed in the Masters thesis report and slides.

[September 2020]: Results have been improved from the initial basis thanks to:

Introduction of residual connections in the architecture
Inclusion of further consecutive steps in the loss with different weighting schemes to reduce the loss at long term predictions

Model	Z500 (6h)	t850 (6h)	Z500 (120h)	t850 (120h)
Weyn et al	103.17	1.0380	611.33	2.957
Iciar June 2020	67.46	0.7172	861.7	3.432
Ours Sep 2020	61.58	0.7110	680.024	2.901

Results can be checked at plot_results.ipynb

Resources:

Installation

For a local installation, follow the below instructions.

Clone this repository.

git clone git@github.com:ownzonefeng/weather_prediction.git
cd weather_prediction

Install the dependencies.
```
conda env create -f environment.yml
```

Create the data folders

mkdir data/equiangular/5.625deg/ data/healpix/5.625deg/

Download the WeatherBench data on the data/equiangular/5.625deg/ folder by following instructions on the WeatherBench repository.
Interpolate the WeatherBench data onto the HEALPix grid. Modify the paremeters in scripts/config_data_interpolation.yml as desired.
```
python -m scripts.data_iterpolation -c scripts/config_data_interpolation.yml
```

Attention:

If it does not find the module SphereHealpix from pygsp, install the development branch using:

conda activate weather_modelling
pip install git+https://github.com/Droxef/pygsp@new_sphere_graph

Reproduce latest results

Train model:

The model listed as "Ours 2020" is trained using the module full_pipeline_multiple_steps.py. An example of how to use it can be found on the notebook Restarting_weights_per_epoch.ipynb.

The config file to be used is configs/config_residual_multiple_steps.json. You may want to modify the model name and/or the data paths if the data has been relocated.

Evaluate model:

You can generate the model predictions using the notebook generate_evaluate_predictions.ipynb. The parameters to be modified are:

model name (third cell)
epochs to be evaluated (you can define a range or a single one)

In order to evaluate the performance of the model, you only need to run up to "Generate plots for evaluation". This sencond part will generate the skill and climatology plots (you may be interested in generate them for a single epoch usually, not all of them)

Compare models:

In order to compare the performance of different models, or the same model at different epochs or simply a model against different baselines, you can use the notebook plot_results.ipynb. Depending on the purpose of the comparison, you may want to run a different section of the notebook. An explanation of each section and its use case can be found under the heading of the notebook.

Modules

full_pipeline_evaluation.py

Allows to train, test, generate predictions and evaluate them for a model trained with a loss function that includes 2 steps. All parameters, except GPU configuration, are defined in a config file such as the ones stored on the folder configs/ .

full_pipeline_multiple_steps.py

Allows to to train and test a model
with a loss function that includes multiple steps that can be defined by the user. It saves the model after every epoch but does not generate the predictions (to save time since it can be done in parallel using the notebook generate_evaluate_predictions.ipynb ). The parameters are defined inside the main function, although it can be adapted to use a config file as in full_pipeline_evalution.py

It is important to remark that the update function that takes care of the weight's update is defined on top of the file and should be adapted to the number of lead steps taken into account in the loss function.

architecture.py

Contains pytorch models used for both full_pipeline_multiple_steps.py and full_pipeline_evaluation.py Previous architectures used can be found in the folder modules/old_architectures/

plotting.py

Contains different functions to generate evaluation plots.

train_last_model.py

Contains code to train model with 2step-ahead prediction such as the one used for Iciar2020 results.

Notebooks

The main notebooks to explore are:

Restarting_weights_per_epoch.ipynb Contains an example of how to use the functions that train the model that reported the best results mentioned earlier.
generate_evaluate_predictions.ipynb Generate values on validation set using the weights of the desired saved model
plot_results.ipynb Generate loss plots and comparison plots against different benchmark models
healpix_resampling.ipynb Generate healpix data from equiangular data
generate_observations.ipynb Generate ground-truth data for evaluation of the models

License

The content of this repository is released under the terms of the MIT license.