Coupled Mode Problems for Bottom Interacting Sound and Coupled Mode Problems for Bottom Interacting Sound: Student Support (Assert) Robert I. Odom Applied Physics Laboratory College of Ocean and Fisheries Sciences University of Washington 1013 N.E. 40th Street phone: (206)685-3788 fax: (206)543-6785 email: [email protected] Award #’s: N00014-96-1-0233 & N00014-97-1-0675 (AASERT) LONG-TERM GOAL The long-term goal of this research is to improve our ability to model and predict VLF acoustic propagation in shallow water with particular emphasis on the range dependence of the medium and the geoacoustic properties of the bottom, and to quantify the various factors affecting the overall acoustic energy budget in shallow water propagation. OBJECTIVES Our scientific objectives are to incorporate the effects of sediment anisotropy, strong sediment attenuation, and the effects of both deterministic and stochastic medium properties into a local coupled mode propagation model, and to develop accurate theory and robust numerical algorithms for the shallow water propagation problem. APPROACH We are using an approach based on coupled local modes to carry out a systematic study of the effects of scattering, normal dispersion, anisotropy and intrinsic attenuation on a propagating shallow water acoustic signal with strong bottom interaction. The coupled mode theory is developed from the first order equations of motion for the stress and displacement rather than from the second order equations for a velocity or displacement potential. The later approach introduces coupling coefficients depending on the second-order derivatives with respect to the range coordinate of the local mode functions. These second-order coupling coefficients are an artifact of the formulation, and not present in the coupled mode theory based on the first order equations of motion. WORK COMPLETED We have derived the proper low p (slowness) mixed boundary condition for the spectral diffusion equation of Odom and Mercer (1996). We have incorporated a very general elastic mode code into our mode coupling program. The elastic mode code (Park, 1996) is capable of computing the elastic (seismo-acoustic) modes of a layered medium comprising anisotropic layers with hexagonal symmetry and arbitrary symmetry axis orientation. Two papers were completed and submitted to refereed journals (Park and Odom, 1998; Park and Odom, in press 1999). Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number. 1. REPORT DATE 3. DATES COVERED 1998 2. REPORT TYPE 00-00-1998 to 00-00-1998 4. TITLE AND SUBTITLE 5a. CONTRACT NUMBER Coupled Mode Problems for Bottom Interacting Sound and Coupled 5b. GRANT NUMBER Mode Problems for Bottom Interacting Sound: Student Support (Assert) 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR(S) 5d. PROJECT NUMBER 5e. TASK NUMBER 5f. WORK UNIT NUMBER 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION University of Washington,Applied Physics Laboratory,Seattle,WA,98105 REPORT NUMBER 9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR’S ACRONYM(S) 11. SPONSOR/MONITOR’S REPORT NUMBER(S) 12. DISTRIBUTION/AVAILABILITY STATEMENT Approved for public release; distribution unlimited 13. SUPPLEMENTARY NOTES See also ADM002252. 14. ABSTRACT 15. SUBJECT TERMS 16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF 18. NUMBER 19a. NAME OF ABSTRACT OF PAGES RESPONSIBLE PERSON a. REPORT b. ABSTRACT c. THIS PAGE Same as 4 unclassified unclassified unclassified Report (SAR) Standard Form 298 (Rev. 8-98) Prescribed by ANSI Std Z39-18 RESULTS Derivation of the low p mixed boundary condition for the spectral diffusion equation requires the use of improper modes representation of the elastic continuum. We employed the formulation of Maupin (1996), and have found that frequency dependent function multiplying the spectral energy flux varies as (1/frequency). This provides the proper total loss limit as frequency increases, and also indicates that the approximate boundary condition assuming total loss is likely to be a poor approximation at even moderate frequencies. Transverse isotropy (TI), which by definition has hexagonal symmetry and a vertical symmetry axis, is an almost ubiquitous feature of marine sediments. There is, however, no reason to expect that the symmetry axis remains vertical in regions with non-planar bottom topography. Bottom topography will cause the symmetry axis to rotate, and when this occurs, bottom interacting modes can develop significant quasi-SH (quasi-Love) particle motion components. Figure 1 is a plot of phase velocity vs. symmetry axis angle. Notice that for symmetry axis tilt of zero degrees (TI), the phase velocities of the quasi-SH (labeled qL) and the quasi-P-SV (labeled qR), are distinct and well separated. But that as the symmetry axis rotates, the qL and qR phase velocities can approach very closely. This has a dramatic effect on the mode coupling in a range dependent medium. Figure 2 is the coupling matrix for a symmetry axis tilt angle of 30 degrees. Whereas the mode coupling matrix for a range dependent isotropic or TI medium is strongly diagonally dominant indicating mostly nearest-neighbor interaction (Odom et al., 1996), the mode coupling matrix shown in Figure 2 has large terms far from the diagonal. This indicates coupling to SH in the bottom, and indicates that coupling to SH may be an important loss mechanism in regions of strong bottom anisotropy. The frequency in the example is 15 Hz. IMPACT/APPLICATIONS The development of additional code for more general anisotropic bottom structures and the derivation of the proper boundary condition for the spectral diffusion equation indicate steady progress towards meeting our scientific objectives. This research is directly applicable to predicting the effect of a complicated shallow water environment on the acoustic field. TRANSITIONS Modal methods for modeling in random range dependent shallow water waveguides should provide important constraints on the most significant waveguide properties affecting propagation at low frequencies. RELATED PROJECTS Our research is directly related to other programs studying surface, volume and bottom interaction effects, including 6.2 and 6.3 efforts to quantify bottom backscatter and bottom loss effects in littoral regions. REFERENCES Maupin, “The radiation modes of a vertically varying half-space: a new representation of the complete Green’s function in terms of modes,” Geophys J. Int., 126, 762-780, 1996. Odom, R.I., M. Park, J.A. Mercer, R.S. Crosson and P. Paik, “Effects of transverse isotropy on modes and mode coupling in shallow water,” J. Acoust. Soc. Am., 100, 2079-2092, 1996. Odom, R.I. and J.A. Mercer, “Spectral diffusion of elastic wave energy,” Geophys. J. Int., 127, 111- 124, 1996. Park, J., “Surface waves in layered anisotropic structures,” Geophys. J. Int., 126, 173-183, 1996. PUBLICATIONS Park, M. and R.I. Odom, “Effects of elastic heterogeneities and anisotropy on mode and coupling and signals in shallow water,” J. Acoust. Soc. Am, 104, 747-757, 1998. Park M. and R.I. Odom, “The effects of stochastic rough interfaces on coupled-mode elastic waves,” Geophys. J. Int., in press, 1998 or 1999. Figure 1. Modal phase velocities for a layered anisotropic hexagonally symmetric sediment model as a function of symmetry axis tilt angle. Note that the quasi-SH modes (qL) and quasi-P-SV (qR) modes attract and repel as the symmetry axis is rotated. Figure 2. Mode-mode coupling matrix for the same sediment model as in Figure 1 with a symmetry axis tilt of 30 degrees. The mode coupling exhibits significant non-nearest neighbor interactions among the modesi due to coupling between qL and qR modes. This is in contrast to the heterogeneity induced mode coupling in isotropic and transversely isotropic media in which the coupling is nearest neighbor dominant.