Jean-Marc Bonnefond

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.

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.

Preface paper to the Semi-Arid Land-Surface-Atmosphere (SALSA) Program special issue.

The Semi-Arid Land-Surface-Atmosphere Program (SALSA) is a multi-agency, multi-national research effort that seeks to evaluate the consequences of natural and human-induced environmental change in semi-arid regions. The ultimate goal of SALSA is to advance scientific understanding of the semi-arid portion of the hydrosphere-biosphere interface in order to provide reliable information for environmental decision making. SALSA approaches this goal through a program of long-term, integrated observations, process research, modeling, assessment, and information management that is sustained by cooperation among scientists and information users. In this preface to the SALSA special issue, general program background information and the critical nature of semi-arid regions is presented. A brief description of the Upper San Pedro River Basin, the initial location for focused SALSA research follows. Several overarching research objectives under which much of the interdisciplinary research contained in the special issue was undertaken are discussed. Principal methods, primary research sites and data collection used by numerous investigators during 1997-1999 are then presented. Scientists from about 20 US, five European (four French and one Dutch), and three Mexican agencies and institutions have collaborated closely to make the research leading to this special issue a reality. The SALSA Program has served as a model of interagency cooperation by breaking new ground in the approach to large scale interdisciplinary science with relatively limited resources.

Estimation of heat and momentum fluxes over complex terrain using a large aperture scintillometer.

A comprehensive experimental plan has been designed to further investigate the potential and the limitations associated with the use of a large aperture scintillometer (LAS) to infer path average sensible and momentum fluxes over complex surfaces as part of the Semi-Arid Land-Surface-Atmosphere (SALSA) Program. The complexity of the terrain is associated with the type and the cover of the vegetation canopy as well as with changes in topography. Scintillometer based estimates of sensible heat flux and friction velocity are compared to those measured by eddy correlation systems over a grassland patch, a mesquite patch, and over a transect spanning both patches. The results show that considering the complexity of the surface, the overall performance of the scintillometer is relatively good.

A long-term study of soil heat flux under a forest canopy

International programmes such as EUROFLUX focus on the analysis of long-term fluxes and energy budgets in the biosphere. Reliable estimates of hourly energy budgets require an accurate estimation of soil heat flux, that is often non-negligible even in a forest, and can be predominant during the night. Over long periods of time such as one to several months, its contribution can also be significant. The present work has been carried out to get good estimates of the soil heat flux in a maritime pine stand in the southwest of France, one of the 15 EUROFLUX sites. Using a whole year’s worth of data, soil heat flux was estimated by a two-step version of the null-alignment method using soil temperature, water content and bulk density measurements between the soil surface and a depth of 1 m. A data subset was firstly used to estimate and model the soil thermal conductivity at various depths. The full data set was then used with the modelled conductivity to estimate heat storage between the surface and a reference depth, and calculate the heat flux at the soil surface. Throughout the investigated year and at a 30 min time scale, the soil heat flux represents 5‐10% of the incident net radiation, i.e. 30‐50% of the net radiation over the understorey. Cumulative values from September 1997 to March 1998 reach a maximum of 70 MJ m 2 , which represents nearly 50% of the cumulative values of transmitted net radiation (140 MJ m 2 ) over the same period. These estimates of soil heat flux allowed the energy budgets of the whole stand and the understorey to be closed, and showed that the storage terms are significant not only at a 30 min time scale but also at longer time scales (a few weeks). An attempt was finally made to model soil heat flux from meteorological data, which has rarely been done for a forest soil and over a long-term data set. In most of the existing models, soil heat flux is taken as a fraction of net radiation or sensible heat flux. Here, the litter acts as a mulch at the soil surface so that the only significant terms of the energy balance at this level are soil heat flux, transmitted net radiation and turbulent sensible heat flux. Soil heat flux is shown to be a linear combination of (1) net radiation above the understorey with a clear dependence of the coefficient on the soil cover fraction, and (2) the difference between the air and litter temperatures, with little influence of soil water content or wind speed on the coefficient.

Production and carbon allocation in monocultures and mixed-species plantations of Eucalyptus grandis and Acacia mangium in Brazil

