Impact of stratospheric variability on tropospheric climate change

An improved stratospheric representation has been included in simulations with the Hadley Centre HadGEM1 coupled ocean atmosphere model with natural and anthropogenic forcings for the period 1979–2003. An improved stratospheric ozone dataset is employed that includes natural variations in ozone as well as the usual anthropogenic trends. In addition, in a second set of simulations the quasi biennial oscillation (QBO) of stratospheric equatorial zonal wind is also imposed using a relaxation towards ERA-40 zonal wind values. The resulting impact on tropospheric variability and trends is described. We show that the modelled cooling rate at the tropopause is enhanced by the improved ozone dataset and this improvement is even more marked when the QBO is also included. The same applies to warming trends in the upper tropical troposphere which are slightly reduced. Our stratospheric improvements produce a significant increase of internal variability but no change in the positive trend of annual mean global mean near-surface temperature. Warming rates are increased significantly over a large portion of the Arctic Ocean. The improved stratospheric representation, especially the QBO relaxation, causes a substantial reduction in near-surface temperature and precipitation response to the El Chichón eruption, especially in the tropical region. The winter increase in the phase of the northern annular mode observed in the aftermath of the two major recent volcanic eruptions is partly captured, especially after the El Chichón eruption. The positive trend in the southern annular mode (SAM) is increased and becomes statistically significant which demonstrates that the observed increase in the SAM is largely subject to internal variability in the stratosphere. The possible inclusion in simulations for future assessments of full ozone chemistry and a gravity wave scheme to internally generate a QBO is discussed.

Transport in the low-latitude tropopause zone diagnosed using particle trajectories

Recent literature has described a “transition zone” between the average top of deep convection in the Tropics and the stratosphere. Here transport across this zone is investigated using an offline trajectory model. Particles were advected by the resolved winds from the European Centre for Medium-Range Weather Forecasts reanalyses. For each boreal winter clusters of particles were released in the upper troposphere over the four main regions of tropical deep convection (Indonesia, central Pacific, South America, and Africa). Most particles remain in the troposphere, descending on average for every cluster. The horizontal components of 5-day trajectories are strongly influenced by the El Niño–Southern Oscillation (ENSO), but the Lagrangian average descent does not have a clear ENSO signature. Tropopause crossing locations are first identified by recording events when trajectories from the same release regions cross the World Meteorological Organization lapse rate tropopause. Most crossing events occur 5–15 days after release, and 30-day trajectories are sufficiently long to estimate crossing number densities. In a further two experiments slight excursions across the lapse rate tropopause are differentiated from the drift deeper into the stratosphere by defining the “tropopause zone” as a layer bounded by the average potential temperature of the lapse rate tropopause and the profile temperature minimum. Transport upward across this zone is studied using forward trajectories released from the lower bound and back trajectories arriving at the upper bound. Histograms of particle potential temperature (θ) show marked differences between the transition zone, where there is a slow spread in θ values about a peak that shifts slowly upward, and the troposphere below 350 K. There forward trajectories experience slow radiative cooling interspersed with bursts of convective heating resulting in a well-mixed distribution. In contrast θ histograms for back trajectories arriving in the stratosphere have two distinct peaks just above 300 and 350 K, indicating the sharp change from rapid convective heating in the well-mixed troposphere to slow ascent in the transition zone. Although trajectories slowly cross the tropopause zone throughout the Tropics, all three experiments show that most trajectories reaching the stratosphere from the lower troposphere within 30 days do so over the west Pacific warm pool. This preferred location moves about 30°–50° farther east in an El Niño year (1982/83) and about 30° farther west in a La Niña year (1988/89). These results could have important implications for upper-troposphere–lower-stratosphere pollution and chemistry studies.

Sting jets in simulations of a real cyclone by two mesoscale models

The existence of sting jets as a potential source of damaging surface winds during the passage of extratropical cyclones has recently been recognized However, there are still very few published studies on the subject Furthermore, although ills known that other models are capable of reproducing sting jets, in the published literature only one numerical model [the Met Office Unified Model (MetUM)] has been used to numerically analyze these phenomena This article alms to improve our understanding of the processes that contribute to the development of sting jets and show that model differences affect the evolution of modeled sting jets A sting jet event during the passage of a cyclone over the United Kingdom on 26 February 2002 has been simulated using two mesoscale models namely the MetUM and the Consortium for Small Scale Modeling (COSMO) model to compare their performance Given the known critical importance of vertical resolution in the simulation of sting jets the vertical resolution of both models has been enhanced with respect to their operational versions Both simulations have been verified against surface measurements of maximum gusts, satellite imagery and Met Office operational synoptic analyses, as well as operational analyses from the ECMWF It is shown that both models are capable of reproducing sting jets with similar, though not identical. features Through the comparison of the results from these two models, the relevance of physical mechanisms, such as evaporative cooling and the release of conditional symmetric instability, in the generation and evolution of sting jets is also discussed

