“…Previously, the two-way integral equation coupled-mode (IECM) approach has proved successful, when an adiabatic approach is not applicable, in producing benchmark solutions of the total pressure field as a function of range in the forward direction for large scale, i.e., smooth, range-dependent boundaries that result in forward-and back-scattering from mode coupling. [3][4][5][6] The goal of this paper is to establish a benchmark solution that is, in principle, an exact numerical solution for the reverberation time series for an environment with a range-dependent fine-scale rough bottom boundary that induces mode coupling and generates a scattered field. The solution includes scattering effects to all orders, i.e., sums the infinite series of forward and backward contributions at each range point and maintains energy conservation.…”