Links between tropical Pacific seasonal, interannual and orbital variability during the Holocene

The El Niño/Southern Oscillation (ENSO) is the leading mode of interannual climate variability. However, it is unclear how ENSO has responded to external forcing, particularly orbitally induced changes in the amplitude of the seasonal cycle during the Holocene. Here we present a reconstruction of seasonal and interannual surface conditions in the tropical Pacific Ocean from a network of high-resolution coral and mollusc records that span discrete intervals of the Holocene. We identify several intervals of reduced variance in the 2 to 7 yr ENSO band that are not in phase with orbital changes in equatorial insolation, with a notable 64% reduction between 5,000 and 3,000 years ago. We compare the reconstructed ENSO variance and seasonal cycle with that simulated by nine climate models that include orbital forcing, and find that the models do not capture the timing or amplitude of ENSO variability, nor the mid-Holocene increase in seasonality seen in the observations; moreover, a simulated inverse relationship between the amplitude of the seasonal cycle and ENSO-related variance in sea surface temperatures is not found in our reconstructions. We conclude that the tropical Pacific climate is highly variable and subject to millennial scale quiescent periods. These periods harbour no simple link to orbital forcing, and are not adequately simulated by the current generation of models.

Holocene lake-level trends in the Rocky Mountains, USA

The availability of water shapes life in the western United States, and much of the water in the region originates in the Rocky Mountains. Few studies, however, have explicitly examined the history of water levels in the Rocky Mountains during the Holocene. Here, we examine the past levels of three lakes near the Continental Divide in Montana and Colorado to reconstruct Holocene moisture trends. Using transects of sediment cores and sub-surface geophysical profiles from each lake, we find that mid-Holocene shorelines in the small lakes (4–110 ha) were as much as ~10 m below the modern lake surfaces. Our results are consistent with existing evidence from other lakes and show that a wide range of settings in the region were much drier than today before 3000–2000 years ago. We also discuss evidence for millennial-scale moisture variation, including an abruptly-initiated and -terminated wet period in Colorado from 4400 to 3700 cal yr BP, and find only limited evidence for low-lake stands during the past millennium. The extent of low-water levels during the mid-Holocene, which were most severe and widespread ca. 7000–4500 cal yr BP, is consistent with the extent of insolation-induced aridity in previously published regional climate model simulations. Like the simulations, the lake data provide no evidence for enhanced zonal flow during the mid-Holocene, which has been invoked to explain enhanced mid-continent aridity at the time. The data, including widespread evidence for large changes on orbital time scales and for more limited changes during the last millennium, confirm the ability of large boundary-condition changes to push western water supplies beyond the range of recent natural variability.

The Global Paleomonsoon as seen through speleothem records from Asia and the Americas

The regional monsoons of the world have long been viewed as seasonal atmospheric circulation reversal-analogous to a thermally-driven land-sea breeze on a continental scale. This conventional view of monsoons is now being integrated at a global scale and accordingly, a new paradigm has emerged which considers regional monsoons to be manifestations of global-scale seasonal changes in response to overturning of atmospheric circulation in the tropics and subtropics, and henceforth, interactive components of a singular Global Monsoon (GM) system. The paleoclimate community, however, tends to view 'paleomonsoon' (PM), largely in terms of regional circulation phenomena. In the past decade, many high-quality speleothem oxygen isotope (delta O-18) records have been established from the Asian Monsoon and the South American Monsoon regions that primarily reflect changes in the integrated intensities of monsoons on orbital-to-decadal timescales. With the emergence of these high-resolution and absolute-dated records from both sides of the Equator, it is now possible to test a concept of the 'Global-Paleo-Monsoon' (GPM) on a wide-range of timescales. Here we present a comprehensive synthesis of globally-distributed speleothem delta O-18 records and highlight three aspects of the GPM that are comparable to the modern GM: (1) the GPM intensity swings on different timescales; (2) their global extent; and (3) an anti-phased inter-hemispheric relationship between the Asian and South American monsoon systems on a wide range of timescales.

Evidenze micropaleontologiche e sedimentologiche di cicli deposizionali e climatici circa millenari nei Depositi Tardoquaternari della Piana dell'Arno e del Delta del Po