Report of the 2010 assessment of the scientific assessment panel: future ozone and its impact on surface UV

SCIENTIFIC SUMMARY Globally averaged total column ozone has declined over recent decades due to the release of ozone-depleting substances (ODSs) into the atmosphere. Now, as a result of the Montreal Protocol, ozone is expected to recover from the effects of ODSs as ODS abundances decline in the coming decades. However, a number of factors in addition to ODSs have led to and will continue to lead to changes in ozone. Discriminating between the causes of past and projected ozone changes is necessary, not only to identify the progress in ozone recovery from ODSs, but also to evaluate the effectiveness of climate and ozone protection policy options. Factors Affecting Future Ozone and Surface Ultraviolet Radiation • At least for the next few decades, the decline of ODSs is expected to be the major factor affecting the anticipated increase in global total column ozone. However, several factors other than ODS will affect the future evolution of ozone in the stratosphere. These include changes in (i) stratospheric circulation and temperature due to changes in long-lived greenhouse gas (GHG) abundances, (ii) stratospheric aerosol loading, and (iii) source gases of highly reactive stratospheric hydrogen and nitrogen compounds. Factors that amplify the effects of ODSs on ozone (e.g., stratospheric aerosols) will likely decline in importance as ODSs are gradually eliminated from the atmosphere. • Increases in GHG emissions can both positively and negatively affect ozone. Carbon dioxide (CO2)-induced stratospheric cooling elevates middle and upper stratospheric ozone and decreases the time taken for ozone to return to 1980 levels, while projected GHG-induced increases in tropical upwelling decrease ozone in the tropical lower stratosphere and increase ozone in the extratropics. Increases in nitrous oxide (N2O) and methane (CH4) concentrations also directly impact ozone chemistry but the effects are different in different regions. • The Brewer-Dobson circulation (BDC) is projected to strengthen over the 21st century and thereby affect ozone amounts. Climate models consistently predict an acceleration of the BDC or, more specifically, of the upwelling mass flux in the tropical lower stratosphere of around 2% per decade as a consequence of GHG abundance increases. A stronger BDC would decrease the abundance of tropical lower stratospheric ozone, increase poleward transport of ozone, and could reduce the atmospheric lifetimes of long-lived ODSs and other trace gases. While simulations showing faster ascent in the tropical lower stratosphere to date are a robust feature of chemistry-climate models (CCMs), this has not been confirmed by observations and the responsible mechanisms remain unclear. • Substantial ozone losses could occur if stratospheric aerosol loading were to increase in the next few decades, while halogen levels are high. Stratospheric aerosol increases may be caused by sulfur contained in volcanic plumes entering the stratosphere or from human activities. The latter might include attempts to geoengineer the climate system by enhancing the stratospheric aerosol layer. The ozone losses mostly result from enhanced heterogeneous chemistry on stratospheric aerosols. Enhanced aerosol heating within the stratosphere also leads to changes in temperature and circulation that affect ozone. • Surface ultraviolet (UV) levels will not be affected solely by ozone changes but also by the effects of climate change and by air quality change in the troposphere. These tropospheric effects include changes in clouds, tropospheric aerosols, surface reflectivity, and tropospheric sulfur dioxide (SO2) and nitrogen dioxide (NO2). The uncertainties in projections of these factors are large. Projected increases in tropospheric ozone are more certain and may lead to reductions in surface erythemal (“sunburning”) irradiance of up to 10% by 2100. Changes in clouds may lead to decreases or increases in surface erythemal irradiance of up to 15% depending on latitude. Expected Future Changes in Ozone Full ozone recovery from the effects of ODSs and return of ozone to historical levels are not synonymous. In this chapter a key target date is chosen to be 1980, in part to retain the connection to previous Ozone Assessments. Noting, however, that decreases in ozone may have occurred in some regions of the atmosphere prior to 1980, 1960 return dates are also reported. The projections reported on in this chapter are taken from a recent compilation of CCM simulations. The ozone projections, which also form the basis for the UV projections, are limited in their representativeness of possible futures since they mostly come from CCM simulations based on a single GHG emissions scenario (scenario A1B of Emissions Scenarios. A Special Report of Working Group III of the Intergovernmental Panel on Climate Change, Cambridge University Press, 2000) and a single ODS emissions scenario (adjusted A1 of the previous (2006) Ozone Assessment). Throughout this century, the vertical, latitudinal, and seasonal structure of the ozone distribution will be different from what it was in 1980. For this reason, ozone changes in different regions of the atmosphere are considered separately. • The projections of changes in ozone and surface clear-sky UV are broadly consistent with those reported on in the 2006 Assessment. • The capability of making projections and attribution of future ozone changes has been improved since the 2006 Assessment. Use of CCM simulations from an increased number of models extending through the entire period of ozone depletion and recovery from ODSs (1960–2100) as well as sensitivity simulations have allowed more robust projections of long-term changes in the stratosphere and of the relative contributions of ODSs and GHGs to those changes. • Global annually averaged total column ozone is projected to return to 1980 levels before the middle of the century and earlier than when stratospheric halogen loading returns to 1980 levels. CCM projections suggest that this early return is primarily a result of GHG-induced cooling of the upper stratosphere because the effects of circulation changes on tropical and extratropical ozone largely cancel. Global (90°S–90°N) annually averaged total column ozone will likely return to 1980 levels between 2025 and 2040, well before the return of stratospheric halogens to 1980 levels between 2045 and 2060. • Simulated changes in tropical total column ozone from 1960 to 2100 are generally small. The evolution of tropical total column ozone in models depends on the balance between upper stratospheric increases and lower stratospheric decreases. The upper stratospheric increases result from declining ODSs and a slowing of ozone destruction resulting from GHG-induced cooling. Ozone decreases in the lower stratosphere mainly result from an increase in tropical upwelling. From 1960 until around 2000, a general decline is simulated, followed by a gradual increase to values typical of 1980 by midcentury. Thereafter, although total column ozone amounts decline slightly again toward the end of the century, by 2080 they are no longer expected to be affected by ODSs. Confidence in tropical ozone projections is compromised by the fact that simulated decreases in column ozone to date are not supported by observations, suggesting that significant uncertainties remain. • Midlatitude total column ozone is simulated to evolve differently in the two hemispheres. Over northern midlatitudes, annually averaged total column ozone is projected to return to 1980 values between 2015 and 2030, while for southern midlatitudes the return to 1980 values is projected to occur between 2030 and 2040. The more rapid return to 1980 values in northern midlatitudes is linked to a more pronounced strengthening of the poleward transport of ozone due to the effects of increased GHG levels, and effects of Antarctic ozone depletion on southern midlatitudes. By 2100, midlatitude total column ozone is projected to be above 1980 values in both hemispheres. • October-mean Antarctic total column ozone is projected to return to 1980 levels after midcentury, later than in any other region, and yet earlier than when stratospheric halogen loading is projected to return to 1980 levels. The slightly earlier return of ozone to 1980 levels (2045–2060) results primarily from upper stratospheric cooling and resultant increases in ozone. The return of polar halogen loading to 1980 levels (2050–2070) in CCMs is earlier than in empirical models that exclude the effects of GHG-induced changes in circulation. Our confidence in the drivers of changes in Antarctic ozone is higher than for other regions because (i) ODSs exert a strong influence on Antarctic ozone, (ii) the effects of changes in GHG abundances are comparatively small, and (iii) projections of ODS emissions are more certain than those for GHGs. Small Antarctic ozone holes (areas of ozone <220 Dobson units, DU) could persist to the end of the 21st century. • March-mean Arctic total column ozone is projected to return to 1980 levels two to three decades before polar halogen loading returns to 1980 levels, and to exceed 1980 levels thereafter. While CCM simulations project a return to 1980 levels between 2020 and 2035, most models tend not to capture observed low temperatures and thus underestimate present-day Arctic ozone loss such that it is possible that this return date is biased early. Since the strengthening of the Brewer-Dobson circulation through the 21st century leads to increases in springtime Arctic column ozone, by 2100 Arctic ozone is projected to lie well above 1960 levels. Uncertainties in Projections • Conclusions dependent on future GHG levels are less certain than those dependent on future ODS levels since ODS emissions are controlled by the Montreal Protocol. For the six GHG scenarios considered by a few CCMs, the simulated differences in stratospheric column ozone over the second half of the 21st century are largest in the northern midlatitudes and the Arctic, with maximum differences of 20–40 DU between the six scenarios in 2100. • There remain sources of uncertainty in the CCM simulations. These include the use of prescribed ODS mixing ratios instead of emission fluxes as lower boundary conditions, the range of sea surface temperatures and sea ice concentrations, missing tropospheric chemistry, model parameterizations, and model climate sensitivity. • Geoengineering schemes for mitigating climate change by continuous injections of sulfur-containing compounds into the stratosphere, if implemented, would substantially affect stratospheric ozone, particularly in polar regions. Ozone losses observed following large volcanic eruptions support this prediction. However, sporadic volcanic eruptions provide limited analogs to the effects of continuous sulfur emissions. Preliminary model simulations reveal large uncertainties in assessing the effects of continuous sulfur injections. Expected Future Changes in Surface UV. While a number of factors, in addition to ozone, affect surface UV irradiance, the focus in this chapter is on the effects of changes in stratospheric ozone on surface UV. For this reason, clear-sky surface UV irradiance is calculated from ozone projections from CCMs. • Projected increases in midlatitude ozone abundances during the 21st century, in the absence of changes in other factors, in particular clouds, tropospheric aerosols, and air pollutants, will result in decreases in surface UV irradiance. Clear-sky erythemal irradiance is projected to return to 1980 levels on average in 2025 for the northern midlatitudes, and in 2035 for the southern midlatitudes, and to fall well below 1980 values by the second half of the century. However, actual changes in surface UV will be affected by a number of factors other than ozone. • In the absence of changes in other factors, changes in tropical surface UV will be small because changes in tropical total column ozone are projected to be small. By the middle of the 21st century, the model projections suggest surface UV to be slightly higher than in the 1960s, very close to values in 1980, and slightly lower than in 2000. The projected decrease in tropical total column ozone through the latter half of the century will likely result in clear-sky surface UV remaining above 1960 levels. Average UV irradiance is already high in the tropics due to naturally occurring low total ozone columns and high solar elevations. • The magnitude of UV changes in the polar regions is larger than elsewhere because ozone changes in polar regions are larger. For the next decades, surface clear-sky UV irradiance, particularly in the Antarctic, will continue to be higher than in 1980. Future increases in ozone and decreases in clear-sky UV will occur at slower rates than those associated with the ozone decreases and UV increases that occurred before 2000. In Antarctica, surface clear-sky UV is projected to return to 1980 levels between 2040 and 2060, while in the Arctic this is projected to occur between 2020 and 2030. By 2100, October surface clear-sky erythemal irradiance in Antarctica is likely to be between 5% below to 25% above 1960 levels, with considerable uncertainty. This is consistent with multi-model-mean October Antarctic total column ozone not returning to 1960 levels by 2100. In contrast, by 2100, surface clear-sky UV in the Arctic is projected to be 0–10% below 1960 levels.

