Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesoscale science”

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 811 records · Page 45

Kinetics and Thermodynamics of Sr Permeation in CeO 2 -Based Barrier Layers for Solid-Oxide Electrolyzer Cells

Solid-oxide electrolyzer cells (SOECs) convert steam to hydrogen efficiently at high temperatures. However, during operation, the diffusion of cations or impurities through the cells due to electrode degradation can cause unwanted secondary phases to form, which may degrade device performance. Here, in this study, we use atomistic and mesoscale simulations coupled with experimental analysis to study the diffusion of Sr through the Gd-doped CeO 2 (GDC) barrier layer used to protect the yttria-stabilized zirconia (YSZ) electrolyte in SOECs. From our atomistic calculations, we find Sr diffusion to be negligibly slow in bulk GDC; however, surface diffusion is much more favorable. Subsequent mesoscale simulations show that Sr diffusion is activated when the porosity of GDC exceeds ∼10% and significantly exceeds diffusion in bulk and grain boundary regions. We also find that SrO-based species can accumulate at GDC surfaces; however, SrO aggregation and coarsening will be limited by the large lattice mismatch between GDC and SrO. Energy-dispersive X-ray spectroscopy (EDS) and electron diffraction confirm that Sr can accumulate within GDC pores and form disperse Sr-containing secondary phases. Altogether, Sr diffusion in dense GDC is unlikely to give rise to thick SrO layers, which would severely limit device performance. The formation of Sr-containing secondary phases can largely be avoided by restricting the porosity of the GDC layer as much as possible.

36 MATERIALS SCIENCE↗

The Impacts of Rotational Mixing on the Precipitation Simulated by a Convection Permitting Model

With increased availability of computational resources, regional and global scale convection-permitting model (CPM, Δx ~ 1–10 km) simulations are becoming more common. CPMs have improved accuracy in their representation of deep convection and mesoscale convective systems (MCSs) compared to coarser resolution models. However, CPMs still exhibit convective cloud and precipitation biases relative to observations, notably a lesser frequency of light precipitation rates and greater frequency of heavy precipitation rates. In this work we hypothesize that these CPM biases are related to under-resolved mixing between convective updrafts and their surrounding environment. To test this hypothesis, we introduce a parameterization to the Weather Research and Forecasting model (WRF) that adds a small angular rotation of the grid-scale flow about the axis perpendicular to the plane of convective drafts. This rotated flow is then allowed to alter advection of moisture and hydrometeors. The effects of such mixing on precipitation characteristics are evaluated in month-long 4-km grid spacing simulations over the Amazon. The enhanced mixing transports moisture and condensate from convective cores to other areas including downdrafts. This increases the frequency of low-precipitable water and light precipitation. It also decreases the frequency of intense precipitation from isolated deep convection and MCSs, increases cloud top temperatures, reduces radar echo-top heights, and increases overall precipitation by altering the relationship of precipitation with precipitable water, in better agreement with observations. The results suggest when optimized using multiple observations, such an approach may provide a path toward more accurate representation of convection and precipitation statistics in convection-permitting simulations.

54 ENVIRONMENTAL SCIENCES↗

Lattice disorder effect on magnetic ordering of iron arsenides

This study investigates magnetic ordering temperature in nano- and mesoscale structural features in an iron arsenide. Although magnetic ground states in quantum materials can be theoretically predicted from known crystal structures and chemical compositions, the ordering temperature is harder to pinpoint due to potential local lattice variations that calculations may not account for. In this work we find surprisingly that a locally disordered material can exhibit a significantly larger Néel temperature (T N ) than an ordered material of precisely the same chemical stoichiometry. Here, a EuFe 2 As 2 crystal, which is a ‘122’ parent of iron arsenide superconductors, is found through synthesis to have ordering below T N = 195 K (for the locally disordered crystal) or T N = 175 K (for the ordered crystal). In the higher T N crystals, there are shorter planar Fe-Fe bonds [2.7692(2) Å vs. 2.7745(3) Å], a randomized in-plane defect structure, and diffuse scattering along the [00 L] crystallographic direction that manifests as a rather broad specific heat peak. For the lower T N crystals, the a-lattice parameter is larger and the in-plane microscopic structure shows defect ordering along the antiphase boundaries, giving a larger T N and a higher superconducting temperature (T c ) upon the application of pressure. First-principles calculations find a strong interaction between c-axis strain and interlayer magnetic coupling, but little impact of planar strain on the magnetic order. Neutron single-crystal diffraction shows that the low-temperature magnetic phase transition due to localized Eu moments is not lattice or disorder sensitive, unlike the higher-temperature Fe sublattice ordering. This study demonstrates a higher magnetic ordering point arising from local disorder in 122.

