Forest composition change and biophysical climate feedbacks across boreal North America.

Deciduous tree cover is expected to increase in North American boreal forests with climate warming and wildfire. This shift in composition has the potential to generate biophysical cooling via increased land surface albedo. Here we use Landsat-derived maps of continuous tree canopy cover and deciduous fractional composition to assess albedo change over recent decades. We find, on average, a small net decrease in deciduous fraction from 2000 to 2015 across boreal North America and from 1992 to 2015 across Canada, despite extensive fire disturbance that locally increased deciduous vegetation. We further find near-neutral net biophysical change in radiative forcing associated with albedo when aggregated across the domain. Thus, while there have been widespread changes in forest composition over the past several decades, the net changes in composition and associated post-fire radiative forcing have not induced systematic negative feedbacks to climate warming over the spatial and temporal scope of our study.

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.

Influence of atmospheric teleconnections on interannual variability of Arctic-boreal fires.

Fires across the Arctic-boreal zone (ABZ) play an important role in the boreal forest succession, permafrost thaw, and the regional and global carbon cycle and climate. These fires occur mainly in summer with large interannual variability. Previous studies primarily focused on the impacts of local surface climate and tropical El Niño-Southern Oscillation (ENSO). This study, for the first time, comprehensively investigates the influence of summer leading large-scale atmospheric teleconnection patterns in the Northern Hemisphere extra-tropics on interannual variability of ABZ fires. We use correlation and regression analysis of 1997-2019 multiple satellite-based products of burned area and observed/reanalyzed climate data. Results show that eight leading teleconnection patterns significantly affect 63 ± 2 % of burned areas across the ABZ. Western North America is affected by the East Pacific/North Pacific pattern (EP/NP) and the West Pacific pattern (WP); boreal Europe by the Scandinavia pattern (SCA); eastern North America, western and central Siberia, and southeastern Siberia by the North Atlantic Oscillation (NAO); and eastern Siberia /Russian Far East by the East Atlantic pattern (EA). NAO/EA induces lower-tropospheric drier northwesterly/northerly airflow passing through the east of boreal North America/Eurasia, which decreases surface relative humidity. Other teleconnections trigger a high-pressure anomaly, forcing downward motion that suppresses cloud formation and increases solar radiation reaching the ground to warm the surface air as well as brings drier air downward to reduce surface relative humidity. The drier and/or warmer surface air can decrease fuel wetness and thus increase burned area. Our study highlights the important role of the extra-tropical teleconnection patterns on ABZ fires, which is much stronger than ENSO that was thought to control interannual variability of global fires. It also establishes a theoretical foundation for ABZ fire prediction based on extra-tropical teleconnections, and has the potential to facilitate ABZ fire prediction and management.

Escalating carbon emissions from North American boreal forest wildfires and the climate mitigation potential of fire management.

Wildfires in boreal forests release large quantities of greenhouse gases to the atmosphere, exacerbating climate change. Here, we characterize the magnitude of recent and projected gross and net boreal North American wildfire carbon dioxide emissions, evaluate fire management as an emissions reduction strategy, and quantify the associated costs. Our results show that wildfires in boreal North America could, by mid-century, contribute to a cumulative net source of nearly 12 gigatonnes of carbon dioxide, about 3% of remaining global carbon dioxide emissions associated with keeping temperatures within the Paris Agreement's 1.5°C limit. With observations from Alaska, we show that current fire management practices limit the burned area. Further, the costs of avoiding carbon dioxide emissions by means of increasing investment in fire management are comparable to or lower than those of other mitigation strategies. Together, our findings highlight the climate risk that boreal wildfires pose and point to fire management as a cost-effective way to limit emissions.

Climate change, fire return intervals and the growing risk of permanent forest loss in boreal Eurasia.

Climate change has driven an increase in the frequency and severity of fires in Eurasian boreal forests. A growing number of field studies have linked the change in fire regime to post-fire recruitment failure and permanent forest loss. In this study we used four burned area and two forest loss datasets to calculate the landscape-scale fire return interval (FRI) and associated risk of permanent forest loss. We then used machine learning to predict how the FRI will change under a high emissions scenario (SSP3-7.0) by the end of the century. We found that there are currently 133,000 km2 forest at high, or extreme, risk of fire-induced forest loss, with a further 3 M km2 at risk by the end of the century. This has the potential to degrade or destroy some of the largest remaining intact forests in the world, negatively impact the health and economic wellbeing of people living in the region, as well as accelerate global climate change.

