A Moment-Based Polarimetric Radar Forward Operator for Rain Microphysics

There is growing interest in combining microphysical models and polarimetric radar observations to improve our understanding of storms and precipitation. Mapping model-predicted variables into the radar observational space necessitates a forward operator, which requires assumptions that introduce uncertainties into model-observation comparisons. These include uncertainties arising from the microphysics scheme a priori assumptions of a fixed drop size distribution (DSD) functional form, whereas natural DSDs display far greater variability. To address this concern, this study presents a moment-based polarimetric radar forward operator with no fundamental restrictions on the DSD form by linking radar observables to integrated DSD moments. The forward operator is built upon a dataset of > 200 million realistic DSDs from one-dimensional bin microphysical rain shaft simulations, and surface disdrometer measurements from around the world. This allows for a robust statistical assessment of forward operator uncertainty and quantification of the relationship between polarimetric radar observables and DSD moments. Comparison of "truth" and forward-simulated vertical profiles of the polarimetric radar variables are shown for bin simulations using a variety of moment combinations. Higher-order moments (especially those optimized for use with the polarimetric radar variables: the 6th and 9th) perform better than the lower-order moments (0th and 3rd) typically predicted by many bulk microphysics schemes

Understanding and Visualizing Droplet Distributions in Simulations of Shallow Clouds

