Authors: Sainsbury, BL; Sainsbury, DP; Osorio, A; Carroll, D; Lett, J

Open access courtesy of:


Cite As:
Sainsbury, BL, Sainsbury, DP, Osorio, A, Carroll, D & Lett, J 2022, 'A coupled modelling approach for discontinuous subsidence at the Cadia East mine', in Y Potvin (ed.), Caving 2022: Fifth International Conference on Block and Sublevel Caving, Australian Centre for Geomechanics, Perth, pp. 797-804,

Download citation as:   ris   bibtex   endnote   text   Zotero

It is known that in addition to the production schedule, a number of orebody, geologic and topographic features can influence the nature and extent of surface subsidence related to a cave mine. Some of these factors include the shape and depth of the orebody/extraction level in plan, strength of the rock mass including the intact blocks, joints and geological features such as faults and dykes, ground surface profile, including historical surface mining and the presence of nearby underground excavations (Sainsbury 2012). Brady and Brown (2006) classify subsidence into two types, continuous and discontinuous. Continuous (or trough) subsidence refers to the formation of a smooth (gradual) change in profile that does not have step changes. Discontinuous subsidence involves large surface displacements (both tension and shear) and the formation of steps/breaks in the surface profile. Discontinuous subsidence is strongly affected by the presence of pre-existing features within the rock mass. As such, it cannot be accurately assessed without the inclusion of such features within a three-dimensional numerical model. A coupled continuous–discontinuous approach has been developed to predict caving and subsidence efficiently and accurately at the Cadia East mine over an eight-year production life (or 150 MT drawn). Numerically simulated results with the coupled approach are compared to in situ observations. It can be seen that in the case of sub-vertical structures within the near-surface rock mass, traditional strain criteria applied alongside continuum models are unable to capture the true shape or extent of cracking and that discontinuum models provide greater accuracy.

Keywords: caving, subsidence, continuum-discontinuum

Abel, JF & Lee, FT 1980, Subsidence Potential in Shale and Crystalline Rocks, open-file report, U.S. Geological Survey.
Brady, BHG & Brown, ET 2006, Rock Mechanics for Underground Mining, 3rd edn, Springer, Dordrecht.
Cavieres, P, Gaete, S, Lorig, L & Gomez, P 2003, ‘Three-dimensional analysis of fracturing limits induced by large scale underground mining at El Teniente mine’, in P Culligan, H Einstein & A Whittle (eds), Proceedings of the 39th U.S. Rock Mechanics Symposium, pp. 893–900
Clark, I 2006, ‘Simulation of rock mass strength using ubiquitous joints’, in R Hart, P Varona (eds), Numerical Modeling in Geomechanics 2006: Proceedings of the 4th International FLAC Symposium, Itasca Consulting Group, Minneapolis, paper no. 08-07.
Davy, P, Darcel, C, Le Goc, R, Munier, R, Selroos, JO & Mas Ivars, D 2018, ‘DFN, why, how and what for, concepts, theories and issues’, ARMA Virtual Conference: 2nd International Discrete Fracture Network Engineering Conference, American Rock Mechanics Association, Alexandria.
Itasca Consulting Group 2022a, 3DEC – Distinct-Element Modeling of Jointed and Blocky Material in 3D, computer software, Itasca Consulting Group, Minneapolis,
Itasca Consulting Group 2022b, FLAC3D – Explict Continuum Modeling of Non-Linear Material Behavior in 3D, version 7.0, computer software, Minneapolis,
Lambe, TW 1973, ‘Predictions in soil engineering’, Géotechnique, vol. 23, pp. 151–202,
Sainsbury, B 2012, A Model for Cave Propagation and Subsidence Assessment in Jointed Rock Masses, PhD thesis, UNSW Sydney, Kensington.
Sainsbury, B 2018, ‘Consideration of the volumetric changes that accompany rock mass failure’, Rock Mechanics and Rock Engineering, vol. 52,
Sainsbury, B, Sainsbury, D & Carroll, D 2018, ‘Back-analysis of PC1 cave propagation and subsidence behaviour at the Cadia East mine’, in Y Potvin & J Jakubec (eds), Caving 2018: Proceedings of the Fourth International Symposium on Block and Sublevel Caving, Australian Centre for Geomechanics, Perth, pp. 167–178,
Sainsbury, B & Stokel, B-M 2012, ‘Historical assessment of caving induced subsidence at the Kiirunavaara Lake orebody’, MassMin 2012: Proceedings of the 6th International Conference & Exhibition on Mass Mining, Canadian Institute of Mining, Metallurgy and Petroleum, Westmount.
Sainsbury, B & Sainsbury, D 2017, ‘Practical use of the ubiquitous-joint constitutive model for the simulation of anisotropic rock masses’, Rock Mechanics and Rock Engineering, vol. 50, pp. 1507–1528,
Sainsbury, D, Sainsbury, B, Paetzold, H-D, Lourens, P & Vakili, A 2016a, ‘Caving-induced subsidence behaviour of Lift 1 at the Palabora block cave mine’, MassMin 2016: Proceedings of the Seventh International Conference & Exhibition On Mass Mining, Australasian Institute of Mining and Metallurgy, Melbourne.
Sainsbury, D, Sainsbury, B & Lorig, L 2010, ‘Investigation of caving induced subsidence at the abandoned Grace Mine’, Transactions of the Institutions of Mining and Metallurgy, Section A: Mining Technology, vol. 119, pp. 151–161,
Sainsbury, D, Sainsbury, B & Sweeney, E 2016b, ‘Three-dimensional analysis of complex anisotropic slope instability at MMG’s Century Mine’, Transactions of the Institutions of Mining and Metallurgy, Section A: Mining Technology, vol. 125,
Vyazmensky, A, Elmo, D, Stead, D & Rance JR 2008, ‘Numerical analysis of the influence of geological structures on development of surface subsidence associated with block caving mining’, MassMin 2008: Proceedings of the 5th International Conference and Exhibition on Mass Mining, Luleå, pp. 467–475.

© Copyright 2022, Australian Centre for Geomechanics (ACG), The University of Western Australia. All rights reserved.
Please direct any queries or error reports to