Front. Earth Sci. Frontiers in Earth Science Front. Earth Sci. 2296-6463 Frontiers Media S.A. 10.3389/feart.2017.00007 Earth Science Original Research Evaluation of the Use of Moist Potential Vorticity and Moist Potential Vorticity Vector in Describing Annual Cycles of Rainfall over Different Regions in Tanzania Luhunga Philbert M. 1 2 * Djolov George 1 1Department of Geography, Geo-Informatics and Meteorology, University of Pretoria Pretoria, South Africa 2Tanzania Meteorological Agency, Research Section Dar es Salaam, Tanzania

Edited by: Tercio Ambrizzi, University of São Paulo, Brazil

Reviewed by: Eduardo Zorita, Helmholtz-Zentrum Geesthacht Centre for Materials and Coastal Research (HZ), Germany; Meiry Sayuri Sakamoto, Fundacao Cearense de Meteorologia e Recursos Hidricos, Brazil

*Correspondence: Philbert M. Luhunga philuhunga@yahoo.com

This article was submitted to Atmospheric Science, a section of the journal Frontiers in Earth Science

 03 02 2017 2017 5 7 10 06 2016 19 01 2017 Copyright © 2017 Luhunga and Djolov. 2017 Luhunga and Djolov

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

The economy of Tanzania heavily depends on agriculture sector which is primarily rain-fed. In this paper, we compute the moist potential vorticity (MPV) and evaluate its usefulness to describe annual cycles of rainfall. We also modify the convective vorticity vector (CVV→) which was defined as the cross product of absolute vorticity and the gradient of equivalent potential temperature to moist potential vorticity vector (MPVV⃗). This vector is calculated as a cross product of absolute vorticity and the gradient of moist air entropic potential temperature. The performance of MPVV⃗ to describe the annual cycles of rainfall over different regions in Tanzania is analyzed. Twenty six years (1976–2001) daily data of air temperature, specific humidity, zonal and meridional components of the wind at 850 and 600 hPa from numerical output generated by the Rossby Center regional climate model version four (RCA4) are used for computation of MPV and MPVV⃗ at 700 hPa. The statistical relationship between MPVV⃗ and MPV against observed rainfall data from 22 synoptic meteorological stations using Pearson correlation coefficient indicates that MPVV⃗ bears a stronger and more statistically significant correlation coefficient to rainfall than MPV suggesting its potential use as predictor of annual cycles of rainfall over different regions in Tanzania.

 moist potential vorticity moist potential vorticity vector moist air entropic potential temperature Introduction

The general circulation models (GCMs) represent the most satisfactory approach for predicting climate change (IPCC, 2013). They describe the relevant physical processes in the atmosphere, hydrosphere, and cryosphere that make up the climate system. However, GCMs have coarse spatial resolutions and cannot resolve small scale features such as orography, and land use land change that characterize the climate of many regions in the world. This makes their climate simulations of limited use in impact studies of climate change on for instance biodiversity, ecosystem services, agricultural systems, species distributions, conservation planning, and other landscape related matters (Villegas and Jarvis, 2010; Daniels et al., 2012; Tumbo et al., 2012; Xiaoduo et al., 2012; Hassan et al., 2013; Vigaud et al., 2013). These types of impact studies require climate information with much finer spatial resolution.

Downscaling of GCMs outputs is a widely applicable technique for obtaining high resolution climate information that takes into account regional patterns and valuable local knowledge. It is defined as a process of making a link between the state of some atmospheric variable representing a large space (henceforth referred to as the “large scale”) and the state of some atmospheric variable representing a smaller space (henceforth referred to as the “small scale”) (Benestad et al., 2007).

There are two broad categories of downscaling techniques (Hewitson and Crane, 1996). The first category is the dynamical downscaling. This is based on nesting a high resolution regional climate model (RCM) within GCM and drives it using boundary condition from GCM (Danis et al., 2002). The second category is statistical downscaling. This is based on establishing statistical links between large scale atmospheric variables (predictors) and local scale atmospheric variables (predictands).

Dynamical downscaling technique has been extensively used to provide high resolution climate simulation over different regions (Danis et al., 2002; Jones et al., 2004; Roux, 2009; Wilby and Fowler, 2011; Xiaoduo et al., 2012). However, this technique suffers to reproduce the spatial and temporal distributions of climate variables with strong spatial and temporal variability such as rainfall. This is due to fact that the contemporary numerical grids of the RCMs are still too coarse to represent all drivers of rainfall at local scales, circulation patterns, small scale topography, thunderstorms and cloud micro-physics processes (Goosse et al., 2010; Wilby and Fowler, 2011).

Statistical downscaling can be employed to better or adjust the output from RCMs (Benestad et al., 2007). The statistical downscaling techniques require the selection of realistic predictors that are relevant to the predictands. Several researchers (Zorita and Von storch, 1999; Chen et al., 2010; Villegas and Jarvis, 2010; Hassan et al., 2013; Vigaud et al., 2013; Muchuru et al., 2014) have used predictors such as, sea level pressure, Sea Surface Temperature (SST), geo-potential height, wind fields, relative humidity, or temperature variables in statistical downscaling to develop predictors-predictands transfer functions.

However, some of the predictors that are used for statistical downscaling are debated in literatures. For instance, the sea surface temperatures (SSTs), which partly depends on the ocean dynamics, is not represented well in ocean models (Benestad et al., 2007). The spatial resolution of these models tends to be too coarse to describe the ocean currents which are important influences on the SSTs. Fung et al. (2011) argued that circulation predictors alone are unlikely to capture precipitation mechanisms linked to thermodynamics and vapor contents. Wilby and Wigley (2000) suggest that atmospheric moisture must be considered as predictors as well as atmospheric circulations. Charles et al. (1999) suggest that inclusion of moisture variables as predictors can lead to convergence between statistical and dynamical downscaling approaches.

Therefore there is no consensus in literature about the most appropriate predictor to be used in statistical downscaling (Fung et al., 2011). In this study we explore the use of MPV to describe annual cycles of rainfall over different regions of Tanzania to see if it can be used as a predictor for downscaling climate change projections. The same task is performed using the moist potential vorticity vector (MPVV⃗) which is a modified version of convective vorticity vector (CVV⃗) defined by Gao et al. (2004b).

