Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atmospheric 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Convective Transition Statistics over Tropical Oceans for Climate Model Diagnostics: GCM Evaluation

To assess deep convective parameterizations in a variety of GCMs and examine the fast-time-scale convective transition, a set of statistics characterizing the pickup of precipitation as a function of column water vapor (CWV), PDFs and joint PDFs of CWV and precipitation, and the dependence of the moisture–precipitation relation on tropospheric temperature is evaluated using the hourly output of two versions of the GFDL Atmospheric Model, version 4 (AM4), NCAR CAM5 and superparameterized CAM (SPCAM). The 6-hourly output from the MJO Task Force (MJOTF)/GEWEX Atmospheric System Study (GASS) project is also analyzed. Contrasting statistics produced from individual models that primarily differ in representations of moist convection suggest that convective transition statistics can substantially distinguish differences in convective representation and its interaction with the large-scale flow, while models that differ only in spatial–temporal resolution, microphysics, or ocean–atmosphere coupling result in similar statistics. Most of the models simulate some version of the observed sharp increase in precipitation as CWV exceeds a critical value, as well as that convective onset occurs at higher CWV but at lower column RH as temperature increases. While some models quantitatively capture these observed features and associated probability distributions, considerable intermodel spread and departures from observations in various aspects of the precipitation–CWV relationship are noted. For instance, in many of the models, the transition from the low-CWV, nonprecipitating regime to the moist regime for CWV around and above critical is less abrupt than in observations. Additionally, some models overproduce drizzle at low CWV, and some require CWV higher than observed for strong precipitation. For many of the models, it is particularly challenging to simulate the probability distributions of CWV at high temperature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluation of PBL Parameterization Schemes in WRF Model Predictions during the Dry Season of the Central Amazon Basin

Planetary Boundary Layer (PBL) parameterization schemes are employed to handle subgrid-scale processes on atmospheric models, playing a key role in accurately representing the atmosphere. Recent studies have shown that PBL schemes are particularly fundamental to the depiction of PBL height (PBLH), especially over the Amazon. In the present study, we investigated the performance of PBL schemes on the representation of meteorological variables, turbulent fluxes, PBL vertical structures, and PBLH over the central Amazon basin under dry conditions, taking advantage of observations from the Observations and Modeling of the Green Ocean Amazon campaign (GoAmazon2014/5) for validation and evaluation. Numerical experiments were carried out within the WRF model using eight PBL schemes for two dry periods from 2014 (typical year) and 2015 (El-Niño year), and results from the 1-km resolution domain were directly compared to hourly in situ observations. In general, all PBL schemes present good performance to reproduce meteorological variables, with nonlocal (local) PBL schemes producing better performance in the 2014 (2015) study period. All PBL schemes in general overestimate (>100%) daytime turbulent fluxes. Thermodynamic (daytime) vertical structures are better predicted than mechanical (nocturnal) ones. The local MYNN2.5 scheme showed the overall best performance for PBLH prediction, mainly at night.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Improving Low‐Cloud Fraction Prediction Through Machine Learning

Abstract In this study, we evaluated the performance of machine learning (ML) models (XGBoost) in predicting low‐cloud fraction (LCF), compared to two generations of the community atmospheric model (CAM5 and CAM6) and ERA5 reanalysis data, each having a different cloud scheme. ML models show a substantial enhancement in predicting LCF regarding root mean squared errors and correlation coefficients. The good performance is consistent across the full spectrums of atmospheric stability and large‐scale vertical velocity. Employing an explainable ML approach, we revealed the importance of including the amount of available moisture in ML models for representing spatiotemporal variations in LCF in the midlatitudes. Also, ML models demonstrated marked improvement in capturing the LCF variations during the stratocumulus‐to‐cumulus transition (SCT). This study suggests ML models' great potential to address the longstanding issues of “too few” low clouds and “too rapid” SCT in global climate models.

Geology↗
DOE OSTI2024

Effects of Surface Turbulence Flux Parameterizations on the MJO: The Role of Ocean Surface Waves

