Satellite Ecology (SATECO)—linking ecology,remote sensing and micrometeorology,from plot to regional scale,for the study of ecosystem structure and function

There is a growing requirement for ecosystem science to help inform a deeper understanding of the effects of global climate change and land use change on terrestrial ecosystem structure and function, from small area (plot) to landscape, regional and global scales. To meet these requirements, ecologists have investigated plant growth and carbon cycling processes at plot scale, using biometric methods to measure plant carbon accumulation, and gas exchange (chamber) methods to measure soil respiration. Also at the plot scale, micrometeorologists have attempted to measure canopy- or ecosystem-scale CO₂ flux by the eddy covariance technique, which reveals diurnal, seasonal and annual cycles. Mathematical models play an important role in integrating ecological and micrometeorological processes into ecosystem scales, which are further useful in interpreting time-accumulated information derived from biometric methods by comparing with CO₂ flux measurements. For a spatial scaling of such plot-level understanding, remote sensing via satellite is used to measure land use/vegetation type distribution and temporal changes in ecosystem structures such as leaf area index. However, to better utilise such data, there is still a need for investigations that consider the structure and function of ecosystems and their processes, especially in mountainous areas characterized by complex terrain and a mosaic distribution of vegetation. For this purpose, we have established a new interdisciplinary approach named ‘Satellite Ecology’, which aims to link ecology, remote sensing and micrometeorology to facilitate the study of ecosystem function, at the plot, landscape, and regional scale. This article was contributed at the invitation of the Editorial Committee.     

新一代遥感技术助力生态系统生态学研究

随着气候变化和人类活动的加剧, 生态系统正处于剧烈变化中, 生态学家需要从更大的时空尺度去理解生态系统过程和变化规律, 应对全球变化带来的威胁和挑战。传统地面调查方法主要获取的是样方尺度、离散的数据, 难以满足大尺度生态系统研究对数据时空连续性的要求。相比于传统地面调查方法, 遥感技术具有实时获取、重复监测以及多时空尺度的特点, 弥补了传统地面调查方法空间观测尺度有限的缺点。遥感通过分析电磁波信息从而识别地物属性和特征, 反演生态系统组成、能量流动和物质循环过程中的关键要素, 已逐渐成为生态学研究中必不可少的数据来源。近年来, 随着激光雷达、日光诱导叶绿素荧光等新型遥感技术以及无人机、背包等近地面遥感平台的发展, 个人化、定制化的近地面遥感观测逐渐成熟, 新一代遥感技术正在推动遥感信息“二维向三维”的转变, 为传统样地观测与卫星遥感之间搭建了尺度推绎桥梁, 这也给生态系统生态学带来了新的机遇, 推动生态系统生态学向多尺度、多过程、多学科、多途径发展。因此, 该文从生态系统生态学角度出发, 重点关注陆地生态系统中生物组分, 并分别从生态系统类型、结构、功能和生物多样性等方面, 结合作者在实际研究工作中的主要成果和该领域国际前沿动态, 阐述遥感技术在生态系统生态学中的研究现状并指出我国生态系统遥感监测领域发展方向及亟待解决的问题。     

Upscaling leaf area index in an Arctic landscape through multiscale observations

Monitoring and understanding global change requires a detailed focus on upscaling, the process for extrapolating from the site‐specific scale to the smallest scale resolved in regional or global models or earth observing systems. Leaf area index (LAI) is one of the most sensitive determinants of plant production and can vary by an order of magnitude over short distances. The landscape distribution of LAI is generally determined by remote sensing of surface reflectance (e.g. normalized difference vegetation index, NDVI) but the mismatch in scales between ground and satellite measurements complicates LAI upscaling. Here, we describe a series of measurements to quantify the spatial distribution of LAI in a sub‐Arctic landscape and then describe the upscaling process and its associated errors. Working from a fine‐scale harvest LAI–NDVI relationship, we collected NDVI data over a 500 m × 500 m catchment in the Swedish Arctic, at resolutions from 0.2 to 9.0 m in a nested sampling design. NDVI scaled linearly, so that NDVI at any scale was a simple average of multiple NDVI measurements taken at finer scales. The LAI–NDVI relationship was scale invariant from 1.5 to 9.0 m resolution. Thus, a single exponential LAI–NDVI relationship was valid at all these scales, with similar prediction errors. Vegetation patches were of a scale of ～0.5 m and at measurement scales coarser than this, there was a sharp drop in LAI variance. Landsat NDVI data for the study catchment correlated significantly, but poorly, with ground‐based measurements. A variety of techniques were used to construct LAI maps, including interpolation by inverse distance weighting, ordinary Kriging, External Drift Kriging using Landsat data, and direct estimation from a Landsat NDVI–LAI calibration. All methods produced similar LAI estimates and overall errors. However, Kriging approaches also generated maps of LAI estimation error based on semivariograms. The spatial variability of this Arctic landscape was such that local measurements assimilated by Kriging approaches had a limited spatial influence. Over scales >50 m, interpolation error was of similar magnitude to the error in the Landsat NDVI calibration. The characterisation of LAI spatial error in this study is a key step towards developing spatio‐temporal data assimilation systems for assessing C cycling in terrestrial ecosystems by combining models with field and remotely sensed data.     

