Publications

Export 149 results:
Sort by: Author [ Title  (Asc)] Type Year
A B [C] D E F G H I J K L M N O P Q R S T U V W X Y Z   [Show ALL]
C
Gawarkiewicz, G, Jan S, Lermusiaux PFJ, McClean JL, Centurioni L, Taylor K, Cornuelle B, Duda TF, Wang J, Yang YJ, Sanford T, Lien RC, Lee C, Lee MA, Leslie W, Haley PJ, Niiler PP, Gopalakrishnan G, Velez-Belchi P, Lee DK, Kim YY.  2011.  Circulation and Intrusions Northeast of Taiwan: Chasing and Predicting Uncertainty in the Cold Dome. Oceanography. 24:110-121. AbstractWebsite

An important element of present oceanographic research is the assessment and quantification of uncertainty. These studies are challenging in the coastal ocean due to the wide variety of physical processes occurring on a broad range of spatial and temporal scales. In order to assess new methods for quantifying and predicting uncertainty, a joint Taiwan-US field program was undertaken in August/September 2009 to compare model forecasts of uncertainties in ocean circulation and acoustic propagation, with high-resolution in situ observations. The geographical setting was the continental shelf and slope northeast of Taiwan, where a feature called the "cold dome" frequently forms. Even though it is hypothesized that Kuroshio subsurface intrusions are the water sources for the cold dome, the dome's dynamics are highly uncertain, involving multiple scales and many interacting ocean features. During the experiment, a combination of near-surface and profiling drifters, broad-scale and high-resolution hydrography, mooring arrays, remote sensing, and regional ocean model forecasts of fields and uncertainties were used to assess mean fields and uncertainties in the region. River runoff from Typhoon Morakot, which hit Taiwan August 7-8, 2009, strongly affected shelf stratification. In addition to the river runoff, a cold cyclonic eddy advected into the region north of the Kuroshio, resulting in a cold dome formation event. Uncertainty forecasts were successfully employed to guide the hydrographic sampling plans. Measurements and forecasts also shed light on the evolution of cold dome waters, including the frequency of eddy shedding to the north-northeast, and interactions with the Kuroshio and tides. For the first time in such a complex region, comparisons between uncertainty forecasts and the model skill at measurement locations validated uncertainty forecasts. To complement the real-time model simulations, historical simulations with another model show that large Kuroshio intrusions were associated with low sea surface height anomalies east of Taiwan, suggesting that there may be some degree of predictability for Kuroshio intrusions.

Ponte, AL, Cornuelle BD.  2013.  Coastal numerical modelling of tides: Sensitivity to domain size and remotely generated internal tide. Ocean Modelling. 62:17-26.   10.1016/j.ocemod.2012.11.007   AbstractWebsite

The propagation of remotely generated superinertial internal tides constitutes a difficulty for the modelling of regional ocean tidal variability which we illustrate in several ways. First, the M2 tidal solution inside a control region located along the Southern California Bight coastline is monitored while the extent of the numerical domain is increased (up to 512 x 512 km). While the amplitude and phase of sea level averaged over the region is quasi-insensitive to domain size, a steady increase of kinetic energy, predominantly baroclinic, is observed with increasing domain size. The increasing flux of energy into the control region suggests that this trend is explained by the growing contribution from remote generation sites of internal tide which can propagate up to the control region. Increasing viscosities confirms this interpretation by lowering baroclinic energy levels and limiting their rate of increase with domain size. Doubling the grid spacing allows consideration of numerical domains 2 times larger. While the coarse grid has lower energy levels than the finer grid, the rate of energy increase with domain size appears to be slowing for the largest domain of the coarse grid simulations. Forcing the smallest domain with depth-varying tidal boundary conditions from the simulation in the largest domain produces energy levels inside the control region comparable to those in the control region for the largest domain, thereby confirming the feasibility of a nested approach. In contrast, simulations forced with a subinertial tidal constituent (K1) show that when the propagation of internal tide is limited, the control region kinetic energy is mostly barotropic and the magnitudes of variations of the kinetic energy with domain size are reduced. (C) 2012 Elsevier Ltd. All rights reserved.

Kim, SY, Cornuelle BD.  2015.  Coastal ocean climatology of temperature and salinity off the Southern California Bight: Seasonal variability, climate index correlation, and linear trend. Progress in Oceanography. 138:136-157.   10.1016/j.pocean.2015.08.001   AbstractWebsite

