Snow and cloud feedbacks modelled by an atmospheric general circulation model

One of the most important parametrizations in general circulation models used for climate change experiments is that of the surface albedo. The results of an albedo feedback experiment carried out under the auspices of the US Department of Energy are presented. An analysis of long and short wave components of the model response shows that short wave response dominates changes in fixed to variable albedo experiments, but that long wave response dominates in clear to cloudy sky changes. Cloud distribution changes are also discussed and are related to changes in global sensitivity. At the surface, the heat balance change for perturbed sea surface temperatures is dominated by changes in latent heat flux and downward long wave radiation. If albedo is freed up however, the major contrast lies in the change in surface reflected short wave radiation, amplified by changes in downward short wave radiation caused by cloud amount changes.     

Climate feedbacks in a general circulation model incorporating prognostic clouds

 This study performs a comprehensive feedback analysis on the Bureau of Meteorology Research Centre General Circulation Model, quantifying all important feedbacks operating under an increase in atmospheric CO₂. The individual feedbacks are analysed in detail, using an offline radiation perturbation method, looking at long- and shortwave components, latitudinal distributions, cloud impacts, non-linearities under 2xCO₂ and 4xCO₂ warmings and at interannual variability. The water vapour feedback is divided into terms due to moisture height and amount changes. The net cloud feedback is separated into terms due to cloud amount, height, water content, water phase, physical thickness and convective cloud fraction. Globally the most important feedbacks were found to be (from strongest positive to strongest negative) those due to water vapour, clouds, surface albedo, lapse rate and surface temperature. For the longwave (LW) response the most important term of the cloud ‘optical property’ feedbacks is due to the water content. In the shortwave (SW), both water content and water phase changes are important. Cloud amount and height terms are also important for both LW and SW. Feedbacks due to physical cloud thickness and convective cloud fraction are found to be relatively small. All cloud component feedbacks (other than height) produce conflicting LW/SW feedbacks in the model. Furthermore, the optical property and cloud fraction feedbacks are also of opposite sign. The result is that the net cloud feedback is the (relatively small) product of conflicting physical processes. Non-linearities in the feedbacks are found to be relatively small for all but the surface albedo response and some cloud component contributions. The cloud impact on non-cloud feedbacks is also discussed: greatest impact is on the surface albedo, but impact on water vapour feedback is also significant. The analysis method here proves to be a␣powerful tool for detailing the contributions from different model processes (and particularly those of the clouds) to the final climate model sensitivity. Received: 15 June 2000 / Accepted: 10 January 2001     

Evaluation of tropical cloud regimes in observations and a general circulation model

Tropical cloud regimes defined by cluster analysis of International Satellite Cloud Climatology Project (ISCCP) cloud top pressure (CTP)–optical thickness distributions and ISCCP-like Goddard Institute for Space Studies (GISS) general circulation model (GCM) output are analyzed in this study. The observations are evaluated against radar–lidar cloud-top profiles from the atmospheric radiation measurement (ARM) Program active remote sensing of cloud layers (ARSCL) product at two tropical locations and by placing them in the dynamical context of the Madden–Julian oscillation (MJO). ARSCL highest cloud-top profiles indicate that differences among some of the six ISCCP regimes may not be as prominent as suggested by ISCCP at the ARM tropical sites. An experimental adjustment of the ISCCP CTPs to produce cloud-top height profiles consistent with ARSCL eliminates the independence between those regimes. Despite these ambiguities, the ISCCP regime evolution over different phases of the MJO is consistent with existing MJO mechanisms, but with a greater mix of cloud types in each phase than is usually envisioned. The GISS Model E GCM produces two disturbed and two suppressed regimes when vertical convective condensate transport is included in the model’s cumulus parameterization. The primary model deficiencies are the absence of an isolated cirrus regime, a lack of mid-level cloud relative to ARSCL, and a tendency for occurrences of specific parameterized processes such as deep and shallow convection and stratiform low cloud formation to not be associated preferentially with any single cloud regime.     

A modelling study of aerosol processing by stratocumulus clouds and its impact on general circulation model parameterisations of cloud and aerosol