A multidisciplinary study was carried out on the Late Quaternary-Holocene subsurface deposits of two Mediterranean coastal areas: Arno coastal plain (Northern Tyrrhenian Sea) and Modern Po Delta (Northern Adriatic Sea). Detailed facies analyses, including sedimentological and micropalaeontological (benthic foraminifers and ostracods) investigations, were performed on nine continuously-cored boreholes of variable depth (ca. from 30 meters to100 meters). Six cores were located in the Arno coastal plain and three cores in the Modern Po Delta. To provide an accurate chronological framework, twenty-four organic-rich samples were collected along the fossil successions for radiocarbon dating (AMS 14C). In order to reconstruct the depositional and palaeoenvironmental evolution of the study areas, core data were combined with selected well logs, provided by local companies, along several stratigraphic sections. These sections revealed the presence of a transgressive-regressive (T-R) sequence, composing of continental, coastal and shallow-marine deposits dated to the Late Pleistocene-Holocene period, beneath the Arno coastal plain and the Modern Po Delta. Above the alluvial deposits attributed to the last glacial period, the post-glacial transgressive succession (TST) consists of back-barrier, transgressive barrier and inner shelf deposits. Peak of transgression (MFS) took place around the Late-Middle Holocene transition and was identified by subtle micropalaeontological indicators within undifferentiated fine-grained deposits. Upward a thick prograding succession (HST) records the turnaround to regressive conditions that led to a rapid delta progradation in both study areas. Particularly, the outbuilding of modern-age Po Delta coincides with mud-belt formation during the late HST (ca. 600 cal yr BP), as evidenced by a fossil microfauna similar to the foraminiferal assemblage observed in the present Northern Adriatic mud-belt. A complex interaction between allocyclic and autocyclic factors controlled facies evolution during the highstand period. The presence of local parameters and the absence of a predominant factor prevent from discerning or quantifying consequences of the complex relationships between climate and deltaic evolution. On the contrary transgressive sedimentation seems to be mainly controlled by two allocyclic key factors, sea-level rise and climate variability, that minimized the effects of local parameters on coastal palaeoenvironments. TST depositional architecture recorded in both study areas reflects a well-known millennial-scale variability of sea-level rising trend and climate during the Late glacial-Holocene period. Repeated phases of backswamp development and infilling by crevasse processes (parasequences) were recorded in the subsurface of Modern Po Delta during the early stages of transgression (ca. 11,000-9,500 cal yr BP). In the Arno coastal plain the presence of a deep-incised valley system, probably formed at OSI 3/2 transition, led to the development of a thick (ca. 35-40 m) transgressive succession composed of coastal plain, bay-head delta and estuarine deposits dated to the Last glacial-Early Holocene period. Within the transgressive valley fill sequence, high-resolution facies analyses allowed the identification and lateral tracing of three parasequences of millennial duration. The parasequences, ca. 8-12 meters thick, are bounded by flooding surfaces and show a typical internal shallowing-upward trend evidenced by subtle micropalaeontological investigations. The vertical stacking pattern of parasequences shows a close affinity with the step-like sea-level rising trend occurred between 14,000-8,000 cal years BP. Episodes of rapid sea-level rise and subsequent stillstand phases were paralleled by changes in climatic conditions, as suggested by pollen analyses performed on a core drilled in the proximal section of the Arno palaeovalley (pollen analyses performed by Dr. Marianna Ricci Lucchi). Rapid shifts to warmer climate conditions accompanied episodes of rapid sea-level rise, in contrast stillstand phases occurred during temporary colder climate conditions. For the first time the palaeoclimatic signature of high frequency depositional cycles is clearly documented. Moreover, two of the three "regressive" pulsations, recorded at the top of parasequences by episodes of partial estuary infilling in the proximal and central portions of Arno palaeovalley, may be correlated with the most important cold events of the post-glacial period: Younger Dryas and 8,200 cal yr BP event. The stratigraphic and palaeoclimatic data of Arno coastal plain and Po Delta were compared with those reported for the most important deltaic and coastal systems in the worldwide literature. The depositional architecture of transgressive successions reflects the strong influence of millennial-scale eustatic and climatic variability on worldwide coastal sedimentation during the Late glacial-Holocene period (ca. 14,000-7,000 cal yr BP). The most complete and accurate record of high-frequency eustatic and climatic events are usually found within the transgressive succession of very high accommodation settings, such as incised-valley systems where exceptionally thick packages of Late glacial-Early Holocene deposits are preserved.

