Search NASA⌕ Search

SEARCH · Search NASA

Results for “updraft modeling”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

[Relationship Between Core Convective Structure and Intensity Change in Tropical Cyclones]

Tropical Storm Chantal during August 2001 was a storm that failed to intensify over the few days prior to making landfall on the Yucatan Peninsula. An observational study of Tropical Storm Chantal is presented using a diverse data set including remote and in situ measurements from the NASA ER-2 and DC-8 and the NOAA WP-3D N42RF aircraft and satellite data. The authors discuss the storm structure from the larger scale environment down to the convective scale. Large vertical shear (850-200 hPa shear magnitude range 8-15 m/s) plays a very important role in preventing Chantal from intensifying. The storm had a poorly defined vortex that only extended up to 5-6 km altitude, and an adjacent intense convective region that comprised an Mesoscale Convective System (MCS). The entire low-level circulation center was in the rain-free western side of the storm, about 80 km to the west-southwest of the MCS. The MCS appears to have been primarily the result of intense convergence between large scale, low-level easterly flow with embedded downdrafts, and the cyclonic vortex flow. The individual cells in the MCS such as Cell 2 during the period of the observations, were extremely intense with reflectivity core diameters of 10 km and peak updrafts exceeding 20 m/s. Associated with this MCS were two broad subsidence (warm) regions both of which had portions over the vortex. The first layer near 700 hPa was directly above the vortex and covered most of it. The second layer near 500 hPa was along the forward and right flanks of Cell 2 and undercut the anvil divergence region above. There was not much resemblance of these subsidence layers to typical upper level warm cores in hurricanes that are necessary to support strong surface winds and a low central pressure. The observations are compared to previous studies of weakly sheared storms and modeling studies of shear effects and intensification. The configuration of the convective updrafts, low-level circulation, and lack of vertical coherence between the upper and low level warming regions, likely inhibited intensification of Chantal. This configuration is consistent with modeling of vortices in sheared environments, which suggest strongest convection and rain in the downshear left quadrant of the storm, and subsidence in the upshear right quadrant. The vertical shear profile is however different from what was assumed in previous modeling in that the winds are strongest in the lowest levels and the deep tropospheric vertical shear is on the order of 10-12 m/s.

Source record↗
NASA NTRS2006

Structure of Highly Sheared Tropical Storm Chantal during CAMEX-4

Tropical Storm Chantal during August 2001 was a storm that failed to intensify over the few days prior to making landfall on the Yucatan Peninsula. An observational study of Tropical Storm Chantal is presented using a diverse dataset including remote and in situ measurements from the NASA ER-2 and DC-8 and the NOAA WP-3D N42RF aircraft and satellite. The authors discuss the storm structure from the larger-scale environment down to the convective scale. Large vertical shear (850-200-hPa shear magnitude range 8-15 m/s) plays a very important role in preventing Chantal from intensifying. The storm had a poorly defined vortex that only extended up to 5-6-km altitude, and an adjacent intense convective region that comprised a mesoscale convective system (MCS). The entire low-level circulation center was in the rain-free western side of the storm, about 80 km to the west-southwest of the MCS. The MCS appears to have been primarily the result of intense convergence between large-scale, low-level easterly flow with embedded downdrafts, and the cyclonic vortex flow. The individual cells in the MCS such as cell 2 during the period of the observations were extremely intense, with reflectivity core diameters of 10 km and peak updrafts exceeding 20 m/s. Associated with this MCS were two broad subsidence (warm) regions, both of which had portions over the vortex. The first layer near 700 hPa was directly above the vortex and covered most of it. The second layer near 500 hPa was along the forward and right flanks of cell 2 and undercut the anvil divergence region above. There was not much resemblance of these subsidence layers to typical upper-level warm cores in hurricanes that are necessary to support strong surface winds and a low central pressure. The observations are compared to previous studies of weakly sheared storms and modeling studies of shear effects and intensification. The configuration of the convective updrafts, low-level circulation, and lack of vertical coherence between the upper- and lower-level warming regions likely inhibited intensification of Chantal. This configuration is consistent with modeled vortices in sheared environments, which suggest the strongest convection and rain in the downshear left quadrant of the storm, and subsidence in the upshear right quadrant. The vertical shear profile is, however, different from what was assumed in previous modeling in that the winds are strongest in the lowest levels and the deep tropospheric vertical shear is on the order of 10-12 m/s.

