Forest demography and biomass accumulation rates are associated with transient mean tree size vs. density scaling relations.

Linking individual and stand-level dynamics during forest development reveals a scaling relationship between mean tree size and tree density in forest stands, which integrates forest structure and function. However, the nature of this so-called scaling law and its variation across broad spatial scales remain unquantified, and its linkage with forest demographic processes and carbon dynamics remains elusive. In this study, we develop a theoretical framework and compile a broad-scale dataset of long-term sample forest stands (n = 1,433) from largely undisturbed forests to examine the association of temporal mean tree size vs. density scaling trajectories (slopes) with biomass accumulation rates and the sensitivity of scaling slopes to environmental and demographic drivers. The results empirically demonstrate a large variation of scaling slopes, ranging from -4 to -0.2, across forest stands in tropical, temperate, and boreal forest biomes. Steeper scaling slopes are associated with higher rates of biomass accumulation, resulting from a lower offset of forest growth by biomass loss from mortality. In North America, scaling slopes are positively correlated with forest stand age and rainfall seasonality, thus suggesting a higher rate of biomass accumulation in younger forests with lower rainfall seasonality. These results demonstrate the strong association of the transient mean tree size vs. density scaling trajectories with forest demography and biomass accumulation rates, thus highlighting the potential of leveraging forest structure properties to predict forest demography, carbon fluxes, and dynamics at broad spatial scales.

Non-symmetric responses of leaf onset date to natural warming and cooling in northern ecosystems.

The northern hemisphere has experienced regional cooling, especially during the global warming hiatus (1998-2012) due to ocean energy redistribution. However, the lack of studies about the natural cooling effects hampers our understanding of vegetation responses to climate change. Using 15,125 ground phenological time series at 3,620 sites since the 1950s and 31-year satellite greenness observations (1982-2012) covering the warming hiatus period, we show a stronger response of leaf onset date (LOD) to natural cooling than to warming, i.e. the delay of LOD caused by 1°C cooling is larger than the advance of LOD with 1°C warming. This might be because cooling leads to larger chilling accumulation and heating requirements for leaf onset, but this non-symmetric LOD response is partially offset by warming-related drying. Moreover, spring greening magnitude, in terms of satellite-based greenness and productivity, is more sensitive to LOD changes in the warming area than in the cooling. These results highlight the importance of considering non-symmetric responses of spring greening to warming and cooling when predicting vegetation-climate feedbacks.

Global water use efficiency saturation due to increased vapor pressure deficit.

The ratio of carbon assimilation to water evapotranspiration (ET) of an ecosystem, referred to as ecosystem water use efficiency (WUEeco), is widely expected to increase because of the rising atmospheric carbon dioxide concentration (Ca). However, little is known about the interactive effects of rising Ca and climate change on WUEeco. On the basis of upscaled estimates from machine learning methods and global FLUXNET observations, we show that global WUEeco has not risen since 2001 because of the asymmetric effects of an increased vapor pressure deficit (VPD), which depressed photosynthesis and enhanced ET. An undiminished ET trend indicates that rising temperature and VPD may play a more important role in regulating ET than declining stomatal conductance. Projected increases in VPD are predicted to affect the future coupling of the terrestrial carbon and water cycles.

Carbon uptake in Eurasian boreal forests dominates the high-latitude net ecosystem carbon budget.

Arctic-boreal landscapes are experiencing profound warming, along with changes in ecosystem moisture status and disturbance from fire. This region is of global importance in terms of carbon feedbacks to climate, yet the sign (sink or source) and magnitude of the Arctic-boreal carbon budget within recent years remains highly uncertain. Here, we provide new estimates of recent (2003-2015) vegetation gross primary productivity (GPP), ecosystem respiration (Reco ), net ecosystem CO2 exchange (NEE; Reco - GPP), and terrestrial methane (CH4 ) emissions for the Arctic-boreal zone using a satellite data-driven process-model for northern ecosystems (TCFM-Arctic), calibrated and evaluated using measurements from >60 tower eddy covariance (EC) sites. We used TCFM-Arctic to obtain daily 1-km2 flux estimates and annual carbon budgets for the pan-Arctic-boreal region. Across the domain, the model indicated an overall average NEE sink of -850 Tg CO2 -C year-1 . Eurasian boreal zones, especially those in Siberia, contributed to a majority of the net sink. In contrast, the tundra biome was relatively carbon neutral (ranging from small sink to source). Regional CH4 emissions from tundra and boreal wetlands (not accounting for aquatic CH4 ) were estimated at 35 Tg CH4 -C year-1 . Accounting for additional emissions from open water aquatic bodies and from fire, using available estimates from the literature, reduced the total regional NEE sink by 21% and shifted many far northern tundra landscapes, and some boreal forests, to a net carbon source. This assessment, based on in situ observations and models, improves our understanding of the high-latitude carbon status and also indicates a continued need for integrated site-to-regional assessments to monitor the vulnerability of these ecosystems to climate change.

