Leaf-level coordination principles propagate to the ecosystem scale.

Fundamental axes of variation in plant traits result from trade-offs between costs and benefits of resource-use strategies at the leaf scale. However, it is unclear whether similar trade-offs propagate to the ecosystem level. Here, we test whether trait correlation patterns predicted by three well-known leaf- and plant-level coordination theories - the leaf economics spectrum, the global spectrum of plant form and function, and the least-cost hypothesis - are also observed between community mean traits and ecosystem processes. We combined ecosystem functional properties from FLUXNET sites, vegetation properties, and community mean plant traits into three corresponding principal component analyses. We find that the leaf economics spectrum (90 sites), the global spectrum of plant form and function (89 sites), and the least-cost hypothesis (82 sites) all propagate at the ecosystem level. However, we also find evidence of additional scale-emergent properties. Evaluating the coordination of ecosystem functional properties may aid the development of more realistic global dynamic vegetation models with critical empirical data, reducing the uncertainty of climate change projections.

ForestClim-Bioclimatic variables for microclimate temperatures of European forests.

Microclimate research gained renewed interest over the last decade and its importance for many ecological processes is increasingly being recognized. Consequently, the call for high-resolution microclimatic temperature grids across broad spatial extents is becoming more pressing to improve ecological models. Here, we provide a new set of open-access bioclimatic variables for microclimate temperatures of European forests at 25 × 25 m2 resolution.

ForestTemp - Sub-canopy microclimate temperatures of European forests.

Ecological research heavily relies on coarse-gridded climate data based on standardized temperature measurements recorded at 2 m height in open landscapes. However, many organisms experience environmental conditions that differ substantially from those captured by these macroclimatic (i.e. free air) temperature grids. In forests, the tree canopy functions as a thermal insulator and buffers sub-canopy microclimatic conditions, thereby affecting biological and ecological processes. To improve the assessment of climatic conditions and climate-change-related impacts on forest-floor biodiversity and functioning, high-resolution temperature grids reflecting forest microclimates are thus urgently needed. Combining more than 1200 time series of in situ near-surface forest temperature with topographical, biological and macroclimatic variables in a machine learning model, we predicted the mean monthly offset between sub-canopy temperature at 15 cm above the surface and free-air temperature over the period 2000-2020 at a spatial resolution of 25 m across Europe. This offset was used to evaluate the difference between microclimate and macroclimate across space and seasons and finally enabled us to calculate mean annual and monthly temperatures for European forest understories. We found that sub-canopy air temperatures differ substantially from free-air temperatures, being on average 2.1°C (standard deviation ± 1.6°C) lower in summer and 2.0°C higher (±0.7°C) in winter across Europe. Additionally, our high-resolution maps expose considerable microclimatic variation within landscapes, not captured by the gridded macroclimatic products. The provided forest sub-canopy temperature maps will enable future research to model below-canopy biological processes and patterns, as well as species distributions more accurately.

The three major axes of terrestrial ecosystem function.

The leaf economics spectrum1,2 and the global spectrum of plant forms and functions3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species2. Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities4. However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability4,5. Here we derive a set of ecosystem functions6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems7,8.

Soil respiration variation along an altitudinal gradient in the Italian Alps: Disentangling forest structure and temperature effects.

On the mountains, along an elevation gradient, we generally observe an ample variation in temperature, with the associated difference in vegetation structure and composition and soil properties. With the aim of quantifying the relative importance of temperature, vegetation and edaphic properties on soil respiration (SR), we investigated changes in SR along an elevation gradient (404 to 2101 m a.s.l) in the southern slopes of the Alps in Northern Italy. We also analysed soil physicochemical properties, including soil organic carbon (SOC) and nitrogen (N) stocks, fine root C and N, litter C and N, soil bulk densities and soil pH at five forest sites, and also stand structural properties, including vegetation height, age and basal area. Our results indicated that SR rates increased with temperature in all sites, and 55-76% of SR variability was explained by temperature. Annual cumulative SR, ranging between 0.65-1.40 kg C m-2 yr-1, decreased along the elevation gradient, while temperature sensitivity (Q10) of SR increased with elevation. However, a high SR rate (1.27 kg C m-2 yr-1) and low Q10 were recorded in the mature conifer forest stand at 1731 m a.s.l., characterized by an uneven-aged structure and high dominant tree height, resulting in a nonlinear relationship between elevation and temperature. Reference SR at 10°C (SRref) was unrelated to elevation, but was related to tree height. A significant negative linear relationship was found between bulk density and elevation. Conversely, SOC, root C and N stock, pH, and litter mass were best fitted by nonlinear relationships with elevation. However, these parameters were not significantly correlated with SR when the effect of temperature was removed (SRref). These results demonstrate that the main factor affecting SR in forest ecosystems along this Alpine elevation gradient is temperature, but its regulating role can be strongly influenced by site biological characteristics, particularly vegetation type and structure, affecting litter quality and microclimate. This study also confirms that high elevation sites are rich in SOC and more sensitive to climate change, being prone to high C losses as CO2. Furthermore, our data indicate a positive relationship between Q10 and dominant tree height, suggesting that mature forest ecosystems characterized by an uneven-age structure, high SRref and moderate Q10, may be more resilient.

