**3. Results**

#### *3.1. Long-Term Soil Moisture–Temperature Coupling*

To know the spatial distribution of strong coupling between soil moisture and temperature in China, we first calculated the metric (Π) by using ET (evapotranspiration) and PET (potential evapotranspiration) from the GLEAM data, and T (2-meter temperature) and Rn (surface net radiation) from the ERA-Interim reanalysis data over the period 1980–2013. Figure 2 illustrates the derived long-term soil moisture–temperature coupling annually and in different seasons.

Figure 2a shows the Π values derived from the annual data. The coupling strength appears to be highest in Northeast China and part of the Tibetan Plateau, where the climate is neither too wet nor too dry, showing that soil moisture has the strongest impact on temperature over the annual average in these regions. The results are consistent with previous studies which have suggested that such hot spots of soil moisture–temperature coupling occur most in transitional regions between wet and dry climate [9,39,40]. In spring, there appears to be strong soil moisture–temperature coupling (Figure 2b) in large areas of China, including Northeast, North, Northwest China and the Tibetan Plateau, and Yunnan province as well.

In summer, the coupling strengths are largely reduced compared to those in spring. There are only coupling signals in North China, especially Inner Mongolia. In autumn, the strongest coupling signals appear in Northeast China and the northern part of the Tibetan Plateau, although these are also significantly weakened compared to those in spring. It is not surprising to find strong signals in Northeast China, where it is neither too wet nor too dry in summer. In the cold season of winter, quite limited coupling signals are found, indicating that the soil moisture generally has quite limited impact on temperature across China.

Over the whole region, as shown in Figure 2, soil moisture–temperature coupling is relatively stronger in spring, followed by summer and autumn, and rather insignificant in winter. The seasonality of the soil moisture–temperature coupling strength has a distinct regional variation. Over Northern China, which is mainly a arid/semi-arid region, the contribution of soil moisture to evapotranspiration is mainly limited by water. Here, the results in Figure 2 depict strong coupling, which suggests a stronger impact of soil moisture anomies on temperature. Given that Southern China is mainly a humid/semi-humid region, the comparatively weaker soil moisture–temperature coupling (as seen in Figure 2) demonstrates that the region is mostly dominated by energy-limited conditions [8,13]. The north of China is here identified as a hot spot in spring, mainly because dry conditions (water-limited) persist throughout the year. In spring, a sufficient energy supply for ET, due to melting ice and snow, increases the amount of water in these regions, as well as causes more evaporation and, thus, increased coupling strength to the atmosphere. This may explain why coupling strength is higher in spring than in summer in North China [41]. The seasonal transition from winter to spring affects the soil moisture thawing and radiation budget over the Tibetan Plateau, which results in more heat transfer into the atmosphere. The heat energy transferred to the atmosphere is used to warm the air, thus showing strong soil moisture–temperature coupling [42]. In Yunnan province, the water supply is not sustainable because of the special climatic conditions, thus the soil is drier during the spring, which leads to a strong soil moisture–temperature coupling [43].

Figure 3 shows the soil moisture–temperature coupling strengths in different seasons based on the four climatic regions. Obviously, the coupling strengths vary greatly in different regions and different seasons. Except for the humid region, the coupling strengths of the other three climatic regions are all the strongest in spring and the weakest in winter. From the perspective of different seasons, the strength differences between the arid region, the transitional region, and the Tibetan Plateau region are relatively small. The seasonal transition from winter to spring influences soil moisture thawing and radiation budgets, with more heat energy being transferred into atmosphere [42,43]. For the typical climate regions divided according to Figure 1, the arid and transitional regions and the Tibetan Plateau zone have analogous soil moisture and atmosphere conditions, thus they are very similar, particularly in spring. However, coupling strengths in the humid region appear to be much smaller than the other regions in all seasons. The transitional region appears to have the strongest coupling strength in terms of the annual average, followed by the arid region and the Tibetan Plateau. The soil moisture–temperature coupling in spring appears to be much stronger than in the other seasons, which is particularly significant in the arid and transitional regions and the Tibetan Plateau. In the cold winter, the soil moisture–temperature coupling strengths are rather weak in all of the four climatic regions. Despite notable differences in the soil moisture–temperature coupling strengths among the four seasons, there are still some features in common. There appears to be hot spots of soil moisture–temperature coupling in the transitional zones between wet and dry regimes in all seasons.

**Figure 2.** Soil moisture–temperature coupling over China during the period 1980–2013. (**a**) Whole year; (**b**) spring (MAM, from March to May); (**c**) summer (JJA, from June to August); (**d**) autumn (SON, from September to November); (**e**) winter (DJF, from December to February).

**Figure 3.** Soil moisture–temperature coupling strengths in different seasons in the typical climate regions over China during the period 1980–2013.

