Understanding the Predictability within Convection-Allowing Ensemble Forecasts in East China: Meteorological Sensitivity, Forecast Error Growth and Associated Precipitation Uncertainties Across Spatial Scales

This study investigates the practical predictability of two simulated mesoscale convective systems (MCS1 and MCS2) within a state-of-the-art convection-allowing ensemble forecast system. The two MCSs are both controlled by the synoptic Meiyu-front but differ in mesoscale orographic forcing. An observation system simulation experiment (OSSE) setup is first built, which includes flow-dependent multiple-scale initial and lateral boundary perturbations and a 12 h 30-member ensemble forecast is thereby created. In combination with the difference total energy, the decorrelation scale and the ensemble sensitivity analysis, both forecast error evolution, precipitation uncertainties and meteorological sensitivity that describe the practical predictability are assessed. The results show large variabilities of precipitation forecasts among ensemble members, indicative of the practical predictability limit. The study of forecast error evolution shows that the error energy in the MCS1 region in which the convection is blocked by the Dabie Mountains exhibits a simultaneous peak pattern for all spatial scales at around 6 h due to strong moist convection. On the other hand, when large-scale flow plays a more important role, the forecast error energy in the MCS2 region exhibits a stepwise increase with increasing spatial scale. As a result of error energy growth, the precipitation uncertainties evolve from small scales and gradually transfer to larger scales, implying a strong relationship between error growth and precipitation across spatial scales, thus explaining the great precipitation variability within ensemble members. These results suggest the additional forcing brought by the Dabie Mountains could regulate the predictability of Meiyu-frontal convection, which calls for a targeted perturbation design in convection-allowing ensemble forecast systems with respect to different forcing mechanisms.


Introduction
Forecasting warm-season convective events is a consistent important and challenging issue during the most recent decades over the Yangtze and Huai River Basin (YHRB) [1,2]. In this area, in addition to the synoptic quasi-stationary Meiyu-front which provides large-scale ascent [3], the initiation and development of convective systems are also strongly influenced by steady surface features such as complex terrain (e.g., the Dabie Mountains, Huang mountain, and Mufu mountain) front is chosen ( Figure. 1b), one of which (MCS1) starts upstream of the Dabie Mountains and is influenced by the mountain-flow interaction (an example of how Dabie Mountains regulate the lowlevel flow is shown in Figure. 1c), while the other one (MCS2) locates over the downstream place of the Dabie Mountains. Since both of the two MCSs are governed by the large-scale Meiyu-front, we assume if the presence of Dabie Mountains could substantially impact the practical predictability of synoptic frontal forced precipitation, thus lead to different error growth dynamics and produces different precipitation uncertainties. To this end, three main questions are addressed: 1) what is the precipitation forecast sensitivity with respect to thermodynamic states at initial/previous times of a typical frontal precipitation system and that influenced by the Dabie Mountains; 2) considering multiple-scale errors induced from initial and lateral boundary condition, what is the evolution characteristics of forecast error? 3) how the forecast error growth influences the distribution of precipitation uncertainties across spatial scales and causing large precipitation variability within ensemble forecasts? The consideration of these issues provides insight into the forecast performance within convection-allowing ensemble forecasts system and potentially aids the optimal design for it especially when facing different types of convective events.
The outline of the article is as follows. The model and ensemble design are illustrated in section 2, along with the details of diagnostic approaches. Section 3 gives the results including case overview, ensemble performance, sensitive analysis, error growth, and associated precipitation uncertainty analysis. Section 4 provides the summary and discussion.

