Magnitude of climate variability during middle Pliocene warmth: a palaeoclimate modelling study

The nature and magnitude of climatic variability during the period of middle Pliocene warmth (ca 3.29–2.97 Ma) is poorly understood. We present a suite of palaeoclimate modelling experiments incorporating an advanced atmospheric general circulation model (GCM), coupled to a Q-flux ocean model for 3.29, 3.12 and 2.97 Ma BP. Astronomical solutions for the periods in question were derived from the Berger and Loutre BL2 astronomical solution. Boundary conditions, excluding sea surface temperatures (SSTs) which were predicted by the slab-ocean model, were provided from the USGS PRISM2 2°×2° digital data set. The model results indicate that little annual variation (0.5°C) in SSTs, relative to a ‘control’ experiment, occurred during the middle Pliocene in response to the altered orbital configurations. Annual surface air temperatures also displayed little variation. Seasonally, surface air temperatures displayed a trend of cooler temperatures during December, January and February, and warmer temperatures during June, July and August. This pattern is consistent with altered seasonality resulting from the prescribed orbital configurations. Precipitation changes follow the seasonal trend observed for surface air temperature. Compared to present-day, surface wind strength and wind stress over the North Atlantic, North Pacific and Southern Ocean remained greater in each of the Pliocene experiments. This suggests that wind-driven gyral circulation may have been consistently greater during the middle Pliocene. The trend of climatic variability predicted by the GCM for the middle Pliocene accords with geological data. However, it is unclear if the model correctly simulates the magnitude of the variation. This uncertainty is derived from, (a) the relative insensitivity of the GCM to perturbation in the imposed boundary conditions, (b) a lack of detailed time series data concerning changes to terrestrial ice cover and greenhouse gas concentrations for the middle Pliocene and (c) difficulties in representing the effects of ‘climatic history’ in snap-shot GCM experiments.

Assessment of a climate model to reproduce rainfall variability and extremes over Southern Africa

It is increasingly accepted that any possible climate change will not only have an influence on mean climate but may also significantly alter climatic variability. A change in the distribution and magnitude of extreme rainfall events (associated with changing variability), such as droughts or flooding, may have a far greater impact on human and natural systems than a changing mean. This issue is of particular importance for environmentally vulnerable regions such as southern Africa. The sub-continent is considered especially vulnerable to and ill-equipped (in terms of adaptation) for extreme events, due to a number of factors including extensive poverty, famine, disease and political instability. Rainfall variability and the identification of rainfall extremes is a function of scale, so high spatial and temporal resolution data are preferred to identify extreme events and accurately predict future variability. The majority of previous climate model verification studies have compared model output with observational data at monthly timescales. In this research, the assessment of ability of a state of the art climate model to simulate climate at daily timescales is carried out using satellite-derived rainfall data from the Microwave Infrared Rainfall Algorithm (MIRA). This dataset covers the period from 1993 to 2002 and the whole of southern Africa at a spatial resolution of 0.1° longitude/latitude. This paper concentrates primarily on the ability of the model to simulate the spatial and temporal patterns of present-day rainfall variability over southern Africa and is not intended to discuss possible future changes in climate as these have been documented elsewhere. Simulations of current climate from the UKMeteorological Office Hadley Centre’s climate model, in both regional and global mode, are firstly compared to the MIRA dataset at daily timescales. Secondly, the ability of the model to reproduce daily rainfall extremes is assessed, again by a comparison with extremes from the MIRA dataset. The results suggest that the model reproduces the number and spatial distribution of rainfall extremes with some accuracy, but that mean rainfall and rainfall variability is underestimated (over-estimated) over wet (dry) regions of southern Africa.

Observed trends and teleconnections of the Siberian high: A recently declining center of action

This study investigates variability in the intensity of the wintertime Siberian high (SH) by defining a robust SH index (SHI) and correlating it with selected meteorological fields and teleconnection indices. A dramatic trend of -2.5 hPa decade(-1) has been found in the SHI between 1978 and 2001 with unprecedented (since 1871) low values of the SHI. The weakening of the SH has been confirmed by analyzing different historical gridded analyses and individual station observations of sea level pressure (SLP) and excluding possible effects from the conversion of surface pressure to SLP. SHI correlation maps with various meteorological fields show that SH impacts on circulation and temperature patterns extend far outside the SH source area extending from the Arctic to the tropical Pacific. Advection of warm air from eastern Europe has been identified as the main mechanism causing milder than normal conditions over the Kara and Laptev Seas in association with a strong SH. Despite the strong impacts of the variability in the SH on climatic variability across the Northern Hemisphere, correlations between the SHI and the main teleconnection indices of the Northern Hemisphere are weak. Regression analysis has shown that teleconnection indices are not able to reproduce the interannual variability and trends in the SH. The inclusion of regional surface temperature in the regression model provides closer agreement between the original and reconstructed SHI.

