Search

Vegetation indices (VIs), which describe remotely sensed vegetation properties such as photosynthetic activity and canopy structure, are widely used to study vegetation dynamics across scales. However, VI-based results can vary between indices, sensors, quality control measures, compositing algorithms, and atmospheric and sun–target–sensor geometry corrections. These variations make it difficult to draw robust conclusions about ecosystem change and highlight the need for consistent VI application and verification. In this Technical Review, we summarize the history and ecological applications of VIs and the linkages and inconsistencies between them. VIs have been used since the early 1970s and have evolved rapidly with the emergence of new satellite sensors with more spectral channels, new scientific demands and advances in spectroscopy. When choosing VIs, the spectral sensitivity and features of VIs and their suitability for target application should be considered. During data analyses, steps must be taken to minimize the impact of artefacts, VI results should be verified with in situ data when possible and conclusions should be based on multiple sets of indicators. Next-generation VIs with higher signal-to-noise ratios and fewer artefacts will be possible with new satellite missions and integration with emerging vegetation metrics such as solar-induced chlorophyll fluorescence, providing opportunities for studying terrestrial ecosystems globally.

Mapping canopy photosynthesis in both high spatial and temporal resolution is essential for carbon cycle monitoring in heterogeneous areas. However, well established satellites in sun-synchronous orbits such as Sentinel-2, Landsat and MODIS can only provide either high spatial or high temporal resolution but not both. Recently established CubeSat satellite constellations have created an opportunity to overcome this resolution trade-off. In particular, Planet Fusion allows full utilization of the CubeSat data resolution and coverage while maintaining high radiometric quality. In this study, we used the Planet Fusion surface reflectance product to calculate daily, 3-m resolution, gap-free maps of the near-infrared radiation reflected from vegetation (NIRvP). We then evaluated the performance of these NIRvP maps for estimating canopy photosynthesis by comparing with data from a flux tower network in Sacramento-San Joaquin Delta, California, USA. Overall, NIRvP maps captured temporal variations in canopy photosynthesis of individual sites, despite changes in water extent in the wetlands and frequent mowing in the crop fields. When combining data from all sites, however, we found that robust agreement between NIRvP maps and canopy photosynthesis could only be achieved when matching NIRvP maps to the flux tower footprints. In this case of matched footprints, NIRvP maps showed considerably better performance than in situ NIRvP in estimating canopy photosynthesis both for daily sum and data around the time of satellite overpass (R 2 = 0.78 vs. 0.60, for maps vs. in situ for the satellite overpass time case). This difference in performance was mostly due to the higher degree of consistency in slopes of NIRvP -canopy photosynthesis relationships across the study sites for flux tower footprint-matched maps. Our results show the importance of matching satellite observations to the flux tower footprint and demonstrate the potential of CubeSat constellation imagery to monitor canopy photosynthesis remotely at high spatio-temporal resolution.

Continuous monitoring of vegetation indices (VIs) the fraction of absorbed photosynthetically active radiation (fPAR) and leaf area index (LAI) through satellite remote sensing has advanced our understanding of biosphere–atmosphere interactions. Substantial efforts have been put into monitoring individual variables in the field, but options to concurrently monitor VIs, fPAR, and LAI in-situ have been lacking. In this paper, we present the Smart Surface Sensing System (4S), which automatically collects, transfers and processes VIs, fPAR and LAI data streams. The 4S consists of a microcomputer, controller and camera, a multi-spectral spectrometer built in with a light-emitting diode (LED) and an internet connection. Lab testing and field observations in a rice paddy site that experiences wet summer monsoon seasons confirmed the linear response of 4S to light intensities in the blue, green, red and near-infrared spectral channels, with wide ranging temperatures and humidity having only a minor impact on 4S throughout the growing season. Applied over an entire rice growing season (day of year [DOY] 120 - 248), VIs and fPAR from 4S were linearly related to corresponding VIs from a reference spectrometer (R2 = 0.98; NDVI, R2 = 0.96; EVI) and the LAI-2200 instrument (R2 = 0.76), respectively. Integration of gap fraction-based LAI from LED sensors and a green index from the micro-camera allowed tracking of the seasonality of green LAI. The continuous and diverse nature of 4S observations highlights its potential for evaluating satellite remote sensing products. We believe that 4S will be useful for the expansion of ecological sensing networks across multiple spatial and temporal scales.