36 MATERIALS SCIENCE↗

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗

Exploring the Meteorological Impacts of Surface and Rooftop Heat Mitigation Strategies Over a Tropical City

Different heat mitigation technologies have been developed to improve the thermal environment in cities. However, the regional impacts of such technologies, especially in the context of a tropical city, remain unclear. The deployment of heat mitigation technologies at city-scale can change the radiation balance, advective flow, and energy balance between urban areas and the overlying atmosphere. We used the mesoscale Weather Research and Forecasting model coupled with a physically based single-layer urban canopy model to assess the impacts of five different heat mitigation technologies on surface energy balance, standard surface meteorological fields, and planetary boundary layer (PBL) dynamics for premonsoon typical hot summer days over a tropical coastal city in the month of April in 2018, 2019, and 2020. Results indicate that the regional impacts of cool materials (CMs), super-cool broadband radiative coolers, green roofs (GRs), vegetation fraction change, and a combination of CMs and GRs (i.e., “Cool city (CC)”) on the lower atmosphere are different at diurnal scale. Results showed that super-cool materials have the maximum potential of ambient temperature reduction of 1.6°C during peak hour (14:00 LT) compared to other technologies in the study. During the daytime hours, the PBL height was considerably lower than the reference scenario with no implementation of strategies by 700 m for super-cool materials and 500 m for both CMs and CC cases; however, the green roofing system underwent nominal changes over the urban area. During the nighttime hours, the PBL height increased by CMs and the CC strategies compared to the reference scenario, but minimal changes were evident for super-cool materials. The changes of temperature on the vertical profile of the heat mitigation implemented city reveal a stable PBL over the urban domain and a reduction of the vertical mixing associated with a pollution dome. This would lead to crossover phenomena above the PBL due to the decrease in vertical wind speed. Therefore, assessing the coupled regional impact of urban heat mitigation over the lower atmosphere at city-scale is urgent for sustainable urban planning.

54 ENVIRONMENTAL SCIENCES↗

Uncertainty Quantification Framework for Predicting Material Response with Large Number of Parameters: Application to Creep Prediction in Ferritic-Martensitic Steels Using Combined Crystal Plasticity and Grain Boundary Models

This paper presents an uncertainty quantification (UQ) framework for the physics-based model prediction of material response with a large number of parameters. The application problem presented in this work is that of predicting creep in Grade 91 steel at 600°C. The material response is defined with a physically based microstructural model with constitutive equations emulating several observed phenomena in Grade 91 and embodied into an explicit geometry mesoscale finite element model for prior austenite grains and grain boundaries. Creep within the grains and in grain boundaries are represented by crystal plasticity for dislocation motion and a physics-based model for cavity growth and nucleation, respectively. The creep behavior of this material is influenced by several parameters, some of which have a wide range of variation based on experimental data. UQ combined with microstructural modeling can discover the core microstructural causes of experimental variability, leading to improved materials with lower variability in critical long-term material properties. In this study, we investigate the model's uncertainty to identify material properties that may be modified during production to increase creep life and analyze different components of the crystal plasticity model for improvements. For this purpose, a quantity of interest is defined as time to minimum creep rate, which correlates well to the creep failure of the material. A deep neural network model was trained and validated to be used as a surrogate for the finite element model. Then, a variance-based sensitivity analysis is performed on the surrogate model to find the Sobol indices of the input parameters in respect to the output quantity of interest. The Sobol indices are used to reduce the dimensionality of the model. Generalized polynomial chaos expansion is used on the reduced basis models to propagate the uncertainty from the input parameters to the quantity of interest using the deep neural network surrogate model. These results are benchmarked against uncertainty propagation using Monte Carlo simulations. In conclusion, the UQ performed through the reduced basis model captures almost all the uncertainty in the model with significantly fewer simulations, making it possible to perform the UQ directly via simulations with the finite element model rather than surrogate machine-learned models.