Paleosol-based stratigraphy of late Quaternary deposits from the Po basin, between Bologna and the Po river (Northern Italy)

Stratigraphic studies carried out over the last decades in Italy and elsewhere testify a growing interest in Quaternary deposits and in the influence of climate change on their architecture. The subsurface of the Po Plain, in its topmost portion, is made up of alluvial deposits organized in depositional cycles at different scales. This PhD thesis provides millennial-scale stratigraphic reconstruction of the Late Pleistocene-Holocene deposits beneath the southern Po Plain, based on basin-scale correlation of laterally-extensive buried soil horizons. Far from the aim of characterizing palaeosols from a mineralogical and geochemical point of view, we focused on the physical and stratigraphic significance of these horizons. In the Bologna urban area, which hosts an abundance of stratigraphic data, the correlation between seventeen continuously-cored boreholes led to the identification of five vertically-stacked palaeosol-bounded sequences within the 14C time window. In a wide portion of the alluvial plain north of Bologna, far away from the Apenninic margin and from the Po River, where subsurface stratigraphic architecture is dominated by markedly lenticular sediment bodies, palaeosols revealed to be the only stratigraphic marker of remarkable lateral continuity. These horizons are characterized by peculiar resistance values, which make them easily identifiable via pocket penetration tests. Palaeosols reveal specific geometric relationships with the associated alluvial facies associations, allowing reliable estimates of soil development as a function of alluvial dynamics. With the aid of sixty new radiocarbon dates, a reliable age attribution and likely time intervals of exposure were assigned to each palaeosol. Vertically-stacked palaeosols delimitate short-term depositional cycles, likely related to the major episodes of climatic change of the last 40 ky. Through integration of stratigraphic data with 750 archaeological reports from the Bologna area, the impact of human settlements on depositional and pedogenic processes during the late Holocene was investigated.

Impact of Holocene climate changes on alpine and treeline vegetation at Sanetsch Pass, Bernese Alps, Switzerland

In order to infer reactions of treeline and alpine vegetation to climatic change, past vegetation changes are reconstructed on the basis of pollen, macrofossil and charcoal analysis. The sampled sediment cores originate from the small pond Emines, located at the Sanetsch Pass (connecting the Valais and Bern, Switzerland) at an altitude of 2288 m a.s.l. Today's treeline is at ca. 2200 m a.s.l. in the area, though due to special pass (saddle) conditions it is locally depressed to ca. 2060 m a.s.l. Our results reveal that the area around Emines was covered by treeless alpine vegetation during most of the past 12,000 years. Single individuals of Betula, Larix decidua and possibly Pinus cembra occurred during the Holocene. Major centennial to millennial-scale responses of treeline vegetation to climatic changes are evident. However, alpine vegetation composition remained rather stable between 11,500 and 6000 cal. BP, showing that Holocene climatic changes of +/− 1 °C hardly influenced the local vegetation at Emines. The rapid warming of 3–4 °C at the Late Glacial/Holocene transition (11,600 cal. BP) caused significant altitudinal displacements of alpine species that were additionally affected by the rapid upward movement of trees and shrubs. Since the beginning of the Neolithic, vegetation changes at Sanetsch Pass resulted from a combination of climate change and human impact. Anthropogenic fire increase and land-use change combined with a natural change from subcontinental to more oceanic climate during the second half of the Holocene led to the disappearance of P. cembra in the study area, but favoured the occurrence of Picea abies and Alnus viridis. The mid- to late-Holocene decline of Abies alba was primarily a consequence of human impact, since this mesic species should have benefitted from a shift to more oceanic conditions. Future alpine vegetation changes will be a function of the amplitude and rapidity of global warming as well as human land use. Our results imply that alpine vegetation at our treeline pass site was never replaced by forests since the last ice-age. This may change in the future if anticipated climate change will induce upslope migration of trees. The results of this study emphasise the necessity of climate change mitigation in order to prevent biodiversity losses as a consequence of unprecedented community and species displacement in response to climatic change.

Impact of suborbital climate changes in the North Atlantic on ice sheet dynamics at the Mid-Pleistocene Transition

