Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixing thermodynamics”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Thermodynamic and kinematic drivers of atmospheric boundary layer stability in the central Arctic during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC)

Observations collected during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) provide a detailed description of the impact of thermodynamic and kinematic forcings on atmospheric boundary layer (ABL) stability in the central Arctic. This study reveals that the Arctic ABL is stable and near-neutral with similar frequencies, and strong stability is the most persistent of all stability regimes. MOSAiC radiosonde observations, in conjunction with observations from additional measurement platforms, including a 10 m meteorological tower, ceilometer, microwave radiometer, and radiation station, provide insight into the relationships between atmospheric stability and various atmospheric thermodynamic and kinematic forcings of ABL turbulence and how these relationships differ by season. We found that stronger stability largely occurs in low-wind (i.e., wind speeds are slow), low-radiation (i.e., surface radiative fluxes are minimal) environments; a very shallow mixed ABL forms in low-wind, high-radiation environments; weak stability occurs in high-wind, moderate-radiation environments; and a near-neutral ABL forms in high-wind, high-radiation environments. Surface pressure (a proxy for synoptic staging) partially explains the observed wind speeds for different stability regimes. Cloud frequency and atmospheric moisture contribute to the observed surface radiation budget. Unique to summer, stronger stability may also form when moist air is advected from over the warmer open ocean to over the colder sea ice surface, which decouples the colder near-surface atmosphere from the advected layer, and is identifiable through observations of fog and atmospheric moisture.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Energy Barriers for H 2 Uptake and Release by Defects in Boron Nitride

Density functional theory simulations (DFT) have been used in combination with the climbing image nudged elastic band method (CI-NEB) to determine the energy barriers required to dissociate H 2 at four types of h-BN defects and evaluate the potential for reacted H species to hop to adjacent available sites. The defects were selected based on previous work in which ab initio thermodynamics was used to identify specific defects with optimal thermodynamics for reversible H 2 storage, that is, the defects would dissociate molecular H 2 while not binding H species too strongly so that hydrogen release is favorable. The four defects investigated are the N and B monovacancies, which only offer either B or N defective sites respectively, as well as larger cluster vacancy defects named 3V(1B2N) and hexagonal 6V(3B3N), which present a mixture of B and N defect sites available for H binding. While homogeneous H 2 dissociation by B and N monovacancies have energy barriers larger than 0.4 eV, we found that mixed BN terminated defects have energy barrier lower than 0.4 eV for H 2 dissociation. In particular, H 2 dissociation by the hexagonal 6V(3B3N) defect, exclusively made of Frustrated Lewis Pairs (FLPs), is almost barrierless. The energy barrier for subsequent H hopping to available adjacent B or N sites of the defect were also calculated. It was found that the hopping of H species generally involves high energy barriers ranging from 0.26 eV to 2.57 eV. In conclusion, this work completes previous thermodynamic theoretical investigations and provides a more comprehensive picture of H 2 reaction with defective h-BN for potential H 2 storage and recovery.

