An assessment of the MODIS collection 5 leaf area index product for a region of mixed coniferous forest

Canopy leaf area index (LAI), defined as the single-sided leaf area per unit ground area, is a quantitative measure of canopy foliar area. LAI is a controlling biophysical property of vegetation function, and quantifying LAI is thus vital for understanding energy, carbon and water fluxes between the land surface and the atmosphere. LAI is routinely available from Earth Observation (EO) instruments such as MODIS. However EO-derived estimates of LAI require validation before they are utilised by the ecosystem modelling community. Previous validation work on the MODIS collection 4 (c4) product suggested considerable error especially in forested biomes, and as a result significant modification of the MODIS LAI algorithm has been made for the most recent collection 5 (c5). As a result of these changes the current MODIS LAI product has not been widely validated. We present a validation of the MODIS c5 LAI product over a 121 km2 area of mixed coniferous forest in Oregon, USA, based on detailed ground measurements which we have upscaled using high resolution EO data. Our analysis suggests that c5 shows a much more realistic temporal LAI dynamic over c4 values for the site we examined. We find improved spatial consistency between the MODIS c5 LAI product and upscaled in situ measurements. However results also suggest that the c5 LAI product underestimates the upper range of upscaled in situ LAI measurements.

Preliminary leaf area index estimates from airborne small footprint full-waveform LiDAR data

This study has compared preliminary estimates of effective leaf area index (LAI) derived from fish-eye lens photographs to those estimated from airborne full-waveform small-footprint LiDAR data for a forest dataset in Australia. The full-waveform data was decomposed and optimized using a trust-region-reflective algorithm to extract denser point clouds. LAI LiDAR estimates were derived in two ways (1) from the probability of discrete pulses reaching the ground without being intercepted (point method) and (2) from raw waveform canopy height profile processing adapted to small-footprint laser altimetry (waveform method) accounting for reflectance ratio between vegetation and ground. The best results, that matched hemispherical photography estimates, were achieved for the waveform method with a study area-adjusted reflectance ratio of 0.4 (RMSE of 0.15 and 0.03 at plot and site level, respectively). The point method generally overestimated, whereas the waveform method with an arbitrary reflectance ratio of 0.5 underestimated the fish-eye lens LAI estimates.

Specific leaf area and leaf area index in developing stands of Fagus sylvatica L. and Picea abies Karst

European beech (Fagus sylvatica L.) and Norway spruce (Picea abies Karst.) are two of the most ecologically and economically important forest tree species in Europe. These two species co-occur in many locations in Europe, leading to direct competition for canopy space. Foliage characteristics of two naturally regenerated pure stands of beech and spruce with fully closed canopies were contrasted to assess the dynamic relationship between foliage adaptability to shading, stand LAI and tree growth. We found that individual leaf size is far more conservative in spruce than in beech. Individual leaf and needle area was larger at the top than at the bottom of the canopy in both species. Inverse relationship was found for specific leaf area (SLA), highest SLA values were found at lowest light availability under the canopy. There was no difference in leaf area index (LAI) between the two stands, however LAI increased from 10.8 to 14.6 m2m-2 between 2009 and 2011. Dominant trees of both species were more efficient in converting foliage mass or area to produce stem biomass, although this relationship changed with age and was species-specific. Overall, we found larger foliage plasticity in beech than in spruce in relation to light conditions, indicating larger capacity to exploit niche openings.

Exploring the spatial variation of take-all (Gaeumannomyces graminis var. tritici) for site-specific management

The fungus Gaeumannomyces graminis var. tritici (Ggt), commonly known as the take-all fungus, causes damage to roots of wheat and barley that limits crop growth and causes loss of yield. There was little knowledge on the within-field spatial variation of take-all and relations with features in the growing crop, selected soil properties and spectral information from remotely sensed imagery. Geostatistical analyses showed that take-all, chlorosis and leaf area index had similar patchy distributions. Many of the spectral bands from a hyperspectral image also had similar spatial patterns to take-all and chlorosis. Relations between take-all and mineral nitrogen, elevation and pH were generally weaker.

Effect of plant density and initial crown size on growth, development and yield in strawberry cultivars Elsanta and Bolero

The effects of density (plant spacing) and initial plant size on vegetative growth, flowering and fruiting were studied in the strawberry cultivars Elsanta and Bolero in their first and second years of cropping. The influence of these factors on light use and dry-matter partitioning was investigated. The size of planting material in 'Elsanta' and 'Bolero' slightly affected plant growth and yield, but this effect was not consistent and radiation use efficiency (RUE) and harvest index were unaltered. Plant spacing did not significantly affect the early stages of crop growth, but was important in determining growth and yield later in the season, this effect being more significant in the second year of cropping. Plant growth and yield per plant increased as plant spacing increased from 20 to 30 cm in both 'Elsanta' and 'Bolero', but the highest harvest index and yield per square metre were obtained at the closest spacing. Increased plant spacing also resulted in a greater leaf area and leaf area index. However, light was used less efficiently resulting in a lower RUE and lower harvest index (HI).

