From 33bf72c1c738c69bed65ae40d54340f31ff73708 Mon Sep 17 00:00:00 2001 From: Naomi Parsons Date: Wed, 13 May 2026 14:29:37 +0100 Subject: [PATCH 1/2] Adding CMIP7 versions Changed end date, variables and test dataset --- .../recipe_seaice_sensitivity_cmip7.yml | 172 ++++++++++++++++++ 1 file changed, 172 insertions(+) create mode 100644 esmvaltool/recipes/recipe_seaice_sensitivity_cmip7.yml diff --git a/esmvaltool/recipes/recipe_seaice_sensitivity_cmip7.yml b/esmvaltool/recipes/recipe_seaice_sensitivity_cmip7.yml new file mode 100644 index 0000000000..68b8dd022a --- /dev/null +++ b/esmvaltool/recipes/recipe_seaice_sensitivity_cmip7.yml @@ -0,0 +1,172 @@ +# ESMValTool +--- +documentation: + title: Sea ice sensitivity + + description: | + Recipe for quantifying the sensitivity of sea ice to global warming. + Siconc data is summed for each hemisphere and then compared to the + change in globally meaned, annually meaned surface air temperature. + In the northern hemisphere, September sea ice data is used. + In the southern hemisphere, annual mean sea ice data is used. + Two plots are produced for each hemisphere, one showing the gradient + of the direct regression of sea ice area over temperature, and the + other showing the two separate trends over time. + + authors: + - parsons_naomi + - sellar_alistair + - blockley_ed + + maintainer: + - parsons_naomi + +defaults: &defaults {activity: CMIP, project: CMIP7, exp: historical, } +datasets: + - {<<: *defaults, dataset: UKCM2-0-LL, institute: MOHC, ensemble: r1i1p1f3, exp: 1pctCO2, grid: g100, label_dataset: True} + + +preprocessors: + extract_test_period: + &extract_test_period + extract_time: + start_day: 1 + start_month: 1 + start_year: 1979 + end_day: 31 + end_month: 12 + end_year: 2021 + + extract_sept: + &extract_sept + extract_month: + month: 9 + + nh_total_area: + &nh_total_area + extract_region: + start_longitude: 0 + end_longitude: 360 + start_latitude: 0 + end_latitude: 90 + area_statistics: + operator: sum + convert_units: + units: 1e6 km2 + + sh_total_area: + &sh_total_area + extract_region: + start_longitude: 0 + end_longitude: 360 + start_latitude: -90 + end_latitude: 0 + area_statistics: + operator: sum + convert_units: + units: 1e6 km2 + + global_mean: + &global_mean + area_statistics: + operator: mean + + annual_mean: + &annual_mean + annual_statistics: + operator: mean + keep_group_coordinates: true + + pp_arctic_sept_sea_ice: + <<: *extract_test_period + <<: *extract_sept + <<: *nh_total_area + + pp_antarctic_avg_ann_sea_ice: + <<: *extract_test_period + <<: *annual_mean + <<: *sh_total_area + + pp_avg_ann_global_temp: + <<: *extract_test_period + <<: *global_mean + <<: *annual_mean + +diagnostics: + arctic: + description: Plots September sea ice sensitivity above 0 latitude in millions of square kilometres + variables: + siconc: &siconc + preprocessor: pp_arctic_sept_sea_ice + realm: &siconc_realm seaIce + mip: *siconc_realm + frequency: mon + region: glb + branding_suffix: tavg-u-hxy-u + tas: &tas + preprocessor: pp_avg_ann_global_temp + realm: &tas_realm atmos + mip: *tas_realm + frequency: mon + region: glb + branding_suffix: tavg-h2m-hxy-u + scripts: + sea_ice_sensitivity_script: + script: seaice/seaice_sensitivity.py + observations: + observation period: 1979-2014 + sea ice sensitivity (Notz-style plot): + mean: -4.01 + standard deviation: 0.32 + plausible range: 1.28 + annual trends (Roach-style plot): + first point: + GMST trend: + SIA trend: + Pearson CC of SIA over GMST: + significance of SIA over GMST: + second point: + GMST trend: + SIA trend: + Pearson CC of SIA over GMST: + significance of SIA over GMST: + third point: + GMST trend: + SIA trend: + Pearson CC of SIA over GMST: + significance of SIA over GMST: + + + antarctic: + description: Plots annual mean sea ice sensitivity below 0 latitude in millions of square kilometres + variables: + siconc: + <<: *siconc + preprocessor: pp_antarctic_avg_ann_sea_ice + tas: + *tas + scripts: + sea_ice_sensitivity_script: + script: seaice/seaice_sensitivity.py + observations: + observation period: + sea ice sensitivity (Notz-style plot): + mean: + standard deviation: + plausible range: + annual trends (Roach-style plot): + first point: + GMST trend: + SIA trend: + Pearson CC of SIA over GMST: + significance of SIA over GMST: + second point: + GMST trend: + SIA trend: + Pearson CC of SIA over GMST: + significance of SIA over GMST: + third point: + GMST trend: + SIA trend: + Pearson CC of SIA over GMST: + significance of SIA over GMST: From 2399cfc11a20dad750d05736ecee4ac5c4351ee8 Mon Sep 17 00:00:00 2001 From: Naomi Parsons Date: Tue, 19 May 2026 16:46:07 +0100 Subject: [PATCH 2/2] Could roll back: Trying to work with observations. Dates are tricky as fake CMIP7 data is only 1900-1910 and observations are later... --- .../recipe_seaice_sensitivity_cmip7.yml | 102 +++++++++--------- 1 file changed, 53 insertions(+), 49 deletions(-) diff --git a/esmvaltool/recipes/recipe_seaice_sensitivity_cmip7.yml b/esmvaltool/recipes/recipe_seaice_sensitivity_cmip7.yml index 68b8dd022a..9a0b594530 100644 --- a/esmvaltool/recipes/recipe_seaice_sensitivity_cmip7.yml +++ b/esmvaltool/recipes/recipe_seaice_sensitivity_cmip7.yml @@ -21,9 +21,22 @@ documentation: maintainer: - parsons_naomi -defaults: &defaults {activity: CMIP, project: CMIP7, exp: historical, } -datasets: - - {<<: *defaults, dataset: UKCM2-0-LL, institute: MOHC, ensemble: r1i1p1f3, exp: 1pctCO2, grid: g100, label_dataset: True} +model_defaults: &model_defaults {activity: CMIP, project: CMIP7, exp: historical, } +model_datasets: &model_datasets + - {<<: *model_defaults, dataset: UKCM2-0-LL, institute: MOHC, ensemble: r1i1p1f3, exp: 1pctCO2, grid: g100, label_dataset: True} + +obs_defaults: &obs_defaults { project: OBS, observation: True, tier: 2} +tasa_obs: &tasa_obs +# - { <<: *obs_defaults, dataset: HadCRUT5, type: ground, version: 5.0.1.0-analysis } + - { <<: *obs_defaults, dataset: NOAAGlobalTemp, type: ground, version: v5.0.0, mip: Amon } + +arctic_siconc_obs: &arctic_siconc_obs + - { <<: *obs_defaults, dataset: OSI-450-nh, type: reanaly, version: v3, mip: OImon, supplementary_variables: [{short_name: areacello, mip: fx}]} + - { <<: *obs_defaults, dataset: HadISST, type: reanaly, version: 1, mip: OImon} + +antarctic_siconc_obs: &antarctic_siconc_obs + - { <<: *obs_defaults, dataset: OSI-450-sh, type: reanaly, version: v3, mip: OImon, supplementary_variables: [{short_name: areacello, mip: fx}]} + - { <<: *obs_defaults, dataset: HadISST, type: reanaly, version: 1, mip: OImon} preprocessors: @@ -32,10 +45,10 @@ preprocessors: extract_time: start_day: 1 start_month: 1 - start_year: 1979 + start_year: &data_start 1900 # TODO: 1979 end_day: 31 end_month: 12 - end_year: 2021 + end_year: &data_end 1909 # TODO: 2021 extract_sept: &extract_sept @@ -75,7 +88,6 @@ preprocessors: &annual_mean annual_statistics: operator: mean - keep_group_coordinates: true pp_arctic_sept_sea_ice: <<: *extract_test_period @@ -103,6 +115,13 @@ diagnostics: frequency: mon region: glb branding_suffix: tavg-u-hxy-u + additional_datasets: + *model_datasets + siconc_obs: + <<: *siconc + short_name: siconc + additional_datasets: + *arctic_siconc_obs tas: &tas preprocessor: pp_avg_ann_global_temp realm: &tas_realm atmos @@ -110,32 +129,27 @@ diagnostics: frequency: mon region: glb branding_suffix: tavg-h2m-hxy-u + additional_datasets: + *model_datasets + tasa_obs: &tasa + <<: *tas + short_name: tasa + additional_datasets: + *tasa_obs scripts: sea_ice_sensitivity_script: script: seaice/seaice_sensitivity.py observations: - observation period: 1979-2014 - sea ice sensitivity (Notz-style plot): + data_period: &data_period + start_year: *data_start + end_year: *data_end + observation_period: + start_year: 1979 + end_year: 2014 + sea_ice_sensitivity: mean: -4.01 - standard deviation: 0.32 - plausible range: 1.28 - annual trends (Roach-style plot): - first point: - GMST trend: - SIA trend: - Pearson CC of SIA over GMST: - significance of SIA over GMST: - second point: - GMST trend: - SIA trend: - Pearson CC of SIA over GMST: - significance of SIA over GMST: - third point: - GMST trend: - SIA trend: - Pearson CC of SIA over GMST: - significance of SIA over GMST: - + standard_deviation: 0.32 + plausible_range: 1.28 antarctic: description: Plots annual mean sea ice sensitivity below 0 latitude in millions of square kilometres @@ -143,30 +157,20 @@ diagnostics: siconc: <<: *siconc preprocessor: pp_antarctic_avg_ann_sea_ice + additional_datasets: + *model_datasets + siconc_obs: + <<: *siconc + short_name: siconc + preprocessor: pp_antarctic_avg_ann_sea_ice + additional_datasets: + *antarctic_siconc_obs tas: - *tas + <<: *tas + tasa_obs: + <<: *tasa scripts: sea_ice_sensitivity_script: script: seaice/seaice_sensitivity.py observations: - observation period: - sea ice sensitivity (Notz-style plot): - mean: - standard deviation: - plausible range: - annual trends (Roach-style plot): - first point: - GMST trend: - SIA trend: - Pearson CC of SIA over GMST: - significance of SIA over GMST: - second point: - GMST trend: - SIA trend: - Pearson CC of SIA over GMST: - significance of SIA over GMST: - third point: - GMST trend: - SIA trend: - Pearson CC of SIA over GMST: - significance of SIA over GMST: + data_period: *data_period