Skip to main content

Net primary productivity of forest stands in New Hampshire estimated from Landsat and MODIS satellite data



A simulation model that relies on satellite observations of vegetation cover from the Landsat 7 sensor and from the Moderate Resolution Imaging Spectroradiometer (MODIS) was used to estimate net primary productivity (NPP) of forest stands at the Bartlett Experiment Forest (BEF) in the White Mountains of New Hampshire.


Net primary production (NPP) predicted from the NASA-CASA model using 30-meter resolution Landsat inputs showed variations related to both vegetation cover type and elevational effects on mean air temperatures. Overall, the highest predicted NPP from the NASA-CASA model was for deciduous forest cover at low to mid-elevation locations over the landscape. Comparison of the model-predicted annual NPP to the plot-estimated values showed a significant correlation of R2 = 0.5. Stepwise addition of 30-meter resolution elevation data values explained no more than 20% of the residual variation in measured NPP patterns at BEF. Both the Landsat 7 and the 250-meter resolution MODIS derived mean annual NPP predictions for the BEF plot locations were within ± 2.5% of the mean of plot estimates for annual NPP.


Although MODIS imagery cannot capture the spatial details of NPP across the network of closely spaced plot locations as well as Landsat, the MODIS satellite data as inputs to the NASA-CASA model does accurately predict the average annual productivity of a site like the BEF.


The capacity of forests to sequester carbon from the increasing pool of atmospheric CO2 is becoming an issue of central importance for land managers and policy makers. Forested areas that consistently add carbon by growth in ecosystem production are potentially important as current and future sinks for industrial CO2 emissions. Conversely, land areas that do not consistently sequester carbon over time may be adding to already rising atmospheric CO2 levels from fossil fuel burning sources, Temperate forests recovering from disturbances, such as harvest for wood products or regrowth on abandoned agricultural lands, may represent important sinks globally for CO2 [1]. Nonetheless, affordable and rapid methods to understand and quantify the factors controlling the productivity of forest stands over most of the United States still await development.

Although regional- to global-scale relationships between forest carbon uptake and major climatic gradients have been demonstrated in several types of vegetation productivity models [2, 3], capturing landscape scale (e.g., less than 1 km2), patterns in forest carbon cycles has proven to be a challenge. At a landscape level, the influence of macroclimate is often less important than other sources of spatial variation, including disturbances, hydrology, and soil nutrient supply [46]. Nevertheless, several studies have addressed sub-landscape level variation in forest production through coupled applications of remote sensing and ecosystem process models [79].

Satellite remote sensing from instruments like Landsat has been applied to forest canopy studies for over a decade [10, 11]. The launch of NASA's Terra satellite platform in 1999 with the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument on-board initiated a new era in remote sensing of the Earth system with promising implications for forest carbon research. Direct input of satellite vegetation index "greenness" data from the MODIS sensor into ecosystem simulation models can now be used to estimate spatial variability in monthly net primary production (NPP), biomass accumulation, and litter fall inputs to soil carbon pools [12]. These global MODIS vegetation data sets are available at no charge from NASA data centers, which makes their application for carbon cycle studies most affordable, once proven to be scientifically robust.

Satellite data

Operational MODIS algorithms were the first to generate the Enhanced Vegetation Index (EVI) [13] as global image coverages from 2000-present. As a successor to the two-channel normalized difference vegetation index (NDVI), EVI represents an optimized vegetation index, whereby the vegetation index isolines in red and near infra-red spectral bands are designed to approximate vegetation biophysical isolines derived from canopy radiative transfer theory and/or measured biophysical and optical relationships. EVI was developed to optimize the greenness signal, or area-averaged canopy photosynthetic capacity, with improved sensitivity in high biomass regions and improved vegetation monitoring through a de-coupling of the canopy background signal and a reduction in atmosphere influences. Houborg and Soegaard [14] found MODIS EVI was able to accurately describe the variation in green biomass, in agriculture areas in Denmark, up to green LAI of 5 (R2 = 0.91). The EVI has been found useful in estimating absorbed PAR related to chlorophyll contents in vegetated canopies [15], and has been shown to be highly correlated with processes that depend on absorbed light, such as gross primary productivity (GPP) [16, 17].

