Como una respuesta parcial que puede mostrar la primera ecuación es válida:
En la parte superior de la superficie de la esfera,
$$\vec r=\langle x,y,z\rangle=\langle x,y,\sqrt{r^2-x^2-y^2}\rangle$$
Así
$$d\vec r=\langle 1,0,\frac{-x}{\sqrt{r^2-x^2-y^2}}\rangle dx+\langle 0,1,\frac{-y}{\sqrt{r^2-x^2-y^2}}\rangle dx$$
Y, a continuación,
$$\begin{align}d\vec A&=\pm\langle 1,0,\frac{-x}{\sqrt{r^2-x^2-y^2}}\rangle dx\times\langle 0,1,\frac{-y}{\sqrt{r^2-x^2-y^2}}\rangle dx\\
&=\pm\langle\frac{x}{\sqrt{r^2-x^2-y^2}},\frac{y}{\sqrt{r^2-x^2-y^2}},1\rangle dx\,dy\\
&=\langle\frac{x}{\sqrt{r^2-x^2-y^2}},\frac{y}{\sqrt{r^2-x^2-y^2}},1\rangle dx\,dy\end{align}$$
Porque en la apariencia normal de los puntos hasta aquí. Entonces
$$d^2A=\left|\left|d^2\vec A\right|\right|=\frac{r\,dx\,dy}{\sqrt{r^2-x^2-y^2}}$$
En la parte inferior de la superficie de la esfera,
$$\vec r=\langle x,y,z\rangle=\langle x,y,-\sqrt{r^2-x^2-y^2}\rangle$$
Y podemos trabajar fuera
$$d^2A=\left|\left|d^2\vec A\right|\right|=\frac{r\,dx\,dy}{\sqrt{r^2-x^2-y^2}}$$
Tal y como estaba en la superficie superior.
Comience con el lado derecho:
$$\begin{align}\int\int_{S_r}\cos(x+y+z)d^2A&=\int_{-r}^r\int_{-\sqrt{r^2-x^2}}^{\sqrt{r^2-x^2}}\left[\cos\left(x+y-\sqrt{r^2-x^2-y^2}\right)\right.\\
&+\left.\cos\left(x+y+\sqrt{r^2-x^2-y^2}\right)\right]\frac{r\,dy}{\sqrt{r^2-x^2-y^2}}dx\end{align}$$
El exterior de la $(x)$ integral tiene un dominio que es un círculo paralelo a la $yz$-plane y el interior de la integral sólo necesita los dos puntos que tienen el mismo $x$$y$. Ahora podemos ampliar la cosenos utilizando el ángulo de la fórmula de la suma para obtener
$$\begin{align}\int\int_{S_r}\cos(x+y+z)d^2A&=\int_{-r}^r\int_{-\sqrt{r^2-x^2}}^{\sqrt{r^2-x^2}}2\cos(x+y)\cos\left(\sqrt{r^2-x^2-y^2}\right)\frac{r\,dy}{\sqrt{r^2-x^2-y^2}}dx\\
&=\int_{-r}^r\int_{-\sqrt{r^2-x^2}}^{\sqrt{r^2-x^2}}\cos(x+y)\cos\left(-\sqrt{r^2-x^2-y^2}\right)\frac{r\,dy}{\sqrt{r^2-x^2-y^2}}dx\\
&+\int_{-r}^r\int_{-\sqrt{r^2-x^2}}^{\sqrt{r^2-x^2}}\cos(x+y)\cos\left(\sqrt{r^2-x^2-y^2}\right)\frac{r\,dy}{\sqrt{r^2-x^2-y^2}}dx\\
&=\int\int_{S_r}\cos(x+y)\cos z\,d^2A\end{align}$$
Vamos a configurar la misma integral con $x$ exterior, $z$ media y $y$ interior:
$$\begin{align}\int\int_{S_r}\cos(x+y)\cos z\,d^2A&=\int_{-r}^r\int_{-\sqrt{r^2-x^2}}^{\sqrt{r^2-x^2}}\left[\cos\left(x-\sqrt{r^2-x^2-z^2}\right)\cos z\right.\\
&+\left.\cos\left(x+\sqrt{r^2-x^2-z^2}\right)\cos z\right]\frac{r\,dz}{\sqrt{r^2-x^2-z^2}}dx\\
&=\int_{-r}^r\int_{-\sqrt{r^2-x^2}}^{\sqrt{r^2-x^2}}2\cos x\cos\left(\sqrt{r^2-x^2-z^2}\right)\cos z\frac{r\,dz}{\sqrt{r^2-x^2-z^2}}dx\\
&=\int_{-r}^r\int_{-\sqrt{r^2-x^2}}^{\sqrt{r^2-x^2}}\cos x\cos\left(-\sqrt{r^2-x^2-z^2}\right)\cos z\frac{r\,dz}{\sqrt{r^2-x^2-z^2}}dx\\
&+\int_{-r}^r\int_{-\sqrt{r^2-x^2}}^{\sqrt{r^2-x^2}}\cos x\cos\left(\sqrt{r^2-x^2-z^2}\right)\cos z\frac{r\,dz}{\sqrt{r^2-x^2-z^2}}dx\\
&=\int\int_{S_r}\cos x\cos y\cos z\,d^2A\end{align}$$
EDIT: Ahora tenemos que terminar la integral. Comenzamos con la mano derecha y transformar a
$$\begin{align}u&=\frac1{\sqrt3}(x+y+z)\\
v&=\frac1{\sqrt6}(-2x+y+z)\\
w&=\frac1{\sqrt2}(z-y)\end{align}$$
Este es sólo un rígido rotación sobre el origen, por lo que
$$\begin{align}\int\int_{S_r}\cos(x+y+x)d^2A&=\int\int_{S_r}\cos\left(\sqrt3u\right)d^2A\\
&=\int_{-r}^r\int_{-\sqrt{r^2-u^2}}^{\sqrt{r^2-u^2}}2\cos\left(\sqrt3u\right)\frac{r\,dv}{\sqrt{r^2-u^2-v^2}}du\\
&=\int_{-r}^r\int_{-\pi/2}^{\pi/2}2\cos\left(\sqrt3u\right)\frac{r\sqrt{r^2-u^2}\cos\theta\,d\theta}{\sqrt{r^2-u^2}\cos\theta}du\\
&=2\pi r\int_{-r}^r\cos\left(\sqrt3u\right)du\\
&=\frac{2\pi r}{\sqrt3}\left.\sin\left(\sqrt3u\right)\right|_{-r}^r\\
&=\frac{4\pi r}{\sqrt3}\sin\left(\sqrt3r\right)\end{align}$$
EDIT: Lástima que el factor de $\frac r{\sqrt{r^2-u^2-v^2}}$ en el área de elemento destruido mi hermosa solución con las funciones de Bessel y nos dejó con esta prosaica respuesta.