Variability of the Atlantic thermohaline circulation described by three-dimensional empirical orthogonal functions

We describe the use of bivariate 3d empirical orthogonal functions (EOFs) in characterising low frequency variability of the Atlantic thermohaline circulation (THC) in the Hadley Centre global climate model, HadCM3. We find that the leading two modes are well correlated with an index of the meridional overturning circulation (MOC) on decadal timescales, with the leading mode alone accounting for 54% of the decadal variance. Episodes of coherent oscillations in the sub-space of the leading EOFs are identified; these episodes are of great interest for the predictability of the THC, and could indicate the existence of different regimes of natural variability. The mechanism identified for the multi-decadal variability is an internal ocean mode, dominated by changes in convection in the Nordic Seas, which lead the changes in the MOC by a few years. Variations in salinity transports from the Arctic and from the North Atlantic are the main feedbacks which control the oscillation. This mode has a weak feedback onto the atmosphere and hence a surface climatic influence. Interestingly, some of these climate impacts lead the changes in the overturning. There are also similarities to observed multi-decadal climate variability.

Bjerknes compensation and the decadal variability of the energy transports in a coupled climate model

In the 1960s, Jacob Bjerknes suggested that if the top-of-the-atmosphere (TOA) fluxes and the oceanic heat storage did not vary too much, then the total energy transport by the climate system would not vary too much either. This implies that any large anomalies of oceanic and atmospheric energy transport should be equal and opposite. This simple scenario has become known as Bjerknes compensation. A long control run of the Third Hadley Centre Coupled Ocean-Atmosphere General Circulation Model (HadCM3) has been investigated. It was found that northern extratropical decadal anomalies of atmospheric and oceanic energy transports are significantly anticorrelated and have similar magnitudes, which is consistent with the predictions of Bjerknes compensation. ne degree of compensation in the northern extratropics was found to increase with increasing, time scale. Bjerknes compensation did not occur in the Tropics, primarily as large changes in the surface fluxes were associated with large changes in the TOA fluxes. In the ocean, the decadal variability of the energy transport is associated with fluctuations in the meridional overturning circulation in the Atlantic Ocean. A stronger Atlantic Ocean energy transport leads to strong warming of surface temperatures in the Greenland-Iceland-Norwegian (GIN) Seas. which results in a reduced equator-to-pole surface temperature gradient and reduced atmospheric baroclinicity. It is argued that a stronger Atlantic Ocean energy transport leads to a weakened atmospheric transient energy transport.

The variation of solar wind correlation lengths over three solar cycles

We present the results of a study of solar wind velocity and magnetic field correlation lengths over the last 35 years. The correlation length of the magnetic field magnitude λ | B| increases on average by a factor of two at solar maxima compared to solar minima. The correlation lengths of the components of the magnetic field λ_{B_{XYZ}} and of the velocity λ_{V_{YZ}} do not show this change and have similar values, indicating a continual turbulent correlation length of around 1.4×106 km. We conclude that a linear relation between λ | B|, VB 2, and Kp suggests that the former is related to the total magnetic energy in the solar wind and an estimate of the average size of geoeffective structures, which is, in turn, proportional to VB 2. By looking at the distribution of daily correlation lengths we show that the solar minimum values of λ | B| correspond to the turbulent outer scale. A tail of larger λ | B| values is present at solar maximum causing the increase in mean value.

A comparison of the variability of biological nutrients against depth and potential density.

The main biogeochemical nutrient distributions, along with ambient ocean temperature and the light field, control ocean biological productivity. Observations of nutrients are much sparser than physical observations of temperature and salinity, yet it is critical to validate biogeochemical models against these sparse observations if we are to successfully model biological variability and trends. Here we use data from the Bermuda Atlantic Time-series Study and the World Ocean Database 2005 to demonstrate quantitatively that over the entire globe a significant fraction of the temporal variability of phosphate, silicate and nitrate within the oceans is correlated with water density. The temporal variability of these nutrients as a function of depth is almost always greater than as a function of potential density, with he largest reductions in variability found within the main pycnocline. The greater nutrient variability as a function of depth occurs when dynamical processes vertically displace nutrient and density fields together on shorter timescales than biological adjustments. These results show that dynamical processes can have a significant impact on the instantaneous nutrient distributions. These processes must therefore be considered when modeling biogeochemical systems, when comparing such models with observations, or when assimilating data into such models.

Seasonal and interannual variability in the Mozambique Channel from moored current observations