In order to understand how the soil moisture–temperature coupling strengths are linked to the soil moisture amount, Figure 4 shows the density scatter plot of coupling strengths against the soil moisture amount in spring. It appears that, principally, the coupling strengths are linearly related to the soil moisture amounts across China; when soil increases, the strength of soil moisture-temperature coupling decreases, and this linear relationship is particularly clear when the soil moisture amount is more than 0.2 m<sup>3</sup>/m3. When soil is too dry, the coupling is not sensitive to soil moisture amount, and where there is too much soil moisture, especially in the humid regions, the density of data points is also high and shows that the coupling strengths are rather low. This result is also consistent with other studies [11,13].

**Figure 4.** The density scatter plot of coupling strengths against the soil moisture amount in spring.

#### *3.2. Coupling Anomalies in Heatwaves*

As seen in Figure 2, the soil moisture–temperature coupling is strongest in spring than in other seasons in China; however, many studies have reported heat wave amplifications through the feedback loops between soil moisture deficit and temperature in summer [18,44]. To understand the detailed processes of soil moisture–temperature coupling in the heatwave events, we conducted two case studies on daily scales.

#### 3.2.1. Case 1: Heatwave of Southeast China in Summer 2013

In the summer of 2013, Southeast China experienced abnormally high temperatures, which broke the heat records for the past 141 years and led to an unprecedented heatwave across China [1]. This unprecedented anomalies reached high values from 23 July to 14 August [3]. This disaster caused about US\$10 billion in crop damage, and a total of 5758 Heatwave-related illness cases were reported [40–47].

Figure 5a illustrates the soil moisture–temperature coupling from 1 June to 30 August, when the heatwave occurred. Figure 5b,c shows the soil moisture and temperature anomalies, referring to the multi-year average of 1980–2013. It appears that the summer temperature was strongly coupled to land surface soil moisture in East China, where the heatwave occurred. Meanwhile, there were significant soil moisture deficits and large-scale positive temperature anomalies, reaching roughly 6 ◦C, particularly in the middle and lower reaches of the Yangtze River basin. Atmospheric circulation anomalies are generally considered to be the main cause of heatwaves in China, e.g., the movement of the Northwest Pacific subtropical high [46,48,49]. However, land surface feedback on the atmosphere have been found to be an important factor for heatwaves over China, which may contribute 30–70% of the high temperature anomalies [22]. Our finding from Figure 5a–c is very likely to tell such a story that soil moisture deficit resulted in, at least partly, a significant heatwave through the coupled processes between land and atmosphere.

To better understand the temporal evolution of the heat wave, we show in Figure 5d–f the changes of the temperature anomalies ( *T* ) and the heat anomalies ( *H* − *Hp* ) with time. Figure 5d shows the positive heat anomalies before the heatwave occurred (12–23 July) in Southern China, which are related to the soil moisture deficit and lead to enhanced evaporative stress. Figure 5e illustrates the mega-heatwave from 24 July to 16 August in Southern China, when both the temperature anomalies (*T* ) and the heat anomalies ( *H* − *Hp* ) reached their maximum values with the largest spatial coverage. Figure 5f shows the temperature anomalies ( *T* ) and the heat anomalies ( *H* − *Hp* ) during 17–28 August, when the heatwave had almost vanished with largely reduced temperature and heat anomalies, as well as their spatial coverages. Further, in Figure 5g, we show the temporal variations of *H* − *Hp* and *T* from 1 June to 13 August, which are averaged within a small region in the epicenter of the heatwave (marked in Figure 5e). The right Y-axis indicates the metric π, and the left Y-axis indicates the anomalies of *H* − *Hp* and *T* . It clearly shows how soil moisture deficit contributed to the enhancement of heat and temperature anomalies; it is obvious that the land–atmosphere coupling was strongest with the largest π values during the mega heatwave (23 July–18 August).

**Figure 5.** The soil moisture–temperature coupling and related processes during the summer heatwave of Southeast China during 24 July–16 August 2013. (**a**) The coupling metric π; (**b**) the soil moisture anomalies and (**c**) the temperature anomalies, referring to the multi-year average of 1980–2013; (**d**) pre-heatwave (12–23 July); (**e**) the occurrence of the heatwave (24 July–16 August); (**f**) the days after the heatwave (17–28 August); (**g**) daily time series of *T* , *H* , *Hp* , and the coupling metric π.

3.2.2. Case 2: Heatwave of North China in Summer 2009

There was a heatwave from 20 June to 4 July in 2009 in North China. Although it was not as severe as the mega-heatwave in Southeast China in 2013, the continuous hot weather in the North China plain was rare since 1949 [50]. At the end of June, an unprecedented heatwave hit Hebei and Shandong provinces, setting new temperature records. By early July, the heatwave gradually dissipated in North China and the high temperature shifted to Southern China [51]. The cause analysis of this event mainly focused on the effect of El Niño and is typical of the high-pressure systems [49]. In addition, the soil moisture deficit of North China would increase the sensible heat flux and influence the atmospheric boundary layer temperature, conducive to strengthening the subtropical high and causing a heatwave [52]. Figure 6a–c depicts the coupling metric π, as well as the soil moisture and temperature anomalies referring to multi-year average of 1980–2013 during the heatwave (20 June–4 July) in North China in 2009. Where there is the strongest coupling, there is significant soil moisture deficit and maximum temperature anomalies.

