CalcNoahmpWatBudg calculates water budget components from WRF-Hydro
(w/NoahMP) model output.
The LDASOUT dataframe (required)
The RTOUT dataframe, if overland or subsurface routing was turned on (default=NULL)
The GW_OUT dataframe, if groundwater model was turned on (default=NULL)
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)
A list of soil layer depths in mm (top to bottom, default=c(100, 300, 600, 1000))
The basin area in square km (necessary only if gwoutDf is provided)
CalcNoahmpWatBudg reads dataframes derived from WRF-Hydro output
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.
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
A new dataframe containing the water budget components in mm.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.