Sub-continental-scale carbon stocks of individual trees in African drylands.

The distribution of dryland trees and their density, cover, size, mass and carbon content are not well known at sub-continental to continental scales1-14. This information is important for ecological protection, carbon accounting, climate mitigation and restoration efforts of dryland ecosystems15-18. We assessed more than 9.9 billion trees derived from more than 300,000 satellite images, covering semi-arid sub-Saharan Africa north of the Equator. We attributed wood, foliage and root carbon to every tree in the 0-1,000 mm year-1 rainfall zone by coupling field data19, machine learning20-22, satellite data and high-performance computing. Average carbon stocks of individual trees ranged from 0.54 Mg C ha-1 and 63 kg C tree-1 in the arid zone to 3.7 Mg C ha-1 and 98 kg tree-1 in the sub-humid zone. Overall, we estimated the total carbon for our study area to be 0.84 (±19.8%) Pg C. Comparisons with 14 previous TRENDY numerical simulation studies23 for our area found that the density and carbon stocks of scattered trees have been underestimated by three models and overestimated by 11 models, respectively. This benchmarking can help understand the carbon cycle and address concerns about land degradation24-29. We make available a linked database of wood mass, foliage mass, root mass and carbon stock of each tree for scientists, policymakers, dryland-restoration practitioners and farmers, who can use it to estimate farmland tree carbon stocks from tablets or laptops.

Extensive global wetland loss over the past three centuries.

Wetlands have long been drained for human use, thereby strongly affecting greenhouse gas fluxes, flood control, nutrient cycling and biodiversity1,2. Nevertheless, the global extent of natural wetland loss remains remarkably uncertain3. Here, we reconstruct the spatial distribution and timing of wetland loss through conversion to seven human land uses between 1700 and 2020, by combining national and subnational records of drainage and conversion with land-use maps and simulated wetland extents. We estimate that 3.4 million km2 (confidence interval 2.9-3.8) of inland wetlands have been lost since 1700, primarily for conversion to croplands. This net loss of 21% (confidence interval 16-23%) of global wetland area is lower than that suggested previously by extrapolations of data disproportionately from high-loss regions. Wetland loss has been concentrated in Europe, the United States and China, and rapidly expanded during the mid-twentieth century. Our reconstruction elucidates the timing and land-use drivers of global wetland losses, providing an improved historical baseline to guide assessment of wetland loss impact on Earth system processes, conservation planning to protect remaining wetlands and prioritization of sites for wetland restoration4.

Wetland emission and atmospheric sink changes explain methane growth in 2020.

Atmospheric methane growth reached an exceptionally high rate of 15.1 ± 0.4 parts per billion per year in 2020 despite a probable decrease in anthropogenic methane emissions during COVID-19 lockdowns1. Here we quantify changes in methane sources and in its atmospheric sink in 2020 compared with 2019. We find that, globally, total anthropogenic emissions decreased by 1.2 ± 0.1 teragrams of methane per year (Tg CH4 yr-1), fire emissions decreased by 6.5 ± 0.1 Tg CH4 yr-1 and wetland emissions increased by 6.0 ± 2.3 Tg CH4 yr-1. Tropospheric OH concentration decreased by 1.6 ± 0.2 per cent relative to 2019, mainly as a result of lower anthropogenic nitrogen oxide (NOx) emissions and associated lower free tropospheric ozone during pandemic lockdowns2. From atmospheric inversions, we also infer that global net emissions increased by 6.9 ± 2.1 Tg CH4 yr-1 in 2020 relative to 2019, and global methane removal from reaction with OH decreased by 7.5 ± 0.8 Tg CH4 yr-1. Therefore, we attribute the methane growth rate anomaly in 2020 relative to 2019 to lower OH sink (53 ± 10 per cent) and higher natural emissions (47 ± 16 per cent), mostly from wetlands. In line with previous findings3,4, our results imply that wetland methane emissions are sensitive to a warmer and wetter climate and could act as a positive feedback mechanism in the future. Our study also suggests that nitrogen oxide emission trends need to be taken into account when implementing the global anthropogenic methane emissions reduction pledge5.

Recent increases in annual, seasonal, and extreme methane fluxes driven by changes in climate and vegetation in boreal and temperate wetland ecosystems.

