CO2 climate sensitivity and snow-sea-ice albedo parameterization in an atmospheric GCM coupled to a mixed-layer ocean model

The snow-sea-ice albedo parameterization in an atmospheric general circulation model (GCM), coupled to a simple mixed-layer ocean and run with an annual cycle of solar forcing, is altered from a version of the same model described by Washington and Meehl (1984). The model with the revised formulation is run to equilibrium for 1 × CO₂ and 2 × CO₂ experiments. The 1 ×CO₂ (control) simulation produces a global mean climate about 1° warmer than the original version, and sea-ice extent is reduced. The model with the altered parameterization displays heightened sensitivity in the global means, but the geographical patterns of climate change due to increased carbon dioxide (CO₂) are qualitatively similar. The magnitude of the climate change is affected, not only in areas directly influenced by snow and ice changes but also in other regions of the globe, including the tropics where sea-surface temperature, evaporation, and precipitation over the oceans are greater. With the less-sensitive formulation, the global mean surface air temperature increase is 3.5 °C, and the increase of global mean precipitation is 7.12%. The revised formulation produces a globally averaged surface air temperature increase of 4.04 °C and a precipitation increase of 7.25%, as well as greater warming of the upper tropical troposphere. Sensitivity of surface hydrology is qualitatively similar between the two cases with the larger-magnitude changes in the revised snow and ice-albedo scheme experiment. Variability of surface air temperature in the model is comparable to observations in most areas except at high latitudes during winter. In those regions, temporal variation of the sea-ice margin and fluctuations of snow cover dependent on the snow-ice-albedo formulation contribute to larger-than-observed temperature variability. This study highlights an uncertainty associated with results from current climate GCMs that use highly parameterized snow-sea-ice albedo schemes with simple mixed-layer ocean models.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Impact of global warming on the interannual and interdecadal climate modes in a coupled GCM

In this study, we investigated the impact of global warming on the variabilities of large-scale interannual and interdecadal climate modes and teleconnection patterns with two long-term integrations of the coupled general circulation model of ECHAM4/OPYC3 at the Max-Planck-Institute for Meteorology, Hamburg. One is the control (CTRL) run with fixed present-day concentrations of greenhouse gases. The other experiment is a simulation of transient greenhouse warming, named GHG run. In the GHG run the averaged geopotential height at 500?hPa is increased significantly, and a negative phase of the Pacific/North American (PNA) teleconnection-like distribution pattern is intensified. The standard deviation over the tropics (high latitudes) is enhanced (reduced) on the interdecadal time scales and reduced (enhanced) on the interannual time scales in the GHG run. Except for an interdecadal mode related to the Southern Oscillation (SO) in the GHG run, the spatial variation patterns are similar for different (interannual?+?interdecadal, interannual, and interdecadal) time scales in the GHG and CTRL runs. Spatial distributions of the teleconnection patterns on the interannual and interdecadal time scales in the GHG run are also similar to those in the CTRL run. But some teleconnection patterns show linear trends and changes of variances and frequencies in the GHG run. Apart from the positive linear trend of the SO, the interdecadal modulation to the El Niño/SO cycle is enhanced during the GHG 2040?～?2099. This is the result of an enhancement of the Walker circulation during that period. La Niña events intensify and El Niño events relatively weaken during the GHG 2070?～?2090. It is interesting to note that with increasing greenhouse gas concentrations the relation between the SO and the PNA pattern is reversed significantly from a negative to a positive correlation on the interdecadal time scales and weakened on the interannual time scales. This suggests that the increase of the greenhouse gas concentrations will trigger the nonstationary correlation between the SO and the PNA pattern both on the interdecadal and interannual time scales.     

The DMS-cloud albedo feedback effect: Greatly underestimated?

There are a number of ways by which the biosphere may counter any impetus for global warming that might be produced by the rising CO₂ content of earth's atmosphere. Evidence for one of these phenomena, the DMS-cloud feedback effect, is discussed in light of recent claims that it is not of sufficient strength to be of much importance.     

