1999,,,low Grade Metamorphsim of Volcanic Rocks ,,Iceland.02Neuhoff
1999,,,low Grade Metamorphsim of Volcanic Rocks ,,Iceland.02Neuhoff
1999,,,low Grade Metamorphsim of Volcanic Rocks ,,Iceland.02Neuhoff
467–501]
467
468 Philip S. Neuhoff and others—Porosity evolution and mineral
GEOLOGICAL SETTING
Geology.—Teigarhorn is located near the eastern edge of the volcanic plateau of
present-day Iceland (fig. 1). Lavas erupted from fissures and central volcanoes during
magmatism along a divergent plate boundary between the North American and Eur-
asian plates (Kristjánsson, Gudmundsson, and Haraldsson, 1995). The oldest exposed
lavas are only about 13 Ma (Watkins and Walker, 1977). Upper crustal stratigraphy in
eastern Iceland consists largely of subaerial basaltic lavas with minor intercalated
sediments, silicic lavas, lignite layers, and pyroclastic units (Robinson and others, 1982;
Gústafsson, 1992; Walker, 1960b; Wood, 1977; Sæmundsson, 1979). Continued volca-
nism and rifting led to pronounced tilting of the lavas that dominate crustal structure
around Teigarhorn (Pálmason, 1973); dips increase from ⬍2° to 6° southwest at the top
of the lava pile to ⬃6° to 9° southwest at sealevel (fig. 1; Walker, 1974).
Disrupting the monoclinal structure of the lava pile are mafic and silicic intrusive
centers and central volcanic complexes (fig. 1; Walker, 1963). Extrusive products of the
extinct central volcanoes are interlayered with fissure-erupted lavas that make up most of
the crust in eastern Iceland. Glacial exhumation has revealed an alignment of mafic dice
swarms, central volcanoes, and intrusive centers that defines the geometry of paleo-rift
zones. Teigarhorn lies within a prominent dike swarm (called the Álftafjördur dike
paragenesis during low-grade metamorphism of basaltic lavas at Teigarhorn 469
Fig. 1. Generalized geological map of a portion of eastern Iceland (after Jóhannesson and Sæmundsson,
1989) showing the distribution of major rock types, locations of intrusions and central volcanoes, and selected
lava orientations from Walker (1963, 1974). Dashed curve outlines the Álftafjördur mafic dike swarm where
dikes make up greater than 8 percent (by volume) of the crust (Walker, 1963).
swarm; Walker, 1963): the Álftafjördur volcanic center and the Austurhorn intrusion
both lie within this swarm (figs. 1 and 2; Walker, 1963; Blake, 1970). The composition
and mineralogy of these dikes are similar to the surrounding basalts discussed below
(Walker, 1963). The dikes are subvertical, trend north-northeast (fig. 2), and vary in
thickness from 0.25 to 3 m.
Lithology.—The volcanic stratigraphy exposed at Teigarhorn includes approx 150 m
of subaerial mafic lava flows with minor intercalated pyroclastic units. The uppermost
Fig. 2. Map of Teigarhorn showing the distribution of individual mafic dikes (determined by field
mapping and from aerial photographs) and the locations of samples and porosity measurements. Also shown is
the location of the Skessa tuff (Walker, 1963).
Philip S. Neuhoff and others 471
lavas belong to the Grænavatn Porphyritic Group (Walker, 1959), and the oldest units
are the Skessa tuff (see below) and surrounding undesignated lavas. These units were
buried to approx 1500m at Teigarhorn (Walker, 1960b). Compositionally, the lavas
exposed at Teigarhorn (fig. 2) are 3 to 25 m thick subaerial low-MgO tholeiites and
ferrobasalts (Flower and others, 1982). Most flows are fine-grained, containing abundant
groundmass plagioclase, olivine, pyroxene, and glass with minor olivine and pyroxene
phenocrysts. Lavas with up to 30 modal percent plagioclase phenocrysts occur infre-
quently. The flows have scoraceous, brecciated flow tops. Vesicles are abundant near the
top and bottom of the flows, but the flow centers are essentially non-vesicular and
massive. A 20 m thick, green, welded, silicic tuffaceous unit (the Skessa tuff of Walker,
1959) is exposed near the southern portion of Teigarhorn (marked ‘‘Skessa tuff’’ on fig.
2). The Skessa tuff is a regionally extensive unit that appears to be the product of an
eruption from an unexposed volcanic center approx 25 km north of Teigarhorn (Walker,
1962). Separating many flows are undulatory weathering surfaces that frequently are
overlain by bright red basaltic tephras that are 1 to 15 cm thick and constitute less than 1
percent of the stratigraphy. The tephras are phreatomagmatic deposits produced by
interaction between basaltic magmas and surface or groundwaters at eruption sites 10’s
to 100’s of km away from Teigarhorn (Heister, ms).
The volcanic stratigraphy in eastern Iceland has been subjected to regional zeolite-
facies metamorphism; higher grades of metamorphism are exposed around intrusive
complexes and central volcanoes (Walker, 1960b). Walker (1960b) described depth-
controlled zeolite-facies mineral zones (fig. 3) that parallel the paleosurface of the lava
pile and transgress the volcanic stratigraphy. This implies that the zones formed after
tilting of the lava pile and do not solely reflect variations in chemistry of the host lavas.
The zeolite zones described by Walker (1960b) are underlain by zones characterized by
stilbite ⫾ heulandite ⫾ mordenite, laumontite, and prehnite ⫹ epidote with increasing
depth (Mehegan, Robinson, and Delany, 1982). Coastal exposures around Berufjördur
(fig. 1) have provided museum-grade zeolite specimens for over 200 yrs (Betz, 1981).
Teigarhorn, in particular, has been an important source for zeolites used in studies of
crystal chemistry (Galli and Rinaldi, 1974; Passaglia, 1975; Passaglia and others, 1978;
Alberti, Pongiluppi, and Vezzalini, 1982), crystal structures (Slaughter and Kane, 1969),
reaction kinetics (Ragnarsdóttir, 1993), optical properties (Bauer and Hrı̂chová, 1962),
and piezoelectricity (Bellanca, 1940).
METHODS
Mapping and sample collection.—Representative samples of altered and unaltered lavas
were collected throughout the coastal exposures at Teigarhorn in conjunction with
mapping of dike distribution and outcrop porosity measurements (fig. 2). Estimates of
total macroscopic porosity (herein defined as the percentage of a given volume of rock
occupied by pore space) prior to mineral precipitation were conducted in situ during
mapping following the methods of Manning and Bird (1991). In the present communica-
tion, a distinction is made between the total porosity at a given point in time (⌽) and the
cumulative porosity that reflects pore abundance (⌽*) but does not include the effects of
pore occlusion by mineral precipitation (Manning and Bird, 1991). The contribution of
primary pores to ⌽* (⌽1*) was estimated by placing a clear mylar sheet imprinted with a
54 by 41.5 cm grid (node spacing 1.5 cm) directly on the outcrop as a basis for point
counting macroscopic pores. Estimates of ⌽1* included vesicles and other primary
structures as well as open space between flow units that could not texturally be linked to
post-eruption processes. Although this technique only identifies macroscopic pores on
the order of 1 mm or larger, the contribution of these pores to ⌽1* is large, and
uncertainties are on the order of 1.5 percent or less for grids of over 1000 points
(Manning and Bird, 1991; Chayes, 1956). The contribution of secondary pore spaces to
472 Philip S. Neuhoff and others—Porosity evolution and mineral
Fig. 3. Depth distribution of regional zeolite facies mineral zones in eastern Iceland. Zones above sealevel
are from mapping in olivine normative tholeiites around Berufjördur by Walker (1960b) and below sealevel are
from studies of the Reydarfjördur drill core located ⬃40 km north-northwest of Teigarhorn (Mehegan,
Robinson, and Delaney, 1982). Note that the outcrops at Teigarhorn are approximately at sealevel. Stage II
alteration at Teigarhorn includes both the mesolite/scolecite and heulandite ⫹ stilbite mineral zones.
Mineral Paragenesis
Alteration at Teigarhorn occurred during three paragenetic stages. The textural and
mineralogical characteristics of these stages are summarized in table 1 and the sequence
of mineral precipitation during each stage is shown in figure 4. Representative composi-
tions and molar formulas for clays, zeolites, and feldspars are given in tables 2 through 5.
Stage I.—Paragenetic stage I is the earliest alteration at Teigarhorn and is character-
ized by linings of primary pore space by celadonite and silica minerals (table 1; figs. 4
and 5B). All other mineral assemblages are either precipitated on top of or cross cut Stage
I alteration (fig. 5B). Stage I alteration is common but not present in all samples.
Celadonite and silica minerals typically form layers along the walls of pore space with a
layered or botryoidal habit. Occasionally Stage I silica forms layered geopetal structures
TABLE 1
Mineralogical and textural characteristics of paragenetic stages observed at Teigarhorn,
eastern Iceland
Mineral
Paragenetic Mineral Assemblages Assemblages
Stage* in Pore Space in Matrix Textures
Stage I ⫾celadonite ⫾ silica minor celadonite pores: lining primary pore
space; geopetal structures
matrix: minor celadonite
replacing groundmass
Stage II C/S2 ⫹ scolecite ⫾ chalcedony ⫾ C/S ⫹ albite ⫹ pores: C/S pore linings
quartz zeolites
C/S ⫹ stilbite ⫹ heulandite ⫾ partial or complete primary
mordenite ⫾ epistilbite ⫾ pore fillings by zeolites
qu artz ⫾ chalcedony
matrix: pseudomorphic
replacement of groundmass
and phenocrysts
StageIII quartz ⫾ clinoptilolite ⫾ stilbite ⫾ quartz ⫹ C/S pores: brecciation, fracturing,
heulandite ⫾ mordenit e ⫾ mineral precipitation
scolecite ⫾ laumontite ⫾ calcite
matrix: discoloration (SiO2
metasomatism)
*Stage I was the earliest stage, Stage III the last. 2Mixed-layer chlorite smectite (see text).
474 Philip S. Neuhoff and others—Porosity evolution and mineral
Fig. 4. Temporal sequence of secondary minerals developed in lavas at Teigarhorn. Events toward the
right of the diagram occurred after events to the left.
on the bottom of vesicles that are oriented roughly parallel to the lava stratigraphy.
Minor celadonite occurs in the volcanic groundmass.
Celadonite typically precedes silica in the crystallization sequence (fig. 5B). A
representative celadonite composition from Teigarhorn is given in table 2, and 6
celadonite analyses are plotted in figure 6. Stage I celadonite is close to the endmember
compositions as defined by Buckley and others (1978), with only slightly lower potas-
sium content and higher aluminum than the ideal endmember composition
[K2(Fe,Al)2(Mg,Fe)2Si8O20 (OH)4; table 2; fig. 6]. The slight deviation from ideal compo-
sition may be due to the presence of a small trioctahedral smectite component (either as
intergrowths or interstratifications) in the celadonites, although smectite is not detected
in the XRD patterns. Celadonite is distinguished from other phyllosilicates by its
distinctive bright green color in both hand specimen and thin section (figs. 5B and 5D),
its layered habit, and by its potassium-rich composition (table 2; fig. 6). Silica minerals
consist of quartz and/or chalcedony. It appears likely that quartz and chalcedony
replaced preexisting Stage I opaline silica during burial, as opaline silica is the only Stage
I silica phase observed near the top of the lava pile above Teigarhorn (Neuhoff,
unpublished data).
Stage II.—Paragenetic Stage II alteration is characterized by continued infilling of
primary pore spaces and partial replacement of groundmass and phenocryst minerals.
As with Stage I, the distribution at the outcrop scale of Stage II alteration is largely a
function of the distribution of primary porosity in the lavas (fig. 5A). Stage II mineral
assemblages and paragenetic sequences are summarized in table 1 and figure 4. Mixed
layer chlorite/smectite (C/S) is the first phase to precipitate in primary pore space during
Stage II (fig. 5B, C). The term ‘‘mixed layer chlorite/smectite’’ and the abbreviation C/S
are used in the present communication without specific inference of ordering among
layers (Reynolds, 1988). Brown to brownish green radiating bundles of C/S for linings of
primary pore space. The bundles nucleated on pore walls or on Stage I minerals where
present and developed toward the center of the pores. As shown in figure 6, Stage II C/S
paragenesis during low-grade metamorphism of basaltic lavas at Teigarhorn 475
compositions (table 2) project near the binary compositional join between ferromagne-
sian trioctahedral smectite (saponite-ferrosaponite; Güven, 1988) to ferromagnesian
chlorite (clinochlore-chamosite; Bailey, 1988). The percentage of chlorite (or smectite)
interlayers was estimated from electron microprobe analysis of the total molar Si ⫹ Al ⫹
Mg ⫹ Fe content of the mineral by linear interpolation between smectite and chlorite
endmembers (Schiffman and Fridleifsson, 1991). This interpolation requires calculation
of molar formulas based on 28 oxygen charge equivalents (O ⫹ 2OH), for which the
molar Si ⫹ Al ⫹ Mg ⫹ Fe contents of endmember chlorite and trioctahedral smectite are
20 and 17.82, respectively. By this analysis, Stage II clays contain between 12 and 85 mol
percent chlorite. The most prominent interlayer cations in C/S clays are Ca⫹2 and Na⫹
(table 2). Potassium typically accounts for less than 5 mol percent of the total interlayer
cations but in rare cases constitutes more than 25 mol percent with one sample as high as
60 mol percent. Mixed-layer chlorite-smectite with high K⫹ contents always occur in
close proximity to celadonite, and textural evidence suggests that local replacement of
Stage I celadonite by Stage II C/S may have occurred.
The compositions of C/S clays become progressively enriched in chlorite during
Stage II alteration. Figure 7 shows the percentage of chlorite in a rim of Stage II C/S as a
function of distance from the vesicle wall. In general, C/S near the wall of the vesicle
precipitated earlier than that closer to the center of the vesicle. It can be seen in figure 7
that the chlorite content of the clay generally increases from ⬃20 percent at the vesicle
wall to about 80 percent near the center of the vesicle, indicating that the chlorite content
of Stage II C/S increased with time.
The last phases to crystallize during Stage II are minor quartz and chalcedony along
with zeolites (table 1 and fig. 4) that serve as index minerals for the regional zeolite zones
shown in figure 3. Zeolites partially or completely fill the pore space remaining after
precipitation of C/S (fig. 5B). As indicated in table 1, two Stage II zeolite assemblages
occur at Teigarhorn; the mesolite ⫹ scolecite zone described by Walker (1960b; fig. 3) is
underlain by a zone characterized by heulandite ⫹ stilbite with occasional mordenite
and epistilbite. The lower zone is similar to zones characterized by heulandite and/or
stilbite identified in Icelandic geothermal systems (Kristmannsdóttir and Tómasson,
1978), East Greenland (Neuhoff and others, 1997), the Faeroe Islands (Jørgensen, 1984),
the North Shore Volcanic Group (Minnesota, USA; Schmidt, 1990), and the Reydarfjör-
dur drillhole (Mehegan, Robinson, and Delaney, 1982).
Mesolite ⫹ scolecite zone alteration at Teigarhorn is characterized by radiating
bundles of near stoichiometric scolecite fibers (table 3) that fill nearly all pore spaces
remaining after C/S precipitation (fig. 5D). Scolecite is white in hand specimen and often
has a brownish hue in plane polarized light. Growth direction is very distinct; the
bundles widen as a function of distance from the apex where nucleation first occurred at
the edge of Stage II C/S pore linings.
Typical vesicle fillings in the heulandite ⫹ stilbite zone consist of intergrowths of
stilbite, heulandite, and mordenite with occasional epistilbite (fig. 5B, C). Frequently
only one or two of these minerals are present in a given pore (usually stilbite or
heulandite). Mordenite is fibrous, occurring alone or as discrete fibers within heulandite
or stilbite crystals (fig. 5C). Stilbite, heulandite, and epistilbite occur as euhedral blocky
to lath shaped crystals. No evidence was found for heulandite ⫹ stilbite zone assem-
blages overprinting (replacing) mesolite ⫹ scolecite zone minerals. Representative
compositions of Stage II stilbite, heulandite, and epistilbite are given in table 3 (Stage II
mordenite is too fine-grained for reliable analysis). Epistilbite compositions are slightly
richer in Al⫹3 and Na⫹ than the ideal formula (CaAl2Si6O16 · 6H2O). The ratio of Si to Al
in Stage II heulandite and stilbite ranges from ⬃2.9 to ⬃3.4 and ⬃2.9 to ⬃3.5,
respectively. Exchangeable cation sites in Stage II heulandite (fig. 8) and stilbite (fig. 9)
476 Philip S. Neuhoff and others—Porosity evolution and mineral
Fig. 5. Stage I and II alteration at Teigarhorn. (A) Outcrop view showing the flow scale distribution
of primary pore spaces that host Stage I and II alteration. Two lava flows are visible separated by an undulating
surface. The top of the lower flow is highly vesicular and contains abundant Stage I celadonite rims and Stage II
scolecite visible as white pore fillings. The massive flow center of the upper flow contains only isolated,
large vesicles filled with Stage II scolecite. (B) Stage I celadonite and silica rims preceding Stage II mixed-
layer chlorite/smectite (C/S) and stilbite in sample ISMP-2. Photomicrograph taken with polars partially
crossed. Analyses in columns 2 and 7 of table 2 were conducted on celadonite and C/S rims shown in photo.
paragenesis during low-grade metamorphism of basaltic lavas at Teigarhorn 477
(C) Stage II heulandite ⫹ stilbite zone assemblages (sample 94-53). Photomicrograph taken with polars
partially crossed. Vesicle in rimmed with celadonite (Stage I) and Stage II C/S. Filling the pore are lath-shaped
crystals of stilbite (birefringement phase) and blocky crystals of heulandite. The stilbite and heulandite crystals
contain abundant mordenite fibers. Blue areas are residual porosity (color due to impregnated epoxy). (D)
Stage II mesolite ⫹ scolecite zone assemblages (sample 94-77). Photomicrographs taken with polars partially
crossed. Vesicles are rimmed with C/S and filled with scolecite. The large plagioclase phenocrysts have been
replaced by scolecite.
478 Philip S. Neuhoff and others—Porosity evolution and mineral
are dominantly occupied by Ca⫹2, although in some samples Na⫹ and K⫹ constitute up
to 30 percent of the exchangeable ions.
Stage II minerals in the matrix of the lavas occur predominantly as pseudomorphic
replacements of primary phenocryst and groundmass phases (table 1). Groundmass
glass, pyroxene, and olivine as well as pyroxene and olivine phenocrysts are partially
replaced by C/S with compositions similar to that lining primary pore space. Ground-
mass and phenocryst plagioclase is pseudomorphically replaced by scolecite and/or
completely albitized (table 4).
Stage III.—The latest alteration at Teigarhorn (Stage III) is developed in secondary
pore space (fractures, breccias, et cetera) and primary pores intersected by this porosity
(figs. 10, 11). Stage I and II minerals filling vesicles and other primary pore spaces are
cross cut by Stage III porosity. Stage III breccias commonly entrain Stage I and II
minerals (fig. 11C). This is frequently visible in outcrop as rafts of vesicle walls coated
with celadonite and Stage II heulandite within Stage III breccias. Development of
fracture porosity during Stage III is associated with intrusion of mafic dikes that cross cut
the lavas. The near-vertical dikes make up more than 8 percent of the volume of rocks at
Teigarhorn (fig. 2) and vary in thickness between 0.5 and 3 m. Spacing between dikes is
uneven, ranging from less than 1 to over 30 m. Lava flows within a few meters of the
dikes are extensively fractured, but the fractures generally do not contain secondary
TABLE 2
Representative phyllosilicate analysis from Teigarhorn, eastern Iceland determined by EPMA
paragenesis during low-grade metamorphism of basaltic lavas at Teigarhorn
Fig. 7. Temporal changes in the compositions of Stage II mixed-layer chlorite-smectite (C/S) depicted as a
function of percent chlorite (calculated using method of Schiffman and Fridleifsson, 1991). Analyses are from a
300 µm linear traverse through a Stage II clay rim of a vesicle in sample 94-50. The compositions of C/S trend
from early smectitic clays (near the vesicle wall) to relatively late chloritic clays (towards the center of the
vesicle).
minerals. At distances greater than 3 m from some dikes, brittle deformation features
host Stage III mineral assemblages. This alteration is only locally developed.
Stage III mineral assemblages fill 0.2 to 2 cm thick veins oriented parallel to the
dikes and irregular fractures and breccias up to 0.8 m across. Veins and fractures can be
traced for up to 10 m. The margins of veins are occasionally lined with 5 micron equant
crystals of clinoptilolite (compositions given in fig. 8 and table 5) and usually marked by
one or more layers of chalcedony. Filling the veins is quartz that is an- to subhedral near
the margin and euhedral in the center with the c axis oriented normal to the vein walls.
Clusters of euhedral quartz crystals (individual crystal size two to four mm) frequently
form stalactiform deposits hanging vertically from the roofs of large cavities in the
breccias. Generally chlorite-rich C/S occurs intergrown with quartz and as discrete
fracture fillings (table 2). The last assemblages to fill the veins and breccias are euhedral
stilbite, heulandite, or stilbite ⫹ heulandite ⫹ mordenite (fig. 11A, B) with local
development of later scolecite, laumontite, and calcite in large, open cavities (fig. 11B).
Stage III zeolites frequently occur as vein fillings and as large (cm-scale) euhedral crystals
for which Teigarhorn is famous. Representative compositions of Stage III zeolites are
given in table 5. The compositions of heulandite, stilbite, and scolecite formed during
Stage III are similar to those formed during Stage II (table 5; figs. 8 and 9).
Fluid inclusion heating/freezing experiments were conducted on a 3 cm crystal of
Iceland spar type calcite that formed late in Stage III along with scolecite and laumontite.
Two-phase inclusions (liquid with vapor bubbles constituting ⬃5 percent of the inclusion
paragenesis during low-grade metamorphism of basaltic lavas at Teigarhorn 481
TABLE 3
Representative compositions of Stage II zeolites from vesicles determined by EPMA
volume) occur throughout the crystal. Most inclusions are up to 0.25 mm in diameter
with ovoid, elongate, or triangular morphologies. Primary inclusions homogenized at
115° ⫾ 8°C and frozen at ⫺0.1° to 0°C.
Stage III fractures and breccias are frequently surrounded by aureoles (1-40 cm
wide) of extensive SiO2 metasomatism (fig. 10). These zones are lighter in color than the
surrounding lavas, and in extreme cases are pale yellow. The lava is extensively replaced
by fine-grained aggregates of quartz with about 5 modal percent C/S (fig. 11C and D).
Stage III C/S (table 2, fig. 6) is more chloritic than in Stage II, ranging from 55 to 90
percent chlorite. Where this alteration is most intense the original texture of the lava is
completely obliterated. The replaced lava is highly porous with a spongy texture (fig.
11D) and has approx 20 percent secondary porosity.
Bulk Rock Chemistry
Bulk rock compositions of samples from two successive lava flows (and an interven-
ing mafic volcanoclastic unit) from Teigarhorn (samples 94-86, 94-87, 94-90, 94-91,
94-96, and 94-101; fig. 2) are given in table 6. Samples 94-86 and 94-101 are from the
overlying (younger) flow, 94-87 is the mafic volcanoclastic horizon, and 94-90, 94-91,
and 94-96 are from the underlying (older) flow (see footnote to table 6 for sample
descriptions). The least altered samples of both flows (94-86 and 94-90) are similar in
composition and belong to the low MgO tholeiite lavas described by Flower and others
(1982). The mafic volcanoclastic unit is compositionally similar to the unaltered lavas
except for higher Fe, Si, K, Rb, Cu, Cr, and Ni and lower Ca, Na, P, Sr, and V contents in
the oxide-rich, more oxidized volcanoclastic horizon.
The effects of alteration on bulk rock composition can be assessed by comparing
these relatively unaltered samples with examples of Stages I and II (94-91 and 94-96) and
Stage III (94-101) alteration from the same flows. This is illustrated by the isocon plots
482 Philip S. Neuhoff and others—Porosity evolution and mineral
Fig. 8. Molar proportions of Ca, Na, and K in heulandite and clinoptilolite from Teigarhorn.
(Grant, 1982) in figure 12, which compare the compositions of sample 94-96 (which
exhibits Stage I and II alteration) and sample 94-101 (extreme Stage III SiO2 metasoma-
tism) to those of the relatively unaltered samples 94-90 and 94-86, respectively, of the
same flows. The linear trends (labeled ‘‘constant mass isocon’’ in fig. 12) represent no
elemental loss or gain during alteration. It can be seen in figure 12A that Stage I and II
alteration together lead to only minor changes in bulk rock composition, as illustrated by
the fact that the major element data plot along the isocon. With the exception of Rb and
La in sample 94-96, the other trace element data listed in table 6 have been omitted for
clarity, but they all plot along the constant mass trend as well. The most significant bulk
chemical signature of Stages I and II is an increase in La, Rb, and to a lesser extent K2O
concentration; Rb and K2O were likely concentrated due to precipitation of celadonite
in 94-96. The fact that most of the data in figure 12A plot along the constant mass isocon
reflects the isochemical nature of the hydration reactions that typify Stage I and II
alteration. In contrast, the extreme metasomatism associated with Stage III is clearly
reflected in sample 94-101. Quartz along with minor stilbite and chlorite have entirely
replaced the lava (94-86) resulting in a bulk rock composition that is ⬃92 percent SiO2.
As can be seen in figure 12B, this sample is clearly enriched in SiO2 and depleted in all
other analyzed elements relative to the unaltered lavas.
Porosity
Lava flows at Teigarhorn exhibit values of ⌽1* (volume percent of primary pores) up
to 22 percent with subsequent processes during Stage III leading to values of ⌽2* of up to
11 percent (table 7). The range of values of ⌽1* given in table 7 is typical of basaltic lavas
paragenesis during low-grade metamorphism of basaltic lavas at Teigarhorn 483
reported by Christensen and Wilkens (1982) from the Reydarfjördur drillhole in eastern
Iceland. Variations in ⌽1* (table 7) reflect differences in pore density both between flows
and within various portions of a given flow. The estimates of ⌽2* (volume percent of
secondary pores) are similar to those obtained by Manning and Bird (1991) for fractured
basalts near the margin of the Skaergaard intrusion in East Greenland.
Percentages of ⌽1* occupied by Stage I celadonite and silica and Stage II clays and
zeolites along with the percent of ⌽1* remaining after precipitation of Stage I and II
minerals are given in table 8. Inspection of the table shows that the extent of porosity
reduction during surficial alteration and burial metamorphism is heterogeneous between
samples and is a strong function of paragenetic stage. Precipitation of Stage I minerals,
where developed, leads to a 10 to 30 percent reduction in ⌽1*. Chlorite/smectite
interlayer clays are more common than Stage I minerals and fill 0 to 80 percent of ⌽1*.
Zeolites, the last minerals to precipitate during Stage II, fill an additional 0 to 90 percent
of ⌽1*. The percent of ⌽1* remaining after Stage II varies between ⬃0 and 35 percent; that
is, ⌽ after Stage II ranges from about 0 percent to nearly 8 percent.
DISCUSSION
Time-space development of mineral parageneses.—Mineral parageneses developed at
Teigarhorn present a time-integrated view of the petrogenesis of zeolite-facies basaltic
aquifers in eastern Iceland. In this context, the temporal development of mineral
assemblages can be directly related to the geologic evolution of the crust in eastern
Iceland. This is summarized in figure 13, which depicts the history of the lavas at
484 Philip S. Neuhoff and others—Porosity evolution and mineral
TABLE 4
Representative compositions of altered and fresh plagioclase
determined by EPMA
Teigarhorn. The horizontal axis in figure 13 represents time after eruption and the
vertical axis is crustal depth. Also illustrated in figure 13 are petrotectonic events (lava
eruption, burial during crustal accretion, dike intrusion, glacial exhumation) and mineral
parageneses experienced by the lavas.
Celadonite and silica rims (Stage I) formed at shallow depths (fig. 13). Stage I
geopetal structures at Teigarhorn have orientations similar to the lava flows, and in other
localities exhibit progressively more shallow dips between silica bands due to progres-
sive tilting (Walker, 1974). This relationship suggests that the geopetal structures formed
prior to tilting. Near surface genesis for Stage I minerals is consistent with the low
temperatures (⬍ ⬃50°C) at which celadonite forms in natural systems (Odin, 1988) and
the supersaturation of Icelandic surface waters with respect to chalcedony (Gı́slason and
Eugster, 1987b; Gı́slason, Arnórsson, and Ármannsson, 1996). Silica and potassium are
among the most mobile elements in near-surface environments in Iceland (Gı́slason,
Arnórsson, and Ármannsson, 1996), providing components necessary for formation of
celadonite and opal.
paragenesis during low-grade metamorphism of basaltic lavas at Teigarhorn 485
Fig. 10. Stage III fractures filled with quartz (outcrop view). Pencil for scale. The dashed white lines are
boundaries of metasomatic zones containing quartz ⫹ minor C/S that replace the basalts (see fig. 11C and D).
The aureoles are visible as areas of yellow alteration marked by the symbol ‘‘Si.’’ Also visible are vesicles
containing Stage II scolecite.
Fig. 11 (continued)
Teigarhorn is not present at elevations ⬎100m in lavas exposed in the mountains around
Berufjördur (Neuhoff, unpublished data). This suggests that albitization occurred concur-
rently with Stage II zeolite formation. The paucity of absolute age dates in eastern
Iceland do not permit concise constraints on the timing and duration of zeolite
488
TABLE 5
Representative compositions of Stage III fracture-filling zeolites determined by EPMA.
Philip S. Neuhoff and others—Porosity evolution and mineral
paragenesis during low-grade metamorphism of basaltic lavas at Teigarhorn 489
TABLE 6
Bulk rock compositions from Teigarhorn, eastern Iceland
‡Sample descriptions:
94-86: compact lava with minor Stage III heulandite veins
94-87: bright red mafic volcanoclastic unit marking boundary between two lava flows
94-90: compact, massive lava with minor small vesicles filled with smectite
94-91: compact basalt with minor vesicles containing Stage II clay minerals
94-96: vesicular lava with abundant Stage I celadonite, Stage II heulandite and smectite
94-101: yellow, SiO2-metasomatized lava from aureole around Stage III breccia
*Weight percent of oxide in anhydrous (calcined) sample. †Total Fe reported as Fe203. §Determined by loss on
ignition at 110°C (H2O⫺) and 900°C (H2O⫹). †Parts per million of element in untreated sample. B/D: below
detection limit.
precipitation during Stage II; however, it seems likely that it did not last more than a few
Ma, as observed in East Greenland (Neuhoff and others, 1997).
Pore structures filled with Stage III alteration cross-cut Stage II mineralization. Stage
III is spatially and, presumably, temporally linked to dike intrusion (fig. 13). The spatial
association of Stage III with late dikes indicates that this last paragenetic stage is a local,
and not a regional, phenomenon. The close association with dike intrusion indicates that
Stage III occurred while the region around Teigarhorn was still undergoing intrusive
490 Philip S. Neuhoff and others—Porosity evolution and mineral
Fig. 12. Isocon plots (Grant, 1982) depicting the variation in bulk rock chemistry accompanying Stages I
and II (A) and Stage III (B) alteration at Teigarhorn. The lines labeled ‘‘constant mass isocon’’ depict the
theoretical distribution of data assuming that the mass and compositions of the lavas do not change during
alteration. The SiO2, Rb, and La analyses have been multiplied (by factors of 0.25, 10000, and 10000,
respectively) to facilitate comparison with data for other major oxide components. Data used to construct this
plot are listed in table 6.
paragenesis during low-grade metamorphism of basaltic lavas at Teigarhorn 491
TABLE 7
Outcrop porosity measurements
1Station
locations given in figure 2. Descriptions: (1) flow top; (2) flow center with abundant Stage III fractures
developed ⬃10 m from dike; (3) flow top; (4) flow center; (5) flow center; (6) flow top; (7) flow center with
abundant Stage III fractures developed ⬃15 m from dike; (8) flow bottom with abundant Stage III fractures
developed about 5 m from dike.
TABLE 8
Percent of pore space occupied by Stage I and II minerals based on thin section area measurement.
All samples are from flow tops
aNumber of individual pores measured. bRims of paragenetic Stage I celadonite or silica on vesicle walls.
cIncludes all paragenetic Stage II trioctahedral smectite and chlorite interlayer clays. dStage II zeolite
assemblages, dominantly mesolite ⫹ scolecite and heulandite ⫹ stilbite ⫹ mordenite. ePrimary pore space (⌽1*)
remaining after paragenetic Stage II.
activity related to rifting. Increased heat flow and permeability due to fracturing during
dike emplacement initiated hydrothermal fluid circulation that led to the extensive
wall-rock replacement and fracture filling characteristic of Stage III mineralization. No
age dating has been conducted on the dikes at Teigarhorn that would permit constraints
on the absolute timing of Stage III alteration.
Thermobarometric conditions during Stages II and III.—Barometric conditions during
paragenetic Stages II and III were complex functions of the burial and hydrogeologic
492 Philip S. Neuhoff and others—Porosity evolution and mineral
Fig. 13. Spatial and temporal development of pore-filling mineral assemblages at Teigarhorn. The vertical
axis depicts depth below land surface at the time of each event depicted in the figure. Time elapsed after
eruption increases to the right. No scale is implied on the horizontal axis. The parallel curves on the figure are
boundaries of the lavas exposed at Teigarhorn in space and time. Note that the timing of dike emplacement is
meant to infer only the timing of dikes associated with Stage III alteration at Teigarhorn.
grades during Stage II in eastern Iceland. The relatively high chlorite contents of C/S
associated with the zeolite assemblages at Teigarhorn may be a consequence of albitiza-
tion of primary plagioclases, which promotes chlorite formation through release of
aluminum (Shau and Peacor, 1992; Schmidt and Robinson, 1997). As noted by Schiff-
man and Fridleifsson (1991), the temperatures at which C/S occurs in geothermal
systems are often higher than in regional diagenetic systems. In contrast, temperature
estimates based on zeolite assemblages and fluid inclusion thermometry are similar (see
below). This suggests that estimation of temperatures from observations of zeolite
assemblages in geothermal systems may be more reliable than similar estimates based on
phyllosilicate distribution.
Assuming a temperature of formation for Stage II zeolite assemblages of 90° ⫾ 10°C
(see above) leads to an estimated thermal gradient of approx 53° ⫾ 7°C/km (given a
surface temperature of 10°C). A thermal conductivity for water-saturated basalt of 1.5 to
2.0 W/(mK) (Oxburgh and Agrell, 1982) gives an estimated heat flow during Stage II
zeolite precipitation of 1.9 to 3.2 heat flow units (HFU ⫽ 10⫺6 cal · cm⫺2 · s⫺1). This value
is consistent with conditions on the flank of the present-day active volcanic zone in
Iceland (Flóvenz and Sæmundsson, 1993).
Mineral assemblages and fluid inclusion measurements indicate that temperatures
during Stage III were slightly higher than in Stage II. Empirically-based chlorite
geothermometry (Cathelineau, 1988; Cathelineau and Nieva, 1985) suggests that the
relatively early, chlorite-rich Stage III C/S formed at temperatures between 140° and
180°C. It should be noted that temperatures estimated using chlorite geothermometry
generally disagree with measured temperatures in Icelandic geothermal systems (Schiff-
man and Fridleifsson, 1991). Late Stage III zeolite assemblages are broadly similar to
those in Stage II but may have formed at slightly higher temperatures based on the
presence of laumontite. Laumontite becomes prominent at temperatures above 120°C in
Icelandic geothermal systems (Kristmannsdóttir and Tómasson, 1978), consistent with
pressure corrected homogenization temperatures for late-Stage III calcite (125° ⫾ 8°C
assuming hydrostatic pressures). The apparent decrease in temperature during Stage III
may reflect dissipation of heat flow associated with crystallization of the dikes or may be
the result of errors associated with the chlorite geothermometer noted above.
Chemical fluxes during alteration.—The widespread precipitation of phyllosilicates,
silica minerals, and zeolites in pore spaces during Stages I and II indicates that Si, Al, Ca,
Na, K, Fe, and Mg were mobile over scales of a millimeter or more. This is further
illustrated by the textural preservation of pore walls and by successive lining of pore
spaces during mineral paragenesis, which require that chemical components necessary
for precipitation of secondary minerals were transported into the pores. All components
necessary to form the mineral assemblages described above were readily available in the
primary phases in the lavas (glass, olivine, plagioclase, pyroxene) with the exception of
H2O, which was provided by the metamorphic fluids. The similarity between the bulk
rock compositions of lavas altered during Stages I and II and those of unaltered lavas
(table 6; fig. 12A) indicates that many elements were relatively immobile (that is,
conserved) on the scale of individual lava flows during regional alteration. An exception
is potassium, which is concentrated during Stage I in celadonite (table 6). As noted
above, potassium is very mobile during near-surface alteration of basaltic lavas (Gı́sla-
son, Arnórsson, and Ármannsson, 1996). The concentration of Rb, K2O, and La
apparent in figure 12A and table 6 is consistent with the observations of Wood (1976)
that these elements were mobilized during zeolite-facies metamorphism.
In contrast, high fluid fluxes through some Stage III (see below) fractures resulted in
widespread leaching at the outcrop scale of all elements except silicon (table 6; fig. 12B).
Extensive replacement of the lavas by quartz (plus minor chlorite-rich C/S) is suggestive
of interaction with acidic fluids. The low pH of these fluids may reflect the presence of
494 Philip S. Neuhoff and others—Porosity evolution and mineral
relatively large amounts of dissolved magmatic gases. The solubilities of primary phases
in the lavas as well as C/S and zeolites are all enhanced under very acidic conditions. In
contrast, quartz solubility is insensitive to pH under acidic conditions, and it formed
early during Stage III from SiO2 released during hydrolysis/dissolution of primary
phases. The large model abundance (⬎90 percent) of quartz in Stage III aureoles
suggests that some SiO2 may have been added during alteration. Aside from minor C/S
precipitated in the aureoles and later in veins, no mineralogical sinks for iron and
magnesium formed during Stage III, suggesting that these elements were transported
over distances of meters (the scale of outcrops at Teigarhorn) away from their sources.
Late precipitation of zeolites and calcite in Stage III veins most likely occurred as the
aqueous solutions became more alkaline as protons were consumed by hydrolysis.
Both Stages II and III exhibit a pronounced temporal change during paragenesis
from phyllosilicates ⫹ quartz to assemblages dominated by Ca-zeolites. In both cases,
this paragenetic trend was probably coincident with changes in temperature and
pressure. However, the abrupt shift in the mineral composition from C/S to zeolites must
reflect the chemical evolution of the aqueous electrolyte solutions from which these
minerals precipitated, as schematically illustrated in the phase diagram in figure 14. The
arrows A and B in figure 14 represent two possible trends in fluid composition that would
result in mineral paragenesis of Ca-zeolites after C/S (figs. 4, 5B and C). The decrease in
log [(aMg⫹2/(aH⫹)2] suggested by paragenetic trends A and B would result from decreasing
the activity of Mg⫹2 in the pore fluid as a consequence of precipitation of C/S
(represented by clinochlore in the phase diagram). Mineralogical sources of Mg⫹2
(olivine, basaltic glass) in basaltic lavas are generally less stable than plagioclase under
low-grade metamorphic conditions (Gı́slason and Eugster, 1987a; Gı́slason and Arnórs-
son, 1993) and are often hydrolyzed before plagioclase. Albitization of plagioclase was
likely contemporaneous with precipitation of zeolites in pores during Stage II (see
above), adding Ca⫹2 to aqueous solution and increasing log ([aCa⫹2)/(aH⫹)2] as suggested
by trend A in figure 14. Instability of zeolites relative to C/S during early Stage II
alteration provides an explanation for the apparent lack of zeolite formation during
burial of the lavas inferred from our petrographic observations that heulandite ⫹ stilbite
zone assemblages do not replace mesolite/scolecite zone assemblages. These observa-
tions, combined with the phase relations depicted in figure 14, suggest that the lavas were
buried before olivine and basaltic glass were exhausted, resulting in fluid compositions
during burial in which C/S was stable relative to zeolites.
The general paragenetic trend from C/S to zeolites, as well as fluid composition
trends A and B in figure 14, is consistent with the results of experimental simulations of
basalt-meteoric water interaction at low temperatures. Ghiara and others (1993) con-
ducted experiments on closed-system dissolution of basaltic glass in deionized water at
200°C that resulted in a mineral paragenesis of early smectite ⫹ phillipsite followed by
analcime. Total Mg concentrations in solution during these experiments initially rose
rapidly, and then declined to a steady state level as Mg-rich smectite precipitated.
Calcium concentrations also rose rapidly during the initial stages of these experiments,
and continued to increase until the last stages of the experiments. Thus, log [(aCa⫹2)/
(aMg⫹2)] increased with time in response to precipitation of Mg-rich smectite and the
fluids evolved toward the stability fields of zeolites as suggested in figure 14. Similar
experiments by Gı́slason and Eugster (1987a) on dissolution of basaltic glass and
crystalline basalt at 25°, 45°, and 65°C resulted in nearly identical solution composition
trends. Subsequent analysis of the solid phases produced in these experiments (Gı́slason,
Veblen, and Livi, 1993) identified several magnesium-bearing phases (kerolite, chryso-
tile, talc, and smectite) along with quartz, calcite, kaolinite, and amorphous precipitates.
Early formation of magnesian minerals in these experiments (including smectite) con-
paragenesis during low-grade metamorphism of basaltic lavas at Teigarhorn 495
Fig. 14. Stabilities of minerals in the system MgO-CaO-Na2O3-Al2O3-SiO 2-H2O-HCl at 100°C, 500 bar
in equilibrium with quartz and albite as functions of log [(aCa⫹2)/(aH⫹)2] and log [(aMg⫹2)/(aH⫹)2] in an aqueous
solution with aH2O ⫽ 1. Diagram calculated using thermodynamic properties and algorithms in the computer
code SUPCRT92 of Johnson, Oelkers, and Helgeson (1992), with addition of data for stilbite
(CaAl2Si7O18 · 7H2O; ⌬Gfo (cal/mol) ⫽ ⫺2381234; ⌬Hfo (cal/mol) ⫽ ⫺25588306; So (cal/molK) ⫽ 191.41; Vo
(cm3/mol) ⫽ 332.61; Cpo (cal/molK) ⫽ 154.12 ⫹ 0.168 T(K) ⫺ 10860000 T(K)⫺2; Neuhoff, unpublished data).
The diagram was constructed for conditions approximating the maximum burial temperature and lithostatic
pressure at Teigarhorn (see text). The right-hand axis gives values of log fco2 in equilibrium with calcite. Phase
boundaries between Al-bearing minerals are shown as solid lines. Saturation surfaces for stoichiometric
diposide, talc, wollastonite, and tremolite are depicted as dotted lines. Two sets of phase boundaries involving
clinochlore are depicted in figure 12: one for endmember clinochlore (aclinochlore ⫽ 1; solid lines) and one
corresponding to aclinochlore ⫽ 0.0068 (dash-dot-dot lines), the activity of clinochlore in sample 94-94 (table 2)
assuming ideal mixing (Holland, Baker, and Powell, 1998). Prehnite is not depicted on the diagram as the
thermodynamic properties reported for this phase in Johnson, Oeklers, and Helgeson (1992) imply that it is
stable relative to all Ca-zeolites, conflicting with experimental and field evidence indicating that Ca-zeolites are
stable relative to prehnite at this temperature and pressure. Trends A and B are discussed in the text.
sumed Mg, causing the value of log [(aCa⫹2)/(a Mg⫹2)] to increase as inferred above from
natural parageneses.
Fluid compositions in local equilibrium with Stage II mineral assemblages were
predicted using Version 8 of the computer code EQ3NR (Wolery, 1991) assuming
heterogeneous equilibrium at 100°C, 500 bar between clinochlore, stilbite, quartz, and
albite. Additional constraints were provided by fixing the concentration of Cl⫺ to 0.001
(representative concentration in meteoric-water dominated geothermal systems in Ice-
496 Philip S. Neuhoff and others—Porosity evolution and mineral
Fig. 15. Variation of bulk porosity (percent) in lavas at Teigarhorn as a function of paragenetic stage
(based on the measurements in tables 7 and 8). Time elapsed after eruption increases to the right. Vertical
brackets illustrate the ranges in total primary and secondary porosity determined by outcrop and thin section
measurements. The gray region illustrates changes in the range of bulk porosity values realized through
geologic time, with the relative height and absolute position of the boundaries of this region illustrating the
variation in and magnitude of porosity.
suggest that the percent of ⌽ occluded by mineral precipitation during Stage III was
smaller than that in Stage II as evidenced by the abundant open cavities remaining in
breccias and fractures that host museum-grade zeolite, quartz, and Iceland spar crystals.
The textures and extent of infilling of primary pore spaces during Stages I and II
suggest that although ⌽ decreased during alteration, the matrix of the lavas was
sufficiently permeable to permit transport of chemical components necessary for precipi-
tation of observed mineral assemblages. As noted above, mass transfer must have
occurred from the matrix of the lavas to primary pore spaces. This apparently was true
even after pore walls were coated with Stage I silica and celadonite and Stage II C/S.
With the exception of scoreaceous zones, where connections between pore spaces are
readily observable, much of the Stage I and II pore filling occurred in effectively isolated
vesicles. Vesicles hosting Stage II minerals were thus only partially sealed by precipita-
tion of Stage I minerals and Stage II C/S. The nearly complete infilling of ⌽1* during
Stage II was particularly acute because of the large molar volumes of zeolites and
smectites. Subsequent alteration was limited by the lack of pore space for storage and
migration of fluids, except where brittle deformation locally increased ⌽, resulting in
Stage III alteration.
Although ⌽ was larger in flow tops at the onset of Stage II than in Stage III (table 7),
fractures conducting fluid flow during Stage III are often laterally continuous in outcrop
for as much as 10 m and associated with irregular porous metasomatized aureoles (figs.
10, 11C and D). Thus, k during Stage III was probably much greater initially than in
porous zones occupied by Stage I and II minerals, despite lower values of ⌽ (table 7).
The extensive metasomatism producing the vuggy SiO2-rich zones near Stage III pores
was likely a consequence of high water/rock ratios possible in large, conductive
fractures. Late precipitation of zeolites closed the fractures, resulting in a decrease of both
⌽ and k. The duration of convective fluid and heat transport during Stage III was
probably limited by declining k due to fracture filling.
498 Philip S. Neuhoff and others—Porosity evolution and mineral
CONCLUDING REMARKS
The mineral parageneses observed in eastern Iceland and other low-grade metaba-
salt localities often result from multiple stages of alteration that can be difficult to
distinguish. However, as shown in figure 13, the textural and mineralogical signatures of
low-grade metamorphism can be interpreted by relatively simple paragenetic schemes
that account for the tectonic history of the lavas. It is also readily apparent that mineral
precipitation and tectonic events accompanying the paragenetic stages depicted in figure
13 profoundly influence the hydrogeologic properties of basaltic crust (fig. 15).
An example of the importance of petrotectonic interpretations of zeolite-facies
mineral parageneses is the ‘‘zone of abundant zeolites,’’ described by Walker (1960b) as
the deepest exposed mineral zone developed in non-olivine normative tholeiitic lavas in
eastern Iceland. The assemblage described by Walker (1960b) in this zone (‘‘magnificent
specimens’’ of heulandite, stilbite, scolecite, epistilbite) is similar to the integrated
mineralogy of Stage II and III alteration found at Teigarhorn, which Walker (1960b)
cited as a type locality for the ‘‘zone of abundant zeolites.’’ Inspection of figures 3 and 4 in
Walker (1960b) suggests that the vertical extent of the ‘‘zone of abundant zeolites’’ is
correlated to dike distribution along the northern shore of Berufjördur just as the
distribution of Stage III alteration at Teigarhorn is associated with dikes. Based on these
comparisons, it appears that the ‘‘zone of abundant zeolites’’ reported by Walker (1960b)
represents hydrothermal aureoles around dike swarms and is not a signature of burial
metamorphism.
ACKNOWLEDGMENTS
The authors would like to thank Herbert Hjörleifsson for graciously granting us
access to the outcrops at Teigarhorn, Elsa Thórey Eysteinsdóttir, Gunnlaugur Fridbjarnar-
son, and Obba, Óskar, and Berta Sandholt for friendship and hospitality, and the citizens
of Djúpivogur for accommodating us during our work at Teigarhorn. Anne Dyreborg
assisted with the fluid inclusion analyses. This work was funded by the Shell and McGee
Funds at Stanford University, the Geological Society of America, the Icelandic Student
Innovation Fund, the Petroleum Research Fund (grant ACS-PRF 31742-AC2 to DKB),
and graduate fellowships from the United States National Science Foundation and
Environmental Protection Agency (to PSN). Discussions with L. Heister, A. Stefánsson,
and S. Jakobsson and reviews by J. Liou and P. Schiffman contributed significantly to our
interpretations.
REFERENCES
Ague, J. J., 1995, Deep crustal growth of quartz, kyanite, and garnet into large-aperture, fluid-filled fractures,
north-eastern Connecticut, USA: Journal of Metamorphic Geology, v. 13, p. 299-314.
Alberti, A., Pongiluppi, D., and Vezzalini, G., 1982, The Crystal Chemistry of Natrolite, Mesolite and
Scolecite: Neues Jarbuch fur Mineralogie, Abhandlungen, v. 143, p. 231-248.
Alt, J. C., 1999, Very low-grade hydrothermal metamorphism of basic igenous rocks, in Frey, M., and
Robinson, D., editors, Low Grade Metamorphism: Oxford, Blackwell Science Ltd., p. 169-201.
Arnórsson, S., Gunnlaugsson, E., and Svavarsson, H., 1983, The chemistry of geothermal waters in Iceland; II,
Mineral equilibria and independent variables controlling water compositions: Geochemica et Cosmo-
chimica Acta, v. 47, p. 547-566.
Arnórsson, S., Sigurdsson, S., and Svavarsson, H., 1982, The chemistry of geothermal waters in Iceland; I,
Calculation of aqueous speciation from 0 degrees to 370 degrees C: Geochemica et Cosmochimica Acta,
v. 46, p. 1513-1532.
Bailey, S. W., 1988, Chlorites: structures and crystal chemistry, In Bailey, S. W., editor, Hydrous Phyllosili-
cates: Mineralogical Society of America Reviews in Mineralogy, v. 19, p. 347-403.
Bauer, J., and Hrı̂chová, R., 1962, Ptilolite (mordenite) from Teigarhorn, Iceland: Sbornik Vysoke Skoly
Chemicko-Technologicke V Praze, Mineralogie, v. 6, p. 27-36.
Bellanca, A., 1940, Piezoelettricita e legge di geminazione della scolecite: Societa de Scienze Naturali Ed
Economiche di Palermo. Bollettino, v. 22, p. 11-18.
Betz, V., 1981, Famous mineral localities; zeolites from Iceland and the Faeroes: The Mineralogical Record,
v. 12, p. 5-26.
Bevins, R. E., Robinson, D., and Rowbotham, G., 1991, Compositional variations in mafic phyllosilicates from
regional low-grade metabasites and application of the chlorite geothermometer: Journal of Metamorphic
Geology, v. 9, p. 711-721.
paragenesis during low-grade metamorphism of basaltic lavas at Teigarhorn 499
Bevins, R. E., Rowbotham, G., and Robinson, D., 1991, Zeolite to prehnite-pumpellyite facies metamorphism
of the late Proterozoic Zig-Zag Dal Basalt Formation, eastern North Greenland: Lithios, v. 27, p. 155-165.
Blake, D. H., 1970, Geology of Álftafjördur volcano, a Tertiary volcanic centre in south-eastern Iceland:
Visindafelag Islendinga [Rit], v. 2, p. 43-63.
Bolton, E. W., Lasaga, A. C., and Rye, D. M., 1999, Long-term flow/chemistry feedback in a porous medium
with heterogeneous permeability: kinetic control of dissolution and precipitation: American Journal of
Science, v. 299, p. 1-68.
Buckley, H. A., Bevan, J. C., Brown, K. M., and Johnson, R., 1978, Glauconite and celadonite: Two separate
mineral species: Mineralogical Magazine, v. 42, p. 373-382.
Cathelineau, M., 1988, Cation site occupancy in chlorites and illites as a function of temperature: Clay
Minerals, v. 23, p. 471-485.
Cathelineau, M., and Nieva, D., 1985, A chlorite solid solution geothermometer. The Los Azulfres geothermal
system (Mexico): Contributions to Mineralogy and Petrology, v. 91, p. 235-244.
Chayes, F., 1956, Petrographic model analysis: an elementary statistical appraisal: New York, John Wiley and
Sons, 113 p.
Christiansen, F. G., Boesen, A., Bojesen-Koefoed, J., Dallhoff, F., Dam, G., Neuhoff, P. S., Pedersen, A. K.,
Pedersen, G. K., Stannius, L. S., and Zinck-Joergensen, K., 1999, Petroleum geological activities onshore
West Greenland in 1997: Geology of Greenland Survey Bulletin, v. 180, p. 10-17.
Christensen, N. I., and Wilkens, R. H., 1982, Seismic properties, density, and composition of the Icelandic
crust near Reydarfjördur: Journal of Geophysical Research, v. 87, p. 6389-6395.
Dewers, T., and Ortoleva, P., 1990, Geochemical self-organization III: Mechano-chemical model of metamor-
phic differentiation: American Journal of Science, v. 290, p. 473-521.
Flóvenz, O. G., and Sæmundsson, K., 1993, Heat flow and geothermal processes in Iceland: Tectonophysics,
v. 225, p. 123-138.
Flower, M. F. J., Pritchard, R. G., Brem, G., Cann, J. R., Delaney, J., Emmerman, R., Gibson, I. L., Oakley, P. J.,
Robinson, P. T., and Schmincke, H.-U., 1982, Chemical stratigraphy, Iceland Research Drilling Project,
Reydarfjordur, eastern Iceland: Journal of Geophysical Research, v. 87, p. 6489-6610.
Galli, R., and Rinaldi, R., 1974, The Crystal Chemistry of Epistilbites: American Mineralogist, v. 59,
p. 1055-1061.
Ghiara, M. R., Franco, E., Petti, C., Stanzione, D., and Valentino, G.M., 1993, Hydrothermal interaction
between basaltic glass, deionized water and seawater: Chemical Geology, v. 104, p. 125-138.
Gı́slason, S. R., and Arnórsson, S., 1993, Dissolution of primary basaltic minerals in natural waters: saturation
state and kinetics: Chemical Geology, v. 105, p. 117-135.
Gı́slason, S. R., Arnórsson, S., and Ármannsson, H., 1996, Chemical weathering of basalt in Southwest Iceland;
effects of runoff, age of rocks and vegetative/glacial cover: American Journal of Science, v. 296,
p. 837-907.
Gı́slason, S. R., and Eugster, H. P., 1987a, Meteoric water-basalt interactions. I: A laboratory study:
Geochemica et Cosmochimica Acta, v. 51, p. 2827-2840.
————— 1987b, Meteroric water-basalt interactions. II: A field study in NE Iceland: Geochemica et Cosmo-
chimica Acta, v. 51, p. 2841-2855.
Gı́slason, S. R., Veblen, D. R., and Livi, K. J. T., 1993, Experimental meteoric water-basalt interactions:
Characterization and interpretation of alteration products: Geochimica et Cosmochimica Acta, v. 57,
p. 1459-1471.
Grant, J. A., 1982, The Isocon Diagram—A Simple Solution to Gresens’ Equation for Metasomatic Alteration:
Economic Geology, v. 81, p. 1976-1982.
Gústafsson, L. E., 1992, Geology and petrography of the Dyrfjöll central volcano, eastern Iceland: Berliner
Geowissenschaftliche Abhandlungen, v. 138, p. 1-101.
Güven, N., 1988, Smectites, In Bailey, S. W., editor, Hydrous Phyllosilicates: Mineralogical Society of America
Reviews in Mineralogy, v. 19, p. 497-559.
Heister, L. E., ms, 1997, The petrogenesis of Tertiary phonolitic volcanism Gronau West Nunatak, East
Greenland: M.S. thesis, Arizona State University, Tempe, Arizona, 12 p.
Holland, T., Baker, J., and Powell, R., 1998, Mixing properties and activity-composition relationships of
chlorites in the system MgO-FeO-Al2O3-SiO2-H2O: European Journal of Mineralogy, v. 10, p. 395-406.
Jóhannesson, H., and Sæmundsson, K., 1989, Geological Map of Iceland. 1:500000. Bedrock Geology:
Reykjavik, Iceland, Icelandic Museum of Natural History and Iceland Geodetic Survey, 1 sheet.
Johnson, J. W., Oelkers, E. H., and Helgeson, H. C., 1992, SUPCRT92: Software package for calculating the
standard molar thermodynamic properties of minerals, gases, aqueous species, and reactions among them
as functions of temperature and pressure: Computers in Geoscience, v. 18, p. 899-947.
Jørgensen, O., 1984, Zeolite zones in the basaltic lavas of the Faeroe Islands. In The deep drilling project
1980-1981 in the Faeroe Islands: Annales Societatia Scientarum Faroensis, Supplementum, v. 9, p. 71-91.
Kristjánsson, L., Gudmundsson, A., and Haraldsson, H., 1995, Stratigraphy and paleomagnetism of a
3-km-thick Miocene lava pile in the Mjoifjordur area, eastern Iceland: Geologische Rundschau, v. 84,
p. 813-830.
Kristmannsdóttir, H., 1979, Alteration of basaltic rocks by hydrothermal activity at 100-300°C, in Mortland,
M., and Farmer, V., editors, Developments in Sedimentology, v. 27: Amsterdam, Elsevier, p. 359-367.
Kristmannsdóttir, H., and Tómasson, J., 1978, Zeolite zones in geothermal areas in Iceland, in Sand, L. B., and
Mumpton, F. A., editors, Natural Zeolites: Oxford, Pergamon Press Ltd., p. 277-284.
Lasaga, A. C., 1986, Metamorphic reaction rate laws and the development of isograds: Mineralogical
Magazine, v. 50, p. 359-373.
————— 1989, Fluid flow and chemical reaction kinetics in metamorphic systems: a new simple model: Earth
and Planetary Science Letters, v. 94, p. 417-424.
500 Philip S. Neuhoff and others—Porosity evolution and mineral
Lasaga, A. C., and Rye, D. M., 1993, Fluid flow and chemical reaction kinetics in metamorphic systems:
American Journal of Science, v. 293, p. 361-404.
Liou, J. G., Seki, Y. S., Guillemette, R. N., and Sakai, H. U., 1985, Compositions and parageneses of secondary
minerals in the Onikobe geothermal system, Japan: Chemical Geology, v. 49, p. 1-20.
Manning, C. E., and Bird, D. K., 1991, Porosity evolution and fluid flow in the basalts of the Skaergaard
magma-hydrothermal system, East Greenland: American Journal of Science, v. 291, p. 201-257.
————— 1995, Porosity, permeability, and basalt metamorphism, in Schiffman, P., and Day, H. W., editors,
Low-Grade Metamorphism of Mafic Rocks: Geological Society of America Special Paper 296, p. 123-140.
Manning, C. E., and Ingebritsen, S. E., 1999, Permeability of the continental crust: Implications of geothermal
data and metamorphic systems: Reviews in Geophysics, v. 37, p. 127-150.
Manning, C. E., Ingebritsen, S. E., and Bird, D. K., 1993, Missing mineral zones in contact metamorphosed
basalts: American Journal of Science, v. 293, p. 894-938.
Mehegan, J. M., Robinson, P. T., and Delaney, J. R., 1982, Secondary mineralization and hydrothermal
alteration in the Reydarfjördur drill core, eastern Iceland: Journal of Geophysical Research, v. 87,
p. 6511-6524.
Murata, K. J., Formoso, M. L. L., and Roisenberg, A., 1987, Distribution of zeolites in lavas of southeastern
Parana Basin, state of Rio Grande do Sul, Brazil: Journal of Geology, v. 95, p. 455-467.
Neuhoff, P. S., Watt, W. S., Bird, D. K., and Pedersen, A. K., 1997, Timing and Structural Relations of Regional
Zeolite Zones in Basalts of the East Greenland Continental Margin: Geology, v. 25, p. 803-806.
Norton, D., 1988, Metasomatism and permeability: American Journal of Science, v. 288, p. 604-618.
Norton, D., and Knapp, R., 1977, Transport phenomena in hydrothermal systems: The nature of porosity:
American Journal of Science, v. 277, p. 913-936.
Odin, G. S., Desprairies, A., Fullagar, P. D., Bellon, H., Decarreau, A., Froehlich, F., and Zelvelder, M., 1988,
Nature and geological significance of celadonite, in Odin, G. S., editor, Green marine clays; oolitic
ironstone facies, verdine facies, glaucony facies and celadonite-bearing facies; a comparative study:
Developments in Sedimentology, v. 45, p. 337-398.
Ortoleva, P., Chadam, J., Merino, E., and Sen, A., 1987, Geochemical self-organization II: The reactive-
infiltration instability: American Journal of Science, v. 287, p. 1008-1040.
Oxburgh, E. R., and Agrell, S. O., 1982, Thermal conductivity and temperature structure of the Reydarfjordur
borehole: Journal of Geophysical Research, v. 87, p. 6423-6428.
Pálmason, G., 1973, Crustal rifting, and related thermo-mechanical processes in the lithosphere beneath
Iceland: Geologische Rundschau, v. 70, p. 244-260.
Passaglia, E., 1975, The crystal chemistry of mordenites: Contributions to Mineralogy and Petrology, v. 50,
p. 65-77.
Passaglia, E., Galli, E., Leoni, L., and Rossi, G., 1978, The Crystal Chemistry of Stilbites and Stellerites: Bulletin
de Mineralogie, v. 101, p. 368-375.
Ragnarsdóttir, K. V., 1993, Dissolution kinetics of heulandite at pH 2-12 and 25°C: Geochimica et Cosmo-
chimica Acta, v. 57, p. 2439-2449.
Reynolds, R. C. J., 1988, Interstratified clay minerals, in Brindley, G. W., and Brown, G., editors, Crystal
Structures of Clay Minerals and Their X-ray Identification: London, Mineralogical Society, p. 249-304.
Robert, C., and Goffé, B., 1993, Zeolitization of basalts in subaqueous freshwater settings: Field observations
and experimental study: Geochimica et Cosmochimica Acta, v. 57, p. 3597-3612.
Robinson, D., and Bevins, R. E., 1999, Patterns of regional low-grade metamorphism in metabasites, in Frey,
M., and Robinson, D., editors, Low Grade Metamorphism: Oxford, Blackwell Science Ltd., p. 143-168.
Robinson, D., Bevins, R. E., and Rowbotham, G., 1993, The characterization of mafic phyllosilicates in
low-grade metabasalts from eastern North Greenland: American Mineralogist, v. 78, p. 377-390.
Robinson, P. T., Mehegan, J., Gibson, I. L., and Schmincke, H.-U., 1982, Lithology and structure of the
volcanic sequence in eastern Iceland: Journal of Geophysical Research, v. 87, p. 6429-6436.
Rose, N. M., 1995, Geochemical consequences of fluid flow in porous basaltic crust containing permeability
contrasts: Geochemica et Cosmochimica Acta, v. 59, p. 4381-4392.
Sæmundsson, K., 1979, Outline of the geology of Iceland: Jökull, v. 29, p. 7-28.
Schiffman, P., and Fridleifsson, G. Ó., 1991, The smectite-chlorite transition in drillhole NJ-15, Nesjavellir
geothermal field, Iceland: XRD, BSE and electron microprobe investigations: Journal of Metamorphic
Geology, v. 9, p. 679-696.
Schmidt, S. Th., 1990, Alteration under conditions of burial metamorphism in the North Shore Volcanic
Group, Minnesota—Mineralogical and geochemical zonation: Heidelberger Geowissenschaftliche Abhand-
lungen, v. 41, 309 pp.
Schmidt, S. Th., and Robinson, D., 1997, Metamorphic grade and porosity and permeability controls on mafic
phyllosilicate distributions in a regional zeolite to greenschist facies transition of the North Shore Volcanic
Group, Minnesota: Geological Society of America Bulletin, v. 109, p. 683-697.
Seki, Y., Liou, J. G., Guillemette, R., Sakai, H., Oki, Y., Hirano, T., and Onuki, H., 1983, Investigation of
geothermal systems in Japan. I. Onikobe Geothermal Area: Saitama University 3, Hydroscience and
Geotechnology Laboratory, Memoir, 206 pp.
Shau, Y.-H., and Peacor, D. R., 1992, Phyllosilicates in hydrothermally altered basalts from DSDP Hole 504B,
leg 83_A TEM and AEM study: Contributions to Mineralogy and Petrology, v. 112, p. 119-133.
Slaughter, M., and Kane, W. T., 1969, The crystal structure of a disordered epistilbite: Zeitschrift für
Kristallographie, v. 130, p. 68-87.
Steefel, C. I., and Lasaga, A. C., 1994, A coupled model for transport of multiple chemical species and kinetic
precipitation/dissolution reactions with applications to reactive flow in single phase hydrothermal
systems: American Journal of Science, v. 294, p. 529-592.
paragenesis during low-grade metamorphism of basaltic lavas at Teigarhorn 501
Steefel, C. I., and Van Cappellen, P., 1990, A new kinetic approach to modeling water-rock interaction: the
role of nucleation, precursors, and Ostwald ripening: Geochimica et Cosmochimica Acta, v. 54, p. 2657-
2677.
Sukheswala, R. N., Avasia, R. K., and Gangopadhyay, M., 1974, Zeolites and associated secondary minerals in
the Deccan Traps of Western India: Mineralogical Magazine, v. 39, p. 658-671.
Tingle, T. N., Neuhoff, P. S., Ostergren, J. D., Jones, R. E., and Donovan, J. J., 1996, The effect of ‘missing’
(unanalyzed) oxygen on quantitative electron microprobe microanalysis of hydrous silicate and oxide
minerals: Geological Society of America, Abstracts with Programs, v. 28, p. A-212.
Walker, G. P. L., 1951, The amygdale minerals in the Tertiary lavas of Ireland. I. The distribution of chabazite
habits and zeolites in the Garron plateau area, County Antrim: Mineralogy Magazine, v. 29, p. 773-791.
————— 1959, Geology of the Reyoarfjörour area, eastern Iceland. The Quarterly Journal of the Geological
Society of London, v. 114, p. 367-393.
————— 1960a, The amygdale minerals of the Tertiary lavas of Ireland. III. Regional distribution: Mineralogical
Magazine, v. 32, p. 515-528.
————— 1960b, Zeolite zones and dike distribution in relation to the structure of the basalts of eastern Iceland:
Journal of Geology, v. 68, p. 515-528.
————— 1962, Tertiary welded tuffs in eastern Iceland: The Quarterly Journal of the Geological Society of
London, v. 118, p. 275-293.
————— 1963, The Breiddalur central volcano, eastern Iceland [with discussion]: The Quarterly Journal of the
Geological Society of London, v. 119, p. 29-63.
————— 1974, The structure of eastern Iceland: in Kristjánsson, l., editor, Geodynamics of Iceland and the North
Atlantic Area: Reykjavik, Iceland, NATO Advanced Study Institute, p. 177-188.
Watkins, N. D., and Walker, G. P., L., 1977, Magnetostratigraphy of eastern Iceland: American Journal of
Science, v. 277, p. 513-584.
Wolery, T. J., 1991, EQ3NR. A computer program for geochemical aqueous speciation solubility calculations:
Theoretical manual, user’s guide and related documentation (Version 7.0). Livermore, California,
Lawrence Livermore National Laboratory report UCRL-MA-110662 Pt III.
Wood, D. A., 1976, Elemental Mobility during Zeolite Facies Metamorphism of the Tertiary Basalts of Eastern
Iceland: Contributions to Mineralogy and Petrology, v. 55, p. 241-254.
————— 1977, Major and trace element variations in the Tertiary lavas of eastern Iceland and their significance
with respect to the Iceland geochemical anomaly: Journal of Petrology, v. 19, p. 393-436.