The MAESTRA Model

This website maintained by Belinda Medlyn: bmedlyn@bio.mq.edu.au
Last updated April 2007
Hosted by Biological Science at Macquarie University

CONTENTS:

Model Description Download Files On-Line Manual Bibliography

Join Mailing List Most Recent Changes Input Checklist

DESCRIPTION OF THE MAESTRA MODEL

Some history:

MAESTRO/MAESTRA is a model of forest canopy radiation absorption and photosynthesis. The model has a long history, going back to the work of Norman & Jarvis in the 1970's and 80's. Ying Ping Wang improved and tested the model for his PhD thesis, and it was published in Wang & Jarvis (1990) Agric For Meteorol 51:257-280. A lot of other people have worked on the model over the years, and as of 1997 there were several versions of the model in existence, most of which were complicated and difficult to understand or modify. In 1997 I obtained a version of the code from Ying-Ping Wang and set about revising it, with the objectives of (i) modularising the code to make the program easier to understand and modify; and (ii) incorporating standard formulations of physiological sub-modules (Farquhar - von Caemmerer, Ball-Berry). This standardised code is available from this website. 

Model description:

• The forest canopy is represented in the model as an array of tree crowns, whose positions and dimensions are specified.
• Calculations are done for just one crown, the 'target crown’.
• The distribution of leaf area within the target crown is specified, as is the leaf angle distribution.
• The main meteorological driving variables are incident radiation, temperature and humidity.
• The target crown is divided into up to120 grid points, and the radiation penetrating to each grid point is calculated for three wavebands (PAR, near infra-red and longwave). Direct, diffuse, and scattered radiation are considered separately.
• The absorbed PAR at each grid point drives the photosynthesis and transpiration routines.
• Photosynthesis is calculated from the Farquhar-von Caemmerer model.
• Stomatal conductance can be calculated from either of the Jarvis, Ball-Berry, or Ball-Berry-Leuning models.
• Transpiration is calculated by applying the Penman-Monteith equation at each gridpoint, and summing over the gridpoints.

The model is implemented in DOS-based Fortran, with text input and output files. At present there are no plans to develop a Windows version.

DOWNLOAD FILES

Please Note! The MAESTRA program is provided freely here. We expect you to use it in good faith. Please contact us if you would like to use the program in your work. Also, we take no responsibility for errors encountered in the model or in outputs generated by the model.

Also, if you publish anything using Maestra, please let me know so I can add it to the bibliography. Thankyou.

Documentation: The on-line manual gives full details of how to run the program. Please see the bibliography for articles with information regarding model operation. 

Executable: There are hourly ( EXE24.ZIP) and half-hourly ( EXE48.ZIP) versions of the program, updated April 07. Each of these files (1.2 MB, zipped) contain:

• MAESTRA.EXE is the executable program. Just type 'MAESTRA' to run.(NB currently only runs on PC).
• MAESHR.EXE is an alternative version of the program which only runs for one hour but gives a very detailed printout of what happens in that hour. Useful for checking the program!
• Sample input files: STR.DAT, PHY.DAT, TREES.DAT, CONFILE.DAT, and MET.DAT.
• Sample output files: DAYFLX.DAT, HRFLUX.DAT, LAYFLX.DAT, PTFLX.DAT, TESTFLX.DAT, DAYTOT.DAT, HRTOT.DAT, MAESERR.DAT, RESP.DAT, RESPHR.DAT
• Additional programs: SUMTREES.EXE, RESP.EXE, MAESTEST.EXE


Code: The file MAESTRA.ZIP (updated April 07; 88 KB, zipped) contains:

• The FORTRAN files MAESTRA.FOR and MAESHR.FOR which are the two alternative main programs, plus the files of code INOUT.FOR,GETMET.FOR, PHYSIOL.FOR, RADN.FOR and UTILS.FOR. These Fortran files contain all the code for the program. I know that MAESTRA can be compiled with Microsoft Fortran on a PC but have not tested it on any other systems.
• The header files MAESTCOM, METCOM - these contain program constants and are referred to in the code.
• The FORTRAN files SUMTREES.FOR, RESP.FOR and MAESTEST.FOR for the additional programs.

MAILING LIST

Would you like to be kept informed of changes to the MAESTRA program? If you write me a note I will add your name to my mailing list. This mailing list will be advised whenever an updated version of MAESTRA becomes available from these pages.

INPUT CHECKLIST

The little red dot indicates that a particular input is necessary. Other inputs are optional.

Plot details:

• Latitude and longitude.
• Stocking.
• Bearing. The bearing from north of the x-axis of the plot.
• Slope. The slope, in degrees, in x and y directions.

  Canopy structure:

