Search NASA⌕ Search

SEARCH · Search NASA

Results for “storm resolving model”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Cloud-to-Ground Lightning Estimates Derived from SSMI Microwave Remote Sensing and NLDN

Lightning observations are collected using ground-based and satellite-based sensors. The National Lightning Detection Network (NLDN) in the United States uses multiple ground sensors to triangulate the electromagnetic signals created when lightning strikes the Earth's surface. Satellite-based lightning observations have been made from 1998 to present using the Lightning Imaging Sensor (LIS) on the NASA Tropical Rainfall Measuring Mission (TRMM) satellite, and from 1995 to 2000 using the Optical Transient Detector (OTD) on the Microlab-1 satellite. Both LIS and OTD are staring imagers that detect lightning as momentary changes in an optical scene. Passive microwave remote sensing (85 and 37 GHz brightness temperatures) from the TRMM Microwave Imager (TMI) has also been used to quantify characteristics of thunderstorms related to lightning. Each lightning detection system has fundamental limitations. TRMM satellite coverage is limited to the tropics and subtropics between 38 deg N and 38 deg S, so lightning at the higher latitudes of the northern and southern hemispheres is not observed. The detection efficiency of NLDN sensors exceeds 95%, but the sensors are only located in the USA. Even if data from other ground-based lightning sensors (World Wide Lightning Location Network, the European Cooperation for Lightning Detection, and Canadian Lightning Detection Network) were combined with TRMM and NLDN, there would be enormous spatial gaps in present-day coverage of lightning. In addition, a globally-complete time history of observed lightning activity is currently not available either, with network coverage and detection efficiencies varying through the years. Previous research using the TRMM LIS and Microwave Imager (TMI) showed that there is a statistically significant correlation between lightning flash rates and passive microwave brightness temperatures. The physical basis for this correlation emerges because lightning in a thunderstorm occurs where ice is first present in the cloud and electric charge separation occurs. These ice particles efficiently scatter the microwave radiation at the 85 and 37 GHz frequencies, thus leading to large brightness temperature depressions. Lightning flash rate is related to the total amount of ice passing through the convective updraft regions of thunderstorms. Confirmation of this relationship using TRMM LIS and TMI data, however, remains constrained to TRMM observational limits of the tropics and subtropics. Satellites from the Defense Meteorology Satellite Program (DMSP) have global coverage and are equipped with passive microwave imagers that, like TMI, observe brightness temperatures at 85 and 37 GHz. Unlike the TRMM satellite, however, DMSP satellites do not have a lightning sensor, and the DMSP microwave data has never been used to derive global lightning. In this presentation, a relationship between DMSP Special Sensor Microwave Imager (SSMI) data and ground-based cloud-to-ground (CG) lightning data from NLDN is investigated to derive a spatially complete time history of CG lightning for the USA study area. This relationship is analogous to the established using TRMM LIS and TMI data. NLDN has the most spatially and temporally complete CG lightning data for the USA, and therefore provides the best opportunity to find geospatially coincident observations with SSMI sensors. The strongest thunderstorms generally have minimum 85 GHz Polarized Corrected brightness Temperatures (PCT) less than 150 K. Archived radar data was used to resolve the spatial extent of the individual storms. NLDN data for that storm spatial extent defined by radar data was used to calculate the CG flash rate for the storm. Similar to results using TRMM sensors, a linear model best explained the relationship between storm-specific CG flash rates and minimum 85 GHz PCT. However, the results in this study apply only to CG lightning. To extend the results to weaker storms, the probability of CG lightning (instead of the flash rate) was calculated for storms having 85 GHz PCT greater than 150 K. NLDN data was used to determine if a CG strike occurred for a storm. This probability of CG lightning was plotted as a function of minimum 85 GHz PCT and minimum 37 GHz PCT. These probabilities were used in conjunction with the linear model to estimate the CG flash rate for weaker storms with minimum 85 GHz PCTs greater than 150 K. Results from the investigation of CG lightning and passive microwave radiation signals agree with the previous research investigating total lightning and brightness temperature. Future work will take the established relationships and apply them to the decades of available DMSP data for the USA to derive a map of CG lightning flash rates. Validation of this method and uncertainty analysis will be done by comparing the derived maps of CG lightning flash rates against existing NLDN maps of CG lightning flash rates.

