Search NASA⌕ Search

SEARCH · Search NASA

Results for “FOSSIL”

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 1,117 records · Page 62

THE EATING HABITS OF MILKY WAY-MASS HALOS: DESTROYED DWARF SATELLITES AND THE METALLICITY DISTRIBUTION OF ACCRETED STARS

We study the mass spectrum of destroyed dwarfs that contribute to the accreted stellar mass of Milky Way (MW)-mass (M{sub vir} ∼ 10{sup 12.1} M{sub ⊙}) halos using a suite of 45 zoom-in dissipationless simulations. Empirical models are employed to relate (peak) subhalo mass to dwarf stellar mass, and we use constraints from z = 0 observations and hydrodynamical simulations to estimate the metallicity distribution of the accreted stellar material. The dominant contributors to the accreted stellar mass are relatively massive dwarfs with M{sub star} ∼ 10{sup 8}–10{sup 10}M{sub ⊙}. Halos with more quiescent accretion histories tend to have lower mass progenitors (10{sup 8}–10{sup 9} M{sub ⊙}), and lower overall accreted stellar masses. Ultra-faint mass (M{sub star} < 10{sup 5} M{sub ⊙}) dwarfs contribute a negligible amount (≪1%) to the accreted stellar mass and, despite having low average metallicities, supply a small fraction (∼2%–5%) of the very metal-poor stars with [Fe/H] < −2. Dwarfs with masses 10{sup 5} < M{sub star}/M{sub ⊙} < 10{sup 8} provide a substantial amount of the very metal-poor stellar material (∼40%–80%), and even relatively metal-rich dwarfs with M{sub star} > 10{sup 8} M{sub ⊙} can contribute a considerable fraction (∼20%–60%) of metal-poor stars if their metallicity distributions have significant metal-poor tails. Finally, we find that the generic assumption of a quiescent assembly history for the MW halo seems to be in tension with the mass spectrum of its surviving dwarfs. We suggest that the MW could be a “transient fossil”; a quiescent halo with a recent accretion event(s) that disguises the preceding formation history of the halo.

79 ASTRONOMY AND ASTROPHYSICS↗

THE ACS LCID PROJECT. XI. ON THE EARLY TIME RESOLUTION OF SFHs OF LOCAL GROUP DWARF GALAXIES: COMPARING THE EFFECTS OF REIONIZATION IN MODELS WITH OBSERVATIONS

The analysis of the early star formation history (SFH) of nearby galaxies, obtained from their resolved stellar populations, is relevant as a test for cosmological models. However, the early time resolution of observationally derived SFHs is limited by several factors. Thus, direct comparison of observationally derived SFHs with those derived from theoretical models of galaxy formation is potentially biased. Here we investigate and quantify this effect. For this purpose, we analyze the duration of the early star formation activity in a sample of four Local Group dwarf galaxies and test whether they are consistent with being true fossils of the pre-reionization era; i.e., if the quenching of their star formation occurred before cosmic reionization by UV photons was completed. Two classical dSph (Cetus and Tucana) and two dTrans (LGS-3 and Phoenix) isolated galaxies with total stellar masses between 1.3×10{sup 6} and 7.2×10{sup 6} M{sub ⊙} have been studied. Accounting for time resolution effects, the SFHs peak as much as 1.25 Gyr earlier than the optimal solutions. Thus, this effect is important for a proper comparison of model and observed SFHs. It is also shown that none of the analyzed galaxies can be considered a true fossil of the pre-reionization era, although it is possible that the outer regions of Cetus and Tucana are consistent with quenching by reionization.

79 ASTRONOMY AND ASTROPHYSICS↗

Chandra and JVLA Observations of HST Frontier Fields Cluster MACS J0717.5+3745

To investigate the relationship between thermal and non-thermal components in merger galaxy clusters, we present deep JVLA and Chandra observations of the HST Frontier Fields cluster MACS J0717.5+3745. The Chandra image shows a complex merger event, with at least four components belonging to different merging subclusters. Northwest of the cluster, ∼0.7 Mpc from the center, there is a ram-pressure-stripped core that appears to have traversed the densest parts of the cluster after entering the intracluster medium (ICM) from the direction of a galaxy filament to the southeast. We detect a density discontinuity north-northeast of this core, which we speculate is associated with a cold front. Our radio images reveal new details for the complex radio relic and radio halo in this cluster. In addition, we discover several new filamentary radio sources with sizes of 100–300 kpc. A few of these seem to be connected to the main radio relic, while others are either embedded within the radio halo or projected onto it. A narrow-angled-tailed (NAT) radio galaxy, a cluster member, is located at the center of the radio relic. The steep spectrum tails of this active galactic nucleus lead into the large radio relic where the radio spectrum flattens again. This morphological connection between the NAT radio galaxy and relic provides evidence for re-acceleration (revival) of fossil electrons. The presence of hot ≳20 keV ICM gas detected by Chandra near the relic location provides additional support for this re-acceleration scenario.