Conventional process-based crop simulation models and agro-land surface models require numerous forcing variables and input parameters. The regional application of these crop simulation models is complicated by factors concerning input data requirements and parameter uncertainty. In addition, the empirical remotely sensed regional scale crop yield estimation method does not enable growth process modeling. In this study, we developed a process-based rice yield estimation model by integrating an assimilate allocation module into the satellite remote sensing-derived and biophysical process-based Breathing Earth System Simulator (BESS). Normalized accumulated gross primary productivity ( G P P n o r m - a c c u ) was used as a scaler for growth development, and the relationships between G P P n o r m - a c c u and dry matter partitioning coefficients were determined from the eddy covariance and biometric measurements at the Cheorwon Rice paddy KoFlux site. Over 95% of the variation in the dry matter allocation coefficients of rice grain could be explained by G P P n o r m - a c c u . The dynamics of dry matter distribution among different rice components were simulated, and the annual grain yields were estimated. BESS-Rice simulated GPP and dry matter partitioning dynamics, and rice yields were evaluated against in-situ measurements at three paddy rice sites registered in KoFlux. The results showed that BESS-Rice performed well in terms of rice productivity estimation, with average root mean square error (RMSE) value of 2.2 g C m−2 d−1 (29.5%) and bias of –0.5 g C m−2 d−1 (–7.1%) for daily GPP, and an average RMSE value of 534.8 kg ha−1 (7.7%) and bias of 242.1 kg ha−1 (3.5%) for the annual yield, respectively. BESS-Rice is much simpler than conventional crop models and this helps to reduce the uncertainty related to the forcing variables and input parameters and can result in improved regional yield estimation. The process-based mechanism of BESS-Rice also enables an agronomic diagnosis to be made and the potential impacts of climate change on rice productivity to be investigated.

Time series of vegetation indices (e.g. normalized difference vegetation index [NDVI]) and color indices (e.g. green chromatic coordinate [GCC]) based on radiometric measurements are now available at different spatial and temporal scales ranging from weekly satellite observations to sub-hourly in situ measurements by means of near-surface remote sensing (e.g. spectral sensors or digital cameras). In situ measurements are essential for providing validation data for satellite-derived vegetation indices. In this study we used a recently developed method to calculate NDVI from near-infrared (NIR) enabled digital cameras (NDVIC) at 17 sites (for a total of 74 year-sites) encompassing six plant functional types (PFT) from the PhenoCam network. The seasonality of NDVIC was comparable to both NDVI measured by ground spectral sensors and by the moderate resolution imaging spectroradiometer (MODIS). We calculated site- and PFT-specific scaling factors to correct NDVIC values and recommend the use of site-specific NDVI from MODIS in order to scale NDVIC. We also compared GCC extracted from red-green-blue images to NDVIC and found PFT-dependent systematic differences in their seasonalities. During senescence, NDVIC lags behind GCC in deciduous broad-leaf forests and grasslands, suggesting that GCC is more sensitive to changes in leaf color and NDVIC is more sensitive to changes in leaf area. In evergreen forests, NDVIC peaks later than GCC in spring, probably tracking the processes of shoot elongation and new needle formation. Both GCC and NDVIC can be used as validation tools for the MODIS Land Cover Dynamics Product (MCD12Q2) for deciduous broad-leaf spring phenology, whereas NDVIC is more comparable than GCC with autumn phenology derived from MODIS. For evergreen forests, we found a poor relationship between MCD12Q2 and camera-derived phenology, highlighting the need for more work to better characterize the seasonality of both canopy structure and leaf biochemistry in those ecosystems. Our results demonstrate that NDVIC is in excellent agreement with NDVI obtained from spectral measurements, and that NDVIC and GCC can complement each other in describing ecosystem phenology. Additionally, NDVIC allows the detection of structural changes in the canopy that cannot be detected by visible-wavelength imagery.