Spectral characterisation of cooled filters for remote sensing

The improvements obtained on cooling atmospheric remote-sensing instruments for space flight applications has promoted research in characterization of the necessary optical filters. By modelling the effects of temperature on the dispersive spectrum of some constituent thin film materials, the cooled performance can be simulated and compared. multilayer filter designs with the measured spectra from actual filters. Two actual filters are discussed, for the 7µm region, one a composite cut-on/cut-off design of 13% HBW and the other an integral narrowband design of 4% HBW.

Performance model for cooled IR filters

A low-temperature model is described for infrared multilayer filters containing PbTe (or other semiconductor) and ZnSe (or other II/VI). The model is based on dielectric dispersion with semiconductor carrier dispersion added. It predicts an improved performance on cooling such as would be useful to avoid erroneous signals from optics in spaceflight radiometers. Agreement with measurement is obtained over the initial temperature range 70-400K and wavelength range 2.5-20µm.

Measurement of infrared multilayer filters at temperatures down to 4 deg K

Measurement is reported at 4 deg K (and blocked transmission below 10-5) of PbTe/ZnS thin-film filters deposited on Ge substrates. The reduced carrier-absorption which is obtained by cooling these PbTe films is found to accord with simple theory. Advantage for various high-performance multilayers by cooling is significant at the longer wavelengths, and has been verified.

