Merge pull request #672 from FelicitasBeier/updateMehta2024

Updated AEI initialization to Mehta (2024)
magpiemodel · May 28, 2024 · 77d5df1 · 77d5df1
2 parents 58cc585 + 01fa7a3
commit 77d5df1
Show file tree

Hide file tree

Showing 14 changed files with 466 additions and 209 deletions.
diff --git a/CHANGELOG.md b/CHANGELOG.md
@@ -7,6 +7,7 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/).
 ## [Unreleased]
 
 ### changed
+- **41_area_equipped_for_irrigation** updated (non-default) AEI data (from Mehta2022 to Mehta2024)
 - **22_land_conservation and default.cfg** Added options for baseline protection
 - **15_food, default.cfg and scenario_config.csv** changed fader setup and introduced new switches for specifying food substitution scenarios and exogeneous food intake scenarios
 - **70_livestock** default.cfg and scenario_config.csv** changed fader setup and introduced new switches for specifying feed substitution with SCP scenarios

diff --git a/config/default.cfg b/config/default.cfg
@@ -1087,7 +1087,8 @@ cfg$gms$area_equipped_for_irrigation <- "endo_apr13"    # def = endo_apr13
 
 # * switch for initialization area
 # * (LUH2v2): area equipped for irrigation based on LUH2v2 irrigated areas
-# * (Mehta2022): area equipped for irrigation from Mehta et al. (2022)
+# * (Mehta2024_Siebert2013): area equipped for irrigation from Mehta et al. (2024) using GMIA (Siebert et al, 2013) base map
+# * (Mehta2024_Meier2018): are equipped for irrigation from Mehta et al. (2024) using Meier et al. (2018) base map
 cfg$gms$c41_initial_irrigation_area  <- "LUH2v2"        # def = LUH2v2
 
 # * Sets the rate of depreciation of irrigation infrastructure in every timestep.

diff --git a/config/scenario_config.csv b/config/scenario_config.csv
diff --git a/literature.bib b/literature.bib
@@ -1023,14 +1023,47 @@ @article{siebert_FAO_2007
   year = {2007},
 }
 