Respiratory loss during late-growing season determines the net carbon dioxide sink in northern permafrost regions.

Warming of northern high latitude regions (NHL, > 50 °N) has increased both photosynthesis and respiration which results in considerable uncertainty regarding the net carbon dioxide (CO2) balance of NHL ecosystems. Using estimates constrained from atmospheric observations from 1980 to 2017, we find that the increasing trends of net CO2 uptake in the early-growing season are of similar magnitude across the tree cover gradient in the NHL. However, the trend of respiratory CO2 loss during late-growing season increases significantly with increasing tree cover, offsetting a larger fraction of photosynthetic CO2 uptake, and thus resulting in a slower rate of increasing annual net CO2 uptake in areas with higher tree cover, especially in central and southern boreal forest regions. The magnitude of this seasonal compensation effect explains the difference in net CO2 uptake trends along the NHL vegetation- permafrost gradient. Such seasonal compensation dynamics are not captured by dynamic global vegetation models, which simulate weaker respiration control on carbon exchange during the late-growing season, and thus calls into question projections of increasing net CO2 uptake as high latitude ecosystems respond to warming climate conditions.

Lowering water table reduces carbon sink strength and carbon stocks in northern peatlands.

Peatlands at high latitudes have accumulated >400 Pg carbon (C) because saturated soil and cold temperatures suppress C decomposition. This substantial amount of C in Arctic and Boreal peatlands is potentially subject to increased decomposition if the water table (WT) decreases due to climate change, including permafrost thaw-related drying. Here, we optimize a version of the Organizing Carbon and Hydrology In Dynamic Ecosystems model (ORCHIDEE-PCH4) using site-specific observations to investigate changes in CO2 and CH4 fluxes as well as C stock responses to an experimentally manipulated decrease of WT at six northern peatlands. The unmanipulated control peatlands, with the WT <20 cm on average (seasonal max up to 45 cm) below the surface, currently act as C sinks in most years (58 ± 34 g C m-2  year-1 ; including 6 ± 7 g C-CH4 m-2  year-1 emission). We found, however, that lowering the WT by 10 cm reduced the CO2 sink by 13 ± 15 g C m-2  year-1 and decreased CH4 emission by 4 ± 4 g CH4 m-2  year-1 , thus accumulating less C over 100 years (0.2 ± 0.2 kg C m-2 ). Yet, the reduced emission of CH4 , which has a larger greenhouse warming potential, resulted in a net decrease in greenhouse gas balance by 310 ± 360 g CO2-eq  m-2  year-1 . Peatlands with the initial WT close to the soil surface were more vulnerable to C loss: Non-permafrost peatlands lost >2 kg C m-2 over 100 years when WT is lowered by 50 cm, while permafrost peatlands temporally switched from C sinks to sources. These results highlight that reductions in C storage capacity in response to drying of northern peatlands are offset in part by reduced CH4 emissions, thus slightly reducing the positive carbon climate feedbacks of peatlands under a warmer and drier future climate scenario.

Understory plant removal counteracts tree thinning effect on soil respiration in a temperate forest.

