The HYCOM (HYbrid Coordinate Ocean Model) data assimilative system

Abstract

This article provides an overview of the effort centered on the HYbrid Coordinate Ocean Model (HYCOM) to develop an eddy-resolving, real-time global and basin-scale ocean hindcast, nowcast, and prediction system in the context of the Global Ocean Data Assimilation Experiment (GODAE). The main characteristics of HYCOM are first presented, followed by a description and assessment of the present near real-time Atlantic forecasting system. Regional/coastal applications are also discussed since an important attribute of the data assimilative HYCOM simulations is the capability to provide boundary conditions to regional and coastal models. The final section describes the steps taken toward the establishment of the fully global eddy-resolving HYCOM data assimilative system and discusses some of the difficulties associated with advanced data assimilation given the size of the problem.

Introduction

A broad partnership of institutions¹ is presently collaborating in developing and demonstrating the performance and application of eddy-resolving, real-time global and basin-scale ocean hindcast, nowcast, and prediction systems using the HYbrid Coordinate Ocean Model (HYCOM). The plan is to transition these systems for operational use by the U.S. Navy at the Naval Oceanographic Office (NAVOCEANO), Stennis Space Center, MS, and the Fleet Numerical Meteorology and Oceanography Center (FNMOC), Monterey, CA; and by NOAA at the National Centers for Environmental Prediction (NCEP), Washington, D.C. The partnership is also the eddy-resolving global ocean data assimilative system development effort that is sponsored by the U.S. component of the Global Ocean Data Assimilation Experiment (GODAE). GODAE is a coordinated international effort envisioning “a global system of observations, communications, modeling, and assimilation that will deliver regular, comprehensive information on the state of the oceans, in a way that will promote and engender wide utility and availability of this resource for maximum benefit to the community”. Three of the GODAE specific objectives are to apply state-of-the-art models and assimilation methods to produce short-range open ocean forecasts, boundary conditions to extend predictability of coastal and regional subsystems, and initial conditions for climate forecast models (GODAE Strategic Plan, International GODAE Steering Team, 2000). HYCOM development is the result of collaborative efforts among the University of Miami, the Naval Research Laboratory (NRL), and the Los Alamos National Laboratory (LANL), as part of the multi-institutional HYCOM Consortium for Data Assimilative Ocean Modeling funded by the National Ocean Partnership Program (NOPP) in 1999 to develop and evaluate a data assimilative hybrid isopycnal-sigma-pressure (generalized) coordinate ocean model (Bleck, 2002, Chassignet et al., 2003, Halliwell, 2004).

Numerical modeling studies over the past several decades have demonstrated advances in both model architecture and the availability of computational resources for the scientific community. Perhaps the most noticeable aspect of this progression has been the evolution from simulations on coarse-resolution horizontal/vertical grids outlining basins of simplified geometry and bathymetry and forced by idealized stresses, to fine-resolution simulations incorporating realistic coastal definition and bottom topography, forced by observational data on relatively short time scales (Hurlburt and Hogan, 2000, Smith et al., 2000, Chassignet and Garraffo, 2001). Traditional Ocean General Circulation Models (OGCMs) use a single coordinate type to represent the vertical, but recent model comparison exercises performed in Europe (DYnamics of North Atlantic MOdels—DYNAMO) (Willebrand et al., 2001) and in the U.S. (Data Assimilation and Model Evaluation Experiment—DAMÉE) (Chassignet et al., 2000) have shown that no single vertical coordinate–depth, density, or terrain-following σ-levels–can by itself be optimal everywhere in the ocean. These and earlier comparison studies (Chassignet et al., 1996, Roberts et al., 1996, Marsh et al., 1996) have shown that the models considered are able to simulate the large-scale characteristics of the oceanic circulation reasonably well, but that the interior water mass distribution and associated thermohaline circulation are strongly influenced by localized processes that are not represented equally by each model's vertical discretization. The choice of the vertical coordinate system is one of the most important aspects of an ocean model's design and practical issues of representation and parameterization are often directly linked to the vertical coordinate choice (Griffies et al., 2000). Currently, there are three main vertical coordinates in use, none of which provides universal utility. Hence, many developers have been motivated to pursue research into hybrid approaches. Isopycnal (density tracking) layers are best in the deep stratified ocean, z-levels (constant fixed depths) are best used to provide high vertical resolution near the surface within the mixed layer, and σ-levels are often the best choice in shallow coastal regions. HYCOM combines all three approaches and the optimal distribution is chosen at every time step. The model makes a dynamically smooth transition between the coordinate types via the continuity equation using the hybrid vertical coordinate generator.

The layout of the paper is as follows. First, in Section 2, an overview of the HYCOM characteristics is presented with model performance illustrated using non-data assimilative basin-scale and regional nested simulations. The near real-time North Atlantic Ocean data assimilative system in then introduced in Section 3 and its hindcast capabilities evaluated. In Section 4, issues associated with regional/coastal applications are introduced and discussed. Future development plans are presented in Section 5.