Diagnosing destabilization risk in global land carbon sinks.

Global net land carbon uptake or net biome production (NBP) has increased during recent decades1. Whether its temporal variability and autocorrelation have changed during this period, however, remains elusive, even though an increase in both could indicate an increased potential for a destabilized carbon sink2,3. Here, we investigate the trends and controls of net terrestrial carbon uptake and its temporal variability and autocorrelation from 1981 to 2018 using two atmospheric-inversion models, the amplitude of the seasonal cycle of atmospheric CO2 concentration derived from nine monitoring stations distributed across the Pacific Ocean and dynamic global vegetation models. We find that annual NBP and its interdecadal variability increased globally whereas temporal autocorrelation decreased. We observe a separation of regions characterized by increasingly variable NBP, associated with warm regions and increasingly variable temperatures, lower and weaker positive trends in NBP and regions where NBP became stronger and less variable. Plant species richness presented a concave-down parabolic spatial relationship with NBP and its variability at the global scale whereas nitrogen deposition generally increased NBP. Increasing temperature and its increasing variability appear as the most important drivers of declining and increasingly variable NBP. Our results show increasing variability of NBP regionally that can be mostly attributed to climate change and that may point to destabilization of the coupled carbon-climate system.

Widespread seasonal compensation effects of spring warming on northern plant productivity.

Climate change is shifting the phenological cycles of plants1, thereby altering the functioning of ecosystems, which in turn induces feedbacks to the climate system2. In northern (north of 30° N) ecosystems, warmer springs lead generally to an earlier onset of the growing season3,4 and increased ecosystem productivity early in the season5. In situ6 and regional7-9 studies also provide evidence for lagged effects of spring warmth on plant productivity during the subsequent summer and autumn. However, our current understanding of these lagged effects, including their direction (beneficial or adverse) and geographic distribution, is still very limited. Here we analyse satellite, field-based and modelled data for the period 1982-2011 and show that there are widespread and contrasting lagged productivity responses to spring warmth across northern ecosystems. On the basis of the observational data, we find that roughly 15 per cent of the total study area of about 41 million square kilometres exhibits adverse lagged effects and that roughly 5 per cent of the total study area exhibits beneficial lagged effects. By contrast, current-generation terrestrial carbon-cycle models predict much lower areal fractions of adverse lagged effects (ranging from 1 to 14 per cent) and much higher areal fractions of beneficial lagged effects (ranging from 9 to 54 per cent). We find that elevation and seasonal precipitation patterns largely dictate the geographic pattern and direction of the lagged effects. Inadequate consideration in current models of the effects of the seasonal build-up of water stress on seasonal vegetation growth may therefore be able to explain the differences that we found between our observation-constrained estimates and the model-constrained estimates of lagged effects associated with spring warming. Overall, our results suggest that for many northern ecosystems the benefits of warmer springs on growing-season ecosystem productivity are effectively compensated for by the accumulation of seasonal water deficits, despite the fact that northern ecosystems are thought to be largely temperature- and radiation-limited10.

Compensatory water effects link yearly global land CO2 sink changes to temperature.

Large interannual variations in the measured growth rate of atmospheric carbon dioxide (CO2) originate primarily from fluctuations in carbon uptake by land ecosystems. It remains uncertain, however, to what extent temperature and water availability control the carbon balance of land ecosystems across spatial and temporal scales. Here we use empirical models based on eddy covariance data and process-based models to investigate the effect of changes in temperature and water availability on gross primary productivity (GPP), terrestrial ecosystem respiration (TER) and net ecosystem exchange (NEE) at local and global scales. We find that water availability is the dominant driver of the local interannual variability in GPP and TER. To a lesser extent this is true also for NEE at the local scale, but when integrated globally, temporal NEE variability is mostly driven by temperature fluctuations. We suggest that this apparent paradox can be explained by two compensatory water effects. Temporal water-driven GPP and TER variations compensate locally, dampening water-driven NEE variability. Spatial water availability anomalies also compensate, leaving a dominant temperature signal in the year-to-year fluctuations of the land carbon sink. These findings help to reconcile seemingly contradictory reports regarding the importance of temperature and water in controlling the interannual variability of the terrestrial carbon balance. Our study indicates that spatial climate covariation drives the global carbon cycle response.