Background information on potential vorticity

The concept of Potential Vorticity (PV) has long history in the study of fluid dynamics. It has been used in meteorology and oceanography for many years back (see Bjerknes, 1898a; Rossby, 1939; Ertel, 1942). This concept has many applications in meteorology, oceanography and aerodynamics (Hoskins et al., 1985). It is mentioned in Hoskins et al. (1985) that “PV can be used to understand the dynamics and thermal conditions of atmospheric flow to the lower limit to the fineness of the structures that may occur, all the way down to the length scales at which molecular diffusion acts.” Many important synoptic scale processes can be understood within the framework of PV. Recently McIntyre (2015) indicated that PV can explain the balanced flows and basic dynamical processes of large scale features such as breaking and propagation of Rossby-wave and its many consequences in the Earth's atmosphere. PV can describe global-scale teleconnections, anti-frictional phenomena such as jet stream self-sharpening, and the genesis of cyclones, anticyclones and storm tracks. PV is conservative and is subject to invertibility principle (Hoskins et al., 1985). The conservative and invertibility properties of PV form the bases of understanding many important atmospheric flow processes.

Ertel (1942) definition of PV see Schubert et al. (2004) is

PV(θ) = ρ−1 ζa · ∇(θ),

where ρ is density of air, ζa is absolute vorticity and ∇(θ) is three dimensional gradient of the potential temperature. This definition emerged from fundamental concepts on circulation and vorticity that have been laid in the works of Bjerknes (1898a,b) and Rossby (1936, 1938, 1940).

Equation (1) is often used to study the thermodynamic properties of the atmosphere (Hoskins and Sardeshmukh, 1987; Stoelinga, 1996; Hoskins, 1997) while it is based only on dry-air potential temperature θ. In moist atmosphere, Equation (1) is not conserved when latent heat release is taken into account (Cao and Cho, 1995; Mofor and Lu, 2008). To avoid this drawback, Bennetts and Hoskins (1979) defined a moist potential vorticity (MPV) by replacing dry-air potential temperature θ with wet bulb potential temperature θw as;

MPV(θw) = ρ-1ζa·∇(θw).

However, θw in Equation (2) cannot fulfill the demand to verify at the same time the conservative property of the moist air and invertibility principle (Marquet, 2014).

Therefore neither Equation (1) nor Equation (2) can be used to study non-uniform saturated atmospheric flow and fulfill the demand to verify, at the same time, a moist and dry air conservative property and an invertibility principle. To overcome this drawback Gao et al. (2004a) defined a new Generalized Moist Potential Vorticity (GMPV) by replacing θ with a Generalized Potential Temperature GMPV(θ*) (GPT) as;

GMPV(θ*) = ρ-1ζa·∇(θ*),

where θ is defined as;

θ*(T,p,q) = θexp(LqsCpT(qqs)k),

where q and qs are specific humidity and saturated specific humidity respectively, (qqs)k is a condensation probability function. In case of absolutely dry atmosphere where q = 0, Equation (4) reduces to dry potential temperature θ*(T, p, q) = θ while in completely saturated atmosphere where q = qs it reduces to equivalent potential temperature θ*(T,p,q)=θe=θexp(LqsCpT). In realistic atmosphere which is non-uniformly saturated, the introduction of condensation probability function fixes the discontinuity of latent heat term due to the impact of water phase changes in the thermodynamic equation. Therefore a smooth transition from completely dry atmosphere and saturated atmosphere is achieved through the change of specific humidity from q to qs. Equation (4) has been used in computation of MPV (Gao et al., 2004a; Mofor and Lu, 2008; Liang et al., 2010; Yang et al., 2014), and found that the solenoidal term does not cancel out in the MPV tendency equation in moist and dry atmosphere. However, Gao et al. (2004a) noted some limitations of applicability of the condensation density function in regions of no condensation or lower relative humidity conditions.

Recently Marquet (2014) defined a new MPV using the specific entropy formulation expressed in terms of moist-air entropy potential temperature (θs) as;

MPV(θs) = 1ρζa·∇(θs),

where θs is defined in Marquet (2011) as;

θs ≡(θs)1(TTr)⋋qt(ppr)-kδqt(rrrv)γqt(1 + ηrv)k(1 + δqt)(1 + ηrr)kδqt,  (θs)1 = θexp (Λrqt) exp(-Lvql + LsqicpdT),

where Λr=(sv0 - sd0)cpd≈5.87 is a key quantity. It is mentioned in Marquet (2011) that Λr depends on the standard entropies of water vapor and dry air (sv0 and sd0). It is also mentioned that (θs)1 is a good approximation of θs. For detailed derivation of θs the reader may consult Marquet (2011, 2014). The main advantage of θs, is that it represent exactly the moist air entropy, it is valid for a general mixing of dry air, water vapor and all possible condensed water species. It is mentioned in Marquet (2014) that “θs verifies the same conservative properties as the moist entropy, even for varying dry air or total water content”. The moist formulation for θs is valid for a general mixing of dry air, water vapor, and all possible condensed water species (Marquet, 2011). In this paper, we compute Marquet (2014)'s MPV formulation and explore its usefulness as a predictor of the annual cycles of rainfall over different regions of Tanzania.

We also suggest a modification of the convective vorticity vector (CVV⃗) proposed by Gao et al. (2004b) by replacing the equivalent potential temperature θe with conservative moist-air entropy potential temperature (θs) to form moist potential vorticity vector (MPVV⃗). This vector is also evaluated as a predictor of annual cycles of rainfall over different regions in Tanzania. We argue that the results from the scalar product of absolute vorticity and gradient of temperature may not explain all atmospheric dynamics over the tropics that contribute to formation and distribution of rainfall events. This is due to (1) the Coriolis parameter that contributes for many dynamical processes in mid and extra tropical regions is very small over the tropics and is zero over the equator. Therefore equatorial flows, especially two dimensional equatorial flows, may not be explained by any version of MPV as Gao et al. (2004b) argued. (2) The vertical gradient of temperature over the tropics is small due to strong convective mixing processes. The scalar MPV derived from dot product of absolute vorticity and gradient of temperature may not be accurate enough to explain convective processes associated with rainfall events over tropics.

 Data and analysis Study area

