Abstract

Studying the variation in the boundary position of the Asian summer monsoon in the Tibetan Plateau (TP) region and its potential drivers is important for understanding the climate in this region. Three sets of mean monthly precipitation data from 1980 to 2019 were sourced from the Global Precipitation Climatology Centre, the Climate Research Unit, and China Meteorological Information Service Centre. Several indicators that represent the Asian summer monsoon boundary (ASMB) were selected to compare their applicability to the TP region and elucidate the changes in the location of the ASMB in the TP over the last four decades. The results showed that the ASMB in the TP region extends in a southwest-northeast direction, with a clear north-south variation. It reaches as far north as the Kunlun Mountains and as far south as the Himalayas. The largest amplitude in spatial fluctuation occurs in the middle of the TP, and the smallest amplitude occurs at both ends of the region. A “small-large-small” fluctuation pattern was observed from west to east. The water vapor mainly originates from the South Asian region. The South Asian summer monsoon can move the ASMB position northward, whereas the westerly wind moves the ASMB position southward. Variation in the ASMB in the TP region is closely associated with the South Asian monsoon and westerly wind.

1. Introduction

Monsoons have long been a popular topic of research in the meteorological community. The modern concept of monsoons was introduced as early as the late 17th century, by Halley [1], who argued that they are caused by differences in the heating of the ocean and land by the sun, which, consequently, lead to differences in atmospheric pressure. Later, Hadley [2] extended the monsoon model, whereby the rotation effect of the earth causes the airflow to be affected by the Coriolis force, deflecting it to the right to form a southwesterly wind. The concept of “global monsoon” (GM) has been introduced and improved through observations at different levels and in-depth studies in recent 20 years. “GM” studies have evolved from evaluating regional climate phenomena to investigating global scale effects [3], opening a new phase of research on monsoon precipitation and regions exposed to monsoons [47].

The Asian summer monsoon (ASM) is an important part of the GM system. It is the leading climate phenomenon and the main cause of precipitation in Asia [810]. Changes to the intensity of the ASM affect this region, including how rainbands are distributed, which may lead to meteorological disasters, such as flooding, under precipitation anomalies [1113]. On the TP, the precipitation is less than that in the core areas of Asia, such as the East and South Asian area; therefore, the precipitation in the region varies considerably, which, in turn, impacts the vegetation and ecology along the monsoon boundary, and it is also more sensitive to the climate at the edge of the ASM [1418]. Thus, there is a need to assess how ASMB in the TP changed in the past four decades to better understand the climatic characteristics of TP and GM.

Extensive efforts have been undertaken to delineate the boundary of the monsoon. Previously, the “boundary” was called the northern edge zone of the monsoon, with thermal and dynamic factors being used to monitor how the position of the monsoon boundary changes. For example, Wang et al. [19] defined the monsoon boundary as areas with southwesterly winds >2.0 ms−1, pseudo-equivalent potential temperature (θse) south of 30°N greater than 340 K, and θse north of 30°N greater than 335 K. Wu et al. [20] defined the monsoon boundary as five-day mean southwesterly winds at 850 hPa and southerly winds >3.0 ms−1, and a five-day mean θse at 850 hPa greater than 335 K. However, the monsoon boundary varies greatly depending on the indicators used to define it. Because the summer monsoon carries a large amount of water vapor, the monsoon boundary could also be defined based on precipitation indicators including annual precipitation of 300–500 mm [21], five-day average precipitation [22], or 20 mm isoline of precipitation variance during summer (June–August) [23]. Additionally, the boundary can be directly delineated from the standardized precipitation indicator, by establishing areas with a difference of 2.0 mm day−1 between summer and winter precipitation and more than 55% of summer precipitation as monsoon zones [24], and using the 2.0 mm day−1 isoline of summer (May–September) precipitation to determine a northern monsoon boundary indicator [25]. However, defining the monsoon boundary index using kinetic or thermodynamic indicators could have limitations, due to the large topographic effect of the TP [26], with precipitation data possibly being more informative. In this study, we compared the monsoon boundary index by using monthly mean precipitation data and explored its applicability in the TP. We expect that our results will provide new insights on the spatial and temporal fluctuations of the monsoon boundary, and possible factors affecting its location.