• Positions of all the individual crowns. The x and y co-ordinates of each crown can be specified, or calculated from the plot size and number of trees.
• Crown dimensions. The radii in the x and y directions and the heights of the trunk and green canopy can be specified for each crown, or it can be assumed that all crowns have the same dimensions. If the simulation is to run for a long period of time, these values should be specified for a series of dates, to allow the crowns to grow.
• Leaf area. The leaf area for each crown (or an average value for all crowns) should be specified for the whole period of the simulation.
• Stem diameters.
• Canopy shape. Conical, ellipsoidal, or paraboloidal.
• Leaf inclination angle distrutibution. The distribution of leaf angles, or the average leaf angle.
• Proportion of leaf area in each leaf age class.
• Distribution of leaf area within the crown, in vertical and/or horizontal directions. This can be specified for up to three leaf age classes.
• Clumping factor: the ratio of shoot projected area : leaf projected area.


Meteorological data: on an hourly or daily basis

• Incident radiation. Either PAR or total radiation.
• Fraction of incident radiation which is direct-beam.
• Temperature. Either hourly temperature or daily minimum and maximum.
• Relative humidity or vapour pressure deficit.
• Wind speed above the canopy.
• Atmospheric CO2 concentration. Average value sufficient.
• Atmospheric pressure. Average value sufficient.


Photosynthesis Parameters:

• Jmax and Vcmax. These can be specified for different canopy layers and different foliage age classes, also on different dates. They can also be specified as functions of leaf nitrogen content, in which case the slope and intercept of the relationship are required.
• Leaf N content. This can be specified for different canopy layers and different foliage age classes, also on different dates.
• Theta. Curvature of light-response curve of Jmax.
• Parameters of temperature-response curves of Vcmax and Jmax.


Stomatal Conductance Parameters:

• Jarvis model parameters. GSREF, GSJI0, GSJD0, GSJA, GSJB, GSJTREF, GSJTMAX, GSJT0
• Ball-Berry model parameters. Slope and intercept of Ball-Berry model.
• Ball-Berry-Leuning model parameters. Slope, intercept, and D0 in Ball-Berry-Leuning model.
• Leaf width.


Respiration Parameters:

• Leaf dark respiration rate.This can be specified for different canopy layers and different foliage age classes, also on different dates. It can also be specified as a function of leaf nitrogen content, in which case the slope and intercept of the relationship are required.
• Q10 for leaf respiration.
• Fraction by which leaf dark respiration is reduced in the light.
• Growth and maintenance respiration coefficients for respiration of other tree components. Stem, branches, coarse roots, fine roots.
• Allometric coefficients for biomass as a function of stem height and diameter. These are used to calculate the amount of respiring biomass.


Reflectance and Transmittance Parameters:

• Soil reflectance. For three wavebands.
• Leaf reflectance and transmittance. For three wavebands. These can be specified separately for different canopy layers and foliage age classes.

VALIDATION DATA

The outputs from MAESTRA are the hourly values of absorbed radiation, photosynthesis and transpiration of the target tree. Useful validation data are intercepted PAR and CO2 and H2O fluxes.

CHANGES TO 31/03/07

Soil moisture effect on stomatal conductance. Soil moisture can now be input as either soil water content or potential. Three different functions are now provided for the effect of soil moisture on stomatal conductance.

Cylinder crown shape. The absorbed radiation was incorrect if you used the cylinder crown shape. This bug has been corrected.

Other Bugs. The calculation of azimuth angle in the southern hemisphere was corrected. The equation for incident thermal radiation was improved. A small bug in the interpolation of tree parameters was corrected.

Other compilers. .I'm attempting to fix the code so that it will be able to be compiled with other compilers (other than Compaq Visual Fortran). This is work in progress - I've made some small changes but more will likely come.

CHANGES TO 22/01/04

Q10 of foliage and wood respiration. Can now be specified by date.

Bugs. A couple of small bugs in the routine that calculates pathlength (DIST) were corrected. There were also some inconsistencies if you had different numbers of foliage age classes for the physiology and leaf area density distributions. 

CHANGES TO 10/09/01

Met data file. The start of the actual met data in the file met.dat must be preceded by the line
DATA STARTS
(note: caps essential; no tabs or extra blank spaces). This ensures the file is read from the correct place, regardless of system.

