This document presents a methodology for modeling mortality rates across Canadian provinces. It reviews literature on single and multi-population mortality models. Mortality indices from 1921-2009 for 9 provinces are analyzed using unit root tests and the Johansen cointegration test, which finds common trends among the indices. Vector error correction (VECM) and vector autoregressive (VAR) models are estimated and show better fit for female mortality. The models are used to project and forecast mortality indices 50 years into the future and to price annuities for cohorts from 1960-2000, with VECM producing more accurate pricing. In conclusion, mortality trends are declining across provinces with common factors, and VECM performs better than VAR or ARIMA models