Anthropogenic greenhouse gas contribution to flood risk in England and Wales in autumn 2000

Interest in attributing the risk of damaging weather-related events to anthropogenic climate change is increasing1. Yet climate models used to study the attribution problem typically do not resolve the weather systems associated with damaging events2 such as the UK floods of October and November 2000. Occurring during the wettest autumn in England and Wales since records began in 17663, 4, these floods damaged nearly 10,000 properties across that region, disrupted services severely, and caused insured losses estimated at £1.3 billion (refs 5, 6). Although the flooding was deemed a ‘wake-up call’ to the impacts of climate change at the time7, such claims are typically supported only by general thermodynamic arguments that suggest increased extreme precipitation under global warming, but fail8, 9 to account fully for the complex hydrometeorology4, 10 associated with flooding. Here we present a multi-step, physically based ‘probabilistic event attribution’ framework showing that it is very likely that global anthropogenic greenhouse gas emissions substantially increased the risk of flood occurrence in England and Wales in autumn 2000. Using publicly volunteered distributed computing11, 12, we generate several thousand seasonal-forecast-resolution climate model simulations of autumn 2000 weather, both under realistic conditions, and under conditions as they might have been had these greenhouse gas emissions and the resulting large-scale warming never occurred. Results are fed into a precipitation-runoff model that is used to simulate severe daily river runoff events in England and Wales (proxy indicators of flood events). The precise magnitude of the anthropogenic contribution remains uncertain, but in nine out of ten cases our model results indicate that twentieth-century anthropogenic greenhouse gas emissions increased the risk of floods occurring in England and Wales in autumn 2000 by more than 20%, and in two out of three cases by more than 90%.

Climate on the Southern Black Sea coast during the Holocene: implications from the Sofular Cave record

We present the updated Holocene section of the Sofular Cave record from the southernBlackSeacoast (northern Turkey); an area with considerably different present-day climate compared to that of the neighboring Eastern Mediterranean region. Stalagmite δ13C, growth rates and initial (234U/238U) ratios provide information about hydrological changes above the cave; and prove to be more useful than δ18O for deciphering Holocene climatic variations. Between ∼9.6 and 5.4 ka BP (despite a pause from ∼8.4 to 7.8 ka BP), the Sofular record indicates a remarkable increase in rainfall amount and intensity, in line with other paleoclimate studies in the Eastern Mediterranean. During that period, enhanced summertime insolation either produced much stronger storms in the following fall and winter through high sea surface temperatures, or it invoked a regional summer monsoon circulation and rainfall. We suggest that one or both of these climatic mechanisms led to a coupling of the BlackSea and the Mediterranean rainfall regimes at that time, which can explain the observed proxy signals. However, there are discrepancies among the Eastern Mediterranean records in terms of the timing of this wet period; implying that changes were probably not always occurring through the same mechanism. Nevertheless, the Sofular Cave record does provide hints and bring about new questions about the connection between regional and large scale climates, highlighting the need for a more extensive network of high quality paleoclimate records to better understand Holoceneclimate.

Support for global climate reorganization during the "Medieval Climate Anomaly"

Widely distributed proxy records indicate that the Medieval Climate Anomaly (MCA; *900–1350 AD) was characterized by coherent shifts in large-scale Northern Hemisphere atmospheric circulation patterns. Although cooler sea surface temperatures in the central and eastern equatorial Pacific can explain some aspects of medieval circulation changes, they are not sufficient to account for other notable features, including widespread aridity through the Eurasian sub-tropics, stronger winter westerlies across the North Atlantic and Western Europe, and shifts in monsoon rainfall patterns across Africa and South Asia. We present results from a full-physics coupled climate model showing that a slight warming of the tropical Indian and western Pacific Oceans relative to the other tropical ocean basins can induce a broad range of the medieval circulation and climate changes indicated by proxy data, including many of those not explained by a cooler tropical Pacific alone. Important aspects of the results resemble those from previous simulations examining the climatic response to the rapid Indian Ocean warming during the late twentieth century, and to results from climate warming simulations—especially in indicating an expansion of the Northern Hemisphere Hadley circulation. Notably, the pattern of tropical Indo-Pacific sea surface temperature (SST) change responsible for producing the proxy-model similarity in our results agrees well with MCA-LIA SST differences obtained in a recent proxy-based climate field reconstruction. Though much remains unclear, our results indicate that the MCA was characterized by an enhanced zonal Indo-Pacific SST gradient with resulting changes in Northern Hemisphere tropical and extra-tropical circulation patterns and hydroclimate regimes, linkages that may explain the coherent regional climate shifts indicated by proxy records from across the planet. The findings provide new perspectives on the nature and possible causes of the MCA—a remarkable, yet incompletely understood episode of Late Holocene climatic change.

