This study focused on producing flash flood hazard susceptibility maps (FFHSM) using frequency ratio (FR) and statistical index (SI) models in the Xiqu Gully (XQG) of Beijing, China. First, a total of 85 flash flood hazard locations (n = 85) were surveyed in the field and plotted using geographic information system (GIS) software. Based on the flash flood hazard locations, a flood hazard inventory map was built. Seventy percent (n = 60) of the flooding hazard locations were randomly selected for building the models. The remaining 30% (n = 25) of the flooded hazard locations were used for validation. Considering that the XQG used to be a coal mining area, coalmine caves and subsidence caused by coal mining exist in this catchment, as well as many ground fissures. Thus, this study took the subsidence risk level into consideration for FFHSM. The ten conditioning parameters were elevation, slope, curvature, land use, geology, soil texture, subsidence risk area, stream power index (SPI), topographic wetness index (TWI), and short-term heavy rain. This study also tested different classification schemes for the values for each conditional parameter and checked their impacts on the results. The accuracy of the FFHSM was validated using area under the curve (AUC) analysis. Classification accuracies were 86.61%, 83.35%, and 78.52% using frequency ratio (FR)-natural breaks, statistical index (SI)-natural breaks and FR-manual classification schemes, respectively. Associated prediction accuracies were 83.69%, 81.22%, and 74.23%, respectively. It was found that FR modeling using a natural breaks classification method was more appropriate for generating FFHSM for the Xiqu Gully.