Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “carbon budget”

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.

At least 145 records · Page 8

Dispersal and fire limit Arctic shrub expansion

Abstract Arctic shrub expansion alters carbon budgets, albedo, and warming rates in high latitudes but remains challenging to predict due to unclear underlying controls. Observational studies and models typically use relationships between observed shrub presence and current environmental suitability (bioclimate and topography) to predict shrub expansion, while omitting shrub demographic processes and non-stationary response to changing climate. Here, we use high-resolution satellite imagery across Alaska and western Canada to show that observed shrub expansion has not been controlled by environmental suitability during 1984–2014, but can only be explained by considering seed dispersal and fire. These findings provide the impetus for better observations of recruitment and for incorporating currently underrepresented processes of seed dispersal and fire in land models to project shrub expansion and climate feedbacks. Integrating these dynamic processes with projected fire extent and climate, we estimate shrubs will expand into 25% of the non-shrub tundra by 2100, in contrast to 39% predicted based on increasing environmental suitability alone. Thus, using environmental suitability alone likely overestimates and misrepresents shrub expansion pattern and its associated carbon sink.

54 ENVIRONMENTAL SCIENCES↗

Is there warming in the pipeline? A multi-model analysis of the Zero Emissions Commitment from CO 2

Abstract. The Zero Emissions Commitment (ZEC) is the change in global mean temperature expected to occur following the cessation of net CO 2 emissions and as such is a critical parameter for calculating the remaining carbon budget. The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) was established to gain a better understanding of the potential magnitude and sign of ZEC, in addition to the processes that underlie this metric. A total of 18 Earth system models of both full and intermediate complexity participated in ZECMIP. All models conducted an experiment where atmospheric CO 2 concentration increases exponentially until 1000 PgC has been emitted. Thereafter emissions are set to zero and models are configured to allow free evolution of atmospheric CO 2 concentration. Many models conducted additional second-priority simulations with different cumulative emission totals and an alternative idealized emissions pathway with a gradual transition to zero emissions. The inter-model range of ZEC 50 years after emissions cease for the 1000 PgC experiment is -0.36 to 0.29°C, with a model ensemble mean of -0.07 °C, median of -0.05°C, and standard deviation of 0.19°C. Models exhibit a wide variety of behaviours after emissions cease, with some models continuing to warm for decades to millennia and others cooling substantially. Analysis shows that both the carbon uptake by the ocean and the terrestrial biosphere are important for counteracting the warming effect from the reduction in ocean heat uptake in the decades after emissions cease. This warming effect is difficult to constrain due to high uncertainty in the efficacy of ocean heat uptake. Overall, the most likely value of ZEC on multi-decadal timescales is close to zero, consistent with previous model experiments and simple theory.

MacDougall, Andrew H.↗

Influence of atmospheric oxygen on leaf structure and starch deposition in Arabidopsis thaliana

Plant culture in oxygen concentrations below ambient is known to stimulate vegetative growth, but apart from reports on increased leaf number and weight, little is known about development at subambient oxygen concentrations. Arabidopsis thaliana (L.) Heynh. (cv. Columbia) plants were grown full term in pre-mixed atmospheres with oxygen partial pressures of 2.5, 5.1, 10.1, 16.2, and 21.3 kPa O2, 0.035 kPa CO2 and the balance nitrogen under continuous light. Fully expanded leaves were harvested and processed for light and transmission electron microscopy or for starch quantification. Growth in subambient oxygen concentrations caused changes in leaf anatomy (increased thickness, stomatal density and starch content) that have also been described for plants grown under carbon dioxide enrichment. However, at the lowest oxygen treatment (2.5 kPa), developmental changes occurred that could not be explained by changes in carbon budget caused by suppressed photorespiration, resulting in very thick leaves and a dwarf morphology. This study establishes the leaf parameters that change during growth under low O2, and identifies the lower concentration at which O2 limitation on transport and biosynthetic pathways detrimentally affects leaf development. Grant numbers: NAG5-3756, NAG2-1020, NAG2-1375.