36 MATERIALS SCIENCE↗

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗

Segmentation of RDX and TNT in X‐Ray Computed Tomography Reconstructions of Melt‐Cast Explosives

ABSTRACT Three‐dimensional mesoscale characterization of heterogeneous melt‐cast high explosives is challenging because of the difficulty differentiating binder from explosive crystals: two functionally different materials which are typically similar in density by design. Here, we report an algorithm which can differentiate hexahydro‐1,3,5‐trinitro‐1,3,5‐triazine (RDX) from 2,4,6‐trinitrotoluene (TNT) in x‐ray computed tomography (CT) volumes with tens of microns resolution. This method allows us to quantify RDX/TNT content, porosity, and RDX domain size. We calibrated the segmentation algorithm using simulated x‐ray CT volumes containing object models of RDX crystals within a TNT matrix. We then segmented and analyzed CT data for Composition B (Comp B), a 60/40 RDX/TNT mixture, and Cyclotol, a 75/25 RDX/TNT mixture. We examined melt‐cast samples fabricated with 100% theoretical maximum density (TMD) and 85% TMD. For the 100% TMD Comp B and Cyclotol samples, the RDX content values calculated by segmentation were 3% and 9% lower, respectively, than the values measured by high‐performance liquid chromatography on material from the same synthesis lots. This result is consistent with the expected underreporting of RDX content resulting from x‐ray CT resolution limits on RDX particles with diameters smaller than 25 µm. The 85% TMD samples were less accurately segmented with our algorithm due to the confounding presence of voids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mesoscale Cellular Convection Detection and Classification Using Convolutional Neural Networks: Insights From Long-Term Observations at ARM Eastern North Atlantic Site

Marine boundary layer clouds are crucial in Earth's climate system. They frequently manifest as closed or open cell mesoscale cellular convection (MCC). MCC clouds are challenging to represent accurately in current climate models, highlighting the need for detailed observational data sets and in-depth analyses. This study utilizes over 8 years of observations from the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) User Facility Eastern North Atlantic (ENA) site at Graciosa Island, Azores, to investigate these clouds. We first apply a convolutional neural network with a U-Net architecture to classify open and closed cells, marking the first application of such an approach for automatically detecting MCC patterns from ground-based radar measurements. This method addresses some observational gaps in satellite data related to low temporal resolution, nighttime challenges, and limited vertical structure capture. The analysis of the MCC cases shows clear differences between closed and open MCCs: Closed MCC clouds are characterized by lower cloud tops and bases, shallower cloud geometrical depth, weaker horizontal wind speeds, stronger atmospheric stability, and a more homogeneous liquid water path than open MCCs. Finally, we demonstrate two potential applications of our radar-based MCC classifications: (a) facilitating the investigation of aerosol-cloud interactions and (b) exploring meteorological factors along with MCC's evolution by integrating satellite imagery and back-trajectory analysis. The identified MCC cases offer a valuable resource for the scientific community to study MCC processes further and improve climate model accuracy.

54 ENVIRONMENTAL SCIENCES↗

Correlative molecular-to-mesoscale evolution in conjugated polymers for intrinsically stretchable organic photovoltaics

Conjugated polymer thin films offer a unique combination of tunable optoelectronic properties and mechanical flexibility, making them as promising materials for intrinsically stretchable optoelectronic devices. However, achieving both mechanical robustness and high device performance remains a key challenge. Addressing this requires a fundamental understanding of how molecular and mesoscale structures evolve under mechanical strain. Here, we employ a comprehensive suite of X-ray spectroscopy and scattering techniques to investigate the multiscale structural evolution of conjugated polymer thin films during uniaxial deformation. We uncover a two-stage morphological response: an initial stage characterized by polymer chain alignment and rapid crystallite disruption, followed by continued chain orientation accompanied by intrachain torsion at higher strains. These correlative structural adaptations govern key material properties, including stress dissipation, optical absorption, and photovoltaic performance. Our findings establish a mechanistic framework for understanding deformation in semiconducting polymers and provide design principles for developing mechanically robust, high-performance stretchable electronics.