Increasing fire and the decline of fire adapted black spruce in the boreal forest.

Intensifying wildfire activity and climate change can drive rapid forest compositional shifts. In boreal North America, black spruce shapes forest flammability and depends on fire for regeneration. This relationship has helped black spruce maintain its dominance through much of the Holocene. However, with climate change and more frequent and severe fires, shifts away from black spruce dominance to broadleaf or pine species are emerging, with implications for ecosystem functions including carbon sequestration, water and energy fluxes, and wildlife habitat. Here, we predict that such reductions in black spruce after fire may already be widespread given current trends in climate and fire. To test this, we synthesize data from 1,538 field sites across boreal North America to evaluate compositional changes in tree species following 58 recent fires (1989 to 2014). While black spruce was resilient following most fires (62%), loss of resilience was common, and spruce regeneration failed completely in 18% of 1,140 black spruce sites. In contrast, postfire regeneration never failed in forests dominated by jack pine, which also possesses an aerial seed bank, or broad-leaved trees. More complete combustion of the soil organic layer, which often occurs in better-drained landscape positions and in dryer duff, promoted compositional changes throughout boreal North America. Forests in western North America, however, were more vulnerable to change due to greater long-term climate moisture deficits. While we find considerable remaining resilience in black spruce forests, predicted increases in climate moisture deficits and fire activity will erode this resilience, pushing the system toward a tipping point that has not been crossed in several thousand years.

Impacts of pre-fire conifer density and wildfire severity on ecosystem structure and function at the forest-tundra ecotone.

Wildfire frequency and extent is increasing throughout the boreal forest-tundra ecotone as climate warms. Understanding the impacts of wildfire throughout this ecotone is required to make predictions of the rate and magnitude of changes in boreal-tundra landcover, its future flammability, and associated feedbacks to the global carbon (C) cycle and climate. We studied 48 sites spanning a gradient from tundra to low-density spruce stands that were burned in an extensive 2013 wildfire on the north slope of the Alaska Range in Denali National Park and Preserve, central Alaska. We assessed wildfire severity and C emissions, and determined the impacts of severity on understory vegetation composition, conifer tree recruitment, and active layer thickness (ALT). We also assessed conifer seed rain and used a seeding experiment to determine factors controlling post-fire tree regeneration. We found that an average of 2.18 ± 1.13 Kg C m-2 was emitted from this fire, almost 95% of which came from burning of the organic soil. On average, burn depth of the organic soil was 10.6 ± 4.5 cm and both burn depth and total C combusted increased with pre-fire conifer density. Sites with higher pre-fire conifer density were also located at warmer and drier landscape positions and associated with increased ALT post-fire, greater changes in pre- and post-fire understory vegetation communities, and higher post-fire boreal tree recruitment. Our seed rain observations and seeding experiment indicate that the recruitment potential of conifer trees is limited by seed availability in this forest-tundra ecotone. We conclude that the expected climate-induced forest infilling (i.e. increased density) at the forest-tundra ecotone could increase fire severity, but this infilling is unlikely to occur without increases in the availability of viable seed.

Overwintering fires in boreal forests.

Forest fires are usually viewed within the context of a single fire season, in which weather conditions and fuel supply can combine to create conditions favourable for fire ignition-usually by lightning or human activity-and spread1-3. But some fires exhibit 'overwintering' behaviour, in which they smoulder through the non-fire season and flare up in the subsequent spring4,5. In boreal (northern) forests, deep organic soils favourable for smouldering6, along with accelerated climate warming7, may present unusually favourable conditions for overwintering. However, the extent of overwintering in boreal forests and the underlying factors influencing this behaviour remain unclear. Here we show that overwintering fires in boreal forests are associated with hot summers generating large fire years and deep burning into organic soils, conditions that have become more frequent in our study areas in recent decades. Our results are based on an algorithm with which we detect overwintering fires in Alaska, USA, and the Northwest Territories, Canada, using field and remote sensing datasets. Between 2002 and 2018, overwintering fires were responsible for 0.8 per cent of the total burned area; however, in one year this amounted to 38 per cent. The spatiotemporal predictability of overwintering fires could be used by fire management agencies to facilitate early detection, which may result in reduced carbon emissions and firefighting costs.

