The rapid decrease in Arctic sea ice cover and thickness not only has a linkage with extreme weather in the mid-latitudes but also brings more opportunities for Arctic shipping routes and polar resource exploration, both of which motivate us to further understand causes of sea-ice variations and to obtain more accurate estimates of sea-ice cover in the future. Here, a novel data-driven method, the causal effect networks algorithm, is applied to identify the direct precursors of September sea-ice extent covering the North Sea Route and Transpolar Sea Route at different lead times so that statistical models can be constructed for sea-ice prediction. The whole study area was also divided into two parts: the north region covered by multiyear ice and the south region covered by seasonal ice. The forecast models of September sea-ice extent in the whole study area (TSIE) and south region (SSIE) at lead times of 1–4 months can explain over 65％ and 79％ of the variances, respectively, but the forecast skill of sea-ice extent in the north region (NSIE) is limited at a lead time of 1 month. At lead times of 1–4 months, local sea-ice concentration and sea-ice thickness have a larger influence on September TSIE and SSIE than other teleconnection factors. When the lead time is more than 4 months, the surface meridional wind anomaly from north Europe in the preceding autumn or early winter is dominant for September TSIE variations but is comparable to thermodynamic factors for NSIE and SSIE. We suggest that this study provides a complementary approach for predicting regional sea ice and is helpful in evaluating and improving climate models.