Abstract

Tibetan Plateau (TP) mesoscale vortex (TPMV) was regarded as one of the most important rain bearing systems in China. Previous studies focused on the mechanisms of the TPMV in the viewpoint of deterministic forecast; however, few studies investigate the predictability of the TPMV using the Observing System Research and Predictability Experiment (THORPEX) Interactive Grand Global Ensemble (TIGGE) from the European Center for Medium Range Weather Forecasts (ECWMF). This paper investigates the location and the intensity of the larger-scale synoptic systems that influenced the development of the TPMV and its associated heavy rainfall by correlation and composite analysis. The case study on 18 July 2013 shows that stronger Balkhash Lake ridge, weaker Baikal Lake trough, and weaker western Pacific subtropical high (WPSH) are favorable to formation of TPMV over the Sichuan basin (SCB); otherwise, weaker Balkhash Lake ridge, stronger Baikal Lake trough, and stronger WPSH result in formation of TPMV to west of the SCB slightly. After the initial time, forecast for next 48 h of the geopotential height over the SCB can be viewed as a precursor of the subsequent time-averaged 90–108 h forecast of TPMV. TPMV had critical contributions to the heavy rainfall over the SCB on 18 July 2013.

1. Introduction

Tibetan Plateau (TP) mesoscale vortex (TPMV), with a typical horizontal scale of several hundred kilometers, frequently occurs over the TP in boreal summer [14]. As one of the major rain bearing systems [57], TPMV can give rise to the extreme rainfall events in East Asia, including China, Korea, and Japan [2, 3]. Based on the locations and heights in the troposphere, generally, the subsynoptic or mesoscale convective systems are named as TPMV or southwest vortex (SWV), respectively. In China, TPMV is even regarded to be the second heavy precipitation causing system after the tropical cyclones [8]. For example, the severe flood in 1998 over the Yangtze River basin was considered to be triggered by the eastward-propagating TPMV [911].

Partly due to the scarcity of meteorological observations over the TP and its adjacent area, our knowledge concerning the formation of TPMV is less deep than that of the tropical cyclones in China. To compensate for the limitation in capturing the TPMV and the SWV by a conventional observational network, two major Tibetan Plateau scientific field experiments were conducted in 1979 and 1998 [12], respectively. Over the past decades, a number of hypotheses have been made to explain the genesis, movement, and structure of the TPMV and SWV from thermodynamical and dynamical environments over the plateau [13, 14]. Takahashi [15] suggested that the lower level cold air mass contributed to the formation of TPMV. Ye and Gao [7] pointed out that the shallow cyclonic vortices developed and moved out from the TP when the favorable conditions occurred at a higher level. The large-scale circulations and the TP topographic effects and the release of convective latent heat were contributed to the development of TPMV during the rainy season (e.g., Wang [3]; Shi et al. [9]; Dell’Osso and Chen [16]; Gao [17]; Shen et al. [18]; Sugimoto and Ueno [19]; Zhou et al. [20]). Recently, the climatology of TPMVs and SWVs were studied by objective tracking approaches [21, 22]; these results provided the geographical distribution of TPMVs and SWVs and further enhanced the understanding of the TPMVs and SWVs.

As ensembles of numerical weather prediction models are more and more used in the daily weather forecasting, some investigators used ensemble forecasts to analyze the dynamics of the weather system. For example, Schumacher [23] used ensemble forecasts to determine the synoptic and mesoscale factors on the development of a vortex and its associated heavy rainfall over the southern plains of the United States during 25 to 30 June 2007. Qian et al. [24] examined the uncertainties in the track forecast of the European Center for Medium Range Weather Forecasts (ECMWF) operational ensemble about super typhoon Megi. Torn [25] also applied the techniques to investigate the dynamics of the tropical cyclones (TCs). Hawblitzel et al. [26] as well as Sippel and Zhang [27, 28] explored the genesis of mesoscale convective vortices (MCVs) and TCs using similar techniques. Recently, Li et al. [29] examined the parameters associated with the SWV movement using an ensemble prediction system to diagnose the dynamical and thermodynamical characteristics about the eastward movement and growth of SWV.

