energy.closure {bigleaf} | R Documentation |
Calculates the degree of the energy balance non-closure for the entire time span based on the ratio of two sums (energy balance ratio), and ordinary least squares (OLS).
energy.closure(
data,
Rn = "Rn",
G = NULL,
S = NULL,
LE = "LE",
H = "H",
instantaneous = FALSE,
missing.G.as.NA = FALSE,
missing.S.as.NA = FALSE
)
data |
Data.frame or matrix containing all required variables. |
Rn |
Net radiation (W m-2) |
G |
Ground heat flux (W m-2); optional |
S |
Sum of all storage fluxes (W m-2); optional |
LE |
Latent heat flux (W m-2) |
H |
Sensible heat flux (W m-2) |
instantaneous |
should the energy balance be calculated at the time step
of the observations ( |
missing.G.as.NA |
if |
missing.S.as.NA |
if |
The energy balance ratio (EBR) is calculated as:
EBR = sum(LE + H)/sum(Rn - G - S)
the sum is taken for all time steps with complete observations (i.e. where all energy balance terms are available).
a named vector containing:
n |
number of complete (all energy balance terms available) observations |
intercept |
intercept of the OLS regression |
slope |
slope of the OLS regression |
r_squared |
r^2 of the OLS regression |
EBR |
energy balance ratio |
if instantaneous = TRUE
, only EBR
is returned.
Wilson K., et al. 2002: Energy balance closure at FLUXNET sites. Agricultural and Forest Meteorology 113, 223-243.
## characterize energy balance closure for DE-Tha in June 2014
energy.closure(DE_Tha_Jun_2014,instantaneous=FALSE)
## look at half-hourly closure
EBR_inst <- energy.closure(DE_Tha_Jun_2014,instantaneous=TRUE)
summary(EBR_inst)