36 MATERIALS SCIENCE↗

Amplified Mesoscale and Submesoscale Variability and Increased Concentration of Precipitation under Global Warming over Western North America

Abstract Cold-season precipitation statistics in simulations from the storm-resolving WRF Model at 6-km and 1-h resolution over western North America are analyzed. Pseudo–global warming future simulations for the 2041–80 period, constrained by GCMs under the RCP8.5 scenario, are compared to the 1981–2020 historical simulation. The analysis focuses on the dynamical properties of precipitation time series at subdaily scales and on the morphology of storms. The statistical distribution of precipitation intensities in each pixel of the simulation domain is characterized through nonparametric statistical indicators: frequency of wet hours, mean wet-hour precipitation intensity, and Gini coefficient as a measure of the temporal concentration of the precipitation volume. Additionally, the temporal and spatial Fourier power spectra of precipitation time series and precipitation fields are analyzed. The half-power period (HPP) and half-power wavelength (HPW) are defined as spectral measures of the characteristic scales of precipitation’s temporal and spatial patterns. The results show statistically significant increases in the mean wet-hour precipitation intensity and in the Gini coefficient in 99% of the pixels, indicating that the seasonal precipitation volume becomes more concentrated within a smaller number of hours with higher precipitation intensity. The statistics of change in the frequency of wet hours are more contrasted across the simulation domain. The changes are also reflected in the power spectra, which show the spatial and temporal variability increasing proportionally more with finer spatial and temporal scales and the HPW and HPP decreasing. These projected changes are expected to have consequences, not only in terms of hydrologic impacts but also in terms of the predictability of precipitation patterns. Significance Statement The precipitation characteristics of winter storms over the western United States and southwestern Canada are analyzed in future climate simulations for the 2041–80 period. As compared to present-day climate, the most intense parts of the storms are projected to produce a higher rainfall volume, with increased concentration over smaller areas and shorter time intervals. The propensity of rainfall intensity to vary rapidly over time will be enhanced in the future according to the simulations. These model predictions imply an increased risk of rapid flooding in small basins. They also suggest that predicting several hours ahead the time and location at which a storm will produce maximum rainfall may become more challenging in the future.

Climate change↗

Mesoscale Organization in Cumulus-Coupled Stratocumulus

Marine cloud systems cover a substantial portion of the world’s oceans. Most of these clouds form relatively close to the ocean surface, typically within one to two kilometers, a region referred to by meteorologists as the marine boundary layer. They are composed predominantly of liquid water, although ice particles can occur in mid- and high-latitude marine clouds during winter. In satellite imagery, these clouds appear bright against the darker ocean surface below, reflecting a large fraction of incoming sunlight back into space that would otherwise warm the ocean. Because marine boundary layer clouds cover such an extensive area of the ocean, they exert a significant influence on Earth’s overall transfer of solar energy absorbed by the surface and thermal energy emitted to space, a balance known as the planetary radiation budget. Marine boundary layer clouds are typically thin, and their formation and dissipation depend on a delicate balance between processes acting at the ocean surface below and the warm, dry air above. They are notoriously difficult to simulate accurately in weather forecast models, which often produce too few marine low clouds in midlatitudes and clouds in tropical regions that are excessively bright, meaning they reflect too much solar radiation. The marine boundary layer is frequently characterized by widespread overcast cloud cover that often transitions from a continuous, single-layer deck to more broken cloud fields toward the tropics. These transitions typically proceed through an intermediate stage in which shallow, broken clouds form beneath the overlying stratiform cloud deck. Once broken clouds develop below the overcast, they frequently self-organize into cloud clusters known as marine boundary layer convective complexes (MBLCCs), although the mechanisms governing the formation and organization of MBLCCs remain poorly understood. Accurately representing these transitions in long-range weather forecast models is essential because they influence the properties of air masses advected over the continental United States and Europe, and they become increasingly important for forecasts on seasonal and longer timescales. We employed two complementary approaches to investigate the processes controlling MBLCCs and their impact on marine cloud cover. Long-term observations from the U.S. Department of Energy’s Eastern North Atlantic (ENA) Observatory provided a unique dataset that allowed us to characterize fundamental properties of MBLCCs, including their typical size and frequency of occurrence. These observations were combined with high-resolution numerical simulations performed on supercomputers to examine the evolution of MBLCCs during cold-air outbreaks over the ENA region.

