ALBERT

Description: A taxonomy of tropical convective and stratiform vertical structures is constructed through cluster analysis of 3 yr of Tropical Rainfall Measuring Mission (TRMM) “warm-season” (surface temperature greater than 10°C) precipitation radar (PR) vertical profiles, their surface rainfall, and associated radar-based classifiers (convective/stratiform and brightband existence). Twenty-five archetypal profile types are identified, including nine convective types, eight stratiform types, two mixed types, and six anvil/fragment types (nonprecipitating anvils and sheared deep convective profiles). These profile types are then hierarchically clustered into 10 similar families, which can be further combined, providing an objective and physical reduction of the highly multivariate PR data space that retains vertical structure information. The taxonomy allows for description of any storm or local convective spectrum by the profile types or families. The analysis provides a quasi-independent corroboration of the TRMM 2A23 convective/stratiform classification. The global frequency of occurrence and contribution to rainfall for the profile types are presented, demonstrating primary rainfall contribution by midlevel glaciated convection (27%) and similar depth decaying/stratiform stages (28%–31%). Profiles of these types exhibit similar 37- and 85-GHz passive microwave brightness temperatures but differ greatly in their frequency of occurrence and mean rain rates, underscoring the importance to passive microwave rain retrieval of convective/stratiform discrimination by other means, such as polarization or texture techniques, or incorporation of lightning observations. Close correspondence is found between deep convective profile frequency and annualized lightning production, and pixel-level lightning occurrence likelihood directly tracks the estimated mean ice water path within profile types.

Description: During its first three years, the Tropical Rainfall Measuring Mission (TRMM) satellite observed nearly six million precipitation features. The population of precipitation features is sorted by lightning flash rate, minimum brightness temperature, maximum radar reflectivity, areal extent, and volumetric rainfall. For each of these characteristics, essentially describing the convective intensity or the size of the features, the population is broken into categories consisting of the top 0.001%, top 0.01%, top 0.1%, top 1%, top 2.4%, and remaining 97.6%. The set of “weakest/smallest” features composes 97.6% of the population because that fraction does not have detected lightning, with a minimum detectable flash rate of 0.7 flashes (fl) min−1. The greatest observed flash rate is 1351 fl min−1; the lowest brightness temperatures are 42 K (85 GHz) and 69 K (37 GHz). The largest precipitation feature covers 335 000 km2, and the greatest rainfall from an individual precipitation feature exceeds 2 × 1012 kg h−1 of water. There is considerable overlap between the greatest storms according to different measures of convective intensity. The largest storms are mostly independent of the most intense storms. The set of storms producing the most rainfall is a convolution of the largest and the most intense storms. This analysis is a composite of the global Tropics and subtropics. Significant variability is known to exist between locations, seasons, and meteorological regimes. Such variability will be examined in Part II. In Part I, only a crude land–ocean separation is made. The known differences in bulk lightning flash rates over land and ocean result from at least two differences in the precipitation feature population: the frequency of occurrence of intense storms and the magnitude of those intense storms that do occur. Even when restricted to storms with the same brightness temperature, same size, or same radar reflectivity aloft, the storms over water are considerably less likely to produce lightning than are comparable storms over land.

Description: A taxonomy of tropical convective and stratiform vertical structures is constructed through cluster analysis of 3 yr of Tropical Rainfall Measuring Mission (TRMM) "warm-season" (surface temperature greater than 10 C) precipitation radar (PR) vertical profiles, their surface rainfall, and associated radar-based classifiers (convective/ stratiform and brightband existence). Twenty-five archetypal profile types are identified, including nine convective types, eight stratiform types, two mixed types, and six anvil/fragment types (nonprecipitating anvils and sheared deep convective profiles). These profile types are then hierarchically clustered into 10 similar families, which can be further combined, providing an objective and physical reduction of the highly multivariate PR data space that retains vertical structure information. The taxonomy allows for description of any storm or local convective spectrum by the profile types or families. The analysis provides a quasi-independent corroboration of the TRMM 2A23 convective/ stratiform classification. The global frequency of occurrence and contribution to rainfall for the profile types are presented, demonstrating primary rainfall contribution by midlevel glaciated convection (27%) and similar depth decaying/stratiform stages (28%-31%). Profiles of these types exhibit similar 37- and 85-GHz passive microwave brightness temperatures but differ greatly in their frequency of occurrence and mean rain rates, underscoring the importance to passive microwave rain retrieval of convective/stratiform discrimination by other means, such as polarization or texture techniques, or incorporation of lightning observations. Close correspondence is found between deep convective profile frequency and annualized lightning production, and pixel-level lightning occurrence likelihood directly tracks the estimated mean ice water path within profile types.

