Search

Accurate representation of the turbulent exchange of carbon, water, and heat between the land surface and the atmosphere is critical for modelling global energy, water, and carbon cycles in both future climate projections and weather forecasts. Evaluation of models' ability to do this is performed in a wide range of simulation environments, often without explicit consideration of the degree of observational constraint or uncertainty and typically without quantification of benchmark performance expectations. We describe a Model Intercomparison Project (MIP) that attempts to resolve these shortcomings, comparing the surface turbulent heat flux predictions of around 20 different land models provided with in situ meteorological forcing evaluated with measured surface fluxes using quality-controlled data from 170 eddy-covariance-based flux tower sites. Predictions from seven out-of-sample empirical models are used to quantify the information available to land models in their forcing data and so the potential for land model performance improvement. Sites with unusual behaviour, complicated processes, poor data quality, or uncommon flux magnitude are more difficult to predict for both mechanistic and empirical models, providing a means of fairer assessment of land model performance. When examining observational uncertainty, model performance does not appear to improve in low-turbulence periods or with energy-balance-corrected flux tower data, and indeed some results raise questions about whether the energy balance correction process itself is appropriate. In all cases the results are broadly consistent, with simple out-of-sample empirical models, including linear regression, comfortably outperforming mechanistic land models. In all but two cases, latent heat flux and net ecosystem exchange of CO2 are better predicted by land models than sensible heat flux, despite it seeming to have fewer physical controlling processes. Land models that are implemented in Earth system models also appear to perform notably better than stand-alone ecosystem (including demographic) models, at least in terms of the fluxes examined here. The approach we outline enables isolation of the locations and conditions under which model developers can know that a land model can improve, allowing information pathways and discrete parameterisations in models to be identified and targeted for future model development.

The O2 : CO2 exchange ratio (ER) between terrestrial ecosystems and the atmosphere is a key parameter for partitioning global ocean and land carbon fluxes. The long-term terrestrial ER is considered to be close to 1.10 mol of O2 consumed per mole of CO2 produced. Due to the technical challenge in measuring directly the ER of entire terrestrial ecosystems (EReco), little is known about variations in ER at hourly and seasonal scales, as well as how different components contribute to EReco. In this modeling study, we explored the variability in and drivers of EReco and evaluated the hypothetical uncertainty in determining ecosystem O2 fluxes based on current instrument precision. We adapted the one-dimensional, multilayer atmosphere–biosphere gas exchange model “CANVEG” to simulate hourly EReco from modeled O2 and CO2 fluxes in a temperate beech forest in Germany. We found that the modeled annual mean EReco ranged from 1.06 to 1.12 mol mol−1 within the 5-year study period. Hourly EReco showed strong variations over diel and seasonal cycles and within the vertical canopy profile. The determination of ER from O2 and CO2 mole fractions in air above and within the canopy (ERconc) varied between 1.115 and 1.15 mol mol−1. CANVEG simulations also indicated that ecosystem O2 fluxes could be derived with the flux-gradient method using measured vertical gradients in scalar properties, as well as fluxes of CO2, sensible heat and latent energy derived from eddy covariance measurements. Owing to measurement uncertainties, however, the uncertainty in estimated O2 fluxes derived with the flux-gradient approach could be as high as 15 µmol m−2 s−1, which represented the 90 % quantile of the uncertainty in hourly data with a high-accuracy instrument. We also demonstrated that O2 fluxes can be used to partition net CO2 exchange fluxes into their component fluxes of photosynthesis and respiration if EReco is known. The uncertainty of the partitioned gross assimilation ranged from 1.43 to 4.88 µmol m−2 s−1 assuming a measurement uncertainty of 0.1 or 2.5 µmol m−2 s−1 for net ecosystem CO2 exchange and from 0.1 to 15 µmol m−2 s−1 for net ecosystem O2 exchange, respectively. Our analysis suggests that O2 measurements at ecosystem scale have the potential to partition net CO2 fluxes into their component fluxes, but further improvement in instrument precision is needed.

Information about forest background reflectance is needed for accurate biophysical parameter retrieval from forest canopies (overstory) with remote sensing. Separating under- and overstory signals would enable more accurate modeling of forest carbon and energy fluxes. We retrieved values of the normalized difference vegetation index (NDVI) of the forest understory with the multi-angular Moderate Resolution Imaging Spectroradiometer (MODIS) bidirectional reflectance distribution function (BRDF)/albedo data (gridded 500 m daily Collection 6 product), using a method originally developed for boreal forests. The forest floor background reflectance estimates from the MODIS data were compared with in situ understory reflectance measurements carried out at an extensive set of forest ecosystem experimental sites across Europe. The reflectance estimates from MODIS data were, hence, tested across diverse forest conditions and phenological phases during the growing season to examine their applicability for ecosystems other than boreal forests. Here we report that the method can deliver good retrievals, especially over different forest types with open canopies (low foliage cover). The performance of the method was found to be limited over forests with closed canopies (high foliage cover), where the signal from understory becomes too attenuated. The spatial heterogeneity of individual field sites and the limitations and documented quality of the MODIS BRDF product are shown to be important for the correct assessment and validation of the retrievals obtained with remote sensing.

