Significant coal reserves in New South Wales lie under areas with steep topography such as valley slopes and cliffs. Subsidence predictions are difficult to make for these locations as existing empirical databases cannot cover all permutations of topography and mining geometry. Two dimensional numerical model alternatives are of limited applicability. Consequently a novel three-dimensional numerical model is developed including the detailed panel layout, realistic rock mass properties including the full stratigraphy, goaf behaviour and actual surface topography. Model results show the same trends as subsidence measurements from Baal Bone Colliery in the Western Coalfield of New South Wales.