Figure 6d–f is analogous to Figure 5d–f, dividing the study period into pre-heatwave (the left panel), the mega heatwave (the middle panel), and post-heatwave (the right panel). It appears that the heat anomalies (*H* − *Hp* ) and the temperature anomalies (*T* ) reached their maximum values during the heatwave, while such anomalies existed in neither the left nor the right panel. In Figure 6f, we show *H* − *Hp* and *T* for the period 1 June–31 August, averaged within a small region in the epicenter of the heatwave (marked in Figure 6e). Similarly, it is obvious that the land–atmosphere coupling is strongest with the largest π values during the mega heatwave 20 during June to 4 July, when there is the largest heat and temperature anomalies associated with a significant soil moisture deficit.

**Figure 6.** The soil moisture–temperature coupling and related processes during the summer heatwave of North China during 20 June–4 July 2009. (**a**) The coupling metric π; (**b**) the soil moisture anomalies and (**c**) the temperature anomalies referring to the multi-year average of 1980–2013; (**d**) pre-heatwave (12–19 June); (**e**) the occurrence of the heatwave (20 June–5 July); (**f**) the days after the heatwave (6–13 July); (**g**) daily time series of *T* , *H* , *Hp* , and the coupling metric π.

#### **4. Conclusions and Discussion**

This study attempted to utilize the GELAM and ERA-Interim datasets to study land–atmosphere coupling in China for the period 1980–2013. The key findings are as follows.

Hot spots of soil moisture–temperature coupling were found in North China and over the Tibetan Plateau, which indicate that the soil moisture–temperature coupling is strongest in the transitional climate zones. These results are in agreemen<sup>t</sup> with [9,14]. The seasonality of soil moisture–temperature coupling strength has marked regional variation, which suggests that soil moisture–temperature coupling strength is stronger in Northern China than in the southern, and coupling is relatively stronger in spring, followed by summer and autumn, and insignificant in winter. In spring, soil moisture–temperature coupling is stronger within dry areas, and these regimes are mainly water-limited regions, where evaporation depends on the supply of water.

Case studies involving the 2013 Southeast China heatwave and the 2009 North China heatwave were conducted to understand the role of soil moisture–temperature coupling and the related heating processes during heatwave events. It was found that enhanced heat and temperature anomalies associated with soil moisture deficit, when the soil moisture–temperature coupling intensifies, could result in enhanced evaporative stress and heat anomalies, which finally leads to enhanced soil moisture coupling with temperature [1,50]. However, there is much debate about the exact physical mechanisms of how heatwaves occur and evolve. Mirelles et al. [11] found that the prevailing persistent synoptic patterns led to warm air advection and clear skies, along with a high atmospheric demand, intensified soil desiccation (causing a strong surface sensible heat flux), causing the mega-heatwaves of 2003 and 2010 in Europe. Zhang et al. [53] supports that soil moisture in spring and early summer may be an important contributor to heatwaves in July via positive subtropical high anomalies. Several studies show that the formation mechanism of heatwaves is not only caused by a certain external force, but may be influenced by circulation systems, external forcing, or local e ffect, and soil moisture deficit may contribute directly and indirectly to all of these processes [10].

Land–atmosphere coupling involves water, energy, and chemical elements, which a ffect di fferent processes in the hydrological cycle and thus play a critical role in the climate system [8]. However, in-situ observations of soil moisture and land surface fluxes are scarce and uncertain at the appropriate scale, which has caused grea<sup>t</sup> di fficulties in the study of land–atmosphere coupling. The recent advances in satellite remote sensing have provided near-real-time datasets for us, and reanalysis data can also provide long-term databases for such studies. The long-term memory of soil moisture could help better understand the land–atmosphere interactions and may provide valuable information in weather forecasts to aid the managemen<sup>t</sup> of extremely warm climates.

Other dynamic relations between soil moisture and atmosphere coupling were ignored in this study, which may also a ffect the participation of sensible heat and latent heat. Furthermore, the ultimate causal relationship was not demonstrated in the diagnostics between soil moisture and temperature coupling. The physical mechanism and causal relationship between soil moisture and temperature still need to be further explored in future work.

**Author Contributions:** Q.Y. and G.W. conceived and designed the overall project, and wrote most of sections of the manuscript. C.Z. and D.F.T.H. performed methodology and data analysis. D.L., X.M. and M.Z. supplied suggestions and comments for the manuscript. All authors have read and agreed the published version of the manuscript.

**Funding:** This study was funded by the National Natural Science Foundation of China (41875094, 41850410492).

**Acknowledgments:** All authors are grateful to the anonymous reviewers and editors for their constructive comments on earlier versions of the manuscript.

**Conflicts of Interest:** The authors declare no conflict of interest.