Long-term effect of volcanic forcing on ocean heat content

Explosive volcanic eruptions cause episodic negative radiative forcing of the climate system. Using coupled atmosphere-ocean general circulation models (AOGCMs) subjected to historical forcing since the late nineteenth century, previous authors have shown that each large volcanic eruption is associated with a sudden drop in ocean heat content and sea-level from which the subsequent recovery is slow. Here we show that this effect may be an artefact of experimental design, caused by the AOGCMs not having been spun up to a steady state with volcanic forcing before the historical integrations begin. Because volcanic forcing has a long-term negative average, a cooling tendency is thus imposed on the ocean in the historical simulation. We recommend that an extra experiment be carried out in parallel to the historical simulation, with constant time-mean historical volcanic forcing, in order to correct for this effect and avoid misinterpretation of ocean heat content changes

Simulations of the London urban heat island

We present simulations of London's meteorology using the Met Office Unified Model with a new, sophisticated surface energy-balance scheme to represent the urban surfaces, called MORUSES. Simulations are performed with the urban surfaces represented and with the urban surfaces replaced with grass in order to calculate the urban increment on the local meteorology. The local urban effects were moderated to some extent by the passage of an onshore flow that propagated up the Thames estuary and across the city, cooling London slightly in the afternoon. Validations of screen-level temperature show encouraging agreement to within 1–2 K, when the urban increment is up to 5 K. The model results are then used to examine factors shaping the spatial and temporal structure of London's atmospheric boundary layer. The simulations reconcile the differences in the temporal evolution of the urban heat island (UHI) shown in various studies and demonstrate that the variation of UHI with time depends strongly on the urban fetch. The UHI at a location downwind of the city centre shows a decrease in UHI during the night, while the UHI at the city centre stays constant. Finally, the UHI at a location upwind of the city centre increases continuously. The magnitude of the UHI by the time of the evening transition increases with urban fetch. The urban increments are largest at night, when the boundary layer is shallow. The boundary layer experiences continued warming after sunset, as the heat from the urban fabric is released, and a weakly convective boundary layer develops across the city. The urban land-use fraction is the dominant control on the spatial structure in the sensible heat flux and the resulting urban increment, although even the weak advection present in this case study is sufficient to advect the peak temperature increments downwind of the most built-up areas. Copyright © 2011 Royal Meteorological Society and British Crown Copyright, the Met Office

