Extract EMEP atmospheric deposition data for a catchment

Usage

get_atmo_dep(
  catchment_boundary_path,
  t_ext = "year",
  start_year = 1990,
  end_year = 2020
)

Arguments

catchment_boundary_path: path to basin boundary shape file.
t_ext: string, which EMEP data to access 'year' for yearly averages, 'month' - monthly averages. Optional (default - 'year').
start_year: integer year to start data extraction. Optional (default - 1990).
end_year: integer year to end data extraction. Optional (default - 2020).

Value

dafaframe with "DATE", "NH4_RF", "NO3_RF" , "NH4_DRY" and "NO3_DRY" columns. Values in SWAT+ units. "NH4_RF" - ammonia in rainfall (mg/l), "NO3_RF" - nitrate in rainfall (mg/l), NH4_RF - ammonia deposition (kg/ha/yr), "NO3_DRY" - nitrate dry deposition (kg/ha/yr)

Examples

if (FALSE) {
basin_path <- system.file("extdata", "GIS/basin.shp", package = "svatools")
df <- get_atmo_dep(basin_path)
##Plot results
ggplot(pivot_longer(df, !DATE, names_to = "par", values_to = "values"), aes(x = DATE, y = values))+
geom_line()+
facet_wrap(~par, scales = "free_y")+
theme_bw()
}