Publications

Export 151 results:
Sort by: Author Title [ Type  (Asc)] Year
Journal Article
Wolfe, CL, Cessi P, Cornuelle BD.  2017.  An intrinsic mode of interannual variability in the Indian Ocean. Journal of Physical Oceanography. 47:701-719.   10.1175/jpo-d-16-0177.1   AbstractWebsite

An intrinsic mode of self-sustained, interannual variability is identified in a coarse-resolution ocean model forced by an annually repeating atmospheric state. The variability has maximumloading in the Indian Ocean, with a significant projection into the South Atlantic Ocean. It is argued that this intrinsic mode is caused by baroclinic instability of the model's Leeuwin Current, which radiates out to the tropical Indian and South Atlantic Oceans as long Rossby waves at a period of 4 yr. This previously undescribed mode has a remarkably narrowband time series. However, the variability is not synchronized with the annual cycle; the phase of the oscillation varies chaotically on decadal time scales. The presence of this internal mode reduces the predictability of the ocean circulation by obscuring the response to forcing or initial condition perturbations. The signature of this mode can be seen in higher-resolution global ocean models driven by high-frequency atmospheric forcing, but altimeter and assimilation analyses do not show obvious signatures of such a mode, perhaps because of insufficient duration.

Muccino, JC, Arango HG, Bennett AF, Chua BS, Cornuelle BD, Di Lorenzo E, Egbert GD, Haidvogel D, Levin JC, Luo H, Miller AJ, Moore AA, Zaron ED.  2008.  The Inverse Ocean Modeling system. Part II: Applications. Journal of Atmospheric and Oceanic Technology. 25:1623-1637.   10.1175/2008jtecho522.1   AbstractWebsite

The Inverse Ocean Modeling (IOM) System is a modular system for constructing and running weak-constraint four-dimensional variational data assimilation (W4DVAR) for any linear or nonlinear functionally, smooth dynamical model and observing array. The IOM has been applied to four ocean models with widely varying characteristics. The Primitive Equations Z-coordinate-Harmonic Analysis of Tides (PEZ-HAT) and the Regional Ocean Modeling System (ROMS) are three-dimensional, primitive equations models while the Advanced Circulation model in 2D (ADCIRC-2D) and Spectral Element Ocean Model in 2D (SEOM-2D) are shallow-water models belonging to the general finite-element family. These models. in conjunction with the IOM, have been used to investigate a wide variety of scientific phenomena including tidal. mesoscale, and wind-driven circulation. In all cases, the assimilation of data using the IOM provides a better estimate of the ocean state than the model alone.

Hammer, PTC, Dorman LM, Hildebrand JA, Cornuelle BD.  1994.  Jasper Seamount Structure - Seafloor Seismic Refraction Tomography. Journal of Geophysical Research-Solid Earth. 99:6731-6752.   10.1029/93jb02170   AbstractWebsite

The velocity structure of Jasper Seamount was modeled using one- and three-dimensional inversions of P wave travel times. The results represent the first detailed seismic images of a submerged, intraplate volcano. Two seismic refraction experiments were completed on Jasper Seamount, incorporating ocean bottom seismometers and navigated seafloor shots. The P wave travel times were first used to compute a one-dimensional velocity profile which served as a starting model for a three-dimensional tomographic inversion. The seamount P velocities are significantly slower than those observed in typical oceanic crust at equivalent subbasement depths. This suggests that Jasper Seamount is constructed predominantly of extrusive lavas with high average porosity. The velocity models confirm morphological predictions: Jasper Seamount is a shield volcano with rift zone development. High seismic velocities were detected beneath the large radial ridges while low velocities characterize the shallow summit and flanks. Comparisons between P velocity models of Jasper Seamount and the island of Hawaii reveal that these two shield volcanoes are not structurally proportional. Jasper Seamount is far smaller than Hawaii, yet both volcanoes exhibit an outer extrusive layer of similar thickness. This suggests that seamount size influences the intrusive/extrusive proportions; density equilibrium between melt and country rock may explain this behavior.