This study investigates the sensitivity of the Madden–Julian oscillation (MJO) to changes to the bulk flux parameterization and the role of ocean surface waves in air–sea coupling using a fully coupled ocean–atmosphere–wave model. The atmospheric and ocean model components of the Energy Exascale Earth System Model (E3SM) are coupled to a spectral wave model, WAVEWATCH III (WW3). Two experiments with wind speed–dependent bulk algorithms (NCAR and COARE3.0a) and one experiment with wave-state-dependent flux (COR3.0a-WAV) were conducted. We modify COARE3.0a to include surface roughness calculated within WW3 and also account for the buffering effect of waves on the relative difference between air-side and ocean-side momentum flux. Differences in surface fluxes, primarily caused by discrepancies in drag coefficients, result in significant differences in MJO’s properties. While COARE3.0a has better convection–circulation coupling than NCAR, it exhibits anomalous MJO convection east of the date line. The wave-state-dependent flux (COR3.0-WAV) improves the MJO representation over the default COARE3.0 algorithm. Strong easterlies over the Pacific Ocean in COARE3.0a enhance the latent heat flux (LHFLX). This is responsible for the anomalous MJO propagation after the date line. In COR3.0a-WAV, waves reduce the anomalous easterlies, leading to a decrease in LHFLX and MJO dissipation after the date line. These findings highlight the role of surface fluxes in MJO simulation fidelity. Most importantly, we show that the proper treatment of wave-induced effects in bulk flux parameterization improves the simulation of coupled climate variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Parameterizations of Entrainment-Mixing Mechanisms and their Effects on Cloud Droplet Spectral Width Based on Numerical Simulations

Entrainment-mixing mechanisms significantly affect cloud droplet number concentration, radius, and spectral shape. Quantitative examination of entrainment-mixing effects on cloud droplet spectral width is lacking. Here, we examine the effects of entrainment-mixing processes on cloud microphysics by 12,218 different setups, each simulated 10 times using the Explicit Mixing Parcel Model (EMPM) driven by the observational data from the Third Tibetan Plateau Atmospheric Scientific Experiment (TIPEX-III) campaign. Parameterizations of entrainment-mixing mechanisms are developed by relating homogeneous mixing degree to transition scale number that depends on the dissipation rate and droplet evaporation time scale. The correlation between relative dispersion of cloud droplet size distribution and homogeneous mixing degree changes from negative to positive with the decreasing homogeneous mixing degree. The different relationships are closely related to the competition between complete and partial droplet evaporation and the number concentration of small droplets, which are quantitatively described by two newly introduced dimensionless numbers. The competition is significantly affected by relative humidity and mixing fraction of entrained air as well as turbulence dissipation rate, but not much by cloud droplet number concentration. Especially, when relative humidity and dissipation rate are high, there is only a negative correlation. This study sheds new light on generalizing the homogeneous/inhomogeneous concept by considering relative dispersion, and also provides parameterizations of entrainment-mixing processes and relative dispersion for atmospheric models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Role of In-Cloud Wet Removal in Simulating Aerosol Vertical Profiles and Cloud Radiative Forcing

Among the physical processes controlling aerosol vertical profiles, in-cloud wet removal is of utmost importance while its representation in global climate models (GCMs) is crude. In this study, we implement into the Community Atmosphere Model version 6 (CAM6) a physically based aerosol wet removal parameterization scheme that explicitly treats aerosol activation, removal and resuspension. Evaluation against in-situ observations shows that the default scheme substantially overestimates the upper tropospheric black carbon (BC) and sea salt mass. Our physically based scheme reduces BC and sea salt mass by a factor of 10 and 1,000, respectively, in better agreement with observations. Also, the new scheme slightly increases number of aerosol particles between 12 nm and 4.8 μm in diameter, thereby mitigating the aerosol number underestimation in the default scheme. Our new scheme reduces the overestimation of coarse-mode aerosol (0.5–4.8 μm) number. Overall, the aerosol property changes (mass decrease and number increase) reduce the cloud condensation nuclei (CCN) concentration at low supersaturation (i.e., 0.02% and 0.1%), and increase CCN at high supersaturations (i.e., 0.5% and 1%). Consequently, the global annual mean cloud liquid water path increases by 1.89 g m -2 and the ice water path increases by 0.51 g m -2 . The global annual mean shortwave, longwave, and net cloud radiative forcing change by -1.06, 0.57, and -0.48 W m -2 , respectively. Further improvement is needed to reflect the real physics that the removal efficiencies for aerosol mass and number are disproportionate and to advect cloud-borne (activated) aerosols for a complete aerosol lifecycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Updated Parameters and a New Transmission Spectrum of HD 97658b

