The relationship between the North American summer monsoon, the Rocky Mountains and the North Pacific subtropical anticyclone in HadAM3

In this study the relationship between the North American monsoon, the Californian sea surface temperature (SST) cold pool, the Rocky Mountains and the North Pacific subtropical anticyclone is investigated using the Hadley Centre's atmospheric climate model, HadAM3. In 1996 Hoskins hypothesized that heating in the North American monsoon might be important for the maintenance of the summertime North Pacific subtropical anticyclone, since the monsoon heating may induce descent to the north-west of the monsoon in the descending eastern flank of the subtropical anticyclone. This descent is further enhanced by radiative cooling and is associated with equatorward surface winds parallel to the western coast of North America. These equatorward winds induce oceanic upwelling of cold water and contribute to the formation of the Californian SST cold pool, which may feed back on the anticyclone by further suppressing convection and inducing descent. More recently, Rodwell and Hoskins also investigated the global summer monsoon–subtropical anticyclone relationship. They examined the role that mountains play in impeding the progress of the low-level mid-latitude westerlies, either deflecting the westerlies northwards where they ascend along the sloping mid-latitude isentropes or deflecting them southwards forcing them to descend along the isentropes. In particular, the introduction of the Rockies into a primitive-equation model adiabatically induces descent in the eastern descending flank of the North Pacific subtropical anticyclone. These hypothesized mechanisms have been investigated using HadAM3, focusing on the possible suppression of convection by the Californian SST cold pool, the response of the North Pacific subtropical anticyclone to the strength of the North American monsoon and the ‘blocking’ of the mid-latitude westerlies by the Rocky Mountains. The role of the Rockies is examined by integrating the model with modified orography for the Rocky Mountains. Changing the height of the Rockies alters the circulation in a way consistent with the mechanism outlined above. Higher Rocky mountains force the westerlies southwards, inducing descent in the eastern flank of the subtropical anticyclone as the air descends along the sloping isentropes. The relationship between the North American monsoon and the North Pacific subtropical anticyclone is investigated by suppressing the monsoon in HadAM3. The suppression of the monsoon is accomplished by increasing the surface albedo over Mexico, which induces anomalous ascent on the eastward flank of the subtropical anticyclone and anomalous polewards surface winds along the western coast of the North American continent, also providing support for the above hypothesis. The removal of the Californian SST cold pool, however, has a statistically insignificant effect on the model, suggesting that in this model the feedback of the SST cold pool on the eastern flank of the anticyclone is weak.

Upper-ocean heat budget and ocean eddy transport in the south-east Pacific in a high-resolution coupled model

We present an analysis of the oceanic heat advection and its variability in the upper 500 m in the southeastern tropical Pacific (100W–75W, 25S–10S) as simulated by the global coupled model HiGEM, which has one of the highest resolutions currently used in long-term integrations. The simulated climatology represents a temperature advection field arising from transient small-scale (<450 km) features, with structures and transport that appear consistent with estimates based on available observational data for the mooring at 20S, 85W. The transient structures are very persistent (>4 months), and in specific locations they generate an important contribution to the local upper-ocean heat budget, characterised by scales of a few hundred kilometres, and periods of over a year. The contribution from such structures to the local, long-term oceanic heat budget however can be of either sign, or vanishing, depending on the location; and, although there appears some organisation in preferential areas of activity, the average over the entire region is small. While several different mechanisms may be responsible for the temperature advection by transients, we find that a significant, and possibly dominant, component is associated with vortices embedded in the large-scale, climatological salinity gradient associated with the fresh intrusion of mid-latitude intermediate water which penetrates north-westward beneath the tropical thermocline

The North Atlantic jet stream under climate change and its relation to the NAO and EA patterns

