[Paper Review] Some applications of Rice formulas to waves

Jean‐Marc Azäis, José R. León+1 more|ArXiv.org|Oct 5, 2009

Ocean Waves and Remote Sensing13 references3 citations

TL;DR

This paper applies Rice's formulas to compute moments—particularly expectation and second moments—of geometric functionals in random wavefields, including specular points, twinkles, dislocations, and level set measures in oceanography and optics. It derives explicit formulas for these functionals under Gaussian random field models and establishes a Central Limit Theorem for certain functionals via higher-order moment computations.

ABSTRACT

We use Rice's formulas in order to compute the moments of some level functionals which are linked to problems in oceanography and optics. For instance, we consider the number of specular points in one or two dimensions, the number of twinkles, the distribution of normal angle of level curves and the number or the length of dislocations in random wavefronts. We compute expectations and in some cases, also second moments of such functionals. Moments of order greater than one are more involved, but one needs them whenever one wants to perform statistical inference on some parameters in the model or to test the model itself. In some case we are able to use these computations to obtain a Central Limit Theorem.

Motivation & Objective

• To compute the expected number and second moments of specular points in one- and two-dimensional random wavefronts.
• To analyze the distribution of twinkles—transient points of maximum reflection intensity—arising from time-varying wavefronts.
• To study dislocations in random wavefronts, defined as loci where wave amplitude vanishes, particularly in 2D and 3D settings.
• To derive explicit formulas for the expected Hausdorff measure of level sets of random fields, including higher-order moments.
• To establish conditions under which a Central Limit Theorem holds for wavefront geometric functionals using higher-order moment computations.

Proposed method

• Uses Rice’s formula to compute the expected number of zeros of a random field and its derivatives, particularly for level sets and critical points.
• Applies the Longuet-Higgins approximation to simplify the condition for specular points to $ W_x(x,t) = kx $, enabling tractable moment computation.
• Derives the joint characteristic function of the gradient and Hessian of a Gaussian random field to compute moments of dislocation and twinkling counts.
• Employs multivariate complex Gaussian field representations and moment-generating functions to compute integrals over the joint distribution of field and derivative values.
• Uses residue calculus and algebraic manipulation to evaluate high-dimensional integrals arising from characteristic function inversion.
• Transforms the problem into a form equivalent to known formulas in the literature (e.g., [6]), validating the derivation through comparison.

Experimental results

Research questions

• RQ1 What is the expected number of specular points in a random wavefront under the Longuet-Higgins approximation?
• RQ2 How can the second moment of the number of twinkles be computed in a time-dependent random wavefield?
• RQ3 What is the expected number of dislocations (phase singularities) in a 2D or 3D random wavefield with independent real and imaginary parts?
• RQ4 Can a Central Limit Theorem be established for the number of level set crossings or dislocations based on higher-order moment computations?
• RQ5 What is the expected Hausdorff measure of the level set $ \{ \mathbf{x} : W(\mathbf{x}) = \mathbf{u} \} $ for a Gaussian random field in $ \mathbb{R}^d $?

Key findings

• The expected number of specular points in a one-dimensional wavefield is computed explicitly using Rice’s formula under the Longuet-Higgins approximation.
• The second moment of the number of twinkles is derived via the joint distribution of the gradient and Hessian of the wavefield, enabling statistical inference on wave dynamics.
• For 2D random wavefronts with independent real and imaginary parts, the expected number of dislocations is shown to be finite and computable via the joint zero set of two Gaussian fields.
• The paper derives a closed-form expression for the characteristic function of the wavefield’s gradient and Hessian, which is used to compute the distribution of dislocations and twinkles.
• The derived formulas for the moment-generating function of the wavefield’s derivatives are shown to be equivalent to known results in the literature, validating the approach.
• A Central Limit Theorem is established for the number of level set crossings under suitable moment and mixing conditions, based on the computed second moments.