# Harmonic Analysis of Precipitation Time Series in Lake Tana Basin, Ethiopia

Authors: Zewdu Alamineh Fetene; Tesfay Mekonnen Weldegerima; Tadesse Terefe Zeleke; Melessew Nigussie
Journal: Advances in Meteorology
 (2018)

---

## Abstract

This study presents harmonic analysis of precipitation observations within the Lake Tana Basin for the periods of 1985–2015. The livelihood of several millions of people within the basin and outside the basin is governed by the precipitation conditions within this basin. Large spatial and temporal variabilities of precipitation can increase the incidence of extreme events such as floods and droughts. It is important to identify the characteristics of these variations, and this study aims at investigating the characteristics of the seasonal and annual cycles of precipitation within the Lake Tana Basin using harmonic analysis. Precipitation data of 31 years from four weather stations were used in the analysis. We then applied harmonic analysis to calculate the amplitude, phase shift, and variance of observation. Detailed characteristics of the first five harmonics are presented and discussed. We found the amplitude of the first harmonic to be173.42,177.93,127.77, and 188.78mm for Debre Tabor, Bahir Dar, Gondar, and Dangila, respectively. This shows that Dangila areas got more rainfall during this fundamental period than others increasing from Gondar to Dangila direction. Also, the variance in the first harmonic is smaller than the variances of other harmonics, and this means that the large variations of the precipitation originate from higher harmonics (short time periods). This shows that precipitation variations are governed mainly by monthly, seasonal, and semiannual variations. The analysis has shown that maximum precipitation for all stations occurred in July and August.


---

## Body

## 1. Introduction

The impacts of climate change and climate variability on human life are widely acknowledged. These led to monitor the behavior of weather and climate variables. Rainfall, as one of the most important climate variables, has a direct and indirect impact on the natural environment and human life. Large spatial and temporal variability of it leads to increased incidence of extreme events such as droughts and floods [1]. Drought has three distinguishing features that depend on the characteristics of precipitation; such as intensity, duration and spatial coverage [2]. Floods include river floods, flash floods, urban floods, and sewer floods and can be caused by intense and/or long-lasting precipitation, dam break, and reduced conveyance due to landslides.Many studies have shown the spatial and temporal changes in the amount, intensity, and duration of rainfall worldwide. For example, a study of rainstorms in China [3] has shown an increase in rainfall in the subtropical areas. Rainfall shows a decreasing trend in Inner Mongolia [4] from 1960 to 2012. Similarly, a study on long-term trend in precipitation of different spatial regions of Bangladesh [5] suggested that, for annual precipitation, a significant increasing monotonic trend was found in whole Bangladesh.One of the advanced techniques of studying the rainfall distributions and variabilities is the harmonic analysis [6]. This analysis technique allows studying different parameters, such as the amplitudes of rainfall, frequency, and phase angle (showing time of maximum rainfall). Harmonic analysis technique has been used by many researchers to investigate the spatial and temporal distribution of rainfall. Hsu and Wallace [7] have investigated seasonal precipitation in global scale using harmonic analysis and found much less seasonal variation of precipitation over the oceans than the land. Horn and Bryson [8] and McGee [9] also studied the seasonal distribution and climatology of precipitation using harmonic analysis. Tarawneh [10] presented rainfall climatology of Saudi Arabia. Tarawneh and Kadioglu [11] and Kadioglu et al. [12] used the harmonic technique to analyze precipitation climatology in Jordan and Turkey, respectively; the results of these two articles agree with each other in terms of the percentages of variance of the first harmonic, which represent more than 90% of variation. This is due to the fact that Jordan and Turkey are located in the Mediterranean region and are affected by similar weather systems. The above studies have demonstrated the efficiency of the harmonic analysis technique in analyzing precipitation time series in terms of the amplitudes of many harmonics and periodicity.Variability and distribution of precipitation within the Upper Blue Nile Basin, where the Lake Tana Basin is found (Figure1), has been investigated by many researchers (e.g., [13, 14]). Similarly, Worku [15] stated that flooding is a familiar event in the Upper Blue Nile Basin and did cause a lot of destruction in the past years. But those all researchers did not explain the features of precipitation in terms of percentage of variance, amplitude, and phase shift since the actual data they used cannot show the hidden features of precipitation distribution at different harmonics. Even though harmonic analysis gives us better understanding of the periodicity of annual, seasonal, and interseasonal rainfall, none of these studies have applied this method in the region. Therefore, this study aims to investigate the seasonal and annual cycles of rainfall in the Lake Tana Basin using harmonic analysis. We have utilized 31-year precipitation data gathered at four separate meteorological stations within the Lake Tana Basin.Figure 1
The location of the Lake Tana Basin in the Ethiopian and the Upper Blue Nile Basin system with meteorological stations ([16]).


