For large $k$ the matrix product $\epsilon\epsilon^\top\rightarrow k I_k$, so $y\rightarrow (A+kBB^\top)^{1/2}\epsilon$ and the matrix elements of $y$ are Gaussian with zero mean and covariance matrix $C=A+kBB^\top$.