It is shown how suitably scaled, order-m moments, Dm±, of the Elsässer vorticity fields in three-dimensional magnetohydrodynamics (MHD) can be used to identify three possible regimes for solutions of the MHD equations with magnetic Prandtl number PM=1. These vorticity fields are defined by ω±=curlz±=ω±j, where z± are Elsässer variables, and where ω and j are, respectively, the fluid vorticity and current density. This study follows recent developments in the study of three-dimensional Navier-Stokes fluid turbulence [Gibbon, Nonlinearity 27, 2605 (2014)NONLE50951-771510.1088/0951-7715/27/10/2605]. Our mathematical results are then compared with those from a variety of direct numerical simulations, which demonstrate that all solutions that have been investigated remain in only one of these regimes which has depleted nonlinearity. The exponents q± that characterize the inertial range power-law dependencies of the z± energy spectra, E±(k), are then examined, and bounds are obtained. Comments are also made on (a) the generalization of our results to the case PM≠1 and (b) the relation between Dm± and the order-m moments of gradients of magnetohydrodynamic fields, which are used to characterize intermittency in turbulent flows. © 2016 American Physical Society.