Search

MetPy is an open-source, Python-based package for meteorology, providing domain-specific functionality built extensively on top of the robust scientific Python software stack, which includes libraries like NumPy, SciPy, Matplotlib, and xarray. The goal of the project is to bring the weather analysis capabilities of GEMPAK (and similar software tools) into a modern computing paradigm. MetPy strives to employ best practices in its development, including software tests, continuous integration, and automated publishing of web-based documentation. As such, MetPy represents a sustainable, long-term project that fills a need for the meteorological community. MetPy’s development is substantially driven by its user community, both through feedback on a variety of open, public forums like Stack Overflow, and through code contributions facilitated by the GitHub collaborative software development platform. MetPy has recently seen the release of version 1.0, with robust functionality for analyzing and visualizing meteorological datasets. While previous versions of MetPy have already seen extensive use, the 1.0 release represents a significant milestone in terms of completeness and a commitment to long-term support for the programming interfaces. This article provides an overview of MetPy’s suite of capabilities, including its use of labeled arrays and physical unit information as its core data model, unit-aware calculations, cross sections, skew T and GEMPAK-like plotting, station model plots, and support for parsing a variety of meteorological data formats. The general road map for future planned development for MetPy is also discussed.

Lightning is frequently initiated within the convective regions of thunderstorms, and so flash rates tend to follow trends in updraft speed and volume. It has been suggested that lightning production is linked to the turbulent flow generated by updrafts as turbulent eddies organize charged hydrometeors into complex charge structures. These complex charge structures consist of local regions of increased charge magnitudes between which flash-initiating electric fields may be generated. How turbulent kinematics influences lightning production, however, remains unclear. In this study, lightning flashes produced in a multicell and two supercell storms simulated using The Collaborative Model for Multiscale Atmospheric Simulation (COMMAS) were examined to identify the kinematic flow structures within which they occurred. By relating the structures of updrafts to thermals, initiated lightning flashes were expected to be located where the rate of strain and rotational flow are equal, or between updraft and eddy flow features. Results showed that the average lightning flash is initiated in kinematic flow structures dominated by vortical flow patterns, similar to those of thermals, and the structures’ kinematics are characterized by horizontal vorticity and vertical shearing. These kinematic features were common across all cases and demonstrated that where flash-initiating electric fields are generated is along the periphery of updrafts where turbulent eddies are produced. Careful consideration of flow structures near initiated flashes is consistent with those of thermals rising through a storm.

This study employed a parallel-plate capacitor model by which the electrostatic energy of lightning flashes could be estimated by considering only their physical dimensions and breakdown electric fields in two simulated storms. The capacitor model has previously been used to approximate total storm electrostatic energy but is modified here to use the geometry of individual lightning flashes to mimic the local charge configuration where flashes were initiated. The energy discharged may then be diagnosed without context of a storm’s entire charge structure. The capacitor model was evaluated using simulated flashes from two storms modeled by the National Severe Storms Laboratory’s Collaborative Model for Multiscale Atmospheric Simulation (COMMAS). Initial capacitor model estimates followed the temporal evolution of the flash discharge energy of COMMAS for each storm but demonstrated the need to account for an adjustment factor μc to represent the fraction of energy a flash dissipates, as this model assumes the entire preflash energy is discharged by a flash. Individual values of μc were obtained simply by using the ratio of the COMMAS flash to capacitor energy. Median values were selected to represent the flash populations for each storm, and were in range of . Application of aligned the magnitudes of the capacitor model discharge energy estimates to those of COMMAS and to those estimated in previous studies. Therefore, by considering a μc within range of , application of the capacitor model for observed lightning datasets is suggested.

A lightning risk assessment for application to human safety was created and applied in 10 west Texas locations from 2 May 2016 to 30 September 2016. The method combined spatial lightning mapping data, probabilistic risk calculation adapted from the International Electrotechnical Commission Standard 62305-2, and weighted average interpolation to produce risk magnitudes that were compared with tolerability thresholds to issue lightning warnings. These warnings were compared with warnings created for the same dataset using a more standard lightning safety approach that was based on National Lightning Detection Network (NLDN) total lightning within 5 n mi (1 n mi = 1.852 km) of each location. Four variations of the calculation as well as different units of risk were tested to find the optimal configuration to calculate risk to an isolated human outdoors. The best-performing risk configuration using risk (10 min)−1 or larger produced the most comparable results to the standard method, such as number of failures, average warning duration, and total time under warnings. This risk configuration produced fewer failures than the standard method but longer total time under warnings and higher false alarm ratios. Median lead times associated with the risk configuration were longer than the standard method for all units considered, whereas median down times were shorter for risk (10 min)−1 and risk (15 min)−1. Overall, the risk method provides a baseline framework to quantify the changing lightning hazard on the storm scale and could be a useful tool to aid in lightning decision support scenarios.

