Chuixiang Yi

Long-Term Carbon Dioxide Fluxes from a Very Tall Tower in a Northern Forest: Flux Measurement Methodology

Methodology for determining fluxes of CO 2 and H2O vapor with the eddy-covariance method using data from instruments on a 447-m tower in the forest of northern Wisconsin is addressed. The primary goal of this study is the validation of the methods used to determine the net ecosystem exchange of CO 2. Two-day least squares fits coupled with 30-day running averages limit calibration error of infrared gas analyzers for CO2 and H2O signals to 2%‐3%. Sonic anemometers are aligned with local streamlines by fitting a sine function to tilt and wind direction averages, and fitting a third-order polynomial to the residual. Lag times are determined by selecting the peak in lagged covariance with an error of 1.5%‐2% for CO2 and 1% for H2O vapor. Theory and a spectral fit method allow determination of the underestimation in CO2 flux ( ,5% daytime, ,12% nighttime) and H2O vapor flux ( ,21%), which is due to spectral degradation induced by long air-sampling tubes. Scale analysis finds 0.5-h flux averaging periods are sufficient to measure all flux scales at 30-m height, but 1 h is necessary at higher levels, and random errors in the flux measurements due to limited sampling of atmospheric turbulence are fairly large ( 15%‐20% for CO2 and 20%‐40% for H2O vapor at lower levels for a 1-h period).

Influence of advection on measurements of the net ecosystem‐atmosphere exchange of CO2 from a very tall tower

In most studies of the net ecosystem-atmosphere exchange of CO2 (NEE) using tower-based eddy covariance (EC) systems it has been assumed that advection is negligible. In this study we use a scalar conservation budget method to estimate the contribution of advection to NEE measurements from a very tall tower in northern Wisconsin. We examine data for June-August 1997. Measured NEE0, calculated as the sum of the EC flux plus the rate of change of storage below the EC measurement level, is expected to be constant with measurement height, and we take the differences between levels as a measure of advection. We find that the average difference in total advection DFCadtot between 30 and 122 m is as large as 6 mmol m 22 s 21 during the morning transition from stable to convective conditions and the average difference DFCadtot between 122 and 396 m is as large as 4 mmol m 22 s 21 during daytime. For the month of July, advection between 30 and 122 m is 27% of the diurnally integrated NEE0 at 122 m, and advection between 122 and 396 m accounts for 5% of the NEE0 observed at 396 m. The observed differences of advection often have significant correlation with the vertical integral of wind speed within the same layer. This indicates that the horizontal advection contribution to NEE could be significant. Direct observations of the vertical gradient in CO2 show that DFCadtot cannot be explained by vertical advection alone. It is hypothesized that differing flux footprints and pooling of CO2 in the heterogeneous landscape causes the advection contribution. The magnitudes of the total advection component FCadtot of NEE at the 30 m level are roughly estimated by a linear extrapolation. A peak in F Cadtot at 30 m of ;3 mmol m 22 s 21 during the morning transition is predicted for all three months. The July integrated F Cadtot is estimated to be 10% of the diurnally integrated NEE0 at 30 m.

Long-Term Observations of the Dynamics of the Continental Planetary Boundary Layer

Time series of mixed layer depth, zi, and stable boundary layer height from March through October of 1998 are derived from a 915-MHz boundary layer profiling radar and CO 2 mixing ratio measured from a 447-m tower in northern Wisconsin. Mixed layer depths from the profiler are in good agreement with radiosonde measurements. Maximum zi occurs in May, coincident with the maximum daytime surface sensible heat flux. Incoming radiation is higher in June and July, but a greater proportion is converted to latent heat by photosynthesizing vegetation. An empirical relationship between zi and the square root of the cumulative surface virtual potential temperature flux is obtained ( r 2 5 0.98) allowing estimates of zi from measurements of virtual potential temperature flux under certain conditions. In fair-weather conditions the residual mixed layer top was observed by the profiler on several nights each month. The synoptic mean vertical velocity (subsidence rate) is estimated from the temporal evolution of the residual mixed layer height during the night. The influence of subsidence on the evolution of the mixed, stable, and residual layers is discussed. The CO 2 jump across the inversion at night is also estimated from the tower measurements.

Modeling and measuring the nocturnal drainage flow in a high‐elevation, subalpine forest with complex terrain

