Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical climate model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Savings in Action: Lessons from Observed and Modeled Residential Solar Plus Storage Systems

The electric grid is rapidly evolving as small-scale, demand-side resources play increasingly important roles in grid operations and decarbonization. Maximizing the potential of demand-side resources involves incentivizing electricity customers to use those resources in ways that benefit the broader electrical grid. These incentives depend largely on the electricity cost savings that customers can realize from demand-side resource adoption. Determining these potential cost savings is a complex task. Cost savings depend on numerous factors, including the characteristics of different technologies, the algorithms that control these devices, system performance, customer behavior, electricity rate structures, and climatic factors. Another challenge is that estimated cost savings are frequently based on modeled rather than observed system performance, particularly in the academic literature. In this study, we begin to fill the gap in empirical research of demand-side resources using data from a new construction residential community equipped with rooftop solar and storage (S+S) in Arizona. We use these data to analyze the factors that determine customer electricity cost savings and emissions impacts of S+S in the real world. We then compare these data to modeled system performance to understand how models deviate from real-world outcomes. Based on these findings, we explore ways to improve such models and, conversely, use modeled results to suggest improvements to actual S+S deployment. The results of these analyses can be summarized in four key findings: 1) rate structures play a central role in the grid and customer value of demand-side resources; 2) certain customers can benefit more from demand-side resource adoption than others; 3) modeled battery dispatch and sizing reveals opportunities for additional cost savings; and 4) optimal dispatches can reduce grid emissions while maximizing bill savings.

14 SOLAR ENERGY↗
DOE OSTI2023

A Machine Learning Correction Model of the Winter Clear-Sky Temperature Bias over the Arctic Sea Ice in Atmospheric Reanalyses

Atmospheric reanalyses are widely used to estimate the past atmospheric near-surface state over sea ice. They provide boundary conditions for sea ice and ocean numerical simulations and relevant information for studying polar variability and anthropogenic climate change. Previous research revealed the existence of large near-surface temperature biases (mostly warm) over the Arctic sea ice in the current generation of atmospheric reanalyses, which is linked to a poor representation of the snow over the sea ice and the stably stratified boundary layer in the forecast models used to produce the reanalyses. These errors can compromise the employment of reanalysis products in support of polar research. Here, we train a fully connected neural network that learns from remote sensing infrared temperature observations to correct the existing generation of uncoupled atmospheric reanalyses (ERA5, JRA-55) based on a set of sea ice and atmospheric predictors, which are themselves reanalysis products. The advantages of the proposed correction scheme over previous calibration attempts are the consideration of the synoptic weather and cloud state, compatibility of the predictors with the mechanism responsible for the bias, and a self-emerging seasonality and multidecadal trend consistent with the declining sea ice state in the Arctic. The correction leads on average to a 27% temperature bias reduction for ERA5 and 7% for JRA-55 if compared to independent in situ observations from the MOSAiC campaign (respectively, 32% and 10% under clear-sky conditions). These improvements can be beneficial for forced sea ice and ocean simulations, which rely on reanalyses surface fields as boundary conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

On CCN Effects upon Convective Cold Pool Timing and Features

Abstract Cold pools produced by deep convection can initiate new convection, and their representation in larger-scale weather and climate models could improve prediction of the extent and timing of upscale growth. Cold pools originate from latent cooling from precipitation changing phase, but little attention has been paid to microphysical influences on cold pool characteristics, particularly CCN effects. Datasets obtained from the CACTI and RELAMPAGO field campaigns, along with idealized numerical modeling, are utilized to investigate the hypothesis that convective storms forming in higher-CCN environments generate their first surface rainfall later, delaying cold pool initiation. Aircraft observations of CCN and shallow convection on 9 days do suggest a CCN effect. Those ingesting more CCN contained fewer drizzle drops, although a decreased cloud depth with increasing CCN was also likely a limiting factor. In three of those cases that later developed into deep convection, the timing of cold pool onset was not ubiquitously delayed in environments with more CCN. Idealized numerical simulations suggest that an ordinary thunderstorm can experience small delays in cold pool onset with increasing CCN due to changes in graupel production, but CCN effects on the cold pool from a supercell thunderstorm can be easily overpowered by its unique dynamics. A strong inverse relationship between cold pool strength, expansion rate, and depth with increasing CCN is suggested by the results of the ordinary thunderstorm simulation. Further consideration of CCN appears warranted for future cold pool parameterization development, but other environmental factors affecting storm morphology and precipitation cannot be ignored.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Orographic enhancement of rainfall over the Congo Basin

