Abstract
Tibito Tephra was first recorded in the central highlands of Papua New Guinea (PNG) in 1971. By the late 1970s, the tephra had been mapped across tens of thousands of square kilometres, traced to its source on Long Island in the Bismarck Sea, and linked to pyroclastic density currents in the Matapun Beds on the island. With a tephra fall volume >10 km3, this eruption was clearly one of the 10 largest globally in the past 600 years. Although almost certainly in the AD 17th or 18th centuries, determining just when this VEI 6 eruption occurred has proved difficult. Whether this eruption occurred before or after William Dampier sailed past, named, described and drew a profile of Long Island has been debated also. Dating the Long Island eruption has implications for assessing the rate of revegetation of the island, recolonisation of the island and caldera lake, Jared Diamond’s theory of ‘supertramp’ birds, correlations between major eruptions and ice core chronologies, the constant rate of supply model of 210Pb concentration in lake sediments, reservoir effects in Lake Kutubu in the southern highlands of PNG, the timing of a prehistoric phase of agriculture in the Western Highlands and the longevity of oral histories recording a taim tudak (time of darkness) when sand fell from the sky, houses collapsed, crops were ruined and people died. This paper reviews early dates based on radiocarbon, 210Pb and paleosecular magnetic variation, historical reports and genealogical dates from oral histories, and speculation based on tree ring and ice core evidence. Analyses of 10 new radiocarbon dates from the Matapun Beds on Long Island are reported. Our best estimate places the eruption between 1651 and 1671 AD with a 95.4% probability and between 1655 and 1665 AD with a 68.2% probability.
Keywords
Introduction
Across the highlands of Papua New Guinea (PNG), local people from a variety of language groups preserved oral histories of a taim tudak (time of darkness) when sand fell from the sky, crops were ruined, a few houses collapsed and people died. Tephrostratigraphic investigations at numerous sites led to the naming of the taim tudak deposit as ‘Tibito Tephra’ (Blong, 1982).
By the late 1970s, Tibito Tephra had been mapped across an area of tens of thousands of square kilometres of the PNG highlands, and had been traced back to its source on Long Island (5.33°S, 147.1°E), an island in the Bismarck Sea to the north of the PNG mainland (Figure 1). At the same time, it became clear that Tibito Tephra was coeval with major pyroclastic flow deposits on the island first briefly described by G. A. M. Taylor (1953), labelled the ‘pyroclastic mantle’ by Ball and Johnson (1976) and formally named the Matapun Beds by Pain et al. (1981). With the tephra fall volume alone exceeding 10 km3, it was clear that the eruption responsible for this tephra could be classified VEI 6.

Study region and important locations. Tibito Tephra isopachs in centimetres, revised with additional data since Blong (1982) in blue. Tephra thicknesses recorded in the field are shown. The heavy dashed line encompasses the area where stratigraphic investigations confirm that Tibito Tephra is the most recent recognisable tephra.
Dating the VEI 6 eruption that deposited Tibito Tephra across the highlands and the associated pyroclastic Matapun Beds on Long Island proved to be difficult, though all estimates lay between early 17th and late 19th century. Radiocarbon dates in both areas demonstrated that the eruption occurred at a time when the relationship between radiocarbon years and calendar years was at its most complex and offering a range of interpretations. The evidence provided by the journals of the few European navigators that had traversed the area since the 16th century proved similarly difficult to interpret, and estimates based on genealogies recorded in oral histories suggested more recent dates.
By the 1980s, evidence that this was a VEI 6 eruption encouraged dendrochronologists and ice core scientists to attempt to match cold years and/or sulphate spikes to the eruption, usually without reference to or an understanding of other information that might constrain estimates of the age. As tree ring counting and ice core analysis improved, eruptions elsewhere were associated with particular tree ring and ice core features, and new radiocarbon calibration curves became available, speculation about the year in which the last major eruption of Long Island occurred became a moveable feast – generally believed to be within the latter half of the 17th century, and without any new evidence based on the age of Tibito Tephra or the Matapun Beds.
This contribution reviews the earlier evidence for the age of the eruption, outlines the reasons that dating the eruption as accurately as possible is important and provides the results of a new program of radiocarbon dating of the Matapun Beds.
Background
Long Island lies 130 km due east of Madang and about 55 km north of the nearest point on the PNG mainland (Figure 1). The maximum N-S length of the island is about 27 km and the maximum width about 25 km, with a land area of about 330 km2 and a large central caldera. Mt Reaumur, the highest point at 1280 m, is located near the northern end of the island. Remarkably little is known about the climate on the island: Ball and Glucksman (1978) estimated the mean annual rainfall at about 2800 mm with the wettest months with common thunderstorms between December and March, and a dry season extending April to November. Thornton (2001) concurs with earlier authors that the vegetation appears to be an arrested sub-climax of fast growing softwood species dispersed by floating seeds or by birds and bats. Above about 300 m elevation, the forest appears more open. Surface streams are inhibited by the deep porous volcanic deposits across much of the island; prolonged dry seasons and droughts are common.
Details of the European discovery and exploration of Long Island (also known as Arup or Ahrup), PNG, are given in Ball (1982b). Although Abel Tasman sailed past the island in late April 1643, he made no comment on it, possibly believing it was a peninsula of the mainland (Sharp, 1968: 228). Long Island was named by the English navigator William Dampier on 31 March 1700, 1 who also drew a profile of the island. Blong et al. (2016) have shown, after determining the likely location from where Dampier’s profile was drawn, that his view of the island is remarkably similar to the present-day profile (Figure 2), suggesting that the last major eruption of Long Island more likely occurred before AD 1700.

