PublishedPaper HiJASE KBasnet
PublishedPaper HiJASE KBasnet
net/publication/347739674
CITATIONS READS
2 245
4 authors, including:
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Keshav Basnet on 14 January 2021.
Abstract
Climate change is considered as the most critical global challenge of the century. Unusual precipitation
pattern and outflanking of hydraulic structures frequently reported these days in the Kaski district of
Nepal. This research aimed to analyze the trend of climate change in the lower Kaski using annual
and seasonal weather data (2001-2017) of five meteorological stations out of which three on the upper
tropical region, one on the sub-tropical region and the other one on the temperate region. Individual
trend analysis was performed for rainfall and temperature patterns using Mann-Kendell test. Trend
analysis for seasonal average precipitation shows an increasing trend for winter and pre-monsoon
seasons and a decreasing trend for monsoon and post-monsoon seasons. In contrast, overall annual
precipitation shows a decreasing trend. An increasing trend was found for maximum seasonal
temperature for all seasons except pre-monsoon season. The seasonal minimum temperature shows a
growing trend and overall annual trends for both the maximum and minimum temperatures were also
increasing. All of these trends are the evidence of climatic changes that are happening over time.
Additionally, a comparative study was conducted between a meteorological station and a locally
established manual station. The differences in the numeric values of cumulative rainfall with
comparable precipitation readings suggest the need to take into account local station data for design,
construction, and planning of hydraulic structures.
Keywords: Global warming, Local rain gauge station, Mann Kendall test, Rainfall trend, Temperature
trend.
11
Himalayan Journal of Applied Science and Engineering (HiJASE), Vol. 1, Issue 1, No. 11, 2020
Northern Hemisphere [1]. The warming trend in required to explain the processes is also limited.
Nepal after 1997 has ranged from 0.06 °C to 0.12 With few Monitoring stations (just 280 across the
°C per year in most of the middle mountains and entire country), the hydro-meteorological data has
Himalayan regions; this trend has been less than been collected only since the late 1960s [5]. The
0.03 °C per year in the Siwalik and Terai regions observed climate trend analysis of Nepal is
[2]. Predictions are constantly being made to available from 1971 to 2014 [6]. However, there
understand what will happen based on observed is no study on climate change trend primarily
climate shifts over the past decades and projected focused on the Kaski district comparing the trend
greenhouse gas emissions. Changing precipitation of various climatic zones.
and temperature analysis trends are useful in
The study area of this research is based on the
engineering, environmental, forestry, and soil
Kaski district of Gandaki Province, Nepal (Fig. 1).
conservation studies. This trend analysis is
It is located at latitude 28°18′19.08″ North and
necessary for proper engineering designs of
longitude 84°04′37.20″ East. The altitude of the
hydrologic structures like dams, bridges,
Kaski district, with Pokhara as headquarter,
irrigation canals, sewerage, etc. Addisu et al. [3]
ranges from lower 450 masl to upper 8091 masl in
conducted a time series trend analysis of rainfall
the Himalaya range. It has a Metropolitan City, 4
and temperature in the Ethiopia's lake Tana Sub-
Rural Municipalities, and 3 Electoral Sectors.
basin. Similarly, Shrestha [4] studied the variation
Basnet and Acharya [7] reported that Kaski is the
of summer monsoon rainfall over Nepal and its
area of high intensity of rainfall throughout the
relation to the southern oscillation index.
year in Nepal, and about 80% of rainfall occurs
Nepal has amazingly diverse geological features typically in four months; June, July, August, and
and climatic conditions throughout the country. September. With a focus on our study area, the
The hydrology of Nepal is mainly fed by the South unusual precipitation patterns are frequently being
Asian Monsoon System (SAM). Still, the reported these days. Besides that, the overflow of
relationship between the timing, volume of the drain, outflanking of hydraulic structures such
monsoon rainfall, and the mountain landscape are as bridges, dams, river training works, etc. are also
poorly understood [5]. The altitude variation over observed recurrently. This incidence is not a mere
a short distance limits our ability to explain the coincidence, but the result of climatic changes
precipitation pattern in Nepal. The data set happening over the course of time.
This study's main objective is to analyze the Based on this study, out of the seven climatic
climate change trend in the Kaski district of zones of Kaski district (see Table 1), only three
Nepal. For that, the trend of rainfall and climatic zones are taken for historical trend
temperature in different climatic zones were analysis. The results are based on Sen's slope and
analyzed, and also the need for local rain gauge the established rain gauge is of manual type.
stations in construction activities was assessed.
K. Basnet et al. 12
Himalayan Journal of Applied Science and Engineering (HiJASE), Vol. 1, Issue 1, No. 11, 2020
This study's secondary objective is to compare the rainfall data relevant to the study area or proposed
rainfall data collected from a manually established site (e.g., while designing storm water drain, see
rain gauge device with the precipitation data of a [8]). For this purpose, a rain gauge was
DHM station. Researchers and hydraulic established manually to collect the local rainfall
engineers in Nepal facing the problem of lack of data [9].
Figure 2: Research framework; a) for trend analysis and b) for need assessment.
K. Basnet et al. 13
Himalayan Journal of Applied Science and Engineering (HiJASE), Vol. 1, Issue 1, No. 11, 2020
Trend analysis is the widespread practice of trends in series of environmental data, climate
collecting information and identifying a pattern. data or hydrological data [12]. The null
The detection, estimation and prediction of trends hypothesis, H0, is that the data come from a
and associated statistical and physical population with independent realizations and are
significance are important aspects of climate identically distributed. In statistics, the Kendall
research [10, 11]. Given a time series of (say) rank correlation coefficient, commonly referred
temperatures, the trend is the rate at which to as Kendall's tau coefficient (after the Greek
temperature changes over a time period which letter τ), is a statistic used to measure the ordinal
may be linear or non-linear. Simple linear association between two measured quantities.
regression is most commonly used to estimate the Being a non-parametric test, no necessity of
linear trend (slope) and statistical significance. normally distributed data and having low
The non-parametric (i.e. distribution free) Mann- sensitivity to abrupt breaks due to
Kendall (M-K) test can also use to assess inhomogeneous time series, this statistical
monotonic trend (linear or non-linear) method is preferred for climatologic trend time
significance. The non-parametric Mann-Kendall series analysis.
test is commonly employed to detect monotonic
Table 2: Hydro-meteorological stations in the lower Kaski district. [Source: DHM, Nepal]
Latitude Longitude Elevation
S. No. Name of station Station type
(°N) (°E) (m)
1 Pokhara Airport Aeronautical 28.13 84.00 827
2 Malepatan Agrometeorology 28.07 84.07 856
3 Lumle Agrometeorology 28.18 83.48 1740
4 Begnas Climatology 28.11 84.08 669
5 Panchase Climatology 28.13 83.47 2493
The Mann-Kendall S Statistic is computed [13] as normally distributed with the mean and variance
follows: (σ2) as:
∑
𝑆 ∑ ∑ 𝑠𝑖𝑔𝑛 𝑇 𝑇 (1) 𝜎
(3)
K. Basnet et al. 14
Himalayan Journal of Applied Science and Engineering (HiJASE), Vol. 1, Issue 1, No. 11, 2020
relationship between the two variables, and is Autocorrelation or serial correlation is defined as
carried out on the ranks of the data. For each the correlation of a variable with itself over
variable separately, the values are put in order and successive time intervals, prior to testing for
numbered, 1 for the lowest value, 2 for the next trends. It increases the chances of detecting
lowest and so on. Kendall's tau takes values of significant trends even if they are absent and vice
range [-1,1]. A positive correlation indicates that versa. Hence, it is essential to be considered in the
the ranks of both variables increase together and time series analysis. Autocorrelation between the
a negative correlation indicates that as the rank of ranks of the data can be calculated as suggested
one variable increases, the other decreases [15]. by Hamed and Rao (1998) after removing the
apparent trend [16]. The adjusted variance is
given as:
𝟏 𝑵
𝑽𝑨𝑹 𝑺 𝑵 𝑵 𝟏 𝟐𝑵 𝟓 (5)
𝟏𝟖 𝑵𝑺∗
𝑵 𝟐
𝟏 ∑𝒑𝒊 𝟏 𝑵 𝒊 𝑵 𝑰 𝟏 𝑵 𝒊 𝟐 𝑷𝒔 𝒊 (6)
𝑵𝑺∗ 𝑵 𝑵 𝟏 𝑵 𝟐
Malepatan 30310.8 7558.7 (July) 80.8 (Dec) 4251.5 (2007) 2982.5 (2006)
Lumle 45348.8 11535.0 (July) 93.3 (Dec) 6310.2 (2003) 4294.7 (2006)
K. Basnet et al. 15
Himalayan Journal of Applied Science and Engineering (HiJASE), Vol. 1, Issue 1, No. 11, 2020
K. Basnet et al. 16
Himalayan Journal of Applied Science and Engineering (HiJASE), Vol. 1, Issue 1, No. 11, 2020
As presented in Tables 5 and 6, in the winter range, 6.71 mm per season in the sub-tropical
season, the average precipitation shows an range and 5.29 mm per season in the temperate
increasing trend of 1.27 mm per season at the range.
upper tropical range, 1.41 mm per season at the
3.2.2. Annual Precipitation Trend:
sub-tropical range and 3.36 mm per season at the
temperaure range. In pre-monsoon season, the As shown in Tables 7 and 8, in the upper tropical
average precipitation shows an increasing trend range, the average annual precipitation is
of 0.01 mm per season at the upper tropical range, decreasing at the rate of 2.22 mm per year, 8.13
15.21 mm per season at the temperature range and mm per the year in sub-tropical range and is
a decreasing trend of 0.60 mm per season in the increasing at the rate of 13.47 mm per year in
sub-tropical range. In the monsoon season, the temperate range. More moisture can be hold by a
average precipitation shows an increasing trend warmer atmosphere, and water vapor increases by
of 2.84 mm per season at the sub-tropical range 7 % for every degree centigrade of warming
and 142.38 mm per season at temperate range but globally [18]. Although, the translation of this
a decreasing trend of 71.58 mm per season at the phenonmenon to global changes in precipitation
upper tropical range. In the post-monsoon season, is quite unclear, the total volume of precipitation
the average precipitation showed a decreasing is likely to increase by 1-2 % per degree of
trend of 23.89 mm per season at the upper tropical warming [18].
3.3. Maximum Temperature Trend: season in the upper tropical range, 0.19 °C per
season in the sub-tropical range and 0.27 °C per
3.3.1. Seasonal Maximum Temperature
season in temperature range. In the monsoon
Trend:
season, the maximum temperature is increasing
Seasonal maximum temperature trend analysis by 0.07 °C per season in the upper tropical range
was performed for the two periods which are but decreasing at the rate of 0.19 °C per season in
shown in Table 9 and Table 10. In the winter sub-tropical range and 0.17 °C per season in
season, maximum seasonal temperature showed temperate range. In the post-monsoon season, the
an increasing trend of 0.17 °C per season in the average maximum temperature is increasing at
upper tropical range and 0.06 °C per season in the rate of 0.1 °C per season in the upper tropical
temperate- range but the decreasing trend of 0.05 range, decreasing at the rate of 0.02 °C per season
°C per season in sub-tropical range. Similarly, in in the sub-tropical range and increasing at the rate
the pre-monsoon season, the maximum of 0.05 °C per season in the temperate range.
temperature decreases at the rate of 0.09 °C per
K. Basnet et al. 17
Himalayan Journal of Applied Science and Engineering (HiJASE), Vol. 1, Issue 1, No. 11, 2020
3.3.2. Annual Maximum Temperature Trend: range. In the sub-tropical range, the average
annual maximum temperature is decreasing at the
As Table 11 and Table 12 represent, the average
rate of 0.06 °C per year. In the temperate range,
annual maximum temperature is increasing at the
the average annual maximum temperature is
rate of 0.082 °C per year in the upper tropical
decreasing at the rate of 0.066 °C per year.
K. Basnet et al. 18
Himalayan Journal of Applied Science and Engineering (HiJASE), Vol. 1, Issue 1, No. 11, 2020
Table 12: Annual maximum temperature trend (2010-2017).
S. No. Name of station τ p-value Sen's Slope Trend
1 Pokhara Airport 0.17 0.39 0.08 Increasing
2 Malepatan 0.23 0.25 0.03 Increasing
3 Begnas 0.32 0.29 0.16 Increasing
4 Lumle -0.31 0.43 -0.17 Decreasing
5 Panchase -0.05 0.64 -0.07 Decreasing
K. Basnet et al. 19
Himalayan Journal of Applied Science and Engineering (HiJASE), Vol. 1, Issue 1, No. 11, 2020
precipitation collected in the local rain gauge was established by the Department of Hydrology was
2129.4 mm (June-August, 2018), and the amount 1963.4 mm (June-August, 2018).
of precipitation collected by the rain gauge
Table 14: Seasonal minimum temperature trend (2010-2017).
S. Name of p- Sen's p- Sen's
τ Trend τ Trend
No. Stations value Slope value Slope
Winter Pre- Monsoon Season
Pokhara
1 0.64 0.04 0.28 Increasing 0.07 0.90 0.05 Increasing
Airport
2 Malepatan 0.29 0.39 0.25 Increasing 0.07 0.90 0.08 Increasing
3 Begnas 0.50 0.11 0.19 Increasing 0.14 0.71 0.10 Increasing
4 Lumle 0.62 0.07 0.36 Increasing 0.33 0.37 0.25 Increasing
5 Panchase 0.29 0.39 0.13 Increasing -0.21 0.54 -0.10 Decreasing
Monsoon Season Post Monsoon Season
Pokhara
1 0.47 0.14 0.08 Increasing 0.50 0.11 0.12 Increasing
Airport
2 Malepatan 0.36 0.27 0.06 Increasing -0.26 0.45 -0.05 Decreasing
3 Begnas 0.14 0.71 0.03 Increasing -0.07 0.90 -0.03 Decreasing
4 Lumle 0.14 0.76 0.03 Increasing 0.68 0.05 0.18 Increasing
5 Panchase 0.29 0.39 0.17 Increasing -0.04 1.00 0.00 No trend
Table 16: Comparision of established and DHM As the stations were not located exactly at the
rain gauge stations. same place, slightly different values of
Cumulative Cumulative precipitations are observed in those stations
precipitation precipitation
though it was similar for the first few days.
Month at DHM rain at local rain
gauge station gauge station Therefore, the local rain gauge station's
(mm) (mm) precipitations values obviously differ with the
June 578.5 671.5 precipitation values of DHM station if one keeps
July 831.5 708.0 recording for a more extended period than the
August 553.4 750.0 three months of record period for this study.
Total
1963.4 2129.4
Precipitation
K. Basnet et al. 20
Himalayan Journal of Applied Science and Engineering (HiJASE), Vol. 1, Issue 1, No. 11, 2020
gauge station can help to obtain in-situ
precipitation data for designing hydraulic
structures. The World Meteorological
Organization (WMO) suggests the minimum area
for one rain gauge station to be less than 25 km2
for the mountainous region of the world's
temperate and Mediterranean and tropical zone.
But the station at Malepatan and Pokhara Airport
is nearer and also the station at Lekhnath and
Begnas are under the area of 25 km2, yet the
precipitation varies spatially. Hence, in the
districts like Kaski, where there is a heavy
downpour and massive variation in precipitation
Figure 3: Comparison of cumulative precipitation (mm)
patterns within such a short distance, the local
of local and DHM rain gauge stations.
rain gauge station readings need to be considered
4. Conclusion: for economic hydrological constructions and to
avoid the failures of hydraulic structures.
The climate change study was performed for three _______________________________________
climatic zones of Kaski, Nepal, namely; upper
tropical region, sub-tropical region and temperate References:
region. The main goal of this research was to
[1] IPCC, "Climate change 2007: impacts,
analyze the climate changing trend by individual
adaptation and vulnerability”,
trend analysis of rainfall and temperature patterns
Contribution of Working Group II to the
by using the Mann-Kendell test. For the lower
Fourth Assessment Report of the
Kaski district, it is found that the seasonal
Intergovernmental Panel on Climate
precipitation in winter and pre-monsoon seasons
Change, Cambridge University Press,
are increasing, and the seasonal precipitation in
Cambridge, UK, 2007.
monsoon and post-monsoon seasons are
decreasing. The overall annual precipitation trend [2] A. B. Shrestha, C. P. Wake, P. A.
is decreasing. The seasonal maximum Mayewski, and J. E. Dibb, “Maximum
temperature trend of winter, monsoon and post- temperature trends in the Himalaya and
monsoon seasons are increasing and that of pre- its vicinity: an analysis based on
monsoon season is decreasing. The overall annual temperature records from Nepal for the
maximum temperature trend is increasing. The period 1971–94”, Journal of Climate,
seasonal minimum temperature trend of temperate vol. 12, pp. 2775-2786, 1999.
region is decreasing whereas other two regions are
[3] S. Addisu, Y. G. Sealassie, G. Fissha, and
increasing and the overall annual minimum
B. Gedif, “Time series trend analysis of
temperature is also increasing. However this study
temperature and rainfall in lake Tana
is limited due to the use of insufficient data for
sub-basin, Ethiopia", Environmental
trend analysis of climate change. Hence, it can be
System Research, DOI: 10.1186/s40068-
further extended in future to use climate data of
015-0051-0, 2015.
more than thirty years. Furthermore detail
assessment can be performed using the daily and [4] M. L. Shrestha, "Interannual variation of
monthly weather data as this study focused only summer monsoon rainfall over Nepal and
on annual and seasonal trends. its relation to southern oscillation index",
Meteorology and Atmospheric Physics,
In addition, a local rain gauge was established
vol. 75, pp. 21-28, 2000.
manually and compared to the rainfall data with a
DHM station. Precipitation pattern from two [5] A. Dixit., " Climate Change in Nepal:
methods was found similar but with the Impacts and Adaptive Strategies".,
difference in numeric values of cumulative Institution for Social and Environmental
rainfall. It suggests that a locally established rain Transition-Nepal, World Resources
Report., https://www.wri.org/our-
K. Basnet et al. 21
Himalayan Journal of Applied Science and Engineering (HiJASE), Vol. 1, Issue 1, No. 11, 2020
work/project/world-resources- International Journal of Environmental
report/climate-change-nepal-impacts- Sciences, vol. 6, no. 6, pp. 954-969,
and-adaptive-strategies. 2016.
[6] Government of Nepal, Ministry of [12] R. Ilaboya, “Non-parametric Mann-
Population and Environment, Kendall test statistics for rainfall trend
Department of Hydrology and analysis in some selected states within
Meteorology, "Observed climate trend the coastal region of Nigeria", Journal of
analysis of Nepal (1971-2014)”, 2017. Civil, Construction and Environmental
Engineering, vol. 3, pp. 17-28, 2018.
[7] K. Basnet and D. Acharya, “Flood
analysis at Ramghat, Pokhara, Nepal [13] H. Motiee, E. McBean., “An Assessment
using HEC-RAS”, Technical Journal, of Long Term Trends in Hydrologic
vol. 1, no. 1, pp. 41-53, 2019. Components and Implications for Water
https://doi.org/10.3126/tj.v1i1.27591. Levels in Lake Superior “, Hydrology
Research, 40.6, 564-579, 2009.
[8] S. Khadka and K. Basnet, “Storm water
management of Barahi Chowk area, [14] K. Drapela, I. Drapelova., "Application
Lakeside, Pokhara, Nepal using of Mann-Kendall test and the Sen’s slope
SWMM”, Proceedings of 2nd estimates for trend detection in
International Conference on Engineering deposition data from Bílý Kříž (Beskydy
& Technology, Lalitpur, Nepal, vol. 2, Mts., the Czech Republic) 1997–2010",
pp. 320-325, 2019. Beskdy Mendel University in Brno 4 (2),
133–146, 2011.
[9] K. Basnet and M. Neupane, “Storm water
drainage design based on hydrological [15] P. Llukan., "Kendall’s Tau".
analysis: a case study on Lamachaur International Encyclopedia of Statistical
catchment area, Pokhara, Nepal”, Science., DOI:
Oodbodhan: A Journal of TUTA, vol. 5, https://doi.org/10.1007/978-3-642-
no. 1, pp. 100-109, 2018. 04898-2_324, 2014.
[10] A. Basarir, H. Arman, S. Hussein, A. [16] K. H. Hamed and A. Rao, “A modified
Murad, A. Aldahan, M. A. Al-Abri, Mann-Kendall trend test for
"Trend detection in annual temperature autocorrelated data”, Journal of
and precipitation using Mann–Kendall Hydrology, vol. 204, pp.182-196, 1998.
test - a case study to assess climate
[17] T. Sinha, K. Cherkauer, "Time series
change in Abu Dhabi, United Arab
analysis of soil freeze and thaw processes
Emirates", Proceedings of 3rd
in Indiana", Journal of Hydrometrology,
International Sustainable Buildings
https://doi.org/10.1175/2008JHM934.1,
Symposium, Lecture Notes in Civil
2008.
Engineering, vol. 7. Springer, Cham,
2017. [18] Will Climate Change Lead To More
Rainfall And Extreme Floods In India?,
[11] S. S. Chinchorkar, D. M. Bhavin Ram,
Scientific India,
Paradava, and M. M. Trivedi, "A case
https://scind.org/1294/Environment/will
study of temperature and rainfall trends
-climate-change-lead-to-more-rainfall-
using Mann-Kendall test in Saurashtra
and-extreme-floods-in-india.html , 2018.
Region (Junagadh) of Gujaart, India",
K. Basnet et al. 22