Jérôme Ogée

Europe-wide reduction in primary productivity caused by the heat and drought in 2003

Future climate warming is expected to enhance plant growth in temperate ecosystems and to increase carbon sequestration. But although severe regional heatwaves may become more frequent in a changing climate, their impact on terrestrial carbon cycling is unclear. Here we report measurements of ecosystem carbon dioxide fluxes, remotely sensed radiation absorbed by plants, and country-level crop yields taken during the European heatwave in 2003. We use a terrestrial biosphere simulation model to assess continental-scale changes in primary productivity during 2003, and their consequences for the net carbon balance. We estimate a 30 per cent reduction in gross primary productivity over Europe, which resulted in a strong anomalous net source of carbon dioxide (0.5 Pg C yr-1) to the atmosphere and reversed the effect of four years of net ecosystem carbon sequestration. Our results suggest that productivity reduction in eastern and western Europe can be explained by rainfall deficit and extreme summer heat, respectively. We also find that ecosystem respiration decreased together with gross primary productivity, rather than accelerating with the temperature rise. Model results, corroborated by historical records of crop yields, suggest that such a reduction in Europes primary productivity is unprecedented during the last century. An increase in future drought events could turn temperate ecosystems into carbon sources, contributing to positive carbon-climate feedbacks already anticipated in the tropics and at high latitudes.

A dynamic global vegetation model for studies of the coupled atmosphere‐biosphere system

This work presents a new dynamic global vegetation model designed as an extension of an existing surface-vegetation-atmosphere transfer scheme which is included in a coupled ocean-atmosphere general circulation model. The new dynamic global vegetation model simulates the principal processes of the continental biosphere influencing the global carbon cycle (photosynthesis, autotrophic and heterotrophic respiration of plants and in soils, fire, etc.) as well as latent, sensible, and kinetic energy exchanges at the surface of soils and plants. As a dynamic vegetation model, it explicitly represents competitive processes such as light competition, sapling establishment, etc. It can thus be used in simulations for the study of feedbacks between transient climate and vegetation cover changes, but it can also be used with a prescribed vegetation distribution. The whole seasonal phenological cycle is prognostically calculated without any prescribed dates or use of satellite data. The model is coupled to the IPSL-CM4 coupled atmosphere-ocean-vegetation model. Carbon and surface energy fluxes from the coupled hydrology-vegetation model compare well with observations at FluxNet sites. Simulated vegetation distribution and leaf density in a global simulation are evaluated against observations, and carbon stocks and fluxes are compared to available estimates, with satisfying results.

Photosynthetic carbon isotope discrimination and its relationship to the carbon isotope signals of stem, soil and ecosystem respiration

• Photosynthetic carbon (C) isotope discrimination (Δ(Α)) labels photosynthates (δ(A) ) and atmospheric CO(2) (δ(a)) with variable C isotope compositions during fluctuating environmental conditions. In this context, the C isotope composition of respired CO(2) within ecosystems is often hypothesized to vary temporally with Δ(Α). • We investigated the relationship between Δ(Α) and the C isotope signals from stem (δ(W)), soil (δ(S)) and ecosystem (δ(E)) respired CO(2) to environmental fluctuations, using novel tuneable diode laser absorption spectrometer instrumentation in a mature maritime pine forest. • Broad seasonal changes in Δ(Α) were reflected in δ(W,) δ(S) and δ(E). However, respired CO(2) signals had smaller short-term variations than Δ(A) and were offset and delayed by 2-10 d, indicating fractionation and isotopic mixing in a large C pool. Variations in δ(S) did not follow Δ(A) at all times, especially during rainy periods and when there is a strong demand for C allocation above ground. • It is likely that future isotope-enabled vegetation models will need to develop transfer functions that can account for these phenomena in order to interpret and predict the isotopic impact of biosphere gas exchange on the C isotope composition of atmospheric CO(2).

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

A single‐substrate model to interpret intra‐annual stable isotope signals in tree‐ring cellulose

The carbon and oxygen stable isotope composition of wood cellulose (delta(13)C(cellulose) and delta(18)O(cellulose), respectively) reveal well-defined seasonal variations that contain valuable records of past climate, leaf gas exchange and carbon allocation dynamics within the trees. Here, we present a single-substrate model for wood growth to interpret seasonal isotopic signals collected in an even-aged maritime pine plantation growing in South-west France, where climate, soil and flux variables were also monitored. Observed seasonal patterns in delta(13)C(cellulose) and delta(18)O(cellulose) were different between years and individuals, and mostly captured by the model, suggesting that the single-substrate hypothesis is a good approximation for tree ring studies on Pinus pinaster, at least for the environmental conditions covered by this study. A sensitivity analysis revealed that the model was mostly affected by five isotopic discrimination factors and two leaf gas-exchange parameters. Modelled early wood signals were also very sensitive to the date when cell wall thickening begins (t(wt)). Our model could therefore be used to reconstruct t(wt) time series and improve our understanding of how climate influences this key parameter of xylogenesis.

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.

The impact of soil microorganisms on the global budget of δ18O in atmospheric CO2

