\documentclass[11pt,twoside]{article}\makeatletter

\IfFileExists{xcolor.sty}%
  {\RequirePackage{xcolor}}%
  {\RequirePackage{color}}
\usepackage{colortbl}
\usepackage{wrapfig}
\usepackage{ifxetex}
\ifxetex
  \usepackage{fontspec}
  \usepackage{xunicode}
  \catcode`⃥=\active \def⃥{\textbackslash}
  \catcode`❴=\active \def❴{\{}
  \catcode`❵=\active \def❵{\}}
  \def\textJapanese{\fontspec{Noto Sans CJK JP}}
  \def\textChinese{\fontspec{Noto Sans CJK SC}}
  \def\textKorean{\fontspec{Noto Sans CJK KR}}
  \setmonofont{DejaVu Sans Mono}
  
\else
  \IfFileExists{utf8x.def}%
   {\usepackage[utf8x]{inputenc}
      \PrerenderUnicode{–}
    }%
   {\usepackage[utf8]{inputenc}}
  \usepackage[english]{babel}
  \usepackage[T1]{fontenc}
  \usepackage{float}
  \usepackage[]{ucs}
  \uc@dclc{8421}{default}{\textbackslash }
  \uc@dclc{10100}{default}{\{}
  \uc@dclc{10101}{default}{\}}
  \uc@dclc{8491}{default}{\AA{}}
  \uc@dclc{8239}{default}{\,}
  \uc@dclc{20154}{default}{ }
  \uc@dclc{10148}{default}{>}
  \def\textschwa{\rotatebox{-90}{e}}
  \def\textJapanese{}
  \def\textChinese{}
  \IfFileExists{tipa.sty}{\usepackage{tipa}}{}
\fi
\def\exampleFont{\ttfamily\small}
\DeclareTextSymbol{\textpi}{OML}{25}
\usepackage{relsize}
\RequirePackage{array}
\def\@testpach{\@chclass
 \ifnum \@lastchclass=6 \@ne \@chnum \@ne \else
  \ifnum \@lastchclass=7 5 \else
   \ifnum \@lastchclass=8 \tw@ \else
    \ifnum \@lastchclass=9 \thr@@
   \else \z@
   \ifnum \@lastchclass = 10 \else
   \edef\@nextchar{\expandafter\string\@nextchar}%
   \@chnum
   \if \@nextchar c\z@ \else
    \if \@nextchar l\@ne \else
     \if \@nextchar r\tw@ \else
   \z@ \@chclass
   \if\@nextchar |\@ne \else
    \if \@nextchar !6 \else
     \if \@nextchar @7 \else
      \if \@nextchar (8 \else
       \if \@nextchar )9 \else
  10
  \@chnum
  \if \@nextchar m\thr@@\else
   \if \@nextchar p4 \else
    \if \@nextchar b5 \else
   \z@ \@chclass \z@ \@preamerr \z@ \fi \fi \fi \fi
   \fi \fi  \fi  \fi  \fi  \fi  \fi \fi \fi \fi \fi \fi}
\gdef\arraybackslash{\let\\=\@arraycr}
\def\@textsubscript#1{{\m@th\ensuremath{_{\mbox{\fontsize\sf@size\z@#1}}}}}
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\def\abbr{}
\def\corr{}
\def\expan{}
\def\gap{}
\def\orig{}
\def\reg{}
\def\ref{}
\def\sic{}
\def\persName{}\def\name{}
\def\placeName{}
\def\orgName{}
\def\textcal#1{{\fontspec{Lucida Calligraphy}#1}}
\def\textgothic#1{{\fontspec{Lucida Blackletter}#1}}
\def\textlarge#1{{\large #1}}
\def\textoverbar#1{\ensuremath{\overline{#1}}}
\def\textquoted#1{‘#1’}
\def\textsmall#1{{\small #1}}
\def\textsubscript#1{\@textsubscript{\selectfont#1}}
\def\textxi{\ensuremath{\xi}}
\def\titlem{\itshape}
\newenvironment{biblfree}{}{\ifvmode\par\fi }
\newenvironment{bibl}{}{}
\newenvironment{byline}{\vskip6pt\itshape\fontsize{16pt}{18pt}\selectfont}{\par }
\newenvironment{citbibl}{}{\ifvmode\par\fi }
\newenvironment{docAuthor}{\ifvmode\vskip4pt\fontsize{16pt}{18pt}\selectfont\fi\itshape}{\ifvmode\par\fi }
\newenvironment{docDate}{}{\ifvmode\par\fi }
\newenvironment{docImprint}{\vskip 6pt}{\ifvmode\par\fi }
\newenvironment{docTitle}{\vskip6pt\bfseries\fontsize{22pt}{25pt}\selectfont}{\par }
\newenvironment{msHead}{\vskip 6pt}{\par}
\newenvironment{msItem}{\vskip 6pt}{\par}
\newenvironment{rubric}{}{}
\newenvironment{titlePart}{}{\par }

\newcolumntype{L}[1]{){\raggedright\arraybackslash}p{#1}}
\newcolumntype{C}[1]{){\centering\arraybackslash}p{#1}}
\newcolumntype{R}[1]{){\raggedleft\arraybackslash}p{#1}}
\newcolumntype{P}[1]{){\arraybackslash}p{#1}}
\newcolumntype{B}[1]{){\arraybackslash}b{#1}}
\newcolumntype{M}[1]{){\arraybackslash}m{#1}}
\definecolor{label}{gray}{0.75}
\def\unusedattribute#1{\sout{\textcolor{label}{#1}}}
\DeclareRobustCommand*{\xref}{\hyper@normalise\xref@}
\def\xref@#1#2{\hyper@linkurl{#2}{#1}}
\begingroup
\catcode`\_=\active
\gdef_#1{\ensuremath{\sb{\mathrm{#1}}}}
\endgroup
\mathcode`\_=\string"8000
\catcode`\_=12\relax

\usepackage[a4paper,twoside,lmargin=1in,rmargin=1in,tmargin=1in,bmargin=1in,marginparwidth=0.75in]{geometry}
\usepackage{framed}

\definecolor{shadecolor}{gray}{0.95}
\usepackage{longtable}
\usepackage[normalem]{ulem}
\usepackage{fancyvrb}
\usepackage{fancyhdr}
\usepackage{graphicx}
\usepackage{marginnote}

\renewcommand{\@cite}[1]{#1}


\renewcommand*{\marginfont}{\itshape\footnotesize}

\def\Gin@extensions{.pdf,.png,.jpg,.mps,.tif}

  \pagestyle{fancy}

\usepackage[pdftitle={Determining Temperature Trends in the Granary Areas of Peninsular Malaysia using Mann-Kendall and Sen's Slope Estimator},
 pdfauthor={}]{hyperref}
\hyperbaseurl{}

	 \paperwidth210mm
	 \paperheight297mm
              
\def\@pnumwidth{1.55em}
\def\@tocrmarg {2.55em}
\def\@dotsep{4.5}
\setcounter{tocdepth}{3}
\clubpenalty=8000
\emergencystretch 3em
\hbadness=4000
\hyphenpenalty=400
\pretolerance=750
\tolerance=2000
\vbadness=4000
\widowpenalty=10000

\renewcommand\section{\@startsection {section}{1}{\z@}%
     {-1.75ex \@plus -0.5ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large\bfseries}}
\renewcommand\subsection{\@startsection{subsection}{2}{\z@}%
     {-1.75ex\@plus -0.5ex \@minus- .2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large}}
\renewcommand\subsubsection{\@startsection{subsubsection}{3}{\z@}%
     {-1.5ex\@plus -0.35ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\large}}
\renewcommand\paragraph{\@startsection{paragraph}{4}{\z@}%
     {-1ex \@plus-0.35ex \@minus -0.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\normalsize}}
\renewcommand\subparagraph{\@startsection{subparagraph}{5}{\parindent}%
     {1.5ex \@plus1ex \@minus .2ex}%
     {-1em}%
     {\reset@font\normalsize\bfseries}}


\def\l@section#1#2{\addpenalty{\@secpenalty} \addvspace{1.0em plus 1pt}
 \@tempdima 1.5em \begingroup
 \parindent \z@ \rightskip \@pnumwidth 
 \parfillskip -\@pnumwidth 
 \bfseries \leavevmode #1\hfil \hbox to\@pnumwidth{\hss #2}\par
 \endgroup}
\def\l@subsection{\@dottedtocline{2}{1.5em}{2.3em}}
\def\l@subsubsection{\@dottedtocline{3}{3.8em}{3.2em}}
\def\l@paragraph{\@dottedtocline{4}{7.0em}{4.1em}}
\def\l@subparagraph{\@dottedtocline{5}{10em}{5em}}
\@ifundefined{c@section}{\newcounter{section}}{}
\@ifundefined{c@chapter}{\newcounter{chapter}}{}
\newif\if@mainmatter 
\@mainmattertrue
\def\chaptername{Chapter}
\def\frontmatter{%
  \pagenumbering{roman}
  \def\thechapter{\@roman\c@chapter}
  \def\theHchapter{\roman{chapter}}
  \def\thesection{\@roman\c@section}
  \def\theHsection{\roman{section}}
  \def\@chapapp{}%
}
\def\mainmatter{%
  \cleardoublepage
  \def\thechapter{\@arabic\c@chapter}
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \pagenumbering{arabic}
  \setcounter{secnumdepth}{6}
  \def\@chapapp{\chaptername}%
  \def\theHchapter{\arabic{chapter}}
  \def\thesection{\@arabic\c@section}
  \def\theHsection{\arabic{section}}
}
\def\backmatter{%
  \cleardoublepage
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \setcounter{secnumdepth}{2}
  \def\@chapapp{\appendixname}%
  \def\thechapter{\@Alph\c@chapter}
  \def\theHchapter{\Alph{chapter}}
  \appendix
}
\newenvironment{bibitemlist}[1]{%
   \list{\@biblabel{\@arabic\c@enumiv}}%
       {\settowidth\labelwidth{\@biblabel{#1}}%
        \leftmargin\labelwidth
        \advance\leftmargin\labelsep
        \@openbib@code
        \usecounter{enumiv}%
        \let\p@enumiv\@empty
        \renewcommand\theenumiv{\@arabic\c@enumiv}%
	}%
  \sloppy
  \clubpenalty4000
  \@clubpenalty \clubpenalty
  \widowpenalty4000%
  \sfcode`\.\@m}%
  {\def\@noitemerr
    {\@latex@warning{Empty `bibitemlist' environment}}%
    \endlist}

\def\tableofcontents{\section*{\contentsname}\@starttoc{toc}}
\parskip0pt
\parindent1em
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\newenvironment{reflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemsep}{0pt}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\itshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{sansreflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\itemsep}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\upshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{specHead}[2]%
 {\vspace{20pt}\hrule\vspace{10pt}%
  \phantomsection\label{#1}\markright{#2}%

  \pdfbookmark[2]{#2}{#1}%
  \hspace{-0.75in}{\bfseries\fontsize{16pt}{18pt}\selectfont#2}%
  }{}
      \def\TheFullDate{2020-01-15 (revised: 15 January 2020)}
\def\TheID{\makeatother }
\def\TheDate{2020-01-15}
\title{Determining Temperature Trends in the Granary Areas of Peninsular Malaysia using Mann-Kendall and Sen's Slope Estimator}
\author{}\makeatletter 
\makeatletter
\newcommand*{\cleartoleftpage}{%
  \clearpage
    \if@twoside
    \ifodd\c@page
      \hbox{}\newpage
      \if@twocolumn
        \hbox{}\newpage
      \fi
    \fi
  \fi
}
\makeatother
\makeatletter
\thispagestyle{empty}
\markright{\@title}\markboth{\@title}{\@author}
\renewcommand\small{\@setfontsize\small{9pt}{11pt}\abovedisplayskip 8.5\p@ plus3\p@ minus4\p@
\belowdisplayskip \abovedisplayskip
\abovedisplayshortskip \z@ plus2\p@
\belowdisplayshortskip 4\p@ plus2\p@ minus2\p@
\def\@listi{\leftmargin\leftmargini
               \topsep 2\p@ plus1\p@ minus1\p@
               \parsep 2\p@ plus\p@ minus\p@
               \itemsep 1pt}
}
\makeatother
\fvset{frame=single,numberblanklines=false,xleftmargin=5mm,xrightmargin=5mm}
\fancyhf{} 
\setlength{\headheight}{14pt}
\fancyhead[LE]{\bfseries\leftmark} 
\fancyhead[RO]{\bfseries\rightmark} 
\fancyfoot[RO]{}
\fancyfoot[CO]{\thepage}
\fancyfoot[LO]{\TheID}
\fancyfoot[LE]{}
\fancyfoot[CE]{\thepage}
\fancyfoot[RE]{\TheID}
\hypersetup{citebordercolor=0.75 0.75 0.75,linkbordercolor=0.75 0.75 0.75,urlbordercolor=0.75 0.75 0.75,bookmarksnumbered=true}
\fancypagestyle{plain}{\fancyhead{}\renewcommand{\headrulewidth}{0pt}}

\date{}
\usepackage{authblk}

\providecommand{\keywords}[1]
{
\footnotesize
  \textbf{\textit{Index terms---}} #1
}

\usepackage{graphicx,xcolor}
\definecolor{GJBlue}{HTML}{273B81}
\definecolor{GJLightBlue}{HTML}{0A9DD9}
\definecolor{GJMediumGrey}{HTML}{6D6E70}
\definecolor{GJLightGrey}{HTML}{929497} 

\renewenvironment{abstract}{%
   \setlength{\parindent}{0pt}\raggedright
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
   \textcolor{GJBlue}{\large\bfseries\abstractname\space}
}{%   
   \vskip8pt
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
}

\usepackage[absolute,overlay]{textpos}

\makeatother 
      \usepackage{lineno}
      \linenumbers
      
\begin{document}

             \author[1]{Danladi Yusuf  Gumel}

             \author[2]{Ahmad Makmom Bn  Abdullah}

             \author[3]{Rasheeda E.  Elhadi}

             \author[4]{Da?u Abba  Umar}

             \affil[1]{  Jigawa State College of Education}

\renewcommand\Authands{ and }

\date{\small \em Received: 14 December 2019 Accepted: 1 January 2020 Published: 15 January 2020}

\maketitle


\begin{abstract}
        


The spatiotemporal dynamics of temperature as well as rainfall have received greater attention from the scientific communities. This study analysed temperature variability in the three granary areas of Peninsular Malaysia using descriptive statistics, parametric (least square regression) and nonparametric (Mann- Kendall and Sen?s slope estimator). The study identified significant warming trend in the annual mean maximum temperature in two of the study areas, i.e. Subang Jaya and Kota Bharu. Also significant warming trend was detected in the annual minimum temperature and significant increasing trend in some of the monthly maximum and minimum temperatures for all the three stations. Also the result reveals spatial and temporal variation in both the maximum and minimum temperature at annual, monthly and seasonal scales. For the annual scale maximum temperature, this study identified a warming trend for the two stations with about 0.014oC per year (1.4oC per 100 years).

\end{abstract}


\keywords{temperature trends, mann- kendall, sen?s slope estimator, peninsular malaysia.}

\begin{textblock*}{18cm}(1cm,1cm) % {block width} (coords) 
\textcolor{GJBlue}{\LARGE Global Journals \LaTeX\ JournalKaleidoscope\texttrademark}
\end{textblock*}

\begin{textblock*}{18cm}(1.4cm,1.5cm) % {block width} (coords) 
\textcolor{GJBlue}{\footnotesize \\ Artificial Intelligence formulated this projection for compatibility purposes from the original article published at Global Journals. However, this technology is currently in beta. \emph{Therefore, kindly ignore odd layouts, missed formulae, text, tables, or figures.}}
\end{textblock*}


\let\tabcellsep& 	 	 		 
\section[{Introduction}]{Introduction}\par
here has been a growing concern from scientific communities across the globe on analyzing spatiotemporal dynamics of rainfall and temperatures. As these two critical weather elements exerts overriding influence on agriculture and other aspects of human society, triggered by the increase in anthropogenic greenhouse gas emission causing the warming of the globe \hyperref[b15]{(IPCC, 2007)}. The global warming phenomenon creates series of feedback mechanisms affecting the natural processes of the hydrologic cycles which alters rainfall patterns in terms of intensity, duration, frequency, and onset and cessation. Although the issue of climate change has become a global issue, yet its impacts are rather deterministic in nature and depends on the specific region of concern \hyperref[b0]{(Amirabadizadeh , Huang, \& Shui Lee, 2015)}. Global temperature studies have clearly reported on the recent long-term warming of the global around the middle of 1970 to 2013 with an average global trend of 0.2 ºC per decade \hyperref[b32]{(Rohde et al., 2013;}\hyperref[b44]{Turco, Palazzi, Hardenberg, \& Provenzale, 2015)}.\par
Though, in recent times the global temperature trends all over the world have either declined or lacks statistical significance in a lot of regional series around the globe around 1997 \hyperref[b20]{(Kaufmann, Kauppi, Mann, \& Stock, 2011)}, although temperatures have consistently hovers above the long-term averages \hyperref[b1]{(Campra \& Morales 2016)}. Recent observations and global averages show a significant decrease in the warming trend from 0.12 °C per decades between 1951-2012 to 0.05 °C per decades between 1998 to 2012  {\ref (Hartmann et al., 2013)}. Moreover, surface air temperatures exhibits greater spatial and temporal variability \hyperref[b23]{(Lovejoy, 2014;}\hyperref[b36]{Steinman, Mann, \& Miller, 2015;}\hyperref[b44]{Turco et al., 2015)}.\par
Previous studies conducted to analyse the trends in the temperature and rainfall used Mann -Kendall statistics and Theil Sen's slope mainly due to the simplicity and versatility of the approach \hyperref[b25]{(Mustapha, 2013)}. For instance studies in India, by \hyperref[b16]{Jain and Kumar (2012)} studied trends in rainfall, rainy days and temperature over India using Sen's non-parametric estimator and Mann-Kendall test. Their findings showed inconsistent rainfall trends amongst the stations under their study, 15 basins indicated decreasing trend with only one station showing statistically significant trend at 95\% confidence level, while the mean maximum temperature series showed a rising trend for most of the stations; it showed a falling trend at some stations. The mean minimum temperature showed a rising as well as a falling trend. Also in the north-eastern United States, \hyperref[b19]{(Karmeshu, 2012)} similarly used Mann-Kendall test on annual temperature and precipitation for the nine states, their findings revealed statistically significant increasing trend in temperatures for seven out of the nine states with annual linear trend ranging from 0.0006 to 0.02 o F per annum.\par
In a more recent studies, Chakraborty et al. (2017) studied changes in mean air temperature in the parts of eastern Himalaya, in the northeast Indian states. They observed spatial variability in trends with statistically significant increase in annual mean temperature for most of their stations. Despite the spatial variability, the overall range of increase in mean temperature is 0.2 °C to 1.6 °C per decade across the study region. Significant rise in average temperature during the winter is experienced by five out of seven places. While in Austria Herath and Sarukkalige (2017) using binning technique have reported variation of rainfall-temperature scaling with location and reported increasing trend for more extreme short span precipitation and a decreasing trend in the average long span precipitation incidences. The Fifth Assessment Report (AR5) of the intergovernmental panel on climate change (IPCC) revealed that the period 2016 to 2035 will witness increase in the mean temperature around the globe \hyperref[b37]{(Stocker et al., 2013)}. As such the need for the determination of local as well sub-regional trends and variability for a better understanding of the pattern of changes in the world climate has been reiterated \hyperref[b1]{(Campra \& Morales 2016)}.\par
In spite of this growing concern and the prevalence of the trends associated with a number of climatic variables all over the world, studies of this nature in Malaysia are limited in number \hyperref[b0]{(Amirabadizadeh et al., 2015)}. A study by \hyperref[b48]{Wai, Camerlengo, Khairi, and Wahab (2005)} reported an upward trend in the average yearly temperature. Their studies predicted temperature changes ranging from 0.99 o C to 3.44 o C for the next one hundred years. They also reported upward warming trend for all their stations for the past three decades. Tangang, Juneng, and Ahmad (2007), also reported a warming trend of 2.7 -4.0 o C per 100 years for all regions in Peninsular Malaysia and the northern Borneo. Studies by (Suhaila, Deni, Zin, \& Jemain, 2010; Varikoden, Preethi, Samah, \& Babu, 2011; Zin, Jamaludin, Deni, \& Jemain, 2010) have equally demonstrated the variability and occurrence of extreme events signalling changes in the pattern of climate in some parts of the country with its attendant consequences on the environment and other activities.\par
Similarly, the study by Amirabadizadeh et al. (2015) reported a warming trend of minimum and maximum temperatures ranging from 3.5 o C to 4.0 o C per 100 years (0.035 o C to 0.04 o C per annum) over the Langat River Basin in Malaysia. Suhaila and Yusop (2017) used Pettit and sequential Mann-Kendall (SQ-MK) tests to examine the annual and seasonal trends, and change point detection associated with the mean, maximum and minimum temperature data series in Peninsular Malaysia. Their findings detected abrupt changes in the data series and observed significant increasing trends in the annual and seasonal mean, maximum and minimum temperatures in over the country ranging from 2 to 5 °C per 100 years during the last 32 years. They detected large increase in magnitudes of the minimum temperature trend greater than that of the maximum temperatures for most the stations under study.\par
A number of possible factors governing spatial and temporal variability of surface temperature warming in a particular region or continental area have been highlighted in the literature. They include those related variability arising from atmosphere-ocean interface. For example North Atlantic Oscillation (NAO), the Pacific Decadal Oscillation (PDO), the Indian Ocean Dipole (IOD) and the El Ni-no-Southern Oscillation (ENSO).\par
The variability of inter-annual temperatures perhaps have largely been ascribed to the effects of El Ni-no -Southern Oscillation (ENSO) and El-Nina \hyperref[b39]{(Suhaila and Yusop (2017)}. In parts of Southeast Asia recent study by  {\ref Malik et al. (2012)} visualized the pattern of extreme events in the rainfall-fields revealed the variability in the moisture sinks over the region. These natural variabilities are believed to have caused both inter-annual and decadal variations of temperatures in some part of the globe and may influence the long term warming trend in any part of the world \hyperref[b42]{(Tangang et al., 2007)}, as well rainfall variability \hyperref[b41]{(Tan, Ibrahim, Cracknell, \& Yusop, 2017)}.\par
Beside these, other anthropogenic ally induced variability, such as the consequences of the Urban Heat Island (UHI) (Li, Zhang, Liu, \& Huang, 2004; Philandras, Metaxas, \& Nastos, 1999), and deforestation effects \hyperref[b47]{(Voldoire \& Royer, 2004}) are capable of causing variation in the local patterns of temperatures. But in the case of Malaysia and perhaps major part of Southeast Asian sub-region studies have demonstrated continuous changes in the inter-annual rainfall variability and other anomalies \hyperref[b18]{(Juneng \& Tangang, 2005;}\hyperref[b42]{Tangang et al., 2007)}. Furthermore, \hyperref[b42]{Tangang et al. (2007)} averred that Malaysia being located amid the Indian and the Pacific Ocean which is the origin of the expanding Walker cell, it is probable that the IOD might have an effect on the country's temperature and rainfall variation. Hence the need for more studies on the rainfall behaviour especially at the micro levels to effectively comprehend the significance of the climatic changes has also been stressed \hyperref[b38]{(Suhaila et al., 2010)}.\par
This study analyzed and examined the trend in the changes associated with minimum and maximum annual, monthly and seasonal scale temperatures specific to the granary areas of MADA, Kedah, IADA, Barat Laut Selangor, and KADA, Kelantan in the Peninsular Malaysia. The study however employed Theil-Sen's Slope and measured the rate of changes of these variables temporally. 
\section[{II.}]{II.} 
\section[{Materials and Methods}]{Materials and Methods} 
\section[{a) Data Sources}]{a) Data Sources}\par
Time series data for mean monthly values of temperature for the period of  Generally, the need to identify outliers is one of the significant step in the data quality control (González-Rouco, Jiménez, Quesada, \& Valero, 2001). Outliers could be error in measurement or might be accurate extreme values. In the context of this dataset, the concern was not identification of erroneous observations especially with the regards to using hydrological parameters, but to reduce the size of the distribution tails \hyperref[b11]{(González-Rouco et al., 2001)}. In this study, possible outliers were examined using Q-Q plots of the individual data sets, the identified possible outliers were corrected by trimming the extreme values relative to the mean value. In this case outliers were considered as those values above a maximum threshold for each of the time series (González-Rouco et al   
\section[{d) Homogeneity Test}]{d) Homogeneity Test}\par
In this study the test of inhomogeneity of datasets (temperature series) was performed by applying four techniques of the standard normal homogeneity test, Buishand range test, Pettitt test, and Von Neumann ratio tests. The results showed that the data series obtained from the three stations were found to be homogeneous at ? = 0.05. 
\section[{e) Linear Regression Test}]{e) Linear Regression Test}\par
Simple linear regression is considered as a conventional approach and one of the simplest methods used in detecting changes in the time series of meteorological variables \hyperref[b1]{(Campra \& Morales 2016)}. Simple linear regression was applied to the temperature and rainfall data series. Trends and their 95\% CIs were estimated by least squares linear regression. Linear trends were estimated in every series from the slopes of the fit using values of monthly, annual and seasonal averages of Tmax, Tmin  {\ref (} Where X 's are the K explanatory variables (years) and y is the dependent variable (maximum temperature, minimum temperature). The slope line is b, and a is the intercept (value of y when x =0). The slope of regression describes the trend whether positive or negative. Linear regression works with the assumption of normal distribution. In this instance, the null hypothesis is that the slope of the line is zero or there is no trend in the temperature data. The significance of the slope shows the probability value (p-value). Microsoft excel was therefore, employed in the plotting of the trend lines and the Xlstat was used in the determination of the statistical values of the linear regression analysis. The P-value from the regression analysis was tested at the significance level ? = 0.05. The R 2 value or the square of the correlation from regression is used to indicate the strength of association and relationship between the variables X and Y. It has ratio between 0 and 1.0. R 2 value of 1.0 indicates stronger correlation and it means all points lie linearly. Whereas when R 2 is 0.0 it means no correlation or linear relation between variables X and Y. 
\section[{f) Mann Kendall Trend Test}]{f) Mann Kendall Trend Test}\par
In this study Mann Kendall test was used to examine the performance of a class of non -parametric trend test and the relative magnitude of the data rather than their measured values \hyperref[b17]{(Juahir et al., 2010;}\hyperref[b21]{Kendall 1975)}. In this context the method was used to detect long term trend of the meteorological variables (i.e. temperature) in the respective study areas.\par
Monthly and annual series were determined for each of the station using the seasonal Mann Kendall Trend Test \hyperref[b17]{(Juahir et al., 2010;}\hyperref[b30]{Río et al., 2011)}. XLSTAT software was used in the graphical presentation of the data sets. Moreover, XLSTAT and MAKESEN software were also used to calculate the statistical significance and estimation of trend using Sen Slope estimator for the variability and trend detection \hyperref[b30]{(Río et al., 2011)}.\par
The underlying principle of this model was based on the statistic (S) which is considered to be zero (0) meaning there was no trend. Each pair of observed values y i, y j (1> j) of the random variable was examined to find out whether y i > y j or y i < y j . The test statistic for the Mann-Kendall test was given as;S = ? ? ?????? ??? ?? ? ?? ?? ? ?? ?? =?? =1 ???1 ???1 (2.3) Where; Sign??? ?? ? ?? ?? ? = ? 1 ???? ?? ?? ? ?? ?? = 0 0 ?????? ?? ? ?? ?? = 0 ?1 ?????? ?? ? ?? ?? = 0 ? (2.4)\par
This means that the number of positive differences minus the number of negative differences. Variance of s is therefore computed by;Var (s) =[??(?? ? 1)(2?? + 5) ? ? (?? ? 1) ?? (2?? + 5)]/18 (2.5)\par
In a situation where n is greater than 10, the standard normal variate z is computed by using the following equation \hyperref[b8]{(Douglas, Vogel, \& Kroll, 2000)}.Z = ? ? ? ? ? ???1 ??????? (??) ???? ?? > 0 0 ???? ?? = 0 ??+1 ??????? (??) ???? ?? < 0 ? ? ? ? ? (2.6)\par
Therefore, the presence of a statistical trend is assessed using z value. A positive or negative z value is an indication of upward or downward trend. This study used Sen Slope estimator which is a non-parametric test procedure discovered by Sen \hyperref[b33]{(Sen 1968}) and advanced by \hyperref[b10]{(Gilbert, 1987)} to measure the actual slope in Mann -Kendall trend tests. It estimates the degree of any significant monotonic increase or decrease in trends examined in the Mann -Kendall S tests. The estimator Sen Slope was used where trend identified in the time series data is considered to be linear, illustrating the measure of the change per unit time \hyperref[b25]{(Mustapha, 2013;}\hyperref[b40]{Tabari \& Talaee, 2011)}. The Sen Slope estimator method is not sensitive to single data error or outliers. It is represented in the following equation;Q = ???? ?? ? ???? ?? (2.7)\par
Where Q represents the value of Sen in slope estimator; xj and xk are data values at time j and k. If there are n values of xj in the time series, the Sen Slope estimator is the median of n(n-1)/2 pairwise slopes, hence the Sen Slope estimator can be determined using;Q = ??? ??+1 2 ????? ?????? ?????? (2.8) Q = 1 2 (?? 2 ?? ) + ??? ??+2 2 ? ???? ?? ???? ???????? (2.9)\par
III. 
\section[{Results \& Discussions a) Descriptive Statistical Analysis}]{Results \& Discussions a) Descriptive Statistical Analysis}\par
Table \hyperref[tab_2]{1} provides a simple descriptive statistics of the annual values for the three variables used in this  From the mean temperature records (Table \hyperref[tab_2]{1}), it is evident that Subang Jaya station recorded the highest mean maximum temperature of 34.5 o C, followed by Alor Setar (34.4 o C) and Kota Bharu station observed the lowest mean maximum temperature (32.9 o C). For the Mean minimum temperature values, Subang Jaya observed the highest mean followed by Alor Setar and Kota Bharu recorded the lowest mean. Alor Setar and Subang Jaya observed average maximum temperature slightly higher than the Malaysia average (33 o C), while Kota Bharu station exhibited average maximum temperature almost the same with that of Malaysia. Similarly, all the three stations observed slightly higher average minimum temperature than the mean minimum temperature value of 22.0 o C for Malaysia (Chee-Wan \& Meng-Chang, 2012).\par
The values of the standard deviation further revealed the absolute variability of the temperatures. Subang Jaya station observed the highest absolute variability as indicated by the coefficients of standard deviation of 0.45 and 0.66 for both minimum and maximum temperatures respectively. Alor Setar recorded second highest standard deviation of 0.39 for the maximum temperature. Kota Bharu station observed the lowest standard deviation of 0.32 for the maximum temperature, but observed higher standard deviation of minimum temperature next to Subang Jaya. The coefficient of variability (CV) relatively indicated that Subang Jaya had the highest coefficient of variation of 0.013 in annual maximum temperature, Alor Setar observed second highest coefficient of variability (0.011) for maximum temperature. While Subang Jaya and Kota Bharu observed highest minimum temperature coefficient of variation of 0.029 each.\par
Table \hyperref[tab_3]{2} indicates the pattern of the temperature variability when the temperature variables were constructed into 10 year intervals to examine their possible changes in the mean, standard deviation and the coefficient of variability over time. In Alor Setar station the changes in mean annual maximum temperature and the mean annual minimum temperature were first steady, and between the periods 2001 to 2010 the mean annual maximum temperature and the mean annual minimum temperature appreciated with about 1.5\% and 0.9\% respectively with the decreasing variability.\par
For the Subang Jaya station, the mean annual maximum temperature decreased with about 6\% within the 30 years period, while the mean annual minimum temperature increased by about 2.8\% within the periods, but there were general increased in the temperature variability. At Kota Bharu station both the mean annual maximum and mean annual minimum temperatures increased by 0.3\% and 2.8\% respectively, with 1.1\% relative variability.    The highest mean monthly minimum temperature also coincided with the month of February (Figure \hyperref[fig_6]{4}). In Kota Bharu station (Figure  {\ref 6}), there is little variation in the mean monthly maximum temperature throughout the year, yet the month of May was the warmest (34.4 o C) followed by the month of June (34.3 o C) and the lowest mean maximum temperature corresponded with the month of December (30.6 o C). Also, average minimum temperature in this station indicates little variability (Figure  {\ref 6}). Characteristically the average minimum values peaked correspondingly along with the months of higher average monthly maximum temperature, which is April with 23 o C; but, the months of January through February were noted to be less warm (Figure  {\ref 7}). For the purpose of this study, the seasonal variability in terms of temperature were analysed based on the two paddy growing seasons, i.e. the main season is considered as the period when paddy is grown without supply of water from irrigation . Though varies usually from August/September to February/Mach the following year. The off season is regarded as the dry period when paddy planting normally depends on an irrigation system, the time mostly span from February/March until July/August (DOA, 2014). In this respect, Figures \hyperref[fig_9]{7-8} presents the descriptive statistics of the seasonal temperatures for the period 1981 to 3014 in the study areas.\par
From Figure \hyperref[fig_9]{8} the mean temperature for the main season was generally lower than the off season temperatures, probably due to the moderating effect of the relatively high amount of rainfall received during the main season. Comparatively, MADA and IADA recorded high mean maximum temperature of 34.2 o C each, while KADA had the lowest mean maximum value (32.3 o C).\par
During the off season (Figure \hyperref[fig_10]{9}),Subang Jaya station representing IADA recorded the highest mean maximum temperature (34.6 o C) followed by Alor Setar representing MADA (34.4 o C) and Kota Bharu station representing KADA also had the lowest mean. There was higher maximum temperature variability in IADA during the main season, in the MADA area showing less variability (Figure \hyperref[fig_9]{8}). But, in Figure \hyperref[fig_10]{9}, during the off-season the maximum temperature is more variable at IADA and KADA than at MADA. The minimum temperature shows greater variability more than the maximum temperature for all the seasons in the three areas. Comparatively, the main season minimum temperature for IADA is more variable than in MADA and KADA (Figure \hyperref[fig_10]{9}). But, minimum temperature for KADA was the least variable during the off-season compared to MADA and the highest variability of minimum temperature was observed in IADA during this season (Figure \hyperref[fig_10]{9}).  
\section[{b) Ordinary Linear Regression Trend Analysis}]{b) Ordinary Linear Regression Trend Analysis}\par
The result of the linear regression of the annual, monthly and seasonal maximum temperature, minimum temperature and rainfall trend analysis for the three study areas are calculated for the34 years independently. 
\section[{c) Maximum Temperature Regression Analysis}]{c) Maximum Temperature Regression Analysis}\par
The maximum temperature linear regression analysis for Alor Setar was conducted on the annual and monthly values, in Alor Setar, the linear trend line of the annual mean maximum temperature shows increasing trend, although there was a weak relationship between the maximum temperature changes and year (R 2 = .002) as illustrated in Figure \hyperref[fig_11]{10a}. Higher mean maximum temp is recorded in the months of February, August and October (Figure \hyperref[fig_11]{10b}). As for the monthly maximum temperature trend, the months of March, April, May, July and November showed decreasing trend, but the other months demonstrated an increasing trend. From the probability value (p-value) of the regression analyses, the coefficients of the monthly and annual trend lines were greater than 95\% confidence level (?=0.05), the null hypothesis (that there is no significant trend in the annual and monthly maximum temperature is therefore retained). This means that the annual and monthly maximum temperature exhibited no statistically significant trend, except for the months of January and February. The regression result also shows that more warming was recorded in the month of February with coefficient of .045 0 C with a unit change in the year.\par
Figure \hyperref[fig_1]{11a} shows the annual maximum temperature increasing linear trend with a moderate relationship between maximum temperature changes and the year (R 2 = .365), while Figure \hyperref[fig_1]{11b} illustrates the monthly maximum temperature pattern for Subang Jaya over the years. The maximum temperature linear regression analysis was further conducted on the annual and monthly values. The result showed that both annual and monthly maximum temperature demonstrated increasing trend. Only the months of January, February and March revealed statistically significant upward trend. There was no statistically significant increasing trend for the rest of months at 95\% confidence level. The maximum temperature for the month of March revealed more significant warming trend with the coefficient of .046 0 C for every unit increase in the year (p = 0.000). Figure \hyperref[fig_13]{13} shows the seasonal mean maximum temperature trends for the three study areas. The off seasons maximum temperature for Alor Setar and Kota Bharu showed decreasing trend, while the maximum temperature for all the seasons in Subang Jaya and main season in Alor Star and Kota Bharu showed increasing trend. All the R 2 shows a weak relationship between seasonal maximum temperature changes and the changes in the year.\par
The result for the linear regression indicated that only main season maximum temperature at Subang Jaya representing IADA showed statistically significant upward trend for the period with a coefficient of .016 (p = .002). Also the off seasons, mean maximum temperatures for Alor Setar representing MADA and Kota Bharu representing KADA revealed statistically significant downward trend with the coefficient of -.031 (p = .004)and -.022( p= .033) respectively.  
\section[{) Minimum Temperature Regression Analysis}]{) Minimum Temperature Regression Analysis}\par
Figure \hyperref[fig_15]{14a} shows the linear trend of the annual minimum temperatures while Figure \hyperref[fig_15]{14b} shows the mean monthly minimum temperature pattern. The result from the regression analysis on the mean annual minimum temperature and mean monthly minimum temperature for Alor Setar showed increasing trends for all the months as well as the mean annual minimum temperature. But only the minimum temperature for the months of February and July were shown to have no statistically significant trend at 95\% confidence level (? = .050). From the coefficients of the regression analyses, the minimum temperature for the month of March recorded highest increasing trend within the period with the coefficient of .054 0 C (p =.000). The lowest increasing trend was recorded in the month of  
\section[{Kota Bharu}]{Kota Bharu}\par
Subang Jaya Alor Setar October with coefficient of .014 0 C (p =.000). The monthly increasing minimum temperature trend for the month of March surpasses the annual increasing trend over the years with coefficient of .018 0 C (p = .000). R 2 shows a weak relationship between annual minimum temperature changes and the year.  \hyperref[fig_7]{15} a shows the annual linear trend for the mean minimum temperature for Subang Jaya with moderate relationship between annual minimum temperature changes and the year (R 2 = .328), while Figure \hyperref[fig_7]{15b} shows the pattern of the mean monthly temperatures. From the linear regression analysis, the result shows statistically significant upward trend for all the monthly and annual mean monthly minimum temperature for this station. The linear trend coefficient for the mean monthly minimum temperature ranges from .037 0 C to .063 0 C. The highest increasing trend was recorded in the month of March, while the months of June and July were the lowest increasing trend. Figure \hyperref[fig_17]{17} shows linear trends of seasonal mean minimum temperature trends for both main and off seasons for the three stations. The regression result for the seasonal changes in the mean minimum temperatures (Appendix F) showed statistically significant increasing trends for all the seasons. Subang Jaya recorded higher upward trend for all the season with the coefficient of .058 (p =.000) each. Alor Setar recorded lowest trend during the main season, while during the off-season Kota Bharu station recorded the lowest seasonal minimum temperature trend. All the R 2 values show that substantial changes in the seasonal temperatures were determined by the changes in the time (year). In Table \hyperref[tab_4]{3} present the results of Mann-Kendall and Sen's slope estimator for the mean annual and mean monthly maximum temperatures for Alor Setar. Based on this result the maximum annual temperature indicated a statistically not significant increasing trend in the data series. The MK test confirmed the result of the regression analysis which shows no statistically significant increase in the annual maximum temperature. Mann-Kendall trend revealed similar pattern of upward and downward trends. Similar to the linear regression analysis, the months of March, April and November have negative sign signifying decreasing trend in the monthly maximum temperature. While the months of May, July and October indicated no trend in their monthly maximum temperature, the rest of the months showed positive sign indicating increasing trend. From the result only the mean maximum temperature for the months of January, August and September shows statistically significant increasing trend. Highest warming was indicated during the months of January (Z= 2.69; Q = 0.05). Table \hyperref[tab_5]{4} presents the result from the Mann-Kendall and Sen's slope estimator for the annual and monthly maximum temperature for Subang Jaya station. All the monthly and annual maximum temperature for this station revealed increasing trend with the exception of the months of October and November exhibiting no trend and decreasing trend in the maximum temperature respectively. Only the months of February, March, April, September, and the annual maximum temperature for this station revealed statistical significant increasing trends. The MK trend test result confirmed the regression analysis where both the two tests revealed increasing trends in all the monthly maximum temperature except for the October maximum temperature. Table \hyperref[tab_6]{5} shows the result of the Mann-Kendall trend test and Sen's slope estimator for the annual and monthly maximum temperature for Kota Bharu station. The result shows a downward trend in the months of January, February, September and November maximum temperature, while the maximum temperature for the month of June showed no trend. All other months including the annual maximum temperature showed increasing trend. There were no statistically significant trends for all the months except for July, December and the annual maximum temperature at 95\% confidence level (? =.05). Table \hyperref[tab_7]{6} presents the result of Mann-Kendall trend test and Sen's slope estimator for the seasonal maximum temperature for the three stations. The result shows increasing trend in the maximum temperature for all the seasons in all the three study areas. Moreover, only the off-season temperature in all the three areas showed statistically significant increasing trend. More warming is observed during the off-season in Alor Setar (Z= 2.38; Q= .032) higher than all the other areas. The MK test result is similar to the regression analysis for seasonal maximum temperatures.  \hyperref[tab_8]{7} revealed statistically significant upward trends in minimum temperature records for all the months as well as for the annual values over the periods in Alor Setar, except for the months of February and July. The rising minimum temperature was observed to be significantly higher in the month of March (Z =3.73; Q = 0.056), whereas lowest significant increasing trend in the monthly minimum temperature was noticed in the month of October (Z= 2.15; Q=.0014). The MK trend statistics confirmed the regression result of the minimum temperature for this station. Table \hyperref[tab_9]{8} shows the result of Mann-Kendall and Sen's slope estimator for the minimum temperature in Subang Jaya stations over the years. The result revealed statistically significant increase in the minimum temperature for all the months as well as the annual minimum temperatures with the exception of the month of July which showed no trend at 95\% confidence level. Statistically significant upward warming corresponded with the month of March, December and the annual value with (.063 o C) each. The month of August recorded the lowest significant minimum temperature warming (.0.028 o C). The result confirmed the regression result for the minimum temperature for this station. Table \hyperref[tab_10]{9} shows the result of Mann-Kendall tests and Sen's slope estimator for monthly and the annual minimum temperature of Kota Bharu station. The result revealed increasing trend in the monthly as well the annual minimum temperature over the periods. But, the minimum temperature for the months of May and November were statistically not significant, while all the rest of the months as well as the annual values showed statistical significant increasing trends at 95\% confidence level.\par
For all the monthly temperatures, January minimum temperature observed the highest upward trends (Z= 3.21; Q= .056), while the month of April recorded the lowest upward trend (Z = 2.15; Q= .018). The result of Mann-Kendall and Sen's slope estimator of seasonal minimum temperature for the three stations are presented in Table \hyperref[tab_11]{10}. The result revealed statistically significant increasing trend for the minimum temperature in all the seasons. The rising seasonal minimum temperature was higher for all the seasons in Subang Jaya station, while Kota Bharu station observed the lowest upward trend for the two seasons.  
\section[{Summary}]{Summary}\par
In summary, the variability of temperature in the three study areas were investigated using descriptive statistics, parametric (least square regression) and nonparametric (Mann-Kendall and Sen's slope estimator). The study identified significant warming trend in the annual mean maximum temperature in two of the study areas, i.e. Subang Jaya and Kota Bharu typifying the climate over Peninsular Malaysia. Also significant warming trend detect in the annual minimum temperature and significant increasing trend in some of the monthly maximum and minimum temperatures for all the three stations. Also, the result reveals spatial and temporal variation in both the maximum and minimum temperature at annual, monthly and seasonal scales.\par
For the annual scale maximum temperature, this study identified a warming trend for the two stations with about 0. V. 
\section[{Conclusion}]{Conclusion}\par
Firstly, it can be concluded that there were spatial as well as temporal variation of temperature across the three granary areas. Secondly, the study identified significant warming trend in the annual mean maximum temperature in two of the study areas, i.e. Subang Jaya and Kota Bharu. Also significant warming trend were detected in the annual and seasonal minimum temperature for all the stations as well as significant increasing trend in some of the monthly maximum and minimum temperatures for all the three stations.\par
The findings from this study have the following implications; firstly, the result from the trend analysis provides an insight to agricultural development agencies in these areas and paddy farmers themselves, to make pre-emptive measure in relation to climate change variability. Timely measures and institutional actions will surely assist in ameliorating the damages that may be caused by the climate variability. This is in view of the fact that 34 years temperatures and rainfall data do not suffices the denial of the occurrence of climate change variability.\par
Funding: This research did not receive any specific grant from funding agencies in the public, commercial, or notfor-profit sectors. 
\section[{Conflict of Interests:}]{Conflict of Interests:}\par
The authors wish to summarily declare that they have no conflict of interest concerning the publication of this paper.\begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-2.png}
\caption{\label{fig_0}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1}\includegraphics[]{image-3.png}
\caption{\label{fig_1}Figure 1 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-4.png}
\caption{\label{fig_2}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-5.png}
\caption{\label{fig_3}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{2}\includegraphics[]{image-6.png}
\caption{\label{fig_4}Figure 2}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{233}\includegraphics[]{image-7.png}
\caption{\label{fig_5}Figure 2 :Figure 3 :Figure 3}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{4}\includegraphics[]{image-8.png}
\caption{\label{fig_6}Figure 4 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{5}\includegraphics[]{image-9.png}
\caption{\label{fig_7}Figure 5 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{67}\includegraphics[]{image-10.png}
\caption{\label{fig_8}Figure 6 :Figure 7 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{8}\includegraphics[]{image-11.png}
\caption{\label{fig_9}Figure 8 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{9}\includegraphics[]{image-12.png}
\caption{\label{fig_10}Figure 9 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{10}\includegraphics[]{image-13.png}
\caption{\label{fig_11}Figure 10 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1112}\includegraphics[]{image-14.png}
\caption{\label{fig_12}Figure 11 :Figure 12 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{13}\includegraphics[]{image-15.png}
\caption{\label{fig_13}Figure 13 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-16.png}
\caption{\label{fig_14}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{14}\includegraphics[]{image-17.png}
\caption{\label{fig_15}Figure 14 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1516}\includegraphics[]{image-18.png}
\caption{\label{fig_16}Figure 15 :Figure 16 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{17}\includegraphics[]{image-19.png}
\caption{\label{fig_17}Figure 17 :}\end{figure}
  \begin{figure}[htbp]
\noindent\textbf{1} \par 
\begin{longtable}{P{0.07858880778588807\textwidth}P{0.21094890510948905\textwidth}P{0.1902676399026764\textwidth}P{0.09306569343065693\textwidth}P{0.0951338199513382\textwidth}P{0.08479318734793187\textwidth}P{0.09720194647201946\textwidth}}
Granary\tabcellsep Variables\tabcellsep Min\tabcellsep Max\tabcellsep Mean\tabcellsep SD\tabcellsep CV\\
Alor Setar\tabcellsep Tmax ( 0 C)\tabcellsep 34.4\tabcellsep 35.0\tabcellsep 34.4\tabcellsep 0.39\tabcellsep 0.011\\
\tabcellsep Tmin( 0 C)\tabcellsep 22.2\tabcellsep 23.2\tabcellsep 22.5\tabcellsep 0.25\tabcellsep 0.011\\
\tabcellsep R/Fall(mm)\tabcellsep 1575.2\tabcellsep 2626.4\tabcellsep 2016.7\tabcellsep 249.3\tabcellsep 0.123\\
Subang Jaya\tabcellsep Tmax ( 0 C)\tabcellsep 34.5\tabcellsep 35.1\tabcellsep 34.5\tabcellsep 0.45\tabcellsep 0.013\\
\tabcellsep Tmin( 0 C)\tabcellsep 22.3\tabcellsep 23.8\tabcellsep 22.4\tabcellsep 0.66\tabcellsep 0.029\\
\tabcellsep R/Fall(mm)\tabcellsep 1944.8\tabcellsep 3210.3\tabcellsep 2551.8\tabcellsep 328.3\tabcellsep 0.017\\
Kota Bahru\tabcellsep Tmax ( 0 C)\tabcellsep 32.9\tabcellsep 33.6\tabcellsep 32.9\tabcellsep 0.32\tabcellsep 0.010\\
\tabcellsep Tmin( 0 C)\tabcellsep 22.1\tabcellsep 24.2\tabcellsep 22.3\tabcellsep 0.65\tabcellsep 0.029\\
\tabcellsep R/Fall(mm)\tabcellsep 1540.5\tabcellsep 3734.5\tabcellsep 2576.5\tabcellsep 595.8\tabcellsep 0.231\\
\tabcellsep \tabcellsep \multicolumn{5}{l}{Data source: Malaysia Meteorological Department}\end{longtable} \par
 
\caption{\label{tab_2}Table 1 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{2} \par 
\begin{longtable}{P{0.03289205702647657\textwidth}P{0.0675152749490835\textwidth}P{0.07097759674134418\textwidth}P{0.0804989816700611\textwidth}P{0.08482688391038697\textwidth}P{0.07184317718940936\textwidth}P{0.12118126272912422\textwidth}P{0.08482688391038697\textwidth}P{0.07184317718940936\textwidth}P{0.0804989816700611\textwidth}P{0.08309572301425662\textwidth}}
\tabcellsep Period\tabcellsep \tabcellsep 1981 -1990\tabcellsep \tabcellsep \tabcellsep 1991 -2000\tabcellsep \tabcellsep \tabcellsep 2001 -2010\tabcellsep \\
Station\tabcellsep Statistics\tabcellsep Tmax ( 0 C)\tabcellsep Tmin ( 0 C)\tabcellsep R/F (mm)\tabcellsep Tmax ( 0 C)\tabcellsep Tmin ( 0 C)\tabcellsep R/F (mm)\tabcellsep Tmax ( 0 C)\tabcellsep Tmin ( 0 C)\tabcellsep R/F (mm)\\
\tabcellsep Lowest\tabcellsep 33.9\tabcellsep 22.3\tabcellsep 1615.3\tabcellsep 33.9\tabcellsep 22.2\tabcellsep 1473.5\tabcellsep 34.5\tabcellsep 22.3\tabcellsep 1928\\
Alor Setar\tabcellsep Highest Mean SD\tabcellsep 34.8 34.2 0.327\tabcellsep 22.7 22.4 0.131\tabcellsep 2469.5 1989.8 284.2\tabcellsep 34.8 34.2 0.327\tabcellsep 22.6 22.4 0.141\tabcellsep 2183.9 1911.6 193.4\tabcellsep 35.0 34.7 0.162\tabcellsep 22.9 22.6 0.164\tabcellsep 2573.5 2248.8 193.5\\
\tabcellsep CV\tabcellsep 0.010\tabcellsep 0.006\tabcellsep 0.143\tabcellsep 0.010\tabcellsep 0.006\tabcellsep 0.101\tabcellsep 0.005\tabcellsep 0.007\tabcellsep 0.086\\
\tabcellsep Lowest\tabcellsep 34.0\tabcellsep 21.3\tabcellsep 1971.3\tabcellsep 34.3\tabcellsep 21.8\tabcellsep 2419.0\tabcellsep 34.1\tabcellsep 22.2\tabcellsep 2292.4\\
Subang Jaya\tabcellsep Highest Mean SD\tabcellsep 34.8 34.3 0.241\tabcellsep 22.1 21.6 0.287\tabcellsep 3331.4 2390.9 405.5\tabcellsep 34.8 34.5 0.171\tabcellsep 23.0 22.3 0.323\tabcellsep 2811.5 2646.5 126.1\tabcellsep 35.2 34.1 0.359\tabcellsep 23.3 22.6 0.294\tabcellsep 3455.2 2908.1 320.7\\
\tabcellsep CV\tabcellsep 0.07\tabcellsep 0.013\tabcellsep 0.170\tabcellsep 0.005\tabcellsep 0.014\tabcellsep 0.048\tabcellsep 0.010\tabcellsep 0.013\tabcellsep 0.110\\
\tabcellsep Min\tabcellsep 32.4\tabcellsep 21.7\tabcellsep 1540.5\tabcellsep 32.7\tabcellsep 21.9\tabcellsep 1689.0\tabcellsep 32.7\tabcellsep 22.0\tabcellsep 1928.6\\
Kota Bharu\tabcellsep Max Mean SD\tabcellsep 33.3 32.8 0.342\tabcellsep 22.2 21.8 0.147\tabcellsep 2859.3 2240.6 469.9\tabcellsep 33.6 33.1 0.273\tabcellsep 22.9 22.3 0.329\tabcellsep 3734.5 2886.2 671.2\tabcellsep 33.3 32.9 0.231\tabcellsep 22.8 22.4 0.238\tabcellsep 3566.2 2547.1 494.0\\
\tabcellsep CV\tabcellsep 0.010\tabcellsep 0.007\tabcellsep 0.210\tabcellsep 0.008\tabcellsep 0.015\tabcellsep 0.233\tabcellsep 0.007\tabcellsep 0.011\tabcellsep 0.194\\
\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \multicolumn{5}{l}{Data Source: Malaysia Meteorological Department}\end{longtable} \par
 
\caption{\label{tab_3}Table 2 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{3} \par 
\begin{longtable}{P{0.07541966426858514\textwidth}P{0.04280575539568345\textwidth}P{0.12026378896882495\textwidth}P{0.1528776978417266\textwidth}P{0.18752997601918464\textwidth}P{0.15491606714628298\textwidth}P{0.11618705035971223\textwidth}}
\tabcellsep \tabcellsep Alor Setar\tabcellsep \tabcellsep \tabcellsep Subang Jaya\tabcellsep \\
\tabcellsep \tabcellsep \tabcellsep Kota Bharu\tabcellsep \tabcellsep \tabcellsep \\
Month\tabcellsep n\tabcellsep Z-Value\tabcellsep Theil Sen's Slope (Q)\tabcellsep Trends\tabcellsep Kendall's Tau\tabcellsep P-Value\\
Jan\tabcellsep 34\tabcellsep 2.69\tabcellsep 0.05\tabcellsep Increasing\tabcellsep 0.329\tabcellsep 0.007\\
Feb\tabcellsep 34\tabcellsep 1.25\tabcellsep 0.025\tabcellsep No Trend\tabcellsep 0.155\tabcellsep 0.211\\
Mar\tabcellsep 34\tabcellsep -1.13\tabcellsep -0.014\tabcellsep No Trend\tabcellsep -0.139\tabcellsep 0.259\\
April\tabcellsep 34\tabcellsep -1.05\tabcellsep -0.02\tabcellsep No Trend\tabcellsep -0.130\tabcellsep 0.291\\
May\tabcellsep 34\tabcellsep 0.03\tabcellsep 0.0\tabcellsep No trend\tabcellsep 0.005\tabcellsep 0.976\\
Jun\tabcellsep 34\tabcellsep 0.78\tabcellsep 0.005\tabcellsep No Trend\tabcellsep 0.098\tabcellsep 0.438\\
Jul\tabcellsep 34\tabcellsep 0.22\tabcellsep 0.0\tabcellsep No trend\tabcellsep 0.029\tabcellsep 0.823\\
Aug\tabcellsep 34\tabcellsep 2.41\tabcellsep 0.029\tabcellsep Increasing\tabcellsep 0.297\tabcellsep 0.016\\
Sep\tabcellsep 34\tabcellsep 2.12\tabcellsep 0.033\tabcellsep Increasing\tabcellsep 0.260\tabcellsep 0.034\\
Oct\tabcellsep 34\tabcellsep 0.43\tabcellsep 0.0\tabcellsep No trend\tabcellsep 0.056\tabcellsep 0.665\end{longtable} \par
 
\caption{\label{tab_4}Table 3 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{4} \par 
\begin{longtable}{P{0.08397177419354838\textwidth}P{0.04627016129032258\textwidth}P{0.1028225806451613\textwidth}P{0.14909274193548389\textwidth}P{0.2090725806451613\textwidth}P{0.13538306451612903\textwidth}P{0.12338709677419356\textwidth}}
Month\tabcellsep n\tabcellsep Z-Value\tabcellsep Theil Sen's Slope (Q)\tabcellsep Trends\tabcellsep Kendall's Tau\tabcellsep P-Value\\
Jan\tabcellsep 34\tabcellsep 1.16\tabcellsep 0.025\tabcellsep No Trend\tabcellsep 0.198\tabcellsep 0.108\\
Feb\tabcellsep 34\tabcellsep 3.27\tabcellsep 0.033\tabcellsep Increasing\tabcellsep 0.405\tabcellsep 0.001\\
Mar\tabcellsep 34\tabcellsep 3.20\tabcellsep 0.050\tabcellsep Increasing\tabcellsep 0.392\tabcellsep 0.001\\
April\tabcellsep 34\tabcellsep 2.08\tabcellsep 0.029\tabcellsep Increasing\tabcellsep 0.255\tabcellsep 0.050\\
May\tabcellsep 34\tabcellsep 1.53\tabcellsep 0.024\tabcellsep No Trend\tabcellsep 0.190\tabcellsep 0.125\\
Jun\tabcellsep 34\tabcellsep 1.30\tabcellsep 0.014\tabcellsep No Trend\tabcellsep 0.161\tabcellsep 0.195\\
Jul\tabcellsep 34\tabcellsep 1.38\tabcellsep 0.020\tabcellsep No trend\tabcellsep 0.170\tabcellsep 0.167\\
Aug\tabcellsep 34\tabcellsep 1.32\tabcellsep 0.030\tabcellsep No Trend\tabcellsep 0.162\tabcellsep 0.186\\
Sep\tabcellsep 34\tabcellsep 2.44\tabcellsep 0.050\tabcellsep Increasing\tabcellsep 0.300\tabcellsep 0.015\\
Oct\tabcellsep 34\tabcellsep 2.48\tabcellsep 0.039\tabcellsep Increasing\tabcellsep 0.304\tabcellsep 0.013\\
Nov\tabcellsep 34\tabcellsep -0.85\tabcellsep -0.009\tabcellsep No Trend\tabcellsep -0.106\tabcellsep 0.396\\
Dec\tabcellsep 34\tabcellsep 1.40\tabcellsep 0.013\tabcellsep No Trend\tabcellsep 0.173\tabcellsep 0.162\\
Annual\tabcellsep 34\tabcellsep 2.18\tabcellsep 0.014\tabcellsep Increasing\tabcellsep 0.275\tabcellsep 0.029\end{longtable} \par
 
\caption{\label{tab_5}Table 4 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{5} \par 
\begin{longtable}{P{0.08397177419354838\textwidth}P{0.04627016129032258\textwidth}P{0.10796370967741935\textwidth}P{0.14909274193548389\textwidth}P{0.19879032258064513\textwidth}P{0.14052419354838708\textwidth}P{0.12338709677419356\textwidth}}
Month\tabcellsep n\tabcellsep Z-Value\tabcellsep Theil Sen's Slope (Q)\tabcellsep Trends\tabcellsep Kendall's Tau\tabcellsep P-Value\\
Jan\tabcellsep 34\tabcellsep -1.04\tabcellsep -0.025\tabcellsep No Trend\tabcellsep -0.129\tabcellsep 0.298\\
Feb\tabcellsep 34\tabcellsep -1.86\tabcellsep -0.03\tabcellsep No Trend\tabcellsep -0.229\tabcellsep 0.063\\
Mar\tabcellsep 34\tabcellsep 1.65\tabcellsep 0.021\tabcellsep No Trend\tabcellsep 0.203\tabcellsep 0.099\\
April\tabcellsep 34\tabcellsep 1.81\tabcellsep 0.025\tabcellsep No Trend\tabcellsep 0.223\tabcellsep 0.070\\
May\tabcellsep 34\tabcellsep 1.11\tabcellsep 0.007\tabcellsep No Trend\tabcellsep 0.139\tabcellsep 0.268\\
Jun\tabcellsep 34\tabcellsep 0.30\tabcellsep 0.0\tabcellsep No trend\tabcellsep 0.038\tabcellsep 0.766\\
Jul\tabcellsep 34\tabcellsep 2.34\tabcellsep 0.028\tabcellsep Increasing\tabcellsep 0.288\tabcellsep 0.020\\
Aug\tabcellsep 34\tabcellsep 1.25\tabcellsep 0.013\tabcellsep No Trend\tabcellsep 0.155\tabcellsep 0.211\\
Sep\tabcellsep 34\tabcellsep -1.01\tabcellsep -0.009\tabcellsep No Trend\tabcellsep -0.127\tabcellsep 0.311\\
Oct\tabcellsep 34\tabcellsep 1.26\tabcellsep 0.013\tabcellsep No Trend\tabcellsep 0.156\tabcellsep 0.206\\
Nov\tabcellsep 34\tabcellsep -1.41\tabcellsep -0.016\tabcellsep No Trend\tabcellsep -0.175\tabcellsep 0.158\\
Dec\tabcellsep 34\tabcellsep 2.69\tabcellsep 0.038\tabcellsep Increasing\tabcellsep 0.331\tabcellsep 0.007\\
Annual\tabcellsep 34\tabcellsep 2.16\tabcellsep 0.014\tabcellsep Increasing\tabcellsep 0.270\tabcellsep 0.031\end{longtable} \par
 
\caption{\label{tab_6}Table 5 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{6} \par 
\begin{longtable}{P{0.41666666666666663\textwidth}P{0.024074074074074074\textwidth}P{0.05740740740740741\textwidth}P{0.0925925925925926\textwidth}P{0.11111111111111112\textwidth}P{0.07962962962962962\textwidth}P{0.06851851851851852\textwidth}}
Max. Temperature\tabcellsep n\tabcellsep Z-value\tabcellsep Theil Sen's Slope(Q)\tabcellsep Trends\tabcellsep Kendall's Tau\tabcellsep P-value\\
Alor Setar\textunderscore MS\tabcellsep 34\tabcellsep 1.36\tabcellsep 0.008\tabcellsep No Trend\tabcellsep 0.170\tabcellsep 0.175\\
Alor Setar\textunderscore OS\tabcellsep 34\tabcellsep 2.38\tabcellsep 0.032\tabcellsep Increasing\tabcellsep 0.294\tabcellsep 0.017\\
Kota Bharu\textunderscore MS\tabcellsep 34\tabcellsep 0.69\tabcellsep 0.004\tabcellsep No Trend\tabcellsep 0.087\tabcellsep 0.492\\
Kota Bharu\textunderscore OS\tabcellsep 34\tabcellsep 2.11\tabcellsep 0.025\tabcellsep Increasing\tabcellsep 0.259\tabcellsep 0.035\\
Subang Jaya\textunderscore MS\tabcellsep 34\tabcellsep 3.02\tabcellsep 0.016\tabcellsep Increasing\tabcellsep 0.377\tabcellsep 0.003\\
Subang Jaya\textunderscore OS\tabcellsep 34\tabcellsep 1.65\tabcellsep 0.017\tabcellsep No Trend\tabcellsep 0.206\tabcellsep 0.098\\
\multicolumn{3}{l}{f) MK Test for Minimum Temperature}\tabcellsep \tabcellsep \tabcellsep \tabcellsep \\
Table\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \end{longtable} \par
 
\caption{\label{tab_7}Table 6 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{7} \par 
\begin{longtable}{P{0.08280318091451291\textwidth}P{0.045626242544731604\textwidth}P{0.09970178926441352\textwidth}P{0.1453280318091451\textwidth}P{0.22306163021868788\textwidth}P{0.13180914512922465\textwidth}P{0.12166998011928429\textwidth}}
Month\tabcellsep n\tabcellsep Z-Value\tabcellsep Theil Sen's Slope (Q)\tabcellsep Trends\tabcellsep Kendall's Tau\tabcellsep P-Value\\
Jan\tabcellsep 34\tabcellsep 2.56\tabcellsep 0.033\tabcellsep Increasing\tabcellsep 0.313\tabcellsep 0.011\\
Feb\tabcellsep 34\tabcellsep 1.31\tabcellsep 0.013\tabcellsep No Trend\tabcellsep 0.163\tabcellsep 0.190\\
Mar\tabcellsep 34\tabcellsep 3.73\tabcellsep 0.056\tabcellsep Increasing\tabcellsep 0.455\tabcellsep 0.000\\
April\tabcellsep 34\tabcellsep 3.36\tabcellsep 0.027\tabcellsep Increasing\tabcellsep 0.417\tabcellsep 0.001\\
May\tabcellsep 34\tabcellsep 4.29\tabcellsep 0.031\tabcellsep Increasing\tabcellsep 0.532\tabcellsep 0.000\\
Jun\tabcellsep 34\tabcellsep 2.13\tabcellsep 0.022\tabcellsep Increasing\tabcellsep 0.264\tabcellsep 0.033\\
Jul\tabcellsep 34\tabcellsep 1.71\tabcellsep 0.012\tabcellsep No Trend\tabcellsep 0.214\tabcellsep 0.088\\
Aug\tabcellsep 34\tabcellsep 2.67\tabcellsep 0.022\tabcellsep Increasing\tabcellsep 0.333\tabcellsep 0.008\\
Sep\tabcellsep 34\tabcellsep 3.02\tabcellsep 0.024\tabcellsep Increasing\tabcellsep 0.377\tabcellsep 0.003\\
Oct\tabcellsep 34\tabcellsep 2.15\tabcellsep 0.014\tabcellsep Increasing\tabcellsep 0.267\tabcellsep 0.032\\
Nov\tabcellsep 34\tabcellsep 3.64\tabcellsep 0.027\tabcellsep Increasing\tabcellsep 0.453\tabcellsep 0.000\\
Dec\tabcellsep 34\tabcellsep 3.44\tabcellsep 0.042\tabcellsep Increasing\tabcellsep 0.425\tabcellsep 0.001\\
Annual\tabcellsep 34\tabcellsep 4.16\tabcellsep 0.018\tabcellsep Increasing\tabcellsep 0.524\tabcellsep 0.000\end{longtable} \par
 
\caption{\label{tab_8}Table 7 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{8} \par 
\begin{longtable}{P{0.0821499013806706\textwidth}P{0.04526627218934911\textwidth}P{0.09891518737672583\textwidth}P{0.14418145956607495\textwidth}P{0.22800788954635107\textwidth}P{0.13076923076923078\textwidth}P{0.12071005917159763\textwidth}}
Month\tabcellsep n\tabcellsep Z-Value\tabcellsep Theil Sen's Slope (Q)\tabcellsep Trends\tabcellsep Kendall's Tau\tabcellsep P-Value\\
Jan\tabcellsep 34\tabcellsep 3.49\tabcellsep 0.052\tabcellsep Increasing\tabcellsep 0.425\tabcellsep 0.000\\
Feb\tabcellsep 34\tabcellsep 3.03\tabcellsep 0.040\tabcellsep Increasing\tabcellsep 0.375\tabcellsep 0.002\\
Mar\tabcellsep 34\tabcellsep 4.75\tabcellsep 0.063\tabcellsep Increasing\tabcellsep 0.584\tabcellsep 0.000\\
April\tabcellsep 34\tabcellsep 4.72\tabcellsep 0.057\tabcellsep Increasing\tabcellsep 0.577\tabcellsep 0.000\\
May\tabcellsep 34\tabcellsep 4.32\tabcellsep 0.047\tabcellsep Increasing\tabcellsep 0.530\tabcellsep 0.000\\
Jun\tabcellsep 34\tabcellsep 4.15\tabcellsep 0.050\tabcellsep Increasing\tabcellsep 0.510\tabcellsep 0.000\\
Jul\tabcellsep 34\tabcellsep 2.75\tabcellsep 0.041\tabcellsep Increasing\tabcellsep 0.337\tabcellsep 0.006\\
Aug\tabcellsep 34\tabcellsep 2.33\tabcellsep 0.028\tabcellsep Increasing\tabcellsep 0.288\tabcellsep 0.020\\
Sep\tabcellsep 34\tabcellsep 3.54\tabcellsep 0.033\tabcellsep Increasing\tabcellsep 0.437\tabcellsep 0.000\\
Oct\tabcellsep 34\tabcellsep 4.33\tabcellsep 0.040\tabcellsep Increasing\tabcellsep 0.532\tabcellsep 0.000\\
Nov\tabcellsep 34\tabcellsep 4.90\tabcellsep 0.062\tabcellsep Increasing\tabcellsep 0.599\tabcellsep 0.000\\
Dec\tabcellsep 34\tabcellsep 3.45\tabcellsep 0.063\tabcellsep Increasing\tabcellsep 0.422\tabcellsep 0.001\\
Annual\tabcellsep 34\tabcellsep 6.13\tabcellsep 0.063\tabcellsep Increasing\tabcellsep 0.751\tabcellsep 0.000\end{longtable} \par
 
\caption{\label{tab_9}Table 8 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{9} \par 
\begin{longtable}{P{0.08263888888888889\textwidth}P{0.045535714285714284\textwidth}P{0.13660714285714287\textwidth}P{0.10962301587301589\textwidth}P{0.22261904761904763\textwidth}P{0.13154761904761905\textwidth}P{0.12142857142857141\textwidth}}
Month\tabcellsep n\tabcellsep \multicolumn{2}{l}{Z-Value Theil Sen's Slope (Q)}\tabcellsep Trends\tabcellsep Kendall's Tau\tabcellsep P-Value\\
Jan\tabcellsep 34\tabcellsep 3.21\tabcellsep 0.056\tabcellsep Increasing\tabcellsep 0.390\tabcellsep 0.001\\
Feb\tabcellsep 34\tabcellsep 2.87\tabcellsep 0.050\tabcellsep Increasing\tabcellsep 0.351\tabcellsep 0.004\\
Mar\tabcellsep 34\tabcellsep 2.41\tabcellsep 0.044\tabcellsep Increasing\tabcellsep 0.295\tabcellsep 0.016\\
April\tabcellsep 34\tabcellsep 2.15\tabcellsep 0.018\tabcellsep Increasing\tabcellsep 0.268\tabcellsep 0.032\\
May\tabcellsep 34\tabcellsep 0.57\tabcellsep 0.004\tabcellsep No Trend\tabcellsep 0.072\tabcellsep 0.571\\
Jun\tabcellsep 34\tabcellsep 2.20\tabcellsep 0.023\tabcellsep Increasing\tabcellsep 0.272\tabcellsep 0.028\\
Jul\tabcellsep 34\tabcellsep 2.41\tabcellsep 0.025\tabcellsep Increasing\tabcellsep 0.299\tabcellsep 0.016\\
Aug\tabcellsep 34\tabcellsep 2.52\tabcellsep 0.024\tabcellsep Increasing\tabcellsep 0.312\tabcellsep 0.012\\
Sep\tabcellsep 34\tabcellsep 2.99\tabcellsep 0.027\tabcellsep Increasing\tabcellsep 0.371\tabcellsep 0.003\\
Oct\tabcellsep 34\tabcellsep 3.28\tabcellsep 0.019\tabcellsep Increasing\tabcellsep 0.409\tabcellsep 0.001\\
Nov\tabcellsep 34\tabcellsep 1.75\tabcellsep 0.010\tabcellsep No Trend\tabcellsep 0.221\tabcellsep 0.079\\
Dec\tabcellsep 34\tabcellsep 3.24\tabcellsep 0.041\tabcellsep Increasing\tabcellsep 0.399\tabcellsep 0.010\\
Annual\tabcellsep 34\tabcellsep 4.56\tabcellsep 0.028\tabcellsep Increasing\tabcellsep 0.569\tabcellsep 0.000\end{longtable} \par
 
\caption{\label{tab_10}Table 9 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{10} \par 
\begin{longtable}{P{0.37447552447552446\textwidth}P{0.025757575757575757\textwidth}P{0.06142191142191142\textwidth}P{0.09906759906759906\textwidth}P{0.13076923076923078\textwidth}P{0.0851981351981352\textwidth}P{0.0733100233100233\textwidth}}
Max. Temperature\tabcellsep n\tabcellsep Z-value\tabcellsep Theil Sen's Slope(Q)\tabcellsep Trends\tabcellsep Kendall's Tau\tabcellsep P-value\\
Alor Setar\textunderscore MS\tabcellsep 34\tabcellsep 3.76\tabcellsep 0.033\tabcellsep Increasing\tabcellsep 0.465\tabcellsep 0.000\\
Alor Setar\textunderscore OS\tabcellsep 34\tabcellsep 3.88\tabcellsep 0.030\tabcellsep Increasing\tabcellsep 0.483\tabcellsep 0.000\\
Kota Bharu\textunderscore MS\tabcellsep 34\tabcellsep 4.12\tabcellsep 0.032\tabcellsep Increasing\tabcellsep 0.515\tabcellsep 0.000\\
Kota Bharu\textunderscore OS\tabcellsep 34\tabcellsep 4.16\tabcellsep 0.023\tabcellsep Increasing\tabcellsep 0.519\tabcellsep 0.000\\
Subang Jaya\textunderscore MS\tabcellsep 34\tabcellsep 5.79\tabcellsep 0.059\tabcellsep Increasing\tabcellsep 0.710\tabcellsep 0.000\\
Subang Jaya\textunderscore OS\tabcellsep 34\tabcellsep 6.24\tabcellsep 0.060\tabcellsep Increasing\tabcellsep 0.764\tabcellsep 0.000\\
IV.\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \end{longtable} \par
 
\caption{\label{tab_11}Table 10 :}\end{figure}
 			\footnote{© 2020 Global Journals Determining Temperature trends in the Granary Areas of Peninsular Malaysia using Mann-Kendall and Sen's Slope Estimator} 		 		\backmatter   			 
\subsection[{Acknowledgement}]{Acknowledgement}\par
The authors wish to express their gratitude to Malaysia Meteorological Department (MMD) for providing the historical weather data used in this study. 			  			  				\begin{bibitemlist}{1}
\bibitem[Hartmann et al.]{b12}\label{b12} 	 		\textit{},  		 			D L Hartmann 		,  		 			A M K Tank 		,  		 			M Rusticucci 		,  		 			L V Alexander 		,  		 			S Brönnimann 		,  		 			Y A R Charabi 		.  		 	 
\bibitem[ Biometrika]{b35}\label{b35} 	 		\textit{},  	 	 		\textit{Biometrika}  		52  (3/4)  p. .  	 
\bibitem[ Geophysical Research Letters]{b45}\label{b45} 	 		\textit{},  	 	 		\textit{Geophysical Research Letters}  		42  (9)  p. .  	 
\bibitem[Amirabadizadeh et al. ()]{b0}\label{b0} 	 		\textit{},  		 			M Amirabadizadeh 		,  		 			Y F Huang 		,  		 			T Lee 		.  		 doi:10.1155/ 2015/ 579437.  	 	 		\textit{Malaysia. Advances in Meteorology}  		2015. 2015.  	 	 (Recent trends in temperature and precipitation in the Langat river basin) 
\bibitem[Rohde et al. ()]{b32}\label{b32} 	 		\textit{A new estimate of the average Earth surface land temperature spanning 1753 to},  		 			R Rohde 		,  		 			R Muller 		,  		 			R Jacobsen 		,  		 			E Muller 		,  		 			S Perlmutter 		,  		 			A Rosenfeld 		,  		 			. . Wickham 		,  		 			C 		.  		2013. 2011.  	 	 (Geoinfor Geostat Overview 1: 1. of, 7, 2) 
\bibitem[Wai et al. ()]{b48}\label{b48} 	 		\textit{A study of global warming in Malaysia},  		 			N M Wai 		,  		 			A Camerlengo 		,  		 			A Khairi 		,  		 			A Wahab 		.  		2005.  	 
\bibitem[Shapiro and Wilk ()]{b34}\label{b34} 	 		\textit{An analysis of variance test for normality},  		 			S S Shapiro 		,  		 			M B Wilk 		.  		1965.  	 	 (complete samples) 
\bibitem[Steinman et al. ()]{b36}\label{b36} 	 		‘Atlantic and Pacific multidecadal oscillations and Northern Hemisphere temperatures’.  		 			B A Steinman 		,  		 			M E Mann 		,  		 			S K Miller 		.  	 	 		\textit{Science}  		2015. 347  (6225)  p. .  	 
\bibitem[Tan et al. ()]{b41}\label{b41} 	 		‘Changes in precipitation extremes over the Kelantan River Basin’.  		 			M L Tan 		,  		 			A L Ibrahim 		,  		 			A P Cracknell 		,  		 			Z Yusop 		.  	 	 		\textit{Malaysia. International Journal of Climatology}  		2017. 37  (10)  p. .  	 
\bibitem[Chee-Wan and Meng-Chang ()]{b3}\label{b3} 	 		 			C Chee-Wan 		,  		 			C Meng-Chang 		.  		\textit{Country Report: Malaysia. Paper presented at the Asia Pacific Economic Cooperation (APEC) Workshop on food security},  				2012.  	 
\bibitem[Philandras et al. ()]{b27}\label{b27} 	 		‘Climate variability and urbanization in Athens’.  		 			C Philandras 		,  		 			D Metaxas 		,  		 			P T Nastos 		.  	 	 		\textit{Theoretical and Applied Climatology}  		1999. 63  (1-2)  p. .  	 
\bibitem[Mustapha ()]{b25}\label{b25} 	 		‘Detecting Surface Water Quality Trends Using Mann-Kendall Tests and Sen's Slope Estimates’.  		 			A Mustapha 		.  	 	 		\textit{International Journal of Advanced and Innovative Research}  		2013. IJAIR. p. .  	 
\bibitem[Sen ()]{b33}\label{b33} 	 		‘Estimates of the regression coefficient based on Kendall's tau’.  		 			P K Sen 		.  	 	 		\textit{Journal of the American Statistical Association}  		1968. 63  (324)  p. .  	 
\bibitem[Herath and Sarukkalige ()]{b14}\label{b14} 	 		‘Evaluation of empirical relationships between extreme rainfall and daily maximum temperature in Australia’.  		 			S M Herath 		,  		 			R Sarukkalige 		.  	 	 		\textit{Journal of Hydrology}  		2017.  	 
\bibitem[Juneng and Tangang ()]{b18}\label{b18} 	 		‘Evolution of ENSO-related rainfall anomalies in Southeast Asia region and its relationship with atmosphere-ocean variations in Indo-Pacific sector’.  		 			L Juneng 		,  		 			F T Tangang 		.  	 	 		\textit{Climate Dynamics}  		2005. 25  (4)  p. .  	 
\bibitem[Juahir et al. ()]{b17}\label{b17} 	 		‘Hydrological trend analysis due to land use changes at Langat River Basin’.  		 			H Juahir 		,  		 			S M Zain 		,  		 			A Z Aris 		,  		 			M K Yusof 		,  		 			M A A Samah 		,  		 			M Mokhtar 		.  	 	 		\textit{Environment Asia}  		2010. 2010. 3 p. .  	 
\bibitem[Ipcc ()]{b15}\label{b15} 	 		 			Ipcc 		.  		\textit{Climate Change: Synthesis Report. Retrieved from Valentia},  				 (Geneva)  		2007.  	 
\bibitem[Kendall ()]{b21}\label{b21} 	 		 			M G Kendall 		.  		\textit{Rank correlation methods},  				 (London)  		1975. Charles Griffin.  	 	 (4th ed.) 
\bibitem[Doane and Seward (201)]{b7}\label{b7} 	 		‘Measuring skewness: a forgotten statistic’.  		 			D P Doane 		,  		 			L E Seward 		.  	 	 		\textit{Journal of Statistics Education}  		201. 19  (2)  p. .  	 
\bibitem[Kaplan ()]{b13}\label{b13} 	 		\textit{Observations: atmosphere and surface Climate Change 2013 the Physical Science Basis: Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change},  		 			A Kaplan 		.  		2013. Cambridge University Press.  	 
\bibitem[Turco et al. ()]{b44}\label{b44} 	 		\textit{Observed climate change hotspots},  		 			M Turco 		,  		 			E Palazzi 		,  		 			J Hardenberg 		,  		 			A Provenzale 		.  		2015.  	 
\bibitem[Doa ()]{b5}\label{b5} 	 		\textit{Paddy Statistics of Malaysia},  		 			Doa 		.  		 BK 15/01-07/400.  		2014. 2013.  	 	 (First Edition ed.) 
\bibitem[Peninsular Malaysia: Department of Agriculture Peninsular Malaysia]{b6}\label{b6} 	 		\textit{Peninsular Malaysia: Department of Agriculture Peninsular Malaysia},  		 	 
\bibitem[Razali and Wah ()]{b28}\label{b28} 	 		‘Power comparisons of shapiro-wilk, kolmogorov-smirnov, lilliefors and anderson-darling tests’.  		 			N M Razali 		,  		 			Y B Wah 		.  	 	 		\textit{Journal of statistical modeling and analytics}  		2011. 2  (1)  p. .  	 
\bibitem[González-Rouco et al. ()]{b11}\label{b11} 	 		‘Quality control and homogeneity of precipitation data in the southwest of Europe’.  		 			J F González-Rouco 		,  		 			J L Jiménez 		,  		 			V Quesada 		,  		 			F Valero 		.  	 	 		\textit{Journal of Climate}  		2001. 14  (5)  p. .  	 
\bibitem[Wan et al. ()]{b49}\label{b49} 	 		‘Recent changes in extreme rainfall events in Peninsular Malaysia: 1971-2005’.  		 			Z Wan 		,  		 			W Z Suhaila 		,  		 			J Deni 		,  		 			S M Jemain 		,  		 			AA 		.  	 	 		\textit{Theoretical and Applied Climatology}  		2010. 99  (3-4)  p. .  	 
\bibitem[Zin et al. ()]{b50}\label{b50} 	 		‘Recent changes in extreme rainfall events in Peninsular Malaysia: 1971-2005’.  		 			W Z W Zin 		,  		 			S Jamaludin 		,  		 			S M Deni 		,  		 			A A Jemain 		.  	 	 		\textit{Theoretical and Applied Climatology}  		2010. 99  (3-4)  p. 303.  	 
\bibitem[Tabari and Talaee ()]{b40}\label{b40} 	 		‘Recent trends of mean maximum and minimum air temperatures in the western half of Iran’.  		 			H Tabari 		,  		 			P H Talaee 		.  	 	 		\textit{Meteorology and Atmospheric Physics}  		2011. 111  (3-4)  p. .  	 
\bibitem[Kaufmann et al. ()]{b20}\label{b20} 	 		‘Reconciling anthropogenic climate change with observed temperature’.  		 			R K Kaufmann 		,  		 			H Kauppi 		,  		 			M L Mann 		,  		 			J H Stock 		.  	 	 		\textit{Proceedings of the National Academy of Sciences}  		2011. 1998-2008. 108  (29)  p. .  	 
\bibitem[Lovejoy ()]{b23}\label{b23} 	 		‘Return periods of global climate fluctuations and the pause’.  		 			S Lovejoy 		.  	 	 		\textit{Geophysical Research Letters}  		2014. 41  (13)  p. .  	 
\bibitem[Muñoz-Díaz and Rodrigo ()]{b24}\label{b24} 	 		‘Seasonal rainfall variations in Spain (1912-2000) and their links to atmospheric circulation’.  		 			D Muñoz-Díaz 		,  		 			F Rodrigo 		.  	 	 		\textit{Atmospheric Research}  		2006. 81  (1)  p. .  	 
\bibitem[Varikoden et al. ()]{b46}\label{b46} 	 		‘Seasonal variation of rainfall characteristics in different intensity classes over Peninsular Malaysia’.  		 			H Varikoden 		,  		 			B Preethi 		,  		 			A Samah 		,  		 			C Babu 		.  	 	 		\textit{Journal of Hydrology}  		2011. 404  (1-2)  p. .  	 
\bibitem[Rodriguez-Puebla et al. ()]{b31}\label{b31} 	 		‘Spatial and temporal patterns of annual precipitationVariability over the Iberian Peninsula’.  		 			C Rodriguez-Puebla 		,  		 			A H Encinas 		,  		 			S Nieto 		,  		 			J Garmendia 		.  	 	 		\textit{International Journal of Climatology}  		1998. 18 p. .  	 
\bibitem[Río et al. ()]{b30}\label{b30} 	 		‘Spatial distribution of recent rainfall trends in Spain’.  		 			S D Río 		,  		 			L Herrero 		,  		 			R Fraile 		,  		 			A Penas 		.  	 	 		\textit{International Journal of Climatology}  		2011. 1961-2006. 31  (5)  p. .  	 
\bibitem[Nyatuame et al. ()]{b26}\label{b26} 	 		‘Statistical analysis of rainfall trend for volta region in Ghana’.  		 			M Nyatuame 		,  		 			V Owusu-Gyimah 		,  		 			F Ampiaw 		.  	 	 		\textit{International Journal of Atmospheric Sciences}  		2014. 2014.  	 
\bibitem[Gilbert ()]{b10}\label{b10} 	 		\textit{Statistical methods for environmental pollution monitoring},  		 			R O Gilbert 		.  		1987. John Wiley \& Sons.  	 
\bibitem[Stocker et al. ()]{b37}\label{b37} 	 		 			T Stocker 		,  		 			D Qin 		,  		 			G Plattner 		,  		 			M Tignor 		,  		 			S Allen 		,  		 			J Boschung 		,  		 			. . Midgley 		,  		 			P 		.  		\textit{Inter governmental panel on climate change, working group I contribution to the IPCC fifth assessment report (AR5)},  				 (New York)  		2013. Cambridge University Press.  	 
\bibitem[Trenberth and Paolino ()]{b43}\label{b43} 	 		‘The Northern Hemisphere sea-level pressure data set: Trends, errors and discontinuities’.  		 			K E Trenberth 		,  		 			D A PaolinoJr 		.  	 	 		\textit{Monthly Weather Review}  		1980. 108  (7)  p. .  	 
\bibitem[Eischeid et al. ()]{b9}\label{b9} 	 		‘The quality control of long-term climatological data using objective data analysis’.  		 			J K Eischeid 		,  		 			Bruce Baker 		,  		 			C Karl 		,  		 			T R Diaz 		,  		 			HF 		.  	 	 		\textit{Journal of Applied Meteorology}  		1995. 34  (12)  p. .  	 
\bibitem[Cramer and Howitt ()]{b4}\label{b4} 	 		\textit{The Sage dictionary of statistics: a practical resource for students in the social sciences},  		 			D Cramer 		,  		 			D L Howitt 		.  		2004. Sage.  	 
\bibitem[Suhaila and Yusop ()]{b39}\label{b39} 	 		‘Trend analysis and change point detection of annual and seasonal temperature series in Peninsular Malaysia’.  		 			J Suhaila 		,  		 			Z Yusop 		.  		 doi:10. 1007/s00703-017-0537-6.  	 	 		\textit{Meteorology and Atmospheric Physics}  		2017.  	 
\bibitem[Chakraborty et al. ()]{b2}\label{b2} 	 		\textit{Trend Analysis and Change Point Detection of Mean Air Temperature: A Spatio-Temporal Perspective of North-Eastern India. Environmental Processes},  		 			D Chakraborty 		,  		 			S Saha 		,  		 			R K Singh 		,  		 			B K Sethy 		,  		 			A Kumar 		,  		 			U S Saikia 		,  		 			. . Daschaudhuri 		,  		 			D 		.  		 \xref{http://dx.doi.org/10.1007/s40710-017-0263-6}{10.1007/s40710-017-0263-6}.  		2017.  	 
\bibitem[Campra and Morales ()]{b1}\label{b1} 	 		‘Trend analysis by a piecewise linear regression model applied to surface air temperatures in Southeastern Spain’.  		 			P Campra 		,  		 			M Morales 		.  		 \xref{http://dx.doi.org/10.5194/npg-2016-29}{10.5194/npg-2016-29}.  	 	 		\textit{Nonlinear Processes in Geophys. Discussions},  				2016. 1973-2014. 2016.  	 
\bibitem[Rimi et al. ()]{b29}\label{b29} 	 		‘Trend analysis of climate change and investigation on its probable impacts on rice production at Satkhira’.  		 			R H Rimi 		,  		 			S H Rahman 		,  		 			S Karmakar 		,  		 			S G Hussain 		.  	 	 		\textit{Bangladesh. Pak. J. Meteorol}  		2009. 6 p. .  	 
\bibitem[Jain and Kumar ()]{b16}\label{b16} 	 		‘Trend analysis of rainfall and temperature data for India’.  		 			S K Jain 		,  		 			V Kumar 		.  	 	 		\textit{Current Science}  		2012. p. .  	 
\bibitem[Tangang et al. ()]{b42}\label{b42} 	 		‘Trend and interannual variability of temperature in Malaysia’.  		 			F Tangang 		,  		 			L Juneng 		,  		 			S Ahmad 		.  	 	 		\textit{Theoretical and Applied Climatology}  		2007. 89  (3-4)  p. .  	 
\bibitem[Karmeshu ()]{b19}\label{b19} 	 		\textit{Trend detection in annual temperature \& precipitation using the Mann Kendall test-a case study to assess climate change on select states in the northeastern United States},  		 			N Karmeshu 		.  		2012.  	 
\bibitem[Douglas et al. ()]{b8}\label{b8} 	 		‘Trends in floods and low flows in the United States: impact of spatial correlation’.  		 			E Douglas 		,  		 			R Vogel 		,  		 			C Kroll 		.  	 	 		\textit{Journal of Hydrology}  		2000. 240  (1)  p. .  	 
\bibitem[Suhaila et al. ()]{b38}\label{b38} 	 		‘Trends in Peninsular Malaysia rainfall data during the southwest monsoon and northeast monsoon seasons’.  		 			J Suhaila 		,  		 			S M Deni 		,  		 			W Z W Zin 		,  		 			A Jemain 		.  	 	 		\textit{Sains Malaysiana}  		2010. 39  (4)  p. .  	 
\bibitem[Voldoire and Royer ()]{b47}\label{b47} 	 		‘Tropical deforestation and climate variability’.  		 			A Voldoire 		,  		 			J Royer 		.  	 	 		\textit{Climate Dynamics}  		2004. 22  (8)  p. .  	 
\bibitem[Li et al. ()]{b22}\label{b22} 	 		‘Urban heat island effect on annual mean temperature during the last 50 years in China’.  		 			Q Li 		,  		 			H Zhang 		,  		 			X Liu 		,  		 			J Huang 		.  	 	 		\textit{Theoretical and Applied Climatology}  		2004. 79  (3-4)  p. .  	 
\end{bibitemlist}
 			 		 	 
\end{document}