### 1.1. Synoptic Climatology of Rainfall in Lake Tana Basin

The climate of Lake Tana Basin is dominated by highland tropical monsoon, seasonal migration of intertropical convergence zone (ITCZ), and associated atmospheric circulations, but the complex topography of the region has also an effect on local climate conditions. Even if the basin is found within the tropics, it is not a hot region since the altitudinal climate controlling effect is significantly higher than its latitudinal controlling effect. In other words, Lake Tana Basin, because of its high altitudes, mostly experiences temperate climate conditions. Thus, the seasonal variation in the climate is low compared to other tropical areas of low altitude [17].Two distinct seasons, dry and wet, are documented. The dry season is from November to May, while the wet season covers the remaining parts of the year, when most of the precipitation takes place. Summer season from June to September is the most important and agriculturally active season, most of the rainfall (70–90% total rainfall) occurs, in the basin. The mean annual rainfall of the area is about 1465 mm with significant spatial variation. The mean annual maximum temperature is 25.5°C, and the mean annual minimum temperature is 10.8°C [18].


## 1.1. Synoptic Climatology of Rainfall in Lake Tana Basin

The climate of Lake Tana Basin is dominated by highland tropical monsoon, seasonal migration of intertropical convergence zone (ITCZ), and associated atmospheric circulations, but the complex topography of the region has also an effect on local climate conditions. Even if the basin is found within the tropics, it is not a hot region since the altitudinal climate controlling effect is significantly higher than its latitudinal controlling effect. In other words, Lake Tana Basin, because of its high altitudes, mostly experiences temperate climate conditions. Thus, the seasonal variation in the climate is low compared to other tropical areas of low altitude [17].Two distinct seasons, dry and wet, are documented. The dry season is from November to May, while the wet season covers the remaining parts of the year, when most of the precipitation takes place. Summer season from June to September is the most important and agriculturally active season, most of the rainfall (70–90% total rainfall) occurs, in the basin. The mean annual rainfall of the area is about 1465 mm with significant spatial variation. The mean annual maximum temperature is 25.5°C, and the mean annual minimum temperature is 10.8°C [18].


## 2. Methodology

### 2.1. The Study Area

The Lake Tana (shown in Figure1) is located in Northwestern Ethiopia (latitude 10.95° and 12.78°N and longitude 36.89° and 38.25°E) with a drainage area of about 15,000 km2 [19]. The Lake Tana Basin (Figure 1), the largest lake in Ethiopia and the third largest in the Nile Basin, is found in this region. It is shared by three administrative zones of the Amhara regional state, namely, North Gondor, South Gondor, and West Gojjam directly. The climate of the area is largely controlled by the movement of the intertropical convergence zone and tropical highland monsoon, which results in a single rainy season between June and September. 70–90% of total rainfall occurs in kiremt season between June and September. The mean annual rainfall varies from 1,200 to 1,600 mm [20], and the mean actual evapotranspiration of the catchments is 773 mm.


### 2.2. Data

