Extrapolating biosphere-atmosphere CO2 flux observations to larger scales in space, part of the so-called “upscaling” problem, is a central challenge for surface-atmosphere exchange research. Upscaling CO2 flux in tundra is complicated by the pronounced spatial variability of vegetation cover. We demonstrate that a simple model based on chamber observations with a pan-Arctic parameterization accurately describes up to 75% of the observed temporal variability of eddy covariance—measured net ecosystem exchange (NEE) during the growing season in an Abisko, Sweden, subarctic tundra ecosystem, and differed from NEE observations by less than 4% for the month of June. These results contrast with previous studies that found a 60% discrepancy between upscaled chamber and eddy covariance NEE sums. Sampling an aircraft-measured normalized difference vegetation index (NDVI) map for leaf area index (L) estimates using a dynamic flux footprint model explained less of the variability of NEE across the late June to mid-September period, but resulted in a lower root mean squared error and better replicated large flux events. Findings suggest that ecosystem structure via L is a critical input for modeling CO2 flux in tundra during the growing season. Future research should focus on quantifying microclimate, namely photosynthetically active radiation and air temperature, as well as ecosystem structure via L, to accurately model growing season tundra CO2 flux at chamber and plot scales.
Introduction
Tundra ecosystems are characterized by pronounced spatial variability in vegetation structure and have the highest coefficient of variation of leaf area index (L) of any global biome (Asner et al., 2003). This spatial heterogeneity complicates efforts to upscale CO2 flux observations from chamber to eddy covariance flux tower to region (Oechel et al., 1998; Soegaard et al., 2000; Fox et al., 2008). Some studies in arctic and subarctic tundra demonstrate excellent agreement between CO2 flux measurements and model estimates at multiple spatial scales (Soegaard et al., 2000), but other efforts find relatively poor agreement (Oechel et al., 1998; Fox et al., 2008), due in part to the challenges posed by the spatial heterogeneity of vegetation structure and thereby ecosystem function. Robust approaches for scaling CO2 flux in tundra must be developed to make progress in our understanding of the terrestrial C cycle and its responses to a changing climate.
Multiple techniques are used to measure the flux of CO2 between biosphere and atmosphere, and estimates do not always agree, especially when extrapolated to different scales in space. Many studies find substantially different estimates of the net ecosystem exchange of carbon (NEE) and/or its components, gross ecosystem productivity (GEP), and ecosystem respiration (ER), when examining chamber and eddy covariance—measured fluxes (Janssens et al., 2001), whereas other studies find good agreement, at least for selected periods (Dore et al., 2003). A solution to the so-called upscaling problem from chamber to tower is required to make chamber-based measurements relevant at larger spatial scales and to link eddy covariance (EC) measurements to the leaf, stem, root, and soil mechanisms that ultimately give rise to the observed ecosystem-scale flux.
Chamber and EC flux estimates may converge or diverge for different reasons, including (1) flux overestimation with chamber-based techniques, especially in windy conditions with chambers that do not account for pressure pumping via the Venturi effect (Bain et al., 2005; Xu et al., 2006); (2) abscission of roots and mycorrhizae when placing chambers in the soil (Subke et al., 2009); (3) bias in the EC measurements due to advection, lack of energy balance closure, or diurnal footprint variability (Grace et al., 1996; Wilson et al., 2002; Oren et al., 2006); (4) uncertainty in gapfilling missing flux measurements for the calculation of seasonal sums (Goulden et al., 1996; Moncrieff et al., 1996); (5) choosing statistically unrepresentative locations within the flux footprint with which to make chamber-based measurements (Fox et al., 2008); (6) uncertainty as to how to partition vegetation within the flux footprint into functional units for modeling, including uncertainty in the flux footprint model itself (Schmid, 2002; Kljun et al., 2004); and (7) complications owing to the spatial aggregation of what is often a nonlinear relationship between surface-atmosphere exchange and micrometeorological drivers (Rastetter et al., 1992; Kustas and Norman, 2000; Hong and Kim, 2008; Stoy et al., 2009b; Hill et al., 2011). These issues can be partially resolved by studying situations where bias terms are likely to be minor and where spatial variability in ecosystem structure and function within the eddy covariance flux footprint has been characterized.
We argue that our tundra study ecosystem near Abisko, Sweden, represents an ideal test case for CO2 flux upscaling. Upscaled NEE estimates using ecosystem models based on a unique mobile chamber designed to minimize issues 1 and 2 (see Williams et al., 2006) can be compared against eddy covariance NEE measurements for periods that minimize issues 3 and 4. The unvented closed chamber system does not disturb the soil environment to account for issues 1 and 2, and the site is characteristically windy (averaging nearly 3.3 m s-1) with a short tower height to minimize uncertainties owing to CO2 storage in the air space below the eddy covariance system (Yang et al., 2007). The implications of the ecosystem model used here, that tundra ecosystem CO2 flux follows a predictable relationship with L (Shaver et al., 2007), simplifies issues 5 and 6 if spatial heterogeneity within the flux footprint is characterized, as has been addressed in our case using high-resolution airborne remote-sensing observations (Stoy et al., 2009b). This leaves only issue 7, which has been explored at the study site using an information theoretic framework to avoid biases due to nonlinearities (Stoy et al., 2009c).
Despite the proposed advantages for scaling CO2 flux at the Abisko tundra site, Fox et al. (2008) found that upscaled 1 m2 chamber flux measurements and EC tower measurements differed by 60% over a 40-day period during the peak growing season. Fox et al. (2008) summed chamber-based measurements of different vegetation classes to the scale of the turbulent eddy covariance flux footprint (Schmid, 1994, 1997) based on a detailed vegetation map, and used error estimates from both measurement techniques to demonstrate that the flux sums did not agree. They suggested that either error in the tower measurements (i.e. due to 3 and/or 4 above), or biased sampling using chambers (5 and/or 6) as opposed to chamber errors (1 and 2) were responsible for the lack of agreement.
A major assumption of the Fox et al. (2008) study was that detailed spatial information on community composition, regardless of their structural properties like L, is required to upscale chamber-based flux measurements to larger spatial scales. However, the implications of chamber studies from the Abisko tundra site and other arctic ecosystem types (Shaver et al., 2007; Street et al., 2007) suggest that only L, air temperature (Ta), and photosynthetic photon flux density (PPFD) are required to model CO2 flux in tundra during the growing season; community composition can be effectively excluded from the model. Specifically, Shaver et al. (2007) parameterized a simple model with such chamber measurements to argue that tundra ecosystems demonstrate functional convergence, meaning that the parameters of simple models for photosynthesis and respiration were similar across vegetation types and regions. Measured CO2 flux followed a predictable relationship with L (Williams et al., 2006; Street et al., 2007) regardless of species type or community composition. If these results are applicable to larger spatial scales than the 1 × 1 m chambers employed by previous studies (Shaver et al., 2007; Street et al., 2007), namely the eddy covariance flux footprint (Loranty et al., 2011), then NEE at the footprint scale should be accurately predicted with an unbiased estimate of L and accurate meteorological inputs.
We hypothesize that a model of ecosystem CO2 flux based on previous chamber measurements in tundra will match eddy covariance NEE measurements during growing season periods for which the model was parameterized (Williams and Rastetter, 1999; Shaver et al., 2007). We further test if footprint modeling coupled with spatially explicit normalized difference vegetation index (NDVI) maps or simple tower-based NDVI observations, as proxies for L, are required for accurate chamber-to-tower upscaling. We discuss the consequences of our findings for upscaling surface CO2 fluxes in tundra.
Methods
SITE DESCRIPTION
The Abisko tundra site (hereafter AT) is 6.2 km SSE of Abisko, Sweden, in mixed tundra vegetation. An eddy covariance tower was located on a moraine plateau at UTM coordinates 411180.486 E, 7577776.226 N, (zone 34W) at an elevation of 752.3 m a.s.l. A meteorological tower was located nearby at 411191.840 E, 7577785.175 N, 751.9 m a.s.l. Sensor location was measured by real-time kinetic GPS using a Trimble 5800 rover with Trimble 5700 base station (Trimble Navigation Ltd., Sunnyvale, California, U.S.A.). Six types of tundra vegetation were previously classified at AT: fell field, open Empetrum heath, closed Empetrum. heath, poor fen, shrub tundra, and grey willow scrub (mostly Salix glauca, noting frequent hybridization) (Shaver et al., 2007; Street et al., 2007; Fox et al., 2008). The same assemblages existed during the study period in 2007.
EDDY COVARIANCE
NEE was measured using an open path eddy covariance system comprised of a LI-COR 7500 infrared gas analyzer (Li-Cor, Lincoln, Nebraska, U.S.A.) and a Gill R3 Sonic Anemometer (Gill Instruments Limited, Lymington, Hampshire, U.K.). Sensors were mounted at 3.0 m on a telescopic mast. Half-hourly fluxes were calculated using FluxView (Centre for Ecology and Hydrology, Wallingford, U.K.). Processing involved correction for the angle of attack (Gash and Dolman, 2003), adjustment for lag times between sonic anemometer and infrared gas analyzer, coordinate rotations to align the horizontal wind vector along the half-hourly mean wind direction and to force the vertical component to zero, correcting observed heat flux for humidity effects (Schotanus et al., 1983), correction for spectral flux losses, the Webb-Pearman-Leuning correction for the effects of air density fluctuations on open path eddy covariance measurements (Webb et al., 1980), and the Burba cor- rection to account for heating of the open-path infrared gas analyzer (Burba et al., 2008).
METEOROLOGICAL MEASUREMENTS
Photosynthetic photon flux density (PPFD) was measured using a SKP 215 Quantum Sensor (Skye Instruments, Llandrindod Wells, U.K.) at a height of 2 m on the meteorological tower. Missing measurements were replaced using the linear relationship with incident shortwave radiation measured by a co-located CM11 Pyranometer (Kipp & Zonen, Delft, Netherlands). Air temperature was measured using a HMP45 sensor (Vaisala, Helsinki, Finland), also at a height of 2 m.
TOWER AND AIRCRAFT NDVI MEASUREMENTS
NDVI was measured from the micrometeorological tower using 20-min averaged output of SKR1800 two-channel reflectance sensors (Skye Instruments) pointing both upward and downward and mounted at 1.5 m on the meteorological tower. The SKR1800 as configured has a field of view of 60° and a footprint of 1.7 m in diameter. NDVI was converted to the formulation used by van Wijk and Williams (2005) and Williams et al. (2008) as discussed in Stoy et al. (2009b). L was estimated using the NDVI conversion for sensors at 1.5 m height following Williams et al. (2008):
Representative daily NDVI values were chosen by selecting periods with solar zenith angles between 55° and 75° to ensure that different days of observation were comparable given the low solar zenith angles experienced in Abisko before mid-June and after early September.
Aircraft-based NDVI measurements were used for the spatially explicit flux footprint analysis. NDVI was calculated for 4 m pixels in a 2 km × 2 km area including the AT tower footprint using observations from an Azimuth Systems AZ-16 Airborne Thematic Mapper (ATM) flown on 17 July 2005 (day of year [DOY] 198) as described in Stoy et al. (2009b).
The tower-based and aircraft-based NDVI measurements were compared against surface reflectance measurements made every 20 m in 200 m transects in the four cardinal directions from the meteorological tower during the peak of the growing season (DOY 185; 2007) using a FieldSpec2 spectroradiometer (ASD Inc., Boulder, Colorado, U.S.A.), henceforth called an “ASD.” The present analysis concentrates on the 200 × 200 m area centered around the AT tower that corresponds to the flux source area during periods with adequate convective and mechanical production of turbulent kinetic energy for accepting the half-hourly eddy covariance flux values (Aubinet et al., 2000). A further check on NDVI surrounding AT was provided by ASD measurements from a University of Edinburgh G-GEOS aircraft flyover on DOY 198, 2008.
NET ECOSYSTEM EXCHANGE MODEL
NEE was modeled using the observed relationship between photosynthesis, L, and irradiation (i.e. PPFD) and between respiration, air temperature, and L (“PLIRTLE”; Shaver et al., 2007). The gross primary productivity (GPP) submodel follows the aggregated canopy photosynthesis model of Rastetter et al. (1992). GPP, in µmol CO2 m-2 s-1, is assumed to follow a saturating response to PPFD, in µmol photons m-2 s-1, and is integrated through the canopy using the Beer-Lambert law for light attenuation:
where k is the Beer's law light extinction coefficient, assumed here to be 0.5 after Shaver et al. (2007), Eo?? (µmol CO2 µmol photons-1) is the initial slope of the light sensitivity response of photosynthesis, and Pmax is the light saturated rate of photosynthesis per unit leaf area (µmol m-2 s-1).
Shaver et al. (2007) found that ER models for tundra that include two sources, one sensitive to L and Ta and one not (the latter presumably from deeper soil horizons), fare better than models that assume a single substrate pool. We use model “ER2” (Shaver et al., 2007):
where R0 (µmol m-2 s-1) is base respiration at 0 °C, Rx (µmol m-2 s-1) is the component of ecosystem respiration that is insensitive to observed L and Ta, and β (°C-1) is the temperature sensitivity of R0. The global parameter set using observations from Sweden and Alaska as described in Shaver et al. (2007) was chosen for this analysis to investigate the proposed pan-Arctic applicability of PLIRTLE and the generality of the chamber-based model upscaling approach (Table 1). The net CO2 flux from atmosphere to biosphere is denoted as negative following the micrometeorological convention.
FLUX FOOTPRINT MODEL
A two-dimensional semi-analytic footprint model after Hsieh et al. (2000), extended to two dimensions by Detto et al. (2006) (see also Novick et al., 2004; Stoy et al., 2006; Oishi et al., 2008) was used to sample NDVI pixels from the ATM-generated NDVI map (see Fig. 3 later herein). AU NDVI pixels that contributed less than 0.01% to the footprint were excluded for efficient computation, and the total remaining footprint was rescaled to equal one. These pixels were weighted by flux footprint source area to derive a probability distribution of flux-weighted NDVI values, which were converted to L following Equation (1). This derived distribution of L entered the PLIRTLE model prior to averaging following Stoy et al. (2009b). Modeled NEE based on L estimates from the tower-mounted NDVI sensor (NEEmod,tower) and modeled NEE based on L estimates from the source-area weighted footprint analysis (NEEmod,footprint) were compared against eddy covariance—measured NEE. We focused the model comparison on the 83-day period between DOY 174 when continuous eddy covariance measurements began and DOY 257 when tower-mounted NDVI measurements became unreliable. We additionally highlight a period between DOY 186 and DOY 200 when PPFD and Ta exhibited characteristic mid-summer variability but tower-mounted NDVI measurements were to a first order constant.
TABLE 1
The pan-Arctic parameter set for the PLIRTLE model with ecosystem respiration model (ER) 2 and a fixed light extinction coefficient (k) of 0.5 following Shaver et al. (2007). Units for Eo are defined in the text.
![t01_275.gif](ContentImages/Journals/aare/45/2/1938-4246-45.2.275/graphic/WebImages/t01_275.gif)
FIGURE 1.
(A) Air temperature (Ta), and (B) photosynthetically active photon flux density (PPFD) between day of year (DOY) 120 and 250 of 2007, with a focus on the model comparison period between DOY 186 and DOY 200 (C: Ta; D: PPFD).
![f01_275.jpg](ContentImages/Journals/aare/45/2/1938-4246-45.2.275/graphic/WebImages/f01_275.jpg)
We note that chamber observations and subsequent modeling work were made at the same site during previous years (Shaver et al., 2007; Street et al., 2007), and that the ATM flyover also took place during previous years (Stoy et al., 2009b). If this previous measurement and modeling work can be used as a basis for upscaling NEE during subsequent years, it may be argued that the ability to upscale models based on chamber measurements to the eddy covariance flux footprint is further simplified.
Results
MICROCLIMATE
The summer of 2007 was characterized by brief warm spells extending from multiple days to over one week (Fig. 1). The highlighted period (DOY 186–200) was characterized by variable PPFD and Ta, which was generally >8 °C.
TOWER AND AIRCRAFT-MEASURED NDVI
Mean tower-measured daily NDVI was ∼0.55 preceding the period of rapid vegetation growth, which began around DOY 150. Mean daily NDVI reached its peak of ∼0.75 around DOY 180 and began a slow decline around DOY 210 (Fig. 2, part A). The mean daily L estimates that resulted from the NDVI observations following Equation (1) ranged from <0.2 m2 m-2 to ∼1 m2 m-2 (Fig. 2, part B).
Mean NDVI from the ASD-sampled transects within the flux footprint (0.683) and NDVI measured by the tower-mounted sensor (0.675) were not significantly different (two-sided t-test, p < 0.05) on DOY 185, 2007, when the ASD measurements were made. Mean tower-measured NDVI during DOY 186–200 was 0.73 and the mean weighted NDVI captured by the dynamic flux footprint (e.g. Fig. 3) over the same period was 0.71. A representative distribution of NDVI and L within the flux footprint for a half-hourly period is shown in Figure 4. Mean ATM-measured NDVI from the Airborne Remote Sensing Facility flyover on 17 July (DOY 198) 2005 in the vicinity of the tower representing the flux footprint was 0.72. Mean ASD-measured NDVI from the G-GEOS flyover in the vicinity of the tower was 0.73 on DOY 198, 2008.
MODEL-MEASUREMENT COMPARISON
Across the observation period, NEEmod,tower explained 61% of the variance of NEE measurements with a modeling efficiency (Loague and Green, 1991; Meyer and Butler, 1993) of 0.58 and root mean squared error (RMSE) of 15.7 µmol CO2 m-2 s-1. Mean NEEmod,tower was -1.74 µmol CO2 m-2 s-1, which is less than the mean measured NEE of -1.27 µmol CO2 m-2 s-1, indicating that the model overestimated ecosystem C uptake by 37%. NEEmod,footprint explained only 53% of the variability of NEE and had a lower modeling efficiency (0.53), but also had a far lower RMSE (1.3 µmol CO2 m-2 s-1). Mean NEEmod,footprint differed from observed NEE by less than 3%, far within the range of uncer-tainty on the order of 10% or more commonly attributed to eddy covariance observations (Moncrieff et al., 1996).
The differences in goodness-of-fit between NEEmod,footprint and NEEmod,tower extended beyond single-day examples; the RMSE of NEEmod,tower was generally lower during the earlier part of the measurement period before DOY 210, and the RMSE of NEEmod, footprint was generally lower during the latter part of the measurement period (Fig. 6). Both approaches replicated observed variability poorly during a period centered around DOY 203. This period was marked by relatively low retrieval of eddy covariance flux data that averaged less than 50% of available measurements per day owing in part to large rain events. Approximately 8.5 mm of rain was measured at the Abisko Research Station during the week beginning DOY 200. The RMSE of NEEmod,footprint tended to be lower during the latter part of the observation period after DOY 210. This period corresponded to larger footprint volumes; the peak of the footprint source weight function averaged 41m before DOY 210 and was over 60% longer (66 m) after DOY 210, on average. The increase in footprint dimension corresponded to a decrease in the sensible heat flux toward the end of the measurement period (r = -0.76; p < 0.001); wind speed was not statistically related to DOY (p = 0.59). Both modeling approaches explained over 60% of the variability of measured NEE when excluding the period between DOY 200 and DOY 204 during the period of low NEE measurement retrieval, and after DOY 240 when L estimates were frequently lower than 0.4 m2 m-2 and vegetation was senescing (Fig. 2, part B).
FIGURE 2.
(A) Daily measurements of the normalized difference vegetation index (NDVI) using a tower-mounted sensor, and (B) subsequent estimates of the leaf area index (L) at the Abisko Tundra flux tower during the 2007 growing season. Circles represent the daily mean of acceptable NDVI measurements and subsequent L estimates. Error bars represent 1 standard deviation about the daily mean.
![f02_275.jpg](ContentImages/Journals/aare/45/2/1938-4246-45.2.275/graphic/WebImages/f02_275.jpg)
When focusing on the highlighted period between DOY 186 and DOY 200, NEEmod,footprint matched the diurnal patterns of observed NEE better during certain days (e.g. DOY 186 and DOY 196; Fig. 5, part A) and tended to better replicate periods with large flux magnitude (Fig. 5, part C). NEEmod,tower better matched the diurnal pattern or observed NEE on other days (e.g. DOY 195 and DOY 197). Both modeling approaches frequently missed large positive and negative excursions of observed NEE (e.g. DOY 190). The local Abisko-area PLIRTLE parameterizations (Shaver et al., 2007) did not improve model fit versus the pan-Arctic parameterization (data not shown).
FIGURE 3.
A normalized difference vegetation index (NDVI) map at 4 m resolution of the area around the Abisko Tundra eddy covariance tower (AT), overlaid with a flux footprint for a half-hourly flux measurement at 10:00 (local time) on DOY 186, 2007. The outermost ellipsoid represents the 80% footprint, and the innermost ellipsoid represents the 20% footprint.
![f03_275.jpg](ContentImages/Journals/aare/45/2/1938-4246-45.2.275/graphic/WebImages/f03_275.jpg)
FIGURE 4.
(A) The probability distribution function of tower-measured normalized difference vegetation index (NDVI) at the Abisko, Sweden, tundra research site with a kernel density estimate of footprintweighted NDVI for the time period shown in Figure 3 (10:00 on DOY 186,2007). (B) The leaf area index (L) distributions that result from converting the NDVI distributions in part A using the transfer functions described in Williams et al. (2008) (see Equation 1).
![f04_275.jpg](ContentImages/Journals/aare/45/2/1938-4246-45.2.275/graphic/WebImages/f04_275.jpg)
FIGURE 5.
(A) The measured net ecosystem exchange of CO2 (NEE) at a tundra ecosystem near Abisko, Sweden, with PLIRTLE (Shaver et al., 2007) model outputs that employ leaf area index (L) estimated from tower-based normalized difference vegetation index (NDVI) measurements (NEEmod,tower, B) and NDVI sampled from a dynamic flux footprint (NEEmod,footprint, C).
![f05_275.jpg](ContentImages/Journals/aare/45/2/1938-4246-45.2.275/graphic/WebImages/f05_275.jpg)
FIGURE 6.
Goodness-of-fit statistics for the PLIRTLE model with leaf area index estimated by a tower-mounted normalized difference vegetation index (NDVI) sensor (“Tower”) and with leaf area index estimated using NDVI from a flux footprint analysis following Hsieh et al. (2000) and Detto et al. (2006) using a 4 m NDVI map (“Footprint”, see Fig. 3). (A) EF: modeling efficiency; (B) r 2: coefficient of determination; (C) RMSE: root mean squared error.
![f06_275.jpg](ContentImages/Journals/aare/45/2/1938-4246-45.2.275/graphic/WebImages/f06_275.jpg)
Discussion
UPSCALING CO2FLUX IN TUNDRA
A model (PLIRTLE) based on microclimate and L (Shaver et al., 2007) that was parameterized using pan-Arctic chamber-based flux measurements (Shaver et al., 2007; Street et al., 2007) explained up to 75% of the variance of NEE for select two-week periods and explained over 60% of the variance of NEE across most of the growing season. The 4% difference in magnitude between NEE and NEEmod,footprint is well within NEE error from eddy covariance observations, which is usually considered to be on the order of 10–15% (Goulden et al., 1996; Moncrieff et al., 1996; Stoy et al., 2006) depending on flux magnitude (Hollinger and Richardson, 2005; Richardson et al., 2006). The maximum percent of variance explained by the models is similar to but slightly less than the explanatory power of the PLIRTLE model based on chamber measurements (80% ± 5%) in pan-Arctic tundra ecosystems (Shaver et al., 2007).
Our results contrast with those of Fox et al. (2008), who found that upscaling the chamber measurements of Shaver et al. (2007) and Street et al. (2007) based on a map of dominant species gave NEE values 60% greater than eddy covariance-based measurements, and suggested that the chamber sampling design was not systematic. This interpretation is correct; the sampling design of Shaver et al. (2007) and Street et al. (2007) was meant to capture representative examples of dominant tundra ecosystem types for the purposes of a pan-Arctic comparison, not to capture characteristic rates of ecosystem CO2 exchange at the eddy covariance footprint scale. In other words, assuming that vegetation type controlled surface fluxes, rather than vegetation characteristics (here via L), resulted in erroneous upscaled fluxes. We recommend upscaling biosphere-atmosphere flux observations through validated models of ecosystem function (Rastetter et al., 1992) rather than species composition.
FOOTPRINT MODELING
Adding a dynamic flux footprint via NEEmod,footprint decreased the explanatory power of the model for many periods. NEEmod,foot print represented an improvement over NEEmodtower in certain instances, including large CO2 uptake events (Fig. 5). On occasion, NEEmod,footprint estimated greater CO2 sequestration than measurements supported. We suggest that this overestimate is likely due to the frequent presence of high NDVI outliers (Figs. 3 and 4) when scaling using the flux footprint. These high NDVI patches correspond largely to Salix shrubs, which were also frequent outliers when parameterizing the PLIRTLE model (Shaver et al., 2007). Further modeling work with chamber observations suggested that the complex light environment near the edges of Salix vegetation patches caused GPP to deviate from a simple, linear relationship with L (Fletcher et al., 2012). NEEmod,footprint often included NDVI pixels that corresponded to the Salix patches to which PLIRTLE does not fit as well. Future research should focus on modeling NEE in Salix, especially if shrub growth and encroachment across the circumpolar tundra biome continue (Forbes et al., 2010; Macias-Fauria et al., 2012).
Our findings suggest that adding complexity via the footprint model, and thereby uncertainty, to an upscaling approach that matches observations reasonably well does not necessarily improve model fit. NEEmod,footprint did not represent an unambiguous improvement over NEEmod,tower (Fig. 6) except during the latter part of the measurement period when the flux footprint tended to be larger. It is unclear if disturbance near the tower or the larger population of NDVI pixels sampled by the larger footprint can be attributed to the improved fit of NEEmod,footprint during the latter part of the growing season. Regardless, fortuitous tower-mounted NDVI sensor placement in an area with an average NDVI that closely matched the NDVI in the area of the flux footprint (Fig. 3) contributed to the success of NEEmod,tower.
OBSERVATIONAL CHALLENGES FOR TUNDRA CO2 FLUX UPSCALING
The nonlinear relationship between NDVI and L creates a nonlinear relationship between NDVI and NEE in PLIRTLE (van Wijk and Williams, 2005; Williams et al., 2008; Stoy et al., 2009b, 2009c). Most growing season NDVI measurements at AT are on the order of 0.7 (Fig. 2, part A), which is near the point where small errors in NDVI create large biases in PLIRTLE-based NEE estimates (see Fig. 7, part A, in Stoy et al., 2009c). This sensitivity, coupled with the sensitivity of NDVI to soil reflectance (Rocha and Shaver, 2009), suggests that including multiple remote-sensing products including the enhanced vegetation index (EVI) may improve multi-scale carbon flux estimates in tundra (Boelman et al., 2003; Rocha and Shaver, 2009). Future work should add multiple observational constraints to improve multi-scale understanding of tundra CO2 flux (Quaife et al., 2008; Hill et al., 2011).
ADDITIONAL APPLICATIONS OF THE PLIRTLE MODEL
PLIRTLE and similar models have been used extensively in recent studies of CO2 exchange in tundra at the chamber (Shaver et al., 2007; Street et al., 2007), tower (Rastetter et al., 2010; Loranty et al., 2011), landscape (Stoy et al., 2009a, 2009b), and regional scales (Loranty et al., 2011). The present analysis demonstrates that it can likewise be used as a strategy for upscaling chamber-scale observations to larger scales in space. It is important to note potential pitfalls of such an approach, including recent findings of lower GPP per unit ground area in tundra patch edges (Fletcher et al., 2012) and the potential role of patterned ground in explaining submeter spatial patterns of ecosystem respiration (Sommerkorn, 2008). We also note that our analysis is restricted to the growing season, and that shoulder and cold-season fluxes from tundra ecosystems may be non-trivial (Larsen et al., 2007a, 2007b; Street et al., 2012); in fact, both NEEmod,tower and NEEmod, footprint performed more poorly later in the growing season when vegetation began to senesce (Fig. 6). Despite potential limitations, we suggest that our upscaling approach may represent an improvement over previous studies that were unable to find convergent estimates of NEE by combining inference from chamber studies and from eddy covariance towers (Fox et al., 2008).
CONCLUSIONS AND FUTURE STUDIES
Our findings emphasize that canopy characteristics (via L) and microclimate (via PPFD and Ta), rather than community composition, can be used for upscaling CO2 flux in tundra. Measurements of plant communities, including the extent of vegetation patches and the distribution of species, are not necessary for explaining or upscaling carbon flux during the growing season (Fox et al., 2008). We note, however, that plant response to regional climate change occurs at the species level. These alterations to plant community composition and the emergent ecosystem consequences, including subsequent changes to L (Jia et al., 2003), must be further studied (Chapin et al., 2005; McGuire et al., 2007; Spadavecchia et al., 2008). Research on present tundra carbon cycling should quantify the magnitude, distribution, and uncertainty of L and microclimate in tundra across scales in space (Loranty et al., 2011). Future research should focus on the response of species distribution and ecosystem attributes to a changing climate (Stow et al., 2004), including the response of climate to changes to vegetation distribution and land cover characteristics (McGuire et al., 2006, 2007; Forbes et al., 2010).
Acknowledgments
We acknowledge funding from the U.S. National Science Foundation (grant numbers OPP-0096523, OPP-0352897, DEB-0087046, and DEB-00895825), University of Edinburgh, and Natural Environment Research Council (NERC) (award NE/D005922/ 1 for the ABACUS consortium; award 03/17 for the NERC ARSF flight that carried the ATM sensor; and support from the NERC Field Spectroscopy Facility [FSF] for the ASD instrument). Stoy acknowledges funding from the Marie Curie European Incoming International Fellowship (project number 237348, TSURF), the National Science Foundation (“Scaling Ecosystem Function: Novel Approaches from MaxEnt and Multiresolution,” DBI #1021095), and Montana State University. We thank Brian Huntley for use of the Trimble dGPS, Dario Papale for providing the flux despiking algorithm, Markus Reichstein and MPI-Jena for research support, Poppy Potter for G-GEOS flight assistance, Terry Callaghan for general support, and Annika Kristofferson at Abisko Scientific Research Station for the provision of meteorological data.