Temperature thresholds of ecosystem respiration at a global scale.

Ecosystem respiration is a major component of the global terrestrial carbon cycle and is strongly influenced by temperature. The global extent of the temperature-ecosystem respiration relationship, however, has not been fully explored. Here, we test linear and threshold models of ecosystem respiration across 210 globally distributed eddy covariance sites over an extensive temperature range. We find thresholds to the global temperature-ecosystem respiration relationship at high and low air temperatures and mid soil temperatures, which represent transitions in the temperature dependence and sensitivity of ecosystem respiration. Annual ecosystem respiration rates show a markedly reduced temperature dependence and sensitivity compared to half-hourly rates, and a single mid-temperature threshold for both air and soil temperature. Our study indicates a distinction in the influence of environmental factors, including temperature, on ecosystem respiration between latitudinal and climate gradients at short (half-hourly) and long (annual) timescales. Such climatological differences in the temperature sensitivity of ecosystem respiration have important consequences for the terrestrial net carbon sink under ongoing climate change.

Ecosystem transpiration and evaporation: Insights from three water flux partitioning methods across FLUXNET sites.

We apply and compare three widely applicable methods for estimating ecosystem transpiration (T) from eddy covariance (EC) data across 251 FLUXNET sites globally. All three methods are based on the coupled water and carbon relationship, but they differ in assumptions and parameterizations. Intercomparison of the three daily T estimates shows high correlation among methods (R between .89 and .94), but a spread in magnitudes of T/ET (evapotranspiration) from 45% to 77%. When compared at six sites with concurrent EC and sap flow measurements, all three EC-based T estimates show higher correlation to sap flow-based T than EC-based ET. The partitioning methods show expected tendencies of T/ET increasing with dryness (vapor pressure deficit and days since rain) and with leaf area index (LAI). Analysis of 140 sites with high-quality estimates for at least two continuous years shows that T/ET variability was 1.6 times higher across sites than across years. Spatial variability of T/ET was primarily driven by vegetation and soil characteristics (e.g., crop or grass designation, minimum annual LAI, soil coarse fragment volume) rather than climatic variables such as mean/standard deviation of temperature or precipitation. Overall, T and T/ET patterns are plausible and qualitatively consistent among the different water flux partitioning methods implying a significant advance made for estimating and understanding T globally, while the magnitudes remain uncertain. Our results represent the first extensive EC data-based estimates of ecosystem T permitting a data-driven perspective on the role of plants' water use for global water and carbon cycling in a changing climate.

Stomatal optimization based on xylem hydraulics (SOX) improves land surface model simulation of vegetation responses to climate.

Land surface models (LSMs) typically use empirical functions to represent vegetation responses to soil drought. These functions largely neglect recent advances in plant ecophysiology that link xylem hydraulic functioning with stomatal responses to climate. We developed an analytical stomatal optimization model based on xylem hydraulics (SOX) to predict plant responses to drought. Coupling SOX to the Joint UK Land Environment Simulator (JULES) LSM, we conducted a global evaluation of SOX against leaf- and ecosystem-level observations. SOX simulates leaf stomatal conductance responses to climate for woody plants more accurately and parsimoniously than the existing JULES stomatal conductance model. An ecosystem-level evaluation at 70 eddy flux sites shows that SOX decreases the sensitivity of gross primary productivity (GPP) to soil moisture, which improves the model agreement with observations and increases the predicted annual GPP by 30% in relation to JULES. SOX decreases JULES root-mean-square error in GPP by up to 45% in evergreen tropical forests, and can simulate realistic patterns of canopy water potential and soil water dynamics at the studied sites. SOX provides a parsimonious way to incorporate recent advances in plant hydraulics and optimality theory into LSMs, and an alternative to empirical stress factors.