The main data used in this study are monthly precipitation for each year in the whole statistical period, which are collected from the National Meteorology Agency (NMA) of Ethiopia. We have used four weather stations (Table1) which spread out over the Lake Tana Basin along with a full statistical period of 31 years (1985–2015). Station selection was based on quality, long range data, and representation of various climate zones in the Basin. Some station selection criteria justify that the study of historical climate variability and change should utilize reliable data that are free of artificial trends or changes. Artifacts of measurement caused by changes in observation practice, equipment, site exposure, and location can lead to misleading results when used in trend analyses [21]. Hence, the homogeneity test of the time series of precipitation data, following reliable procedures (e.g., [22]), have been applied in this study to reduce these artifacts. We have used the total monthly average precipitation for the harmonic analysis.Table 1
Stations of the studying area, locations, elevations, and mean annual rainfall.
Station name
Latitude (degree)
Longitude (degree)
Elevation (m)
Mean annual rainfall (mm)
Bahir Dar
11.57
37.36
1800
1403.17
Debre Tabor
11.86
38.02
2706
1459.40
Gondar
12.60
37.45
1967
1069.40
Dangila
11.25
36.84
2127
1633.10


### 2.3. Harmonic Analysis

Harmonic analysis is commonly used to study periodic variations that can represent the fluctuations or variations in a time series as having arisen from adding together a series of sine and cosine functions [6, 23]. It is based on the mathematical principle that a curve, viewed as a function, may be represented by a series of trigonometric function. These trigonometric functions are “harmonic” in the sense that they are chosen to have frequencies exhibiting integer multiples of the fundamental frequency determined by the sample size (i.e., length) of the data series [24].


#### 2.3.1. Representation of Seasonal and Annual Cycles by Transforming a Cosine Wave

Necessarily, harmonic function laterally shifts in order to have it match the peaks and troughs of a data series. The time series data represented byY, t=0, at origin Y, can be written as a combination of different sine and cosine signals as follows:(1)Y=y´+A1cos2πtT−Ф1+A2cos4πtT−Ф2+⋯+Ancos2nπtT−Фn,where y´ is an arithmetic average of the original data series, while the rest of the terms represent the first,second,…,nth harmonics of the period T and Ф is called the phase angle or phase shift which can then be easily interpreted as corresponding to the time of maximum harmonic function.This time shifting is most conveniently accomplished when the cosine function is used because its maximum value is achieved when the angle on which it operates is zero. Shifting the cosine function to the right by the angleФ results in a new function that is maximized at Ф.Equation (1) can be rewritten as a series:(2)Y=y´+Ccos2πtT−Ф.The argument of cosine function is zero when(3)Ф=2πtT.The function is maximized at time(4)ti=ФiT2π,where ti is the ith month and Фi is the phase shift for the corresponding month of the data series.


#### 2.3.2. Estimation of the Harmonic Parameters

In order to calculate amplitude (C) and phase shift (Ф) of a periodicity for which the data comprise a perfect sine curve for a known period T, direct substitution is being applied. For representing the associated precipitation data reasonably well, better choices for C and Ф can be found [6, 23]. The easiest way to do this is by applying trigonometric identity approach.(5)cosα−Ф=cosФcosα+sinαsinФ.Forα=2πt/T and multiplying both sides by amplitude C reveals the following equation:(6)Ccos2πtT−Ф=CcosФcos2πtT+CsinФsin2πtT+Acos2πtT+Bsin2πtT,where A=CcosФ∧B=CsinФ.(7)Y=y´+Acos2πtT+Bsin2πtT,where(8)C=A2+B21/2.Fornth order of harmonics, (6) can be rewritten as follows:(9)Y=y´+∑i=1naicos2πtiT+bisin2πtiT.Similarly,(10)Ai=ai2+bi21/2,(11)Фi=tan−1biai.For the special situation where the data values are equally spaced in time with no missing values, the properties of the sine and cosine functions allow the same least-squares parameter values to be obtained more easily and efficiently using the following equation:(12)ai=2N∑t=1N−1Ycos2πtiT,bi=2N∑t=1N−1Ysin2πtiT,where i=1,2,3,…,N/2, aiandbi are coefficients of the ith harmonic given by [6, 25], and Y is the monthly mean precipitation of the ith month, and N is the number of observations.The variance of harmonic can be also obtained using Wilks [6]:(13)Vi=Ai22.The percentage of the total month-to-month variations explained by a given harmonic can be determined by forming a ratio of the square of the amplitude of that specific harmonic to the sum of all harmonics. The first harmonic can show the observed data of precipitation as a single annual cycle, whereas the second harmonic can represent the tendency towards a semiannual variation of precipitation that cannot appeared from the observed data [10]. With the same concept, the third harmonic describes more detail of the annual variation of precipitation and so on. In brief, the first harmonics show an annual cycle with the greatest amplitude as compared to other harmonics. The second harmonic indicates a semiannual change, while the third one illustrates four-month seasonal change in details [12]. The phase angle shift indicates the displacements of the maximum along the time axis, which can be represented by the x-axis. The phase angle shift can express the precipitation regime and the boundaries [26].


