Navigation path

LEft

17-May-2022   
Path: RAMI-IV : EXPERIMENTS : ACTUAL CANOPIES : JARVSELJA WINTER BIRCHSTAND
EXPERIMENTSRESULTSMODELSPARTICIPANTS Up Print version Decrease text Increase text Home

Järvselja Birch Stand (Winter): HET15_JBS_WIN

This page provides descriptions of the architectural, spectral and illumination related properties of a 49 year old Betula pendula stand located near 58° 16′ 38.67″ N 27° 20′ 30.38″ E. The stand was inventoried in summer 2007 by Andres Kuusk, Joel Kuusk, Mait Lang, Tõnu Lükk, Matti Mõttus, Tiit Nilson, Miina Rautiainen, and Alo Eenmäe of the Tartu Observatory, in Tõravere, Estonia as well as the Estonian University of Life Sciences, Tartu, Estonia. Potential RAMI participants thus are to treat the information presented on this page as actual 'inventory data', that is, they should identify/extract those parameters and characteristics that are required as input to their canopy reflectance models. In some cases this may mean that simplifications have to be made to the available information, or, that parts of the available information can not be exploited with a given radiative transfer model. Whatever the case may be, all potential RAMI participants should mimic the standard practices that they use when matching actual field measurements to the required set(s) of input parameters of their model(s). If this means that you need more information than provided, please do not hesitate in contacting us. Last but not least, for those 3D models capable of maintaining architectural fidelity down to the individual shoot and branch level a series of ASCII (text) files containing the Cartesian coordinates of various geometric primitives (triangles, spheres and cylinders) and their transformations will be given. This should facilitate the reconstruction of the birch canopy architecture as it is described on this page.

In order to facilitate the generation of the Järvselja birch stand (Winter) the information on this page has been subdivided into four different categories. For each one of these categories the relevant descriptions will be contained within a uniquely coloured text frame and can be accessed by clicking on one of the four links below:

architectural
characteristics
spectral
characteristics
illumination
characteristics
measurements
characteristics

In case of difficulties or missing data on this page please do not hesitate in contacting us so that the problems may be resolved as fast as possible.


Architectural information up

1) General canopy characteristics

The Järvselja birch forest inventory was carried out over a 100×100 m² area placing the origin of the coordinate system at its south-western end. In order to include also the tree crowns of the inventoried tree locations within the RAMI birch Stand representation it was necessary to expand the scene area slightly beyond one hectare. Maintaining the origin of the tree location coordinate system thus resulted in some negative x,y values in the table below. Overall architectural characteristics of the scene are thus as follows:

Scene dimensions:
(ΔX × ΔY × ΔZ)
105.5115 × 106.1535 × 30.5130 [m × m × m]
(Xmin, Ymin, Zmin) −2.6507, −2.9759, 0.0 [m, m, m]
(Xmax, Ymax, Zmax) 102.8608, 103.1776, 30.5130 [m, m, m]
Number of trees in scene 1029 (465 BEPE, 205 TICO, 196 ALGL, POTR 78, PIAB 39, FREX 30, ACPL 16)
Leaf Area Index of scene* 0.0346
Fractional scene coverage** 0.2510
*The LAI of the pine trees is computed using half the total area of the needles in a shoot.
**The fractional cover is defined as 1 - direct transmission at zero solar zenith angle.

2) Foliage structure

The table below provides the structural characteristics of the Norway spruce shoots that are part of the RAMI Järvselja winter birch-stand scene. RT models capable of representing the detailed architecture of individual foliage elements may want to use the information provided in the ASCII (text) files accessible from the last row in each table below.