Evidence for a non-linear carbon accumulation pattern along an Alpine glacier retreat chronosequence in Northern Italy.

BACKGROUND: The glaciers in the Alps, as in other high mountain ranges and boreal zones, are generally retreating and leaving a wide surface of bare ground free from ice cover. This early stage soil is then colonized by microbes and vegetation in a process of primary succession. It is rarely experimentally examined whether this colonization process is linear or not at the ecosystem scale. Thus, to improve our understanding of the variables involved in the carbon accumulation in the different stages of primary succession, we conducted this research in three transects on the Matsch glacier forefield (Alps, N Italy) at an altitude between 2,350 and 2,800 m a.s.l. METHODS: In three field campaigns (July, August and September 2014) a closed transparent chamber was used to quantify the net ecosystem exchange (NEE) between the natural vegetation and the atmosphere. On the five plots established in each of the three transects, shading nets were used to determine ecosystem response function to variable light conditions. Ecosystem respiration (Reco) and gross ecosystem exchange (GEE) was partitioned from NEE. Following the final flux measurements, biometric sampling was conducted to establish soil carbon (C) and nitrogen (N) content and the biomass components for each transect. RESULTS: A clear difference was found between the earlier and the later successional stage. The older successional stages in the lower altitudes acted as a stronger C sink, where NEE, GEE and Reco were significantly higher than in the earlier successional stage. Of the two lower transects, the sink capacity of intermediate-succession plots exceeded that of the plots of older formation, in spite of the more developed soil. Total biomass (above- and belowground) approached its maximum value in the intermediate ecosystem, whilst the later stage of succession predominated in the corresponding belowground organic mass (biomass, N and C). OUTLOOK: We found that the process of carbon accumulation along a glacier retreat chronosequence is not linear, and after a quite rapid increase in carbon accumulation capacity in the first 150 years, in average 9 g C m-2 year-1, it slows down, taking place mainly in the belowground biomass components. Concurrently, the photosynthetic capacity peaks in the intermediate stage of ecosystem development. If confirmed by further studies on a larger scale, this study would provide evidence for a predominant effect of plant physiology over soil physical characteristics in the green-up phase after glacier retreat, which has to be taken into account in the creation of scenarios related to climate change and future land use.

Timescale effects on the environmental control of carbon and water fluxes of an apple orchard.

Model parameterization and validation of earth-atmosphere interactions are generally performed using a single timescale (e.g., nearly instantaneous, daily, and annual), although both delayed responses and hysteretic effects have been widely recognized. The lack of consideration of these effects hampers our capability to represent them in empirical- or process-based models. Here we explore, using an apple orchard ecosystem in the North of Italy as a simplified case study, how the considered timescale impacts the relative importance of the single environmental variables in explaining observed net ecosystem exchange (NEE) and evapotranspiration (ET). Using 6 years of eddy covariance and meteorological information as input data, we found a decay of the relative importance of the modeling capability of photosynthetically active radiation in explaining both NEE and ET moving from half-hourly to seasonal timescale and an increase in the relative importance of air temperature (T) and VPD. Satellite NDVI, used as proxy of leaf development, added little improvement to overall modeling capability. Increasing the timescale, the number of variables needed for parameterization decreased (from 5 to 1), while the proportion of variance explained by the model increased (r2 from 0.56-0.78 to 0.85-0.90 for NEE and ET respectively). The wavelet coherence and the phase analyses showed that the two variables that increased their relative importance when the scale increased (T, VPD) were not in phase at the correlation peak of both ET and NEE. This phase shift in the time domain corresponds to a hysteretic response in the meteorological variables domain. This work confirms that the model parameterization should be performed using parameters calculated at the appropriate scale. It suggests that in managed ecosystems, where the interannual variability is minimized by the agronomic practices, the use of timescales large enough to include hysteretic and delayed responses reduces the number of required input variables and improves their explanatory capacity.

Quantifying deforestation and forest degradation with thermal response.

Deforestation and forest degradation cause the deterioration of resources and ecosystem services. However, there are still no operational indicators to measure forest status, especially for forest degradation. In the present study, we analysed the thermal response number (TRN, calculated by daily total net radiation divided by daily temperature range) of 163 sites including mature forest, disturbed forest, planted forest, shrubland, grassland, savanna vegetation and cropland. TRN generally increased with latitude, however the regression of TRN against latitude differed among vegetation types. Mature forests are superior as thermal buffers, and had significantly higher TRN than disturbed and planted forests. There was a clear boundary between TRN of forest and non-forest vegetation (i.e. grassland and savanna) with the exception of shrubland, whose TRN overlapped with that of forest vegetation. We propose to use the TRN of local mature forest as the optimal TRN (TRNopt). A forest with lower than 75% of TRNopt was identified as subjected to significant disturbance, and forests with 66% of TRNopt was the threshold for deforestation within the absolute latitude from 30° to 55°. Our results emphasized the irreplaceable thermal buffer capacity of mature forest. TRN can be used for early warning of deforestation and degradation risk. It is therefore a valuable tool in the effort to protect forests and prevent deforestation.

