ELA-ICP-MS U-Pb zircon geochronology of regional volcanism hosting the Bajo de la Alumbrera Cu-Au deposit: implications for porphyry-related mineralization

ELA-ICP-MS U-Pb zircon geochronology has been used to show that the porphyritic intrusions related to the formation of the Bajo de la Alumbrera porphyry Cu-Au deposit, NW Argentina, are cogenetic with stratigraphically well-constrained volcanic and volcaniclastic rocks of the Late Miocene Farallon Negro Volcanic Complex. Zircon geochronology for intrusions in this deposit and the host volcanic sequence show that multiple mineralized porphyries were emplaced in a volcanic complex that developed over 1.5 million years. Volcanism occurred in a multivent volcanic complex in a siliciclastic intermontane basin. The complex evolved from early mafic-intermediate effusive phases to a later silicic explosive phase associated with mafic intrusions. Zircons from the basal mafic-intermediate lavas have ages that range from 8.46 +/- 0.14 to 7.94 +/- 0.27 Ma. Regionally extensive silicic explosive volcanism occurred at similar to8.0 Ma (8.05 +/- 0.13 and 7.96 +/- 0.11 Ma), which is co-temporal with intrusion of the earliest mineralized porphyries at Bajo de la Alumbrera (8.02 +/- 0.14 and 7.98 +/- 0.14 Ma). Regional uplift and erosion followed during which the magmatic-hydrothermal system was probably unroofed. Shortly thereafter, dacitic lava domes were extruded (7.95 +/- 0.17 Ma) and rhyolitic diatremes (7.79 +/- 0.13 Ma) deposited thick tuff blankets, across the region. Emplacement of large intermediate composition stocks occurred at 7.37 +/- 0.22 Ma, shortly before renewed magmatism occurred at Bajo de la Alumbrera (7.10 +/- 0.07 Ma). The latest porphyry intrusive event is temporally associated with new ore-bearing magmatic-hydrothermal fluids. Other dacitic intrusions are associated with subeconomic deposits that formed synchronously with the mineralized porphyries at Bajo de la Alumbrera. However, their emplacement continued (from 7.10 +/- 0.06 to 6.93 +/- 0.07 Ma) after the final intrusion at Bajo de al Alumbrera. Regional volcanism had ceased by 6.8 Ma (6.92 +/- 0.07 Ma). The brief history of the volcanic complex hosting the Bajo de la Alumbrera Cu-Au deposit differs from that of other Andean provinces hosting porphyry deposits. For example, at the El Salvador porphyry copper district in Chile, magmatism related to Cu mineralization was episodic in regional igneous activity that occurred over tens of millions of years. Bajo de la Alumbrera resulted from the superposition of multiple porphyry-related hydrothermal systems, temporally separated by a million years. It appears that the metal budget in porphyry ore deposits is not simply a function of their longevity and/or the superposition of multiple porphyry systems. Nor is it a function of the duration of the associated cycle of magmatism. Instead, the timing of processes operating in the parental magma body is the controlling factor in the formation of a fertile porphyry-related ore system.

Volcanic resurfacing and the early terrestrial crust: Zircon U-Pb and REE constraints from the Isua Greenstone Belt, southern West Greenland

The oldest known bona fide succession of elastic metasediments Occurs in the Isua Greenstone Belt. SW Greenland and consists of a variety of mica schists and rare metaconglomerates. The metasediments are in direct contact with a felsic metavolcanic lithology that has previously been dated to 3.71 Ga. Based on trace element geochemical data for 30 metasediments, we selected the six samples with highest Zr concentrations for zircon extraction. These samples all yielded very few or no zircon, Those extracted from mica schists yielded ion probe U/Pb ages between 3.70 and 3,71 Ga. One metaconglomerate sample yielded just a single zircon of 3.74 Ga age. The mica schist hosted zircons have U/Pb ages. Th/U ratios, REE patterns and Eu anomalies indistinguishable from zircon in the adjacent 3.71 Ga felsic metavolcanic unit. Trace element modelling requires the bulk of material in the metasediments to be derived from variably weathered mafic lithologies but some metasediments contain substantial contribution from more evolved source lithologies. The paucity of zircon in the mica schists is thus explained by incorporation of material from largely zircon-free volcanic lithologies. The absence of older zircon in the mica schists and the preponderance of mafic source material imply intense, mainly basaltic resurfacing of the early Earth. The implications of this process are discussed, Thermal considerations suggest that horizontal growth of Hadean crust by addition of mafic ultramafic lavas must have triggered self-reorganisation of the protocrust by remelting. Reworking oft Hadean crust may have been aided by burial of hydrated (weathered) metabasalt due to semi-continuous addition of new voluminous basalt Outpouring,;, This process Causes a bias towards eruption of Zr-saturated partial melts at the surface with O-isotope corn posit ion,, potentially different from the mantle. The oldest zircons hosted in sediments would have been buried to substantial depth or formed in plutons that crystallised at some depth from which it took hundreds of millions of years for them to be exhumed and incorporated into much younger sediments. (C) 2005 Elsevier B.V.All rights reserved.