In this study, we present the results of the NASA-CASA (Carnegie-Ames-Stanford Approach) model to predict net primary productivity (NPP) fluxes using both Landsat 7 and MODIS imagery as a means to infer variability in temperate forests of the eastern United States. Our NASA-CASA model [1820] has been designed to estimate monthly patterns in carbon fixation and plant biomass increments using moderate-to-high spatial resolution (30 m to 250 m) satellite remote sensing of surface vegetation characteristics and driven with spatially interpolated climate. The main objectives of this analysis were to (1) evaluate uncertainties in the use of Landsat 30-m imagery collected once during a summer growing season to predict NPP over a forested landscape and (2) to assess whether MODIS image data can capture landscape-level variability in forest production as well as Landsat image data in an eastern U. S. hardwood forest.

Site description

The Bartlett Experimental Forest – BEF (at approximately Latitude 44.0° N, Longitude 71.3° W) was established in 1932 as a U.S. Forest Service long-term research site. The 25 km2 BEF is located within the White Mountain National Forest, a heavily forested, mountainous region in north-central New Hampshire. The BEF is representative of the larger White Mountain region, having similar vegetative composition (Figure 1), soil types, disturbance histories and topography. As is true of the wider White Mountain region, all but a small portion of the BEF is believed to have a history of logging, although roughly 45% of the land area (particularly at higher elevations) has remained uncut since at least 1890. BEF is an actively managed forest and continues to be subjected to a variety of harvesting practices typical of the region. This is reflected in a range of successional sequences, forest patch sizes, and structural distributions. Natural disturbances affecting portions of the forest include a late 19th century fire, severe wind disturbance resulting from hurricanes in 1938 and 1954, and a damaging ice storm in 1998.

Figure 1
figure 1

National Land Cover Dataset (NLCD) cover types for the BEF (outline boundary shown).

In 1932, a network of permanent forest inventory plots (approximately 0.1 ha each) were established on a regular grid, spaced 200 m by 100 m apart [9]. All trees on most of these plots have been measured by 1-inch diameter classes in at least three time periods, the most recent complete re-measurement at the time of this study being in 2001–2003 [21]. Plot elevations ranged from 220 m to 731 m and represent a range in species composition and successional status. Major tree species encountered were sugar maple (Acer saccharum), American beech (Fagus grandifolia), white ash (Fraxinus americana), paper birch (Betula papyrifera), yellow birch (Betula alleghaniensis), red maple (Acer rubrum), pin cherry (Prunus pennsylvanica), eastern hemlock (Tsuga canadensis), red spruce (Picea rubens), balsam fir (Abies balsamea) and eastern white pine (Pinus strobus). Most plots contained mixtures of two or more species.

Ollinger and Smith [9] hypothesized that there could be a shift in the importance of moisture versus temperature limitations on forest production from low to high elevations within the White Mountains region. Increases in water availability with elevation may result from several factors in the northeast, including increased precipitation, a decrease in transpiration caused by lower temperatures, and the tendency for mid-elevation soils to be derived from deeper deposits of fine-textured glacial till. Supporting this interpretation, Federer [22] concluded that New Hampshire forests in low-elevation areas probably experience water deficits in most years.


A comparison of different modeling methods is presented in the following section. Each of these methods use remote sensing data inputs to characterize forest cover attributes that can influence NPP. The principal aim of this analysis is to determine how well satellite image data from Landsat and MODIS can capture landscape-level variability in forest production.

PnNET-II model for plot-based NPP

Above ground net primary productivity (ANPP, g C m-2 yr-1) for BEF plots has been estimated as the sum of wood production plus foliar production, based on plot measurements of tree diameters, litterfall, and allometeric scaling. Using these ANPP estimates for validation, a complete forest NPP (above- and below-ground) data set was generated from the forest process model PnET-II [23], initialized with high spectral resolution imagery from AVIRIS (Airborne Visible/Infrared Imaging Spectrometer) for foliar N estimation [9]. PnET-II requires a number of other input parameters summarizing vegetation and site characteristics, along with monthly climatic data. Vegetation parameters include foliar N, leaf retention time and growing-degree day variables describing the phenology of leaf production and senescence. Required climatic and soil inputs include temperature, precipitation, photosynthetically-active radiation (PAR), and soil water holding capacity (WHC). For pixel-by-pixel application at BEF, PnET-II was run in conjunction with image-derived foliar N estimates and a 20-m resolution digital elevation model (DEM). For each 20-m pixel estimate of NPP, geographic coordinates and elevation were extracted and used to estimate maximum and minimum temperature, vapor pressure, precipitation, and PAR.