Permafrost carbon feedbacks threaten global climate goals.

Rapid Arctic warming has intensified northern wildfires and is thawing carbon-rich permafrost. Carbon emissions from permafrost thaw and Arctic wildfires, which are not fully accounted for in global emissions budgets, will greatly reduce the amount of greenhouse gases that humans can emit to remain below 1.5 °C or 2 °C. The Paris Agreement provides ongoing opportunities to increase ambition to reduce society's greenhouse gas emissions, which will also reduce emissions from thawing permafrost. In December 2020, more than 70 countries announced more ambitious nationally determined contributions as part of their Paris Agreement commitments; however, the carbon budgets that informed these commitments were incomplete, as they do not fully account for Arctic feedbacks. There is an urgent need to incorporate the latest science on carbon emissions from permafrost thaw and northern wildfires into international consideration of how much more aggressively societal emissions must be reduced to address the global climate crisis.

Statistical upscaling of ecosystem CO2 fluxes across the terrestrial tundra and boreal domain: Regional patterns and uncertainties.

The regional variability in tundra and boreal carbon dioxide (CO2 ) fluxes can be high, complicating efforts to quantify sink-source patterns across the entire region. Statistical models are increasingly used to predict (i.e., upscale) CO2 fluxes across large spatial domains, but the reliability of different modeling techniques, each with different specifications and assumptions, has not been assessed in detail. Here, we compile eddy covariance and chamber measurements of annual and growing season CO2 fluxes of gross primary productivity (GPP), ecosystem respiration (ER), and net ecosystem exchange (NEE) during 1990-2015 from 148 terrestrial high-latitude (i.e., tundra and boreal) sites to analyze the spatial patterns and drivers of CO2 fluxes and test the accuracy and uncertainty of different statistical models. CO2 fluxes were upscaled at relatively high spatial resolution (1 km2 ) across the high-latitude region using five commonly used statistical models and their ensemble, that is, the median of all five models, using climatic, vegetation, and soil predictors. We found the performance of machine learning and ensemble predictions to outperform traditional regression methods. We also found the predictive performance of NEE-focused models to be low, relative to models predicting GPP and ER. Our data compilation and ensemble predictions showed that CO2 sink strength was larger in the boreal biome (observed and predicted average annual NEE -46 and -29 g C m-2  yr-1 , respectively) compared to tundra (average annual NEE +10 and -2 g C m-2  yr-1 ). This pattern was associated with large spatial variability, reflecting local heterogeneity in soil organic carbon stocks, climate, and vegetation productivity. The terrestrial ecosystem CO2 budget, estimated using the annual NEE ensemble prediction, suggests the high-latitude region was on average an annual CO2 sink during 1990-2015, although uncertainty remains high.

Siberian and temperate ecosystems shape Northern Hemisphere atmospheric CO2 seasonal amplification.

The amplitude of the atmospheric CO2 seasonal cycle has increased by 30 to 50% in the Northern Hemisphere (NH) since the 1960s, suggesting widespread ecological changes in the northern extratropics. However, substantial uncertainty remains in the continental and regional drivers of this prominent amplitude increase. Here we present a quantitative regional attribution of CO2 seasonal amplification over the past 4 decades, using a tagged atmospheric transport model prescribed with observationally constrained fluxes. We find that seasonal flux changes in Siberian and temperate ecosystems together shape the observed amplitude increases in the NH. At the surface of northern high latitudes, enhanced seasonal carbon exchange in Siberia is the dominant contributor (followed by temperate ecosystems). Arctic-boreal North America shows much smaller changes in flux seasonality and has only localized impacts. These continental contrasts, based on an atmospheric approach, corroborate heterogeneous vegetation greening and browning trends from field and remote-sensing observations, providing independent evidence for regionally divergent ecological responses and carbon dynamics to global change drivers. Over surface midlatitudes and throughout the midtroposphere, increased seasonal carbon exchange in temperate ecosystems is the dominant contributor to CO2 amplification, albeit with considerable contributions from Siberia. Representing the mechanisms that control the high-latitude asymmetry in flux amplification found in this study should be an important goal for mechanistic land surface models moving forward.

Wildfire combustion and carbon stocks in the southern Canadian boreal forest: Implications for a warming world.

Boreal wildfires are increasing in intensity, extent, and frequency, potentially intensifying carbon emissions and transitioning the region from a globally significant carbon sink to a source. The productive southern boreal forests of central Canada already experience relatively high frequencies of fire, and as such may serve as an analog of future carbon dynamics for more northern forests. Fire-carbon dynamics in southern boreal systems are relatively understudied, with limited investigation into the drivers of pre-fire carbon stocks or subsequent combustion. As part of NASA's Arctic-Boreal Vulnerability Experiment, we sampled 79 stands (47 burned, 32 unburned) throughout central Saskatchewan to characterize above- and belowground carbon stocks and combustion rates in relation to historical land use, vegetation characteristics, and geophysical attributes. We found southern boreal forests emitted an average of 3.3 ± 1.1 kg C/m2 from field sites. The emissions from southern boreal stands varied as a function of stand age, fire weather conditions, ecozone, and soil moisture class. Sites affected by historical timber harvesting had greater combustion rates due to faster carbon stock recovery rates than sites recovering from wildfire events, indicating that different boreal forest land use practices can generate divergent carbon legacy effects. We estimate the 2015 fire season in Saskatchewan emitted a total of 36.3 ± 15.0 Tg C, emphasizing the importance of southern boreal fires for regional carbon budgets. Using the southern boreal as an analog, the northern boreal may undergo fundamental shifts in forest structure and carbon dynamics, becoming dominated by stands <70 years old that hold 2-7 kg C/m2 less than current mature northern boreal stands. Our latitudinal approach reinforces previous studies showing that northern boreal stands are at a high risk of holding less carbon under changing disturbance conditions.

Climate change decreases the cooling effect from postfire albedo in boreal North America.

Fire is a primary disturbance in boreal forests and generates both positive and negative climate forcings. The influence of fire on surface albedo is a predominantly negative forcing in boreal forests, and one of the strongest overall, due to increased snow exposure in the winter and spring months. Albedo forcings are spatially and temporally heterogeneous and depend on a variety of factors related to soils, topography, climate, land cover/vegetation type, successional dynamics, time since fire, season, and fire severity. However, how these variables interact to influence albedo is not well understood, and quantifying these relationships and predicting postfire albedo becomes increasingly important as the climate changes and management frameworks evolve to consider climate impacts. Here we developed a MODIS-derived 'blue sky' albedo product and a novel machine learning modeling framework to predict fire-driven changes in albedo under historical and future climate scenarios across boreal North America. Converted to radiative forcing (RF), we estimated that fires generate an annual mean cooling of -1.77 ± 1.35 W/m2 from albedo under historical climate conditions (1971-2000) integrated over 70 years postfire. Increasing postfire albedo along a south-north climatic gradient was offset by a nearly opposite gradient in solar insolation, such that large-scale spatial patterns in RF were minimal. Our models suggest that climate change will lead to decreases in mean annual postfire albedo, and hence a decreasing strength of the negative RF, a trend dominated by decreased snow cover in spring months. Considering the range of future climate scenarios and model uncertainties, we estimate that for fires burning in the current era (2016) the cooling effect from long-term postfire albedo will be reduced by 15%-28% due to climate change.

Increasing wildfires threaten historic carbon sink of boreal forest soils.

Boreal forest fires emit large amounts of carbon into the atmosphere primarily through the combustion of soil organic matter1-3. During each fire, a portion of this soil beneath the burned layer can escape combustion, leading to a net accumulation of carbon in forests over multiple fire events4. Climate warming and drying has led to more severe and frequent forest fires5-7, which threaten to shift the carbon balance of the boreal ecosystem from net accumulation to net loss1, resulting in a positive climate feedback8. This feedback will occur if organic-soil carbon that escaped burning in previous fires, termed 'legacy carbon', combusts. Here we use soil radiocarbon dating to quantitatively assess legacy carbon loss in the 2014 wildfires in the Northwest Territories of Canada2. We found no evidence for the combustion of legacy carbon in forests that were older than the historic fire-return interval of northwestern boreal forests9. In forests that were in dry landscapes and less than 60 years old at the time of the fire, legacy carbon that had escaped burning in the previous fire cycle was combusted. We estimate that 0.34 million hectares of young forests (<60 years) that burned in the 2014 fires could have experienced legacy carbon combustion. This implies a shift to a domain of carbon cycling in which these forests become a net source-instead of a sink-of carbon to the atmosphere over consecutive fires. As boreal wildfires continue to increase in size, frequency and intensity7, the area of young forests that experience legacy carbon combustion will probably increase and have a key role in shifting the boreal carbon balance.

Expansion of high-latitude deciduous forests driven by interactions between climate warming and fire.

High-latitude regions have experienced rapid warming in recent decades, and this trend is projected to continue over the twenty-first century1. Fire is also projected to increase with warming2,3. We show here, consistent with changes during the Holocene4, that changes in twenty-first century climate and fire are likely to alter the composition of Alaskan boreal forests. We hypothesize that competition for nutrients after fire in early succession and for light in late succession in a warmer climate will cause shifts in plant functional type. Consistent with observations, our ecosystem model predicts evergreen conifers to be the current dominant tree type in Alaska. However, under future climate and fire, our analysis suggests the relative dominance of deciduous broadleaf trees nearly doubles, accounting for 58% of the Alaska ecosystem's net primary productivity by 2100, with commensurate declines in contributions from evergreen conifer trees and herbaceous plants. Post-fire deciduous broadleaf tree growth under a future climate is sustained from enhanced microbial nitrogen mineralization caused by warmer soils and deeper active layers, resulting in taller trees that compete more effectively for light. The expansion of deciduous broadleaf forests will affect the carbon cycle, surface energy fluxes and ecosystem function, thereby modifying important feedbacks with the climate system.

Large loss of CO2 in winter observed across the northern permafrost region.

Recent warming in the Arctic, which has been amplified during the winter1-3, greatly enhances microbial decomposition of soil organic matter and subsequent release of carbon dioxide (CO2)4. However, the amount of CO2 released in winter is highly uncertain and has not been well represented by ecosystem models or by empirically-based estimates5,6. Here we synthesize regional in situ observations of CO2 flux from arctic and boreal soils to assess current and future winter carbon losses from the northern permafrost domain. We estimate a contemporary loss of 1662 Tg C yr-1 from the permafrost region during the winter season (October through April). This loss is greater than the average growing season carbon uptake for this region estimated from process models (-1032 Tg C yr-1). Extending model predictions to warmer conditions in 2100 indicates that winter CO2 emissions will increase 17% under a moderate mitigation scenario-Representative Concentration Pathway (RCP) 4.5-and 41% under business-as-usual emissions scenario-RCP 8.5. Our results provide a new baseline for winter CO2 emissions from northern terrestrial regions and indicate that enhanced soil CO2 loss due to winter warming may offset growing season carbon uptake under future climatic conditions.

Cross-scale controls on carbon emissions from boreal forest megafires.

Climate warming and drying is associated with increased wildfire disturbance and the emergence of megafires in North American boreal forests. Changes to the fire regime are expected to strongly increase combustion emissions of carbon (C) which could alter regional C balance and positively feedback to climate warming. In order to accurately estimate C emissions and thereby better predict future climate feedbacks, there is a need to understand the major sources of heterogeneity that impact C emissions at different scales. Here, we examined 211 field plots in boreal forests dominated by black spruce (Picea mariana) or jack pine (Pinus banksiana) of the Northwest Territories (NWT), Canada after an unprecedentedly large area burned in 2014. We assessed both aboveground and soil organic layer (SOL) combustion, with the goal of determining the major drivers in total C emissions, as well as to develop a high spatial resolution model to scale emissions in a relatively understudied region of the boreal forest. On average, 3.35 kg C m-2 was combusted and almost 90% of this was from SOL combustion. Our results indicate that black spruce stands located at landscape positions with intermediate drainage contribute the most to C emissions. Indices associated with fire weather and date of burn did not impact emissions, which we attribute to the extreme fire weather over a short period of time. Using these results, we estimated a total of 94.3 Tg C emitted from 2.85 Mha of burned area across the entire 2014 NWT fire complex, which offsets almost 50% of mean annual net ecosystem production in terrestrial ecosystems of Canada. Our study also highlights the need for fine-scale estimates of burned area that represent small water bodies and regionally specific calibrations of combustion that account for spatial heterogeneity in order to accurately model emissions at the continental scale.

Detecting early warning signals of tree mortality in boreal North America using multiscale satellite data.

Increasing tree mortality from global change drivers such as drought and biotic infestations is a widespread phenomenon, including in the boreal zone where climate changes and feedbacks to the Earth system are relatively large. Despite the importance for science and management communities, our ability to forecast tree mortality at landscape to continental scales is limited. However, two independent information streams have the potential to inform and improve mortality forecasts: repeat forest inventories and satellite remote sensing. Time series of tree-level growth patterns indicate that productivity declines and related temporal dynamics often precede mortality years to decades before death. Plot-level productivity, in turn, has been related to satellite-based indices such as the Normalized difference vegetation index (NDVI). Here we link these two data sources to show that early warning signals of mortality are evident in several NDVI-based metrics up to 24 years before death. We focus on two repeat forest inventories and three NDVI products across western boreal North America where productivity and mortality dynamics are influenced by periodic drought. These data sources capture a range of forest conditions and spatial resolution to highlight the sensitivity and limitations of our approach. Overall, results indicate potential to use satellite NDVI for early warning signals of mortality. Relationships are broadly consistent across inventories, species, and spatial resolutions, although the utility of coarse-scale imagery in the heterogeneous aspen parkland was limited. Longer-term NDVI data and annually remeasured sites with high mortality levels generate the strongest signals, although we still found robust relationships at sites remeasured at a typical 5 year frequency. The approach and relationships developed here can be used as a basis for improving forest mortality models and monitoring systems.

Vulnerability of eastern US tree species to climate change.

Climate change is expected to alter the distribution of tree species because of critical environmental tolerances related to growth, mortality, reproduction, disturbances, and biotic interactions. How this is realized in 21st century remains uncertain, in large part due to limitations on plant migration and the impacts of landscape fragmentation. Understanding these changes is of particular concern for forest management, which requires information at an appropriately fine spatial resolution. Here we provide a framework and application for tree species vulnerability to climate change in the eastern United States that accounts for influential drivers of future distributions. We used species distribution models to project changes in habitat suitability at 800 m for 40 tree species that vary in physiology, range, and environmental niche. We then developed layers of adaptive capacity based on migration potential, forest fragmentation, and propagule pressure. These were combined into metrics of vulnerability, including an overall index and spatially explicit categories designed to inform management. Despite overall favorable changes in suitability, the majority of species and the landscape were considered vulnerable to climate change. Vulnerability was significantly exacerbated by projections of pests and pathogens for some species. Northern and high-elevation species tended to be the most vulnerable. There were, however, some notable areas of particular resilience, including most of West Virginia. Our approach combines some of the most important considerations for species vulnerability in a straightforward framework, and can be used as a tool for managers to prioritize species, areas, and actions.

Not all droughts are created equal: the impacts of interannual drought pattern and magnitude on grassland carbon cycling.

Climate extremes, such as drought, may have immediate and potentially prolonged effects on carbon cycling. Grasslands store approximately one-third of all terrestrial carbon and may become carbon sources during droughts. However, the magnitude and duration of drought-induced disruptions to the carbon cycle, as well as the mechanisms responsible, remain poorly understood. Over the next century, global climate models predict an increase in two types of drought: chronic but subtle 'press-droughts', and shorter term but extreme 'pulse-droughts'. Much of our current understanding of the ecological impacts of drought comes from experimental rainfall manipulations. These studies have been highly valuable, but are often short term and rarely quantify carbon feedbacks. To address this knowledge gap, we used the Community Land Model 4.0 to examine the individual and interactive effects of pulse- and press-droughts on carbon cycling in a mesic grassland of the US Great Plains. A series of modeling experiments were imposed by varying drought magnitude (precipitation amount) and interannual pattern (press- vs. pulse-droughts) to examine the effects on carbon storage and cycling at annual to century timescales. We present three main findings. First, a single-year pulse-drought had immediate and prolonged effects on carbon storage due to differential sensitivities of ecosystem respiration and gross primary production. Second, short-term pulse-droughts caused greater carbon loss than chronic press-droughts when total precipitation reductions over a 20-year period were equivalent. Third, combining pulse- and press-droughts had intermediate effects on carbon loss compared to the independent drought types, except at high drought levels. Overall, these results suggest that interannual drought pattern may be as important for carbon dynamics as drought magnitude and that extreme droughts may have long-lasting carbon feedbacks in grassland ecosystems.