Tanzania (Figure 1) is located in East Africa between longitudes 29° to 41°E and latitudes 1° and 12°S. The country has an area of 945,000 km2 of which 884,000 km2 is land mass and 61,000 km2 is lakes, rivers and seashore. Tanzania has complex topography that is very heterogeneous. The height of the topography ranges from sea level in the East to 1600 m in the West. In the northeastern highlands is the highest mountain in Africa: Mt Kilimanjaro with an altitude of 5895 m, while in the north is the largest lake in Africa: Lake Victoria. In the south is Lake Nyasa and the Ruvuma river, and in the west is the deepest lake in Africa: Lake Tanganyika. Much of the country lies above 1000 m altitude with many areas over central and northern regions above 1500 m.

Map of Tanzania showing the location of meteorological stations in Unimodal and Bimodal regions separated by the red line.

The climate over Tanzania is mainly controlled by the movement of the Inter-Tropical-Convergence-Zone (ITCZ). However, seasonal interactions within the ITCZ, perturbations in global climate circulation and changes in local circulation systems which are influenced by complex topographical features all contribute to high local climate variability. The seasonal rainfall is modulated by changes in the global sea surface temperatures (SSTs) especially over the equatorial Pacific and Indian Oceans (Black et al., 2003; Black, 2005; Anyah and Semazzi, 2007).

Tanzania is characterized by two rainfall seasons, namely March-April-May (MAM) and October-November-December (OND). These seasons are mainly driven by the migration of the ITCZ, which lags behind the overhead sun by 3–4 weeks over the region (Luhunga et al., 2016). The ITCZ migrates toward southern regions of Tanzania in October-December, reaching southern parts of the country in January-February and reverses northwards in March, April and May (Timiza, 2011). Due to this movement, some areas experience single and double passages of the ITCZ (Luhunga et al., 2016). Regions with a single passage are known as unimodal areas (see Figure 1) and include the southern, southwestern, central and western parts of the country which receive rainfall from October through to April or May (Timiza, 2011). Areas that experience a double passage are known as bimodal areas, and include north, northern coast, northeastern highlands, the Lake Victoria basin and the islands of Zanzibar (Unguja and Pemba). These regions receive two distinct rainfall seasons; the long rain season (known as Masika in Swahili) which starts in March and continues through May (MAM) and the short rainfall season (Vuli in Swahili) which starts in October and continues through December (OND) (Agrawala et al., 2003).

The amount of seasonal rainfall varies significantly in space and time, with higher variation observed during the Vuli season than in Masika. The rainfall falling in these seasons usually ranges from 50 to 200 mm per month but varies greatly between regions and can be as much as 300 mm per month in wettest regions and seasons (McSweeney et al., 2010). Higher amounts of seasonal rainfall are recorded over the southwestern and northeastern highlands, while central Tanzania is semi-arid, receiving seasonal rainfall of less than 50 mm per month. Annual average rainfall over Tanzania ranges from 534 to 1837 mm.

 Model data

Rossby center has produced and made available a very large number of climate simulation that are included in the Coordinated Regional Downscaling Experiment (CORDEX). Rossby Centre Atmospheric (RCA) model is driven by several general circulation models (GCMs) and ERA-Interim global reanalysis datasets (Strandberg et al., 2014). In this study daily data of air temperature, specific humidity, zonal and meridional components of the wind at 850 and 600 hPa from numerical output generated by the Rossby Center regional climate model version four (RCA4) forced by CCCma-CanESM2 (GCM) are used for computation of MPV and MPVV⃗ at 700 hPa. RCA4 has high space resolution of 0.40 by 0.40 corresponding to 50 km by 50 km and has 10 vertical levels (standard pressure levels).

The land surface characteristics used to initialize Soil-Vegetation-Atmosphere Transfer schemes (SVTs) in RCA4 come from a new global physiography data bases ECOCLIMAP for vegetation, lake depth and soil carbon density, and Gtopo30 for orography. The near surface diagnostic quantities: temperature, specific humidity and wind speed are solved using new version of Turbulence Kinetic Energy (TKE) scheme (Strandberg et al., 2014). This scheme combines TKE based on local stability measurement (using Richardson number) and the non-local parcel method. Convection scheme in RCA4 is based on Bechtold-KF (Strandberg et al., 2014). In this scheme the triggering function forcing for convection is the large scale vertical velocity, closure assumption and cloud top is based on CAPE closure (Bechtold et al., 2001).

 Observation data

Monthly observation rainfall data from 22 synoptic meteorological stations over the period of 1976–2001 are obtained from the Tanzania Meteorological Agency (TMA). These data are quality controlled to remove the inhomogeneity and gaps by using the HOMER software package. For detailed descriptions of methodology used for homogeneity tests the reader may consult Luhunga et al. (2014).

 Analysis Marquet (<xref ref-type="bibr" rid="B32">2014</xref>)'s MPV

For simplification of the analysis Equation (5) can be re-written into three terms under hydrostatic balance as;

MPV(θs) = -g(-∂v∂p∂θs∂x + ∂u∂p∂θs∂y+(ζ + f)∂θs∂p), first term = MPVz=-g(ζ + f)∂θs∂p, second term = MPVx=g∂v∂p∂θs∂x, third term = MPVy=-g∂u∂p∂θs∂y. The moist potential vorticity vector (<inline-formula><mml:math id="m28"><mml:mover class="overrightarrow"><mml:mrow><mml:mi>M</mml:mi><mml:mi>P</mml:mi><mml:mi>V</mml:mi><mml:mi>V</mml:mi></mml:mrow><mml:mo>⃗</mml:mo></mml:mover></mml:math></inline-formula>)

Gao et al. (2004b) defined convective vorticity vector as;

CVV⃗=ζa × ∇θeρ.

We modify Equation (12) by replacing θe with θs to form MPVV⃗ as;

MPVV⃗(θs) = ζa × ∇θsρ,

Equation (13) can be re-written into component form as;

MPVV→(θs)=1ρ(−∂v∂z∂u∂z(ζ + f))  ×  (∂θs∂x∂θs∂y∂θs∂z)                            =1ρ(∂u∂z∂θs∂z − (ζ + f)∂θs∂y∂v∂z∂θs∂z + (ζ + f)∂θs∂x−∂v∂z∂θs∂y   −   ∂u∂z∂θs∂x)=(MPVV→xMPVV→yMPVV→z)