2. Data and Methods

2.1. Study Area and Data

The TP (26–40°N and 73–104°E) is the highest plateau in the world, with an average elevation of >4500 m [27]. The TP is located to the south of the southern edge of the Himalayas, north of the Kunlun and Qilian Mountains, west to the Pamir Plateau, and east and northeast to the western section of the Qinling Mountains; it has a total area of approximately 2.5 million km2 [28] (see Figure 1).

Three different monthly precipitation datasets covering the period from 1980 to 2019 were used in this study. The first dataset was obtained from the Global Precipitation Climatology Centre (GPCC) [29], with a horizontal resolution of 0.5° × 0.5°longitude. The second dataset was obtained from the China Meteorological Data Service Centre, providing monthly precipitation data of China’s surface precipitation (CN05) [30], with a horizontal resolution of 0.25° × 0.25°. The third dataset was obtained from the Climate Research Unit (CRU) [31, 32], which provided monthly precipitation data at a horizontal resolution of 0.5° × 0.5°. Monthly wind data were obtained from the monthly reanalyzed data of the National Center for Environmental Prediction/Atmospheric Research (NCEP/NCAR) [33], with a horizontal resolution of 2.5° × 2.5°. The surface wind field was based on reanalyzed monthly values of era5-land monthly averaged data from 1980 to 2019, with a horizontal resolution of 0.1° × 0.1° [34]. Water vapor flux was calculated from era5 monthly averaged pressure data from 1980 to 2019, with a horizontal resolution of 0.25° × 0.25° and a vertical resolution of 37 layers [35]. Land cover type data were obtained from MODerate Resolution Imaging Spectroradiometer (MODIS) data from October 2000 to October 2001, with a horizontal resolution of 0.25° × 0.25°, in which the 17 land cover classification schemes were used based on the International Geosphere-Biosphere Program (IGBP) [36].

2.2. Methods

To assess the interannual and interdecadal variability of the ASMB, we first used time series analysis, discussed its subregional variability, and subsequently employed wavelet analysis to identify its multiple cycles of variability, reflecting the trends at different time scales. To better understand the corresponding dynamic mechanisms or external forcing, composite analysis was used to discuss the wind field at 500 hPa [37, 38] and also calculate the water vapor flux on the TP and its surrounding areas [3941]. The westerly wind index and the South Asian monsoon index were used to represent the intensity variation of the westerly wind and the South Asian monsoon. The cross wavelets were used to analyze the relationship between the westerly wind-ASMB and the South Asian monsoon-ASMB [42, 43]. The functions and their corresponding definitions are detailed below.

2.2.1. Wavelet Analysis

We selected the Morlet wavelet [44] to analyze the relationship between the South Asian monsoon-monsoon boundary-westerly winds. The Morlet wavelet function is expressed aswhere is the wavelet function; is the imaginary symbol of a complex number; and is the nondimensional frequency. The Morlet wavelet transform coefficients are defined aswhere and are the wavelet coefficients of the original and conjugate functions, respectively, of x; a and b represent the time scale and time position factors, respectively. The magnitude of the mode of the Morlet wavelet coefficient of variation indicates the strength of the signal of the characteristic time scale; the real part indicates the information of the distribution and the bit phase of different characteristic time scales at different times. The wavelet variance can be obtained by integrating the squares of all wavelet coefficients for different scales. The wavelet variance can indicate the strength of the fluctuation of that cycle in the time series and is defined as

2.2.2. Cross-Wavelet Analysis

To assess the influence of westerly winds and South Asian monsoon on ASMB, the cross-wavelet power spectrum (XWT) and the cross-wavelet coherence spectrum (WTC) were analyzed [44, 45]. The XWT can be defined using the two times series X and Y, with wavelet transforms and , where n is the time index and s is the scale:where is the complex conjugate of . The XWT can reflect two sequences with the same energy spectral region after wavelet transform, thereby revealing the significance of the interaction between two sequences in different frequency domains.

The WTC is defined aswhere is a smoothing operator. The WTC can reflect the coherence range of two wavelet transform in the time-frequency domain.

2.2.3. Westerly Winds Index

