www.nature.com/scientificreports OPEN Historical and Projected Surface Temperature over India during the th st 20 and 21 century Received: 22 July 2016 1 2 1 1 3 Ghouse Basha , P. Kishore , M. Venkat Ratnam , A. Jayaraman , Amir Agha Kouchak , Taha B. 4,5 2 Accepted: 7 April 2017 M. J. Ouarda & Isabella Velicogna Published: xx xx xxxx Surface Temperature (ST) over India has increased by ~0.055 K/decade during 1860–2005 and follows the global warming trend. Here, the natural and external forcings (e.g., natural and anthropogenic) responsible for ST variability are studied from Coupled Model Inter-comparison phase 5 (CMIP5) th st models during the 20 century and projections during the 21 century along with seasonal variability. Greenhouse Gases (GHG) and Land Use (LU) are the major factors that gave rise to warming during the th 20 century. Anthropogenic Aerosols (AA) have slowed down the warming rate. The CMIP5 projection over India shows a sharp increase in ST under Representative Concentration Pathways (RCP) 8.5 where st it reaches a maximum of 5 K by the end of the 21 century. Under RCP2.6 emission scenarios, ST st increases up to the year 2050 and decreases afterwards. The seasonal variability of ST during the 21 century shows significant increase during summer. Analysis of rare heat and cold events for 2080–2099 relative to a base period of 1986–2006 under RCP8.5 scenarios reveals that both are likely to increase substantially. However, by controlling the regional AA and LU change in India, a reduction in further warming over India region might be achieved. Global Surface Temperature (ST) has increased significantly during the last three decades. The recent report of the Intergovernmental Panel on Climate Change (IPCC) showed the anthropogenic footprint on the global warming of atmosphere and oceans, reduction in snow and ice, and increase in sea level. Climate change and variability comprise complex interactions between natural and external forcings at different temporal and spatial scales. However, quantifying the contribution of each component in the overall observed climate signal is very challeng- ing. Climate models represent important tools to investigate the historical climate change due to anthropogenic forcings. With these models and controlled simulations, it is possible to decompose the contributions from indi- vidual factors or assess their combined effects on regional or global climate. This allows exploring the possible causes and understanding the physical mechanisms behind climate change and variability . The World Climate Research Program (WCRP) has supported the development of the Coupled Model Intercomparison Project (CMIP), which provides simulation from state-of-the-art global climate models. CMIP model simulations have been used in the Assessment Reports of the Intergovernmental Panel on Climate th Change (e.g., IPCC AR4 and AR5). The CMIP5 simulations include both historic (20 century) and future pro- st 2 jections (21 century) under different Representative Concentration Pathways (RCPs) . The RCPs represent different mitigation scenarios that influence the future emissions of greenhouse gasses, aerosols, ozone and land use change. The CMIP5 experiment simulates historical changes in temperature with forcings driven by Anthropogenic Aerosols (AA), Greenhouse Gases (GHGs), Land Use change (LU), Natural (NAT) forces (Solar Irradiance + Volcanic Activity), and Solar radiation (SL) . Comparing these forcings with each other provides a 4–8 better understanding of the contributors to climate change over a particular region or location . India is one of the fastest-growing economies in the world. The annual ST has increased since the 1950’s, 9–11 particularly over India. Although, previous studies have shown a clear increasing trend in mean, minimum, and maximum temperatures over India during the last century, the contributing factors to this have not been 1 2 National Atmospheric Research Laboratory, Gadanki, Tirupati, India. Department of Earth System Science, University of California, Irvine, California, 92697, USA. Department of Civil and Environmental Engineering, University of California, Irvine, California, 92697, USA. Institute Center for Water and Environment (iWATER), Masdar Institute of Science and Technology, P.O. Box 54224, Abu Dhabi, UAE. INRS-ETE, National Institute of Scientific Research, Quebec City (QC), G1K9A9, Canada. Correspondence and requests for materials should be addressed to G.B. (email: firstname.lastname@example.org) Scientific Repo R ts | 7: 2987 | DOI:10.1038/s41598-017-02130-3 1 www.nature.com/scientificreports/ Figure 1. Taylor diagram for mean surface temperature between CRU and CMIP5 individual model during the period 1901–2005 (a), IMD and CMIP5 individual model for the years 1969–2005 (b) over India. Solid lines correspond to standard deviation; dotted lines for correlations (Figure was created using the Interactive Data Language (IDL) version 8.2 software http://www.harrisgeospatial.com/docs/whats_new_in_82.html). explored in detail. The observed climate change in India is due to the complex interactions between natural and anthropogenic forcings at different temporal and spatial scales. However, it is not well understood to what extent climate change in India during the past century, particularly the past several decades, can be attributed to human and natural forcings. In this study, we have utilized the CMIP5 simulation from the years 1860–2100 along with ground-based observational data sets (Climate Research Unit, CRU and India Meteorological Department, IMD gridded temperature) to quantify the historical (Natural + Anthropogenic), NAT, and external forcings (AA, GHGs, LU and SL) in the observed warming over India. We evaluate the 17 CMIP5 historical ST data with obser- vational data sets along with future projections. We investigate the contribution of different natural and anthro- pogenic forcings in the observed and expected future warming. We also evaluate trends in different simulations during different periods between the years 2006–2100. In addition, we discuss changes in ST extremes over India. Results Evaluation of CMIP5 ST with CRU and IMD. The ST simulation from the individual CMIP5 model is evaluated against the CRU and IMD data sets as shown in Figures S1 and S2, respectively. These simulations are compared for different periods: 1901–2005 (against CRU data) and 1969–2005 (against IMD data). From both figures, it is clear that the North Eastern part and Himalayan regions show a large bias for all models, possibly due to complexities in higher elevations also mentioned in other studies . Further, the overall bias is even larger when CIMP5 simulations are compared with IMD relative to the CRU data set, mainly on the western coast of India. The models CSIRO-Mk3-6-0, GISS-E2-H, GISS-E2-R, HadGEM2-CC, HadGEM2-ES, MIRCO-ESM, MRI-CGCM3, and inmcm4 show larger biases compared to other models. For further evaluation, we show the Taylor diagrams for representing the simulated mean STs relative to CRU (Fig. 1a) and IMD (Fig. 1b). The former covers 1901–2005, whereas the latter covers 1969–2005 over India. In these figures, the solid lines correspond to standard deviations, while the dotted lines represent the corre- lations between simulations and observations. As shown, the correlation between 8 CMIP5 models (CCSM4, CSIRO-Mk3-6-0, GISS-E2-H, GISS-E2-R, HadGEM2-CC, HadGEM2-ES, GFDL-CM3, and MICRO-ESM) with CRU and IMD data sets is less than 0.5. To continue the analysis, we have discarded these models, and considered the ones that show high correlations (i.e., CNRM-CM5, CanESM2, GFDL-CM3, IPSL-CM5A-LR, MIRCOC5, MPI-ESM-LR, NorESM1-M, and bcc-csm1-1). The Taylor diagram presents the spread of the models in terms of their spatial correlation. The correlations between models and the two observational data sets are quite similar (Fig. 1a and b). The highest correlation score of about of 0.78 is observed for MIROC5 and IPSL-CM5-AL. The climatological mean biases between CMIP5 simulations and the two observational data sets (i.e., IMD, CR) are shown in Fig. 2. Note that this bias and trends are estimated for the period of 1969–2005 when the two observational datasets overlap. The mean bias between IMD and CRU is less than ± 1 K, as shown in Fig. 2a, except for the Northeastern part and the Himalayan regions. The climatological bias between IMD and the ensemble mean of the selected 8 CMIP5 models is shown in Fig. 2b. The bias over India varies around ± 1 K, except over the Northeastern and Himalayan regions where the bias exceeds 4 K. Similar bias is observed between CRU and CMIP5 models. Higher bias in the Northeastern and Himalayan regions may be due to the fact that there are few observations available in these regions. In addition, we computed the spatial trends based on the historical climate model simulations and observations (IMD, CRU) - see Fig. 2(d–f ). In Fig. 2(d–f ), the solid black mark indicates the 95% confidence level, as determined at each individual grid point and was computed by using student’s t-test. Over India, there is a warming trend at the rate of ~0.2 K/decade from the year 1969–2005. Few parts (Northeastern and Himalayan regions) of India exhibit decreasing trends in IMD data (Fig. 2d). Recent studies also indicate relatively suppressed warming over India, with some cooling over the central regions, during the past century . This feature is not shown in the CRU data and CMIP5 simulations. The central part of India shows the highest warming compared to other regions as shown by both CRU and CMIP5 simulations (~0.36 K/ Scientific Repo R ts | 7: 2987 | DOI:10.1038/s41598-017-02130-3 2 www.nature.com/scientificreports/ Figure 2. Spatial temperature difference between (a) IMD and CRU, (b) IMD and ensemble mean of CMIP5 models, and (c) CRU and ensemble mean CMIP5 models. Spatial surface temperature trends over India for (d) IMD, (e) CRU and Ensemble mean of CMIP5 models for the years 1969–2005 with 95% confidence level (hatched regions) (Figure was created using the Interactive Data Language (IDL) version 8.2 software http:// www.harrisgeospatial.com/docs/whats_new_in_82.html). Seasonal trends (K/decade) Winter Spring Summer Fall All years Name 1860–2005 Historical 0.056 0.053 0.071 0.040 0.055 Historical_AA −0.029 −0.025 −0.019 −0.033 −0.027 Historical_GHG 0.121 0.092 0.115 0.099 0.106 Historical_LU 0.086 0.054 0.079 0.074 0.073 Historical_NAT 0.019 0.025 0.012 0.004 0.006 Historical_SL −0.014 −0.004 −0.015 0.005 −0.007 1901–2005 Historical 0.077 0.069 0.100 0.489 0.074 CRU 0.084 0.022 0.064 0.072 0.060 1969–2005 Historical 0.225 0.183 0.217 0.199 0.207 IMD 0.137 0.112 0.068 0.089 0.081 CRU 0.252 0.116 0.099 0.202 0.168 Table 1. Historical, CRU and IMD seasonal trends at 95% confidence level over the Indian region. decade). Overall, there are discrepancies in the warming trends. However, CRU observations and CMIP5 simula- tions seem to be more consistent with one another compared to IMD data. We have also performed the seasonal trend analysis in ST over India as shown in Table 1. Scientific Repo R ts | 7: 2987 | DOI:10.1038/s41598-017-02130-3 3 www.nature.com/scientificreports/ Figure 3. Temporal variation of annual surface temperature anomalies during the years 1860–2005. ‘All’ indicates the historical forcing, which contains all types of forcings. Individual forcings include Anthropogenic Aerosol (AA), Greenhouse Gases (GHG), Land Use change (LU), Natural forcings (solar radiation + volcanic eruptions), and Solar radiation (SL). In each forcing the observed CRU temperature data is also plotted for the years 1901–2005. Robust regression analysis was performed at the 95% confidence level for trend analysis which is over plotted in each panel (Figure was created using the Interactive Data Language (IDL) version 8.2 software http://www.harrisgeospatial.com/docs/whats_new_in_82.html). th Impact of different forcings on ST. The ensemble means of different historical forcings for the 19 and th 20 centuries (shown in Table S1) are plotted in Fig. 3. e Th forcing responses are computed and compared with the observed ST based on ground-based data. Figure 3 show temporal variation of mean ST over India from CRU, and historical simulation (GHG, NAT, AA, LU, and SL forcings) data. These temporal variations are calculated as the mean of all years over India (using data from Jan. 1860 to Dec. 2005) subtracted from the yearly ST values of each simulation. In each panel, the ensemble bounds (maximum and minimum values) are shown in shaded color. The results of the robust regression fitting, along with CRU ST are shown in Fig. 3 for each forcing. The increase in temperature due to GHGs forcing from 1900 onwards is clearly observed in Fig. 3b. The turning point (changing from negative to positive temperature anomaly) in the yearly time series of ST variation matches very well with the historical ensemble mean of the different forcings. The observed CRU ST is superimposed on the ensemble of historical model simulations (Fig. 3a). Both historical simulations and CRU ST show similar vari- th ability throughout the 20 century (Fig. 3a). Except for NAT, SL and AA, all other forcings show a remarkable increase in ST, some even higher than the observations (e.g., 3c and 3d). This shows that considering LU and GHG alone leads to an overestimation of the surface temperature warming in the historical period. However, increase in ST from the combination of all forcings closely follows the observed historical temperatures. The warming due to GHG shows a sharp increase after1960 over India, which is very much consistent with LU forcings. The combination of both direct and indirect AA forcings depicts the significant decreasing from 1960 onwards. The impacts of NAT and SL forcings on warming in ST are comparatively much smaller than those from GHG, and LU. These simulations with individual forcings are useful to compare the role of different forcings. Individually, they are not expected to reproduce historical observations. However, we have included the historical observations as a common reference for evaluation. Overall, the results show that AA, NAT and SL provided cooling effect and slow the warming rate, while the other forcings (GHG and LU) contributed significantly to the observed warming. The results show that AA has negative contribution to the warming when we considered both the direct and indirect forcings. When direct forcings alone is consider, there is a weak cooling trend. However, with both direct and indirect AA the cooling trend becomes highly significant. Under GHG forcings, only the warming Scientific Repo R ts | 7: 2987 | DOI:10.1038/s41598-017-02130-3 4 www.nature.com/scientificreports/ Figure 4. (a) Box-and-whisker plots for temperature of different forcings along with observed CRU data sets from 1901–2005. The boxes indicate the interquartile model spread (range between the 25th and 75th quantiles), the black/grey solid marks within the boxes show the multimodel median and the whiskers indicate the total intermodel range. (b) Probability density distribution of surface temperature over the Indian region relative to the years 1910–1940 for different forcings (Figure was created using the Interactive Data Language (IDL) version 8.2 software http://www.harrisgeospatial.com/docs/whats_new_in_82.html). trend of 0.16 K/decade is observed which is very strong compared to all other forcings. The seasonal variations of different forcings are depicted in Table 1. The largest seasonal variations in GHG and LU are observed during winter followed by summer, fall and spring. GHG variability is larger during winter and summer compared to the other seasons. The observed increasing trend is larger for GHG compared to LU. In case of AA forcings, the seasonal variation of decreasing trend is large during fall and minimum during summer. A small positive trend is noticed from NAT forcing in all the seasons. The SL radiation shows a significant negative trend, except during the fall season. Figure 4 shows the box and whisker plot along with probability density functions of ST from different forcings th th relative to the base period of 1910–1940. The boxes indicate the interquartile model spread (i.e., 25 and 75 th quartiles). The 75 quartile of the GHG forcing is highest among others and reaches 298.2 K. Figure 4 shows that the difference between GHG and NAT forcings is around ~3 K, which suggests the role of GHG alone relative to the natural variability. The highest ST values are found in GHGs followed by LU, AA, SL and NAT, respectively. e p Th robability density functions of the historical temperature anomalies (relative to the period 1910–1940) show clear warming shifts of about ~0.15 K in GHGs, followed by LU. The AA forcings shows the negative shift in the probability distribution. The distributions of the remaining forcings fall between 0–0.2 K temperature changes in the historical period . Showed the competing roles of GHG and AA with regard to Indian Ocean warming. e Th spatial distribution of ST trends over India is estimated using 146 years of simulations (1860–2005) with different forcings (Fig. 5). Here, the ST trends are computed at each grid point using the robust regression tech- nique . In this process, we have tested the statistical significance at each grid cell and only considered trends that are significant at the 0.05 significance level (95% confidence) marked with ‘x’ symbol. Figure 5a shows the warming response over India at a 0.18 K/decade rate based on all forcings. The NAT forcings show a slight warm- ing trend over central India and cooling trends in other parts, which, as expected, is not consistent with the Scientific Repo R ts | 7: 2987 | DOI:10.1038/s41598-017-02130-3 5 www.nature.com/scientificreports/ Figure 5. Geographical distribution of surface temperature trends (K/decade) for historical data along with different forcings. (a) Historical surface temperature, (b) AA, (c) GHG, (d) LU, (e) NAT, (f) SL. Trends were estimated at the 95% confidence level (hatched regions) for each grid (Figure was created using the Interactive Data Language (IDL) version 8.2 software http://www.harrisgeospatial.com/docs/whats_new_in_82.html). observations. The response of GHGs and LU forcing is very large and a warming trend is observed throughout the region. The warming (cooling) rate response to GHG (AA) forcing in North India is large compared to the South India (Fig. 5b and c). Sensitivity experiments, such as correlation analysis with all forcings, are also performed in the present work. e co Th rrelation values are presented in Table S2. The observed ST from CRU shows good correlation with GHG forcings while historic data shows a very high correlation with GHG and LU. A negative correlation is detected between AA and SL with observed data. GHG emissions contributed to a warming of about 0.106 K/decade over the period 1860 to 2005 whereas the contributions LU to the warming trend are around 0.073 K/decade. The contribution from natural forcings is about 0.006 K/decade, very much less than external forcings. Solar radiation shows a negative trend, which means a cooling effect on ST. The contribution of natural and external forcings contributed to a mean surface warming of about 0.055 K/decade. It is extremely likely that more than half of the observed increase in average ST from 1950 onwards was caused by external forcings. st Future projections. For projecting climate change impacts on ST during the 21 century, we considered low and high RCPs (i.e., RCP2.6 and RCP8.5) emission scenarios. The evolution of ST changes relative to 1901–1960 over India from IMD, CRU, historic simulations, and ensemble mean of 8 CMIP5 simulations are shown in Fig. 6. RCP2.6 and RCP8.5 show around 3.2 degrees temperature difference by the end of the century. During the last century, the change in surface temperature is well captured by all data sets. Temperature projections show consist- ent increasing trends in RCP8.5 whereas the projected mean temperature from RCP2.6 increases up to 2050, and then slightly decreases, consistent with the radiative forcing under RCP2.6. The mean temperature over India is st st expected to increase over the 21 century under all the RCPs. Around the middle of 21 century, the rate of warm- ing increases sharply depending on the scenario (e.g., RCP 8.5). On average, the projected change in temperature st by the end of the 21 century relative to 1901–1960 is expected to be 1.8 K, 2 K, 3.5 K, and 5 K for RCP2.6, RCP4.5, RCP6.0 and RCP8.5 emission scenarios, respectively. st Figure 7 illustrates the spatial pattern of ST trends over India during the 21 century under RCP2.6 and RCP8.5 emission scenarios. The ST trends were estimated for different periods such as near future (2006–2035), mid-future (2046–2065), and far future (2080–2099) as well as for the entire time series of model simulations (2006–2099). Under RCP2.6, the temperature increases in near future and decreases in both mid- and far-future, but the entire time series shows an overall increase in ST of about ~0.15 K/decade. The far-future projections show Scientific Repo R ts | 7: 2987 | DOI:10.1038/s41598-017-02130-3 6 www.nature.com/scientificreports/ Figure 6. Projected annual mean surface temperature from multimodel average CMIP5 data for the years 1860–2100 relative to the base period 1901–1960. The period 2006–2099 represents the future projection scenarios from different RCPs. The shaded region represents the standard deviation of the ensemble mean of CMIP5 models. The observed data sets CRU (1901–2005) and IMD (1969–2005) are also plotted (Figure was created using the Interactive Data Language (IDL) version 8.2 software http://www.harrisgeospatial.com/docs/ whats_new_in_82.html). Figure 7. Geographical distribution of surface temperature trends (K/decade) at the 95% confidence level st (hatched regions) from RCP2.6 and RCP 8.5 for different periods over India for the 21 century. Different periods correspond to near (2016–2035), mid (2046–2065) and far (2080–2099) and for the whole period (2006–2099) (Figure was created using the Interactive Data Language (IDL) version 8.2 software http://www. harrisgeospatial.com/docs/whats_new_in_82.html). cooling all over India except for southern India (Fig. 7c). Continuous increases in ST are expected in near, mid-, st and far future under the RCP8.5 scenario. Significantly, higher temperatures are expected in the late 21 century compared to near and middle future temperatures. The long-term (2006–2099) trends show warming of about st 0.2 K/decade over India. In addition, the warming increases rapidly toward the end of the 21 century, reaching a peak value of about 0.4 K/decade. We have also analyzed the seasonal variations of ST trends over India from RCP2.6 and RCP8.5 as shown in Table S3 for different time periods. In the near future, under RCP2.6 the temperature increases sharply and is highest during the summer season, while ST decreases in mid- and far-future. Under RCP8.5, the ST increases during all periods. The far-future ST trend increases drastically (nearly 3 times the near and mid-future trends). The winter ST shows significant increase followed by spring over India. The entire time series indicates that ST st can increase by ~4.5 K by the end of the 21 century. Extreme Temperatures. Recent studies have shown the variability and mechanisms for occurrence of heat 15, 16 waves over India . Increasing trends in frequency, total duration and maximum duration of heat waves are Scientific Repo R ts | 7: 2987 | DOI:10.1038/s41598-017-02130-3 7 www.nature.com/scientificreports/ observed over the central and northwestern parts of India . showed that the increase in heat wave activity is due to global warming over India. Here, the extremely high and low monthly temperatures (i.e. rare heat and rare cold events) from CMIP5 simulations are analyzed. The 20-year difference value of the annual maximum or minimum 17, 18 temperature is a measure of changes in rare temperature extremes . Figure S3 indicates the projected change st in the 20-yr return value of the annual maximum and minimum daily ST over India by the end of the 21 century (2081–2100), relative to the recent past (1986–2005) from RCP2.6 and RCP8.5 emission scenarios. The annual st maximum temperature extremes under RCP2.6 are projected to decrease by the end of the 21 century, except for the Eastern part of central India and the Himalayan region. Under RCP8.5, the Himalayan region is pro- st jected to face higher extreme temperatures by the end of the 21 century. Moderate increases are expected in the Northeastern, Central and Northwestern parts of India (Figure S3 top right panel). Figure S3 (bottom) shows the minimum temperature extremes (rare cold events) under RCP2.6 and RCP8.6 emission scenarios. The projected changes in rare cold events are higher under RCP8.5 than RCP 2.6. The northern part of India can expect larger changes than the other parts of the country. Summary and Conclusions The surface temperature (ST) changes over India are studied by using 109 simulations from global coupled cli - th mate models with different external forcings during the 20 and future projections under different RCP emission st scenarios in the 21 century. Based on the CMIP5 multi-model simulations, the relative contribution of external forcings such as natural (solar radiation, volcanic emissions) and human activities (GHGs, AA, and LU) to the observed ST during the last 100 years were computed. In addition, we have also considered the individual and combined effects of direct and indirect AA forcings. The evaluation of CMIP5 historic data with observational st data sets and future projections during the 21 century under different RCP scenarios is also presented. The long-term trends were estimated by using robust regression analysis. The main findings of the present study are summarized below: th • During the 20 century, the major contributors for the increase in ST are GHGs and LU over India whereas the AA has slowed the warming rate over India. The NAT and SL forcings shows slight negative warm- ing trend. The warming rate from GHGs and Lu are estimated to be 0.14 ± 0.53 and 0.06 ± 0.2 K/centur y, respectively. • The AA forcings show cooling effect when we considered both the direct and indirect forcings. The cooling effect is weekend when only direct effect of AA is considered. • Human activities are highly responsible for the increase in warming over India. The warming is predomi- nantly attributed to GHGs followed by LU. A sharp increase in warming is noticed from the 1960s onward 19, 20 from GHGs, whereas for LU the increase is observed from 1980 . • The spatial variation of trends from die ff rent forcings over India indicates that the impact of GHGs is more severe over the Northern and Western parts of India. th • The trends estimated during the 20 century show highest warming in historic data followed by CRU and IMD data. The central part of India warmed severely compared to other regions. The observed warming is more pronounced during the summer followed by winter, monsoon season, and post-monsoon. • The multi-model mean of CMIP5 captures the basic features of observed ST. The future protections under different RCP emission scenarios are very distinct. Rapid increase in warming is noticed under RCP8.5 where st it reaches a maximum temperature of ~5 K by the end of the 21 century. In RCP2.6 emission scenarios the warming increases until 2050 and decreases aer ft wards. • Future warming is larger during the winter season followed by summer, monsoon and post-monsoon seasons st over India. By the end of the 21 century, both rare cold and heat events are expected to increase over India under RCP 8.5 emission scenario. GHGs (water vapor, CO , CH , N O) act to make the surface warmer by absorbing and emitting heat energy 2 4 2 in all the directions. Addition of more GHGs to the atmosphere make more effective in preventing the heat from escaping into space. The GHG’s concentration increased significantly since the industrial revolution began. The atmospheric concentrations of CO , CH , N O has increased by 40%, 150% and 20%, respectively. This increase 2 4 2 in GHGs is due to emissions by human activities, which alter the Earth’s energy balance resulting in significant increase in ST thus global warming. The anthropogenic forcings plays a significant role in climate change hence global warming. Humans alter aerosols not only at the surface (e.g. industrial emissions) but also at the higher altitudes. Due to increase in populations, the natural vegetation decreases for thousands of years by deforestation. th e Th largest changes occurred in 20 century. Hence changing land cover ae ff cts the climate by modifying the sur - face reflectivity and the hydrological cycle. In addition to the direct radiative effects, LU can ae ff ct the sources and sinks of GHG’s and the amount of dust lifted into the atmosphere by the wind. For example, atmospheric CO , CH , and N O concentrations are increased by deforestation, cultivation of rice and by use of nitrogen fertilizers 4 2 in agriculture activity. Natural variability of aerosols, particularly due to volcanoes eruptions (e.g. Mt. Pinatubo in 1991) are recognized as a significant climate change indication, which alters the Earth’s radiation balance and thus tends to cause the increase in temperature. These and other lines of evidence point conclusively to the fact that the elevated CO concentration in our atmosphere is the result of human activities. It is well known that net global warming is very sensitive to the region of focus. Dier ff ent local and regional topographical and meteor - ological factors play an important role in determining global warming and climate change. India will observe changes in precipitation, cyclone occurrence, human health, agriculture, extreme events . Although a sudden increase in ST was observed during the last century, the annual precipitation over India showed only moder- 22, 23 ate increase. However, the extreme precipitation frequency and magnitude increased rapidly . Nevertheless, newly introduced forcing agents such as historical emissions of carbonaceous aerosols should be included in Scientific Repo R ts | 7: 2987 | DOI:10.1038/s41598-017-02130-3 8 www.nature.com/scientificreports/ future studies for a better understanding of global warming impacts. Uncertainties related to future warming in model analysis also need further investigation. CMIP5 coordinated climate models experiments had leads to better understanding of past, present, and future climate change and variability. However, the science gaps and outstanding questions have promoted to design new experiments known as CMIP6, which address the science aspects such as Earth system respond to forcings, sources and consequences of systematic biases and assessment of future climate changes given climate variability, predictability and uncertainties in scenarios. In recent years, the asymmetry of the forcing has been noticed and attracted more and more attention of scientists. For example, the asymmetric volcanic eruptions at different latitudes can ae ff ct the global climate differently . In CMIP6, the latitude-dependent GHG concentration will also considered, that should have different effect on the Indian sur - face temperature compared to CMIP5. Data and Methodology Historical and Future projections from CMIP5 models. In this article, we have analyzed 109 simula- tions from 17 CMIP5 models, which include different emission scenarios, and historical simulations as shown in Table S1. Note that not all CMIP5 models simulate all the forcings. 17 CMIP5 models under different RCP 2.6, 4.5, 6.5, 8.5 emission scenarios are used along with NAT forcings from 12 simulations, GHGs forcings from 15 simulations, AA forcings from 5 simulations, LU forcings from 5 simulations and SL forcings from 4 simulations. e Th period 1860–2005, which is designated as the historic period, has been chosen for coordinated climate model experiments under CMIP5 . In this simulation, varying atmospheric composition (e.g., CO ) resulting from solar forcing, anthropogenic, volcanic, emissions of short-lived species and natural and anthropogenic aerosols are forced to reproduce the historical climate. These experiments are part of the CMIP5 coordinated experiments, used in the fifth Assessment Report by IPCC . In this study, we used multiple ensembles from 17 CMIP5 models with monthly mean ST outputs available for all four scenarios (RCP2.6, RCP4.5, RCP6.0, and RCP8.5) , and his- torical simulations (NAT, AA, GHGs, LU, and SL) with different forcing (see Table S1). The RCP4.5 corresponds −2 to radiative forcing of 4.5 Wm aer t ft he year 2100, which means that the CO emissions may exceed 650 ppm. −2 Similarly, the RCP8.5 corresponds to 8.5 Wm aer t ft he year 2100, which means that the equivalent CO exceeds 1370 ppm. Observational data sets. CRU data. The monthly gridded ST data from Climate Research Unit (CRU) version 3.22 is used in the present study for the years 1901–2005. This data is available with 0.5° × 0.5° latitude/ longitude grid over land regions only over the whole globe. This data is produced by CRU at the University of East Anglia making use of more than 4000 weather station data globally . IMD gridded data. e Th India Meteorological Department (IMD) developed gridded ST data at 1° × 1° resolu- tion . The observed data from 395 stations were selected for the development of gridded data. The observed data was subjected to quality checks such as rejecting values, which exceed known extreme levels, maximum temper- ature less than minimum temperature, same temperature values for many consecutive days etc. The 395 station data were interpolated into grids with the modified version of Shepard angular distance weighting algorithm . e d Th ata period spans between the years 1969–2005 and it is cross validated by estimating the root mean square error. The errors associated with the interpolation scheme were used in preparing gridded data over the plains, which was found to be 0.5 °C at the maximum. Methodology. The ST data from different forcings and all the CMIP5 model data sets under various RCPs are brought to a common grid of 1° × 1° longitude and latitude by using cubic spline interpolation. For trend estima- tion, the robust regression technique is used. This method is an important tool for analyzing the data ae ff cted by outliers and it is based on Iteratively Reweighted Least Squares Regression (IRLS) . Robust regression is used to 30, 31 detect and provide results that are resistant to outliers . The t-test analysis is employed to calculate the statis- tical significance of the temperature trends. In this process, we have estimated the confidence level at each trend grid point and only those spectral peaks with greater than 95% confidence level are considered. References 1. Pachauri, R. K. et al. Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC, 2014). 2. Taylor, K. E., Stouffer, R. J. & Meehl, G. A. An overview of CMIP5 and the experiment design. Bull. Amer. Meteor. Soc. 93, 485–498, doi:10.1175/BAMS-D-11-00094.1 (2012). 3. Santer, B. D. et al. Identifying human influences on atmospheric temperature. P. Natl. Acad. Sci. USA 110, 26–33, doi:10.1073/ pnas.1210514109 (2013). 4. Christidis, N., Stott, P. A., Zwiers, F. W., Shiogama, H. & Nozawa, T. The contribution of anthropogenic forcings to regional changes in temperature during the last decade. Clim. Dyn. 39(6), 1259–1274, doi:10.1007/s00382-011-1184-0 (2011). 5. Polson, D., Hegerl, G. C., Allan, R. P. & Sarojini, B. Have greenhouse gases intensified the contrast between wet and dry Regions? Geophys Res Lett. 40, 4783–4787, doi:10.1002/grl.50923 (2013). 6. Wilcox, L. J., Highwood, E. J. & Dunstone, N. J. The influence of anthropogenic aerosol on multi-decadal variations of historical global climate. Environ. Res. Lett. 8, 024033, doi:10.1088/1748-9326/8/2/024033 (2013). 7. Wu, P., Christidis, N. & Stott, P. Anthropogenic impact on Earth’s hydrological cycle. Nature Clim. Change, doi:10.1038/nclimate1932 (2013, 34). 8. Cheng, L., AghaKouchak, A. & Phillips, T. Non-stationary Return Levels of CMIP5 Multi-Model Temperature Extremes. Clim. Dyn. 44, 2947–2963, doi:10.1007/s00382-015-2625-y (2015). 9. Kothawale, D. R. & Kumar, K. R. On the recent changes in surface temperature trends over India. Geophys. Res. Lett. 32, L18714–n/a, doi:10.1029/2005GL023528 (2005). 10. Kothawale, D. R., Munot, A. A. & Krishna Kumar, K. Surface air temperature variability over India during 1901–2007, and its association with ENSO. Climate Res. 42, 89–104, doi:10.3354/cr00857 (2010). 11. Murari, K. K., Ghosh, S., Patwardhan, A., Daly, E. & Salvi, K. Intensification of future severe heat waves in India and their effect on heat stress and mortality. Reg. Environ Change. 15, 569–579, doi:10.1007/s10113-014-0660-6 (2015). Scientific Repo R ts | 7: 2987 | DOI:10.1038/s41598-017-02130-3 9 www.nature.com/scientificreports/ 12. Liu, Z., Mehran, A., Phillips, T. J. & AghaKouchak, A. Seasonal and regional biases in CMIP5 precipitation simulations. Clim Res. 60, 35–50, 10.3354/cr01221 (2014). 13. Roxy, M. K. et al. Drying of Indian subcontinent by rapid Indian Ocean warming and a weakening land-sea thermal gradient. Nature Comm. 6, 7423, doi:10.1038/ncomms8423 (2015). 14. Dong, L. & Zhou, T. The Indian Ocean sea surface temperature warming simulated by CMIP5 models during the twentieth century: Competing forcing roles of GHGs and anthropogenic aerosols. Jour. of Clim. 27(9), 3348–3362, doi:10.1175/JCLI-D-13-00396.1 (2014). 15. Ratnam, J. V. et al. Anatomy of Indian heat waves. Sci. Rep. 6, doi:10.1038/srep24395 (2016). 16. Rohini, P., Rajeevan, M. & Srivastava, A. K. On the Variability and Increasing Trends of Heat Waves over India. Sci. Rep. 6, 26153, doi:10.1038/srep26153 (2016). 17. Cheng, L., AghaKouchak, A., Gilleland, E. & Katz, R. W. Non-stationary Extreme Value Analysis in a Changing Climate. Clim. Change. 127(2), 353–369, doi:10.1007/s10584-014-1254-5 (2014). 18. Wuebbles, D. et al. CMIP5 climate model analyses: climate extremes in the United States. Bull. Am. Meteorol. Soc. 95, 571–583, doi:10.1175/BAMS-D-12-00172.1 (2014). 19. Tett, S. F. B., Stott, P. A., Allen, M. R., Ingram, W. J. & Mitchell, J. F. B. Causes of twentieth-century temperature change near the Earth’s surface. Nature 399, 569–572, doi:10.1038/21164 (1999). 20. Tett, S. F. B. et al. Estimation of natural and anthropogenic contributions to 20th century temperature change. J. Geophys. Res. 107, doi:10.1029/2000JD000028 (2002). 21. Chandrappa, R., Kurshrestha, U. C. & Gupta, S. Chapter 17: Changing Livelihood Due to Climate Change: Asian Perspective. In coping with Climate Change (pp. 245–269). Berlin, Heidelberg: Springer Berlin Heidelberg, doi:10.1007/978-3-642-19674-4 (2011). 22. Goswami, B. N., Venugopal, V., Sengupta, S., Madhusoodanan, M. S. & Xavier, P. K. Increasing trend of extreme rain events over India in a warming environment. Science 314, 1442–1444, doi:10.1126/science.1132027 (2006). 23. Kishore, P. et al. Precipitation climatology over India: validation with observations and reanalysis datasets and spatial trends. Clim. Dyn, doi:10.1007/s00382-015-2597-y (2015). 24. Haywood, J. M., Jones, A., Bellouin, N. & Stephenson, D. Asymmetric forcing from stratospheric aerosols impacts Sahelian rainfall. Nature Clim. Change 3, 660–665, doi:10.1038/nclimate1857 (2013). 25. Meehl, G. A. et al. The WCRP CMIP3 multimodel dataset. A new era in climate change research. Bull. Am. Meteorol. Soc. 88, 1383–1394, doi:10.1175/BAMS-88-9-1383 (2007). 26. Jones, P. D. CRU TS3.22: Climatic Research Unit (CRU) Time-Series (TS) Version 3.22 of High Resolution Gridded Data of Month- by-month Variation in Climate (Jan. 1901- Dec. 2013). NCAS British Atmospheric Data Centre, 24 September 2014. doi:10.5285/18BE23F8-D252-482D-8AF9-5D6A2D40990C (2014). 27. Srivastava, A. K., Rajeevan, M. & Kshirsagar, S. R. Development of a high resolution daily gridded temperature data set (1969–2005) for the Indian region. Atmos. Science Lett. 6, n/a–n/a, doi:10.1002/asl.232 (2009). 28. Shepard, D. A two-dimensional function for irregularly spaced data. In 23rd ACM National Conference, pp. 517–524 (1968). 29. Holland, P. E. & Welsch, R. E. Robust regression using iteratively reweighted least-squares. Commun Stat A6, 813–827 (1997). 30. Kutner, M., Nachsteim, C. J. & Neter, J. Applied Linear Regression Models 4th Ed. McGraw-Hill, New York, 2004. 31. Kishore, P. et al. Evaluating CMIP5 models using GPS radio occultation COSMIC temperature in UTLS region during 2006–2013, twenty-first century projection and trends. Clim. Dyn. 47, 3253–3270, doi:10.1007/s00382-016-3024-8 (2016). Acknowledgements We acknowledge the GCM modeling groups, the Program for Climate Model Diagnosis and Inter-comparison (PCMDI), and the WCRP’s Working Group on Coupled Modeling for their roles in making available the WCRP CMIP5 multi-model datasets. The authors would like to thank the National Center for Atmospheric Research (NCAR) for providing the CRU data set. The authors are thankful to the National Climate Centre, India Meteorological Department for providing gridded data. All the figures were created using Interactive Data Language (IDL) version 8.2 (http://www.harrisgeospatial.com/docs/whats_new_in_82.html). The authors would like to thank three anonymous reviewers whose comments helped considerably in improving the quality of this paper. Author Contributions G.B. conceived and designed the idea. P.K. carried out the analysis. M.V.R., A.J. A.A., T.B.M.J.O. and I.V. provided the guidance for the analysis and all the authors have contributed to writing the paper. Additional Information Supplementary information accompanies this paper at doi:10.1038/s41598-017-02130-3 Competing Interests: The authors declare that they have no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2017 Scientific Repo R ts | 7: 2987 | DOI:10.1038/s41598-017-02130-3 10
Scientific Reports – Springer Journals
Published: Jun 7, 2017
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.
Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.
All the latest content is available, no embargo periods.
“Hi guys, I cannot tell you how much I love this resource. Incredible. I really believe you've hit the nail on the head with this site in regards to solving the research-purchase issue.”Daniel C.
“Whoa! It’s like Spotify but for academic articles.”@Phil_Robichaud
“I must say, @deepdyve is a fabulous solution to the independent researcher's problem of #access to #information.”@deepthiw
“My last article couldn't be possible without the platform @deepdyve that makes journal papers cheaper.”@JoseServera