Heymsfield, G. M.↗
NASA NTRS2014

Physical and Dynamical Linkages Between Lightning Jumps and Storm Conceptual Models

The presence and rates of total lightning are both correlated to and physically dependent upon storm updraft strength, mixed phase precipitation volume and the size of the charging zone. The updraft modulates the ingredients necessary for electrification within a thunderstorm, while the updraft also plays a critical role in the development of severe and hazardous weather. Therefore utilizing this relationship, the monitoring of lightning rates and jumps provides an additional piece of information on the evolution of a thunderstorm, more often than not, at higher temporal resolution than current operational radar systems. This correlation is the basis for the total lightning jump algorithm that has been developed in recent years. Currently, the lightning jump algorithm is being tested in two separate but important efforts. Schultz et al. (2014; this conference) is exploring the transition of the algorithm from its research based formulation to a fully objective algorithm that includes storm tracking, Geostationary Lightning Mapper (GLM) Proxy data and the lightning jump algorithm. Chronis et al. (2014) provides context for the transition to current operational forecasting using lightning mapping array based products. However, what remains is an end-to-end physical and dynamical basis for coupling total lightning flash rates to severe storm manifestation, so the forecaster has a reason beyond simple correlation to utilize the lightning jump algorithm within their severe storm conceptual models. Therefore, the physical basis for the lightning jump algorithm in relation to severe storm dynamics and microphysics is a key component that must be further explored. Many radar studies have examined flash rates and their relationship to updraft strength, updraft volume, precipitation-sized ice mass, etc.; however, their relationship specifically to lightning jumps is fragmented within the literature. Thus the goal of this study is to use multiple Doppler and polarimetric radar techniques to resolve the physical and dynamical storm characteristics specifically around the time of the lightning jump. This information will help forecasters anticipate lightning jump occurrence, or even be of use to determine future characteristics of a given storm (e.g., development of a mesocyclone, downdraft, or hail signature on radar), providing additional lead time/confidence in the severe storm warning paradigm.

Schultz, Christopher J.↗
NASA NTRS2014

Integration of the Total Lightning Jump Algorithm into Current Operational Warning Environment Conceptual Models

The presence and rates of total lightning are both correlated to and physically dependent upon storm updraft strength, mixed phase precipitation volume and the size of the charging zone. The updraft modulates the ingredients necessary for electrification within a thunderstorm, while the updraft also plays a critical role in the development of severe and hazardous weather. Therefore utilizing this relationship, the monitoring of lightning rates and jumps provides an additional piece of information on the evolution of a thunderstorm, more often than not, at higher temporal resolution than current operational radar systems. This correlation is the basis for the total lightning jump algorithm that has been developed in recent years. Currently, the lightning jump algorithm is being tested in two separate but important efforts. Schultz et al. (2014; AMS 10th Satellite Symposium) is exploring the transition of the algorithm from its research based formulation to a fully objective algorithm that includes storm tracking, Geostationary Lightning Mapper (GLM) Proxy data and the lightning jump algorithm. Chronis et al. (2014; this conference) provides context for the transition to current operational forecasting using lightning mapping array based products. However, what remains is an end to end physical and dynamical basis for relating lightning rates to severe storm manifestation, so the forecaster has a reason beyond simple correlation to utilize the lightning jump algorithm within their severe storm conceptual models. Therefore, the physical basis for the lightning jump algorithm in relation to severe storm dynamics and microphysics is a key component that must be further explored. Many radar studies have examined flash rates and their relation to updraft strength, updraft volume, precipitation‐sized ice mass, etc.; however, relation specifically to lightning jumps is fragmented within the literature. Thus the goal of this study is to use multiple Doppler techniques to resolve the physical and dynamical storm characteristics specifically around the time of the lightning jump. This information will help forecasters anticipate lightning jump occurrence, or even be of use to determine future characteristics of a given storm (e.g., development of a mesocyclone, downdraft, or hail signature on radar), providing additional lead time/confidence in the severe storm warning paradigm.

