Observed Soil Moisture Impact on Strong Convection over Mountainous Tibetan Plateau

: Convectionover the TibetanPlateau(TP) has beenlinked to heavyrain and ﬂoodingin downstream parts of China. Understanding processes which inﬂuence the development of convection on the TP could contribute to better forecastingof these extreme events.TP scale( ; 1000km) soilmoisture gradientshave beenshown to inﬂuence formation of convective systems over the eastern TP. The importance of smaller-scale ( ; 10km) variability has been identiﬁed in other regions (including the Sahel and Mongolia) but has yet to be investigated for the TP. In addition, compared to studies over ﬂat terrain, much less is known about soil moisture–convection feedbacks above complex topography. In this study we use satellite observations of cold cloud, land surface temperature, and soil moisture to analyze the effect of mesoscale soil moisture heterogeneity on the initiation of strong convection in the complex TP environment. We ﬁnd that strong convection is favored over negative (positive) land surface temperature (soil moisture) gradients. The signal is strongest for less vegetation and low topographic complexity, though still signiﬁcant up to a local standard deviation of 300m in elevation,accountingfor65%ofcases.Inaddition,thesignalisdependentonbackgroundwind.Strongconvectiveinitiationisonlysensitivetolocal(tensofkilometers)soilmoistureheterogeneityforlightwindspeeds,thoughlarge-scale(hundreds of kilometers) gradientsmay still be important for strong wind speeds. Our results demonstrate that, even in the presence of complex local soil moisture variability plays an important role in storm development.


Introduction
The Tibetan Plateau (TP) is the highest and most extensive plateau in the world, profoundly affecting climate and weather in the region (Yang et al. 2013;Li et al. 2014). Due to its average elevation of more than 4000 m, it provides strong thermal and dynamical forcing in the midtroposphere during the summer months, fostering the frequent development of intense storms (e.g., Yang et al. 2004). Tibetan convective systems (TCSs) can be associated with particularly extreme rainfall events and contribute up to ;70% of rainfall over the TP and adjacent areas (Hu et al. 2016). In particular, propagation of TCSs can bring heavy rain and flooding to downstream parts of China, affecting millions of people (e.g., Li et al. 2008;Zhao et al. 2019). A better understanding of the processes that impact TCS genesis over the TP could contribute to better forecasting of these extreme events. Furthermore, there is strong evidence for accelerated climate warming on the TP (e.g., Guo and Wang 2012;Duan and Xiao 2015). This warming may affect convection on the TP, and hence hazardous weather in China, which makes the identification of factors for TCS development even more important.
Numerous works have analyzed characteristics of TCSs using various satellite products (e.g., Guo et al. 2006;Li et al. 2008;Hu et al. 2016;Li et al. 2019;Gao et al. 2019). These studies have revealed two core regions for convective activity during the monsoon season, the southern-central plateau and the eastern plateau. The high frequency of TCSs over the southern-central plateau has been attributed to topography (Kurosaki and Kimura 2002;Barros et al. 2004;Fujinami et al. 2005). Large-scale (.100 km) mountain ranges have an organizing effect on afternoon strong convection. However, the correlation between cloud distribution and topography with horizontal scales of less than 100 km is less obvious (Kurosaki and Kimura 2002). Indeed, that study suggested that, at the mesoscale, the spatial distribution of cloud may be controlled by local thermally induced circulations. Sugimoto and Ueno (2010), using a combination of satellite observations and model simulations, have demonstrated that the high frequency of TCSs over the eastern plateau is linked to a longitudinal surface wetness gradient. Sub-plateau-scale convergences are established by plateau scale heating contrasts, Denotes content that is immediately available upon publication as open access.
Supplemental information related to this paper is available at the Journals Online website: https://doi.org/10.1175/JHM-D-20-0129.s1. effectively forming TCSs in the eastern plateau during the monsoon season. The authors commented that further work was needed to assess sensitivity to sub-plateau-scale surface wetness patterns.
The impact of soil moisture (SM) on moist convection has been recognized by numerous works (e.g., see reviews by Seneviratne et al. 2010;Santanello et al. 2018). SM affects the partitioning of surface fluxes into sensible and latent heat, which control planetary boundary layer (PBL) growth and moisture availability, respectively. Depending on the atmospheric stability profile, either can have a direct effect on the formation of clouds and deep convection (e.g., Findell and Eltahir 2003;Ek and Holtslag 2004). On the TP the importance of these (one dimensional) thermodynamic triggers of convection has been investigated for the Nam Co Lake Basin. A high-resolution cloud-resolving simulation revealed a clear positive correlation between SM, latent heat flux and cloud cover (Gerken et al. 2015). That study also identified a convective optimum at intermediate SM, where moisture availability was balanced with sufficient sensible heat flux to drive vertical motion.
Spatial variations in near-surface air pressure and low-level convergence, linked to SM, can also influence convective initiation as observed by Sugimoto and Ueno (2010) for the TP at large scales, though this effect can also be important at smaller scales (e.g., Pielke 2001;Taylor et al. 2007). Observational studies have shown that horizontal circulations driven by mesoscale land surface heterogeneity enhance the probability of convective initiation in various regions, including the Sahel, Europe, and Mongolia (Taylor et al. 2011, hereafter T11;Taylor 2015, hereafter T15;Teramura et al. 2019). The importance of this (two dimensional) dynamical process has also been identified over complex terrain in southwest Germany (Kottmeier et al. 2008;Corsmeier et al. 2011;Hauck et al. 2011) but has yet to be investigated for the TP.
In this study we use satellite observations of cloud-top temperature (CTT), land surface temperature (LST), and SM to examine the influence of mesoscale SM variability on the initiation of strong convection during summer months June-September (JJAS). This paper is organized as follows. Section 2 describes the study region, data, and methods. Section 3 details our results including subset analyses where we consider the role of topographic environment, background wind conditions, and vegetation. This is followed by a discussion in section 4 and a list of our main conclusions in section 5.