Colonization of a Deglaciated Moraine: Contrasting Patterns of Carbon Uptake and Release from C3 and CAM Plants.

INTRODUCTION: Current glacier retreat makes vast mountain ranges available for vegetation establishment and growth. As a result, carbon (C) is accumulated in the soil, in a negative feedback to climate change. Little is known about the effective C budget of these new ecosystems and how the presence of different vegetation communities influences CO2 fluxes. METHODS: On the Matsch glacier forefield (Alps, Italy) we measured over two growing seasons the Net Ecosystem Exchange (NEE) of a typical grassland, dominated by the C3 Festuca halleri All., and a community dominated by the CAM rosettes Sempervivum montanum L. Using transparent and opaque chambers, with air temperature as the driver, we partitioned NEE to calculate Ecosystem Respiration (Reco) and Gross Ecosystem Exchange (GEE). In addition, soil and vegetation samples were collected from the same sites to estimate the Net Ecosystem Carbon Balance (NECB). RESULTS: The two communities showed contrasting GEE but similar Reco patterns, and as a result they were significantly different in NEE during the period measured. The grassland acted as a C sink, with a total cumulated value of -46.4±35.5 g C m-2 NEE, while the plots dominated by the CAM rosettes acted as a source, with 31.9±22.4 g C m-2. In spite of the different NEE, soil analysis did not reveal significant differences in carbon accumulation of the two plant communities (1770±130 for F. halleri and 2080±230 g C m-2 for S. montanum), suggesting that processes often neglected, like lateral flows and winter respiration, can have a similar relevance as NEE in the determination of the Net Ecosystem Carbon Balance.

Potential and limitations of inferring ecosystem photosynthetic capacity from leaf functional traits.

The aim of this study was to systematically analyze the potential and limitations of using plant functional trait observations from global databases versus in situ data to improve our understanding of vegetation impacts on ecosystem functional properties (EFPs). Using ecosystem photosynthetic capacity as an example, we first provide an objective approach to derive robust EFP estimates from gross primary productivity (GPP) obtained from eddy covariance flux measurements. Second, we investigate the impact of synchronizing EFPs and plant functional traits in time and space to evaluate their relationships, and the extent to which we can benefit from global plant trait databases to explain the variability of ecosystem photosynthetic capacity. Finally, we identify a set of plant functional traits controlling ecosystem photosynthetic capacity at selected sites. Suitable estimates of the ecosystem photosynthetic capacity can be derived from light response curve of GPP responding to radiation (photosynthetically active radiation or absorbed photosynthetically active radiation). Although the effect of climate is minimized in these calculations, the estimates indicate substantial interannual variation of the photosynthetic capacity, even after removing site-years with confounding factors like disturbance such as fire events. The relationships between foliar nitrogen concentration and ecosystem photosynthetic capacity are tighter when both of the measurements are synchronized in space and time. When using multiple plant traits simultaneously as predictors for ecosystem photosynthetic capacity variation, the combination of leaf carbon to nitrogen ratio with leaf phosphorus content explains the variance of ecosystem photosynthetic capacity best (adjusted R2 = 0.55). Overall, this study provides an objective approach to identify links between leaf level traits and canopy level processes and highlights the relevance of the dynamic nature of ecosystems. Synchronizing measurements of eddy covariance fluxes and plant traits in time and space is shown to be highly relevant to better understand the importance of intra- and interspecific trait variation on ecosystem functioning.

Joint control of terrestrial gross primary productivity by plant phenology and physiology.

