R/calcEmisNitrogenOceans.R

In pik-piam/mrcommons: MadRat commons Input Data Library

#' @title calcEmisNitrogenOceans
#' @description Conputes (nitrogen) emissions of Oceans.
#' @param method deposition inventory
#' @return List of magpie objects with results on country level, weight on country level, unit and description.
#' @author Benjamin Leon Bodirsky
#' @seealso
#' [calcAtmosphericDeposition()],
#' [calcNitrogenBudgetCropland()]
#' @examples
#' \dontrun{
#' calcOutput("EmisNitrogenOceans")
#' }
#'
calcEmisNitrogenOceans <- function(method = "ACCMIP") {

  if (method == "ACCMIP") {
    emi <- collapseNames(calcOutput("AtmosphericDeposition", datasource = "ACCMIP", glo_incl_oceans = FALSE,
                                    cellular = FALSE, emission = TRUE, aggregate = FALSE))
    emiGlo <- collapseNames(calcOutput("AtmosphericDeposition", datasource = "ACCMIP", glo_incl_oceans = TRUE,
                                       cellular = FALSE, emission = TRUE, aggregate = FALSE))
    out <- emi
    out[, , ] <- 0
    out["ATA", , ] <- emiGlo - dimSums(emi, dim = 1)
    out <- add_columns(out, dim = 3.1, addnm = c("n2o_n_direct", "accumulation"))
    out[, , c("n2o_n_direct", "accumulation")] <- 0
    # add global estimates for N2O and N burial from Voss et al (2013)
    # Voss, M., H. W. Bange, J. W. Dippner, J. J. Middelburg, J. P. Montoya, and B. Ward. 2013.
    # "The Marine Nitrogen Cycle: Recent Discoveries, Uncertainties and the Potential Relevance of Climate Change."
    # Philosophical Transactions of the Royal Society B: Biological Sciences 368 (1621): 20130121–20130121.
    # doi:10.1098/rstb.2013.0121. #nolint: commented_code_linter.
    vcat(2, "Values for oceans are inserted into the non-exisiting country Antarctica")
    out["ATA", , "n2o_n_direct"] <- 3.8
    out["ATA", , "accumulation"] <- 22
  } else if (method %in% c("Nsurplus", "Nsurplus2")) {
    baseyear <- "y2000"
    if (method == "Nsurplus2") {
      budget <- collapseNames(calcOutput("NitrogenBudgetOcean", deposition = "Nsurplus",
                                         aggregate = FALSE)[, , "surplus"])
    } else {
      budget <- collapseNames(calcOutput("NitrogenBudgetOcean", deposition = "CEDS", leaching = "Nsurplus",
                                         aggregate = FALSE)[, , "surplus"])
    }
    emi <- collapseNames(calcOutput("EmisNitrogenOceans", aggregate = FALSE, method = "ACCMIP")[, baseyear, ])
    out <- setYears(emi / dimSums(budget, dim = 1)[, baseyear, ], NULL) * dimSums(budget, dim = 1)
    out <- mbind(out, setNames(budget - dimSums(out, dim = 3), "n2_n"))
    out <- add_columns(out, dim = 3.1, addnm = c("no3_n"))
    out[, , "no3_n"] <- 0
  }

  out <- add_dimension(out, dim = 3.1, nm = "oceans")
  return(list(
    x = out,
    weight = NULL,
    unit = "Mt Nr, NH3N and NO2N",
    description = "Nitrogen emissions of oceans"))
}