Lien, RC, Ma B, Lee CM, Sanford TB, Mensah V, Centurioni LR, Cornuelle BD, Gopalakrishnan G, Gordon AL, Chang MH, Jayne SR, Yang YJ.  2015.  The Kuroshio and Luzon undercurrent east of Luzon Island. Oceanography. 28:54-63.   10.5670/oceanog.2015.81   AbstractWebsite

Current structure, transport, and water mass properties of the northward-flowing Kuroshio and the southward-flowing Luzon Undercurrent (LU) were observed for nearly one year, June 8, 2012-June 4, 2013, across the Kuroshio path at 18.75 degrees N. Observations were made from four platforms: an array of six subsurface ADCP moorings, two Seagliders, fivepressure inverted echo sounders (PIES), and five horizontal electric field (HEF) sensors, providing the most detailed time series of the Kuroshio and Luzon Undercurrent water properties to date. Ocean state estimates of the western boundary current system were performed using the MIT general circulation model-four-dimensional variational assimilation (MITgcm-4D-Var) system. Prominent Kuroshio features from observations are simulated well by the numerical model. Annual mean Kuroshio transport, averaged over all platforms, is similar to 16 Sv with a standard deviation similar to 4 Sv. Kuroshio and LU transports and water mass pathways east of Luzon are revealed by Seaglider measurements. In a layer above the salinity maximum associated with North Pacific Tropical Water (NPTW), Kuroshio transport is similar to 7 Sv and contains North Equatorial Current (NEC) and Western Philippine Sea (WPS) waters, with an insignificant amount of South China Sea water on the shallow western flank. In an intermediate layer containing the core of the NPTW, Kuroshio transport is similar to 10 Sv, consisting mostly of NEC water. In the lower layer of the Kuroshio, transport is similar to 1.5 Sv of mostly North Pacific Intermediate Water (NPIW) as a part of WPS waters. Annual mean Luzon Undercurrent southward transport integrated to 1,000 m depth is similar to 2.7 Sv with a standard deviation similar to 2 Sv, carrying solely WPS waters below the salinity minimum of the NPIW. The transport of the western boundary current integrated over the full ocean depth east of Luzon Island is similar to 14 +/- 4.5 Sv. Sources of the water masses in the Kuroshio and Luzon Undercurrent are confirmed qualitatively by the numerical model.

Subramanian, AC, Hoteit I, Cornuelle B, Miller AJ, Song H.  2012.  Linear versus Nonlinear Filtering with Scale-Selective Corrections for Balanced Dynamics in a Simple Atmospheric Model. Journal of the Atmospheric Sciences. 69:3405-3419.   10.1175/JAS-D-11-0332.1   AbstractWebsite

This paper investigates the role of the linear analysis step of the ensemble Kalman filters (EnKF) in disrupting the balanced dynamics in a simple atmospheric model and compares it to a fully nonlinear particle-based filter (PF). The filters have a very similar forecast step but the analysis step of the PF solves the full Bayesian filtering problem while the EnKF analysis only applies to Gaussian distributions. The EnKF is compared to two flavors of the particle filter with different sampling strategies, the sequential importance resampling filter (SIRF) and the sequential kernel resampling filter (SKRF). The model admits a chaotic vortical mode coupled to a comparatively fast gravity wave mode. It can also be configured either to evolve on a so-called slow manifold, where the fast motion is suppressed, or such that the fast-varying variables are diagnosed from the slow-varying variables as slaved modes. Identical twin experiments show that EnKF and PF capture the variables on the slow manifold well as the dynamics is very stable. PFs, especially the SKRF, capture slaved modes better than the EnKF, implying that a full Bayesian analysis estimates the nonlinear model variables better. The PFs perform significantly better in the fully coupled nonlinear model where fast and slow variables modulate each other. This suggests that the analysis step in the PFs maintains the balance in both variables much better than the EnKF. It is also shown that increasing the ensemble size generally improves the performance of the PFs but has less impact on the EnKF after a sufficient number of members have been used.

Orsi, AJ, Cornuelle BD, Severinghaus JP.  2012.  Little Ice Age cold interval in West Antarctica: Evidence from borehole temperature at the West Antarctic Ice Sheet (WAIS) Divide. Geophysical Research Letters. 39   10.1029/2012gl051260   AbstractWebsite

