Estimating Forest fAPAR from Multispectral Landsat-8 Data Using the Invertible Forest Reflectance Model INFORM

The estimation of the Fraction of Absorbed Photosynthetically Active Radiation in forests (forest fAPAR) from multi-spectral Landsat-8 data is investigated in this paper using a physically based radiative transfer model (Invertible Forest Reflectance Model, INFORM) combined with an inversion strategy based on artificial neural nets (ANN). To derive the forest fAPAR for the Dabie mountain test site in China in 30 m spatial resolution (size approximately 3000 km2), a database of forest canopy spectral reflectances was simulated with INFORM taking into account structural variables such as leaf area index (LAI), crown coverage and stem density as well as leaf composition. To establish the relationship between forest fAPAR and the reflectance modeled by INFORM, a logarithmic relationship between LAI and fAPAR was used previously established using on-site field measurements. On this basis, predictive models between Landsat-8 reflectance and fAPAR were established using an artificial neural network. After calibrating INFORM for the test site, forty-two forest stands were used to validate the performance of the method. The results show that spectral signatures modeled by INFORM correspond reasonably well with the forest canopy reflectance spectra derived from Landsat data. Deviations increase with increasing angle between surface normal of the hilly terrain and OPEN ACCESS Remote Sens. 2015, 7 7426 sun incidence. The comparison of estimated and measured fAPAR (R2 = 0.47, RMSE = 0.11) demonstrates that INFORM can be inverted using neural nets to provide acceptable estimates of forest fAPAR. The accuracy of the predictions increased significantly when excluding pixels located in very steep terrain. This demonstrates that the applied topographic correction was not sufficiently accurate and should be improved for making optimum use of radiative transfer models such as INFORM.