Shultz, Christopher J.↗
NASA NTRS2014

Physical and Dynamical Linkages between Lightning Jumps and Storm Conceptual Models

The presence and rates of total lightning are both correlated to and physically dependent upon storm updraft strength, mixed phase precipitation volume and the size of the charging zone. The updraft modulates the ingredients necessary for electrification within a thunderstorm, while the updraft also plays a critical role in the development of severe and hazardous weather. Therefore utilizing this relationship, the monitoring of lightning rates and jumps provides an additional piece of information on the evolution of a thunderstorm, more often than not, at higher temporal resolution than current operational radar systems. This correlation is the basis for the total lightning jump algorithm that has been developed in recent years. Currently, the lightning jump algorithm is being tested in two separate but important efforts. Schultz et al. (2014; this conference) is exploring the transition of the algorithm from its research based formulation to a fully objective algorithm that includes storm tracking, Geostationary Lightning Mapper (GLM) Proxy data and the lightning jump algorithm. Chronis et al. (2014; this conference) provides context for the transition to current operational forecasting using lightning mapping array based products. However, what remains is an end-to-end physical and dynamical basis for coupling total lightning flash rates to severe storm manifestation, so the forecaster has a reason beyond simple correlation to utilize the lightning jump algorithm within their severe storm conceptual models. Therefore, the physical basis for the lightning jump algorithm in relation to severe storm dynamics and microphysics is a key component that must be further explored. Many radar studies have examined flash rates and their relationship to updraft strength, updraft volume, precipitation-sized ice mass, etc.; however, their relationship specifically to lightning jumps is fragmented within the literature. Thus the goal of this study is to use multiple Doppler and polarimetric radar techniques to resolve the physical and dynamical storm characteristics specifically around the time of the lightning jump. This information will help forecasters anticipate lightning jump occurrence, or even be of use to determine future characteristics of a given storm (e.g., development of a mesocyclone, downdraft, or hail signature on radar), providing additional lead time/confidence in the severe storm warning paradigm.

Schultz, Christopher J.↗
NASA NTRS2013

Integration of the Total Lightning Jump Algorithm into Current Operational Warning Environment Conceptual Models

The presence and rates of total lightning are both correlated to and physically dependent upon storm updraft strength, mixed phase precipitation volume and the size of the charging zone. The updraft modulates the ingredients necessary for electrification within a thunderstorm, while the updraft also plays a critical role in the development of severe and hazardous weather. Therefore utilizing this relationship, the monitoring of lightning rates and jumps provides an additional piece of information on the evolution of a thunderstorm, more often than not, at higher temporal resolution than current operational radar systems. This correlation is the basis for the total lightning jump algorithm that has been developed in recent years. In order to become a viable option for operational forecasters to incorporate into their severe storm monitoring process, the total lightning jump must be placed into the framework of several severe storm conceptual models (e.g., radar evolution, storm morphology) which forecasters have built through training and experience. Thus, one of the goals of this study is to examine and relate the lightning jump concept to often used radar parameters (e.g., dBZ vertical structure, VIL, MESH, MESO/shear) in the warning environment. Tying lightning trends and lightning jump occurrences to these radar based parameters will provide forecasters with an additional tool that they can use to build an accurate realtime depiction as to what is going on in a given environment. Furthermore, relating the lightning jump concept to these parameters could also increase confidence in a warning decision they have already made, help tip the scales on whether or not to warn on a given storm, or to draw the forecaster s attention to a particular storm that is rapidly developing. Furthermore the lightning information will add vital storm scale information in regions that are not well covered by radar, or when radar failures occur. The physical basis for the lightning jump algorithm in relation to severe storm dynamics and microphysics is a key component that must be further explored. Many radar studies have examined flash rates and their relation to updraft strength, updraft volume, precipitation -sized ice mass, etc.; however, very few have related the concept of the lightning jump and manifestation of severe weather to storm dynamics and microphysics using multi -Doppler and polarimetric radar techniques. Therefore, the second half of this study will combine the lightning jump algorithm and these radar techniques in order to place the lightning jump concept into a physical and dynamical framework. This analysis includes examining such parameters as mixed phase precipitation volume, charging zone, updraft strength and updraft volume. Such a study should provide increased understanding of and confidence in the strengths and limitations of the lightning jump algorithm in the storm warning process.

