Abstract

In accordance with the China Meteorological Administration definition, this study considered a weather process with a maximum surface temperature of ≥35°C for more than three consecutive days as a heatwave event. Based on a dataset of daily maximum temperatures from meteorological stations on the North China Plain, including ordinary and national basic/reference surface stations, the intensity-area-duration method was used to analyze the spatiotemporal distribution characteristics of heatwave events on the North China Plain (1961–2017). Moreover, based on demographic data from the Statistical Yearbook and Greenhouse Gas Initiative (GGI) Population Scenario Database of the Austrian Institute for International Applied Systems Analysis (IIASA), population exposure to heatwave events was also studied. The results showed that the frequency, intensity, and area of impact of heatwave events on the North China Plain initially decreased (becoming weaker and less extensive) and then increased (becoming stronger and more extensive). Similarly, the trend of population exposure to heatwave events initially decreased and then increased, and the central position of exposure initially moved southward and then returned northward. Population exposure in the eastern Taihang Mountains was found significantly higher than in the western Taihang Mountains. In relation to the change of population exposure to heatwave events on the North China Plain, the influence of climatic factors was found dominant with an absolute contribution rate of >75%. Except for 2011–2017, increase in population also increased the exposure to heatwaves, particularly in the first half of the study period. Interaction between climatic and population factors generally had less impact on population exposure than either climatic factors or population factors alone. This study demonstrated a method for assessing the impact of heatwave events on population exposure, which could form a scientific basis for the development of government policy regarding adaption to climate change.

1. Introduction

Climate change includes not only changes in mean climate but also variation in weather extremes [1]. Over the past century, because of the combined effects of human activities and natural factors, Earth’s climate has undergone significant change in terms of the primary features of warming and this change has shown an increasing trend in recent years [2, 3]. The trend has led to a wide range of climate anomalies, especially an increased frequency of occurrence of regional extreme high-temperature events [4]. The Fifth Assessment Report issued by the United Nations Intergovernmental Panel on Climate Change in 2014 showed that the global surface temperature increased by 0.85°C from 1880 to 2012 [5]. In different regions of the world, the annual and daily maximum temperatures have both shown a significant upward trend since the latter part of the 20th century [6, 7]. For most of China, the trend of increase in air temperature has accelerated since the 1980s [8]. From a regional perspective, the characteristics of high temperatures are complex. Overall, significant positive trends in the frequency of occurrence of extreme maximum temperature events have been observed in most of China, and the North China Plain is one area that has experienced the most significant warming.

Increased temperature will lead to an increase in the risk of heat-related diseases and death in humans [9]. The risk of morbidity and accidental death among urban vulnerable populations and among urban and rural outdoor workers during a period of heat stress has been identified as one of the key risks of climate change in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [5]. Changes in heatwaves are an important aspect of climate change, and their impact on both natural ecosystems and human society might be more severe than the effects of the rise in average temperature [1013]. The latest Hadley Centre Global Historical Climatology Network Daily observations show that the intensity, frequency, and duration of global heatwave events are increasing [14]. In recent years, heatwaves have surpassed other natural disasters such as hurricanes, lightning, tornadoes, floods, and earthquakes to become the primary cause of fatalities in America related to natural disasters [1517]. Moreover, heatwaves have killed more people in Australia than any other natural hazard since 1900 [18]. It is estimated that around 40,000 people in central and western parts of continental Europe died during a record heatwave that occurred in the summer of 2003 [19, 20]. In 2013, heatwaves occurred frequently in most parts of China. For example, Changsha, Chongqing, and Hangzhou recorded 57, 49, and 47 hot days (daily maximum temperature ≥35°C), respectively, during which the highest recorded temperatures reached 38.2, 40.9, and 40.5°C, respectively. These heatwave incidents led to a direct economic loss of 83.56 billion RMB [21]. Research on high-temperature disasters has attracted widespread attention from government departments and the scientific community, and it has become the focus of climate research under the background context of global warming.