Integration of ground survey and remote sensing derived data: Producing robust indicators of habitat extent and condition

The availability of suitable habitat is a key predictor of the changing status of biodiversity. Quantifying habitat availability over large spatial scales is, however, challenging. Although remote sensing techniques have high spatial coverage, there is uncertainty associated with these estimates due to errors in classification. Alternatively, the extent of habitats can be estimated from ground‐based field survey. Financial and logistical constraints mean that on‐the‐ground surveys have much lower coverage, but they can produce much higher quality estimates of habitat extent in the areas that are surveyed. Here, we demonstrate a new combined model which uses both types of data to produce unified national estimates of the extent of four key habitats across Great Britain based on Countryside Survey and Land Cover Map. This approach considers that the true proportion of habitat per km² (Z_i) is unobserved, but both ground survey and remote sensing can be used to estimate Z_i. The model allows the relationship between remote sensing data and Z_i to be spatially biased while ground survey is assumed to be unbiased. Taking a statistical model‐based approach to integrating field survey and remote sensing data allows for information on bias and precision to be captured and propagated such that estimates produced and parameters estimated are robust and interpretable. A simulation study shows that the combined model should perform best when error in the ground survey data is low. We use repeat surveys to parameterize the variance of ground survey data and demonstrate that error in this data source is small. The model produced revised national estimates of broadleaved woodland, arable land, bog, and fen, marsh and swamp extent across Britain in 2007.     

遥感主体图的准确度对景观生态学研究的影响

用各种案例系统地解释了遥感数据分类误差对景观指数误差的必然影响。一方面 ,遥感数据在各种时间和空间尺度上为景观生态学研究提供必需的土地类型数据 ;另一方面 ,遥感技术的灵活性和复杂性可以产生出各种质量的土地类型数据。但景观生态学方面的用户对土地类型数据基本上是没有选择地使用 ,甚至是不知好坏地使用 ,所以景观生态学的发现和结论具有不可避免的任意性。总结了在各种情况下景观指数的变动区间 ,指出了现实较低的遥感数据的分类准确度会引起更低的景观指数的准确度 ,当进行景观变化分析时 ,这种误差的放大效应将更加明显。当前 ,人们对除面积以外的景观指数的误差仍然束手无策 ,尽可能地提高遥感数据的分类准确度是唯一力所能及的办法。     

Model–data synthesis in terrestrial carbon observation: methods, data requirements and data uncertainty specifications

Systematic, operational, long‐term observations of the terrestrial carbon cycle (including its interactions with water, energy and nutrient cycles and ecosystem dynamics) are important for the prediction and management of climate, water resources, food resources, biodiversity and desertification. To contribute to these goals, a terrestrial carbon observing system requires the synthesis of several kinds of observation into terrestrial biosphere models encompassing the coupled cycles of carbon, water, energy and nutrients. Relevant observations include atmospheric composition (concentrations of CO₂ and other gases); remote sensing; flux and process measurements from intensive study sites; in situ vegetation and soil monitoring; weather, climate and hydrological data; and contemporary and historical data on land use, land use change and disturbance (grazing, harvest, clearing, fire). A review of model–data synthesis tools for terrestrial carbon observation identifies ‘nonsequential’ and ‘sequential’ approaches as major categories, differing according to whether data are treated all at once or sequentially. The structure underlying both approaches is reviewed, highlighting several basic commonalities in formalism and data requirements. An essential commonality is that for all model–data synthesis problems, both nonsequential and sequential, data uncertainties are as important as data values themselves and have a comparable role in determining the outcome. Given the importance of data uncertainties, there is an urgent need for soundly based uncertainty characterizations for the main kinds of data used in terrestrial carbon observation. The first requirement is a specification of the main properties of the error covariance matrix. As a step towards this goal, semi‐quantitative estimates are made of the main properties of the error covariance matrix for four kinds of data essential for terrestrial carbon observation: remote sensing of land surface properties, atmospheric composition measurements, direct flux measurements, and measurements of carbon stores.     

