Abstract. This paper presents the ECCO v4 non-linear inverse modeling framework and its baseline solution for the evolving ocean state over the period 1992–2011. Both components are publicly available and subjected to regular, automated regression tests. The modeling framework includes sets of global conformal grids, a global model setup, implementations of data constraints and control parameters, an interface to algorithmic differentiation, as well as a grid-independent, fully capable Matlab toolbox. The baseline ECCO v4 solution is a dynamically consistent ocean state estimate without unidentified sources of heat and buoyancy, which any interested user will be able to reproduce accurately. The solution is an acceptable fit to most data and has been found to be physically plausible in many respects, as documented here and in related publications. Users are being provided with capabilities to assess model–data misfits for themselves. The synergy between modeling and data synthesis is asserted through the joint presentation of the modeling framework and the state estimate. In particular, the inverse estimate of parameterized physics was instrumental in improving the fit to the observed hydrography, and becomes an integral part of the ocean model setup available for general use. More generally, a first assessment of the relative importance of external, parametric and structural model errors is presented. Parametric and external model uncertainties appear to be of comparable importance and dominate over structural model uncertainty. The results generally underline the importance of including turbulent transport parameters in the inverse problem.

Abstract. This paper presents the ECCO v4 non-linear inverse modeling framework and its baseline solution for the evolving ocean state over the period 1992–2011. Both components are publicly available and highly integrated with the MITgcm. They are both subjected to regular, automated regression tests. The modeling framework includes sets of global conformal grids, a global model setup, implementations of model-data constraints and adjustable control parameters, an interface to algorithmic differentiation, as well as a grid-independent, fully capable Matlab toolbox. The reference ECCO v4 solution is a dynamically consistent ocean state estimate (ECCO-Production, release 1) without un-identified sources of heat and buoyancy, which any interested user will be able to reproduce accurately. The solution is an acceptable fit to most data and has been found physically plausible in many respects, as documented here and in related publications. Users are being provided with capabilities to assess model-data misfits for themselves. The synergy between modeling and data synthesis is asserted through the joint presentation of the modeling framework and the state estimate. In particular, the inverse estimate of parameterized physics was instrumental in improving the fit to the observed hydrography, and becomes an integral part of the ocean model setup available for general use. More generally, a first assessment of the relative importance of external, parametric and structural model errors is presented. Parametric and external model uncertainties appear to be of comparable importance and dominate over structural model uncertainty. The results generally underline the importance of including turbulent transport parameters in the inverse problem.

[1] The spatial distribution and fate of riverine dissolved organic carbon (DOC) in the Arctic may be significant for the regional carbon cycle but are difficult to fully characterize using the sparse observations alone. Numerical models of the circulation and biogeochemical cycles of the region can help to interpret and extrapolate the data and may ultimately be applied in global change sensitivity studies. Here we develop and explore a regional, three-dimensional model of the Arctic Ocean in which, for the first time, we explicitly represent the sources of riverine DOC with seasonal discharge based on climatological field estimates. Through a suite of numerical experiments, we explore the distribution of DOC-like tracers with realistic riverine sources and a simple linear decay to represent remineralization through microbial degradation. The model reproduces the slope of the DOC-salinity relationship observed in the eastern and western Arctic basins when the DOC tracer lifetime is about 10 years, consistent with published inferences from field data. The new empirical parameterization of riverine DOC and the regional circulation and biogeochemical model provide new tools for application in both regional and global change studies.

In recent years, high-resolution ("eddying") global three-dimensional ocean general circulation models have begun to include astronomical tidal forcing alongside atmospheric forcing. Such models can carry an internal tide field with a realistic amount of nonstationarity, and an internal gravity wave continuum spectrum that compares more closely with observations as model resolution increases. Global internal tide and gravity wave models are important for understanding the three-dimensional geography of ocean mixing, for operational oceanography, and for simulating and interpreting satellite altimeter observations. Here we describe the most important technical details behind such models, including atmospheric forcing, bathymetry, astronomical tidal forcing, self-attraction and loading, quadratic bottom boundary layer drag, parameterized topographic internal wave drag, shallow-water tidal equations, and a brief summary of the theory of linear internal gravity waves. We focus on simulations run with two models, the HYbrid Coordinate Ocean Model (HYCOM) and the Massachusetts Institute of Technology general circulation model (MITgcm). We compare the modeled internal tides and internal gravity wave continuum to satellite altimeter observations, moored observational records, and the predictions of the Garrett-Munk (1975) internal gravity wave continuum spectrum. We briefly examine specific topics of interest, such as tidal energetics, internal tide nonstationarity, and the role of nonlinearities in generating the modeled internal gravity wave continuum. We also describe our first attempts at using a Kalman filter to improve the accuracy of tides embedded within a general circulation model. We discuss the challenges and opportunities of modeling stationary internal tides, non-stationary internal tides, and the internal gravity wave continuum spectrum for satellite altimetry and other applications. Introductionhis book chapter is about global modeling of oceanic internal tides and the oceanic internal gravity wave continuum. The chapter focuses on hydrodynamical modeling, rather than empirical modeling, of such motions. Due to the operational oceanography theme of the book in which this chapter resides, we focus on high-spatial-resolution numerical models run over relatively short time scales-i.e., simulations that could form the dynamical backbone of operational models-rather than on lower-resolution models run over decades or centuries for climate forecasting purposes. In this introductory section, after defining internal gravity waves and internal tides, we discuss the motivation for, requirements for, and history of global modeling of internal tides and the internal gravity wave continuum. A subsequent section focuses on the technical details underlying such models, such as atmospheric forcing, bathymetry, astronomical tidal forcing, self-attraction and loading, quadratic bottom boundary layer drag, parameterized topographic internal wave drag, shallow-water tidal equations, and a brief synopsis of internal wave theor...

We revisit the challenges and prospects for ocean circulation models following Griffies et al. (2010). Over the past decade, ocean circulation models evolved through improved understanding, numerics, spatial discretization, grid configurations, parameterizations, data assimilation, environmental monitoring, and process-level observations and modeling. Important large scale applications over the last decade are simulations of the Southern Ocean, the Meridional Overturning Circulation and its variability, and regional sea level change. Submesoscale variability is now routinely resolved in process models and permitted in a few global models, and submesoscale effects are parameterized in most global models. The scales where nonhydrostatic effects become important are beginning to be resolved in regional and process models. Coupling to sea ice, ice shelves, and high-resolution atmospheric models has stimulated new ideas and driven improvements in numerics. Observations have provided insight into turbulence and mixing around the globe and its consequences are assessed through perturbed physics models. Relatedly, parameterizations of the mixing and overturning processes in boundary layers and the ocean interior have improved. New diagnostics being used for evaluating models Fox-Kemper et al. Ocean Circulation Models alongside present and novel observations are briefly referenced. The overall goal is summarizing new developments in ocean modeling, including: how new and existing observations can be used, what modeling challenges remain, and how simulations can be used to support observations.