# Trend characteristics of rainy days and evaporation at a tropical rainforest region in East Malaysia, Borneo.

IntroductionEfforts to understand the impacts of climate change at a regional or even local scale, in this case a river basin, has acquired greater importance during the last few decades due to its increased adverse effects on hydrological, agricultural, industrial as well as economic sectors. The effects of the Earth's changing climate are demonstrated through varying hydrometeorological characteristics at a spatial and temporal level (Frich et al., 2002). Spatial and temporal variability of the impacts of climate change are not only limited to the amount of rain falland an increase in temperature but also include increased rates of rainfall in some regions, frequent fash foods, continuous droughts (meteorological as well as agriculture), decreasing groundwater levels, cyclones, melting of glaciers and so on (Raziei et al., 2005; IPCC, 2007; Sharad & Vijay 2012; Zafor et al., 2016a). All these phenomena have been well studied during the last three decades in an effort to reduce the crumbling effects of climate change (Crisci et al., 2002; Walther et al., 2002; Benestad, 2013; Paparrizos & Matzarakis, 2017; Sobral et al., 2019; Parry et al., 2020). This data collection has been improved due to the development of sophisticated climatic parameter measuring and monitoring instruments as well as the increased availability of satellite-based observations and measurements (Schloss et al., 1996; Rodell et al., 2004; Cretauxet et al., 2011; Rhee & Im, 2017; Asfaw et al., 2018). Along with this increased availability of higher resolution gridded climate data sets increased the ease of regional scale statistical analysis of climatic parameters (Abatzoglou, 2013; Cecil et al., 2014; Dumitrescu & Birsan, 2015; Akinsanola et al., 2017; Serrano-Notivoli et al., 2018; Simpson & McCarthy, 2019; Singh & Xiaosheng, 2019; Hassan et al., 2020; Xu et al., 2020; Yao et al., 2020). Among the observations and research carried out to understand the effects of climate change, statistical analysis of rainfall trends has attained wider importance at a global as well as local scale (Yue & Wang, 2004; Shadmani et al., 2012; Ambuan et al.,2013; Nyatuame et al., 2014; Ahmad et al.,2015; Zhang et al., 2015; Teodoro et al., 2016; Gavrilov et al., 2018; Sayyad et al., 2019; Sobral et al., 2020; Wang et al., 2020). Along with rainfall trend analysis, the climatic parameters such as air temperature, evaporation and humidity have also been analyzed at various scales through statistical measurements (Santhe et al.,2012; Mikkonen, et al., 2015; Gebremicael et al., 2017; Livada et al., 2019; Behzadi et al., 2020). Global scale studies on the effect of climate change in rainfall indicate an increase of 0.5 - 1% in decadal rainfall and a 2% increase in annual average regional rainfall during the 20 (th)century (Mosmann et al., 2004; Xu et al., 2003; 2005; Yu et al., 2006; IPCC, 2007).