Direct observations from an array of current meter moorings across the Mozambique Channel in the south-west Indian Ocean are presented covering a period of more than 4 years. This allows an analysis of the volume transport through the channel, including the variability on interannual and seasonal time scales. The mean volume transport over the entire observational period is 16.7 Sv poleward. Seasonal variations have a magnitude of 4.1 Sv and can be explained from the variability in the wind field over the western part of the Indian Ocean. Interannual variability has a magnitude of 8.9 Sv and is large compared to the mean. This time scale of variability could be related to variability in the Indian Ocean Dipole (IOD), showing that it forms part of the variability in the ocean-climate system of the entire Indian Ocean. By modulating the strength of the South Equatorial Current, the weakening (strengthening) tropical gyre circulation during a period of positive (negative) IOD index leads to a weakened (strengthened) southward transport through the channel, with a time lag of about a year. The relatively strong interannual variability stresses the importance of long-term direct observations.

Effects of solar wind magnetosphere coupling recorded at different geomagnetic latitudes: Separation of directly-driven and storage/release systems

The effect on geomagnetic activity of solar wind speed, compared with that of the strength of the interplanetary magnetic field, differs with geomagnetic latitude. In this study we construct a new index based on monthly standard deviations in the H-component of the geomagnetic field for all geomagnetic latitudes. We demonstrate that for this index the response at auroral regions correlates best with interplanetary coupling functions which include the solar wind speed while mid- and low-latitude regions respond to variations in the interplanetary magnetic field strength. These results are used to isolate the responsible geomagnetic current systems.

The impact of two arable field margin management schemes on litter decomposition

A primary objective of agri-environment schemes is the conservation of biodiversity; in addition to increasing the value of farmland for wildlife, these schemes also aim to restore natural ecosystem functioning. The management of scheme options can influence their value for delivering ecosystem services by modifying the composition of floral and faunal communities. This study examines the impact of an agri-environment scheme prescription on ecosystem functioning by testing the hypothesis that vegetation management influences decomposition rates in grassy arable field margins. The effects of two vegetation management practices in arable field margins - cutting and soil disturbance (scarification) - on litter decomposition were compared using a litterbag experimental approach in early April 2006. Bags had either small mesh designed to restrict access to soil macrofauna, or large mesh that would allow macrofauna to enter. Bags were positioned on the soil surface or inserted into the soil in cut and scarified margins, retrieved after 44, 103 and 250 days and the amount of litter mass remaining was calculated. Litter loss from the litterbags with large mesh was greater than from the small mesh bags, providing evidence that soil macrofauna accelerate rates of litter decomposition. In the large mesh bags, the proportion of litter remaining in bags above and belowground in the cut plots was similar, while in the scarified plots, there was significantly more litter left in the aboveground bags than in the belowground bags. This loss of balance between decomposition rates above and belowground in scarified margins may have implications for the development and maintenance of grassy arable field margins by influencing nutrient availability for plant communities. (C) 2008 Elsevier B.V. All rights reserved.

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 importance of interplanetary magnetic field (By)on polar cap patch formation

A number of poleward moving events were observed between 1130 and 1300 UT on 11 February 2004, during periods of southward interplanetary magnetic field (IMF), while the steerable antenna of the European Incoherent Scatter (EISCAT) Svalbard radar (ESR)and the Tromsø VHF radar pointed nearly northward at low elevation. In this interval, simultaneous SuperDARN CUTLASS Finland radar measurements showed poleward moving radar aurora forms (PMRAFs) which appeared very similar to the density enhancements observed by the ESR northward pointing antenna. These events appeared quasiperiodically with a period of about 10 min. Comparing the observations from the above three radars, it is inferred that there is an almost one‐to‐one correspondence between the poleward moving plasma concentration enhancements (PMPCEs) observed by the ESR and the VHF radar and the PMRAFs measured by the CUTLASS Finland radar. These observations are consistent with the interpretation that the polar cap patch material was generated by photoionization at subauroral latitudes and that the plasma was structured by bursts of magnetopause reconnection giving access to the polar cap. There is clear evidence that plasma structuring into patches was dependent on the variability in IMF |By|. The duration of these events implies that the average evolution time of the newly opened flux tubes from the subauroral region to the polar cap was about 33 min.

High frequency variability of the Atlantic meridional overturning circulation

