Paul Berbigier

Energy balance closure at FLUXNET sites

A comprehensive evaluation of energy balance closure is performed across 22 sites and 50 site-years in FLUXNET, a network of eddy covariance sites measuring long-term carbon and energy fluxes in contrasting ecosystems and climates. Energy balance closure was evaluated by statistical regression of turbulent energy fluxes (sensible and latent heat (LE)) against available energy (net radiation, less the energy stored) and by solving for the energy balance ratio, the ratio of turbulent energy fluxes to available energy. These methods indicate a general lack of closure at most sites, with a mean imbalance in the order of 20%. The imbalance was prevalent in all measured vegetation types and in climates ranging from Mediterranean to temperate and arctic. There were no clear differences between sites using open and closed path infrared gas analyzers. At a majority of sites closure improved with turbulent intensity (friction velocity), but lack of total closure was still prevalent under most conditions. The imbalance was greatest during nocturnal periods. The results suggest that estimates of the scalar turbulent fluxes of sensible and LE are underestimated and/or that available energy is overestimated. The implications on interpreting long-term CO2 fluxes at FLUXNET sites depends on whether the imbalance results primarily from general errors associated

Estimates of the annual net carbon and water exchange of forests: The EUROFLUX methodology

Publisher Summary The chapter has described the measurement system and the procedure followed for the computation of the fluxes and the procedure of flux summation, including data gap filling strategy, night flux corrections and error estimation. It begins with the introduction of estimates of the annual net carbon and water exchange of forests using the EUROFLUX methodology. The chapter then provides us with the theory and moves on to discuss the eddy covariance system and its sonic anemometer, temperature fluctuation measurements, infrared gas analyser, air transport system, and tower instrumentation. Additional measurements are also given in the chapter. Data acquisition and its computation and correction is discussed next in the chapter by giving its general procedure, half-hourly means (co-)variances and uncorrected fluxes, intercomparison of software, and correction for frequency response losses. The chapter has also discussed about quality control and four criteria are investigated here for the same. Spatial representativeness of measured fluxes and summation procedure are reviewed. The chapter then moves on to the discussion of data gap filling through interpolation and parameterization and neural networks. Corrections to night-time data and error estimation are also explored in the chapter. Finally, the chapter closes with conclusions.

Respiration as the main determinant of carbon balance in European forests

Carbon exchange between the terrestrial biosphere and the atmosphere is one of the key processes that need to be assessed in the context of the Kyoto Protocol. Several studies suggest that the terrestrial biosphere is gaining carbon, but these estimates are obtained primarily by indirect methods, and the factors that control terrestrial carbon exchange, its magnitude and primary locations, are under debate. Here we present data of net ecosystem carbon exchange, collected between 1996 and 1998 from 15 European forests, which confirm that many European forest ecosystems act as carbon sinks. The annual carbon balances range from an uptake of 6.6 tonnes of carbon per hectare per year to a release of nearly 1 t C ha -1 yr-1, with a large variability between forests. The data show a significant increase of carbon uptake with decreasing latitude, whereas the gross primary production seems to be largely independent of latitude. Our observations indicate that, in general, ecosystem respiration determines net ecosystem carbon exchange. Also, for an accurate assessment of the carbon balance in a particular forest ecosystem, remote sensing of the normalized difference vegetation index or estimates based on forest inventories may not be sufficient.

The human footprint in the carbon cycle of temperate and boreal forests

Temperate and boreal forests in the Northern Hemisphere cover an area of about 2 × 107 square kilometres and act as a substantial carbon sink (0.6–0.7 petagrams of carbon per year). Although forest expansion following agricultural abandonment is certainly responsible for an important fraction of this carbon sink activity, the additional effects on the carbon balance of established forests of increased atmospheric carbon dioxide, increasing temperatures, changes in management practices and nitrogen deposition are difficult to disentangle, despite an extensive network of measurement stations. The relevance of this measurement effort has also been questioned, because spot measurements fail to take into account the role of disturbances, either natural (fire, pests, windstorms) or anthropogenic (forest harvesting). Here we show that the temporal dynamics following stand-replacing disturbances do indeed account for a very large fraction of the overall variability in forest carbon sequestration. After the confounding effects of disturbance have been factored out, however, forest net carbon sequestration is found to be overwhelmingly driven by nitrogen deposition, largely the result of anthropogenic activities. The effect is always positive over the range of nitrogen deposition covered by currently available data sets, casting doubts on the risk of widespread ecosystem nitrogen saturation under natural conditions. The results demonstrate that mankind is ultimately controlling the carbon balance of temperate and boreal forests, either directly (through forest management) or indirectly (through nitrogen deposition).

CO2 and water vapour fluxes for 2 years above Euroflux forest site

