Economic Geology; January 2005; v. 100; no. 1;
p. 131-148; DOI: 10.2113/100.1.0131
© 2005 Society of Economic Geologists
Alunite in the Pascua-Lama High-Sulfidation Deposit: Constraints on Alteration and Ore Deposition Using Stable Isotope Geochemistry
C. L. Deyell
Centre for Ore Deposit Research, University of Tasmania, Private Bag 79, Hobart, Tasmania 7001, Australia
R. Leonardson
Barrick Gold Exploration, 293 Spruce Road, Elko, Nevada 89801
R. O. Rye
U.S. Geological Survey, Mail Stop 963, Denver Federal Center, Denver Colorado 80225
J. F. H. Thompson
Teck-Cominco, 200 Burrard Street, Vancouver, British Columbia, Canada V6C 3L9
T. Bissig
Universidad Católica del Norte, Depto. Ciencias Geológicas, Av. Angamos 0610, Antofagasta, Chile
D. R. Cooke
Centre for Ore Deposit Research, University of Tasmania, Private Bag 79, Hobart, Tasmania, 7001, Australia
Corresponding author: e-mail, cdeyell{at}utas.edu.au
 |
Abstract
|
|---|
The Pascua-Lama high-sulfidation system, located in the El Indio-Pascua belt
of Chile and Argentina, contains over 16 million ounces (Moz) Au and 585 Moz Ag.
The deposit is hosted primarily in granite rocks of Triassic age with
mineralization occurring in several discrete Miocene-age phreatomagmatic
breccias and related fracture networks. The largest of these areas is Brecha
Central, which is dominated by a mineralizing assemblage of
alunite-pyrite-enargite and precious metals. Several stages of hydrothermal
alteration related to mineralization are recognized, including all types of
alunite-bearing advanced argillic assemblages (magmatic-hydrothermal,
steam-heated, magmatic steam, and supergene). The occurrence of alunite
throughout the paragenesis of this epithermal system is unusual and detailed
radiometric, mineralogical, and stable isotope studies provide constraints on
the timing and nature of alteration and mineralization of the
alunite-pyrite-enargite assemblage in the deposit.
Early (preore) alteration occurred prior to ca. 9 Ma and consists of intense
silicic and advanced argillic assemblages with peripheral argillic and
widespread propylitic zones. Alunite of this stage occurs as fine intergrowths
of alunite-quartz ± kaolinite, dickite, and pyrophyllite that selectively
replaced feldspars in the host rock. Stable isotope systematics suggest a
magmatic-hydrothermal origin with a dominantly magmatic fluid source. Alunite is
coeval with the main stage of Au-Ag-Cu mineralization (alunite-pyrite-enargite
assemblage ore), which has been dated at approximately 8.8 Ma. Ore-stage alunite
has an isotopic signature similar to preore alunite, and
34Salun-py
data indicate depositional temperatures of 245° to 305°C. The
D
and
18O data exclude significant involvement of
meteoric water during mineralization and indicate that the assemblage formed
from H2S-dominated magmatic fluids. Thick steam-heated alteration
zones are preserved at the highest elevations in the deposit and probably formed
from oxidation of H2S during boiling of the magmatic ore fluids.
Coarsely crystalline magmatic steam alunite (8.4 Ma) is restricted to the
near-surface portion of Brecha Central. Postmineral alunite ± jarosite were
previously interpreted to be supergene crosscutting veins and overgrowths,
although stable isotope data suggest a mixed magmatic-meteoric origin for this
late-stage alteration. Only late jarosite veinlets (8.0 Ma) associated with
fine-grained pseudocubic alunite have a supergene isotopic signature.
The predominance of magmatic fluids recorded throughout the paragenesis of
the Pascua system is atypical for high-sulfidation deposits, which typically
involve significant meteoric water in near-surface and peripheral alteration
and, in some systems, even ore deposition. At Pascua, the strong magmatic
signature of both alteration and main-stage (alunite-pyrite-enargite assemblage)
ore is attributed to limited availability of meteoric fluids. This is in
agreement with published data for the El Indio-Pascua belt, indicating an event
of uplift and subsequent pediment incision, as well as a transition from
semiarid to arid climatic conditions, during the formation of the deposit in the
mid to late Miocene.
 |
Introduction
|
|---|
THE PASCUA-LAMA Au-Ag-Cu deposit is situated in the Cordillera Principal of
Chile and Argentina at the north end of the El Indio-Pascua belt (Fig. 1).
Widespread zones of hydrothermal alteration occur throughout the district and
are visible as distinct color anomalies and topographic features at elevations
of 4,500 to 5,100 m above sea level (asl). Recognition of these zones led to the
discovery of high-sulfidation gold mineralization at Nevada (now Esperanza) in
the Chilean region in 1977 by Companía Minera San José S.A. The property was
explored by a number of different operators, including Lac Minerals Ltd., until
the acquisition of this company by Barrick Gold Corporation in 1994. Continued
exploration by Barrick Gold has uncovered several zones of both sulfide-rich and
oxidized mineralization to the east of the original discovery and across the
international border into Argentina. Estimated reserves for the deposit total
16.8 million ounces (Moz) Au and 585 Moz Ag.

View larger version (38K):
[in this window]
[in a new window]
|
FIG. 2. Simplified geology and major faults of the El Indio-Pascua belt,
showing locations of major ore deposits in the region. Geologic information is
taken from Martin et al. (1995) and Ramos et al. (1989), as summarized in Bissig
et al. (2002). The upper Miocene Pascua Formation is not shown due to its
restricted occurrence (see text). BdTF = Baños del Toro fault.
|
|
Pascua-Lama is classified as a high-sulfidation system because of the
dominant enargite-pyritehosted Au-Ag-Cu mineralization and associated
advanced argillic alteration. Several distinct types of advanced argillic
alteration and alunite can be recognized, including near-surface steam-heated
zones, widespread magmatic-hydrothermal alunite-kaolinite (±
dickite-pyrophyllite) assemblages, coarse magmatic steam alunite as defined by
Rye et al. (1992), and late, apparently supergene, alunite-jarosite veinlets and
disseminations. Alunite also commonly formed contemporaneously with, and is host
to, precious metal mineralization. This style of mineralization occurs as an
assemblage of alunite-pyrite-enargite and contains most of the identified gold
resource in the region.
The coprecipitation of alunite with precious metal-bearing sulfides is an
unusual feature of the Pascua-Lama system. In most high-sulfidation deposits,
advanced argillic alteration and alunite deposition clearly precede
mineralization, regardless of how the ore was emplaced (e.g., Stoffregen, 1987;
Rye 1993; Cooke and Simmons, 2000, and references therein). The occurrence of
significant alunite as gangue to mineralization is seldom documented, with the
exception of epithermal veins at la Mejicana, Nevados del Famatina district,
Argentina (Losada-Calderon and McPhail, 1996), late alunite-enargite veins at El
Indio (Jannas et al., 1999), and breccia- and vein-hosted mineralization at
Tambo (Jannas et al., 1999; Deyell et al., in press). In this study, we describe
the occurrence of alunite-sulfide-gold mineralization at Pascua-Lama and
associated advanced argillic alteration assemblages. Stable isotope data and 40Ar-39Ar
alunite ages are used to constrain the origin of both alteration and ore fluids
and detail the multistage evolution of this richly mineralized epithermal
system.
 |