In recent years, studies on population exposure to natural disasters have achieved many results [2227]. Research on extreme climate events has developed from value analysis at a single site [2830] to related investigations of intensity, area of impact, and duration [26, 31]. Based on the intensity-area-duration (IAD) method, Zhai et al. [31] developed an approach for identification of a regional extreme event and analysis of its spatial coverage and duration. Extreme heatwave events can be more serious than a single high-temperature event [10]. However, few studies have considered regional persistent extreme heatwave events and the variation of associated population exposure.

The population of the North China Plain accounts for approximately 30% of the total population of China. Furthermore, the North China Plain represents the center of China’s political, economic, cultural, transportation, education, and tourism activities. Therefore, it is of great scientific and practical importance to study the evolution of heatwave events and associated population exposure on the North China Plain. This study was based on daily maximum temperature data (1961–2017) from national meteorological stations on the North China Plain, in combination with demographic data and Greenhouse Gas Initiative (GGI) Population Scenario Database of the Austrian Institute for International Applied Systems Analysis (IIASA) [32, 33]. First, the frequency, intensity, area of impact, and spatial variation characteristics of heatwave events in North China are analyzed by grid unit scale. Second, the temporal and spatial evolution characteristics of population exposure to heatwave events are studied. Finally, attribution analysis is performed on the factors affecting population exposure. The main objectives of this study were to reveal the contributions of various factors to population exposure to past heatwave events on the North China Plain and to provide a scientific basis for the development of strategies effective in helping humans cope with the adverse effects of such heatwave events.

2. Data and Methods

2.1. Study Area

The area of the North China Plain considered in this study comprised the region bounded by 33°–42°N, 110°–120°E, which includes Beijing, Tianjin, Inner Mongolia, Hebei, Henan, Shanxi, Shandong, Jiangsu, Anhui, and Shaanxi (Figure 1). This region covers an area of approximately 9.0 × 105 km2 and it contains approximately 30% of the total population of China. The North China Plain has a warm temperate semiarid monsoon climate, which has the characteristics of cold and snow in winter, a dry and windy spring, hot rainy summer, and cool dry autumn.

2.2. Data Description

The climate dataset used in this study was derived from daily maximum temperature data (1961–2017) recorded at national meteorological stations on the North China Plain, which were integrated by the National Meteorological Information Center of the Chinese Meteorological Administration. Following quality control, the daily maximum temperature data were interpolated to a 0.5° × 0.5° grid using the inverse distance weighting method.

The North China Plain population data were obtained from the Chinese Demographic Yearbook. Using the GGI Population Scenario Database of the IIASA (which has spatial resolution of 0.5° × 0.5° and temporal resolution of 10a) as the background of population distribution, the population data of the North China Plain provinces from the Chinese Statistical Yearbook were interpolated to a 0.5° × 0.5° grid using the grid weight calculation method.

2.3. Method
2.3.1. Heatwave Event Identification

In accordance with the definition of heatwaves by the China Meteorological Administration, this study considered a weather process with a maximum temperature ≥35°C for more than three consecutive days as a heatwave event. This definition focuses on the hazard of persistently high air temperature. When a grid met this definition, it was assumed that there was a heatwave at that grid. Thus, we analyzed the frequency of heatwave events (i.e., number of heatwave events), intensity of heatwave events (i.e., maximum daily maximum temperature of a heatwave event), and areas affected by heatwave events (i.e., maximum area covered by the heatwave).

2.3.2. Intensity-Area-Duration (IAD) Method

This study adopted the IAD method, which links the three most important features of extreme events [26, 31, 34, 35]. This method has been applied successfully in analyses of drought and extreme precipitation events in China [26, 31, 34]. The principle of the IAD method is to identify grids that neighbor the daily maximum temperature ≥35°C over a given timescale and to select continuous areas as an extreme heatwave event. The area of impact is determined by the summation of all grids, while the intensity is determined by the maximum temperature of all grids included in one event. The IAD method can be used to study the simultaneous changes in intensity and area of impact over a given duration. The steps required in the process are as follows.

