Patrick Heimbach

Bio: Patrick Heimbach is an academic researcher from University of Texas at Austin. The author has contributed to research in topics: Ocean current & Thermohaline circulation. The author has an hindex of 47, co-authored 153 publications receiving 7775 citations. Previous affiliations of Patrick Heimbach include Massachusetts Institute of Technology & Max Planck Society.

ECCO version 4: an integrated framework for non-linear inverse modeling and global ocean state estimation

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.

ECCO version 4: an integrated framework for non-linear inverse modeling and global ocean state estimation

Abstract: United States. National Aeronautics and Space Administration. Physical Oceanography Program

An Eddy-Permitting Southern Ocean State Estimate

Abstract: An eddy-permitting general circulation model of the Southern Ocean is fit by constrained least squares to a large observational dataset during 2005–06. Data used include Argo float profiles, CTD synoptic sections, Southern Elephant Seals as Oceanographic Samplers (SEaOS) instrument-mounted seal profiles, XBTs, altimetric observations [Envisat, Geosat, Jason-1, and Ocean Topography Experiment (TOPEX)/Poseidon], and infrared and microwave radiometer observed sea surface temperature. An adjoint model is used to determine descent directions in minimizing a misfit function, each of whose elements has been weighted by an estimate of the observational plus model error. The model is brought into near agreement with the data by adjusting its control vector, here consisting of initial and meteorological boundary conditions. Although total consistency has not yet been achieved, the existing solution is in good agreement with the great majority of the 2005 and 2006 Southern Ocean observations and better repr...

Global ocean circulation during 1992-1997, estimated from ocean observations and a general circulation model

Abstract: [1] A three-dimensional oceanic state is estimated for the period 1992–1997 as it results from combining large-scale ocean data sets with a general circulation model. At the cost of increased computational load, the estimation (assimilation) method is chosen specifically so that the resulting state estimate is consistent with the model equations, having no artificial sources or sinks. To bring the model into close agreement with observations, its initial temperature and salinity conditions are permitted to change, as are the time-dependent surface fluxes of momentum, heat and freshwater. Resulting changes of these “control vectors” are largely consistent with accepted uncertainties in the hydrographic climatology and meteorological analyses. The assimilation procedure is able to correct for many of the traditional shortcomings of the flow field by changing the surface boundary conditions. Changes in the resulting flow field are predominantly on the gyre scale and affect many features that are often poorly simulated in traditional numerical simulations, such as the strengths of the Gulf Stream and its extension, the Azores Current and the anticyclonic circulation associated with the Labrador Sea. Tests of the results and their consistency with prior error assumptions show that the constrained model has moved considerably closer to the observations imposed as constraints, but has also moved closer to independent data from the World Ocean Circulation Experiment not used in the assimilation procedure. In some regions where the comparisons remain indeterminate, not enough ocean observations are available, and it is difficult to ascribe the residuals to either the model or the observations. Although problems remain, a useful first solution to the global time-dependent ocean state estimation problem has been found. The estimates will continue to improve through the evolution of numerical models, computer power increases, more data, and more efficient estimation methods.

Climate Change 2007: The Physical Science Basis.

Abstract: Cause, conseguenze e strategie di mitigazione Proponiamo il primo di una serie di articoli in cui affronteremo l’attuale problema dei mutamenti climatici. Presentiamo il documento redatto, votato e pubblicato dall’Ipcc - Comitato intergovernativo sui cambiamenti climatici - che illustra la sintesi delle ricerche svolte su questo tema rilevante.

The Ensemble Kalman Filter: Theoretical formulation and practical implementation

Abstract: The purpose of this paper is to provide a comprehensive presentation and interpretation of the Ensemble Kalman Filter (EnKF) and its numerical implementation. The EnKF has a large user group, and numerous publications have discussed applications and theoretical aspects of it. This paper reviews the important results from these studies and also presents new ideas and alternative interpretations which further explain the success of the EnKF. In addition to providing the theoretical framework needed for using the EnKF, there is also a focus on the algorithmic formulation and optimal numerical implementation. A program listing is given for some of the key subroutines. The paper also touches upon specific issues such as the use of nonlinear measurements, in situ profiles of temperature and salinity, and data which are available with high frequency in time. An ensemble based optimal interpolation (EnOI) scheme is presented as a cost-effective approach which may serve as an alternative to the EnKF in some applications. A fairly extensive discussion is devoted to the use of time correlated model errors and the estimation of model bias.

Sea-Level Rise and Its Impact on Coastal Zones

Abstract: Global sea levels have risen through the 20th century. These rises will almost certainly accelerate through the 21st century and beyond because of global warming, but their magnitude remains uncertain. Key uncertainties include the possible role of the Greenland and West Antarctic ice sheets and the amplitude of regional changes in sea level. In many areas, nonclimatic components of relative sea-level change (mainly subsidence) can also be locally appreciable. Although the impacts of sea-level rise are potentially large, the application and success of adaptation are large uncertainties that require more assessment and consideration.

The development and evaluation of the Earth Gravitational Model 2008 (EGM2008)

Abstract: [1] EGM2008 is a spherical harmonic model of the Earth's gravitational potential, developed by a least squares combination of the ITG-GRACE03S gravitational model and its associated error covariance matrix, with the gravitational information obtained from a global set of area-mean free-air gravity anomalies defined on a 5 arc-minute equiangular grid This grid was formed by merging terrestrial, altimetry-derived, and airborne gravity data Over areas where only lower resolution gravity data were available, their spectral content was supplemented with gravitational information implied by the topography EGM2008 is complete to degree and order 2159, and contains additional coefficients up to degree 2190 and order 2159 Over areas covered with high quality gravity data, the discrepancies between EGM2008 geoid undulations and independent GPS/Leveling values are on the order of ±5 to ±10 cm EGM2008 vertical deflections over USA and Australia are within ±11 to ±13 arc-seconds of independent astrogeodetic values These results indicate that EGM2008 performs comparably with contemporary detailed regional geoid models EGM2008 performs equally well with other GRACE-based gravitational models in orbit computations Over EGM96, EGM2008 represents improvement by a factor of six in resolution, and by factors of three to six in accuracy, depending on gravitational quantity and geographic area EGM2008 represents a milestone and a new paradigm in global gravity field modeling, by demonstrating for the first time ever, that given accurate and detailed gravimetric data, asingle global model may satisfy the requirements of a very wide range of applications