Abstract

Turbulent kinetic energy (TKE) budget variations during thunderstorm days (TD) and nonthunderstorm days (NTD) of premonsoon seasons of 2007, 2009, and 2010 have been investigated at a tropical station Kharagpur (22°30′N, 87°20′E) using the surface layer turbulence data obtained during severe thunderstorms-observations and regional modeling (STORM) experiment. Significant variations in the contributions of the TKE budget parameters with respect to stability are observed on these contrasting days of weather activity. In highly unstable conditions, smaller dissipation rates are seen on TD compared to NTD, while approaching near neutral conditions, higher dissipation rates are found in TD. New relationships between TKE dissipation rates with respect to atmospheric stability are proposed at Kharagpur for TD and NTD.

1. Introduction

Evolution of turbulence in atmospheric surface layer (ASL) assumes special importance, as it plays major role in the transportation of heat and moisture from near the surface to higher levels in the atmosphere. Turbulent kinetic energy (TKE), which is a measure of turbulence in the atmosphere, is directly related to the transport of momentum, heat, and moisture through the boundary layer [1]. Thus, understanding the variation of the individual budget components is crucial for energy exchange mechanism from atmospheric surface layer to upper atmosphere and vice versa. TKE budget has been studied during various field campaigns based on Monin-Obukhov similarity theory (MOST) [26]. TKE budget equation is still an uncertain relation in the quantitative case and is commonly used to parameterize turbulent properties of the surface layer in atmospheric models of larger scale [7]. TKE budget associates the local storage of turbulence to the shear production, buoyancy production, dissipation, and the transport processes and has numerous applications in both empirical and computational modeling in boundary-layer meteorology [810].

Study of TKE budget under unstable conditions and convective conditions is important for understanding the timely response of ASL. It was found in [11] that in the unstable surface layer, all terms in the TKE budget are considerable, and TKE generated through buoyancy forces is transported out of the layer, while the dissipation rate can be regarded as an approximation of the sum of mechanical production and the residual, with the unmeasured pressure term. Results at different sites about TKE budget and dissipation are emphasizing the importance of various TKE terms during convective situations and the necessity to be included in numerical weather prediction models, particularly in models that resolve mesoscale structures in storms (e.g., [1216]). Dissipation rate of the TKE is an important nondimensional parameter to study the kinetic energy evaluation in atmospheric boundary layer [17, 18] and is reported by researchers over the years in various parts of the world [1, 10, 18].

However, detailed studies of TKE budget variations are absent over Gangetic West Bengal region due to lack of surface layer observational network, where the premonsoon (March–May) thunderstorms occur every year [19]. To overcome the scarcity of available data over the region, a multi-institutional and international co-coordinated observational program known as severe thunderstorm-observations and regional modeling (STORM) has been launched by Department of Science and Technology (DST), Government of India [19]. Kharagpur (22°30′N, 87°20′E) is one of the experimental sites for STORM experiment, and data collected during premonsoon seasons of 2007, 2009, and 2010 as a part of this experiment are used in the present study. Highly convective conditions existing over this study region conduce to the occurrence of thunderstorms existing during premonsoon months over the study region. But these storms occur over a few days only. Hence, the present study aims to discern the variations and differences in the energy transport mechanisms in TD and NTD by studying the contributions of various physical forcings to total TKE during the premonsoon period over Kharagpur. This study is not looking into the dynamics of thunderstorm activity but rather looking into the turbulence transport during the whole day of thunderstorm.

2. Site Description, Data, and Quality Check of Data

The data for the present study are obtained from a 50 m instrumented tower, situated in an agriculture farm at the Indian Institute of Technology, Kharagpur region of West Midnapore, Gangetic West Bengal, India. The site is flat and grassy. The soil is sandy loam with a mixture of sand (64.1%), silt (20.1%), and clay (15.8%) and with a bulk density of 1.65 Mg m−3, volumetric heat capacity of  J m−2 K−1, field capacity of 26.7%, and wilting point of 9.3% [20, 21]. The mean sea level height of the station is 39 m. The data comprise air temperature (K), relative humidity (%), wind speed (ms−1), and wind direction (degrees) at six different levels (namely, 2, 4, 8, 16, 32, and 50 m) obtained from the 50 m micrometeorological tower. At 10 m level a sonic anemometer (Manufactured by R. M. Young) has been installed to measure the turbulent components of wind and temperature with 10 Hz frequency. For the present study, the fast response data averaged for 30 min are used. Complete description of sensors, with their manufacturer and model, is given in [22]. In the present study, both slow response data (SRD) of 1 Hz and fast response data (FRD) of 10 Hz (sonic data) of atmospheric variables for TD and NTD have been used for 2007, 2009, and 2010. Eight thunderstorm days (TD) and twelve clear nonthunderstorm days (NTD) are selected. The information about selected TD cases (time of thunderstorm, and associated rainfall) and NTD cases is given in Table 1.