Introducing nitrogen-fixing tree species in fast-growing eucalypt plantations has the potential to improve soil nitrogen availability compared with eucalypt monocultures. Whether or not the changes in soil nutrient status and stand structure will lead to mixtures that out-yield monocultures depends on the balance between positive interactions and the negative effects of interspecific competition, and on their effect on carbon (C) uptake and partitioning. We used a C budget approach to quantify growth, C uptake and C partitioning in monocultures of Eucalyptus grandis (W. Hill ex Maiden) and Acacia mangium (Willd.) (treatments E100 and A100, respectively), and in a mixture at the same stocking density with the two species at a proportion of 1 : 1 (treatment MS). Allometric relationships established over the whole rotation, and measurements of soil CO(2) efflux and aboveground litterfall for ages 4-6 years after planting were used to estimate aboveground net primary production (ANPP), total belowground carbon flux (TBCF) and gross primary production (GPP). We tested the hypotheses that (i) species differences for wood production between E. grandis and A. mangium monocultures were partly explained by different C partitioning strategies, and (ii) the observed lower wood production in the mixture compared with eucalypt monoculture was mostly explained by a lower partitioning aboveground. At the end of the rotation, total aboveground biomass was lowest in A100 (10.5 kg DM m(-2)), intermediate in MS (12.2 kg DM m(-2)) and highest in E100 (13.9 kg DM m(-2)). The results did not support our first hypothesis of contrasting C partitioning strategies between E. grandis and A. mangium monocultures: the 21% lower growth (ΔB(w)) in A100 compared with E100 was almost entirely explained by a 23% lower GPP, with little or no species difference in ratios such as TBCF/GPP, ANPP/TBCF, ΔB(w)/ANPP and ΔB(w)/GPP. In contrast, the 28% lower ΔB(w) in MS than in E100 was explained both by a 15% lower GPP and by a 15% lower fraction of GPP allocated to wood growth, thus partially supporting our second hypothesis: mixing the two species led to shifts in C allocations from above- to belowground, and from growth to litter production, for both species.

Turbulent Structures in a Pine Forest with a Deep and Sparse Trunk Space: Stand and Edge Regions

Forested landscapes often exhibit large spatial variability in vertical and horizontal foliage distributions. This variability may affect canopy-atmosphere exchanges through its action on the development of turbulent structures. Here we investigate in neutral stratification the turbulent structures encountered in a maritime pine forest characterized by a high, dense foliated layer associated with a deep and sparse trunk space. Both stand and edge regions are considered. In situ measurements and the results of large-eddy simulations are used and analyzed together. In stand conditions, far from the edge, canopy-top structures appear strongly damped by the dense crown layer. Turbulent wind fluctuations within the trunk space, where the momentum flux vanishes, are closely related to these canopy-top structures through pressure diffusion. Consequently, autocorrelation and spectral analyses are not quite appropriate to characterize the vertical scale of coherent structures in this type of canopy, as pressure diffusion enhances the actual scale of structures. At frequencies higher than those associated with canopy-top structures, wind fluctuations related to wake structures developing behind tree stems are observed within the trunk space. They manifest themselves in wind velocity spectra as secondary peaks in the inertial subrange region, confirming the hypothesis of spectral short-cuts in vegetation canopies. In the edge region specific turbulent structures develop just below the crown layer, in addition to canopy-top structures. They are generated by the wind shear induced by the sub-canopy wind jet that forms at the edge. These structures provide a momentum exchange mechanism similar to that observed at the canopy top but in the opposite direction and with a lower magnitude. They may develop as in plane mixing-layer flows, with some perturbations induced by canopy-top structures. Wake structures are also observed within the trunk space in the edge region.

Integrated Sensible Heat Flux Measurements of a Two-Surface Composite Landscape using Scintillometry

The potential of the LAS (large aperture scintillometry) method for measuring sensible heat flux (H) directly integrated over a two-field composite surface is evaluated. We describe a field experiment performed within the Alpilles/ReSeDa project in the south-east of France over a composite surface made up of wheat and bare soil (451 and 216 m long respectively) using two 0.15-m diameter scintillometers mounted at heights of 2.05 and 4.54 m. When compared against reference values obtained by the eddy correlation technique, LAS-measured sensible heat flux reveals a systematic overestimation of about 10%. A simple model describing the integration of the scintillometer signal along the beam for a two-field composite surface is described. A simulation of the experiment confirms that the bias observed isrelated to non-linearities in the integration process in relation with thebell-shape sensitivity curve of the instrument to the structure parameter for the refractive index it measures. The model is used to test the sensitivity of the LAS-derived H values to the composition of the pathlength (ratio of both surfaces) and to the contrast in sensible heat flux and roughness length between the two fields. Sensitivity tests to the aggregation scheme for roughness length (two of them are tested) and to the measurement height are also presented. The composition of the surface in combination with the contrast in sensible heat flux (in direct relation with the contrast in latent heat flux) explains most of the bias, with possible deviations ranging from -50 up to 80 W m-2. A tentative semi-empirical method for correcting the bias is suggested, which only requires a crude estimate of the contrast in component sensible heat fluxes along the pathlength.