Satellite image fusion methods that improve spatial and temporal resolution have significant potential to advance understanding of ecosystem dynamics in space and time. However, systematic evaluations of image fusion methods against in situ spectral data are lacking. Here, we used a suite of in situ spectral data collected at 60 elementary sampling units (10 × 10 m) covering 15 Landsat pixel (30 × 30 m) plots and one Moderate Resolution Imaging Spectroradiometer (MODIS) pixel (250 × 250 m) throughout the entire growing season in a heterogeneous rice paddy landscape to evaluate four state-of-the-art image fusion NDVI products. They include the Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model (ESTARFM), Flexible Spatiotemporal DAta Fusion (FSDAF), SaTellite dAta IntegRation (STAIR), and the CubeSat Enabled Spatio-Temporal Enhancement Method (CESTEM); the former three blended Landsat and MODIS data, whereas the latter combined CubeSats, Landsat, and MODIS observations. All fusion products showed strong linear relationships against in situ data when combining all spatial and temporal observations (R2: 0.73 to 0.93) although there were partly negative biases (–1% to –9%). These biases resulted from forcing data to image fusion algorithms, such as Landsat (–4%) and MODIS (–7%). Performance difference between fusion methods were considerably larger for spatial than for temporal variation. Furthermore, Landsat NDVI explained only 17–22% of spatial variation against in situ spectral data, which can be translated into weak performance of image fusion products to predict spatial variability in NDVI. Image fusion products that relied on spatial interpolation showed large biases (–15% to –30%) for a vegetation plot surrounded by mixed land cover plots. Our results highlight key sources of uncertainty and will be instrumental in improving satellite image fusion methods to monitor land surface phenology in space and time.

It is assumed that the ratio of actual evapotranspiration (AET) to potential evapotranspiration (PET) is mostly controlled by the soil water content available for ET. This control is formulated using the soil moisture stress function (SSF), where the evaporative fraction (EF) or the fraction of the AET to PET (fPET) is assumed to be either a linear or a non-linear function of soil moisture. We examine the effectiveness of the soil moisture stress function to quantify soil moisture control on EF or fPET over a dryland wheat field in Victoria, Australia. Micrometeorological observations from two cropping seasons were used for the analysis. The efficacy of a root-density-weighted soil moisture estimate in predicting EF and fPET was investigated as against the commonly assumed fixed-depth root zone soil moisture. However, results indicate a strong relationship between EF and available soil water fraction (AWF) in the root zone only when solar radiation is higher than 5 MJ/m2/day. As the rooting depth increases with vegetation growth, SSF exhibits the strongest correlation with AWF for increasing soil profile depth. In the early and harvesting crop growth stages, ET is constrained mostly by surface soil moisture (0–5 cm). In the mid-growth stages, ET is strongly influenced by soil moisture in the root zone (0–60 cm). The shape of SSF, however, changes significantly between the two years (2012 and 2013). It is inferred that different temporal rainfall patterns between the years caused wheat’s different response to water stress.

Digital cover photography (DCP) has emerged as an indirect method to measure gap fraction of vegetation canopies. However, as with other photographic methods, determining camera relative exposure value (REV) and threshold for pixel classification, cause substantial uncertainties in gap fraction estimates. Here we propose a new method to improve the measurement of gap fraction under various solar zenith angles (SZAs), sky conditions, and canopy structures. This method computes gap fractions of ambiguous vegetation or sky pixels using an unsaturated raw image from DCP and a reconstructed sky image from the raw image, thus taking full advantage of the potential of raw image processing. This is combined with pre-classification of pixels that are unambiguously canopy and sky to greatly reduce light scattering effects that are likely to be present within the canopy. To test the sensitivity of the new method, we acquired images at one-hour intervals between 20 and 85° of SZAs under closed, half-closed, and open canopies with REV settings from 0 to −5. The new method showed little variation in gap fractions across the diverse SZAs in closed, half-closed, and open canopies. A perforated panel experiment, which was used to test the accuracy of the estimated gap fractions, confirmed that the new method accurately estimated gap fractions across a range of hole size, gap fractions and SZAs. We conclude that the new method opens new opportunities to estimate gap fractions accurately regardless of solar positions from open to closed canopies, and is a significant advance for accurate and precise monitoring of canopy cover and leaf area index (LAI), and for calibrating and evaluating satellite remote sensing LAI products.