Effects of cultural system (organic and conventional) on growth, photosynthesis and yield components of sweet corn (zea mays L.) under semi-arid environment

Organic sweet maize consists of a new industrial crop product. Field experiment was conducted to determine the effects of cultural systems on growth, photosynthesis and yield components of sweet maize crop (Zea mays L. F-1 hybrid 'Midas'). A randomized complete block design was employed with four replicates per treatment (organic fertilization: cow manure (5, 10 and 20 t ha(-1)), poultry manure (5, 10 and 20 t ha(-1)) and barley mulch (5, 10 and 20 t ha(-1)), synthetic fertilizer (240 kg N ha(-1)): 21-0-0 and control). The lowest dry weight, height and leaf area index and sod organic matter were measured in the control treatment. Organic matter content was proportionate to the amount of manure applied. The control plots had the lowest yield (1593 kg ha(-1)) and the double rate cow manure plots the had,greatest one. (6104 kg ha(-1)). High correlation between sweet corn yield and organic matter was registered. Moreover, the lowest values of 1000-grain weight were obtained with control plot. The fertilizer plot gave values which were similar to the full rate cow manure treatment. The photosynthetic race of the untreated control was significantly lower than that of the other treatments. The phorosynthetic rate increased as poultry manure and barley mulch ram decreased and as cow manure increased. Furthermore the untreated control had the lowest stomatal conductance and chlorophyll content. Our results indicated that sweet corn growth and yield in the organic plots was significantly higher than those in the conventional plots.

Effect of plant density and initial crown size on growth, development and yield in strawberry cultivars Elsanta and Bolero

The effects of density (plant spacing) and initial plant size on vegetative growth, flowering and fruiting were studied in the strawberry cultivars Elsanta and Bolero in their first and second years of cropping. The influence of these factors on light use and dry-matter partitioning was investigated. The size of planting material in 'Elsanta' and 'Bolero' slightly affected plant growth and yield, but this effect was not consistent and radiation use efficiency (RUE) and harvest index were unaltered. Plant spacing did not significantly affect the early stages of crop growth, but was important in determining growth and yield later in the season, this effect being more significant in the second year of cropping. Plant growth and yield per plant increased as plant spacing increased from 20 to 30 cm in both 'Elsanta' and 'Bolero', but the highest harvest index and yield per square metre were obtained at the closest spacing. Increased plant spacing also resulted in a greater leaf area and leaf area index. However, light was used less efficiently resulting in a lower RUE and lower harvest index (HI).

Characterization of a highly biodiverse floodplain meadow using hyperspectral remote sensing within a plant functional trait framework

We assessed the potential for using optical functional types as effective markers to monitor changes in vegetation in floodplain meadows associated with changes in their local environment. Floodplain meadows are challenging ecosystems for monitoring and conservation because of their highly biodiverse nature. Our aim was to understand and explain spectral differences among key members of floodplain meadows and also characterize differences with respect to functional traits. The study was conducted on a typical floodplain meadow in UK (MG4-type, mesotrophic grassland type 4, according to British National Vegetation Classification). We compared two approaches to characterize floodplain communities using field spectroscopy. The first approach was sub-community based, in which we collected spectral signatures for species groupings indicating two distinct eco-hydrological conditions (dry and wet soil indicator species). The other approach was “species-specific”, in which we focused on the spectral reflectance of three key species found on the meadow. One herb species is a typical member of the MG4 floodplain meadow community, while the other two species, sedge and rush, represent wetland vegetation. We also monitored vegetation biophysical and functional properties as well as soil nutrients and ground water levels. We found that the vegetation classes representing meadow sub-communities could not be spectrally distinguished from each other, whereas the individual herb species was found to have a distinctly different spectral signature from the sedge and rush species. The spectral differences between these three species could be explained by their observed differences in plant biophysical parameters, as corroborated through radiative transfer model simulations. These parameters, such as leaf area index, leaf dry matter content, leaf water content, and specific leaf area, along with other functional parameters, such as maximum carboxylation capacity and leaf nitrogen content, also helped explain the species’ differences in functional dynamics. Groundwater level and soil nitrogen availability, which are important factors governing plant nutrient status, were also found to be significantly different for the herb/wetland species’ locations. The study concludes that spectrally distinguishable species, typical for a highly biodiverse site such as a floodplain meadow, could potentially be used as target species to monitor vegetation dynamics under changing environmental conditions.