A model of the aqueous phase processing of an aerosol population undergoing multiple cycling through a stratocumulus (Sc) cloud layer is presented. Results indicate that a significant modification of the aerosol properties is achieved following the first cycle through cloud. In a polluted atmosphere, further modification in subsequent cycles is seen to be hydrogen peroxide limited unless there is a flux of ammonia entering the system through cloud base (CB). The modification of the aerosol population is seen to have little effect on the microphysics (specifically the cloud droplet concentration and effective radius) of the processing cloud. However, it enables processed aerosols to subsequently act as efficient cloud condensation nuclei (CCN) in less vigorous clouds (as a result of reducing the critical supersaturation required to activate them). The effects of variations in the internal mixture of soluble components of aerosols on the microphysics of clouds forming on them are also investigated using the cloud model. A (K2) parameterisation of the effects of variations in internally mixed nitrate loadings on the cloud droplet number concentration is presented. The effects of applying this K2 correction to the droplet number (derived from a parameterisation based on sulphate) for the presence of nitrate in aerosol have been investigated using the HadAM3 version of the Hadley Centre General Circulation Model (GCM). The effect on global annual mean simulations of the indirect forcing and effective radius is small, but more pronounced regionally. Suggestions (based on model results and observations) for parameterising the size distribution and in-cloud growth of aerosols for use in GCMs are presented.     

Meridional overturning circulation: stability and ocean feedbacks in a box model

A box model of the inter-hemispheric Atlantic meridional overturning circulation is developed, including a variable pycnocline depth for the tropical and subtropical regions. The circulation is forced by winds over a periodic channel in the south and by freshwater forcing at the surface. The model is aimed at investigating the ocean feedbacks related to perturbations in freshwater forcing from the atmosphere, and to changes in freshwater transport in the ocean. These feedbacks are closely connected with the stability properties of the meridional overturning circulation, in particular in response to freshwater perturbations. A separate box is used for representing the region north of the Antarctic circumpolar current in the Atlantic sector. The density difference between this region and the north of the basin is then used for scaling the downwelling in the north. These choices are essential for reproducing the sensitivity of the meridional overturning circulation observed in general circulation models, and therefore suggest that the southernmost part of the Atlantic Ocean north of the Drake Passage is of fundamental importance for the stability of the meridional overturning circulation. With this configuration, the magnitude of the freshwater transport by the southern subtropical gyre strongly affects the response of the meridional overturning circulation to external forcing. The role of the freshwater transport by the overturning circulation (M _ov ) as a stability indicator is discussed. It is investigated under which conditions its sign at the latitude of the southern tip of Africa can provide information on the existence of a second, permanently shut down, state of the overturning circulation in the box model. M _ov will be an adequate indicator of the existence of multiple equilibria only if salt-advection feedback dominates over other processes in determining the response of the circulation to freshwater anomalies. M _ov is a perfect indicator if feedbacks other than salt-advection are negligible.     

Design and performance of a new cloud microphysics scheme developed for the ECHAM general circulation model

A new cloud microphysics scheme including a prognostic treatment of cloud ice (PCI) is developed to yield a more physically based representation of the components of the atmospheric moisture budget in the general circulation model ECHAM. The new approach considers cloud water and cloud ice as separate prognostic variables. The precipitation formation scheme for warm clouds distinguishes between maritime and continental clouds by considering the cloud droplet number concentration, in addition to the liquid water content. Based on several observational data sets, the cloud droplet number concentration is derived from the sulfate aerosol mass concentration as given from the sulfur cycle simulated by ECHAM. Results obtained with the new scheme are compared to satellite observations and in situ measurements of cloud physical and radiative properties. In general, the standard model ECHAM4 and also PCI capture the overall features, and the simulated results usually lie within the range of observed uncertainty. As compared to ECHAM4, only slight improvements are achieved with the new scheme. For example, the overestimated liquid water path and total cloud cover over convectively active regions are reduced in PCI. On the other hand, some shortcomings of the standard model such as underestimated shortwave cloud forcing over the extratropical oceans of the respective summer hemisphere are more pronounced in PCI.This paper was presented at the Third International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 4–9 Sept. 1995 under the auspices of the Max Planck Institute for Meteorology, Hamburg. Editor for these papers is L. Dümenil.     

