Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earth Systems”

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,045 records · Page 58

Multirate partitioned Runge–Kutta methods for coupled Navier–Stokes equations

Earth system models are complex integrated models of atmosphere, ocean, sea ice, and land surface. Coupling the components can be a significant challenge due to the difference in physics, temporal, and spatial scales. Further, this study explores multirate partitioned Runge-Kutta methods for the fluid-fluid interaction problem and demonstrates its parallel performance by using the PETSc library. We consider compressible Navier-Stokes equations with gravity coupled through a rigid-lid interface. Our large-scale numerical experiments reveal that multirate partitioned Runge-Kutta coupling schemes (1) can conserve total mass; (2) have second-order accuracy in time; and (3) provide favorable strong- and weak-scaling performance on modern computing architectures. We also show that the speedup factors of multirate partitioned Runge-Kutta methods match theoretical expectations over their base (single-rate) method.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the reductions in mortality from air-pollution by cancelling new coal power plants

Deep decarbonization paths to the 1.5°C or 2°C temperature stabilization futures require a rapid reduction in conventional coal-fired power plants, but countries are currently building 223 GW of new coal power capacity and plan to build 377 GW more in the next decade. Coal-fired plants are also a major contributor to air pollution related health impacts. Here, we couple an integrated human-earth system model (GCAM) with an air quality model (TM5-FASST) to examine regional health co-benefits from cancelling new coal-fired plants worldwide. Additionally, we find that cancelling all new proposed projects would decrease air pollution related premature mortality between 101,388 - 213,205 deaths (2-5%) in 2030, and 213,414 - 373,054 (5-8%) in 2050, globally, but heavily concentrated in developing Asia. Furthermore, we estimate that strengthening the climate target from 2°C to 1.5°C would avoid 326,351 additional mortalities in 2030, of which 251,011 (75%) are attributable to the incremental coal plant shutdown.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

The short-term comprehensive impact of the phase-out of global coal combustion on air pollution and climate change

With the continuous intensification of global warming, the reduction and ultimate phase-out of coal combustion is an inevitable trend in the future global energy transformation. This study comprehensively analyzed the impact of phasing out coal combustion on global emissions and concentrations of air pollutants, radiative fluxes, meteorology and climate using Community Earth System Model 2 (CESM2). The results indicate that after the global phase-out of coal combustion, there is a marked decrease in the concentrations of sulfur dioxide (SO 2 ), nitrogen oxides (NO x ) and fine particulate matter (PM 2.5 ), with some regions experiencing a reduction of exceeding 50%. There is no significant change in global ozone (O 3 ) concentration. There are decreasing AOD and positive radiative fluxes globally in the short term, though the cloud contributes minor negative radiative fluxes. The global air temperature may increase by approximately (0.02 ± 0.15) °C on average with regional and seasonal variations, and the precipitation may potentially increase by approximately (2.7 ± 40.6) mm yr −1 globally and over 20% in equatorial regions in the short term. But combined with the decreasing trend of cloud water content in the Northern Hemisphere, it indicates a potential increase in the extremity of precipitation events. In conclusion, this study provides references for global control of air pollution, mitigation strategies of climate change, and transformation of energy structures under the objective of “carbon neutrality”, such as focusing on the negative climate impacts of exacerbating regional warming and increasing extreme precipitation resulting from the rapid reduction of aerosols in the short term.

01 COAL, LIGNITE, AND PEAT↗

A Comprehensive Chemistry Evaluation and Diagnostics Package for E3SM – ChemDyg Version 1.1.0

The Chemistry Evaluation and Diagnostics Package (ChemDyg) is an open-source tool designed for the Energy Exascale Earth System Model (E3SM) developed by the U.S. Department of Energy. ChemDyg facilitates routine evaluation, tailored development, and in-depth analysis of atmospheric chemistry through its modular architecture, allowing users to compare model outputs with observational data. Version 1.1.0 introduces a robust set of diagnostic capabilities, including climatology, time evolution of key tracers, diurnal and annual cycle analyses, and extensive budget diagnostics. These features help identify model discrepancies and enhance the representation of atmospheric chemistry in E3SM. Each self-contained diagnostic set includes dedicated scripts and documentation for ease of use. The interactive HTML output improves data accessibility, accelerating chemistry model development. Additionally, ChemDyg's flexible framework allows for customization, enabling users to create unique diagnostic sets for specific scientific contributions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dissolved pyrogenic carbon leaching in soil: Effects of soil depth and pyrolysis temperature