We demonstrate the utility of transient polarimetric signatures (ZDR and KDP columns, a proxy for surges in a thunderstorm updraft) to explain variability in lightning flash rates in a tornadic supercell. Observational data from a WSR-88D and the Oklahoma lightning mapping array are used to map the temporal variance of polarimetric signatures and VHF sources from lightning channels. It is shown, via three-dimensional and cross-sectional analyses, that the storm was of inverted polarity resulting from anomalous electrification. Statistical analysis confirms that mean flash area in the ZDR column region was 10 times smaller than elsewhere in the storm. On an average, 5 times more flash initiations occurred within ZDR column regions, thereby supporting existing theory of an inverse relationship between flash initiation rates and lightning channel extent. Segmentation and object identification algorithms are applied to gridded radar data to calculate metrics such as height, width, and volume of ZDR and KDP columns. Variability in lightning flash rates is best explained by the fluctuations in ZDR column volume with a Spearman’s rank correlation coefficient value of 0.72. The highest flash rates occur in conjunction with the deepest ZDR columns (up to 5 km above environmental melting level) and largest volumes of ZDR columns extending up to the −20°C level (3 km above the melting level). Reduced flash rates toward the end of the analysis are indicative of weaker updrafts manifested as low ZDR column volumes at and above the −10°C level. These findings are consistent with recent studies linking lightning to the interplay between storm dynamics, kinematics, thermodynamics, and precipitation microphysics.

The existence of mesoscale lightning discharges on the order of 100 km in length has been known since the radar-based findings of Ligda in the mid-1950s. However, it took the discovery of sprites in 1989 to direct significant attention to horizontally extensive “megaflashes” within mesoscale convective systems (MCSs). More recently, 3D Lightning Mapping Arrays (LMAs) have documented sprite-initiating lightning discharges traversing several hundred kilometers. One such event in a 2007 Oklahoma MCS having an LMA-derived length of 321 km, has been certified by the WMO as the longest officially documented lightning flash. The new Geostationary Lightning Mapper (GLM) sensor on GOES-16/17 now provides an additional tool suited to investigating mesoscale lightning. On 22 October 2017, a quasi-linear convective system moved through the central United States. At 0513 UTC, the GLM indicated a lightning discharge originated in northern Texas, propagated north-northeast across Oklahoma, fortuitously traversed the Oklahoma LMA (OKLMA), and finally terminated in southeastern Kansas. This event is explored using the OKLMA, the National Lightning Detection Network (NLDN), and the GLM. The NLDN reported 17 positive cloud-to-ground flashes (+CGs), 23 negative CGs (−CGs), and 37 intracloud flashes (ICs) associated with this massive discharge, including two +CGs capable of inducing sprites, with others triggering upward lightning from tall towers. Combining all available data confirms the megaflash, which illuminated 67,845 km2, was at least 500 km long, greatly exceeding the current official record flash length. Yet even these values are being superseded as GLM data are further explored, revealing that such vast discharges may not be all that uncommon.

Cleaning in a flash Hydroxyl radicals (OH) are the most important oxidizing species in the atmosphere and provide much of its ability to cleanse itself. It is known that nitric oxide production by lightning leads to the formation of OH and other atmospheric oxidants, such as ozone and hydroperoxyl radicals (HO 2 ), through a variety of chemical reactions. Brune et al. used airborne measurements of OH and HO 2 to show that lightning also produces them directly and in amounts much greater than expected. They found that this mechanism may be responsible for as much as one-sixth of the oxidizing capacity of Earth's atmosphere. Science , this issue p. 711