A coastal ocean climatology of temperature and salinity in the Southern California Bight is estimated from conductivity-temperature-depth (CTD) and bottle sample profiles collected by historical California Cooperative Oceanic Fisheries Investigation (CalCOFI) cruises (1950-2009; quarterly after 1984) off southern California and quarterly/monthly nearshore CTD surveys (within 30 km from the coast except for the surfzone; 1999-2009) off San Diego and Los Angeles. As these fields are sampled regularly in space, but not in time, conventional Fourier analysis may not be possible. The time dependent temperature and salinity fields are modeled as linear combinations of an annual cycle and its five harmonics, as well as three standard climate indices (El Nino-Southern Oscillation (ENSO), Pacific Decadal Oscillation (PDO), North Pacific Gyre Oscillation (NPGO)), the Scripps Pier temperature time series, and a mean and linear trend without time lags. Since several of the predictor indices are correlated, the indices are successively orthogonalized to eliminate ambiguity in the identification of the contributed variance of each component. Regression coefficients are displayed in both vertical transects and horizontal maps to evaluate (1) whether the temporal and spatial scales of the two data sets of nearshore and offshore observations are consistent and (2) how oceanic variability at a regional scale is related to variability in the nearshore waters. The data-derived climatology can be used to identify anomalous events and atypical behaviors in regional-scale oceanic variability and to provide background ocean estimates for mapping or modeling. (C) 2015 Elsevier Ltd. All rights reserved.

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.

Gopalakrishnan, G, Hoteit I, Cornuelle BD, Rudnick DL.  2019.  Comparison of 4DVAR and EnKF state estimates and forecasts in the Gulf of Mexico. Quarterly Journal of the Royal Meteorological Society. 145:1354-1376.   10.1002/qj.3493   AbstractWebsite

An experiment is conducted to compare four-dimensional variational (4DVAR) and ensemble Kalman filter (EnKF) assimilation systems and their predictability in the Gulf of Mexico (GoM) using the Massachusetts Institute of Technology general circulation model (MITgcm). The quality of the ocean-state estimates, forecasts, and the contribution of ensemble prediction are evaluated. The MITgcm-Estimating the Circulation and Climate of the Ocean (ECCO) 4DVAR (MITgcm-ECCO) and the MITgcm-Data Assimilation Research Testbed (DART) EnKF (MITgcm-DART) systems were used to compute two-month hindcasts (March-April, 2010) by assimilating satellite-derived along-track sea-surface height (SSH) and gridded sea-surface temperature (SST) observations. The estimates from both methods at the end of the hindcast period were then used to initialize forecasts for two months (May-June, 2010). This period was selected because a loop current (LC) eddy (Eddy Franklin: Eddy-F) detachment event occurred at the end of May 2010, immediately after the Deepwater Horizon (DwH) oil spill. Despite some differences between the setups, both systems produce analyses and forecasts of comparable quality and both solutions significantly outperformed model persistence. A reference forecast initialized from the 1/12 degrees Hybrid Coordinate Ocean Model (HYCOM)/NCODA global analysis also performed well. The EnKF experiments for sensitivity to filter parameters showed enhanced predictability when using more ensemble members and stronger covariance localization, but not for larger inflation. The EnKF experiments varying the number of assimilation cycles showed enhanced short-term (long-term) predictability with fewer (more) assimilation cycles. Additional hindcast and forecast experiments at other times of significant LC evolution showed mixed performance of both systems, which depends strongly on the background state of the GoM circulation. The present work demonstrates a practical application of both assimilation methods for the GoM and compares them in a limited number of realizations. The overall conclusion showing improved short-term (long-term) predictability for EnKF (4DVAR) carries an important caveat that the results from this study are specific to a few 4DVAR and EnKF LC eddy separation experiments in the GoM and cannot be generalized to conclude the relative performance of both methods, especially in other applications. However, some of the concepts and methods should carry over to other applications.

Worcester, PF, Cornuelle BD, Hildebrand JA, Hodgkiss WS, Duda TF, Boyd J, Howe BM, Mercer JA, Spindel RC.  1994.  A Comparison of Measured and Predicted Broad-Band Acoustic Arrival Patterns in Travel Time-Depth Coordinates at 1000-Km Range. Journal of the Acoustical Society of America. 95:3118-3128.   10.1121/1.409977   AbstractWebsite