Terrestrial gross primary productivity (GPP) varies greatly over time and space. A better understanding of this variability is necessary for more accurate predictions of the future climate-carbon cycle feedback. Recent studies have suggested that variability in GPP is driven by a broad range of biotic and abiotic factors operating mainly through changes in vegetation phenology and physiological processes. However, it is still unclear how plant phenology and physiology can be integrated to explain the spatiotemporal variability of terrestrial GPP. Based on analyses of eddy-covariance and satellite-derived data, we decomposed annual terrestrial GPP into the length of the CO2 uptake period (CUP) and the seasonal maximal capacity of CO2 uptake (GPPmax). The product of CUP and GPPmax explained >90% of the temporal GPP variability in most areas of North America during 2000-2010 and the spatial GPP variation among globally distributed eddy flux tower sites. It also explained GPP response to the European heatwave in 2003 (r(2) = 0.90) and GPP recovery after a fire disturbance in South Dakota (r(2) = 0.88). Additional analysis of the eddy-covariance flux data shows that the interbiome variation in annual GPP is better explained by that in GPPmax than CUP. These findings indicate that terrestrial GPP is jointly controlled by ecosystem-level plant phenology and photosynthetic capacity, and greater understanding of GPPmax and CUP responses to environmental and biological variations will, thus, improve predictions of GPP over time and space.

Terrestrial gross primary production inferred from satellite fluorescence and vegetation models.

Determining the spatial and temporal distribution of terrestrial gross primary production (GPP) is a critical step in closing the Earth's carbon budget. Dynamical global vegetation models (DGVMs) provide mechanistic insight into GPP variability but diverge in predicting the response to climate in poorly investigated regions. Recent advances in the remote sensing of solar-induced chlorophyll fluorescence (SIF) opens up a new possibility to provide direct global observational constraints for GPP. Here, we apply an optimal estimation approach to infer the global distribution of GPP from an ensemble of eight DGVMs constrained by global measurements of SIF from the Greenhouse Gases Observing SATellite (GOSAT). These estimates are compared to flux tower data in N. America, Europe, and tropical S. America, with careful consideration of scale differences between models, GOSAT, and flux towers. Assimilation of GOSAT SIF with DGVMs causes a redistribution of global productivity from northern latitudes to the tropics of 7-8 Pg C yr(-1) from 2010 to 2012, with reduced GPP in northern forests (~3.6 Pg C yr(-1) ) and enhanced GPP in tropical forests (~3.7 Pg C yr(-1) ). This leads to improvements in the structure of the seasonal cycle, including earlier dry season GPP loss and enhanced peak-to-trough GPP in tropical forests within the Amazon Basin and reduced growing season length in northern croplands and deciduous forests. Uncertainty in predicted GPP (estimated from the spread of DGVMs) is reduced by 40-70% during peak productivity suggesting the assimilation of GOSAT SIF with models is well-suited for benchmarking. We conclude that satellite fluorescence augurs a new opportunity to quantify the GPP response to climate drivers and the potential to constrain predictions of carbon cycle evolution.

Thermal optimality of net ecosystem exchange of carbon dioxide and underlying mechanisms.

• It is well established that individual organisms can acclimate and adapt to temperature to optimize their functioning. However, thermal optimization of ecosystems, as an assemblage of organisms, has not been examined at broad spatial and temporal scales. • Here, we compiled data from 169 globally distributed sites of eddy covariance and quantified the temperature response functions of net ecosystem exchange (NEE), an ecosystem-level property, to determine whether NEE shows thermal optimality and to explore the underlying mechanisms. • We found that the temperature response of NEE followed a peak curve, with the optimum temperature (corresponding to the maximum magnitude of NEE) being positively correlated with annual mean temperature over years and across sites. Shifts of the optimum temperature of NEE were mostly a result of temperature acclimation of gross primary productivity (upward shift of optimum temperature) rather than changes in the temperature sensitivity of ecosystem respiration. • Ecosystem-level thermal optimality is a newly revealed ecosystem property, presumably reflecting associated evolutionary adaptation of organisms within ecosystems, and has the potential to significantly regulate ecosystem-climate change feedbacks. The thermal optimality of NEE has implications for understanding fundamental properties of ecosystems in changing environments and benchmarking global models.

Does canopy mean nitrogen concentration explain variation in canopy light use efficiency across 14 contrasting forest sites?

