Next Article in Journal
Can Citizen Science Promote Flood Risk Communication?
Previous Article in Journal
Urban Hydroinformatics: Past, Present and Future
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial and Temporal Variations in the Rainy Season Onset over the Qinghai–Tibet Plateau

1
Institute of Remote Sensing and Earth Science, College of Science, Hangzhou Normal University, Hangzhou 311132, China
2
State Key Laboratory of Estuarine and Coastal Research, East China Normal University, Shanghai 200062, China
3
Institute of Agricultural Remote Sensing and Information Technology Application, College of Environmental and Resource Sciences, Zhejiang University, Hangzhou 310058, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(10), 1960; https://doi.org/10.3390/w11101960
Submission received: 23 July 2019 / Revised: 6 September 2019 / Accepted: 18 September 2019 / Published: 20 September 2019
(This article belongs to the Section Water, Agriculture and Aquaculture)

Abstract

:
Precipitation on the Qinghai–Tibet Plateau (TP) in southwestern China is subject to interactions between the complex and variable terrain and the sensitive climate. The regional climate is mainly affected by three circulations: westerlies, the South Asian monsoon, and the East Asian monsoon. Spatial and temporal variations in the rainy season onset were characterised based on daily precipitation from 106 meteorological stations on the TP from 1971 to 2015. Using the Theil–Sen Median trend analysis, Mann–Kendall test and mutation detection, the characteristics and reasons for the variations during the rainy season over the plateau over the past 45 years were investigated. The following results were obtained from the analysis: (1) There were obvious regional differences in the rainy season onset over the TP, and the rainy season began on the southeastern plateau and moved northwestward. (2) The TP rainy season underwent a significant mutation in approximately 1997, and following this mutation, the area affected by the delayed rainy season increased. (3) Against the background of global warming, the rainy season trend over the TP was advanced; however, there were still several multiple contiguous concentrated areas on the plateau. (4) Before the rainy season mutation, there were two centres of delayed precipitation on the plateau, which existed primarily due to their location at the end of the plateau water vapour transport channel. After the mutation, the number of delayed precipitation centres on the plateau increased to three and presented a spatially expanding trend, which may be related to the weakening trend in atmospheric circulation.

1. Introduction

Precipitation is a key part of the global atmospheric and hydrological cycles and an important factor that affects the distribution of water resources [1]. Climate warming has accelerated the water cycle by speeding the processes of precipitation, evapotranspiration, and surface runoff and has, thereby, changed the spatial and temporal distributions of water resources, which widely impact the distribution of freshwater resources and the stability of natural ecosystems [2,3]. Subtle changes in precipitation have profound effects on plateau ecosystems, energy cycles, agricultural production, food security, and so on [4,5,6,7]. With global warming, the temperatures of both continents and oceans have increased; however, the temperature has increased more slowly in the ocean due to its specific heat, reducing the sea–land thermal gradient, which has changed the characteristics of atmospheric circulation and caused variations in the timing and intensity of precipitation events in the region [1], [8,9,10,11,12].
As the “third pole” of the world, the unique ecological geography unit of the Qinghai–Tibet Plateau (TP) has been extensively studied in the climate response field, and the sensitivity of the climate of the TP is well documented. Thus, this plateau has been called the “driver” and “amplifier” of climate change [13,14,15]. Currently, precipitation characteristics, such as seasonal precipitation, precipitation intensity, and extreme precipitation changes on the plateau have been widely studied by researchers [16,17]. The rainy season is the most concentrated period of precipitation during the year on the TP, and the precipitation during this season can reach 90% of the annual total precipitation [18]. The onset of the rainy season serves as a signal that the plateau climate has changed from dry to wet. The onset of the rainy season is a direct reflection of the modal changes in atmospheric circulation. Studying the spatial and temporal changes in the rainy season onset plays an important role in informing agricultural and water security planning.
The climate of the TP (especially precipitation) is affected by a combination of the South Asian monsoon, westerlies, and the East Asian monsoon [19,20], and subtle changes in atmospheric circulation intensity can cause significant interannual variations in the rainy season onset over the TP. Several investigations have indicated a clear trend in the intensities of the westerlies and monsoon circulation. The TP climate tends towards warm and humid [11,21]. In the context of global climate and circulation changes, the temporal and spatial trends of the rainy season onset over the TP are still unknown. Therefore, in the present study, the temporal and spatial variation characteristics of the rainy season onset over the Tibetan Plateau were investigated, and the results may reveal the variation characteristics in the monsoon circulation system and shed light on the potential impacts of climate change on the global water cycle, which is also meaningful for environmental ecological risk assessments.
The daily precipitation data from 106 meteorological stations from 1971 to 2015 were analysed (Figure 1). Based on GIS spatial analysis, the Theil–Sen Median analysis method combined with the Mann–Kendall test was used to study the relationship between atmospheric circulation changes and the onset of the rainy season on the plateau against a global warming background. The findings can provide a reference for understanding the characteristics of the global water cycle and inform planning to avoid ecological risks.

2. Materials and Methods

2.1. Study Area and the Partitioned Regions

The Qinghai–Tibet Plateau (TP), with an area of 2.5 million km2, lies in southwestern China (73° E–105° E, 25° N–40° N). The TP is the highest-altitude region in the world. The terrain of the plateau is complex, and the climate over the plateau has obvious regional characteristics. In this paper, the study area covers Tibet in Qinghai Province and the counties of Sichuan, Yunnan, Xinjiang, and Gansu. The water vapour transport on the TP is mainly derived from three major water vapour channels [20]. The first is located in the middle of the Himalayas in the southwestern part of the TP, the second is in the Brahmaputra River valley in the southeastern part of the TP, and the third is in the Tarim Basin and Qaidam Basin in the northern part of the TP. The water vapour flux of the water vapour channel on the southern TP is mainly affected by the Indian summer monsoon. The water vapour flux of the northern water vapour channel is affected by westerlies. The atmospheric precipitation from the southeast water vapour channel of the plateau is obviously greater than that of the other channels, and there are three high-precipitation centres, which are located in the southwestern, southeastern, and southern parts of the TP [22]. There are obvious spatial and temporal differences in the plateau precipitation. In terms of time, there are obvious dry and wet seasons on the plateau, with precipitation mainly concentrated in the rainy season. In terms of space, the TP precipitation gradually decreases from southwest to northeast. On the western Sichuan Plateau, on the Qinghai Plateau, and south of Bayan Kala Mountain, the annual precipitation is greater than 400 mm, and this area is humid and semihumid. In contrast, the annual precipitation is less than 400 mm in the western and northern parts of the plateau, which are arid and semiarid regions.
As mentioned above, there are regional differences in precipitation over the plateau. In this paper, based on the climate division of the Chinese Academy of Sciences and by combining the topography, atmospheric circulation, watershed division, and regional precipitation on the plateau, we divided the plateau into 11 regions (Figure 1) to better understand the regional rainy season onset and trend over the TP (A. Northern Tibet, NT; B. Central Tibet, CT; C. Southern Tibet, ST; D. Qaidam Basin, QB; E. Southern Qinghai, SQ; F. Changdu, CD; G. Dawang–Chayu, DC; H. Qilian–Qinghai Lake, QQ; I. Bomi–Western Sichuan, BW; J. Northern Yunnan, NY; K. Sichuan, SC).

