This paper considers the phenomenon of distinct regional frequencies recently observed in some power systems. We first study the causes through dynamic simulations and then develop the mathematical model describing this behaviour. Then, techniques to solve the model are discussed, demonstrating that the post-fault frequency evolution in any given region is equal to the frequency evolution of the Centre Of Inertia plus certain inter-area oscillations. This finding leads to the deduction of conditions for guaranteeing frequency stability in all regions of a power system, a deduction performed using a mixed analyticalnumerical approach that combines mathematical analysis with regression methods on simulation samples. The proposed stability conditions are linear inequalities that can be implemented in any optimisation routine allowing the co-optimisation of all existing ancillary services for frequency support: inertia, multi-speed frequency response, load damping and an optimised largest power infeed. This is the first reported mathematical framework with explicit conditions to maintain frequency stability in a power system exhibiting inter-area oscillations in frequency.