First, for a given timescale (here, ≥3 d), calculate the intensities for the different periods for all grids. The grid point with the highest intensity is regarded as the “center with highest intensity” of a regional extreme daily maximum temperature event. Second, among the surrounding eight grids, find the one with the second highest intensity to establish the “center with second highest intensity.” Third, among the grids surrounding the “center with second highest intensity,” find the one with the third highest intensity. Determine all those grids with daily maximum temperature greater than the threshold and combine them all as a regional extreme daily maximum temperature event. Fourth, determine another “center with highest intensity” and repeat the above steps, deducing all the regional extreme daily maximum temperature events over a given timescale.

2.3.3. Spatialization of Population

Population is an imperative factor in exposure assessment in relation to heatwave disasters. The spatial distribution of the population was established using the IIASA GGI Population Scenario Database (resolution: 0.5° × 0.5°, time: 1990, 2000, and 2010) based on the total number of provinces in the Chinese Statistical Yearbook. In this process, the weight coefficient is interpolated to each grid point and calculated as follows:where denotes the population at grid i in province , is the population of the GGI Population Scenario Database at grid i in province j, is the total population of the GGI Population Scenario Database in province , is the total population of province j, m is the total number of provinces, and n is the total number of grids in province j.

In this study, the population distribution of the North China Plain during 1961–1990, 1991–2000, and 2001–2017 was based on the GGI Population Scenario Database in 1990, 2000, and 2010, respectively.

2.3.4. Definition of Population Exposure and Analysis of Contribution Rate

The number of people exposed to a heatwave event was defined as the population exposure [36]. Population exposure to heatwave events is affected by the frequency of heatwaves and the number of people. Consequently, population exposure is affected by the combined effects of climatic factors (change in frequency of heatwave events and constant population size), population factors (constant frequency of heatwave events and change in population size), and climatic-population factors (changes in both frequency of heatwaves events and population size). Changes in population exposure to heatwaves can be expressed as [37]

For ease of comprehension, it is assumed that the contribution rate of increased population exposure is positive and that the contribution rate of different influencing factors in the change of population exposure to heatwaves is as follows:

Rate of contribution of climatic factors:

Rate of contribution of population factors:

Rate of contribution of combined climatic-population factors:where F and P are the frequency of heatwave events and population size, respectively, F× ∆P is the influence of population factors, P× ∆F is the influence of climate factors, and ∆F× ∆P is the combined influence of climatic-population factors.

3. Results

3.1. Frequency, Intensity, Area of Impact, and Spatial Changes of Heatwave Events

The frequency of heatwaves reflects the number of local heatwave events, where a higher frequency usually means greater impact. Figure 2 shows the frequency of heatwave events on the North China Plain during 1961–2017. It can be seen that the decadal mean frequency of heatwave events on the North China Plain shows a significant initial trend of decrease followed by an increase. During the period with a significant downward trend (1961–1990), the decadal mean frequency of heatwave events during 1961–1970, 1971–1980, and 1981–1990 decreased from 30.0 to 17.6 and 10.7 events, respectively. During the period with an upward trend (1990–2017), the decadal mean frequency of heatwave events during 1991–2000, 2001–2010, and 2011–2017 increased to 19.4, 23.1, and 38.1 events, respectively. The highest decadal mean frequencies of 60, 59, and 58 events occurred in 2017, 2013, and 2014, respectively. Based on the daily maximum temperature data, we found that changes in air temperature caused these changes in the frequency of occurrence of heatwave events.

