I would like to show that $$ g_{\sigma}(x, y) * g_{\tau}(x, y)=g_{\sqrt{\sigma^{2}+\tau^{2}}}\left(x, y\right) $$ where $g_{\sigma}(x, y)=\frac{1}{2 \pi \sigma^{2}} e^{-\frac{x^{2}+y^{2}}{2 \sigma^{2}}}$. My attempt so far:
Since the Fourier transform of $g_{\sigma}$ is given by $$ \hat{g}_{\sigma}(x, y) = e^{-(\alpha^2 + \beta^2)2\pi\sigma^2} $$ However I simply cannot simplify/get the product $$ \hat{g}_{\sigma}(x, y) \hat{g}_{\tau}(x, y)= e^{-(\alpha^2 + \beta^2)2\pi\sigma^2} e^{-(\alpha^2 + \beta^2)2\pi\tau^2} $$ to $g_{\sqrt{\sigma^{2}+\tau^{2}}}\left(x, y\right)$ - am I on the wrong track here, or is it just some manipulation that needs to be done?
You are almost there. Use $$e^Ae^B=e^{A+B}$$ In your case $$e^{-(u^2+v^2)2\pi\sigma^2}e^{-(u^2+v^2)2\pi\tau^2}=e^{-(u^2+v^2)2\pi(\sigma^2+\tau^2)}=e^{-(u^2+v^2)2\pi(\sqrt{\sigma^2+\tau^2})^2}=\mathcal{F}\{G(x, y, \sqrt{\sigma^2+\tau^2})\}$$