Description: During its first three years, the Tropical Rainfall Measuring Mission (TRMM) satellite observed nearly six million precipitation features. The population of precipitation features is sorted by lightning flash rate, minimum brightness temperature, maximum radar reflectivity. areal extent, and volumetric rainfall. For each of these characteristics, essentially describing the convective intensity or the size of the features, the population is broken into categories consisting of the top 0.001%, top 0.01%, top 0.1%, top 1%, top 2.4%. and remaining 97.6%. The set of weakest/smallest features composes 97.6% of the population because that fraction does not have detected lightning, with a minimum detectable flash rate of 0.7 flashes (fl) per minute. The greatest observed flash rate is 1351 fl per minute; the lowest brightness temperatures are 42 K (85 GHz) and 69 K (37 GHz). The largest precipitation feature covers 335 000 square kilometers and the greatest rainfall from an individual precipitation feature exceeds 2 x 10 kg per hour of water. There is considerable overlap between the greatest storms according to different measures of convective intensity. The largest storms are mostly independent of the most intense storms. The set of storms producing the most rainfall is a convolution of the largest and the most intense storms. This analysis is a composite of the global Tropics and subtropics. Significant variability is known to exist between locations. seasons, and meteorological regimes. Such variability will be examined in Part II. In Part I, only a crude land-ocean separation is made. The known differences in bulk lightning flash rates over land and ocean result from at least two differences in the precipitation feature population: the frequency of occurrence of intense storms and the magnitude of those intense storms that do occur. Even when restricted to storms with the same brightness temperature, same size, or same radar reflectivity aloft, the storms over water are considerably less likely to produce lightning than are comparable storms over land.

Description: A "dimensional reduction" (DR) method is introduced for analyzing lightning field changes whereby the number of unknowns in a discrete two-charge model is reduced from the standard eight to just four. The four unknowns are found by performing a numerical minimization of a chi-squared goodness-of-fit function. At each step of the minimization, an Overdetermined Fixed Matrix (OFM) method is used to immediately retrieve the best "residual source". In this way, all 8 parameters are found, yet a numerical search of only 4 parameters is required. The inversion method is applied to the understanding of lightning charge retrievals. The accuracy of the DR method has been assessed by comparing retrievals with data provided by the Lightning Detection And Ranging (LDAR) instrument. Because lightning effectively deposits charge within thundercloud charge centers and because LDAR traces the geometrical development of the lightning channel with high precision, the LDAR data provides an ideal constraint for finding the best model charge solutions. In particular, LDAR data can be used to help determine both the horizontal and vertical positions of the model charges, thereby eliminating dipole ambiguities. The results of the LDAR-constrained charge retrieval method have been compared to the locations of optical pulses/flash locations detected by the Lightning Imaging Sensor (LIS).

Description: We describe the clustering algorithm used by the Lightning Imaging Sensor (LIS) and the Optical Transient Detector (OTD) for combining the lightning pulse data into events, groups, flashes, and areas. Events are single pixels that exceed the LIS/OTD background level during a single frame (2 ms). Groups are clusters of events that occur within the same frame and in adjacent pixels. Flashes are clusters of groups that occur within 330 ms and either 5.5 km (for LIS) or 16.5 km (for OTD) of each other. Areas are clusters of flashes that occur within 16.5 km of each other. Many investigators are utilizing the LIS/OTD flash data; therefore, we test how variations in the algorithms for the event group and group-flash clustering affect the flash count for a subset of the LIS data. We divided the subset into areas with low (1-3), medium (4-15), high (16-63), and very high (64+) flashes to see how changes in the clustering parameters affect the flash rates in these different sizes of areas. We found that as long as the cluster parameters are within about a factor of two of the current values, the flash counts do not change by more than about 20%. Therefore, the flash clustering algorithm used by the LIS and OTD sensors create flash rates that are relatively insensitive to reasonable variations in the clustering algorithms.

