Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “environmental fate”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High methane concentrations in tidal salt marsh soils: Where does the methane go?

Abstract Tidal salt marshes produce and emit CH 4 . Therefore, it is critical to understand the biogeochemical controls that regulate CH 4 spatial and temporal dynamics in wetlands. The prevailing paradigm assumes that acetoclastic methanogenesis is the dominant pathway for CH 4 production, and higher salinity concentrations inhibit CH 4 production in salt marshes. Recent evidence shows that CH 4 is produced within salt marshes via methylotrophic methanogenesis, a process not inhibited by sulfate reduction. To further explore this conundrum, we performed measurements of soil–atmosphere CH 4 and CO 2 fluxes coupled with depth profiles of soil CH 4 and CO 2 pore water gas concentrations, stable and radioisotopes, pore water chemistry, and microbial community composition to assess CH 4 production and fate within a temperate tidal salt marsh. We found unexpectedly high CH 4 concentrations up to 145,000 μmol mol −1 positively correlated with S 2− (salinity range: 6.6–14.5 ppt). Despite large CH 4 production within the soil, soil–atmosphere CH 4 fluxes were low but with higher emissions and extreme variability during plant senescence (84.3 ± 684.4 nmol m −2 s −1 ). CH 4 and CO 2 within the soil pore water were produced from young carbon, with most Δ 14 C‐CH 4 and Δ 14 C‐CO 2 values at or above modern. We found evidence that CH 4 within soils was produced by methylotrophic and hydrogenotrophic methanogenesis. Several pathways exist after CH 4 is produced, including diffusion into the atmosphere, CH 4 oxidation, and lateral export to adjacent tidal creeks; the latter being the most likely dominant flux. Our findings demonstrate that CH 4 production and fluxes are biogeochemically heterogeneous, with multiple processes and pathways that can co‐occur and vary in importance over the year. This study highlights the potential for high CH 4 production, the need to understand the underlying biogeochemical controls, and the challenges of evaluating CH 4 budgets and blue carbon in salt marshes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A Comprehensive Greenhouse Gas Assessment of Biomass-Based Carbon Dioxide Removal and Storage

Biomass-based Carbon Removal and Storage (BiCRS) is a suite of technological pathways that provide a relatively affordable and high-potential way to remove carbon dioxide from the atmosphere and simultaneously re-purpose residual biomass. BiCRS technologies are individually optimized to different biomass feedstocks, temperatures and oxygen, with the primary goal of capturing a high percentage of biomass-derived carbon dioxide and storing it in permanent below-ground storage (i.e. geologic storage), and a secondary goal of producing renewable energy. The current standard for BiCRS life cycle assessments does not account for greenhouse gas (carbon dioxide, methane and nitrous oxide) fluxes from surface soil amendment of BiCRS byproducts (e.g. char, ash) and their alternate fates. This is partly due to a paucity of empirical data due to the novelty of BiCRS conversion technologies, the range in feedstocks, and the heterogeneity of soils to which they might be amended. To fill this gap in in the life-cycle assessment and provide a realistic range of parameters for a soil amendment component of the BiCRS life cycle assessment, we conducted an incubation experiment to measure soil carbon changes, microbial respiration, methane fluxes, and nitrous oxide fluxes from two different soils amended with char- and ash- byproducts from biomass after gasification to hydrogen, fast pyrolysis to bio-oil, and torrefaction. Our results suggest that nitrous oxide and methane emission or consumption from BiCRS char amendments to soil are context dependent. Nitrous oxide emissions from amendments are higher in soils with higher pH, whereas some char amendments to low pH soil can reduce N 2 O emissions relative to control. In contrast, methane is emitted from BiCRS char amended to low pH soil, but consumed or neutral relative to control when amended to high pH soil. We present carbon and nitrogen mass balance throughout the experiment to help understand short-term durability of BiCRS chars and their counterfactuals, and preliminary suggestions for revisions to the broader BiCRS Measurement, Reporting, and Verification framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Modeling framework for flocculated cohesive sediment transport in the current bottom boundary layer

