New EAT Lancet Diet Implementation

CHANGELOG.md

@@ -8,6 +8,8 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/).
 ## [Unreleased]
 
 ### changed
+- **15_food** Added improved EAT Lancet diet implementation (EAT 2p0)
+- **default.cfg** update input data to rev 4.99 (new validation data, new EATLancet recommendations, new semi-natural vegtation cropland data)
 - **start_functions** Check if cfg$recalibrate is consistent with cfg$gms$s14_use_yield_calib
 - **31_past** in grasslands_apr22 realization: changed structure of f31_pastr_suitability to align with ssp-rcp specific input data formulation. Changed input filename from cs3 to cs2. Added `cc`, `nocc` and `nocc_hist` options for `i31_manpast_suit` and changed input gams code from table to parameter. Climate scenario assignment moved from preloop.gms to input.gms. Removed pastSuit set in sets.gms as not needed anymore. Adjusted not_used.txt in both grasslands_apr22 and static realizations.
 - **default.cfg** Default for cfg$gms$peatland changed from `on` to `v2`
@@ -16,7 +18,7 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/).
 - **scenario_config.csv** settings for cfg$gms$s35_secdf_distribution in `ForestryEndo` and `ForestryExo` changed from `2` to `0`
 - **scenario_config.csv** removed erroneous setting `cc` from column `input`
 - **default.cfg** update default `cfg$gms$c56_pollutant_prices` and `cfg$gms$c60_2ndgen_biodem` to `R32M46-SSP2EU-NPi`
-- **default.cfg** update input data to rev 4.96 
+- **default.cfg** update input data to rev 4.96
 - **default.cfg** changed default for `cfg$gms$s56_buffer_aff` from 0.2 to 0.5
 - **default.cfg** changed default for `cfg$gms$s32_aff_prot` from 0 to 1
 - **21_trade** s21_trade_bal_damper for roundwood changed from 0.75 to 0.65

config/default.cfg

@@ -22,10 +22,10 @@ cfg$model <- "main.gms"   #def = "main.gms"
 #### input settings ####
 
 # which input data sets should be used?
