Search Results
You are looking at 1 - 10 of 39 items for :
- Author or Editor: Jochem Marotzke x
- Article x
- Refine by Access: All Content x
Abstract
Boundary mixing is implemented in an ocean general circulation model such that the vertical mixing coefficient k
υ
is nonzero only near side boundaries and in convection regions. The model is used in a highly idealized configuration with no wind forcing and very nearly fixed surface density to investigate the three-dimensional dynamics of the thermohaline circulation. For k
υ
= 20 × 10−4 m2 s−1 and lower, the meridional overturning strength to great accuracy is proportional to
For k υ = 30 × 10−4 m2 s−1 and higher, no steady solutions have been found; the meridional overturning oscillates on a timescale of about 25 years. A time-averaged thermally direct overturning cell is not supported dynamically because convection extends longitudinally across the entire basin, and upwelling near the western boundary does not lead to densities higher than at the eastern boundary.
Assuming uniform upwelling in the west, level isopycnals near the equator, and level isopycnals along the eastern boundary south of the outcropping latitude permits the analytic determination of convection depth at the eastern wall and hence the density difference between the eastern and western walls. This difference is at most one-quarter the surface density difference between high and low latitudes, and agrees in magnitude and latitudinal dependence with the numerical experiments. Scaling arguments estimate overturning strength as of the order of 10 × 106 m3 s−1 and confirm the 2/3 power dependence on k υ . The derivation also gives a dependence of overturning strength with latitude that agrees qualitatively with the numerical results. The scaling for the dependence of meridional heat transport on latitude agrees well with the model results; scaling for heat transport amplitude agrees less well but correctly predicts a weaker dependence on k υ than maximum overturning.
Abstract
Boundary mixing is implemented in an ocean general circulation model such that the vertical mixing coefficient k
υ
is nonzero only near side boundaries and in convection regions. The model is used in a highly idealized configuration with no wind forcing and very nearly fixed surface density to investigate the three-dimensional dynamics of the thermohaline circulation. For k
υ
= 20 × 10−4 m2 s−1 and lower, the meridional overturning strength to great accuracy is proportional to
For k υ = 30 × 10−4 m2 s−1 and higher, no steady solutions have been found; the meridional overturning oscillates on a timescale of about 25 years. A time-averaged thermally direct overturning cell is not supported dynamically because convection extends longitudinally across the entire basin, and upwelling near the western boundary does not lead to densities higher than at the eastern boundary.
Assuming uniform upwelling in the west, level isopycnals near the equator, and level isopycnals along the eastern boundary south of the outcropping latitude permits the analytic determination of convection depth at the eastern wall and hence the density difference between the eastern and western walls. This difference is at most one-quarter the surface density difference between high and low latitudes, and agrees in magnitude and latitudinal dependence with the numerical experiments. Scaling arguments estimate overturning strength as of the order of 10 × 106 m3 s−1 and confirm the 2/3 power dependence on k υ . The derivation also gives a dependence of overturning strength with latitude that agrees qualitatively with the numerical results. The scaling for the dependence of meridional heat transport on latitude agrees well with the model results; scaling for heat transport amplitude agrees less well but correctly predicts a weaker dependence on k υ than maximum overturning.
Abstract
The length of time an ocean model and its adjoint should be integrated in determining a steady state compatible with observed data is investigated. The starting point is based upon a suggestion that only one time step is required. This method fails to converge to an acceptable solution when applied to a general circulation model (GCM) of the North Atlantic. Using a very coarse resolution GCM in an idealized geometry, the problem is traced to the interplay of convective adjustment and the very short integration time.
The general assimilation technique is explored using a very simple model, a linear first-order equation with forcing and damping. The model is unable to provide a dynamical coupling between the forcing and the model response, owing to a mismatch of integration time and adjustment time scale. Coupling can be enforced in the simple linear model through a careful choice of weighting factors, a strategy excluded in the GCM due to the presence of very fast processes like convective adjustment. An integration over a sufficiently long time can avoid the problems encountered. Experiments with the idealized GCM prove successful for longer integrations, and a tentative upper limit of 50 years is given for inversions aiming at the main thermocline structure.
Abstract
The length of time an ocean model and its adjoint should be integrated in determining a steady state compatible with observed data is investigated. The starting point is based upon a suggestion that only one time step is required. This method fails to converge to an acceptable solution when applied to a general circulation model (GCM) of the North Atlantic. Using a very coarse resolution GCM in an idealized geometry, the problem is traced to the interplay of convective adjustment and the very short integration time.
The general assimilation technique is explored using a very simple model, a linear first-order equation with forcing and damping. The model is unable to provide a dynamical coupling between the forcing and the model response, owing to a mismatch of integration time and adjustment time scale. Coupling can be enforced in the simple linear model through a careful choice of weighting factors, a strategy excluded in the GCM due to the presence of very fast processes like convective adjustment. An integration over a sufficiently long time can avoid the problems encountered. Experiments with the idealized GCM prove successful for longer integrations, and a tentative upper limit of 50 years is given for inversions aiming at the main thermocline structure.
Abstract
A general circulation model of the Indian Ocean is fitted to monthly averaged climatological temperatures, salinities, and surface fluxes using the adjoint method. Interannual variability is minimized by penalizing the temporal drift from one seasonal cycle to another during a two-year integration. The resultant meridional overturning and heat transport display large seasonal variations, with maximum amplitudes of 18 and 22 (× 106 m3 s−1) for the overturning and 1.8 and 1.4 (× 1015 W) for heat transport near 10°S and 10°N, respectively. A dynamical decomposition of the overturning and heat transport shows that the time-varying Ekman flow plus its barotropic compensation can explain a large part of the seasonal variations in overturning and heat transport. The maximum variations at 10°N and 10°S are associated with monsoon reversal over the northern Indian Ocean and changes of the easterlies over the southern Indian Ocean. An external mode with variable topography has a moderate contribution where the Somali Current and the corresponding gyre reverse direction seasonally. Contribution from vertical shear (thermal wind and ageostrophic shear) is dominant near the southern boundary and large near the Somali Current latitudes. The dominant balance in the zonally integrated heat budget is between heat storage change and heat transport convergence except south of 15°S.
Optimization with seasonal forcings improves estimates of sea surface temperatures, but the annual average overturning and heat transport are very similar to previous results with annual mean forcings. The annual average heat transport consists of roughly equal contributions from time-mean and time-varying fields of meridional velocities and temperatures in the northern Indian Ocean, indicating a significant rectification to the heat transport due to the time-varying fields. The time-mean and time-varying contributions are primarily due to the overturning and horizontal gyre, respectively.
Inclusion of TOPEX data enhances the seasonal cycles of the estimated overturning and heat transport in the central Indian Ocean significantly and improves the estimated equatorial zonal flows but leads to unrealistic estimates of the velocity structure near the Indonesian Throughflow region, most likely owing to the deficiencies in the lateral boundary conditions.
Abstract
A general circulation model of the Indian Ocean is fitted to monthly averaged climatological temperatures, salinities, and surface fluxes using the adjoint method. Interannual variability is minimized by penalizing the temporal drift from one seasonal cycle to another during a two-year integration. The resultant meridional overturning and heat transport display large seasonal variations, with maximum amplitudes of 18 and 22 (× 106 m3 s−1) for the overturning and 1.8 and 1.4 (× 1015 W) for heat transport near 10°S and 10°N, respectively. A dynamical decomposition of the overturning and heat transport shows that the time-varying Ekman flow plus its barotropic compensation can explain a large part of the seasonal variations in overturning and heat transport. The maximum variations at 10°N and 10°S are associated with monsoon reversal over the northern Indian Ocean and changes of the easterlies over the southern Indian Ocean. An external mode with variable topography has a moderate contribution where the Somali Current and the corresponding gyre reverse direction seasonally. Contribution from vertical shear (thermal wind and ageostrophic shear) is dominant near the southern boundary and large near the Somali Current latitudes. The dominant balance in the zonally integrated heat budget is between heat storage change and heat transport convergence except south of 15°S.
Optimization with seasonal forcings improves estimates of sea surface temperatures, but the annual average overturning and heat transport are very similar to previous results with annual mean forcings. The annual average heat transport consists of roughly equal contributions from time-mean and time-varying fields of meridional velocities and temperatures in the northern Indian Ocean, indicating a significant rectification to the heat transport due to the time-varying fields. The time-mean and time-varying contributions are primarily due to the overturning and horizontal gyre, respectively.
Inclusion of TOPEX data enhances the seasonal cycles of the estimated overturning and heat transport in the central Indian Ocean significantly and improves the estimated equatorial zonal flows but leads to unrealistic estimates of the velocity structure near the Indonesian Throughflow region, most likely owing to the deficiencies in the lateral boundary conditions.
Abstract
A general circulation model with a highly idealized geometry is used to investigate which fundamentally different equilibria of the global thermohaline circulation may exist. The model comprises two identical basins representing the Atlantic and Pacific oceans, which are connected by a circumpolar channel in the south. The model circulation is driven, in addition to wind forcing by restoring the sea surface temperature to prescribed values and specified freshwater fluxes in the surface salinity budget (mixed boundary conditions). The boundary conditions are symmetric with respect to the equator and identical for both oceans.
Four fundamentally different, stable steady states are found under the same set of boundary conditions. Two of the equilibria show both oceans in the same state, with high-altitude deep-water formation occuring either in both northern or in both southern oceans, respectively. Two additional equilibria exist in which the thermohaline circulations of the basins differ fundamentally from each other: one ocean forms deep water at northern high latitudes, while the other has a much weaker circulation with sinking in the Southern Hemisphere. One of these equilibria qualitatively corresponds to today's global thermohaline circulation pattern (conveyor belt).
It is demonstrated that a transition from one equilibrium to another can be accomplished by relatively small differences in the freshwater fluxes. The preference and sensitivity of the steady states depends critically on the freshwater forcing applied.
Abstract
A general circulation model with a highly idealized geometry is used to investigate which fundamentally different equilibria of the global thermohaline circulation may exist. The model comprises two identical basins representing the Atlantic and Pacific oceans, which are connected by a circumpolar channel in the south. The model circulation is driven, in addition to wind forcing by restoring the sea surface temperature to prescribed values and specified freshwater fluxes in the surface salinity budget (mixed boundary conditions). The boundary conditions are symmetric with respect to the equator and identical for both oceans.
Four fundamentally different, stable steady states are found under the same set of boundary conditions. Two of the equilibria show both oceans in the same state, with high-altitude deep-water formation occuring either in both northern or in both southern oceans, respectively. Two additional equilibria exist in which the thermohaline circulations of the basins differ fundamentally from each other: one ocean forms deep water at northern high latitudes, while the other has a much weaker circulation with sinking in the Southern Hemisphere. One of these equilibria qualitatively corresponds to today's global thermohaline circulation pattern (conveyor belt).
It is demonstrated that a transition from one equilibrium to another can be accomplished by relatively small differences in the freshwater fluxes. The preference and sensitivity of the steady states depends critically on the freshwater forcing applied.
Abstract
An idealized three-dimensional model of buoyancy-driven flow in a single hemisphere is used to investigate the relationship between the meridional overturning and the efficiency by which convective mixing eliminates static instability. In the “fast” limit (mixing timescale hours to weeks), the meridional overturning is not rate limited by the efficiency of convective mixing. If convective mixing is made less efficient, the model’s meridional overturning strength increases. Moreover, the dominant downwelling occurs not at the highest surface density;hence the deep ocean is relatively buoyant. The numerical results are explained by the different influences of convective mixing and downward advection on the deep-ocean heat budget; they underscore the fundamentally three-dimensional nature of the meridional overturning. In addition, the narrowness of deep downwelling is related to the geostrophic dynamics of deep temperature anomalies near the eastern wall. The model results presented here are in contrast to the expectation that deep-water formation by convective mixing is a necessary, if not rate-limiting, ingredient to the existence of a thermohaline circulation.
Abstract
An idealized three-dimensional model of buoyancy-driven flow in a single hemisphere is used to investigate the relationship between the meridional overturning and the efficiency by which convective mixing eliminates static instability. In the “fast” limit (mixing timescale hours to weeks), the meridional overturning is not rate limited by the efficiency of convective mixing. If convective mixing is made less efficient, the model’s meridional overturning strength increases. Moreover, the dominant downwelling occurs not at the highest surface density;hence the deep ocean is relatively buoyant. The numerical results are explained by the different influences of convective mixing and downward advection on the deep-ocean heat budget; they underscore the fundamentally three-dimensional nature of the meridional overturning. In addition, the narrowness of deep downwelling is related to the geostrophic dynamics of deep temperature anomalies near the eastern wall. The model results presented here are in contrast to the expectation that deep-water formation by convective mixing is a necessary, if not rate-limiting, ingredient to the existence of a thermohaline circulation.
Abstract
The wind-driven circulation adds a significant contribution to poleward meridional heat transport. Numerical models indicate that equatorward of ϕ 0, the zero wind stress latitude (30° lat), most of the wind-induced heat transport is due to the meridional overturning circulation known as the subtropical cell. The volume transport of this overturning is approximately given by the surface Ekman transport. By combining this fact with the assumption that Ekman-downwelled water approximately follows isotherms except near the equator, the authors derive an expression for the meridional heat transport that depends only on wind stress and surface temperature. The expression is confirmed in numerical models with simplified geometry and forcing. Numerical results indicate that peak heat transport due to the subtropical cell is about 0.1 × 1015 W for the North Atlantic and 0.3 × 1015 W for the North Pacific.
Abstract
The wind-driven circulation adds a significant contribution to poleward meridional heat transport. Numerical models indicate that equatorward of ϕ 0, the zero wind stress latitude (30° lat), most of the wind-induced heat transport is due to the meridional overturning circulation known as the subtropical cell. The volume transport of this overturning is approximately given by the surface Ekman transport. By combining this fact with the assumption that Ekman-downwelled water approximately follows isotherms except near the equator, the authors derive an expression for the meridional heat transport that depends only on wind stress and surface temperature. The expression is confirmed in numerical models with simplified geometry and forcing. Numerical results indicate that peak heat transport due to the subtropical cell is about 0.1 × 1015 W for the North Atlantic and 0.3 × 1015 W for the North Pacific.
Abstract
A coarse resolution, three-dimensional numerical model is used to study how external parameters control the existence and strength of equatorially asymmetric thermohaline overturning in a large-scale, rotating ocean basin. Initially, the meridional surface density gradient is directly set to be larger in a “dominant” hemisphere than in a “subordinate” hemisphere. The two-hemisphere system has a broader thermocline and weaker upwelling than the same model with the dominant hemisphere only. This behavior is in accord with classical scaling arguments, providing that the continuity equation is employed, rather than the linear vorticity equation.
The dominant overturning cell, analogous to North Atlantic Deep Water formation, is primarily controlled by the surface density contrast in the dominant hemisphere, which in turn is largely set by temperature. Consequently, in experiments with mixed boundary conditions, the dominant cell strength is relatively insensitive to the magnitude Q S of the salinity forcing. However, Q S strongly influences subordinate hemisphere properties, including the volume transport of a shallow overturning cell and the meridional extent of a tongue of low-salinity intermediate water reminiscent of Antarctic Intermediate Water.
The minimum Q S is identified for which the steady, asymmetric flow is stable; below this value, a steady, equatorially symmetric, temperature-dominated overturning occurs. For high salt flux, the asymmetric circulation becomes oscillatory and eventually gives way to an unsteady, symmetric, salt-dominated overturning. For given boundary conditions, it is possible to have at least three different asymmetric states, with significantly different large-scale properties. An expression for the meridional salt transport allows one to roughly predict the surface salinity and density profile and stability of the asymmetric state as a function of Q S and other external parameters.
Abstract
A coarse resolution, three-dimensional numerical model is used to study how external parameters control the existence and strength of equatorially asymmetric thermohaline overturning in a large-scale, rotating ocean basin. Initially, the meridional surface density gradient is directly set to be larger in a “dominant” hemisphere than in a “subordinate” hemisphere. The two-hemisphere system has a broader thermocline and weaker upwelling than the same model with the dominant hemisphere only. This behavior is in accord with classical scaling arguments, providing that the continuity equation is employed, rather than the linear vorticity equation.
The dominant overturning cell, analogous to North Atlantic Deep Water formation, is primarily controlled by the surface density contrast in the dominant hemisphere, which in turn is largely set by temperature. Consequently, in experiments with mixed boundary conditions, the dominant cell strength is relatively insensitive to the magnitude Q S of the salinity forcing. However, Q S strongly influences subordinate hemisphere properties, including the volume transport of a shallow overturning cell and the meridional extent of a tongue of low-salinity intermediate water reminiscent of Antarctic Intermediate Water.
The minimum Q S is identified for which the steady, asymmetric flow is stable; below this value, a steady, equatorially symmetric, temperature-dominated overturning occurs. For high salt flux, the asymmetric circulation becomes oscillatory and eventually gives way to an unsteady, symmetric, salt-dominated overturning. For given boundary conditions, it is possible to have at least three different asymmetric states, with significantly different large-scale properties. An expression for the meridional salt transport allows one to roughly predict the surface salinity and density profile and stability of the asymmetric state as a function of Q S and other external parameters.
Abstract
A one-dimensional radiative–convective equilibrium model is used to investigate the influence of clouds on the onset of a runaway greenhouse under strong solar forcing. By comparing experiments with clear-sky conditions (clouds are transparent to radiation) to experiments with full-sky conditions (clouds are radiatively active), the authors find that the critical solar irradiance that is necessary to trigger a runaway greenhouse is increased from around 1.15–1.20 times the present-day total solar irradiance (TSI) on Earth S 0 for clear-sky conditions to around 1.40–1.45S 0 for full-sky conditions. Cloud thickness increases with TSI, leading to a substantially higher albedo, which in turn allows the climate to remain in equilibrium for markedly higher values of TSI. The results suggest that steady states with sea surface temperatures higher than 335 K exist for a large range of TSI. The thickening clouds in these states do not reduce the outgoing longwave radiation any more, implying that the thickening of clouds increases only their shortwave effect. This mechanism allows the column to remain in balance even at high sea surface temperatures. The authors find double equilibria for both clear-sky and full-sky conditions, but the range for which they occur extends to considerably higher values of TSIs for full-sky conditions. Moreover, when clouds are included in the radiative transfer calculations, climate instabilities are no longer caused by longwave effects but by the cloud albedo effect.
Abstract
A one-dimensional radiative–convective equilibrium model is used to investigate the influence of clouds on the onset of a runaway greenhouse under strong solar forcing. By comparing experiments with clear-sky conditions (clouds are transparent to radiation) to experiments with full-sky conditions (clouds are radiatively active), the authors find that the critical solar irradiance that is necessary to trigger a runaway greenhouse is increased from around 1.15–1.20 times the present-day total solar irradiance (TSI) on Earth S 0 for clear-sky conditions to around 1.40–1.45S 0 for full-sky conditions. Cloud thickness increases with TSI, leading to a substantially higher albedo, which in turn allows the climate to remain in equilibrium for markedly higher values of TSI. The results suggest that steady states with sea surface temperatures higher than 335 K exist for a large range of TSI. The thickening clouds in these states do not reduce the outgoing longwave radiation any more, implying that the thickening of clouds increases only their shortwave effect. This mechanism allows the column to remain in balance even at high sea surface temperatures. The authors find double equilibria for both clear-sky and full-sky conditions, but the range for which they occur extends to considerably higher values of TSIs for full-sky conditions. Moreover, when clouds are included in the radiative transfer calculations, climate instabilities are no longer caused by longwave effects but by the cloud albedo effect.
Abstract
This study analyzes the response of the Atlantic meridional overturning circulation (AMOC) to different CO2 concentrations and two ice sheet configurations in simulations with the coupled climate model MPI-ESM. With preindustrial (PI) ice sheets, there are two different AMOC states within the studied CO2 range: one state with a strong and deep upper overturning cell at high CO2 concentrations and one state with a weak and shallow upper cell at low CO2 concentrations. Changes in AMOC variability with decreasing CO2 indicate two stability thresholds. The strong state is stable above the first threshold near 217 ppm, and the weak state is stable below the second threshold near 190 ppm. Between the two thresholds, both states are marginally unstable, and the AMOC oscillates between them on millennial time scales. The weak AMOC state is stable when Antarctic Bottom Water becomes dense and salty enough to replace North Atlantic Deep Water (NADW) in the deep North Atlantic and when the density gain over the North Atlantic becomes too weak to sustain continuous NADW formation. With Last Glacial Maximum (LGM) ice sheets, the density gain over the North Atlantic and the northward salt transport are enhanced with respect to the PI ice sheet case. This enables active NADW formation and a strong AMOC for the entire range of studied CO2 concentrations. The AMOC variability indicates that the simulated AMOC is far away from a stability threshold with LGM ice sheets. The nonlinear relationship among AMOC, CO2, and prescribed ice sheets provides an explanation for the large intermodel spread of AMOC states found in previous coupled LGM simulations.
Abstract
This study analyzes the response of the Atlantic meridional overturning circulation (AMOC) to different CO2 concentrations and two ice sheet configurations in simulations with the coupled climate model MPI-ESM. With preindustrial (PI) ice sheets, there are two different AMOC states within the studied CO2 range: one state with a strong and deep upper overturning cell at high CO2 concentrations and one state with a weak and shallow upper cell at low CO2 concentrations. Changes in AMOC variability with decreasing CO2 indicate two stability thresholds. The strong state is stable above the first threshold near 217 ppm, and the weak state is stable below the second threshold near 190 ppm. Between the two thresholds, both states are marginally unstable, and the AMOC oscillates between them on millennial time scales. The weak AMOC state is stable when Antarctic Bottom Water becomes dense and salty enough to replace North Atlantic Deep Water (NADW) in the deep North Atlantic and when the density gain over the North Atlantic becomes too weak to sustain continuous NADW formation. With Last Glacial Maximum (LGM) ice sheets, the density gain over the North Atlantic and the northward salt transport are enhanced with respect to the PI ice sheet case. This enables active NADW formation and a strong AMOC for the entire range of studied CO2 concentrations. The AMOC variability indicates that the simulated AMOC is far away from a stability threshold with LGM ice sheets. The nonlinear relationship among AMOC, CO2, and prescribed ice sheets provides an explanation for the large intermodel spread of AMOC states found in previous coupled LGM simulations.
Abstract
A theoretical analysis of the interactions between atmospheric meridional transports and the thermohaline circulation is presented, using a four-box ocean-atmosphere model in one hemisphere. The model is a simplified version of that developed by Nakamura Stone, and Marotzke and is amenable to analytical solutions. The ocean model is Stommel's; the atmospheric model gives the surface heat and freshwater fluxes as residuals of the atmospheric energy and moisture budgets, assumed in balance. Radiation at the top of the atmosphere depends linearly on surface temperature; atmospheric meridional heat and moisture transports are proportional to the meridional temperature gradient.
A Newtonian cooling law is derived for differential surface heat flux. The restoring coefficient is proportional to the efficiency of atmospheric transports and inversely proportional to the relative ocean area compared to total surface area. Surface freshwater flux increases with increasing temperature gradient and is inversely proportional to the ratio of ocean area to catchment area. The range of stable solutions with high-latitude sinking is smaller than in related, uncoupled box models due to the dependence of freshwater flux on the temperature gradient, which leads to a positive feedback with the thermohaline circulation. A strong control of the temperature gradient by atmospheric transports enhances the positive feedback between the salinity gradient and thermohaline Circulation simultaneously, it weakens the positive feedback between atmospheric moisture transport and the thermohaline circulation.
Overestimating the atmospheric moisture transport and underestimating oceanic mass transport both artificially destabilize the high-latitude sinking state. Overestimating the atmospheric heat transport and hence the Newtonian restoring coefficient can be artificially stabilizing or destabilizing. These erroneous sensitivities ate not alleviated if flux adjustments are added to obtain the correct mean climate, and then held fixed in climate change experiments. We derive alternate flux adjustment schemes, which do preserve the model's stability properties for particular sources of error.
Abstract
A theoretical analysis of the interactions between atmospheric meridional transports and the thermohaline circulation is presented, using a four-box ocean-atmosphere model in one hemisphere. The model is a simplified version of that developed by Nakamura Stone, and Marotzke and is amenable to analytical solutions. The ocean model is Stommel's; the atmospheric model gives the surface heat and freshwater fluxes as residuals of the atmospheric energy and moisture budgets, assumed in balance. Radiation at the top of the atmosphere depends linearly on surface temperature; atmospheric meridional heat and moisture transports are proportional to the meridional temperature gradient.
A Newtonian cooling law is derived for differential surface heat flux. The restoring coefficient is proportional to the efficiency of atmospheric transports and inversely proportional to the relative ocean area compared to total surface area. Surface freshwater flux increases with increasing temperature gradient and is inversely proportional to the ratio of ocean area to catchment area. The range of stable solutions with high-latitude sinking is smaller than in related, uncoupled box models due to the dependence of freshwater flux on the temperature gradient, which leads to a positive feedback with the thermohaline circulation. A strong control of the temperature gradient by atmospheric transports enhances the positive feedback between the salinity gradient and thermohaline Circulation simultaneously, it weakens the positive feedback between atmospheric moisture transport and the thermohaline circulation.
Overestimating the atmospheric moisture transport and underestimating oceanic mass transport both artificially destabilize the high-latitude sinking state. Overestimating the atmospheric heat transport and hence the Newtonian restoring coefficient can be artificially stabilizing or destabilizing. These erroneous sensitivities ate not alleviated if flux adjustments are added to obtain the correct mean climate, and then held fixed in climate change experiments. We derive alternate flux adjustment schemes, which do preserve the model's stability properties for particular sources of error.