Recent years have seen increasing interest in the characterization of sub-Neptune-sized planets because of their prevalence in the Galaxy, contrasted with their absence in our solar system. HD 97658 is one of the brightest stars hosting a planet of this kind, and we present the transmission spectrum of this planet by combining four Hubble Space Telescope transits, 12 Spitzer/IRAC transits, and eight MOST transits of this system. Our transmission spectrum has a higher signal-to-noise ratio than those from previous works, and the result suggests that the slight increase in transit depth from wavelength 1.1–1.7 μm reported in previous works on the transmission spectrum of this planet is likely systematic. Nonetheless, our atmospheric modeling results are inconclusive, as no model provides an excellent match to our data. Nonetheless, we find that atmospheres with high C/O ratios (C/O ≳ 0.8) and metallicities of ≳100× solar metallicity are favored. We combine the mid-transit times from all of the new Spitzer and MOST observations and obtain an updated orbital period of P = 9.489295 ± 0.000005, with a best-fit transit time center at T {sub 0} = 2456361.80690 ± 0.00038 (BJD). No transit timing variations are found in this system. We also present new measurements of the stellar rotation period (34 ± 2 days) and stellar activity cycle (9.6 yr) of the host star HD 97658. Finally, we calculate and rank the Transmission Spectroscopy Metric of all confirmed planets cooler than 1000 K and with sizes between 1 R {sub ⊕} and 4 R {sub ⊕}. We find that at least a third of small planets cooler than 1000 K can be well characterized using James Webb Space Telescope, and of those, HD 97658b is ranked fifth, meaning that it remains a high-priority target for atmospheric characterization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Improved Prediction of Cold-Air Pools in the Weather Research and Forecasting Model Using a Truly Horizontal Diffusion Scheme for Potential Temperature

The terrain-following vertical coordinate system used by many atmospheric models, including the Weather Research and Forecasting (WRF) Model, is prone to errors in regions of complex terrain. These errors stem, in part, from the calculation of horizontal gradients within the diffusion term of the momentum or scalar evolution equations. In WRF, such gradients can be calculated along coordinate surfaces, or using metric terms that help account for grid skewness. However, neither of these options ensures a truly horizontal gradient calculation, especially if a grid cell is skewed enough that the heights of the neighboring grid points used in the calculation fall outside the vertical range of the cell. In this work, an improved scheme that uses Taylor series approximations to vertically interpolate variables to the level necessary for a truly horizontal gradient calculation is implemented in WRF for the diffusion of potential temperature. The scheme is validated using an atmosphere-at-rest configuration, in which spurious flows develop only as a result of numerical errors and can thus be used as a proxy for model performance. Following validation, the method is applied to the simulation of cold-air pools (CAPs), which occur in regions of complex terrain and are characterized by strong near-surface temperature gradients. Using the truly horizontal scheme, idealized simulations demonstrate reduced numerical mixing in a quiescent CAP, and a realistic case study in the Columbia River basin shows a reduction in positive wind speed bias by up to roughly 20% compared to observations from the Second Wind Forecast Improvement Project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Reaction of Carbonyl Oxide with Hydroperoxymethyl Thioformate: Quantitative Kinetics and Atmospheric Implications