-cfg$input <- c(regional    = "rev4.96_h12_magpie.tgz",
-               cellular    = "rev4.96_h12_fd712c0b_cellularmagpie_c200_MRI-ESM2-0-ssp370_lpjml-8e6c5eb1.tgz",
-               validation  = "rev4.96_h12_validation.tgz",
-               additional  = "additional_data_rev4.47.tgz",
+cfg$input <- c(regional    = "rev4.99_h12_magpie.tgz",
+               cellular    = "rev4.99_h12_fd712c0b_cellularmagpie_c200_MRI-ESM2-0-ssp370_lpjml-8e6c5eb1.tgz",
+               validation  = "rev4.99_h12_validation.tgz",
+               additional  = "additional_data_rev4.46.tgz",
                calibration = "calibration_H12_per_ton_fao_may22_glo_08Aug23.tgz")
 
 # NOTE: It is recommended to recalibrate the model when changing cellular input data
@@ -56,8 +56,8 @@ cfg$force_download <- FALSE
 cfg$force_replace <- FALSE
 
 # Settings for the yield calibration
-# For this yield calibration (that uses results from a special MAgPIE calibration run 
-# triggered by the recalibrate switch) to be activated the switch s14_use_yield_calib needs 
+# For this yield calibration (that uses results from a special MAgPIE calibration run
+# triggered by the recalibrate switch) to be activated the switch s14_use_yield_calib needs
 # to be turned on. (Default is off)
 # This switch should only be activated for penalty_apr22 crop realization.
 # For other realizations, it is recommended not to use this yield calibration.
@@ -325,7 +325,7 @@ cfg$gms$s14_degradation <- 0      # def = 0
 
 # Switch to toggle the use of yield calibration factors (that resulted from a calibration run).
 # If 0, no yield calibration factors are used, meaning all calibration factors are set to 1.
-# If 1, yield calibration factors are used. (For this option to function, either an existing calibration file 
+# If 1, yield calibration factors are used. (For this option to function, either an existing calibration file
 # must be supplied in the input directory, or yields must be recalibrated during preprocessing.)
 # This switch should only be activated for penalty_apr22 crop realization.
 # For other realizations, it is recommended not to use this yield calibration.
@@ -421,16 +421,25 @@ cfg$gms$s15_exo_waste <- 0               # def = 0
 # * (1.2):  corresponds to 20% food waste ~ half waste of HIC
 cfg$gms$s15_waste_scen <- 1.2          # def = 1.2
 
-# * switch for transition to exogenous diet scenarios (EAT Lancet and National Institute of Nutrition (NIN))
-# * (2):      transition towards exogenous diets (NIN for India and EAT for other regions)
-# *           Please use this switch in combination with cfg$gms$c15_EAT_scen <- "FLX_hmilk"
-# *           If only India's diet needs to be changed, then set cfg$gms$scen_countries15 <- "IND"
-# * (1):      transition towards exogenous diets and food demand
-# * (0):      regression-based estimation of diets and food demand
+# * Switch for transition to exogenous diet scenarios
+# * (EAT Lancet and National Institute of Nutrition (NIN))
+# * (3):  MAgPIE-specific realization of the EAT-Lancet diet, where model-internal
+# *       (regression-based) food demand projections are constrained by ranges for
+# *       intake targets of food groups to ensure healthy and sustainable diets as
+# *       recommended by the EAT-Lancet Commission
+# * (2):  transition towards exogenous diets (NIN for India and EAT for other regions)
+# *       Please use this switch in combination with cfg$gms$c15_EAT_scen <- "FLX_hmilk"
+# *       If only India's diet needs to be changed, then set cfg$gms$scen_countries15 <- "IND"
+# * (1):  transition towards exogenous diets and food demand parametrized
+# *       to a food-specific data set published by the EAT-Lancet Commission
+# *       (Willett et al., 2019). This EATLancet implementation is deprecated.
+# *       Please use MAgPIE-specific realization of the EAT-Lancet diet (3)
+# *       for EATLancet diet.
+# * (0):  regression-based estimation of diets and food demand
 cfg$gms$s15_exo_diet <- 0               # def = 0
 
 # * exogenous calorie scenario (EAT Lancet diet scenarios)
-# * only activated if s15_exo_diet is set to 1
+# * only activated if s15_exo_diet is > 0
 # *   options:   healthy_BMI, 2100kcal, 2500kcal,
 # *              endo, no_underweight, no_overweight
 # *              half_overweight, no_underweight_half_overweight
@@ -445,35 +454,36 @@ cfg$gms$s15_exo_diet <- 0               # def = 0
 cfg$gms$c15_kcal_scen <- "healthy_BMI"       # def = healthy_BMI
 #
 # * exogenous food-specific diet scenario (EAT Lancet diet scenarios)
-# * only activated if s15_exo_diet is set to 1
+# * only activated if s15_exo_diet is > 0
 # *   options:   BMK, FLX, PSC, VEG, VGN, FLX_hmilk, FLX_hredmeat
 cfg$gms$c15_EAT_scen <- "FLX"                # def = FLX
 
-# * Subspecifications: Which commodities shall be included in the diet shift
+# * Sub-specifications: Which commodities shall be included in the diet shift
 # * towards cfg$gms$c15_EAT_scen ? Only selected (1) commodities will be
 # * included, while the others will remain to be endogenous.
 # * Any shift will be compensated to reach the calorie target by adjusting
 # * staple calories.
-cfg$gms$s15_exo_monogastric <- 1   # def = 1, options: 0,1
-cfg$gms$s15_exo_ruminant    <- 1   # def = 1, options: 0,1
-cfg$gms$s15_exo_fish        <- 1   # def = 1, options: 0,1
-cfg$gms$s15_exo_fruitvegnut <- 1   # def = 1, options: 0,1
-cfg$gms$s15_exo_pulses      <- 1   # def = 1, options: 0,1
-cfg$gms$s15_exo_sugar       <- 1   # def = 1, options: 0,1
-cfg$gms$s15_exo_oils        <- 1   # def = 1, options: 0,1
-cfg$gms$s15_exo_brans       <- 1   # def = 1, options: 0,1
-cfg$gms$s15_exo_scp         <- 1   # def = 1, options: 0,1
+cfg$gms$s15_exo_monogastric  <- 1   # def = 1, options: 0,1
+cfg$gms$s15_exo_ruminant     <- 1   # def = 1, options: 0,1
+cfg$gms$s15_exo_fish         <- 1   # def = 1, options: 0,1
+cfg$gms$s15_exo_fruitvegnutroots  <- 1   # def = 1, options: 0,1
+cfg$gms$s15_exo_pulses       <- 1   # def = 1, options: 0,1
+cfg$gms$s15_exo_sugar        <- 1   # def = 1, options: 0,1
+cfg$gms$s15_exo_oils         <- 1   # def = 1, options: 0,1
+cfg$gms$s15_exo_brans        <- 0   # def = 0 (as no target defined in EATLancet. with 1: brans are set to 0), options: 0,1
+cfg$gms$s15_exo_scp          <- 1   # def = 1, options: 0,1
 
 # * Scenario target for the inclusion of alcohol in the EAT-Lancet diet
-# * The EAT-Lancet diet only allows for added sugars, but does not include processed food or
-# * alcohol.
-# * only activated if s15_exo_diet is set to 1
-cfg$gms$s15_exo_alcohol <- 1 # default: 1
+# * The EAT-Lancet diet only allows for added sugars,
+# * but does not include processed food or alcohol.
+# * Only activated if s15_exo_diet is set to 1
+cfg$gms$s15_exo_alcohol <- 1        # def = 1, options: 0,1
 # * Via 's15_alc_scen' a maximum target for alcohol consumption can be defined.
+# * Only activated if s15_exo_alcohol = 1
 # * (0):      no alcohol consumption, as in the original version of the EAT-Lancet diet
 # * (0.014):  maximum target for alcohol consumption is 1.4% of total calorie consumption
 # *           (see Lassen et al., 2020)
-cfg$gms$s15_alc_scen <- 0.014                # def = 0.014
+cfg$gms$s15_alc_scen <- 0               # def = 0
 #
 #
 # * Switch and specification of countries for which exogenous food scenarios
@@ -772,11 +782,11 @@ cfg$gms$s32_max_aff_area <- Inf              # def = Inf
 # Protection of afforested areas (endogenous CO2 price driven)
 # 0=until end of planning horizon (see s32_planing_horizon)
 # 1=forever
-# Note: Without a price on CO2 emissions from land-use change in module 32_forestry, 
-#       which is the current default (c56_emis_policy <- "reddnatveg_nosoil"), 
-#       the recommended setting for s32_aff_prot is 1. 
-#       Otherwise, the model can do re/afforestation on existing re/afforestation areas 
-#       after the end of planning horizon without penalty for CO2 emissions. 
+# Note: Without a price on CO2 emissions from land-use change in module 32_forestry,
+#       which is the current default (c56_emis_policy <- "reddnatveg_nosoil"),
+#       the recommended setting for s32_aff_prot is 1.
+#       Otherwise, the model can do re/afforestation on existing re/afforestation areas
+#       after the end of planning horizon without penalty for CO2 emissions.
 cfg$gms$s32_aff_prot <- 1              # def = 1
 
 # Type of afforestation constraint
@@ -1271,7 +1281,7 @@ cfg$gms$c50_scen_neff_noselect <- "baseeff_add3_add5_add10_max65"   # def = base
 
 # * For pastures, only two scenarios exist: a constant scenario, and one where
 # * the minimum NUE increases to 55,60 and 65% for 2020,2050 and 2100, to avoid
-# * that pastures which dont have a nutrient deficit are fertilized
+# * that pastures which do not have a nutrient deficit are fertilized
 # * (historical NUE may be low despite no fertilization)
 # * Options:
 # * constant, constant_min55_min60_min65
@@ -1364,7 +1374,7 @@ cfg$gms$s56_cprice_red_factor <- 1    # def = 1
 # * Note: For best consistency it is recommended to use trajectories from the most recent
 # *     coupled REMIND-MAgPIE runs. Currently, this is R32M46.
 # * Available options:
-# *     none: no 2nd generation bioenergy demand 
+# *     none: no 2nd generation bioenergy demand
 # * R32M46: Coupled REMIND-MAgPIE runs with REMIND 3.2.1.dev412 and MAgPIE 4.6.11
 # *   NPi: Current policies; limit peak warming to 3.0°C
 # *   NDC: Nationally determined contributions; limit peak warming to 2.5°C
@@ -1463,9 +1473,9 @@ cfg$gms$c56_cprice_aff <- "secdforest_vegc"   # def = secdforest_vegc
 # * should not be higher than s32_planing_horizon
 cfg$gms$s56_c_price_exp_aff <- 50   # def = 50
 
-# * Share of carbon credits for re/afforestation projects set aside in a buffer to 
+# * Share of carbon credits for re/afforestation projects set aside in a buffer to
 # * account for potential failure of permanent carbon storage.
-# * Alternatively, this parameter can be interpreted as a reduction factor 
+# * Alternatively, this parameter can be interpreted as a reduction factor
 # * of the carbon price for non-permanent CDR.
 # * Values > 0 will reduce the incentive for c-price-induced re/afforestation
 cfg$gms$s56_buffer_aff <- 0.5   # def = 0.5
@@ -1548,10 +1558,10 @@ cfg$gms$s57_maxmac_ch4_awms  <- -1   # def = -1
 
 # ***---------------------    58_peatland   ------------------------------------
 # * (off):  Peatland area and associated GHG emissions are assumed zero
-# * (v2):   Peatland area is initialized with degraded and intact peatland area for the year 2022,  
+# * (v2):   Peatland area is initialized with degraded and intact peatland area for the year 2022,
 # *         and projected into the future based on changes in managed land.
-# *         For boreal and tropical climates, GHG emissions are calculated using emission 
-# *         factors from the IPCC 2013 Wetland supplement. For temperate climates, more recent 
+# *         For boreal and tropical climates, GHG emissions are calculated using emission
+# *         factors from the IPCC 2013 Wetland supplement. For temperate climates, more recent
 # *         estimates from Tiemeyer et al 2020 are used.
 cfg$gms$peatland  <- "v2"                     # def = v2
 
@@ -1572,8 +1582,8 @@ cfg$gms$s58_cost_drain_rewet_onetime  <- 0   # def = 0
 cfg$gms$s58_cost_drain_recur  <- 0   # def = 0
 
 
-# * Switch for fixing peatland area until the year given by s58_fix_peatland to historic levels (not available in `off`). 
-# * Note: In the realisation `v2`, historic peatland area is only available for the year 2022. 
+# * Switch for fixing peatland area until the year given by s58_fix_peatland to historic levels (not available in `off`).
+# * Note: In the realisation `v2`, historic peatland area is only available for the year 2022.
 # *     Due to the 5-year time step limitation in MAgPIE, the suggested default is 2020.
 # *     Fixing  peatland area in previous time steps provides a better
 # *     proxy for GHG emissions from peatlands than assuming no peatland area.
@@ -1625,7 +1635,7 @@ cfg$gms$c60_1stgen_biodem <- "const2020"               # def = const2020
 # * Note: For best consistency it is recommended to use trajectories from the most recent
 # *     coupled REMIND-MAgPIE runs. Currently, this is R32M46.
 # * Available options:
-# *     none: no 2nd generation bioenergy demand 
+# *     none: no 2nd generation bioenergy demand
 # * R32M46: Coupled REMIND-MAgPIE runs with REMIND 3.2.1.dev412 and MAgPIE 4.6.11
 # *   NPi: Current policies; limit peak warming to 3.0°C
 # *   NDC: Nationally determined contributions; limit peak warming to 2.5°C

config/scenario_config.csv

@@ -23,6 +23,7 @@ gms$s15_exo_oils;;;;;;;;;;1;0;1;1;;;;;;;1;1;1;1;1;1;1;;;;;;;;;;;;;;;;;;;
 gms$s15_exo_brans;;;;;;;;;;1;0;1;1;;;;;;;1;1;1;1;1;1;1;;;;;;;;;;;;;;;;;;;
 gms$s15_exo_scp;;;;;;;;;;1;0;1;1;;;;;;;1;1;1;1;1;1;1;;;;;;;;;;;;;;;;;;;
 gms$s15_exo_alcohol;;;;;;;;;;1;0;1;1;;;;;;;1;1;1;1;1;1;1;;;;;;;;;;;;;;;;;;;
+gms$s15_alc_scen;0.014;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 gms$c21_trade_liberalization;;;;l908080r807070;l909090r808080;l909090r808080;l909595r809090;l908080r807070;l908080r807070;l908080r807070;l908080r807070;l909595r809090;l908080r807070;;;;;;;;;l908080r807070;l908080r807070;l908080r807070;l908080r807070;l908080r807070;;;;;;;;;;;;;;;;;;;
 gms$c22_protect_scenario;;;;none;none;none;none;none;none;BH;none;BH_IFL;BH;;;;;;;BH;;BH_IFL;BH_IFL;BH_IFL;BH_IFL;BH_IFL;;;;;;;;;;;;;;;;;;;
 gms$c22_protect_scenario_noselect;;;;;;;;;;;;;;;;;;;;;;none;none;BH_IFL;none;BH_IFL;;;;;;;;;;;;;;;;;;;

config/scenario_fsec.csv

@@ -19,6 +19,7 @@ gms$s15_exo_sugar;;;;;;;;0;0;0;0;0;0;1;0;0;;;;;;;;;;;;;;;;;1;1;;;;;;;;;;;;;;;
 gms$s15_exo_oils;;;;;;;;0;0;0;0;0;0;1;0;0;;;;;;;;;;;;;;;;;1;1;;;;;;;;;;;;;;;
 gms$s15_exo_scp;;;;;;;;0;0;0;0;0;0;1;0;0;;;;;;;;;;;;;;;;;1;1;;;;;;;;;;;;;;;
 gms$s15_exo_alcohol;;;;;;;;0;0;0;0;0;0;1;0;0;;;;;;;;;;;;;;;;;1;1;;;;;;;;;;;;;;;
+gms$s15_alc_scen;0.014;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 gms$c15_food_scenario;;;SSP1;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 gms$c21_trade_liberalization;;;;;;;;;;;;;;;;;;;;;;;l908080r807070;;;;;;;;;;;;;;;;;;;;;;;;;;
 gms$c22_protect_scenario;;;;;;;;;;;;;;;;;;;;;;;;;;;;;BH_IFL;;;;;;;;;;;;;;;;;;;;
@@ -65,9 +66,9 @@ gms$s62_max_dem_bioplastic;0;;;;400;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 gms$c70_fac_req_regr;reg;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 gms$c70_feed_scen;;;;;;;;;;;;;;;;;;ssp1;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 gms$c73_build_demand;;;;;;;;;;;;;;;;;;;;;;;;50pc;;;;;;;;;;;;;;;;;;;;;;;;;
-input['cellular'];rev4.96_FSEC_3c888fa5_cellularmagpie_c200_MRI-ESM2-0-ssp460_lpjml-8e6c5eb1.tgz;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;rev4.96_FSEC_0bd54110_cellularmagpie_c200_MRI-ESM2-0-ssp119_lpjml-8e6c5eb1.tgz;rev4.96_FSEC_6819938d_cellularmagpie_c200_MRI-ESM2-0-ssp126_lpjml-8e6c5eb1.tgz;;rev4.96_FSEC_1b5c3817_cellularmagpie_c200_MRI-ESM2-0-ssp245_lpjml-8e6c5eb1.tgz;rev4.96_FSEC_3c888fa5_cellularmagpie_c200_MRI-ESM2-0-ssp460_lpjml-8e6c5eb1.tgz;rev4.96_FSEC_fd712c0b_cellularmagpie_c200_MRI-ESM2-0-ssp370_lpjml-8e6c5eb1.tgz;rev4.96_FSEC_09a63995_cellularmagpie_c200_MRI-ESM2-0-ssp585_lpjml-8e6c5eb1.tgz;
-input['regional'];rev4.96_FSEC_magpie.tgz;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-input['validation'];rev4.96_FSEC_validation.tgz;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+input['cellular'];rev4.99_FSEC_3c888fa5_cellularmagpie_c200_MRI-ESM2-0-ssp460_lpjml-8e6c5eb1.tgz;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;rev4.99_FSEC_0bd54110_cellularmagpie_c200_MRI-ESM2-0-ssp119_lpjml-8e6c5eb1.tgz;rev4.99_FSEC_6819938d_cellularmagpie_c200_MRI-ESM2-0-ssp126_lpjml-8e6c5eb1.tgz;;rev4.99_FSEC_1b5c3817_cellularmagpie_c200_MRI-ESM2-0-ssp245_lpjml-8e6c5eb1.tgz;rev4.99_FSEC_3c888fa5_cellularmagpie_c200_MRI-ESM2-0-ssp460_lpjml-8e6c5eb1.tgz;rev4.99_FSEC_fd712c0b_cellularmagpie_c200_MRI-ESM2-0-ssp370_lpjml-8e6c5eb1.tgz;rev4.99_FSEC_09a63995_cellularmagpie_c200_MRI-ESM2-0-ssp585_lpjml-8e6c5eb1.tgz;
+input['regional'];rev4.99_FSEC_magpie.tgz;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+input['validation'];rev4.99_FSEC_validation.tgz;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 input['additional'];additional_data_rev4.47.tgz;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 input['calibration'];calibration_FSEC_24Mar23.tgz;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 magicc_emis_scen;bjoernAR6_C_SSP2-NDC.mif;;;bjoernAR6_C_SSP2-PkBudg900.mif;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;bjoernAR6_C_SSP1-NDC.mif;;;;;;;;;;;;bjoernAR6_C_SSP1-PkBudg900.mif