Broadband acoustic signals were transmitted from a moored 250-Hz source to a 3-km-long vertical line array of hydrophones 1000 km distant in the eastern North Pacific Ocean during July 1989. The sound-speed field along the great circle path connecting the source and receiver was measured directly by nearly 300 expendable bathythermograph (XBT), conductivity-temperature-depth (CTD), and air-launched expendable bathythermograph (AXBT) casts while the transmissions were in progress. This experiment is unique in combining a vertical receiving array that extends over much of the water column, extensive concurrent environmental measurements, and broadband signals designed to measure acoustic travel times with 1-ms precision. The time-mean travel times of the early raylike arrivals, which are evident as wave fronts sweeping across the receiving array, and the time-mean of the times at which the acoustic reception ends (the final cutoffs) for hydrophones near the sound channel axis, are consistent with ray predictions based on the direct measurements of temperature and salinity, within measurement uncertainty. The comparisons show that subinertial oceanic variability with horizontal wavelengths shorter than 50 km, which is not resolved by the direct measurements, significantly (25 ms peak-to-peak) affects the time-mean ray travel times. The final cutoffs occur significantly later than predicted using ray theory for hydrophones more than 100-200 m off the sound channel axis. Nongeometric effects, such as diffraction at caustics, partially account for this observation.

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

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.

Mazloff, MR, Cornuelle BD, Gille ST, Verdy A.  2018.  Correlation lengths for estimating the large-scale carbon and heat content of the Southern Ocean. Journal of Geophysical Research-Oceans. 123:883-901.   10.1002/2017jc013408   AbstractWebsite

The spatial correlation scales of oceanic dissolved inorganic carbon, heat content, and carbon and heat exchanges with the atmosphere are estimated from a realistic numerical simulation of the Southern Ocean. Biases in the model are assessed by comparing the simulated sea surface height and temperature scales to those derived from optimally interpolated satellite measurements. While these products do not resolve all ocean scales, they are representative of the climate scale variability we aim to estimate. Results show that constraining the carbon and heat inventory between 35 degrees S and 70 degrees S on time-scales longer than 90 days requires approximately 100 optimally spaced measurement platforms: approximately one platform every 20 degrees longitude by 6 degrees latitude. Carbon flux has slightly longer zonal scales, and requires a coverage of approximately 30 degrees by 6 degrees. Heat flux has much longer scales, and thus a platform distribution of approximately 90 degrees by 10 degrees would be sufficient. Fluxes, however, have significant subseasonal variability. For all fields, and especially fluxes, sustained measurements in time are required to prevent aliasing of the eddy signals into the longer climate scale signals. Our results imply a minimum of 100 biogeochemical-Argo floats are required to monitor the Southern Ocean carbon and heat content and air-sea exchanges on time-scales longer than 90 days. However, an estimate of formal mapping error using the current Argo array implies that in practice even an array of 600 floats (a nominal float density of about 1 every 7 degrees longitude by 3 degrees latitude) will result in nonnegligible uncertainty in estimating climate signals.

Rudnick, DL, Gopalakrishnan G, Cornuelle BD.  2015.  Cyclonic eddies in the Gulf of Mexico: Observations by underwater gliders and simulations by numerical model. Journal of Physical Oceanography. 45:313-326.   10.1175/jpo-d-14-0138.1   AbstractWebsite

Circulation in the Gulf of Mexico (GoM) is dominated by the Loop Current (LC) and by Loop Current eddies (LCEs) that form at irregular multimonth intervals by separation from the LC. Comparatively small cyclonic eddies (CEs) are thought to have a controlling influence on the LCE, including its separation from the LC. Because the CEs are so dynamic and short-lived, lasting only a few weeks, they have proved a challenge to observe. This study addresses that challenge using underwater gliders. These gliders' data and satellite sea surface height (SSH) are used in a four-dimensional variational (4DVAR) assimilation in the Massachusetts Institute of Technology (MIT) general circulation model (MITgcm). The model serves two purposes: first, the model's estimate of ocean state allows the analysis of four-dimensional fields, and second, the model forecasts are examined to determine the value of glider data. CEs have a Rossby number of about 0.2, implying that the effects of flow curvature, cyclostrophy, to modify the geostrophic momentum balance are slight. The velocity field in CEs is nearly depth independent, while LCEs are more baroclinic, consistent with the CEs origin on the less stratified, dense side of the LCE. CEs are formed from water in the GoM, rather than the Atlantic water that distinguishes the LCE. Model forecasts are improved by glider data, using a quality metric based on satellite SSH, with the best 2-month GoM forecast rivaling the accuracy of a global hindcast.

