Geospatial Modeling for Investigating Spatial Pattern and Change Trend of Temperature and Rainfall

Bangladesh has been experiencing increased temperature and change in precipitation regime, which might adversely affect the important ecosystems in the country differentially. The river flows and groundwater recharge over space and time are determined by changes in temperature, evaporation and crucially precipitation. These again have a spatio-temporal dimension. This geospatial modeling research aimed at investigating spatial patterns and changing trends of temperature and rainfall within the geographical boundary of Bangladesh. This would facilitate better understanding the change pattern and their probable impacts on the ecosystem. The southeastern region, which is one of the most important forest ecosystem zones in the country, is experiencing early onset and withdrawal of rain but increasing trends in total rainfall except in the Monsoon season. This means that the region is experiencing a lower number of rainy days. However, total rainfall has not changed significantly. The differential between maximum and minimum showed an increasing trend. This changing pattern in average max and min temperature along with precipitation might cause a situation in which the species that are growing now may shift to suitable habitats elsewhere in the future. Consequently, the biodiversity, watersheds and fisheries, productivity of land, agriculture and food security in the region will be affected by these observed changes in climate.


Introduction
Bangladesh is the lowest riparian area in the Ganges, Brahmaputra and Meghna (GBM) River Basins.Nearly 93% of the basin area lies outside Bangladesh in India, Nepal, Bhutan and China [1].The climate in Bangladesh is mostly governed by its geographical location in the mouth of the Bay of Bengal and the deltaic flood plain in the foot hill of the Himalayas [2].Ocean circulations in the Indian Ocean, Bay of Bengal and the Tropical Channel have great influence on Bangladesh climate.Bangladesh has a humid, warm, tropical climate.Its climate is influenced primarily by the Monsoon and partly by Pre-Monsoon and Post-Monsoon circulations.The southwest Monsoon originates over the Indian Ocean and carries warm, moist, and unstable air [3].The Monsoon has its onset during the first week of June and ends in the first week of October, with some inter-annual variability in dates.Besides the Monsoon, the easterly trade winds are also active, providing warm and relatively drier circulation.In Bangladesh, there are four prominent seasons, namely, Winter (December to February), Pre-Monsoon (March to May), Monsoon (June to August), Post-Monsoon (September to November).
The Germanwatch Long-Term Climate Risk Index (CRI) 2011 has identified Bangladesh as the country most affected by extreme climate events during 1990-2009 [4].Many projected climate change (CC) impacts including sea level rise, higher temperatures, erratic Monsoon precipitation and run-off, potentially reduced dry season precipitation, and increase in cyclone frequency and intensity would in fact reinforce many other baseline stresses that already pose a serious impediment to the economic development of Bangladesh.In the past, a number of studies were carried out on trends in climate parameters over Bangladesh and South Asia [5].Chowdhury and Debsharma [6] and Mia [7] pointed out that temperature has been changed (by analyzing historical data of some selected meteorological stations).References [8,9] reported that mean annual temperatures of Bangladesh have increased during the period of 1895-1980 at 0.31 ˝C over the past two decades.Karmakar and Shrestha [10] using the 1961-1990 data for Bangladesh have projected that annual mean maximum temperature would increase to 0.4 ˝C and 0.73 ˝C by the year of 2050 and 2100, respectively.In terms of precipitation, most models show increased annual precipitation for Bangladesh.Some models also showed the possibility of increasing Monsoon precipitation (May to September), and decreasing dry season precipitation (December to February).The mean annual rainfall is about 2320 mm, but there are places with a mean annual rainfall of 6000 mm or more [11].A long duration of heavy rainfall associated with "norwester" thunderstorms is very common in Bangladesh [11,12].More than 80% of the 2,300 mm of annual precipitation that falls on Bangladesh comes during the Monsoon period [13].Annual rainfall ranges from 1200 mm in the extreme west to over 5000 mm in the east and northeast [14].Changes in the distribution of river flows and groundwater recharge over space and time are determined by changes in temperature, evaporation and, crucially, precipitation [15].Some CC impacts on hydrological processes have been observed already [16,17], and further changes are projected.This geospatial research aimed at investigating spatial patterns and changing trends of temperature and rainfall within the geographical boundary of Bangladesh and the interlinkage between them.This would facilitate better understanding of changing patterns, and their probable impacts on the ecosystem that were discussed in another author's paper [18] and a book chapter [19].
In short, geospatial is a methodology or set of analytical procedures used to derive information about spatial relationships among geographic phenomena.In most cases, geospatial models are two-dimensional, which take parameters on time and space to simulate a process and visualize a phenomena at given time and space.The evaluation of high resolution Regional Climate Models (RCMs) is seriously constrained by the lack of representative high-resolution observational datasets [18].Gaps in well distributed and quality datasets on climate parameters are also big challenges in carrying out dynamical three-dimensional experiments covering Bangladesh.It is also important to keep in mind that, while horizontal resolution is readily increased down to about 10 km, real vertical resolution in the surface layers is more a function of the boundary layer scheme than of numbers, which makes the computational time and resources prohibitively high.This study utilized geospatial techniques and tools for generating national to landscape level high-resolution geospatial scenarios of climate parameters.