Cohesive sediment transport, where its settling velocity is controlled by the flocculation process, is a crucial component in determining biochemical cycles, fate of pollutants, and morphodynamics in many aquatic ecosystems. In this study, a modeling framework is presented to investigate how flocculation influences cohesive sediment transport in the current bottom boundary layer in dilute conditions, consistent with the calibration range of the flocculation model. From a local analysis of floc dynamics in homogenous turbulence, we identify that the floc size distribution is mainly controlled by floc cohesion and yield strength. The uncertainty in fractal dimension plays a minor role for the floc size but it influences the resulting floc density and settling velocity. The transport analysis in the current boundary layer shows that the flocculation process alters the vertical distribution of the settling velocity and hence the sediment concentration with a strong dependence on cohesion, floc yield strength, and floc structure. When the flocs are more susceptible to breaking, a well-mixed concentration profile is obtained. In contrast, for flocs with higher cohesion or yield strength, higher concentration with a sharp gradient is observed close to the bed. Overall, the settling velocity exhibits a low vertical variability within 20 % of the depth-averaged value except near the bed. Further, this suggests that using a depth-averaged settling velocity yields acceptable predictions of the sediment concentration profiles, especially for flocs with lower cohesion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Coordination of rooting, xylem, and stomatal strategies explains the response of conifer forest stands to multi-year drought in the southern Sierra Nevada of California

Abstract. Extreme droughts are a major determinant of ecosystem disturbance that impacts plant communities and feeds back into climate change through changes in plant functioning. However, the complex relationships between aboveground and belowground plant hydraulic traits and their role in governing plant responses to drought are not fully understood. In this study, we use a model, the Functionally Assembled Terrestrial Ecosystem Simulator in a configuration that includes plant hydraulics (FATES-Hydro), to investigate ecosystem responses to the 2012–2015 California drought in comparison with observations at a site in the southern Sierra Nevada that experienced widespread tree mortality during this drought. We conduct a sensitivity analysis to explore how different plant water sourcing and hydraulic strategies lead to differential responses during normal and drought conditions. The analysis shows the following. Deep roots that sustain productivity through the dry season are needed for the model to capture observed seasonal cycles of evapotranspiration (ET) and gross primary productivity (GPP) in normal years, and deep-rooted strategies are nonetheless subject to large reductions in ET and GPP when the deep soil reservoir is depleted during extreme droughts, in agreement with observations. Risky stomatal strategies lead to greater productivity during normal years as compared to safer stomatal control, but they also lead to a high risk of xylem embolism during the 2012–2015 drought. For a given stand density, stomatal and xylem traits have a stronger impact on plant water status than on ecosystem-level fluxes. Our study highlights the significance of resolving plant water sourcing strategies to represent drought impacts on plants and consequent feedbacks in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Ambient Aerosol Is Physically Larger on Cloudy Days in Bondville, Illinois

Particle chemical composition affects aerosol optical and physical properties in ways important for the fate, transport, and impact of atmospheric particulate matter. For example, hygroscopic constituents take up water to increase the physical size of a particle, which can alter the extinction properties and atmospheric lifetime. At the collocated AERosol RObotic NETwork (AERONET) and Interagency Monitoring of PROtected Visual Environments (IMPROVE) network monitoring stations in rural Bondville, Illinois, we employ a novel cloudiness determination method to compare measured aerosol physicochemical properties on predominantly cloudy and clear sky days from 2010 to 2019. On cloudy days, aerosol optical depth (AOD) is significantly higher than on clear sky days in all seasons. Measured angstrom ngstro''m exponents are significantly smaller on cloudy days, indicating physically larger average particle size for the sampled populations in all seasons except winter. Mass concentrations of fine particulate matter that include estimates of aerosol liquid water (ALW) are higher on cloudy days in all seasons but winter. More ALW on cloudy days is consistent with larger particle sizes inferred from angstrom ngstro''m exponent measurements. Aerosol chemical composition that affects hygroscopicity plays a determining impact on cloudy versus clear sky differences in AOD, angstrom ngstro''m exponents, and ALW. This work highlights the need for simultaneous collocated, high-time-resolution measurements of both aerosol chemical and physical properties, in particular at cloudy times when quantitative understanding of tropospheric composition is most uncertain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Facet-dependent electron transfer induces distinct arsenic reallocations on hematite