Interactions between mass balance, atmospheric circulation, and recent climate change on the Djankuat Glacier, Caucasus Mountains, Russia

This paper reports recent changes in the mass balance record from the Djankuat Glacier, central greater Caucasus, Russia, and investigates possible relationships between the components of mass balance, local climate, and distant atmospheric forcing. The results clearly show that a strong warming signal has emerged in the central greater Caucasus, particularly since the 1993/1994 mass balance year, and this has led to a significant increase in the summer ablation of Djankuat. At the same time, there has been no compensating consistent increase in winter precipitation and accumulation leading to the strong net loss of mass and increase in glacier runoff. Interannual variability in ablation and accumulation is partly associated with certain major patterns of Northern Hemisphere climatic variability. The positive phase of the North Pacific (NP) teleconnection pattern forces negative geopotential height and temperature anomalies over the Caucasus in summer and results in reduced summer melt, such as in the early 1990s, when positive NP extremes resulted in a temporary decline in ablation rates. The positive phase of the NP is related to El Nino-Southern Oscillation, and it is possible that a teleconnection between the tropical Pacific sea surface temperatures and summer air temperatures in the Caucasus is bridged through the NP pattern. More recently, the NP pattern was predominantly negative, and this distant moderating forcing on summer ablation in the Caucasus was absent. Statistically significant correlations are observed between accumulation and the Scandinavian (SCA) teleconnection pattern. The frequent occurrence of the positive SCA phase at the beginning of accumulation season results in lower than average snowfall and reduced accumulation. The relationship between the North Atlantic Oscillation (NAO), Arctic Oscillation, and accumulation is weak, although positive precipitation anomalies in the winter months are associated with the negative phase of the NAO. A stronger positive correlation is observed between accumulation on Djankuat and geopotential height over the Bay of Biscay unrelated to the established modes of the Northern Hemisphere climatic variability. These results imply that the mass balance of Djankuat is sensitive to the natural variability in the climate system. Distant forcing, however, explains only 16% of the variance in the ablation record and cannot fully explain the recent increase in ablation and negative mass balance.

Fire regimes during the last glacial

Sedimentary charcoal records document changes in fire regime. We have identified 67 sites (30 sites with better than millennial resolution) which have records for some part of the Last Glacial to analyse changes in global fire regimes. Fire was consistently lower during the glacial than during the Eemian and Holocene. Within the glacial, Marine Isotope Stage (MIS) 3 is characterised globally by more fire than MIS 2. The signal for MIS 4 is less clear: there is more fire in the Northern Hemisphere and less fire in the Southern Hemisphere than during MIS 2 and 3. The records, most particularly records from the northern extratropics, show millennial-scale variability in fire regimes corresponding to the rapid climate changes associated with Dansgaard–Oeschger (D-O) cycles. Most of the D-O cycles during the Last Glacial and all of the Heinrich stadials are apparent in the composite global record of fire regime: fire increases during D-O warming events and decreases during intervals of rapid cooling. Our analyses show that fire regimes show a lagged response to rapid climate changes of ca 100–200 years in the case of D-O warming events, ca 0–100 years in the case of D-O cooling events and ca 200 years in the case of Heinrich Stadials. The Strong climatic variability experienced during the glacial resulted in important changes in fire regimes even though the base level of biomass burning was less than today.

Postglacial vegetation, fire and climate dynamics at Central Chilean Patagonia (Lake Shaman, 44ºS)

