Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ecosystem model”

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,063 records · Page 59

Contrasting Regional Carbon Cycle Responses to Seasonal Climate Anomalies Across the East-West Divide of Temperate North America

Across temperate North America, interannual variability (IAV) in gross primary production (GPP) and net ecosystem exchange (NEE) and their relationship with environmental drivers are poorly understood. Here, we examine IAV in GPP and NEE and their relationship to environmental drivers using two state-of-the-science flux products: NEE constrained by surface and space-based atmospheric CO 2 measurements over 2010-2015 and satellite up-scaled GPP from FluxSat over 2001-2017. We show that the arid western half of temperate North America provides a larger contribution to IAV in GPP (104% of east) and NEE (127% of east) than the eastern half, in spite of smaller magnitude of annual mean GPP and NEE. This occurs because anomalies in western ecosystems are temporally coherent across the growing season leading to an amplification of GPP and NEE. In contrast, IAV in GPP and NEE in eastern ecosystems is dominated by seasonal compensation effects, associated with opposite responses to temperature anomalies in spring and summer. Terrestrial biosphere models in the MsTMIP ensemble generally capture these differences between eastern and western temperate North America, although there is considerable spread between models.

54 ENVIRONMENTAL SCIENCES↗

Short-term variation in leaf-level water use efficiency in a tropical forest

Here we report the representation of stomatal regulation of transpiration and CO 2 assimilation is key to forecasting terrestrial ecosystem responses to global change. Given its importance in determining the relationship between forest productivity and climate, accurate and mechanistic model representation of the relationship between stomatal conductance (g s ) and assimilation is crucial. We assess possible physiological and mechanistic controls on the estimation of the g 1 (stomatal slope, inversely proportional to water use efficiency) and g 0 (stomatal intercept) parameters, using diurnal gas exchange surveys and leaf level response curves of six tropical broadleaf evergreen tree species. g 1 estimated from ex-situ response curves averaged 50% less than g 1 estimated from survey data. While g 0 and g 1 varied between leaves of different phenological stages, the trend was not consistent among species. We identified a diurnal trend associated with g 1 and g 0 that significantly improved model projections of diurnal trends in transpiration. The accuracy of modelled g s can be improved by accounting for variation in stomatal behavior across diurnal periods, and between measurement approaches, rather than focusing on phenological variation in stomatal behavior. Additional investigation into the primary mechanisms responsible for diurnal variation in g 1 will be required to account for this phenomena in land surface models.

54 ENVIRONMENTAL SCIENCES↗

A Novel Freshwater to Marine Evolutionary Transition Revealed within Methylophilaceae Bacteria from the Arctic Ocean

Bacteria inhabiting polar oceans, particularly the Arctic Ocean, are less studied than those at lower latitudes. Discovering bacterial adaptations to Arctic Ocean conditions is essential for understanding responses to the accelerated environmental changes occurring in the North. The Methylophilaceae are emerging as a model for investigating the genomic basis of habitat adaptation, because related lineages are widely distributed across both freshwater and marine ecosystems. Here, we investigated Methylophilaceae diversity in the salinity-stratified surface waters of the Canada Basin, Arctic Ocean. In addition to a diversity of marine OM43 lineages, we report on the genomic characteristics and evolution of a previously undescribed Methylophilaceae clade (BS01) common to polar surface waters yet related to freshwater sediment Methylotenera species. BS01 is restricted to the lower-salinity surface waters, while OM43 is found throughout the halocline. An acidic proteome supports a marine lifestyle for BS01, but gene content shows increased metabolic versatility compared to OM43 and evidence for ongoing genome-streamlining. Phylogenetic reconstruction shows that BS01 colonized the pelagic ocean independently of OM43 via convergent evolution. Salinity adaptation and differences in one-carbon and nitrogen metabolism may play a role in niche differentiation between BS01 and OM43. In particular, urea utilization by BS01 is predicted to provide an ecological advantage over OM43 given the limited amount of inorganic nitrogen in the Canada Basin. These observations provide further evidence that the Arctic Ocean is inhabited by distinct bacterial groups and that at least one group (BS01) evolved via a freshwater to marine environmental transition.