We report Pyrogenic organic carbon (PyC) is found in soils as a heterogeneous mixture of thermally altered plant residues that range in their susceptibility to losses via mineralization, leaching, and erosion. Leaching of PyC as DOC (DPyC) within the soil profile is likely influenced by the chemical composition of solid PyC and the dominant soil processes and properties at a particular depth. Here we report the results of a 2-year laboratory decomposition and leaching study designed to investigate the interactive effects of pyrolysis temperature levels (no pyrolysis, 300 °C, and 450 °C) and soil depth (surface soil, 0–10 cm depth; subsurface soil, 50–70 cm depth) on losses of 13 C-labeled jack pine (Pinus banksiana) wood, wood pyrolyzed at 300 °C (PyC300) and wood pyrolyzed at 450 °C (PyC450) as dissolved organic carbon (DOC). Losses of wood and PyC in the 13 C- DOC pool were measured in leachates drained from soils once a month over the course of 1 year, and at one single leaching event after 2 years. We found that pyrolysis temperature levels interacted with time and soil depth to affect losses of wood and PyC as DOC throughout the 1-year incubation and leaching study, with greater DOC losses from wood in surface soil (0.73 ± 0.076% of added C) than in subsurface soil (0.40 ± 0.063% of added C) averaged across sampling time. Monthly DOC losses from PyC were not affected by soil depth. Cumulative data indicated a small contribution of wood (2.86 ± 0.07%), PyC300 (0.40 ± 0.04%), and PyC450 (0.16 ± 0.01%) to total DOC leached from soils. DOC losses from wood and PyC300 (as a proportion of added C) were greater in surface soil than in subsurface soil, whereas DOC losses from PyC450 were unaffected by differences between surface and subsurface soil. Losses of DOC were greater from PyC300 than from PyC450 in surface soil, with no significant differences in DOC losses between PyC300 and PyC450 in subsurface soil. One single leaching event at the end of the 2-year decomposition study resulted in higher DOC losses from wood in subsurface soil than in surface soil likely due to desorption, and no differences in DOC losses between PyC300 and PyC450 regardless of soil depth. Our results suggest strong interactions between the initial physicochemical composition of organic C inputs and soil properties (soil depth as a proxy) that control the mobility and transport of PyC and should be better represented in Earth system models.

54 ENVIRONMENTAL SCIENCES↗

Spatiotemporal dynamics of ecosystem fires and biomass burning-induced carbon emissions in China over the past two decades

Fire is a major type of disturbance that has important influences on ecosystem dynamics and carbon cycles. Yet our understanding of ecosystem fires and their carbon cycle consequences is still limited, largely due to the difficulty of large-scale fire monitoring and the complex interactions between fire, vegetation, climate, and anthropogenic factors. Here, using data from satellite-derived fire observations and ecosystem model simulations, we performed a comprehensive investigation of the spatial and temporal dynamics of China's ecosystem fire disturbances and their carbon emissions over the past two decades (1997–2016). Satellite-derived results showed that on average about 3.47 - 4.53 x 10 4 km 2 of the land was burned annually during the past two decades, among which annual burned forest area was about 0.81 - 1.25 × 10 4 km 2 , accounting for 0.33-0.51% of the forest area in China. Biomass burning emitted about 23.02 TgC per year. Compared to satellite products, simulations from the Energy Exascale Earth System Land Model (ELM) strongly overestimated China's burned area and fire-induced carbon emissions. Annual burned area and fire-induced carbon emissions were high for boreal forest in Northeast China's Daxing'anling region and subtropical dry forest in South Yunnan, as revealed by both the satellite product and the model simulations. Our results suggest that climate and anthropogenic factors play critical roles in controlling the spatial and seasonal distribution of China's ecosystem fire disturbances. Our findings highlight the importance of multiple complementary approaches in assessing ecosystem fire disturbance and its carbon consequences. Further studies are required to improve the methods of observing and modelling China's ecosystem fire disturbances, which will provide valuable information for fire management and ecosystem sustainability in an era when both human activities and the natural environment are rapidly changing.