where ρ is density which is defined as

ρ=pθv(p,T,qv)Rd(p0p)k,

where p is atmospheric pressure (in Pa) at different level, p0 is atmospheric pressure at reference level, θv(p, T, qv) is virtual potential temperature, Rd is specific gas constant for dry air and k=Rdcp, cp is specific heat capacity at constant pressure. The absolute vorticity ζa is the is defined as;

ζa = -∂v∂zi⃗ + ∂u∂zj⃗ + (ζ + f)k⃗,

where ζ is the relative vorticity defined as ζ=(∂v∂x-∂u∂y)k⃗+uatan(φ), where a is the radius of the earth and φ is the latitude, f is the coriolis parameter defined as f = 2Ωsin(φ).

Considering the hydrostatic equilibrium ∂/∂z = −ρg∂/∂p, Equation (14) can be re-written as;

first component MPVV⃗x=ρg2∂u∂p∂θs∂p-(ζ + f)ρ∂θs∂y, second component MPVV⃗y=ρg2∂v∂p∂θs∂p + (ζ+f)ρ∂θs∂x, third component MPVV⃗z=g(∂v∂p∂θs∂y+∂u∂p∂θs∂x), magnitude of MPVV⃗is written as|MPVV⃗|=(MPVV⃗x)2+(MPVV⃗y)2+(MPVV⃗z)2.

MPVV⃗x, MPVV⃗y and MPVV⃗z are the x, y and z component of MPVV⃗ respectively and |MPVV⃗| is the magnitude of MPVV⃗, u and v are zonal and meridional winds. It is important to note that the vertical component of the wind w is neglected in the computation of MPVV⃗ due to fact that it is smaller than the horizontal components of the winds, therefore ∂w∂x and ∂w∂y≈0.

 MPV and <inline-formula><mml:math id="m51"><mml:mover class="overrightarrow"><mml:mrow><mml:mi>M</mml:mi><mml:mi>P</mml:mi><mml:mi>V</mml:mi><mml:mi>V</mml:mi></mml:mrow><mml:mo>⃗</mml:mo></mml:mover></mml:math></inline-formula> interpolation and statistical analysis

The MPV and MPVV⃗ are calculated at each grid point. In order to be compared with rainfall data at different meteorological stations the grid values of points MPV and MPVV⃗ are interpolated to the location of meteorological stations using arithmetic mean technique. The daily values of MPV and MPVV⃗ interpolated at each station are used to calculate monthly averages of MPV and MPVV⃗. The Pearson correlation coefficient between observed station rainfall and the MPV and between rainfall and MPVV⃗ is computed at each meteorological station using Equation (21) and Equation (22).

r(R,MPV)=∑ i= 1N(Ri−R¯)(MPVi−MPV¯)∑ i= 1N(Ri−R¯)2∑i=0N(MPVi-MPV¯)2, r(R,MPVV)=∑ i= 1N(Ri−R¯)(MPVVi−MPVV¯)∑ i= 1N(Ri−R¯)2∑i=0N(MPVVi-MPVV¯)2,

where R, MPV and MPVV⃗ are the observed rainfall, moist potential vorticity and moist potential vorticity vector respectively, while i refers to the observed rainfall and MPV or MPVV⃗ pairs and N is the total number of such pairs. The value of r ranges between −1 for the perfect negative relationship to +1 for the perfect positive relationship between two variables. The statistical analysis for significant testing of Pearson correlation coefficient used in this study is called a test of the statistical significance of a regressor which is well documented in many statistical text books (e.g., Rangaswamy, 2006).

 Results

First we motivate, why our computation of MPV and MPVV⃗ is done at 700 hPa. The reason is that over the tropics, especially over equatorial regions, moist air thermodynamics become active within the lower troposphere, thus why weather system that determine the day to day forecast are frequently diagnosed at 850 and 700 hPa. The 850 hPa level is the lower level close to the boundary layer and is usable for diagnosis of weather triggering systems along the coastal regions where boundary layer clouds produced by large amount of moisture flux convergence may determine the formation of rainfall. Away from the coastal regions over high grounds, weather systems that are triggering the formation of rainfall are normally diagnosed at 700 hPa. Generally over the tropics rainfall formation is dominated by the low level clouds and shallow convergence of moist air as large amount of moisture is found over low levels. This is different from mid and extra tropics where rainfall is determined by deep convection triggered by the movement of cold fronts and cut off lows.

Statistical analyses

The Pearson correlation coefficient is the measure of relationship between two variables. It is performed to measure the strength of relationship between MPV and rainfall and between MPVV⃗ and rainfall. The Pearson correlation coefficient above 0.4 is considered relatively strong and correlation between 0.2 and 0.4 is considered moderate and those below 0.2 are considered weak (Mayor and Mesquita, 2015). Further statistical test is carried computing the statistical significance level (p) and the coefficient of variation (R2).

Table 1, indicates the Pearson correlation coefficient between annual cycles of MPVz against rainfall computed as monthly average from 1976 to 2001. It is clear that MPVz indicate relatively strong correlation with observed rainfall at 4 meteorological stations, moderate correlated with rainfall at 6 meteorological stations and weakly correlated with rainfall at 12 meteorological stations. The highest correlation coefficient between MPVz and rainfall is observed at Igeri (r = 0.62, p = 0.032). On the other hand, Table 2 indicates the Pearson correlation coefficient between annual cycles MPVV⃗x against rainfall computed as monthly average from 1976 to 2001. MPVV⃗x, has relative strong correlation with rainfall at 12 meteorological stations and moderately correlated with rainfall at 6 meteorological stations and weakly correlated with rainfall at 4 meteorological stations. As shown in Table 2, Kibaha and Morogoro have the highest correlation coefficients of (r = 0.8, p = 0.002) and (r = 0.79, p = 0.002) respectively.

Indicates p-values and Coefficients of determination (R2) for MPV, values in bold are statistically significance at alpha = 0.05.