Non-NASA Center↗

Airborne Lidar-Based Estimates of Tropical Forest Structure in Complex Terrain: Opportunities and Trade-Offs for REDD+

Background: Carbon stocks and fluxes in tropical forests remain large sources of uncertainty in the global carbon budget. Airborne lidar remote sensing is a powerful tool for estimating aboveground biomass, provided that lidar measurements penetrate dense forest vegetation to generate accurate estimates of surface topography and canopy heights. Tropical forest areas with complex topography present a challenge for lidar remote sensing. Results: We compared digital terrain models (DTM) derived from airborne lidar data from a mountainous region of the Atlantic Forest in Brazil to 35 ground control points measured with survey grade GNSS receivers. The terrain model generated from full-density (approx. 20 returns/sq m) data was highly accurate (mean signed error of 0.19 +/-0.97 m), while those derived from reduced-density datasets (8/sq m, 4/sq m, 2/sq m and 1/sq m) were increasingly less accurate. Canopy heights calculated from reduced-density lidar data declined as data density decreased due to the inability to accurately model the terrain surface. For lidar return densities below 4/sq m, the bias in height estimates translated into errors of 80-125 Mg/ha in predicted aboveground biomass. Conclusions: Given the growing emphasis on the use of airborne lidar for forest management, carbon monitoring, and conservation efforts, the results of this study highlight the importance of careful survey planning and consistent sampling for accurate quantification of aboveground biomass stocks and dynamics. Approaches that rely primarily on canopy height to estimate aboveground biomass are sensitive to DTM errors from variability in lidar sampling density.

Airborne lidar↗

flat10MIP: an emissions-driven experiment to diagnose the climate response to positive, zero and negative CO2 emissions

Abstract. The proportionality between global mean temperature and cumulative emissions of CO2 predicted in Earth system models (ESMs) is the foundation of carbon budgeting frameworks. Deviations from this behavior could impact estimates of required net-zero timings and negative emissions requirements to meet the Paris Agreement climate targets. However, existing ESM diagnostic experiments do not allow for direct estimation of these deviations as a function of defined emissions pathways. Here, we perform a set of climate model diagnostic experiments for the assessment of transient climate response to cumulative CO2 emissions (TCRE), the Zero Emissions Commitment (ZEC), and climate reversibility metrics in an emissions-driven framework. The emissions-driven experiments provide consistent independent variables simplifying simulation, analysis and interpretation, with emissions rates more comparable to recent levels than existing protocols using model-specific compatible emissions from the CMIP DECK 1pctCO2 experiment, where emissions rates tend to increase during the experiment, such that at the time of CO2 doubling in year 70, emissions are much greater than present-day values. A base experiment, “esm-flat10”, has constant emissions of CO2 of 10 GtC per year (near-present-day values), and initial results show that the TCRE estimated in this experiment is about 0.1 K less than that obtained using 1pctCO2. A subset of ESMs exhibit land carbon sinks that saturate during this experiment. A branch experiment, esm-flat10-zec, illustrates that both positive and negative ZEC effects are less pronounced under esm-flat10 than under 1pctCO2 – the magnitude of ZEC50 in ESMs is, on average, reduced by 30 % compared with 1pctCO2 branch experiments. A final experiment, esm-flat10-cdr, assesses climate reversibility under negative emissions, where we find that peak warming may occur before or after net zero and that the asymmetry in temperature at a given level of cumulative emissions between the positive and negative emissions phases is well described by ZEC in most models. Further, we find that existing probabilistic simple climate model (SCM) ensembles tend to overestimate temperature reversibility compared with ESMs, highlighting the need for additional constraints. We propose a set of climate diagnostic indicators to quantify various aspects of climate reversibility. These experiments were suggested as potential candidates in CMIP7 and have since been adopted as “fast track” simulations.

Sanderson, Benjamin M↗

Estimating CO 2 fluxes through integrating spatial and temporal input layers via deep learning algorithms