54 ENVIRONMENTAL SCIENCES↗

Towards the next generation of Geospatial Artificial Intelligence

Geospatial Artificial Intelligence (GeoAI), as the integration of geospatial studies and AI, has become one of the fastest-developing research directions in spatial data science and geography. This rapid change in the field calls for a deeper understanding of the recent developments and envision where the field is going in the near future. In this work, we provide a quantitative analysis of the GeoAI literature from the spatial, temporal, and semantic aspects. We briefly discuss the history of AI and GeoAI by highlighting some pioneering work. Then we discuss the current landscape of GeoAI by selecting five representative subdomains including remote sensing, urban computing, Earth system science, cartography, and geospatial semantics. Finally, we highlight several unique future research directions of GeoAI which are classified into two groups: GeoAI method development challenges and GeoAI Ethics challenges. Topics include heterogeneity-aware GeoAI, knowledge-guided GeoAI, spatial representation learning, geo-foundation models, fairness-aware GeoAI, privacy-aware GeoAI, as well as interpretable and explainable GeoAI. We hope our review of GeoAI’s past, present, and future is comprehensive and can enlighten the next generation of GeoAI research.

58 GEOSCIENCES↗

Choices in land representation materially affect modeled biofuel carbon intensity estimates

Estimates of biofuel carbon intensity are uncertain and depend on modeled land use change (LUC) emissions. While analysts have focused on economic and agronomic assumptions affecting the quantity of land converted, researchers have paid less attention to how models classify land into broad categories and designate some categories as ineligible for LUC. To explore the effect of these land representation attributes, we use three versions of a global human and Earth systems model, GCAM, and compute the “carbon intensity of land-use change” (CI-LUC) from increased U.S. corn ethanol production. We consider uncertainty in model parameters along with the choice of land representation and find the latter is one of the most influential parameters on estimated CI-LUC. A version of the model that protects 90% of non-commercial land reduced estimated CI-LUC by an average of 32% across Monte Carlo trials compared to our baseline model. Another version that mimics the GTAP-BIO-ADV land representation, which protects all non-commercial land, reduced CI-LUC by an average of 19%. The results of this experiment suggest that analyses of LUC that exclude non-commercial land categories or assume more land is protected than empirically justified underestimate biofuel induced LUC GHG emissions.

09 BIOMASS FUELS↗

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport

In a computational fluid model of the atmosphere, the advective transport of trace species, or tracers, can be computationally expensive. For efficiency, models often use semi-Lagrangian advection methods. High-order interpolation semi-Lagrangian (ISL) methods, in particular, can be extremely efficient, if the problem of property preservation specific to them can be addressed. Atmosphere models often use geometrically and logically nonuniform grids for efficiency and, as a result, element-based discretizations. Such grids and discretizations make stability a particular problem for ISL methods. Generally, high-order, element-based ISL methods that use the natural polynomial interpolant associated with a nodal finite-element discretization are unstable. Here, we derive new bases having order of accuracy up to nine, with positive nodal weights, that stabilize the element-based ISL method. We use these bases to construct the linear advection operator in the property-preserving Interpolation Semi-Lagrangian Element-based Transport (Islet) method. Then we discuss key software implementation details. Finally, we show performance results for the Energy Exascale Earth System Model's atmosphere dynamical core, comparing the original and new transport methods. These simulations used up to 27,600 Graphical Processing Units (GPU) on the Oak Ridge Leadership Computing Facility's Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗

Impacts of climate change on future hurricane induced rainfall and flooding in a coastal watershed: A case study on Hurricane Harvey