Climate warming is expected to increase global methane (CH4 ) emissions from wetland ecosystems. Although in situ eddy covariance (EC) measurements at ecosystem scales can potentially detect CH4 flux changes, most EC systems have only a few years of data collected, so temporal trends in CH4 remain uncertain. Here, we use established drivers to hindcast changes in CH4 fluxes (FCH4 ) since the early 1980s. We trained a machine learning (ML) model on CH4 flux measurements from 22 [methane-producing sites] in wetland, upland, and lake sites of the FLUXNET-CH4 database with at least two full years of measurements across temperate and boreal biomes. The gradient boosting decision tree ML model then hindcasted daily FCH4 over 1981-2018 using meteorological reanalysis data. We found that, mainly driven by rising temperature, half of the sites (n = 11) showed significant increases in annual, seasonal, and extreme FCH4 , with increases in FCH4 of ca. 10% or higher found in the fall from 1981-1989 to 2010-2018. The annual trends were driven by increases during summer and fall, particularly at high-CH4 -emitting fen sites dominated by aerenchymatous plants. We also found that the distribution of days of extremely high FCH4 (defined according to the 95th percentile of the daily FCH4 values over a reference period) have become more frequent during the last four decades and currently account for 10-40% of the total seasonal fluxes. The share of extreme FCH4 days in the total seasonal fluxes was greatest in winter for boreal/taiga sites and in spring for temperate sites, which highlights the increasing importance of the non-growing seasons in annual budgets. Our results shed light on the effects of climate warming on wetlands, which appears to be extending the CH4 emission seasons and boosting extreme emissions.

COVID-19 lockdowns drive decline in active fires in southeastern United States.

Fire is a common ecosystem process in forests and grasslands worldwide. Increasingly, ignitions are controlled by human activities either through suppression of wildfires or intentional ignition of prescribed fires. The southeastern United States leads the nation in prescribed fire, burning ca. 80% of the country's extent annually. The COVID-19 pandemic radically changed human behavior as workplaces implemented social-distancing guidelines and provided an opportunity to evaluate relationships between humans and fire as fire management plans were postponed or cancelled. Using active fire data from satellite-based observations, we found that in the southeastern United States, COVID-19 led to a 21% reduction in fire activity compared to the 2003 to 2019 average. The reduction was more pronounced for federally managed lands, up to 41% below average compared to the past 20 y (38% below average compared to the past decade). Declines in fire activity were partly affected by an unusually wet February before the COVID-19 shutdown began in mid-March 2020. Despite the wet spring, the predicted number of active fire detections was still lower than expected, confirming a COVID-19 signal on ignitions. In addition, prescribed fire management statistics reported by US federal agencies confirmed the satellite observations and showed that, following the wet February and before the mid-March COVID-19 shutdown, cumulative burned area was approaching record highs across the region. With fire return intervals in the southeastern United States as frequent as 1 to 2 y, COVID-19 fire impacts will contribute to an increasing backlog in necessary fire management activities, affecting biodiversity and future fire danger.

Societal shifts due to COVID-19 reveal large-scale complexities and feedbacks between atmospheric chemistry and climate change.

The COVID-19 global pandemic and associated government lockdowns dramatically altered human activity, providing a window into how changes in individual behavior, enacted en masse, impact atmospheric composition. The resulting reductions in anthropogenic activity represent an unprecedented event that yields a glimpse into a future where emissions to the atmosphere are reduced. Furthermore, the abrupt reduction in emissions during the lockdown periods led to clearly observable changes in atmospheric composition, which provide direct insight into feedbacks between the Earth system and human activity. While air pollutants and greenhouse gases share many common anthropogenic sources, there is a sharp difference in the response of their atmospheric concentrations to COVID-19 emissions changes, due in large part to their different lifetimes. Here, we discuss several key takeaways from modeling and observational studies. First, despite dramatic declines in mobility and associated vehicular emissions, the atmospheric growth rates of greenhouse gases were not slowed, in part due to decreased ocean uptake of CO2 and a likely increase in CH4 lifetime from reduced NO x emissions. Second, the response of O3 to decreased NO x emissions showed significant spatial and temporal variability, due to differing chemical regimes around the world. Finally, the overall response of atmospheric composition to emissions changes is heavily modulated by factors including carbon-cycle feedbacks to CH4 and CO2, background pollutant levels, the timing and location of emissions changes, and climate feedbacks on air quality, such as wildfires and the ozone climate penalty.

