This paper discusses the errors in surface tracer and flux fields in ocean
models induced by using approximate surface boundary conditions involving r
elaxation toward observed values rather than more physically realistic cond
itions that involve (often inaccurate) surface fluxes. The authors show the
oretically and with a global model example that where there is a net annual
surface flux of tracer (balanced by advection), (i) the annual mean surfac
e tracer held is biased compared with the observations and (ii) the annual
mean tracer Flux is also biased if the surface tracer field has a feedback
on the surface tracer advection or diffusion. As previously shown, the ampl
itude of the annual cycle of tracers is also decreased. The global model in
dicates that temperature offsets of 1 degrees-2 degrees C (or even greater)
and heat flux errors of 30 W m(-2) occur in regions of strong advection, s
uch as the equatorial upwelling zone, western boundary currents, and the An
tarctic Circumpolar Current. These are all areas crucial for the thermohali
ne circulation, so that the use of such boundary conditions is likely to yi
eld incorrect estimates for climate simulation models. Zonally integrated m
eridional hear fluxes may be in error by up to 25%.