A 19 cal ka BP pollen and charcoal record from Lake Shaman (44°S; 71°W, Chile) was analyzed to establish vegetation, fire and climate dynamics of the forest-steppe ecotone in Central Chilean Patagonia. Lake Shaman record indicates that the upper Río Cisnes valley was free of ice at around 19 cal ka BP. From this date and until 14.8 cal ka BP, a grass steppe with high proportions of shrubs associated to colder and drier conditions than present developed in this area. A continuous increase of Nothofagus accompanied by a decline in the steppe shrubs and sudden dominance of paludal over aquatic plants from 11 cal ka BP was associated to effective moisture increase but still under modern values. The replacement of the cold-dry grass-shrub steppe by a similar-than-present forest-steppe ecotone suggests an increase in temperature indicating the onset of the Holocene. At the same time, moderate fire activity suggested by the charcoal record could be related to major fuel availability as consequence of Nothofagus forest expansion. Between 8 and 3 cal ka BP, the record indicates the easternmost position of the forest-steppe ecotone suggesting the highest effective moisture with the establishment of seasonality between 5 and 3 cal ka BP. From 3 cal ka BP, the record indicates a retraction of the forest-steppe ecotone accompanied by a high pollen record variability and an increased fire activity. These late changes suggest decreased effective moisture associated with a high climatic variability. At regional and extra-regional scale, climatic changes at Lake Shaman's record are mostly associated to changes (latitudinal shifts and/or strengthening/weakening) of past Southern Westerlies that were previously recorded along Patagonia from the Lateglacial to the mid-Holocene. During the Late Holocene, a regional pattern characterized by high record variability emerges throughout Central Chilean Patagonia. This variability would be related to (1) low magnitude Southern Westerlies changes probably associated to ENSO and/or SAM or (2) the complex relationships between vegetation, fire and human occupations during the last 3 cal ka.

Examining the relationships between Holocene climate change, hydrology, and human society in Ireland

This thesis explores human-environment interactions during the Mid-Late Holocene in raised bogs in central Ireland. The raised bogs of central Ireland are widely-recognised for their considerable palaeoenvironmental and archaeological resources: research over the past few decades has established the potential for such sites to preserve sensitive records of Holocene climatic variability expressed as changes in bog surface wetness (BSW); meanwhile archaeological investigations over the past century have uncovered hundreds of peatland archaeological features dating from the Neolithic through to the Post-Medieval period including wooden trackways, platforms, and deposits of high-status metalwork. Previous studies have attempted to explore the relationship between records of past environmental change and the occurrence of peatland archaeological sites reaching varying conclusions. More recently, environmentally-deterministic models of human-environment interaction in Irish raised bogs at the regional scale have been explicitly tested leading to the conclusion that there is no relationship between BSW and past human activity. These relationships are examined in more detail on a site-by-site basis in this thesis. To that end, testate amoebae-derived BSW records from nine milled former raised bogs in central Ireland were produced from sites with known and dated archaeological records. Relationships between BSW records and environmental conditions within the study area were explored through both the development of a new central Ireland testate amoebae transfer function and through comparisons between recent BSW records and instrumental weather data. Compilation of BSW records from the nine fossil study sites show evidence both for climate forcing, particularly during 3200-2400 cal BP, as well as considerable inter-site variability. Considerable inter-site variability was also evident in the archaeological records of the same sites. Whilst comparisons between BSW and archaeological records do not show a consistent linear relationship, examination of records on a site-by-site basis were shown to reveal interpretatively important contingent relationships. It is concluded therefore, that future research on human-environment interactions should focus on individual sites and should utilise theoretical approaches from the humanities in order to avoid the twin pitfalls of masking important local patterns of change, and of environmental determinism.

Explaining patterns of avian diversity and endemicity: climate and biomes of southern Africa over the last 140,000 years

Aim Test hypotheses that present biodiversity and endemic species richness are related to climatic stability and/or biome persistence.Location Africa south of 15° S. Methods Seventy eight HadCM3 general circulation model palaeoclimate experiments spanning the last 140,000 years, plus a pre-industrial experiment,were used to calculate measures of climatic variability for 0.5° grid cells. Models were fitted relating distributions of the nine biomes of South Africa,Lesotho and Swaziland to present climate. These models were used to simulate potential past biome distribution and extent for the 78 palaeoclimate experiments, and three measures of biome persistence. Climatic response surfaces were fitted for 690 bird species regularly breeding in the region and used to simulate present species richness for cells of the 0.5° grid. Species richness was evaluated for residents, mobile species (nomadic or partially/altitudinally migrant within the region), and intra-African migrants, and also separately for endemic/near-endemic (hereafter ‘endemic’) species as a whole and those associated with each biome. Our hypotheses were tested by analysing correlations between species richness and climatic variability or biome persistence. Results The magnitude of climatic variability showed clear spatial patterns. Marked changes in biome distributions and extents were projected, although limited areas of persistence were projected for some biomes. Overall species richness was not correlated with climatic variability, although richness of mobile species showed a weak negative correlation. Endemic species richness was significantly negatively correlated with climatic variability. Strongest correlations, however, were positive correlations between biome persistence and richness of endemics associated with individual biomes. Main conclusions Low climatic variability, and especially a degree of stability enabling biome persistence, is strongly correlated with species richness of birds endemic to southern Africa. This probably principally reflects reduced extinction risk for these species where the biome to which they are adapted persisted

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.