Geology of the Pascua-Lama Region
|
|---|
Regional framework
The Pascua-Lama Au-Cu deposit is located at the northern end of the El
Indio-Pascua belt in the Central Andean Cordillera Principal, straddling the
Chile-Argentina border (Fig. 1). The belt occurs in the middle of the flat-slab
segment of the Andes between two other well-known mineralized areas: the
Maricunga porphyry Au and high-sulfidation Au-Ag belt to the north (26°28°
S) and the Los Pelambres-El Teniente porphyry Cu-Mo belt (32°34° S) to the
south (e.g., Kay et al., 1999). In addition to Pascua-Lama, the El Indio-Pascua
belt hosts widespread zones of hydrothermal alteration and several world-class
epithermal deposits and prospects (Fig. 2), including El Indio, which has
recently closed, and Veladero, with estimated resources of 9.4 Moz Au.
Mineralized systems in the El Indio-Pascua belt formed in the late Miocene with
all significant high-sulfidation mineralization occurring between 6.2 and 9.4
Ma, although only minor associated magmatism is recognized (Bissig et al.,
2001).
The El Indio-Pascua belt occurs within a north-southstriking tectonic
depression bound by the steeply dipping, east-verging, high-angle reverse Baños
del Toro-Chollay fault to the west (Fig. 2) and opposing west-verging structures
(informally named the El Indio high-angle reverse fault zone) farther to the
east. Within this block, an Upper Paleozoic to Lower Jurassic basement composed
predominantly of calc-alkaline felsic intrusions and volcanic rocks (Martin et
al., 1999) is overlain by up to 1,500 m of Tertiary subaerial volcanic strata
(Fig. 2). The latter are extensively preserved in the southern part of the belt,
where the Tambo and El Indio mines are located, but only limited exposures are
preserved in the northern extremity of the district, near Pascua-Lama and
Veladero (Maksaev et al., 1984; Martin et al., 1995). The Tertiary volcanic
sequence comprises predominantly volcanic and volcaniclastic rocks of the Tilito
Formation (2723 Ma), the Escabroso Group (2117.5 Ma), and the 17 to 14 Ma
Cerro de las Tortolas Formation (Martin et al., 1995; Bissig et al., 2001).
Magmatism decreased markedly after the deposition of these units. Isolated
dacitic tuffs of the Vacas Heladas Formation erupted between 12.7 and 11 Ma. A
single 7.84 Ma dacite dike, the only igneous rock coeval with the mineralization
in the district, is reported from Pascua (Bissig et al., 2001). The Vallecito
Formation (6.25.5 Ma) is only exposed in the Valle de Cura region (Ramos et
al., 1989) and the El Indio-Tambo area. Volcanism ceased in the upper Pliocene
after the eruption of the 2 Ma Cerro de Vidrio rhyolite dome in the northeast
part of the belt (Bissig et al., 2002a).
Local geology
The geology of the immediate Pascua area (Fig. 3) is dominated by large
volumes of Mesozoic granitic rocks that include fine- to coarse-grained
porphyritic granite and small stocks of dacite porphyry and granodiorite. These
granitic units have been dated at approximately 190 Ma (Bissig et al., 2001) and
may be related to the regionally extensive and slightly older (212 Ma) Colorado
Group intrusive rocks identified by Martin et al. (1999). Crystal and lithic ash
flow tuffs of probable Lower Permian age (265 Ma; Martin et al., 1995) are the
oldest exposed rocks in the region and occur on the western edge of the
Esperanza sector and the southern and southeastern margins of the property.
Small leucocratic diorite stocks and a larger dacite porphyry of Tertiary age,
as well as various small younger dikes of andesitic to rhyolitic composition
(not shown), have intruded the older intrusive and extrusive units.

