The MHD model enables derivation and analysis of the rich structure of geodesic acoustic modes (GAMs) and zonal modes in axisymmetric magnetic confined plasmas. The modes are identifiable from a single dispersion relation as two branches of slow magnetosonic continua. The lower frequency branch can be identified as a zonal flow, which in the simplified limit of static plasmas, has vanishing magnetic component. It is shown in this contribution that axisymmetric, and lesser known nonaxisymmetric, zonal modes can be derived from MHD and kinetic models. The work provides a comprehensive derivation of the GAMs and zonal flow continua in static and stationary toroidally rotating plasmas, and investigates the exact solution and structure of a generalised family of zonal modes.