Publications

Export 58 results:
Sort by: Author Title Type [ Year  (Desc)]
2019
Subramanian, AC, Balmaseda MA, Centurioni L, Chattopadhyay R, Cornuelle BD, DeMott C, Flatau M, Fujii Y, Giglio D, Gille ST, Hamill TM, Hendon H, Hoteit I, Kumar A, Lee JH, Lucas AJ, Mahadevan A, Matsueda M, Nam S, Paturi S, Penny SG, Rydbeck A, Sun R, Takaya Y, Tandon A, Todd RE, Vitart F, Yuan DL, Zhang CD.  2019.  Ocean observations to improve our understanding, modeling, and forecasting of subseasonal-to-seasonal variability. Frontiers in Marine Science. 6   10.3389/fmars.2019.00427   AbstractWebsite

Subseasonal-to-seasonal (S2S) forecasts have the potential to provide advance information about weather and climate events. The high heat capacity of water means that the subsurface ocean stores and re-releases heat (and other properties) and is an important source of information for S2S forecasts. However, the subsurface ocean is challenging to observe, because it cannot be measured by satellite. Subsurface ocean observing systems relevant for understanding, modeling, and forecasting on S2S timescales will continue to evolve with the improvement in technological capabilities. The community must focus on designing and implementing low-cost, high-value surface and subsurface ocean observations, and developing forecasting system capable of extracting their observation potential in forecast applications. S2S forecasts will benefit significantly from higher spatio-temporal resolution data in regions that are sources of predictability on these timescales (coastal, tropical, and polar regions). While ENSO has been a driving force for the design of the current observing system, the subseasonal time scales present new observational requirements. Advanced observation technologies such as autonomous surface and subsurface profiling devices as well as satellites that observe the ocean-atmosphere interface simultaneously can lead to breakthroughs in coupled data assimilation (CDA) and coupled initialization for S2S forecasts.

Boas, ABV, Ardhuin F, Ayet A, Bourassa MA, Brandt P, Chapron B, Cornuelle BD, Farrar JT, Fewings MR, Fox-Kemper B, Gille ST, Gommenginger C, Heimbach P, Hell MC, Li Q, Mazloff MR, Merrifield ST, Mouche A, Rio MH, Rodriguez E, Shutler JD, Subramanian AC, Terrill EJ, Tsamados M, Ubelmann C, van Sebille E.  2019.  Integrated observations of global surface winds, currents, and waves: Requirements and challenges for the next decade. Frontiers in Marine Science. 6   10.3389/fmars.2019.00425   AbstractWebsite

Ocean surface winds, currents, and waves play a crucial role in exchanges of momentum, energy, heat, freshwater, gases, and other tracers between the ocean, atmosphere, and ice. Despite surface waves being strongly coupled to the upper ocean circulation and the overlying atmosphere, efforts to improve ocean, atmospheric, and wave observations and models have evolved somewhat independently. From an observational point of view, community efforts to bridge this gap have led to proposals for satellite Doppler oceanography mission concepts, which could provide unprecedented measurements of absolute surface velocity and directional wave spectrum at global scales. This paper reviews the present state of observations of surface winds, currents, and waves, and it outlines observational gaps that limit our current understanding of coupled processes that happen at the air-sea-ice interface. A significant challenge for the coming decade of wind, current, and wave observations will come in combining and interpreting measurements from (a) wave-buoys and high-frequency radars in coastal regions, (b) surface drifters and wave-enabled drifters in the open-ocean, marginal ice zones, and wave-current interaction "hot-spots," and (c) simultaneous measurements of absolute surface currents, ocean surface wind vector, and directional wave spectrum from Doppler satellite sensors.

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.

Heimbach, P, Fukumori I, Hills CN, Ponte RM, Stammer D, Wunsch C, Campin JM, Cornuelle B, Fenty I, Forget G, Kohl A, Mazloff M, Menemenlis D, Nguyen AT, Piecuch C, Trossman D, Verdy A, Wang O, Zhang H.  2019.  Putting it all together: Adding value to the global ocean and climate observing systems with complete self-consistent ocean state and parameter estimates. Frontiers in Marine Science. 6   10.3389/fmars.2019.00055   AbstractWebsite

