Publications

Export 142 results:
Sort by: Author Title Type [ Year  (Desc)]
2005
Colosi, JA, Baggeroer AB, Cornuelle BD, Dzieciuch MA, Munk WH, Worcester PF, Dushaw BD, Howe BM, Mercer JA, Spindel RC, Birdsall TG, Metzger K, Forbes AMG.  2005.  Analysis of multipath acoustic, field variability and coherence in the finale of broadband basin-scale transmissions in the North Pacific Ocean. Journal of the Acoustical Society of America. 117:1538-1564.   10.1121/1.1854615   AbstractWebsite

The statistics of low-frequency, long-range acoustic transmissions in the North Pacific Ocean are presented. Broadband signals at center frequencies of 28, 75, and 84 Hz are analyzed at propagation ranges of 3252 to 5171 km, and transmissions were received on 700 and 1400 in long vertical receiver arrays with 35 in hydrophone spacing. In the analysis we focus on the energetic "finale" region of the broadband time front arrival pattern, where a multipath interference pattern exists. A Fourier analysis of 1 s regions in the finale provide narrowband data for examination as well. Two-dimensional (depth and time) phase unwrapping is employed to study separately the complex field phase and intensity. Because data sampling occured in 20 or 40 min intervals followed by long gaps, the acoustic fields are analyzed. in terms of these 20 and 40 min and multiday observation times. An analysis of phase, intensity, and complex envelope variability as a function of depth and time is presented in terms of mean fields, variances, probability density functions (PDFs), covariance, spectra, and coherence. Observations are compared to a random multipath model of frequency and vertical wave number spectra for phase and log intensity, and the observations are compared to a broadband multipath model of scintillation index and coherence. 2005 Acoustical Society of America.

Vera, MD, Heaney KD, Grp N.  2005.  The effect of bottom interaction on transmissions from the North Pacific Acoustic Laboratory Kauai source. Journal of the Acoustical Society of America. 117:1624-1634.   10.1121/1.1854491   AbstractWebsite

Acoustic signals transmitted from a 75-Hz broadband source near Kauai as part of the North Pacific Acoustic Laboratory (NPAL), experiment were recorded on an array of receivers near California at a range of 3890 km, and on a vertical line array at a range of 3336 km in the Gulf of Alaska. Because the source is approximately 2 m above the seafloor, and the bottom depth at the receivers near California is approximately 1800 m, acoustic interaction with the bathymetry complicates the identification of the recorded arrivals with those present in numerical simulations of the experiment. Ray methods were used to categorize acoustic energy according to interactions with the sea bottom and surface and to examine the significance of seafloor geometry. A modal decomposition was also used to examine the role of range-dependent bathymetry and to associate the effects on the acoustic field with seafloor features at specific ranges. Parabolic-equation simulations were performed in order to investigate the sensitivity of the received signal to geoacoustic parameters; shear excitations within the seafloor were modeled using a complex-density, equivalent-fluid technique. Incorporation of bottom interaction into models of the propagation enables an identification between experimental and simulated arrivals. (c) 2005 Acoustical Society of America.

Andrew, RK, Howe BM, Mercer JA, Group NPAL, Cornuelle B, Colosi J.  2005.  Transverse horizontal spatial coherence of deep arrivals at megameter ranges. Journal of the Acoustical Society of America. 117:1511-1526.   10.1121/1.1854851   AbstractWebsite

Predictions of transverse horizontal spatial coherence from path integral theory are compared with measurements for two ranges between 2000 and 3000 km. The measurements derive from a low-frequency (75 Hz) bottom-mounted source at depth 810 m near Kauai that transmitted m-sequence signals over several years to two bottom-mounted horizontal line arrays in the North Pacific. In this paper we consider the early arriving portion of the deep acoustic field at these arrays. Horizontal coherence length estimates, on the order of 400 m, show good agreement with lengths calculated from theory. These lengths correspond to about 1 degrees in horizontal arrival angle variability using a simple, extended, spatially incoherent source model, Estimates of scintillation index, log-amplitude variance, and decibel intensity variance indicate that the fields were partially saturated. There was no significant seasonal variability in these measures. The scintillation index predictions agree quite well with the dataset estimates; nevertheless, the scattering regime predictions (fully saturated) vary from the regime classification (partially saturated) inferred from observation. This contradictory result suggests that a fuller characterization of scattering regime metrics may be required. (c) 2005 Acoustical Society of America.

Hoteit, I, Cornuelle B, Kohl A, Stammer D.  2005.  Treating strong adjoint sensitivities in tropical eddy-permitting variational data assimilation. Quarterly Journal of the Royal Meteorological Society. 131:3659-3682.: Royal Meteorological Society, 104 Oxford Rd. Reading Berks RG1 7LL UK   10.1256/qj.05.97   AbstractWebsite

