PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 7350.
Published online 2017 August 4. doi:  10.1038/s41598-017-07627-5
PMCID: PMC5544675

Unravelling Diurnal Asymmetry of Surface Temperature in Different Climate Zones

Abstract

Understanding the evolution of Diurnal Temperature Range (DTR), which has contradicting global and regional trends, is crucial because it influences environmental and human health. Here, we analyse the regional evolution of DTR trend over different climatic zones in India using a non-stationary approach known as the Multidimensional Ensemble Empirical Mode Decomposition (MEEMD) method, to explore the generalized influence of regional climate on DTR, if any. We report a 0.36 °C increase in overall mean of DTR till 1980, however, the rate has declined since then. Further, arid deserts and warm-temperate grasslands exhibit negative DTR trends, while the west coast and sub-tropical forest in the north-east show positive trends. This transition predominantly begins with a 0.5 °C increase from the west coast and spreads with an increase of 0.25 °C per decade. These changes are more pronounced during winter and post-monsoon, especially in the arid desert and warm-temperate grasslands, the DTR decreased up to 2 °C, where the rate of increase in minimum temperature is higher than the maximum temperature. We conclude that both maximum and minimum temperature increase in response to the global climate change, however, their rates of increase are highly local and depend on the underlying climatic zone.

Introduction

Global warming due to the rise in greenhouse gas emissions is usually quantified through the long-term changes in the mean temperature (T mean) of the earth’s surface, expressed in yearly, decadal or even centurial timescales. However, regional scale impact assessment, adaptation, and mitigation of climate change need an evaluation of suitable climate change indicators at daily time scales. Diurnal Temperature Range (DTR), defined as the difference between daily maximum and minimum temperatures, is a key driving factor in most of the climatological processes13 and is an important index of diurnal variations4. The intensity, frequency, and duration of extreme precipitation events are increasing, on average, with a decrease in DTR due to increase in night time temperatures than daytime temperatures59. Further, any change in DTR would increase the risks associated with drought10, heat stress11, crop failure12, 13, human health14 and mortality rate15.

Numerous studies1, 16, 17 reported a global decrease in DTR at a rate of around 0.07 °C per decade18, with the decrease being evident during 1950–198019, 20 and attributed this to the higher rate of increase of minimum temperature (T min) than maximum temperature (T max). However, a few regional studies, mainly in the tropics and subtropics, reported an increasing trend21, 22 in DTR and attributed it solely to the increase in T max, which may be ascribed to the influence of regional climate on the behaviour of DTR. Further, ref. 20 highlighted that the acuteness of DTR is uncertain, especially in the south Asian regions such as India and China, arising mainly due to the unavailability of reliable regional data products. In this study, we analysed the evolution of T mean, T max, T min and DTR over India using the recently available 1°×1° gridded daily temperature data provided by Indian Metrological Department (IMD)23. Detailed analysis of DTR is crucial over this region owing to the associated risks12, 13, 15, aggravated by an agriculture-dependent economy and dense population. Spatio-temporal exploration of DTR in different heterogeneous climatic zones24 in India, i.e., arid/semi-arid in north-west, warm-temperate in central India, sub-tropical humid in north-east and foothills of Himalayas, semi-arid and equatorial grassland in southern India and warm-humid coastal regions, is expected to reveal the generalized influence of regional climate on DTR.

Past studies over India2, 17, 25 attributed an increasing trend in DTR to an increase in T max and no/insignificant change in T min. On the contrary, few studies26, 27 observed a decreasing trend in both T min and T max over western Himalayas. Further, ref. 28 used high-resolution temperature data and observed no significant trend in DTR, except significant decrease over northernmost India. These contradictions can be attributed to the limitations in traditional trend analysis. Moreover, traditional statistical trend approaches employed in the past consider either a constant linear trend or a shape established a priori 29. Such approaches are sensitive to the boundaries of the time series and hence fail to capture the multi-decadal variability30. To eliminate these disadvantages and to understand the non-stationary behaviour, we use Multidimensional Ensemble Empirical Mode Decomposition (MEEMD)31, which is an improved method for scientific data analysis (See Methods). Considering the availability of data and in order to compare the present study with the previous studies27, 32, we restrict the analysis to a period of 60 years from 1951 to 2010.

