This paper discusses the prediction of subsidence over partial extraction systems in thick coal seams. A model based on empirical techniques using influence functions is developed to calculate the surface ground settlements above extraction panels. Independent of the influence function, two field parameters, the subsidence factor and the angle of draw (limit angle) are required to define the model. Both these parameters are best observed from critical or super critical subsidence areas. However in thick seams, extraction widths are typically sub critical and therefore it was necessary to refine an estimation process for sub critical conditions. Partial extraction systems in NZ are characterised by variability with respect to extraction thickness and depth. As extraction proceeds along a panel, conditions deteriorate once the critical span width is exceeded and therefore percentage extraction decreases. A further variable was introduced to the modelling which allows for varying extraction sequences. Two comparisons with field observations, one over a partial extraction panel, the other over abandoned workings, showed the method gave good agreement to measured subsidence values.