Engineering Papers⌕ Search

Engineering topics

Morrison, Hugh

Publications and source records attributed to Morrison, Hugh.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Bayesian Approach for Statistical–Physical Bulk Parameterization of Rain Microphysics. Part I: Scheme Description

A novel framework is proposed for the bulk parameterization of rain microphysics: the Bayesian Observationally Constrained Statistical–Physical Scheme (BOSS). It is designed to facilitate direct constraint by observations using Bayesian inference. BOSS combines existing process-level microphysical knowledge with flexible process rate formulations and parameters constrained by observations within a Bayesian framework. Furthermore, using a raindrop size distribution (DSD) normalization method that relates DSD moments to one another via generalized power series, generalized multivariate power expressions are derived for the microphysical process rates as functions of a set of prognostic DSD moments. The scheme is flexible and can utilize any number and combination of prognostic moments and any number of terms in the process rate formulations. This means that both uncertainty in parameter values and structural uncertainty associated with the process rate formulations can be investigated systematically, which is not possible using traditional schemes. In this paper, BOSS is compared to two- and three-moment versions of a traditional bulk rain microphysics scheme (denoted as MORR). It is demonstrated that some process formulations in MORR are analytically equivalent to the generalized power expressions in BOSS using one or two terms, while others are not. BOSS is able to replicate the behavior of MORR in idealized one-dimensional rainshaft tests, but with a much more flexible and systematic design. Part II of this study describes the application of BOSS to derive rain microphysical process rates and posterior parameter distributions in Bayesian experiments using Markov chain Monte Carlo sampling constrained by synthetic observations.

58 GEOSCIENCES↗
DOE OSTI2020

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 introduced 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 further benefit of rigorously quantified uncertainty of these estimates.

58 GEOSCIENCES↗
DOE OSTI2020

Surprising similarities in model and observational aerosol radiative forcing estimates

Abstract. The radiative forcing from aerosols (particularly through their interaction with clouds) remains one of the most uncertain components of the human forcing of the climate. Observation-based studies have typically found a smaller aerosol effective radiative forcing than in model simulations and were given preferential weighting in the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5). With their own sources of uncertainty, it is not clear that observation-based estimates are more reliable. Understanding the source of the model and observational differences is thus vital to reduce uncertainty in the impact of aerosols on the climate. These reported discrepancies arise from the different methods of separating the components of aerosol forcing used in model and observational studies. Applying the observational decomposition to global climate model (GCM) output, the two different lines of evidence are surprisingly similar, with a much better agreement on the magnitude of aerosol impacts on cloud properties. Cloud adjustments remain a significant source of uncertainty, particularly for ice clouds. However, they are consistent with the uncertainty from observation-based methods, with the liquid water path adjustment usually enhancing the Twomey effect by less than 50 %. Depending on different sets of assumptions, this work suggests that model and observation-based estimates could be more equally weighted in future synthesis studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Sensitivity of Simulated Deep Convection to a Stochastic Ice Microphysics Framework

Ice microphysics parameterizations in models must make major simplifications relative to observations, typically employing empirical relationships to represent average functional properties of particles. However, previous studies have established that ice particle properties vary even in similar cloud types and thermodynamic environments, and it remains unclear how this so-called “natural variability” impacts simulated deep convection. This uncertainty is addressed by implementing a stochastic framework into the Predicted Particle Properties microphysics scheme in the Weather Research and Forecasting model. The approach stochastically varies the coefficients of the mass-size (m-D) relationship (m=aD b ) for unrimed and partially rimed ice. Using guidance from aircraft in situ measurements obtained during the Midlatitude Continental Convective Clouds Experiment (MC3E), the scheme samples from distributions of the prefactor (a) and the exponent (b) of the m-D relationship. Simulations of two MC3E deep convective cases indicate that the stochastic m-D scheme produces considerable variability of anvil cirrus cloud optical depth (τ) distributions, even for the same ice water path (IWP). Thus, the stochastic scheme produces variable cloud radiative forcing that is independent of IWP. This τ-IWP relationship variability is nonexistent using the deterministic m-D ensemble. Additional sensitivity tests are performed in which the fallspeed-size relationship (V=cD d ) is stochastically varied, resulting in variable precipitation amounts and rain rate distributions. Results are presented in the context of satellite and precipitation observations and include comparison with other ensemble configurations using perturbed initial and lateral boundary conditions and small-amplitude noise added to the potential temperature field.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

