update

config/default.cfg

@@ -1760,6 +1760,11 @@ cfg$gms$s58_rewet_exo_start_value <- 0.3   # def = 0.3
 # * Target value for exogenous peatland rewetting as share of drained peatland in reference period
 cfg$gms$s58_rewet_exo_target_value <- 0.5   # def = 0.5
 
+# * Switch for exogenous protection of intact peatland for selected countries (0=off, 1=on)
+cfg$gms$s58_intact_prot_exo <- 0   # def = 0
+# * Switch for exogenous protection of intact peatland for all other countries (0=off, 1=on)
+cfg$gms$s58_intact_prot_exo_noselect <- 0   # def = 0
+
 # * One-time and recurring costs for peatland rewetting (USD17MER per ha)
 # * (inflated from default originally in USD05 using USD05 --> USD17 inflation rate:1.23)
 cfg$gms$s58_cost_rewet_onetime  <- 1230   # def = 1000 * 1.23

modules/58_peatland/v2/declarations.gms

@@ -22,7 +22,8 @@ parameters
  p58_peatland_ref(j,land58)          Peatland area in reference period (mio. ha)
  p58_country_dummy(iso)              Dummy parameter indicating whether country is affected by exogenous peatland rewetting (1)
  p58_country_weight(i)               Policy country weight per region (1)
- i58_rewetting_exo(t,j)              Share of drained peatland in reference period to be rewetted (1)
+ i58_rewetting_exo(t,j)              Cellular share of drained peatland in reference period to be rewetted (1)
+ i58_intact_protection_exo(j)        Cellular share of intact peatland that should be protected (1)
 ;
 
 equations

modules/58_peatland/v2/input.gms

@@ -16,10 +16,12 @@ scalars
   s58_balance_penalty            Penalty for technical peatland balance term (USD17MER) / 1e+06 /
   s58_rewetting_exo              Switch for exogenous peatland rewetting for selected countries (1) / 0 /
   s58_rewetting_exo_noselect     Switch for exogenous peatland rewetting for all other countries (1) / 0 /
-  s58_rewet_exo_start_year       Start year for exogenous peatland rewetting / 2030 /
-  s58_rewet_exo_target_year      Target year for exogenous peatland rewetting / 2050 /
+  s58_rewet_exo_start_year       Start year for exogenous peatland rewetting (1) / 2030 /
+  s58_rewet_exo_target_year      Target year for exogenous peatland rewetting (1) / 2050 /
   s58_rewet_exo_start_value      Start value for exogenous peatland rewetting as share of drained peatland in reference period (1) / 0.3 /
   s58_rewet_exo_target_value     Target value for exogenous peatland rewetting as share of drained peatland in reference period (1) / 0.5 /
+  s58_intact_prot_exo            Switch for exogenous protection of intact peatland for selected countries (1) / 0 /
+  s58_intact_prot_exo_noselect   Switch for exogenous protection of intact peatland for all other countries (1) / 0 /
 ;
 
 * Set-switch for countries affected by certain policies

modules/58_peatland/v2/preloop.gms

@@ -23,6 +23,9 @@ i58_rewetting_exo(t,j) = i58_peatland_rewetting_fader(t) *
   (s58_rewetting_exo * sum(cell(i,j), p58_country_weight(i))
   + s58_rewetting_exo_noselect * sum(cell(i,j), 1-p58_country_weight(i)));
 
+i58_intact_protection_exo(j) = s58_intact_prot_exo * sum(cell(i,j), p58_country_weight(i))
+                      + s58_intact_prot_exo_noselect * sum(cell(i,j), 1-p58_country_weight(i));
+
 * fix uncontrolled pollutants to zero
 vm_emissions_reg.fx(i,"peatland",pollutants) = 0;
 vm_emissions_reg.lo(i,"peatland",poll58) = -Inf;

modules/58_peatland/v2/presolve.gms

@@ -38,6 +38,7 @@ else
   v58_peatland.l(j,land58) = pc58_peatland(j,land58);
   v58_peatland.up(j,drained58) = Inf;
   v58_peatland.up(j,"rewetted") = s58_rewetting_switch;
+  v58_peatland.lo(j,"intact") = pc58_peatland(j,"intact") * i58_intact_protection_exo(j);
   v58_peatland.up(j,"intact") = pc58_peatland(j,"intact");
   v58_peatland.fx(j,"peatExtract") = pc58_peatland(j,"peatExtract");
   v58_balance.lo(j,manPeat58) = 0;

scripts/start/projects/paper_peatlandTax.R

@@ -23,7 +23,7 @@ source("scripts/start_functions.R")
 source("config/default.cfg")
 
 # create additional information to describe the runs
-cfg$info$flag <- "PTax43"
+cfg$info$flag <- "PTax44"
 
 cfg$results_folder <- "output/:title:"
 cfg$results_folder_highres <- "output"
@@ -153,6 +153,7 @@ cfg$title <- .title(cfg, paste(ssp, "Ref", sep = "-"))
 cfg$gms$c56_mute_ghgprices_until <- "y2150"
 cfg$gms$c56_pollutant_prices <- "T0-CO2"
 cfg$gms$s58_rewetting_exo <- 0
+cfg$gms$s58_intact_prot_exo <- 0
 start_run(cfg, codeCheck = FALSE)
 
 ## GHG policy scenarios
@@ -166,6 +167,7 @@ for (tax in c("T25-CO2",
   cfg$title <- .title(cfg, paste(ssp, tax, sep = "-"))
   cfg$gms$c56_mute_ghgprices_until <- "y2025"
   cfg$gms$c56_pollutant_prices <- tax
+  cfg$gms$s58_intact_prot_exo <- 1
   start_run(cfg, codeCheck = FALSE)
 }
 
@@ -174,16 +176,19 @@ cfg$title <- .title(cfg, paste(ssp, "NRL33", sep = "-"))
 cfg$gms$c56_mute_ghgprices_until <- "y2150"
 cfg$gms$c56_pollutant_prices <- "T0-CO2"
 cfg$gms$s58_rewetting_exo <- 0.33
+cfg$gms$s58_intact_prot_exo <- 1
 start_run(cfg, codeCheck = FALSE)
 
 cfg$title <- .title(cfg, paste(ssp, "NRL50", sep = "-"))
 cfg$gms$c56_mute_ghgprices_until <- "y2150"
 cfg$gms$c56_pollutant_prices <- "T0-CO2"
 cfg$gms$s58_rewetting_exo <- 0.5
+cfg$gms$s58_intact_prot_exo <- 1
 start_run(cfg, codeCheck = FALSE)
 
 cfg$title <- .title(cfg, paste(ssp, "NRL100", sep = "-"))
 cfg$gms$c56_mute_ghgprices_until <- "y2150"
 cfg$gms$c56_pollutant_prices <- "T0-CO2"
 cfg$gms$s58_rewetting_exo <- 1
+cfg$gms$s58_intact_prot_exo <- 1
 start_run(cfg, codeCheck = FALSE)