54 ENVIRONMENTAL SCIENCES↗

Reconciliation of asynchronous satellite-based $\mathrm{NO_2}$ and $\mathrm{XCO_2}$ enhancements with mesoscale modeling over two urban landscapes

Fossil fuel carbon dioxide (CO 2ff ), the main driver of global warming and climate change, is often co-emitted with nitrogen oxides (NO x ) and precursors to ground-level ozone from anthropogenic sources like power plants or vehicles. In urban and suburban areas, satellite-based NO 2 can be used as a proxy to track the emissions of CO 2ff . Because of NO 2 ’s shorter lifetime, urban NO 2 plumes are more distinguishable from backgrounds and more sensitive to variations in emissions. However, the combination of these two gases is limited by the asynchrony among NO 2 and CO 2 monitoring satellites. We used CO 2ff simulated by the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) model to reconcile the tropospheric NO 2 vertical column density (VCD) from the Tropospheric Monitoring Instrument (TROPOMI) and column-averaged dry-air mole fractions of carbon dioxide enhancements (ΔXCO 2 ) from Orbiting Carbon Observatory 3 (OCO-3) Snapshot Area Maps (SAMs) over a multicity area, Washington D.C.-Baltimore (DC-Balt), and a basin city, Mexico City. NO 2 /CO 2ff ratios over DC-Balt are smaller than Mexico City, indicative of stricter emission restrictions, a more combustion-efficient vehicle fleet, and higher combustion efficiency due to lower altitude in DC-Balt. For single-track cases, the spatial correlations between NO 2 and ΔXCO 2 over Mexico City are stronger than DC-Balt because the NO 2 and CO 2 are mostly trapped in the valley of Mexico City, while DC-Balt is severely affected by distant sources (i.e., US East Coast cities). Using multi-track averaging, spatial correlation coefficients increase with the number of days used for averaging. The correlations reached a maximum when averaging >12 continuous images for DC-Balt and > 10 continuous images for Mexico City. This finding indicates that multi-track averaging using modeled CO 2ff as a proxy is helpful to filter the noise in single-track images, to cancel the interference from distant sources, and to magnify correlations between NO 2 and CO 2ff . Mexico City showed stronger spatial correlations but weaker temporal correlations than DC-Balt due to biomass burning hot spots and large transport errors caused by the trapping effects of the surrounding mountains. Tracking the 20-day moving average of CO 2ff emissions using TROPOMI NO 2 seems technically feasible, considering the relationship between correlation coefficients and the number of available satellite images.

54 ENVIRONMENTAL SCIENCES↗

Develop a new integrated macro→micro←nano (MMN) multiscale modeling framework to optimize high strength aluminum alloys and processes for vehicle light-weighting​

Bending tests provide a means to study plane strain fracture performance of 6000 series aluminum alloys. Metrics from bending tests have been correlated with self-pierce riveting (SPR) performance of a high strength AA6111 automotive aluminum alloy in previous works. Using the ORNL HPC resources, this project developed an innovative macro→micro←nano (MMN) multiscale microstructure-based finite element (FE) code to further understanding of the relationship between microstructure and fracture properties of high-strength 6000-series alloys. This work started with microstructural characterization in both mesoscale and nanoscale and bending performance characterization of AA6111 HS2-T6 alloy at Ford, the MMN framework was applied to this alloy to simulate 3-point VDA bending. From the results of the macro-modeling of 3-point VDA bending, the critical region of fracture was identified, and the region geometry was used to construct the micro-model. The fracture criterion of micron-scale precipitates and aluminum matrix which contains submicron and nano particles (AL-SMP-NP) in the micro-model was calibrated and validated by comparing simulated and measured bending results. With the AL-SMP-NP fracture strain obtained, the fracture strain of Al-matrix containing nano particles (ALNP) will similarly be determined by a submicron scale-model using an edge-constrained FE modeling approach developed by Hu et al. With the ALNP fracture strain obtained, the fracture strain of Al-matrix containing no particles will similarly be determined by a nano scale-model using an edge-constrained FE modeling approach. After the MMN framework is built and fracture criterion calibrated, nano-model FE simulations with virtual microstructures was performed to obtain a reduced order model (ROM) of the fracture criterion of the ALNP as a function of volume fraction, size, and distribution of the nanoparticle. This nano→submicron→micro modeling part allows exploration of the influence of different material nanostructures from different process conditions on the bending properties within the multiscale bending simulation framework and the ROM of material bendability as a function of nanoparticle size and shape was established. This obtained reduced order model (ROM) could help guide the design and selection of materials to improve existing SPR process models that could replace trial-and-error rivet/die selection and help to design new rivet/die combinations capable of robustly joining new higher strength 6000 5 series alloys in automotive body structures. This would enable lightweighting of Ford vehicles leading to greater fuel efficiency and reduce manufacturing time and energy.

