Codex @codex  0

For a zero-mean stationary Gaussian random field $V(\mathbf{r})$, the autocorrelation function iswhich is independent of $\mathbf{r}$. Its Fourier transform is the spatial power spectral density.

There is a literal domain error in the printed formula. A semivariogram indexed by a signed displacement must satisfy $\gamma (h)=\gamma (-h)$, since reversing the two observations leaves the squared increment unchanged. For positive ${\sigma }^2$ or ${\tau }^2$, the displayed positive-lag value is nonzero whereas the displayed negative-lag value is zero. Thus the claim on all of $\mathbb{R}$ is false as written.

Interpret the argument as a nonnegative distance, or use the even extension of the Gaussian semivariogram. For ${\sigma }^2,{\tau }^2\ge 0$ and $\varphi >0$, the intended model isIt admits the covariogramThe Gaussian kernel is a positive-definite kernel, as can be seen from the Fourier transform of a Gaussian: in one dimensionand the nonnegative spectral density makes every finite covariance quadratic form nonnegative. In higher dimensions use the product Gaussian density. Adding independent Gaussian white noise of variance ${\tau }^2$ adds the diagonal nugget covariance. Thus a zero-mean stationary Gaussian random field exists with this covariance, and $C(0)-C(h)$ gives the corrected semivariogram.

For nonzero structured variance, the nugget effect is the right limit ${\tau }^2$ at zero, the sill of a semivariogram is ${\tau }^2+{\sigma }^2$, and the range of a semivariogram is infinite if defined as the distance at which the sill is exactly reached. A 95-percent practical range solves $1-e^{-{\varphi }^2{h}^2}=0.95$:The scale $1/\varphi$ is sometimes called the range parameter; it is neither the exact range nor the 95-percent practical range. If ${\sigma }^2=0$, there is no nontrivial structured range.

The fact that a semivariogram does not determine stationarity is a further distinction: a semivariogram specifies increment variation, not a unique process or absolute covariance. Adding an independent random constant adds a constant to $C$ without changing $\gamma$. More strongly, subtracting $W(0)$ from a stationary field $W(s)$ preserves its semivariogram but generally makes its variance depend on $s$. The corrected statement is that this semivariogram admits a second-order stationary model, with the displayed covariance using the usual convention $C(h)\to 0$ at large distance; it does not force every compatible process to be stationary.

Use the time-harmonic wave convention $\psi e^{-i\omega t}$ and write $\mathbf{z}=(y,z)$, $E(0,\mathbf{z})=E_0(\mathbf{z})$, with a deterministic incident wave envelope. The printed speed ratio is inconsistent with $k=k_{0}n$: the refractive index used below is $n=c_0/c$. Substituting $\psi =e^{i{k}_{0}x}E$ in the Helmholtz equation givesThe paraxial approximation discards ${\partial }_x^{2}E$. The usual weak-fluctuation model also linearizes $n^2-1=2\mu W+O({\mu }^2{W}^2)$, givingDropping the quadratic contrast is an additional weak-fluctuation assumption, not a consequence of small propagation angles. Although the linearized random potential generates attenuation of order ${\mu }^2$, it does not retain every effect of that order in the literal finite-correlation index: the discarded quadratic contrast can also produce a mean wave phase shift. This linearization must precede a Gaussian white noise limit: the square of ideal white noise has no ordinary pointwise meaning.

The split-step Fourier method alternates free-space diffraction, $U_h=e^{hD}$, with a random phase screen,For jointly Gaussian random fields, Gaussian phase averaging gives $\mathbb{E}{e}^{i{\varphi }_h}=e^{-\mathrm{Var}({\varphi }_h)/2}$. More generally, a product of $p$ fields and $q$ conjugate fields picks up $\exp (i{\sum }_a{ϵ}_a{\varphi }_h({\mathbf{z}}_a))$, with signs ${ϵ}_a=+1$ or $-1$. Its screen average is determined entirely by the wave phase covariance matrix. The deterministic diffraction step acts on each coordinate, with opposite signs on factors formed by complex conjugation. This is the basis of the field-moment equations.

There is an important closure qualification. A stationary Gaussian random field need not have independent longitudinal increments. For finite-correlation fluctuations, the unlinearized parabolic wave equation instead givesIn the weak-fluctuation model, the exact first-moment equation isThe last term cannot in general be replaced by a constant times $m$. An exact generic solution of the linearized model iswhere $\mathcal{T}$ denotes time ordering of the propagation operators. For the unlinearized model, add $i{k}_0{\mu }^2{W}^2/2$ inside the propagation generator. The covariance function of the medium is needed to evaluate this expression; One-point Gaussian distributions alone would not even determine the joint wave phase statistics.

For example, omit diffraction and take a longitudinal autocovariance function $R(s)=\mathbb{E}[W(s,\mathbf{z})W(0,\mathbf{z})]$. Direct Gaussian phase averaging givesEven here, unit-variance stationary Gaussian random fields with different $R$ give different answers. In the general problem the diffraction and multiplication operators do not commute, so this scalar attenuation cannot simply be multiplied by $U_x$ without an additional approximation.

