Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cloud parameterizations”

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.

807 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Mechanical and Thermal Forcing for Upslope Flows and Cumulus Convection over the Sierras de Córdoba

Abstract The upslope flow processes affecting the vertical extent of orographic cumulus convection are examined using observations from the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. Specifically, clear air returns from the U.S. Department of Energy (DOE) second-generation C-band scanning Atmospheric Radiation Measurement (ARM) precipitation radar (CSAPR2) are used to characterize the structure and variability of the ridge-normal (i.e., up/downslope) flow components, which transport mass to the crest of Argentina’s Sierras de Córdoba and contribute to convective initiation. Data are compiled for the entire CACTI period (October–April), including days with clear skies, shallow cumuli, cumulus congestus, and deep convection. To examine shared variability among >70 000 radar scans, we use (i) a principal component analysis (PCA) to isolate modes of variability in the upslope flow and (ii) composite analysis based on convective outcomes, determined from GOES-16 satellite observations. These data are contextualized with observed surface sensible heat fluxes, thermodynamic profiles, and synoptic-scale analysis. Results indicate distinct thermally and mechanically forced upslope flow modes, modulated by diurnal heating and synoptic-scale variations, respectively. In some instances, there is a superposition of thermal and mechanical forcing, yielding either deeper or shallower upslope flow. The composite analyses based on satellite data show that successively deeper convective outcomes are associated with successively deeper upslope flow layers that more readily transport mass to the ridge crest in conjunction with lower lifting condensation levels, facilitating convective initiation. These results help to isolate the forcing mechanisms for orographic convection and thus provide a foundation for parameterizing orographic convective processes in coarse resolution models.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Impacts of Suppressing Excessive Light Rain on Aerosol Radiative Effects and Health Risks

Global climate models (GCMs) have been used widely to study radiative forcing and health risks of aerosols. A recent study using two GCMs found that light rain plays a dominant role in controlling aerosol loading. However, “too much light rain and too little heavy rain” is a longstanding bias in GCMs. It is unclear how much light rain affects aerosol-cloud-radiation interactions and health risks from air pollution. Here, in this work, we show that, with the correction of the rainfall intensity spectrum in the National Center for Atmospheric Research Community Atmosphere Model version 5.3 by introducing a stochastic deep convection scheme, the reduced frequency of light rain (1–20 mm d -1 ) results in changes of aerosol direct radiative effects (DRE) of up to -0.5 ± 0.03 W/m 2 and aerosol cloud radiative effects (CRE) of up to -0.9 ± 0.03 W/m 2 . The total (CRE + DRE) radiative effects of light rain-mediated aerosol changes exceed the present-day anthropogenic forcing of aerosols relative to preindustrial levels from the Coupled Model Intercomparison Project (CMIP5&6) models. However, the correction of the rainfall intensity spectrum has little effect on anthropogenic aerosol forcing (defined as the radiative perturbation due to changes in aerosol concentrations between the industrial era and preindustrial levels). Due to increased exposure to fine particulates (PM 2.5 ), the estimated global total premature mortality is much higher than previously estimated, by 300,000 ± 60,000 deaths per year, and is more severe in populous regions such as India and China. The findings in this study highlight the need to understand uncertainties in radiative effects and health risks of aerosols due to simulation biases of precipitation in GCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

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

