Preventing burst fracture is critical in maintaining quality of life in patients affected by metastatic cancer to the spine. The objective of this study was to use parametric finite element modeling to improve burst fracture risk assessment and treatment options for patients with spinal metastases. The model was used to evaluate the effects of complex loading and geometric factors in the metastatic spine. Axial loading was found to be the primary load type contributing to burst fracture initiation. Upper thoracic vertebrae were found to be at increased risk of burst fracture. Fracture risk was reduced by presence of the ribcage and increased kyphotic angle. Transcortical tumours were found to have a lower risk of burst fracture as compared to contained tumour scenarios. Protocols were developed for pathologic burst fracture risk assessment and prophylactic treatment based on these findings. This work is important in directing treatment for patients with spinal metastases.