Greening drylands despite warming consistent with carbon dioxide fertilization effect.

The rising atmospheric CO2 concentration leads to a CO2 fertilization effect on plants-that is, increased photosynthetic uptake of CO2 by leaves and enhanced water-use efficiency (WUE). Yet, the resulting net impact of CO2 fertilization on plant growth and soil moisture (SM) savings at large scale is poorly understood. Drylands provide a natural experimental setting to detect the CO2 fertilization effect on plant growth since foliage amount, plant water-use and photosynthesis are all tightly coupled in water-limited ecosystems. A long-term change in the response of leaf area index (LAI, a measure of foliage amount) to changes in SM is likely to stem from changing water demand of primary productivity in water-limited ecosystems and is a proxy for changes in WUE. Using 34-year satellite observations of LAI and SM over tropical and subtropical drylands, we identify that a 1% increment in SM leads to 0.15% (±0.008, 95% confidence interval) and 0.51% (±0.01, 95% confidence interval) increments in LAI during 1982-1998 and 1999-2015, respectively. The increasing response of LAI to SM has contributed 7.2% (±3.0%, 95% confidence interval) to total dryland greening during 1999-2015 compared to 1982-1998. The increasing response of LAI to SM is consistent with the CO2 fertilization effect on WUE in water-limited ecosystems, indicating that a given amount of SM has sustained greater amounts of photosynthetic foliage over time. The LAI responses to changes in SM from seven dynamic global vegetation models are not always consistent with observations, highlighting the need for improved process knowledge of terrestrial ecosystem responses to rising atmospheric CO2 concentration.

Causes of slowing-down seasonal CO2 amplitude at Mauna Loa.

Changing amplitude of the seasonal cycle of atmospheric CO2 (SCA) in the northern hemisphere is an emerging carbon cycle property. Mauna Loa (MLO) station (20°N, 156°W), which has the longest continuous northern hemisphere CO2 record, shows an increasing SCA before the 1980s (p < .01), followed by no significant change thereafter. We analyzed the potential driving factors of SCA slowing-down, with an ensemble of dynamic global vegetation models (DGVMs) coupled with an atmospheric transport model. We found that slowing-down of SCA at MLO is primarily explained by response of net biome productivity (NBP) to climate change, and by changes in atmospheric circulations. Through NBP, climate change increases SCA at MLO before the 1980s and decreases it afterwards. The effect of climate change on the slowing-down of SCA at MLO is mainly exerted by intensified drought stress acting to offset the acceleration driven by CO2 fertilization. This challenges the view that CO2 fertilization is the dominant cause of emergent SCA trends at northern sites south of 40°N. The contribution of agricultural intensification on the deceleration of SCA at MLO was elusive according to land-atmosphere CO2 flux estimated by DGVMs and atmospheric inversions. Our results also show the necessity to adequately account for changing circulation patterns in understanding carbon cycle dynamics observed from atmospheric observations and in using these observations to benchmark DGVMs.

State of the science in reconciling top-down and bottom-up approaches for terrestrial CO2 budget.

Robust estimates of CO2 budget, CO2 exchanged between the atmosphere and terrestrial biosphere, are necessary to better understand the role of the terrestrial biosphere in mitigating anthropogenic CO2 emissions. Over the past decade, this field of research has advanced through understanding of the differences and similarities of two fundamentally different approaches: "top-down" atmospheric inversions and "bottom-up" biosphere models. Since the first studies were undertaken, these approaches have shown an increasing level of agreement, but disagreements in some regions still persist, in part because they do not estimate the same quantity of atmosphere-biosphere CO2 exchange. Here, we conducted a thorough comparison of CO2 budgets at multiple scales and from multiple methods to assess the current state of the science in estimating CO2 budgets. Our set of atmospheric inversions and biosphere models, which were adjusted for a consistent flux definition, showed a high level of agreement for global and hemispheric CO2 budgets in the 2000s. Regionally, improved agreement in CO2 budgets was notable for North America and Southeast Asia. However, large gaps between the two methods remained in East Asia and South America. In other regions, Europe, boreal Asia, Africa, South Asia, and Oceania, it was difficult to determine whether those regions act as a net sink or source because of the large spread in estimates from atmospheric inversions. These results highlight two research directions to improve the robustness of CO2 budgets: (a) to increase representation of processes in biosphere models that could contribute to fill the budget gaps, such as forest regrowth and forest degradation; and (b) to reduce sink-source compensation between regions (dipoles) in atmospheric inversion so that their estimates become more comparable. Advancements on both research areas will increase the level of agreement between the top-down and bottom-up approaches and yield more robust knowledge of regional CO2 budgets.