Model Configureuration
Following Zhuang et al. [12], the Advanced Research core of the Weather Research and Forecasting (WRF) Model version 3.7.1 [39] is used with an inner convection-allowing domain of 258 × 258 horizontal grid points at 3-km grid spacing, nested within an outer parameterizing domain of 180 × 180 grid points at 18-km grid spacing, and 41 terrain-following hydrostatic-pressure vertical levels are applied in all domains. The outer domain covers East China, while the inner domain targets the YHRB area, including Anhui and Jiangsu provinces and surrounding areas ( Figure. 1a).
The initial conditions (IC) and lateral boundary conditions (LBC) are derived from the U.S. National Center for Environmental Protection Global Forecast System (GFS) (0.5°) analysis data. The physical parameterization schemes include the WSM6 microphysics [40], Grell-3 cumulus [41] (which is only applied to the outer domain), YSU boundary layer [42], RRTM longwave radiation [43], and Goddard shortwave radiation schemes [44]. The inner domain is used as the analysis domain. These schemes were chosen for their good implements over the YHRB region [45,46].  Figure.1 a. The blue square marks the general range of synoptic low-level jet while the orange square marks the range of boundary layer low-level jet.

Ensemble Generation
To generate IC and LBC perturbations in a flow-dependent manner [47], a modified version of ensemble data assimilation and forecast system is constructed based on Wang et al. [48]. The data assimilation method is based on the Ensemble Square Root Filter approach proposed by [49]. For an OSSE, a true state is first created from a 36-h nature run initialized from the GFS analysis 24 h before the analysis time for each case with a 3-km resolution over the outer domain (TRUE). Thereafter, the mesoscale sounding datasets for the outer domain are created at randomly chosen locations within China's mainland and interpolated from TRUE with observation errors of 2.5 m/s, 1.2 K, and 0.005 kg/kg for wind velocity, temperature, and water vapor, respectively [50]. For the inner domain, TRUE is interpolated to seven real CINRAD-SA radars (see Figure. 1 in Zhuang et al. [12]), a Chinese version of Weather Surveillance Radar-1988 Doppler (WSR-88D) [51], with observation errors of 2 m/s and 5 dBZ for radial velocity and reflectivity [47].
For the outer domain, the initial ensemble is initialized by adding analysis perturbations to U, V, T, Qv in IC and LBC with a supersaturation correction for Qv. The perturbations are acquired from the first 30-member European Centre for Medium-Range Weather Forecasts global ensemble prediction products (0.5°, https://apps.ecmwf.int/datasets/) that well represent the synoptic-scale uncertainty. For the present case study, a 3-h spin-up ensemble is start at 12:00 UTC (2000 LST) 6 July 2013, along with the data assimilation of conventional observation datasets in the outer domain every 3-h over a 24-h period. The analysis perturbations at the seventh cycle are downscaled to the inner domain using the WRF ndown tool [52]. For the inner domain, the radar datasets (including reflectivity and radial velocity) are assimilated for 3-h with a 10-min interval, while the outer domain provides larger-scale LBC perturbations every 15-min during the assimilation process. Details of the data assimilation parameters applied can be found in Zhuang et al. [12]. Thereafter, a 30-member 12h ensemble forecast from 12:00 UTC 6 July 2013 to 20:00 UTC 7 July 2013 is performed with perturbed IC and LBC (CTRL).

Ensemble Sensitivity Analysis
The ensemble sensitivity analysis (ESA) proposed by Ancell and Hakim [53] is an efficient tool for estimating the meteorological forecast sensitivity to the thermodynamics variables in earlier times. Recently, studies have demonstrated that ESA could be used in convection-allowing ensemble forecast systems to estimate the sensitivity of non-gaussian distribution convective metrics such as precipitation and composite reflectivity to earlier thermodynamics variables [54,55]. In the present study, we use ESA to explore the meteorological sensitivity of Meiyu-frontal precipitation to the thermodynamic conditions at earlier forecast lead times. Details of information about ESA can be seen in Hill et al. [55].

Measurement of Forecast Error
Following Nielsen and Schumacher [19], the forecast error is quantified using a modified version of DTE. In addition to the typical DTE, a vertical function is applied to calculate the two-dimensional root mean vertically integrated DTE (RMDTE): where , , and are the differences of zonal wind, meridional wind, and temperature subtracted the ensemble mean from each ensemble member, respectively ( = 1004.9 JKg K and = 270 K). The subscripts i, j, k, t, and m represent the x-direction, y-direction, vertical level, forecast time, and ensemble member. To further investigate the scale-dependent error energy characteristics, we decompose RMDTE into three scale ranges: the small-scale (36 km ≥ wavelength), the transitionscale (120 km ≥ wavelength > 36 km), and the larger-scale (wavelength > 120 km). The small-scale represents the flow-dependent initial errors at the convective scale and can only be resolved in the inner domain; the transition-scale encompasses the smallest scale that can be resolved by the outer domain. To obtain variable fields at different scales, the discrete cosine transform (DCT) [27] method is used. DCT is similar to the fast Fourier transform and can avoid discontinuity problems at domain boundaries. The RMDTE at different scale ranges can be acquired in two steps: 1) decomposing , , and into different scale ranges , , and ; 2) computing RMDTE( ) using equation 1.

