Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SYSTEMS MODELING”

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 649 records · Page 36

Greenland Surface Mass Balance as Simulated by the Community Earth System Model. Part II: Twenty-First-Century Changes

This study presents the first twenty-first-century projections of surface mass balance (SMB) changes for the Greenland Ice Sheet (GIS) with the Community Earth System Model (CESM), which includes a new ice sheet component. For glaciated surfaces, CESM includes a sophisticated calculation of energy fluxes, surface albedo, and snowpack hydrology (melt, percolation, refreezing, etc.). To efficiently resolve the high SMB gradients at the ice sheet margins and provide surface forcing at the scale needed by ice sheet models, the SMB is calculated at multiple elevations and interpolated to a finer 5-km ice sheet grid. During a twenty-first-century simulation driven by representative concentration pathway 8.5 (RCP8.5) forcing, the SMB decreases from 372 ± 100 Gt yr -1 in 1980–99 to -78 ± 143 Gt yr -1 in 2080–99. The 2080–99 near-surface temperatures over the GIS increase by 4.7 K (annual mean) with respect to 1980–99, only 1.3 times the global increase (+3.7 K). Snowfall increases by 18%, while surface melt doubles. The ablation area increases from 9% of the GIS in 1980–99 to 28% in 2080–99. Over the ablation areas, summer downward longwave radiation and turbulent fluxes increase, while incoming shortwave radiation decreases owing to increased cloud cover. The reduction in GIS-averaged July albedo from 0.78 in 1980–99 to 0.75 in 2080–99 increases the absorbed solar radiation in this month by 12%. Summer warming is strongest in the north and east of Greenland owing to reduced sea ice cover. In the ablation area, summer temperature increases are smaller due to frequent periods of surface melt.

54 ENVIRONMENTAL SCIENCES↗

Combining Earth System Modeling and Machine Learning to Investigate Volcanic Sulfate Deposition in Polar Ice Cores

Volcanic eruptions emit large amounts of sulfur dioxide (SO2), water, and other chemicals into the atmosphere, both in the troposphere and the stratosphere. Most of the SO2 is converted to sulfate aerosol, which is eventually deposited following long-range transport. The deposits from large eruptions are potentially detectable in ice cores, but there are many cases in which sulfate layers have not been linked to their source volcanoes. As volcanoes can act as significant shocks to the global climate system, we are interested in locating these eruptions in order to increase understanding of the volcanic record. To narrow down the search, we performed 140 simulations of volcanic eruptions using the GISS ModelE Earth system model. We varied the latitude, longitude, Julian day, plume top, plume bottom, and injected SO2 and H2O amounts using a Latin hypercube sampling approach, and analyzed correlations between these parameters and sulfate depositions at ice core sites in Antarctica and Greenland. Using machine learning and parameter estimation, we generated probability distributions and maximum likelihood estimates for the parameters given sulfate deposition data, which can predict latitude with some skill. We find that the volcano latitude and SO2 content are best correlated with sulfate depositions at each pole, while longitude, Julian day, and H2O have small or insignificant effects. Plume altitude and thickness are important because they determine how much of the SO2 is injected into the stratosphere, which has implications for sulfur transport and lifetimes.

Earth system models↗

A Macroalgal Cultivation Modeling System (MACMODS): Evaluating the Role of Physical-Biological Coupling on Nutrients and Farm Yield

Offshore aquaculture has the potential to expand the macroalgal industry. However, moving into deeper waters requires suspended structures that will present novel farm-environment interactions. Here, we present a computational modeling framework, the Macroalgal Cultivation Modeling System (MACMODS), to explore within-farm modifications to light, seawater flow, and nutrient fields across time and space scales relevant to macroalgae. A regional ocean model informs the site-specific setting, the Santa Barbara Channel in the Southern California Bight. A fine-scale hydrodynamic model predicts modified flows and turbulent mixing within the farm. A spatially resolved macroalgal growth model, parameterized for giant kelp, Macrocystis pyrifera , predicts kelp biomass. Key findings from model integration are that regional ocean conditions set overall farm performance, while fine-scale within-farm circulation and nutrient delivery are important to resolve variation in within-farm macroalgal performance. Therefore, we conclude that models resolving within-farm dynamics can provide benefit to farmers with insight on how farm design and regional ocean conditions interact to influence overall yield. Here, the presence of repeating longlines aligned with the mean current generate flow diversions around the farm as well as attached Langmuir circulations and increased turbulence intensity. These flow-induced phenomena lead to less biomass in the interior portion of the farm relative to the edges. We also find that there is an effluent “footprint” that extends as much as 20 km beyond the farm. In this regard, MACMODS can be used to not only evaluate farm design and cultivation practices that maximize yield but also explore interactions between the farm and ecosystem in order to minimize impacts.

