Subsidence hazard has never been considered seriously until recent yews in Korea, although its socioeconomic impact on Korea becomes more and more enormous. There have been a few studies for the application of GIS analysis technique to the prediction of subsidence hazard. For GIS analysis, several factors, which are represented by coverage, are considered and selected for building a GIS model. Numerical method was used to quantify the importance of each factor in GIS model and the result from numerical modeling using FLAC was compared with that from previous research based on empirical methods. Analysis in 3-D needs more computer resources (i.e. memory). Therefore that in 2.5-D was considered to overcome the problem. Not only maximum vertical subsidence but also maximum horizontal strain and maximum slope have been considered for the assessment of subsidence hazard. The model can be easily modified for the purpose of applications in any subsidence area, especially cavern or abandoned mines under thick soil layer.