North Atlantic decadal regimes in a coupled GCM simulation

 The non-stationarity of the North Atlantic atmosphere-ocean coupling is investigated utilizing a long time integration of a coupled atmosphere-ocean general circulation model (GCM) and a consistent atmospheric experiment forced by the climatological sea surface temperature (SST) of the coupled GCM. The temporal behavior of the North Atlantic Oscillation (NAO) is non-stationary with two different decadal regimes being identified: (a) phases with enhanced (active) low-frequency variability of the NAO index are characterized by regional modes with a baroclinic Pacific-North America (PNA) and a dominant barotropic North Atlantic pattern; (b) in phases with reduced (passive) low-frequency variability a global mode connects tropics and midlatitudes. The characteristic space scales are similar in the coupled and the consistent atmospheric experiment; the time scales of the atmospheric eigenmodes are modified by ocean dynamics. In the active (passive) phase the corresponding atmospheric mode is reinforced by the North Atlantic (tropical Pacific) SST. Received: 15 September 2000 / Accepted: 30 March 2001     

Tropical Pacific internal atmospheric dynamics and resolution in a coupled GCM

A noise reduction technique, namely the interactive ensemble (IE) approach is adopted to reduce noise at the air–sea interface due to internal atmospheric dynamics in a state-of-the-art coupled general circulation model (CGCM). The IE technique uses multiple realization of atmospheric general circulation models coupled to a single ocean general circulation model. The ensembles mean fluxes from the atmospheric simulations are communicated to the ocean component. Each atmospheric simulation receives the same SST coming from the ocean component. The only difference among the atmospheric simulations comes from perturbed initial conditions, thus the atmospheric states are, in principle synoptically independent. The IE technique can be used to better understand the importance of weather noise forcing of natural variability such as El Niño Southern Oscillation (ENSO). To study the impact of weather noise and resolution in the context of a CGCM, two IE experiments are performed at different resolutions. Atmospheric resolution is an important issue since the noise statistics will depend on the spatial scales resolved. A simple formulation to extract atmospheric internal variability is presented. The results are compared to their respective control cases where internal atmospheric variability is left unchanged. The noise reduction has a major impact on the coupled simulation and the magnitude of this effect strongly depends on the horizontal resolution of the atmospheric component model. Specifically, applying the noise reduction technique reduces the overall climate variability more effectively at higher resolution. This suggests that “weather noise” is more important in sustaining climate variability as resolution increases. ENSO statistics, dynamics, and phase asymmetry are all modified by the noise reduction, in particular ENSO becomes more regular with less phase asymmetry when noise is reduced. All these effects are more marked for the higher resolution case. In contrast, ENSO frequency is unchanged by the reduction in the weather noise, but its phase-locking to the annual cycle is strongly dependent on noise and resolution. At low resolution the noise structure is similar to the signal, whereas the spatial structure of the noise deviates from the spatial structure of the signal as resolution increases. It is also suggested that event-to-event differences are largely driven by atmospheric noise as opposed to chaotic dynamics within the context of the large-scale coupled system, suggesting that there is a well-defined “canonical” event.     

Surface currents and equatorial thermocline in a coupled upper ocean-atmosphere GCM

