Aircraft‐derived regional scale CO2 fluxes from vegetated drained thaw‐lake basins and interstitial tundra on the Arctic Coastal Plain of Alaska

The landscape surface of the Barrow Peninsula of Alaska is a mosaic of small ponds, thaw lakes, different aged vegetated drained thaw‐lake basins (VDTLBs), and interstitial tundra which have been dynamically formed by both short‐ and long‐term processes. We used a combination of tower‐ and aircraft‐based eddy covariance measurements to characterize the spatial and temporal patterns of CO₂, latent, and sensible heat fluxes along with MODIS NDVI, and were able to scale the aircraft‐based CO₂ fluxes to the 1802 km² Barrow Peninsula region. During typical 2006 summer conditions, the midday hourly CO₂ flux over the region was ?2.04 × 10⁵ kg CO₂ h^?1. The CO₂ fluxes among the interstitial tundra, Ancient, and Old VDTLBs, as well as between the Medium and Young VDTLBs were not significantly different. Combined, the interstitial tundra and Old and Ancient VDTLBs represent～67% of the Barrow Peninsula surface area, accounting for ～59% of the regional flux signal. Although the Medium and Young VDTLBs represent ～11% of the surface area, they account for a large portion, ～35%, of the total regional flux. The remaining ～22% of the surface area are lakes and contributed the remaining ～6% of the total regional flux. Previous studies treated vegetated areas of the region as a single surface type with measurements from a few study sites; doing so could underestimate the regional flux by ～22%. Here, we demonstrate that aircraft‐based systems have the ability to cover large spatial scales while measuring the turbulent fluxes across a number of surfaces and combined with ground‐ and satellite‐based measurements provide a valuable tool for both scaling and validation of regional‐scale fluxes.     

Terrestrial carbon balance in a drier world: the effects of water availability in southwestern North America

Global modeling efforts indicate semiarid regions dominate the increasing trend and interannual variation of net CO₂ exchange with the atmosphere, mainly driven by water availability. Many semiarid regions are expected to undergo climatic drying, but the impacts on net CO₂ exchange are poorly understood due to limited semiarid flux observations. Here we evaluated 121 site‐years of annual eddy covariance measurements of net and gross CO₂ exchange (photosynthesis and respiration), precipitation, and evapotranspiration (ET) in 21 semiarid North American ecosystems with an observed range of 100 – 1000 mm in annual precipitation and records of 4–9 years each. In addition to evaluating spatial relationships among CO₂ and water fluxes across sites, we separately quantified site‐level temporal relationships, representing sensitivity to interannual variation. Across the climatic and ecological gradient, photosynthesis showed a saturating spatial relationship to precipitation, whereas the photosynthesis–ET relationship was linear, suggesting ET was a better proxy for water available to drive CO₂ exchanges after hydrologic losses. Both photosynthesis and respiration showed similar site‐level sensitivity to interannual changes in ET among the 21 ecosystems. Furthermore, these temporal relationships were not different from the spatial relationships of long‐term mean CO₂ exchanges with climatic ET. Consequently, a hypothetical 100‐mm change in ET, whether short term or long term, was predicted to alter net ecosystem production (NEP) by 64 gCm^?2 yr^?1. Most of the unexplained NEP variability was related to persistent, site‐specific function, suggesting prioritization of research on slow‐changing controls. Common temporal and spatial sensitivity to water availability increases our confidence that site‐level responses to interannual weather can be extrapolated for prediction of CO₂ exchanges over decadal and longer timescales relevant to societal response to climate change.     

FLUXNET and modelling the global carbon cycle