The interfacial electron transfer (ET) between electron shuttling compounds and iron (Fe) oxyhydroxides plays a crucial role in the reductive dissolution of Fe minerals and the fate of surface-bound arsenic (As). However, the impact of exposed facets of highly crystalline hematite on reductive dissolution and As immobilization is poorly understood. In this study, we systematically investigated the interfacial processes of the electron shuttling compound cysteine (Cys) on various facets of hematite and the reallocations of surface-bound As(III) or As(V) on the respective surfaces. Our results demonstrate that the ET process between Cys and hematite generates Fe(II) and leads to reductive dissolution, with more Fe(II) generated on {001} facets of exposed hematite nanoplates (HNPs). Reductive dissolution of hematite leads to significantly enhanced As(V) reallocations on hematite. Nevertheless, upon the addition of Cys, a raipd release of As(III) can be halted by its prompt re-adsorption, leaving the extent of As(III) immobilization on hematite unchanged throughout the course of reductive dissolution. This is due to that Fe(II) can form new precipitates with As(V), a process that is facet-dependent and influenced by water chemistry. Electrochemical analysis reveals that HNPs exhibit higher conductivity and ET ability, which is beneficial for reductive dissolution and As reallocations on hematite. In conclusion, these findings highlight the facet-dependent reallocations of As(III) and As(V) facilitated by electron shuttling compounds and have implications for the biogeochemical processes of As in soil and subsurface environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Modeling the impact of extreme weather events and future climate on the radiologically contaminated sites of Enewetak Atoll

Enewetak Atoll underwent 43 historical nuclear tests from 1948 to 1958, including the first hydrogen bomb test, resulting in a substantial nuclear material fallout contaminating the Atoll and the lagoon waters. The radionuclide fallout material deposited in lagoon sediments and land soil will remain for decades to come. With intensifying climate and extreme weather events, the possibility of redistribution of deposited radionuclide material has become a great concern. This study uses a numerical modeling approach to estimate the potential elevated radionuclide concentrations that can be distributed during storm events under current and future climates. We simulated three historical storm scenarios that are most likely to impact Atoll’s environment and remobilize the radionuclide-bound sediments. WRF-ARW was used to reconstruct these storm scenarios under current year (2015) and future year (2090) climates. Storm-induced ocean hydrodynamics conditions were generated using FVCOM. FVCOM-ICM was externally coupled to simulate the fate and transport of radionuclides. Given that the 239 Pu is the largest inventory of the lagoon and Atoll islands, the model results show the highest average 239 Pu concentration that an island may be exposed to is 3.25E-4 Bq/m 3 (becquerel per cubic meters), which is an increase of 84 times the average baseline/existing 239 Pu concentration without the storm conditions. The overall increase in 239 Pu average over all the islands of Atoll is about 20 folds relative to the baseline concentration. Despite the high relative increase ratios, the significantly low activity concentrations may not pose an immediate exposure risk. However, due to the limitations of the study and uncertainties/biases in the historical data used, further research supported by field surveys to better characterize the current contamination level may be needed to make more accurate predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Resolving dynamic mineral-organic interactions in the rhizosphere by combining in-situ microsensors with plant-soil reactive transport modeling