The Oregon State University coupled upper ocean-atmosphere GCM is evaluated in terms of the simulated winds, ocean currents and thermocline depth variations. Although the zonal wind velocities in the model are underestimated by a factor of about three and the zonal current velocities are underestimated by a factor of about five, the model is seen to qualitatively simulate the major features of the gyral scale currents, and the phases of the seasonal variation of the principal equatorial currents are in reasonable agreement with observations. The simulated tropical currents are dominated by Ekman transport and the eastern boundary currents do not penetrate far enough equatorward, while the western boundary currents do not penetrate far enough poleward. The subtropical trade wind belt and the mid-latitude westerlies are displaced equatorward of observations; hence, the mid-latitude eastward currents, principally the Kuroshio-North Pacific Drift and the Gulf Stream-North Atlantic Current are displaced equatorward. In spite of these shortcomings the surface current simulation of this two-layer upper ocean model is comparable with that of other ocean GCMs of coarse resolution. The coupled model successfully simulates the deepening of the thermocline westward across Pacific as a consequence of the prevailing Walker circulation. The region of most intense simulated surface forcing is located in the western Pacific due to a southwestward displacement of the northeast trade winds relative to observations; hence the equatorial Pacific is dominated by eastward propagation of thermocline depth variations. The excessively strong Ekman divergence and upwelling in the western Pacific cools the local warm pool, while incorrectly simulated westerlies in the eastern Pacific suppress upwelling and inhibit cooling from below. These features reduce the simulated trans-Pacific sea-surface temperature gradient, weakening the Walker circulation and the anomalies associated with the simulated Southern Oscillation. Offprint requests to: KR Sperber     

Polar amplification in a coupled climate model with locked albedo

In recent years, a substantial reduction of the sea ice in the Arctic has been observed. At the same time, the near-surface air in this region is warming at a rate almost twice as large as the global average—this phenomenon is known as the Arctic amplification. The role of the ice-albedo feedback for the Arctic amplification is still a matter of debate. Here the effect of the surface-albedo feedback (SAF) was studied using a coupled climate model CCSM3 from the National Center for Atmospheric Research. Experiments, where the SAF was suppressed by locking the surface albedo in the entire coupled model system, were conducted. The results reveal polar temperature amplification when this model, with suppressed albedo, is forced by a doubling of the atmospheric CO₂ content. Comparisons with variable albedo experiments show that SAF amplifies the surface-temperature response in the Arctic area by about 33%, whereas the corresponding value for the global-mean surface temperature is about 15%. Even though SAF is an important process underlying excessive warming at high latitudes, the Arctic amplification is only 15% larger in the variable than in the locked-albedo experiments. It is found that an increase of water vapour and total cloud cover lead to a greenhouse effect, which is larger in the Arctic than at lower latitudes. This is expected to explain a part of the Arctic surface–air-temperature amplification.     

The effect of Arabian Sea optical properties on SST biases and the South Asian summer monsoon in a coupled GCM

This study examines the effect of seasonally varying chlorophyll on the climate of the Arabian Sea and South Asian monsoon. The effect of such seasonality on the radiative properties of the upper ocean is often a missing process in coupled general circulation models and its large amplitude in the region makes it a pertinent choice for study to determine any impact on systematic biases in the mean and seasonality of the Arabian Sea. In this study we examine the effects of incorporating a seasonal cycle in chlorophyll due to phytoplankton blooms in the UK Met Office coupled atmosphere-ocean GCM HadCM3. This is achieved by performing experiments in which the optical properties of water in the Arabian Sea—a key signal of the semi-annual cycle of phytoplankton blooms in the region—are calculated from a chlorophyll climatology derived from Sea-viewing Wide Field-of-View Sensor (SeaWiFS) data. The SeaWiFS chlorophyll is prescribed in annual mean and seasonally-varying experiments. In response to the chlorophyll bloom in late spring, biases in mixed layer depth are reduced by up to 50% and the surface is warmed, leading to increases in monsoon rainfall during the onset period. However when the monsoons are fully established in boreal winter and summer and there are strong surface winds and a deep mixed layer, biases in the mixed layer depth are reduced but the surface undergoes cooling. The seasonality of the response of SST to chlorophyll is found to depend on the relative depth of the mixed layer to that of the anomalous penetration depth of solar fluxes. Thus the inclusion of the effects of chlorophyll on radiative properties of the upper ocean acts to reduce biases in mixed layer depth and increase seasonality in SST.     

Processes governing the predictability of the Atlantic meridional overturning circulation in a coupled GCM