Materials
Maximum and minimum monthly temperature and rainfall observation data from 34 meteorological stations of BMD (Bangladesh Meteorological Department) located all over Bangladesh were collected from the BMD database for a thirty-year period of 1978-2007 for this study.For data analysis spreadsheet packages of MS Excel and SPSS (Statistical Package for Social Science) were used, whereas ArcGIS 10.1 was used for geo-spatial analysis.A Spatial Analysis Tool in the ArcGIS 10.1 package was used for investigating spatial patterns in data.

Study Area
Bangladesh is located between 20 ˝to 26 ˝N and 88 ˝to 92 ˝E.It is bordered on the west, north and east by India, on the southeast by Myanmar, and on the south by the Bay of Bengal.The northern most area of the country is a few kilometers south of the Himalayans, which has profound influence on the climate in Bangladesh.Most of the country is low-lying land comprising mainly the delta of the Ganges and Brahmaputra rivers.Floodplains occupy 80% of the country.Mean elevations range from less than 1 meter on tidal floodplains, 1 to 3 meters on the main river and estuarine floodplains, and up to 6 meters in the Sylhet basin in the northeast [19].Only in the extreme northwest are elevations greater than 30 meters above the mean sea level.The northeast and southeast parts of the country are undulating and hilly (compared to the mostly deltaic floodplain), with some tertiary hills over 1000 meters above mean sea level [20].These two regions are historically home to the country's terrestrial tropical forest ecosystems.on the climate in Bangladesh.Most of the country is low-lying land comprising mainly the delta of the Ganges and Brahmaputra rivers.Floodplains occupy 80% of the country.Mean elevations range from less than 1 meter on tidal floodplains, 1 to 3 meters on the main river and estuarine floodplains, and up to 6 meters in the Sylhet basin in the northeast [19].Only in the extreme northwest are elevations greater than 30 meters above the mean sea level.The northeast and southeast parts of the country are undulating and hilly (compared to the mostly deltaic floodplain), with some tertiary hills over 1000 meters above mean sea level [20].These two regions are historically home to the country's terrestrial tropical forest ecosystems.

Methodology
Figure 1 presents conceptual and sequential analytical process flows of research.It is important to note that this paper is part of a larger piece of research and only the geospatial modeling part of the research is presented here.

Regression Model
Regression analysis is widely used for prediction and forecasting, where its use has substantial overlap with the field of machine learning.Regression analysis is also used to understand which among the independent variables are related to the dependent variable, and to explore the forms of these relationships.In restricted circumstances, regression analysis can be used to infer causal relationships between the independent and dependent variables.However this can lead to illusions or false relationships, so caution is advisable [21].
In statistics, regression analysis includes many techniques for and analyzing several variables in discerning the relationship between a dependent variable and one or more independent variables.It helps one understand how the value of the dependent variable changes when any one of the independent variables is varied, while the other independent variables are held fixed.
A large body of techniques for carrying out regression analysis has been developed.Familiar methods such as linear regression and ordinary least squares regression are parametric, in that the regression function is defined in terms of a finite number of unknown parameters that are estimated

Regression Model
Regression analysis is widely used for prediction and forecasting, where its use has substantial overlap with the field of machine learning.Regression analysis is also used to understand which among the independent variables are related to the dependent variable, and to explore the forms of these relationships.In restricted circumstances, regression analysis can be used to infer causal relationships between the independent and dependent variables.However this can lead to illusions or false relationships, so caution is advisable [21].
In statistics, regression analysis includes many techniques for and analyzing several variables in discerning the relationship between a dependent variable and one or more independent variables.It helps one understand how the value of the dependent variable changes when any one of the independent variables is varied, while the other independent variables are held fixed.
A large body of techniques for carrying out regression analysis has been developed.Familiar methods such as linear regression and ordinary least squares regression are parametric, in that the regression function is defined in terms of a finite number of unknown parameters that are estimated from the data.Nonparametric regression refers to techniques that allow the regression function to lie in a specified set of functions, which may be infinite-dimensional.
The performance of regression analysis methods in practice depends on the form of the data generating process, and how it relates to the regression approach being used.Since the true form of the data-generating process is generally not known, regression analysis often depends to some extent on making assumptions about this process.These assumptions are sometimes testable if a large amount of data is available.Regression models for prediction are often useful even when the assumptions are moderately violated, although they may not perform optimally.However, in many applications, especially with small effects or questions of causality based on observational data, regression methods give misleading results [22,23].