Vegetation and environmental dynamics in the southern Black Sea region since 18 kyr BP derived from the marine core 22-GC3

Sediments from the Black Sea, a region historically dominated by forests and steppe landscapes, are a valuable source of detailed information on the changes in regional terrestrial and aquatic environments at decadal to millennial scales. Here we present multi-proxy environmental records (pollen, dinoflagellate cysts, Ca, Ti and oxygen isotope data) from the uppermost 305 cm of the core 22-GC3 (42°13.53′N, 36°29.55′E) collected from a water depth of 838 m in the southern part of the Black Sea in 2007. The records span the last ~ 18 kyr (all ages are given in cal kyr BP). The pollen data reveal the dominance of the Artemisia-steppe in the region, suggesting rather dry/cold environments ~ 18–14.5 kyr BP. Warming/humidity increase during melt-water pulses (~ 16.1–14.5 kyr BP), indicated by δ18O records from the 22-GC3 core sediment and from the Sofular Cave stalagmite, is expressed in more negative δ13C values from the Sofular Cave, usually interpreted as the spreading of C3 plants. The records representing the interstadial complex (~ 14.5–12.9 kyr BP) show an increase in temperature and moisture, indicated by forest development, increased primary productivity and reduced surface run-off, whereas the switch from primary terrigenous to primary authigenic Ca origin occurs ~ 500 yr later. The Younger Dryas cooling is clearly demonstrated by more negative δ13C values from the Sofular Cave and a reduction of pines. The early Holocene (11.7–8.5 kyr BP) interval reveals relatively dry conditions compared to the mostly moist and warm middle Holocene (8.5–5 kyr BP), which is characterized by the establishment of the species-rich warm mixed and temperate deciduous forests in the low elevation belt, temperate deciduous beech-hornbeam forests in the middle and cool conifer forest in upper mountain belt. The border between the early and middle Holocene in the vegetation records coincides with the opening of the Mediterranean corridor at ~ 8.3 kyr BP, as indicated by a marked change in the dinocyst assemblages and in the sediment lithology. Changes in the pollen assemblages indicate a reduction in forest cover after ~ 5 kyr BP, which was likely caused by increased anthropogenic pressure on the regional vegetation.

Progressive palaeoenvironmental change during the Late Barremian-Early Aptian as prelude to Oceanic Anoxic Event 1a: Evidence from the Gorgo a Cerbara section (Umbria-Marche basin, central Italy)