Parameter estimation in surface exchange models using nonlinear inversion: how many parameters can we estimate and which measurements are most useful?

Models of mass and energy exchanges between the biosphere and the atmosphere generally contain a nonlinear dependence between fluxes and model parameters, and thus estimation of these parameters from measurements in a heterogeneous landscape depends on the scale of the observations. The scale‐dependence of a typical surface‐exchange model (the CSIRO Biospheric Model, CBM) is examined using the diurnal variation of hourly fluxes of CO₂, latent heat, sensible heat and soil heat. The fluxes were measured using micrometeorological techniques over six sites in a grazing/pasture system in SE Australia during a period of three weeks in 1995. Nonlinear parameter inversion was used to determine model parameters. Analysis of the covariance of the estimates of the parameters and the unexplained residuals of the model showed that a maximum of three or four parameters could be determined independently from the observations for all six sites. Estimates of a key model parameter, j_max, the mean of maximum potential electron transport rate of all leaves within the canopy, was best determined by the measurements of net CO₂ flux at all sites examined. Measurements of ground heat flux provide little information about any of the model parameters in CBM. Because of nonlinearities in the surface exchange model, calculated fluxes will be in error if parameters for the component vegetation types are simply averaged in proportion to their areal fraction. The magnitude of these errors was examined for CBM using a hypothetical land surface consisting of two surface types, each with different parameter values. Predictions of net CO₂, latent heat and ground heat fluxes using a linear combination of model parameters for the two surface types were quite similar with those found using optimal estimates of the parameters for the landscape, but were significantly poorer for sensible heat fluxes.     

A cross-biome comparison of daily light use efficiency for gross primary production

Vegetation light use efficiency is a key physiological parameter at the canopy scale, and at the daily time step is a component of remote sensing algorithms for scaling gross primary production (GPP) and net primary production (NPP) over regional to global domains. For the purposes of calibrating and validating the light use efficiency ( ε _g) algorithms, the components of ε _g– absorbed photosynthetically active radiation (APAR) and ecosystem GPP – must be measured in a variety of environments. Micrometeorological and mass flux measurements at eddy covariance flux towers can be used to estimate APAR and GPP, and the emerging network of flux tower sites offers the opportunity to investigate spatial and temporal patterns in ε _g at the daily time step. In this study, we examined the relationship of daily GPP to APAR, and relationships of ε _g to climatic variables, at four micrometeorological flux tower sites – an agricultural field, a tallgrass prairie, a deciduous forest, and a boreal forest. The relationship of GPP to APAR was close to linear at the tallgrass prairie site but more nearly hyperbolic at the other sites. The sites differed in the mean and range of daily ε _g, with higher values associated with the agricultural field than the boreal forest. ε_g decreased with increasing APAR at all sites, a function of mid‐day saturation of GPP and higher ε _g under overcast conditions. ε _g was generally not well correlated with vapor pressure deficit or maximum daily temperature. At the agricultural site, a ε _g decline towards the end of the growing season was associated with a decrease in foliar nitrogen concentration. At the tallgrass prairie site, a decline in ε _g in August was associated with soil drought. These results support inclusion of parameters for cloudiness and the phenological status of the vegetation, as well as use of biome‐specific parameterization, in operational ε _g algorithms.     

Ground and remote sensing-based measurements of leaf area index in a transitional forest and seasonal flooded forest in Brazil