Large-eddy-resolving simulations using the Collaborative Model for Multiscale Atmospheric Simulation (COMMAS), which contains microphysical charging and branched-lightning parameterizations, produce much more complex net charge structures than conventionally visualized from previous observations, simulations, and conceptual diagrams. Many processes contribute to the hydrometeor charge budget within a thunderstorm, including advection, hydrometeor differential sedimentation, subgrid turbulent mixing and diffusion, ion drift, microphysical separation, and the attachment of ion charge deposited by the lightning channel. The lightning deposition, sedimentation, and noninductive charging tendencies contribute the most overall charge at relatively large scales, while the advection tendency, from resolved turbulence, provides the most “texture” at small scales to the net charge density near the updraft region of the storm. The scale separation increases for stronger storm simulations. In aggregate, lightning deposition and sedimentation resemble the smoother distribution of the electric potential, while evidence suggests individual flashes could be responding to the fine texture in the net charge. The clear scale separation between the advection and other net charge tendencies suggest the charge advection is most capable of providing net charge texture; however, a clear-cut causality is not obtained from this study.

This study examined 34 lightning flashes within four separate thundersnow events derived from lightning mapping arrays (LMAs) in northern Alabama, central Oklahoma, and Washington DC. The goals were to characterize the in-cloud component of each lightning flash, as well as the correspondence between the LMA observations and lightning data taken from national lightning networks like the National Lightning Detection Network (NLDN). Individual flashes were examined in detail to highlight several observations within the dataset. The study results demonstrated that the structures of these flashes were primarily normal polarity. The mean area encompassed by this set of flashes is 375 km(2), with a maximum flash extent of 2300 km(2), a minimum of 3 km(2), and a median of 128 km(2). An average of 2.29 NLDN flashes were recorded per LMA-derived lightning flash. A maximum of 11 NLDN flashes were recorded in association with a single LMA-derived flash on 10 January 2011. Additionally, seven of the 34 flashes in the study contain zero NLDN identified flashes. Eleven of the 34 flashes initiated from tall human-made objects (e.g., communication towers). In at least six lightning flashes, the NLDN detected a return stroke from the cloud back to the tower and not the initial upward leader. This study also discusses lightning’s interaction with the human built environment and provides an example of lightning within heavy snowfall observed by GOES-16’s Geostationary Lightning Mapper.

High-resolution wind measurements at 2.25 m in height were used to investigate the mean and turbulence properties of an extreme thunderstorm wind event in West Texas. These data were combined with single Doppler scans from the Texas Tech University Ka-band mobile Doppler radars systems (TTUKa) to provide meteorological context over the surface measurement stations for portions of the outflow. Several features characteristic of a severe wind event were noted in the radar data, including a bowing portion of the thunderstorm complex and a small circulation on the leading edge. These features were reflected in the surface wind time histories and provided natural separation between various regions of the outflow. These features also contributed to the peak 1-s gust at all measurement stations. The turbulence characteristics of each outflow region were also investigated and compared. Reduced values of running turbulence intensity and elevated values of longitudinal integral scales were noted during the period of peak wind speed. Larger scales of turbulence within the outflow were also suggested via spectral analysis.

Convective clouds play an important role in the Earth's climate system as a driver of large-scale circulations and a primary mechanism for the transport of heat, moisture, aerosols, and momentum throughout the troposphere. Despite their climatic importance, multi-scale models continue to have persistent biases produced by insufficient representation of convective clouds. This is the result of an incomplete understanding of key processes such as convective initiation, updraft and downdraft dynamics, cloud and precipitation microphysics, and aerosol-convection interactions.The Aerosol-Cloud-Precipitation-Climate Initiative, an international research group dedicated to advancing understanding of aerosol impacts on clouds relevant to climate, has identified the Houston, Texas region as an optimal location for targeted studies of aerosol-convection interactions within frequently developing isolated deep convection. Houston lies within a humid subtropical climate regime, where onshore flow and sea-breeze convection interact with a range of aerosol conditions associated with Houston's urban and industrial emissions. Pilot studies have suggested that convective clouds in this region are potentially significantly impacted by the varying aerosol conditions.

