
We present a lattice calculation of the nucleon isovector axial and induced pseudoscalar form factors on the CLS ensembles using $N_{\rm f}=2$ dynamical flavours of nonperturbatively $\mathcal{O}(a)$improved Wilson fermions and an $\mathcal{O}(a)$improved axial current together with the pseudoscalar density. Excitedstate effects in the extraction of the form factors are treated using a variety of methods, with a detailed discussion of their respective merits. The chiral and continuum extrapolation of the results is performed both using formulae inspired by Heavy Baryon Chiral Perturbation Theory (HBChPT) and a global approach to the form factors based on a chiral effective theory (EFT) including axial vector mesons. Our results indicate that careful treatment of excitedstate effects is important in order to obtain reliable results for the axial form factors of the nucleon, and that the main remaining error stems from the systematic uncertainties of the chiral extrapolation. As final results, we quote $g_{\rm A} = 1.278 \pm 0.068\genfrac{}{}{0pt}{1}{+0.000}{0.087}$, $\langle r_{\rm A}^2\rangle = 0.360 \pm 0.036\genfrac{}{}{0pt}{1}{+0.080}{0.088}~\mathrm{fm}^2$, and $g_{\rm P} = 7.7 \pm 1.8 \genfrac{}{}{0pt}{1}{+0.8}{2.0}$ for the axial charge, axial charge radius and induced pseudoscalar charge, respectively, where the first error is statistical and the second is systematic.