ALBERT

Description: A diagnostic analysis of the VVP (volume velocity processing) retrieval method is presented, with emphasis on understanding the technique as a linear, multivariate regression. Similarities and differences to the velocity-azimuth display and extended velocity-azimuth display retrieval techniques are discussed, using this framework. Conventional regression diagnostics are then employed to quantitatively determine situations in which the VVP technique is likely to fail. An algorithm for preparation and analysis of a robust VVP retrieval is developed and applied to synthetic and actual datasets with high temporal and spatial resolution. A fundamental (but quantifiable) limitation to some forms of VVP analysis is inadequate sampling dispersion in the n space of the multivariate regression, manifest as a collinearity between the basis functions of some fitted parameters. Such collinearity may be present either in the definition of these basis functions or in their realization in a given sampling configuration. This nonorthogonality may cause numerical instability, variance inflation (decrease in robustness), and increased sensitivity to bias from neglected wind components. It is shown that these effects prevent the application of VVP to small azimuthal sectors of data. The behavior of the VVP regression is further diagnosed over a wide range of sampling constraints, and reasonable sector limits are established.

Description: The Lightning Imaging Sensor (LIS) is a NASA Earth Observing System (EOS) instrument on the Tropical Rainfall Measuring Mission (TRMM) platform designed to acquire and investigate the distribution and variability of total lightning (i.e., cloud-to-ground and intracloud) between q35' in latitude. Since lightning is one of the responses of the atmosphere to thermodynamic and dynamic forcing, the LIS data is being used to detect deep convection without land-ocean bias, estimate the precipitation mass in the mixed phased region of thunderclouds, and differentiate storms with strong updrafts from those with weak vertical motion.

Description: The Optical Transient Detector (OTD) is a space-based instrument specifically designed to detect and locate lightning discharges (intracloud and cloud-to-ground) as it orbits the Earth. A statistical examination of OTD lightning data reveals that nearly 1.2 billion flashes occurred over the entire earth during the one year period from September 1995 through August 1996. This translates to an average of 37 lightning flashes occurring around the globe every second, which is well below the traditional estimate of 100 flashes per second. An average of 75% of the global lightning activity during the year occurs between 30' S and 30' N. An analysis of the annual lightning distribution reveals that an average of 82% of the lightning flashes occur over the continents and 18% over the oceans, which translates to an average land-ocean flash density ratio of nearly 11.

Description: The mapping of the lightning optical pulses detected by the Lightning Imaging Sensor (LIS) is compared with the radiation sources by Lightning Detection and Ranging (LDAR) and the National Lightning Detection Network (NLDN) for three thunderstorms observed during and overpasses on 15 August 1998. The comparison involves 122 flashes including 42 ground and 80 cloud flashes. For ground flash, the LIS recorded the subsequent strokes and changes inside the cloud. For cloud flashes, LIS recorded those with higher sources in altitude and larger number of sources. The discrepancies between the LIS and LDAR flash locations are about 4.3 km for cloud flashes and 12.2 km for ground flashes. The reason for these differences remain a mystery.

Description: Since April 1995, lightning activity around the globe has been monitored with the Optical Transient Detector (OTD). The OTD observations acquired during the one year period from September 1995 through August 1996 have been used to statistically determine the number of flashes that occur over the Earth during each hour of the diurnal cycle, expressed both as a function of local time and universal time. The globally averaged local [il,htnina activity displays a peak in late afternoon (1500-1800 local time) and a minimum in the morning hours (0600- 1000 local time) consistent with convection associated with diurnal heating. No diurnal variation is found for oceanic storms. The diurnal lightning distribution (universal time) for the globe displays a variation of about 35% about its mean as compared to the Carnegie curve which has a variation of only 15% above and below the mean.

Description: In the framework of the project POLINAT 2 (Pollution in the North Atlantic Flight Corridor) we measured NO(x) (NO and NO2) and ozone on 98 flights through the North Atlantic Flight Corridor (NAFC) with a fully automated system permanently installed aboard an in-service Swissair B-747 airliner in the period of August to November 1997. The averaged NO, concentrations both in the NAFC and at the U.S. east coast were similar to that measured in autumn 1995 with the same system. The patchy occurrence of NO(x), enhancements up to 3000 pptv over several hundred kilometers (plumes), predominately found over the U.S. east coast lead to a log-normal NO(x) probability density function. In three case-studies we examine the origins of such plumes by combining back-trajectories with brightness temperature enhanced (IR) satellite imagery, with lightning observations from the U.S. National Lightning Detection Network (NLDN) or with the Optical Transient Detector (OTD) satellite. For frontal activity above the continental U.S., we demonstrate that the location of NO(x) plumes can be well explained with maps of convective influence. For another case we show that the number of lightning flashes in a cluster of marine thunderstorms is proportional to the NO(x) concentrations observed several hundred kilometers downwind of the anvil outflows and suggest that lightning was the dominant source. From the fact that in autumn the NO, maximum was found several hundred kilometers off the U.S. east coast, it can be inferred that thunderstorms triggered over the warm Gulf Stream current are an important source for the regional upper tropospheric NO(x) budget in autumn.

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.