Sassi, Michel J. [Pacific Northwest National Labor↗
DOE OSTI2020

Mirostructure Characterization of Friction Consolidated Copper-Nickel using a Machine Learning Approach: Developing Process to Microstructure Associations

Friction consolidation (FC) is a solid phase processing approach where discrete material forms such as powders, chips, nuggets, etc. are densified via shear deformation. The precursors are placed in a billet container and brought in contact with a rotating tool that applying the desirable amount of normal force. Under the combined action of the rotation and normal pressure, the discrete precursor is consolidated through porosity reduction and shear deformation. FC is increasingly being studied as an attractive approach to manufacturing fully dense parts from powder forms owing to its ability to mix, alloy and consolidate difficult-to-process precursors in minimal number of process steps. Material consolidation and deformation in shear consolidation processes have been studied extensively previously for different material combinations previously. However, despite the extensive research in this area, understanding of the mechanistic processes in pore consolidation, deformation-induced mixing and material solubility during FC is still evolving. Material development using solid phase processing approaches such as FC is often performed based on research experience/education, which can be biased. Conventional analysis and simulation tools in this area tend to be successful only when material thermodynamic pathways and microstructural evolution sequences resulting from processing are clearly defined or known. They are not as effective for emerging advanced manufacturing technologies where material evolution pathways are not well established. The ability to predict optimal process parameters based on material chemistry and bulk properties is essential to accelerate materials design and processing, as are an understanding of the relevant structure-processing-property relationships. These structure-processing-property-performance relationships are at the core of materials science research. Microstructure characterization provides the link to these four core areas, often through visualizing material microstructure using imaging techniques. However, linking microstructure image data (i.e., micrographs) to variables of interest (e.g., processing parameters, material chemistry) in a reproducible, generalizable, and quantitative manner is a significant challenge. Typically, quantitatively linking image data to processing history relies on significant domain knowledge and manual or subject matter expert (SME)-heuristic based image analysis. Such an approach to image analysis has the potential to be biased, inefficient, and difficult to replicate.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Physics-Based Model to Represent Membrane-Electrode Assemblies of Solid-Oxide Fuel Cells Based on Gadolinium-Doped Ceria

This paper reports a physics-based model that predicts membrane-electrode assembly (MEA) performance of solid-oxide fuel cells (SOFCs) with Ce 0.9 Gd 0.1 O 2− δ (GDC10) electrolyte membranes. The paper derives self-consistent thermodynamic and transport properties for GDC1o mobile charged defects (oxide vacancies and reduced-ceria small polarons) by fitting published measurements of oxygen non-stoichiometry and conductivity over ranges of temperature and O 2 partial pressures. The button-cell model is applied to evaluate how mixed ionic-electronic conductivity influences the performance of an SOFC MEA with a GDC10 electrolyte sandwiched between a porous, composite Ni-GDC10 anode and a porous, composite cathode of Sm 0.5 Sr 0.5 CoO 3− δ (i.e., SSC) and GDC10. SSC properties are also derived by fitting published conductivity and oxygen non-stoichiometry measurements. Mixed conductivity of GDC10 and competing charge transfer reactions at both electrodes reduce open circuit voltages due to leakage current and buildup of defect concentrations at electrode-electrolyte interfaces. To fit polarization data, the button-cell model includes heterogeneous reaction rates for defect incorporation on the GDC10 surface along with Butler–Volmer expressions derived for competing charge transfer reaction rates from rigorous analyses assuming rate-limiting, elementary charge transfer reactions for each electrode. The calibrated MEA model can support rigorous SOFC modeling with GDC10 electrolytes over the range of conditions within a fully operating cell.

Electrochemistry↗
DOE OSTI2022

High entropy alloys with hexagonal close-packed structure derived from thin film combinatorial approach

A thin film combinatorial synthesis approach was utilized to develop a compositional and structural library for the OsRuWCo alloy system. The survey of this alloy space produced 24 unique compositions and their corresponding crystal structures, which represent a phase map where amorphous and single-phase hexagonal close-packed alloys were identified. Based on the identification of a promising candidate film region, a new high entropy alloy based on OsRuWCo, but with non-equiatomic composition, was synthesized in bulk form. The alloy exhibited a single-phase hexagonal close-packed structure in the as-cast state. Three derivatives from this system were also developed by considering the heat of mixing, atomic size, and binary solubility. These new alloys are based on OsRuWCoIr and OsRuWCoFe, both of which exhibit single-phase hexagonal close-packed as-cast structures, as determined by X-ray diffraction and electron microscopy. Additionally, this large compositional space was utilized to evaluate conventional parameters that describe high entropy alloys (ΔS mix , ΔH mix , Ω, δ, and VEC), and these were found to elucidate the criteria for solid solution formation in this system. Additionally, trends illustrating the evolution from amorphous to crystalline phases are demonstrated and discussed. It is concluded that in this particular case, the dominant factors promoting crystallinity are thermodynamic or electronic, rather than geometric.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Locking water molecules via ternary O–H⋯O intramolecular hydrogen bonds in perhydroxylated closo -dodecaborate

A multitude of applications related to perhydroxylated closo-dodecaborate B 12 (OH) 12 2− in the condensed phase are inseparable from the fundamental mechanisms underlying the high water orientation selectivity based on the base B 12 (OH) 12 2− . Herein, we directly compare the structural evolution of water clusters, ranging from monomer to hexamer, oriented by functional groups in the bases B 12 H 12 2− , B 12 H 11 OH 2− and B 12 (OH) 12 2− using multiple theoretical methods. A significant revelation is made regarding B 12 (OH) 12 2− : each additional water molecule is locked into the intramolecular hydrogen bond B–O–H ternary ring in an embedded form. This new pattern of water cluster growth suggests that B–(H–O)⋯H–O interactions prevail over the competition from water–hydrogen bonds (O⋯H–O), distinguishing it from the behavior observed in B 12 H 12 2− and B 12 H 11 OH 2− bases, in which competition arises from a mixed competing model involving dihydrogen bonds (B–H⋯H–O), conventional hydrogen bonds (B–(H–O)⋯H–O) and water hydrogen bonds (O⋯H–O). Through aqueous solvation and ab initio molecular dynamics analysis, we further demonstrate the largest water clusters in the first hydrated shell with exceptional thermodynamic stability around B 12 (OH) 12 2− . These findings provide a solid scientific foundation for the design of boron cluster chemistry incorporating hydroxyl-group-modified borate salts with potential implications for various applications.

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

The dynamics, mixing, and thermonuclear burn of compressed foams with varied gas fills

Inertial confinement fusion (ICF) implosions involve highly coupled physics and complex hydrodynamics that are challenging to model computationally. Due to the sensitivity of such implosions to small features, detailed simulations require accurate accounting of the geometry and dimensionality of the initial conditions, including capsule defects and engineering features such as fill tubes used to insert gas into the capsule, yet this is computationally prohibitive. It is therefore difficult to evaluate whether discrepancies between the simulation and experiment arise from inadequate fidelity to the capsule geometry and drive conditions, uncertainties in physical data used by simulations, or inadequate physics. We present results from detailed high-resolution three-dimensional simulations of ICF implosions performed as part of the MARBLE campaign on the National Ignition Facility [Albright et al., Phys. Plasmas 29, 022702 (2022)]. These experiments are foam-filled separated-reactant experiments, where deuterons reside in the foam and tritons reside in the capsule gas fill and deuterium–tritium (DT) fusion reactions only occur in the presence of mixing between these materials. Material mixing in these experiments is primarily seeded by shock interaction with the complex geometry of the foam and gas fill, which induces the Richtmyer–Meshkov instability. We compare results for experiments with two different gas fills (ArT and HT), which lead to significant differences in the hydrodynamic and thermodynamic developments of the materials in the implosion. Our simulation results show generally good agreement with experiments and demonstrate a substantial impact of hydrodynamic flows on measured ion temperatures. The results suggest that viscosity, which was not included in our simulations, is the most important unmodeled physics and qualitatively explains the few discrepancies between the simulation and experiment. The results also suggest that the hydrodynamic treatment of shocks is inadequate to predict the heating and yield produced during shock flash, when the shock converges at the center of the implosion. Alternatively, underestimation of the level of radiative preheat from the shock front could explain many of the differences between the experiment and simulation. Nevertheless, simulations are able to reproduce many experimental observables within the level of experimental reproducibility, including most yields, time-resolved X-ray self-emission images, and an increase in burn-weighted ion temperature and neutron down-scattered ratio in the line of sight that includes a jet seeded by the glue spot that joins capsule hemispheres.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Layered Transition Metal Oxides as Ca Intercalation Cathodes: A Systematic First-Principles Evaluation

Finding high-voltage Ca cathode materials is a critical step to unleashing the full potential of high-energy-density Ca-ion batteries. Therefore, first-principles calculations are used to demonstrate that P-type layered calcium transition metal (TM) oxide materials (CaTM 2 O 4 ) with a range of TM substitutions (TM = Ti, V, Cr, Mn, Fe, Co, and Ni) have excellent battery-related properties including thermodynamic stability, average voltage, energy density, synthesizability, ionic mobility, and electronic structure. However, the thermodynamic stability of the charged phase and TM redox activity are shown to be sensitive to TM selection, with CaCo 2 O 4 having the best balance of all considered properties. The utility of combining multiple TMs to expand the chemical search space for TM substitutions is demonstrated by mixing Co and Ni in layered CaTM 2 O 4 .

25 ENERGY STORAGE↗
DOE OSTI2026

Understanding the causes of satellite–model discrepancies in aerosol–cloud interactions using near-LES simulations of marine boundary layer clouds

Aerosol–cloud interactions (ACI) remain the largest source of uncertainty in model estimates of anthropogenic radiative forcing, primarily because of deficiencies in representing aerosol–cloud microphysical processes that lead to inconsistent cloud liquid water path (LWP) responses to aerosol perturbations between observations and models. To investigate this discrepancy, we conducted a series of large-eddy-scale simulations driven by realistic meteorology over the eastern North Atlantic, and evaluated LWP susceptibility, precipitation processes, and boundary layer thermodynamics using satellite and ground-based observations. Simulated LWP responses show a strong dependence on cloud state. Non-precipitating thin clouds exhibit a modest LWP decrease with increasing cloud droplet number concentration (N d ), consistent in sign but weaker in magnitude than satellite estimates, reflecting enhanced turbulent mixing and evaporation. The largest model-observation discrepancy occurs in non-precipitating thick clouds, where simulated LWP susceptibilities are strongly positive (+0.32) while observations indicate large negative values (−0.69). This discrepancy stems from excessive precipitation driven by underestimated entrainment, overly active accretion, and overly broad drop-size distributions in polluted conditions. While our high-resolution setup mitigates the excessive drizzling common in coarser models and captures key regime transitions, these biases persist – highlighting that improved parameterizations of cloud-top processes, precipitation, and aerosol effects are needed beyond simply increasing model resolution. Additionally, misrepresented moisture inversions in reanalysis introduce a moist bias in cloud-top relative humidity, further amplifying positive LWP susceptibility. Our results also suggest that large negative N d –LWP relationships in observations may reflect internal cloud processes rather than true ACI effects.

Aerosol-cloud interaction↗
DOE OSTI2022

Thermodynamic assessment of lithium halide reciprocal salt systems for energy applications

Lithium halides are of interest because of their unique applications as fuel or coolants in molten salt reactors (MSRs) and electrolytes in thermally activated batteries. Because iodine is an important fission product for MSRs, the Molten Salt Thermal Properties Database-Thermochemical (MSTDB-TC) is being expanded to include relevant iodide reciprocal salt systems. In this work, we assess the thermochemical properties of lithium halide pseudobinary systems based on reported experimental phase diagrams, eutectic temperatures, and enthalpies of mixing. Within the framework of the modified quasi-chemical model in the quadruplet approximation (MQMQA), extrapolations to the related pseudoternary and pseudoquaternary systems have been performed. The pseudoquaternary representations of phase equilibria, enthalpy, entropy, and heat capacity are represented in the resultant models. Of interest for generating source terms for accident analysis is the ability to compute the vapor pressures of the dominant iodide species for the systems at any composition. Here, an example of such calculations is provided for the LiF-LiI system at x LiI = 0.01 and x LiI = 0.70, in this case indicating that the LiI vapor species are predominant.

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

Predicting Multicomponent Aqueous Phase Adsorption Equilibria of Organic Acids Using the Generalized Brunauer–Emmett–Teller Isotherm Model

Here, to support process development of adsorptive separation of organic acids, this work presents a thermodynamic modeling methodology to predict multicomponent aqueous phase organic acid adsorption equilibria from single-component adsorption isotherms using the generalized Brunauer–Emmett–Teller isotherm model (gBET). With the organic acid fugacities rigorously accounting for the aqueous phase organic acid solution chemistry and solution nonideality, gBET precisely represents the single-component adsorption isotherms and accurately predicts the binary and ternary mixed-acid adsorption equilibria for the ranges of initial pH (∼3–7), acid concentration (100–400 mmol/L), and temperature (298.15–328.15 K) with less than 10% average absolute relative deviation. In addition, gBET with pH-independent parameters provides insights into the underlying adsorption phenomena, including the adsorbate loadings and compositions in the monolayer and subsequent layers under varying initial pH, temperature, concentration, and composition. The gBET model predictions outperform the predictions from the overloading model and the Ideal Dilute Solution Theory.

09 BIOMASS FUELS↗
DOE OSTI2024

Mixed-phase Convective Clouds in the Polar Marine Boundary Layer

This proposal, a collaboration between two universities, firstly aims to describe the high-latitude CAO cloud regime using the array of instruments at the NSA (Barrow, Alaska) and at two mobile ARM sites (MARCUS 2017-18, and COMBLE 2020). These ARM data, and the multi-sensor variables we propose to derive, allow detailed descriptions of mesoscale organization of precipitation, thermodynamic profiles, vertical velocity, cloud depth, and cloud and precipitation properties. These observations are key for our second goal, which is to use output from validated cloud- and eddy-resolving simulations to gain insights into the linkages between microphysical processes, the MBL vertical structure, and the MBL dynamics that control cloud macrostructure.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Pulsed Light Synthesis of High Entropy Nanocatalysts with Enhanced Catalytic Activity and Prolonged Stability for Oxygen Evolution Reaction

The ability to synthesize compositionally complex nanostructures rapidly is a key to high-throughput functional materials discovery. In addition to being time-consuming, a majority of conventional materials synthesis processes closely follow thermodynamics equilibria, which limit the discovery of new classes of metastable phases such as high entropy oxides (HEO). Herein, a photonic flash synthesis of HEO nanoparticles at timescales of milliseconds is demonstrated. By leveraging the abrupt heating and cooling cycles induced by a high-power-density xenon pulsed light, mixed transition metal salt precursors undergo rapid chemical transformations. Hence, nanoparticles form within milliseconds with a strong affinity to bind to the carbon substrate. Oxygen evolution reaction (OER) activity measurements of the synthesized nanoparticles demonstrate two orders of magnitude prolonged stability at high current densities, without noticeable decay in performance, compared to commercial IrO 2 catalyst. This superior catalytic activity originates from the synergistic effect of different alloying elements mixed at a high entropic state. It is found that Cr addition influences surface activity the most by promoting higher oxidation states, favoring optimal interaction with OER intermediates. The proposed high-throughput method opens new pathways toward developing next-generation functional materials for various electronics, sensing, and environmental applications, in addition to renewable energy conversion.

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

Breaking the Brønsted–Evans–Polanyi Relation with Dual-Metal Sites

Linear scaling relationships impose inherent limitations on catalyst activity; the Brønsted−Evans−Polanyi (BEP) relation, which correlates activation and reaction energies, is a prominent example. Here we report a dual-metal site catalyst (DMSC) on ceria that breaks the BEP relation for C−C coupling of methyl intermediates an elementary step in methane coupling to form ethane. The DMSC structure on CeO 2 (111) was discovered by density-functional theory (DFT) structural exploration and confirmed to be stable via ab initio thermodynamics and ab initio molecular dynamics. Homonuclear and heteronuclear DMSCs of Ni, Pd, Pt, Fe, Ru, Os, Co, Rh, and Ir (45 pairs in total) were examined for methyl affinity and methyl−methyl coupling activation energy. We found that many heteronuclear DMSCs break the BEP linear scaling due to a mixed low-affinity/high-affinity coadsorption of the two methyl groups, decoupling the step responsible for the activation energy (E a ) at the low-affinity site from the overall reaction energy (ΔE) determined by both sites. This mechanism of breaking the BEP relationship via the DMSCs offers a catalyst design principle for C−C coupling reactions.

Chen, Yiming [Vanderbilt Univ., Nashville, TN (Uni↗
DOE OSTI2026

High entropy oxides prediction and discovery by the Mixed Enthalpy-Entropy Descriptor

The vast, high-dimensional composition space of high-entropy oxides (HEOs) offers exceptional opportunities for functional materials discovery, yet it also poses a fundamental challenge: the rational and efficient prediction of stable, synthesizable compositions and the corresponding structure–property relationships. Despite growing interest, the field still lacks broadly applicable, physically grounded descriptors capable of navigating various large chemical spaces. Here, we introduce a Mixed Enthalpy–Entropy Descriptor (MEED) that enables rapid, first-principles–based prediction of HEOs synthesizability across diverse chemistries. Using MEED, we perform high-throughput screening of two distinct HEO families: rocksalt oxides and perovskite oxides. The predicted top candidates in each family were experimentally validated. MEED reveals unifying thermodynamic and structural principles governing stability across both chemical compositions and polymorphs, providing mechanistic insight into the formation of high-entropy phases. This work significantly broadens the accessible chemical design space for HEOs and establishes a data-efficient framework for accelerating the discovery of next-generation functional materials.

Yu, Liping [University of Central Florida]↗
DOE OSTI2020

Inorganic Halide Double Perovskites with Optoelectronic Properties Modulated by Sublattice Mixing

All-inorganic halide double perovskites have emerged as a promising class of materials that are potentially more stable and less toxic than lead-containing hybrid organic-inorganic perovskite optoelectronic materials. In this work, 311 cesium chloride double perovskites (Cs 2 BB'Cl 6 ) were selected from a set of 903 compounds as likely being stable on the basis of a statistically learned tolerance factor (t) for perovskite stability. First-principles calculations on these 311 double perovskites were then performed to assess their stability and identify candidates with band gaps appropriate for optoelectronic applications. We predict that 261 of the 311 Cs 2 BB'Cl 6 compounds are likely synthesizable on the basis of a thermodynamic analysis of their decomposition to competing compounds (decomposition enthalpy <0.05 eV/atom). Of these 261 likely synthesizable compounds, 47 contain no toxic elements and have direct or nearly direct (within 100 meV) band gaps between 1 and 3 eV, as computed with hybrid density functional theory (HSE06). Within this set, we identify the triple-alkali perovskites Cs 2 [Alk] + [TM] 3+ Cl 6 , where Alk is a group 1 alkali cation and TM is a transition-metal cation, as a class of Cs 2 BB'Cl 6 double perovskites with remarkable optical properties, including large and tunable exciton binding energies as computed by the GW-Bethe-Salpeter equation (GW-BSE) method. We attribute the unusual electronic structure of these compounds to the mixing of the Alk-Cl and TM-Cl sublattices, leading to materials with small band gaps, large exciton binding energies, and absorption spectra that are strongly influenced by the identity of the transition metal. The role of the double-perovskite structure in enabling these unique properties is probed through an analysis of the electronic structures and chemical bonding of these compounds in comparison with other transition-metal and alkali transition-metal halides.

14 SOLAR ENERGY↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2023

Asymmetries in cloud microphysical properties ascribed to sea ice leads via water vapour transport in the central Arctic

To investigate the influence of sea ice openings like leads on wintertime Arctic clouds, the air mass transport is exploited as a heat and humidity feeding mechanism which can modify Arctic cloud properties. Cloud microphysical properties in the central Arctic are analysed as a function of sea ice conditions during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in 2019–2020. The Cloudnet classification algorithm is used to characterize the clouds based on remote sensing observations and the atmospheric thermodynamic state from the observatory on board the research vessel (RV) Polarstern. To link the sea ice conditions around the observational site with the cloud observations, the water vapour transport (WVT) being conveyed towards RV Polarstern has been utilized as a mechanism to associate upwind sea ice conditions with the measured cloud properties. This novel methodology is used to classify the observed clouds as coupled or decoupled to the WVT based on the location of the maximum vertical gradient of WVT height relative to the cloud-driven mixing layer. Only a conical sub-sector of sea ice concentration (SIC) and the lead fraction (LF) centred on the RV Polarstern location and extending up to 50 km in radius and with an azimuth angle governed by the time-dependent wind direction measured at the maximum WVT is related to the observed clouds. We found significant asymmetries for cases when the clouds are coupled or decoupled to the WVT and selected by LF regimes. Liquid water path of low-level clouds is found to increase as a function of LF, while the ice water path does so only for deep precipitating systems. Clouds coupled to WVT are found to generally have a lower cloud base and larger thickness than decoupled clouds. Thermodynamically, for coupled cases the cloud-top temperature is warmer and accompanied by a temperature inversion at the cloud top, whereas the decoupled cases are found to be closely compliant with the moist adiabatic temperature lapse rate. The ice water fraction within the cloud layer has been found to present a noticeable asymmetry when comparing coupled versus decoupled cases. This novel approach of coupling sea ice to cloud properties via the WVT mechanism unfolds a new tool to study Arctic surface–atmosphere processes. With this formulation, long-term observations can be analysed to enforce the statistical significance of the asymmetries. Furthermore, our results serve as an opportunity to better understand the dynamic linkage between clouds and sea ice and to evaluate its representation in numerical climate models for the Arctic system.

54 ENVIRONMENTAL SCIENCES↗