A variational data assimilation system has been implemented for the tropical Pacific Ocean using an eddy-permitting regional implementation of the MITgcm. The adjoint assimilation system was developed by the Estimation of the Circulation and the Climate of the Ocean consortium, and has been extended to deal with open boundaries. This system is used to adjust the model to match observations in the tropical Pacific region using control parameters which include initial conditions, open boundaries and time-dependent surface fluxes. This paper focuses on problems related to strong adjoint sensitivities that may impede the model fit to the observations. A decomposition of the velocities at the open boundaries into barotropic and baroclinic modes is introduced to deal with very strong sensitivities of the model sea surface height to the barotropic component of the inflow. Increased viscosity and diffusivity terms are used in the adjoint model to reduce exponentially growing sensitivities in the backward run associated with nonlinearity of the forward model. Simplified experiments in which the model was constrained with Levitus temperature and salinity data, Reynolds sea surface temperature data and TOPEX/POSEIDON altimeter data were performed to demonstrate the controllability of this assimilation system and to study its sensitivity to the starting guesses for forcing and initial conditions.

2004
Skarsoulis, EK, Cornuelle BD.  2004.  Travel-time sensitivity kernels in ocean acoustic tomography. Journal of the Acoustical Society of America. 116:227-238.   10.1121/1.1753292   AbstractWebsite

Wave-theoretic ocean acoustic propagation modeling is combined with the peak arrival approach for tomographic travel-time observables to derive the sensitivity kernel of travel times with respect to sound-speed variations. This is the Born-Frechet kernel relating the three-dimensional spatial distribution of sound-speed variations with the induced travel-time variations. The derivation is based on the first Born approximation of the Green's function. The application of the travel-time sensitivity kernel to an ocean acoustic waveguide gives a picture close to the ray-theoretic one in the case of high frequencies. However, in the low-frequency case, of interest in ocean acoustic tomography, for example, there are significant deviations. Low-frequency travel times are sensitive to sound-speed changes in Fresnel-zone-scale areas surrounding the eigenrays, but not on the eigenrays themselves, where the sensitivity is zero. Further, there are areas of positive sensitivity, where, e.g., a sound-speed increase results in an increase of arrival times, i.e., a further delay of arrivals, in contrast with the common expectation. These findings are confirmed by forward acoustic predictions from a coupled-mode code. (C) 2004 Acoustical Society of America.

Hursky, P, Porter MB, Cornuelle BD, Hodgkiss WS, Kuperman WA.  2004.  Adjoint modeling for acoustic inversion. Journal of the Acoustical Society of America. 115:607-619.   10.1121/1.1636760   AbstractWebsite

The use of adjoint modeling for acoustic inversion is investigated. An adjoint model is derived from a linearized forward propagation model to propagate data-model misfit at the observation points back through the medium to the medium perturbations not being accounted for in the model. This adjoint model can be used to aid in inverting for these unaccounted medium perturbations. Adjoint methods are being applied to a variety of inversion problems, but have not drawn much attention from the underwater acoustic community. This paper presents an application of adjoint methods to acoustic inversion. Inversions are demonstrated in simulation for both range-independent and range-dependent sound speed profiles using the adjoint of a parabolic equation model. Sensitivity and error analyses are discussed showing how the adjoint model enables calculations to be performed in the space of observations, rather than the often much larger space of model parameters. Using an adjoint model enables directions of steepest descent in the model parameters (what we invert for) to be calculated using far fewer modeling runs than if a forward model only were used. (C) 2004 Acoustical Society of America.

Willis, JK, Roemmich D, Cornuelle B.  2004.  Interannual variability in upper ocean heat content, temperature, and thermosteric expansion on global scales. Journal of Geophysical Research-Oceans. 109   10.1029/2003jc002260   AbstractWebsite

[1] Satellite altimetric height was combined with approximately 1,000,000 in situ temperature profiles to produce global estimates of upper ocean heat content, temperature, and thermosteric sea level variability on interannual timescales. Maps of these quantities from mid-1993 through mid-2003 were calculated using the technique developed by Willis et al. [ 2003]. The time series of globally averaged heat content contains a small amount of interannual variability and implies an oceanic warming rate of 0.86 +/- 0.12 watts per square meter of ocean (0.29 +/- 0.04 pW) from 1993 to 2003 for the upper 750 m of the water column. As a result of the warming, thermosteric sea level rose at a rate of 1.6 +/- 0.3 mm/yr over the same time period. Maps of yearly heat content anomaly show patterns of warming commensurate with ENSO variability in the tropics, but also show that a large part of the trend in global, oceanic heat content is caused by regional warming at midlatitudes in the Southern Hemisphere. In addition to quantifying interannual variability on a global scale, this work illustrates the importance of maintaining continuously updated monitoring systems that provide global coverage of the world's oceans. Ongoing projects, such as the Jason/TOPEX series of satellite altimeters and the Argo float program, provide a critical foundation for characterizing variability on regional, basin, and global scales and quantifying the oceans' role as part of the climate system.