Schultz, Christopher J.↗
NASA NTRS1985

A simple model of evaporatively driven downdraft - Application to microburst downdraft

It is pointed out that the physics of a downdraft are very different from that of an updraft. An updraft is virtually saturated, while the downdraft can be unsaturated. The present investigation is concerned with the dynamics of a simple one-dimensional, time-dependent model of a downdraft. Earlier treatments are extended on the basis of a number of considerations, taking into account a distribution of drop sizes, modification of the drop size distribution by evaporation, entrainment of environmental air into the downdraft, and, most importantly, explicit calculation of the intensity of the downdraft. A large number of calculations are performed to study the effects of conditions at the top of the downdraft and in the environment.

Srivastava, R. C.↗
NASA NTRS2014

A Wrf-Chem Flash Rate Parameterization Scheme and LNO(x) Analysis of the 29-30 May 2012 Convective Event in Oklahoma During DC3

The Deep Convective Clouds and Chemistry (DC3) field campaign in 2012 provided a plethora of aircraft and ground-based observations (e.g., trace gases, lightning and radar) to study deep convective storms, their convective transport of trace gases, and associated lightning occurrence and production of nitrogen oxides (NOx). Based on the measurements taken of the 29-30 May 2012 Oklahoma thunderstorm, an analysis against a Weather Research and Forecasting Chemistry (WRF-Chem) model simulation of the same event at 3-km horizontal resolution was performed. One of the main objectives was to include various flash rate parameterization schemes (FRPSs) in the model and identify which scheme(s) best captured the flash rates observed by the National Lightning Detection Network (NLDN) and Oklahoma Lightning Mapping Array (LMA). The comparison indicates how well the schemes predicted the timing, location, and number of lightning flashes. The FRPSs implemented in the model were based on the simulated thunderstorms physical features, such as maximum vertical velocity, cloud top height, and updraft volume. Adjustment factors were added to each FRPS to best capture the observed flash trend and a sensitivity study was performed to compare the range in model-simulated lightning-generated nitrogen oxides (LNOx) generated by each FRPS over the storms lifetime. Based on the best FRPS, model-simulated LNOx was compared against aircraft measured NOx. The trace gas analysis, along with the increased detail in the model specification of the vertical distribution of lightning flashes as suggested by the LMA data, provide guidance in determining the scenario of NO production per intracloud and cloud-to-ground flash that best matches the NOx mixing ratios observed by the aircraft.

Cloud physics↗
NASA NTRS2014

A WRF-Chem Flash Rate Parameterization Scheme and LNOx Analysis of the 29-30 May 2012 Convective Event in Oklahoma During DC3

The Deep Convective Clouds and Chemistry (DC3) field campaign in 2012 provided a plethora of aircraft and ground-based observations (e.g., trace gases, lightning and radar) to study deep convective storms, their convective transport of trace gases, and associated lightning occurrence and production of nitrogen oxides (NOx). Based on the measurements taken of the 29-30 May 2012 Oklahoma thunderstorm, an analysis against a Weather Research and Forecasting Chemistry (WRF-Chem) model simulation of the same event at 3-km horizontal resolution was performed. One of the main objectives was to include various flash rate parameterization schemes (FRPSs) in the model and identify which scheme(s) best captured the flash rates observed by the National Lightning Detection Network (NLDN) and Oklahoma Lightning Mapping Array (LMA). The comparison indicates how well the schemes predicted the timing, location, and number of lightning flashes. The FRPSs implemented in the model were based on the simulated thunderstorms physical features, such as maximum vertical velocity, cloud top height, and updraft volume. Adjustment factors were applied to each FRPS to best capture the observed flash trend and a sensitivity study was performed to compare the range in model-simulated lightning-generated nitrogen oxides (LNOx) generated by each FRPS over the storms lifetime. Based on the best FRPS, model-simulated LNOx was compared against aircraft measured NOx. The trace gas analysis, along with the increased detail in the model specification of the vertical distribution of lightning flashes as suggested by the LMA data, provide guidance in determining the scenario of NO production per intracloud and cloud-to-ground flash that best matches the NOx mixing ratios observed by the aircraft.

