CalcNoahmpWatBudg: Calculate water balance from WRF-Hydro (w/NoahMP) output

In mccreigh/rwrfhydro: R tools for the WRF Hydro Model

Description

CalcNoahmpWatBudg calculates water budget components from WRF-Hydro (w/NoahMP) model output.

Usage

CalcNoahmpWatBudg(ldasoutDf, rtoutDf = NULL, gwoutDf = NULL,
  sfcrt = FALSE, soildeps = c(100, 300, 600, 1000), basarea = NULL)

Arguments

ldasoutDf	The LDASOUT dataframe (required)
rtoutDf	The RTOUT dataframe, if overland or subsurface routing was turned on (default=NULL)
gwoutDf	The GW_OUT dataframe, if groundwater model was turned on (default=NULL)
sfcrt	A flag whether surface overland flow routing was active. All other routing options are determined based on the input dataframes, as needed (e.g., if gwoutDf is provided, it is assumed that the groundwater model was active). (default=FALSE)
soildeps	A list of soil layer depths in mm (top to bottom, default=c(100, 300, 600, 1000))
basarea	The basin area in square km (necessary only if gwoutDf is provided)

Details

CalcNoahmpWatBudg reads dataframes derived from WRF-Hydro output (i.e., using GetMultiNcdf) and calculates water budget partitioning (e.g., surface runoff, evaporation, groundwater). Assumes WRF-Hydro output dataframes have already been masked to the desired basin. See GetMultiNcdf documentation for examples of how to do this. NOTE: Currently only works for model runs using NoahMP as the LSM.

REQUIRED variables (these terms must be in your input dataframes):

• LDASOUT: ACCPRCP, ACCECAN, ACCETRAN, ACCEDIR, SFCRNOFF, UGDRNOFF, SOIL_M (all layers), SNEQV, CANICE, CANLIQ

• RTOUT (optional, use if overland or subsurface routing were activated): QSTRMVOLRT, SFCHEADSUBRT, QBDRY

• GWOUT (optional, use if groundwater bucket model was activated): q_cms, POSIXct

OUTPUT water budget terms (may vary depending on model configuration):

• LSM_PRCP: Total precipitation (mm)

• LSM_ECAN: Total canopy evaporation (mm)

• LSM_ETRAN: Total transpiration (mm)

• LSM_EDIR: Total surface evaporation (mm)

• LSM_DELSWE: Change in snowpack snow water equivalent (mm)

• LSM_DELCANWAT: Change in canopy water storage (liquid + ice) (mm)

• LSM_SFCRNOFF: Surface runoff from LSM (for an LSM-only run) (mm)

• LSM_UGDRNOFF: Subsurface runoff from LSM (for an LSM-only run) (mm)

• LSM_DELSOILM: Change in total soil moisture storage (mm)

• HYD_QSTRMVOL: Total runoff into channel from land (routing model only) (mm)

• HYD_DELSFCHEAD: Change in surface storage (routing model only) (mm)

• HYD_QBDRY: Total flow outside of domain (routing model only) (mm)

• HYD_GWOUT: Total groundwater outflow (routing model only) (mm)

• WB_SFCRNOFF: Total surface runoff used in the water budget calculation (either LSM_SFCRNOFF or HYD_QSTRMVOL) (mm)

• WB_GWOUT: Total groundwater outflow used in the water budget calculation (either LSM_UGDRNOFF or HYD_GWOUT) (mm)

• WB_DELGWSTOR: Change in groundwater storage (adjusted by surface runoff when no surface routing active) (mm)

• ERROR: Remainder in water budget (mm)

• RUN_FRAC: Runoff fraction, runoff/precipitation

• EVAP_FRAC: Evaporative fraction, evapotranspiration/precipitation

• STOR_FRAC: Change in storage fraction, storagechange/precipitation

Value

A new dataframe containing the water budget components in mm.

See Also

Other modelEvaluation: CalcModPerfMulti, CalcModPerf, CalcNoahmpFluxes, CalcStatCategorical, CalcStatCont

Examples

## Take a NoahMP LDASOUT dataframe for a model run of Fourmile Creek with no routing
## options turned on and return a water budget summary.

## Not run: 
wb.nort.fc <- CalcNoahmpWatBudg(modLDASOUT1d.nort.fc)
wb.nort.fc
##
## Take NoahMP LDASOUT, HYDRO model RTOUT, and GW_outflow dataframes for a model
## run of Fourmile Creek with subsurface, overland, and groundwater routing options
## turned on and return a water budget summary. The default soil depths were used
## and the basin is 63.1 km2. NOTE: We MUST specify with the sfcrt flag that overland
## flow routing was turned on. Otherwise the LSM surface runoff term is used.

wb.allrt.fc <- CalcNoahmpWatBudg(modLDASOUT1d.allrt.fc, 
                                 modRTOUT1h.allrt.fc, 
                                 modGWOUT1h.allrt.fc, 
                                 sfcrt=TRUE, basarea=63.1)
wb.allrt.fc

## End(Not run)

mccreigh/rwrfhydro documentation built on Feb. 28, 2021, 1:53 p.m.