D
Subramanian, AC, Miller AJ, Cornuelle BD, Di Lorenzo E, Weller RA, Straneo F.  2013.  A data assimilative perspective of oceanic mesoscale eddy evolution during VOCALS-REx. Atmospheric Chemistry and Physics. 13:3329-3344.   10.5194/acp-13-3329-2013   AbstractWebsite

Oceanic observations collected during the VOCALS-REx cruise time period, 1-30 November 2008, are assimilated into a regional ocean model (ROMS) using 4DVAR and then analyzed for their dynamics. Nonlinearities in the system prevent a complete 30-day fit, so two 15-day fits for 1-15 November and 16-30 November are executed using the available observations of hydrographic temperature and salinity, along with satellite fields of SST and sea-level height anomaly. The fits converge and reduce the cost function significantly, and the results indicated that ROMS is able to successfully reproduce both large-scale and smaller-scale features of the flows observed during the VOCALS-REx cruise. Particular attention is focused on an intensively studied eddy at 76 degrees W, 19 degrees S. The ROMS fits capture this eddy as an isolated rotating 3-D vortex with a strong subsurface signature in velocity, temperature and anomalously low salinity. The eddy has an average temperature anomaly of approximately -0.5 degrees C over a depth range from 50-600 m and features a cold anomaly of approximately -1 degrees C near 150 m depth. The eddy moves northwestward and elongates during the second 15-day fit. It exhibits a strong signature in the Okubo-Weiss parameter, which indicates significant non-linearity in its evolution. The heat balance for the period of the cruise from the ocean state estimate reveals that the horizontal advection and the vertical mixing processes are the dominant terms that balance the temperature tendency of the upper layer of the ocean locally in time and space. Areal averages around the eddies, for a 15-day period during the cruise, suggest that vertical mixing processes generally balance the surface heating. Although, this indicates only a small role for lateral advective processes in this region during this period, this quasi-instantaneous heat budget analysis cannot be extended to interpret the seasonal or long-term upper ocean heat budget in this region.

Dushaw, BD, Worcester PF, Munk WH, Spindel RC, Mercer JA, Howe BM, Metzger K, Birdsall TG, Andrew RK, Dzieciuch MA, Cornuelle BD, Menemenlis D.  2009.  A decade of acoustic thermometry in the North Pacific Ocean. Journal of Geophysical Research-Oceans. 114   10.1029/2008jc005124   AbstractWebsite

Over the decade 1996-2006, acoustic sources located off central California (1996 1999) and north of Kauai (1997-1999, 2002-2006) transmitted to receivers distributed throughout the northeast and north central Pacific. The acoustic travel times are inherently spatially integrating, which suppresses mesoscale variability and provides a precise measure of ray-averaged temperature. Daily average travel times at 4-day intervals provide excellent temporal resolution of the large-scale thermal field. The interannual, seasonal, and shorter-period variability is large, with substantial changes sometimes occurring in only a few weeks. Linear trends estimated over the decade are small compared to the interannual variability and inconsistent from path to path, with some acoustic paths warming slightly and others cooling slightly. The measured travel times are compared with travel times derived from four independent estimates of the North Pacific: (1) climatology, as represented by the World Ocean Atlas 2005 (WOA05); (2) objective analysis of the upper-ocean temperature field derived from satellite altimetry and in situ profiles; (3) an analysis provided by the Estimating the Circulation and Climate of the Ocean project, as implemented at the Jet Propulsion Laboratory (JPL-ECCO); and (4) simulation results from a high-resolution configuration of the Parallel Ocean Program (POP) model. The acoustic data show that WOA05 is a better estimate of the time mean hydrography than either the JPL-ECCO or the POP estimates, both of which proved incapable of reproducing the observed acoustic arrival patterns. The comparisons of time series provide a stringent test of the large-scale temperature variability in the models. The differences are sometimes substantial, indicating that acoustic thermometry data can provide significant additional constraints for numerical ocean models.