Dampier’s profile of Long Island, reproduced from Williamson (1939), compared with a modern view from approximately the same location. The shaded portion in the lower profile shows Dampier’s view after allowing for curvature of the earth (Blong et al., 2016). North is to the right side of the figure.
The central freshwater lake (Figure 3) that occupies most of the island was evidently not known to Europeans until 1928 (Ball, 1982b). The lake has a surface area of a little over 90 km2, with a surface elevation of about 196 m a.s.l. 2 Ball and Glucksman (1978) used an echosounder to produce a bathymetric map of the lake; a similar survey by Rabaul Volcanological Survey staff in 1976 resulted in a reasonably similar map (Chris McKee, personal communication, 2014). Our lake bathymetry, based on a DEM most likely derived from the Ball and Glucksman (1978) map, shows that more than 60% of the lake floor (i.e. more than 55 km2) lies below modern sea level with the deepest point at −173 m to the east of the inner island, Motmot (Figure 4). Eruptions have been recorded seven or eight times from Motmot since 1933 and once (1993) from a linear vent system beneath the east central part of the lake (Global Volcanism Program, 2013b).

A one arc-second (approximately 30 m) resolution SRTM image of Long Island. The elevation of the water surface of Lake Wisdom was 196 m above mean sea level in the SRTM data.

The bathymetry of Lake Wisdom purchased from https://www.tcarta.com in July 2014. The island in the south-central part of Lake Wisdom is Motmot, which has erupted seven or eight times since the first eruption recorded from the lake in 1933.
Based on the bathymetry in Figure 4, the volume of the lake is estimated at about 19.7 km3. The highest point on the caldera rim is approximately 630 m a.s.l. at the northern end; this is almost certainly a ridge of Mt Reaumur. Most of the 40 km long caldera rim lies between elevations of 375 and 500 m a.s.l. If we assume the average height of the rim is 400 m, the caldera volume can be estimated at 40.7 km3; at an average elevation of 500 m, the volume becomes 51.5 km3. Thus, between 30% and 48% of the caldera volume is water-filled.
Immediately after the most recent major eruption, the surface area of the caldera was probably smaller as parts of the caldera walls appear to be eroding. However, as little or no sediment can escape from the lake, the caldera was probably deeper immediately post-eruption. Additionally, Motmot Island and other intra-caldera features have been produced by a number of post-caldera eruptions.
Only the broad details of the volcanic history of Long Island are known. Mt Reaumur (1280 m) at the northern end (the right hand end in Figure 2) and Cerisy Peak (1112 m) at the southern end are ancient basaltic stratovolcanoes, named by the navigator Dumont D’Urville in 1827 (Ball, 1982b; Ball and Johnson, 1976). Prominent satellite cones lie on the southwest flank of Reaumur and east and northwest of Cerisy Peak (Figure 2). Most of the rest of the island is mantled with thick pyroclastic deposits that slope away from the central caldera.
Three major eruptions from the caldera have been identified (Blong et al., 1982; Pain et al., 1981). The earliest of these eruptions is represented by the Kiau Beds, which comprise a thick sequence of airfall ash and lapilli beds overlain by a firm ignimbrite. The outer portion of a completely carbonised log horizontal in the basal beds gave a radiocarbon age of 16,040 ± 270 yr BP (SUA-623). The Kiau Beds have been correlated on geochemical grounds with Ep Tephra widespread in the highlands of mainland PNG (Blong et al., 2017). Two radiocarbon dates from near Kuk (Mt Hagen) provide an age range for Ep tephra between 12,000 and 15,000 radiocarbon years. 3
The Kiau Beds are overlain by the Biliau Beds, which form a similar sequence of airfall deposits followed by pyroclastic density currents (PDCs), ranging 1 to 8 m in thickness. A carbonised log in the basal unit provided a radiocarbon age of 3990 ± 110 yr BP (SUA-624) (Pain et al., 1981). The Biliau Beds are believed to be correlated with Kim Tephra widespread in highland PNG and dated at between 3500 and 3800 radiocarbon years. 4
The uppermost thick eruption sequence on Long Island, the Matapun Beds (Pain et al. 1981), are better exposed and better known. The lowermost unit lies directly on a palaeosol at the top of the Biliau Beds and consists of as much as 2 m of well-bedded coarse airfall ash and pumiceous lapilli. This is overlain by unwelded, poorly bedded PDCs composed of scoriaceous lapilli and non-vesicular blocks in a fine sandy matrix. Accretionary lapilli are common. The PDCs are overlain by up to 1.5 m of discontinuous airfall-bedded ash and lapilli, including accretionary lapilli. On mainland PNG and in the highlands, this eruption deposited Tibito Tephra, mapped across an area of about 94,000 km2 (Figure 1), with a volume of more than 10 km3 and geochemically linked to a Long Island source (Blong, 1982).
The relative ‘sizes’ of these three major eruptions are not known with any precision but Ep Tephra appears to be the thickest of the three tephras on the mainland, with Kim Tephra possibly the thinnest and/or least widespread. All three eruptions have been assigned magnitudes M ≥ 6.3 in the LaMeve (Large Magnitude Explosive Volcanic Eruptions) database (http://www.bgs.ac.uk/vogripa/view/controller.cfc?method=lameve). It is likely that all three eruptions contributed to the caldera volume of 40–50 km3 or more, though it is also probable that minor eruptions built new cones and refilled portions of the caldera at various times between the major events. It is not possible to relate caldera size to the volume of each eruption.
While Long Island has produced three major eruptions of global significance in the last 15,000 years, it is the timing of the eruption that produced the Matapun Beds and Tibito Tephra that is the focus of the current contribution.
Significance of the most recent major eruption
Dating of the most recent major eruption of Long Island has implications for several areas of research:
The airfall tephra volume of more than 10 km3 identifies the latest major eruption of Long Island as a VEI 6 eruption. As only 10 or so VEI 6 or greater eruptions have been recorded globally in the last 600 years, the eruption of Long Island is clearly of global significance. In fact, only 56 out of the 7742 in the Smithsonian catalogue (Siebert et al., 2010) of Holocene eruptions have been of this magnitude or greater (roughly 0.7% of all eruptions in the Holocene). Our recent estimates of airfall tephra volume, based on updated and improved methods of volume estimation, suggest the Long Island eruption was larger than the 1883 Krakatau eruption. 5 A more accurate date will lead to improved estimates of the return periods of such large magnitude events.
Surveys of bird populations on islands in the eastern Bismarck Archipelago show that Long Island, Crown Island and Tolokiwa (Figure 1) have spectacularly high population densities of birds of ‘supertramp’ species, compared with other islands in the region (Diamond, 1974). Supertramps are birds of high overwater dispersal ability but probably low competitive ability. Supertramps are confined to islands placing a premium on dispersal ability – that is, small islands subject to frequent population extinctions, remote islands and islands with recent large volcanic eruptions resulting in defaunation. In fact, Long Island has populations of more supertramp bird species than any other known island in the world (Jared Diamond, 2015, personal communication).
Long Island’s biotic communities appear less developed than expected when compared with islands in the Krakatau group after the 1883 eruption, with the vegetation in an arrested stage of development. Thornton (2001) attributes this at Long Island to the highly porous volcanic soils, poor water retention, low rainfall and frequent extreme seasonal aridity producing slow growth of trees and arrested vegetational development which in turn inhibits maturation of the island’s animal communities.The biota in Lake Wisdom are less developed than that in Lake Dakataua in West New Britain, which formed following a caldera-forming eruption in the AD 7th century (McKee et al., 2011). Ball and Glucksman (1980) hypothesise this results from the shorter time since the creation of the lake and the greater distance of Long Island from sources of colonists (Thornton, 2001).
Tibito Tephra is readily identifiable in numerous mainland PNG lakes and alpine bogs where its use as a tephrostratigraphic marker has contributed to discussions about rates of sedimentation and infilling, nutrient accession and human modification of catchments (Hope et al., 1988). The presence and estimated age of Tibito Tephra in Lakes Egari and Ipea in the highlands has also been crucial to the prominence of the constant rate of supply (c.r.s.) model versus the constant initial concentration (c.i.c.) model of 210Pb concentrations in lake sediments (Oldfield et al., 1978). Tibito Tephra has also been used to correlate between lake cores measuring recent palaeomagnetic secular variation (Thompson and Oldfield, 1978).
Recent radiocarbon dating of lake sediments in Lake Kutubu in the southern highlands of mainland PNG (see Figure 1) indicates substantial reservoir effects. Geochemical identification of Tibito and Olgaboli Tephras in the stratigraphic sequence suggests that the reservoir effect amounts to more than 2000 years at the level of Tibito Tephra (Schneider et al., in preparation).
At the Kuk UNESCO World Heritage Site (Golson, 2017) in the Western Highlands (where Tibito Tephra was first identified in 1971), Tibito and other tephra were vital tephrostratigraphic markers elucidating the various phases of agricultural activity at the site. Tibito provides an important marker at the end of Phase 5 of the agricultural history (Golson et al., 2017). The abandonment of wetland cultivation at Kuk took place just before the fall of Tibito Tephra, whereas the most dramatic period of forest disturbance at numerous highland sites occurs immediately after Tibito Tephra fell with evidence of increased erosion (Corlett, 1984). Dating Tibito Tephra also provides crucial evidence for the development of agroforestry based on Casuarina oligodon, increased forest disturbance, increased soil erosion and the appearance of raised-bed garden forms and ditching features in the archaeological record (Haberle, 1998a).
During the 1970s, it became apparent that widespread stories of a taim tudak (time of darkness) across the highlands of PNG were indigenous recognition of the fall of tephra; in most locations across at least 25 language groups, Tibito Tephra was the youngest identifiable tephra (Figure 1). Blong (1982) showed that many of the legends were essentially accurate accounts of the probable consequences of a 10- to 100-mm-thick fall of volcanic ash, except in those aspects relating to the likely duration and timing of the ash fall. Although many accounts included a genealogical basis for the timing of the ash fall (Blong, 1982; Lacey, 1990; Mai, 1981), and the ash fall and legend had earlier been attributed to the 1883 eruption of Krakatau (Watson, 1962), more thorough examination revealed that the oral histories of the tephra fall had survived more or less intact for far longer than the few generations proposed by many of the story tellers. Accurate dating of the fall of Tibito Tephra improves our understanding of the longevity of these oral histories, perhaps with implications about the timing of events portrayed in oral histories elsewhere.
Across much of the highlands, the only historical time that can be accurately determined is the first arrival of European man (usually in or after the 1930s). However, some legends refer to events before or after the taim tudak; for example, several accounts state that the people concerned lived elsewhere when the ash fall occurred or that they migrated to or from another place before or after the ash fall. Several of the oral histories collected about the taim tudak refer to sweet potato vines or tubers (e.g. Nelson, 1971), implying the presence of the food crop at the time of the eruption. 6
Along the Madang coast, oral histories refer to a major eruption and disappearance of a high island west of Long Island. These legends consistently place the Yomba eruption (Mennis, 2007; Global Volcanism Program, 2013a) at least two generations before the Long Island eruption (Mennis, 2007). While some authorities believe that Yomba was located at or near Hankow Reef (location on Figure 1) northwest of Long Island, an eruption in Holocene times from this source seems very unlikely (Eli Silver, 2014, personal communication; Silver et al., 2009). Although Blong (1982) implied that Yomba may have been the source volcano for Olgaboli Tephra, widespread across the highlands of PNG, this tephra, dated at 980–1180 yr cal BP, probably erupted from Karkar Island (Schneider et al., 2017).
In each of the above examples, improved dating of the most recent major eruption of Long Island provides crucial constraints on the timing, on the rate of operation and/or rate of change of environmental processes, and on our understanding of the human environment in PNG.
Age estimates for the most recent major eruption of Long Island
The timing of the most recent major eruption of Long Island has been of interest for more than a hundred years. Figure 5 summarises views about the age of the eruption based on previous publications (see also Appendix 1). Most of the estimates lie in the range early–mid-17th century to 1883, the latter because the ash fall that produced the taim tudak story in part of the Eastern Highlands was assumed to result from the cataclysmic eruption of Krakatau in 1883 (Watson, 1962).

The number of published age estimates for the most recent major eruption of Long Island and/or Tibito Tephra in the period AD 1630–1910; G = estimates based on geological evidence; O = oral histories, both from Long Island and from the highlands; RC = radiocarbon dates; V = various lines of evidence; and VF = estimates appearing in the ice core and tree ring literature (see Appendix 1). The dashed line represents the time of Dampier’s voyage.
Age estimates made before the mid-1970s are based entirely on oral histories (O in Figure 5) and nearly all suggest AD post-1700. Bamler (1912), Nurton (1932a) and Carey’s (1938) accounts are based on stories about the eruption of Long Island itself. Vicedom and Tischner (1943, Volume 3, p. 50), early Catholic missionaries among the recently ‘discovered’ highland populations around Mt Hagen, recorded the taim tudak story in the 1930s, the first record of the very widespread accounts of the fall of Tibito Tephra from highland PNG. 7 These recorders heard a local story about a ‘time of darkness’ when sand fell from the sky, crops, animals and many people died and houses collapsed. These oral accounts frequently included genealogies which referred to an ancestor who was alive at the time of the event; an age estimate can then be made by assigning 25–30 years to each generation before that of the story teller (Blong, 1979b; 1982). The dashed line on Figure 1 encloses the area in which we are certain Tibito Tephra is the youngest identifiable volcanic ash.
Geological (G in Figure 5) or volcanological evidence of a ‘recent’ major eruption on Long Island was first recognised by G. A. M. Taylor in 1953, although the two samples he later submitted for radiocarbon dating (NZ R332 and R517) produced ‘Modern’ 8 results and it has not been possible to identify accurately the locations near the Lake Wisdom crater wall from where these samples were obtained. Three samples collected by Ian Hughes (ANU-1125-1127) from the coast southwest of Poin Kiau (Figure 3) in the early 1970s were also not stratigraphically secure, but the stratigraphy along this coast exposes only Matapun Beds (Pain et al., 1981).
Blong (1982) used four radiocarbon age estimates from Long Island and six from the Kuk prehistoric site in the Western Highlands of PNG. Kuk samples were corrected to allow for (constant) sedimentation rates in the Kuk swamp, the means of the various samples were pooled and the nascent correction curves were examined, to decide that the Stuiver (1978) curves were likely to be the most reliable. 9 It was concluded: These calibrated radiocarbon dates place the time of darkness, the fall of Tibito Tephra and the eruption of Long Island fairly firmly in the early–mid 17th century. However, consideration of 2σ each side of the mean ages indicates a possibility that the eruption occurred in the early 18th century or, for that matter, in the late 16th century (Blong 1982: 193). Polach (1981) reached essentially the same conclusion using a slightly different route based on the same initial samples.
Blong also examined genealogical estimates of the timing of the occurrence of the time of darkness based on Mai’s 24 interviews (Mai, 1981) in Enga Province and a further 29 accounts collected from across the highlands (Blong, 1981), concluding that most oral histories suggested that the eruption occurred in the mid-18th century or more recently, even when the ‘gap’ between generations was assumed to be an unlikely long 30 years.
Oldfield and co-workers reported an estimated age for Tibito Tephra of ~AD 1860 based on 210Pb concentration versus depth curves in Lakes Egari and Pipiak (see locations on Figure 1), based on a constant initial concentration model. Using the more realistic constant rate of supply model provided an age range of AD 1800–1840 (with a mean of AD 1814) (Oldfield et al., 1978). Subsequently, improved measurement techniques for 210Pb and 226Ra using direct gamma assay altered this estimate to AD 1680–1690 (Oldfield et al., 1980), by curve fitting through points fixed at 130–160 years by 210Pb and 14C dates >1000 years (210Pb has a half-life of 22.6 years and this updated estimate is beyond the range of the 210Pb technique). Thompson and Oldfield (1978) also compared palaeomagnetic secular variation in the Lake Ipea sediments with AD post-1650 secular variation based on all observatory records, concluding that Tibito Tephra was deposited AD pre-1700.
Other evidence of the possible age of the eruption included Dampier’s profile of Long Island drawn in AD 1700 (Figure 2) and his description that: ‘The Trees appeared very green and flourishing …’ (Williamson, 1939). The two lines of evidence suggest that the eruption must have occurred at least 20–30 years prior to his visit (i.e. before AD 1670–1680 to explain the unchanged shape of the island and to allow time for revegetation). Limited historical evidence provided by subsequent European navigators suggested that the eruption could not have taken place after about AD 1800 (Ball, 1982b; Blong, 1982: 193).
After weighing all the available evidence and favouring the radiocarbon dates, Blong (1982: 194) concluded that in the mid-17th century, AD 1630–1670 was the most likely time for the last major eruption of Long Island and the associated events to have occurred. This age estimate was recorded in the Smithsonian Institution’s Volcanoes of the World (1981) as 1660 ± 20 years based on uncorrected radiocarbon years and remains recorded that way today (Global Volcanism Program, 2013b).
Interpreting the Long Island radiocarbon results, and those dating the fall of Tibito Tephra across the highlands of PNG, has not been straightforward as illustrated in Figures 6 and 7. For example, ANU-1126 is a sample of the outer rings of a tree trunk carbonised in a PDC of the Matapun Beds near Poin Kiau (one of the samples collected by Ian Hughes). When the sample was analysed in 1973, the result, 230 ± 75 yr BP (Before Present = 1950), implied a calendar year age of AD 1720, with a one standard deviation range of 1645–1795 and a two standard deviation range of 1570–1870. Figure 6 illustrates that there is not a simple one-to-one correspondence between radiocarbon years and calendar years and that one of the periods most difficult to interpret is AD 1660–1950.

Radiocarbon calibration curve for ANU-1126 (collected from the northwest coast of Long Island in the early 1970s) illustrating the difficulties in converting radiocarbon years to calendar years.

All dates available for the Long Island eruption, calibrated using OxCal 4.2 and the IntCalnorthern hemisphere calibration curve.
During the last few decades, each new version of the calibration curves gave slightly different results, with none of them resolving satisfactorily one of the central issues – whether the last major eruption occurred before or after William Dampier’s description of Long Island in AD 1700. As Figure 2 shows, there is excellent agreement between Dampier’s 1700 profile of Long Island and a modern profile, but this similarity was less convincing until recently when it became possible to pinpoint the location from where Dampier’s profile had been drawn (Blong et al., 2016). As almost all the slopes on Long Island facing east and between the peaks of Reaumur and Cerisy are mantled with thick deposits of Matapun Beds (Pain et al., 1981), it is now reasonably certain that Dampier’s profile postdates the last major eruption; that is, the eruption occurred at a time before AD 1700.
Dampier also provided a description of the vegetation on Long Island in 1700:
The 31st in the Forenoon we shot in between 2 Islands, lying about 4 Leagues asunder; with Intention to pass between them. The Southernmost is a long Island, with a high Hill at each End; this I named Long Island. The Northernmost is a round high Island towering up with several Heads or Tops, something resembling a Crown; This I named Crown-Isle, from its Form. Both these Islands appear’d very pleasant, having Spots of green Savannahs mixt among the Wood-land: The Trees appeared very green and flourishing, and some of them looked white and full of Blossoms. (Dampier, 1729, in Williamson, 1939)
Although there has been subsequent discussion about the state of the vegetation (Ball, 1982a; Ball and Johnson, 1976), it is clear by analogy with the 1883 Krakatau eruption and vegetation recovery that at least two or three decades must have passed between the eruption and Dampier’s voyage (see, for example, Thornton, 2000, 2001).
By the early 1980s, it was evident that sulphate spikes in the Greenland and Antarctic ice caps recorded a sequence of major volcanic eruptions (e.g. Hammer et al., 1980). As ice cores record more or less annual snow deposition and it is reasonably straightforward to identify the signal of large tropical eruptions such as Krakatau (in 1883) and Tambora (in 1815), the search began for large eruptions that might match other sulphate spikes. As large eruptions also produce global cooling, it was possible to match (particularly European) tree ring patterns to eruptions. The early efforts of Zielinski et al. (1994) and Briffa et al. (1998) were certainly cognisant of the arguments limiting the date of Long Island’s most recent major eruption presented above. However, some subsequent ‘volcanic forcing’ publications (VF in Figure 5) failed to recognise that at least 20 or 30 years were necessary after the eruption for revegetation of the island before Dampier’s AD 1700 description; the authors were seemingly unaware of any of the discussions about the age or arguments for and against timing of the eruption. In the light of Dampier’s description of the vegetation, an AD 1675 date for the eruption seems marginally possible and an AD 1695 date quite unlikely.
These issues were compounded by the problem that ice core records were not always annual as seasonal accumulation layers could not always be easily identified, and erosion sometimes removed other layers so that matches between cores were less than perfect, requiring revision. Similarly, the match between tree rings and ice cores was not always error-free (Baillie, 2010). The efforts of Gao et al. (2008) in comparing numerous ice cores from both the Arctic and Antarctic appears to have reduced errors inherent in reconstructions based on small numbers of cores. However, as Figure 6 and Appendix 1 show, the Long Island eruption has been matched to a number of sulphate spikes and to calendar years. While the VEI 6 eruption of Long Island was certainly large enough, sulphurous enough (Schonwalder et al., 2017) and close enough to the equator to produce a recognisable sulphate spike in ice accumulations in both Northern and Southern Hemispheres, it is not yet possible to match the eruption with any particular spike. 10
New 14C dates
In June 2014, fresh samples of charcoalised logs from the Matapun Beds on Long Island were collected. Several logs were found in multiple locations around the island. All logs were small and did not contain visible rings preventing any wiggle matching attempts. Charcoal analysis determined the logs were from multiple taxa, Cocos nucifera, Pandanus cf. tectorius, Terminalia cf. catappa and Calophyllum inophyllum (M. Prebble, personal communication, 2017). Samples were returned to the Australian National University Radiocarbon Laboratory for subsequent analysis. A total of 10 samples were selected for radiocarbon dating. Subsamples weighing ~5mg were dried, then pretreated using standard acid:base:acid techniques, 1M HCl for 1 hr, multiple 1M NaOH rinses until clear, 1M HCl, then repeated rinses with high purity water (milli-Q). Approximately 1.5 mg of material was converted to CO2 at 900°C in a sealed quartz tube containing ~60mg CuO and silver. All charcoal samples were >50% carbon indicating well-preserved samples. The resulting CO2 was purified and converted to graphite in the presence of hydrogen with powdered Fe as a catalyst. The graphite was then measured on the ANU single stage accelerator mass spectrometer (Fallon et al., 2010). All data were corrected using online AMS d13C, normalised to Oxalic Acid I and background subtracted using 14C free coal treated in the same ABA manner. Data are presented according to the recommendations of Stuiver and Polach (1977).
OxCal 4.2 (Bronk Ramsey, 2009a, 2009b) was used to calibrate the radiocarbon ages to calendar ages. Consistent with previous studies (McKee et al., 2015), we used the northern hemisphere IntCal dataset for calibration as PNG lies close to the northern hemisphere and the SHCal dataset is primarily made up from wood south of 35°S. Figure 6 shows that the eruption date falls in a very difficult period for radiocarbon calibration, with multiple calendar ages possible for a conventional radiocarbon age. Performing a standard calibration on the available Long Island material results in a calendar age eruption estimate of ~1400–present (Table 1, Figure 7).
All radiocarbon samples available for Long Island. Includes OxCal 4.2 calibrated calendar ages and outlier modelled ages.
In order to refine the estimate of the eruption age, Bayesian modelling was performed in OxCal (Bronk Ramsey, 2009b). We used a standard phase model with outlier analysis to determine a date for the eruption at Long Island (Bronk Ramsey, 2009a). The general outlier model was used with a 0.05 probability of any sample being an outlier (Bronk Ramsey, 2009b). We only used the recently collected samples (those with SANU sample numbers) in the phase model as the providence is robust and we are confident that the outer edge of the logs corresponds to a single event. The modelled calibrated ages for each sample are shown in Figure 8. The modelled date for the eruption has a 1 sigma calendar date range of 1649–1665 AD (68.2%) and a 2 sigma range of 1640–1675 AD (92.5%) and 1785–1793 AD (2.9%).

(a) Phase model results from OxCal 4.2 for the recently collected carbon logs from Long Island. (b) Calibrated calendar age estimate of the Long Island eruption from the radiocarbon measurements of the recently collected carbon logs.
If we add a priori knowledge to the OxCal phase model using the evidence from Dampier’s observations in 1700 and allowing a minimum 20 years for post-eruption vegetation growth (i.e. the eruption had to occur before 1680 AD), we can further constrain the eruption date to 1651–1671 (95.4% probability, Figure 9). 11

Long Island eruption date estimated from the OxCal phase model using the AD 1643 date of Abel Tasman’s visit and AD 1680 (to allow 20 years for revegetation of the island) before Dampier’s observations in AD 1700 as constraints for the eruption period.
Conclusion
For nearly the last five decades, the date of the most recent major eruption of Long Island has been of interest, contributing to debates of historical, archaeological, biological, volcanological and climatic interest. In the 1970s, estimates of the timing of the eruption ranged from the mid-17th century to the very early 20th century with all of the post-1800 AD estimates relying on genealogies reported in oral histories (Figure 5). Estimates based on 210Pb concentrations and palaeomagnetic secular variations in lake sediments, supported by radiocarbon dates, then moved estimates into the late 17th or early 18th century. William Dampier’s AD 1700 profile of Long Island (Figure 2) and description of the vegetation added to the confusion about the timing of the eruption.
Resolution of the age of the eruption was complicated by the remarkably non-linear relationship between radiocarbon ages and calibrated dates (Figure 6) and the minor changes engendered by each new calibration. Nonetheless, in the 1980s, the Smithsonian Institution recorded the age of the eruption as 1660 ± 20 years (described as ‘uncorrected’ though that was not the case). Undoubtedly, this age range for a VEI 6 eruption stirred the interest of those researching tree rings and ice cores and the suggested matching to a range of dates between 1640 and 1695, usually without knowledge of the constraints imposed by Dampier’s observations.
The 10 new radiocarbon dates reported here further constrain the possible timing, with our best estimate producing a 95.4% probability that the eruption occurred between AD 1651 and 1671.
Further refinement of this estimate will depend on wiggle matching of charcoalised logs exhibiting annual rings in the Matapun Beds or a convincing match to ice core and/or tree ring records.
Footnotes
Appendix 1
Published estimates of the timing of the last major eruption of Long Island and/or the fall of Tibito Tephra. Estimates are presented in order of publication or, where known, the year in which the information was collected.a Calendar year estimates are often less precise than they appear – often central estimates ±20 years. Not infrequently, there are gaps of 5 or more years between record and publication. Calendar year estimates are then grouped by decades in Figure 5.
Acknowledgements
We gratefully acknowledge the following individuals: Jack Golson who initiated the search for tephra in archaeological sites in the Western Highlands in 1970 and his continuing interest in our efforts to refine the dating; Ian Hughes who first proposed, in the mid-1970s, that Long Island might be the source of Tibito Tephra; Colin Pain who accompanied CMcK and RB to Long Island in 1976 and extended far and wide our understanding of the distribution of Tibito Tephra; Geoff Hope who over more than 40 years cored far-flung lakes and bogs that contained Tibito across the highlands; Mike Prebble for his identification of the tree species in our 14C sample collection in 2014; and Paul Grasso for his manipulation of spreadsheets. We are extremely grateful also to the staff of the Madang Resort Hotel, Andrew Barter and the crew of the Kalibobo Spirit for making the visit by RB, SF and CMcK to Madang and Long Island in 2014 so pleasant and rewarding. We also want to thank Frank Oldfield and two anonymous reviewers for their valuable comments on the manuscript.
Funding
The author(s) received no financial support for the research, authorship and/or publication of this article.