Elucidating the response mechanism of soil respiration (Rs) to silvicultural practices is pivotal to evaluating the effects of management practices on soil carbon cycling in planted forest ecosystems. However, as common management practices, how thinning, understory plant removal, and their interactions affect Rs and its autotrophic and heterotrophic components (Ra and Rh) remains unclear. Therefore, we investigated Rs, Ra and Rh by the trenching method from 2011 to 2015 in a Pinus tabuliformis plantation in northern China, subjecting to four treatments (intact control plots [CK], thinning [T], understory removal [UR], and thinning with understory removal [TUR]). Mean annual Rs was significantly increased by thinning (by 15.3%), whereas decreased by UR (by 17.4%), compared with CK. These variations in Rs were mainly attributed to changes in Ra. The increments of Ra were caused by the enhanced growth of fine root biomass after thinning. However, UR led to lower Ra compared with CK (p < .05), indicating that understory growth is inadequate to compensate for the decreased respiring root biomass induced by understory removal. Rs was unchanged between TUR and the intact control plot due to the opposite effects of thinning and UR on the Ra. Changes in Rh exhibited no significant differences among the treatments, partly because of the stable microbial biomass carbon (MBC) and forest floor mass (litter and fine woody debris). No interaction effect between thinning and understory removal was detected on Rs, Ra, and Rh. The lowest temperature sensitivity (Q10 ) value of Ra was found in CK. This study highlights the necessity of incorporating understory plant effects on soil CO2 efflux in assessing forest management practices on soil carbon cycling.

Field-based tree mortality constraint reduces estimates of model-projected forest carbon sinks.

Considerable uncertainty and debate exist in projecting the future capacity of forests to sequester atmospheric CO2. Here we estimate spatially explicit patterns of biomass loss by tree mortality (LOSS) from largely unmanaged forest plots to constrain projected (2015-2099) net primary productivity (NPP), heterotrophic respiration (HR) and net carbon sink in six dynamic global vegetation models (DGVMs) across continents. This approach relies on a strong relationship among LOSS, NPP, and HR at continental or biome scales. The DGVMs overestimated historical LOSS, particularly in tropical regions and eastern North America by as much as 5 Mg ha-1 y-1. The modeled spread of DGVM-projected NPP and HR uncertainties was substantially reduced in tropical regions after incorporating the field-based mortality constraint. The observation-constrained models show a decrease in the tropical forest carbon sink by the end of the century, particularly across South America (from 2 to 1.4 PgC y-1), and an increase in the sink in North America (from 0.8 to 1.1 PgC y-1). These results highlight the feasibility of using forest demographic data to empirically constrain forest carbon sink projections and the potential overestimation of projected tropical forest carbon sinks.

Surface temperatures reveal the patterns of vegetation water stress and their environmental drivers across the tropical Americas.

Vegetation is a key component in the global carbon cycle as it stores ~450 GtC as biomass, and removes about a third of anthropogenic CO2 emissions. However, in some regions, the rate of plant carbon uptake is beginning to slow, largely because of water stress. Here, we develop a new observation-based methodology to diagnose vegetation water stress and link it to environmental drivers. We used the ratio of remotely sensed land surface to near surface atmospheric temperatures (LST/Tair ) to represent vegetation water stress, and built regression tree models (random forests) to assess the relationship between LST/Tair and the main environmental drivers of surface energy fluxes in the tropical Americas. We further determined ecosystem traits associated with water stress and surface energy partitioning, pinpointed critical thresholds for water stress, and quantified changes in ecosystem carbon uptake associated with crossing these critical thresholds. We found that the top drivers of LST/Tair , explaining over a quarter of its local variability in the study region, are (1) radiation, in 58% of the study region; (2) water supply from precipitation, in 30% of the study region; and (3) atmospheric water demand from vapor pressure deficits (VPD), in 22% of the study region. Regions in which LST/Tair variation is driven by radiation are located in regions of high aboveground biomass or at high elevations, while regions in which LST/Tair is driven by water supply from precipitation or atmospheric demand tend to have low species richness. Carbon uptake by photosynthesis can be reduced by up to 80% in water-limited regions when critical thresholds for precipitation and air dryness are exceeded simultaneously, that is, as compound events. Our results demonstrate that vegetation structure and diversity can be important for regulating surface energy and carbon fluxes over tropical regions.

Empirical estimates of regional carbon budgets imply reduced global soil heterotrophic respiration.

Resolving regional carbon budgets is critical for informing land-based mitigation policy. For nine regions covering nearly the whole globe, we collected inventory estimates of carbon-stock changes complemented by satellite estimates of biomass changes where inventory data are missing. The net land-atmospheric carbon exchange (NEE) was calculated by taking the sum of the carbon-stock change and lateral carbon fluxes from crop and wood trade, and riverine-carbon export to the ocean. Summing up NEE from all regions, we obtained a global 'bottom-up' NEE for net land anthropogenic CO2 uptake of -2.2 ± 0.6 PgC yr-1 consistent with the independent top-down NEE from the global atmospheric carbon budget during 2000-2009. This estimate is so far the most comprehensive global bottom-up carbon budget accounting, which set up an important milestone for global carbon-cycle studies. By decomposing NEE into component fluxes, we found that global soil heterotrophic respiration amounts to a source of CO2 of 39 PgC yr-1 with an interquartile of 33-46 PgC yr-1-a much smaller portion of net primary productivity than previously reported.