Abstract The Congo rainforest located in central equatorial Africa is an important, yet understudied part of the globe surrounded by complex orographic features. A primitive understanding of precipitation processes such as mesoscale convective dynamics magnifies uncertainties in the future climate projections of the hydrological cycle over the Congo. Furthermore, the effects of orography, which is an important forcing for convection and precipitation, are poorly resolved by climate models, and ill‐conceptualized over the Congo. To address this knowledge gap, perturbed orographic forcing experiments are conducted using the Weather Research and Forecasting (WRF) mesoscale numerical model in a high‐resolution convection‐permitting model setup. The model simulated selected dates in November 2014. The thunderstorms and rainfall simulated in the control run for the case study analyzed in this article compared reasonably well to satellite‐derived brightness temperature and rainfall data. The results from this case study show that the dynamical impact of increasing the height of the East African Highlands is the blocking of the lower‐tropospheric tropical easterlies. This weakening of the lower‐tropospheric zonal winds increases the windshear over the Congo Basin resulting in slower propagating, more intense mesoscale convective systems with enhanced rainfall.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Clouds and Convective Self-Aggregation in a Multimodel Ensemble of Radiative-Convective Equilibrium Simulations

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climate sensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique among intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than unaggregated simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Micro-Pulse DIAL for Atmospheric Water Vapor Profiling (Phase I Final Report)

Water vapor plays a key role in Earth’s climate and weather patterns. Accurate, high resolution measurements of the water vapor vertical profiles in a region have been found to significantly improve the numerical weather prediction models’ short-term forecasts, precipitation forecasts, as well as severe weather forecasts. However, current measurement techniques provide very coarse resolution in space, time, and/or altitude. A network of continuously sampling water vapor profilers would greatly expand the scientific knowledge of the atmosphere and serve as a valuable input into numerical weather models. This would be of interest not only to the Department of Energy’s Environmental Sciences Division but to other institutions such as the National Weather Service, the National Center for Atmospheric Research, the Federal Emergency Management Agency, and the military.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Long‐Term Impacts of Global Solid Biofuel Emissions on Ambient Air Quality and Human Health for 2000–2019

Globally, solid biofuels (SB) have been widely used for household cooking and energy production for decades due to electricity shortages and socio-economic barriers to adopting renewable energy alternatives. This has detrimental effects on air quality, human health, and climate through trace gas and aerosol emissions. Despite numerous studies, the long-term consequences of SB emissions remain poorly understood. Here, we use the Community Earth System Model and the Community Emissions Data System emission inventory to investigate the SB emission impacts on air quality and human health for 2000–2019. Global SB emission increased the ambient PM 2.5 (particulate matter with aerodynamic diameters ≤2.5 μm) and ozone (O 3 ) concentrations up to 23.61 μg/m 3 and 13.69 ppbv, with significant effects found in India, China, and the Rest of Asia (ROA). Our study estimates total annual premature deaths (APDs) associated with global SB-attributable PM 2.5 and O 3 exposure as 1.11 million [95% confidence interval (95% CI): 1.00–1.22 million] in 2000 up to 1.43 million (95% CI: 1.30–1.56 million) in 2019. China's SB emissions and associated APDs have reduced substantially, whereas India and ROA had a major leap in both estimates in 2019 compared to 2000. China's progress in cutting residential SB emissions accounts for its improvements. Our study urges the reduction of SB usage and emissions to potentially improve overall air quality and human health conditions, especially in highly populated, low- and middle-income countries, where the poor air quality and associated health burden attributable to SB emissions are estimated to be higher.

O 3↗
DOE OSTI2024

Upland Yedoma taliks are an unpredicted source of atmospheric methane

Landscape drying associated with permafrost thaw is expected to enhance microbial methane oxidation in arctic soils. Here we show that ice-rich, Yedoma permafrost deposits, comprising a disproportionately large fraction of pan-arctic soil carbon, present an alternate trajectory. Field and laboratory observations indicate that talik (perennially thawed soils in permafrost) development in unsaturated Yedoma uplands leads to unexpectedly large methane emissions (35–78 mg m -2 d -1 summer, 150–180 mg m -2 d -1 winter). Upland Yedoma talik emissions were nearly three times higher annually than northern-wetland emissions on an areal basis. Approximately 70% emissions occurred in winter, when surface-soil freezing abated methanotrophy, enhancing methane escape from the talik. Remote sensing and numerical modeling indicate the potential for widespread upland talik formation across the pan-arctic Yedoma domain during the 21 st and 22 nd centuries. Contrary to current climate model predictions, these findings imply a positive and much larger permafrost-methane-climate feedback for upland Yedoma.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Deep Learning Modeling Framework to Capture Mixing Patterns in Reactive-Transport Systems

