Bruce D. Cook

Evaluation of remote sensing based terrestrial productivity from MODIS using regional tower eddy flux network observations

The Moderate Resolution Spectroradiometer (MODIS) sensor has provided near real-time estimates of gross primary production (GPP) since March 2000. We compare four years (2000 to 2003) of satellite-based calculations of GPP with tower eddy CO2 flux-based estimates across diverse land cover types and climate regimes. We examine the potential error contributions from meteorology, leaf area index (LAI)/fPAR, and land cover. The error between annual GPP computed from NASAs Data Assimilation Offices (DAO) and tower-based meteorology is 28%, indicating that NASAs DAO global meteorology plays an important role in the accuracy of the GPP algorithm. Approximately 62% of MOD15-based estimates of LAI were within the estimates based on field optical measurements, although remaining values overestimated site values. Land cover presented the fewest errors, with most errors within the forest classes, reducing potential error. Tower-based and MODIS estimates of annual GPP compare favorably for most biomes, although MODIS GPP overestimates tower-based calculations by 20%-30%. Seasonally, summer estimates of MODIS GPP are closest to tower data, and spring estimates are the worst, most likely the result of the relatively rapid onset of leaf-out. The results of this study indicate, however, that the current MODIS GPP algorithm shows reasonable spatial patterns and temporal variability across a diverse range of biomes and climate regimes. So, while continued efforts are needed to isolate particular problems in specific biomes, we are optimistic about the general quality of these data, and continuation of the MOD17 GPP product will likely provide a key component of global terrestrial ecosystem analysis, providing continuous weekly measurements of global vegetation production

Canopy nitrogen, carbon assimilation, and albedo in temperate and boreal forests: Functional relations and potential climate feedbacks

The availability of nitrogen represents a key constraint on carbon cycling in terrestrial ecosystems, and it is largely in this capacity that the role of N in the Earths climate system has been considered. Despite this, few studies have included continuous variation in plant N status as a driver of broad-scale carbon cycle analyses. This is partly because of uncertainties in how leaf-level physiological relationships scale to whole ecosystems and because methods for regional to continental detection of plant N concentrations have yet to be developed. Here, we show that ecosystem CO2 uptake capacity in temperate and boreal forests scales directly with whole-canopy N concentrations, mirroring a leaf-level trend that has been observed for woody plants worldwide. We further show that both CO2 uptake capacity and canopy N concentration are strongly and positively correlated with shortwave surface albedo. These results suggest that N plays an additional, and overlooked, role in the climate system via its influence on vegetation reflectivity and shortwave surface energy exchange. We also demonstrate that much of the spatial variation in canopy N can be detected by using broad-band satellite sensors, offering a means through which these findings can be applied toward improved application of coupled carbon cycle–climate models.

Amazon forests maintain consistent canopy structure and greenness during the dry season

The seasonality of sunlight and rainfall regulates net primary production in tropical forests. Previous studies have suggested that light is more limiting than water for tropical forest productivity, consistent with greening of Amazon forests during the dry season in satellite data. We evaluated four potential mechanisms for the seasonal green-up phenomenon, including increases in leaf area or leaf reflectance, using a sophisticated radiative transfer model and independent satellite observations from lidar and optical sensors. Here we show that the apparent green up of Amazon forests in optical remote sensing data resulted from seasonal changes in near-infrared reflectance, an artefact of variations in sun-sensor geometry. Correcting this bidirectional reflectance effect eliminated seasonal changes in surface reflectance, consistent with independent lidar observations and model simulations with unchanging canopy properties. The stability of Amazon forest structure and reflectance over seasonal timescales challenges the paradigm of light-limited net primary production in Amazon forests and enhanced forest growth during drought conditions. Correcting optical remote sensing data for artefacts of sun-sensor geometry is essential to isolate the response of global vegetation to seasonal and interannual climate variability.

Dissolved organic carbon in old field soils: Total amounts as a measure of available resources for soil mineralization

Dissolved organic carbon (DOC) and C and N mineralization were measured during a 210 day regulated in vitro incubation of soils from an old field successional sequence at Cedar Creek Natural History Area. The objective of the study was to evaluate the hypothesis that soil DOC constitutes a readily-available microbial resource, and that DOC concentrations are related to rates of biological decomposition and associated nutrient release from soil organic matter. Soils from five previously cultivated old fields undergoing secondary succession and an oak savanna were selected because they had demonstrated different patterns of C and N cycling. Although amounts of total C differed dramatically (496–1371 μmol g−1), DOC concentrations of all soils at the time of collection were between 0.70 and 1.30 μmol g−1. During the incubation, total and relative DOC concentrations generally remained constant or increased while mineralization rates decreased. When all soils and incubation intervals were considered, there was no obvious relationship between DOC and instantaneous rates of mineralization. Asymptotic exponential response curves did describe positive associations between DOC and CO2-C mineralization rates at early incubation times (R2 = 0.98 for 14 and 35 days), but not later. Similar models did not show a strong relationship between DOC and net-N mineralization rates. By the end of the incubation, the DOC pool could potentially supply 1.5–3.4 days of total C mineralization, but the instantaneous C mineralization rate at any given DOC concentration was 3–10 times lower than at 14 days. These results reflect decreased DOC utilization relative to supply, and could be caused by the accumulation of recalcitrant DOC.