Tree species: Picea abies
foliage shape:
Max. foliage length o 14.2 cm
Max. radial extent of foliage+ 3.56 cm
one-sided foliage areax 83.2 cm²
twig length 3 - 6 cm#
twig diameter 0.3 cm
structural description file (ASCII) file
o The foliage length is measured from the bottom of the twig to the tip of the upper needles, i.e., it is the maximum vertical distance on the image displayed in the first row of this table
+ The maximum radial foliage extension is defined as the maximum distance from the axis connecting the bottom of the twig to the tip of the elementary foliage unit, i.e., along the vertical in the pictures displayed in the first row of this table).
x In the case of the Picea Abies shoots this one-sided foliage area value is half the total needle surface area (i.e., the value obtained when summing up the silhouette areas of all individual needles in the shoot). Here needles are assumed to be elongated spheres. If individual needles are represented as cylinders (with discs as endcaps) then the total needle area of the shoot will be different and the number of shoots per pine tree should be adjusted accordingly.
# A detailed description of the structural characteristics of the Picea abies shoot structure can be found in the header of the structural description file accessible via the link at the bottom of this table column.

3) Tree structure

The Järvselja Birch forest is generated on the basis of 18 individual tree representations. These are distributed among the following tree species: ACPL=maple (Acer Platenoides), BEPE1-4=birch (Betula Pendula), ALGL1-4=alder (Almus Glutinosa), TICO1-5=linden (Tilio Cordata), POTR1-2=poplar (Populus Tremuloides), PIAB=Norway spruce (Picea Abies), and FREX=ash (Fraxinus Exelsior). Only the Norway spruce trees contain foliage. The two tables below provides an overview of some structural characteristics of these 18 tree representations. NOTE: all footnotes are listed at the end of the second table. For those RT models capable of representing the 3D architecture of a given tree through a series of geometric primitives, the last lines of each table contain links to data files with detailed specifications of the foliage and wood structural properties of the Järvselja Birch forest (Winter) trees.

TABLE 1:
tree identifierACPLBEPE1BEPE2BEPE3BEPE4ALGL1ALGL2ALGL3ALGL4
tree height [m]15.3719.8625.4927.9930.5119.3522.5825.7627.09
crown radiusx
mean [m]:
maximum [m]:
picture
0.87
2.07
picture
0.71
1.94
picture
0.78
2.18
picture
1.09
3.10
picture
1.21
3.52
picture
0.53
1.87
picture
0.81
2.31
picture
1.10
2.71
picture
1.51
3.54
total foliage area of tree [m²]0.00.00.00.00.00.00.00.00.0
trunk DBH+ [cm]11.511.517.526.031.015.021.028.035.0
total wood area of tree [m²]14.969917.915537.247263.594667.116534.187147.25150.27973.082
vertical profile of wood areao [m²]graph
data
graph
data
graph
data
graph
data
graph
data
graph
data
graph
data
graph
data
graph
data
tree shape image
wood structure (ASCII file)woodstem
branch
stem
branch
stem
branch
stem
branch
woodwoodwoodwood

