energy_closure: Energy Balance Closure
In grahamstewart12/tidyflux:

Description Usage Arguments Details Value References Examples

Calculates the degree of the energy balance non-closure for individual records, or the entire time span based on the ratio of two sums (energy balance ratio), and ordinary least squares (OLS).

1
2
3

energy_closure(data, Rn = "Rn", G = NULL, S = NULL, LE = "LE",
  H = "H", instantaneous = TRUE, name_out = "EBR",
  missing.G.as.NA = FALSE, missing.S.as.NA = FALSE)

`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 (`TRUE`), or over the entire time period provided as input (`FALSE`).
`missing.G.as.NA`	If `TRUE`, missing G are treated as `NA`s, otherwise set to 0.
`missing.S.as.NA`	If `TRUE`, missing S are treated as `NA`s, otherwise set to 0.
`x`	A data.frame or matrix containing all required variables.
`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.

The energy balance ratio (EBR) is calculated as:

EBR = sum(LE + H)/sum(Rn - G - S)

where the sum is taken for all time steps with complete observations (i.e. all energy balance terms are available).

A named vector.

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)