Vegetation Model User Guide

This document describes rundeck specifications required to run different configurations of the Ent Terrestrial Biosphere Model coupled to ModelE2. This includes different initialization options for different levels of dynamics. These options include:

(1) Ent biophysics (suface fluxes) only, with Matthews (1983) and Rosenzweig & Abramopoulos (1997) prescribed seasonality of leaf area index (LAI) by vegetation type. This requires only cover initilization (VEG and CROPS) and other vegetation structure is prescribed.

(2) (AVAILABLE TO DEVELOPERS) Ent biophysics and Ent prognostic phenology (interactive seasonal plant growth). This configuration can be initialized as in (1), but it is best to initialize additionally with geographically varying files for initial LAI, canopy height, and canopy maximum LAI. Data sets for these values can be obtained from Nancy Kiang.

Sample rundecks in templates directory:

E1M20iEnt.R
E4F40iEnt.R

Components:

Ent giss_LSM

Component Options:

OPTS_Ent = ONLINE=YES PS_MODEL=FBB PFT_MODEL=ENT

Description:
OPTS_Ent: compile options for Ent physics modules in model/Ent/
   PS_MODEL (FBB): FBB=Farquhar/Ball-Berry biophysics for Ent ; this will soon be the default. If not set, the Friend & Kiang (2005) biophysics scheme is still available in the code.
   PFT_MODEL (ENT): Ent 16 plant functional types are run. If not set, the Matthews (1983) 8 biome types are still available, but the Ent 16 PFTs have current support.
OPTS_giss_LSM: compile options for land hydrology physics modules in model/giss_LSM/

INPUT_FILES required include all the ground hydrology boundary conditions required for any ModelE runs:

TOPO - Fractions of ocean, lake, ground, ice, thicknesses, topography. GISS layer format.
SOIL - ModelE soil boundary conditions file. Binary file of grid arrays, no titles.
TOP_INDEX = topographic index if for TOP_MODEL run-off. GISS layer format.
VEG - global arrays for cover fraction by each plant functional type (PFT). GISS layer format.
CROPS - global arrays of historical crop cover by year. GISS layer format.
CD_coef - global array of surface roughness length. GISS layer format.

Additional INPUT_FILES required for Ent

soil_textures - global arrays of soil textures (sand, silt, clay fractions) for top 30 cm of soil for soil biogeochemistry.
SOILCARB_global (optional) - global array of total soil carbon per m2 for intialization of soil biogeochemistry. Used if PARAMETER do_soilinit=.true.
LAImax - global arrays by PFT of geographic maximum annual canopy LAI. Required for setting plant population density. Used if PARAMETER do_phenology_activegrowth=TRUE or do_init_geo=TRUE.
LAIinit - global arrays by PFT of initial leaf area index (LAI) if running prognostic seasonality. Required if PARAMETER do_phenology_activegrowth=TRUE or do_init_geo=TRUE.
HITEent - global arrays by PFT of geographic canopy heights. Used if PARAMETER do_phenology_activegrowth=TRUE or do_init_geo=TRUE.

PARAMETERS

force_VEG - logical whether to force LAI from a file (set to .false. if running prognostic phenology)
do_soilinit - logical, initialize soil carbon from file (.true. (default) or .false)
do_phenology_activegrowth - logical, simulate seasonal growth (leaf, fine roots, sapwood) and LAI (default=.true.)
do_structuralgrowth - logical, simulate woody structure growth (stem, coarse root) and height change (default=.false.)
do_init_geo - logical, initialize from file for geographic variation in canopy heights and maximum LAI (default=.true.). These are used to calculate canopy structure (plant geometry and density). Otherwise, the Rosenzweig & Abrampoulos (1983) prescriptions of PFT canopy heights and maximum LAI are used (no geographic variation) .

Diagnostics

The diagnostics for Ent include:

Carbon fluxes (mass/area/time)
GPP - gross primary productivity, gross uptake of carbon
IPP - isoprene emission from veg
RAUTO - autotrophic respiration, respiration by vegetation
SOILRESP - soil respiration

Carbon stocks (mass/area)
CLAB - plant labile carbon, carbon reserves in the plants
SOILCPOOLSUM - soil organic carbon

LAI - leaf area index (area foliage/area ground)

References

Friend, A. D. and N. Y. Kiang (2005). "Land Surface Model Development for the GISS GCM: Effects of Improved Canopy Physiology on Simulated Climate." Journal of Climate 18(15): 2883-2902.

Rosenzweig, C. and F. Abramopoulos (1997). "Land-surface model development for the GISS GCM." Journal of Climate 10: 2040-2054.