Measurement of Precipitation Uncertainties
A novelty of this work is the detailed investigation of scale-dependent characteristics between error energy growth and associated precipitation uncertainties. Surcel et al. [27] suggested a quantitative method to evaluate the scale-dependent features of precipitation uncertainties, namely, the decorrelation scale (λ0), which allows for displacements in the forecasts. In this method, for the scale λ ≤ λ0, there is no predictability of precipitation forecasts from ensemble members. For the scale λ ≥ λ0, the ensemble members are correlated, indicative of some predictability. The power ratio in terms of decorrelation scale is defined as: where ( ) refers to the variance of the precipitation field ( = 1, … , ) at scale λ, ( ) is the uncertainty value varying from 1/n to one, where a larger ( ) corresponds to higher uncertainty and lower precipitation forecast predictability at scale λ. Theoretically, a threshold value of ( ) = 1 at a given scale means the complete loss of predictability. In this study, the threshold is 0.9 to eliminate noise without introducing any significant bias [10,23,35,36].

Results
This section presents results on the practical predictability within convection-allowing ensembles of two MCSs along the same Meiyu-front over YHRB. We first introduce the general features of the events and analyze the ensemble performance. In a step forward, we present the scaledependent forecast error growth and explore how error energy evolution could impact the diagnostic precipitation variable, which could explain the precipitation variability within ensemble forecasts. Finally, the forecast sensitivities for the two convective systems are analyzed and compared using subjective analysis in combination with the ESA method.

Case Overview and Ensemble Performance
As frontal forced convections, the synoptic conditions combined to produce ingredients conducive to the formation and development of the convective events. During the study period, there is a southwest-northeast-oriented quasi-stationary synoptic cold front (indicated by the 850-hPa wind shear line in Figure. 2) dominating convections in the analysis domain which is known as the Meiyufront. A strong southwest low-level (850 hPa ≈1500 m height) jet is evident over southern Anhui, advecting moisture northward toward the frontal areas. During the mid-night (Figureure. 2b-d, gi), the low-level jet gradually strengthens and sustains the development of convections (Figureure. 3a). Therefore, the synoptic and subsynoptic forcing benefit the convections we focus on. Besides, for the Dabie Mountains region (red boxes), the 925-hPa boundary layer low-level (925 hPa ≈ 700 m height) jet [56] (Figure 2f-j) could interact with the Dabie Mountains [6] and impact the precipitation upstream the mountain region.  Figure. 4, it is clear that the area-averaged values for the MCS2 region have a more evident semidiurnal cycle than that for the MCS1 region. As a short summary, the two MCSs have similar large-scale but different mesoscale forcing mechanisms, leading to different precipitation distribution and semidiurnal features. The 30-member convection-allowing ensemble forecasts (CTRL) show a great divergence of precipitation within ensemble members ( Figure. 4), indicative of the practical predictability limit [23,57] induced by errors originated from initial and lateral boundary conditions. Most of the members in CTRL produce rainfall of similar magnitudes to those in TRUE but have errors in location and organizational structure to different degrees (Not shown). Figureure 4 shows the time evolution of area-averaged 30 min accumulated precipitation for CTRL, MEAN (the ensemble mean forecast) and TRUE. Most members reflect the semidiurnal feature of both MCS1 and MCS2. Specifically, almost all members generate a wet bias between 17:00-21:00 UTC (5-9 h) for the MCS1 region and a dry bias during 17:00-22:00 UTC (5-10 h) for the MCS2 region, highlighting a challenge for optimal ensemble design over YHRB region, especially for events controlled by different forcing mechanisms. Besides, the variability within CTRL is evident and the large values are generally corresponding to the peak pattern of precipitation (e.g., 4 h for the MCS1 region and 9 h for the MCS2 region), which potentially implies the importance of moist convection on predictability.

