The climate response to CO2 of the Hadley Centre coupled AOGCM with and without flux adjustment

Coupled atmosphere‐ocean general circulation models have a tendency to drift away from a realistic climatology. The modelled climate response to an increase of CO2 concentration may be incorrect if the simulation of the current climate has significant errors, so in many models, including ours, the drift is counteracted by applying prescribed fluxes of heat and fresh water at the ocean‐atmosphere interface in addition to the calculated surface exchanges. Since the additional fluxes do not have a physical basis, the use of this technique of “flux adjustment” itself introduces some uncertainty in the simulated response to increased CO2. We find that the global‐average temperature response of our model to CO2 increasing at 1% per year is about 30% less without flux adjustment than with flux adjustment. The geographical patterns of the response are similar, indicating that flux adjustment is not causing any gross distortion. The reduced size of the response is due to more effective vertical transport of heat into the ocean, and a somewhat smaller climate sensitivity. Although the response in both cases lies within the generally accepted range for the climate sensitivity, systematic uncertainties of this size are clearly undesirable, and the best strategy for future development is to improve the climate model in order to reduce the need for flux adjustment.

A simple relationship between volcanic sulfate aerosol optical depth and surface temperature change simulated in an atmosphere-ocean general circulation model.

In this study we quantify the relationship between the aerosol optical depth increase from a volcanic eruption and the severity of the subsequent surface temperature decrease. This investigation is made by simulating 10 different sizes of eruption in a global circulation model (GCM) by changing stratospheric sulfate aerosol optical depth at each time step. The sizes of the simulated eruptions range from Pinatubo‐sized up to the magnitude of supervolcanic eruptions around 100 times the size of Pinatubo. From these simulations we find that there is a smooth monotonic relationship between the global mean maximum aerosol optical depth anomaly and the global mean temperature anomaly and we derive a simple mathematical expression which fits this relationship well. We also construct similar relationships between global mean aerosol optical depth and the temperature anomaly at every individual model grid box to produce global maps of best‐fit coefficients and fit residuals. These maps are used with caution to find the eruption size at which a local temperature anomaly is clearly distinct from the local natural variability and to approximate the temperature anomalies which the model may simulate following a Tambora‐sized eruption. To our knowledge, this is the first study which quantifies the relationship between aerosol optical depth and resulting temperature anomalies in a simple way, using the wealth of data that is available from GCM simulations.

The distribution of solar wind speeds during solar minimum: calibration for numerical solar wind modeling constraints on the source of the slow solar wind

It took the solar polar passage of Ulysses in the early 1990s to establish the global structure of the solar wind speed during solar minimum. However, it remains unclear if the solar wind is composed of two distinct populations of solar wind from different sources (e.g., closed loops which open up to produce the slow solar wind) or if the fast and slow solar wind rely on the superradial expansion of the magnetic field to account for the observed solar wind speed variation. We investigate the solar wind in the inner corona using the Wang-Sheeley-Arge (WSA) coronal model incorporating a new empirical magnetic topology–velocity relationship calibrated for use at 0.1 AU. In this study the empirical solar wind speed relationship was determined by using Helios perihelion observations, along with results from Riley et al. (2003) and Schwadron et al. (2005) as constraints. The new relationship was tested by using it to drive the ENLIL 3-D MHD solar wind model and obtain solar wind parameters at Earth (1.0 AU) and Ulysses (1.4 AU). The improvements in speed, its variability, and the occurrence of high-speed enhancements provide confidence that the new velocity relationship better determines the solar wind speed in the outer corona (0.1 AU). An analysis of this improved velocity field within the WSA model suggests the existence of two distinct mechanisms of the solar wind generation, one for fast and one for slow solar wind, implying that a combination of present theories may be necessary to explain solar wind observations.

Quantifying the response of the ORAC aerosol optical depth retrieval for MSG SEVIRI to aerosol model assumptions

We test the response of the Oxford-RAL Aerosol and Cloud (ORAC) retrieval algorithm for MSG SEVIRI to changes in the aerosol properties used in the dust aerosol model, using data from the Dust Outflow and Deposition to the Ocean (DODO) flight campaign in August 2006. We find that using the observed DODO free tropospheric aerosol size distribution and refractive index increases simulated top of the atmosphere radiance at 0.55 µm assuming a fixed erosol optical depth of 0.5 by 10–15 %, reaching a maximum difference at low solar zenith angles. We test the sensitivity of the retrieval to the vertical distribution f the aerosol and find that this is unimportant in determining simulated radiance at 0.55 µm. We also test the ability of the ORAC retrieval when used to produce the GlobAerosol dataset to correctly identify continental aerosol outflow from the African continent and we find that it poorly constrains aerosol speciation. We develop spatially and temporally resolved prior distributions of aerosols to inform the retrieval which incorporates five aerosol models: desert dust, maritime, biomass burning, urban and continental. We use a Saharan Dust Index and the GEOS-Chem chemistry transport model to describe dust and biomass burning aerosol outflow, and compare AOD using our speciation against the GlobAerosol retrieval during January and July 2006. We find AOD discrepancies of 0.2–1 over regions of intense biomass burning outflow, where AOD from our aerosol speciation and GlobAerosol speciation can differ by as much as 50 - 70 %.

