Understanding and forecasting the spatial and temporal distributions of extreme precipitation over urban areas is crucial for effective planning and mitigation efforts. However, this task remains challenging as accurate forecasting depends on properly representing urban surfacees and their interactions with the planetary boundary layer (PBL). We examined the hindcast of an extreme precipitation event over Beijing on 21–22 July 2012. The primary focus was assessing its sensitivity to two widely used PBL parameterizations (MYJ and YSU), two urban parameterizations (SLUCM and BEP_BEM), and two different land-use and land-cover (LULC) datasets.Sensitivity experiments were initialized at different times to explore the model dependence on initial conditions. The analyses were conducted over three selected regions: the entire model domain covering the Beijing metropolitan area, an upwind region of Beijing, and the entire urban area of Beijing. The results show that the MYJ PBL scheme performs better than the YSU PBL scheme in capturing near-surface air temperature as well as the location and timing of the heaviest precipitation. The variability in simulated precipitation among the chosen PBL schemes is lower compared to that among different time of initializations. The LULC impacted the spatial distribution of precipitation but its effect on the amount of precipitation was minimal. Overall, using a combination of the MYJ PBL scheme, SLUCM urban parameterization, and locally-enhanced Beijing LULC, and initializing the model simulations at 0000 UTC July 20, 2012, demonstrated superior performance in capturing precipitation levels, despite some spatial discrepancies in the precipitation distribution. The performance of BEP_BEM urban parameterization is similar to SLUCM across various factors such as average rain rate, maximum rain rate, and rain volume. These findings offer valuable insights towards better simulations of extreme precipitation and flooding in rapidly urbanizing areas such as Beijing.