Di Lorenzo, E, Miller AJ, Neilson DJ, Cornuelle BD, Moisan JR.  2004.  Modelling observed California Current mesoscale eddies and the ecosystem response. International Journal of Remote Sensing. 25:1307-1312.   10.1080/01431160310001592229   AbstractWebsite

Satellite and in situ observations are used to test model dynamics for the California Current System (CCS). The model and data are combined to reconstruct the mesoscale ocean structure during a given three-week period. The resulting physical flow field is used to drive a 3D ecosystem model to interpret SeaWiFS and in situ chlorophyll-a (chl-a) variations. With this approach a more complete and consistent picture of the physical and ecosystem processes of the CCS is obtained, providing the basis for addressing fundamental questions about dynamics and predictability of the coastal ocean.

Moore, AM, Arango HG, Di Lorenzo E, Cornuelle BD, Miller AJ, Neilson DJ.  2004.  A comprehensive ocean prediction and analysis system based on the tangent linear and adjoint of a regional ocean model. Ocean Modelling. 7:227-258.   10.1016/j.ocemod.2003.11.001   AbstractWebsite

The regional ocean modelling system (ROMS) is a new generation ocean general circulation model that is rapidly gaining favour in the ocean modelling community. The tangent linear and adjoint versions of ROMS have recently been developed, and a new suite of tools that utilize these models for a variety of applications are now available to the ocean modelling community. In this paper we will describe the tangent linear and adjoint components of ROMS, and present examples from the tools that are currently available to ROMS users. In particular we will consider the finite time eigenmodes and the adjoint finite time eigenmodes of the tangent linear propagator, the singular vectors of the propagator, and its forcing singular vectors and stochastic optimals. The pseudospectra of the tangent linear resolvent matrix are also considered. Examples of each type of calculation will be presented for a time evolving double gyre ocean circulation in a rectangular ocean basin. (C) 2003 Elsevier Ltd. All rights reserved.

2003
Willis, JK, Roemmich D, Cornuelle B.  2003.  Combining altimetric height with broadscale profile data to estimate steric height, heat storage, subsurface temperature, and sea-surface temperature variability. Journal of Geophysical Research-Oceans. 108   10.1029/2002jc001755   AbstractWebsite

A new technique is demonstrated for combining altimetric height (AH) and sea-surface temperature (SST) with in situ data to produce improved estimates of 0/800 m steric height (SH), heat content, and temperature variability. The technique uses a linear regression onto AH to construct an initial guess for the subsurface quantity. This guess is then corrected toward the in situ data creating an estimate with substantially less error than could be achieved using either data set alone. Inclusion of the SST data further improves the estimates and illustrates how the procedure can be generalized to allow inclusion of additional data sets. The technique is demonstrated over a region in the southwestern Pacific enclosing the Tasman Sea. Nine-year time series of heat storage and temperature variability, averaged over 4degrees latitude and longitude and 1 year in time, are calculated. The estimates have RMS errors of approximately 4.6 W/m(2) in heat storage, 0.10degreesC in subsurface temperature and 0.11degreesC in surface temperature, and fractional errors of 20, 28, and 18%, respectively, relative to the total variance overall spatial and temporal scales considered. These represent significant improvements over previous estimates of these quantities. All the time series show strong interannual variability including the El Nino event of 1997. Application of these techniques on a global scale could provide new insight into the variability of the general circulation and heat budget of the upper ocean.

2002
Send, U, Worcester PF, Cornuelle BD, Tiemann CO, Baschek B.  2002.  Integral measurements of mass transport and heat content in the Strait of Gibraltar from acoustic transmissions. Deep-Sea Research Part Ii-Topical Studies in Oceanography. 49:4069-4095.   10.1016/s0967-0645(02)00143-1   AbstractWebsite