a. Study area and climatology
The TP covers 2.5 million km 2 between 258 and 458N and between 708 and 1058E (Fig. 1). Elevations on the plateau range from 2500 to over 8000 m (Fig. 1a) with local topographic height variations exceeding 400 m in places (Fig. 1b). Known as the water tower of Asia (Immerzeel et al. 2010), it contains nearly 400 lakes with aerial extent greater than 10 km 2 ( Fig. 1a; Wan et al. 2017). In the summer there exists a clear northwest-southeast gradient in vegetation linked to precipitation, with highest normalized difference vegetation index (NDVI), SM, and rainfall total occurring in the southeastern plateau (Figs. 1c,d).

b. Datasets
The datasets used in the current work are described in detail in this section and summarized in Table S1 in the online supplemental material. Topographic variability is computed from the Shuttle Radar Topography Mission 1 arc second global dataset (SRTM; Farr et al. 2007), vegetation state is interpreted from the MODIS Terra/Aqua monthly NDVI product at 0.058 (MOD13C2/MYD13C2; Didan 2015a,b) and the location and size of lakes (with area . 10 km 2 ) are taken from Wan et al. (2017).
Strong convection is identified from the Fengyun-2 (FY2) CTT product (NSMC 2013). FY2 is a series of geostationary satellites located at 08, 798-1128E. The CTT product is available from 2013 to present and provides hourly (half-hourly for 2019 onward) images at ;5-km resolution. We use data from FY2F (2013-15), FY2G (2016-18), and FY2H (2019). For each year we selected data from the most recently launched satellite because of (i) the improved data continuity of FY2G compared to FY2F and (ii) the higher temporal resolution of FY2H compared to FY2G. To ensure consistency between satellites we inspect images from years with more than one satellite. For FY2H/FY2G, coincident observations produce the same results in our identification of strong convective initiation (see section 2c). For FY2G/FY2F it is not possible to do a direct comparison due to a change in observation time (from 30 min past the hour to on the hour). However, the evolution of the cloud field in multisatellite simulations appears consistent. Due to the large viewing angles of the satellites it is necessary to correct for parallax. For the correction, cloud-top height is estimated using lapse rates from ERA5 reanalysis (0.258, 1-h time step; Hersbach et al. 2019). Wind vectors at 10 m from ERA5 are also used to determine the direction and speed of background wind. Winds on the TP are highly variable in space, and the observational network is sparse. We therefore test the sensitivity of our results to the choice of wind reanalysis product, with parts of our analysis repeated using MERRA2 (0.58 3 0.6258, 1-h time step; Gelaro et al. 2017) and ERA-Interim (0.758, 6-h time step; Berrisford et al. 2011). To assess the statistical relationship between CTT and rainfall intensity we employ instantaneous radar measurements of surface rainfall rate from TRMM (2013, TRMM_2A25; TRMM 2011a) and GPM (2014-19;GPM_2ADPR;Iguchi and Meneghini 2017).
Prestorm surface wetness conditions are characterized from satellite observations of LST and SM. Due to the higher quality of LST observations (lower uncertainties, less missing data) compared to SM, the bulk of our analysis uses LST. Previous works (e.g., T11; T15) have used LST anomalies (LSTAs) as a proxy for surface wetness conditions. Locally positive (negative) LSTAs are indicative of locally drier (wetter) conditions compared to the surroundings. This study employs LSTA observations from the 10-km AMSR2 downscaled LPRM Level 2 product (overpass time 1530 LT, UTC 1 8; de Jeu and Owe 2014; Owe et al. 2008). The microwave radiometer provides all-sky observations; however, the 36.5-GHz channel used for LST retrieval is sensitive to rain droplets in the atmosphere Holmes et al. 2009). To minimize rainfall contamination, the dataset is filtered using GPM IMERG Level 3 high quality precipitation based on microwave-only data (including AMSR2 images; Huffman et al. 2019). All pixels with rainfall rates greater than 1 mm h 21 between 1500 and 1700 LT (UTC 1 8) are excluded. A monthly mean LST climatology per pixel is constructed for the period 2013-19 on the 10-km grid from which daily LSTAs are computed. The AMSR2 LPRM product also includes SM observations; however, this dataset contains a discontinuity in 2016 affecting SM values globally (note this discontinuity is not present in the LST time series). Instead we utilize independent SM observations contained within the ESA Climate Change Initiative (ESA-CCI) SM version 03.3 daily 25-km product (ESA-CCI C3S 2018). Although there are uncertainties and a large number of data gaps in the ESA-CCI product over the TP, there is a significant correlation between positive SM anomalies and observed precipitation in summer (Meng et al. 2018). However, the current application requires prestorm, rather than daily, SM values. Therefore pixels containing only midday (1230 LT, UTC 1 8) ASCAT observations are extracted and filtered. Only pixels with the highest quality flag of 0 and no rainfall contamination are retained. The resulting dataset is relatively sparse but does contain sufficient information to test our interpretation of results based on LST. As for LST, we compute a 7-yr SM climatology from which daily SM anomalies (SMA) are calculated. An evaluation study of remotely sensed and reanalysis SM products found that variations in ASCAT SM correlated well with in situ measurements (Zeng et al. 2015). Positive (negative) ASCAT SMAs should therefore be representative of wetter (drier) soils.