During Oceanic Anoxic Event 1a (OAE 1a, 120 Ma; Li et al., 2008), organic carbon-rich layers were deposited in marine environments under anoxic conditions on a global scale. In this study, palaeoenvironmental conditions leading to this event are characterised by studying the Upper Barremian to the Lower Aptian succession of the Gorgo a Cerbara section (central Italy). For this, an integrated multi-proxy approach (δ13Ccarb; δ13Corg; δ18O; phosphorus; Total Organic Carbon, TOC; bulk-rock mineralogy, as well as redox-sensitive trace elements — RSTEs) has been applied. During the LateBarremian, thin organic-rich layers occur episodically, and associated Corg:Ptot ratios indicate the presence of intermittent dysoxic to anoxic conditions. Coarse correlations are observed between TOC, P and biogenic silica contents, indicating links between P availability, productivity, and TOC preservation. However, the corresponding δ13Ccarb and δ18O records remain quite stable, indicating that these brief periods of enhanced TOC preservation did not have sufficient impact on the marine carbon reservoir to deviate δ13C records. Around the Barremian–Aptian boundary, TOC-enriched layers become more frequent. These layers correlate with negative excursions in the δ13Ccarb and δ13Corg records, possibly due to a warming period as indicated by the δ18O record. During the earliest Aptian, this warming trend is reverted into a cooling trend, which is then followed by an important warming step near the onset of Oceanic Anoxic Event 1a (OAE 1a). During this time period, organic-rich intervals occur, which are characterised by the progressive increase in RSTE. The warming step prior the onset of OAE 1a is associated with the well-known negative spike in δ13Ccarb and δ13Corg records, an important peak in P accumulation, RSTE enrichments and Corg:Ptot ratios indicating the prevalence of anoxic conditions. The Selli Level itself may document a cooling phase. RSTE enrichments and Corg:Ptot ratios confirm the importance of anoxic conditions during OAE 1a at this site. The Gorgo a Cerbara section is interpreted to reflect the progressive impact of palaeoenvironmental change related to the formation of the Ontong-Java plate-basalt plateau, which started already around the Barremian–Aptian boundary and culminated into OAE 1a.

Oxygen isotope ratios in the shell of Mytilus edulis: archives of glacier meltwater in Greenland?

Melting of the Greenland Ice Sheet (GrIS) is accelerating and will contribute significantly to global sea level rise during the 21st century. Instrumental data on GrIS melting only cover the last few decades, and proxy data extending our knowledge into the past are vital for validating models predicting the influence of ongoing climate change. We investigated a potential meltwater proxy in Godthåbsfjord (West Greenland), where glacier meltwater causes seasonal excursions with lower oxygen isotope water (δ18Ow) values and salinity. The blue mussel (Mytilus edulis) potentially records these variations, because it precipitates its shell calcite in oxygen isotopic equilibrium with ambient seawater. As M. edulis shells are known to occur in raised shorelines and archaeological shell middens from previous Holocene warm periods, this species may be ideal in reconstructing past meltwater dynamics. We investigate its potential as a palaeo-meltwater proxy. First, we confirmed that M. edulis shell calcite oxygen isotope (δ18Oc) values are in equilibrium with ambient water and generally reflect meltwater conditions. Subsequently we investigated if this species recorded the full range of δ18Ow values occurring during the years 2007 to 2010. Results show that δ18Ow values were not recorded at very low salinities (< ~ 19), because the mussels appear to cease growing. This implies that Mytilus edulis δ18Oc values are suitable in reconstructing past meltwater amounts in most cases, but care has to be taken that shells are collected not too close to a glacier, but rather in the mid-region or mouth of the fjord. The focus of future research will expand on the geographical and temporal range of the shell measurements by sampling mussels in other fjords in Greenland along a south–north gradient, and by sampling shells from raised shorelines and archaeological shell middens from prehistoric settlements in Greenland.

Evidence for global climate reorganization during medieval times

A synthesis of global climate model results and inferences from proxy records suggests an increased sea surface temperature gradient between the tropical Indian and Pacific Oceans during medieval times.

Accurate analysis of noble gas concentrations in small water samples and its application to fluid inclusions in stalagmites

The concentrations of dissolved noble gases in water are widely used as a climate proxy to determine noble gas temperatures (NGTs); i.e., the temperature of the water when gas exchange last occurred. In this paper we make a step forward to apply this principle to fluid inclusions in stalagmites in order to reconstruct the cave temperature prevailing at the time when the inclusion was formed. We present an analytical protocol that allows us accurately to determine noble gas concentrations and isotope ratios in stalagmites, and which includes a precise manometrical determination of the mass of water liberated from fluid inclusions. Most important for NGT determination is to reduce the amount of noble gases liberated from air inclusions, as they mask the temperature-dependent noble gas signal from the water inclusions. We demonstrate that offline pre-crushing in air to subsequently extract noble gases and water from the samples by heating is appropriate to separate gases released from air and water inclusions. Although a large fraction of recent samples analysed by this technique yields NGTs close to present-day cave temperatures, the interpretation of measured noble gas concentrations in terms of NGTs is not yet feasible using the available least squares fitting models. This is because the noble gas concentrations in stalagmites are not only composed of the two components air and air saturated water (ASW), which these models are able to account for. The observed enrichments in heavy noble gases are interpreted as being due to adsorption during sample preparation in air, whereas the excess in He and Ne is interpreted as an additional noble gas component that is bound in voids in the crystallographic structure of the calcite crystals. As a consequence of our study's findings, NGTs will have to be determined in the future using the concentrations of Ar, Kr and Xe only. This needs to be achieved by further optimizing the sample preparation to minimize atmospheric contamination and to further reduce the amount of noble gases released from air inclusions.