The Influence of Vertical Wind Shear on Moist Thermals

Although it is well established that vertical wind shear helps to organize and maintain convective systems, there is a longstanding colloquial notion that it inhibits the development of deep convection. To investigate this idea, the vertical momentum budgets of sheared and unsheared moist thermals were compared in idealized cloud model simulations. Consistent with the idea of vertical wind shear inhibiting convective development, convection generally deepened at a slower rate in sheared simulations than in unsheared simulations, and the termination heights of thermals in sheared runs were correspondingly lower. These differences in deepening rates resulted from weaker vertical acceleration of thermals in the sheared compared to the unsheared runs. Downward-oriented dynamic pressure acceleration was enhanced by vertical wind shear, which was the primary reason for relatively weak upward acceleration of sheared thermals. This result contrasts with previous ideas that entrainment or buoyant perturbation pressure accelerations are the primary factors inhibiting the growth of sheared convection. A composite thermal analysis indicates that enhancement of dynamic pressure acceleration in the sheared runs is caused by asymmetric aerodynamic lift forces associated with shear-driven cross flow perpendicular to the direction of the thermals’ ascent. These results provide a plausible explanation for why convection is slower to deepen in sheared environments and why slanted convection tends to be weaker than upright convection in squall lines.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Modeling of Cloud Microphysics: Can We Do Better?

Representation of cloud microphysics is a key aspect of simulating clouds. From the early days of cloud modeling, numerical models have relied on an Eulerian approach for all cloud and thermodynamic and microphysics variables. Over time the sophistication of microphysics schemes has steadily increased, from simple representations of bulk masses of cloud and rain in each grid cell, to including different ice particle types and bulk hydrometeor concentrations, to complex schemes referred to as bin or spectral schemes that explicitly evolve the hydrometeor size distributions within each model grid cell. As computational resources grow, there is a clear trend toward wider use of bin schemes, including their use as benchmarks to develop and test simplified bulk schemes. We argue that continuing on this path brings fundamental challenges difficult to overcome. The Lagrangian particle-based probabilistic approach is a practical alternative in which the myriad of cloud and precipitation particles present in a natural cloud is represented by a judiciously selected ensemble of point particles called superdroplets or superparticles. Here, the advantages of the Lagrangian particle-based approach when compared to the Eulerian bin methodology are explained, and the prospects of applying the method to more comprehensive cloud simulations—for instance, targeting deep convection or frontal cloud systems—are discussed.

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

Assessing the Influence of Microphysical and Environmental Parameter Perturbations on Orographic Precipitation

Microphysical (MP) schemes contain parameters whose values can impact the amount and location of forecasted precipitation, and sensitivity is typically explored by perturbing one parameter at a time while holding the rest constant. Although much can be learned from these “one-at-a-time” studies, the results are limited as these methods do not allow for nonlinear interactions of multiple perturbed parameters. This study applies the Morris one-at-a-time (MOAT) method, a robust statistical tool allowing for simultaneous perturbation of numerous parameters, to explore orographic precipitation sensitivity to changes in microphysical and environmental parameters within an environment characteristic of an atmospheric river. Results show parameters associated with snow fall speed coefficient A s and density ρ s , ice-cloud water collection efficiency (ECI), rain accretion (WRA), relative humidity, zonal wind speed, and surface potential temperature cause the largest influence on simulated precipitation. MP and environmental parameter perturbations can cause precipitation changes of similar magnitude, but results vary by location on the mountain. Different environments are also tested, with A s being the most influential MP parameter regardless of environment. Fewer MP parameters influence precipitation in a faster-wind-speed environment, possibly due to the stronger dynamical forcing upwind and different wave dynamics downwind, compared to a slower-wind-speed environment. Finally, perturbing MP parameters within a single scheme can result in precipitation variations of similar magnitude compared to using entirely different microphysics schemes. MOAT results presented in this study have implications for Bayesian parameter estimation methods and stochastic parameterization within ensemble forecasting.

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

A novel approach for characterizing the variability in mass–dimension relationships: results from MC3E