We compare the variability of the Atlantic meridional overturning circulation (AMOC) as simulated by the coupled climate models of the RAPID project, which cover a wide range of resolution and complexity, and observed by the RAPID/MOCHA array at about 26N. We analyse variability on a range of timescales, from five-daily to interannual. In models of all resolutions there is substantial variability on timescales of a few days; in most AOGCMs the amplitude of the variability is of somewhat larger magnitude than that observed by the RAPID array, while the time-mean is within about 10% of the observational estimate. The amplitude of the simulated annual cycle is similar to observations, but the shape of the annual cycle shows a spread among the models. A dynamical decomposition shows that in the models, as in observations, the AMOC is predominantly geostrophic (driven by pressure and sea-level gradients), with both geostrophic and Ekman contributions to variability, the latter being exaggerated and the former underrepresented in models. Other ageostrophic terms, neglected in the observational estimate, are small but not negligible. The time-mean of the western boundary current near the latitude of the RAPID/MOCHA array has a much wider model spread than the AMOC does, indicating large differences among models in the simulation of the wind-driven gyre circulation, and its variability is unrealistically small in the models. In many RAPID models and in models of the Coupled Model Intercomparison Project Phase 3 (CMIP3), interannual variability of the maximum of the AMOC wherever it lies, which is a commonly used model index, is similar to interannual variability in the AMOC at 26N. Annual volume and heat transport timeseries at the same latitude are well-correlated within 15--45N, indicating the climatic importance of the AMOC. In the RAPID and CMIP3 models, we show that the AMOC is correlated over considerable distances in latitude, but not the whole extent of the north Atlantic; consequently interannual variability of the AMOC at 50N, where it is particularly relevant to European climate, is not well-correlated with that of the AMOC at 26N, where it is monitored by the RAPID/MOCHA array.

Resolving the spatial variability of soil N using fractions of soil organic matter

The spatial variability of soil nitrogen (N) mineralisation has not been extensively studied, which limits our capacity to make N fertiliser recommendations. Even less attention has been paid to the scale-dependence of the variation. The objective of this research was to investigate the scale-dependence of variation of mineral N (MinN, N–NO3− plus N–NH4+) at within-field scales. The study was based on the spatial dependence of the labile fractions of SOM, the key fractions for N mineralisation. Soils were sampled in an unbalanced nested design in a 4-ha arable field to examine the distribution of the variation of SOM at 30, 10, 1, and 0.12 m. Organic matter in free and intra-aggregate light fractions (FLF and IALF) was extracted by physical fractionation. The variation occurred entirely within 0.12 m for FLF and at 10 m for IALF. A subsequent sampling on a 5-m grid was undertaken to link the status of the SOM fractions to MinN, which showed uncorrelated spatial dependence. A uniform application of N fertiliser would be suitable in this case. The failure of SOM fractions to identify any spatial dependence of MinN suggests that other soil variables, or crop indicators, should be tested to see if they can identify different N supply areas within the field for a more efficient and environmentally friendly N management.

Seasonality and interannual variability of the westerly jet in the Tibetan Plateau region

In this study, 40-yr ECMWF Re-Analysis (ERA-40) data are used for the description of the seasonal cycle and the interannual variability of the westerly jet in the Tibetan Plateau region. To complement results based on the analysis of monthly mean horizontal wind speeds, an occurrence-based jet climatology is constructed by identifying the locations of the jet axes at 6-hourly intervals throughout 1958–2001. Thus, a dataset describing the highly transient and localized features of jet variability is obtained. During winter and summer the westerly jet is located, respectively, to the south and north of the Tibetan Plateau. During the spring and autumn seasons there are jet transitions from south to north and vice versa. The median dates for these transitions are 28 April and 12 October. The spring transition is associated with large interannual variations, while the fall transition occurs more reliably within a 3-week period. The strength of the jet exhibits a peculiar seasonal cycle. During northward migration in April/May, the jet intensity weakens and its latitudinal position varies largely. In some springs, there are several transitions and split configurations occur before the jet settles in its northern summer position. In June, a well-defined and unusually strong jet reappears at the northern flanks of the Tibetan Plateau. In autumn, the jet gradually but reliably recedes to the south and is typically more intense than in spring. The jet transitions between the two preferred locations follow the seasonal latitudinal migration of the jet in the Northern Hemisphere. An analysis of interannual variations shows the statistical relationship between the strength of the summer jet, the tropospheric meridional temperature gradient, and the all-India rainfall series. Both this analysis and results from previous studies point to the particular dynamical relevance of the onsetting Indian summer monsoon precipitation and the associated diabatic heating for the formation of the strong summer jet. Finally, an example is provided that illustrates the climatological significance of the jet in terms of the covariation between the jet location and the spatial precipitation distribution in central Asia.

