Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earth System Model Development”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Collaborative Research: Improved Efficiency and Coupling of the Radiation Code in the ACME Earth System Model. Final Report

This final report details all work performed on the project by both project partners. This project provided support to properly couple RTE+RRTMGP, a high-performance broadband radiation code, within DOE’s Energy Exascale Earth System Model (E3SM). RTE+RRTMGP is a successor to the RRTMG radiation code, which has been widely accepted for its speed and accuracy by the global modeling community, and has been in use in the NCAR CESM for many years and was implemented in the initial version of E3SM. However, the computational cost of RRTMG remains high relative to other components in part due to its complexity and to its inefficient use of modern optimization strategies, issues that were rectified by the development of RTE+RRTMGP. Many of the accomplishment in this project necessitated significant collaboration with the E3SM development team. One focus of the project was to enhance the code’s optimization on the limited number of emerging computing systems on which the model is expected be used, including Many Integrated Core (MIC) architectures and Graphics Processing Unit (GPU) hardware. We also developed additional capabilities for RTE+RRTMGP that E3SM scientists identified as important for the planned applications of the model. The result of our project was optimization of a key physical component (radiative transfer calculations) of E3SM, directly supporting E3SM’s overarching global modeling objectives. More broadly, this project provided overall advancements in the use of radiative transfer calculations in atmospheric modeling and simulation, particularly for climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Characterization of carbonaceous aerosols during TRACER-CAT

Absorbing aerosols (AA) have an important impact on the global radiation budget and cloud properties. The composition and properties of AA can vary substantially throughout the atmosphere, depending on the particle source and the influence of chemical aging. Uncertainties associated with the radiative effects of AA remain substantial. A key contributor to this uncertainty is understanding the extent to which coatings in general, and water uptake especially, alters absorption by AA particles and how this depends on particle composition. We deployed new and existing experimental tools during the Tracking Aerosol Convection Interactions Experiment (TRACER) campaign in Houston, TX as part of the Carbonaceous Aerosols Thrust (CAT) to provide detailed characterization of aerosol optical, chemical, and physical properties. Our TRACER-CAT measurements complemented and expanded on the planned TRACER instrumentation, allowing for more detailed characterization of aerosol properties of relevance to cloud development (a core focus of TRACER), such as the composition of particles that can act as cloud condensation nuclei, than would otherwise be available. Our measurements have allowed for assessment of the relationship(s) between AA optical properties (with a focus on absorption) and the chemical and physical characteristics (including the mixing state of black carbon (BC) containing particles). These field observations occurred in collaboration with Los Alamos National Laboratory in summer 2022 during the TRACER intensive operating period. The instrumentation we co-deployed provided for measurement of (i) multi-wavelength dry aerosol absorption, scattering, and extinction, (ii) the size-dependent composition and abundance of sub-micron aerosol, differentiating between those particles that do and do not contain BC, (iii) BC-specific concentrations and size distributions, (iv) particle size, and (v) the first field measurements at an ARM site of the influence of RH on multi-wavelength absorption by ambient AA. We have leveraged the natural variability of the atmosphere and of aerosol sources in the Houston region to (i) specifically disentangle contributions to light absorption from BC, absorbing organic carbon (brown carbon), and coatings on BC, (ii) characterize the mixing state of BC and assess the factors that give rise to compositional differences between BC-containing and BC-free aerosol, and (iii) establish how water uptake influences absorption and how any such effect depends on particle composition and BC mixing state. Overall, our study contributed to the mission of the Atmospheric System Research program in multiple ways. Through the deployment of complementary, advanced instrumentation for characterization of a wide range of aerosol properties our work helped to maximize the scientific impact of the TRACER campaign. Our work also allowed for development of new insights into the relationship(s) between aerosol composition, hygroscopicity, and the mixing state of BC with aerosol optical properties. Through this, our work has provided knowledge that can improve understanding and model representation of aerosol processes as they affect the Earth’s radiation budget.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

PAVC: The foundation for a Pan-Arctic Vegetation Cover database

