<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AJCC</journal-id><journal-title-group><journal-title>American Journal of Climate Change</journal-title></journal-title-group><issn pub-type="epub">2167-9495</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ajcc.2014.35038</article-id><article-id pub-id-type="publisher-id">AJCC-52628</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  A Study of Decadal December Temperature Variability in Pakistan
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>omana</surname><given-names>Ambreen</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Iftikhar</surname><given-names>Ahmad</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Shahzad</surname><given-names>Sultan</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Zhaobo</surname><given-names>Sun</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Muhammad</surname><given-names>Nawaz</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff3"><addr-line>Nanjing University of Information Science and Technology, Nanjing, China</addr-line></aff><aff id="aff1"><addr-line>Department of Geography, University of Balochistan, Quetta, Pakistan</addr-line></aff><aff id="aff2"><addr-line>Pakistan Meteorological Department, Islamabad, Pakistan</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>iftigeog@gmail.com(IA)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>25</day><month>12</month><year>2014</year></pub-date><volume>03</volume><issue>05</issue><fpage>429</fpage><lpage>437</lpage><history><date date-type="received"><day>19</day>	<month>September</month>	<year>2014</year></date><date date-type="rev-recd"><day>21</day>	<month>October</month>	<year>2014</year>	</date><date date-type="accepted"><day>15</day>	<month>November</month>	<year>2014</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The spatial distribution of December temperature in Pakistan has been assessed by statistical method based on mean monthly data from 51 ground stations. The analysis is performed at decadal scale over the period of 1950-2000. December is one of the representative months of winter season in Pakistan, the country with subtropical location and complex rugged terrains, plateaus and plains. The results support a slight rising temperature trend in December. However, this change in temperature varies from region to region as well as from decade to decade and reflects a complicated spatial-temporal structure of temperature anomalies. The assessment shows that the temperature anomalies in different national territories at local scale do not follow the assumption that winter months are warming in northern hemisphere. Both the isothermal shift and temperature anomalies confirm that the mountainous areas of Pakistan face more temperature variability than plains.
 