Origin of the Archean Sask craton and its extent within the Trans-Hudson orogen: evidence from Pb and Nd isotopic compositions of basement rocks and post-orogenic intrusions

U-Pb zircon ages from the exposed Sask craton are 2450-3100 Ma, from the Peter Lake Domain 2575-2640 Ma, and from rocks of the Trans-Hudson orogen 1840-1880 Ma. U-Pb monazite and zircon ages of post-orogenic pegmatites and aplites are 1770-1800 Ma. Common Pb and Sm-Nd isotopic compositions of post-orogenic intrusions, as probes of crust beneath the orogen, were compared to Sask craton rocks and ca. 1850 Ma orogenic rocks to infer the origin and subsurface distribution of the Sask craton within the internides of the Trans-Hudson orogen. Results show that post-orogenic intrusions within most of the Glennie Domain and Hanson Lake block were derived, at least in part, from Archean source materials, demonstrating that the Sask craton lies beneath Paleoproterozoic orogenic rocks present at the surface. In contrast, common Pb and Sm-Nd isotopic compositions from pegmatites and aplites of the La Ronge Domain are essentially identical with those of the Paleoproterozoic orogenic rocks into which they are intruded, indicating derivation by partial melting of similar rocks. Thus, if the Sask craton extended to the west beneath the La Ronge Domain, it was beneath the zone of melting that produced the post-orogenic intrusions, making it unlikely that the Sask craton is a detached part of the Hearne craton. Many samples from the Sask craton have elevated Pb-208/Pb-204 ratios, unlike Superior craton or Hearne craton rocks, suggesting that the Sask craton was derived from an exotic source, such as the Wyoming craton, which shares similar elevated Pb-208/Pb-204 ratios.

Applications of accurate, high-precision Pb isotope ratio measurement by multi-collector ICP-MS

The isotope composition of Ph is difficult to determine accurately due to the lack of a stable normalisation ratio. Double and triple-spike addition techniques provide one solution and presently yield the most accurate measurements. A number of recent studies have claimed that improved accuracy and precision could also be achieved by multi-collector ICP-MS (MC-ICP-MS) Pb-isotope analysis using the addition of Tl of known isotope composition to Pb samples. In this paper, we verify whether the known isotope composition of Tl can be used for correction of mass discrimination of Pb with an extensive dataset for the NIST standard SRM 981, comparison of MC-ICP-MS with TIMS data, and comparison with three isochrons from different geological environments. When all our NIST SRM 981 data are normalised with one constant Tl-205/Tl-203 of 2.38869, the following averages and reproducibilities were obtained: Pb-207/Pb-206=0.91461+/-18; Pb-208/Ph-206 = 2.1674+/-7; and (PbPh)-Pb-206-Ph-204 = 16.941+/-6. These two sigma standard deviations of the mean correspond to 149, 330, and 374 ppm, respectively. Accuracies relative to triple-spike values are 149, 157, and 52 ppm, respectively, and thus well within uncertainties. The largest component of the uncertainties stems from the Ph data alone and is not caused by differential mass discrimination behaviour of Ph and Tl. In routine operation, variation of sample introduction memory and production of isobaric molecular interferences in the spectrometer's collision cell currently appear to be the ultimate limitation to better reproducibility. Comparative study of five different datasets from actual samples (bullets, international rock standards, carbonates, metamorphic minerals, and sulphide minerals) demonstrates that in most cases geological scatter of the sample exceeds the achieved analytical reproducibility. We observe good agreement between TIMS and MC-ICP-MS data for international rock standards but find that such comparison does not constitute the ultimate. test for the validity of the MC-ICP-MS technique. Two attempted isochrons resulted in geological scatter (in one case small) in excess of analytical reproducibility. However, in one case (leached Great Dyke sulphides) we obtained a true isochron (MSWD = 0.63) age of 2578.3 +/- 0.9 Ma, which is identical to and more precise than a recently published U-Pb zircon age (2579 3 Ma) for a Great Dyke websterite [Earth Planet. Sci. Lett. 180 (2000) 1-12]. Reproducibility of this age by means of an isochron we regard as a robust test of accuracy over a wide dynamic range. We show that reliable and accurate Pb-isotope data can be obtained by careful operation of second-generation MC-ICP magnetic sector mass spectrometers. (C) 2002 Elsevier Science B.V. All rights reserved.