[1] Early and Mid-Pleistocene climate, ocean hydrography and ice sheet dynamics have been reconstructed using a high-resolution data set (planktonic and benthicδ18O time series, faunal-based sea surface temperature (SST) reconstructions and ice-rafted debris (IRD)) record from a high-deposition-rate sedimentary succession recovered at the Gardar Drift formation in the subpolar North Atlantic (Integrated Ocean Drilling Program Leg 306, Site U1314). Our sedimentary record spans from late in Marine Isotope Stage (MIS) 31 to MIS 19 (1069–779 ka). Different trends of the benthic and planktonic oxygen isotopes, SST and IRD records before and after MIS 25 (∼940 ka) evidence the large increase in Northern Hemisphere ice-volume, linked to the cyclicity change from the 41-kyr to the 100-kyr that occurred during the Mid-Pleistocene Transition (MPT). Beside longer glacial-interglacial (G-IG) variability, millennial-scale fluctuations were a pervasive feature across our study. Negative excursions in the benthicδ18O time series observed at the times of IRD events may be related to glacio-eustatic changes due to ice sheets retreats and/or to changes in deep hydrography. Time series analysis on surface water proxies (IRD, SST and planktonicδ18O) of the interval between MIS 31 to MIS 26 shows that the timing of these millennial-scale climate changes are related to half-precessional (10 kyr) components of the insolation forcing, which are interpreted as cross-equatorial heat transport toward high latitudes during both equinox insolation maxima at the equator.

Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

Correct estimation of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice core studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: outputs of a firn densification model, and measurements of δ15N of N2 in air trapped in ice core, assuming that δ15N is only affected by gravitational fractionation in the firn column. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi-coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available ice core air-δ15N measurements from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rates and temperature conditions. Our δ15N profiles reveal a heterogeneous response of the firn structure to glacial–interglacial climatic changes. While firn densification simulations correctly predict TALDICE δ15N variations, they systematically fail to capture the large millennial-scale δ15N variations measured at BI and the δ15N glacial levels measured at JRI and EDML – a mismatch previously reported for central East Antarctic ice cores. New constraints of the EDML gas–ice depth offset during the Laschamp event (~41 ka) and the last deglaciation do not favour the hypothesis of a large convective zone within the firn as the explanation of the glacial firn model–δ15N data mismatch for this site. While we could not conduct an in-depth study of the influence of impurities in snow for firnification from the existing datasets, our detailed comparison between the δ15N profiles and firn model simulations under different temperature and accumulation rate scenarios suggests that the role of accumulation rate may have been underestimated in the current description of firnification models.

Estimated strength of the Atlantic overturning circulation during the last deglaciation

The Atlantic meridional overturning circulation affects the latitudinal distribution of heat, and is a key component of the climate system. Proxy reconstructions, based on sedimentary Pa-231/Th-230 ratios and the difference between surface-and deep-water radiocarbon ages, indicate that during the last glacial period, the overturning circulation was reduced during millennial-scale periods of cooling(1-5). However, much debate exists over the robustness of these proxies(6-8). Here e combine proxy reconstructions of sea surface and air temperatures and a global climate model to quantitatively estimate changes in the strength of the Atlantic meridional overturning circulation during the last glacial period. We find that, relative to the Last Glacial Maximum, the overturning circulation was reduced by approximately 14 Sv during the cold Heinrich event 1. During the Younger Dryas cold event, the overturning circulation was reduced by approximately 12 Sv, relative to the preceding warm interval. These changes are consistent with qualitative estimates of the overturning circulation from sedimentary Pa-231/Th-230 ratios. In addition, we find that the strength of the overturning circulation during the Last Glacial Maximum and the Holocene epoch are indistinguishable within the uncertainty of the reconstruction.

High-resolution glacial and deglacial record of atmospheric methane by continuous-flow and laser spectrometer analysis along the NEEM ice core

The Greenland NEEM (North Greenland Eemian Ice Drilling) operation in 2010 provided the first opportunity to combine trace-gas measurements by laser spectroscopic instruments and continuous-flow analysis along a freshly drilled ice core in a field-based setting. We present the resulting atmospheric methane (CH4) record covering the time period from 107.7 to 9.5 ka b2k (thousand years before 2000 AD). Companion discrete CH4 measurements are required to transfer the laser spectroscopic data from a relative to an absolute scale. However, even on a relative scale, the high-resolution CH4 data set significantly improves our knowledge of past atmospheric methane concentration changes. New significant sub-millennial-scale features appear during interstadials and stadials, generally associated with similar changes in water isotopic ratios of the ice, a proxy for local temperature. In addition to the midpoint of Dansgaard–Oeschger (D/O) CH4 transitions usually used for cross-dating, sharp definition of the start and end of these events brings precise depth markers (with ±20 cm uncertainty) for further cross-dating with other palaeo- or ice core records, e.g. speleothems. The method also provides an estimate of CH4 rates of change. The onsets of D/O events in the methane signal show a more rapid rate of change than their endings. The rate of CH4 increase associated with the onsets of D/O events progressively declines from 1.7 to 0.6 ppbv yr−1 in the course of marine isotope stage 3. The largest observed rate of increase takes place at the onset of D/O event #21 and reaches 2.5 ppbv yr−1.