Physically consistent responses of the global atmospheric hydrological cycle in models and observations

Robust and physically understandable responses of the global atmospheric water cycle to a warming climate are presented. By considering interannual responses to changes in surface temperature (T), observations and AMIP5 simulations agree on an increase in column integrated water vapor at the rate 7 %/K (in line with the Clausius­Clapeyron equation) and of precipitation at the rate 2-­3 %/K (in line with energetic constraints). Using simple and complex climate models, we demonstrate that radiative forcing by greenhouse gases is currently suppressing global precipitation (P) at ~ -0.15 %/decade. Along with natural variability, this can explain why observed trends in global P over the period 1988-2008 are close to zero. Regional responses in the global water cycle are strongly constrained by changes in moisture fluxes. Model simulations show an increased moisture flux into the tropical wet region at 900 hPa and an enhanced outflow (of smaller magnitude) at around 600 hPa with warming. Moisture transport explains an increase in P in the wet tropical regions and small or negative changes in the dry regions of the subtropics in CMIP5 simulations of a warming climate. For AMIP5 simulations and satellite observations, the heaviest 5-day rainfall totals increase in intensity at ~15 %/K over the ocean with reductions at all percentiles over land. The climate change response in CMIP5 simulations shows consistent increases in P over ocean and land for the highest intensities, close to the Clausius-Clapeyron scaling of 7 %/K, while P declines for the lowest percentiles, indicating that interannual variability over land may not be a good proxy for climate change. The local changes in precipitation and its extremes are highly dependent upon small shifts in the large-scale atmospheric circulation and regional feedbacks.

Response of the northern stratospheric polar vortex to the seasonal alignment of QBO phase transitions

This study considers the strength of the Northern Hemisphere Holton-Tan effect (HTE) in terms of the phase alignment of the quasi-biennial oscillation (QBO) with respect to the annual cycle. Using the ERA-40 Reanalysis, it is found that the early winter (Nov–Dec) and late winter (Feb–Mar) relation between QBO phase and the strength of the stratospheric polar vortex is optimized for subsets of the 44-year record that are chosen on the basis of the seasonality of QBO phase transitions at the 30 hPa level. The timing of phase transitions serves as a proxy for changes in the vertical structure of the QBO over the whole depth of the tropical stratosphere. The statistical significance of the Nov–Dec (Feb–Mar) HTE is greatest when 30 hPa QBO phase transitions occur 9–14 (4–9) months prior to the January of the NH winter in question. This suggests that there exists for both early and late winter a vertical structure of tropical stratospheric winds that is most effective at influencing the interannual variability of the polar vortex, and that an early (late) winter HTE is associated with an early (late) progression of QBO phase towards that structure. It is also shown that the seasonality of QBO phase transitions at 30 hPa varies on a decadal timescale, with transitions during the first half of the calendar year being relatively more common during the first half of the tropical radiosonde wind record. Combining these two results suggests that decadal changes in HTE strength could result from the changing seasonality of QBO phase transitions. Citation: Anstey, J. A., and T. G. Shepherd (2008), Response of the northern stratospheric polar vortex to the seasonal alignment of QBO phase transitions, Geophys. Res. Lett., 35, L22810, doi:10.1029/2008GL035721.

Influence of the quasi-biennial oscillation on the extratropical winter stratosphere in an atmospheric general circulation model and in reanalysis data