+@article{siebert_FAO_2007,
+  author = {Stefan Siebert and Petra Döll and Sebastian Feick and Jippe Hoogeveen and Karen Frenken},
+  year = {2007},
+  title = {Global Map of Irrigation Areas version 4.0.1.},
+  journal = {Johann Wolfgang Goethe University, Frankfurt am Main, Germany / Food and Agriculture Organization of the United Nations, Rome, Italy.},
+  url = {http://www.fao.org/nr/water/aquastat/irrigationmap/index10.stm},
+}
+
+@article{siebert_FAO_2013,
+  author = {Stefan Siebert, Verena Henrich, Karen Frenken and Jacob Burk},
+  year = {2013},
+  title = {Global Map of Irrigation Areas version 5},
+  journal = {Rheinische Friedrich-Wilhelms-University, Bonn, Germany / Food and Agriculture Organization of the United Nations, Rome, Italy},
+  url = {https://data.apps.fao.org/catalog/iso/f79213a0-88fd-11da-a88f-000d939bc5d8},
+}
+
+@article{meier_global_2018,
+  author = {Meier, Jonas; Zabel, Florian; Mauser, Wolfram},
+  year = {2018},
+  title = {Global Irrigated Areas dataset},
+  journal = {PANGAEA},
+  doi = {https://doi.org/10.1594/PANGAEA.884744},
+  url = {https://doi.pangaea.de/10.1594/PANGAEA.884744},
+}
+
+@article{mehta_half_2024,
+  author = {Mehta, P., Siebert, S., Kummu, M. et al.},
+  year = {2024},
+  title = {Half of twenty-first century global irrigation expansion has been in water-stressed regions},
+  journal = {Nature Water},
+  doi = {https://doi.org/10.1038/s44221-024-00206-9},
+}
+
 @book{fao_aquastat_2016,
   address = {Rome},
   title = {{AQUASTAT} core database},
   publisher = {Food and Agriculture Organization of the United Nations (FAO)},
   author = {{FAO}},
   url = {{https://data.harvestportal.org/de/dataset/fao-aquastat/resource/c4668555-eb76-4882-83b1-230038e24f02?inner_span=True}},
   year = {2016},
-  note = {Database accessed on 2023/02/06}
+  note = {Database accessed on 2023/02/06},
 }
 
 @techreport{worldbank_irrigation_1995,

diff --git a/modules/37_labor_prod/exo/sets.gms b/modules/37_labor_prod/exo/sets.gms
@@ -6,7 +6,7 @@
 *** |  Contact: [email protected]
 
 sets
- rcp37 climate change impact sceanrios
+ rcp37 climate change impact scenarios
     / rcp119, rcp585 /
 
  metric37 labor productivity impact metric

diff --git a/modules/41_area_equipped_for_irrigation/endo_apr13/realization.gms b/modules/41_area_equipped_for_irrigation/endo_apr13/realization.gms
@@ -7,8 +7,9 @@
 
 
 *' @description This realization allows the model to endogenously decide on investments to deploy additional
-*' irrigation infrastructure, i.e. to increase the area equipped for irrigation (AEI). Initial values for AEI
-*' in 1995 are taken from @siebert_FAO_2007. Contraction of AEI is not possible.
+*' irrigation infrastructure, i.e. to increase the area equipped for irrigation (AEI).
+*' Initial values for AEI are areas that are irrigated in the initial time step in LUH2v2.
+*' Depreciation of AEI can be set via `s41_AEI_depreciation`.
 *' Irrigated crop production can only take place where irrigation infrastructure is present.
 *'
 *' Unit costs per hectare for AEI expansion are derived from a World Bank study (@worldbank_irrigation_1995)

diff --git a/modules/41_area_equipped_for_irrigation/endo_apr13/sets.gms b/modules/41_area_equipped_for_irrigation/endo_apr13/sets.gms
@@ -10,5 +10,5 @@ sets
        / y1995, y2000, y2005, y2010, y2015 /
 
   aei41    Data source for area equipped for irrigation
-       / LUH2v2, Mehta2022 /
+       / LUH2v2, Mehta2024_Siebert2013, Mehta2024_Meier2018 /
 ;
diff --git a/modules/41_area_equipped_for_irrigation/static/realization.gms b/modules/41_area_equipped_for_irrigation/static/realization.gms
@@ -7,7 +7,7 @@
 
 
 *' In this realization, area equipped for irrigation is fixed to input data (around the year 2000)
-*' for all time steps. The source of the input data is @siebert_FAO_2007.
+*' for all time steps. The sources of the input data are @siebert_FAO_2013 and @mehta_half_2024.
 *'
 *' @limitations No irrigation is possible on areas that have not been equipped for irrigation in the past.
 *'

diff --git a/modules/41_area_equipped_for_irrigation/static/sets.gms b/modules/41_area_equipped_for_irrigation/static/sets.gms
@@ -10,5 +10,5 @@ sets
        / y1995, y2000, y2005, y2010, y2015 /
 
   aei41    Data source for area equipped for irrigation
-       / LUH2v2, Mehta2022 /
+       / LUH2v2, Mehta2024_Siebert2013, Mehta2024_Meier2018 /
 ;
diff --git a/modules/42_water_demand/agr_sector_aug13/realization.gms b/modules/42_water_demand/agr_sector_aug13/realization.gms
@@ -33,6 +33,7 @@
 *'   efficiency times management factor") from @PIK_report104_2007.
 *'   Here, irrigated area from @siebert_FAO_2007 has been used as aggregation
 *'   weight.
+*'   Contraction of AEI happens if a depreciation rate is set in the switch `s41_AEI_depreciation`. 
 *'
 *' * A regression of country values of the "conveyance efficiency times
 *'   management factor" from @PIK_report104_2007 on GDP.

diff --git a/scripts/output/projects/FSEC_water.R b/scripts/output/projects/FSEC_water.R
@@ -59,3 +59,7 @@ watStressViolations <- watStress
 # mark violations in different color
 watStressViolations[efvViolation == 1] <- 100
 write.magpie(watStressViolations, file_name = file.path(outputdir, "watStressViolations.mz"))
+
+# Water EFV ratio (EFV to EFR)
+watEFVratio <- waterEFVratio(gdx, level = "grid", dir = outputdir)
+write.magpie(watEFVratio, file_name = file.path(outputdir, "watEFVratio.mz"))
diff --git a/scripts/output/projects/inms_reporting_reg.R b/scripts/output/projects/inms_reporting_reg.R
@@ -24,15 +24,15 @@ print("Start inms reporting reg runscript")
 
 if(!exists("source_include")) {
 
-  title       <- "inms_SSP2_RCP4p5_PolicyLow_v4"
-  outputdir       <- "output/inms_SSP2_RCP4p5_PolicyLow_v4_2020-07-13_15.37.07"
+  title     <- "inms_SSP2_RCP4p5_PolicyLow_v4"
+  outputdir <- "output/inms_SSP2_RCP4p5_PolicyLow_v4_2020-07-13_15.37.07"
 
   ###Define arguments that can be read from command line
-  readArgs("outputdir","title")
+  readArgs("outputdir", "title")
 }
 #########################################################################################
 
-print(paste0("script started for output directory",outputdir))
+print(paste0("script started for output directory", outputdir))
 
 withr::local_dir(outputdir)
 
@@ -41,24 +41,26 @@ title <- cfg$title
 print("generating INMS output for the run: ")
 print(title)
 
-filename=paste0("report_",title,".mif")
-gdx=paste0("fulldata.gdx")
-a=getReportINMS(gdx,file=filename,scenario = title,dir=".")
+filename <- paste0("report_", title, ".mif")
+gdx <- paste0("fulldata.gdx")
+a <- getReportINMS(gdx, file = filename, scenario = title, dir = ".")
 
 print(filename)
-mif=read.report(filename)
+mif <- read.report(filename)
 
 
-missingyears=function(x){
-  history=paste0("y",1965+((0:5)*5))
-  x[[1]][[1]] <- time_interpolate(x[[1]][[1]],interpolated_year = c(history,paste0("y",2005+((0:9)*10))),integrate_interpolated_years = TRUE)
-  x[[1]][[1]][,history,] = 0
+missingyears <- function(x) {
+  history <- paste0("y", 1965 + ((0:5) * 5))
+  x[[1]][[1]] <- time_interpolate(x[[1]][[1]],
+                                  interpolated_year = c(history,paste0("y",2005+((0:9)*10))),
+                                  integrate_interpolated_years = TRUE)
+  x[[1]][[1]][, history, ] <- 0
   return(x)
 }
 
 #a=c(missingyears(ssp1),missingyears(ssp2))
-a=missingyears(mif)
+a <- missingyears(mif)
 
-write.reportProject(a,mapping=paste0(wdbefore,"/mapping_inms.csv"),file="report_inms.mif")
+write.reportProject(a, mapping = paste0(wdbefore, "/mapping_inms.csv"), file = "report_inms.mif")
 #write.report(a,file="magpie_results_nov2019.mif")
 warnings()