Spatial gradients of temperature, accumulation and delta18O-ice in Greenland over a series of Dansgaard-Oeschger events

Air and water stable isotope measurements from four Greenland deep ice cores (GRIP, GISP2, NGRIP and NEEM) are investigated over a series of Dansgaard–Oeschger events (DO 8, 9 and 10), which are representative of glacial millennial scale variability. Combined with firn modeling, air isotope data allow us to quantify abrupt temperature increases for each drill site (1σ = 0.6 °C for NEEM, GRIP and GISP2, 1.5 °C for NGRIP). Our data show that the magnitude of stadial–interstadial temperature increase is up to 2 °C larger in central and North Greenland than in northwest Greenland: i.e., for DO 8, a magnitude of +8.8 °C is inferred, which is significantly smaller than the +11.1 °C inferred at GISP2. The same spatial pattern is seen for accumulation increases. This pattern is coherent with climate simulations in response to reduced sea-ice extent in the Nordic seas. The temporal water isotope (δ18O)–temperature relationship varies between 0.3 and 0.6 (±0.08) ‰ °C−1 and is systematically larger at NEEM, possibly due to limited changes in precipitation seasonality compared to GISP2, GRIP or NGRIP. The gas age−ice age difference of warming events represented in water and air isotopes can only be modeled when assuming a 26% (NGRIP) to 40% (GRIP) lower accumulation than that derived from a Dansgaard–Johnsen ice flow model.

State of the Antarctic and Southern Ocean Climate System

This paper reviews developments in our understanding of the state of the Antarctic and Southern Ocean climate and its relation to the global climate system over the last few millennia. Climate over this and earlier periods has not been stable, as evidenced by the occurrence of abrupt changes in atmospheric circulation and temperature recorded in Antarctic ice core proxies for past climate. Two of the most prominent abrupt climate change events are characterized by intensification of the circumpolar westerlies (also known as the Southern Annular Mode) between similar to 6000 and 5000 years ago and since 1200-1000 years ago. Following the last of these is a period of major trans-Antarctic reorganization of atmospheric circulation and temperature between A. D. 1700 and 1850. The two earlier Antarctic abrupt climate change events appear linked to but predate by several centuries even more abrupt climate change in the North Atlantic, and the end of the more recent event is coincident with reorganization of atmospheric circulation in the North Pacific. Improved understanding of such events and of the associations between abrupt climate change events recorded in both hemispheres is critical to predicting the impact and timing of future abrupt climate change events potentially forced by anthropogenic changes in greenhouse gases and aerosols. Special attention is given to the climate of the past 200 years, which was recorded by a network of recently available shallow firn cores, and to that of the past 50 years, which was monitored by the continuous instrumental record. Significant regional climate changes have taken place in the Antarctic during the past 50 years. Atmospheric temperatures have increased markedly over the Antarctic Peninsula, linked to nearby ocean warming and intensification of the circumpolar westerlies. Glaciers are retreating on the peninsula, in Patagonia, on the sub-Antarctic islands, and in West Antarctica adjacent to the peninsula. The penetration of marine air masses has become more pronounced over parts of West Antarctica. Above the surface, the Antarctic troposphere has warmed during winter while the stratosphere has cooled year-round. The upper kilometer of the circumpolar Southern Ocean has warmed, Antarctic Bottom Water across a wide sector off East Antarctica has freshened, and the densest bottom water in the Weddell Sea has warmed. In contrast to these regional climate changes, over most of Antarctica, near-surface temperature and snowfall have not increased significantly during at least the past 50 years, and proxy data suggest that the atmospheric circulation over the interior has remained in a similar state for at least the past 200 years. Furthermore, the total sea ice cover around Antarctica has exhibited no significant overall change since reliable satellite monitoring began in the late 1970s, despite large but compensating regional changes. The inhomogeneity of Antarctic climate in space and time implies that recent Antarctic climate changes are due on the one hand to a combination of strong multidecadal variability and anthropogenic effects and, as demonstrated by the paleoclimate record, on the other hand to multidecadal to millennial scale and longer natural variability forced through changes in orbital insolation, greenhouse gases, solar variability, ice dynamics, and aerosols. Model projections suggest that over the 21st century the Antarctic interior will warm by 3.4 degrees +/- 1 degrees C, and sea ice extent will decrease by similar to 30%. Ice sheet models are not yet adequate enough to answer pressing questins about the effect of projected warming on mass balance and sea level. Considering the potentially major impacts of a warming climate on Antarctica, vigorous efforts are needed to better understand all aspects of the highly coupled Antarctic climate system as well as its influence on the Earth's climate and oceans.