Background Accurate estimation of net ecosystem exchange of CO 2 fluxes (Fc) is essential for understanding carbon cycle processes and assessing ecosystem carbon budgets. However, conventional modeling approaches often emphasize temporal dynamics while overlooking the pronounced spatial heterogeneity within the footprint of eddy covariance (EC) towers, potentially limiting predictive accuracy and interpretability of Fc estimates. To address this challenge, we developed a spatiotemporal model that integrates high-resolution footprint-weighted spatial information with sequential environmental drivers. Results The integrated model combines a deeper graph convolutional network to characterize fine-scale spatial variability within EC footprints and a gated recurrent unit network to capture temporal dependencies in biophysical conditions. Using multi-year flux tower observations, remote sensing vegetation indices and footprint modeling, we evaluate the proposed method across three land cover types. This spatiotemporal model consistently outperforms temporal-only and spatial-only baselines, achieving the highest overall accuracy (R 2 = 0.9569) and the lowest RMSE (1.8128 μmol m −2 s −1 ) and MAE (1.1939 μmol m −2 s −1 ). Performance gains are particularly evident in ecosystems with strong vegetation heterogeneity, where spatial structure substantially modulates Fc variability. Conclusions This study demonstrates the importance of joint modeling spatial heterogeneity and temporal dynamics for improving Fc estimation and provides a robust method for advancing footprint-based Fc estimates across diverse ecosystems, supporting refined assessments of terrestrial carbon fluxes, and enhancing scientific foundations for carbon studies.

CO2 flux estimate↗

Damage to living trees contributes to almost half of the biomass losses in tropical forests

Abstract Accurate estimates of forest biomass stocks and fluxes are needed to quantify global carbon budgets and assess the response of forests to climate change. However, most forest inventories consider tree mortality as the only aboveground biomass (AGB) loss without accounting for losses via damage to living trees: branchfall, trunk breakage, and wood decay. Here, we use ~151,000 annual records of tree survival and structural completeness to compare AGB loss via damage to living trees to total AGB loss (mortality + damage) in seven tropical forests widely distributed across environmental conditions. We find that 42% (3.62 Mg ha −1 year −1 ; 95% confidence interval [CI] 2.36–5.25) of total AGB loss (8.72 Mg ha −1 year −1 ; CI 5.57–12.86) is due to damage to living trees. Total AGB loss was highly variable among forests, but these differences were mainly caused by site variability in damage‐related AGB losses rather than by mortality‐related AGB losses. We show that conventional forest inventories overestimate stand‐level AGB stocks by 4% (1%–17% range across forests) because assume structurally complete trees, underestimate total AGB loss by 29% (6%–57% range across forests) due to overlooked damage‐related AGB losses, and overestimate AGB loss via mortality by 22% (7%–80% range across forests) because of the assumption that trees are undamaged before dying. Our results indicate that forest carbon fluxes are higher than previously thought. Damage on living trees is an underappreciated component of the forest carbon cycle that is likely to become even more important as the frequency and severity of forest disturbances increase.

Zuleta, Daniel↗

Estimation of coarse dead wood stocks in intact and degraded forests in the Brazilian Amazon using airborne lidar

Coarse dead wood is an important component of forest carbon stocks, but it is rarely measured in Amazon forests and is typically excluded from regional forest carbon budgets. Our study is based on line intercept sampling for fallen coarse dead wood conducted along 103 transects with a total length of 48 km matched with forest inventory plots where standing coarse dead wood was measured in the footprints of larger areas of airborne lidar acquisitions. We developed models to relate lidar metrics and Landsat time series variables to coarse dead wood stocks for intact, logged, burned, or logged and burned forests. Canopy characteristics such as gap area produced significant individual relations for logged forests. For total fallen plus standing coarse dead wood (hereafter defined as total coarse dead wood), the relative root mean square error for models with only lidar metrics ranged from 33 % in logged forest to up to 36 % in burned forests. The addition of historical information improved model performance slightly for intact forests (31 % against 35 % relative root mean square error), not justifying the use of a number of disturbance events from historical satellite images (Landsat) with airborne lidar data. Lidar-derived estimates of total coarse dead wood compared favorably with independent ground-based sampling for areas up to several hundred hectares. The relations found between total coarse dead wood and variables quantifying forest structure derived from airborne lidar highlight the opportunity to quantify this important but rarely measured component of forest carbon over large areas in tropical forests.