Thunderstorm events over Kharagpur during premonsoon season are delineated using the log book information at the tower site and information obtained from cyclone detection radar at Kolkata, India Meteorological Department for the region. The classification of a thunderstorm day at Kharagpur is based on occurrence of a thunderstorm event at any time of that particular day [23] which says that if any thunderstorm event occurred within 24 hours of a day that day is termed as TD. And if there is no activity for whole day 24-hour period, that day is termed as NTD. Data sets have been subjected to quality check before analysis. Steps of quality checks employed in the present study are the same as those employed in [24]. For the present study, after employing steady state test in quality check of data, only those data sets which are satisfying the steady state assumption are used in analysis. This is essential as MOST can only be valid for the steady state. A total of 960 half-hourly data sets of FRD have been analysed for NTD (576) and TD (384), out of which 608 data sets have been finally considered into analysis (394 for NTD and 214 for TD) after quality check. The rejected data sets also include the rain hour’s data sets because of the possible effects of water droplets on transducer heads of sonic anemometers significantly influencing the data quality.

3. Methodology

Assuming a steady state and horizontally homogeneous conditions and choosing a coordinate system aligned to mean wind, the turbulent kinetic energy (TKE) budget equation can be written as (e.g., [3, 25]):

From the left hand side of (1), the first term is buoyant production term which acts as a source (sink) term during day- (night-) time. The second term is shear production term, which is always a source term. The third term is a divergence of vertical TKE flux, which can be a source or sink term depending on whether there is a flux convergence or divergence. The fourth term is a pressure correlation term that describes the redistribution of TKE by pressure perturbations. The fifth and last term is viscous dissipation of TKE, which is always a sink term, where are zonal, meridional, and vertical components of wind fluctuations; is the mean wind speed; and are turbulent parts of pressure and virtual potential temperature; is acceleration due to gravity; is log mean height; and are mean density and mean virtual potential temperature; is turbulent kinetic energy per unit mass; is dissipation; , , , and are vertical kinematic eddy fluxes of heat, momentum, energy, and pressure, respectively. will be computed as

During the pilot experiment at Kharagpur, the fast response turbulence measurements using sonic anemometer are available at 10 m height only, making it not possible to estimate the flux divergence following the methodology suggested in [3]. We make estimates of the magnitudes of flux divergence following [25] under the guidelines of Monin-Obukhov similarity theory (MOST). MOST is proven to be a good tool for understanding and explanation of the turbulent characteristics of lower atmosphere at any observational site. Buoyancy and shear terms are computed from slow response as well as fast response data.

Dissipation term is calculated from the high-frequency end of spectra obtained by available fast response turbulence measurements of sonic anemometer. The following equations are used in calculating dissipation (e.g., [26]): where is natural frequency (), is spectral density, , and are zonal, meridional, and vertical velocity components, is mean wind speed, and is a universal constant for the inertial subrange (0.55 chosen for and 0.73 for and ), following [27]. Pressure transport term of the budget equation (1) is estimated as a residual, which is likely to contain any possible errors in all the remaining terms of TKE budget equation.

The nondimensional form of dissipation rate of TKE is assumed to follow MOST and therefore can be expressed as a function of the stability parameter, [17]; that is, where is von Kármán constant (0.40), and is frictional velocity (ms−1).

4. Results and Discussion

4.1. Mean TKE Variation

Mean TKE () has been computed using (2) for all TD and NTD cases at the site. It is noticed that diurnal variation of is different during TD and NTD at the site. Diurnal variations of mean TKE at Kharagpur during 26 April 2007 (TD) and 6 May 2007 (NTD) are presented in Figure 1. Mean TKE values are generally higher on TD than NTD. The gap in the curve for TD is due to data removal by the quality checks. It is noticed that even after rainfall, there exists a high value peak in TD TKE curve at around 17:30 h, which may be attributed to the presence of high winds. Even in nighttime, value of TKE remains higher on TD (Figure 1). This shows presence of highly turbulent atmosphere on TD compared to NTD. Similar kind of variation has been noticed with higher values of TKE during individual TD cases compared to those of NTD at Kharagpur.