Atmospheric chemistry↗
NASA NTRSDate not supplied

Development of a Unified SHOC-EDMF PBL Parameterization in the NASA GEOS Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies.

SHOC-EDMF↗
NASA NTRSDate not supplied

A Unified SHOC-MF Boundary Layer Parameterization in the NASA Goddard Earth Observing System Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC; Bogenschutz and Krueger, 2013). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies. Results indicate that many higher order moments are underestimated when the mass flux contribution is neglected.

Nathan Arnold↗
NASA NTRSDate not supplied

A Unified SHOC-MF Boundary Layer Parameterization in the NASA Goddard Earth Observing System Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC; Bogenschutz and Krueger, 2013). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies. Results indicate that many higher order moments are underestimated when the mass flux contribution is neglected.

Nathan Arnold↗
NASA NTRS1998

Pulsing Inertial Oscillation, Supercell Storms, and Surface Mesonetwork Data

The pulsing inertial oscillation (PIO) model is a nonlinear, time-dependent, translating vortex solution of the inviscid, compressible fluid dynamic equations in the middle troposphere. The translation of this vortex during a pulse is strikingly similar to that of a supercell storm -- a rotating thunderstorm that can generate tornadoes and hail. Two studies were performed to test the hypothesis that some supercell storms are manifestations of a PIO pulse. The first study applied the model to an intense interior draft whose buoyancy was bounded by a temperature excess of +/- 12 K. The peak updraft speed achieved was 41.5 m/ s and the peak Rossby number was 92.9. The study also pointed to an advanced concept for attaining higher values. The second study applied the PIO model to a supercell storm as a whole and succeeded in replicating its bulk properties, such as mesocyclonic circulation, net mass and moisture influxes, and time track. This study also identified a critical feature of the PIO model that could be tested against storm data: The average vertical draft is downward before the turn in the storm track and upward afterwards. In the conventional theory, the average vertical draft is upward from storm inception until dissipation. These differing draft predictions were compared with the best available data, which are surface mesonetwork data. These data were found to support the PIO model. However, surface data alone are not conclusive, and further measurements are warranted.

Costen, R. C.↗
NASA NTRS1985

Lidar observations of the planetary boundary layer

The application of an airborne downward-looking lidar to the study of organized cellular convection in the planetary boundary layer (PBL) over the ocean is described. The lidar consisted of a frequency doubled Nd-YAG 530 mm-wavelength laser whose axis was aligned colinearly with the optical axis of an all-reflecting 40 mm-diameter Newtonian telescope. The airborne lidar provided a unique observation of both microscale and mesoscale variations of the PBL top. The lidar data, presented as constant backscatter isopleth soundings, provide a visual indication of the presence of vertically organized convection cells. Comparisons of the lidar-derived PBL structure with both a conceptual model of the PBL and laboratory simulations of Deardorf et al. (1980) of a developing convective PBL showed that the observations are consistent with a model of mixing in the PBL, which involves a field of organized updrafts separated by downdrafts.

Melfi, S. H.↗
NASA NTRS1987

The terminal area simulation system. Volume 2: Verification cases