Associations between minerals and organic matter represent one of the most important carbon storage mechanisms in soils. Plant roots are major sources of soil carbon, and resolving the dynamics and dominance of microbial consumption versus mineral sorption of root-derived carbon is critical to understanding soil carbon storage. Here we integrate in-situ rhizosphere microsensor and plant physiological measurements with a 3-D plant-soil reactive transport model to explore the fate of dissolved organic carbon (DOC) in the rhizosphere, particularly its microbial consumption and interaction with Fe oxide minerals. Over several days, a microdialysis probe sampling pore water at the root-soil interface of growing Vicia faba roots in live soil, revealed clear diel patterns of DOC concentration. Daytime DOC spikes coincided with peaks in leaf-level photosynthesis rates and were accompanied by declining redox potential and dissolved oxygen as well as increasing pH in the rhizosphere. Incorporating microsensor data into our modeling framework showed that the measured rapid loss of DOC after each mid-day spike could not be explained by consumption via aerobic respiration, nor via anaerobic respiration dominated by Fe oxide reduction. Rather, in the model, a large fraction of rhizosphere DOC was rapidly immobilized each day by adsorption to Fe oxides. Further, modeled microbial Fe reduction (fueled by DOC) did not mobilize significant organic carbon from Fe oxides during the day. Instead, the model predicted equilibrium desorption of organic carbon from Fe oxides at night. This new mechanistic modeling framework, which couples aboveground plant physiological measurements with non-destructive high-resolution monitoring of rhizosphere processes, has great potential for exploring the dynamics and balance of the various microbial reactions and mineral interactions controlling carbon transformations and storage in soils.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Dinitrogen Emissions Dominate Nitrogen Gas Emissions From Soils With Low Oxygen Availability in a Moist Tropical Forest

Lowland tropical forest soils are relatively N rich and are the largest global source of N 2 O (a powerful greenhouse gas) to the atmosphere. Despite the importance of tropical N cycling, there have been few direct measurements of N 2 (an inert gas that can serve as an alternate fate for N 2 O) in tropical soils, limiting our ability to characterize N budgets, manage soils to reduce N 2 O production, or predict the future role that N limitation to primary productivity will play in buffering against climate change. We collected soils from across macro‐ and micro‐topographic gradients that have previously been shown to differ in O 2 availability and trace gas emissions. We then incubated these soils under oxic and anoxic headspaces to explore the relative effect of soil location versus transient redox conditions. No matter where the soils came from, or what headspace O 2 was used in the incubation, N 2 emissions dominated the flux of N gas losses. In the macrotopography plots, production of N 2 and N 2 O were higher in low O 2 valleys than on more aerated ridges and slopes. In the microtopography plots, N 2 emissions from plots with lower mean soil O 2 (5%–10%) were greater than in plots with higher mean soil O 2 (10%–20%). We estimate an N gas flux of ∼37 kg N/ha/yr from this forest, 99% as N 2 . These results suggest that N 2 fluxes may have been systematically underestimated in these landscapes, and that the measurements we present call for a reevaluation of the N budgets in lowland tropical forest ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A robust method of nuclei isolation for single-cell RNA sequencing of solid tissues from the plant genus Populus

Single-cell transcriptome analysis has been extensively applied in humans and animal models to uncover gene expression heterogeneity between the different cell types of a tissue or an organ. It demonstrated its capability to discover key regulatory elements that determine cell fate during developmental programs. Single-cell analysis requires the isolation and labeling of the messenger RNA (mRNA) derived from each cell. These challenges were primarily addressed in mammals by developing microfluidic-based approaches. For plant species whose cells contain cell walls, these approaches have generally required the generation of isolated protoplasts. Many plant tissues’ secondary cell wall hinders enzymatic digestion required for individual protoplast isolation, resulting in an unequal representation of cell types in a protoplast population. This limitation is especially critical for cell types located in the inner layers of a tissue or the inner tissues of an organ. Consequently, single-cell RNA sequencing (scRNA-seq) studies using microfluidic approaches in plants have mainly been restricted to Arabidopsis roots, for which well-established procedures of protoplast isolation are available. Here we present a simple alternative approach to generating high-quality protoplasts from plant tissue by characterizing the mRNA extracted from individual nuclei instead of whole cells. We developed the protocol using two different plant materials with varying cellular complexity levels and cell wall structure, Populus shoot apices, and more lignified stems. Using the 10× Genomics Chromium technology, we show that this procedure results in intact mRNA isolation and limited leakage, with a broad representation of individual cell transcriptomes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Distribution of Infiltration in the 216-U-10 and 216-B-3 Pond Systems 1944-1997

