Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reflectivity”

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 1,063 records · Page 59

Thermodynamics shapes the in vivo enzyme burden of glycolytic pathways

Thermodynamically constrained reactions and pathways are hypothesized to impose greater protein demands on cells, requiring higher enzyme amounts to sustain a given flux compared to those with stronger thermodynamics. To test this, we quantified the absolute concentrations of glycolytic enzymes in three bacterial species— Zymomonas mobilis, Escherichia coli , and Clostridium thermocellum —which employ distinct glycolytic pathways with varying thermodynamic driving forces. By integrating enzyme concentration data with corresponding in vivo metabolic fluxes and ΔG measurements, we found that the highly favorable Entner-Doudoroff pathway in Z. mobilis requires only one-fourth the amount of enzymatic protein to sustain the same flux as the thermodynamically constrained pyrophosphate-dependent glycolytic pathway in C. thermocellum , with the Embden-Meyerhof-Parnas pathway in E. coli exhibiting intermediate thermodynamic favorability and enzyme demand. Across all three pathways, early reactions with stronger thermodynamic driving forces generally required lower enzyme investment than later, less favorable steps. Additionally, reflecting differences in glycolytic strategies, the highly reversible ethanol fermentation pathway in C. thermocellum requires 10-fold more protein to maintain the same flux as the irreversible, forward-driven ethanol fermentation pathway in Z. mobilis . Thus, protein investment across glycolytic pathways reflects differences in their thermodynamic favorability.

Clostridium thermocellum↗

Atypical Divergence of SARS-CoV-2 Orf8 from Orf7a within the Coronavirus Lineage Suggests Potential Stealthy Viral Strategies in Immune Evasion

Orf8, one of the most puzzling genes in the SARS lineage of coronaviruses, marks a unique and striking difference in genome organization between SARS-CoV-2 and SARS-CoV-1. Here, using sequence comparisons, we unequivocally reveal the distant sequence similarities between SARS-CoV-2 Orf8 with its SARS-CoV-1 counterparts and the X4-like genes of coronaviruses, including its highly divergent “paralog” gene Orf7a, whose product is a potential immune antagonist of known structure. Supervised sequence space walks unravel identity levels that drop below 10% and yet exhibit subtle conservation patterns in this novel superfamily, characterized by an immunoglobulin-like beta sandwich topology. We document the high accuracy of the sequence space walk process in detail and characterize the subgroups of the superfamily in sequence space by systematic annotation of gene and taxon groups. While SARS-CoV-1 Orf7a and Orf8 genes are most similar to bat virus sequences, their SARS-CoV-2 counterparts are closer to pangolin virus homologs, reflecting the fine structure of conservation patterns within the SARS-CoV-2 genomes. The divergence between Orf7a and Orf8 is exceptionally idiosyncratic, since Orf7a is more constrained, whereas Orf8 is subject to rampant change, a peculiar feature that may be related to hitherto-unknown viral infection strategies. Despite their common origin, the Orf7a and Orf8 protein families exhibit different modes of evolutionary trajectories within the coronavirus lineage, which might be partly attributable to their complex interactions with the mammalian host cell, reflected by a multitude of functional associations of Orf8 in SARS-CoV-2 compared to a very small number of interactions discovered for Orf7a.

60 APPLIED LIFE SCIENCES↗

Large $$N_c$$ QCD phase diagram at $$\mu _B=0$$