TABLE 2:
tree identifierFREXTICO1TICO2TICO3TICO4TICO5PIABPOTR1POTR2
tree height [m]13.725.9111.2714.4118.3420.7010.9025.2730.49
Foliage normal distribution: zenith angle=NANANANANANAgraph
data
NANA
Foliage normal distribution: azimuth angle=NANANANANANAgraph
data
NANA
height to live/green crown [m]NANANANANANA1.82NANA
crown radiusx
mean [m]:
maximum [m]:
picture
0.85
1.99
picture
1.07
2.48
picture
0.66
1.86
picture
0.83
2.27
picture
0.89
2.44
picture
0.96
2.57
picture
0.50
1.60
picture
1.33
2.60
picture
1.18
3.11
vertical profile of crown radius* [m]NANANANANANAgraph
data
NANA
half-total foliage area of tree# [m²]0.00.00.00.00.00.09.94130.00.0
vertical profile of leaf area o [m²]NANANANANANAgraph
data
NANA
trunk DBH+ [cm]9.57912151710.51728.0
total wood area of tree [m²]11.2675.83048.94016.52223.81741.62889.052141.3079103.829
vertical profile of wood areao [m²]graph
data
graph
data
graph
data
graph
data
graph
data
graph
data
graph
data
graph
data
graph
data
tree shape image
foliage structure (ASCII file)NANANANANANAfileNANA
wood structure (ASCII file)woodwoodwoodwoodwoodwoodwoodstem
branch
stem
branch
= For shoots the zenith angle of the foliage normal is defined as the angle between the vertical and the normal of the inner twig of the shoot (for a shoot axis aligned along the z-axis the normal was arbitrarily chosen to lie along the y-axis). Rather than spanning the full range of possible zenith angles (i.e., from 0 to 180 degree) as could be expected for non-flat asymmetric objects, it was chosen to follow the convention of foliage normals pointing only into the upper hemisphere. This is because RAMI participants, that make use of this foliage normal distribution information, will in all likelihood have models where scatterers are represented as flat (disc or equilateral triangle shaped) objects. However, should your model require a description of the foliage normal zenith angle distribution up to 180 degrees then please do not hesitate in contacting us and we will provide this information to you. For both the zenith and azimuth angle distributions the 'graph' link shows an image of the normalised foliage normal distribution versus zenith (or azimuth) angle of the foliage normal. The 'data' files for the zenith and azimuth angle distribution have three columns indicating 1) the upper value of the zenith (or azimuth) angle in a given bin, 2) the fraction of foliage area having a normal in this zenith (or azimuth) angle range, and 3) the fraction of wood area having a normal that falls in this zenith (or azimuth) angle range. Bin angle widths were chosen to be 5 degrees and 10 degrees for zenith and azimuth angles, respectively.
x The crown radius of actual trees is varying with azimuth angle. This can be seen in the various pictures showing a perspective-free nadir view of a given tree located at x=0,y=0 (concentric circles indicate the distance from the origin in steps of 0.25m). The mean and maximum values were computed from the triangle objects making up the 3D trees depicted in the picture in the the fourth-last row of each table column.
* The graphs show the maximum radial distance of foliage elements in a given height interval plotted against the upper height level of that height interval. The data files have five columns: lower-height-of-bin-in-units-of-meters     upper-height-of-bin-in-units-of-meters     minimum_radial-distance_of_foliage-in-units-of-m     maximum_radial-distance_of_foliage-in-units-of-m.     mean_radial-distance_of_foliage-in-units-of-m
# This value corresponds to the one-sided leaf area for flat leaves. For Norway spruce trees it corresponds to the sum of the (maximum) silhouettes of all the individual needles in the tree (i.e., half the total needle area per tree).
o The data files have 3 columns: lower-height-of-bin-in-units-of-meters     upper-height-of-bin-in-units-of-meters     area-of-wood-or-foliage-in-units-of-m2.
+ This is the nominal value derived from the inventory data for a tree of this height. The actual value of the tree representations provided in the ASCII files at the bottom of this table might be slightly different.

4) Stand structure

The Järvselja Birch forest (winter) is composed of 1120 individual trees. The following two tables indicates how these trees are distributed among the above tree classes and specifies their respective x,y locations of the tree centers of each tree class in the forest stand. The last row of each table contains an ASCII file with tree rotation and translation information for those RT models capable of ingesting the detailed 3D architecture of the tree models specified in the previous section.

TABLE 1:
tree identifierACPLBEPE1BEPE2BEPE3BEPE4ALGL1ALGL2ALGL3ALGL4
tree number per class 16115246921250973811
x,y coordinates of tree centers [m,m] datadatadatadatadatadatadatadatadata
tree rotations and translation (ASCII file)x datadatadatadatadatadatadatadatadata

TABLE 2:
tree identifierFREXTICO1TICO2TICO3TICO4TICO5PIABPOTR1POTR2
tree number per class 301743636220396117
x,y coordinates of tree centers [m,m] datadatadatadatadatadatadatadatadata
tree rotations and translation (ASCII file)x datadatadatadatadatadatadatadatadata
x These files contain pseudo code to rotate individual trees around their z axis and translate them from the origin to the x,y locations specified in the data files of the previous row of this table. Positive rotation angles in these files indicate that when looking down from the positive Z axis towards the origin of the coordinate system a counterclockwise rotation will result in moving the positive x axis towards the positive y axis. The angle of rotation is in the 7th column of these data files (starting from the count from 1).