The complexity of radiative transfer, its importance to the exchange of energy in the climate system, and its high computational cost establishes importance of an accurate and efficient radiative transfer parameterization for climate simulation. Previous work by the proposing team at AER led to the development of the radiation code, RRTMG, which has been widely accepted by the global modeling community as a fast and accurate advancement over the previous generation of radiation codes. It has been in use in the NCAR CESM for many years, and it has been implemented in the initial version the DOE E3SM model. However, its computational cost remains high relative to other components in part due to its complexity and to its inefficient use of modern optimization strategies, and this project helped address this limitation for the code’s application in E3SM. Under other funding, the Investigators of this project led an effort to develop a high-performance broadband radiation code, called RTE+RRTMGP, which is a completely restructured code that will take advantage of modern computational capabilities to enhance its performance while retaining the strengths and accuracy of the original code. Designed to perform over a range of computer architectures, RTE+RRTMGP makes extensive use of Fortran 2003 features to improve both its efficiency and its use of memory. RTE+RRTMGP is expected to be adopted in the next generation of RTE+RRTMGP. This project allowed for advancements to the code’s computational capabilities and adding new and enhanced features. This included revising RTE+RRTMGP to run on GPU processors, analysis of and improvements to the code’s timing, modifying the gas optics in RRTMGP to increase the code’s accuracy, extensive validation of computed fluxes, heating rates and forcings, generating a cloud optical property and vertical sampling capabilities for RTE+RRTMGP, implementing a fast and accurate longwave scattering capability, and groundwork for the inclusion of a capability to specify solar variability. Many of the accomplishment in this project necessitated significant collaboration with the E3SM development team. The result of this project was optimization of a key physical component (radiative transfer calculations) of E3SM, directly supporting E3SM’s overarching global modeling objectives. More broadly, this project provided overall advancements in the use of radiative transfer calculations in atmospheric modeling and simulation, particularly for climate.

58 GEOSCIENCES↗
DOE OSTI2020

Identifying a regional aerosol baseline in the eastern North Atlantic using collocated measurements and a mathematical algorithm to mask high-submicron-number-concentration aerosol events

Abstract. High-time-resolution measurements of in situ aerosol and cloud properties provide the ability to study regional atmospheric processes that occur on timescales of minutes to hours. However, one limitation to this approach is that continuous measurements often include periods when the data collected are not representative of the regional aerosol. Even at remote locations, submicron aerosols are pervasive in the ambient atmosphere with many sources. Therefore, periods dominated by local aerosol should be identified before conducting subsequent analyses to understand aerosol regional processes and aerosol–cloud interactions. Here, we present a novel method to validate the identification of regional baseline aerosol data by applying a mathematical algorithm to the data collected at the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) user facility in the eastern North Atlantic (ENA). The ENA central facility (C1) includes an aerosol observing system (AOS) for the measurement of aerosol physical, optical, and chemical properties at time resolutions from seconds to minutes. A second temporary supplementary facility (S1), located ∼0.75 km from C1, was deployed for ∼1 year during the Aerosol and Cloud Experiments (ACE-ENA) campaign in 2017. First, we investigate the local aerosol at both locations. We associate periods of high submicron number concentration (Ntot) in the fine-mode condensation particle counter (CPC) and size distributions from the Ultra-High Sensitivity Aerosol Spectrometer (UHSAS) as a function of wind direction using a meteorology sensor with local sources. Elevated concentrations of Aitken-mode (< 100 nm diameter) particles were observed in correspondence with the wind directions associated with airport operations. At ENA, the Graciosa Airport and its associated activities were found to be the main sources of high-concentration aerosol events at ENA, causing peaks in 1 min Ntot that exceeded 8000 and 10 000 cm−3 at C1, in summer and winter, respectively, and 5000 cm−3 at S1 in summer. Periods with high Ntot not associated with these wind directions were also observed. As a result, the diverse local sources at ENA yielded a poor relationship between Ntot measurements collected at C1 and S1 (R2=0.03 with a slope =0.05±0.001). As a first approach to mask these events, the time periods when the wind direction was associated with the airport operations (west to northwest and southeast to south at C1 and east to south at S1) were applied. The meteorological masks removed 38.9 % of the data at C1 and 43.4 % at S1, and they did not significantly improve the relationship between the two sites (R2=0.18 with a slope =0.06±0.001). Due to the complexity of high-Ntot events observed at ENA, we develop and validate a mathematical ENA Aerosol Mask (ENA-AM) to identify high-Ntot events using 1 min resolution data from the AOS CPC at C1 and S1. After its parameterization and application, ENA-AM generated a high correlation between Ntot in the summer at C1 and S1 (R2=0.87 with a slope =0.84±0.001). We identified the regional baseline at ENA to be 428±228 cm−3 in the summer and 346±223 cm−3 in the winter. Lastly, we compared masked measurements from the AOS with the ARM Aerial Facility (AAF) during flights over C1 in the summer to understand submicron aerosol vertical mixing over C1. The high correlation (R2=0.71 with a slope of 1.04±0.01) observed between C1 and the AAF Ntot collected within an area of 10 km surrounding ENA and at altitudes < 500 m indicated that the submicron aerosol at ENA was well mixed within the first 500 m of the marine boundary layer during the month of July during ACE-ENA. Our novel method for determining a regional aerosol baseline at ENA can be applied to other time periods and at other locations with validation by a secondary site or additional collocated measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Microphysical and Dynamical Effects of Mixed-Phase Hydrometeors in Convective Storms Using a Bin Microphysics Model: Melting