Prediction and control of chemical mixing are vital for many scientific areas such as subsurface reactive transport, climate modeling, combustion, epidemiology, and pharmacology. Due to the complex nature of mixing in heterogeneous and anisotropic media, the mathematical models related to this phenomenon are not analytically tractable. Numerical simulations often provide a viable route to predict chemical mixing accurately. However, contemporary modeling approaches for mixing cannot utilize available spatial-temporal data to improve the accuracy of the future prediction and can be compute-intensive, especially when the spatial domain is large and for long-term temporal predictions. To address this knowledge gap, in this work we will present in this paper a deep learning (DL) modeling framework applied to predict the progress of chemical mixing under fast bimolecular reactions. This framework uses convolutional neural networks (CNN) for capturing spatial patterns and long short-term memory (LSTM) networks for forecasting temporal variations in mixing. By careful design of the framework—placement of non-negative constraint on the weights of the CNN and the selection of activation function, the framework ensures non-negativity of the chemical species at all spatial points and for all times. Our DL-based framework is fast, accurate, and requires minimal data for training. The time needed to obtain a forecast using the model is a fraction (≈ O(-6)) of the time needed to obtain the result using a high-fidelity simulation. To achieve an error of 10% (measured using the infinity norm) for capturing local-scale mixing features such as interfacial mixing, only 24% to 32% of the sequence data for model training is required. To achieve the same level of accuracy for capturing global-scale mixing features, the sequence data required for model training is 64% to 70% of the total spatial-temporal data. Hence, the proposed approach—a fast and accurate way to forecast long-time spatial-temporal mixing patterns in heterogeneous and anisotropic media—will be a valuable tool for modeling reactive-transport in a wide range of applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Comparative analysis of thermal management systems in electric vehicles at extreme weather conditions: Case study on Nissan Leaf 2019 Plus, Chevrolet Bolt 2020 and Tesla Model 3 2020

With the surge in electric vehicle (EV) adoption and the need for extended driving ranges, optimizing energy efficiency, particularly through thermal management, is critical, especially in extreme weather. Managing the substantial energy needed for cabin climate control and battery temperature regulation can increase energy demands by over 50 %, severely limiting range. This study conducts a comparative analysis of thermal management systems (TMS) in three popular EV vehicles, 2020 Chevrolet Bolt, 2019 Nissan Leaf Plus, and 2020 Tesla Model 3, evaluating their distinct TMS configurations and performance under varied weather conditions. Using both numerical simulations and experimental data collected on a controlled test bench at Argonne National Laboratory, we assess how TMS architecture and operational modes influence energy consumption and range. A comprehensive TMS model was developed, integrating cabin and battery thermal sub-models in the Autonomie software platform, to simulate temperature fluctuations and range impacts. Cabin climate was modeled using a mono-zonal approach, while battery cell temperature distribution was estimated through a 2D nodal structure. Each vehicle's distinct TMS setup was evaluated: the Chevrolet Bolt and Tesla Model 3 use a dual evaporator vapor compression cycle with a PTC heater for the cabin and a coolant loop for battery thermal management; the Nissan Leaf Plus employs a heat pump with a PTC heater for the cabin and air-cooling for the battery. Tests conducted at ambient temperatures of 35°C, 22°C, -7°C, and -18°C reveal significant differences in energy use and range reduction across both configurations and conditions. At 35°C, the Tesla Model 3, Chevrolet Bolt, and Nissan Leaf Plus have a range reduction of 8%, 9%, and 13%, respectively, due to air conditioning. In winter, heating technology is paramount; at -7°C, the Nissan Leaf's heat pump configuration achieves a lower range reduction (19.3%) compared to the Tesla and Chevrolet Bolt PTC heaters, which reduce range by 28.3% and 31%, respectively. Further, this study provides valuable insights for automotive engineers, EV technology researchers, and thermal management system designers aiming to enhance electric vehicle performance by understanding how different weather conditions and TMS architectures impact energy consumption and driving range.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Sensitivity of Fine‐Resolution Urban Heat Island Simulations to Soil Moisture Parameterization