54 ENVIRONMENTAL SCIENCES↗

Binder model system to be used for determination of prepolymer functionality

Development of a method for determining the functionality distribution of prepolymers used for rocket binders is discussed. Research has been concerned with accurately determining the gel point of a model polyester system containing a single trifunctional crosslinker, and the application of these methods to more complicated model systems containing a second trifunctional crosslinker, monofunctional ingredients, or a higher functionality crosslinker. Correlations of observed with theoretical gel points for these systems would allow the methods to be applied directly to prepolymers.

Martinelli, F. J.↗

Climate, Variability, and Climate Sensitivity of “Middle Atmosphere” Chemistry Configurations of the Community Earth System Model Version 2, Whole Atmosphere Community Climate Model Version 6 (CESM2(WACCM6))

Abstract Simulating whole atmosphere dynamics, chemistry, and physics is computationally expensive. It can require high vertical resolution throughout the middle and upper atmosphere, as well as a comprehensive chemistry and aerosol scheme coupled to radiation physics. An unintentional outcome of the development of one of the most sophisticated and hence computationally expensive model configurations is that it often excludes a broad community of users with limited computational resources. Here, we analyze two configurations of the Community Earth System Model Version 2, Whole Atmosphere Community Climate Model Version 6 (CESM2(WACCM6)) with simplified “middle atmosphere” chemistry at nominal 1 and 2° horizontal resolutions. Using observations, a reanalysis, and direct model comparisons, we find that these configurations generally reproduce the climate, variability, and climate sensitivity of the 1° nominal horizontal resolution configuration with comprehensive chemistry. While the background stratospheric aerosol optical depth is elevated in the middle atmosphere configurations as compared to the comprehensive chemistry configuration, it is comparable among all configurations during volcanic eruptions. For any purposes other than those needing an accurate representation of tropospheric organic chemistry and secondary organic aerosols, these simplified chemistry configurations deliver reliable simulations of the whole atmosphere that require 35% and 86% fewer computational resources at nominal 1 and 2° horizontal resolution, respectively.

54 ENVIRONMENTAL SCIENCES↗

Reducing Long‐Standing Surface Ozone Overestimation in Earth System Modeling by High‐Resolution Simulation and Dry Deposition Improvement

The overestimation of surface ozone concentration in low‐resolution global atmospheric chemistry and climate models has been a long‐standing issue. We first update the ozone dry deposition scheme in both high‐ (0.25°) and low‐resolution (1°) Community Earth System Model (CESM) version 1.3 runs, by adding the effects of leaf area index and correcting the sunlit and shaded fractions of stomatal resistances. With this update, 5‐year‐long summer simulations (2015–2019) using the low‐resolution CESM still exhibit substantial ozone overestimation (by 6.0–16.2 ppbv) over the U.S., Europe, eastern China, and ozone pollution hotspots. The ozone dry deposition scheme is further improved by adjusting the leaf cuticle conductance, reducing the mean ozone bias by 19%, and increasing the model resolution further reduces the ozone overestimation by 43%. We elucidate the mechanism by which model grid spacing influences simulated ozone, revealing distinctive pathways in urban versus rural areas. In rural areas, grid spacing mainly affects daytime ozone levels, where additional NO x emissions from nearby urban areas result in an ozone boost and overestimation in low‐resolution simulations. In contrast, over urban areas, daytime ozone overestimation follows a similar mechanism due to the influence of volatile organic compounds from surrounding rural areas. However, nighttime ozone overestimation is closely linked to weakened NO titration owing to the redistribution of urban NO x to rural areas. Additionally, stratosphere‐troposphere exchange may also contribute to reducing ozone bias in high‐resolution simulations, warranting further investigation. This optimized high‐resolution CESM may enhance understanding of ozone formation mechanisms, sources, and changes in a warming climate.

