Next Article in Journal
A Comprehensive Comparison of Machine Learning and Feature Selection Methods for Maize Biomass Estimation Using Sentinel-1 SAR, Sentinel-2 Vegetation Indices, and Biophysical Variables
Next Article in Special Issue
Complementarity Characteristics of Actual and Potential Evapotranspiration and Spatiotemporal Changes in Evapotranspiration Drought Index over Ningxia in the Upper Reaches of the Yellow River in China
Previous Article in Journal
An Algorithm for Retrieving the 2-D Distribution of Moderate Rain by X-SAR
Previous Article in Special Issue
PS-InSAR Based Monitoring of Land Subsidence by Groundwater Extraction for Lahore Metropolitan City, Pakistan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inter-Comparison of Diverse Heatwave Definitions in the Analysis of Spatiotemporally Contiguous Heatwave Events over China

1
Department of Atmospheric Science, School of Environmental Studies, China University of Geosciences, Wuhan 430074, China
2
Centre for Severe Weather and Climate and Hydro-Geological Hazards, Wuhan 430074, China
3
Laboratory of Critical Zone Evolution, School of Geography and Information Engineering, China University of Geosciences, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(16), 4082; https://doi.org/10.3390/rs14164082
Submission received: 18 July 2022 / Revised: 12 August 2022 / Accepted: 18 August 2022 / Published: 20 August 2022

Abstract

:
A heatwave (HW) is a spatiotemporally contiguous event that is spatially widespread and lasts many days. HWs impose severe impacts on many aspects of society and terrestrial ecosystems. Here, we systematically investigate the influence of the selected threshold method (the absolute threshold method ( A B S ), quantile-based method ( Q T L ), and moving quantile-based method ( Q T L m o v )) and selected variables (heat index ( H I ), air temperature ( T a i r )) on the change patterns of spatiotemporally contiguous heatwave (STHW) characteristics over China from 1961–2017. Moreover, we discuss the different STHW change patterns among different HW severities (mild, moderate, and severe) and types (daytime and nighttime). The results show that (1) all threshold methods show a consistent phenomenon in most regions of China: STHWs have become longer-lasting (6.42%, 66.25%, and 148.58% HW days ( HWD ) increases were found from 1991–2017 compared to 1961–1990 corresponding to A B S , Q T L , and Q T L m o v , respectively, as below), more severe (14.83%, 89.17%, and 158.92% increases in HW severity ( HWS ) increases), and more spatially widespread (14.92%, 134%, and 245.83% increases in the summed HW area ( HWA sum )). However, the HW frequency ( HWF ) of moderate STHWs in some regions decreased as mild and moderate STHWs became extreme; (2) for threshold methods that do not consider seasonal variations (i.e., A B S and Q T L ), the spatial H I exceedance continuity was relatively weak, thus resulting in underestimated STHW characteristics increase rates; (3) for different variables defining STHWs, the relative changing ratio of the H I -based STHW was approximately 20% higher than that of the T a i r -based STHW for all STHW characteristics, under the Q T L m o v threshold; (4) for different STHW types, the nighttime STHW was approximately 60% faster than the daytime STHW increase considering the Q T L threshold and approximately 120% faster for the Q T L m o v method. This study provides a systematic investigation of different STHW definition methods and will benefit future STHW research.

Graphical Abstract

1. Introduction

Heatwaves (HWs) severely impact many aspects of society and terrestrial ecosystems [1,2,3,4]. For example, HWs can directly cause organ damage and a range of cardiovascular diseases (e.g., myocardial infarction and heat exhaustion) and indirectly increase the risk of spreading viruses [5,6,7,8]. Additionally, HWs can reduce food and water supplies and increase potential conflicts, such as the mortality caused by HW and other environmental and social problems. For example, the 2012 extreme HW in the United States led to a maize production decrease of 13% compared to the 2011 level [9]; the 2013 extreme HW in southern China led to a net reduction of 101.54  Tg C in regional carbon sequestration in two months and the most extensive crop failures since 1960 [10].
In the past few decades, many researchers have devoted themselves to revealing the change patterns of HWs in the warming climate. However, different researchers have used different methods to define HW events. It is generally accepted that an HW is an episode when the temperature variable exceeds a certain threshold. However, when selecting a more particular definition, divergences arise regarding which threshold method and temperature variable to use [11,12,13]. For the threshold, there are three prevalent methods: the absolute threshold method ( A B S ), the quantile-based method ( Q T L ), and the moving quantile-based method ( Q T L m o v ) [12,13]. For example, a HW is recognized by the Royal Netherlands Meteorological Institute as a period when the daily maximum temperature is greater than 25 °C for more than five consecutive days. Following the World Meteorological Organization’s (WMO) suggestions, HWs occur when the daily maximum temperature surpasses the average maximum temperature by 5 °C for five consecutive days. In China, HWs are defined as events exceeding 35 °C and lasting at least three consecutive days, according to the Chinese National Meteorological Administration. These commonly used threshold features are summarized in Table 1. Regarding the temperature variable, some studies have used air temperature ( T a i r ) [14,15,16,17]. However, many studies have argued that the heat index ( H I ) is more relevant to human-perceived temperature [18] and is thus more suitable to define HWs [19,20,21,22]. Although some significant progress has been made, the impacts of variable selection on the spatiotemporal change pattern of HWs are not well understood.
On the other hand, HW is a kind of spatiotemporally continuous event that are spatially widespread and last many days, transferring from one region to another, and merging or splitting over time during their development and diminution. However, most previous studies have ignored these HW characteristics and studied HWs simply at the grid or site scale. With advances in the spatiotemporal clustering technique [23,24], it has become feasible to track spatiotemporally contiguous heatwaves (STHWs) and study their change patterns on the 3-D event scale [12,25,26,27]. However, HW studies performed at this 3-D event scale are much fewer than those performed at the grid and site scales. Consequently, the problem of the effects of the variable and threshold method selection on the STHW changing pattern is increasingly prominent.
Therefore, in this study, we mainly attempt to illustrate the influence of threshold method selection (the absolute threshold method ( A B S ), quantile-based method ( Q T L ), and moving quantile-based method ( Q T L m o v )), variable selection (heat index ( H I ), air temperature ( T a i r )) on the change patterns of spatiotemporally contiguous heatwaves (STHWs) characteristics over China from 1961 to 2017. In addition, we provide an in-depth discussion on the drivers of differences in the intercomparison of HW definitions.

2. Data and Methods

2.1. Data

The daily maximum air temperature ( T m a x ), minimum air temperature ( T m i n ), and relative humidity ( R H ) are used in this study. The T m a x and T m i n data are derived from CN05.1 [28] at a spatial resolution of 0.25° covering 1961–2021 period. The CN05.1 dataset is interpolated from 2416 stations in mainland China with strict quality control [28]. Previous studies have revealed that unhomogenized R H could lead to a fake decreasing trend of H I in South China [29,30]. To avoid this problem, we used ChinaRHv1.0, a homogenized dataset containing site observations from 746 stations on the Chinese mainland from 1960–2017 [30]. We interpolated ChinaRHv1.0 into 0.25° by the same method as that applied to CN05.1 [28], where the thin-plate spline method was proposed for the climatology and the angular distance weighting method was proposed for the anomaly. The interpolated climatology-plus-anomaly result forms the final gridded R H dataset. Limited by the R H coverage period, we analyze STHW only from 1961 to 2017.

