This paper describes and evaluates a satellite rainfall estimation technique that combines infrared and lightning information to estimate precipitation in deep convective systems. The algorithm is developed and tested using seven years (2002-08) of TRMM measurements over the southern
The Geostationary Operational Environmental Satellite R series (GOES-R, to be launched in 2015; Gurka et al. 2006) will carry an Advanced Baseline Imager (Schmit et al. 2005) and the first-ever Geostationary Lightning Mapper (GLM; Christian 2008). GLM is modeled after the
Most IR-based techniques employ empirical regressions to derive surface rainfall rates from cloud-top brightness temperatures (Tb; Arkin 1979; Negri and Adler 1981; Arkin and Meisner 1987; Griffith et al. 1978; Scofield 1987; Adler and Negri 1988; Goodman et al. 1994; Kuligowski 2002). Although some IR techniques consider additional factors (e.g., cloud growth and cloud-top structure; Adler and Negri 1988; Vicente et al. 1998), combine multiple spectral channels (Sorooshian et al. 2000; Ba and
Previous rainfall retrieval studies have demonstrated lightning information as having a great potential in improving the estimation of convective rainfall (
This study aims to develop and evaluate an IR- lightning combined rainfall estimation algorithm using spaceborne lightning and IR measurements from the
2. Data and methodology
a. TRMM dataset
This study applies three years (2002-04) of TRMM measurements to develop and tune an IR-lightning combined rainfall algorithm. We then use both the training (2002-04) and an independent (2005-08) dataset of TRMMrainfall retrievals to evaluate this particular technique. To be fairly evaluated by the independent data (4 yr), both datasets are selected to be after the TRMM orbit altitude boost in 2001. Our interests in this study focus on the warm-season (May-August) precipitation over the southern continental
The TRMM version-7 dataset is used in this study. TRMM measurements are from four instruments (Kummerow et al. 1998), including the precipitation radar (PR), TRMM Microwave Imager (TMI), Lightning Image Sensor (LIS; Christian 1999), and Visible and IR Scanner (VIRS). The PR provides threedimensional storm information of hydrometeors (radar reflectivity) at approximately 5-km horizontal and 250-m vertical resolution. In this study, we mainly use PR-based convective/stratiform rainfall separation (2A23; Kozu et al. 2001; Awaka et al. 2009) and rainfall rates (2A25) (
The VIRS is a five-channel imaging spectroradiometer scanning the earth through a 720-km swath with 2.11-3.02-km footprints (Kummerow et al. 1998). This study applies brightness temperatures from the 10.8-mm window channel to develop the IR-based rainfall estimation technique. The LIS is an optical sensor, which detects total lightning (intracloud and cloud-to-ground lightning). LIS monitors total lightning activity in a 600 3 600km2 area during a period between 40 and 90 s with a footprint of ;5 km and temporal resolution of 2ms (Christian 1999). The LIS algorithm separates lightning events (resolution at 5km and 2ms) with strong irradiance from the background and defines lightning flashes as clusters of adjacent events that occur within 330ms. LIS detects 70%-90% of all lightning flashes within its field of view (FOV). The GOES-R GLM is modeled after TRMM LIS and is thus very similar to LIS, except that GLM is a staring sensor with a broad FOV (13 000 3 13 000km2). GLM's scheduled resolution will range between 8 and 12km over
b. Algorithm characteristics
Two major steps are involved in this algorithm development: the first step is to develop an IR-only rainfall estimation technique while the second step is to couple lightning information into that IR-based algorithm. The convective-stratiform technique (Adler and Negri 1988) is an IR-based rainfall estimation technique that has been applied to geostationary data. Because the CST already separately estimates convective core rainfall and anvil stratiform rainfall, it is a good candidate for the base algorithm, since the lightning information is meant to improve only the convective rain estimation. The flowchart of the algorithms is given in Fig. 1. First of all, we apply a version of the IR-based CST, with parameters determined from three years of TRMM VIRS IR data and PR and TMI rain retrievals. The CST convective portion (second column in Fig. 1) includes finding local minima in the Tb field (Tmin), conducting a slope test (Tb difference between the Tmin and surrounding points) on identified Tmin for final decision of location of convective cores, defining the convective area (as a function of Tmin and centered on the Tmin location), and assigning the convective rain rate (also as a function of Tmin over the convective area).
Lightning information is then combined with the IRbased CST to develop a lightning-enhanced CST (third column in Fig. 1). Lightning information is used to confirm or eliminate IR-based convective cores from Tmin features that already passed the slope test, identify convective areas from lightning flash rate features, even if no IR Tmin is collocated with the lightning, and assign convective rain rate (as a function of both Tmin and lightning flash density). Stratiform rain areas are determined by a background or mode Tb (Tmode) in the anvil, which is close in value to the tropopause temperature (Adler and Negri 1988). Therefore, the stratiform rain areas are similar for both CST and CSTL (first column in Fig. 1). In the order of defining raining areas, convective core areas take precedence over stratiform areas. Any convective area defined by the IR data that fails to pass the lightning test will go through the stratiform test again. Therefore, if an IR-defined convective core is eliminated by the lightning test but reaches the stratiform threshold, it will be reassigned as stratiform. Otherwise, the IR-based convective core area will be assigned as a no-rain area. Finally, the stratiform areas are assigned a median stratiform rain rate (;2.5mmh21) calculated from the training dataset. It should be noted that all the functions used in CST and CSTL are empirically defined and tuned by the training dataset used in this study.
c. Functions and parameters applied in the algorithm
First of all, the CST identifies local minima (Tmin) of Tb. These Tmin points are considered as the center points of potential convective cores. Certain empirical relationships are derived from the training dataset and used as part of the CST and CSTL. These functions include the slope test and relations between convective area and Tmin, convective rain rate and Tmin, stratiform rain rate and Tb, and convective rain rate and lightning. They are specifically defined as follows:
1) A slope function is defined as a function that linearly discriminates convective and nonconvective cores through linear discriminant analysis (Fig. 2a). A slope is defined as the difference between the local minimum of Tb (Tmin) and the mean Tb of the surrounding pixels within a 10-km radius [similar to Adler and Negri (1988)]. Convective cores are defined as Tmin cores that have at least two PR-defined convective pixels located within 10-km radius, while nonconvective cores have no PR convective pixels within that distance.
2) Convective area and rain rate over that area are further connected to Tmin, and finally derived from the regression of Tmin (Figs. 3a,b). The convective area is derived as the total number of PR convective pixels that are within a 25-km radius of the Tmin location (Fig. 3a). On the other hand, the convective rain rate of each pixel within the convective core area is derived from the relationship betweenPR maximum rain rate of the convective core and Tmin (Fig. 3b).
3) The stratiform rainfall area is defined by the background anvil stratiform threshold temperature (or Tmode) of the mature convective systems in the image (Adler and Negri 1988). Tmode has been used previously by Adler et al. (1985) as a background Tb and an indicator of the thick portion of mature anvils where precipitation is present. Tmode is calculated for each Tmin with slope is less than or equal to 4K, since most Tmin reaching this slope criterion are mature thunderstorms found embedded in extensive cirrus background (Adler et al. 1985). Specifically, here this background Tb, or Tmode, is defined as the weighted mean of Tb surrounding ;80km of the Tmin: Tmode 5 E[W(i) 3 Tb(i)], where E is the expected value and the weight W(i) is the number of IR pixels at Tb(i). The 80-km surrounding area is selected because this area is large enough to obtain a stable statistic, yet small enough to be representative of a particular Tmin point. The final Tmode within the satellite image is derived from the average of all Tmin mode temperatures. The same Tmode temperature is applied to all pixels within a single IR image. Pixels with Tb at and below the Tmode value that have not been assigned a convective rain rate are given a fixed stratiform rain rate. Using the training dataset this value is derived from the PR data as 2.5mmh21. Figure 3c shows that the distribution of PR rain rate in these stratiform areas as a function of VIRS Tb is nearly constant, indicating the significance of applying the constant rain rate value over a range of Tb.
4) For locations with lightning flashes, the convective rain rate is derived from lightning flash density. Based on Xu et al. (2013), there is no clear regression relationship between lightning and the PR-defined rain rate, probably related to location differences between the surface radar echo and lightning locations at fine spatial resolution. This might also be due to efficient autoconversion and/or coalescence (warm rain) processes in storm cells under certain environments. On the other hand, the close correlation between lightning flash density and TMI rain rate is probably due to both the microwave measurements and the lightning occurrence having a physical connection to mixed-phase precipitation microphysics, but with only indirect connections to surface precipitation. Therefore, the TMI rain-rate estimate is used here to connect with lightning. The passive microwave (TMI) rain cores are better collocated with locations of lightning flash density cores. The rain rate-lightning function is derived from the mean TMI rain rate and lightning flash density within a 10-km radius (Fig. 3d).
All constants in the above empirical functions are based on sensitivity tests, validation, and statistics. For example, multiple tests have been done for choosing a window (including 25-, 50-, 75-, and 100-km radius) to define convective core area. These tests have been applied to the 3-yr training dataset and evaluated by the PR convective area. The validation results show that the 25-km radius window is the best window for this purpose. A 10-km area for lightning area is also chosen in this way. Sensitivity tests have also been done for using a specific lightning density for the restriction of convective core. In this study, a density of two flashes per minute within 10-km radius ends up to be the best proxy. Sensitivity tests show that a five flashes per minute criterion would eliminate 50% of the convective cores and would not allow the algorithm to correctly represent the convective core population or area. We recognize that a true convective core, especially a wide convective core, might be separated into multiple "convective cores" in the CST algorithm. This might make the criterion of lightning density lower for identifying convective cores. Similarly, Xu et al. (2013) showed that the threshold of two flashes per minute restricts 90% of the precipitation cores to be convective cores (having convective elements).
There are two reasons that TMI rain rates are used to derive the relationships between rain rate and lightning in item 4 above. First, our main goal of using lightning is to improve the IR rain estimates up to the accuracy of microwave estimates, as lightning measurements have close relationships with the microwave measurements. If this goal is fulfilled, the IR rain estimates already have benchmark improvements. Second, when the TMI rain rates are used for the lightning- rainfall regression, then it will also be easy to do the validation/calibration by microwave estimates in the future. This is also straightforward, since calibration by passive microwave measurements is proposed for the final version of CSTL that will be applied toGOES-Rdata.
d. Evaluation methods and definitions
Dependent (2002-04) and independent (2005-08) TRMM PR and TMI estimates are used to evaluate the characteristics and accuracy of the CST and CSTL techniques. Functions and coefficients used in the algorithm are derived from the dependent data. The algorithm is then applied back to the dependent data and further to the independent data. The CST and CSTL are evaluated in terms of convective-stratiform rain identification, instantaneous rain-rate estimates, and total rain volume. Accuracy statistics of correlation coefficient (CC), bias, and root-mean-square error (RMSE) are calculated. This study also examines contingency score parameters including the probability of detection (POD), false alarm rate (FAR), and critical success index (CSI). Specifically, the success-failure parameters are defined as POD 5 Nsuccess/(Nsuccess 1 Nfailure), FAR 5 Nfalse alarm/(Nsuccess 1 Nfalse alarm), and CSI 5 Nsuccess/(Nsuccess 1 Nfailure 1 Nfalse alarm).
3. Details of the CST and CSTL algorithms
a. IR-based convective/stratiform technique
The basic idea of the CST is that cloud top above the convective core is usually colder than the surrounding clouds as convective updraftdevelops and pushes more ice particles to the upper level (Zipser 1977; Houze 1989). As a result, local minima of brightness temperatures (Tmin) can represent the location of convective cores (and therefore heavy precipitation regions) better than assigning rain rate as a function of brightness temperatures through a regressionmethod (the colder the brightness temperature is, the heavier is the rain rate). The CST algorithm first searches for local Tmin cores. Tmin cores are then examined by the slope test. If a Tmin core passes the slope test, it is identified as a convective core.
A function used for the slope test is first defined using the statistics of the training dataset. The slope parameter is defined as the difference between the mean of the Tb values surrounding the Tmin location and the Tmin value itself. This slope value is then compared to whether or not the PR has convective rain identified at that location. The slope function is then defined as a function that linearly discriminates the PR-defined convective and nonconvective Tmin features (cloud) based on VIRS and
If the measured slope value of a Tmin core is larger than the corresponding value in the slope function (Fig. 2a) and absolutely greater than 2K, then it passes the slope test and is accepted as a location of a convective core; otherwise it fails. Given the sensor sensitivity, noise of the IR brightness field, and oscillation of cloudtop field due to factors such as turbulence, a 2-K gradient seems to be a reasonable lower value. Furthermore, some sensitivity tests have been done for this selection, and 2Kworks out to be the lowest required gradient for a convective core using the CST technique.As shown in Fig. 2a, convective and nonconvective (relative to PR) Tmin features generally fall into two groups. These two groups are mainly discriminated at the temperature of 210K for low-slope Tmin cores (,6K) and at the temperature of 220K for high-gradient (.12 K) Tmin cores. The performance of the slope function defined here is quite stable when applied to different years or months in either the training or independent dataset (Fig. 2b), with both the POD and FAR identification of convective cores being relatively invariant month by month. In addition, the detection of convective cores is significantly high (POD 5 0.6-0.8), while the false alarm is at a reasonable level (FAR 5 0.3-0.5). However, the technique obviously misses some convective cores and allows room for additional information to improve the detection algorithm.
Convective cores (Tmin features passing the slope test) are further assigned convective rain area (Fig. 3a) and rain rate (Fig. 3b) as a function of their Tmin values. These functions are derived from the median values of convective area and convective rain rate of PR estimates as a function of Tmin, respectively. In the original CST (Adler and Negri 1988), similar functions (slope, convective area-Tmin, and rain rate-Tmin) were derived from one-dimensional cloud model results. In the CST stratiform rain areas are determined as those with Tb values lower than the mode or background brightness temperature, which has been shown to be coincident with the outline of stratiform rainfall and approximately equal to the tropopause temperature (Adler and Negri 1988). The stratiform rain rate is assigned as the median stratiform rainfall rate (2.5mmh21) derived from the PR rain estimates in the training dataset (Fig. 3c). As shown in Figs. 3a-c, convective rain area and rain rate have strong relationships with Tmin, but the stratiform rain rate is independent of Tb. These relationships also indicate the importance of separating convective and stratiform precipitation, even for IR-based techniques.
Figure 4 shows an example of the CST applied to a young thunderstorm complex from the step of finding local Tb minima (Fig. 4a) to assigning convective/ stratiform areas (Fig. 4c) and to deriving rain-rate estimates (Fig. 4d). Generally, the CST can catch the convective cores and define the heavy precipitation (convective) areas very well for young convective cells (when compared with radar; Fig. 4b). In such new developing systems, cloud-top overshooting is very evident, and the convective area is closely related to theminimum cloud-top temperature.However, as the convective storm develops and becomes mature, the relationships between IR Tb and rain rate in the CST (and most IR techniques) become more complicated (Fig. 5). As shown by the example in Fig. 5c and previous studies (Adler and Negri 1988; Xu et al. 2013), the CST may incorrectly identify cold cloud debris of mesoscale convective storms (MCS) or thick cirrus as convective areas and define too many convective cores in the broad, cold Tb field. Some of the Tmin coresmight result fromthe oscillation of the Tb field due to forced tropospheric mesoscale gravity waves generated by overshooting tops (Fovell et al. 1992) rather than true overshooting cloud tops. As shown in the next section, this situation could be corrected through combining lightning information with the CST.
b. Lightning-enhanced CST
The IR-based CST, described in the last section, is combined with lightning measurements to develop the CSTL. The use of lightning information is based on the lightning-cloud-rainfall relationships quantified by Xu et al. (2013). Xu et al. (2013) showed that the lightning information is useful in aiding the identification of convective cores in thick anvils missed by the IR technique, eliminating misidentified convective cores in cold cloud debris or thick cirrus, and improving convective rainfall volumeestimates. In addition, the effective lightning flash area also has a very good correlation with the convective rainfall area (Grecu et al. 2000; Xu et al. 2013).
In the lightning-enhanced algorithm, lightning information is incorporated into the CST after Tmin cores are identified and pass the slope test (Fig. 1) using just IR data. In the CSTL, Tmin cores that pass the CST slope test are further tested and must have at least two lightning flashes (in a 1-min period) within 10-km radius of the Tmin to remain as convective cores. Areas of these convective cores are derived from the area-Tmin function (Fig. 3a). In addition, lightning-associated areas (areas within a 10-kmradius of lightning flash centers) are also considered as convective areas, since lightning flash areas are highly correlated to convective area (Xu et al. 2013). No areas survive as being convective after failing these two kinds of lightning test. The lightning clusters add an increment to the IR convective cores (with more than two flashes) rather than dominating the numbers. Grecu et al. (2000) showed that lightning clusters only defined 50% of the true convective area. In most cases lightning clusters are much smaller than the IR-defined convective cores. For example, even for a convective core having 10 lightning flashes, these lightning flashes might overlap over a relatively smaller area than the large convective core. On the other hand, lightning-defined areas help to redefine real convective areas missed by CST in the situation when CST is misled by large uniform cold cloud cover over the convective core. Although stratiform or anvil lightning flashes do occur, they are only a very small fraction of the total lighting and most of them are related to a convective core nearby (more details in section 5b).
The Tmin-based convective areas (with nearby lightning) plus those 10-km areas with lightning (even with no IR Tmin) together make up the estimated convective rain areas. The example in Fig. 5 shows obvious improvement (CSTL over CST) as compared to the passive microwave (Fig. 5b) in identifying convective areas through the use of the lightning information. The lightning information helps to remove convective regions incorrectly defined by the CST (Fig. 5d). This incorrect definition might be due to cold temperatures and Tmin features appearing due to fluctuations and oscillations in thick anvil cirrus (Xu et al. 2013). Since we need an overall large view of the case, TMI retrievals instead of PR are shown here. There is some difference even between PR and TMI as seen from the heavy precipitation areas in Fig. 6. Based on Fig. 6, CSTL (or TMI) catches most convective precipitation (e.g., .10mmh21) that has been detected by the PR, except for the narrow line over the most western part of the system. In comparison with PR, CST still generates much more convective precipitation.
After defining convective areas, convective rain rate is assigned as a function of Tmin and lightning flash density. For convective areas without lightning occurrence, rain rate is derived from the rain rate-Tmin function as shown in Fig. 3b. In addition, rain rate is also derived as a function of lightning flash density in convective areas having lightning flashes (Fig. 3d). The final estimates of instantaneous rainfall are degraded to 10- or 20-km resolution (Fig. 6). As shown in the example, both CST and CSTL generally reproduce the distribution of heavy and light (stratiform) rain even under the large cold cloud shield. CST artificially overestimates the rainfall in some no-rain or very light rain regions (Fig. 6a) where CST falsely defines convective features (Fig. 5c). In contrast, the CSTL helps to solve this problem (Figs. 5d and 6d) and produces more accurate rain-rate patterns. As also shown in this example, CSTL estimates even catch those extremely heavy (.30mmh21) precipitation locations, especially when they are compared with TMI estimations. On the other hand, CST and CSTL also show skill in identifying the light rain (stratiform) area or delineating the rain and no-rain boundary in this mature MCS as shown in Figs. 5 and 6.
4. Evaluations of CST and CSTL rain estimates
This section evaluates CST and CSTL rain estimates by using TRMM PR and TMI rain retrievals dependently (using the TRMM 2002-04 training set) and independently (using 2005-08 TRMM data). Both the dependent and independent datasets include over a million PR- and TMI-defined precipitation pixels (Table 1). To fully evaluate the modified CST and the value of incorporating lightning information, the CST and CSTL estimates are compared with PR and TMI estimations on the basis of convective-stratiform separation, instantaneous rainfall rates, and total rainfall volume. Ground-based rain rates are not used for validation in this study for two major reasons: 1) the current version of the algorithm is a prototype and is for testing the value of adding lightning data (ground-based measurements will be included in later validation of products) and 2) time- space matching is easy and accurate when all measurements are on the same platform (TRMM satellite), whereas matching to ground-based measurements would be much more difficult.
a. Convective and stratiform identification
As mentioned previously, cloud-top Tb-rainfall and lightning-rainfall relationships are totally different in the convective and stratiform regions. One of the advantages of CST and CSTL is that they define convective and stratiform precipitation type before assigning specific rain rates. Figure 7 shows their POD, FAR, and CSI in identifying convective-type rainfall (with reference to PR-based identification). All comparisons are done at 10-km resolution; the contingent performance scores are even higher for the 20-km resolution. In general, both the CST and CSTL show significant accuracy in identifying convective rain features and the level of accuracy is invariant between dependent and independent estimates. First of all, the CST already has substantial ability (POD 5 0.7) to determine convective precipitation regions even before lightning information is used (Fig. 7a). This is reasonable, since overshooting cloud tops and associated minimum Tb features (defined as convective areas in CST) often correspond to heavy precipitation near the surface. Lightning information slightly improves (0.05 or 8% increase in POD) the convective detection after it is incorporated into CST. The combination of both IR and lightning should be superior to IRonly or lightning-only methods (Grecu et al. 2000). In addition, the improvement by using total lightning versusCGis quite obvious. For example, Grecu et al. (2000) reported that CG lightning cluster method detects only 36%(POD) of the convective areas. Therefore, only using CG lightning area as a convective proxy could largely underestimate the convective areas. The relatively small increase in POD from CST to CSTL is due to the relatively few, but probably important, convective features that are missed just using IR under a dense anvil cirrus deck. In fact, the most important role of lightning information seems to be to screen out convective areas incorrectly defined by CST (Fig. 7b). Specifically, lightning data help to reduce the false alarms in convective rain by about 30% (i.e., the FAR drops from 0.5 for CST to 0.35 in CSTL). Because of the increase of convective detection and false alarm decrease, CSTL improves its CSI over that of CST by about 25% (Fig. 7c). These statistics clearly support our hypothesis that lightning information can be used to identify convective cores (areas) that are missed by CST and remove CST false alarms on convective cores (areas). These results are critical in understanding how the lightning information can improve convective system rain estimation.
Figure 8 shows convective areas defined by CST and CSTL versus those derived from PR (general statistics are in red text). Each point in the figure is the total convective area over the 800 3 800km2 TRMM sector image, with a number of individual convective features contributing to the total area for that image. In general, both CST (CC50.86) and CSTL (CC50.92) show great skill in representing convective features over a region. The CST has a fairly small overall bias (15%) but does generate larger (up to 100%) convective precipitation areas than the corresponding PR estimates for large precipitating system (Figs. 8a,b). This area overestimation at high values is due to the high convective false alarm rate of CST at the latter stage of storms such as MCSs when cold and complicated cloud debris mislead the IR-based technique searching for Tmin features. This is a common issue for IR-only rainfall estimation techniques. Since lightning information helps to lower the convective identification FAR in the CST (Fig. 7), the CSTL produces more accurate estimates of convective area (Figs. 8c,d). For example, CSTL lowers the overall bias on independent PR convective area estimation from 16% to 22%, while reducing the RMSE from 0.95 to 0.66 (3104km2). Both the CST and CSTL display very small changes in statistics between the dependent and independent periods (2002-04 vs 2005- 08) in deriving convective area. The correlation coefficient between CST (or CSTL) and PR estimates of convective area is very high (0.86-0.92). The addition of the lightning information, however, does improve the bias and RMSE of the identification of the convective areas (statistics shown in Fig. 8).
The stratiform area is estimated primarily from the IR data using a thresholding technique in both the CST and CSTL. Figure 9 shows that IR techniques in this study have significant skill in estimating stratiform rain area (CC 5 0.78-0.82). However, both CST and CSTL estimates have a somewhat larger scatter and higher bias in the identification of stratiform areas than convective areas. This reflects the difficulty of defining light rain boundaries with IR information. The lightning information improves the stratiform area estimation (e.g., reduces the bias by 15%-20%), although the scatter (as seen in CC and RMSE values) is not significantly affected by the lightning information. The bias reduction is due primarily to improving the convective area estimation (e.g., by eliminating false convective signatures), and thereby in turn improving the remaining stratiform area definition. As expected, the primary impact of the inclusion of (total) lightning data in this algorithm lies in improving estimations in convective areas, not stratiform areas.
b. Estimates of instantaneous rainfall rates
Figures 10 and 11 display the CST and CSTL rain estimates in comparison with PR-TMI estimates at 20- km resolution in
The results above are also consistent for different months (May-August) and for different years (see Fig. 12). Table 2 and Fig. 12 show the complete set of validation statistics (i.e., correlation coefficient, bias, and RMSE). In general, the CSTL improves the statistics over CST (Table 2), with the correlation coefficient increasing by more than 30% when evaluating by PR (CC increases from 0.36 to 0.49) or TMI (CC increases from 0.53 to 0.66) retrievals. CSTL reduces the bias over CST by 30% (50%) and RMSE by 20% (25%) when compared with PR (TMI) estimates. Furthermore, lightning information also helps to improve the POD and lowers the FAR. For example, in the validation by TMI estimates with rainfall rates . 5mmh21, CSTL increases the POD slightly (from 0.67 to 0.70) and greatly reduces the FAR (from 0.64 to 0.52) over CST. These statistics further support the conclusion in the earlier section that lightning does help to improve detection and eliminate false alarms in convective (or heavy) precipitation.
All results (Figs. 10-12 and Table 2) show that both the CST and CSTL have significant skill in instantaneous rainfall estimation and the CSTL has better statistical performance due to the addition of the lightning information. CST estimates show a correlation coefficient of 0.53, bias of 0.8mmh21, and RMSE of 6.5mmh21, when evaluated by comparing to passive microwave estimates (Table 2). Compared with several IR techniques at similar time and space scales reported by Kuligowski (2002), the performance score of CST falls into the best group. Lightning-enhanced CSTL shows even more accurate rainfall estimates, with a correlation coefficient of 0.66 and bias of only 0.4mmh21. These improvements on IR-only estimates by lightning information are comparable to previous studies (Grecu et al. 2000; Morales and Anagnostou 2003; Chronis et al. 2004), although this study focuses on a much larger dataset.
c. Total rainfall volume estimations
This section evaluates the overall performance of the CST and CSTL by considering both the precipitation area and intensity. Accurate estimates of the rainfall total (or volumetric rainfall) over a region are important for hydrological applications and even climate studies. Figure 13 shows the instantaneous rainfall volume estimations by CST and CSTL as a function of TRMM PR-TMI estimates in each 800 3 800 km2 box or image over the southern
5. Summary and discussion
This study develops and evaluates an IR-lightning combined rainfall estimation algorithm using TRMM measurements over the southern
Major results from this study include the following: 1) Both CST and CSTL show significant skill and stable performance in warm-season rainfall estimates over the southern
The use of lightning shows significant improvement over IR-based rainfall estimates as also shown in previous studies (Grecu et al. 2000;Morales andAnagnostou 2003; Chronis et al. 2004).Asmentioned previously,most IR-based techniques employ empirical regressions to derive surface rainfall rate from cloud-top brightness temperatures. However, the relationship between cloudtop Tb and surface rain rate is subject to great variability. Precipitating systems with efficient warm rain processes are great examples, whose surface rain rates are high but updrafts are too weak to transport ice particles to the upper troposphere; therefore the cloud tops are warm. On the storm scale, the Tb-rain rate relationship changes at different stages of the storm life cycle and in different parts of the storm. If the same Tb-rain rate relationship is applied throughout the storm life cycle, IR techniques could frequently underestimate rainfall at early stages when the system is dominated in terms of area by convective cores at high rain rates and overestimate rainfall at later stage, when the system is mostly stratiform by area. Similarly, the IR technique would overestimate stratiform precipitation and underestimate convective rainfall volume when the same Tb-rain rate function is derived from the whole precipitation area. The CST IR approach takes these variations into account, with the result that a CST pixel with the same Tb could have different rain rates depending on whether it is defined as convective or stratiform. In addition, most IR techniques could falsely treat thick cirrus clouds or cloud anvils as active convection, and miss convective cores under large uniform cold cloud shields or sheared cloud tops. The CST itself will improve things in this regard, but lightning information can help significantly more to solve these issues, as lightning provides indirect information about convective intensity and microphysics. The CST approach of separating the convective estimation from the stratiform estimation is well suited for the inclusion of lightning information since the lightning will help in only the estimation in areas of convective cores. For example, many IR-only techniques might have high false alarm rates in heavy precipitation due to overestimation in areas of stratiform precipitation under extremely cold cloud shields. This false alarm rate could be lowered significantly by certain techniques (still based on only IR), but at the great expense of detection possibility (e.g., Fig. 2 in Kuligowski 2002). However, the lightning information will help to solve this dilemma as shown in this study.
The general performance of the lightning-enhanced IR technique is comparable to previous studies in using IR-lightning or lightning-only data for rainfall estimation (
There are, of course, limitations in the lightningenhanced rainfall estimation technique developed in this study. First of all, lightning information in our technique is only connected to identifying convective areas and in estimating convective rain intensity, but lightning does occur in the stratiform or anvil cloud region of convective systems. However, lightning flashes initiating in stratiform or anvil areas are less than 10% of the total lightning (Dye et al. 2007;
In terms of practical application, the approaches described and tested in this paper seem well matched with the data from the forthcoming GOES-R satellite with its IR and lightning information. This approach, or others like it, seemed poised to take rain estimation from geosynchronous satellite information into a new era of improved results.
Acknowledgments. This research was supported by the
Adler, R. F., and
Arkin, P. A., 1979: The relationship between fractional coverage of high cloud and rainfall accumulations during GATE over the B-scale array. Mon. Wea. Rev., 107, 1382-1387.
Awaka, J., T.
Ba, M. B., and A.
Boccippio,D. J.,K. L.Cummins,H.
Carey, L. D., and
Christian, H. J., 1999: Optical detection of lightning from space. Proc. 11th Int. Conf. on Atmospheric Electricity (available as NASA Conference Publication 209261),
_____, 2008: The geostationary lightning mapper instrument for GOES-R. Extended Abstracts, Third Conf. on Meteorological Applications of Lightning Data,
Chronis, T. G.,
Dye, J. E., and Coauthors, 2007: Electric fields, cloud microphysics, and reflectivity in anvils of
Fierro, A. O.,
Gopalan, K., N.-Y. Wang,
Griffith, C. G.,
Gurka, J. J.,
Houze, R. A., 1989: Observed structure of mesoscale convective systems and implications for large-scale heating. Quart.
Jacobson, A. R., and
Joyce, R. J.,
Kozu, T., and Coauthors, 2001: Development of precipitation radar onboard the
Kuligowski, R. J., 2002: A self-calibrating real-time GOES rainfall algorithm for short-term rainfall estimates. J. Hydrometeor., 3, 112-130.
Kummerow, C., W.
_____, and Coauthors, 2001: The evolution of the Goddard profiling algorithm (GPROF) for rainfall estimation from passive microwave sensors.
Lang, T. J., and
MacGorman, D. R., and
McCaul, E. W.,
Morales, C., and
Negri, A. J., and
Petersen, W. A., and
Reap, R. M., and
Schmit, T. J.,
Schultz, C. J.,
Scofield, R. A., 1987: The NESDIS operational convective precipitation technique. Mon. Wea. Rev., 115, 1773-1792.
Seity, Y., S. Soula, and
Sheridan, S. C.,
Soriano, L. R.,
Soula, S., and S. Chauzy, 2001: Some aspects of the correlation between lightning and rain activities in severe storms.
Suszcynsky, D. M., and
Takahashi, T., 1978: Riming electrification as a charge generation mechanism in thunderstorms. J.
Vicente, G. A.,
Wiens, K. C.,
Williams, E. R., and S. Stanfill, 2002: The physical origin of the land-ocean contrast in lightning activity. C.R. Phys., 3, 1277-1292.
Xu, W., and
Zipser, E. J., 1977: Mesoscale and convective-scale downdrafts as distinct components of squall-line structure. Mon. Wea. Rev., 105, 1568-1589.
_____, 1994: Deep cumulonimbus cloud systems in the tropics with and without lightning. Mon. Wea. Rev., 122, 1837-1851.
Earth System Science Interdisciplinary Center,
Corresponding author address: Dr.
Most Popular Stories
- Major Phone Makers Sign Anti-Phone-Theft Pledge
- College Board Offers a Sneak Peak at New SAT
- 'Beige Book' Federal Reserve Survey, April 2014: Full Text
- Earthlike Planet Found in Red Dwarf's Goldilocks Zone
- Salsa Legend Cheo Feliciano Killed in Car Crash
- Twitter Offers App Install Ads
- Natural Gas Shoots Up on Bullish Stockpile Report
- Hiring and Weekly Jobless Claims Both Edge Up
- Chipotle Plans First Price Rise in 3 Years
- Generating Solar Power in the Dark