Generation of inertia-gravity waves in the rotating thermal annulus by a localised boundary layer instability

Waves with periods shorter than the inertial period exist in the atmosphere (as inertia-gravity waves) and in the oceans (as Poincaré and internal gravity waves). Such waves owe their origin to various mechanisms, but of particular interest are those arising either from local secondary instabilities or spontaneous emission due to loss of balance. These phenomena have been studied in the laboratory, both in the mechanically-forced and the thermally-forced rotating annulus. Their generation mechanisms, especially in the latter system, have not yet been fully understood, however. Here we examine short period waves in a numerical model of the rotating thermal annulus, and show how the results are consistent with those from earlier laboratory experiments. We then show how these waves are consistent with being inertia-gravity waves generated by a localised instability within the thermal boundary layer, the location of which is determined by regions of strong shear and downwelling at certain points within a large-scale baroclinic wave flow. The resulting instability launches small-scale inertia-gravity waves into the geostrophic interior of the flow. Their behaviour is captured in fully nonlinear numerical simulations in a finite-difference, 3D Boussinesq Navier-Stokes model. Such a mechanism has many similarities with those responsible for launching small- and meso-scale inertia-gravity waves in the atmosphere from fronts and local convection.

Centennial changes in the heliospheric magnetic field and open solar flux: the consensus view from geomagnetic data and cosmogenic isotopes and its implications

Svalgaard and Cliver (2010) recently reported a consensus between the various reconstructions of the heliospheric field over recent centuries. This is a significant development because, individually, each has uncertainties introduced by instrument calibration drifts, limited numbers of observatories, and the strength of the correlations employed. However, taken collectively, a consistent picture is emerging. We here show that this consensus extends to more data sets and methods than reported by Svalgaard and Cliver, including that used by Lockwood et al. (1999), when their algorithm is used to predict the heliospheric field rather than the open solar flux. One area where there is still some debate relates to the existence and meaning of a floor value to the heliospheric field. From cosmogenic isotope abundances, Steinhilber et al. (2010) have recently deduced that the near-Earth IMF at the end of the Maunder minimum was 1.80 ± 0.59 nT which is considerably lower than the revised floor of 4nT proposed by Svalgaard and Cliver. We here combine cosmogenic and geomagnetic reconstructions and modern observations (with allowance for the effect of solar wind speed and structure on the near-Earth data) to derive an estimate for the open solar flux of (0.48 ± 0.29) × 1014 Wb at the end of the Maunder minimum. By way of comparison, the largest and smallest annual means recorded by instruments in space between 1965 and 2010 are 5.75 × 1014 Wb and 1.37 × 1014 Wb, respectively, set in 1982 and 2009, and the maximum of the 11 year running means was 4.38 × 1014 Wb in 1986. Hence the average open solar flux during the Maunder minimum is found to have been 11% of its peak value during the recent grand solar maximum.

Remote forcing of the Antarctic Circumpolar Current by diapycnal mixing

We show that diapycnal mixing can drive a significant Antarctic Circumpolar Current (ACC) volume transport, even when the mixing is located remotely in northern-hemisphere ocean basins. In the case of remote forcing, the globally-averaged diapycnal mixing coefficient is the important parameter. This result is anticipated from theoretical arguments and demonstrated in a global ocean circulation model. The impact of enhanced diapycnal mixing on the ACC during glacial periods is discussed.

Multimodel climate and variability of the stratosphere