Abstract. Mass–dimension (m–D) relationships determining bulkmicrophysical properties such as total water content (TWC) and radarreflectivity factor (Z) from particle size distributions are used in bothnumerical models and remote sensing retrievals. The a and b coefficientsrepresenting m=aD b relationships, however, can vary significantlydepending on meteorological conditions, particle habits, the definition ofparticle maximum dimension, the probes used to obtain the data, techniquesused to process the cloud probe data, and other unknown reasons. Thus,considering a range of a,b coefficients may be more applicable for use innumerical models and remote sensing retrievals. Microphysical data collectedby two-dimensional optical array probes (OAPs) installed on the University ofNorth Dakota (UND) Citation aircraft during the Mid-latitude Continental ConvectiveClouds Experiment (MC3E) were used in conjunction with TWC data from aNevzorov probe and ground-based S-band radar data to determine a and busing a technique that minimizes the chi-square difference between the TWC andZ derived from the OAPs and those directly measured by a TWC probe andradar. All a and b values within a specified tolerance were regarded as equallyplausible solutions. Of the 16 near-constant-temperature flight legs analyzedduring the 25 April, 20 May, and 23 May 2011 events, the derived surfaces ofsolutions on the first 2 days where the aircraft-sampled stratiform cloudhad a larger range in a and b for lower temperature environments thatcorrespond to less variability in N(D), TWC, and Z for a flight leg.Because different regions of the storm were sampled on 23 May, differences inthe variability in N(D), TWC, and Z influenced the distribution ofchi-square values in the (a,b) phase space and the specified tolerance in away that yielded 2.8 times fewer plausible solutions compared to the flightlegs on the other dates. These findings show the importance of representingthe variability in a,b coefficients for numerical modeling and remotesensing studies, rather than assuming fixed values, as well as the need tofurther explore how these surfaces depend on environmental conditions inclouds containing ice hydrometeors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

A General N -Moment Normalization Method for Deriving Raindrop Size Distribution Scaling Relationships

A general drop size distribution (DSD) normalization method is formulated in terms of generalized power series relating any DSD moment to any number and combination of reference moments. This provides a consistent framework for comparing the variability of normalized DSD moments using different sets of reference moments, with no explicit assumptions about the DSD functional form (e.g., gamma). It also provides a method to derive any unknown moment plus an estimate of its uncertainty from one or more known moments, which is relevant to remote sensing retrievals and bulk microphysics schemes in weather and climate models. Herein, the approach is applied to a large dataset of disdrometer-observed and bin microphysics-modeled DSDs. As expected, the spread of normalized moments decreases as the number of reference moments is increased, quantified by the logarithmic standard deviation of the normalized moments, σ. Averaging σ for all combinations of reference moments and normalized moments of integer order 0–10, 42.9%, 81.3%, 93.7%, and 96.9% of spread are accounted for applying one-, two-, three-, and four-moment normalizations, respectively. Thus, DSDs can be well characterized overall using three reference moments, whereas adding a fourth reference moment contributes little independent information. The spread of disdrometer-observed DSD moments from uncertainty associated with drop count statistics generally lies between values of σ using two- and three-moment normalizations. However, this uncertainty has little impact on the derived DSD scaling relationships or σ when considered.

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

Broadening of Modeled Cloud Droplet Spectra Using Bin Microphysics in an Eulerian Spatial Domain

This study investigates droplet size distribution (DSD) characteristics from condensational growth and transport in Eulerian dynamical models with bin microphysics. A hierarchy of modeling frameworks is utilized, including parcel, one-dimensional (1D), and three-dimensional large-eddy simulation (LES). The bin DSDs from the 1D model, which includes only vertical advection and condensational growth, are nearly as broad as those from the LES and in line with observed DSD widths for stratocumulus clouds. These DSDs are much broader than those from Lagrangian microphysical calculations within a parcel framework that serve as a numerical benchmark for the 1D tests. In contrast, the bin-modeled DSDs are similar to the Lagrangian microphysical benchmark for a rising parcel in which Eulerian transport is not considered. These results indicate that numerical diffusion associated with vertical advection is a key contributor to broadening DSDs in the 1D model and LES. This DSD broadening from vertical numerical diffusion is unphysical, in contrast to the physical mixing processes that previous studies have indicated broaden DSDs in real clouds. It is proposed that artificial DSD broadening from vertical numerical diffusion compensates for underrepresented horizontal variability and mixing of different droplet populations in typical LES configurations with bin microphysics, or the neglect of other mechanisms that broaden DSDs such as growth of giant cloud condensation nuclei. These results call into question the ability of Eulerian dynamical models with bin microphysics to investigate the physical mechanisms for DSD broadening, even though they may reasonably simulate overall DSD characteristics.

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