Decadal variability in land carbon sink efficiency.

BACKGROUND: The climate mitigation target of limiting the temperature increase below 2 °C above the pre-industrial levels requires the efforts from all countries. Tracking the trajectory of the land carbon sink efficiency is thus crucial to evaluate the nationally determined contributions (NDCs). Here, we define the instantaneous land sink efficiency as the ratio of natural land carbon sinks to emissions from fossil fuel and land-use and land-cover change with a value of 1 indicating carbon neutrality to track its temporal dynamics in the past decades. RESULTS: Land sink efficiency has been decreasing during 1957-1990 because of the increased emissions from fossil fuel. After the effect of the Mt. Pinatubo eruption diminished (after 1994), the land sink efficiency firstly increased before 2009 and then began to decrease again after 2009. This reversal around 2009 is mostly attributed to changes in land sinks in tropical regions in response to climate variations. CONCLUSIONS: The decreasing trend of land sink efficiency in recent years reveals greater challenges in climate change mitigation, and that climate impacts on land carbon sinks must be accurately quantified to assess the effectiveness of regional scale climate mitigation policies.

Reconciling carbon-cycle processes from ecosystem to global scales.

Understanding carbon (C) dynamics from ecosystem to global scales remains a challenge. Although expansion of global carbon dioxide (CO2) observatories makes it possible to estimate C-cycle processes from ecosystem to global scales, these estimates do not necessarily agree. At the continental US scale, only 5% of C fixed through photosynthesis remains as net ecosystem exchange (NEE), but ecosystem measurements indicate that only 2% of fixed C remains in grasslands, whereas as much as 30% remains in needleleaf forests. The wet and warm Southeast has the highest gross primary productivity and the relatively wet and cool Midwest has the highest NEE, indicating important spatial mismatches. Newly available satellite and atmospheric data can be combined in innovative ways to identify potential C loss pathways to reconcile these spatial mismatches. Independent datasets compiled from terrestrial and aquatic environments can now be combined to advance C-cycle science across the land-water interface.

Asymmetric responses of ecosystem productivity to rainfall anomalies vary inversely with mean annual rainfall over the conterminous United States.

The CONterminous United States (CONUS) presents a large range of climate conditions and biomes where terrestrial primary productivity and its inter-annual variability are controlled regionally by rainfall and/or temperature. Here, the response of ecosystem productivity to those climate variables was investigated across different biomes from 2010 to 2018 using three climate datasets of precipitation, air temperature or drought severity, combined with several proxies of ecosystem productivity: a remote sensing product of aboveground biomass, an net primary productivity (NPP) remote sensing product, an NPP model-based product and four gross primary productivity products. We used an asymmetry index (AI) where positive AI indicates a greater increase of ecosystem productivity in wet years compared to the decline in dry years, and negative AI indicates a greater decline of ecosystem productivity in dry years compared to the increase in wet years. We found consistent spatial patterns of AI across the CONUS for the different products, with negative asymmetries over the Great Plains and positive asymmetries over the southwestern CONUS. Shrubs and, to a lesser extent, evergreen forests show a persistent positive asymmetry, whilst (natural) grasslands appear to have transitioned from positive to negative anomalies during the last decade. The general tendency of dominant negative asymmetry response for ecosystem productivity across the CONUS appears to be influenced by the negative asymmetry of precipitation anomalies. AI was found to be a function of mean rainfall: more positive AIs were found in dry areas where plants are adapted to drought and take advantage of rainfall pulses, and more negative AIs were found in wet areas, with a threshold delineating the two regimes corresponding to a mean annual rainfall of 200-400 mm/year.

Evolution of woodcutting behaviour in Early Pliocene beaver driven by consumption of woody plants.