Winesett, Thomas↗
NASA NTRS2013

Diagnosing Aircraft Icing Potential from Satellite Cloud Retrievals

The threat for aircraft icing in clouds is a significant hazard that routinely impacts aviation operations. Accurate diagnoses and forecasts of aircraft icing conditions requires identifying the location and vertical distribution of clouds with super-cooled liquid water (SLW) droplets, as well as the characteristics of the droplet size distribution. Traditional forecasting methods rely on guidance from numerical models and conventional observations, neither of which currently resolve cloud properties adequately on the optimal scales needed for aviation. Satellite imagers provide measurements over large areas with high spatial resolution that can be interpreted to identify the locations and characteristics of clouds, including features associated with adverse weather and storms. This paper describes new techniques for interpreting cloud products derived from satellite data to infer the flight icing threat to aircraft. For unobscured low clouds, the icing threat is determined using empirical relationships developed from correlations between satellite imager retrievals of liquid water path and droplet size with icing conditions reported by pilots (PIREPS). For deep ice over water cloud systems, ice and liquid water content (IWC and LWC) profiles are derived by using the imager cloud properties to constrain climatological information on cloud vertical structure and water phase obtained apriori from radar and lidar observations, and from cloud model analyses. Retrievals of the SLW content embedded within overlapping clouds are mapped to the icing threat using guidance from an airfoil modeling study. Compared to PIREPS and ground-based icing remote sensing datasets, the satellite icing detection and intensity accuracies are approximately 90% and 70%, respectively, and found to be similar for both low level and deep ice over water cloud systems. The satellite-derived icing boundaries capture the reported altitudes over 90% of the time. Satellite analyses corresponding to the time and location of several recent aviation accidents and with icing PIREPS are also presented that reveal skill in identifying severe icing conditions. These results demonstrate the utility of satellite cloud retrievals for quantitatively diagnosing the potential for icing conditions on temporal and spatial scales that should be useful to the aviation community. Plans are being developed to deliver these new satellite products to the GOES-R Proving Ground in the near future so that they can be evaluated in operational applications.

Smith, William L., Jr.↗
NASA NTRS2004

Sensitivity of a Cloud-Resolving Model to the Bulk and Explicit Bin Microphysical Schemes

A cloud-resolving model is used to study sensitivities of two different microphysical schemes, one is the traditional bulk type, and the other is an explicit bin scheme, in simulating a mid-latitude squall line case (PRE-STORM, June 10-1 1,1985). Simulations using different microphysical schemes are compared with each other and also with the observations. Both the bulk and bin models reproduce the general features during the developing and mature stage of the system. Furthermore, the observations and the well-proven bulk scheme simulation serve as validations for the newly incorporated bin scheme. However, it is also shown that the bulk and bin simulations have distinct differences, most notably in the stratiform region of the squall line system. Weak convective cells exist in the stratiform region in the bulk simulation, but not in the bin simulation. These weak convective cells in the stratiform region simulated in the bulk scheme model are remnants of the stronger convections previously at the leading edge of the system, sustained by horizontal vorticity generated by its own cool pool near the surface. The bin simulation, on the other hand, has a horizontally homogeneous stratiform cloud structure, which agrees better with the observations. Examinations of the downdraft core strength, the potential temperature perturbation, and the evaporative cooling rate show that the differences between the bulk and bin models are due mainly to the stronger low-level evaporative cooling in the convective zone simulated in the bulk microphysical scheme, which is unrealistic because of the assumptions made in raindrop size distribution. Further sensitivity tests that reduce the evaporation rate in bulk scheme artificially produce more upright convective core and less weak cores in stratiform region. However, they produce weaker upper level outflow and consequently less stratiform rain area. The addition of a more realistic raindrop breakup scheme in the bin scheme results more realistic radar reflectivity and stronger surface rainfall. Despite the increase of the rain evaporation and strengthening of the near surface cool pool, bin scheme with rain breakup shows homogeneous stratiform rain. These sensitivity tests prove the robustness of the bin microphysical scheme and the difficulty of tuning the limited parameters in the bulk microphysical scheme to realistically reproduce detail structures in a mid-latitude squall line case study.

Li, Xiao-Wen↗
NASA NTRS1999