Fluctuation dissipation in a general circulation model

This paper considers the climate response to step function changes in the solar constant in two versions of a general circulation model with simplified geography. The NCAR CCM0 model is applied to an all-land planet with no topography (Terra Blanda). In one version there is moisture in the air (as well as self-generated clouds) as evaporated from an ideal surface at a fixed 80% of saturation. In the other version there is no moisture in the atmosphere. We examine the decay of natural anomalies in the large-scale temperature field in each model and compare the time dependence of the ensemble average with the average temporal behavior of the response to step function changes in the solar constant. The fluctuation-dissipation theorem of statistical mechanics makes specific predictions about the relationship between the two curves. We conduct the experiments for both versions of the model since the sensitivity is quite different for each. The theorem is found to hold reasonably well in each case.     

Cloud-radiation feedback and atmosphere-ocean coupling in a stochastic multicloud model

Despite recent advances in supercomputing, current general circulation models (GCMs) have significant problems in representing the variability associated with organized tropical convection. Furthermore, due to high sensitivity of the simulations to the cloud radiation feedback, the tropical convection remains a major source of uncertainty in long-term weather and climate forecasts. In a series of recent studies, it has been shown, in paradigm two-baroclinic-mode systems and in aquaplanet GCMs, that a stochastic multicloud convective parameterization based on three cloud types (congestus, deep and stratiform) can be used to improve the variability and the dynamical structure of tropical convection, including intermittent coherent structures such as synoptic and mesoscale convective systems. Here, the stochastic multicloud model is modified with a parameterized cloud radiation feedback mechanism and atmosphere-ocean coupling. The radiative convective feedback mechanism is shown to increase the mean and variability of the Walker circulation. The corresponding intensification of the circulation is associated with propagating synoptic scale systems originating inside of the enhanced sea surface temperature area. In column simulations, the atmosphere ocean coupling introduces pronounced low frequency convective features on the time scale associated with the depth of the mixed ocean layer. However, in the presence of the gravity wave mixing of spatially extended simulations, these features are not as prominent. This highlights the deficiency of the column model approach at predicting the behavior of multiscale spatially extended systems. Overall, the study develops a systematic framework for incorporating parameterized radiative cloud feedback and ocean coupling which may be used to improve representation of intraseasonal and seasonal variability in GCMs.     

Multi-year predictability in a coupled general circulation model

Multi-year to decadal variability in a 100-year integration of a BMRC coupled atmosphere-ocean general circulation model (CGCM) is examined. The fractional contribution made by the decadal component generally increases with depth and latitude away from surface waters in the equatorial Indo-Pacific Ocean. The relative importance of decadal variability is enhanced in off-equatorial “wings” in the subtropical eastern Pacific. The model and observations exhibit “ENSO-like” decadal patterns. Analytic results are derived, which show that the patterns can, in theory, occur in the absence of any predictability beyond ENSO time-scales. In practice, however, modification to this stochastic view is needed to account for robust differences between ENSO-like decadal patterns and their interannual counterparts. An analysis of variability in the CGCM, a wind-forced shallow water model, and a simple mixed layer model together with existing and new theoretical results are used to improve upon this stochastic paradigm and to provide a new theory for the origin of decadal ENSO-like patterns like the Interdecadal Pacific Oscillation and Pacific Decadal Oscillation. In this theory, ENSO-driven wind-stress variability forces internal equatorially-trapped Kelvin waves that propagate towards the eastern boundary. Kelvin waves can excite reflected internal westward propagating equatorially-trapped Rossby waves (RWs) and coastally-trapped waves (CTWs). CTWs have no impact on the off-equatorial sub-surface ocean outside the coastal wave guide, whereas the RWs do. If the frequency of the incident wave is too high, then only CTWs are excited. At lower frequencies, both CTWs and RWs can be excited. The lower the frequency, the greater the fraction of energy transmitted to RWs. This lowers the characteristic frequency (reddens the spectrum) of variability off the equator relative to its equatorial counterpart. At low frequencies, dissipation acts as an additional low pass filter that becomes more effective, as latitude increases. At the same time, ENSO-driven off-equatorial surface heating anomalies drive mixed layer temperature responses in both hemispheres. Both the eastern boundary interactions and the accumulation of surface heat fluxes by the surface mixed layer act to low pass filter the ENSO-forcing. The resulting off-equatorial variability is therefore more coherent with low pass filtered (decadal) ENSO indices [e.g. NINO3 sea-surface temperature (SST)] than with unfiltered ENSO indices. Consequently large correlations between variability and NINO3 extend further poleward on decadal time-scales than they do on interannual time-scales. This explains why decadal ENSO-like patterns have a broader meridional structure than their interannual counterparts. This difference in appearance can occur even if ENSO indices do not have any predictability beyond interannual time-scales. The wings around 15–20°S, and sub-surface variability at many other locations are predictable on interannual and multi-year time-scales. This includes westward propagating internal RWs within about 25° of the equator. The slowest of these take up to 4 years to reach the western boundary. This sub-surface predictability has significant oceanographic interest. However, it is linked to only low levels of SST variability. Consequently, extrapolation of delayed action oscillator theory to decadal time-scales might not be justified.     