c. Identification of strong convective initiation locations
We define the development of strong convection as a newly emergent, rapidly cooling contiguous cloud area less than 30 km across, with CTT less than 2548C. The coordinate of the coldest pixel within the cloud area (corrected for parallax) is taken as the initiation point. Previous works have used rapidly cooling CTTs to identify convective initiation (e.g., Morel and Senesi 2002;T15), and cloud areas with spatial scales smaller than 35 km have been identified as having an increased probability of being associated with deep convection (Klein et al. 2018). There is some disparity in the literature on the temperature threshold for strong convection over the TP. Even very recent studies (e.g., Liu et al. 2019;Li et al. 2019), which use the same observational data (FY2), have opted for subtly different thresholds. The choices are 2548C (e.g., Liu et al. 2019;Sugimoto and Ueno 2010;Guo et al. 2006), 2538C (e.g., Li et al. 2008;Hu et al. 2016), and 2528C (e.g., Li et al. 2019).
A comparison of contiguous cold cloud (FY2) larger than 350 km 2 (following Klein et al. 2018) and coincident radar precipitation (TRMM_2A25 and GPM_2APDR) indicates that the probability of intense rain increases with decreasing CTT (Fig. 2). We therefore opt for the coldest threshold from the literature in our definition of strong convection. To minimize ''false initiations,'' in other words the detection of preexisting cloud areas associated with propagating cold cloud, for each case we examine surrounding pixels within 50 km of the detected cloud area in the previous time step. We reject cases where the decrease in minimum CTT is less than 108 over an hour (following T15).
On the TP it is vital to consider the impact of topography and lakes on the development of strong convection. To minimize the impact of lake triggers on our results, we exclude cases that occur within 20 km of water bodies greater than 10 km 2 (;20% of cases are removed by this filter). The thresholds were chosen to maximize the number of retained cases while reducing the impact of lakes on the analysis. Increasing the proximity threshold to 50 km (as was done in T15) results in a large (.50%) reduction in cases. Considering topography on the other hand, it is known from previous studies (e.g., Barros et al. 2004;Sugimoto and Ueno 2010;Li et al. 2019) that there is a higher frequency of convective activity over the southerncentral plateau, where topographic complexity is high. Rather than confine our analysis to relatively flat areas of the plateau, which would exclude this convective hotspot, we choose not to initially filter out cases in areas of strong topographic complexity. Instead we perform an additional analysis (described in section 3b) to explicitly examine the impact of topographic variability on the SM-convection coupling. Figure 3 displays the 7-yr season mean spatial distribution and diurnal cycle of strong convective initiations (excluding lake cases) identified using our method for the entire season (JJAS). The same figures by month are provided in the supplemental material (Fig. S1). We observe the same core regions of convective activity as previous works, the southern-central and eastern plateau, although we find the former to be slightly stronger (Fig. 3a). This is similar to the findings of Li et al. (2019), who also used FY2 data to identify strong convection over the TP. We observe a peak in strong convective initiations in the afternoon (Fig. 3b), around 50% of all storms form between the hours of 1400-2000 LT (0600-1200 UTC). This is similar to the diurnal cycle of summer TCS initiations reported by Li et al. (2008). To ensure we are sampling strong convection that initiated after our observations of LST (1530 LT, UTC 1 8), we choose to focus our analysis on 1700-2000 LT (0900-1200 UTC), around 25% of all storms form between these hours (Fig. 3b). We assume these later storms have similar characteristics and formation mechanisms as those that develop earlier in the afternoon.