A multi-satellite framework to rapidly evaluate extreme biosphere cascades: The Western US 2021 drought and heatwave.

The increasing frequency and intensity of climate extremes and complex ecosystem responses motivate the need for integrated observational studies at low latency to determine biosphere responses and carbon-climate feedbacks. Here, we develop a satellite-based rapid attribution workflow and demonstrate its use at a 1-2-month latency to attribute drivers of the carbon cycle feedbacks during the 2020-2021 Western US drought and heatwave. In the first half of 2021, concurrent negative photosynthesis anomalies and large positive column CO2 anomalies were detected with satellites. Using a simple atmospheric mass balance approach, we estimate a surface carbon efflux anomaly of 132 TgC in June 2021, a magnitude corroborated independently with a dynamic global vegetation model. Integrated satellite observations of hydrologic processes, representing the soil-plant-atmosphere continuum (SPAC), show that these surface carbon flux anomalies are largely due to substantial reductions in photosynthesis because of a spatially widespread moisture-deficit propagation through the SPAC between 2020 and 2021. A causal model indicates deep soil moisture stores partially drove photosynthesis, maintaining its values in 2020 and driving its declines throughout 2021. The causal model also suggests legacy effects may have amplified photosynthesis deficits in 2021 beyond the direct effects of environmental forcing. The integrated, observation framework presented here provides a valuable first assessment of a biosphere extreme response and an independent testbed for improving drought propagation and mechanisms in models. The rapid identification of extreme carbon anomalies and hotspots can also aid mitigation and adaptation decisions.

Spatial heterogeneity of global forest aboveground carbon stocks and fluxes constrained by spaceborne lidar data and mechanistic modeling.

Forest carbon is a large and uncertain component of the global carbon cycle. An important source of complexity is the spatial heterogeneity of vegetation vertical structure and extent, which results from variations in climate, soils, and disturbances and influences both contemporary carbon stocks and fluxes. Recent advances in remote sensing and ecosystem modeling have the potential to significantly improve the characterization of vegetation structure and its resulting influence on carbon. Here, we used novel remote sensing observations of tree canopy height collected by two NASA spaceborne lidar missions, Global Ecosystem Dynamics Investigation and ICE, Cloud, and Land Elevation Satellite 2, together with a newly developed global Ecosystem Demography model (v3.0) to characterize the spatial heterogeneity of global forest structure and quantify the corresponding implications for forest carbon stocks and fluxes. Multiple-scale evaluations suggested favorable results relative to other estimates including field inventory, remote sensing-based products, and national statistics. However, this approach utilized several orders of magnitude more data (3.77 billion lidar samples) on vegetation structure than used previously and enabled a qualitative increase in the spatial resolution of model estimates achievable (0.25° to 0.01°). At this resolution, process-based models are now able to capture detailed spatial patterns of forest structure previously unattainable, including patterns of natural and anthropogenic disturbance and recovery. Through the novel integration of new remote sensing data and ecosystem modeling, this study bridges the gap between existing empirically based remote sensing approaches and process-based modeling approaches. This study more generally demonstrates the promising value of spaceborne lidar observations for advancing carbon modeling at a global scale.

Modeled production, oxidation, and transport processes of wetland methane emissions in temperate, boreal, and Arctic regions.

Wetlands are the largest natural source of methane (CH4 ) to the atmosphere. The eddy covariance method provides robust measurements of net ecosystem exchange of CH4 , but interpreting its spatiotemporal variations is challenging due to the co-occurrence of CH4 production, oxidation, and transport dynamics. Here, we estimate these three processes using a data-model fusion approach across 25 wetlands in temperate, boreal, and Arctic regions. Our data-constrained model-iPEACE-reasonably reproduced CH4 emissions at 19 of the 25 sites with normalized root mean square error of 0.59, correlation coefficient of 0.82, and normalized standard deviation of 0.87. Among the three processes, CH4 production appeared to be the most important process, followed by oxidation in explaining inter-site variations in CH4 emissions. Based on a sensitivity analysis, CH4 emissions were generally more sensitive to decreased water table than to increased gross primary productivity or soil temperature. For periods with leaf area index (LAI) of ≥20% of its annual peak, plant-mediated transport appeared to be the major pathway for CH4 transport. Contributions from ebullition and diffusion were relatively high during low LAI (<20%) periods. The lag time between CH4 production and CH4 emissions tended to be short in fen sites (3 ± 2 days) and long in bog sites (13 ± 10 days). Based on a principal component analysis, we found that parameters for CH4 production, plant-mediated transport, and diffusion through water explained 77% of the variance in the parameters across the 19 sites, highlighting the importance of these parameters for predicting wetland CH4 emissions across biomes. These processes and associated parameters for CH4 emissions among and within the wetlands provide useful insights for interpreting observed net CH4 fluxes, estimating sensitivities to biophysical variables, and modeling global CH4 fluxes.