In 1999, the consortium on Estimating the Circulation and Climate of the Ocean (ECCO) set out to synthesize the hydrographic data collected by the World Ocean Circulation Experiment (WOCE) and the satellite sea surface height measurements into a complete and coherent description of the ocean, afforded by an ocean general circulation model. Twenty years later, the versatility of ECCO's estimation framework enables the production of global and regional ocean and sea-ice state estimates, that incorporate not only the initial suite of data and its successors, but nearly all data streams available today. New observations include measurements from Argo floats, marine mammal-based hydrography, satellite retrievals of ocean bottom pressure and sea surface salinity, as well as ice-tethered profiled data in polar regions. The framework also produces improved estimates of uncertain inputs, including initial conditions, surface atmospheric state variables, and mixing parameters. The freely available state estimates and related efforts are property-conserving, allowing closed budget calculations that are a requisite to detect, quantify, and understand the evolution of climate-relevant signals, as mandated by the Coupled Model Intercomparison Project Phase 6 (CMIP6) protocol. The solutions can be reproduced by users through provision of the underlying modeling and assimilation machinery. Regional efforts have spun off that offer increased spatial resolution to better resolve relevant processes. Emerging foci of ECCO are on a global sea level changes, in particular contributions from polar ice sheets, and the increased use of biogeochemical and ecosystem data to constrain global cycles of carbon, nitrogen and oxygen. Challenges in the coming decade include provision of uncertainties, informing observing system design, globally increased resolution, and moving toward a coupled Earth system estimation with consistent momentum, heat and freshwater fluxes between the ocean, atmosphere, cryosphere and land.

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

2015
Edwards, CA, Moore AM, Hoteit I, Cornuelle BD.  2015.  Regional ocean data assimilation. Annual Review of Marine Science, Vol 7. 7:21-42.   10.1146/annurev-marine-010814-015821   AbstractWebsite

This article reviews the past 15 years of developments in regional ocean data assimilation. A variety of scientific, management, and safety-related objectives motivate marine scientists to characterize many ocean environments, including coastal regions. As in weather prediction, the accurate representation of physical, chemical, and/or biological properties in the ocean is challenging. Models and observations alone provide imperfect representations of the ocean state, but together they can offer improved estimates. Variational and sequential methods are among the most widely used in regional ocean systems, and there have been exciting recent advances in ensemble and four-dimensional variational approaches. These techniques are increasingly being tested and adapted for biogeochemical applications.

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

Gopalakrishnan, G, Cornuelle BD, Hoteit I.  2013.  Adjoint sensitivity studies of loop current and eddy shedding in the Gulf of Mexico. Journal of Geophysical Research: Oceans. 118:3315-3335.   10.1002/jgrc.20240   AbstractWebsite

Adjoint model sensitivity analyses were applied for the loop current (LC) and its eddy shedding in the Gulf of Mexico (GoM) using the MIT general circulation model (MITgcm). The circulation in the GoM is mainly driven by the energetic LC and subsequent LC eddy separation. In order to understand which ocean regions and features control the evolution of the LC, including anticyclonic warm-core eddy shedding in the GoM, forward and adjoint sensitivities with respect to previous model state and atmospheric forcing were computed using the MITgcm and its adjoint. Since the validity of the adjoint model sensitivities depends on the capability of the forward model to simulate the real LC system and the eddy shedding processes, a 5 year (2004–2008) forward model simulation was performed for the GoM using realistic atmospheric forcing, initial, and boundary conditions. This forward model simulation was compared to satellite measurements of sea-surface height (SSH) and sea-surface temperature (SST), and observed transport variability. Despite realistic mean state, standard deviations, and LC eddy shedding period, the simulated LC extension shows less variability and more regularity than the observations. However, the model is suitable for studying the LC system and can be utilized for examining the ocean influences leading to a simple, and hopefully generic LC eddy separation in the GoM. The adjoint sensitivities of the LC show influences from the Yucatan Channel (YC) flow and Loop Current Frontal Eddy (LCFE) on both LC extension and eddy separation, as suggested by earlier work. Some of the processes that control LC extension after eddy separation differ from those controlling eddy shedding, but include YC through-flow. The sensitivity remains stable for more than 30 days and moves generally upstream, entering the Caribbean Sea. The sensitivities of the LC for SST generally remain closer to the surface and move at speeds consistent with advection by the high-speed core of the current, while sensitivities to SSH generally extend to deeper layers and propagate more slowly. The adjoint sensitivity to relative vorticity deduced from the sensitivities to velocity fields suggests that advection of cyclonic (positive) relative vorticity anomalies from the YC or the LCFEs accelerate the LC eddy separation. Forward model perturbation experiments were performed to complement and check the adjoint sensitivity analysis as well as sampling the predictability and nonlinearity of the LC evolution. The model and its adjoint can be used in four-dimensional variational assimilation (4D-VAR) to produce dynamically consistent ocean state estimates for analysis and forecasts of the circulation of the GoM.

