skip to main content


Title: Discovery of the Youngest Toba Tuff in South Africa; Challenges of Extremely Low Abundance Cryptotephra
The 74 ka Youngest Toba Tuff (YTT) was discovered as cryptotephra in South African archaeological sites at Pinnacle Point (PP) 5-6N, Vleesbaai [1] and Klasies River on the Indian Ocean, and the Diepkloof Rock Shelter on the Atlantic coast nearly 750 km west of PP. The YTT eruption distributed tephra across eastern and southern Africa and provides a widespread isochron useful for dating archaeological deposits, testing age models, and precisely determining the timing of changes in human behavior. At PP, we demonstrated that the MIS 4-5 transition began just before the YTT eruption. Humans thrived both through the YTT event and the changing climate, and important changes in technology occurred just after the Toba eruption [1]. Controversy related to trapped charge age models at the Diepkloof rock shelter [2,3] were resolved by identifying YTT at a location in the stratigraphic section that confirmed the Jacobs et al. [2] model for the site, confirming that technological changes similar to those observed at PP occurred synchronously to those at Diepkloof, not substantially before as suggested by a prior published age model [3]. Processing samples with very low abundance cryptotephra, such as that found in South Africa, is a challenge and requires revision of standard laboratory techniques. Samples with high organic or clay content benefit by being treated with 10% HCl and 3% hydrogen peroxide (H2O2). This step does not degrade shard integrity or affect chemistry and is useful in separating shards from clay and organic particles allowing better recovery in heavy liquids. We also modified the heavy liquid density range from 1.95 - 2.55 to 2.2 - 2.5 g/cm3 to more effectively remove quartz and feldspar as well as biogenic silica. This density range captures shards ranging from rhyolite to dacite in composition. [1] Smith, E. I. et al. Nature 555, 511, 2018 [2] Jacobs, Z. & Roberts, R. G. Journal of Archaeological Science 63, 175-192, 2015. [3] Tribolo, C. et al. Journal of Archaeological Science 40, 3401-3411 2013.  more » « less
Award ID(s):
1917173
NSF-PAR ID:
10278334
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Annual VM Goldschmidt Conference
ISSN:
1945-7022
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Site description. This data package consists of data obtained from sampling surface soil (the 0-7.6 cm depth profile) in black mangrove (Avicennia germinans) dominated forest and black needlerush (Juncus roemerianus) saltmarsh along the Gulf of Mexico coastline in peninsular west-central Florida, USA. This location has a subtropical climate with mean daily temperatures ranging from 15.4 °C in January to 27.8 °C in August, and annual precipitation of 1336 mm. Precipitation falls as rain primarily between June and September. Tides are semi-diurnal, with 0.57 m median amplitudes during the year preceding sampling (U.S. NOAA National Ocean Service, Clearwater Beach, Florida, station 8726724). Sea-level rise is 4.0 ± 0.6 mm per year (1973-2020 trend, mean ± 95 % confidence interval, NOAA NOS Clearwater Beach station). The A. germinans mangrove zone is either adjacent to water or fringed on the seaward side by a narrow band of red mangrove (Rhizophora mangle). A near-monoculture of J. roemerianus is often adjacent to and immediately landward of the A. germinans zone. The transition from the mangrove to the J. roemerianus zone is variable in our study area. An abrupt edge between closed-canopy mangrove and J. roemerianus monoculture may extend for up to several hundred meters in some locations, while other stretches of ecotone present a gradual transition where smaller, widely spaced trees are interspersed into the herbaceous marsh. Juncus roemerianus then extends landward to a high marsh patchwork of succulent halophytes (including Salicornia bigellovi, Sesuvium sp., and Batis maritima), scattered dwarf mangrove, and salt pans, followed in turn by upland vegetation that includes Pinus sp. and Serenoa repens. Field design and sample collection. We established three study sites spaced at approximately 5 km intervals along the western coastline of the central Florida peninsula. The sites consisted of the Salt Springs (28.3298°, -82.7274°), Energy Marine Center (28.2903°, -82.7278°), and Green Key (28.2530°, -82.7496°) sites on the Gulf of Mexico coastline in Pasco County, Florida, USA. At each site, we established three plot pairs, each consisting of one saltmarsh plot and one mangrove plot. Plots were 50 m^2 in size. Plots pairs within a site were separated by 230-1070 m, and the mangrove and saltmarsh plots composing a pair were 70-170 m apart. All plot pairs consisted of directly adjacent patches of mangrove forest and J. roemerianus saltmarsh, with the mangrove forests exhibiting a closed canopy and a tree architecture (height 4-6 m, crown width 1.5-3 m). Mangrove plots were located at approximately the midpoint between the seaward edge (water-mangrove interface) and landward edge (mangrove-marsh interface) of the mangrove zone. Saltmarsh plots were located 20-25 m away from any mangrove trees and into the J. roemerianus zone (i.e., landward from the mangrove-marsh interface). Plot pairs were coarsely similar in geomorphic setting, as all were located on the Gulf of Mexico coastline, rather than within major sheltering formations like Tampa Bay, and all plot pairs fit the tide-dominated domain of the Woodroffe classification (Woodroffe, 2002, "Coasts: Form, Process and Evolution", Cambridge University Press), given their conspicuous semi-diurnal tides. There was nevertheless some geomorphic variation, as some plot pairs were directly open to the Gulf of Mexico while others sat behind keys and spits or along small tidal creeks. Our use of a plot-pair approach is intended to control for this geomorphic variation. Plot center elevations (cm above mean sea level, NAVD 88) were estimated by overlaying the plot locations determined with a global positioning system (Garmin GPS 60, Olathe, KS, USA) on a LiDAR-derived bare-earth digital elevation model (Dewberry, Inc., 2019). The digital elevation model had a vertical accuracy of ± 10 cm (95 % CI) and a horizontal accuracy of ± 116 cm (95 % CI). Soil samples were collected via coring at low tide in June 2011. From each plot, we collected a composite soil sample consisting of three discrete 5.1 cm diameter soil cores taken at equidistant points to 7.6 cm depth. Cores were taken by tapping a sleeve into the soil until its top was flush with the soil surface, sliding a hand under the core, and lifting it up. Cores were then capped and transferred on ice to our laboratory at the University of South Florida (Tampa, Florida, USA), where they were combined in plastic zipper bags, and homogenized by hand into plot-level composite samples on the day they were collected. A damp soil subsample was immediately taken from each composite sample to initiate 1 y incubations for determination of active C and N (see below). The remainder of each composite sample was then placed in a drying oven (60 °C) for 1 week with frequent mixing of the soil to prevent aggregation and liberate water. Organic wetland soils are sometimes dried at 70 °C, however high drying temperatures can volatilize non-water liquids and oxidize and decompose organic matter, so 50 °C is also a common drying temperature for organic soils (Gardner 1986, "Methods of Soil Analysis: Part 1", Soil Science Society of America); we accordingly chose 60 °C as a compromise between sufficient water removal and avoidance of non-water mass loss. Bulk density was determined as soil dry mass per core volume (adding back the dry mass equivalent of the damp subsample removed prior to drying). Dried subsamples were obtained for determination of soil organic matter (SOM), mineral texture composition, and extractable and total carbon (C) and nitrogen (N) within the following week. Sample analyses. A dried subsample was apportioned from each composite sample to determine SOM as mass loss on ignition at 550 °C for 4 h. After organic matter was removed from soil via ignition, mineral particle size composition was determined using a combination of wet sieving and density separation in 49 mM (3 %) sodium hexametaphosphate ((NaPO_3)_6) following procedures in Kettler et al. (2001, Soil Science Society of America Journal 65, 849-852). The percentage of dry soil mass composed of silt and clay particles (hereafter, fines) was calculated as the mass lost from dispersed mineral soil after sieving (0.053 mm mesh sieve). Fines could have been slightly underestimated if any clay particles were burned off during the preceding ignition of soil. An additional subsample was taken from each composite sample to determine extractable N and organic C concentrations via 0.5 M potassium sulfate (K_2SO_4) extractions. We combined soil and extractant (ratio of 1 g dry soil:5 mL extractant) in plastic bottles, reciprocally shook the slurry for 1 h at 120 rpm, and then gravity filtered it through Fisher G6 (1.6 μm pore size) glass fiber filters, followed by colorimetric detection of nitrite (NO_2^-) + nitrate (NO_3^-) and ammonium (NH_4^+) in the filtrate (Hood Nowotny et al., 2010,Soil Science Society of America Journal 74, 1018-1027) using a microplate spectrophotometer (Biotek Epoch, Winooski, VT, USA). Filtrate was also analyzed for dissolved organic C (referred to hereafter as extractable organic C) and total dissolved N via combustion and oxidation followed by detection of the evolved CO_2 and N oxide gases on a Formacs HT TOC/TN analyzer (Skalar, Breda, The Netherlands). Extractable organic N was then computed as total dissolved N in filtrate minus extractable mineral N (itself the sum of extractable NH_4-N and NO_2-N + NO_3-N). We determined soil total C and N from dried, milled subsamples subjected to elemental analysis (ECS 4010, Costech, Inc., Valencia, CA, USA) at the University of South Florida Stable Isotope Laboratory. Median concentration of inorganic C in unvegetated surface soil at our sites is 0.5 % of soil mass (Anderson, 2019, Univ. of South Florida M.S. thesis via methods in Wang et al., 2011, Environmental Monitoring and Assessment 174, 241-257). Inorganic C concentrations are likely even lower in our samples from under vegetation, where organic matter would dilute the contribution of inorganic C to soil mass. Nevertheless, the presence of a small inorganic C pool in our soils may be counted in the total C values we report. Extractable organic C is necessarily of organic C origin given the method (sparging with HCl) used in detection. Active C and N represent the fractions of organic C and N that are mineralizable by soil microorganisms under aerobic conditions in long-term soil incubations. To quantify active C and N, 60 g of field-moist soil were apportioned from each composite sample, placed in a filtration apparatus, and incubated in the dark at 25 °C and field capacity moisture for 365 d (as in Lewis et al., 2014, Ecosphere 5, art59). Moisture levels were maintained by frequently weighing incubated soil and wetting them up to target mass. Daily CO_2 flux was quantified on 29 occasions at 0.5-3 week intervals during the incubation period (with shorter intervals earlier in the incubation), and these per day flux rates were integrated over the 365 d period to compute an estimate of active C. Observations of per day flux were made by sealing samples overnight in airtight chambers fitted with septa and quantifying headspace CO_2 accumulation by injecting headspace samples (obtained through the septa via needle and syringe) into an infrared gas analyzer (PP Systems EGM 4, Amesbury, MA, USA). To estimate active N, each incubated sample was leached with a C and N free, 35 psu solution containing micronutrients (Nadelhoffer, 1990, Soil Science Society of America Journal 54, 411-415) on 19 occasions at increasing 1-6 week intervals during the 365 d incubation, and then extracted in 0.5 M K_2SO_4 at the end of the incubation in order to remove any residual mineral N. Active N was then quantified as the total mass of mineral N leached and extracted. Mineral N in leached and extracted solutions was detected as NH_4-N and NO_2-N + NO_3-N via colorimetry as above. This incubation technique precludes new C and N inputs and persistently leaches mineral N, forcing microorganisms to meet demand by mineralizing existing pools, and thereby directly assays the potential activity of soil organic C and N pools present at the time of soil sampling. Because this analysis commences with disrupting soil physical structure, it is biased toward higher estimates of active fractions. Calculations. Non-mobile C and N fractions were computed as total C and N concentrations minus the extractable and active fractions of each element. This data package reports surface-soil constituents (moisture, fines, SOM, and C and N pools and fractions) in both gravimetric units (mass constituent / mass soil) and areal units (mass constituent / soil surface area integrated through 7.6 cm soil depth, the depth of sampling). Areal concentrations were computed as X × D × 7.6, where X is the gravimetric concentration of a soil constituent, D is soil bulk density (g dry soil / cm^3), and 7.6 is the sampling depth in cm. 
    more » « less
  2. IODP Expedition 379 to the Amundsen Sea continental rise recovered latest Miocene-Holocene sediments from two sites on a drift in water depths >3900m. Sediments that are dominated by clay and silty clay host pebbles and cobbles of ice-rafted detritus (IRD) (Gohl et al. 2021, doi: 10.14379/iodp.proc.379.2021). Cobble-sized dropstones also appear as fall-in, in the top of cores recovered from sediments >5.3 Ma. The principle means to deposit abundant IRD and sparse dropstones in deep sea sediment is through melting of icebergs released by Antarctic ice-sheet calving events. We use petrological and age characteristics of clasts from the Exp379 sites to fingerprint their bedrock provenance to extend knowledge of subglacial bedrock, and with the intention to illuminate changes in icesheet extent between 7 – 3 Ma that lend credence to forecasts of extensive future change. Mapped onshore geology shows pronounced distinctions in bedrock age between tectonic provinces of West or East Antarctica (e.g. Cox et al. 2020, doi:10.21420/7SH7-6K05; Jordan et al. 2020, doi.org/10.1038/s43017-019-0013-6). This allows us to use geochronology and thermochronology of rock clasts and minerals for tracing their provenance, and hence ascertain whether IRD deposited at the 379 drillsites originated from proximal or distal Antarctic sources. We here report zircon and apatite U-Pb dates from several sand samples and dropstones taken from latest Miocene, early Pliocene, and Plio-Pleistocene-boundary sediments. Additional Hf isotope data, and apatite fission track and 40Ar/39Ar Kfeldspar ages for some of the same samples help to strengthen provenance interpretations. The study revealed three distinct zircon age populations at ca. 100, 175, and 250 Ma. Using Kolmogorov-Smirnov (K-S) statistical tests to compare our new igneous and detrital zircon (DZ) U-Pb results with previously published data, we found strong similarities to West Antarctic bedrock, but low correspondence to prospective sources in East Antarctica, implying a role for icebergs calved from the West Antarctic Ice Sheet (WAIS). The ~100 Ma age resembles plutonic ages from Marie Byrd Land and islands in Pine Island Bay. The ~250 and 175 Ma populations match published mineral dates from shelf sediments in the eastern Amundsen Sea Embayment as well as granite ages from the Antarctic Peninsula and the Ellsworth-Whitmore Mountains (EWM). The different derivation of coarse sediment sources requires changes in iceberg origin through the latest Miocene, early Pliocene, and Plio/Pleistocene, likely the result of changes in WAIS extent. One unique dropstone recovered from Exp379 Site U1533B is green quartz arenite, which yielded mostly 500-625 Ma detrital zircons. In visual appearance and dominant U-Pb age population, it resembles a sandstone dropstone recovered from Exp382 Site U1536 in the Scotia Sea (Hemming et al. 2020, https://gsa.confex.com/gsa/2020AM/meetingapp.cgi/Paper/357276). K-S tests yield high values (P ≥ 0.6), suggesting a common provenance for both dropstones recovered from late Miocene to Pliocene sediments, despite the 3270 km distance separating the sites. Comparisons to published data, in progress, narrow the group of potential on-land sources to exposures in the EWM or isolated ranges at far south latitudes in the Antarctic interior. If both dropstones originated from the same source area, they signify that dramatic shifts in the WAIS grounding line position do occur, along with periodic opening of a seaway connecting the Amundsen and Weddell Seas. 
    more » « less
  3. Abstract Caves and rock shelters contribute important records to local, regional and sub-continental reconstructions of environment and climate change through the southern African Quaternary. Against a backdrop of pronounced climate change, the archaeological record of the Marine Isotope Stage 6 to 1 period in southern Africa documents a remarkable time in the behavioural and technological evolution of anatomically modern humans. Significant evidence of this evolution is represented in diverse components of the sedimentary record in caves and rock shelters in the region. We present a catalogue of published caves and rock shelters in southern Africa that preserve temporally-relevant clastic and chemical palaeoclimatic proxies in order to: (1) facilitate the integration of cave and rock shelter sedimentary data into broader, regional chronostratigraphically-correlated palaeoclimatic sequences; and (2) identify possible areas and proxies that require focused research in the future. To demonstrate the complexity of the Marine Isotope Stage 6 to 1 stratigraphic record and use of palaeoenvironmental proxies, we present three case studies representing interior and coastal contexts: Border Cave, Klasies River Mouth and Pinnacle Point. These examples aptly demonstrate the challenges of these contexts, but also the opportunities for palaeoenvironmental research in southern Africa when conducted through integrated, multidisciplinary approaches. Published records of palaeoenvironmental research from cave and rock shelter sequences in southern Africa are heavily biased to the South African coastal areas and the record is temporally and spatially fragmented. However, there are interesting patterns in the chronostratigraphic record and in the distribution of sites within the context of the geology and vegetation ecology of southern Africa that require further exploration. There are also promising techniques in stable isotope analysis that can be applied to abundant sedimentary components found in the region’s caves and rock shelters, and in its museums. 
    more » « less
  4. The Toba eruption ∼74,000 y ago was the largest volcanic eruption since the start of the Pleistocene and represents an important test case for understanding the effects of large explosive eruptions on climate and ecosystems. However, the magnitude and repercussions of climatic changes driven by the eruption are strongly debated. High-resolution paleoclimate and archaeological records from Africa find little evidence for the disruption of climate or human activity in the wake of the eruption in contrast with a controversial link with a bottleneck in human evolution and climate model simulations predicting strong volcanic cooling for up to a decade after a Toba-scale eruption. Here, we use a large ensemble of high-resolution Community Earth System Model (CESM1.3) simulations to reconcile climate model predictions with paleoclimate records, accounting for uncertainties in the magnitude of Toba sulfur emissions with high and low emission scenarios. We find a near-zero probability of annual mean surface temperature anomalies exceeding 4 °C in most of Africa in contrast with near 100% probabilities of cooling this severe in Asia and North America for the high sulfur emission case. The likelihood of strong decreases in precipitation is low in most of Africa. Therefore, even Toba sulfur release at the upper range of plausible estimates remains consistent with the muted response in Africa indicated by paleoclimate proxies. Our results provide a probabilistic view of the uneven patterns of volcanic climate disruption during a crucial interval in human evolution, with implications for understanding the range of environmental impacts from past and future supereruptions.

     
    more » « less
  5. IODP Expedition 379 to the Amundsen Sea continental rise recovered latest Miocene-Holocene sediments from two sites on a drift in water depths >3900m. Sediments are dominated by clay and silty clay with coarser-grained intervals and ice-rafted detritus (IRD) (Gohl et al. 2021, doi:10.14379/iodp.proc.379.2021). Cobble-sized dropstones appear as fall-in, in cores recovered from sediments >5.3 Ma. We consider that abundant IRD and the sparse dropstones melted out of icebergs formed due to Antarctic ice-sheet calving events. We are using petrological and age characteristics of the clasts from the Exp379 sites to fingerprint their bedrock provenance. The results may aid in reconstruction of past changes in icesheet extent and extend knowledge of subglacial bedrock. Mapped onshore geology shows pronounced distinctions in bedrock age between tectonic provinces of West or East Antarctica (e.g. Cox et al. 2020, doi:10.21420/7SH7-6K05; Jordan et al. 2020, doi.org/10.1038/s43017-019-0013-6). This allows us to use geochronology and thermochronology of rock clasts and minerals for tracing their provenance, and ascertain whether IRD deposited at IODP379 drillsites originated from proximal or distal Antarctic sources. We here report zircon and apatite U-Pb dates from four sand samples and five dropstones taken from latest Miocene, early Pliocene, and Plio-Pleistocene-boundary sediments. Additional Hf isotope data, and apatite fission track and 40Ar/39Ar Kfeldspar ages for some of the same samples help to strengthen provenance interpretations. The study revealed three distinct zircon age populations at ca. 100, 175, and 250 Ma. Using Kolmogorov-Smirnov (K-S) statistical tests to compare our new igneous and detrital zircon (DZ) U-Pb results with previously published data, we found strong similarities to West Antarctic bedrock, but low correspondence to prospective sources in East Antarctica, implying a role for icebergs calved from the West Antarctic Ice Sheet (WAIS). The ~100 Ma age resembles plutonic ages from Marie Byrd Land and islands in Pine Island Bay. The ~250 and 175 Ma populations match published mineral dates from shelf sediments in the eastern Amundsen Sea Embayment as well as granite ages from the Antarctic Peninsula and the Ellsworth-Whitmore Mountains (EWM). The different derivation of coarse sediment sources requires changes in iceberg origin through the latest Miocene, early Pliocene, and Plio/Pleistocene, likely the result of changes in WAIS extent. One unique dropstone recovered from Exp379 Site U1533B is green quartz arenite, which yielded mostly 500-625 Ma detrital zircons. In visual appearance and dominant U-Pb age population, it resembles a sandstone dropstone recovered from Exp382 Site U1536 in the Scotia Sea (Hemming et al. 2020, https://gsa.confex.com/gsa/2020AM/meetingapp.cgi/Paper/357276). K-S tests yield high values (P ≥ 0.6), suggesting a common provenance for both dropstones recovered from late Miocene to Pliocene sediments, despite the 3270 km distance separating the sites. Comparisons to published data, in progress, narrow the group of potential on-land sources to exposures in the EWM or isolated ranges at far south latitudes in the Antarctic interior. If both dropstones originated from the same source area, they could signify dramatic shifts in the WAIS grounding line position, and the possibility of the periodic opening of a seaway connecting the Amundsen and Weddell Seas. https://meetingorganizer.copernicus.org/EGU21/EGU21-9151.html 
    more » « less