The interannual variability of the stratospheric polar vortex during winter in both hemispheres is observed to correlate strongly with the phase of the quasi-biennial oscillation (QBO) in tropical stratospheric winds. It follows that the lack of a spontaneously generated QBO in most atmospheric general circulation models (AGCMs) adversely affects the nature of polar variability in such models. This study examines QBO–vortex coupling in an AGCM in which a QBO is spontaneously induced by resolved and parameterized waves. The QBO–vortex coupling in the AGCM compares favorably to that seen in reanalysis data [from the 40-yr ECMWF Re-Analysis (ERA-40)], provided that careful attention is given to the definition of QBO phase. A phase angle representation of the QBO is employed that is based on the two leading empirical orthogonal functions of equatorial zonal wind vertical profiles. This yields a QBO phase that serves as a proxy for the vertical structure of equatorial winds over the whole depth of the stratosphere and thus provides a means of subsampling the data to select QBO phases with similar vertical profiles of equatorial zonal wind. Using this subsampling, it is found that the QBO phase that induces the strongest polar vortex response in early winter differs from that which induces the strongest late-winter vortex response. This is true in both hemispheres and for both the AGCM and ERA-40. It follows that the strength and timing of QBO influence on the vortex may be affected by the partial seasonal synchronization of QBO phase transitions that occurs both in observations and in the model. This provides a mechanism by which changes in the strength of QBO–vortex correlations may exhibit variability on decadal time scales. In the model, such behavior occurs in the absence of external forcings or interannual variations in sea surface temperatures.

Stable isotope and modelling evidence for CO2 as a driver of glacial–interglacial vegetation shifts in southern Africa

Atmospheric CO2 concentration is hypothesized to influence vegetation distribution via tree–grass competition, with higher CO2 concentrations favouring trees. The stable carbon isotope (δ13C) signature of vegetation is influenced by the relative importance of C4 plants (including most tropical grasses) and C3 plants (including nearly all trees), and the degree of stomatal closure – a response to aridity – in C3 plants. Compound-specific δ13C analyses of leaf-wax biomarkers in sediment cores of an offshore South Atlantic transect are used here as a record of vegetation changes in subequatorial Africa. These data suggest a large increase in C3 relative to C4 plant dominance after the Last Glacial Maximum. Using a process-based biogeography model that explicitly simulates 13C discrimination, it is shown that precipitation and temperature changes cannot explain the observed shift in δ13C values. The physiological effect of increasing CO2 concentration is decisive, altering the C3/C4 balance and bringing the simulated and observed δ13C values into line. It is concluded that CO2 concentration itself was a key agent of vegetation change in tropical southern Africa during the last glacial–interglacial transition. Two additional inferences follow. First, long-term variations in terrestrial δ13Cvalues are not simply a proxy for regional rainfall, as has sometimes been assumed. Although precipitation and temperature changes have had major effects on vegetation in many regions of the world during the period between the Last Glacial Maximum and recent times, CO2 effects must also be taken into account, especially when reconstructing changes in climate between glacial and interglacial states. Second, rising CO2 concentration today is likely to be influencing tree–grass competition in a similar way, and thus contributing to the "woody thickening" observed in savannas worldwide. This second inference points to the importance of experiments to determine how vegetation composition in savannas is likely to be influenced by the continuing rise of CO2 concentration.

A high resolution late Holocene palaeo environmental record from the central Adriatic Sea