The Dependence on Grid Resolution of Numerically Simulated Convective Cloud Systems Using Ice Microphysics

Mesoscale research and forecast models are increasingly being used at horizontal resolutions of 1-8 km to simulate a variety of precipitating systems. When the model is used to simulate convective systems, it is uncertain to what extent the dynamics and microphysics of convective updrafts can be resolved with grids larger than 1 km. In this study, two- and three-dimensional versions of the Goddard Cumulus Ensemble model are used to determine the impact of horizontal grid resolution on the behavior of the simulated storms and on the characteristics of the cloud microphysical fields. It will be shown that as resolution decreases from about 1 km to greater than 3 km, there is a fairly rapid degradation of the storm structure in the form of reduced convective mass fluxes, updraft tilts, and cloud microphysics. A high-resolution simulation of hurricane outer rainbands using the MM5 mesoscale model shows also that there can be a substantial modification of the key microphysical processes that contribute to rainfall as a result of reducing the horizontal resolution.

Braun, Scott A.↗
DOE OSTI2026

Revealing the Formation of the <20 MeV Inner Proton Radiation Belt at L ∼ 2 During the 10–11 May 2024 Superstorm

A new proton radiation belt was identified during the geomagnetic superstorm of 10–11 May 2024. To investigate its origin, we use an MHD-test particle simulation to model solar energetic proton (SEP) trapping and the evolution of the initial trapped proton population during the storm. The simulation weights the injected SEP population with interplanetary proton measurements and the initial trapped proton radial profile with energy-resolved Weather System Follow-on-Microwave (WSF-M) proton data, enabling quantitative comparison with observations. Results show that SEP contributions are minimal; the new ∼2–18 MeV belt at L ∼ 2 forms mainly through redistribution and energization of the initial trapped population within hours of the CME shock, driven by subsequent electric field impulses, rather than the initial shock alone. These results provide quantitative evidence based on realistic, plasma-dependent field dynamics and significantly advance understanding of inner radiation belt formation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2015

Kinematic and Microphysical Control of Lightning Flash Rate over Northern Alabama

The Deep Convective Clouds and Chemistry (DC3) experiment seeks to examine the relationship between deep convection and the production of nitrogen oxides (NO (sub x)) via lightning (LNO (sub x)). A critical step in estimating LNO (sub x) production in a cloud-resolving model (CRM) without explicit lightning is to estimate the flash rate from available model parameters that are statistically and physically correlated. As such, the objective of this study is to develop, improve and evaluate lightning flash rate parameterizations in a variety of meteorological environments and storm types using radar and lightning mapping array (LMA) observations taken over Northern Alabama from 2005-2012, including during DC3. UAH's Advanced Radar for Meteorological and Operational Research (ARMOR) and the Weather Surveillance Radar - 1988 Doppler (WSR 88D) located at Hytop (KHTX) comprises the dual-Doppler and polarimetric radar network, which has been in operation since 2004. The northern Alabama LMA (NA LMA) in conjunction with Vaisala's National Lightning Detection Network (NLDN) allow for a detailed depiction of total lightning during this period. This study will integrate ARMOR-KHTX dual Doppler/polarimetric radar and NA LMA lightning observations from past and ongoing studies, including the more recent DC3 results, over northern Alabama to form a large data set of 15-20 case days and over 20 individual storms, including both ordinary multicell and supercell convection. Several flash rate parameterizations will be developed and tested, including those based on 1) graupel/small hail volume; 2) graupel/small hail mass, and 3) convective updraft volume. Sensitivity of the flash rate parameterizations to storm intensity, storm morphology and environmental conditions will be explored.

Carey, Lawrence D.↗
NASA NTRS1990

Impulsive microwave burst and solar noise storm emission resolved with the VLA

The evolution of a microwave burst at 20.7 cm wavelength and a type I noise storm at 91.6 cm wavelength are examined using VLA. The magnetic loops in the two spectral regions are studied. The sizes and brightness temperatures of the 20.7 cm burst sources are compared with those predicted by multithermal and nonthermal models of microwave burst emission of Dulk and Dennis (1982). The data reveal that: (1) the precursor, impulsive, and post-impulsive phases of the 20.7 cm burst are located in spatially separated sources; (2) the gradual enhancement of a 91.6 cm noise storm continuum source suggests a feed-back mechanism exists between activity in higher and lower lying loops; and (3) the derived magnetic field strengths at 20.7 cm for the nonthermal models are 35-160 G with core field strengths of 125-220 G.