[1] The nocturnal drainage flow of air causes significant uncertainty in ecosystem CO2, H2O, and energy budgets determined with the eddy covariance measurement approach. In this study, we examined the magnitude, nature, and dynamics of the nocturnal drainage flow in a subalpine forest ecosystem with complex terrain. We used an experimental approach involving four towers, each with vertical profiling of wind speed to measure the magnitude of drainage flows and dynamics in their occurrence. We developed an analytical drainage flow model, constrained with measurements of canopy structure and SF6 diffusion, to help us interpret the tower profile results. Model predictions were in good agreement with observed profiles of wind speed, leaf area density, and wind drag coefficient. Using theory, we showed that this one-dimensional model is reduced to the widely used exponential wind profile model under conditions where vertical leaf area density and drag coefficient are uniformly distributed. We used the model for stability analysis, which predicted the presence of a very stable layer near the height of maximum leaf area density. This stable layer acts as a flow impediment, minimizing vertical dispersion between the subcanopy air space and the atmosphere above the canopy. The prediction is consistent with the results of SF6 diffusion observations that showed minimal vertical dispersion of nighttime, subcanopy drainage flows. The stable within-canopy air layer coincided with the height of maximum wake-to-shear production ratio. We concluded that nighttime drainage flows are restricted to a relatively shallow layer of air beneath the canopy, with little vertical mixing across a relatively long horizontal fetch. Insight into the horizontal and vertical structure of the drainage flow is crucial for understanding the magnitude and dynamics of the mean advective CO2 flux that becomes significant during stable nighttime conditions and are typically missed during measurement of the turbulent CO2 flux. The model and interpretation provided in this study should lead to research strategies for the measurement of these advective fluxes and their inclusion in the overall mass balance for CO2 at this site with complex terrain.

Regional carbon dioxide fluxes from mixing ratio data

We examine the atmospheric budget of CO2 at temperate continental sites in the Northern Hemisphere. On a monthly time scale both surface exchange and atmospheric transport are important in determining the rate of change of CO2 mixing ratio at these sites. Vertical differences between the atmospheric boundary layer and free troposphere over the continent are generally greater than large-scale zonal gradients such as the difference between the free troposphere over the continent and the marine boundary layer. Therefore, as a first approximation we parametrize atmospheric transport as a vertical exchange term related to the vertical gradient of CO2 and the mean vertical velocity from the NCEP Reanalysis model. Horizontal advection is assumed to be negligible in our simple analysis. We then calculate the net surface exchange of CO2 from CO2 mixing ratio measurements at four tower sites. The results provide estimates of the surface exchange that are representative of a regional scale (i.e. ∼106 km2). Comparison with direct, local-scale (eddy covariance) measurements of net exchange with the ecosystems around the towers are reasonable after accounting for anthropogenic CO2 emissions within the larger area represented by the mixing ratio data. A network of tower sites and frequent aircraft vertical profiles, separated by several hundred kilometres, where CO2 is accurately measured would provide data to estimate horizontal and vertical advection and hence provide a means to derive net CO2 fluxes on a regional scale. At present CO2 mixing ratios are measured with sufficient accuracy relative to global reference gas standards at only a few continental sites. The results also confirm that flux measurements from carefully sited towers capture seasonal variations representative of large regions, and that the midday CO2 mixing ratios sampled in the atmospheric surface layer similarly capture regional and seasonal variability in the continental CO2 budget.

Momentum Transfer within Canopies

Abstract To understand the basic characteristics of the observed S-shaped wind profile and the exponential flux profile within forest canopies, three hypotheses are postulated. The relationship between these fundamental profiles is well established by combining the postulated hypotheses with momentum equations. Robust agreements between theoretical predictions and observations indicate that the nature of momentum transfer within canopies can be well understood by combining the postulated hypotheses and momentum equations. The exponential Reynolds stress profiles were successfully predicted by the leaf area index (LAI) profile alone. The characteristics of the S-shaped wind profile were theoretically explained by the plant morphology and local drag coefficient distribution. Predictions of maximum drag coefficient were located around the maximum leaf area level for most forest canopies but lower than the maximum leaf area level for a corn canopy. A universal relationship of the Reynolds stress between the t...

Tipping point of a conifer forest ecosystem under severe drought