Result and Discussion

The regional average of T mean, T max, T min and DTR over India is shown in Fig. 1a. T mean over India shows an increasing trend since 1951. An identical pattern is observed for T max. However, T min decreases slightly till 1980 with a magnitude of 0.19 °C, which results in subsequent increase (0.36 °C) in DTR. Thereafter, insignificant trend in DTR can be attributed to the rapid increase in T min. Figure 1b shows the 60-year average of DTR, which ranges from 7 °C to 16 °C, with its values gently increasing from humid coastal regions to north-western arid regions. To understand the temporal evolution of temperature, we evaluated the spatial EEMD-trend (See Methods for Definition) of annual DTR, T max, and T min for four distinct time windows (1951–1980, 1951–1990, 1951–2000 and 1951–2010), as shown in Fig. 2. DTR (Fig. 2a, top panel) shows a positive trend over southern west coast, foothills of the Himalayas and some parts of north-west, which is analogous to results from a previous study32. The spatial extent of positive trend is found to gradually increase over the Equatorial and north-eastern regions, with an approximate increase of 0.25 °C per decade. On the contrary, significant negative trend is observed over some parts of northern west coast during the initial 30 years, which has gradually extended towards central and northern India during the subsequent decades. The spatial extent and magnitude of the positive trend over the north-west also decreases over the decades. Negative trends in DTR over north-west and central India, known for its arid/semi-arid climate, well agree with the behaviour of climatologically similar regions taken up by previous studies33, 34. A similar analysis of T max reveals the substantial increase in the spatial extent and the magnitude of the positive trend over the southern west coast (Fig. 2b, middle panel). T max of the entire coastal region is increasing at a rate of 0.25 °C per decade. Similar intensification of warming is visible in the north-east during the recent decades. This explains the increase in DTR over these regions. However, a reversal of trend pattern is observed for T min (Fig. 2c, bottom-panel). While, the negative trend of T min is manifested in the foot-hills of Himalayas26 in the initial years, it has shrunk gradually over the decades, along with the concurrent evolution of its positive trend in the north-west after 1990. Interestingly, positive and negative trend patterns of DTR are consistent with the positive trends in T max and T min, respectively27.

Figure 1
(a) Annual variation of spatially averaged annual temperature over India and (b) Average DTR from 1951 to 2010. Grids with inconsistent data are shown in grey colour (extreme bottom of colour bar), which is represented as NA (Not Analysed). The maps were ...
Figure 2
Spatial evolution of Ensemble Empirical Mode Decomposition trend of annual (a) DTR, (b) Maximum and (c) Minimum temperatures. Each sub panel has four time windows 1) 1951–1980, 2) 1951–1990, 3) 1951–2000 and 4) 1951–2010. ...

We also explored the DTR trends in different seasons to understand the influence of seasons on DTR in various regions. MEEMD analysis is repeated for four prominent seasons, i.e., winter (DJF: December to February), pre-monsoon (MAM: March to May), summer monsoon (JJA: June to August) and post-monsoon (SON: September to November) as shown in Fig. 3. During winter, DTR shows negative trend in central India and few parts of the north-west; and positive trend in west coast and the north-east32 (Fig. 3a, top panel). The spatial extent of this negative trend has moved towards the entire north-west and central India over the decades, with a decrease of at most 2 °C in most of these regions by 2010. Likewise, the positive trend of DTR in west coast has increased by at most 3 °C covering the entire Peninsular India. Pre-monsoon season shows a constant positive change in DTR in most of the grids (Fig. 3, 2nd panel). It is evident that there is no significant change, except approximately 0.5 °C/decade increase in a few grids of Deccan Plateau. Positive trend dominates during monsoon season too, except over central India. However, a manifestation of considerable evolution is absent over southern India, east coast and north-east (Fig. 3, 3rd panel). Moreover, a reversal of trend is seen over the entire northern India by 2010. The trend in the post-monsoon period (Fig. 3, bottom panel), follows the same pattern as that of the winter with more dominance of positive trend over the region.