Three 2 kHz acoustic transceivers were deployed on either side of the eastern entrance of the Strait of Gibraltar during April-May 1996 to determine the feasibility of using acoustic methods to make routine, rapidly repeated, horizontally integrated measurements of flow and temperature in straits. Reciprocal transmissions between the transceivers were used to test the feasibility of using traditional ray differential travel times to monitor the component of flow along the acoustic paths. Transmissions directly across the Strait were used to test the feasibility of using horizontal arrival angle fluctuations and acoustic intensity scintillations to monitor the flow perpendicular to the acoustic path. The geometry was selected to provide ray paths that only sample the lower-layer Mediterranean water, so that the feasibility of monitoring the Mediterranean outflow using the various methods could be evaluated. The acoustic scintillation method did not yield useful current estimates, but the experimental parameters were not optimized for this approach. Since the low-frequency variability in log-amplitude was found to be highly correlated at receivers 228 m apart, it is possible that acoustic scintillation measurements using different receiver spacings and more rapid sampling might yield better results. The horizontal deflection method gave encouraging results at the time of neap tides, but less so during spring tides. For this approach, both theoretical estimates and measured phase differences between the horizontally separated receivers suggest that internal-wave-induced horizontal arrival angle fluctuations may fundamentally limit the precision with which arrival angles can be measured. Further work is needed to determine if a smaller horizontal spacing and higher signal-to-noise ratios would yield better results. Reciprocal travel time measurements diagonally across the Strait performed the best of the three methods, giving absolute flow estimates consistent with those derived from current-meter data. The fractional uncertainty variance for the lower layer tidal transport from a single tomographic path was estimated to be 0.017 (i.e. 98% of the a priori tidal transport variance was resolved). The spatial scales of the sub-tidal flow are thought to be significantly shorter than those of the tidal flow, however, which means that a more elaborate monitoring network is required to achieve the same performance for sub-tidal variability. Finally, sum travel times from the reciprocal transmissions were found to provide good measurements of the temperature and heat content in the lower layer. (C) 2002 Published by Elsevier Science Ltd.

2001
Roemmich, D, Gilson J, Cornuelle B, Weller R.  2001.  Mean and time-varying meridional transport of heat at the tropical subtropical boundary of the North Pacific Ocean. Journal of Geophysical Research-Oceans. 106:8957-8970.   10.1029/1999jc000150   AbstractWebsite

Ocean heat transport near the tropical/subtropical boundary of the North Pacific during 1993-1999 is described, including its mean and time variability. Twenty-eight trans-Pacific high-resolution expendable bathythermograph (XBT)/expendable conductivity-temperature-depth (XCTD) transects are used together with directly measured and operational wind estimates to calculate the geostrophic and Ekman transports. The mean heat transport across the XBT transect was 0.83 +/- 0.12 pW during the 7 year period. The large number of transects enables a stable estimate of the mean field to be made, with error bars based on the known variability. The North Pacific heat engine is a shallow meridional overturning circulation that includes warm Ekman and western boundary current components flowing northward, balanced by a southward flow of cool thermocline waters (including Subtropical Mode Waters). A near-balance of geostrophic and Ekman transports holds in an interannual sense as well as for the time mean. Interannual variability in geostrophic transport is strikingly similar to the pattern of central North Pacific sea level pressure variability (the North Pacific Index). The interannual range in heat transport was more than 0.4 pW during 1993-1999, with maximum northward values about 1 pW in early 1994 and early 1997. The ocean heat transport time series is similar to that of European Centre for Medium-Range Weather Forecasts air-sea heat flux integrated over the Pacific north of the XBT line. The repeating nature of the XBT/XCTD transects, with direct wind measurements, allows a substantial improvement over previous heat transport estimates based on one-time transects. A global system is envisioned for observing the time-varying ocean heat transport and its role in the Earth's heat budget and climate system.

Tiemann, CO, Worcester PF, Cornuelle BD.  2001.  Acoustic scattering by internal solitary waves in the Strait of Gibraltar. Journal of the Acoustical Society of America. 109:143-154.   10.1121/1.1329624   AbstractWebsite

High-freqnency underwater acoustic transmissions across the Strait of Gibraltar were used to examine acoustic scattering caused by the unique internal wave field in the Strait. Internal solitary waves of 100 m in amplitude propagate along the interface between an upper layer of Atlantic water and a lower layer of Mediterranean water. The interface is also strongly modulated by internal tides of comparable amplitude. As internal solitary waves cross the acoustic path, they cause sharp soundspeed gradients which intermittently refract acoustic rays away from normal sound channels. Internal tides vertically shift soundspeed profiles for additional travel time variability. Although the acoustic scattering is quite complicated, it is also surprisingly robust, making it a good candidate for modeling. Key features of the acoustic arrival pattern can be accounted for in some detail by a model description of the complex hydraulics in the Strait. (C) 2001 Acoustical Society of America. [DOI: 10.1121/1.1329624].

Tiemann, CO, Worcester PF, Cornuelle BD.  2001.  Acoustic remote sensing of internal solitary waves and internal tides in the Strait of Gibraltar. Journal of the Acoustical Society of America. 110:798-811.   10.1121/1.1382617   AbstractWebsite