The warming climate is likely to increase hurricane-associated extreme rainfall and lead to sea-level rise (SLR). Thus, how the floods induced by intense hurricanes respond to these potential changes is of great concern. This study investigates the future warmer climate impacts on hurricane-induced extreme rainfall, and—more importantly—the subsequent compound flooding at the watershed scale (from an event-based analysis perspective). To this goal, a modeling framework is designed based on the Distributed Hydrology Soil Vegetation Model (DHSVM), the Two-Dimensional fvand the Regional Community Earth System Model (R-CESM). The framework was applied to Hurricane Harvey (2017) at the Clear Creek watershed (a coastal watershed in the southern Houston) as a case study. The results show that the projected maximum rainfall totals over the watershed would be exacerbated by 17.7 % and 49.7 % in the 2050s and 2090s (respectively) under Representative Concentration Pathway 8.5 (RCP 8.5). This means a 16.1 % increase in Harvey rainfall over the watershed per degree Celsius increase in Mean Surface Temperature over the Gulf of Mexico region (18°~31° N, 77°~98° W). Meanwhile, the increases in maximum inundation extent would be 11.0 % (2050s) and 19.5 % (2090s). Furthermore, considerable increases in maximum inundation depth and duration in regions along the middle and downstream of Clear Creek (and also those around Clear Lake) are expected. The projected SLR will have little effect on the maximum inundation depth and extent if storm surge changes are not taken into account; meanwhile, it will influence the inundation duration at downstream locations. In conclusion, this modeling framework can be also applied at other coastal watersheds to evaluate the projected climate change impacts on the compound flooding induced by extreme climate events.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the long-term changes of terrestrial water storage and their driving factors

Global warming is expected to cause changes in terrestrial water storage (TWS) across the land surface, with widespread impacts on ecosystems and society. Although extensive research has been performed to analyze TWS changes and possible drivers during the post-2000 period, longer-term evolution of TWS and associated environmental forcings remain relatively unexplored. In this study, we evaluated the performance of the Energy Exascale Earth System model (E3SM) land model ELM version 1 (ELM v1) in simulating global TWS, and used factorial simulations of ELMv1 to quantify global TWS changes and their drivers during 1948–2012. We found that ELM’s agreed best with existing satellites and reconstruction datasets in temperate regions unaffected by irrigation. Biome- and climate zone-averaged TWS mainly increased at rates between 0 and 10 mm/year over 1948–2012, but the second half of that period saw smaller positive trends than the first half or even negative trends. Climate change explained >80 % of the TWS trends across most biomes and climate zones, followed by land use and land cover change. The physiological and phenological effects of CO2 primarily induced noticeable TWS trends in the more humid biomes and climate zones across different latitudes. In contrast, nitrogen deposition and aerosol deposition generally had smaller and negative impacts across the biomes and climate regions. Among the meteorological drivers analyzed, the long-term average imbalance between precipitation (P), evapotranspiration (E), and runoff (Q) contributed >50 % of the TWS trends in most biomes and climate zones, with nonlinearity being induced by spatially heterogenous changes in E/P and Q/P ratios. The accumulated detrended anomalies in P, E, and Q also often contributed substantially, while the trends difference between P, E, and Q contributed little. Together, these findings unveiled an intensification of the global TWS and its diverse patterns of climate change and different non-withdrawal human-induced alterations, contributing to a more comprehensive understanding and projection of the global water cycle.

54 ENVIRONMENTAL SCIENCES↗

A hybrid Penman-Monteith and machine learning model for simulating evapotranspiration and its components