54 ENVIRONMENTAL SCIENCES↗

Forced Patterns of Sea Level Rise in the Community Earth System Model Large Ensemble From 1920 to 2100

To provide context for observed sea level rise, the forced response (FR) in dynamic sea level (DSL) during the 20th and 21st centuries is examined in the Community Earth System Model Large Ensemble (LE). After accounting for simulation drift, which in the LE is sizable, the DSL FR is found to be complex, both in space and time. Its evolving character is suggested to arise from both the diversity and transient evolution of climate forcing agents and the slow adjustment timescales of the intermediate and deep oceans. Nonetheless, various intervals of spatially coherent change simulated for the recent past and near future are identified, and their characteristics and associated driving mechanisms are identified and discussed. The mid- to late-20th century DSL FR is characterized by an hemispherically asymmetric pattern of change, with depressed rates of rise in the northern oceans. There is also a dipole of change in the Southern Ocean caused by changes in near-surface zonal winds. Through the late 20th and early 21st centuries, a different pattern of rise emerges with elevated rates in the tropics and depressed rates at high latitudes. Zonal and interbasin variations characterize both intervals and involve the pattern and depth of anomalous ocean heat content storage and spatial contrasts in the expansion coefficient tied mainly to base state temperature, with greater rates of rise in warm regions per unit warming. Finally, the relative roles of surface flux and ocean convergence anomalies are examined.

54 ENVIRONMENTAL SCIENCES↗

Thermodynamic versus Dynamic Controls on Extreme Precipitation in a Warming Climate from the Community Earth System Model Large Ensemble

The moisture budget is evaluated as a function of the probability distribution of precipitation for the end of the twentieth century and projected end of the twenty-first century in the Community Earth System Model Large Ensemble. For a given precipitation percentile, a conditional moisture budget equation relates precipitation minus evaporation ( P − E) to vertical moisture transport, horizontal moisture advection, and moisture storage. At high percentiles, moisture advection and moisture storage cancel and evaporation is negligible, so that precipitation is approximately equal to vertical moisture transport, and likewise for projected changes. Therefore, projected changes to extreme precipitation are approximately equal to the sum of thermodynamic and dynamic tendencies, representing changes to the vertical profiles of moisture content and mass convergence, respectively. The thermodynamic tendency is uniform across percentiles and regions as an intensification of the hydrological cycle, but the dynamic tendency is more complex. For extreme events, per degree of warming, in the mid-to-high latitudes the dynamic tendency is small, so that precipitation approximately scales by the Clausius–Clapeyron 7% K −1 increase. In the subtropics, a drying tendency originating from dynamics offsets the thermodynamic wetting tendency, with the net effect on precipitation varying among regions. The effect of this dynamic drying decreases with increasing percentile. In the deep tropics, a positive dynamic tendency occurs with magnitude similar to or greater than the positive thermodynamic tendency, resulting in generally a 10%–15% K −1 precipitation increase, and with a >25% K −1 increase over the tropical east Pacific. This reinforcing dynamical tendency increases rapidly for high percentiles.

54 ENVIRONMENTAL SCIENCES↗

Relative Controls of Vapor Pressure Deficit and Soil Water Stress on Canopy Conductance in Global Simulations by an Earth System Model

Climate change impacts drought and vapor pressure deficit (VPD) that critically influence land-atmosphere water vapor and CO 2 exchange. There is a lack of model-based global evaluation of the relative roles of the two types of water stresses in limiting canopy conductance (G c ). Using the Energy Exascale Earth System Model (E3SM), we conducted four land-atmosphere coupled simulations for the historical period with and without a plant hydraulics scheme (PHS) and under present-day and quadrupling of CO 2 concentrations. By fitting an empirical model of G c to the E3SM simulations, we evaluated the sensitivity of G c to water stresses caused by VPD and soil water (SW) in the growing season. The empirical model is based on a multiplicative algorithm that adjusts G c according to VPD, soil water saturation, and their interaction. Our results show positive trends of VPD over land during the historical period, with large differences between simulations with and without PHS. Despite the differences, a majority of land pixels that experience both SW and VPD stresses are in semi-arid and tropical savanna ecosystems. G c is more sensitive to VPD than SW in less than 20% of land pixels where the empirical model provides a good fit for the simulated G c . Ignoring SW and VPD interaction may underestimate the relative role of VPD on G c by as much as 28%. Including PHS accentuates SW stress limitation on G c by simulating wetter soil and higher ET that lower VPD. Elevated CO 2 reduces G c by increasing water use efficiency and ameliorates SW stress on G c .