For over a century, the electric Potential Gradient (PG) of the atmosphere has been measured and studied. The local vertical electric field (Ez) is strongly influenced by the presence of lightning, electrified clouds, rainfall, aerosols, and many others. The One Year Electric Field Study-North Slope of Alaska (OYES-NSA) field campaign was established in the summer of 2017 to measure the vertical electric field at the ARM site in Barrow, Alaska alongside a wide array of supplementary instrumentation, including a Micro-Pulse Lidar and upward facing Ka-band radar. Two years of observations (072017-062019) have shown the possibility to quantify the local effects from aerosols and clouds observed by the Lidar and Radar on the measured EZ. Throughout the manuscript, the physics convention of negative downward fair-weather electric fields is used. Three cases (convective clouds, high concentration of near surface aerosols, and blowing snow) are used to demonstrate the localized effects on the measured EZ. Utilizing the relationships between EZ and backscatter/reflectivity, we have developed a methodology to distinguish samples with local influences. A fair-weather (FW) condition is determined to be associated with a low Lidar backscattering (less than 15 km−1sr−1), in the presence of no significant cloud activity (radar reflectivity less than −10dBZ). The samples satisfying these criteria are found with a 5-min averaged standard deviation of less than 15 V/m, and EZ between −250 V/m to −50 V/m. Using only properties of the EZ measurements allows for the simultaneous comparisons of FW at multiple sites, without the need for supplementary information of local weather conditions. Simultaneous EZ measurements from 8 FW cases are shown between Barrow, AK and Corpus Christi, TX on the timescale of minutes to hours. Similar variation patterns in the FW EZ are shown at both sites, providing evidence of the global nature of the atmospheric electric system. Furthermore, the seasonal-diurnal variability of FW at multiple sites shows similar distributions of the PG.

Approximately 63 million lightning flashes have been identified and analyzed from multiple years of Washington, D. C., northern Alabama, and northeast Colorado lightning mapping array (LMA) data using an open-source flash-clustering algorithm. LMA networks detect radiation produced by lightning breakdown processes, allowing for high-resolution mapping of lightning flashes. Similar to other existing clustering algorithms, the algorithm described herein groups lightning-produced radiation sources by space and time to estimate total flash counts and information about each detected flash. Various flash characteristics and their sensitivity to detection efficiency are investigated to elucidate biases in the algorithm, detail detection efficiencies of various LMAs, and guide future improvements. Furthermore, flash density values in each region are compared to corresponding satellite estimates. While total flash density values produced by the algorithm in Washington, D. C. (~20 flashes km−2 yr−1), and Alabama (~35 flashes km−2 yr−1) are within 50% of satellite estimates, LMA-based estimates are approximately a factor of 3 larger (50 flashes km−2 yr−1) than satellite estimates in northeast Colorado. Accordingly, estimates of the ratio of in-cloud to cloud-to-ground flashes near the LMA network (~20) are approximately a factor of 3 larger than satellite estimates in Colorado. These large differences between estimates may be related to the distinct environment conducive to intense convection, low-altitude flashes, and unique charge structures in northeast Colorado.

A database consisting of approximately 4000 storm observations has been objectively analyzed to determine environmental characteristics that produce high radar reflectivities above the freezing level, large total lightning flash rates on the order of 10 flashes per minute, and anomalous vertical charge structures (most notably, dominant midlevel positive charge). The storm database is drawn from four regions of the United States featuring distinct environments, each with coinciding Lightning Mapping Array (LMA) network data. LMAs are able to infer total lightning flash rates using flash clustering algorithms, such as the one implemented in this study. Results show that anomalous charge structures inferred from LMA data, significant lightning flash rates, and increased radar reflectivities above the freezing level tend to be associated with environments that have high cloud base heights (approximately 3 km above ground level) and large atmospheric instability, quantified by normalized convective available potential energy (NCAPE) near 0.2 m s−2. Additionally, we infer that aerosols may affect storm intensity. Maximum flash rates were observed in storms with attributed aerosol concentrations near 1000 cm−3, while total flash rates decrease when aerosol concentrations exceed 1500 cm−3, consistent with previous studies. However, this effect is more pronounced in regions where the NCAPE and cloud base height are low. The dearth of storms with estimated aerosol concentrations less than 700 cm−3 (approximately 1% of total sample) does not provide a complete depiction of aerosol invigoration.

The 2009 Redoubt Volcano eruption produced a series of explosive events, the largest of which produced episodes of volcanic lightning similar to thunderstorms. Flash energy spectra were calculated from lightning mapping data collected during the eruption. The spectra were compared to the turbulence characteristics expected from each stage of plume development. Small flash length scales present at early times were associated with the gas thrust and initial convective stages. Increases in flash length scales and flash energy as the explosive events progressed were associated with an increase in volume of the plume. Spectra with a large range in flash length scales (0.2–10 km) and a spectral peak at small flash length scales (< 1 km) were associated with a superposition of spectra from gas thrust, convective, and umbrella/ash cloud regions. An approximate 5/3 power law slope was observed on the order of 10 min after gas thrust forcing had ended and the transition to a drifting ash cloud was underway, which may reflect the dissipating state of the drifting ash cloud.