Field-measured Arctic vegetation cover data is essential for creating accurate, high-quality vegetation structure and composition maps. Extrapolating field data into high-resolution cover maps provides detailed, function-specific information for use in Earth System Models, vegetation classifications, and monitoring vegetation change over time and space. However, field campaigns that collect plant cover vary substantially in scope, method, and purpose, which makes them difficult to unify across data stores, and they are often not designed to meet remote sensing needs. In this work, we synthesized and harmonized field-based fractional cover data from various data stores to create a high-quality, consistent repository schema for remote sensing-based vegetation cover mapping applications. We developed a reproducible workflow for synthesizing visual estimate and point-intercept fractional cover data. The resultant Pan-Arctic Vegetation Cover (PAVC) database contains synthesized fractional cover at both the species and plant functional type levels. The latter includes absolute foliar cover for deciduous shrubs and trees, evergreen shrubs and trees, forbs, graminoids, lichen, bryophytes, and “other” vegetation, as well as absolute cover for litter and top cover for water and bare ground.

Steckler, Morgan R. [Oak Ridge National Laboratory↗
DOE OSTI2020

Causal Interactions between Southern Ocean Polynyas and High-Latitude Atmosphere–Ocean Variability

Abstract Weddell Sea open-ocean polynyas have been observed to occasionally release heat from the deep ocean to the atmosphere, indicating that their sporadic appearances may be an important feature of high-latitude atmosphere–ocean variability. Yet, observations of the phenomenon are sparse and many standard-resolution models represent these features poorly, if at all. We use a fully coupled, synoptic-scale preindustrial control simulation of the Energy Exascale Earth System Model (E3SMv0-HR) to effectively simulate open-ocean polynyas and investigate their role in the climate system. Our approach employs statistical tests of Granger causality to diagnose local and remote drivers of, and responses to, polynya heat loss on interannual to decadal time scales. First, we find that polynya heat loss Granger causes a persistent increase in surface air temperature over the Weddell Sea, strengthening the local cyclonic wind circulation. Along with responding to polynyas, atmospheric conditions also facilitate their development. When the Southern Ocean experiences a rapid poleward shift in the circumpolar westerlies following a prolonged negative phase of the southern annular mode (SAM), Weddell Sea salinity increases, promoting density destratification and convection in the water column. Finally, we find that the reduction of surface heat fluxes during periods of full ice cover is not fully compensated by ocean heat transport into the high latitudes. This imbalance leads to a buildup of ocean heat content that supplies polynya heat loss. These results disentangle the complex, coupled climate processes that both enable the polynya’s existence and respond to it, providing insights to improve the representation of these highly episodic sea ice features in climate models.

Kaufman, Zachary S.↗
DOE OSTI2022

Advances in hexagon mesh-based flow direction modeling

Watershed delineation and flow direction representation are the foundations of streamflow routing in spatially distributed hydrologic modeling. A recent study showed that hexagon-based watershed discretization has several advantages compared to the traditional Cartesian (latitude–longitude) discretization, such as uniform connectivity and compatibility with other Earth system model components based on unstructured mesh systems (e.g., oceanic models). Despite these advantages, hexagon-based discretization has not been widely adopted by the current generation of hydrologic models. One major reason is that there is no existing model that can delineate hexagon-based watersheds while maintaining accurate representations of flow direction across various spatial resolutions. In this study, we explored approaches such as spatial resampling and hybrid breaching-filling stream burning techniques to improve watershed delineation and flow direction representation using a newly developed hexagonal mesh watershed delineation model (HexWatershed). We applied these improvements to the Columbia River basin and performed 16 simulations with different configurations. The results show that (1) spatial resampling modulates flow direction around headwaters and provides an opportunity to extract subgrid information; and (2) stream burning corrects the flow directions in mountainous areas with complex terrain features.

58 GEOSCIENCES↗
DOE OSTI2024

Seeing the forest for the trees: implementing dynamic representation of forest management and forest carbon in a long-term global multisector model

Abstract Studies have found that understanding forest management is critical in understanding the interaction between the carbon cycle and the integrated human-Earth system. This makes effectively representing forest management decisions such as planting and harvesting important. Here, we implement a novel dynamic forest harvest model in a global state of the art multi-sector dynamics model, namely the Global Change Analysis Model (GCAM). We implement an approach that explicitly tracks forest age and generates rotation ages for forest harvest that are responsive to changes in wood prices, changes in forest age and regional preferences for forest rotation. Furthermore, the forest sector in GCAM competes for investment with other land use types in the future years based on expected profit. Our baseline scenario results indicate that with the new forest harvest model, the current global wood product demand in GCAM can be met with minimal loss of old growth forest through the age-based harvest decisions. We find that economic pressure for deforestation and consequent loss of forest carbon is a bigger driver of global forest change than wood harvests, especially in developing regions. Under alternative scenarios where an economic value is placed on carbon across the terrestrial and energy systems, while there is an increase in forest plantations, there can be corresponding decreases in forest cover in some regions as forest land competes with land for bio-energy crops. When the carbon in forests is assigned a price, we find that the average rotation age for wood harvests can be reduced across regions to harvest forests in a more carbon efficient manner.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Overlooked cooling effects of albedo in terrestrial ecosystems

Radiative forcing (RF) resulting from changes in surface albedo is increasingly recognized as a significant driver of global climate change but has not been adequately estimated, including by Intergovernmental Panel on Climate Change (IPCC) assessment reports, compared with other warming agents. Here, we first present the physical foundation for modeling albedo-induced RF and the consequent global warming impact (GWI Δα ). We then highlight the shortcomings of available current databases and methodologies for calculating GWI Δα at multiple temporal scales. There is a clear lack of comprehensive in situ measurements of albedo due to sparse geographic coverage of ground-based stations, whereas estimates from satellites suffer from biases due to the limited frequency of image collection, and estimates from earth system models (ESMs) suffer from very coarse spatial resolution land cover maps and associated albedo values in pre-determined lookup tables. Field measurements of albedo show large differences by ecosystem type and large diurnal and seasonal changes. As indicated from our findings in southwest Michigan, GWI Δα is substantial, exceeding the RF Δα values of IPCC reports. Inclusion of GWI Δα to landowners and carbon credit markets for specific management practices are needed in future policies. We further identify four pressing research priorities: developing a comprehensive albedo database, pinpointing accurate reference sites within managed landscapes, refining algorithms for remote sensing of albedo by integrating geostationary and other orbital satellites, and integrating the GWI Δα component into future ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

ARM Aerial Instrument Workshop Report

The mission of the U.S. Department of Energy’s (DOE) Biological and Environmental Research (BER) program is to “support transformative science and scientific user facilities to achieve a predictive understanding of complex biological, earth, and environmental systems for energy and infrastructure security, independence, and prosperity.” (https://science.osti.gov/ber) Aligned with the BER central mission, the Earth and Environmental Systems Sciences Division (EESSD) plays a vital role in supporting the fundamental research to understand and predict Earth’s climate and environmental systems, and is also in a unique position to inform the development of sustainable solutions to the nation’s energy and environmental challenges. Specifically, EESSD manages two scientific user facilities: the Atmospheric Radiation Measurement (ARM) user facility and the Environmental Molecular Sciences Laboratory (EMSL). These facilities provide the broader scientific community with scientific expertise, technical capabilities, and unique data sets to facilitate science in areas of importance to DOE. As a multi-platform scientific user facility, ARM aims to fulfill the needs predominantly within the EESSD Atmospheric System Research (ASR) and the Earth and Environmental System Modeling (EESM) mission areas, and provide the critical measurements required to improve understanding of aerosol and cloud life cycles and their interactions, and their coupling with the Earth’s surface. Over the years, ARM has carried out piloted and unmanned aircraft campaigns under different organizational and operational paradigms (Schmid et al. 2014, 2016). Building on its success, the ARM Aerial Facility (AAF) continues to complement the ground-based observations with airborne in situ cloud, aerosol, and trace gas observations as well as measurements of atmospheric state and atmospheric radiation. During the past three years, ARM has managed field campaigns using unmanned aerial systems (UAS) and tethered balloon systems (TBS) at Oliktok Point in Alaska to improve understanding of atmospheric processes in the Arctic. In 2019, following a careful evaluation of scientific community needs, ARM acquired a Bombardier Challenger 850 regional jet to replace the vintage Grumman Gulfstream-159 turboprop aircraft previously used by AAF. With this new “laboratory in the sky”, AAF is evaluating its current and future aerial observation capabilities to continue satisfying the needs of the research community.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A More La Niña–Like Response to Radiative Forcing after Flux Adjustment in CESM2

Abstract In response to greenhouse gas forcing, most coupled global climate models project the tropical Pacific SST trend toward an “El Niño–like” state, with a reduced zonal SST gradient and a weakened Walker circulation. However, observations over the last five decades reveal a trend toward a more “La Niña–like” state with a strengthening zonal SST gradient. Recent research indicates that the identified trend differences are unlikely to be entirely due to internal variability and probably result, at least in part, from systematic model biases. In this study, Community Earth System Model, version 2 (CESM2), is used to explore how mean-state biases within the model may influence its forced response to radiative forcing in the tropical Pacific. The results show that using flux adjustment to reduce the mean-state bias in CESM2 over the tropical regions results in a more La Niña–like trend pattern in the tropical Pacific, with a strengthening of the tropical Pacific zonal SST gradient and a relatively enhanced Walker circulation, as hypothesized to occur if the ocean thermostat mechanism is stronger than the atmospheric mechanisms which by themselves would weaken the Walker circulation. We also find that the historical strengthening of the tropical Pacific zonal gradient is transient but persists into the near term in a high-emissions future warming scenario. These results suggest the potential of flux adjustment as a method for developing alternative projections that represent a wider range of possible future tropical Pacific warming scenarios, especially for a better understanding of regional patterns of climate risk in the near term.

Zhuo, Jing-Yi [Lamont-Doherty Earth Observatory, C↗
DOE OSTI2021

Water balance response of permafrost-affected watersheds to changes in air temperatures

Abstract Observations show increases in river discharge to the Arctic Ocean especially in winter over the last decades but the physical mechanisms driving these changes are not yet fully understood. We hypothesize that even in the absence of a precipitation increase, permafrost degradation alone can lead to increased annual river runoff. To test this hypothesis we perform 12 millennium-long simulations over an idealized hypothetical watershed (1 km 2 ) using a distributed, physically based water balance model (Water flow and Balance Simulation Model, WaSiM). The model is forced by both a hypothetical warming defined by an air temperature increase of 7.5 ∘ C over 100 years, and a corresponding cooling scenario. To assess model sensitivity we vary soil saturated hydraulic conductivity and lateral subsurface flow configuration. Under the warming scenario, changes in subsurface water transport due to ground temperature changes result in a 7%–14% increase in annual runoff accompanied by a 6%–20% decrease in evapotranspiration. The increase in runoff is most pronounced in winter. Hence, the simulations demonstrate that changes in permafrost characteristics due to climate warming and associated changes in evapotranspiration provide a plausible mechanism for the observed runoff increases in Arctic watersheds. In addition, our experiments show that when lateral subsurface moisture transport is not included, as commonly done in global-scale Earth System Models, the equilibrium water balance in response to the warming or cooling is similar to the water balance in simulations where lateral subsurface transport is included. However, the transient changes in water balance components prior to reaching equilibrium differ greatly between the two. For example, for high saturated hydraulic conductivity only when lateral subsurface transport is considered, a period of decreased runoff occurs immediately after the warming. This period is characterized by a positive change in soil moisture storage caused by the soil moisture deficit developed during prior cooling.

Debolskiy, Matvey V. (ORCID:0000000296343627)↗
DOE OSTI2025

Signatures of aerosol-cloud interactions in GiOcean: a coupled global reanalysis with two-moment cloud microphysics

Aerosols influence the Earth's radiative balance through direct interactions with radiation and by affecting cloud properties. Anthropogenic aerosols have led to cooling during the industrial era through aerosol–cloud interactions (ACI), including aerosol effects on cloud microphysical properties and the subsequent adjustments. However, large uncertainties remain in Earth system models (ESMs) regarding the magnitude of this cooling. In part, ESMs substantially disagree on cloud properties, thermodynamics, the hydrological cycle, and general circulation. Reanalysis provides a useful avenue for exploring the impact of ACI on clouds and radiation because its atmosphere is forced to match realistic conditions through the assimilation of observations. Here, we explore the impact of ACI on clouds in the GiOcean reanalysis – the first to incorporate aerosol-cloud interactions. We contrast variables important for ACI between GiOcean and satellite observations and develop 2-dimensional lookup tables of ACI for both using a source-sink budget perspective to attribute the changes in cloud droplet number (Nd) and liquid water path (LWP) to aerosol and meteorology. A compositing analysis using lookup tables shows that GiOcean captures key aspects of aerosol–cloud–precipitation interactions, including (1) activation of aerosol into cloud droplets, (2) effective precipitation scavenging of Nd, (3) suppression of precipitation by high Nd in regions with heavy aerosol emissions. In contrast, satellite observations do not exhibit clear patterns for processes (2) and (3). Random Forest analysis shows that interannual variability in Nd and LWP over the Northern Hemisphere ocean in GiOcean is primarily driven by precipitation, consistent with satellite observations.

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

Continental United States may lose 1.8 petagrams of soil organic carbon under climate change by 2100

Abstract Aims High‐resolution information on soils’ vulnerability to climate‐induced soil organic carbon (SOC) loss can enable environmental scientists, land managers, and policy makers to develop targeted mitigation strategies. This study aims to estimate baseline and decadal changes in continental US surface SOC stocks under future emission scenarios. Location Continental United States. Time period 2014–2100. Methods We used recent SOC field observations ( n = 6,213 sites), environmental factors ( n = 32), and an ensemble machine learning (ML) approach to estimate baseline SOC stocks in surface soils across the continental United States at 100‐m spatial resolution, and decadal changes under the projected climate scenarios of Coupled Model Intercomparison Project Phase Six (CMIP6) earth system models (ESMs). Results Baseline SOC projections from ML approaches captured more than 50% of variability in SOC observations, whereas ESMs represented only 6–16% of observed SOC variability. ML estimates showed a mean total loss of 1.8 Pg C from US surface soils under the high‐emission scenario by 2100, whereas ESMs showed no significant change in SOC stocks with wide variation among ESMs. Both ML and ESM predictions agree on the direction of SOC change (net emissions or sequestration) across 46–51% of continental US land area. These differences are attributable to the high‐resolution site‐specific data used in the ML models compared to the relatively coarse grid represented in CMIP6 ESMs. Main conclusions Our high‐resolution estimates of baseline SOC stocks, identification of key environmental controllers, and projection of SOC changes from US land cover types under future climate scenarios suggest the need for high‐resolution simulations of SOC in ESMs to represent the heterogeneity of SOC. We found that the SOC change is sensitive to key soil related factors (e.g. soil drainage and soil order) that have not been historically considered as input parameters in ESMs, because currently more than 95% variability in the SOC of CMIP6 ESMs is controlled by net primary productivity, temperature, and precipitation. Using additional environmental factors to estimate the baseline SOC stocks and predict the future trajectory of SOC change can provide more accurate results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Bridging New Observational Capabilities and Process-Level Simulation: Insights into Aerosol Roles in the Earth System

The spatial distribution of ambient aerosol particles significantly impacts aerosol–radiation–cloud interactions, which contribute to the largest uncertainty in global anthropogenic radiative forcing estimations. However, the atmospheric boundary layer and lower free troposphere have not been adequately sampled in terms of spatiotemporal resolution, hindering a comprehensive characterization of various atmospheric processes and impeding our understanding of the Earth system. To address this research data gap, we have leveraged the development of uncrewed aerial systems (UAS) and advanced measurement techniques to obtain mesoscale spatial data on aerosol microphysical and optical properties around the U.S. Southern Great Plains (SGP) atmospheric observatory. Our study also benefits from state-of-the-art laboratory facilities that include three-dimensional molecular imaging techniques enabled by secondary ion mass spectrometry and nanogram-level chemical composition analysis via micronebulization aerosol mass spectrometry. Through our study, we have developed a framework for observation–modeling integration, enabling an examination of how various assumptions about the organic–inorganic components mixing state, inferred from chemical analysis, affect clouds and radiation in observation-constrained model simulations. By integrating observational constraints (derived from offline chemical analysis of the aerosol surface using collected samples) with in situ UAS observations, we have identified a prominent role of organic-enriched nanometer layers located at the surface of aerosol particles in determining profiles of aerosol optical and hygroscopic properties over the SGP observatory. Furthermore, we have improved the agreement between predicted clouds and ground-based cloud lidar measurements. This UAS–model–laboratory integration exemplifies how these new advanced capabilities can significantly enhance our understanding of aerosol–radiation–cloud interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Life-Cycle Inventory of Critical Materials: Nickel, Copper, Titanium, and Rare-Earth Elements

The United States has developed a list of critical minerals/materials whose sustained and reliable supply is pivotal to the robust functioning of critical industrial sectors. A key concern with these minerals is their environmental effects as a function of their production location. This requires material and energy flow details for their processing steps. This report provides a life-cycle inventory (LCI) for producing four critical minerals and/or material systems (nickel, copper, titanium, and rare-earth elements) incorporated in the updated GREET® (Greenhouse gases, Regulated Emissions, and Energy use in Technologies) model. For these systems, we provide an LCI as a function of production location – domestic (within the United States) and international (geographies from where the US imports these minerals). Our LCI also considers variations in ore grades for nickel and copper. This report also provides an LCI for all intermediate materials used to produce these critical materials.

58 GEOSCIENCES↗
DOE OSTI2025

The Capacity Expansion Regional Feasibility (CERF) Model: High-Resolution Power Plant Siting

Abstract This presentation gives an overview of the geospatial power plant siting model CERF. CERF (Capacity Expansion Regional Feasibility) is an open source Python package developed under the Integrated Multisector Multiscale Modeling (IM3) Project at PNNL. This presentation covers an overview of how the CERF model works, walks through various power plant siting analyses, and discusses future research opportunities for the model. The CERF model can be accessed at https://github.com/IMMM-SFA/cerf. PNNL Information Release Number: PNNL-SA-207336 Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

Mongird, Kendall [Pacific Northwest National Labor↗
DOE OSTI2021

Systemic effects of rising atmospheric vapor pressure deficit on plant physiology and productivity

Abstract Earth is currently undergoing a global increase in atmospheric vapor pressure deficit (VPD), a trend which is expected to continue as climate warms. This phenomenon has been associated with productivity decreases in ecosystems and yield penalties in crops, with these losses attributed to photosynthetic limitations arising from decreased stomatal conductance. Such VPD increases, however, have occurred over decades, which raises the possibility that stomatal acclimation to VPD plays an important role in determining plant productivity under high VPD. Furthermore, evidence points to more far‐ranging and complex effects of elevated VPD on plant physiology, extending to the anatomical, biochemical, and developmental levels, which could vary substantially across species. Because these complex effects are typically not considered in modeling frameworks, we conducted a quantitative literature review documenting temperature‐independent VPD effects on 112 species and 59 traits and physiological variables, in order to develop an integrated and mechanistic physiological framework. We found that VPD increase reduced yield and primary productivity, an effect that was partially mediated by stomatal acclimation, and also linked with changes in leaf anatomy, nutrient, and hormonal status. The productivity decrease was also associated with negative effects on reproductive development, and changes in architecture and growth rates that could decrease the evaporative surface or minimize embolism risk. Cross‐species quantitative relationships were found between levels of VPD increase and trait responses, and we found differences across plant groups, indicating that future VPD impacts will depend on community assembly and crop functional diversity. Our analysis confirms predictions arising from the hydraulic corollary to Darcy's law, outlines a systemic physiological framework of plant responses to rising VPD, and provides recommendations for future research to better understand and mitigate VPD‐mediated climate change effects on ecosystems and agro‐systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A global anthropogenic emission inventory of atmospheric pollutants from sector- and fuel-specific sources (1970–2017): an application of the Community Emissions Data System (CEDS)

Global anthropogenic emission inventories remain vital for understanding the sources of atmospheric pollution and the associated impacts on the environment, human health, and society. Rapid changes in today's society require that these inventories provide contemporary estimates of multiple atmospheric pollutants with both source sector and fuel type information to understand and effectively mitigate future impacts. To fill this need, we have updated the open-source Community Emissions Data System (CEDS) (Hoesly et al., 2019) to develop a new global emission inventory, CEDSGBD-MAPS. This inventory includes emissions of seven key atmospheric pollutants (NOx; CO; SO2; NH3; non-methane volatile organic compounds, NMVOCs; black carbon, BC; organic carbon, OC) over the time period from 1970–2017 and reports annual country-total emissions as a function of 11 anthropogenic sectors (agriculture; energy generation; industrial processes; on-road and non-road transportation; separate residential, commercial, and other sectors (RCO); waste; solvent use; and international shipping) and four fuel categories (total coal, solid biofuel, the sum of liquid-fuel and natural-gas combustion, and remaining process-level emissions). The CEDSGBD-MAPS inventory additionally includes monthly global gridded (0.5°?×?0.5°) emission fluxes for each compound, sector, and fuel type to facilitate their use in earth system models. CEDSGBD-MAPS utilizes updated activity data, updates to the core CEDS default scaling procedure, and modifications to the final procedures for emissions gridding and aggregation. Relative to the previous CEDS inventory (Hoesly et al., 2018), these updates extend the emission estimates from 2014 to 2017 and improve the overall agreement between CEDS and two widely used global bottom-up emission inventories. The CEDSGBD-MAPS inventory provides the most contemporary global emission estimates to date for these key atmospheric pollutants and is the first to provide global estimates for these species as a function of multiple fuel types and source sectors. Dominant sources of global NOx and SO2 emissions in 2017 include the combustion of oil, gas, and coal in the energy and industry sectors as well as on-road transportation and international shipping for NOx. Dominant sources of global CO emissions in 2017 include on-road transportation and residential biofuel combustion. Dominant global sources of carbonaceous aerosol in 2017 include residential biofuel combustion, on-road transportation (BC only), and emissions from the waste sector. Global emissions of NOx, SO2, CO, BC, and OC all peak in 2012 or earlier, with more recent emission reductions driven by large changes in emissions from China, North America, and Europe. In contrast, global emissions of NH3 and NMVOCs continuously increase between 1970 and 2017, with agriculture as a major source of global NH3 emissions and solvent use, energy, residential, and the on-road transport sectors as major sources of global NMVOCs. Due to similar development methods and underlying datasets, the CEDSGBD-MAPS emissions are expected to have consistent sources of uncertainty as other bottom-up inventories. The CEDSGBD-MAPS source code is publicly available online through GitHub: https://github.com/emcduffie/CEDS/tree/CEDS_GBD-MAPS (last access: 1 December 2020). The CEDSGBD-MAPS emission inventory dataset (both annual country-total and monthly global gridded files) is publicly available under https://doi.org/10.5281/zenodo.3754964 (McDuffie et al., 2020c).

Mcduffie, Erin↗
DOE OSTI2023

Comparing the Performance of Julia on CPUs versus GPUs and Julia-MPI versus Fortran-MPI: a case study with MPAS-Ocean (Version 7.1)

Abstract. Some programming languages are easy to develop at the cost of slow execution, while others are fast at runtime but much more difficult to write. Julia is a programming language that aims to be the best of both worlds – a development and production language at the same time. To test Julia's utility in scientific high-performance computing (HPC), we built an unstructured-mesh shallow water model in Julia and compared it against an established Fortran-MPI ocean model, the Model for Prediction Across Scales–Ocean (MPAS-Ocean), as well as a Python shallow water code. Three versions of the Julia shallow water code were created: for single-core CPU, graphics processing unit (GPU), and Message Passing Interface (MPI) CPU clusters. Comparing identical simulations revealed that our first version of the Julia model was 13 times faster than Python using NumPy, where both used an unthreaded single-core CPU. Further Julia optimizations, including static typing and removing implicit memory allocations, provided an additional 10–20× speed-up of the single-core CPU Julia model. The GPU-accelerated Julia code was almost identical in terms of performance to the MPI parallelized code on 64 processes, an unexpected result for such different architectures. Parallelized Julia-MPI performance was identical to Fortran-MPI MPAS-Ocean for low processor counts and ranges from 2× faster to 2× slower for higher processor counts. Our experience is that Julia development is fast and convenient for prototyping but that Julia requires further investment and expertise to be competitive with compiled codes. We provide advice on Julia code optimization for HPC systems.

54 ENVIRONMENTAL SCIENCES↗