Integrating physical processes with machine learning has advanced evapotranspiration (ET) simulation, yet most hybrid models fail to partition total ET into its components: soil evaporation (E) and vegetation transpiration (T). This study introduces Residual Neural Network–Penman–Monteith (RNN-PM), a novel hybrid dual-source ET model designed to overcome this limitation. The model synergizes the physically-based Penman–Monteith framework with three specialized residual neural networks trained to estimate key conductance parameters (canopy conductance, soil surface conductance, and aerodynamic conductance). Furthermore this explicit parameterization allows for the direct partitioning of total ET. Validation at National Ecological Observatory Network (NEON) flux sites using high-frequency partitioned E and T shows that RNN-PM reliably reproduces ET and the transpiration fraction (T/ET). For ET, the model achieves an average Kling–Gupta efficiency (KGE) of 0.89 and a root-mean-square error (RMSE) of 0.55 mm/day; for T/ET, the KGE is 0.87 with an RMSE of 0.06. Furthermore, RNN-PM demonstrates robust generalization, accurately simulating ET and its components well beyond the initial training dataset, even under extreme climatic conditions. This study extended the analysis by comparing the RNN-PM model with seven established dual-source ET models. The results indicate that RNN-PM outperforms both conventional machine learning models and purely physical process-based models in simulating ET components in most cases. Among the purely physical process-based dual-source models, those based on surface temperature decomposition showed improved performance as the leaf area index (LAI) decreased when evaluated against high-frequency ET component datasets. In contrast, the performance of conductance-based dual-source models declined with decreasing LAI. Although purely machine learning-based models can produce relatively accurate simulations of ET components, they often exhibit limited generalization capability, an issue that the RNN-PM model effectively overcomes. Ultimately, the RNN-PM model represents a significant advance in simulating ET components, offering a novel and scalable approach for improving the representation of land–atmosphere interactions in Earth system models.

54 ENVIRONMENTAL SCIENCES↗

Storylines for the 1997 New Year’s Flood: The role of watershed antecedent conditions and future warming in shaping discharge in the Truckee River watershed

The 1997 New Year’s flood was among the most devastating floods in the Truckee River watershed located in western Nevada. This event resulted from complex interactions of flood drivers, such as extreme precipitation, wet antecedent watershed conditions, warm temperatures and rapid snowmelt. We leveraged simulated forcings from the regionally refined mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) and a process-based hydrological model to recreate the 1997 New Year’s flood for the Truckee River watershed across four climate warming levels ranging from the current temperatures to + 4° C. For each scenario, we conducted ensemble simulations with the same forcing but with 100 different seasonal watershed antecedent conditions, which were randomly sampled from long-term hydrological simulations. The results show that the 1997 New Year’s flood can be reproduced or exceeded consistently only when the antecedent watershed conditions are wet, specifically when streamflows are above the 75th percentile of the climatological value. There is negligible change in ensemble mean peakflows for Truckee River near Reno; however, there are increases of 18% and 14% under the warming levels of + 3° C and + 4° C, respectively. The increases in peakflows under future climate warming are attributed to wetter antecedent watershed conditions and enhanced snowmelt. Furthermore, the largest increases in peakflows occur at small, high-elevation headwater basins along the Sierra Nevada crest. This study highlights that changes in extreme flood events will result from the complex interplay of multiple flood drivers. It also demonstrates the potential of storyline approaches to analyze future realizations of these extreme events under different climate scenarios.

Climate change impact study↗

Scalable self attraction and loading calculations for unstructured ocean tide models

Self attraction and earth-loading effects are important for accurately modeling global tides. A common approach of handling this forcing is to expand mass anomalies into spherical harmonics, which are scaled by load Love numbers to account for elastic earth deformation. We investigate two different approaches to perform these calculations for ocean models that employ unstructured meshes and distributed memory parallelization. The first approach leverages a highly efficient spherical harmonics library, but requires all-to-one and one-to-all communications and interpolation operations between the unstructured and a structured mesh. This approach is compared to a parallel algorithm that computes the spherical harmonic transformations directly on the unstructured mesh with an all-reduce communication. Here, our results show that although the unstructured mesh calculations are more expensive, the scalability of the unstructured mesh approach allows for more efficient spherical harmonics transforms for high-resolution meshes and large processor counts. This methodology enables the efficient inclusion of tidal dynamics large-scale Earth system model simulations.

54 ENVIRONMENTAL SCIENCES↗

Distinct Assembly Processes and Microbial Communities Constrain Soil Organic Carbon Formation

Soil stores more carbon (C) than all vegetation and the atmosphere combined. Soil C stocks are broadly shaped by temperature, moisture, soil physical characteristics, vegetation, and microbial-mediated metabolic processes. The efficiency with which microorganisms use soil C regulates the balance between C storage in soil and the atmosphere. In this review, we discuss how microbial physiology and community assembly processes determine microbial growth rate and efficiency and, in turn, soil organic C cycling through the lens of community ecology. We introduce a conceptual framework cataloging life history (i.e., growth rate, resource acquisition, and stress tolerance) and assembly traits (i.e., competition, facilitation, and dispersal) that correspond with different growth efficiencies. We also compare how dominant mycorrhizal fungal type affects growth efficiency. We propose that further development and inclusion of specific community parameters in microbial-explicit Earth system models are needed for accurately predicting soil organic C responses to global change.