Brazilian Amazon↗

CO2 partial pressure and fluxes in the Amazon River Plume using in situ and remote sensing data

Estimations of the global carbon budget include a quantitative understanding of the evolving processes that occur along river-to-ocean gradients. However, high spatiotemporal resolution observations of these processes are limited. Here we present in situ measurements of the partial pressure of CO2 (pCO2) made through the Amazon River plume (ARP) during different discharge seasons, from 2010-2012. We evaluated the spatiotemporal distribution of pCO2 using Soil Moisture and Ocean Salinity (SMOS) satellite observations for each hydrologic period in the ARP. Regression models were used to estimate pCO2 at the ARP for the period of 2010-2014. From these distributions we calculated sea-air gas exchange of CO2 between the plume waters and the atmosphere (F_co2^sea). Intra-annual variability of F_co2^sea was related to discharge at the river mouth and ocean currents as well as trade winds in the plume. Climatic events during the study period had a significant impact on the F_co2^sea. Including the plume area closer to the river mouth makes the ARP a net source of CO2 with an annual net sea-air flux of 8.6 ±7.1 Tg C y-1 from 2011-2014.

Amazon↗

Substantial hysteresis in emergent temperature sensitivity of global wetland CH4 emissions

Abstract Wetland methane (CH 4 ) emissions ( $${F}_{{{CH}}_{4}}$$ F C H 4 ) are important in global carbon budgets and climate change assessments. Currently, $${F}_{{{CH}}_{4}}$$ F C H 4 projections rely on prescribed static temperature sensitivity that varies among biogeochemical models. Meta-analyses have proposed a consistent $${F}_{{{CH}}_{4}}$$ F C H 4 temperature dependence across spatial scales for use in models; however, site-level studies demonstrate that $${F}_{{{CH}}_{4}}$$ F C H 4 are often controlled by factors beyond temperature. Here, we evaluate the relationship between $${F}_{{{CH}}_{4}}$$ F C H 4 and temperature using observations from the FLUXNET-CH 4 database. Measurements collected across the globe show substantial seasonal hysteresis between $${F}_{{{CH}}_{4}}$$ F C H 4 and temperature, suggesting larger $${F}_{{{CH}}_{4}}$$ F C H 4 sensitivity to temperature later in the frost-free season (about 77% of site-years). Results derived from a machine-learning model and several regression models highlight the importance of representing the large spatial and temporal variability within site-years and ecosystem types. Mechanistic advancements in biogeochemical model parameterization and detailed measurements in factors modulating CH 4 production are thus needed to improve global CH 4 budget assessments.

54 ENVIRONMENTAL SCIENCES↗

Towards a Circular Economy for PET Bottles in the U.S. - 4P Model