Abstract The dynamics of convective systems are inherently linked to microphysical processes through phase changes that result in warming or cooling. This is especially true of near-surface cooling via evaporation and melting of falling hydrometeors. In most numerical simulations, the melting of frozen hydrometeors (e.g., hail, graupel, snow) is computed within parameterized bulk microphysics schemes, many of which lack the ability to accurately represent mixed-phase hydrometeors (i.e., partially melted ice), which can affect hydrometeor sedimentation, melting, and evaporation of shed drops. To better understand the microphysical and dynamical effects of melting in convective storms, a bin microphysics scheme was used in the Weather Research and Forecasting Model for two idealized cases: a supercell storm and a squall line. Physically based predicted liquid fraction, instantaneous melting, and instantaneous shedding schemes were used to examine the role and importance of melting hydrometeors for these two storm modes. The results suggest that the amount of precipitation is dependent on the representation of melting. Moreover, the dynamic and thermodynamic characteristics of the simulated storms are found to differ substantially between the melting scenarios, resulting in varied storm system evolution; these differences are found to be dependent on the ambient aerosol concentration, although the differences induced by changing the representation of melting generally outweigh those of changing the aerosol loading. The results highlight the large role of melting in convective storm characteristics and suggest that further model improvements are needed in the near future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

The Fifth International Workshop on Ice Nucleation phase 2 (FIN-02): laboratory intercomparison of ice nucleation measurements