Leaf area index (LAI) is a key driver of forest productivity and evapotranspiration; however, it is a difficult and labor-intensive variable to measure, making its measurement impractical for large-scale and long-term studies of tropical forest structure and function. In contrast, satellite estimates of LAI have shown promise for large-scale and long-term studies, but their performance has been equivocal and the biases are not well known. We measured total, overstory, and understory LAI of an Amazon-savanna transitional forest (ASTF) over 3 years and a seasonal flooded forest (SFF) during 4 years using a light extinction method and two remote sensing methods (LAI MODIS product and the Landsat-METRIC method), with the objectives of (1) evaluating the performance of the remote sensing methods, and (2) understanding how total, overstory and understory LAI interact with micrometeorological variables. Total, overstory and understory LAI differed between both sites, with ASTF having higher LAI values than SFF, but neither site exhibited year-to-year variation in LAI despite large differences in meteorological variables. LAI values at the two sites have different patterns of correlation with micrometeorological variables. ASTF exhibited smaller seasonal variations in LAI than SFF. In contrast, SFF exhibited small changes in total LAI; however, dry season declines in overstory LAI were counteracted by understory increases in LAI. MODIS LAI correlated weakly to total LAI for SFF but not for ASTF, while METRIC LAI had no correlation to total LAI. However, MODIS LAI correlated strongly with overstory LAI for both sites, but had no correlation with understory LAI. Furthermore, LAI estimates based on canopy light extinction were correlated positively with seasonal variations in rainfall and soil water content and negatively with vapor pressure deficit and solar radiation; however, in some cases satellite-derived estimates of LAI exhibited no correlation with climate variables (METRIC LAI or MODIS LAI for ASTF). These data indicate that the satellite-derived estimates of LAI are insensitive to the understory variations in LAI that occur in many seasonal tropical forests and the micrometeorological variables that control seasonal variations in leaf phenology. While more ground-based measurements are needed to adequately quantify the performance of these satellite-based LAI products, our data indicate that their output must be interpreted with caution in seasonal tropical forests.     

Comparison and Assessment of Aerial and Ground Estimates of Waterbird Colonies

Abstract: Aerial surveys are often used to quantify sizes of waterbird colonies; however, these surveys would benefit from a better understanding of associated biases. We compared estimates of breeding pairs of waterbirds, in colonies across southern Louisiana, USA, made from the ground, fixed-wing aircraft, and a helicopter. We used a marked-subsample method for ground-counting colonies to obtain estimates of error and visibility bias. We made comparisons over 2 sampling periods: 1) surveys conducted on the same colonies using all 3 methods during 3–11 May 2005 and 2) an expanded fixed-wing and ground-survey comparison conducted over 4 periods (May and Jun, 2004–2005). Estimates from fixed-wing aircraft were approximately 65% higher than those from ground counts for overall estimated number of breeding pairs and for both dark and white-plumaged species. The coefficient of determination between estimates based on ground and fixed-wing aircraft was ≤0.40 for most species, and based on the assumption that estimates from the ground were closer to the true count, fixed-wing aerial surveys appeared to overestimate numbers of nesting birds of some species; this bias often increased with the size of the colony. Unlike estimates from fixed-wing aircraft, numbers of nesting pairs made from ground and helicopter surveys were very similar for all species we observed. Ground counts by one observer resulted in underestimated number of breeding pairs by 20% on average. The marked-subsample method provided an estimate of the number of missed nests as well as an estimate of precision. These estimates represent a major advantage of marked-subsample ground counts over aerial methods; however, ground counts are difficult in large or remote colonies. Helicopter surveys and ground counts provide less biased, more precise estimates of breeding pairs than do surveys made from fixed-wing aircraft. We recommend managers employ ground counts using double observers for surveying waterbird colonies when feasible. Fixed-wing aerial surveys may be suitable to determine colony activity and composition of common waterbird species. The most appropriate combination of survey approaches will be based on the need for precise and unbiased estimates, balanced with financial and logistical constraints. (JOURNAL OF WILDLIFE MANAGEMENT 72(3):697–706; 2008)     

基于地面调查与遥感调查协同的城市绿地规模和分布研究——以青岛市为例

城市绿地调查作为园林绿化规划建设管理的基础工作，其数据准确程度直接影响着城市园林绿化的科学发展。提出了基于遥感调查与地面调查协同的城市绿地调查的技术方法，通过该方法获得数据与单独使用遥感调查或地面调查技术获得数据相比提升了城市绿地，尤其是公园绿地和附属绿地调查的时效性和准确性，为城市绿地调查提供了可操作的技术方法。     

Agricultural land-use in China: a comparison of area estimates from ground-based census and satellite-borne remote sensing