Introduction
The estimation of the Fraction of Absorbed Photosynthetically Active Radiation (fAPAR) provides key information for estimating vegetation net primary productivity (NPP) [1].Consequently, fAPAR has been identified by the Global Climate Observing System (GCOS) as one of the Essential Climate Variables (ECVs).Forest ecosystems play a critical role in the global carbon cycle but the accurate retrieval of forest fAPAR is hindered by the complex interaction of electromagnetic radiation with forest canopies [2].Additional problems arise in hilly terrain.
Most studies use empirical models to estimate forest fAPAR [3].The functional relationships (linear or nonlinear) between satellite-derived vegetation indices and field-measured fAPAR established by regression analysis are easy to understand and computationally efficient to operate over large scales.They also constitute an efficient means to generate reference information from field samples useful for testing alternative approaches [4,5].Empirical models, however, do not exploit the complex transfer of solar radiation within the forest stands and can only be applied to the conditions where they were developed.
A potentially more robust approach consists in using radiative transfer models.Physically-based models can handle many different types of conditions and account for the effect of biochemical, structural or observational variables on canopy reflectance [6].After model construction and calibration, it is possible to invert canopy reflectance models (CRM) to retrieve variables of interest from the measured reflectance.Canopy reflectance models act thus as the interface between the remotely sensed images and forest variables of interest.A prerequisite for using CRMs is that the satellite observations have to be corrected for atmospheric and terrain effects so that the resulting image corresponds to the Bottom of Atmosphere (BoA) reflectance [7][8][9].
Since the 1970s, a variety of CRMs have been developed to simulate forest canopy reflectance.In early CRMs such as SUITS (a radiative transfer model named by the developer' name) [10] and SAIL(Scattering by Arbitrarily Inclined Leaves) [11], vegetation canopies have been modeled assuming a homogenous canopy layer composed of infinite small scattering elements (the leaves).The models do not mimic the influence of clumpiness, open spots, and shadowing effects, but assume a homogenous distribution of scattering and absorbing particles (e.g., leaves) in the canopy.This over-simplifies forest canopies, where leaves are generally not homogeneously distributed, and where shadows play an important role in the scene reflectance [12,13].The simplification of canopy layer in the mentioned models also neglects the effects of understory which is known to have a strong influence on canopy reflectance [14].Sometimes the sensitivity of some bands onto the overstory leaf area index (LAI) even disappears when the understory reaches an LAI higher than one [15].
To deal with shadowing effects resulting from tree crowns, geometrical-optical models such as the Li-Strahler model [16] describe the forest as a composite of opaque geometrical shapes.These models, however, mostly ignore the non-zero transmissivity of tree crowns.
In recent years, the spectral invariants theory (p-theory) [17] was applied in forests and a PARAS (a simple parameterization model using photon recollision probability, a spectral invariant which describes the structure of a canopy, for calculating forest stand reflectance.)model family was developed [18][19][20][21].PARAS models require only a limited set of input parameters and the canopy structure is described using only one parameter, the photon recollision probability p.Therefore, the PARAS models are also called semi-physical forest reflectance models, neglecting crown level clumping and the associated shading patterns on the background, as well as multiple interactions of photons between the tree layer and the understory layer [20].
Among existing CRMs for forests (for an overview see [22]), INFORM [23,24] provides a simple while still realistic description of the interaction of the electromagnetic radiation with forest canopies.This computationally simple model takes into account the effects of shadowing and crown transmittance as well as the effects of leaf geometry and leaf optical properties on the crown transmittance.The impact of leaf geometry and leaf optical properties were neglected by the forest-light interaction model (FLIM) [25] on which INFORM builds.
Whatever CRM is used, it is worth noting that fAPAR usually cannot be derived directly by inverting CRMs as fAPAR is usually not a build-in variable of CRMs.The common way to retrieve fAPAR from the modeled reflectance is utilizing the modeled spectral reflectance or its derivatives (e.g., VIs) to establish a linear or non-linear relationship with fAPAR [14,26,27].This leads to so-called semi-empirical models.A close relationship between canopy reflectance and fAPAR can be expected because both are functionally linked to leaf area index (LAI) [28].For example, based on the continuous observations of fAPAR and LAI during the growing period, linear or logarithmic relationships were identified between the two biophysical variables [29][30][31][32].Equally well known is also the fact that LAI has a close relationship with the reflectance modeled from CRM, which is the foundation for estimates of LAI by inverting the CRM.Therefore, the relationship between LAI and fAPAR, as well as between LAI and canopy reflectance, presents a possible link between the modeled canopy reflectance and fAPAR.This link can be exploited using various techniques such as look-up-tables (LUT) and numerical optimization procedures [33][34][35].Compared with numerical optimization techniques and LUT, a well calibrated artificial neural network (ANN) can be applied to estimate variables with higher computational efficiency and more accurate approximations of complex non-linear functions [36].
In this context, this paper will investigate the estimation of forest fAPAR from multispectral Landsat-8 data using the INFORM model.Important sub-tasks include: (a) to calibrate the INFORM parameters for adapting the CRM to local conditions, (b) to evaluate its performance in modeling the forest canopy reflectance by sensitivity analysis and comparison with observed spectral signatures, and (c) to estimate and evaluate forest fAPAR by inverting INFORM based on ANN.
As the study is conducted in a hilly environment, special attention will be given to possible terrain effects.

Study Area and Field Measurements
The Dabie study area (115°24′-116°3′E, 31°13′-31°44′N) is located on the northern slope of the Dabie mountain with steep topography and complex terrain.Dabie is located in southern China and covers an area of about 3814 km 2 , including a forest area of 2940 km 2 .Twenty-one relatively homogenous broad-leaved forest stands with poplar, chestnut and beech as main species were identified, located between 125 m to 1210 m above mean sea level (Figure 1).Within each of these stands, two quadratic plots of 30 m length were established leading to a total number of forty-two reference plots.Terrain effects represent one of the most important perturbations in remote sensing over mountainous terrain [37].To study terrain effects, several reference plots with low cosine of the incidence angle (angle i between the direction of the sun and the local surface normal) [38] were chosen.Information regarding elevation, slope, aspect and cos i is reported in Table 1 for each of the 42 reference plots.As shown in Figure 2, according to the location of the reference plots and the sun, two groups are distinguished: (1) plots with 'high' cos i > 0.51 (n = 33; shown in black), and (2) plots with "low" cos i ≤ 0.51 (n = 9; shown in bold).
The central location of each reference plot was determined with an accuracy of about ±5 m using a differential GPS device (Trimble Pro XR).At each plot, measurements of forest structural variables were conducted in October 2012 including canopy fAPAR, LAI, stem density (SD), crown diameter (CD) and stand height (H) (Table 2).
Table 1.Information regarding elevation, slope and aspect of the 42 individual field plots.The variable cosi "i" is the angle between the direction of the sun and the local surface normal.The 42 reference plots are split into two groups: (1) reference plots with high cosi (e.g., cosi > 0.51) in black, and (2) reference plots with low cosi (e.g., cosi ≤ 0.51) in bold.The latter reference plots represent non-favorable observations conditions (orientation of the surface) and will therefore be analyzed separately.On clear days, fAPAR measurements were taken using a Tracing Radiation and Architecture of Canopies (TRAC) instrument which was developed at Canada Centre for Remote Sensing [39].A fish-eye lens instrument was used under overcast condition (Regent Instruments Inc., Tokyo, Japan).Compared to TRAC showing fAPAR directly, the fish-eye lens provided images requiring further processing by CAN-EYE software [40].LAI was estimated using a Li-Cor LAI-2200 Plant Canopy Analyzer under overcast sky conditions (Lincoln Inc., Nebraska) [41].As shown in Figure 3, the red dots are the locations of LAI-2200 measurements and TRAC measurements, and the box respects one plot of 30 m length, which corresponds to one Landsat pixel.LAI and fAPAR measurements were taken at eight regularly spaced points within each plot, from which the average was calculated.Stem density (SD) was obtained by counting the number of trees in a plot.Stand height (H) was calculated from the mean height of all trees.The height of each tree (from the ground to the top) was estimated from angular measurements.Crown diameter (CD) was the average of the crown diameter of the trees.Stem density, stand height and crown diameter were all measured from the trees whose diameter at breast height (DBH) were greater than 5 cm within a plot.In a forest stand, the average soil reflectance was determined from three measurements using a FieldSpec 4 Hi-Res spectroradiometer (Analytical Spectra Devices., Inc., Boulder, CO, USA).All canopy variables were measured between 10:00 and 16:00 to minimize sun angle effects.

Image Data and Image Pre-Processing
For modelling and mapping LAI, a Landsat-8 image acquired on 12 October 2013 was used.The Landsat-8 image was radiometrically corrected by utilizing the Radiometric Calibration module and the FLAASH Atmospheric Correction module integrated in ENVI 5.0 with the spectral response function [9].
As most of the existing CRMs are not suitable for areas with steep topography and complex terrains, the Landsat-8 image was topographically corrected with the C-correction model [42,43] (Figure 4).Topographic corrections with this simple model remove major illumination induced brightness changes.However, residual effects are to be expected [38,44] and were evident in our data set (Figure 5).No further attempts were made to remove the remaining terrain effects, as such data quality is to be expected in typical remote sensing application studies.

Method
The estimation of forest fAPAR from multispectral remote sensing data and the INFORM canopy reflectance model consisted of five steps (Figure 6): (a) establishing an empirical relationships between field measured LAI and in situ fAPAR, (b) modeling visible and near-infrared bi-directional spectral reflectance of forest canopies with variable structure, leaf bio-chemistry and canopy background reflectance using INFORM, (c) using INFORM to generate a synthetic fAPAR dataset through application of the empirical relationship between LAI and fAPAR to the LAI values used as input to the synthetic reflectance data base, (d) determining the predictive relationship between spectral reflectance and fAPAR of the synthetic data base based on back-propagation (BP) neural networks, and (e) applying the trained ANN to remote sensing data and analyzing the simulated results.Figure 6.General procedure for estimating forest fAPAR from Landsat data using a forest reflectance model.

INFORM
INFORM was used to simulate the canopy reflectance of forest stands between 400 and 2500 nm [23,24].The model is based on FLIM [25], combined with SAILH (SAIL model with hot) [11,45], and PROSPECT (a radiative transfer model based on Alle's generalized "plate model" that represents the optical properties of plant leaves from 400 nm to 2500 nm) [46].
In INFORM, a forest is considered as a discontinuous canopy layer with tree crowns and gaps [25].The average crown transmittance in sun direction and in observation direction (Ts and To) is calculated from SAIL.The required leaf transmittance and reflectance (τ and ρ) are calculated from PROSPECT.The bi-directional reflectance of the forest canopy in a given spectral band is expressed as: where RC is the crown reflectance at infinite crown depth and RG is the background reflectance.C is the "crown factor": where Ts and To are the average crown transmittance in the sun direction and in observation direction, respectively; co and cs are the ground coverage by crowns in observation direction and by shadow in sun direction, respectively.
G is the "ground factor": where Fcd, Fcs, Fod and Fos are the ground fractions, representing crowns with shadowed background, crowns with sunlit background, shadowed open space and sunlit open space respectively.The four fractions are calculated as follows: where rho represents the positive correlation between co and cs as: where g is a geometrical factor depending on the illumination and observing geometry, θo and θs are observation and solar zenith angles, and Ψ is the relative azimuth between sun and observer.H is the tree height.CD is the tree crown diameter.
To derive C and G, several variables, such as co and cs, are defined.
where SD is the stem density and k is the average crown horizontal area depending on CD.
Regarding the average crown transmittance, the SAILH model is used to estimate Ts and To with the leaf area index (LAI), average leaf inclination angle (ALA), leaf reflectance (ρleaf), leaf transmittance (τleaf), observation zenith angle (θo), solar zenith angle (θs), the relative azimuth angle between sun and observer (Ψ), the fraction of diffuse radiation (skyl) and the hot spot (hot) [45] In addition, RG and RC are also derived from the SAILH model.
where LAIU and ALAU are LAI and ALA of the understory, LAIinf is the LAI at infinite crown depth and ρsoil is the soil reflectance (from field measurements).The leaf transmittance and reflectance (τleaf and ρleaf) in Equation ( 12) are calculated using the PROSPECT model [46] with the leaf structure parameter (N), chlorophyll content (a+b) (Cab), cellulose and lignin content (Cm) and equivalent water thickness (Cw) as input parameters.
As shown in Table 3, INFORM is parameterized based on twelve synthetic variables.Consequently, the reflectance of forest stands canopy can be expressed as: Except the forest stand variables (SD, CD, LAI and H, Table 2), the remaining input parameters (geometrical information and the leaf-level parameters) were fixed (Table 4).The geometrical information (observation and sun zenith angles; azimuth angle) was taken from the corresponding Landsat-8 orbit characteristics, location and date.The leaf-level parameters including leaf chlorophyll content (a+b), dry matter content and equivalent water thickness were taken from LOPEX'93 experiment [47] representing poplar, chestnut and beech leaves.During the processing of the field-measured data, it was found that stem density (SD) and crown diameter (CD) were closely related.This permitted the constraining of the INFORM parameterization by using their fitted relation: SD = 8388CD −1.486 (R 2 of 0.68, n = 42) (Figure 7).

Training the Artificial Neural Network
To invert the CRM, Artificial Neural Networks (ANN) were employed.ANN "learn" the relation between a set of input variables (here the INFORM simulated forest reflectances) and the target variable (here: fAPAR) [48].In the present work, a simple three layer feed-forward backpropagation network was used with tan-sigmoidal transfer functions in hidden and output layers.The ANN was trained using a training database generated from INFORM.To assess forest fAPAR, a previously established relation between LAI and fAPAR was used (see next section and Figure 8).For successful training, the distribution of inputs and outputs was important as well as a compact network architecture and efficient training algorithms.

Generating the Training Database
The training database consisted of the input variables (i.e., the modeled spectra from INFORM) and the corresponding target variables (i.e., the canopy structural and bio-chemical variables used to simulate the spectra and with fAPAR linked to LAI).
To ensure that the training data was consistent with the real Landsat-8 data, the continuous reflectance data modeled from CRM were spectrally resampled to multispectral wavebands using the appropriate spectral response functions, following: where ρ is the resampled reflectance, ρs is the modeled reflectance, and B(λ) is the spectral response function.To permit using alternative satellite sensors, only the first four Landsat channels were simulated: blue (band 1), green (band 2), red (band 3) and near infrared (band 4).
Several studies found that fAPAR had a significant correlation with LAI at the canopy-level.Canopy LAI (LAIc) can be expressed as the product of single tree LAI (LAIS) and crown closure (Co).The crown closure Co was defined as the ground coverage by crowns in nadir direction and was calculated as a function of stem density (SD) and crown diameter (CD) (Equations ( 6) and ( 8)).Consequently, the canopy LAI (LAIc) was computed as: Using field data (n = 42), the logarithmic function relating fAPAR with LAI was determined (Figure 8): On this basis, a fAPAR dataset was generated corresponding with the simulated spectral reflectances.Both the values of fAPAR and spectral reflectance fall between 0 and 1, so no data normalization was applied.Our empirical relation corresponds well with similar models reported in the literature [49].

Establishing the Network Architecture and Training
The backpropagation ANN was realized within MATLAB software (www.mathworks.com).The chosen neural network architecture consists of one input layer with four linear neurons for input of the three visible bands and the near-infrared reflectance of Landsat-8 plus one hidden layer with 20 neurons and one output layer with one neuron (fAPAR).The study was restricted to the four mentioned bands as these bands are included in most sensors, such as Quickbird, SPOT and two small Chinese optical satellites (HJ-1A and HJ-1B).Hence, the developed methodology can potentially be applied to map fAPAR from other sensors as well.Both the hidden layer and output layer use the sigmoid transfer function.The train epochs were set as 100.
Once the ANN was established, the training data were divided into three subsets [36].Half of the simulated datasets were randomly selected for training, and the remaining half were divided into two parts: validation datasets and testing datasets.The validation datasets were used during the learning process of the network to prevent over-fitting.The testing datasets serve as internal evaluation of the ANN estimation performance.
The training process to calibrate ANN coefficients relied on the minimization of a misfit function by a back-propagation algorithm [50].During this process, the values of the weights and biases in the network were iteratively adjusted to minimize the error between the network outputs and the measured outputs (fAPAR).The coefficient of determination (R 2 ) and the root mean squared error (RMSE) between the network outputs and the measured outputs were calculated to measure accuracy of the model inversion.The optimal network architecture was determined after several trials to provide the best estimation performance.Once calibrated, a neural net can be considered as an approximation of the surface response of the inverse radiative transfer model by means of a simple parametric model with known coefficients.A calibrated ANN can be applied easily to real image data.In case the retrieved fAPAR value fell outside the definition range (0-1), the value of the closest bound was used.

Sensitivity Analysis of INFORM Parameters
Sensitivity analysis of model parameters is an important step in the process of model adaptation to regional conditions for later generation of synthetic reflectance data using CRMs.It is also an efficient way to investigate the effects of forest canopy variables on reflectance signatures and helps understanding how canopy structure and leaf bio-chemicals influence the model simulations.
In the sensitivity analysis, the variable parameters to be analyzed were classified into five ranks and the corresponding reflectance spectra were simulated from INFORM.During the simulation, other parameters in INFORM were kept constant (shown in Table 4).The analyzed variables were varied within the range of measured data (Table 2) providing information on the sensitivity of INFORM generated forest canopy reflectance for different stem densities, crown diameters, tree heights, dry matter contents and chlorophyll contents (a+b), respectively, in the wavelength range between 400 and 2500 nm (Figure 9).When one variable (e.g., LAI) was analyzed, other variables (e.g., SD, CD, H, Cab and Cm) were set as the average of the minimum and maximum.
Stem density (SD), crown diameter (CD) and LAI had very strong influence in all bands, in particular in the near-infrared.With increasing stem density (SD), the modeled reflectance gradually became saturated.The chlorophyll contents (a+b) and dry matter content had little influence in all bands due to very small parameter ranges.As a result, the chlorophyll contents (a+b) and dry matter content were kept constant for the remaining study (at 0.0035 and 44, respectively) (see Table 4).

Validation of the Modelled Canopy Reflectance
For the forty-two reference plots of our study area, forest reflectance spectra from 400 to 2500 nm were simulated with INFORM and resampled to Landsat-8 wavebands (visible and near-infrared bands).For the simulations, the observed forest structural and bio-chemical parameters were used.Modelled reflectances in the visible and near-infrared wavelength regions of samples with relatively high cosi (e.g., above 0.51-see Table 1) mostly compared well to the measured values (Figure 10).Most points were relatively closely distributed along the 1:1 line (e.g., within ±0.01 reflectance units in the visible channels and ±0.1 in the nIR) thus confirming other experimental results by [24].
On the contrary, for the nine samples taken in steep terrain (e.g., with low cosi-see Table 1), it was not possible to achieve a good forward simulation accuracy with INFORM.RMSE between measured and simulated reflectances were substantially higher compared to the samples with high cosi (Table 5).This finding clearly shows that the topographic correction was not well performing in very steep terrain.For this reason, it was decided to clearly distinguish between these two samples (e.g., high/low cosi) when analysing the fAPAR modeling results, as one cannot expect good CRM inversion results if the observed reflectance spectra are perturbed by terrain effects, etc.
As the improvement of existing topographic correction algorithms was out of the scope of the present study, we preferred to work with the imperfect Landsat data and instead distinguished the two groups of reference samples (see Table 1).Obviously, the imperfections related to the terrain correction would need to be resolved if similar methods are to be applied in hilly terrain.

Validation of fAPAR Inversion Results
A scatterplot (n = 33) of field measured and Landsat-8 based forest fAPAR values is shown in Figure 11 excluding the nine samples from very steep terrain.Most points are slightly below or scatter around the 1:1 line with a RMSE of 0.11.The explained variance (R 2 ) amounts to 0.47.As stems are "misinterpreted" as leaves in optical devices used for taking reference measurements, a slight overestimation of field measurements can be expected.
We noticed again a strong impact of topography (e.g., imperfect topographic correction) on modelling results.Regarding the fAPAR estimations obtained from the nine samples with low cosi, only 3% of the variance was explained with RMSE of 0.16 (Table 6).The main reason for the strong errors are related to the topography of the study area.

Estimation of Forest fAPAR in the Study Area
The trained ANN was applied to the Landsat-8 imagery of the entire study test site.No pixels were excluded for terrain orientation.As shown in Figures 12 and 13, estimated values of fAPAR are generally in a reasonable range.fAPAR values are stretched in Figure 12 along a color ramp for cartographic reasons.Although most of the forest stands in the fAPAR map appear rather homogenous, certain variations can be observed.Gradual variations in forest fAPAR are expected according to changing growing conditions.On the other hand, extreme values appearing in some regions with low cosi (e.g., high uncertainty) are possibly artefacts related to terrain effects.

Conclusions
In this paper, the forest fAPAR was estimated from Landsat-8 data through inversion of the physically-based INFORM radiative transfer model based on ANN.Three main conclusions can be drawn from our study: (1) INFORM seems moderately well suited for modeling the bi-directional spectral reflectance of forest canopies in the visible and near-infrared bands as a function of structural and bio-chemical forest characteristics.The parameters most strongly influencing the simulated reflectances of the model are LAI, crown diameter, and stem density.The modelled reflectance showed a consistent underestimation in red band and overestimation in near-infrared band.Probably, part of the differences could be reduced by using higher precision field measurement data.For example, in-situ leaf characteristics and other species-related factors were neglected in this study or unavailable.Obviously, being a relatively simple CRM, INFORM only captures the main factors of variability.More advanced CRMs need to be tested to reach a closer agreement between forward simulations and Landsat observations.(2) The forest fAPAR was modeled for the Dabie mountain test site by inverting INFORM through an artificial neural network approach.Results suggest that this method can successfully estimate forest fAPAR from multispectral images and achieve an acceptable accuracy with RMSE of 0.11 (14% of average) and R 2 of 0.47.In further studies, efforts will be made to estimate and evaluate time series of forest fAPAR based on the method.(3) The mountainous terrain posed the biggest challenge for successful retrieval of fAPAR.The employed simple topographic C-correction method created several overcorrections and artifacts negatively influencing the fAPAR retrieval.The insufficient topographic correction already became visible when running INFORM in forward mode (e.g., with field measured canopy characteristics entered into the CRM).In doing so, it was noted that the simulated canopy reflectance spectra strongly deviated from the observed (Landsat-8) spectral profiles.

Figure 1 .
Figure 1.Location of the study area in China (Left) and the field-measured sampling plots (red dots) shown in a digital elevation model (DEM) (Right).At each of the twenty-one field sampling plots, two reference plots were established.The total number of reference plots thus equals 42.

Figure 2 .
Figure 2. Location relationship between the reference plots and the sun.

Figure 3 .
Figure 3. Sampling scheme for TRAC transects and LAI-2200 measurement points.The distance between two LAI measurement points was ten meters in each cardinal direction.

Figure 4 .
Figure 4. False color composite image (bands 4-3-2) of the study area after atmospheric and terrain related preprocessing.

Figure 5 .
Figure 5. Observed impact of cosi (local terrain orientation with respect to the sun) on Landsat-8 reflectance in the nIR band.The linear regression line between near-infrared band reflectance and cosi (n = 42) is also shown.

Figure 7 .
Figure 7. Empirical relationship between crown diameter and stem density used in INFORM to link stem density (SD) to crown size (CD) (n = 42).

Figure 8 .
Figure 8. fAPAR as a function of LAI for the forest stands at three test sites (n = 42).Both, fAPAR and LAI were measured in the field.

Figure 9 .
Figure 9. INFORM generated forest reflectance spectra (0.4-2.5 µm) for different values of stem density, crown diameter, LAI, height, dry matter content and leaf chlorophyll content (a+b), respectively.The ranges of simulated structural parameter were chosen according to observed ranges.Ranges of bio-chemical parameters were chosen according to LOPEX experiment.

Figure 10 .
Figure 10.Modeled and measured forest canopy reflectance (33 plots from 'normal' terrain-high cosi) in different wavebands: (a) blue band; (b) green band; (c) red band and (d) near-infrared band of Landsat-8.The 1:1 line is shown as continuous line.The dotted lines in (a), (b) and (c), are described as y = x ± 0.01, and in (d), y = x ± 0.1.

Figure 12 .
Figure 12.Distribution of fAPAR in the study area Dabie mountain in (a) and distribution of cosi in (b).Note that non-forest pixels are shown in white.

Figure 13 .
Figure 13.Frequency distribution of the forest fAPAR map of the study area.Only pixels covered with forest are taken into account.No pixel is excluded for terrain orientation.

Table 2 .
Summary statistics for field measured fraction of absorbed photosynthetically active radiation in forests (fAPAR) and forest stand variables (n = 42).
as input parameters.

Table 3 .
Summary of main invertible forest reflectance model (INFORM) parameters.Variable meaning is given in the main text.

Table 4 .
INFORM parameters kept constant during model simulations.

Table 5 .
Accuracy of INFORM reflectance simulations with field measured canopy characteristics as inputs.The RMSE is between INFORM simulations and Landsat-8 observations.We distinguish steep terrain (e.g., low cosi) and normal terrain (e.g., high cosi).

Table 6 .
Accuracy of forest fAPAR from INFORM against field measured values.We distinguish steep terrain (e.g., low cosi) and normal terrain (e.g., high cosi).