Description: Robust description of the diurnal cycle from TRMM observations is complicated by the limitations of Low Earth Orbit (LEO) sampling; from a 'climatological' perspective, sufficient sampling must exist to control for both spatial and seasonal variability, before tackling an additional diurnal component (e.g., with 8 additional 3-hourly or 24 1-hourly bins). For documentation of vertical structure, the narrow sample swath of the TRMM Precipitation Radar limits the resolution of any of these components. A neural-network based 'virtual radar" retrieval has been trained and internally validated, using multifrequency / multipolarization passive microwave(TM1) brightness temperatures and textures parameters and lightning (LIS) observations, as inputs, and PR volumetric reflectivity as targets (outputs). By training the algorithms (essentially highly multivariate, nonlinear regressions) on a very large sample of high-quality co-located data from the center of the TRMM swath, 3D radar reflectivity and derived parameters (VIL, IWC, Echo Tops, etc.) can be retrieved across the entire TMI swath, good to 8-9% over the dynamic range of parameters. As a step in the retrieval (and as an output of the process), each TMI multifrequency pixel (at 85 GHz resolution) is classified into one of the 25 archetypal radar profile vertical structure "types", previously identified using cluster analysis. The dynamic range of retrieved vertical structure appears to have higher fidelity than the current (Version 6) experimental GPROF hydrometeor vertical structure retrievals. This is attributable to correct representation of the prior probabilities of vertical structure variability in the neural network training data, unlike the GPROF cloud-resolving model training dataset used in the V6 algorithms. The LIS lightning inputs are supplementary inputs, and a separate offline neural network has been trained to impute (predict) LIS lightning from passive-microwave-only data. The virtual radar retrieval is thus, in principle, extensible to Aqua/AMSR-E and NPOESS/CMIS passive microwave instruments. The virtual radar approach yields a threefold increase in effective sampling from the mission, albeit of lower-quality "retrieved" data, reducing the variance of local estimates by one third (or the standard deviation by-0.57). In this talk, the variance reduction is leveraged to more finely resolve global diurnal variability in both space and time (local hour).

Description: Lightning optical flash parameters (e.g., radiance, area, duration, number of optical groups, and number of optical events) derived from almost 5 yrs of Optical Transient Detector (OTD) data are compared with peak current and multiplicity observations derived from the US National Lightning Detection Networkm (NLDN). Despite the relatively low lightning geolocation accuracy afforded by OTD, a total of 48,870 NLDN cloud-to-ground (CG) flashes were correlated with OTD flashes, or about 10,000 CGs per year. The median values of the above OTD flash parameters for the 48,870 CGs were, respectively: 0.137 J/square meters/sr/micrometers, 313.7 square kilometers, 0.189 s, 4 optical groups per CG, and 8 optical events per CG. Invoking the multiplicity data, the median number of optical groups per stroke was 2.5, and the median number of optical events per stroke was 5.0. Median values of peak current for negative and positive CGs were -21.6 kA and 17.8 kA, respectively, and as expected, the negative CGs had a larger average multiplicity than the positive CGs. A statistical summary is provided for all CGs, for positive and negative CGs, and for CGs from different seasons. Standard two-distribution hypothesis tests were perfonned to intercompare the population means of the various lightning parameters. In particular, and to greater than the 99% confidence level, it was found that positive CGs are on average more radiant, of greater areal extent, and are longer lasting than negative CGs. Rankings from a complete set of hypothesis tests between CGs of different polarities and from different seasons are also provided. Most notably, wintertime positive CGs tend to be more radiant, of greater areal extent, and longer lasting than any other group of CGs (i.e., negative springtime CGs, positive summertime CGs, etc.).

Description: Joint observations from the TRMM mission are used to examine the contribution which lightning observations can make towards estimation of vertical structure and rainfall from passive microwave data. Passive microwave (TMI) observations are binned into distinct vertical structure categories based on cluster analysis of radar (PR) vertical profiles and convective/stratiform (C/S) classifiers. TMI rain estimates are high-biased relative to radar for stratiform and mixed convective/stratiform vertical structures and low-biased relative to radar for convective structures. Significant ambiguity exists in the TMI brightness temperature space between midlevel stratiform and convective profile types which are primary contributors to tropical rainfall. The ambiguity is worst for convective profile types with radar echo tops within and just above the mixed-phase region. The ability of TMI data (including all low and high frequency polarized brightness temperatures; 19V, 19H, 21V, 37V, 37H, 85V, 85H) to predict vertical structure is assessed. The incremental benefit of including TMI convective/stratiform (C/S) classifiers (both 85 GHz polarization and 19, 37 and 85 GHz texture-based) and LZS lightning observations is then considered. The use of all parameters (brightness temperatures, TMI C/S classifiers, lightning data, predicted vertical structure) to reduce TMI/PR rainfall estimate scatter is demonstrated. For both vertical structure and rainfall estimation, inclusion of C/S classifiers and lightning observations provide statistical skill improvements but do not completely alleviate critical midlevel profile C/S ambiguity and related rain estimate errors.