Abstract Lattice studies suggest that at zero baryon chemical potential and increasing temperature there are three characteristic regimes in QCD that are connected by smooth analytical crossovers: a hadron gas regime at$$T < T_{ch}\sim 155$$ T < T ch ∼ 155 MeV, an intermediate regime, called stringy fluid, at$$T_{ch}< T < \sim 3 T_{ch}$$ T ch < T < ∼ 3 T ch , and a quark-gluon plasma regime at higher temperatures. These regimes have been interpreted to reflect different approximate symmetries and effective degrees of freedom. In the hadron gas the effective degrees of freedom are hadrons and the approximate chiral symmetry of QCD is spontaneously broken. The intermediate regime has been interpreted as lacking spontaneous chiral symmetry breaking along with the emergence of new approximate symmetry, chiral spin symmetry, that is not a symmetry of the Dirac Lagrangian, but is a symmetry of the confining part of the QCD Lagrangian. While the high temperature regime is the usual quark-gluon plasma which is often considered to reflect “deconfinement” in some way. This paper explores the behavior of these regimes of QCD as the number of colors in the theory,$$N_c$$ N c , gets large. In the large$$N_c$$ N c limit the theory is center-symmetric and notions of confinement and deconfinement are unambiguous. The energy density is$$\mathcal{O}(N_c^0)$$ O ( N c 0 ) in the meson gas,$${{\mathcal {O}}}(N_c^1)$$ O ( N c 1 ) in the intermediate regime and$${{\mathcal {O}}}(N_c^2)$$ O ( N c 2 ) in the quark-gluon plasma regime. In the large$$N_c$$ N c limit these regimes may become distinct phases separated by first order phase transitions. The intermediate phase has the peculiar feature that glueballs should exist and have properties that are unchanged from what is seen in the vacuum (up to$$1/N_c $$ 1 / N c corrections), while the ordinary dilute gas of mesons with broken chiral symmetry disappears and approximate chiral spin symmetry should emerge.

Physics↗

Shape phase changes with N in 72-84 Kr isotopes

The Kr isotopes lying in between the lighter isotopes of (Zn, Ge, and Se) and the heavier isotopes of (Sr and Zr) in the A=70-80 region exhibit very interesting spectral features. The spectra of 72-84 Kr isotopes display a striking contrast from those of Zn, Ge and Se isotopes. The role of spherical and oblate and prolate deformed subshell gaps at specific Z and N and the resulting re-inforcement are strikingly evident in these contrasting features, with variation in Z or N. The evolution of the spectral features in Kr isotopes with N as reflected in the quadrupole deformation, K-band structures, E0 decay, B(E2) values, β-softness of the nuclear core, and odd–even staggering in K π = 2 + γ-bands is studied to explore the role of the under lying nuclear interactions. Further, the correlations with odd–A isotopes are explored. The shape co-existence displayed in some Kr isotopes is studied. The large deformation observed in the ground bands of 74,76 Kr, as exhibited in the B(E2) values is especially interesting.

74 ATOMIC AND MOLECULAR PHYSICS↗

Perspective on Tsallis statistics for nuclear and particle physics

This is a concise introduction to the topic of nonextensive Tsallis statistics meant especially for those interested in its relation to high-energy proton–proton, proton–nucleus and nucleus–nucleus collisions. The three types of Tsallis statistics are reviewed. Only one of them is consistent with the fundamental hypothesis of equilibrium statistical mechanics. The single-particle distributions associated with it, namely Boltzmann, Fermi–Dirac and Bose–Einstein, are derived. These are not equilibrium solutions to the conventional Boltzmann transport equation which must be modified in a rather nonintuitive manner for them to be so. Nevertheless, the Boltzmann limit of the Tsallis distribution is extremely efficient in representing a wide variety of single-particle distributions in high-energy proton–proton, proton–nucleus and nucleus–nucleus collisions with only three parameters, one of them being the so-called nonextensitivity parameter [Formula: see text]. This distribution interpolates between an exponential at low transverse energy, reflecting thermal equilibrium, to a power law at high transverse energy, reflecting the asymptotic freedom of Quantum Chromodynamics (QCD). It should not be viewed as a fundamental new parameter representing nonextensive behavior in these collisions.

Physics↗

PVwindow (Photovoltaic Window Simulator) [SWR-22-01]

PVwindow simulates photovoltaic window properties by allowing the user build stacks of thin film materials that compose the window layers. The only needed inputs are the complex refractive indices and thicknesses of the materials. The software will simulate the optical transmission, front reflection, and back reflection of the stack using the transfer matrix method. Properties specific to PV window optics such as visible transmittance and color are also calculated. The optics code extracts the number of absorbed photons in the PV absorber layer and solves the single diode equation to yield theoretical photovoltaic metrics such as power conversion efficiency, fill factor, open circuit voltage, and short circuit current of the PV window when illuminated at an arbitrary angle or from either side of the stack. The software tool ultimately allows researchers to design photovoltaic windows and simulate their properties. The software also exports files that can be read into LBL Window software or into the EnergyPlus and Open Studio building energy modeling software to determine larger scale energy impacts of PV window integration.

Wheeler, Lance↗

BCS (CSP Beam Characterization System) [SWR-23-13]