Improved global estimates of terrestrial photosynthesis and respiration are critical for predicting the rate of change in atmospheric CO2. The oxygen isotopic composition of atmospheric CO2 can be used to estimate these fluxes because oxygen isotopic exchange between CO2 and water creates distinct isotopic flux signatures. The enzyme carbonic anhydrase (CA) is known to accelerate this exchange in leaves, but the possibility of CA activity in soils is commonly neglected. Here, we report widespread accelerated soil CO2 hydration. Exchange was 10–300 times faster than the uncatalyzed rate, consistent with typical population sizes for CA-containing soil microorganisms. Including accelerated soil hydration in global model simulations modifies contributions from soil and foliage to the global CO18O budget and eliminates persistent discrepancies existing between model and atmospheric observations. This enhanced soil hydration also increases the differences between the isotopic signatures of photosynthesis and respiration, particularly in the tropics, increasing the precision of CO2 gross fluxes obtained by using the δ18O of atmospheric CO2 by 50%.

Interactive effects of nocturnal transpiration and climate change on the root hydraulic redistribution and carbon and water budgets of southern United States pine plantations

Deep root water uptake and hydraulic redistribution (HR) have been shown to play a major role in forest ecosystems during drought, but little is known about the impact of climate change, fertilization and soil characteristics on HR and its consequences on water and carbon fluxes. Using data from three mid-rotation loblolly pine plantations, and simulations with the process-based model MuSICA, this study indicated that HR can mitigate the effects of soil drying and had important implications for carbon uptake potential and net ecosystem exchange (NEE), especially when N fertilization is considered. At the coastal site (C), characterized by deep organic soil, HR increased dry season tree transpiration (T) by up to 40%, and such an increase affected NEE through major changes in gross primary productivity (GPP). Deep-rooted trees did not necessarily translate into a large volume of HR unless soil texture allowed large water potential gradients to occur, as was the case at the sandy site (S). At the Piedmont site (P) characterized by a shallow clay-loam soil, HR was low but not negligible, representing up to 10% of T. In the absence of HR, it was predicted that at the C, S and P sites, annual GPP would have been diminished by 19, 7 and 9%, respectively. Under future climate conditions HR was predicted to be reduced by up to 25% at the C site, reducing the resilience of trees to precipitation deficits. The effect of HR on T and GPP was predicted to diminish under future conditions by 12 and 6% at the C and P sites, respectively. Under future conditions, T was predicted to stay the same at the P site, but to be marginally reduced at the C site and slightly increased at the S site. Future conditions and N fertilization would decrease T by 25% at the C site, by 15% at the P site and by 8% at the S site. At the C and S sites, GPP was estimated to increase by 18% and by >70% under future conditions, respectively, with little effect of N fertilization. At the P site, future conditions would stimulate GPP by only 12%, but future conditions plus N fertilization would increase GPP by 24%. As a consequence, in all sites, water use efficiency was predicted to improve dramatically with future conditions. Modeling the effect of reduced annual precipitation indicated that limited water availability would decrease all carbon fluxes, including NEE and respiration. Our simulations highlight the interactive effects of nutrients and elevated CO(2), and showed that the effect of N fertilization would be greater under future climate conditions.

Validation of eddy flux measurements above the understorey of a pine forest

Measurements of turbulent exchange in the understorey of a forest canopy are necessary to understand and model both the functioning of the lower stratum and its contribution to turbulent exchange at the whole canopy scale. Eddy covariance measurements of eddy fluxes just above the floor of a forest canopy must be thoroughly validated, given the particular conditions prevailing there. Our objective is two-fold: (i) check the overall quality of such eddy flux measurements through the analysis of the understorey energy balance closure and (ii) define quality criteria for each half-hourly sample, based on the residual term of the energy balance. A subset of the EUROFLUX data base, collected within a pine forest canopy in south-west France, was used for this purpose. During this experiment, all heat storage terms were carefully measured, which allowed the closure of the understorey energy balance to be rigorously tested. As in most experiments storage term measurements are not available, we also developed a method to estimate them, in order to apply the above-mentioned data selection method. The energy balance closure was found to be quite satisfactory (the slope of the sum of eddy fluxes and storage terms versus transmitted net radiation is 0.99, the intercept is less than 1 W m 2 , r 2 is 0.94, there is no deviation from a linear trend). The data selection procedure allows a fair description of the daily and day-to-day variation of turbulent fluxes while rejecting the most dubious data, whether experimental or estimated storage terms are used. This analysis proves the validity of eddy flux measurements in the lower part of the forest and offers tools for flux data selection, depending on the type of studies such data are intended for.

Improvement of water and light availability after thinning at a xeric site: which matters more? A dual isotope approach

Thinning fosters individual tree growth by increasing the availability of water, light and nutrients. At sites where water rather than light is limiting, thinning also enhances soil evaporation and might not be beneficial. Detailed knowledge of the short- to long-term physiological response underlying the growth responses to thinning is crucial for the management of forests already suffering from recurrent drought-induced dieback. We applied a dual isotope approach together with mechanistic isotope models to study the physiological processes underlying long-term growth enhancement of heavily thinned Pinus sylvestris in a xeric forest in Switzerland. This approach allowed us to identify and disentangle thinning-induced changes in stomatal conductance and assimilation rate. At our xeric study site, the increase in stomatal conductance far outweighed the increase in assimilation, implying that growth release in heavily thinned trees is primarily driven by enhanced water availability rather than increased light availability. We conclude that in forests with relatively isohydric species (drought avoiders) that are growing close to their physiological limits, thinning is recommended to maintain a less negative water balance and thus foster tree growth, and ultimately the survival of forest trees under drought.