54 ENVIRONMENTAL SCIENCES↗

Response of Simple, Model Systems to Extreme Conditions

The focus of the research was on the application of high-pressure/high-temperature techniques, together with intense energetic ion beams, to the study of the behavior of simple oxide systems (e.g., SiO 2 , GeO 2 , CeO 2 , TiO 2 , HfO 2 , SnO 2 , ZnO and ZrO 2 ) under extreme conditions. These simple stoichiometries provide unique model systems for the analysis of structural responses to pressure up to and above 1 Mbar, temperatures of up to several thousands of kelvin, and the extreme energy density generated by energetic heavy ions (tens of keV/atom). The investigations included systematic studies of radiation- and pressure-induced amorphization of high P-T polymorphs. By studying the response of simple stoichiometries that have multiple structural “outcomes”, we have established the basic knowledge required for the prediction of the response of more complex structures to extreme conditions. We especially focused on the amorphous state and characterized the different non-crystalline structure-types that result from the interplay of radiation and pressure. For such experiments, we made use of recent technological developments, such as the perforated diamond-anvil cell and in situ investigation using synchrotron x-ray sources. We have been particularly interested in using extreme pressures to alter the electronic structure of a solid prior to irradiation. We expected that the effects of modified band structure would be evident in the track structure and morphology, information which is much needed to describe theoretically the fundamental physics of track-formation. Finally, we investigated the behavior of different simple-oxide, composite nanomaterials (e.g., uncoated nanoparticles vs. core/shell systems) under coupled, extreme conditions. This provided insight into surface and boundary effects on phase stability under extreme conditions.

36 MATERIALS SCIENCE↗

Photovoltaic System Modeling: Uncertainty and Sensitivity Analyses

We report an uncertainty and sensitivity analysis for modeling AC energy from photovoltaic systems. Output from a PV system is predicted by a sequence of models. We quantify uncertainty in the output of each model using empirical distributions of each model’s residuals. We propagate uncertainty through the sequence of models by sampling these distributions to obtain an empirical distribution of a PV system’s output. We consider models that: (1) translate measured global horizontal, direct and global diffuse irradiance to plane-of-array irradiance; (2) estimate effective irradiance; (3) predict cell temperature; (4) estimate DC voltage, current and power; (5) reduce DC power for losses due to inefficient maximum power point tracking or mismatch among modules; and (6) convert DC to AC power. Our analysis considers a notional PV system comprising an array of FirstSolar FS-387 modules and a 250 kW AC inverter; we use measured irradiance and weather at Albuquerque, NM. We found the uncertainty in PV system output to be relatively small, on the order of 1% for daily energy. We found that uncertainty in the models for POA irradiance and effective irradiance to be the dominant contributors to uncertainty in predicted daily energy. Our analysis indicates that efforts to reduce the uncertainty in PV system output predictions may yield the greatest improvements by focusing on the POA and effective irradiance models.

14 SOLAR ENERGY↗

Workshop on the Initialization of High-Resolution Earth System Models