The Beam Characterization System (BCS) is an image-processing tool for estimating heliostat tracking errors. A large mirror, or heliostat, reflects the sun onto a distant target. A carefully calibrated camera then takes an image of this target to capture the shape of the reflected beam. The software kit hosts a series of scripts to identify key features and run analysis on the captured images.

Kesseli, Devon↗

LATTE: Los Alamos TravelTime package based on Eikonal equation

This Fortran code focuses on traveltime computation and tomography based on eikonal equation. Specifically, the package provides three major functionalities: (1) forward modeling of traveltime from single-point or ensemble source based on factorized eikonal equation, (2) adjoint-state first-arrival traveltime tomography based on picked first arrival traveltime using steepest descent, conjugate gradient, or limited-memory BFGS inversion scheme, and (3) adjoint-state joint transmission-reflection tomography based on picked first-arrival and reflection traveltimes. The package applies to forward modeling and tomography based on traveltime in 2D and 3D isotropic regular-grid models. We name this package LATTE – Los Alamos TravelTime package based on Eikonal equation. * The code is for accompanying a journal paper under preparation. The paper will be submitted via LA-UR separately later.

Gao, Kai↗

OpenCSP v.1.0

SAND2024-08987O OpenCSP analyzes and improves the performance of optical mirrors that reflect sunlight and concentrate it onto a receiver. The algorithms in this library analyze the flow of light from an original source—such as the sun—reflecting off one or more mirrors, and onto a receiver, optical screen, or camera. OpenCSP also analyzes optical systems using other techniques, such as shape or slope comparison. The software calculates aggregate properties of the optical system, such as its efficiency, manufacturing error, pointing control error, expected performance, and other useful parameters. OpenCSP also includes high-level tools, such as Open SOFAST, a program that obtains high-resolution measurement of CSP mirrors, full heliostats, and other CSP collectors. OpenCSP supports producing meaningful output data to help CSP designers and field operators improve the performance of their systems. OpenCSP also includes classes and functions that support these calculations, including input utilities, analysis output libraries, geometric constraint analysis for drone flight planning, and more. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Brost, Randolph↗

Fireworks on Weather Radar and Camera

High-sensitivity weather radars easily detect nonmeteorological phenomena characterized by weak radar returns. Fireworks are the example presented here. To understand radar observations, an experiment was conducted in which the National Severe Storms Laboratory (NSSL)’s research (3-cm wavelength) dual-polarization radar and a video camera were located at 1 km from fireworks in Norman, Oklahoma. The fireworks from the 4 July 2017 celebration were recorded by both instruments. The experiment is described. Few bursts recorded by the camera are analyzed to obtain the height of the explosion, its maximum diameter, number of stars, and the duration of the visible image. Radar volume scans are examined to characterize the height of the observation, the maximum reflectivity, and its distribution with height. The fireworks location is close to the Terminal Doppler Weather Radar (TDWR) that operates in single polarization at a 5-cm wavelength and monitors hazardous weather over the Oklahoma City airport. A third radar with data from the event is the Weather Surveillance Radar-1988 Doppler (WSR-88D) located in Norman. It has a wavelength of 10 cm and supports technical developments at the Radar Operation Center. Reflectivity factors measured by the three radars are compared to infer the size of dominant scatterers. The polarimetric characteristics of fireworks returns are analyzed. Although these differ from those of precipitation, they are indistinguishable from insect returns. Radar observation of larger fireworks in Fort Worth, Texas, with a WSR-88D is included and compared with the observations of the smaller fireworks in Norman. We expect the detectability of explosions would be similar as of fireworks. Pinpointing locations would be useful to first responders, or air quality forecasters. We note a benefit of fireworks recognition in weather radar data is that it can prevent contamination of precipitation accumulations.

54 ENVIRONMENTAL SCIENCES↗

Effective Visualization of Radar Data for Users Impacted by Color Vision Deficiency