4.2. Variation of TKE Budget Terms

The contribution of different terms of TKE budget during TD as well as NTD is analysed. Diurnal variations of TKE budget parameters on one TD case (26 April 2007) and one NTD case (6 May 2007) are depicted in Figure 2, with general variations observed on all other TD and NTD discussed. The daily daytime (0600 to 1800 h LST) and nighttime (0000 to 0600 h LST and 1800 to 2330 h LST) averages of all TKE budget terms during each of the TD and NTD cases as well as respective cumulative averages are presented in Tables 2 and 3. The variation and significance of each term of (1) have been explained for the present study.

4.2.1. Shear Term

Higher value of shear production is observed on TD compared to NTD, as seen in Figure 2(a). Shear production values in the early morning and late evening compared to daytime are more for TD, but less for NTD. The reason for such high shear production values on TD can be attributed to high winds throughout the day. During NTD, until 06:30 h the winds are of the order of 1 to 2.5 ms−1 and then picked up to the order of 4.5 to 6.5 ms−1 during the rest of the day. This can be attributed to the higher shear production values from 06:30 h onwards in NTD. On the close examination of the shear term (Table 2), maximum shear production is observed on 26 and 27 April 2007. As expected, the shear term is the main source term of the TKE budget.

4.2.2. Buoyancy Term

The daily averages for 26 April (Figure 2(b)) as well as other days (Tables 2 and 3) reveal higher buoyancy generation during TD compared to NTD. It acts as source during daytime (which shows active thermal convection up into the atmosphere) and sink during nighttime (tends to suppress or consume TKE) to the generation of TKE, during TD as well as NTD cases. Negative nighttime buoyancy values can be explained on the basis of static stability present in the atmosphere. Day to day variability of buoyancy term is noticed during TD.

4.2.3. Flux Divergence Term

Flux divergence and buoyancy terms vary opposite to each other. Variation of this term for 26 April (TD) and 6 May (NTD) is shown in Figure 2(c). Daytime negative flux divergence means upward transportation of TKE from the surface layer, whereas positive flux divergence in nighttime indicates downward transportation of TKE into the surface layer. In general, the magnitude of the term is more during daytime as well as nighttime in TD compared to NTD except on 16 April, where high values of flux divergence have been observed. Table 2 shows more upward transportation of TKE in the TD cases compared to the NTD cases (Table 3). On 3 May 2007 (NTD), the value of flux divergence is very low, and daytime value is showing the downward transport of TKE while nighttime value shows upward transport.

4.2.4. Dissipation Term

Dissipation acts opposite to shear production during the present study. Dissipation values are higher on TD than NTD (Figure 2(d)) during both daytime and nighttime. Higher dissipation in TD varies in accordance with higher shear production in those days. Dissipation of TKE can be a significant source of heat in storms, and it increases the efficiency of the storm [16]. Existence of higher dissipation values on TD can be attributed to availability of high heat availability in ASL, making the atmosphere suitable for thunderstorm happening.

4.2.5. Pressure Transport Term

Negative pressure transport term with higher magnitude is seen during the entire TD than NTD (Figure 2(e)). Negative values of pressure transport term during daytime as well as nighttime were found during TD days whereas during NTD cases, they are positive during daytime and negative during nighttime except on 16 April 2007 (NTD) when the negative pressure transport term value is observed in daytime and during nighttime. It was found in [4] that pressure transport is greater than flux divergence term, and pressure transport is always negative. We have also got similar results during TD and with change of magnitudes in nighttime in case of NTD. It is difficult to interpret this result, since this term is estimated as a residual which possibly contains the errors of all the other terms of the budget equation. The analysis of all the terms of the TKE budget equation reveals that shear term is source for both daytime and nighttime; buoyancy term is source during daytime and sink during nighttime; flux divergence term is sink during daytime and source during nighttime; dissipation term is sink during both daytime and nighttime; and finally pressure transportation term is sink during both daytime and nighttime for TD and sink during daytime and source during nighttime for NTD.

