The spatial structure of the magnetohydrodynamic shear Alfvén and cusp continuum modes in axisymmetric toroidal geometry is analyzed for an incompressible plasma. The normal component $\xi_{\psi}$ of the plasma displacement is found to have an oscillatory type of singularity, $\xi_{\psi} \sim (\psi - \psi_0)^{\delta} \sim \sin(|\delta| \ln|\psi-\psi_0| + \mbox{const})$, or, for example, in the case of up/down symmetry, a purely logarithmic behavior. This result was also obtained in the more general case of a compressible plasma (Salat, A. and Tataronis, J. A. 1999 Phys. Plasmas6, 3207). However, in that work the incompressible limit could not be taken because the continuum equations were written in terms of a variable that becomes ill-defined in the limit. The present ab initio derivation shows that, in general, the incompressible limit in the previous work, $\gamma \rightarrow\infty$, where $\gamma$ is the ratio of the specific heats, does give the correct spatial dependence. Furthermore, we prove that it is possible to recover the oscillatory type of singularity for the cylindrical screw pinch from the logarithmic singularity of an up/down symmetric torus.