Addressing the impact of vegetation architecture on shortwave radiation transfer in land surface models is important for accurate weather forecasting, carbon budget estimates, and climate predictions. This paper investigates to what extent it is possible to retrieve structural parameters of two different parameterization schemes from direct transmittance derived from digital hemispherical photography and 3D radiative transfer modeling for two study sites with different vegetation canopy architectures. Neglecting the representation of 3D canopy structure in radiative transfer schemes leads to significant errors in shortwave radiation partitioning (up to 3.5 times more direct transmittance in the 3D model). Structural parameters, referred to as whole canopy ‘clumping indices’, were obtained in order to evaluate the impact of angular variation in clumping on shortwave radiation transfer. Impacts on photosynthesis were evaluated at site level with the UKESM land surface model, JULES. A comparison between flux tower derived and modeled photosynthesis indicates that considering zenith angular variations of structural parameters in the radiative transfer scheme of the UKESM land surface model significantly improves photosynthesis prediction in light limited ecosystems (from RMSE = 2.91 μmol CO2.m − 2.s − 1 to RMSE = 1.51 μmol CO2.m − 2.s − 1, 48% smaller), typically with enhanced photosynthesis from bottom layers.

Rice is a major crop that feeds more than half of the global population. Rice paddies play a complex role in the carbon cycle by emitting CH4 to the atmosphere while either sequestering or releasing CO2 to the atmosphere. Here, we present 3.5 years of eddy covariance measurements of CO2 and CH4 fluxes over an intermittently irrigated, single crop rice paddy in South Korea in tandem with carbon stock measurements for leaves, stems, grains, and roots. The rice paddy acted as a slight sink or neutral in CO2 (−47 ± 51 g C m−2 y−1; mean ± 95% confidence interval) with 941 ± 130 g C m−2 y−1 and 894 ± 108 g C m−2 y−1 as the gross primary production (GPP) and ecosystem respiration (Reco), respectively. The mean annual net primary productivity (NPP) was 666 ± 31 g C m−2 y−1, which was partitioned to leaves (9%), stems (30%), grains (51%), and roots (10%), for a harvest index of 57%. The paddy emitted 20.6 ± 1.5 g C m−2 y−1 in the form of CH4, of which between 82 and 96% was emitted during the cultivation period (~132 days) . The CH4 flux showed a mid-season decline caused by mid-season drainage, which produced aerobic soil conditions that inhibited CH4 emissions. The proportion of estimated ebullitive flux to total CH4 flux ranged from 10 to 17% during the cultivation period. GPP, Reco, and CH4 fluxes varied by 12%, 11%, and 7%, respectively, over the 3-year period. When considering the harvest and the global warming potential of CH4, the paddy shifted from a greenhouse gas sink to a strong source (2,741 ± 468 g CO2 eq m−2 y−1). These findings will improve our understanding of carbon budgets and cycles, and help develop greenhouse gas and water management strategies for rice paddies.

a b s t r a c t Estimates of the canopy gap fraction, on which calculations of leaf area index (LAI) are based, are sensi- tive to photographic exposure in upward-facing images. In this article we describe a simple, automated method of image acquisition and processing that eliminates both subjectivity and the need for the oper- ator to consider photographic exposure in the field. A key strength of our methodology was the use of a test apparatus (perforated aluminum screen) with a precisely known gap fraction; this allowed us to separate the confounding effects of gap size and gap fraction on the optimal photographic exposure for a canopy. We took photographs of the test apparatus at different photographic exposures; we varied the gap fraction by covering a proportion of the holes in the screen, and also varied gap size by varying the distance of the camera from the screen. We acquired both raw images and JPEG images. We found that the optimal exposure of JPEG images varied with both gap size and gap fraction, not just gap fraction as previously assumed. Underexposing by one stop yielded raw data that were never clipped resulting in data loss, but that used most of the 14-bit range of the raw file. We also found that it was easily possible to standardize photographic exposure during image processing by acquiring raw images in the field; thus eliminating the variation in estimated gap fraction and LAI associated with exposure variations. This result was replicated in both fisheye images and cover images that we acquired in real canopies. We recommend the following protocol for acquiring canopy images in future studies

