How well can we simulate past climates? Evaluating the models using global palaeoenvironmental datasets.

Global syntheses of palaeoenvironmental data are required to test climate models under conditions different from the present. Data sets for this purpose contain data from spatially extensive networks of sites. The data are either directly comparable to model output or readily interpretable in terms of modelled climate variables. Data sets must contain sufficient documentation to distinguish between raw (primary) and interpreted (secondary, tertiary) data, to evaluate the assumptions involved in interpretation of the data, to exercise quality control, and to select data appropriate for specific goals. Four data bases for the Late Quaternary, documenting changes in lake levels since 30 kyr BP (the Global Lake Status Data Base), vegetation distribution at 18 kyr and 6 kyr BP (BIOME 6000), aeolian accumulation rates during the last glacial-interglacial cycle (DIRTMAP), and tropical terrestrial climates at the Last Glacial Maximum (the LGM Tropical Terrestrial Data Synthesis) are summarised. Each has been used to evaluate simulations of Last Glacial Maximum (LGM: 21 calendar kyr BP) and/or mid-Holocene (6 cal. kyr BP) environments. Comparisons have demonstrated that changes in radiative forcing and orography due to orbital and ice-sheet variations explain the first-order, broad-scale (in space and time) features of global climate change since the LGM. However, atmospheric models forced by 6 cal. kyr BP orbital changes with unchanged surface conditions fail to capture quantitative aspects of the observed climate, including the greatly increased magnitude and northward shift of the African monsoon during the early to mid-Holocene. Similarly, comparisons with palaeoenvironmental datasets show that atmospheric models have underestimated the magnitude of cooling and drying of much of the land surface at the LGM. The inclusion of feedbacks due to changes in ocean- and land-surface conditions at both times, and atmospheric dust loading at the LGM, appears to be required in order to produce a better simulation of these past climates. The development of Earth system models incorporating the dynamic interactions among ocean, atmosphere, and vegetation is therefore mandated by Quaternary science results as well as climatological principles. For greatest scientific benefit, this development must be paralleled by continued advances in palaeodata analysis and synthesis, which in turn will help to define questions that call for new focused data collection efforts.

The neolithisation of Liguria (NW Italy): an environmental archaeological and palaeoenvironmental perspective

The archaeological evidence compiled for Liguria has enabled the formulation of a comprehensive model of Neolithic social, technological and economic development (∼7800–5700 cal yrs BP). The model indicates that during the Early and Middle Neolithic (∼7800–6300 cal yrs BP; ‘Impressed Ware’ and ‘Square Mouthed’ pottery cultures) human activity mainly focussed on low (coastal) and mid-altitude areas. By the Late Neolithic (∼6300–5700 cal yrs BP; ‘Chassey’ culture) farming practices were taking place over a wider range of altitudes and involved transhumant pastoralism. Complementary environmental archaeological and palaeoecological records from caves, open-air sites, lakes and mires indicate that human activities had a more significant impact on the environment than previously thought. This included clearance, especially Abies, Ulmus, Fraxinus and Tilia, and woodland utilisation and management (e.g. leaf foddering), as well as cereal cultivation and animal husbandry. The influence of Middle Holocene climatic changes, especially from ∼7800 cal yrs BP, on the direction of vegetation changes and socio-economic developments during the Neolithic remain uncertain.

Extreme rainfall variability in Australia: Patterns, drivers and predictability

Leading patterns of observed monthly extreme rainfall variability in Australia are examined using an Empirical Orthogonal Teleconnection (EOT) method. Extreme rainfall variability is more closely related to mean rainfall variability during austral summer than in winter. The leading EOT patterns of extreme rainfall explain less variance in Australia-wide extreme rainfall than is the case for mean rainfall EOTs. We illustrate that, as with mean rainfall, the El Niño-Southern Oscillation (ENSO) has the strongest association with warm-season extreme rainfall variability, while in the cool-season the primary drivers are atmospheric blocking and the subtropical ridge. The Indian Ocean Dipole and Southern Annular Mode also have significant relationships with patterns of variability during austral winter and spring. Leading patterns of summer extreme rainfall variability have predictability several months ahead from Pacific sea surface temperatures (SSTs) and as much as a year in advance from Indian Ocean SSTs. Predictability from the Pacific is greater for wetter than average summer months than for months that are drier than average, whereas for the Indian Ocean the relationship has greater linearity. Several cool-season EOTs are associated with mid-latitude synoptic-scale patterns along the south and east coasts. These patterns have common atmospheric signatures denoting moist onshore flow and strong cyclonic anomalies often to the north of a blocking anti-cyclone. Tropical cyclone activity is observed to have significant relationships with some warm season EOTs. This analysis shows that extreme rainfall variability in Australia can be related to remote drivers and local synoptic-scale patterns throughout the year.