Modern beavers (Castor) are prolific ecosystem engineers and dramatically alter the landscape through tree harvesting and dam building. Little is known, however, about the evolutionary drivers of their woodcutting behaviour. Here we investigate if early woodcutting behaviour in Castoridae was driven by nutritional needs. We measured stable carbon and nitrogen isotopes (δ13C and δ15N) of coeval subfossil plants and beaver collagen (Dipoides sp.) from the Early Pliocene, High Arctic Beaver Pond fossil locality (Ellesmere Island), in order to reconstruct Dipoides sp. diet. Isotopic evidence indicates a diet of woody plants and freshwater macrophytes, supporting the hypothesis that this extinct semiaquatic beaver engaged in woodcutting behaviour for feeding purposes. In a phylogenetic context, the isotopic evidence implies that woodcutting and consumption of woody plants can be traced back to a small-bodied, semiaquatic Miocene castorid, suggesting that beavers have been consuming woody plants for over 20 million years. We propose that the behavioural complex (swimming, woodcutting, and consuming woody plants) preceded and facilitated the evolution of dam building. Dam building and food caching behaviours appear to be specializations for cold winter survival and may have evolved in response to late Neogene northern cooling.

Increased high-latitude photosynthetic carbon gain offset by respiration carbon loss during an anomalous warm winter to spring transition.

Arctic and boreal ecosystems play an important role in the global carbon (C) budget, and whether they act as a future net C sink or source depends on climate and environmental change. Here, we used complementary in situ measurements, model simulations, and satellite observations to investigate the net carbon dioxide (CO2 ) seasonal cycle and its climatic and environmental controls across Alaska and northwestern Canada during the anomalously warm winter to spring conditions of 2015 and 2016 (relative to 2010-2014). In the warm spring, we found that photosynthesis was enhanced more than respiration, leading to greater CO2 uptake. However, photosynthetic enhancement from spring warming was partially offset by greater ecosystem respiration during the preceding anomalously warm winter, resulting in nearly neutral effects on the annual net CO2 balance. Eddy covariance CO2 flux measurements showed that air temperature has a primary influence on net CO2 exchange in winter and spring, while soil moisture has a primary control on net CO2 exchange in the fall. The net CO2 exchange was generally more moisture limited in the boreal region than in the Arctic tundra. Our analysis indicates complex seasonal interactions of underlying C cycle processes in response to changing climate and hydrology that may not manifest in changes in net annual CO2 exchange. Therefore, a better understanding of the seasonal response of C cycle processes may provide important insights for predicting future carbon-climate feedbacks and their consequences on atmospheric CO2 dynamics in the northern high latitudes.

Biophysical feedback of global forest fires on surface temperature.

The biophysical feedbacks of forest fire on Earth's surface radiative budget remain uncertain at the global scale. Using satellite observations, we show that fire-induced forest loss accounts for about 15% of global forest loss, mostly in northern high latitudes. Forest fire increases surface temperature by 0.15 K (0.12 to 0.19 K) one year following fire in burned area globally. In high-latitudes, the initial positive climate-fire feedback was mainly attributed to reduced evapotranspiration and sustained for approximately 5 years. Over longer-term (> 5 years), increases in albedo dominated the surface radiative budget resulting in a net cooling effect. In tropical regions, fire had a long-term weaker warming effect mainly due to reduced evaporative cooling. Globally, biophysical feedbacks of fire-induced surface warming one year after fire are equivalent to 62% of warming due to annual fire-related CO2 emissions. Our results suggest that changes in the severity and/or frequency of fire disturbance may have strong impacts on Earth's surface radiative budget and climate, especially at high latitudes.

Data on beetle-killed and surviving lodgepole pine (Pinus contorta) radial growth from the Beaverhead-Deerlodge National Forest, MT prior to a severe mountain pine beetle outbreak.

This article contains measurements of raw radial growth, distance to pith, and calculated basal area increments (BAI) from 444 5-mm increment cores (237 trees) collected in July 2016 from the Beaverhead-Deerlodge National Forest, MT. These data were used for the study presented in "Mountain pine beetle attack faster growing lodgepole pine at low elevations in western Montana, USA" [1]. Plot locations where increment cores were taken as well as code to calculate BAI are also included. Cores were collected from lodgepole pine (Pinus contorta) trees that were killed during a recent bark beetle outbreak (220 cores; 117 trees) as well as trees that survived the outbreak (210 cores; 113 trees) in twelve stands spanning north and south aspects and three elevational bands along a 600-m gradient. 14 additional cores were collected from 7 strip-attacked trees. Increment cores were prepared and measured using standard dendrochronological techniques, "An Introduction to Tree-Ring Dating" [2]. Master chronologies for each aspect-elevation combination were created using approximately ten cores from surviving trees at each location. Cores were cross-dated, then scanned at 2400 dpi. Annual ring widths were measured using CooRecorder 7.7, "Cybis Electronic, CDendro and CooRecorder V.7.7" [3], and final chronologies were quantitatively validated in COFECHA, "Computer-assisted quality control in tree-ring dating and measurement, Tree-Ring Society" [4].