A multi-proxy study of a Holocene sediment core (RF 93-30) from the western flank of the central Adriatic, in 77 m of water, reveals a sequence of changes in terrestrial vegetation, terrigenous sediment input and benthic fauna, as well as evidence for variations in sea surface temperature spanning most of the last 7000 yr. The chronology of sedimentation is based on several lines of evidence, including AMS 14C dates of foraminifera extracted from the core, palaeomagnetic secular variation, pollen indicators and dated tephra. The temporal resolution increases towards the surface and, for some of the properties measured, is sub-decadal for the last few centuries. The main changes recorded in vegetation, sedimentation and benthic foraminiferal assemblages appear to be directly related to human activity in the sediment source area, which includes the Po valley and the eastern flanks of the central and northern Appenines. The most striking episodes of deforestation and expanding human impact begin around 3600 BP (Late Bronze Age) and 700 BP (Medieval) and each leads to an acceleration in mass sedimentation and an increase in the proportion of terrigenous material, reflecting the response of surface processes to widespread forest clearance and cultivation. Although human impact appears to be the proximal cause of these changes, climatic effects may also have been important. During these periods, signs of stress are detectable in the benthic foram morphotype assemblages. Between these two periods of increased terrigeneous sedimentation there is smaller peak in sedimentation rate around 2400BP which is not associated with evidence for deforestation, shifts in the balance between terrigenous and authigenic sedimentation, or changes in benthic foraminifera. The mineral magnetic record provides a sensitive indicator of changing sediment sources: during forested periods of reduced terrigenous input it is dominated by authigenic bacterial magnetite, whereas during periods of increased erosion, anti-ferromagetic minerals (haematite and/or goethite) become more important, as well as both paramagnetic minerals and super-paramagnetic magnetite. Analysis of the alkenone, U37k′, record provides an indication of possible changes in sea surface temperature during the period, but it is premature to place too much reliance on these inferred changes until the indirect effects of past changes in the depth of the halocline and in circulation have been more fully evaluated. The combination of methods used and the results obtained illustrate the potential value of such high resolution near-shore marine sedimentary sequences for recording wide-scale human impact, documenting the effects of this on marine sedimentation and fauna and, potentially, disentangling evidence for human activities from that for past changes in climate.

Simulating the interannual variability of major dust storms on Mars using variable lifting thresholds

The redistribution of a finite amount of martian surface dust during global dust storms and in the intervening periods has been modelled in a dust lifting version of the UK Mars General Circulation Model. When using a constant, uniform threshold in the model’s wind stress lifting parameterisation and assuming an unlimited supply of surface dust, multiannual simulations displayed some variability in dust lifting activity from year to year, arising from internal variability manifested in surface wind stress, but dust storms were limited in size and formed within a relatively short seasonal window. Lifting thresholds were then allowed to vary at each model gridpoint, dependent on the rates of emission or deposition of dust. This enhanced interannual variability in dust storm magnitude and timing, such that model storms covered most of the observed ranges in size and initiation date within a single multiannual simulation. Peak storm magnitude in a given year was primarily determined by the availability of surface dust at a number of key sites in the southern hemisphere. The observed global dust storm (GDS) frequency of roughly one in every 3 years was approximately reproduced, but the model failed to generate these GDSs spontaneously in the southern hemisphere, where they have typically been observed to initiate. After several years of simulation, the surface threshold field—a proxy for net change in surface dust density—showed good qualitative agreement with the observed pattern of martian surface dust cover. The model produced a net northward cross-equatorial dust mass flux, which necessitated the addition of an artificial threshold decrease rate in order to allow the continued generation of dust storms over the course of a multiannual simulation. At standard model resolution, for the southward mass flux due to cross-equatorial flushing storms to offset the northward flux due to GDSs on a timescale of ∼3 years would require an increase in the former by a factor of 3–4. Results at higher model resolution and uncertainties in dust vertical profiles mean that quasi-periodic redistribution of dust on such a timescale nevertheless appears to be a plausible explanation for the observed GDS frequency.

Past and recent changes in the North Atlantic oscillation

The North Atlantic oscillation (NAO) is under current climate conditions the leading mode of atmospheric circulation variability over the North Atlantic region. While the pattern is present during the entire year, it is most important during winter, explaining a large part of the variability of the large-scale pressure field, being thus largely determinant for the weather conditions over the North Atlantic basin and over Western Europe. In this study, a review of recent literature on the basic understanding of the NAO, its variability on different time scales and driving physical mechanisms is presented. In particular, the observed NAO variations and long-term trends are put into a long term perspective by considering paleo-proxy evidence. A representative number of recently released NAO reconstructions are discussed. While the reconstructions agree reasonably well with observations during the instrumental overlapping period, there is a rather high uncertainty between the different reconstructions for the pre-instrumental period, which leads to partially incoherent results, that is, periods where the NAO reconstructions do not agree even in sign. Finally, we highlight the future need of a broader definition of the NAO, the assessment of the stability of the teleconnection centers over time, the analysis of the relations to other relevant variables like temperature and precipitation, as well as on the relevant processes involved