Uso del espacio de cazadores recolectores y paleoambiente Holoceno en el Valle del Río Cisnes de Aisén, Chile

This paper presents results on archaeological research conducted at the Cisnes river basin (~44° S), valley which passes through several environments in western Patagonia, from the westernmost limits of the steppe to the Pacific channels. These are assessed in light of a palaeoenvironmental reconstruction spanning from the Late Pleistocene to the Holocene. Results allow interpreting different ways for approaching the environment; these are exposed both spatially and chronologically. Cultural units defined occupied the space discontinuously and ranking it differentially attending to environmental variability. We correlate our data with available archaeological information for Aisén and adjacent areas. Additionally, we discuss cultural continuity and eventual interaction of steppe hunter-gatherers with peoples from the western archipelagic area.

The surface temperatures of Earth: steps towards integrated understanding of variability and change

Surface temperature is a key aspect of weather and climate, but the term may refer to different quantities that play interconnected roles and are observed by different means. In a community-based activity in June 2012, the EarthTemp Network brought together 55 researchers from five continents to improve the interaction between scientific communities who focus on surface temperature in particular domains, to exploit the strengths of different observing systems and to better meet the needs of different communities. The workshop identified key needs for progress towards meeting scientific and societal requirements for surface temperature understanding and information, which are presented in this community paper. A "whole-Earth" perspective is required with more integrated, collaborative approaches to observing and understanding Earth's various surface temperatures. It is necessary to build understanding of the relationships between different surface temperatures, where presently inadequate, and undertake large-scale systematic intercomparisons. Datasets need to be easier to obtain and exploit for a wide constituency of users, with the differences and complementarities communicated in readily understood terms, and realistic and consistent uncertainty information provided. Steps were also recommended to curate and make available data that are presently inaccessible, develop new observing systems and build capacities to accelerate progress in the accuracy and usability of surface temperature datasets.

Mechanisms for AMOC variability simulated by the NEMO model

We have investigated mechanisms for the Atlantic Meridional Overturning Circulation (AMOC) variability at 26.5° N (other than the Ekman component) that can be related to external forcings, in particular wind variability. Resolution dependence is studied using identical experiments with 1° and 1/4° NEMO model runs over 1960–2010. The analysis shows that much of the variability in the AMOC at 26° N can be related to the wind strength over the North Atlantic, through mechanisms lagged on different timescales. At ~ 1-year lag the January–June difference of mean sea level pressure between high and mid-latitudes in the North Atlantic explains 35–50% of the interannual AMOC variability (with negative correlation between wind strength and AMOC). At longer lead timescales ~ 4 years, strong (weak) winds over the northern North Atlantic (specifically linked to the NAO index) are followed by higher (lower) AMOC transport, but this mechanism only works in the 1/4° model. Analysis of the density correlations suggests an increase (decrease) in deep water formation in the North Atlantic subpolar gyre to be the cause. Therefore another 30% of the AMOC variability at 26° N can be related to density changes in the top 1000 m in the Labrador and Irminger seas occurring ~ 4 years earlier.

Contrasting interannual and multidecadal NAO variability

Decadal and longer timescale variability in the winter North Atlantic Oscillation (NAO) has considerable impact on regional climate, yet it remains unclear what fraction of this variability is potentially predictable. This study takes a new approach to this question by demonstrating clear physical differences between NAO variability on interannual-decadal (<30 year) and multidecadal (>30 year) timescales. It is shown that on the shorter timescale the NAO is dominated by variations in the latitude of the North Atlantic jet and storm track, whereas on the longer timescale it represents changes in their strength instead. NAO variability on the two timescales is associated with different dynamical behaviour in terms of eddy-mean flow interaction, Rossby wave breaking and blocking. The two timescales also exhibit different regional impacts on temperature and precipitation and different relationships to sea surface temperatures. These results are derived from linear regression analysis of the Twentieth Century and NCEP-NCAR reanalyses and of a high-resolution HiGEM General Circulation Model control simulation, with additional analysis of a long sea level pressure reconstruction. Evidence is presented for an influence of the ocean circulation on the longer timescale variability of the NAO, which is particularly clear in the model data. As well as providing new evidence of potential predictability, these findings are shown to have implications for the reconstruction and interpretation of long climate records.

