Skip to content

embodiedWater

Source: R/embodiedWater.R
embodiedWater.Rd

Calculates production-based and consumption-based (embodied) water footprint accounting using bilateral trade flows. Water use is allocated to traded products based on production ratios and bilateral trade patterns.

Usage

embodiedWater(
  gdx,
  file = NULL,
  level = "reg",
  type = "all",
  waterType = "consumption",
  bilateral = FALSE,
  disaggLivestock = FALSE
)

Arguments

gdx

GDX file

file

a file name the output should be written to using write.magpie

level

Level of regional aggregation; "reg" (regional), "glo" (global), "regglo" (regional and global) or any other aggregation level defined in superAggregate. Only used when bilateral=FALSE.

type

Type of accounting: "production" (production-based), "consumption" (consumption-based), "trade" (export, import, and net-trade), "all" (all five), or "flows" (bilateral flows, requires bilateral=TRUE)

waterType

Type of water to report: "withdrawal" (water withdrawal), "consumption" (consumptive water use)

bilateral

Logical; if TRUE, returns bilateral flows with dimensions (exporter.importer, year, product) instead of regional totals (default FALSE)

disaggLivestock

Logical; if TRUE, the feed pathway retains the livestock product dimension, so water is attributed per animal product × feed crop combination. Passes disaggLivestock to tradedPrimariesBilateral. Use dimSums(x[kli_items], dim=3.1) to collapse to feed crops, or dimSums(x[kli_items], dim=3.2) to collapse to animal products. Default is FALSE (current behaviour: feed attributed to crops).

Value

Embodied water use as MAgPIE object (unit depends on waterType). When bilateral=FALSE: dimensions are (region, year, accounting.product). When bilateral=TRUE: dimensions are (exporter.importer, year, pathway.product). When disaggLivestock=TRUE: feed items have dimensions kli_product.kve instead of feed.kve; prim/secd items and all aggregated accounting outputs are unchanged.

See also

water_usage, trade

Author

David M Chen

Examples

if (FALSE) { # \dontrun{
  x <- embodiedWater(gdx, type = "all", waterType = "consumption")
  # Bilateral flows
  xBilat <- embodiedWater(gdx, type = "flows", bilateral = TRUE)
} # }