While the previous studies mostly focused on the genesis, movement, and structure of the TPMV in viewpoint of determinable prediction, few studies [30] have investigated the TPMV using the ECWMF Ensemble Prediction product [31], e.g., using the ECMWF ensemble data and NCEP reanalysis data, Wang et al. [30] found that the plateau vortex moving eastward and coupling with SWV is the key factor to the persistent heavy rainfall over the SCB on 30 June 2013. In this study, we attempt to investigate the weather systems that favor for or not favor for the development of TPMV and the associated heavy rainfall, using the Observing System Research and Predictability Experiment (THORPEX) Interactive Grand Global Ensemble (TIGGE; Bougeault et al. [32]). The following questions will be answered: why do some ensemble members produce a TPMV-associated heavy rainfall, while others do not? What are the precursors to the development or nondevelopment of TPMV in the ensemble? How long is leading time to forecast the TPMV and the associated heavy rainfall? In addition to answering these questions, another important purpose of this study is to evaluate the utility of the TIGGE in the one of the most complicated topography areas in the world.

The remainder of this article is organized as follows. Section 2 briefly introduces the datasets and methodology. In Section 3, the episode selection and the synoptic overview are presented. Section 4 presents the results, including differences in the strength of TPMV of the ensemble members, lager-scale factors influencing the development of the TPMV, and detailed processes that are favorable or detrimental to the development of TPMV and its associated heavy rainfall over the Sichuan basin (SCB). Finally, the conclusions and discussion are provided in Section 5.

2. Data and Analysis Methods

Six hourly ERA-Interim reanalysis data with a horizontal 0.75° × 0.75° resolution and 37 pressure levels from 1000 to 1 hPa [33] are used for synoptic analyses. To determine the synoptic mesoscale factors that lead to the development of the TPMV and its associated heavy rainfall, the ensemble forecast data used here are TIGGE from the ECWMF Ensemble Prediction System [32].

The ensemble includes a control run and 50 members, and the initial perturbation method of members is adopted by singular vectors in pairs and by stochastic physics. A 51-member ensemble is coming from the spectral truncation of T399 with 62 vertical levels and through 240 forecast hours from the TIGGE archive. The horizontal fields were bilinearly interpolated from global Gaussian grid to 0.75° × 0.75° to coincide with the ERA-Interim reanalysis data. The daily gauged rainfall data (1800–1200 UTC) are provided by Meteorological Information Center of China Meteorological Administration (CMA).

To elucidate the favorable or detrimental processes to the development of TPMV and heavy rainfall on 18 July 2013 (1800 UTC 17 July–1200 UTC 18 July), the value of area-averaged 500 hPa geopotential height (GPH) and original gauged rainfall over the SCB (26°N–34.25°N, 97.25°E–108.5°E) of 18 July (4 times average) are defined as the TPMV index (TPMVI) and P, respectively.

Following the method used by Sippel and Zhang [27], linear correlations between the TPMVI and 500 hPa GPH within the entire ensemble members are calculated at each grid point and different forecast time. Correlations are used here for simplicity to illustrate studying physical relationships between several atmospheric variables and TPMV within the ensemble, although the linear correlations do not necessarily interpret causality (Schumacher [23]; Sippel and Zhang [27]; Hakim and Torn [34]). To physically interpret the difference between the stronger and weaker members for the predictability of TPMV, the composite (averaged) method is also used in this paper.

3. Case Overview

Located in the east of TP, the SCB (red dashed box in Figure 1) experienced consecutive heavy rainfall events (Figure 2) in July 2013. The continuous precipitation caused the deadly and destructive flooding for Sichuan Province [35]. The total economic loss brought by the flood was viewed as the most seriously events since the heavy rainfall occurred over the SCB in 11–15 July 1981 [35]. As the heaviest rainfall occurred on 18 July 2013 (Figure 2), this episode is selected for study in this paper.