Global vegetation biomass production efficiency constrained by models and observations.

Plants use only a fraction of their photosynthetically derived carbon for biomass production (BP). The biomass production efficiency (BPE), defined as the ratio of BP to photosynthesis, and its variation across and within vegetation types is poorly understood, which hinders our capacity to accurately estimate carbon turnover times and carbon sinks. Here, we present a new global estimation of BPE obtained by combining field measurements from 113 sites with 14 carbon cycle models. Our best estimate of global BPE is 0.41 ± 0.05, excluding cropland. The largest BPE is found in boreal forests (0.48 ± 0.06) and the lowest in tropical forests (0.40 ± 0.04). Carbon cycle models overestimate BPE, although models with carbon-nitrogen interactions tend to be more realistic. Using observation-based estimates of global photosynthesis, we quantify the global BP of non-cropland ecosystems of 41 ± 6 Pg C/year. This flux is less than net primary production as it does not contain carbon allocated to symbionts, used for exudates or volatile carbon compound emissions to the atmosphere. Our study reveals a positive bias of 24 ± 11% in the model-estimated BP (10 of 14 models). When correcting models for this bias while leaving modeled carbon turnover times unchanged, we found that the global ecosystem carbon storage change during the last century is decreased by 67% (or 58 Pg C).

Climate-Driven Variability and Trends in Plant Productivity Over Recent Decades Based on Three Global Products.

Variability in climate exerts a strong influence on vegetation productivity (gross primary productivity; GPP), and therefore has a large impact on the land carbon sink. However, no direct observations of global GPP exist, and estimates rely on models that are constrained by observations at various spatial and temporal scales. Here, we assess the consistency in GPP from global products which extend for more than three decades; two observation-based approaches, the upscaling of FLUXNET site observations (FLUXCOM) and a remote sensing derived light use efficiency model (RS-LUE), and from a suite of terrestrial biosphere models (TRENDYv6). At local scales, we find high correlations in annual GPP among the products, with exceptions in tropical and high northern latitudes. On longer time scales, the products agree on the direction of trends over 58% of the land, with large increases across northern latitudes driven by warming trends. Further, tropical regions exhibit the largest interannual variability in GPP, with both rainforests and savannas contributing substantially. Variability in savanna GPP is likely predominantly driven by water availability, although temperature could play a role via soil moisture-atmosphere feedbacks. There is, however, no consensus on the magnitude and driver of variability of tropical forests, which suggest uncertainties in process representations and underlying observations remain. These results emphasize the need for more direct long-term observations of GPP along with an extension of in situ networks in underrepresented regions (e.g., tropical forests). Such capabilities would support efforts to better validate relevant processes in models, to more accurately estimate GPP.

Quantifying the biophysical and socioeconomic drivers of changes in forest and agricultural land in South and Southeast Asia.

South and Southeast Asia (SSEA) has been a hotspot for land use and land cover change (LULCC) in the past few decades. The identification and quantification of the drivers of LULCC are crucial for improving our understanding of LULCC trends. So far, the biophysical and socioeconomic drivers of forest change have not been quantified at the regional scale, particularly for SSEA. In this study, we quantify the biophysical and socioeconomic drivers of forest change on a country-by-country basis in SSEA using an integrated quantitative methodology, which systematically accounts for previously published driver information and regional datasets. We synthesize more than 200 publications to identify the drivers of the forest change at different spatial scales in SSEA. Subsequently, we collect spatially explicit proxy data to represent the identified drivers. We quantify the dynamics of forest and agricultural land from 1992 to 2015 using the Climate Change Initiative (CCI) land cover data developed by the European Space Agency (ESA). A geographically weighted regression method is employed to quantify the spatially heterogeneous drivers of forest change. Our results show that socioeconomic drivers are more important than biophysical drivers for the conversion of forest to agricultural land in South Asia and maritime Southeast Asia. In contrast, biophysical drivers are more important than socioeconomic drivers for the conversion of agricultural land to forest in maritime Southeast Asia and less important in South Asia. Both biophysical and socioeconomic drivers contribute approximately equally to both changes in the mainland Southeast Asia region. By quantifying the dynamics of forest and agricultural land and the spatially explicit drivers of their changes in SSEA, this study provides a solid foundation for LULCC modeling and projection.

