We consider a topological insulator (TI) of spherical geometry and numerically investigate the influence of disorder on the density of surface states. The energy spectrum of the spherical TI surface is discrete, for a sphere of finite radius, and can be truncated by imposing a high-energy cutoff at the scale of the bulk band gap. To this clean system we add a surface disorder potential of the most general hermitian form, V = V 0 (θ, φ)1 1 + V(θ, φ) · σ, where V 0 describes the spin-independent part of the disorder and the three components of V describe the spin-dependent part. We expand these four disorder functions in spherical harmonics and draw the expansion coefficients randomly from a four-dimensional, zero-mean gaussian distribution. Different strengths and classes of disorder are realized by specifying the 4 × 4 covariance matrix. For each instantiation of the disorder, we solve for the energy spectrum via exact diagonalization. Then we compute the disorder-averaged density of states, ρ(E), by averaging over 200,000 different instantiations. Disorder broadens the Landau-level delta functions of the clean density of states into peaks that decay and merge together. If the spin-dependent term is dominant, these peaks split due to the breaking of the degeneracy between time-reversed partner states. Increasing disorder strength pushes states closer and closer to zero energy (the Dirac point), resulting in a low-energy density of states that becomes nonzero for sufficient disorder, typically approaching an energy-independent saturation value, for most classes of disorder. But for purely spin-dependent disorder with V either entirely out-of-surface or entirely in-surface, we identify intriguing disorder-induced features in the vicinity of the Dirac point. In the out-of-surface case, a new peak emerges at zero energy. In the in-surface case, we see a symmetryprotected zero at zero energy, with ρ(E) increasing linearly toward nonzero-energy peaks. These striking features are explained in terms of the breaking (or not) of two chiral symmetries of the clean Hamiltonian.