We compare estimates of total cropland area, paddy rice area, and irrigated cropland area in China from land cover maps derived from optical remote sensing in 1992–93 (1-km resolution NOAA AVHRR) and county-level agricultural census data for 1990. At national, regional, provincial, and county scales, the total cropland area estimated by remote sensing is 50–100% higher than reported in the agricultural census. For topographically flat North and Central China, there is a high correlation between county-level cropland area estimates by the two approaches. For other regions, the correlation between remote sensing and agricultural census cropland area is much weaker. Throughout China, there is only moderate to weak correlation between remote sensing-based and census-bases estimates of paddy rice area and total irrigated cropland area. It is likely that the census data underestimates and the remote sensing overestimates the actual cropland area. These uncertainties in agricultural land cover estimates will contribute to uncertainty in any large-scale biogeochemical analyses.     

Satellite-Based Analysis of Evapotranspiration and Water Balance in the Grassland Ecosystems of Dryland East Asia

The regression tree method is used to upscale evapotranspiration (ET) measurements at eddy-covariance (EC) towers to the grassland ecosystems over the Dryland East Asia (DEA). The regression tree model was driven by satellite and meteorology datasets, and explained 82% and 76% of the variations of ET observations in the calibration and validation datasets, respectively. The annual ET estimates ranged from 222.6 to 269.1 mm yr⁻¹ over the DEA region with an average of 245.8 mm yr⁻¹ from 1982 through 2009. Ecosystem ET showed decreased trends over 61% of the DEA region during this period, especially in most regions of Mongolia and eastern Inner Mongolia due to decreased precipitation. The increased ET occurred primarily in the western and southern DEA region. Over the entire study area, water balance (the difference between precipitation and ecosystem ET) decreased substantially during the summer and growing season. Precipitation reduction was an important cause for the severe water deficits. The drying trend occurring in the grassland ecosystems of the DEA region can exert profound impacts on a variety of terrestrial ecosystem processes and functions.     

Spectral response of the seagrass Zostera noltii with different sediment backgrounds

The efficiency of vegetation indices (VIs) to estimate the above-ground biomass of the seagrass species Zostera noltii Hornem. from remote sensing was tested experimentally on different substrata, since terrestrial vegetation studies have shown that VIs can be adversely influenced by the spectral properties of soils and background surfaces. Leaves placed on medium sand, fine sand and autoclaved fine sand were incrementally removed, and the spectral reflectance was measured in the 400–900 nm wavelength range. Several VIs were evaluated: ratios using visible and near infrared wavelengths, narrow-band indices, indices based on derivative analysis and continuum removal. Background spectral reflectance was clearly visible in the leaf reflectance spectra, showing marked brightness and spectral contrast variations for the same amount of vegetation. Paradoxically, indices used to minimize soil effects, such as the Soil-Adjusted Vegetation Index (SAVI) and the Modified second Soil-Adjusted Vegetation Index (MSAVI₂) showed a high sensitivity to background effects. Similar results were found for the widely used Normalized Difference Vegetation Index (NDVI) and for Pigment Specific Simple Ratios (PSSRs). In fact, background effects were most reduced for VIs integrating a blue band correction, namely the modified specific ratio (mSR₍₇₀₅₎), the modified Normalized Difference (mND₍₇₀₅₎), and two modified NDVIs proposed in this study. However, these indices showed a faster saturation for high seagrass biomass. The background effects were also substantially reduced using Modified Gaussian Model indices at 620 and 675 nm. The blue band corrected VIs should now be tested for air-borne or satellite remote sensing applications, but some require sensors with a hyperspectral resolution. Nevertheless, this type of index can be applied to analyse broad band multispectral satellite images with a blue band.     

Inversion of net ecosystem CO2 flux measurements for estimation of canopy PAR absorption