Observational constraints reduce model spread but not uncertainty in global wetland methane emission estimates.

The recent rise in atmospheric methane (CH4 ) concentrations accelerates climate change and offsets mitigation efforts. Although wetlands are the largest natural CH4 source, estimates of global wetland CH4 emissions vary widely among approaches taken by bottom-up (BU) process-based biogeochemical models and top-down (TD) atmospheric inversion methods. Here, we integrate in situ measurements, multi-model ensembles, and a machine learning upscaling product into the International Land Model Benchmarking system to examine the relationship between wetland CH4 emission estimates and model performance. We find that using better-performing models identified by observational constraints reduces the spread of wetland CH4 emission estimates by 62% and 39% for BU- and TD-based approaches, respectively. However, global BU and TD CH4 emission estimate discrepancies increased by about 15% (from 31 to 36 TgCH4 year-1 ) when the top 20% models were used, although we consider this result moderately uncertain given the unevenly distributed global observations. Our analyses demonstrate that model performance ranking is subject to benchmark selection due to large inter-site variability, highlighting the importance of expanding coverage of benchmark sites to diverse environmental conditions. We encourage future development of wetland CH4 models to move beyond static benchmarking and focus on evaluating site-specific and ecosystem-specific variabilities inferred from observations.

Uncovering the 3d and 4d Electronic Interactions in Solvated Ru Complexes with 2p3d Resonant Inelastic X-ray Scattering.

The electronic structure and dynamics of ruthenium complexes are widely studied given their use in catalytic and light-harvesting materials. Here we investigate three model Ru complexes, [RuIII(NH3)6]3+, [RuII(bpy)3]2+, and [RuII(CN)6]4-, with L3-edge 2p3d resonant inelastic X-ray scattering (RIXS) to probe unoccupied 4d valence orbitals and occupied 3d orbitals and to gain insight into the interactions between these levels. The 2p3d RIXS maps contain a higher level of spectral information than the L3 X-ray absorption near edge structure (XANES). This study provides a direct measure of the 3d spin-orbit splittings of 4.3, 4.0, and 4.1 eV between the 3d5/2 and 3d3/2 orbitals of the [RuIII(NH3)6]3+, [RuII(bpy)3]2+, and [RuII(CN)6]4- complexes, respectively.

Regional trends and drivers of the global methane budget.

The ongoing development of the Global Carbon Project (GCP) global methane (CH4 ) budget shows a continuation of increasing CH4 emissions and CH4 accumulation in the atmosphere during 2000-2017. Here, we decompose the global budget into 19 regions (18 land and 1 oceanic) and five key source sectors to spatially attribute the observed global trends. A comparison of top-down (TD) (atmospheric and transport model-based) and bottom-up (BU) (inventory- and process model-based) CH4 emission estimates demonstrates robust temporal trends with CH4 emissions increasing in 16 of the 19 regions. Five regions-China, Southeast Asia, USA, South Asia, and Brazil-account for >40% of the global total emissions (their anthropogenic and natural sources together totaling >270 Tg CH4  yr-1 in 2008-2017). Two of these regions, China and South Asia, emit predominantly anthropogenic emissions (>75%) and together emit more than 25% of global anthropogenic emissions. China and the Middle East show the largest increases in total emission rates over the 2000 to 2017 period with regional emissions increasing by >20%. In contrast, Europe and Korea and Japan show a steady decline in CH4 emission rates, with total emissions decreasing by ~10% between 2000 and 2017. Coal mining, waste (predominantly solid waste disposal) and livestock (especially enteric fermentation) are dominant drivers of observed emissions increases while declines appear driven by a combination of waste and fossil emission reductions. As such, together these sectors present the greatest risks of further increasing the atmospheric CH4 burden and the greatest opportunities for greenhouse gas abatement.