View larger version (39K):
[in this window]
[in a new window]
|
FIG. 3. Generalized geology of the Pascua-Lama deposit (based on regional
mapping by Barrick geologists and data presented in Bissig et al. (2001) and
Chouinard (2003); see Figure 2 for location. Mesozoic granitic rocks (the Pascua-Lama
Complex) include several bodies of similar composition, ranging from porphyritic
granite to crowded granite porphyry. Several generations of Tertiary
hydrothermal breccias have been grouped together and include Brecha Central, the
largest breccia body in the region. Also shown are the trace of the Alex Tunnel
(4,680 m asl) and line of section CA-00 (section A-A') illustrated in Figure 5.
Locations of 40Ar/39Ar samples described in the text and
Table 2 are given (projected to surface; abbreviations as given in Table 1).
|
|

View larger version (57K):
[in this window]
[in a new window]
|
FIG. 5. Pascua-Lama property map showing the simplified distribution of
alteration zones at surface (using unpub. data from D. Heberlein, 1999).
|
|
View this table:
[in this window]
[in a new window]
|
TABLE 1 TABLE 1. Paragenetic Sequence and Characteristics of the Different Stages of Advanced Argillic Alteration and Alunite Deposition at Pascua
|
|
Several generations of breccias are recognized in the Pascua-Lama area,
including large, premineral phreatomagmatic breccia pipes and small postmineral
pebble dikes and alunite-jarosite matrix breccias (Chouinard, 2003). The most
significant breccia in the region is Brecha Central, which hosts roughly 25
percent of total gold mineralization. This breccia is sited on and elongate
parallel to east-southeast fractures of the Quebrada de Pascua regional
lineament and contains both matrix- and clast-supported components. The breccia
formed prior to the main stage of mineralization and is inferred to be a relict
maar-diatreme complex, with tuff rings that slumped back into the vents and
produced blocks of subhorizontal volcaniclastic sediments within the breccia-filled
vent (Chouinard, 2003). Mineralization followed this main stage of brecciation,
with significant Au-Ag-Cu ore deposition occurring in several discrete centers.
 |
Mineralization
|
|---|
Several styles of Au (+ Ag-Cu) mineralization are recognized at Pascua. The
greatest volume of gold mineralization is contained within an assemblage of
alunite, pyrite, and enargite that is concentrated in the Brecha Central area.
Other styles of mineralization, described in detail by Chouinard (2003), include
either pyrite, pyrite-szomolnokite (FeSO4 H2O), or
native gold-dominated assemblages, which typically occur above and lateral to
alunite-pyrite-enargite ore. High-grade Ag-rich mineralization also occurs in
the upper parts of the deposit (>4,700 m) but formed late and overprinted the
main-stage Au-Ag-Cu ore (Chouinard, 2003). Although these peripheral mineralized
zones contain significant metal accumulations, this paper focuses exclusively on
alunite-pyrite-enargitetype mineralization, where the precious metal
deposition is directly associated with coeval alunite.
Alunite-pyrite-enargite mineralization
Alunite-pyrite-enargite mineralization occurs as disseminations and
open-space filling in the Brecha Central pipe and as banded alunite-sulfide
veins (Fig. 4A). Textural evidence from mineral intergrowths and inclusions
indicates that alunite was deposited synchronously with the sulfides and
precious metals. Gold and silver are hosted primarily in enargite and pyrite,
predominantly bound within the sulfide crystal lattices (Chouinard, 2003) but
also as inclusions of native gold, calaverite (AuTe2), and minor
electrum. Alunite is intimately intergrown with, and locally host to, precious
metal-bearing sulfides (Fig. 4B) and locally contains native gold. Accessory
phases in the alunite-pyrite-enargite assemblage include quartz, diaspore,
anglesite, pyrophyllite, stibnite, cassiterite, goldfieldite, covellite, galena,
and trace chalcopyrite. Native sulfur was also deposited locally with the
alunite-pyrite-enargite assemblage although most is paragenetically late and was
precipitated near the end of the mineralizing event (Chouinard, 2003).

View larger version (115K):
[in this window]
[in a new window]
|
FIG. 4. Examples of alunite-pyrite-enargite mineralization. A. Banded veins
in the stockwork zone surrounding Brecha Central. Alunite (white bands)
alternates with pyrite-enargite (dark bands). B. Backscattered electron image
showing alunite intergrown with pyrite and enargite. The latter hosts inclusions
of calaverite (AuTe2). DDH-111, 200.5m.
|
|
 |
Alteration
|
|---|
Widespread zones of hydrothermal alteration are recognized in the Pascua-Lama
region. Alteration is typically pervasive and texturally destructive, and the
relationships between alteration assemblages are complex. Subtle changes in
mineralogy, particularly where clay and sericitic minerals are present, are not
easily identified in the field. A detailed alteration study using samples
systematically selected from surface and underground exposures, drill core, and
reverse-circulation drill holes has clarified the distribution of alteration
assemblages (Deyell, 2001). The character and mineralogy of the assemblages were
determined by short-wave infrared spectroscopy (SWIR) techniques with
supplemental petrography and X-ray diffraction analyses.
At a broad scale, intense silicic and advanced argillic alteration zones are
associated with major breccias and mineralized centers (Fig. 5). Silicic
alteration is dominated by patchy vuggy silica zones, generally above about
4,750 m asl, consisting of fine-grained secondary quartz, rare primary quartz
phenocrysts, and relict rutile grains (Chouinard, 2003). Local silicification is
also evident as fine-grained granular quartz that occurs as veinlets, infilling
voids and vugs, or as a pervasive replacement of advanced argillic-altered
rocks. These zones of intense alteration grade outward to peripheral argillic
and propylitic assemblages. Thick blankets of steam-heated advanced argillic
alteration (Table 1) are also preserved at upper elevations, particularly above
the Brecha Central zone (Fig. 6).

View larger version (49K):
[in this window]
[in a new window]
|
FIG. 6. Section CA-00 (see Fig. 3 for location). A. Distribution of major
lithological units. B. Distribution of alteration assemblages and Au-grade
contours. C. Distribution of phyllosilicates (kaolinite, dickite, pyrophyllite)
within the advanced argillic alteration assemblage. Mineral abbreviations: alun
= alunite, dick = dickite, ill = illite, jar = jarosite, kao = kaolinite,
pyrophyllite (pyl), quartz (qtz).
|
|
In detail, several stages of advanced argillic alteration and all of the
major types of alunite deposition (cf. Rye et al., 1992) are recognized in the
Pascua region, in addition to the alunite-pyrite-enargite stage of
mineralization. These assemblages and their respective ages (Table 1) are
described below, followed by the stable isotope systematics, which are used to
constrain the origin of alunite and alunite-bearing mineral assemblages. All age
determinations were carried out at the Queens University (Kingston, Canada) 40Ar-39Ar
laboratory, equipped with a Mass Analyzer Products MAP 216 mass spectrometer.
About 10 mg of each sample was irradiated for 7.5 h at the McMaster nuclear
reactor in Hamilton, Canada, using biotite standard Mac-83 as a radiation-flux
monitor (24.36 ± 0.17 Ma; Sandeman et al., 1999). Samples were step heated
using a defocused LEXEL 3500 Ar laser beam. Ages were calculated using the decay
constants suggested by Steiger and Jäger (1977), and all errors are given at 2
.
Analytical details are provided in Table 2.
Preore advanced argillic alteration
The earliest recognized advanced argillic alteration in the Pascua area is
characterized by fine-grained intergrowths of alunite ± quartz, kaolinite
and/or dickite and/or pyrophyllite, and pyrite. This assemblage is widespread
and extends from surface to the maximum depth investigated by drilling (between
5,200- and 4,100-m elev). Alteration formed prior to brecciation at about 9.1 Ma
(Bissig et al., 2001) and is superimposed by the main mineralizing event,
although no reliable 40Ar/39Ar ages have been obtained for this alteration.
Alunite in the preore advanced argillic assemblage occurs as aggregates that
replaced feldspar phenocrysts, as disseminations in the wall-rock matrix, and
rarely as irregular alunite veinlets and stringers. Only primary quartz
phenocrysts and minor rutile remain from the original wall rock. Alunite is the
main constituent of the assemblage and comprises 10 to 40 vol percent of the
altered rocks. Preore alunite contains rare inclusions of
aluminum-phosphate-sulfate minerals (primarily walthierite, svanbergite,
woodhouseite, and florencite), pyrite, and diaspore (Fig. 7A). Zunyite also
occurs in association with alunite but exclusively at depths below a 4,500-m
elevation and only where pyrophyllite is present.

View larger version (161K):
[in this window]
[in a new window]
|
FIG. 7. Photographs of alteration minerals. A. Backscattered electron image
showing irregular alumino-phospho-sulfate grains (bright zones) in cores of
alunite associated with advanced argillic alteration (alun). DDH-116, 289m. B.
Steam-heated alteration in drill core sample DDH 119-47m above the Frontera
zone. Alunite + quartz replaces feldspar phenocrysts. C. Backscattered electron
image of ore-stage alunite containing REE-rich alumino-phospho-sulfate (APS)
inclusions (florencite). Sampled from DDH-111, 189.9m. D. Backscattered electron
image showing oscillatory PO4 ± Sr-enriched bands (light bands) in
coarse-grained magmatic steam alunite (sample PS-26c). E. Backscattered electron
image of pseudocubic alunite (gray) with overgrowing jarosite (white) in
late-stage vein. F. Backscattered electron image showing supergene alunite,
jarosite, and intermediate alunite-jarosite solid solution. Sample PM-33, Alex
Tunnel.
|
|
A vertical zonation in silicate mineralogy associated with this assemblage is
evident in the Brecha Central area (Fig. 6C). Pyrophyllite is dominant below
about a 4,450-m elevation, whereas dickite typically occurs between a 4,450- to
4,850-m elevation. Kaolinite is restricted to areas peripheral to Brecha Central
and to the steam-heated alteration zone (see below).
Steam-heated alteration
Blanketlike zones of silica minerals (cristobalite, opal, chalcedony) ±
kaolinite, alunite, and native sulfur occur at or near the present-day surface
at elevations greater than approximately 4,900 m (Fig. 6B). This alteration
extends to greater depths along structures, locally crosscutting preore advanced
argillic and peripheral argillic assemblages. The tabular, near-surface
distribution of this assemblage, in addition to the alteration mineralogy, is
consistent with the steam-heated environment (e.g., Schoen et al., 1974). Alteration is typically pervasive, occurring as powdery, friable masses, and
original textures are only locally preserved. Alunite is a minor constituent in
the surficial blanket assemblage and typically occurs toward the base of the
alteration zones. Where present, alunite occurs as extremely fine grained
(<5-µm), pseudocubic, tabular, and anhedral crystals in
irregular patches and disseminations in veinlets and open spaces (Fig. 7B).
Fine-grained jarosite occurs locally as overgrowths on alunite.
The timing of steam-heated alteration is unresolved. Attempts to date alunite
from this zone were mostly unsuccessful and only an approximate 40Ar/39Ar
age of 9.14 ± 1.98 Ma was obtained (Table 2). Textural relationships suggest
that steam-heated assemblages formed after preore advanced argillic wall-rock
alteration, but these zones are never observed in contact with
alunite-pyrite-enargitestage mineralization. Based on the overlap in age data
with this event (Fig. 8), steam-heated alteration may have been, at least in
part, coeval with mineralization.

View larger version (23K):
[in this window]
[in a new window]
|
FIG. 8. Range of 40Ar/39Ar age data (with 2
error bars) for alunite from the alunite-pyrite-enargite (APE)-stage ore,
steam-heated alteration (SH), magmatic steam (MS), late-stage veins (LV), and
one sample of supergene jarosite (SP). Also indicated is the approximate minimum
age for preore advanced argillic alteration (AA1), based on age dates for
advanced argillic alteration in the Frontera and Lama areas from Bissig et al.
(2001). The age of the Pascua Formation rhyodacite dike is also reported in
Bissig et al. (2001). The cross-hatched area represents the possible age range
of APE mineralization, given constraints from preore AA1 alteration and postore
MS alunite.
|
|
Alunite-pyrite-enargite ore
Alunite associated with the alunite-pyrite-enargite assemblage occurs as
coarse-grained clusters and rosettes in open vugs of the breccia matrix and as
densely intergrown euhedral white, gray, or pink crystals in the banded veins.
Grains are typically tabular to bladed, ranging from 30 µm to
>2 cm in length. They contain inclusions of pyrite, enargite, galena, native
sulfur, barite, anglesite, aluminium-phosphate-sulfate minerals (walthierite and
florencite; Fig. 7C), rare diaspore, and pyrophyllite. Compositional zoning
between alunite (K-rich) and walthierite (Ba-rich) end members is common.
Minor alunite + quartz also occur in restricted alteration halos peripheral
to veins hosting alunite-pyrite-enargite mineralization. This is thought to be a
wall-rock alteration envelope related to the alunite-pyrite-enargite veins,
although the extent of these envelopes is difficult to determine because they
grade into the more pervasive preore advanced argillic assemblages. For this
reason, we restricted our sampling and analysis to alunite clearly related (as
gangue) to alunite-pyrite-enargite ore.
The 40Ar/39Ar ages for alunite range from 8.78 ± 0.63
(this study, Table 2) to 8.1 ± 0.1 Ma (Bissig et al., 2001). We consider the
8.78 Ma age to be a better representation of the alunite-pyrite-enargite
assemblage because the younger samples reported in Bissig et al. (2001) were
taken from a smaller oxidized breccia center that is now interpreted to have
formed after the peak of the mineralizing event. However, due to the relatively
large analytical error on the reported alunite sample, the age of mineralization
can only be constrained to between about 9.1 and 8.2 Ma (Fig. 8).
Magmatic steam alunite
Another style of alunite deposition is recognized in isolated surface
exposures of coarse vug- and matrix-filling alunite above Brecha Central. This
alunite is medium to coarse grained and clear to golden yellow in color.
Electron microprobe data indicate that this alunite is nearly stoichiometric
[KAl3(SO4)2(OH)6] with minor Na and
P substitution in oscillatory growth zones (Fig. 7D). Rare inclusions of barite
and jarosite are present. The coarse-grained nature of this alunite and lack of
coeval sulfides (in addition to stable isotope characteristics described below)
are consistent with a magmatic steam origin. The magmatic steam environment was
initially defined by Rye et al. (1992) on the basis of data for veins of
coarsely crystalline alunite at Marysvale, Utah (Cunningham et al., 1984).
The reported age of 8.38 ± 0.17 Ma (Table 2) indicates that magmatic steam
alunite may be younger than the alunite-pyrite-enargite mineralizing event,
although the two stages overlap within analytical error.
Late-stage alunite veins
A late stage of alunite occurs as veins, veinlets, and local disseminations
throughout the Pascua district. Alunite is typically fine to medium grained (520
µm) and occurs as tabular to anhedral crystals locally intermixed
with silica. A pseudocubic crystal habit is observed locally in finer grained
(<5-µm) alunite (Fig. 7E). Jarosite is also common in the
veins, but in all examples it has overprinted alunite. This assemblage occurs as
either a fine-grained, granular to powdery pervasive alteration,
cryptocrystalline alunite (±jarosite, kaolinite) veins, or as medium- to
coarse-grained, euhedral crystals in veins and disseminations. Crosscutting
relationships indicate that this stage of alteration postdates both preore
advanced argillic alteration and alunite-pyrite-enargite mineralization.
Repeated attempts to obtain a reliable 40Ar/39Ar age for
late-stage alunite were unsuccessful due to insufficient release of argon,
although a rough age of 7.97 ± 1.59 Ma is reported (Table 2).
Supergene jarosite (± alunite)
Jarosite is present as a late mineral in association with advanced argillic,
silicic, and argillic (quartz-illite) alteration assemblages. It occurs
interstitial to alunite, quartz, and illite as fine-grained disseminations,
granular crystals in vugs and open spaces, and locally as veinlets that dissect
alunite grains or crosscut preore advanced argillic and argillic alteration. The
distribution of jarosite is strongly fault and fracture controlled and extends
to the maximum depths investigated by drilling (4,100-m elevation). Accessory
hematite is present locally. Alunite rarely occurs intergrown with jarosite, in
thin veinlets and fine-grained aggregates in the rock matrix. The two cannot be
distinguished in thin section, but SEM-EDS analysis shows pseudocubic alunite
grains (<5 µm) overgrown by successive bands of jarosite and
compositional zones intermediate to alunite-jarosite end members (Fig. 7F). An
age of 7.98 ± 0.43 Ma (this study, Table 2) for jarosite of this type is the
youngest alteration event recorded in the Pascua-Lama district.
 |
Stable Isotope Study
|
|---|
Stable isotope data for alunite and associated alteration minerals are used
to characterize the type and source of alteration and alunite-pyrite-enargite
mineralizing fluids. Alunite from each of the assemblages described above and
two samples of jarosite was analyzed for
34S and a
subset of samples was selected for
D and
18O
analysis. Sulfur data were also obtained for ore-stage pyrite, enargite, and
barite. All analyses were conducted at the U.S. Geological Survey isotope
laboratory in Denver. The
34S analyses were carried
out by an online method, using an elemental analyzer coupled to a Micromass
Optima mass spectrometer, following the method of Giesemann et al. (1994).
Hydrogen isotopes were determined by a step-heated and cryogenic technique
modified from Godfrey (1962), as described in Wasserman et al. (1992). Oxygen
isotope data were collected for alunite (both sulfate and hydroxyl groups) and
barite, using the conventional BrF5 method described in Wasserman et
al. (1992). Several samples were also analyzed by pyrolysis with a Finnigan TC/EA
coupled to a Finnigan Delta Plus XL mass spectrometer, using continuous-flow
methods modified from Kornexl et al. (1999). Details of sample preparation and
analytical methods are reported in Deyell (2001). Brief sample descriptions and
all isotope data are listed in Table 3.
Sulfur isotopes
All sulfur isotope analyses are summarized in Figure 9. The
34S
value of bulk sulfur in the magma (
34S
S)
cannot be determined directly since no large bodies of synore volcanic rocks
have been identified in the Pascua district, with the exception of the 7.84 Ma
rhyodacite dike noted previously. Estimation of
34S
S
(2.8 ± 2
) is therefore based on data for magmatic steam alunite (see below).
Preore advanced argillic alunite has
34S values
that range from 13 to 25 per mil (Fig. 9). Pyrite could not be separated in
sufficient quantity for analysis from preore advanced argillic alteration, but
R. Beane (unpub. report for Barrick Chile, Ltda., 1988) reported a
34S
value of 3.5 per mil for disseminated pyrite intergrown with alunite in
quartz monzonite from the Brecha Central area. This result is plotted for
comparison in Figure 9. The range of
34S values for
alunite associated with the alunite-pyrite-enargite assemblage is identical to
that for alunite in preore advanced argillic alteration, and data for ore-stage
barite (22.1
) also fall within this range. Ore-stage enargite and pyrite have
34S values that cluster around cluster around 3 to 6 per
mil, consistent with the value for pyrite from preore advanced argillic
alteration.
Alunite
34S data for the preore advanced argillic- and
alunite-pyrite-enargite assemblages are significantly higher than all other
stages of alunite and jarosite in the deposit (Fig. 9). The 14 to 30 per mil
difference between these alunites and associated sulfides is consistent with
sulfate derived from the disproportionation of SO2 (Rye et al., 1992;
Rye, 1993) and is interpreted to reflect a magmatic-hydrothermal origin. The
large variation in
34S values for alunite is
attributed to deposition over a range of temperatures. These temperatures can be
estimated using the isotopic fractionation equation of Ohmoto and Lasaga (1982)
for coexisting alunite and pyrite. Data for the one alunite-pyrite pair from a
preore advanced argillic sample reported by R. Beane (unpub. report, 1988) give
an equilibrium
34Salun-py temperature of
380°C. Using an average
34S value for pyrite of -4
per mil, depositional temperatures for other preore advanced argillic alunite
samples are estimated to be 190° to 350°C. Likewise, equilibrium temperatures
calculated for three alunite-pyrite-enargite assemblage alunite-pyrite pairs
(Table 4) range from 245° to 305° C. The H2S/SO4
ratios for these same pairs average about 2.
Steam-heated alunite has
34S values of 0 to 6
per mil, which is slightly higher than
34S values
of sulfides in the deposit (Fig. 9). These results are consistent with a
steam-heated origin for this stage of alteration (e.g., Rye et al., 1992),
although the range of
34S data indicates that
limited sulfur isotope exchange occurred between aqueous SO4 and H2S
during this alteration event.
Magmatic steam alunite was sampled from coarse-grained, euhedral crystals
obtained near surface in the matrix of Brecha Central. This alunite has a
34S
value of 2.8 per mil, which is slightly higher than sulfides in the underlying
deposit but significantly lower than alunite from either the preore advanced
argillic or alunite-pyrite-enargite assemblages (Fig. 9).
The
34S values for late-stage alunite range from
4 to 7 per mil. These values are higher than the
34S
values of sulfides in the deposit and suggest that the sulfate was derived from
either the oxidation of sulfides or H2S degassed during the collapse
of the hydrothermal system. In either case, limited sulfur isotope exchange
between sulfate and sulfide species is inferred based on the positive
34S
values for this stage of alteration.
Only two samples of jarosite were separated for analysis. One sample, taken
from a series of fine-grained veinlets approximately 300 m below the present-day
surface, has the lowest
34S value (2.5
) of
all sulfates in the deposit. This value approaches the
34S
values for pyrite associated with the alunite-pyrite-enargite assemblage and is
consistent with sulfate derived from the oxidation of precursor sulfides. A
second sample of jarosite, from veinlets that crosscut preore advanced argillic-altered
wall rock, has a
34S value (3.3
), similar to
that of late-stage alunite.
Oxygen and hydrogen isotopes
Oxygen and hydrogen isotope compositions for each stage of alteration are
given in Figures 10 through 12. Fluid compositions in equilibrium with alunite
of each stage are calculated from the fractionation equations of Stoffregen et
al. (1994) over a range of temperatures, as specified. These temperatures are
estimated from
34Salun-py values or
18OSO4-OHisotope fractionations for each stage, where available. The
D
of local meteoric water is estimated to have been 110 ± 10 per mil based on
paleoelevation models (B. Taylor, pers. commun.) and previous studies of the El
Indio-Pascua belt (Jannas et al., 1999; Deyell, 2001).

View larger version (20K):
[in this window]
[in a new window]
|
FIG. 10. D, 18OSO4,
and 18OOH values for alunite in the
preore advanced argillic assemblage (AA1). Fluid compositions ( DH2O
and 18OH2O) in equilibrium with alunite
are calculated from equations of Stoffregen et al. (1994) at 200° to 380°C
(based on range of 34Salun-py for AA1
alteration). Reference lines and fields shown include: meteoric water line
(Craig, 1961), kaolinite line (Savin and Epstein, 1970), typical fluids
dissolved in felsic magmas (Taylor, 1988), the range of water compositions
discharged from high-temperature fumaroles (volcanic vapor; Giggenbach, 1992),
and the composition of paleometeoric waters in the El Indio-Pascua belt
(estimated from paleotopography; B. Taylor, pers. commun., 2001).
|
|
Preore advanced argillic alteration: Alunite from this stage of
alteration is characterized by high
18OSO4
values of 14 to 21 per mil (Fig. 10). The composition of preore alunite fluids
overlaps with that of typical felsic magmatic fluids (Taylor, 1988) and
indicates a dominant magmatic component in the source fluids (cf. Rye et al,
1992). The trend of the data in Figure 10, which shows gradual decreasing
Dalun
nearly parallel to the meteoric water line, suggests some mixing of these
magmatic fluids with exchanged meteoric water, resulting from either water-rock
interaction and/or boiling or evaporation. Depositional temperatures calculated
from
18OSO4-OH range from 50° to
+140°C and are attributed to retrograde isotopic exchange in the OH site.
Alunite-pyrite-enargite ore: Calculated fluid compositions for ore-stage
alunite are nearly identical to preore alunite, with
DH2O
ranging between 30 and 47 per mil (Fig. 11). Again, these data overlap
with the
DH2O of magmatic water (Taylor, 1988) and
indicate a dominantly magmatic source for the mineralizing fluid. The
18OSO4-OH
temperatures (50° to +790°C) indicate disequilibrium, as noted above for
preore advanced argillic alunite.
Steam-heated alunite: Fluid compositions in equilibrium with steam-heated
alunite (Fig. 11) were calculated over a temperature range of 60° to 140°C.
These temperatures are based on
18OSO4-OH
data for three steam-heated alunite samples, which correspond to 60°, 140°,
and 270°C. The lower temperatures are reasonable for the steam-heated
environment, in which aqueous sulfate forms at shallow levels from the oxidation
of H2S by atmospheric oxygen (Rye et al., 1992; Ebert and Rye, 1997),
and suggest that oxygen isotope equilibrium was obtained between aqueous sulfate
and water in this fluid. The third sample gives an unreasonably high
18OSO4-OH
temperature, given the shallow depth of formation, although
18OOH
data are in the range of the previous two samples. These data indicate that
significant sulfur isotope exchange occurred between aqueous SO4 and
H2S, but
18OSO4-OH may have
been affected by postdepositional exchange for this sample. Calculated fluid
compositions for steam-heated alunite, as shown in Figure 11, have
DH2O
values of 38 to 54 per mil and
18OH2O
values between 5 and +4 per mil. These values overlap with those calculated
for ore-stage alunite and are significantly higher than those for meteoric
water.
Magmatic steam alunite: Depositional temperatures for alunite of this
stage are poorly constrained. A calculated
18OSO4-OH
temperature of 110°C is lower than expected for magmatic steam alunite, given
the coarse-grained nature of the sample and its inferred magmatic source, and
may indicate equilibration with a low
DH2O fluid.
Similar processes have been inferred for magmatic steam alunite from Marysvale,
Utah (Rye et al., 1992), and Tambo, Chile (Deyell et al., in press). Calculation
of fluid compositions in equilibrium with this alunite was therefore made over a
range of assumed temperatures from 150° to 250°C. These fluid compositions
(Fig. 12) overlap with those for preore advanced argillic alteration and
alunite-pyrite-enargite mineralization and indicate a major magmatic component
in the source fluid. Combined with a low
34S value,
lack of coeval sulfides, and the coarse-grained nature of the alunite, these
data are consistent with a magmatic steam origin (cf. Rye et al., 1992). Alunite
of this type is considered to have
34S values close
to that of the bulk sulfur in the magma, based on relationships determined by
Cunningham et al. (1984) and Rye et al. (1992) for alunite at Marysvale, Utah.
Late-stage alunite veins: Data for late-stage alunite are plotted in
Figure 12, which includes reference fields for supergene alunite and jarosite as
given in Rye et al. (1992) and Rye and Alpers (1997), respectively. Late-stage
alunite
18OOH values range from 5.5 to
13.5 per mil, and
18OSO4 values range
from 9.7 to 20.7 per mil. Two of the three
18OSO4
values fall outside of the supergene alunite sulfate field and have large
positive
18OSO4-OH values that are not
consistent with a supergene origin (Rye et al., 1992).
Calculated
DH2O values for late-stage alunite
range between 35 and 75 per mil, with
18OH2O
less than 0 per mil. The data fall on a trend that projects, at low
D
values, close to the meteoric water line, and suggests mixing of magmatic fluids
with varying amounts of meteoric water.
Late jarosite: Oxygen isotope data are available only for one sample of
jarosite. The
18OOH data for this
jarosite (Fig. 12B) is slightly heavier than the supergene jarosite OH ("SJOH")
zone as defined by Rye and Alpers (1997), but
18OSO4
data (4.5
) fall within the supergene jarosite sulfate field ("SJSF,"
Rye and Alpers, 1997). These values may indicate limited OH exchange with low pH
water (Rye and Alpers, 1997) or formation from partly exchanged meteoric water.
The composition of fluid calculated in equilibrium with jarosite is close to the
predicted composition of meteoric water and is consistent with formation in the
supergene environment.
 |
Constraints on Alteration and Ore Deposition
|
|---|
The Pascua deposit is a large, complex system that formed from episodic
magmatic-hydrothermal activity recorded over a period of about 1 Myr. The
occurrence of alunite throughout the paragenesis of this deposit is atypical of
most high-sulfidation systems, where advanced argillic alteration and alunite
deposition clearly precede mineralization (e.g., Stoffregen, 1987; Rye 1993;
Cooke and Simmons, 2000, and references therein). Constraints on the evolution
of this system from preore alteration through main-stage ore deposition to
postore processes can be deduced from stable isotope systematics for alunite, in
combination with mineralogical data and 40Ar/39Ar dating
(Fig. 13).

View larger version (35K):
[in this window]
[in a new window]
|
FIG. 13. Schematic section showing the evolution of the Pascua deposit and
associated alteration assemblages over time. A. ~9.1 to 8.8 Ma: widespread,
early hydrothermal alteration in the greater Pascua region created patchy vuggy
silica (silicic) alteration zones and widespread advanced argillic assemblages
(AA1) grading outward to argillic and propylitic (not shown) assemblages.
Sericite at depth is inferred from old core logs and company reports. Minor
steam-heated alteration zones may have developed at or near surface at this
time. B. 8.8 to ~8.4 Ma: main-stage Au-Ag-Cu mineralization followed brecciation
in the Brecha Central area and in several smaller satellite breccia bodies.
Alunite-pyrite-enargite (APE) mineralization was deposited in open spaces in the
breccia matrices and surrounding vein networks. Upwelling magmatic fluids
(condensed magmatic vapors) overwhelmed and displaced meteoric water. Ore
deposition was coincident with lowering of the water table due to regional
uplift and erosion. C. 8.4 to 7.9 Ma: postore alteration and late-stage
processes coincided with the waning stages of the magmatic-hydrothermal system.
Pulses of magmatic vapor deposited magmatic steam alunite near surface above
Brecha Central. Late-stage alunite ± jarosite formed thin veinlets and
disseminations from mixed magmatic-meteoric fluids down to depths of ~600 m
below the present-day surface. Isolated disseminations and veinlets of supergene
jarosite ± alunite (not shown) formed from cooler meteoric waters that
penetrated back into the system.
|
|
Preore alteration
Preore advanced argillic alteration in the Pascua area occurred prior to
about 9.1 Ma. Variations in the clay mineralogy associated with this assemblage
(kaolinite-dickite-pyrophyllite) are interpreted to record a progressive cooling
of the hydrothermal fluids, with increasing elevation and with distance away
from Brecha Central. Similarly, the zonation from advanced argillic and silicic
assemblages outward to surrounding argillic and propylitic alteration indicates
neutralization of hypogene hydrothermal fluids by wall-rock reaction away from
the main fluid conduits. This zonation is typical of many high-sulfidation
systems (e.g., Steven and Ratté, 1960; Stoffregen, 1987; Arribas, 1995). Stable
isotope data for preore alunite are consistent with a magmatic-hydrothermal
origin and suggest that alteration resulted from SO2
disproportionation following the condensation of a magmatic vapor plume (Rye et
al., 1992; Rye, 1993). These vapors were absorbed or condensed into fluids with
magmatic H and O isotope signatures, although limited mixing with exchanged
meteoric water is indicated by the trend of alunite data nearly parallel to the
meteoric water line (see Fig. 10B).
Steam-heated alteration
Steam-heated alunite forms from aqueous sulfate created by the oxidation of H2S
at or above the water table (Schoen et al., 1974). The fluids responsible for
this alunite typically have a meteoric isotopic signature (e.g., Rye et al.,
1992; Ebert and Rye, 1997), however steam-heated alunite at Pascua has
DH2O
values up to 70 per mil higher than those of local meteoric water. Deuterium
enrichments are reported in steam-heated fluids from active geothermal systems,
due to the effects of subsurface boiling and/or evaporation (Truesdell et al.,
1977; Henley and Stewart, 1983), but these are typically much smaller (1020
).
This enrichment is typically on the order of 10 to 20 per mil, however, and is
much less than that recognized at Pascua. The Pascua data, therefore, cannot be
attributed to boiling or evaporation alone, and we interpret these results to
indicate that there was little meteoric water involved in steam-heated
alteration. In fact, much of the ground water seems to be composed of condensed
magmatic fluids that effectively displaced meteoric ground water.
The timing of steam-heated alteration above Brecha Central cannot be
constrained by 40Ar/39Ar methods due to the large overlap
in analytical errors (see Fig. 8) but may have been coeval with either the
development of preore advanced argillic alteration or the main ore stage. We
suggest the latter is more likely because only during the main ore stage would
renewed heat and fluid input, combined with decompression related to the
eruption of multiple breccias, have caused boiling and release of sufficient
volatiles to account for such widespread alteration in the steam-heated zone.
Alunite-pyrite-enargite assemblage mineralization
Deposition of the alunite-pyrite-enargite assemblage with associated Au and
Cu mineralization occurred at about 8.8 Ma. Alunite-sulfide stable isotopes
indicate that alunite-pyrite-enargite mineralizing fluids were of magmatic
origin with a wide range of temperatures from 200° to 350°C. Although abundant
sulfate minerals are recognized during this stage, sulfur isotope data indicate
that the Pascua system, as a whole, was H2S dominant with H2S/SO4
averaging about 1 to 3. This dominance of H2S is typical of many
magmatic-hydrothermal systems (Rye, 1993, in press; Bethke et al., in press). At
Summitville, there is evidence that rock buffering locally produced H2S/SO4
<1 (Bethke et al., in press), however similar shifts in the ratio of sulfur
species have not been recognized at Pascua during alunite-pyrite-enargite
deposition. The bulk of alunite-pyrite-enargite ore precipitated over a narrow
range of redox and pH conditions (Fig. 14), although local variations in the
mineralogy indicate that a broader range of depositional conditions was achieved
locally.

View larger version (42K):
[in this window]
[in a new window]
|
FIG. 14. Log (O2)-pH diagram illustrating condition in
the Pascua alunite-pyrite-enargite (APE) mineralizing event at 275°C and vapor
saturation (59 bars). The stability fields of alunite (alun), enargite (en),
hematite (hem), magnetite (mag), muscovite (musc), pyrrhotite (po), pyrite (py),
tennantite (tenn), and the predominance fields of aqueous sulfur-bearing species
are shown. Shaded and stippled gray areas represent total gold solubility
contours for Au(HS)2, AuHS(aq), and
AuCl2 at 1 ppb, 10 ppb, and 1 ppm, resepectively.
Calculations for the distribution of species and gold concentrations are
modified after Cooke et al. (1996). The white area labeled "APE ore
fluids" and the black area contained within represent the chemical
characteristics of alunite-pyrite-enargite mineralizing fluids. The black area
represents the bulk of alunite-pyrite-enargite ore and is contained within the
stability fields of alunite, pyrite, and enargite and is constrained by H2S/SO4
ratios from 1 to 2.5 (based on 34S results from
this study). Local mineralogical variability (e.g., monomineralic veins of
alunite, enargite, and/or pyrite and the presence of native sulfur) suggests
that the composition of the ore fluid varied over a wider range locally,
schematically represented by the larger white area.
|
|
The isotopic signature of ore-stage alunite and associated sulfides is
similar to preore advanced argillic alteration and indicates that alunite formed
from the disproportionation of SO2 in magmatic vapor (cf. Rye et al.,
1992). During the alunite-pyrite-enargite event, these vapors were focused into
structures and open spaces in the breccia matrices and continued wall-rock
reaction is recorded as restricted alteration halos around veins. Although the
exact conditions of alunite-pyrite-enargite formation are difficult to constrain
due to a lack of fluid inclusion data, the mineralogical assemblages are
consistent with gold transport as AuHS(aq). Continued influx
of heat, volatiles, and fluids into previously condensed magmatic vapors were
most likely responsible for alteration and metal deposition during the
alunite-pyrite-enargite event.
Postore alteration
Late-stage alteration events record the gradual collapse of the Pascua
magmatic-hydrothermal system after about 7.9 Ma. Results from the stable isotope
study indicate that late alunite veins were not generated under typical
supergene conditions, although they have a field appearance suggestive of a
supergene origin. Instead, alunite of this stage formed from mixed magmatic-meteoric
fluids at moderate temperatures of about 100°C. The aqueous sulfate required to
precipitate alunite of this stage was likely derived from the oxidation of
precursor sulfides, given their abundance within the Pascua-Lama system.
Significant silver enrichment may have accompanied this stage of alteration,
although the exact nature and timing of this event remain unclear (Chouinard,
2003).
The occurrence of jarosite with a supergene isotopic signature suggests that
most magmatic activity in the region had ceased by 7.5 Ma, at the latest.
Evidence for the supergene remobilization of gold is observed locally throughout
the Pascua district in oxidized breccias such as Esperanza and Penelope, where
visible gold occurs as coatings on high-level, oxidized fracture surfaces.
 |
Conclusions
|
|---|
The unusual occurrence of coeval alunite-pyrite-enargite with the major stage
of Au, Ag, and Cu mineralization at Pascua constrains the nature and source of
the ore-forming fluids. In most other high-sulfidation systems, ore fluids are
considered to be of mixed magmatic-meteoric origin, resulting from the
adsorption of magmatic vapors or brines by shallow meteoric water (e.g., Rye,
1993, Arribas, 1995; Hedenquist et al., 1998, Cooke and Simmons, 2000, and
references therein; Muntean and Einaudi, 2001). However, at Pascua, there is no
evidence for the involvement of meteoric fluid in the main ore-forming stage.
Alunite-pyrite-enargite ore precipitated from fluids that were acidic, H2S
dominant, and had a wide temperature range from 200° to 350°C. Gold was likely
transported as a bisulfide complex, suggesting that boiling was a viable
depositional process.
The predominance of magmatic fluids is recorded throughout the paragenesis of
the Pascua deposit. Significant meteoric input is only recognized in the
precipitation of supergene jarosite (+ alunite) in the final stages of the
hydrothermal system. Similar conditions have been recorded at the Tambo high-sulfidation
deposit, which formed at nearly the same time in the southern end of the El
Indio-Pascua belt (Deyell et al., in press), suggesting that the same hydrologic
regime was affecting epithermal processes at a regional scale. These findings
are in agreement with published data from Bissig et al. (2002b) that recorded
regional uplift and subsequent pediment erosion throughout this region during
the mid to late Miocene. These events, combined with a transition from semiarid
to arid climatic conditions (e.g., Bissig et al., 2002b, and references
therein), would have caused significant disturbances to the local water table
and contributed to the limited availability of meteoric fluids in the epithermal
environment. As a result, magmatic fluids were largely undiluted for the
duration of magmatic-hydrothermal activity in the Pascua-Lama region.
 |
APPENDIX
|
|---|
All 40Ar-39Ar alunite analyses reported in this study
were carried out at the Queens University laboratory. All samples were step
heated and dates and errors (all given at 2
) are calculated
using formulae given by Dalrymple et al. (1981) and the constants recommended by
Steiger and Jäger (1977). The interlaboratory standard (flux-monitor) used is
Mac-83 biotite (Sandeman et al., 1997), with the most recent published age of
24.36 ±