Several field campaigns since the year 2000 have focused on anomalously electrified or “inverted polarity” thunderstorms. This study synthesizes these recent results, and considers how variability in the non-inductive relative-growth rate electrification mechanism might clarifying the meaning of “inverted polarity”. Instead of falling into two polarity classes, electrification and charge structure in strong updrafts vary continuously, as expected if depletion of supercooled water is a primary control on electrification. Two- or three-dimensional storm flows or other electrification mechanisms are required to combine one or more of these electrification regimes into “inverted” or otherwise complicated local charge sequences. Cloud flashes that result from these local charge sequences should be termed “positive” and “negative” instead of “normal” and “inverted” because cloud flashes of either polarity can occur at any altitude in thunderstorms.

Previous analyses of very high frequency (VHF) Lightning Mapping Array (LMA) observations relative to the location of deep convective updrafts have noted a systematic pattern in flash characteristics. In and near strong updrafts, flashes tend to be smaller and more frequent, while flashes far from strong vertical drafts exhibit the opposite tendency. This study quantitatively tests these past anecdotal observations using LMA data for two supercell storms that occurred in Oklahoma in 2004. The data support a prediction from electrostatics that frequent breakdown and large flash extents are opposed. An energetic scaling that combines flash rate and flash area exhibits a power-law scaling regime on scales of a few kilometers and a maximum in flash energy at about 10 km. The spectral shape is surprisingly consistent across a range of moderate to large flash rates. The shape of this lightning flash energy spectrum is similar to that expected of turbulent kinetic energy spectra in thunderstorms. In line with the hypothesized role of convective motions as the generator of thunderstorm electrical energy, the correspondence between kinematic and electrical energy spectra suggests that advection of charge-bearing precipitation by the storm’s flow, including in turbulent eddies, couples the electrical and kinematic properties of a thunderstorm.

[1] The temporal and spatial development of sprite-producing lightning flashes is examined with coordinated observations over an asymmetric mesoscale convective system (MCS) on 29 June 2011 near the Oklahoma Lightning Mapping Array (LMA). Sprites produced by a total of 26 lightning flashes were observed simultaneously on video from Bennett, Colorado and Hawley, Texas, enabling a triangulation of sprites in comparison with temporal development of parent lightning (in particular, negatively charged stepped leaders) in three-dimensional space. In general, prompt sprites produced within 20 ms after the causative stroke are less horizontally displaced (typically 30 km). However, both prompt and delayed sprites are usually centered within 30 km of the geometric center of relevant LMA sources (with affinity to negative stepped leaders) during the prior 100 ms interval. Multiple sprites appearing as dancing/jumping events associated with a single lightning flash could be produced either by distinct strokes of the flash, by a single stroke through a series of current surges superposed on an intense continuing current, or by both. Our observations imply that sprites elongated in one direction are sometimes linked to in-cloud leader structure with the same elongation, and sprites that were more symmetric were produced above the progression of multiple negative leaders. This suggests that the large-scale structure of sprites could be affected by the in-cloud geometry of positive charge removal. Based on an expanded dataset of 39 sprite-parent flashes by including more sprites recorded by one single camera over the same MCS, the altitude (above mean sea level, MSL) of positively charged cloud region tapped by sprite-producing strokes declined gradually from ~10 km MSL (−35°C) to around 6 km MSL (−10°C) as the MCS evolved through the mature stage. On average, the positive charge removal by causative strokes of sprites observed on 29 June is centered at 3.6 km above the freezing level or at 7.9 km above ground level.

