ELSEVIER Palaeogeography, Palaeoclimatology, Palaeoecology 152 (1999) 101±127
Environmental and sea-level changes on Banks Peninsula (Canterbury, New Zealand) through three glaciation±interglaciation cycles
James Shulmeister a,Ł, Jane M. Soons b,GlennW.Bergerc, Margaret Harper a, Sarah Holt a, Neville Moar d, John A. Carter a a School of Earth Sciences, Victoria University, P.O. Box 600, Wellington, New Zealand b Department of Geography, University of Canterbury, Private Bag 4800, Christchurch, New Zealand c Quaternary Sciences Centre, Desert Research Institute, Reno, NE, USA d Landcare Research, P.O. Box 69, Lincoln, New Zealand Received 4 November 1997; accepted 7 January 1999
Abstract
A greater than 200 ka record of marine transgressions and regressions is recorded from a 75 m core from Banks Peninsula, Canterbury, New Zealand. This record comprises thick suites of muddy sediments attributed to back barrier, lake and lagoonal environments alternating with thin soil and loess complexes. These deposits have been dated using radio- carbon and thermoluminescence (TL) techniques supported by proxy data (diatoms, phytoliths, pollen and sedimentology). The aqueous deposits are attributed to three interglacials and an interstadial (Marine Isotope Stages 1, 5a, 5c, 6, and 7). The loesses and paleosols date to the intervening stadials (Isotope Stages 2, 5d (or 6?) and probably 8). On the basis of transgressive beach facies, back barrier swamps and barrier-blocked lake deposits, a partial sea-level curve including data from Isotope Stage 5 is presented. Our data indicate that Banks Peninsula has been tectonically stable over that period and we provide sea-level points that support the existing isotope curve during Stages 5 and 6. Detailed diatom records are limited to Isotope Stage 1 and the latter part of Stage 5. Diatom histories recorded from these stages are remarkably consistent. Both indicate a progressive ¯oral change from marine types through freshwater colonising species to freshwater planktonic assemblages. These re¯ect parallel histories of coastal evolution during the two interglacials. In both cases, marine transgression in the early part of the isotope phase was followed by lagoon development implying that a gravel spit extended across the embayment from the west. This was succeeded by lake development when the lagoon was cut off by the juncture of the spit with Banks Peninsula. This lake deepened as the coast rotated into swash alignment and the spit was converted into a gravel barrier. The vegetation history of the site indicates that mixed podocarp broadleaf forests, similar to the pre-European ¯ora of Banks Peninsula, occupied the region during Isotope Stages 1 and 7. This contrasts with the palynological interpretation of a marine record (DSDP Site 594) from off the Canterbury coast which suggested that Isotope Stage 7 was markedly cooler than the Holocene. During glacial periods, forest was eliminated and replaced by a tall shrubland of mixed montane and coastal af®nities. 1999 Elsevier Science B.V. All rights reserved.
Keywords: late Quaternary; coastal evolution; sea-level change; New Zealand; micropalaeontology; thermoluminescence dating
Ł Corresponding author. Tel.: C64-4-495-5233; Fax: C64-4-495-5186; E-mail: [email protected]
0031-0182/99/$ ± see front matter 1999 Elsevier Science B.V. All rights reserved. PII: S0031-0182(99)00035-8 102 J. Shulmeister et al. / Palaeogeography, Palaeoclimatology, Palaeoecology 152 (1999) 101±127
1. Introduction marine record to terrestrial New Zealand is fraught with potential pitfalls and good terrestrial records Despite an extensive history of Quaternary re- are required to verify the marine story. Terrestrial search in New Zealand, there is a paucity of long, records which can be correlated directly to the late Quaternary records. This re¯ects both the highly oceanic record are particularly valuable. Estuarine dynamic New Zealand landscape, where most de- and lagoonal systems provide the possibility of such positional environments are ephemeral on glacial± a correlation through direct evidence of marine trans- interglacial timeframes, and a tendency for research gressions and regressions. Barrier-blocked lagoonal to focus on the glacial record in mountainous ar- systems are potentially very useful sites for sea-level eas, where long sedimentary records are rare. In the investigations, as the fresh to brackish water environ- absence of these records, knowledge of the overall ments behind the barrier are sea-level-controlled and late Quaternary history of terrestrial New Zealand largely depositional. is remarkably fragmentary. By contrast, considerable The Kaitorete `Spit'±Lake Ellesmere complex is progress has been made in deciphering climate his- a large barrier-blocked lagoon system that occurs on tories from the marine record in the New Zealand the southern ¯ank of Banks Peninsula, Canterbury region (e.g. Nelson et al., 1985; Heusser and van (Fig. 1). Water levels in Lake Ellesmere are currently de Geer, 1994). These records con®rm the general managed to prevent inundation of agricultural land correspondence of glacial±interglacial cycles in New but in historical times the lake extended into the val- Zealand with the Northern Hemisphere pattern but leys around the southwestern ¯ank of the peninsula. are of low temporal resolution. Direct transfer of the These valleys are steep sided and narrow mouthed
Fig. 1. Location map showing Banks Peninsula and the location of the Gebbies Valley drill hole. The drill site is marked with an asterisk. J. Shulmeister et al. / Palaeogeography, Palaeoclimatology, Palaeoecology 152 (1999) 101±127 103 and, even on open marine sections of the Banks during the early history of Banks Peninsula. Well- Peninsula coast, provide loci for ®ne sediment depo- man (1979) suggested that subsidence occurred at a sition. Thick sediment ®lls exist in all these valleys. rate of 0.2 mm=year, and about 600 m of sediment A drilling programme in Gebbies Valley along the is recorded in the topographic low on the adjacent southwest ¯ank of peninsula, behind the Kaitorete Canterbury Plains, now occupied by Lake Ellesmere `Spit', recovered 75 m of largely ®ne-grained sed- (Talbot et al., 1986). However, Lawrie (1993) points iments. These sediments extend through the last to the presence of buried shore platforms at C5to >200 ka and are believed to represent three glacial± C8 m AMSL (above mean sea level), of probable interglacial cycles, and record several marine trans- Last Interglacial age, around the southwestern side gressions and regressions. In this paper, we present of the peninsula and similar features have now been an environmental history and the coastal evolution identi®ed elsewhere on the peninsula (Bal, 1997). along the southwestern ¯ank of Banks Peninsula for These imply local tectonic stability for at least the the last 200 ka. last 125,000 years.
1.1. Physiographic and geological setting 1.2. General late Quaternary history
Banks Peninsula is the remnant of a complex Two major types of Quaternary sediments are of Miocene shield volcanos (Sewell, 1988). It is recorded on Banks Peninsula: loess and valley ®lls. located along the eastern seaboard of the South Is- Loess is widespread, characteristically forming a thin land, New Zealand, immediately south of the city cover on upper valley sides and spurs, but increas- of Christchurch (Fig. 1). The volcanic complex is ing to several metres in depth on lower slopes. It is largely basaltic and contains two large calderas, now derived mainly from the Southern Alps greywackes, occupied by Akaroa and Lyttelton Harbours. The but the occurrence of sponge spicules in deposits complex has been heavily eroded and steep-sided on the eastern side of the peninsula indicates a valleys have been incised and=or exhumed along the subsidiary source on the exposed former sea¯oor outer ¯anks of the peninsula. These valleys proba- to the south and east, during low sea-stands (Rae- bly ®rst developed during the early history of the side, 1964). Around the southwestern slopes of the complex (Shelley, 1989), and the two calderas were peninsula loess overlies wave-cut platforms at C5 formed before the deposition of the youngest lava m AMSL, indicating a pre-Holocene age, and these ¯ows 7±5.8 million years ago (Sewell et al., 1992). are interpreted as Last Interglacial by Lawrie (1993). Initially an offshore island, Banks Peninsula has At least three loess layers separated by colluvium been alternately attached to, or separated from, the and palaeosols have been identi®ed in this area and mainland of the South Island of New Zealand, in Akaroa Harbour (Grif®ths, 1973). A C-14 age of depending on interglacial±glacial variations in sea 17,450 š 2070 yr BP was obtained for the upper- level, and on the extent to which the alluvial fans most layer in the Barry's Bay site in Akaroa Harbour of the Canterbury Plains have been built out from (Grif®ths, 1973), but Goh et al. (1977) demonstrated the eastern front of the Southern Alps (Brown and that this and other ages which they obtained from Wilson, 1988). These fans are the source of material untreated material in the palaeosols gave minimum for the construction of barrier beaches, dunes and ages. More recent studies yielded ages of over 140 ka lagoons along their seaward edge which have linked in the lower loess units (e.g. Berger et al., 1994a) and the peninsula to the mainland. it is clear that many of the shore platforms around On the southwest side of the peninsula, from the southwest ¯ank of Banks Peninsula are at least Christchurch to Lake Forsyth, basement rock is 80± Last Interglacial in age, but could be considerably 100 m below present sea level (Brown and Weeber, older. 1992, 1994), and the lower parts of the valleys are The valley ®lls on Banks Peninsula are domi- ®lled with sediment. The considerable depth of these nantly composed of muds and very ®ne sands (e.g. valleys below present sea level may be related to in- Dingwall, 1974; Soons et al., 1997). This study is the cision during low glacial sea levels, or to subsidence ®rst detailed investigation of a valley ®ll complex. 104 J. Shulmeister et al. / Palaeogeography, Palaeoclimatology, Palaeoecology 152 (1999) 101±127
2. Methods lar sampling at about 1 m intervals was undertaken. Foraminifera were found in only one sample and 2.1. Core and sample recovery work in this area is unlikely to be pursued. Sam- ples were examined at less than one metre intervals The core was located at the mouth of Gebbies for both phytoliths and diatoms. Diatom preserva- Valley on a relict foredune system of historical Lake tion through the core was intermittent. Zones of Ellesmere (see Fig. 1). The hole was cored using a good preservation occurred between 7andC1.5 truck-mounted, cable tool, water well drilling sys- m AMSL and between 21.5 and 15.5 m AMSL. tem. The hole was cased using standard 8 inch water Some preservation also occurred between 63.5 and well steel casing, which was advanced behind the 57.5 m AMSL. From the initial examination, 50 corer every 3 m. This eliminated signi®cant down- samples were counted from these zones. A mini- hole contamination. The top 42 m of material was mum of 300 diatoms per level were counted. Three cored using a 1 m long, 75 mm internal diameter intervals of phytolith preservation were recognised push corer mounted on the drilling bit. Below 42 roughly matching those for diatoms. m, samples were recovered using a 100 mm internal Initial samples for pollen analyses were recovered diameter hollow ¯ight piston corer. These samples at 1 m intervals but most proved barren. Subse- were recovered in 0.25±0.50 m increments and came quently, samples for pollen analysis were collected out as cohesive core except in loose gravels near the at irregular intervals from those sections of the core bottom of the hole. that offered the greatest concentration of organic material. Samples were obtained at relatively close 2.2. Sedimentary descriptions intervals between 44.77 and 50.25 m. It is certain that important changes have been missed due to non- All samples were described visually either in the preservation of pollen. Damaged and unidenti®able ®eld (bagged material) or the laboratory (core sec- pollen were found, and are believed to be reworked tion). Attributes such as the presence of worm casts, from the Late Cretaceous Gebbies Pass Plant Beds, organics, concretions and coarse lithics were record- presently exposed at the head of the valley (Sewell et ed. Munsell colours were determined for all samples. al., 1992). Depending on the visible grain size distribution, 20 g to about 250 g of material was extracted for 2.4. Luminescence dating particle size analyses, using a sedigraph for the ®ner than 63 µm fraction and an automated settling tube In ideal circumstances, sediments up to at least (Rapid Sediment Analyser) for the coarser material. 800 ka old can now be aged by speci®c thermolu- The results were combined and composite grain size minescence (TL) procedures (Berger et al., 1994b; distributions and summary statistics using Folk and Berger, 1995). Therefore we chose twelve samples Ward (1957) were calculated. Grain size information for TL dating, representing most of the recovered is integrated into Fig. 2. core length. These samples were extracted from the Based on the gross characteristics, informal sed- centre of the split cores (above 39 m AMSL) or the imentary units were identi®ed. Sixteen subsamples interior of large cohesive samples recovered during from these sedimentary units were recovered for XRD auguring (below 39 m AMSL). Parallel samples (X-ray diffraction) analyses of both the clay and, were recovered for dose-rate determination. Thick- where appropriate, sand fractions. Samples for lumi- source alpha-particle counting was used to measure nescence dating were selected on basis of sediment U and Th concentrations, while K concentration was homogeneity and inferred origins (e.g. Berger, 1988). determined from atomic absorption spectrophotom- etry. Cosmic-ray dose rates were estimated and the 2.3. Micropalaeontology average water content was assumed to be equiva- lent to the present saturation values for the samples. Samples were recovered for pollen, phytolith, di- Alpha ef®ciency values (b-values) were determined atom, and foraminiferal analyses. Initially sub-regu- from a comparison of extrapolated alpha-dose re- J. Shulmeister et al. / Palaeogeography, Palaeoclimatology, Palaeoecology 152 (1999) 101±127 105
Fig. 2. Stratigraphy of the Gebbies Valley Borehole (G2). 106 J. Shulmeister et al. / Palaeogeography, Palaeoclimatology, Palaeoecology 152 (1999) 101±127 0.96 .For 1.7 2.9 9.1 ¦ except š 19 12) 21 13 12 11 13 10 31 33 37 š š š š š š š š š š š š š š 21, 25, and 27 TL age ., 1994a). These g TB1 TB2 PB 111 TB1 2.8 PB PB PB PB PB PB PB ADD (111 PB PB 9.18 TB2TB2 93.5 154 TB2 114 ADD 115 ADD 99 TB2 119 TB2 136 TB2 180 PB 172 TB2 204 PB 24.5 ======Temperature 2.7 250±340 3.0 220±280 5.4 230±350 3.3 250±320 9.6 230-310 64 220±340 69 250±340 130 260±320 59 230±340 58 240±350 66 220±320 34 250±340 160 210±300 40 230±300 200 250±350 32 240±320 68 250±320 41 230±300 40 220±310 40 230±300 41 220±380 25 240±320 91 200±290 110 210±300 140 250±380 š š š š š š š š š š š š š š š š š š š š š š š š š f E D 28.5 38.1 80.0 400 357 570 407 417 455 348 590 830 e throughout table. ¦ Preheat 1 š 0.096 120 32.6 0.18) 417 0.13 130 85.8 0.16 1400.23 1350.22 365 140 571 434 0.17 140 413 0.22 140 417 0.17 135 429 0.12 140 378 0.17 130 555 0.13 130 560 0.38 140 880 š š š š š š š š š š š š š ka) (ëC) (Gy) (Gy) (ka) = R D )(Gy 2 0.11 3.552 0.11 3.50 0.17 3.90 0.23 3.58 0.26 4.20 0.20 3.60 0.16 2.79 0.21 3.09 0.14 3.28 0.48 4.31 0.03 for others. Errors are d 0.2) (3.74 0.30.3 3.71 3.82 š š š š š š š š š š š š š š -value b (1.0 0.73 1.04 0.92 1.07 1.1 0.73 1.2 1.1 0.98 1.47 1.1 0.71 1.1 1.46 1.1 0.84 1.2 1.1 1.1 1.65 1.0 1.12 1.2 0.63 š š š š š š š š š š š š 0.02 for ®rst sample and š 0.34 9.8 0.28 9.69 0.22 10.26 0.33 12.8 0.35 10.5 0.37 8.4 0.34 9.7 0.30 7.12 0.31 6.0 0.35 10.0 0.37 8.8 0.34 11.4 š š š š š š š š š š š š plot. The lower limit is 10ëC below lowest data point because each point represents a 10ëC slice of the readout curves. The value having the lowest error was chosen, when comparing PB and either TB2 or ADD results, when these results agree with 1 UTh T E ± D E c D O 2 K (%) (ppm) (ppm) (pGy cm b mass dry sample. Uncertainty is = 0.05. Water š a 11.48 0.42 3.32 2.18 17.38 0.24 2.28 2.59 19.06 0.4421.90 0.38 3.10 3.06 3.28 3.21 32.68 0.34 3.01 2.82 27.02 0.30 2.6733.28 3.38 0.29 2.78 3.00 ======42.47 0.28 2.18 3.48 44.89 0.34 1.91 3.10 55.62 0.37 2.28 3.30 56.87 0.34 2.94 2.72 65.62 0.35 2.78 3.79 = = = = = -measuring method is indicated after the slash. Partial bleach (PB), total-bleach1 (TB1), total-bleach2 (TB2) and additive dose (ADD) (Berger et al E Ratio mass water Alpha ef®ciency factor. This was not measured for samples 15 and 13, butSelected only estimated temperature from range measured range in for other samples. After slash: depth (m) within core. Uncertainty is For 4 days. Weighted (by inverse variance) mean and average error over temperature range in next column. D methods are illustrated inand the ®gures. TB1 TB2 used used thesample bleaching 435±700 15. for range For 24 for age h these calculation, with samples. the a PB full-spectrum used Hg-vapour the lamp. 435±700 PB nm used range 550±700 for nm the wavelengths for other samples samples, 3, and 5, TB1 used the range 550±700 nm for all of the others Table 1 TL dating results Sample BP-27 Please note thata 4.5 m should be addedb to all core depthc data to convert tod m AMSL. e f g sample 27, the TB2 result is ignored (see Fig. 3C). BP-25 BP-21 BP-23 BP-13 BP-15 BP-11 BP-9 BP-7 BP-5 BP-3 BP-1 J. Shulmeister et al. / Palaeogeography, Palaeoclimatology, Palaeoecology 152 (1999) 101±127 107 sponse curves with beta-dose response curves (see 3.1. TL dating Berger, 1988). Alpha doses were administered by sealed-foil 241Am sources. The above data were used Age determinations from TL are presented in to calculate a dose rate .Dr/ (equations in Berger, Table 1. The TL results indicate that about 12.5 1988). m of Holocene muds and silts overlie a yellowish A luminescence age is equal to the equivalent silt (loess) of Last Glacial Maximum age (BP-25: dose .DE/ divided by the dose rate Dr. DE values 22:8 š 1:7 ka) (Fig. 3). The major feature of the were determined from extrapolations of beta-dose TL results is the concentration of ages around 90± response curves. Beta doses were administered by 130 ka indicating a thick sequence of deposits of a 125 mCi sealed-foil 90Sr±90Y source. The partial- Marine Isotope Stage (MIS) 5 (and possibly 6) be- bleach and total-bleach procedures (Aitken, 1985; tween 12.5 and 43.75 m. Of these ages, BP-15 at Berger, 1988) were used to measure DE values, 22.52 m AMSL is anomalous because of its larger which were then plotted against readout temperature age and greater error range (though it still overlaps to provide a plateau test. To remove thermally unsta- surrounding samples at 2¦). The basal three TL sam- ble and anomalous-fading TL components from the ples between 51 and 61 m yield ages averaging laboratory-irradiated subsamples prior to TL read- around 190 ka and may represent either MIS 7 or out, elevated-temperature storage treatments were MIS7andMIS8deposits. employed (see refs. in Berger, 1994). The chosen preheating temperatures are listed 3.2. Radiocarbon dating in Table 1. Subsamples of polymineralic grains 4± 11 µm in diameter were heated at 5ëC=s. The TL Radiocarbon dating results are presented in Ta- was recorded through a 5-mm-thick KOPP optical ble 2. Only one completely satisfactory age (NZA glass ®lter CS-560, which transmits 360±495 nm 5262) was obtained from a wood sample recovered wavelengths at 5% cut (345±505 nm at 1% cut). at a depth of 7.08 m AMSL. This yielded an early Thus most ultraviolet emissions are blocked, as are Holocene age of 7558 š 89 BP (8271 š 89 CAL BP: those from several plagioclase feldspars, and the Stuiver and Reimer (1993)) which is consistent with detected signal is dominated by the ca. 410 nm the results from the luminescence dating. Two of the emissions common to most K-rich feldspars (e.g. remaining samples yielded ®nite ages. NZA 5264 Berger, 1995). ( 21.54 m AMSL) yielded an age of 43,900 š 1700 BP. This is at the limit of the radiocarbon technique 2.5. Radiocarbon (the sample has only 0.42% of the modern C-14 content) and is considerably younger than the TL Six samples were submitted for AMS dating at ages from these levels. The age of 37,230 š 620 the Institute of Geological and Nuclear Sciences BP derived from NZA 5663 at 54.75 m AMSL is laboratory in Wellington, New Zealand. clearly erroneous as there are three overlying radio- carbon samples which yield minimum ages that are signi®cantly older and occur up to 30 m higher in the 3. Results core. In addition, two consistent TL (BP-3 and BP-5) ages of around 180 ka were recorded from within a A summary stratigraphic column is presented in few metres of this sample. Fig. 2. Of the 75 m of material recovered, most (about 60 m) are ®ne silts and clays. Gravels are 3.3. Diatoms common only in the basal 10 m. In the clay frac- tion, XRD indicates that kaolinite, mica and chlorite Over 130 species of diatoms were identi®ed from are ubiquitous with smectite and vermiculite locally samples in the core. Of these about 30 species oc- present. curred commonly or in signi®cant numbers in sec- tions of the core. They were preserved in three zones. The lowest of the zones occurs from 63.5 to 108 J. Shulmeister et al. / Palaeogeography, Palaeoclimatology, Palaeoecology 152 (1999) 101±127
Fig. 3. Luminescence growth curves for samples BP-25 (A, B), BP-27 (C) and BP-21 (D).