If $X$ and $Y$ are independent random variables that are normally distributed, then their sum is also normally distributed. i.e., if
$X\sim N(\mu_X, \sigma_X^2)$
$Y \sim N(\mu_Y, \sigma_Y^2)$
$Z=X+Y$
and $X$ and $Y$ are independent, then
$Z \sim N(\mu_X + \mu_Y, \sigma_X^2 + \sigma_Y^2).$
You can prove it using both characteric functions and convolutions. I will refer you to read this page for the proofs.
In the event that the variables $X$ and $Y$ are jointly normally distributed correlated random variables, then $Z=X + Y$ is still normally distributed
and the mean is the sum of the means. However, the variances are not additive
due to the correlation. Indeed,
$\sigma_{X+Y} = \sqrt{\sigma_X^2+\sigma_Y^2+2\rho\sigma_X \sigma_Y}$
This one is also shown in the previous link.
You can also find this using the fact that $Q^{'}=[X \; Y]$ is bivariate normal with variables $X$ and $Y$ and distributed as $N_{2}(\mu,\Sigma).$ Then $Z=AQ$ with $A=[1\;1]$ is distributed as $N_{q}(A\mu,A \Sigma A^{'})$, where $q$ is the number of linear combinations (or the number of rows in $A$). Each row in $A$ corresponds to a linear combination of the variables $X$ and $Y$. Here we have only 1 row. So, $q=1$. Put $\sigma_X^2$, $\sigma_Y^2$ and $\sigma_{XY}$ in $\Sigma_{(2\times2)}$ and find the result for the correlated case.