## 2.1. The Study Area

The Lake Tana (shown in Figure1) is located in Northwestern Ethiopia (latitude 10.95° and 12.78°N and longitude 36.89° and 38.25°E) with a drainage area of about 15,000 km2 [19]. The Lake Tana Basin (Figure 1), the largest lake in Ethiopia and the third largest in the Nile Basin, is found in this region. It is shared by three administrative zones of the Amhara regional state, namely, North Gondor, South Gondor, and West Gojjam directly. The climate of the area is largely controlled by the movement of the intertropical convergence zone and tropical highland monsoon, which results in a single rainy season between June and September. 70–90% of total rainfall occurs in kiremt season between June and September. The mean annual rainfall varies from 1,200 to 1,600 mm [20], and the mean actual evapotranspiration of the catchments is 773 mm.


## 2.2. Data

The main data used in this study are monthly precipitation for each year in the whole statistical period, which are collected from the National Meteorology Agency (NMA) of Ethiopia. We have used four weather stations (Table1) which spread out over the Lake Tana Basin along with a full statistical period of 31 years (1985–2015). Station selection was based on quality, long range data, and representation of various climate zones in the Basin. Some station selection criteria justify that the study of historical climate variability and change should utilize reliable data that are free of artificial trends or changes. Artifacts of measurement caused by changes in observation practice, equipment, site exposure, and location can lead to misleading results when used in trend analyses [21]. Hence, the homogeneity test of the time series of precipitation data, following reliable procedures (e.g., [22]), have been applied in this study to reduce these artifacts. We have used the total monthly average precipitation for the harmonic analysis.Table 1
Stations of the studying area, locations, elevations, and mean annual rainfall.
Station name
Latitude (degree)
Longitude (degree)
Elevation (m)
Mean annual rainfall (mm)
Bahir Dar
11.57
37.36
1800
1403.17
Debre Tabor
11.86
38.02
2706
1459.40
Gondar
12.60
37.45
1967
1069.40
Dangila
11.25
36.84
2127
1633.10


## 2.3. Harmonic Analysis

Harmonic analysis is commonly used to study periodic variations that can represent the fluctuations or variations in a time series as having arisen from adding together a series of sine and cosine functions [6, 23]. It is based on the mathematical principle that a curve, viewed as a function, may be represented by a series of trigonometric function. These trigonometric functions are “harmonic” in the sense that they are chosen to have frequencies exhibiting integer multiples of the fundamental frequency determined by the sample size (i.e., length) of the data series [24].


### 2.3.1. Representation of Seasonal and Annual Cycles by Transforming a Cosine Wave

Necessarily, harmonic function laterally shifts in order to have it match the peaks and troughs of a data series. The time series data represented byY, t=0, at origin Y, can be written as a combination of different sine and cosine signals as follows:(1)Y=y´+A1cos2πtT−Ф1+A2cos4πtT−Ф2+⋯+Ancos2nπtT−Фn,where y´ is an arithmetic average of the original data series, while the rest of the terms represent the first,second,…,nth harmonics of the period T and Ф is called the phase angle or phase shift which can then be easily interpreted as corresponding to the time of maximum harmonic function.This time shifting is most conveniently accomplished when the cosine function is used because its maximum value is achieved when the angle on which it operates is zero. Shifting the cosine function to the right by the angleФ results in a new function that is maximized at Ф.Equation (1) can be rewritten as a series:(2)Y=y´+Ccos2πtT−Ф.The argument of cosine function is zero when(3)Ф=2πtT.The function is maximized at time(4)ti=ФiT2π,where ti is the ith month and Фi is the phase shift for the corresponding month of the data series.


### 2.3.2. Estimation of the Harmonic Parameters