The processes that govern the predictability of decadal variations in the North Atlantic meridional overturning circulation (MOC) are investigated in a long control simulation of the ECHO-G coupled atmosphere?Cocean model. We elucidate the roles of local stochastic forcing by the atmosphere, and other potential ocean processes, and use our results to build a predictive regression model. The primary influence on MOC variability is found to come from air?Csea heat fluxes over the Eastern Labrador Sea. The maximum correlation between such anomalies and the variations in the MOC occurs at a lead time of 2?years, but we demonstrate that the MOC integrates the heat flux variations over a period of 10?years. The corresponding univariate regression model accounts for 74.5% of the interannual variability in the MOC (after the Ekman component has been removed). Dense anomalies to the south of the Greenland-Scotland ridge are also shown to precede the overturning variations by 4?C6?years, and provide a second predictor. With the inclusion of this second predictor the resulting regression model explains 82.8% of the total variance of the MOC. This final bivariate model is also tested during large rapid decadal overturning events. The sign of the rapid change is always well represented by the bivariate model, but the magnitude is usually underestimated, suggesting that other processes are also important for these large rapid decadal changes in the MOC.     

Interannual and interdecadal variabilities in the Pacific in an MRI coupled GCM

Interannual and interdecadal variabilities in the Pacific are investigated with a coupled atmosphere-ocean GCM developed at MRI, Japan. The model is run for 70 years with flux adjustments. The model shows interannual variability in the tropical Pacific which has several typical characteristics shared with the observed ENSO. A basin-scale feature of the principal SST variation for the ENSO time scale shows negative correlation in the central North Pacific with the tropical SST, similar to that of the observed one. Associated variation of the model atmosphere indicates an intensification of the Aleutian Low and a PNA-like teleconnection pattern as a response to the tropical warm SST anomaly. The ENSO time scale variability in the midlatitude ocean consists of the westward propagation of the subsurface temperature signal and the temperature variation within the shallow mixed layer forced by the anomalous atmospheric heat fluxes. For the interdecadal time scale, variation of the SST is simulated realistically with a geographical pattern similar to that for the ENSO time scale, but it has a larger relative amplitude in the northern Pacific. For the atmosphere, spatial structure of the variation in the interdecadal time scale is also similar to that in the ENSO time scale, but has smaller amplitude in the northern Pacific. Long oceanic spin-up time (>∼10 y) in the mid-high latitude, however, makes oceanic response in the interdecadal time scale larger than that in the ENSO time scale. The lagged-regression analysis for the ocean temperature variation relative to the wind stress variation indicates that interdecadal variation of the ocean subsurface at the mid-high latitudes is considered as enhanced ocean gyre spin-up process in response to the atmospheric circulation change at the mid-high latitudes, remotely forced by the interdecadal variation of the tropical SST. Received: 6 November 1995 / Accepted: 19 April 1996     

The surface heat flux feedback. Part II: direct and indirect estimates in the ECHAM4/OPA8 coupled GCM

The surface heat flux feedback in the Atlantic Ocean is estimated in the ECHAM4/OPA8 coupled model. The net heat flux feedback is negative everywhere, mostly ranging between 15 and 35 W m^-2 K^-1, but reaching up to 50 W m^-2 K^-1 in the tropics, so that it damps existing sea surface temperature anomalies. The bulk of it is due to the turbulent flux, although in the tropics the radiation feedback also strongly contributes. The turbulent heat flux feedback is strongest in fall and winter at extra-tropical latitudes, and in spring and summer near the equator. At mid-latitudes, the radiation feedback remains small in each season, but it can be strongly negative in parts of the tropics. At extra-tropical latitudes the model feedback compares rather well with estimates derived in Part I from the COADS observations and the NCEP reanalysis, but in the tropical Atlantic the negative heat flux feedback is much too strong. An indirect estimation of the model heat flux feedback is also attempted in regions of small mean surface current, based on the difference in decay time of sea surface temperature and salinity anomalies. The inferred negative heat flux feedback is qualitatively correct, but the seasonal changes in the mixed-layer depth are too large for the method to be accurate at high latitudes.     