ABSTRACT Urban areas experience the impact of natural disasters, such as heatwaves and flash floods, disparately in different neighbourhoods across a city. The demand for precise urban hydrometeorological and hydroclimatological modelling to examine this disparity, and the interacting challenges posed by climate change and urbanisation, has thus surged. The Weather Research and Forecasting (WRF) model has served such operational and research purposes for decades. Recent advancements in WRF, including enhanced numerical schemes and sophisticated urban atmospheric‐hydrological parameterizations, have empowered the simulation of urban geophysical processes at high resolution (~1 km), but even this resolution misses significant urban microclimate variability. This study applies the large‐eddy simulations (LES) mode within WRF, coupled with single‐layer urban canopy models (SLUCM), to enable even finer‐scale modelling (150 m) of the Urban Heat Island (UHI) effect in the Baltimore metropolitan area. We run nine scenarios to evaluate various methods of initializing soil moisture and various spinup lead times, and to assess the impact of WRF's Mosaic approach in depicting subgrid‐scale processes. We evaluate the scenarios by comparing the WRF simulated land surface temperature (LST) against Landsat LST and the WRF simulated hourly 2‐m air temperatures (AT) with observations from eight weather stations across the domain. Results underscore the paramount influence of the lead spinup time on the spatiotemporal distribution of simulated soil moisture, consequently shaping WRF's efficacy in predicting the UHI. Furthermore, interpolating soil moisture‐related parameters from the parent for child domain initialization yields a notable reduction in mean and root‐mean‐squared errors. This improvement was particularly evident in simulations with the longest spinup time, affirming the importance of carefully designing the initialization of soil moisture for improved urban temperature predictions.

Talebpour, Mahdad↗
DOE OSTI2024

Are turbulence effects on droplet collision–coalescence a key to understanding observed rain formation in clouds?

Rain formation is a critical factor governing the lifecycle and radiative forcing of clouds and therefore it is a key element of weather and climate. Cloud microphysics–turbulence interactions occur across a wide range of scales and are challenging to represent in atmospheric models with limited resolution. Based on past experiments and idealized numerical simulations, it has been postulated that cloud turbulence accelerates rain formation by enhancing drop collision–coalescence. We provide substantial evidence for significant impacts of turbulence on the evolution of cloud droplet size distributions and rain formation by comparing high-resolution observations of cumulus congestus clouds with state-of-the-art large-eddy simulations coupled with a Lagrangian particle-based microphysics scheme. Turbulent coalescence must be included in the model to accurately represent the observed drop size distributions, especially for drizzle drop sizes at lower heights in the cloud. Turbulence causes earlier rain formation and greater rain accumulation compared to simulations with gravitational coalescence only. The observed rain size distribution tail just above cloud base follows a power law scaling that deviates from theoretical scalings considering either a purely gravitation collision kernel or a turbulent kernel neglecting droplet inertial effects, providing additional evidence for turbulent coalescence in clouds. In contrast, large aerosols acting as cloud condensation nuclei (“giant CCN”) do not significantly impact rain formation owing to their long timescale to reach equilibrium wet size relative to the lifetime of rising cumulus thermals. Overall, turbulent drop coalescence exerts a dominant influence on rain initiation in warm cumulus clouds, with limited impacts of giant CCN.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

High-Resolution Regional Wave Hindcast for the U.S. Alaska Coast

Regional long-term wave hindcast with state-of-art third-generation models is essential for characterization of regional wave resources. This report summarizes modeling efforts for the simulation of the wave climate along southern Alaska, which used an unstructured, nested-grid modeling approach that incorporates a global-regional nested grid using WAVEWATCH III ® and the high-resolution SWAN (Simulating Waves Nearshore) model. Wave resource and wave bulk parameters were simulated for a 32-year period from 1979 to 2010 and were subsequently validated with wave buoy data in the model domain. Error statistics of model skills were calculated at all buoy stations. Overall, model results, both resource and bulk parameters, match observations well. The model hindcast was able to reproduce the seasonal variation of the sea state with large waves that occur in the winter and early spring months when wind forcing is strong and the calm sea state during the summer when wind forcing is weak. The hindcast also considers the effect of sea ice on wave growth, which affects the southern Bering Sea. The hindcast results, including resource and bulk parameters, were archived at every model grid point for the entire model domain at three-hour intervals. In addition, spectral results were saved at National Data Buoy Center buoy stations and selected virtual buoy locations at hourly intervals. The nested-grid modeling framework employed in this study provides a powerful and efficient modeling approach for accurately simulating wave climate at regional and long-term temporal scales with sufficiently fine resolutions in the nearshore region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A Bottom-up Approach to Improve the Representation of Deep Convective Clouds in Weather and Climate Models (Final Report)