79 ASTRONOMY AND ASTROPHYSICS↗

Deep Very Large Array Observations of the Merging Cluster CIZA J2242.8+5301: Continuum and Spectral Imaging

Despite the progress that has been made in understanding radio relics, there are still open questions regarding the underlying particle acceleration mechanisms. In this paper, we present deep 1–4 GHz Very Large Array (VLA) observations of CIZA J2242.8+5301 ( z = 0.1921), a double radio relic cluster characterized by small projection on the plane of the sky. Our VLA observations reveal, for the first time, the complex morphology of the diffuse sources and the filamentary structure of the northern relic. Additionally, we discover new, faint, diffuse radio emission extending north of the main northern relic. Our Mach number estimates for the northern and southern relics, based on the radio spectral index map obtained using the VLA observations and existing LOw Frequency ARray (LOFAR) and Giant Metrewave Radio Telescope data, are consistent with previous radio and X-ray studies ( and ). However, color–color diagrams and models suggest a flatter injection spectral index than the one obtained from the spectral index map, indicating that projection effects might be not entirely negligible. The southern relic consists of five "arms." Embedded in it, we find a tailed radio galaxy that seems to be connected to the relic. A spectral index flattening, where the radio tail connects to the relic, is also measured. We propose that the southern relic may trace AGN fossil electrons that are reaccelerated by a shock, with an estimated strength of . High-resolution mapping of other tailed radio galaxies also supports a scenario where AGN fossil electrons are revived by the merger event and could be related to the formation of some diffuse cluster radio emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Creating Accurate Methane Emission Inventories through Data-Driven Airborne Survey Strategies

Because natural gas emits less carbon than other fossil fuels, it holds promise as a green energy transition fuel. However, the overall carbon footprint of natural gas is significantly elevated by methane emissions that occur during its production and transmission (Cusworth et al. 2022). Methane “super-emitters,” while comprising only about 1% of sites, are responsible for the majority of oil- and gas-sourced methane emissions, making their detection and mitigation critical in reducing the climate impact of natural gas and in meeting national and global sustainability goals (Sherwin et al. 2024). Yet, despite advancements in detection, significant uncertainties remain regarding the size, frequency, and duration distributions of methane emissions (e.g., Frankenberg et al. 2016, Cusworth et al. 2022, Chen, Sherwin et al. 2022, Conrad et al. 2023, Johnson et al. 2023, Sherwin et al. 2024) underscoring the need for comprehensive emissions inventories segmented by basin across the US. Airborne surveys are well-suited for collecting data to build these comprehensive, basin-level inventories because they allow for extensive spatial coverage, and have the spatial resolution, and the sensitivity to pinpoint individual methane sources. As remote sensing technologies enable rapid basin-scale surveys, it is imperative to establish scientifically and statistically robust standards to generate reliable and actionable emissions inventories. Recent work has shown that differences in airborne sampling strategies, detection technologies, and analysis can lead to large differences between survey conclusions if not correctly accounted for (Chen et al. 2024). This elevates the importance of incorporating proper sampling and analysis techniques when designing a methane emissions monitoring campaign to produce accurate results and facilitate cross-study comparisons. In this paper, we describe a survey strategy designed using the latest conclusions from the literature to align results from different aerial surveys. We identify several sampling and analysis principles, including large sample sizes, balanced sampling across oil and gas production, careful survey area definition, and a unified protocol for analysis, to be vital to producing an unbiased estimate of basin-scale emissions. We present results from a Department of Energy-funded project that deployed this survey strategy in two understudied oil and gas- producing regions in the United States: the Haynesville Basin in Texas and Louisiana, and the Woodford Shale in the Anadarko Basin in Oklahoma.

03 NATURAL GAS↗

Radiative forcing of the direct aerosol effect from AeroCom Phase II simulations

