Two seasons of Darwin, Australia, C-band polarimetric (CPOL) research radar, radiosoundings, and lightning data are examined to study the relative influence of the large-scale atmospheric regimes and the underlying surface types on tropical convective cloud properties and their diurnal evolution. The authors find that in the “deep westerly” regime, which corresponds to the monsoon period, the convective cloud occurrence rate is highest, consistent with its highest relative humidity. However, these convective clouds have relatively low cloud-top heights, smaller-than-average cell volumes, and are electrically least active. In this regime, the cloud cell volume does not vary significantly across different underlying surfaces and afternoon convective activity is suppressed. Thus, the picture emerging is that the convective cloud activity in the deep westerly regime is primarily regulated by the large-scale conditions. The remaining regimes (“easterly,” “shallow westerly,” and “moist easterly”) also demonstrate strong dependence on the large-scale forcing and a secondary dependence on the underlying surface type. The easterly regime has a small convective cloud occurrence rate and low cloud heights but higher lightning counts per convective cloud. The other two regimes have moderate convective cloud occurrence rates and larger cloud sizes. The easterly, shallow westerly, and moist easterly regimes exhibit a strong, clearly defined semidiurnal convective cloud occurrence pattern, with peaks in the early morning and afternoon periods. The cell onset times in these three regimes depend on the combination of local time and the underlying surface.
Convection patterns in the vicinity of Darwin, Australia, a site typical of the monsoon climate of northern Australia, have been investigated using ground remote sensing observations (e.g., Keenan and Carbone 1992; Rutledge et al. 1992; Williams et al. 1992; May et al. 2008; Protat et al. 2011). The main reasons for focusing on Darwin are that 1) the site has one of the most comprehensive long-term meteorological observational networks anywhere in the tropics; and 2) it experiences a wide variety of convective systems, and therefore should have important implications for the wider tropical Asia–Pacific region. Furthermore, the Darwin site combines seasonally varying meteorological conditions with distinct dry, wet, and transition seasons, with a complex topography of coastlines, islands, and oceanic areas. This makes Darwin an ideal location to investigate the relative roles of large-scale meteorology and surface types.
Past studies using data around Darwin explored the statistical characteristics of convection, where meteorological regimes were broadly separated into two categories; the buildup–break periods with low-level easterly winds and monsoon periods with low-level westerly winds (e.g., Keenan and Carbone 1992; Rutledge et al. 1992; Williams et al. 1992; May and Ballinger 2007). In break periods, cloud cells were reported to be more intense, taller, and electrically more active compared to monsoon periods. However, the daily total rain accumulation is higher during monsoon periods, with a ~50% contribution from stratiform rain (May et al. 2012).
Recent cluster analysis of thermodynamic sounding data using 49 wet seasons (defined as October–April) of radiosonde measurements showed that the Darwin wet season can be subdivided into five objective regimes, which have significantly different synoptic environments (Pope et al. 2009a). These regimes have been shown to be associated with significantly different properties of ice clouds for the Darwin region (Protat et al. 2011). Consequently it is worthwhile to investigate how convective cloud properties may change when data are separated into the five regimes instead of using the simple monsoon–break separation. Such a separation can aid the evaluation and development of convective parameterizations in models (e.g., Jakob 2003, 2010) as it can help better identify the relationship between the large-scale state (as defined by the cluster regimes) and small-scale cloud properties.
Several convective cloud properties will be considered in this study, including convective cloud occurrence and convective cloud-top heights, volume, kinematics, cell onset times, and electrical properties. Another important element affecting the growth of convective cloud systems is the merging of individual clouds since this leads to formation of larger cloud systems (Westcott 1994; Simpson et al. 1993). Previous research efforts in this area generally focused on a single convective cloud property. For example, Westcott (1994) considered case studies of convective cell merging and proposed that merging occurs because of horizontal expansion. Carbone et al. (2000) studied the Hector storms over Tiwi Islands and found that they formed mostly due to sea-breeze convergence. Pope et al. (2009b), using six wet seasons of satellite observations, found that in the north Australian region, mesoscale convective systems (MCSs) during the westerly (easterly) flow generally first formed over the western (eastern) side of Australia and then move across the continent. Building on these previous studies, a unified study of several convective cloud properties as to be carried out here will provide a more complete understanding of convective cloud properties and competing factors that regulate cloud growth.
May and Ballinger (2007) considered a small subset of aforementioned convective cloud properties for the Darwin region, which were identified using the automated Thunderstorm Identification, Tracking, Analysis and Nowcasting (TITAN) radar analysis tool (Dixon and Wiener 1993). Their convective radar echo-top height (ETH) statistics showed little evidence for a multimodal distribution as hypothesized by early observations (Johnson et al. 1999) and models (Liu and Moncrieff, 1998). Instead they found a continuous distribution of ETH with the peak of the distribution shifting toward the tropical tropopause layer (~15 km) as the distributions are conditioned on higher reflectivity (May and Ballinger 2007). But they did not provide any information on the variations with respect to underlying surface and local time or the variability with respect to recently identified large-scale atmospheric regimes.
The present paper aims to extend the May and Ballinger (2007) study. The specific objectives of this study are 1) to assess how the large-scale atmospheric regime affects the distribution of the convective clouds, ETH, and associated electrical activity by analyzing the diurnal and spatial variability; 2) to examine the variability of convective cell volume, kinematics, and cell onset times during the respective large-scale atmospheric regimes; and 3) attempt to ascertain the significance of the large-scale regime against other competing factors such as underlying surface and diurnal cycle, in the production of tropical convective clouds. This paper is organized as follows: the datasets, together with the techniques employed to extract convective cloud properties from radar reflectivities are described in section 2. The basic characteristics and spatial variability of convective clouds properties as a function of the large-scale atmospheric regimes is described in section 3a, followed by an analysis of the diurnal variability in section 3b. Finally, the results are summarized and discussed in section 4.
2. Datasets and method
The study makes use of two wet seasons (October 2005–April 2006 and October 2006–April 2007) of data from the Darwin C-band polarimetric research radar (CPOL; Keenan et al. 1998), the Australian GPATS (http://www.gpats.com.au/) lightning products and radiosoundings at Darwin airport. The sounding data are from the daily 2300 UTC (0830 LT) operational observations. The 2300 UTC data are selected to avoid modification of the environment by strong diurnal convection.
The CPOL radar (12.25°S, 131.04°E) as shown in Fig. 1, collects a three-dimensional volume of data out to a range of 150 km once every 10 min. Each volume consists of a series of 16 conical sweeps at elevations ranging from 0.5° to 42°. The radar transmits alternate linear horizontal and vertical polarization pulses of wavelength 5.3 cm. The main data source used in the present paper is the three-dimensional radar reflectivity after attenuation by rain is corrected for using the method developed by Bringi and Chandrasekar (2001). Other important polarimetric radar retrievals, such as drop size distribution and precipitating water contents are analyzed separately in a paper in preparation.
Figure 1 shows the extent of the domain sampled by the CPOL radar. Only data from the highlighted gray region (radar ranges of 20–120 km) are analyzed in this paper. This is done to reduce errors resulting from limited sample size at close ranges caused by the “cone of silence” occurring at elevation angles greater than 42° and at large ranges due to beam spreading. We also found that the mean radar ETH near maximum range of 150 km is ~1 km higher than the mean ETH within 120 km of radar center. The radar ETH statistics has a small range bias due to the radar scanning geometry; however, this effect is quite small within our radar sampling domain.
Reflectivity data are gridded by constructing a series of the constant altitude plan position indicator (CAPPI) at every 0.5 km in height (with a horizontal bin size of 2.5 km × 2.5 km) extending up to 20 km, using the Sorted Position Radar Interpolation (SPRINT) software. The gridded reflectivity data at a CAPPI level of altitude 2.5 km are processed using the “Steiner” convective–stratiform classification algorithm (Steiner et al. 1995) to determine the occurrence of the convective and stratiform precipitation at individual radar pixels. The Steiner algorithm classifies the gridded reflectivity as convective if the reflectivity value is at least 40 dBZ or greater than a fluctuating threshold depending on the area-averaged background reflectivity (within a radius of 11 km around the grid point). Each convective center has a radius of influence (ranging from 1 to 5 km) also depending on the surrounding background reflectivity (Steiner et al. 1995).
For each identified convective pixel at 2.5-km CAPPI level, the maximum height of the 5-dBZ echoes is computed to provide an estimate of the ETH. Specifically, the ETH corresponded to radar echo height whose reflectivity is the closest to 5 dBZ, but with a reflectivity value within the range of 0–10 dBZ, and provided there are continuous (in the vertical) reflectivity fields between the 2.5-km CAPPI level and this ETH. This procedure filtered out any possible effects of detached cloud layers situated above the convective towers. The 5-dBZ radar ETH definition has been previously used by May and Ballinger (2007).
In most cases, the true cloud-top height will extend higher than the 5-dBZ ETH; however, using CloudSat data the difference between cloud-top heights and radar 0- or 10-dBZ ETH has been found to often be within 2 km (Casey et al. 2012). Selecting the lowest available reflectivity per convective column might appear to be a better proxy of cloud-top height. However, this will introduce artifacts because the radar sensitivity drops with range, leading to fewer signals detected at longer ranges. The radar detection sensitivity is 0 dBZ near its maximum range of 150 km, so the choice of 5-dBZ threshold is sufficiently high to allow for detection of echoes at any radar range considered in this study.
This study also makes use of cell-based analysis, such as cell lifetime, speed, direction of movement, and volume. These parameters are derived using the TITAN radar analysis tool (Dixon and Wiener 1993). TITAN identifies convective cloud volumes based on radar reflectivity and volume thresholds. It then tracks these cloud volumes (hereafter referred as simply as “cells”) in space at discrete times (every 10 min in this case). Here a minimum volume requirement of 30 km3 and a reflectivity threshold of 35 dBZ are used to identify convective cells (e.g., May and Ballinger 2007). To reduce noise, filters are applied to the data. We only use information from cells that could be tracked over at least two consecutive radar scans. Thus, the analyzed cells had a minimum lifetime of 10 min. Moreover, only cells that formed and decayed within the radar sampling domain are used in the analysis. This is achieved by rejecting any track that passed beyond a 140-km radius (the maximum radar coverage radius is 150 km). Similar TITAN cell selection criteria have been used elsewhere (Goudenhoofdt et al. 2010). Overall, from a total of 50 485 cells that were detected by TITAN during the two seasons, these filters rejected ~56% of the cells, leaving just over 22 000 cells in our analysis. However, if one chooses to restrict the maximum radius to 120 km, as has been done for the Steiner method, 4500 more TITAN cells are discarded. Importantly, the TITAN analysis tool does not require gridded radar data, so the interpolation of the observed conical scans into CAPPIs is not a concern, which allows for an investigation of up to the 140-km range.
A crucial difference between the Steiner and TITAN methods of convective cell identification is that the former is likely to capture small cells such as those in the early growth or decay phase as well as mature cells, while TITAN has been designed to find mostly mature, intense cells. This is because the Steiner method does not have a minimum volume or lifetime requirement and permits lower reflectivities in the analysis (Steiner et al. 1995).
Finally, the electrical properties of the convective clouds are estimated using GPATS lightning data. Similar to other lightning detection networks, the GPATS network uses GPS-synchronized time stamps of the observed lightning sferics signals from each station and locates the strokes using the time of arrival method. To study the response of lightning associated with convective clouds, a lightning stroke was only used for the subsequent analysis provided there was at least one convective pixel occurring within a radius 10 km in distance and 10 min in time of this lightning stroke. These criteria rejected ~6% of strokes, from a total of 153 125 strokes detected within the radar domain over the two seasons.
In the subsequent analysis, the lightning occurrences are expressed in units of flashes per minute per pixel–cell. Several thousands of convective pixels had no lightning stroke associated with them and these “0” flash rates are retained during the calculation.
a. Basic convective cloud characteristics during the different large-scale atmospheric regimes
1) Mean regime characteristics
This section provides an account of the average cloud characteristics and associated electrical properties (Table 1), together with horizontal wind vectors, vertical shear of horizontal winds (hereafter, vertical wind shear), and vertical profile of humidity profiles (Fig. 2) for the five large-scale atmospheric regimes identified by Pope et al. (2009a). The long-term thermodynamic profiles and the large-scale environment are described in Pope et al. (2009a). Note that their details, but not their broadscale characteristics, will differ somewhat from our results because of interannual variability. The vertical wind shear profile complement results from Pope et al. (2009a), as wind shear has long been known to have an impact on convective organization, strength, and propagation properties (e.g., Rotunno et al. 1988).
Table 1 shows the 95% confidence interval range of the number of Steiner-identified convective pixels and TITAN cells, together with the lightning flash rate per pixel/cell as a function of the large-scale atmospheric regime. As TITAN keeps track of cell splits and mergers during successive radar volume scans, cell identification can be complex. Here, all cells that had the same “complex track identification number” are treated as one cell. The complex track number remains the same even if the cell splits or merges during its lifetime. For the study period, 78% of the detected TITAN tracks have a simple structure free of any splitting or merging events. The remaining 22% of cells have a complex structure, with a majority of them undergoing cell mergers. Merged cells are typically taller and larger than simple cells (Westcott 1994). The Steiner method treats each individual radar pixel independently.
Note that the 95% confidence interval of the convective occurrence frequency for the respective regimes does not overlap when using the Steiner method (Table 1). This is an initial indication that significantly different convective occurrence patterns do occur during the five large-scale regimes. Differences in convective cloud properties such as cell volume, propagation speed, and lifetime are significant when one compares the results of the deep westerly (DW) regime (corresponds to the active monsoon period) against that of the other regimes.
The dry easterly (DE) regime may be viewed as the trade wind regime. It mainly occurs in October and November (Table 1). The winds are southeasterly in this regime at low altitudes (Fig. 2a), reversing to westerly at ~8 km and back to easterly above a 16-km height. The upper-level (>15 km) easterly winds occur persistently in all regimes and are due to the presence of an upper-level jet. The DE regime has cells that typically lasted longer than other regimes. This could be due to the strongest low-level (0–3 km) and midlevel vertical wind shear (Fig. 2b). Robe and Emanuel (2001) and several earlier studies indeed suggested that strong shear in the lower levels produced more organized and longer lived convection. Both the Steiner and TITAN methods show that this regime has the lowest rate of convective activity; however, the lightning flash rate per convective pixel or cell is the highest. Low convective cloud activity is consistent with the lowest relative humidity (Fig. 2c), which is due to a dry continental air mass being advected over Darwin (Pope et al. 2009a). The existence of higher lightning flash rates during premonsoon (and monsoon break) conditions, than during the active monsoon period, has been previously documented over Darwin using lightning data from a separate lightning network (Höller et al. 2009; Labrador et al. 2009). Overall, the DE regime occurs ~11% of the time in our two-season sample and contains only very few detectable radar convective pixels (on average 312 pixels per day). We, therefore, choose not to show any further results from this regime from hereon.
The easterly (E) regime is typically seen as the transition between the trade wind regime and monsoon onset. It occurs mainly in the early and late part of the wet season. For the study period, this regime is the least frequent and accounts for only 7% of our total sample. The E regime has a higher average number of both convective pixels and cells than the DE regime, but is still smaller compared to the other regimes. The large-scale synoptic environment advects an air mass from the Coral Sea over Darwin (Pope et al. 2009a), which creates a moister environment than that of the DE regime (Fig. 2c). The horizontal wind vectors and vertical shear wind profile are similar to the DE regime except in the midtroposphere (8–15 km) where they are much weaker in the E regime (Figs. 2a,b). The lightning flash rate per pixel is moderately high in this regime and is consistent with premonsoonal lightning features (Höller et al. 2009; Labrador et al. 2009).
The DW regime is associated with typical monsoon conditions, and its occurrence peaks between January and March (Table 1). It accounts for 18% of the total sample. The large-scale synoptic environment indicates the presence of northwesterly winds at low levels (Fig. 2a) transporting an air mass of equatorial origin into the region (Pope et al. 2009a). They also found that this regime produced the highest amount of rainfall consistent with the highest relative humidity of all regimes (Fig. 2c). Both the Steiner and TITAN methods reveal that the DW regime generates the highest convective area and cell counts per day, respectively. However, the mean volume of the convective cells in the DW regime is relatively small, ~68 km3, compared to the other regimes with a mean cell volume of close to ~100 km3. This may be partly because ~90% of cells in this regime have a simple track structure, whereas the other two convective activity regimes [i.e., the shallow westerly (SW) and the moist easterly (ME) regimes described below] have only 67%–70% cells as simple. Also, the low and midlevel vertical wind shear is weakest in this regime, so the convection is predicted to be relatively short lived (e.g., Table 1) and less organized (e.g., Rotunno et al. 1988; Robe and Emanuel 2001). Consistent with previous studies, the DW regime is found to have the least amount of lightning discharges (Höller et al. 2009; Labrador et al. 2009).
The SW regime has previously been found to occur when the active monsoon region moves to the east of Darwin (Pope et al. 2009a). They found this regime to be associated with the largest mean convective available potential energy (CAPE) values of about 1100 J kg−1 and potentially stronger updrafts. The SW regime occurs 16% of the sample time. Table 1 shows that the SW regime has the second highest convective area per day and similar number of convective cells as the ME regime (described next). The SW regime is found to have the highest percentage of cells undergoing merger. Also, the electrical activity is consistently higher than in the other two frequently occurring regimes (DW and ME) regardless of the data processing procedure. The wind vectors in the SW regime change fairly rapidly in the first 2 km, veering from westerly near the surface to southerly at ~2 km (Fig. 2a). Between 2–8 km, the winds in the SW regime continued to be southerly and then strongly easterly above 15-km height. The rapid changes in the near-surface winds caused the low-level shear in the SW regime to be approximately 7 times more than in the DW regime (Fig. 2b). The relative humidity level is slightly less than that observed during the DW regime (Fig. 2c).
The ME regime can be viewed as the typical break monsoon period. This regime is the most frequent, occurring 48% of the sample time, and could be interpreted as the “default” state of the Darwin wet season. The convective area and cells numbers are similar to the SW regime but electrical activity seems to be slightly lower. The large-scale synoptic environment indicates the presence of easterly wind anomalies transporting an air mass of equatorial origin, together with large region of convergence over Darwin (Pope et al. 2009a). The sounding data (Fig. 2a) highlight the presence of easterly winds extending throughout the troposphere with the lowest wind magnitude near the ground level and at ~10 km. The low and midlevel vertical wind shear is moderately high in this regime, and therefore favors more organized convection compared to the DW regime.
2) Convective 5-dbz echo-top heights and associated lightning
This section shows the overall variation of the 5-dBZ ETH extracted using Steiner convective pixels and associated lightning as a function of the large-scale atmospheric regime. The diurnal features of these two cloud properties are in section 3b(2). Figure 3a shows the probability distribution function (PDFs) for ETH using 1-km bins in height for the four regimes with sufficient samples: E, DW, SW, and ME. The shaded gray region in this figure and all subsequent figures is the average distribution obtained using data from all days, regardless of regime classification. Figure 3b shows the vertical profile of the lightning flash rate per pixel.
The ETH distribution for all convective pixels (gray shaded region) shows a broad peak between 8 and 14 km (Fig. 3a). Each large-scale atmospheric regime shows a single peak occurrence in the ETH. There is no clear evidence of a multimodal distribution of convective ETH as reported by previous studies (e.g., Liu and Moncrieff 1998; Johnson et al. 1999) even though a significant amount of cumulus congestus cloud is present in our analysis. Our analysis is unable to reproduce the trimodal distribution of Johnson et al. (1999) because 1) each individual convective cloud could have several ETH, which will smear out the less-dominant peak occurring near the tropopause layer [~15 km; our main goal here is to study the convective fractions, so ETH data are considered more suitable than cloud-top height (CTH)]; and 2) the shallow cumulus clouds with peak heights within 1–2 km are usually missed because they are typically nonprecipitating and so cannot be captured with our C-band radar due to minimum detectable signal and sampling issues. However, there seems to be some evidence of a multimodal peak in convective ETH when it is presented as a function of diurnal cycle (see Fig. 11).
Figure 3a shows that the E regime (black) has the peak occurrence at the lowest height of all regimes (~8 km), followed by the DW regime (~11 km). The deepest convective clouds form in the SW regime with a peak occurrence at ~14 km and could be associated with stronger updrafts (e.g., Pope et al. 2009a). The mean distribution (gray shaded region) and the ME regime are mostly similar, since the ME regime is by far the most frequent (see Table 1). The TITAN method also produces the same dependence of ETH on large-scale regime except that the occurrence peak height is higher by 1–2 km (not shown).
The lightning occurrence profiles (Fig. 3b) show that the lightning rates increase strongly with convective ETH, with all but the DW regime showing the most lightning for the deepest clouds. The convective clouds in the E regime produce about 2–3 times more lightning than other regimes for all ETH up to a height of ~15 km. Notably, the SW regime has a secondary peak in lightning production rate associated with ETH around 10 km (Fig. 3b). In general, lightning is believed to be triggered when there is interaction between the upward flux of supercooled liquid water and the downward flux of graupel in the mixed phase (−10° and −40°C) of thunderstorms (Deierling et al. 2008). To maintain this process, sufficient CAPE to support vertical motions in excess of 6–7 m s−1 is required to supply supercooled liquid water in the mixed phase (van den Broeke et al. 2005). Large CAPE values potentially lead to stronger updrafts and higher ETH, so it is logical to expect the lightning flash rates to increase with ETH.
Figure 4 provides the spatial distribution of the average ETH, convective cloud occurrence frequency, and associated lightning flash rates. All data in this figure are interpolated to a 5 km × 5 km grid. We notice that the average ETH (top panels) is slightly higher beyond the ranges of 120 km (not shown) due to the beam spreading effect. Small size convective cells (which are typically shallow in height) with a narrow horizontal cross-sectional area become less frequent as the horizontal distance between adjacent beams widen at farther ranges because they are likely to be missed during the SPRINT interpolation. As a result mostly wider, taller cells contribute to the mean ETH near the maximum sampling range.
During the E regime (left column) a maximum in convective clouds is found over the ocean. The oceanic clouds in this regime generally have a higher mean ETH of ~10.5 km, compared to those occurring over land, whose average height is ~8 km. The lightning occurrence peaks tend to be collocated with regions of higher average ETH, with a significant proportion occurring along the coastline. Data from lightning networks have also shown significant lightning along the top end coastline of Darwin (e.g., Labrador et al. 2009).
During the DW regime the convective cloud occurrences are found to be larger over the western half of the domain, with the majority of them occurring in the Beagle Gulf (see Fig. 1 for location) and its coastal boundary regions. This region has been shown to have maximum precipitation during Darwin monsoon periods (May et al. 2012). The mean ETH is ~10 km, which is low compared to the other convectively active regimes (SW and ME), but convective cloud occurrence rate, especially over the ocean, is highest in this regime. This may be because during the DW regime, convection is embedded in a large-scale ascending region associated with the monsoon trough (May et al. 2012). The lightning locations are generally widespread and low in occurrence, with a maximum lightning occurrence being collocated with the maximum occurrence of convective cells.
A comparison of the spatial maps of the SW regime against the DW regime suggests that the peak convective occurrence locations show some tendency to shift eastward, from the western half in the DW regime to the region within 50 km surrounding the radar center. This is consistent with the conjecture that during the SW regime, the active monsoon region has moved to the east of Darwin (Pope et al. 2009a). The mean ETH is clearly the highest of all the regimes. The lightning occurrence rate is also the higher in this regime compared to the DW regime, with the maximum lightning occurrence located mainly over the ocean. Possible reasons for the higher lightning occurrence over the ocean than the land are discussed in section 3b. A closer examination of radar reflectivity loops and lightning occurrence reveals that the observed lightning occurrence peak is due to a significant number of events, not just a few extreme events.
During the most common regime (ME), maximum convective cloud occurrences are on the western part of the Tiwi Islands, consistent with the frequent occurrence of Hector storms (e.g., Carbone et al. 2000). Early storms typically occur over the eastern part of the Tiwi Islands and propagate westward during the break monsoon conditions. Carbone et al. (2000) explains that these storms intensify as they approach the west coast due to cell merger, and so more convective pixels are detected by the radar on the western part of Tiwi Islands. However, these Hector storms do not seem to be as electrically active as storms forming along the top end coastline. Focusing only over the Tiwi Islands, the lightning flash rate per convective pixel seems to the highest along the west coast region where the cell merger is most likely to occur. This is consistent with the electrical activity associated with typical Hector storms (Carey and Rutledge 2000). They found no significant lightning during the developing stage of the Hector storms and the maximum flash intensity was associated with the cell merger during the mature phase. The mean ETH shows moderate dependence on the underlying surface, with ETH slightly higher over the mainland than over ocean and Tiwi Islands. The convective activity is minimum northeast of Darwin in this ME regime.
Overall, the results shown in Figs. 3 and 4 indicate that convective cloud occurrence, ETH, and associated lightning depend both on the large-scale atmospheric conditions (as exemplified by the Pope et al. 2009a regimes) and the underlying surface. In the next sections, we investigate the effects of these two factors on other properties of convective cells such as cell lifetime, propagation parameters, volume, and cell genesis time. These cell properties are derived using the TITAN analysis tool.
3) Convective cell kinematics
The aim of this section is to examine the variation of convective cell kinematics (i.e., cell lifetime, speed, direction, and displacement) obtained using the TITAN tool and their spatial distribution in the four large-scale atmospheric regimes. Results for the DE regime are again not presented because on average this regime had two TITAN cell tracks per day. In all four regimes, the cells are mostly short lived with a mode occurrence lifetime of 20 min and a strongly positively skewed duration frequency (Fig. 5a). Longer-lived cells, such as those with lifetime exceeding 100 min (~5% of all TITAN cells), are found to be least frequent in the DW regime.
In contrast, the cell speed varied significantly during the respective regimes (Fig. 5b). The easterly regimes (E and ME) exhibit a much narrower distribution of cell speed with a peak occurrence near 3 m s−1. However, the westerly regimes (SW and DW) are characterized by a broader distribution, with 30% (15%) of the cells in the DW (SW) regime having cell speeds exceeding 10 m s−1. This greater cell speed in the westerly regime, particularly the DW regime, is because the steering flow speeds (wind speed at 700 hPa or ~3 km, see Fig. 2a) are larger in those regimes.
Figure 6 shows spatial maps of the cell track distribution and their average displacement, lifetime, and speed, as a function of regime. The cell displacement is calculated as follows. First, the coordinates of the cell center at first detection (t = 0 h) are grouped into 20 km × 20 km bins with respect to radar center. A 20 km × 20 km bin size is chosen to give at least 5 TITAN tracks per bin. Then for all cells in a bin the average location of the cell center at decay (t = termination of cell) is calculated. The average displacement vector is then defined as the position of cell decay relative to its onset and is shown as an arrow for each bin in the third panels of Fig. 6.
The spatial distribution of the TITAN tracks (top panels in Fig. 6) is similar to the distribution of convective pixels (second panels in Fig. 4). The most noticeable difference occurs in the ME regime, with the western part of Tiwi Islands showing comparatively less TITAN tracks than convective pixels. This can be explained since the TITAN occurrence maps show a given track only once at cell onset. As indicated above, the western island maximum found by the Steiner method represents Hector storms, which are usually born on the eastern part of the Tiwi Islands and then they propagate westward where the sea-breeze interaction makes them more intense (Carbone et al. 2000).
According to the bottom three panels in Fig. 6, cells tend to propagate for larger distances in regions located on the windward side of the incoming large-scale atmospheric circulation. For example, cells located in the northwest half of the domain in the DW regime and those in the southeast half in ME regime propagated for longer distances since they last longer and/or propagate faster. The steering flow mainly controls the direction of propagation of the TITAN cell, but it cannot explain the gradual drop in the cell propagation distance as they move from the windward side to the leeward side. The hypothesis that this gradual drop in cell propagation is an artifact because fast-moving and long-lived cells are more likely to be filtered out from the leeward side by our cell selection criteria (since they are more likely to propagate beyond 140 km from the radar center) was investigated and rejected. A similar result is obtained when we used all cells, even those that extended beyond 140 km from the radar center.
To further investigate this, we calculated spatial variation of the percentage of cells rejected compared to all cells when using a 140-km maximum radius requirement (provided the cell lifetime was at least 10 min). The results of this analysis are shown in the second panels of Fig. 6. It shows that our filters rejected less than 2% of cells in the circular region of radius 100 km bounded by the second concentric ring. Importantly this region does not show any spatial gradient in the cell-rejection frequency, but we still observed longer-propagating cells on the windward side compared to the leeward side in this inner region. Analysis of the spatial variation of the ratio of merged cells to all cells (results not shown) indicated that cell mergers on the windward side tend to be higher than on the leeward side. This result supports the observations of longer-lived cells (e.g., Westcott 1994). High-resolution 3D winds and gridded thermodynamic profiles for the region around Darwin would be needed to further understand the salient cloud physics causing this effect, which will be the subject of further investigations.
4) Convective cell volume
The aim of this section is to examine the variation of convective cell volume and its spatial distribution with a large-scale atmospheric regime.
The SW and ME regime show a similar distribution of TITAN cell volumes (Fig. 7), with both the DW and E regimes deviating from the mean distribution more significantly. The proportion of cells with a small volume of 30 km3 is ~15% for the SW and ME regimes, while it is much larger (~22%) for the DW regime and smaller (~11%) for the E regime. Bigger volume cells are most frequent in the E regime, though results are drawn from a smaller number of events. Within the convectively active regimes (DW, SW, and ME), cells with a large volume are more frequent in the SW regime (55% of the cells had volume >60 km3) and ME (51%) regimes compared to the DW regime (37%). An interesting feature of the SW and ME regimes is that the cells over land have a larger volume compared to those occurring over ocean (Fig. 8). In contrast, in the DW regime the cell volume shows little dependency on the underlying surface. The drop in cell volume at the far southeast of Darwin could be an artifact associated with increase in the rejection of TITAN cells by our filters (second panels Fig. 6). Overall, this points out that the convective clouds in the DW regime are embedded within the large-scale monsoon trough.
Overall, the variability in cell volume is linked to both the large-scale atmospheric circulations and the nature of the underlying surface. For example, cell volume is largest in the E regime, smallest in the DW regime, and intermediate in the SW and ME regimes. Comparing the three most frequent regimes, they all, except for the DW regime, have larger cells over the continent than over the ocean. Since cell volume (Fig. 8) reveals a similar response as the cell area (results not shown here) and to some extend as the ETH (top panels Fig. 4), it is fair to assume that cells with larger volume will have a greater mean ETH and a wider horizontal extent.
b. Effects of the large-scale regime on the diurnal cycle of convection
Having identified significant differences in basic cloud cell characteristics for the four large-scale regimes used in this study, this section focuses on the diurnal cycle of cell characteristics, in particular convective ETH occurrence and associated lightning, as they are indicative of the intensity and microphysical characteristics of the convective systems.
1) Convective cell onset time
In this section we examine variation in cell onset time by binning the onset times with respect to the Darwin local time (LT = UTC + 9 h, 30 min). The distribution of cell onset times (Fig. 9) shows that most of the cells are triggered during the day with a secondary peak occurring in the early morning period. For the DW regime, the daytime peak of the cell onset occurs the earliest, around midday, followed by the SW regime at 1400 LT and around 1500 LT for the easterly regimes.
The spatial maps of the cell onset times (Fig. 10) show well-defined differences in the dominant local time of the onset of convective cells with respect to the underlying surface. Some caution must be exercised when interpreting the results shown in Fig. 10, as the colors only represent the modal local time of the onset of convective cell development. Obviously some cells will be born outside the modal local time period for a given underlying surface. Over the ocean, the cells are triggered mainly in the early morning and in some cases around midnight, regardless of the regime type. Over land, the cells are predominantly triggered in the afternoon except for the DW regime. In the DW regime, the triggering of the cells within ~60 km from the coastline happens around midday, while for the remaining land region it still occurs in the afternoon. These features in the diurnal cycle of cell onset time with respect to different underlying surface types are consistent with earlier research (Liu and Zipser 2008 and references therein).
In all regimes except the DW regime, convective cells over land are likely initiated by sea breezes whereas ocean cells are predominantly triggered by the land breeze. Thus, the cell onset times are strongly dependent on diurnal cycle and on the underlying surface, in at least three out of four regimes. In contrast, in the DW regime (or monsoon period) with extensive cloud cover, radiative heating of the land is less effective resulting in changes to the mechanisms that trigger convection (May et al. 2012).
2) The diurnal cycle of convective ETH
Figure 11 shows the evolution of ETH occurrence frequency as a function of time of day and height for each of the large-scale regimes. The ETH occurrences are calculated separately for each bin of 1 h in local time and 1 km in height, and then normalized by the number of days in each regime. For clarity, the counts are then further divided by the peak occurrence value in each panel (peak values given on the bottom right-hand corner). The density of points as a percentage of the maximum occurrence is presented using a color scale with white indicating that no data is recorded in this bin.
In the E regime, the convective echo occurrence is highest in the afternoon and in the early morning period (Fig. 11). It appears that, especially in the afternoon period, the clouds are generally shallow during the early growth phase and progressively develop into deeper clouds in the mature stage. This diurnal cycle is consistent with that of the nonprecipitating ice clouds over Darwin during that same regime, as characterized in Protat et al. (2011). This consistency suggests that in the E regime, nonprecipitating ice clouds are predominantly convectively generated. At all times, except for the afternoon period, mean ETH (black curve) during the E regime is lower than the mean values for all regimes (black–white dashed curve). The electrical activity in the E regime is semidiurnal and follows the convective echo occurrence frequency, with the lightning flash rate peaks occurring fewer hours prior to peaks in convective ETH occurrence (white curve).
The DW regime shows a prolonged period of occurrence of convective clouds from midnight through the morning with a peak around midday, and a clear occurrence minimum in the evening (Fig. 11). Typically during monsoon conditions, which the DW regime represents, there is a large proportion of stratiform clouds (May and Ballinger 2007). Hence, the convective ETH diurnal cycle is expected to deviate from that of rainfall, which often shows a maximum in the afternoon and evening. Overall, the average ETH of ~10 km is generally lower than in the all-regime average. The DW regime is the least active in terms of lightning and this could be due to insufficient updraft speeds within the convective core to produce lightning (e.g., van den Broeke et al. 2005). Unlike the E regime, the frequency of occurrence of nonprecipitating ice clouds in the DW regime in Protat et al. (2011) is very different from the convective ETH statistics obtained here. The maximum in nonprecipitating ice cloud occurrence occurs later than the convective ETH occurrence maximum, between 1500–2000 LT (Fig. 2d in Protat et al. 2011). This comparison suggests that during the DW regime, thick nonprecipitating anvils and cirrus decks produced by deep convection are much longer lived than during other regimes. During the DW regime, the diurnal variation in atmospheric temperature is weak due to widespread cloud cover reducing the daytime heating of the land (May et al. 2012). This largely explains the lack of a strong evening peak in the occurrence of convection during this regime.
During the SW regime the average ETH is higher than the mean values for all regimes at all times of the diurnal cycle, with two peaks: one in the morning and one in the afternoon. We previously have shown that the SW regime also contains the tallest convective ETH (Fig. 3a) and with moderate cell volume (Fig. 7) possibly due to stronger updrafts and increased occurrence of cell merging. The peak in nonprecipitating ice cloud occurrence (Protat et al. 2011) is shifted to a later time (2000–2400 LT), suggesting again the production of extended anvils by deep convection associated with the SW regime, as is the case for the DW regime as well. The SW regime is found to have the second highest lightning activity, with the majority of lightning strokes generated by the early morning storms. Again the peak in lightning flash rates tends to occur few hours ahead of the peak in convective ETH occurrence.
During the most frequent ME regime, the results reveal that the early phase of storm development occurs at ~1500 LT with a peak height of 9 km (Fig. 11). These cells mature within a few hours, becoming towering cumulonimbus clouds with a peak occurrence height of 14 km. This diurnal cycle is consistent with that of the nonprecipitating ice clouds (Protat et al. 2011) in this regime. This suggests that thick anvils and cirrus decks produced by deep convection are shorter lived than during the DW and SW regimes. From the evening through the night the convective systems gradually decay causing a gradual drop in ETH. This drop is also found in the nonprecipitating ice cloud statistics (Protat et al. 2011). The infrared satellite observations analyzed by Pope et al. (2009b) confirm a mesoscale convective system genesis time near 1500 LT during monsoon break periods, and these usually decay within approximately 3 h. The lightning flash rates are highest when the ETH reached the peak heights in the evening period.
A shortfall in Fig. 11 is that the responses in ETH may be affected by the complex topographic environment around Darwin. We attempt to rectify this by further splitting the time–height pdfs of ETH into three groups of different underlying surface types, namely, oceanic, coastal, and continental. The results are shown in Fig. 12 and the area covered by the three surface types is shown in Fig. 1. In the DW regime (Fig. 12a), the peak occurrence in convective clouds occurs earliest over the oceanic surface in the early morning period and progressively shifts inland, peaking over the continental surface near midday. This progression of convective cloud activity from the oceanic region through coastal and then over land is consistent with the picture that convection in the DW regime is embedded in the large-scale forcing by the monsoon trough. In contrast, convection in the ME regimes seems to be primarily dependent on conditioning of the atmosphere by land and sea-breeze processes. For example, the majority of convective cloud activity occurs above the oceanic region in the early morning period, with peak heights at 14 km. In contrast, during the afternoon and evening periods, the convective cloud occurrence is highest above the coastal (bimodal peak height of 9 and 14 km) and continent (peak height of 14 km) regions, respectively. There is little evidence to suggest that the early born coastal convection is progressing over the continent since storms in the ME regime mainly propagate toward the ocean (see Fig. 6). Results for the E and SW regimes are not shown because they exhibited less noticeable differences in the convective cloud occurrence over the three underlying surfaces.
Overall, the results shown in Figs. 11 and 12 indicate that the diurnal cycle of convective cloud occurrence and their top heights, and spatial location of the convective clouds, contrast considerably among the four large-scale atmospheric regimes. First, all the regimes, except the DW regime, show intense convective activity in the late afternoon, presumably initiated by the sea-breeze circulation that forms on the top end coastline. The results also indicate that sea-breeze effects are less important during the DW regime. Second, the DW regime clearly shows oceanic characteristics, while the ME regime demonstrates much more continental characteristics. Third, the SW regime (and the E regime, though results are drawn from a smaller number of events) show high convective activity after midnight and in the early morning, thus showing that convection in this regime exhibits somewhat oceanic characteristics. Finally, the comparison of the diurnal cycle nonprecipitating ice clouds and convective cloud towers indicate that the thick anvils and cirrus decks produced by deep convection are shorter lived during easterly regimes (E, ME) and longer lived during the westerly regimes (DW, SW).
Higher lightning flash rates after midnight (vs afternoon or evening), particularly in the E and SW regime, and over the coastal boundary region (vs continent) do not seem to be consistent with the traditional picture of having more lightning over land and in the afternoon period. The complex topography of coastlines, islands, and oceanic areas within our sample area combined with the distinct wet regimes may be contributing toward this discrepancy. On the other hand, since the Darwin site with its Doppler radar pair can provide higher-resolution 3D wind data, it will offer an opportunity in the future to derive upward mass fluxes and to check consistency with lightning activity (e.g., Deierling et al. 2008). The question is as follows: for a given mass flux rate, do convective cells produce more lightning when located over land (vs sea) or in the afternoon period (vs the early morning period)? This is the subject of ongoing investigations.
4. Conclusions and summary
Polarimetric weather radar data collected over two wet seasons (October 2005–April 2006; October 2006–April 2007) at the tropical low-latitude station of Darwin, northern Australia, are used to study the variability of convective cloud properties with both the large-scale state of the atmosphere, the diurnal cycle, and the underlying surface type. The properties of convection studied here include the frequency of convective cloud occurrence, 5-dBZ echo-top heights (ETH), kinematics (lifetime, speed, and direction of propagation), cell structures, and volumes. Both the spatial and diurnal variability of these tropical convective cloud properties are studied as a function of the identified main large-scale atmospheric states in this area.
A summary of the key findings is as follows:
The most frequent ME (break) regime shows the highest convective activity from afternoon to midnight and a secondary occurrence peak in the early morning. These convective clouds occur most frequently on the western part of Tiwi Islands, which is consistent with the signature of the well-known Hector storms. In the afternoon the convective clouds are initially shallow with a modal height of ~9 km, and within a few hours grow into deeper convective towers with a modal height of ~14 km. The ETHs are higher and the cloud cell volumes are larger over land than sea. It is also very clear from the results that the land cells in this regime are predominantly initiated in the afternoon by sea-breeze processes whereas ocean cells pop up in the early morning due to land breeze effects. Overall, the convection in the ME regimes seems to be well organized and shows characteristics similar to continental convection. Since this regime occurred for nearly 48% of the wet season, its convection patterns could be a fair representation of the default climatology of Darwin.
In contrast, the DW regime, which corresponds to the active monsoon period, exhibits the highest overall probability of generating convective cells. It has a peak convective cloud occurrence over the coastal boundary region from midnight to early afternoon. The evening convective activity is least frequent in this regime and is thought to be due to the presence of continuous cloud cover reducing daytime heating that prevents the establishment of sea-breeze convergence. The vertical wind shear in the low levels, convective ETH, cloud cell volumes, and lightning activity are all smaller in this regime compared to the other convectively activity regimes (SW and ME regimes). Also, the effect of the underlying surface types on most convective cloud properties is the weakest in the DW regime. Overall, clouds in this regime exhibit oceanic characteristics, with convection being embedded in the large-scale forcing of the monsoon trough.
In the SW regime, the peak convective occurrence location shifts eastward compared to the DW regime. This observation supports the hypothesis that these two regimes are connected to the eastward propagation of the monsoon trough. Another feature in the SW regime that matches with the DW regime is the increase in occurrence of convective clouds in the early morning period. However, unlike the DW regime, the effect of the underlying surface on the convective cloud properties is somewhat strong in the SW regime. For example, the land cells predominantly initiate in the afternoon and have a larger volume compared to those that form in the early morning over the ocean. Another contrasting feature is that the convective cloud activity in the SW regime is moderately high in the afternoon. Overall, this indicates that the SW regime are regulated by a mixture of large-scale forcing that are important for the DW regime and the sea-breeze effects that dominate the ME regime.
The E regime behaves in a similar manner to the SW regime. Like the DW regime, the E regime has the highest convective cloud activity in the early morning period. While the observed secondary peak in convective cloud activity in the evening period can be attributed to the sea-breeze effects, the effect of the underlying surface on the convective cloud properties is moderate. Contrary to the SW and ME regimes, the E regime has somewhat higher ETHs and larger cloud volumes over ocean than land. The convective clouds in this regime have one of the highest tendencies of producing lightning flashes, and most of these electrically active clouds are located at the top end of the Darwin coastline.
The main purpose of the study was to use the complex meteorological and topographic environment around Darwin to study the relative influence of the large-scale atmospheric conditions, as represented by a set of synoptic regimes, and the underlying surface types on the basic characteristics of convective systems and their diurnal evolution. The picture emerging from this study shows an intricate interplay between the large-scale regime and surface-type influences on the properties of convection. To first order, the large-scale regime determines much of the convective evolution, as exemplified by the rare occurrence of convection in the E and DE regimes, and the widespread occurrence of relatively weak convection in the DW regime. However, complex topography, such as the presence of coastlines, is a major secondary factor in determining the structural characteristics of convection. For example, during the ME regime, much of the convection is triggered along sea-breeze fronts either over the Tiwi Islands or the mainland. This indicates that the large-scale state does not allow convection to occur spontaneously over the ocean, but does allow for more organized forms of convection. This picture is likely typical not only for the north of Australia, but the entire Maritime Continent, where the existence of numerous islands of varying size can trigger sea-breeze convection even in large-scale conditions unfavorable for widespread convection over oceanic areas. In contrast, during the DW regime, the surface influence becomes negligible, as the large-scale upward motion associated with the monsoon trough provides sufficient forcing to allow widespread convection with large areas of long-lived stratiform cloud, which in turn suppresses the daytime heating of the land.
This work has been supported by the U.S. Department of Energy ARM. We would like to acknowledge the contributions of Brad Atkinson and Michael Whimpey in supporting the Darwin observatory and data management. Rodney Potts and Kevin Cheong are thanked for providing the TITAN data and for their discussion about it. James Sofra is thanked for GPATs lightning data. Special thanks to Susan Rennie and Surendra Rauniyar for their discussion and comments.
The Centre for Australian Weather and Climate Research is a partnership between the Bureau of Meteorology and the Commonwealth Scientific and Industrial Research Organisation.