The effect of stochastic fluctuations of the interface boundaries has been
incorporated into the elastic coupled-mode equations for 2-D range-dependen
t media. We assume the medium to be characterized by some deterministic ran
ge-dependent layered structure superposed by stochastic boundary fluctuatio
ns. The deterministic range-dependent structure defines the reference struc
ture, and the reference structure with superposed stochastic fluctuations d
efines the true layer boundaries in the medium. The boundary conditions on
the true boundary are expanded about the reference boundary, and first-orde
r perturbation theory is applied to the boundary conditions, as well as to
the equation of motion. The O(1) system and the O(epsilon) system represent
wave propagation in the deterministic model and in the stochastic model, r
espectively, with a first-order perturbation of roughness. The solution of
the O(1) system is referred to as the primary field. The solutions to the O
(epsilon) system relate the coherent field and the scattered field, which a
re represented as the superposition of local modes multiplied by the stocha
stic modal amplitudes. Propagation of coupled-mode elastic waves in a 2-D d
eterministic range-dependent medium is represented by a unitary coupled-mod
e propagator. The evolution equation for the stochastic medium and the stoc
hastic coupling matrices, which acts to convert the coherent field to the s
cattered field, is derived. Enforcing energy conservation on the O(epsilon)
system leads to the requirement that the propagator for the stochastic med
ium must satisfy a Lippmann-Schwinger-type integral equation, whose solutio
n can be represented by a formal perturbation series for the multiply scatt
ered wavefields. The integral equations for the mean field and the covarian
ce of the field are also presented. These formal theoretical results are va
lid to all orders of multiple scattering. In the second part of the paper t
he Born approximation to the Lippmann-Schwinger integral equation is used t
o extract information about attenuation due to rough surface scattering and
in the design of an inverse problem for the boundary roughness variance an
d correlation length. By defining the modal scattering cross-section, the f
ormula for the scattering Q(s)(-1) from the Born approximation for the scat
tered field is derived. The modal scattering cross-section and scattering Q
(s) for a range-dependent model with stochastic roughness are computed for
both exponential and Gaussian correlation functions. Finally, a formula for
the power spectrum of the coherent field is derived from the Born approxim
ation, and an inversion for the roughness variance and correlation length i
s designed by power spectrum fitting.