We report on the AeroCom Phase II direct aerosol effect (DAE) experiment where 16 detailed global aerosol models have been used to simulate the changes in the aerosol distribution over the industrial era. All 16 models have estimated the radiative forcing (RF) of the anthropogenic DAE, and have taken into account anthropogenic sulphate, black carbon (BC) and organic aerosols (OA) from fossil fuel, biofuel, and biomass burning emissions. In addition several models have simulated the DAE of anthropogenic nitrate and anthropogenic influenced secondary organic aerosols (SOA). The model simulated all-sky RF of the DAE from total anthropogenic aerosols has a range from -0.58 to -0.02 Wm -2 , with a mean of -0.27 Wm -2 for the 16 models. Several models did not include nitrate or SOA and modifying the estimate by accounting for this with information from the other AeroCom models reduces the range and slightly strengthens the mean. Modifying the model estimates for missing aerosol components and for the time period 1750 to 2010 results in a mean RF for the DAE of -0.35 Wm -2 . Compared to AeroCom Phase I (Schulz et al., 2006) we find very similar spreads in both total DAE and aerosol component RF. However, the RF of the total DAE is stronger negative and RF from BC from fossil fuel and biofuel emissions are stronger positive in the present study than in the previous AeroCom study. We find a tendency for models having a strong (positive) BC RF to also have strong (negative) sulphate or OA RF. This relationship leads to smaller uncertainty in the total RF of the DAE compared to the RF of the sum of the individual aerosol components. The spread in results for the individual aerosol components is substantial, and can be divided into diversities in burden, mass extinction coefficient (MEC), and normalized RF with respect to AOD. We find that these three factors give similar contributions to the spread in results.

54 ENVIRONMENTAL SCIENCES↗

Modeling the formation and aging of secondary organic aerosols in Los Angeles during CalNex 2010

Four different literature parameterizations for the formation and evolution of urban secondary organic aerosol (SOA) frequently used in 3-D models are evaluated using a 0-D box model representing the Los Angeles metropolitan region during the California Research at the Nexus of Air Quality and Climate Change (CalNex) 2010 campaign. We constrain the model predictions with measurements from several platforms and compare predictions with particle- and gas-phase observations from the CalNex Pasadena ground site. That site provides a unique opportunity to study aerosol formation close to anthropogenic emission sources with limited recirculation. The model SOA that formed only from the oxidation of VOCs (V-SOA) is insufficient to explain the observed SOA concentrations, even when using SOA parameterizations with multi-generation oxidation that produce much higher yields than have been observed in chamber experiments, or when increasing yields to their upper limit estimates accounting for recently reported losses of vapors to chamber walls. The Community Multiscale Air Quality (WRF-CMAQ) model (version 5.0.1) provides excellent predictions of secondary inorganic particle species but underestimates the observed SOA mass by a factor of 25 when an older VOC-only parameterization is used, which is consistent with many previous model–measurement comparisons for pre-2007 anthropogenic SOA modules in urban areas. Including SOA from primary semi-volatile and intermediate-volatility organic compounds (P-S/IVOCs) following the parameterizations of Robinson et al. (2007), Grieshop et al. (2009), or Pye and Seinfeld (2010) improves model–measurement agreement for mass concentration. The results from the three parameterizations show large differences (e.g., a factor of 3 in SOA mass) and are not well constrained, underscoring the current uncertainties in this area. Our results strongly suggest that other precursors besides VOCs, such as P-S/IVOCs, are needed to explain the observed SOA concentrations in Pasadena. All the recent parameterizations overpredict urban SOA formation at long photochemical ages (≈ 3 days) compared to observations from multiple sites, which can lead to problems in regional and especially global modeling. However, reducing IVOC emissions by one-half in the model to better match recent IVOC measurements improves SOA predictions at these long photochemical ages. Among the explicitly modeled VOCs, the precursor compounds that contribute the greatest SOA mass are methylbenzenes. Measured polycyclic aromatic hydrocarbons (naphthalenes) contribute 0.7% of the modeled SOA mass. The amounts of SOA mass from diesel vehicles, gasoline vehicles, and cooking emissions are estimated to be 16–27, 35–61, and 19–35%, respectively, depending on the parameterization used, which is consistent with the observed fossil fraction of urban SOA, 71(±3) %. The relative contribution of each source is uncertain by almost a factor of 2 depending on the parameterization used. In-basin biogenic VOCs are predicted to contribute only a few percent to SOA. A regional SOA background of approximately 2.1 μg m −3 is also present due to the long-distance transport of highly aged OA, likely with a substantial contribution from regional biogenic SOA. The percentage of SOA from diesel vehicle emissions is the same, within the estimated uncertainty, as reported in previous work that analyzed the weekly cycles in OA concentrations (Bahreini et al., 2012; Hayes et al., 2013). However, the modeling work presented here suggests a strong anthropogenic source of modern carbon in SOA, due to cooking emissions, which was not accounted for in those previous studies and which is higher on weekends. Lastly, this work adapts a simple two-parameter model to predict SOA concentration and O/C from urban emissions. This model successfully predicts SOA concentration, and the optimal parameter combination is very similar to that found for Mexico City. This approach provides a computationally inexpensive method for predicting urban SOA in global and climate models. We estimate pollution SOA to account for 26 Tg yr −1 of SOA globally, or 17% of global SOA, one-third of which is likely to be non-fossil.