Challenging terrestrial biosphere models with data from the long-term multifactor Prairie Heating and CO2 Enrichment experiment.

Multifactor experiments are often advocated as important for advancing terrestrial biosphere models (TBMs), yet to date, such models have only been tested against single-factor experiments. We applied 10 TBMs to the multifactor Prairie Heating and CO2 Enrichment (PHACE) experiment in Wyoming, USA. Our goals were to investigate how multifactor experiments can be used to constrain models and to identify a road map for model improvement. We found models performed poorly in ambient conditions; there was a wide spread in simulated above-ground net primary productivity (range: 31-390 g C m-2  yr-1 ). Comparison with data highlighted model failures particularly with respect to carbon allocation, phenology, and the impact of water stress on phenology. Performance against the observations from single-factors treatments was also relatively poor. In addition, similar responses were predicted for different reasons across models: there were large differences among models in sensitivity to water stress and, among the N cycle models, N availability during the experiment. Models were also unable to capture observed treatment effects on phenology: they overestimated the effect of warming on leaf onset and did not allow CO2 -induced water savings to extend the growing season length. Observed interactive (CO2  × warming) treatment effects were subtle and contingent on water stress, phenology, and species composition. As the models did not correctly represent these processes under ambient and single-factor conditions, little extra information was gained by comparing model predictions against interactive responses. We outline a series of key areas in which this and future experiments could be used to improve model predictions of grassland responses to global change.

Assessing uncertainties in land cover projections.

Understanding uncertainties in land cover projections is critical to investigating land-based climate mitigation policies, assessing the potential of climate adaptation strategies and quantifying the impacts of land cover change on the climate system. Here, we identify and quantify uncertainties in global and European land cover projections over a diverse range of model types and scenarios, extending the analysis beyond the agro-economic models included in previous comparisons. The results from 75 simulations over 18 models are analysed and show a large range in land cover area projections, with the highest variability occurring in future cropland areas. We demonstrate systematic differences in land cover areas associated with the characteristics of the modelling approach, which is at least as great as the differences attributed to the scenario variations. The results lead us to conclude that a higher degree of uncertainty exists in land use projections than currently included in climate or earth system projections. To account for land use uncertainty, it is recommended to use a diverse set of models and approaches when assessing the potential impacts of land cover change on future climate. Additionally, further work is needed to better understand the assumptions driving land use model results and reveal the causes of uncertainty in more depth, to help reduce model uncertainty and improve the projections of land cover.

Dynamics and determinants of land change in India: integrating satellite data with village socioeconomics.

We examine the dynamics and spatial determinants of land change in India by integrating decadal land cover maps (1985-1995-2005) from a wall-to-wall analysis of Landsat images with spatiotemporal socioeconomic database for ~630,000 villages in India. We reinforce our results through collective evidence from synthesis of 102 case studies that incorporate field knowledge of the causes of land change in India. We focus on cropland-fallow land conversions, and forest area changes (excludes non-forest tree categories including commercial plantations). We show that cropland to fallow conversions are prominently associated with lack of irrigation and capital, male agricultural labor shortage, and fragmentation of land holdings. We find gross forest loss is substantial and increased from ~23,810 km2 (1985-1995) to ~25,770 km2 (1995-2005). The gross forest gain also increased from ~6000 km2 (1985-1995) to ~7440 km2 (1995-2005). Overall, India experienced a net decline in forest by ~18,000 km2 (gross loss-gross gain) consistently during both decades. We show that the major source of forest loss was cropland expansion in areas of low cropland productivity (due to soil degradation and lack of irrigation), followed by industrial development and mining/quarrying activities, and excessive economic dependence of villages on forest resources.