In order to calculate amplitude (C) and phase shift (Ф) of a periodicity for which the data comprise a perfect sine curve for a known period T, direct substitution is being applied. For representing the associated precipitation data reasonably well, better choices for C and Ф can be found [6, 23]. The easiest way to do this is by applying trigonometric identity approach.(5)cosα−Ф=cosФcosα+sinαsinФ.Forα=2πt/T and multiplying both sides by amplitude C reveals the following equation:(6)Ccos2πtT−Ф=CcosФcos2πtT+CsinФsin2πtT+Acos2πtT+Bsin2πtT,where A=CcosФ∧B=CsinФ.(7)Y=y´+Acos2πtT+Bsin2πtT,where(8)C=A2+B21/2.Fornth order of harmonics, (6) can be rewritten as follows:(9)Y=y´+∑i=1naicos2πtiT+bisin2πtiT.Similarly,(10)Ai=ai2+bi21/2,(11)Фi=tan−1biai.For the special situation where the data values are equally spaced in time with no missing values, the properties of the sine and cosine functions allow the same least-squares parameter values to be obtained more easily and efficiently using the following equation:(12)ai=2N∑t=1N−1Ycos2πtiT,bi=2N∑t=1N−1Ysin2πtiT,where i=1,2,3,…,N/2, aiandbi are coefficients of the ith harmonic given by [6, 25], and Y is the monthly mean precipitation of the ith month, and N is the number of observations.The variance of harmonic can be also obtained using Wilks [6]:(13)Vi=Ai22.The percentage of the total month-to-month variations explained by a given harmonic can be determined by forming a ratio of the square of the amplitude of that specific harmonic to the sum of all harmonics. The first harmonic can show the observed data of precipitation as a single annual cycle, whereas the second harmonic can represent the tendency towards a semiannual variation of precipitation that cannot appeared from the observed data [10]. With the same concept, the third harmonic describes more detail of the annual variation of precipitation and so on. In brief, the first harmonics show an annual cycle with the greatest amplitude as compared to other harmonics. The second harmonic indicates a semiannual change, while the third one illustrates four-month seasonal change in details [12]. The phase angle shift indicates the displacements of the maximum along the time axis, which can be represented by the x-axis. The phase angle shift can express the precipitation regime and the boundaries [26].


## 2.3.1. Representation of Seasonal and Annual Cycles by Transforming a Cosine Wave

Necessarily, harmonic function laterally shifts in order to have it match the peaks and troughs of a data series. The time series data represented byY, t=0, at origin Y, can be written as a combination of different sine and cosine signals as follows:(1)Y=y´+A1cos2πtT−Ф1+A2cos4πtT−Ф2+⋯+Ancos2nπtT−Фn,where y´ is an arithmetic average of the original data series, while the rest of the terms represent the first,second,…,nth harmonics of the period T and Ф is called the phase angle or phase shift which can then be easily interpreted as corresponding to the time of maximum harmonic function.This time shifting is most conveniently accomplished when the cosine function is used because its maximum value is achieved when the angle on which it operates is zero. Shifting the cosine function to the right by the angleФ results in a new function that is maximized at Ф.Equation (1) can be rewritten as a series:(2)Y=y´+Ccos2πtT−Ф.The argument of cosine function is zero when(3)Ф=2πtT.The function is maximized at time(4)ti=ФiT2π,where ti is the ith month and Фi is the phase shift for the corresponding month of the data series.


## 2.3.2. Estimation of the Harmonic Parameters