Correct estimates of gap fraction are essential for quantifying canopy architectural variables, such as leaf area and clumping indices, which modify land–atmosphere interactions. However, gap fraction measurements from optical sensors are contaminated by radiation that is scattered by plant elements and ground surfaces. In this study, we propose a simple one-dimensional, invertible, bidirectional transmission model to remove scattering effects from gap fraction measurements. To evaluate how well the proposed model computes scattered radiance under a variety of ecosystem conditions, we compared simulated scattered radiance by the proposed model to a more sophisticated three-dimensional model in four ecosystem types (oak-grass savanna, birch, pine, and spruce stands). The simple model showed good agreement with the three-dimensional model in the scattering factor (scattered radiation from leaves normalized by sky radiation), except for highly reflective stems such as birch. The simple model showed that the scattering factor is highest when the leaf area index (LAI) is low (1–2 m 2 m −2 ) in a non-clumped canopy, potential errors in estimating the LAI increase with an increase in LAI, and bright land surfaces (e.g., snow and bright soil) and bright stems (e.g., birch) can contribute significantly to scattering effects. By applying the simple model with LAI-2200 data collected in an oak-grass savanna woodland, we found that the scattering factor causes significant underestimation of the LAI (up to 26% for sunny conditions, 7.7% for diffuse sky conditions) and significant overestimation of the apparent clumping index (up to 14% for sunny conditions, 4.3% for diffuse sky conditions). The LAI is underestimated because of the effect of scattered radiation on gap fraction estimates, which cause overestimation of the clumping index. Even under highly diffuse sky conditions, errors in LAI estimates due to scattering effects are not always negligible (up to 7.7% underestimation). The proposed inversion scheme provides an opportunity to quantify gap fractions, LAI, and apparent clumping index even under sunny conditions.

a b s t r a c t Vegetation phenology is highly sensitive to climate change. Phenology also controls many feedbacks of vegetation to the climate system by influencing the seasonality of albedo, surface roughness length, canopy conductance, and fluxes of water, energy, CO2 and biogenic volatile organic compounds. In this review, we first discuss the environmental drivers of phenology, and the impacts of climate change on phenology, in different biomes. We then examine the vegetation-climate feedbacks that are mediated by phenology, and assess the potential impact on these feedbacks of shifts in phenology driven by climate change. We finish with an overview of phenological modeling and we suggest ways in which models might be improved using existing data sets. Several key weaknesses in our current understanding emerge from this analysis. First, we need a better understanding of the drivers of phenology, particularly in under-studied biomes (e.g. tropical forests). We do not have a mechanistic understanding of the role of photoperiod, even in well-studied biomes. In all biomes, the factors controlling senescence and dormancy are not well-documented. Second, for the most part (i.e. with the exception of phenology impacts on CO2 exchange) we have only a qualitative understanding of the feedbacks between vegetation and climate that are mediated by phenology. We need to quantify the magnitude of these feedbacks, and ensure that they are accurately reproduced by models. Third, we need to work towards a new understanding of phenological processes that enables progress beyond the modeling paradigms currently in use. Accurate representation of phenological processes in models that couple the land surface to the climate system is particularly important, especially when such models are being used to predict future climate.