Thorough analysis of local droplet-level interactions is crucial to better understand the microphysical processes in clouds and their effect on the global climate. High-accuracy simulations of relevant droplet size distributions from Large Eddy Simulations (LES) of bin microphysics challenge current analysis techniques due to their high dimensionality involving three spatial dimensions, time, and a continuous range of droplet sizes. Utilizing the compact latent representations from Variational Autoencoders (VAEs), we produce novel and intuitive visualizations for the organization of droplet sizes and their evolution over time beyond what is possible with clustering techniques. This greatly improves interpretation and allows us to examine aerosol-cloud interactions by contrasting simulations with different aerosol concentrations. We find that the evolution of the droplet spectrum is similar across aerosol levels but occurs at different paces. This similarity suggests that precipitation initiation processes are alike despite variations in onset times.Comment: 4 pages, 3 figures, accepted at NeurIPS 2023 (Machine Learning and the Physical Sciences Workshop

Use of Polarimetric Radar Measurements to Constrain Simulated Convective Cell Evolution: A Pilot Study with Lagrangian Tracking

To probe the potential value of a radar-driven field campaign to constrain simulation of isolated convection subject to a strong aerosol perturbation, convective cells observed by the operational KHGX weather radar in the vicinity of Houston, Texas, are examined individually and statistically. Cells observed in a single case study of onshore flow conditions during July 2013 are first examined and compared with cells in a regional model simulation. Observed and simulated cells are objectively identified and tracked from observed or calculated positive specific differential phase (K(sub DP)) above the melting level, which is related to the presence of supercooled liquid water. Several observed and simulated cells are subjectively selected for further examination. Below the melting level, we compare sequential cross sections of retrieved and simulated raindrop size distribution parameters. Above the melting level, we examine time series of KDP and radar differential reflectivity (Z(sub DR)) statistics from observations and calculated from simulated supercooled rain properties, alongside simulated vertical wind and supercooled rain mixing ratio statistics. Results indicate that the operational weather radar measurements offer multiple constraints on the properties of simulated convective cells, with substantial value added from derived K(sub DP) and retrieved rain properties. The value of collocated three-dimensional lightning mapping array measurements, which are relatively rare in the continental US, supports the choice of Houston as a suitable location for future field studies to improve the simulation and understanding of convective updraft physics. However, rapid evolution of cells between routine volume scans motivates consideration of adaptive scan strategies or radar imaging technologies to amend operational weather radar capabilities. A 3-year climatology of isolated cell tracks, prepared using a more efficient algorithm, yields additional relevant information. Isolated cells are found within the KHGX domain on roughly 40 % of days year-round, with greatest concentration in the northwest quadrant, but roughly 5-fold more cells occur during June through September. During this enhanced occurrence period, the cells initiate following a strong diurnal cycle that peaks in the early afternoon, typically follow a south-to-north flow, and dissipate within 1 h, consistent with the case study examples. Statistics indicate that 150 isolated cells initiate and dissipate within 70 km of the KHGX radar during the enhanced occurrence period annually, and roughly 10 times as many within 200 km, suitable for multi-instrument Lagrangian observation strategies. In addition to ancillary meteorological and aerosol measurements, robust vertical wind speed retrievals would add substantial value to a radar-driven field campaign

Inverse Modeling to Quantify Microphysical Parameterization Uncertainty

This dissertation document details work completed in the field of cloud microphysical model error and uncertainty quantification. The general motivation is a need to account for error and uncertainty in microphysical parameterization schemes, which remain a large source numerical weather prediction error. The techniques used have been drawn from the fields of nonlinear parameter estimation, data assimilation and uncertainty quantification. First, microphysical parameter uncertainty has been quantified in a manner identical to that of Posselt and Vukicevic (2010), except that a vertically resolved simulated radar reflectivity observation was used as observation constraint on parameter uncertainty, rather than column-integral observations (such as IWP, LWP, etc.). This necessitated work on estimation of the error characteristic of the vertically resolved radar reflectivity observations, in particular, the estimation of a vertically correlated radar error model. Additionally, work has been conducted on estimating uncertainty in microphysics schemes by using perturbations on the hydrometeor time tendency equations as control parameters within a probabilistic Monte Carlo inversion. In order to further facilitate probabilistic inverse modeling studies, and to provide an educational tool for future students, a simplified, Matlab-based framework has been developed for testing advanced Monte Carlo techniques on simple models (e.g. damped harmonic oscillator, Lorenz 1963 model, etc). Finally, an ensemble Kalman transform smoother (EnTKS or ETKS, Posselt and Bishop (2012)) is used to confim that new choices of uncertain model control variables may aid in the use of ensemble Kalman techniques for non-Gaussain model uncertainty analysis

CloudBOSS: Version used for CloudBOSS introductory papers Part 2

This is the version of the CloudBOSS Markov Chain Monte Carlo codes used for part 2 of the introductory papers for CloudBOSS, as well as an offline driver and various scripts used to generate plots and additional data sets.This version has the "ab" calculation bug that affected results in the paper

Dual-Polarization Radar Fingerprints of Precipitation Physics: A Review

This article reviews how precipitation microphysics processes are observed in dual-polarization radar observations. These so-called “fingerprints” of precipitation processes are observed as vertical gradients in radar observables. Fingerprints of rain processes are first reviewed, followed by processes involving snow and ice. Then, emerging research is introduced, which includes more quantitative analysis of these dual-polarization radar fingerprints to obtain microphysics model parameters and microphysical process rates. New results based on a detailed rain shaft bin microphysical model are presented, and we conclude with an outlook of potentially fruitful future research directions

A Bayesian Approach for Statistical–Physical Bulk Parameterization of Rain Microphysics. Part II: Idealized Markov Chain Monte Carlo Experiments

Observationally informed development of a new framework for bulk rain microphysics, the Bayesian Observationally Constrained Statistical-Physical Scheme (BOSS; described in Part I of this study), is demonstrated. This scheme's development is motivated by large uncertainties in cloud and weather simulations associated with approximations and assumptions in existing microphysics schemes. Here, a proof-of-concept study is presented using a Markov chain Monte Carlo sampling algorithm with BOSS to probabilistically estimate microphysical process rates and parameters directly from a set of synthetically generated rain observations. The framework utilized is an idealized steady-state one-dimensional column rainshaft model with specified column-top rain properties and a fixed thermodynamical profile. Different configurations of BOSS-flexibility being a key feature of this approach-are constrained via synthetic observations generated from a traditional three-moment bulk microphysics scheme. The ability to retrieve correct parameter values when the true parameter values are known is illustrated. For cases when there is no set of true parameter values, the accuracy of configurations of BOSS that have different levels of complexity is compared. It is found that addition of the sixth moment as a prognostic variable improves prediction of the third moment (proportional to bulk rain mass) and rain rate. In contrast, increasing process rate formulation complexity by adding more power terms has little benefit-a result that is explained using further-idealized experiments. BOSS rainshaft simulations are shown to well estimate the true process rates from constraint by bulk rain observations, with the additional benefit of rigorously quantified uncertainty of these estimates.6 month embargo; published online: 4 March 2020This item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at [email protected]