Bigleaf—An R package for the calculation of physical and physiological ecosystem properties from eddy covariance data

We present the R package bigleaf (version 0.6.5), an open source toolset for the derivation of meteorological, aerodynamic, and physiological ecosystem properties from eddy covariance (EC) flux observations and concurrent meteorological measurements. A 'big-leaf' framework, in which vegetation is represented as a single, uniform layer, is employed to infer bulk ecosystem characteristics top-down from the measured fluxes. Central to the package is the calculation of a bulk surface/canopy conductance (Gs/Gc) and a bulk aerodynamic conductance (Ga), with the latter including formulations for the turbulent and canopy boundary layer components. The derivation of physical land surface characteristics such as surface roughness parameters, wind profile, aerodynamic and radiometric surface temperature, surface vapor pressure deficit (VPD), potential evapotranspiration (ET), imposed and equilibrium ET, as well as vegetation-atmosphere decoupling coefficients, is described. The package further provides calculation routines for physiological ecosytem properties (stomatal slope parameters, stomatal sensitivity to VPD, bulk intercellular CO2 concentration, canopy photosynthetic capacity), energy balance characteristics (closure, biochemical energy), ancillary meteorological variables (psychrometric constant, saturation vapor pressure, air density, etc.), customary unit interconversions and data filtering. The target variables can be calculated with a different degree of complexity, depending on the amount of available site-specific information. The utilities of the package are demonstrated for three single-level (above-canopy) eddy covariance sites representing a temperate grassland, a temperate needle-leaf forest, and a Mediterranean evergreen broadleaf forest. The routines are further tested for a two-level EC site (tree and grass layer) located in a Mediterranean oak savanna. The limitations and the ecophysiological interpretation of the derived ecosystem properties are discussed and practical guidelines are given. The package provides the basis for a consistent, physically sound, and reproducible characterization of biometeorological conditions and ecosystem physiology, and is applicable to EC sites across vegetation types and climatic conditions with minimal ancillary data requirements.