diff --git a/esmvaltool/recipes/ref/recipe_ref_annual_cycle_region_cmip7.yml b/esmvaltool/recipes/ref/recipe_ref_annual_cycle_region_cmip7.yml
new file mode 100644
index 0000000000..64118d65a8
--- /dev/null
+++ b/esmvaltool/recipes/ref/recipe_ref_annual_cycle_region_cmip7.yml
@@ -0,0 +1,227 @@
+# ESMValTool
+---
+documentation:
+  title: Benchmarking of a single model for a region through extract_shape.
+  description: >
+    Benchmarking: annual cycle.
+  authors:
+    - lauer_axel
+    - bock_lisa
+    - hassler_birgit
+    - lindenlaub_lukas
+    - schlund_manuel
+    - weigel_katja
+    - andela_bouwe
+  maintainer:
+    - lauer_axel
+  references:
+    - lauer25gmd
+  projects:
+    - dlrmabak
+
+regions: &regions
+  - "Arabian-Peninsula"
+  - "Arabian-Sea"
+  - "Arctic-Ocean"
+  - "Bay-of-Bengal"
+  - "C.Australia"
+  - "C.North-America"
+  - "Caribbean"
+  - "Central-Africa"
+  - "E.Antarctica"
+  - "E.Asia"
+  - "E.Australia"
+  - "E.C.Asia"
+  - "E.Europe"
+  - "E.North-America"
+  - "E.Siberia"
+  - "E.Southern-Africa"
+  - "Equatorial.Atlantic-Ocean"
+  - "Equatorial.Indic-Ocean"
+  - "Equatorial.Pacific-Ocean"
+  - "Greenland/Iceland"
+  - "Madagascar"
+  - "Mediterranean"
+  - "N.Atlantic-Ocean"
+  - "N.Australia"
+  - "N.Central-America"
+  - "N.E.North-America"
+  - "N.E.South-America"
+  - "N.Eastern-Africa"
+  - "N.Europe"
+  - "N.Pacific-Ocean"
+  - "N.South-America"
+  - "N.W.North-America"
+  - "N.W.South-America"
+  - "New-Zealand"
+  - "Russian-Arctic"
+  - "Russian-Far-East"
+  - "S.Asia"
+  - "S.Atlantic-Ocean"
+  - "S.Australia"
+  - "S.Central-America"
+  - "S.E.Asia"
+  - "S.E.South-America"
+  - "S.Eastern-Africa"
+  - "S.Indic-Ocean"
+  - "S.Pacific-Ocean"
+  - "S.South-America"
+  - "S.W.South-America"
+  - "Sahara"
+  - "South-American-Monsoon"
+  - "Southern-Ocean"
+  - "Tibetan-Plateau"
+  - "W.Antarctica"
+  - "W.C.Asia"
+  - "W.North-America"
+  - "W.Siberia"
+  - "W.Southern-Africa"
+  - "West&Central-Europe"
+  - "Western-Africa"
+
+datasets:
+  # Dataset to be benchmarked
+  - dataset: DUMMY-MODEL
+    drs_specs: MIP-DRS7
+    project: CMIP7
+    exp: historical
+    ensemble: r1i1p1f2
+    grid: g999
+    activity: CMIP
+    institute: MOHC
+    mip: atmos
+    frequency: mon
+    region: glb
+
+preprocessors:
+  pp_abs:
+    regrid:
+      target_grid: 1x1
+      scheme: linear
+    extract_shape:
+      shapefile: ar6
+      decomposed: true
+      ids:
+        Name: *regions
+    climate_statistics:
+      period: monthly
+    area_statistics:
+      operator: mean
+
+  pp_abs_200:
+    extract_levels:
+      levels: 20000
+      scheme: linear
+      coordinate: air_pressure
+    regrid:
+      target_grid: 1x1
+      scheme: linear
+    extract_shape:
+      shapefile: ar6
+      decomposed: true
+      ids:
+        Name: *regions
+    climate_statistics:
+      period: monthly
+    area_statistics:
+      operator: mean
+
+
+diagnostics:
+  near-surface_air_temperature_annual_cycle:
+    description: Plot global mean mean annual cycle for the given period including a reference dataset.
+    variables:
+      tas: &var_default_abs
+        start_year: 1992
+        end_year: 2021
+        preprocessor: pp_abs
+        branding_suffix: tavg-h2m-hxy-u
+        mip: Amon
+        additional_datasets:
+          - {dataset: HadCRUT5, project: OBS, type: ground, version: 5.0.1.0-analysis, tier: 2}
+    scripts:
+      figures: &figures
+        script: monitor/multi_datasets.py
+        plot_folder: '{plot_dir}'
+        group_variables_by: shape_id
+        plot_filename: "{plot_type}_{real_name}_{dataset}_{mip}_{exp}_{ensemble}_{shape_id}"
+        plots:
+          annual_cycle:
+            plot_kwargs:
+              default:
+                color: red
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 4
+              HadCRUT5:
+                color: black
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+              ERA5:
+                color: black
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+              ERA-5:
+                color: black
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+              GPCP-V3.3:
+                color: black
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+            pyplot_kwargs:
+              title: '{long_name} {shape_id}'
+              ylabel: '{short_name} [{units}]'
+  precipitation_annual_cycle:
+    description: Plot global mean mean annual cycle for the given period including a reference dataset.
+    variables:
+      pr:
+        <<: *var_default_abs
+        branding_suffix: tavg-u-hxy-u
+        additional_datasets:
+          - {dataset: GPCP-V3.3, project: obs4MIPs, tier: 1}
+    scripts:
+      figures: *figures
+  sea_level_pressure_annual_cycle:
+    description: Plot global mean mean annual cycle for the given period including a reference dataset.
+    variables:
+      psl:
+        <<: *var_default_abs
+        branding_suffix: tavg-u-hxy-u
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1}
+    scripts:
+      figures: *figures
+  specific_humidity_annual_cycle:
+    description: Plot global mean mean annual cycle for the given period including a reference dataset.
+    variables:
+      hus_200:
+        <<: *var_default_abs
+        short_name: hus
+        preprocessor: pp_abs_200
+        branding_suffix: tavg-p19-hxy-u
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1}
+    scripts:
+      figures: *figures
+  eastward_wind_annual_cycle:
+    description: Plot global mean mean annual cycle for the given period including a reference dataset.
+    variables:
+      ua_200:
+        <<: *var_default_abs
+        short_name: ua
+        preprocessor: pp_abs_200
+        branding_suffix: tavg-p19-hxy-air
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1}
+    scripts:
+      figures: *figures
diff --git a/esmvaltool/recipes/ref/recipe_ref_timeseries_region_cmip7.yml b/esmvaltool/recipes/ref/recipe_ref_timeseries_region_cmip7.yml
new file mode 100644
index 0000000000..1b4f2bfded
--- /dev/null
+++ b/esmvaltool/recipes/ref/recipe_ref_timeseries_region_cmip7.yml
@@ -0,0 +1,338 @@
+# ESMValTool
+---
+documentation:
+  title: Benchmarking of a single model for a region through extract_shape.
+  description: >
+    Benchmarking: time series.
+  authors:
+    - lauer_axel
+    - bock_lisa
+    - hassler_birgit
+    - lindenlaub_lukas
+    - schlund_manuel
+    - weigel_katja
+    - andela_bouwe
+  maintainer:
+    - lauer_axel
+  references:
+    - lauer25gmd
+  projects:
+    - dlrmabak
+
+regions: &regions
+  - "Arabian-Peninsula"
+  - "Arabian-Sea"
+  - "Arctic-Ocean"
+  - "Bay-of-Bengal"
+  - "C.Australia"
+  - "C.North-America"
+  - "Caribbean"
+  - "Central-Africa"
+  - "E.Antarctica"
+  - "E.Asia"
+  - "E.Australia"
+  - "E.C.Asia"
+  - "E.Europe"
+  - "E.North-America"
+  - "E.Siberia"
+  - "E.Southern-Africa"
+  - "Equatorial.Atlantic-Ocean"
+  - "Equatorial.Indic-Ocean"
+  - "Equatorial.Pacific-Ocean"
+  - "Greenland/Iceland"
+  - "Madagascar"
+  - "Mediterranean"
+  - "N.Atlantic-Ocean"
+  - "N.Australia"
+  - "N.Central-America"
+  - "N.E.North-America"
+  - "N.E.South-America"
+  - "N.Eastern-Africa"
+  - "N.Europe"
+  - "N.Pacific-Ocean"
+  - "N.South-America"
+  - "N.W.North-America"
+  - "N.W.South-America"
+  - "New-Zealand"
+  - "Russian-Arctic"
+  - "Russian-Far-East"
+  - "S.Asia"
+  - "S.Atlantic-Ocean"
+  - "S.Australia"
+  - "S.Central-America"
+  - "S.E.Asia"
+  - "S.E.South-America"
+  - "S.Eastern-Africa"
+  - "S.Indic-Ocean"
+  - "S.Pacific-Ocean"
+  - "S.South-America"
+  - "S.W.South-America"
+  - "Sahara"
+  - "South-American-Monsoon"
+  - "Southern-Ocean"
+  - "Tibetan-Plateau"
+  - "W.Antarctica"
+  - "W.C.Asia"
+  - "W.North-America"
+  - "W.Siberia"
+  - "W.Southern-Africa"
+  - "West&Central-Europe"
+  - "Western-Africa"
+
+datasets:
+  # Dataset to be benchmarked
+  - dataset: DUMMY-MODEL
+    drs_specs: MIP-DRS7
+    project: CMIP7
+    exp: historical
+    ensemble: r1i1p1f2
+    grid: g999
+    activity: CMIP
+    institute: MOHC
+    mip: atmos
+    frequency: mon
+    region: glb
+
+
+preprocessors:
+  pp_abs:
+    regrid:
+      target_grid: 1x1
+      scheme: linear
+    extract_shape:
+      shapefile: ar6
+      crop: true
+      decomposed: true
+      ids:
+        Name: *regions
+    area_statistics:
+      operator: mean
+
+  pp_abs_200:
+    extract_levels:
+      levels: 20000
+      scheme: linear
+      coordinate: air_pressure
+    regrid:
+      target_grid: 1x1
+      scheme: linear
+    extract_shape:
+      shapefile: ar6
+      crop: true
+      decomposed: true
+      ids:
+        Name: *regions
+    area_statistics:
+      operator: mean
+
+  pp:
+    regrid:
+      target_grid: 1x1
+      scheme: linear
+    extract_shape:
+      shapefile: ar6
+      decomposed: true
+      ids:
+        Name: *regions
+    anomalies:
+      period: month
+      reference:
+        start_year: 1992
+        start_month: 1
+        start_day: 1
+        end_year: 2021
+        end_month: 12
+        end_day: 31
+    area_statistics:
+      operator: mean
+
+  pp_200:
+    extract_levels:
+      levels: 20000
+      scheme: linear
+      coordinate: air_pressure
+    regrid:
+      target_grid: 1x1
+      scheme: linear
+    extract_shape:
+      shapefile: ar6
+      decomposed: true
+      ids:
+        Name: *regions
+    anomalies:
+      period: month
+      reference:
+        start_year: 1992
+        start_month: 1
+        start_day: 1
+        end_year: 2021
+        end_month: 12
+        end_day: 31
+    area_statistics:
+      operator: mean
+
+
+diagnostics:
+
+  near-surface_air_temperature_anomalies:
+    description: Plot "classical" time series of global mean anomalies including a reference dataset.
+    variables:
+      tas: &var_default
+        start_year: 1992
+        end_year: 2021
+        preprocessor: pp
+        mip: Amon
+        branding_suffix: tavg-h2m-hxy-u
+        additional_datasets:
+          - {dataset: HadCRUT5, project: OBS, type: ground, version: 5.0.1.0-analysis, tier: 2}
+    scripts:
+      figures: &figures
+        script: monitor/multi_datasets.py
+        plot_folder: '{plot_dir}'
+        group_variables_by: shape_id
+        plot_filename: "{plot_type}_{real_name}_{dataset}_{mip}_{exp}_{ensemble}_{shape_id}"
+        plots:
+          timeseries:
+            plot_kwargs: &plot_kwargs
+              default:
+                color: red
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 4
+              HadCRUT5:
+                color: black
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+              ERA5:
+                color: black
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+              ERA-5:
+                color: black
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+              GPCP-V3.3:
+                color: black
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+            pyplot_kwargs:
+              title: '{long_name} Anomalies {shape_id}'
+              ylabel: '{short_name} [{units}]'
+  precipitation_anomalies:
+    description: Plot "classical" time series of global mean anomalies including a reference dataset.
+    variables:
+      pr:
+        <<: *var_default
+        branding_suffix: tavg-u-hxy-u
+        additional_datasets:
+          - {dataset: GPCP-V3.3, project: obs4MIPs, tier: 1}
+    scripts:
+      figures: *figures
+  sea_level_pressure_anomalies:
+    description: Plot "classical" time series of global mean anomalies including a reference dataset.
+    variables:
+      psl:
+        <<: *var_default
+        branding_suffix: tavg-u-hxy-u
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1}
+    scripts:
+      figures: *figures
+  specific_humidity_anomalies:
+    description: Plot "classical" time series of global mean anomalies including a reference dataset.
+    variables:
+      hus_200:
+        <<: *var_default
+        short_name: hus
+        preprocessor: pp_200
+        branding_suffix: tavg-p19-hxy-u
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1}
+    scripts:
+      figures: *figures
+  eastward_wind_anomalies:
+    description: Plot "classical" time series of global mean anomalies including a reference dataset.
+    variables:
+      ua_200:
+        <<: *var_default
+        short_name: ua
+        preprocessor: pp_200
+        branding_suffix: tavg-p19-hxy-air
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1}
+    scripts:
+      figures: *figures
+
+  near-surface_air_temperature:
+    description: Plot "classical" time series of global mean including a reference dataset.
+    variables:
+      tas: &var_default_abs
+        <<: *var_default
+        preprocessor: pp_abs
+        branding_suffix: tavg-h2m-hxy-u
+        additional_datasets:
+          - {dataset: HadCRUT5, project: OBS, type: ground, version: 5.0.1.0-analysis, tier: 2}
+    scripts:
+      figures: &figures_abs
+        script: monitor/multi_datasets.py
+        plot_folder: '{plot_dir}'
+        group_variables_by: shape_id
+        plot_filename: "{plot_type}_{real_name}_{dataset}_{mip}_{exp}_{ensemble}_{shape_id}"
+        plots:
+          timeseries:
+            plot_kwargs: *plot_kwargs
+            pyplot_kwargs:
+              title: '{long_name} {shape_id}'
+  precipitation:
+    description: Plot "classical" time series of global mean including a reference dataset.
+    variables:
+      pr:
+        <<: *var_default_abs
+        branding_suffix: tavg-u-hxy-u
+        additional_datasets:
+          - {dataset: GPCP-V3.3, project: obs4MIPs, tier: 1}
+    scripts:
+      figures: *figures_abs
+  sea_level_pressure:
+    description: Plot "classical" time series of global mean including a reference dataset.
+    variables:
+      psl:
+        <<: *var_default_abs
+        branding_suffix: tavg-u-hxy-u
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1}
+    scripts:
+      figures: *figures_abs
+  specific_humidity:
+    description: Plot "classical" time series of global mean including a reference dataset.
+    variables:
+      hus_200:
+        <<: *var_default_abs
+        short_name: hus
+        preprocessor: pp_abs_200
+        branding_suffix: tavg-p19-hxy-u
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1}
+    scripts:
+      figures: *figures_abs
+  eastward_wind:
+    description: Plot "classical" time series of global mean including a reference dataset.
+    variables:
+      ua_200:
+        <<: *var_default_abs
+        short_name: ua
+        preprocessor: pp_abs_200
+        branding_suffix: tavg-p19-hxy-air
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1}
+    scripts:
+      figures: *figures_abs
diff --git a/esmvaltool/recipes/ref/recipe_ref_trend_regions_cmip7.yml b/esmvaltool/recipes/ref/recipe_ref_trend_regions_cmip7.yml
new file mode 100644
index 0000000000..4317e54c38
--- /dev/null
+++ b/esmvaltool/recipes/ref/recipe_ref_trend_regions_cmip7.yml
@@ -0,0 +1,334 @@
+# ESMValTool
+---
+documentation:
+  description:
+    Box plot, linear trends for
+    the 46 IPCC AR6 reference regions for land (Iturbide et al., 2020).
+    The diagnostic called in this recipe uses the seaborn.boxplot function.
+    See seaborn manual for detailed information
+    (https://seaborn.pydata.org/generated/seaborn.boxplot.html).
+
+  authors:
+    - schlund_manuel
+    - weigel_katja
+  maintainer:
+    - weigel_katja
+  title: Linear trends.
+
+datasets:
+  # - {dataset: MIROC6,      project: CMIP6, exp: historical, mip: Amon, ensemble: r1i1p1f1, grid: gn}
+  - dataset: DUMMY-MODEL
+    drs_specs: MIP-DRS7
+    project: CMIP7
+    exp: historical
+    ensemble: r1i1p1f2
+    grid: g999
+    activity: CMIP
+    institute: MOHC
+    mip: atmos
+    frequency: mon
+    region: glb
+
+preprocessors:
+  trend_tas: &trend
+    custom_order: true
+    extract_shape:
+      shapefile: ar6
+      decomposed: true
+      ids:
+        Acronym:
+          - 'NWN'
+          - 'NEN'
+          - 'GIC'
+          - 'WNA'
+          - 'CNA'
+          - 'ENA'
+          - 'NCA'
+          - 'SCA'
+          - 'CAR'
+          - 'NWS'
+          - 'NSA'
+          - 'NES'
+          - 'SAM'
+          - 'SWS'
+          - 'SES'
+          - 'SSA'
+          - 'NEU'
+          - 'WCE'
+          - 'EEU'
+          - 'MED'
+          - 'SAH'
+          - 'WAF'
+          - 'CAF'
+          - 'NEAF'
+          - 'SEAF'
+          - 'WSAF'
+          - 'ESAF'
+          - 'MDG'
+          - 'RAR'
+          - 'WSB'
+          - 'ESB'
+          - 'RFE'
+          - 'WCA'
+          - 'ECA'
+          - 'TIB'
+          - 'EAS'
+          - 'ARP'
+          - 'SAS'
+          - 'SEA'
+          - 'NAU'
+          - 'CAU'
+          - 'EAU'
+          - 'SAU'
+          - 'NZ'
+          - 'WAN'
+          - 'EAN'
+    area_statistics:
+      operator: mean
+    linear_trend:
+    convert_units:
+      units: K yr-1
+
+  trend_pr:
+    <<: *trend
+    convert_units:
+      units: mm kg m-3 day-1 yr-1
+
+  trend_psl:
+    <<: *trend
+    convert_units:
+      units: m-1 kg s-2 yr-1
+
+
+  trend_ua200: &trend200
+    custom_order: true
+    extract_levels:
+      levels: 20000
+      scheme: linear
+      coordinate: air_pressure
+    extract_shape:
+      shapefile: ar6
+      decomposed: true
+      ids:
+        Acronym:
+          - 'NWN'
+          - 'NEN'
+          - 'GIC'
+          - 'WNA'
+          - 'CNA'
+          - 'ENA'
+          - 'NCA'
+          - 'SCA'
+          - 'CAR'
+          - 'NWS'
+          - 'NSA'
+          - 'NES'
+          - 'SAM'
+          - 'SWS'
+          - 'SES'
+          - 'SSA'
+          - 'NEU'
+          - 'WCE'
+          - 'EEU'
+          - 'MED'
+          - 'SAH'
+          - 'WAF'
+          - 'CAF'
+          - 'NEAF'
+          - 'SEAF'
+          - 'WSAF'
+          - 'ESAF'
+          - 'MDG'
+          - 'RAR'
+          - 'WSB'
+          - 'ESB'
+          - 'RFE'
+          - 'WCA'
+          - 'ECA'
+          - 'TIB'
+          - 'EAS'
+          - 'ARP'
+          - 'SAS'
+          - 'SEA'
+          - 'NAU'
+          - 'CAU'
+          - 'EAU'
+          - 'SAU'
+          - 'NZ'
+          - 'WAN'
+          - 'EAN'
+    area_statistics:
+      operator: mean
+    linear_trend:
+    convert_units:
+      units: m s-1 yr-1
+
+  trend_hus200:
+    <<: *trend200
+    convert_units:
+      units: year-1
+
+
+diagnostics:
+  tas_trends:
+    variables:
+      tas:
+        preprocessor: trend_tas
+        mip: Amon
+        timerange: 1992/2021
+        branding_suffix: tavg-h2m-hxy-u
+        additional_datasets:
+          - {dataset: HadCRUT5, project: OBS, type: ground,
+             version: 5.0.1.0-analysis, tier: 2}
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1,
+             reference_for_metric: true}
+    scripts:
+      plot:
+        script: seaborn_diag.py
+        seaborn_func: barplot
+        reset_index: true
+        write_netcdf: true
+        seaborn_kwargs:
+          x: shape_id
+          y: tas
+          hue: dataset
+        add_aux_coords: true
+        facets_as_columns: [dataset]
+        legend_title: Data
+        plot_object_methods:
+          set_xlabel: region
+          set_ylabel: tas trend [K/yr]
+        seaborn_settings:
+          style: ticks
+          rc:
+            figure.figsize: [28, 6]
+        suptitle: Trend in Near-Surface Air Temperature
+
+  pr_trends:
+    variables:
+      pr:
+        preprocessor: trend_pr
+        mip: Amon
+        timerange: 1992/2021
+        branding_suffix: tavg-u-hxy-u
+        additional_datasets:
+          - {dataset: GPCP-V3.3, project: obs4MIPs, tier: 1,
+             reference_for_metric: true}
+          - {dataset: ERA5, project: native6, type: reanaly, version: v1,
+             tier: 3, reference_for_metric: true}
+    scripts:
+      plot:
+        script: seaborn_diag.py
+        seaborn_func: barplot
+        reset_index: true
+        write_netcdf: true
+        seaborn_kwargs:
+          x: shape_id
+          y: pr
+          hue: dataset
+        add_aux_coords: true
+        facets_as_columns: [dataset]
+        legend_title: Data
+        plot_object_methods:
+          set_xlabel: region
+          set_ylabel: pr trend [mm kg m-3 day-1 yr-1]
+        seaborn_settings:
+          style: ticks
+          rc:
+            figure.figsize: [28, 6]
+        suptitle: Trend in Precipitation
+
+
+  psl_trends:
+    variables:
+      psl:
+        preprocessor: trend_psl
+        mip: Amon
+        timerange: 1992/2021
+        branding_suffix: tavg-u-hxy-u
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1,
+             reference_for_metric: true}
+    scripts:
+      plot:
+        script: seaborn_diag.py
+        seaborn_func: barplot
+        reset_index: true
+        write_netcdf: true
+        seaborn_kwargs:
+          x: shape_id
+          y: psl
+          hue: dataset
+        add_aux_coords: true
+        facets_as_columns: [dataset]
+        legend_title: Data
+        plot_object_methods:
+          set_xlabel: region
+          set_ylabel: psl trend [m-1 kg s-2 yr-1]
+        seaborn_settings:
+          style: ticks
+          rc:
+            figure.figsize: [28, 6]
+        suptitle: Trend in Sea Level Pressure
+
+  ua200_trends:
+    variables:
+      ua:
+        preprocessor: trend_ua200
+        mip: Amon
+        timerange: 1992/2021
+        branding_suffix: tavg-p19-hxy-air
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1,
+             reference_for_metric: true}
+    scripts:
+      plot:
+        script: seaborn_diag.py
+        seaborn_func: barplot
+        reset_index: true
+        write_netcdf: true
+        seaborn_kwargs:
+          x: shape_id
+          y: ua
+          hue: dataset
+        add_aux_coords: true
+        facets_as_columns: [dataset]
+        plot_object_methods:
+          set_xlabel: region
+          set_ylabel: ua trend [m s-1 yr-1]
+        seaborn_settings:
+          style: ticks
+          rc:
+            figure.figsize: [28, 6]
+        suptitle: Trend in Zonal Wind at 200 hPa
+
+  hus200_trends:
+    variables:
+      hus:
+        preprocessor: trend_hus200
+        mip: Amon
+        timerange: 1992/2021
+        branding_suffix: tavg-p19-hxy-u
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1,
+             reference_for_metric: true}
+    scripts:
+      plot:
+        script: seaborn_diag.py
+        seaborn_func: barplot
+        reset_index: true
+        write_netcdf: true
+        seaborn_kwargs:
+          x: shape_id
+          y: hus
+          hue: dataset
+        add_aux_coords: true
+        facets_as_columns: [dataset]
+        plot_object_methods:
+          set_xlabel: region
+          set_ylabel: hus trend [yr-1]
+        seaborn_settings:
+          style: ticks
+          rc:
+            figure.figsize: [28, 6]
+        suptitle: Trend in Specific Humidity at 200 hPa