High-frequency underwater acoustic transmissions across the Strait of Gibraltar are used to examine the feasibility of acoustically measuring several physical processes in the Strait, a difficult area to sample with conventional instruments. Internal undular bores propagate along the interface between an upper layer of Atlantic water and a lower layer of Mediterranean water. As they cross the acoustic path they are recognized by their scattering effects in the acoustic record. The time between internal bore crossings is influenced more by the tidal phase of the bore release at the Camarinal Sill than by variability in the bore's propagation time to the acoustic path. When internal bores were present, the acoustic arrival patterns could be classified as one of three types with different internal bore and internal tide amplitudes. The arrival types alternate during spring to neap tide transitions, suggesting that internal bore amplitude is not linearly related to tidal height. The sensitivity of acoustic observables to several physical parameters is investigated using a forward model, and a demonstration of inverse techniques provides estimates of several physical parameters from spring tidal cycles. (C) 2001 Acoustical Society of America.

Dushaw, B, Bold G, Chiu CS, Colosi J, Cornuelle B, Desaubies Y, Dzieciuch M, Forbes A, Gaillard F, Gavrilov A, Gould J, Howe BM, Lawrence M, Lynch J, Menemenlis D, Mercer J, Mikhalevsky PN, Munk W, Nakano I, Schott F, Send U, Spindel R, Terre T, Worcester P, Wunsch C.  2001.  Observing the ocean in the 2000's: A strategy for the role of acoustic tomography in ocean climate observation. Observing the Oceans in the 21st Century. ( Koblinsky C, Smith NR, Eds.).:391-418., Melbourne: Bureau of Meteorology Abstract
n/a
2000
Miller, AJ, Di Lorenzo E, Neilson DJ, Cornuelle BD, Moisan JR.  2000.  Modeling CalCOFI observations during El Nino: Fitting physics and biology. California Cooperative Oceanic Fisheries Investigations Reports. 41:87-97. AbstractWebsite

Surveys of temperature, salinity, and velocity from CalCOFI, altimetric measurements of sea level, and drifter observations of temperature and velocity during the 1997-98 El Nino are now being fit with an eddy-resolving ocean model of the Southern California Eight region to obtain dynamically consistent estimates of eddy variability. Skill is evaluated by the model-data mismatch (rms error) during the fitting interval and eventually by forecasting independent data. Preliminary results of fitting July 1997 physical fields are discussed. The physical fields are used to drive a three-dimensional NPZD-type model to be fit to subsurface chlorophyll a (chl a), nitrate, and bulk zooplankton from CalCOFI surveys, and surface chi a from SeaWiFS. Preliminary results of testing the ecosystem model in one-dimensional and three-dimensional form are discussed.

Cornuelle, BD, Chereskin TK, Niiler PP, Morris MY, Musgrave DL.  2000.  Observations and modeling of a California undercurrent eddy. Journal of Geophysical Research-Oceans. 105:1227-1243.   10.1029/1999jc900284   AbstractWebsite

A deep, nonlinear warm eddy advecting water that was also anomalously saltier, lower in oxygen, and higher in nutrients relative to surrounding waters was observed in moored current and temperature measurements and in hydrographic data obtained at a site similar to 400 km off the coast of northern California. The eddy was reproduced using a nonlinear quasi-geostrophic model, initialized by an iterative procedure using time series of 2-day averaged moored current measurements. The procedure demonstrates how a data assimilative technique synthesizes and enhances the resolution of a relatively sparse data set by incorporating time-dependence and model physics. The model forecast showed significant skill above persistence or climatology for 40 days. Our hypothesis, that the eddy was generated at the coast in winter and subsequently moved 400 km offshore by May, is consistent with the eddy movement diagnosed by the model and with the observations and coastal climatology. The model evolution significantly underpredicted the temperature anomaly in the eddy owing in part to unmodeled salinity compensation in trapped California Undercurrent water. Together, observations and model results show a stable nonlinear eddy in the California Current System that transported water and properties southwestward through the energetic eastern boundary region. Coherent features such as this one may be a mechanism for property transfer between the eddy-rich coastal zone and the eddy desert of the eastern North Pacific Ocean.

Pinkel, R, Munk W, Worcester P, Cornuelle BD, Rudnick D, Sherman J, Filloux JH, Dushaw BD, Howe BM, Sanford TB, Lee CM, Kunze E, Gregg MC, Miller JB, Merrifield MA, Luther DS, Firing E, Brainard R, Flament PJ, Chave AD, Moum JM, Caldwell DR, Levine MD, Boyd T, Egbert GD.  2000.  Ocean mixing studied near Hawaiian Ridge. Eos, Transactions American Geophysical Union. 81:545-553.   10.1029/EO081i046p00545-02   AbstractWebsite