Heatwave intensity reflects the strength of heatwave events. Figure 3 shows the intensity of heatwave events on the North China Plain. It can be seen that the most and least intense heatwave events on the North China Plain occurred during 1961–1970 and 1971–1980, respectively. Since the 1970s, there has been a trend of increase in heatwave intensity. The average intensity of heatwaves during 1971–1980 was approximately 36.8°C; during 1981–2000, the average intensity of heatwave events remained around 36.9°C, and the mean intensity during 2001–2010 and 2011–2017 reached 37.1 and 37.3°C, respectively. During 1961–2017, the average intensity of heatwaves on the North China Plain always remained >36°C, and there were four years with intensity >38°C.

The area of impact of heatwave events on the North China Plain during 1961–2017 is shown in Figure 4. The chronological trend of the mean area of impact of heatwave events shows a decrease during 1961–1990, i.e., it decreased from 0.65 × 106 km2 in 1961–1970 to 0.19 × 106 km2 in 1981–1990. During 1961–1990, the area of impact of heatwave events on the North China Plain decreased by an annual average of 2.2 × 104 km2. During 1981–2017, the area of impact of heatwave events increased significantly. For example, the average annual area of impact in 2011–2017 was 4.5 times that during 1981–1990. During 1991–2017, the area of impact of heatwave events on the North China Plain increased by an annual average of approximately 2.9 × 104 km2. The maximum area of impact of >1.6 × 106 km2 appeared in 2014.

The spatial distribution of the annual average frequency of heatwaves on the North China Plain shows that the frequency of heatwaves decreased during 1961–1990 and that the area of impact gradually diminished (Figures 5(a)5(c)). The center of high frequency gradually narrowed from southeast Hebei, central and eastern Henan, western Shandong and northwest Anhui to small areas of southern Hebei, central Henan, and northern Anhui. Conversely, the frequency and area of impact of heatwave events on the North China Plain increased significantly during 1981–2017 (Figures 5(c)5(f)). The center of high frequency gradually expanded from southeastern Hebei and central Henan to southern Hebei, southern Shandong, and northern-central Henan. According to the spatial distribution of the annual average frequency during 1961–2017 (Figure 5(g)), heatwave events have occurred in most parts on the North China Plain, and the center of high frequency (frequency: >20 times) has been located mainly in southeastern Hebei, eastern Henan, and northern Anhui. It is evident that topography determines the spatial distribution of the frequency of heatwaves, i.e., to the east of the Taihang Mountains, low terrain means high frequency, and to the west of the Taihang Mountains, high terrain means low frequency.

The spatial distribution of the frequency of heatwaves is related to population exposure. For a population of a certain size, more frequent heatwave events will mean greater population exposure. This means that if the current size of the population is maintained, an increase in the frequency of future heatwave events will lead to an increase in population exposure.

3.2. Population Exposure to Heatwave Events

The severity of the impact of a heatwave event depends not only on the event itself but also on the levels of exposure and vulnerability of the population. The level of exposure is the extent to which the adverse effects of heatwave events affect the population, economy, and other factors [39, 40]. Population exposure is defined as the size of the population within an area affected by a heatwave. The variation in population exposure to heatwave events on the North China Plain shows a downward trend during 1961–1990 and a significant upward trend during 1981–2017 (Figure 6). During 1981–1990, population exposure decreased by approximately 58.9% (in comparison with 1961–1970) at an annual average rate of 6.8 million. During 2011–2017, population exposure increased by approximately 331.6% (in comparison with 1981–1990) at an average annual rate of around 12 million. According to current trends, population exposure on the North China Plain will exceed 1 billion in the next decade.