Iron Fertilization of the Subantarctic Ocean During the Last Ice Age

John H. Martin, who discovered widespread iron limitation of ocean productivity, proposed that dust-borne iron fertilization of Southern Ocean phytoplankton caused the ice age reduction in atmospheric carbon dioxide (CO2). In a sediment core from the Subantarctic Atlantic, we measured foraminifera-bound nitrogen isotopes to reconstruct ice age nitrate consumption, burial fluxes of iron, and proxies for productivity. Peak glacial times and millennial cold events are characterized by increases in dust flux, productivity, and the degree of nitrate consumption; this combination is uniquely consistent with Subantarctic iron fertilization. The associated strengthening of the Southern Ocean’s biological pump can explain the lowering of CO2 at the transition from mid-climate states to full ice age conditions as well as the millennial-scale CO2 oscillations.

Stacking of discontinuous regional palaeoclimate records: Chironomid-based summer temperatures from the Alpine region

Since multi-site reconstructions are less affected by site-specific climatic effects and artefacts, regional palaeotemperature reconstructions based on a number of sites can provide more robust estimates of centennial- to millennial-scale temperature trends than individual, site-specific records. Furthermore, reconstructions based on multiple records are necessary for developing continuous climate records over time scales longer than covered by individual sequences. Here, we present a procedure for developing such reconstructions based on relatively short (centuries to millennia), discontinuously sampled records as are typically developed when using biotic proxies in lake sediments for temperature reconstruction. The approach includes an altitudinal correction of temperatures, an interpolation of individual records to equal time intervals, a stacking procedure for sections of the interval of interest that have the same records available, as well as a splicing procedure to link the individual stacked records into a continuous reconstruction. Variations in the final, stacked and spliced reconstruction are driven by variations in the individual records, whereas the absolute temperature values are determined by the stacked segment based on the largest number of records. With numerical simulations based on the NGRIP δ18O record, we demonstrate that the interpolation and stacking procedure provides an approximation of a smoothed palaeoclimate record if based on a sufficient number of discontinuously sampled records. Finally, we provide an example of a stacked and spliced palaeotemperature reconstruction 15000–90 calibrated 14C yr BP based on six chironomid records from the northern and central Swiss Alps and eastern France to discuss the potential and limitations of this approach.

Combining charcoal sediment and molecular markers to infer a Holocene fire history in the Maya Lowlands of Petén, Guatemala

Vegetation changes in the Maya Lowlands during the Holocene are a result of changing climate conditions, solely anthropogenic activities, or interactions of both factors. As a consequence, it is difficult to assess how tropical ecosystems will cope with projected changes in precipitation and land-use intensification over the next decades. We investigated the role offire during the Holocene by combining macroscopic charcoal and the molecular fire proxies levoglucosan, mannosan and galactosan. Combining these two different fire proxies allows a more robust understanding of the complex history of fire re- gimes at different spatial scales during the Holocene. In order to infer changes in past biomass burning, we analysed a lake sediment core from Lake Peten Itza, Guatemala, and compared our results with millennial-scale vegetation and climate change available in the area. We detected three periods of high fire activity during the Holocene: 9500 e 6000 cal yr BP, 3700 cal yr BP and 2700 cal yr BP. We attribute the first maximum mostly to climate conditions and the last maximum to human activities. The rapid change between burned vegetation types at the 3700 cal yr BP fire maximum may result from human activity.