2.2. Datasets

The ground observations of 106 meteorological stations covering the TP were obtained from the National Meteorological Information Center of China [23]. In the time series, daily precipitation data with strict quality control from 1971 to 2015 (45 years) were selected to ensure the integrity of the datasets and the accuracy of the research results. The datasets were processed in a standard manner to eliminate invalid and erroneous data. We rearranged the original data according to the year and site. Suspicious data were compared to the missing value for quality control, and the missing data accounted for fewer than 5 parts per million, too low to influence the research results.
Shuttle Radar Topography Mission Digital Elevation Model (SRTM DEM) data was obtained from the United States Geological Survey (USGS) website, and the data has a spatial resolution of 30 m [24].

2.3. Method

2.3.1. Determining the Rainy Season Onset

The onset of the rainy season was determined according to the relative precipitation coefficient model [25,26]. The formula is shown as follows:
C n = A n A ,
where An is the average precipitation of day n, A is the average annual precipitation of the meteorological station, and Cn is the relative precipitation coefficient of day n. Generally, the precipitation coefficient is related to the local annual precipitation [26]. To determine the onset of the rainy season, we need to calculate the 5-day, 10-day, and 15-day relative precipitation coefficients. The relative precipitation coefficient was set as 1.5 to avoid the impact of extreme precipitation on the results.
Thus, the criteria for the rainy season onset are as follows: After a moderate rain day (i.e., daily precipitation greater than 5.0 mm), the three indicators C5, C10, and C15 for a relative precipitation coefficient are simultaneously greater than 1.5; that is, this rainy day can be regarded as the rainy season onset.

2.3.2. Theil–Sen Median Trend Analysis

The average annual precipitation clearly differed among the partitioned regions of the TP, and the frequency of extreme precipitation (precipitation events greater than the 95th percentile) on the plateau has increased over the past 50 years [17,27]. The Qaidam Basin and its surrounding area, which is between the Kunlun Mountains and the Qilian Mountains, is a typical arid and semiarid region. The annual average precipitation was less than 400 mm, and a slight change in precipitation may affect the determination of the rainy season onset. The Theil–Sen Median analysis is a nonparametric trend analysis method that is less susceptible to extreme values and is more robust than a traditional parametric analysis [28]. The slope of the rainy season onset was calculated by the following formula:
S ONSET = Median ( ONSET j ONSET i j i ) , 1971 i j 2015 ,
where i and j represent the observation year, ONSET i and ONSET j represent the rainy season onset (Julian Day) for the corresponding year from 1971 to 2015, and S ONSET represents the median value of the slope of any two years in n years (a total of C n 2 sets of data). If S ONSET was less than 0, the trend of the rainy season onset was delayed; if S ONSET was equal to 0, the trend of the rainy season onset had not changed at the corresponding temporal scale; if S ONSET was greater than 0, the trend of the rainy season onset was advanced. In practice, this value of S ONSET is usually magnified by 10-fold, which is called the ten-year rate of change and is in units of d/10 a.

2.3.3. Mann–Kendall Test