NASA Goddards LiDAR, Hyperspectral and Thermal (G-LiHT) Airborne Imager

The combination of LiDAR and optical remotely sensed data provides unique information about ecosystem structure and function. Here, we describe the development, validation and application of a new airborne system that integrates commercial off the shelf LiDAR hyperspectral and thermal components in a compact, lightweight and portable system. Goddard’s LiDAR, Hyperspectral and Thermal (G-LiHT) airborne imager is a unique system that permits simultaneous measurements of vegetation structure, foliar spectra and surface temperatures at very high spatial resolution (~1 m) on a wide range of airborne platforms. The complementary nature of LiDAR, optical and thermal data provide an analytical framework for the development of new algorithms to map plant species composition, plant functional types, biodiversity, biomass and carbon stocks, and plant growth. In addition, G-LiHT data enhance our ability to validate data from existing satellite missions and support NASA Earth Science research. G-LiHT’s data processing and distribution system is designed to give scientists open access to both low- and high-level data products (http://gliht.gsfc.nasa.gov), which will stimulate the community development of synergistic data fusion algorithms. G-LiHT has been used to collect more than 6,500 km2 of data for NASA-sponsored studies across a broad range of ecoregions in the USA and Mexico. In this paper, we document G-LiHT design considerations, physical specifications, instrument performance and calibration and acquisition parameters. In addition, we describe the data processing system and higher-level data products that are freely distributed under NASA’s Data and Information policy.

Dissolved organic carbon in old field soils: Compositional changes during the biodegradation of soil organic matter

Abstract The quality of soil dissolved organic carbon (DOC) was examined using the Leenheer DOC fractionation scheme, which separates soluble organic compounds into well-defined functional groups that exhibit similar reactive properties. DOC fractions were measured for five previously cultivated old fields undergoing secondary succession, and an undisturbed oak savanna. Despite differences in field age (time since abandonment), plant community composition, distribution and amounts of phytomass, C and N storage, and potential amounts of CO 2 -C and net-N mineralization, the quality of soil DOC did not appear to differ. Nearly all DOC occurred in acid fractions (77%); hydrophilic acids alone constituted 50%. The fractionation procedure was also performed at four different times during a 210-day regulated in vitro incubation of the soils. Despite decreasing mineralization response to soil DOC concentrations, the fractional composition of the DOC remained relatively constant throughout the incubation. Although we could not evaluate DOC utilization, the results demonstrated that soil DOC was altered during the decomposition of soil organic matter; both total amounts and the relative N content of the hydrophobic acid fraction increased during the incubation period.

Integrating Solar Induced Fluorescence and the Photochemical Reflectance Index for Estimating Gross Primary Production in a Cornfield

The utilization of remotely sensed observations for light use efficiency (LUE) and tower-based gross primary production (GPP) estimates was studied in a USDA cornfield. Nadir hyperspectral reflectance measurements were acquired at canopy level during a collaborative field campaign conducted in four growing seasons. The Photochemical Reflectance Index (PRI) and solar induced chlorophyll fluorescence (SIF), were derived. SIF retrievals were accomplished in the two telluric atmospheric oxygen absorption features centered at 688 nm (O2-B) and 760 nm (O2-A). The PRI and SIF were examined in conjunction with GPP and LUE determined by flux tower-based measurements. All of these fluxes, environmental variables, and the PRI and SIF exhibited diurnal as well as day-to-day dynamics across the four growing seasons. Consistent with previous studies, the PRI was shown to be related to LUE (r2 = 0.54 with a logarithm fit), but the relationship varied each year. By combining the PRI and SIF in a linear regression model, stronger performances for GPP estimation were obtained. The strongest relationship (r2 = 0.80, RMSE = 0.186 mg CO2/m2/s) was achieved when using the PRI and SIF retrievals at 688 nm. Cross-validation approaches were utilized to demonstrate the robustness and consistency of the performance. This study highlights a GPP retrieval method based entirely on hyperspectral remote sensing observations.

Airborne Lidar-Based Estimates of Tropical Forest Structure in Complex Terrain: Opportunities and Trade-Offs for REDD+

BackgroundCarbon stocks and fluxes in tropical forests remain large sources of uncertainty in the global carbon budget. Airborne lidar remote sensing is a powerful tool for estimating aboveground biomass, provided that lidar measurements penetrate dense forest vegetation to generate accurate estimates of surface topography and canopy heights. Tropical forest areas with complex topography present a challenge for lidar remote sensing.ResultsWe compared digital terrain models (DTM) derived from airborne lidar data from a mountainous region of the Atlantic Forest in Brazil to 35 ground control points measured with survey grade GNSS receivers. The terrain model generated from full-density (~20 returns m−2) data was highly accurate (mean signed error of 0.19 ± 0.97 m), while those derived from reduced-density datasets (8 m−2, 4 m−2, 2 m−2 and 1 m−2) were increasingly less accurate. Canopy heights calculated from reduced-density lidar data declined as data density decreased due to the inability to accurately model the terrain surface. For lidar return densities below 4 m−2, the bias in height estimates translated into errors of 80–125 Mg ha−1 in predicted aboveground biomass.ConclusionsGiven the growing emphasis on the use of airborne lidar for forest management, carbon monitoring, and conservation efforts, the results of this study highlight the importance of careful survey planning and consistent sampling for accurate quantification of aboveground biomass stocks and dynamics. Approaches that rely primarily on canopy height to estimate aboveground biomass are sensitive to DTM errors from variability in lidar sampling density.

DART: Recent Advances in Remote Sensing Data Modeling With Atmosphere, Polarization, and Chlorophyll Fluorescence

To better understand the life-essential cycles and processes of our planet and to further develop remote sensing (RS) technology, there is an increasing need for models that simulate the radiative budget (RB) and RS acquisitions of urban and natural landscapes using physical approaches and considering the three-dimensional (3-D) architecture of Earth surfaces. Discrete anisotropic radiative transfer (DART) is one of the most comprehensive physically based 3-D models of Earth-atmosphere radiative transfer, covering the spectral domain from ultraviolet to thermal infrared wavelengths. It simulates the optical 3-D RB and optical signals of proximal, aerial, and satellite imaging spectrometers and laser scanners, for any urban and/or natural landscapes and for any experimental and instrumental configurations. It is freely available for research and teaching activities. In this paper, we briefly introduce DART theory and present recent advances in simulated sensors (LiDAR and cameras with finite field of view) and modeling mechanisms (atmosphere, specular reflectance with polarization and chlorophyll fluorescence). A case study demonstrating a novel application of DART to investigate urban landscapes is also presented.

Structural Dynamics of Tropical Moist Forest Gaps.

Gap phase dynamics are the dominant mode of forest turnover in tropical forests. However, gap processes are infrequently studied at the landscape scale. Airborne lidar data offer detailed information on three-dimensional forest structure, providing a means to characterize fine-scale (1 m) processes in tropical forests over large areas. Lidar-based estimates of forest structure (top down) differ from traditional field measurements (bottom up), and necessitate clear-cut definitions unencumbered by the wisdom of a field observer. We offer a new definition of a forest gap that is driven by forest dynamics and consistent with precise ranging measurements from airborne lidar data and tall, multi-layered tropical forest structure. We used 1000 ha of multi-temporal lidar data (2008, 2012) at two sites, the Tapajos National Forest and Ducke Reserve, to study gap dynamics in the Brazilian Amazon. Here, we identified dynamic gaps as contiguous areas of significant growth, that correspond to areas > 10 m2, with height <10 m. Applying the dynamic definition at both sites, we found over twice as much area in gap at Tapajos National Forest (4.8 %) as compared to Ducke Reserve (2.0 %). On average, gaps were smaller at Ducke Reserve and closed slightly more rapidly, with estimated height gains of 1.2 m y-1 versus 1.1 m y-1 at Tapajos. At the Tapajos site, height growth in gap centers was greater than the average height gain in gaps (1.3 m y-1 versus 1.1 m y-1). Rates of height growth between lidar acquisitions reflect the interplay between gap edge mortality, horizontal ingrowth and gap size at the two sites. We estimated that approximately 10 % of gap area closed via horizontal ingrowth at Ducke Reserve as opposed to 6 % at Tapajos National Forest. Height loss (interpreted as repeat damage and/or mortality) and horizontal ingrowth accounted for similar proportions of gap area at Ducke Reserve (13 % and 10 %, respectively). At Tapajos, height loss had a much stronger signal (23 % versus 6 %) within gaps. Both sites demonstrate limited gap contagiousness defined by an increase in the likelihood of mortality in the immediate vicinity (~6 m) of existing gaps.