2.2. Heat Index

In previous research, either T a i r or H I was used to evaluate HW. However, there is no definite conclusion as to which of these is more suitable. In this paper, both variables are used for comparison. H I (in °C), also known as the apparent temperature, is more correlated with the human-perceived temperature [31,32] and is calculated through Rothfusz regression [21,22,33]:
H I = 8.784695 + 1.61139411 T a i r 2.339549 R H 0.14611605 T a i r R H 1.2308094 × 10 2 T a i r 2 1.6424828 × 10 2 R H 2 + 2.211732 × 10 3 T a i r 2 R H + 7.2546 × 10 4 T a i r R H 2 + 2.211732 × 10 3 T a i r 2 R H
The general H I is calculated when T a i r is above 4.4 °C (40 °F); otherwise, H I is equal to T a i r , where T a i r can be T m a x or T m i n . When R H and T a i r encounter the following situations, further adjustments are applied to H I [22,33]:
H I = H I ( 13 R H ) 4 × 1 1.8 T a i r 63 17 , R H < 13 %   a n d   26.7 ° C < T a i r < 44.4 ° C H I = H I + R H 85 10 × ( 11 0.36 T a i r ) , R H > 85 %   a n d   26.7 ° C < T a i r < 30.6 ° C

2.3. HW Thresholds

If the temperature of a grid cell exceeds its corresponding threshold, we consider that a HW has occurred here. The determination of thresholds consists of two main components: quantiles and calculation methods. The detailed classification is described below:
The quantile can also refer to the HW severity levels depending on its magnitude: (1) 90th (i.e., mild level), (2) 95th (i.e., moderate level), and (3) 99th (i.e., severe level). When the quantile is greater than the 99th percentile, only a few HW events are identified. Such a small sample size would considerably impair the reliability of the estimated trends. This could lead to an uncertainty factor that is too large to reflect the variation trend patterns of most HW events.
Regarding the calculation methods, three methods are selected, including (1) A B S , (2) Q T L , and (3) Q T L m o v .
(1)
For the A B S method, all grid cells share the same threshold and eastern China (longitude ≥ 105 °E) is used to compute the threshold. First, the regional area-weighted average of T a i r or H I is calculated at each time step. The thresholds are determined in accordance with certain percentiles (i.e., 90th, 95th and 99th corresponding to mild, moderate and severe HWs, respectively) during the 1961–1990 reference period from the time series.
(2)
For the Q T L method, each grid has a unique threshold. For each grid, the threshold is calculated as a certain percentile during the 1961–1990 reference period.
(3)
For the Q T L m o v method, the threshold varies among each day of the year ( d o y ) from 1 to 366 in the Julian calendar. For each grid, a 15-day moving window (i.e., 7 days prior and posterior to the corresponding calendar day sampled from each year) for each d o y is used to create a subset, and a specific quantile calculated during the 1961–1990 reference period is used as the threshold corresponding to the center date. For instance, all 1 January to 15 January data from 1961–1990 are selected, and the threshold for 8 January is the quantile of these 450 days. Due to the dramatic temperature changes that have occurred in China since the 1990s, it is appropriate to take 1961–1990 as the reference period for a historical comparison and climate change monitoring; this application is also recommended by the WMO [34].

2.4. STHW Detection

In spatiotemporal detecting, the HW spatial-state array (i.e., whether the temperature is greater than or equal to the threshold) is used as the input, and temporally and spatially contiguous grid cells are connected to the same STHW event. This idea has been applied in drought studies previously [35]. The 3-D connectivity algorithm is relatively common but is rarely used in HW identification tasks. It is often used for lesion tracking in medical and 3D image denoising, and similar functions can be found in the Python package scikit-image, MATLAB tool P i x e l   C o n n e c t i v i t y , R package n e u r o i m , and so on [36].
Here, we rewrite the spatiotemporal clustering algorithm of Samaniego [24] to detect STHWs. Compared to the original version, ours has an excellent computing efficiency (e.g., 30~50 times faster, thanks to the speed of the Julia language), low memory consumption (e.g., uses the disjoint-set data algorithm) and can also connect contiguous temporally related HW events (e.g., all merging, splitting, expanding, and receding HW events associated with a specific event) into the same STHW event. Figure 1 briefly illustrates the main steps of this process, and the computational implementations are as follows:
(1)
The original data are converted to the HW spatial-state array. Initially, the original data have three dimensions (i.e., longitude × latitude × time). For each voxel, if the value is greater than the corresponding threshold, it is marked as “True” and considered as HW grid; otherwise, it is marked as “False”. The threshold is calculated as described in Section 2.3. Therefore, the data are converted to an HW spatial-state array that is stored in a “True/False” Boolean format.
(2)
Spatially adjacent HW grids are connected to spatial HW events. In this step, the HW state array is cyclically traversed moment-by-moment in a top-to-bottom, left-to-right order to search for spatially adjacent voxels in each time layer. If an HW grid is detected, it is given a unique c l u s t e r I d . After that, according to the 4-connected neighborhood principle (i.e., two adjoining cells are part of the same object if they are both on and are connected along the horizontal or vertical direction), all recursively searched edge-adjacent HW grids are given the same c l u s t e r I d . Here, a merged spatial HW event with a grid number less than n c o n n e c t is eliminated, because a small isolated HW event is usually unable to develop into a regional or impactful STHW event; n c o n n e c t is set to 16.
(3)
The spatial HW events are connected to spatiotemporally contiguous HW events. If there are at least n o v e r l a p grids at the same location in the preceding and following time steps from two spatial HW events, the events are merged into one STHW event and given a new c l u s t e r I d ( c l u s t e r I d n e w = y e a r × 1000 + c l u s t e r I d o l d , where c l u s t e r I d o l d is the c l u s t e r I d of the previous spatial HW event). Here, n o v e r l a p is set to correspond to n c o n n e c t . Because HW events usually have a small scale when appearing and dying, this threshold can be set to better distinguish different STHW events from the onset and demise. Finally, the c l u s t e r I d values are stored in three dimensions as original data (i.e., longitude × latitude × time). Due to the thresholds, the identified STHW can ensure that the spatial coverage is not too small and lasts at least 2 days.
Figure 1. The proposed method for detecting STHWs. (a) The identified STHW events from 28 June to 1 July 2017, in China; the features of merging and splitting events are labeled in the third step. (b) Explanation of the STHW clustering process, including three steps: (1) identifying HWs; (2) spatially connecting HWs; and (3) temporally connecting HWs.
Figure 1. The proposed method for detecting STHWs. (a) The identified STHW events from 28 June to 1 July 2017, in China; the features of merging and splitting events are labeled in the third step. (b) Explanation of the STHW clustering process, including three steps: (1) identifying HWs; (2) spatially connecting HWs; and (3) temporally connecting HWs.
Remotesensing 14 04082 g001
Based on the selection of the temperature and data variables, the STHW types are separated into two groups: (1) daytime STHWs identified by the daily maximum temperature and (2) nighttime STHWs identified by the daily minimum temperature. Similarly, T a i r -based STHWs and H I -based STHWs are identified by the air temperature and heat index, respectively.
To date, these diverse HW definitions can be derived from different threshold calculation methods ( A B S , Q T L and Q T L m o v ), severity levels (mild, moderate and severe), types (daytime and nighttime), and data variables ( T a i r and H I ). The rest of this study is based on the use of control variable methods to compare these definitions.