Figure 3
Spatial evolution of Ensemble Empirical Mode Decomposition trend of seasonal DTR for (a) DJF, (b) MAM, (c) JJA and (d) SON. Each sub panel has four time windows 1) 1951–1980, 2) 1951–1990, 3) 1951–2000 and 4) 1951–2010. ...

The seasonal variations in DTR are further understood with an investigation of the corresponding variation of T max and T min (Figures S2 and S3). This reveals that the EEMD-trend is similar to that of annual pattern, with different spatial extent and magnitude (for a detailed description see Supplementary Section-2). Especially, the increase in T min (~3 °C after 1950, Figure S3) during winter, particularly in the agriculturally rich areas of north-west, will severely affect the yield of winter crops12 and the warming phase of T min, may also affect the yield of paddy in other parts of India, if prolonged13. Intensification of warming, irrespective of seasons, is apparent from the spatial pattern of various aspects of T max and T min extremes, computed using extreme event indices defined by CLIMDEX (see Table S1). High (low) extremes of T max and T min are found to be following an increasing (decreasing) trend over the decades (Figures S4 and S5). This suggests that both T max and T min are increasing, however, their respective rate of increase determines the DTR trend. This relationship between DTR and the temperature extremes is further analysed through scatter plots of DTR against the respective T max and T min for the grids with positive and negative trends separately (Fig. 4). It shows that the positive trend in DTR is because of a higher rate of increase in T max than T min, while the reverse is also true for the negative trend of DTR.

Figure 4
Scatter plot of Ensemble Empirical Mode Decomposition trend in a) positive trend of DTR and T max (T min) and b) negative trend of DTR and T max (T min) for the last time window (1951–2010). Unit is in °C. The maps were generated using ...

The influence of regional climate on the relationship between DTR, T max, and T min is investigated by comparing their EEMD-trends for different climatic zones and for different time windows (Fig. 5). For this analysis, representative grids were taken from these climatic zones. We observed that the north-western arid and semi-arid regions; and foothills of Himalayas showed a strengthening of negative DTR trends with an increase in the time window, which is because of a strengthening of positive trends in T min compared to T max, likewise the reverse is also true. Regions with tendency of DTR towards positive trends (north-east, semi-arid Deccan plateau, Equatorial grasslands and west coast) owe it to the tendency of T max to move towards positive trends compared to T min. The increase in DTR caused by the increase in T max can result in precipitation deficit6, 10. Especially in north-east and the Western Ghats, which are high rainfall receiving regions, the positive trend of DTR may signal the increase in dry days30. This suggests that the difference in the regional evolution of DTR is influenced by the different regional response of T max and T min. Clearly, T min is increasing at a faster pace in the north-west. Even the magnitude of the negative trend at other places has significantly reduced over the decades. Traditional methods of linear trend analysis fail to capture this increasing phase2, 26, 27. For statistical justification, the significance of trends are checked at 10% significance level and the grids showing significant trends for all variables are shown in Figure S6.

Figure 5
Regional Ensemble Empirical Mode Decomposition trend of DTR, T max and T min for heterogeneous climatic zones (a) 1951–1980, (b) 1951–1990, (c) 1951–2000 and (d) 1951–2010, (e) representative grids in different climatic ...