The westerly winds index (WWI) in the Northern Hemisphere is based on the Rossby [46] westerly wind index, and it is calculated based on the following equation:where indicates 35°N at 500 hPa height and indicates 55°N at 500 hPa height.

2.2.4. South Asian Monsoon Index

The South Asian monsoon index is calculated using the South Asian monsoon index (SASMI) based on Li and Zeng [47, 48]:where is the climate average wind vector for January, is the climate average wind vector for January and July, and is the one month (n) in a certain year (m).

3. Results

3.1. Compare the ASMB Definition in the TP

Till date, a considerable amount of research has been dedicated to defining and elucidating the variation of the ASMB. However, there are relatively few relevant studies on TP, and the applicability of these definitions of the ASMB in the TP needs to be discussed further. To this end, we first selected three representative definitions of the ASMB and compared their applicability in the TP (Figure 2). The selected representative boundary definitions include the following: (1) the summer (May–September) 2.0 mm day−1 isoline proposed by Chen et al. [25] as the monsoon boundary (“Chen” in Figure 2); (2) the global monsoon index proposed by Wang et al. [24], which is the region in which the difference between the summer (May–September) and winter (November–March of the next year) precipitation exceeded 2.0 mm day−1, and summer precipitation exceeded 55% of total annual precipitation (“Wang” in Figure 2); and (3) the monsoon boundary proposed by Tang et al. [23], as the 20 mm isoline of variance in summer precipitation from June to August (“Tang” in Figure 2).

Figure 2 shows that the three indicators are roughly in the southwest-northeast direction. However, there are also significant differences among these three indicators. For instance, the position of the indicator “Chen” is more northward than the other two indicators, the indicator “Tang” is the more southerly, and the indicator “Wang” is located between the other two. The land cover types in northern TP are barren and sparsely vegetated, with central and southeast TP being dominated by open shrublands and grasslands, respectively (Figure 2(a)). The position of the ASMB indicator “Chen” is located between sparsely vegetated areas, open shrublands, and grasslands; the indicator “Wang” crosses from open shrublands and grasslands; and the indicator “Tang” is located to the south of open shrublands and grasslands. The land cover types are more influenced by climatic conditions and also more sensitive in monsoon boundary areas [5052]. The fluctuations of the indicator “Chen” generally better reflect the spatial distribution of modern land cover types in the TP. In the arid-humid zoning map, which is divided based on annual average precipitation (Figure 2(b)) [49], the entire TP is roughly divided into arid (north), semiarid and semihumid (center), and humid areas (south). Figure 2(b) shows that the fluctuations of the indicator “Chen” essentially match the spatial distribution of the dividing line of arid and semiarid areas. However, as the indicators “Wang” and “Tang” are located in the semihumid and humid areas, respectively, such a monsoon boundary is unrealistic for the ASM. As for the surface wind field (Figure 2(c)), the ASMB depicted by indicator “Chen” is located in the northernmost position of the southerly wind, which is more compatible than the fluctuating situation presented by the other two indicators in the TP. Overall, the indicator “Chen” is better at portraying the fluctuations of ASMB in the TP and more meaningful from a climatic-ecological perspective [52]. Thus, the “Chen” 2 mm day−1 precipitation isoline during summer (May–September) was considered for further discussions on the ASMB.

Additionally, we analyzed the differences of the three indicators under different datasets using the CRU (Figures 2(d)2(f)) and CN05 datasets (Figures 2(g)2(i)). The results indicated that the three datasets (GPCC, CRU, and CN05) performed differently for the TP. The ASMB generated by the GPCC (Figure 2(a)) and CN05 (Figure 2(g)) datasets is more ecologically related and located at the boundary of sparsely vegetated areas-open shrublands-grasslands. In contrast, the ASMB generated by the CRU dataset (Figure 2(d)) does not correspond well to the modern land cover type map. In the arid-humid zoning map, the ASMB generated by the GPCC (Figure 2(b)), CRU (Figure 2(e)), and CN05 (Figure 2(h)) datasets are at the junction of semiarid and arid areas, and all exhibit good climatological significance. In the surface wind field map, the ASMB generated by the GPCC (Figure 2(c)) and CN05 (Figure 2(i)) datasets has a northerly position but is more southerly for the CRU (Figure 2(f)) dataset.