The terrestrial biosphere as a net source of greenhouse gases to the atmosphere.

The terrestrial biosphere can release or absorb the greenhouse gases, carbon dioxide (CO2), methane (CH4) and nitrous oxide (N2O), and therefore has an important role in regulating atmospheric composition and climate. Anthropogenic activities such as land-use change, agriculture and waste management have altered terrestrial biogenic greenhouse gas fluxes, and the resulting increases in methane and nitrous oxide emissions in particular can contribute to climate change. The terrestrial biogenic fluxes of individual greenhouse gases have been studied extensively, but the net biogenic greenhouse gas balance resulting from anthropogenic activities and its effect on the climate system remains uncertain. Here we use bottom-up (inventory, statistical extrapolation of local flux measurements, and process-based modelling) and top-down (atmospheric inversions) approaches to quantify the global net biogenic greenhouse gas balance between 1981 and 2010 resulting from anthropogenic activities and its effect on the climate system. We find that the cumulative warming capacity of concurrent biogenic methane and nitrous oxide emissions is a factor of about two larger than the cooling effect resulting from the global land carbon dioxide uptake from 2001 to 2010. This results in a net positive cumulative impact of the three greenhouse gases on the planetary energy budget, with a best estimate (in petagrams of CO2 equivalent per year) of 3.9 ± 3.8 (top down) and 5.4 ± 4.8 (bottom up) based on the GWP100 metric (global warming potential on a 100-year time horizon). Our findings suggest that a reduction in agricultural methane and nitrous oxide emissions, particularly in Southern Asia, may help mitigate climate change.

Role of forest regrowth in global carbon sink dynamics.

Although the existence of a large carbon sink in terrestrial ecosystems is well-established, the drivers of this sink remain uncertain. It has been suggested that perturbations to forest demography caused by past land-use change, management, and natural disturbances may be causing a large component of current carbon uptake. Here we use a global compilation of forest age observations, combined with a terrestrial biosphere model with explicit modeling of forest regrowth, to partition the global forest carbon sink between old-growth and regrowth stands over the period 1981-2010. For 2001-2010 we find a carbon sink of 0.85 (0.66-0.96) Pg year-1 located in intact old-growth forest, primarily in the moist tropics and boreal Siberia, and 1.30 (1.03-1.96) Pg year-1 located in stands regrowing after past disturbance. Approaching half of the sink in regrowth stands would have occurred from demographic changes alone, in the absence of other environmental changes. These age-constrained results show consistency with those simulated using an ensemble of demographically-enabled terrestrial biosphere models following an independent reconstruction of historical land use and management. We estimate that forests will accumulate an additional 69 (44-131) Pg C in live biomass from changes in demography alone if natural disturbances, wood harvest, and reforestation continue at rates comparable to those during 1981-2010. Our results confirm that it is not possible to understand the current global terrestrial carbon sink without accounting for the sizeable sink due to forest demography. They also imply that a large portion of the current terrestrial carbon sink is strictly transient in nature.

Forest responses to last-millennium hydroclimate variability are governed by spatial variations in ecosystem sensitivity.

Forecasts of future forest change are governed by ecosystem sensitivity to climate change, but ecosystem model projections are under-constrained by data at multidecadal and longer timescales. Here, we quantify ecosystem sensitivity to centennial-scale hydroclimate variability, by comparing dendroclimatic and pollen-inferred reconstructions of drought, forest composition and biomass for the last millennium with five ecosystem model simulations. In both observations and models, spatial patterns in ecosystem responses to hydroclimate variability are strongly governed by ecosystem sensitivity rather than climate exposure. Ecosystem sensitivity was higher in models than observations and highest in simpler models. Model-data comparisons suggest that interactions among biodiversity, demography and ecophysiology processes dampen the sensitivity of forest composition and biomass to climate variability and change. Integrating ecosystem models with observations from timescales extending beyond the instrumental record can better understand and forecast the mechanisms regulating forest sensitivity to climate variability in a complex and changing world.

Beyond ecosystem modeling: A roadmap to community cyberinfrastructure for ecological data-model integration.

