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.sup.-2 yr.sup.-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.sup.-2 yr.sup.-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 (SR.sub.ref) 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 (SR.sub.ref). 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 CO.sub.2 . 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 SR.sub.ref and moderate Q10, may be more resilient.
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.
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 CO ₂ uptake period (CUP) and the seasonal maximal capacity of CO ₂ uptake (GPP ₘₐₓ). The product of CUP and GPP ₘₐₓ 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 ² = 0.90) and GPP recovery after a fire disturbance in South Dakota ( r ² = 0.88). Additional analysis of the eddy–covariance flux data shows that the interbiome variation in annual GPP is better explained by that in GPP ₘₐₓ than CUP. These findings indicate that terrestrial GPP is jointly controlled by ecosystem-level plant phenology and photosynthetic capacity, and greater understanding of GPP ₘₐₓ and CUP responses to environmental and biological variations will, thus, improve predictions of GPP over time and space.
Significance Terrestrial gross primary productivity (GPP), the total photosynthetic CO ₂ fixation at ecosystem level, fuels all life on land. However, its spatiotemporal variability is poorly understood, because GPP is determined by many processes related to plant phenology and physiological activities. In this study, we find that plant phenological and physiological properties can be integrated in a robust index—the product of the length of CO ₂ uptake period and the seasonal maximal photosynthesis—to explain the GPP variability over space and time in response to climate extremes and during recovery after disturbance.
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.
•Seven light use efficiency models were compared at global eddy covariance towers.•Performance of seven models differed substantially among ecosystem types.•It is needed to improve LUE models by ...integrating impacts of diffuse radiation and reliable water stress equations.
Simulating gross primary productivity (GPP) of terrestrial ecosystems has been a major challenge in quantifying the global carbon cycle. Many different light use efficiency (LUE) models have been developed recently, but our understanding of the relative merits of different models remains limited. Using CO2 flux measurements from multiple eddy covariance sites, we here compared and assessed major algorithms and performance of seven LUE models (CASA, CFix, CFlux, EC-LUE, MODIS, VPM and VPRM). Comparison between simulated GPP and estimated GPP from flux measurements showed that model performance differed substantially among ecosystem types. In general, most models performed better in capturing the temporal changes and magnitude of GPP in deciduous broadleaf forests and mixed forests than in evergreen broadleaf forests and shrublands. Six of the seven LUE models significantly underestimated GPP during cloudy days because the impacts of diffuse radiation on light use efficiency were ignored in the models. CFlux and EC-LUE exhibited the lowest root mean square error among all models at 80% and 75% of the sites, respectively. Moreover, these two models showed better performance than others in simulating interannual variability of GPP. Two pairwise comparisons revealed that the seven models differed substantially in algorithms describing the environmental regulations, particularly water stress, on GPP. This analysis highlights the need to improve representation of the impacts of diffuse radiation and water stress in the LUE models.
We upscaled FLUXNET observations of carbon dioxide, water, and energy fluxes to the global scale using the machine learning technique, model tree ensembles (MTE). We trained MTE to predict site‐level ...gross primary productivity (GPP), terrestrial ecosystem respiration (TER), net ecosystem exchange (NEE), latent energy (LE), and sensible heat (H) based on remote sensing indices, climate and meteorological data, and information on land use. We applied the trained MTEs to generate global flux fields at a 0.5° × 0.5° spatial resolution and a monthly temporal resolution from 1982 to 2008. Cross‐validation analyses revealed good performance of MTE in predicting among‐site flux variability with modeling efficiencies (MEf) between 0.64 and 0.84, except for NEE (MEf = 0.32). Performance was also good for predicting seasonal patterns (MEf between 0.84 and 0.89, except for NEE (0.64)). By comparison, predictions of monthly anomalies were not as strong (MEf between 0.29 and 0.52). Improved accounting of disturbance and lagged environmental effects, along with improved characterization of errors in the training data set, would contribute most to further reducing uncertainties. Our global estimates of LE (158 ± 7 J × 1018 yr−1), H (164 ± 15 J × 1018 yr−1), and GPP (119 ± 6 Pg C yr−1) were similar to independent estimates. Our global TER estimate (96 ± 6 Pg C yr−1) was likely underestimated by 5–10%. Hot spot regions of interannual variability in carbon fluxes occurred in semiarid to semihumid regions and were controlled by moisture supply. Overall, GPP was more important to interannual variability in NEE than TER. Our empirically derived fluxes may be used for calibration and evaluation of land surface process models and for exploratory and diagnostic assessments of the biosphere.
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⁻¹ from 2010 to 2012, with reduced GPP in northern forests (~3.6 Pg C yr⁻¹) and enhanced GPP in tropical forests (~3.7 Pg C yr⁻¹). 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.
► Mean energy balance closure at 173 FLUXNET sites is 0.84. ► Mean forest and non-forest closure does not differ. ► Significant differences in closure were found among plant functional types. ► ...Landscape-level vegetation variability should not be excluded from the interpretation.
The energy balance at most surface-atmosphere flux research sites remains unclosed. The mechanisms underlying the discrepancy between measured energy inputs and outputs across the global FLUXNET tower network are still under debate. Recent reviews have identified exchange processes and turbulent motions at large spatial and temporal scales in heterogeneous landscapes as the primary cause of the lack of energy balance closure at some intensively-researched sites, while unmeasured storage terms cannot be ruled out as a dominant contributor to the lack of energy balance closure at many other sites. We analyzed energy balance closure across 173 ecosystems in the FLUXNET database and explored the relationship between energy balance closure and landscape heterogeneity using MODIS products and GLOBEstat elevation data. Energy balance closure per research site (CEB,s) averaged 0.84±0.20, with best average closures in evergreen broadleaf forests and savannas (0.91–0.94) and worst average closures in crops, deciduous broadleaf forests, mixed forests and wetlands (0.70–0.78). Half-hourly or hourly energy balance closure on a percent basis increased with friction velocity (u*) and was highest on average under near-neutral atmospheric conditions. CEB,s was significantly related to mean precipitation, gross primary productivity and landscape-level enhanced vegetation index (EVI) from MODIS, and the variability in elevation, MODIS plant functional type, and MODIS EVI. A linear model including landscape-level variability in both EVI and elevation, mean precipitation, and an interaction term between EVI variability and precipitation had the lowest Akaike's information criterion value. CEB,s in landscapes with uniform plant functional type approached 0.9 and CEB,s in landscapes with uniform EVI approached 1. These results suggest that landscape-level heterogeneity in vegetation and topography cannot be ignored as a contributor to incomplete energy balance closure at the flux network level, although net radiation measurements, biological energy assimilation, unmeasured storage terms, and the importance of good practice including site selection when making flux measurements should not be discounted. Our results suggest that future research should focus on the quantitative mechanistic relationships between energy balance closure and landscape-scale heterogeneity, and the consequences of mesoscale circulations for surface-atmosphere exchange measurements.
The respiratory release of carbon dioxide (CO₂) from the land surface is a major flux in the global carbon cycle, antipodal to photosynthetic CO₂ 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₁₀) across 60 FLUXNET sites with the use of a methodology that circumvents confounding effects. Contrary to previous findings, our results suggest that Q₁₀ 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.