Mining induced seismic activity and rockbursting are critical concerns for many underground operations. Seismic activity may arise from the crushing of highly stressed volumes of rock around mine openings or from shear motion on planes of weakness. Shear-slip on major planes of weakness such as faults, shear zones and weak contacts has long been recognized as a dominant mode of failure in underground mines. In certain circumstances, it can generate large seismic events and induce substantial damage to mine openings. The Big Bell Gold mine began experiencing major seismic activity and resultant damage in 1999. Several seismic events were recorded around the second graphitic shear between April 2000 and February 2002. It is likely that the seismic activity occurred as a result of the low strength of the shear structure combined with the high level of mining induced stresses. The stability of the second graphitic shear was examined in order to gain a better understanding of the causes and mechanisms of the seismic activity recorded in the vicinity of the shear structure as mining advanced. The data were derived from the observation of the structure exposures, numerical modelling and seismic monitoring. The numerical modelling predictions and the interpreted seismic monitoring data were subsequently compared in order to identify potential relationships between the two. This thesis proposes the Incremental Work Density (IWD) as a measure to evaluate the relative likelihood of shear-slip induced seismic activity upon major planes of weakness. IWD is readily evaluated using numerical modelling and is calculated as the product of the average driving shear stress and change in inelastic shear deformation during a given mining increment or step. IWD is expected to correlate with shear-slip induced seismic activity in both space and time. In this thesis, IWD was applied to the case study of the second graphitic shear at the Big Bell mine. Exposures of the second graphitic shear yielded information about the physical characteristics of the structure and location within the mine. Numerical modelling was used to examine the influence of mining induced stresses on the overall behaviour of the shear structure. A multi-step model of the mine was created using the three- dimensional boundary element code of Map3D. The shear structure was physically incorporated into the model in order to simulate inelastic shear deformation. An elasto-plastic Mohr-Coulomb material model was used to describe the structure behaviour. The structure plane was divided into several elements in order to allow for the comparison of the numerical modelling predictions and the interpreted seismic data. Stress components, deformation components and IWD values were calculated for each element of the shear structure and each mining step. The seismic activity recorded in the vicinity of the second graphitic shear was back analysed. The seismic data were also gridded and smoothed. Gridding and smoothing of individual seismic moment and seismic energy values resulted in the definition of indicators of seismic activity for each element and mining step. The numerical model predicted inelastic shear deformation upon the second graphitic shear as mining advanced. The distribution of modelled IWD suggested that shear deformation was most likely seismic upon a zone below the stopes and most likely aseismic upon the upper zone of the shear structure. The distribution of seismic activity recorded in the vicinity of the shear structure verified the above predictions. The seismic events predominantly clustered upon the zone below the stopes. The results indicated that the seismic activity recorded in the vicinity of the second graphitic shear was most likely related to both the change in inelastic shear deformation and the level of driving shear stress during mechanical shearing. Time distribution of the seismic events also indicated that shear deformation and accompanying seismic activity were strongly influenced by mining and were time-dependant. Seismic activity in the vicinity of the second graphitic shear occurred as a result of the overall inelastic shear deformation of the shear structure under mining induced stresses. A satisfactory relationship was found between the spatial distribution of modelled IWD upon the shear structure and the spatial distribution of interpreted seismic activity (measured as either smoothed seismic moment or smoothed seismic energy). Seismic activity predominantly clustered around a zone of higher IWD upon the second graphitic shear as mining advanced. However, no significant statistical relationship was found between the modelled IWD and the interpreted seismic activity. The lack of statistical relationship between the modelled and seismic data may be attributed to several factors including the limitations of the techniques employed (e.g. Map3D modelling, seismic monitoring) and the complexity of the process involved.