The fractional absorption of photosynthetically active radiation (f_PAR) is frequently a key variable in models describing terrestrial ecosystem–atmosphere interactions, carbon uptake, growth and biogeochemistry. We present a novel approach to the estimation of the fraction of incident photosynthetically active radiation absorbed by the photosynthetic components of a plant canopy (f_Chl). The method uses micrometeorological measurements of CO₂ flux and incident radiation to estimate light response parameters from which canopy structure is deduced. Data from two Ameriflux sites in Oklahoma, a tallgrass prairie site and a wheat site, are used to derive 7‐day moving average estimates of f_Chl during three years (1997–1999). The inverse estimates are compared to long‐term field measurements of PAR absorption. Good correlations are obtained when the field‐measured f_PAR is scaled by an estimate of the green fraction of total leaf area, although the inverse technique tends to be lower in value than the field measurements. The inverse estimates of f_Chl using CO₂ flux measurements are different from measurements of f_PAR that might be made by other, more direct, techniques. However, because the inverse estimates are based on observed canopy CO₂ uptake, they might be considered more biologically relevant than direct measurements that are affected by non‐physiologically active components of the canopy. With the increasing number of eddy covariance sites around the world the technique provides the opportunity to examine seasonal and inter‐annual variation in canopy structure and light harvesting capacity at individual sites. Furthermore, the inverse f_Chl provide a new source of data for development and testing of f_PAR retrieval using remote sensing. New remote sensing algorithms, or adjustments to existing algorithms, might thus become better conditioned to ‘biologically significant’ light absorption than currently possible.     

北京市农业景观生态与美学质量空间评价

选择了重点反映农业景观生态与美学功能的生态系统功能、自然性、开阔与多样性、污染概率和整洁度等10个单一空间显性指数,基于田间调查和专家系统赋予其不同权重,构建了农业景观质量综合评价指标,并基于地理信息系统支持下的土地利用数据和利用遥感获得的植被指数,对北京市农业景观质量及其空间差异进行了评价.结果表明：北京市农业景观质量的区域差别较大,城区及小城镇边缘的农业景观质量最差,近郊外围及远郊区较好;北京市农业景观的优劣主要与地形和人为压力有关.基于空间信息技术和景观指数,在大尺度上构建的综合指标体系,对评价农业景观质量的空间差异具有重要意义.     

Airborne hyperspectral data predict Ellenberg indicator values for nutrient and moisture availability in dry grazed grasslands within a local agricultural landscape

Species-based ecological indices, such as Ellenberg indicators, reflect plant habitat preferences and can be used to describe local environment conditions. One disadvantage of using vegetation data as a substitute for environmental data is the fact that extensive floristic sampling can usually only be carried out at a plot scale within limited geographical areas. Remotely sensed data have the potential to provide information on fine-scale vegetation properties over large areas. In the present study, we examine whether airborne hyperspectral remote sensing can be used to predict Ellenberg nutrient (N) and moisture (M) values in plots in dry grazed grasslands within a local agricultural landscape in southern Sweden. We compare the prediction accuracy of three categories of model: (I) models based on predefined vegetation indices (VIs), (II) models based on waveband-selected VIs, and (III) models based on the full set of hyperspectral wavebands. We also identify the optimal combination of wavebands for the prediction of Ellenberg values. The floristic composition of 104 (4 m × 4 m grassland) plots on the Baltic island of Öland was surveyed in the field, and the vascular plant species recorded in the plots were assigned Ellenberg indicator values for N and M. A community-weighted mean value was calculated for N (mN) and M (mM) within each plot. Hyperspectral data were extracted from an 8 m × 8 m pixel window centred on each plot. The relationship between field-observed and predicted mean Ellenberg values was significant for all three categories of prediction models. The performance of the category II and III models was comparable, and they gave lower prediction errors and higher R² values than the category I models for both mN and mM. Visible and near-infrared wavebands were important for the prediction of both mN and mM, and shortwave infrared wavebands were also important for the prediction of mM. We conclude that airborne hyperspectral remote sensing can detect spectral differences in vegetation between grassland plots characterised by different mean Ellenberg N and M values, and that remote sensing technology can potentially be used to survey fine-scale variation in environmental conditions within a local agricultural landscape.     

Site-level evaluation of satellite-based global terrestrial gross primary production and net primary production monitoring