CASA modeling methods for NPP

As documented in by Potter [19], the monthly NPP flux, defined as net fixation of CO2 by vegetation, is computed in NASA-CASA on the basis of light-use efficiency [24]. Monthly production of plant biomass is estimated as a product of time-varying surface solar irradiance, Sr, and EVI from the MODIS or Landsat satellite data, plus a constant light utilization efficiency term (emax) that is modified by time-varying stress scalar terms for temperature (T) and moisture (W) effects (Equation 1).

NPP = S r EVI e max T W

The default emax term is set uniformly at 0.39 g C MJ-1 PAR, a value that derives from calibration of predicted annual NPP to previous field estimates [2]. This model calibration has been validated globally by comparing predicted annual NPP to more than 1900 field measurements of NPP (Figure 2a). Interannual NPP fluxes from the CASA model have been reported [25] and validated against multi-year estimates of NPP from field stations and tree rings [26]. Our NASA-CASA model has been validated against field-based measurements of NPP fluxes and carbon pool sizes at multiple boreal forest sites in North America [18, 2628] and against atmospheric inverse model estimates of global carbon fluxes [20].

Figure 2
figure 2

Comparison of annual observed NPP to predicted values from the NASA-CASA model (driven by 0.5° data inputs). (a) Inputs from AVHRR-FPAR 1982 and climate means from New et al. (2000), (b) Inputs from MODIS-EVI 2001 and climate from NCEP reanalysis products for 2001. Both figures include the 1:1 regression line. The data set of more than 1900 observed NPP points was compiled for the Ecosystem Model-Data Intercomparison (EMDI) activity by the Global Primary Productivity Data Initiative (GPPDI) working groups of the International Geosphere Biosphere Program Data and Information System [34,35].

The T stress scalar is computed with reference to derivation of an optimal seasonal temperature (Topt) for plant production. Over large areas, the Topt setting will vary by latitude and longitude, ranging from near 0°C in the Arctic to the middle thirties in low latitude deserts. The W stress scalar is estimated from monthly water deficits, based on a comparison of moisture supply (precipitation) to potential evapotranspiration (PET) demand using the method of Thornthwaite and Mather [29]. The Moderate Resolution Imaging Spectroradiometer (MODIS) 1-km land cover map [30] was used to specify the predominant land cover class for the W term in each pixel as either forest or non-forest classes.

Whereas previous versions of the NASA-CASA model [2, 18] used NDVI bands to approximate FPAR, the current model version instead has been calibrated to use MODIS EVI datasets as direct inputs to Equation 1 above. In long-term (1982 to 2004) simulations, continuity between AVHRR and MODIS sensor data for inputs to NASA-CASA is an issue that must be addressed by recalibration of annual NPP results post 2000. NASA-CASA model predictions with 2001 monthly MODIS EVI inputs have been adjusted using the same set of field measurements of NPP shown in Figure 2a, to which the model was previously calibrated for a best linear fit to AVHRR inputs [20]. To best match predictions with previously measured NPP estimates at the global scale (Figure 2b), the model emax term for 2001 MODIS EVI inputs should be reset to 0.55 g C MJ-1 PAR, a value that is globally 42% higher than previously used in the model for AVHRR-driven NPP predictions from 1982–1998 [20]. The regression coefficient (with line intercept forced through zero) of R2 = 0.91 for this NPP recalibration to 2001 MODIS EVI inputs was statistically significant (p < 0.01).