Spatio-temporal surface soil heat flux estimates from satellite data: results for the AMMA experiment at the Fakara (Niger) supersite

This paper describes a method that employs Earth Observation (EO) data to calculate spatiotemporal estimates of soil heat flux, G, using a physically-based method (the Analytical Method). The method involves a harmonic analysis of land surface temperature (LST) data. It also requires an estimate of near-surface soil thermal inertia; this property depends on soil textural composition and varies as a function of soil moisture content. The EO data needed to drive the model equations, and the ground-based data required to provide verification of the method, were obtained over the Fakara domain within the African Monsoon Multidisciplinary Analysis (AMMA) program. LST estimates (3 km × 3 km, one image 15 min−1) were derived from MSG-SEVIRI data. Soil moisture estimates were obtained from ENVISAT-ASAR data, while estimates of leaf area index, LAI, (to calculate the effect of the canopy on G, largely due to radiation extinction) were obtained from SPOT-HRV images. The variation of these variables over the Fakara domain, and implications for values of G derived from them, were discussed. Results showed that this method provides reliable large-scale spatiotemporal estimates of G. Variations in G could largely be explained by the variability in the model input variables. Furthermore, it was shown that this method is relatively insensitive to model parameters related to the vegetation or soil texture. However, the strong sensitivity of thermal inertia to soil moisture content at low values of relative saturation (<0.2) means that in arid or semi-arid climates accurate estimates of surface soil moisture content are of utmost importance, if reliable estimates of G are to be obtained. This method has the potential to improve large-scale evaporation estimates, to aid land surface model prediction and to advance research that aims to explain failure in energy balance closure of meteorological field studies.

The REFLEX project: Comparing different algorithms and implementations for the inversion of a terrestrial ecosystem model against eddy covariance data

We describe a model-data fusion (MDF) inter-comparison project (REFLEX), which compared various algorithms for estimating carbon (C) model parameters consistent with both measured carbon fluxes and states and a simple C model. Participants were provided with the model and with both synthetic net ecosystem exchange (NEE) of CO2 and leaf area index (LAI) data, generated from the model with added noise, and observed NEE and LAI data from two eddy covariance sites. Participants endeavoured to estimate model parameters and states consistent with the model for all cases over the two years for which data were provided, and generate predictions for one additional year without observations. Nine participants contributed results using Metropolis algorithms, Kalman filters and a genetic algorithm. For the synthetic data case, parameter estimates compared well with the true values. The results of the analyses indicated that parameters linked directly to gross primary production (GPP) and ecosystem respiration, such as those related to foliage allocation and turnover, or temperature sensitivity of heterotrophic respiration, were best constrained and characterised. Poorly estimated parameters were those related to the allocation to and turnover of fine root/wood pools. Estimates of confidence intervals varied among algorithms, but several algorithms successfully located the true values of annual fluxes from synthetic experiments within relatively narrow 90% confidence intervals, achieving >80% success rate and mean NEE confidence intervals <110 gC m−2 year−1 for the synthetic case. Annual C flux estimates generated by participants generally agreed with gap-filling approaches using half-hourly data. The estimation of ecosystem respiration and GPP through MDF agreed well with outputs from partitioning studies using half-hourly data. Confidence limits on annual NEE increased by an average of 88% in the prediction year compared to the previous year, when data were available. Confidence intervals on annual NEE increased by 30% when observed data were used instead of synthetic data, reflecting and quantifying the addition of model error. Finally, our analyses indicated that incorporating additional constraints, using data on C pools (wood, soil and fine roots) would help to reduce uncertainties for model parameters poorly served by eddy covariance data.

A data assimilation method for using low-resolution Earth observation data in heterogeneous ecosystems

We present an approach for dealing with coarse-resolution Earth observations (EO) in terrestrial ecosystem data assimilation schemes. The use of coarse-scale observations in ecological data assimilation schemes is complicated by spatial heterogeneity and nonlinear processes in natural ecosystems. If these complications are not appropriately dealt with, then the data assimilation will produce biased results. The “disaggregation” approach that we describe in this paper combines frequent coarse-resolution observations with temporally sparse fine-resolution measurements. We demonstrate the approach using a demonstration data set based on measurements of an Arctic ecosystem. In this example, normalized difference vegetation index observations are assimilated into a “zero-order” model of leaf area index and carbon uptake. The disaggregation approach conserves key ecosystem characteristics regardless of the observation resolution and estimates the carbon uptake to within 1% of the demonstration data set “truth.” Assimilating the same data in the normal manner, but without the disaggregation approach, results in carbon uptake being underestimated by 58% at an observation resolution of 250 m. The disaggregation method allows the combination of multiresolution EO and improves in spatial resolution if observations are located on a grid that shifts from one observation time to the next. Additionally, the approach is not tied to a particular data assimilation scheme, model, or EO product and can cope with complex observation distributions, as it makes no implicit assumptions of normality.