Effect of seasonal forcing on global circulation in a world ocean general circulation model

 Effects of the seasonal variation in thermohaline and wind forcing on the abyssal circulation are investigated by using an ocean general circulation model. To isolate effects of the seasonality in the thermohaline forcing from those in the wind forcing, we carry out three experiments with (1) annual-mean wind forcing and perpetual-winter thermohaline forcing, (2) annual-mean wind forcing and seasonal thermohaline forcing, and (3) seasonal wind forcing and seasonal thermohaline forcing. The deep water under the seasonal thermohaline forcing becomes warmer than under the perpetual-winter thermohaline forcing. Although the perpetual-winter thermohaline forcing is widely used and believed to reproduce the deep water better than the annual-mean forcing, the difference between the results of the perpetual-winter and the seasonal thermohaline forcing is significant. The seasonal variation of the Ekman convergence and divergence produces meridional overturning cells extending to the bottom because the period of seasonal cycle is shorter than the adjustment timescale by baroclinic Rossby waves. The heat transport owing to those Ekman flows and temperature anomalies makes the upper water (0–200 m) colder at low to mid-latitudes (40S–40N) and warmer at high latitudes. Also the deep water becomes warmer owing to the warming of the northern North Atlantic, the main source region of North Atlantic Deep Water. The model is also synchronously (i.e., without acceleration) integrated with seasonal forcing for 5400 y. A past study suggested that under seasonal forcing, a sufficient equilibrium state can be achieved after only decades of synchronous integration following more than 10 000 y of accelerated integration. Here, the result so obtained is compared with that of the 5400-y synchronous integration. The difference in the global average temperature is as small as 0.12 °C, and most of the difference is confined to the Southern Ocean. Received: 1 May 1998 / Accepted: 5 January 1999     

Walker circulation controls ENSO atmospheric feedbacks in uncoupled and coupled climate model simulations

Many climate models strongly underestimate the two most important atmospheric feedbacks operating in El Niño/Southern Oscillation (ENSO), the positive (amplifying) zonal surface wind feedback and negative (damping) surface-heat flux feedback (hereafter ENSO atmospheric feedbacks, EAF). This hampers a realistic representation of ENSO dynamics in these models. Here we show that the atmospheric components of climate models participating in the 5th phase of the Coupled Model Intercomparison Project (CMIP5) when forced by observed sea surface temperatures (SST), already underestimate EAF on average by 23%, but less than their coupled counterparts (on average by 54%). There is a pronounced tendency of atmosphere models to simulate stronger EAF, when they exhibit a stronger mean deep convection and enhanced cloud cover over the western equatorial Pacific (WEP), indicative of a stronger rising branch of the Pacific Walker Circulation (PWC). Further, differences in the mean deep convection over the WEP between the coupled and uncoupled models explain a large part of the differences in EAF, with the deep convection in the coupled models strongly depending on the equatorial Pacific SST bias. Experiments with a single atmosphere model support the relation between the equatorial Pacific atmospheric mean state, the SST bias and the EAF. An implemented cold SST bias in the observed SST forcing weakens deep convection and reduces cloud cover in the rising branch of the PWC, causing weaker EAF. A warm SST bias has the opposite effect. Our results elucidate how biases in the mean state of the PWC and equatorial SST hamper a realistic simulation of the EAF.     