Planktonic dimethylsulfide and cloud albedo: An estimate of the feedback response

Partial control of climate by the biosphere may be possible through a chain of processes that ultimately links marine plankton production of dimethylsulfide (DMS) with changes in cloud albedo (Charlson et al., 1987). Changes in cloud optical properties can have profound impacts on atmospheric radiation transfer and, hence, the surface environment. In this study, we have developed a simple model that incorporates empirically based parameterizations to account for the biological control of cloud droplet concentration in a first attempt to estimate the strength of the DMS-cloud albedo feedback mechanism. We find that the feedback reduces the global climatic response to imposed perturbations in solar insolation by less than 7%. Likewise, it modifies the strength of other feedbacks affecting surface insolation over oceans by roughly the same amount. This suggests that the DMS-cloud albedo mechanism will be unable to substantially reduce climate sensitivity, although these results should be confirmed with less idealized models when more is known about the net production of DMS by the marine biosphere and its relation to aerosol/cloud microphysics and climate.     

The influence of systematic errors in the Southeast Pacific on ENSO variability and prediction in a coupled GCM

The impact of the warm SST bias in the Southeast Pacific (SEP) on the quality of seasonal and interannual variability and ENSO prediction in a coupled GCM is investigated. The reduction of this bias is achieved by means of empirical heat flux correction that is constant in time. It leads to a wide range of changes in the tropical Pacific climate including enhanced southeast trades, well-defined dry zone in the SEP, better simulation of the South Pacific Convergence Zone and stronger cross-equatorial asymmetry of the mean state in the eastern Pacific. As a result of the mean climate correction, significant improvements in the simulation of the seasonal cycle of the oceanic and atmospheric states are also observed both at the equator and basin-wide. Due to more realistic simulation of the seasonal evolution of the cold tongue, tropical convection and surface winds in the corrected version of the model, phase-lock of ENSO to the annual cycle looses its strong semi-annual component and becomes quite similar to the observed, although the amplitude of ENSO is reduced. Zonal wind stress response to the SST anomalies in the central-eastern Pacific also becomes more realistic. ENSO retrospective forecast experiments conducted with the directly coupled and the flux-corrected versions of the model demonstrate that deficiencies in the seasonal evolution of the cold tongue/Inter-Tropical Convergence Zone complex (that were largely due to the SEP bias in this model) and the related errors in the ENSO phase-lock to the annual cycle can seriously degrade ENSO prediction. By reducing these errors, ENSO predictive skill in the coupled model was substantially enhanced.     

Decadal and centennial variability of the southern semiannual oscillation simulated in the GFDL coupled GCM

 We investigate the behavior of the semiannual oscillation (SAO) in surface pressure and 500 hPa baroclinicity at high southern latitudes in a 1000-year GFDL coupled ocean-atmosphere GCM run. The model represents this feature but is shown to underestimate its amplitude and percentage of variance explained in the midlatitudes. South of 60 °S the simulation of the pressure oscillation, although somewhat too weak, is considerably better. Our analysis reveals significant interannual, decadal and centennial variability in the modeled SAO. While there is only a short historical record of observational data in the middle and high southern latitudes, existing studies suggest that the strength of the SAO does show significant variability on at least the first two of these time scales. Strong relationships between the semiannual cycles of surface pressure and baroclinicity are apparent in the model output, reinforcing the findings in earlier studies that the differing annual march of temperature between the midlatitudes and the Antarctic coast leads to a semiannual component in the baroclinicity and thence the surface pressure. We draw attention to extended periods when the model SAO is weak and strong, and have investigated the nature of the circulation during each period. The GFDL model results suggest that a significant proportion of the SAO variation was associated more with variations in the strength of the winter pressure maximum rather than the springtime minimum. The extent to which this and other aspects of the modeled longterm variability are related to actual atmospheric structure must await the availability of longer data records. Received: 11 November 1996 / Accepted: 28 July 1997     