A modeling system to simulate fate and transport of contaminant plumes in groundwater beneath the Central Plateau to support the Cumulative Impact Evaluation (CIE) and Composite Analysis (CA) is being developed. The CIE will evaluate the cumulative effect of multiple Central Plateau radiological and chemical sources on groundwater quality (DOE/RL-2018-69, Cumulative Impact Evaluation Technical Approach Document ). This will enable cleanup decisions to be evaluated in the context of all Central Plateau waste sites and existing groundwater conditions. The CA evaluates the cumulative effect of all sources of radioactive materials on radiological dose to future members of the public (DOE/RL-2018-60, Annual Status Report (FY 2018): Composite Analysis of Low Level Waste Disposal in the Central Plateau of the Hanford Site ). Fate and transport modeling will be used to forecast future groundwater concentrations for both the CIE and CA. The CIE/CA subsurface modeling system is divided into source (inventory), vadose zone, and saturated zone facets. Source inventory is a necessary input to the vadose zone facet, and releases from the vadose zone to groundwater are inputs to the saturated zone facet. The vadose zone facet simulates historical and future releases from the vadose zone to groundwater for both liquid waste and solid waste disposal sites. Effluent release volumes and radionuclide inventory estimates for liquid waste sites are contained in ECFHANFORD- 17-0079, Hanford Soil Inventory Model (SIM-v2) Calculated Radionuclide Inventory of Direct Liquid Discharges to Soil in the Hanford Site’s 200 Areas , herein referred to as SIM-v2. The liquid waste sites in SIM-v2 include the 216-B-3 Pond near 200 East Area and the 216-U-10 Pond in 200 West Area. The 216-B-3 Pond system consisted of a main pond and 3 expansion lobes that were operational during overlapping time periods. It also included influent ditches, but these are not listed as liquid waste sources in SIM-v2 and are not included in the CIE/CA modeling. The total effluent volume released to the pond system was measured, but the partitioning of this volume to the various expansion lobes was not determined and is not accounted for in SIM-v2. Regarding the 216-U-10 Pond, it received water from several influent ditches. The pond and the influent ditches are listed as liquid waste sites in SIM-v2 and are included as sources in the CIE/CA modeling. Because water in the ditches flowed into the pond, only a portion of the effluent volumes and contaminant inventories assigned to the ditches in SIM-v2 infiltrated from the ditches. Thus, partitioning of infiltration and contaminant inventories between the ditches and the pond is not fully accounted for in SIM-v2. Further, some of the water in the 216-U-10 Pond entered overflow ditches and the amount of this overflow was not recorded and is not accounted for in SIM-v2. Thus, the purpose of this Environmental Calculation File (ECF) is to partition the infiltration of wastewater and contaminant inventory between the main pond and each expansion lobe of the 216-B-3 Pond system, and to partition the infiltration of wastewater and contaminant inventory between the ditches and pond of the 216-U-10 Pond system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Nutrient Dynamics in a Coupled Terrestrial Biosphere and Land Model (ELM-FATES-CNP)