Spatiotemporal Characteristics of Scale-Dependent Forecast Error Growth
In a convection-allowing NWP, the IC errors could nonlinearly grow at small scales within precipitation regions and break the predictability afterward, which leads to the large precipitation variability within CTRL ( Figure. 4). In this subsection, we tend to understand the error growth dynamics through the scale-dependent analysis of spatiotemporal characteristics of forecast error energy evolution. Figureure 5 shows the distribution of RMDTE at different valid times across different spatial scales. As flow-dependent perturbations derived from a multiscale ensemble data assimilation and forecasting system, the initial small-scale RMDTE clearly corresponds to the precipitation areas south of the Meiyu-front ( Figure. 5a). During the next few hours, the small-scale RMDTE increases in magnitude, spreads out and keeps following the convection regions along the Meiyu-front zone. For the transition-scale ( Figure. 5b), the initial RMDTE exhibits a different pattern compared with that in the small-scale but still partly corresponds to the precipitation region. After 2 h (the second column in Figure. 5b), the transition-scale RMDTE gradually evolves into a similar pattern as the small-scale RMDTE, albeit with a broader area, illustrating an upscale transfer of error [16,58]. For the largerscale, the initial pattern of RMDTE ( Figure. 5a3) generally follows the synoptic shear line ( Figure. 1b), depicting flow-dependent uncertainties at larger scales. Similar to the transition-scale, we also find the larger-scale RMDTE gradually evolves into the precipitation zone, but with a slower speed. This feature further confirms the upscale transfer of errors is clear in this study.
Additionally, a comparison between the two MCS regions shows the RMDTEs at small-and transition-scale increase simultaneously, which is due to error growth caused by moist convection. However, for the larger-scale, the increase of error magnitude is faster in the MCS1 region ( Figure. 9d3) than that in the MCS2 region. It's speculated this is due to the larger-scale LBC errors imposed from the western boundary which directly promotes the forecast error through regulating the westerly flow at the larger-scale and impact the practical predictability afterward. This result indirectly suggests the importance of LBC perturbations for the design of convection-allowing ensembles. To further illustrate the error evolution characteristics for the two MCSs regions, we next perform a time series analysis of RMDTE ( Figure. 6). In the MCS1 region, the total RMDTE (black curve in Figure. 6a) generally shows two peaks at around 1 and 7 h, while for the MCS2 region the total RMDTE (black curve in Figure. 6b) exhibits a moderate trend after the rapid increase before 4 h. The difference of error energy growth between the two MCS regions suggests different predictability for these two events and thereby imply the complexity for forecasting Meiyu-front forced convective events.
We further decompose the total RMDTE into different scale ranges. The time evolution of RMDTE at different scales clearly shows the upscale transfer processes for both MCS regions. In general, the small-scale RMDTE could rapidly increase and reaches a saturation (1 h for the MCS1 region and 3 h for the MCS2 region) (blue curves in Figure. 6) under the influence of moist instability and latent heat release [16][17][18][59][60][61][62], followed with lagging increase and saturation of errors at the transition-and larger-scale, respectively. This result implies that the three-stage error growth conceptual model [16,58] is generally applicable for convective events controlled by Meiyu-front.
Though followed with a similar error growth model, it is important to notice that the scaledependent evolution trends are also different between the MCS1 and MCS2 regions. The difference can be linked to the meteorological nature (forcing mechanisms) of the two MCSs. Specifically, for the MCS2 region in which the convective system is controlled by large-scale flow and moving eastward along the frontal zone, the evolution of scale-dependent errors strictly follows the stepwise three-stage error growth mode ( Figure. 6b) (i.e., the errors increasingly transfer into larger scales at lagging times): first, the small-scale errors rapidly increase and reach saturation at around 3 h, after that, the transition-scale errors keep another two-hour increasing and saturate at 5.5 h, while for the larger-scale, the errors could keep a steadier increasing throughout the 12-h range. This growth mode can be explained by the displacement of individual cells driven by large-scale forced ascent [16,63,64] around the Meiyu-front. Inconsistent with MCS2, in the MCS1 region, the small-scale RMDTE performs a secondary increase within 4-6 h after the first saturation at 1 h ( Figure. 6a). Concurrently, the RMDTEs at transition-and larger-scale practically exhibit a simultaneously mark-up during 4-7 h which is consistent with the pronounced strengthening of precipitation ( Figure. 4a). According to previous work, Nielsen and Schumacher [19] proposed a theoretical assumption that the trend of error growth can be generally controlled by the large-scale flow or quasi-stationary convection diurnal cycle, while the magnitude (peak or other big fluctuations) can be modified by strong moist convection. We think this theory can be slightly extended to explain the scale-dependent error growth for the MCSs regions in this study. In other words, if the large-scale flow is more important, the error growth mode in the MCS2 region is common. For the MCS1 region in which the convection is strongly restricted by the Dabie Mountains, the forecast errors at small scales are able to grow consistently under the control of moist convection, impact larger scales via upscale transfer, and muddle the stepwise growth trend ( Figure. 6b). Hohenegger et al. [59] also suggested a theory if the errors could propagate against the mean flow in the unstable area, the errors may reach a sufficient amplitude and seriously disrupt predictability, which corroborates the explanation for the MCS1 region. In conclusion, though under similar large-scale forcing mechanisms (Meiyu-front and southwest low-level jet), the additional mesoscale orographic forcing exerted by the Dabie Mountains could largely influence the atmospheric predictability of frontal precipitation. This explains the recently found challenge of convection-allowing NWP in forecasting convective events adjunct to the Dabie Mountain during Meiyu seasons [65], especially for that near the LBC of the model domain.