The stratospheric climate and variability from simulations of sixteen chemistry‐climate models is evaluated. On average the polar night jet is well reproduced though its variability is less well reproduced with a large spread between models. Polar temperature biases are less than 5 K except in the Southern Hemisphere (SH) lower stratosphere in spring. The accumulated area of low temperatures responsible for polar stratospheric cloud formation is accurately reproduced for the Antarctic but underestimated for the Arctic. The shape and position of the polar vortex is well simulated, as is the tropical upwelling in the lower stratosphere. There is a wide model spread in the frequency of major sudden stratospheric warnings (SSWs), late biases in the breakup of the SH vortex, and a weak annual cycle in the zonal wind in the tropical upper stratosphere. Quantitatively, “metrics” indicate a wide spread in model performance for most diagnostics with systematic biases in many, and poorer performance in the SH than in the Northern Hemisphere (NH). Correlations were found in the SH between errors in the final warming, polar temperatures, the leading mode of variability, and jet strength, and in the NH between errors in polar temperatures, frequency of major SSWs, and jet strength. Models with a stronger QBO have stronger tropical upwelling and a colder NH vortex. Both the qualitative and quantitative analysis indicate a number of common and long‐standing model problems, particularly related to the simulation of the SH and stratospheric variability.

On the response of the Antarctic Circumpolar Current transport to climate change in coupled climate models

The CMIP3 (IPCC AR4) models show a consistent intensification and poleward shift of the westerly winds over the Southern Ocean during the 21st century. However, the responses of the Antarctic Circumpolar Currents (ACC) show great diversity in these models, with many even showing reductions in transport. To obtain some understanding of diverse responses in the ACC transport, we investigate both external atmospheric and internal oceanic processes that control the ACC transport responses in these models. While the strengthened westerlies act to increase the tilt of isopycnal surfaces and hence the ACC transport through Ekman pumping effects, the associated changes in buoyancy forcing generally tend to reduce the surface meridional density gradient. The steepening of isopycnal surfaces induced by increased wind forcing leads to enhanced (parameterized) eddy-induced transports that act to reduce the isopycnal slopes. There is also considerable narrowing of the ACC that tends to reduce the ACC transport, caused mainly by the poleward shifts of the subtropical gyres and to a lesser extent by the equatorward expansions of the subpolar gyres in some models. If the combined effect of these retarding processes is larger than that of enhanced Ekman pumping, the ACC transport will be reduced. In addition, the effect of Ekman pumping on the ACC is reduced in weakly stratified models. These findings give insight into the reliability of IPCC-class model predictions of the Southern Ocean circulation, and into the observed decadal-scale steady ACC transport.

How is open solar magnetic flux lost over the solar cycle?

The Sun's open magnetic field, magnetic flux dragged out into the heliosphere by the solar wind, varies by approximately a factor of 2 over the solar cycle. We consider the evolution of open solar flux in terms of a source and loss term. Open solar flux creation is likely to proceed at a rate dependent on the rate of photospheric flux emergence, which can be roughly parameterized by sunspot number or coronal mass ejection rate, when available. The open solar flux loss term is more difficult to relate to an observable parameter. The supersonic nature of the solar wind means open solar flux can only be removed by near-Sun magnetic reconnection between open solar magnetic field lines, be they open or closed heliospheric field lines. In this study we reconstruct open solar flux over the last three solar cycles and demonstrate that the loss term may be related to the degree to which the heliospheric current sheet (HCS) is warped, i.e., locally tilted from the solar rotation direction. This can account for both the large dip in open solar flux at the time of sunspot maximum as well as the asymmetry in open solar flux during the rising and declining phases of the solar cycle. The observed cycle-to-cycle variability is also well matched. Following Sheeley et al. (2001), we attribute modulation of open solar flux by the degree of warp of the HCS to the rate at which opposite polarity open solar flux is brought together by differential rotation.

Bistability of the Atlantic overturning circulation in a global climate model and links to ocean freshwater transport

The possibility of a rapid collapse in the strength of the Atlantic meridional overturning circulation (AMOC), with associated impacts on climate, has long been recognized. The suggested basis for this risk is the existence of two stable regimes of the AMOC (‘on’ and ‘off’), and such bistable behaviour has been identified in a range of simplified climate models. However, up to now, no state-of-the-art atmosphere-ocean coupled global climate model (AOGCM) has exhibited such behaviour, leading to the interpretation that the AMOC is more stable than simpler models indicate. Here we demonstrate AMOC bistability in the response to freshwater perturbations in the FAMOUS AOGCM - the most complex AOGCM to exhibit such behaviour to date. The results also support recent suggestions that the direction of the net freshwater transport at the southern boundary of the Atlantic by the AMOC may be a useful physical indicator of the existence of bistability. We also present new estimates for this net freshwater transport by the AMOC from a range of ocean reanalyses which suggest that the Atlantic AMOC is currently in a bistable regime, although with large uncertainties. More accurate observational constraints, and an improved physical understanding of this quantity, could help narrow uncertainty in the future evolution of the AMOC and to assess the risk of a rapid AMOC collapse.