A fundamental yet outstanding question about deep cumulus convection regards how convective updrafts, downdrafts, and surface-based cold pools are inter-related within organized storms, and how these three convective components are modulated by external and internal factors. We considered this question through numerical model simulations as well as data collected during the MC3E field campaign. One of our basic hypotheses was that deep cold pools are forced by wide downdrafts which arise from wide updrafts. This “bottom-up” perspective required consideration of the external controls on updraft width, and then the internal controls on cold-pool characteristics including depth. Model experiments showed that environmental vertical wind shear (an external factor) exerts a particularly strong control on updraft width; analyses of linear and nonlinear dynamic pressure forcing readily explain this control. The experiments also showed that increases in updraft width are accompanied by corresponding increases in downdraft width and cold-pool depth. The updraft–downdraft connection is through precipitation: the wide area of precipitation that will be initiated in a wide updraft will in turn lead to a wide downdraft. The downdraft–cold pool connection is through mass continuity: at the downdraft terminus, the large downdraft mass flux through a wide downdraft must lead to a deep cold pool. Model experiments also showed that the cold-pool depth and other characteristics can be significantly modulated by microphysical processes (internal factors). Specifically, the cold-pool expansion rate, depth, and intensity was found to depend most strongly on the amount of rain evaporation near the ground. However, the sustenance of the cold pool itself depended upon the downdrafts strengthened by sublimation of graupel higher in the storm. Finally, field campaign data revealed compelling inter-relationships between the cold-pool characteristics, and downdraft transience in terms of the microphysics. For example, pulses of hail and graupel (also internal factors) within downdrafts resulted in short-lived increases in the cold-pool depth. However, the cold-pool speed appeared to be relatively independent of the cold-pool depth and intensity, suggesting limitations of density current theory to convective cold pools in nature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Impact of vegetation physiology and phenology on watershed hydrology in a semiarid watershed in the Pacific Northwest in a changing climate

Concurrent with changes in climate conditions, increases in atmospheric composition and transport can affect land surface processes by regulating physiological and phenological responses of plants, and consequently perturb water, carbon, and nitrogen (N) dynamics at multiple scales. Howerver, the relationships bwtween plants physiological and phenological and land surface processes in a changing climate are less well known due to the limitation of representing palnts properties in traditional hydrological models. In this study, we use the Community Land Model version 5 (CLM5) to investigate how projected changes in large scale atmospheric conditions can affect hydrologic and biogeochemical dynamics in the Upper Columbia-Priest Rapids (UCPR) watershed, a typical semi-arid watershed located in the north-western United States dominated by croplands. Six numerical experiments considering climate change, CO2 concentration, and N deposition consistent with the Representative Concentration Pathways (RCPs) 4.5 and 8.5 are performed at 1-km resolution over the period of 2015-2099. Our results show that, compared with historical results (1985-2014), the elevated CO2 directly increases plants N uptake rates that provide additional leaf N to enhance photosynthesis by the end of 21st century (2070-2099). The higher GPP leads to changes in stomatal resistance and increased leaf area index (LAI) that affect evaporatranspiration (ET) and its components. As a result of stomatal responses to a higher CO2 concentration, ET over corn fields declines by 24.6 mm yr-1 under RCP4.5, and by 64.0 mm yr-1 under RCP8.5, respectively by the end of the century from its historical mean of 428 mm yr-1. Meanwhile, ET over soybean fields increases by 27.7 mm yr-1 under RCP4.5 and 9.98 mm yr-1 under RCP8.5 from its historical mean of 507 mm yr-1 as a result of CO2 fertilization as soybeans are able to uptake more N to increase LAI significantly. Furtherfore, the changes in physiological and phenological in corn and soybean can strongly modulate water budget terms over the watershed in response to expected changes in greenhouse gas concentrations and atmospheric composition. Additionally, increases in CO2 concentration and N deposition are beneficial for improving water use efficiency of corn and soybean, therefore lead to less water consumption per unit gram of carbon produced from these crops. Our study elucidates the need to jointly consider effects of surface climate and atmospheric composition and transport in understanding and predicting watershed hydrologic and biogeochemical responses, especially for an agricultural watershed.