Measurements of the net CO₂ flux between terrestrial ecosystems and the atmosphere using the eddy covariance technique have the potential to underpin our interpretation of regional CO₂ source–sink patterns, CO₂ flux responses to forcings, and predictions of the future terrestrial C balance. Information contained in FLUXNET eddy covariance data has multiple uses for the development and application of global carbon models, including evaluation/validation, calibration, process parameterization, and data assimilation. This paper reviews examples of these uses, compares global estimates of the dynamics of the global carbon cycle, and suggests ways of improving the utility of such data for global carbon modelling. Net ecosystem exchange of CO₂ (NEE) predicted by different terrestrial biosphere models compares favourably with FLUXNET observations at diurnal and seasonal timescales. However, complete model validation, particularly over the full annual cycle, requires information on the balance between assimilation and decomposition processes, information not readily available for most FLUXNET sites. Site history, when known, can greatly help constrain the model‐data comparison. Flux measurements made over four vegetation types were used to calibrate the land‐surface scheme of the Goddard Institute for Space Studies global climate model, significantly improving simulated climate and demonstrating the utility of diurnal FLUXNET data for climate modelling. Land‐surface temperatures in many regions cool due to higher canopy conductances and latent heat fluxes, and the spatial distribution of CO₂ uptake provides a significant additional constraint on the realism of simulated surface fluxes. FLUXNET data are used to calibrate a global production efficiency model (PEM). This model is forced by satellite‐measured absorbed radiation and suggests that global net primary production (NPP) increased 6.2% over 1982–1999. Good agreement is found between global trends in NPP estimated by the PEM and a dynamic global vegetation model (DGVM), and between the DGVM and estimates of global NEE derived from a global inversion of atmospheric CO₂ measurements. Combining the PEM, DGVM, and inversion results suggests that CO₂ fertilization is playing a major role in current increases in NPP, with lesser impacts from increasing N deposition and growing season length. Both the PEM and the inversion identify the Amazon basin as a key region for the current net terrestrial CO₂ uptake (i.e. 33% of global NEE), as well as its interannual variability. The inversion's global NEE estimate of −1.2 Pg [C] yr⁻¹ for 1982–1995 is compatible with the PEM‐ and DGVM‐predicted trends in NPP. There is, thus, a convergence in understanding derived from process‐based models, remote‐sensing‐based observations, and inversion of atmospheric data. Future advances in field measurement techniques, including eddy covariance (particularly concerning the problem of night‐time fluxes in dense canopies and of advection or flow distortion over complex terrain), will result in improved constraints on land‐atmosphere CO₂ fluxes and the rigorous attribution of mechanisms to the current terrestrial net CO₂ uptake and its spatial and temporal heterogeneity. Global ecosystem models play a fundamental role in linking information derived from FLUXNET measurements to atmospheric CO₂ variability. A number of recommendations concerning FLUXNET data are made, including a request for more comprehensive site data (particularly historical information), more measurements in undisturbed ecosystems, and the systematic provision of error estimates. The greatest value of current FLUXNET data for global carbon cycle modelling is in evaluating process representations, rather than in providing an unbiased estimate of net CO₂ exchange.     

Soil Carbon Dioxide Flux in Antarctic Dry Valley Ecosystems

The Antarctic dry valleys of southern Victoria Land are extreme desert environments where abiotic factors, such as temperature gradients, parent material, and soil water dynamics, may have a significant influence on soil carbon dioxide (CO₂) flux. Previous measurements of soil respiration have demonstrated very low rates of CO₂ efflux, barely above detection limits. We employed a modified infrared gas-analyzer system that enabled detection of smaller changes in CO₂ concentration in the field than previously possible. We measured diel CO₂ fluxes and monitored soil microclimate at three sites in Taylor Valley. Soil CO₂ flux ranged from –0.1 to 0.15 mol m^–2 s^–1. At two of the three sites, we detected a physically driven flux associated with diel variability in soil temperature. At these sites, CO₂ uptake (negative flux) was associated with dropping soil temperatures, whereas CO₂ evolution (positive flux) was associated with increases in soil temperature. These observations are corroborated by laboratory experiments that suggest that CO₂ flux is influenced by physically driven processes. We discuss four potential mechanisms that may contribute to physically driven gas exchange. Our results suggest there are strong interactions between biological and abiotic controls over soil CO₂ flux in terrestrial ecosystems of the Antarctic dry valleys, and that the magnitude of either may dominate depending on the soil environment and biological activity.     

Towards a coupled paradigm of NH3‐CO2 biosphere–atmosphere exchange modelling

