Publications

Export 144 results:
Sort by: Author Title Type [ Year  (Asc)]
2000
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.

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.

2001
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
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.

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].

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.

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.

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.

2004
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.

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.

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.

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.

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.

2005
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.

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.

Schneider, N, Cornuelle BD.  2005.  The forcing of the Pacific decadal oscillation. Journal of Climate. 18:4355-4373.   10.1175/jcli3527.1   AbstractWebsite

The Pacific decadal oscillation (PDO), defined as the leading empirical orthogonal function of North Pacific sea surface temperature anomalies, is a widely used index for decadal variability. It is shown that the PDO can be recovered from a reconstruction of North Pacific sea surface temperature anomalies based on a first-order autoregressive model and forcing by variability of the Aleutian low, El Nino-Southern Oscillation (ENSO), and oceanic zonal advection anomalies in the Kuroshio-Oyashio Extension. The latter results from oceanic Rossby waves that are forced by North Pacific Ekman pumping. The SST response patterns to these processes are not orthogonal, and they determine the spatial characteristics of the PDO. The importance of the different forcing, processes is frequency dependent. At interannual time scales, forcing from ENSO and the Aleutian low determines the response in equal parts. At decadal time scales, zonal advection in the Kuroshio-Oyashio Extension, ENSO, and anomalies of the Aleutian low each account for similar amounts of the PDO variance. These, results support the hypothesis that the PDO is not a dynamical mode, but arises from the superposition of sea surface temperature fluctuations with different dynamical origins.

2007
Di Lorenzo, E, Moore AM, Arango HG, Cornuelle BD, Miller AJ, Powell B, Chua BS, Bennett AF.  2007.  Weak and strong constraint data assimilation in the inverse Regional Ocean Modeling System (ROMS): Development and application for a baroclinic coastal upwelling system. Ocean Modelling. 16:160-187.   10.1016/j.ocemod.2006.08.002   AbstractWebsite

We describe the development and preliminary application of the inverse Regional Ocean Modeling System (ROMS), a four dimensional variational (4DVAR) data assimilation system for high-resolution basin-wide and coastal oceanic flows. Inverse ROMS makes use of the recently developed perturbation tangent linear (TL), representer tangent linear (RP) and adjoint (AD) models to implement an indirect representer-based generalized inverse modeling system. This modeling framework is modular. The TL, RP and AD models are used as stand-alone sub-models within the Inverse Ocean Modeling (IOM) system described in [Chua, B.S., Bennett, A.F., 2001. An inverse ocean modeling system. Ocean Modell. 35 137-165.]. The system allows the assimilation of a wide range of observation types and uses an iterative algorithm to solve nonlinear assimilation problems. The assimilation is performed either under the perfect model assumption (strong constraint) or by also allowing for errors in the model dynamics (weak constraints). For the weak constraint case the TL and RP models are modified to include additional forcing terms on the right hand side of the model equations. These terms are needed to account for errors in the model dynamics. Inverse ROMS is tested in a realistic 3D baroclinic upwelling system with complex bottom topography, characterized by strong mesoscale eddy variability. We assimilate synthetic data for upper ocean (0-450 m) temperatures and currents over a period of 10 days using both a high resolution and a spatially and temporally aliased sampling array. During the assimilation period the flow field undergoes substantial changes from the initial state. This allows the inverse solution to extract the dynamically active information from the synthetic observations and improve the trajectory of the model state beyond the assimilation window. Both the strong and weak constraint assimilation experiments show forecast skill greater than persistence and climatology during the 10-20 days after the last observation is assimilated. Further investigation in the functional form of the model error covariance and in the use of the representer tangent linear model may lead to improvement in the forecast skill. (c) 2006 Elsevier Ltd. All rights reserved.

Miller, AJ, Neilson DJ, Luther DS, Hendershott MC, Cornuelle BD, Worcester PF, Dzieciuch MA, Dushaw BD, Howe BM, Levin JC, Arango HG, Haidvogel DB.  2007.  Barotropic Rossby wave radiation from a model Gulf Stream. Geophysical Research Letters. 34   10.1029/2007gl031937   AbstractWebsite

