Abstract
Despite the importance of precipitation and moisture transport over the Tibetan Plateau for glacier mass balance, river runoff and local ecology, changes in these quantities remain highly uncertain and poorly understood. Here we use observational data and model simulations to explore the close relationship between summer rainfall variability over the southwestern Tibetan Plateau (SWTP) and that over central-eastern India (CEI), which exists despite the separation of these two regions by the Himalayas. We show that this relationship is maintained primarily by ‘up-and-over’ moisture transport, in which hydrometeors and moisture are lifted by convective storms over CEI and the Himalayan foothills and then swept over the SWTP by the mid-tropospheric circulation, rather than by upslope flow over the Himalayas. Sensitivity simulations confirm the importance of up-and-over transport at event scales, and an objective storm classification indicates that this pathway accounts for approximately half of total summer rainfall over the SWTP.
Similar content being viewed by others
Introduction
The Tibetan Plateau (TP), sometimes referred to as ‘the world’s water tower’1, is a critical region for water resources in Asia. Nearly one sixth of the global population depends on rivers with headwaters on the TP, and changes in precipitation over the TP have great impacts on glaciers, river discharge, plant phenology and natural hazards1,2,3,4,5. Summer (June–September) precipitation, which occurs mainly in convective storms (CSs), accounts for 60–90% of total annual precipitation over the TP and surrounding areas6,7,8. Trends and variability in summertime precipitation over the TP are therefore likely to exert strong influences on climate and socio-economic conditions throughout many parts of Asia.
The spatiotemporal variation of summer precipitation over the TP is complex, and is mainly regulated by a combination of summer monsoons and regional moisture recycling2,3,9,10. Rainfall over the eastern TP has been extensively studied and well-documented, with a relative abundance of ground measurements, and is closely linked to the South Asian monsoon via a moisture corridor that extends along the Brahmaputra River from the Bay of Bengal9,11,12. However, despite recent sharp declines in glacier extent2,3, the southwestern TP (SWTP) has received considerably less attention. The SWTP is a semiarid region with a vulnerable ecosystem, and precipitation systems over the SWTP often serve as precursors of synoptic systems downstream4,6,11. Some convective instabilities that develop over the SWTP cause extreme rainfall and severe flooding in East Asia as they develop and move eastward9,13.
The factors that control precipitation variations over the SWTP remain poorly understood. Here we establish that summer rainfall over the SWTP and CEI is closely related via a previously unreported ‘up-and-over’ moisture transport path, in which moist, cloud-laden air lifted by convective systems over CEI is swept over the SWTP by the mid-tropospheric circulation. This relationship is evident in both observational data sets and global climate model simulations, and depends on the coincident occurrence of convective systems over CEI and southwesterlies in the middle and upper troposphere over the Himalayas. Several lines of evidence confirm both the existence of the up-and-over moisture transport path and its importance in the water budget of the SWTP, including satellite measurements, isotopic ratios in precipitation, numerical sensitivity simulations and backward trajectory calculations. The results consistently indicate that up-and-over transport accounts for at least half of total summer rainfall over the SWTP, and show that this pathway dominates moisture supply in this region. The SWTP will be substantially drier without the up-and-over transport associated with convective systems over CEI.
Results
Linkage of summer rainfall between SWTP and CEI
Evidence of a strong connection between rainfall over the SWTP and rainfall over the Indian subcontinent emerges from historical records of precipitation (Fig. 1; see Methods for details). In particular, despite a large difference in rainfall magnitude across the Himalayas (Fig. 1a), significant decreasing trends have been observed over both the SWTP and central-eastern India (CEI; Fig. 1b). Decreases in rainfall over the Indian subcontinent over the past six decades have been attributed both to natural variability7,14 and to anthropogenic forcing15,16,17. Summer rainfall variability over the SWTP is highly correlated with summer rainfall variability over CEI (Fig. 1c). Indeed, summer precipitation over CEI is a reliable predictor for summer precipitation over the SWTP (Fig. 2a), with correlations of r=0.68 for APHRODITE18 data (based on rain gauge measurements during 1951–2007; P<0.001) and r=0.64 for TRMM19 data (during 1998–2013; P<0.01). The correlation between TRMM and APHRODITE estimates of precipitation exceeds 0.85 during their period of overlap (1998–2007; P<0.01) with a difference of <8% between their mean magnitudes. Strong correlations between precipitation variability over these two regions are likewise evident in two other observational data sets, as well as in simulations of historical climate and scenario-based projections of future climate by many global climate models (Supplementary Fig. 1). These results indicate a close and robust relationship between summer rainfall over the SWTP and summer rainfall over CEI that transcends data sources and model configurations. Quantitative differences in the correlation coefficients may be tied to differences in model physics, dynamics or resolutions, and are not explored here.
The strong relationship between summer rainfall over the SWTP and summer rainfall over CEI is somewhat surprising, as these two regions are separated by the Himalayan mountain range. The Himalayas are on average more than 5,000 m above mean sea level between 75 and 95° E, and act as a natural barrier for northward water vapour transport into the SWTP. Northward flow impinging on the Himalayas is mechanically lifted by the orography, with the southern foothills of the Himalayas among the most prominent regions of orographic precipitation on Earth20,21. This situation creates a strong precipitation gradient in the lee of the Himalayas. The sharp decrease in precipitation on the leeside of the mountain range suggests that the Himalayas strictly limit upslope moisture transport under stably stratified conditions. In other words, these mountains are sufficiently high that condensation during upslope flow ‘wrings out’ most of the water, leaving little moisture available for precipitation in the interior of the SWTP. However, the interior of the SWTP experiences more precipitation than other similar places blocked by high mountain ranges, such as desert areas in the southwestern United States and South America. Furthermore, numerical modelling experiments9 and analyses of the isotopic composition of rainfall22,23,24,25,26,27 suggest that much of the rainfall over the SWTP during boreal summer can be traced back to the nearby ocean. The strong connection between the two regions must therefore be associated with alternative routes for moisture transport, rather than direct upslope flow. These results call for a re-evaluation of moisture transport into the SWTP and a careful analysis of how the relationship between the SWTP and CEI is established.
An up-and-over route for moisture transport into SWTP
Summertime rainfall over the Indian subcontinent is dominated by deep CSs10,21,28,29,30. These CSs originate over the surrounding oceans, and are then observed to move northward or northwestward in company with strong onshore moisture transport31. Individual CSs may organize and develop into monsoon depressions, which in turn contribute significantly to monsoon rainfall over the northern part of the Indian subcontinent32,33. CS development preferentially occurs when monsoon southwesterlies and high-latitude northwesterlies converge and ascend along the southwestern flank of the Himalayas13.
Almost every large precipitation event over the SWTP occurs in tandem with a corresponding event over CEI. For example, precipitation over these two regions during summer 2002 (Fig. 2b) exhibits similar variations at both daily and weekly time scales (r=0.72, P<0.01). Variations at these time scales correspond to episodes of strong rainfall; for example, between June and September 2002, the same 34 days accounted for >60% of the total rainfall over both regions. The co-variability of precipitation over the two regions is thus primarily determined by the correspondence of individual convective precipitation events. We propose that this link manifests as strong transport of moisture and hydrometeors at middle-to-upper levels of the troposphere: moist, cloud-laden air is elevated by deep CSs over the Indian subcontinent and then advected over the Himalayas into the SWTP by the mid-tropospheric circulation. We term this transport pathway the ‘up-and-over’ route.
Intrusive and non-intrusive deep CSs
Not all precipitation events over CEI are accompanied by strong precipitation over the SWTP. In some cases, convective precipitation over CEI is confined to the south side of the Himalayas. We therefore categorize CSs over CEI into two types based on the characteristics of the mid-tropospheric flow near the Himalayas (see Methods). We define each CS as intrusive if the mid-tropospheric circulation permits it to influence the SWTP, and non-intrusive if its effects are confined to the south side of the Himalayas. Supplementary Figure 2 shows examples of typical intrusive and non-intrusive CS events. Both examples were characterized by a closed low-level (850 hPa) centre of low pressure south of the Himalayas; however, strong rainfall extended northward into the SWTP for the intrusive CS and did not for the non-intrusive CS. The northward extension of rainfall in the intrusive case was associated with strong cross-barrier (southwesterly) winds of up to 10 m s–1 at 500 hPa and a region of persistent large-scale ascent that extended from the Indian subcontinent to the SWTP. At 300 hPa, strong southwesterlies ahead of an approaching upper level trough pushed the upper tropospheric high southeastward. These middle and upper tropospheric circulation patterns promoted the advection of convectively lifted moisture and hydrometeors from the Indian subcontinent northeastward over the SWTP, in stark contrast to the circulation patterns observed during the non-intrusive CS. In the non-intrusive case, the low-pressure system extended upward from the surface to 500 hPa, with predominant easterly flow along the mountain range. Meanwhile, the centre of high pressure at 300 hPa straddled the mountain range, with northwesterly flow towards the SWTP and northeasterly flow towards the Indian Subcontinent. The absence of cross-barrier flow in the middle and upper troposphere (500 and 300 hPa) effectively prevented the northward transport of moisture and hydrometeors, so that rainfall was mostly confined to the Indian subcontinent. The moist static energy at 300 hPa was likewise much larger during the intrusive case than during the non-intrusive case, indicating that upper levels of the troposphere experienced greater moistening and diabatic heating during the intrusive case.
WRF simulations and backward trajectory calculations
It is difficult to distinguish the relative contributions of up-and-over transport, upslope transport and local processes using observations alone. For example, the pattern of precipitation we attribute to up-and-over transport might instead be generated by two separate rainfall events over the SWTP and CEI. We therefore conduct a series of sensitivity simulations using the Weather Research and Forecasting (WRF) model (see Methods), in which we estimate the quantitative contributions of up-and-over transport, upslope transport and local recycling processes to the precipitation over the SWTP by reducing the advection of moisture and hydrometeors across the southern boundary into the SWTP within certain layers (Supplementary Fig. 3). A sensitivity simulation with hydrometeor advection entirely disallowed gives negligible precipitation change over the SWTP as compared with the control simulation (Figure not shown); however, this does not necessarily indicate that hydrometeor transport is unimportant given the large uncertainties in model simulated hydrometeor fields. In light of these uncertainties, we focus on the moisture transport within different layers for the following sensitivity simulations. We assume that water vapour advection across the southern boundary of the Plateau over CEI initially located more than 5 km above ground level (AGL) is mainly related to up-and-over transport, while advection of CEI water vapour initially located below 2.5 km AGL is mostly associated with upslope transport. Advection of CEI water vapour between 2.5 and 5 km AGL is assumed to be a mixture of the up-and-over and upslope transport. We use transport in this layer to represent uncertainties associated with partitioning between the two pathways.
The distribution of precipitation and the upper level circulation produced by the control simulation agree well with TRMM observations and ERA-Interim reanalyses (Fig. 3). Simulated precipitation over the SWTP decreases by ∼56% when advection of upper layer (>5 km AGL) water vapour is reduced by half. Halving the transport of water vapour at levels >2.5 km, AGL further reduces the precipitation, with simulated rainfall only ∼18% of that produced in the control simulation. Simulated precipitation over the SWTP essentially disappears (∼2% of the control simulation) when water vapour transport across the southern boundary is halved through the entire vertical column. These results indicate that the vast majority of the precipitation in this case (∼98%) depends on transport of moisture from the Indian subcontinent, with large contributions (56–82%) attributable to up-and-over transport. Additional simulations for other intrusive CSs show similar results (Supplementary Fig. 4), indicating that these contributions are robust.
The topographic features of the Himalayas, including both valleys and ridges, are well-represented by the 9-km grid used for the WRF simulations (Supplementary Fig. 5). Moisture transport can be found both in the valleys and at middle-to-upper levels, indicating that the upslope and up-and-over routes coexist. However, the reduction of moisture transport in the sensitivity simulations is most pronounced when advection of upper layer moisture is eliminated, indicating that moisture supply during these events depended most critically on up-and-over transport.
Our conclusions regarding the importance of up-and-over transport during the 12 September 2002 intrusive case are further supported by backward trajectory calculations using the FLEXPART particle dispersion model34 (Supplementary Fig. 6). Trajectory-derived transport pathways preceding the rain event are qualitatively consistent with the results of the WRF simulations described above. The quantitative partitioning of ‘up-and-over’ and ‘local convection’ trajectories is very sensitive to the criteria used to define these subsets; however, reasonable criteria assign 50–75% of these particles to the up-and-over route, indicating that transport via the up-and-over route was at the very least comparable to local convection (including that dependent on upslope moisture transport) during this event. Overall, the evidence consistently supports both the existence and the critical importance of the up-and-over transport to SWTP precipitation at event scales.
Moisture transport and isotopic observations
Both satellite measurements and reanalyses confirm the existence of mid-tropospheric moisture corridors extending from CEI into the SWTP associated with intrusive CSs (for example, Supplementary Figs 7 and 8). By contrast, the same moisture fields show distinct breaks near the Himalayas for non-intrusive cases (not shown). The consistency between AIRS retrievals and ERA-Interim estimates of water vapour in this case indicates that the reanalysis is able to adequately capture the circulation and moisture fields over this area, as AIRS radiances were not assimilated by ERA-Interim until April 2003 (ref. 35). The vertical structure of the reanalysis moisture fields and circulation along the 80° E longitude transect (Supplementary Fig. 8) during typical intrusive and non-intrusive cases shows significant differences in the vicinity of the Himalayas and the SWTP. During the intrusive case, moist air and condensate lofted by strong upward motion were advected northward to the SWTP in the 400–600 hPa layer. By contrast, advection at these levels was directed southward during the non-intrusive case. Upslope transport also tends to be enhanced during intrusive events (Supplementary Fig. 8), but the sensitivity simulations described above indicate that rainfall during this event depended primarily on simultaneous increases in up-and-over transport. Contrasts in upper level circulation and rainfall between these two cases are also reproduced by the high-Asia refined analysis (Supplementary Fig. 9).
Additional evidence for the moisture source associated with intrusive CSs over the SWTP can be gleaned from the stable oxygen isotope ratio (δ18O) in precipitation. This tracer has been widely used to infer key characteristics of precipitation, including moisture source regions, transport history and the degree of rainout24,25,26,27. We use event-based observations of δ18O in precipitation from four stations in the Tibetan Network for Isotopes in Precipitation (TNIP) and one station operated by the Institute of TP Research (ITP) during 1998–2007 to evaluate the contributions of different moisture sources and transport pathways to summertime rainfall over the SWTP. A total of 98 summertime rainfall events (59 intrusive cases and 39 non-intrusive cases) were selected corresponding to the identified CSs (see Supplementary Table 3 for details). As shown by previous studies, precipitation amount and δ18O are inversely related in this region (Fig. 4), with a lower ratio of heavy isotopes during events with higher precipitation amounts27. Moreover, we find that intrusive CSs were consistently associated with larger precipitation amounts and lower values of δ18O than non-intrusive CSs. This difference is consistent with differences in moisture source and transport effects, but is also potentially consistent with the well-known amount effect, which also manifests as an inverse relationship between δ18O and precipitation amount36,37. To distinguish between these two possibilities, we generate theoretical estimates of the amount effect in local convection (dashed curves in Fig. 4; see Methods for details) using a model intended to represent rain events due to either upslope transport or local recycling (which both deliver moisture to the atmospheric surface layer above the SWTP). The theoretical estimates are consistent with the observations collected during non-intrusive cases (when local convection is expected to dominate), but diverge sharply from the observations collected during intrusive cases. Ratios of δ18O in summer precipitation over the SWTP for intrusive cases are much larger than expected based on the theoretical curves. Previous studies have likewise shown that water vapour detrained from tropical convective clouds is more enriched in heavy isotopes than expected based on theoretical estimates38. The isotopic observations collected during intrusive cases are therefore consistent with efficient moisture transport from the tropical ocean via the up-and-over pathway, and inconsistent with heavy convective precipitation that depends on moisture sources that deliver water vapour primarily to the surface layer (that is, upslope flow and local evapotranspiration).
Composite and mean precipitation profiles
To extend our analysis beyond event scales, we construct composites for all intrusive and non-intrusive CSs during the period 1998–2005 (Fig. 5). Intrusive CSs account for approximately half of total summertime precipitation over both CEI (46%) and the SWTP (47%), with small interannual fluctuations. Precipitation occurs over the SWTP during 90% of the intrusive CSs, while rainfall is mostly confined to the Indian subcontinent during non-intrusive CSs. Non-intrusive CSs account for ∼35% of total summer rainfall over CEI, but with relatively little rainfall over the SWTP (17% of total summertime precipitation). Similar to the example cases discussed above, both composites have a low-level cyclonic circulation over CEI. CSs are more intense and extend further northwestward in the intrusive composite (Fig. 5a), while CSs and the monsoon trough are typically shifted slightly eastward towards the northwestern Bay of Bengal in the non-intrusive composite (Fig. 5d). The area of heavy rainfall during intrusive CSs extends towards the northeast in a band that stretches from CEI to the SWTP, but the area of heavy rainfall during non-intrusive CSs is smaller and mostly confined to CEI. Northward cross-barrier flow from CEI to the SWTP at 500 hPa is apparent in the intrusive composite (Fig. 5b), and contrasts with southward flow from the SWTP to CEI in the non-intrusive composite (Fig. 5e). An eastward centre of high pressure at 300 hPa straddles the southern periphery of the Plateau in the intrusive composite (Fig. 5c), while anticyclonic winds associated with the Iran High are much stronger and extend further eastward in the non-intrusive composite (Fig. 5f), suggesting that intrusive and non-intrusive CSs preferentially occur during the TP and Iranian Plateau modes of the upper tropospheric anticyclone, respectively39. This difference means that the upper level circulation over the SWTP during non-intrusive CSs is dominated by flow from the mid-latitudes rather than by flow from the Indian subcontinent. Moist static energy, which reflects the effects of convective heating and moistening, is likewise larger in the intrusive composite, with a southwest–northeast tilt that extends further over the SWTP30,40.
We apply the same method to create composite distributions of moisture anomalies in the middle and upper troposphere using twice-daily AIRS retrievals during 2002–2005 (Fig. 6). The intrusive composite features a tongue of positive moisture anomalies that extends northeastward from CEI to the SWTP in the middle and upper troposphere (Fig. 6a–c), consistent with the prevailing cross-barrier (southwesterly) flow. By contrast, positive moisture anomalies are confined mainly to the southern side of the Himalayas during non-intrusive CSs (Fig. 6d–f), consistent with the anomalous eastward extension of the Iran High and the associated anticyclonic flow. In summary, strong southwesterly flow promotes moisture transport over the SWTP during intrusive CSs, while strong northeasterly flow inhibits moisture transport over the SWTP during non-intrusive CSs.
Figure 7 shows mean precipitation along a transect perpendicular to the Plateau boundary for the intrusive and non-intrusive CS composites. Intrusive and non-intrusive CSs both result in large amounts of precipitation over the Indian subcontinent, but intrusive CSs produce much more precipitation over the SWTP than non-intrusive CSs. In particular, while orographic precipitation along the south slope of the plateau is apparent during both intrusive and non-intrusive CSs, precipitation extends northward into the Plateau during intrusive CSs but drops sharply at altitudes >2.5 km during non-intrusive CSs. This difference highlights the importance of up-and-over transport to precipitation over the SWTP. The up-and-over route is open during intrusive CSs, and enhances rainfall by delivering moisture and condensate to the inner part of the SWTP. By contrast, the up-and-over route is effectively closed during non-intrusive CSs, and precipitation rarely penetrates beyond the flank of the plateau.
Discussion
Our results emphasize that mid-tropospheric advection of convectively lofted hydrometeors and moisture from CSs in CEI northward over the Himalayas (that is, ‘up-and-over’ transport) is critically important for summertime precipitation variability over the SWTP. Relative to conventional upslope moisture transport, which is frequent but inefficient, up-and-over transport is more efficient but less frequent. Up-and-over transport occurs on ∼23% of days during summer and accounts for approximately half of total summer rainfall over the SWTP. Local convection also occurs frequently over the SWTP in summer, and follows a clear diurnal cycle with an afternoon peak. Episodes of local convection often follow intrusive CSs (Fig. 2b), when up-and-over transport is particularly strong. This lag may be explained by the importance of precipitation recycling processes13,20 in these episodes. Up-and-over transport, complemented by upslope transport, and associated rainfall represent the ultimate source of much of the moisture that participates in these recycling processes, and therefore contribute to total precipitation over the SWTP even on days when this transport is weak or non-existent. In this sense, up-and-over transport associated with intrusive CSs over India fundamentally determines the overall wetness of the SWTP: without this transport, the SWTP would be substantially drier. Although most precipitation over the SWTP occurs during summer, future studies should investigate the role of precipitation variability in other seasons (especially that associated with mid-latitude systems during winter) in the regional water budget.
We present strong evidence that intrusive CSs over the Indian subcontinent directly and significantly influence rainfall over the SWTP at both seasonal and event scales. This result has important implications for understanding and evaluating projections of future changes in regional precipitation to climate forcings, and for assessing the resilience of the biosphere, hydrosphere and cryosphere in this region. Confidence in projected changes in regional rainfall is notoriously difficult to achieve, but is vitally important for planning adaptation strategies, managing water resources and establishing environmental protections. Summer precipitation is projected to increase over both CEI and the SWTP under a warmer climate41, which would benefit ‘summer accumulation’ type glacier growth in this region3. Although changes in both regions remain uncertain, the importance and robust nature of the relationship between these two regions establishes a framework for evaluating model performance in this region, eliminating unrealistic simulations and ultimately bolstering confidence in projections of future change.
Methods
Core data
APHRODITE18 is a long-term continental-scale daily gridded precipitation data set over Asia. These data have previously been used in studies of precipitation characteristics over both India7 and the TP. In addition to APHRODITE, we use observational estimates of precipitation from the TRMM Multi-satellite Precipitation Analysis19, GPCP42, CMAP43 and model output from CMIP5 historical simulations and future projections. Moisture fields are based on satellite retrievals from AIRS44,45 and reanalysis output from ERA-Interim35. ERA-Interim reanalysis and forecast fields are also used for dynamic and thermodynamic conditions in both the case study and the construction of composites. Additional details are provided in Supplementary Tables 1 and 2.
Time series correlation coefficient
We calculate correlation coefficients using normalized anomalies (anomalies from the mean over the entire record divided by the s.d. over the entire record) of summer rainfall over CEI and the SWTP.
CS identification and categorization
We identify and classify CSs using 3-h CLAUS46 and TRMM19 data. First, CSs over CEI are detected based on a conventional definition13 as contiguous areas larger than 10,000 km2 with CLAUS brightness temperatures (Tb) colder than 219 K. Second, large rainfall events8 are identified as days with rainfall amounts larger than 6 mm per day. CS days are defined as days that satisfy both of the above criteria over CEI. We then further classify CS days into two types based on mid-tropospheric (600–400 hPa) wind directions within a region between 78 and 82° E longitude and 26 and 32° N latitude. CS days with northward flow stronger than seasonal mean (1 m s–1) in the mid-troposphere over the mountains are classified as intrusive, while CS days with southward flow in the mid-troposphere over the mountains are classified as non-intrusive. Finally, we use 850 hPa relative vorticity to track the location and intensity of each CS.
WRF simulations
WRF model simulations are conducted using one-way nesting approach via the ndown approach (Detailed model setup is listed in Supplementary Table 4). A high-resolution (9 km) inner domain is nested within a coarse-resolution (27 km) domain (Supplementary Fig. 3), with the lateral boundary conditions of the inner domain obtained from hourly outputs from the outer domain simulations. This approach enables us to alter hydrometeor and water vapour fluxes at the boundaries of the inner domain. In addition to the control simulation, three other sensitivity simulations are conducted to differentiate and quantify the contributions of moisture fluxes in different vertical layers to the rainfall amount over the SWTP. Advection of water vapour is reduced by 50% in the whole air column, above 5 km AGL and above 2.5 km AGL, respectively, in these three sensitivity simulations. Radiosonde measurements indicate that the middle and upper troposphere over CEI are nearly saturated (relative humidities of 70–100%) on CS days, whereas non-convective days are rather dry with relative humidities of 20–50% above the boundary layer (Supplementary Fig. 10). For simplicity, we reduce the specific humidity above 5 km AGL by 50% to mimic the situation without convection. This simulation is intended to quantify the contribution of up-and-over transport. The second simulation, in which we reduce moisture advection by 50% above 2.5 km AGL, is intended to quantify the contribution of upslope transport. The final simulation, in which we reduce moisture advection throughout the column by 50%, is intended to quantify the total contribution of transport across the southern boundary.
Isotope observations
The TNIP includes more than 30 observational stations on the TP, with the longest records dating back to the 1960s and the most recent records collected in 2007 (ref. 27). We select event-based precipitation data from four TNIP stations and one ITP station that measured δ18O ratios in precipitation over the southern or SWTP during 1998–2007. We select the station at Lulang, located in the southeastern TP, to compare differences in moisture sources for rainfall over the SWTP and the southeastern TP. A total of 98 summer rainfall events corresponding to our categorized CSs were observed at the five selected stations (Supplementary Table 3). Additional details regarding these measurements have been provided by Yao et al27.
Isotope theoretical model
The theoretical amount effect is modelled by treating local convection as a Rayleigh distillation process operating on a single parcel, with condensate immediately removed as precipitation after it forms36. Isotopic fractionation is assumed to occur at thermodynamic equilibrium. The liquid–ice transition is neglected for simplicity; its inclusion does not affect the qualitative nature of the results. The model includes five free parameters, the surface pressure (set to be 560 hPa) and the initial temperature (ranges from 260 to 280 K), relative humidity (set to 80%), vapour isotopic composition (δ18Ov=–15‰) and mass of the parcel. In the absence of suitable constraints, the last two parameters are used to fit the model to the observations. With the exception of δ18Ov, which acts to shift the model output left and right to lower or higher values of δ18O, variation of these parameters only affects the quantitative end points and does not affect the qualitative output of the model. The shapes of the theoretical curves, which are the basis for our conclusion that moisture source and transport effects govern the differences between intrusive and non-intrusive cases, are unaffected by these variations. The most important source of uncertainty in this model and its interpretation is the role of post-condensation exchange, as equilibration and partial evaporation of raindrops falling through an unsaturated cloud layer can have important impacts on the isotopic composition of precipitation at the ground37. These effects are dependent on raindrop size, relative humidity below the cloud base and the depth of the unsaturated layer, and are a key component of the overall amount effect. Reanalysis data indicate that the relative humidity of near-surface air over the SWTP is generally higher during intrusive cases than during non-intrusive cases. Given similar raindrop sizes, post-condensation exchange will therefore tend to increase δ18O in precipitation more during non-intrusive cases than during intrusive cases; however, the opposite would need to be true for post-condensation exchange to explain the deviation of the intrusive observations from the theoretical curve. Although the effects of post-condensation exchange likely cause us to overestimate δ18Ov relative to its actual value, they are unlikely to affect our main conclusion.
Backward trajectory analysis
Trajectories are calculated using version 9.02 of the FLEXPART particle dispersion model34 driven by ERA-Interim reanalysis and forecast fields35 at 1° × 1° horizontal resolution and 3-h temporal resolution. About 50,000 particles filling the domain 30–33° N, 80–84° E and 100–8,000 m AGL are initialized at 18UTC 12 September 2002 and tracked backward in time until 12UTC 8 September 2002 (102 h), with particle locations output hourly. To properly account for the complex topography and relative frequency of moist convection in this region, we conduct these simulations with parametrizations of moist convection, sub-grid scale turbulence and mesoscale velocity fluctuations enabled.
Additional information
How to cite this article: Dong, W. et al. Summer rainfall over the southwestern Tibetan Plateau controlled by deep convection over the Indian subcontinent. Nat. Commun. 7:10925 doi: 10.1038/ncomms10925 (2016).
References
Xu, X., Lu, C., Shi, X. & Gao, S. World water tower: An atmospheric perspective. Geophys. Res. Lett. 35, L20815 (2008).
Yao, T. et al. Different glacier status with atmospheric circulations in Tibetan Plateau and surroundings. Nat. Clim. Change 2, 663–667 (2012).
Mölg, T., Maussion, F. & Scherer, D. Mid-latitude westerlies as a driver of glacier variability in monsoonal High Asia. Nat. Clim. Change 4, 68–79 (2014).
Tao, S.-Y. & Ding, Y.-H. Observational evidence of the influence of the Qinghai-Xizang (Tibet) Plateau on the occurrence of heavy rain and severe convective storms in China. Bull. Am. Meterol. Soc. 62, 23–30 (1981).
Dorji, T. et al. Plant functional traits mediate reproductive phenology and success in response to experimental warming and snow addition in Tibet. Global Change Biol. 19, 459–472 (2013).
Maussion, F. et al. Precipitation seasonality and variability over the Tibetan Plateau as resolved by the High Asia Reanalysis. J. Clim. 27, 1910–1927 (2014).
Turner, A. G. & Annamalai, H. Climate change and the South Asian summer monsoon. Nat. Clim. Change 2, 587–595 (2012).
Singh, D., Tsiang, M., Rajaratnam, B. & Diffenbaugh, N. S. Observed changes in extreme wet and dry spells during the South Asian summer monsoon season. Nat. Clim. Change 4, 456–461 (2014).
Chen, B., Xu, X., Yang, S. & Zhang, W. On the origin and destination ofatmospheric moisture and air mass over the Tibetan Plateau. Theoret. Appl. Climatol. 110, 423–435 (2012).
Annamalai, H. & Slingo, J. M. Active/break cycles: diagnosis of the intraseasonal variability of the asian summer monsoon. Clim. Dynam. 18, 85–102 (2001).
Wu, G. & Zhang, Y. Tibetan Plateau forcing and the timing of the monsoon onset over South Asia and the South China Sea. Mon. Weather Rev. 126, 913–927 (1998).
Yang, K. et al. Recent climate change over the Tibetan Plateau and their impacts on energy and water cycle: A review. Glob. Planet. Change 112, 79–91 (2014).
Sugimoto, S. & Ueno, K. Formation of mesoscale convective systems over the eastern Tibetan Plateau affected by plateau-scale heating contrast. J. Geophys. Res. 115, D16105 (2010).
Wang, B. et al. Northern hemisphere summer monsoon intensified by mega-El nino/southern oscillation and atlantic multidecadal oscillation. Proc. Natl Acad. Sci. USA 110, 5347–5352 (2013).
Bollasina, M. A., Ming, Y. & Ramaswamy, V. Anthropogenic aerosols and the weakening of the South Asian summer monsoon. Science 334, 502–505 (2011).
Lau, K. M., Kim, M. K. & Kim, K. M. Asian summer monsoon anomalies induced by aerosol direct forcing: the role of the Tibetan Plateau. Clim. Dynam. 26, 855–864 (2006).
Ramanathan, V. et al. Atmospheric brown clouds: Impacts on South Asian climate and hydrological cycle. Proc. Natl Acad. Sci. USA 102, 5326–5333 (2005).
Yatagai, A. et al. APHRODITE: constructing a long-term daily gridded precipitation data set for Asia based on a dense network of rain gauges. Bull. Am. Meteorol. Soc. 93, 1401–1415 (2012).
Huffman, G. J. et al. The TRMM Multisatellite Precipitation Analysis (TMPA): Quasi-global, multiyear, combined-sensor precipitation estimates at fine scales. J. Hydrometeorol. 8, 38–55 (2007).
Liu, X., Bai, A. & Liu, C. Diurnal variations of summertime precipitation over the Tibetan Plateau in relation to orographically-induced regional circulations. Environ. Res. Lett. 4, 045203 (2009).
Xie, S., Xu, H., Saji, N. H., Wang, Y. & Liu, W. Role of narrow mountains in large-scale organization of Asian monsoon convection. J. Clim. 19, 3420–3429 (2005).
Numaguti, A. Origin and recycling processes of precipitating water over the Eurasian continent: experiments using an atmospheric general circulation model. J. Geophys. Res. 104, 1957–1972 (1999).
Yang, M., Yao, T., Wang, H., Gou, X. & Tian, L. Estimating the criterion for determining water vapour sources of summer precipitation on the northern Tibetan Plateau. Hydrol. Process. 20, 505–513 (2005).
Tian, L. et al. Oxygen-18 concentrations in recent precipitation and ice cores on the Tibetan Plateau. J. Geophys. Res. 108, 4293 (2003).
Tian, L. et al. Stable isotopic variations in west China: A consideration of moisture sources. J. Geophys. Res. 112, D10112 (2007).
Kurita, N. & Yamada, H. The role of local recycling evaluated using stable isotope data from over the middle of the Tibetan Plateau during the monsoon season. J. Hydrometeor. 9, 760–765 (2013).
Yao, T. et al. A review of climatic controls on δ18O in precipitation over the Tibetan Plateau: Observations and simulations. Rev. Geophys. 51, 525–548 (2013).
Gadgil, S. The Indian monsoon and its variability. Annu. Rev. Earth Pl. Sc. 31, 429–467 (2003).
Romatschke, U. & Houze, R. A. Jr Characteristics of precipitating convective systems in the South Asian monsoon. J. Hydrometeor. 12, 3–26 (2011).
Chen, T. & Weng, S. Interannual and intraseasonal variations in monsoon depressions and their westward-propagating predecessors. Mon. Wea. Rev. 127, 1005–1020 (1999).
Krishnamurthy, V. & Shukla, J. Intraseasonal and seasonally persisting patterns of Indian monsoon rainfall. J. Clim. 20, 3–20 (2007).
Prajeesh, A. G., Ashok, K. & Rao, D. V. B. Falling monsoon depression frequency: A Gray-Sikka conditions perspective. Sci. Rep. 3, 2989 (2013).
Molnar, P., Boos, W. R. & Battisti, D. S. Orographic controls on climate and paleoclimate of Asia: thermal and mechanical roles for the Tibetan Plateau. Annu. Rev. Earth Planet. Sci. 38, 77–102 (2010).
Stohl, A. et al. Technical note: the Lagrangian particle dispersion model FLEXPART version 6.2. Atmos. Chem. Phys. 5, 2461–2474 (2005).
Dee, D. P. et al. The ERA-Interim reanalysis: configuration and performance of the data assimilation system. Quart. J. R. Meteorol. Soc. 137, 553–597 (2011).
Dansgaard, W. Stable isotopes in precipitation. Tellus 16, 436–468 (1964).
Lee, J.-E. & Fung, I. ‘Amount effect’ of water isotopes and quantitative analysis of post-condensation processes. Hydrol. Process. 22, 1–8 (2007).
Webster, C. R. & Heymsfield, A. J. Water isotope ratios D/H, 18O/16O, 17O/16O in and out of clouds map dehydration pathways. Science 302, 1742–1745 (2003).
Zhang, Q., Wu, G. & Qian, Y. The bimodality of the 100 hPa South Asia high and its relationship to the climate anomaly over East Asia in summer. J. Meteorol. Soc. Japan. Ser. II 80, 733–744.
Boos, W. R. & Kuang, Z. Dominant control of the South Asian monsoon by orographic insulation versus plateau heating. Nature 463, 218–223 (2010).
Collins, M. et al. in: Climate Change 2013: The Physical Science Basis eds Stocker T. F.et al. Cambridge Univ. Press (2013).
Adler, R. F. et al. The version 2 Global Precipitation Climatology Project (GPCP) monthly precipitation analysis (1979-Present). J. Hydrometeor. 4, 1147–1167 (2003).
Xie, P. & Arkin, P. A. Global precipitation: a 17-year monthly analysis based on gauge observations, satellite estimates, and numerical model outputs. Bull. Am. Meteor. Soc. 78, 2539–2558 (1997).
Aumann, H. H. et al. AIRS/AMSU/HSB on the Aqua mission: Design, science objectives, data products, and processing systems. IEEE T. Geosci. Remote 41, 253–264 (2003).
Tian, B. et al. AIRS/AMSU/HSB version 6 level 3 product user guide. tech. rep., jet propulsion laboratory, Pasadena. Available at http://disc.sci.gsfc.nasa.gov/AIRS/documentation/v6_docs/v6releasedocs-1/V6_L3_User_Guide.pdf (2013).
Hodges, K. I., Chappell, D. W., Robinson, G. J. & Yang, G. An improved algorithm for generating global window brightness temperatures from multiple satellite infrared imagery. J. Atmos. Oceanic Technol. 17, 1296–1312 (2000).
Acknowledgements
We gratefully acknowledge the APHRODITE, TRMM, GPCP, CMAP, CLAUS, AIRS and ERA-Interim projects for preparing the precipitation, brightness temperature and other atmospheric data, and the associated data repositories (listed in Supplementary Table 2) for providing public access to these data. We further thank the climate modelling groups for producing the model output and making it available. This work was supported by the Ministry of Science and Technology of China (Grant 2013CBA01805 and 2014CB441303). L.W. is supported by the National Natural Science Foundation of China (Grant 41322001) and the Hundred Talents Program of Chinese Academy of Sciences. L.T. is supported by the National Natural Science Foundation of China (Grant 41530748).
Author information
Authors and Affiliations
Contributions
W.D. and Y.Lin. developed the essential research idea. W.D. performed the analysis and wrote the initial draft of the paper. Y.Lin., Y.X. and J.S.W. contributed editorial suggestions and supporting analyses. L.T. and L.W. provided the isotope measurements and their interpretation. Y.M., B.W., Y.Luo, W.H., J.H., Y.P., and F.X. contributed to the interpretation of the results and the preparation of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Figures 1-10, Supplementary Tables 1-4 and Supplementary References. (PDF 1684 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Dong, W., Lin, Y., Wright, J. et al. Summer rainfall over the southwestern Tibetan Plateau controlled by deep convection over the Indian subcontinent. Nat Commun 7, 10925 (2016). https://doi.org/10.1038/ncomms10925
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms10925
This article is cited by
-
Variations of Stable Isotopic Composition in Precipitation and their Controlling Factors, a Case Study in Dongying, Yellow River Delta
Asia-Pacific Journal of Atmospheric Sciences (2024)
-
Westerlies Affecting the Seasonal Variation of Water Vapor Transport over the Tibetan Plateau Induced by Tropical Cyclones in the Bay of Bengal
Advances in Atmospheric Sciences (2024)
-
Controls of the recent precipitation anomalies in the southeastern Tibetan Plateau: from the perspective of Indian summer monsoon activities and moisture sources
Climate Dynamics (2024)
-
Anthropogenic forcing and Pacific internal variability-determined decadal increase in summer precipitation over the Asian water tower
npj Climate and Atmospheric Science (2023)
-
The enhanced relationship between summer rainfall over the eastern Tibetan Plateau and sea surface temperature in the tropical Indo-Pacific Ocean
Climate Dynamics (2023)