2.5. STHW Characteristics

To evaluate STHW, STHW characteristics are defined separately at temporal and spatial scales. The temporal characteristics are all counted at the event scale (Table 2). Here, we focus only on the most impactful STHW event of the year, corresponding to the largest HW severity ( HWS ) each year. Because the number of small STHW events is much larger than that of typical HW events, tracking only the annual largest STHW event can prevent a large number of small events from offsetting the real temporal variants. Here, HW days ( HWD ) expresses the duration of STHW events, HW intensity ( HWI ) and HWS are used to describe the intensity and cumulative intensity, respectively, and HW area ( HWA ) is the occurrence area, which is further divided into three aspects: the mean, maximum, and total area.
In contrast to our focus only on the severest HW event temporally, all HW events that occur over the year are counted in the spatial analysis (Table 3). This is because the spatial HW characteristics should reflect the overall spatial HW distribution, while the largest annual event cannot mask the entire country. Here, HWD is calculated based on the number of days the HW occurs; HW frequency ( HWF ) is added to express the number of HW events; and HWI   and HWS are still calculated by obtaining the anomaly between the temperature and threshold.

3. Results

3.1. STHW Characteristics of Different Threshold Methods

Figure 2 shows the annual variations in the characteristics of the most severe STHW events that occurred in China from 1961–2017, defined by three different threshold methods: A B S , Q T L and Q T L m o v . The magnitudes of the STHW characteristics vary substantially among different threshold methods. For HWD , HWI , and HWS , the magnitudes obtained with the A B S method are significantly larger than those obtained with Q T L and Q T L m o v , with the rank of A B S > Q T L > Q T L m o v . The mean values obtained with A B S are 100.00 days, 6.33 °C, and 632.71 °C∙days for HWD , HWI , and HWS , respectively. In comparison, these mean values are only 58.00 days, 2.50 °C, and 145.82 °C·days, respectively, for Q T L and 18.75 days, 1.86 °C and 35.88 °C∙days, respectively, for Q T L m o v . For HWA avg and HWA max , the situation is reversed. The rank becomes Q T L > Q T L m o v > A B S , and the mean values of HWA avg are 2.38 × 106 km2, 1.58 × 106 km2 and 1.20 × 106 km2, respectively, while the mean HWA max values are 5.31 × 106 km2, 3.77 × 106 km2 and 2.51 × 106 km2, respectively, for those three threshold methods.
Regarding the trends, for HWD , the Q T L threshold exhibits an increasing trend of 7.42 days/decade, much higher than the trends obtained using the other methods (1.08 days/decade and 2.88 days/decade) but not significantly so. HWI   and HWS exhibit similar results. The A B S threshold always indicates the highest slope, reaching 0.24 °C/decade and 28.24 °C·day/decade. However, for HWI , only the A B S threshold shows significance, while all methods show significance for HWS . For all area-related characteristics, the Q T L threshold identifies the highest increasing rate. Especially for HWA sum , the Q T L threshold shows an amazing trend of 25.3 × 106 km2/decade, approximately 4–5 times higher than those obtained with the other methods. In comparison, the A B S threshold has the lowest increasing rates of 0.03 km2/decade, 0.03 km2/decade, and 4.04 km2/decade corresponding to HWA avg , HWA max , and HWA sum , respectively.
Regarding the spatial magnitude (Figure 3), except for HWF , the magnitudes obtained with Q T L and Q T L m o v are similar for all STHW characteristics. Because the Q T L method focuses only on the highest temperature period, while Q T L m o v   considers all periods in a year, Q T L m o v detects more STHWs than Q T L , resulting in a higher HWF (Figure 3b2,b3). Using the A B S method, the derived spatial distribution of the multiannual average STHW characteristics is substantially different from those obtained with Q T L and Q T L m o v . Most STHWs are detected in the North China Plain and South China, while no STHWs are detected on the Tibetan Plateau because of the low climatological temperatures there.
Regarding the spatial trend, the three analyzed methods show consistently increasing trends for all STHW spatial characteristics except HWI (Figure 4). The STHWs become longer in duration, more widespread in impact, and more intense in intensity and severity. Regarding HWD , HWA avg , HWA max and HWA sum , the ranking of the increase rates is Q T L > Q T L m o v > A B S . For HWD , HWI , and HWS , there is a prevalent and consistent increasing trend among Q T L , Q T L m o v , and A B S . However, for HWF , the increasing trend is not significant in most regions of China (Figure 4a2,b2) and even becomes negative in North China (Figure 4b2). This is mainly because mild and moderate STHWs gradually turn into severe STHWs under the warming climate conditions. The frequency of mild and moderate STHWs decreases, while that of severe STHWs increased (Figures S1 and S2). Additionally, for HWD , HWS , and HWI , noticeable decreasing trends arise in the middle of the North China Plain and across large parts of Northeast China.

3.2. HI-Based and T a i r -Based STHWs

Many researchers utilize T a i r   rather than H I to measure HWs, though this may produce errors when attempting to accurately reflect HWs’ impacts on humankind. Moreover, the results may differ to a relatively great extent due to the addition of relative humidity. Hence, we compare the difference between H I -based and T a i r -based STHWs.
Figure 5 presents the difference between the H I -based and T a i r -based moderate daytime STHWs from the perspective of the relative change ratio between the present period of 1991–2017 and the historical period of 1961–1990 in China. Regarding the Q T L m o v threshold, the relative change ratio of the H I -based STHW is approximately 20% higher than that of the T a i r -based STHW for all STHW characteristics (Figure 5). For the three threshold methods, the differences in the relative change ratio between the H I -based and T a i r -based STHWs are extremely small, especially when using the A B S threshold. Moreover, regarding the Q T L method, the result is even reversed; the relative change ratio of the H I -based STHW is lower than that of the T a i r -based STHW. For threshold methods that do not consider any seasonal variation (i.e., A B S and Q T L ), the spatial continuity of the H I -based exceedance is weaker than that of the T a i r -based exceedance, thus contributing to the weak and reversed phenomena observed when using the A B S and Q T L threshold methods.

3.3. STHWs with Different Severity Levels