The presence of vegetation alters hydrological cycles of ecosystems. Complex plant–soil interactions govern the fate of precipitation input and water transitions through ecosystem compartments. Disentangling these interactions is a major challenge in the field of ecohydrology and a pivotal foundation for understanding the carbon cycle of semi-arid ecosystems. Stable water isotopes can be used in this context as tracer to quantify water movement through soil–vegetation–atmosphere interfaces. The aim of this study is to disentangle vegetation effects on soil water infiltration and distribution as well as dynamics of soil evaporation and grassland water use in a Mediterranean cork oak woodland during dry conditions. An irrigation experiment using δ18O labelled water was carried out in order to quantify distinct effects of tree and herbaceous vegetation on the infiltration and distribution of event water in the soil profile. Dynamic responses of soil and herbaceous vegetation fluxes to precipitation regarding event water use, water uptake depth plasticity, and contribution to ecosystem soil evaporation and transpiration were quantified. Total water loss to the atmosphere from bare soil was as high as from vegetated soil, utilizing large amounts of unproductive evaporation for transpiration, but infiltration rates decreased. No adjustments of main root water uptake depth to changes in water availability could be observed during the experiment. This forces understorey plants to compete with adjacent trees for water in deeper soil layers at the onset of summer. Thus, understorey plants are subjected to chronic water deficits faster, leading to premature senescence at the onset of drought. Despite this water competition, the presence of cork oak trees fosters infiltration and reduces evapotranspirative water losses from the understorey and the soil, both due to altered microclimatic conditions under crown shading. This study highlights complex soil–plant–atmosphere and inter-species interactions controlling rain pulse transitions through a typical Mediterranean savannah ecosystem, disentangled by the use of stable water isotopes.

Savannah-type ecosystems account for 26–30% of global gross primary productivity GPP, with water being one of the major driving factors. In Europe, savannah-type woodlands cover an area of about 1.5 million ha. Here, the recent past has shown a significant decrease in precipitation P in winter and spring as well as a decrease in total annual precipitation. Strong effects on local water balance and carbon sink strength have thus been reported due to changes in precipitation regime. The objective of this study is to quantify the impact of the extreme drought event in 2012 on the water balance, gross primary productivity and carbon sink strength of a typical Portuguese cork-oak woodland (montado) compared to the wet year of 2011. Physiological responses of the dominant tree species Quercus suber (L.) are disentangled employing combined photosynthesis and stomatal conductance modelling. Precipitation effectiveness ET/P increased from 86% in 2011 to 122% in the 2012 dry year due to deep soil or groundwater access of the Q. suber trees leaving no water for groundwater replenishment. Understorey and overstorey GPP were strongly reduced, by 53 and 28%, respectively, in 2012 compared to 2011, due to the late onset of the autumn rains in 2011 and an additional severe winter/spring drought. However, the ecosystem was still a carbon sink in both years, but with a 38% reduced sink strength under extreme drought in 2012 compared to 2011. The combined photosynthesis–stomatal conductance model yielded the best results if it was allowed to adjust photosynthetic and stomatal parameters simultaneously. If stomatal response was modelled with the Leuning approach, which allows for a different sensitivity to vapour pressure deficit, the stomatal model parameters were highly coupled. A change in either of the parameters needed to be compensated by the other to guarantee a stable sensitivity of stomatal conductance to assimilation, independent of variations in vapour pressure deficit. The Q. suber trees showed a 37% reduced stomatal conductance during the drought period of 2012 compared to 2011, due to water supply limitations. In response to reduced leaf-internal CO2 availability, the trees strongly reduced the apparent maximum carboxylation rate by 43% in 2012 compared to 2011. Unexpectedly, the optimum temperature Topt of the maximum electron transport rate decreased during the drought period, enhancing the susceptibility of the trees to high temperature stress during the summer. Our results suggest that, if the trend of decreasing annual precipitation and changed precipitation patterns on the Iberian Peninsula continues, sustained effects on local groundwater reservoirs, understorey species composition and tree mortality have to be expected in the long term. To model the effect of drought on the montado ecosystem successfully, variable apparent maximum carboxylation rate Vc,max, stomatal conductance parameter m and vapour pressure deficit sensitivity parameter D0 need to be incorporated into photosynthesis–stomatal conductance modelling.

Stable isotope analysis is a powerful tool for assessing plant carbon and water relations and their impact on biogeochemical processes at different scales. Our process-based understanding of stable isotope signals, as well as technological developments, has progressed significantly, opening new frontiers in ecological and interdisciplinary research. This has promoted the broad utilisation of carbon, oxygen and hydrogen isotope applications to gain insight into plant carbon and water cycling and their interaction with the atmosphere and pedosphere. Here, we highlight specific areas of recent progress and new research challenges in plant carbon and water relations, using selected examples covering scales from the leaf to the regional scale. Further, we discuss strengths and limitations of recent technological developments and approaches and highlight new opportunities arising from unprecedented temporal and spatial resolution of stable isotope measurements.