Color vision deficiency (CVD) is a decreased ability to discern between particular colors. Eight percent of genetic males and half a percent of genetic females have some form of CVD, with many in the radar community falling into this group. When presenting data on a two-dimensional plane, it is common to use colors to represent values via a colormap. Colormap choice in the radar community is influenced by the ability to highlight scientifically interesting features in data, institutional choices, and domain dominance of legacy colormaps. The problem with these current colormaps is that many do not project well for those with CVD (i.e., green next to red). In working with the CVD community to address this problem, multiple colormaps for moments such as equivalent reflectivity factor and Doppler velocity were created for users with forms of CVD such as deuteranomaly, protanomaly, protanopia, and deuteranopia using Python tools such as colorspacious and viscm. We show how these colormaps can improve interpretability for four cases: a mesoscale convective system, a pyrocumulonimbus storm, a wintertime midlatitude cyclone, and widespread storms with a large bird migration. These new radar equivalent reflectivity factor, Doppler velocity, and polarization colormaps are designed to highlight rain, frozen precipitation, nonmeteorological targets, and velocity-based items, are perceptually uniform, and are visually friendly for those with CVD.

54 ENVIRONMENTAL SCIENCES↗

Synergetic Use of the WSR-88D Radars, GOES-R Satellites, and Lightning Networks to Study Microphysical Characteristics of Hurricanes

Abstract This study analyzes the microphysics and precipitation pattern of Hurricanes Harvey (2017) and Florence (2018) in both the eyewall and outer rainband regions. From the retrievals by a satellite red–green–blue scheme, the outer rainbands show a strong convective structure while the inner eyewall has less convective vigor (i.e., weaker upper-level reflectivities and electrification), which may be related to stronger vertical wind shear that hinders fast vertical motions. The WSR-88D column-vertical profiles further confirm that the outer rainband clouds have strong vertical motion and large ice-phase hydrometeor formation aloft, which correlates well with 3D Lightning Mapping Array source counts in height and time. From the results from this study, it is determined that the inner eyewall region is dominated by warm rain, whereas the external rainband region contains intense mixed-phase precipitation. External rainbands are defined here as those that reside outside of the main hurricane circulation, associated with surface tropical storm wind speeds. The synergy of satellite and radar dual-polarization parameters is instrumental in distinguishing between the key microphysical features of intense convective rainbands and the warm-rain-dominated eyewall regions within the hurricanes. Substantial amounts of ice aloft and intense updrafts in the external rainbands are indicative of heavy surface precipitation, which can have important implications for severe weather warnings and quantitative precipitation forecasts. The novel part of this study is to combine ground-based radar measurement with satellite observations to study hurricane microphysical structure from surface to cloud top so as to fill in the gaps between the two observational techniques.

Meteorology & Atmospheric Sciences↗

WRF Modelling of Deep Convection and Hail for Wind Power Applications

Deep convection and the related occurrence of hail, intense precipitation and wind gusts represent a hazard to a range of energy infrastructure including wind turbine blades. Wind turbine blade leading edge erosion (LEE) is caused by the impact of falling hydrometeors onto rotating wind turbine blades. It is a major source of wind turbine maintenance costs and energy losses from wind farms. In the United States Southern Great Plains (SGP), where there is widespread wind energy development, deep convection and hail events are common, increasing the potential for precipitation-driven LEE. A 25-day Weather Research and Forecasting (WRF) model simulation conducted at convection-permitting resolution and using a detailed microphysics scheme is carried out for the SGP to evaluate the effectiveness in modelling the wind and precipitation conditions relevant to LEE potential. WRF output for these properties is evaluated using radar observations of precipitation (including hail) and reflectivity, in situ wind speed measurements and wind power generation. This research demonstrates some skill for the primary drivers of LEE. Wind speeds, rainfall rates and precipitation totals show good agreement with observations. Here, the occurrence of precipitation during power producing wind speeds is also shown to exhibit fidelity. Hail events frequently occur during periods when wind turbines are rotating and are especially important to LEE in the SGP. The presence of hail is modelled with a mean proportion correct of 0.77 and an odds ratio 4.55. Further research is needed to demonstrate sufficient model performance to be actionable for the wind energy industry and there is evidence for positive model bias in cloud reflectivity.

17 WIND ENERGY↗

Radar Retrieval Evaluation and Investigation of Dendritic Growth Layer Polarimetric Signatures in a Winter Storm

