The accuracy prediction for the crop yield is conducive to the food security in regions and/or nations. To some extent, the prediction model for crop yields combining the crop mechanism model with statistical regression model (SRM) can improve the timeliness and robustness of the final yield prediction. In this study, the accumulated biomass (AB) simulated by the Agricultural Production Systems sIMulator (APSIM) model and multiple climate indices (e.g., climate suitability indices and extreme climate indices) were incorporated into SRM to predict the wheat yield in the North China Plain (NCP). The results showed that the prediction model based on the random forest (RF) algorithm outperformed the prediction models using other regression algorithms. The prediction for the wheat yield at SM (the period from the start of grain filling to the milky stage) based on RF can obtain a higher accuracy (r = 0.86, RMSE = 683 kg ha⁻¹ and MAE = 498 kg ha⁻¹). With the progression of wheat growth, the performances of yield prediction models improved gradually. The prediction of yield at FS (the period from flowering to the start of grain filling) can achieve higher precision and a longer lead time, which can be viewed as the optimum period providing the decent performance of the yield prediction and about one month’s lead time. In addition, the precision of the predicted yield for the irrigated sites was higher than that for the rainfed sites. The APSIM-simulated AB had an importance of above 30% for the last three prediction events, including FIF event (the period from floral initiation to flowering), FS event (the period from flowering to the start of grain filling) and SM event (the period from the start of grain filling to the milky stage), which ranked first in the prediction model. The climate suitability indices, with a higher rank for every prediction event, played an important role in the prediction model. The winter wheat yield in the NCP was seriously affected by the low temperature events before flowering, the high temperature events after flowering and water stress. We hope that the prediction model can be used to develop adaptation strategies to mitigate the negative effects of climate change on crop productivity and provide the data support for food security.

Global climate change has had a significant impact on crop production and agricultural water use. Investigating different future climate scenarios and their possible impacts on crop production and water consumption is critical for proposing effective responses to climate change. In this study, based on daily downscaled climate data from 22 Global Climate Models (GCMs) provided by Coupled Model Intercomparison Project Phase 6 (CMIP6), we applied the well-validated Agricultural Production Systems sIMulator (APSIM) to simulate crop phenology, yield, and water use of the rice–wheat rotation at four representative stations (including Hefei and Shouxian stations in Anhui province and Kunshan and Xuzhou stations in Jiangsu province) across the Huang-Huai-Hai Plain, China during the 2041–2070 period (2050s) under four Shared Socioeconomic Pathways (i.e., SSP126, SSP245, SSP370, and SSP585). The results showed a significant increase in annual mean temperature (Temp) and solar radiation (Rad), and annual total precipitation (Prec) at four investigated stations, except Rad under SSP370. Climate change mainly leads to a consistent advance in wheat phenology, but inconsistent trends in rice phenology across four stations. Moreover, the reproductive growth period (RGP) of wheat was prolonged while that of rice was shorted at three of four stations. Both rice and wheat yields were negatively correlated with Temp, but positively correlated with Rad, Prec, and CO₂ concentration ([CO2]). However, crop ET was positively correlated with Rad, but negatively correlated with [CO2], as elevated [CO2] decreased stomatal conductance. Moreover, the water use efficiency (WUE) of rice and wheat was negatively correlated with Temp, but positively correlated with [CO2]. Overall, our study indicated that the change in Temp, Rad, Prec, and [CO2] have different impacts on different crops and at different stations. Therefore, in the impact assessment for climate change, it is necessary to explore and analyze different crops in different regions. Additionally, our study helps to improve understanding of the impacts of climate change on crop production and water consumption and provides data support for the sustainable development of agriculture.

Climate change has and will continue to exert significant effects on social economy, natural environment, and human life. Research on the climatic suitability of crops is critical for mitigating and adapting to the negative impacts of climate change on crop production. In the study, we developed the climate suitability model of maize and investigated the climate suitability of summer maize during the base period (1981–2010) and two future periods of 2031–2060 (2040s) and 2071–2100 (2080s) in the North China Plain (NCP) based on BCC-CSM2-MR model (BCC) from the Coupled Model Comparison Program (CMIP6) under two Shared Socioeconomic Pathways (SSP) 245 and SSP585. The phenological shift of maize under future climate scenarios was simulated by the Agricultural Production Systems Simulator (APSIM). The results showed that the root mean square errors (RMSE) between observations and projections for sunshine suitability (S_S), temperature suitability (S_T), precipitation suitability (S_P), and integrated climate suitability (S_Z) during the whole growth period were 0.069, 0.072, 0.057, and 0.040, respectively. Overall, the BCC projections for climate suitability were in suitable consistency with the observations in the NCP. During 1981–2010, the S_P, S_T, and S_Z were high in the north of the NCP and low in the south. The S_P, S_T, and S_Z showed a downward trend under all the future climate scenarios in most areas of NCP while the S_S increased. Therein, the change range of S_P and S_S was 0–0.1 under all the future climate scenarios. The S_T declined by 0.1–0.2 in the future except for the decrease of more than 0.3 under the SSP585 scenario in the 2080s. The decrease in S_Z in the 2040s and 2080s under both SSP scenarios varied from 0 to 0.2. Moreover, the optimum area decreases greatly under future scenarios while the suitable area increases significantly. Adjusting sowing data (SD) would have essential impacts on climate suitability. To some extent, delaying SD was beneficial to improve the climate suitability of summer maize in the NCP, especially under the SSP585 scenario in the 2080s. Our findings can not only provide data support for summer maize production to adapt to climate change but also help to propose agricultural management measures to cope with future climate change.

Extreme temperature events are becoming more frequent due to global warming, and have critical effects on natural ecosystems, social and economic spheres, human production and life. Predicting changes in temperature extremes and trends under future climate scenarios helps to assess the impact of climate change accurately. Based on climate observations from 54 meteorological stations in the North China Plain and the projection data from seven general circulation models (GCMs) from the Coupled Model Intercomparison Project phase 6 (CMIP6), this paper researches nine representative extreme temperature indices under four typical climate scenarios. The aim is to reveal the temporal and spatial variations in extreme temperature indices in the North China Plain during the past (1971–2010) and the future (2061–2100). The results show that: using a support vector machine (SVM) to perform regression analysis on the multi-GCMs prediction results, the root mean square error (RMSE) and relative standard deviation (RSD) of the multi-model ensemble simulations obtained by the SVM method are lower than those of the arithmetic mean method and can better match the trend of the historical extreme temperature index; the extreme high temperature index is predicted to show a significant upward trend in the future, while the extreme low temperature index will decrease significantly; and there are significant spatial differences in the extreme temperature index in both historical and future periods, with the extreme temperature index under the high radiation forcing scenario (SSP585) showing the most considerable variation and the most significant spatial differences.