Precipitation thresholds regulate net carbon exchange at the continental scale.

Understanding the sensitivity of ecosystem production and respiration to climate change is critical for predicting terrestrial carbon dynamics. Here we show that the primary control on the inter-annual variability of net ecosystem carbon exchange switches from production to respiration at a precipitation threshold between 750 and 950 mm yr-1 in the contiguous United States. This precipitation threshold is evident across multiple datasets and scales of observation indicating that it is a robust result and provides a new scaling relationship between climate and carbon dynamics. However, this empirical precipitation threshold is not captured by dynamic global vegetation models, which tend to overestimate the sensitivity of production and underestimate the sensitivity of respiration to water availability in more mesic regions. Our results suggest that the short-term carbon balance of ecosystems may be more sensitive to respiration losses than previously thought and that model simulations may underestimate the positive carbon-climate feedbacks associated with respiration.

Limited evidence for CO2 -related growth enhancement in northern Rocky Mountain lodgepole pine populations across climate gradients.

Forests sequester large amounts of carbon annually and are integral in buffering against effects of global change. Increasing atmospheric CO2 may enhance photosynthesis and/or decrease stomatal conductance (gs ) thereby enhancing intrinsic water-use efficiency (iWUE), having potential indirect and direct benefits to tree growth. While increasing iWUE has been observed in most trees globally, enhanced growth is not ubiquitous, possibly due to concurrent climatic constraints on growth. To investigate our incomplete understanding of interactions between climate and CO2 and their impacts on tree physiology and growth, we used an environmental gradient approach. We combined dendrochronology with carbon isotope analysis (δ13 C) to assess the covariation of basal area increment (BAI) and iWUE over time in lodgepole pine. Trees were sampled at 18 sites spanning two climatically distinct elevation transects on the lee and windward sides of the Continental Divide, encompassing the majority of lodgepole pine's northern Rocky Mountain elevational range. We analyzed BAI and iWUE from 1950 to 2015, and explored correlations with monthly climate variables. As expected, iWUE increased at all sites. However, concurrent growth trends depended on site climatic water deficit (CWD). Significant growth increases occurred only at the driest sites, where increases in iWUE were strongest, while growth decreases were greatest at sites where CWD has been historically lowest. Late summer drought of the previous year negatively affected growth across sites. These results suggest that increasing iWUE, if strong enough, may indirectly benefit growth at drier sites by effectively extending the growing season via reductions in gs . Strong growth decreases at high elevation windward sites may reflect increasing water stress as a result of decreasing snowpack, which was not offset by greater iWUE. Our results imply that increasing iWUE driven by decreasing gs may benefit tree growth in limited scenarios, having implications for future carbon uptake potential of semiarid ecosystems.

Future global productivity will be affected by plant trait response to climate.

Plant traits are both responsive to local climate and strong predictors of primary productivity. We hypothesized that future climate change might promote a shift in global plant traits resulting in changes in Gross Primary Productivity (GPP). We characterized the relationship between key plant traits, namely Specific Leaf Area (SLA), height, and seed mass, and local climate and primary productivity. We found that by 2070, tropical and arid ecosystems will be more suitable for plants with relatively lower canopy height, SLA and seed mass, while far northern latitudes will favor woody and taller plants than at present. Using a network of tower eddy covariance CO2 flux measurements and the extrapolated plant trait maps, we estimated the global distribution of annual GPP under current and projected future plant community distribution. We predict that annual GPP in northern biomes (≥45 °N) will increase by 31% (+8.1 ± 0.5 Pg C), but this will be offset by a 17.9% GPP decline in the tropics (-11.8 ± 0.84 Pg C). These findings suggest that regional climate changes will affect plant trait distributions, which may in turn affect global productivity patterns.