Operational monitoring of global terrestrial gross primary production (GPP) and net primary production (NPP) is now underway using imagery from the satellite‐borne Moderate Resolution Imaging Spectroradiometer (MODIS) sensor. Evaluation of MODIS GPP and NPP products will require site‐level studies across a range of biomes, with close attention to numerous scaling issues that must be addressed to link ground measurements to the satellite‐based carbon flux estimates. Here, we report results of a study aimed at evaluating MODIS NPP/GPP products at six sites varying widely in climate, land use, and vegetation physiognomy. Comparisons were made for twenty‐five 1 km² cells at each site, with 8‐day averages for GPP and an annual value for NPP. The validation data layers were made with a combination of ground measurements, relatively high resolution satellite data (Landsat Enhanced Thematic Mapper Plus at ～30 m resolution), and process‐based modeling. There was strong seasonality in the MODIS GPP at all sites, and mean NPP ranged from 80 g C m^?2 yr^?1 at an arctic tundra site to 550 g C m^?2 yr^?1 at a temperate deciduous forest site. There was not a consistent over‐ or underprediction of NPP across sites relative to the validation estimates. The closest agreements in NPP and GPP were at the temperate deciduous forest, arctic tundra, and boreal forest sites. There was moderate underestimation in the MODIS products at the agricultural field site, and strong overestimation at the desert grassland and at the dry coniferous forest sites. Analyses of specific inputs to the MODIS NPP/GPP algorithm – notably the fraction of photosynthetically active radiation absorbed by the vegetation canopy, the maximum light use efficiency (LUE), and the climate data – revealed the causes of the over‐ and underestimates. Suggestions for algorithm improvement include selectively altering values for maximum LUE (based on observations at eddy covariance flux towers) and parameters regulating autotrophic respiration.     

Assimilation exceeds respiration sensitivity to drought: A FLUXNET synthesis

The intensification of the hydrological cycle, with an observed and modeled increase in drought incidence and severity, underscores the need to quantify drought effects on carbon cycling and the terrestrial sink. FLUXNET, a global network of eddy covariance towers, provides dense data streams of meteorological data, and through flux partitioning and gap filling algorithms, estimates of net ecosystem productivity (F_NEP), gross ecosystem productivity (P), and ecosystem respiration (R). We analyzed the functional relationship of these three carbon fluxes relative to evaporative fraction (EF), an index of drought and site water status, using monthly data records from 238 micrometeorological tower sites distributed globally across 11 biomes. The analysis was based on relative anomalies of both EF and carbon fluxes and focused on drought episodes by biome and climatic season. Globally P was ≈50% more sensitive to a drought event than R. Network‐wide drought‐induced decreases in carbon flux averaged ?16.6 and ?9.3 g C m^?2 month^?1 for P and R, i.e., drought events induced a net decline in the terrestrial sink. However, in evergreen forests and wetlands drought was coincident with an increase in P or R during parts of the growing season. The most robust relationships between carbon flux and EF occurred during climatic spring for F_NEP and in climatic summer for P and R. Upscaling flux sensitivities to a global map showed that spatial patterns for all three carbon fluxes were linked to the distribution of croplands. Agricultural areas exhibited the highest sensitivity whereas the tropical region had minimal sensitivity to drought. Combining gridded flux sensitivities with their uncertainties and the spatial grid of FLUXNET revealed that a more robust quantification of carbon flux response to drought requires additional towers in all biomes of Africa and Asia as well as in the cropland, shrubland, savannah, and wetland biomes globally.     

Application of two remote sensing GPP algorithms at a semiarid grassland site of North China

Aims Estimation of gross primary production (GPP) from remote sensing data is an important approach to study regional or global carbon cycle. However, for a given algorithm, it usually has its limitation on applications to a wide range of vegetation types and/or under diverse environmental conditions. This study was conducted to compare the performance of two remote sensing GPP algorithms, the MODIS GPP and the vegetation photosynthesis model (VPM), in a semiarid temperate grassland ecosystem.Methods The study was conducted at a typical grassland site in Ujimuqin of Inner Mongolia, North China, over 2 years in 2006 and 2007. Environmental controls on GPP measured by the eddy covariance (EC) technique at the study site were first investigated with path analysis of meteorological and soil moisture data at a daily and 8-day time steps. The estimates of GPP derived from the MODIS GPP and the VPM with site-specific inputs were then compared with the values of EC measurements as ground truthing at the site. Site-specific ? max (α) was estimated by using rectangular hyperbola function based on the 7-day flux data at 30-min intervals over the peak period of the growing season (May to September).Important findings Between the two remote sensing GPP algorithms and various estimates of the fraction of absorbed photosynthetic active radiation (FPAR), the VPM based on FPAR derived from the enhanced vegetation index (EVI) works the best in predicting GPP against the ground truthing of EC GPP. A path analysis indicates that the EC GPP in this semiarid temperate grassland ecosystem is controlled predominantly by both soil water and temperature. The site water condition is slightly better simulated by the moisture multiplier in the VPM than in the MODIS GPP algorithm, which is a most probable explanation for a better performance of the VPM than MODIS GPP algorithm in this semiarid grassland ecosystem.