Station name MPVz MPVx MPVy MPV
p-values r R2 p-values r R2 p-values r R2 p-values r R2
ARUSHA 0.912 0.04 0.001 0.755 −0.10 0.010 0.491 −0.22 0.049 0.921 0.03 0.001
DIA 0.606 −0.17 0.028 0.977 0.01 0.0001 0.828 −0.07 0.005 0.674 −0.14 0.018
BUKOBA 0.557 0.19 0.036 0.160 −0.43 0.187 0.416 0.26 0.067 0.566 0.18 0.034
DODOMA 0.819 −0.07 0.006 0.472 −0.23 0.053 0.006 −0.74 0.542 0.761 −0.10 0.010
IRINGA 0.360 −0.29 0.084 0.138 −0.45 0.207 0.013 0.69 0.478 0.371 −0.28 0.081
MBEYA 0.048 0.58 0.336 0.747 −0.10 0.011 0.338 0.30 0.092 0.073 0.54 0.287
MOROGORO 0.874 0.05 0.003 0.287 −0.33 0.112 0.332 −0.31 0.094 0.936 0.03 0.001
SONGEA 0.768 −0.10 0.009 0.559 −0.19 0.035 0.197 0.40 0.160 0.770 −0.09 0.009
TABORA 0.100 0.50 0.248 0.319 0.31 0.099 0.007 −0.73 0.537 0.101 0.50 0.246
TANGA 0.404 0.27 0.071 0.815 0.08 0.006 0.485 0.22 0.050 0.411 0.26 0.068
IGERI 0.032 0.62 0.382 0.006 0.74 0.550 0.020 0.66 0.433 0.028 0.63 0.398
ILONGA 0.994 0.00 0.0000001 0.253 −0.36 0.128 0.244 −0.36 0.133 0.970 −0.01 0.0001
KIBAHA 0.517 −0.21 0.043 0.824 −0.07 0.005 0.386 −0.28 0.076 0.526 −0.20 0.041
KIGOMA 0.471 0.23 0.053 0.029 −0.63 0.392 0.114 −0.48 0.231 0.640 0.15 0.023
LYAMUNGO 0.629 −0.16 0.024 0.076 0.53 0.281 0.441 −0.25 0.060 0.663 −0.14 0.020
MLINGANO 0.699 −0.12 0.016 0.683 0.13 0.017 0.343 0.30 0.090 0.771 −0.09 0.009
MOSHI 0.267 0.35 0.121 0.202 0.40 0.157 0.834 −0.07 0.005 0.262 0.35 0.124
MTWARA 0.017 −0.67 0.447 0.367 −0.29 0.082 0.768 −0.10 0.009 0.013 −0.69 0.479
MUSOMA 0.677 −0.13 0.018 0.471 −0.23 0.053 0.216 −0.39 0.149 0.625 −0.16 0.025
MWANZA 0.934 −0.03 0.001 0.190 −0.41 0.165 0.044 −0.59 0.348 0.807 −0.08 0.006
SAME 0.353 −0.29 0.087 0.353 −0.29 0.086 0.232 −0.37 0.140 0.350 −0.30 0.088
ZANZIBAR 0.754 −0.10 0.010 0.778 0.09 0.008 0.983 0.01 0.0001 0.776 −0.09 0.008

Indicates p-values and Coefficients of determination (R2) for MPVV⃗ values in bold are statistically significance at alpha = 0.05.

Station name MPVV⃗x MPVV⃗y MPVV⃗z |MPVV⃗|
p-values r R2 p-values r R2 p-values r R2 p-values r R2
ARUSHA 0.228 0.38 0.142 0.414 0.26 0.068 0.104 −0.49 0.243 0.183 0.41 0.170
DIA 0.006 0.740 0.547 0.303 0.32 0.105 0.131 0.46 0.213 0.006 0.74 0.546
BUKOBA 0.406 0.264 0.070 0.504 0.21 0.046 0.727 −0.11 0.013 0.352 0.30 0.087
DODOMA 0.009 0.712 0.507 0.062 0.55 0.306 0.346 0.30 0.089 0.002 0.79 0.627
IRINGA 0.206 0.393 0.154 0.015 0.68 0.465 0.080 0.52 0.275 0.002 0.80 0.637
MBEYA 0.635 −0.153 0.023 0.363 0.29 0.083 0.211 0.39 0.152 0.503 0.21 0.046
MOROGORO 0.002 0.786 0.618 0.120 0.47 0.224 0.967 −0.01 0.0002 0.003 0.78 0.610
SONGEA 0.111 −0.484 0.234 0.876 0.05 0.003 0.245 0.36 0.132 0.308 −0.32 0.103
TABORA 0.127 0.466 0.217 0.056 0.56 0.319 0.944 −0.02 0.001 0.001 0.83 0.682
TANGA 0.269 0.347 0.120 0.264 0.35 0.123 0.165 0.43 0.183 0.057 0.56 0.316
IGERI 0.457 −0.237 0.056 0.684 0.13 0.017 0.041 −0.59 0.354 0.896 −0.04 0.002
ILONGA 0.007 0.728 0.530 0.116 0.48 0.228 0.996 0.001 0.0001 0.010 0.71 0.501
KIBAHA 0.002 0.801 0.642 0.216 0.39 0.149 0.359 0.29 0.085 0.001 0.83 0.682
KIGOMA 0.352 0.295 0.087 0.195 0.40 0.162 0.036 0.61 0.371 0.005 0.75 0.558
LYAMUNGO 0.677 −0.134 0.018 0.729 −0.11 0.012 0.208 −0.39 0.154 0.628 −0.16 0.024
MLINGANO 0.027 0.633 0.401 0.378 0.28 0.079 0.618 0.16 0.026 0.060 0.56 0.311
MOSHI 0.791 0.086 0.007 0.924 0.03 0.001 0.095 −0.50 0.254 0.801 0.08 0.007
MTWARA 0.123 0.470 0.221 0.643 0.15 0.022 0.974 −0.01 0.0001 0.398 0.27 0.072
MUSOMA 0.984 −0.007 0.00004 0.192 0.40 0.163 0.062 −0.55 0.306 0.268 0.35 0.121
MWANZA 0.423 0.256 0.065 0.031 0.62 0.386 0.083 −0.52 0.270 0.013 0.69 0.478
SAME 0.045 0.588 0.345 0.357 0.29 0.085 0.387 −0.28 0.076 0.050 0.58 0.331
ZANZIBAR 0.019 0.663 0.439 0.398 0.27 0.072 0.119 0.47 0.225 0.021 0.65 0.426