We present a representation of nitrogen and phosphorus cycling in the Functionally Assembled Terrestrial Ecosystem Simulator, a demographic vegetation model within the Energy Exascale Earth System land model. This representation is modular, and designed to allow testing of multiple hypothetical approaches for carbon-nutrient coupling in plants. Novel model hypotheses introduced in this work include, (a) the controls on plant acquisition of aqueous mineralized nutrients in the soil and (b) fairly straight forward methods of allocating nutrients to specific plant organs and their losses through live plant turnover as well as litter fluxes generated through plant mortality. This combines the new with pre-existing hypotheses (such as nitrogen fixation and soil decomposition) into a system that can accommodate plant-soil dynamics for a large number of size- and functional-type-resolved plant cohorts within a time-since-disturbance-resolved ecosystem. Root uptake of nutrients is governed by fine root biomass, and plants vary in their fine root biomass allocation in order to balance carbon and nutrient limitations to growth. We test the sensitivity of the model to a wide range of parameter variations and structural representations, and in the context of observations at Barro Colorado Island, Panama. A key model prediction is that plants in the high-light-availability canopy positions allocate more carbon to fine roots than plants in low-light understory environments, given the widely different carbon versus nutrient constraints of these two niches within a given ecosystem. This model provides a basis for exploring carbon-nutrient coupling with vegetation demography within Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Vadose zone flushing of fertilizer tracked by isotopes of water and nitrate

Abstract A substantial fraction of nitrogen (N) fertilizer applied in agricultural systems is not incorporated into crops and moves below the rooting zone as nitrate (NO 3 − ). Understanding mechanisms for soil N retention below the rooting zone and leaching to groundwater is essential for our ability to track the fate of added N. We used dual stable isotopes of nitrate (δ 15 N–NO 3 − and δ 18 O–NO 3 − ) and water (δ 18 O–H 2 O and δ 2 H–H 2 O) to understand the mechanisms driving nitrate leaching at three depths (0.8, 1.5, and 3.0 m) of an irrigated corn field sampled every 2 weeks from 2016 to 2020 in the southern Willamette Valley, Oregon, USA. Distinct periods of high nitrate concentrations with lower δ 15 N–NO 3 − values indicated that a portion of that nitrate was from recent fertilizer applications. We used a mixing model to quantify nitrate fluxes associated with recently added fertilizer N versus older, legacy soil N during these “fertilizer signal periods.” Nitrate leached below 3.0 m in these periods made up a larger proportion of the total N leached at that depth (∼52%) versus the two shallower depths (∼13%–16%), indicating preferential movement of recently applied fertilizer N through the deep soil into groundwater. Further, N associated with recent fertilizer additions leached more easily when compared to remobilized legacy N. A high volume of fall and winter precipitation may push residual fertilizer N to depth, potentially posing a larger threat to groundwater than legacy N. Optimizing fertilizer N additions could minimize fertilizer losses and reduce nitrate leaching to groundwater.

Agriculture↗
DOE OSTI2022

Degradation and fate of 2,4-dinitroanisole (DNAN) and its intermediates treated with Mg/Cu bimetal: Surface examination with XAS, DFT, and LDI-MS

In this article, a novel Mg-based bimetal reagent (Mg/Cu) was used as an enhanced reductive system to degrade insensitive munition 2,4-dinitroanisole (DNAN), a contaminant found in energetic-laden waste. Degradation of DNAN was significantly impacted by dissolved oxygen and studied in anoxic and oxic bimetal systems (i.e., purging with N 2 , air, or O 2 gas). Degradation occurred through sequential nitroreduction: first one nitro group was reduced (ortho or para) to form short-lived intermediates 2-amino-4-nitroanisole or 4-amino-2-nitroanisole (2-ANAN or 4-ANAN), and then subsequent reduction of the other nitro group formed 2,4-diaminoanisole (DAAN). The nitro-amino intermediates demonstrated regioselective reduction in the ortho position to 2-ANAN; Regioselectivity was also impacted by the anoxic/oxic environment. Under O 2 -purging DNAN degradation rate was slightly enhanced, but most notably O 2 significantly accelerated DAAN generation. DAAN also further degraded only in the oxygenated Mg/Cu system. Adsorption of DNAN byproducts to the reagent occurred regardless of anoxic/oxic condition, resulting in a partition of carbon mass between the adsorbed phase (27%-35%) and dissolved phase (59%-72%). Additional surface techniques were applied to investigate contaminant interaction with Cu. Density functional theory (DFT) calculations identified preferential adsorption structures for DNAN on Cu with binding through two O atoms of one or both nitro groups. X-ray absorption spectroscopy (XAS) measurements determined the oxidation state of catalytic metal Cu and formation of a Cu-O-N bond during treatment. Laser desorption ionization mass spectrometry (LDI-MS) measurements also identified intermediate 2-ANAN adsorbed to the bimetal surface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Matrix Approach to Land Carbon Cycle Modeling