Stratospheric temperature trends: our evolving understanding

We review the scientific literature since the 1960s to examine the evolution of modeling tools and observations that have advanced understanding of global stratospheric temperature changes. Observations show overall cooling of the stratosphere during the period for which they are available (since the late 1950s and late 1970s from radiosondes and satellites, respectively), interrupted by episodes of warming associated with volcanic eruptions, and superimposed on variations associated with the solar cycle. There has been little global mean temperature change since about 1995. The temporal and vertical structure of these variations are reasonably well explained bymodels that include changes in greenhouse gases, ozone, volcanic aerosols, and solar output, although there are significant uncertainties in the temperature observations and regarding the nature and influence of past changes in stratospheric water vapor. As a companion to a recent WIREs review of tropospheric temperature trends, this article identifies areas of commonality and contrast between the tropospheric and stratospheric trend literature. For example, the increased attention over time to radiosonde and satellite data quality has contributed to better characterization of uncertainty in observed trends both in the troposphere and in the lower stratosphere, and has highlighted the relative deficiency of attention to observations in the middle and upper stratosphere. In contrast to the relatively unchanging expectations of surface and tropospheric warming primarily induced by greenhouse gas increases, stratospheric temperature change expectations have arisen from experiments with a wider variety of model types, showingmore complex trend patterns associated with a greater diversity of forcing agents.

The temperature response to stratospheric water vapour changes