The MPVx has shown correlation coefficient of greater than or equal to 0.4 at 6 meteorological stations (Table 1), while MPVV⃗y has high correlation coefficient of greater or equal to 0.4 at 8 stations (Table 2). The MPVy has correlated with rainfall with correlation coefficient of greater or equal to 0.4 at 7 stations (Table 1), while MPVV⃗z has correlation coefficient of greater or equal to 0.4 at 11 meteorological stations (Table 2). The MPV has relatively strong correlated with rainfall at 4 meteorological stations (Table 1). On the other hand |MPVV⃗| has relatively strongly correlated with rainfall at 14 meteorological stations (Table 2). Overall, the findings implied that MPVV⃗ correlate better with rainfall than the MPV. Also MPVV⃗x and |MPVV⃗| correlated better with rainfall than the other components.

The bivariate regression analysis to examine the predictive power between annual cycle rainfall as independent variable and the MPV and its terms are also presented in (Table 1). This table indicates that MPVz explains better the variation of rainfall over Mbeya (R2 = 0.336, p = 0.048), while MPVx explain better the variation of rainfall over Igeri (R2 = 0.550, p = 0.006). The MPVy, explain better the variation of rainfall over Dodoma (R2 = 0.542, p = 0.006), Iringa (R2 = 0.478, p = 0.013), Tabora (R2 = 0.537, p = 0.007), and Mwanza (R2 = 0.348, p = 0.044). The MPV explain better the variation in rainfall only over Mtwara (R2 = 0.479, p = 0.013).

On the other hand, the bivariate regression analysis is conducted to examine the predictive power between annual cycle rainfall as independent variable and the magnitude and the components of MPVV⃗ are presented in Table 2. Results showed that over Dar es Salaam (DIA), MPVV⃗x explains the variation of rainfall with (R2= 0.547, p = 0.006), followed by |MPVV⃗| with (R2 = 0.546, p = 0.006). MPVV⃗x, explain better the variation of rainfall over Morogoro with (R2 = 0.618, p = 0.002), Ilonga with (R2 = 0.530, p = 0.007), and over Mlingano with (R2 = 0.401, p = 0.027).

|MPVV⃗| explain better the variation in rainfall over Dodoma (R2 = 0.627, p = 0.002), Iringa (R2 = 0.637, p = 0.002), Tabora (R2 = 0.682, p = 0.001), Kibaha (R2 = 0.682, p = 0.001), Kigoma (R2 = 0.558, p = 0.005), Mwanza (R2 = 0.478, p = 0.013) and Zanzibar (R2 = 0.426, p = 0.021). Term3 better explain the variation of rainfall at Igeri (R2 = 0.354, p = 0.041).

The variations of annual cycles in rainfall and MPVV⃗x over unimodal and bimodal regions represented by Tabora and Kibaha stations respectively are presented in Figures 2, 3. It is clear that MPVV⃗x catches the annual cycles of rainfall.

Annual cycle of rainfall and MPVVx⃗ (102*PV-units) at 700 hPa calculated from 1976-2001 at Tabora meteorological station.

Annual cycle of rainfall and MPVVx⃗ (102*PV-units) at 700 hPa calculated from 1976-2001 at Kibaha meteorological station.

 Summary and recommendations

The aim of this study was to compute the moist potential vorticity (MPV) and moist potential vorticity vector MPVV⃗ and compare their performance in describing annual cycles of rainfall over different regions of Tanzania. Results indicated that MPVV⃗ perform better in explaining the variation of rainfall than MPV. MPVV⃗x, and |MPVV⃗| provide strong correlation coefficient with rainfall at significant level less than 0.05 at many stations when compared to other components. This suggests that MPVV⃗x, and |MPVV⃗| can be used as predictors of rainfall over different regions where they have shown strong correlation coefficient, and high coefficient of determination at significant level of 0.05. For instance, in Figures 2, 3 MPVV⃗x captured the annual cycle of rainfall over Tabora in unimodal region and Kibaha in bimodal region. One can construct a transfer function between MPVV⃗x and monthly rainfall and use it to improve climate projections for rainfall or for seasonal climate prediction. This is important especially for Tanzania where seasonal climate forecasting is based on analyzing analog years. This is the simple seasonal prediction technique based on analyzing weather maps that resemble other weather maps for different years within the season in the historical record and normally is based on subjective judgments from human eyes. This method produces fairly inaccuracy seasonal climate prediction (Huijun et al., 2015).

In this study we recommend the use of MPVV⃗ as predictor of annual cycles of rainfall over different regions in Tanzania. MPVV⃗, show the observed pattern of rainfall, that in MAM where there is higher rainfall total, MPVV⃗ is also higher and in OND where there is lower rainfall total, MPVV⃗ is lower too. Therefore a transfer function constructed based on MPVV⃗ may accurately predict the seasonal variation of rainfall over different regions in Tanzania.

Furthermore MPVV⃗ need to be explored more on its application to better seasonal prediction in East Africa where seasonal climate prediction depends among other techniques on simulations from the general circulation models (GCMs) forced by sea surface temperatures. However the GCMs have coarse space resolution to reproduce climate details at different regions. The results presented in this study contribute on the existing predictors that are used for development of empirical models for seasonal climate prediction. Most statistical model are constructed using predictors such as relative humidity, geo potential height, upper level wind speed (e.g., at 500 hPa). These predictors suffer to reproduce at the same time the dynamics and thermodynamics of the atmosphere. However MPVV⃗ reproduces at the same time the dynamics of atmospheric flows (through vorticity) and thermodynamics of atmospheric flows (through the gradient of moist entropy potential temperature).

It is important to note that the MPV and MPVV⃗ presented in this study were computed using data from the RCM driven by GCM. Therefore further studies are recommended to explore the performance of MPV and MPVV⃗ in describing rainfall events in Tanzania using data from RCM driven by ERA-Interim data. Moreover, it is recommended that MPVV⃗ should be tested on ability to reproduce interannual variability of rainfall to have more confidence to use it as predictor of rainfall events.

 Author contributions

