R/reportSDG6.R

In pik-piam/magpie4: MAgPIE outputs R package for MAgPIE version 4.x

#' @title reportSDG6
#' @description reports all SDG indicators relevant for SDG6 - Access to Water
#'
#' @export
#'
#' @param gdx GDX file
#' @param level level of aggregation (cluster: "cell", regional: "regglo")
#' @param outputdir output directory
#'
#' @return MAgPIE object
#' @author Felicitas Beier, Isabelle Weindl
#' @import magclass
#' @examples
#'
#'   \dontrun{
#'     x <- reportSDG6(gdx)
#'   }
#'

reportSDG6 <- function(gdx, level = "regglo", outputdir = ".") {
  x <- NULL
  #cfg <- NULL
  #load(paste0(outputdir, "/config.Rdata"))

  indicatorname <- "SDG|SDG06|Safe sanitation"
  unit          <- "fraction"
  #missing
  # getNames(out) <- paste0(indicatorname, " (",unit,")")
  # x <- mbind(x,out)

  indicatorname <- "SDG|SDG06|Safe wastewater"
  unit          <- "fraction"
  #missing (retrieve from moinput?)
  # getNames(out) <- paste0(indicatorname, " (",unit,")")
  # x <- mbind(x,out)

  indicatorname <- "SDG|SDG06|Water quality"
  unit          <- "fraction"
  #missing
  # getNames(out) <- paste0(indicatorname, " (",unit,")")
  # x <- mbind(x,out)

  indicatorname <- "SDG|SDG06|N water loading"
  unit          <- "Mt N/yr"
  #missing
  # getNames(out) <- paste0(indicatorname, " (",unit,")")
  # runoff
  # x <- mbind(x,out)

  indicatorname <- "SDG|SDG06|P water loading"
  unit          <- "Mt P/yr"
  #missing
  # getNames(out) <- paste0(indicatorname, " (",unit,")")
  # x <- mbind(x,out)

  budget <- NitrogenBudget(gdx, level = "regglo")

  indicatorname <- "SDG|SDG06|Fertilizer use"
  unit          <- "Mt N/yr"
  # Def.: Nitrogen fertilizer use
  out <- budget[,,"fertilizer"]
  getNames(out) <- paste0(indicatorname, " (", unit, ")")
  x <- mbind(x,out)

  indicatorname <- "SDG|SDG06|Nitrogen surplus on cropland"
  unit          <- "Mt N/yr"
  # Def.: Nitrogen surplus on cropland
  out <- budget[,,"surplus"]
  getNames(out) <- paste0(indicatorname, " (", unit, ")")
  x <- mbind(x, out)

  indicatorname <- "SDG|SDG06|Nitrate concentration in water"
  unit          <- "tN/km3"
  #missing
  # getNames(out) <- paste0(indicatorname, " (",unit,")")
  # x <- mbind(x,out)

  indicatorname <- "SDG|SDG06|Water use efficiency"
  unit          <- "rate"
  #missing
  # Def.: a country's total gross domestic product (GDP) divided by total freshwater withdrawals (OurWorldInData: sdg-tracker.org????)
  # getNames(out) <- paste0(indicatorname, " (",unit,")")
  # x <- mbind(x,out)

  # indicatorname="SDG|SDG06|Water stress"
  # unit="fraction"
  # # Def.: total quantity of freshwater withdrawals (agriculture, manufacturing, domestic; km^3) as a share of total available freshwater resources (km^3)
  # # water use from MAgPIE
  # wateruse  <- water_usage(gdx,level=level,users=c("agriculture", "manufacturing", "electricity", "domestic"),digits=10) # unit: km^3/yr
  # nonaguses <- dimSums(wateruse[,,c("manufacturing","electricity","domestic")],dim=3)
  # wateruse  <- dimSums(wateruse,dim=3)
  # # total water availability (km^3)
  # waterav   <- read.magpie(paste0(outputdir,"/lpj_watavail_total_c200.mz"))/1000
  # years     <- intersect(getYears(wateruse),getYears(waterav))
  # if (cfg$gms$c43_watavail_scenario=="nocc") {
  #   waterav[,years,] <- waterav[,"y1995",]
  # }
  # waterav   <- gdxAggregate(gdx, waterav, weight="sum", to=level, absolute=TRUE)
  # # water scarcity
  # out           <- wateruse[,years,]/waterav[,years,]
  # # where non-agricultural uses exceed water availability (special rule in MAgPIE): scarcity indicator capped to 1
  # out[nonaguses[,years,]>waterav[,years,]] <- 1
  # getNames(out) <- paste0(indicatorname, " (",unit,")")
  # x <- mbind(x,out)
  #
  # indicatorname="SDG|SDG06|Agricultural water stress"
  # unit="fraction"
  # # Def.: quantity of agricultural freshwater withdrawals as a share of total available freshwater resources (km^3)
  # # water use from MAgPIE
  # wateruse  <- water_usage(gdx,level=level,users="agriculture",digits=10,sum=TRUE) # unit: km^3/yr
  # # total water availability (km^3)
  # waterav   <- read.magpie(paste0(outputdir,"/lpj_watavail_total_c200.mz"))/1000
  # years     <- intersect(getYears(wateruse),getYears(waterav))
  # if (cfg$gms$c43_watavail_scenario=="nocc") {
  #   waterav[,years,] <- waterav[,"y1995",]
  # }
  # waterav   <- gdxAggregate(gdx, waterav, weight="sum", to=level, absolute=TRUE)
  # # water scarcity
  # out           <- wateruse[,years,]/waterav[,years,]
  # getNames(out) <- paste0(indicatorname, " (",unit,")")
  # x <- mbind(x,out)

  indicatorname <- "SDG|SDG06|People under water stress"
  unit          <- "million"
  #missing (Def.: number of people living in water stressed region)
  # getNames(out) <- paste0(indicatorname, " (",unit,")")
  # x <- mbind(x,out)

  indicatorname <- "SDG|SDG06|Environmental flow exceedance"
  unit          <- "ratio"
  # Def.: ratio of violated environmental flows over total environmental flows
  #EFV <- water_EFexceedance(gdx,level=level,outputdir=outputdir,users=c("agriculture", "manufacturing", "electricity", "domestic"))

  #tmp <- EFV
  #tmp[EFV>0] <-0 # Cells where EFRs are not exceeded
  #tmp[EFV<0] <-(-1) # Cells that exceed EFRs

  #EFV_volume <- EFV*tmp # area affected by environmental flow violation
  #ratio_EFV  <- EFV_volume[,years,]/EFR[,years,]  # ratio of violation volume to total EFV
  #global_ratio <- dimSums(EFV_volume[,years,],dim=1)/dimSums(EFR[,years,],dim=1)
  #EFV_area <- superAggregate(EFV_area, aggr_type="sum",level=level,crop_aggr=TRUE) # area affected by environmental flow violation at reporting level


  #getNames(out) <- paste0(indicatorname, " (",unit,")")
  #x <- mbind(x,out)

  indicatorname <- "SDG|SDG06|Agricultural water use"
  unit          <- "km3/yr"
  # Def.: water usage in agriculture in the growing period
  out <- water_usage(gdx, level = "regglo", users = "agriculture", seasonality = "grper",
                     sum = TRUE, digits = 3)
  getNames(out) <- paste0(indicatorname, " (", unit, ")")
  x <- mbind(x,out)

  indicatorname <- "SDG|SDG06|Water-related ecosystems"
  unit          <- "million ha"
  #missing
  # getNames(out) <- paste0(indicatorname, " (",unit,")")
  # x <- mbind(x,out)

  #x <- x[,,sort(getNames(x))]
  return(x)
}