No abstract

Seismic images of oceanic thermohaline finestructure record vertical displacements from internal waves and turbulence over large sections at unprecedented horizontal resolution. Where reflections follow isopycnals, their displacements can be used to estimate levels of turbulence dissipation, by applying the Klymak-Moum slope spectrum method. However, many issues must be considered when using seismic images for estimating turbulence dissipation, especially sources of random and harmonic noise. This study examines the utility of seismic images for estimating turbulence dissipation in the ocean, using synthetic modeling and data from two field surveys, from the South China Sea and the eastern Pacific Ocean, including the first comparison of turbulence estimates from seismic images and from vertical shear. Realistic synthetic models that mimic the spectral characteristics of internal waves and turbulence show that reflector slope spectra accurately reproduce isopycnal slope spectra out to horizontal wavenumbers of ;0.04 cpm, corresponding to horizontal wavelengths of 25 m. Using seismic reflector slope spectra requires recognition and suppression of shot-generated harmonic noise and restriction of data to frequency bands with signal-to-noise ratios greater than about 4. Calculation of slope spectra directly from Fourier transforms of the seismic data is necessary to determine the suitability of a particular dataset to turbulence estimation from reflector slope spectra. Turbulence dissipation estimated from seismic reflector displacements compares well to those from 10-m shear determined by coincident expendable current profiler (XCP) data, demonstrating that seismic images can produce reliable estimates of turbulence dissipation in the ocean, provided that random noise is minimal and harmonic noise is removed.

One of the major problems in numerically simulating waves traveling in the Earth is that an artificial boundary must be introduced to produce unique solutions. To eliminate the spurious reflections introduced by this artificial boundary, we use a damping expression based on analogies to shock absorbers. This method can reduce the amplitude of the reflected wave to any prespecified value and is successful for waves at any angle of incidence. The method can eliminate unwanted reflections from the surface, reflections at the corners of the model, and waves reflected off an interface that strike the artificial boundary.Many of the boundary conditions currently used in the numerical solution of waves are approximations to perfectly absorbing boundary conditions and depend upon the angle of incidence of the incoming wave at the artificial boundary. Stability problems often occur with these boundary conditions. The method we use at the artificial boundary allows use of stable Dirichlet or von Neumann conditions.Since the surface motion is easily measured when the waves are induced by normal seismic techniques, approximations of surface waves are needed to obtain information on Rayleigh waves. An implicit finitedifference scheme that is relatively easy to incorporate into existing numerical simulators is used to obtain the surface data for the forward finite-difference approximations.

No abstract

No abstract