36 MATERIALS SCIENCE↗

Co-Occurring Atmospheric Features and Their Contributions to Precipitation Extremes

Object-based identification algorithms for atmospheric features are commonly utilized to attribute global precipitation. This study employs a systematic approach to examine feature co-occurrences and their relationships to mean and extreme precipitation. Four features are identified using existing data sets for atmospheric rivers (ARs), mesoscale convective systems (MCSs), low-pressure systems (LPSs), and fronts (FTs). Often, a single atmospheric phenomenon satisfies the criteria set by multiple feature identification algorithms, yielding an association between precipitation and multiple features. Over the extra-tropics, the number of features attributed to a single event typically increases with precipitation intensity. Over two-thirds of the precipitation is from co-occurring features, with a considerable fraction related to AR-FT co-occurrences. Over the tropics, about one-quarter of precipitation is associated with co-occurring features, with LPS-MCS co-occurrences contributing substantially in monsoon regions. MCSs are the leading single-feature contributors over tropical land and oceans. In the extra-tropics, FTs, ARs, and their co-occurrences account for over half of the total precipitation over oceans. AR-FT-MCS and FT-MCS co-occurrences contribute to extremes (precipitation exceeding the 95th percentile) over both oceans (over 30%) and land (over 20%). Any combination of features involving MCSs shows a larger contribution to high percentiles of precipitation intensity. A case analysis indicates that AR-FT-MCS co-occurrences exhibit convective instability and deep vertical motion, suggesting that the feature trackers and reanalysis are capturing physics relevant to both convective and frontal systems. The results here emphasize the need for simultaneous identifications of multiple features when attributing precipitation to atmospheric phenomena.

54 ENVIRONMENTAL SCIENCES↗

Monsoonal MCS Initiation, Rainfall, and Diurnal Gravity Waves over the Bay of Bengal: Observation and a Linear Model

Abstract Previous observational studies have indicated that mesoscale convective systems (MCSs) contribute the majority of precipitation over the Bay of Bengal (BoB) during the summer monsoon season, yet their initiation and propagation remain incompletely understood. To fill this knowledge gap, we conducted a comprehensive study using a combination of 20-yr satellite observations, MCS tracking, reanalysis data, and a theoretical linear model. Satellite observations reveal clear diurnal propagation signals of MCS initiation frequency and rainfall from the west coast of the BoB toward the central BoB, with the MCS rainfall propagating slightly slower than the MCS initiation frequency. Global reanalysis data indicate a strong association between the offshore-propagating MCS initiation frequency/rainfall and diurnal low-level wind perturbations, implying the potential role of gravity waves. To verify the hypothesis, we developed a 2D linear model that can be driven by realistic meteorological fields from reanalysis. The linear model realistically reproduces the characteristics of offshore-propagating diurnal wind perturbations. The wind perturbations, as well as the offshore propagation signals of MCS initiation frequency and rainfall, are associated with diurnal gravity waves emitted from the coastal regions, which in turn are caused by the diurnal land–sea thermal contrast. The ambient wind speed and vertical wind shear play crucial roles in modulating the timing, propagation, and amplitude of diurnal gravity waves. Using the linear model and satellite observations, we further show that the stronger monsoonal flows lead to faster offshore propagation of diurnal gravity waves, which subsequently control the offshore propagation signals of MCS initiation and rainfall. Significance Statement Rainfall over the Bay of Bengal (BoB) is primarily contributed by large and organized rainfall systems in the summer monsoon season. During this season, these systems are commonly observed over the east coast of India around midnight, the western BoB in the morning, and the central BoB in the afternoon. This eastward rainfall propagation is confirmed by observations, reanalysis data, and a theoretical model to have a strong association with atmospheric diurnal gravity waves. These waves are caused by land–sea thermal contrast and can trigger rainfall systems over the offshore regions. We also found that the diurnal gravity waves, as well as their triggered rainfall systems, can be greatly modulated by the large-scale monsoonal flows. All the above findings improve our understanding of diurnal rainfall cycle over the BoB and may contribute to the future improvement of rainfall forecast over the region.