carbon sequestration↗

Unveiling dark brown carbon: The hidden light-absorbing aerosols from biomass burning

The global rise in wildfire severity highlights the significant warming impact of brown carbon aerosols. Further, new findings reveal water-insoluble “dark brown carbon” tar balls that absorb far more light than the previously studied water-soluble brown carbon, representing a major unaccounted factor in global Earth System Models.

54 ENVIRONMENTAL SCIENCES↗

Leaf isoprene and monoterpene emission distribution across hyperdominant tree genera in the Amazon basin

Tropical forests are acknowledged to be the largest global source of isoprene (C 5 H 8 ) and monoterpenes (C 10 H 16 ) emissions, with current synthesis studies suggesting few tropical species emit isoprenoids (20-38%) and do so with highly variable emission capacities, including within the same genera. This apparent lack of a clear phylogenetic thread has created difficulties both in linking isoprenoid function with evolution and for the development of accurate biosphere-atmosphere models. In this work, we present a systematic emission study of "hyperdominant" tree species in the Amazon Basin. Across 162 individuals, distributed among 25 botanical families and 113 species, isoprenoid emissions were widespread among both early and late successional species (isoprene: 61.9% of the species; monoterpenes: 15.0%; both isoprene and monoterpenes: 9.7%). The hyperdominant species (69) across the top five most abundant genera, which make up about 50% of all individuals in the Basin, had a similar abundance of isoprenoid emitters (isoprene: 63.8%; monoterpenes: 17.4%; both 11.6%). Among the abundant genera, only Pouteria had a low frequency of isoprene emitting species (15.8% of 19 species). In contrast, Protium, Licania, Inga, and Eschweilera were rich in isoprene emitting species (83.3% of 12 species, 61.1% of 18 species, 100% of 8 species, and 100% of 12 species, respectively). Light response curves of individuals in each of the five genera showed light-dependent, photosynthesis-linked emission rates of isoprene and monoterpenes. Importantly, in every genus, we observed species with light-dependent isoprene emissions together with monoterpenes including β-ocimene. These observations support the emerging view of the evolution of isoprene synthases from β-ocimene synthases. Our results have important implications for understanding isoprenoid function-evolution relationships and the development of more accurate Earth System Models.

59 BASIC BIOLOGICAL SCIENCES↗

A mechanistic understanding of oxygen isotopic changes in the Western United States at the Last Glacial Maximum

At the Last Glacial Maximum (LGM), records suggest drier conditions in the northwest United States and wetter conditions in the southwest United States relative to present-day as well as widespread changes in the isotopic composition of water. However, the mechanisms responsible for these changes remain ambiguous. Here, we explore differences in western United States hydroclimate between the LGM and preindustrial with a water isotope tracer enabled Earth System Model. We then use proxy forward models to compare simulated and recorded δ 18 O in speleothems. We find that the pattern of hydroclimate response in the western United States at the LGM relates to a combination of 1) increased frequency and southward shifted wintertime extratropical cyclones in the North Pacific, 2) greater rainout of moisture as it moves over the continent, and 3) reduced evaporation in the cooler LGM climate. The simulated lower δ 18 O of precipitation at the LGM relates predominantly to an increase in cool season moisture removal efficiency, with a secondary contribution from relatively more cool season precipitation. Both surface temperatures and North American ice sheets contribute to these hydroclimate changes at the LGM. Comparisons between δ 18 O from proxy forward models and speleothem records in the western United States show general agreement at the LGM, with increasing depletion moving towards the continental interior. Finally, this study highlights the similarities and differences between hydrologic and δ 18 O changes at the LGM and emphasizes the utility of model-proxy comparison for interpretation.

54 ENVIRONMENTAL SCIENCES↗