Global change pressures on soils from land use and management.

Soils are subject to varying degrees of direct or indirect human disturbance, constituting a major global change driver. Factoring out natural from direct and indirect human influence is not always straightforward, but some human activities have clear impacts. These include land-use change, land management and land degradation (erosion, compaction, sealing and salinization). The intensity of land use also exerts a great impact on soils, and soils are also subject to indirect impacts arising from human activity, such as acid deposition (sulphur and nitrogen) and heavy metal pollution. In this critical review, we report the state-of-the-art understanding of these global change pressures on soils, identify knowledge gaps and research challenges and highlight actions and policies to minimize adverse environmental impacts arising from these global change drivers. Soils are central to considerations of what constitutes sustainable intensification. Therefore, ensuring that vulnerable and high environmental value soils are considered when protecting important habitats and ecosystems, will help to reduce the pressure on land from global change drivers. To ensure that soils are protected as part of wider environmental efforts, a global soil resilience programme should be considered, to monitor, recover or sustain soil fertility and function, and to enhance the ecosystem services provided by soils. Soils cannot, and should not, be considered in isolation of the ecosystems that they underpin and vice versa. The role of soils in supporting ecosystems and natural capital needs greater recognition. The lasting legacy of the International Year of Soils in 2015 should be to put soils at the centre of policy supporting environmental protection and sustainable development.

Hotspots of uncertainty in land-use and land-cover change projections: a global-scale model comparison.

Model-based global projections of future land-use and land-cover (LULC) change are frequently used in environmental assessments to study the impact of LULC change on environmental services and to provide decision support for policy. These projections are characterized by a high uncertainty in terms of quantity and allocation of projected changes, which can severely impact the results of environmental assessments. In this study, we identify hotspots of uncertainty, based on 43 simulations from 11 global-scale LULC change models representing a wide range of assumptions of future biophysical and socioeconomic conditions. We attribute components of uncertainty to input data, model structure, scenario storyline and a residual term, based on a regression analysis and analysis of variance. From this diverse set of models and scenarios, we find that the uncertainty varies, depending on the region and the LULC type under consideration. Hotspots of uncertainty appear mainly at the edges of globally important biomes (e.g., boreal and tropical forests). Our results indicate that an important source of uncertainty in forest and pasture areas originates from different input data applied in the models. Cropland, in contrast, is more consistent among the starting conditions, while variation in the projections gradually increases over time due to diverse scenario assumptions and different modeling approaches. Comparisons at the grid cell level indicate that disagreement is mainly related to LULC type definitions and the individual model allocation schemes. We conclude that improving the quality and consistency of observational data utilized in the modeling process and improving the allocation mechanisms of LULC change models remain important challenges. Current LULC representation in environmental assessments might miss the uncertainty arising from the diversity of LULC change modeling approaches, and many studies ignore the uncertainty in LULC projections in assessments of LULC change impacts on climate, water resources or biodiversity.

The Interplay Between Bioenergy Grass Production and Water Resources in the United States of America.

We apply a land surface model to evaluate the interplay between potential bioenergy grass (Miscanthus, Cave-in-Rock, and Alamo) production, water quantity, and nitrogen leaching (NL) in the Central and Eastern U.S. Water use intensity tends to be lower where grass yields are modeled to be high, for example in the Midwest for Miscanthus and Cave-in-Rock and the upper southeastern U.S. for Alamo. However, most of these regions are already occupied by crops and forests and substitution of these biome types for ethanol production implies trade-offs. In general, growing Miscanthus consumes more water, Alamo consumes less water, and Cave-in-Rock consumes approximately the same amount of water as existing vegetation. Bioenergy grasses can maintain high productivity over time, even in water limited regions, because their roots can grow deeper and extract the water from the deep, moist soil layers. However, this may not hold where there are frequent and intense drought events, particularly in regions with shallow soil depths. One advantage of bioenergy grasses is that they mitigate nitrogen leaching relative to row crops and herbaceous plants when grown without applying N fertilizer; and bioenergy grasses, especially Miscanthus, generally require less N fertilizer application than row crops and herbaceous plants.

Global patterns and controls of soil organic carbon dynamics as simulated by multiple terrestrial biosphere models: Current status and future directions.