Tree species and locations for the Järvselja birch stand (winter). The origin of the coordinate system is in the lower left hand side corner of the image.

RAMI participants with 3D RT models capable of representing objects using geometric primitives can download a single compressed ZIP archive with all the tree architectural ASCII information that is listed in the above tables by clicking HERE. Note: The size of the compressed archive is about 22 megabytes. It contains 44 ASCII files and can be unzipped using 'WINZIP' on windows or 'unzip' on linux/unix operating systems. Beware that the inflated archive will take up 184.9 Megabytes of storage.

spectral canopy characteristics up

Only the foliage and woody components in the Järvselja birch stand (Winter) scene feature LAMBERTIAN scattering properties. The background properties on the other hand are NON-LAMBERTIAN. To capture the directional variability of the hemispherical directional reflectance factor (HDRF) of snow, the RPV model has been fitted to a series of actual goniometer observations. The tables below contains the magnitudes of the reflectance and transmission characteristics of the various canopy components for nineteen different spectral bands as well as the RPV parameters describing the anisotropy of the background HDRF. The experimental identifier for the Järvselja Birch Stand (Winter) scene is given by HET15_JBS_WIN_B**_54 where B** relates to the spectral bands (B01, B02, …, B19). An ASCII (text) file that resumes all of the information in this table can be found here.

Trees: Lambertian scattering laws

spectral identifierB01B02B03B04B05B06B07B08B09B10
ACPL wood reflectance o0.124740.138060.157650.166890.173900.197470.206270.209900.239190.25201
BEPE stem reflectanceo0.355090.379720.402570.410540.414960.432510.436340.438050.463450.47211
BEPE branch reflectanceo0.071840.077380.089100.094540.098430.108110.100830.098800.161660.21918
ALGL wood reflectanceo0.124740.138060.157650.166890.173900.197470.206270.209900.239190.25201
FREX wood reflectanceo0.124740.138060.157650.166890.173900.197470.206270.209900.239190.25201
TICO wood reflectanceo0.124740.138060.157650.166890.173900.197470.206270.209900.239190.25201
PIAB foliage reflectance0.026380.029120.069620.078290.066680.042830.033950.032470.075020.12841
PIAB foliage transmittance0.010330.014680.044110.053130.045250.027470.020060.018020.052660.10298
PIAB wood reflectanceo0.116620.123530.135680.141250.143750.152660.154790.157910.185380.19673
POTR stem reflectanceo0.168550.161940.181240.188970.192330.209850.214510.218960.246310.27118
POTR branch reflectanceo0.111110.094380.102650.107460.110300.117380.118320.120150.140010.15722
 
spectral identifieroB11B12B13B14B15B16B17B18B19
ACPL wood reflectanceo0.258910.282830.289870.309010.363400.375820.382930.432591.0
BEPE stem reflectanceo0.479370.490290.495740.509270.538980.546020.548970.566291.0
BEPE branch reflectanceo0.254360.350290.366680.399450.464380.474120.478720.501681.0
ALGL wood reflectanceo0.258910.282830.289870.309010.363400.375820.382930.432591.0
FREX wood reflectanceo0.258910.282830.289870.309010.363400.375820.382930.432591.0
TICO wood reflectance0.258910.282830.289870.309010.363400.375820.382930.432591.0
PIAB foliage reflectance0.173230.392670.426210.446740.458550.458520.457720.434220.5
PIAB foliage transmittanceo0.146050.321750.343880.357210.369930.371300.371430.365700.5
PIAB wood reflectanceo0.202630.221770.227450.244260.298820.311540.318580.370081.0
POTR stem reflectanceo0.289830.361930.374680.399130.441650.446740.449700.463831.0
POTR branch reflectanceo0.168760.211050.218080.233840.266480.271520.273760.287231.0