Paired comparison of water, energy and carbon exchanges over two young maritime pine stands (Pinus pinaster Ait.): effects of thinning and weeding in the early stage of tree growth.

The effects of management practices on energy, water and carbon exchanges were investigated in a young pine plantation in south-west France. In 2009-10, carbon dioxide (CO(2)), H(2)O and heat fluxes were monitored using the eddy covariance and sap flow techniques in a control plot (C) with a developed gorse layer, and an adjacent plot that was mechanically weeded and thinned (W). Despite large differences in the total leaf area index and canopy structure, the annual net radiation absorbed was only 4% lower in plot W. We showed that higher albedo in this plot was offset by lower emitted long-wave radiation. Annual evapotranspiration (ET) from plot W was 15% lower, due to lower rainfall interception and transpiration by the tree canopy, partly counterbalanced by the larger evaporation from both soil and regrowing weedy vegetation. The drainage belowground from plot W was larger by 113 mm annually. The seasonal variability of ET was driven by the dynamics of the soil and weed layers, which was more severely affected by drought in plot C. Conversely, the temporal changes in pine transpiration and stem diameter growth were synchronous between sites despite higher soil water content in the weeded plot. At the annual scale, both plots were carbon sinks, but thinning and weeding reduced the carbon uptake by 73%: annual carbon uptake was 243 and 65 g C m(-2) on plots C and W, respectively. Summer drought dramatically impacted the net ecosystem exchange: plot C became a carbon source as the gross primary production (GPP) severely decreased. However, plot W remained a carbon sink during drought, as a result of decreases in both GPP and ecosystem respiration (R(E)). In winter, both plots were carbon sources, plots C and W emitting 67.5 and 32.4 g C m(-2), respectively. Overall, this study highlighted the significant contribution of the gorse layer to mass and energy exchange in young pine plantations.

Analysis of the limits of the CT2-profile method for sensible heat flux measurements in unstable conditions

We present a test of the CT 2 -profile method described by Hill et al. [J. Atmos. Ocean. Technol. 9 (5) (1992) 526] to estimate the surface sensible heat flux over an homogeneous surface. A comparison with traditional eddy correlation measurements performed over a pasture (during the SALSA-Mexico experiment) using three identical large aperture scintillometers (LASs) along a 330 m propagation path and placed at heights 2.50, 3.45 and 6.45 m is first given. Scintillometer derived fluxes using the classical method at one level [Agric. For. Meteorol. 76 (1995) 149] reveal that the three scintillometers provide consistent measurements but underestimate by 15% the flux obtained with the 3D sonic anemometer. This is attributed to spatial non-homogeneities of the experimental site. Considerable scatter (and even the impossibility of performing computations) is found when using the CT 2 -profile method which is particularly prone to errors in nearly neutral and highly unstable conditions. The sensitivity of these errors to the accuracy of scintillometer measurements, the calibration errors and the measurement heights is investigated numerically. Simulations are made assuming a normal distribution of the relative error for CN 2 with standard deviations between 2 and 5% and no calibration error in a first step. Only calibration errors (up to 4% between instruments) are simulated in a second step. They confirm that the profile method degrades very rapidly with the accuracy of CN 2 : for instance the RMS error for H reaches 68 W m 2 (and the cases of impossible computation 28%) for a realistic D 5% value, with heights 2.50 and 3.45 m. Results appear slightly less sensitive to small calibration errors. The choice of the measurement heights z1 and z2 is also analysed: a ratio z2=z1 3 or 4 with z1 > 2 m seems the best compromise to minimise errors in H. Nevertheless the accuracy of the profile method is always much lower than that given by the classical method using measurements at one level, provided a good estimate of roughness length is available. We conclude that the CT 2 -profile method is not suitable for routine applications.