The maximum light use efficiency (LUE = gross primary production (GPP)/absorbed photosynthetic photon flux density (aPPFD)) of plant canopies has been reported to vary spatially and some of this variation has previously been attributed to plant species differences. The canopy nitrogen concentration [N] can potentially explain some of this spatial variation. However, the current paradigm of the N-effect on photosynthesis is largely based on the relationship between photosynthetic capacity (A(max)) and [N], i.e., the effects of [N] on photosynthesis rates appear under high PPFD. A maximum LUE-[N] relationship, if it existed, would influence photosynthesis in the whole range of PPFD. We estimated maximum LUE for 14 eddy-covariance forest sites, examined its [N] dependency and investigated how the [N]-maximum LUE dependency could be incorporated into a GPP model. In the model, maximum LUE corresponds to LUE under optimal environmental conditions before light saturation takes place (the slope of GPP vs. PPFD under low PPFD). Maximum LUE was higher in deciduous/mixed than in coniferous sites, and correlated significantly with canopy mean [N]. Correlations between maximum LUE and canopy [N] existed regardless of daily PPFD, although we expected the correlation to disappear under low PPFD when LUE was also highest. Despite these correlations, including [N] in the model of GPP only marginally decreased the root mean squared error. Our results suggest that maximum LUE correlates linearly with canopy [N], but that a larger body of data is required before we can include this relationship into a GPP model. Gross primary production will therefore positively correlate with [N] already at low PPFD, and not only at high PPFD as is suggested by the prevailing paradigm of leaf-level A(max)-[N] relationships. This finding has consequences for modelling GPP driven by temporal changes or spatial variation in canopy [N].

Recent decline in the global land evapotranspiration trend due to limited moisture supply.

More than half of the solar energy absorbed by land surfaces is currently used to evaporate water. Climate change is expected to intensify the hydrological cycle and to alter evapotranspiration, with implications for ecosystem services and feedback to regional and global climate. Evapotranspiration changes may already be under way, but direct observational constraints are lacking at the global scale. Until such evidence is available, changes in the water cycle on land−a key diagnostic criterion of the effects of climate change and variability−remain uncertain. Here we provide a data-driven estimate of global land evapotranspiration from 1982 to 2008, compiled using a global monitoring network, meteorological and remote-sensing observations, and a machine-learning algorithm. In addition, we have assessed evapotranspiration variations over the same time period using an ensemble of process-based land-surface models. Our results suggest that global annual evapotranspiration increased on average by 7.1 ± 1.0 millimetres per year per decade from 1982 to 1997. After that, coincident with the last major El Niño event in 1998, the global evapotranspiration increase seems to have ceased until 2008. This change was driven primarily by moisture limitation in the Southern Hemisphere, particularly Africa and Australia. In these regions, microwave satellite observations indicate that soil moisture decreased from 1998 to 2008. Hence, increasing soil-moisture limitations on evapotranspiration largely explain the recent decline of the global land-evapotranspiration trend. Whether the changing behaviour of evapotranspiration is representative of natural climate variability or reflects a more permanent reorganization of the land water cycle is a key question for earth system science.

Influence of spring and autumn phenological transitions on forest ecosystem productivity.

We use eddy covariance measurements of net ecosystem productivity (NEP) from 21 FLUXNET sites (153 site-years of data) to investigate relationships between phenology and productivity (in terms of both NEP and gross ecosystem photosynthesis, GEP) in temperate and boreal forests. Results are used to evaluate the plausibility of four different conceptual models. Phenological indicators were derived from the eddy covariance time series, and from remote sensing and models. We examine spatial patterns (across sites) and temporal patterns (across years); an important conclusion is that it is likely that neither of these accurately represents how productivity will respond to future phenological shifts resulting from ongoing climate change. In spring and autumn, increased GEP resulting from an 'extra' day tends to be offset by concurrent, but smaller, increases in ecosystem respiration, and thus the effect on NEP is still positive. Spring productivity anomalies appear to have carry-over effects that translate to productivity anomalies in the following autumn, but it is not clear that these result directly from phenological anomalies. Finally, the productivity of evergreen needleleaf forests is less sensitive to phenology than is productivity of deciduous broadleaf forests. This has implications for how climate change may drive shifts in competition within mixed-species stands.

Global convergence in the temperature sensitivity of respiration at ecosystem level.

The respiratory release of carbon dioxide (CO(2)) from the land surface is a major flux in the global carbon cycle, antipodal to photosynthetic CO(2) uptake. Understanding the sensitivity of respiratory processes to temperature is central for quantifying the climate-carbon cycle feedback. We approximated the sensitivity of terrestrial ecosystem respiration to air temperature (Q(10)) across 60 FLUXNET sites with the use of a methodology that circumvents confounding effects. Contrary to previous findings, our results suggest that Q(10) is independent of mean annual temperature, does not differ among biomes, and is confined to values around 1.4 +/- 0.1. The strong relation between photosynthesis and respiration, by contrast, is highly variable among sites. The results may partly explain a less pronounced climate-carbon cycle feedback than suggested by current carbon cycle climate models.