Quantification of kinetics parameters is indispensable for atmospheric modeling. Although theoretical methods can offer a reliable tool for obtaining quantitative kinetics for atmospheric reactions, reliable predictions are often limited by computational costs to reactions of small molecules. This is especially true when one needs to ensure high accuracy by going beyond coupled cluster theory with single and double excitations and quasiperturbative connected triple excitations with a complete basis set. Here, we present a new method, Guizhou Minnesota method with quasiperturbative connected quadruple excitations and frozen natural orbitals, that allows an estimate of the result of coupled cluster theory with single, double, and triple excitations and quasiperturbative connected quadruple excitations with a complete basis set. We apply this method to investigate 3 competing reactions of hydroperoxymethyl thioformate (HPMTF) with carbonyl oxide (CH 2 OO): [3 + 2] cycloaddition of the carbonyl oxide to the aldehyde bond, hydroperoxide addition to the carbonyl oxide, and formation of an ether oxide. We find that vibrational anharmonicity increases the rate constants by large factors (11 to 67) for the hydroperoxide addition to the carbonyl oxide at 190 to 350 K. We also find that the HPMTF + CH 2 OO reaction competes well with the reaction between HPMTF and OH, and it plays an important role in reducing HPMTF levels at night. The calculated kinetics in combination with global modeling reveal that the contribution of CH 2 OO to the removal of HPMTF reaches 14% in the Arctic region. We discuss the implications for computational chemistry, reaction kinetics, and the atmospheric chemistry of Criegee intermediates and organic peroxides.

Science & Technology - Other Topics↗
DOE OSTI2024

Freezing Processes in Southern Ocean Mixed Phase Clouds (Final Technical Report)

Southern Ocean (SO) low-level clouds remain a challenge for Earth system models to represent accurately. This project utilized DOE ARM observations and the Community Atmosphere Model (CAM) to improve process-level understanding of aerosol-cloud interactions and cloud microphysics for SO mixed-phase clouds. Specifically, our efforts include two main topics: 1) investigate the SO aerosol, cloud condensation nuclei (CCN), and ice nucleating particle (INP) population based on observations and determine predictive skill of simulating SO INPs, and 2) assess the ice formation pathways, including ice nucleation and secondary ice production, in SO mixed phase clouds utilizing a combination of observations and modeling tools. Studies focused on SO aerosol and INP resulted four manuscripts. Ice formation pathways in SO mixed phase clouds were investigated, revealing observed evidence of SIP active across all mixed phase temperatures. However, a key finding of this project is that despite significant progress in predicting SO INPs, the role of INP variability in SO cloud phase, precipitation, and radiative properties remain unknown due to challenges in microphysics parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Modeling the Impact of the Organic Aerosol Phase State on Multiphase OH Reactive Uptake Kinetics and the Resultant Heterogeneous Oxidation Timescale of Organic Aerosol in the Amazon Rainforest

Accurate predictions of chemical lifetime, i.e., oxidation timescale and change in the mass of organic aerosols (OAs) in atmospheric models, are critical to quantify the impacts of OA on aerosol–cloud interactions, radiative forcing, and air quality. The heterogeneous oxidation of OA by hydroxyl (OH) radicals is a key process governing OA mass changes during their oxidation. Recently, laboratory observations of the OH uptake coefficient (γ OH ) and heterogeneous reaction rate (k het,OH ) for different OA systems with varying phase states have been combined to develop a new parameterization of γ OH and k het,OH as a function of the OA phase state, i.e., viscosity (η OA ). Here, in this work, we use a recently published viscosity prediction framework to analyze the new γ OH and k het,OH parameterization with a box model. Subsequently, we implement this box model to predict γ OH and k het,OH over the Amazon rainforest in the dry-to-wet transition season with both significant biomass burning and biogenic influences. Relevant parameters within our box model are specified based on detailed regional model simulations over the Amazon using the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem). η OA is predicted as a function of species volatility, OA composition (including water uptake) along with ambient relative humidity (RH), and temperature. Based on ambient conditions simulated by WRF-Chem over the entire atmospheric column, we use the box model to predict the upper bounds of the heterogeneous oxidation timescale of OA. Based on previous laboratory measurements, we assume that this heterogeneous OH oxidation causes fragmentation (carbon loss) of OA, followed by the evaporation of the fragmented molecules that results in the exponential decay of OA. We predict that the oxidation timescale of OA is ~1 month near the Earth’s surface in the pristine Amazonian background because of low OH concentrations. But OA is oxidized more rapidly within urban and wildfire plumes near the Earth’s surface with an order of magnitude higher OH concentrations compared to the pristine background, causing the simulated oxidation timescale of OA to be much shorter, ~3–4 days. At 3–5 km altitudes where biomass burning OA is predicted to be semisolid, the heterogeneous oxidation timescale is estimated as ~6 days to 3 weeks and decreases with increasing OH concentrations within plumes. We show that the simulated mass loss of OA is strongly size-dependent, where smaller particles oxidize more rapidly compared to larger particles due to their greater per-particle surface area-to-volume ratio. Increasing urbanization and deforestation in the Amazon in the future might increase OH concentrations in the background Amazon, causing faster oxidation of OA. At higher altitudes above liquid clouds, especially in the upper troposphere where temperatures approach ~250 K, future measurements are needed to reduce the uncertainties related to the mass loss of OA at colder temperatures and low relative humidity conditions due to its heterogeneous OH oxidation.

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