In order to calculate amplitude (C) and phase shift (Ф) of a periodicity for which the data comprise a perfect sine curve for a known period T, direct substitution is being applied. For representing the associated precipitation data reasonably well, better choices for C and Ф can be found [6, 23]. The easiest way to do this is by applying trigonometric identity approach.(5)cosα−Ф=cosФcosα+sinαsinФ.Forα=2πt/T and multiplying both sides by amplitude C reveals the following equation:(6)Ccos2πtT−Ф=CcosФcos2πtT+CsinФsin2πtT+Acos2πtT+Bsin2πtT,where A=CcosФ∧B=CsinФ.(7)Y=y´+Acos2πtT+Bsin2πtT,where(8)C=A2+B21/2.Fornth order of harmonics, (6) can be rewritten as follows:(9)Y=y´+∑i=1naicos2πtiT+bisin2πtiT.Similarly,(10)Ai=ai2+bi21/2,(11)Фi=tan−1biai.For the special situation where the data values are equally spaced in time with no missing values, the properties of the sine and cosine functions allow the same least-squares parameter values to be obtained more easily and efficiently using the following equation:(12)ai=2N∑t=1N−1Ycos2πtiT,bi=2N∑t=1N−1Ysin2πtiT,where i=1,2,3,…,N/2, aiandbi are coefficients of the ith harmonic given by [6, 25], and Y is the monthly mean precipitation of the ith month, and N is the number of observations.The variance of harmonic can be also obtained using Wilks [6]:(13)Vi=Ai22.The percentage of the total month-to-month variations explained by a given harmonic can be determined by forming a ratio of the square of the amplitude of that specific harmonic to the sum of all harmonics. The first harmonic can show the observed data of precipitation as a single annual cycle, whereas the second harmonic can represent the tendency towards a semiannual variation of precipitation that cannot appeared from the observed data [10]. With the same concept, the third harmonic describes more detail of the annual variation of precipitation and so on. In brief, the first harmonics show an annual cycle with the greatest amplitude as compared to other harmonics. The second harmonic indicates a semiannual change, while the third one illustrates four-month seasonal change in details [12]. The phase angle shift indicates the displacements of the maximum along the time axis, which can be represented by the x-axis. The phase angle shift can express the precipitation regime and the boundaries [26].


## 3. Result and Discussion