Figure 6 shows the relative change ratios of STHW characteristics between the present and historical periods under different HW severity levels in China. As the severity level rises, most characteristics show increased growth trends together with significance. Specifically, under extreme severity levels, the HWA max values obtained with the Q T L m o v and Q T L threshold methods reach 0.51 and 0.53, respectively, with increases of nearly 50% and at least threefold compared to the mild and moderate levels (0.17 and 0.18 for the Q T L m o v threshold, 0.08 and 0.18 for the Q T L threshold). Generally, compared to mild STHWs, the average relative change ratios of moderate and severe STHWs among different STHW characteristics are approximately 1.2 and 2.4 times higher, respectively, indicating that extreme STHWs increase much faster than mild and moderate STHWs; this finding is consistent with the results of previous HW studies performed at the grid scale [22]. However, when using the A B S threshold method, these values are approximately 0.87 and 0.93 times smaller than those of mild STHWs, further confirming that the spatial continuity of the results obtained using the A B S threshold method is weak and may be unsuitable for STHW detection. Because most nonsignificant trends occur in the HWI results, we mainly focus on comparing the other five indexes. Among those five indexes, the relative changing ratio of HWA max is the highest, while that of HWA sum ranks second, indicating that the impacted spatial areas of extreme STHWs increase much faster than their duration, intensity, and severity.
Figure 3. Spatial distributions of the multiannual average daytime, moderate, H I -based STHW spatial characteristics from 1961–2017: the combination results of using (a) A B S , (b) Q T L and (c) Q T L m o v method and (1) HWD , (2) HWF , (3) HWI and (4) HWS characteristics. The frequency distribution is shown in the lower-left corner of each subplot.
Figure 3. Spatial distributions of the multiannual average daytime, moderate, H I -based STHW spatial characteristics from 1961–2017: the combination results of using (a) A B S , (b) Q T L and (c) Q T L m o v method and (1) HWD , (2) HWF , (3) HWI and (4) HWS characteristics. The frequency distribution is shown in the lower-left corner of each subplot.
Remotesensing 14 04082 g003
Figure 4. Spatial distribution of the annual trend of H I -based, moderate, daytime STHW characteristics from 1961–2017: the combination results of using (a) A B S , (b) Q T L and (c) Q T L m o v method and (1) HWD , (2) HWF , (3) HWI and (4) HWS characteristics. The filled color represents the value of Sen’s slope. The hatches denote trends that are significant at the 0.05 level. For each subplot, the percentages of positively significant (PS) and negatively significant (NS) trends are shown in the top-center, and a histogram of trends is plotted in the lower-left corner. The frequency distribution is shown in the lower-left corner of each subplot. In b4, the highlighted area is the focus area used to explain the negative trends in Section 4.1.
Figure 4. Spatial distribution of the annual trend of H I -based, moderate, daytime STHW characteristics from 1961–2017: the combination results of using (a) A B S , (b) Q T L and (c) Q T L m o v method and (1) HWD , (2) HWF , (3) HWI and (4) HWS characteristics. The filled color represents the value of Sen’s slope. The hatches denote trends that are significant at the 0.05 level. For each subplot, the percentages of positively significant (PS) and negatively significant (NS) trends are shown in the top-center, and a histogram of trends is plotted in the lower-left corner. The frequency distribution is shown in the lower-left corner of each subplot. In b4, the highlighted area is the focus area used to explain the negative trends in Section 4.1.
Remotesensing 14 04082 g004
Figure 5. The change ratios of temporal STHW characteristics between the present period of 1991–2017 and the historical period of 1961–1990: for (a) HWD , (b) HWI , (c) HWS , (d) HWA avg , (e) HWA max and (f) HWA sum . Daytime, severe STHWs are used here. The change ratio equals ( x n o w x h i s ) / x h i s , where x h i s and x n o w are the mean values of the considered STHW characteristics during the historical and present periods, respectively. The “*” symbol denotes that the Sen’s slopes of the STHW characteristics from 1961–2017 are significant at the 0.05 level.
Figure 5. The change ratios of temporal STHW characteristics between the present period of 1991–2017 and the historical period of 1961–1990: for (a) HWD , (b) HWI , (c) HWS , (d) HWA avg , (e) HWA max and (f) HWA sum . Daytime, severe STHWs are used here. The change ratio equals ( x n o w x h i s ) / x h i s , where x h i s and x n o w are the mean values of the considered STHW characteristics during the historical and present periods, respectively. The “*” symbol denotes that the Sen’s slopes of the STHW characteristics from 1961–2017 are significant at the 0.05 level.
Remotesensing 14 04082 g005
Figure 6. Same as Figure 5, but comparing different HW severity levels. The H I -based, daytime STHW results are used here.
Figure 6. Same as Figure 5, but comparing different HW severity levels. The H I -based, daytime STHW results are used here.
Remotesensing 14 04082 g006

3.4. Daytime and Nighttime STHWs

Figure 7 compares the differences in the relative change ratios of STHW characteristics between daytime and nighttime STHWs. Surprisingly, daytime STHWs show higher ratios only for HWI but lower ratios for HWD , HWS , HWA avg , HWA max , and HWA sum . For instance, for the HWI values identified by the Q T L threshold, the ratios corresponding to daytime and nighttime STHWs are 0.10 and 0.07, respectively. For the Q T L m o v threshold method, daytime STHWs possess at least twofold higher HWD (1.84 and 0.87), HWA avg   (0.36 and 0.16), HWA max   (0.56 and 0.18), and HWA sum   (3.11 and 1.22) values than nighttime STHWs. Overall, the results show that the average increase in nighttime STHWs is approximately 60% faster than that of daytime STHWs when using the Q T L threshold, and approximately 120% faster when using the Q T L m o v method.

4. Discussion

4.1. Explanation of Negative STHW Characteristic Trends in North China and Northeast China

Although prevalent increasing trends in spatial STHW characteristics (e.g., HWD , HWI , and HWS ) are observed in most regions of China, significant decreasing trends are presented in parts of North China and Northeast China (Figure 4). Taking one representative region, North China, as an example, we further illustrate the underlying reasons for these decreasing trends (Figure 8 and Figure 9).
We found that negative STHW trends can be attributed to the following two aspects: (1) On the one hand, temperature warming is spatially asymmetric. Although increasing temperatures are prevalent in most regions of the globe [37], some regions are experiencing the inverse trend, with significant decreasing temperature trends. North China is one of these cases [38,39]. Undoubtedly, the decreasing local temperatures contribute to the decreasing STHW characteristics at the spatial level. (2) On the other hand, temperature warming is asymmetric over time. As shown in Figure 9, strong increasing temperature trends are observed in spring and autumn; a weak increasing trend is observed in winter; and reverse decreasing trends are observed in summer and at the end of spring. We also found that decreasing R H primarily contributes to the decreasing H I in summer (Figure 8) and, hence, to the asymmetric warming observed among different methods. This phenomenon directly leads to the decreasing trends of Q T L threshold-based STHW characteristics, as the Q T L threshold method mainly focuses on the hottest period each year (namely, summer), when significant decreasing trends occur. Moreover, regarding the Q T L m o v threshold-based STHWs, this phenomenon weakens the increasing STHW trends. However, whether the trend direction becomes negative depends on the increasing STHW magnitudes in other seasons, as Q T L m o v   focuses on H I anomalies throughout the whole year.
Figure 9 further investigates the offsetting effect of summertime cooling and warming in other seasons regarding the Q T L m o v   -based STHW results in North China. The figure shows that although the annual averaged H I and T a i r have positive trends, the 90%, 95% and 99% quantiles of their Q T L m o v   -based anomalies in the current period from 1991–2017 apparently decrease compared to those in the historical period from 1961–1990 (Figure 9). This indicates that summertime cooling has a stronger impact on the negative STHW trends in the spatial component than warming in other seasons in North China and explains the negative STHW characteristics trends derived in North China.

4.2. Possible Reasons for the STHW Characteristic Differences among Different STHW Definitions

