Publications

Export 43 results:
Sort by: Author Title Type [ Year  (Desc)]
2017
Sagen, H, Worcester PF, Dzieciuch MA, Geyer F, Sandven S, Babiker M, Beszczynska-Moller A, Dushaw BD, Cornuelle B.  2017.  Resolution, identification, and stability of broadband acoustic arrivals in Fram Strait. Journal of the Acoustical Society of America. 141:2055-2068.   10.1121/1.4978780   AbstractWebsite

An ocean acoustic tomography system consisting of three moorings with low frequency, broad-band transceivers and a moored receiver located approximately in the center of the triangle formed by the transceivers was installed in the central, deep-water part of Fram Strait during 2010-2012. Comparisons of the acoustic receptions with predictions based on hydrographic sections show that the oceanographic conditions in Fram Strait result in complex arrival patterns in which it is difficult to resolve and identify individual arrivals. In addition, the early arrivals are unstable, with the arrival structures changing significantly over time. The stability parameter a suggests that the instability is likely not due to small-scale variability, but rather points toward strong mesoscale variability in the presence of a relatively weak sound channel as being largely responsible. The estimator-correlator [Dzieciuch, J. Acoust. Soc. Am. 136, 2512-2522 (2014)] is shown to provide an objective formalism for generating travel-time series given the complex propagation conditions. Because travel times obtained from the estimator-correlator are not associated with resolved, identified ray arrivals, inverse methods are needed that do not use sampling kernels constructed from geometric ray paths. One possible approach would be to use travel-time sensitivity kernels constructed for the estimator-correlator outputs. (C) 2017 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license.

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

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.

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.

Dzieciuch, MA, Cornuelle BD, Skarsoulis EK.  2013.  Structure and stability of wave-theoretic kernels in the ocean. Journal of the Acoustical Society of America. 134:3318-3331.   10.1121/1.4818846   AbstractWebsite

Wave-theoretic modeling can be applied to obtain travel-time sensitivity kernels (TSKs) representing the amount ray travel times are affected by sound-speed variations anywhere in the medium. This work explores the spatial frequency content of the TSK compared to expected ocean variability. It also examines the stability of the TSK in environments that produce strong sensitivity of ray paths to initial conditions. The conclusion is that the linear TSK model is an effective predictor of travel-time changes and that the rays perform nearly as well as the full-wave kernel. The TSK is examined in physical space and in wavenumber space, and it is found that this is the key to understanding how the travel time reacts to ocean perturbations. There are minimum vertical and horizontal length scales of ocean perturbations that are required for the travel time to be affected. The result is that the correspondence between true travel times and those calculated from the kernel is high for large-scale perturbations and somewhat less for the small scales. This demonstrates the validity of ray-based inversion of travel time observations for the cases under study. (C) 2013 Acoustical Society of America.

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

Todd, RE, Rudnick DL, Mazloff MR, Cornuelle BD, Davis RE.  2012.  Thermohaline structure in the California Current System: Observations and modeling of spice variance. Journal of Geophysical Research-Oceans. 117   10.1029/2011jc007589   AbstractWebsite

Upper ocean thermohaline structure in the California Current System is investigated using sustained observations from autonomous underwater gliders and a numerical state estimate. Both observations and the state estimate show layers distinguished by the temperature and salinity variability along isopycnals (i.e., spice variance). Mesoscale and submesoscale spice variance is largest in the remnant mixed layer, decreases to a minimum below the pycnocline near 26.3 kg m(-3), and then increases again near 26.6 kg m(-3). Layers of high (low) meso-and submesoscale spice variance are found on isopycnals where large-scale spice gradients are large (small), consistent with stirring of large-scale gradients to produce smaller scale thermohaline structure. Passive tracer adjoint calculations in the state estimate are used to investigate possible mechanisms for the formation of the layers of spice variance. Layers of high spice variance are found to have distinct origins and to be associated with named water masses; high spice variance water in the remnant mixed layer has northerly origin and is identified as Pacific Subarctic water, while the water in the deeper high spice variance layer has southerly origin and is identified as Equatorial Pacific water. The layer of low spice variance near 26.3 kg m(-3) lies between the named water masses and does not have a clear origin. Both effective horizontal diffusivity, kappa(h), and effective diapycnal diffusivity, kappa(v), are elevated relative to the diffusion coefficients set in the numerical simulation, but changes in kappa(h) and kappa(v) with depth are not sufficient to explain the observed layering of thermohaline structure.

2011
Song, H, Miller AJ, Cornuelle BD, Di Lorenzo E.  2011.  Changes in upwelling and its water sources in the California Current System driven by different wind forcing. Dynamics of Atmospheres and Oceans. 52:170-191.   10.1016/j.dynatmoce.2011.03.001   AbstractWebsite