This study characterizes Lightning Mapping Array performance for networks that participated in the Deep Convective Clouds and Chemistry field program using new Monte Carlo and curvature matrix model simulations. These open‐source simulation tools are readily adapted to real‐time operations or detailed studies of performance. Each simulation accounted for receiver threshold and location, as well as a reference distribution of source powers and flash sizes based on thunderstorm observations and the mechanics of station triggering. Source and flash detection efficiency were combined with solution bias and variability to predict flash area distortion at long ranges. Location errors and detection efficiency were highly dependent on the station configuration and thresholds, especially at longer ranges, such that performance varied more than expected across different networks and with azimuth within networks. Error characteristics matched prior studies, which led to an increase in flash distortion with range. Predicted flash detection efficiency exceeded 95% within 100 km of all networks., Key Points LMA location errors and detection efficiency are highly dependent on the station configuration and thresholds, especially at longer rangesPerformance varied greatly across different DC3 networks and with azimuth in each, but overall error characteristics matched prior studiesPredicted flash detection efficiency exceeded 95% within 100 km of all DC3 networks, and these flashes are distorted more at larger ranges

The reported range for global production of nitrogen oxides (NOx = NO + NO2) by lightning remains large (e.g., 32 to 664 mol NOx flash−1), despite incorporating results from over 30 individual laboratory, theoretical, and field studies since the 1970s. Airborne and ground-based observations from the Deep Convective Clouds and Chemistry experiment in May and June 2012 provide a new data set for calculating moles of NOx produced per lightning flash, P(NOx), in thunderstorms over the United States Great Plains. This analysis utilizes a combination of in situ observations of storm inflow and outflow from three instrumented aircraft, three-dimensional spatial information from ground-based radars and satellite observations, and spatial and temporal information for intracloud and cloud-to-ground lightning flashes from ground-based lightning mapping arrays. Evaluation of two analysis methods (e.g., a volume-based approach and a flux-based approach) for converting enhancements in lightning-produced NOx from volume-based mixing ratios to moles NOx flash−1 suggests that both methods equally approximate P(NOx) for storms with elongated anvils, while the volume-based approach better approximates P(NOx) for storms with circular-shaped anvils. Results from the more robust volume-based approach for three storms sampled over Oklahoma and Colorado during DC3 suggest a range of 142 to 291 (average of 194) moles NOx flash−1 (or 117–332 mol NOx flash−1 including uncertainties). Although not vastly different from the previously reported range for storms occurring in the Great Plains (e.g., 21–465 mol NOx flash−1), results from this analysis of DC3 storms offer more constrained upper and lower limits for P(NOx) in this geographical region.

This study examined lightning activity relative to the rapidly evolving kinematics of a hail-producing storm on 15 August 2006. Data were provided by the National Weather Radar Testbed Phased-Array Radar, the Oklahoma Lightning Mapping Array, and the National Lightning Detection Network. This analysis is the first to compare the electrical characteristics of a hail-producing storm with the reflectivity and radial velocity structure at temporal resolutions of less than 1 min. Total flash rates increased to approximately 220 min−1 as the storm’s updraft first intensified, leveled off during its first mature stage, and then decreased for 2–3 min despite the simultaneous development of another updraft surge. This reduction in flash rate occurred as wet hail formed in the new updraft and was likely related to the wet growth; wet growth is not conducive to hydrometeor charging and probably contributed to the formation of a “lightning hole” without a mesocyclone. Total flash rates subsequently increased to approximately 450 min−1 as storm volume and inferred graupel volume increased, and then decreased as the storm dissipated. The vertical charge structure in the storm initially formed a positive tripole (midlevel negative charge between upper and lower positive charges). The charge structure in the second updraft surge consisted of a negative charge above a deep midlevel positive charge, a reversal consistent with the effects of large liquid water contents on hydrometeor charge polarity in laboratory experiments. Prior to the second updraft surge, the storm produced two cloud-to-ground flashes, both lowering the usual negative charge to ground. Shortly before hail likely reached ground, the storm produced four cloud-to-ground flashes, all lowering the positive charge. Episodes of high singlet VHF sources were observed at approximately 13–15 km during the initial formation and later intensification of the storm’s updraft.

Lightning mapping, electric field, and radar data from the 26 May 2004 supercell in central Oklahoma are used to examine the storm’s charge structure. An initial arc-shaped maximum in lightning activity on the right flank of the storm’s bounded weak echo region was composed of an elevated normal polarity tripole in the region of precipitation lofted above the storm’s weak echo region. Later in the storm, two charge structures were associated with precipitation that reached the ground. To the left of the weak echo region, six charge regions were inferred, with positive charge carried on hail at the bottom of the stack. Farther forward in the storm’s precipitation region, four charge regions were inferred, with negative charge at the bottom of the stack. There were different charge structures in adjacent regions of the storm at the same time, and regions of opposite polarity charge were horizontally adjacent at the same altitude. Flashes occasionally lowered positive charge to ground from the forward charge region. A conceptual model is presented that ties charge structure in different regions of the storm to storm structure inferred from radar reflectivity.