The impact of resolution on the adjustment and decadal variability of the Atlantic Meridional Overturning Circulation in a coupled climate model

Variations in the Atlantic Meridional Overturning Circulation (MOC) exert an important influence on climate, particularly on decadal time scales. Simulation of the MOC in coupled climate models is compromised, to a degree that is unknown, by their lack of fidelity in resolving some of the key processes involved. There is an overarching need to increase the resolution and fidelity of climate models, but also to assess how increases in resolution influence the simulation of key phenomena such as the MOC. In this study we investigate the impact of significantly increasing the (ocean and atmosphere) resolution of a coupled climate model on the simulation of MOC variability by comparing high and low resolution versions of the same model. In both versions, decadal variability of the MOC is closely linked to density anomalies that propagate from the Labrador Sea southward along the deep western boundary. We demonstrate that the MOC adjustment proceeds more rapidly in the higher resolution model due the increased speed of western boundary waves. However, the response of the Atlantic Sea Surface Temperatures (SSTs) to MOC variations is relatively robust - in pattern if not in magnitude - across the two resolutions. The MOC also excites a coupled ocean-atmosphere response in the tropical Atlantic in both model versions. In the higher resolution model, but not the lower resolution model, there is evidence of a significant response in the extratropical atmosphere over the North Atlantic 6 years after a maximum in the MOC. In both models there is evidence of a weak negative feedback on deep density anomalies in the Labrador Sea, and hence on the MOC (with a time scale of approximately ten years). Our results highlight the need for further work to understand the decadal variability of the MOC and its simulation in climate models.

Seasonal temperature variability of the Neoglacial (3300–2500 BP) and Roman Warm Period (2500–1600 BP) reconstructed from oxygen isotope ratios of limpet shells (Patella vulgata), Northwest Scotland

Seasonal sea-surface temperaturevariability for the Neoglacial (3300–2500 BP) and Roman WarmPeriod (RWP; 2500–1600 BP), which correspond to the Bronze and Iron Ages, respectively, was estimated using oxygen isotope ratios obtained from high-resolution samples micromilled from radiocarbon-dated, archaeological limpet (Patella vulgata) shells. The coldest winter months recorded in Neoglacial shells averaged 6.6 ± 0.3 °C, and the warmest summer months averaged 14.7 ± 0.4 °C. One Neoglacial shell captured a year without a summer, which may have resulted from a dust veil from a volcanic eruption in the Katla volcanic system in Iceland. RWP shells record average winter and summer monthly temperatures of 6.3 ± 0.1 °C and 13.3 ± 0.3 °C, respectively. These results capture a cooling transition from the Neoglacial to RWP, which is further supported by earlier studies of pine history in Scotland, pollen type analyses in northeast Scotland, and European glacial events. The cooling transition observed at the boundary between the Neoglacial and RWP in our study also agrees with the abrupt climate deterioration at 2800–2700 BP (also referred to as the Subboreal/Subatlantic transition) and therefore may have been driven by decreased solar radiation and weakened North Atlantic Oscillation conditions.

Model studies of the interannual variability of the northern-hemisphere stratospheric winter circulation: the role of the quasi-biennial oscillation

A series of experiments are described that examine the sensitivity of the northern-hemisphere winter evolution to the equatorial quasi-biennial oscillation (QBO). The prime tool for the experiments is a stratosphere-mesosphere model. The model is integrated over many years with the modelled equatorial winds relaxed towards observed values in order to simulate a realistic QBO. In experiment A the equatorial winds are relaxed towards Singapore radiosonde observations in the height region 16-32 km. In contrast to previous modelling studies, the Holton-Tan relationship (warm/cold winters associated with easterly/westerly QBO winds in the lower stratosphere) is absent. However, in a second experiment (run B) in which the equatorial winds are relaxed towards rocketsonde data over the extended height range 16-58 km, a realistic Holton-Tan relationship is reproduced. A series of further studies are described that explore in more detail the sensitivity to various equatorial height regions and to the bottom-boundary forcing. The experiments suggest that the evolution of the northern-hemisphere winter circulation is sensitive to equatorial winds throughout the whole depth of the stratosphere and not just to the lower-stratospheric wind direction as previously assumed.