As reported in [25], our results show buoyancy and flux divergence terms are of opposite signs but the magnitudes of both terms are not equal. It is seen from Tables 2 and 3 that the magnitude of flux divergence term is higher than that of buoyancy term during TD, but it is less during NTD, with opposite sign. This is one of the characteristic differences observed between TD and NTD cases. Higher dissipation values during TD provide more heat availability for the thunderstorm initiation [16].

4.3. Kinetic Energy Dissipation Rates

The nondimensional dissipation functions for TKE using eddy correlation method and following MOST are computed. Thus, can be expressed as a function of the stability parameter [17], where , is reference height, and is Obukhov length. It is suggested by various researchers [2, 3, 28] that in the nearly neutral conditions the mechanical production term equals the dissipation one. that is, is equal to 1. In [29] a relationship of with during unstable conditions was proposed

It was proposed in [28] that relationships based on Kansas results are as follows:

And for stable conditions,

Relationship of given by [29] has been widely used for unstable conditions and have been further modified by [30] for marine ASL. Relationships for both unstable and stable conditions provided in [30] were as follows: where for [30], and thus, the relationship for stable conditions comes out the same as (7). It was suggested in [17] that (9) can be written in a general form as

Equation (10) has been followed in [31] for . For very high stability, the is negligibly small and in [31], and 3.7 in [32]. Hence it was proposed in [17] that during stable conditions (nighttime), when the ground is colder than the air, TKE dissipation rate can adopt the -less stratification behavior, while under unstable conditions (during daytime) dissipation rates are not following any specific behavior.

For the present study, normalized dissipation rates of vertical wind with respect to for stable and unstable conditions during TD and NTD are estimated. Bin quartiles are computed for equal log spaced bins from data points to fit the empirical relationships in both TD and NTD. It has been found that the earlier proposed relationships for unstable or stable conditions are not satisfying the present study at both sites, and hence, there is a need for new relations in the present study for TD and NTD. Relations are different for TD and NTD at both sites, which is showing different dissipation during TD compared to NTD cases. Relations of dissipation rates at Kharagpur are based on relations of [28] and are as follows.

During NTD:

During TD:

These empirical relations, along with relations of EF and TG, with respect to , bin quartile points for unstable and stable cases during both TD and NTD, are shown in Figure 3. Interestingly, the dissipation rates of TD are lower than that of NTD in the unstable region of but higher than that of NTD for . These results are showing different behaviors of two dissipation rates in the free convection sublayer as proposed in [33]. However, in the stable region, the dissipation rates for NTD are higher than TD. The effect of highly convective atmosphere and higher production rates can be attributed to the higher dissipation during unstable conditions of TD.

4.4. Variation of Budget Parameters with Stability

Based on the stability , the period of the study is divided into three stability groups: near neutral , unstable , and stable during TD as well as NTD. Table 4 provides the information of the stability ranges chosen, TKE budget parameters during TD and NTD separated individually into different stability groups, and their averages computed under each stability group. All the TD and NTD cases (of 2007, 2009, and 2010) have been taken into account for dividing them into stability groups. The variation of TKE budget parameters with respect to stability is shown in Table 4 and discussed here under for Kharagpur.

4.4.1. Unstable Case

Mean stability is showing that TD cases are more unstable (−0.172) compared to NTD cases (−0.154). In unstable conditions, shear production, buoyancy, and pressure transport terms act as source terms, while flux divergence and dissipation terms act as sink terms in both TD and NTD. Flux divergence and dissipation together compensate the other terms of the budget equation. Magnitude of terms in TD is higher than that of NTD during unstable conditions. The buoyancy production is always a strong function of stability. Availability of sufficient moisture in the lower layers of the atmosphere and rapid increase of convection in terms of sensible heat flux play an important role in occurrence of the thunderstorm, and hence, more buoyancy production during TD compared to NTD can be explainable. Higher wind speeds during daytime in TD can be attributed to higher values of shear production. Flux divergence is showing that the transportation of TKE is more up in the atmosphere during NTD. Because of higher production, high values of dissipation in TD during unstable conditions are explainable.

4.4.2. Near-Neutral Case