59 BASIC BIOLOGICAL SCIENCES↗

Scaling carbon fluxes from eddy covariance sites to globe: synthesis and evaluation of the FLUXCOM approach

FLUXNET comprises globally distributed eddy-covariance-based estimates of carbon fluxes between the biosphere and the atmosphere. Since eddy covariance flux towers have a relatively small footprint and are distributed unevenly across the world, upscaling the observations is necessary to obtain global-scale estimates of biosphere–atmosphere exchange. Based on cross-consistency checks with atmospheric inversions, sun-induced fluorescence (SIF) and dynamic global vegetation models (DGVMs), here we provide a systematic assessment of the latest upscaling efforts for gross primary production (GPP) and net ecosystem exchange (NEE) of the FLUXCOM initiative, where different machine learning methods, forcing data sets and sets of predictor variables were employed. Spatial patterns of mean GPP are consistent across FLUXCOM and DGVM ensembles (R 2 >0.94 at 1° spatial resolution) while the majority of DGVMs show, for 70 % of the land surface, values outside the FLUXCOM range. Global mean GPP magnitudes for 2008–2010 from FLUXCOM members vary within 106 and 130 PgC yr -1 with the largest uncertainty in the tropics. Seasonal variations in independent SIF estimates agree better with FLUXCOM GPP (mean global pixel-wise R 2 ~0.75) than with GPP from DGVMs (mean global pixel-wise R 2 ~0.6). Seasonal variations in FLUXCOM NEE show good consistency with atmospheric inversion-based net land carbon fluxes, particularly for temperate and boreal regions (R 2 >0.92). Interannual variability of global NEE in FLUXCOM is underestimated compared to inversions and DGVMs. The FLUXCOM version which also uses meteorological inputs shows a strong co-variation in interannual patterns with inversions (R 2 =0.87 for 2001–2010). Mean regional NEE from FLUXCOM shows larger uptake than inversion and DGVM-based estimates, particularly in the tropics with discrepancies of up to several hundred grammes of carbon per square metre per year. These discrepancies can only partly be reconciled by carbon loss pathways that are implicit in inversions but not captured by the flux tower measurements such as carbon emissions from fires and water bodies. We hypothesize that a combination of systematic biases in the underlying eddy covariance data, in particular in tall tropical forests, and a lack of site history effects on NEE in FLUXCOM are likely responsible for the too strong tropical carbon sink estimated by FLUXCOM. Furthermore, as FLUXCOM does not account for CO 2 fertilization effects, carbon flux trends are not realistic. Overall, current FLUXCOM estimates of mean annual and seasonal cycles of GPP as well as seasonal NEE variations provide useful constraints of global carbon cycling, while interannual variability patterns from FLUXCOM are valuable but require cautious interpretation. Exploring the diversity of Earth observation data and of machine learning concepts along with improved quality and quantity of flux tower measurements will facilitate further improvements of the FLUXCOM approach overall.

54 ENVIRONMENTAL SCIENCES↗

Modeling the Sources and Transport Processes During Extreme Ammonia Episodes in the U.S. Corn Belt