Worcester, P, Dushaw BD, Andrew RK, Howe BM, Mercer JA, Spindel RC, Cornuelle B, Dzieciuch M, Birdsall TG, Metzger K, Menemenlis D.  2008.  A decade of acoustic thermometry in the North Pacific Ocean: Using long-range acoustic travel times to test gyre-scale temperature variability derived from other observations and ocean models. Journal of the Acoustical Society of America. 123 AbstractWebsite

Large-scale, range- and depth-averaged temperatures in the North Pacific Ocean were measured by long-range acoustic transmissions over the decade 1996-2006. Acoustic sources off central California and north of Kauai transmitted to receivers throughout the North Pacific. Even though acoustic travel times are spatially integrating, suppressing mesoscale variability and providing a precise measure of large-scale temperature, the travel times sometimes vary significantly on time scales of only a few weeks. The interannual variability is large, with no consistent warming or cooling trends. Comparison of the measured travel times with travel times derived from (i) the World Ocean Atlas 2005 (WOA05), (ii) an upper ocean temperature estimate derived from satellite altimetry and in situ profiles, (iii) an analysis provided by the Estimating the Circulation and Climate of the Ocean (ECCO) project, and (iv) simulation results from a high-resolution configuration of the Parallel Ocean Program (POP) show similarities, but also reveal substantial differences. The differences suggest that the data can provide significant additional constraints for numerical ocean simulations. The acoustic data show that WOA05 is a much better estimate of the time-mean hydrography than either the ECCO or POP estimates and provide significantly better time resolution for large-scale ocean variability than can be derived from satellite altimetry and in situ profiles.

Kim, SY, Cornuelle BD, Terrill EJ.  2010.  Decomposing observations of high-frequency radar-derived surface currents by their forcing mechanisms: Decomposition techniques and spatial structures of decomposed surface currents. Journal of Geophysical Research-Oceans. 115   10.1029/2010jc006222   AbstractWebsite

Surface current observations from a high-frequency radar network deployed in southern San Diego are decomposed according to their driving forces: pure tides and their neighboring off-band energy, local winds, and low frequency. Several superposed ocean responses are present as a result of the complicated bottom topography and relatively weak winds off southern San Diego, as opposed to coastal regions where circulation can be explained by a dominant forcing mechanism. This necessitates an application of a statistical decomposition approach. Surface currents coherent with pure tides are calculated using harmonic analysis. Locally wind-driven surface currents are estimated by regression of observed winds on observed surface currents. The dewinded and detided surface currents are filtered by weighted least-squares fitting assuming white noise and three colored signal bands: low-frequency band (less than 0.4 cycles per day) and near-tidal peaks at the diurnal (K-1) and semidiurnal (M-2) frequencies. The spatial and temporal variability of each part of the decomposed surface currents is investigated in terms of ocean response to the driving forces. In addition, the spatial correlations of individual components exhibit Gaussian and exponential shapes with varying decorrelation length scales.

Kim, SY, Cornuelle BD, Terrill EJ.  2010.  Decomposing observations of high-frequency radar-derived surface currents by their forcing mechanisms: Locally wind-driven surface currents. Journal of Geophysical Research-Oceans. 115   10.1029/2010jc006223   AbstractWebsite

The wind impulse response function and transfer function for high-frequency radar-derived surface currents off southern San Diego are calculated using several local wind observations. The spatial map of the transfer function reflects the influence of the coast on wind-current dynamics. Near the coast (within 20 km from the shoreline), the amplitudes of the transfer function at inertial and diurnal frequencies are reduced due to effects of coastline and bottom bathymetry. Meanwhile, the amplitude of low-frequency currents increases near the coast, which is attributed to the local geostrophic balance between cross-shore pressure gradients against the coast and currents. Locally wind-driven surface currents are estimated from the data-derived response function, and their power spectrum shows a strong diurnal peak superposed on a red spectrum, similar to the spectra of observed winds. Current magnitudes and veering angles to a quasi-steady wind are typically 2-5% of the wind speed and vary 50 degrees-90 degrees to the right of the wind, respectively. A wind skill map is introduced to present the fractional variance of surface currents explained by local winds as a verification tool for wind data quality and relevance. Moreover, the transfer functions in summer and winter are presented to examine the seasonal variation in ocean surface current response to the wind associated with stratification change.