On the relationship between decadal buoyancy anomalies and variability of the Atlantic meridional overturning circulation

Owing to the role of the Atlantic meridional overturning circulation (AMOC) in ocean heat transport, AMOC variability is thought to play a role in climate variability on a wide range of time scales. This paper focuses on the potential role of the AMOC in climate variability on decadal time scales. Coupled and ocean-only general circulation models run in idealized geometries are utilized to study the relationships between decadal AMOC and buoyancy variability and determine whether the AMOC plays an active role in setting sea surface temperature on decadal time scales.DecadalAMOC variability is related to changes in the buoyancy field along the western boundary according to the thermal wind relation. Buoyancy anomalies originate in the upper ocean of the subpolar gyre and travel westward as baroclinic Rossby waves. When the buoyancy anomalies strike the western boundary, they are advected southward by the deep western boundary current, leading to latitudinally coherent AMOC variability. The AMOC is observed to respond passively to decadal buoyancy anomalies: although variability of the AMOC leads to meridional ocean heat transport anomalies, these transports are not responsible for creating the buoyancy anomalies in the subpolar gyre that drive AMOC variability.

Mean climate and variability of the atmosphere and ocean on an aquaplanet

Numerical experiments are described that pertain to the climate of a coupled atmosphere–ocean–ice system in the absence of land, driven by modern-day orbital and CO2 forcing. Millennial time-scale simulations yield a mean state in which ice caps reach down to 55° of latitude and both the atmosphere and ocean comprise eastward- and westward-flowing zonal jets, whose structure is set by their respective baroclinic instabilities. Despite the zonality of the ocean, it is remarkably efficient at transporting heat meridionally through the agency of Ekman transport and eddy-driven subduction. Indeed the partition of heat transport between the atmosphere and ocean is much the same as the present climate, with the ocean dominating in the Tropics and the atmosphere in the mid–high latitudes. Variability of the system is dominated by the coupling of annular modes in the atmosphere and ocean. Stochastic variability inherent to the atmospheric jets drives variability in the ocean. Zonal flows in the ocean exhibit decadal variability, which, remarkably, feeds back to the atmosphere, coloring the spectrum of annular variability. A simple stochastic model can capture the essence of the process. Finally, it is briefly reviewed how the aquaplanet can provide information about the processes that set the partition of heat transport and the climate of Earth.

The importance of vertical velocity variability for estimates of the indirect aerosol effects

The activation of aerosols to form cloud droplets is dependent upon vertical velocities whose local variability is not typically resolved at the GCM grid scale. Consequently, it is necessary to represent the subgrid-scale variability of vertical velocity in the calculation of cloud droplet number concentration. This study uses the UK Chemistry and Aerosols community model (UKCA) within the Hadley Centre Global Environmental Model (HadGEM3), coupled for the first time to an explicit aerosol activation parameterisation, and hence known as UKCA-Activate. We explore the range of uncertainty in estimates of the indirect aerosol effects attributable to the choice of parameterisation of the subgrid-scale variability of vertical velocity in HadGEM-UKCA. Results of simulations demonstrate that the use of a characteristic vertical velocity cannot replicate results derived with a distribution of vertical velocities, and is to be discouraged in GCMs. This study focuses on the effect of the variance (σw2) of a Gaussian pdf (probability density function) of vertical velocity. Fixed values of σw (spanning the range measured in situ by nine flight campaigns found in the literature) and a configuration in which σw depends on turbulent kinetic energy are tested. Results from the mid-range fixed σw and TKE-based configurations both compare well with observed vertical velocity distributions and cloud droplet number concentrations. The radiative flux perturbation due to the total effects of anthropogenic aerosol is estimated at −1.9 W m−2 with σw = 0.1 m s−1, −2.1 W m−2 with σw derived from TKE, −2.25 W m−2 with σw = 0.4 m s−1, and −2.3 W m−2 with σw = 0.7 m s−1. The breadth of this range is 0.4 W m−2, which is comparable to a substantial fraction of the total diversity of current aerosol forcing estimates. Reducing the uncertainty in the parameterisation of σw would therefore be an important step towards reducing the uncertainty in estimates of the indirect aerosol effects. Detailed examination of regional radiative flux perturbations reveals that aerosol microphysics can be responsible for some climate-relevant radiative effects, highlighting the importance of including microphysical aerosol processes in GCMs.