There is evidence that H I increases faster than T a i r regardless of the area-weighted original value or area-weighted anomaly (Figure 10a1–d1). This result is consistent with previous findings that the joint impact of reduced R H and increased temperatures leads to more unprecedented HWs [3,13,36]. However, only for the Q T L m o v threshold do the H I -based STHW characteristics always present higher ratios than the T a i r -based characteristics (Figure 5). This result implies that the Q T L m o v threshold is more suitable for STHW detection. For threshold methods that do not consider seasonal variations (i.e., A B S and Q T L ), the spatial continuity of the H I exceedance is relatively weak, thus resulting in the increase rates of STHW characteristics being underestimated.
Regarding different HW severity levels, the averaged relative change ratios of moderate and severe STHWs among different STHW characteristics are approximately 1.2 and 2.4 times higher than those of the mild STHW for the Q T L m o v and Q T L threshold methods, respectively, implying that extreme STHWs will increase much faster. Due to relatively the low numbers of STHW events identified as the severity level increases. The occurrence of fewer events will cause the extreme values to be more prominent, and more STHW events may “average” the results, resulting in relatively low change ratios.
Regarding daytime and nighttime STHWs, the increase in nighttime STHWs is approximately 60% faster than that in daytime STHWs using the Q T L threshold and approximately 120% faster using the Q T L m o v method. In addition, the latter method usually results in a higher change ratio than the former. This finding is consistent with previous studies [17,40,41] and is mainly the result of intensified heating at nighttime compared to daytime (0.46 °C∙decade−1 at nighttime and 0.42 °C∙decade−1 at daytime) [17,42,43], further leading to an increasing tendency for nighttime HWs to replace the previously dominating independent daytime events [17]. As a further explanation, Luo [44] revealed that increased cloud cover can reduce the solar radiation received by the land surface, thus potentially increasing the downward longwave radiation emitted by the atmosphere and clouds at the surface at night and leading to reduced cooling of longwave radiation and increased nighttime temperatures. This may be the root cause of such rapid increases in nighttime HWs.

4.3. Significance and Limitations

Recently, Luo [25] proposed a 3D perspective to analyze the spatiotemporal patterns of HWs. In his result, the characteristic trends are 1.14 days/decade for the HW lifetime, 0.04 °C/decade for the HW intensity, 0.18 × 106 km2/decade for the maximum area and 1.30 × 106 km2/decade for the total area, while the corresponding metrics obtained in our study are 2.88 days/decade for HWD , 0.06 °C/decade for HWI , 0.19 × 106 km2/decade for HWA max and 5.77 × 106 km2/decade for HWA sum . Compared to his approach, we focused on the most impactful STHW events regarding the temporal characteristics, resulting in relatively high increasing growth in HWD and HWA sum . Apart from these indexes, our results show a high degree of consistency. Despite the divergences among the analyzed methods, both results showed significantly positive trends, demonstrating the robustness of our results.
The study offers some important insights, outlined as follows:
(1)
An improved algorithm is used to detect STHWs. Instead of simply merging all temporally and spatially adjacent HW grids into a single STHW event, we set two thresholds: n c o n n e c t and n o v e r l a p . Compared to the ordinary 3D connection algorithm, this method considers the merging and splitting behaviors of STHW events during their development and diminution.
(2)
Past articles rarely used spatiotemporal clustering algorithms to compare the spatiotemporal patterns under different HW definitions. The majority of investigations of spatial or temporal variations focused only on the grid scale and ignored joint spatial and temporal evolution pattern. In addition, H I -based HWs are closely relevant to human perception, and while most past studies used either T a i r or H I , the difference among these variables remains largely unexplored. In our study, we compare these HW definitions as comprehensively as possible.
Due to practical constraints, this paper cannot provide a comprehensive review of the following considerations:
(1)
The data selection method may influence the results, especially for R H . According to our early study, using unhomogenized R H data can produce entirely different results. To ensure the objectivity of the results, we replaced the original R H with a homogenized R H . However, the differences introduced by different homogenization methods were not assessed. It is thus better to use multiple datasets to demonstrate the robustness of the results. In addition, the selection of temperature data may produce differences, as expressed in the trend magnitudes and the different distributions of spatial characteristics, but not to the same extent as the R H selection.
(2)
Different selections of the n c o n n e c t and n o v e r l a p thresholds were not investigated herein. Comparing the similarity of the results derived with different thresholds may enhance the robustness of our result.
Overall, our study can provide a more comprehensive understanding of the change patterns of STHW events in China and help researchers select a suitable HW definition in China in future research. However, the influence of the chosen data and n c o n n e c t and n o v e r l a p values remains unknown in this work, and this will be an important issue in future research.

5. Conclusions

In this study, we used a modern spatiotemporal clustering algorithm to detect spatiotemporally contiguous heatwaves (STHWs) from a 3-D perspective in China from 1961–2017. We systematically investigated the influence of various threshold methods ( A B S , Q T L , and Q T L m o v ), variables ( T a i r and H I ), severity levels (mild, moderate, and severe), and time (daytime and nighttime) selections on the STHW change characteristics. The results suggest the following:
(1)
Although the trends’ magnitudes were divergent, all threshold methods showed a consistent phenomenon in most regions of China: STHWs have become longer-lasting, more severe, and have impacted larger areas spatially. However, the HWI increasing trends identified in most cases are insignificant. Moreover, the HWF of the moderate STHWs has decreased mainly because the transfer of mild and moderate STHWs to extreme STHWs decreased the frequency of the former events. In small parts of central-eastern and northeast China, some negative STHW trends were observed due to asymmetric warming in these areas, i.e., these regions experience negative H I trends in summer but positive trends in other seasons. For threshold methods that do not consider any seasonal variation (i.e., A B S and Q T L ), the spatial continuity of the H I exceedance is relatively weak, resulting in the increase rates of STHW characteristics being underestimated. In addition, the A B S method uses eastern China to calculate the threshold, and it is difficult to identify HWs for the lower-elevation areas in western China (e.g., the Tibetan Plateau). In comparison, Q T L m o v performs much better and is more suitable for use across China.
(2)
Regarding H I -based and T a i r -based STHWs, when using the Q T L m o v threshold, the relative change ratios of all H I -based STHW characteristics are approximately 20% higher than the T a i r -based STHW characteristics, suggesting that the contribution of R H to HW variability is visible and causes human-perceived HWs to increase rapidly, as hot temperatures are even more unbearable as the presence of humidity in the environment reduces the ability of the human body to cool itself. Therefore, H I -based STHW should be priorly used. T a i r   alone without considering R H may underestimate the increasing STHW magnitudes. However, this phenomenon is weak when considering the Q T L threshold and is even reversed when using the A B S method due to the weak spatial continuity of these two methods.
(3)
Regarding different HW severity levels, the averaged relative change ratios of moderate and severe STHWs among different STHW characteristics are approximately 1.2 and 2.4 times higher than those of mild STHWs when using Q T L m o v and Q T L threshold methods, respectively, implying that extreme STHWs will increase much faster than mild or moderate STHWs; this finding is inextricably related to climate warming. Therefore, this phenomenon requires preparation and mitigation programs to avoid or mitigate the effects of the rapid growth of extreme weather events.
(4)
Regarding daytime and nighttime STHWs, the increase in nighttime STHWs is approximately 60% faster than that of daytime STHWs when applying the Q T L threshold and approximately 120% faster when applying the Q T L m o v method. Increasing nighttime exposures may increase the risk of heat-related illnesses.
Our study provides a systematic investigation of different STHW definition methods and will benefit HW studies in the future.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/rs14164082/s1 and https://zenodo.org/record/7012026. Figure S1. Same as Figure 3, but the H I -based, daytime, mild STHW results are used here. Figure S2. Same as Figure 3, but the H I -based, daytime, severe STHW results are used here.