The ‘Le Bray’ site, a 28-year old plantation of maritime pine (Pinus pinaster Ait.) with an understorey of graminae, is part of the European network Euroflux. Carbon dioxide and water vapour fluxes were measured continuously for 2 years above the canopy with an eddy-covariance system. The 2-year accumulated evapotranspiration amounts to 3184 ± 440 MJ m −2 (1331 ± 186 mm of water) compared to 1860 mm rainfall or 6011 MJ m −2 net radiation. The overall gap in the energy balance is 1322 MJ m −2 , partly due to the poor estimation of evapotranspiration during rainy periods. The 2-year sum of carbon sequestration is 11.5 ± 0.8 t of carbon per hectare. The 2-year total ecosystem respiration, estimated by adjusting statistically night-time CO2 fluxes against the soil–litter interface temperature for turbulent nights only, is about 33.6 ± 2.0 t of carbon. A fairly good statistical relationship is found by multiple regression between daily gross primary production (GPP) on the one hand, and direct and diffuse absorbed PAR on the other hand. The water use efficiency (defined as the ratio of GPP to evapotranspiration) is inversely related to air saturation deficit.

Energy Partitioning Between Latent And Sensible Heat Flux During The Warm Season At Fluxnet Sites

The warm season (mid-June through late August) partitioning between sensible (H) and latent (LE) heat flux, or the Bowen ratio (beta=H/LE), was investigated at 27 sites over 66 site years within the international network of eddy covariance sites (FLUXNET). Variability in beta across ecosystems and climates was analyzed by quantifying general climatic and surface characteristics that control flux partitioning. The climatic control on beta was quantified using the climatological resistance (R-i), which is proportional to the ratio of vapor pressure deficit (difference between saturation vapor pressure and atmospheric vapor pressure) to net radiation (large values of R-i decrease beta). The control of flux partitioning by the vegetation and underlying surface was quantified by computing the surface resistance to water vapor transport (R-c, with large values tending to increase beta). There was a considerable range in flux partitioning characteristics (R-c, R-i and beta) among sites, but it was possible to define some general differences between vegetation types and climates. Deciduous forest sites and the agricultural site had the lowest values of R-c and beta (0.25-0.50). Coniferous forests typically had a larger R-c and higher beta (typically between 0.50 and 1.00 but also much larger). However, there was notable variability in R-c and R-i between coniferous site years, most notably differences between oceanic and continental climates and sites with a distinct dry summer season (Mediterranean climate). Sites with Mediterranean climates generally had the highest net radiation, R-c, R-i, and beta. There was considerable variability in beta between grassland site years, primarily the result of interannual differences in soil water content and R-c

Phase and amplitude of ecosystem carbon release and uptake potentials as derived from FLUXNET measurements