Stomatal conductance, one of the major plant physiological controls within NH₃ biosphere–atmosphere exchange models, is commonly estimated from semi‐empirical multiplicative schemes or simple light‐ and temperature‐response functions. However, due to their inherent parameterization on meteorological proxy variables, instead of a direct measure of stomatal opening, they are unfit for the use in climate change scenarios and of limited value for interpreting field‐scale measurements. Alternatives based on H₂O flux measurements suffer from uncertainties in the partitioning of evapotranspiration at humid sites, as well as a potential decoupling of transpiration from stomatal opening in the presence of hygroscopic particles on leaf surfaces. We argue that these problems may be avoided by directly deriving stomatal conductance from CO₂ fluxes instead. We reanalysed a data set of NH₃ flux measurements based on CO₂‐derived stomatal conductance, confirming the hypothesis that the increasing relevance of stomatal exchange with the onset of vegetation activity caused a rapid decrease of observed NH₃ deposition velocities. Finally, we argue that developing more mechanistic representations of NH₃ biosphere–atmosphere exchange can be of great benefit in many applications. These range from model‐based flux partitioning, over deposition monitoring using low‐cost samplers and inferential modelling, to a direct response of NH₃ exchange to climate change.     

Assessing the eddy covariance technique for evaluating carbon dioxide exchange rates of ecosystems: past, present and future

The eddy covariance technique ascertains the exchange rate of CO₂ across the interface between the atmosphere and a plant canopy by measuring the covariance between fluctuations in vertical wind velocity and CO₂ mixing ratio. Two decades ago, the method was employed to study CO₂ exchange of agricultural crops under ideal conditions during short field campaigns. During the past decade the eddy covariance method has emerged as an important tool for evaluating fluxes of carbon dioxide between terrestrial ecosystems and the atmosphere over the course of a year, and more. At present, the method is being applied in a nearly continuous mode to study carbon dioxide and water vapor exchange at over a hundred and eighty field sites, worldwide. The objective of this review is to assess the eddy covariance method as it is being applied by the global change community on increasingly longer time scales and over less than ideal surfaces. The eddy covariance method is most accurate when the atmospheric conditions (wind, temperature, humidity, CO₂) are steady, the underlying vegetation is homogeneous and it is situated on flat terrain for an extended distance upwind. When the eddy covariance method is applied over natural and complex landscapes or during atmospheric conditions that vary with time, the quantification of CO₂ exchange between the biosphere and atmosphere must include measurements of atmospheric storage, flux divergence and advection. Averaging CO₂ flux measurements over long periods (days to year) reduces random sampling error to relatively small values. Unfortunately, data gaps are inevitable when constructing long data records. Data gaps are generally filled with values produced from statistical and empirical models to produce daily and annual sums of CO₂ exchange. Filling data gaps with empirical estimates do not introduce significant bias errors because the empirical algorithms are derived from large statistical populations. On the other hand, flux measurement errors can be biased at night when winds are light and intermittent. Nighttime bias errors tend to produce an underestimate in the measurement of ecosystem respiration. Despite the sources of errors associated with long‐term eddy flux measurements, many investigators are producing defensible estimates of annual carbon exchange. When measurements come from nearly ideal sites the error bound on the net annual exchange of CO₂ is less than ±50 g C m^?2 yr^?1. Additional confidence in long‐term measurements is growing because investigators are producing values of net ecosystem productivity that are converging with independent values produced by measuring changes in biomass and soil carbon, as long as the biomass inventory studies are conducted over multiple years.     

A scaling approach for quantifying the net CO2 flux of the Kuparuk River Basin,Alaska