Willson, Robert F.↗
NASA NTRSDate not supplied

Mars Global Climate Modeling.

Scientists use Global Climate Models (GCMs) to better under the current and past climate states of terrestrial (solid surface) bodies in our solar system and beyond, and the physical processes that control them. GCMs are complex, multi-dimensional computer codes that can generally be divided into two parts: the geophysical fluid dynamics (GFD) framework, which represents accelerations and spatially resolved processes, and the physics routines, which provide the forcing functions for the circulation. Producing a GCM that is appropriate for a particular body—Mars, for example—requires implementing the appropriate physics routines (e.g., radiative transfer, planetary boundary layer physics, dust lifting physics to generate dust storms, etc.) for that body. In this talk, I will give an overview of the components of the NASA Ames Mars GCM and discuss some of the scientific questions we address with this state-of-the-art numerical model.

Melinda A. Kahre↗
NASA NTRS2003

Cloud-Resolving Model and GPM

Over the past twenty years, rainfall retrieval algorithms have been developed to retrieve rainfall and vertical hydrometeor structures from passive microwave observations by making use of the fact that weighting functions for various frequencies peak at different levels within a rainy atmosphere. GPROF is one of two TMI rainfall algorithms. It is physically based retrieval that finds the vertical hydrometeor profile that best fits the brightness temperatures in the available passive radiometer channels. Matching is achieved using a library of hydrometeor profiles generated by cloud-resolving models (CRMs). The hydrometeor profiles have a corresponding surface precipitation rate. The algorithm retrieves the hydrometeor profiles and associated surface rainfall using a Bayesian approach that gives the estimated expected values. The ability of CRMs to produce cloud structures that are reliable and representative of observed storms is crucial for the success of GPROF. The cloud mycrophysics are one of the keys to achieving this. In addition, CRMs have been a very useful tool for GPM-algorithm developers through Cloud-Radiation Simulations (CRS), one of the nine GPM disciplinary research themes. This paper will discuss how to generate consistent and comprehensive 4D cloud datasets from an improved (i.e., in regard to bulk and multi-moment microphysics) CRM for TRMM and GPM rainfall retrieval algorithm developers. These cloud datasets include CRM-simulated clouds and cloud systems from different geographic locations in the tropics and midlatitudes. By linking the CRM with a passive microwave radiative-transfer model and using satellite and airborne data, the performance of the "cloud physics" can be assessed and in turn modified and improved. This paper will also address how to assess and improve the performance of various latent and diabatic heating algorithms and develop an algorithm to retrieve the vertical structure of apparent moistening (Q2). Considering that the GPM will produce high (temporal and spatial) resolution heating and rainfall data, these algorithms will be used to obtain the temporal and spatial distributions of surface rainfall and the associated vertical heating and moistening profiles throughout the subtropical and midlatitudes.

Tao, Wei-Kuo↗
NASA NTRS1999

Wave Dynamics and Transport in the Stratosphere

The report discusses: (1) Gravity waves generated by tropical convection: A study in which a two-dimensional cloud-resolving model was used to examine the possible role of gravity waves generated by a simulated tropical squall line in forcing the quasi-biennial oscillation was completed. (2) Gravity wave ray tracing studies:It was developed a linear ray tracing model of gravity wave propagation to extend the nonlinear storm model results into the mesosphere and thermosphere. (3) tracer filamentation: Vertical soundings of stratospheric ozone often exhibit laminated tracer structures characterized by strong vertical tracer gradients. (4) Mesospheric gravity wave modeling studies: Although our emphasis in numerical simulation of gravity waves generated by convection has shifted from simulation of idealized two-dimensional squall lines to the most realistic (and complex) study of wave generation by three-dimensional storms. (5) Gravity wave climatology studies: Mr. Alexander applied a linear gravity wave propagation model together with observations of the background wind and stability fields to compute climatologies of gravity wave activity for comparison to observations. (6) Convective forcing of gravity waves: Theoretical study of gravity wave forcing by convective heat sources has completed. (7) Gravity waves observation from UARS: The objective of this work is to apply ray tracing, and other model technique, in order to determine to what extend the horizontal and vertical variation in satellite observed distribution of small-scale temperature variance can be attributed to gravity waves from particular sources. (8) The annual and interannual variations in temperature and mass flux near the tropical tropopause. and (9) Three dimensional cloud model.