The Hawaii Ocean Mixing Experiment (HOME) is a grassroots program to study turbulent mixing processes near the Hawaiian Ridge. The HOME is motivated by the desire to understand diffusive aspects of the advective-diffusive balance that mediates the general circulation of the oceans. HOME is focused on tidally driven mixing, given the ubiquity of the tide as a deep-sea energy source. As the sea surface cools at high latitude, surface waters sink. Subsidence rate is sufficient to fill the worlds ocean with cold bottom water in approximately 3,000 years. Diffusive processes that transfer heat into the abyssal ocean are required to maintain a steady-state thermal structure. An effective eddy diffusivity of order Kp=10−4 m2 s−1, 700 times the molecular diffusivity of heat, is necessary [Munk, 1966]. Such a diffusivity might be supported by either mechanical mixing (turbulent transport) or thermodynamic (so-called doubly diffusive) processes.

1999
Colosi, JA, Scheer EK, Flatte SM, Cornuelle BD, Dzieciuch MA, Munk WH, Worcester PF, Howe BM, Mercer JA, Spindel RC, Metzger K, Birdsall TG, Baggeroer AB.  1999.  Comparisons of measured and predicted acoustic fluctuations for a 3250-km propagation experiment in the eastern North Pacific Ocean. Journal of the Acoustical Society of America. 105:3202-3218.   10.1121/1.424650   AbstractWebsite

During the Acoustic Engineering Test (AET) of the Acoustic Thermometry of Ocean Climate (ATOC) program, acoustic signals were transmitted from a broadband source with 75-Hz center frequency to a 700-m-long vertical array of 20 hydrophones at a distance of 3252 km receptions occurred over a period of-six days. Each received pulse showed early identifiable timefronts, followed by about 2 s of highly variable energy. For the identifiable timefronts, observations of travel-time variance, average pulse shape, and the probability density function (PDF) of intensity are presented, and calculations of internal-wave contributions to those fluctuations are compared to the observations. Individual timefronts have rms travel time fluctuations of 11 to 19 ms, with time scales of less than 2 h. The pulse time spreads are between 0 and 5.3 ms rms, which suggest that internal-wave-induced travel-time biases are of the same magnitude. The PDFs of intensity for individual ray arrivals are compared to log-normal and exponential distributions. The observed PDFs are closer to the log-normal distribution, and variances of log intensity are between (3.1 dB)(2) (with a scintillation index of 0.74) for late-arriving timefronts and (2.0 dB)(2) (with a scintillation index of 0.2) for the earliest timefronts. Fluctuations of the pulse termination time of the transmissions are observed to be 22 ms rms. The intensity PDF of nonidentified peaks in the pulse crescendo are closer to a log-normal distribution than an exponential distribution, but a Kolmogorov-Smimov test rejects both distributions. The variance of the nonidentified peaks is (3.5 dB)(2) land the-scintillation index is 0.92. As a group, the observations suggest that the propagation is on the border of the unsaturated and partially saturated regimes. After improving the specification of the. ray weighting function, predictions of travel-time variance using the Garrett-Munk (GM) internal-wave spectrum at one-half the reference energy are in good agreement with the observations, and the one-half GM energy level compares well with XBT data taken along the transmission path. Predictions of pulse spread and wave propagation regime are in strong disagreement with the observations. Pulse time spread estimates are nearly two orders of magnitude too large, and Lambda-Phi methods for predicting the wave propagation regime predict full saturation. (C) 1999 Acoustical Society of America. [S0001-4966(99)04606-8].

Worcester, PF, Cornuelle BD, Dzieciuch MA, Munk WH, Howe BM, Mercer JA, Spindel RC, Colosi JA, Metzger K, Birdsall TG, Baggeroer AB.  1999.  A test of basin-scale acoustic thermometry using a large-aperture vertical array at 3250-km range in the eastern North Pacific Ocean. Journal of the Acoustical Society of America. 105:3185-3201.   10.1121/1.424649   AbstractWebsite

Broadband acoustic signals were transmitted during November 1994 from a 75-Hz source suspended near the depth of the sound-channel axis to a 700-m long vertical receiving array approximately 3250 km distant in the eastern North Pacific Ocean. The early part of the arrival pattern consists of raylike wave fronts that are resolvable, identifiable, and stable. The later part of the arrival pattern does not contain identifiable raylike arrivals, due to scattering from internal-wave-induced sound-speed fluctuations. The observed ray travel times differ from ray predictions based on the sound-speed field constructed using nearly concurrent temperature and salinity measurements by more than a priori variability estimates, suggesting that the equation used to compute sound speed requires refinement. The range-averaged ocean sound speed can be determined with an uncertainty of about 0.05 m/s from the observed ray travel times together with the time at which the near-axial acoustic reception ends, used as a surrogate for the group delay of adiabatic mode 1. The change in temperature over six days can be estimated with an uncertainty of about 0.006 degrees C. The sensitivity of the travel times to ocean variability is concentrated near the ocean surface and at the corresponding conjugate depths, because all of the resolved ray arrivals have upper turning depths within a few hundred meters of the surface. (C) 1999 Acoustical Society of America. [S0001-4966(99)04506-3].