The largest climate anomaly of the last 1000 years in the Northern Hemisphere was the Little Ice Age (LIA) from 1400-1850 C. E., but little is known about the signature of this event in the Southern Hemisphere, especially in Antarctica. We present temperature data from a 300 m borehole at the West Antarctic Ice Sheet (WAIS) Divide. Results show that WAIS Divide was colder than the last 1000-year average from 1300 to 1800 C.E. The temperature in the time period 1400-1800 C.E. was on average 0.52 +/- 0.28 degrees C colder than the last 100-year average. This amplitude is about half of that seen at Greenland Summit (GRIP). This result is consistent with the idea that the LIA was a global event, probably caused by a change in solar and volcanic forcing, and was not simply a seesaw-type redistribution of heat between the hemispheres as would be predicted by some ocean-circulation hypotheses. The difference in the magnitude of the LIA between Greenland and West Antarctica suggests that the feedbacks amplifying the radiative forcing may not operate in the same way in both regions. Citation: Orsi, A. J., B. D. Cornuelle, and J. P. Severinghaus (2012), Little Ice Age cold interval in West Antarctica: Evidence from borehole temperature at the West Antarctic Ice Sheet (WAIS) Divide, Geophys. Res. Lett., 39, L09710, doi: 10.1029/2012GL051260.

Skarsoulis, EK, Cornuelle BD, Dzieciuch MA.  2013.  Long-range asymptotic behavior of vertical travel-time sensitivity kernels. Journal of the Acoustical Society of America. 134:3201-3210.   10.1121/1.4818785   AbstractWebsite

Vertical travel-time sensitivity kernels (VTSKs) describe the effect of horizontally uniform sound-speed changes on travel times in range-independent ocean environments. Wave-theoretic VTSKs can be obtained either analytically, through perturbation of the normal-mode representation, or numerically, as horizontal marginals of the corresponding two-dimensional and three-dimensional travel-time sensitivity kernels. In previous works, it has been observed that wave-theoretic finite-frequency VTSKs approach the corresponding ray-theoretic sensitivity kernels as the propagation range increases. The present work is an attempt to explain this behavior. A stationary-phase approach is used to obtain a long-range asymptotic expression for the wave-theoretic VTSKs. The resulting asymptotic VTSKs are very close to the corresponding ray-theoretic ones. The smoothness condition, required for the stationary-phase approximation to hold, is used to obtain an estimate for the range beyond which the asymptotic behavior sets in. (C) 2013 Acoustical Society of America.

Orsi, AJ, Cornuelle BD, Severinghaus JP.  2014.  Magnitude and temporal evolution of Dansgaard-Oeschger event 8 abrupt temperature change inferred from nitrogen and argon isotopes in GISP2 ice using a new least-squares inversion. Earth and Planetary Science Letters. 395:81-90.   10.1016/j.epsl.2014.03.030   AbstractWebsite

Polar temperature is often inferred from water isotopes in ice cores. However, non-temperature effects on 3180 are important during the abrupt events of the last glacial period, such as changes in the seasonality of precipitation, the northward movement of the storm track, and the increase in accumulation. These effects complicate the interpretation of 8180 as a temperature proxy. Here, we present an independent surface temperature reconstruction, which allows us to test the relationship between delta O-18(ice) and temperature, during Dansgaard-Oeschger event 8, 38.2 thousand yrs ago using new delta N-15 and delta Ar-40 data from the GISP2 ice core in Greenland. This temperature reconstruction relies on a new inversion of inert gas isotope data using generalized least-squares, and includes a robust uncertainty estimation. We find that both temperature and delta O-18 increased in two steps of 20 and 140 yrs, with an overall amplitude of 11.80 +/- 1.8 degrees C between the stadial and interstadial centennial-mean temperature. The coefficient alpha = d delta O-18/dT changes with each time-segment, which shows that non-temperature sources of fractionation have a significant contribution to the delta O-18 signal. When measured on century-averaged values, we find that alpha = d delta O-18/dT = 0.32 +/- 0.06%(0)/degrees C, which is similar to the glacial/Holocene value of 0.328%(o)/degrees C. (C) 2014 Elsevier B.V. All rights reserved.