d. Characterization of prestorm surface wetness
To calculate mean prestorm surface wetness conditions we adopt the method used in T11 and T15. First, the initiation location for each strong convection case is matched to the closest grid point in the surface data (LSTA or SMA). A surface domain (6200 km relative to the initiation location) is then extracted and rotated in the direction of the background wind (using 10-m wind vectors from reanalysis). Finally, a composite of mean surface conditions is generated from an average of rotated surface domains for all cases. To ensure a fair comparison, the spatial mean LSTA/SMA is subtracted from the individual surface domains before averaging. It is necessary to consider the potential influence of missing data (particularly relevant for the SM data) on our mean surface composites. Therefore, to reduce the impact of sampling issues, we only include cases where the surface domain contains at least 50% (75%) valid AMSR2 (ASCAT) observations. Our resulting sample size is 1435 (172) strong convection initiation cases for LSTA (SMA).

a. Mean conditions
First, we consider our full initiation sample to analyze the mean surface conditions associated with afternoon strong convection, presented in Fig. 4. The two-dimensional mean LSTA composite depicts a relatively warmer area upwind of the initiation location, with cooler conditions downwind (Fig. 4a). The preference for convection to develop on the downwind edge of dry patches, close to dry-wet boundaries, has been identified in modeling studies (e.g., Garcia-Carreras et al. 2011) and observed previously in the Sahel (T11) and Europe (T15). The two dimensional patterns of LSTA observed here show particular resemblance to the Sahel result, with positive anomalies over and upwind of the initiation location and negative anomalies around the edge of the domain.
Taking a downwind cross section through the center of the domain, we find the one dimensional LSTA gradient is strongest at and downwind of the initiation location (Fig. 4b).
To assess the significance of the mean negative downwind LSTA gradient, we adopt the method from T15. For strong convection cases where there is sufficient local LSTA data, we sample an initiation gradient between 210 and 20 km (solid black rectangle in Fig. 4a, shading in Fig. 4b), corresponding to the scale around which we expect the effect of surface heterogeneity on initiations to be maximized (T11). Then, for the same cases, we sample up to eight non-initiation gradients (solid gray rectangles in Fig. 4a). These gradients are similarly computed over 30 km downwind using a linear regression (following T11/T15). Overall this gives us a set of n i 5 1194 initiation gradients, with a mean of 20.68 K (100 km) 21 , and n ni 5 4071 non-initiation gradients with a mean of 20.02 K (100 km) 21 (5265 gradients in total). To test whether the mean initiation gradient is significantly more negative than the mean noninitiation gradient, we perform a single-tail p test on a distribution of random D values, where D is the difference between mean initiation and non-initiation gradients. This distribution is computed from 100 000 random combinations of n i and n ni gradients from the pool of 5265 gradients (initiation plus non-initiation). Based on this test we find the mean initiation gradient to be statistically significant (p value of 0.02). Therefore, we can conclude that strong convection on the TP occurs preferentially over negative LSTA gradients.
To test the role of soil moisture conditions in creating observed LSTA patterns, we generate a mean SMA composite. A downwind cross section of mean composite SMA is presented in Fig. 4b. For the 30-km section over which initiation gradients are computed, downwind gradients in LSTA are anticorrelated with gradients in SMA, although the correlation is not consistent along the whole transect. This could be attributed to the limited number of valid pixels (only 172 cases with sufficient local SM data compared to 1435 cases for LST) and the lower spatial resolution of the SM product (25 km for SM versus 10 km for LST). If we consider only the 172 cases with sufficient local SM data, the two dimensional composites of SMA and LSTA have broadly similar structures (Fig. S2).
We also consider the potential role of preinitiation cloud cover on LSTA patterns. This is necessary since, in contrast to T11 and T15, the present study makes use of microwave LST observations, which provide data for both cloudy and cloud-free days. The LSTA signal may therefore be affected by variations in incoming shortwave radiation at the surface. Here we consider the 651 cases with FY2 cloud observations at time of the AMSR2 overpass (1530 LT, UTC 1 8). This is less than the total number of cases due to a change in satellite observation time (see section 2). For each pixel in the 400 km 3 400 km domain we determine the fraction of cases where cloud is present (according to FY2) and examine the resulting downwind cross section with comparison to mean composite LSTA (Fig. S3a). We observe slightly less cloud (;4%) upwind of the initiation location coincident with warmer LSTA, which may be contributing to the observed LSTA gradient. To quantify the potential impact of this difference in cloudiness and surface insolation on our LSTA gradients, we first estimate relationships between LSTA and cloud cover fraction (CF) and SMA separately using linear regressions. Variables are sampled over the 50 km 3 50 km area of composite maximum positive LSTA (purple box with dotted outline in Fig. 4a). Mean LSTA and CF are sampled for the 651 cases with suitably timed cloud observations, while mean LSTA and mean SMA are sampled for the 172 cases with sufficient SM observations. In this instance CF is defined as the fraction of cloudy pixels in the 50 km 3 50 km area. As expected (Fig. S4), both drivers are negatively correlated with LSTA, though the variance explained by SMA (r 2 5 27%) is more than twice as large as for CF (r 2 5 12%). Using the slopes of these regression lines, we predict the magnitude of downwind LSTA variations induced by SM and cloudiness (Fig. S3b). The magnitude of the SMA derived signal (from 172 cases) is more than twice as large as the effect due to cloud variability and is consistent in both amplitude and location with the observed LSTA, which is based on 1435 cases. Overall these tests indicate that SM, rather than cloudiness, is the dominant driver of our observed gradients in LSTA, though the impact of cloud variability likely reinforces the observed gradient.