Land ecosystems contribute to climate change mitigation by taking up approximately 30% of anthropogenically emitted carbon. However, estimates of the amount and distribution of carbon uptake across the world's ecosystems or biomes display great uncertainty. The latter hinders a full understanding of the mechanisms and drivers of land carbon uptake, and predictions of the future fate of the land carbon sink. The latter is needed as evidence to inform climate mitigation strategies such as afforestation schemes. To advance land carbon cycle modeling, we have developed a matrix approach. Land carbon cycle models use carbon balance equations to represent carbon exchanges among pools. Our approach organizes this set of equations into a single matrix equation without altering any processes of the original model. The matrix equation enables the development of a theoretical framework for understanding the general, transient behavior of the land carbon cycle. While carbon input and residence time are used to quantify carbon storage capacity at steady state, a third quantity, carbon storage potential, integrates fluxes with time to define dynamic disequilibrium of the carbon cycle under global change. The matrix approach can help address critical contemporary issues in modeling, including pinpointing sources of model uncertainty and accelerating spin-up of land carbon cycle models by tens of times. The accelerated spin-up liberates models from the computational burden that hinders comprehensive parameter sensitivity analysis and assimilation of observational data to improve model accuracy. Such computational efficiency offered by the matrix approach enables substantial improvement of model predictions using ever-increasing data availability. Overall, the matrix approach offers a step change forward for understanding and modeling the land carbon cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Simulating the Transport and Rupture of Pollen in the Atmosphere

Pollen, one type of primary biological aerosol particle (PBAP), is emitted from the terrestrial biosphere and can undergo physical changes in the atmosphere via particle rupture. To examine the fate of pollen and its atmospheric processing, a pollen emission and transport scheme is coupled to the Weather Research and Forecasting Model with Chemistry (WRF-Chem). We simulate the emission of pollen and its impacts on the cloud properties and precipitation in the Southern Great Plains from 12 to 19 April 2013, a period with both high pollen emissions and convective activity. We conduct a suite of ensemble runs that simulate primary pollen and three different pollen rupture mechanisms that generate subpollen particles, including (a) high humidity-induced surface rupture, (b) high humidity-induced in-atmosphere plus surface rupture, and (c) lightning-induced rupture, where in-cloud and cloud-to-ground lightning strikes trigger pollen rupture events. When relative humidity is high (>80%), coarse primary pollen (~1 μg m –3 ) is converted into fine subpollen particles (~1.2e –4 μg m –3 ), which produces 80% more subpollen particles than lightning-induced rupture. The in-atmosphere humidity-driven rupture predominantly produces subpollen particles, which is further enhanced during a frontal thunderstorm. During strong convection, vertical updrafts lift primary pollen and subpollen particles (~0.5e –4 μg m –3 ) to the upper troposphere (~12 km) and laterally transports the ruptured pollen in the anvil top outflow. In regions of high pollen and strong convection, ruptured pollen can influence warm cloud formation by decreasing low cloud (<4 km) cloud water mixing ratios and increasing ice phase hydrometeors aloft (>10 km)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Wildfire‐Induced Losses of Soil Particulate and Mineral‐Associated Organic Carbon Persist for Over 4 Years in a Chaparral Ecosystem