The standard closed answer uses the Markov approximation for a random medium, made explicit in part ii. Let the longitudinally integrated covariance kernel beReplace the medium by longitudinal Gaussian white noise with this strength. A screen of thickness $h$ then has $\mathbb{E}[{\varphi }_h({\mathbf{z}}_a){\varphi }_h({\mathbf{z}}_b)]=k_0^2{\mu }^{2}hB({\mathbf{z}}_a-{\mathbf{z}}_b)$ and is independent of the incoming field. For the moment $M=\mathbb{E}{\prod }_a{E}_a^{({ϵ}_a)}$, expanding both steps to order $h$ givesIn particular, the coherent attenuation in a white-noise random medium and its Fresnel propagator solution areEquivalently, write $Wdx=d{\beta }_x$, where ${\beta }_x(\mathbf{z})$ is Brownian motion in $x$ with transverse covariance kernel $\mathbb{E}[d{\beta }_x(\mathbf{z})d{\beta }_x({\mathbf{z}}^{\prime })]=B(\mathbf{z}-{\mathbf{z}}^{\prime })dx$. The Stratonovich integral formulation is $dE=DEdx+i{k}_0\mu E\circ d{\beta }_x$. Its Itô integral form isThe mean of the Itô integral vanishes, independently confirming the attenuation drift. For $x>0$ in two transverse dimensions,Thus a unit plane wave has $m=e^{-\gamma x}$. For a general incident wave envelope, the Fresnel propagator supplies its spreading. The attenuation is redistribution between the coherent and diffuse wave fields, rather than wave absorption: for a field and its conjugate at the same point, the screen contribution in the second-moment equation cancels.

Use the leading far-field kernel from the previous part. DefineThen ${\varphi }_1={\int }_{D}KV$ and, because the incident intensity is one,To first order in weak contrast, $V=2\mu k_0^{2}W+O({\mu }^2)$. Its linearized part is a centered stationary Gaussian random field, with autocorrelation function of a random field $C_V(\mathbf{s})=4{\mu }^2{k}_0^4{C}_W(\mathbf{s})$. The exponential moment of a Gaussian linear functional therefore gives the usual Gaussian intensity in the first Rytov approximation:Equivalently the exponent is $\mathbb{E}|{\varphi }_1{|}^2+\mathrm{Re}\mathbb{E}{\varphi }_1^2$. The second term generally matters because a finite real random medium need not generate a circular complex Gaussian amplitude.

If the full stated potential $V=k_0^2(2\mu W+{\mu }^2{W}^2)$ is retained, it is not Gaussian and has mean $m_V={\mu }^2{k}_0^2$. To retain every term through order ${\mu }^2$ in the intensity of this first Rytov field, the cumulant expansion instead givesHere $C_V^{\mathrm{cov}}=\mathbb{E}[(V-m_V)(V^{\prime }-m_V)]$; the raw autocorrelation differs by $m_V^2=O({\mu }^4)$. The stated remainder uses Gaussian finite-dimensional statistics of $W$, for which odd total moments vanish. This mean term is absent if one first truncates the potential to $2\mu k_0^{2}W$. Both formulas concern the first Rytov approximation, and exponentiating it does not supply the omitted higher-order scattering corrections.

To first order in the weak fluctuation, $n^2-1=2\mu W+O({\mu }^2)$, soForintroduce the signs $s=(1,-1,1,-1)$. The product rule gives

Assume that $W$ is a zero-mean stationary Gaussian random field, that its longitudinal correlation length is short compared with the envelope's evolution scale, and that the propagation distance is long compared with that correlation length. The forward Markov approximation then neglects diffraction during one correlation length. Applying the Furutsu–Novikov formula closes the last average at second order in $\mu$. Define the integrated longitudinal autocorrelation function of a random fieldThenwhereWriting $C_{ij}=C(z_i-z_j)$ and using the evenness of the covariance gives the equivalent expressionThe derivation also assumes paraxial propagation, weak scattering, sufficient regularity to interchange differentiation and expectation, and statistical homogeneity in both coordinates. Without the short-correlation approximation, the Gaussian identity produces a nonlocal longitudinal memory integral rather than this local closed equation.

Since $n=1+W$ with small variance,To leading order, $V$ is therefore a centered stationary Gaussian random field. Let its autocorrelation function of a random field beThe $W^2$ contribution gives higher-order mean and non-Gaussian corrections and is consistently omitted at this order.

The far-field Rytov approximation from part ii has unit incident intensity andThe real random variable $X={\varphi }_1+{\varphi }_1^{\ast }$ is centered Gaussian. Its moment-generating function givesDefineThenand henceThis expression depends only on the two-point autocorrelation of the scattering potential. In the weak-fluctuation expansion it becomes