Author Contributions

Conceptualization, D.K.; software, H.S. and D.K.; data curation, L.X.; writing—original draft preparation, H.S. and D.K.; writing—review and editing, H.S., D.K., X.G. and J.L.; visualization, H.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Natural Science Foundation of Hubei Province (Grant 2021CFB198), the National Natural Science Foundation of China (Grants 4210011820, 41901041, and 42001042), the Fundamental Research Funds for the Central Universities, China University of Geosciences (Wuhan) (Grant CUG2106107).

Data Availability Statement

The CN05.1 dataset can be available at http://ccrc.iap.ac.cn/resource/detail?id=228 (accessed on 28 March 2022). The homogenized R H data is available at https://www.scidb.cn/en/detail?dataSetId=633694461334913025 (accessed on 5 April 2022).

Acknowledgments

The authors thank Wu, J. for providing CN05.1 and Li, Z. for providing homogenized R H data.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this paper:
HWHeatwave
STHWSpatiotemporally continuous heatwave
A B S Absolute threshold method
Q T L Quantile-based threshold method
Q T L m o v Moving quantile-based threshold method
T a i r Air temperature
T m a x Daily maximum air temperature
T m i n Daily minimum air temperature
R H Relative humidity
H I Heat index
d o y Day of year
HWD Heatwave duration
HWI Heatwave intensity
HWS Heatwave severity
HWA avg Heatwave mean occurrence area
HWA max Heatwave maximum occurrence area
HWA sum Heatwave total occurrence area
HWF Heatwave frequency