Most land surface and ecosystem models assume that a vegetated canopy can be abstracted as a turbid medium when such models compute mass, energy, and carbon exchange. However, those models fail to simulate radiation environments in heterogeneous landscapes. This study aims to couple a spatially explicit three-dimensional (3D) shortwave and longwave radiative transfer model with a soil and canopy energy balance and canopy physiology model (CANOAK-FLiES), to investigate how well the 3D model performs in a heterogeneous landscape compared to a 1D model. The canopy structural parameters were extracted using airborne-based Light Detection and Ranging (LiDAR) data and digital cover photographs. The developed model was compared with a wide variety of field and remote sensing data (e.g., hyper-spectral remote sensing images, overstory and understory radiation, and tower-based energy and CO2 fluxes) in a woody oak savanna in CA, USA. Overall, the simulated spatial and diurnal patterns of the radiation environment were consistent with the measurements of shortwave and longwave wavebands. The 3D approach worked better than the 1D approach from the wet mild spring to the dry hot summer while the relative importance of using the 3D approach depends on climate and canopy physiological conditions as well as the canopy structure. With low leaf area index in the oak woodland, the woody elements absorbed 12%–39% of the PAR and 20%–52% of the NIR radiation. Consequently, 12% of the daytime energy flux was used as woody heat storage at our study site. This novel 3D model has the potential to serve as a useful tool for analyzing the spatio-temporal variability of radiation and energy fluxes in evaluating footprints of radiation sensors and eddy covariance fluxes, and serve as a standard in evaluating the performance of a hierarchy of simpler land surface models to compute mass and energy exchange of heterogeneous landscapes.

The continued spread of invasive weeds is threatening ecosystem health throughout North America. Understanding the relationships between invasive weeds’ key phenological phases and structural and/or functional canopy development is an essential step for making informed decisions regarding their management. We analyzed a three-year image archive obtained from an inexpensive webcam overlooking a perennial pepperweed (Lepidium latifoliumL.) infestation in California to explore the ability of red (R)‐green (G)‐blue (B) color space information to track the structural and functional development of the pepperweed.Wecharacterizedstructuralandfunctionalcanopydevelopmentthroughsurfaceroughness length (z 0m ; a proxy for canopy height and leaf area index) and canopy photosynthesis (F A ), respectively, both of which we derived from eddy covariance measurements. Here we demonstrate the use of cross-correlation functions to determine the temporal lags between chromatic coordinates and two color indices, all calculated fromRGBbrightness levels, withz 0m andF A . We found that these color metrics fail to represent the structural and/or functional state of the canopy. In contrast,relative luminance(CIE Y) appears to be a better indicator forz 0m and especially forF A . We calculatedCIE Yfrom pepperweedRGB brightness levels in relation to hypothetical horizontal referenceRGBbrightness levels. We obtained the latter by applying the ratio between horizontally measured and hypothetical incoming solar radiation on a vertical surface toRGBbrightness levels of a vertically oriented reference of invariant light-grey color. We conclude that webcam image archives may provide an inexpensive tool for making informed decisions regarding the timing but not for assessing the effectiveness of invasive plant control measures such as mowing.

We developed and tested the use of light emitting diodes (LEDs) to monitor vegetation reflectance in narrow spectral bands as a tool suitable for quantifying and monitoring ecosystem structure, function and metabolism. LEDs are appealing because they are inexpensive, small and reliable light sources that, when used in reverse mode, can measure spectrally selective radiation. We selected LEDs in red and near-infrared wavebands as they are used to calculate the normalized difference vegetation index (NDVI). The lab experiments revealed that the LEDs showed linear relation with a hyper-spectral spectroradiometer (R2 > 0.94 and 0.99 for red and NIR, respectively) and marginal sensitivity to temperature. To test the efficacy of this novel sensor, we measured spectral reflectance with LEDs over an annual grassland in California over 3.5 years. The LED-sensor captured daily to interannual variation of the spectral reflectance at the two bands with reliable and stable performance. The spectral reflectance in the two bands and NDVI proved to be useful to identify the leaf-on and leaf-off dates (mean bias errors of 5.3 and 4.2 days, respectively) and to estimate canopy photosynthesis (R2 = 0.91). We suggest that this novel instrument can monitor other structural and functional (e.g. leaf area index, leaf nitrogen) variables by employing LEDs that have other specific wavelengths bands. Considering that off-the-shelf LEDs cover a wide range of wavebands from the ultraviolet to near-infrared regions, we believe that the research community could explore a range of similar instruments across a range of bands for a variety of ecological applications.

