In recent years, extreme weather events such as heavy rainfall are common and have contributed to loss of lives, damage of properties, destruction of environment and other social economic livelihood of people in many countries   . Tanzania has witnessed many heavy rainfall events in recent years. Table 1 depicts some of the more recent heavy rainfall events which have occurred over different regions in Tanzania. Forecasting of these events is of considerable benefit to different sectors (agriculture, hydropower, health, water resources) and community at large. However, it is difficult to forecast these events accurately and reliably   . In Tanzania, the difficulty to forecast heavy rainfall is associated with the country’s complex topographical landscapes, numerous large inland water bodies, variation in vegetation types and land-ocean contrast  .
Several studies  -  have analysed the drivers of heavy rainfall events over East African region in general and Tanzania in particular.  analyzed the link between El Nino events in the equatorial Indian Ocean and the catastrophic rainfall of 1961-1962 in East Africa. They found that, heavy rainfall over East African region is associated with large Sea Surface Temperatures (SSTs) and wind anomalies over equatorial Indian Ocean.  analysed the atmospheric circulation patterns over the Indian Ocean that was linked to heavy rainfall event that occurred over northern Tanzania during the short rain of October-November-December (OND) season in 2006. They found that heavy rainfall was associated with strong warming over the Indian Ocean, coupled with convective zones over the western Indian Ocean and East African region.
The above mentioned studies have analysed the impacts of large scale atmospheric circulation, atmosphere-ocean interactions and SSTs on heavy rainfall events. However, relatively less attention has been directed towards analysing the impact of local scale induced circulations from orography, land ocean contrast and difference in vegetation
Table 1. Regions received heavy rainfall from 2011to 2015 in Tanzania.
type on formation of heavy rainfall events. Furthermore the combined impact of local and large scale circulations to the atmospheric environment at which was the main cause of heavy rainfall events have not been analysed. In this study, we modify the convective vorticity vector () which is defined as a cross product of absolute vorticity and gradient of equivalent potential temperature to moist potential vorticity vector by replacing the equivalent potential temperature with the moist-air entropy potential temperature (). The usefulness of in diagnosis of heavy rainfall events is analysed. The ultimate goal is to highlight the possibilities of using to help prediction of heavy rainfall events over different regions in Tanzania.
Motivation of the Use of New Moist Potential Vorticity Vector to Describe Rainfall Events
Heavy rainfall events, among other factors are triggered by transport and convection of large volume of water vapour. Physical variable such as moist potential vorticity (MPV), that describe both dynamics and thermodynamics aspect of water vapour can be used to diagnose and potentially predict the occurrence and distribution of heavy rainfall  . Indeed, MPV combines the two general characteristics responsible for the rainfall formation: the dynamics via the vorticity component and the thermodynamics via the gradient of the moist-air temperature.  used MPV as a fundamental dynamical variable to examine a four dimensional dataset obtained from an 18-hr high resolution 25 km simulation of a squall line that occurred during 1200UCT 10 June to 0600UCT 11 June 1985. They found that the stratiform region is characterised by negative MPV and suggest that trailing stratiform precipitation could be viewed as an end of slantwise convection with pronounced moist symmetric instability (MSI).  used MPV as a diagnostic variable of heavy rainfall event that occurred in MEIYU and found that MPV is an important physical variable to reveal heavy rainfall structure and dynamic mechanisms. They also found that negative MPV correspond to the MEIYU front-wind shear line system and the negative center corresponds to the heavy rain center.
However, the scalar MPV which is derived from the dot product of absolute vorticity and gradient of moist-air temperature may not fully explain the complex atmospheric circulations patterns over the tropics. The Coriolis parameter over the tropics is very small and is zero over the equator. Furthermore the vertical gradient of temperature over the tropics is small due to strong mixing processes. Thus the scalar MPV may not represents some dynamical patterns over the tropics, particularly close to equatorial region. For instance,  demonstrated that the vertical component of planetary vorticity in MPV computation is zero and cannot explain two dimension (2-D) equatorial atmospheric flows.  argued that over the tropics, the moist potential vorticity computed from the dot product is mainly in the horizontal direction since the absolute vorticity and the gradient of temperature are small and cannot represent the atmospheric motions reasonably.  proposed the use of convective vorticity vector () to study deep convection particularly over tropical regions. However, is computed using the gradient of equivalent potential temperature. It is conservative only in moist adiabatic and frictionless processes and its generalization to be used in dry atmosphere lead to annihilation of solenoidal term in its tendency equation  . Since the atmosphere is never completely dry or saturated, but non-uniformly saturated. Thus cannot be used to study non-uniform saturated atmospheric flow and fulfil the demand to verify, at the same time, a moist and dry air conservative property and an invertibility principle. In this study we modify to moist potential vorticity vector by replacing the equivalent potential temperature with moist-air entropy potential temperature (). This temperature is valid for a general mixing of dry air, water vapour and all possible condensed water species. is analysed on its ability to diagnose heavy rainfall events over different regions in Tanzania.
2. Data and Methodology
2.1. Model Description Experimental Design and Data from the Model
The Weather Research and Forecasting (WRF) model version 3.3.1, jointly developed by the National Oceanic and Atmospheric Administration (NOAA) and National Centre for Atmospheric Research (NCAR) is used in this study. This model is a non-hy- drostatic mesoscale NWP model, fully compressible and has terrain following sigma coordinates   . The WRF is chosen in this study because it has been used by different researchers at government institution and universities and is currently running at Tanzania Meteorological Agency (TMA) Numerical Weather Prediction (NWP) modelling section. It features multiple dynamical cores, a 3-dimensional variation and 4-dimensional (3DVAR and 4DVAR) data assimilation system, and has software architecture that allow for computational parallelism. WRF is suitable for a broad spectrum of applications across scales ranging from meters to thousands of kilometres.
In this study the domain of WRF is set-up over East Africa region (8˚N to 18˚S; 25˚E to 52˚E) (Figure 1), on a horizontal resolution of 15 km, and the model has 28 vertical levels. The model was set to simulate heavy rainfall events that occurred over different regions in Tanzania from 20th to 22nd December 2011 and 5th to 8th May 2015. The boundary condition from the National Centre for Environmental Prediction (NCEP) global forecast system (GFS) dataset with a 0.5˚ × 0.5˚ spatial resolution were used to force the WRF. The input data has a time interval of three hours, while the output data are served at hourly interval. The simulations were initiated at 00h UTC and run for 48 hours ahead. Simulated meteorological variables: wind speed (zonal, meridional and vertical components), pressure, temperature and relative humidity at different pressure levels are used to compute.
2.2. Rainfall Data
Rainfall as cumulated total for every 3 hours was accessed from the Tropical Rainfall Measuring Mission (TRMM). This is an international project led by the National Aeronautics and Space Administration (NASA) and Japan Aerospace Exploration Agency (JAXA). The project provides improved estimate of rainfall over tropical regions using remote sensing technique through satellite observations. The TRMM rainfall
Figure 1. The topographical map in the domain where WRF is set indicating surface land cover that can be resolved by the model in the study region.
estimates used here are from TRMM 3B42V7, available at spatial resolution of 0.25˚ by 0.25˚ which is approximately 27.8 km by 27.8 km. These data were accessed from the website:
http://giovanni.sci.gsfc.nasa.gov/giovanni/#service=TmAvMp&starttime=&endti and was used to map the spatial distribution of rainfall (mm/day) during 20th to 22nd December 2011 and 5th to 8th May 2015. These maps are compared with the spatial patterns of.
2.3. Observed Rainfall Data
Daily accumulated rainfall (mm/day) during 20th to 22nd December 2011 and 5th to 8th May 2015 were acquired from the Tanzania Meteorological Agency (TMA). In order to compare the spatial patterns of MPVV with rainfall, the Inverse Distance Weighting (IDW) interpolation technique was used to interpolate the patterns of observed rainfall using data from 22 weather stations. For detailed description about the IDW interpolation technique a reader may consult  . The spatial rainfall maps generated using IDW interpolation techniques are compared with spatial patterns of.
2.4. Descriptions of Rainfall Events
On 20st to 22nd December 2011, there was very heavy rainfall over the city of Dar es Salaam which triggered a worst flooding event where forty three people were reported dead and many were missing. The major transportation networks of the city were destroyed by the flood. Some bridges were completely swept out and others were covered with water. The heavy rainfall event were caused by warming over the eastern, North and northwestern Australia and the prevailed surface Westerly wind, that pushed warm surface water towards western Indian Ocean. This warm surface water spread around western equatorial Indian Ocean coupled with significant warming developed over western Indian Ocean (coast of Dar es Salaam) on 20-22 December enhanced more convection along the coastal regions. Moreover, low level convergence dominated over most regions that enhanced heavy rainfall over much parts of the country. Figure 2 indicates the cloud band over the coastal region detected from satellite image on 21-12- 2011, at 0900UCT.
On 6th May 2015, the city of Dar es Salaam received heavy rainfall of 111.4 mm in 24-hours. Kibaha region received 88.3 mm in 24-hours, while Zanzibar reported 45.5 mm in 24-hours. On the 7th May 2015, rainfall of 60.2 mm, 52.5 mm, 49.7 mm, 41.4 mm and 37.4 mm were reported at Mahenge, Dar es Salaam, Kibaha, Zanzibar and
Figure 2. Satellite image (a) visible band; (b) infra- red band on 21-12-2011, at 0900UCT.
Mtwara regions respectively. These rainfalls were influenced by the enhancement of the inter-tropical convergence zone (ITCZ) over Tanzania, and significant moisture flux over the coastal belt coupled with low level convergence that triggered deep convection over most regions.
2.5. Modification of CVV to MPVV
We start by defining the moist-air entropy potential temperature. The moist-air entropy potential temperature is defined in  as
where is a key quantity. It is mentioned in  that de-
pends on the standard entropies of water vapour and dry air and. It is also mentioned that is a good approximation of. For detailed derivation of please refer to  -  .
Another parameter to use in computation of is density which is defined from the equation of state as
where p is atmospheric pressure (in Pa) at different level, is atmospheric pressure at reference level, is potential temperature, is specific gas constant for
dry air and, is specific heat capacity at constant pressure. The last pa-
rameter used for computing is the absolute vorticity which is defined as
where is the relative vorticity defined as, where a is
the radius of the earth and is the latitude, f is the coriolis parameter defined as.
is defined as
Considering the hydrostatic equilibrium, Equation (5) can be re- written as
First component, (6)
Second component, (7)
The third component, (8)
The magnitude of is written as
3. Results and Discussion
In this study we analyse whether can serve as diagnostic variable of heavy rainfall events over different regions of Tanzania. This is achieved by comparing the spatial patterns of and rainfall. Starting with the first case, Figure 3 shows the distribution of rainfall estimates derived from satellite, where subplot 3(a) and 3(b) respectively represent accumulated rainfall estimates on 20th and 21st of December, 2011. While subplot 3(c) is 48-hour accumulated rainfall estimates from 20th to 21st December, 2011. It is seen that the distribution of rainfall almost covered the whole domain of Tanzania. However, there are differences in rainfall intensity across the regions. The coastal regions and southwestern highlands received larger amount of rainfall compared to Northern regions and northeastern highlands. The maximum 24-hours accumulated rainfall is seen over Dar es Salaam (Figure 3(a) and Figure 3(b)). These rainfall estimates from satellite collocate with observed rainfall in Figure 4, where maximum rainfall of about 156 mm was measured over Dar es Salaam weather station on 20th December, 2011. Further, Figure 4 also shows that most regions received rainfall above 22 mm on 20th December, 2011.
The distribution of the share similar patterns to that of rainfall areas,
Figure 3. Rainfall in mm/day (a) 20/12/2011; (b) 21/12/2011; (c) 20-21/12/2011.
Figure 4. Distribution of observed rainfall over different regions in Tanzania on 20-12-2011.
Figure 5 shows the distribution of the first component of, where subplot 5(a) and 5(b) presents the first component of at 700 hPa on 20th and 21st December 2011 respectively, while subplot 5(c) is the 48-hour average of at 700 hPa. The reason why we calculate at 700 hPa level and compare with the rainfall areas is due to the fact that, this level is commonly used for diagnosis of rainfall triggering systems over the tropics, particularly in Tanzania. Thus the behaviour of computed at this level is expected to better represent the pattern of rainfall over different regions compared to at other levels. The positive values of are seen distributed almost over the whole domain of Tanzania (Figure 5(a)). This collocates with the distribution of rainfall in Figure 3(a). Positive values of extended from Congo to western part of Tanzania, where they extended again to northern and southern regions. This distribution may have contributed to the observed high amount of rainfall over the southwestern highlands, where high ground topography may be enhancing convection processes and over southern part of Lake Victoria, where lake-land contrast may be triggering convection activity. However, a negative band of is seen along the coast (Figure 5(a)). This band is surrounded by positive values of, where maximum values are seen over the Indian Ocean. Figure 5(b) shows the band of negative extends from the coast to inland, this indicates convective processes over the coast regions reduced on 21 of December 2011.
Figure 6 shows the distribution of the second component of, where subplot 6(a) and 6(b) present the patterns of on 20th and 21st, December 2011 respec-
Figure 5. Spatial distribution of the first component of the moist potential vorticity vector (in PV-units) on (a) 20/12/2011; (b) 21/12/2011; (c) 20-21/12/2011.
Figure 6. Similar to Figure 5 but for the second component of the moist potential vorticity vector (in PV-units) on (a) 20/12/2011; (b) 21/12/2011; (c) 20-21/12/2011.
tively. Subplot 6(c) is the average of the second component of from 20th to 21st December 2011. The distributions of the second component of collocate with the distribution of rainfall locations. Positive values of are distributed almost over the whole of Tanzania. Similar to what was observed with the first component, the coastal region is dominated with negative values of which are surrounded with positive values of. From Figure 6(a), maximum positive values of are seen over southwestern highland and southern part of the lake Victoria region and parts of central regions this correspond to high rainfall over those areas.
The distribution of the third component of is shown in Figure 7. From this figure it can be seen that subplot 7(a) and subplot 7(b) which represents the patterns of the third component of on 20th and 21st December 2011 respectively are dominated with positive values. The 48-hour average of the third component of shows that the country was dominated with small positive values (Figure 7(c)).
The magnitude of the shown in Figure 8 indicates that the band of maximum values of is located over the coastal regions. This band bear different shapes, on 20th December, it is oriented in south-eastern direction touching the coast, with maximum values of located close to Dar es Salaam. This collocated with the observed heavy rainfall over that region that triggered catastrophic flooding event. On 21st December 2011, the band of maximum magnitude of is oriented
Figure 7. Similar like Figure 6 but for third component of the moist potential vorticity vector (in PV-units) on (a) 20/12/2011; (b) 21/12/2011; (c) 20-21/12/2011.
Figure 8. The magnitude of the moist potential vorticity vector (in PV-units) on (a) 20/12/2011; (b) 21/12/2011; (c) 20-21/12/2011.
parallel to the coast but also touches most regions along the coast. This collocated with the observed rainfall over the coastal regions. Generally it is seen that the area close to maximum values of the magnitude of experienced heavy rainfall (Figure 3(a), Figure 4 and Figure 8(a)). To compare the strength of convective systems along the coast regions between day 20th and 21st it can be seen that convection processes was stronger along the coastal regions on 20th than on 21st of December, 2011.
The second case analysed in this study is that occurred on 5th to 7th May, 2015. Figure 9 present distribution of rainfall estimates from the satellite, where subplot 9(a), 9(b) and 9(c) presents rainfall distribution on 5th, 6th and 7th May 2015 respectively. It can be seen that most areas received rainfall on the 5th where the coastal region received more rainfall compared to the other regions (Figure 9(a)). On the 6th the amount of rainfall has increased over coastal regions, where the city of Dar es Salaam received heavy rainfall which triggered catastrophic flooding. On the 7th rainfall has increased over the entire coastal regions.
Figures 10-12 represents the distribution of observed rainfall over different regions. These figures present similar patterns as the satellite rainfall estimates. However, the amount of rainfall estimated by satellite is lower than that from the weather stations. The coastal region received more rainfall than other regions and the city of Dar es Salaam received about 111mm of rain on 6th May 2015.
The patterns of follow similar patterns to that of the rainfall, Figure 13
Figure 9. Rainfall in mm/day (a) 5/5/2015; (b) 6/5/2015; (c) 7/5/2015; (d) 5-7/5/2015.
Figure 10. Distribution of observed rainfall over different regions in Tanzania on 5-5-2015.
Figure 11. Distribution of observed rainfall over different regions in Tanzania on 6-5-2015.
Figure 12. Distribution of observed rainfall over different regions in Tanzania on 7-5-2011.
shows the distribution of the first component of; where on 5th positive values of dominate almost over the entire domain of Tanzania with maximum values located over the northern coast. However, on 6th the positive values of squeezed to form a kidney shape, extending from northern to parts of southern coastal regions. This band of positive values of contributed to heavy rainfall over Dar es Salaam. On the 7th the band of positive values of MPVV over the coast relaxed and moved to parts of southwestern high grounds. This might contributed to heavy rainfall events observed over Mahenge (Figure 13).
The distribution of second component of is shown on Figure 14. Similar pattern to that of the first component is seen. However, the second component shows that on 5th the maximum positive values of were located close to northern Tanzania Kenya boarder. It moved over northern coast on the 6th and influenced heavy rainfall over those areas. On the 7th the entire coast was dominated with positive values of the second component of. The third component of shows little contribution to observed rainfall as it indicates the country was dominated with small positive values of on the 6th without showing the regions with maximum values especially to the area where received maximum rainfall (Figure 15).
The magnitude of shows that maximum values of were located on the northern Tanzania to Kenya boarder on the 5th, it then moved to touch the entire northern coast of Tanzania on the 6th (Figure 16). In general the areas with maximum values of the magnitude of collocate with areas received heavy rainfall.
Figure 13. The first component of the moist potential vorticity vector (in PV-units) (a) at 2015:5:5:00-2015:5:5:22; (b) 2015:5:5:22-2015:5:6:21; (c) 2015:5:6:21-2015:5:7:20; (d) 2015:5:5:00- 2015:5:7:20.
Figure 14. The second component of the moist potential vorticity vector (in PV-units) (a) at 2015:5:5:00-2015:5:5:22; (b) 2015:5:5:22-2015:5:6:21; (c) 2015:5:6:21-2015:5:7:20; (d) 2015:5:5:00- 2015:5:7:20.
Figure 15. The third component of the moist potential vorticity vector (in PV-units) (a) at 2015:5:5:00-2015:5:5:22; (b) 2015:5:5:22-2015:5:6:21; (c) 2015:5:6:21-2015:5:7:20; (d) 2015:5:5:00- 2015:5:7:20.
Figure 16. The magnitude of moist potential vorticity vector (in PV-units) (a) at 2015:5:5:00- 2015:5:5:22; (b) 2015:5:5:22-2015:5:6:21; (c) 2015:5:6:21-2015:5:7:20; (d) 2015:5:5:00-2015:5:7:20.
This gives more confidence that can be used to diagnose the occurrences of heavy rainfall.
4. Conclusion and Recommendations
In this study, we presented a first step of new paradigm, to use the Moist Potential Vorticity Vector as diagnostic variable of heavy rainfall events in Tanzania. The main purpose of the article was to compute and compare the patterns of derived from the gradient of the moist-air entropic potential temperature and heavy rainfall events that occurred over different regions in Tanzania on 20th to 22nd December, 2011 and on 5th to 8th May, 2015. Moreover, the article aimed at assessing the relative contributions of horizontal and vertical components of MPVV detecting on the location of the observed heavy rainfall events. Dynamic and Thermodynamic variables: wind speed, temperature, Atmospheric pressure and relative humidity at 600 hPa and 800 hPa levels generated by the Weather Research and Forecasting (WRF) Model running at Tanzania Meteorological Agency (TMA) were used to compute MPVV at 700 hPa level. It is found that can serve as a diagnostic variable of heavy rainfall events over different regions in Tanzania. The first and the second components of contribute almost equally to locate the areas which received heavy rainfall. Maximum values of are found to coincide with the areas received heavy rainfall. Although we got good relation between with heavy rainfall but only five days data were used, therefore we recommended more studies to be done to validate the usefulness of. Furthermore, it might be interesting if can be tested to be used as a predictor in climate change studies to help in downscaling climate change projections.
Authors are grateful to the Tanzania Meteorological Agency for provision of observed meteorological data, and the output from WRF model which have been used in this study. Special thanks to Pascal Marquet from the Météo-France, CNRM/GMAP/PROC for the useful discussion on computation of his new novelty moist air entropic potential temperature.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
 Hally, A., Richard, E., Fresnay, S. and Lambert, D. (2014) Ensemble Simulations with Perturbed Physical Parametrizations: Pre-HyMeX Case Studies. Quarterly Journal of the Royal Meteorological Society, 140, 1900-1917.
 Global Weather Hazards Summary (GWHS) (2016) Flooding Risks Elevated as Heavy Rains Forecast to Continue in Southern Tanzania.
 CNRS (Délégation Paris Michel-Ange) (2013) Elucidating Heavy Precipitation Events. Science Daily.
 Mahmud, M. and Ijaya Kumar, T.S.V.V. (2008) Forecasting Severe Rainfall in the Equatorial Southeast Asia. Geofizika, 25, 109-128.
 Kijazi, A.L. and Reason, C.J.C. (2009) Analysis of the 2006 Floods over Northern Tanzania. International Journal of Climatology, 29, 955-970.
 Kapala, A., Born, K. and Flohn, H. (1994) Monsoon Anomaly or an El Niño Event in the Equatorial Indian Ocean? Catastrophic Rains 1961/62 in East Africa and Their Teleconections. WMO/TD, 619, 119-126.
 Anyah, R.O. and Semazzi, F.H. (2006) Climate Variability over the Greater Horn of Africa Based on NCAR AGCM Ensemble. Theoretical and Applied Climatology, 86, 39-62.
 Liang, Z., Lu, C. and Tollerud, E.I. (2010) Diagnostic Study of Generalized Moist Potential Vorticity in a Non-Uniformly Saturated Atmosphere with Heavy Precipitation. Quarterly Journal of the Royal Meteorological Society, 136, 1275-1288.
 Zhang, D.-L. and Cho, H.-R. (1992) The Development of Negative Moist Potential Vorticity in the Stratiform Region of a Simulated Squall Line. Monthly Weather Review, 120, 1322-1341.
 Gao, S.T., Wang, X.R. and Zhou, Y.S. (2004) Generation of Generalised Moist Potential Vorticity in a Frictionless and Moist Adiabatic Flow. Geophysical Research Letters, 31, L12113.
 Schubert, H.W., Hausman, S.A., Garcia, M., et al. (2001) Potential Vorticity in a Moist Atmosphere. Journal of the Atmospheric Sciences, 58, 3148-3157.
 Lo, J.C., Yang, Z.L. and Pielke, R.A. (2008) Assessment of Three Dynamical Climate Down-scaling Methods Using the Weather Research and Forecasting (WRF) Model. Journal of Geophysical Research, 113.
 Ratana, S.B., Ratnam, J.V., Behera, A.K., Rautenbach, C.J., Ndarana, T., Takahashi, K. and Yamagata, T. (2013) Performance Assessment of Three Convective Parameterization Schemes in WRF for Downscaling Summer Rainfall over South Africa. Climate Dynamics, 42, 2931-2953.
 Marquet, P. (2011). Definition of a Moist Entropic Potential Temperature. Application to FIRE-I Data Fights. Quarterly Journal of the Royal Meteorological Society, 137, 768-791.
 Marquet, P. (2014) On the Definition of a Moist-Air Potential Vorticity. Quarterly Journal of the Royal Meteorological Society, 140, 917-929.
 Marquet, P. (2015). An Improved Approximation for the Moist-Air Entropy Potential Temperature θ_s.