Behringer, D, Birdsall T, Brown M, Cornuelle B, Heinmiller R, Knox R, Metzger K, Munk W, Spiesberger J, Spindel R, Webb D, Worcester P, Wunsch C.  1982.  A demonstration of ocean acoustic tomography. Nature. 299:121-125.   10.1038/299121a0   AbstractWebsite

Over the past decade oceanographers have become increasingly aware of an intense and compact ocean ‘mesoscale’ eddy structure (the ocean weather) that is superimposed on a generally sluggish large-scale circulation (the ocean climate). Traditional ship-based observing systems are not adequate for monitoring the ocean at mesoscale resolution. A 1981 experiment mapped the waters within a 300 × 300 km square south-west of Bermuda, using a peripheral array of moored midwater acoustic sources and receivers. The variable acoustic travel times between all source–receiver pairs were used to construct the three-dimensional (time-variable) eddy fields, using inverse theory. Preliminary results from inversions are consistent with the shipborne and airborne surveys.

Verlinden, CMA, Sarkar J, Cornuelle BD, Kuperman WA.  2017.  Determination of acoustic waveguide invariant using ships as sources of opportunity in a shallow water marine environment. Journal of the Acoustical Society of America. 141:EL102-EL107.   10.1121/1.4976112   AbstractWebsite

The waveguide invariant (WGI) is a property that can be used to localize acoustic radiators and extract information about the environment. Here the WGI is determined using ships as sources of opportunity, tracked using the Automatic Identification System (AIS). The relationship between range, acoustic intensity, and frequency for a ship in a known position is used to determine the WGI parameter beta. These b values are interpolated and a map of b is generated. The method is demonstrated using data collected in a field experiment on a single hydrophone in a shallow water environment off the coast of Southern California. (C) 2017 Acoustical Society of America

Roemmich, D, Cornuelle B.  1987.  Digitization and Calibration of the Expendable Bathythermograph. Deep-Sea Research Part a-Oceanographic Research Papers. 34:299-307.   10.1016/0198-0149(87)90088-4   AbstractWebsite