In the following section, an evaluation procedure at the landscape scale is presented for the NASA-CASA model at BEF using growing season EVI computed from Landsat 7 Enhanced Thematic Mapper (ETM+) imagery at a 30-m gridded spatial resolution. We hypothesize that to match estimated NPP across BEF sample plots, the spatial detail provided by Landsat imagery collected during the summer period when leaf area is near maximum for the year will compensate for the lack of repeated EVI sampling during each month of the year (that, for example, is provided by the MODIS sensor, even though MODIS data is at a coarser 250-m minimum resolution). The rationale behind this hypothesis is that monthly climate inputs to the NASA-CASA model that uses a single Landsat image may be sufficient to capture the canopy phenology patterns throughout any given year that MODIS monthly satellite observations are used to capture. Furthermore, 30-m elevation and Landsat data inputs to NASA-CASA may serve as proxies for local weather variation effects on NPP.

In all the NASA-CASA simulation results presented for the BEF in this study, 4-km resolution spatial grids from PRISM (Parameter-elevation Regressions on Independent Slopes Model) [31] were used for precipitation, average maximum temperature, and average minimum temperature. These 4-km climatologies were derived from U. S. weather stations records interpolated first into 30 arc-second data sets. PRISM is unique in that it incorporates a spatial climate knowledge base that accounts for topographic influences such as rain shadows, temperature inversions, and coastal effects, in the climate mapping process. An example of the PRISM monthly temperature product for the BEF area is shown in Figure 3, wherein higher elevation areas in the western portion of the area typically have cooler summertime temperatures than the lower elevation portions of the area.

Figure 3
figure 3

PRISM [31] mean monthly temperature grid (20 km × 16 km) for the BEF area (outlined at center) showing the gradient of relatively cooler to warmer conditions moving from higher elevation to lower elevation, predominantly west to east (left to right).

CASA model NPP using Landsat EVI inputs

The Landsat 7 imagery used for this NASA-CASA modeling was collected on August 29, 2001. Predicted NPP from the NASA-CASA model for the BEF (Figure 4) shows variations related to both vegetation cover type (Figure 1) and elevational effects on mean air temperatures (Figure 3). Overall, the highest predicted NPP from the NASA-CASA model was for deciduous forest cover at low to mid-elevation locations over the landscape.

Figure 4
figure 4

CASA model NPP for the BEF area using Landsat 7 EVI and monthly PRISM climate inputs from 2001. Plot measurement sites are shown as a point grid.

The Landsat 30-m resolution grid coverage was used to stratify 3 × 3 pixel areas around each BEF plot location into comparable predictions of annual NPP from the CASA model. Comparison of the NASA-CASA model-predicted annual NPP to the BEF-estimated values (Figure 5) showed a significant correlation of R2 = 0.5 (p < 0.05) following the removal of 90 outlier plot values that were determined to have suffered heavy damage to forest stands during winter ice storms of January, 1998. The AVIRIS imagery used as an input to BEP plot NPP estimation was collected in 1997 and much of the Forest Service field data reflects pre-1997 tree status as well. Because we have used a more recent (i.e., 2001) Landsat 7 image for this study, NASA-CASA predicted NPP should be a more updated representation of BEF plot productivity.

Figure 5
figure 5

Plot measurement NPP for the BEF area versus CASA model NPP using Landsat 7 EVI and monthly PRISM climate inputs from 2001.

In an attempt to improve upon these initial results, we subdivided the BEF plots into three different forest types: Deciduous, Evergreen, and Mixed, defined by the National Land Cover Dataset (NLCD) [32], 1999; (Figure 1) at 30-m resolution. Separate correlation results for NASA-CASA predictions versus the measured plot NPP were generated for each NLCD forest type. Previously published studies for the BEF [9, 21] suggest that Deciduous stands have higher nitrogen in their tree leaves and hence higher average annual NPP even at the same leaf cover level as the Evergreen stands. We find that the simple breakout by NLCD forest type did not improve correlation statistic values, but did result in reduced variations around the mean and hence the newly stratified correlation line falls closer to a 1:1 line for the Evergreen forest stands (data not shown). This closer 1:1 match was not observed for the Deciduous stands, which presumably grow on a wider range of soil nutrient conditions.

