Description: ABSTRACT Bankfull discharge is a key parameter in the context of river engineering and geomorphology, as an indicator of flood discharge capacity in alluvial rivers, and varying in response to the incoming flow and sediment regimes. Bankfull channel dimensions have significantly adjusted along the Lower Yellow River (LYR) due to recent channel degradation, caused by the operation of the Xiaolangdi Reservoir, which has led to longitudinal variability in cross-sectional bankfull discharges. Therefore, it is more representative to describe the flood discharge capacity of the LYR, using the concept of reach-averaged bankfull discharge. Previous simple mean methods to estimate reach-scale bankfull discharge cannot meet the condition of flow continuity or account for the effect of different spacing between two sections. In this study, a general method to calculate cross-sectional bankfull discharge using the simulated stage-discharge relation is firstly outlined briefly, and an integrated method is then proposed for estimating reach-scale bankfull discharge. The proposed method integrates a geometric mean based on the log-transformation with a weighted average based on the spacing between two consecutive sections, which avoids the shortcomings of previous methods. The post-flood reach-scale bankfull discharges in three different channel-pattern reaches of the LYR were estimated annually during the period from 1999 to 2010 using the proposed method, based on surveyed post-flood profiles at 91 sedimentation sections and the measured hydrological data at seven hydrometric sections. The calculated results indicate that: (i) the estimated reach-scale bankfull discharges can effectively represent the flood discharge capacity of different reaches, with their ranges of variation being less than those of typical cross-sectional bankfull discharges; and (ii) the magnitude of the reach-scale bankfull discharge in each reach can respond well to the accumulative effect of incoming flow and sediment conditions. Finally, empirical relationships for different reaches in the LYR were developed between the reach-scale bankfull discharge and the previous four-year average discharge and incoming sediment coefficient during flood seasons, with relatively high correlation coefficients between them being obtained, and the reach-scale bankfull discharges in different reaches predicted by the delayed response model were also presented for a comparison. These relations for the prediction of reach-scale bankfull discharges were validated using the cross-sectional profiles and hydrological data measured in 2011. This article is protected by copyright. All rights reserved.