A study was undertaken of signal digitization and temperature calibration in expendable bathythermographs (XBT's) to learn how to minimize temperature errors in that system. An XBT digitizer was built into a PC-type microcomputer and used to calibrate 24 XBT probes at 5 temperature points, and later, to calibrate 72 probes at a single temperature. Twenty of the first set of probes were fastened rigidly in pairs and dropped in the ocean as a field test of the calibrations. Calibration of individual probes reduced the standard deviation of temperature calibration errors from around 0.05°C to <0.01°C. The calibration procedure is simple and nondestructive, so the probes can be used normally after calibration. Errors in the temperature digitizer can be held to about 0.01°C by periodic adjustment. An advantage of the PC-based digitizer is the ease with which the calibrations are accomplished and applied to the ocean tracers. Two substantial sources of transient systematic error in XBT temperatures are mentioned: an electronic transient lasting about 0.1 s occurs on entry of the probe into seawater, and a longer transient is due to the thermal mass of the XBT nose.

Ubelmann, C, Cornuelle B, Fu LL.  2016.  Dynamic mapping of along-track ocean altimetry: Method and performance from observing system simulation experiments. Journal of Atmospheric and Oceanic Technology. 33:1691-1699.   10.1175/jtech-d-15-0163.1   AbstractWebsite

Simulated along-track ocean altimetry data were used to implement the use of a nonlinear dynamic propagator to perform three-dimensional (time and 2D space) interpolation of mesoscale sea surface height (SSH). The method is an inverse approach to processing altimetry data unevenly sampled in time and space into high-level gridded altimetry maps. The inverse approach, similar to the standard objective mapping, contains some correction terms to the innovation vectors to account for nonlinear dynamics. Another key improvement is to solve for the covariance functions through a Green's function approach. From the Observing System Simulation Experiments carried out to simulate a three-satellite constellation over the Gulf Stream region, the new method can significantly reduce mapping errors and improve the resolving capabilities compared to the standard linear objective analysis such as that used by the AVISO gridding.

E
Hoteit, I, Cornuelle B, Heimbach P.  2010.  An eddy-permitting, dynamically consistent adjoint-based assimilation system for the tropical Pacific: Hindcast experiments in 2000. Journal of Geophysical Research-Oceans. 115   10.1029/2009jc005437   AbstractWebsite

An eddy-permitting adjoint-based assimilation system has been implemented to estimate the state of the tropical Pacific Ocean. The system uses the Massachusetts Institute of Technology's general circulation model and its adjoint. The adjoint method is used to adjust the model to observations by controlling the initial temperature and salinity; temperature, salinity, and horizontal velocities at the open boundaries; and surface fluxes of momentum, heat, and freshwater. The model is constrained with most of the available data sets in the tropical Pacific, including Tropical Atmosphere and Ocean, ARGO, expendable bathythermograph, and satellite SST and sea surface height data, and climatologies. Results of hindcast experiments in 2000 suggest that the iterated adjoint-based descent is able to significantly improve the model consistency with the multivariate data sets, providing a dynamically consistent realization of the tropical Pacific circulation that generally matches the observations to within specified errors. The estimated model state is evaluated both by comparisons with observations and by checking the controls, the momentum balances, and the representation of small-scale features that were not well sampled by the observations used in the assimilation. As part of these checks, the estimated controls are smoothed and applied in independent model runs to check that small changes in the controls do not greatly change the model hindcast. This is a simple ensemble-based uncertainty analysis. In addition, the original and smoothed controls are applied to a version of the model with doubled horizontal resolution resulting in a broadly similar "downscaled'' hindcast, showing that the adjustments are not tuned to a single configuration (meaning resolution, topography, and parameter settings). The time-evolving model state and the adjusted controls should be useful for analysis or to supply the forcing, initial, and boundary conditions for runs of other models.

Rasmussen, LL, Cornuelle BD, Levin LA, Largier JL, Di Lorenzo E.  2009.  Effects of small-scale features and local wind forcing on tracer dispersion and estimates of population connectivity in a regional scale circulation model. Journal of Geophysical Research-Oceans. 114   10.1029/2008jc004777   AbstractWebsite

A small-scale model of the Southern California-Northern Baja California coastline has been developed to explore dispersion over the continental shelf, with specific attention to physical parameters pertinent to simulations of larval dispersal and population connectivity. The ROMS simulation employs a nested grid system, with an inner domain resolution of 600 m and an outer domain resolution of 1.5 km. Realistic bathymetry and forcing were employed to investigate the effects of passive transport of tracers introduced at locations with known communities of mytilid mussels along the coastline. The effects of topographic resolution, boundary conditions, and choice of meteorological forcing products on dispersion rates, tracer trajectories, and the subsequent measures of population connectivity were examined. In particular, the choice of wind forcing product resulted in different circulation patterns and tracer trajectories and had especially important consequences on measures of larval connectivity such as self-seeding, potential for larval settlement ( import), and contribution to the pool of available larvae ( export). While some forcing products performed better when model data were compared to field measurements, no product was clearly superior. The uncertainty in results, which may appear minor in larger-scale temperature or surface velocity fields, is significant when examining a sensitive passive tracer. This modeling uncertainty needs to be addressed when interpreting connectivity results.

Gemba, KL, Sarkar J, Cornuelle B, Hodgkiss WS, Kuperman WA.  2018.  Estimating relative channel impulse responses from ships of opportunity in a shallow water environment. The Journal of the Acoustical Society of America. 144:1231-1244.   10.1121/1.5052259   Abstract

The uncertainty of estimating relative channel impulse responses (CIRs) obtained using the radiated signature from a ship of opportunity is investigated. The ship observations were taken during a 1.4 km (11 min) transect in a shallow water environment during the Noise Correlation 2009 (NC09) experiment. Beamforming on the angle associated with the direct ray-path yields an estimate of the ship signature, subsequently used in a matched filter. Relative CIRs are estimated every 2.5 s independently at three vertical line arrays (VLAs). The relative arrival-time uncertainty is inversely proportional to source bandwidth and CIR signal-to-noise ratio, and reached a minimum standard deviation of 5 μs (equivalent to approximately 1 cm spatial displacement). Time-series of direct-path relative arrival-times are constructed for each VLA element across the 11 min observation interval. The overall structure of these time-series compares favorably with that predicted from an array element localization model. The short-term standard deviations calculated on the direct-path (7 μs) and bottom-reflected-path (17 μs) time-series are in agreement with the predicted arrival-time accuracies. The implications of these observed arrival-time accuracies in the context of estimating sound speed perturbations and bottom-depth are discussed.

Verdy, A, Cornuelle B, Mazloff MR, Rudnick DL.  2017.  Estimation of the tropical Pacific Ocean state 2010-13. Journal of Atmospheric and Oceanic Technology. 34:1501-1517.   10.1175/jtech-d-16-0223.1   AbstractWebsite

A data-assimilating 1/38 regional dynamical ocean model is evaluated on its ability to synthesize components of the Tropical Pacific Ocean Observing System. The four-dimensional variational data assimilation (4DVAR) method adjusts initial conditions and atmospheric forcing for overlapping 4-month model runs, or hindcasts, that are then combined to give an ocean state estimate for the period 2010-13. Consistency within uncertainty with satellite SSH and Argo profiles is achieved. Comparison to independent observations from Tropical Atmosphere Ocean (TAO) moorings shows that for time scales shorter than 100 days, the state estimate improves estimates of TAO temperature relative to an optimally interpolated Argo product. The improvement is greater at time scales shorter than 20 days, although unpredicted variability in the TAO temperatures implies that TAO observations provide significant information in that band. Larger discrepancies between the state estimate and independent observations from Spray gliders deployed near the Galapagos, Palau, and Solomon Islands are attributed to insufficient model resolution to capture the dynamics in strong current regions and near coasts. The sea surface height forecast skill of the model is assessed. Model forecasts using climatological forcing and boundary conditions are more skillful than climatology out to 50 days compared to persistence, which is a more skillful forecast than climatology out to approximately 20 days. Hindcasts using reanalysis products for atmospheric forcing and open boundary conditions are more skillful than climatology for approximately 120 days or longer, with the exact time scale depending on the accuracy of the state estimate used for initializing and on the reanalysis forcing. Estimating the model representational error is a goal of these experiments.

Fujii, Y, Cummings J, Xue Y, Schiller A, Lee T, Balmaseda MA, Remy E, Masuda S, Brassington G, Alves O, Cornuelle B, Martin M, Oke P, Smith G, Yang XS.  2015.  Evaluation of the Tropical Pacific Observing System from the ocean data assimilation perspective. Quarterly Journal of the Royal Meteorological Society. 141:2481-2496.   10.1002/qj.2579   AbstractWebsite

The drastic reduction in the number of observation data from the Tropical Atmospheric Ocean (TAO)/Triangle Trans-Ocean Buoy Network (TRITON) array since 2012 has given rise to a need to assess the impact of those data in ocean data assimilation (DA) systems. This article provides a review of existing studies evaluating the impacts of data from the TAO/TRITON array and other components of the Tropical Pacific Observing System (TPOS) on current ocean DA systems used for a variety of operational and research applications. It can be considered as background information that can guide the evaluation exercise of TPOS. Temperature data from TAO/TRITON array are assimilated in most ocean DA systems which cover the tropical Pacific in order to constrain the ocean heat content, stratification, and circulation. It is shown that the impacts of observation data depend considerably on the system and application. The presence of model error often makes the results difficult to interpret. Nevertheless there is consensus that the data from TAO/TRITON generally have positive impacts complementary to Argo floats. In the equatorial Pacific, the impacts are generally around the same level or larger than those of Argo. We therefore conclude that, with the current configuration of TPOS, the loss of the TAO/TRITON data is having a significant detrimental impact on many applications based on ocean DA systems. This conclusion needs to be kept under review because the equatorial coverage by Argo is expected to improve in the future.

Worcester, PF, Lynch JF, Morawitz WML, Pawlowicz R, Sutton PJ, Cornuelle BD, Johannessen OM, Munk WH, Owens WB, Shuchman R, Spindel RC.  1993.  Evolution of the Large-Scale Temperature-Field in the Greenland Sea During 1988-89 from Tomographic Measurements. Geophysical Research Letters. 20:2211-2214.   10.1029/93gl02373   AbstractWebsite

The Greenland Sea Ocean Acoustic Tomography Experiment was conducted during 1988-89, as one component of the international Greenland Sea Project, to study deep water formation and the response of the gyre to variations in wind stress and ice cover. Six acoustic transceivers moored in an array 200-km across transmitted to one another at four hour intervals. Near the end of February, 1989, a sub-surface temperature maximum at several hundred meters depth disappeared over a suprisingly large area of the central Greenland Sea . While the water column was modified to about 1000 m depth over much of the gyre, the surface remained colder than the deeper water, contrary to what might be expected from simple models of convective renewal.