Future Community Efforts in Understanding and Modeling Atmospheric Processes

Atmospheric processes are important for weather, climate, and water prediction, as well as climate change projection, but their representation through subgrid parameterizations remains challenging. Understanding and representing these processes are therefore major goals of international research programs. The Global Atmospheric System Studies (GASS) Panel (Zeng and Klocke 2017), as part of the Global Energy and Water Exchanges Project (GEWEX), is dedicated to developing and improving the representation of the atmosphere in weather and climate models by facilitating and supporting the international community to carry out and use observations, process studies, and numerical model experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Impacts of Ice Particle Shape and Density Evolution on the Distribution of Orographic Precipitation

Abstract An IMPROVE-2 orographic precipitation case is simulated using the Ice-Spheroids Habit Model with Aspect-Ratio Evolution (ISHMAEL) microphysics. In ISHMAEL, the evolution of ice particle properties such as mass, shape, size, density, and fall speed are predicted. These ice particle properties along with the ice size distributions from ISHMAEL and model-derived spatial distribution of accumulated precipitation are compared to observations. ISHMAEL predicts planar and columnar particles at spatial locations that agree with observations. Sensitivity simulations are used to explore the impact of predicting ice particle shape evolution on orographic cloud properties and precipitation compared to the traditional approach of representing snow and graupel using separate categories with conversion from snow to graupel during riming. High biases in both IWCs aloft and surface precipitation accumulation occur in the Umpqua River valley using separate snow and graupel categories because snow that does not convert to graupel is advected over the Coast Range and precipitates out in the valley. Improvements in IWCs aloft and surface precipitation using ISHMAEL occur from both predicting various vapor-grown habits and predicting the impact of partial riming on ice particle properties. Compared to traditional microphysics schemes, ISHMAEL also produces less spatial variability in accumulated precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Effects of Scavenging, Entrainment, and Aqueous Chemistry on Peroxides and Formaldehyde in Deep Convective Outflow Over the Central and Southeast United States

Deep convective transport of gaseous precursors to ozone (O3) and aerosols to the upper troposphere is affected by liquid phase and mixed-phase scavenging, entrainment of free tropospheric air and aqueous chemistry. The contributions of these processes are examined using aircraft measurements obtained in storm inflow and outflow during the 2012 Deep Convective Clouds and Chemistry (DC3) experiment combined with high-resolution (dx = 3 km) WRF-Chem simulations of a severe storm, an air mass storm, and a mesoscale convective system (MCS). The simulation results for the MCS suggest that formaldehyde (CH2O) is not retained in ice when cloud water freezes, in agreement with previous studies of the severe storm. By analyzing WRF-Chem trajectories, the effects of scavenging, entrainment, and aqueous chemistry on outflow mixing ratios of CH2O, methyl hydroperoxide (CH3OOH), and hydrogen peroxide (H2O2) are quantified. Liquid phase microphysical scavenging was the dominant process reducing CH2O and H2O2 outflow mixing ratios in all three storms. Aqueous chemistry did not significantly affect outflow mixing ratios of all three species. In the severe storm and MCS, the higher than expected reductions in CH3OOH mixing ratios in the storm cores were primarily due to entrainment of low-background CH3OOH. In the air mass storm, lower CH3OOH and H2O2 scavenging efficiencies (SEs) than in the MCS were partly due to entrainment of higher background CH3OOH and H2O2. Overestimated rain and hail production in WRF-Chem reduces the confidence in ice retention fraction values determined for the peroxides and CH2O.

Bela, Megan↗
DOE OSTI2018

Sensitivity of Mountain Hydroclimate Simulations in Variable-Resolution CESM to Microphysics and Horizontal Resolution