</p></abstract><kwd-group><kwd>Temperature Anomaly</kwd><kwd> December</kwd><kwd> Pakistan</kwd><kwd> Decadal Scale</kwd><kwd> Spatial Distribution</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The true essence of climate change is hard to be understood without the temperature analysis at different temporal and spatial scales. Pakistan faces extreme weather events like drought, floods and heat waves [<xref ref-type="bibr" rid="scirp.52628-ref1">1</xref>] -[<xref ref-type="bibr" rid="scirp.52628-ref4">4</xref>] . The variability of temperature over complex rugged terrains remains great attraction for searchers in climate change sciences e.g. [<xref ref-type="bibr" rid="scirp.52628-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.52628-ref6">6</xref>] and surface air temperature is the most important climate element that helps understanding of climate variability [<xref ref-type="bibr" rid="scirp.52628-ref7">7</xref>] . In the current national scenario, economic development is not sustainable in the agro- based developing society of Pakistan without understanding of climate change and its impacts on agriculture [<xref ref-type="bibr" rid="scirp.52628-ref8">8</xref>] - [<xref ref-type="bibr" rid="scirp.52628-ref10">10</xref>] where the faster growing population already has faced meteorological hazards [<xref ref-type="bibr" rid="scirp.52628-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.52628-ref12">12</xref>] . The monthly diagnostic of temperature variability is instrumental in future climate policy [<xref ref-type="bibr" rid="scirp.52628-ref13">13</xref>] - [<xref ref-type="bibr" rid="scirp.52628-ref16">16</xref>] . More than half of the total area of Pakistan is dominated by rugged mountains including Hindukush Karakoram Himalayas (HKH). The Himalayas were already declared sensitive area in term of high glacial ablation by Working Group II in the Asian Chapter of IPCC [<xref ref-type="bibr" rid="scirp.52628-ref17">17</xref>] .</p><p>This paper is focused on the analysis of December temperature variability with emphasis on spatial distribution and its regional detail in Pakistan. Most of climatic studies are supported by “time series figures of temperature”, while rarer studies map the anomalies and also neglect the spatial distribution of the observations. In this paper, we mapped the isotherms through which we learned about the shifts of average temperature at regional scale and recognized their anomalies by simulating temperature coefficient at decadal scale.</p></sec><sec id="s2"><title>2. Data and Methods</title><p>The monthly temperature data of surface air temperature were obtained from Pakistan Meteorological Department (PMD). The 51 stations (<xref ref-type="fig" rid="fig1">Figure 1</xref>, <xref ref-type="table" rid="table1">Table 1</xref>) were selected from all over the country. The mean monthly data of December temperature were analyzed separately for each decade from 1950s to 2000s.</p><p>The isothermal maps (<xref ref-type="fig" rid="fig2">Figure 2</xref>) for each decade were prepared to check the shift of average temperatures. The distribution of spatial anomalies is shown with the help of maps (<xref ref-type="fig" rid="fig3">Figure 3</xref>) that configure the spatial distribution of temperature coefficient. This simulation is instrumental in understanding the temperature variability in different parts of the country.</p><p>The following technique was adopted to calculate the temperature coefficient trend,</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x5.png" xlink:type="simple"/></inline-formula>.</p><p>The above relation establishes linear regression between, the time series t<sub>i</sub> and variable x<sub>i</sub> (temperature) for the specified period.</p><p>By taking into account t<sub>i</sub> as independent and x<sub>i</sub> dependent variable, regression coefficient “b” and the regression constant “a” of least-squares estimation have been calculated, respectively by using the following principle.</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> The geographical distribution of the meteorological stations utilized in the study</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x6.png"/></fig><table-wrap-group id="1"><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> The 51 stations data were utilized in the analysis that almost covers the whole national territories</title></caption><table-wrap id="1_1"><table><tbody><thead><tr><th align="center" valign="middle" >Serial No.</th><th align="center" valign="middle" >Stations</th><th align="center" valign="middle" >Latitude</th><th align="center" valign="middle" >Longitude</th><th align="center" valign="middle" >Region</th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >Astor</td><td align="center" valign="middle" >35<sup>˚</sup>20'</td><td align="center" valign="middle" >74<sup>˚</sup>54'</td><td align="center" valign="middle" >Gilgit-Baltistan</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >Badin</td><td align="center" valign="middle" >24<sup>˚</sup>38'</td><td align="center" valign="middle" >68<sup>˚</sup>54'</td><td align="center" valign="middle" >Sindh</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >Bahawal Nagar</td><td align="center" valign="middle" >29<sup>˚</sup>57'</td><td align="center" valign="middle" >73<sup>˚</sup>51'</td><td align="center" valign="middle" >Punjab</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >Bahawal Pur</td><td align="center" valign="middle" >29<sup>˚</sup>24'</td><td align="center" valign="middle" >71<sup>˚</sup>47'</td><td align="center" valign="middle" >Punjab</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >Balakot</td><td align="center" valign="middle" >34<sup>˚</sup>23'</td><td align="center" valign="middle" >73<sup>˚</sup>21'</td><td align="center" valign="middle" >KPK</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >Bar Khan</td><td align="center" valign="middle" >29<sup>˚</sup>53'</td><td align="center" valign="middle" >69<sup>˚</sup>43'</td><td align="center" valign="middle" >Balcohistan</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >Bunji</td><td align="center" valign="middle" >35<sup>˚</sup>67'</td><td align="center" valign="middle" >74<sup>˚</sup>63'</td><td align="center" valign="middle" >Gilgit-Baltistan</td></tr><tr><td align="center" valign="middle" >8</td><td align="center" valign="middle" >Cherat</td><td align="center" valign="middle" >33<sup>˚</sup>49'</td><td align="center" valign="middle" >71<sup>˚</sup>53'</td><td align="center" valign="middle" >KPK</td></tr><tr><td align="center" valign="middle" >9</td><td align="center" valign="middle" >Chhor</td><td align="center" valign="middle" >25<sup>˚</sup>31'</td><td align="center" valign="middle" >69<sup>˚</sup>47'</td><td align="center" valign="middle" >Sindh</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle" >Chilas</td><td align="center" valign="middle" >35<sup>˚</sup>25'</td><td align="center" valign="middle" >74<sup>˚</sup>06'</td><td align="center" valign="middle" >Gilgit-Baltistan</td></tr><tr><td align="center" valign="middle" >11</td><td align="center" valign="middle" >Chitral</td><td align="center" valign="middle" >35<sup>˚</sup>51'</td><td align="center" valign="middle" >71<sup>˚</sup>50'</td><td align="center" valign="middle" >KPK</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >Dalbandin</td><td align="center" valign="middle" >28<sup>˚</sup>53'</td><td align="center" valign="middle" >64<sup>˚</sup>24'</td><td align="center" valign="middle" >Balochistan</td></tr><tr><td align="center" valign="middle" >13</td><td align="center" valign="middle" >Dera Ismail Khan</td><td align="center" valign="middle" >31<sup>˚</sup>55'</td><td align="center" valign="middle" >70<sup>˚</sup>52'</td><td align="center" valign="middle" >KPK</td></tr><tr><td align="center" valign="middle" >14</td><td align="center" valign="middle" >Dir</td><td align="center" valign="middle" >35<sup>˚</sup>12'</td><td align="center" valign="middle" >71<sup>˚</sup>51'</td><td align="center" valign="middle" >KPK</td></tr><tr><td align="center" valign="middle" >15</td><td align="center" valign="middle" >Drosh</td><td align="center" valign="middle" >35<sup>˚</sup>34'</td><td align="center" valign="middle" >71<sup>˚</sup>47'</td><td align="center" valign="middle" >KPK</td></tr><tr><td align="center" valign="middle" >16</td><td align="center" valign="middle" >Faisalabad</td><td align="center" valign="middle" >31<sup>˚</sup>26'</td><td align="center" valign="middle" >73<sup>˚</sup>08'</td><td align="center" valign="middle" >Punjab</td></tr><tr><td align="center" valign="middle" >17</td><td align="center" valign="middle" >Garhi Dupatta</td><td align="center" valign="middle" >34<sup>˚</sup>13'</td><td align="center" valign="middle" >73<sup>˚</sup>37</td><td align="center" valign="middle" >AJK</td></tr><tr><td align="center" valign="middle" >18</td><td align="center" valign="middle" >Gilgit</td><td align="center" valign="middle" >35<sup>˚</sup>55'</td><td align="center" valign="middle" >74<sup>˚</sup>20'</td><td align="center" valign="middle" >Gilgit-Baltistan</td></tr><tr><td align="center" valign="middle" >19</td><td align="center" valign="middle" >Hyderabad</td><td align="center" valign="middle" >25<sup>˚</sup>23'</td><td align="center" valign="middle" >68<sup>˚</sup>25'</td><td align="center" valign="middle" >Sindh</td></tr><tr><td align="center" valign="middle" >20</td><td align="center" valign="middle" >Islamabad</td><td align="center" valign="middle" >33<sup>˚</sup>43'</td><td align="center" valign="middle" >73<sup>˚</sup>05'</td><td align="center" valign="middle" >Punjab</td></tr><tr><td align="center" valign="middle" >21</td><td align="center" valign="middle" >Jacobabad</td><td align="center" valign="middle" >28<sup>˚</sup>18'</td><td align="center" valign="middle" >68<sup>˚</sup>28'</td><td align="center" valign="middle" >Sindh</td></tr><tr><td align="center" valign="middle" >22</td><td align="center" valign="middle" >Jhang</td><td align="center" valign="middle" >31<sup>˚</sup>27'</td><td align="center" valign="middle" >73˚32'</td><td align="center" valign="middle" >Punjab</td></tr><tr><td align="center" valign="middle" >23</td><td align="center" valign="middle" >Jiwani</td><td align="center" valign="middle" >25<sup>˚</sup>04'</td><td align="center" valign="middle" >61<sup>˚</sup>48'</td><td align="center" valign="middle" >Balochistan</td></tr><tr><td align="center" valign="middle" >24</td><td align="center" valign="middle" >Kakul</td><td align="center" valign="middle" >34<sup>˚</sup>11'</td><td align="center" valign="middle" >73<sup>˚</sup>15'</td><td align="center" valign="middle" >KPK</td></tr><tr><td align="center" valign="middle" >25</td><td align="center" valign="middle" >Kalat</td><td align="center" valign="middle" >29<sup>˚</sup>02'</td><td align="center" valign="middle" >66<sup>˚</sup>35'</td><td align="center" valign="middle" >Balochistan</td></tr><tr><td align="center" valign="middle" >26</td><td align="center" valign="middle" >Karachi</td><td align="center" valign="middle" >24<sup>˚</sup>54'</td><td align="center" valign="middle" >66<sup>˚</sup>56</td><td align="center" valign="middle" >Sindh</td></tr><tr><td align="center" valign="middle" >27</td><td align="center" valign="middle" >Khanpur</td><td align="center" valign="middle" >28<sup>˚</sup>39'</td><td align="center" valign="middle" >70<sup>˚</sup>41'</td><td align="center" valign="middle" >Punjab</td></tr><tr><td align="center" valign="middle" >28</td><td align="center" valign="middle" >Kohat</td><td align="center" valign="middle" >33<sup>˚</sup>57'</td><td align="center" valign="middle" >71<sup>˚</sup>43'</td><td align="center" valign="middle" >KPK</td></tr><tr><td align="center" valign="middle" >29</td><td align="center" valign="middle" >Kotli</td><td align="center" valign="middle" >33<sup>˚</sup>31'</td><td align="center" valign="middle" >73<sup>˚</sup>54'</td><td align="center" valign="middle" >Sindh</td></tr><tr><td align="center" valign="middle" >30</td><td align="center" valign="middle" >Lahore</td><td align="center" valign="middle" >3<sup> </sup>1<sup>˚</sup>35'</td><td align="center" valign="middle" >74<sup>˚</sup>24'</td><td align="center" valign="middle" >Punjab</td></tr><tr><td align="center" valign="middle" >31</td><td align="center" valign="middle" >Larkana</td><td align="center" valign="middle" >27<sup>˚</sup>32'</td><td align="center" valign="middle" >68<sup>˚</sup>14'</td><td align="center" valign="middle" >Sindh</td></tr><tr><td align="center" valign="middle" >32</td><td align="center" valign="middle" >Mianwali</td><td align="center" valign="middle" >32<sup>˚</sup>55'</td><td align="center" valign="middle" >71<sup>˚</sup>52'</td><td align="center" valign="middle" >Punjab</td></tr><tr><td align="center" valign="middle" >33</td><td align="center" valign="middle" >Multan</td><td align="center" valign="middle" >30<sup>˚</sup>12'</td><td align="center" valign="middle" >71<sup>˚</sup>26'</td><td align="center" valign="middle" >Punjab</td></tr><tr><td align="center" valign="middle" >34</td><td align="center" valign="middle" >Murree</td><td align="center" valign="middle" >33<sup>˚</sup>55'</td><td align="center" valign="middle" >73<sup>˚</sup>23'</td><td align="center" valign="middle" >Punjab</td></tr><tr><td align="center" valign="middle" >35</td><td align="center" valign="middle" >Muzaffarabad</td><td align="center" valign="middle" >34<sup>˚</sup>22'</td><td align="center" valign="middle" >73<sup>˚</sup>29'</td><td align="center" valign="middle" >AJK</td></tr><tr><td align="center" valign="middle" >36</td><td align="center" valign="middle" >Nawabshah</td><td align="center" valign="middle" >26<sup>˚</sup>15'</td><td align="center" valign="middle" >68<sup>˚</sup>22'</td><td align="center" valign="middle" >Sindh</td></tr><tr><td align="center" valign="middle" >37</td><td align="center" valign="middle" >Ormara</td><td align="center" valign="middle" >25<sup>˚</sup>12'</td><td align="center" valign="middle" >64<sup>˚</sup>40'</td><td align="center" valign="middle" >Balochistan</td></tr><tr><td align="center" valign="middle" >38</td><td align="center" valign="middle" >Padidan</td><td align="center" valign="middle" >26<sup>˚</sup>51'</td><td align="center" valign="middle" >68<sup>˚</sup>08'</td><td align="center" valign="middle" >Sindh</td></tr><tr><td align="center" valign="middle" >39</td><td align="center" valign="middle" >Panjgur</td><td align="center" valign="middle" >26<sup>˚</sup>58'</td><td align="center" valign="middle" >64<sup>˚</sup>06'</td><td align="center" valign="middle" >Balochistan</td></tr></tbody></table></table-wrap><table-wrap id="1_2"><table><tbody><thead><tr><th align="center" valign="middle" >40</th><th align="center" valign="middle" >Parachinar</th><th align="center" valign="middle" >33<sup>˚</sup>52'</th><th align="center" valign="middle" >70<sup>˚</sup>05'</th><th align="center" valign="middle" >FATA</th></tr></thead><tr><td align="center" valign="middle" >41</td><td align="center" valign="middle" >Pasni</td><td align="center" valign="middle" >25<sup>˚</sup>16'</td><td align="center" valign="middle" >63<sup>˚</sup>29'</td><td align="center" valign="middle" >Balochistan</td></tr><tr><td align="center" valign="middle" >42</td><td align="center" valign="middle" >Peshawar</td><td align="center" valign="middle" >34<sup>˚</sup>01</td><td align="center" valign="middle" >71<sup>˚</sup>35</td><td align="center" valign="middle" >KPK</td></tr><tr><td align="center" valign="middle" >43</td><td align="center" valign="middle" >Quetta</td><td align="center" valign="middle" >30<sup>˚</sup>05'</td><td align="center" valign="middle" >66<sup>˚</sup>57'</td><td align="center" valign="middle" >Balochistan</td></tr><tr><td align="center" valign="middle" >44</td><td align="center" valign="middle" >Risalpur</td><td align="center" valign="middle" >71<sup>˚</sup>98'</td><td align="center" valign="middle" >34<sup>˚</sup>07'</td><td align="center" valign="middle" >KPK</td></tr><tr><td align="center" valign="middle" >45</td><td align="center" valign="middle" >Rohri</td><td align="center" valign="middle" >27<sup>˚</sup>40'</td><td align="center" valign="middle" >68<sup>˚</sup>54'</td><td align="center" valign="middle" >Sindh</td></tr><tr><td align="center" valign="middle" >46</td><td align="center" valign="middle" >Saidu Sharif</td><td align="center" valign="middle" >34<sup>˚</sup>44'</td><td align="center" valign="middle" >72<sup>˚</sup>33'</td><td align="center" valign="middle" >KPK</td></tr><tr><td align="center" valign="middle" >47</td><td align="center" valign="middle" >Sarghoda</td><td align="center" valign="middle" >32<sup>˚</sup>05'</td><td align="center" valign="middle" >72<sup>˚</sup>67'</td><td align="center" valign="middle" >Punjab</td></tr><tr><td align="center" valign="middle" >48</td><td align="center" valign="middle" >Sialkot</td><td align="center" valign="middle" >32<sup>˚</sup>31'</td><td align="center" valign="middle" >74<sup>˚</sup>32'</td><td align="center" valign="middle" >Punjab</td></tr><tr><td align="center" valign="middle" >49</td><td align="center" valign="middle" >Sibbi</td><td align="center" valign="middle" >29<sup>˚</sup>33'</td><td align="center" valign="middle" >67<sup>˚</sup>53'</td><td align="center" valign="middle" >Balochistan</td></tr><tr><td align="center" valign="middle" >50</td><td align="center" valign="middle" >Skardu</td><td align="center" valign="middle" >35<sup>˚</sup>18'</td><td align="center" valign="middle" >75<sup>˚</sup>41'</td><td align="center" valign="middle" >Gilgit-Baltistan</td></tr><tr><td align="center" valign="middle" >51</td><td align="center" valign="middle" >Zhob</td><td align="center" valign="middle" >31<sup>˚</sup>21'</td><td align="center" valign="middle" >69<sup>˚</sup>28'</td><td align="center" valign="middle" >Balochistan</td></tr></tbody></table></table-wrap></table-wrap-group><disp-formula id="scirp.52628-formula1048"><graphic  xlink:href="http://html.scirp.org/file/2-2360206x7.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3"><title>3. Results and Discussion</title><p>Based on results <xref ref-type="fig" rid="fig2">Figure 2</xref> and <xref ref-type="fig" rid="fig3">Figure 3</xref> show the geography of average temperature marked by isotherms and observed decadal temperature anomalies marked by temperature coefficient respectively. The comparison of various decades configures that the isotherms were shifted from decade to decade. It means spatial temperature changes were there in the study area. Therefore, the temperature anomalies <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x8.png" xlink:type="simple"/></inline-formula> mark these changes by statistical significance.</p><p>Pakistan have variety of temperatures from south (north) to north (south) and east (west) to west (east) in December. The reason for this obvious difference in temperature in the same month is latitudinal extent and landforms. Obviously, the temperature divides Pakistan into parts of high and low temperatures in Indus Plains and north/western rugged parts respectively. The said reasons are (were) only true for the distribution of average temperatures in normal atmospheric condition but the distribution of spatial anomalies was quite complex and varied from area to area. The comparison of various decades is evident of the isothermal shift for a selected temperature (e.g. 15˚C or 12˚C or any other isotherm), it means temperature change occurs at local scale in the study period. It has been observed that where the isotherms were closely spaced the temperature anomalies were more obvious and vice versa. This mostly happened in mountains and confirmed that rugged lofty portions of the country were more susceptible to temperature variability.</p><sec id="s3_1"><title>3.1. The Decade of 1950s</title><p>Reference to <xref ref-type="fig" rid="fig2">Figure 2</xref>(a) the temperature above 18˚C prevailed over southern coastal area, while the northern most areas showed the average temperature below 0˚C. The comparison of 1950s (<xref ref-type="fig" rid="fig2">Figure 2</xref>(a)) and 1960s (<xref ref-type="fig" rid="fig2">Figure 2</xref>(b)) shows that the 12˚C isotherms has been moved down towards east in Federally Administered Tribal Areas (FATA) and the belt of temperature which was over the mountains in 1950s was shifted to the piedmonts in 1960s. The shift of isotherms is clear indication that there was change in average temperature in Pakistan at local scale in various decades.</p><p>Reference to <xref ref-type="fig" rid="fig3">Figure 3</xref>(a) most of Balochistan shows temperature in the ranges below average (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x9.png" xlink:type="simple"/></inline-formula>&amp; −2) except Kalat and Khuzdar region where the temperature is above average with the temperature coefficient range of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x10.png" xlink:type="simple"/></inline-formula> &amp; 1 (above average), <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x11.png" xlink:type="simple"/></inline-formula>represents temperature change. In most parts of the Indus</p><fig-group id="fig2"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> (a)-(f): The spatial distribution of average temperatures of December their shift in various decades in Pakistan over the decades of 1950s to 2000s.</title></caption><fig id ="fig2_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x12.png"/></fig><fig id ="fig2_2"><label>(c)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x13.png"/></fig><fig id ="fig2_3"><label> (d)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x14.png"/></fig><fig id ="fig2_4"><label>(e)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x15.png"/></fig><fig id ="fig2_5"><label> (f)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x16.png"/></fig><fig id ="fig2_6"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x17.png"/></fig></fig-group><p>plains (in both Punjab and Sindh) and Khyber-Pakhtoonkhwa (KPK), the observed anomalies were</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x18.png" xlink:type="simple"/></inline-formula>. Nevertheless, Bahawalpur region, FATA, southeastern parts of Sindh falls with in the anomalous temperature coefficient of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x18.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x19.png" xlink:type="simple"/></inline-formula>. The Potwar Plateau reflects slight cooling tendency acknowledged by<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x18.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x20.png" xlink:type="simple"/></inline-formula>.</p></sec><sec id="s3_2"><title>3.2. The Decade of 1960s</title><p>The area in Indus plains between the isotherms of 15˚C and 12˚C is wider (<xref ref-type="fig" rid="fig2">Figure 2</xref>(b)) than in case of 1950s</p><fig-group id="fig3"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> (a)-(f): The spatial distribution of December temperature coefficient in the decades of 1950s to 2000s, the figure depicts the temperature variability at local scale in Pakistan.</title></caption><fig id ="fig3_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x21.png"/></fig><fig id ="fig3_2"><label>(c)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x22.png"/></fig><fig id ="fig3_3"><label> (d)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x23.png"/></fig><fig id ="fig3_4"><label>(e)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x24.png"/></fig><fig id ="fig3_5"><label> (f)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x25.png"/></fig><fig id ="fig3_6"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2360206x26.png"/></fig></fig-group><p>(<xref ref-type="fig" rid="fig2">Figure 2</xref>(a)) the same is true for Balochistan also while the area between 18˚C and 15˚C is wider in Sindh in 1950s than in case of 1960s, this confirms the changes in average temperatures at local levels.</p><p>Now consider the temperature anomalies (<xref ref-type="fig" rid="fig3">Figure 3</xref>(b)) most of Pakistan was under cooling temperatures <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x27.png" xlink:type="simple"/></inline-formula> including most of Sindh, Punjab, Azad Kashmir and Balochistan except the Gawadar area. Most of KPK and FATA were characterized by slight negative coefficient signatures <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x28.png" xlink:type="simple"/></inline-formula> except its upper parts in Hindukush Mountains. In the extreme northeast of the country in HKH mountains temperature was found below average<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x29.png" xlink:type="simple"/></inline-formula>. The Quetta valley and surroundings experienced temperature</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x30.png" xlink:type="simple"/></inline-formula>well below reference point in 1960s, that belt of temperature change extends up to Noshki region close to Pak-Afghan border.</p></sec><sec id="s3_3"><title>3.3. The Decade of 1970s</title><p>The comparison between 1960s and 1970s acknowledged that in 1970s (<xref ref-type="fig" rid="fig2">Figure 2</xref>(c)), the isotherm of 18˚C has been pushed northward and the same is true for 15˚C isotherm also. The 12˚C, 9˚C, 6˚C, 3˚C and 0˚C isotherms also recorded up slope shift in the mountains. It confers the warming condition especially in rugged parts during the decade of 1970s.</p><p>In the decade of 1970s (<xref ref-type="fig" rid="fig3">Figure 3</xref>(c)), many parts of the country like KPK, FATA, Jammu &amp; Azad Kashmir (AJK), lower Indus Plains in southern Punjab and upper Sindh experienced temperature change above average acknowledged by<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x31.png" xlink:type="simple"/></inline-formula>. The highest warming <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x32.png" xlink:type="simple"/></inline-formula> of the decade in December was found in Kalat and Surroundings. The KKH, upper Indus Plains, southwestern Balochistan and southern Sindh were dominated by temperature change below average with confidence level of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x33.png" xlink:type="simple"/></inline-formula>. The Quetta Valley and surrounding were found with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x34.png" xlink:type="simple"/></inline-formula> anomalous temperatures and noticed to be the coolest area of the decade.</p></sec><sec id="s3_4"><title>3.4. The Decade of 1980s</title><p>The coastal areas reflected variation of temperature averages almost in all decades, here one thing can not be ignored that is the impact of changes in sea surface temperature (SST) of Arabian Sea [<xref ref-type="bibr" rid="scirp.52628-ref18">18</xref>] . The 12˚C isotherm has climbed slightly up slopes in piedmonts areas (<xref ref-type="fig" rid="fig2">Figure 2</xref>(d)). However, the spacing of isothermal lines is more or less analogous in 1970s and 1980s.</p><p>The 1980s (<xref ref-type="fig" rid="fig3">Figure 3</xref>(d)) registered no obvious warming surface temperatures and most of the country was dominated by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x35.png" xlink:type="simple"/></inline-formula> that shows temperature slightly below average. The Chitral in KPK and an elongated patch from Mastung-Kalat area to Makran Hills in Balochistan experienced temperature decrease <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x36.png" xlink:type="simple"/></inline-formula> and the temperature was obviously low than their surrounding territories. The only temperature change of the decade with coefficient value below −1 was registered in Sulaiman Ranges.</p></sec><sec id="s3_5"><title>3.5. The Decade of 1990s</title><p>The comparison of 1980s (<xref ref-type="fig" rid="fig2">Figure 2</xref>(d)) and 1990s (<xref ref-type="fig" rid="fig2">Figure 2</xref>(e)) supports that a northward as well as northwestward shift of isotherms has occurred in 1990s. The northward shift was not only true for mountains but also obvious in plains and piedmonts which is clear indication of warming. Both the isothermal shift and temperature anomalies confirmed that the mountainous areas of Pakistan face more temperature changes (warming) than plains.</p><p>The temperature anomaly in 1990s (<xref ref-type="fig" rid="fig3">Figure 3</xref>(e)) was not as uniform over the country as in the decade of 1980s. Abrupt temperature changes were noticed in 1990s where almost all Pakistan was dominated by warming condition. The Sulaiman Ranges and surrounding was below average in 1980s but in 1990s opposite situation persisted with obvious warming tendency<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x37.png" xlink:type="simple"/></inline-formula>. Maximum portion of the national territory configures the change of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x38.png" xlink:type="simple"/></inline-formula>. The Hindukush Mountains in Chitral, parts of Himalayas, parts of eastern Punjab along the Pak-Indian border southern most Sindh, central and northern Balochistan were subjected to warming and have shown changes in temperature within the range of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x39.png" xlink:type="simple"/></inline-formula>.</p></sec><sec id="s3_6"><title>3.6. The Decade of 2000s</title><p>The shift of isothermal lines in 2000s (<xref ref-type="fig" rid="fig2">Figure 2</xref>(f)) towards north and northwest was even more than the case of 1990s. The shift of high temperature averages towards north is devastating for the big valley glaciers in HKH and will badly affect the agriculture-based economy of the country based on the Indus Rivers System.</p><p>Reference to the decade of 2000s (<xref ref-type="fig" rid="fig3">Figure 3</xref>(f)), except FATA, southern KPK, Sulaiman Ranges and Noshki with surroundings the rest of Pakistan is under warming condition. The amount of warming observed was significant as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x40.png" xlink:type="simple"/></inline-formula> in HKH, KPK, area surrounding Sulaiman Ranges lobe, eastern parts of Sindh, Kharan region in Balochistan and upper Indus plains while temperature changes of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x41.png" xlink:type="simple"/></inline-formula> was observed in southeastern Punjab and Kirthar Ranges. The upper Indus plains, parts of lower Indus plains, most of Balochistan especially its piedmonts and plains exhibits the change of temperature in December with the range of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2360206x42.png" xlink:type="simple"/></inline-formula>. Generally, based on results, as one proceeds down from mountainous territories to the plains, the temperature coefficient values are increasing. Therefore, in 2000s warming was much obvious in the plains than mountainous regions while in 1990s the mountains regions have shown more warming than the plains.</p></sec></sec><sec id="s4"><title>4. Conclusions</title><p>In Pakistan, the temperature anomalies in December vary from decade to decade and region to region within the same decade. It was noticed that after 1980s, the warming tendency was obvious in most of the national territories. Nevertheless, the regional detail of temperature anomalies at local scale were complex and did not comply with the assumption totally that winter months were warming. Generally, by shift from mountainous territories to the plains, the temperature coefficient values are generally increasing. Therefore, in the 2000 decade warming was much obvious in the plains than that in mountainous regions while in 1990s the mountains regions showed more warming than the plains. This sort of study could be useful for local agriculture especially in various enclosed valleys in rugged parts of the country.</p><p>It has been observed that wherever the isotherms were closely spaced, the temperature anomalies found were with high statistical significance on the contrary where the isotherms were widely spaced and the temperature change per unit area was less. Generally, it is clear that the shift of isotherm was obvious in the northern and western parts of the study domain acknowledged by high level of statistical significance pertaining to temperature anomalies there. Both the isothermal shift and temperature anomalies confirm that the mountainous areas of Pakistan face more temperature variability (warming) than plains.</p></sec><sec id="s5"><title>Acknowledgements</title><p>The help and cooperation of PMD is highly recognized. The comments and suggestions of anonymous reviewers substantially improved the paper.</p></sec></body><back><ref-list><title>References</title><ref id="scirp.52628-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">UNFCCC (2007) Climate Change: Impacts, Vulnerabilities and Adaptation in Developing Countries. http://unfccc.int/resource/docs/publications/impacts</mixed-citation></ref><ref id="scirp.52628-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Pachauri, R.K. (2009) Regional Conference on “Climate Change: Challenges and Opportunities for South Asia”. 13-14 January 2009, Islamabad, Pakistan.</mixed-citation></ref><ref id="scirp.52628-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Pakistan Meteorological Department (PMD) (2009) Climate Change Indicators of Pakistan, Technical Report No. PMD 22/2009 Published by PMD, Islamabad Pakistan. http://www.pakmet.com.pk</mixed-citation></ref><ref id="scirp.52628-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Easterling, D.R., Evans, J.L., Groisman, P.Y., Karl, T.R., Kunkel, K.E. and Ambenje, P. (2000) Observed Variability and Trends in Extreme Climate Events: A Brief Review. Bulletin of the American Meteorological Society, 81, 417-425. http://dx.doi.org/10.1175/1520-0477(2000)081&lt;0417:OVATIE&gt;2.3.CO;2</mixed-citation></ref><ref id="scirp.52628-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Mahrt, L. (2006) Variation of Surface Air Temperature in Complex Terrain. Journal of Applied Meteorology and Climatology, 45, 1481-1493. http://dx.doi.org/10.1175/JAM2419.1</mixed-citation></ref><ref id="scirp.52628-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Ahmad, I., Ambreen, R., Sultan, S., Sun, Z. and Nawaz, M. (2014) Regional Characteristics of Temperature Anomalies in Pakistan with Emphasis on Spatial Distribution at Decadal Scale: A Case Study of August (1950s-2000s). Atmospheric and Climate Sciences, 4, 721-726. http://dx.doi.org/10.4236/acs.2014.44065</mixed-citation></ref><ref id="scirp.52628-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Buhairi, M. (2010) Analysis of Monthly, Seasonal and Annual Air Temperature Variability and Trends in Taiz City— Republic of Yemen. Journal of Environmental Protection, 1, 401-409. http://dx.doi.org/10.4236/jep.2010.14046</mixed-citation></ref><ref id="scirp.52628-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Ministry of Finance Government of Pakistan (2012) 2012 Pakistan Economic Survey, 2011-12. http://www.finance.gov.pk/survey/chapter_12/16-Environment</mixed-citation></ref><ref id="scirp.52628-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Ministry of Finance Government of Pakistan (2013) 2013 Pakistan Economic Survey, 2012-13. http://www.finance.gov.pk/survey_1213.html</mixed-citation></ref><ref id="scirp.52628-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Hussain, S.S. and Mudasser, M. (2007) Prospects for Wheat Production under Changing Climate in Mountain Areas of Pakistan—An Econometric Analysis. Agricultural Systems, 94, 494-501. http://dx.doi.org/10.1016/j.agsy.2006.12.001</mixed-citation></ref><ref id="scirp.52628-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Himayatullah, K. and Abuturab, K. (2008) Natural Hazards and Disaster Management in Pakistan. Munich Personal RePEc Archive.</mixed-citation></ref><ref id="scirp.52628-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">PADMU (1983) Country Report Pakistan—Desertification Problems, Extent and Remedial Measures. Pakistan Desertification Monitoring Unit (PADMU), Islamabad.</mixed-citation></ref><ref id="scirp.52628-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Griggs, D.J. and Noguer, M. (2002) Climate Change 2001: The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change. Weather, 57, 267-269.http://dx.doi.org/10.1256/004316502320517344</mixed-citation></ref><ref id="scirp.52628-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Nordhaus, W.D. (2007) A Review of the Stern Review on the Economics of Climate Change. Journal of Economic Literature, 45, 686-702. http://dx.doi.org/10.1257/jel.45.3.686</mixed-citation></ref><ref id="scirp.52628-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Ahmad, I., Zhaobo, S., Deng, W. and Ambreen, R. (2010) Trend Analysis of January Temperature in Pakistan over the Period of 1961-2006: Geographical Perspective. Pakistan Journal of Meteorology, 13, 11-22.</mixed-citation></ref><ref id="scirp.52628-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Ahmad, I., Ambreen, R., Sultan, S., Sun, Z.B. and Deng, W.T. (2014) Spatial-Temporal Variations in January Temperature in Pakistan and Their Possible Links with SLP and 500-hPa Levels over the Period of 1950-2000: A Geographical Approach. Atmospheric and Climate Sciences, 4, 524-533. http://dx.doi.org/10.4236/acs.2014.44048</mixed-citation></ref><ref id="scirp.52628-ref17"><label>17</label><mixed-citation publication-type="book" xlink:type="simple">Cruz, V., Harasawa, H., Lal, M., Wu, S., Anokhin, Y., Punsalmaa, B., Honda, Y., Jafari, M., Liand, C. and Huu Ninh, N. (2007) Asia. Climate Change 2007: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Parry, M.L., Canziani, O.F., Palutikof, J.P., vander Linden, P.J. and Hanson, C.E., Eds., Cambridge University Press, Cambridge, UK, 471.</mixed-citation></ref><ref id="scirp.52628-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Khan, T.M.A., Khan, F.A. and Jilani, R. (2008) Sea Surface Temperature Variability along Pakistan Coast and Its Relation to El Nino-Southern Oscillation. Journal of Basic and Applied Sciences, 4, 67-72.</mixed-citation></ref></ref-list></back></article>