Relationship Between Error Energy Growth And Precipitation Forecasts
A novelty of this study is to build a quantitatively and scale-dependent relationship between error energy growth and associated precipitation forecasts. Since the evolution of forecast error energy is clearly linked to the precipitation forcing mechanisms and varies between the two MCS regions, we further inspect the corresponding characteristics of precipitation uncertainties, which directly relate to the precipitation variability found in section 3.1 (Figure. 4).
Figureure 7a,b shows the power ratio (abbreviated as R) as a function of spatial scale for the two MCS regions at different lead times, which quantitatively illustrates the evolution of precipitation uncertainties in terms of spatial scales regardless of the selection of a certain threshold. It is clear that R increases at small scales and gradually transfer to larger scales with lead time, corresponding to the upscale transfer processes that revealed in the error energy growth processes ( Figure. 6). For both regions, with the rapid increase of small-scale error energy ( Figure. 6), R first increases at scales within 20 km at 0.5 h (black curves in Figure. 7). This finding also indicates that the precipitation could only be monotonously influenced from small to larger scales even multiple-scale information is contained in the initial perturbations. After the first saturation of small-scale RMDTE, R exceeds the threshold 0.9 at scales within 20 km (blue curves in Figure. 7), implying a full loss of predictability at the meso-γ-scale. Meanwhile, along with the upscale transfer of errors, R reflects a slower increase at larger scales at 2 h. After that, with the constant upscale error transfer from small to larger scales of RMDTE, R exhibit increases at all scales within the meso-β-scale, revealing the predictability loss of precipitation forecasts at relevant scales.
Considering the dissimilarity of RMDTE growth between the two MCS regions ( Figure. 6), we also find corresponding differences revealed on the precipitation uncertainties. The decorrelation scale λ0, which is further calculated based on R, is displayed in Figure. 7 (c, d). For both MCS regions, the precipitation forecasts of ensemble members at scales less than λ0 are fully uncorrelated, indicating totally loss of predictability (i.e., untrusted forecasting results). For both MCS regions, the precipitation forecasts in CTRL are more or less correlated at all spatial scales at 0.5 h and start at 6 km at 1 h. After that, for the MCS2 region, the consistent increase of area-averaged small-scale RMDTE (1-2.5 h) leads to an earlier growth of λ0 ( Figure. 7b) than that in the MCS1 region. Next, along with the subsequent saturation of small-scale RMDTE, the growth rate of λ0 gradually becomes slower in the MCS2 region than that in the MCS1 region. This feature is owing to the secondary increase of small-scale RMDTE in the MCS1 (Figure. 6a). Besides, in the MCS1 region, λ0 reaches a peak at around 7 h, which is also consistent with the peak phase of RMDTE showed in Figure. 6a. The fit lines computed with the least square regression also reveal the strong relationship between error energy growth and precipitation uncertainties.
In conclusion, the precipitation uncertainties first arise from small scales and gradually transfer to larger scales even impacted by multiscale initial errors. Specifically, the argument of small-scale error energy magnitude first improves R start from small scales and reaches a maximum (predictability loss) with the saturation of small-scale errors. Then with the upscale transfer of smallscale errors, R gradually increases from smaller to larger scales with the peak shifting to larger scales simultaneously.