Kim, SY, Terrill EJ, Cornuelle BD.  2008.  Mapping surface currents from HF radar radial velocity measurements using optimal interpolation. Journal of Geophysical Research-Oceans. 113   10.1029/2007jc004244   AbstractWebsite

An optimal interpolation (OI) method to compute surface vector current fields from radial velocity measurements derived from high-frequency (HF) radars is presented. The method assumes a smooth spatial covariance relationship between neighboring vector currents, in contrast to the more commonly used un-weighted least-squares fitting (UWLS) method, which assumes a constant vector velocity within a defined search radius. This OI method can directly compute any quantities linearly related to the radial velocities, such as vector currents and dynamic quantities (divergence and vorticity) as well as the uncertainties of those respective fields. The OI method is found to be more stable than the UWLS method and reduces spurious vector solutions near the baselines between HF radar installations. The OI method produces a covariance of the uncertainty of the estimated vector current fields. Three nondimensional uncertainty indices are introduced to characterize the uncertainty of the vector current at a point, representing an ellipse with directional characteristics. The vector current estimation using the OI method eliminates the need for multiple mapping steps and optimally fills intermittent coverage gaps. The effects of angular interpolation of radial velocities, a commonly used step in the preprocessing of radial velocity data prior to vector current computation in the UWLS method, are presented.

Kim, SY, Terrill EJ, Cornuelle BD, Jones B, Washburn L, Moline MA, Paduan JD, Garfield N, Largier JL, Crawford G, Kosro PM.  2011.  Mapping the U.S. West Coast surface circulation: A multiyear analysis of high-frequency radar observations. Journal of Geophysical Research-Oceans. 116   10.1029/2010jc006669   AbstractWebsite

The nearly completed U. S. West Coast (USWC) high-frequency radar (HFR) network provides an unprecedented capability to monitor and understand coastal ocean dynamics and phenomenology through hourly surface current measurements at up to 1 km resolution. The dynamics of the surface currents off the USWC are governed by tides, winds, Coriolis force, low-frequency pressure gradients (less than 0.4 cycles per day (cpd)), and nonlinear interactions of those forces. Alongshore surface currents show poleward propagating signals with phase speeds of O(10) and O(100 to 300) km day(-1) and time scales of 2 to 3 weeks. The signals with slow phase speed are only observed in southern California. It is hypothesized that they are scattered and reflected by shoreline curvature and bathymetry change and do not penetrate north of Point Conception. The seasonal transition of alongshore surface circulation forced by upwelling-favorable winds and their relaxation is captured in fine detail. Submesoscale eddies, identified using flow geometry, have Rossby numbers of 0.1 to 3, diameters in the range of 10 to 60 km, and persistence for 2 to 12 days. The HFR surface currents resolve coastal surface ocean variability continuously across scales from submesoscale to mesoscale (O(1) km to O(1000) km). Their spectra decay with k(-2) at high wave number (less than 100 km) in agreement with theoretical submesoscale spectra below the observational limits of present-day satellite altimeters.

Cornuelle, BD, Malanotterizzoli P.  1986.  A Maximum-Gradient Inverse for the Gulf-Stream System. Journal of Geophysical Research-Oceans. 91:566-580.   10.1029/JC091iC09p10566   AbstractWebsite

Acoustic tomography uses integrating measurements which require inverse methods to resolve the averages into estimates of spatial structure. Statistical inverse methods have been extensively used to solve the reconstruction problem over different tomographic ranges and configurations. These inverses become very difficult to apply in frontal regions like the Gulf Stream (GS) system, where the statistics are acutely inhomogeneous and anisotropic and the mean is not a likely representation of the GS front at any time. In this paper we propose an alternative inverse which asks for the solution which gives a front instead of asking for the smoothest solution. The inverse solution minimizes the errors in the fit to the data while simultaneously maximizing the sum of the squares of the gradients observed in the reconstructed section and minimizing the absolute value norm for stability. The inverse is aimed at detecting changes in the GS front, thus the data are used to estimate the perturbations to a previous estimate of the frontal structure, instead of reconstructing the entire front as a perturbation from some average state. This approach is intended to merge well with eventual dynamic updating schemes and can be used with various types of data, given a proper model. Several examples have been run intercomparing the traditional linear least squares (LLSI) with the maximum gradient inverse (MGI), from very idealized cases to a real Gulf Stream section reconstructed from hydrographic data. Different transceiver configurations were also compared and mid-depth instruments were found to be superior to bottom mounted instruments. The simplest cases show a significant improvement in the estimate of the Gulf Stream front by the MGI compared to the weighted least squares inverse (LLSI). As the cases became more complicated (and more realistic), the differences between inverse methods become less pronounced, although the strength and location of the perturbation maxima were always determined more accurately by the MGI. The decline is at least partially due to the numerical algorithm which lumps data misfits and external constraints (the maximum gradient) into a single penalty criterion which is minimized. The most immediate way to overcome this limitation is to break up the problem into a two-step procedure, first a least squares inverse to fit the data and second an iterative, nonlinear optimization maximizing the gradient and minimizing the absolute value norm.

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.