Abstract. The second phase of the Fifth International Ice Nucleation Workshop (FIN-02)involved the gathering of a large number of researchers at the KarlsruheInstitute of Technology's Aerosol Interactions and Dynamics of the Atmosphere(AIDA) facility to promote characterization and understanding of icenucleation measurements made by a variety of methods used worldwide.Compared to the previous workshop in 2007, participation was doubled,reflecting a vibrant research area. Experimental methods involved sampling ofaerosol particles by direct processing ice nucleation measuring systems fromthe same volume of air in separate experiments using different ice nucleatingparticle (INP) types, and collections of aerosol particle samples ontofilters or into liquid for sharing amongst measurement techniques thatpost-process these samples. In this manner, any errors introduced bydifferences in generation methods when samples are shared across laboratorieswere mitigated. Furthermore, as much as possible, aerosol particle sizedistribution was controlled so that the size limitations of different methodswere minimized. The results presented here use data from the workshop toassess the comparability of immersion freezing measurement methods activatingINPs in bulk suspensions, methods that activate INPs in condensation and/orimmersion freezing modes as single particles on a substrate, continuous flowdiffusion chambers (CFDCs) directly sampling and processing particles wellabove water saturation to maximize immersion and subsequent freezing ofaerosol particles, and expansion cloud chamber simulations in which liquidcloud droplets were first activated on aerosol particles prior to freezing.The AIDA expansion chamber measurements are expected to be the closestrepresentation to INP activation in atmospheric cloud parcels in thesecomparisons, due to exposing particles freely to adiabatic cooling. The different particle types used as INPs included the minerals illite NX andpotassium feldspar (K-feldspar), two natural soil dusts representative of arable sandy loam(Argentina) and highly erodible sandy dryland (Tunisia) soils, respectively,and a bacterial INP (Snomax ® ). Consideredtogether, the agreement among post-processed immersion freezing measurementsof the numbers and fractions of particles active at different temperaturesfollowing bulk collection of particles into liquid was excellent, withpossible temperature uncertainties inferred to be a key factor in determiningINP uncertainties. Collection onto filters for rinsing versus directly intoliquid in impingers made little difference. For methods that activatedcollected single particles on a substrate at a controlled humidity at orabove water saturation, agreement with immersion freezing methods was good inmost cases, but was biased low in a few others for reasons that have not beenresolved, but could relate to water vapor competition effects. AmongstCFDC-style instruments, various factors requiring (variable) highersupersaturations to achieve equivalent immersion freezing activation dominatethe uncertainty between these measurements, and for comparison with bulkimmersion freezing methods. When operated above water saturation to includeassessment of immersion freezing, CFDC measurements often measured at orabove the upper bound of immersion freezing device measurements, but oftenunderestimated INP concentration in comparison to an immersion freezingmethod that first activates all particles into liquid droplets prior tocooling (the PIMCA-PINC device, or Portable Immersion Mode Cooling chAmber–Portable Ice Nucleation Chamber), and typically slightly underestimated INPnumber concentrations in comparison to cloud parcel expansions in the AIDAchamber; this can be largely mitigated when it is possible to raise therelative humidity to sufficiently high values in the CFDCs, although this isnot always possible operationally. Correspondence of measurements of INPs among direct sampling andpost-processing systems varied depending on the INP type. Agreement was bestfor Snomax ® particles in the temperature regimecolder than −10 ∘ C, where their ice nucleation activity is nearlymaximized and changes very little with temperature. At temperatures warmer than −10 ∘ C, Snomax ® INP measurements (allvia freezing of suspensions) demonstrated discrepancies consistent withprevious reports of the instability of its protein aggregates that appear tomake it less suitable as a calibration INP at these temperatures. ForArgentinian soil dust particles, there was excellent agreement across allmeasurement methods; measures ranged within 1 order of magnitude for INPnumber concentrations, active fractions and calculated active site densitiesover a 25 to 30 ∘ C range and 5 to 8 orders of correspondingmagnitude change in number concentrations. This was also the case for alltemperatures warmer than −25 ∘ C in Tunisian dust experiments. Incontrast, discrepancies in measurements of INP concentrations or active sitedensities that exceeded 2 orders of magnitude across a broad range of temperaturemeasurements found at temperatures warmer than −25 ∘ C in a previous study werereplicated for illite NX. Discrepancies also exceeded 2 orders of magnitude attemperatures of −20 to −25 ∘ C for potassium feldspar (K-feldspar), but these coincidedwith the range of temperatures at which INP concentrations increase rapidly atapproximately an order of magnitude per 2 ∘ C cooling forK-feldspar. These few discrepancies did not outweigh the overall positive outcomes of theworkshop activity, nor the future utility of this data set or future similarefforts for resolving remaining measurement issues. Measurements of the samematerials were repeatable over the time of the workshop and demonstratedstrong consistency with prior studies, as reflected by agreement of databroadly with parameterizations of different specific or general (e.g., soildust) aerosol types. The divergent measurements of the INP activity of illiteNX by direct versus post-processing methods were not repeated for otherparticle types, and the Snomax ® datademonstrated that, at least for a biological INP type, there is no expectedmeasurement bias between bulk collection and direct immediately processedfreezing methods to as warm as −10 ∘ C. Since particle size rangeswere limited for this workshop, it can be expected that for atmosphericpopulations of INPs, measurement discrepancies will appear due to thedifferent capabilities of methods for sampling the full aerosol sizedistribution, or due to limitations on achieving sufficient watersupersaturations to fully capture immersion freezing in direct processinginstruments. Overall, this workshop presents an improved picture of presentcapabilities for measuring INPs than in past workshops, and providesdirection toward addressing remaining measurement issues.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Chemical composition, microstructure, and hygroscopic properties of aerosol particles at the Zotino Tall Tower Observatory (ZOTTO), Siberia, during a summer campaign