The CRU and GPCC global gridded land surface precipitation datasets are both widely used. Both datasets are based on observation data collected by global meteorological stations and have long time series (1901–present) and high horizontal resolution (0.5° × 0.5°) [52, 53]. The poorer performance of the CRU dataset can be attributed to the method used to process the data. For instance, in the early years, there was a lack of observation values in the TP; consequently, the climatic average was used instead of actual precipitation data [31, 32, 54]. The CN05 precipitation dataset is based on the latest precipitation data of 2472 stations on the ground in China, compiled by the National Meteorological Information Center. The thin-disk spline method of ANUSPLIN software was used for spatial interpolation to generate monthly precipitation data in China with the latest horizontal resolution of 0.5° × 0.5° from 1961 to 2013 [55]. The CN05 dataset is based on the regional dataset of China and shows only specific changes in the TP of the Chinese region. In general, the position of the ASMB generated by the GPCC dataset was superior in terms of delineating and detecting changes, indicating that is a better indicator of climatic-ecological significance. Subsequent analyses focused on the GPCC dataset.

3.2. Variation in the ASMB in the TP

Based on the GPCC dataset, the ASMB was illustrated from 1980 to 2019 using the indicator of summer precipitation of 2.0 mm day-1 to analyze interannual variability. Figure 3 shows that the ASMB generally exhibits a southwest-northeast trend, and any changes were mainly concentrated in the central region of the TP near the Tanggula Mountains and northeastern portion of the TP near the Qilian Mountains, extending as far north as 36N and as far south as 29N in the center part. Along the western part of the plateau, the ASMB changed minimally, especially at 80°E and primarily fluctuated between 32N and 34N. In the eastern part of the TP, the ASMB fluctuated between 36N and 40N. The fluctuations in the amplitude from west to the east suggested a “small-large-small” fluctuation situation.

When analyzing interdecadal variation in the ASM, Ding et al. [56] proposed that it exhibited interdecadal evolution that was first strong, then weak, and then strong after the 1970s. Figure 4 shows that the main fluctuations in interdecadal variation are concentrated in the central part of the TP over the four decades. Lower fluctuations were observed on the western and eastern sides of the TP, indicating an overall southwest-northeast trend. In the western side of the TP, the position of the ASMB extends as far north as 34°N. In contrast, in central TP, the ASMB extends as far north as 36°N. However, interdecadal variation in the eastern part of the TP is not obvious. Overall, interdecadal variability was low in the 1980s and 1990s (Figures 4(a) and 4(b)), but was large in the central and western parts of the TP in the 2000s and 2010s, where it was primarily concentrated (Figures 4(c) and 4(d)).

On analyzing the interannual and interdecadal variations, it can be observed that the ASMB shows different fluctuation situations in different regions of TP; to better reflect the fluctuation magnitude of ASMB and explore how it fluctuated in the TP from 1980 to 2019, three regions (68°E–80°E, 80°E–95°E, and 95°E–105°E) were delineated (Figure 5), based on the results of the interannual and interdecadal “small-large-small” fluctuations shown from west to east (Figures 3 and 4). For these regions, the temporal variation in the average annual latitude of the ASMB was calculated, and the results indicate that the fluctuation of ASMB was relatively smooth before 2000, gradually increasing after 2000 (Figure 5). In the easternmost region (95°E–105°E), the position of the ASMB was the most northerly as compared to the other regions (Figure 5(a)). In the central region (80°E–95°E), the position of the ASMB fluctuated the most, by up to 2.5° latitude in some years (Figure 5(b)). In the westernmost region (68°E–80°E), the position of the ASMB was the most southerly and had the lowest fluctuation (Figure 5(c)). The overall range in fluctuation was noticeably smaller in the western part of the TP (68°E–80°E). In summary, the fluctuation was primarily concentrated in the central part of the TP (80°E–95°E).