Previously published studies also suggested that low elevation stands at BEF tend to have lower NPP than mid- to higher elevation stands [9], a trend that runs counter to most other forested locations with rugged terrain. Forest soils at low elevation at the BEF can have lower moisture availability compared to mid-elevation soils. We hypothesized therefore that including elevation as an independent variable in the regression analysis would explain a significant additional portion of the unexplained variance of estimated NPP in Figure 5. This hypothesis was tested using stepwise multiple linear regression, from which it was determined that addition of 30-meter elevation data values increased the R2 coefficient by 0.17. The addition of 30-meter estimates for aspect direction (0–360°) in the stepwise regression increased the R2 coefficient by another 0.02, leading to a final R2 = 0.69 for correlation of NASA-CASA predicted NPP plus topographic variables with the measured NPP patterns at BEF plots.

Because we can compute both NDVI and EVI from the 2001 Landsat 7 bands, it is possible to compare NASA-CASA NPP results to the BEF NPP estimates using either vegetation index from the satellite imagery. It is evident that NDVI is saturating at the high end of the greenness (i.e., actual leaf area index – LAI) range, whereas EVI is saturating to a lower degree in this manner and instead spreads out the greenness inputs to NASA-CASA more linearly over the entire range of NPP estimates at the BEF (Figure 6). The result is a closer fit to the predicted 1:1 line using the EVI compared to the NDVI as the NASA-CASA model input variable (with the same climate and topographic settings) over the BEF coverage area.

Figure 6
figure 6

Plot measurement NPP for the BEF area versus CASA model NPP using Landsat 7 EVI (blue diamond symbols) or NDVI (red square symbols) and monthly PRISM climate inputs from 2001. Best-fit regression lines are shown for both EVI and NDVI results.

CASA Model NPP Comparison using MODIS EVI

MODIS EVI composite images for all months of 2001 were used in place of the single Landsat 7 EVI imagery to predict annual NPP at the BEF. These continuous gridded observations of canopy greenness cover from the MODIS instrument have been collected to improve the tracking of canopy green-up and green-down over a growing season, particularly for Deciduous and Mixed forest stand types. MODIS time series profiles for EVI at three 250-m pixel locations within the BEF could be identified at locations of nearly uniform coverage (> 90%) of single forest classes from the NLCD map. These profiles showed that Deciduous and Mixed forest types reached maximum greenness cover levels by end of July, whereas Coniferous forest types increased in greenness cover levels until mid-August (Figure 7). Compared to the other two forest types, greenness cover in Coniferous forest types consistently remained about 20% lower from May through December and then increased to about 10% higher during January and February.

Figure 7
figure 7

MODIS time series profiles for EVI at three 250-m pixel locations selected as representative of nearly uniform coverage (> 90%) of Deciduous, Coniferous, and Mixed forest classes at the BEF.

Nevertheless, using fully gridded 250-m MODIS EVI inputs to the NASA-CASA model (with no sub-sampling bias), comparison of predicted annual NPP to the BEF-estimated plot NPP values showed no significant correlation result (at p < 0.1), regardless of whether the entire collection of plot estimates were considered together, or were separated into predominantly Deciduous, Coniferous, and Mixed forest classes (data not shown). This implies that the relatively coarse spatial resolution of the MODIS 250-m EVI imagery cannot capture landscape-level variability in tree production in the same way that 30-m Landsat image data can for this eastern hardwood forest location.

On the other hand, both the MODIS and Landsat 7 derived mean annual NPP predictions for the BEF plot locations were within ± 2.5% of the mean of PnET-based plot estimates for annual NPP. The mean annual NPP prediction from the NASA-CASA model using MODIS 250-m EVI imagery for the BEF was 503 g C m-2 yr-1, and therefore was slightly closer to the mean of all PnET-based plot measurements for annual NPP of 492 g C m-2 yr-1 than was the mean annual NPP prediction from the NASA-CASA model using Landsat 30-m EVI imagery of 479 g C m-2 yr-1. This finding suggests that, although MODIS 250-m EVI imagery cannot capture the spatial details of NPP across the network of closely spaced plot locations, it does accurately predict the average annual productivity of a site like the BEF.