This paper describes recent variations of the North Atlantic eddy-driven jet stream and analyzes the mean response of the jet to anthropogenic forcing in climate models. Jet stream changes are analyzed both using a direct measure of the near-surface westerly wind maximum and using an EOF-based approach. This allows jet stream changes to be related to the widely used leading patterns of variability: the North Atlantic Oscillation (NAO) and East Atlantic (EA) pattern. Viewed in NAO–EA state space, isolines of jet latitude and speed resemble a distorted polar coordinate system, highlighting the dependence of the jet stream quantities on both spatial patterns. Some differences in the results of the two methods are discussed, but both approaches agree on the general characteristics of the climate models. While there is some agreement between models on a poleward shift of the jet stream in response to anthropogenic forcing, there is still considerable spread between different model projections, especially in winter. Furthermore, the model responses to forcing are often weaker than their biases when compared to a reanalysis. Diagnoses of jet stream changes can be sensitive to the methodologies used, and several aspects of this are also discussed.

A mechanism for the effect of tropospheric jet structure on the annular mode-like response to stratospheric forcing

For many climate forcings the dominant response of the extratropical circulation is a latitudinal shift of the tropospheric mid-latitude jets. The magnitude of this response appears to depend on climatological jet latitude in general circulation models (GCMs): lower latitude jets exhibit a larger shift. The reason for this latitude dependence is investigated for a particular forcing, heating of the equatorial stratosphere, which shifts the jet poleward. Spin-up ensembles with a simplified GCM are used to examine the evolution of the response for five different jet structures. These differ in the latitude of the eddy-driven jet, but have similar sub-tropical zonal winds. It is found that lower latitude jets exhibit a larger response due to stronger tropospheric eddy-mean flow feedbacks. A dominant feedback responsible for enhancing the poleward shift is an enhanced equatorward refraction of the eddies, resulting in an increased momentum flux, poleward of the low latitude critical line. The sensitivity of feedback strength to jet structure is associated with differences in the coherence of this behaviour across the spectrum of eddy phase speeds. In the configurations used, the higher latitude jets have a wider range of critical latitude locations. This reduces the coherence of the momentum flux anomalies associated with different phase speeds, with low phase speeds opposing the effect of high phase speeds. This suggests that, for a given sub-tropical zonal wind strength, the latitude of the eddy driven jet affects the feedback through its influence on the width of the region of westerly winds and the range of critical latitudes on the equatorward flank of the jet.

The effect of nonlinearity on the variational assimilation of satellite observations using a simple column model

Cloud imagery is not currently used in numerical weather prediction (NWP) to extract the type of dynamical information that experienced forecasters have extracted subjectively for many years. For example, rapidly developing mid-latitude cyclones have characteristic signatures in the cloud imagery that are most fully appreciated from a sequence of images rather than from a single image. The Met Office is currently developing a technique to extract dynamical development information from satellite imagery using their full incremental 4D-Var (four-dimensional variational data assimilation) system. We investigate a simplified form of this technique in a fully nonlinear framework. We convert information on the vertical wind field, w(z), and profiles of temperature, T(z, t), and total water content, qt (z, t), as functions of height, z, and time, t, to a single brightness temperature by defining a 2D (vertical and time) variational assimilation testbed. The profiles of w, T and qt are updated using a simple vertical advection scheme. We define a basic cloud scheme to obtain the fractional cloud amount and, when combined with the temperature field, we convert this information into a brightness temperature, having developed a simple radiative transfer scheme. With the exception of some matrix inversion routines, all our code is developed from scratch. Throughout the development process we test all aspects of our 2D assimilation system, and then run identical twin experiments to try and recover information on the vertical velocity, from a sequence of observations of brightness temperature. This thesis contains a comprehensive description of our nonlinear models and assimilation system, and the first experimental results.

The FAMOUS climate model (versions XFXWB and XFHCC): description update to version XDBUA

FAMOUS is an ocean-atmosphere general circulation model of low resolution, based on version 4.5 of the UK MetOffice Unified Model. Here we update the model description to account for changes in the model as it is used in the CMIP5 EMIC model intercomparison project (EMICmip) and a number of other studies. Most of these changes correct errors found in the code. The EMICmip version of the model (XFXWB) has a better-conserved water budget and additional cooling in some high latitude areas, but otherwise has a similar climatology to previous versions of FAMOUS. A variant of XFXWB is also described, with changes to the dynamics at the top of the model which improve the model climatology (XFHCC).