The barotropic Rossby wave field in the North Atlantic Ocean is studied in an eddy-resolving ocean model simulation. The meandering model Gulf Stream radiates barotropic Rossby waves southward through preferred corridors defined by topographic features. The smoother region between the Bermuda Rise and the mid-Atlantic Ridge is a particularly striking corridor of barotropic wave radiation in the 20-50 day period band. Barotropic Rossby waves are also preferentially excited at higher frequencies over the Bermuda Rise, suggesting resonant excitation of topographic Rossby normal modes. The prevalence of these radiated waves suggests that they may be an important energy sink for the equilibrium state of the Gulf Stream.

Kohl, A, Stammer D, Cornuelle B.  2007.  Interannual to decadal changes in the ECCO global synthesis. Journal of Physical Oceanography. 37:313-337.   10.1175/jpo3014.1   AbstractWebsite

An estimate of the time-varying global ocean circulation for the period 1992 - 2002 was obtained by combining most of the World Ocean Circulation Experiment ( WOCE) ocean datasets with a general circulation model on a 1 horizontal grid. The estimate exactly satisfies the model equations without artificial sources or sinks of momentum, heat, and freshwater. To bring the model into agreement with observations, its initial temperature and salinity conditions were permitted to change, as were the time-dependent surface fluxes of momentum, heat, and freshwater. The estimation of these "control variables" is largely consistent with accepted uncertainties in the hydrographic climatology and meteorological analyses. The estimated time-mean horizontal transports of volume, heat, and freshwater, which were largely underestimated in the previous 2 optimization performed by Stammer et al., have converged with time-independent estimates from box inversions over most parts of the World Ocean. Trends in the model's heat content are 7% larger than those reported by Levitus and correspond to a global net heat uptake of about 1.1 W m(-2) over the model domain. The associated model trend in sea surface height over the estimation period resembles the observations from Ocean Topography Experiment ( TOPEX)/Poseidon over most of the global ocean. Sea surface height changes in the model are primarily steric but show contributions from mass redistributions from the subpolar North Atlantic Ocean and the Southern Ocean to the subtropical Pacific Ocean gyres. Steric contributions are primarily temperature based but are partly compensated by salt variation. However, the North Atlantic and the Southern Ocean reveal a clear contribution of salt to large-scale sea level variations.

Kim, SY, Terrill E, Cornuelle B.  2007.  Objectively mapping HF radar-derived surface current data using measured and idealized data covariance matrices. Journal of Geophysical Research-Oceans. 112   10.1029/2006jc003756   AbstractWebsite

Surface currents measured by high-frequency radars are objectively mapped using covariance matrices computed from hourly surface current vectors spanning two years. Since retrievals of surface radial velocities are inherently gappy in space and time, the irregular density of surface current data leads to negative eigenvalues in the sample covariance matrix. The number and the magnitude of the negative eigenvalues depend on the degree of data continuity used in the matrix computation. In a region of 90% data coverage, the negative eigenvalues of the sample covariance matrix are small enough to be removed by adding a noise term to the diagonal of the matrix. The mapping is extended to regions of poorer data coverage by applying a smoothed covariance matrix obtained by spatially averaging the sample covariance matrix. This approach estimates a stable covariance matrix of surface currents for regions with the intermittent radar coverage. An additional benefit is the removal of baseline errors that often exist between two radar sites. The covariance matrices and the correlation functions of the surface currents are exponential in space rather than Gaussian, as is often assumed in the objective mapping of oceanographic data sets. Patterns in the decorrelation length scale provide the variabilities of surface currents and the insights on the influence of topographic features (bathymetry and headlands). The objective mapping approach presented herein lends itself to various applications, including the Lagrangian transport estimates, dynamic analysis through divergence and vorticity of current vectors, and statistical models of surface currents.