The results of changes in the position of the monsoon boundary were assessed periodically using wavelet analysis (Figure 6), and it was confirmed that the characteristics of the ASMB varied for these three regions. Figure 6(a) shows the existence of a 13a–16a periodic signal during 1998–2016 and a 3a periodic signal during 2003–2015, corresponding to the wavelet variance; Figure 6(b) shows the 3a as the main cycle. Figure 6(c) shows that there are clear periodic signals at 5a, 10a–15a, and 25a. Near the 5a scale, the ASMB exhibits a small cycle change in north-south fluctuations and a clear alternating north-south change. At 10a, it exhibits a more stable periodic variation of the signal and a significant interdecadal variation. At 25a, although there are obvious periodic signal variations, the contours are not closed and are not subject to discussion. The corresponding variance (Figure 6(d)) indicates that 4a and 10a are the first and second cycles, respectively. Figure 6(e)) shows 5a–6a and 14a–15a periodic signals, which mainly occurred in 1993–2019 and 1983–2015. The variance (Figure 6(f)) shows two distinct peaks, which correspond to the 6a and 15a scale periods, with the main cycle occurring at 15a and the second cycle occurring at 6a.

3.3. Cause of the ASMB in the TP

As discussed previously, the ASMB exhibits clear variation in the central part of the TP (80°E–95°E). In this area, there were large fluctuations in the position of the ASMB, with significant interannual variability. Consequently, variation in the position of the ASMB in this region might be a better representative of changes in the entire TP. To elucidate the mechanism driving this variation, five years, representing the most northerly and southerly years of the ASMB in the TP, were selected in the longitude range from 80°E to 95°E (based on Figure 5). The selected northerly years were 1999, 2000, 2003, 2008, and 2010. The selected southerly years were 1983, 1990, 2006, 2015, and 2019. The 500 hPa wind fields from July to August in the northerly and southerly years were then analyzed individually using composite analysis (Figure 7).

The southern areas of the TP were mainly controlled by the South Asian monsoon (Figures 7(a) and 7(c)), and in northerly years, more areas of the TP were affected, resulting in fluctuations in the ASMB in the TP. The summer wind intensity reflects the characteristics of the main body, and the location of the ASMB reflects the characteristics of the edge [25, 5759]. In particular, water vapor transport in South Asian monsoons strongly impacts precipitation in the TP (Figure 8). For instance, in a northerly year (Figure 7(b)), there was an anticyclone in the Bay of Bengal, which allowed the South Asian monsoon to affect the regions falling at about 35°N and helped transport warm and humid air to more areas of the TP. This phenomenon caused precipitation to increase in the region and caused the ASMB to shift northward. In a southerly year (Figure 7(d)), there was a strong anticyclone on the west side and a cyclone on the northeast side of the TP. This circulation situation made the westerly winds stronger, which allowed the westerly winds to propagate in the TP and created a situation wherein the South Asian monsoon could only affect the regions falling about 29°N; this reduced water vapor transport, and caused the ASMB in the TP to shift southward.

As shown in Figure 7, the relationship between the water vapor transport of the South Asian monsoon and the location of the ASMB is relatively clear. To examine the water vapor transport problem of the monsoon, two years that had the most northerly and southerly ASMB positions (2003 and 2015, respectively) were selected. Water vapor transport of the entire summer layer in the TP and its surrounding areas was then analyzed (Figure 8): water vapor transport in the TP in summer primarily involved northward zonal transport. There was a strong water vapor transport belt in the Indian Ocean–Bay of Bengal, with a maximum value of 64 × 10−5 kg/(m−2s). Moisture flowed from the Indian Ocean, through the Indian Peninsula to the Bay of Bengal, and then northward to the TP. It then converted to meridional transport under the influence of the westerly winds. In the northerly years (Figure 8(a)), the center of the water vapor has larger values on the west side of the Indian peninsula. The monsoon transported water vapor along the transport channel up to 32°N or even further north, which might have caused the position of the ASMB to exhibit a significant northward shift. In the southerly years (Figure 8(b)), the meridional transport is deflected at 28°N under the action of strong westerly winds, which might cause a southerly shift in the position of the ASMB.

By analyzing the water vapor flux transport of the monsoon, the effect on the position of the ASMB was elucidated. However, the climate of TP is also influenced by the westerly winds [60, 61]. Therefore, we examined how the westerly winds and South Asian monsoon impact the position of the ASMB (Figure 9). For our investigation, we selected the westerly winds index (WWI) in the Northern Hemisphere, which is based on the index proposed by Rossby [46], and the South Asian monsoon index, which is calculated using the index based on Li and Zeng [47, 48]. The anomaly of the ASMB was then calculated in the TP from 80°E to 95°E; moreover, the anomalies of the westerly wind index (WWI) and SASMI were also calculated (Figure 9).