Inheritance of early Archaean Pb-isotope variability from long-lived Hadean protocrust

Comparison of initial Pb-isotope signatures of several early Archaean (3.65-3.82 Ga) lithologies (orthogneisses and metasediments) and minerals (feldspar and galena) documents the existence of substantial isotopic heterogeneity in the early Archaean, particularly in the Pb-207/Pb-204 ratio. The magnitude of isotopic variability at 3.82-3.65 Ga requires source separation between 4.3 and 4.1 Ga, depending on the extent of U/Pb fractionation possible in the early Earth. The isotopic heterogeneity could reflect the coexistence of enriched and depleted mantle domains or the separation of a terrestrial protocrust with a U-238/Pb-204 (mu) that was ca. 20-30% higher than coeval mantle. We prefer this latter explanation because the high-p signature is most evident in metasediments (that formed at the Earth's surface). This interpretation is strengthened by the fact that no straightforward mantle model can be constructed for these high-mu lithologies without violating bulk silicate Earth constraints. The Pb-isotope evidence for a long-lived protocrust complements similar Hf-isotope data from the Earth's oldest zircons, which also require an origin from an enriched (low Lu/Hf) environment. A model is developed in which greater than or equal to3.8-Ga tonalite and monzodiorite gneiss precursors (for one of which we provide zircon U-Pb data) are not mantle-derived but formed by remelting or differentiation of ancient (ca. 4.3 Ga) basaltic crust which had evolved with a higher U/Pb ratio than coeval mantle in the absence of the subduction process. With the initiation of terrestrial subduction at, we propose, ca. 3.75 Ga, most of the greater than or equal to3.8-Ga basaltic shell (and its differentiation products) was recycled into the mantle, because of the lack of a stabilising mantle lithosphere. We argue that the key event for preservation of all greater than or equal to3.8-Ga terrestrial crust was the intrusion of voluminous granitoids immediately after establishment of global subduction because of complementary creation of a lithospheric keel. Furthermore, we argue that preservation of !3.8-Ga material (in situ rocks and zircons) globally is restricted to cratons with a high U/Pb source character (North Atlantic, Slave, Zimbabwe, Yilgarn, and Wyoming), and that the Pb-isotope systematics of these provinces are ultimately explained by reworking of material that was derived from ca. 4.3 Ga (i.e. Hadean) basaltic crust.

Chemical characterization of earth's most ancient clastic metasediments from the Isua Greenstone Belt, southern West Greenland

New and published major and trace element abundances of elastic metasediments (mainly garnet-biotite-plagioclase schists) from the similar to 3.8 Ga Isua Greenstone Belt (IGB), southern West Greenland, are used in an attempt to identify the compositional characteristics of the protoliths of these sediments. Compositionally, the metasediments are heterogeneous with enrichment of LREE (La/Sm-chord = 1.1-3.9) and variable enrichment and depletion of HREE (Gd/Yb-chord = 0.8-4.3). Chondrite-normalized Eu is also variable, spanning a range from relative Eu depletion to enrichment (Eu/Eu* = 0.6-1.3). A series of geochemical and geological criteria provides conclusive evidence for a sedimentary origin, in disagreement with some previous studies that questioned the presence of genuine elastic metasediments. In particular, trace element systematics of IGB metasediments show strong resemblance to other well-documented Archaean clastic sediments, and are consistent with a provenance consisting of ultramafic, malic and felsic igneous rocks. Two schists, identified as metasomatized mafic igneous rocks from petrographic and field evidence, show distinct compositional differences to the metasediments. Major element systematics document incipient-to-moderate source weathering in the majority of metasediments, while signs of secondary K-addition are rare. Detailed inspection of Eu/Eu*, Fe2O3 and CIW (chemical index of weathering) relationships reveals that elevated iron contents (when compared to averages for continental crust) and strong relative enrichment in Eu may be due to precipitation of marine Fe-oxyhydroxides during deposition or diagenesis on the seafloor. Some of the IGB metasediments have yielded anomalous Nd-142 and W-182 isotopic compositions that were respectively interpreted in terms of early mantle differentiation processes and the presence of a meteorite component. Alternatively, W and possibly Nd isotopes could have been affected by thermal neutron capture on the Hadean surface. The latter process was tested in this study by analysis of Sm isotope compositions, which serve as an effective monitor for neutron capture effects. As no anomalous variation from terrestrial values was detected, we infer that isotope systematics (including W-182 and Nd-142) of IGB metasediments were not affected by neutron capture, but reflect decay of radioactive parent isotopes. Copyright (c) 2005 Elsevier Ltd.

40Ar/39Ar dating of Mn-iron ore of the Sapecado Mine, Qadrilatero Ferrifero, Brazil

The SW region of Amazonian craton presents policyclic evolution between 1.80-1.00 Ga and is comprised of the Rio Negro-Juruena, Rondoniana and Sunsas Provinces. The evolution of this region has being characterized by four orogens: Alto Jauru (1.79-1.74 Ga), Cachoeirinha (1.58-1.52 Ga), Suíte Santa Helena (1.45-1.42 Ga) e Sunsas/Aguapeí (1.0-0.9 Ga). The Alto Jauru orogen consists of TTG gneissic associations, greenstone sequences and intrusive granitoids origined in volcanic arc setting. Eight 40Ar/39Ar step-heating analyses were carried out in minerals (biotiteand hornblende) to investigate the thermal history and crustal evolution of this region. From the Alto Jauru orogen was sampled the gneiss banded and two biotite grains provide large dispersion of apparent ages, suggesting heterogenity in reservoir of the argon. Apparent age diagram yielded integrated ages of 1472 ± 6 Ma, interpreted as minimum ages of regional cooling episode. Three analyses of hornblende present ages varing from 1310 to 1400 Ma, possibly because smaller grain size become more susceptible to argon loss. 40Ar/39Ar step-heating methodology applied on biotite of pyroclastic tuff (U-Pb age about 1758 ± 7 Ma) presented integrated age of 1507 ± 7 Ma. The results found for this terrane demonstrated a geochronological correlation with metamorphic process linked Cachoeirinha orogen. Biotite and hornblende grains separates from granite and a tonalite origined during Cachoeirinha orogen were analyzed and the apparent age diagrams indicated well-defined plateau ages of 1520-1540 Ma. Biotite grains from a granitic sample were analized, and integrated ages about 1526 ± 2 Ma were obtained due argon loss in the initial steps. Thermochronologic history of SW region Amazonian craton is coherent with regional policyclic events and 40Ar/39Ar ages here presented probably correspond to regional cooling period of Cachoeirinha orogen.

Lamproites from Gaussberg, Antarctica: Possible transition zone melts of Archaean subducted sediments

Petrogenetic models for the origin of lamproites are evaluated using new major element, trace element, and Sr, Nd, and Pb isotope data for Holocene lamproites from the Gaussberg volcano in the East Antarctic Shield. Gaussberg lamproites exhibit very unusual Pb isotope compositions (Pb-206/Pb-204 = 17.44-17.55 and Pb-207/Pb-204 = 15.56-15.63), which in common Pb isotope space plot above mantle evolution lines and to the left of the meteorite isochron. Combined with very unradiogenic Nd, such compositions are shown to be inconsistent with an origin by melting of sub-continental lithospheric mantle. Instead, a model is proposed in which late Archaean continent-derived sediment is subducted as K-hollandite and other ultra-high-pressure phases and sequestered in the Transition Zone (or lower mantle) where it is effectively isolated for 2-3 Gyr. The high Pb-207/Pb-204 ratio is thus inherited from ancient continent-derived sediment, and the relatively low Pb-206/Pb-204 ratio is the result of a single stage of U/Pb fractionation by subduction-related U loss during slab dehydration. Sr and Nd isotope ratios, and trace element characteristics (e.g. Nb/Ta ratios) are consistent with sediment subduction and dehydration-related fractionation. Similar models that use variable time of isolation of subducted sediment can be derived for all lamproites. Our interpretation of lamproite sources has important implications for ocean island basalt petrogenesis as well as the preservation of geochemically anomalous reservoirs in the mantle.

On the overabundance of light rare earth elements in terrestrial zircons and its implication for Earth's earliest magmatic differentiation

We present whole-rock and zircon rare earth element (REE) data from two early Archaean gneisses (3.81 Ga and 3.64 Ga) from the Itsaq gneiss complex, south-west Greenland. Both gneisses represent extremely rare examples of unaltered, fresh and relatively undeformed igneous rocks of such antiquity. Cathodoluminescence imaging of their zircons indicates a single crystallisation episode with no evidence for either later metamorphic and/or anatectic reworking or inheritance of earlier grains. Uniform, single-population U/Pb age data confirm the structural simplicity of these zircons. One sample, a 3.64 Ga granodioritic gneiss from the Gothabsfjord, yields a chondrite-normalised REE pattern with a positive slope from La to Lu as well as substantial positive Ce and slight negative Eu anomalies, features generally considered to be typical of igneous zircon. In contrast, the second sample, a 3.81 Ga tonalite from south of the Isua Greenstone Belt, has variable but generally much higher light REE abundances, with similar middle to heavy REE. Calculation of zircon/melt distribution coefficients (D-REE(zircon/melt)) from each sample yields markedly different values for the trivalent REE (i.e. Ce and Eu omitted) and simple application of one set of D-REE(zircon/melt) to model the melt composition for the other sample yields concentrations that are in error by up to two orders of magnitude for the light REE (La-Nd). The observed light REE overabundance in the 3.81 Ga tonalite is a commonly observed feature in terrestrial zircons for which a number of explanations ranging from lattice strain to disequilibrium crystallisation have been proposed and are further investigated herein. Regardless of the cause of light REE overabundance, our study shows that simple application of zircon/melt distribution coefficients is not an unambiguous method for ascertaining original melt composition. In this context, recent studies that use REE data to claim that > 4.3 Ga Hadean detrital zircons originally crystallised from an evolved magma, in turn suggesting the operation of geological processes in the early Earth analogous to those of the present day (e.g. subduction and melting of hydrated oceanic crust), must be regarded with caution. Indeed, comparison of terrestrial Hadean and > 3.9 Ga lunar highland zircons shows remarkable similarities in the light REE, even though subduction processes that have been used to explain the terrestrial zircons have never operated on the Moon. (C) 2002 Elsevier Science B.V. All rights reserved.

Zinc deposits and related mineralization of the Burketown mineral field, including the world-class Century Deposit, Northern Australia: Fluid inclusion and stable isotope evidence for basin fluid sources

The stratiform Century Zn-Pb deposit and the discordant Zn-Pb lode deposits of the Burketown mineral field, northern Australia, host ore and gangue minerals with primary fluid inclusions that have not been affected by the Isan orogeny, thus providing a unique opportunity to investigate the nature of the ore-forming brines. All of the deposits are hosted in shales and siltstones belonging to the Isa superbasin and comprise sphalerite, pyrite, carbonate, quartz, galena, minor chalcopyrite, and minor illite. According to Pb model ages, the main ore stage of mineralization at Century formed at I575 Ma, some 20 m.y. after deposition of the host shale sequence. Microthermometry on undeformed, primary fluid inclusions hosted in porous sphalerite shows that the Zn at Century was transported to the deposit by a homogeneous, Ca2+- and Na+-bearing brine with a salinity of 21.6 wt percent NaCl equiv. delta D-fluid of the fluid inclusion water ranges from -89 to -83 per mil, consistent with a basinal brine that evolved from meteoric water. Fluid inclusion homogenization temperatures range between 74 degrees and 125 degrees C, which are lower than the 120 degrees to 160 degrees C range calculated from vitrinite reflectance and illite crystallinity data from the deposit. This discrepancy indicates that mineralization likely formed at 50 to 85 Mpa, corresponding to a depth of 1,900 to 3,100 m. Transgressive galena-sphalerite veins that cut stratiform mineralization at Century and breccia-filled quartz-dolomite-sphalerite-galena veins in the discordant Zn-Pb lodes have Pb model ages between 1575 and 1485 Ma. Raman spectroscopy and microthermometry reveal that the primary fluid inclusions in these veins contain Ca2+, Na+. but they have lower salinities between 23 and 10 wt percent NaCl equiv and higher delta D-fluid values ranging from -89 to -61 per mil than fluid inclusions in porous sphalerite from Century. Fluid inclusion water from sphalerite in one of the lode deposits has delta O-18(fluid) values of 1.6 and 2.4 per mil, indistinguishable from delta O-18(fluid) values between -0.3 to +7.4 per mil calculated from the isotopic composition of coexisting quartz, dolomite, and illite. The trend toward lower salinities and higher delta D-fluid values relative to the earlier mineralizing fluids is attributed to mixing between the fluid that formed Century and a seawater-derived fluid from a different source. Based on seismic data from the Lawn Hill platform and paragenetic and geochemical results from the Leichhardt River fault trough to the south, diagenetic aquifers in the Underlying Calvert superbasin appear to have been the most likely sources for the fluids that formed Century and the discordant lode deposits. Paragenetically late sphalerite and calcite cut sphalerite, quartz, and dolomite in the lode deposits and contain Na+-dominated fluid inclusions with much lower salinities than their older counterparts. The isotopic composition of calcite also indicates delta O-18(fluid) from 3.3 to 10.7 per mil, which is larger than the range obtained from synmineralization minerals, supporting the idea that a unique fluid source was involved. The absolute timing of this event is unclear, but a plethora of Pb model, K-Ar, and Ar-40/Ar-39 ages between 1440 and 1300 Ma indicate that a significant volume of fluid was mobilized at this time. The deposition of the Roper superbasin from ca. 1492 +/- 4 Ma suggests that these late veins formed from fluids that may have been derived from aquifers in overlying sediments of the Roper superbasin. Clear, buck, and drusy quartz in veins unrelated to any form of Pb-Zn mineralization record the last major fluid event in the Burketown mineral field and form distinct outcrops and ridges in the district. Fluid inclusions in these veins indicate formation from a low-salinity, 300 degrees +/- 80 degrees C fluid. Temperatures approaching 300 degrees C recorded in organic matter adjacent to faults and at sequence boundaries correspond to K-Ar ages spanning 1300 to 1100 Ma, which coincides with regional hydrothermal activity in the northern Lawn Hill platform and the emplacement of the Lakeview Dolerite at the time of assemblage of the Rodinia supercontinent.

Late Quaternary uplift and subsidence of the west coast of Tanna, south Vanuatu, southwest Pacific: U-Th ages of raised coral reefs in the Median Sedimentary Basin

Twelve Late Quaternary TIMS U-Th ages are reported here from 10 coral samples collected in situ from five transgressive coral/algal raised reefs (height: max. 113 m, min. 8 m) and two raised lagoonal deposits (height: max. 18 m, min. 8 m) along and near the west coast of Tanna, which lies in the Median Sedimentary Basin of South Vanuatu, southwest Pacific. These reefs and raised lagoonal deposits represent several age groups: (i) 215 ka (marine oxygen-isotope stage 7) penultimate interglacial (highest elevation and oldest); (ii) one lagoonal deposit of ca 127 ka (marine oxygen-isotope stage 5e); (iii) three last interglacial reefs with ages 102, 89 and 81 ka (representing marine oxygen-isotope stages 5c, 5b and 5a, respectively, of the latter part of the last interglacial); (iv) a lagoonal deposit with a 92 ka age (5b); and (v) a Holocene reef (age >5.7-5.0 ka) (lowest elevation and youngest). A ca 4.9 ka regressive reef (at elevation of 1.5 m above sea-level) is consistent with an island-wide 6.5 m uplift (probably largely coseismic), and a probable further island-wide uplift occurred in the late Holocene. The U-series ages taken together with the heights of transgressive reefs show that uplift since 215 ka was, on average, at similar to0.52 mm/y; although since 5 ka the uplift rate was, on average, similar to1.6 mm/y (the assumption being that a 1.5 m above sea-level reef has a coseismic origin). Elevation of transgressive reefs 5a, 5b and 5c and their ages indicates an island-wide subsidence during the period ?124-89 ka (i.e. Late Quaternary uplift/subsidence was jerky). Late Quaternary uplift/subsidence on the northwest coast of Tanna is considered to be due to irregular thicknesses of crust being subducted beneath Tanna.