This study projected the future rainfall (2046–2065 and 2081–2100) for the North China Plain (NCP) using two stochastic statistical downscaling models, the non-homogeneous hidden Markov model and the generalized linear model for daily climate time series, conditioned by the large-scale atmospheric predictors from six general circulation models for three emission scenarios (A1B, A2 and B1). The results indicated that the annual total rainfall, the extreme daily rainfall and the maximum length of consecutive wet/dry days would decline, while the number of annual rainfall days would slightly increase (correspondingly rainfall intensity would decrease) in the NCP, in comparison with the base period (1961–2010). Moreover, the summer monsoon rainfall, which accounted for 50–75 % of the total annual rainfalls in NCP, was projected to decrease in the latter half of twenty-first century. The spatial patterns of change showed generally north–south gradients with relatively larger magnitude decrease in the northern NCP and less decrease (or even slightly increase) in the southern NCP. This could result in decline of the annual runoff by −5.5 % (A1B), −3.3 % (A2) and −4.1 % (B1) for 2046–2065 and −5.3 % (A1B), −4.6 % (A2) and −1.9 % (B1) decrease for 2081–2100. These rainfall changes, combined with the warming temperature, could lead to drier catchment soil profiles and further reduce runoff potential, would hence provide valuable references for the water availability and related climate change adaption in the NCP.