Zhu, Bowen↗
DOE OSTI2025

Huge ensembles – Part 1: Design of ensemble weather forecasts using spherical Fourier neural operators

Abstract. Simulating low-likelihood high-impact extreme weather events in a warming world is a significant and challenging task for current ensemble forecasting systems. While these systems presently use up to 100 members, larger ensembles could enrich the sampling of internal variability. They may capture the long tails associated with climate hazards better than traditional ensemble sizes. Due to computational constraints, it is infeasible to generate huge ensembles (comprised of 1000–10 000 members) with traditional, physics-based numerical models. In this two-part paper, we replace traditional numerical simulations with machine learning (ML) to generate hindcasts of huge ensembles. In Part 1, we construct an ensemble weather forecasting system based on spherical Fourier neural operators (SFNOs), and we discuss important design decisions for constructing such an ensemble. The ensemble represents model uncertainty through perturbed-parameter techniques, and it represents initial condition uncertainty through bred vectors, which sample the fastest-growing modes of the forecast. Using the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (IFS) as a baseline, we develop an evaluation pipeline composed of mean, spectral, and extreme diagnostics. With large-scale, distributed SFNOs with 1.1 billion learned parameters, we achieve calibrated probabilistic forecasts. As the trajectories of the individual members diverge, the ML ensemble mean spectra degrade with lead time, consistent with physical expectations. However, the individual ensemble members' spectra stay constant with lead time. Therefore, these members simulate realistic weather states during the rollout, and the ML ensemble passes a crucial spectral test in the literature. The IFS and ML ensembles have similar extreme forecast indices, and we show that the ML extreme weather forecasts are reliable and discriminating. These diagnostics ensure that the ensemble can reliably simulate the time evolution of the atmosphere, including low-likelihood high-impact extremes. In Part 2, we generate a huge ensemble initialized each day in summer 2023, and we characterize the simulations of extremes.

Mahesh, Ankur↗
DOE OSTI2019

Evaluation of Implicit‐Explicit Additive Runge‐Kutta Integrators for the HOMME‐NH Dynamical Core

Abstract The nonhydrostatic High‐Order Method Modeling Environment (HOMME‐NH) atmospheric dynamical core supports acoustic waves that propagate significantly faster than the advective wind speed, thus greatly limiting the time step size that can be used with standard explicit time integration methods. Resolving acoustic waves is unnecessary for accurate climate and weather prediction. This numerical stiffness is addressed herein by considering implicit‐explicit additive Runge‐Kutta (ARK IMEX) methods that can treat the acoustic waves in a stable manner without requiring implicit treatment of nonstiff modes. Various ARK IMEX methods are evaluated for their efficiency in producing accurate solutions, ability to take large time step sizes, and sensitivity to grid cell length ratio. Both the gravity wave test and baroclinic instability test from the 2012 Dynamical Core Model Intercomparison Project are used to recommend 5 of the 27 ARK IMEX methods tested for use in HOMME‐NH.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Local Lagrangian reduced-order modeling for the Rayleigh-Taylor instability by solution manifold decomposition

The Rayleigh-Taylor instability is a classical hydrodynamic instability of great interest in various disciplines of science and engineering, including astrophysics, atmospheric sciences and climate, geophysics, and fusion energy. Analytical methods cannot be applied to explain the long-time behavior of the Rayleigh-Taylor instability, and therefore, numerical simulation of the full problem is required. However, in order to capture the growth of amplitude of perturbations accurately, both the spatial and temporal discretizations need to be extremely fine for traditional numerical methods, and long-time simulation may become prohibitively expensive. In this paper, we propose efficient reduced order model techniques to accelerate the simulation of the Rayleigh-Taylor instability in compressible gas dynamics. Here, we introduce a general framework for decomposing the solution manifold to construct the temporal domain partition and temporally-local reduced order model construction with varying Atwood number. We propose two practical approaches in this framework, namely decomposition by physical time and by penetration distance. Numerical results are presented to examine the performance of the proposed approaches.

97 MATHEMATICS AND COMPUTING↗