2008
Roux, P, Cornuelle BD, Kuperman WA, Hodgkiss WS.  2008.  The structure of raylike arrivals in a shallow-water waveguide. Journal of the Acoustical Society of America. 124:3430-3439.   10.1121/1.2996330   AbstractWebsite

Acoustic remote sensing of the oceans requires a detailed understanding of the acoustic forward problem. The results of a shallow-water transmission experiment between a vertical array of sources and a vertical array of receivers are reported. The source array is used to provide additional degrees of freedom to isolate and track raylike arrivals by beamforming over both source and receiver arrays. The coordinated source-receiver array processing procedure is presented and its effectiveness in an example of tracking raylike arrivals in a fluctuating ocean environment is shown. Many of these arrivals can be tracked over an hour or more and show slowly varying amplitude and phase. The use of a double-beamforming algorithm lays the foundation for shallow-water acoustic remote sensing using travel time and source and receive angles of selected eigenrays. (C) 2008 Acoustical Society of America. [DOI: 10.1121/1.2996330]

Hoteit, I, Cornuelle B, Thierry V, Stammer D.  2008.  Impact of resolution and optimized ECCO forcing on Simulations of the tropical pacific. Journal of Atmospheric and Oceanic Technology. 25:131-147.   10.1175/2007jtecho528.1   AbstractWebsite

The sensitivity of the dynamics of a tropical Pacific Massachusetts Institute of Technology (MIT) general circulation model (MITgcm) to the surface forcing fields and to the horizontal resolution is analyzed. During runs covering the period 1992-2002, two different sets of surface forcing boundary conditions are used, obtained 1) from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis project and 2) from the Estimating the Circulation and Climate of the Ocean (ECCO) assimilation consortium. The "ECCO forcing" is the "NCEP forcing" adjusted by a state estimation procedure using the MITgcm with a 1 degrees x 1 degrees global grid and the adjoint method assimilating a multivariate global ocean dataset. The skill of the model is evaluated against ocean observations available in situ and from satellites. The model domain is limited to the tropical Pacific, with open boundaries located along 26 degrees S, 26 degrees N, and in the Indonesian throughflow. To account for large-scale changes of the ocean circulation, the model is nested in the global time-varying ocean state provided by the ECCO consortium on a 1 grid. Increasing the spatial resolution to 1/3 degrees and using the ECCO forcing fields significantly improves many aspects of the circulation but produces overly strong currents in the western model domain. Increasing the resolution to 1/6 degrees does not yield further improvements of model results. Using the ECCO heat and freshwater fluxes in place of NCEP products leads to improved time-mean model skill (i.e., reduced biases) over most of the model domain, underlining the important role of adjusted heat and freshwater fluxes for improving model representations of the tropical Pacific. Combinations of ECCO and NCEP wind forcing fields can improve certain aspects of the model solutions, but neither ECCO nor NCEP winds show clear overall superiority.

Raghukumar, K, Cornuelle BD, Hodgkiss WS, Kuperman WA.  2008.  Pressure sensitivity kernels applied to time-reversal acoustics. Journal of the Acoustical Society of America. 124:98-112.   10.1121/1.2924130   AbstractWebsite

Sensitivity kernels for receptions of broadband sound transmissions are used to study the effect of the transmitted signal on the sensitivity of the reception to environmental perturbations. A first-order Born approximation is used to obtain the pressure sensitivity of the received signal to small changes in medium sound speed. The pressure perturbation to the received signal caused by medium sound speed changes is expressed as a linear combination of single-frequency sensitivity kernels weighted by the signal in the frequency domain. This formulation can be used to predict the response of a source transmission to sound speed perturbations. The stability of time-reversal is studied and compared to that of a one-way transmission using sensitivity kernels. In the absence of multipath, a reduction in pressure sensitivity using time reversal is only obtained with multiple sources. This can be attributed both to the presence of independent paths and to cancellations that occur due to the overlap of sensitivity kernels for different source-receiver paths. The sensitivity kernel is then optimized to give a new source transmission scheme that takes into account knowledge of the medium statistics and is related to the regularized inverse filter. (c) 2008 Acoustical Society of America.