Effective leaf area index is routinely quantified with optical instruments that measure gap fraction through the probability of beam penetration of sunlight through the vegetation. However, there have been few efforts to obtain theoretically consistent effective leaf area indices from those measurements. To apply the Beer–Lambert law, multiple gap fraction measurements may be averaged in two ways: (1) by taking the mean of the logarithms of the individual gap fraction values or (2) by taking the logarithm of the mean gap fraction. Based on a theoretical model and gap fraction measurements from 41 sites, we report that effective leaf area index must be quantified using the second approach. The first approach implemented in the LAI-2000 instrument considers clumping effects at scales larger than shoots. Thus, the combination of the first approach with an independent clumping index overestimates leaf area index up to 30% at the investigated sites. Clumping effects accounted for by the LAI-2000 instrument, called the “apparent” clumping index, were dependent on canopy cover, crown shape, and canopy height. A forest gap fraction model showed that short canopy height, vertically prolonged crown shape and higher canopy cover are associated with the lowest apparent clumping indices. We show that the apparent clumping index is a useful quantity to constrain the true clumping index and to investigate spatial and temporal variation of clumping effects. Such information would be useful to evaluate a coarse global clumping index map and improve land surface models.

Savannas are spatially heterogeneous, open ecosystems, thus efforts to quantify canopy structure with methods developed for homogeneous, closed canopies are prone to failure. We examine the applicability of two direct (litterfall, allometry) and five indirect (LAI-2000, TRAC, digital hemispheric photography, digital cover photography, traversing radiometer system) methods to determine leaf area index across a 9 ha domain in an oak-savanna ecosystem in California, USA. Interpretation of the leaf area index measurements is supported by two gap-fraction models. We recommend that leaf inclination angle distribution should be characterized first. For this purpose, we propose a simple, reliable and reproducible method using a digital camera. We show that the combination of digital cover photography and LAI-2000 could provide spatially representative leaf area index, gap fraction and element clumping index. Based on these two indirect methods, we quantify spatially representative element clumping index and leaf area index at ecosystem scale as 0.49 � 0.10 (mean � 95% confidence interval) and 0.77 � 0.27, respectively. In contrast to previous studies in northern ecosystems, measurement and modeling results suggest that element clumping index decreases with view zenith angle, most likely due to apparent changes of tree distribution pattern with the view zenith angle. Our results highlight the importance of ecosystem-scale clumping effects for the adequate quantification of tree leaf area index in savannas. Finally, we suggest a protocol to quantify leaf area index and its associated canopy structure variables in open canopy ecosystems.

Many studies on land surface radiation balances have relied on geostationary satellites. These satellites have provided data with high temporal resolution (less than 3 h); however, the spatial resolution was too coarse (20–250 km scale) to investigate local-scale land surface radiation balances. Moderate resolution imaging spectroradiometer (MODIS) – onboard both the Terra and Aqua satellites – yields a tradeoff with regard to this problem by providing higher spatial resolution (∼1 km scale) and sensing all over the earth nearly twice a day during daytime; this provides a potential tool for the periodical monitoring of the land surface energy balance. The reliability of MODIS-derived estimates is, however, affected by the presence of multiple error sources, such as those related to heterogeneous land cover and complex topography. In this study, we have used atmospheric (5 and 10 km scale) and land (1 km scale) products obtained from both the Terra and Aqua MODIS devices as inputs in order to estimate the radiation components (1 km scale) under clear daytime conditions over a heterogeneous farmland area and a rugged deciduous forest in the Korea Flux Network (KoFlux). The reliability of these estimates and the associated errors were evaluated by comparing against field measurements taken for 41 and 26 clear days with regard to the farmland and forest sites, respectively. Solar radiation was successfully retrieved with a root mean square error (RMSE) of ∼20 W m−2 for both the Terra and Aqua devices over the flat farmland site, whereas the rugged forest site exhibited corresponding values of 40 and 65 W m−2 RMSE values with consistent positive biases (presumably caused by topographic effects). The RMSE values of the downward longwave radiation were ∼20 W m−2 for both the Terra and Aqua devices for both these sites. The sensitivities of the upward components of the shortwave and longwave radiations varied with the RMSE values to the scale of the spatial heterogeneity of both the sites. Consequently, the RMSE values of the net radiation ranged from 33 to 61 W m−2 for both the devices at both the sites. Our results suggest that the scales of the patch mosaics within the landscapes need to be quantified for proper retrieval of the MODIS-derived radiation products. More extensive validation efforts are required to identify and account for major error sources across diverse land surface conditions.

