Diagnostic for calculating Lamb Weathertypes

doc/sphinx/source/recipes/index.rst

@@ -64,6 +64,7 @@ Atmosphere
    recipe_validation
    recipe_radiation_budget
    recipe_aod_aeronet_assess
+   recipe_weathertyping
 
 Climate metrics
 ^^^^^^^^^^^^^^^

doc/sphinx/source/recipes/recipe_weathertyping.rst

@@ -0,0 +1,113 @@
+.. _recipes_weathertyping:
+
+Lamb Weathertypes
+=====
+
+Overview
+--------
+
+A diagnostic to calculate Lamb weathertypes over a given region. Furthermore, 
+correlations between weathertypes and precipitation patterns over a given area can be calculated 
+and 'combined' or 'simplified' weathertypes can be derived. Additionally, mean fields, as well as
+anomalies and standard deviations can be plotted.
+
+
+Available recipes and diagnostics
+---------------------------------
+
+Recipes are stored in esmvaltool/recipes/
+
+    * recipe_weathertyping.yml
+
+Diagnostics are stored in esmvaltool/diag_scripts/weathertyping/
+
+    * weathertyping.py: calculate lamb and simplified WT, plot mean, anomalies and std for each WT for psl, tas, pr
+
+
+User settings in recipe
+-----------------------
+
+#. weathertyping.py
+
+   *Required settings for script*
+
+   *Optional settings for script*
+
+   * correlation_threshold: correlation threshold for selecting similar WT pairs, only needed if automatic_slwt==True and predefined_slwt==False. default=0.9
+   * rmse_threshold: rmse threshold for selecting similar WT pairs, only needed if automatic_slwt==True and predefined_slwt==False. default=0.002
+   * plotting: if true, create plots of means, anomalies and std for psl, tas, prcp
+   * automatic_slwt: if true, automatically combine WT with similar precipitation patterns over specified area (via thresholds of correlation and rmse OR via predefined_slwt)
+   * predefined_slwt: dictionary of mappings between weathertypes
+
+.. note::
+
+  predefined_slwt can be a dictionary where keys are slwt and the values are arrays of lwt OR where keys are lwt and values are slwt
+
+   *Required settings for variables*
+
+   *Optional settings for variables*
+
+   *Required settings for preprocessor*
+
+   *Optional settings for preprocessor*
+
+   *Color tables*
+
+
+Variables
+---------
+
+* psl (atmos, day, time longitude latitude)
+* tas (atmos, day, time longitude latitude)
+* tp (atmos, day, time longitude latitude)
+* pr (atmos, day, time longitude latitude)
+
+
+Observations and reformat scripts
+---------------------------------
+
+*Note: (1) obs4MIPs data can be used directly without any preprocessing;
+(2) see headers of reformat scripts for non-obs4MIPs data for download
+instructions.*
+
+* E-OBS: European Climate Assessment & Dataset gridded daily precipitation sum
+* ERA5: ECMWF reanalysis
+
+References
+----------
+
+* Maraun, D., Truhetz, H., & Schaffer, A. (2021). Regional climate model biases, their dependence on synoptic circulation biases and the potential for bias adjustment: A process-oriented evaluation of the Austrian regional climate projections. Journal of Geophysical Research: Atmospheres, 126, e2020JD032824. https://doi.org/10.1029/2020JD032824 
+* Jones, P.D., Hulme, M. and Briffa, K.R. (1993), A comparison of Lamb circulation types with an objective classification scheme. Int. J. Climatol., 13: 655-663. https://doi.org/10.1002/joc.3370130606
+
+Example plots
+-------------
+
+.. _fig_weathertyping_1:
+.. figure::  /recipes/figures/weathertyping/lwt_1_ERA5_psl_mean.png
+   :align:   center
+
+   PSL mean map of Lamb WT 1 for ERA5.
+
+.. _fig_weathertyping_2:
+.. figure::  /recipes/figures/weathertyping/lwt_1_TaiESM1_psl_mean.png
+   :align:   center
+
+   PSL mean map of Lamb WT 1 for TaiESM1.
+
+.. _fig_weathertyping_3:
+.. figure::  /recipes/figures/weathertyping/slwt_EOBS_4_ERA5_psl_mean.png
+   :align:   center
+
+   PSL mean map of slwt_EOBS 4 for ERA5 (in this case combined Lamb WT 24 and 23).
+
+.. _fig_weathertyping_4:
+.. figure::  /recipes/figures/weathertyping/correlation_matrix_E-OBS.png
+   :align:   center
+
+   Heatmap of correlation values for Lamb WTs 1-27.
+
+.. _fig_weathertyping_5:
+.. figure::  /recipes/figures/weathertyping/ERA5__lwt_rel_occurrence.png
+   :align:   center
+
+   Stackplot of seasonal relative occureences of each WT for ERA5.

esmvaltool/config-references.yml

@@ -286,6 +286,11 @@ authors:
     name: Juckes, Martin
     institute: BADC, UK
     orcid:
+  jury_martin:
+    name: Jury, Martin
+    institute: WEGC, Austria
+    e-mail: [email protected]
+    orcid: https://orcid.org/0000-0003-0590-7843
   kadygrov_nikolay:
     name: Kadygrov, Nikolay
     institute: IPSL, France
@@ -315,6 +320,12 @@ authors:
     name: Krasting, John
     institute: NOAA, USA
     orcid: https://orcid.org/0000-0002-4650-9844
+  kroissenbrunner_thomas:
+    name: Kroissenbrunner, Thomas
+    institute: WEGC, Austria
+    email: [email protected]
+    github: thomaskroi1996
+    orcid: 
   kuehbacher_birgit:
     name: Kuehbacher, Birgit
     institute: DLR, Germany
@@ -793,6 +804,7 @@ projects:
   trr181: DFG Project TRR-181, Energy transfers in Atmosphere and Ocean
   ukesm: UKESM, UK Earth System Model project (NERC)
   usmile: ERC Synergy Grant USMILE
+  preval: PREVAL ÖKS
 
 
 realms:

esmvaltool/diag_scripts/climate_metrics/feedback_parameters.py

@@ -365,7 +365,7 @@ def _create_regression_plot(tas_cube,
     y_reg = reg.slope * x_reg + reg.intercept
 
     # Plot data
-    title = (f'{FEEDBACK_PARAMETERS.get(var,var)} TOA radiance for '
+    title = (f'{FEEDBACK_PARAMETERS.get(var, var)} TOA radiance for '
              f'{dataset_name}')
     filename = f'{var}_regression_{dataset_name}'
     if description is not None: