Ice gouging problem is a significant challenge threatening the integrity of subsea pipelines in the Arctic (e.g., Beaufort Sea) and even non-Arctic (e.g., Caspian Sea) offshore regions. Determining the seabed response to ice scour through the subgouge soil deformations and the keel reaction forces are important aspects for a safe and cost-effective design. In this study, the subgouge soil deformations and the keel reaction forces were simulated by the extreme learning machine (ELM) for the first time. Nine ELM models (ELM 1–ELM 9) were developed using the key parameters governing the ice–seabed interaction. The number of neurons in the hidden layer was optimized and the best activation function for the ELM network was identified. The premium ELM model, resulting in the lowest level of inaccuracy and complexity and the highest level of correlation with experimental values was identified by performing a sensitivity analysis. The gouge depth ratio and the shear strength of the seabed soil were found to be the most influential input parameters affecting the subgouge soil deformations and the keel reaction forces. A set of the ELM-based equations were proposed to approximate the ice gouging parameters. The uncertainty analysis showed that the premium ELM model slightly underestimated the subgouge soil deformation.