The regional patterns of DTR depend on the rate of increase in T min and T max. In India, the regional climate affects this rate of increase through many different factors like cloud cover1, 2, 12, 28, precipitation1, 2, 34, soil moisture1, 2, vegetation cover34, etc. All the above processes are inversely correlated with DTR. However, these feedback processes are local in nature, which are consequently stimulated by some global factors. For instance, DTR over northern India decreases with increase in cloud cover formed due to aerosol loading2, 12, 28, 35. However, a contradictory pattern of DTR is observed over southern India under high cloud cover12, 28. This behaviour can be attributed to the proximity to the sea and the rising trend of Sea Surface Temperature (SST)12, as significant positive correlation between SST and near surface air temperature exists35, 36. Moreover, we observe an increasing trend of near surface air temperature in coastal regions, where the effect of urbanization is not pertinent32. To understand this behaviour of diurnal temperature variations, we also analysed the SST (see supplementary information). Since SST directly depends on the ocean heat content, which is affected by the global mean temperature37, it is pertinent to suggest that the effect of the global climate on SST is transferred to the observed changes in the land surface temperature. Figure S7 reveals that the rate of warming of the Arabian Sea is higher than the Bay of Bengal, which could explain the evolution of DTR from west to east over time. The difference in the SST of Arabian Sea and Bay of Bengal can impact the surface temperature behaviours over land, while it is also affected by the local factors like vegetation cover of the region38. The increase in humidity in the arid deserts and warm-temperate grasslands due to the warming of the Arabian Sea and increase in cropland cover39 can explain the significant increase in T min, which leads to a decrease in DTR in these regions. In case of west coast, the warmed up sea causes both T max and T min to rise, however the higher rate of increase in T max than T min can be attributed to the decline in forest cover. The diurnal asymmetry pattern over the sub-tropical forest in the north-east can be attributed to the warming up of Bay of Bengal. From the above discussion, it is evident that the rise in both T min and T max in all climatic zones is in line with global climate change, whereas the rate at which they increase is strongly influenced by the regional factors, with clear distinction among the climate zones. However, a detailed study investigating the factors affecting DTR changes is needed to better understand this behaviour of diurnal temperature variations over India.

Conclusion

This study presents the analysis of Diurnal Temperature Range (DTR) and temperature employing a non-stationary and non-linear approach known as Multi-dimensional Ensemble Empirical Mode Decomposition (MEEMD) method, using gridded daily temperature data of past 60 years (1951–2010) over India. To further understand the influence of regional climate on DTR, spatio-temporal exploration of DTR over different heterogeneous climatic zones are performed. The following conclusions are derived from this study:

  1. The regions conventionally showing high DTR such as arid, semi-arid and warm temperate grasslands of India, show negative trend in DTR, because the rate of increase in minimum temperature is higher than that of the maximum temperature.
  2. Sub-tropical forests, equatorial grasslands down south and the west coast show a higher rate of increase in maximum temperature with no/negligible change in minimum temperature, resulting in positive trend of DTR.
  3. The season-wise analysis reveals that the evolution in trend of DTR is predominantly evident during winter (DJF) and post-monsoon (SON) seasons, indicating warmer winters.
  4. Analysis of extremes shows an increasing (decreasing) warm (cold) days and nights. In addition, the regional trend of minimum temperature suggests a shift towards global pattern of increasing night time temperature, contrasting the pattern reported in the recent literature.
  5. The increase in minimum, maximum and mean temperatures can be attributed to global climate change. However, the changes in DTR mainly depend on the relative rate of increase in maximum and minimum temperatures, which is chiefly dominated by the regional climate.

While previous studies mainly attributed DTR changes to local factors1, 2, 28, our analysis on DTR over India with SST reveals that both global and local factors are working in tandem. We argue that the significant decrease in DTR in arid desert and warm-temperate grasslands can be attributed to the global influence on the Arabian Sea and subsequent increase in humidity over these regions. However, the rate of increase of T min, which is a factor of local micro-climate, also plays a dominating role leading to a negative trend in DTR over these regions. In case of the west coast, the Arabian Sea warming has a direct influence on the observed positive trend in DTR since the warmed up sea causes the T max to rise faster than the T min. However, plausible the above reasons may sound, further detail analysis is required to obtain a clear picture of the factors influencing diurnal asymmetry, both local and global.

Methods

Data Description