As length and timing of the growing season are major factors explaining differences in carbon exchange of ecosystems, we analyzed seasonal patterns of net ecosystem carbon exchange (FNEE) using eddy covariance data of the FLUXNET data base (http://www-eosdis.ornl.gov/FLUXNET). The study included boreal and temperate, deciduous and coniferous forests, Mediterranean evergreen systems, rainforest, native and managed temperate grasslands, tundra, and C3 and C4 crops. Generalization of seasonal patterns are useful for identifying functional vegetation types for global dynamic vegetation models, as well as for global inversion studies, and can help improve phenological modules in SVAT or biogeochemical models. The results of this study have important validation potential for global carbon cycle modeling. The phasing of respiratory and assimilatory capacity differed within forest types: for temperate coniferous forests seasonal uptake and release capacities are in phase, for temperate deciduous and boreal coniferous forests, release was delayed compared to uptake. According to seasonal pattern of maximum nighttime release (evaluated over 15-day periods, Fmax) the study sites can be grouped in four classes: (1) boreal and high altitude conifers and grasslands; (2) temperate deciduous and temperate conifers; (3) tundra and crops; (4) evergreen Mediterranean and tropical forests. Similar results are found for maximum daytime uptake (Fmin) and the integral net carbon flux, but temperate deciduous forests fall into class 1. For forests, seasonal amplitudes of Fmax and Fmin increased in the order tropical C3-crops>temperate deciduous forests>temperate conifers>boreal conifers>tundra ecosystems. Due to data restrictions, our analysis centered mainly on Northern Hemisphere temperate and boreal forest ecosystems. Grasslands, crops, Mediterranean ecosystems, and rainforests are under-represented, as are savanna systems, wooded grassland, shrubland, or year-round measurements in tundra systems. For regional or global estimates of carbon sequestration potentials, future investigations of eddy covariance should expand in these systems.

Interception loss, throughfall and stemflow in a maritime pine stand. I. Variability of throughfall and stemflow beneath the pine canopy

Patterns of spatial variability of throughfall and stemflow were determined in a maritime pine (Pinus pinaster Ait.) stand for two consecutive years. Data were obtained from 52 fixed rain gauges and 12 stemflow measuring devices located in a 50m × 50m plot at the centre of an 18-year-old stand. The pine trees had been sown in rows 4m apart and had reached an average height of 12.6m. The spatial distribution of stems had a negligible effect on the throughfall partitioning beneath the canopy. Variograms of throughfall computed for a sample of storms did not reveal any spatial autocorrelation of throughfall for the sampling design used. Differences in throughfall, in relation to the distance from the rows, were not consistently significant. In addition, the distance from the tree stem did not influence the amount of throughfall. The confidence interval on the amount of throughfall per storm was between 3 and 8%. The stemflow was highly variable between trees. The effect of individual trees on stemflow was significant but the amount of stemflow per tree was not related to tree size (i.e. height, trunk diameter, etc.). The cumulative sampling errors on stemflow and throughfall for a single storm created a confidence interval of between ±7 and ±51% on interception. This resulted mainly from the low interception rate and sampling error on throughfall.

Implications of the carbon cycle steady state assumption for biogeochemical modeling performance and inverse parameter retrieval

We analyze the impacts of the steady state assumption on inverse model parameter retrieval from biogeochemical models. An inverse model parameterization study using eddy covariance CO 2 flux data was performed with the Carnegie Ames Stanford Approach (CASA) model under conditions of strict and relaxed carbon cycle steady state assumption (CCSSA) in order to evaluate both the robustness of the models structure for the simulation of net ecosystem carbon fluxes and the assessment of the CCSSA effects on simulations and parameter estimation. Net ecosystem production (NEP) measurements from several eddy covariance sites were compared with NEP estimates from the CASA model driven by local weather station climate inputs as well as by remotely sensed fraction of photosynthetically active radiation absorbed by vegetation and leaf area index. The parameters considered for optimization are directly related to aboveground and belowground modeled responses to temperature and water availability, as well as a parameter (η) that relaxed the CCSSA in the model, allowing for site level simulations to be initialized either as net sinks or sources. A robust relationship was observed between NEP observations and predictions for most of the sites through the range of temporal scales considered (daily, weekly, biweekly, and monthly), supporting the conclusion that the model structure is able to capture the main processes explaining NEP variability. Overall, relaxing CCSSA increased model efficiency (21%) and decreased normalized average error (-92%). Intersite variability was a major source of variance in model performance differences between fixed (CCSSA f ) and relaxed (CCSSA r ) CCSSA conditions. These differences were correlated with mean annual NEP observations, where an average increase in modeling efficiency of 0.06 per 100 g Cm -2 a -1 (where a is years) of NEP is observed (α < 0.003). The parameter η was found to be a key parameter in the optimization exercise, generating significant model efficiency losses when removed from the initial parameter set and parameter uncertainties were significantly lower under CCSSA r . Moreover, modeled soil carbon stocks were generally closer to observations once the steady state assumption was relaxed. Finally, we also show that estimates of individual parameters are affected by the steady state assumption. For example, estimates of radiation-use efficiency were strongly affected by the CCSSA f indicating compensation effects for the inadequate steady state assumption, leading to effective and thus biased parameters. Overall, the importance of model structural evaluation in data assimilation approaches is thus emphasized.

Partitioning net ecosystem carbon exchange into net assimilation and respiration with canopy‐scale isotopic measurements: An error propagation analysis with 13CO2 and CO18O data

13 C and 18 O isotopic signature of FR (dR), and for both isotopes the a priori uncertainty associated with this term is estimated to be around 2% at our site. Using d 13 C-CO2 and [CO2] measurements, we then show that the uncertainty on instantaneous values of FA and FR can be as large as 4 mmol m � 2 s � 1 . Even if we could get more accurate estimates of the net CO2 flux, the isoflux, and the isotopic signatures of FA and FR ,t his uncertainty would not be significantly reduced because the isotopic disequilibrium between FA and FR is too small, around 2–3%. With d 18 O-CO2 and [CO2] measurements the uncertainty associated with the gross fluxes lies also around 4 mmol m � 2 s � 1 but could be dramatically reduced if we were able to get more accurate estimates of the CO 18 O isoflux and the associated discrimination during photosynthesis. This is because the isotopic disequilibrium between FA and FR is large, of the order of 12–17%. The isotopic disequilibrium between FA and FR and the uncertainty on dR vary among ecosystems and over the year. Our approach should help to choose the best strategy to study the carbon budget of a given ecosystem using stable isotopes. INDEX TERMS: 4806 Oceanography: Biological and Chemical: Carbon cycling; 0315 Atmospheric Composition and Structure: Biosphere/atmosphere interactions; 1615 Global Change: Biogeochemical processes (4805); KEYWORDS: carbon cycle, carbon 13, oxygen 18, CO2 assimilation, respiration