Miller, AJ, Cornuelle BD.  1999.  Forecasts from fits of frontal fluctuations. Dynamics of Atmospheres and Oceans. 29:305-333.   10.1016/s0377-0265(99)00009-3   AbstractWebsite

A primitive equation ocean model is fit with strong constraints to non-synoptic hydrographic surveys in an unstable frontal current region, the Iceland-Faeroe Front. The model is first initialized from a time-independent objective analysis of non-synoptic data (spanning 2 to 6 days). A truncated set of eddy-scale basis functions is used to represent the initial error in temperature, salinity, and velocity. A series of model integrations, each perturbed with one basis function for one dependent variable in one layer, is used to determine the sensitivity to the objective-analysis initial state of the match to the non-synoptic hydrographic data. A new initial condition is then determined from a generalized inverse of the sensitivity matrix and the process is repeated to account for non-linearity. The method is first tested in 'identical twin' experiments to demonstrate the adequacy of the basis functions in representing initial condition error and the convergence of the method to the true solution. The approach is then applied to observations gathered in August 1993 in the Iceland-Faeroe Front. Model fits are successful in improving the match to the true data, leading to dynamically consistent evolution scenarios. However, the forecast skill (here defined as the variance of the model-data differences) of the model runs from the optimized initial condition is not superior to less sophisticated methods of initialization, probably due to inadequate initialization data. The limited verification data in the presence of strong frontal slopes may not be sufficient to establish Forecast skill, so that it must be judged subjectively or evaluated by other quantitative measures. (C) 1999 Elsevier Science B.V. All rights reserved.

Colosi, JA, Grp A.  1999.  A review of recent results on ocean acoustic wave propagation in random media: Basin scales. Ieee Journal of Oceanic Engineering. 24:138-155.   10.1109/48.757267   AbstractWebsite

Measurements of basin-scale acoustic transmissions made during the last four years by the Acoustic Thermometry of Ocean Climate (ATOC) program have allowed for the study of acoustic fluctuations of low-frequency pulse propagation at ranges of 1000 to 5000 km, Analysis of data from the ATOC Acoustic Engineering Test conducted in November 1994 has revealed new and unexpected results for the physics of ocean acoustic wave propagation in random media, In particular, use of traditional Lambda, Phi methods (using the Garrett-Munk (GM) internal wave model) to identify the wave propagation regime for early identifiable wavefronts predict the saturated regime, whereas observations of intensity probability density functions, intensity variance, and pulse time spread and wander suggest that the propagation is more likely near the border between the unsaturated and partially saturated regimes. Calculations of the diffraction parameter Lambda are very sensitive to the broad-band nature of the transmitted pulse, with CW calculations differing from a simplistic broad-band calculation by 10(3)! A simple model of pulse propagation using the Born approximation shows that CW and broad-band cases are sensitive to a random medium very differently and a theoretical description of broad-band effects for pulse propagation through a random media remains a fundamental unsolved problem in ocean acoustics. The observations show that, at 75-Hz center frequency, acoustic normal mode propagation is strongly nonadiabatic due to random media effects caused by internal waves. Simulations at a lower frequency of 28 Hz suggest that the first few modes might be treated adiabatically even in a random ocean. This raises the possibility of using modal techniques for ocean acoustic tomography, thereby increasing the vertical resolution of thermometry. Finally, the observation of unsaturated or partially saturated propagation for 75-Hz broad-band transmissions, like those of ATOC, suggests that ray-based tomography will be robust at basin-scales. This opens up the possibility of ray-based internal wave tomography using the observables of travel time variance, and vertical and temporal coherence, Using geometrical optics and the GM internal wave spectrum, internal wave tomography for an assortment of parameters of the GM model can be formulated in terms of a mixed linear/nonlinear inverse, This is a significant improvement upon a Monte Carlo approach presented in this paper which is used to infer average internal wave energies as a function of depth for the SLICE89 experiment. However, this Monte Carlo approach demonstrated, for the SLICE89 experiment, that the GM model failed to render a consistent inverse for acoustic energy which sampled the upper 100 m of the ocean, Until a new theory for the forward problem is advanced, internal wave tomography utilizing the signal from strong mode coupling can only be carried out using time-consuming Monte Carlo methods.

Dushaw, BD, Howe BM, Mercer JA, Spindel RC, Baggeroer AB, Menemenlis D, Wunsch C, Birdsall TG, Clark KMC, Colosi JA, Cornuelle BD, Dzieciuch M, Munk W.  1999.  Multimegameter-Range Acoustic Data Obtained by Bottom-Mounted Hydrophone Arrays for Measurement of Ocean Temperature. Ieee Journal of Oceanic Engineering. 24:202-214.   10.1109/48.757271   AbstractWebsite