Electrification and lightning are simulated for a small continental multicell storm. The results are consistent with observations and thus provide additional understanding of the charging processes and evolution of this storm. The first six observed lightning flashes were all negative cloud-to-ground (CG) flashes, after which intracloud (IC) flashes also occurred between middle and upper levels of the storm. The model simulation reproduces the basic evolution of lightning from low and middle levels to upper levels. The observed lightning indicated an initial charge structure of at least an inverted dipole (negative charge above positive). The simulations show that noninductive charge separation higher in the storm can enhance the main negative charge sufficiently to produce negative CG flashes before upper-level IC flashes commence. The result is a “bottom-heavy” tripole charge structure with midlevel negative charge and a lower positive charge region that is more significant than the upper positive region, in contrast to the traditional tripole structure that has a less significant lower positive charge region. Additionally, the occurrence of cloud-to-ground lightning is not necessarily a result of excess net charge carried by the storm, but it is primarily caused by the local potential imbalance between the lowest charge regions. The two-moment microphysics scheme used for this study predicted mass mixing ratio and number concentration of cloud droplets, rain, ice crystals, snow, and graupel. Bulk particle density of graupel was also predicted, which allows a single category to represent a greater range of particle characteristics. (An additional hail category is available but was not needed for the present study.) The prediction of hydrometeor number concentration is particularly critical for charge separation at higher temperatures (−5° < T < −20°C) in the mixed phase region, where ice crystals are produced by rime fracturing (Hallett–Mossop process) and by splintering of freezing drops. Cloud droplet concentration prediction also affected the rates of inductive charge separation between graupel and droplets.

Measurements during TELEX by a lightning mapping array, polarimetric and mobile Doppler radars, and balloon-borne electric-field meters and radiosondes show how lightning and other electrical properties depend on storm structure, updrafts, and precipitation formation.

Data from a three-dimensional lightning mapping array (LMA) and from two soundings by balloon-borne electric field meters (EFMs) were used to analyze the electrical structures of a multicell storm observed on 25 June 2000 during the Severe Thunderstorm Electrification and Precipitation Study (STEPS). This storm had a complex, multicell structure with four sections, each of whose electrical structure differed from that of the others during all or part of the analyzed period. The number of vertically stacked charge regions in any given section of the storm ranged from two to six. The most complex charge and lightning structures occurred in regions with the highest reflectivity values and the deepest reflectivity cores. Intracloud flashes tended to concentrate in areas with large radar reflectivity values, though some propagated across more than one core of high reflectivity or into the low-reflectivity anvil. Intracloud lightning flash rates decreased as the vertical extent and maximum value of reflectivity cores decreased. Cloud-to-ground flash rates increased as cores of high reflectivity descended to low altitudes. Most cloud-to-ground flashes were positive. All observed positive ground flashes initiated between the lowest-altitude negative charge region and a positive charge region just above it. The storm’s complexity makes it hard to classify the vertical polarity of its overall charge structure, but most of the storm had a different vertical polarity than what is typically observed outside the Great Plains. The electrical structure can vary considerably from storm to storm, or even within the same storm, as in the present case.

On 28–29 June 2004 a multicellular thunderstorm west of Oklahoma City, Oklahoma, was probed as part of the Thunderstorm Electrification and Lightning Experiment field program. This study makes use of radar observations from the Norman, Oklahoma, polarimetric Weather Surveillance Radar-1988 Doppler, three-dimensional lightning mapping data from the Oklahoma Lightning Mapping Array (LMA), and balloon-borne vector electric field meter (EFM) measurements. The storm had a low flash rate (30 flashes in 40 min). Four charge regions were inferred from a combination of LMA and EFM data. Lower positive charge near 4 km and midlevel negative charge from 4.5 to 6 km MSL (from 0° to −6.5°C) were generated in and adjacent to a vigorous updraft pulse. Further midlevel negative charge from 4.5 to 6 km MSL and upper positive charge from 6 to 8 km (from −6.5° to −19°C) were generated later in quantity sufficient to initiate lightning as the updraft decayed. A negative screening layer was present near the storm top (8.5 km MSL, −25°C). Initial lightning flashes were between lower positive and midlevel negative charge and started occurring shortly after a cell began lofting hydrometeors into the mixed phase region, where graupel was formed. A leader from the storm’s first flash avoided a region where polarimetric radar suggested wet growth and the resultant absence of noninductive charging of those hydrometeors. Initiation locations of later flashes that propagated into the upper positive charge tracked the descending location of a polarimetric signature of graupel. As the storm decayed, electric fields greater than 160 kV m−1 exceeded the minimum threshold for lightning initiation suggested by the hypothesized runaway breakdown process at 5.5 km MSL, but lightning did not occur. The small spatial extent (≈100 m) of the large electric field may not have been sufficient to allow runaway breakdown to fully develop and initiate lightning.