Antarctic air over New Zealand following vortex breakdown in 1998

An ozonesonde profile over the Network for Detection of Stratospheric Change (NDSC) site at Lauder (45.0° S, 169.7° E), New Zealand, for 24 December 1998 showed atypically low ozone centered around 24 km altitude (600 K potential temperature). The origin of the anomaly is explained using reverse domain filling (RDF) calculations combined with a PV/O3 fitting technique applied to ozone measurements from the Polar Ozone and Aerosol Measurement (POAM) III instrument. The RDF calculations for two isentropic surfaces, 550 and 600 K, show that ozone-poor air from the Antarctic polar vortex reached New Zealand on 24–26 December 1998. The vortex air on the 550 K isentrope originated in the ozone hole region, unlike the air on 600 K where low ozone values were caused by dynamical effects. High-resolution ozone maps were generated, and their examination shows that a vortex remnant situated above New Zealand was the cause of the altered ozone profile on 24 December. The maps also illustrate mixing of the vortex filaments into southern midlatitudes, whereby the overall mid-latitude ozone levels were decreased.

Response of the North Atlantic storm track to climate change shaped by ocean–atmosphere coupling

A poleward shift of the mid-latitude storm tracks in response to anthropogenic greenhouse-gas forcing has been diagnosed in climate model simulations1, 2. Explanations of this effect have focused on atmospheric dynamics3, 4, 5, 6, 7. However, in contrast to storm tracks in other regions, the North Atlantic storm track responds by strengthening and extending farther east, in particular on its southern flank8. These adjustments are associated with an intensification and extension of the eddy-driven jet towards western Europe9 and are expected to have considerable societal impacts related to a rise in storminess in Europe10, 11, 12. Here, we apply a regression analysis to an ensemble of coupled climate model simulations to show that the coupling between ocean and atmosphere shapes the distinct storm-track response to greenhouse-gas forcing in the North Atlantic region. In the ensemble of simulations we analyse, at least half of the differences between the storm-track responses of different models are associated with uncertainties in ocean circulation changes. We compare the fully coupled simulations with both the associated slab model simulations and an ocean-forced experiment with one climate model to establish causality. We conclude that uncertainties in the response of the North Atlantic storm track to anthropogenic emissions could be reduced through tighter constraints on the future ocean circulation.

Controlling the computational modes of the arbitrarily structured C-grid

The arbitrarily structured C-grid, TRiSK (Thuburn, Ringler, Skamarock and Klemp, 2009, 2010) is being used in the ``Model for Prediction Across Scales'' (MPAS) and is being considered by the UK Met Office for their next dynamical core. However the hexagonal C-grid supports a branch of spurious Rossby modes which lead to erroneous grid-scale oscillations of potential vorticity (PV). It is shown how these modes can be harmlessly controlled by using upwind-biased interpolation schemes for PV. A number of existing advection schemes for PV are tested, including that used in MPAS, and none are found to give adequate results for all grids and all cases. Therefore a new scheme is proposed; continuous, linear-upwind stabilised transport (CLUST), a blend between centred and linear-upwind with the blend dependent on the flow direction with respect to the cell edge. A diagnostic of grid-scale oscillations is proposed which gives further discrimination between schemes than using potential enstrophy alone and indeed some schemes are found to destroy potential enstrophy while grid-scale oscillations grow. CLUST performs well on hexagonal-icosahedral grids and unrotated skipped latitude-longitude grids of the sphere for various shallow water test cases. Despite the computational modes, the hexagonal icosahedral grid performs well since these modes are easy and harmless to filter. As a result TRiSK appears to perform better than a spectral shallow water model.

Simulated anthropogenic changes in the Brewer–Dobson circulation, including its extension to high latitudes