Canopy characteristics of contrasting clones of cacao (Theobroma cacao)

Canopy characteristics (leaf area index, fractional light interception, extinction coefficient) of mature trees of ten clonally propagated cacao cultivars were measured over a period of 14 months at an experimental site in Bahia, Brazil. Differences in leaf area index between clones became more pronounced over time. When an approximately constant leaf area index was reached (after about nine months), LAI varied between clones from 2.8 to 4.5. Clonal differences in the relationship between leaf area index and fractional light interception implied differences in canopy architecture, as reflected by the range of extinction coefficients (mean values ranged from 0.63 for the clone TSH-565 to 0.82 for CC-10). The results demonstrate the potential for breeding more photosynthetically efficient cacao canopies.

The boreal forest and the polar front

The analysis presented in this paper suggests that the larger heating over the boreal forest in the spring and summer, as contrasted with weaker heating over the adjacent tundra, results in a preferred position of the polar front along the northern edge of the boreal forest. This positioning is well documented in the literature (see, for example, Bryson, 1966; Barry and Hare, 1974; Kreps and Barry, 1970). This heating results from the lower albedo of the boreal forest which is not compensated by an increase in transpiration, even with the larger leaf area index of the forest. The warmer temperatures are mixed upward by the deep boundary layer over the forest and mesoscale circulations which result from the patchiness of heating associated with the heterogeneous landscapes of the forest. Thus in contrast to previous assumptions in which the arctic front position in the summer determines the northern limit of the boreal tree line, our study suggests the boreal forest itself significantly influences the preferred position of the front. This conclusion reinforces the findings of Bonan et al. (1992) and Foley et al. (1994) on the important role of boreal forest-tundra interactions with climate.

The simulation of the opposing fluxes of latent heat and CO2 over various land-use types: coupling a gas exchange model to a mesoscale atmospheric model

A mesoscale meteorological model (FOOT3DK) is coupled with a gas exchange model to simulate surface fluxes of CO2 and H2O under field conditions. The gas exchange model consists of a C3 single leaf photosynthesis sub-model and an extended big leaf (sun/shade) sub-model that divides the canopy into sunlit and shaded fractions. Simulated CO2 fluxes of the stand-alone version of the gas exchange model correspond well to eddy-covariance measurements at a test site in a rural area in the west of Germany. The coupled FOOT3DK/gas exchange model is validated for the diurnal cycle at singular grid points, and delivers realistic fluxes with respect to their order of magnitude and to the general daily course. Compared to the Jarvis-based big leaf scheme, simulations of latent heat fluxes with a photosynthesis-based scheme for stomatal conductance are more realistic. As expected, flux averages are strongly influenced by the underlying land cover. While the simulated net ecosystem exchange is highly correlated with leaf area index, this correlation is much weaker for the latent heat flux. Photosynthetic CO2 uptake is associated with transpirational water loss via the stomata, and the resulting opposing surface fluxes of CO2 and H2O are reproduced with the model approach. Over vegetated surfaces it is shown that the coupling of a photosynthesis-based gas exchange model with the land-surface scheme of a mesoscale model results in more realistic simulated latent heat fluxes.

Linking urban water balance and energy balance models to analyse urban design options

Using a water balance modelling framework, this paper analyses the effects of urban design on the water balance, with a focus on evapotranspiration and storm water. First, two quite different urban water balance models are compared: Aquacycle which has been calibrated for a suburban catchment in Canberra, Australia, and the single-source urban evapotranspiration-interception scheme (SUES), an energy-based approach with a biophysically advanced representation of interception and evapotranspiration. A fair agreement between the two modelled estimates of evapotranspiration was significantly improved by allowing the vegetation cover (leaf area index, LAI) to vary seasonally, demonstrating the potential of SUES to quantify the links between water sensitive urban design and microclimates and the advantage of comparing the two modelling approaches. The comparison also revealed where improvements to SUES are needed, chiefly through improved estimates of vegetation cover dynamics as input to SUES, and more rigorous parameterization of the surface resistance equations using local-scale suburban flux measurements. Second, Aquacycle is used to identify the impact of an array of water sensitive urban design features on the water balance terms. This analysis confirms the potential to passively control urban microclimate by suburban design features that maximize evapotranspiration, such as vegetated roofs. The subsequent effects on daily maximum air temperatures are estimated using an atmospheric boundary layer budget. Potential energy savings of about 2% in summer cooling are estimated from this analysis. This is a clear ‘return on investment’ of using water to maintain urban greenspace, whether as parks distributed throughout an urban area or individual gardens or vegetated roofs.