References

  1. Li, Y.; Ding, Y.; Li, W. Observed Trends in Various Aspects of Compound Heat Waves across China from 1961 to 2015. J. Meteorol. Res. 2017, 31, 455–467. [Google Scholar] [CrossRef]
  2. Meehl, G.A.; Tebaldi, C. More Intense, More Frequent, and Longer Lasting Heat Waves in the 21st Century. Science 2004, 305, 994–997. [Google Scholar] [CrossRef] [PubMed]
  3. Fischer, E.M.; Schär, C. Consistent Geographical Patterns of Changes in High-Impact European Heatwaves. Nat. Geosci. 2010, 3, 398–403. [Google Scholar] [CrossRef]
  4. Habeeb, D.; Vargo, J.; Stone, B. Rising Heat Wave Trends in Large US Cities. Nat. Hazards 2015, 76, 1651–1665. [Google Scholar] [CrossRef]
  5. Kenney, W.L.; Craighead, D.H.; Alexander, L.M. Heat Waves, Aging, and Human Cardiovascular Health. Med. Sci. Sports Exerc. 2014, 46, 1891–1899. [Google Scholar] [CrossRef]
  6. Mora, C.; Counsell, C.W.W.; Bielecki, C.R.; Louis, L.V. Twenty-Seven Ways a Heat Wave Can Kill You. Circ. Cardiovasc. Qual. Outcomes 2017, 10, e004233. [Google Scholar] [CrossRef]
  7. Yin, Q.; Wang, J. The Association between Consecutive Days’ Heat Wave and Cardiovascular Disease Mortality in Beijing, China. BMC Public Health 2017, 17, 223. [Google Scholar] [CrossRef]
  8. Carlson, C.J.; Albery, G.F.; Merow, C.; Trisos, C.H.; Zipfel, C.M.; Eskew, E.A.; Olival, K.J.; Ross, N.; Bansal, S. Climate Change Increases Cross-Species Viral Transmission Risk. Nature 2022, 607, 555–562. [Google Scholar] [CrossRef]
  9. Chung, U.; Gbegbelegbe, S.; Shiferaw, B.; Robertson, R.; Yun, J.I.; Tesfaye, K.; Hoogenboom, G.; Sonder, K. Modeling the Effect of a Heat Wave on Maize Production in the USA and Its Implications on Food Security in the Developing World. Weather Clim. Extrem. 2014, 5–6, 67–77. [Google Scholar] [CrossRef]
  10. Yuan, W.; Cai, W.; Chen, Y.; Liu, S.; Dong, W.; Zhang, H.; Yu, G.; Chen, Z.; He, H.; Guo, W.; et al. Severe Summer Heatwave and Drought Strongly Reduced Carbon Uptake in Southern China. Sci. Rep. 2016, 6, 18813. [Google Scholar] [CrossRef]
  11. Perkins, S.E. A Review on the Scientific Understanding of Heatwaves-Their Measurement, Driving Mechanisms, and Changes at the Global Scale. Atmos. Res. 2015, 164–165, 242–267. [Google Scholar] [CrossRef]
  12. Vogel, M.M.; Zscheischler, J.; Fischer, E.M.; Seneviratne, S.I. Development of Future Heatwaves for Different Hazard Thresholds. J. Geophys. Res. Atmos. 2020, 125, e2019JD032070. [Google Scholar] [CrossRef]
  13. You, Q.; Jiang, Z.; Kong, L.; Wu, Z.; Bao, Y.; Kang, S.; Pepin, N. A Comparison of Heat Wave Climatologies and Trends in China Based on Multiple Definitions. Clim. Dyn. 2016, 48, 3975–3989. [Google Scholar] [CrossRef]
  14. Brooke Anderson, G.; Bell, M.L. Heat Waves in the United States: Mortality Risk during Heat Waves and Effect Modification by Heat Wave Characteristics in 43 U.S. Communities. Environ. Health Perspect. 2011, 119, 210–218. [Google Scholar] [CrossRef] [PubMed]
  15. Su, Q.; Dong, B.W. Recent Decadal Changes in Heat Waves over China: Drivers and Mechanisms. J. Clim. 2019, 32, 4215–4234. [Google Scholar] [CrossRef]
  16. Wang, J.; Chen, Y.; Tett, S.F.B.; Yan, Z.; Zhai, P.; Feng, J.; Xia, J. Anthropogenically-Driven Increases in the Risks of Summertime Compound Hot Extremes. Nat. Commun. 2020, 11, 582. [Google Scholar] [CrossRef]
  17. Chen, Y.; Li, Y. An Inter-Comparison of Three Heat Wave Types in China during 1961–2010: Observed Basic Features and Linear Trends. Sci. Rep. 2017, 7, srep45619. [Google Scholar] [CrossRef]
  18. Li, J.; Chen, Y.D.; Gan, T.Y.; Lau, N.-C. Elevated Increases in Human-Perceived Temperature under Climate Warming. Nat. Clim. Chang. 2018, 8, 43–47. [Google Scholar] [CrossRef]
  19. Mora, C.; Dousset, B.; Caldwell, I.R.; Powell, F.E.; Geronimo, R.C.; Bielecki, C.R.; Counsell, C.W.W.; Dietrich, B.S.; Johnston, E.T.; Louis, L.V.; et al. Global Risk of Deadly Heat. Nat. Clim. Chang. 2017, 7, 501–506. [Google Scholar] [CrossRef]
  20. Russo, S.; Sillmann, J.; Sterl, A. Humid Heat Waves at Different Warming Levels. Sci. Rep. 2017, 7, 7477. [Google Scholar] [CrossRef]
  21. Luo, M.; Lau, N. Increasing Heat Stress in Urban Areas of Eastern China: Acceleration by Urbanization. Geophys. Res. Lett. 2018, 45, 13060–13069. [Google Scholar] [CrossRef]
  22. Kong, D.; Gu, X.; Li, J.; Ren, G.; Liu, J. Contributions of Global Warming and Urbanization to the Intensification of Human-Perceived Heatwaves Over China. J. Geophys. Res. Atmos. 2020, 125, e2019JD032175. [Google Scholar] [CrossRef]
  23. Diaz, V.; Corzo Perez, G.A.; Van Lanen, H.A.J.; Solomatine, D.; Varouchakis, E.A. An Approach to Characterise Spatio-Temporal Drought Dynamics. Adv. Water Resour. 2020, 137, 103512. [Google Scholar] [CrossRef]
  24. Samaniego, L.; Kumar, R.; Zink, M. Implications of Parameter Uncertainty on Soil Moisture Drought Analysis in Germany. J. Hydrometeorol. 2013, 14, 47–68. [Google Scholar] [CrossRef]
  25. Luo, M.; Lau, N.; Liu, Z.; Wu, S.; Wang, X. An Observational Investigation of Spatiotemporally Contiguous Heatwaves in China From a 3D Perspective. Geophys. Res. Lett. 2022, 49, e2022GL097714. [Google Scholar] [CrossRef]
  26. Lyon, B.; Barnston, A.G.; Coffel, E.; Horton, R.M. Projected Increase in the Spatial Extent of Contiguous US Summer Heat Waves and Associated Attributes. Environ. Res. Lett. 2019, 14, 114029. [Google Scholar] [CrossRef]
  27. Lo, S.H.; Chen, C.T.; Russo, S.; Huang, W.R.; Shih, M.F. Tracking Heatwave Extremes from an Event Perspective. Weather Clim. Extrem. 2021, 34, 100371. [Google Scholar] [CrossRef]
  28. Wu, J.; Gao, X.-J. A Gridded Daily Observation Dataset over China Region and Comparison with the Other Datasets. Chin. J. Geophys. 2013, 56, 1102–1111. [Google Scholar]
  29. Freychet, N.; Tett, S.F.B.; Yan, Z.; Li, Z. Underestimated Change of Wet-Bulb Temperatures over East and South China. Geophys. Res. Lett. 2020, 47, e2019GL086140. [Google Scholar] [CrossRef]
  30. Li, Z.; Yan, Z.; Zhu, Y.; Freychet, N.; Tett, S. Homogenized Daily Relative Humidity Series in China during 1960–2017. Adv. Atmos. Sci. 2020, 37, 318–327. [Google Scholar] [CrossRef]
  31. De Bono, A.; Peduzzi, P.; Kluser, S.; Giuliani, G. Impacts of Summer 2003 Heat Wave in Europe. Environ. Alert Bull. 2004, 2, 4. [Google Scholar]
  32. Steadman, R.G. A Universal Scale of Apparent Temperature. J. Clim. Appl. Meteorol. 1984, 23, 1674–1687. [Google Scholar] [CrossRef]
  33. Rothfusz, L.P.; NWS Southern Region Headquarters. The Heat Index Equation (or, More than You Ever Wanted to Know about Heat Index). Atmos. Clim. Sci. 1990, 4, 9023. [Google Scholar]
  34. Qi, L.; Wang, Y. Changes in the Observed Trends in Extreme Temperatures over China around 1990. J. Clim. 2012, 25, 5208–5222. [Google Scholar] [CrossRef]
  35. Lloyd-Hughes, B. A Spatio-Temporal Structure-Based Approach to Drought Characterisation. Int. J. Climatol. 2012, 32, 406–418. [Google Scholar] [CrossRef]
  36. Van der Walt, S.; Schönberger, J.L.; Nunez-Iglesias, J.; Boulogne, F.; Warner, J.D.; Yager, N.; Gouillart, E.; Yu, T. Scikit-Image: Image Processing in Python. PeerJ 2014, 2014, e453. [Google Scholar] [CrossRef]
  37. IPCC Summary for Policymakers. In Climate Change 2022: Mitigation of Climate Change. Contribution of Working Group III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Shukla, P.R.; Skea, J.; Slade, R.; al Khourdajie, A.; van Diemen, R.; McCollum, D.; Pathak, M.; Some, S.; Vyas, P.; Fradera, R.; et al. (Eds.) Cambridge University Press: Cambridge, UK; New York, NY, USA, 2022. [Google Scholar]
  38. Fu, D.; Ding, Y. The Study of Changing Characteristics of the Winter Temperature and Extreme Cold Events in China over the Past Six Decades. Int. J. Climatol. 2021, 41, 2480–2494. [Google Scholar] [CrossRef]
  39. Qian, W.; Lin, X. Regional Trends in Recent Temperature Indices in China. Clim. Res. 2004, 27, 119–134. [Google Scholar] [CrossRef]
  40. Zhang, Y.; Mao, G.; Chen, C.; Lu, Z.; Luo, Z.; Zhou, W. Population Exposure to Concurrent Daytime and Nighttime Heatwaves in Huai River Basin, China. Sustain. Cities Soc. 2020, 61, 102309. [Google Scholar] [CrossRef]
  41. An, N.; Zuo, Z. Changing Structures of Summertime Heatwaves over China during 1961–2017. Sci. China Earth Sci. 2021, 64, 1242–1253. [Google Scholar] [CrossRef]
  42. Du, Z.; Zhao, J.; Liu, X.; Wu, Z.; Zhang, H. Recent Asymmetric Warming Trends of Daytime versus Nighttime and Their Linkages with Vegetation Greenness in Temperate China. Environ. Sci. Pollut. Res. 2019, 26, 35717–35727. [Google Scholar] [CrossRef] [PubMed]
  43. Song, Z.; Yang, H.; Huang, X.; Yu, W.; Huang, J.; Ma, M. The Spatiotemporal Pattern and Influencing Factors of Land Surface Temperature Change in China from 2003 to 2019. Int. J. Appl. Earth Obs. Geoinf. 2021, 104, 102537. [Google Scholar] [CrossRef]
  44. Luo, M.; Lau, N.-C.; Liu, Z. Different Mechanisms for Daytime, Nighttime, and Compound Heatwaves in Southern China. Weather Clim. Extrem. 2022, 36, 100449. [Google Scholar] [CrossRef]