The United States generates the most plastic waste of any country. Along with that GHG emissions from the global plastic economy are expected to increase to 15% of the global carbon budget by 2050. It is imperative that plastic recycling is made a reality to reduce both plastic pollution in the environment and GHG emissions. A portfolio of end-of-life strategies must be implemented to minimize environmental impacts and retain valuable plastic material, but it is challenging to compare options that generate products with different utility and lifetime. Plastic use reduction, reuse and recycling are thus increasingly important, but making informed policy and research decisions within this space can be challenging given the diverse range of available solutions. The novel analysis framework, Plastic Parallel Pathways Platform (4P) has been equipped with consequential life cycle assessment, techno-economic analysis, and a plastic circularity indicator to estimate the greenhouse gas (GHG) emissions, circularity, and cost of polyethylene terephthalate (PET) down-cycling to lower-quality resin, closed-loop recycling to food-grade PET bottles, up-cycling to fiber-reinforced plastic (FRP), and conversion to non-plastic products (electricity, oil) on a United States economy-wide basis. Integrating system dynamics into this robust plastics model that already incorporates techno-economics, circularity, and environmental impacts will enable identification of key bottlenecks between manufacturers, waste sorters, and reclaimers that currently prevent rapid decarbonization of the plastics economy. System dynamics (SD) explore the evolution of activities and technologies based on changed macro parameters such as plastic demand and supply, market shifts, and cross-sectoral interactions. This project particularly aims to explore the interplay between, waste collection, plastic waste sorting, recycling, and manufacturing, as well as the effect of plastic bale quality and plastic reuse initiatives on the surrounding process stages. This functionality will facilitate combinatory analysis in which a portfolio of end-of-life pathways are assessed simultaneously, with the exact makeup of that portfolio affected by parameters such as technology scales, resource constraints, and waste mitigation efforts. Integrating SD with the 4P framework enables analyzing the effect of increased revenue and reinvestment into improving process efficiencies, sorting and collection quantities. Through that, market effects of increased recycled resin availability can be studied for the plastics systems model for the US. The results will help identify technical or economic bottlenecks that currently limit efforts to decarbonize the U.S. plastics economy.

carbon↗

Diel variation in CO 2 flux is substantial in many lakes

Lakes play a significant role in the global carbon cycle, acting as sources and sinks of carbon dioxide (CO 2 ). In situ measurements of CO 2 flux (FCO 2 ) from lakes have generally been collected during daylight, despite indications of significant diel variability. This introduces bias when scaling up to whole-lake annual aquatic carbon budgets. We conducted an international sampling program to ascertain the extent of diel variation in FCO 2 across lakes. We sampled 21 lakes over 41 campaigns and measured FCO 2 at 4-h intervals over a full diel cycle. Rates of FCO 2 ranged from −3.16 to 4.39 mmol m −2 h −1 . Integrated over a day, FCO 2 ranged from −381.68 to 878.49 mg C m −2 d −1 (mean = 76.54) across campaigns. We identified three characteristic diel patterns in FCO 2 related to trophic status and show that for half of the campaigns, daily flux estimates were biased by > 50% if based on a single (daytime) measurement.

de Eyto, Elvira [Marine Institute, Mayo (Ireland)]↗

Influence of sun zenith angle on canopy clumping and the resulting impacts on photosynthesis

Addressing the impact of vegetation architecture on shortwave radiation transfer in land surface models is important for accurate weather forecasting, carbon budget estimates, and climate predictions. This paper investigates to what extent it is possible to retrieve structural parameters of two different parameterization schemes from direct transmittance derived from digital hemispherical photography and 3D radiative transfer modeling for two study sites with different vegetation canopy architectures. Neglecting the representation of 3D canopy structure in radiative transfer schemes leads to significant errors in shortwave radiation partitioning (up to 3.5 times more direct transmittance in the 3D model). Structural parameters, referred to as whole canopy ‘clumping indices’, were obtained in order to evaluate the impact of angular variation in clumping on shortwave radiation transfer. Impacts on photosynthesis were evaluated at site level with the UKESM land surface model, JULES. A comparison between flux tower derived and modeled photosynthesis indicates that considering zenith angular variations of structural parameters in the radiative transfer scheme of the UKESM land surface model significantly improves photosynthesis prediction in light limited ecosystems (from RMSE = 2.91 μmol CO 2 .m -2 .s -1 to RMSE = 1.51 μmol CO 2 .m -2 .s -1 , 48% smaller), typically with enhanced photosynthesis from bottom layers.

54 ENVIRONMENTAL SCIENCES↗

Decoupling of greenness and gross primary productivity as aridity decreases