Drought-induced tree mortality has recently received considerable attention. Questions have arisen over the necessary intensity and duration thresholds of droughts that are sufficient to trigger rapid forest declines. The values of such tipping points leading to forest declines due to drought are presently unknown. In this study, we have evaluated the potential relationship between the level of tree growth and concurrent drought conditions with data of the tree growth-related ring width index (RWI) of the two dominant conifer species (Pinus edulis and Pinus ponderosa) in the Southwestern United States (SWUS) and the meteorological drought-related standardized precipitation evapotranspiration index (SPEI). In this effort, we determined the binned averages of RWI and the 11 month SPEI within the month of July within each bin of 30 of RWI in the range of 0–3000. We found a significant correlation between the binned averages of RWI and SPEI at the regional-scale under dryer conditions. The tipping point of forest declines to drought is predicted by the regression model as SPEItp = −1.64 and RWItp = 0, that is, persistence of the water deficit (11 month) with intensity of −1.64 leading to negligible growth for the conifer species. When climate conditions are wetter, the correlation between the binned averages of RWI and SPEI is weaker which we believe is most likely due to soil water and atmospheric moisture levels no longer being the dominant factor limiting tree growth. We also illustrate a potential application of the derived tipping point (SPEItp = −1.64) through an examination of the 2002 extreme drought event in the SWUS conifer forest regions. Distinguished differences in remote-sensing based NDVI anomalies were found between the two regions partitioned by the derived tipping point.

A nonparametric method for separating photosynthesis and respiration components in CO2 flux measurements

[1] Future climate change is expected to affect ecosystematmosphere CO2 exchange, particularly through the influence of temperature. To date, however, few studies have shown that differences in the response of net ecosystem CO2 exchange (NEE) to temperature among ecosystems can be explained by differences in the photosynthetic and respiratory processes that compose NEE. Using a new nonparametric statistical model, we analyzed data from four forest ecosystems. We observed that differences among forests in their ability to assimilate CO2 as a function of temperature were attributable to consistent differences in the temperature dependence of photosynthesis and respiration. This observation provides empirical validation of efforts to develop models of NEE from the first-principle relationships between photosynthetic and respiratory processes and climate. Our results also showed that models of seasonal dynamics in NEE that lack specific consideration of the temperature dependence of respiration and photosynthesis are likely to carry significant uncertainties. INDEX TERMS: 0315 Atmospheric Composition and Structure: Biosphere/atmosphere interactions; 1615 Global Change: Biogeochemical processes (4805); 3307 Meteorology and Atmospheric Dynamics: Boundary layer processes; 3322 Meteorology and Atmospheric Dynamics: Land/atmosphere interactions; 4806 Oceanography: Biological and Chemical: Carbon cycling. Citation: Yi, C., et al. (2004), A nonparametric method for separating photosynthesis and respiration components in CO2 flux measurements, Geophys. Res. Lett., 31, L17107, doi:10.1029/2004GL020490.

Focus on extreme events and the carbon cycle

Climate physics indicates that warming climate is a likely cause of extreme weather and more frequent and intense climate events. These extreme events can disrupt terrestrial carbon dynamics dramatically by triggering ecological disturbances and potentially forcing climate–carbon feedbacks. In this paper we synthesize the findings of 26 papers that focus on collecting evidence and developing knowledge of how extreme events disturb terrestrial carbon dynamics.

Climate extremes and grassland potential productivity

The considerable interannual variability (IAV) (?5 PgC yr?1) observed in atmospheric CO2 is dominated by variability in terrestrial productivity. Among terrestrial ecosystems, grassland productivity IAV is greatest. Relationships between grassland productivity IAV and climate drivers are poorly explained by traditional multiple-regression approaches. We propose a novel method, the perfect-deficit approach, to identify climate drivers of grassland IAV from observational data. The maximum daily value of each ecological or meteorological variable for each day of the year, over the period of record, defines the ?perfect? annual curve. Deficits of these variables can be identified by comparing daily observational data for a given year against the perfect curve. Links between large deficits of ecosystem activity and extreme climate events are readily identified. We applied this approach to five grassland sites with 26 site-years of observational data. Large deficits of canopy photosynthetic capacity and evapotranspiration derived from eddy-covariance measurements, and leaf area index derived from satellite data occur together and are driven by a local-dryness index during the growing season. This new method shows great promise in using observational evidence to demonstrate how extreme climate events alter yearly dynamics of ecosystem potential productivity and exchanges with atmosphere, and shine a new light on climate?carbon feedback mechanisms.