In this study we describe the hygroscopic properties of accumulation- and coarse-mode aerosol particles sampled at the Zotino Tall Tower Observatory (ZOTTO) in Central Siberia (61° N; 89° E) from 16 to 21 June 2013. The hygroscopic growth measurements were supplemented with chemical analyses of the samples, including inorganic ions and organic/elemental carbon. In addition, the microstructure and chemical composition of aerosol particles were analyzed by X-ray micro-spectroscopy (STXM-NEXAFS) and transmission electron microscopy (TEM). A mass closure analysis indicates that organic carbon accounted for 61 and 38% of PM in the accumulation mode and coarse mode, respectively. The water soluble fraction of organic matter was estimated to be 52 and 8% of PM in these modes. Sulfate, predominantly in the form of ammoniated sulfate, was the dominant inorganic component in both size modes: ~ 34% in the accumulation vs. ~ 47% in the coarse mode. The hygroscopic growth measurements were conducted with a filter-based differential hygroscopicity analyzer (FDHA) over the range of 5–99.4% RH in the hydration and dehydration operation modes. The FDHA study indicates that both accumulation and coarse modes exhibit pronounced water uptake approximately at the same RH, starting at ~ 70%, while efflorescence occurred at different humidities, i.e., at ~ 35% RH for submicron particles vs. ~ 50% RH for supermicron particles. This ~ 15% RH difference was attributed to higher content of organic material in the submicron particles, which suppresses water release in the dehydration experiments. The kappa mass interaction model (KIM) was applied to characterize and parameterize non-ideal solution behavior and concentration-dependent water uptake by atmospheric aerosol samples in the 5–99.4% RH range. Based on KIM, the volume-based hygroscopicity parameter, &kappa; v , was calculated. The &kappa; v, ws value related to the water soluble (ws) fraction was estimated to be ~ 0.15 for the accumulation mode and ~ 0.36 for the coarse mode, respectively. The obtained &kappa; v, ws for the accumulation mode is in good agreement with earlier data reported for remote sites in the Amazon rain forest (&kappa; v &approx; 0.15) and a Colorado boreal forest (&kappa; v &approx; 0.16). We used the Zdanovskii–Stokes–Robinson (ZSR) mixing rule to predict the chemical composition dependent hygroscopicity, &kappa; v, p . The obtained &kappa; v, p values overestimate the experimental FDHA-KIM-derived &kappa; v, ws by factors of 1.8 and 1.5 for the accumulation and coarse modes, respectively. This divergence can be partly explained by incomplete dissolution of the hygroscopic inorganic compounds resulting from kinetic limitations due to a sparingly soluble organic coating. The TEM and STXM-NEXAFS results indicate that aged submicron (>300 nm) and supermicron aerosol particles possess core-shell structures with an inorganic core, and are enriched in organic carbon at the mixed particle surface. The direct FDHA kinetic studies provide a bulk diffusion coefficient of water of ~ 10 &minus;12 cm 2 s &minus;1 indicating a semi-solid state of the organic-rich phase leading to kinetic limitations of water uptake and release during hydration and dehydration cycles. Overall the present ZOTTO data set, obtained in the growing season, has revealed a strong influence of organic carbon on the hygroscopic properties of the ambient aerosols. The sparingly soluble organic coating controls hygroscopic growth, phase transitions, and microstructural rearrangement processes. The observed kinetic limitations can strongly influence the outcome of experiments performed on multi-second time scales, such as the commonly applied HTDMA (Hygroscopicity Tandem Differential Mobility Analyzer) and CCNC (Cloud Condensation Nuclei Counter) measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

A comprehensive laboratory study on the immersion freezing behavior of illite NX particles. A comparison of 17 ice nucleation measurement techniques

