Ballooning instabilities are investigated in three-dimensional magnetic toroidal plasma confinement systems with low global magnetic shear. The lack of any continuous symmetry in the plasma equilibrium can lead to these modes being localized along the field lines by a process similar to Anderson localization. This produces a multibranched local eigenvalue dependence, where each branch corresponds to a different unit cell of the extended covering space in which the eigenfunction peak resides. These phenomena are illustrated numerically for the three-field-period heliac H-1, and contrasted with an axisymmetric s-α tokamak model. The localization allows a perturbative expansion about zero shear, enabling the effects of shear to be investigated. Ballooning instabilities are pressure-driven ideal magnetohydrodynamic (MHD) instabilities which limit the maximum β (plasma pressure/magnetic pressure) that can be obtained in a plasma. They are localized about regions where the field lines are concave to the plasma, which are known as unfavourable regions of curvature. Another localizing influence is the magnetic shear, which measures the rate at which neighboring field lines at different minor radii separate as they wind their way around the torus. Large shear helps stabilize these modes, thereby playing an important role in the MHD stability. In this paper however we consider the effects of very small or zero shear, such as occurs in the heliac class of stellarators or in the shear-reversal layers of an advanced tokamak.We begin by making the usual assumption that the magnetic field lines map out nested flux surfaces, or magnetic surfaces. These are labeled using a normalizedtoroidal-flux variable s, which varies between zero at the center of the plasma and unity at the plasma edge. Within each surface the straight-field-line poloidal θ and toroidal ζ angle variables are defined such that the field lines appear as straight lines in the (θ, ζ) plane. The magnetic field may then be written B = ∇ζ×∇ψ − q∇θ×∇ψ ≡ ∇α×∇ψ, where the field-line label α ≡ ζ − qθ. Here, 2πψ represents the poloidal magnetic flux, while q = q(s) is the safety factor (inverse of rotational transform), which is equal to the average number of toroidal circuits traversed by a field line per poloidal circuit traversed around the torus.Ballooning modes can be characterized as having a long parallel and short perpendicular wavelength with respect to the field lines. By ordering the perpendicular wavelength to be small and expanding to lowest order in an asymptotic series the local mode behavior can be expressed by a one-dimensional equation along a field line [1]. Taking the plasma to be incompressible, the ballooning equation may be written [2]where the eigenfunction ξ is related to the mode displacement while the eigenvalue λ is equal to the mode growth rate squared. This represents the local stability, local to a field line. In forming global modes, ray tracing must be performed in the three-dimensional λ phase space to determine which of these local sol...