Net CO₂ flux measurements conducted during the summer and winter of 1994–96 were scaled in space and time to provide estimates of net CO₂ exchange during the 1995–96 (9 May 1995–8 May 1996) annual cycle for the Kuparuk River Basin, a 9200 km² watershed located in NE Alaska. Net CO₂ flux was measured using dynamic chambers and eddy covariance in moist‐acidic, nonacidic, wet‐sedge, and shrub tundra, which comprise 95% of the terrestrial landscape of the Kuparuk Basin. CO₂ flux data were used as input to multivariate models that calculated instantaneous and daily rates of gross primary production (GPP) and whole‐ecosystem respiration (R) as a function of meteorology and ecosystem development. Net CO₂ flux was scaled up to the Kuparuk Basin using a geographical information system (GIS) consisting of a vegetation map, digital terrain map, dynamic temperature and radiation fields, and the models of GPP and R. Basin‐wide estimates of net CO₂ exchange for the summer growing season (9 May?5 September 1995) indicate that nonacidic tundra was a net sink of ?31.7 ± 21.3 GgC (1 Gg = 10⁹ g), while shrub tundra lost 32.5 ± 6.3 GgC to the atmosphere (negative values denote net ecosystem CO₂ uptake). Acidic and wet sedge tundra were in balance, and when integrated for the entire Kuparuk River Basin (including aquatic surfaces), whole basin summer net CO₂ exchange was estimated to be in balance (?0.9 ± 50.3 GgC). Autumn to winter (6 September 1995–8 May 1996) estimates of net CO₂ flux indicate that acidic, nonacidic, and shrub tundra landforms were all large sources of CO₂ to the atmosphere (75.5 ± 8.3, 96.4 ± 11.4, and 43.3 ± 4.7 GgC for acidic, nonacidic, and shrub tundra, respectively). CO₂ loss from wet sedge surfaces was not substantially different from zero, but the large losses from the other terrestrial landforms resulted in a whole basin net CO₂ loss of 217.2 ± 24.1 GgC during the 1995–96 cold season. When integrated for the 1995–96 annual cycle, acidic (66.4 + 25.25 GgC), nonacidic (64.7 ± 29.2 GgC), and shrub tundra (75.8 ± 8.4 GgC) were substantial net sources of CO₂ to the atmosphere, while wet sedge tundra was in balance (0.4 + 0.8 GgC). The Kuparuk River Basin as a whole was estimated to be a net CO₂ source of 218.1 ± 60.6 GgC over the 1995–96 annual cycle. Compared to direct measurements of regional net CO₂ flux obtained from aircraft‐based eddy covariance, the scaling procedure provided realistic estimates of CO₂ exchange during the summer growing season. Although winter estimates could not be assessed directly using aircraft measurements of net CO₂ exchange, the estimates reported here are comparable to measured values reported in the literature. Thus, we have high confidence in the summer estimates of net CO₂ exchange and reasonable confidence in the winter net CO₂ flux estimates for terrestrial landforms of the Kuparuk river basin. Although there is larger uncertainty in the aquatic estimates, the small surface area of aquatic surfaces in the Kuparuk river basin (≈ 5%) presumably reduces the potential for this uncertainty to result in large errors in basin‐wide CO₂ flux estimates.     

Terrestrial biosphere model performance for inter‐annual variability of land‐atmosphere CO2 exchange

Interannual variability in biosphere‐atmosphere exchange of CO₂ is driven by a diverse range of biotic and abiotic factors. Replicating this variability thus represents the ‘acid test’ for terrestrial biosphere models. Although such models are commonly used to project responses to both normal and anomalous variability in climate, they are rarely tested explicitly against inter‐annual variability in observations. Herein, using standardized data from the North American Carbon Program, we assess the performance of 16 terrestrial biosphere models and 3 remote sensing products against long‐term measurements of biosphere‐atmosphere CO₂ exchange made with eddy‐covariance flux towers at 11 forested sites in North America. Instead of focusing on model‐data agreement we take a systematic, variability‐oriented approach and show that although the models tend to reproduce the mean magnitude of the observed annual flux variability, they fail to reproduce the timing. Large biases in modeled annual means are evident for all models. Observed interannual variability is found to commonly be on the order of magnitude of the mean fluxes. None of the models consistently reproduce observed interannual variability within measurement uncertainty. Underrepresentation of variability in spring phenology, soil thaw and snowpack melting, and difficulties in reproducing the lagged response to extreme climatic events are identified as systematic errors, common to all models included in this study.     

Remote sensing of forest gas exchange: Considerations derived from a tomographic perspective