Abstract Atmospheric ammonia (NH 3 ) is the primary form of reactive nitrogen (N r ) and a precursor of ammonium (NH 4 + ) aerosols. Ammonia has been linked to adverse impacts on human health, the loss of ecosystem biodiversity, and plays a key role in aerosol radiative forcing. The midwestern United States is the major NH 3 source in North America because of dense livestock operations and the high use of synthetic nitrogen fertilizers. Here, we combine tall‐tower (100 m) observations in Minnesota and Weather Research and Forecasting model coupled with Chemistry (WRF‐Chem) modeling to investigate high and low NH 3 emission episodes within the U.S. Corn Belt to improve our understanding of the distribution of emission sources and transport processes. We examined observations and performed model simulations for cases in February through November of 2017 and 2018. The results showed the following: (1) Peak emissions in November 2017 were enhanced by above‐normal air temperatures, implying a Q 10 (i.e., the change in NH 3 emissions for a temperature increase of 10°C) of 2.5 for emissions. (2) The intensive livestock emissions rom northern Iowa, approximately 400 km away from the tall tower, accounted for 17.6% of the abundance in tall‐tower NH 3 mixing ratios. (3) Ammonia mixing ratios in the innermost domain 3 frequently (i.e., 336 hr, 48% of November 2017) exceeded 5.3 ppb, an important air quality health standard. (4) In November 2017, simulated NH 3 net ecosystem exchange (the difference between NH 3 emissions and dry deposition) accounted for 60–65% of gross NH 3 emissions for agricultural areas and was 2.8–3.1 times the emissions of forested areas. (5) We estimated a mean annual NH 3 net ecosystem exchange of 1.60 ± 0.06 nmol · m −2 · s −1 for agricultural lands and −0.07 ± 0.02 nmol · m −2 · s −1 for forested lands. These results imply that future warmer fall temperatures will enhance agricultural NH 3 emissions, increase the frequency of dangerous NH 3 episodes, and enhance dry NH 3 deposition in adjacent forested lands.

Hu, Cheng↗

The Kokkos Ecosystem [Brief]

In 2016/2017, the field of High-Performance Computing (HPC) entered a new era driven by fundamental physics challenges to produce ever more energy and cost-efficient processors. Since the convergence on the Message-Passing Interface (MPI) standard in the mid-1990s, application developers enjoyed a seemingly static view of the underlying machine — that of a distributed collection of homogeneous nodes executing in collaboration. However, after almost two decades of dominance, the sole use of MPI to derive parallelism acted as a limiter to improved future performance. While MPI is widely expected to continue to function as the basic mechanism for communication between compute nodes for the immediate future, additional parallelism is required on the computing node itself if high performance and efficiency goals are to be realized. When reviewing the architectures of the top HPC systems today, the change in paradigm is clear: the compute nodes of the leading machines in the world are either powered by many-core chips with a few dozen cores each, or use heterogeneous designs, where traditional CPUs marshal work to massively parallel compute accelerators which has as many as 200,000 processing threads in flight simultaneously. Complicating matters further for application developers, each processor vendor has its own preferred way of writing code for their architecture.The Kokkos EcoSystem was released by Sandia in 2017 to address this new era in HPC system design by providing a vendor independent performance portable programming system for scientific, engineering, and mathematical software applications written in the C++ programming language. Using Kokkos, application developers can be more productive because they will not have to create and maintain separate versions of their software for each architecture, nor will they have to be experts in each architecture's peculiar requirements. Instead, they will have a single method of programming for the diverse set of modern HPC architectures. While Kokkos started in 2011 as a programming model only, it soon became clear that complex applications needed more. It is also critical to have a portable mathematical functions and developers need tools to debug their applications, gain insight into the performance characteristics of their codes and tune algorithm performance parameters through automated processes. The Kokkos EcoSystem addresses those needs through its three main components: the Kokkos Core programming model, the Kokkos Kernels math library, and the Kokkos Tools project.

97 MATHEMATICS AND COMPUTING↗

Hydraulic Conductivity Measurements, Utqiagvik (Barrow), Alaska, 2014