Holton, James R.↗
NASA NTRS2011

Recent THEMIS and Coordinated GBO Measurements of Substorm Expansion Onset: Do We Finally Have an Answer?

For nearly 30 years an often-times heated debate has engaged the substorm community: Do substorms begin with the formation of a new reconnection site in the midtail plasma sheet (the Near-Earth Neutral Line model) or do they begin near the transition region between stretched tail and dipolar field lines (the Current Disruption model). The THEMIS mission, with a coordinated suite of five in-situ spacecraft and ground observatories, has greatly extended our understanding of how substorms initiate and evolve. But have the new data resolved the fundamental question? In this talk I review the last few year's of sub storm research, with an emphasis of how the THEMIS data have revolutionized our understanding.

Kepko, L.↗
NASA NTRS2016

Empirical Modeling of the Storm Time Innermost Magnetosphere Using Van Allen Probes and THEMIS Data: Eastward and Banana Currents

The structure of storm time currents in the inner magnetosphere, including its innermost region inside 4R(sub E), is studied for the first time using a modification of the empirical geomagnetic field model TS07D and new data from Van Allen Probes and Time History of Events and Macroscale Interactions during Substorms missions. It is shown that the model, which uses basis-function expansions instead of ad hoc current modules to approximate the magnetic field, consistently improves its resolution and magnetic field reconstruction with the increase of the number of basis functions and resolves the spatial structure and evolution of the innermost eastward current. This includes a connection between the westward ring current flowing largely at R > or approx. 3R(sub E) and the eastward ring current concentrated at R < or approx. 3R(sub E) resulting in a vortex current pattern. A similar pattern coined 'banana current' was previously inferred from the pressure distributions based on the energetic neutral atom imaging and first-principles ring current simulations. The morphology of the equatorial currents is dependent on storm phase. During the main phase, it is complex, with several asymmetries forming banana currents. Near SYM-H minimum, the banana current is strongest, is localized in the evening-midnight sector, and is more structured compared to the main phase. It then weakens during the recovery phase resulting in the equatorial currents to become mostly azimuthally symmetric.

Stephens, G. K.↗
DOE OSTI2026

Improving Process Level Understanding of Boundary Layer Winds over the Northeast U.S. Shelf: The Third Wind Forecast Improvement Project (WFIP3)

The third Wind Forecast Improvement Project (WFIP3), a U.S. Department of Energy and National Oceanic and Atmospheric Administration sponsored investigation, sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the marine atmospheric boundary layer. WFIP3 focused on mesoscale and submesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high-use coastal zone, using a 3D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the marine atmospheric boundary layer over the ocean was done from an air-sea interaction flux tower and extended deployments of a large autonomous barge platform. Numerous critical forecasting phenomena were observed that are being evaluated within regional coupled and uncoupled modeling systems, including the National Oceanic and Atmospheric Administration's foundational High-Resolution Rapid Refresh forecast model.

Kirincich, Anthony↗
NASA NTRS2000

The Estimation of Hydrometeor Profiles from Wideband Microwave Observations

Profiles of the microphysical properties of clouds and raincells are essential in many areas of atmospheric research and operational meteorology. In order to enhance the understanding of the nonlinear and underconstrained relationships between cloud and hydrometeor microphysical profiles and passive microwave brightness temperatures, estimations of cloud profiles for an anvil, a convective, and an updraft region of an oceanic squall were performed. The estimations relied on comparisons between radiative transfer calculations of incrementally estimated microphysical profiles and concurrent dual-altitude wideband brightness temperatures from the 22 February 1993 flight during TOGA-COARE. The wideband observations (10--220 GHz) are necessary for estimating cloud profiles reaching up to 20 km. The low frequencies enhance the rain and cloud water profiles, while the high frequencies are required to detail the higher altitude ice microphysics. A microphysical profile was estimated for each of the three regions of the storm. Each of the three estimated profiles produced calculated brightness temperatures within approximately 10 K of the observations. A majority, of the total iterative adjustment were to the estimated profile's frozen hydrometeor characteristics and were necessary to match the high frequency calculations with the observations. This indicates a need to validate cloud resolving models using high frequencies. Some difficulties matching the 37 GHz observation channels on the DC-8 and ER-2 aircrafts with the calculations simulated at the two aircraft heights (approximately 11 km and 20 km, respectively) were noted and potential causes presented.