b. Dependence on topographic complexity
On the TP, topography can help to mechanically lift unstable near-surface air, potentially triggering deep convection independently from SM states. During nighttime and morning hours, the fraction of initiation cases associated with high topographic complexity ranges between 39% and 45% (Fig. S5b). This decreases to 25%-37% in the afternoon, suggesting an increased importance of initiation mechanisms linked to surface heating.
To investigate the sensitivity of the all case signal to topographic complexity, we divide our sample into subsets according to the standard deviation in elevation computed over a 40 km 3 40 km box (SD) centered on the initiation location (dashed box in Fig. 4a). The scale at which topographic complexity is considered was chosen to encompass the area within which initiation gradients are computed. Downwind cross sections are presented in Figs. 5a-d. Statistically significant (p , 0.05) negative downwind LSTA gradients are found for the two subsets with relatively modest topographic variability (SD less than 300 m), corresponding to 65% of the total dataset (Fig. S5a). Changing the box size for the calculation of topographic complexity alters the magnitudes of SD and redistributes some of the cases. We find our chosen box size is more effective than larger box sizes at identifying cases where soil moisture is playing a significant role (Table S2). A weakening of the soil moisture signal with increasing topographic complexity is consistent with observations over East Mongolia (Teramura et al. 2019) and results from idealized model experiments (Imamovic et al. 2017).