As we can see in Figure 9, the ASMB changes with the alternation of westerly winds and South Asian monsoon. During strong westerly winds and weak South Asian monsoon years (e.g., 1982, 1990, 1994, and 2016), the ASMB was positioned further south in the TP (Figure 9). During weak westerly winds and strong South Asian monsoon years (e.g., 1985, 1994, 2001, 2006, and 2017), the ASMB was positioned further north in the TP. During years when the westerly winds and South Asian monsoon were both weak (e.g., 1983, 1992, and 2013) or both strong (1990, 1994), three-position patterns of ASMB were exhibited: a stable state, or a southward or northward deflection, indicating the combined influence of the westerly winds and South Asian monsoon during these periods. During 1980–1995, the South Asian monsoon and westerly winds changed alternately, with this pattern being reflected in the position of the ASMB (Figure 9). During 1996–2008, the ASMB anomaly was positive, and its position was more northerly; during this period, the SASMI anomaly was mostly positive and strongly fluctuated, whereas the WWI anomaly weakly fluctuated. However, there were still such periods when the South Asian monsoon (westerly winds) was strong, and yet changes to the westerly winds (South Asian monsoon) affected the position of the ASMB, which shifted to the south when the westerly winds (South Asian monsoon) became stronger, and to the north when the westerly winds (South Asian monsoon) weakened.

To improve our current understanding of the interactions between the westerly winds, the ASMB, and the South Asian monsoon, cross-wavelet analysis was used. According to the power of the cross-wavelet spectrum (Figure 10(a)), the spectral energy of the ASMB-westerly winds was concentrated in 1-2 a (2005–2008) and 9–11 a (1994–2003) cycles. The spectral energy of the ASMB-South Asian monsoon was concentrated in a 3–5 a (2013–2019) cycle (Figure 10(c)). However, we estimated the significance level for each scale only using values outside the COI [62]. According to the cross-wavelet coherence spectrogram, there were positive (1995–2000) and negative (2003–2006) correlations between the westerly winds and the ASMB (Figure 10(b)). The phase angle indicates that the westerly winds lagged behind the ASMB. For the cross-wavelet coherence spectrogram of the South Asian monsoon-ASMB (Figure 10(d)), negative (1995–2004) and positive (2010–2014) correlations were also observed. In summary, the phase difference between the westerly winds and the South Asian monsoon may indicate the impact of a synergistic effect, which could be amplified by interactions between the westerly winds and the South Asian monsoon.

Based on the calculation results of the cross-wavelet analysis, there may exist a synergistic effect between the westerly winds and South Asian monsoon. To examine this synergistic effect, precipitation anomaly maps for positive South Asian monsoon and neutral westerly winds events, neutral South Asian monsoon and positive westerly winds events, positive South Asian monsoon years, and positive westerly winds events were studied (Figure 11); the events with variability exceeding the +1 standard deviation of the SASMI are defined as the positive South Asian monsoon years. Events within +1 and −1 standard deviations for SASMI are neutral South Asian monsoon years [63]. Positive westerly winds and neutral westerly wind years are also defined in this manner. Thus, there are years with positive South Asian monsoon and neutral westerly winds events (1991, 2001, 2006, and 2017), neutral South Asian monsoon and positive westerly winds events (1998, 1999, 2002, and 2016), and positive South Asian monsoon and westerly winds events (1990 and 1994), which impact the position of the monsoon belt differently. In years with positive South Asian monsoon and neutral westerly winds events, the southern part of the TP received more precipitation and the ASMB was more northerly (Figure 11(a)). This positioning reflected the positive phase relationship between the South Asian monsoon and ASMB. During years with neutral South Asian monsoon and westerly winds events (Figure 11(b)), only the central-eastern part of the TP experienced more precipitation (and the other regions had less precipitation), with the ASMB being more southerly. This phenomenon reflected the inverse phase relationship between the westerly winds and ASMB. However, most of the regions, especially in the center part of TP, did not pass the 95% confident level; consequently, it can be inferred that the precipitation anomalies were not influenced by a single westerly wind or South Asian monsoon. In contrast, in years with positive South Asian monsoon and westerly winds events, a region of low precipitation occurred in the central part of TP, with the ASMB being more southerly (Figure 11(c)). This result indicates that while South Asian monsoon and westerly winds have positive and negative effects on the ASMB, the positioning of the ASMB was influenced by a synergistic effect between the South Asian monsoon and westerly winds.