Emergence of Chenopodium album and Stellaria media of different origins under different climatic conditions

The emergence behaviour of weed species in relation to cultural and meteorological events was studied. Dissimilarities between populations in dormancy and germination ecology, between-year maturation conditions and seed quality and burial site climate all contribute to potentially unpredictable variability. Therefore, a weed emergence data set was produced for weed seeds of Stellaria media and Chenopodium album matured and collected from three populations (Italy, Sweden and UK). The seeds were collected in two consecutive seasons (1999 and 2000) and subsequently buried in the autumn of the same year of maturation in eight contrasting climatic locations throughout Europe and the USA. The experiment sought to explore and explain differences between the three populations in their emergence behaviour. Evidence was demonstrated of synchrony in the timing of the emergence of different populations of a species at a given burial site. The relative magnitudes of emergence from the three populations at a given burial site in a given year were generally similar across all the burial sites in the study. The resulting data set was also used to construct a simple weed emergence model, which was tested for its application to the range of different burial environments and populations. The study demonstrated the possibility of using a simple thermal time-based model to describe part of the emergence behaviour across different burial sites, seed populations and seasons, and a simple winter chilling relationship to adjust for the magnitude of the flush of emergence at a given burial site. This study demonstrates the possibility of developing robust generic models for simple predictions of emergence timing across populations.

Adaptation to climate change and variability: farmer responses to intra-seasonal precipitation trends in South Africa

We describe the nature of recent (50 year) rainfall variability in the summer rainfall zone, South Africa, and how variability is recognised and responded to on the ground by farmers. Using daily rainfall data and self-organising mapping (SOM) we identify 12 internally homogeneous rainfall regions displaying differing parameters of precipitation change. Three regions, characterised by changing onset and timing of rains, rainfall frequencies and intensities, in Limpopo, North West and KwaZulu Natal provinces, were selected to investigate farmer perceptions of, and responses to, rainfall parameter changes. Village and household level analyses demonstrate that the trends and variabilities in precipitation parameters differentiated by the SOM analysis were clearly recognised by people living in the areas in which they occurred. A range of specific coping and adaptation strategies are employed by farmers to respond to climate shifts, some generic across regions and some facilitated by specific local factors. The study has begun to understand the complexity of coping and adaptation, and the factors that influence the decisions that are taken.

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. 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 amplitude of the simulated annual cycle is similar to observations. 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. 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 is not well-correlated with the AMOC at 26N.

Climatic and oceanic associations with daily rainfall extremes over southern Africa

Changes in climate variability and, in particular, changes in extreme climate events are likely to be of far more significance for environmentally vulnerable regions than changes in the mean state. It is generally accepted that sea-surface temperatures (SSTs) play an important role in modulating rainfall variability. Consequently, SSTs can be prescribed in global and regional climate modelling in order to study the physical mechanisms behind rainfall and its extremes. Using a satellite-based daily rainfall historical data set, this paper describes the main patterns of rainfall variability over southern Africa, identifies the dates when extreme rainfall occurs within these patterns, and shows the effect of resolution in trying to identify the location and intensity of SST anomalies associated with these extremes in the Atlantic and southwest Indian Ocean. Derived from a Principal Component Analysis (PCA), the results also suggest that, for the spatial pattern accounting for the highest amount of variability, extremes extracted at a higher spatial resolution do give a clearer indication regarding the location and intensity of anomalous SST regions. As the amount of variability explained by each spatial pattern defined by the PCA decreases, it would appear that extremes extracted at a lower resolution give a clearer indication of anomalous SST regions.

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.