c. Sensitivity to background wind speed and dataset
Background wind affects the persistence of surface-induced variability in the PBL and may change the propagation speed of initiated convection (Froidevaux et al. 2014). At the same time, wind speed and direction are variables difficult to correctly capture in reanalysis data for a complex region such as the TP (Yu et al. 2019). In the following, we therefore assess the sensitivity of our signal to wind conditions using three different reanalysis datasets.
To assess the impact of background wind speed on the SM-convection coupling signal, we consider only the 925 cases with relatively modest local topographic complexity (SD , 300 m), for which we found a significant relationship. This sample is stratified into approximately equally sized subsets according to 1200 LT (0400 UTC) low-level wind speed using 10-m wind vectors from ERA5. Downwind cross sections are presented in Figs. 6a-c. Statistically significant (p # 0.05) negative downwind LSTA gradients are found for all subsets; however, the length scale of the upwind warmer area appears to increase with increasing wind speed. A local (30-50 km) feature is only apparent for the two subsets with wind speeds lower than 2.5 m s 21 . This is in agreement with results from large-eddy simulations which showed that wind speeds greater than 2.5 m s 21 suppressed the impact of surface heterogeneity on length scales up to 40 km (Avissar and Schmidt 1998). For high wind cases our results indicate that, although atmospheric sensitivity to small-scale heterogeneity (tens of kilometers) is suppressed, there is still a preference for convective initiation over larger scale (hundreds of kilometers) negative LSTA gradients.
To test the sensitivity of our results to the choice of reanalysis product, we recomputed the background wind speed subset (high topographic complexity cases removed) and all case LSTA composites using 10-m wind vectors from MERRA2 [0400 UTC (1200 LT)] and ERA-Interim [0600 UTC (1400 LT)]. Downwind cross sections are presented in Figs. 7a-d and mean 30-km initiation gradients are presented in Table 1. The one dimensional all case LSTA composites based on all three products are very similar (Fig. 7a) and the two dimensional composites also have consistent patterns (not shown). This gives us confidence that our mean result is not significantly impacted by uncertainties in modeled wind direction. For the background wind speed subset analysis, LSTA composites based on more recent products MERRA2 and ERA5 display similar behavior (Figs. 7b-d), particularly around the initiation location, although the initiation gradients derived from MERRA2 results are weaker and less significant (Table 1). In the case of ERA-Interim, the LSTA composites for moderate and high wind speeds display similar behavior to the MERRA2 and ERA5 results; however, this is not true for low wind speeds (Fig. 7b). This may be due to higher uncertainties in wind direction in the older model for this wind class. Stronger winds are more likely to be the result of persistent large-scale dynamical forcing, while lighter winds are more likely to be controlled by more rapidly evolving local drivers. The lower spatial and temporal resolution of ERA-Interim could mean the model is unable to accurately capture the latter. Indeed a global evaluation of daily averaged surface wind speeds found that ERA-Interim (along with two other reanalysis products with a 6-hourly time step) showed poor performance compared to an average of multiple reanalysis products (including ERA5 and MERRA2; Ramon et al. 2019). A comparison of wind direction between the products for our cases reveals a larger spread for lower wind speeds (Fig. S6), as expected. Both ERA-Interim and MERRA2 display a similar disagreement with ERA5, although we find the composite LSTA gradient to be more significant for MERRA2 than ERA-Interim (p 5 0.07 and p 5 0.6 respectively, Table 1) suggesting there may be other factors contributing to the worse skill of ERA-Interim. Wind direction between the three products is most consistent for the high wind speed class, where we find significant largescale negative downwind LSTA gradients regardless of product used (p # 0.05, Table 1).

d. Impact of vegetation
The presence of vegetation can suppress the impact of SM on surface flux partitioning and potentially weaken the coupling between SM and convection (Williams and Torn 2015). To determine the effect of vegetation on the local SMconvection coupling signal, we consider the 647 cases with relatively modest local topographic complexity (SD , 300 m) and lower background wind speeds (s , 2.5 m s 21 ). This sample is divided into low and high vegetation classes of approximately equal size according to MODIS Terra/Aqua monthly NDVI. Due to the northwest-southeast gradient in vegetation on the plateau, cases are broadly separated into groups of southeastern (high vegetation) and northwestern (low vegetation) storms (Fig. S7). Downwind cross sections are presented in Figs. 8a and 8b. The LSTA composites for both vegetation classes have similar characteristics around the initiation location, although the mean initiation gradient is 0.24 K (100 km) 21 weaker for the high NDVI subset. This is consistent with SM controls on convection being more important in climatologically drier (less vegetated) environments. For example, numerical experiments have found that vegetation exerts a ''homogenizing'' effect on domain SM due to increased moisture uptake from the root zone in wetter areas (Ivanov et al. 2010). However, the significant LSTA gradient  Table 1.  Cui et al. (2020) using laboratory measurements. It should also be noted that for a given SM induced gradient in sensible heat flux, denser (aerodynamically rougher) vegetation would produce a weaker LST gradient than a surface with more bare soil. Therefore we cannot conclude that the steeper LSTA gradient found in Fig. 8a is necessarily associated with a stronger sensible heat flux gradient. It is also worth noting that vegetation contrasts on the TP may have an effect on the timing of local convection (Babel et al. 2014), which will not be captured by the current analysis.

Discussion
This work reveals the influence of mesoscale SM heterogeneity on the initiation of strong convection over the TP. Using satellite observations we show a preference for strong convective development on the downwind side of dry surfaces, close to wetter areas, similar to results from Europe (T15) and West Africa (T11). The spatial relationship between surface gradients and convective initiation in these studies suggests that the preference for afternoon rain to occur over locally drier soils globally (Taylor et al. 2012 First, over the mountainous TP, there are known orographic effects on strong convection (e.g., Kurosaki and Kimura 2002;Barros et al. 2004;Fujinami et al. 2005). In general, orographic effects on convective initiation are considered to be stronger than land surface heterogeneity effects (Houze 2012). For this reason, unlike the Europe study, we chose to explicitly analyze the impact of topography on our result. Consistent with expectations, our observed impact of local SM variability on the development of strong convection is moderated by local topographic complexity, as demonstrated by our subset analysis. The downwind LSTA gradient for the low topographic complexity subset [21.23 K (100 km) 21 ] is nearly double the mean gradient for all cases [20.68 K (100 km) 21 ].
Having an average elevation of 4000 m and ''on plateau'' mountain ranges, the atmospheric dynamics on the TP are significantly more complex than the Sahel. This is highlighted by the high sensitivity of our signal to background wind speed. It must be acknowledged that our subset analysis also displays a sensitivity to the choice of reanalysis product. However, all three sets of results show a change in surface gradient length scale with increasing wind speed, from tens of kilometers for lighter winds to hundreds of kilometers for stronger winds. This suggests that strong wind speeds have the effect of suppressing the impact of small-scale heterogeneity on convective initiation (consistent with numerical simulations; e.g., Avissar and Schmidt 1998), but larger-scale gradients may still be important in these cases. This result is different to both previous studies as no sensitivity to wind speed was found for the Sahel, and no significant relationship between gradients and initiations for high wind speeds was found in Europe. There will likely be some uncertainty in the composites due to satellite product limitations. The TP is a challenging region of the world for remote sensing (e.g., Meng et al. 2018;Zeng et al. 2015). Although we have made special effort to retain only the highest quality observations by filtering the datasets (see section 2b), some errors may remain. This will include uncertainties in the measurements and potential misidentification of rain contaminated pixels. This work makes use of microwave LST observations that provide data for all-sky conditions, hence the measurements may be influenced by cloud cover. We observe a small sensitivity of our LSTA signal to variations in cloud cover, which will have increased the derived gradients. Also, our initiation locations are derived from hourly CTT data, whereas both T11 and T15 had 15 min observations available to them. The cloud field can evolve significantly in 1 h, particularly if storms are moving quickly. This will undoubtedly compromise our ability to accurately map the location where rapid growth in convection first occurs, effectively introducing some spatial noise. In addition the LSTA gradients may be weaker for Tibet due to the limited time window within the diurnal cycle of convective initiations considered in our analysis (limited by our observations of LST). The fraction of initiations occurring over complex topography (SD . 300 m) is smallest around 1400 LT (28% compared to the average 35% in our time window, Fig. S5b). This may suggest the impact of soil moisture (verse topography) on convective triggering is stronger earlier in the afternoon.

Conclusions
Our main conclusions may be summarized as follows: 1) As found in the Sahel, strong convection over the Tibetan Plateau is favored over negative (positive) land surface temperature (soil moisture) gradients.
New findings for Tibetan Plateau: 2) The signal is strongest for low topographic complexity, but still significant for local standard deviation in topography up to 300 m. 3) The atmosphere is only sensitive to local (tens of kilometers) soil moisture heterogeneity under calm synoptic conditions (light background wind speeds). 4) Increased vegetation moderates the surface temperature gradients; however, the soil moisture-convective initiation coupling signal is still prevalent in these cases.
Our results demonstrate that, even in the complex TP environment, mesoscale variations in SM still play an important role in storm development. Models will not only need to capture the SM-convection feedback, but also its dependence on the highly varying TP surface and synoptic conditions.