Duda, TF, Flatte SM, Colosi JA, Cornuelle BD, Hildebrand JA, Hodgkiss WS, Worcester PF, Howe BM, Mercer JA, Spindel RC.  1992.  Measured Wave-Front Fluctuations in 1000-Km Pulse-Propagation in the Pacific-Ocean. Journal of the Acoustical Society of America. 92:939-955.   10.1121/1.403964   AbstractWebsite

A 1000-km acoustical transmission experiment has been carried out in the North Pacific, with Pulses broadcast between a moored broadband source (250-Hz center frequency) and a moored sparse vertical line of receivers. Two data records are reported: a period of 9 days at a pulse rate of one per hour, and a 21 -h period on the seventh day at six per hour. Many wave-front segments were observed at each hydrophone depth, and arrival times were tracked and studied as functions of time and depth. Arrivals within the final section of the pulse are not trackable in time or space at the chosen sampling rates, however. Broadband fluctuations, which are uncorrelated over 10-min sampling and 60-m vertical spacing, are observed with about 40 (ms)2 variance. The variance of all other fluctuations (denoted as low-frequency) is comparable or smaller than the broadband value; this low-frequency variance can be separated into two parts: a wave-front segment displacement (with vertical correlation length greater than 1 km) that varies substantially between rays with different ray identifiers, and a distortion (with vertical correlation length between 60 m and 1 km) of about 2 (ms)2 variance. The low-frequency variance may be explained as the effect of internal waves, including internal tides. The variance of the broadband fluctuations is reduced somewhat but not eliminated if only high-intensity peaks are selected; this selection does not affect the statistics of the low-frequency fluctuations.

Schonau, MC, Rudnick DL, Cerovecki I, Gopalakrishnan G, Cornuelle BD, McClean JL, Qiu B.  2015.  The Mindanao Current mean structure and connectivity. Oceanography. 28:34-45.   10.5670/oceanog.2015.79   AbstractWebsite

The Mindanao Current (MC), a low-latitude western boundary current in the Pacific Ocean, plays an important role in heat and freshwater transport to the western Pacific warm pool and the Indian Ocean. However, there have been relatively few comprehensive studies of the structure and variability of the MC and its connectivity to regional circulation. The Origins of the Kuroshio and Mindanao Current (OKMC) initiative combines four years of glider observations of the MC, a historical conductivity-temperature-depth (CTD)/float climatology, and results from a global strongly eddying forward ocean general circulation model simulation and a regional ocean state estimate. The MC is resolved as a strong southward current primarily within the upper 200 m, approaching 1 m s(-1), and extending roughly 300 km offshore of Mindanao. Observations and model simulations show a persistent northward Mindanao Undercurrent (MUC) below the thermocline. The MC transports water masses of North Pacific origin southward, while the MUC carries water with South Pacific characteristics northward. The subthermocline transport of the MC and the MUC is connected to other undercurrents in the Philippine Sea. The variability of this transport is a topic of continuing research.

Hoteit, I, Hoar T, Gopalakrishnan G, Collins N, Anderson J, Cornuelle B, Kohl A, Heimbach P.  2013.  A MITgcm/DART ensemble analysis and prediction system with application to the Gulf of Mexico. Dynamics of Atmospheres and Oceans. 63:1-23.   10.1016/j.dynatmoce.2013.03.002   AbstractWebsite