Immersion freezing is the most relevant heterogeneous ice nucleation mechanism through which ice crystals are formed in mixed-phase clouds. In recent years, an increasing number of laboratory experiments utilizing a variety of instruments have examined immersion freezing activity of atmospherically relevant ice-nucleating particles. However, an intercomparison of these laboratory results is a difficult task because investigators have used different ice nucleation (IN) measurement methods to produce these results. A remaining challenge is to explore the sensitivity and accuracy of these techniques and to understand how the IN results are potentially influenced or biased by experimental parameters associated with these techniques. Within the framework of INUIT (Ice Nuclei Research Unit), we distributed an illite-rich sample (illite NX) as a representative surrogate for atmospheric mineral dust particles to investigators to perform immersion freezing experiments using different IN measurement methods and to obtain IN data as a function of particle concentration, temperature ( T ), cooling rate and nucleation time. A total of 17 measurement methods were involved in the data intercomparison. Experiments with seven instruments started with the test sample pre-suspended in water before cooling, while 10 other instruments employed water vapor condensation onto dry-dispersed particles followed by immersion freezing. The resulting comprehensive immersion freezing data set was evaluated using the ice nucleation active surface-site density, n s , to develop a representative n s ( T ) spectrum that spans a wide temperature range (-37 °C < T < -11 °C) and covers 9 orders of magnitude in n s . In general, the 17 immersion freezing measurement techniques deviate, within a range of about 8 °C in terms of temperature, by 3 orders of magnitude with respect to n s . In addition, we show evidence that the immersion freezing efficiency expressed in n s of illite NX particles is relatively independent of droplet size, particle mass in suspension, particle size and cooling rate during freezing. A strong temperature dependence and weak time and size dependence of the immersion freezing efficiency of illite-rich clay mineral particles enabled the n s parameterization solely as a function of temperature. We also characterized the n s ( T ) spectra and identified a section with a steep slope between -20 and -27 °C, where a large fraction of active sites of our test dust may trigger immersion freezing. This slope was followed by a region with a gentler slope at temperatures below -27 °C. While the agreement between different instruments was reasonable below ~ -27 °C, there seemed to be a different trend in the temperature-dependent ice nucleation activity from the suspension and dry-dispersed particle measurements for this mineral dust, in particular at higher temperatures. For instance, the ice nucleation activity expressed in n s was smaller for the average of the wet suspended samples and higher for the average of the dry-dispersed aerosol samples between about -27 and -18 °C. Only instruments making measurements with wet suspended samples were able to measure ice nucleation above -18 °C. A possible explanation for the deviation between -27 and -18 °C is discussed. Multiple exponential distribution fits in both linear and log space for both specific surface area-based n s ( T ) and geometric surface area-based n s ( T ) are provided. These new fits, constrained by using identical reference samples, will help to compare IN measurement methods that are not included in the present study and IN data from future IN instruments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Decentralized Approach for Modeling Organized Convection Based on Thermal Populations on Microgrids

Abstract In this study, a spectral model for convective transport is coupled to a thermal population model on a two‐dimensional horizontal “microgrid,” covering the typical gridbox size of general circulation models. The goal is to explore new ways of representing impacts of spatial organization in cumulus cloud fields. The thermals are considered the smallest building block of convection, with thermal life cycle and movement represented through binomial functions. Thermals interact through two simple rules, reflecting pulsating growth and environmental deformation. Long‐lived thermal clusters thus form on the microgrid, exhibiting scale growth and spacing that represent simple forms of spatial organization and memory. Size distributions of cluster number are diagnosed from the microgrid through an online clustering algorithm, and provided as input to a spectral multiplume eddy‐diffusivity mass flux scheme. This yields a decentralized transport system, in that the thermal clusters acting as independent but interacting nodes that carry information about spatial structure. The main objectives of this study are (a) to seek proof of concept of this approach, and (b) to gain insight into impacts of spatial organization on convective transport. Single‐column model experiments demonstrate satisfactory skill in reproducing two observed cases of continental shallow convection. Metrics expressing self‐organization and spatial organization match well with large‐eddy simulation results. We find that in this coupled system, spatial organization impacts convective transport primarily through the scale break in the size distribution of cluster number. The rooting of saturated plumes in the subcloud mixed layer plays a key role in this process.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Convection-Permitting Simulations With the E3SM Global Atmosphere Model

This paper describes the first implementation of the Δx = 3.25 km version of the Energy Exascale Earth System Model (E3SM) global atmosphere model and its behavior in a 40-day prescribed-sea-surface-temperature simulation (January 20 through February 28, 2020). This simulation was performed as part of the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) Phase 2 model intercomparison. Effective resolution is found to be ~ 6x the horizontal dynamics grid resolution despite using a coarser grid for physical parameterizations. Despite this new model being in an immature and untuned state, moving to 3.25 km grid spacing solves several long-standing problems with the E3SM model. In particular, Amazon precipitation is much more realistic, the frequency of light and heavy precipitation is improved, agreement between the simulated and observed diurnal cycle of tropical precipitation is excellent, and the vertical structure of tropical convection and coastal stratocumulus look good. In addition, the new model is able to capture the frequency and structure of important weather events (e.g., tropical cyclones, extratropical cyclones including atmospheric rivers, and cold air outbreaks). Interestingly, this model does not get rid of the erroneous southern branch of the intertropical convergence zone nor the tendency for strongest convection to occur over the Maritime Continent rather than the West Pacific, both of which are classic climate model biases. Several other problems with the simulation are identified, underscoring the fact that this model is a work in progress.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Process-Oriented Diagnostics: Principles, Practice, Community Development, and Common Standards