Meteorology & Atmospheric Sciences↗

Using ARM Observations to Evaluate Process-Interactions in MCS Simulations Across Scales (Final Progress Report)

This project, funded by DOE Atmospheric System Research (DE-SC0020050), focused on improving the representation of mesoscale convective systems (MCSs) in numerical weather and climate models by leveraging high-resolution observations from the DOE Atmospheric Radiation Measurement (ARM) program. The research aimed to evaluate model sensitivities to grid spacing, microphysics, and planetary boundary layer (PBL) schemes, with a particular emphasis on improving convection parameterization for high-resolution modeling. Findings from this work highlight several key advancements. Model validation against ARM radar wind profiler data from the Southern Great Plains (SGP) and Manaus (MAO) sites revealed systematic biases in simulated convective mass flux profiles, leading to the development of an observationally constrained evaluation framework for diagnosing and improving model performance. Sensitivity analyses demonstrated that the representation of Amazonian MCSs was highly dependent on PBL scheme selection, while mid-latitude MCSs were more strongly influenced by microphysics parameterizations. A series of high-resolution WRF simulations, ranging from 4 km to 125 m grid spacing, provided insight into the behavior of convective drafts across scales. While updraft properties converged at sub-kilometer resolutions, biases in downdraft intensity persisted even at the finest resolution tested, emphasizing the need for further refinements in model physics. Additionally, comparisons of MCS vertical structures between mid-latitude and tropical environments revealed stronger updrafts and larger mass flux in mid-latitude MCSs, providing critical insights for improving climate model representations of storm-scale dynamics. The project’s findings have already contributed to advancing numerical modeling capabilities, particularly in WRF, MPAS, ICON, and DOE’s SCREAM model, by refining how convective processes are represented in high-resolution climate simulations. Results were disseminated through peer-reviewed publications, conference presentations, and ARM/ASR Research Highlights, engaging the broader scientific community. The project also provided valuable training opportunities for two postdoctoral researchers, who played central roles in model development, analysis, and dissemination of results. Their work contributed to several publications and conference presentations, helping prepare them for careers in atmospheric modeling. By improving the simulation of MCSs, this research directly supports the development of next-generation climate models capable of more accurately representing extreme precipitation and convective processes. The insights gained will inform future improvements in convective parameterization and guide the design of high-resolution weather and climate simulations, ultimately enhancing the reliability of climate projections and weather forecasts.

54 ENVIRONMENTAL SCIENCES↗

Collaborative Research: Advancing Arctic Climate Projection Capability at Seasonal to Decadal Scales (Final Technical Report)

