Add recipe for seasonal cycle and time series of Arctic/Antarctic sea ice area (REF) (#3891)

Co-authored-by: Manuel Schlund <[email protected]>

Author: axel-lauer

doc/sphinx/source/recipes/figures/ref/annual_cycle_sea_ice_area_nh_ambiguous_dataset_ambiguous_mip_historical_r1i1p1f1.png

@@ -25,6 +25,9 @@ Recipes are stored in recipes
   Maps and zonal means of longwave and shortwave cloud radiative effect
 * ref/recipe_ref_scatterplot.yml:
   2D histograms with focus on clouds
+* ref/recipe_ref_sea_ice_area_basic.yml:
+  Seasonal cycle of Arctic (NH) and Antarctic (SH) sea ice area, time series
+  of Arctic September (NH) and Antarctic February (SH) sea ice area
 
 
 Example plots:
@@ -45,9 +48,24 @@ Example plots:
    effect from CERES-EBAF Ed4.2 (solid black), ESACCI-CLOUD (dashed black),
    ISCCP-FH (dotted black) and the MPI-ESM1-2-LR model (blue).
 
-.. _fig_ref_3:
 .. figure::  /recipes/figures/ref/seaborn_jointplot.png
    :align:   center
 
    2D histogram of total cloud fraction (ESACCI-CLOUD) and shortwave cloud radiative
    effect (CERES-EBAF) for the years 2001-2016 with 1D histograms attached.
+
+.. _fig_ref_4:
+.. figure::  /recipes/figures/ref/annual_cycle_sea_ice_area_nh_ambiguous_dataset_ambiguous_mip_historical_r1i1p1f1.png
+   :align:   center
+   :width:   8cm
+
+   Average seasonal cycle of the Arctic (NH) sea ice area from MPI-ESM1-2-LR
+   (red line) compared with OSISAF/CCI (blue line). Created with recipe_ref_sea_ice_area_basic.yml.
+
+.. _fig_ref_5:
+.. figure::  /recipes/figures/ref/timeseries_sea_ice_area_nh_sep_ambiguous_dataset_ambiguous_mip_historical_r1i1p1f1.png
+   :align:   center
+   :width:   8cm
+
+   Time series of Arctic (NH) September (NH) sea ice area from MPI-ESM1-2-LR
+   (red line) compared with OSISAF/CCI (blue line). Created with recipe_ref_sea_ice_area_basic.yml.

doc/sphinx/source/recipes/figures/ref/timeseries_sea_ice_area_nh_sep_ambiguous_dataset_ambiguous_mip_historical_r1i1p1f1.png

@@ -0,0 +1,154 @@
+# ESMValTool
+---
+documentation:
+  title: Seasonal cycle of NH/SH sea ice area
+  description: >
+    Plots of the seasonal cycle and time series of NH and SH sea ice area.
+  authors:
+    - lauer_axel
+  maintainer:
+    - lauer_axel
+
+
+datasets:
+  - {dataset: MPI-ESM1-2-LR, grid: gn, benchmark_dataset: true}  # example model
+
+preprocessors:
+
+  pp_nh:
+    extract_region:
+      start_longitude: 0
+      end_longitude: 360
+      start_latitude: 30
+      end_latitude: 90
+    climate_statistics:
+      period: month
+    area_statistics:
+      operator: sum
+    convert_units:
+      units: 1e6 km2
+
+  pp_nh_sep:
+    extract_region:
+      start_longitude: 0
+      end_longitude: 360
+      start_latitude: 30
+      end_latitude: 90
+    extract_month:
+      month: 9
+    area_statistics:
+      operator: sum
+    convert_units:
+      units: 1e6 km2
+
+  pp_sh:
+    extract_region:
+      start_longitude: 0
+      end_longitude: 360
+      start_latitude: -90
+      end_latitude: -30
+    climate_statistics:
+      period: month
+    area_statistics:
+      operator: sum
+    convert_units:
+      units: 1e6 km2
+
+  pp_sh_feb:
+    extract_region:
+      start_longitude: 0
+      end_longitude: 360
+      start_latitude: -90
+      end_latitude: -30
+    extract_month:
+      month: 2
+    area_statistics:
+      operator: sum
+    convert_units:
+      units: 1e6 km2
+
+
+diagnostics:
+
+  siarea_seas:
+    description: Create seasonal cycle plot of sea ice area.
+    variables:
+      sea_ice_area_nh: &var_basic
+        short_name: siconc
+        timerange: '1995/2014'
+        preprocessor: pp_nh
+        project: CMIP6
+        mip: SImon
+        exp: historical
+        ensemble: r1i1p1f1
+        grid: gr
+        reference_dataset: OSI-450-nh
+        additional_datasets:
+        - {dataset: OSI-450-nh, project: OBS, type: reanaly, version: v3, mip: OImon, tier: 2,
+           supplementary_variables: [{short_name: areacello, mip: fx}]}
+      sea_ice_area_sh:
+        <<: *var_basic
+        preprocessor: pp_sh
+        reference_dataset: OSI-450-sh
+        additional_datasets:
+        - {dataset: OSI-450-sh, project: OBS, type: reanaly, version: v3, mip: OImon, tier: 2,
+           supplementary_variables: [{short_name: areacello, mip: fx}]}
+    scripts:
+      allplots:
+        script: monitor/multi_datasets.py
+        plot_folder: '{plot_dir}'
+        group_variables_by: variable_group
+        facet_used_for_labels: dataset
+        plots:
+          annual_cycle: &plot_basic
+            annual_mean_kwargs: False
+            pyplot_kwargs:
+              title: 'Sea ice area (million km2)'
+            plot_kwargs:
+              OSI-450-nh:
+                color: blue
+                label: 'OSISAF/CCI'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+              OSI-450-sh:
+                color: blue
+                label: 'OSISAF/CCI'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+              default:
+                color: red
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 4
+
+  siarea_min:
+    description: Create time series of September/February sea ice area.
+    variables:
+      sea_ice_area_nh_sep:
+        <<: *var_basic
+        timerange: '1979/2014'
+        preprocessor: pp_nh_sep
+        reference_dataset: OSI-450-nh
+        additional_datasets:
+        - {dataset: OSI-450-nh, project: OBS, type: reanaly, version: v3, mip: OImon, tier: 2,
+           supplementary_variables: [{short_name: areacello, mip: fx}], timerange: 1979/2020}
+      sea_ice_area_sh_feb:
+        <<: *var_basic
+        timerange: '1979/2014'
+        preprocessor: pp_sh_feb
+        reference_dataset: OSI-450-sh
+        additional_datasets:
+        - {dataset: OSI-450-sh, project: OBS, type: reanaly, version: v3, mip: OImon, tier: 2,
+           supplementary_variables: [{short_name: areacello, mip: fx}], timerange: 1979/2020}
+    scripts:
+      allplots:
+        script: monitor/multi_datasets.py
+        plot_folder: '{plot_dir}'
+        group_variables_by: variable_group
+        facet_used_for_labels: dataset
+        plots:
+          timeseries:
+            <<: *plot_basic

doc/sphinx/source/recipes/recipe_ref.rst

@@ -1,10 +1,10 @@
 @misc{osi-450,
-  doi = {10.15770/EUM_SAF_OSI_0008},
-  url = {https://navigator.eumetsat.int/product/EO:EUM:DAT:MULT:OSI-450},
+  doi = {10.15770/EUM_SAF_OSI_0013},
+  url = {https://user.eumetsat.int/catalogue/EO:EUM:DAT:0826},
   author = {{OSI SAF}},
   keywords = {Climate,  Thematic Climate Data Record,  Ocean,  Sea Ice},
   language = {en},
-  title = {Global Sea Ice Concentration Climate Data Record v2.0 - Multimission},
+  title = {Global Sea Ice Concentration Climate Data Record v3.0 - Multimission},
   publisher = {OSI SAF},
-  year = {2017}
+  year = {2022}
 }