The scientific contribution of both authors is significant to the manuscript, the computation and data search was done by PL. The validation of the model was done by GD. Both authors participated fully in writing and analyzing the results.

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Authors are grateful to Tanzania Meteorological Agency, Rossby center for regional climate modeling, and the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR), for provision of data 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.

 References Agrawala S. Moehder A. Hemp A. Van Aalst M. Hitz S. Meena H. . (2003). Development and Climate Change in Tanzania: Focus on Mount Kilimanjaro. Paris: OECD. Anyah R. O. Semazzi F. H. M. (2007). Variability of East african rainfall based on multiyearRegCM3 simulations. Int. J. Climatol. 27, 357–371. 10.1002/joc.1401 Bechtold P. Bazile E. Guichard F. Mascart P. Richard E. (2001). A mass-flux convection scheme for regional and global models. Q. R. J. Meteorol. Soc. Vol. 127, 869–886. 10.1002/qj.49712757309 Benestad R. E. Hanssen-Bauer I. Førland E. J. (2007). An evaluation of statistical models for downscaling precipitation and their ability to capture long-term trends. Int. J. Climsyol. 27, 649–665. 10.1002/joc.1421 Bennetts D. A. Hoskins B. J. (1979). Conditional symmetric instability - a possible explanation for frontal rainbands. Q. J. R. Meteorol. Soc. 105, 945–962. 10.1002/qj.49710544615 Bjerknes V. (1898a). Uber die Bildung von Circulationsbewegung und Wirbeln in reibungslosen Fliissigkeiten. Videnskabsselskapets Skrifter. I Math. Naturu. Klasse. 29. Bjerknes V. (1898b). Uber einen hydrodynamischen Fundamentalsatz und seine Anwendung besonders auf die Mechanik der Atmosphare und des Weltmeeres. Kgl. Suenska Vetenskapsakad. Handl. 31, 35. Black E. (2005). The relationship between indian ocean sea-surface temperature and east african rainfall. Philos. Trans. R. Soc. 1826, 43–47. 10.1098/rsta.2004.147415598619 Black E. Slingo J. Sperber K. R. (2003). An observational study of the relationship between excessively strong short rains in coastal east Africa and Indian ocean SST. Mon. Weather Rev. 131, 74–94. 10.1175/1520-0493(2003)131<0074:AOSOTR>2.0.CO;2 Cao Z. Cho H.-R. (1995). Generation of moist potential vorticity in extratropical cyclones. J. Atmos. Sci. 52, 3263–3281. Charles S. P. Bates B. C. Hughes J. P. (1999). A spatio-temporal model for downscaling precipitation occurrence and amounts. J. Geophys. Res. 104, 31657–31669. 10.1029/1999JD900119 Chen S. T. Yu P. S. Tang Y. H. (2010). Statistical downscaling of daily precipitation using support vector machines and multivariate analysis J. Hydrol. 385, 13–22. 10.1016/j.jhydrol.2010.01.021 Daniels A. E. Morrison J. F. Joyce L. A. Crookston N. L. Chen S. C. McNulty S. G. (2012). Climate projections FAQ. General Technical Report RMRS-GTR-277WWW. Fort Collins, CO: U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station, 32. Danis B. Laprise R. Caya D. Cot^e J. (2002). Downscaling ability of one-way nested regional climate models: the Big-Brother Experiment. Clim. Dyn. 18, 627–646 10.1007/s00382-001-0201-0 Ertel H. (1942). Ein Neuer hydrodynamischer Wirbelsatz. Met. Z. 59, 271–281. Fung F. Lopez A. New M. (2011). Modelling the Impact of Climate Change on Water Resources. Oxford: Wiley-Blackwell. Gao S. Ping F. Li X. Tao W.-K. (2004b). A convective vorticity vector associated with tropical convection: a two-dimensional cloud-resolving modeling study. J. Geophys. Res. 109:D14106. 10.1029/2004JD004807 Gao S. Wang X. Zhou Y. (2004a). Generation of generalised moist potential vorticity in a frictionless and moist adiabatic flow. Geophys. Res. Lett. 31:L12113. 10.1029/2003GL019152 Goosse H. Barriat P. Y. Lefebvre W. Loutre M. F. Zunz V. (2010). Introduction to Climate Dynamics and Climate Modeling. Available online at: http://www.climate.be/textbook Hassan Z. Shamsudin S. Harun S. (2013). Application of SDSM and LARS-WG for simulating and downscaling of rainfall and temperature. Theor. Appl. Climatol. 116, 243–257. 10.1007/s00704-013-0951-8 Hewitson B. C. Crane R. G (1996). Climate downscaling: techniques and application. Clim. Res. 7, 85–95. Hoskins B. (1997). A potential vorticity view of synoptic development. Meteorol. Appl. 325–334. Hoskins B. J. McIntyre M. E. Robertson A. W (1985). On the use and significance of isentropic potential vorticity maps. Q. J. R. Meteorol. Soc. 111, 877–946. Hoskins B. J. Sardeshmukh P. D. (1987). A diagnostic study of the dynamics of the northern hemisphere winter of 1985/86. Q. J. R. Meteorol. Soc. 113, 759–778. 10.1002/qj.49711347705 Huijun W. Ke F. Jianqi S. Shuanglin L. Zhaohui L. Guangqing Z. . (2015). A review of seasonal climate prediction research in china. Adv. Atmos. Sci. 32, 149–168. 10.1007/s00376-014-0016-7 IPCC (2013). Summary for Policy, Makersand In: Climate, Change (2013). The Physical Science Basis, in Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, eds Stocker T. F. Qin D. Plattner G. K. Tignor M. Allen S. K. Boschung J. (New York, NY: Cambridge University Press). Jones R. G. Noguer M. Hassell D. C. Hudson D. Wilson S. S. Jenkins G. J. . (2004). Generating High Resolution Climate Change Scenarios Using PRECIS. Exeter: Met Office Hadley Centre. Liang Z. Lu C. Tollerud E. I. (2010). Diagnostic study of generalized moist potential vorticity in a non-uniformly saturated atmosphere with heavy precipitation. Q. J. R. Meteorol. Soc. 136, 1275–1288. 10.1002/qj.636 Luhunga P. Botai J. Kahimba F. (2016). Evaluation of the performance of CORDEX regional climate models in simulating present climate conditions of Tanzania. J. South Hemisphere Earth Syst. Sci. 66, 32–54. 10.22499/3.6601.005 Luhunga P. Mutayoba E. Ng'ongolo H (2014). Homogeneity of monthly mean air temperature of the United Republic of Tanzania with HOMER. Atmospher. Clim. Sci. 4, 70–77. 10.4236/acs.2014.41010 Marquet P. (2011). Definition of a moist entropic potential temperature. Application to FIRE-I data flights. Q. J. R. Meteorol. Soc. 137, 768–791. 10.1002/qj.787 Marquet P. (2014). On the de?nition of a moist-air potential vorticity. Q. J. R. Meteorol. Soc. 140, 917–929. 10.1002/qj.2182 Mayor Y. G. Mesquita M. D. S. (2015). Numerical simulations of the 1 May 2012 deep convection event over cuba: sensitivity to cumulus and microphysical schemes in a high-resolution model. Adv. Meteorol. 2015:973151. 10.1155/2015/973151 McIntyre M. E. (2015). Potential Vorticity, in Encyclopedia of Atmospheric Sciences, Vol 2, 2nd Edn., eds North G. R. editor-in-chief Pyle J. Zhang F. (New York, NY: Elsevier Science; Academic Press), 375–383. McSweeney C. New M. Lizcano G. (2010). UNDP Climate Change Country Profiles: Tanzania. Available online at: http://www.geog.ox.ac.uk/research/climate/projects/undp-cp/UNDP_reports/Tanzania/Tanzania.lowres.report.pdf (Accessed May 10, 2013). Mofor L. A. Lu C. (2008). Generalized moist potential vorticity and its application in the analysis of atmospheric flows. Prog. Nat. Sci. 19, 285–289. 10.1016/j.pnsc.2008.07.009 Muchuru S. Landman W. A. DeWitt D. Lötter D. (2014). Seasonal Rainfall Predictability over the Lake Kariba Catchment Area, Water SA Vol. 40. Available online at: http://www.wrc.org.za ISSN 0378-4738 Rangaswamy R. (2006). Agricultural Statistics. New Delhi: New Age International (P) Ltd Publishers. Rossby C. G. (1936). Dynamics of steady ocean currents in the light of experimental fluid mechanics. Pap. Phys. Oceanogr. 5, 1–43. 10.1575/1912/1088 Rossby C. G. (1938). On the mutual adjustment of pressure and velocity distributions in certain simple current systems, II. J. Marine Res. 2, 239–263. Rossby C. G. (1939). Relation between variations in the intensity of the zonal circulation of the atmosphere and the displacements of the semi-permanent centers of action. J. Marine Res. 2, 38–55. 10.1357/002224039806649023 Rossby C. G. (1940). Planetary flow patterns in the atmosphere. Q. J. R. Meteorol. Soc. 66, 68–87. Roux B. (2009). Ultra High-Resolution Climate Simulations over the Stellenbosch Wine Producing Region Using a Variable-Resolution Mode, MSc Dissertation, University of Pretoria. Available online at: http://upetd.up.ac.za/thesis/available/etd-11302009-185214/unrestricted/dissertation.pdf Schubert W. Ruprecht E. Hertenstein R. Nieto-Ferreira R. Taft R. Rozo? C. . (2004). English translations of twenty-one of Ertel's papers on geophysical ?uid dynamics. Meteorolo. Zeitschrift. 13, 527–576. 10.1127/0941-2948/2004/0013-0527 Strandberg G. Bärring L. Hansson U. Jansson C. Jones C. Kjellström E. . (2014). CORDEX Scenarios for Europe from the Rossby Centre Regional Climate Model RCA4, Reports Meteorology and Climatology, 116SMHI, SE-60176 Norrköping, Sverige (2014) Stoelinga M. (1996). A potential vorticity-based study of the role of diabatic heating and friction in a numerically simulated cyclone. Mon. Wea. Rev. 124:849–874. 10.1175/1520-0493 Timiza W. (2011). Climate Variability and Satellite – Observed Vegetation Responses in Tanzania. Master thesis Physical Geography and Ecosystem Analysis, Lund University, Seminar series 205, 30 ECTS. Tumbo S. D. Mpeta E. Tadross M. Kahimba F. C. Mbillinyi B. P. Mahoo H. F. (2012). Application of self-organizing maps technique in downscaling GCMs climate change projections for Same, Tanzania. J. Phys. Chem. Earth 35, 608–617. 10.1016/j.pce.2010.07.023 Vigaud N. Vrac M. Caballero Y. (2013). Probabilistic downscaling of GCM scenarios over southern India. Int. J. Climatol. 33, 1248–1263 10.1002/joc.3509 Villegas J. R. Jarvis A. (2010). Downscaling Global Circulation Model Outputs: The Delta Method Decision and Policy Analysis Working Paper No. 1, Centro internacional de agricultura Tropical Wilby L. R. Fowler J. H (2011). Regional Climate Downscaling. Available online at: http://www.pages-perso-julie-carreau.univ-montp2.fr/UM2/Packages_and_Tutorial_files/downscaling.pdf Wilby R. L. Wigley T. M.L (2000). Precipitation predictors for downscaling: observed and general circulation model relationships. Int. J. Climatol. 20, 641–661. 10.1002/(SICI)1097-0088(200005)20:6<641::AID-JOC501>3.0.CO;2-1 Xiaoduo P. Xin L.i, Xiaokang, S. H. I. Xujun H. A. N. Lihui L. U. O. Liangxu W. A. N. G. (2012). Dynamic downscaling of near-surface air temperature at the basin scale using WRF–a case study in the Heihe River Basin,China. Front. Earth Sci. 6, 314–323. 10.1007/s11707-012-0306-2 Yang S. Gao S. T. Lu C. G. (2014). A generalized frontogenesis function and its application. Adv. Atmos. Sci. 31, 1065–1078. 10.1007/s00376-014-3228-y Zorita E. Von storch H. (1999). The Analog method as a simple statistical downscaling technique: comparison with more complicated methods. J. Clim. 12, 2474–2489. 10.1175/1520-0442(1999)012<2474:TAMAAS>2.0.CO;2