From weather to climate timescales, weeks to decades, societal needs require detailed predictions of the atmosphere, ocean, and land-surface to facilitate informed planning. Central to this scientific endeavor is the ability to accurately initialize the earth system models (ESMs), since frequently the skill in forecasts, predictions, and projections depends on the quality of the initialization. However, initialization of high-resolution ESMs remains an outstanding challenge, particularly for the climate modeling community where initialization has received relatively less attention and where longer-term information is needed. This need, after a review of previous workshop outcomes, led the U.S. Department of Energy (DOE) and National Oceanic and Atmospheric Administration (NOAA) to host a workshop that focused exclusively on the opportunities that data assimilation (DA) and, specifically, coupled data assimilation (CDA) might provide to the ESM enterprise with particular attention toward initialization. The workshop was held April 9-10, 2018, at the Rockville Hilton in Rockville, Maryland, with 41 participants. It was composed of seven sessions, each with a distinct topical theme that related to the broad challenge of initializing high-resolution ESMs. The first session was used to frame the workshop challenge. Two of the sessions (2 and 7) focused the current state of the art of CDA within laboratory and operational centers at both the short-term, subseasonal and long-term, climate time scales. Three of the sessions (3, 4, and 5) addressed new algorithmic approaches to CDA at high resolution. Session 6 focused on the computational science solutions for the management of large input/output streams required to carry out the assimilation procedure.

54 ENVIRONMENTAL SCIENCES↗

STITCHES: a Python package to amalgamate existing Earth system model output into new scenario realizations