High resolution, 1° × 1° gridded, daily maximum, minimum and mean temperature prepared by National Climatic Centre (NCC), Indian Meteorological Department (IMD) is used in this study23. The spatial domain of this gridded data set is 7.5°N to 37.5°N and 67.5°E to 97.5°E, covering the main land region of India. After proper quality checks, the daily temperature record obtained from 395 stations are selected for the period of 1951–2013. Further, the most appropriate interpolation technique for irregular spaced data40, which is, modified version of Shepard’s angular distance weighing algorithm41, is used to convert the station data into gridded data. The accuracy of the gridded dataset evaluated using cross-validation and the RMSE (Root Mean Square Error) was found to be lesser than 0.5 °C. Moreover, the present dataset is compared with the monthly mean temperature dataset prepared University of Delaware42 and correlation coefficient is noted to be around 0.8 in most of the grids23. This dataset is found to be useful for extreme temperature analysis, since the occurrence of historical heat waves exactly associate with the present data set. The all-India average of maximum, minimum and mean temperature varies from 20.48 °C to 34.18 °C, 7.7 °C to 24.4 °C and 14.1 °C to 28.9 °C respectively. This data set is extensively used in the literature, especially for analysis of extremes11, 43.

Extreme Event Indices

The extreme event indices adopted in the analysis have been defined based on the ‘Expert Team on Climate Change Detection and Indices’ (ETCCDI) recommendation44 (available from http://www.climdex.org/indices.html). Totally 9 indices are analyzed out of 16 indices proposed44 by ETCCDI. The other indices such as GSL (Growing Season Length), FD (Number of Frost Days), ID (Number of Icing Days) are not suitable indices for Indian weather. Instead of taking the extremes of maximum and minimum temperature, we calculated the average annual/seasonal maximum and minimum temperature to analyze the changes in DTR. The percentile threshold value is defined as the calendar day 90th percentile centered on a 5-day window for the base period 1961-1990. Bootstrap approach is used to calculate the number of days in the base period to avoid homogeneity in the in-base and out-base period45. Indices used in this study are shown in Table S1.

Multidimensional Ensemble Empirical Mode Decomposition

Multidimensional Ensemble Empirical Mode Decomposition (MEEMD)31 is a method for eliminating the oscillatory component of a time series and revealing while preserving the slow varying component46. This approach does not need a functional form a priori and is capable of extracting the hidden nonlinear and non-stationary nature of the time series29, 47. In MEEMD, a time series, Y(t), at a grid point, is decomposed using EEMD46, 48 into oscillatory components known as Intrinsic Mode Functions (IMF) as shown in equation (1). The sequential elimination of IMFs (I j) produces residual (R n), which is either monotonic or contains only one extremum, and cannot be further decomposed into an oscillatory component.

Y(t)=j=1nIj(t)+Rn(t)
1

A detailed methodology is presented in the Supplementary Section-1. MEEMD is less sensitive to the boundary compared to the other commonly used methods, hence the physical meaning of the time series will be retained even after appendage of data49. MEEMD has been successfully implemented in various climate studies29, 49 due to its effectiveness and the spatial-temporal locality. Since the trend varies with time, as described above, the averaged slope of the trend within a time window cannot appropriately capture the evolution of the trend. Hence, the trend at a particular time(t) is represented as the increment of the EEMD-trend (Residual)29 from that of the reference time (1951) as shown in equation (2).

TrendEEMD(t) = Rn(t) − Rn(1951)
2

Test for Significance

The statistical significance of these EEMD trends at a given spatio-temporal locations is evaluated by Monte Carlo Simulation adopted by ref. 29. Here, 10000 samples of red noise series were generated, having the same temporal length (60) and lag-1 auto-correlation as that of the original series. The slow varying component of each of the generated series is computed using MEEMD approach and the probability density function (PDF) of the EEMD. The trend at any time is calculated. Then the trend value (ratio between EEMD and standard deviation of the respective variable) of the particular grid is checked whether it is within 10% significance level of the empirical PDF. If it is satisfied, then the EEMD trend is considered to be statistically significant. However, statistically non-significant trend does not mean a complete rejection of the null hypothesis, but only indicates that there is no clear evidence for a trend.

Data Availability statement

The data that support the findings of this study are available from Indian Meteorological Department (IMD). More details are available in http://www.imd.gov.in/advertisements/20170320_advt_34.pdf.

Electronic supplementary material

Author Contributions

Author Contributions

C.T.D. designed the study. R.V. developed the software code, carried out the analysis and plotted the figures. R.V. and A.C. interpreted the results. All the authors wrote the manuscript.

Notes

Competing Interests

The authors declare that they have no competing interests.

Footnotes

Electronic supplementary material

Supplementary information accompanies this paper at doi:10.1038/s41598-017-07627-5

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Dai A, Trenberth KE, Karl TR. Effects of Clouds, Soil Moisture, Precipitation, and Water Vapor on Diurnal Temperature Range. J. Clim. 1999;12:2451–2473. doi: 10.1175/1520-0442(1999)012<2451:EOCSMP>2.0.CO;2. [Cross Ref]
2. Rai A, Joshi MK, Pandey AC. Variations in diurnal temperature range over India: Under global warming scenario. J. Geophys. Res. Atmos. 2012;117:1–12. doi: 10.1029/2011JD016697. [Cross Ref]
3. Davy, R., Esau, I., Chernokulsky, A., Outten, S. & Zilitinkevich, S. Diurnal asymmetry to the observed global warming. Int. J. Climatol., doi:10.1002/joc.4688 (2016).
4. Cook, T. Toward a Reassessment of Daily Temperature Range Trends. Eos (Washington. DC). 97 (2016).
5. Donat MG, et al. Updated analyses of temperature and precipitation extreme indices since the beginning of the twentieth century: The HadEX2 dataset. J. Geophys. Res. Atmos. 2013;118:2098–2118. doi: 10.1002/jgrd.50150. [Cross Ref]
6. Miao C, Sun Q, Borthwick AGL, Duan Q. Linkage Between Hourly Precipitation Events and Atmospheric Temperature Changes over China during the Warm Season. Sci. Rep. 2016;6:22543. doi: 10.1038/srep22543. [PMC free article] [PubMed] [Cross Ref]
7. Shrestha, A. B., Bajracharya, S. R., Sharma, A. R., Duo, C. & Kulkarni, A. Observed trends and changes in daily temperature and precipitation extremes over the Koshi river basin 1975–2010. Int. J. Climatol., doi:10.1002/joc.4761 (2016).
8. Singh D, Tsiang M, Rajaratnam B, Diffenbaugh NS. Observed changes in extreme wet and dry spells during the South Asian summer monsoon season. Nat. Clim. Chang. 2014;4:456–461. doi: 10.1038/nclimate2208. [Cross Ref]
9. Klein Tank AMG, et al. Changes in daily temperature and precipitation extremes in central and south Asia. J. Geophys. Res. 2006;111:D16105. doi: 10.1029/2005JD006316. [Cross Ref]
10. He B, Huang L, Wang Q. Precipitation deficits increase high diurnal temperature range extremes. Sci. Rep. 2015;5:12004. doi: 10.1038/srep12004. [PMC free article] [PubMed] [Cross Ref]
11. Dash SK, Kjellstrom T. Workplace heat stress in the context of rising temperature in India. Curr. Sci. 2011;101:496–503.
12. Jaswal AK. Recent winter warming over India - spatial and temporal characteristics of monthly maximum and minimum temperature trends for January to March. Mausam. 2010;61:163–174.
13. Welch JR, et al. Rice yields in tropical/subtropical Asia exhibit large but opposing sensitivities to minimum and maximum temperatures. Proc. Natl. Acad. Sci. 2010;107:14562–14567. doi: 10.1073/pnas.1001222107. [PubMed] [Cross Ref]
14. Cheng J, et al. Impact of diurnal temperature range on human health: a systematic review. Int. J. Biometeorol. 2014;58:2011–24. doi: 10.1007/s00484-014-0797-5. [PubMed] [Cross Ref]
15. Briga M, Verhulst S. Large diurnal temperature range increases bird sensitivity to climate change. Sci. Rep. 2015;5:16600. doi: 10.1038/srep16600. [PMC free article] [PubMed] [Cross Ref]
16. Karl TR, et al. A New Perspective on Recent Global Warming: Asymmetric Trends of Daily Maximum and Minimum Temperature. Bull. Am. Meteorol. Soc. 1993;74:1007–1023. doi: 10.1175/1520-0477(1993)074<1007:ANPORG>2.0.CO;2. [Cross Ref]
17. Alexander LV, et al. Global observed changes in daily climate extremes of temperature and precipitation. J. Geophys. Res. 2006;111:D05109.
18. Trenberth, K. E. et al. In Climate Change 2007: The Physical Science Basis (eds. Solomon, S. et al.) 235–336 (Cambridge University Press, 2007).
19. Vose RS, Easterling DR, Gleason B. Maximum and minimum temperature trends for the globe: An update through 2004. Geophys. Res. Lett. 2005;32:L23822. doi: 10.1029/2005GL024379. [Cross Ref]
20. Thorne PW, et al. Reassessing changes in diurnal temperature range: Intercomparison and evaluation of existing global data set estimates. J. Geophys. Res. Atmos. 2016;121:5138–5158. doi: 10.1002/2015JD024584. [Cross Ref]
21. Peralta-Hernandez AR, Balling RC, Jr, Barba-Martinez LR. Analysis of near-surface diurnal temperature variations and trends in southern Mexico. Int. J. Climatol. 2009;29:205–209. doi: 10.1002/joc.1715. [Cross Ref]
22. Shahid S, Harun S, Bin, Katimon A. Changes in diurnal temperature range in Bangladesh during the time period 1961–2008. Atmos. Res. 2012;118:260–270. doi: 10.1016/j.atmosres.2012.07.008. [Cross Ref]
23. Srivastava AK, Rajeevan M, Kshirsagar SR. Development of a high resolution daily gridded temperature data set (1969-2005) for the Indian region. Atmos. Sci. Lett. 2009;10:249–254.
24. Peel MC, Finlayson BL, McMahon TA. Updated world map of the Köppen-Geiger climate classification. Hydrol. Earth Syst. Sci. 2007;11:1633–1644. doi: 10.5194/hess-11-1633-2007. [Cross Ref]
25. Kothawale DR, Rupa Kumar K. On the recent changes in surface temperature trends over India. Geophys. Res. Lett. 2005;32:1–4. doi: 10.1029/2005GL023528. [Cross Ref]
26. Yadav RR, Park W-K, Singh J, Dubey B. Do the western Himalayas defy global warming? Geophys. Res. Lett. 2004;31:1–5. doi: 10.1029/2004GL020201. [Cross Ref]
27. Kothawale DR, Kumar KK, Srinivasan G. Spatial asymmetry of temperature trends over India and possible role of aerosols. Theor. Appl. Climatol. 2012;110:263–280. doi: 10.1007/s00704-012-0628-8. [Cross Ref]
28. Roy SS, Balling RC. Analysis of trends in maximum and minimum temperature, diurnal temperature range, and cloud cover over India. Geophys. Res. Lett. 2005;32:1–4. doi: 10.1029/2004GL022201. [Cross Ref]
29. Ji F, Wu Z, Huang J, Chassignet EP. Evolution of land surface air temperature trend. Nat. Clim. Chang. 2014;4:462–466. doi: 10.1038/nclimate2223. [Cross Ref]
30. Vinnarasi R, Dhanya CT. Changing characteristics of extreme wet and dry spells of Indian monsoon rainfall. J. Geophys. Res. Atmos. 2016;121:2146–2160. doi: 10.1002/2015JD024310. [Cross Ref]
31. Wu Z, Huang NE, Chen X. The Multi-Dimensional Ensemble Empirical Mode Decomposition Method. Adv. Adapt. Data Anal. 2009;01:339–372. doi: 10.1142/S1793536909000187. [Cross Ref]
32. Kumar KR, Kumar KK, Pant GB. Diurnal asymmetry of surface temperature trends over India. Geophys. Res. Lett. 1994;21:677–680. doi: 10.1029/94GL00007. [Cross Ref]
33. Zhou L, Dickinson RE, Tian Y, Vose RS, Dai Y. Impact of vegetation removal and soil aridation on diurnal temperature range in a semiarid region: Application to the Sahel. Proc. Natl. Acad. Sci. 2007;104:17937–17942. doi: 10.1073/pnas.0700290104. [PubMed] [Cross Ref]
34. Zhou L, et al. Spatial dependence of diurnal temperature range trends on precipitation from 1950 to 2004. Clim. Dyn. 2009;32:429–440. doi: 10.1007/s00382-008-0387-5. [Cross Ref]
35. Kothawale DR, Munot AA, Krishna Kumar K. Surface air temperature variability over India during 1901–2007, and its association with ENSO. Clim. Res. 2010;42:89–104. doi: 10.3354/cr00857. [Cross Ref]
36. Revadekar JV, Varikoden H, Murumkar PK, Ahmed SA. On the relationship between sea surface temperatures, circulation parameters and temperatures over west coast of India. Sci. Total Environ. 2016;551–552:175–185. doi: 10.1016/j.scitotenv.2016.02.016. [PubMed] [Cross Ref]
37. Levitus S, Antonov JI, Boyer TP, Stephens C. Warming of the World. Ocean. Science (80-.). 2000;287:2225–2229. doi: 10.1126/science.287.5461.2225. [Cross Ref]
38. Collatz GJ, et al. A mechanism for the influence of vegetation on the response of the diurnal temperature range to changing climate. Geophys. Res. Lett. 2000;27:3381–3384. doi: 10.1029/1999GL010947. [Cross Ref]
39. Tian H, Banger K, Bo T, Dadhwal VK. History of land use in India during 1880–2010: Large-scale land transformations reconstructed from satellite data and historical archives. Glob. Planet. Change. 2014;121:78–88. doi: 10.1016/j.gloplacha.2014.07.005. [Cross Ref]
40. New M, Hulme M, Jones P. Representing twentieth-century space-time climate variability. Part II: Development of 1901-96 monthly grids of terrestrial surface climate. J. Clim. 2000;13:2217–2238. doi: 10.1175/1520-0442(2000)013<2217:RTCSTC>2.0.CO;2. [Cross Ref]
41. Shepard, D. A two-dimensional interpolation function for irregularly-spaced data. in Proceedings of the 1968 23rd ACM national conference on - 517–524, doi:10.1145/800186.810616 (ACM Press, 1968).
42. Willmott, C. & Matsuura, K. T A Temperature and Precipitation: Monthly and Annual Time Series (1950–1999), at http://www.esrl.noaa.gov/psd/data/gridded/data.UDel_AirT_Precip.html (2001).
43. Deshpande, N. R., Kothawale, D. R. & Kulkarni, A. Changes in climate extremes over major river basins of India. Int. J. Climatol., doi:10.1002/joc.4651 (2016).
44. Weather and Climate Extremes., doi:10.1007/978-94-015-9265-9 (Springer Netherlands, 1999).
45. Zhang X, Hegerl G, Zwiers FW, Kenyon J. Avoiding Inhomogeneity in Percentile-Based Indices of Temperature Extremes. J. Clim. 2005;18:1641–1651. doi: 10.1175/JCLI3366.1. [Cross Ref]
46. Wang YH, Yeh CH, Young HWV, Hu K, Lo MT. On the computational complexity of the empirical mode decomposition algorithm. Phys. a-Statistical Mech. Its Appl. 2014;400:159–167. doi: 10.1016/j.physa.2014.01.020. [Cross Ref]
47. Wu Z, Huang NE, Long SR, Peng C-K. On the trend, detrending, and variability of nonlinear and nonstationary time series. Proc. Natl. Acad. Sci. 2007;104:14889–14894. doi: 10.1073/pnas.0701020104. [PubMed] [Cross Ref]
48. Wu Z, Huang NE. Ensemble Empirical Mode Decomposition: A noise-assisted data analysis method. Adv. Adapt. Data Anal. 2009;01:1–41. doi: 10.1142/S1793536909000047. [Cross Ref]
49. Wu Z, Huang NE, Wallace JM, Smoliak BV, Chen X. On the time-varying trend in global-mean surface temperature. Clim. Dyn. 2011;37:759–773. doi: 10.1007/s00382-011-1128-8. [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group