Background: Anisotropic scattering

band IDB01B02B03B04B05B06B07B08B09B10
RPV ρ00.8741540.8527920.8496580.8483040.8478620.8451800.8370760.8361750.8307650.827047
RPV k1.0090201.0029301.0038001.0044501.0049801.0070001.0056701.0054101.0048001.003840
RPV Θ0.0561820.0562460.0587700.0603190.0612630.0639160.0640640.0640170.0648870.064917
RPV ρc0.8004550.7490300.7338450.7254980.7217170.7212580.7096990.7106240.7071160.702685
RPV data#filefilefilefilefilefilefilefilefilefile
 
band IDB11B12B13B14B15B16B17B18B19
RPV ρ00.8237840.8153200.8148550.7953300.7737740.7562810.7614870.6568671.0
RPV k1.0032201.0023001.0034100.9998770.9996220.9959971.0006700.9820861.0
RPV Θ0.0652370.0666150.0673500.0706380.0723730.0729840.0767450.0891450.0
RPV ρc0.6966060.6861880.6978750.6763580.7034110.7133580.7306270.6671591.0
RPV data#filefilefilefilefilefilefilefilefile
o The transmittance for woody elements (stem and branches) is equal to zero.
# These files are the BRF output generated by the RPV code when fed with the above input parameters. Note that the solar zenith angle in these files takes account of the direction of the incident light. It thus lies in the range [90-180] and is computed as 180 - SZA, where SZA is the solar zenith angle defined below. A relative azimuth of 0 (180) degree relates to forward (backward) scattering conditions. For more information on RPV please look here. Participants who wish to fit their own anisotropic background model to the RPV-simulated BRF data should inform the RAMI coordinators of this via the report files.

illumination characteristics up

The illumination conditions for the Järvselja Birch stand (winter) relate to the 1st April at GMT 9:45. More specifically the illumination features both direct and isotropic diffuse components. Direct solar light is characterised by a solar zenith angle (SZA) of 54.0 degree and a solar azimuth angle of 291.3 degree. The table below indicates the ratio of isotropically diffuse to total incident radiation for the nineteen different spectral bands:

spectral identifierB01B02B03B04B05B06B07B08B09B10
diffuse/total solar flux ratiox0.26600.18920.15680.14210.12910.09310.08280.07950.07670.0756
 
spectral identifierB11B12B13B14B15B16B17B18B19
diffuse/total solar flux ratiox0.07480.07130.07000.06650.05550.05380.05270.04460.0
x The 'direct/total solar flux ratio' is thus equal to 1 - (diffuse/total solar flux ratio).

The figure below shows perspective-free view of the Järvselja birch stand with the Cartesian coordinate system and direction of the incident solar radiation (blue arrow) superimposed. Azimuth angles are counted in an anti-clockwise direction from the positive X-axis towards the positive Y-axis as indicated by the (dotted blue) arc around the origin.

measurement characteristics up

The experimental identifier <EXP> that is needed in the naming of the various measurement results files (see file naming and formatting conventions) for the Järvselja Birch Stand (Winter) scene is given by HET15_JBS_WIN_B**_54 where B** relates to the spectral bands (B01, B02, …, B19). For each one of these spectral bands a series of radiative measurements have to be performed. In addition a lidar experiment and a fisheye experiment are proposed.

The following are the prescribed measurements for the Järvselja birch stand (Winter):


Prior to the performing of any RT model simulations, please refer to the 'definitions' pages for detailed instructions regarding the angular sign conventions for BRF simulations, as well as other RT model technicalities. Also read the relevant file naming and formatting conventions that must be adhered to by all participants. In addition, RAMI-IV offers participants the possibility to test the compliance of their model-generated results files with these file-naming and formatting convention, prior to their submission via ftp: To do so follow the on-line format checker link that appears in the top navigation bar during the active submission period.

Right navigation