Global-scale studies will help to improve a broader understanding of the effects of climate change for the planet in general but to better understand and develop management strategies, more detailed local as well as regional studies are necessary. In order to estimate the effects of climate change, most researchers have studied trends in monthly, seasonal and annual rainfall using historical datasets and by applying a range of scientific methods irrespective of climatic conditions (Lau & Wu, 2007; Allan et al., 2010; Uduak et al., 2012; Zafor et al., 2016b; Asfaw et al., 2018; Sa'adi et al., 2019; Wang et al., 2020). Most recent studies have revealed increasing, decreasing or no trend in rainfall amounts indicating the influence of climate change over rainfall (in the case of regions showing increasing or decreasing trends) or no influence (regions showing no trends). A few studies have also reported statistical trends in rainy days, evaporation rates, air temperature and extreme rainfall events correlated with climate change scenarios (Xu et al., 2006; Burn & Hesch, 2007; Bandyopadhya et al., 2009; Fu et al., 2009; Kumar & Jain, 2010, 2011; Abtew et al., 2011; Jhajharia et al., 2012; Soltani et al., 2012; Kumar et al., 2012; Madhu et al., 2015;Soro et al., 2016; Djaman et al., 2017; Oliveira et al., 2017; Toride et al., 2018; Yang et al., 2019; Behzadi et al., 2020). The statistical trend results provide evidence for climate change effects on parameters analyzed and its contribution to regional climatic changes.

Numerous studies have been reported from Malaysia related to determining statistical rainfall trends using Mann Kendall, Spearman's Rho, Sen's slope and linear regression techniques and have been concentrated in river basins in West Malaysia (Suhaila et al., 2008, 2010; Deni et al., 2010; Zin et al., 2010; Suhaila & Jemain, 2012; Palizdan et al., 2014, 2015, 2017; CheRos et al., 2016; Suhaila & Yusop, 2018;Yatim et al., 2019; Pour et al., 2020). When looking at the East Malaysian region, especially Sarawak, very few local or regional reports are available regarding the statistical trend determination of rainfall (Manton et al., 2001; Hua et al., 2013; Krishnan et al., 2018, 2019a,b; Sa'adi et al., 2017a, 2019). Further regional scale climatic model for the southeast Asian region (including the study area) indicates an increasing trend in rainfall along with temperature and also suggesting extreme rainfall in very localized nature and increase in dry areas in upper river basins in East Malaysia (Christensen et al., 2007; Chotamonsak et al., 2011; Sa'adi et al., 2017b; Tahir et al., 2018). Though these studies reported a significant increase in rainfall in the region, no attempts have been made to assess the characteristics of distribution and trend patterns of annual rainy or non-rainy days and evaporation from the region. Therefore, the present study has been framed with the objective of analyzing statistical trends in annual rainy days, non-rainy days and monthly evaporation rates at a pre-determined region in northern Borneo. The Limbang River Basin (LRB), which possesses the unique characteristics of an equatorial tropical river basin, was selected for a pilot study to determine trends in rainy as well as non-rainy days and evaporation by applying well-defined and widely used statistical techniques. The findings of the present study will help researchers to detect the impacts of climate change by applying various hydro-meteorological parameters, especially those from tropical regions. This is why these regions are significantly influenced by regional as well as local climate conditions and which often result in low agricultural productivity, forest deterioration and loss of biodiversity (Thomas et al., 2004; Thompson et al., 2009; Piao et al., 2010; Balint et al., 2011). Currently, most researchers have concentrated on determining trends in rainfall and temperature (Krishnan et al., 2019a, b; Sa'adi et al., 2019). However, information related to trends in evaporation rates as well as rainy and non-rainy days is a basic pre-requisite for the planning and development of major hydroelectric projects and large-scale agricultural schemes (Jana et al., 2016; Eymen & Koylu, 2019; Liu et al., 2019; Mandale et al., 2019; Ruiz-Alvarez et al., 2019; Ghumman et al., 2020).

Study area

Borneo, the third largest island in the world with a tropical equatorial rainfall regime is located in South East Asia. The Limbang River Basin (LRB) is one of the largest river basins in the state of Sarawak, Malaysia and covers a total area of 3,687 [km.sup.2] , located between N latitudes [03.sup.0] 47' 7" to [04.sup.0] 54' and E longitudes [114.sup.0] 46' 7" to [115.sup.0] 31' 30" (Figure 1). The Limbang River flows a distance of 196 km before it discharges into the South China Sea, at Brunei Bay. The LRB exhibits varied terrain characteristics and a highly undulating nature with an elevation ranging from sea level to 2,400 m above mean sea level. The majority of the LRB is covered with primary and secondary forests (more than 90%) but also includes oil palm plantations and mixed agriculture. Sarawak receives rainfall of various amounts throughout the year in two monsoon seasons, which determine the climatic conditions of the study area. The monsoon seasons in this region last from May to September (Southwest Monsoon) and from November to March (Northeast Monsoon) but these are separated with two short inter-monsoon months i.e. April and October. The Northeast monsoon is the wettest period and the Southwest monsoon comparatively drier (DID, 2015). The study area (LRB) receives an annual average rainfall of 3,851 mm with monthly mean rainfall exceeding 250 mm. Due to the high amount of rainfall and the peculiar relief characteristics of the research site, lower areas in the LRB experience repeated severe fash fooding during the peak monsoon seasons. The average daily temperature ranges from 25[degrees]C to 35[degrees]C during the day and 18[degrees]C to 21[degrees]C at night with the humidity varying from 70 to 90 percent (Tuh et al., 2013).The major human development activity in the study region is timber harvesting, which alters the terrain conditions and induces severe soil erosion during heavy rainfall as well as affecting the climatic conditions of this region (Vijith et al., 2018).

Materials and methods

Before collecting the necessary hydrometeorological information, the available data were obtained from the Department of Irrigation and Drainage (DID), Sarawak, Malaysia. It was noted that, though a number of rain gauging stations had been installed in the LRB, only three rain gauge stations, the Limbang DID, Ukong and Long Napir possessed long-term measurements, which is essential for statistical trend determination. The DID measured rainfall using manual as well as automated weather stations and the rate of evaporation was measured from a pan by assessing the change in the level of its free water from the evaporation pan with a height of 254 mm and a diameter of 1,206 mm (DID, 2017). These stations represent different terrain (physiographic) and conditions (coastal, midland and highland). Daily rainfall and evaporation data were collected from the DID and processed. The general characteristics of the rain gauge stations and the data sets collected are given in Table 1 and the locations of the rain gauges as well as the evaporation measuring stations are shown in Figure 1. The time span of the meteorological data collected ranged from 35 to 46 years (1970 to 2015) and the percentage of the missing data ranged from 4 to 17%. The elevation of the rain gauge stations varied from 7 m (Limbang DID), 22 m (Ukong) and 147 m (Long Napir). Among the rain gauge stations, only Limbang DID and Ukong had evaporation measurements (Figure 1 and Table 2). Before proceeding to the calculation of the missing data and the statistical trend determination, a lag-1 autocorrelation test was performed to determine the method of preprocessing and to test for homogeneity and correlation (Tabari et al., 2012; Zhaoet al., 2015). This test classified data as homogenous but non-correlated and suggests the suitability of the collected data for statistical analysis without applying any pre-whitening technique to determine the trend. Then the missing data in the daily evaporation were filled through a normal ratio method. Three different statistical methods, the linear regression model (parametric), Mann Kendall and Sen's slope (non-parametric) methods were implemented to determine the trends in rainy and non-rainy days and evaporation in the LRB. Statistical techniques used in the present study are discussed in the following sections.

Linear regression model

Linear regression model is one of the simplest method to calculate the trend of time series data. Linear regression model considers normally distributed data and analyzes the linear relationship between time and variables, which are temporally and spatially independent (Helsel & Hirsch, 1992; Wang & Zhou, 2005; Soro et al., 2016). The linear regression equation is:

Y= a + bX (1)

where, Y is the dependent variable, X is the explanatory variable, b is the slope and a is the intercept. The slope of regression can describe the positive and negative nature of the time series data. An assumption of normal distribution is necessary for linear regression model. The null hypothesis is that the slope of the line is zero or no trend in the data. The probability value shows the significance of the slope. A positive value indicates an increasing linear trend and a negative value a decreasing linear trend (Brito et al., 2017).

Mann Kendall Test

The Mann Kendall test is the non-parametric test used to statistically assess the trend in a series of monotonic variables such as hydrological, climatic and environmental data which may be serially independent (Mann, 1945; Kendall, 1948). Mann Kendall test statistics are widely used to detect the trend characteristics of stream flow, rainfall and temperature in connection with the assessment for the impact of climate change on these variables (Partal & Kahya, 2006; Dindang et al, 2013; Diop et al.,2016).

Mann Kendall statistics for time series data of [X.sub.1], [X.sub.2], [X.sub.3], . . .[X.sub.n] can be determined as:

[mathematical expression not reproducible] (2)

where, [X.sub.i] and [X.sub.j] are ith and jth value of data, n is the number of observations and sgn (Xj - Xi) is the sign function calculated as :

[mathematical expression not reproducible] (3)

If the number of observation n> 10, the mean ([micro](S)) and the variance ([sigma] (2) (S)) are calculated using the following equation:

[micro] (s) = 0 (4)

[mathematical expression not reproducible] (5)

where, ti is the number of observations in the ith tied group and m is the number tied groups. If the number of ties are absent in the observation then the variance is calculated as:

[mathematical expression not reproducible] (6)

The normalised standard test statistics is calculated as:

[mathematical expression not reproducible] (7)

The [Z.sub.s] value will determine whether the observed data possess a statistically significant trend. If the [Z.sub.s] value is greater than zero this indicates a positive trend and a value, less than zero signifies a negative trend. To test the monotonic trend, either positive or negative, at [alpha] level of significance, the null hypothesis [H.sub.0] (absence of monotonic trend/ identical distribution) is rejected if the absolute value of [Z.sub.s] is greater than [Z.sub.s(1-[alpha]/2], and the alternative hypothesis [H.sub.1] (follows a monotonic trend) is accepted (Turgay & Ercan, 2006; Karpouzos et al., 2010; Tabari et al., 2015).

Sen's slope

The Sen's slope method helps to predict the magnitude of the trend of the variable under consideration by estimating the slope of the trend (Yue et al., 2002). The Sen's slope estimator is considered an alternative to the parametric least-squares regression line, which estimates the slope using a weighted mean, whereas the Sen's method uses median values to estimate the true slope of the existing trend (Sen, 1968; Gilbert, 1987). The slope T. can be calculated as:

[mathematical expression not reproducible] (8)

where, Xi and Xk are the data values at the time j and k (j>k) and with i varies 1, 2, 3, ...,N. The Sen's estimator of slope is the median of N values of T. and can be calculated as:

[mathematical expression not reproducible] (9)

A two-sided non-parametric test at 100 (1- [alpha]) % confidence interval is used to compute the value of [Q.sub.med]. A positive value if [Q.sub.i] represents an upward trend, a negative value demonstrates a downward trend and no trend is represented by zero.

The annual distribution of rainy days and non-rainy days were used to identify trends in rainy and non-rainy days whereas monthly evaporation rates were used to calculate evaporation trends. Results of the trend patterns of rainy, non-rainy days and evaporation are given in Tables 3 and 4 and are discussed in detail below.

Results and Discussion

General characteristics of rain and evaporation gauging stations

The rainfall analyzed at three rain gauge stations in the LRB showed a range of characteristics (Table 1). The monthly mean evaporation rate measured by the two stations also showed fluctuation during each month and varied widely between the stations (Table 2). The rain gauge station Limbang DID, located in the lower alluvial plains and peat land of the LRB showed an annual rainfall in the range of 2,033 mm to 4,554 mm with a mean of 3,093 mm. In the case of rainy and non-rainy days, the Limbang DID station recorded a mean total of 203 and 153 rainy and non-rainy days, respectively. The mean monthly evaporation rate measured varied between 125 mm (February) to 149 mm (August). Regarding the Ukong station, which is located in the center of the LRB, the annual rainfall ranged from 3,152 mm to 5,233 mm with a mean rainfall of 3,885 mm. This station recorded 183 mean rainy days and 165 mean non-rainy days and the mean monthly evaporation varied between 94 mm (February) to 115 mm (August) during the period of analysis. At the Long Napir station, the annual rainfall amount ranged from 1,833 mm to 3,582 mm with a mean rainfall of 2,480 mm and indicated 195 mean rainy and 123 non-rainy days. Evaporation measurements were not present at this station. Comparing the three rain gauge stations, the Ukong station exhibited the longest rainfall record (45 years) and recorded the highest annual rainfall and mean rainfall of all three stations. However, when analyzing the mean rainy days, the Limbang DID station recorded the maximum rainy days followed by the Long Napir station, whereas the Ukong station logged more non-rainy days. The mean monthly evaporation rate was comparatively higher at the Limbang DID station than at the Ukong station. These results indicate the ideal location of the Ukong station (located in the central catchment region of the Limbang River Basin) which received the maximum rainfall recorded during the least number of days and where less water evaporated despite the daily temperature range being more or less similar throughout the region.

Statistical trend characteristics of rainy and non- rainy days and evaporation Rainy day and non-rainy day trend

Trend of rainy days and non-rainy days recorded at gauging stations in the LRB indicate an increasing trend in rainy days and decreasing trend in non-rainy days with varying statistical significance.The results of the trend analysis are given in Table 3 and the linear relationship of rainy and non-rainy days during assessed periods are shown in Figure 2. Results show an increasing trend in rainy days at 99% confidence at all three stations. The linear regression model shows an increase of 2.08 rainy days/year at the Limbang DID station (Figure 2A) whereas the Mann Kendall and the Sen's slope estimates show the increase of rainy days as 4.16 and 2.27 days/year, respectively. Regarding the rain gauge station at Ukong, the data showed the rainy days increased by 2.03 days/year using linear regression (Figure 2B), 3.53 days/year by means of Mann Kendall and 2.75 via Sen's slope methods of trend estimation. Among the stations, Long Napir showed the maximum increase in rainy days i.e. 4.91 days/year employing the Mann Kendall test, 3.8 utilizing the Sen's slope estimation and 2.89 days/year applying linear regression (Figure 2C) methods.

The data for non-rainy days' showed a direct correlation with the rainy days data as seen in the linear trend graphs. In the present analysis, the representative rain gauge stations selected in the LRB showed significant and a non-significant decreasing trend in non-rainy days (Table 3 and Figure 2D-F). Non-rainy days recorded at the Limbang DID station showed a statistically significant (99% confidence) decreasing trend in all three methods (Figure 2D), used at -1.48 days/year in linear regression, Mann Kendall (-2.92 days/year) and Sen's slope (-1.7 days/year). The Ukong station showed the maximum decreasing trend of -3.49 days/year using the Mann Kendall followed by Sen's slope with -2.76 and -2.03 days/year in the linear regression test (Figure 2E). The Long Napir station showed a non-significant decreasing trend in non-rainy days except for the linear regression model which showed a decreasing trend with a value of -0.65 days/year (Figure 2F), whereas the Mann Kendall and Sen's slope methods showed -0.69 days/year and -0.15 days/year, respectively.

Overall, the results point towards increasing trend characteristics of rainy days. At the same time, regional spatial differences in the trend patterns were also found in the LRB. This is clearly demonstrated by the trend characteristics of the Long Napir and Ukong stations, which show a comparatively higher increasing trend in rainy days and a decreasing trend in non-rainy days compared to the Limbang DID station. These results may be influenced by the geographical location of the rain gauging stations. The Long Napir and Ukong stations are located in comparatively higher elevations and are surrounded by mountains of varying elevation, whereas the Limbang DID station is situated near the coastal area and will experience the direct influence of wind blowing from the South China Sea. The increasing trend characteristics of rainy days can be correlated with the amount of annual rainfall received in the region during the last few decades. Studies on rainfall trend characteristics carried out in Sarawak suggest an increasing trend in rainfall in the region, which supports the results of the present study (Tuh et al., 2013, Sa'adi et al.,2017a). Moreover, a recent study conducted by Krishnan et al., (2019 a,b) on the monthly and annual rainfall trend characteristics of the Limbang River Basin identified a statistically significant increase in monthly rainfall in the river basin.

Monthly evaporation trends

The results of monthly evaporation trend characteristics analyzed through linear regression, Mann Kendall and Sen's slope statistics are given in table 4 and the monthly linear trend characteristics are shown in Figures 3 and 4. During the month of January, the Limbang DID station showed a nonsignificant increasing trend in evaporation (0.55 mm/year, 1.02 mm/year and 0.75mm/year, respectively for linear regression model, Mann Kendall and Sen's slope methods). The Ukong station also showed an increasing trend in evaporation with a 95% confidence level in the Mann Kendall test (1.75 mm/year) compared to the other two methods, which showed a non-significant increase (0.53 mm/year and 1.02 mm/year, respectively for linear regression and Sen's slope methods). In the month of February the Limbang DID station showed a statistically significant increasing trend (95% confidence) in the linear regression (0.74 mm/year) Mann Kendall' (2.04 mm/year) and Sen's slope (0.98) tests. The Ukong station continued to show non-significant increasing trends in the linear regression and Sen's slope methods, whereas in the Mann Kendall test, it showed a statistically significant (90% confidence) increasing trend of 1.8mm/year. In the months of March and April, the Limbang DID and Ukong stations showed a non-significant increasing trend in evaporation.

The Limbang DID and Ukong stations continued to show a non-significant increasing trend in evaporation during the month of May while a maximum increasing trend was shown using the Mann Kendall test at 0.94 mm/ year and 0.85 mm/year for the Limbang DID and Ukong stations, respectively. In June, the Limbang DID station continued to show a non-significant increasing trend in evaporation with a maximum of 1.25 mm/year. A change in trend characteristics was observed at the Ukong station, which showed a non-significant decreasing trend in the monthly evaporation rate in the Mann Kendall and Sen's slope tests (-0.40 mm/year and -0.07 mm/year, respectively) whereas in linear regression tests it showed a slight increasing trend (0.21 mm/year). In the month of July, evaporation measured at both stations showed significant as well as non-significant increasing trends. The Limbang DID station showed an increase of 1.75 mm/year in the Mann Kendall test and 0.79 mm/year in the linear regression with 90% confidence. The Sen's slope method also showed an increasing trend (0.77). During the same month, the Ukong station showed a maximum increase (non-significant) of 1.37 mm/year by means of the Mann Kendall test, whereas the linear regression method showed a statistically significant (90% confidence) increase of 0.63 mm/year. The Sen's slope method also showed an increasing trend (0.34) during this month. In August, both stations showed a non-significant increasing trend. The Limbang DID station showed a maximum increase of 0.99 mm/year using a linear regression, whereas the Ukong station showed an increase of 0.45 mm/year (Mann Kendall test).

During September, both stations showed distinct characteristics of a trend in evaporation rates. The Limbang DID station showed a statistically significant trend whereas, the Ukong station showed a non-significant decreasing trend as well as no trend. The Limbang DID station data using the Mann Kendall and linear regression tests, revealed an increase in evaporation rates of 2.55 mm/year (95% confidence) and 0.79 mm/year (99% confidence), respectively. The Sen's slope trend method also showed an increasing nature (0.53). While the Ukong station under the linear regression method showed a non-significant decreasing trend (-0.10 mm/year), the Mann Kendall test showed no trend (0) and the Sen's slope showed a slightly increasing trend (0.20). In the month of October, both stations showed increasing trends. The Limbang DID station showed an increasing trend in the evaporation rate by means of the linear regression (1.00 mm/year at 99% confidence), the Mann Kendall test (2.55 mm/year at 95%) and a comparatively higher increase of 1.12 for the Sen's slope test. The Ukong station showed a maximum increase of 1.83 mm/year (90% confidence) using the Mann Kendall test and a non-significant increase of 0.28 mm/year via a linear regression and 0.82 mm/year via a Sen's slope test. During the months of November and December, both stations showed an increasing trend in monthly evaporation rates, but a statistically significant increase was noted during the month of December only. During November, the Limbang DID station showed a maximum increase of 0.75 mm/year in the Mann Kendall test, whereas the Ukong station showed a maximum of 0.60 mm/year. During the month of December, the Limbang DID station showed a statistically significant (95% confidence) increasing trend in linear regression (0.78 mm/year) and a non-significant increasing trend in Mann Kendall (1.54 mm/year) test. The Ukong station showed a statistically significant (99% confidence) increasing trend in the linear regression (0.65 mm/year) and Mann Kendall (2.61 mm/year) tests. The Sen's slope method also showed increasing trends for the Limbang DID and Ukong stations (1.01 and 0.86, respectively).

Overall, the monthly evaporation rate in selected rain gauge stations in the LRB showed an increasing trend with varying statistical significance. In the Limbang DID station, statistically significant increasing trends were associated with the months of February, July, September, October and December. At the Ukong station a significant increase in the evaporation rate was observed during the months of January, February, July, October and December and a non-significant decreasing trend during June and September. Being located in an equatorial tropical rainforest region with very dynamic climatic characteristics and varying terrain characteristics likely contributed to comparatively higher evaporation rates. Daily variation in temperature, wind speed, relative humidity and rainfall is the major factors controlling the rate of evaporation for this region. Further, it was also noted that all the regional scale climate models suggest the presence of extreme rainfall with increasing trend in temperature in these regions (Christensen et al., 2007; Chotamonsak et al., 2011). Comparing the results of evaporation trends (present research) with monthly rainfall trend pattern in the region carried out by Krishnan et al., (2019a), both shows very correlating characteristics. It was noted that monthly rainfall shows increasing trend characteristics in most months similar to the evaporation trends particularly for Limbang DID and Ukong station. Both the stations record maximum monthly rainfall higher than 750mm and average annual rainfall > 3000mm. This suggests the increased rate of evaporation associated with heavy rainfall in the region. If the rate of rainfall is lower, then the evaporation rate less and vice versa and also indicate the close relationship between local rainfall and evaporation rate (Yang et al., 2016).This might be due to the increase in difference of maximum and minimum temperatures over these regions and also due to sunshine duration, wind speed and saturation demand of the atmosphere (Mondal et al., 2013; Madhu et al., 2015; Houngue et al., 2019).

In order to make a comprehensive analysis of the relationship between mean evaporation rates with rainfall, temperature and relative humidity, linear regression and correlation analyses were performed (Figure 5). It was noted that, there was a negative correlation between the monthly mean rainfall and the monthly mean evaporation (r = -0.43 and [R.sup.2] = 0.18), whereas a positive correlation was observed between the monthly mean temperature and the mean evaporation (r = 0.46 and [R.sup.2] = 0.21). Though the monthly mean rainfall and evaporation shown a negative correlation, as mentioned earlier, the evaporation rate is highly depends on the local rainfall characteristics. However, relative humidity and evaporation showed a comparatively higher negative correlation (r = -0.57 and [R.sup.2] = 0.32). Though the correlation and linear regression was moderate, the combined effects of these parameters likely influenced the higher rate of evaporation in the LRB. In order to establish and understand the exact co-existence of hydro-meteorological variables more detailed measurements of rainfall, daily temperature, humidity and wind speed along with evaporation is necessary.

Conclusion

Statistical trend analysis for annual rainy days, non-rainy days and evaporation data from the Limbang River Basin analyzed through linear regression model, Mann Kendall and Sen's slope methods indicates significant increasing trends in annual rainy days and evaporation rates and a decreasing trend in non-rainy days. All the three stations showed an increasing trend in annual rainy days with highest increase in Long Napir (4.91 days/year). Similarly in the case of non-rainy days trend, the Ukong station showed the maximum decrease (-3.49 days/year). Increase in the number of rainy days in the study area suggests a positive impact on the rainfall amount in the region. Monthly evaporation trends showed an increasing trend throughout the year at the Limbang DID station whereas Ukong station showed varying trend characteristics with a maximum increase in December and a decreasing trend in June and September. Correlation analysis between the mean evaporation rate with the monthly mean rainfall, monthly mean temperature and relative humidity suggests a combined influence on increased evaporation rates. Further, the increased evaporation rate in the region is linked to the local precipitation characteristics, diurnal variation in temperature, wind characteristics and humidity. Moreover, the LRB suffers from a lack of detailed observation records of evaporation from different locations, which are critical for deriving more comprehensive evaporation characteristics of the region. Overall, the Limbang River Basin is now experiencing a greater number of rainy days with an increasing evaporation rate. The present research provided a new applied theory in the field of statistical trend determination of hydro-meteorological parameters, especially for those who are working on hydroelectric projects and the impacts of climate change at a local/regional level.

References

Abatzoglou, J. T. (2013). Development of gridded surface meteorological data for ecological applications and modelling. International Journal of Climatology, 33(1), 121-131.

Abtew, W., Obeysekera, J., & Iricanin, N., (2011). Pan evaporation and potential evapotranspiration trends in South Florida. Hydrological Process, 25(6), 958-969. DOI:10.1002/hyp.7887.

Akinsanola, A. A., Ogunjobi, K. O., Ajayi, V. O., Adefsan, E. A., Omotosho, J. A. & Sanogo, S. (2017). Comparison of five gridded precipitation products at climatological scales over West Africa. Meteorology and Atmospheric Physics, 129(6), 669-689.

Allan, R.P., Soden, B. J., John, V. O., Ingram, W., & Good, P. (2010). Current changes in tropical precipitation. Environmental Research Letters, 5(2), 025205.

Asfaw, A., Simane, B., Hassen, A., & Bantider, A., (2018). Variability and time series trend analysis of rainfall and temperature in northcentral Ethiopia: A case study in Woleka sub-basin. Weather and Climate Extremes, 19, 29-41.

Balint, M., Domisch, S., Engelhardt, C. H. M., Haase, P. , Lehrian, S., ... Nowak, C. (2011). Cryptic biodiversity loss linked to global climate change. Nature Climate Change, 1(6), 313.

Bandyopadhyay, A., Bhadra, A., Raghuwanshi, N. S., & Singh, R. (2009). Temporal trends in estimates of reference evapotranspiration over India. Journal of Hydrological Engineering, 14(5), 508-515. DOI:10.1061/(ASCE)HE.1943-5584.0000006.

Behzadi, F., Wasti, A., Rahat, S. H., Tracy, J. N., & Ray, P. A. (2020). Analysis of the climate change signal in Mexico City given disagreeing data sources and scattered projections. Journal of Hydrology: Regional Studies, 27, 100662.

Benestad, R. E. (2013). Association between trends in daily rainfall percentiles and the global mean temperature. Journal of Geophysical Research: atmospheres, 118(19), 10802-10810. DOI: 10.1002/jgrd.50814.

Brito, T. T., Oliveira-Junior, J. F., Lyra, G. B., Gois, G., & Zeri, M. (2017). Multivariate analysis applied to monthly rainfall over Rio de Janeiro state, Brazil. Meteorology and Atmospheric Physics, 129(5), 469-478.

Burn, D. H., & Hesch, N. M. (2007). Trends in evaporation for the Canadian Prairies. Journal of Hydrology, 336(1-2), 61-73. DOI:10.1016/j.jhydrol.2006.12.011.

Cecil, D. J., Buechler, D. E., & Blakeslee, R. J. (2014). Gridded lightning climatology from TRMM-LIS and OTD: Dataset description. Atmospheric Research, 135, 404-414.

CheRos, F., Tosaka, H., Sidek, L.M., & Basri, H., (2016). Homogeneity and trends in long-term rainfall data, Kelantan River Basin, Malaysia. International Journal of River Basin Management, 14(2), 151-163.

Chotamonsak, C., Salathe Jr., E. P. , Kreasuwan, J., Chantara, S., & Siriwitayakorn, K. (2011). Projected climate change over Southeast Asia simulated using a WRF regional climate model. Atmospheric Science Letters, 12(2), 213-219.

Christensen, J.H., Hewitson, B., Busuioc, A., Chen, A., Gao, X., Held, I. (...), Whetton, P. (2007). Regional Climate Projections. S. Solomon, D. Qin, M. Manning, Z. Chen, M. Marquis, K. B. Averyt, M. Tignor and H. L. Miller (Eds). In: Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA

Cretaux, J. F., Jelinskia, W., Calmantb, S., Kouraevc, A., Vuglinskid, V. , (...), Maisongrande, P. (2011). SOLS: A lake database to monitor in the Near Real Time water level and storage variations from remote sensing data. Advances in Space Research, 47(9), 1497-1507. DOI: 10.1016/j.asr.2011.01.004

Crisci, A., Gozzini, B., Meneguzzo, F., Pagliara, S., & Maracchi, G. (2002). Extreme rainfall in a changing climate: regional analysis and hydrological implications in Tuscany. Hydrological Processes, 16(6), 1261-1274. doi:10.1002/hyp.1061.

Deni, S. M., Suhaila, J., Zin, W. Z. W., & Jemain, A. A. (2010). Spatial trends of dry spells over Peninsular Malaysia during monsoon seasons. Theoretical and Applied Climatology, 99(3-4), 357.

DID. (2015). The Sarawak hydrological yearbook 2015. Vol 42. Department of Irrigation and Drainage, Sarawak, Malaysia.

DID. (2017). Hydrological station register (revised and updated 2017). Department of Irrigation and Drainage Malaysia. http://[H.sub.2]O.water.gov.my/man_hp1/HP7_2017.pdf (last accessed on 20.06.2018)

Dindang, A., Taat, A. B., Beng, P. E., Alwi, A. B. M., Mandai, A. A., (...) Lah, D. (2013). Statistical and trend analysis of rainfall data in Kuching, Sarawak from 1968 -- 2010. Research Publication No. 6/2013, Malaysian Meteorological Department, Kuching.

Diop, L., Ansoumana, B., & Dior., (2016). Spatio temporal Trend Analysis of the Mean Annual Rainfall in Senegal. European Science Journal, 12(12), 231-245. DOI: 10.19044/esj.2016.v12n12p231.

Djaman, K., Koudahe, K., & Ganyo, K. K. (2017). Trend analysis in annual and monthly pan evaporation and pan coefficient in the context of climate change in Togo. Journal of Geoscience and Environment Protection, 5(12), 41-56.

Dumitrescu, A., & Birsan, M. V. (2015). ROCADA: a gridded daily climatic dataset over Romania (1961 -- 2013) for nine meteorological variables. Natural Hazards, 78(2), 1045-1063.

Eymen, A., & Koylu, U. (2019). Seasonal trend analysis and ARIMA modeling of relative humidity and wind speed time series around Yamula Dam. Meteorology and Atmospheric Physics, 131(3), 601-612.

Frich, P. , Alexander, L., Della-Marta, P. M., Gleason, B., Haylock, M., (...), Peterson, T. C. (2002). Observed coherent changes in climatic extremes during the second half of the twentieth century. Climate Research, 19(3), 193-212. DOI:10.3354/cr019193.

Fu, G., Charles, S. P. , & Yu , J., (2009). A critical overview of pan evaporation trends over the last 50 years. Climatic Change, 97(1), 193 - 214. DOI:10.1007/s10584-009-9579-1

Gavrilov, M. B., Markovic, S., Janc, N., Nikolic, M., Valjarevic, A. D., (...), Bacevic, N. (2018). Assessing average annual air temperature trends using the Mann--Kendall test in Kosovo. Acta Geographica Slovenica, 58(1).

Gebremicael, T. G., Mohamed, Y. A., & Hagos, E. Y. , (2017). Temporal and spatial changes of rainfall and streamflow in the upper Tekeze--Atbara river basin, Ethiopia. Hydrology and Earth System Sciences, 21(4), 2127-2142.

Ghumman, A. R., Ghazaw, Y. M., Alodah, A., Rauf, A., Shafquzzaman, M., & Haider, H. (2020). Identifcation of Parameters of Evaporation Equations Using an Optimization Technique Based on Pan Evaporation. Water, 12(1), 228.

Gilbert, R. O., (1987). Statistical methods for environmental pollution monitoring. John Wiley & Sons.

Hassan, I., Kalin, R. M., White, C. J., & Aladejana, J. A. (2020). Evaluation of daily gridded meteorological datasets over the Niger Delta region of Nigeria and implication to water resources management. Atmospheric and Climate Sciences, 10(1), 21-39.

Helsel, D. R., & Hirsch, R. M., (1992). Statistical Methods in Water Resources. In: Studies in Environmental Science 49, Elsevier: New York, N Y, USA.

Houngue, R., Lawin, A. E., Moumouni, S., & Afouda, A. A., (2019). Change in Climate Extremes and Pan Evaporation Infuencing Factors over Oueme Delta in Benin. Climate, 7(1), 2.

Hua, T. M., Hui, R. C. Y. , & Husen, R., (2013). Trends of rainfall in Sarawak from 1999 to 2008. In: Proceedings of the International Conference on Social Science Research.

IPCC. (2007). Climate change 2007: Synthesis report. Valencia, Spain.

Jana, C., Sharma, G. C., Alam, N. M., Mishra, P. K., Dubey, S. K., & Kumar, R. (2016). Trend analysis of rainfall and rainy days of Acra. International Journal of Agriculture Statistics and Science, 12(1), 263-270.

Jhajharia, D., Yadav, B. K., Maske, S., Chattopadhyay, S., & Kar, A. K. (2012). Identifcation of trends in rainfall, rainy days and 24h maximum rainfall over subtropical Assam in Northeast India. Comptes Rendus Geoscience 344(1), 1-13. DOI:10.1016/j.crte.2011.11.002.

Karpouzos, D. K., Kavalieratou, S., & Babajimopoulos, C. (2010). Trend analysis of precipitation data in Pieria Region (Greece). European Water, 30, 31-40.

Kendall, M. G. (1948). Rank Correlation Methods. Griffn, London.

Krishnan, M. V. N., Prasanna, M. V. , & Vijith, H. (2019a). Statistical analysis of trends in monthly precipitation at the Limbang River Basin, Sarawak (NW Borneo), Malaysia. Meteorology and Atmospheric Physics, 131(4), 883-896.

Krishnan, M. V. N., Prasanna, M. V. , & Vijith, H., (2019b). Annual and Seasonal Rainfall Trends in an Equatorial Tropical River Basin in Malaysian Borneo. Environmental Modeling & Assessment, 24(5), 569-584.

Krishnan, M. V. N., Prasanna, M. V. , & Vijith, H. (2018). Fluctuations in Monthly and Annual Rainfall Trend in the Limbang River Basin, Malaysia: A Statistical Assessment to Detect the Influence of Climate Change. Journal of Climate Change, 4(2), 15-29.

Kumar, N., Yadav, P. B., Tyagi, A., & Jaswal, A. (2012). Trend and spatial distribution of rainfall & rainy days over Andaman & Nicobar Islands. Natural Hazards, 63, 575-587. DOI:10.1007/s11069-012-0173-x.

Kumar, V. , & Jain, S. K., (2010). Trends in seasonal and annual rainfall and rainy days in Kashmir Valley in the last century. Quaternary International, 212(1), 64-69. DOI:10.1016/j.quaint.2009.08.006.

Kumar, V. , & Jain, S. K. (2011). Trends in rainfall amount and number of rainy days in river basins of India (1951-2004). Hydrological Research, 42(4), 290-306. DOI:10.2166/nh.2011.067.

Lau, K. M., & Wu, H. T. (2007). Detecting trends in tropical rainfall characteristics, 1979 -- 2003. International Journal of Climatology, 27(8), 979-988.

Liu, Y. J., Chen, J., & Pan, T. (2019). Analysis of changes in reference evapotranspiration, pan evaporation, and actual evapotranspiration and their influencing factors in the North China Plain during 1998 -- 2005. Earth and Space Science, 6(8), 1366-1377.

Livada, I., Synnefa, A., Haddad, S., Paolini, R., Garshasbi, S., (...), Santamouris, M. (2019). Time series analysis of ambient air-temperature during the period 1970 -- 2016 over Sydney, Australia. Science of the Total Environment, 648, 1627-1638.

Madhu, S., Kumar, T. V. L., Barbosa, H., Rao, K. K., & Bhaskar, V. V. (2015). Trend analysis of evapotranspiration and its response to droughts over India. Theoretical and Applied Climatology, 121(1-2), 41-51.

Mandale, V. P., Jedhe, S. H., & Khadtare, M. Y., (2019). Spatio-temporal Trends of Rainfall and Rainy Days in the Marathwada Region of Maharashtra State. Climate Change, 5(17), 55-61.

Mann, H. B., (1945). Nonparametric tests against trend. Econometrica, 13, 245-259.

Manton, M. J., Della-Marta, P. M., Haylock, M. R., Hennessy, K. J., Nicholls, N., (...), Yee, D. (2001). Trends In extreme daily rainfall and temperature in Southeast Asia and the South Pacific: 1961-1998. International Journal of Climatology, 21, 269-284.

Mikkonen, S., Laine, M., Makela, H. M., Gregow, H., Tuomenvirta, H., (...), Laaksonen, A. (2015). Trends in the average temperature in Finland, 1847 -- 2013. Stochastic Environmental Research and Risk Assessment, 29(6), 1521-1529.

Mosmann, V. , Castro, A., Fraile, R., Dessens, J., & Sanchez, J. L. (2004). Detection of statistically significant trends in the summer precipitation of mainland Spain. Atmospheric Research, 70(1), 43-53. DOI:10.1016/j.atmosres.2003.11.002.

Nyatuame, M., Owusu-Gyimah, V. , & Ampiaw, F., (2014). Statistical analysis of rainfall trend for Volta Region in Ghana. International Journal of Atmospheric Science, 1- 11. doi:10.1155/2014/203245.

Oliveira, P. T., Silva, C. S., & Lima, K. C. (2017). Climatology and trend analysis of extreme precipitation in subregions of Northeast Brazil. Theoretical and Applied Climatology, 130(1-2), 77-90.

Palizdan, N., Falamarzi, Y. , Huang, Y. F., Lee, T. S., & Ghazali, A. H. (2014). Regional precipitation trend analysis at the Langat River Basin, Selangor, Malaysia. Theoretical and Applied Climatology, 11 7 (3-4), 589-606.

Palizdan, N., Falamarzi, Y. , Huang, Y. F., Lee T. S., & Ghazali, A. H. (2015). Temporal precipitation trend analysis at the Langat River Basin, Selangor, Malaysia. Journal of Earth System Sciences, 124(8), 1623-1638.

Palizdan, N., Falamarzi, Y. , Huang, Y. F., & Lee, T.S., (2017). Precipitation trend analysis using discrete wavelet transform at the Langat River Basin, Selangor, Malaysia. Stochastic Environmental Research and Risk Assessment, 31(4), 853-877.

Paparrizos, S., & Matzarakis, A., (2017). Assessment of future climate change impacts on the hydrological regime of selected Greek areas with different climate conditions. Hydrology Research, 48(5), 1327-1342

Parry, L., Harrison, S., Betts, R., Shannon, S., Jones, D. B., & Knight, J. (2020). Impacts of Climate Change on Himalayan Glaciers: Processes, Predictions and Uncertainties. In: A. P. Dimri, B. Bookhagen, M. Stoffel, & T. Yasunari (Eds). Himalayan Weather and Climate and their Impact on the Environment (pp. 331-349). Springer, Cham.

Partal, T., & Kahya, E. (2006). Trend analysis in Turkish precipitation data. Hydrological Processes, 20(9), 2011-2026. DOI: 10.1002/hyp.5993.

Piao, S., Ciais, P. , Huang, Y. , Shen, Z., Peng, S., (...), Fang, J. (2010). The impacts of climate change on water resources and agriculture in China. Nature, 467(7311), 43.

Pour, S. H., Wahab, A. K. A., & Shahid, S. (2020). Physical-empirical models for prediction of seasonal rainfall extremes of Peninsular Malaysia. Atmospheric Research, 233, 104720.

Raziei, T., Arasteh, P. D., & Saghafan, B. (2005). Annual rainfall trend in arid & semi-arid regions of Iran. ICID 21st ( )European regional Conference, 15-19 May 2005--Frankfurt (Oder) and Slubice--Germany and Poland, 20-28.

Rhee, J., & Im, J., (2017). Meteorological drought forecasting for ungauged areas based on machine learning: Using long-range climate forecast and remote sensing data. Agriculture and Forest Meteorology, 237-238, 105-122. DOI: 10.1016/j.agrformet.2017.02.011.

Rodell, M, Houser, P. R., Jambor, U., Gottschalck, J., Mitchell, K. (2004). The global land data assimilation system. Bulletin of American Meteorological Society, 85(3), 381-394. DOI: 10.1 I75/BAMS-85-3-38I.

Ruiz-Alvarez, O., Singh, V. P. , Enciso-Medina, J., Munster, C., Kaiser, R., (...), Dos Santos, C. A. C. (2019). Spatio-temporal trends in monthly pan evaporation in Aguascalientes, Mexico. Theoretical and Applied Climatology, 136(1-2), 775-789.

Sa'adi, Z., Shahid, S., Ismail, T., Chung, E. S., & Wang, X. J. (2019). Trends analysis of rainfall and rainfall extremes in Sarawak, Malaysia using modified Mann-Kendall test. Meteorology and Atmospheric Physics, 131(3), 263-277.

Sa'adi, Z., Shahid, S., Ismail, T., Chung, E. S., & Wang, X. J. (2017a). Distributional changes in rainfall and river flow in Sarawak, Malaysia. Asia-Pacific Journal of Atmospheric Sciences, 53(4), 489-500.

Sa'adi, Z., Shahid, S., Chung, E. S., & Bin Ismail, T. (2017b). Projection of spatial and temporal changes of rainfall in Sarawak of Borneo Island using statistical downscaling of CMIP5 models. Atmospheric research, 197, 446-460.

Sayyad, R. S., Dakhore, K. K., & Phad, S.V. (2019). Analysis of rainfall trend of Parbhani, Maharshtra using Mann-Kendall test. Econometrica, 13, 245-259.

Schloss, A. L., Vorosmarty, C., Willmott, C. J., & Choudhury, B. J. (1996). Analyzing the discharge regime of a large tropical river through remote sensing, ground-based climatic data, and modelling. Water Resources Research, 32(10), 3137-3150. DOI: 10.1029/96WR01333.

Sen, P. K., (1968). Estimates of the regression coefficient based on Kendall's Tau. American Statistical Journal, 63 (324), 1379-1389. DOI:10.2307/2285891.

Serrano-Notivoli, R., Martin-Vide, J., Saz, M. A., Longares, L. A., Begueria, S., (...), Luis, M. (2018). Spatio-temporal variability of daily precipitation concentration in Spain based on a high-resolution gridded data set. International Journal of Climatology, 38, e518-e530.

Shadmani, M., Marof, S., & Roknian, M. (2012). Trend analysis in reference evapotranspiration using Mann-Kendall and Spearman's Rho tests in arid regions of Iran. Water Resources Management, 26(1), 211-224. DOI: 10.1007/s11269-011-9913-z.

Sharad, K. J., & Vijay, K., (2012). Trend analysis of rainfall and temperature data for India. Current Science, 102(1), 37 -49.

Simpson, I. R., & McCarthy, M. P. (2019). Structural and sampling uncertainty in observed UK daily precipitation extremes derived from an intercomparison of gridded data sets. International Journal of Climatology, 39(1), 128-142.

Singh, V., & Xiaosheng, Q. (2019). Data assimilation for constructing long-term gridded daily rainfall time series over Southeast Asia. Climate Dynamics, 53(5-6), 3289-3313.

Sobral, B. S., Oliveira-Junior, J. F., Alecrim, F., Gois, G., Muniz-Junior, J. G., (...), Zeri, M. (2020). PERSIANN-CDR based characterization and trend analysis of annual rainfall in Rio De Janeiro State, Brazil. Atmospheric Research, 238, 104873.

Sobral, B. S., Oliveira-Junior, J. F., Gois, G., Pereira-Junior, E. R., Terassi, P. M. B. (2019). Drought characterization for the state of Rio de Janeiro based on the annual SPI index: trends, statistical tests and its relation with ENSO. Atmospheric research, 220, 141-154.

Soltani, S., Saboohi, R., & Yaghmaei, L., (2012). Rainfall and rainy days trend in Iran. Climatic Change, 110(1-2), 187-213. DOI: 10.1007/s10584-011-0146-1.

Soro, G. E., Noufe, D., Goula Bi, T. A., & Shorohou, B. (2016). Trend analysis for extreme rainfall at sub-daily and daily timescales in Cote d'Ivoire. Climate, 4(3), 37.

Suhaila, J., & Yusop, Z. (2018). Trend analysis and change point detection of annual and seasonal temperature series in Peninsular Malaysia. Meteorology and Atmospheric Physics, 130(5), 565-581.

Suhaila, J., & Jemain, A.A. (2012). Spatial analysis of daily rainfall intensity and concentration index in Peninsular Malaysia. Theoretical and Applied Climatology, 108(1-2), 235-245.

Suhaila, J., Deni, S. M., & Jemain, A. A. (2008). Detecting inhomogeneity of rainfall series in Peninsular Malaysia. Asia-Pacific Journal of Atmospheric Science, 44(4), 369-380.

Suhaila, J., Zin, W. Z. W., Deni, S. M., & Jemain, A. A. (2010). Trends in peninsular Malaysia rainfall data during the southwest monsoon and northeast monsoon seasons: 1975 -- 2004. Sains Malaysiana, 39(4), 533-542.

Tabari, H., Abghari, H., & Hosseinzadeh Talaee, P. (2012). Temporal trends and spatial characteristics of drought and rainfall in arid and semiarid regions of Iran. Hydrological Processes, 26(22), 3351-3361.

Tabari, H., Taye, M. T., & Willems, P. (2015). Statistical assessment of precipitation trends in the upper Blue Nile River basin. Stochastic Environmental Research and Risk Assessment, 29(7), 1751-1761.

Tahir, T., Hashim, A. M., & Yusof, K. W., (2018). Statistical downscaling of rainfall under transitional climate in Limbang River Basin by using SDSM. In: IOP Conference Series: Earth and Environmental Science (Vol. 140, No. 1, p. 012037). IOP Publishing.

Teodoro, P.E., Oliveira-Junior, J. F., Cunha, E. R., Correa, C., Torres, F. E., (...), Ribeiro, L. P. (2016). Cluster analysis applied to the spatial and temporal variability of monthly rainfall in Mato Grosso do Sul State, Brazil. Meteorology and Atmospheric Physics, 128(2), 197-209.

Thomas, C. D., Cameron, A., Green, R., Bakkenes, M., Beaumont, L. J., (...), Williams, S. (2004). Extinction risk from climate change. Nature, 427(6970), 145.

Thompson, I., Mackey, B., McNulty, S., Mosseler, A. (2009). Forest resilience, biodiversity, and climate change. In: Secretariat of the Convention on Biological Diversity, Montreal. Technical Series no. 43. 1-67. (Vol. 43, pp. 1-67).

Toride, K., Cawthorne, D. L., Ishida, K., Kavvas, M. L., Anderson, M. L. (2018). Long-term trend analysis on total and extreme precipitation over Shasta Dam watershed. Science of the Total Environment, 626, 244-254.

Tuh, M. H., Robin, C. Y. H., & Rafdah, H. (2013). Trends of rainfall in Sarawak from 1999 to 2008. Proceeding of the International Conference on Social Science Research, ICSSR 2013, Penang, Malaysia. Organized by WorldConferences.net.

Uduak, C., UdoInyang., & Edem, I. (2012). Analysis of rainfall trends in AkwaIbom state, Nigeria. Journal of Environmental Earth Science, 2(8), 60-70.

Vijith, H., Hurmain, A., & Dodge-Wan, D. (2018) Impacts of land use changes and land cover alteration on soil erosion rates and vulnerability of tropical mountain ranges in Borneo. Remote Sensing Applications: Society and Environment, 12, 57-69.

Walther, G. R., Post, E., Convey, P. , Menzel, A., Parmesan, A., (...), Bairlein, F. (2002). Ecological responses to recent climate change. Nature, 416(6879), 389-395. DOI: 10.1038/416389a.

Wang, Y., Xu, Y. , Tabari, H., Wang, J., Wang, Q., Song, S., & Hue, Z. (2020). Innovative trend analysis of annual and seasonal rainfall in the Yangtze River Delta, eastern China. Atmospheric Research, 231, 104673.

Wang, Y. Q. & Zhou, L. (2005). Observed trends in extreme precipitation events in China during 1961-2001 and the associated changes in large-scale circulation. Geophysical Research Letters, 32, L09707.

Xu, C.Y., Gong, L., Jiang, T., Chen, D., & Singh, V. P. (2006). Analysis of spatial distribution and temporal trend of reference evapotranspiration and pan evaporation in Changjiang (Yangtze River) catchment. Journal of Hydrology, 327(1-2), 81-93. DOI: 10.1016/j.jhydrol.2005.11.029.

Xu, X., Yang, X., Zhu, B., Tang, Z., Wu, H., & Xie, L. (2020). Characteristics of MERRA-2 black carbon variation in east China during 2000-2016. Atmospheric Environment, 222, 117140.

Xu, Z.X., Takeuchi, K., Ishidaira, H., & Li, J. Y. (2005). Long-term trend analysis for precipitation in Asia Pacific Friend river basin. Hydrological Processes, 19(18), 3517-3532. DOI: 10.1002/hyp.5846.

Xu, Z. X., Takeuchi, K., & Ishidaira, H. (2003). Monotonic trend and step changes in Japanese precipitation. Journal of Hydrology, 279(1-4), 144-150. DOI: 10.1016/S0022-1694(03)00178-1.

Yang, T., Sun, F., Liu, W., Wang, H., Wang, T., & Liu, C. (2019). Using Geo-detector to attribute spatio-temporal variation of pan evaporation across China in 1961-2001. International Journal of Climatology, 39(5), 2833-2840.

Yang, Z., Zhang, Q., & Hao, X. (2016). Evapotranspiration trend and its relationship with precipitation over the loess plateau during the last three decades. Advances in Meteorology, 2016.

Yatim, A. N. M., Latif, M. T., Ahamad, F., Khan, M. F., Nadzir, M. S. M., & Juneng, L. (2019). Observed Trends in Extreme Temperature over the Klang Valley, Malaysia. Advances in Atmospheric Sciences, 36(12), 1355-1370.

Yu, P. S., Yang, T. C., & Kuo, C. C. (2006). Evaluating long-term trends in annual and seasonal precipitation in Taiwan. Water Resources Management, 20(1), 1007-1023. DOI: 10.1007/s11269-006-9020-8.

Yue, S., Pilon, P., Phinney, B., & Cavadias, G. (2002). The influence of autocorrelation on the ability to detect trend in hydrological series. Hydrological Processes, 16(9), 1807-1829. DOI: 10.1002/hyp.1095.

Yue, S., & Wang, C. (2004). The Mann-Kendall test modified by effective sample size to detect trend in serially correlated hydrological series. Water Resources Management, 18(3), 201-218. DOI: 10.1023/B:WA RM.0000043140.61082.60.

Zafor, M. A., Farzana, S. Z., Chakraborty, A., & Rahman, A. (2016). Analysis of rainfall trends and variability at Sylhet region in Bangladesh. ARPN Journal of Engineering and Applied Sciences, 11(11), 6836 - 6846.

Zafor, M. A., Chakraborty, A., Muniruzzaman, S. M., & Mojumdar, S. R. (2016). Rainfall forecasting in northeastern part of Bangladesh using time series ARIMA Model. Research Journal of Engineering and Science, E-ISSN, 2278, 9472.

Zhang, X.L., Wang, S., Zhang, J. M., Wang, G., & Tang, X. Y. (2015). Temporal and spatial variability in precipitation trends in the Southeast Tibetan Plateau during 1961-2012. Climate of the Past Discussion, 11(1), 447-487. DOI: 10.5194/cpd-11-447-2015.

Zhao, H., Yang, S., Wang, Z., Zhou, X., Luo, Y., & Wu, L. (2015). Evaluating the suitability of TRMM satellite rainfall data for hydrological simulation using a distributed hydrological model in the Weihe River catchment in China. Journal of Geographical Science, 25(2), 177-195.

Ninu Krishnan MV (*), Prasanna M V, Vijith H

Department of Applied Geology, Faculty of Engineering and Science, Curtin University Malaysia, Malaysia

(*) Corresponding author: mvninukrishna@gmail.com

Record

Manuscript received: 05/02/2019

Accepted for publication: 28/04/2020

How to cite item

Ninu Krishnan, M. V. , Prasanna, M. V. , & Vijith, H. (2020). Trend characteristics of rainy days and evaporation at a tropical rainforest region in East Malaysia, Borneo. Earth Sciences Research Journal, 24(3), 305-315. DOI: https://doi.org/10.15446/esrj.v24n3.77670

https://doi.org/10.15446/esrj.v24n3.77670

Table 1. Details of rainfall, rainy days and non-rainy days at the rain gauging stations considered in the present study. Stations Data availability Elevation (m) Data missing (%) From To Limbang DID 1980 2015 7 4.67 Ukong 1970 2015 22 17.83 Long Napir 1977 2015 147 13.51 Stations Accumulated rainfall (mm) Rainy days Non rainy days Minimum Maximum Mean (mean) (mean) Limbang DID 2033.00 4554.20 3093.63 203 153 Ukong 3152.00 5233.00 3885.54 183 165 Long Napir 1855.00 3582.00 2480.71 195 123 Table 2. Mean monthly evaporation measured at evaporation gauging stations in LRB. Stations Data availability Mean monthly evaporation (mm) From To Jan Feb Mar Apr May Limbang DID 1983 2015 129.29 125.65 148.98 139.54 138.62 Ukong 1973 2004 106.40 94.14 119.47 111.80 109.78 Stations Mean monthly evaporation (mm) Jun Jul Aug Sep Oct Nov Dec Limbang DID 138.72 144.67 149.23 139.20 143.15 133.25 134.59 Ukong 111.22 113.02 115.69 110.99 108.90 101.92 104.92 Table 3. Statistical trend characteristics of rainy days and non-rainy days in LRB. Variables Analysis Methods Rain gauging stations Limbang Ukong Long DID Napir Linear regression 2.08 (***) 2.03 (***) 2.89 (***) model Rainy days MannKendall 4.16 (***) 3.53 (***) 4.91 (***) Sen's slope 2.27 2.75 3.8 Linear regression -1.48 (***) -2.03 (***) -0.65 model Non-rainy MannKendall -2.92 (***) -3.49 (***) -0.69 days Sen's slope -1.7 -2.76 -0.15 Note: bold values indicate statistically significant trends; (***) 99% confdence level. Table 4 statistical trend characteristics of monthly evaporation in LRB. Station Name Analysis Methods Months Jan Feb Mar Apr Linear regression model 0.55 0.74 (**) 0.06 0.12 Limbang DID MannKendall 1.02 2.04 (**) 0.12 0.21 Sen's slope 0.75 0.98 0.39 0.2 Linear regression model 0.53 0.58 -0.05 0.30 Ukong MannKendall 1.75 (*) 1.8 (*) 0.45 1.07 Sen's slope 1.02 0.9 0.22 0.58 Station Name Analysis Methods Months May Jun Jul Aug Linear regression model 0.49 0.40 0.79 (*) 0.99 Limbang DID MannKendall 0.94 1.25 1.75 (*) 0.58 Sen's slope 0.81 0.34 0.77 0.61 Linear regression model 0.23 0.21 0.63 (*) 0.17 Ukong MannKendall 0.85 -0.40 1.37 0.45 Sen's slope 0.33 -0.07 0.34 0.30 Station Name Analysis Methods Months Sep Oct Nov Linear regression model 0.79 (***) 1.00 (***) 0.49 Limbang DID MannKendall 2.55 (**) 2.47 (**) 0.75 Sen's slope 0.53 1.12 0.33 Linear regression model -0.10 0.28 0.05 Ukong MannKendall 0 1.83 (*) 0.6 Sen's slope 0.20 0.82 0.49 Station Name Analysis Methods Months Dec Linear regression model 0.789 (**) Limbang DID MannKendall 1.54 Sen's slope 1.01 Linear regression model 0.65 (***) Ukong MannKendall 2.61 (***) Sen's slope 0.86 Note: bold values: statistically significant trends; italics--decreasingtrends; (***) Confdence level 99%; (**) Confdence level 95%; (*) Confdence level 90%.

Printer friendly Cite/link Email Feedback | |

Author: | M.V., Ninu Krishnan; V., Prasanna M.; H., Vijith |
---|---|

Publication: | Earth Sciences Research Journal |

Geographic Code: | 9INDO |

Date: | Sep 1, 2020 |

Words: | 9448 |

Previous Article: | Study on the Evolution Monitoring of Topographic and Hydrological Characteristics of Small Watershed Based on Remote Sensing and GIS. |

Next Article: | Establishment of Evaluation Model for Shallow Geothermal Energy Resource Development Potential Based on Characteristic of Geotemperature. |

Topics: |