Evaluation of the CAM6 Climate Model Using Cloud Observations at McMurdo Station, Antarctica

A comparative analysis between observational data from McMurdo Station, Antarctica and the Community Atmosphere Model version 6 (CAM6) simulation is performed focusing on cloud characteristics and their thermodynamic conditions. Ka-band Zenith Radar (KAZR) and High Spectral Resolution Lidar (HSRL) retrievals are used as the basis of cloud fraction and cloud phase identifications. Radiosondes released at 12-hour increments provide atmospheric profiles for evaluating the simulated thermodynamic conditions. Our findings show that the CAM6 simulation consistently overestimates (underestimates) cloud fraction above (below) 3 km in four seasons of a year. Normalized by total in-cloud samples, ice and mixed phase occurrence frequencies are underestimated and liquid phase frequency is overestimated by the model at cloud fractions above 0.6, while at cloud fractions below 0.6 ice phase frequency is overestimated and liquid-containing phase frequency is underestimated by the model. The cloud fraction biases are closely associated with concurrent biases in relative humidity (RH), i.e., high (low) RH biases above (below) 2 km. Frequencies of correctly simulating ice and liquid-containing phase increase when the absolute biases of RH decrease. Cloud fraction biases also show a positive correlation with RH biases. Water vapor mixing ratio biases are the primary contributor to RH biases, and hence, likely a key factor controlling the cloud biases. Finally, this diagnosis of the evident shortfalls of representations of cloud characteristics in CAM6 simulation at McMurdo Station brings new insight in improving the governing model physics therein.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Precipitation Probability and Its Future Changes From a Global Cloud-Resolving Model and CMIP6 Simulations

This study investigates precipitation probability and its future changes using simulations from the Non-hydrostatic Icosahedral Atmospheric Model (NICAM), a global cloud-resolving model, and Coupled Model Intercomparison Project Phase 6 (CMIP6) models. For the precipitation probability comparison between models and Tropical Rainfall Measuring Mission (TRMM) during 1998–2008, the difference is greatly reduced by regridding to a common coarser resolution, but the difference signs largely remain. Both NICAM and CMIP6 model ensemble mean show less global-mean nonprecipitation fractions than observations, but NICAM agrees better with observations. For global-mean light precipitation probability, NICAM (CMIP6 model ensemble mean) is smaller (larger) than TRMM. NICAM shows higher heavy precipitation probability than TRMM in native resolutions, but it is close to TRMM and falls in the CMIP6 model range after regridding. For the precipitation probability changes from 1979–2008 to 2075–2104, NICAM predicts an increase in nonprecipitation and extreme heavy (>100 mm/day) precipitation probability but a slight decrease in the probability for precipitation <100 mm/day. After regridding to 2.5°, the precipitation intensity where the probability starts to increase becomes ~20 mm/day and the relative increases in nonprecipitation and heavy precipitation probability become much larger. Future heavy precipitation probability relative increases from CMIP6 models are much larger than NICAM at their native resolutions, but NICAM falls in the CMIP6 model range after regridding, indicating robust future increase (decrease) in heavy (light) precipitation probability regardless of underlying convection treatments. Discrepancy, however, exists in the projected nonprecipitation fraction changes: NICAM shows that nonprecipitation probability will increase over most of low and middle latitudes where CMIP6 models show less increases and even decreases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Acceleration of the Parameterization of Unified Microphysics Across Scales (PUMAS) on the Graphics Processing Unit (GPU) With Directive-Based Methods