Soil water content and its effect on stomatal conductance. Soil water content may be specified in the met.dat file by adding a column 'SW'. This will be converted into a soil moisture deficit using the two additional parameters SWMAX and SWMIN which must be given in the namelist ENVIRON in met.dat. The soil moisture deficit is given by SMD = (SWMAX - SW) / (SWMAX - SWMIN). Thus SW should be in the same units as SWMAX and SWMIN.
The effect of soil moisture on stomatal condutance is calculated as f(SMD) = 1 - SMD1*exp(SMD2*SMD) where SMD1, SMD2 are parameters. This function is implemented for all three stomatal conductance models. The parameters SMD1, SMD2 should be added to the relevant namelists in phy.dat.
In the code you will note that work is also under way to add back in the water balance routines used by Barton & Massheder - however this has not yet been completed.

Quantum yield of electron transport. Previously this was specified as AJQ in namelist JMAXPARS in phy.dat. Alternatively, you can now specify different values according to date, layer, and age class, as you can for JMAX, VCMAX, etc. This is done by adding namelists
&AJQCON    NODATES, NOLAYERS, NOAGES \
&AJQ            DATES, VALUES \
to phy.dat.

Alternative formulations of the temperature dependences of Km and Gamma*.  Use the temperature dependences found by Bernacchi et al (2001) Plant Cell Environ. 24:253-260 by adding the parameter IECO = 0 in namelist JMAXPARS in file phy.dat. The default IECO = 1 is to use the temperature dependences opted for by the ECOCRAFT group.

Estimating incident radiation. The program can use the formula of Bristow & Campbell (1984) Ag For Met 31: 159-166 to calculate incident PAR from air temperature. To use this option it is necessary to put the namelist
& BRISTO    DELTAT    \
in the met.dat file. DELTAT should have 12 values representing the mean monthly daily temperature amplitude for each month of the year.

Stem maintenance respiration on an area basis. To specify stem maintenance respiration on an area basis rather than a mass basis, change the namelist WRESP in phy.dat to have the following parameters:
RMA - maintenance respiration per unit stem surface area (instead of RM, the rate per unit mass)
STEMFORM - used to calculate total stem surface area. Stem SA = STEMFORM * PI * HT * DBH^2
Parameters RTEMP, Q10W, EFFY stay the same.

Outputs. Mean canopy temperature (TCAN) is output on an hourly basis. There is also a column for water evaporation but this is non-operational at the moment.

Bugs. In the version of 14/3/99, a cosine was added to the calculation of WNS1 in TRANSD. This change should not have been made and has been corrected in the current version.

Thermal radiation. The formula used to estimate incident thermal radiation was changed to that of Monteith & Unsworth (1990) Principles of Environmental Physics V2 p53 - this formula takes into account higher radiation from cloudy skies.

Understorey. A program has been written to calculate radiation absorption and photosynthesis by a grass understorey. This program is not available on this web site but may be obtained by writing to Belinda Medlyn.

CHANGES TO 12/1/00

Fixed millennium bug. Dates should still be specified as ‘DD/MM/YY’. YY = 50 - 99, years 1950-1999 assumed; YY = 00 - 49, years 2000 – 2049 assumed. Someone else can change the program come 2050 …

Powerstation Version: Changes made so that program compiles in Microsoft Fortran PowerStation. IMPORTANT:  in input files, arrays involving strings (e.g. dates) must come at the end of the namelist, after all numbers – this means a change to the input file.

Added incident total short-wave radiation to hourly input met data options. Column ‘RAD’ in W m-2.

Added new VPD function for Jarvis stomatal conductance.
    F(VPD) = min(1,1/(VK1*VPD^VK2))
with parameters VK1 and VK2 in namelist JARGS in phy.dat. To use, set j = 3 in namelist MODELGS in confile.dat.

Bugs Fixed: (in order of importance)
(1) The photosynthesis subroutine PHOTOSYN was returning gross photosynthesis, not net, for the Ball-Berry-Leuning case.
(2) A correction to the subroutine DIST in radn.for to handle the case where the ray is completely outwith the crown.
(3) Handle correctly the cases where Jmax or Vcmax <= 0, where photosynthesis is below the light compensation point, and to calculate Gamma* when T < 0.
(4) If physiology is only specified for one age class, but LAD distributions are specified for several, the physiology data must be copied into arrays for the other age classes.
(5) In the calculation of the beam fraction of incident radiation, a more accurate estimate of the zenith angle is now used.

New Parameters Added:
(1) The initial slope of the light-response curve of electron transport is now the parameter AJQ (namelist: JMAXPARS in phy.dat). Default value: ALPHAQ = 0.425 mol mol-1 (defined in maestcom).
(2) Added parameters to force Jmax and Vcmax to go to zero at low temperatures. They begin to decline linearly to zero at T = TVJUP °C and reach zero at T = TVJDN °C. TVJUP, TVJDN are in namelist VCMAXPARS in phy.dat. Default values: -100.