The global exchange of gas (CO₂, H₂O) and energy (sensible and latent heat) between forest ecosystems and the atmosphere is often assessed using remote sensing (RS) products. Although these products are essential in quantifying the spatial variability of forest–atmosphere exchanges, large uncertainties remain from a measurement bias towards top of canopy fluxes since optical RS data are not sensitive for the vertically integrated forest canopy. We hypothesize that a tomographic perspective opens new pathways to advance upscaling gas exchange processes from leaf to forest stands and larger scales. We suggest a 3D modelling environment comprising principles of ecohydrology and radiative transfer modelling with measurements of micrometeorological variables, leaf optical properties and forest structure, and assess 3D fields of net CO₂ assimilation (A_n) and transpiration (T) in a Swiss temperate forest canopy. 3D simulations were used to quantify uncertainties in gas exchange estimates inherent to RS approaches and model assumptions (i.e. a big‐leaf approximation in modelling approaches). Our results reveal substantial 3D heterogeneity of forest gas exchange with top of canopy A_n and T being reduced by up to 98% at the bottom of the canopy. We show that a simplified use of RS causes uncertainties in estimated vertical gas exchange of up to 300% and that the spatial variation of gas exchange in the footprint of flux towers can exceed diurnal dynamics. We also demonstrate that big‐leaf assumptions can cause uncertainties up to a factor of 10 for estimates of A_n and T. Concluding, we acknowledge the large potential of 3D assessments of gas exchange to unravelling the role of vertical variability and canopy structure in regulating forest–atmosphere gas and energy exchange. Such information allows to systematically link canopy with global scale controls on forest functioning and eventually enables advanced understanding of forest responses to environmental change.     

Interannual variability in carbon dioxide fluxes and flux–climate relationships on grazed and ungrazed northern mixed-grass prairie

The annual carbon (C) budget of grasslands is highly dynamic, dependent on grazing history and on effects of interannual variability (IAV) in climate on carbon dioxide (CO₂) fluxes. Variability in climatic drivers may directly affect fluxes, but also may indirectly affect fluxes by altering the response of the biota to the environment, an effect termed ‘functional change’. We measured net ecosystem exchange of CO₂ (NEE) and its diurnal components, daytime ecosystem CO₂ exchange (P_D) and night‐time respiration (R_E), on grazed and ungrazed mixed‐grass prairie in North Dakota, USA, for five growing seasons. Our primary objective was to determine how climatic anomalies influence variability in CO₂ exchange. We used regression analysis to distinguish direct effects of IAV in climate on fluxes from functional change. Functional change was quantified as the improvement in regression on fitting a model in which slopes of flux–climate relationships vary among years rather than remain invariant. Functional change and direct effects of climatic variation together explained about 20% of variance in weekly means of NEE, P_D, and R_E. Functional change accounted for more than twice the variance in fluxes of direct effects of climatic variability. Grazing did not consistently influence the contribution of functional change to flux variability, but altered which environmental variable best explained year‐to‐year differences in flux–climate slopes, reduced IAV in seasonal means of fluxes, lessened the strength of flux–climate correlations, and increased NEE by reducing R_E relatively more than P_D. Most of these trends are consistent with the interpretation that grazing reduced the influence of plants on ecosystem fluxes. Because relationships between weekly values of fluxes and climatic regulators changed annually, year‐to‐year differences in the C balance of these ecosystems cannot be predicted from knowledge of IAV in climate alone.     

Carbon dioxide fluxes in a spatially and temporally heterogeneous temperate grassland

Landscape position, grazing, and seasonal variation in precipitation and temperature create spatial and temporal variability in soil processes, and plant biomass and composition in grasslands. However, it is unclear how this variation in plant and soil properties affects carbon dioxide (CO₂) fluxes. The aim of this study is to explore the effect of grazing, topographic position, and seasonal variation in soil moisture and temperature on plant assimilation, shoot and soil respiration, and net ecosystem CO₂ exchange (NEE). Carbon dioxide fluxes, vegetation, and environmental variables were measured once a month inside and outside long-term ungulate exclosures in hilltop (dry) to slope bottom (mesic) grassland throughout the 2004 growing season in Yellowstone National Park. There was no difference in vegetation properties and CO₂ fluxes between the grazed and the ungrazed sites. The spatial and temporal variability in CO₂ fluxes were related to differences in aboveground biomass and total shoot nitrogen content, which were both related to variability in soil moisture. All sites were CO₂ sinks (NEE>0) for all our measurments taken throughout the growing season; but CO₂ fluxes were four- to fivefold higher at sites supporting the most aboveground biomass located at slope bottoms, compared to the sites with low biomass located at hilltops or slopes. The dry sites assimilated more CO₂ per gram aboveground biomass and stored proportionally more of the gross-assimilated CO₂ in the soil, compared to wet sites. These results indicate large spatio-temporal variability of CO₂ fluxes and suggest factors that control the variability in Yellowstone National Park.     