Figure 2. Annual variations in daytime H I -based STHW temporal characteristics from 1961–2017: for (a) HWD , (b) HWI , (c) HWS , (d) HWA avg , (e) HWA max and (f) HWA sum . The STHW characteristics were calculated at the event scale, and only the most serious annual STHW event was selected. The dashed lines represent the annual variations in STHW characteristics; the solid lines represent the corresponding Sen’s slope trends. In each subplot, Sen’s slope for each threshold method is labeled on the top right. The symbols “*”, “**” and “***” indicate that slopes are significant at the 0.05, 0.01 and 0.001 levels, respectively, while “ns” represents slopes that are not significant.
Figure 2. Annual variations in daytime H I -based STHW temporal characteristics from 1961–2017: for (a) HWD , (b) HWI , (c) HWS , (d) HWA avg , (e) HWA max and (f) HWA sum . The STHW characteristics were calculated at the event scale, and only the most serious annual STHW event was selected. The dashed lines represent the annual variations in STHW characteristics; the solid lines represent the corresponding Sen’s slope trends. In each subplot, Sen’s slope for each threshold method is labeled on the top right. The symbols “*”, “**” and “***” indicate that slopes are significant at the 0.05, 0.01 and 0.001 levels, respectively, while “ns” represents slopes that are not significant.
Remotesensing 14 04082 g002
Figure 7. Same as Figure 5 but comparing daytime and nighttime HWs. The H I -based, severe STHW results are used here.
Figure 7. Same as Figure 5 but comparing daytime and nighttime HWs. The H I -based, severe STHW results are used here.
Remotesensing 14 04082 g007
Figure 8. Sen’s slope of H I from 1961–2017 in North China, exhibiting one representative region with significant decreasing STHW characteristics and as depicted in Figure 4b4. For each month, the H I values are averaged from all grid cells in the analyzed region. The lower and upper whiskers correspond to the 10th and 90th quantiles, respectively, and outliers are removed. β on the top right is the annual H I trend averaged across all grid cells.
Figure 8. Sen’s slope of H I from 1961–2017 in North China, exhibiting one representative region with significant decreasing STHW characteristics and as depicted in Figure 4b4. For each month, the H I values are averaged from all grid cells in the analyzed region. The lower and upper whiskers correspond to the 10th and 90th quantiles, respectively, and outliers are removed. β on the top right is the annual H I trend averaged across all grid cells.
Remotesensing 14 04082 g008
Figure 9. The empirical cumulative distribution function (ECDF) results of the Q T L threshold-based anomalies of (a) H I and (b) T a i r and the Q T L m o v threshold-based anomalies of (c) H I and (d) T a i r   in North China, a representative region with significantly decreasing STHW characteristics as depicted in Figure 4b4. In each subplot, the right panel shows a zoomed image of the gray window in the left panel.
Figure 9. The empirical cumulative distribution function (ECDF) results of the Q T L threshold-based anomalies of (a) H I and (b) T a i r and the Q T L m o v threshold-based anomalies of (c) H I and (d) T a i r   in North China, a representative region with significantly decreasing STHW characteristics as depicted in Figure 4b4. In each subplot, the right panel shows a zoomed image of the gray window in the left panel.
Remotesensing 14 04082 g009
Figure 10. The seasonal (left panel) and monthly (right panel) Sen’s slope values of the area-weighted original temperature variables (a1,a2) and different threshold-based area-weighted anomalies (b1–d2) from 1991–2017 in China. The annual trend β and significance level are presented on the top in each subplot of the right panel. The symbols “**” and “***” indicate that the slopes are significant at the 0.01 and 0.001 levels, respectively.
Figure 10. The seasonal (left panel) and monthly (right panel) Sen’s slope values of the area-weighted original temperature variables (a1,a2) and different threshold-based area-weighted anomalies (b1–d2) from 1991–2017 in China. The annual trend β and significance level are presented on the top in each subplot of the right panel. The symbols “**” and “***” indicate that the slopes are significant at the 0.01 and 0.001 levels, respectively.
Remotesensing 14 04082 g010
Table 1. The features of commonly used thresholds to define HWs.
Table 1. The features of commonly used thresholds to define HWs.
ThresholdMathematicsSpatially VariableTemporally VariableDescription
A B S CustomNoNoCustom defined, the thresholds are static both with different grids and on different dates
Q T L QuantileYesNoCalculated by quantile, the thresholds vary spatially, but are temporally static
Q T L m o v QuantileYesYesCalculated by rolling quantile, the thresholds vary both spatially and temporally
Table 2. Temporal metrics used to characterize STHW events. In each STHW event, t 0 + 1 is the start date, t 0 + d is the end date, and k refers to each grid cell.
Table 2. Temporal metrics used to characterize STHW events. In each STHW event, t 0 + 1 is the start date, t 0 + d is the end date, and k refers to each grid cell.
MetricEquationUnitDescription
Duration HWD = t = t 0 + 1 t = t 0 + d 1 dayFor an STHW event, the number of days from onset to demise
Intensity HWI = 1 d t = t 0 + 1 t = t 0 + d k S k , t ( T k , t T k t r s ) k S k , t °CFor an STHW event, the average daily area-weighted average anomaly over all grid cells
Severity HWS = t = t 0 + 1 t = t 0 + d k k S k , t ( T k , t T k t r s ) k S k , t °C·dayFor an STHW event, the aggregate daily area-weighted average anomaly over all grid cells
Mean Area HWA avg = t = t 0 + 1 t = t 0 + d k S k , t t = t 0 + 1 t = t 0 + d t 106 km2For an STHW event, the average daily areas over all grid cells
Maximum Area HWA max = max ( k S k , t 0 + 1 , , k S k , t 0 + d ) 106 km2For an STHW event, the maximum daily areas over all grid cells
Total Area HWA sum = t = t 0 + 1 t = t 0 + d k S k , t 106 km2For an STHW event, the aggregate daily areas over all grid cells
Table 3. Spatial metrics used to characterize STHW events.
Table 3. Spatial metrics used to characterize STHW events.
MetricEquationUnitDescription
Duration HWD = 1 n i j = 1 n i t i , j T i , j 1 dayFor each grid i , the number of days T i , j over which STHW event j occurred; n is the number of HW events that take place annually
Frequency HWF = n i timeFor each grid i , the number of STHW events that take place annually
Intensity HWI = 1 n i j = 1 n i 1 T i , j t j = 0 T i , j a j , t j °CFor each grid, i , the average mean anomaly across all STHW events; n is the number of STHW events that take place annually
Severity HWS = 1 n i j = 1 n i t j = 0 T i , j a j , t j °C·dayFor each grid i , the average total anomaly across all STHW events; n is the number of STHW events that take place annually
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Song, H.; Kong, D.; Xiong, L.; Gu, X.; Liu, J. Inter-Comparison of Diverse Heatwave Definitions in the Analysis of Spatiotemporally Contiguous Heatwave Events over China. Remote Sens. 2022, 14, 4082. https://doi.org/10.3390/rs14164082

AMA Style

Song H, Kong D, Xiong L, Gu X, Liu J. Inter-Comparison of Diverse Heatwave Definitions in the Analysis of Spatiotemporally Contiguous Heatwave Events over China. Remote Sensing. 2022; 14(16):4082. https://doi.org/10.3390/rs14164082

Chicago/Turabian Style

Song, Heyang, Dongdong Kong, Li Xiong, Xihui Gu, and Jianyu Liu. 2022. "Inter-Comparison of Diverse Heatwave Definitions in the Analysis of Spatiotemporally Contiguous Heatwave Events over China" Remote Sensing 14, no. 16: 4082. https://doi.org/10.3390/rs14164082

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop