Ages of 24 widespread tephras erupted since 30,000 years ago in New Zealand, with re-evaluation of the timing and palaeoclimatic implications of the Lateglacial cool episode recorded at Kaipo bog

Tephras are important for the NZ-INTIMATE project because they link all three records comprising the composite inter-regional stratotype developed for the New Zealand climate event stratigraphy (NZ-CES). Here we firstly report new calendar ages for 24 widespread marker tephras erupted since 30,000 calendar (cal.) years ago in New Zealand to help facilitate their use as chronostratigraphic dating tools for the NZ-CES and for other palaeoenvironmental and geological applications. The selected tephras comprise 12 rhyolitic tephras from Taupo, nine rhyolitic tephras from Okataina, one peralkaline rhyolitic tephra from Tuhua, and one andesitic tephra each from Tongariro and Egmont/Taranaki volcanic centres. Age models for the tephras were obtained using three methods: (i) C-based wiggle-match dating of wood from trees killed by volcanic eruptions (these dates published previously); (ii) flexible depositional modelling of a high-resolution C-dated age-depth sequence at Kaipo bog using two Bayesian-based modelling programs, Bacon and OxCal's P_Sequence function, and the IntCal09 data set (with SH offset correction-44±17yr); and (iii) calibration of C ages using OxCal's Tau_Boundary function and the SHCal04 and IntCal09 data sets. Our preferred dates or calibrated ages for the 24 tephras are as follows (youngest to oldest, all mid-point or mean ages of 95% probability ranges): Kaharoa AD 1314±12; Taupo (Unit Y) AD 232±10; Mapara (Unit X) 2059±118cal.yrBP; Whakaipo (Unit V) 2800±60cal.yrBP; Waimihia (Unit S) 3401±108cal.yrBP; Stent (Unit Q) 4322±112cal.yrBP; Unit K 5111±210cal.yrBP; Whakatane 5526±145cal.yrBP; Tuhua 6577±547cal.yrBP; Mamaku 7940±257cal.yrBP; Rotoma 9423±120cal.yrBP; Opepe (Unit E) 9991±160cal.yrBP; Poronui (Unit C) 11,170±115cal.yrBP; Karapiti (Unit B) 11,460±172cal.yrBP; Okupata 11,767±192cal.yrBP; Konini (bed b) 11,880±183cal.yrBP; Waiohau 14,009±155cal.yrBP; Rotorua 15,635±412cal.yrBP; Rerewhakaaitu 17,496±462cal.yrBP; Okareka 21,858±290cal.yrBP; Te Rere 25,171±964cal.yrBP; Kawakawa/Oruanui 25,358±162cal.yrBP; Poihipi 28,446±670cal.yrBP; and Okaia 28,621±1428cal.yrBP.Secondly, we have re-dated the start and end of the Lateglacial cool episode (climate event NZce-3 in theNZ-CES), previously referred to as the Lateglacial climate reversal, as defined at Kaipo bog in eastern North Island, New Zealand, using both Bacon and OxCal P_Sequence modelling with the IntCal09 data set. The ca1200-yr-long cool episode, indicated by a lithostratigraphic change in the Kaipo peat sequence to grey mudwith lowered carbon content, and a high-resolution pollen-derived cooling signal, began 13,739±125cal.yrBP and ended 12,550±140cal.yrBP (mid-point ages of the 95% highest posterior density regions, Bacon modelling). The OxCal modelling, generating almost identical ages, confirmed these ages. The Lateglacial cool episode (ca 13.8-12.6cal.kaBP) thus overlaps a large part of the entire Antarctic Cold Reversal chronozone (ca 14.1-12.4cal.kaBP or ca 14.6-12.8cal.kaBP), and an early part of the Greenland Stadial-1 (Younger Dryas) chronozone (ca 12.9-11.7cal.kaBP). The timing of the Lateglacial cool episode at Kaipo is broadly consistent with the latitudinal patterns in the Antarctic Cold Reversal signal suggested for the New Zealand archipelago from marine and terrestrial records, and with records from southern South America. © 2012 Elsevier Ltd.

A composite pollen-based stratotype for inter-regional evaluation of climatic events in New Zealand over the past 30,000 years (NZ-INTIMATE project)

Our review of paleoclimate information for New Zealand pertaining to the past 30,000 years has identified a general sequence of climatic events, spanning the onset of cold conditions marking the final phase of the Last Glaciation, through to the emergence to full interglacial conditions in the early Holocene. In order to facilitate more detailed assessments of climate variability and any leads or lags in the timing of climate changes across the region, a composite stratotype is proposed for New Zealand. The stratotype is based on terrestrial stratigraphic records and is intended to provide a standard reference for the intercomparison and evaluation of climate proxy records. We nominate a specific stratigraphic type record for each climatic event, using either natural exposure or drill core stratigraphic sections. Type records were selected on thebasis of having very good numerical age control and a clear proxy record. In all cases the main proxy of the type record is subfossil pollen. The type record for the period from ca 30 to ca 18 calendar kiloyears BP (cal. ka BP) is designated in lake-bed sediments from a small morainic kettle lake (Galway tarn) in western South Island. The Galway tarn type record spans a period of full glacial conditions (Last Glacial Coldest Period, LGCP) within the Otira Glaciation, and includes three cold stadials separated by two cool interstadials. The type record for the emergence from glacial conditions following the termination of the Last Glaciation (post-Termination amelioration) is in a core of lake sediments from a maar (Pukaki volcanic crater) in Auckland, northern North Island, and spans from ca 18 to 15.64±0.41 cal. ka BP. The type record for the Lateglacial period is an exposure of interbedded peat and mud at montane Kaipo bog, eastern North Island. In this high-resolution type record, an initial mild period was succeeded at 13.74±0.13 cal. ka BP by a cooler period, which after 12.55±0.14 cal. ka BP gave way to a progressive ascent to full interglacial conditions that were achieved by 11.88±0.18 cal. ka BP. Although a type section is not formally designated for the Holocene Interglacial (11.88±0.18 cal. ka BP to the present day), the sedimentary record of Lake Maratoto on the Waikato lowlands, northwestern North Island, is identified as a prospective type section pending the integration and updating of existing stratigraphic and proxy datasets, and age models. The type records are interconnected by one or more dated tephra layers, the ages of which are derived from Bayesian depositional modelling and OxCal-based calibrations using the IntCal09 dataset. Along with the type sections and the Lake Maratoto record, important, well-dated terrestrial reference records are provided for each climate event. Climate proxies from these reference records include pollen flora, stable isotopes from speleothems, beetle and chironomid fauna, and glacier moraines. The regional composite stratotype provides a benchmark against which to compare other records and proxies. Based on the composite stratotype, we provide an updated climate event stratigraphic classification for the New Zealand region. © 2013 Elsevier Ltd.

The lateritic profile of Balkouin, Burkina Faso: geochemistry, mineralogy and genesis

This study reports on the geochemical and mineralogical characterization of a lateritic profile cropping out in the Balkouin area, Central Burkina Faso, aimed at obtaining a better understanding of the processes responsible for the formation of the laterite itself and the constraints to its development. The lateritic profile rests on a Paleoproterozoic basement mostly composed of granodioritic rocks related to the Eburnean magmatic cycle passing upwards to saprolite and consists of four main composite horizons (bottom to top): kaolinite and clay-rich horizons, mottled laterite and iron-rich duricrust. In order to achieve such a goal, a multi-disciplinary analytical approach was adopted, which includes inductively coupled plasma (ICP) atomic emission and mass spectrometries (ICP-AES and ICP-MS respectively), X-ray powder diffraction (XRPD), scanning electron microscopy with energy dispersive spectrometry (SEM-EDS) and micro-Raman spectroscopy.

The geochemical data, and particularly the immobile elements distribution and REE patterns, show that the Balkouin laterite is the product of an in situ lateritization process that involved a strong depletion of the more soluble elements (K, Mg, Ca, Na, Rb, Sr and Ba) and an enrichment in Fe; Si was also removed, particularly in the uppermost horizons. All along the profile the change in composition is coupled with important changes in mineralogy. In particular, the saprolite is characterized by occurrence of abundant albitic plagioclase, quartz and nontronite; kaolinite is apparently absent. The transition to the overlying lateritic profile marks the breakdown of plagioclase and nontronite, thus allowing kaolinite to become one of the major components upwards, together with goethite and quartz. The upper part of the profile is strongly enriched in hematite (+ kaolinite). Ti oxides (at least in part as anatase) and apatite are typical accessory phases, while free aluminum hydroxides are notably absent. Mass change calculations emphasize the extent of the mass loss, which exceeds 50 wt% (and often 70 wt%) for almost all horizons; only Fe was significantly concentrated in the residual system.

The geochemical and mineralogical features suggest that the lateritic profile is the product of a continuous process that gradually developed from the bedrock upwards, in agreement with the Schellmann classic genetic model. The laterite formation must have occurred at low pH (? 4.5) and high Eh (? 0.4) values, i.e., under acidic and oxidizing environments, which allowed strongly selective leaching conditions. The lack of gibbsite and bohemite is in agreement with the compositional data: the occurrence of quartz (± amorphous silica) all along the profile was an inhibiting factor for the formation of free aluminum hydroxides.

Using geostatistical Bayesian Updating to integrate airborne radiometrics and soil geochemistry to improve mapping for mineral exploration

Mineral exploration programmes around the world use data from remote sensing, geophysics and direct sampling. On a regional scale, the combination of airborne geophysics and ground-based geochemical sampling can aid geological mapping and economic minerals exploration. The fact that airborne geophysical and traditional soil-sampling data are generated at different spatial resolutions means that they are not immediately comparable due to their different sampling density. Several geostatistical techniques, including indicator cokriging and collocated cokriging, can be used to integrate different types of data into a geostatistical model. With increasing numbers of variables the inference of the cross-covariance model required for cokriging can be demanding in terms of effort and computational time. In this paper a Gaussian-based Bayesian updating approach is applied to integrate airborne radiometric data and ground-sampled geochemical soil data to maximise information generated from the soil survey, to enable more accurate geological interpretation for the exploration and development of natural resources. The Bayesian updating technique decomposes the collocated estimate into a production of two models: prior and likelihood models. The prior model is built from primary information and the likelihood model is built from secondary information. The prior model is then updated with the likelihood model to build the final model. The approach allows multiple secondary variables to be simultaneously integrated into the mapping of the primary variable. The Bayesian updating approach is demonstrated using a case study from Northern Ireland where the history of mineral prospecting for precious and base metals dates from the 18th century. Vein-hosted, strata-bound and volcanogenic occurrences of mineralisation are found. The geostatistical technique was used to improve the resolution of soil geochemistry, collected one sample per 2 km2, by integrating more closely measured airborne geophysical data from the GSNI Tellus Survey, measured over a footprint of 65 x 200 m. The directly measured geochemistry data were considered as primary data in the Bayesian approach and the airborne radiometric data were used as secondary data. The approach produced more detailed updated maps and in particular maximized information on mapped estimates of zinc, copper and lead. Greater delineation of an elongated northwest/southeast trending zone in the updated maps strengthened the potential to investigate stratabound base metal deposits.

North Atlantic marine radiocarbon reservoir ages through Heinrich event H4: A new method for marine age model construction

Cooling and sinking of dense saline water in the Norwegian–Greenland Sea is essential for the formation of North Atlantic Deep Water. The convection in the Norwegian–Greenland Sea allows for a northward flow of warm surface water and southward transport of cold saline water. This circulation system is highly sensitive to climate change and has been shown to operate in different modes. In ice cores the last glacial period is characterized by millennial-scale Dansgaard–Oeschger (D–O) events of warm interstadials and cold stadials. Similar millennial-scale variability (linked to D–O events) is evident from oceanic cores, suggesting a strong coupling of the atmospheric and oceanic circulations system. Particularly long-lasting cold stadials correlate with North Atlantic Heinrich events, where icebergs released from the continents caused a spread of meltwater over the northern North Atlantic and Nordic seas. The meltwater layer is believed to have caused a stop or near-stop in the deep convection, leading to cold climate. The spreading of meltwater and changes in oceanic circulation have a large influence on the carbon exchange between atmosphere and the deep ocean and lead to profound changes in the 14C activity of the surface ocean. Here we demonstrate marine 14C reservoir ages (R) of up to c. 2000 years for Heinrich event H4. Our R estimates are based on a new method for age model construction using identified tephra layers and tie-points based on abrupt interstadial warmings.

Rapid magmatic processes accompany arc-continent collision: the Western Bismarck arc, Papua New Guinea

New U-Th-Ra, major and trace element, and Sr-Nd-Pb isotope data are presented for young lavas from the New Britain and Western Bismarck arcs in Papua New Guinea. New Britain is an oceanic arc, whereas the latter is the site of an arc-continent collision. Building on a recent study of the Manus Basin, contrasts between the two arcs are used to evaluate the processes and timescales of magma generation accompanying arc-continent collision and possible slab detachment. All three suites share many attributes characteristic of arc lavas that can be ascribed to the addition of a regionally uniform subduction component derived from the subducting altered oceanic crust and sediment followed by dynamic melting of the modified mantle. However, the Western Bismarck arc lavas diverge from the Pb isotope mixing array formed by the New Britain and the Manus Basin lavas toward elevated Pb-208/Pb-204. We interpret this to reflect a second and subsequent addition of sediment melt at crustal depth during collision. U-238 and Ra-226 excesses are preserved in all of the lavas and are greatest in the Western Bismarck arc. High-Mg andesites with high Sr/Y ratios in the westernmost arc are attributed to recent shallow mantle flux melting at the slab edge. Data for two historical rhyolites are also presented. Although these rhyolites formed in quite different tectonic settings and display different geochemical and isotopic compositions, both formed from mafic parents within millennia.

Air pollutant contamination and acidification of surface waters in the North York Moors, UK: Multi-proxy evidence from the sediments of a moorland pool

Despite the extensive geographical range of palaeolimnological studies designed to assess the extent of surface water acidification in the United Kingdom during the 1980s, little attention was paid to the status of surface waters in the North York Moors (NYM). In this paper, we present sediment core data from a moorland pool in the NYM that provide a record of air pollution contamination and surface water acidification. The 41-cm-long core was divided into three lithostratigraphic units. The lower two comprise peaty soils and peats, respectively, that date to between approximately 8080 and 6740 cal. BP. The uppermost unit comprises peaty lake muds dating from between approximately ad 1790 and the present day (ad 2006). The lower two units contain pollen dominated by forest taxa, whereas the uppermost unit contains pollen indicative of open landscape conditions similar to those of the present. Heavy metal, spheroidal carbonaceous particle, mineral magnetics and stable isotope analysis of the upper sediments show clear evidence of contamination by air pollutants derived from fossil-fuel combustion over the last c. 150years, and diatom analysis indicates that the naturally acidic pool became more acidic during the 20th century. We conclude that the exceptionally acidic surface waters of the pool at present (pH=c. 4.1) are the result of a long history of air pollution and not because of naturally acidic local conditions. We argue that the highly acidic surface waters elsewhere in the NYM are similarly acidified and that the lack of evidence of significant recovery from acidification, despite major reductions in the emissions of acidic gases that have taken place over the last c. 30years, indicates the continuing influence of pollutant sulphur stored in catchment peats, a legacy of over 150years of acid deposition.

Young, Old, and Weathered Carbon—Part 1: Using Radiocarbon and Stable Isotopes to Identify Carbon Sources in an Alkaline, Humic Lake

This article presents a case study of Lower Lough Erne, a humic, alkaline lake in northwest Ireland, and uses the radiocarbon method to determine the source and age of carbon to establish whether terrestrial carbon is utilized by heterotrophic organisms or buried in sediment. Stepped combustion was used to estimate the degree of the burial of terrestrial carbon in surface sediment. ∆14C, δ13C, and δ15N values were measured for phytoplankton, dissolved inorganic carbon (DIC), dissolved organic carbon (DOC), and particulate organic carbon (POC). ∆14C values were used to indicate the presence of different sources of carbon, including bedrock-derived inorganic carbon, “modern,” “recent,” “subsurface,” and “subfossil” terrestrial carbon in the lake. The use of 14C in conjunction with novel methods (e.g. stepped combustion) allows the determination of the pathway of terrestrial carbon in the system, which has implications for regional and global carbon cycling.