Balloon soundings were made through two supercell storms during the Severe Thunderstorm Electrification and Precipitation Study (STEPS) in summer 2000. Instruments measured the vector electric field, temperature, pressure, relative humidity, and balloon location. For the first time, soundings penetrated both the strong updraft and the rainy downdraft region of the same supercell storm. In both storms, the strong updraft had fewer vertically separated charge regions than found near the rainy downdraft, and the updraft’s lowest charge was elevated higher, its bottom being near the 40-dBZ boundary of the weak-echo vault. The simpler, elevated charge structure is consistent with the noninductive graupel–ice mechanism dominating charge generation in updrafts. In the weak-echo vault, the amount of frozen precipitation and the time for particle interactions are too small for significant charging. Inductive charging mechanisms and lightning may contribute to the additional charge regions found at lower altitudes outside the updraft. Lightning mapping showed that the in-cloud channels of a positive ground flash could be in any one of the three vertically separated positive charge regions found outside the updraft, but were in the middle region, at 6–8 km MSL, for most positive ground flashes. The observations are consistent with the electrical structure of these storms having been inverted in polarity from that of most storms elsewhere. It is hypothesized that the observed inverted-polarity cloud flashes and positive ground flashes were caused by inverted-polarity storm structure, possibly due to a larger than usual rime accretion rate for graupel in a strong updraft.

Balloon-borne electric field soundings and lightning mapping data have been analyzed for three of the storms that occurred in the Severe Thunderstorm Electrification and Precipitation Study field program in 2000 to determine if the storms had inverted-polarity electrical structures. The polarities of all or some of the vertically stacked charge regions in such storms are opposite to the polarities observed at comparable heights in normal storms. Analyses compared the charge structures inferred from electric field soundings in the storms with charges inferred from three-dimensional lightning mapping data. Charge structures were inferred from electric field profiles by combining the one-dimensional approximation of Gauss's law with additional information from three-dimensional patterns in the electric field vectors. The three different ways of inferring the charge structure in the storms were found to complement each other and to be consistent overall. Charge deposition by lightning possibly occurred and increased the charge complexity of one of the storms. Many of the cloud flashes in each case were inverted-polarity flashes. Two storms produced ground flash activity comprised predominantly of positive ground flashes. One storm, which was an isolated thunderstorm, produced inverted-polarity cloud flashes, but no flashes to ground. The positive and negative thunderstorm charge regions were found at altitudes where, respectively, negative and positive charge would be found in normal-polarity storms. Thus, we conclude that these storms had anomalous and inverted-polarity electrical structures. Collectively, these three cases (along with the limited cases in the refereed literature) provide additional evidence that thunderstorms can have inverted-polarity electrical structures.

[1] From a dataset of about 250 soundings of electric field (E), nine were adversely affected by lightning. These soundings are interpreted as ending near lightning initiation locations. Scaled to standard pressure, the largest observed E was 626 kV m−1 and the largest estimated E was 929 kV m−1. E exceeded runaway breakdown threshold, RBth, by factors of 1.1–3.3 before each flash, and overvoltages were 1.4–4.3. Seven cases had rapid E increases (rates of 11–100 kV m−1 s−1) in the few seconds before the flash, and in three the maximum E occurred 3 s or more before the flash. A tenth sounding with E > RBth for 38 s had subsequent lightning initiate 2 km from the balloon; one channel came within 400 m, but the flash and large E did not adversely affect the instruments. The findings suggest that E > RBth is a necessary condition for lightning initiation, but it is not sufficient.