The Regional Arctic System Model (RASM) at process resolving configurations has been used to (i) advance understanding of physical processes and feedbacks involved in Arctic amplification and (ii) understand and potentially reduce uncertainty in prediction of arctic climate change at seasonal to decadal scales. RASM consists the atmosphere (Weather and Research Forecasting model, WRF), ocean (Parallel Ocean Program, POP), sea ice (CICE), land hydrology (Variable Infiltration Capacity model, VIC), river routing scheme (RVIC), marine biogeochemistry components and the coupling framework (CPL7). Its domain is pan-Arctic, with the atmosphere and land components configured on a 50-km or 25-km grid and four configurations of the ocean and sea ice components: 1/12°(~9.3km) or 1/48°(~2.4km) and 45 or 60 vertical layers. These RASM configurations have been motivated by the emerging exascale capability for high performance computing to improve model fidelity. The dynamical downscaling of reanalysis allows comparison of RASM results with observations in place and time to: (i) advance system level understanding of physical processes and coupling involved in an event, (ii) optimize model parameter space, (iii) diagnose and reduce model biases and (iv) produce realistic and consistent across all the components initial conditions for predictions and predictability studies, which are all unique capabilities not available in global Earth System Models (ESMs). An evaluation of RASM 1.0 (Cassano et al. 2017) revealed that it had a cold bias over the oceans and a warm bias over land areas due largely to cloud and radiation biases in the model, with too little cloud cover simulated over land and too much cloud cover simulated over sub-polar oceans. This study has motivated an upgrade to WRF version 3.7.1 in RASM and allowed for the inclusion of the radiative impact of convective clouds. A variety of atmospheric physics parameterizations were evaluated against observations (e.g. data from the Arctic Clouds in Summer Experiment (ACSE); Sedlar et al. 2020) to identify an optimal suite of WRF physics options in RASM. The RASM with the optimized WRF physics were used to study the impact of strong mesoscale winds over the ocean around the southern tip of Greenland (DuVivier and Cassano 2016) and their impact on oceanic convection (DuVivier et al. 2017a). Data from the PolarWinds field campaign were used to evaluate WRF boundary layer physics and resolution impacts on the simulation of a Greenland barrier wind event (DuVivier et al. 2017b). The RVIC streamflow routing model has been implemented in RASM to realistically represent high-resolution streamflow processes (Hamman et al. 2017) and to couple the land buoyancy fluxes to the ocean. The RASM-RVIC high-resolution data set of all coastal freshwater fluxes for the Arctic drainage basin and surrounding areas for 1979-2014 was published as a separate product (https://doi.org/10.5281/zenodo.293037). The fidelity of atmospheric momentum transfer to and the response of polar marine Ekman layer in RASM and Community Earth System Model (CESM) was investigated by Roberts et al. (2015). The increased frequency of oceanic flux exchange in CESM, following the RASM guidance, caused a considerable increase in the median inertial ice speed across the Southern Ocean and parts of the Arctic. A comprehensive evaluation of the RASM1.0 atmosphere-ocean-sea ice-land interface was completed by Brunke et al. (2018). RASM was also demonstrated for its capability to simulate extreme events in agreement with observations in space and time (Lee et al. submitted). In particular, the development of three open water events, known as polynyas, have been simulated north of Greenland in February of 2011, 2017 and 2018, in agreement with satellite observations for the past four decades. The optimized RASM sea ice results have been favorably evaluated against satellite observations and a subset of eleven CMIP6 models (Watts et al. submitted). In a complementary project, Jin et al. (2018) have shown that RASM with higher-resolution and new sea-ice processes contributed to lower model errors in sea-ice conditions, concentrations of nutrients and ice algae, in comparison to results from the coarse-resolution (1°) CESM. In two other complementary studies, RASM results were used (i) to explain areas of concentrated use by bowhead whales, the seasonal progression in the use, and the physical environment within those areas (Citta et al. 2015) and (ii) for a synthesis of fall bowhead whales distribution and migration in the Bering-Chukchi-Beaufort (BCB) Sea to investigate whale movements and feeding to the local ocean hydrography and currents (Citta et al. 2018). However, the multi-decadal output from the CESM Large Ensemble yielded unrealistic forcing. Instead, the shorter NCEP CFSv2 9-month forecasts were successfully tested and afforded an increased ensemble size (~30) to demonstrate gains of dynamical downscaling at sub-seasonal to intra-annual time scales. The improved model physics and coupling among RASM model components have yielded more realistic representation of the sea ice cover and consistent across all model components initial conditions. Consequently, RASM demonstrates significant gains compared to simulation of sea ice in the NCEP reanalysis. In addition, RASM 6-month ensemble forecasts yield very realistic sea ice distribution, which demonstrates both significant gains of dynamical downscaling and the continued impact of the initial conditions on forecasts out to 6 months (https://nps.edu/web/rasm/predictions). A follow up study using RASM for dynamical downscaling of the more realistic CESM initialized Decadal Prediction Large Ensemble output is currently ongoing as part of the DOE RGMA HiLAT-RASM project.

54 ENVIRONMENTAL SCIENCES↗