a b s t r a c t The regular monitoring of evapotranspiration from satellites has been limited because of discontinu- ous temporal coverage, resulting in snapshots at a particular point in space and time. We developed a temporal upscaling scheme using satellite-derived instantaneous estimates of evapotranspiration to produce a daily-sum evapotranspiration averaged over an 8-day interval. We tested this scheme against measured evapotranspiration data from 34 eddy covariance flux towers covering seven plant functional types from boreal to tropical climatic zones. We found that the ratio of a half-hourly-sum of potential solar radiation (extraterrestrial solar irradiance on a plane parallel to the Earth's surface) between 10:00 hh and 14:00 hh to a daily-sum of potential solar radiation provides a robust scaling factor to convert a half-hourly measured evapotranspiration to an estimate of a daily-sum; the estimated and measured daily sum evapotranspiration showed strong linear relation (r 2 = 0.92) and small bias (−2.7%). By com- parison, assuming a constant evaporative fraction (the ratio of evapotranspiration to available energy) during the daytime, although commonly used for temporal upscaling, caused 13% underestimation of evapotranspiration on an annual scale. The proposed temporal upscaling scheme requires only latitude, longitude and time as input. Thus it will be useful for developing continuous evapotranspiration estimates in space and time, which will improve continuous monitoring of hydrological cycle from local to global scales.

International audience; We propose the Breathing Earth System Simulator (BESS), an upscaling approach to quantify global gross primary productivity and evapotranspiration using MODIS with a spatial resolution of 1-5 km and a temporal resolution of 8 days. This effort is novel because it is the first system that harmonizes and utilizes MODIS Atmosphere and Land products on the same projection and spatial resolution over the global land. This enabled us to use the MODIS Atmosphere products to calculate atmospheric radiative transfer for visual and near infrared radiation wave bands. Then we coupled atmospheric and canopy radiative transfer processes, with models that computed leaf photosynthesis, stomatal conductance and transpiration on the sunlit and shaded portions of the vegetation and soil. At the annual time step, the mass and energy fluxes derived from BESS showed strong linear relations with measurements of solar irradiance (r(2) = 0.95, relative bias: 8%), gross primary productivity (r(2) = 0.86, relative bias: 5%) and evapotranspiration (r(2) = 0.86, relative bias: 15%) in data from 33 flux towers that cover seven plant functional types across arctic to tropical climatic zones. A sensitivity analysis revealed that the gross primary productivity and evapotranspiration computed in BESS were most sensitive to leaf area index and solar irradiance, respectively. We quantified the mean global terrestrial estimates of gross primary productivity and evapotranpiration between 2001 and 2003 as 118 +/- 26 PgC yr(-1) and 500 +/- 104 mm yr(-1) (equivalent to 63,000 +/- 13,100 km(3) yr(-1)), respectively. BESS-derived gross primary productivity and evapotranspiration estimates were consistent with the estimates from independent machine-learning, data-driven products, but the process-oriented structure has the advantage of diagnosing sensitivity of mechanisms. The process-based BESS is able to offer gridded biophysical variables everywhere from local to the total global land scales with an 8-day interval over multiple years.

[1] We report on the interannual variability of evapotranspiration (E) and energy exchange of an annual grassland in the Mediterranean climate zone of California. They were measured directly with the eddy covariance technique over a 6-year period that spanned between July 2001 and June 2007 and experienced a large range in precipitation (376 mm to 888 mm). Despite a two-fold range in precipitation, annual E ranged much less, between 266 mm and 391 mm. We found that pronounced energy-limited and water-limited periods occurred within the same year. In the water-limited period, monthly integrated E scaled negatively with solar radiation and was restrained by precipitation. In the energy-limited period, on the other hand, the majority of E scaled positively with solar radiation (Rg) and was confined by potential E (Ep). E was most sensitive to the availability of soil moisture during the transition to the senescence period rather than onset of the greenness period, causing annual E to be strongly modulated by growing season length. Bulk surface conductance scaled consistently with Priestley-Taylor α coefficient regardless of interannual and seasonal variability of precipitation, E, and solar radiation.