Ecosystem primary productivity is a key ecological process influencing many ecosystem services, including carbon storage. Thus, clarifying how primary productivity in terrestrial ecosystems responds to climatic variability can reveal key mechanisms that will drive future changes in the global carbon budget. Satellite products of canopy greenness are widely used as proxies for vegetation productivity to evaluate how ecosystems respond to climate variability. However, to what degree inter-annual variations in productivity are consistent with greenness and how this relationship varies spatially remains unclear. Here we investigated the strength of the coupling between inter-annual variations in leaf area index (LAI, a measure of greenness) and ecosystem gross primary productivity (GPP) derived from eddy covariance towers, i.e., the r2 of the LAI-GPP relationship. Overall, inter-annual GPP and LAI were highly coupled (i.e., high r2) in arid grasslands, but were fully decoupled in mesic evergreen broadleaf forests, indicating that this relationship varies strongly along aridity gradients. A possible mechanism of the spatial variation in the LAI-GPP relationship is that the tradeoff between ecosystem structure (LAI) and physiology (photosynthesis per unit leaf area) becomes stronger in more humid climates. Land models overestimated the r2 of LAI-GPP correlation for most ecosystem types and failed to capture the spatial pattern along aridity gradients. We conclude that relying on greenness products for evaluating inter-annual changes in vegetation productivity may bias assessments, especially in tropical rainforest ecosystems. Our findings may also reconcile observed disparities between responses in greenness and GPP during drought in Amazon forests.

54 ENVIRONMENTAL SCIENCES↗

Lagged Wetland CH 4 Flux Response in a Historically Wet Year

While a stimulating effect of plant primary productivity on soil carbon dioxide (CO 2 ) emissions has been well documented, links between gross primary productivity (GPP) and wetland methane (CH 4 ) emissions are less well investigated. Determination of the influence of primary productivity on wetland CH 4 emissions (FCH 4 ) is complicated by confounding influences of water table level and temperature on CH 4 production, which also vary seasonally. Here, in this work, we evaluate the link between preceding GPP and subsequent FCH 4 at two fens in Wisconsin using eddy covariance flux towers, Lost Creek (US-Los) and Allequash Creek (US-ALQ). Both wetlands are mosaics of forested and shrub wetlands, with US-Los being larger in scale and having a more open canopy. Co-located sites with multi-year observations of flux, hydrology, and meteorology provide an opportunity to measure and compare lag effects on FCH4 without interference due to differing climate. Daily average FCH4 from US-Los reached a maximum of 47.7 ηmol CH 4 m -2 s -1 during the study period, while US-ALQ was more than double at 117.9 ηmol CH 4 m -2 s -1 . The lagged influence of GPP on temperature-normalized FCH4 (T air -FCH 4 ) was weaker and more delayed in a year with anomalously high precipitation than a following drier year at both sites. FCH 4 at US-ALQ was lower coincident with higher stream discharge in the wet year (2019), potentially due to soil gas flushing during high precipitation events and lower water temperatures. Better understanding of the lagged influence of GPP on FCH 4 due to this study has implications for climate modeling and more accurate carbon budgeting.

59 BASIC BIOLOGICAL SCIENCES↗

Near‐Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub‐Arctic Watershed

Abstract Discontinuous permafrost environments exhibit strong spatial heterogeneity at scales too small to be driven by weather forcing or captured by Earth System Models. Here we analyze effects of observed spatial heterogeneity in soil and vegetation properties, hydrology, and thermal dynamics on ecosystem carbon dynamics in a watershed on the Seward Peninsula in Alaska. We apply a Morris global sensitivity analysis to a process‐rich, successfully tested terrestrial ecosystem model (TEM), ecosys , varying soil properties, boundary conditions, and weather forcing. We show that landscape heterogeneity strongly impacts soil temperatures and vegetation composition. Snow depth, O‐horizon thickness, and near‐surface water content, which vary at scales of O(m), control the soil thermal regime more than an air temperature gradient corresponding to a 140 km north–south distance. High shrub productivity is simulated only in talik (perennially unfrozen) soils with high nitrogen availability. Through these effects on plant and permafrost dynamics, landscape heterogeneity impacts ecosystem productivity. Simulations with near‐surface taliks have higher microbial respiration (by 78.0 gC m −2 yr −1 ) and higher net primary productivity (by 104.9 gC m −2 yr −1 ) compared to runs with near‐surface permafrost, and simulations with high shrub productivity have outlying values of net carbon uptake. We explored the prediction uncertainty associated with ignoring observed landscape heterogeneity, and found that watershed net carbon uptake is 60% larger when heterogeneity is accounted for. Our results highlight the complexity inherent in discontinuous permafrost environments and demonstrate that missing representation of subgrid heterogeneity in TEMs could bias predictions of high‐latitude carbon budget.

