Sound at 85 to 450 Hz propagating in approximately 80-m depth water from fixed sources to a joint horizontal/vertical line array (HLA/VLA) is analyzed. The data are from a continental shelf area east of Delaware Bay (USA) populated with tidally generated long-and short-wavelength internal waves. Sound paths are 19 km in the along-shore (along internal-wave crest) direction and 30 km in the cross-shore direction. Spatial statistics of HLA arrivals are computed as functions of beam steering angle and time. These include array gain, horizontally lagged spatial correlation function, and coherent beam power. These quantities vary widely in magnitude, and vary over a broad range of time scales. For example, correlation scale can change rapidly from forty to five wavelengths, and correlation-scale behavior is anisotropic. In addition, the vertical array can be used to predict correlation expected for adiabatic propagation with cylindrical symmetry, forming a benchmark. Observed variations are in concert with internal-wave activity. Temporal variations of three coherence measures, horizontal correlation length, array gain, and ratio of actual correlation length to predicted adiabatic-mode correlation length, are very strong, varying by almost a factor of ten as internal waves pass.

An alternating direction implicit (ADI) three-dimensional fluid parabolic equation solution method with enhanced accuracy is presented. The method uses a square-root Helmholtz operator splitting algorithm that retains cross-multiplied operator terms that have been previously neglected. With these higher-order cross terms, the valid angular range of the parabolic equation solution is improved. The method is tested for accuracy against an image solution in an idealized wedge problem. Computational efficiency improvements resulting from the ADI discretization are also discussed.

A tank experiment was carried out to investigate underwater sound propagation over an elastic bottom in flat and sloping configurations. The purpose of the experiment was to evaluate range-dependent propagation models with high-quality experimental data. The sea floor was modeled as an elastic medium by a polyvinyl chloride slab. The relatively high rigidity of the slab requires accounting for shear waves in this environment. Acoustic measurements were obtained along virtual arrays in the water column using a robotic apparatus. Elastic parabolic equation solutions are in excellent agreement with data.

Using a moored source and horizontal/vertical line array combination, horizontal coherence properties of high signal to noise ratio (≥20dB) 100–1600Hz signals have been measured. Internal waves in the area of the measurement created moving episodic sound-speed anomaly structures, influencing coherence length. Measured horizontal coherence scales for 100Hz ranged from 5 to 20 acoustic wavelengths, and were inversely related to the sound-speed anomaly strength. Horizontal field properties were compared with fields computed using modal decompositions of the vertical signals. The comparison allows azimuthal field coherence properties to be studied apart from normal-mode interference effects.

Recent improvements in the parabolic equation method are combined to extend this approach to a larger class of seismo-acoustics problems. The variable rotated parabolic equation [J. Acoust. Soc. Am. 120, 3534-3538 (2006)] handles a sloping fluid-solid interface at the ocean bottom. The single-scattering solution [J. Acoust. Soc. Am. 121, 808-813 (2007)] handles range dependence within elastic sediment layers. When these methods are implemented together, the parabolic equation method can be applied to problems involving variations in bathymetry and the thickness of sediment layers. The accuracy of the approach is demonstrated by comparing with finite-element solutions. The approach is applied to a complex scenario in a realistic environment.