Soil is the largest organic carbon (C) pool of terrestrial ecosystems, and C loss from soil accounts for a large proportion of land-atmosphere C exchange. Therefore, a small change in soil organic C (SOC) can affect atmospheric carbon dioxide (CO2) concentration and climate change. In the past decades, a wide variety of studies have been conducted to quantify global SOC stocks and soil C exchange with the atmosphere through site measurements, inventories, and empirical/process-based modeling. However, these estimates are highly uncertain, and identifying major driving forces controlling soil C dynamics remains a key research challenge. This study has compiled century-long (1901-2010) estimates of SOC storage and heterotrophic respiration (Rh) from 10 terrestrial biosphere models (TBMs) in the Multi-scale Synthesis and Terrestrial Model Intercomparison Project and two observation-based data sets. The 10 TBM ensemble shows that global SOC estimate ranges from 425 to 2111 Pg C (1 Pg = 1015 g) with a median value of 1158 Pg C in 2010. The models estimate a broad range of Rh from 35 to 69 Pg C yr-1 with a median value of 51 Pg C yr-1 during 2001-2010. The largest uncertainty in SOC stocks exists in the 40-65°N latitude whereas the largest cross-model divergence in Rh are in the tropics. The modeled SOC change during 1901-2010 ranges from -70 Pg C to 86 Pg C, but in some models the SOC change has a different sign from the change of total C stock, implying very different contribution of vegetation and soil pools in determining the terrestrial C budget among models. The model ensemble-estimated mean residence time of SOC shows a reduction of 3.4 years over the past century, which accelerate C cycling through the land biosphere. All the models agreed that climate and land use changes decreased SOC stocks, while elevated atmospheric CO2 and nitrogen deposition over intact ecosystems increased SOC stocks-even though the responses varied significantly among models. Model representations of temperature and moisture sensitivity, nutrient limitation, and land use partially explain the divergent estimates of global SOC stocks and soil C fluxes in this study. In addition, a major source of systematic error in model estimations relates to nonmodeled SOC storage in wetlands and peatlands, as well as to old C storage in deep soil layers.

Managing Multiple Mandates: A System of Systems Model to Analyze Strategies for Producing Cellulosic Ethanol and Reducing Riverine Nitrate Loads in the Upper Mississippi River Basin.

Implementing public policies often involves navigating an array of choices that have economic and environmental consequences that are difficult to quantify due to the complexity of multiple system interactions. Implementing the mandate for cellulosic biofuel production in the Renewable Fuel Standard (RFS) and reducing hypoxia in the northern Gulf of Mexico by reducing riverine nitrate-N loads represent two such cases that overlap in the Mississippi River Basin. To quantify the consequences of these interactions, a system of systems (SoS) model was developed that incorporates interdependencies among the various subsystems, including biofuel refineries, transportation, agriculture, water resources and crop/ethanol markets. The model allows examination of the impact of imposing riverine nitrate-N load limits on the biofuel production system as a whole, including land use change and infrastructure needs. The synergies of crop choice (first versus second generation biofuel crops), infrastructure development, and environmental impacts (streamflow and nitrate-N load) were analyzed to determine the complementarities and trade-offs between environmental protection and biofuel development objectives. For example, the results show that meeting the cellulosic biofuel target in the RFS using Miscanthus x giganteus reduces system profits by 8% and reduces nitrate-N loads by 12% compared to the scenario without a mandate. However, greater water consumption by Miscanthus is likely to reduce streamflow with potentially adverse environmental consequences that need to be considered in future decision making.

Where does the carbon go? A model-data intercomparison of vegetation carbon allocation and turnover processes at two temperate forest free-air CO2 enrichment sites.