Six individual ice cores were collected from the Barrow Environmental Observatory in Barrow, Alaska, in May of 2013 as part of the Next Generation Ecosystem Experiment (NGEE). Each core was drilled at a different location to varying depths. After drilling, the cores were stored in coolers packed with dry ice and flown to Lawrence Berkeley National Laboratory (LBNL) in Berkeley, CA. 3-dimensional images of the cores were constructed using medical X-ray computed tomography (CT) scanner at 120kV. Hydraulic conductivity samples were extracted from these cores at LBNL Richmond Field Station in Richmond, CA, in February 2014 by cutting 5 to 8 inch segments using a chop saw. Samples were packed individually and stored at -20C freezing temperatures to minimize any changes in structure or loss of ice content prior to analysis. Hydraulic conductivity was determined through falling head tests using a permeameter [ELE International, Model #: K-770B] (Appendix A). Samples were placed in a latex membrane via a membrane stretcher while frozen. Use of a membrane stretcher made the membranes easier to secure and minimized contact with the sample. A clear polycarbonate sleeve, fabricated with a stainless steel ring at the bottom to keep the sleeve from floating, was placed around the sample inside the permeameter to minimize deformation during analysis. The permeameter was filled with water and 1.0 PSI of air was applied for confining pressure during sample defrost. Outflow valves were left open to allow for incremental thawing and samples were left to thaw for approximately 12 hours. After approximately 12 hours of thaw, initial falling head tests were performed. When the flow was significantly too fast or too slow, the analysis was stopped and the burette size was adjusted accordingly (i.e. a larger diameter burette was used for flows that were faster than desired or a smaller diameter burette was used for flows that were slower than desired). Two to four measurements were collected on each sample and collection stopped when the applied head load exceeded 25% change from the original load. Analyses were performed between 2 to 3 times for each sample. The final hydraulic conductivity calculations were computed using methodology of Das et al., 1985.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Impacts of Alternative Operations and Renewable Energy Deployment on Columbia River Hydropower

The Columbia River Treaty Tribes in the Pacific Northwest - the Nez Perce, Umatilla, Warm Springs, and Yakama - hold treaty-reserved fishing rights for the Columbia River, one of the world's most productive salmon rivers and a critical resource for these tribes. However, the tribes have expressed that current operating regimes of hydropower dams throughout the Columbia River Basin (CRB) do not fully account for tribal fishing rights and have negatively impacted fish populations. Four tribes acting together through the Columbia River Inter-Tribal Fish Commission (CRITFC)'s 2022 Energy Vision stated that current power system models often do not fully account for the many constraints faced by hydropower facilities in the CRB. As the deployment of variable renewable energy (VRE) like solar and wind continues to increase, power system flexibility will become increasingly important. As a result, there is a growing need to better understand the true capabilities of the hydropower generation fleet while accurately accounting for ecological constraints (Northwest Power and Conservation Council 2022). Understanding hydropower's role in a future grid with a higher share of VRE can inform water resources planning to address ecological needs. The goal of this study is to examine the impacts of alternative hydropower operation rules and weather variability on hydropower generation and grid operations in VRE and transmission infrastructure deployment scenarios. The study evaluates how hydropower scheduling practices can reduce ecological impacts to the region's salmon populations. More specifically, the study examines the impact of today's dam water release rules (called "adjusted base water rules" in this study) and new, ecologically informed water release rules (called "ecosystem water rules" in this study) on two VRE scenarios (current renewable energy levels and higher renewable energy levels) using six weather years for each scenario (2008-2013). CRITFC developed the ecosystem water rules, which capture a portion of the changes they recommend for CRB hydropower operations. These analyses use a water resource planning model (OASIS) and a production cost model (PLEXOS) to simulate CRB reservoir cascade and Western Interconnection power grid operation.

13 HYDRO ENERGY↗

Editorial: Advances in observations and modeling of snow, forest-snow processes and snow hydrology

Seasonal snowpack plays a pivotal role in the regional water cycle. In Western mountainous regions, where precipitation predominantly occurs as snow during winter and is limited in the summer months when water demand peaks, snowpack serves as a critical water source essential for sustaining water supply, ecosystems, agriculture, and industries. Additionally, snowpack dynamically influences extreme events, such as floods, droughts, and wildfires, through intricate land-atmosphere feedback. Nevertheless, significant knowledge gaps persist within this field due to data limitations and the complexity of the physics involved in modeling snow and snow hydrology (i.e., snow acts as a critical upstream hydrological process). Here, this research topic presents a collection of manuscripts that address some of these important research gaps.

54 ENVIRONMENTAL SCIENCES↗

Weakened Orographic Influence on Cool-Season Precipitation in Simulations of Future Warming Over the Western US

High-resolution regional climate model (RCM) simulations of global warming consistently predict larger percentage increases in precipitation in the lee of midlatitude mountain ranges than on their windward slopes, indicating a weakening of the orographic rain shadow. This redistribution of precipitation could have profound consequences for water resources and ecosystems, but its underlying mechanisms are unknown. Here we show that rain-shadow weakening is just one manifestation of a more general decrease in the influence of orography on precipitation under global warming. We introduce a simple model of precipitation change based on this principle, and find that it agrees well with an ensemble of high-resolution simulations performed over the western United States. We argue that diminished orographic influence can be explained by the unique vertical structure of orographically forced ascent, which tends to maximize in the lower atmosphere where condensation is thermodynamically less sensitive to warming.

54 ENVIRONMENTAL SCIENCES↗

Manuscript model outputs, model source code, and figure scripts: The role of geomorphology in mediating biomass allocation impacts on salt-marsh resilience and carbon accumulation

This data package provides model source code (C++), model outputs, and figure generation scripts (R) needed to reproduce the following manuscript: Bruns, Nicholas E., Genevieve L. Noyce, and Matthew L. Kirwan. "The Role of Geomorphology in Mediating Biomass Allocation Impacts on Salt-Marsh Resilience and Carbon Accumulation." Estuarine, Coastal and Shelf Science 327 (December 2025): 109549. https://doi.org/10.1016/j.ecss.2025.109549. This manuscript investigates how geomorphology mediates the impact of biomass allocation shifts on salt marsh persistence and carbon (C) sequestration under sea level rise. We use a 1-D soil-column model (Kirwan and Mudd 2012) to perform experiments across a range of static root:shoot ratios (RSR = 1-4) spanning observed values. The model explicitly simulates interactions between tidal inundation, productivity, inorganic sediment deposition, and organic matter accumulation. Experiments determine whether geomorphic feedbacks amplify, dampen, or leave unchanged the ecosystem response to biomass allocation shifts. A first experiment uses constant sea level rise (2.5 mm/yr) to examine equilibrium dynamics and their influence on carbon accumulation rates across different suspended sediment concentrations (SSC). A second set of experiments calculates threshold sea level rise rates for marsh drowning across RSR and SSC combinations. Final experiments apply accelerating sea level rise scenarios (2000-2200) derived from NOAA projections (Sweet et al. 2022) to generate an envelope of expected responses, quantifying the importance of biomass allocation shifts on marsh carbon accumulation and persistence. All experiments are repeated across SSC ranging from 5-50 mg/L to investigate how these interactions vary in micro-tidal marshes with different sediment supplies. Package contents: * README.txt with detailed description of package contents and instructions for reproducing manuscript figures and model outputs * R scripts for generating all manuscript figures * C++ baseline model code from Kirwan and Mudd (2012) * C++ source code for 4 experimental model variants used in the manuscript, extending above baseline code * Model input files (.csv, .txt) including sea level rise scenarios * Model output files (.txt) used in manuscript analyses Temporal coverage: Model simulations span years 1900-2200, with accelerating sea level rise scenarios for 2000-2200. Key variables: Root:shoot ratio, suspended sediment concentration, carbon accumulation rate, vertical accretion rate, marsh elevation, inundation depth, threshold sea level rise rate.

54 ENVIRONMENTAL SCIENCES↗

Aerobic bacterial methane synthesis

Significance Observations summarized herein contribute to an ongoing paradigm shift in microbial ecology, documenting an emergent property of ecosystem function that further challenges the perception that biogenic methane (CH 4 ) production is strictly an anaerobic process. Relevant metabolites, a model bacterial isolate, gene, and enzyme are identified, and we show how this property can conceivably be broadly distributed in the biosphere and contribute to global CH 4 emissions. Scientifically, this study will enable lines of investigation that will expand our understanding of CH 4 synthesis and emission in nature and illustrates how CH 4 synthesis may actually serve as a nexus for the C and N cycles in nature.

Wang, Qian↗

Non-growing season plant nutrient uptake controls Arctic tundra vegetation composition under future climate

Abstract Plant growth and distribution in high-latitude tundra ecosystems is strongly limited by nutrient availability and is critical for quantifying centennial-scale carbon-climate interactions. However, land model representations of plant–nutrient interactions are uncertain, leading to poor comparisons with high-latitude observations. Although it has been recognized for decades in the observational community that plants continue to acquire nutrients well past when aboveground activity has ceased, most large-scale land models ignore this process. Here we address the role tundra plant nutrient acquisition during the non-growing season (NGS) has on centennial-scale vegetation growth and dynamics, with a focus on shrub expansion. We apply a well-tested mechanistic model of coupled plant, microbial, hydrological, and thermal dynamics that explicitly represents nutrient acquisition based on plant and microbial traits, thereby allowing a prognostic assessment of NGS nutrient uptake. We first show that the model accurately represents observed seasonality of NGS plant nutrient uptake in a northern Alaskan tundra site. Applying the model across the North America tundra indicates that NGS nutrient uptake is consistent with observations and ranges between ∼5% and 50% of annual uptake, with large spatial variability and dependence on plant functional type. We show that NGS plant nutrient acquisition strongly enhances modeled 21 st century tundra shrub growth and expansion rates. Our results suggest that without NGS nutrient uptake, total shrub aboveground dominance would be ∼50% lower, limited primarily by their inability to grow tall enough to maximize their inherent capacity for light competition. Evergreen shrubs would be more strongly affected because of their relatively lower capacity for nutrient remobilization and acquisition compared to deciduous shrubs. Our results highlight the importance of NGS plant and soil processes on high-latitude biogeochemistry and vegetation dynamics and motivate new observations and model structures to represent these dynamics.

54 ENVIRONMENTAL SCIENCES↗

Seabird distribution is better predicted by abundance of prey than oceanography. A case study in the Gulf of Cadiz (SW, Iberian Peninsula)

Abstract Quantifying factors that influence marine predator distributions is essential to understanding the current and future change in marine biodiversity. Here, we test whether marine predator distribution relates to prey, or is better predicted by other specific habitat features. We examine the correlation between spatial distribution of three seabird species and their prey, as well as environmental proxies (oceanographic characteristics) in the Gulf of Cadiz, NE Atlantic. We modeled the at-sea distribution of Cory's shearwater, Balearic shearwater and Northern gannet, based on: (i) pelagic fish abundance according to acoustic surveys, and (ii) a forecast-model of remotely sensed environmental variables (productivity, sea surface temperature, and salinity). In general, seabird distributions were better predicted by abundance of fish than by environmental variables at the habitat scale. We obtained consistent correlations between seabird presence and the abundance of medium-sized (10–20 cm) sardines, anchovies and Mediterranean horse mackerel, providing information on their preferred prey. Additionally, oceanographic productivity variables moderately contributed to seabird distribution models, with better predictive value for the critically endangered Balearic shearwater and Northern gannet in the summer, whilst the model for Cory's shearwater's produced poorer predictions. Predator–prey combined studies may represent essential tools for an efficient ecosystem-based management of marine environments.

de la Cruz, Andrés (ORCID:0000000203632498)↗

Multiscale Land-Atmosphere Interaction in Tropical Ecosystems

A leading source of uncertainty in the magnitude of future climate change arises from feedback between the physical climate system and the carbon cycle. Output from 15 Earth System Models (ESMs) in the Coupled Model Intercomparison Project (CMIP5) shows that when forced by identical emissions, predicted atmospheric CO 2 concentrations in 2100 ranged from 780 ppm to 1120 ppm due to variations in carbon cycle feedback among the models. The resulting difference in radiative forcing of climate is almost 2 W m -2 , comparable to the total uncertainty in the physical climate system for the representative concentration pathway (RCP) 8.5 scenario in 2100. Scaling future carbon sources and sinks to match contemporary atmospheric CO 2 reduced this uncertainty by a factor of five, suggesting that contemporary observations might dramatically constrain future predictions, but their scaling had no mechanistic basis, so it provides little insight on carbon cycle processes.

54 ENVIRONMENTAL SCIENCES↗

Factors controlling marine aerosol size distributions and their climate effects over the northwest Atlantic Ocean region

Abstract. Aerosols over Earth's remote and spatially extensive ocean surfaces have important influences on planetary climate. However, these aerosols and their effects remain poorly understood, in part due to the remoteness and limited observations over these regions. In this study, we seek to understand factors that shape marine aerosol size distributions and composition in the northwest Atlantic Ocean region. We use the GEOS-Chem model with the TwO-Moment Aerosol Sectional (TOMAS) microphysics algorithm model to interpret measurements collected from ship and aircraft during the four seasonal campaigns of the North Atlantic Aerosols and Marine Ecosystems Study (NAAMES) conducted between 2015 and 2018. Observations from the NAAMES campaigns show enhancements in the campaign-median number of aerosols with diameters larger than 3 nm in the lower troposphere (below 6 km), most pronounced during the phytoplankton bloom maxima (May/June) below 2 km in the free troposphere. Our simulations, combined with NAAMES ship and aircraft measurements, suggest several key factors that contribute to aerosol number and size in the northwest Atlantic lower troposphere, with significant regional-mean (40–60∘ N and 20–50∘ W) cloud-albedo aerosol indirect effect (AIE) and direct radiative effect (DRE) processes during the phytoplankton bloom. These key factors and their associated simulated radiative effects in the region include the following: (1) particle formation near and above the marine boundary layer (MBL) top (AIE: −3.37 W m−2, DRE: −0.62 W m−2); (2) particle growth due to marine secondary organic aerosol (MSOA) as the nascent particles subside into the MBL, enabling them to become cloud-condensation-nuclei-sized particles (AIE: −2.27 W m−2, DRE: −0.10 W m−2); (3) particle formation and growth due to the products of dimethyl sulfide, above and within the MBL (−1.29 W m−2, DRE: −0.06 W m−2); (4) ship emissions (AIE: −0.62 W m−2, DRE: −0.05 W m−2); and (5) primary sea spray emissions (AIE: +0.04 W m−2, DRE: −0.79 W m−2). Our results suggest that a synergy of particle formation in the lower troposphere (particularly near and above the MBL top) and growth by MSOA contributes strongly to cloud-condensation-nuclei-sized particles with significant regional radiative effects in the northwest Atlantic. To gain confidence in radiative effect magnitudes, future work is needed to understand (1) the sources and temperature dependence of condensable marine vapors forming MSOA, (2) primary sea spray emissions, and (3) the species that can form new particles in the lower troposphere and grow these particles as they descend into the marine boundary layer.

54 ENVIRONMENTAL SCIENCES↗

Multi-physics modeling of injected nanoparticles effect on remediation of CO 2 leakage through cracks

One potential risk in CO 2 sequestration is the leakage of carbon dioxide, which can result in contamination of underground water, creating potential threats to existing ecosystems. The common leakage pathway is through the pre-existing fractures or discontinuities within cement in the wellbore, incurred by the environmental conditions imposed on the cement. Injecting nanoparticles into pre-existing cracks is one of the most recently proposed ideas for mitigating fracture propagation in cement CO 2 sequestration. To demonstrate the feasibility of this new technology, a numerical approach was taken in this work, as it is challenging to investigate it in a laboratory setting. We proposed a coupled ALE (Arbitrary Lagrangian–Eulerian)–DEM (Discrete Element Method)–peridynamic modeling strategy within the LS-Dyna package to investigate the intertwined interaction among the CO 2 fluid flow, native fluid (e.g., brine), particle clusters, and cracks within the cement. The numerical results demonstrate that injected nanoparticles can effectively reduce the pressure exerted on the crack surface. Accordingly, the potential fracture propagation at the crack tip would be reduced compared to corresponding cases without nanoparticels as pressurized by fluid flow. This result verifies the effectiveness of proposed nanoparticle injection technology. Finally, using this established modeling strategy, the effect of filling particles on the fracture mitigation for different crack geometries (e.g. particle cluster patterns, aspect ratio of crack aperture and length) and CO 2 reservoir pressure are examined. The result shows that injected particles successfully reduce the fracture propagation in these scenarios.

42 ENGINEERING↗

Effects of Rapid Permafrost Thaw on CO2 and CH4 Fluxes in a Warmer and Wetter Future (Final Technical Report)

When ice-rich permafrost thaws, the ground subsides, creating thermokarst landscapes with dramatically different soil conditions and carbon fluxes than the original ecosystem. Roughly 20% of the northern permafrost region is susceptible to thermokarst formation (Olefeldt et al., 2016). Thermokarst formation is often rapid; tens of meters of permafrost can thaw within a few years (Schuur et al., 2015). In topographically low areas, thermokarst thaw converts boreal forest or tundra dry shrub ecosystems into sedge or Sphagnum moss wetlands (Olefeldt et al., 2016). While this type of landscape transformation releases carbon stored in permafrost into the atmosphere, on longer time scales, it facilitates sequestration of atmospheric carbon in plant biomass because permafrost thaw releases plant-available nutrients and wetlands are highly productive (M. C. Jones et al., 2017). However, wetlands also generate methane, which is a potent greenhouse gas. Methane emissions from thermokarst wetlands can cause these carbon-sequestering systems to have a positive global warming potential (Johansson et al., 2006; Turetsky et al., 2007). Our project objective was to improve Earth System and environmental predictability by advancing understanding of how CO 2 and CH 4 flux in permafrost thaw-induced wetlands (thermokarst) will change in the future as temperatures and climate conditions shift. Northern latitudes are expected to get warmer and wetter (IPCC 2013), and initiation and expansion of thermokarst thaw is expected to increase (Jorgenson et al. 2006; Zhang et al. 2017). Given these expected changes, our work sought to address three broad questions: Q1) How will northern latitude CO2 and CH4 emissions respond to warming temperatures? Q2) What is the impact of precipitation on permafrost thaw and carbon emissions? Q3) How do CO2 and CH4 emissions change as wetlands age after permafrost thaw? To answer these questions, we took both a modeling and measurement approach. Modeling work was conducted with DOE’s Energy Exascale Earth System Model (E3SM) land model (ELM). Empirical work took place in two primary locations. The first was a thermokarst site near Fairbanks, AK. The site is part of the Bonanza Creek Long Term Ecological Research program and is well instrumented (Neumann et al., 2019). The second was an isolated thawing permafrost wetland located on Kenai Peninsula — Brown’s Lake bog (B. M. Jones et al., 2016) — where the current climate is representative of what is expected at higher latitudes in the future. This site provides an ideal opportunity to test our hypotheses about how thermokarst wetland dynamics will respond to future environmental conditions. In addition, the project collaborated with researchers asking similar questions who were collecting measurements in high latitude post-glacial lakes located in collaborators tackling similar questions in Sweden (Stordalen Mire) (Emerson et al., 2021). Project efforts directly aligned with the stated goal of the funding opportunity announcement (FOA), which was “to improve the understanding and representation of terrestrial ecosystems in ways that advance Earth system model parameterizations and capabilities... thereby improving the quality of Earth and environmental model projections and providing the scientific foundation needed to support DOE’s science and energy missions.” Specifically, the project improved sophistication and accuracy of the Energy Exascale Earth System Model (E3SM), which is being developed primarily at DOE National Laboratories to support scientific research and decision-making.

54 ENVIRONMENTAL SCIENCES↗