The University of Southampton
University of Southampton Institutional Repository

The use of adjoint models for determining the sensitivity of integral quantities in an eddy resolving Oean General Circulation Model

The use of adjoint models for determining the sensitivity of integral quantities in an eddy resolving Oean General Circulation Model
The use of adjoint models for determining the sensitivity of integral quantities in an eddy resolving Oean General Circulation Model
Adjoint models calculate the exact sensitivity of an output function of a model to infinitesimal perturbations in the forcing or initial conditions. In eddy resolving ocean models the presence of chaotic eddies is expected to lead to sensitivities to infinitesimal perturbations that are very different from the sensitivity to large perturbations and that no longer contain useful information. Previous studies disagree as to whether adjoint models can be used with eddy resolving ocean models on timescales longer than a few months.

Here the MIT ocean general circulation model and its adjoint are used to look at the sensitivity of the time mean heat content, kinetic energy, available potential energy and thermocline depth to the sea surface temperature, zonal wind stress, and vertical diffusivity in an eddy resolving model of a zonally reentrant channel. Using the tangent linear model the non linear timescale of the eddy resolving model is estimated at around 200 days. The adjoint model is integrated over 278 days and 690 days to see whether useful information remains in the sensitivities calculated by the adjoint model for longer than the non linear timescale of the system. The usefulness of the information in the sensitivities calculated by the adjoint model is assessed by comparison with integrations of the full non linear forward model with large spatial scale perturbations to the forcing, finite difference gradient checks, and sensitivities calculated by an adjoint model in a non eddy resolving channel where the adjoint method is known to provide useful information.

Finite difference gradients are found to be unsuitable for calculating sensitivities of time averaged climate quantities in an eddy resolving ocean model as they are also affected by chaos. Comparison of the sensitivities calculated by the adjoint model in the eddy and non eddy resolving models shows that information remains in the spatial structure of the adjoint model results in the eddy resolving model on a time scale of 278 days.

In the non eddy resolving case the adjoint model results agree well with the perturbed forward model experiments, and are clearly climatically relevant on a timescale of 690 days. Use of a parameterisation scheme that reduces the eddy kinetic energy gives adjoint sensitivities that agree with well the perturbed forward model experiments after 690 days, although there are areas of extremely high adjoint sensitivity that may not be physically realistic. Without this parameterisation scheme, adjoint sensitivities involving dynamic variables grow exponentially with time as expected in a chaotic system, but at the end of the integration time of 690 days there is some agreement between the adjoint and forward model results for sensitivities involving thermodynamic variables only.

These results show that even in the presence of chaotic eddies some useful information is retained in the adjoint model solution beyond the nonlinear timescale of the system.
McLay, F.
a476a0cc-1bfc-45b3-86f1-680ef86cbe2c
McLay, F.
a476a0cc-1bfc-45b3-86f1-680ef86cbe2c

McLay, F. (2006) The use of adjoint models for determining the sensitivity of integral quantities in an eddy resolving Oean General Circulation Model. University of Southampton, Faculty of Engineering Science and Mathematics, School of Ocean and Earth Sciences, Doctoral Thesis, 92pp.

Record type: Thesis (Doctoral)

Abstract

Adjoint models calculate the exact sensitivity of an output function of a model to infinitesimal perturbations in the forcing or initial conditions. In eddy resolving ocean models the presence of chaotic eddies is expected to lead to sensitivities to infinitesimal perturbations that are very different from the sensitivity to large perturbations and that no longer contain useful information. Previous studies disagree as to whether adjoint models can be used with eddy resolving ocean models on timescales longer than a few months.

Here the MIT ocean general circulation model and its adjoint are used to look at the sensitivity of the time mean heat content, kinetic energy, available potential energy and thermocline depth to the sea surface temperature, zonal wind stress, and vertical diffusivity in an eddy resolving model of a zonally reentrant channel. Using the tangent linear model the non linear timescale of the eddy resolving model is estimated at around 200 days. The adjoint model is integrated over 278 days and 690 days to see whether useful information remains in the sensitivities calculated by the adjoint model for longer than the non linear timescale of the system. The usefulness of the information in the sensitivities calculated by the adjoint model is assessed by comparison with integrations of the full non linear forward model with large spatial scale perturbations to the forcing, finite difference gradient checks, and sensitivities calculated by an adjoint model in a non eddy resolving channel where the adjoint method is known to provide useful information.

Finite difference gradients are found to be unsuitable for calculating sensitivities of time averaged climate quantities in an eddy resolving ocean model as they are also affected by chaos. Comparison of the sensitivities calculated by the adjoint model in the eddy and non eddy resolving models shows that information remains in the spatial structure of the adjoint model results in the eddy resolving model on a time scale of 278 days.

In the non eddy resolving case the adjoint model results agree well with the perturbed forward model experiments, and are clearly climatically relevant on a timescale of 690 days. Use of a parameterisation scheme that reduces the eddy kinetic energy gives adjoint sensitivities that agree with well the perturbed forward model experiments after 690 days, although there are areas of extremely high adjoint sensitivity that may not be physically realistic. Without this parameterisation scheme, adjoint sensitivities involving dynamic variables grow exponentially with time as expected in a chaotic system, but at the end of the integration time of 690 days there is some agreement between the adjoint and forward model results for sensitivities involving thermodynamic variables only.

These results show that even in the presence of chaotic eddies some useful information is retained in the adjoint model solution beyond the nonlinear timescale of the system.

Text
McLay_2006_PhD.pdf - Other
Download (3MB)

More information

Published date: September 2006
Organisations: University of Southampton

Identifiers

Local EPrints ID: 46000
URI: http://eprints.soton.ac.uk/id/eprint/46000
PURE UUID: f97ba8cd-e061-4a28-90ca-af785eed7498

Catalogue record

Date deposited: 09 May 2007
Last modified: 15 Mar 2024 09:15

Export record

Contributors

Author: F. McLay

Download statistics

Downloads from ePrints over the past year. Other digital versions may also be available to download e.g. from the publisher's website.

View more statistics

Atom RSS 1.0 RSS 2.0

Contact ePrints Soton: eprints@soton.ac.uk

ePrints Soton supports OAI 2.0 with a base URL of http://eprints.soton.ac.uk/cgi/oai2

This repository has been built using EPrints software, developed at the University of Southampton, but available to everyone to use.

We use cookies to ensure that we give you the best experience on our website. If you continue without changing your settings, we will assume that you are happy to receive cookies on the University of Southampton website.

×