54 ENVIRONMENTAL SCIENCES↗

Implementation and Evaluation of Emission‐Driven Land‐Atmosphere Coupled Simulation in E3SMv2.1

Emissions-driven (prognostic CO 2 ) simulations are essential for representing two-way carbon-climate feedback in Earth System Models. We present an emissions-driven land–atmosphere coupled biogeochemistry (BGC) configuration (BGCLNDATM_progCO2) in version 2.1 of the Energy Exascale Earth System Model (E3SMv2.1). This is the first E3SM configuration that performs land-atmosphere emission-hindcasts. Here, we document its implementation, evaluate the model's performance against observations and other models, and propose a structured evaluation protocol for such emissions-driven simulations. We conducted transient historical simulations (1850–2014) with BGCLNDATM_progCO2 and compare them to reference simulations—a land-atmosphere coupled simulation without BGC and a standalone land simulation with BGC, both using prescribed CO 2 concentrations—and to observations. BGCLNDATM_progCO2 overestimates atmospheric CO 2 concentrations by 11–23 ppm yet stays within the 40-ppm spread CMIP6 emission-driven models and retains physical climate properties comparable to the reference runs. The CO 2 biases are partly attributed to underrepresented oceanic CO 2 uptake and inadequate representations of some terrestrial processes. In general, introducing prognostic CO 2 did not change physical climate metrics at the global scale but had larger regional effects, particularly over land where spatially heterogeneous CO 2 and prognostic leaf area index influenced surface energy balance. Finally, we propose a general evaluation protocol including spin-up assessment, atmospheric CO 2 benchmarking, physical climate evaluation, and land biogeochemical analysis to support scientific rigor and facilitate inter-model comparisons. The new configuration lays the groundwork for future enhancements, including improved terrestrial biogeochemical processes, integrated marine biogeochemistry, and additional human–Earth system interactions. These developments advance E3SM toward fully coupled emissions-driven simulations, enabling more accurate carbon–climate feedback projections and informing mitigation policy by providing physically consistent carbon-budget metrics for mitigation scenarios.

54 ENVIRONMENTAL SCIENCES↗

Historically inconsistent productivity and respiration fluxes in the global terrestrial carbon cycle

The terrestrial carbon cycle is a major source of uncertainty in climate projections. Its dominant fluxes, gross primary productivity (GPP), and respiration (in particular soil respiration, R S ), are typically estimated from independent satellite-driven models and upscaled in situ measurements, respectively. We combine carbon-cycle flux estimates and partitioning coefficients to show that historical estimates of global GPP and R S are irreconcilable. When we estimate GPP based on R S measurements and some assumptions about R S :GPP ratios, we found the resulted global GPP values (bootstrap mean ${149}^{+29}_{-23}$ Pg C yr ⁻1 ) are significantly higher than most GPP estimates reported in the literature (${113}^{+18}_{-18}$ Pg C yr ⁻1 ). Similarly, historical GPP estimates imply a soil respiration flux (Rs GPP , bootstrap mean of ${68}^{+10}_{-8}$ Pg C yr ⁻1 ) statistically inconsistent with most published R S values (${87}^{+9}_{-8}$ Pg C yr ⁻1 ), although recent, higher, GPP estimates are narrowing this gap. Furthermore, global R S :GPP ratios are inconsistent with spatial averages of this ratio calculated from individual sites as well as CMIP6 model results. This discrepancy has implications for our understanding of carbon turnover times and the terrestrial sensitivity to climate change. Future efforts should reconcile the discrepancies associated with calculations for GPP and Rs to improve estimates of the global carbon budget.

54 ENVIRONMENTAL SCIENCES↗