The heavy rainfall is associated with the favorable synoptic situations [36]. At 500 hPa weather charts on 17 July 2013 in Figure 3(a), two low pressure centers were observed in mid-high latitudes. One was located in the Ural Mountains, and the other was in Baikal Lake; meanwhile, another weak ridge was located over the Sea of Okhotsk. Kuo et al. [37] had proposed that the unusually strong Baikal Lake trough in summer provided a favorable environment for the organization of a mesoscale vortex over the SCB. In middle latitudes, the flow between the 40°N–50°N was zonally distributed, although the short wave troughs over the Balkhash Lake and the southwest of the Baikal Lake were still observed. Meanwhile, a northeast-southwest trough was located in the north of the SCB. Notably, western Pacific subtropical high (WPSH) meridionally oriented along a northeast-southwest axis was over the west of the Pacific Ocean. After 24 hours in Figure 3(b), on 18 July, the trough over Ural Mountains was steady and varied a little, and the Baikal Lake trough moved northward and became more meridional. And at this time, a closed midlevel circulation (TPMV) had developed, and the TPMV brought about heavy rainfall exceeding 100 mm during the period from 1200 UTC 17 July to 1200 UTC 18 July.

4. Results

4.1. Composites Analysis of the Stronger and Weaker Members

ECWMF Ensemble Prediction System [29] forecasts were initialized at 0000 UTC 14 July 2013. Figure 4 shows the reason why we chose this initialization time. As shown in Figure 4, the initialization time at 0000 UTC 14 July not only has enough lead time for studying the forecast signal but also has more accurate forecast signals of the TPMV than 0000 UTC 15 July and 0000 UTC 16 July. The TPMV is one of the important weather systems producing the extreme rainfall events in SCB [30]. To measure the differences among the ensemble results about TPMV, we divide the 50 members into stronger (weaker) groups with 6 minimum (maximum) TPMVI values. So the “stronger” group consists of 6 members, and the “weaker” group consists of 6 members.

The composites of the 500 hPa GPH fields and its associated wind fields on 18 July for these two groups are shown in Figure 5. A common feature of both stronger group and weaker group is the quasizonal alignment of low GPH centers between 50°N and 60°N. Two troughs are located over the southeast of the Ural Mountains (50°N–60°N, 60°E–80°E) and the Baikal Lake (50°N–60°N, 100°E–120°E, Baikal Lake trough hereafter), respectively. The striking difference is that the Baikal Lake trough in stronger group is weaker and tilts less than that in weaker group while the 5820 contours are much closer to the SCB in stronger group than that in weaker group. In middle latitudes (40°N–50°N, 100°E–120°E), the flow is essentially zonal over Eurasia. The southern end of the trough over the south of Baikal Lake moves to the north of the SCB, and the flow pattern forms the “North trough-South vortex,” which is a typical flow pattern for heavy rain over the SCB [38]. Baikal Lake trough is slightly further south and deeper in the stronger group than that in the weaker group. Note that the TPMV is just located over the SCB in Figure 5(a) and slightly further east in the stronger group than that in the weaker group in Figure 5(b). WPSH is more meridional and its western boundary extends further east in the stronger group than that in the weaker group.

Figure 5(c) further shows the differences between two groups. The distinction is that a significant “- + -” pattern of GPH anomalies dominates over the Eurasian continent in the meridional direction. The negative values cover the northwest of the Sea of Okhotsk (55°N–65°N, 82°E–140°E) and the low latitudes from the east of the TP to the west Pacific Ocean, with positive values over the Eurasian continent. According to an objective criterion of the trough and ridge, the negative GPH differences imply stronger troughs in the stronger group and the positive GPH differences imply stronger ridges in the stronger group. These features imply that Baikal Lake trough and the WPSH are weaker in the stronger group than that in the weaker group; meanwhile, Balkhash Lake ridge (40°N–50°N, 82°E–100°E) and TPMV in the stronger group are stronger than that in the weaker group. These features are consistent with the discussions above.

4.2. Lager-Scale Factors Influencing the Development of the TPMV