The spatial distribution of population exposure to heatwaves on the North China Plain is shown in Figure 7. Comparison with Figure 5 reveals that the spatial distribution of population exposure is related to the spatiotemporal distribution of heatwave events. During 1961–2017, population exposure to heatwaves occurred in Shandong, northern Jiangsu, northeast and southern Hebei, Beijing, Tianjin, Henan, southern Shanxi, and other areas. The most densely populated area affected by heatwave events was southeastern Henan (Figure 7). Population exposure to heatwaves on the North China Plain decreased during 1961–1990 (Figures 7(a)7(c)) and the extent of the area of high values gradually retreated from north to south. The population affected by heatwave events showed an increasing trend during 1981–2017 (Figures 7(c)7(f)) and the extent of the area of high values gradually expanded from south to north. The largest period of population exposure occurred in 2011–2017 (Figure 7(f)). Grid values with population exposure of >40 million occurred in >20% of regions. Figure 7(g) shows the spatial distribution of population exposure on the North China Plain during 1961–2017. It can be seen that population exposure to heatwaves decreased from south to north and that the center of the highest values was located in southeastern Henan. Population exposure is greatly affected by elevation, i.e., to the east of the Taihang Mountains, population exposure is large, and to the west of the Taihang Mountain, population exposure is small. Although Beijing and Tianjin have high population density, population exposure is low in these regions because of the occurrence of fewer heatwave events.

3.3. Analysis of Factors Influencing Population Exposure to Heatwaves

As can be seen from Figure 8, the population of the North China Plain is gradually increasing. The population east of the Taihang Mountains is significantly larger than to the west. The high-density population areas are Beijing, Tianjin, and southeastern parts of the North China Plain.

Changes in population exposure to heatwave events depend on changes in the spatial distributions of climate and population. Table 1 shows the rate of contribution of each factor that influences population exposure, where a negative contribution indicates that the contribution is opposite to the trend of population exposure and where the absolute value indicates the rate of that contribution. The change of population exposure to heatwave events on the North China Plain has been dominated by climatic factors, although the rate of contribution decreased gradually and changed from a negative contribution to a positive contribution during 1961–2010. During 2011–2017, in comparison with 2001–2010, the population influenced by climatic factors increased by up to 0.02 billion (180%) because of the sharp increase in the frequency of occurrence of heatwaves. The most influential period of population factors was during 1971–1980 when the population influenced increased by 0.54 billion (60.7%) in comparison with 1961–1970. Subsequently, the population influenced by population factors gradually decreased as the contribution of population factors changed from positive to negative. In general, interaction between climatic and population factors had less impact on population exposure than either climatic factors or population factors alone.

Climatic factors have had the greatest impact on population exposure, whereas the effect of population factors has had a decreasing trend, which has resulted ultimately in a negative rate of contribution. This is because of the significant increase in heatwaves, while some provinces have experienced slow or even no population growth in recent years.

4. Conclusions and Discussion

Since the 1990s, global warming has had a significant effect on heatwaves in China and this is expected to continue. Using the IAD method based on gridded data, which combines the three most important elements (intensity, area of impact, and duration) of heatwave events, this study investigated the spatiotemporal characteristics of regional heatwaves on the North China Plain. This approach, which differs from previous analyses based solely on individual stations, allows greater understanding of heatwave events. In addition, gridded temperature data can be more reasonably combined with population data, which is helpful in assessment of the impact of heatwave events on a population.