Process-oriented diagnostics (PODs) aim to provide feedback for model developers through model analysis based on physical hypotheses. However, the step from a diagnostic based on relationships among variables, even when hypothesis driven, to specific guidance for revising model formulation or parameterizations can be substantial. The POD may provide more information than a purely performance-based metric, but a gap between POD principles and providing actionable information for specific model revisions can remain. Furthermore, in coordinating diagnostics development, there is a trade-off between freedom for the developer, aiming to capture innovation, and near-term utility to the modeling center. Best practices that allow for the former, while conforming to specifications that aid the latter, are important for community diagnostics development that leads to tangible model improvements. Promising directions to close the gap between principles and practice include the interaction of PODs with perturbed physics experiments and with more quantitative process models as well as the inclusion of personnel from modeling centers in diagnostics development groups for immediate feedback during climate model revisions. Examples are provided, along with best-practice recommendations, based on practical experience from the NOAA Model Diagnostics Task Force (MDTF). Common standards for metrics and diagnostics that have arisen from a collaboration between the MDTF and the Department of Energy’s Coordinated Model Evaluation Capability are advocated as a means of uniting community diagnostics efforts.

58 GEOSCIENCES↗
DOE OSTI2021

Evaluating the Conservation of Energy Variables in Simulations of Deep Moist Convection

Abstract It is often assumed in parcel theory calculations, numerical models, and cumulus parameterizations that moist static energy (MSE) is adiabatically conserved. However, the adiabatic conservation of MSE is only approximate because of the assumption of hydrostatic balance. Two alternative variables are evaluated here: MSE − IB and MSE + KE, wherein IB is the path integral of buoyancy ( B ) and KE is kinetic energy. Both of these variables relax the hydrostatic assumption and are more precisely conserved than MSE. This article quantifies the errors that result from assuming that the aforementioned variables are conserved in large-eddy simulations (LES) of both disorganized and organized deep convection. Results show that both MSE − IB and MSE + KE better predict quantities along trajectories than MSE alone. MSE − IB is better conserved in isolated deep convection, whereas MSE − IB and MSE + KE perform comparably in squall-line simulations. These results are explained by differences between the pressure perturbation behavior of squall lines and isolated convection. Errors in updraft B diagnoses are universally minimized when MSE − IB is assumed to be adiabatically conserved, but only when moisture dependencies of heat capacity and temperature dependency of latent heating are accounted for. When less accurate latent heat and heat capacity formulae were used, MSE − IB yielded poorer B predictions than MSE due to compensating errors. Our results suggest that various applications would benefit from using either MSE − IB or MSE + KE instead of MSE with properly formulated heat capacities and latent heats.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

CCN activity and organic hygroscopicity of aerosols downwind of an urban region in central Amazonia: Seasonal and diel variations and impact of anthropogenic emissions