Simulation of thermohaline circulation with a twenty-layer oceanic general circulation model

Summary This paper presents the basic configuration and preliminary performance of a twenty-layer oceanic general circulation model which represents a portion of the recent progress in developing coupled ocean-atmosphere general circulation models made by the authors. The model uses latitude/depthdependent thermohaline-stratification subduction, -coordinate, three-dimensional implicit diffusion, complete convective adjustment, separating and coupling of external and internal modes and Asselin temporal filter, and thermodynamic sea-ice calculation. With seasonally varying climatological forcing at the surface and enhanced surface salinities in the region adjacent Antarctica, the model has been integrated for one thousand years to reach a quasiequilibrium state. Preliminary verification shows that the model is capable of simulating successfully not only many aspects of the upper ocean circulation but also an acceptable thermohaline circulation. The modelled overturning rate of the North Atlantic Deep Water (NADW) is greater than 15Sv. The simulated overturning rate of the Antarctic Bottom Water (AABW) is about 20Sv. The southward outflow of NADW can be identified from not only the meridional overturning streamfunction but also the current fields at four deeper levels from 1455m to 2475m. The AABW northward outflow exists at some bottom levels below 2600m, and mainly flows towards the Pacific basin.Major problems in the present simulation include the underestimate of the NADW outflow, the failure to simulate the Antarctic Intermediate Water (AAIW), the too fresh bottom water and the too diffuse thermocline of the model. A sensitivity experiment has revealed that the model diffusion process has an important impact on the simulation of both the thermocline and the NADW outflow.With 16 Figures     

Climate predictability experiments with a general circulation model

The atmospheric response to the evolution of the global sea surface temperatures from 1979 to 1992 is studied using the Max-Planck-Institut 19 level atmospheric general circulation model, ECHAM3 at T 42 resolution. Five separate 14-year integrations are performed and results are presented for each individual realization and for the ensemble-averaged response. The results are compared to a 30-year control integration using a climate monthly mean state of the sea surface temperatures and to analysis data. It is found that the ECHAM3 model, by and large, does reproduce the observed response pattern to El Nino and La Niña. During the El Nino events, the subtropical jet streams in both hemispheres are intensified and displaced equatorward, and there is a tendency towards weak upper easterlies over the equator. The Southern Oscillation is a very stable feature of the integrations and is accurately reproduced in all experiments. The inter-annual variability at middle- and high-latitudes, on the other hand, is strongly dominated by chaotic dynamics, and the tropical SST forcing only modulates the atmospheric circulation. The potential predictability of the model is investigated for six different regions. Signal to noise ratio is large in most parts of the tropical belt, of medium strength in the western hemisphere and generally small over the European area. The ENSO signal is most pronounced during the boreal spring. A particularly strong signal in the precipitation field in the extratropics during spring can be found over the southern United States. Western Canada is normally warmer during the warm ENSO phase, while northern Europe is warmer than normal during the ENSO cold phase. The reason is advection of warm air due to a more intense Pacific low than normal during the warm ENSO phase and a more intense Icelandic low than normal during the cold ENSO phase, respectively.     

Ocean-atmosphere interactions and climate drift in a coupled general circulation model

 We have analysed numerical simulations performed with a global 3D coupled atmosphere-ocean model to focus on the role of atmospheric processes leading to sea surface temperature (SST) drift in the tropics. Negative SST errors occur coherently in space and time with large positive errors in latent heat and momentum fluxes at the tropical air-sea interface, as diagnosed from forced SST simulations. The warm pool in the western Pacific disappears after a few years of simulation. Strong SST gradients enforce regions of high precipitation that are thin and stationary north of the equator. We detail the implications for the ocean-atmosphere system of such upheaval in the deep convection location. A sensitivity experiment to empirically formulate air-sea drag coefficient shows that the rapid warm pool erosion is not sensitive to changes in the formulation of the surface drag coefficient over the oceans because the corresponding changes in turbulent heat fluxes and LW cooling approximately cancel one another. In the eastern Pacific, the improvement in SST is striking and caused by feedbacks between SST, surface turbulent fluxes and boundary layer cloud fraction, which decreases as SST warms. Received: 8 December 1998 / Accepted: 6 January 2000     

