Retrieving Chave's environmental index

Share:

Description

Extract the Chave et al. 2014 environmental index thanks to the coordinates of the data. The function is time-consuming at its first use as it downloads a raster in the working directory. However, as soon as the raster is downloaded once, the function then runs fast (if the working directory is not changed or if the raster is copied in the new working directory).

Usage

1
computeE(coord)

Arguments

coord

Coordinates of the site(s), a matrix/dataframe with two columns (e.g. cbind(longitude, latitude)) (see examples).

Details

The Chave's environmental index, E, has been shown to be an important covariable of the diameter-height relationship of trees. It is calculated as:

E = 1.e-3 * (0.178 * TS - 0.938 * CWD - 6.61 * PS)

where TS is temperature seasonality as defined in the Worldclim dataset (bioclimatic variable 4), CWD is the climatic water deficit (in mm/yr, see Chave et al. 2014) and PS is the precipitation seasonality as defined in the Worldclim dataset (bioclimatic variable 15).

The E index is extracted from a raster file (2.5 arc-second resolution, or ca. 5 km) downloadable at http://chave.ups-tlse.fr/pantropical_allometry.htm

Value

The function returns E, the environmental index computed thanks to the Chave et al 2014 formula.

Author(s)

Jerome CHAVE, Maxime REJOU-MECHAIN, Ariane TANGUY

References

Chave et al. (2014) Improved allometric models to estimate the above-ground biomass of tropical trees, Global Change Biology, 20 (10), 3177-3190

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
# One study site
lat <- 4.08
long <- -52.68
coord <- cbind(long, lat)
## Not run: E <- computeE(coord)

# Several study sites (here three sites)
long <- c(-52.68, -51.12, -53.11)
lat <- c(4.08, 3.98, 4.12)
coord <- cbind(long, lat)
## Not run: E <- computeE(coord)