Comparison of CO<Subscript>2</Subscript> and H<Subscript>2</Subscript>O fluxes over grassland vegetations measured by the eddy-covariance technique and by open system chamber

Measurements of CO₂ and H₂O fluxes were carried out using two different techniques—eddy-covariance (EC) and open system gas exchange chamber (OC)—during two-years’ period (2003–2004) at three different grassland sites. OC measurements were made during fourteen measurement campaigns. We found good agreement between the OC and EC CO₂ flux values (n = 63, r ² = 0.5323, OC F_CO2 = −0.6408+0.9508 EC F_CO2). The OC F_H2O values were consistently lower than those measured by the EC technique, probably caused by the air stream difference inside and outside the chamber. Adjusting flow rate within the chamber to the natural conditions would be necessary in future OC measurements. In comparison with EC, the OC proved to be a good tool for gas exchange measurements in grassland ecosystems.     

Small-scale hydrological variation determines landscape CO2 fluxes in the high Arctic

We explored the influence of small-scale spatial variation in soil moisture on CO₂ fluxes in the high Arctic. Of five sites forming a hydrological gradient, CO₂ was emitted from the three driest sites and only the wettest site was a net sink of CO₂. Soil moisture was a good predictor of net ecosystem exchange (NEE). Higher gross ecosystem photosynthesis (GEP) was linked to higher bryophyte biomass and activity in response to the moisture conditions. Ecosystem respiration (R _e) rates increased with soil moisture until the soil became anaerobic and then R _e decreased. At well-drained sites R _e was driven by GEP, suggesting substrate and moisture limitation of soil respiration. We propose that spatial variability in soil moisture is a primary driver of NEE.     

Multiscale model intercomparisons of CO2 and H2O exchange rates in a maturing southeastern US pine forest

We compared four existing process‐based stand‐level models of varying complexity (physiological principles in predicting growth, photosynthesis and evapotranspiration, biogeochemical cycles, and stand to ecosystem carbon and evapotranspiration simulator) and a new nested model with 4 years of eddy‐covariance‐measured water vapor (LE) and CO₂ (Fc) fluxes at a maturing loblolly pine forest. The nested model resolves the ‘fast’ CO₂ and H₂O exchange processes using canopy turbulence theories and radiative transfer principles whereas slowly evolving processes were resolved using standard carbon allocation methods modified to improve leaf phenology. This model captured most of the intraannual variations in leaf area index (LAI), net ecosystem exchange (NEE), and LE for this stand in which maximum LAI was not at a steady state. The model comparisons suggest strong linkages between carbon production and LAI variability, especially at seasonal time scales. This linkage necessitates the use of multilayer models to reproduce the seasonal dynamics of LAI, NEE, and LE. However, our findings suggest that increasing model complexity, often justified for resolving faster processes, does not necessarily translate into improved predictive skills at all time scales. Additionally, none of the models tested here adequately captured drought effects on water and CO₂ fluxes. Furthermore, the good performance of some models in capturing flux variability on interannual time scales appears to stem from erroneous LAI dynamics and from sensitivity to droughts that injects unrealistic flux variability at longer time scales.     

Soil moisture dynamics of calcareous grassland under elevated CO2

Water relations of nutrient-poor calcareous grassland under long-term CO₂ enrichment were investigated. Understanding CO₂ effects on soil moisture is critical because productivity in these grasslands is water limited. In general, leaf conductance was reduced at elevated CO₂, but responses strongly depended on date and species. Evapotranspiration (measured as H₂O gas exchange) revealed only small, non-significant reductions at elevated CO₂, indicating that leaf conductance effects were strongly buffered by leaf boundary layer and canopy conductance (leaf area index was not or only marginally increased under elevated CO₂). However, these minute and non-significant responses of water vapour loss accumulated over time and resulted in significantly higher soil moisture in CO₂-enriched plots (gravimetric spot measurements and continuous readings using a network of time-domain reflectometry probes). Differences strongly depended on date, with the smallest effects when soil moisture was very high (after heavy precipitation) and effects were largest at intermediate soil moisture. Elevated CO₂ also affected diurnal soil moisture courses and rewetting of soils after precipitation. We conclude that ecosystem-level controls of the water balance (including soil feedbacks) overshadow by far the physiological effects observed at the leaf level. Indirect effects of CO₂ enrichment mediated by trends in soil moisture will have far-ranging consequences on plant species composition, soil bacterial and faunal activity as well as on soil physical structure and may indirectly also affect hydrology and trace gas emissions and atmospheric chemistry. Received: 21 December 1997 / Accepted: 3 August 1998     