Two climatic states and feedbacks on thermohaline circulation in an Earth system model of intermediate complexity

The McGill Paleoclimate Model-2 (MPM-2) is employed to study climate–thermohaline circulation (THC) interactions in a pre -industrial climate, with a special focus on the feedbacks on the THC from other climate system components. The MPM-2, a new version of the MPM, has an extended model domain from 90S to 90N, active winds and no oceanic heat and freshwater flux adjustments. In the MPM-2, there are mainly two stable modes for the Atlantic meridional overturning circulation (MOC) under the ‘present-day’ forcing (present-day solar forcing and the pre-industrial atmospheric CO₂ level of 280 ppm). The ‘on’ mode has an active North Atlantic deep water formation, while the ‘off’ mode has no such deep water formation. By comparing the ‘off’ mode climate state with its ‘on’ mode analogue, we find that there exist many large differences between the two climate states, which originate from large changes in the oceanic meridional heat transports. By suppressing or isolating each process associated with a continental ice sheet over North America, sea ice, the atmospheric hydrological cycle and vegetation, feedbacks from these components on the Atlantic MOC are investigated. Sensitivity studies investigating the role of varying continental ice growth and sea ice meridional transport in the resumption of the Atlantic MOC are also carried out. The results show that a fast ice sheet growth and an enhanced southward sea ice transport significantly favor the resumption of the Atlantic MOC in the MPM-2. In contrast to this, the feedback from the atmospheric hydrological cycle is a weak positive one. The vegetation-albedo feedback could enhance continental ice sheet growth and thus could also favor the resumption of the Atlantic MOC. However, before the shut-down of the Atlantic MOC, feedbacks from these components on the Atlantic MOC are very weak.     

Decadal climate variability simulated in a coupled general circulation model

A 1000 year integration of the CSIRO coupled ocean-atmosphere general circulation model is used to study low frequency (decadal to centennial) climate variability in precipitation and temperature. The model is shown to exhibit sizeable decadal variability for these fields, generally accounting for approximately 20 to 40% of the variability (greater than one year) in precipitation and up to 80% for temperature. An empirical orthogonal function (EOF) analysis is applied to the model output to show some of the major statistical modes of low frequency variability. The first EOF spatial pattern looks very much like that of the interannual ENSO pattern. It bears considerable resemblance to observational estimates and is centred in the Pacific extending into both hemispheres. It modulates both precipitation and temperature globally. The EOF has a time evolution that appears to be more than just red noise. Finally, the link between SST in the Pacific with Australian rainfall variability seen in observations is also evident in the model. Received: 29 August 1998 / Accepted: 31 July 1999     

A PDF-based hybrid prognostic cloud scheme for general circulation models

A new cloud parameterization based on prognostic equations for the subgrid-scale fluctuations in temperature and total water content is introduced for global climate models. The proposed scheme, called hybrid prognostic cloud (HPC) parameterization, employs simple probability density functions (PDFs) to the horizontal subgrid-scale inhomogeneity, allowing them to vary in shape in response to small-scale processes such as cumulus detrainment and turbulent mixing. Simple tests indicate that the HPC scheme is highly favorable as compared to a diagnostic scheme in terms of the cloud fraction and cloud water content under either uniform or non-uniform forcing. The relevance of the HPC scheme is investigated by implementing it in an atmospheric component model of the climate model MIROC with a coarse resolution of T42. A comparison of the short-term integrations between the T42 model and a global cloud resolving model (GCRM) reveals that the HPC scheme can reproduce, to a certain degree, the subgrid-scale variance and skewness of temperature and total water content simulated in the GCRM. It is also found that the HPC scheme significantly alters the climatological distributions in cloud cover, precipitation, and moisture, which are all improved from the model using a conventional diagnostic cloud scheme.