Abstract This study evaluates ice particle size distribution and aspect ratio φ Multi-Radar Multi-Sensor (MRMS) dual-polarization radar retrievals through a direct comparison with two legs of observational aircraft data obtained during a winter storm case from the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) campaign. In situ cloud probes, satellite, and MRMS observations illustrate that the often-observed K dp and Z DR enhancement regions in the dendritic growth layer can either indicate a local number concentration increase of dry ice particles or the presence of ice particles mixed with a significant number of supercooled liquid droplets. Relative to in situ measurements, MRMS retrievals on average underestimated mean volume diameters by 50% and overestimated number concentrations by over 100%. IWC retrievals using Z DR and K dp within the dendritic growth layer were minimally biased relative to in situ calculations where retrievals yielded −2% median relative error for the entire aircraft leg. Incorporating φ retrievals decreased both the magnitude and spread of polarimetric retrievals below the dendritic growth layer. While φ radar retrievals suggest that observed dendritic growth layer particles were nonspherical (0.1 ≤ φ ≤ 0.2), in situ projected aspect ratios, idealized numerical simulations, and habit classifications from cloud probe images suggest that the population mean φ was generally much higher. Coordinated aircraft radar reflectivity with in situ observations suggests that the MRMS systematically underestimated reflectivity and could not resolve local peaks in mean volume diameter sizes. These results highlight the need to consider particle assumptions and radar limitations when performing retrievals. significance statement Developing snow is often detectable using weather radars. Meteorologists combine these radar measurements with mathematical equations to study how snow forms in order to determine how much snow will fall. This study evaluates current methods for estimating the total number and mass, sizes, and shapes of snowflakes from radar using images of individual snowflakes taken during two aircraft legs. Radar estimates of snowflake properties were most consistent with aircraft data inside regions with prominent radar signatures. However, radar estimates of snowflake shapes were not consistent with observed shapes estimated from the snowflake images. Although additional research is needed, these results bolster understanding of snow-growth physics and uncertainties between radar measurements and snow production that can improve future snowfall forecasting.

Meteorology & Atmospheric Sciences↗

Frozen Hydrometeor Terminal Fall Velocity Dependence on Particle Habit and Riming as Observed by Vertically Pointing Radars

Vertically pointing Ka-band radar measurements are used to derive fall velocity–reflectivity factor ($V$ t = $aZ$$^{b}_{e}$) relations for frozen hydrometeor populations of different habits during snowfall events observed at Oliktok Point, Alaska, and at the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC). Case study events range from snowfall with highly rimed particles observed during periods with large amounts of supercooled liquid water path (LWP > 320 g m –2 ) to unrimed snowflakes including instances when pristine planar crystals were the dominant frozen hydrometeor habit. The prefactor a and the exponent b in the observed $V$ t –$Z$ e relations scaled to the sea level vary in the approximate ranges 0.5–1.4 and 0.03–0.13, respectively (reflectivities are in mm 6 m –3 and velocities are in m s –1 ). The coefficient a values are the smallest for planar crystals (a ~ 0.5) and the largest (a > 1.2) for particles under severe riming conditions with high LWP. There is no clear distinction between b values for high and low LWP conditions. The range of the observed $V$ t –$Z$ e relation coefficients is in general agreement with results of modeling using fall velocity–size (υ t = αD β ) relations for individual particles found in literature for hydrometeors of different habits, though there is significant variability in α and β coefficients from different studies even for a same particle habit. Correspondences among coefficients in the $V$ t –$Z$ e relations for particle populations and in the individual particle υ t –$D$ relations are analyzed. Furthermore, these correspondences and the observed $V$ t –$Z$ e relations can be used for evaluating different frozen hydrometeor fall velocity parameterizations in models.

54 ENVIRONMENTAL SCIENCES↗

Improved Diagnosis of Precipitation Type with LightGBM Machine Learning