Evidence from this modeling study indicates that the satellite-observed canopy greenness represented by the EVI is useful as a variable to help account for carbon sinks in northeastern forest ecosystems. When controlling for elevation, the NASA-CASA model using inputs of 30-m resolution Landsat 7 EVI for the month of peak growing season temperatures captures a substantial portion (> 67%) of landscape scale variation in the controls (such as foliar N content and moisture stress effects) on annual NPP at BEF. Nonetheless, a major limitation of Landsat imagery as a routine input to forest productivity models is the relatively low number of cloud-free scenes that can be collected over the course of a growing season in most temperate and humid climate zones. With a return schedule of approximately 16 days, the probability of collecting one cloud-free Landsat image per month of the year is low.

In comparison, MODIS imagery can be used in the NASA-CASA model to generate accurate predictions of the average forest NPP for a site the size of the BEF. Because MODIS provides global daily coverage of the entire Earth surface, these 250-meter resolution images can be composited over several weeks to provide an economical and rapid method to capture the regional-scale factors controlling seasonal productivity over most forested areas of the United States. NASA-CASA model predictions from MODIS 250-meter resolution data inputs are best suited to assessments of large forested tracts of land where stand ages are relatively uniform.

A next logical step in this methodology of using satellite image data as inputs to models of forest production at the landscape scale is to combine Landsat and MODIS data sets and enhance the advantages of both sensor measurements. We are investigating techniques to interpolate between dates of 30-meter resolution Landsat EVI data values by an appropriate sampling of monthly MODIS EVI values at selected locations that can represent much of the spatial variability in forest stand types and canopy properties detected in the higher resolution Landsat imagery.


  1. Dixon RK, Brown S, Houghton RA, Solomon AM, Trexler MC, Wisniewski J: Carbon pools and flux of global forest ecosystems. Science 1994, 263: 185–190. 10.1126/science.263.5144.185

    Article  CAS  Google Scholar 

  2. Potter CS, Randerson JT, Field CB, Matson PA, Vitousek PM, Mooney HA, Klooster SA: Terrestrial ecosystem production: A process model based on global satellite and surface data. Global Biogeochem, Cycles 1993, 7: 811–841.

    Article  Google Scholar 

  3. Melillo JM, McGuire AD, Kicklighter DW, Moore B III, Vorosmarty CJ, Schloss AL: Global climate change and terrestrial net primary production. Nature 1993, 363: 234–240. 10.1038/363234a0

    Article  CAS  Google Scholar 

  4. Pastor J, Aber JD, McClaugherty CA, Melillo JM: Aboveground production and N and P cycling along a nitrogen mineralization gradient on Blackhawk Island, Wisconsin. Ecology 1984, 65: 256–68. 10.2307/1939478

    Article  CAS  Google Scholar 

  5. Reich PB, Grigal DF, Aber JD, Gower ST: Nitrogen mineralization and productivity in 50 hardwood and conifer stands on diverse soils. Ecology 1997, 78: 335–47.

    Article  Google Scholar 

  6. Reed RA, Finley ME, Romme WH, Turner MG: Aboveground net primary production and leaf-area index in early postfire vegetation in Yellowstone National Park. Ecosystems 1999,2(1):88–94. 10.1007/s100219900061

    Article  Google Scholar 

  7. Ollinger SV, Aber JD, Federer CA: Estimating regional forest productivity and water yield using an ecosystem model linked to a GIS. Landscape Ecol 1998, 13: 323–34. 10.1023/A:1008004423783

    Article  Google Scholar 

  8. Kimball JS, Keyser AR, Running SW, Saatchi SS: Regional assessment of boreal forest productivity using an ecological process model and remote sensing parameter maps. Tree Physiol 2000,20(11):761–75.

    Article  Google Scholar 

  9. Ollinger SV, Smith ML: Net primary production and canopy nitrogen in a temperate forest landscape: An analysis using imaging spectroscopy, modeling, and field data. Ecosystems 2005, 8: 760–778. 10.1007/s10021-005-0079-5

    Article  CAS  Google Scholar 

  10. Cohen WB, Spies TA: Estimating structural attributes of Douglas Fir/Western Hemlock forest stands from Landsat and SPOT imagery. Remote Sensing of Environment 1992, 41: 1–17. 10.1016/0034-4257(92)90056-P

    Article  Google Scholar 

  11. Curran PJ, Dungan JL, Gholz HL: Seasonal LAI in slash pine estimated with Landsat TM. Remote Sensing of Environment 1992, 39: 3–13. 10.1016/0034-4257(92)90136-8

    Article  Google Scholar 

  12. Potter C, Klooster S, Genovese V, Myneni R: Satellite data helps predict terrestrial carbon sinks. Eos – Transactions of the American Geophysical Union 2003, 84: 502–508.

    Article  Google Scholar 

  13. Huete A, Didan K, Miura T, Rodriguez EP, Gao X, Ferreira LG: Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sensing of Environment 2002, 83: 195–213. 10.1016/S0034-4257(02)00096-2

    Article  Google Scholar 

  14. Houborg RM, Soegaard H: Regional simulation of ecosystem CO 2 and water vapor exchange for agricultural land using NOAA AVHRR and Terra MODIS satellite data. Application to Zealand, Denmark. Remote Sensing of Environment 2004, 93: 150–167. 10.1016/j.rse.2004.07.001

    Article  Google Scholar 

  15. Zhang Q, Xiao X, Braswell B, Linder E, Baret F, Moore B III: Estimating light absorption by chlorophyll, leaf and canopy in a deciduous broadleaf forest using MODIS data and a radiative transfer model. Remote Sensing of Environment 2005, 99: 357–371. 10.1016/j.rse.2005.09.009

    Article  Google Scholar 

  16. Xiao X, Zhang Q, Braswell B, Urbanski S, Boles S, Wofsy S, Moore B III, Ojima D: Modeling gross primary production of temperate deciduous broadleaf forest using satellite images and climate data. Remote Sensing of Environment 2004, 91: 256–270. 10.1016/j.rse.2004.03.010

    Article  Google Scholar 

  17. Rahman AF, Sims DA, Cordova VD, El-Masri BZ: Potential of MODIS EVI and surface temperature for directly estimating per-pixel ecosystem C fluxes. Geophys Res Lett 2005, 32: L19404. 10.1029/2005GL024127

    Article  Google Scholar 

  18. Potter CS, Klooster SA, Brooks V: Interannual variability in terrestrial net primary production: Exploration of trends and controls on regional to global scales. Ecosystems 1999, 2: 36–48. 10.1007/s100219900056

    Article  Google Scholar 

  19. Potter CS: Terrestrial biomass and the effects of deforestation on the global carbon cycle. BioScience 1999, 49: 769–778. 10.2307/1313568

    Article  Google Scholar 

  20. Potter C, Klooster S, Myneni R, Genovese V, Tan P, Kumar V: Continental scale comparisons of terrestrial carbon sinks estimated from satellite data and ecosystem modeling 1982–98. Global and Planetary Change 2003, 39: 201–213. 10.1016/j.gloplacha.2003.07.001

    Article  Google Scholar 

  21. Smith ML, Ollinger SV, Martin ME, Aber JD, Hallett RA, Goodale CL: Direct estimation of aboveground forest productivity through hyperspectral remote sensing of canopy nitrogen. Ecological Applications 2002, 12: 1286–1302. 10.1890/1051-0761(2002)012[1286:DEOAFP]2.0.CO;2

    Article  Google Scholar 

  22. Federer CA: Frequency and intensity of drought in New Hampshire forests: evaluation by the BROOK model. In Applied Modeling in Catchment Hydrology. Water Resour Pub, Littleton, CO; 1982:459–470.

    Google Scholar 

  23. Aber JD, Ollinger SV, Feder CA, Reich PB, Goulden ML, Kicklighter DW, Melillo JM, Lathrop RG Jr: Predicting the effects of climate change on water yield and forest production in Northeastern U.S. Climate Research 1995, 5: 207–222. 10.3354/cr005207

    Article  Google Scholar 

  24. Monteith JL: Solar radiation and productivity in tropical ecosystems. Journal of Applied Ecology 1972, 9: 747–766. 10.2307/2401901

    Article  Google Scholar 

  25. Behrenfeld MJ, Randerson JT, McClain CR, Feldma GC, Los SQ, Tucker CI, Falkowski PG, Field CB, Frouin R, Esaias WE, Kolber DD, Pollack NH: Biospheric primary production during an ENSO transition. Science 2001, 291: 2594–2597. 10.1126/science.1055071

    Article  CAS  Google Scholar 

  26. Malmström CM, Thompson MV, Juday GP, Los SO, Randerson JT, Field CB: Interannual variation in global scale net primary production: Testing model estimates. Global Biogeochem Cycles 1997, 11: 367–392. 10.1029/97GB01419

    Article  Google Scholar 

  27. Amthor JS, Chen JM, Clein JS, Frolking SE, Goulden ML, Grant RF, Kimball JS, King AW, McGuire AD, Nikolov NT, Potter CS, Wang S, Wofsy SC: Boreal forest CO 2 exchange and evapotranspiration predicted by nine ecosystem process models: Inter-model comparisons and relations to field measurements. J Geophys Res 2001, 106: 33,623–33,648. 10.1029/2000JD900850

    Article  CAS  Google Scholar 

  28. Hicke JA, Asner GP, Randerson JT, Tucker CJ, Los SO, Birdsey R, Jenkins JC, Field CB, Holland EA: Satellite-derived increases in net primary productivity across North America, 1982–1998. Geophys Res Lett 2002,29(10):1427. doi:10.1029/2001GL013578 doi:10.1029/2001GL013578 10.1029/2001GL013578

    Article  Google Scholar 

  29. Thornthwaite CW, Mather JR: The water balance. In Publications in Climatology. Drexel Institute of Technology. Volume VIII. Centerton, New Jersey; 1955:86.

    Google Scholar 

  30. Friedl MA, McIver DK, Hodges JCF, Zhang XY, Muchoney D, Strahler AH, Woodcock CE, Gopal S, Schneider A, Cooper A, Baccini A, Gao F, Schaaf C: Global land cover mapping from MODIS: algorithms and early results. Remote Sensing of Environment 2002, 83: 287–302. 10.1016/S0034-4257(02)00078-0

    Article  Google Scholar 

  31. Daly C, Gibson WP, Doggett , Smith MJ, Taylor G: Up-to-date monthly climate maps for the conterminous United States. In Proc., 14th AMS Conf. on Applied Climatology, 84th AMS Annual Meeting Combined Preprints. Amer. Meteorological Soc., Seattle, WA, January 13–16, Paper P5.1, CD-ROM; 2004.

    Google Scholar 

  32. U.S. Geological Survey (USGS) and U.S. Environmental Protection Agency (USEPA): National Land Cover Data (NLCD).1999. []

    Google Scholar 

  33. NASA CQUEST website[]

  34. Barrett DJ: NPP Multi-Biome: VAST Calibration Data, 1965–1998. Data set Oak Ridge National Laboratory Distributed Active Archive Center, Oak Ridge, Tennessee, U.S.A; 2001. []

    Google Scholar 

  35. Zheng D, Prince S, Wright R: Terrestrial net primary production estimates for 0.5° grid cells from field observations – A contribution to global biogeochemical modeling. Global Change Biology 2003, 9: 46–64. 10.1046/j.1365-2486.2003.00534.x

    Article  Google Scholar 

Download references


This work was supported by grants from NASA programs in Earth Observing System (EOS) Interdisciplinary Science and Applied Sciences. NASA-CASA model data sets used in this report are available online [33], as part of the Carbon Query and Evaluation Support Tools (CQUEST) project. Portions of this research are based upon data generated from long-term studies at the Bartlett Experimental Forest, funded by the U.S. Forest Service, Northern Research Station. We thank Scott Ollinger and Richard Birdsey for valuable comments on this and related carbon modeling studies.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Christopher Potter.

Additional information

Competing interests

The author(s) declare that they have no competing interests.

Authors' contributions

All authors have made substantial contributions to the analysis and interpretation of data, have been involved in drafting the manuscript and revising it critically for important intellectual content, and have given final approval of the version to be published.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Potter, C., Gross, P., Genovese, V. et al. Net primary productivity of forest stands in New Hampshire estimated from Landsat and MODIS satellite data. Carbon Balance Manage 2, 9 (2007).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: