SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Abstract The spatial dispersion of photoelements within a vegetation canopy, quantified by the clumping index (CI), directly regulates the within‐canopy light environment and photosynthesis rate, but is not commonly implemented in terrestrial biosphere models to estimate the ecosystem carbon cycle. A few global CI products have been developed recently with remote sensing measurements, making it possible to examine the global impacts of CI. This study deployed CI in the radiative transfer scheme of the Community Land Model version 5 (CLM5) and used the revised CLM5 to quantitatively evaluate the extent to which CI can affect canopy absorbed radiation and gross primary production (GPP), and for the first time, considering the uncertainty and seasonal variation of CI with multiple remote sensing products. Compared to the results without considering the CI impact, the revised CLM5 estimated that sunlit canopy absorbed up to 9%–15% and 23%–34% less direct and diffuse radiation, respectively, while shaded canopy absorbed 3%–18% more diffuse radiation across different biome types. The CI impacts on canopy light conditions included changes in canopy light absorption, and sunlit–shaded leaf area fraction related to nitrogen distribution and thus the maximum rate of Rubisco carboxylase activity ( V cmax ), which together decreased photosynthesis in sunlit canopy by 5.9–7.2 PgC year −1 while enhanced photosynthesis by 6.9–8.2 PgC year −1 in shaded canopy. With higher light use efficiency of shaded leaves, shaded canopy increased photosynthesis compensated and exceeded the lost photosynthesis in sunlit canopy, resulting in 1.0 ± 0.12 PgC year −1 net increase in GPP. The uncertainty of GPP due to the different input CI datasets was much larger than that caused by CI seasonal variations, and was up to 50% of the magnitude of GPP interannual variations in the tropical regions. This study highlights the necessity of considering the impacts of CI and its uncertainty in terrestrial biosphere models.
The continued warming of the Arctic could release vast stores of carbon into the atmosphere from high-latitude ecosystems, especially from thawing permafrost. Increasing uptake of carbon dioxide (CO 2 ) by vegetation during longer growing seasons may partially offset such release of carbon. However, evidence of significant net annual release of carbon from site-level observations and model simulations across tundra ecosystems has been inconclusive. To address this knowledge gap, we combined top-down observations of atmospheric CO 2 concentration enhancements from aircraft and a tall tower, which integrate ecosystem exchange over large regions, with bottom-up observed CO 2 fluxes from tundra environments and found that the Alaska North Slope is not a consistent net source nor net sink of CO 2 to the atmosphere (ranging from ₋6 to +6 Tg C yr -1 for 2012–2017). Our analysis suggests that significant biogenic CO 2 fluxes from unfrozen terrestrial soils, and likely inland waters, during the early cold season (September–December) are major factors in determining the net annual carbon balance of the North Slope, implying strong sensitivity to the rapidly warming freeze-up period. At the regional level, we find no evidence of the previously reported large late-cold-season (January–April) CO 2 emissions to the atmosphere during the study period. Despite the importance of the cold-season CO 2 emissions to the annual total, the interannual variability in the net CO 2 flux is driven by the variability in growing season fluxes. During the growing season, the regional net CO 2 flux is also highly sensitive to the distribution of tundra vegetation types throughout the North Slope. This study shows that quantification and characterization of year-round CO 2 fluxes from the heterogeneous terrestrial and aquatic ecosystems in the Arctic using both site-level and atmospheric observations are important to accurately project the Earth system response to future warming.
The main objective of our research was to characterize biogeochemical cycles at continental and global scales in both terrestrial and aquatic ecosystems. This characterization applied to both natural ecosystems and those disturbed by human activity. The primary elements of interest were carbon and nitrogen and the analysis sought to quantify standing stocks and dynamic cycling processes. The translocation of major nutrients from the terrestrial landscape to the atmosphere (via trace gases) and to fluvial systems (via leaching, erosional losses, and point source pollution) were of particular importance to this study. Our aim was to develop the first generation of Earth System Models. Our research was organized around the construction and testing of component biogeochemical models which treated terrestrial ecosystem processes, aquatic nutrient transport through drainage basins, and trace gas exchanges at the continental and global scale. A suite of three complementary models were defined within this construct. The models were organized to operate at a 1/2 degree latitude by longitude level of spatial resolution and to execute at a monthly time step. This discretization afforded us the opportunity to understand the dynamics of the biosphere down to subregional scales, while simultaneously placing these dynamics into a global context.
Changes in rainfall amounts and patterns have been observed and are expected to continue in the near future with potentially significant ecological and societal consequences. Modelling vegetation responses to changes in rainfall is thus crucial to project water and carbon cycles in the future. In this study, we present the results of a new model-data intercomparison project, where we tested the ability of 10 terrestrial biosphere models to reproduce the observed sensitivity of ecosystem productivity to rainfall changes at 10 sites across the globe, in nine of which, rainfall exclusion and/or irrigation experiments had been performed. The key results are as follows: (a) Inter-model variation is generally large and model agreement varies with timescales. In severely water-limited sites, models only agree on the interannual variability of evapotranspiration and to a smaller extent on gross primary productivity. In more mesic sites, model agreement for both water and carbon fluxes is typically higher on fine (daily–monthly) timescales and reduces on longer (seasonal–annual) scales. (b) Models on average overestimate the relationship between ecosystem productivity and mean rainfall amounts across sites (in space) and have a low capacity in reproducing the temporal (interannual) sensitivity of vegetation productivity to annual rainfall at a given site, even though observation uncertainty is comparable to inter-model variability. (c) Most models reproduced the sign of the observed patterns in productivity changes in rainfall manipulation experiments but had a low capacity in reproducing the observed magnitude of productivity changes. Models better reproduced the observed productivity responses due to rainfall exclusion than addition. (d) All models attribute ecosystem productivity changes to the intensity of vegetation stress and peak leaf area, whereas the impact of the change in growing season length is negligible. Finally, the relative contribution of the peak leaf area and vegetation stress intensity was highly variable among models.
The overall goals of CD-08 team in Phase I were to quantify the contributions of different components of the carbon cycle to overall ecosystem carbon balance in Amazonian tropical forests and to undertake process studies at a number of sites along the eastern LBA transect to understand how and why these fluxes vary with site, season, and year. We divided this work into a number of specific tasks: (1) determining the average rate (and variability) of tree growth over the past 3 decades; (2) determining age demographics of tree populations, using radiocarbon to determine tree age; (3) assessing the rate of production and decomposition of dead wood debris; (4) determining turnover rates for organic matter in soils and the mean age of C respired from soil using radiocarbon measurements; and (5) comparing our results with models and constructing models to predict the potential of tropical forests to function as sources or sinks of C. This report summarizes the considerable progress made towards our original goals, which have led to increased understanding of the potential for central Amazon forests to act as sources or sinks of carbon with altered productivity. The overall picture of tropical forest C dynamics emerging from our Phase I studies suggests that the fraction of gross primary production allocated to growth in these forests is only 25-30%, as opposed to the 50% assumed by many ecosystem models. Consequent slow tree growth rates mean greater mean tree age for a given diameter, as reflected in our measurements and models of tree age. Radiocarbon measurements in leaf and root litter suggest that carbon stays in living tree biomass for several years up to a decade before being added to soils, where decomposition is rapid. The time lags predicted from 14C, when coupled with climate variation on similar time scales, can lead to significant interannual variation in net ecosystem C exchange.
Multiyear trends in Normalized Difference Vegetation Index (NDVI) have been used as metrics of high latitude ecosystem change based on the assumption that NDVI change is associated with ecological change, generally as changes in green vegetation amount (green leaf area index [LAI] or plant cover). Further, no change in NDVI is often interpreted as no change in these variables. Three canopy reflectance models including linear mixture model, the SAIL (Scattering from Arbitrarily Inclined Leaves) model, and the GeoSail model were used to simulate scenarios representing high latitude landscape NDVI responses to changes in LAI and plant cover. The simulations showed inconsistent NDVI responses. Clear increases in NDVI are generally associated with increases in LAI and plant cover. At higher values of LAI, the change in NDVI per unit change in LAI decreases, with very little change in spruce forest NDVI where crown cover is >50% and at the tundra–taiga ecotone with transitions from shrub tundra to spruce woodland. These lower responses may bias the interpretation of greening/browning trends in boreal forests. Variations in water or snow coverage were shown to produce outsized nonbiological NDVI responses. Inconsistencies in NDVI responses exemplify the need for care in the interpretation of NDVI change as a metric of high latitude ecosystem change, and that landscape characteristics in terms of the type of cover and its characteristics, such as the initial plant cover, must be taken into account in evaluating the significance of any observed NDVI trends.
Tropical cyclones are an important cause of forest disturbance, and major storms caused severe structural damage and elevated tree mortality in coastal tropical forests. Model capabilities that can be used to understand post-hurricane forest recovery are still limited. We use a vegetation demography model, the Functionally Assembled Terrestrial Ecosystem Simulator, coupled with the Energy Exascale Earth System Model Land Model (ELM-FATES) to study the processes and the key factors regulating post-hurricane forest recovery. We implemented hurricane-induced forest damage, including defoliation, structural biomass reduction, and tree mortality, performed ensemble model simulations, and used random forest feature importance. For the simulation in the Luquillo Experimental Forest, Puerto Rico, we identified factors controlling the post-hurricane forest recovery, and quantified the sensitivity of key model parameters to the post-hurricane forest recovery. The results indicate a tendency for the Bisley forests to shift toward the light demanding plant functional type (PFT) when the pre-hurricane biomass between the light demanding and shade tolerant PFTs is nearly equal and forests experience hurricane disturbance with mortality >60% for both the two PFTs. Under more realistic conditions where the shade tolerant PFT is initially dominant, mortality >80% is required for a shift toward dominance of the light demanding PFT at Bisley. Hurricane mortality and background mortality are the two major factors regulating post-hurricane forest recovery in simulations. This research improves understanding of the ELM-FATES model behavior associated with hurricane disturbance and provides guidance for dynamic vegetation model development in representing hurricane induced forest damage with varied intensities.
Coastal systems in the northern Gulf of Mexico, including the Mobile Bay, AL estuary, are subject to increasing pressure from a variety of activities including climate change. Climate changes have a direct effect on the discharge of rivers that drain into Mobile Bay and adjacent coastal water bodies. The outflows change water quality (temperature, salinity, and sediment concentrations) in the shallow aquatic areas and affect ecosystem functioning. Mobile Bay is a vital ecosystem that provides habitat for many species of fauna and flora. Historically, submerged aquatic vegetation (SAV) and seagrasses were found in this area of the northern Gulf of Mexico; however the extent of vegetation has significantly decreased over the last 60 years. The objectives of this research are to determine: how climate changes affect runoff and water quality in the estuary and how these changes will affect habitat suitability for SAV and seagrasses. Our approach is to use watershed and hydrodynamic modeling to evaluate the impact of climate change on shallow water aquatic ecosystems in Mobile Bay and adjacent coastal areas. Remotely sensed Landsat data were used for current land cover land use (LCLU) model input and the data provided by Intergovernmental Panel on Climate Change (IPCC) of the future changes in temperature, precipitation, and sea level rise were used to create the climate scenarios for the 2025 and 2050 model simulations. Project results are being shared with Gulf coast stakeholders through the Gulf of Mexico Data Atlas to benefit coastal policy and climate change adaptation strategies.
Heatwaves threaten ecosystem carbon balances, yet the mechanisms driving short-term carbon flux responses remain poorly understood. Here, integrating high-frequency eddy covariance (EC) data from 140 global flux tower sites (872 site-years) with detailed process-based modeling, we examine ecosystem responses during and immediately after heatwaves. We show that heatwaves caused a −40% (range [−29%, −128%]) reduction in net ecosystem productivity (NEP) compared to pre-heatwave values, with this reduction persisting over the following two weeks (−38% range [+3%, −154%]). We attributed NEP decreases to photosynthesis decreases more than to ecosystem respiration (RE) increases. Forest sites had greater NEP decreases during heatwaves than non-forest sites, but remained carbon sinks afterwards, indicating resilience. Our modeling analysis of extreme heatwaves at selected EC sites shows that decreased photosynthesis, increased maintenance respiration, and decreased plant non-structural carbon reserves during heatwaves drive carbon cycle changes that persist for weeks. Consistent with phenocam observations, we modeled a reduction in leaf area index caused by reduced non-structural carbon reserves, leading to early leaf senescence and longer-term impacts. Ongoing increases in heatwaves are therefore likely to reduce NEP across a range of ecosystems, exacerbating carbon cycle feedback.
Abstract Land surface models require continuous validation against observations to improve and reduce simulation uncertainty. However, inferred model performance can be heavily influenced by subjective choices made in the selection and application of observational data products. A key area often misrepresented by models is the Arctic–Boreal region, which is a potential tipping point region in Earth’s climate system due to large permafrost carbon stocks that are vulnerable to release with climate warming. We use the International Land Model Benchmarking (ILAMB) framework to evaluate how the model skill of TRENDY-v9 models varies based on the choice of observational-based benchmark and how benchmarks are applied in model evaluation. This analysis uses global datasets integrated into ILAMB and new, regionally-specific observational products from the Arctic–Boreal Vulnerability Experiment. Our results cover the overall time period of 1979–2019 and show that model scores can vary substantially depending on the data product applied, with higher model scores indicating better model performance against observations. The lowest model scores occur when benchmarked against regional, compared to global, datasets. We also evaluate observed and modeled functional relationships between ecosystem respiration and air temperature and between gross primary production and precipitation. Here, we find that the magnitude and shape of the responses are strongly impacted by the choice of observational dataset and the approach used to construct the functional relationship benchmark. These results suggest that model evaluation studies could conclude a false sense of model skill if only using a single benchmark data product or if not applying regional data products when performing a regional model analysis. Collectively, our findings highlight the influence of benchmarking choices on model evaluation and point to the need for benchmarking guidelines when assessing model skill.
The Pearl River estuary-coastal system in the Northern South China Sea is frequently affected by tropical cyclones (TCs) every year. Previous in-situ and remote sensing studies have found typhoons can enhance phytoplankton biomass and induce blooms in this region. However, the mechanistic links between phytoplankton blooms and typhoons have not been well elucidated due to the interplay of multiple processes along the land-ocean-atmosphere interfaces. Unravelling these interactions will have to rely on an integrated modelling system. In this work, we constructed a realistic, 3-dimensional, land-ocean-atmosphere modelling system with the marine ecosystem and sediment components for the China Great Bay Area. By using the integrated modelling system, we quantitatively investigated phytoplankton response to hydrological conditions variations under Typhoon Hato (2017), a strong typhoon case. Passive tracer experiments showed that with high river discharge induced by heavy rainfall, the residence time of Lingding Bay is as short as 15-day, less than half of that under the climatological discharge. The increase in freshwater pulse washed out the phytoplankton biomass within Lingding Bay. While for the offshore region, the source and sink terms analysis showed that the increase of phytoplankton biomass in the first week was because of the uplift of nutrient-rich subsurface water, while in the second week was because of the seaward propagated nearshore high phytoplankton biomass water. While riverine nutrients support phytoplankton growth in the third week, a large part of phytoplankton biomass was lost to zooplankton grazing, showing the system shifted from the bottom-up control to the top-down control.
Models of photosynthetic production at ecosystem and global scales require multiple input parameters specifying physical and physiological surface features. While certain physical parameters (e.g., absorbed photosynthetically active radiation) can be derived from current satellite sensors, other physiologically relevant measures (e.g., vegetation type, water status, carboxylation capacity, or photosynthetic light-use efficiency), are not generally directly available from current satellite sensors at the appropriate geographic scale. Consequently, many model parameters must be assumed or derived from independent sources, often at an inappropriate scale. An abundance of ecophysiological studies at the leaf and canopy scales suggests strong physiological control of vegetation-atmosphere CO2 and water vapor fluxes, particularly in evergreen vegetation subjected to diurnal or seasonal stresses. For example hot, dry conditions can lead to stomatal closure, and associated "downregulation" of photosynthetic biochemical processes, a phenomenon often manifested as a "midday photosynthetic depression". A recent study with the revised simple biosphere (SiB2) model demonstrated that photosynthetic downregulation can significantly impact global climate. However, at the global scale, the exact significance of downregulation remains unclear, largely because appropriate physiological measures are generally unavailable at this scale. Clearly, there is a need to develop reliable ways of extracting physiologically relevant information from remote sensing. Narrow-band spectrometers offer many opportunities for deriving physiological parameters needed for ecosystem and global scale photosynthetic models. Experimental studies on the ground at the leaf- to stand-scale have indicated that several narrow-band features can be used to detect plant physiological status. One physiological signal is caused by xanthophyll cycle pigment activity, and is often expressed as the Photochemical Reflectance Index (PRI). Because the xanthophyll cycle pigments are photoregulatory pigments closely linked to photosynthetic function, this index can be used to derive relative photosynthetic rates. An additional signal with physiological significance is the 970 nm water absorption band, which provides a measure of liquid water content. This feature has been quantified both using a simple 2-band ratio (900/970 nm, here referred to as the "Water Band Index" or WBI;), and using the "continuum removal" method. Current atmospheric correction methods for AVIRIS imagery also obtain quantitative expressions of surface liquid water absorption based on the 970 nm water band and may be comparable to ground-based estimates of water content using this feature. However, physiological interpretations of both the PRI and the WBI are best understood at the leaf and canopy scales, where complications of atmospheric interference and complex stand and landscape features can be minimized, and where experimental manipulations can be readily applied. Currently it is not known whether these physiological indices can be used to derive meaningful physiological information from AVIRIS imagery. In addition to the problem of atmospheric interference, another challenge is that any simple physiological index can be confounded by multiple factors unrelated to physiology, and this problem can become more severe at progressively larger spatial scales. For example, previous work has suggested that both the PRI and the WBI, are strongly correlated with other optical measures of canopy structure (e.g., the Normalized Difference Vegetation Index or green vegetation fraction), indicating a confounding effect of structure on physiological signals at the larger, landscape scale. Furthermore, the normal operating mode of most imaging spectrometers does not allow simultaneous, ground truthing at a level of detail needed for physiological sampling. Additionally, manipulative experiments of physiology are difficult to apply at a geographic scale suitable for comparison with remote imagery, which often works at spatial scales that are several orders of magnitude larger than those typically used for physiological studies. These limitations require the consideration of alternative approaches to validating physiological information derived from AVIRIS data. In this report, we present a multi-scale sampling approach to detecting physiologically significant signals in narrow-band spectra. This approach explores the multi-dimensional data space provided by narrow-band spectrometry, and combines AVIRIS imagery at a large scale, with ground spectral sampling at an intermediate scale, and detailed ecophysiological measurements at a fine scale, to examine seasonally and spatially changing relationships between multiple structural and physiological variables. Examples of this approach are provided by simultaneous sampling of the Normalized Difference Vegetation Index (NDVI), an index of fractional PAR interception and green vegetation cover, the Water Band Index (WBI, an index of liquid water absorption, and the Photochemical Reflectance Index (PRI, an index of xanthophyll cycle pigment activity and photosynthetic light-use efficiency. By directly linking changing optical properties sampled on the ground with measurable physiological states, we hope to develop a basis for interpreting similar signals in AVIRIS imagery.
Abstract Biomass allocation in plants is the foundation for understanding dynamics in ecosystem carbon balance, species competition, and plant–environment interactions. However, existing work on plant allometry has mainly focused on trees, with fewer studies having developed allometric equations for grasses. Grasses with different life histories can vary in their carbon investment by prioritizing the growth of specific organs to survive, outcompete co‐occurring plants, and ensure population persistence. Further, because grasses are important fuels for wildfire, the lack of grass allocation data adds uncertainty to process‐based models that relate plant physiology to wildfire dynamics. To fill this gap, we conducted a greenhouse experiment with 11 common California grasses varying in photosynthetic pathway and growth form. We measured plant sizes and harvested above‐ and belowground biomass throughout the life cycle of annual species, while for the establishment stage of perennial grasses to quantify allometric relationships for leaf, stem, and root biomass, as well as plant height and canopy area. We used basal diameter as a reference measure of plant size. Overall, basal diameter is the best predictor for leaf and stem biomass, height, and canopy area. Including height as another predictor can improve model accuracy in predicting leaf and stem biomass and canopy area. Fine root biomass is a function of leaf biomass alone. Species vary in their allometric relationships, with most variation occurring for plant height, canopy area, and stem biomass. We further explored potential trade‐offs in biomass allocation across species between leaf and fine root, leaf and stem, and allocation to reproduction. Consistent with our expectation, we found that fast‐growing plants allocated a greater fraction to reproduction. Additionally, plant height and specific leaf area negatively influenced the leaf‐to‐stem ratio. However, contrary to our hypothesis, there were no differences in root‐to‐leaf ratio between perennial and annual or C 4 and C 3 plants. Our study provides species‐specific and functional‐type‐specific allometry equations for both above‐ and belowground organs of 11 common California grass species, enabling nondestructive biomass assessment in California grasslands. These allometric relationships and trade‐offs in carbon allocation across species can improve ecosystem model predictions of grassland species interactions and environmental responses through differences in morphology.
The terrestrial biosphere absorbs about 20% of fossil-fuel CO2 emissions. The overall magnitude of this sink is constrained by the difference between emissions, the rate of increase in atmospheric CO2 concentrations, and the ocean sink. However, the land sink is actually composed of two largely counteracting fluxes that are poorly quantified: fluxes from land-use change andCO2 uptake by terrestrial ecosystems. Dynamic global vegetation model simulations suggest that CO2 emissions from land-use change have been substantially underestimated because processes such as tree harvesting and land clearing from shifting cultivation have not been considered. As the overall terrestrial sink is constrained, a larger net flux as a result of land-use change implies that terrestrial uptake of CO2 is also larger, and that terrestrial ecosystems might have greater potential to sequester carbon in the future. Consequently, reforestation projects and efforts to avoid further deforestation could represent important mitigation pathways, with co-benefits for biodiversity. It is unclear whether a larger land carbon sink can be reconciled with our current understanding of terrestrial carbon cycling. Our possible underestimation of the historical residual terrestrial carbon sink adds further uncertainty to our capacity to predict the future of terrestrial carbon uptake and losses.
Almost half of the global terrestrial soil carbon (C) is stored in the northern circumpolar permafrost region, where air temperatures are increasing two times faster than the global average. As climate warms, permafrost thaws and soil organic matter becomes vulnerable to greater microbial decomposition. Long-term soil warming of ice-rich permafrost can result in thermokarst formation that creates variability in environmental conditions. Consequently, plant and microbial proportional contributions to ecosystem respiration may change in response to long-term soil warming. Natural abundance δ 13 C and Δ 14 C of aboveground and belowground plant material, and of young and old soil respiration were used to inform a mixing model to partition the contribution of each source to ecosystem respiration fluxes. We employed a hierarchical Bayesian approach that incorporated gross primary productivity and environmental drivers to constrain source contributions. We found that long-term experimental permafrost warming introduced a soil hydrology component that interacted with temperature to affect old soil C respiration. Old soil C loss was suppressed in plots with warmer deep soil temperatures because they tended to be wetter. When soil volumetric water content significantly decreased in 2018 relative to 2016 and 2017, the dominant respiration sources shifted from plant aboveground and young soil respiration to old soil respiration. The proportion of ecosystem respiration from old soil C accounted for up to 39% of ecosystem respiration and represented a 30-fold increase compared to the wet-year average. Our findings show that thermokarst formation may act to moderate microbial decomposition of old soil C when soil is highly saturated. However, when soil moisture decreases, a higher proportion of old soil C is vulnerable to decomposition and can become a large flux to the atmosphere. Furthermore, as permafrost systems continue to change with climate, we must understand the thresholds that may propel these systems from a C sink to a source.
Abstract Arctic warming has increased vegetation growth and soil respiration during recent decades. The rate of Arctic warming will likely amplify over the 21st century. Previous studies have revealed that the most severe Arctic warming occurred during the cold season (September to May). The cold-season warming has posited significant CO 2 emissions to the atmosphere via respiration, possibly offsetting warm-season (June to August) net CO 2 uptake. However, prevailing Earth system land models poorly represent cold-season CO 2 emissions, making estimates of Arctic tundra annual CO 2 budgets highly uncertain. Here, we demonstrate that an improved version of the energy exascale Earth system model (E3SM) land model (ELMv1-ECA) captures the large amount of cold-season CO 2 emissions over Alaskan Arctic tundra as reported by two independent, observationally-constrained datasets. We found that the recent seven-decades warming trend of cold-season soil temperature is three times that of the warm-season. The climate sensitivity of warm-season net CO 2 uptake, however, is threefold higher than for the cold-season net CO 2 loss, mainly due to stronger plant resilience than microbial resilience to hydroclimatic extremes. Consequently, the modeled warm-season net CO 2 uptake has a larger positive trend (0.74 ± 0.14 gC m −2 yr −1 ) than that of cold-season CO 2 emissions (0.64 ± 0.11 gC m −2 yr −1 ) from 1950 to 2017, supported by enhanced plant nutrient uptake and increased light- and water-use efficiency. With continued warming and elevated CO 2 concentrations under the representative concentration pathway (RCP) 8.5 scenario, the increasing rate of warm-season net CO 2 uptake is more than twice the rate of cold-season emissions (1.33 ± 0.32 gC m −2 yr −1 vs 0.50 ± 0.12 gC m −2 yr −1 ), making the modeled Alaskan Arctic tundra ecosystem a net CO 2 sink by 2100. However, other geomorphological and ecological disturbances (e.g. abrupt permafrost thaw, thermokarst development, landscape-scale hydrological changes, wildfire, and insects) that are not considered here might alter our conclusion.
The time required for an ecosystem to recover from severe drought is a key component of ecological resilience. The phenology effects on drought recovery are, however, poorly understood. These effects centre on how phenology variations impact biophysical feedbacks, vegetation growth and, ultimately, recovery itself. Using multiple remotely sensed datasets, we found that more than half of ecosystems in mid- and high-latitudinal Northern Hemisphere failed to recover from extreme droughts within a single growing season. Earlier spring phenology in the drought year slowed drought recovery when extreme droughts occurred in mid-growing season. Delayed spring phenology in the subsequent year slowed drought recovery for all vegetation types (with importance of spring phenology ranging from 46% to 58%). The phenology effects on drought recovery were comparable to or larger than other well-known postdrought climatic factors. These results strongly suggest that the interactions between vegetation phenology and drought must be incorporated into Earth system models to accurately quantify ecosystem resilience.