In order to obtain the dominant harmonics, we have calculated different parameters. Consequently, Table2 shows the calculations necessary to obtain least-square estimates for the parameters of the annual harmonic representing the Lake Tana Basin’s mean monthly precipitation, using (12). The precipitation data are shown in the column labeled Y in Table 2.Table 2
Illustration of the mechanics of using (2), (3), and (5)–(13) to estimate the parameters of a fundamental harmonic for each station.
(a) Debre Tabor
Serial number
Y
cos
2
π
t
/
12
sin
2
π
t
/
12
Y
cos
2
π
t
/
12
Y
sin
2
π
t
/
12
1
406.88
−
0.5000
−
0.8660
−203.44
−352.36
2
389.11
−
0.8660
−
0.5000
−336.97
−194.56
3
183.86
−
0.0000
−
1.0000
0.00
−183.86
4
165.34
−
1.0000
0.0000
−165.34
0.00
5
76.86
0.5000
−
0.8660
38.43
−66.56
Sum
1222.05
−1.866
−3.232
−667.32
−798.34
a
i
b
i
A
i
Ф
i
rad
V
i
P
V
i
1
−
33.91
−
58.56
67.67
1.05
2289.6
4.47
2
−90.07
−
90.99
128.03
0.79
8195.9
15.99
3
−90.07
−
121.63
151.35
0.93
11,453.2
22.35
4
−117.63
−121.63
169.21
0.80
14,315.3
27.93
5
−111.22
−132.72
173.16
0.87
14,992.2
29.26
(b) Bahir Dar
Serial number
Y
cos
2
π
t
/
12
sin
2
π
t
/
12
Y
cos
2
π
t
/
12
Y
sin
2
π
t
/
12
1
417.04
−
0.867
−
0.498
−361.57
−207.69
2
381.08
−
0.502
−
0.865
−191.30
−329.63
3
195.84
−
0.000
−
1.000
0.00
−195.84
4
184.91
−
1.000
0.000
−184.91
0.00
5
92.16
0.498
−
0.867
45.90
−79.90
Sum
1271.39
−1.871
−3.23
−691.88
−813.06
a
i
b
i
A
i
Ф
i
rad
V
i
P
V
i
1
−
60.26
−
34.62
69.50
0.52
2414.91
4.55
2
−91.14
−
89.56
127.78
0.78
8163.75
15.38
3
−91.14
−
122.20
154.26
0.93
11897.59
22.41
4
−121.95
−122.20
172.64
0.79
14902.32
28.07
5
−114.30
−135.52
177.29
0.87
15715.08
29.60
(c) Gondar
Serial number
Y
cos
2
π
t
/
12
sin
2
π
t
/
12
Y
cos
2
π
t
/
12
Y
sin
2
π
t
/
12
1
288.34
−0.8660
−0.5000
−249.70
−144.17
2
282.9
−0.5000
−0.8660
−141.45
−244.99
3
154.35
−1
0.00
−154.35
0.00
4
107.51
−0.00
−1
0.00
−107.51
5
82.02
−0.8660
0.5000
−71.03
41.01
Sum
915.12
−3.232
−1.866
−616.53
−455.66
a
i
b
i
A
i
Ф
i
rad
V
i
P
V
i
1
−41.62
−24.03
48.06
0.17
1154.83
4.22
2
−65.20
−64.86
91.97
0.78
4228.93
15.46
3
−90.93
−64.86
111.69
0.62
6237.54
22.81
4
−90.93
−82.78
122.97
0.74
7560.40
27.65
5
−102.77
−75.94
127.78
0.64
8164.28
29.86
(d) Dangila
Serial number
Y
cos
2
π
t
/
12
sin
2
π
t
/
12
Y
cos
2
π
t
/
12
Y
sin
2
π
t
/
12
1
382.32
−0.5000
−0.8660
−191.16
−331.09
2
362.28
−0.8660
−0.5000
−313.73
−181.14
3
261.57
−1
0.00
−261.57
0.00
4
254.58
0.00
−1
0.00
−254.58
5
150.99
−0.8660
0.5000
−130.76
75.50
Sum
1411.74
−3.232
−1.866
−897.22
−691.31
a
i
b
i
A
i
Ф
i
rad
V
i
P
V
i
1
−31.86
−55.18
63.72
1.05
2029.95
3.68
2
−84.15
−85.37
119.87
0.79
7184.63
13.03
3
−127.75
−85.37
153.65
0.59
11,804.05
21.40
4
−127.75
−127.80
180.70
0.79
16,326.45
29.60
5
−149.46
−115.22
188.72
0.66
17,806.97
32.29Figure2 illustrates the foregoing procedure using the mean monthly precipitation (mm) for 1985–2015 in the Lake Tana Basin. This figure is simply a plot of the 12 data points, with t=1 indicating January, t=2 indicating February, and so on. The monthly annual average precipitations, which are represented by dotted points, are 122.06,116.93,89.11, and 136.08mm for Debre Tabor, Bahir Dar, Gondar, and Dangila, respectively. The lowest value of precipitation over Gondar station is due to the surrounding areas’ dryness which implies that the wind flow from the dry land does not bring enough amount of moisture. The consecutive occurrence of wet and dry sequences can alter the variability of precipitation distribution which depends on regional-scale predictors such as humidity and atmospheric pressure as Wilby et al. [27] found. In addition to this, we suppose that, land cover of the area will also be the possible reason, but it needs farther study. So, the rainfall of Gondar has been (highly) variable in time and concentrated in summer season (at t=7 (July) and t=8 (August), Figure 2), considering the above point. Other studies also conclude that the rainfall of mountainous areas is highly variable from time to time and from station to station [28, 29]. Those all bring high interannual and semiannual variability of precipitation as the percentage of variation as shown in Figure 4.Figure 2
The monthly precipitation (mm) averaged over 1985–2015 for the four stations: (a) Bahir Dar, (b) Dangila, (c) Debre Tabor, and (d) Gondar. Circular points show climatic mean, and horizontal dotted lines show annual mean.The monthly maximum precipitation values are406.88 (August), 417.04 (July), 288.34 (July), and 382.32 (August) for Debre Tabor, Bahir Dar, Gondar, and Dangila, respectively. These data appear to be at least approximately sinusoidal, executing a single full-cycle over the course of the 12 months. The precipitation at Bahir Dar station shows the annual peak of 417.04 mm in July and the annual valley value of 1.66 mm in January (Figure 2 for all the other stations). Generally, the annual maximum precipitation for all stations occurred in July and August which are the months of main rainy season for the region and most parts of Ethiopia [20].The blue curve at Figure3 shows the function lifted to the level of the average annual precipitation and stretched so that its range is similar to that of the data series (2), with ϕ=0. The stretching has been done only approximately, by choosing the amplitude C to be the maximum precipitation value between the January and December. Simply it is a cosine function with the argument transformed so that it executes one full cycle in 12 months. It is obviously a poor representation of the data because the representation of the curve is opposite to the normal data distribution which is represented by the circular points. Hence, the cosine curve needs to be shifted downward to line up well with the data. The five harmonics plotted downward are used to overcome this poor representation of the data. The maximum in the curve can be made to occur at t=7 months (July) and t=8 months (August) by introducing the phase shift, using (3), of ϕ=72π/12=7/6π and ϕ=82π/12=4/3π.Figure 3
Illustration of the approximate matching of a cosine function to a data detail about the lines. (a) Bahir Dar, (b) Dangila, (c) Debre Tabor, and (d) Gondar.From Figure3, we can also see that the first five harmonics represent the data distribution appropriately. Accordingly, the first two harmonics green (1st) and red (2nd) curves show the dominant harmonics which occur in July (for Bahir Dar and Gondar) and August (for Debre Tabor and Dangila), respectively. The amplitudes for the first harmonic are 173.42,177.93,127.77, and 188.78mm for Debre Tabor, Bahir Dar, Gondar, and Dangila, respectively.Figure4 shows the percentage of variance for each harmonic. The percentages of variance for the first harmonics are 3.68, 4.22, 4.48, and 4.55 in Dangila, Gondar, Debre Tabor, and Bahir Dar, respectively, which represents low variations that can assures that there is difference in climatic regime in the basin. The percentages of variance for the second harmonic are 13.03, 15.38, 15.46, and 15.99 in Dangila, Bahir Dar, Gondar, and Debre Tabor, respectively. Similarly, the percentage of variance for third harmonic are 21.40, 22.35, 22.41, and 22.81 in Dangila, Debre Tabor, Bahir Dar, and Gondar, respectively, which exactly agrees with the definitions of percentage of variance; low percentage of variance for first harmonic necessarily increase the percentage of other harmonics [6, 10]. From the percentage of variance, the second and third harmonics have higher percentage than the first harmonic (Table 2 and Figure 4) which demonstrates that different parts of the basin do have different rainfall regimes.Figure 4
Percentage of variance for all stations (x-axis: 1 = Bahir Dar, 2 = Dangila, 3 = Debre Tabor, and 4 = Gondar).As it is mentioned in the above, low percentages of the first harmonic variance require more contribution from high percentages of other harmonics (like semiannual or seasonal variations). The last harmonics indicates poor representation of the time series data.