This paper describes the development of an advanced ensemble Kalman filter (EnKF)-based ocean data assimilation system for prediction of the evolution of the loop current in the Gulf of Mexico (GoM). The system integrates the Data Assimilation Research Testbed (DART) assimilation package with the Massachusetts Institute of Technology ocean general circulation model (MITgcm). The MITgcm/DART system supports the assimilation of a wide range of ocean observations and uses an ensemble approach to solve the nonlinear assimilation problems. The GoM prediction system was implemented with an eddy-resolving 1/10th degree configuration of the MITgcm. Assimilation experiments were performed over a 6-month period between May and October during a strong loop current event in 1999. The model was sequentially constrained with weekly satellite sea surface temperature and altimetry data. Experiments results suggest that the ensemble-based assimilation system shows a high predictive skill in the GoM, with estimated ensemble spread mainly concentrated around the front of the loop current. Further analysis of the system estimates demonstrates that the ensemble assimilation accurately reproduces the observed features without imposing any negative impact on the dynamical balance of the system. Results from sensitivity experiments with respect to the ensemble filter parameters are also presented and discussed. (C) 2013 Elsevier B.V. All rights reserved.

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.

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.

The AMODE-MST Group, Birdsall TG, Boyd J, Cornuelle BD, Howe BM, Knox R, Mercer JA, Metzger Jr. K, Spindel RC, Worcester PF.  1994.  Moving ship tomography in the North Atlantic. EOS Trans. AGU. 75:17,21,23. Abstract
n/a
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.

Yoo, JG, Kim SY, Cornuelle BD, Kosro PM, Kurapov AL.  2017.  A Noninterpolated Estimate of Horizontal Spatial Covariance from Nonorthogonally and Irregularly Sampled Scalar Velocities. Journal of Atmospheric and Oceanic Technology. 34:2407-2430.   10.1175/jtech-d-17-0100.1   AbstractWebsite

This paper presents a least squares method to estimate the horizontal (isotropic or anisotropic) spatial covariance of two-dimensional orthogonal vector components, without introducing an intervening mapping step and biases, from the spatial covariance of the nonorthogonally and irregularly sampled raw scalar velocities. The field is assumed to be locally homogeneous in space and sampled in an ensemble so the unknown spatial covariance is a function of spatial lag only. The transformation between the irregular grid on which nonorthogonal scalar projections of the vector are sampled and the regular orthogonal grid on which they will be mapped is created using the geometry of the problem. The spatial covariance of the orthogonal velocity components of the field is parameterized by either the energy (power) spectrum in the wavenumber domain or the lagged covariance in the spatial domain. The energy spectrum is constrained to be nonnegative definite as part of the solution of the inverse problem. This approach is applied to three example sets of data, using nonorthogonally and irregularly sampled radial velocity data obtained from 1) a simple spectral model, 2) a regional numerical model, and 3) an array of high-frequency radars. In tests where the true covariance is known, the proposed direct approaches fitting to parameterization of the nonorthogonally and irregularly sampled raw data in the wavenumber domain and spatial domain outperform methods that map the data to a regular grid before estimating the covariance.

Cornuelle, BD, Morris MY, Roemmich DH.  1993.  An Objective Mapping Method for Estimating Geostrophic Velocity from Hydrographic Sections Including the Equator. Journal of Geophysical Research-Oceans. 98:18109-18118.   10.1029/93jc01729   AbstractWebsite