4. Conclusion and Discussion

Changes in the monsoon affect more than two-thirds of the global population and are of great scientific and social importance [3]. The monsoon system generally has its own variable characteristics [5], the strengths and weaknesses of which are reflected in the variations of their boundaries [1921]. However, the definitions of different monsoon boundary indices vary considerably, and the performances of these indices are highly uncertain [2224]. The majority of previous studies have not accounted for the TP, owing to its pronounced topographic effect [26]. In this study, we compared three monsoon indicators and found that the summer (May–September) 2.0 mm day−1 isoline proposed by Chen et al. [25] better reflects variations in the ASMB on the TP. We found that this indicator has both climatic and ecological significance by combining modern land cover types, an arid-humid zone map, and the surface wind fields. In addition, we analyzed the effects of water vapor transport from the South Asian monsoon and the action of the westerly winds, for which we observed synergistic effects between the two.

With respect to climate change on the TP, the strength of the monsoon has an impact on precipitation on the plateau [60]. Compared with previous studies [52], we confirmed that the indicator proposed by Chen et al. [25] has climatic and ecological significance, and better reflects the variations in ASMB on the TP. Based on comparisons with the other two datasets, we found that the GPCC dataset is more applicable for the TP [3032]. The TP is also influenced by the westerly winds [61]. Recent studies have shown that the modulation of the westerlies-monsoon interaction could affect the climate along the monsoon boundary zone [57, 58], and our findings are in accordance with those of previous studies [59]. These findings improve our understanding of the monsoon activity patterns and climate change on the TP. However, this study only examined data from four decades, owing to the limited timescales of the datasets. Modern land cover types were used to verify the ASMB. It is also difficult to obtain reliable vegetation data that spans longer timescales. Our analysis indicates that the ASMB was affected by a synergistic effect between the South Asian monsoon and the westerly winds. The nature of these interrelated mechanisms should be explored further in future research.

Variations in ASMB on the TP were examined using datasets that span four decades. As demonstrated by the interannual and interdecadal results, these variations are concentrated in the central part of the TP region and exhibit a SW-NE trend. Spatial variations in the amplitudes of these fluctuations exhibited a small-large-small trend from south to north. The ASMB is mainly affected by the synergistic effect of the South Asian monsoon and the action of the westerly winds. Further research is required to determine how the monsoon boundary is affected by the relationship between these parameters. The physical mechanisms involved in this relationship, as well as changes that occur over longer timescales, also require further investigation.

Data Availability

The GPCC dataset can be downloaded from the Global Precipitation Climatology Centre (https://www.cgd.ucar.edu). The CRU dataset can be downloaded from the Climate Research Unit (https://crudata.uea.ac.uk/cru/data). The CN05 dataset can be obtained from the China Meteorological Data Service Centre (https://data.cma.cn). Monthly wind data were obtained from the monthly reanalyzed data of the National Center for Environmental Prediction/Atmospheric Research (https://psl.noaa.gov/data/gridded/reanalysis). The ear5 dataset can downloaded from the Copernicus Climate Change Service (https://www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era5). The modern land cover data can be downloaded from the Earth Observing System Data and Information System data center (https://daac.ornl.gov). Boundary data of the Tibetan Plateau were obtained from the National Tibetan Plateau Data Center (https://data.tpdc.ac.cn/).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the Special Survey Project sponsored by the Ministry of Science and Technology of the People’s Republic of China (Grant No. 2017FY100903) and the Sichuan Central Government Guides Local Science and Technology Development Project (No. 2021ZYD0019) and the Second Tibetan Plateau Scientific Expedition and Research (STEP) Program (No. 2019QZKK0103).