In an era of rapid global change, our ability to understand and predict Earth's natural systems is lagging behind our ability to monitor and measure changes in the biosphere. Bottlenecks to informing models with observations have reduced our capacity to fully exploit the growing volume and variety of available data. Here, we take a critical look at the information infrastructure that connects ecosystem modeling and measurement efforts, and propose a roadmap to community cyberinfrastructure development that can reduce the divisions between empirical research and modeling and accelerate the pace of discovery. A new era of data-model integration requires investment in accessible, scalable, and transparent tools that integrate the expertise of the whole community, including both modelers and empiricists. This roadmap focuses on five key opportunities for community tools: the underlying foundations of community cyberinfrastructure; data ingest; calibration of models to data; model-data benchmarking; and data assimilation and ecological forecasting. This community-driven approach is a key to meeting the pressing needs of science and society in the 21st century.

Plant responses to rising vapor pressure deficit.

Recent decades have been characterized by increasing temperatures worldwide, resulting in an exponential climb in vapor pressure deficit (VPD). VPD has been identified as an increasingly important driver of plant functioning in terrestrial biomes and has been established as a major contributor in recent drought-induced plant mortality independent of other drivers associated with climate change. Despite this, few studies have isolated the physiological response of plant functioning to high VPD, thus limiting our understanding and ability to predict future impacts on terrestrial ecosystems. An abundance of evidence suggests that stomatal conductance declines under high VPD and transpiration increases in most species up until a given VPD threshold, leading to a cascade of subsequent impacts including reduced photosynthesis and growth, and higher risks of carbon starvation and hydraulic failure. Incorporation of photosynthetic and hydraulic traits in 'next-generation' land-surface models has the greatest potential for improved prediction of VPD responses at the plant- and global-scale, and will yield more mechanistic simulations of plant responses to a changing climate. By providing a fully integrated framework and evaluation of the impacts of high VPD on plant function, improvements in forecasting and long-term projections of climate impacts can be made.

A Bornean peat swamp forest is a net source of carbon dioxide to the atmosphere.

Tropical peat forests are a globally important reservoir of carbon, but little is known about CO2 exchange on an annual basis. We measured CO2 exchange between the atmosphere and tropical peat swamp forest in Sarawak, Malaysia using the eddy covariance technique over 4 years from 2011 to 2014. The CO2 fluxes varied between seasons and years. A small carbon uptake took place during the rainy season at the beginning of 2011, while a substantial net efflux of >600 g C/m2 occurred over a 2 month period in the middle of the dry season. Conversely, the peat ecosystem was a source of carbon during both the dry and rainy seasons in subsequent years and more carbon was lost during the rainy season relative to the dry season. Our results demonstrate that the forest was a net source of CO2 to the atmosphere during every year of measurement with annual efflux ranging from 183 to 632 g C m-2  year-1 , noting that annual flux values were sensitive to gap filling methodology. This is in contrast to the typical view of tropical peat forests which must have acted as net C sinks over time scales of centuries to millennia to create the peat deposits. Path analyses revealed that the gross primary productivity (GPP) and ecosystem respiration (RE) were primarily affected by vapour pressure deficit (VPD). Results suggest that future increases in VPD could further reduce the C sink strength and result in additional net CO2 losses from this tropical peat swamp forest in the absence of plant acclimation to such changes in atmospheric dryness.

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.

Critical land change information enhances the understanding of carbon balance in the United States.

Large-scale terrestrial carbon (C) estimating studies using methods such as atmospheric inversion, biogeochemical modeling, and field inventories have produced different results. The goal of this study was to integrate fine-scale processes including land use and land cover change into a large-scale ecosystem framework. We analyzed the terrestrial C budget of the conterminous United States from 1971 to 2015 at 1-km resolution using an enhanced dynamic global vegetation model and comprehensive land cover change data. Effects of atmospheric CO2 fertilization, nitrogen deposition, climate, wildland fire, harvest, and land use/land cover change (LUCC) were considered. We estimate annual C losses from cropland harvest, forest clearcut and thinning, fire, and LUCC were 436.8, 117.9, 10.5, and 10.4 TgC/year, respectively. C stored in ecosystems increased from 119,494 to 127,157 TgC between 1971 and 2015, indicating a mean annual net C sink of 170.3 TgC/year. Although ecosystem net primary production increased by approximately 12.3 TgC/year, most of it was offset by increased C loss from harvest and natural disturbance and increased ecosystem respiration related to forest aging. As a result, the strength of the overall ecosystem C sink did not increase over time. Our modeled results indicate the conterminous US C sink was about 30% smaller than previous modeling studies, but converged more closely with inventory data.