Land surface subsidence is an environmental hazard resulting from the extraction of underground resources. In underground mining, when mineral materials are extracted deep within the ground, the emptying or caving of the mined spaces leads to vertical displacement of the ground, known as subsidence. This subsidence can extend to the surface as troughs subsidence, as the movement and deformation of the hanging-wall rocks of the mining stope propagate upwards. Accurately predicting subsidence is crucial for estimating damage and protecting surface buildings and structures in mining areas. Therefore, developing a model that considers all relevant parameters for subsidence estimation is essential. In this article, we discuss the prediction of land subsidence caused by the caving of a stop's roof, focusing on coal mining using the longwall method. We consider a total of 11 parameters related to coal mining, including mining thickness and depth (related to the deposit), as well as density, cohesion, internal friction angle, elasticity modulus, bulk modulus, shear modulus, Poisson's ratio, uniaxial compressive strength, and tensile strength (related to the overburden). We utilize information collected from 14 coal mines regarding mining and subsidence to achieve this. We then explore the prediction of subsidence caused by mining using the gene expression programming (GEP) algorithm, optimized through a combination of the artificial bee colony (ABC) and ant lion optimizer (ALO) algorithms. Modeling results demonstrate that combining the GEP algorithm with optimization based on the ABC algorithm yields the best subsidence prediction, achieving a correlation coefficient of 0.96. Furthermore, sensitivity analysis reveals that mining depth and density have the greatest and least effects, respectively, on land surface subsidence resulting from coal mining using the longwall method.