Abstract Existing precipitation-type algorithms have difficulty discerning the occurrence of freezing rain and ice pellets. These inherent biases are not only problematic in operational forecasting but also complicate the development of model-based precipitation-type climatologies. To address these issues, this paper introduces a novel light gradient-boosting machine (LightGBM)-based machine learning precipitation-type algorithm that utilizes reanalysis and surface observations. By comparing it with the Bourgouin precipitation-type algorithm as a baseline, we demonstrate that our algorithm improves the critical success index (CSI) for all examined precipitation types. Moreover, when compared with the precipitation-type diagnosis in reanalysis, our algorithm exhibits increased F1 scores for snow, freezing rain, and ice pellets. Subsequently, we utilize the algorithm to compute a freezing-rain climatology over the eastern United States. The resulting climatology pattern aligns well with observations; however, a significant mean bias is observed. We interpret this bias to be influenced by both the algorithm itself and assumptions regarding precipitation processes, which include biases associated with freezing drizzle, precipitation occurrence, and regional synoptic weather patterns. To mitigate the overall bias, we propose increasing the precipitation cutoff from 0.04 to 0.25 mm h −1 , as it better reflects the precision of precipitation observations. This adjustment yields a substantial reduction in the overall bias. Finally, given the strong performance of LightGBM in predicting mixed precipitation episodes, we anticipate that the algorithm can be effectively utilized in operational settings and for diagnosing precipitation types in climate model outputs. Significance Statement Freezing rain can have significant impacts on transportation and infrastructure, making accurate prediction of precipitation types crucial. In this study, we use a machine learning method known as LightGBM to predict precipitation types. We show that the new algorithm performs better than the existing methods for all precipitation types examined. Additionally, we compute a freezing-rain climatology over the eastern United States. Although the resulting climatology pattern corresponds well to observations, the algorithm overpredicts freezing-rain occurrence. We argue that this bias can be substantially reduced by increasing the precipitation cutoff from 0.04 to 0.25 mm h −1 . Overall, this work highlights the potential of the LightGBM algorithm for both weather forecasting and diagnosing precipitation types in climate models.

Meteorology & Atmospheric Sciences↗

A New Method for Ice-Ice Aggregation in the Adaptive Habit Model

In this paper, a novel methodology for modeling ice-ice aggregation is presented. This methodology combines a modified hydrodynamic collection algorithm with bulk aggregate characteristic information from an offline simulator that collects ice particles, namely, the Ice Particle and Aggregate Simulator, and has been implemented into the Adaptive Habit Microphysics scheme in the Weather Research and Forecasting Model. Aggregates, or snow, are formed via collection of cloud ice particles, where initial ice characteristics and the resulting geometry determine aggregate characteristics. Upon implementation, idealized squall-line simulations are performed to examine the new methodology in comparison with commonly used bulk microphysics schemes. It is found that the adaptive habit aggregation parameterization develops snow and reduces ice mass and number concentrations compared to other schemes. The development of aggregates through the new methodology cascades into other interesting effects, including enhancements in ice and snow growth, as well as homogeneous freezing. Further microphysical analyses reveal varying sensitivities, where snow processes are most sensitive to the new parameterization, followed by ice, then cloud, rain, and graupel processes. Further, the new scheme results in enhancements in surface precipitation due to the persistence of snow at lower altitudes. This persistence is a result of shape-dependent melting and sublimation, increasing the residence time. Moreover, these low-level enhancements are reflected in increases in radar reflectivity at the surface and its spatial distribution. Finally, the ability to predict snow shape and density allows for the simulation of polarimetric radar quantities, resulting in signature enhancements compared to schemes that do not consider spatial and temporal variations in snow shape and density.

54 ENVIRONMENTAL SCIENCES↗

How Snow Aggregate Ellipsoid Shape and Orientation Variability Affects Fall Speed and Self-Aggregation Rates

Snow aggregate shapes and orientations have long been known to exhibit substantial variability. Despite this observed variability, most weather and climate prediction models use fixed power-law functions that deterministically map particle size to mass and fall speed. As such, integrated quantities like precipitation and self-aggregation rates currently ignore nonlinear effects resulting from variation in shape and orientation for aggregates of the same size. This study therefore develops an analytic framework that couples an empirically based bivariate distribution of ellipsoid shapes to classical hydrodynamic theory so as to capture an appropriate dispersion of masses, projected areas, and fall speeds for an assumed size distribution. For a fixed aggregate size, shape variations produce approximately ±0.13 m s -1 standard deviation of fall speed which increases the mass flux fall speed dispersion by more than 100% over traditional microphysics models. This increased fall speed dispersion results predominantly from shape-induced mass dispersion whereas orientation and drag dispersion play a lesser role. Shape variations can increase mass- and reflectivity-weighted fall speeds by up to 60% of traditional models whereas self-aggregation rates can increase by a factor of 100 for very small slope parameters. This implies that aggregate shape variations effectively forestall the theorized onset of fall speed distribution narrowing and subsequent quenching of the aggregation process. As a result, it is likely that secondary ice formation is necessary to prevent an ever decreasing slope parameter. The mathematical theory presented in this study is used to develop simple correction factors for snow forecast and climate models.

54 ENVIRONMENTAL SCIENCES↗