The Mann–Kendall (MK) test is a nonparametric statistical method [29,30]. The MK test is a time series method that is commonly used in meteorology and hydrology; this method has the advantages of not being affected by sample values, distribution types, and so on [31,32,33]. When the MK method is used to obtain a significant determinant of the change trend in the rainy season onset, the rainy season onset of the long-term sequence can be regarded as a set of independent sample data. When n is greater than 8, the standardised statistical quantity Z obeys the normal distribution. The calculation is as follows:
Z = {   S 1 var ( S ) ,     S > 0     0 ,     S = 0   S + 1 var ( S ) ,     S < 0   .
In Equation (3), S is calculated with Equation (4), and the variance var(s) is computed with Equation (5).
S = i = 1 n 1 j = i + 1 n sign ( ONSET j ONSET i ) ,
var ( S ) = 1 18 [ n ( n 1 ) ( 2 n + 5 ) p = 1 q ( t p 1 ) ( 2 t p + 5 ) ] ,
sign ( ONSET j ONSET i ) = { 1   ONSET j ONSET i > 0 0   ONSET j ONSET i = 0 1   ONSET j ONSET i < 0 .
In Equation (6), the sign is a function symbol, q represents the number of ties (repeated data) in the current time series, p represents a tie, and t p is the width of a tie (the number of times the same data appears). For a given α , if | Z | Z 1 α / 2 , there is a significant trend in the rainy season onset. In this paper, the confidence level α is taken as 0.1, and the corresponding |Z| is 1.645.

2.3.4. Moving t-test

Due to the mutation detection limitations of the MK test [33,34,35], a moving t-test was applied to test the time point of precipitation mutation over the TP, and the formula is shown as follows:
q = ( y 1 ¯ y 2 ¯ ) / ( s 1 n 1 + 1 n 2 ) ,
s = ( n 1 s 1 2 + n 2 s 2 2 ) / ( n 1 + n 2 2 ) ,
During the inspection process, a reference time point needs to be set. The two subsequences before and after this time are y 1 and y 2 , respectively. The sample sizes of the two subsequences are n 1 and n 2 , and the variances in the two subsequences are s 1 2 and s 2 2 . In addition, q is the statistic that follows the t distribution of degrees of freedom v = n 1 + n 2 2 . In this study, the sliding method was used, and the statistical array q i can be obtained by sequentially recursing the reference point. For a given level of significance α = 0.05 , there was a corresponding q α . If | q i | is greater than | q α | , the set reference time point can be regarded as a mutation point.

2.3.5. Cumulative Anomaly Model

For precipitation data, the cumulative anomaly model can be used as a detector for the mutation in the time series. For the cumulative anomaly value at a certain time t in the sequence, the curve formed by time is the cumulative anomaly curve. The formula is shown as follows:
D P i = i = 1971 i = t ( P i P ¯ ) ( 1971 < t 2015 ) .
In Equation (9), P i represents the annual precipitation in year i and P ¯ represents the mean annual precipitation from 1971 to 2015. The cumulative anomaly curve reflects the change in annual precipitation. The rising part is the precipitation increasing stage, and the falling part is the precipitation decreasing stage. The turning point of the curve can be regarded as the mutation year of precipitation.

3. Results

3.1. Spatial Distribution Characteristics of the Rainy Season Onset over the TP

The mean rainy season onset (1971–2015) of each meteorological station was calculated based on the daily precipitation data (Table 1). The average rainy season onset at most stations on the plateau was in the range of 145 ± 25 days, i.e., the onset was between late May and mid-June. At this time, it was summer on the TP, and the precipitation was related to the Indian summer monsoon. During this period, the plateau precipitation was mainly obtained through the water vapour channel of the Brahmaputra River in the southeastern plateau, which originated from the Indian Ocean [6,36].
An onset date contour was generated using the spatial interpolation model in ArcGIS 10.5 (Figure 2). The onset of the rainy season of the TP increased from south to north and from east to west with the advance of the rain belt, and the overall pattern showed a trend from southeast to the central and western regions.
The annual average rainy season onset on the plateau and its surrounding areas first appeared at the Gongshan station in Yunnan (the 62nd day) on the west side of the Hengduan Mountains. The last appeared at the Shiquanhe Station on the Ali Plateau in the western TP (the 194th day) and the Lenghu station (the 194th day) on the northwestern edge of the northern Qaidam Basin, a difference of 132 days. The rainy season started in March based on six stations (Gongshan, Weixi, Nielamu, Chayu, Pulan, and Bomi). Among these stations, four are located in the Hengduan Mountain region. The precipitation in this area is controlled by winter westerlies and the Indian summer monsoon, and the annual precipitation has a bimodal pattern, i.e., the region has a peak of precipitation in spring and a second peak in summer, and the precipitation is more abundant in summer [37]. The southwestern part of the Hengduan Mountains is a vast plains region. The low-altitude westerlies in the early spring are blocked by the TP after passing through the plains, which causes an early precipitation centre at the end of February and beginning of March in the area west of the Hengduan Mountain region. The rainy season onset over the plateau is closely related to the terrain and mountain ranges. The main mountain ranges in the northern part of the TP are east–west oriented, which provides the main water vapour channel via the midlatitude westerlies. The southern part of the TP is mainly the valley plain of the Brahmaputra River basin, which is located between the Gangdese Mountains and the Himalayas and is mainly influenced by the Indian summer monsoon. The Hengduan Mountains are in the southeastern part of the TP, and the main mountain ranges are north–south, which provides a south–north water vapour channel for the Indian summer monsoon [37].
The existence of the Brahmaputra River valley provides a favourable water vapour transport channel for the Indian summer monsoon. The monsoon carries sufficient water vapour upstream along the valley, making the rainy season onset in the middle and lower reaches of the river begin earlier than in the upper reaches of the river. The coastal areas of the river valley are extremely obvious, e.g., the rainy season onsets in Jiacha (the 150th day) and Milin (the 124th day) were obviously earlier than those in Lazi (the 164th day) and Rikaze (the 167th day).
The southeastern part of the TP is restricted by the Hengduan Mountains, and the water vapour transportation is limited. The area west of the Hengduan Mountains is mainly affected by the Indian summer monsoon. The humid airflow released heat and caused precipitation when climbing the windward slope, causing the atmosphere in the central region of the Hengduan Mountains to have a higher temperature but low humidity, which is a unique regional geographical phenomenon—a dry, hot valley. Thus, the rainy day onsets in the central Hengduan Mountains, e.g., Mangkang (the 150th day), Batang (the 146th day) and Derong (the 163rd day), were significantly later than those in the western Hengduan Mountains, e.g., Chayu (the 78th day) and Bomi (the 86th day).

3.2. The Multiyear Trend of the Rainy Season Onset over the TP

3.2.1. Characteristics of Onset Variations in the Rainy Season on the TP

The annual rainy season onset on the TP was calculated by all meteorological stations on the plateau (Figure 3). The year of the latest onset of the rainy season on the plateau was 1979 (the 159th day), and the year of the earliest was 2011 (the 122nd day). The results of the Theil–Sen Median trend method showed that the rainy season onset over the TP showed a relatively significant oscillating advancing trend, with an amplitude advance of 3.68 d/10 a.
The trend in the rainy season onset on the plateau was expressed on a time scale: in the 1980s, the rainy season onset was significantly later than the plateau average, and there was a delayed trend; in the late 1980s, there was a large step advance in the rainy season onset, after which the trend was relatively flat and showed an early trend.

3.2.2. Characteristics of Regional Variations in Rainy Season Onset

On a 45-year time scale, the rainy season change trend in each climate region was similar to the overall trend on the plateau (Figure 4). The difference in the rainy season onset of each region was as follows: A. Northern Tibet (advanced 7.85 d/10 a) > E. Southern Qinghai (advanced 6.76 d/10 a) > D. Qaidam Basin (advanced 5.28 d/10 a) > K. Sichuan (advanced 5 d/10 a) > B. Central Tibet (advanced 4.94 d/10 a) > C. Southern Tibet (advanced 4.1 d/10 a) > H. Qilian–Qinghai Lake (advanced 3.91 d/10 a) > I. Bomi–Western Sichuan (advanced 2.98 d/10 a) > J. Northern Yunnan (advanced 2.81 d/10 a) > F. Changdu (advanced 2.45 d/10 a) > H. Dawang–Chayu (delayed 4.14 d/10 a). The advancing trend of the rainy season onset in the southeastern part of the TP was lower than the average level of the TP (advanced 3.68 d/10 a), while the advancing trend of the rainy season onset in the central and western parts of the TP was more obvious. The occurrence of this phenomenon has a certain relationship with the topography and circulation of the plateau. The plateau monsoon region represented by the southeastern plateau was mainly affected by the Indian summer monsoon, where the advancing trend was slight. While the other regions were controlled by westerlies, the advancing trend was significantly greater than that of the plateau overall.
In general, the variation trend in the onset of the rainy season on the plateau was consistent with the overall trend on the plateau. The advanced trend of the rainy season onset in the southern monsoon region was moderate, but the interannual variation was more obvious. The variation in the rainy season onset in the westerly region was more obvious, but the interannual variation was slight. The more obvious variation in the monsoon circulation than in the westerlies might explain this phenomenon to some extent; however, there was a delayed trend in the monsoon circulation onset, which caused the advancing trend in the rainy season onset in the monsoon region to slow down.

3.3. Differences in the Variations in the Rainy Season Onset before and after Mutation

3.3.1. Mutation Analysis

The process in which the climate changes from one stable mode to another over a short period of time is called climate change. If the onset of the rainy season over the TP mutated over the course of 45 years, a key part of understanding the water cycle on the plateau is determining whether there was a difference before and after the mutation. The present study showed that there were significant differences in the rainy season onset in different regions, and the maximum difference was more than three-fold.
Annual precipitation is the most important climatic factor and represents the intuitive measure of regional precipitation. The change in precipitation in a long-term sequence is an indicator of regional climate change. Therefore, annual precipitation was used as an indicator in this paper, and the moving t-test model and cumulative anomaly model were used to evaluate the precipitation mutation (Figure 5).
The results of the moving t-test showed that when the step size was set to 6/7/8 (i.e., n1 = n2 = 6/7/8), the plateau precipitation was changed in 1997 and passed the 0.05 significance level test; the cumulative anomaly method showed that the cumulative anomaly of the plateau precipitation exhibited a downward trend during the period 1971–1997, reaching a minimum in 1997, while the cumulative anomaly exhibited an upward trend after 1997, which indicated that the precipitation had mutated in 1997. Therefore, 1997 could be regarded as the mutation year for precipitation and the onset of the rainy season on the TP. In the fields of marine science and atmospheric science, 1997 is also considered an important year of mutation. In marine science, 1997 is regarded as a year of abrupt change in the upper ocean temperature in the north–central Pacific [38]. With respect to the El Niño [39], the El Niño event in 1997 is considered to have been the strongest in the past 100 years and to have substantially changed the rainfall pattern in the eastern Pacific. Research on decadal changes in East Asian summer monsoon circulation also revealed that the monsoon underwent an abrupt change in the mid-1990s [40].

3.3.2. Differences in the Regional Variations in the Rainy Season Onset before and after Mutation

The differences in the rainy season onset over the TP before and after mutation were analysed with the Theil–Sen Median method (Table 2). Before the mutation, the trend in the rainy season onset at most meteorological stations over the TP was advancing (70.1%); 62.3% of stations exhibiting an advancing trend showed a significant trend, and the stations with a delaying trend mainly appeared in the Brahmaputra Valley of the southern Tibet region and the eastern Bomi–Western Sichuan region. After the precipitation mutation on the plateau, the trend in the rainy season onset on the plateau changed significantly; the stations with an advancing trend decreased to 50.5%, and only 18.8% of the stations showed an advancing trend. The proportions of stations with a delaying trend and unchanged stations increased significantly, which indicated that the trends at some stations had changed from the advanced mode to the invariant or delayed mode. Most of those stations were located in the Qaidam Basin region, Qilian–Qinghai Lake region, and Bomi–Western Sichuan region.

3.3.3. Differences in the Station Variations in the Rainy Season Onset before and after Mutation

According to the mutation test results, the study period was divided into two parts (before the mutation and after the mutation), and the trend in the rainy season onset for each station was also different before and after the mutation (Figure 6). Before the plateau precipitation mutation, the trend in the rainy season onset over most of the plateau was advanced and mainly located in the A/B/D/E/F/H/I region. Almost all of the stations were located in the plateau hinterland, while the stations in the southeastern part of the plateau (east of F region, central I region, and K region) and southwest of the plateau (C region) showed a delayed trend. The spatial and temporal variations in the trend in the rainy season onset had changed significantly after the mutation. The proportions and difference between advanced stations and delayed stations were reduced (before mutation: 70.93% and 26.74%, respectively; after mutation: 50.53% and 34.74%, respectively), while most of the stations at the edge of the plateau showed a delayed trend. The stations in the southwestern part of the TP (C region) still showed a delayed trend; most stations in the southeastern, eastern and northeastern parts of the TP showed delayed trends, and only four stations (I region) showed an unchanged trend.
Whether before or after mutation, the delayed stations were mostly at the edge of the TP; there was some spatial overlap, but differences remained. After the mutation, the proportion of delayed stations increased, and stations in the southern TP were located in the upper reaches of the Brahmaputra River (C region). However, after the mutation, the area of delayed onset expanded to the Central Tibet region (B region) and eastward to the southern Hengduan Mountains (J region and K region). The delayed stations in the area east of the TP were mainly located in the area north of the Hengduan Mountains (I region) before the mutation; however, the delayed centre underwent a significant northward trend after the mutation. In addition, after the mutation, a new delayed centre appeared in the area northeast of the TP (D region and H region), and the station trends in this region changed significantly from advanced (10–20 d/10 a) to delayed (within 10 d/10 a).

4. Discussion

The global and regional atmospheric circulation was changed by climatic warming, resulting in significant changes in the temporal and spatial patterns of climatic factors, such as the rainy season and precipitation on the plateau. The plateau climate was mainly controlled by westerlies and the South Asian monsoon, changes in the strength of which could directly lead to changes in climatic factors. The southern and southeastern parts of the TP are mainly controlled by the South Asian monsoon, the northern and central parts of the plateau are mainly affected by westerlies, and the eastern margin of the plateau is affected by the East Asian monsoon. In April, all of the water vapour on the TP mainly originated from the zonal transportation of westerlies. With the onset of summer, the influence of the South Asian monsoon expanded from south to north, and a large amount of water vapour from the Indian Ocean and Bay of Bengal entered the plateau from the main water vapour channel of the plateau. During the main rainy season on the plateau, the precipitation in most of the southeastern and southwestern regions was regulated by the South Asian monsoon, and the northern part was still controlled by westerlies.

4.1. The Trend in the Rainy Season Onset of the Main Plateau was Advanced

Whether before or after mutation, the rainy season onset in most regions of the TP showed an advancing trend. Before mutation, the region of advanced onset over the plateau was mainly located in the central and northern parts of the plateau (north of 30° N) and the southern part of the Hengduan Mountains, which is in the southeastern part of the plateau. After the mutation, the trend in the rainy season onset was still advanced over most of the plateau, but the area was significantly reduced, mainly in regions E, F, and I.
The trend in the rainy season onset was mainly related to the type of atmospheric circulation and its interannual variation on different parts of the plateau. The main mountain ranges in the central and northern parts of the plateau are east–west oriented. In the area north of the Gangdese Mountains and Nyainqentanglha Mountains, the main water vapour channel for the westerlies was precipitation and the rainy season, which were both related to the strength of the westerlies. The Brahmaputra Valley was the main channel for the Indian summer monsoon, and in this region, the rainy season onset was related to the strength of the Indian summer monsoon.
In addition, the abundant precipitation carried by westerlies in the winter was blocked by the mountains on the southwestern margin of the plateau. Therefore, an obvious precipitation peak appeared in these areas from January to April, which can explain the phenomenon in which the rainy season onset in Pulan (the 78th day) and Nielamu (the 73rd day) was significantly earlier than in other regions of the plateau.

4.2. The Trend in the Rainy Season Onset was Delayed in Some Regions over the TP

Though the advancing trend in the rainy season onset was apparent over the plateau, several concentrated stations showed a delayed trend. Before the mutation, the concentrated stations appeared in two different regions: One was in the southern part of the TP (C region), which was mainly distributed in the Brahmaputra Valley, where the precipitation is mainly affected by the South Asian monsoon, and the second was located at the northern end of the Hengduan Mountains and south of Bayan Hara Mountain (I region and F region), which was the boundary of the zone affected by the monsoon. After the mutation, the number of delayed regions increased to three. Among these regions, the delayed region in the Brahmaputra Valley extended northward, reaching the area north of the B region, and the delayed region also extended eastward to the southeastern margin of the plateau. The spatial variation characteristics of the delayed region in the Hengduan Mountains showed meridional and northward extensions (covering the central and northern area of the I region). In addition, in the northeastern part of the TP, a new region of delayed onset appeared (D region and H region), and the trend in this region was advanced before mutation but delayed after the mutation.

4.2.1. The Delayed Region South of Tibet

Before the mutation, the region of delayed onset in the southern part of Tibet was mainly located in the valley region of the middle reaches of the Brahmaputra (i.e., south of the Gangdese Mountains–Nyinqing Tanggula Mountains and north of the Himalayas). The annual precipitation in this region was between 250 mm and 450 mm, contributing to a semiarid and semihumid region. However, after the mutation occurred, the delayed region expanded across all of southern Tibet.
Before mutation, the delayed region in the Brahmaputra Valley was mainly related to the geographical location. The region was located between the two large east–west-oriented mountains on the TP, and the valley south of the plateau serves as the largest water vapour channel of the South Asian monsoon. The South Asian monsoon moves along the Brahmaputra, bringing abundant precipitation to the interior of the plateau. The southerly wind can carry water vapour to approximately the 30° N region on the TP, which was the centre of the delayed region. The delayed region was located at the end of the South Asian monsoon water vapour channel, and the atmospheric precipitation in this area was lacking, which led to the trend in which the onset of the rainy season differed in this region from that in other regions on the TP.
After the mutation, the delayed region had obviously extended south of the TP. The boundary of the extension region reached the B region in the north and the Hengduan Mountains in the east, including the G region and the area south of the I region. The expansion of the delayed region after the mutation was mainly related to the weakening of the Indian summer monsoon, which led the annual precipitation in this region to drop sharply, and the maximum decrease was up to 20 mm/year. Therefore, the main reason for the expansion of the delayed region south of the TP was that under the overall pattern of the weakening Indian summer monsoon, the water vapour carried by the monsoon was reduced, resulting in a decrease in the atmospheric precipitation in the region, which directly led to a delayed trend in those regions.

4.2.2. The Delayed Region in the Hengduan Mountains

Another region of delayed onset was concentrated in the I region. Before the mutation, the delayed region was mainly located north of the Hengduan Mountains and south of the Bayan Hara Mountains, which represented the northern boundary of the monsoon region. The delayed region tended to move northward and was located north of the Hengduan Mountains, which was mainly affected by westerlies.
Before the mutation, the rainy season in the delayed region began during summer, and the precipitation during this period was related to the geographical location and monsoon circulation. In the global atmospheric circulation pattern, the Hengduan Mountains served as the natural boundary between the East Asian monsoon and the South Asian monsoon (Indian summer monsoon), and the climate was affected by both the East Asian monsoon and the South Asian monsoon. The Indian summer monsoon has been weakening in recent decades. The river canyon between the Hengduan Mountains acted as a water vapour transport channel, but the shape of the mountains results in a complex folding of this channel (like a tie), which could block and restrain the transport of water vapour to the north. Therefore, at the end of the water vapour channel in the Hengduan Mountains, the precipitation was less abundant than in the southern part, and the onset of the rainy season showed a delayed trend.
After the mutation, the centre of the delayed region was obviously northward and mainly located to the east of the TP, which was at the end of the westerlies and monsoon. After 1998, the precipitation in this region decreased significantly, the annual onset of the rainy season in this region was on the 128th day, and the precipitation from March to May had a decrease of 0 to 5 mm/a [16], which suggested that the reduction in precipitation during the early rainy season may be an important reason for the delayed trend in the rainy season in the region.

4.2.3. The Delayed Region Northeast of the TP

The third delayed region was located in the Qaidam Basin and its northeastern part (D region and part of the H region); the annual precipitation in this region was lower than that in any other region, and the climate was dry. The rainy season onset was in June, which was later than that on the TP.
The change in the rainy season onset in this region from before to after the mutation was clear. Before the mutation, the rainy season onset in the region showed an advancing trend, and the advance range was 10 to 20 d/10 a. However, after the mutation, the rainy season onset showed a delaying trend, and the delay was within 10 d/10 a. This phenomenon was mainly related to the significant reduction in regional precipitation. The annual precipitation in this area was mainly concentrated from June to September, which coincides with the rainy season onset. The precipitation showed a decreasing trend after 1998, which led to a new delayed region northeast of the TP after mutation.

5. Conclusions

The Theil–Sen Median trend analysis and Mann–Kendall test could provide accurate trends in the onset of the rainy season. The temporal and spatial characteristics of the rainy season onset over the TP were analysed from 1971 to 2015, and some of the conclusions are as follows:
There were obvious spatial characteristics in the rainy season onset; the onset time increased with a trend from south to north along the meridional direction and from east to west with latitude. Due to the influence of the monsoon, the rainy season onset on the southwestern edge of the TP was earlier than in other regions of the TP.
The rainy season onset experienced a significant mutation in 1997. Whether before or after mutation, the trend in the rainy season onset over the TP showed an overall advance and partial delay, but the region with a delayed onset was larger after the mutation.
Before the mutation, the delayed region was mainly located at the end of the water vapour channel (Brahmaputra Valley and Hengduan Mountains), where there was less water vapour flux. After the mutation, the area of delayed onset south of the TP expanded significantly, and the weakening Indian summer monsoon was the main reason. The delayed region east of the TP moved northward, which was related to the weakening westerlies. The new delayed region northeast of the TP was related to the decrease in precipitation.

Author Contributions

J.X. and Y.H. (Yuekai Hu) conceived the study, and all authors contributed to writing the paper.

Funding

This study was supported by the Zhejiang Provincial Natural Science Foundation of China (No. LY16D010007), the Hangzhou Science and Technology Plan Project of China (No. 20170533B01) and the National Natural Science Foundation of China (No. 41571339).

Acknowledgments

The authors are grateful to editors and anonymous reviewers for their constructive comments and suggestion, which highly improving the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Li, L.; Schmitt, R.W.; Ummenhofer, C.C. The role of the subtropical North Atlantic water cycle in recent US extreme precipitation events. Clim. Dyn. 2018, 50, 1291–1305. [Google Scholar] [CrossRef]
  2. Tang, L.; Duan, X.; Kong, F.; Zhang, F.; Zheng, Y.; Li, Z.; Hu, S. Influences of climate change on area variation of Qinghai Lake on Qinghai-Tibetan Plateau since 1980s. Sci. Rep. 2018, 8, 3–9. [Google Scholar] [CrossRef]
  3. Zhang, W.; Yi, Y.; Kimball, J.S.; Kim, Y.; Song, K. Climatic controls on spring onset of the Tibetan Plateau grasslands from 1982 to 2008. Remote Sens. 2015, 7, 16607–16622. [Google Scholar] [CrossRef]
  4. Shen, M.; Tang, Y.; Chen, J.; Zhu, X.; Zheng, Y. Influences of temperature and precipitation before the growing season on spring phenology in grasslands of the central and eastern Qinghai-Tibetan Plateau. Agric. For. Meteorol. 2011, 151, 1711–1722. [Google Scholar] [CrossRef]
  5. Sun, J.; Qin, X. Precipitation and temperature regulate the seasonal changes of NDVI across the Tibetan Plateau. Environ. Earth Sci. 2016, 75, 1–9. [Google Scholar] [CrossRef]
  6. Yang, K.; Wu, H.; Qin, J.; Lin, C.; Tang, W.; Chen, Y. Recent climate changes over the Tibetan Plateau and their impacts on energy and water cycle: A review. Glob. Planet. Chang. 2014, 112, 79–91. [Google Scholar] [CrossRef]
  7. Ye, J.S.; Reynolds, J.F.; Sun, G.J.; Li, F.M. Impacts of increased variability in precipitation and air temperature on net primary productivity of the Tibetan Plateau: A modeling analysis. Clim. Chang. 2013, 119, 321–332. [Google Scholar] [CrossRef]
  8. Bollasina, M.A.; Ming, Y.; Ramaswamy, V. Anthropogenic Aerosols and the Summer Monsoon. Science 2011, 334, 502–505. [Google Scholar] [CrossRef]
  9. Duan, K.; Yao, T.; Thompson, L.G. Response of monsoon precipitation in the Himalayas to global warming. J. Geophys. Res. Atmos. 2006, 111, 1–8. [Google Scholar] [CrossRef]
  10. Manabe, S.; Terpstra, T.B. The effect of moutains on the general circulation of the Atmosphere as Identified by Numerical Experiments. J. Atmos. Sci. 1974, 31, 3–42. [Google Scholar] [CrossRef]
  11. Roxy, M.K.; Ritika, K.; Terray, P.; Murtugudde, R.; Ashok, K.; Goswami, B.N. Drying of Indian subcontinent by rapid Indian ocean warming and a weakening land-sea thermal gradient. Nat. Commun. 2015, 6, 7423. [Google Scholar] [CrossRef]
  12. Wu, B. Weakening of Indian summer monsoon in recent decades. Adv. Atmos. Sci. 2005, 22, 21–29. [Google Scholar] [CrossRef]
  13. Liu, X.; Chen, B. Climatic warming in the Tibetan Plateau during recent decades. Int. J. Climatol. 2000, 20, 1729–1742. [Google Scholar] [CrossRef]
  14. Liu, X.; Kutzbach, J.E.; Liu, Z.; An, Z.; Li, L. The Tibetan Plateau as amplifier of orbital-scale variability of the East Asian monsoon. Geophys. Res. Lett. 2003, 30, 1–4. [Google Scholar] [CrossRef]
  15. Pan, B.; Li, J. Qinghai-Tibetan Plateau: A Driver and Amplifier of the global climatic change. J. Lanzhou Univ. (Nat. Sci.) 1996, 32, 108–115. [Google Scholar]
  16. Ma, Z.; Xu, Y.; Peng, J.; Chen, Q.; Wan, D.; He, K.; Li, H. Spatial and temporal precipitation patterns characterized by TRMM TMPA over the Qinghai-Tibetan plateau and surroundings. Int. J. Remote Sens. 2018, 39, 3891–3907. [Google Scholar] [CrossRef]
  17. Zhang, K.; Pan, S.; Cao, L.; Wang, Y.; Zhao, Y.; Zhang, W. Spatial distribution and temporal trends in precipitation extremes over the Hengduan Mountains region, China, from 1961 to 2012. Quat. Int. 2014, 349, 346–356. [Google Scholar] [CrossRef]
  18. Shi, Y.; Song, L. Spatial Downscaling of Monthly TRMM Precipitation Based on EVI and Other Geospatial Variables Over the Tibetan Plateau From 2001 to 2012. Mt. Res. Dev. 2015, 35, 180–194. [Google Scholar] [CrossRef]
  19. An, Z.; Kutzbach, J.E.; Prell, W.L.; Porter, S.C. Evolution of Asian monsoons an phased uplift of the Himalaya-Tibetan Plateau since Late Miocene times. Nature 2001, 411, 62–66. [Google Scholar] [CrossRef]
  20. Yao, T.; Thompson, L.; Yang, W.; Yu, W.; Gao, Y.; Guo, X.; Joswiak, D. Different glacier status with atmospheric circulations in Tibetan Plateau and surroundings. Nat. Clim. Chang. 2012, 2, 663–667. [Google Scholar] [CrossRef]
  21. Li, L.; Yang, S.; Wang, Z.; Zhu, X.; Tang, H. Evidence of Warming and Wetting Climate over the Qinghai-Tibet Plateau. Arct. Antarct Alp. Res. 2010, 42, 449–457. [Google Scholar] [CrossRef]
  22. Zhou, S.; Wu, P.; Wang, C.; Han, J. Spatial distribution of atmospheric water vapour and its relationship with precipitation in summer over the Tibetan Plateau. J. Geogr. Sci. 2012, 22, 795–809. [Google Scholar] [CrossRef]
  23. National Meteorological Information Center of China. Available online: https://earthexplorer.usgs.gov (accessed on 20 September 2019).
  24. United States Geological Survey. Available online: http://data.cma.cn/data/cdcdetail/dataCode/SURF_CLI_CHN_MUL_DAY_V3.0.html (accessed on 20 September 2019).
  25. Fan, H.; Hu, J.; He, D. Trends in precipitation over the low latitude highlands of Yunnan, China. J. Geogr. Sci. 2013, 23, 1107–1122. [Google Scholar] [CrossRef]
  26. Thomas, A. The onset of the rainy season in Yunnan province, PR China and its significance for agricultural operations. Int. J. Biometeorol. 1993, 37, 170–176. [Google Scholar] [CrossRef]
  27. Cao, L.; Pan, S. Changes in precipitation extremes over the Three-River Headwaters region, hinterland of the Tibetan Plateau, during 1960–2012. Quat. Int. 2014, 321, 105–115. [Google Scholar] [CrossRef]
  28. Theil, H. A Rank-Invariant Method of Linear and Polynomial Regression Analysis, I, II, III. Nederl. Akad. Wetensch. Proc. 1950, 53, 386–392. [Google Scholar]
  29. Mann, H.B. Non-Parametric Test Against Trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
  30. Kendall, M.G. A new measure of rank correlation. Biometrika. 1938, 30, 91–93. [Google Scholar] [CrossRef]
  31. Shahid, S. Rainfall variability and the trends of wet and dry periods in Bangladesh. Int. J. Climatol. 2010, 30, 2299–2313. [Google Scholar] [CrossRef]
  32. Yue, S.; Pilon, P.; Cavadias, G. Power of the Mann-Kendall and Spearman’s rho tests for detecting monotonic trends in hydrological series. J. Hydrol. 2002, 259, 254–271. [Google Scholar] [CrossRef]
  33. Yue, S.; Wang, C.Y. The Mann-Kendall test modified by effective sample size to detect trend in serially correlated hydrological series. Water Resour. Manag. 2004, 18, 201–218. [Google Scholar] [CrossRef]
  34. Hamed, K.H.; Rao, A.R. A modified Mann-Kendall trend test for autocorrelated data. J. Hydrol. 1998, 204, 182–196. [Google Scholar] [CrossRef]
  35. Hamed, K.H. Trend detection in hydrologic data: The Mann-Kendall trend test under the scaling hypothesis. J. Hydrol. 2008, 349, 350–363. [Google Scholar] [CrossRef]
  36. Chen, B.; Xu, X.-D.; Yang, S.; Zhang, W.-D. On the origin and destination of atmospheric moisture and air mass over the Tibetan Plateau. Theor. Appl. Climatol. 2012, 110, 423–435. [Google Scholar] [CrossRef]
  37. Yu, H.; Wang, L.; Yang, R.; Yang, M.; Gao, R. Temporal and spatial variation of precipitation in the Hengduan Mountains region in China and its relationship with elevation and latitude. Atmos. Res. 2018, 213, 1–16. [Google Scholar] [CrossRef]
  38. Kim, S.B.; Lee, T.; Fukumori, I. The 1997–1999 abrupt change of the upper ocean temperature in the north central Pacific. Geophys. Res. Lett. 2004, 31, 1–4. [Google Scholar] [CrossRef]
  39. Li, R.; Min, Q.; Fu, Y. 1997/98 El Niño-induced changes in rainfall vertical structure in the east Pacific. J. Clim. 2011, 24, 6373–6391. [Google Scholar] [CrossRef]
  40. Kwon, M.H.; Jhun, J.G.; Ha, K.J. Decadal change in east Asian summer monsoon circulation in the mid-1990s. Geophys. Res. Lett. 2007, 34, 1–6. [Google Scholar] [CrossRef]
Figure 1. Partitioned regions and meteorological stations over the Qinghai–Tibet Plateau (TP).
Figure 1. Partitioned regions and meteorological stations over the Qinghai–Tibet Plateau (TP).
Water 11 01960 g001
Figure 2. Spatial distribution of rainy season onset on the TP (1971–2015). (Stations are sparse in the western TP; therefore, the five stations in the west were not part of the contour interpolation calculations.).
Figure 2. Spatial distribution of rainy season onset on the TP (1971–2015). (Stations are sparse in the western TP; therefore, the five stations in the west were not part of the contour interpolation calculations.).
Water 11 01960 g002
Figure 3. The trend in the rainy season onset on the TP (1971–2015).
Figure 3. The trend in the rainy season onset on the TP (1971–2015).
Water 11 01960 g003
Figure 4. The differences in the trend in the rainy season onset for each region (1971–2015).
Figure 4. The differences in the trend in the rainy season onset for each region (1971–2015).
Water 11 01960 g004
Figure 5. Mutation of the annual precipitation over the TP (a) moving t-test; (b) cumulative anomaly.
Figure 5. Mutation of the annual precipitation over the TP (a) moving t-test; (b) cumulative anomaly.
Water 11 01960 g005aWater 11 01960 g005b
Figure 6. The differences in the trend in the rainy season onset for each station (a) before mutation; (b) after mutation.
Figure 6. The differences in the trend in the rainy season onset for each station (a) before mutation; (b) after mutation.
Water 11 01960 g006aWater 11 01960 g006b
Table 1. List of meteorological stations in the region and the annual average rainy season onset (Julian Day).
Table 1. List of meteorological stations in the region and the annual average rainy season onset (Julian Day).
Station (Region)OnsetStation (Region)OnsetStation (Region)OnsetStation (Region)Onset
Tashikuergan (A)149Pali (C)135Changdu (F)150Ganzi (I)132
Wudaoliang (A)156Jiacha (C)150Dege (F)146Banma (I)134
Gaize (A)168Mangya (D)171Batang (F)146Seda (I)143
Anduo (A)158Lenghu (D)194Mangkang (F)150Daofu (I)139
Tuotuohe (A)162Dachaidan (D)140Milin (G)124Aba (I)133
Shiquanhe (B)194Delingha (D)130Chayu (G)76Maerkang (I)127
Bange (B)163Nuomuhong (D)166Tuole (H)146Hongyuan (I)128
Naqu (B)154Wulan (D)136Yeniugou (H)152Songpan (I)114
Shenzha (B)169Dulan (D)120Qilian (H)144Luolong (I)127
Dangxiong (B)155Chaka (D)140Gangcha (H)149Bomi (I)86
Jiali (B)144Xiaozaohuo (E)187Menyuan (H)134Basu (I)130
Pulan (C)78Geermu (E)171Wushaoling (H)147Xinlong (I)140
Lazi (C)164Zhiduo (E)157Qiabuqia (H)130Litang (I)147
Nanmulin (C)164Qumalai (E)152Xining (H)130Qianning (I)136
Rikaze (C)167Maduo (E)148Guizhou (H)135Linzhi (I)130
Nimu (C)159Qingshuihe (E)144Xinghai (H)137Zuogong (I)152
Gongga (C)163Zhongxinzhan (E)151Guinan (H)129Daocheng (I)157
Lasa (C)156Guoluo (E)136Tongde (H)138Kangding (I)118
Mozhugongka (C)149Dari (E)141Zeku (I)144Derong (I)163
Qiongjie (C)157Zaduo (F)146Tongren (I)124Deqin (I)103
Zedang (C)158Yushu (F)142Henan (I)134Jiulong (I)134
Nielamu (C)73Shiqu (F)144Jiuzhi (I)134Zhongdian (I)126
Dingri (C)177Suoxian (F)146Maqu (I)140Xiaojin (J)118
Jiangzi (C)163Biru (F)144Langmusi (I)140Muli (K)148
Langkazi (C)163Nangqian (F)143Ruoergai (I)134Gongshan (K)62
Cuona (C)111Leiwuqi (F)144Hezuo (I)130
Longzi (C)160Shiquluoxu (F)141Dingqing (I)148
Table 2. The differences in the regional trends in the rainy season onset before and after the mutation.
Table 2. The differences in the regional trends in the rainy season onset before and after the mutation.
Number of StationsRegionBefore the MutationAfter the Mutation
AdvancedUnchangedDelayedAdvancedUnchangedDelayed
5A. Northern Tibet5/5 3/5 2/5
6B. Central Tibet4/6 2/63/62/61/6
18C. Southern Tibet5/12 7/127/185/186/18
8D. Qaidam Basin5/6 1/63/6 3/6
9E. Southern Qinghai7/7 5/7 2/7
12F. Changdu4/81/83/87/113/111/11
2G. Dawang–Chayu1/1 1/2 1/2
12H. Qilian–Qinghai Lake10/11 1/116/111/114/11
31I. Bomi–Western Sichuan18/271/278/2711/263/2612/26
2J. Northern Yunnan2/2 1/2 1/2
1K. Sichuan 1/11/1
106Total61/862/8623/8648/9514/9533/95
* Due to a data shortage and site relocation, the number of stations before and after the mutation is not equal to the number of all stations.

Share and Cite

MDPI and ACS Style

Hu, Y.; Xu, J.; Huang, Y.; Zhou, Y.; Pang, Y.; Shi, Z.; Chen, X. Spatial and Temporal Variations in the Rainy Season Onset over the Qinghai–Tibet Plateau. Water 2019, 11, 1960. https://doi.org/10.3390/w11101960

AMA Style

Hu Y, Xu J, Huang Y, Zhou Y, Pang Y, Shi Z, Chen X. Spatial and Temporal Variations in the Rainy Season Onset over the Qinghai–Tibet Plateau. Water. 2019; 11(10):1960. https://doi.org/10.3390/w11101960

Chicago/Turabian Style

Hu, Yuekai, Junfeng Xu, Yuxin Huang, Yinying Zhou, Yuwen Pang, Zhou Shi, and Xiaojun Chen. 2019. "Spatial and Temporal Variations in the Rainy Season Onset over the Qinghai–Tibet Plateau" Water 11, no. 10: 1960. https://doi.org/10.3390/w11101960

APA Style

Hu, Y., Xu, J., Huang, Y., Zhou, Y., Pang, Y., Shi, Z., & Chen, X. (2019). Spatial and Temporal Variations in the Rainy Season Onset over the Qinghai–Tibet Plateau. Water, 11(10), 1960. https://doi.org/10.3390/w11101960

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

Article Metrics

Back to TopTop