Mann Kendall Test
Mann-Kendall tests are non-parametric tests for detection of trend in a time series.These tests are widely used in environmental science because they are simple, robust and can cope with missing values and values below a detection limit.Since the first proposals of the test by [24,25], covariance between Mann-Kendall statistics were proposed by [26] and the test was extended in order to include seasonality [27], multiple monitoring sites [28] and covariates representing natural fluctuations [29].The present program (MS Excel and it's Add in XLSTAT) has incorporated all of these features.
By the Mann-Kendall test, this study tested the null hypothesis (H0) of no trend, i.e., the observations x i are randomly ordered in time, against the alternative hypothesis (H1), where there is an increasing or decreasing monotonic trend.The data values were evaluated as an ordered time series.Each data value was compared with all subsequent data values.If a data value from a later time period was higher than a data value from an earlier time period, the statistic S is incremented by 1.On the other hand, if the data value from a later time period was found to be lower than a data value sampled earlier, S had been decremented by 1.The net result of all such increments and decrements yields the final value of S [28].The M-K test statistic S is calculated using the formula: where x j and x k are the annual values in years j and k, j > k, respectively.If n < 10, the value of |S| was compared directly to the theoretical distribution of S derived by Mann and Kendall [30].The two-tailed test was used.At a certain probability level, H0 is rejected in favor of H1 if the absolute value of S equals or exceeds a specified value Sα/2, where Sα/2 is the smallest S, which has a probability less than α/2 to appear in case of no trend.A positive (negative) value of S indicates an upward (downward) trend [31,32].For n ě 10, the statistic S is approximately normally distributed with the mean and variance.

Mathematical Interpolation Methods
The geostatistical interpolation technique kriging was used because of its robustness in meteorological data applications (trend surface analysis, inverse distance weighting) or for their perceived advantages for climatic and phonological interpolations at a variety of spatial and temporal scales (ordinary kriging, partial thin plate splines).This method produces interpolation values by relying on an unknown mean value, allowing local influences due to nearby neighboring values (Figure 2).When the mean is unknown, there are few assumptions.This makes ordinary kriging particularly flexible.Temperature/Rainfall trends for all the weather stations were analyzed based on 30 years (1978-2007) of observation.To make a continuous scenario for the whole country, data interpolation was carried out in two steps.Step 1 Tessellated the whole study area into equal hexagons of 1200 sq.km.Imported a value (Temperature/Rainfall Trend) for all the hexagons from the nearest (Arial Distance) weather station using Spatial Joint Tool in ArcGIS.
Step 2 In the second step, the Kriging Geostatistical Technique was used to interpolate (Figure 2) the values using the Spatial Analyst tool in ArcGIS.Step 1 Tessellated the whole study area into equal hexagons of 1200 sq.km.Imported a value (Temperature/Rainfall Trend) for all the hexagons from the nearest (Arial Distance) weather station using Spatial Joint Tool in ArcGIS.

Bay of Bengal
Step 2 In the second step, the Kriging Geostatistical Technique was used to interpolate (Figure 2) the values using the Spatial Analyst tool in ArcGIS.
Among the 34 stations, four stations with unique climatic conditions were selected as representatives of the whole country for detailed analysis purposes.For regional study of temperature variability, Rajshahi, Barisal, Chittagong and Sylhet meteorological stations were selected as representatives of their own regional climatic conditions.Rajshahi is in the drought prone northwest region of Bangladesh.Barisal, the coastal flood plain is in the southwest while Chittagong is in the southeast hilly coastal region and Sylhet is in the northeast region (foot hills of Meghalaya hill range) region of Bangladesh.The whole year was sub-grouped into four seasons (i) Pre-Monsoon (March, April and May (MAM)), (ii) Monsoon (June, July and August (JJA)), (iii) Post-Monsoon (September, October and November (SON)) and (iv) Winter (December, January and February (DJF)).Trend analysis has been conducted on seasonal average data for all the stations.A basic linear regression model was adopted to analyze time trend and interrelation between temperature and rainfall using SPSS.The findings on trends in temperature and rainfall from 34 stations were then spread within the whole country to make a continuous scenario interactively using Geostatistical Kriging Technique and Spatial Analyst tools in the ArcGIS platform.
Patterns of meteorological parameters such as temperature and rainfall are greatly influenced by global as well as regional geography e.g., presence of mountain, distance from sea and location in relation to the sea.In climate or meteorological data analysis, spatial distribution and pattern bear utmost importance because many other environmental variables [33] that govern ecosystems and their services depend on climate [34].
Geostatistics, which is based on the theory of regionalized variables, is increasingly preferred by the scientific community because it allows for capitalizing on the spatial correlation among neighboring observations to predict attribute values at unsampled locations [35].Geostatistical methods have been shown to be superior to other interpolation/extrapolation methods, such as Thiessen polygon, polynomial interpolation, and inverse distance methods [36,37].One of the important advantages of the geostatistical method is to use additional information.

Kriging-Stochastic/Deterministic Interpolation
Stochastic methods are ones that incorporate randomness.The imaginary interpolated surface is conceptualized as one of many that might have been observed, all of which could have produced the known data points.Kriging interpolation procedures produce trend surface analysis that allows the statistical significance of the surface and uncertainty of the predicted value to be calculated.Daily weather conditions are known to influence the growth and development of many biological organisms.It is therefore surprising that relatively few ecological or environmental studies focus upon the interpolation of continuous national coverage of daily temperature using now well-developed geostatistical or splining techniques (e.g., [38]).
It has been postulated that the actual interpolation algorithm used is less critical than the incorporation of environmentally determined gridded variables to augment the interpolation procedure [39,40], although multiple guiding variables are rarely used in either applied studies or research comparing interpolator performance in studies at a daily time step.Certainly, as reported by Jarvis and Stuart [41], this issue, both [42] and [38] are able to relate patterns of uncertainty in their interpolated temperature surfaces with environmental influences such as proximity to the coast or urban heat islands.
Their conclusions suggest that polynomial regression (trend surface plus covariates) is preferable to other local techniques such as inverse distance weighting.In a study of greater depth for the same area [42], second-degree trend models (including Kriging with External Drift (KED) is a non-stationary geostatistical method.Here, we focus on the use of secondary information from a model to obtain better prediction.
In the case of KED, predictions at new locations are made by: for n ÿ i"1 w KED i pS 0 q ˚qk pS i q " q k pS 0 q ; k " 1, . . ., p where z is the target variable, q k s are the predictor variables, i.e., values at new location pS 0 q , B 0 is the vector of KED weights (w KED i q, P is the number of predictors and z is the vector of n observations at primary locations, [43] External drift kriging algorithm is explained in-depth, which has been used in this study.

Changing Trends of Seasonal Maximum and Minimum Temperature
The observed climatic data from 1978 to 2007 indicate that the average temperature in Bangladesh is increasing for all the four seasons.In Pre-Monsoon (March-May), average maximum and minimum temperature (Figures 3-5) were found to be increasing at a rate of 0.016 ˝C and 0.015 ˝C per year, respectively.Similarly, the rate of increase in the average maximum and minimum were found to be 0.034 ˝C and 0.014 ˝C for Monsoon (June-August), 0.018 ˝C and 0.010 ˝C for Post-Monsoon (September-November) and 0.015 ˝C and 0.009 ˝C for Winter (December-February).
where z is the target variable, q k s are the predictor variables, i.e., values at new location (S 0 ), ∂ 0 is the vector of KED weights ( w i KED ), P is the number of predictors and z is the vector of n observations at primary locations, [43] External drift kriging algorithm is explained in-depth, which has been used in this study.

Changing Trends of Seasonal Maximum and Minimum Temperature
The observed climatic data from 1978 to 2007 indicate that the average temperature in Bangladesh is increasing for all the four seasons.In Pre-Monsoon (March-May), average maximum and minimum temperature (Figures 3-5) were found to be increasing at a rate of 0.016 °C and 0.015 °C per year, respectively.Similarly, the rate of increase in the average maximum and minimum were found to be 0.034 °C and 0.014 °C for Monsoon (June-August), 0.018° C and 0.010 °C for Post-Monsoon (September-November) and 0.015 °C and 0.009 °C for Winter (December-February).It was also revealed that 2006 was the warmest year and 1984 was the coldest year in the last 30 years with a yearly average temperature of 31.1 °C and 20.7 °C.Among the last 30 years, Bangladesh experienced the warmest Pre-Monsoon in 1995 (seasonal average 33.8 °C) and the coldest Winter in the year 1989 (seasonal average 13.1 °C).Reference [10] has studied surface climatological data on monthly and annual mean maximum and minimum temperature, and monthly and annual rainfall for the period of 1961-1990.The study showed an increasing trend of mean maximum and minimum temperature in some seasons and a decreasing trend in some others.The overall trend of the annual mean maximum temperature has shown a significant increase over the period of 1961-1990.Regional variations have been observed around the average trend [7].
Observed data of the Rajshahi Station (from 1978 to 2007) (Figures 4 and 5) indicate that overall annual maximum and minimum temperature are increasing at the rate of 0.008 °C and 0.003 °C , respectively.It was also revealed that 2006 was the warmest year and 1984 was the coldest year in the last 30 years with a yearly average temperature of 31.1 ˝C and 20.7 ˝C.Among the last 30 years, Bangladesh experienced the warmest Pre-Monsoon in 1995 (seasonal average 33.8 ˝C) and the coldest Winter in the year 1989 (seasonal average 13.1 ˝C).Reference [10] has studied surface climatological data on monthly and annual mean maximum and minimum temperature, and monthly and annual rainfall for the period of 1961-1990.The study showed an increasing trend of mean maximum and minimum temperature in some seasons and a decreasing trend in some others.The overall trend of the annual mean maximum temperature has shown a significant increase over the period of 1961-1990.Regional variations have been observed around the average trend [7].
Observed data of the Rajshahi Station (from 1978 to 2007) (Figures 4 and 5) indicate that overall annual maximum and minimum temperature are increasing at the rate of 0.008 ˝C and 0.003 ˝C, respectively.However, rate of change in the Monsoon season is much higher than annual rate of change.The average maximum temperature in the Pre-Monsoon season in Rajshahi exhibited an increasing trend while (Figure 4) Average Maximum Temperature in Pre-Monsoon (Bangladesh) decreasing trend at a rate of 0.014 ˝C per year.This decreasing trend also prevailed for the average maximum temperature of the Post-Monsoon in the region.In Monsoon season (June, July and August), the average maximum and minimum temperature show an increasing trend annually at the rate of 0.032 ˝C and 0.013 ˝C, respectively, which means Monsoon season is becoming warmer.On the other hand, average minimum temperature in Winter season (December, January and February) shows almost no change (0.004 ˝C/year) while maximum temperature shows a slight increasing trend annually at the rate of 0.010 ˝C (Figure 3).However, rate of change in the Monsoon season is much higher than annual rate of change.The average maximum temperature in the Pre-Monsoon season in Rajshahi exhibited an increasing trend while (Figure 4) Average Maximum Temperature in Pre-Monsoon (Bangladesh) decreasing trend at a rate of 0.014 °C per year.This decreasing trend also prevailed for the average maximum temperature of the Post-Monsoon in the region.In Monsoon season (June, July and August), the average maximum and minimum temperature show an increasing trend annually at the rate of 0.032 °C and 0.013 °C, respectively, which means Monsoon season is becoming warmer.On the other hand, average minimum temperature in Winter season (December, January and February) shows almost no change (0.004 °C/year) while maximum temperature shows a slight increasing trend annually at the rate of 0.010 °C (Figure 3).Observations from Barisal meteorological station revealed both maximum and minimum average temperature in this region increasing annually (Figures 4 and 5).This region experienced comparatively hotter Monsoon and Winter seasons than in the past.The highest rate, 0.030 °C/year was recorded from 1978 to 2007 observations (30 years) for the average maximum temperature in the Monsoon season with a similar increasing rate in Winter 0.026 °C/year and 0.027 °C/year for average maximum and minimum temperature, respectively In Chittagong, the increasing trend in temperature is much higher than the other regions of Bangladesh.In Pre-Monsoon, the increase rate is 0.049 °C/year and 0.038 °C/year for average Observations from Barisal meteorological station revealed both maximum and minimum average temperature in this region increasing annually (Figures 4 and 5).This region experienced comparatively hotter Monsoon and Winter seasons than in the past.The highest rate, 0.030 ˝C/year was recorded from 1978 to 2007 observations (30 years) for the average maximum temperature in the Monsoon season with a similar increasing rate in Winter 0.026 ˝C/year and 0.027 ˝C/year for average maximum and minimum temperature, respectively In Chittagong, the increasing trend in temperature is much higher than the other regions of Bangladesh.In Pre-Monsoon, the increase rate is 0.049 ˝C/year and 0.038 ˝C/year for average maximum and minimum temperature, respectively.In Monsoon and Post-Monsoon, the rate of increase rate is relatively lower than Pre-Monsoon and Winter season.The Winter season of Chittagong is getting warmer at a rate of 0.040 ˝C/year.The average maximum temperature is increasing at a rate of 0.067 ˝C/year and the average minimum temperature is increasing at 0.031 ˝C/year.During the last 30 years, the highest average maximum temperature in Winter was 29.1 ˝C in the year 2002, and the lowest average minimum temperature was 10 ˝C in the year 1989.
Sylhet, the northeast region of Bangladesh is also experiencing the impact of global warming all year round.The increasing trend in temperature is a comparatively higher (0.061 ˝C/year) average minimum temperature in Winter season.The remaining observations exhibited a similar warming trend (0.031 ˝C/year to 0.05 ˝C/year) with little exception in average minimum temperature in the Monsoon (0.02 ˝C/year).
All four regions of Bangladesh (Rajshahi, Barisal, Chittagong and Sylhet) recorded a positive trend in temperature over the past 30 years (1978-2008) of observations except little decrease in average maximum temperature in Pre-and Post-Monsoon seasons in the Rajshahi region (Figures 5-7).The past 30 years of observations showed that average temperature is increasing comparatively sharply in the Chittagong and Sylhet regions in comparison with the Barisal and Rajshahi regions.On average, the temperature increase rate in Winter is much more rapid than other seasons in Bangladesh.maximum and minimum temperature, respectively.In Monsoon and Post-Monsoon, the rate of increase rate is relatively lower than Pre-Monsoon and Winter season.The Winter season of Chittagong is getting warmer at a rate of 0.040 °C/year.The average maximum temperature is increasing at a rate of 0.067 °C/year and the average minimum temperature is increasing at 0.031 °C/year.During the last 30 years, the highest average maximum temperature in Winter was 29.1 °C in the year 2002, and the lowest average minimum temperature was 10 °C in the year 1989.
Sylhet, the northeast region of Bangladesh is also experiencing the impact of global warming all year round.The increasing trend in temperature is a comparatively higher (0.061 °C/year) average minimum temperature in Winter season.The remaining observations exhibited a similar warming trend (0.031 °C/year to 0.05 °C/year) with little exception in average minimum temperature in the Monsoon (0.02 °C/year).
All four regions of Bangladesh (Rajshahi, Barisal, Chittagong and Sylhet) recorded a positive trend in temperature over the past 30 years (1978-2008) of observations except little decrease in average maximum temperature in Pre-and Post-Monsoon seasons in the Rajshahi region (Figures 5-7).The past 30 years of observations showed that average temperature is increasing comparatively sharply in the Chittagong and Sylhet regions in comparison with the Barisal and Rajshahi regions.On average, the temperature increase rate in Winter is much more rapid than other seasons in Bangladesh.

Changing Trends in Total and Monthly Average Rainfall
Bangladesh has a subtropical monsoon climate characterized by wide seasonal variations in rainfall.Bangladesh gets maximum rainfall in the monsoon season and minimum rainfall in the dry and cool Winter.The past 30 years (1978-2007) of observations showed that the mean rainfall during the Monsoon is about 1391 ˘189 mm having maximum rainfall 1808 mm in 1984 and 1076 mm in 1992.In Winter, maximum rainfall was recorded as 93 mm in 1981 with a minimum about 3 mm in 2004.Recorded rainfall data revealed that average rainfall in Bangladesh is increasing at a rate of 3.346 mm/year in Monsoon and 8.084 mm/year in Post-Monsoon.In Winter and the Pre-Monsoon, average rainfall trend is negative.
In Rajshahi, the rainfall pattern for the last 30 years was similar in Pre-Monsoon, Monsoon and Post-Monsoon seasons.In Winter season, rainfall is showing a negative trend (´1.85 mm/year).The Rajshahi meteorology station recorded an average of 1510 ˘306 mm rainfall per year.In 1981, this region received 2241 mm (maximum) and 843 mm (minimum) in 1992.
Rainfall records from Barisal illustrated that the total rainfall is decreasing in Pre-Monsoon (´3.978 mm/year) (Figure 8a), Monsoon (´7.0049 mm/year) (Figure 8b) and in Winter (´0.3519 mm/year) (Figure 9b).In the Post-Monsoon, rainfall is found to be increasing at a rate of 8.6912 mm/year (Figures 8 and 9).Chittagong, the hilly coastal region of Bangladesh, is experiencing an increasing trend in total rainfall except in the Monsoon season.It may be due to shifting of Monsoon onset/withdrawal periods of this region (need detailed study).Among the considered last 30 years observations, monsoon in Chittagong got the highest rainfall (2878 mm) in 1983 and the lowest (623 mm) in 2001.

Changing Trends in Total and Monthly Average Rainfall
Bangladesh has a subtropical monsoon climate characterized by wide seasonal variations in rainfall.Bangladesh gets maximum rainfall in the monsoon season and minimum rainfall in the dry and cool Winter.The past 30 years (1978-2007) of observations showed that the mean rainfall during the Monsoon is about 1391 ± 189 mm having maximum rainfall 1808 mm in 1984 and 1076 mm in 1992.In Winter, maximum rainfall was recorded as 93 mm in 1981 with a minimum about 3 mm in 2004.Recorded rainfall data revealed that average rainfall in Bangladesh is increasing at a rate of 3.346 mm/year in Monsoon and 8.084 mm/year in Post-Monsoon.In Winter and the Pre-Monsoon, average rainfall trend is negative.
In Rajshahi, the rainfall pattern for the last 30 years was similar in Pre-Monsoon, Monsoon and Post-Monsoon seasons.In Winter season, rainfall is showing a negative trend (−1.85 mm/year).The Rajshahi meteorology station recorded an average of 1510 ± 306 mm rainfall per year.In 1981, this region received 2241 mm (maximum) and 843 mm (minimum) in 1992.
Rainfall records from Barisal illustrated that the total rainfall is decreasing in Pre-Monsoon (−3.978 mm/year) (Figure 8a), Monsoon (−7.0049 mm/year) (Figure 8b) and in Winter (−0.3519 mm/year) (Figure 9b).In the Post-Monsoon, rainfall is found to be increasing at a rate of 8.6912 mm/year (Figures 8 and 9).Chittagong, the hilly coastal region of Bangladesh, is experiencing an increasing trend in total rainfall except in the Monsoon season.It may be due to shifting of Monsoon onset/withdrawal periods of this region (need detailed study).Among the considered last 30  Rainfall records from Barisal illustrated that the total rainfall is decreasing in Pre-Monsoon (−3.978 mm/year) (Figure 8a), Monsoon (−7.0049 mm/year) (Figure 8b) and in Winter (−0.3519 mm/year) (Figure 9b).In the Post-Monsoon, rainfall is found to be increasing at a rate of 8.6912 mm/year (Figures 8 and 9).Chittagong, the hilly coastal region of Bangladesh, is experiencing an increasing trend in total rainfall except in the Monsoon season.It may be due to shifting of Monsoon onset/withdrawal periods of this region (need detailed study).Among the considered last 30  Among the last five years (2003)(2004)(2005)(2006)(2007), this region received less rainfall in three years (2004, 2005, and 2006) than the 30-year Monsoon average.Though the Monsoon in Sylhet is showing an increasing trend in total rainfall (4.117 mm/year), Pre-Monsoon and Post-Monsoon recorded negative trends, and, in Winter, the trend is almost flat.In year-based analysis, rainfall pattern in this region showed a decreasing trend, at a rate of 6.252 mm/year with an average rainfall at about 4073 (˘598) mm/year.
Among the four regions (Rajshahi, Barisal, Chittagong and Sylhet), only Chittagong showed an increasing trend in rainfall in the Winter season.In Post-Monsoon, a decreasing trend in rainfall is observed only in the Sylhet region.Analysis showed that in both Pre-Monsoon and Monsoon seasons, rainfall is increasing in Rajshahi and decreasing in Barisal.

Temperature-Rainfall Inter-Linkage
Temperature and rainfall are the most prominent parameters in climate study.With global climate change, temperature and rainfall are changing their pattern and spatial distribution.This present research finds a significant linkage between temperature and rainfall in Pre-Monsoon and Winter seasons of Bangladesh.However, in Monsoon and Post-Monsoon, the relationship was found to be insignificant (Table 1).For an increase of 1 ˝C, the average maximum temperature in Pre-Monsoon rainfall will decrease by 153 mm.Similarly, a 1 ˝C increase in average maximum temperature may cause a decrease of 23 mm in rainfall in Winter season.For average minimum temperature, an increase of 1 ˝C may reduce by 127 mm in rainfall in the Pre-Monsoon season.Similarly, for seasonal average temperature, an increase of 1 ˝C may cause a decrease of 170 mm rainfall in Pre-Monsoon and 17 mm in Winter season.

Discussion
The geospatial modeling results revealed that the northeastern Sylhet region has clearly shown a discontinuity in the Monsoon progression (southeast to northwest) across Bangladesh.On the contrary, the northeast region compared to the usual Monsoon progression from the southeast to the northwest showed a peculiar early onset.This peculiarity in Monsoon in the northeast region has also been reported [44,45].Reference [46] has also postulated that this peculiar early onset clearly depicts a smaller scale location specific Monsoon mechanism that influences the rain and climate in northeast Bangladesh.This less understood location specific micro climate process could be attributed to the region's close proximity to highly complex topography of the Himalayan mountain system, where the height of mountaintops changes by about 2000 m over a distance of less than 50 km [45].This abrupt change in topography creates Monsoons through northeast Bangladesh and adjacent India.This strength of the Mann-Kendall test in finding the seasonal and locational processes in geospatial pattern of climate parameters were also shown by [47] and [48].Please see Section 2.3.2 for further deliberation on these.
Several previous studies have shown trends of climate change in climatic parameters over Bangladesh.Chowdhury and Debsharma [6] and Mia [7] found out that temperature has been increasing using historical data of some selected meteorological stations.Reference [7] and [8] reported that mean annual temperature of Bangladesh has increased during the period of 1895-1980 at 0.31 ˝C over the past two decades, and the trends show a close match to the findings of this research.Karmakar and Shrestha [10] projected that annual mean maximum temperature will increase to 0.4 ˝C and 0.73 ˝C by the year of 2050 and 2100, respectively.The differences in trends could be attributed to different time periods used by different research and approaches as well.
The Mann-Kendall test showed significant linkage between temperature and rainfall in Pre-Monsoon and Winter seasons, which is substantiated by the studies in [3,46,49,50].The southeastern region, which is one the most important forest ecosystem zone in the country is experiencing early onset and withdrawal of rain, yet increasing trends in total rainfall except in the Monsoon season [30].This means that the region is experiencing a lower number of rainy days.However, total rainfall has not changed significantly.This changing pattern in average maximum and minimum temperature along with precipitation might cause a situation in which the species which are growing now may shift to suitable habitats elsewhere in the future [51].

Conclusions
The geospatial modeling and perception survey results from Sylhet clearly show a discontinuity in the Monsoon progression (southeast to northwest) across Bangladesh.This divergence of Monsoon progression and a peculiar early onset in the northeast region implies that there is a smaller scale Monsoon mechanism influencing the rain and climate in northeast Bangladesh, and this is peculiar and less understood.Understanding this small mechanism is important for better analysis of GCM (Global Circulation/Climate Model) and RAM (Regional Atmospheric Model) results.For example, if these climate models do not resolve these small mechanisms, early Monsoon onset in this region, RAM/RCM (Regional Climate Model) may give results that will be difficult to explain.Understanding the micro process is also important because it would have direct implication on the tropical ecosystem with rich forests and wetland habitats.More importantly, the region has been historically gifted with the highest amount of rainfall in Bangladesh.These climatic conditions are imperative for sustainable tea production, which is one of the major industries and employment sources.In addition, the model scenarios for the eastern part of the country, much of which is a hilly and coastal region, experienced the impacts of warming all the year round, which is here further substantiated by these findings from a community level perception survey.The region is experiencing early onset and withdrawal of rain, yet increasing trends in total rainfall except in the Monsoon season.This means the region is experiencing fewer numbers of rainy days.However, total rainfall has not changed significantly.This implies that at the beginning of the Monsoon, the region is experiencing more rainfall and hence flooding, which may be associated with some accompanying landslides in the region.In addition to that incremental gap between Tmax and Tmin, temperatures may appear as inhibiting factors for germination and natural propagation of different plant species in the southeast region of Bangladesh.The impact of this changing climate would be varied on NPP (Net Primary Production) according to vegetation types.

Figure 1
Figure 1 presents conceptual and sequential analytical process flows of research.It is important to note that this paper is part of a larger piece of research and only the geospatial modeling part of the research is presented here.

Figure 1 .
Figure 1.Schematic diagram shows the conceptual framework and sequential process of the research.

Figure 1 .
Figure 1.Schematic diagram shows the conceptual framework and sequential process of the research.

Climate 2016, 4 , 21 5 of 16 Figure 2 .
Figure 2. Geostatistical Kriging Mess for interpolation of missing values for the study area.

Figure 2 .
Figure 2. Geostatistical Kriging Mess for interpolation of missing values for the study area.

Figure 4 .
Figure 4. Average mean annual minimum temperature in different seasons.

Figure 5 .
Figure 5. Geospatial distribution changing tends at seasonal average minimum temperature.

Figure 4 .
Figure 4. Average mean annual minimum temperature in different seasons.

Figure 4 .
Figure 4. Average mean annual minimum temperature in different seasons.

Figure 5 .
Figure 5. Geospatial distribution changing tends at seasonal average minimum temperature.

Figure 5 .
Figure 5. Geospatial distribution changing tends at seasonal average minimum temperature.

Figure 6 .
Figure 6.Geospatial distribution changing tends in seasonal average maximum temperature.

Figure 6 .
Figure 6.Geospatial distribution changing tends in seasonal average maximum temperature.

Figure 7 .
Figure 7. Geospatial distribution of changing trends in seasonal average rainfall.

Figure 7 .
Figure 7. Geospatial distribution of changing trends in seasonal average rainfall.

Table 1 .
Interrelation between temperature and rainfall in different seasons in Bangladesh.