Publications by Year: 2009

2009
MO Roman, CB Schaaf, CE Woodcock, AH Strahler, XY Yang, RH Braswell, PS Curtis, KJ Davis, D Dragoni, ML Goulden, LH Gu, DY Hollinger, TE Kolb, TP Meyers, JW Munger, JL Privette, AD Richardson, TB Wilson, and SC Wofsy. 11/16/2009. “The MODIS (Collection V005) BRDF/albedo product: Assessment of spatial representativeness over forested landscapes.” Remote Sensing of Environment, 113, 11, Pp. 2476-2498. Publisher's VersionAbstract
A new methodology for establishing the spatial representativeness of tower albedo measurements that are routinely used in validation of satellite retrievals from global land surface albedo and reflectance anisotropy products is presented. This method brings together knowledge of the intrinsic biophysical properties of a measurement site, and the surrounding landscape to produce a number of geostatistical attributes that describe the overall variability, spatial extent, strength of the spatial correlation, and spatial structure of surface albedo patterns at separate seasonal periods throughout the year. Variogram functions extracted from Enhanced Thematic Mapper Plus (ETM+) retrievals of surface albedo using multiple spatial and temporal thresholds were used to assess the degree to which a given point (tower) measurement is able to capture the intrinsic variability of the immediate landscape extending to a satellite pixel. A validation scheme was implemented over a wide range of forested landscapes, looking at both deciduous and coniferous sites, from tropical to boreal ecosystems. The experiment focused on comparisons between tower measurements of surface albedo acquired at local solar noon and matching retrievals from the MODerate Resolution Imaging Spectroradiometer (MODIS) (Collection V005) Bidirectional Reflectance Distribution Function (BRDF)/albedo algorithm. Assessments over a select group of field stations with comparable landscape features and daily retrieval scenarios further demonstrate the ability of this technique to identify measurement sites that contain the intrinsic spatial and seasonal features of surface albedo over sufficiently large enough footprints for use in modeling and remote sensing studies. This approach, therefore, improves our understanding of product uncertainty both in terms of the representativeness of the field data and its relationship to the larger satellite pixel.
WP Yuan, YQ Luo, AD Richardson, R Oren, S Luyssaert, IA Janssens, R Ceulemans, XH Zhou, T Grunwald, M Aubinet, C Berhofer, DD Baldocchi, JQ Chen, AL Dunn, JL Deforest, D Dragoni, AH Goldstein, E Moors, JW Munger, RK Monson, AE Suyker, G Star, RL Scott, J Tenhunen, SB Verma, T Vesala, and SC Wofsy. 11/2009. “Latitudinal patterns of magnitude and interannual variability in net ecosystem exchange regulated by biological and environmental variables.” Global Change Biology, 15, 12, Pp. 2905-2920. Publisher's VersionAbstract
Over the last two and half decades, strong evidence showed that the terrestrial ecosystems are acting as a net sink for atmospheric carbon. However the spatial and temporal patterns of variation in the sink are not well known. In this study, we examined latitudinal patterns of interannual variability (IAV) in net ecosystem exchange (NEE) of CO2 based on 163 site-years of eddy covariance data, from 39 northern-hemisphere research sites located at latitudes ranging from ∼29°N to ∼64°N. We computed the standard deviation of annual NEE integrals at individual sites to represent absolute interannual variability (AIAV), and the corresponding coefficient of variation as a measure of relative interannual variability (RIAV). Our results showed decreased trends of annual NEE with increasing latitude for both deciduous broadleaf forests and evergreen needleleaf forests. Gross primary production (GPP) explained a significant proportion of the spatial variation of NEE across evergreen needleleaf forests, whereas, across deciduous broadleaf forests, it is ecosystem respiration (Re). In addition, AIAV in GPP and Re increased significantly with latitude in deciduous broadleaf forests, but AIAV in GPP decreased significantly with latitude in evergreen needleleaf forests. Furthermore, RIAV in NEE, GPP, and Re appeared to increase significantly with latitude in deciduous broadleaf forests, but not in evergreen needleleaf forests. Correlation analyses showed air temperature was the primary environmental factor that determined RIAV of NEE in deciduous broadleaf forest across the North American sites, and none of the chosen climatic factors could explain RIAV of NEE in evergreen needleleaf forests. Mean annual NEE significantly increased with latitude in grasslands. Precipitation was dominant environmental factor for the spatial variation of magnitude and IAV in GPP and Re in grasslands.
RF Grant, LR Hutyra, RC de Oliveira, JW Munger, SR Saleska, and SC. Wofsy. 8/1/2009. “Modeling the carbon balance of Amazonian rain forests: resolving ecological controls on net ecosystem productivity.” Publisher, 79, 3, Pp. 445-463. Publisher's VersionAbstract
There is still much uncertainty about ecological controls on the rate and direction of net CO2 exchange by tropical rain forests, in spite of their importance to global C cycling. These controls are thought to arise from hydrologic and nutrient constraints to CO2 fixation caused by seasonality of precipitation and adverse chemical properties of some major tropical soil types. Using the ecosystem model ecosys, we show that water uptake to a depth of 8 m avoids constraints to CO2 and energy exchange from soil drying during five-month dry seasons typical for eastern Amazonian forests. This avoidance in the model was tested with eddy covariance (EC) measurements of CO2 and energy fluxes during 2003 and 2004 over an old-growth forest on an acidic, nutrient-poor oxisol in the Tapajós National Forest (TNF) in Pará, Brazil. Modeled CO2 fixation was strongly constrained by slow phosphorus (P) uptake caused by low soil pH. Daytime CO2 influxes in the model were in close agreement with EC measurements (R2 > 0.8) during both wet and dry seasons. Both modeled and measured fluxes indicated that seasonality of precipitation affected CO2 and energy exchange more through its effect on radiation and air temperature than on soil water content. When aggregated to a yearly scale, modeled and gap-filled EC CO2 fluxes indicated that old-growth forest stands in the TNF remained within 100 g C·m−2·yr−1 of C neutrality in the absence of major disturbance. Annual C transformations in ecosys were further corroborated by extensive biometric measurements taken in the TNF and elsewhere in the Amazon basin, which also indicated that old-growth forests were either small C sources or small C sinks. Long-term model runs suggested that rain forests could be substantial C sinks for several decades while regenerating after stand-replacing disturbances, but would gradually decline toward C neutrality thereafter. The time course of net ecosystem productivity (NEP) in the model depended upon annual rates of herbivory and tree mortality, which were based on site observations as affected by weather (e.g., El Niño Southern Oscillation [ENSO] events). This dependence suggests that rain forest NEP is strongly controlled by disturbance as well as by weather.
S Luyssaert, M Reichstein, ED Schulze, IA Janssens, BE Law, D Papale, D Dragoni, ML Goulden, A Granier, WL Kutsch, S Linder, G Matteucci, E Moors, JW Munger, K Pilegaard, M Saunders, and EM Falge. 7/18/2009. “Toward a consistency cross-check of eddy covariance flux-based and biometric estimates of ecosystem carbon balance.” Global Biogeochemical Cycles, 23, GB3009. Publisher's VersionAbstract
Quantification of an ecosystem's carbon balance and its components is pivotal for understanding both ecosystem functioning and global cycling. Several methods are being applied in parallel to estimate the different components of the CO2 balance. However, different methods are subject to different sources of error. Therefore, it is necessary that site level component estimates are cross-checked against each other before being reported. Here we present a two-step approach for testing the accuracy and consistency of eddy covariance–based gross primary production (GPP) and ecosystem respiration (Re) estimates with biometric measurements of net primary production (NPP), autotrophic (Ra) and heterotrophic (Rh) respiration. The test starts with closing the CO2 balance to account for reasonable errors in each of the component fluxes. Failure to do so within the constraints will classify the flux estimates on the site level as inconsistent. If the CO2 balance can be closed, the test continues by comparing the closed site level Ra/GPP with the Rh/GPP ratio. The consistency of these ratios is then judged against expert knowledge. Flux estimates of sites that pass both steps are considered consistent. An inconsistent ratio is not necessarily incorrect but provides a signal for careful data screening that may require further analysis to identify the possible biological reasons of the unexpected ratios. We reviewed the literature and found 16 sites, out of a total of 529 research forest sites, that met the data requirements for the consistency test. Thirteen of these sites passed both steps of the consistency cross-check. Subsequently, flux ratios (NPP/GPP, Rh/NPP, Rh/Re, and Re/GPP) were calculated for the consistent sites. Similar ratios were observed at sites which lacked information to check consistency, indicating that the flux data that are currently used for validating models and testing ecological hypotheses are largely consistent across a wide range of site productivities. Confidence in the output of flux networks could be further enhanced if the required fluxes are independently estimated at all sites for multiple years and harmonized methods are used.
Y. Wang, J Hao, M.B. McElroy, JW Munger, H Ma, D Chen, and C Nielsen. 7/11/2009. “Ozone air quality during the 2008 Beijing Olympics: effectiveness of emission restrictions.” Atmospheric Chemistry and Physics, 9, 14, Pp. 5237-5251. Publisher's VersionAbstract
A series of aggressive measures was launched by the Chinese government to reduce pollutant emissions from Beijing and surrounding areas during the Olympic Games. Observations at Miyun, a rural site 100 km downwind of the Beijing urban center, show significant decreases in concentrations of O3, CO, NOy, and SO2 during August 2008, relative to August 2006–2007. The mean daytime mixing ratio of O3 was lower by about 15 ppbv, reduced to 50 ppbv, in August 2008. The relative reductions in daytime SO2, CO, and NOy were 61%, 25%, and 21%, respectively. Changes in SO2 and in species correlations from 2007 to 2008 indicate that emissions of SO2, CO, and NOx were reduced at least by 60%, 32%, and 36%, respectively, during the Olympics. Analysis of meteorological conditions and interpretation of observations using a chemical transport model suggest that although the day-to-day variability in ozone is driven mostly by meteorology, the reduction in emissions of ozone precursors associated with the Olympic Games had a significant contribution to the observed decrease in O3 during August 2008, accounting for 80% of the O3 reduction for the month as a whole and 45% during the Olympics Period (8–24 August). The model predicts that emission restrictions such as those implemented during the Olympics can affect O3 far beyond the Beijing urban area, resulting in reductions in boundary layer O3 of 2–10 ppbv over a large region of the North China Plain and Northeastern China.
E van Gorsel, N Delpierre, R Leuning, A Black, JW Munger, S Wofsy, M Aubinet, C Feigenwinter, J Beringer, D Bonal, BZ Chen, JQ Chen, R Clement, KJ Davis, AR Desai, D Dragoni, S Etzold, T Grunwald, LH Gu, B Heinesch, LR Hutyra, WWP Jans, W Kutsch, BE Law, MY Leclerc, I Mammarella, L Montagnani, A Noormets, C Rebmann, and S. Wharton. 6/29/2009. “Estimating nocturnal ecosystem respiration from the vertical turbulent flux and change in storage of CO2.” Agricultural and Forest Meteorology, 149, 11, Pp. 1919-1930. Publisher's VersionAbstract
Micrometeorological measurements of night time ecosystem respiration can be systematically biased when stable atmospheric conditions lead to drainage flows associated with decoupling of air flow above and within plant canopies. The associated horizontal and vertical advective fluxes cannot be measured using instrumentation on the single towers typically used at micrometeorological sites. A common approach to minimize bias is to use a threshold in friction velocity, u*, to exclude periods when advection is assumed to be important, but this is problematic in situations when in canopy flows are decoupled from the flow above. Using data from 25 flux stations in a wide variety of forest ecosystems globally, we examine the generality of a novel approach to estimating nocturnal respiration developed by van Gorsel et al. (2007, Tellus, 59B, 307-403). The approach is based on the assumption that advection is small relative to the vertical turbulent flux (Fc) and change in storage (Fs) of CO2 in the few hours after sundown. The sum of Fc and Fs reach a maximum during this period which is used to derive a temperature response function for ecosystem respiration. Measured hourly soil temperatures are then used with this function to estimate respiration RRmax. The new approach yielded excellent agreement with (1) independent measurements using respiration chambers, (2) with estimates using ecosystem light-response curves of Fc + Fs extrapolated to zero light, RLRC, and (3) with a detailed process-based forest ecosystem model, Rcast. At most sites respiration rates estimated using the u*-filter, Rust, were smaller than RRmax and RLRC. Agreement of our approach with independent measurements indicates that RRmax provides an excellent estimate of night time ecosystem respiration.
AL Dunn, SC Wofsy, and AV Bright. 4/1/2009. “Landscape heterogeneity, soil climate, and carbon exchange in a boreal black spruce forest.” Ecological Applications, 19, 2, Pp. 495-504. Publisher's VersionAbstract
We measured soil climate and the turbulent fluxes of CO2, H2O, heat, and momentum on short towers (2 m) in a 160-yr-old boreal black spruce forest in Manitoba, Canada. Two distinct land cover types were studied: a Sphagnum-dominated wetland, and a feathermoss (Pleurozium and Hylocomium)-dominated upland, both lying within the footprint of a 30-m tower, which has measured whole-forest carbon exchange since 1994. Peak summertime uptake of CO2, was higher in the wetland than for the forest as a whole due to the influence of deciduous shrubs. Soil respiration rates in the wetland were approximately three times larger than in upland soils, and 30% greater than the mean of the whole forest, reflecting decomposition of soil organic matter. Soil respiration rates in the wetland were regulated by soil temperature, which was in turn influenced by water table depth through effects on soil heat capacity and conductivity. Warmer soil temperatures and deeper water tables favored increased heterotrophic respiration. Wetland drainage was limited by frost during the first half of the growing season, leading to high, perched water tables, cool soil temperatures, and much lower respiration rates than observed later in the growing season. Whole-forest evapotranspiration increased as water tables dropped, suggesting that photosynthesis in this forest was rarely subject to water stress. Our data indicate positive feedback between soil temperature, seasonal thawing, heterotrophic respiration, and evapotranspiration. As a result, climate warming could cause covariant changes in soil temperature and water table depths that may stimulate photosynthesis and strongly promote efflux of CO2 from peat soils in boreal wetlands.
AD Richardson, DY Hollinger, DB Dail, JT Lee, JW Munger, and J O'Keefe. 3/1/2009. “Influence of spring phenology on seasonal and annual carbon balance in two contrasting New England forests.” Tree Physiology, 29, 3, Pp. 321-331. Publisher's VersionAbstract
Spring phenology is thought to exert a major influence on the carbon (C) balance of temperate and boreal ecosystems. We investigated this hypothesis using four spring onset phenological indicators in conjunction with surface–atmosphere CO2 exchange data from the conifer-dominated Howland Forest and deciduous-dominated Harvard Forest AmeriFlux sites. All phenological measures, including CO2 source–sink transition dates, could be well predicted on the basis of a simple two-parameter spring warming model, indicating good potential for improving the representation of phenological transitions and their dynamic responsiveness to climate variability in land surface models. The date at which canopy-scale photosynthetic capacity reached a threshold value of 12 μmol m−2 s−1 was better correlated with spring and annual flux integrals than were either deciduous or coniferous bud burst dates. For all phenological indicators, earlier spring onset consistently, but not always significantly, resulted in higher gross primary productivity (GPP) and ecosystem respiration (RE) for both seasonal (spring months, April–June) and annual flux integrals. The increase in RE was less than that in GPP; depending on the phenological indicator used, a one-day advance in spring onset increased springtime net ecosystem productivity (NEP) by 2–4 g C m−2 day−1. In general, we could not detect significant differences between the two forest types in response to earlier spring, although the response to earlier spring was generally more pronounced for Harvard Forest than for Howland Forest, suggesting that future climate warming may favor deciduous species over coniferous species, at least in this region. The effect of earlier spring tended to be about twice as large when annual rather than springtime flux integrals were considered. This result is suggestive of both immediate and lagged effects of earlier spring onset on ecosystem C cycling, perhaps as a result of accelerated N cycling rates and cascading effects on N uptake, foliar N concentrations and photosynthetic capacity.
MS Zahniser, DD Nelson, JB McManus, S Herndon, E Wood, JH Shorter, BW Lee, GH Santoni, R Jimenez, BC Daube, S Park, EA Kort, and SC Wofsy. 1/26/2009. “Infrared QC laser applications to field measurements of atmospheric trace gas sources and sinks in environmental research: enhance capabilities using continuous wave QCLs.” SPIE, 7222. Publisher's VersionAbstract
The advent of continuous wave quantum cascade lasers operating at near room temperature has greatly expanded the capability of spectroscopic detection of atmospheric trace gases using infrared absorption at wavelengths from 4 to 12 μm. The high optical power, narrow line width, and high degree of single mode purity result in minimal fractional absorptions of 5x10-6 Hz-1/2 detectable in direct absorption with path lengths up to 210 meters. The Allan plot minima correspond to a fractional absorbance of 1x10-6 or a minimum absorption per unit path length 5x10-11 cm-1 in 50s. This allows trace gas mixing ratio detection limits in the low part-per-trillion (1 ppt = 10-12) range for many trace gases of atmospheric interest. We present ambient measurements of NO2 with detection precision of 10 ppt Hz-1/2. The detection precision for the methane isotopologue 13CH4 is 25 ppt Hz-1/2 which allows direct measurements of ambient ratios of 13CH4/12CH4 with a precision of 0.5‰ in 100 s without pre-concentration. Projections are given for detection limits for other gases including COS, HONO and HCHO as CWRT lasers become available at appropriate wavelengths.
D Medvigy, SC Wofsy, JW Munger, DY Hollinger, and PR Moorcroft. 1/17/2009. “Mechanistic scaling of ecosystem function and dynamics in space and time: Ecosystem Demography model version 2.” Journal of Geophysical Research: Biogeosciences, 114, G01002. Publisher's VersionAbstract
Insights into how terrestrial ecosystems affect the Earth's response to changes in climate and rising atmospheric CO2 levels rely heavily on the predictions of terrestrial biosphere models (TBMs). These models contain detailed mechanistic representations of biological processes affecting terrestrial ecosystems; however, their ability to simultaneously predict field-based measurements of terrestrial vegetation dynamics and carbon fluxes has remained largely untested. In this study, we address this issue by developing a constrained implementation of a new structured TBM, the Ecosystem Demography model version 2 (ED2), which explicitly tracks the dynamics of fine-scale ecosystem structure and function. Carbon and water flux measurements from an eddy-flux tower are used in conjunction with forest inventory measurements of tree growth and mortality at Harvard Forest (42.5°N, 72.1°W) to estimate a number of important but weakly constrained model parameters. Evaluation against a decade of tower flux and forest dynamics measurements shows that the constrained ED2 model yields greatly improved predictions of annual net ecosystem productivity, carbon partitioning, and growth and mortality dynamics of both hardwood and conifer trees. The generality of the model formulation is then evaluated by comparing the model's predictions against measurements from two other eddy-flux towers and forest inventories of the northeastern United States and Quebec. Despite the markedly different composition throughout this region, the optimized model realistically predicts observed patterns of carbon fluxes and tree growth. These results demonstrate how TBMs parameterized with field-based measurements can provide quantitative insight into the underlying biological processes governing ecosystem composition, structure, and function at larger scales.
HR da Rocha, AO Manzi, OM Cabral, SD Miller, ML Goulden, SR Saleska, NR Coupe, SC Wofsy, LS Borma, P. Artaxo, G Vourlitis, JS Nogueira, FL Cardoso, AD Nobre, B Kruijt, HC Freitas, C von Randow, RG Aguiar, and JF Maia. 1/17/2009. “Patterns of water and heat flux across a biome gradient from tropical forest to savanna in Brazil.” Journal of Geophysical Research: Biogeosciences, 114, G00B12. Publisher's VersionAbstract
We investigated the seasonal patterns of water vapor and sensible heat flux along a tropical biome gradient from forest to savanna. We analyzed data from a network of flux towers in Brazil that were operated within the Large-Scale Biosphere-Atmosphere Experiment in Amazonia (LBA). These tower sites included tropical humid and semideciduous forest, transitional forest, floodplain (with physiognomies of cerrado), and cerrado sensu stricto. The mean annual sensible heat flux at all sites ranged from 20 to 38 Wm−2, and was generally reduced in the wet season and increased in the late dry season, coincident with seasonal variations of net radiation and soil moisture. The sites were easily divisible into two functional groups based on the seasonality of evaporation: tropical forest and savanna. At sites with an annual precipitation above 1900 mm and a dry season length less than 4 months (Manaus, Santarem and Rondonia), evaporation rates increased in the dry season, coincident with increased radiation. Evaporation rates were as high as 4.0 mm d−1 in these evergreen or semidecidous forests. In contrast, ecosystems with precipitation less than 1700 mm and a longer dry season (Mato Grosso, Tocantins and São Paulo) showed clear evidence of reduced evaporation in the dry season. Evaporation rates were as low as 2.5 mm d−1 in the transitional forests and 1 mm d−1 in the cerrado. The controls on evapotranspiration seasonality changed along the biome gradient, with evaporative demand (especially net radiation) playing a more important role in the wetter forests, and soil moisture playing a more important role in the drier savannah sites.