Estimating the uncertainty in annual net ecosystem carbon exchange: spatial variation in turbulent fluxes and sampling errors in eddy-covariance measurements

Above forest canopies, eddy covariance (EC) measurements of mass (CO₂, H₂O vapor) and energy exchange, assumed to represent ecosystem fluxes, are commonly made at one point in the roughness sublayer (RSL). A spatial variability experiment, in which EC measurements were made from six towers within the RSL in a uniform pine plantation, quantified large and dynamic spatial variation in fluxes. The spatial coefficient of variation (CV) of the scalar fluxes decreased with increasing integration time, stabilizing at a minimum that was independent of further lengthening the averaging period (hereafter a ‘stable minimum’). For all three fluxes, the stable minimum (CV=9–11%) was reached at averaging times (τ_p) of 6–7 h during daytime, but higher stable minima (CV=46–158%) were reached at longer τ_p (>12 h) during nighttime. To the extent that decreasing CV of EC fluxes reflects reduction in micrometeorological sampling errors, half of the observed variability at τ_p=30 min is attributed to sampling errors. The remaining half (indicated by the stable minimum CV) is attributed to underlying variability in ecosystem structural properties, as determined by leaf area index, and perhaps associated ecosystem activity attributes. We further assessed the spatial variability estimates in the context of uncertainty in annual net ecosystem exchange (NEE). First, we adjusted annual NEE values obtained at our long‐term observation tower to account for the difference between this tower and the mean of all towers from this experiment; this increased NEE by up to 55 g C m^?2 yr^?1. Second, we combined uncertainty from gap filling and instrument error with uncertainty because of spatial variability, producing an estimate of variability in annual NEE ranging from 79 to 127 g C m^?2 yr^?1. This analysis demonstrated that even in such a uniform pine plantation, in some years spatial variability can contribute ～50% of the uncertainty in annual NEE estimates.     

Spatial variation in regional CO2 exchange for the Kuparuk River Basin, Alaska over the summer growing season

The spatial and temporal patterns in CO₂ flux for the Kuparuk River Basin, a 9200‐km² watershed located in NE Alaska were estimated using the Regional Arctic CO₂ Exchange Simulator (RACES) for the 1994–1995 growing seasons. RACES uses non‐linear models and a Geographical Information System database (GIS) consisting of the normalized difference vegetation index (NDVI) and dynamic temperature and radiation maps. The spatial and temporal patterns in the NDVI during both growing seasons suggest that ecosystem development occurred 2–4 weeks earlier and was relatively more rapid in the southern portion of the Kuparuk River Basin. Rates of gross primary production (GPP) and whole‐ecosystem respiration (R) were 2–4 fold higher in the southern basin than along the arctic coastal plain depending on time of year. The higher rate of GPP estimated for the southern basin was primarily due to higher NDVI values, while the higher R estimated for the southern basin was due in part to higher temperature and the NDVI. While GPP and R showed strong latitudinal trends, spatial and temporal trends in net ecosystem CO₂ exchange (NEE) were much more variable. Thus, while spatial trends in carbon gain (GPP) and loss (R) were highly correlated, small spatial and temporal differences in these large fluxes (GPP and/or R) lead to corresponding large spatial variations in the NEE.     

Coupling between carbon cycling and climate in a high-elevation,subalpine forest: a model-data fusion analysis