Gopalakrishnan, G, Cornuelle BD, Hoteit I, Rudnick DL, Owens BW.  2013.  State estimates and forecasts of the loop current in the Gulf of Mexico using the MITgcm and its adjoint. Journal of Geophysical Research: Oceans. 118:3292-3314.   10.1002/jgrc.20239   AbstractWebsite

An ocean state estimate has been developed for the Gulf of Mexico (GoM) using the MIT general circulation model and its adjoint. The estimate has been tested by forecasting loop current (LC) evolution and eddy shedding in the GoM. The adjoint (or four-dimensional variational) method was used to match the model evolution to observations by adjusting model temperature and salinity initial conditions, open boundary conditions, and atmospheric forcing fields. The model was fit to satellite-derived along-track sea surface height, separated into temporal mean and anomalies, and gridded sea surface temperature for 2 month periods. The optimized state at the end of the assimilation period was used to initialize the forecast for 2 months. Forecasts explore practical LC predictability and provide a cross-validation test of the state estimate by comparing it to independent future observations. The model forecast was tested for several LC eddy separation events, including Eddy Franklin in May 2010 during the deepwater horizon oil spill disaster in the GoM. The forecast used monthly climatological open boundary conditions, atmospheric forcing, and run-off fluxes. The model performance was evaluated by computing model-observation root-mean-square difference (rmsd) during both the hindcast and forecast periods. The rmsd metrics for the forecast generally outperformed persistence (keeping the initial state fixed) and reference (forecast initialized using assimilated Hybrid Coordinate Ocean Model 1/12° global analysis) model simulations during LC eddy separation events for a period of 1̃2 months.

Roux, P, Kuperman WA, Cornuelle BD, Aulanier F, Hodgkiss WS, Song HC.  2013.  Analyzing sound speed fluctuations in shallow water from group-velocity versus phase-velocity data representation. Journal of the Acoustical Society of America. 133:1945-1952.   10.1121/1.4792354   AbstractWebsite

Data collected over more than eight consecutive hours between two source-receiver arrays in a shallow water environment are analyzed through the physics of the waveguide invariant. In particular, the use of vertical arrays on both the source and receiver sides provides source and receiver angles in addition to travel-times associated with a set of eigenray paths in the waveguide. From the travel-times and the source-receiver angles, the eigenrays are projected into a group-velocity versus phase-velocity (Vg-Vp) plot for each acquisition. The time evolution of the Vg-Vp representation over the 8.5-h long experiment is discussed. Group speed fluctuations observed for a set of eigenrays with turning points at different depths in the water column are compared to the Brunt-Vaisala frequency. (C) 2013 Acoustical Society of America.

Song, HJ, Hoteit I, Cornuelle BD, Luo XD, Subramanian AC.  2013.  An adjoint-based adaptive ensemble Kalman filter. Monthly Weather Review. 141:3343-3359. AbstractWebsite

A new hybrid ensemble Kalman filter/four-dimensional variational data assimilation (EnKF/4D-VAR) approach is introduced to mitigate background covariance limitations in the EnKF. The work is based on the adaptive EnKF (AEnKF) method, which bears a strong resemblance to the hybrid EnKF/three-dimensional variational data assimilation (3D-VAR) method. In the AEnKF, the representativeness of the EnKF ensemble is regularly enhanced with new members generated after back projection of the EnKF analysis residuals to state space using a 3D-VAR [or optimal interpolation (OI)] scheme with a preselected background covariance matrix. The idea here is to reformulate the transformation of the residuals as a 4D-VAR problem, constraining the new member with model dynamics and the previous observations. This should provide more information for the estimation of the new member and reduce dependence of the AEnKF on the assumed stationary background covariance matrix. This is done by integrating the analysis residuals backward in time with the adjoint model. Numerical experiments are performed with the Lorenz-96 model under different scenarios to test the new approach and to evaluate its performance with respect to the EnKF and the hybrid EnKF/3D-VAR. The new method leads to the least root-mean-square estimation errors as long as the linear assumption guaranteeing the stability of the adjoint model holds. It is also found to be less sensitive to choices of the assimilation system inputs and parameters.

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

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

2010
Raghukumar, K, Cornuelle BD, Hodgkiss WS, Kuperman WA.  2010.  Experimental demonstration of the utility of pressure sensitivity kernels in time-reversal. Journal of the Acoustical Society of America. 128:989-1003.   10.1121/1.3466858   AbstractWebsite