ABSTRACT Wildfires can lower soil carbon (C) stocks directly through combustion, but also indirectly during post‐fire recovery if microbial C demands outpace photosynthetic C inputs. However, how much C is respired by soil microorganisms post‐fire may depend on wildfire effects on particulate organic carbon (POC; mostly plant material accessible to microbes) and/or mineral‐associated organic carbon (MAOC; considered C protected by minerals from decomposers), meaning assessment of wildfire impacts on these pools is necessary to predict microbial decomposition rates and, thus, the fate of soil C. Here, we measured POC, MAOC, pyrogenic organic matter C, plant cover, extracellular enzyme activity (EEA), and microbial community abundance and composition 17 days, and 1, 3, and 4 years after the Holy Fire burned 94 km 2 of fire‐adapted chaparral. The wildfire immediately decreased POC by 50% (from 51 ± 21 to 26 ± 6 g C kg −1 ) and MAOC by 33% (from 9.3 ± 0.9 to 6.3 ± 0.9 g C kg −1 ), consistent with MAOC being less vulnerable to loss than POC. POC decreased by another 38% 1 year post‐fire, consistent with increases in microbial abundance and EEA suggesting increased microbial decomposition. Between 1 and 4 years after the fire, cover of the dominant shrub (Arctostaphylos glandulosa) increased from 3.9% ± 1.6% to 16% ± 5.4% (compared to 58% ± 4.6% in unburned plots), marking the end of net soil C losses. Still, soil C did not increase between 1 and 4 years post‐fire, suggesting plant C inputs did not outpace microbial respiration, a finding consistent with isotopically heavier C from microorganisms raising bulk soil δ 13 C values. As global changes favor increases in wildfire frequency and severity, C losses via combustion and decomposition may outpace plant C inputs during the first 4 years post‐fire in chaparral, slowing the replenishment of soil C stocks.

Biodiversity & Conservation↗
DOE OSTI2023

Quantifying airborne fraction trends and the destination of anthropogenic CO 2 by tracking carbon flows in a simple climate model

Abstract Atmospheric carbon dioxide (CO 2 ) concentrations have increased as a direct result of human activity and are at their highest level over the last 2 million years, with profound impacts on the Earth system. However, the magnitude and future dynamics of land and ocean carbon sinks are not well understood; therefore, the amount of anthropogenic fossil fuel emissions that remain in the atmosphere (the airborne fraction) is poorly constrained. This work aims to quantify the sources and controls of atmospheric CO 2 , the fate of anthropogenic CO 2 over time, and the likelihood of a trend in the airborne fraction. We use Hector v3.0, a coupled simple climate and carbon cycle model with the novel ability to explicitly track carbon as it flows through the Earth system. We use key model parameters in a Monte Carlo analysis of 15 000 model runs from 1750 to 2300. Results are filtered for physical realism against historical observations and CMIP6 projection data, and we calculate the relative importance of parameters controlling how much anthropogenic carbon ends up in the atmosphere. Modeled airborne fraction was roughly 52%, consistent with observational studies. The overwhelming majority of model runs exhibited a negative trend in the airborne fraction from 1960–2020, implying that current-day land and ocean sinks are proportionally taking up more carbon than the atmosphere. However, the percentage of atmospheric CO 2 derived from anthropogenic origins can be much higher because of Earth system feedbacks. We find it peaks at over 90% between 2010–2050. Moreover, when looking at the destination of anthropogenic fossil fuel emissions, only a quarter ends up in the atmosphere while more than half of emissions are taken up by the land sink on centennial timescales. This study evaluates the likelihood of airborne fraction trends and provides insights into the dynamics of anthropogenic CO 2 in the Earth system.

54 ENVIRONMENTAL SCIENCES↗