Showing papers by "Patrick Heimbach published in 2009"

Global General Circulation of the Ocean Estimated by the ECCO-Consortium

Abstract: Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use. The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.

Greenland ice-sheet volume sensitivity to basal, surface and initial conditions derived from an adjoint model

Abstract: We extend the application of control methods to a comprehensive three-dimensional thermomechanical ice-sheet model, SICOPOLIS (SImulation COde for POLythermal Ice Sheets). Lagrange multipliers, i.e. sensitivities, are computed with an exact, efficient adjoint model that has been generated from SICOPOLIS by rigorous application of automatic differentiation. The case study uses the adjoint model to determine the sensitivity of the total Greenland ice volume to various control variables over a 100 year period. The control space has of the order 1.2 × 10 6 elements, consisting of spatial fields of basal flow parameters, surface and basal forcings and initial conditions. Reliability of the adjoint model was tested through finite-difference perturbation calculations for various control variables and perturbation regions, ascertaining quantitative inferences of the adjoint model. As well as confirming qualitative aspects of ice-sheet sensitivities (e.g. expected regional variations), we detect regions where model sensitivities are seemingly unexpected or counter-intuitive, albeit 'real' in the sense of actual model behavior. An example is inferred regions where sensitivities of ice-sheet volume to basal sliding coefficient are positive, i.e. where a local increase in basal sliding parameter increases the ice- sheet volume. Similarly, positive (generally negative) ice temperature sensitivities in certain parts of the ice sheet are found, the detection of which seems highly unlikely if only conventional perturbation experiments had been used. The object of this paper is largely a proof of concept. Available adjoint-code generation tools now open up a variety of novel model applications, notably with regard to sensitivity and uncertainty analyses and ice-sheet state estimation or data assimilation.

Toward adjoinable MPI

Abstract: Automatic differentiation is the primary means of obtaining analytic derivatives from a numerical model given as a computer program. Therefore, it is an essential productivity tool in numerous computational science and engineering domains. Computing gradients with the adjoint (also called reverse) mode via source transformation is a particularly beneficial but also challenging use of automatic differentiation. To date only ad hoc solutions for adjoint differentiation of MPI programs have been available, forcing automatic differentiation tool users to reason about parallel communication dataflow and dependencies and manually develop adjoint communication code. Using the communication graph as a model we characterize the principal problems of adjoining the most frequently used communication idioms. We propose solutions to cover these idioms and consider the consequences for the MPI implementation, the MPI user and MPI-aware program analysis. The MIT general circulation model serves as a use case to illustrate the viability of our approach.

The Global Zonally Integrated Ocean Circulation, 1992–2006: Seasonal and Decadal Variability

Abstract: The zonally integrated meridional and vertical velocities as well as the enthalpy transports and fluxes in a least squares adjusted general circulation model are used to estimate the top-to-bottom oceanic meridional overturning circulation (MOC) and its variability from 1992 to 2006. A variety of simple theories all produce time scales suggesting that the mid- and high-latitude oceans should respond to atmospheric driving only over several decades. In practice, little change is seen in the MOC and associated heat transport except very close to the sea surface, at depth near the equator, and in parts of the Southern Ocean. Variability in meridional transports in both volume and enthalpy is dominated by the annual cycle and secondarily by the semiannual cycle, particularly in the Southern Ocean. On time scales longer than a year, the solution exhibits small trends with complicated global spatial patterns. Apart from a net uptake of heat from the atmosphere (forced by the NCEP–NCAR reanalysis, which...

Observed and simulated estimates of the meridional overturning circulation at 26.5° N in the Atlantic

Abstract: Daily timeseries of the meridional overturning circulation (MOC) estimated from the UK/US RAPID/MOCHA array at 26.5° N in the Atlantic are used to evaluate the MOC as simulated in two global circulation models: (I) an 8-member ensemble of the coupled climate model ECHAM5/MPI-OM, and (II) the ECCO-GODAE state estimate. In ECHAM5/MPI-OM, we find that the observed and simulated MOC have a similar variability and time-mean within the 99% confidence interval. In ECCO-GODAE, we find that the observed and simulated MOC show a significant correlation within the 99% confidence interval. To investigate the contribution of the different transport components, the MOC is decomposed into Florida Current, Ekman and mid-ocean transports. In both models, the mid-ocean transport is closely approximated by the residual of the MOC minus Florida Current and Ekman transports. As the models conserve volume by definition, future comparisons of the RAPID/MOCHA mid-ocean transport should be done against the residual transport in the models. The similarity in the variance and the correlation between the RAPID/MOCHA, and respectively ECHAM5/MPI-OM and ECCO-GODAE MOC estimates at 26.5° N is encouraging in the context of estimating (natural) variability in climate simulations and its use in climate change signal-to-noise detection analyses. Enhanced confidence in simulated hydrographic and transport variability will require longer observational time series.

Observed and simulated estimates of the meridional overturning circulation at 26.5 degrees N in the Atlantic

Abstract: Daily timeseries of the meridional overturning circulation (MOC) estimated from the UK/US RAPID/MOCHA array at 26.5 ◦ N in the Atlantic are used to evaluate the MOC as simulated in two global circulation models: (I) an 8-member ensemble of the coupled climate model ECHAM5/MPI-OM, and (II) the ECCO-GODAE state estimate. In ECHAM5/MPI-OM, we find that the observed and simulated MOC have a similar variability and time-mean within the 99% confidence interval. In ECCOGODAE, we find that the observed and simulated MOC show a significant correlation within the 99% confidence interval. To investigate the contribution of the different transport components, the MOC is decomposed into Florida Current, Ekman and mid-ocean transports. In both models, the midocean transport is closely approximated by the residual of the MOC minus Florida Current and Ekman transports. As the models conserve volume by definition, future comparisons of the RAPID/MOCHA mid-ocean transport should be done against the residual transport in the models. The similarity in the variance and the correlation between the RAPID/MOCHA, and respectively ECHAM5/MPI-OM and ECCO-GODAE MOC estimates at 26.5 ◦ N is encouraging in the context of estimating (natural) variability in climate simulations and its use in climate change signal-to-noise detection analyses. Enhanced confidence in simulated hydrographic and transport variability will require longer observational time series. Correspondence to: J. Baehr (johanna.baehr@zmaw.de)