To elucidate the favorable or detrimental processes to the development of the TPMV, composites of the “stronger” group and the “weaker” group and linear correlations between 500 hPa GPHs and TPMVI within the entire ensemble members are calculated at each grid point and different forecast time. Forecasted atmospheric fields leading up to the chosen time can be regarded as precursors to the development of the TPMV over the SCB. Due to the movement of the mid-high latitude synoptic systems over time, the possible precursors of synoptic systems are masked by asterisk in the mid-high latitudes. Because of the difficulty to indicate the jump of the WPSH at 500 hPa, the area in the dashed blue box where the correlations coefficient is roughly beyond the 95% confidence level near WPSH depicts the influence of the WPSH. To measure the WPSH, the regions used to be calculated GPH in the dashed blue box from Figure 6(a) to 6(i) were (22.5°N–31.5°N, 113.5°E–125°E), (22.5°N–32.5°N, 113.5°E–123.5°E), (24.5°N–32.5°N, 112.5°E–120.5°E), (24.5°N–33°N, 112.5°E–123.5°E), (25.5°N–32°N, 110.5°E–124.5°E), (25.5°N–32.5°N, 110.5°E–119.5°E), (26.5°N–34.5°N, 110.5°E–119.5°E), and (24.5°N–33.5°N, 108.5°E–117.5°E), respectively. The TPMV is shown in the black box over the SCB in Figure 6. As shown in Figure 6, area with continuous correlations is obvious from the beginning of the forecast for next 48 h. The negative correlations are evident over Ural Mountains at the forecast for next 48 h in Figure 6(e) and then moved to the southwest of the Baikal Lake in Figure 6(f). This indicates that there is a stronger TPMV when there is a weaker Baikal Lake trough. The shallower Baikal Lake trough also suggests that weaker cold air flow coming from the mid-high latitudes benefits the formation of the TPMV. This feature has some differences from the proposal by Kuo et al. [37], which stated that colder air flows would benefit the formation and development of the TPMV. The difference of the east winds located in the north of the SCB in Figure 7 and the weaker Baikal Lake trough in the stronger members (Figures 7(f)7(i)) further support the viewpoint mentioned above.

Besides the feature of a negative correlation from the Ural Mountains to the southwest of the Baikal Lake mentioned above, another positive correlation center located in the south of the Balkhash Lake (38°N–42°N, 75°E–85°E) is also observed in Figure 6. This positive correlation center over Balkhash Lake area is apparent in the forecast for next 60 h (Figure 6(f)). As the time progressed, the positive correlation area mentioned above becomes wider in Figure 6(g), and finally it merges with another positive center over the northwest of the SCB (22°N–50°N, 75°E–100°E) in the forecast for next 84 h (Figure 6(h)). Corresponding to the positive center over the northwest of the SCB in Figure 6(f), a stronger ridge in the south of the Balkhash Lake is observed in the stronger members with respect to the weak members, where positive difference of GPH is also shown in Figure 7(f). This suggests that the stronger Balkhash Lake ridge in the south of the Balkhash Lake favors for the formation of the TPMV over the SCB. Otherwise, the weaker Balkhash Lake ridge in the south of the Balkhash Lake results in formation of TPMV to west of the SCB slightly. Figure 7 also demonstrates that the westerly air flow is weaker in the stronger members than that in the weaker members in the north of the SCB. Combining with the shallower Baikal Lake trough in the Figure 6(i), we can deduce that a weaker cold flow coming from the mid-high latitudes favors the formation and the development of the TPMV.

Meanwhile, in the low-mid latitudes, the previous GPHs over the SCB and WPSH are also important for the formation of TPMV. As shown in Figure 6, a persistent and steady positive correlation center between the TPMVI and the 500 hPa GPH near the WPSH is observed from the initial time. To further illustrate the relationship between the WPSH and TPMV, Figure 8 shows the time series of the correlations between the WPSH and TPMVI. The area-averaged GPHs at 500 hPa in the dashed blue box from Figures 6(a) to 6(i) represent WPSH. The significant positive correlation values increase since the forecast for next 24 h (valid 0000 UTC 15 July), and most correlation values exceed the horizontal line of the 95% confidence level with R = 0.36 after the forecast for next 24 h in Figure 8. This suggests that there is a stronger TPMV when WPSH is weaker. Figure 6 also indicates that TPMV is obstructed by the stronger WPSH in the weaker members, and TPMV is easily located over the TP and slightly further west than that in the stronger members. It is obvious that the TPMV located even slightly further west is not favorable to the heavy rainfall over the SCB.