Mountains are natural dams that impede atmospheric moisture transport and water towers that cool, condense, and store precipitation. They are essential in the western United States where precipitation is seasonal, and snowpack is needed to meet water demand. With anthropogenic climate change increasingly threatening mountain snowpack, there is a pressing need to better understand the driving climatological processes. However, the coarse resolution typical of modern global climate models renders them largely insufficient for this task, and signals a need for an advanced strategy. This paper continues the assessment of variable-resolution in the Community Earth System Model (VR-CESM) in modeling mountain hydroclimatology to understand the role of grid-spacing at 55, 28, 14, and 7 km and microphysics, specifically the Morrison and Gettelman (2008, MG1, https://doi.org/10.1175/2008JCLI2105.1) scheme versus the Gettelman and Morrison (2015, MG2, https://doi.org/10.1175/JCLI-D-14-00102.1) scheme. Eight VR-CESM simulations were performed from 1999 to 2015 with the F_AMIP_CAM5 component set, which couples the atmosphere-land models and prescribes ocean data. Refining horizontal grid-spacing from 28 to 7 km with the MG1 scheme did not improve the simulated mountain hydroclimatology. Substantial improvements occurred with the use of MG2 at grid-spacings ≤28 km compared to MG1 as shown with subsequent statistics. Average SWE bias diminished by 9.4X, 4.9X, and 3.5X from 55 to 7 km. The range in minimum (maximum) DJF spatial correlations increased by 0.1–0.2 in both precipitation and SWE. Mountain windward/leeward distributions and elevation profiles improved across hydroclimate variables, however not always with model resolution alone. Disconcertingly, all VR-CESM simulations exhibited a systemic mountain cold bias that worsened with elevation and will require further examination.

58 GEOSCIENCES↗
DOE OSTI2018

Orographic Precipitation Response to Microphysical Parameter Perturbations for Idealized Moist Nearly Neutral Flow

Abstract This study explores the sensitivity of clouds and precipitation to microphysical parameter perturbations using idealized simulations of moist, nearly neutral flow over a bell-shaped mountain. Numerous parameters are perturbed within the Morrison microphysics scheme. The parameters that most affect cloud and precipitation characteristics are the snow fall speed coefficient As, snow particle density ρs, rain accretion (WRA), and ice–cloud water collection efficiency (ECI). Surface precipitation rates are affected by As and ρs through changes to the precipitation efficiency caused by direct and indirect impacts on snow fall speed, respectively. WRA and ECI both affect the amount of cloud water removed, but the precipitation sensitivity differs. Large WRA results in increased precipitation efficiency and cloud water removal below the freezing level, indirectly decreasing cloud condensation rates; the net result is little precipitation sensitivity. Large ECI removes cloud water above the freezing level but with little influence on overall condensation rates. Two environmental experiments are performed to test the robustness of the results: 1) reduction of the wind speed profile by 30% (LowU) and 2) decreasing the surface potential temperature to induce a freezing level below the mountain top (LowFL). Parameter perturbations within LowU result in similar mechanisms acting on precipitation, but a much weaker sensitivity compared to the control. The LowFL case shows ρs is no longer a dominant parameter and As now induces changes to cloud condensation, since more of the cloud depth is present above the freezing level. In general, perturbations to microphysical parameters affect the location of peak precipitation, while the total amount of precipitation is more sensitive to environmental parameter perturbations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Theoretical Expressions for the Ascent Rate of Moist Deep Convective Thermals

An approximate analytic expression is derived for the ratio λ of the ascent rate of moist deep convective thermals and the maximum vertical velocity within them; λ is characterized as a function of two nondimensional buoyancy-dependent parameters y and h and is used to express the thermal ascent rate as a function of the buoyancy field. The parameter y characterizes the vertical distribution of buoyancy within the thermal, and h is the ratio of the vertically integrated buoyancy from the surface to the thermal top and the vertical integral of buoyancy within the thermal. Theoretical λ values are calculated using values of y and h obtained from idealized numerical simulations of ascending moist updrafts and compared to λ computed directly from the simulations. The theoretical values of [Formula: see text] 0.4–0.8 are in reasonable agreement with the simulated λ (correlation coefficient of 0.86). These values are notably larger than the [Formula: see text] from Hill’s (nonbuoyant) analytic spherical vortex, which has been used previously as a framework for understanding the dynamics of moist convective thermals. The relatively large values of λ are a result of net positive buoyancy within the upper part of thermals that opposes the downward-directed dynamic pressure gradient force below the thermal top. These results suggest that nonzero buoyancy within moist convective thermals, relative to their environment, fundamentally alters the relationship between the maximum vertical velocity and the thermal-top ascent rate compared to nonbuoyant vortices. Implications for convection parameterizations and interpretation of the forces contributing to thermal drag are discussed.

54 ENVIRONMENTAL SCIENCES↗