This study uses an analytical model, based on the cooling-to-space approximation, and a fixed dynamical heating model to investigate the structure of the stratospheric cooling that occurs in response to a uniform increase in stratospheric water vapour (SWV). At all latitudes, the largest cooling occurs in the lower stratosphere and decreases in magnitude with height. The cooling is strongly enhanced in the Extratropics compared to the Tropics. This is markedly different to the case of an increase in CO2, which causes maximum cooling near the stratopause and a small warming in the tropical lower stratosphere. The qualitative differences in the structure of the cooling can be explained by the smaller opacity of water vapour bands in the stratosphere compared to CO2. The small opacity means that the magnitude of the initial heating rate perturbation only decreases by a factor of four between the upper and lower stratosphere for a SWV perturbation. Therefore, to balance the heating rate perturbation, the largest temperature change is required in the lower stratosphere. Increasing the background concentration of SWV causes the water vapour bands to become more opaque. For a SWV perturbation applied to a background SWV concentration ≥30 ppmv, the heating rate perturbation and temperature change structurally resemble those from an increase in CO2. In the Extratropics, the lower height of the tropopause is found to cause the enhancement in the cooling at those latitudes. By controlling the depth of atmosphere which adjusts to the SWV perturbation, the tropopause height affects the net exchange of radiation between the layers in the stratosphere as they cool. The latitudinal gradient in upwelling infrared radiation at the tropopause and variations in the background temperature are found to have only a minor effect on the structure of the stratospheric temperature response to a change in SWV.

Combining satellite data and models to estimate cloud radiative effect at the surface and in the atmosphere

Satellite measurements and numerical forecast model reanalysis data are used to compute an updated estimate of the cloud radiative effect on the global multi-annual mean radiative energy budget of the atmosphere and surface. The cloud radiative cooling effect through reflection of shortwave radiation dominates over the longwave heating effect, resulting in a net cooling of the climate system of –21 Wm-2. The shortwave radiative effect of cloud is primarily manifest as a reduction in the solar radiation absorbed at the surface of -53 Wm-2. Clouds impact longwave radiation by heating the moist tropical atmosphere (up to around 40 Wm-2 for global annual means) while enhancing the radiative cooling of the atmosphere over other regions, in particular higher latitudes and sub-tropical marine stratocumulus regimes. While clouds act to cool the climate system during the daytime, the cloud greenhouse effect heats the climate system at night. The influence of cloud radiative effect on determining cloud feedbacks and changes in the water cycle are discussed.

Stratospheric heating by potential geoengineering aerosols

A fixed dynamical heating model is used to investigate the pattern of zonal-mean stratospheric temperature change resulting from geoengineering with aerosols composed of sulfate, titania, limestone and soot. Aerosol always heats the tropical lower stratosphere, but at the poles the response can be either heating, cooling, or neutral. The sign of the change in stratospheric Pole-Equator temperature difference depends on aerosol type, size and season. This has implications for modelling geoengineering impacts and the response of the stratospheric circulation.

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.

The last glacial cycle: transient simulations with an AOGCM

A number of transient climate runs simulating the last 120kyr have been carried out using FAMOUS, a fast atmosphere-ocean general circulation model (AOGCM). This is the first time such experiments have been done with a full AOGCM, providing a three-dimensional simulation of both atmosphere and ocean over this period. Our simulation thus includes internally generated temporal variability over periods from days to millennia, and physical, detailed representations of important processes such as clouds and precipitation. Although the model is fast, computational restrictions mean that the rate of change of the forcings has been increased by a factor of 10, making each experiment 12kyr long. Atmospheric greenhouse gases (GHGs), northern hemisphere ice sheets and variations in solar radiation arising from changes in the Earth's orbit are treated as forcing factors, and are applied either separately or combined in different experiments. The long-term temperature changes on Antarctica match well with reconstructions derived from ice-core data, as does variability on timescales longer than 10 kyr. Last Glacial Maximum (LGM) cooling on Greenland is reasonably well simulated, although our simulations, which lack ice-sheet meltwater forcing, do not reproduce the abrupt, millennial scale climate shifts seen in northern hemisphere climate proxies or their slower southern hemisphere counterparts. The spatial pattern of sea surface cooling at the LGM matches proxy reconstructions reasonably well. There is significant anti-correlated variability in the strengths of the Atlantic Meridional Overturning Circulation (AMOC) and the Antarctic Circumpolar Current (ACC) on timescales greater than 10kyr in our experiments. We find that GHG forcing weakens the AMOC and strengthens the ACC, whilst the presence of northern hemisphere ice-sheets strengthens the AMOC and weakens the ACC. The structure of the AMOC at the LGM is found to be sensitive to the details of the ice-sheet reconstruction used. The precessional component of the orbital forcing induces ~20kyr oscillations in the AMOC and ACC, whose amplitude is mediated by changes in the eccentricity of the Earth's orbit. These forcing influences combine, to first order, in a linear fashion to produce the mean climate and ocean variability seen in the run with all forcings.