Fundamental questions exist about the effects of climate on terrestrial net ecosystem CO₂ exchange (NEE), despite a rapidly growing body of flux observations. One strategy to clarify ecosystem climate–carbon interactions is to partition NEE into its component fluxes, gross ecosystem CO₂ exchange (GEE) and ecosystem respiration (R _E), and evaluate the responses to climate of each component flux. We separated observed NEE into optimized estimates of GEE and R _E using an ecosystem process model combined with 6 years of continuous flux data from the Niwot Ridge AmeriFlux site. In order to gain further insight into the processes underlying NEE, we partitioned R _E into its components: heterotrophic (R _H) and autotrophic (R _A) respiration. We were successful in separating GEE and R _E, but less successful in accurately partitioning R _E into R _A and R _H. Our failure in the latter was due to a lack of adequate contrasts in the assimilated data set to distinguish between R _A and R _H. We performed most model runs at a twice-daily time step. Optimizing on daily-aggregated data severely degraded the model’s ability to separate GEE and R _E. However, we gained little benefit from using a half-hourly time step. The model-data fusion showed that most of the interannual variability in NEE was due to variability in GEE, and not R _E. In contrast to several previous studies in other ecosystems, we found that longer growing seasons at Niwot Ridge were correlated with less net CO₂ uptake, due to a decrease of available snow-melt water during the late springtime photosynthetic period. Warmer springtime temperatures resulted in increased net CO₂ uptake only if adequate moisture was available; when warmer springtime conditions led into mid-summer drought, the annual net uptake declined.     

Spatial variation in landscape‐level CO2 and CH4 fluxes from arctic coastal tundra: influence from vegetation,wetness, and the thaw lake cycle

Regional quantification of arctic CO₂ and CH₄ fluxes remains difficult due to high landscape heterogeneity coupled with a sparse measurement network. Most of the arctic coastal tundra near Barrow, Alaska is part of the thaw lake cycle, which includes current thaw lakes and a 5500‐year chronosequence of vegetated thaw lake basins. However, spatial variability in carbon fluxes from these features remains grossly understudied. Here, we present an analysis of whole‐ecosystem CO₂ and CH₄ fluxes from 20 thaw lake cycle features during the 2011 growing season. We found that the thaw lake cycle was largely responsible for spatial variation in CO₂ flux, mostly due to its control on gross primary productivity (GPP). Current lakes were significant CO₂ sources that varied little. Vegetated basins showed declining GPP and CO₂ sink with age (R² = 67% and 57%, respectively). CH₄ fluxes measured from a subset of 12 vegetated basins showed no relationship with age or CO₂ flux components. Instead, higher CH₄ fluxes were related to greater landscape wetness (R² = 57%) and thaw depth (additional R² = 28%). Spatial variation in CO₂ and CH₄ fluxes had good satellite remote sensing indicators, and we estimated the region to be a small CO₂ sink of ?4.9 ± 2.4 (SE) g C m^?2 between 11 June and 25 August, which was countered by a CH₄ source of 2.1 ± 0.2 (SE) g C m^?2. Results from our scaling exercise showed that developing or validating regional estimates based on single tower sites can result in significant bias, on average by a factor 4 for CO₂ flux and 30% for CH₄ flux. Although our results are specific to the Arctic Coastal Plain of Alaska, the degree of landscape‐scale variability, large‐scale controls on carbon exchange, and implications for regional estimation seen here likely have wide relevance to other arctic landscapes.     

Linking Annual N2O Emission in Organic Soils to Mineral Nitrogen Input as Estimated by Heterotrophic Respiration and Soil C/N Ratio

Organic soils are an important source of N₂O, but global estimates of these fluxes remain uncertain because measurements are sparse. We tested the hypothesis that N₂O fluxes can be predicted from estimates of mineral nitrogen input, calculated from readily-available measurements of CO₂ flux and soil C/N ratio. From studies of organic soils throughout the world, we compiled a data set of annual CO₂ and N₂O fluxes which were measured concurrently. The input of soil mineral nitrogen in these studies was estimated from applied fertilizer nitrogen and organic nitrogen mineralization. The latter was calculated by dividing the rate of soil heterotrophic respiration by soil C/N ratio. This index of mineral nitrogen input explained up to 69% of the overall variability of N₂O fluxes, whereas CO₂ flux or soil C/N ratio alone explained only 49% and 36% of the variability, respectively. Including water table level in the model, along with mineral nitrogen input, further improved the model with the explanatory proportion of variability in N₂O flux increasing to 75%. Unlike grassland or cropland soils, forest soils were evidently nitrogen-limited, so water table level had no significant effect on N₂O flux. Our proposed approach, which uses the product of soil-derived CO₂ flux and the inverse of soil C/N ratio as a proxy for nitrogen mineralization, shows promise for estimating regional or global N₂O fluxes from organic soils, although some further enhancements may be warranted.