Projections of heat stress and associated work performance over India in response to global warming

Summertime heat stress future projections from multi-model mean of 18 CMIP5 models show unprecedented increasing levels in the RCP 4.5 and RCP 8.5 emission scenarios over India. The estimated heat stress is found to have more impact on the coastal areas of India having exposure to more frequent days of extreme caution to danger category along with the increased probability of occurrence. The explicit amount of change in temperature, increase in the duration and intensity of warm days along with the modulation in large scale circulation in future are seemingly connected to the increasing levels of heat stress over India. A decline of 30 to 40% in the work performance is projected over India by the end of the century due to the elevated heat stress levels which pose great challenges to the country policy makers to design the safety mechanisms and to protect people working under continuous extreme hot weather conditions.

www.nature.com/scientificreports/ risks and adversely impacts the work productivity 20,35 . Study of 370 onsite data sets in Hong Kong reported that the work duration of the construction labour has been decreased by 0.33% with one degree rise in WBGT 36 . Analysis of 8076 workers of 11 different studies revealed the 30% loses in work productivity due to rise in WBGT in the people who are working in single shift 37 . The loss in labour productivity due to heat stress over Australia has been reported 38 and a decrease in productivity globally upto 20%. In tropical lands, the agricultural activities of working labor which involves higher levels of physical exertion had been badly affected by the heat stress 39 . It is reported that the estimated labour productivity in relation to heat stress index for 21 regions of the world infers a reduction in labour productivity due to warming of climate and resultant heat stress 40 . Also for the workers with low income in the tropical regions, heat stress is going to be the greatest health hazard where the facilities such as health surveillance are not available 41 . Venugopal 42 studied the health and work productivity issues of the workers from 18 places in India and found that they are vulnerable to heat stress. The impact of heat stress on working people in hot environments is measured using WBGT by ISO 7243 standard 43 . Also, a few reports infer that the work ability estimated from the metrics of WBGT over New Delhi, India show better results compared to the other metrics such as Predictive Heat Stress, Universal Thermal Climate Index 44 . However, the studies of Brode et al. 44 applied limited criteria in simulating the work ability during the day time working under light, moderate and heavy conditions and found the work ability change with "work load, time of the day and climate type". Also, the results of the Brode et al. 44 were suggested to be representatives of prevailing in agriculture,health and industrial sectors. The present unusual warming is expected to increase faster and becomes vulnerable to physiological acclimatization among the workers carrying out activities in outdoor and indoor work places without thermal neutralities. Report from International Labour Organization (ILO) says the decline in global work productivity will be same as the loss of 80 million full time jobs by 2030 due to the increasing levels of heat stress (https ://unfcc c.int/ news/inter natio nal-labou r-organ izati on-warns -of-heat-relat ed-job-losse s). As India is more vulnerable to the heat extremes such as heat waves, heat stress due to its projected increasing temperatures 39,[45][46][47] , knowing the impact of heat stress on the work performance helps the policy makers to take actions in order to protect the workers from the dangerous exposure of heat stress and to maintain the sustained work performance. Hence, it  www.nature.com/scientificreports/ is important to assess the future changes in the heat stress index over India where in the outdoor activities are significant. The current study is built upon the assessment of heat stress along with its associated decline in work performance over India in the present and future climate under different emission scenarios. The results of the present work may be useful to the policy makers to design a better policy tools for taking necessary steps for the welfare of the working people due to the increased heat stress in future under changing climate conditions. www.nature.com/scientificreports/

Materials and methods
Here we used Steadman's Heat Index (HI) 48 approach which combines ambient air temperature and relative humidity to estimate the heat stress. Steadman's field observations and measurements were translated into a table that showed heat stress values and the corresponding risk levels (Table 1). On the basis of the table generated by Steadman, Rothfusz (1990) 49 developed the following equation for heat index.  www.nature.com/scientificreports/ A relationship is developed based on the experimental data sets of working performance 51 , thermal comfort 52 and the decrement of work performance (P in %) has been found with the rise in temperature 53 . It is reported that the temperature and heat index values do not vary much in many areas 54 and also the sensitivity of heat-health effect estimates are not showed larger differences to measure the exposure with heat index versus temperature 55,56 . Though there is a variation, particularly in summer time, a substantial correlation between these two infer the same quantitative results 54 . Also, when using the mean temperature, the indices such as heat stress index, humidex and relative strain index, the magnitudes of all do not vary much. Average temperature performed similarly to the composite indices, but minimum and maximum temperatures performed relatively poorer. Thus, average temperature may be suitable for the development of weather-health warming systems 55 . Hence in the present study, we have used the heat stress index instead of temperature since the heat stress is a better indicator than temperature in heat-health research.
Thus, the equation for estimating the decline in work performance can be written as: The mean of three reanalysis datasets viz, the ERA-Interim reanalysis from ECMWF 57 , NCEP-DOE-Rea-nalysis2 58 and NOAA-CIRES 20th Century reanalysis 59 for daily mean surface air temperature and relative humidity have been used for the estimation of HI. Daily simulations of 18 general circulation models (GCMs) data from suite of the Coupled Model Intercomparison Project 5 (CMIP5) 60 (Table 2)   www.nature.com/scientificreports/ to an uniform grid resolution of 1° × 1° using the bilinear interpolation method 61-63 for making them consistent to compute the multi model mean and which will be helpful in studying the explicit diagnosis of heat stress on space and time scale. Matthews et al. 14 and Lee and Brenner 64 studied the global heat index, estimated from the temperature and relative humidity while Opitz-Stapleton et al. 65 and Sylla et al. 19 studied the same regionally. Lee and Berner 64 reported that these heat stress indices are used by International Organization for Standardization (ISO) and National Institute for Occupational Safety and Health (NIOSH) to measure the heat loads in order to prevent the heat illness. These formulae need to be weighted by population density in order to study the impact on human health and productivity when they are subjected to study for different countries over the globe. During summer, the spatial variations of temperature are very less over India except in the regions of Himalayas and Northeast region of India, the heat stress levels and the decline in work performance have been taken into account over entire India except these two regions. Hence, the results obtained over Himalayan and Northeastern parts need not be considered in the present study.

Results and discussions
Evaluation of summertime mean climatology of heat stress index for the baseline period (1986 to 2005) from CMIP5 Multi Model Mean (MMM) and observed reanalysis data sets along with their biases (model minus observed climatology) revealed the better performance of CMIP5 MMM in capturing the heat stress characteristics over India. Figure 1 indicates that the higher percentage of heat stress days fall in caution and extreme caution categories in the baseline period. During the four categories of heat stress days (caution, extreme caution, danger and extreme danger), the bias is less than 5% which provides a fair confidence on the reliability of model data as these biases are inevitable due to the different physics schemes used by the different GCMs that are being used as multi model mean 61 . Uncertainties in heat stress estimations that include temperature and humidity together  (Fig. 4), while the number of days with heat stress are predominantly increasing over west coast (> 32 days) under extreme caution category in both the scenarios (Fig. 5). Compared to east coast, west coast is not much vulnerable but both emission scenarios show an increase in extreme caution days and project an increase in danger days under RCP 8.5 by end of the century. It is reported that the anomalous westerlies over the Indian land mass which reduce the land sea thermal contrast may cause the hot conditions over the eastern coastal regions of India 68 Significant rise in maximum temperature over the west coast 69  www.nature.com/scientificreports/ curve for the period 2080-2099 is flatter than the other two curves, under both the scenarios, they show the shift in location towards right. Also the probabilities of low heat stress index do not show much difference for all the three time epochs and both the scenarios (left tails of the distribution). However the probability for higher values of heat stress is substantially more towards the end of the century (right tails). The PDFs look almost symmetric but they are positively skewed which clearly reveals that the heat stress index would be more severe (in danger and extreme danger category) in terms of duration and intensity towards the end of the century. Flatness of the curve implies that the kurtosis of the curves are less towards the end of the century i.e., variance or spreads of the distributions are more. Similar features can be observed over west coast region as well. The important point brought out by this analysis is substantial shift of location to the right and increase in spread of the distribution towards the higher values towards the end of the century under both the scenarios. It is clearly seen from the figures (Figs. 6 and 7) that probability of occurrence of danger and extreme danger days over the two coastal regions increase towards the end of the century because of the increase in the duration and intensity of warm days in future. This lengthening of heat stress period and increasing chance of occurrence may result in the occupational health hazards and productivity loses. A similar increase in the heat stress is being reported over different geographical regions of the globe and it has been documented that the trends in heat stress have been more detectable than those temperature 71 . The increased scenarios of heat stress are mainly due to the increasing temperatures along with the changes in the humidity. The CMIP5 models show warming over India exceeding 5 °C during summer time and this increase in temperature is not consistent within the three time slices under the two climate scenarios (under RCP 4.5 and RCP8.5) (Supp Fig. 1). A slight decrement from the mean has been observed in relative humidity which accounts for more discomfort with the increase in temperature during the twenty-first century (Supp Fig. 2). The mean surface winds at 1000 hPa and sea level pressure (SLP) during summertime for future scenarios have been analyzed to understand the surface circulation during the elevated levels of heat stress. Future surface circulation pattern under RCP 8.5 scenario shows the strengthening of easterlies over Bay of Bengal and southwesterlies over Arabian Sea reaching Indian landmass due to the increased surface pressure gradients which may have possible connection with the increased temperatures in combination with the humid conditions that probes the higher heat stress over the coastal regions (Supp. Fig. 3). As the Bay of Bengal and the Arabian Sea are continued to warm in the recent decades, the higher SSTs cause the evaporative cooling mechanisms of the oceans which transports the winds with the subdued moisture as they reach the Indian land mass 72 . This subdued www.nature.com/scientificreports/ moist winds in combination with the already raised temperatures over land, together create the higher values of heat stress. The increase in the heat stress is also attributed to the increased anthropogenic forcing, increasing trends of temperature and decreasing trends in humidity during summer time 70 . Also, the air which is heated due to the absorbing aerosols in central and northeastern parts of India sinks in northwestern and southern Indian regions that cause more heating in the recent decades 73 . The CMIP5 MMM projections for the summertime furthermore suggest anomalous south easterlies from Bay of Bengal associated with high pressure may cause high temperatures over the east and west coasts of India. This might favor the enhancement in the severity of heat stress. In other words, the background atmospheric conditions are favorable for the enhanced severity of heat stress particularly by the end of the century. Prominent features are seen in the upper troposphere like anti-cyclonic circulation over Pakistan region which is intensified towards the end of the century. The future minus present climatology suggests the positive height anomalies and anomalous anti-cyclonic flow at 500 hPa level over northern portions of India (Supp Fig. 4). Anticyclonic conditions over the study area also provide dry moist weather that increases the solar radiation flux due to subsidence of air masses. There is an indication that geopotential height (GPH) gradually strengthens during the future climate scenario. Future climate projections of amplified GPH anomalies at 500-hPa GPH indicate the significant increase of heat waves over India. Positive GPH anomalies at 500 hPa dynamically cause subsidence, clear skies, light winds, warm-air advection, and prolonged hot weather conditions at the surface 4 . Here we also observe a strong positive GPH in the upper troposphere (at 200-hpa) which might be also a causal factor for the surface heat wave and other heat extremes (supp. Figure 5). The possible link between circulation and extreme heat occurrence infers the role of horizontal temperature advection of warm dry air from the interior continental parts to coastal regions as off shore wind direction.
Climate change leaves the workers in lurch who work under extreme climates 74 . The heat stress is projected to intensify health risks of population and decrease the labour capacity, particularly over tropical and sub-tropical areas 75,76 . The decline in work performance over India is conspicuous during 2046-2065 and 2080-2099 than www.nature.com/scientificreports/ during 2016-2035 (Fig. 8). The cognitive fatigue, difficulties in mental concentration which reduces the work efficiency due to rise in temperature may result in the productivity losses 77 . In addition, the sweat evaporation stops at higher relative humidity values which hinder the cooling effect of the body leading to heat illness 78 . The climate models have shown the heat stress vulnerability in the regions of Southeast Asia, southeastern US and Northern Australia 6 . Dunne et al. 75 proposed a fit for estimating labor capacity and reported a reduction in labor work capacity upto 40% under RCP 8.5 scenario over tropical and mid latitudinal region and also found that India will undergo heat stress more quantitatively than Eurasia and greater Caribbean region. The simulated levels of moist heat stress obtained from RCP 8.5 show the substantial rise in wet bulb temperature extremes over India, China, Southeast Asia, interior South America, and western Africa 79 . The urban heat island effect caused by the heat and moisture from the concrete structures that affects the incoming and outgoing solar radiation also has an impact on the heat stress in urban areas 80 . The impacts of climate change on working performance in India further worsen in view of RCP 8.5 scenario where India is likely to experience 3 to 4 degrees of warming 81 . The decrement in the work performance is highly susceptible in the southeastern coastal regions of India under the RCP 4.5 scenario while the same extends to total country except in the belt of Himalayas under the RCP 8.5 scenario. Our assessment showed a quantitative decline in work performance to up to 40% under the RCP8.5 and 35% under the RCP 4.5 scenario. The initial decline in work performance is conspicuous over the east coastal regions of India and later expanded to most of the regions in the country by end of the century. As India witnesses the notorious climate change impacts over the coastal regions, where 250 million people live within 50 km of the coast line 82 , the decline in work performance will be a major threat to the livelihood of the coastal communities' . However, a rigorous evaluation of work performance needs to be done based on individual and population-based characteristics. We provide here a common understanding in the projected decreasing levels of work performance as an overall assessment though the parameters such as local work place conditions 83 , construction of buildings 84 , clothing 85 and physical work intensity 86 have not been considered as input parameters. As suggested by Brode et al. 44 , Dunne et al. 75 and Buzan and Huber 79 , estimation of the work performance needs to be contextualized and should elevate the role of changing the duty hours of labor, type of the labor. Our effort in the present study is to report the general features of declining work performance under future warm and humid conditions under global warming scenario. Also, the present study considers the mean temperature for estimation of heat stress that will provide a fair understanding of work performance. However, the diurnal cycle of the temperature may also play a key role under changing work hours of the labor. Finally, the practical implications of these elevated heat stress 'need the support of engineering, social and policy-making decisions' 79 .

conclusions
Heat stress is expected to intensify, due to unprecedented warming during the twenty-first century. The increasing levels of heat stress pose a great challenge to mitigate and to adapt proper mechanisms to confront it. Decline in work performance associated with the detectable trends of heat stress over India point out the dire need to frame the sage work procedures, preventing the heat-based illness to reduce the health hazards as well as to maintain the work performance. The coastal regions of India (east and west coast) are found to be more vulnerable to heat stress impacts by showing a perceptible increase in the notorious impact days and a decline of 30 to 40% in the work performance, particularly in east coast region.