Elevated atmospheric CO2 concentration (eCO2) has the potential to increase vegetation carbon storage if increased net primary production causes increased long-lived biomass. Model predictions of eCO2 effects on vegetation carbon storage depend on how allocation and turnover processes are represented. We used data from two temperate forest free-air CO2 enrichment (FACE) experiments to evaluate representations of allocation and turnover in 11 ecosystem models. Observed eCO2 effects on allocation were dynamic. Allocation schemes based on functional relationships among biomass fractions that vary with resource availability were best able to capture the general features of the observations. Allocation schemes based on constant fractions or resource limitations performed less well, with some models having unintended outcomes. Few models represent turnover processes mechanistically and there was wide variation in predictions of tissue lifespan. Consequently, models did not perform well at predicting eCO2 effects on vegetation carbon storage. Our recommendations to reduce uncertainty include: use of allocation schemes constrained by biomass fractions; careful testing of allocation schemes; and synthesis of allocation and turnover data in terms of model parameters. Data from intensively studied ecosystem manipulation experiments are invaluable for constraining models and we recommend that such experiments should attempt to fully quantify carbon, water and nutrient budgets.

Climate-driven uncertainties in modeling terrestrial energy and water fluxes: a site-level to global-scale analysis.

We used a land surface model constrained using data from flux tower sites, to analyze the biases in ecosystem energy and water fluxes arising due to the use of meteorological reanalysis datasets. Following site-level model calibration encompassing major vegetation types from the tropics to the northern high-latitudes, we repeated the site and global simulations using two reanalysis datasets: the NCEP/NCAR and the CRUNCEP. In comparison with the model simulations using observed meteorology from sites, the reanalysis-driven simulations produced several systematic biases in net radiation (Rn ), latent heat (LE), and sensible heat (H) fluxes. These include: (i) persistently positive tropical/subtropical biases in Rn using the NCEP/NCAR, and gradually transitioning to negative Rn biases in the higher latitudes; (ii) large positive H biases in the tropics/subtropics using the NCEP/NCAR; (iii) negative LE biases using the NCEP/NCAR above 40°N; (iv) high tropical LE using the CRUNCEP in comparison with observationally derived global estimates; and (v) flux-partitioning biases from canopy and ground components. Across vegetation types, we investigated the role of the meteorological drivers (shortwave and longwave radiation, atmospheric humidity, temperature, precipitation) and their seasonal biases in controlling these reanalysis-driven uncertainties. At the global scale, our site-level analysis explains several model-data differences in the LE and H fluxes when compared with observationally derived global estimates of these fluxes. Using our results, we discuss the implications of site-level model calibration on subsequent regional/global applications to study energy and hydrological processes. The flux-partitioning biases presented in this study have potential implications on the couplings among terrestrial carbon, energy, and water fluxes, and for the calibration of land-atmosphere parameterizations that are dependent on LE/H partitioning.

Climate-driven uncertainties in modeling terrestrial gross primary production: a site level to global-scale analysis.

We used a land surface model to quantify the causes and extents of biases in terrestrial gross primary production (GPP) due to the use of meteorological reanalysis datasets. We first calibrated the model using meteorology and eddy covariance data from 25 flux tower sites ranging from the tropics to the northern high latitudes and subsequently repeated the site simulations using two reanalysis datasets: NCEP/NCAR and CRUNCEP. The results show that at most sites, the reanalysis-driven GPP bias was significantly positive with respect to the observed meteorology-driven simulations. Notably, the absolute GPP bias was highest at the tropical evergreen tree sites, averaging up to ca. 0.45 kg C m(-2)  yr(-1) across sites (ca. 15% of site level GPP). At the northern mid-/high-latitude broadleaf deciduous and the needleleaf evergreen tree sites, the corresponding annual GPP biases were up to 20%. For the nontree sites, average annual biases of up to ca. 20-30% were simulated within savanna, grassland, and shrubland vegetation types. At the tree sites, the biases in short-wave radiation and humidity strongly influenced the GPP biases, while the nontree sites were more affected by biases in factors controlling water stress (precipitation, humidity, and air temperature). In this study, we also discuss the influence of seasonal patterns of meteorological biases on GPP. Finally, using model simulations for the global land surface, we discuss the potential impacts of site-level reanalysis-driven biases on the global estimates of GPP. In a broader context, our results can have important consequences on other terrestrial ecosystem fluxes (e.g., net primary production, net ecosystem production, energy/water fluxes) and reservoirs (e.g., soil carbon stocks). In a complementary study (Barman et al., ), we extend the present analysis for latent and sensible heat fluxes, thus consistently integrating the analysis of climate-driven uncertainties in carbon, energy, and water fluxes using a single modeling framework.