1. Introduction
The observed global average temperatures have been increasing since the midtwentieth century, and the rate of change has been amplified by anthropogenic forcing (Meehl et al. 2004; IPCC 2013; Q. Li et al. 2017). Impacts of climate change include changes not only in mean temperatures but also in the variability and extreme values of the diurnal temperature range (DTR) (e.g., Zhai and Pan 2003) and precipitation (e.g., Zhai et al. 2005; Guo et al. 2014, 2016; Li et al. 2016; Guo et al. 2017). It has long been known that the observed global warming is generally characterized by a greater increase in minimum temperatures than maximum temperature, which induces a decline in the DTR represented by the difference between the maximum temperature and minimum temperature during the past decades (e.g., Karl et al. 1993; Easterling et al. 1997; Zhai and Ren 1999; Durre and Wallace 2001; Vose et al. 2005; Lauritsen and Rogers 2012). The decline in the DTR has been identified as an index of global climate change (Braganza et al. 2004). In addition, the changes in the DTR may also affect the crop yield (Lobell 2007) and the net exchange of carbon dioxide in the ecosystems strongly depends upon diurnal temperature changes (e.g., Yi et al. 2010). It is therefore of great importance to elucidate the DTR changes under global warming. Many previous studies have sought to identify the controlling factors of DTR changes (e.g., Stenchikov and Robock 1995; Dai et al. 1999, 2006; Stone and Weaver 2002, 2003; L. Zhou et al. 2007; Liu et al. 2016). Cloud cover and precipitation are widely recognized to play an important role in the variations of DTR due to their significant influences on the surface energy and hydrological balance (e.g., Karl et al. 1993; Dai et al. 1999; Xia 2013). The clouds can cool Earth’s surface by reflecting downward solar shortwave radiation in daytime and emitting the downwelling longwave radiation back to surface in nighttime, thereby exerting a significant damping effect on DTR. The precipitation can influence the DTR not only by inducing the increase in the cloud cover, but also by strengthening the evaporative cooling attributed to more soil moisture in rainy days. Therefore, the negative correlation between DTR and precipitation variations largely results from the close association of precipitation with cloud cover and soil moisture (e.g., Dai et al. 1997; Zhou et al. 2009a). Up to 80% of the variance of DTR could be explained by cloud and precipitation records (Dai et al. 1999). In addition, other factors can contribute to the DTR variability, including changes in greenhouse gases (GHGs), aerosols, soil moisture, and surface properties (e.g., Hansen et al. 1995; Mitchell et al. 1995; Stenchikov and Robock 1995; Collatz et al. 2000; Stone and Weaver 2003; Zhou et al. 2004; L. Zhou et al. 2007; Dai et al. 2006; Huang et al. 2006; Liu et al. 2016; Xu et al. 2017).
Although some climate factors, such as cloud cover, precipitation, and anthropogenic aerosols, are regarded as the causes for the variations in DTR (e.g., Dai et al. 1999; Zhou et al. 2009a; Xia 2013; Liu et al. 2016; Z. Li et al. 2017), limited studies have been carried out to analyze the atmospheric circulation related to the changes in DTR (e.g., Wu 2010; Ionita et al. 2012). Since cloud cover and precipitation are both significantly influenced by atmospheric circulation, the changes in atmospheric circulation may play a nonnegligible role in the variations of DTR. Previous studies have demonstrated that the climate anomaly in East Asia (EA) shows some potential relationship to that in Australia (AUS) as a result of the coupling of atmospheric circulation between EA and AUS during boreal winter (e.g., Chen et al. 2000; Wang et al. 2009a; Zhang and Zhang 2010; Gong et al. 2015). Although there have been some studies on the DTR change in EA (e.g., Zhai and Ren 1999; Xia 2013; Shen et al. 2014) and AUS (e.g., Plummer et al. 1995) separately, the possible covariations of DTR between EA and AUS have not been investigated yet. Therefore, an interesting question is raised: Are there coherent variations of DTR between EA and AUS during boreal winter? If so, what roles do cloud cover and precipitation play in the variations of DTR in EA and AUS and what are the dynamic mechanisms associated with the changes of DTR in these regions?
In this study, the empirical orthogonal function (EOF) method is at first employed to identify the large-scale pattern of the DTR change in the EA–AUS region. The roles of atmospheric circulation and related dynamic mechanism responsible for the DTR changes are further investigated. The rest of the paper is organized as follows. Section 2 describes data and methods used in this study. The large-scale pattern of boreal winter DTR variations and its relationship to the cloud cover and precipitation in EA–AUS are presented in section 3. In section 4, we discuss the dynamic mechanisms of atmospheric circulation responsible for the variations of wintertime DTR in EA–AUS. A summary of the key findings and discussion are given in section 5.
2. Data and methods
a. Data
The observational monthly DTR, total cloud cover (CLT), and precipitation (Pre) used in this study are provided by the Climatic Research Unit (CRU) of the University of East Anglia [Time Series version 4.00 (TS4.00)], with a horizontal resolution of 0.5° × 0.5° covering the period from 1901 to 2015 (Harris et al. 2013; Harris and Jones 2017). To investigate the atmospheric processes associated with DTR variations, the present study uses the observational proxies of atmospheric circulation and air temperature data from the National Centers for Environmental Prediction (NCEP)–National Center for Atmospheric Research (NCAR) reanalysis (Kalnay et al. 1996) covering the period from January 1948 to the present with a horizontal resolution of 2.5° × 2.5°. To verify the results from the NCEP–NCAR reanalysis, a parallel analysis based on JRA-55 (Kobayashi et al. 2015; Harada et al. 2016) has to be performed, although this reanalysis covers a shorter time period (1958–2015) at a horizontal resolution of 1.25° × 1.25°. The CRU dataset has been bilinearly interpolated to a resolution of 2.5° × 2.5° (latitude × longitude) for convenience of analyzing the relationship to the atmospheric changes based on the reanalysis data. The sea surface temperature (SST) data from the monthly mean Extended Reconstructed Sea Surface Temperature, version 3 (ERSST.v3) dataset (Xue et al. 2003; Smith et al. 2008) were employed to investigate the possible linkage between SST and DTR variations; the data have a 2° × 2° horizontal resolution and cover the period from January 1854 to the present. The present analysis is based on the time period from 1948 to 2015, during which the above variables are available.
b. Methods
The EOF analysis has been a widely used statistical method in the climate research since it was first introduced to meteorology (e.g., Lorenz 1956; von Storch and Zwiers 1999). In particular, it is a useful tool to extract the dominant patterns with coherent temporal variations in a geophysical field (e.g., Hannachi et al. 2007; Monahan et al. 2009). In this study, the EOF method was applied to reveal the dominant patterns of DTR variations over EA–AUS. The Niño-3.4 index defined as SST anomaly averaged over the equatorial east-central Pacific region (5°S–5°N, 170°–120°W) was adopted to represent El Niño–Southern Oscillation (ENSO) (Rasmusson and Carpenter 1982). In this study, the boreal winter was considered and the winter means were constructed by averaging monthly mean data of December, January, and February (DJF). Here, our convention is that the winter of 1948 refers to the 1947/48 winter. Regression and correlation analyses were used to investigate the linkage of CLT and Pre as well as the associated atmospheric processes to the leading EOF patterns of DTR variations. The significance of a 95% confidence level (α = 0.05 level) of the results was evaluated with the two-tailed Student’s t test.
3. The leading modes of the winter DTR variations in EA–AUS
The EOF analysis is applied to winter DTR in the region of 40°S–50°N, 90°–155°E to reveal the characteristics of the DTR variations in the EA–AUS region. The EOF analysis was carried out by constructing an area-weighted covariance matrix. In an area-weighted covariance matrix, a cosθ weighting, where θ refers to latitude, was applied to the analyzed variable, taking into account of the convergence of the meridians with latitude. In the following discussion, the first two leading modes are analyzed, both of which account for about 45% of the total variances of DTR in EA–AUS. These leading modes are well separated from the rest of the modes, according to the criteria proposed by North et al. (1982).
a. EOF1 of the DTR: Trend mode
The spatial pattern of the first EOF mode (EOF1) is characterized by the same-sign loading over almost all the regions of EA–AUS and an opposite loading in part of southwestern Australia (Fig. 1a). The EOF1 explains 27.3% of the total variance of the DTR in the EA–AUS region. The temporal evolution of the corresponding principal component (PC1) exhibits a sharp drop during the whole time period of 1948–2015 (Fig. 1c). The power spectrum of PC1 shows that the PC1 is dominated by extreme low-frequency variations (Fig. 1e). The linear trend of PC1 reaches −0.43 standard deviation per decade. This pronounced trend of PC1 implies that the EOF1 actually reflects the trend of DTR in the EA–AUS region. To verify our hypothesis, the linear trend of wintertime DTR is displayed in Fig. 2. It can be clearly seen that the winter DTR shows a significantly decreasing trend in most areas of EA–AUS, as opposed to an increasing trend in part of southwestern Australia. The spatial distribution of the DTR trend is almost the same as the EOF1 pattern. This confirms that the trend is the dominant characteristic of DTR variability in EA–AUS, which agrees with previous studies that revealed a decreasing DTR trend over most of the globe (e.g., Karl et al. 1993; Easterling et al. 1997; Vose et al. 2005).
Also, the DTR has been widely recognized to be significantly associated with CLT and Pre because of their damping effects on the solar radiation (e.g., Dai et al. 1999; Zhou et al. 2009a; Xia 2010). An increase in cloud cover reduces the downward solar radiation, contributing to surface cooling during the daytime, and enhances downward longwave radiation, contributing to surface warming during the nighttime, which results in a net effect of decrease in DTR (Dai et al. 1999). To illustrate the possible relationship between EOF1 and CLT, we display in Fig. 3a the correlation map that is obtained by the temporal correlation coefficient between PC1 time series and the CLT time series at each grid point. It shows that the correlation between PC1 and CLT is significantly negative in most areas of the EA–AUS region, and small region with positive correlation is located in southwestern AUS and north of 45°N in EA. This spatial pattern of correlations is quite similar to that of EOF1 (Figs. 1a and 3a), suggesting that the CLT indeed has a good association with DTR in EA–AUS. Moreover, the correlation between PC1 and Pre in EA–AUS is somewhat similar to that of CLT, but with a smaller magnitude in EA (Fig. 3b). This suggests a weaker influence of Pre on DTR associated with EOF1 in EA.
We estimate the spatial distribution of linear trend for CLT and Pre to further investigate the relationship between long-term DTR and CLT and Pre variations. The spatial pattern of the linear trend of CLT is quite similar to that of the correlation map between CLT and PC1, as shown in Fig. 3a, with positive trends in most parts of EA–AUS, and negative trend in southwestern AUS and north of 45°N in EA. The trend of Pre is much weaker and insignificant in most areas of EA (Fig. 3d) compared with that in CLT. The significant negative correlations between CLT and DTR in most areas of EA–AUS suggest that the cloud cover is a crucial factor in the changes of DTR in EA–AUS. This confirms the previous findings that the damping effect of cloud cover on solar radiation plays an important role in the trend of DTR (Dai et al. 1997, 1999; Xia 2013; Shen et al. 2014).
It is noted that although the trend of CLT in most areas of eastern China is positive, the amplitude of CLT trend is relatively weak, which implies that, in addition to CLT, there may be other factors that lead to the decline in DTR. L. Zhou et al. (2007) reported that the DTR was reduced maximally in dry areas of the western Sahel during a period of intensive drought. They proposed that the decrease in DTR is largely attributed to the reduction in vegetation cover or the reduction in soil emissivity. Liu et al. (2016) pointed out that the anthropogenic aerosol plays an important role in the trend of DTR decrease in eastern China due to its negative radiation forcing in daytime. In consideration of the distinct climatic background and properties of surface, the factors responsible for the decline in DTR over different regions may be complicated. Therefore, although the CLT is a very important factor contributing to the decrease of DTR in EA–AUS, other factors, such as anthropogenic aerosol and land processes, may also influence the DTR to a certain extent.
b. EOF2 of the DTR: Interannual variability mode
As illustrated in Fig. 1b, the spatial pattern of EOF2 displays an opposite loading between EA and AUS. Overall, the EOF2 accounts for 17.6% of the total variance of DTR in the EA–AUS domain. The corresponding time series (PC2) shows pronounced fluctuation on interannual time scales with a spectrum peak around 2–3 yr (Figs. 1d,f). During the positive phases of PC2, the DTR has significant negative anomaly in most areas of EA and positive anomaly in the whole AUS region, and vice versa. Therefore, the EOF2 actually represents the out-of-phase variations of DTR over EA–AUS on interannual time scale; that is, during the positive phases of PC2, there are negative DTR anomalies in the EA and positive DTR anomalies in the AUS and vice versa.
4. Physical processes for the leading variations of wintertime DTR in EA–AUS
The above results indicate that the dominant patterns of the winter DTR variations in EA–AUS are closely tied to the changes in CLT. In this section, we investigate the atmospheric dynamic processes associated with the variation of the cloud and DTR based on the two leading modes. Regression analyses are employed to analyze the change of atmospheric circulation corresponding to the change of the specified indices (e.g., PCs and Niño-3.4). This change is regarded as the anomaly in the following section.
a. Role of nonuniform warming in the trends of DTR
The atmospheric circulation plays an important role in the formation and maintenance of CLT and Pre. For example, in the tropics, the organized deep convective clouds form in the ascending branches of the Hadley and Walker circulations, while the low clouds cover the ocean in anticyclonic areas (e.g., Bony et al. 2015). The changes of atmospheric circulation related to the EOF1 of DTR variations are first investigated. Because the EOF1 actually represents the trend of DTR in the EA–AUS region, we mainly analyze the trend of atmospheric circulation over 1948–2015. Figure 5a shows the trends of winter sea level pressure (SLP) and 1000-hPa winds in EA–AUS during 1948–2015.
The SLP shows a pronounced negative trend in most part of China and a positive trend in Japan and the South China Sea (SCS). The reduced contrast of sea–land pressure differences in EA induces the southerly wind responses along coastal China, weakening the East Asian winter monsoon (EAWM). This result is consistent with the previous findings that the EAWM has been weakening during the past several decades (e.g., W. Zhou et al. 2007; Wang et al. 2009b; Li and Yang 2017). The southerly anomalies reduce the transport of cold and dry air from higher latitudes over EA, which is beneficial to the increase of atmospheric moisture. Meanwhile, the changes in winds feature cyclonic anomalies over EA that are conducive to the ascending motion of atmosphere. These processes contribute to an increase in the cloud amount but a reduction in the DTR in EA (Figs. 2 and 3c).
Meanwhile, there are positive SLP trends mostly along the coastal AUS, whereas the trend of SLP is negative in central AUS (Fig. 5a). The resultant change in winds shows an anomalous cyclone response in the central AUS, accompanied by northerly anomalies in the central-eastern AUS and southerly anomalies in the southwestern AUS (Fig. 5a). The northerly wind anomalies bring more warm and wet air from low latitudes to eastern AUS. The anomalous cyclonic winds in central AUS induce anomalous ascending motion, thereby facilitating the increase in cloud and decrease in the DTR in eastern-central AUS. Meanwhile, weak southerly anomalies are present over southwestern AUS and they transport the relatively cold air from midlatitudes to southwestern AUS, leading to a small region of decrease in clouds (Fig. 3c), hence leading to increase in DTR in the southwestern AUS (Fig. 2).
Since the changes of atmospheric circulation are closely related to the changes of thermal condition of atmosphere itself, the trends of 1000-hPa air temperature (T1000) during 1948–2015 in the EA–AUS region are examined. Although the global average temperature is increasing over the past several decades due to increased GHGs, the trends of temperature on different regions are not uniform because of the different properties of surface (IPCC 2013). The trends of T1000 indeed show distinct differences in EA–AUS (Fig. 5b). Over EA, there is large warming north of 30°N in China, but small warming in the western North Pacific and even cooling in Japan and southern China. Consequently, the land–sea thermal gradient is weakened in EA. This reduces the land–sea pressure difference, thereby generating anomalous southerly wind responses along coastal China (Fig. 5a). Over AUS, there is cooling along the west coast of AUS and weak warming in central AUS (Fig. 5b). This thermal condition induces an anomalous cyclone response in central AUS, accompanied by northerly anomalies in central-eastern AUS and southerly anomalies in southwestern AUS (Fig. 5a). It is noted that the trend of wind shows an anomalous anticyclone along the east coast of AUS centered around 160°E, which also contributes to the northerly anomalies in eastern AUS. However, the center of the anomalous anticyclone does not match well with the cooling center of temperature trend located east of 160°E. It implies that although the local thermodynamic response of atmosphere to nonuniform surface warming appears to be an important reason for the changes of circulation, other processes may also contribute to the circulation changes in AUS.
The above results suggest that anomalous water vapor transportation and vertical motions associated with anomalous circulation may significantly influence the trends of cloud cover, thereby affecting the trends of DTR in EA–AUS. Meanwhile, the trends of atmospheric circulation are mainly attributed to the atmospheric thermodynamic responses to the nonuniform trend of temperature in EA–AUS, especially in EA. It is noted that the nonuniform warming in the EA–AUS may be contributed by the anthropogenic forcing. Although the radiation forcing associated with the increasing GHG emission is relatively uniform on globe, the response of regional temperature could be different as a result of the distinct climatic background and surface properties (IPCC 2013). Therefore, the decline in DTR in the EA–AUS region may also be influenced by the anthropogenic forcing. Previous studies found that the climate models can simulate a trend of decrease in the DTR, but with a much smaller magnitude than that found in observations (e.g., Karoly et al. 2003; Braganza et al. 2004; Zhou et al. 2009b, 2010). Therefore, the present results may serve to provide insight into plausible reasons for the weak trend of DTR in the climate models; that is, the climate models may not well reproduce the observational nonuniform warming.
b. Role of ENSO in the interannual variability of DTR
To investigate the factors governing the variations of EOF2, the atmospheric circulation and SST anomalies associated with EOF2 are examined. Figure 6a reveals the regressed 850-hPa wind and SST anomalies in tropical Indo-Pacific Ocean on the normalized time series of PC2. During the positive phases of PC2, an anomalous lower-tropospheric anticyclone is located in the western North Pacific, which is accompanied by anomalous southerly winds along the northwest of the anticyclone (Fig. 6a). The southerly anomalies over the EA not only weaken the EAWM, but also bring warm and wet air from low latitudes to EA and induce the increase in CLT and Pre, thereby reducing the DTR in EA (Figs. 4a,b). AUS is controlled by a large anomalous anticyclone over the east Indian Ocean to the entire AUS region. The anomalous descending motion associated with the anomalous anticyclone can significantly reduce the CLT and Pre in the western-central AUS (Fig. 4). Meanwhile, the southerly anomalies in the eastern part of the anomalous anticyclone bring the relatively cold and dry air to eastern AUS and reduce the CLT and Pre there. Therefore, both processes lead to the increase in DTR in AUS.
As shown in Fig. 6a, the SST anomaly pattern associated with PC2 manifests an evident warming in the tropical eastern-central Pacific and Indian Ocean during the positive phases of PC2. This pattern of SST anomalies is quite similar to that of ENSO that is a dominant interannual mode of climate variability. Figure 6b shows the regression map of SLP, 200-hPa velocity potential, and divergent wind anomalies relative to PC2. There are significant positive and negative SLP anomalies in the western and eastern Pacific, respectively. This pattern is also similar to the classic Southern Oscillation (e.g., Ji et al. 2015). In the upper troposphere, an anomalous divergence center is located in the eastern-central Pacific and convergence centers are mainly located in western Pacific and North Pacific (Fig. 6b). The locations of the dipole structure of convergence and divergence centers in the tropical Pacific well match the dipole structure of SLP anomaly pattern, which indicates that the Walker circulation is weakened in association with the significant warming in the tropical eastern-central Pacific, with a complementary contribution from the tropical Indian Ocean warming. Therefore, the anomalous descending motion associated with the weakening Walker circulation forms a pair of lower-tropospheric anomalous anticyclones over the tropical western Pacific (Fig. 6a). The center of the anomalous anticyclone on the north side of the equatorial western Pacific shifts eastward slightly compared to that in the Southern Hemisphere because of the strong positive local air–sea feedback in the western North Pacific (Wang et al. 2000). These results suggest that the out-of-phase interannual variations of the winter DTR between EA and AUS are mainly attributed to the opposite effects of the atmospheric circulation caused by the SST anomalies in tropical Pacific. It is noted that ENSO represents the natural variability in the climate system, and therefore this opposite-sign pattern of the DTR changes in the EA–AUS region may be attributed to the natural change and be independent of the anthropogenic forcings.
The above analysis indicates that the EOF2 may be largely influenced by ENSO. To further examine the influence of ENSO, the Niño-3.4 index defined as the SST anomaly averaged over the equatorial east-central Pacific region (5°S–5°N, 170°–120°W) was employed. Figure 7 displays the normalized time series of the Niño-3.4 index and PC2. There are consistent variations between PC2 and the Niño-3.4 index and the correlation coefficient between the two time series reaches 0.37 during 1948–2015, exceeding the 99% confidence level. This confirms that the interannual variability of DTR in EA–AUS represented by EOF2 is largely influenced by ENSO.
To delineate the influence of ENSO on the circulation anomalies over EA–AUS in detail, the regressed SST, 850-hPa wind, SLP, 200-hPa velocity potential, and corresponding divergent wind anomalies onto normalized Niño-3.4 index over tropical Indo-Pacific are presented in Fig. 8. It can be seen that the spatial pattern of ENSO-related SST and circulation anomalies are similar to those of PC2 (Fig. 6). During the positive phases of the Niño-3.4 index, the pattern of SST anomalies associated with Niño-3.4 index shows evident warming in the tropical eastern-central Pacific (Fig. 8a). The resultant weakened Walker circulation (Fig. 8b) triggers a pair of the lower-level anticyclone north and south of the equatorial western Pacific (Fig. 8a). Southerly anomalies in the west side of the anomalous anticyclone near the Philippine Sea bring warm and wet air from low latitudes to EA and induce the increase in CLT and Pre, thereby reducing the DTR in EA (Fig. 8a). The AUS region is also controlled by a large anomalous anticyclone, extending from the central Indian Ocean to the South Pacific. The anomalous descending motion associated with the anomalous anticyclone can significantly reduce the CLT and Pre in AUS (Fig. 8a).
To further elucidate the influences of ENSO on the interannual variations of DTR, CLT, and Pre, the correlation maps between the Niño-3.4 index and DTR, CLT, and Pre in EA–AUS are presented in Fig. 9. It shows that the correlation between the Niño-3.4 index and DTR is basically opposite over EA and AUS (Fig. 9a). Meanwhile, the correlations between the Niño-3.4 index and CLT and Pre are broadly opposite in EA and AUS. The results indicate that ENSO induces opposite cloud and precipitation changes over EA and AUS, which in turn contribute to opposite variations of DTR between EA and AUS on the interannual time scale. It is noted that correlations are weak in the Maritime Continent (Fig. 9). This phenomenon may be attributed to the sparse observational stations and relatively low data quality in this region.
5. Summary and discussion
This study attempts to unveil the large-scale patterns of the winter DTR changes in EA–AUS and to investigate the physical processes responsible for these variations of DTR based on the observational and reanalysis data over the past 68 years (1948–2015). The EOF analysis reveals two leading modes with distinct winter DTR variations in EA–AUS. These two modes have different time scales and dynamic mechanisms.
The first mode is characterized by a same-sign variation in most areas of EA–AUS. The corresponding PC1 shows a decreasing trend during 1948–2015, which suggests that the EOF1 represents the trends of DTR in EA–AUS. The trends of DTR are mainly related to the changes of CLT in EA–AUS. The influence of Pre is relatively weak in EA–AUS, especially in EA. The changes of CLT in EA–AUS are closely tied to the changes of atmospheric circulation. The pronounced southerly trends over EA induce the northward transport of moisture from the low latitudes to EA. Meanwhile, the cyclonic wind anomalies in EA are also conducive to the ascending motion, thereby increasing the CLT and reducing the DTR in most areas of EA during the past decades. In addition, the change of wind shows an anomalous cyclone in central AUS accompanied by the northerly wind anomalies in eastern AUS and southerly wind anomalies in southwestern AUS. The anomalous cyclone induces a significant increase in CLT with the ascending motion, thus decreasing the DTR in central AUS during the past decades. The northerly wind anomalies bring the warm wet air from low latitudes to eastern AUS, leading to the increase in clouds and decrease in DTR in eastern AUS. The weak southerly anomalies transport the relative cold air from midlatitudes to southwestern AUS, inducing a decrease in clouds and hence an increase in DTR in the southwestern AUS during the past decades.
Furthermore, the changes of atmospheric circulation are largely due to the thermodynamic responses of atmosphere to the nonuniform warming of temperature in EA–AUS. The strongest warming in most areas of China and relatively weak warming in western North Pacific and cooling in Japan reduce the land–sea thermal contrast, thereby generating a southerly response in EA. The warming in AUS is also not uniform, with cooling in most of the coast of AUS but weak warming in central AUS. This thermal condition may trigger an anomalous cyclone response in central AUS, accompanied by the northerlies and southerlies anomalies in eastern-central and southwestern AUS, thereby increasing CLT in eastern-central AUS and decreasing CLT in southwestern AUS. Since the response of regional temperature to the radiative forcing associated with the increasing GHGs emission could be different because of the distinct climatic background and surface properties (IPCC 2013), the decline in DTR caused by the nonuniform warming in the EA–AUS region may result from anthropogenic forcing.
The second leading mode presents an opposite variation in DTR between EA and AUS, and shows fluctuations on interannual time scales. This mode is closely tied to CLT and Pre change, and the former plays a more important role in the variations of DTR. The second mode of DTR variations is found to be largely forced by the ENSO variability. Meanwhile, the remarkably opposite variations of DTR in EA and AUS are primarily due to the opposite effect of ENSO on the CLT and Pre in EA and AUS. During the positive phases of ENSO, the weakened Walker circulation associated with the warming in tropical eastern-central Pacific triggers a pair of anomalous low-level anticyclone over western Pacific. The CLT decreases in AUS owing to the descending motion related to the anomalous anticyclone which covers the entire AUS, thereby leading to increase the DTR in AUS. In contrast, the anomalous southerlies along the west of the anomalous anticyclone around the Philippine Sea transport warm and wet air to EA and lead to the increase in CLT, thus reducing the DTR in EA. Since ENSO is the natural variability in the climate system, this opposite-sign pattern of the DTR changes in the EA–AUS region may be largely attributed to the natural change.
Meanwhile, it is also noted that although the EOF2 is largely influenced by the ENSO, there are some years when the PC2 and Niño-3.4 index do not match well and are even opposite (Fig. 7). This indicates that some other factors may also contribute to the interannual DTR variations in the EA–AUS region. For instance, the Indian Ocean basinwide SST anomaly, generally regarded as a response to ENSO (e.g., Xie et al. 2009), can in turn affect the ENSO, including its amplitude, evolution, and climate impacts (e.g., Wu and Kirtman 2004; Kug and Kang 2006; Terray et al. 2016). In particular, the Indian Ocean warming (cooling) at the peak of El Niño (La Niña) in boreal winter was found to be able to modulate the effect of El Niño on air temperature over the African, Asian, and Australian landmasses via cloud forcing (Herold and Santoso 2017). Previous modeling studies (e.g., Terray et al. 2016) showed that without feedback from the Indian Ocean, the period of ENSO would be longer compared to the irregular ENSO period ranging from 2 to 7 yr. The short period (2–3 yr) ENSO signals as seen in the power spectrum of EOF2 (Fig. 1f), therefore, potentially highlight the role of the Indian Ocean. The detailed dynamic mechanisms regarding the influence of the Indian Ocean on DTR change need to be further investigated in future studies.
We perform a parallel analysis using JRA-55 for the period 1958–2015. The dominant modes of DTR variations in the EA–AUS region in the period of 1958–2015 show almost the same patterns as those in the period of 1948–2015. The trends of winter SLP, 1000-hPa winds, and 1000-hPa air temperature during 1958–2015 based on JRA-55 display major features similar to those during 1948–2015 based on the NCEP–NCAR reanalysis (Fig. S1 in the supplemental material). Corresponding to the EOF2, the El Niño–like SST pattern, a pair of lower-level anticyclones north and south of the equatorial western Pacific, the weakened Walker circulation represented by the dipole structure of SLP anomaly, and the upper-level convergence and divergence anomaly centers in the tropical Pacific are well reproduced in the JRA-55 dataset (Fig. S2 in the supplemental material).
The present study focuses on the DTR changes in EA–AUS region during boreal winter when the ENSO signal is the strongest. A question arises as to the status of DTR changes in the other seasons and the underlying plausible factors. Our preliminary analysis indicates that the leading patterns of DTR changes in the other seasons are not the same as those in boreal winter (figure not shown). This discrepancy may be related to the climatic background and the surface properties. For example, the SST anomaly patterns are quite different in different seasons. During boreal summer, the ENSO SST patterns have decayed and the SST anomalies are neutral in the eastern-central Pacific, whereas the Indian Ocean has the strongest warming during this season and significantly impacts the climate of surrounding regions (e.g., Yang et al. 2007; Xie et al. 2009). Hence, the change of atmospheric circulation forced by the SST anomalies is also significantly different compared with that in the boreal winter. More systematic investigation is needed to unravel the main features of DTR changes in other seasons and the factors and associated processes.
Acknowledgments
This work was supported by the Ministry of Science and Technology of China under Grants 2017YFC1501401 (PI: Jianping Guo) and 2017YFA0603501 (PI: Hua Zhang, co-PI: Jian Li), the National Natural Science Foundation of China under Grants 41771399, 41605060, 41661144016, and 41530425, and the Chinese Academy of Meteorological Sciences under Grant 2017Z005 (PI: Jianping Guo). The authors extend sincere gratitude to the Climatic Research Unit of the University of East Anglia for providing diurnal temperature range, total cloud cover, and precipitation data publicly accessible. Last but not least, the authors are grateful to the editor and three anonymous reviewers for their insightful comments that helped improve the manuscript.
REFERENCES
Bony, S., and Coauthors, 2015: Clouds, circulation and climate sensitivity. Nat. Geosci., 8, 261–268, https://doi.org/10.1038/ngeo2398.
Braganza, K., D. J. Karoly, and J. M. Arblaster, 2004: Diurnal temperature range as an index of global climate change during the twentieth century. Geophys. Res. Lett., 31, L13217, https://doi.org/10.1029/2004GL019998.
Chen, W., H. F. Graf, and R. Huang, 2000: The interannual variability of East Asian winter monsoon and its relation to the summer monsoon. Adv. Atmos. Sci., 17, 48–60, https://doi.org/10.1007/s00376-000-0042-5.
Collatz, G. J., L. Bounoua, S. O. Los, D. A. Randall, I. Y. Fung, and P. J. Sellers, 2000: A mechanism for the influence of vegetation on the response of the diurnal temperature range to changing climate. Geophys. Res. Lett., 27, 3381–3384, https://doi.org/10.1029/1999GL010947.
Dai, A., A. Del Genio, and I. Y. Fung, 1997: Clouds, precipitation, and temperature range. Nature, 386, 665–666, https://doi.org/10.1038/386665b0.
Dai, A., K. E. Trenberth, and T. R. Karl, 1999: Effects of clouds, soil moisture, precipitation, and water vapor on diurnal temperature range. J. Climate, 12, 2451–2473, https://doi.org/10.1175/1520-0442(1999)012<2451:EOCSMP>2.0.CO;2.
Dai, A., T. R. Karl, B. Sun, and K. E. Trenberth, 2006: Recent trends in cloudiness over the United States: A tale of monitoring inadequacies. Bull. Amer. Meteor. Soc., 87, 597–606, https://doi.org/10.1175/BAMS-87-5-597.
Durre, I., and J. M. Wallace, 2001: Factors influencing the cold-season diurnal temperature range in the United States. J. Climate, 14, 3263–3278, https://doi.org/10.1175/1520-0442(2001)014<3263:FITCSD>2.0.CO;2.
Easterling, D. R., and Coauthors, 1997: Maximum and minimum temperature trends for the globe. Science, 277, 364–367, https://doi.org/10.1126/science.277.5324.364.
Gong, H., L. Wang, W. Chen, D. Nath, G. Huang, and W. Tao, 2015: Diverse influences of ENSO on the East Asian–western Pacific winter climate tied to different ENSO properties in CMIP5 models. J. Climate, 28, 2187–2202, https://doi.org/10.1175/JCLI-D-14-00405.1.
Gong, H., L. Wang, W. Chen, X. Chen, and D. Nath, 2017: Biases of the wintertime Arctic Oscillation in CMIP5 models. Environ. Res. Lett., 12, 014001, https://doi.org/10.1088/1748-9326/12/1/014001.
Guo, J., and Coauthors, 2014: Diurnal variation and the influential factors of precipitation from surface and satellite measurements in Tibet. Int. J. Climatol., 34, 2940–2956, https://doi.org/10.1002/joc.3886.
Guo, J., and Coauthors, 2016: Delaying precipitation and lightning by air pollution over the Pearl River Delta. Part I: Observational analyses. J. Geophys. Res. Atmos., 121, 6472–6488, https://doi.org/10.1002/2015JD023257.
Guo, J., and Coauthors, 2017: Declining frequency of summertime local-scale precipitation over eastern China from 1970–2010 and its potential link to aerosols. Geophys. Res. Lett., 44, 5700–5708, https://doi.org/10.1002/2017GL073533.
Hannachi, A., I. Joliffe, and D. Stephenson, 2007: Empirical orthogonal functions and related techniques in atmospheric science: A review. Int. J. Climatol., 27, 1119–1152, https://doi.org/10.1002/joc.1499.
Hansen, J., M. Sato, and R. Ruedy, 1995: Long-term changes of the diurnal temperature cycle: Implications about mechanisms of global climate change. Atmos. Res., 37, 175–209, https://doi.org/10.1016/0169-8095(94)00077-Q.
Harada, Y., and Coauthors, 2016: The JRA-55 Reanalysis: Representation of atmospheric circulation and climate variability. J. Meteor. Soc. Japan, 94, 269–302, https://doi.org/10.2151/jmsj.2016-015.
Harris, I., and P. D. Jones, 2017: CRU TS4.00: Climatic Research Unit (CRU) Time-Series (TS) version 4.00 of high resolution gridded data of month-by-month variation in climate (Jan. 1901–Dec. 2015). Centre for Environmental Data Analysis, University of East Anglia Climatic Research Unit, accessed 25 August 2017, https://doi.org/10.5285/edf8febfdaad48abb2cbaf7d7e846a86.
Harris, I., P. D. Jones, T. J. Osborn, and D. H. Lister, 2013: Updated high-resolution grids of monthly climatic observations—The CRU TS3.10 dataset. Int. J. Climatol., 34, 623–642, https://doi.org/10.1002/joc.3711.
Herold, N., and A. Santoso, 2017: Indian Ocean warming during peak El Niño cools surrounding land masses. Climate Dyn., https://doi.org/10.1007/s00382-017-4001-6, in press.
Huang, Y., R. E. Dickinson, and W. L. Chameides, 2006: Impact of aerosol indirect effect on surface temperature over East Asia. Proc. Natl. Acad. Sci. USA, 103, 4371–4376, https://doi.org/10.1073/pnas.0504428103.
IPCC, 2013: Climate Change 2013: The Physical Science Basis. Cambridge University Press, 1535 pp., https://doi.org/10.1017/CBO9781107415324.
Ionita, M., G. Lohmann, N. Rimbu, and P. Scholz, 2012: Dominant modes of diurnal temperature range variability over Europe and their relationships with large-scale atmospheric circulation and sea surface temperature anomaly patterns. J. Geophys. Res., 117, D15111, https://doi.org/10.1029/2011JD016669.
Ji, X., J. D. Neelin, and C. R. Mechoso, 2015: El Niño–Southern Oscillation sea level pressure anomalies in the western Pacific: Why are they there? J. Climate, 28, 8860–8872, https://doi.org/10.1175/JCLI-D-14-00716.1.
Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437–471, https://doi.org/10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.
Karl, T. R., and Coauthors, 1993: A new perspective on recent global warming: Asymmetric trends of daily maximum and minimum temperature. Bull. Amer. Meteor. Soc., 74, 1007–1023, https://doi.org/10.1175/1520-0477(1993)074<1007:ANPORG>2.0.CO;2.
Karoly, D. J., K. Braganza, P. A. Stott, J. M. Arblaster, G. A. Meehl, A. J. Broccoli, and K. W. Dixon, 2003: Detection of a human influence on North American climate. Science, 302, 1200–1203, https://doi.org/10.1126/science.1089159.
Kobayashi, S., and Coauthors, 2015: The JRA-55 Reanalysis: General specifications and basic characteristics. J. Meteor. Soc. Japan, 93, 5–48, https://doi.org/10.2151/jmsj.2015-001.
Kug, J. S., and I.-S. Kang, 2006: Interactive feedback between ENSO and the Indian Ocean. J. Climate, 19, 1784–1801, https://doi.org/10.1175/JCLI3660.1.
Lauritsen, R. G., and J. C. Rogers, 2012: U.S. diurnal temperature range variability and regional causal mechanism, 1901–2002. J. Climate, 25, 7216–7231, https://doi.org/10.1175/JCLI-D-11-00429.1.
Li, Q., L. Zhang, W. Xu, T. Zhou, J. Wang, P. Zhai, and P. Jones, 2017: Comparisons of time series of annual mean surface air temperature for China since the 1900s: Observations, model simulations, and extended reanalysis. Bull. Amer. Meteor. Soc., 98, 699–711, https://doi.org/10.1175/BAMS-D-16-0092.1.
Li, Y., and S. Yang, 2017: Feedback attributions to the dominant modes of East Asian winter monsoon variations. J. Climate, 30, 905–920, https://doi.org/10.1175/JCLI-D-16-0275.1.
Li, Z., and Coauthors, 2016: Aerosol and monsoon climate interactions over Asia. Rev. Geophys., 54, 866–929, https://doi.org/10.1002/2015RG000500.
Li, Z., and Coauthors, 2017: Aerosol and boundary-layer interactions and impact on air quality. Natl. Sci. Rev., 4, 810–833, https://doi.org/10.1093/nsr/nwx117.
Liu, L., Z. Li, X. Yang, H. Gong, C. Li, and A. Xiong, 2016: The long-term trend in the diurnal temperature range over Asia and its natural and anthropogenic causes. J. Geophys. Res. Atmos., 121, 3519–3533, https://doi.org/10.1002/2015JD024549.
Lobell, D. B., 2007: Changes in diurnal temperature range and national cereal yields. Agric. For. Meteor., 145, 229–238, https://doi.org/10.1016/j.agrformet.2007.05.002.
Lorenz, E. N., 1956: Empirical orthogonal functions and statistical weather prediction. Statistical Forecasting Project Rep. 1, MIT Department of Meteorology, 49 pp.
Meehl, G. A., W. M. Washington, C. M. Ammann, J. M. Arblaster, T. M. L. Wigley, and C. Tebaldi, 2004: Combinations of natural and anthropogenic forcings in twentieth-century climate. J. Climate, 17, 3721–3727, https://doi.org/10.1175/1520-0442(2004)017<3721:CONAAF>2.0.CO;2.
Mitchell, J., R. Davis, W. Ingram, and C. Senior, 1995: On surface temperature, greenhouse gases, and aerosols: Models and observations. J. Climate, 8, 2364–2386, https://doi.org/10.1175/1520-0442(1995)008<2364:OSTGGA>2.0.CO;2.
Monahan, A. H., J. C. Fyfe, M. H. P. Ambaum, D. B. Stephenson, and G. R. North, 2009: Empirical orthogonal functions: The medium is the message. J. Climate, 22, 6501–6514, https://doi.org/10.1175/2009JCLI3062.1.
North, G. R., T. L. Bell, R. F. Cahalan, and F. J. Moeng, 1982: Sampling errors in the estimation of empirical orthogonal functions. Mon. Wea. Rev., 110, 699–706, https://doi.org/10.1175/1520-0493(1982)110<0699:SEITEO>2.0.CO;2.
Plummer, N., Z. Lin, and S. Torok, 1995: Trends in diurnal temperature range over Australia since 1951. Atmos. Res., 37, 79–86, https://doi.org/10.1016/0169-8095(94)00070-T.
Rasmusson, E. M., and T. H. Carpenter, 1982: Variations in tropical sea surface temperature and surface wind fields associated with the Southern Oscillation/El Niño. Mon. Wea. Rev., 110, 354–384, https://doi.org/10.1175/1520-0493(1982)110<0354:VITSST>2.0.CO;2.
Shen, X., B. Liu, G. Li, Z. Wu, Y. Jin, P. Yu, and D. Zhou, 2014: Spatiotemporal change of diurnal temperature range and its relationship with sunshine duration and precipitation in China. J. Geophys. Res. Atmos., 119, 13 163–13 179, https://doi.org/10.1002/2014JD022326.
Smith, T. M., R. W. Reynolds, T. C. Peterson, and J. Lawrimore, 2008: Improvements to NOAA’s historical merged land–ocean surface temperature analysis (1880–2006). J. Climate, 21, 2283–2296, https://doi.org/10.1175/2007JCLI2100.1.
Stenchikov, G. L., and A. Robock, 1995: Diurnal asymmetry of climatic response to increased CO2 and aerosols: Forcings and feedbacks. J. Geophys. Res., 100, 26 211–26 227, https://doi.org/10.1029/95JD02166.
Stone, D. A., and A. J. Weaver, 2002: Daily maximum and minimum temperature trends in a climate model. Geophys. Res. Lett., 29, 1356, https://doi.org/10.1029/2001GL014556.
Stone, D. A., and A. J. Weaver, 2003: Factors contributing to diurnal temperature range trends in twentieth and twenty-first century simulations of the CCCma coupled model. Climate Dyn., 20, 435–445, https://doi.org/10.1007/s00382-002-0288-y.
Terray, P., S. Masson, C. Prodhomme, M. Roxy, and K. Sooraj, 2016: Impacts of Indian and Atlantic oceans on ENSO in a comprehensive modeling framework. Climate Dyn., 46, 2507–2533, https://doi.org/10.1007/s00382-015-2715-x.
von Storch, H., and F. W. Zwiers, 1999: Statistical Analysis in Climate Research. Cambridge University Press, 494 pp.
Vose, R. S., D. R. Easterling, and B. Gleason, 2005: Maximum and minimum temperature trends for the globe: An update through 2004. Geophys. Res. Lett., 32, L23822, https://doi.org/10.1029/2005GL024379.
Wang, B., R. Wu, and X. Fu, 2000: Pacific–East Asian teleconnection: How does ENSO affect East Asian climate? J. Climate, 13, 1517–1536, https://doi.org/10.1175/1520-0442(2000)013<1517:PEATHD>2.0.CO;2.
Wang, L., W. Chen, W. Zhou, and R. Huang, 2009a: Interannual variations of East Asian trough axis at 500 hPa and its association with the East Asian winter monsoon pathway. J. Climate, 22, 600–614, https://doi.org/10.1175/2008JCLI2295.1.
Wang, L., R. Huang, L. Gu, W. Chen, and L. Kang, 2009b: Interdecadal variations of the East Asian winter monsoon and their association with quasi-stationary planetary wave activity. J. Climate, 22, 4860–4872, https://doi.org/10.1175/2009JCLI2973.1.
Wu, Q., 2010: Associations of diurnal temperature range change with the leading climate variability modes during the Northern Hemisphere wintertime and their implication on the detection of regional climate trends. J. Geophys. Res., 115, D19101, https://doi.org/10.1029/2010JD014026.
Wu, R., and B. P. Kirtman, 2004: Understanding the impacts of the Indian Ocean on ENSO variability in a coupled GCM. J. Climate, 17, 4019–4031, https://doi.org/10.1175/1520-0442(2004)017<4019:UTIOTI>2.0.CO;2.
Xia, X., 2010: Spatiotemporal changes in sunshine duration and cloud amount as well as their relationship in China during 1954–2005. J. Geophys. Res., 115, D00K06, https://doi.org/10.1029/2009JD012879.
Xia, X., 2013: Variability and trend of diurnal temperature range in China and their relationship to total cloud cover and sunshine duration. Ann. Geophys., 31, 795–804, https://doi.org/10.5194/angeo-31-795-2013.
Xie, S.-P., K. Hu, J. Hafner, H. Tokinaga, Y. Du, G. Huang, and T. Sampe, 2009: Indian Ocean capacitor effect on Indo–western Pacific climate during the summer following El Niño. J. Climate, 22, 730–747, https://doi.org/10.1175/2008JCLI2544.1.
Xu, H., and Coauthors, 2017: Warming effect of dust aerosols modulated by overlapping clouds below. Atmos. Environ., 166, 393–402, https://doi.org/10.1016/j.atmosenv.2017.07.036.
Xue, Y., T. M. Smith, and R. W. Reynolds, 2003: Interdecadal changes of 30-yr SST normals during 1871–2000. J. Climate, 16, 1601–1612, https://doi.org/10.1175/1520-0442-16.10.1601.
Yang, J., Q. Liu, S.-P. Xie, Z. Liu, and L. Wu, 2007: Impact of the Indian Ocean SST basin mode on the Asian summer monsoon. Geophys. Res. Lett., 34, L02708, https://doi.org/10.1029/2006GL028571.
Yi, C., and Coauthors, 2010: Climate control of terrestrial carbon exchange across biomes and continents. Environ. Res. Lett., 5, 034007, https://doi.org/10.1088/1748-9326/5/3/034007.
Zhai, P., and F. Ren, 1999: Changes of China’s maximum and minimum temperatures in 1951–1990. Acta Meteor. Sin., 13, 278–290.
Zhai, P., and X. Pan, 2003: Trends in temperature extremes during 1951–1999 in China. Geophys. Res. Lett., 30, 1913, https://doi.org/10.1029/2003GL018004.
Zhai, P., X. Zhang, H. Wan, and X. Pan, 2005: Trends in total precipitation and frequency of daily precipitation extremes over China. J. Climate, 18, 1096–1108, https://doi.org/10.1175/JCLI-3318.1.
Zhang, C., and H. Zhang, 2010: Potential impacts of East Asian winter monsoon on climate variability and predictability in the Australian summer monsoon region. Theor. Appl. Climatol., 101, 161–177, https://doi.org/10.1007/s00704-009-0246-2.
Zhou, L., and Coauthors, 2004: Evidence for a significant urbanization effect on climate in China. Proc. Natl. Acad. Sci. USA, 101, 9540–9544, https://doi.org/10.1073/pnas.0400357101.
Zhou, L., R. E. Dickinson, Y. Tian, R. S. Vose, and Y. Dai, 2007: Impact of vegetation removal and soil aridation on diurnal temperature range in a semiarid region—Application to the Sahel. Proc. Natl. Acad. Sci. USA, 104, 17 937–17 942, https://doi.org/10.1073/pnas.0700290104.
Zhou, L., A. Dai, Y. Dai, R. Vose, C.-Z. Zou, Y. Tian, and H. Chen, 2009a: Spatial dependence of diurnal temperature range trends on precipitation from 1950 to 2004. Climate Dyn., 32, 429–440, https://doi.org/10.1007/s00382-008-0387-5.
Zhou, L., R. E. Dickinson, P. Dirmeyer, A. Dai, and S.-K. Min, 2009b: Spatiotemporal patterns of changes in maximum and minimum temperatures in multi-model simulations. Geophys. Res. Lett., 36, L02702, https://doi.org/10.1029/2008GL036141.
Zhou, L., R. E. Dickinson, A. Dai, and P. Dirmeyer, 2010: Detection and attribution of anthropogenic forcing to diurnal temperature range changes from 1950 to 1999: Comparing multi-model simulations with observations. Climate Dyn., 35, 1289–1307, https://doi.org/10.1007/s00382-009-0644-2.
Zhou, W., C. Li, and X. Wang, 2007: Possible connection between Pacific oceanic interdecadal pathway and East Asian winter monsoon. Geophys. Res. Lett., 34, L01701, https://doi.org/10.1029/2006GL027809.