Although land surface influences atmospheric processes significantly, insufficient studies have been conducted on the ensemble forecasts using the breeding of growing modes (BGM) with perturbed land surface variables. To investigate the practicability of perturbed land variables for ensemble forecasting, we used the ARWv3 mesoscale model to generate ensembles for an event of 24 h heavy rainfall with perturbed atmospheric and land variables by the BGM method. Results show that both atmospheric and land variables can generate initial perturbations with BGM, except that they differ in time and saturation characteristics, e.g., saturation is generally achieved in approximately 30 h with a growth rate of ~1.30 for atmospheric variables versus 102 h and growth rate of 1.02 for land variables. With the increase in precipitation, the importance of the perturbations of land variables also increases as compared to those of atmospheric variables. Moreover, the influence of the perturbations of land variables on simulated precipitation is still relatively large, although smaller than that of atmospheric variables, e.g., the spreads of perturbed atmospheric and land subsets were 7.3 and 3.8 mm, respectively. The benefits of perturbed initialisation can also be observed in terms of probability forecast. All findings indicate that the BGM method with perturbed land variables has the potential to ensemble forecasts for precipitation.