The numerical simulation of five case studies are presented and are compared with available data in order to verify the three-dimensional version of the Terminal Area Simulation System (TASS). A spectrum of convective storm types are selected for the case studies. Included are: a High-Plains supercell hailstorm, a small and relatively short-lived High-Plains cumulonimbus, a convective storm which produced the 2 August 1985 DFW microburst, a South Florida convective complex, and a tornadic Oklahoma thunderstorm. For each of the cases the model results compared reasonably well with observed data. In the simulations of the supercell storms many of their characteristic features were modeled, such as the hook echo, BWER, mesocyclone, gust fronts, giant persistent updraft, wall cloud, flanking-line towers, anvil and radar reflectivity overhang, and rightward veering in the storm propagation. In the simulation of the tornadic storm a horseshoe-shaped updraft configuration and cyclic changes in storm intensity and structure were noted. The simulation of the DFW microburst agreed remarkably well with sparse observed data. The simulated outflow rapidly expanded in a nearly symmetrical pattern and was associated with a ringvortex. A South Florida convective complex was simulated and contained updrafts and downdrafts in the form of discrete bubbles. The numerical simulations, in all cases, always remained stable and bounded with no anomalous trends.

Proctor, F. H.↗
NASA NTRS1994

Roles of drizzle in a one-dimensional third-order turbulence closure model of the nocturnal stratus-topped marine boundary layer

This study focuses on the effects of drizzle in a one-dimensional third-order turbulence closure model of the nocturnal stratus-topped marine boundary layer. When the simulated drizzle rate is relatively small (maximum approximately equal to 0.6 mm/day), steady-state solutions are obtained. The boundary layer stabilizes essentially because drizzle causes evaporative cooling of the subcloud layer. This stabilization considerably reduces the buoyancy flux and turbulence kinetic energy below the stratus cloud. Thus, drizzle tends to decouple the cloud from the subcloud layer in the model, as suggested by many observational studies. In addition, the evaporation of drizzle in the subcloud layer creates small scattered clouds, which are likely to represent cumulus clouds, below the solid stratus cloud in the model. The sensitivity experiments show that these scattered clouds help maintain a coupled boundary layer. When the drizzle rate is relatively large (maximum approximately equal to 0.9 mm/day), the response of the model becomes transient with bursts in turbulent fluxes. This phenomenon is related to the formation of the scattered cloud layer below the solid stratus cloud. It appears that the model is inadequate to represent the heat and moisture transport by strong updrafts covering a small fractional area in cumulus convection.

Wang, Shouping↗
NASA NTRS1986

Thunderstorm cloud top dynamics as inferred from satellite observations and a cloud top parcel model

Thunderstorm top temperature-height structure and temperature-height conversion relations in the convective overshooting top region are examined using data from GOES IR, AVHRR, and aircraft overflights, and GOES stereoscopic observations. Using a Lagrangian one-dimensional cloud parcel model applied in the overshooting region, three classes of storm tops are described and the causes of their unique temperature-height structures are explored. Although the parcel model used does not fully reproduce the three-dimensional structure and dynamics of an evolving thunderstorm top, it is shown to crudely reproduce the parcel dynamics of updraft core parcels. The locations of cold points, high points, warm points, and the magnitude of cold-high offsets compare favorably between the model and the satellite observations.

Adler, Robert F.↗
NASA NTRS2009

Forecasting Lightning Threat using Cloud-resolving Model Simulations

As numerical forecasts capable of resolving individual convective clouds become more common, it is of interest to see if quantitative forecasts of lightning flash rate density are possible, based on fields computed by the numerical model. Previous observational research has shown robust relationships between observed lightning flash rates and inferred updraft and large precipitation ice fields in the mixed phase regions of storms, and that these relationships might allow simulated fields to serve as proxies for lightning flash rate density. It is shown in this paper that two simple proxy fields do indeed provide reasonable and cost-effective bases for creating time-evolving maps of predicted lightning flash rate density, judging from a series of diverse simulation case study events in North Alabama for which Lightning Mapping Array data provide ground truth. One method is based on the product of upward velocity and the mixing ratio of precipitating ice hydrometeors, modeled as graupel only, in the mixed phase region of storms at the -15\dgc\ level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domainwide statistics of the peak values of simulated flash rate proxy fields against domainwide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. A blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Weather Research and Forecast Model simulations of selected North Alabama cases show that this model can distinguish the general character and intensity of most convective events, and that the proposed methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because models tend to have more difficulty in correctly predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models, the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of cloud-allowing forecasts become available.

McCaul, E. W., Jr.↗