Acoustic signals transmitted from the ATOC source on Pioneer Seamount off the coast of California have been received at various sites around the Pacific Basin since January 1996. We describe data obtained using bottom-mounted receivers, including U.S. Navy Sound Surveillance System arrays, at ranges up to 5 Mm from the Pioneer Seamount source. Stable identifiable ray arrivals are observed in several cases, but some receiving arrays are not well suited to detecting the direct ray arrivals. At 5-Mm range, travel-time variations at tidal frequencies (about 50 ms peak to peak) agree well with predicted values, providing verification of the acoustic measurements as well as the tidal model. On the longest and northernmost acoustic paths, the time series of resolved ray travel times show an annual cycle peak-to-peak variation of about 1 s and other fluctuations caused by natural oceanic variability. An annual cycle is not evident in travel times from shorter acoustic paths in the eastern Pacific, though only one realization of the annual cycle is available. The low-pass-filtered travel times are estimated to an accuracy of about 10 ms. This travel-time uncertainty corresponds to errors in range- and depth-averaged temperature of only a few millidegrees, while the annual peak-to-peak variation in temperature averaged horizontally over the acoustic path and vertically over the upper 1 km of ocean is up to 0.5 degree C.

1998
Gilson, J, Roemmich D, Cornuelle B, Fu LL.  1998.  Relationship of TOPEX/Poseidon altimetric height to steric height and circulation in the North Pacific. Journal of Geophysical Research-Oceans. 103:27947-27965.   10.1029/98jc01680   AbstractWebsite

TOPEX/Poseidon altimetric height is compared with 20 transpacific eddy-resolving realizations of steric height. The latter are calculated from temperature (expendable bathythermograph (XBT)) and salinity (expendable conductivity and temperature profiler (XCTD)) profiles along a precisely repeating ship track over a period of 5 years. The overall difference between steric height and altimetric height is 5.2 cm RMS. On long wavelengths (lambda < 500 km), the 3.5 cm RMS difference is due mainly to altimetric measurement errors but also has a component from steric variability deeper than the 800 m limit of the XBT. The data sets are very coherent in the long wavelength band, with coherence amplitude of 0.89. This band contains 65% of the total variance in steric height. On short wavelengths (lambda > 500 km), containing 17% of the steric height variance, the 3.0 cm RMS difference and lowered coherence are due to the sparse distribution of altimeter ground tracks along the XBT section. The 2.4 cm RMS difference in the basin-wide spatial mean appears to be due to fluctuations in bottom pressure. Differences between steric height and altimetric height increase near the western boundary, but data variance increases even more, and so the signal-to-noise ratio is highest in the western quarter of the transect. Basin-wide integrals of surface geostrophic transport from steric height and altimetric height are in reasonable agreement. The 1.9 x 10(4) m(2) s(-1) RMS difference is mainly because the interpolated altimetric height lacks spatial resolution across the narrow western boundary current. A linear regression is used to demonstrate the estimation of subsurface temperature from altimetric data. Errors diminish from 0.8 degrees C at 200 m to 0.3 degrees C at 400 m. Geostrophic volume transport, 0-800 m, shows agreement that is similar to surface transport, with 4.8 Sverdrup (Sv) (10(6) m(3) s(-1)) RMS difference. The combination of altimetric height with subsurface temperature and salinity profiling is a powerful tool for observing variability in circulation and transport of the upper ocean. The continuing need for appropriate subsurface data for verification and for statistical estimation is emphasized. This includes salinity measurements, which significantly reduce errors in specific volume and steric height.

Baggeroer, AB, Birdsall TG, Clark C, Colosi JA, Cornuelle BD, Costa D, Dushaw BD, Dzieciuch M, Forbes AMG, Hill C, Howe BM, Marshall J, Menemenlis D, Mercer JA, Metzger K, Munk W, Spindel RC, Stammer D, Worcester PF, Wunsch C.  1998.  Ocean climate change; comparison of acoustic tomography, satellite altimetry, and modeling. Science. 281:1327-1332., Washington, DC, United States (USA): American Association for the Advancement of Science, Washington, DC   10.1126/science.281.5381.1327   AbstractWebsite

Comparisons of gyre-scale acoustic and direct thermal measurements of heat content in the Pacific Ocean, satellite altimeter measurements of sea surface height, and results from a general circulation model show that only about half of the seasonal and year-to-year changes in sea level are attributable to thermal expansion. Interpreting climate change signals from fluctuations in sea level is therefore complicated. The annual cycle of heat flux is 150 ± 25 watts per square meter (peak-to-peak, corresponding to a 0.2°C vertically averaged temperature cycle); an interannual change of similar magnitude is also detected. Meteorological estimates of surface heat flux, if accurate, require a large seasonal cycle in the advective heat flux.