There are plenty of ways they can be random and correlated. In fact, that’s the whole problem: We don’t know how they’re correlated, so we can’t answer the question.
There’s a closed form for the pdf in terms of the Bessel function BesselI(0,r) [that’s a capital-i, not a small-L], for what it’s worth: Diagonalize the covariance and rescale so that you have uncorrelated Gaussian random variables U=N(0,1) and V=N(0,s[sup]2[/sup]) with s<=1, and Z=sqrt(U[sup]2[/sup]+V[sup]2[/sup]). (The smaller of the two principal-axis variances s[sup]2[/sup] will be our extra parameter.)
Converting to polar coordinates and doing the angular integral gives the pdf,
f(z) = (z/s) Exp[-(z[sup]2[/sup]/4)(1/s[sup]2[/sup]+1)] BesselI(0,(z[sup]2[/sup]/4)(1/s[sup]2[/sup]-1)) .
Note that when s=1 you have the isotropic case and the distribution becomes Rayleigh; when s<<1 the distribution converges to the right half of the normal distribution.