Do we (need to) care about canopy radiation schemes in DGVMs? Caveats and potential impacts

Dynamic global vegetation models (DGVMs) are an essential part of current state-of-the-art Earth system models. In recent years, the complexity of DGVMs has increased by incorporating new important processes like, e.g., nutrient cycling and land cover dynamics, while biogeophysical processes like surface radiation have not been developed much further. Canopy radiation models are however very important for the estimation of absorption and reflected fluxes and are essential for a proper estimation of surface carbon, energy and water fluxes. The present study provides an overview of current implementations of canopy radiation schemes in a couple of state-of-the-art DGVMs and assesses their accuracy in simulating canopy absorption and reflection for a variety of different surface conditions. Systematic deviations in surface albedo and fractions of absorbed photosynthetic active radiation (faPAR) are identified and potential impacts are assessed. The results show clear deviations for both, absorbed and reflected, surface solar radiation fluxes. FaPAR is typically underestimated, which results in an underestimation of gross primary productivity (GPP) for the investigated cases. The deviation can be as large as 25% in extreme cases. Deviations in surface albedo range between −0.15 ≤ Δα ≤ 0.36, with a slight positive bias on the order of Δα ≈ 0.04. Potential radiative forcing caused by albedo deviations is estimated at −1.25 ≤ RF ≤ −0.8 (W m−2), caused by neglect of the diurnal cycle of surface albedo. The present study is the first one that provides an assessment of canopy RT schemes in different currently used DGVMs together with an assessment of the potential impact of the identified deviations. The paper illustrates that there is a general need to improve the canopy radiation schemes in DGVMs and provides different perspectives for their improvement

Recovery of forest canopy parameters by inversion of multispectral LiDAR data

We describe the use of Bayesian inference techniques, notably Markov chain Monte Carlo (MCMC) and reversible jump MCMC (RJMCMC) methods, to recover forest structural and biochemical parameters from multispectral LiDAR (Light Detection and Ranging) data. We use a variable dimension, multi-layered model to represent a forest canopy or tree, and discuss the recovery of structure and depth profiles that relate to photochemical properties. We first demonstrate how simple vegetation indices such as the Normalized Differential Vegetation Index (NDVI), which relates to canopy biomass and light absorption, and Photochemical Reflectance Index (PRI) which is a measure of vegetation light use efficiency, can be measured from multispectral data. We further describe and demonstrate our layered approach on single wavelength real data, and on simulated multispectral data derived from real, rather than simulated, data sets. This evaluation shows successful recovery of a subset of parameters, as the complete recovery problem is ill-posed with the available data. We conclude that the approach has promise, and suggest future developments to address the current difficulties in parameter inversion

Latitudinal gradient of spruce forest understory and tundra phenology in Alaska as observed from satellite and ground-based data

The latitudinal gradient of the start of the growing season (SOS) and the end of the growing season (EOS) were quantified in Alaska (61°N to 71°N) using satellite-based and ground-based datasets. The Alaskan evergreen needleleaf forests are sparse and the understory vegetation has a substantial impact on the satellite signal. We evaluated SOS and EOS of understory and tundra vegetation using time-lapse camera images. From the comparison of three SOS algorithms for determining SOS from two satellite datasets (SPOT-VEGETATION and Terra-MODIS), we found that the satellite-based SOS timing was consistent with the leaf emergence of the forest understory and tundra vegetation. The ensemble average of SOS over all satellite algorithms can be used as a measure of spring leaf emergence for understory and tundra vegetation. In contrast, the relationship between the ground-based and satellite-based EOSs was not as strong as that of SOS both for boreal forest and tundra sites because of the large biases between those two EOSs (19 to 26 days). The satellite-based EOS was more relevant to snowfall events than the senescence of understory or tundra. The plant canopy radiative transfer simulation suggested that 84–86% of the NDVI seasonal amplitude could be a reasonable threshold for the EOS determination. The latitudinal gradients of SOS and EOS evaluated by the satellite and ground data were consistent and the satellite-derived SOS and EOS were 3.5 to 5.7 days degree− 1 and − 2.3 to − 2.7 days degree− 1, which corresponded to the spring (May) temperature sensitivity of − 2.5 to − 3.9 days °C− 1 in SOS and the autumn (August and September) temperature sensitivity of 3.0 to 4.6 days °C− 1 in EOS. This demonstrates the possible impact of phenology in spruce forest understory and tundra ecosystems in response to climate change in the warming Artic and sub-Arctic regions

The Laegeren site: an augmented forest laboratory combining 3-D reconstruction and radiative transfer models for trait-based assessment of functional diversity

Given the increased pressure on forests and their diversity in the context of global change, new ways of monitoring diversity are needed. Remote sensing has the potential to inform essential biodiversity variables on the global scale, but validation of data and products, particularly in remote areas, is difficult. We show how radiative transfer (RT) models, parameterized with a detailed 3-D forest reconstruction based on laser scanning, can be used to upscale leaf-level information to canopy scale. The simulation approach is compared with actual remote sensing data, showing very good agreement in both the spectral and spatial domains. In addition, we compute a set of physiological and morphological traits from airborne imaging spectroscopy and laser scanning data and show how these traits can be used to estimate the functional richness of a forest at regional scale. The presented RT modeling framework has the potential to prototype and validate future spaceborne observation concepts aimed at informing variables of biodiversity, while the trait-based mapping of diversity could augment in situ networks of diversity, providing effective spatiotemporal gap filling for a comprehensive assessment of changes to diversity

