Impact of the Westerly Jet on Rainfall/Runoff in the Source Region of the Yangtze River during the Flood Season

School of Atmospheric Sciences, Plateau Atmosphere and Environment Key Laboratory of Sichuan Province, Joint Laboratory of Climate and Environment Change, Chengdu University of Information Technology, Chengdu 610225, China Key Laboratory for Land Surface Process and Climate Change in Cold and Arid Regions, Chinese Academy of Sciences, Lanzhou 730000, China Chengdu Institute of Plateau Meteorology, China Meteorology Administration, Heavy Rain and Drought-Flood Disasters in Plateau and Basin Key Laboratory of Sichuan Province, Chengdu 610072, China Bazhong Meteorological Office of Sichuan Province, Bazhong 636000, China


Introduction
e Yangtze River originates at the terminuses of the glaciers on Geladaindong Mountain in the middle section of the Tanggula Mountain Range in Qinghai Province. e headwaters consist mainly of three rivers, namely, the Dangqu River (southern source), the Tuotuo River (main source), and the Chumar River (northern source). e Tuotuo River converges with the Dangqu River to form the Tongtian River, which flows southeastward and converges with the Batang River near Yushu County in Qinghai Province to form the Jinsha River. Hydrologically, the region upstream of the Zhimenda hydrological station is generally treated as the source region of the Yangtze River (SRYR). Changes in runoff in the SRYR not only affect the regional ecological environment but also significantly affect the water resources in the middle and lower reaches of the Yangtze River. erefore, a number of researchers have studied the characteristics of runoff in the SRYR [1][2][3][4][5]. e runoff in the SRYR is mainly concentrated between May and September. From the 1960s to the 1990s, the runoff in the SRYR exhibited an overall decreasing trend with four main stages, namely, a wet stage, followed by a dry stage, another wet stage, and another dry stage. e 1990s was a relatively particularly dry stage [6][7][8]. e SRYR entered a wet stage in approximately 2005 [9]. However, these studies only examined relatively short periods of time. By reconstructing the runoff of the Tongtian River since 1485, Qin et al. [10] analyzed the long-term characteristics of its changes. Many factors could influence runoff in the SRYR. For example, a precipitation increase could cause the runoff in the SRYR to increase, whereas an evaporation increase could cause the runoff in the SRYR to decrease. e heating of the Tibetan Plateau (TP), enhanced plateau monsoon, and enhanced glacier melting could also increase the runoff in the SRYR [2,5,[11][12][13]. However, precipitation is the principal factor that affects the runoff [14,15]. e relationship between precipitation and periodic fluctuations in runoff in the SRYR is much more significant than that between the temperature and runoff; in addition, the annual runoff also corresponds to annual precipitation in terms of interdecadal changes and years when abrupt changes take place [16,17]. us, precipitation is the principal factor that affects the runoff in the SRYR. e upper-level westerly jet (WJ) is a planetary-scale atmospheric circulation system in the upper troposphere and is an important circulation system that affects the weather and climate in East Asia [18][19][20][21][22][23]. Hence, numerous researchers have studied the WJ [24][25][26][27][28][29][30]. Changes in the position and intensity of the WJ significantly affect the weather and climate in China. With regard to the influence of changes in the position of the WJ, most research has emphasized the influence of meridional movement. On the one hand, the meridional movement of the WJ can influence the Meiyu occurrence. For example, Tao et al. [31] found that the start and end of East Asian Meiyu are related to two northward jumps of the southern branch of the WJ over Asia between June and July. e second northward jump of the southern branch of the WJ occurs before Meiyu in the Jianghuai region of China [32]. On the other hand, the southward and northward jumps of the subtropical WJ in summer also affect the location of the rainband. In June, precipitation anomalies mainly occur over the sea surface south of Japan. In July, when the WJ jumps northward, the rainband moves northward from the middle and lower reaches of the Yangtze River to the Huai River and the region north of it [33,34]. Additionally, the meridional movement of the subtropical WJ can significantly influence summer precipitation in China. In summer, the subtropical WJ moves southward, resulting in increased precipitation in the Yangtze-Huai River basin and the Jiangnan region [18,20,23,26,35] and decreased precipitation in Northeast China [36,37]. e zonal movement of the subtropical WJ can also influence Meiyu. Zhang et al. noted that the eastward and westward movements of the subtropical WJ affect the start and end of the Meiyu period [38]. In addition, changes in the shape of the subtropical WJ in the eastern and western directions also affect the spatial distribution of precipitation in the middle and lower reaches of the Yangtze River during the Meiyu period [39]. If the center of the East Asian subtropical WJ is located relatively westward in January and its axis is situated relatively southward between April and May, the year may be a rich Meiyu year [40]. Furthermore, many researchers have studied the influence of changes in the strength of the subtropical WJ. Liang and Liu [41] found that when the southern branch of the WJ in winter is relatively strong (weak), there will be relatively greater precipitation in South China in the following spring and in North China in the following summer. Li and Zhang pointed out that a strong subtropical WJ can increase the precipitation in the Meiyu period [42]. Wang and Zuo further pointed out that a strong subtropical WJ decreases the precipitation in the Yangtze River basin in midsummer (16 July-14 August) [19]. Evidently, the subtropical WJ has a significant impact on precipitation in China, especially East China. e SRYR is located in the region impacted by the westerlies, but relatively few studies have been conducted to examine the influence of the WJ on the precipitation in the SRYR. In view of this, this study aims to investigate the impact of the WJ on rainfall in the SRYR and analyze its impact on runoff in the SRYR, with the goal of providing a certain basis for protection of the ecological environment of the SRYR.
is paper is organized as follows. Brief descriptions of the data and method are provided in Section 2. Characteristics of changes in runoff in the SRYR are analyzed in Section 3. Correlation analyses of runoff and rainfall in the SRYR during the flood season are provided in Section 4. e impact of the WJ on rainfall in the SRYR is analyzed in Section 5. e impact of the intensity of the WJ over the key region on atmospheric circulation patterns is investigated in Section 6. e impact of the intensity of the WJ over the key region on water vapor transport is investigated in Section 7. Section 8 discusses the major conclusions of this study.

Data.
e hydrological data used in this study were monthly runoff observation data collected from 1956 to 2012 at the Zhimenda station. In addition, the runoff at the Zhimenda station ( Figure 1) was treated as the runoff in the SRYR. e rainfall data used in this study were complete daily observation data collected at 109 stations ( Figure 1) in China between 85°E and 105°E and between 25°N and 40°N from 1961 to 2012 and were provided by the National Meteorological Information Center of the China Meteorological Administration. e daily reanalysis data, including surface pressure, zonal wind (u), meridional wind (v), vertical velocity (w), and specific humidity (q), are derived from the National Centers of Environmental Prediction/ National Centers of Atmospheric Research (NCEP/NCAR) reanalysis data, with a horizontal resolution of 2.5°. Fang et al. [43,44] noted that the NCEP/NCAR data for East Asia for the period 1958-1967 have substantial quality issues. erefore, the period 1971-2012 was selected for analysis in this study.

Methods.
e vertically integrated water vapor flux vector can be expressed as

Advances in Meteorology
Zonal component: Meridional component: In the above equations, g is the gravity acceleration, with a value of 9.80665 m·s −2 in this study; Ps is the surface pressure; Pt is the top level pressure (300 hPa in this study); q is the specific humidity; and u and v are the zonal and meridional wind, respectively. e Pearson correlation coefficient [45] is one of the most commonly used statistical tools to measure the degree of linear correlation between two given variables. e correlation coefficient between two variables-x 1 , x 2 , . . ., x n and y 1 , y 2 , . . ., y n -can be calculated as follows: (4) e overbars represent the climatological mean. In this correlation calculation, the runoff in the flood season is the sum of monthly runoff from July to September, the rainfall is the sum of daily rainfall from 16 June to 15 September, and the circulation is the average of the daily data from 16 June to 15 September. e correlation coefficient value is located in the range from −1 to 1. A positive correlation with a value greater than 0 means that the two variables change in a similar direction (i.e., both variables increase or decrease). In addition, a negative correlation with a value less than 0 implies that the two variables exhibit opposite changes (i.e., as one variable increases, the other decreases).
e station data were interpolated to the grid data by the Cressman interpolation technique [46]. Multiple passes are made through the grid with increasingly smaller radii of influence. At each pass, a new value is calculated for each grid point based on a correction factor that is determined by looking at each station within the radius of influence. For each such station, an error is defined as the difference between the station value and a value arrived by interpolation from the grid to that station. A distance-weighted formula is then applied to all such errors within the radius of influence of the grid point to arrive at a correction value for that grid point. e correction factors are applied to all grid points before the next pass is made. Observations nearest to the grid point carry the most weight. As the distance increases, the observations carry less weight. e Cressman function in Ingrid calculates the weights as follows: where W is the influence radius and r is the distance between the station and the grid point. Figure 2 shows the monthly variation in the runoff at the Zhimenda station. As shown in Figure 2, relatively high runoff occurs between June and October. In particular, the highest runoff occurs between July and September. e total runoff in these three months accounts for 59.6% of the total annual runoff. In contrast, the runoff at the Zhimenda station is notably lower between January and March and in December. Changes in runoff at the Zhimenda station during the flood season (i.e., from July and September) directly affect the annual changes, and the changes also have a significant impact on the regional ecological environment. Hence, this study focuses mainly on analyzing the runoff in the SRYR between July and September. Figure 3 shows the time series of the runoff at Zhimenda station in the flood season during 1971-2012. Overall, the runoff at the Zhimenda station during the flood season displays a slight increasing, albeit insignificant, trend. e runoff at the Zhimenda station during the flood season increased in the early 1970s, decreased in the late 1970s, and was generally high between the 1980s and the mid-1990s. e runoff significantly increased between the mid-1990s (approximately 1994) and 2012 and has been high since 2005, reaching its highest level in 2009. Evidently, there are notable trends of interannual and interdecadal changes in runoff in the SRYR.

Correlation Analysis of Runoff and Rainfall in the SRYR during the Flood Season
Some research shows that climate change in the eastern TP has a significant impact on water resources in the upper reaches of the Yangtze River and that the annual runoff in the upper reaches of the Yangtze River increases as the precipitation in the eastern TP increases [47]. us, the rainfall has a significant impact on runoff. erefore, this section focuses mainly on examining the impact of rainfall on runoff in the SRYR during the flood season. e runoff is not only affected by simultaneous rainfall but is also affected by antecedent rainfall [48]. us, a lag correlation between rainfall and runoff with the daily rainfall leading runoff by 15 days (for June 16 to September 15) is calculated (Figure 4). From Figure 4, it is seen that there is a significant positive correlation over central and southern Qinghai Province, including from the city of Golmud eastward to the large region to the west of the Golog and Hainan Tibetan Autonomous Prefectures, as well as northwestern Sichuan Province. e correlation regions with correlation coefficients greater than 0.6 are mainly distributed in the coastal zone of the SRYR upstream of the Zhimenda station and the two centers over the Tuotuo River region and Qumarleb. Note that lag correlations between rainfall and runoff with other leading times (e.g., 30, 25, 20, 10, and 5 days) were also calculated (not shown). e results show that the significant positive correlation region is highly similar to Figure 4, but the regional ranges of correlation coefficients greater than 0.7 are obviously smaller than that in Figure 4. e above evidence suggests that the rainfall leading runoff by 15 days is reasonable. erefore, more attention should be paid to the 15-day antecedent rainfall when examining the impact of rainfall on runoff in the SRYR during the flood season. In the following sections of the analysis of the impact of atmospheric circulation patterns on rainfall in the SRYR, the data for the period between June 16 and September 15 (a 15-day antecedent period of runoff) are examined.

Impact of the WJ on Rainfall in the SRYR
e aforementioned analysis shows a significant correlation between changes in runoff in the SRYR and the 15-day antecedent rainfall in the coastal zone of the SRYR. e SRYR is situated in the WJ-affected region. Does the WJ have an impact on rainfall in the coastal zone of the SRYR?
is section focuses mainly on analyzing the impact of the WJ on rainfall in the SRYR.
To analyze the relationship between the WJ and rainfall in the coastal zone of the SRYR, five stations (i.e., Wudaoliang, Tuotuohe, Qumalai, Yushu, and Qingshuihe) were selected (the black circular points in Figure 4). We define a rainfall index (I r for short) as the mean flood season rainfall of the five stations, and the formula is shown as follows: where n is the number of stations and x is the flood season rainfall, which is defined as the sum of the rainfall from 16 June to 15 September.  Siberian plains, which result in increased (decreased) rainfall in the coastal zone of the SRYR. erefore, the upper-level WJ has a close relationship with rainfall in the coastal zone of the SRYR. e northward and southward jumps of the upper-level WJ have a significant impact on rainfall [23,34,35]. Lin [49] noted that the jumps of an upper-level jet are related to the intensity of the westerly winds north of the jet. erefore, a WJ intensity index (I W ) is defined as the area-mean 200 hPa zonal winds over 45-55°N and 55-85°E. In addition, the jet position index (I L ) is defined as the latitude averaged within 55-85°E, and the latitude is the location of the maximum 200 hPa zonal winds every 2.5°of longitude. A small (large) I L indicates that the jet is at a relatively low (high) position. e correlation coefficient between I W and I L reached 0.904, which was significant at the 99% confidence level. Note that the results obtained are very similar for the reasonable change in the regions used to define I W and I L . Clearly, the magnitude of I W can satisfactorily reflect the southward and northward movements of the WJ. In other words, when the WJ over the key region is relatively strong (weak), the WJ is located relatively northward (southward). us, the impact of I W on rainfall is mainly analyzed in the following.
Will changes in the intensity of the WJ over the selected key region affect rainfall anomalies in the coastal zone of the SRYR? Figure 6 displays the normalized time series of I W and I r . e temporal correlation coefficient between them is as high as −0.6 and is significant at the 99% confidence level. Clearly, when I W is relatively low, I r is relatively high, indicating relatively heavy rainfall, particularly in 1981, 1989 and 2009. In contrast, when I W is relatively high, I r is relatively low, indicating abnormally light rainfall, particularly in 1973, 1978, 1984, 1990, and 1997. e simultaneous correlation coefficients between I W and rainfall were further calculated (Figure 7(a)). From Figure 7(a), the significant negative correlation is distributed along the coast of the upper reaches of the Yangtze River, with a center (with a correlation coefficient less than −0.5) over the region near the Tuotuo River. In addition, high-I W (low-I W ) years are defined as years when the normalized I W is greater (less) than 1 (−1). Based on these criteria, six years, including 1973, 1978, 1984, 1990, 1994, and 1997, are identified as high-I W years, and another six years, i.e., 1981,1982,1989,2003,2009, and 2010, are identified as low-I W years. Figure 7(b) shows the composite anomalies of rainfall between the high-I W and low-I W years. Pronounced negative rainfall anomalies are observed over the coastal zone in the upper reaches of the Yangtze River, with an amplitude of approximately 90 mm. e results are similar to those in Figure 7(a). e above evidence suggests that when abovenormal (below-normal) zonal winds occur over the key region, the rainfall in the coastal zone of the SRYR decreases (increases). us, changes in the intensity of the WJ over the key region indeed have an impact on rainfall in the coastal zone of the SRYR.

Impact of the Intensity of the WJ over the Key
Region on Atmospheric Circulation Patterns e analysis in the previous section shows that changes in the intensity of the WJ over the key region have a significant impact on rainfall anomalies in the coastal zone of the SRYR. What is the mechanism of impact?
is section focuses mainly on analyzing the impact of changes in the intensity of the WJ over the key region on the zonal and meridional vertical circulation patterns. Figure 8(a) shows the simultaneous correlation between I W and the zonal vertical circulation averaged over 30-37.5°N. As shown in Figure 8(a), the SRYR is dominated by descending air motion. Additionally, there are significant easterly winds from the surface to 100 hPa over the region to the west of the SRYR. However, no significant correlation is observed over the region to the  Advances in Meteorology east of the SRYR. In addition, Figure 8(b) shows composite anomalies of the zonal vertical circulation averaged over 30-37.5°N between the high-I W and low-I W years. From Figure 8(b), pronounced descending motion anomalies over the SRYR and significant easterly wind anomalies over the region to the west of the SRYR are observed. e results are similar to those in Figure 8(a). e above evidence implies that when the I W is greater (less) than the normal value, the WJ position moves northward (southward), which results in significant easterly (westerly) wind anomalies over the region to the west of the SRYR and significant descending (ascending) motion anomalies over the SRYR. Figure 9(a) shows the simultaneous correlation between I W and the meridional vertical circulation averaged over 90-100°E. As shown in Figure 9(a), there is significant meridional circulation over the region between 25°N and 35°N, and its descending branch is situated over the SRYR, which would result in a significant descending motion of northerly winds over the SRYR.

Impact of the Intensity of the WJ over the Key Region on the Meridional Vertical Circulation.
ere is another enormous meridional circulation over the equatorial region, and its ascending branch coincides with the ascending branch of the meridional circulation over the region between 25°N and 35°N, resulting in significant ascending motion in the air flow over the region between 5°N and 27°N. Consequently, the majority of the southerly winds in the lower meridional circulation over the equatorial region transform to ascending air currents when advancing northward and thus can no longer continue to move northward; only a small portion of these southerly winds continue to move northward. However, because of the presence of a small meridional circulation in the lower atmosphere between 20°N and 25°N, these northward-moving southerly winds transform into descending air currents and thus can no longer move northward when reaching approximately 25°N. is circulation pattern is unfavorable to continued northward movement of the southerly wind toward the SRYR. In addition, Figure 9(b) shows composite anomalies of the meridional vertical circulation averaged over 90-100°E between the high-I W and low-I W years, and     Advances in Meteorology the results are similar to those in Figure 9(a). us, when the I W is greater (less) than the normal value, the southerly wind over the equatorial region cannot (can) easily advance northward, and significant descending (ascending) motion anomalies occur over the SRYR. Figure 10(a) shows the correlation between the I W and the 200 hPa winds and divergence field. From Figure 10(a), there is significant cyclonic circulation over the TP, with significant convergence over the TP. Meanwhile, in the lower level (i.e., 600 hPa, Figure 10 200 hPa zonal wind anomalies over the region from 30°N to 40°N, and positive (negative) anomalies over the area south of the sources of three rivers (i.e., Lancang River, Yellow River and Yangtze River), which result in upper-level cyclonic (anticyclonic) shear anomalies over the TP. In addition, for the 600 hPa zonal winds, positive (negative) anomalies over the area north of the sources of three rivers and negative (positive) anomalies over area south of the sources of the three rivers result in low-level anticyclonic (cyclonic) shear anomalies over the TP. Hence, the convergence (divergence) anomalies in the upper level and the divergence (convergence) anomalies in the lower level are induced by the cyclonic (anticyclonic) circulation and anticyclonic (cyclonic) circulation anomalies, respectively. ese processes result in descending (ascending) motion anomalies over the SRYR, which decrease (increase) rainfall. Hence, these circulation patterns also impact runoff. Because runoff has a close relationship with rainfall, circulation anomalies indirectly impact runoff by impacting rainfall. When the antecedent rainfall decreases (increases) in the SRYR during the flood season, runoff also decreases (increases).

Impact of the Intensity of the WJ over the Key Region on Water Vapor Transport
Water vapor is the source of rainfall. Water vapor is transported to a rainfall region by large-scale air movement and, when facilitated by a specific circulation pattern, ascends and cools to form clouds, which in turn generate rainfall. e SRYR is situated in the eastern TP. Water vapor transport is particularly vital to rainfall in the SRYR. e analysis in the previous section shows that changes in the intensity of the WJ over the key region affect rainfall in the coastal zone of the SRYR by affecting the atmospheric circulations. How do the changes in this circulation pattern affect water vapor? is section focuses mainly on analyzing the impact of the intensity of the WJ over the key region on water vapor transport. Figure 11(a) shows the simultaneous correlation between I W and vertically integrated water vapor flux. From Figure 11(a), there is a significant anticyclonic water vapor transport over the region extending from the Caspian Sea eastward to the western TP, and the southerly water vapor transport on its western side is unfavorable to the southward transport of the water vapor from the Caspian Sea region. e majority of the westerly water vapor is transported eastward and transforms into southwesterly water vapor on the northern side of the anticyclonic water vapor transport, which is unfavorable to the southward transport of water vapor originating from high latitudes. e rest of the westerly water vapor transforms into northwesterly water vapor, which continues to be transported eastward. During the eastward transport process, some of the northwesterly water vapor is transported to Northeast China along the northern edge of the TP via the Hetao region. e remainder transforms into northeasterly water vapor and converges over the southern edge of the TP with the easterly water vapor transported from Taiwan Province to the south of the

Advances in Meteorology
TP via South China and then is transported westward, which obviously is unfavorable for the westerly water vapor to be transported to the SRYR. ere is also another significant westerly water vapor transport over the region extending from the low-latitude Indian Ocean to the region west of the Indochina peninsula, which is also unfavorable to the transport of water vapor from the low-latitude ocean to the SRYR. As a result, there is a decrease in the water vapor content over the SRYR. Figure 11(b) shows composite anomalies of the water vapor fluxes between the high-I W and low-I W years. As shown in Figure 11(b), the significant anticyclonic water transport anomalies over the region  Advances in Meteorology extending from the Caspian Sea eastward to the western TP are unfavorable to the southward transport of water vapor originating from high-latitude regions. e significant westerly water vapor transport anomalies from the Arabian Sea via the Indian subcontinent are unfavorable to the northward transport of water vapor originating from the low-latitude ocean.
e results are similar to those in Figure 11(a). us, changes in the intensity of the WJ over the key region have a significant impact on water vapor transport. When there are positive (negative) I W anomalies, the water vapor content of the SRYR decreases (increases). e above evidence indicates that when the I W is greater (less) than the normal value, the WJ shifts northward (southward). In response, the easterly (westerly) wind anomalies over the region to the west of the SRYR become unfavorable (favorable) to the transport of water vapor from high-latitude regions to the SRYR, and the westerly (easterly) water vapor transport anomalies over the Arabian Sea also become unfavorable (favorable) to the northward transport of water vapor originating from the low-latitude ocean.
us, the water vapor content over the SRYR decreases (increases). In addition, the anomalies descending (ascending) motion occur over the SRYR. ese factors decrease (increase) the rainfall in the coastal zone of the SRYR. Additionally, these factors indirectly impact runoff by impacting rainfall; for example, a decrease (increase) in antecedent rainfall decreases (increases) the runoff in the SRYR during the flood season.

Conclusions and Discussion
Using runoff data from the Zhimenda station, this study analyzes the characteristics of changes in runoff in the SRYR during the flood season (from July to September) and its relationship with antecedent rainfall in the SRYR and further analyzes the impact of the WJ on rainfall in the coastal zone of the SRYR. e main conclusions derived from this study are summarized as follows.
(1) Changes in runoff in the SRYR during the flood season display notable interdecadal characteristics. e runoff in the SRYR during the flood season increased in the early 1970s, decreased in the late 1970s, slowly decreased between the 1980s and the mid-1990s, and significantly increased between the mid-1990s (approximately 1994) and 2012.
(2) Fifteen-day antecedent rainfall in the coastal zone of the SRYR (i.e., June 16 to September 15) is the principal factor that affects the runoff in the SRYR during the flood season. (3) ere is a close relationship between rainfall in the coastal zone of the SRYR and the intensity of the WJ over the key region (55-85°E, 45-55°N). When the WJ over the key region is relatively strong (weak), rainfall in the coastal zone of the SRYR decreases (increases). (4) When there are positive (negative) anomalies of the I W , the jet position moves northward (southward) and the easterly (westerly) wind anomalies over the region to the west of the SRYR become unfavorable (favorable) to the transport of water vapor from high-latitude regions to the SRYR. In addition, the southerly wind over the equatorial region cannot (can) easily advance northward, which results in a westerly (easterly) water vapor transport anomalies over the Arabian Sea, which is unfavorable (favorable) to the northward transport of water vapor from the low-latitude ocean. Hence, these conditions decrease (increase) the water vapor content over the SRYR. (5) When the I W is greater (less) than the normal value, convergence (divergence) anomalies in the upper level (i.e., 200 hPa) and divergence (convergence) anomalies in the lower level (i.e., 600 hPa) are induced by the cyclonic (anticyclonic) and anticyclonic (cyclonic) circulation anomalies, respectively. ese circulation patterns result in descending (ascending) motion anomalies over the SRYR. In addition, the water vapor content over the SRYR decreases (increases). ese conditions decrease (increase) rainfall, resulting in less (more) runoff in the SRYR during the flood season.
is study analyzes the impact of the WJ on rainfall in the SRYR from only a statistical perspective. e results require further validation through numerical simulation. In addition, there are a myriad of factors that affect the rainfall in the SRYR, and this study only focuses on analyzing the impact of the WJ. Other factors, such as the plateau monsoon and South Asian High, also, to a certain extent, affect rainfall in the SRYR. erefore, changes in the intensity of the WJ do not coincide with changes in rainfall in the SRYR in some years with abnormally heavy or light rainfall. Hence, the interactions between various factors will be fully taken into consideration in future work.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare no conflicts of interest.