Sensitivity of the tropical Pacific to a change of orbital forcing in two versions of a coupled GCM

 The changes of the variability of the tropical Pacific ocean forced by a shift of six months in the date of the perihelion are studied using a coupled tropical Pacific ocean/global atmosphere GCM. The sensitivity experiments are conducted with two versions of the atmospheric model, varied by two parametrization changes. The first one concerns the interpolation scheme between the atmosphere and ocean models grids near the coasts, the second one the advection of water vapor in the presence of downstream negative temperature gradients, as encountered in the vicinity of mountains. In the tropical Pacific region, the parametrization differences only have a significant direct effect near the coasts; but coupled feedbacks lead to a 1 °C warming of the equatorial cold tongue in the modified (version 2) model, and a widening of the western Pacific large-scale convergence area. The sensitivity of the seasonal cycle of equatorial SST is very different between the two experiments. In both cases, the response to the solar flux forcing is strongly modified by coupled interactions between the SST, wind stress response and ocean dynamics. In the first version, the main feedback is due to anomalous upwelling and leads to westward propagation of SST anomalies; whereas the version 2 model is dominated by an eastward-propagating thermocline mode. The main reason diagnosed for these different behaviors is the atmospheric response to SST anomalies. In the warmer climate simulated by the second version, the wind stress response in the western Pacific is enhanced, and the off-equatorial curl is reduced, both effects favoring eastward propagation through thermocline depth anomalies. The modifications of the simulated seasonal cycle in version 2 lead to a change in ENSO behavior. In the control climate, the interannual variability in the eastern Pacific is dominated by warm events, whereas cold events tend to be the more extreme ones with a shifted perihelion. Received: 14 December 1999 / Accepted: 24 May 2000     

Time-dependent greenhouse warming computations with a coupled ocean-atmosphere model

Climate changes during the next 100 years caused by anthropogenic emissions of greenhouse gases have been simulated for the Intergovernmental Panel on Climate Change Scenarios A (business as usual) and D (accelerated policies) using a coupled ocean-atmosphere general circulation model. In the global average, the near-surface temperature rises by 2.6 K in Scenario A and by 0.6 K in Scenario D. The global patterns of climate change for both IPCC scenarios and for a third step-function 2 x CO₂ experiment were found to be very similar. The warming delay over the oceans is larger than found in simulations with atmospheric general circulation models coupled to mixed-layer models, leading to a more pronounced land-sea contrast and a weaker warming (and in some regions even an initial cooling) in the Southern Ocean. During the first forty years, the global warming and sea level rise due to the thermal expansion of the ocean are significantly slower than estimated previously from box-diffusion-upwelling models, but the major part of this delay can be attributed to the previous warming history prior to the start of present coupled ocean-atmosphere model integration (cold start).     

Sensitivity of the thermohaline circulation in coupled oceanic GCM — atmospheric EBM experiments

We analyze the sensitivity of the oceanic thermohaline circulation (THC) regarding perturbations in fresh water flux for a range of coupled oceanic general circulation — atmospheric energy balance models. The energy balance model (EBM) predicts surface air temperature and fresh water flux and contains the feedbacks due to meridional transports of sensible and latent heat. In the coupled system we examine a negative perturbation in run-off into the southern ocean and analyze the role of changed atmospheric heat transports and fresh water flux. With mixed boundary conditions (fixed air temperature and fixed surface fresh water fluxes) the response is characterized by a completely different oceanic heat transport than in the reference case. On the other hand, the surface heat flux remains roughly constant when the air temperature can adjust in a model where no anomalous atmospheric transports are allowed. This gives an artificially stable system with nearly unchanged oceanic heat transport. However, if meridional heat transports in the atmosphere are included, the sensitivity of the system lies between the two extreme cases. We find that changes in fresh water flux are unimportant for the THC in the coupled system.