Efficient approach to computing travel time with the parabolic equation model.

An efficient method for computing the travel time of an acoustic wave using the parabolic equation model is presented. The frequency derivative of the acoustic phase is the differential travel time associated with a propagation in range. By taking this difference across closely spaced frequencies this method computes the acoustic travel time. This method requires the computation of the field at two frequencies rather than over the full band. The method compares well with other travel time methods for four different cases, including deep water, upslope and shallow water, and a three-dimensional propagation environment.

Computational Investigation of a Boundary-Layer Ingesting Propulsion System for the Common Research Model.

The present paper examines potential propulsive and aerodynamic benefits of integrating a Boundary-Layer Ingestion (BLI) propulsion system into the Common Research Model (CRM) geometry and the NASA Tetrahedral Unstructured Software System (TetrUSS). The Numerical Propulsion System Simulation (NPSS) environment is used to generate engine conditions for Computational Fluid Dynamics (CFD) analyses. Improvements to the BLI geometry are made using the Constrained Direct Iterative Surface Curvature (CDISC) design method. Potential benefits of the BLI system relating to cruise propulsive power are quantified using a power balance method, and a comparison to the baseline case is made. Iterations of the BLI geometric design are shown, and improvements between subsequent BLI designs are presented. Simulations are conducted for a cruise flight condition of Mach 0.85 at an altitude of 38,500 feet, with Reynolds number of 40 million based on mean aerodynamic chord and an angle of attack of 2° for all geometries. Results indicate an 8% reduction in engine power requirements at cruise for the BLI configuration compared to the baseline geometry. Small geometric alterations of the aft portion of the fuselage using CDISC has been shown to marginally increase the benefit from boundary-layer ingestion further, resulting in an 8.7% reduction in power requirements for cruise, as well as a drag reduction of approximately twelve counts over the baseline geometry.

Bathymetric diffraction of basin-scale hydroacoustic signals.

The ocean is nearly transparent to low frequency sound permitting the observation of distant events such as earthquakes or explosions at fully basin scales. For very low frequency the ocean acts as a shallow-water waveguide and lateral variability in bathymetry can lead to out-of-plane effects. In this paper, data from the International Monitoring System of the Comprehensive Test Ban Treaty Organization (CTBTO) is used to present two cases where robustly localized seismic events in locations clearly within the two-dimensional (2-D) shadow of a continent or large island generate T-phase signals that are received on a hydro-acoustic station. A fully three- dimensional parabolic equation model is used to demonstrate that lateral variability of the bathymetry can lead to diffraction, explaining both observations. The implications of this are that the CTBTO network has greater coverage than predicted by 2-D models and that inclusion of diffraction in future processing can improve the automatic global association of hydroacoustic events.

Three-dimensional parabolic equation modeling of mesoscale eddy deflection.

The impact of mesoscale oceanography, including ocean fronts and eddies, on global scale low-frequency acoustics is examined using a fully three-dimensional parabolic equation model. The narrowband acoustic signal, for frequencies from 2 to 16 Hz, is simulated from a seismic event on the Kerguellen Plateau in the South Indian Ocean to an array of receivers south of Ascension Island in the South Atlantic, a distance of 9100 km. The path was chosen for its relevance to seismic detections from the HA10 Ascension Island station of the International Monitoring System, for its lack of bathymetric interaction, and for the dynamic oceanography encountered as the sound passes the Cape of Good Hope. The acoustic field was propagated through two years (1992 and 1993) of the eddy-permitting ocean state estimation ECCO2 (Estimating the Circulation and Climate of the Ocean, Phase II) system. The range of deflection of the back-azimuth was 1.8° with a root-mean-square of 0.34°. The refraction due to mesoscale oceanography could therefore have significant impacts upon localization of distant low-frequency sources, such as seismic or nuclear test events.

Deep water towed array measurements at close range.

During the North Pacific Acoustic Laboratory Philippine Sea 2009 experiment, towed array receptions were made from a towed source as the two ships transited from a separation of several Convergence Zones through a Closest Point of Approach at 3 km. A combination of narrowband tones and broadband pulses were transmitted covering the frequency band 79-535 Hz. The received energy arrives from two general paths-direct path and bottom bounce. Bearing-time records of the narrowband arrivals at times show a 35° spread in the angle of arrival of the bottom bounce energy. Doppler processing of the tones shows significant frequency spread of the bottom bounce energy. Two-dimensional modeling using measured bathymetry, a geoacoustic parameterization based upon the geological record, and measured sound-speed field was performed. Inclusion of the effects of seafloor roughness and surface waves shows that in-plane scattering from rough interfaces can explain much of the observed spread in the arrivals. Evidence of out-of-plane scattering does exist, however, at short ranges. The amount of out-of-plane scattering is best observed in the broadband impulse-beam response analysis, which in-plane surface roughness modeling cannot explain.

Estimating uncertainty in subsurface glider position using transmissions from fixed acoustic tomography sources.