Recent studies using comprehensive middle atmosphere models predict a strengthening of the Brewer-Dobson circulation in response to climate change. To gain confidence in the realism of this result it is important to quantify and understand the contributions from the different components of stratospheric wave drag that cause this increase. Such an analysis is performed here using three 150-yr transient simulations from the Canadian Middle Atmosphere Model (CMAM), a Chemistry-Climate Model that simulates climate change and ozone depletion and recovery. Resolved wave drag and parameterized orographic gravity wave drag account for 60% and 40%, respectively, of the long-term trend in annual mean net upward mass flux at 70 hPa, with planetary waves accounting for 60% of the resolved wave drag trend. Synoptic wave drag has the strongest impact in northern winter, where it accounts for nearly as much of the upward mass flux trend as planetary wave drag. Owing to differences in the latitudinal structure of the wave drag changes, the relative contribution of resolved and parameterized wave drag to the tropical upward mass flux trend over any particular latitude range is highly sensitive to the range of latitudes considered. An examination of the spatial structure of the climate change response reveals no straightforward connection between the low-latitude and high-latitude changes: while the model results show an increase in Arctic downwelling in winter, they also show a decrease in Antarctic downwelling in spring. Both changes are attributed to changes in the flux of stationary planetary wave activity into the stratosphere.

Mechanisms linking volcanic aerosols to the Atlantic meridional overturning circulation

This study examines the sensitivity of the climate system to volcanic aerosol forcing in the third climate configuration of the Met Office Unified Model (HadCM3). The main test case was based on the 1880s when there were several volcanic eruptions, the well-known Krakatau being the largest. These eruptions increased atmospheric aerosol concentrations and induced a period of global cooling surface temperatures. In this study, an ensemble of HadCM3 has been integrated with the standard set of radiative forcings and aerosols from the Intergovernmental Panel on Climate Change Fourth Assessment Report simulations, from 1860 to present. A second ensemble removes the volcanic aerosols from 1880 to 1899. The all-forcings ensemble shows an attributable 1.2-Sv (1 Sv ≡ 106 m3 s−1) increase in the Atlantic meridional overturning circulation (AMOC) at 45°N—with a 0.04-PW increase in meridional heat transport at 40°N and increased northern Atlantic SSTs—starting around 1894, approximately 11 years after the first eruption, and lasting a further 10 years at least. The mechanisms responsible are traced to the Arctic, with suppression of the global water cycle (high-latitude precipitation), which leads to an increase in upper-level Arctic and Greenland Sea salinities. This then leads to increased convection in the Greenland–Iceland–Norwegian (GIN) Seas, enhanced Denmark Strait overflows, and AMOC changes with density anomalies traceable southward along the western Atlantic boundary. The authors investigate whether a similar response to the Pinatubo eruption in 1991 could still be ongoing, but do not find strong evidence.

The Aqua-Planet Experiment (APE): Response to changed meridional SST profile

This paper explores the sensitivity of Atmospheric General Circulation Model (AGCM) simulations to changes in the meridional distribution of sea surface temperature (SST). The simulations are for an aqua-planet, a water covered Earth with no land, orography or sea-ice and with specified zonally symmetric SST. Simulations from 14 AGCMs developed for Numerical Weather Prediction and climate applications are compared. Four experiments are performed to study the sensitivity to the meridional SST profile. These profiles range from one in which the SST gradient continues to the equator to one which is flat approaching the equator, all with the same maximum SST at the equator. The zonal mean circulation of all models shows strong sensitivity to latitudinal distribution of SST. The Hadley circulation weakens and shifts poleward as the SST profile flattens in the tropics. One question of interest is the formation of a double versus a single ITCZ. There is a large variation between models of the strength of the ITCZ and where in the SST experiment sequence they transition from a single to double ITCZ. The SST profiles are defined such that as the equatorial SST gradient flattens, the maximum gradient increases and moves poleward. This leads to a weakening of the mid-latitude jet accompanied by a poleward shift of the jet core. Also considered are tropical wave activity and tropical precipitation frequency distributions. The details of each vary greatly between models, both with a given SST and in the response to the change in SST. One additional experiment is included to examine the sensitivity to an off-equatorial SST maximum. The upward branch of the Hadley circulation follows the SST maximum off the equator. The models that form a single precipitation maximum when the maximum SST is on the equator shift the precipitation maximum off equator and keep it centered over the SST maximum. Those that form a double with minimum on the equatorial maximum SST shift the double structure off the equator, keeping the minimum over the maximum SST. In both situations only modest changes appear in the shifted profile of zonal average precipitation. When the upward branch of the Hadley circulation moves into the hemisphere with SST maximum, the zonal average zonal, meridional and vertical winds all indicate that the Hadley cell in the other hemisphere dominates.