In the California Current System (CCS), upwelling is one of the most important features that enrich the coastal ecosystem. It is highly dependent on both wind stress and wind stress curl, because they contribute to the upwelling system through Ekman transport away from the coast and Ekman pumping as a result of the surface divergence, respectively. Various wind stress products are known to contain sharply different patterns of wind stress, and well-resolved wind forcing products have been shown to drive stronger upwelling due to their better-resolved wind stress curl in previous studies. However, sensitivities of upwelling to changes in wind stress patterns, and each of their control to the source waters and paths of the upwelling cells, are not yet well known for the CCS. Here we study these effects using the Regional Ocean Modeling System (ROMS) and its adjoint model under idealized wind stress forcing patterns representing three widely-used products in addition to a constant wind stress field (no curl): the NCEP/NCAR Reanalysis, the QuikSCAT satellite observations, and the Regional Spectral Model (RSM) downscaling. Changes in currents and isopycnal patterns during the upwelling season are first studied in ROMS under the four different wind stress fields. The model simulations show that the locations of the core of the equatorward flow and the gradient of the cross-shore isopycnals are controlled by the wind stress curl field. The core of the equatorward flow is found under negative wind stress curl, and a deeper upwelling cell is found as the gradient from positive and negative wind stress curl increases. Source waters for the upwelling in each of the four wind stress patterns are investigated using the ROMS adjoint model. The simulations follow a passive tracer backward in time and track the source waters for upwelling in two key areas of interest: inshore and offshore of the Point Sur region of California. The upwelling source waters depend strongly on the depth of the upwelling cell and the alongshore current location. We further relate these results to recent studies of the observed trends in upwelling favorable winds and consequent wind stress curl changes in the CCS. (C) 2011 Elsevier B.V. All rights reserved.

Skarsoulis, EK, Cornuelle BD, Dzieciuch MA.  2011.  Second-Order Sensitivity of Acoustic Travel Times to Sound Speed Perturbations. Acta Acustica United with Acustica. 97:533-543.   10.3813/aaa.918434   AbstractWebsite

The second-order sensitivity of finite-frequency acoustic travel times to sound speed perturbations in range-independent environments is studied. Using the notion of peak arrivals and the normal-mode representation of the Green's function first- and second-order perturbation expressions are derived for the travel times in terms of the underlying perturbations in the Green's function and finally in the sound speed profile. The resulting theoretical expressions are numerically validated. Assuming small and local perturbations the non-linear effects appear to be strongest for sound speed perturbations taking place close to the lower turning depths of the corresponding eigenrays. At the upper turning depths - in the case of temperate propagation conditions - the effects are much weaker due to the larger sound speed gradients. The magnitude of the second-order sensitivity of travel times relative to the first-order sensitivity can be used to obtain an estimate for the limits of linearity.

Todd, RE, Rudnick DL, Mazloff MR, Davis RE, Cornuelle BD.  2011.  Poleward flows in the southern California Current System: Glider observations and numerical simulation. Journal of Geophysical Research-Oceans. 116   10.1029/2010jc006536   AbstractWebsite

Three years of continuous Spray glider observations in the southern California Current System (CCS) are combined with a numerical simulation to describe the mean and variability of poleward flows in the southern CCS. Gliders provide upper ocean observations with good across-shore and temporal resolution along two across-shore survey lines while the numerical simulation provides a dynamically consistent estimate of the ocean state. Persistent poleward flows are observed in three areas: within 100 km of the coast at Point Conception, within the Southern California Bight (SCB), and offshore of the SCB and the Santa Rosa Ridge (SRR). Poleward transport by the flows within the SCB and offshore of the SRR exceeds the poleward transport off Point Conception, suggesting that the poleward flows are not continuous over the 225 km between observation lines. The numerical simulation shows offshore transport between the survey lines that is consistent with some of the poleward flow turning offshore before reaching Point Conception. The poleward current offshore of the SRR is unique in that it is strongest at depths greater than 350 m and it is observed to migrate westward away from the coast. This westward propagation is tied to westward propagating density anomalies originating in the SCB during the spring-summer upwelling season when wind stress curl is most strongly positive. The across-shore wave number, frequency, and phase speed of the westward propagation and the lack of across-shore transport of salinity along isopycnals are consistent with first-mode baroclinic Rossby dynamics.

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.

2009
Skarsoulis, EK, Cornuelle BD, Dzieciuch MA.  2009.  Travel-time sensitivity kernels in long-range propagation. Journal of the Acoustical Society of America. 126:2223-2233.   10.1121/1.3224835   AbstractWebsite