Figure 9 shows that 500 hPa GPH in the forecast for the next 48 h over the SCB plays an important role on the TPMV on 18 July 2013. The correlation coefficient between the averaged GPHs and TPMVI is 0.33, which passed the 95% confidence level. Figures 7(h)7(i) also illustrate that the negative value of the GPH over SCB is corresponding to the stronger TPMV in the subsequent time-averaged 90–108 h forecast of TPMV. The GPH in the forecast for next 48 h over the SCB can be viewed as a precursor of the subsequent time-averaged 90–108 h forecast of TPMV in the east of TP.

4.3. Relationship between the TPMV and the Heavy Rainfall over the SCB

Similar instantaneous linear relationships are also evident between the averaged 500 hPa GPH and precipitation (P) over the SCB (Figure 10). The correlation coefficient between the averaged 500 hPa GPH and is −0.348, which passed the 95% confidence level. This further demonstrates that TPMV may lead to more rainfall over the SCB. For example, higher (lower) heights usually correspond to less (more) precipitation. Since the averaged 500 hPa GPH can be regarded as an indication of the strength of the TPMV, the above results imply a positive correlation between the strength of TPMV and the total precipitation there.

5. Concluding Remarks

The heavy rainfall event was previously shown to be very difficult to be predicted at the medium term forecast [35]. The heavy rainfall was mainly influenced by the TPMV over the SCB. To investigate larger-scale factors influencing the development of the TPMV and associated heavy rainfall on 18 July 2013 over the SCB, statistical relationships were calculated to examine the atmospheric conditions influencing the development of the TPMV. The above results are summarized in Figure 11, which illustrates the weather situations in favor or not for the development of the TPMV and associated heavy rainfall.

TPMV, as a mesoscale baroclinic vortex, is influenced by the interaction of the cold and warm flows significantly [7]. As shown in Figure 11, the TPMV is found to be closely related to the strength of two mid-high latitudes weather systems (e.g., Balkhash Lake ridge and Baikal Lake trough) and previous GPHs over the SCB and the WPSH area. Relatively weaker cold flow with a stronger ridge over the eastern of Balkhash Lake and a weaker trough over the southwest of Baikal Lake, meeting with the weaker warm flow coming from the WPSH, is favorable for the just location of TPMV over SCB. On the other hand, a relatively stronger cold flow coming from a weaker ridge over the eastern of Balkhash Lake and a stronger trough over the southwest of Baikal Lake, due to more cold flow into the western SCB as well as the resistance of stronger WPSH, results in the formation of TPMV over the west of SCB. The slightly further western location of the TPMV leads to less precipitation in the SCB. The GPHs in forecast for next 48 h over the SCB can be viewed as a precursor of subsequent time-averaged 90–108 h forecast of TPMV over the east of TP. These relatively small differences in the wind and weather situations early in the ensemble forecast, in conjunction with modifications of TPMV, will lead to very large spread in the resulting precipitation forecasts.

Compared with the deterministic forecast, by using correlation and composite methods, ensemble forecasts can demonstrate the mechanism of the TPMV formation more easily, thus helping to better understand the dynamics and predictability of the TPMV. For example, it is well known that cold air intrusion is an important forcing mechanism of the TPMV formation [3]. This case study shows that a weaker cold flow coming from the mid-high latitudes and a weaker WPSH are very important to the formation of the TPMV over the SCB. More TPMV cases can be studied using these methods to further verify the conclusion mentioned above, which will broaden the understanding of TPMV dynamic process and why this system can be predicted accurately.

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was jointly supported by the State Key Program of the National Natural Science of China (91337215 and 41475051), the Starting Foundation of the Civil Aviation University of China (2016QD05X), the Special Fund for Climate Change (CCSF201706), and the Special Fund for Development of Weather Forecasting Key Technologies (YBGJXM(2017)03-13).