Understanding the interaction between humans and the Earth system is a computationally daunting task, with many possible approaches depending on resources available and questions of interest. For example, state-of-the-art impact models require decade-long time series of relatively high frequency, spatially resolved and often multiple variables representing climatic impact-drivers (Ruane et al., 2022). Most commonly these are derived from the outputs of detailed, computationally expensive Earth System Models (ESMs) run according to a standard, limited set of future scenarios, the latest being the SSP-RCPs run under CMIP6/ScenarioMIP (Eyring et al., 2016; O’Neill et al., 2016). At the time of writing, O’Neill et al. (2016) has been cited more than 1750 times and Eyring et al. (2016) more than 5000 times, highlighting the broad, general applications of this data. Often, however, impact modeling seeks to explore new scenarios that were not part of the ScenarioMIP protocol, and/or needs a larger set of initial condition ensemble members than are typically available to quantify the effects of ESM internal variability. In addition, the recognition that the human and Earth systems are fundamentally intertwined, and may feature potentially significant feedback loops, is making integrated, simultaneous modeling of the coupled human-Earth system increasingly necessary, if computationally challenging with most existing tools (Thornton et al., 2017). For impact modelers, climate model emulators can be the answer to meet both the needs of: 1) creating realizations for novel scenarios and 2) achieving a simplified, computationally tractable representation of ESM behavior in a coupled human-Earth system modeling framework. We proposed a new, comprehensive approach to such emulation of gridded, multivariate ESM outputs for novel scenarios without the computational cost of a full ESM, STITCHES (Tebaldi et al., 2022). The approach outlined in Tebaldi et al. (2022) should be extensible to future CMIP eras, although the STITCHES software at present is strictly focused on CMIP6/ScenarioMIP data hosted on Pangeo (https://gallery.pangeo.io/repos/pangeo-gallery/cmip6/). The corresponding STITCHES Python package uses existing archives of ESMs’ scenario experiments from CMIP6/ScenarioMIP to construct gridded, multivariate realizations of new scenarios provided by reduced complexity climate models (Hartin et al., 2015; Meinshausen et al., 2011; Smith et al., 2018), or to enrich existing initial condition ensembles. Its output provides the same characteristics as the emulated ESM output: multivariate (spanning potentially all variables that the ESM has saved), spatially resolved (down to the native grid of the ESM), and preserving the same high frequency as the original data. A new realization of multiple variables can be generated on the order of minutes with STITCHES, rather than the hours or sometimes days that ESMs require.

97 MATHEMATICS AND COMPUTING↗

Important ice processes are missed by the Community Earth System Model (CESM) in Southern Ocean mixed–phase clouds: Bridging SOCRATES observations to model developments

Global climate models (GCMs) are challenged by difficulties in simulating cloud phase and cloud radiative effect over the Southern Ocean (SO). Some of the new generation GCMs predict too much liquid and too little ice in mixed-phase clouds. This misrepresentation of cloud phase in GCMs results in weaker negative cloud feedback over the SO and a higher climate sensitivity. Based on a model comparison with observational data obtained during the Southern Ocean Cloud Radiation and Aerosol Transport Experimental Study (SOCRATES), this study addresses a key uncertainty in the Community Earth System Model version 2 (CESM2) related to cloud phase, namely ice formation in pristine remote SO clouds. It is found that sea spray organic aerosols (SSOA) are the most important type of ice nucleating particles (INPs) over the SO with concentrations one order of magnitude higher than those of dust INPs based on measurements and CESM2 simulations. Secondary ice production (SIP) which includes riming splintering, rain droplet shattering, and ice-ice collisional fragmentation as implemented in CESM2 is the dominant ice production process in moderately cold clouds with cloud temperatures greater than –20°C. SIP enhances the in-cloud ice number concentrations (Ni) by 1-3 orders of magnitude and predicts more mixed-phase (with percentage occurrence increased from 15% to 21%), in better agreement with the observations. Finally, this study highlights the importance of accurately representing cloud phase over the pristine remote SO by considering the ice nucleation of SSOA and SIP processes, which are currently missing in most GCM cloud microphysics parameterizations.

54 ENVIRONMENTAL SCIENCES↗

Climate response to off-equatorial stratospheric sulfur injections in three Earth system models – Part 2: Stratospheric and free-tropospheric response

The paper constitutes Part 2 of a study performing a first systematic inter-model comparison of the atmospheric responses to stratospheric aerosol injection (SAI) at various single latitudes in the tropics, as simulated by three state-of-the-art Earth system models – CESM2-WACCM6, UKESM1.0, and GISS-E2.1-G. Building on Part 1 (Visioni et al., 2023) we demonstrate the role of biases in the climatological circulation and specific aspects of the model microphysics in driving the inter-model differences in the simulated sulfate distributions. We then characterize the simulated changes in stratospheric and free-tropospheric temperatures, ozone, water vapor, and large-scale circulation, elucidating the role of the above aspects in the surface SAI responses discussed in Part 1. We show that the differences in the aerosol spatial distribution can be explained by the significantly faster shallow branches of the Brewer–Dobson circulation in CESM2, a relatively isolated tropical pipe and older tropical age of air in UKESM, and smaller aerosol sizes and relatively stronger horizontal mixing (thus very young stratospheric age of air) in the two GISS versions used. We also find a large spread in the magnitudes of the tropical lower-stratospheric warming amongst the models, driven by microphysical, chemical, and dynamical differences. These lead to large differences in stratospheric water vapor responses, with significant increases in stratospheric water vapor under SAI in CESM2 and GISS that were largely not reproduced in UKESM. For ozone, good agreement was found in the tropical stratosphere amongst the models with more complex microphysics, with lower stratospheric ozone changes consistent with the SAI-induced modulation of the large-scale circulation and the resulting changes in transport. In contrast, we find a large inter-model spread in the Antarctic ozone responses that can largely be explained by the differences in the simulated latitudinal distributions of aerosols as well as the degree of implementation of heterogeneous halogen chemistry on sulfate in the models. The use of GISS runs with bulk microphysics demonstrates the importance of more detailed treatment of aerosol processes, with contrastingly different stratospheric SAI responses to the models using the two-moment aerosol treatment; however, some problems in halogen chemistry in GISS are also identified that require further attention. Overall, our results contribute to an increased understanding of the underlying physical mechanisms as well as identifying and narrowing the uncertainty in model projections of climate impacts from SAI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluating the WRF-Hydro Modeling System in Alaska

National Water Model (NWM) implemented operationally in August 2016 to improve hydrological prediction (OWP, 2017). (1) Four operational configurations (2) Only covers contiguous United States (US). NWM is instantiation of Weather Research and Forecasting model hydrological extension package (WRF-Hydro)(Gochis et al., 2013) coupled with Noah Land Surface Model with Multi-Parameterization options (Noah-MP)(Niu et al., 2011). WRF-Hydro is extensible, high-resolution hydrologic routing and streamflow modeling framework, coupling column land surface, terrain routing, and channel routing modules (NCAR, 2017). This project uses experimental version of WRF-Hydro in Alaska mimicking the NWM to: (1) Identify modeling challenges for NWM development in Alaska (2) Assess WRF-Hydro and NWM ability to represent unique hydrological processes of arctic regions and accurately predict high and low flow events (3) Examine impacts of assimilating Surface Water Ocean Topography (SWOT) (Biancamaria et al., 2016) observations to improve model initialization.

Hydrology↗