In this thesis I study a variety of two-dimensional turbulent systems using a mixed analytical, phenomenological and numerical approach. The systems under consideration are governed by the two-dimensional Navier-Stokes (2DNS), surface quasigeostrophic (SQG), alpha-turbulence and magnetohydrodynamic (MHD) equations. The main analytical focus is on the number of degrees of freedom of a given system, defined as the least value $N$ such that all $n$-dimensional ($n$ ≥ $N$) volume elements along a given trajectory contract during the course of evolution. By equating $N$ with the number of active Fourier-space modes, that is the number of modes in the inertial range, and assuming power-law spectra in the inertial range, the scaling of $N$ with the Reynolds number $Re$ allows bounds to be put on the exponent of the spectrum. This allows the recovery of analytic results that have until now only been derived phenomenologically, such as the $k$[superscript(-5/3)] energy spectrum in the energy inertial range in SQG turbulence. Phenomenologically I study the modal interactions that control the transfer of various conserved quantities. Among other results I show that in MHD dynamo triads (those converting kinetic into magnetic energy) are associated with a direct magnetic energy flux while anti-dynamo triads (those converting magnetic into kinetic energy) are associated with an inverse magnetic energy flux. As both dynamo and anti-dynamo interacting triads are integral parts of the direct energy transfer, the anti-dynamo inverse flux partially neutralises the dynamo direct flux, arguably resulting in relatively weak direct energy transfer and giving rise to dynamo saturation. These theoretical results are backed up by high resolution numerical simulations, out of which have emerged some new results such as the suggestion that for alpha turbulence the generalised enstrophy spectra are not closely approximated by those that have been derived phenomenologically, and new theories may be needed in order to explain them.