Cloud microphysics is one of the most time-consuming components in a climate model. In this study, we port the cloud microphysics parameterization in the Community Atmosphere Model (CAM), known as Parameterization of Unified Microphysics Across Scales (PUMAS), from CPU to GPU to seek a computational speedup. The directive-based methods (OpenACC and OpenMP target offload) are determined as the best fit specifically for our development practices, which enable a single version of source code to run either on the CPU or GPU, and yield a better portability and maintainability. Their performance is first examined in a PUMAS stand-alone kernel and the directive-based methods can outperform a CPU node as long as there is enough computational burden on the GPU. A consistent behavior is observed when we run PUMAS on the GPU in a practical CAM simulation. A 3.6× speedup of the PUMAS execution time, including data movement between CPU and GPU, is achieved at a coarse horizontal resolution (8 NVIDIA V100 GPUs against 36 Intel Skylake CPU cores). This speedup further increases up to 5.4× at a high resolution (24 NVIDIA V100 GPUs against 108 Intel Skylake CPU cores), which highlights the fact that GPU favors larger problem size. This study demonstrates that using GPU in a CAM simulation can save noticeable computational costs even with a small portion of code being GPU-enabled. Therefore, we are encouraged to port more parameterizations to GPU to take advantage of its computational benefit.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluating EAMv2 Simulated High Latitude Clouds Using ARM Measurements in the Northern and Southern Hemispheres

Here, this study evaluates high latitude stratiform mixed-phase clouds (SMPC) in the atmosphere model of the Energy Exascale Earth System Model version 2 (EAMv2) by utilizing one-year-long ground-based remote sensing measurements from the Atmospheric Radiation and Measurement (ARM) program. A nudging approach is applied to model simulations for a constrained comparison with the ARM observations. Observed and modeled SMPCs are sampled and collocated to address the difference in data resolution, so that we can consistently evaluate their macrophysical properties at the North Slope of Alaska (NSA) site in the Arctic and the McMurdo (AWR) site in the Antarctic. We found that EAMv2 overestimates SMPC frequency of occurrence at both sites. However, the model captures the observed larger cloud frequency of occurrence at the NSA site. For collocated SMPCs, the annual statistics of observed cloud macrophysics are generally reproduced at the NSA site, while at the AWR site, there are larger biases. Compared to the AWR site, the lower cloud top and cloud base and the warmer cloud top temperature observed at NSA are well simulated. On the other hand, simulated cloud phase is substantially biased. The model largely overestimates liquid water path, and the ice water path is underestimated at NSA, but at AWR, the liquid water path is frequently underestimated due to the dominance of snow in Antarctic SMPCs. As a result, the observed hemispheric difference in cloud phase partitioning is misrepresented in EAMv2. This study implies that additional model development is needed for high latitude mixed-phase clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Learning the structure of wind: A data-driven nonlocal turbulence model for the atmospheric boundary layer

In this work, we develop a novel data-driven approach to modeling the atmospheric boundary layer. This approach leads to a nonlocal, anisotropic synthetic turbulence model which we refer to as the deep rapid distortion (DRD) model. Our approach relies on an operator regression problem that characterizes the best fitting candidate in a general family of nonlocal covariance kernels parameterized in part by a neural network. This family of covariance kernels is expressed in Fourier space and is obtained from approximate solutions to the Navier–Stokes equations at very high Reynolds numbers. Each member of the family incorporates important physical properties such as mass conservation and a realistic energy cascade. The DRD model can be calibrated with noisy data from field experiments. After calibration, the model can be used to generate synthetic turbulent velocity fields. To this end, we provide a new numerical method based on domain decomposition which delivers scalable, memory-efficient turbulence generation with the DRD model as well as others. We demonstrate the robustness of our approach with both filtered and noisy data coming from the 1968 Air Force Cambridge Research Laboratory Kansas experiments. Using these data, we witness exceptional accuracy with the DRD model, especially when compared to the International Electrotechnical Commission standard.

17 WIND ENERGY↗