Temporal and spatial variability of auroral forms in the 10–14 MLT sector: Relationship to plasma convection and solar wind-magnetosphere coupling

Ground-based observations of dayside auroral forms and magnetic perturbations in the arctic sectors of Svalbard and Greenland, in combination with the high-resolution measurements of ionospheric ion drift and temperature by the EISCAT radar, are used to study temporal/spatial structures of cusp-type auroral forms in relation to convection. Large-scale patterns of equivalent convection in the dayside polar ionosphere are derived from the magnetic observations in Greenland and Svalbard. This information is used to estimate the ionospheric convection pattern in the vicinity of the cusp/cleft aurora. The reported observations, covering the period 0700-1130 UT, on January 11, 1993, are separated into four intervals according to the observed characteristics of the aurora and ionospheric convection. The morphology and intensity of the aurora are very different in quiet and disturbed intervals. A latitudinally narrow zone of intense and dynamical 630.0 nm emission equatorward of 75 degrees MLAT, was observed during periods of enhanced antisunward convection in the cusp region. This (type 1 cusp aurora) is considered to be the signature of plasma entry via magnetopause reconnection at low magnetopause latitudes, i.e. the low-latitude boundary layer (LLB I,). Another zone of weak 630.0 nm emission (type 2 cusp aurora) was observed to extend up to high latitudes (similar to 79 degrees MLAT) during relatively quiet magnetic conditions, when indications of reverse (sunward) convection was observed in the dayside polar cap. This is postulated to be a signature of merging between a northward directed IMF (B-z > 0) and the geomagnetic field poleward of the cusp. The coexistence of type 1 and 2 auroras was observed under intermediate circumstances. The optical observations from Svalbard and Greenland were also used to determine the temporal and spatial evolution of type 1 auroral forms, i.e. poleward-moving auroral events occurring in the vicinity of a rotational convection reversal in the early post-noon sector. Each event appeared as a local brightening at the equatorward boundary of the pre-existing type 1 cusp aurora, followed by poleward and eastward expansions of luminosity. The auroral events were associated with poleward-moving surges of enhanced ionospheric convection and F-layer ion temperature as observed by the EISCAT radar in Tromso. The EISCAT ion flow data in combination with the auroral observations show strong evidence for plasma flow across the open/closed field line boundary.

Variability of dayside high latitude convection associated with a sequence of auroral transients

10 second resolution ionospheric convection data covering the invariant latitude range from 71° to 76°, obtained by using the EISCAT UHF and VHF radars, are combined with optical data from Ny Ålesund during a sequence of auroral transients in the post-noon sector (∼ 15 MLT). Satellite observations of polar cap convection patterns suggest negative BZ and BY components of the interplanetary magnetic field. Burst-like enhancements of westward (sunward) post-noon convection were accompanied by eastward moving auroral forms at higher latitudes, above the convection reversal boundary. In this case the background convection was weak, whereas the integrated potential drop across the radar field-of-view associated with the westward flow bursts was typically ∼ 20-35 kV. The auroral phenomenon consists of a series of similar events with a mean repetition period of 8 min. A close correlation between the auroral activity and convection enhancements in the cleft ionosphere is demonstrated.

Variability of dayside convection and motions of the cusp/cleft aurora

We present measurements of the ionospheric plasma flow over the range of invariant latitudes 71–76°, observed at 10-second resolution using both the EISCAT radars, with simultaneous observations of the 630 nm cusp/cleft aurora made by a meridian-scanning photometer at Ny Ålesund, Svalbard. A major increase in the trans-auroral voltage from 5 to 40 kV (associated with sunward convection in the early afternoon sector) is found to follow a southward motion of the aurora and coincide with the onset of regular transient auroral breakup events. It is shown that these observations are consistent with recent theoretical work on how ionospheric flows are excited by time-dependent reconnection at the dayside magnetopause.