54 ENVIRONMENTAL SCIENCES↗

Variability and quasi-decadal changes in the methane budget over the period 2000–2012

Following the recent Global Carbon Project (GCP) synthesis of the decadal methane (CH 4 ) budget over 2000–2012 (Saunois et al., 2016), we analyse here the same dataset with a focus on quasi-decadal and inter-annual variability in CH 4 emissions. The GCP dataset integrates results from top-down studies (exploiting atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up models (including process-based models for estimating land surface emissions and atmospheric chemistry), inventories of anthropogenic emissions, and data-driven approaches. The annual global methane emissions from top-down studies, which by construction match the observed methane growth rate within their uncertainties, all show an increase in total methane emissions over the period 2000-2012, but this increase is not linear over the 13 years. Despite differences between individual studies, the mean emission anomaly of the top-down ensemble shows no significant trend in total methane emissions over the period 2000-2006, during the plateau of atmospheric methane mole fractions, and also over the period 2008-2012, during the renewed atmospheric methane increase. However, the top-down ensemble mean produces an emission shift between 2006 and 2008, leading to 22 [16-32]Tg CH 4 yr -1 higher methane emissions over the period 2008-2012 compared to 2002-2006. This emission increase mostly originated from the tropics, with a smaller contribution from mid-latitudes and no significant change from boreal regions. The regional contributions remain uncertain in top-down studies. Tropical South America and South and East Asia seem to contribute the most to the emission increase in the tropics. However, these two regions have only limited atmospheric measurements and remain therefore poorly constrained. The sectorial partitioning of this emission increase between the periods 2002-2006 and 2008-2012 differs from one atmospheric inversion study to another. However, all top-down studies suggest smaller changes in fossil fuel emissions (from oil, gas, and coal industries) compared to the mean of the bottom-up inventories included in this study. This difference is partly driven by a smaller emission change in China from the top-down studies compared to the estimate in the Emission Database for Global Atmospheric Research (EDGARv4.2) inventory, which should be revised to smaller values in a near future. We apply isotopic signatures to the emission changes estimated for individual studies based on five emission sectors and find that for six individual top-down studies (out of eight) the average isotopic signature of the emission changes is not consistent with the observed change in atmospheric 13 CH 4 . However, the partitioning in emission change derived from the ensemble mean is consistent with this isotopic constraint. At the global scale, the top-down ensemble mean suggests that the dominant contribution to the resumed atmospheric CH 4 growth after 2006 comes from microbial sources (more from agriculture and waste sectors than from natural wetlands), with an uncertain but smaller contribution from fossil CH 4 emissions. In addition, a decrease in biomass burning emissions (in agreement with the biomass burning emission databases) makes the balance of sources consistent with atmospheric 13 CH 4 observations. In most of the top-down studies included here, OH concentrations are considered constant over the years (seasonal variations but without any inter-annual variability). As a result, the methane loss (in particular through OH oxidation) varies mainly through the change in methane concentrations and not its oxidants. For these reasons, changes in the methane loss could not be properly investigated in this study, although it may play a significant role in the recent atmospheric methane changes as briefly discussed at the end of the paper.

58 GEOSCIENCES↗

Aerosols at the poles: an AeroCom Phase II multi-model evaluation

Abstract. Atmospheric aerosols from anthropogenic and natural sources reach the polar regions through long-range transport and affect the local radiation balance. Such transport is, however, poorly constrained in present-day global climate models, and few multi-model evaluations of polar anthropogenic aerosol radiative forcing exist. Here we compare the aerosol optical depth (AOD) at 550 nm from simulations with 16 global aerosol models from the AeroCom Phase II model intercomparison project with available observations at both poles. We show that the annual mean multi-model median is representative of the observations in Arctic, but that the intermodel spread is large. We also document the geographical distribution and seasonal cycle of the AOD for the individual aerosol species: black carbon (BC) from fossil fuel and biomass burning, sulfate, organic aerosols (OAs), dust, and sea-salt. For a subset of models that represent nitrate and secondary organic aerosols (SOAs), we document the role of these aerosols at high latitudes.The seasonal dependence of natural and anthropogenic aerosols differs with natural aerosols peaking in winter (sea-salt) and spring (dust), whereas AOD from anthropogenic aerosols peaks in late spring and summer. The models produce a median annual mean AOD of 0.07 in the Arctic (defined here as north of 60° N). The models also predict a noteworthy aerosol transport to the Antarctic (south of 70° S) with a resulting AOD varying between 0.01 and 0.02. The models have estimated the shortwave anthropogenic radiative forcing contributions to the direct aerosol effect (DAE) associated with BC and OA from fossil fuel and biofuel (FF), sulfate, SOAs, nitrate, and biomass burning from BC and OA emissions combined. The Arctic modelled annual mean DAE is slightly negative (−0.12 W m−2), dominated by a positive BC FF DAE in spring and a negative sulfate DAE in summer. The Antarctic DAE is governed by BC FF. We perform sensitivity experiments with one of the AeroCom models (GISS modelE) to investigate how regional emissions of BC and sulfate and the lifetime of BC influence the Arctic and Antarctic AOD. A doubling of emissions in eastern Asia results in a 33 % increase in Arctic AOD of BC. A doubling of the BC lifetime results in a 39 % increase in Arctic AOD of BC. However, these radical changes still fall within the AeroCom model range.

54 ENVIRONMENTAL SCIENCES↗

Detection and variability of combustion-derived vapor in an urban basin

Abstract. Water emitted during combustion may comprise a significant portion of ambient humidity (&gt; 10 %) in urban areas, where combustion emissions are strongly focused in space and time. Stable water vapor isotopes can be used to apportion measured humidity values between atmospherically transported and combustion-derived water vapor, as combustion-derived vapor possesses an unusually negative deuterium excess value (d-excess, d = δ 2 H − 8 δ 18 O). We investigated the relationship between the d-excess of atmospheric vapor, ambient CO 2 concentrations, and atmospheric stability across four winters in Salt Lake City, Utah. We found a robust inverse relationship between CO 2 excess above background and d-excess on sub-diurnal to seasonal timescales, which was most prominent during periods of strong atmospheric stability that occur during Salt Lake City winter. Using a Keeling-style mixing model approach, and assuming a molar ratio of H 2 O to CO 2 in emissions of 1.5, we estimated the d-excess of combustion-derived vapor in Salt Lake City to be −179 ± 17 ‰, consistent with the upper limit of theoretical estimates. Based on this estimate, we calculate that vapor from fossil fuel combustion often represents 5–10 % of total urban humidity, with a maximum estimate of 16.7 %, consistent with prior estimates for Salt Lake City. Moreover, our analysis highlights that changes in the observed d-excess during periods of high atmospheric stability cannot be explained without a vapor source possessing a strongly negative d-excess value. Further refinements in this humidity apportionment method, most notably empirical validation of the d-excess of combustion vapor or improvements in the estimation of the background d-excess value in the absence of combustion, can yield more certain estimates of the impacts of fossil fuel combustion on urban humidity and meteorology.

54 ENVIRONMENTAL SCIENCES↗

Online characterization of primary and secondary emissions of particulate matter and acidic molecules from a modern fleet of city buses

The potential impact of transitioning from conventional fossil fuel to a non-fossil-fuel vehicle fleet was investigated by measuring primary emissions via extractive sampling of bus plumes and assessing secondary mass formation using the Gothenburg Potential Aerosol Mass (Go:PAM) reactor from 76 in-use transit buses. Online chemical characterization of gaseous and particulate emissions from these buses was conducted using chemical ionization mass spectrometry (CIMS) with acetate as the reagent ion, coupled with the Filter Inlet for Gases and AEROsols (FIGAERO). Acetate reagent ion chemistry selectively ionizes acidic compounds, including organic and inorganic acids, as well as nitrated and sulfated organics. A significant reduction (48 %–98 %) in fresh particle emissions was observed in buses utilizing compressed natural gas (CNG), biodiesels like rapeseed methyl ester (RME) and hydrotreated vegetable oil (HVO), and hybrid-electric HVO (HVO HEV ) compared to diesel (DSL). However, secondary particle formation from photooxidation of emissions was substantial across all the fuel types. The median ratio of particle mass emission factors of aged to fresh emissions increased in the following order: DSL buses at 4.0, HVO buses at 6.7, HVO HEV buses at 10.5, RME buses at 10.8, and CNG buses at 84. Of the compounds that can be identified by CIMS, fresh gaseous emissions from all Euro V/EEV (Enhanced Environmentally friendly Vehicle) buses, regardless of fuel type, were dominated by nitrogen-containing compounds such as nitrous acid (HONO), nitric acid (HNO 3 ), and isocyanic acid (HNCO), alongside small monoacids (C 1 –C 3 ). Notably, the emission of nitrogen-containing compounds was lower in Euro VI buses equipped with more advanced emission control technologies. Secondary gaseous organic acids correlated strongly with gaseous HNO 3 signals (R 2 =0.85–0.99) in Go:PAM, but their moderate to weak correlations with post-photooxidation secondary particle mass suggest that they are not reliable tracers of secondary organic aerosol formation from bus exhaust. Our study highlights that non-regulated compounds and secondary pollutant formation, not currently addressed in legislation, are crucial considerations in the evaluation of environmental impacts of future fuel and engine technology shifts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantification of regional net CO 2 flux errors in the Orbiting Carbon Observatory-2 (OCO-2) v10 model intercomparison project (MIP) ensemble using airborne measurements

Inverse model intercomparison projects (MIPs) provide a chance to assess the uncertainties in inversion estimates arising from various sources. However, accurately quantifying ensemble CO 2 flux errors remains challenging and often relies on the ensemble spread. This study proposes a method for quantifying the errors in regional net surface–atmosphere CO 2 flux estimates from models taken from the Orbiting Carbon Observatory-2 (OCO-2) v10 MIP by using independent airborne CO 2 measurements for the period 2015–2017. We first calculate the root mean square error (RMSE) between the ensemble mean of posterior CO 2 concentrations and airborne observations and then isolate the CO 2 concentration errors caused solely by the ensemble mean of posterior net fluxes by subtracting the observation, representation, and transport errors from seven regions. Our analysis reveals that the flux errors projected onto CO 2 space account for 55 %–85 % of the regional average RMSE over the 3 years, ranging from 0.88 to 1.91 ppm. In five regions, the error estimates based on observations exceed those computed from the ensemble spread of posterior fluxes by a factor of 1.33–1.93, implying an underestimation of the actual flux errors, while their magnitudes are comparable in two regions. The adjoint sensitivity analysis identifies that the underestimation of flux errors is prominent where the magnitudes of fossil fuel emissions exceed those of terrestrial-biosphere fluxes by a factor of 3–31 over the 3 years. This suggests the presence of systematic biases in the inversion estimates associated with errors in the prescribed fossil fuel emissions common to all models. Our study emphasizes the value of airborne measurements for quantifying regional errors in ensemble net CO 2 flux estimates.

54 ENVIRONMENTAL SCIENCES↗

A model of the methane cycle, permafrost, and hydrology of the Siberian continental margin

A two-dimensional model of a sediment column, with Darcy fluid flow, biological and thermal methane production, and permafrost and methane hydrate formation, is subjected to glacial–interglacial cycles in sea level, alternately exposing the continental shelf to the cold atmosphere during glacial times and immersing it in the ocean in interglacial times. The glacial cycles are followed by a "long-tail" 100 kyr warming due to fossil fuel combustion. The salinity of the sediment column in the interior of the shelf can be decreased by hydrological forcing to depths well below sea level when the sediment is exposed to the atmosphere. There is no analogous advective seawater-injecting mechanism upon resubmergence, only slower diffusive mechanisms. This hydrological ratchet is consistent with the existence of freshwater beneath the sea floor on continental shelves around the world, left over from the last glacial period. The salt content of the sediment column affects the relative proportions of the solid and fluid H 2 O-containing phases, but in the permafrost zone the salinity in the pore fluid brine is a function of temperature only, controlled by equilibrium with ice. Ice can tolerate a higher salinity in the pore fluid than methane hydrate can at low pressure and temperature, excluding methane hydrate from thermodynamic stability in the permafrost zone. The implication is that any methane hydrate existing today will be insulated from anthropogenic climate change by hundreds of meters of sediment, resulting in a response time of thousands of years. The strongest impact of the glacial–interglacial cycles on the atmospheric methane flux is due to bubbles dissolving in the ocean when sea level is high. When sea level is low and the sediment surface is exposed to the atmosphere, the atmospheric flux is sensitive to whether permafrost inhibits bubble migration in the model. If it does, the atmospheric flux is highest during the glaciating, sea level regression (soil-freezing) part of the cycle rather than during deglacial transgression (warming and thawing). The atmospheric flux response to a warming climate is small, relative to the rest of the methane sources to the atmosphere in the global budget, because of the ongoing flooding of the continental shelf. The increased methane flux due to ocean warming could be completely counteracted by a sea level rise of tens of meters on millennial timescales due to the loss of ice sheets, decreasing the efficiency of bubble transit through the water column. The model results give no indication of a mechanism by which methane emissions from the Siberian continental shelf could have a significant impact on the near-term evolution of Earth's climate, but on millennial timescales the release of carbon from hydrate and permafrost could contribute significantly to the fossil fuel carbon burden in the atmosphere–ocean–terrestrial carbon cycle.

Archer, D.↗

On a simplified solution of climate-carbon dynamics in idealized flat10MIP simulations

Abstract. Idealized experiments with coupled climate-carbon Earth system models (ESMs) provide a basis for understanding the response of the carbon cycle to external forcing and for quantifying climate-carbon feedbacks. Here, we analyze globally-averaged results from idealized esm-flat10 experiments and show that most models exhibit a quasi-linear relationship between cumulative carbon uptake on land and in the ocean during a period of constant fossil fuel emissions of 10 Pg C yr−1. We hypothesize that this relationship does not depend on emission pathways. Further, as a simplification, we quantify the relationship between cumulative ocean carbon uptake and changes in ocean heat content using a linear approximation. In this way, changes in oceanic heat content and atmospheric CO2 concentration become interdependent variables, reducing the coupled temperature-CO2 system to just one differential equation. The equation can be solved analytically or numerically for the atmospheric CO2 concentration as a function of fossil fuel emissions. This approach leads to a simplified description of global carbon and climate dynamics, which could be used for applications beyond existing analytical frameworks.

Brovkin, Victor↗

Global Carbon Budget 2017

Here an accurate assessment of anthropogenic carbon dioxide (CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the global carbon budget – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. CO 2 emissions from fossil fuels and industry ( E FF ) are based on energy statistics and cement production data, respectively, while emissions from land-use change ( E LUC ), mainly deforestation, are based on land-cover change data and bookkeeping models. The global atmospheric CO 2 concentration is measured directly and its rate of growth ( G ATM ) is computed from the annual changes in concentration. The ocean CO 2 sink ( S OCEAN ) and terrestrial CO 2 sink ( S LAND ) are estimated with global process models constrained by observations. The resulting carbon budget imbalance ( B IM ), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1 σ . For the last decade available (2007–2016), E FF was 9.4 ± 0.5 GtC yr –1 , E LUC 1.3 ± 0.7 GtC yr –1 , G ATM 4.7 ± 0.1 GtC yr –1 , S OCEAN 2.4 ± 0.5 GtC yr –1 , and S LAND 3.0 ± 0.8 GtC yr –1 , with a budget imbalance B IM of 0.6 GtC yr –1 indicating overestimated emissions and/or underestimated sinks. For year 2016 alone, the growth in E FF was approximately zero and emissions remained at 9.9 ±\ 0.5 GtC yr –1 . Also for 2016, E LUC was 1.3 ± .7 GtC yr –1 , G ATM was 6.1 ± 0.2 GtC yr –1 , S OCEAN was 2.6 ± 0.5 GtC yr –1 , and S LAND was 2.7 ± 1.0 GtC yr –1 , with a small B IM of –0.3 GtC. G ATM continued to be higher in 2016 compared to the past decade (2007–2016), reflecting in part the high fossil emissions and the small S LAND consistent with El Niño conditions. The global atmospheric CO 2 concentration reached 402.8 ± 0.1 ppm averaged over 2016. For 2017, preliminary data for the first 6–9 months indicate a renewed growth in E FF of +2.0 % (range of 0.8 to 3.0 %) based on national emissions projections for China, USA, and India, and projections of gross domestic product (GDP) corrected for recent changes in the carbon intensity of the economy for the rest of the world. This living data update documents changes in the methods and data sets used in this new global carbon budget compared with previous publications of this data set.

54 ENVIRONMENTAL SCIENCES↗

Global Carbon Budget 2021

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize datasets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the first time, an approach is shown to reconcile the difference in our ELUC estimate with the one from national greenhouse gas inventories, supporting the assessment of collective countries' climate progress. For the year 2020, EFOS declined by 5.4 % relative to 2019, with fossil emissions at 9.5 ± 0.5 GtC yr−1 (9.3 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 0.9 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission of 10.2 ± 0.8 GtC yr−1 (37.4 ± 2.9 GtCO2). Also, for 2020, GATM was 5.0 ± 0.2 GtC yr−1 (2.4 ± 0.1 ppm yr−1), SOCEAN was 3.0 ± 0.4 GtC yr−1, and SLAND was 2.9 ± 1 GtC yr−1, with a BIM of −0.8 GtC yr−1. The global atmospheric CO2 concentration averaged over 2020 reached 412.45 ± 0.1 ppm. Preliminary data for 2021 suggest a rebound in EFOS relative to 2020 of +4.8 % (4.2 % to 5.4 %) globally. Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2020, but discrepancies of up to 1 GtC yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use changes emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and datasets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this dataset (Friedlingstein et al., 2020, 2019; Le Quéré et al., 2018b, a, 2016, 2015b, a, 2014, 2013). The data presented in this work are available at https://doi.org/10.18160/gcp-2021 (Friedlingstein et al., 2021).

Friedlingstein, Pierre (ORCID:0000000333094739)↗

Global Carbon Budget 2022

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2021, EFOS increased by 5.1 % relative to 2020, with fossil emissions at 10.1 ± 0.5 GtC yr−1 (9.9 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 1.1 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 10.9 ± 0.8 GtC yr−1 (40.0 ± 2.9 GtCO2). Also, for 2021, GATM was 5.2 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEAN was 2.9 ± 0.4 GtC yr−1, and SLAND was 3.5 ± 0.9 GtC yr−1, with a BIM of −0.6 GtC yr−1 (i.e. the total estimated sources were too low or sinks were too high). The global atmospheric CO2 concentration averaged over 2021 reached 414.71 ± 0.1 ppm. Preliminary data for 2022 suggest an increase in EFOS relative to 2021 of +1.0 % (0.1 % to 1.9 %) globally and atmospheric CO2 concentration reaching 417.2 ppm, more than 50 % above pre-industrial levels (around 278 ppm). Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2021, but discrepancies of up to 1 GtC yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extratropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set. The data presented in this work are available at https://doi.org/10.18160/GCP-2022 (Friedlingstein et al., 2022b).

54 ENVIRONMENTAL SCIENCES↗

National CO 2 budgets (2015–2020) inferred from atmospheric CO 2 observations in support of the global stocktake

Accurate accounting of emissions and removals of CO 2 is critical for the planning and verification of emission reduction targets in support of the Paris Agreement. Here, we present a pilot dataset of country-specific net carbon exchange (NCE; fossil plus terrestrial ecosystem fluxes) and terrestrial carbon stock changes aimed at informing countries' carbon budgets. These estimates are based on “top-down” NCE outputs from the v10 Orbiting Carbon Observatory (OCO-2) modeling intercomparison project (MIP), wherein an ensemble of inverse modeling groups conducted standardized experiments assimilating OCO-2 column-averaged dry-air mole fraction (X CO 2 ) retrievals (ACOS v10), in situ CO 2 measurements or combinations of these data. The v10 OCO-2 MIP NCE estimates are combined with “bottom-up” estimates of fossil fuel emissions and lateral carbon fluxes to estimate changes in terrestrial carbon stocks, which are impacted by anthropogenic and natural drivers. These flux and stock change estimates are reported annually (2015–2020) as both a global 1° x 1° gridded dataset and a country-level dataset and are available for download from the Committee on Earth Observation Satellites' (CEOS) website. Across the v10 OCO-2 MIP experiments, we obtain increases in the ensemble median terrestrial carbon stocks of 3.29–4.58 Pg CO2 yr -1 (0.90–1.25 Pg C yr -1 ). This is a result of broad increases in terrestrial carbon stocks across the northern extratropics, while the tropics generally have stock losses but with considerable regional variability and differences between v10 OCO-2 MIP experiments. We discuss the state of the science for tracking emissions and removals using top-down methods, including current limitations and future developments towards top-down monitoring and verification systems.

54 ENVIRONMENTAL SCIENCES↗