Meteorological Sensitivities
Since the two MCSs are found to be different in both error energy growth and precipitation uncertainties, what are the detail meteorological forecast sensitivities that are responsible for the dissimilarity in predictability, especially for the particular time when forecast bifurcations occur? Understanding this issue could help to identify the sensitivity region and variables which in turn improve the perturbation design for the ensemble forecast system.
To address this issue, a subjective analysis of the ensemble members compared with TRUE is further performed with the determination of the spatial phasing and temporal phasing of the precipitation forecasts. Members 11 and 3 are then selected as representations of forecast success (GOOD) and forecast failure (POOR) for further analysis. Although there are very small differences of IC between these two representative members (e.g., T < 0.5 K; Wind-925 hPa < 0.5m/s) and very subtle difference in precipitation intensity is therefore caused within 0.5 h (compare the red curves with the blue curves in Figure. 4), the precipitation location and intensity vary widely over the next few hours (compare Figure. 3a2-e2 with Figure.3 a3-e3). The good member, member 11 (GOOD), could generally reproduce the intensity and location of the two precipitation bands in the MCS1 region. The central location and intensity for MCS2 can also be well represented by member 11 ( Figure. 3a2-e2). In contrast, in the MCS1 region, the convective systems for member 3 (POOR) only develops in the western lee of the Dabie Mountains at early 2 h lead hour ( Figure.   Next, two valid lead times when great forecast divergence initially occurs (regard as an approximately beginning of forecast bifurcation, Figure. 4) in the MCS1 (2 h) and MCS2 (7 h) regions are selected, respectively. Figureure 8 illustrates the precipitation and environmental conditions at 1.5 h forecast lead hour (0.5 h previous to the divergence time) for the MCS1 region. The nearboundary convection is barely decayed in POOR ( Figure. 8b) but evident in GOOD at this time, causing large precipitation difference ( Figure. 8c) and reflecting the early expansion of ensemble variability (Figure.4a). The second row in Figure. 8 shows that GOOD has a stronger 925-hPa southerly wind south of the MCS1 region than that in POOR (marked by the arrows in Figure. 8f), benefiting the convergence difference upstream the Dabie Mountains (marked by the red circle in Figure. 8f). This finding reveals that the nighttime strengthened boundary layer low-level jet ( Figures.  2f-j) is important for MCS1 by bringing in more warm, moist air (larger ) and causing greater upward vertical motion in GOOD (Figure. 8g). We also find the convergence can be strengthened by the substantial stronger 925-hPa westerly wind that appears in POOR ( Figure. 8e) than GOOD ( Figure. 8d), resulting in the obvious precipitation location error in POOR (Figure. 8b). As a short summary, the location and intensity of precipitation in the MCS1 region are closely related to the boundary layer low-level jet and the westerly wind upstream of the Dabie Mountains. For the MCS2 region, evident expansion of ensemble variability is found to start at around 7 h ( Figure. 4b). Figureure 9 displays the precipitation and environmental conditions for GOOD and POOR. The 1-h accumulated precipitation shows POOR is gradually becoming disorganized ( Figure. 9b) and has large location and intensity errors ( Figure. 9c) compared with GOOD. At the 850-hPa, GOOD has a stronger 850-hPa southwest wind than POOR ( Figure. 9f), leading to more moist and warm air ( Figure. 9g) from the Bay of Bengal and a local convergence difference (red circle in Figure. 9f). Besides, the stronger westerly wind in GOOD is associated with warmer air that POOR at around 118.5° (Fig. 9i), causing a greater upward motion accordingly. These results reveal that the accurate forecasting of the strength of the synoptic low-level jet ( Figures. 2a-e) is responsible for the precipitation forecast in the MCS2 region. Note that the simulation of MCS1 could be greatly impacted by the western lateral boundary, but we mainly focus on the meteorological sensitivities themselves not that lead to the sensitivities. On the other hand, in the next sections, we care about the total error growth evolution and the associated impacts on precipitation which are comprehensive results of both initial and lateral boundary errors. Thus, the sensitivity and relative weight of errors from each source to forecast is beyond the scope and will not influence the reliability of our analysis.

Ensemble Sensitivity Analysis
In this subsection, ESA is introduced to quantitatively estimate the meteorological forecast sensitivity as a supplement of the subjective analysis in section 3.4 and shed more light on the understanding of practical predictability of the Meiyu-frontal events. The area-averaged 1-h precipitation in each MCS region is selected as the response metric. Figureure 10 shows the ESA results in the MCS1 region valid at 14:00 UTC (2 h) to 925-hPa horizontal wind field at the beginning hour (12:00 UTC). The precipitation forecast is negatively sensitive to the 925-hPa ( Figure. 10a) zonal wind, indicating the weakening of 925-hPa westerly flow upstream the Dabie Mountains produces higher simulated precipitation in the response region. This finding also indirectly implies the Dabie Mountains can strengthen the upstream precipitation by the physical block of the westerly flow. We also find a dipolar distribution of precipitation forecast sensitivity to meridional wind ( Figure. 10b) at the lower-left quarter of the domain, implying a higher southerly wind speed southern of response area and a decrease of southerly wind within the response area could produce a larger precipitation amount. In short, the boundary layer low-level jet is responsible for the accuracy precipitation forecast in the MCS1 region, consistent with the subjective analysis in Section 3.4 ( Figure. 8). For the MCS2 region, the precipitation forecast is positively sensitive to the upstream 850-hPa zonal wind ( Figure. 11a) and meridional wind ( Figure. 11b), indicating the nighttime strengthening of the southwesterly low-level jet is clearly responsible for the simulated precipitation amount. Therefore, the ESA results confirm the subjective GOOD versus POOR analysis shown in section 3.4 in a more objective aspect. Improvements of the model physical processes for the nighttime strengthened ( Figure. 2) southwest low-level jet, boundary layer low-level jet as well as its interaction with the Dabie Mountains is able to promote the practical predictability of Meiyu-front forced convective events. In the meantime, since the low-level dynamic conditions are crucial to the frontal convection present in our study, targeted perturbations can be added to the low-level jet area to improve the ensemble spread performance. Figure 11. Sensitivity of hourly precipitation in the MCS2 region (red rectangle) at 19:00 UTC 6 July (7 h) to (a) 850-hPa zonal wind (shaded, mm (m s −1 ) −1 ) and (b) 850-hPa meridional wind (shaded, mm (m s −1 ) −1 ) at 17:00 UTC (5 h). Black shading indicates regions where the sensitivity is statistically significant at the 95% confidence level. The black vectors are 850-hPa ensemble mean wind fields.

Conclusion and Discussion
The mechanism of warm-season convective events dominated by the Meiyu-front in east China has been well studied during the last three decades [2,3,6,29,32,33,65], whereas the relevant researches on predictability are still limited [1], especially on the small scales in which the errors could rapidly growth nonlinearly across spatial scales and destroy the predictability. In this study, the practical predictability of two simulated mesoscale convective systems (MCS1 and MCS2) along the Meiyu-front is investigated using a state-of-the-art convection-allowing ensemble forecast system. The two events are distinguished by the additional orographic forcing presented in MCS1. Then on this basis, both the ensemble performance, scale-dependent forecast error growth, precipitation uncertainties and forecast sensitivity are evaluated.
An OSSE setup based on EnKF is first built that includes flow-dependent multiple-scale initial and lateral boundary perturbations. Therefore, a 12-h 30-member convection-allowing ensemble forecast from 12:00 UTC to 20:00 UTC 6 July 2013 is created. The forecast error growth dynamics are evaluated by investigating RMDTE [19] at different spatial scales and the associated precipitation uncertainties caused by forecast error growth are further assessed using the decorrelation scale method [27,28]. The meteorological forecast sensitivity is evaluated using a subjective good versus poor member comparison [15,23] with an objective ESA method [23,55].
The 30-member forecasts show a large precipitation difference for the two mesoscale convective systems along the Meiyu-front, implying a chaotic divergence among the members. The forecast error growth which leads to the large variability is then examined. Results show that the error energy evolution in the MCS1 region in which the convection is influenced by the Dabie Mountains shows a simultaneous peak pattern for all spatial scales at around 6 h with the presence of strong moisture dynamic. On the other hand, it is showing a stepwise increase with increasing spatial scale of forecast error in the MCS2 region where the precipitation is purely controlled by the large-scale flow. This dissimilarity is strongly linked to the forcing mechanisms of the convection. With the evolution of forecast error energy, the precipitation uncertainties develop and change correspondingly. Unlike the error energy, the precipitation uncertainties can only evolve from small scales and gradually transfer to larger scales, explaining the great precipitation variability within ensemble members. For different convective systems, this feature is also different, indicating different practical predictability limits.
Further investigations on the meteorological sensitivity demonstrate that the location and intensity of precipitation in the MCS1 region are closely related to the nighttime strengthened boundary layer low-level jet and its interaction with the Dabie Mountains. While the strength of synoptic southwest low-level jet is responsible for the accuracy of precipitation in the MCS2 region.
To conclude, these results suggest the predictability of Meiyu-front events can be impacted by the Dabie Mountains which locate over the low-level southwesterly winds. Therefore, the perturbation design for a convection-allowing ensemble forecast system covering the YHRB should be treated accordingly with respect to different forcing mechanisms.
In this work, the general features of meteorological sensitivity and error growth for Meiyufrontal forcing events are revealed, and the impact of mesoscale orographic (the Dabie Mountains) forcing on such cases is preliminarily discussed. In a step forward, more on the impacts of the Dabie Mountains on practical predictability will be investigated through terrain modification experiments [66] and scale-based error sensitivity experiments [67]. In addition, it's noted that the ensemble forecast for MCS1 presents a generally wet bias, while in the MCS2 region there is a dry bias ( Figure.  4). This feature may link to the diurnal nature of convection in this region. The influence of the Dabie Mountains and their associated mountain-valley circulation on the MCSs may lead to the differences. Thus, the diurnal nature of the MCSs and relevant predictabilities will also be studied in our next works.
Author Contributions: conceptualization, X. Z., N. W. and J. M.; methodology, X. Z. and N. W.; formal analysis, X. Z. and N. W.; data curation, X. Z. and Y. X.; investigation, X. Z., N. W. and Y. X.; visualization, X. Z.; writingreview and editing, N. W.; funding acquisition, J. M.; supervision, J. M. All authors have read and agreed to the published version of the manuscript.