Wave-theoretic travel-time sensitivity kernels (TSKs) are calculated in two-dimensional (2D) and three-dimensional (3D) environments and their behavior with increasing propagation range is studied and compared to that of ray-theoretic TSKs and corresponding Fresnel-volumes. The differences between the 2D and 3D TSKs average out when horizontal or cross-range marginals are considered, which indicates that they are not important in the case of range-independent sound-speed perturbations or perturbations of large scale compared to the lateral TSK extent. With increasing range, the wave-theoretic TSKs expand in the horizontal cross-range direction, their cross-range extent being comparable to that of the corresponding free-space Fresnel zone, whereas they remain bounded in the vertical. Vertical travel-time sensitivity kernels (VTSKs)-one-dimensional kernels describing the effect of horizontally uniform sound-speed changes on travel-times-are calculated analytically using a perturbation approach, and also numerically, as horizontal marginals of the corresponding TSKs. Good agreement between analytical and numerical VTSKs, as well as between 2D and 3D VTSKs, is found. As an alternative method to obtain wave-theoretic sensitivity kernels, the parabolic approximation is used; the resulting TSKs and VTSKs are in good agreement with normal-mode results. With increasing range, the wave-theoretic VTSKs approach the corresponding ray-theoretic sensitivity kernels. (C) 2009 Acoustical Society of America. [DOI: 10.1121/1.3224835]

Moore, AM, Arango HG, Di Lorenzo E, Miller AJ, Cornuelle BD.  2009.  An Adjoint Sensitivity Analysis of the Southern California Current Circulation and Ecosystem. Journal of Physical Oceanography. 39:702-720.   10.1175/2008jpo3740.1   AbstractWebsite

Adjoint methods of sensitivity analysis were applied to the California Current using the Regional Ocean Modeling Systems (ROMS) with medium resolution, aimed at diagnosing the circulation sensitivity to variations in surface forcing. The sensitivities of coastal variations in SST, eddy kinetic energy, and baroclinic instability of complex time-evolving flows were quantified. Each aspect of the circulation exhibits significant interannual and seasonal variations in sensitivity controlled by mesoscale circulation features. Central California SST is equally sensitive to wind stress and surface heat flux, but less so to wind stress curl, displaying the greatest sensitivity when upwelling-favorable winds are relaxing and the least sensitivity during the peak of upwelling. SST sensitivity is typically 2-4 times larger during summer than during spring, although larger variations occur during some years. The sensitivity of central coast eddy kinetic energy to surface forcing is constant on average throughout the year. Perturbations in the wind that align with mesoscale eddies to enhance the strength of the circulation by local Ekman pumping yield the greatest sensitivities. The sensitivity of the potential for baroclinic instability is greatest when nearshore horizontal temperature gradients are largest, and it is associated with variations in wind stress concentrated along the core of the California Current. The sensitivity varies by a factor of similar to 1.5 throughout the year. A new and important aspect of this work is identification of the complex flow dependence and seasonal dependence of the sensitivity of the ROMS California Current System (CCS) circulation to variations in surface forcing that was hitherto not previously appreciated.

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.

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.

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

Haidvogel, DB, Arango H, Budgell WP, Cornuelle BD, Curchitser E, Di Lorenzo E, Fennel K, Geyer WR, Hermann AJ, Lanerolle L, Levin J, McWilliams JC, Miller AJ, Moore AM, Powell TM, Shchepetkin AF, Sherwood CR, Signell RP, Warner JC, Wilkin J.  2008.  Ocean forecasting in terrain-following coordinates: Formulation and skill assessment of the Regional Ocean Modeling System. Journal of Computational Physics. 227:3595-3624.   10.1016/j.jcp.2007.06.016   AbstractWebsite

Systematic improvements in algorithmic design of regional ocean circulation models have led to significant enhancement in simulation ability across a wide range of space/time scales and marine system types. As an example, we briefly review the Regional Ocean Modeling System, a member of a general class of three-dimensional, free-surface, terrain-following numerical models. Noteworthy characteristics of the ROMS computational kernel include: consistent temporal averaging of the barotropic mode to guarantee both exact conservation and constancy preservation properties for tracers; redefined barotropic pressure-gradient terms to account for local variations in the density field; vertical interpolation performed using conservative parabolic splines; and higher-order, quasi-monotone advection algorithms. Examples of quantitative skill assessment are shown for a tidally driven estuary, an ice-covered high-latitude sea, a wind- and buoyancy-forced continental shelf, and a mid-latitude ocean basin. The combination of moderate-order spatial approximations, enhanced conservation properties, and quasi-monotone advection produces both more robust and accurate, and less diffusive, solutions than those produced in earlier terrain-following ocean models. Together with advanced methods of data assimilation and novel observing system technologies, these capabilities constitute the necessary ingredients for multi-purpose regional ocean prediction systems. (c) 2007 Elsevier Inc. All rights reserved.

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

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.

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

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

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

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

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

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

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
2000
Miller, AJ, Di Lorenzo E, Neilson DJ, Cornuelle BD, Moisan JR.  2000.  Modeling CalCOFI observations during El Nino: Fitting physics and biology. California Cooperative Oceanic Fisheries Investigations Reports. 41:87-97. AbstractWebsite

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

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

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