Four acoustic Seagliders were deployed in the Philippine Sea November 2010 to April 2011 in the vicinity of an acoustic tomography array. The gliders recorded over 2000 broadband transmissions at ranges up to 700 km from moored acoustic sources as they transited between mooring sites. The precision of glider positioning at the time of acoustic reception is important to resolve the fundamental ambiguity between position and sound speed. The Seagliders utilized GPS at the surface and a kinematic model below for positioning. The gliders were typically underwater for about 6.4 h, diving to depths of 1000 m and traveling on average 3.6 km during a dive. Measured acoustic arrival peaks were unambiguously associated with predicted ray arrivals. Statistics of travel-time offsets between received arrivals and acoustic predictions were used to estimate range uncertainty. Range (travel time) uncertainty between the source and the glider position from the kinematic model is estimated to be 639 m (426 ms) rms. Least-squares solutions for glider position estimated from acoustically derived ranges from 5 sources differed by 914 m rms from modeled positions, with estimated uncertainty of 106 m rms in horizontal position. Error analysis included 70 ms rms of uncertainty due to oceanic sound-speed variability.

Long range acoustic measurements of an undersea volcano.

A seamount 8 km southeast of Sarigan Island erupted on 29 May 2010 and was visually observed. The recordings on two sets of hydrophones, operated by International Monitoring System (IMS) of the Comprehensive Test Ban Treaty Organization (CTBTO) are analyzed. Each array is a triplet of axial single hydrophones deployed as a 2 km triangle. Measurements of acoustic intensity for the path to the southern triplet are on the order of 6 dB lower than those received on the northern triplet. Temporal cross-correlation beamforming estimation is performed and the estimated arrival angles for the two arrays, 265° and 267° were consistent with the predicted geodesic arrival of 264.6° and 267.8°, respectively. Cross-correlation between single phones on the northern and southern arrays reveals a peak at 266°, with a cross-correlation of 0.1. Nx2D parabolic equation modeling predicts complete blockage due to seamount interaction along the geodesic path. Overprediction of the seamount blockage indicates that the 2D approximation is incorrect, and three-dimensional propagation must be used to explain the observations. This is demonstrated by the computation of the Adiabatic Mode Parabolic Equation Transmission Loss, which predicts a 5-10 dB lower reception at the southern site.

Site specific probability of passive acoustic detection of humpback whale calls from single fixed hydrophones.

Passive acoustic monitoring of marine mammal calls is an increasingly important method for assessing population numbers, distribution, and behavior. A common mistake in the analysis of marine mammal acoustic data is formulating conclusions about these animals without first understanding how environmental properties such as bathymetry, sediment properties, water column sound speed, and ocean acoustic noise influence the detection and character of vocalizations in the acoustic data. The approach in this paper is to use Monte Carlo simulations with a full wave field acoustic propagation model to characterize the site specific probability of detection of six types of humpback whale calls at three passive acoustic monitoring locations off the California coast. Results show that the probability of detection can vary by factors greater than ten when comparing detections across locations, or comparing detections at the same location over time, due to environmental effects. Effects of uncertainties in the inputs to the propagation model are also quantified, and the model accuracy is assessed by comparing calling statistics amassed from 24,690 humpback units recorded in the month of October 2008. Under certain conditions, the probability of detection can be estimated with uncertainties sufficiently small to allow for accurate density estimates.

Comparison of hybrid three-dimensional modeling with measurements on the continental shelf.

During the CALOPS 2007 experiment, off the coast of Fort Lauderdale, Florida, three-dimensional (3D) multipath was observed using a bottom mounted horizontal line array during source tows along the 200 m isobath [Kevin D. Heaney and James J. Murray, J. Acoust. Soc. Am. 125(4), 1394-1402 (2008)]. In this paper a hybrid modeling approach is presented to model the 3D sound on the Florida shelf, nearly shaped like the canonical wedge. The hybrid approach combines vertical acoustic normal modes with the parabolic equation solution (in range/cross-range). The approach is shown to satisfy the 3D Cartesian-coordinate wave equation in the limit of adiabatic mode propagation. In the adiabatic mode parabolic equation (AMPE) approach modal phase speeds vs position are used as the input to the parabolic equation computation with dimensions of easting (km) and northing (km). Vertical adiabatic modes and horizontal rays are also computed to illustrate the 3D multipath arrival. The AMPE field is computed for all the modes for each element of the horizontal array. Beamforming vs source range is then conducted and excellent agreement with data is achieved.

Experimental observations of active invariance striations in a tank environment.

The waveguide invariant in shallow water environments has been widely studied in the context of passive sonar. The invariant provides a relationship between the frequency content of a moving broadband source and the distance to the receiver, and this relationship is not strongly affected by small perturbations in environment parameters such as sound speed or bottom features. Recent experiments in shallow water suggest that a similar range-frequency structure manifested as striations in the spectrogram exists for active sonar, and this property has the potential to enhance the performance of target tracking algorithms. Nevertheless, field experiments with active sonar have not been conclusive on how the invariant is affected by the scattering kernel of the target and the sonar configuration (monostatic vs bistatic). The experimental work presented in this paper addresses those issues by showing the active invariance for known scatterers under controlled conditions of bathymetry, sound speed profile and high SNR. Quantification of the results is achieved by introducing an automatic image processing approach inspired on the Hough transform for extraction of the invariant from spectrograms. Normal mode simulations are shown to be in agreement with the experimental results.