Objective mapping can remove the equatorial singularity from the problem of estimating geostrophic shear from noisy density measurements. The method uses the complete thermal wind relation, so it is valid uniformly on and off the equator. Errors in the thermal wind balance are due to neglected terms in the momentum balance, which are treated as noise in the inverse problem. The question of whether the geostrophic balance holds near the equator is restated as a need to estimate the size of the ageostrophic noise in the thermal wind equation. Objective mapping formalizes the assumptions about the magnitudes and scales of the geostrophic currents and about the magnitudes and scales of the ageostrophic terms and measurement errors. The uncertainty of the velocity estimates is calculated as part of the mapping and depends on the signal to noise ratio (geostrophic density signal to ageostrophic ''noise'') in the data, as well as the station spacing and the scales assumed for the geostrophic velocities. The method is used to map zonal velocity from a mean Hawaii-Tahiti Shuttle density section. These are compared with previous velocity estimates for the same dataset calculated using other techniques. By choosing appropriate scales, the objective map can duplicate previous results. New temperature data are presented from a repeating, high-resolution expendable bathythermograph section crossing the equator at about 170-degrees-W with four cruises a year between 1987-1991. There appear to be significant differences between this mean temperature and the shuttle mean temperature. Temperature is converted to density with the aid of a mean T-S relation and geostrophic velocity maps are calculated for the 4-year mean. The mean geostrophic undercurrent obtained from our sections is weaker than in the shuttle estimate and is centered slightly north of the equator. Enforcing symmetry about the equator removes the offset of the current, giving a stronger, but narrow undercurrent. The density field apparently includes significant (O(0.5 kg M-3)) large-scale ageostrophic variability which makes velocity estimates from single cruises poorly determined near the equator.

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.

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.

Colosi, JA, Van Uffelen LJ, Cornuelle BD, Dzieciuch MA, Worcester PF, Dushaw BD, Ramp SR.  2013.  Observations of sound-speed fluctuations in the western Philippine Sea in the spring of 2009. Journal of the Acoustical Society of America. 134:3185-3200.   10.1121/1.4818784   AbstractWebsite

As an aid to understanding long-range acoustic propagation in the Philippine Sea, statistical and phenomenological descriptions of sound-speed variations were developed. Two moorings of oceanographic sensors located in the western Philippine Sea in the spring of 2009 were used to track constant potential-density surfaces (isopycnals) and constant potential-temperature surfaces (isotherms) in the depth range 120-2000 m. The vertical displacements of these surfaces are used to estimate sound-speed fluctuations from internal waves, while temperature/salinity variability along isopycnals are used to estimate sound-speed fluctuations from intrusive structure often termed spice. Frequency spectra and vertical covariance functions are used to describe the space-time scales of the displacements and spiciness. Internal-wave contributions from diurnal and semi-diurnal internal tides and the diffuse internal-wave field [related to the Garrett-Munk (GM) spectrum] are found to dominate the sound-speed variability. Spice fluctuations are weak in comparison. The internal wave and spice frequency spectra have similar form in the upper ocean but are markedly different below 170-m depth. Diffuse internal-wave mode spectra show a form similar to the GM model, while internal-tide mode spectra scale as mode number to the minus two power. Spice decorrelates rapidly with depth, with a typical correlation scale of tens of meters.

Colosi, JA, Duda TF, Lin YT, Lynch JF, Newhall AE, Cornuelle BD.  2012.  Observations of sound-speed fluctuations on the New Jersey continental shelf in the summer of 2006. Journal of the Acoustical Society of America. 131:1733-1748.   10.1121/1.3666014   AbstractWebsite

Environmental sensors moored on the New Jersey continental shelf tracked constant density surfaces (isopycnals) for 35 days in the summer of 2006. Sound-speed fluctuations from internal-wave vertical isopycnal displacements and from temperature/salinity variability along isopycnals (spiciness) are analyzed using frequency spectra and vertical covariance functions. Three varieties of internal waves are studied: Diffuse broadband internal waves (akin to waves fitting the deep water Garrett/Munk spectrum), internal tides, and, to a lesser extent, nonlinear internal waves. These internal-wave contributions are approximately distinct in the frequency domain. It is found that in the main thermocline spicy thermohaline structure dominates the root mean square sound-speed variability, with smaller contributions coming from (in order) nonlinear internal waves, diffuse internal waves, and internal tides. The frequency spectra of internal-wave displacements and of spiciness have similar form, likely due to the advection of variable-spiciness water masses by horizontal internal-wave currents, although there are technical limitations to the observations at high frequency. In the low-frequency, internal-wave band the internal-wave spectrum follows frequency to the -1.81 power, whereas the spice spectrum shows a -1.73 power. Mode spectra estimated via covariance methods show that the diffuse internal-wave spectrum has a smaller mode bandwidth than Garrett/Munk and that the internal tide has significant energy in modes one through three. (C) 2012 Acoustical Society of America. [DOI: 10.1121/1.3666014]