An Optical Sensor Network for Vegetation Phenology Monitoring and Satellite Data Calibration

We present a network of sites across Fennoscandia for optical sampling of vegetation properties relevant for phenology monitoring and satellite data calibration. The network currently consists of five sites, distributed along an N-S gradient through Sweden and Finland. Two sites are located in coniferous forests, one in a deciduous forest, and two on peatland. The instrumentation consists of dual-beam sensors measuring incoming and reflected red, green, NIR, and PAR fluxes at 10-min intervals, year-round. The sensors are mounted on separate masts or in flux towers in order to capture radiation reflected from within the flux footprint of current eddy covariance measurements. Our computations and model simulations demonstrate the validity of using off-nadir sampling, and we show the results from the first year of measurement. NDVI is computed and compared to that of the MODIS instrument on-board Aqua and Terra satellite platforms. PAR fluxes are partitioned into reflected and absorbed components for the ground and canopy. The measurements demonstrate that the instrumentation provides detailed information about the vegetation phenology and variations in reflectance due to snow cover variations and vegetation development. Valuable information about PAR absorption of ground and canopy is obtained that may be linked to vegetation productivity

Validation practices for satellite based earth observation data across communities

Assessing the inherent uncertainties in satellite data products is a challenging task. Different technical approaches have been developed in the Earth Observation (EO) communities to address the validation problem which results in a large variety of methods as well as terminology. This paper reviews state-of-the-art methods of satellite validation and documents their similarities and differences. First the overall validation objectives and terminologies are specified, followed by a generic mathematical formulation of the validation problem. Metrics currently used as well as more advanced EO validation approaches are introduced thereafter. An outlook on the applicability and requirements of current EO validation approaches and targets is given

Validating canopy clumping retrieval methods using hemispherical photography in a simulated Eucalypt forest

The so-called clumping factor (Ω) quantifies deviation from a random 3D distribution of material in a vegetation canopy and therefore characterises the spatial distribution of gaps within a canopy. Ω is essential to convert effective Plant or Leaf Area Index into actual LAI or PAI, which has previously been shown to have a significant impact on biophysical parameter retrieval using optical remote sensing techniques in forests, woodlands, and savannas. Here, a simulation framework was applied to assess the performance of existing in situ clumping retrieval methods in a 3D virtual forest canopy, which has a high degree of architectural realism. The virtual canopy was reconstructed using empirical data from a Box Ironbark Eucalypt forest in Eastern Australia. Hemispherical photography (HP) was assessed due to its ubiquity for indirect LAI and structure retrieval. Angular clumping retrieval method performance was evaluated using a range of structural configurations based on varying stem distribution and LAI. The CLX clumping retrieval method (Leblanc et al., 2005) with a segment size of 15° was the best performing clumping method, matching the reference values to within 0.05 Ω on average near zenith. Clumping error increased linearly with zenith angle to > 0.3 Ω (equivalent to a 30% PAI error) at 75° for all structural configurations. At larger zenith angles, PAI errors were found to be around 25–30% on average when derived from the 55–60° zenith angle. Therefore, careful consideration of zenith angle range utilised from HP is recommended. We suggest that plot or site clumping factors should be accompanied by the zenith angle used to derive them from gap size and gap size distribution methods. Furthermore, larger errors and biases were found for HPs captured within 1 m of unrepresentative large tree stems, so these situations should be avoided in practice if possible

The fourth phase of the radiative transfer model intercomparison (RAMI) exercise : Actual canopy scenarios and conformity testing

The RAdiative transfer Model Intercomparison (RAMI) activity focuses on the benchmarking of canopy radiative transfer (RT) models. For the current fourth phase of RAMI, six highly realistic virtual plant environments were constructed on the basis of intensive field data collected from (both deciduous and coniferous) forest stands as well as test sites in Europe and South Africa. Twelve RT modelling groups provided simulations of canopy scale (directional and hemispherically integrated) radiative quantities, as well as a series of binary hemispherical photographs acquired from different locations within the virtual canopies. The simulation results showed much greater variance than those recently analysed for the abstract canopy scenarios of RAMI-IV. Canopy complexity is among the most likely drivers behind operator induced errors that gave rise to the discrepancies. Conformity testing was introduced to separate the simulation results into acceptable and non-acceptable contributions. More specifically, a shared risk approach is used to evaluate the compliance of RI model simulations on the basis of reference data generated with the weighted ensemble averaging technique from ISO-13528. However, using concepts from legal metrology, the uncertainty of this reference solution will be shown to prevent a confident assessment of model performance with respect to the selected tolerance intervals. As an alternative, guarded risk decision rules will be presented to account explicitly for the uncertainty associated with the reference and candidate methods. Both guarded acceptance and guarded rejection approaches are used to make confident statements about the acceptance and/or rejection of RT model simulations with respect to the predefined tolerance intervals. (C) 2015 The Authors. Published by Elsevier Inc.Peer reviewe