The Aqua-Planet Experiment (APE): CONTROL SST Simulation

Climate simulations by 16 atmospheric general circulation models (AGCMs) are compared on an aqua-planet, a water-covered Earth with prescribed sea surface temperature varying only in latitude. The idealised configuration is designed to expose differences in the circulation simulated by different models. Basic features of the aqua-planet climate are characterised by comparison with Earth. The models display a wide range of behaviour. The balanced component of the tropospheric mean flow, and mid-latitude eddy covariances subject to budget constraints, vary relatively little among the models. In contrast, differences in damping in the dynamical core strongly influence transient eddy amplitudes. Historical uncertainty in modelled lower stratospheric temperatures persists in APE. Aspects of the circulation generated more directly by interactions between the resolved fluid dynamics and parameterized moist processes vary greatly. The tropical Hadley circulation forms either a single or double inter-tropical convergence zone (ITCZ) at the equator, with large variations in mean precipitation. The equatorial wave spectrum shows a wide range of precipitation intensity and propagation characteristics. Kelvin mode-like eastward propagation with remarkably constant phase speed dominates in most models. Westward propagation, less dispersive than the equatorial Rossby modes, dominates in a few models or occurs within an eastward propagating envelope in others. The mean structure of the ITCZ is related to precipitation variability, consistent with previous studies. The aqua-planet global energy balance is unknown but the models produce a surprisingly large range of top of atmosphere global net flux, dominated by differences in shortwave reflection by clouds. A number of newly developed models, not optimised for Earth climate, contribute to this. Possible reasons for differences in the optimised models are discussed. The aqua-planet configuration is intended as one component of an experimental hierarchy used to evaluate AGCMs. This comparison does suggest that the range of model behaviour could be better understood and reduced in conjunction with Earth climate simulations. Controlled experimentation is required to explore individual model behaviour and investigate convergence of the aqua-planet climate with increasing resolution.

Timing and structure of the 8.2 kyr B.P. event inferred from δ18O records of stalagmites from China, Oman, and Brazil

Oxygen isotope records of stalagmites from China and Oman reveal a weak summer monsoon event, with a double-plunging structure, that started 8.21 ± 0.02 kyr B.P. An identical but antiphased pattern is also evident in two stalagmite records from eastern Brazil, indicating that the South American Summer Monsoon was intensified during the 8.2 kyr B.P. event. These records demonstrate that the event was of global extent and synchronous within dating errors of <50 years. In comparison with recent model simulations, it is plausible that the 8.2 kyr B.P. event can be tied in changes of the Atlantic Meridional Overturning Circulation triggered by a glacial lake draining event. This, in turn, affected North Atlantic climate and latitudinal position of the Intertropical Convergence Zone, resulting in the observed low-latitude monsoonal precipitation patterns.

Complex latitudinal variation in the morphology of the kleptoparasitic spider Argyrodes kumadai associated with host use and climatic conditions

We examined complex geographical patterns in the morphology of a kleptoparasitic spider, Argyrodes kumadai, across its distributional range in Japan. To disentangle biotic and abiotic factors underlying morphological variation, latitudinal trends were investigated in two traits, body size and relative leg length, across separate transition zones for host use and voltinism. Statistical analyses revealed complex sawtooth clines. Adult body size dramatically changed at the transition zones for host use and voltinism, and exhibited a latitudinal decline following the converse to Bergmann’s cline under the same host use and voltinism in both sexes. A similar pattern was observed for relative leg length in females but not in males. A genetic basis for a part of observed differences in morphology was supported by a common-garden experiment. Our data suggest that local adaptation to factors other than season length such as resource availability (here associated with host use) obscures underlying responses to latitude.