Pressure sensitivity kernels were recently applied to time-reversal acoustics in an attempt to explain the enhanced stability of the time-reversal focal spot [Raghukumar et al., J. Acoust. Soc. Am. 124, 98-112 (2008)]. The theoretical framework developed was also used to derive optimized source functions, closely related to the inverse filter. The use of these optimized source functions results in an inverse filter-like focal spot which is more robust to medium sound speed fluctuations than both time-reversal and the inverse filter. In this paper the theory is applied to experimental data gathered during the Focused Acoustic Fields experiment, conducted in 2005, north of Elba Island in Italy. Sensitivity kernels are calculated using a range-independent sound-speed profile, for a geometry identical to that used in the experiment, and path sensitivities are identified with observed arrivals. The validity of the kernels in tracking time-evolving Green's functions is studied, along with limitations that result from a linearized analysis. An internal wave model is used to generate an ensemble of sound speed profiles, which are then used along with the calculated sensitivity kernels to derive optimized source functions. Focal spots obtained using the observed Green's functions with these optimized source functions are then compared to those obtained using time-reversal and the inverse-filter. It is shown that these functions are able to provide a focal spot superior to time-reversal while being more robust to sound speed fluctuations than the inverse filter or time-reversal. (C) 2010 Acoustical Society of America. [DOI: 10.1121/1.3466858]

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.

Song, H, Hoteit I, Cornuelle BD, Subramanian AC.  2010.  An Adaptive Approach to Mitigate Background Covariance Limitations in the Ensemble Kalman Filter. Monthly Weather Review. 138:2825-2845.   10.1175/2010mwr2871.1   AbstractWebsite

A new approach is proposed to address the background covariance limitations arising from undersampled ensembles and unaccounted model errors in the ensemble Kalman filter (EnKF). The method enhances the representativeness of the EnKF ensemble by augmenting it with new members chosen adaptively to add missing information that prevents the EnKF fromfully fitting the data to the ensemble. The vectors to be added are obtained by back projecting the residuals of the observation misfits from the EnKF analysis step onto the state space. The back projection is done using an optimal interpolation (OI) scheme based on an estimated covariance of the subspace missing from the ensemble. In the experiments reported here, the OI uses a stationary background covariance matrix, as in the hybrid EnKF-three-dimensional variational data assimilation (3DVAR) approach, but the resulting correction is included as a new ensemble member instead of being added to all existing ensemble members. The adaptive approach is tested with the Lorenz-96 model. The hybrid EnKF-3DVAR is used as a benchmark to evaluate the performance of the adaptive approach. Assimilation experiments suggest that the new adaptive scheme significantly improves the EnKF behavior when it suffers from small size ensembles and neglected model errors. It was further found to be competitive with the hybrid EnKF-3DVAR approach, depending on ensemble size and data coverage.

2009
Hoteit, I, Cornuelle B, Kim SY, Forget G, Kohl A, Terrill E.  2009.  Assessing 4D-VAR for dynamical mapping of coastal high-frequency radar in San Diego. Dynamics of Atmospheres and Oceans. 48:175-197.   10.1016/j.dynatmoce.2008.11.005   AbstractWebsite

The problem of dynamically mapping high-frequency (HF) radar radial velocity observations is investigated using a three-dimensional hydrodynamic model of the San Diego coastal region and an adjoint-based assimilation method. The HF radar provides near-real-time radial velocities from three sites covering the region offshore of San Diego Bay. The hydrodynamical model is the Massachusetts Institute of Technology general circulation model (MITgcm) with 1 km horizontal resolution and 40 vertical layers. The domain is centered on Point Loma, extending 117 km offshore and 120 km alongshore. The reference run (before adjustment) is initialized from a single profile of T and S and is forced with wind data from a single shore station and with zero heat and fresh water fluxes. The adjoint of the model is used to adjust initial temperature, salinity, and velocity, hourly temperature, salinity and horizontal velocities at the open boundaries, and hourly surface fluxes of momentum, heat and freshwater so that the model reproduces hourly HF radar radial velocity observations. Results from a small number of experiments suggest that the adjoint method can be successfully used over 10-day windows at coastal model resolution. It produces a dynamically consistent model run that fits HF radar data with errors near the specified uncertainties. In a test of the forecasting capability of the San Diego model after adjustment, the forecast skill was shown to exceed persistence for up to 20 h. (C) 2008 Elsevier B.V. All rights reserved.

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.

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.

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

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

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

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

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

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

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.

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.