Voy a seguir el trabajo iniciado en David H la respuesta.
Se ha demostrado que hay que
$$
I_\alpha(n) = \int_{-\frac{\pi}{2}}^{\frac{\pi}{2}} (\alpha + \sin x)^n \cos^2 x\,dx = \int_{\alpha -1}^{\alpha+1} u^n\sqrt{1-(u-\alpha)^2}\,du.
$$
Se consideran los casos $\alpha \geq 1$, $0 < \alpha < 1$, $\alpha = 0$, $-1 < \alpha < 0$, y $\alpha \leq -1$.
Caso 1: $\alpha \geq 1$. Comenzar haciendo la sustitución de $u = (\alpha+1)e^{-t}$ para obtener
$$
I_\alpha(n) = (\alpha+1)^{n+1} \int_0^{-\log\left(\frac{\alpha-1}{\alpha+1}\right)} e^{nt} e^{-t} \sqrt{1-\left[(\alpha+1)e^{-t}-\alpha\right)^2}\,dt.
$$
Desde
$$
\begin{align}
&e^{-t} \sqrt{1-\left[(\alpha+1)e^{-t}-\alpha\right]^2}\,dt \\
&= \sqrt{2+2\alpha}\, t^{1/2} \left[1-\frac{\alpha+6}{4}t - \frac{3\alpha^2-36\alpha-116}{96}t^2 - \frac{\alpha^3-6\alpha^2+36\alpha+88}{128}t^3 + O(t^4)\right]
\end{align}
$$
como $t \to 0^+$, podemos aplicar Watson Lema para obtener la expresión asintótica
$$
\begin{align}
I_\alpha(n) &= \sqrt{\frac{\pi}{2}}\frac{(\alpha+1)^{n+3/2}}{n^{3/2}}\left[1 -\frac{3 (\alpha+6)}{8 n} - \frac{5 \left(3\alpha^2-36\alpha-116\right)}{128 n^2} \right. \\
&\qquad\qquad \left. - \frac{105 \left(\alpha^3-6\alpha^2+36\alpha+88\right)}{1024 n^3} + O\left(\frac{1}{n^{4}}\right)\right].
\end{align}
$$
El uso de más términos en la expansión de la serie de el integrando dará más términos de la expansión asintótica de $I_\alpha(n)$.
Caso 2: $0 < \alpha < 1$. Hemos dividido la integral en dos piezas
$$
I_\alpha(n) = \int_{0}^{\alpha+1} u^n\sqrt{1-(u-\alpha)^2}\,du + \int_{\alpha -1}^{0} u^n\sqrt{1-(u-\alpha)^2}\,du = J_1(n) + J_2(n).
$$
En el primero hacemos la sustitución de $u = (\alpha+1)e^{-t}$ como antes, y en el segundo hacemos la sustitución de $u=(\alpha-1)e^{-t}$. Las integrales se convierten
$$
J_1(n) = (\alpha+1)^{n+1} \int_0^{\infty} e^{nt} e^{-t} \sqrt{1-\left[(\alpha+1)e^{-t}-\alpha\right)^2}\,dt,
$$
$$
J_2(n) = (-1)^n (1-\alpha)^{n+1} \int_0^{\infty} e^{nt} e^{-t} \sqrt{1-\left[(1-\alpha)e^{-t}+\alpha\right)^2}\,dt.
$$
Mediante la aplicación de Watson lema a cada uno de estos de forma individual nos encontramos con que
$$
J_1(n) \sim \sqrt{\frac{\pi}{2}}\frac{(\alpha+1)^{n+3/2}}{n^{3/2}}
$$
y
$$
J_2(n) \sim (-1)^n \sqrt{\frac{\pi}{2}}\frac{(1-\alpha)^{n+3/2}}{n^{3/2}},
$$
por tanto, los términos procedentes de la primera integral dominar a aquellos que vienen de la segunda. Por lo tanto $I_\alpha(n)$ tiene el mismo asintótica de expansión como en el Caso 1.
Caso 3: $\alpha = 0$. En este caso,$J_2(n) = (-1)^n J_1(n)$, por lo que
$$
I_0(n) = \Bigl(1+(-1)^n\Bigr)\int_0^{\infty} e^{nt} e^{-t} \sqrt{1-e^{-2t}}\,dt.
$$
Aquí
$$
e^{-t} \sqrt{1-e^{-2t}} = \sqrt{2}\,t^{1/2} \left[1-\frac{3}{2}t+\frac{29}{24}t^2-\frac{11}{16}t^3 + O(t^4)\right],
$$
por lo que la aplicación de Watson lema de los rendimientos
$$
I_0(n) = \sqrt{\frac{\pi}{2}} \frac{1+(-1)^n}{n^{3/2}} \left[1-\frac{9}{4 n}+\frac{145}{32 n^2}-\frac{1155}{128 n^3} + O\left(\frac{1}{n^4}\right)\right].
$$
Caso 4: $-1 < \alpha < 0$. Este caso es similar a la del Caso 2. La única diferencia es que los términos de $J_2(n)$ dominar. Tenemos
$$
\begin{align}
I_\alpha(n) &= (-1)^n \sqrt{\frac{\pi}{2}} \frac{(1-\alpha)^{n+3/2}}{n^{3/2}} \left[1+\frac{3(\alpha-6)}{8 n} - \frac{5\left(3\alpha^2+36\alpha-116\right)}{128 n^2} \right. \\
&\qquad\qquad \left. + \frac{105 \left(\alpha^3+6\alpha^2+36\alpha-88\right)}{1024 n^3} + O\left(\frac{1}{n^4}\right)\right].
\end{align}
$$
Caso 5: $\alpha < -1$. Este caso es similar a la del Caso 1 (y su relación con el Caso 2). Si hacemos la sustitución de $u=(\alpha-1)e^{-t}$ $I_\alpha(n)$ y aplicada Watson lema nos encontraríamos con que $I_\alpha(n)$ tiene el mismo asintótica de expansión como en el Caso 4.
En resumen,
Si $a > 0$
$$
\begin{align}
I_\alpha(n) &= \sqrt{\frac{\pi}{2}}\frac{(\alpha+1)^{n+3/2}}{n^{3/2}}\left[1 -\frac{3 (\alpha+6)}{8 n} - \frac{5 \left(3\alpha^2-36\alpha-116\right)}{128 n^2} \right. \\
&\qquad\qquad \left. - \frac{105 \left(\alpha^3-6\alpha^2+36\alpha+88\right)}{1024 n^3} + O\left(\frac{1}{n^{4}}\right)\right]
\end{align}
$$
como $n \to \infty$.
Si $a = 0$
$$
I_0(n) = \sqrt{\frac{\pi}{2}} \frac{1+(-1)^n}{n^{3/2}} \left[1-\frac{9}{4 n}+\frac{145}{32 n^2}-\frac{1155}{128 n^3} + O\left(\frac{1}{n^4}\right)\right]
$$
como $n \to \infty$.
Si $a < 0$
$$
\begin{align}
I_\alpha(n) &= (-1)^n \sqrt{\frac{\pi}{2}} \frac{(1-\alpha)^{n+3/2}}{n^{3/2}} \left[1+\frac{3(\alpha-6)}{8 n} - \frac{5\left(3\alpha^2+36\alpha-116\right)}{128 n^2} \right. \\
&\qquad\qquad \left. + \frac{105 \left(\alpha^3+6\alpha^2+36\alpha-88\right)}{1024 n^3} + O\left(\frac{1}{n^4}\right)\right]
\end{align}
$$
como $n \to \infty$.
A continuación están algunas de las parcelas de los tres regímenes de $5 \leq n \leq 40$. Los puntos azules son los valores de $I_\alpha(n)$ obtenido numéricamente y los puntos rojos son los asintótica aproximaciones obtenidas anteriormente. Los dos primeros términos de la asymptotics se utilizan.
Figura 1: $\alpha = 1/10$.
Figura 2: $\alpha = 0$.
Figura 3: $\alpha = -1/10$.