Based on daily maximum temperature data and spatiotemporal data of population distribution on the North China Plain during 1961–2017, we calculated and analyzed the spatiotemporal evolution characteristics, population exposure change, and factors influencing heatwave events on the North China Plain over the past 57 years. The main conclusions are as follows:(1)Heatwave events on the North China Plain exhibited obvious spatiotemporal change. From the perspective of the decadal mean, the frequency and area of impact of heatwave events both declined initially and then increased, with the lowest values during 1981–1990. The intensity of heatwave events showed an initial downward trend, remained reasonably stable during 1981–2000, and then increased during 2001–2017. From the perspective of spatial change, the area of impact and frequency of heatwave events both decreased at first and then increased, and the center of heatwave events receded southward and then extended northward. From the spatial distribution of heatwave events during 1961–2017, it was found that areas with high incidence of heatwave events were located mainly in southeastern Henan and southern Hebei. The frequency of occurrence of heatwaves was found to be affected by elevation, i.e., heatwaves occurred more frequently to the east of the Taihang Mountains in comparison with areas to the west. Thus, the findings of this study in relation to the characteristics of heatwaves on the North China Plain are similar to other earlier studies [14, 41, 42].(2)Under the background of global warming and the given spatial distribution of population, it was found that the trend of population exposure on the North China Plain showed significant decrease during 1961–1990 followed by an increase during 1991–2017. Spatially, during 1961–1990, the area of population exposure gradually narrowed and the center of high values moved southward. Conversely, during 1991–2017, the center of high values moved northward and the number of people affected by heatwaves increased. During 1961–2017, the region with the highest population exposure was located in southeastern parts of Henan.(3)Population exposure on the North China Plain was found to be influenced by climatic factors, population factors, and combined climatic-population factors. In the change of population exposure to heatwave events, climatic factors played the leading role, population factors increased the population exposure to heatwaves in all decades except for a marginal decrease from 2001 to 2010 to 2011–2017, but there was a decreasing trend in the magnitude of this contribution. During 1961–2010, there was a decreasing trend in the interaction between climatic factors and population factors. The population exposure to climatic factors first decreased and then increased, consistent with the change of frequency of occurrence of heatwave events. The increase in the population in heatwave event areas has led to an increase in the number of people exposed to heatwave events. According to the climatic and demographic trends, it was established that climatic factors will continue to contribute positively and will remain the most important driver in the future, whereas population factors will contribute at a low or even negative rate, consistent with the conclusions of Huang et al. [37].

The effect of a heatwave event on population exposure is not only related to its strength, area of impact, and duration but also to the spatiotemporal distribution of the population. Limited by the temporal resolution of the IIASA GGI Population Scenario Database, the population distribution data used in this study were the background fields of the spatial distribution of population in 1990, 2000, and 2010. Thus, we did not consider population distributions with higher temporal resolution. This will be a focus of future research.

The mechanism maintaining high temperatures on the North China Plain has not yet been fully elucidated. Xie et al. [43] suggested that the continental subtropical high is the main synoptic system responsible for high temperatures on the North China Plain. Conversely, Wei et al. [43] proposed that the occurrence of summer temperature anomalies on the North China Plain is determined primarily by east-west anomalies of the western Pacific subtropical high. However, the question of how the western Pacific subtropical high and continental high pressure might interact to produce high temperatures on the North China Plain will require further in-depth research.

In addition to population exposure, the disaster risk associated with heatwaves is related to factors such as high temperature adaptability of different populations and to the regional level of defense against the effects of heatwaves. Age, gender, occupation, disease status, medical resources, and economic level all affect the response of a population to the effects of heatwaves. Therefore, it will be necessary to undertake further study on the spatiotemporal evolution of heatwave disaster vulnerability under future scenarios, based on the findings of this study, to improve the level of risk management in relation to heatwave events.

Our results provide valuable evidence of the suitability of the proposed method for assessing the impact of heatwave events on population exposure, and they could provide a scientific basis for the development of government policy regarding adaption to climate change.

Data Availability

The daily maximum temperature data (1961–2017) used to support the findings of this study were supplied by http://data.cma.cn/under account and so cannot be made freely available. The North China Plain population data used to support the findings of this study have been deposited in the Chinese Demographic Yearbook (http://www.stats.gov.cn/tjsj/ndsj/). GGI Population Scenario Database of the IIASA used to support the findings of this study were supplied by http://www.iiasa.ac.at/web-apps/ggi/GgiDb/dsd?Action=spatial under account and can be freely available. All other data are available from the authors upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the CDM Fund “Jiangxi Climate Change Adaptation Program” (2014102) and bilateral cooperation project between the Natural Science Foundation of China and the Pakistan Science Foundation (41661144027). The CMA Climate Change Science Fund (CCSF (201722, 201810)) provides a policy-oriented training course for PhD students. The authors are thankful for the support of the High-level Talent Recruitment Program of the Nanjing University of Information Science and Technology (NUIST).