Gail M Skofronick-Jackson↗
NASA NTRS2020

An Examination of the Impact of Grid Spacing on WRF Simulations of Wintertime Precipitation in the Mid-Atlantic United States

A large set of deterministic and ensemble forecasts was produced to identify the optimal spacing for forecasting U.S. East Coast snowstorms. WRF forecasts were produced on cloud-allowing (~1-km grid spacing) and convection-allowing (3–4 km) grids, and compared against forecasts with parameterized convection (>~10 km). Performance diagrams were used to evaluate 19 deterministic forecasts from the winter of 2013–14. Ensemble forecasts of five disruptive snowstorms spanning the years 2015–18 were evaluated using various methods to evaluate probabilistic forecasts. While deterministic forecasts using cloud-allowing grids were not better than convection-allowing forecasts, both had lower bias and higher success ratios than forecasts with parameterized convection. All forecasts were underdispersive. Nevertheless, forecasts on the higher-resolution grids were more reliable than those with parameterized convection. Forecasts on the cloud-allowing grid were best able to discriminate areas that received heavy snow and those that did not, while the forecasts with parameterized convection were least able to do so. It is recommended to use convection-resolving and (if computationally possible) to use cloud-allowing forecast grids when predicting East Coast winter storms.

Snowstorms↗
NASA NTRS1984

Clearing the Martian air - The troubled history of dust storms

This note is an attempt to resolve some misconceptions regarding the historical record of the Martian atmospheric phenomena referred to as 'dust storms,' but often called yellow storms, yellow clouds, planetwide dust storms, global dust storms, great dust storms, etc. The known frequency of planet-encircling storms will be specifically addressed. Better knowledge of the sizes, frequencies, and locations of Martian dust storms is needed for atmospheric modeling and for future mission planning.

Martin, L. J.↗
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.↗
NASA NTRS2005

Lightning NO(x) Production and Its Consequences for Tropospheric Chemistry

Cloud-resolving case-study simulations of convective transport and lightning NO production have yielded results which are directly applicable to the design of lightning parameterizations for global chemical transport models. In this work we have used cloud-resolving models (the Goddard Cumulus Ensemble Model (GCE) and MMS) to drive an off-line cloud-scale chemical transport model (CSCTM). The CSCTM, in conjunction with aircraft measurements of NO(x) in thunderstorms and ground-based lightning observations, has been used to constrain the amount of NO produced per flash. Cloud and chemistry simulations for several case studies of storms in different environments will be presented. Observed lightning flash rates have been incorporated into the CSCTM, and several scenarios of NO production per intracloud (IC) and per cloud-to-ground (CG) flash have been tested for each storm. The resulting NO(x) mixing ratios are compared with aircraft measurements taken within the storm (typically the anvil region) to determine the most likely NO production scenario. The range of values of NO production per flash (or per meter of lightning channel length) that have been deduced from the model will be shown and compared with values of production in the literature that have been deduced from observed NO spikes and from anvil flux calculations. Results show that on a per flash basis, IC flashes are nearly as productive of NO as CG flashes. This result simplifies the lightning parameterization for global models (ie., an algorithm for estimating the ICKG ratio is not necessary). Vertical profiles of lightning NOx mass at the end of the 3-D storm simulations have been summarized to yield suggested profiles for use in global models. Estimates of mean NO production per flash vary by a factor of three from one simulated storm to another. When combined with the global flash rate of 44 flashes per second from NASA's Optical Transient Detector (OTD) measurements, these estimates and the results from other techniques yield global NO production rates of 2-9 TgN/year. Simulations of the photochemistry over the 24 hours following a storm has been performed to determine the additional ozone production which can be attributed to lightning NO. Convective transport of HO(x), precursors leads to the generation of a HO(x), plume which substantially aids the downstream ozone production.

Pickering, Kenneth E.↗