## 4. Conclusion

We have presented the rainfall climatology of Lake Tana Basin using the harmonic analysis technique. This method is commonly used to study periodic variations that can represent the fluctuations or variations in a time series for diverse topographical features as having arisen from adding together a series of sine and cosine functions. Precipitation features are explained in terms of percentage of variances, amplitudes, and phase angle.Climatologically, these regions are affected by various well-known weather conditions, such as ITCZ, ENSO, and so on, which affect different regions at specific times of the year. The percentages of variations are related to the effect of these conditions that usually affect each region separately. The maximum precipitation for all stations occurred in July and August. The annual average precipitation of Dangila was 136.08 mm which was the highest record than the other stations, and the average maximum monthly precipitation (417.04 mm) was recorded in Bahir Dar during July. The amplitudes for the first harmonic are173.42,177.93,127.77, and 188.78mm for Debre Tabor, Bahir Dar, Gondar, and Dangila, respectively. The percentage of variance shows lower value in the first harmonic than others, implying that there is difference in the climatic regime in the Lake Tana Basin. Low percentages of the first harmonic variance mean that it requires more contribution from high percentages of other harmonics (like semiannual or seasonal variations).In all stations, the original time series precipitation has been represented fairly by the first five harmonics. But the number of stations used in this study is less in number to represent the whole basin. Hence, we recommend extension of the number of stations and spatial coverage to incorporate the Upper Blue Nile Basin and beyond, performing the proposed data quality analysis and calculation of the harmonic parameters for further studies.



---

*Source: 1598195-2018-06-07.xml*