Generally, near-neutral conditions are associated with overcast sky and strong winds resulting in more shear production in TD cases compared to NTD cases. Mean stability is almost the same for both TD and NTD cases (−0.005 for TD and −0.002 for NTD). In the near neutral stability, it is seen that shear production and flux divergence terms are source terms, while buoyancy and dissipation terms are sink terms during both TD and NTD. Pressure transport term is source term during TD and sink term during NTD. Magnitude of terms in TD is higher than that in NTD. During premonsoon season, higher sensible heat flux and high winds can be attributed to the higher order of buoyancy and shear production during both TD and NTD.

4.4.3. Stable Case

Mean stability values are higher during TD (0.187) compared to that of NTD (0.170) for stable conditions. For stable region, shear production and buoyancy terms have higher magnitudes in NTD compared to TD. As expected, buoyancy is negative and acts as sink term. In the stable case, buoyancy and dissipation together compensate all other terms.

Higher shear production, buoyancy, and dissipation term values for TD compared to NTD during unstable and neutral conditions are showing highly convective atmosphere supporting the thunderstorm activity. During stable conditions also, dissipation values are higher for TD, but shear production and buoyancy values are high for NTD. This shows the availability of more heat in the atmosphere during TD cases. Flux divergence values are almost of the same order.

4.5. Comparison with TKE Budget Variations during Monsoon Season

TKE budget terms variations for monsoon season over the same site have been studied in [14] using data sets recorded during 1990. Their results showed that dissipation term and flux divergence term together compensate other terms of the TKE budget equation during active as well as nonactive days of summer monsoon at Kharagpur. In the present study during premonsoon season, shear production term and buoyancy term together compensate the other terms of the budget equation. This variation is in agreement with the findings of [3, 14] which found that flux divergence term is significant close to ground and the term is getting insignificant while moving to higher heights. Higher values of flux divergence term have been found on the days of nonactive phase of monsoon. In the present study, the flux divergence term is significant during days of TD and negligibly small during NTD. For the monsoon season, it has been found in [14] that mechanically generated turbulence dominates the thermally induced turbulence in active phase. For the present study, it has been found that during TD, both mechanically and thermally induced turbulences are higher than those of NTD. These results are bringing out the major difference in the contribution of the production/sink terms in generating the TKE during summer monsoon and premonsoon season over Kharagpur.

5. Summary

The present study aims to discern the variations in the atmospheric surface layer turbulence transport and the processes that contribute to the total TKE during the days of thunderstorm and nonthunderstorm (no weather activity) within the same premonsoon period at Kharagpur. In the present study the hours during thunderstorm event are discarded because of not following steady state condition and hence MOST. The study is useful over the site as the transport of energy such as fluxes of heat, moisture, and momentum from surface layer to above is through turbulence mechanism only. A higher magnitude of TKE is noticed during TD than NTD. TKE budget analysis reveals buoyancy and shear production processes (mechanical and thermal forcing) together contributing to the total TKE in TD and NTD, but more in the former case. It is found that the transportation of TKE from the surface layer to upper layers of the atmosphere in TD is almost twice stronger than that in NTD. Highly unstable conditions were noticed in TD, while during nighttime, highly stable stratification was noticed in NTD. Shear production is main source term which compensates other budget terms along with buoyancy during unstable conditions and with flux divergence during near neutral and stable conditions in both TD and NTD. In TD higher dissipation is noticed. TKE dissipation rates for unstable and stable conditions on TD and NTD are different, with higher rates in TD during unstable conditions and in NTD during stable conditions. New relationships for surface layer dissipation rate with respect to stability for TD and NTD are proposed. Under stable conditions, dissipation rate presents a -less stratification behavior. This study reveals the difference in contribution of various TKE budget terms during monsoon (reported in the literature) and premonsoon period at the same site. During the monsoon season, mechanically generated turbulence dominates the thermally induced turbulence in active phase, whereas for the present study both mechanically and thermally induced turbulences are higher in TD than in NTD.

Acknowledgments

The authors gratefully acknowledge the Department of Science & Technology, Government of India, New Delhi, for providing funding for conducting a pilot experiment under STORM program during premonsoon months at Kharagpur. The authors thank Dr. M. Mandal for providing tower data of Kharagpur. They are also thankful to Mr. D. Pradhan, DDGM, India Meteorological Department, Kolkata, for providing the thunderstorm information based on radar imageries. The authors acknowledge with thanks the STORM Advisory Committee and its chairman Professor U. C. Mohanty, IIT Delhi, for providing necessary information about thunderstorm events.