During the Observations and Modeling of the Green Ocean Amazon (GoAmazon2014/5) campaign, size-resolved cloud condensation nuclei (CCN) spectra were characterized at a research site (T3) 60km downwind of the city of Manaus, Brazil, in central Amazonia for one year (12 March 2014 to 3 March 2015). Particle hygroscopicity (κCCN) and mixing state were derived from the size-resolved CCN spectra, and the hygroscopicity of the organic component of the aerosol (κorg) was then calculated from κCCN and concurrent chemical composition measurements. The annual average κCCN increased from 0.13 at 75 nm to 0.17 at 171&thinsp;nm, and the increase was largely due to an increase in sulfate volume fraction. Also, during both wet and dry seasons, κCCN, κorg, and particle composition under background conditions exhibited essentially no diel variations. The constant κorg of ~0.15 is consistent with the largely uniform and high O:C value (~0.8), indicating that the aerosols under background conditions are dominated by the aged regional aerosol particles consisting of highly oxygenated organic compounds. For air masses strongly influenced by urban pollution and/or local biomass burning, lower values of κorg and organic O:C atomic ratio were observed during night, due to accumulation of freshly emitted particles, dominated by primary organic aerosol (POA) with low hygroscopicity, within a shallow nocturnal boundary layer. The O:C, κorg, and κCCN increased from the early morning hours and peaked around noon, driven by the formation and aging of secondary organic aerosol (SOA) and dilution of POA emissions into a deeper boundary layer, while the development of the boundary layer, which leads to mixing with aged particles from the residual layer aloft, likely also contributed to the increases. The hygroscopicities associated with individual organic factors, derived from PMF analysis of AMS spectra, were estimated through multi-variable linear regression. For the SOA factors, the variation of the κ value with O:C agrees well with the linear relationship reported from earlier laboratory studies of SOA hygroscopicity. On the other hand, the variation in O:C of ambient aerosol organics is largely driven by the variation in the volume fractions of POA and SOA factors, which have very different O:C values. As POA factors have hygroscopicity values well below the linear relationship between SOA hygroscopicity and O:C, mixtures with different POA and SOA fractions exhibit a steeper slope for the increase of κorg with O:C, as observed during this and earlier field studies. Our finding helps better understand and reconcile the differences in the relationships between κorg and O:C observed in laboratory and field studies, therefore providing a basis for improved parameterization in global models, especially in a tropical context.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

CCN activity and organic hygroscopicity of aerosols downwind of an urban region in central Amazonia: seasonal and diel variations and impact of anthropogenic emissions

Abstract. During the Observations and Modeling of the Green Ocean Amazon (GoAmazon2014/5) campaign, size-resolved cloud condensation nuclei (CCN) spectra were characterized at a research site (T3) 60 km downwind of the city of Manaus, Brazil, in central Amazonia for 1 year (12 March 2014 to 3 March 2015). Particle hygroscopicity (κCCN) and mixing state were derived from the size-resolved CCN spectra, and the hygroscopicity of the organic component of the aerosol (κorg) was then calculated from κCCN and concurrent chemical composition measurements. The annual average κCCN increased from 0.13 at 75 nm to 0.17 at 171 nm, and the increase was largely due to an increase in sulfate volume fraction. During both wet and dry seasons, κCCN, κorg, and particle composition under background conditions exhibited essentially no diel variations. The constant κorg of ∼ 0. 15 is consistent with the largely uniform and high O : C value (∼ 0. 8), indicating that the aerosols under background conditions are dominated by the aged regional aerosol particles consisting of highly oxygenated organic compounds. For air masses strongly influenced by urban pollution and/or local biomass burning, lower values of κorg and organic O : C atomic ratio were observed during night, due to accumulation of freshly emitted particles, dominated by primary organic aerosol (POA) with low hygroscopicity, within a shallow nocturnal boundary layer. The O : C, κorg, and κCCN increased from the early morning hours and peaked around noon, driven by the formation and aging of secondary organic aerosol (SOA) and dilution of POA emissions into a deeper boundary layer, while the development of the boundary layer, which leads to mixing with aged particles from the residual layer aloft, likely also contributed to the increases. The hygroscopicities associated with individual organic factors, derived from PMF (positive matrix factorization) analysis of AMS (aerosol mass spectrometry) spectra, were estimated through multivariable linear regression. For the SOA factors, the variation of the κ value with O : C agrees well with the linear relationship reported from earlier laboratory studies of SOA hygroscopicity. On the other hand, the variation in O : C of ambient aerosol organics is largely driven by the variation in the volume fractions of POA and SOA factors, which have very different O : C values. As POA factors have hygroscopicity values well below the linear relationship between SOA hygroscopicity and O : C, mixtures with different POA and SOA fractions exhibit a steeper slope for the increase in κorg with O : C, as observed during this and earlier field studies. This finding helps better understand and reconcile the differences in the relationships between κorg and O : C observed in laboratory and field studies, therefore providing a basis for improved parameterization in global models, especially in a tropical context.

Thalman, Ryan↗