W - M - W H I T e G e o C H e M I S T R y Chapter 12: The Crust Chapter 12: Geochemistry of The Solid Earth II: The Crust

W. M.
Wh i t e G eoc h emi s t ry
Chapter 12: The Crust
Chapter 12: Geochemistry of the Solid Earth II: The

Crust
Introduction
We now turn our attention to the crust. Though the crust forms only a small fraction of the mass of
the Earth (about 0.5%), it is arguably the most varied and interesting fraction. Further, itÕs the frac-
tion we can examine directly and therefore know most about. The crust has formed through igneous
processes from the mantle over geologic time. There are two fundamental kinds of crust: oceanic and
continental. Oceanic crust, created by magmatism at mid-ocean ridges, is basaltic in composition,
thin, ephemeral, and relatively uniform. It is important in a number of respects, however. First, its
composition tells us much about the composition of the mantle from which it is derived. Second, oce-
anic crust may be, at least some times, the raw material from which continental crust is formed. Thus
we begin by briefly examining the composition of oceanic crust and the factors that control this compo-
sition. We then turn our attention to the continental crust, which is much thicker, essentially perma-
nent, and on average andesitic in composition. The continental crust is also much more varied in com-
position. Although it too has formed by magma-
tism, its evolution is far more complex than t h a t
of oceanic crust. Though we have an excellent un-
1 derstanding of how oceanic crust forms, our under-
Sediment standing of the processes that have led to t h e
2A Basaltic present continental crust is far from complete.
Pillows and Subduction-related, or Ôisland arcÕ volcanism ap-
pears to play a particularly important role in t h e
Sheet Flows formation of the continental crust, so we will pay
special attention to processes in islands arcs. W e
Sheeted Dikes will then consider the problem of interaction of
mantle-derived magma with the crust, then t h e
2B problems of the composition and evolution and
the continental crust, and finally differentiation
of the crust through melting and metamorphism.
The Oceanic Crust
Isotropic The crust beneath the oceans differs from t h e
Gabbro continental crust in a number of important re-
spects. First, it is thinner, with a typical thick-
ness of 6 km or, compared to an average thickness
3 of 35 km or so for continental crust. Second, it is
more mafic, i.e., richer in Mg and Fe and poorer in
Si that the continental crust. Third, it is tempo-
Layered rary, on average, the time between its creation a t
Gabbro mid-ocean ridges and destruction at subduction
zones is 100 million years or less, compared to an
average age of about 2 billion years for the conti-
nental crust. Finally, it is essentially monoge-
netic; the vast majority of oceanic crust is created
Figure 12.1. Schematic cross-section of the oce- at mid-ocean ridges. As a result, it is much more
anic crust. Numbers on the left indicate t h e uniform in composition than the continental crust.
designation of seismically identifiable layers. In this section, we will focus entirely on the igne-
512 January 25, 1998

W. M. Wh i t e G eoc h emi s t ry
ous part of the oceanic crust. The geochemistry of
Table 12.1: Composition of MORB ocean sediments and hydrothermal interaction be-
EPR MORB1 Primitive MORB2 tween seawater and oceanic crust are discussed in
Average Std Dev. a Chapter 15.
In a 1962 paper that he called Òan essay in
SiO 2 50.39 1.89 49.10
geopoetry,Ó Harry Hess summarized his radical
TiO 2 1.72 0.47 0.6
views on seafloor spreading. He speculated t h a t
Al 2 O 3 14.93 1.13 16.4
mid-ocean ridges were produced by rising mantle
Σ FeO 10.20 1.52 8.8
convention currents, and that these convection cur-
MnO 0.18 0.04
rents then moved laterally away from the mid-
MgO 7.34 1.30 10.3
ocean ridges, producing the phenomenon of conti-
CaO 11.29 1.38 12.4
nental drift. This concept now forms the basis of
Na 2 O 2.86 0.46 1.9
plate tectonics, the fundamental paradigm of ge-
K 2O 0.25 0.47 0.1
ology. Hess did miss one detail, however. He
P 2O 5 0.35 0.48
thought the oceanic crust was hydrated mantle,
99.52 99.60 consisting of Òserpentinized peridotite, hydrated
1
Average of 1266 analyses of basalts from t h e by release of water from the mantle over the ris-
East Pacific Rise compiled by C. Langmuir. ing limb of a [convection] current.Ó However,
2
A primitive MORB composition from B a s a l - when mantle decompresses as it rises, it does not
tic Volcanism on the Terrestrial Planets merely dehydrate, it melts. This melting gener-
ates the basalts that form the oceanic crust. In
Table 12.2. Concentrations of some respects, though, HessÕs mistake is a very minor in-
Trace Elements in Average MORB deed. Oceanic crust is very ephemeral, and for this rea-
Element ppm Element ppm son, it is sometimes better to think of it as part of t h e
K 883.7 Nd 11.18 mantle reservoir than the crustal one. Nevertheless, ig-
Sc 41.37 Sm 3.752 neous processes at mid-ocean ridges have fascinated
Co 17.07 Eu 1.335 many geochemists and much has been learned about them
Ni 149.5 Gd 5.077 in the past several decades.
Cu 74.4 Tb 0.885 Seismic studies show that the oceanic crust has a l a y -
Rb 1.262 Dy 6.304 ered structure (Figure 12.1). The uppermost layer, which
Sr 113.2 Ho 1.342 is not present at mid-ocean ridges, consists of sediments
(Seismic Layer 1). Beneath this lies Layer 2, composed
Y 35.82 Er 4.143
of basaltic lava flows and the dikes that fed their erup-
Zr 104.2 Tm 0.621
tion (the Òsheeted dike complexÓ), and finally gabbros
Nb 3.507 Yb 3.90
(Layer 3). The gabbros apparently consist both of basal-
Sn 1.382 Lu 0.589 tic magmas that crystallized in place (isotropic gabbros)
Cs 0.0141 Hf 2.974 and accumulations of minerals that crystallized as t h e
Ba 13.87 Ta 0.192 basaltic magma was held in crustal magma chambers
La 3.895 Pb 0.489 (layered gabbros). Because of the latter, the gabbros are
Ce 12.00 Th 0.1871 probably somewhat more mafic on average than are t h e
Pr 2.074 U 0.0711 basalts. Layer 2 is often divided into Layer 2A and
From Hofmann (1988).
Layer 2B, with the latter having slightly higher seis-
mic velocities. For many years it was thought t h a t
boundary between the two was the boundary between the lava flows and the dikes. Based on the re-
sults the Ocean Drilling Program, however, it appears the seismic boundary reflects instead a change
in porosity due to filling of voids and fractures by secondary minerals in Layer 2B. It is thus a meta-
morphic boundary, with the transition to the sheeted dike complex actually occurring within Layer
2B.
The average and standard deviations of the major oxides in basalts from the East Pacific Rise are
listed in Table 12.1. An average of trace element concentration in mid-ocean ridge basalt (MORB) is
listed in Table 12.2. The incompatible trace element abundances and isotope ratios are controlled
513 January 25, 1998

mainly by mantle chemistry, which we discussed in
400 SLO the preceding chapter. Hence here we will focus
< 50 mm/yr mainly on the major element composition of MORB,
300 Ave Mg# = 57.1 which is controlled mainly by igneous processes.
There was a considerable debate in the l a t e
200 Indian 1960Õs and throughout the 1970Õs whether mid-
Pacific ocean ridge basalts were ÔprimaryÕ (or close to pri-
mary), i.e., whether they were direct mantle melts
100 Atlantic or whether they had experienced extensive frac-
tionation crystallization before eruption. In t h e
1980Õs, this debate was resolved in favor of t h e
150 INTERMEDIATE view that most MORB had experienced extensive
50-80 mm/yr fractional crystallization. A critical observation is
Frequency
Ave Mg# = 52.7 that all primary mantle-derived magmas must be

100 in equilibrium with olivine of the composition
found in the mantle. Mantle olivine is typically
90-92% forsterite and according to experimental
50 studies such as that of Roedder and Emslie (1972)
(see Chapter 4), a melt would have to have an
atomic Mg/(Mg+Fe2+) ratio of around 0.72 to be in
300 FAST equilibrium with such olivine. The average compo-
250 > 80 mm/yr

200 Ave Mg# = 52.8 H AAD
4 E Tamayo
150 H
H
H H J
Kolbeinsey
100 H H H
Na2O %
50 3 Na8.0 =2.31 E EEE

E
H H HHH
H HH H
H H HHH
E E E EE EE
E E
E E
E EE E
0 10 20 30 40 50 60 70 J
EE EE
E
Mg # J JJ
J J
J JJ
E
E
E
Figure 12.2. Mg# in basalts from mid-ocean 2 JJ J

J
J JJ J
JJ
E
J
JJ
J
ridges of differing spreading rates. The Mg#
is used as an index of the extent of fractional 5 6 7 8 9 10
crystallization. Because olivine is the pri- MgO %
mary crystallizing phase, and because t h e
Figure 12.3. Na2O vs. MgO in MORB from three ar-
Mg# of olivine is much higher than that of
eas of the mid-ocean ridge system: the Australian-
the liquid, Mg# will decrease during frac-
Antarctic Discordance, the Tamayo Fracture Zone
tional crystallization. The figure shows
area of the East Pacific Rise, and the Kolbeinsey
that basalts from slow spreading ridges ex-
Ridge, just north of Iceland. Thick arrows show
perience somewhat less fractional crystal-
calculated fractional crystallization paths
lization on average than basalts from fast
(Ôliquid line of descentÕ) for each data set. Na 8.0
spreading ridges. The most fractionated ba-
values are the intersection of the fractional
salts occur on intermediate spreading rate
crystallization path with MgO concentration of
ridges where small discontinuous melt lens
8%. This is illustrated for the Tamayo data
can occur. After Sinton and Detrick (1992).
(dashed lines), which has Na 8.0 of 2.31. After
Klein and Langmuir (1987).
This ratio, multiplied by 100, is called the ÔMg numberÕ, generally written Mg#.
514 January 25, 1998

sition listed in Table 12.1 has an Mg# of 59, and therefore could not be in equilibrium with mantle
olivine. From this we can conclude that the oceanic crust, including the gabbroic section, is probably
somewhat more mafic than the average MORB listed here. The primitive MORB composition in Ta-
ble 12.1 has an Mg# of 70 and is close to a possible ÔprimaryÕ composition.
One of the remarkable features of basalts erupted at mid-ocean ridges is their uniform chemistry.
Nevertheless, the small variations in chemistry that do occur are significant in that they reveal t h e
details of the magmatic processes at mid-ocean ridges. As the discussion above suggests, fractional
crystallization is an important control on MORB composition. Olivine, augitic clinopyroxene, calcic
plagioclase, and spinel, are, with rare exceptions, the only minerals to crystallize from MORB before
eruption. Spreading rate, which is closely related to magma supply rates, seems to be an important
factor in the extent of fractional crystallization: basalts erupted on fast spreading ridges are gener-
ally more fractionated than those erupted on slow spreading ridges (Figure 12.2). On fast spreading
ridges, such as the East Pacific Rise, magma supplies rates are generally sufficient to maintain a
small steady-state magma chamber, perhaps a few hundred meters deep and 1-2 km wide, at a depth
of a few kilometers beneath
the ridge axis (Sinton and
Detrick, 1992). Melt trapped
in this chamber will cool and
3.5 F CAYMAN & WOODLARK A
B BACK ARC BASINS
crystallize. New magma ris- A F F
J ATLANTIC
ing from the mantle will mix
with the fractionated magma 3.0 Z PACIFIC FA
AA A A
in the chamber before erup- A INDIAN J B J B
AAJJJ JJ J J
tion. In contrast, magma sup- JZJZ
J A A JJB
2..5 ZZZZAJ Z Z J
Na8.0
ply rates on slow spreading ZZA

Z ZJA A ZZ Z
Z
ridges, such as the Mid-Atlan- J JJ
Z JJZZ
tic Ridge, are too low to main- J Z
JB JJ
tain a steady-state magma
2.0 J J J
Z J
J
J JJ Z
chamber. As new magma in- JJ J
jected into the crust can erupt
without mixing with older,
1.5
more fractionated magma.
The other factors that con- 12 J
trol the major element compo-
J JJ J
sition of MORB are the degree J J
JJ JJ
and depth of melting. The ef- JZ J
fects of these factors were in- 10 J Z
Z
ZZ J ZZ
Z
Z Z Z ZZ J J
vestigated by Klein and Lang- Z A A
Z
J Z
A Z
J
Fe8.0
ZZ
J
J J
muir (1987). The difficulty Z
ZJJ AJ A AJ JJ
with understanding the melt- J BJJ ZJ A A
8 J A AJ A J
ing process is that all MORB J J
J A A
AB A
have suffered some fractional FF F
crystallization, and this tends B
J B
to obscure the melting effects. 6
Klein and Langmuir found that
MORB major oxide data from 0 1000 2000 3000 4000 5000 6000
each locality on the mid-ocean Depth, m
ridge system formed a coherent
Figure 12.4. Relationship between Na8.0 and Fe8.0 in MORB and ax i al
and distinct array when oxide
depth of mid-ocean ridges. Data from regions near hot spots are
abundances or ratios of oxides
circled. Data from ridge segments near the Gal‡pagos, Azores, and
were plotted against MgO
Jan Mayen hot spots, which are enclosed in the dashed field, appear
(Figure 12.3). Calculated frac-
to have anomalously low Fe. From Klein and Langmuir (1987).
515 January 25, 1998

tional crystallization paths, illustrated in
Figure 12.3, passed through these arrays, so Maximum
g
they concluded that the trends observed in
ltin
25 Melting
Me
these plots reflected fractional crystalliza-
15 kb
um
tion. They then argued that differences be-
xim
tween the trends reflected differences in t h e
Ma
melting process. 10 kb
MgO (mole %)
Rather than attempt to solve the difficult Paths of
problem of determining the composition of
20
parental magmas, Klein and Langmuir sim-
Adiabatic
ply corrected regionally averaged data to a
5 kb Melting
common MgO concentration of 8%. They did
so simply by projecting each regional array 15 fractional
on an oxide versus MgO plot to 8% MgO. For crystallization
example, a line drawn through Na 2O data of olivine
from the Tamayo Fracture Zone region of t h e
East Pacific Rise intersects 8% MgO at 2.31% Kolbeinsey
AAD Tamayo
Na 2O (Figure 12.3). They called this t h e 10
Na 8.0 value; thus Tamayo has an Na 8.0 of
2.31. So instead of comparing primary 4 6 8 10 12
magma compositions, they compared compo- FeO (mole %)
sitions that had experienced similar extents Figure 12.5. Variation of MgO and FeO in partial
for fractional crystallization. melts of mantle peridotite. Grayed fields show t h e
Klein and Langmuir found that the re- compositions of experimental produced partial
gional corrected averages correlated strongly melts of peridotite at 3 different pressures. Hashed
with the depth of the ridge axis (Figure fields show compositions of high MgO basalts from
12.4). Comparing the variations in concen- the AAD, Tamayo, and Kolbeinsey regions. Curves
trations they observed with those in experi- for the calculated compositions at maximum and
mentally produced partial melts, they found minimum melting are also shown. Dashed arrows
that both variations in the mean pressure of show the path of melt composition produced by
melting and in the mean extent of melting melting of adiabatically rising mantle. Curved ar-
were required to explain the data. For ex- row shows how the compositional a 15 kbar melt
ample, sodium concentrations in partial will evolve due to fractional crystallization of oli-
melts appears to be controlled only by the ex- vine. After Klein and Langmuir (1987).
tent of partial melting. Iron concentrations,
on the other hand, are only weakly affected
by degree of partial melting, but are strong functions of the pressure at which melting occurs (Figure
12.5).
Klein and Langmuir concluded mantle temperature was probably the key factor in accounting for
both depth of the ridge axis and the composition of melts erupted, because mantle temperature affects
both degree of melting and the mean depth of melting. Shallow segments of the mid-ocean ridge sys-
tem overlie relatively hot mantle. The hot mantle intersects the solidus at greater depth and ulti-
mately melts to a greater degree (Figure 12.6). Hotter mantle is less dense and therefore more buoy-
ant, so that ridges overlying hotter mantle will be more elevated. Cooler mantle will not begin to
melt until it reaches shallower depth, and total extent of melting will be more limited. Klein and
Langmuir concluded that a range in degree of melting of 8-20% and in mean pressure of melting of 0.5 to
1.6 GPa were required to produced the range in compositions observed. The hottest regions of the man-
tle occur near mantle plumes such Iceland. The coolest region occurs at the Australian-Antarctic Dis-
cordance, a region where the ridge is particularly deep and isotope studies have suggested is a bound-
ary between mantle convection cells. Overall, the data suggested a range in mantle temperature of
some 250¡ C.
516 January 25, 1998

Processes in Subduction Zones and Island
a Arc Geochemistry
20%
40% 30
Me 10% M Solid Though the evolution of the continental crust has undoubt-
Me % M
ltin el us edly been complex and involved a number of processes, there is
ltin elt
g t in
g
virtual unanimous agreement among scientists that magma-
tism is primary way in which the continental crust has
formed. Today, and throughout the Phanerozoic, additions to
g
ing
the continental crust occur primarily by magmatism associated
with subduction. As we shall see, there are also good chemi-
cal arguments that subduction zone volcanism has been t h e
most important, though not necessarily exclusive, mechanism
by which the continents have formed. Volcanos of the Andes
Y X Mountains of South America, which overlie the subducting
Nazca Plate, are perhaps the best example of this process.
Not all subduction zones are located along continental mar-
gins; indeed, most are not. The Marianas are a good example
Y Ridge Axes X b of an intra-oceanic subduction zone. Some old island arcs,
however, such as Japan, parts of Indonesia, and the Aleutians
Crust have crustal seismic structures that are intermediate between
those of continental and oceanic crust. This suggests that is-
Melting 40% land arcs may eventually transform to continental crust. Plate
10% Columns motion may eventually result in these arcs accreting to conti-
30%
% Melting
nents. Thus intra-oceanic arcs may also contribute to continen-

solidus 15 kb tal growth. Thus in attempting to understand how the conti-
20% nental crust has formed, it is well worth while to consider
subduction zone processes.
10% Subduction zones are, of course, the place where oceanic
crust and lithosphere (often referred to as Òthe slabÓ) are re-
solidus 40 kb turned to the mantle. However, oceanic crust is not chemically
identical to the basalt produced at mid-ocean ridges. Two im-
Figure 12.6. a.) Cartoon of t h e portant things have happened to the crust between its crea-
pressure-temperature relationship tion and subduction. First, it has reacted with seawater at a
of adiabatically rising mantle variety of temperatures. This process, which we discuss in de-
undergoing melting. Hotter mantle tail in Chapter 15, hydrates the oceanic crust, adds some ele-
(X) intersects the solidus at higher ments to it from seawater (e.g., Mg, U), and extracts others.
pressure and ultimately melts to a Isotopic exchange with seawater affects the isotopic composi-
higher degree than cooler mantle tion of Li, B, and Sr. The second thing that happens is t h a t
(Y). The break in slope occurs the oceanic crust acquires a sedimentary veneer. Elements
because energy is consumed by gained through hydrothermal alteration and the elements in
melting (enthalpy of fusion). b.) sediments are derived almost entirely from the continental
Cartoon illustrating the re- crust. Thus the subducting oceanic crust carries with it a cer-
lationship between axial depth, tain amount of continental crustal material (particularly
crustal thickness, melting, and when sediment is subducted) and subduction zones are there-
mantle temperature. Hotter fore the principle sites of crustÑtoÑmantle mass transfer.
mantle (X) maintains the ridge a t
higher elevation because of its Major Element Composition
buoyancy. It also has a deeper melt Magmas found in island arcs (we'll use the term island arc
column and melts to a greater for all subduction zone magmatism, including continental mar-
degree, producing thicker crust gin type) appear to be predominantly andesitic. It seems un-
than cooler mantle (Y). After likely that andesite is the principle magma produced in arcs,
Klein and Langmuir (1989).
517 January 25, 1998

however. Generally, we don't see the lower
parts of arc volcanic edifices, which may be F
basaltic. There is also considerable doubt
whether an andesite can be produced by par-
tial melting of the mantle, particularly a t
depth. Most arcs sit about 100 km above t h e
3.0 0
4.
Benioff zone, and magmas may be generated
close to this depth. A safer bet is that the pri- S 2.0 5
mary magma is actually basaltic, of which an- P R 1.
desites are fractional crystallization products. I
In any case, basalt is not uncommon in intra- S
oceanic arcs. HR
In major element composition, island-arc
1.0
volcanics (IAV) are not much different from g O=
other volcanic rocks. Compared to MORB, t h e */M
principal difference is perhaps simply that si- FeO
liceous compositions are much more common
among the island-arc volcanics. Most IAV sil-
ica saturated or oversaturated; silica under- A M
saturated magmas (alkali basalts) are rare. In Figure 12.7 AFM (A=K 2O+Na 2O, F=FeO+MnO,
that sense, we might call them tholeiitic. M=MgO) diagram illustrating the difference be-
However, in the context of island arc magmas, tween tholeiitic and calc-alkaline lava series of
the term tholeiite has a more restrictive island arcs. Calc-alkaline rocks plot below t h e
meaning. Two principal magma series are rec- heavy line, tholeiites above. PRS is Kuno's pi-
ognized, one called tholeiitic, the other called geonite rock series; HRS is Kuno's hypersthene
calc-alkaline. The principle difference is that rock series. I is the Thingmuli series of Iceland.
the tholeiites differentiate initially toward After Carmichael, et al. (1973).
higher Fe and higher Fe/Mg than the calc-al-
kaline lavas. This is illustrated in Figure 4 J J
J
12.7. In addition, the tholeiites tend to be poorer in K and J
some other incompatible elements than the calc-alkaline
3 J JJ J J
Na6.0
J J J J J JJ
J JJJ
rocks. Kay et al. (1983) argued that the difference, at least for 2 JJ J
the Aleutians, relates to tectonic environment and depth of J
magma stagnation and crystallization. Tholeiites occur in ex- 1

tensional environments in the arc where magma can ascend 14
relatively rapidly and cool less. Fractional crystallization 12 J JJ J J
occurs at shallow levels. Calc-alkaline lavas ascend more J J J J
J
slowly and undergo crystallization at greater depth. 10 JJ
J
Ca6.0
J J J
J
J J
A third series is sometimes defined, the high-alumina se- 8 J
ries. IAV in general tend to have slightly higher Al contents J
6
on average than MORB or OIB, but there is very considerable 10 20 30 40 50 60 70
overlap with the OIB and MORB fields. Crustal Thickness (km)
IAV also tend to be somewhat poorer in Ti than MORB and Figure 12.8 Correlation of Ca and
6.0
OIB, though this is not necessarily a primary feature. Perfit Na with crustal thickness in is-
6.0
et al. (1980) argued that the difference in Ti content is due to land arc basalts. Ca and Na
6.0 6.0
early crystallization of Fe-Ti oxides (e.g., magnetite and i l - are the CaO and Na O concentra-
2
menite) in IAV, which buffers the Ti concentration. Thus t h e tions after correction for fractional
higher Ti concentrations in MORB may reflect the way t h e crystallization to 6.0 MgO. From
magmas differentiate rather than higher Ti concentrations of Plank and Langmuir (1988).
MORB primary melts or MORB sources. The same is true for
iron enrichment. Fe-Ti oxide precipitation may in turn depend
518 January 25, 1998

on oxygen fugacity. Oxygen fugacity is higher in IAV because of the higher water contents. IAV are
enriched in water as well as other volatiles compared to MORB and at least some OIB (e.g., H a -
waii). Island arc volcanics may also have higher CO2/H 2O ratios than MORB (Perfit et al., 1980).
Plank and Langmuir (1988) investigated the factors that control the major element composition of
island arc basalts. They treated the data in a manner analogous to Klein and Langmuir (1987), cor-
recting regional data sets to a common MgO content, but they used 6% MgO rather than the 8% used by
Klein and Langmuir. They found that Na6.0 and Ca6.0 (i.e., Na 2O and CaO concentrations corrected to
6% MgO) correlated well with crustal thickness (Figure 12.8). They argued that crustal thickness de-
termines the height of the mantle column available for melting. Most island arc volcanoes are lo-
cated above the point where the subducting lithosphere reaches a depth of 100-120 km. This suggests
that melting begins at a relatively constant depth in all island arcs. If this is so, then over which
mantle can rise and undergo decompressional melting will be less if the arc crust is thick, leading to
smaller extents of melting beneath arcs with thick crust, and higher Na 6.0 and Ca6.0 in the parental
magmas.
Trace Element Composition
The differences in trace elements between island arc volcanics and those from other tectonic envi-
ronments are probably more significant than the differences in major elements. Rare earths, however,
are not particularly distinctive. There is a very considerable range in rare earth patterns: from LRE
depleted to LRE enriched (Figure 12.9). IAV are virtually never as LRE depleted as MORB, but abso-
lute REE concentrations are, however, often low, and it is not unusual for the middle and heavy rare
earths to be present at lower concentrations than in MORB. One other aspect is of interest. Ce anoma-
lies occur in some IAV, whereas they are never seen in MORB or OIB, though they have been observed
in continental carbonatites and kimberlites. The significance of the Ce anomalies remains uncertain.
Mariana-Izu 40
20 Aleutians
20
10
40 10
New Britain 80
Sample/Chondrite
20 40
Sunda
10 20
40 10
20 L. Antilles South 20 L. Antilles North
10 10
5 5
La Pr Nd Eu Dy Er Yb Lu La Pr Nd Eu Dy Er Yb Lu
Figure 12.9. Rare earth patterns of some typical island arc volcanics. From
White and Patchett (1984).
519 January 25, 1998

Ce anomalies occur in sediment, so
11 E there is the immediate suspicion t h a t
the anomalies in IAV could be inher-
10 ited from subducted sediment.
Island-arc volcanics are richer in
9 E
the incompatible alkalis and alkaline
8 E
E earths (K, Rb, Cs, Sr, and Ba) relative
to other incompatible elements when
7 compared with MORB or OIB. This is
(Ba/La)N
6 E EE E illustrated in Figure 12.10, using t h e

EE E E Ba/La ratio. Though both IAV and
E
5 E EEEEE E
E
E E
E
E
E oceanic basalts can have a large range
E
E
E EE in rare earth patterns, as illustrated
4 EE EE E
EE EE E E
E E EEE
EE
E EE E
by the range in La/Sm ratios, t h e
E E
EE
EE
EE
E E
E EE Oceanic
3 E E EEE E E E E
E E
E EE E E
EE
E E EE E E
E EEE Basalts Ba/La ratios of IAV are generally
EEEE EE E E higher.
2 E
E EE E E
E
E
E EE E E
E
EEE E
EE
EE
E Island-arc volcanics tend to be poor
E E E E E
EE
E
E E E in the so-called high-field-strength
1 E E
E E
E
elements (HFS), i.e., those elements
with charge of +4 or +5 (Chapter 7).
1 2 3 4 5 6 While Nb and Ta are almost invaria-
(La/Sm)N bly depleted in arcs, Zr, Hf, and T h
Figure 12.10. Relative alkali-alkaline earth enrichment are not. It would be better to speak of
of IAV illustrated by plotting the (Ba/La)N ratio vs. t h e Nb and Ta depletion. This is illus-
(La/Sm)N ratio. The subscript N denotes normalization trated in Figure 12.11.
to chondritic values. From Perfit et al. (1980). Figure 12.12 is a spider diagram
summarizing the trace element differ-
ences between island arc volcanics and
MORB. To summarize, IAV
may be either light-rare-
earth enriched or light-rare- 40 B Lesser Antilles
Haw ritic
earth depleted. They gener- C SW Pacific

d
J
Sunda
Chon
ally are enriched relative to

aii
MORB in the alkalis and a l -

E
Scotia
MOR
kaline earths. Significantly, 30 G

Aleutians
d
Icelan
these are the elements most

Nb, ppm
enriched in sediments, and t h e

most soluble elements. Rela-
tive to MORB, IAV tend to be
20
depleted in Nb, and Ta, and B
sometimes in Ti, Hf, and Zr. J
JB
B B B B
J
B
These are the elements least 10 JJ BB B
B
CE
J
B
B
B
BB
B
B J
J
J
J
enriched in sediments, and are B B
BJCB BB
B BC
B B BB
BBB J JJ J JJ C C
BJ B C B B C B B G
generally highly insoluble in B BBBBE J BBB J CC C
C CB J J BCJE JJB JCJ JJ GC B B J B
J CB B
B C B
B BJ J JCBCJCBJ CCCC GCC C BJ BBB J
aqueous solutions because of JBEECBBBJG
CC
J
CCCC CC
E
G CC
CC E CCG C CJ B CC BC
G CCGCCC C C BC BC
J B
CCJ J BBB
CCJJBB G C C J
their high ratios of ionic
charge to radius (Z/r). 200 400 600 800 1000 1200
Ba, ppm
Figure 12.11. Low Nb concentrations of IAV compared to MORB
and OIB illustrated by plotting Nb vs. Ba. From Perfit et a l .
(1980).
520 January 25, 1998

100 Isotopic Composition and
Sediment Subduction
50 high-Al basalt
Island arcs overlie subduc-
Normalized Concentration
calc-alkaline basalt tion zones, which raises t h e

J J alkaline basalt obvious question of the degree
J
▲▲▲ to which subducting oceanic
J ▲ ▲ ▲ ▲▲
▲ crust and sediment might con-
10 ▲▲ MORB
J ▲ tribute to island arc magmas.
▲ J J J J J J
5 J These questions have been
J
▲ J J most successfully addressed
▲▲
Arc Tholeiite through isotope geochemistry.
J Sr isotope ratios are gener-
ally higher, and Nd isotope
1 ▲ ratios generally lower than in
Cs Rb Ba Th U N K La Ce Nd Sm Eu Gd Dy Er Yb Lu MORB, with 87Sr/86Sr ratios
Ta Sr Hf Zr Ti Y around 0.7033 and εNd of +8 be-
Figure 12.12. Spider diagram comparing typical incompatible ele- ing fairly typical of intra-
ment contents of island arc volcanics and MORB. Note the relative oceanic IAV (Hawkesworth et
enrichment in alkalis and alkaline earths and the depletion in Ta in al., 1991). Though there is
the island arc volcanics. After Sun (1980). considerable overlap with
oceanic basalts (MORB and
15.9
Sediments
15.8 les
. Antil
Sunda L
15.7 Marianas a
S. Sandwich Band
15.6 Taiwan
Pb/ Pb
ns
B
MOR
204
tia
u
Ale
ines
15.5 Philipp
207
Indian MORB
15.4
15.3
17 17.5 18 18.5 19 19.5 20
206 204
Pb/ Pb
Figure 12.13. Pb isotope ratios in island arc volcanics. Fields for the South Sandwich, Lesser An-
tilles, Aleutians, Marianas, Philippines, Taiwan, Banda and Sunda arcs are shown and compared
with fields for Atlantic and Pacific MORB (field labeled MORB) and Indian Ocean MORB
(IMORB), and modern marine sediments.
521 January 25, 1998

OIB), there is a slight tendency for IAV to have
higher Sr isotope ratios for a given Nd isotope ratio 18.95 18.9 206Pb/204Pb
L es s
and hence plot to the right of the oceanic basalt ar- 18.91
ray on a NdÐSr isotope ratio plot. This shift to 19.2 19.0
er Ant il l e s
higher Sr isotope ratios appears to result from a con- 15° 19.3
19.5
tribution of subducted oceanic crust to IAV magma 19.5
sources. This conclusion is based on the observation
19.2 19.1
19.3
that during weathering and hydrothermal altera- 19.3
tion of the oceanic crust, isotopic exchange with
19.5
seawater shifts Sr isotope ratios of the oceanic crust
to higher values. Nd isotope ratios, are virtually 10° South
unaffected by these processes because of the ex- America
Or i noc o R .
tremely low concentration of Nd in seawater.
206
Pb/204Pb isotope ratios overlap values of oceanic 65° 60° 55°
basalts, but generally having a more restricted range.
207 Figure 12.14. Pb isotope ratios in volcanics of
Pb/204Pb ratios are typically higher IAV than most
the Lesser Antilles island arc and on the At-
oceanic basalts. They tend to form steeper arrays on
207 204 206 204 lantic Plate subducting beneath it. Contours
Pb/ PbÐ Pb/ Pb plots, and overlap the field of
show the 206Pb/204Pb ratios in the sediment,
marine sediments (Figure 12.13). For most island
which increases from north to south due to
arcs, Pb isotope ratios in the arc volcanics lie be-
the increasing contribution of Orinoco River
tween sediment local to the arc and the MORB field 206 204
(Karig and Kay, 1980). We pointed out in the last sediment. The average Pb/ Pb for each
chapter that Indian Ocean MORB have lower island, which is shown beside the island,
206
Pb/204Pb and higher 207Pb/204Pb than do MORB from also increases from south to north. From
the Atlantic and Pacific. This presumably reflects a White and DuprŽ (1986).
difference in the isotopic composition between t h e
upper mantle beneath the Atlantic and Pacific on the one hand and the Indian on the other. Interest-
ingly, the Pb isotope arrays for Indian and southwestern Pacific arcs (Sunda, Banda, Philippines, and
Taiwan) are elongate toward the Indian Ocean MORB field rather than Atlantic and Pacific MORB.
On the whole then, Pb in island arc magmas appears to be a mixture Pb from local sediment and local
upper mantle.
The Lesser Antilles arc illustrates this particularly well. The Lesser Antilles arc lies on the Car-
ibbean Plate just to the north of South America. The Orinoco River, which drains the Archean
Guiana Highland, delivers sediment that contains particularly radiogenic Pb to the front of the arc.
As a result, the Pb isotope ratios in sediments of the Demerara Abyssal Plain decrease systematically
from south to north (Figure 12.14). A similar pattern of decreasing Pb isotope ratios can be seen in t h e
volcanics of the arc, apparently reflecting the changing isotopic composition of the sediment being
subducted.
Pb isotope ratios are particularly sensitive to the presence of subducted sediment because the con-
centration of Pb in sediment is so much higher, well over two orders of magnitude, than in the mantle.
In a mixture of sediment and mantle, the Pb isotope ratios of the mixture are virtually identical to
those of the sediment when as little as 2 to 3% sediment is present in the mixture. This is not true of
Sr and Nd isotope ratios.
10
Be is an even more sensitive indicator of the presence of sediment in IAV magma sources than Pb,
because it is present in young sediment, but entirely absent from the mantle. As we found in Chapter 8,
10
Be is a cosmogenic isotope; it is produced in the atmosphere by cosmic ray spallation of 14N. It has a
half-life of only 1.5 million years, so we would not expect to find significant amounts of 10Be in t h e
Earth's interior; any created before the solar system formed has long since decayed away. 10Be cre-
ated in the atmosphere is purged by rainfall and is strongly absorbed by clays of sediment and soil.
522 January 25, 1998

Workers at Carnegie Institution of
A
Japan Peru Island Arcs Washington turned to 10Be when they were
A
H
A seeking a means of determining subduction
F
HH H H A A
rates. They decided to attempt to deter-
mine the amount of 10Be in island arc vol-
canics on the assumption that sediments
Central America containing 10Be were subducted, and by
H
H H F H H
H H F H H H H H H H comparing the amount in arc volcanics with
24 the amount of sediment, they hoped to de-
A
termine how much had decayed away, and
AA
AH A A
Aleutians thereby how much time had elapsed since
AA AA AAH H H the sediment left the surface. The problem
0 5 10 15 proved to be too complex to determine sub-
10Be (106 atoms/g) duction rates quantitatively. However, in
the course of this attempt, they produced
F by far the best evidence that sediment is
F All Others indeed subducted in island arcs. Figure
12.15 compares the amount of 10Be in arc
F
lavas with that in non-arc lavas. 10Be has
H
Non-arc Volcanos Flood basalts not been found in measurable quantities (106
H
MORB atoms per gram) in non-arc lavas, but has
H
H
n = 26 been found in lavas of some arcs. The inter-
Average = 0.3 x 106 g–1 pretation is that the 10Be originates from
F F
H sediment subducted to the magma genesis
H zone. The absence of 10Be in other arc lavas
H does not mean sediment is not subducted in
0 5 10 15 those arcs. Only very young sediment, <
10Be (106 atoms/g) 10-15 Ma, will contain appreciable
amounts of 10Be, so if the uppermost sedi-
Figure 12.15. 10Be in arc and non-arc lavas. ÒAÓ indi-
ment is removed in an accretionary wedge,
cates sample from active volcano, ÒHÓ a historic erup-
no 10Be will be delivered to the magma
tion, and ÒFÓ (fresh) indicates collected during erup-
genesis zone even though deeper, older
tion. From Tera et al., 1986.
sediment is being subducted.
Figure 8.36 showed that IAV have
higher and more variable δ18O than either MORB or OIB. This may also be attributed to a contribu-
tion of subducted oceanic crust and sediment. Indeed, as we have said before, oxygen isotope ratios can
be changed only at or near the surface of the Earth, so this is the only way the δ18O of IAV can be
readily explained.
Magma Genesis in Subduction Zones
Now that we have an overview of the composition of arc magmas, letÕs consider in more detail t h e
processes that lead to the unique geochemistry of island arc magmas. Figure 12.16 summarizes petro-
genesis in the subduction zone environment.
Most geochemists and petrologists believe that arc magmas are produced primiarily within t h e
mantle wedgeà (e.g., Hawkesworth et al., 1991; Kay and Kay, 1994). The evidence for this is: ( a )
Primary arc magmas differ only slightly in their major element chemistry from oceanic basalts,
which are definitely mantle-derived. The andesitic nature of many arc magmas probably results
from fractionation crystallization in a crustal or subcrustal magma chamber. It is therefore most
àThe mantle wedge is the part of the mantle overlying the subducting slab; because the slab descends at
an angle, this region in triangular, or wedge-shaped.
523 January 25, 1998

arc crust sediment
Mantle
magma chamber
desecending lithospheric slab
magma generation
Slab dehydration
Figure 12.16. Cross-section of a subduction zone illustrating island arc magma genesis.
likely they are partial melts of peridotite rather than subducted basalt or sediment. (b) Radiogenic
isotopic and trace element systematics generally allow only small fraction of sediment (generally a
few percent or less) to be present in arc magma sources. Relatively high 3He/4He ratios in arc lavas
confirm this. (c) Rare earth patterns of island arc magmas are consistent with these magmas being
generated by partial melting of peridotite, but not by partial melting of eclogite, which would be t h e
stable form of subducted basalt at 100 km depth (the subducting lithosphere is typically located a t
about 100 to 120 km depth beneath island arc volcanos). Because the heavy rare earths partition
strongly into garnet (e.g., Figure 6.15), melts of eclogite should show steep rare earth patterns, with
low concentrations of the heavy rare earths. This is not generally the case. Rare high magnesium an-
desites, sometimes called ÒadakitesÓ (after a well documented occurrence on Adak Island in the Aleu-
tians) with steep rare earths patterns may represent exceptions to this rule and may indeed by gener-
ated by small extents of melting of subducted oceanic crust (Kay, 1978; Defant and Drummond, 1990).
It is possible that such Òslab meltsÓ were more common several billion years ago. We discuss this pos-
sibility later in the chapter. Finally, we should point out that some scientists believe arc magmas
are indeed generated in the subducting slab (e.g., Brophy and Marsh, 1986).
We have seen that isotope systematics, particularly Pb and 10Be, provide evidence of subducted
sediment in the sources of most arcs. Subducted sediment appears to influence the trace element com-
positions of arc lavas as well. Plank and Langmuir (1993) carried out careful study of the composition
of volcanics from 8 arcs and the sediments being subducted beneath them. By analyzing representa-
tive samples from the sediments and considering the proportions of sediment types being carried be-
neath the arc, they estimated the flux of elements being carried by sediment beneath the arc. They
found they could relate the degree of enrichment of most incompatible elements to the sediment flux
of that element. For example, the Ba/Na and Th/Na ratio (after correction for fractional crystalli-
zation to 6% MgO) correlate strongly with the Ba and Th sediment fluxes (Figure 12.17). Different
arcs are enriched to different degrees in these elements: for example, the Lesser Antilles arc has mod-
erate Th/Na ratios but low Ba/Na ratios. The difference appears to be due to the difference in t h e
sediment flux.
The trace element geochemistry of island arc magmas cannot, however, be explained solely by par-
tial melting of simple bulk mixtures of mantle, sediment and subducted oceanic crust. In particular,
the characteristic enrichment of the alkalis and alkaline earths cannot be accounted for in this way.
This is illustrated in Figure 12.18, which compares the observed enrichment of Lesser Antilles low-K
basalts in incompatible elements with the enrichment predicted assuming the source was a mixture of
depleted mantle and sediment. White and DuprŽ (1986) calculated the fraction of sediment in t h e
mixture from Nd isotope ratios of the basalts, and assumed simple batch melting. This simple model
524 January 25, 1998

predicted the concentrations of the rare earths and T h
140 a reasonably well, but the enrichment of Pb, Cs, Rb, U, K ,
120 G Ba, and Sr was greater than predicted. The elements
that are over-enriched are those most soluble in aqueous
Ba6.0/Na6.0
100 J solution, so the excess abundance of these ÒmobileÓ ele-

Al
80 ments suggests they are enriched due to preferential
T Mex transport in aqueous fluids. Particularly strong enrich-
60 Mar ment of these elements is characteristic of virtual all arc
40 Ant magmas, not just those of the Lesser Antilles. Thus it is
R2 = 0.924 widely believed that water released by dehydration of
20
0 200 400 600 800 the subducting oceanic crust and sediment transports mo-
Ba sediment flux bile elements from the slab to the overlying mantle
1.2 wedge. This idea is supported by the positive correlation
b .
observed between concentrations of these elements and
1.0 water contents of submarine basalts from the Marianas
J
0.8 (Stolper and Newman, 1994).
Th6.0/Na6.0
A long standing and key question is why magmas are

0.6 Al generated at all in an area where cold lithosphere is de-
scending. Decades of experimental work has shown t h a t
0.4 Ant
Mex water lowers the solidus of rock and leads to enhanced
0.2 Mar melting at any given temperature compared to dry condi-
G T R2 = 0.914 tions, hence water released by the subducting slab may
0.0 induce melting in the wedge. The best evidence of this
0 1 2 3
Th sediment flux comes from work of Stolper and Newman, 1994) on subma-
Figure 12.17. Relationship between
10
Ba/Na6.0 and Th/Na6.0 in volcanics from
8 arcs and Ba and Th sediment flux be-
neath those arcs. Horizontal bars
5
represent uncertainty in the amount of
sediment subducted; vertical bars reflect
the variance of the ratio in the arc vol-
canics. G is Guatemala, M: Mexico, J: 1
Java, T: Tonga, Mar: Marianas, Al:
Aleutians, Ant: Lesser Antilles. From 0.5
Plank and Langmuir (1993).
rine basalts from the Marianas. As mentioned

above, water concentrations correlated with
those of mobile trace elements (as well as
0.1
Sr Ba K URbCsPbTh La Ce Nd SmEu Gd Dy Er Yb Lu
those of several less mobile, but incompatible
trace elements, such as La). However, water Figure 12.18. Ratio of observed enrichment of incom-
concentrations correlated inversely with t h e patible elements in Lesser Antilles low-K basalts
concentrations of moderately incompatible compared to enrichment calculated from a sediment-
elements, such as Ti, Zr, and Na. These obser- depleted mantle mixing model based on Nd isotopic
vations at first seem contradictory. To explain composition. Enrichment of the rare earths and T h
them, Stolper and Newman (1994) developed can be accounted for by this model, but the enrich-
a model in which a mantle source consisting of ment of Pb, Cs, Rb, U, K, Ba, and Sr is greater than
a mixture of depleted mantle (i.e., a composi- predicted, probably due to preferential transport of
tion that would melt to produce average these elements in fluids. Gray area is 1 standard
MORB) and an H 2O-rich component melts to deviation of the basalts used in the average. From
varying degrees. Their calculated H 2O-rich White and DuprŽ (1986).
525 January 25, 1998

component was strongly enriched in incompatible and mobile elements, including Na. To explain t h e
inverse correlation between Na and H2O in the basalts, they found that the extent of melting must
vary inversely with the amount of H2O-rich component in the source mixture. The smallest extents of
melting (about 5%) occur in H2O-poor sources and give rise to incompatible element-poor basalts,
while the highest extents (over 20%) give rise to H2O and incompatible element-rich basalts. The
large extents of melting producing that latter leads to their being poor in elements such as Na and Ti,
which are only moderately incompatible. Thus Stolper and NewmanÕs results suggest water induces
melting in the mantle wedge.
The origin of the Nb-Ta depletion remains unclear, though there are a variety of ideas: (a) In
part, this depletion may reflect the low abundance of these elements in sediments, but it seems un-
likely this alone can account for this depletion. (b) The high charge to radius ratios of these ele-
ments makes them relatively insoluble in the aqueous solutions, and hence they may be transported to
the mantle wedge less efficiently. (c) They may be retained in specific phases, such as ilmenite, sta-
bilized by unusual conditions or composition of the slab or the web. (d) Finally, these elements may
have the same abundances in the mantle wedge as in the source of MORB, but higher degrees of melt-
ing involved in the generation of IAV compared to MORB lead to lower abundances in the former.
These possibilities are not mutually exclusive; several factors may contribute to Nb-Ta depletion.
Why island arc magmas are more silicic and appear to have experience more fractional crystalli-
zation than oceanic basalts is also unclear. Again, there are several possibilities: (a) Retardation of
plagioclase crystallization and early onset of oxide crystallization due to high concentrations of wa-
ter in IAV may drive island arc magmas to more silicic compositions than water-poor magmas such as
MORB. (b) Island arc crust is generally thicker than oceanic crust, providing a greater opportunity for
stagnation and fractional crystallization. This will certainly be true of continental margins such as
the Andes. It would be particularly difficult for basaltic magmas to ascend through such thick, low
density continental crust. (c) Subduction zones are fundamentally compressional environments, which
may make it more difficult for magmas to ascend to the surface, again causing magmas to stagnate and
crystallize. (d) Primary magmas may be more silicic and less mafic, though not necessarily andesitic,
due to generation under ÒwetÓ conditions. Again, these possibilities are not mutually exclusive.
Crust—Mantle Interaction: Assimilation
Some of the characteristics of island arc magmas that have been ascribed to sediment subduction
can also be produced when mantle-derived magmas assimilate continental crust. When mantle-de-
rived magmas encounter continental crust, the basalt heats the crust. If temperatures are high
enough, the crust may melt. The temperature reached will depend on the size of the magma body: a
relatively thin dike will cause partial melting only very locally; a larger body will cause partial
melting at greater distances (Figure 12.19). The melt produced can be assimilated by the intruding
magma. Assimilation of crust is a pervasive phenomenon and affects most, though certainly not a l l ,
magmas erupted through continental crust, both in subduction and other tectonic environments. For
this reason, we need to devote special attention to the problem.
Geochemistry of Two-Component Mixtures
In the simplest case, the resulting magma becomes a two component mixture of the crustal and man-
tle melts. If homogenization in not complete, there will be some compositional variation. If we plot
the concentration of any two elements in different samples of this mixture against each other, they
must lie on a straight line between the two end members. However, if we plot ratios of elements, or
isotope ratios, they need not lie on a straight line. Indeed, in the general case they do not; rather
they will define a curve whose equation is:
Ax + Bxy + Cy + D = 0 12.1
where x and y are the variables of the abscissa and ordinate respectively (Langmuir et al., 1978). I f
end members are designated 1 and 2 and have ratios x1 and y1, and x2 and y2 respectively, then
526 January 25, 1998

A = a 2 b 1 y 2 – a1 b 2 y 1 12.2
after after
1 yr 2 yr B = a1 b 2 – a 2 b 1 12.3
10
10 m C = a2b2x1 –a1b2x2 12.4
20 10 m
D = a 1 b 2 x 2 y 2 – a2 b 1 x 1 y 1 12.5
30 where a i is the denominator of y i and bi is t h e
B MG BM G denominator of xi. The curvature of the mixing
line will depend on the ratio r:
after after r = a1b2/a2b1 12.6
10 yr
10 100 yr
The greater the value of r is, the greater the cur-
100 m 100 m vature. Only in the special case were r=1 is t h e
20 line straight. This is illustrated in Figure 12.20.
This result is completely general and applies to
30 mixing of river water and seawater, etc. as well
B MG BMG as mixing of magmas.
after Taking the example of 143Nd/144Nd versus
2000 yr after 87
Sr/86Sr, the curvature depends on the ratio of
10 10000 yr ( Nd186Sr2)/(144Nd286Sr1). Since to a very good
144
1 km 1 km
approximation the amount of 144Nd and 86Sr are
20 proportional to total Nd and Sr respectively, r is
approximated by Nd1Sr2/Nd2Sr1. If we express
30 this ratio as r = (Nd/Sr)1/(Nd/Sr)2 we see t h a t
B MG B MG the curvature depends on the ratio of the Nd/Sr
ratio in the two end members. In mantle-derived
Figure 12.19. Temperatures as a function of time rocks Sr/Nd does not deviate greatly from 10, so
and distance surrounding basaltic dikes of mixing curves will be close to straight lines. In
various thickness. Lines B, M, and G are t h e crustal rocks and sediments, deviations from r = 1
temperatures of biotite breakdown, muscovite are more likely and curved mixing lines therefore
breakdown and granite minimum melting re-
spectively. Gray area is region of partial
melting. From Patchett (1980).
( Qa )1 100
207 204 206 204
more common. Note that on a Pb/ PbÑ Pb/ Pb plot, 10
mixing curves will always be straight lines because t h e
denominators are the same (i.e., a = b = 204Pb). 2
r=
Two component mixtures will also form straight lines
on isochron plots (e.g., 87Sr/86SrÑ87Rb/86Sr such as Figure
Q 0.5 1
7.6), because the denominators are the same. Thus mixing a 0.1
lines can be mistaken for isochrons and visa versa. One
way to distinguish the two is a ratio-element plot. A ra-
0.01
Q
tio-element plot, for example 87Sr/86Sr vs. Sr, will also be
a curve described by equation 12.1, but a ratio plotted a 2 ( )
P
against the inverse of the denominator, for example
87 86
Sr/ SrÑ1/Sr, will be a straight line. Such a plot can be b 1 P ( ) P
b 2 ( )
a useful discriminator between isochrons and mixing lines
b
because only in the latter case will 87Sr/86SrÑ1/Sr neces- Figure 12.20. Plots of ratios of elements or
sarily define a straight line (Figure 12.21). Again, this isotopes, Q/a versus P/b for mixing of end
result is completely general, and while the general prin- members 1 and 2. The numbers along t h e
curves are the values for r. From
Langmuir et al. (1978).
527 January 25, 1998

cipals have been illus-
trated with isotope ratios,
.725 ◆ A (200, 0.725) ◆ A they apply equally well to
1.0 ratios of concentrations.
.720 When a magma or series
0.8 of magmas appear to reflect
◆ ◆ 0.8 mixing, we are often faced
.715
87Sr/86Sr
with having to decide

◆ 0.6 ◆ 0.6
whether (1) two mantle-
.710 derived magmas are mix-
◆
0.4 ◆ 0.4
ing, (2) two distinct mantle
0.2 sources are mixing before
◆ 0 ◆ melting, or (3) a mantle-de-
.705 ◆ 0.2
◆ rived magma is mixing
B (450, 0.704) B with assimilated crust. In
.700 case (2), plots involving an
200 300 400 500 2 3 4 5 × 10–3 elemental concentration
Sr, ppm 1/Sr, ppm–1 will not fall on mixing lines
because partial melting and
Figure 12.21. Mixing hyperbola formed by components A and B. Af- fractional crystallization
ter Faure (1986). will change element concen-
trations. Isotope ratios
will not be changed by magma genesis so a plot of two isotope ratios will describe a mixing line in case
(2) as well as cases (1) and (3). As we pointed out in Chapter 9, stable isotopes are particularly useful
in recognizing assimilation. This is so because mantle materials have comparatively uniform stable
isotope ratios and crustal rocks very often have stable isotope ratios that are different from mantle
values.
Recognizing crustal assimilation in subduction
zones magmas can be particularly difficult be- CRUSTAL
cause many of geochemical effects characteristic CONTAMINATION
of crustal assimilation can also result from t h e 12.0 C
presence of subducted sediment component in such x=5 5 : 1
magmas. For example, the δ18OSMOW of sediment is x = 2 Sr c=
generally greater than +15 whereas the mantle 10.0 / :1
2 :1
m
value is about +5.6. Island arc magmas derived δ O x = 1

Sr
18 1 2
from a source having a component of subducted x = .5 1:
sediment will thus have slightly δ O higher
18 8.0
x = .2
1:
10
than mantle values. However, high δ18O can also

result from assimilation of crust during magma as- 6.0 SOURCE
cent. Similarly, elevated values of 87Sr/86Sr could M CONTAMINATION
result from either the subducted sediment in t h e
magma source of assimilation of continental crust. 0.703 0.705 0.707 0.709 0.711
87Sr/86Sr
By combining O isotope analyses with radio-
genic isotope analyses, it is possible to distinguish Figure 12.22. OÑSr isotope plot showing the dif-
between these two processes. Virtually all sili- ference is mixing curves produced by contaminat-
cate rocks, including sediments, magmas, and ing magma with crust (Òcrustal contaminationÓ)
peridotite, have similar O concentrations as opposed to contaminating the magma source
(generally close to 50%), whereas concentrations with subducted material (Òsource contamina-
of the radiogenic elements Sr, Nd, Hf, Pb, etc. tionÓ). x is the fraction of end member ÒCÓ (crust
vary widely. In particular, since these elements or subducted sediment) in the mixture. After
(except Os) are incompatible to some degree, their James (1981).
528 January 25, 1998

concentrations in magmas will be higher than in the mantle. Many island arc magmas have Sr con-
centrations that are greater than those typical of continental crust, but the mantle has much lower
concentrations of Sr than sediment. This means mixing curves for sediment-mantle mixing will look
very different from for magma-crust mixing.
In the case of sediment-mantle mixing, the value of r = SrM/OM/SrS/OS (where the subscripts M
and S denote the concentrations in mantle and sediment respectively) will be typically much less
than 1. In the case of magma-crust mixing, the value of r = SrM/OM/SrC/OC (where the subscripts M
and C denote the concentrations in magma and crust respectively) will typically be 1 or greater. On a
plot of δ18O vs. 87Sr/86Sr, this produces a convex curve in the case of assimilation and a concave curve in
the case of subducted sediment mixing with mantle (Figure 12.22).
Assimilation–Fractional Crystallization
As we pointed out in Chapter 9, assimilation will inevitably be accompanied by fractionation
crystallization, because heat released during crystallization provides the energy for melting sur-
rounding crust. Where fractional crystallization also occurs, assimilation will not produce simple
mixing curves. LetÕs consider this in more detail.
The change in concentration of an element, C,
100 in a magma undergoing assimilation and frac-
Ma/Mc = 0.2 tional crystallization (AFC) is given by:
R Ca
100 D = 0.1 C m/C om = F –z +
R – 1 zC om
(1 – F –z) 12.7
10
Relative Concentration in the magma (C/C0)
where Cm is the concentration in the magma, C 0m

10 is the original concentration in the magma, F is
1 the mass of magma remaining as a fraction of
0 10 the original magma mass, R is the ratio of mass
1 D=2
assimilated to mass crystallized, D is the solid-
liquid partition coefficient of the element, and
1 Ca is the concentration in the assimilated mate-
10 0.1 rial, and
z= R+D–1 12.8
0.1 D = 10 R–1
(this treatment follows DePaolo, 1981a). These
1 equations are invalid if R = 1, but as we found in
Chapter 9, R will generally be less than 1. The
variation of concentration as a function of F for R
0.01 0 = 0.2 and various values of D is shown in Figure
12.23.
0.1 The variation of an isotope ratio in assimila-
tion-fractional crystallization (AFC) is given
0.001 by:
1.0 0.8 0.6 0.4 0.2 0
F (fraction of magma remaining) R C a 1 – F –z ε + C 0 F –zε
R–1 z a m 0
εm = 12.9
Figure 12.23. Variation of concentration during R C a 1 – F –z + C 0 F –z
assimilation-fractional crystallization. Short R–1 z m
dashed red lines are for simple fractional crys- where εm is the isotope ratio in the magma, εa is
tallization. Numbers on the curves refer to v a l - the isotope ratio in the assimilated material, ε 0
ues of Ca/C om , the ratio of the concentration in is the isotope ratio in the original magma and
the assimilant to the original concentration in other variables are as defined above. We de-
the magma. After DePaolo (1981a). rived an equation for the behavior δ18O during
assimilation in Chapter 9. That equation (8.67)
529 January 25, 1998

is somewhat simpler because the concentration of
oxygen does not vary. AFC will produce curves on
18 10 0
200
plots of δ18O vs. radiogenic isotope ratios that can
ing
) 30 be quite different from the simple mixing curves in
16 60 Figure 12.23. Figure 12.24 shows some examples of
D le Mix
.2
14 = 1 1.5 such curves for δ18O vs. 87Sr/86Sr.
p
D=
(Sim
Composition of the Continental

δ18O ‰
12 5:1 Crust
D = 2:1
D= 1
9 The continental crust is extremely heterogene-

10 ous, thus task of estimating its overall composition
8 [Cm0 /Ca]Sr = 10 is a difficult one. Furthermore, only the upper

part of the continental crust is exposed to direct
sampling: the deepest borehole, drilled by t h e
6 Russians in the Kola Peninsula has reached only
.705 .710 .715 .720 .725 .730 .735 12 km and the average thickness of the continental
87Sr/86Sr crust is about 35 km. Therefore, geochemists must
rely heavily on inferences made from indirect ob-
Figure 12.24. Variation of δ18O with 87Sr/86Sr servations to estimate the composition of the con-
during AFC for a magma with an initial δ18O = tinental crust. Beginning with Clarke (1924) and
5.7 and 87Sr/86Sr = 0.703, and an assimilant with Goldschmidt (1933), a number of such estimates of
87
Sr/86Sr = 0.735 and δ18O = +19. All curves are the composition of the continental crust have been
for R = 0.2 (5:1), except for one with D = 2 for made. These have become increasingly sophisti-
which R = 0.11 (labeled 9:1). dashed lines are cated with time. Among the most widely cited
calculated Sr concentrations (ppm) assuming an works are those of Taylor and McLennan (1985,
initial concentration of Sr in the magma of 500 1995), Weaver and Tarney (1984), and Shaw et a l .
ppm. In the case where D = 1, the problem sim- (1986). Two very recent estimates are those of
plifies to one of simple mixing. From Taylor Rudnick and Fountain (1995) and Wedepohl
(1980). (1995). These estimates are not entirely independ-
ent. For example, Weaver and Tarney (1984) rely
on an early version of Taylor and McLennanÕs upper crustal estimate; Rudnick and Fountain (1995) also
rely on Taylor and McLennanÕs upper crustal estimate. Taylor and McLennan in turn rely on the work
of Shaw (1967) for many elements, as does Wedepohl (1995). In the following, we will focus particu-
larly on the estimates of Taylor and McLennan (1985, 1995), Rudnick and Fountain (1995), and Wede-
pohl (1995). In doing so, we want both to learn of the composition of the crust and to understand how
these estimates are made.
We can divide the problem of estimating crustal composition into two parts. The first is to esti-
mate the composition of the upper, accessible parts of the crust. This is referred to as the Òupper
crust.Ó Direct observations provide the most important constraints on the composition of this part of
the crust. The second problem is the composition of the deeper, less accessible part of the crust. For
this part of the crust, indirect observations, particularly geophysical ones such as seismic velocity
Table 12.3. Abundance of Igneous and Sedimentary Rocks in the Upper

Crust
Plutonic Rocks Volume % Sedimentary Rocks Volume %
Granite, Granodiorite 77 Shales 72
Quartz Diorite 8 Carbonates 15
Diorite 1 Sandstones 11
Gabbros 13 Evaporites 2
Syenites, anorthosites, peridotites 1
From Taylor and McLennan (1985).
530 January 25, 1998

and heat flow, provide key constraints on
composition. As we shall see, these observa-
SiO2 X
tions indicate that the continental crust is 10.0
Upper Continental Crust 'Shaw et al.'

X
compositionally stratified, with the lower CaOX Al2O3
K2O XX
X
MgO XNaFeO
part being distinctly more mafic (i.e., richer
in Mg and Fe and poorer in SiO 2 and incom- 1.0 2O
%
patible elements). Some workers divide t h e X
TiO2
deep crust into a ÒmiddleÓ and Òlower crustÓ,
while others consider only a single entity 1000 Ba X X P2O5
X MnO
that they refer to as the Òlower crustÓ.
Zr X X Sr
100
ppm
The Upper Crust
V XXZn
Historically, three approaches to estim- X Cr
NiX
X
X
X Y La
ating the composition of the upper continen-
tal crust have been used. The first is to esti-
10 Th X X
XPb
Cu
mate the volume of various rock types and
U Sc
X
then use typical or average compositions of
each to derive a compositional estimate. Ta- 10 100 1000 1.0 10.0
ble 12.3 gives such an estimate of the rela- ppm %
tive volumes of various igneous and sedimen- Upper Continental Crust 'Eade & Fahrig'
tary rocks in the upper crust. Continental Figure 12.25. Comparison of estimates of the average
sediments constitute about 8% of the mass of composition of the upper crust by Eade and Fahrig
the crust; if pelagic sediments are added t h e (1971) and Shaw (1967). Both were based on studies
total sediment mass is about 11%. of the Canadian shield, but different areas. From
A second approach is to average analyses Taylor and McLennan (1985).
of samples taken over a large area. An alter-
native is to mix sample powders to form com-
posites of various rock types and thus reduce
Al 2O3
10 SiO2
the number of analyses to be made (e.g., Na2O FeO
Shaw, 1967; Eade and Fahrig, 1971; Shaw et CaO
Upper Continental Crust
al., 1986; Wedepohl, 1995). Studies of these

1.0 MgO
kinds consistently produce an average upper
crustal composition similar to that of grano- TiO 2
diorite. This is encouraging since grano-
diorite is the most common igneous rock in the
1000 Mn Ba
crust. Such estimates also tend to produce Sr
Zr
V ZnRb
relatively similar average concentrations for 100
minor and trace elements, as can be seen in
Cu Cr
Figure 12.25. Ga PbNi
A third method is to let the Earth make 10 Co Th
the composites for us. A couple of kinds of Sn UCs Hf
such materials are available. Goldschmidt W
(1933) suggested the use of glacial clays in
10 100 1000 1.0 % 10
melt-water lakes adjacent to the Pleistocene ppm
ice front. An alternative but similar ap- Loess
proach is to use loess, which is fine-grained Figure 12.26. Comparison of elemental concentra-
aeolian material of Pleistocene age. The tions in loess with estimated upper crust. From Tay-
most readily available of this kind of natu- lor and McLennan (1985).
ral composite, however, is simply sediments.
531 January 25, 1998

One advantage of sedi-
ments over glacial mate-
rial is that whereas most
glacial deposits are of B • PAAS
100 J▲ JB ▲ NASC
Pleistocene age (but there ▲ B
J
are various deposits of ▲ B B ES
Sample/Chondrites
J
▲
ages ranging up to 2.3 Ga),
sediments of all ages are B
▲
J
available so that secular
B
variations in crustal com- B ▲ B
▲
▲ J ▲ B B
▲ B
J J ▲ ▲ B
▲
position can be deter-
mined.
10 J J J J J
The advantages of us-

ing geologic composites
should be obvious, but
there are disadvantages
as well. The primary
La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb
problem is that chemical Figure 12.27. Rare earth patterns of Post-Archean Australian S h a l e
fractionations are in- (PAAS) composite, the North American Shale Composite (NASC) and
volved in producing sedi- the European Shale (EC) composite. From Taylor and McLennan (1985).
ments from their parents.
Weathering of rock typically produces three fractions: sands consisting of resistant minerals, clays,
and a solution. These products are transported with varying degrees of efficiency away from the site
of production. Since elements tend to be concentrated in one of these three fractions, none of the frac-
tions will have a composition representative of the parent rock. Because it is produced primarily by
physical, rather than chemical, action, glacial loess is less susceptible to this kind of chemical frac-
tionation, though some fractionation nevertheless occurs. This is illustrated in Figure 12.26. Loess is
enriched in SiO2, Hf, and Zr relative to estimated concentrations in the upper crust. This reflects t h e
enrichment of loess in mechanically and chemically stable minerals, such as quartz and zircon, be-
cause the lighter clays are carried further from their site of origin by wind and water. Loess is also
depleted in Na and Ca, reflecting loss by leaching.
Numerous studies have shown that when rock weathers to produce a sediment, the rare earth pat-
tern of the parent is usually preserved in the sediment. This is because all the rare earths are concen-
trated in the clay fraction, which ultimately form shales. Other Group 3 elements (Sc and Y), as well
as Th, behave similarly to the rare earths during weathering. Furthermore, rare earth patterns are
remarkably similar in different shales, suggesting shales are indeed good averages of crustal compo-
sition. This is illustrated in Figure 12.27, which compares three shale composites from three conti-
nents. Because of these properties of the rare earths, S. R. Taylor and colleagues at the Australian
National University used them as a point of departure for estimating the composition of the upper
continental crust.
Though rare earth patterns in shales are representative of their sources, their absolute concentra-
tions are not. Because some elements are lost to the sand fraction and others to solution (and ultimate-
ly to chemical sediments such as evaporites and carbonates), the REE are enriched in shale by about
20%. Thus Taylor and McLennan (1985, 1995) reduced REE (and Sc, Y, and Th) concentrations by 20%,
in making their estimate of upper crustal composition. Elements other than the REE, Sc, Th, and Y
were estimated either from known ratios of elements to one of these, or were borrowed from the esti-
mates based on analysis of large numbers of samples. For example, Taylor and McLennan (1985) ob-
tained the K concentration from the Th concentrations by assuming a Th/U of 3.8 and a K/U ratio of
104. Their Sr concentration, on the other hand, is an average of Eade and Fahrig (1968) and Shaw
(1967).
532 January 25, 1998

Table 12.4 lists the upper crustal compositions estimated by Taylor and McLennan (1985, 1995) and
Wedepohl (1995). The estimates are broadly similar, and they agree on the concentrations of many
elements within 10%. Both indicate an upper crust of ÒgranodioriticÓ or ÒtonaliticÓ composition. For
some elements, however, the good agreement simply reflects the reliance of both papers on the previ-
ous work of Shaw (1967), but this is not true of all elements. For example, the two estimates agree
closely on the concentrations of the light rare earths despite having used completely different ap-
proaches in estimating these concentrations. For a few elements, the agreement is poor. For example,
Taylor and McLennanÕs estimate of the copper concentration is 75% greater than that of WedepohlÕs.
Plank and Langmuir (in press) suggested several revisions to Taylor and McLennanÕs values for Cs,
Nb, Ta, and TiO2 based on their study of marine sediments. Their values are shown in parentheses in
Table 12.4. Notice that these are substantial changes: Cs increases by over 50%, while Nb and Ta de-
Table 12.4. Composition of the Upper Continental Crust
Major Oxides ( wt %) T& M Wedepohl Normative Mineralogy (T & M)
SiO 2 66.0 64.9 Quartz 15.7
TiO 2 0.5 (0.76) 0.52 Orthoclase 20.1
Al 2O 3 15.2 14.6 Albite 13.6
FeO 4.5 3.97 Diopside 6.1
MgO 2.2 2.24 Hypersthene 9.9
CaO 4.2 4.12 Il 0.95
Na 2O 3.9 3.46
K 2O 3.4 4.04
T & M Wedepohl T & M Wedepohl T&M Wedepohl
Li ppm 20 22 Ga ppm 17 14 Nd ppm 26 25.9
Be ppm 3 3.1 Ge ppm 1.6 1.4 Sm ppm 4.5 4.7
B ppm 15 17 As ppm 1.5 2 Eu ppm 0.88 0.95
C ppm 3240 Se ppm 0.05 0.083 Gd ppm 3.8 2.8
N ppm 83 Br ppm 1.6 Tb ppm 0.64 0.5
F ppm 611 Rb ppm 112 110 Dy ppm 3.5 2.9
Na % 2.89 2.57 Sr ppm 350 316 Ho ppm 0.8 0.62
Mg % 1.33 1.35 Y ppm 22 20.7 Er ppm 2.3
Al % 8.04 7.74 Zr ppm 190 237 Tm ppm 0.33
Si % 30.8 30.35 Nb ppm 25 (13.7) 26 Yb ppm 2.2 1.5
P ppm 700 665 Mo ppm 1.5 1.4 Lu ppm 0.32 0.27
S 953 Pd ppb 0.5 Hf ppm 5.8 5.8
Cl 640 Ag ppb 50 55 Ta ppm 2.2 (0.96) 1.5
K% 2.8 2.87 Cd ppb 98 102 W ppm 2 1.4
Ca % 3 2.94 In ppb 50 61 Re ppb 0.4
Sc ppm 11 7 Sn ppm 5.5 2.5 Os ppb 0.05
Ti ppm 3000 (4560) 3117 Sb ppm 0.2 0.31 Ir ppb 0.02
V ppm 60 53 Te ppb Au ppb 1.8
Cr ppm 35 35 I ppm 1.4 Hg ppb 56
Mn ppm 600 527 Cs ppm 3.7 (7.3) 5.8 Tl ppb 750 750
Fe % 3.5 3.09 Ba ppm 550 668 Pb ppm 20 17
Co ppm 10 11.6 La ppm 30 32.3 Bi ppm 127 123
Ni ppm 20 18.6 Ce ppm 64 65.7 Th ppm 10.7 10.3
Cu ppm 25 14.3 Pr ppm 7.1 6.3 U ppm 2.8 2.5
Zn ppm 71 52
T&M: Taylor and McLennan (1985,1995), Wedepohl: Wedepohl (1995). Values shown in parentheses are revisions of
Plank and Langmuir (in press) to Taylor and McLennanÕs estimates.
533 January 25, 1998

crease by over 50%.Ê Taylor and McLennanÕs values for these elements are poorly constrained so Plank
and LangmuirÕs values should be considered superior. WedepohlÕs concentrations for these elements
are intermediate between the original Taylor and McLennan values and Plank and LangmuirÕs revised
values. The values in Table 12.4 will undoubtedly be further refined in the future.
The Middle and Lower Crust
Rocks metamorphosed deep within the crust are called granulites; an important characteristic of
such rocks is that they are anhydrous, with pyroxene replacing amphibole and biotite. They are
sometimes exposed at the surface by tectonic processes and hence can provide insights into the nature
of the lower crust. However, these granulite terranes have often been subjected to retrograde meta-
morphism (metamorphism occurring while temperatures and pressure decrease), which compromises
their value. Furthermore, questions have been raised as to how typical they are of lower continental
crust. These questions arise because granulite terranes are generally significantly less mafic than
xenoliths from the lower crust. Xenoliths perhaps provide a better direct sample of the lower crust,
but they are rare. The point is, any estimate of the composition of the middle and lower crust will
have to depend heavily on indirect inference and geophysical constraints. There are two principal
geophysical constraints:
• Heat flow in the continental crust. A portion of the heat flowing out the crust is produced by
radioactive decay of K, U, and Th within the crust (other radioactive elements do not contrib-
ute significantly to heat generation because of their long half-lives are low abundances). The
concentrations of these elements can be related to rock type, as indicated in Table 12.2.
• Seismic velocities in the continental crust. Seismic velocities depend on density, compressibil-
ity and the shear modulus, which can in turn be related to composition.
Both tell us something of first order importance about the nature of the continental crust: it is verti-
cally zoned, becoming more mafic (i.e., richer in Fe and Mg and poorer in Si and incompatible ele-
ments) with depth. LetÕs consider them in greater detail.
Heat is transported conductively through the lithosphere. The equation governing heat conduc-
tion is identical to that governing diffusive chemical transport (FickÕs First Law, equ. 5.90):
∂T 12.10
Q(z) = – k
∂z
where Q is the heat flow at some depth, z, k is the thermal conductivity, and ¶T/¶z is the change in
temperature with depth. The thermal conductivity of rocks can be measured in the laboratory and
averages about 2 WmÐ1KÐ1. Thus by measuring the thermal gradient, in boreholes and mines for exam-
ple, the heat flow can be calculated. The average heat flow of the continents is about 60 mW/m2.
This heat has two components: heat conducted out of the mantle, which is about 20 mW/m2, and h e a t
generated by radioactive decay within the continents. The concentrations of K, U, and Th observed a t
the surface of the crust would produce more heat than is observed to be leaving the continental crust i f
these concentrations were uniform through the crust. Thus, the concentrations of these elements must
Table 12.5. U, Th, and K Concentrations and Heat Production in Various Rock
Types.
Igneous Rock Type U Th K Th/U K/U Density Heat Production
(ppm) (ppm) (%) g/cm3 10-6 W/m-3
Granite/Rhyolite 3.9 16.0 3.6 4.1 0.9× 104 2.67 2.5

Granodiorite/Dacite 2.3 9.0 2.6 3.9 1.1× 104 2.72 1.5
Diorite/Andesite 1.7 7.0 1.1 4.1 0.7× 104 2.82 1.1
Gabbro/Basalt 0.5 1.6 0.4 3.2 0.8× 104 2.98 0.3
Peridotite 0.02 0.06 0.006 3.0 0.3× 104 3.28 0.01
Continental Crust 1.25 4.8 1.25 3.8 1.0× 104 Ñ 0.8
534 January 25, 1998

decrease with depth. The concentrations of these elements are related to rock type, as is illustrated
in Table 12.5. The problem is complicated, however, by variations in the ÒmantleÓ heat flow. Heat
flow varies significantly with tectonic age, as is illustrated in Figure 12.28. If, as we believe, the con-
tinental crust is created by magmatism, it will be
initially hot and then cool over time. Subsequent
C
episodes of magmatism may also heat the crust. In 100 M
addition, the variation seen in Figure 12.28 may
represent different thickness of the lithosphere LPa
80
Heat Flow (mW/m 2)

(Vitorello and Pollack, 1990; Nyblade and Pollack, EPa
1993). The lithosphere is a conductive boundary B
60
B LPr
layer, so that the thicker the lithosphere, t h e B EPr A
lower the mantle heat flow out the top of it. Ny- B B B
blade and Pollack (1993) have argued that regions 40 I B
of old Archean crust are underlain by particularly II
thick mantle lithosphere, an argument supported
by geochronological and thermobarometric studies 20 III
of mantle xenoliths from these regions.
Seismic velocities increase with depth in t h e
crust: P-wave velocities increase from about 6 0 0.5 1.0 1.5 2.0 2.5
km/sec in the upper crust to about 7 km/sec in t h e Tectonic Age (Ga)
lower crust. The dependence of the velocity of P Figure 12.28. Heat flow as a function of tectonic
waves (compressional-waves) on physical proper- age. Component I is radiogenic heat produced in
ties is: the crust, II is heat from a transient thermal
K + 43 µ perturbation associated with tectonism, and III
vP = ρ 12.11 is background heat from deeper sources. A is Ar-
chean, EPr is Early Proterozoic, LPr is Late Pro-
and that for velocity of shear waves (S-waves) is:
terozoic, EPa is early Paleozoic, LPa is late Pa-
µ 12.12
vS = ρ leozoic, M is Mesozoic, and C is Cenozoic. From
Vitorello and Pollack (1980).
Paleozoic Paleozoic
Shields & Orogens: Orogens: Mes./Cen. Rifted Mes./Cen. Active
Platforms European other Extensional Arcs ForeArcs Margins Contractional Rifts
10 6.3±0.2 10
6.5±0.2 6.5±0.1 6.5±0.4 6.5±0.3 6.7±0.2
20 6.6±0.1 6.4±0.1 6.6±0.4 20
7.1±0.1 7.2±0.3 6.5±0.1 7.0±0.4
6.9±0.2 6.9±0.1
30 7.0±0.2 7.1±0.1 n =14 n = 30 n=8 30
7.1±0.1 n = 10 n=8
40 n = 19 6.9±0.3 40
km n = 30 n=5 km
n = 10
Key to Velocities
≤ 6.2 km/sec 6.2–6.5 km/sec 6.5–6.9 km/sec ≥ 6.9 km/sec
Figure 12.29. Seismic velocity structure of the continental crust, illustrating its 3-layered
nature. Velocity structure falls into 9 types. The number of profiles used to construct each
type is shown below each type. From Rudnick and Fountain (1995).
535 January 25, 1998

where K is the bulk modulus (inverse of β, com- 8.0
pressibility), µ is the shear modulus, and ρ is den- Measured Vp at 600 MPa
sity. Thus the increase in seismic velocity with
7.5 XXX X
depth implies increasing density and decreasing X
XXX
X XX
Vp (km/sec)
compressibility with depth. Some, but not all, of X
XX X
XX XX
the increase in density results merely from com- 7.0 X
XX
X
pression. The remainder represents changing X X
X X
composition. The real problem is to understand XX X
precisely how seismic velocity depends on compo- 6.5 X XXXXX X
XXXXX
X XX
X
sition. X
Seismic velocity profiles vary widely from 6.0
place to place, as does crustal thickness. Rudnick 40 50 60 70 80
and Fountain (1995) examined a global database SiO2 (wt. %)
of seismic cross sections and found that they can be
divided into 9 classes, which are illustrated in Figure 12.30. Correlation between measured
Figure 12.29. One must next relate seismic veloc- seismic velocity (v P ) and SiO 2 concentration.
ity to composition. This can be done by making From Rudnick and Fountain (1995).
measurements of seismic velocity in the laboratory on samples of known composition. For example,
Figure 12.30 shows the relationship between SiO2 and seismic velocity in a variety of rock types.
To produce an estimate of crustal composition, Rudnick and Fountain assigned an average lithology
to the seismic sections shown in Figure 12.29. They then assigned a composition to each lithology us-
ing a database of the composition of lower crustal xenoliths. Then by estimating the aerial extent of
each type of crustal section, they produced the compositional estimate in Table 12.6. This table
shows that the composition of the lower crust corresponds to that of tholeiitic basalt; in metamorphic
terminology it would be a mafic granulite. The composition of the middle crust corresponds to that of
an andesite. At the prevailing pressures and temperatures this rock would be an amphibolite, con-
sisting mainly of amphibole and plagioclase.
Wedepohl (1995) used the European Geotraverse as a model of the seismic structure of the crust.
This seismic cross-section runs from northern Scandinavia to Tunisia and crosses a great variety of tec-
tonic provinces, ranging from the Archean Fennoscandian Shield to the young fold belts to the young
Alpine orogen. He assigned 3 lithologies to 3 ranges of seismic velocities: sediments, granites, and
gneisses (VP < 6.5 km/s) corresponding to the upper crust, felsic granulites (6.5 < V P < 6.9 km/s), and
mafic granulites (6.9 < VP < 7.5 km/s). He used a database of compositions of felsic and mafic granu-
lites from both xenoliths and exposed terranes to calculate an average composition for each of the lat-
ter two. He then computed a lower crustal composition by weighting felsic and mafic granulites in the
proportions their characteristic seismic velocities were observed in the European Geotraverse. His
estimate of the composition of the lower crustal is also listed in Table 12.6.
Rare earth patterns of upper, middle and lower crust as estimated by Rudnick and Fountain (1995)
are compared in Figure 12.31a. The negative Eu anomaly in the upper crust and slight positive
anomalies in the middle and lower crust (such positive anomalies are typical of many granulites) are
an interesting features of these patterns. Eu is strongly held in plagioclase (Chapter 7). The presence
of plagioclase in the melting residue would produce a negative Eu anomaly in the melt and a positive
one in the residue. Thus these anomalies suggest that crustal has differentiated to form distinct l a y -
ers at least partially through partial melting, with granitic melts forming the upper crust and granu-
litic residues forming the deeper crust. Figure 12.31b displays the estimated composition of the mid-
dle and lower crust relative to the upper crust. The lower and middle crust are depleted in incompati-
ble elements and enriched in compatible elements relative to the upper crust. This is also consistent
with the idea that magmatic processes have been important in creating the compositional layering
observed in the crust.
536 January 25, 1998

Table 12.6. Composition of the Middle and Lower Continental Crust
Major Oxides, %
R&F R & F Wedepohl
Lower Middle Lower
SiO 2 52.3 60.6 58.05
TiO 2 0.8 0.7 0.84
Al 2O 3 16.6 15.5 15.52
FeO 8.4 6.4 7.34
MnO 0.1 0.1 0.12
MgO 7.1 3.4 5.23
CaO 9.4 5.1 6.80
Na 2O 2.6 3.2 2.86
K 2O 0.6 2.01 1.85
P2O5 0.1 0.1 0.20
Trace Elements, ppm
R&F R & F Wedepohl R&F R & F Wedepohl
Lower Middle Lower Lower Middle Lower
Li 6 7 13 Sn 2.1
Be 1.7 Sb 0.3
B 5 I 0.14
C 588 Cs 0.3 2.4 0.8
N 34 Ba 259 402 568
F 429 La 8 17 26.8
S 408 Ce 20 45 53.1
Cl 278 Pr 2.6 5.8 7.4
Sc 31 22 25.3 Nd 11 24 28.1
V 196 118 149 Sm 2.8 4.4 6.0
Cr 215 83 228 Eu 1.1 1.5 1.6
Co 38 25 38 Gd 3.1 4 5.4
Ni 88 33 99 Tb 0.48 0.58 0.81
Cu 26 20 37.4 Dy 3.1 3.8 4.7
Zn 78 70 79 Ho 0.68 0.82 0.99
Ga 13 17 5.4 Er 1.9 2.3
As 1.3 Tm 0.81
Se 0.17 Yb 1.5 2.3 2.5
Br 0.28 Lu 0.25 4.1 0.43
Rb 11 62 41 Hf 1.9 4 4.0
Sr 348 281 352 Ta 0.6 0.6 0.84
Y 16 22 27.2 W 0.6
Zr 68 125 165 Tl 0.26
Nb 5 8 11.3 Hg 0.21
Mo 0.6 Pb 4.2 15.3 12.5
Ag 0.08 Bi 0.037
Cd 0.101 Th 1.2 6.1 6.6
In 0.052 U 0.2 1.6 0.93
R & F: Rudnick and Fountain (1995), Wedepohl: Wedepohl (1995).
537 January 25, 1998

a The Total Continental Crust
100 H
H
H Upper Crust
The approach used by most workers to
B B H
B H B Middle Crust
B estimate the composition of the total
Sample/Chondrite
H B
B J J J Lower Crust continental crust is simply to calculate an
J J
J J H HB average of two or three crustal sections,
H J B H B B B B
10 J J H J J
H B
H H weighting each by its mass. This was
J J
done, for example, by Rudnick and Foun-
tain (1995) and Weaver and Tarney
(1984), both of whom divided the crust
into an upper, lower, and middle section.
1 Both Weaver and Tarney (1984) and
La Ce Pr Nd Sm Eu Gd TbDy Ho Er Tm Yb Lu Rudnick and Fountain (1995) relied on
10
b Mg E Cr
E Taylor and McLennanÕs upper crustal es-
Fe
Ratio to Upper Crust
TiE E E EJ JNiJ J
E E timate. Weaver and Tarney (1984) used
Y Si Al
Ba JNa J J J J E EJ E J Sc V
J J J average Lewisian* amphibolite as their
1 J J J J J
E E E E J Eu
J
Ca middle crust composition and average Le-
E E
J
J J J J
Nb J Ce E E Sm Yb Ga Sr wisian granulite as their lower crust
KEPb
E E La Zr
J E E
Th composition. Shaw et al. (1986) and
0.1 Cs uERbE E J Lower/Upper Wedepohl (1995) used a similar ap-
E E Middle/Upper proach, but divided the crust only into
upper and lower parts. An important
step in this approach is estimating t h e
Figure 12.31. (a). Comparison of chondrite-normalized thickness of the various sections. All use
in upper, middle and lower crust. (b). Elemental geophysical constraints for this.
enrichment or depletion of the middle and lower crust Taylor and McLennan (1985, 1995)
relative to the upper crust. From Rudnick and Fountain used an entirely different approach to es-
(1995). timating total crustal composition, one
based on the Òandesite modelÓ of Taylor (1967). Taylor (1967) noted the role played by subduction-re-
lated volcanism in creation of the continental crust and assumed that on average the crust consisted of
island arc andesite. Thus average island arc andesite was used as the estimated composition of t h e
continental crust. This approach was modified in subsequent work, as Taylor concluded that while
post-Archean crust was created at subduction zones, Archean crust was not and is compositionally dif-
ferent. Taylor and McLennan (1985) essentially modify the Taylor (1967) andesite model for their es-
timate of Archean crustal composition.
Estimates of the major element composition of the continental crust by Weaver and Tarney (1984),
Shaw et al. (1986), Taylor and McLennan (1995), Wedepohl (1995), and Rudnick and Fountain (1995)
are given in Table 12.7. Also listed are estimates of trace element concentrations by Taylor and
McLennan (1995), Rudnick and Fountain (1995), and Wedepohl (1995). Since Rudnick and Fountain
(1995) rely on the Taylor and McLennan upper crustal estimates in calculating the total crust composi-
tion, the revisions to the Taylor and McLennan upper crustal values for TiO 2, Cs, Nb, and Ta of Plank
and Langmuir (in press) affect the Rudnick and Fountain total crustal values of these elements. These
revisions have been made in the Table.
The ranges of estimates for SiO 2 and Al2O3 in Table 12.7 are about 10% and 8% respectively; t h e
range in Mg# (52 to 57) is similarly only about 10%. Interestingly, earlier estimates of crustal SiO 2
and Al2O3, going back to Goldschmidt (1933) also fall within this range. Thus we can conclude with
some confidence that the continental crust on the whole is similar to that of siliceous andesite.
The details of the composition of the crust are less certain, however. Ranges for the other oxides
are substantially larger: 75% for FeO, 68% for MgO, and 100% for MnO. Of these estimates, the com-
* The Lewisian, which outcrops in Northwest Scotland, is perhaps the classic exposure of lower crust.
538 January 25, 1998

position of Taylor and McLen- Table 12.7. Composition of the Continental Crust
nan is the most mafic, and that Major Oxides, wt. %
of Weaver and Tarney t h e R&F T&M W & T We Shaw
least mafic (ranges for FeO SiO 2 59.1 57.3 63.2 61.5 63.2
and MnO decrease to 30% and TiO 2 0.75 0.9 0.6 0.68 0.7
21% respectively if the esti- Al 2O 3 15.8 15.9 16.1 15.1 14.8
mates of Taylor and McLennan FeO 6.6 9.1 4.9 5.67 5.60
are excluded). MnO 0.1 0.18 0.08 0.10 0.09
The Th/U ratio in both MgO 4.4 5.3 2.8 3.7 3.15
Taylor and McLennanÕs and CaO 6.4 7.4 4.7 5.5 4.66
Rudnick and FountainÕs esti- Na 2O 3.2 3.1 4.2 3.2 3.29
mate (3.85 and 3.94 respec- K 2O 1.88 1.1 2.1 2.4 2.34
tively) is probably too low. P2O5 0.2 0.19 0.18 0.14
The Th/U ratio of the de- Trace Elements (in ppm unless otherwise noted)
pleted mantle is well con- R&F T&M We R&F T&M We
strained by 230Th/232Th ratios Li 11 13 18 Sb 0.2 0.3
in MORB to be <2.5. If t h e Be 1.5 2.4 Te, ppb 5
bulk Earth value is 3.9Ð4.2 (as B 10 11 I, ppb 800
would be the case if the terres- C 1990 Cs 3 1 3.4
trial Th/U is chondritic), then N 60 Ba 390 250 584
mass balance requires the crus- F 525 La 18 16 30
tal ratio to be higher than t h e S 697 Ce 42 33 60
bulk Earth value, yet both Cl 472 Pr 5 3.9 6.7
Taylor and McLennanÕs and Sc 22 30 16 Nd 20 16 27
Rudnick and FountainÕs Th/U V 151 230 98 Sm 3.9 3.5 5.3
ratios are close to bulk Earth. Cr 119 185 126 Eu 1.2 1.1 1.3
WedepohlÕs Th/U ratio of 5 Co 25 29 24 Gd 3.6 3.3 4.0
may be more appropriate. Ni 51 105 56 Tb 0.56 0.6 .65
Figure 12.32 shows t h e Cu 24 75 25 Dy 3.5 3.7 3.8
primitive mantle-normalized Zn 73 80 65 Ho 0.76 0.78 0.8
incompatible element concen- Ga 16 18 15 Er 2.2 2.2 2.1
trations in the total continen- Ge 1.6 1.4 Tm 0.32 0.3
tal crust. All these estimates As 1 1.7 Yb 2 2.2 2.0
show the continental crust to Se 0.05 0.12 Lu 0.33 0.3 0.35
be enriched in incompatible Br 1.0 Hf 3.7 3 4.9
elements. This is consistent Rb 58 32 78 Ta 0.7 1 1.1
with the view that the crust Sr 325 260 333 W 1 1.0
was created by partial melting Y 20 20 24 Re, ppb 0.4 0.4
of the mantle. The elemental Zr 123 100 203 Os, ppb 0.005 0.05
enrichment pattern is not Nb 8.5 11 19 Ir, ppb 0.1 0.05
smooth, however. First, Nb Mo 1 1.1 Pt, ppb 0.4
(and Ta, not shown), and to a Ru, ppb 0.1 Au, ppb 3 2.5
lesser degree Ti, are not as en- Rh, ppb 0.06 Hg, ppb 40
riched as elements of similar Pd, ppb 1 0.4 Tl, ppb 360 520
incompatibility. Second, Pb is Ag, ppb 80 70 Pb 12.6 8 14.8
anomalously enriched. As we Cd, ppb 98 100 Bi, ppb 60 85
found earlier in the chapter, In, ppb 50 50 Th 5.6 3.5 8.5
these are features of island arc Sn 2.5 2.3 U 1.42 0.91 1.7
magmas and their presence R & F: Rudnick and Fountain (1995) as revised by Plank and Langmuir (in
here strengthens the notion press), T & M: Taylor and McLennan (1985, 1995), We: Wedepohl (1995),
that island arc volcanism has Shaw: Shaw et al. (1986), W & T: Weaver and Tarney (1984).
539 January 25, 1998

played an important role in creation of
the continental crust.
Rudnick and FountainÕs (1995) or
100 B B
Sample/Primitive Mantle
B WedepohlÕs (1995) estimated composi-

BBB
JB tions are significantly more enriched in in-
J JJJ J
J
B
JB
compatible elements than the estimated
B
J J BB composition of Taylor and McLennan
JJ B J BB (1985, 1995). Since the former two esti-
10 Others J JJB
B
J mates are based on observations while t h e
J Taylor and McLennan (1985) J
BB JJ
BB J
Taylor and McLennan estimate is based on
B Rudnick and Fountain (1995) assumptions about how the crust formed,
the former estimates should probably be
1 preferred. This strong enrichment of t h e
Cs RbBaThU KNb LaCePb Pr Sr NdZrHfSmEuTi YHo Yb
continental crust in incompatible elements
Figure 12.32. Comparison of incompatible element has important implications for the evolu-
enrichment of estimated bulk crustal composition of tion of both the crust and the mantle.
Taylor and McLennan (1985), Rudnick and Fountain Based on Rudnick and FountainÕs (1995) or
(1995) with other estimates (Weaver and Tarney, WedepohlÕs (1995) estimates, the crust is
1984; Shaw et al., 1986; Wedepohl, 1995). From more enriched in incompatible elements
Rudnick and Fountain (1995). than typical island arc andesite. Thus
subduction-related magmatism alone is
inadequate as a mechanism for creation of the continental crust. Based on these estimates, a substan-
tial proportion of the EarthÕs total inventory of lithophile incompatible elements are in the crust.
For example, using the primitive mantle values listed in Chapter 11, we can calculate that 50% or
67% of the EarthÕs total Rb is in the EarthÕs crust using Rudnick and FountainÕs and WedepohlÕs crus-
tal estimates respectively. Similar fractions of other highly incompatible elements such as Cs, B a ,
and Th are also concentrated in the crust. This implies that the melt must have been extracted from
over half the mantle to create the crust, and therefore that at least half the mantle is incompatible
element depleted. Thus differentiation of
the EarthÕs has been remarkably effi-
Relative Mass of Contiental Crust, %
cient.
Growth of the Continental 100 F AM
Crust
The composition of the crust gives us
important clues as to how it was created.
R&S
For example, one possible model for t h e 50 O'N V&J H&R
formation of the continental crust is t h a t D&W
it was produced by late accretion of a M&T
volatile-rich veneer when the Earth
formed. But the composition given in Ta-
ble 12.7 is clearly inconsistent with this 0
view: the crust is not systematically en-
riched in volatile elements. The crust i s
4 3 2 1 0
systematically enriched in incompatible Time, Ga
elements; this leads to the hypothesis Figure 12.33. Models of the rates of crustal growth. AM:
that the crust was created by partial Armstrong (1981a), R&S: Rymer and Schubert (1984), F:
melting of the mantle. However, even Fyfe (1978), D & W: DePaolo and Wasserburg (1979), M
here there are some inconsistencies, since & T: McLennan and Taylor (1982), OÕN: OÕNions and
a partial melt of the mantle is unlikely to Hamilton (1981), V & J: Veizer and Jansen (1979), H & R:
have as much SiO 2 as the crust (~60%). Hurley and Rand (1969). Adapted from Taylor and
McLennan (1985).
540 January 25, 1998

Also, the trace element composition of the
crust is not like that of any single type of
mantle-derived magma. This suggests
the growth of the crust has been a com-
plex process. The details of this process
are not yet fully understood. In this sec-
°W
140
tion, we examine what is known about
crustal growth.
The Pace of Crustal Growth 40°
The first question we might ask is how
N
has the crust grown with time? A variety
of answers to this question have been sug-
gested. These are illustrated in Figure
12.33. They can be broken into three 12 0
types: (1) growth rate increasing through °W
time, e.g., curve V & J (=Vezier and Jan-
sen) and H & R (Hurley and Rand), (2)
approximately linear growth through
time, such as the curve marked O'N (= 100°W 80°W 60° W
O'Nions), and (3) early rapid growth fol-
lowed by later slow growth or no growth,
e.g., curves Am (=Armstrong) and F
<440 Ma 800-1700 Ma
(Fyfe). 800-1300 Ma 1700-2350 Ma
Early studies favored accelerating 2350-2700 Ma
crustal growth through time. Age prov-
1300-1700 Ma
inces for North America based on radio- Figure 12.34. Age provinces of the North American con-
metric dates compiled by Hurley and tinent (after Hurley and Rand, 1969).
Rand (1969) are shown in Figure 12.34. Based on this, Hurley suggested most of the crust was produced
in the last 1 or 2 Ga. However, subsequent work utilizing U-Pb dating of zircons and Sm-Nd model
ages has substantially changed this view. U-Pb zircon ages are not so easily reset by metamorphism
as are the Rb-Sr ages relied on by Hurley. Zircon dating has revealed large areas of the continental
crust that were created more than 2.7 Ga ago. Zircon dating as also identified limited areas of crust
that are old than 3.5 Ga. The oldest dates are from two localities in Australia. Zircons from Mt. Nar-
ryer, Western Australia give ages of around 4.15 Ga (Frounde et al., 1983). Slightly older ages (4.2-
4.3 Ga) were determined on zircons from a second, nearby locality, the Jack Hills. These zircon ana-
lyzes were done by ion probe rather than conventional mass spectrometry and were initially contro-
versial for that reason, but they ages are now generally accepted. In both cases, the zircons are found
in metasedimentary rocks of much younger age (3Ð3.5 Ga), and only a handful of the zircons in these
rocks are this old. While these ancient zircons have been preserved, the rocks in which they crystal-
lized apparently have not. Nevertheless, the zircons are relatively U-rich, suggesting they come
from silica-rich rocks typical of true continental crust rather than more mafic rocks.
More recently, zircons in the Acasta gneisses from the Slave Province in Canada have been dated
at 3.96 Ga by this method (Bowring et al., 1989). These ages are interpreted as the age of crystalliza-
tion of protoliths of these gneisses. Interestingly, their initial εNd are negative, indicating deriva-
tion from a light rare earth-enriched source, which could be even older crust. Thus there is evidence
that at least some continental crust formed very early in EarthÕs history.
In Chapter 11, we found that the many very ancient rocks have positive initial εNd ratios. For ex-
ample, the Isua rocks of East Greenland, which have crystallization ages around 3.8 Ga, have initial
εNd around +1 to +3; rocks of similar age from India have similar initial εNd values. This implies
their sources had high Sm/Nd ratios, and therefore were LRE-depleted for a substantial amount of
541 January 25, 1998

time, hundreds of millions of
years, before the crystalliza- τDM > 2.7 Ga
tion of these rocks. This re- F F J
quires that a complimentary B . 3 G a JJJ
B 2
2.0–
B
LRE enriched reservoir, pos- B B J
B τ DM =
B J
sibly an early continental J
San A
crust, formed before 4.0 Ga.
B J J
How large this crust reser- 1 2
0 Ga
ndre
voir was, however, depends 2 .
–
= 1.8
as F
on how large the volume of
τ
H
M
aul
depleted mantle was, which B JJJ D a
1.8 G
t
is not constrained (McCulloch B B
B BB
1 . 7 –
and Bennett, 1994). Never- B B =
τ
B H
theless these data provide Pac BB BB BB HH
DM 3
ific B HHHH
indirect evidence that a t Oc BB
least crust formation began ean
well before 4.0 Ga. H HH H
H H
Sm-Nd model ages also H
H
have the power to Òsee 0 100 200 km
throughÓ metamorphism and
establish Òcrust formation Figure 12.35. Isotopic provinces, based on crustal residence times
agesÓ, as we found in Chap- (τDM) of the Western U.S. From Bennett and DePaolo (1987).
ter 8. The example of t h e
western U. S. illustrates this point. Figure 12.35 is a map of the Western U. S. showing contours of Nd
crustal residence times (τDM). The data define 3 distinct provinces and suggest the existence of several
others. There is a general similarity to Hurley's map (Figure 12.34), but there is greater detail, and
the ages are generally older.
Figure 12.36 shows the initial εNd values of the granites from the three numbered provinces of Fig-
ure 12.35 plotted as a function of their crystallization age. Despite the variations in crustal residence
times, the crystallization ages indicate Provinces 1-3 all formed between 1.65 and 1.8 Ga. Only t h e
εNd from Province 3 plot close to the depleted mantle evolution curve. From this we can conclude t h a t
only Province 3 was a completely new addition to the crustal mass at that time. Initial εNd for t h e
remaining provinces plot below the depleted mantle evolution curve, suggesting they are mixtures of
new mantle-derived material and older crust. The crustal residence ages of these provinces are older
than the crystallization ages because they are mixtures of mantle and older crust.
In each province there have been subsequent episodes of magmatism. However, the initial εNd l i e
along the same growth trajectory as the older rocks. This suggests that magmatism in these subse-
quent episodes simply recycled pre-existing crustal material and there were no new additions to crus-
tal mass from the mantle. Thus HurleyÕs map (Figure 12.34) must reflect orogeny in which the radio-
genic clocks are reset rather than new crustal additions.
At the other extreme of the continental growth question, Armstrong (1968, 1981a) argued that t h e
mass of the crust has remained nearly constant of the past 4 Ga or so. Armstrong recognized that new
crust has been continually created through time, but he argued that the rate of crustal creation was
balanced by an equal rate of crustal destruction throughÊerosion and subduction of sediment. Three
lines of evidence support this view. First, Armstrong cited the absence of ancient pelagic sediment
anywhere on the EarthÕs surface and argued that most sediment must therefore be subducted. In sup-
port of this, he cited estimates of sediment subduction rates by several workers that are sufficient to
balance crustal growth. Second, he pointed out that neither sea level nor the average thickness of
stable continental cratons have changed with time, and therefore that continental volume also must
not have changed. Finally, he pointed out that positive εNd in the earliest rocks requires a very early
depletion of the mantle (Figure 10.19), most likely through generation of continental crust. Further-
542 January 25, 1998

+10 more, the subsequent nearly linear growth of
Deplete εNd is consistent with constant continental mass
d Mantle and inconsistent with continental growth
1.65-1.8 (Armstrong, 1981b, DePaolo, 1981b).
+5 JJ
J 2.0-2.3 Because sediments effectively sample rocks
JG
J G
J GG of a wide variety of ages, isotope ratios in
G
G
G
JJ G
JG
JJ
sediments provide a means of studying crustal
0 JG
G
G JJ
JJJ evolution. By studying sediments of various
JJ
JG
GJ
J ages, we should be able to draw some infer-
J
J ences about the rates of continental growth.
-5 J Goldstein et al. (1984) found that the mean
Early Paleozoic Rocks

3 crustal residence time (τ DM ) of the modern
εNd 2 river particulates was 1.7 Ga, which they in-
Anorogenic Rocks
E 1 terpreted as the mean age of the crust now be-

-10 J
E
E ing eroded. However, they estimated t h e
J
E
J mean crustal residence time of the entire sedi-
mentary mass to be about 1.9 Ga. Figure 12.37
-15 compares the stratigraphic age* of sediments
J
E
EE
J
E No New Input (a)
-20 4.0 Decreasing
Rate
3.0 Uniform
τCR(Ga)
0.0 0.5 1.0 1.5 2.0 Rate

τCR = τSTRAT
Crystallization Age, Ma 2.0 (first cycle)
Figure 12.36. εNd (initial) as a function of crystal- Increasing
lization age of Western U.S. Groupings 1 2, and 3 Systematics Of Rate
refer to provinces shown in Figure 12.35. From 1.0 Sediment Source Evolution
Bennett and DePaolo (1987). 4.0 3.0 τ 2.0 1.0 0
with their crustal residence ages. We expect that crus- STRAT (Ga)
tal residence ages will be somewhat older than t h e
stratigraphic ages. Only when a rock is eroded into t h e 4.0 (b)
sedimentary mass immediately after its derivation ●
from the mantle will its stratigraphic (τST) and crustal A

τCR(Ga)
●
3.0
residence age (τCR) be equal. ▲
▲
▲
▲
B
τCR = τSTRAT
▲
The top diagram illustrates the relationships be- ●
tween τST and τCR that we would expect to see for various
●
2.0 ●● ●●
crustal growth scenarios, assuming there is a relation- Metasediments ●● ●
●●
● ●●
▲ Canadian Shield ▲ ●
ship between the amount of new material added to t h e Composites ●
1.0 ●
continents and the amount of new material added to t h e ● Sediments
sedimentary mass. If the continents had been created 4.0

Ga ago and if there had been no new additions to conti- 4.0 3.0τ 2.0 1.0 0
nental crust since that time, then the crustal residence STRAT (Ga)
time of all sediments should be 4.0 Ga regardless of Figure 12.37. Relationship between strati-
stratigraphic age. This is illustrated by the line l a - graphic age of sediments and the crustal
beled ÒNo New InputÓ. If, on the other hand, the rate of residence age of material in sediments.
continent growth through time has been uniform since 4.0 From Goldstein et al., 1984).
*The stratigraphic age is the age of deposition of the sediment determined by conventional
geochronological or geological means.
543 January 25, 1998

Ga, then τST and τCR of the sedimentary mass should lie along a line with slope of 1/2, which is t h e
line labeled ÒUniform RateÓ. The reason for this is as follows. If the sedimentary mass at any given
time samples the crust in a representative fashion, then τCR of the sedimentary mass at the time of its
deposition (at τST) should be (4.0-τST)/2 , i.e., the mean time between the start of crustal growth
(which we arbitrarily assume to be 4.0 Ga) and τST. A scenario where the rate of crustal growth de-
creases with time is essentially intermediate between the one-time crust creation at 4.0 and t h e
uniform growth rate case. Therefore, we would expect the decreasing rate scenario to follow a trend
intermediate between these two, for example, the line labeled 'Decreasing Rate'. On the other hand,
if the rate has increased with time, the τCR of the sedimentary mass would be younger than in t h e
case of uniform growth rate, but still must be older than τST, so this scenario should follow a path be-
tween the uniform growth rate case and the line τST = τCR, for example, the line labeled ÒIncreasing
RateÓ.
Line A in Figure 12.37b is the uniform growth rate line with a slope of 1/2. Thus the data seem to
be compatible with a nearly uniform rate of growth of the continental crust, such as the line labeled
ÒOÕNÓ in Figure 12.33 or even with an increasing growth rate. However, Goldstein et al. noted sedi-
mentary mass is cannibalistic: sediments are eroded and redeposited. Line B represents the evolution
of the source of sediments where crustal growth is constant but erosion and re-deposition of old sedi-
ments occurs. The situation is further complicated when crust is destroyed through erosion and subduc-
tion of marine sediments. Goldstein et al. (1984) concluded their data could be consistent with either
nearly constant continental growth rate or a nearly constant continental mass if the rate rate of conti-
nent-to-mantle recycling decreased through time.
In summary, formation of the continental crust began 4.0 Ga ago or earlier. The rarity of rocks of
this age, however, suggests either that crustal growth was initially quite slow, or that much of this
early continental crust may have been destroyed and recycled into the mantle. Since then, there have
been continued additions to the continental crust. The degree to which destruction of crust through
erosion and subduction may have balanced these additions is unclear and remains highly controver-
sial. The average age of the presently existing continental crust is around 2.2 Ga. Models, such as
that of Hurley, in which the rate of crustal growth is increasing cannot be correct. Reconnaissance
Sm-Nd studies of the sort illustrated in Figure 12.35 show that the crust is created in large blocks, and
that the rate has been somewhat episodic. Judging from these kinds of studies, the late Archean and
early Proterozoic (the period between roughly 3.0 and 1.8 Ga) appears to have been a time of rela-
tively rapid crustal growth (e.g., Patchett and Arndt, 1986; McCulloch, 1987). On the other hand,
crustal growth rates in the Phanerozoic (the last 550 Ma) appear to be lower than during earlier
times.
Mechanisms of Crustal Growth
We have seen how the composition of the continental crust is qualitatively consistent with t h e
crust having been created by partial melting of the mantle. But what are the details of this process?
In what tectonic setting has the crust been produced? We can identify a number of possible mecha-
nisms for the creation of the continental crust. These are as follows:
One way to rationalized this equation is to think of newly deposited sediment at τST as a 50-50 mixture
of material derived from the mantle at 4.0 Ga and τST. The equation for the τCR of this mixture would be:
4.0Ê+ÊτST
τCR = .
2
4.0Ê+ÊτST 4.0ÊÐÊτST
At time of deposition, its crustal residence age would have been: τCR = Ð τST = .
2 2
You could satisfy yourself that a mixture of material having τCR of all ages between 4.0 Ga and τST
would have the same τCR as given by this equation.
544 January 25, 1998

¥ Subduction-related volcanism. The Andes Moun-
Table 12.8. Comparison of tains with their magnificent and very active andesitic
Continental Crust with Andesite volcanos represent a modern example of how the crust
Continental Average grows above a subduction zone.
Crust Andesite ¥ Accretion of oceanic crust and oceanic plateaus.
SiO 2 57.3-63.2 59.7 Oceanic crust is usually subducted and returned to t h e
TiO 2 0.6-0.9 0.7 mantle, but anomalously thick crust, such as that of t h e
Al 2O 3 14.8-16.1 17.1 Ontong-Java Plateau of the Western Pacific may be dif-
FeO 5.6-9.1 6.3 ficult to subduct. The alternative is that is accretes to
MnO 0.08-0.18 0.12 continental margins, becoming part of the continental
MgO 2.8-5.3 3.2 mass.
¥ Continental volcanism. Rifting and mantle plumes
CaO 4.7-7.4 6.6
are two causes of volcanism unrelated to subduction.
Na 2O 3.1-3.2 3.3
Flood basalts, probably produced by the start-up phase
K 2O 1.1-2.4 1.5 of mantle plumes, have occasionally been erupted in
P2O5 0.14-0.20 0.19 tremendous volumes.
Continental Crust from Table 12.7. Average ¥ Underplating. Magmas erupted on continents may
medium-K orogenic andesite from Gill (1981). represent only the tip of the iceberg so to speak. Be-
cause of the low density of the continental crust, basalt magma may stagnate at the crust-mantle
boundary. In this scenario, the crust primarily grows from the bottom down.
¥ Intrusion by small-degree melts. Small degree melts of shallow mantle under the continents t h a t
are highly enriched in incompatible elements, such as kimberlites and lamprophyres, may not con-
tribute significant volumes to the continents, but may contribute a disproportionate fraction of t h e
incompatible element inventory.
Of course, all these mechanisms have contributed to crustal growth to varying degrees. The real
question is which mechanism is most important?
There are a number of cogent arguments favoring subduction-related magmatism as the principal
mechanism of crustal growth. First, in a qualitative way, the continental crust has a major element
composition similar to andesites erupted on continental margins and in island arcs. This point is made
in Table 12.8, which compares average ÒorogenicÓ (i.e., island arc or continental margin) andesite to
the estimates of crustal composition in Table 12.7. The resemblance is strong; only the Al concentra-
tion lies outside the range of estimated concentrations in the crust.
The similarity in composition also extends to incompatible elements. This point is made in Figure
12.38, which compares the composition of a siliceous andesite from the Banda arc (Indonesia) with
the range of estimates of crustal composition. While the match is not an exact one, the continental
crust shares the strong enrichment in highly incompatible elements, the relative depletion in Nb and
Ta, and the excess enrichment in Pb observed in most island arc magmas. We should point out, how-
ever, that to some extent, the incompatible element enrichment of island arc magmas results from a
Òcontinental componentÓ in them, originating either through sediment subduction or assimilation of
crust. Thus in trying to understand the incompatible element enrichment of the crust, postulating an
origin by subduction zone magmatism begs the question in a certain sense.
Finally, as we noted earlier, subduction-related magmatism is the most important mechanism of
crustal growth at present, and probably throughout the Phanerozoic as well. Many igneous and meta-
igneous rocks of Proterozoic age also have chemical features suggesting they were produced in subduc-
tion settings. Bennett and DePaolo (1987) concluded the provinces in the Southwest U.S. shown in
Figure 12.35 formed by successive accretion or growth of island arcs on the edge of the pre-existing Ar-
chean Wyoming Craton to the north. Judging from their older Sm-Nd model ages, the northeastern
arcs contain a substantial component of older crust derived from the craton. This could have occurred
through erosion and subduction, or, if the volcanos were built directly on the continent, through as-
similation of crust. As new Proterozoic crust was built outward from the continent, it screened subse-
545 January 25, 1998

quent arcs from the contribu-
200 tion of material from the Ar-
Primitive Mantle Normalized Concentration
chean crust. A similar effect

100 H H has been observed in the Pro-
H
HH terozoic provinces of Canada.
H H Accretion of intra-oceanic
arcs, modern examples of
H which include the Marianas
H and the Aleutians, may also
HH
10 H H H H contribute to crustal growth.
H H H HH H H HH Indeed, major segments of con-
H HH tinental crust appear to have
grown so rapidly that t h e
could not have been created
in a single arc (Rymer and
Schubert, 1986); they there-
1 fore may have originated as
CsRb BaTh U NbTa KLa CePb Pr SrNdSmZrHf EuGdTbDyYHo ErTmYbLu several simultaneously ac-
Figure 12.38. Comparison of incompatible element concentrations in tive arcs subsequently ac-
a siliceous andesite from the Banda arc (red triangles) with t h e creted to a continental mar-
range of estimated concentrations in the continental crust (gray gin. Even in this case, how-
field). Both share a relative depletion in Nb and Ta and a relative ever, global magma produc-
enrichment in Pb. tion rates at subduction zones
during at least some parts of
the Proterozoic must have been 6 to 10 times higher than during the Phanerozoic (Rymer and
Schubert, 1986). These high growth rates may mean some mechanism other than subduction-related
volcanism was responsible for growth of much of the crust.
Though subduction-related volcanism may have been the dominant mode of crustal growth in t h e
Proterozoic and Phanerozoic, other mechanisms have played a role. The Wrangalia Terrane in N W
British Columbia and Alaska, is widely considered to consist in part of oceanic plateaus. The pla-
teaus were produced over a mantle plume in Paleozoic time and later accreted to the North American
continent by plate tectonic processes. The Coast Ranges of Oregon represent an example of more re-
cently accreted oceanic crust. Mantle plumes surfacing beneath continents also produce magmas t h a t
add mass to the continents. The most voluminous eruptions occur in the initial stages of the plume,
when the large buoyant plume head approaches the surface. Under these circumstances, enormous
volumes of basalt erupt. Examples of such flood basalts include the Siberian Traps, the Karoo of
South Africa, the Deccan of India, the Parana of Brazil, and the Columbia River of the NW U. S .
Gravity anomalies suggest even greater volumes of basaltic magma were trapped at deep crustal lev-
els.
Continental rifts can also be sites of voluminous eruption of basaltic magma. A well-documented
example is the Proterozoic Keweenawan or Mid-Continent Rift of the U. S., which formed some 1 to
1.2 Ga ago. Though now mostly covered by Phanerozoic sediments, where it is exposed the rift con-
sists of a trough 150 km wide and 1500 km long filled with up to 15 km of volcanics, primarily basalt,
and clastic sediments derived from them. Modern examples of continental rifts include the Rio
Grande Rift of New Mexico and the East African Rift.
Crustal Growth in the Archean
The case for crustal growth through subduction-related volcanism is much less strong in the Ar-
chean. While andesites dominate island-arc volcanism at present, andesites seem to be relatively
rare in the Archean. Taylor and McLennan (1995) remarked that Òthe composition of the Archean
upper crust stands in marked contrast to that of the post-Archean crust.Ó This suggests the Archean-
546 January 25, 1998

Proterozoic boundary marks a
change in the manner in
which crust was generated.
Since nearly half the conti-
nental crust was created by 100
the end of the Archean, un- Finland
derstanding the genesis of Yilgarn
Sample/Chondrites
Archean crust is important.
There are two principle
types of Archean terranes.
The first is greenstone belts,
of which the Abitibi Belt in 10
the Superior Province of
Canada is a classic example. Amitsoq
They consist of thick se- Ancient Gneiss Complex
quences of volcanic and sedi- Shaw Batholith, Pilbara
mentary rocks in elongate ba-
sins punctuated by circular or
Mt. Edgar Batholitn, Pilbara
elongate granitic batholiths.
They have experienced only La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb
low grade metamorphism Figure 12.39. REE patterns of typical Archean tonalites and trond-
(hence the term greenstone: hjemites. From Taylor and McLennan (1985).
chlorite-serpentinite is t h e
typical metamorphic grade).
Tholeiitic basalt and komatiite are the predominant volcanics, but more siliceous ones are also com-
mon. The combination of mafic volcanics and granitic intrusives makes greenstone belts distinctly bi-
modal in composition.
The presence of komatiites is significant. They are ultramafic lavas erupted at temperatures of
1400 to 1600¡ C, much hotter than basalts, whose eruption temperatures are typically 1100 to 1200¡ C.
Their chemistry indicates they are products of large extents of melting, perhaps up to perhaps 40%.
Komatiites are largely restricted to the Archean; there are a few early Proterozoic occurrences and
only one documented Phaner-ozoic occurrence. The absence of komatiites in the latter part of EarthÕs
history undoubtedly reflects secular cooling of the mantle.
The other type of Archean crust, the high-grade gneiss terranes, such as West Greenland, is much
different. They typically consist of felsic gneisses and sedimentary and volcanic rocks metamor-
phosed at amphibolite to granulite grade. These terranes do not seem to be simply a highly meta-
morphosed version of the greenstone belt terranes as they differ in structure and in sedimentary f a -
cies. The felsic gneisses consist of metamorphosed plutonic rocks of the so-called ÒTTGÓ (Tonalite*Ð
Trondhjemite ÐGranodiorite) suite, characterized by higher Na/K ratios than are found in most post-
Archean granitoids. Many have steep REE patterns (Figure 12.39) and no Eu anomalies. These rare
earth characteristics indicate they formed by partial melting at great depth (>60 km), where garnet,
rather than plagioclase, was present in the residua. The presence of garnet would account for their
steep heavy rare-earth patterns and the absence of plagioclase would explain the lack of Eu-anoma-
lies. Positive initial εNd values of many of these gneisses, particularly the oldest of them, indicates
they formed directly from the mantle, or more likely, by partial melting of basalts that had them-
* Tonalite, or quartz diorite, is a plutonic rock of more or less andesitic composition, having a pre-
dominance of plagioclase over orthoclase (i.e., K-poor diorite).
Trondhjemite is a plutonic rock of more or less granitic composition, but is poorer in orthoclase and
richer in plagioclase than granite senso stricto (i.e., K-poor granite).
547 January 25, 1998

selves formed only shortly before. The oldest preserved crust (>3.5 Ga) invariably consists of high-
grade gneiss terranes.
A number of tectonic models have been suggested for the origin of these two terranes. The green-
stones themselves are clearly of mantle derivation, though there is evidence in a few cases of crustal
contamination. Hypotheses for their origin include continental rifts, island-arcs or continental mar-
gins, back-arc environments, and oceanic plateaus formed over mantle plumes. Many have argued
that the presence of komatiites precludes formation of greenstone belts in island arcs. Because
komatiites require extremely high temperatures for their generation, the relatively cool subduction
zone environment seems an unlikely setting for their production. It seems more likely they formed
over a zone of mantle upwelling, such as a mid-ocean ridge or mantle plume. In the Abitibi Belt,
which formed between 2730 Ma and 2690 Ma, the volcanics lack Nb-Ta anomalies characteristic of
subduction related magmas. This, and the rapid production of great volumes of magma, led Vervoort
et al. (1993) to argue this belt originally formed as an oceanic plateau produced by a mantle plume.
The granitic intrusions, which are generally slightly younger, however, do possess Nb-Ta anomalies,
and hence may have formed near a subduction zone. Vervoort et al. (1993) suggested these granites
were generated over a subduction zone and intruded the greenstone belt as it was accreted to the grow-
ing Superior craton. Early Proterozoic greenstone belts of the West African Craton may also represent
accreted oceanic plateaus (Abouchami et al., 1990). Volcanics of other greenstone belts do, however,
have Nb-Ta anomalies, suggesting they formed in an island arc setting (Condie, 1989).
The origin of the high-grade gneiss terranes is perhaps more problematic. Some have argued t h a t
the TTG suites are mantle melts. While it is possible to create siliceous melts of peridotite under con-
ditions of shallow depth and high water pressure, there are a number of problems with this idea, and
it seems more likely the TTG suite represents melts of a basaltic precursor (Taylor and McLennan,
1985; Ellam and Hawkesworth, 1988; Rudnick, 1995). There are two basic variants of this idea. First,
they could represent melts of basaltic lower crust. The basaltic precursor could have formed in any
number of tectonic environments. Alternatively, they could represent melts of subducted oceanic crust
(Martin, 1986; Drummond and Defant, 1990). Higher temperatures and more rapid plate movements
in the Archean may have meant that melting of subducting oceanic crust, which is rare today, was
common then. In either case, melting must have occurred under sufficient pressure that the basaltic
source was metamorphosed to amphibolite or eclogite.
Refining the Continental Crust
We can summarize by saying that the continental crust has been created by partial melting of t h e
mantle. Subduction-related volcanism appears to have been the principal environment in which new
crust has formed. Rudnick (1995) estimates that 65 to 90% of the crust has been produced in this man-
ner. Archean crust is compositionally distinct. This may reflect a greater proportion of crust produced
by some other mechanism (e.g., melting of mantle plumes), a difference in the composition of subduc-
tion-related magmas due to higher temperatures, or both.
One serious problem remains, however, namely that the composition of the continental crust does
not match that of a primary mantle melt. Primary mantle melts are, with rare exceptions, basaltic.
If the continental crust were simply a melt of the mantle, it should have a basaltic composition, as
does the oceanic crust. Compared to basalt, which is the principal melting product of the mantle, it is
too rich in SiO 2, too poor in MgO and CaO (and probably FeO), has too low a Mg# and too high a
Rb/Sr ratio, and is generally too rich in incompatible elements. This is true irrespective of the envi-
ronment in which these melts are produced. Compared to typical island arc basalts, its Al2O3 is also
too low, as is its Sr/Nd ratio. As we stated earlier, basalt, not andesite, is the primary melt produced
in subduction zones, and intra-oceanic arcs have basaltic bulk compositions. Andesites, while abun-
dant in island arcs and even dominant in continental margin settings such as the Andes, are products of
intracrustal fractional crystallization of basaltic parents, which in the case of continental margins, is
often accompanied by assimilation of pre-existing crust (e.g., Hawkesworth et al., 1982). Further-
more, andesites produced in this way have concentrations of highly compatible elements, such as N i
548 January 25, 1998

and Cr, that are much lower than those of the continental crust. Similarly, mantle-derived magmas
produced at rifts or over mantle plumes are predominantly basaltic. While some andesitic magmas
can be produced by melting of mantle peridotite or by reaction of basaltic magma with peridotite
(e.g., boninites), such melts are rare and require very special conditions for their generation. Fur-
thermore, their compositions do not match that of continental crust. These observations have led to
the suggestion that the continental crust has been ÒrefinedÓ in some way; that is, some other process or
processes have operated to transform the composition of the crust.
One possible explanation for this problem is that lower continental crust, consisting of mafic cumu-
lates or partial melting residues, detaches, or ÒdelaminatesÓ, along with underlying subcontinental
mantle lithosphere (Kay and Kay, 1991). The idea is that when the crust is thickened the 50 km or
more in collisional zones such as the Andes or the Himalayas, mafic lower crust will be transformed
to garnet granulite and eclogite, greatly increasing its density. Under these circumstances, it may
simply detach from the overlying crust and be swept away in asthenospheric circulation along with
the underlying mantle lithosphere. The lithosphere would be replaced by hotter, lower density as-
thenosphere, which would produce uplift and basaltic volcanism. Kay and Kay argued that pre-
cisely these phenomena are observed in the Puna of Argentina, where the crust was thickened to as
much as 65 km in the Miocene. As we noted in the previous chapter, McKenzie and OÕNions (1983)
suggested that detached subcontinental lithosphere could later form the source of mantle plumes.
Since the lower continental crust is clearly more mafic and poorer in incompatible elements than
the bulk crust (compare Tables 12.6 and 12.7), loss of lower crust to the mantle clearly provides a
mechanism for adjusting the composition of the crust in the desired direction. On the other hand, t h e
great age of many xenoliths from the subcontinental lithosphere attests to it stability (it would be
difficult to detach the lower crust without also detaching the underlying mantle lithosphere). As we
noted in Chapter 11, peridotite xenoliths derived from the mantle lithosphere often have the same
age as the overlying crust (however, none of these ancient xenoliths have been found in areas where
lithospheric detachment is thought to have occurred). Furthermore, the distinctive OsÐNd isotope
systematics of subcontinental lithosphere have not been found in any mantle plumes yet. Neverthe-
less, detachment of the lower crust is a viable and interesting hypothesis, and one that deserves fur-
ther study.
Yet another possibility is that the composition of the crust has been refined through erosion and
subduction of weathering products. During weathering, certain elements, most notably Mg, preferen-
tially go into solution. Mg is carried to the oceans where it is removed from solution by reaction with
the basaltic oceanic crust during hydrothermal activity at mid-ocean ridges (discussed in Chapter
15). Depending on the fluxes assumed, it would appear that this process has the potential for remov-
ing substantial amounts of Mg from the crust and hence make it less mafic (see Problem 12.8 at the end
of the chapter). In the same way, other soluble elements, such as Ca and Na, would be removed from
the crust.
While this process probably contributes to refining the composition of the crust, its effect is proba-
bly quite limited. If it were important, we would expect the major element chemistry of the crust to be
enriched in the least soluble elements (Al, Si, and Fe) and systematically depleted in the most solu-
ble ones (Na, Mg, Ca) compared to typical igneous rocks. This is not the case. Among trace elements,
it would not increase the extent of LRE enrichment, nor is it clear that this mechanism could increase
the Rb/Sr ratio of the crust as required (though it would decrease Sr/Nd).
OÕNions and McKenzie (1988) pointed out that the incompatible element abundances in the conti-
nental crust resemble that of a very small (~1%) degree melt of the mantle. They argued that such
small degree melts will be generated during episodes of continental extension. The melts would mi-
grate into the subcontinental lithospheric mantle or lower crust and freeze there. Subsequent melting
of this lower crust would produce strongly incompatible element enriched granites, which are common
components of the upper crust. Interestingly, if this mechanism does indeed have a significant influ-
ence on the incompatible element budget of the continents, it would mean that Sm-Nd model ages sig-
nificantly underestimate the age of the continents. Relatively little testing of this interesting h y -
549 January 25, 1998

pothesis has been carried out, so it is difficult to judge the importance of this process. Even if the in-
compatible element abundances in the crust can be explained this way, they hypothesis does not ad-
dress the problems with major element composition listed above.
References and Suggestions for Further Reading
Abouchami, W., M. Boher, A. Michard and F. Albar•de. 1990. A major 2.1 Ga event of mafic magma-
tism in West Africa Ñ an early stage of crustal accretion. J. Geophys. Res. 95: 17605-17629.
Armstrong, R. L. 1968. A model for the evolution of strontium and lead isotopes in a dynamic Earth.
Rev. Geophys. 6: 175-199.
Armstrong, R. L. and J. A. Cooper. 1971. Lead isotopes in island arcs. Earth Planet. Sci. Lett. 35: 27-37.
Armstrong, R. L.1981a. Radiogenic isotopes: the case for crustal recycling on a near-steady-state no-
continental-growth Earth. in The Origin and Evolution of the Earth's Continental Crust, ed. S. M.
and B. F. Windley. 259-287. London: The Royal Society.
Armstrong, R. L. 1981b. Comment on "Crustal growth and mantle evolution: inferences from models of
element transport and Nd and Sr isotopesÓ. Geochim. Cosmochim. Acta. 45: 1251.
Bennett, V. C. and D. J. DePaolo, Proterozoic crustal history of the western United States as deter-
mined by neodymium isotope mapping. Bull Geol. Soc. Amer., 99: 674-685, 1987.
Bowring, S. A., I. S. Williams, and W. Compston, 1989. 3.96 Ga gneisses from the Slave province,
Northwest Territories, Canada. Geology, 17: 971-975.
Bowring, S. A., J. E. King, T. B. Housh, C. E. Isachsen, and F. A. Podosek, 1989. Neodymium and lead
isotopic evidence for enriched early Archaean crust in North America. Nature, 340: 222-225.
Brophy, J. G. and B. D. Marsh. 1986. On the origin of high-alumina arc basalt and the mechanics of
melt extraction. J. Petrol. 27: 763-789.
Carmichael, I. S. E., F. J. Turner and J. Verhoogen. 1973. Igneous Petrology. New York: McGraw-Hill.
Clarke, F. W. 1924. The Data of Geochemistry. U. S. Geol. Surv. Bull. 770:
Condie, K. C. 1989. Plate Tectonics and Crustal Evolution. Oxford: Pergamon.
Defant, M. J. and M. S. Drummond. 1990. Derivation of some modern arc magmas by melting of young
subducted lithosphere. Nature. 347: 662-665.
DePaolo, D. J. 1980. Crustal growth and mantle evolution: inferences from models of element transport
and Nd and Sr isotopes. Geochim Cosmochim Acta. 44: 1185-1196.
DePaolo, D. J. 1981a. Trace element and isotopic effects of combined wallrock assimilation and frac-
tional crystallization, Earth. Planet. Sci. Lett., 53: 189-202.
DePaolo, D. J. 1981b. Crustal growth and mantle evolution: inferences from models of element trans-
port and Nd and Sr isotopes (reply to a comment by R. L. Armstrong). Geochim. Cosmochim. Acta.
45: 1253-1254.
DePaolo, D. J. and G. J. Wasserburg. 1979. Sm-Nd age of the Stillwater complex and the mantle evo-
lution curve for neodymium. Geochim Cosmochim Acta. 43: 999-1008.
Drummond, M. S. and M. J. Defant. 1990. A model for trondhjemite-tonalite-dacite genesis and crustal
growth via slab melting: Archean to modern comparisons. J. Geophys. Res. 95: 21503-21521.
Eade, K. E. and W. F. Fahirg, Chemical evolution trends of continental platesÑa preliminary study
of the Canadian Shield, Geol. Serv. Canada Bull., 179, 1971.
Ellam, R. M. and C. J. Hawkesworth. 1988. Is average continental crust generated at subduction zones?
Geology. 16: 314-317.
Faure, G., 1986. Principles of Isotope Geology, 2nd ed., Wiley & Sons, New York, 589p.
Frounde, D. O., T. R. Ireland, P. D. Kinney, R. S. Williams, W. Compston, A. R. Williams, and J. S .
Myers, 1983. Ion microprobe identification of 4,100-4,200 Myr-old detrital zircons. Nature, 304: 616-
618.
Fyfe, W. S. 1978. Evolution of the Earth's crust: modern plate tectonics to ancient hot spot tectonics?
Chem. Geol. 23: 89.
Gill, J. B. 1981. Orogenic Andesites and Plate Tectonics. Berlin: Springer Verlag.
550 January 25, 1998

Goldschmidt, V. M. 1933. Grundlagen der quantitativen Geochemie. Fortschr. Mineral. Kristall.
Petrolgr. 17: 112-156.
Goldstein, S. J., R. K. O'Nions and P. J. Hamilton, A Sm-Nd isotopic study of atmospheric dusts and
particulates from major river systems. Earth Planet. Sci. Lett., 70, 221-236, 1984.
Hawkesworth, C. J., M. Hammill, A. R. Gledhill, P. van Calsteren and G. Rogers. 1982. Isotope and
trace element evidence for late-stage intra-crustal melting in the High Andes. Earth Planet. Sci.
Lett. 58: 240-254.
Hawkesworth, C. J., J. M. Hergt, R. M. Ellam and F. McDermott. 1991. Element fluxes associated with
subduction related magmatism. Phil. Trans. R. Soc. Lond. A. 335: 393-405.
Hofmann, A. W. 1988. Chemical differentiation of the Earth: the relationship between mantle, con-
tinental crust, and oceanic crust. Earth Plan Sci Lett. 90: 297-314.
Hurley, P. M., and J. R. Rand, Pre-drift continental nuclei, Science, 164, 1229-1242, 1969.Gill, J. B.,
1981, Orogenic Andesites and Plate Tectonics, Springer-Verlag, Berlin, 390 p.
James, D. E., 1981. The combined use of oxygen and radiogenic isotopes as indicators of crustal contam-
ination, Ann. Rev. Earth. Planet. Sci., 9: 311-344.
Karig, D. E. and R. W. Kay, Fate of sediments on the descending plate at convergent margins, P h i l .
Trans. R. Soc. Lond. Ser. A, 301, 233-251, 1981.
Kay, R. W. 1978. Aleutian magnesian andesites: metls from subducted Pacific Ocean crust. J. Volcanol.
Geotherm. Res. 4: 117-132.
Kay, R. W. and S. M. Kay. 1991. Creation and destruction of lower continental crust. Geol. Rund. 80:
259-278.
Kay, S. M. and R. W. Kay. 1994. Aleutian magmas in space and time. In The Geology of Alaska, T h e
Geology of North America Vol. G-1, G. Plafker and H. C. Berg. ed., pp. 687-722. Boulder: Geol. Soc.
Am.
Kay, S. M., R. W. Kay, and G. P. Citron. 1983. Tectonic controls on tholeiitic and calc-alkaline mag-
matism. J. Geophys. Res. 87: 4051-4072.
Klein, E. M. and C. H. Langmuir. 1987. Ocean ridge basalt chemistry, axial depth, crustal thickness
and temperature variations in the mantle. J. Geophys. Res. 92, 8089-8115.
Klein, E. M. and C. H. Langmuir. 1989. Local versus global variations in ocean ridge basalt composi-
tion: a reply. J. Geophys. Res. 94: 4241-4252.
Langmuir, C. H., R. D. Vocke, G. N. Hanson, and S. R. Hart, 1978. A general mixing equation with
applications to Icelandic basalts, Earth Planet. Sci. Lett., 37: 380-392.
Martin, H. 1986. Effect of steeper Archean geothermal gradient on geochemistry of subduction-zone
magmas. Geology. 14: 753-756.
McCulloch, M. T. 1987. Sm-Nd isotopic constraints on the evolution of Precambrian crust in the Aus-
tralian continent. In Proterozoic Lithosphere Evolution, Geodyn. Ser. Vol. 17, A. Kršner. ed., pp.
115-130. Washington: AGU.
McCulloch, M. T. and V. C. Bennett. 1994. Progressive growth of the Earth's continental crust and de-
pleted mantle: geochemical constraints. Geochim. Cosmochim. Acta. 58: 4717-4738.
McKenzie, D. P. and R. K. O'Nions. 1983. Mantle reservoirs and ocean island basalts. Nature. 301:
229-231.
Nyblade, A. A. and H. N. Pollack. 1993. A global analysis of heat flow from Precambrian terrains:
implications for the thermal structure of Archean and Proterozoic lithosphere. J. Geophys. Res. 98:
12207-12218.
O'Nions, R. K. and P. J. Hamilton. 1981. Isotope and trace element models of crustal evolution. P h i l .
Trans. Roy. Soc. Lond. A. 301: 473-487.
O'Nions, R. K. and D. P. McKenzie. 1988. Melting and continent generation. Earth Plan Sci Lett. 90:
449-456.
Patchett, P. J. 1980. Thermal effects of basalt on continental crust and crustal contamination of mag-
mas. Nature. 283, 559-561.
551 January 25, 1998

Patchett, P. J. and N. T. Arndt. 1986. Nd isotopes and tectonics of 1.9Ð1.7 Ga crustal genesis. E a r t h
Planet. Sci. Lett. 78: 329-338.
Perfit, M. R., D. A. Gust, A. E. Bence, R. J. Arculus and S. R. Taylor. 1980. Chemical characteristics of
island-arc basalts: implications for mantle sources. Chem. Geol. 30: 227-256.
Plank, T. and C. H. Langmuir. 1988. An evaluation of the global variations in the major element
chemistry of arc basalts. Earth Planet. Sci. Lett. 90, 349-370.
Plank, T. and C. H. Langmuir. 1993. Tracing trace elements from sediment input ot volcanic output a t
subduction zones. Nature. 362: 739-742.
Plank, T. and C. H. Langmuir. in press. The geochemical composition of subducting sediment and its
consequences for the crust and mantle. Chem. Geol.
Reymer, A. and G. Schubert. 1984. Phanerozoic addition rates to the cotinental crust and crustal
growth. Tectonics. 3: 63-77.
Reymer, A. and G. Schubert. 1986. Rapid growth of some major segments of continental crust. Geology.
14: 299-302.
Roeder, P. L. and R. F. Emslie. 1970. Olivine-liquid equilibrium. Contrib. Mineral. Petrol. 29: 275-289.
Rudnick, R. 1995. Making continental crust. Nature. 378: 571-578.
Rudnick, R. L. and D. M. Fountain. 1995. Nature and composition of the continental crust: a lower crus-
tal perspective. Rev. Geophys. 33: 267-309.
Shaw, D. M. 1967. An estimate of of the chemical composition of the Canadian Shield, Can. J. E a r t h
Sci. 4, 82.
Shaw, D. M., J. J. Gramer, M. D. Higgins and M. G. Truscott. 1986. Composition of the Candadian Pre-
cambrian shield and the continental crust of the earth. In The Nature of the Lower Continental
Crust, Vol. J. B. Dawson, C. D. A., J. Hall and K. H. Wedepohl. ed., pp. 275-282. Oxford: Blackwell
Scientific Publ.
Sinton, J. M. and R. Detrick. 1992. Mid-ocean ridge magma chambers, J. Geophys. Res. 97, 197-216.
Stolper, E. and S. Newman. 1994. the role of water in the petrogenesis of Mariana trough basalts.
Earth Planet. Sci. Lett. 121: 293-325.
Sun, S. S. 1980. Lead isotope study of young volcanic rocks from mid-ocean ridges, oceanic islands, and
island arcs, Philos. Trans. R. Soc. London, Ser. A, 297: 409-445.
Taylor, H. P. 1980. The effects of assimilation of country rocks by magmas on 18O/16O and 87Sr/86Sr sys-
tematics in igneous rocks. Earth Planet. Sci. Lett. 47, 243-254.
Taylor, S. R. 1967. The origin and growth of continents, Tectonophysics. 4, 17-34.
Taylor, S. R. and S. M. McLennan. 1985. The Continental Crust: its composition and evolution. Oxford:
Blackwell Scientific Publishers.
Taylor, S. R. and S. M. McLennan. 1995. The geochemical evolution of the continental crust. Rev. Geo-
phys. 33: 241-265.
Tera, F., L. Brown, J. Morris, I. S. Sacks, J. Klein, R. Klein and R. Middleton. 1986. sediment incorpora-
tion in island-arc magmas: inferences from 10Be, Geochim. Cosmochim. Acta. 50: 535-550.
Thorpe, R.S. (ed.). 1982. Orogenic Andesites and Related Rocks. Chichester: Wiley & Sons.
Veizer, J. and S. L. Jansen. 1979. Basement and sedimentary recycling and continental evolution. J.
Geol. 87: 341-370.
Vervoort, J. D., W. M. White, R. I. Thorpe and J. M. Franklin. 1993. Postmagmatic thermal activity
in the Abitibi Greenstone Belt, Noranda and Matagami Districts: evidence from whole rock Pb iso-
tope data. Econ. Geol. 88: 1598-1614.
Vitorello, I. and H. N. Pollack, On the variation of continental heat flow with age and the thermal
evolution of continents, J. Geophys. Res., 85, 983-995, 1980.
Weaver, B. L. and J. Tarney. 1984. Major and trace element composition of the continental litho-
sphere. In Structure and Evolution of the Continental Lithosphere, Physics and Chemistry of t h e
Earth, vol 15, Vol. H. N. Pollack and V. R. Murthy. ed., pp. 39-68. Oxford: Pergamon Press.
552 January 25, 1998

Wedepohl, K. H. 1995. The composition of the continental Data for Problem 12.1
crust. Geochim. Cosmochim. Acta. 59: 1217-1232. A B C
White, W. M. and B. DuprŽ, Sediment subduction and magma
Cs 1.07 0.015
genesis in the Lesser Antilles: Isotopic and trace element
constraints, J. Geophys. Res., 91, 5927-5941, 1986. Rb 34.00 17.96 1.33
White, W. M. and P. J. Patchett, 1984, Hf-Nd-Sr isotopes Ba 389.47 353.83 12.30
and incompatible element abundances in island arcs: impli- Th 7.80 14.50 0.29
cations for magma origins and crust-mantle evolution, U 1.89 4.99 0.11
Earth Planet. Sci. Lett., 67: 167-185. Nb 45.67 9.65 5.93
Problems Ta 2.51 0.45 0.35
1. The adjacent table shows trace element data from three K 25921.00 5130.00 1160.00
volcanic rocks, each erupted in a different tectonic environ- La 42.78 31.38 6.93
ment. Plot the data on a Òspider diagramÓ. Make the y-axis Ce 93.76 66.73 23.10
logarithmic and order the elements as in Figure 12.32. Pb 4.14 7.86 0.68
ÒNormalizeÓ the data to ÒPrimitive MantleÓ using the con- Pr 11.96 7.41 4.31
centrations given in Table 10.3. Based on your plot, identify
Sr 717.87 798.75 92.60
the tectonic setting in which each erupted.
Nd 48.20 28.21 22.80
2. In a sedimentary core, the surface sample had a 10Be con- Sm 10.36 5.29 7.89
centration of 100 × 106 atoms/g. Assuming steady state accu- Zr 300.91 121.92 227.00
mulation of 10Be at the surface, what concentration of 10B Hf 6.97 3.86 7.15
would you predict for a 10 million year old sample from Eu 3.06 1.45 2.50
lower in the same core? (HINT: use the decay constant given
Gd 8.90 4.83 10.90
in Table 7.5).
Tb 1.36 0.69
3. The following Sr and Nd isotope data from St. Lucia show Dy 6.37 3.51 13.00
that these volcanics are mixtures of two components. Assum- Y 30.34 19.73 69.10
ing that component A has 87Sr/86Sr, εNd, Sr, and Nd concentra- Ho 1.06 0.68 2.89
tions of 0.706, -0.5, 200 and 8.5 ppm respectively, and compo- Er 2.33 1.71 8.46
nent B has 87Sr/86Sr and εNd of 0.710 and -13.4, estimate the Sr
Tm 0.30 0.27
and Nd concentrations of component B. (HINT: Plot the d a t a
and superimpose a 2 component mixing model on the plot; a d - Yb 1.63 1.58 7.94
just the concentrations by trial and error until you find t h e Lu 0.22 0.23 1.19
best fit).
Sr and Nd Isotope Data
4. The adjacent table below shows from St. Lucia
concentrations of Sr and 87Sr/86Sr Sample 87
Sr/86Sr εN d Data for Problem 4
87
measured on river samples. These Sr, µg/l Sr/86Sr
STL220 0.70655 -3.02
samples were taken just below t h e SL124 0.70933 -11.76 45 0.71571
point where two major rivers join. 60 0.71305
SL125 0.70782 -8.74
Are data consistent with simple mix-
SL51 0.70845 -10.26 78 0.71162
ing between the waters of these two
SL52 0.70850 -10.09 100 0.71073
rivers? Why or why not.
SL121 0.70907 -11.26 123 0.71013
5. Plot oxygen and strontium isotope SL48 0.70885 -10.90 136 0.70969
ratios measured on lavas from SL42 0.70829 -9.71 161 0.70935
Galeras Volcano (Colombia) in t h e SL40 0.70837 -9.99 179 0.70909
following table, then decide
whether mixing curve they define a result of sediment subduction (Òsource contaminationÓ) or assimi-
lation (Òcrustal contaminationÓ)? Justify your answer.
553 January 25, 1998

6. Construct an AFC model of the data from Galeras Volcano

Data from Galeras Volcano
in the accompanying table. Assume that the 87Sr/86Sr of t h e
87
Sr/86Sr δ18O Sr, ppm
assimilant and original magma are 0.712 and 0.70375, respec-
0.70414 6.79 503 tively, the δ18O of the assimilant and original magma are
0.70409 6.65 504 +10.5 and +5.8, respectively, and Sr concentrations of the as-
0.70453 7.65 545 similant and original magma were 100 and 470 ppm respec-
0.70423 6.99 521 tively. Find values of R (ratio of mass assimilated to mass
0.70434 7.51 521 crystallized) and DSr (bulk partition coefficient for Sr) t h a t
0.70435 7.11 534 best fit the data (assume the partition coefficient for O is 1).
0.70419 6.72 514 (HINT: plot both 87Sr/ 86Sr vs. Sr and δ18O vs 87Sr/ 86Sr, and
adjust values of R and D until the line describing your m o d e l
passes through the data).
7. (a.) Assuming the Bulk Silicate Earth (BSE) has the Òprimitive mantleÓ Rb concentration given in
Table 10.3, what fraction of the Rb in the BSE is in the crust according to the Taylor and McLennan,
Rudnick and Fountain, and Wedepohl crustal compositions (Table 12.7)? (HINT: use the volumes o f
crust and mantle given in the Appendix).
(b.) Assume that the silicate Earth consists of only three reservoirs: continental crust, depleted man-
tle and primitive mantle, that the depleted mantle has a Rb concentration that is 10 times lower
than the average MORB concentration (Table 12.2), and that the primitive mantle has the Rb concen-
tration given in Table 10.3. What fraction of the mantle would consist of depleted mantle when you
use the Taylor and McLennan, Rudnick and Fountain, and Wedepohl crustal Rb concentrations (Table
12.7)?
8. Weathering and erosion remove Mg from the continents and transport it in dissolved form to t h e
oceans. Dissolved Mg is removed from the oceans by hydrothermal activity at mid-ocean ridges. As-
suming that the flux of Mg into the oceanic crust given in Chapter 15 has been constant through time
over the past 2.5 Ga and that all this Mg is derived from the continental crust, how much has t h e
crustal Mg concentration decreased as a result of this process? How might this result change if we
make other ÒreasonableÓ assumptions about the flux of Mg into the oceanic crust through time?
554 January 25, 1998

W - M - W H I T e G e o C H e M I S T R y Chapter 12: The Crust Chapter 12: Geochemistry of The Solid Earth II: The Crust

Uploaded by

Copyright:

Available Formats

W - M - W H I T e G e o C H e M I S T R y Chapter 12: The Crust Chapter 12: Geochemistry of The Solid Earth II: The Crust

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

W - M - W H I T e G e o C H e M I S T R y Chapter 12: The Crust Chapter 12: Geochemistry of The Solid Earth II: The Crust

Uploaded by

Copyright:

Available Formats

W. M.

Chapter 12: Geochemistry of the Solid Earth II: The

512 January 25, 1998

513 January 25, 1998

Ave Mg# = 52.7 that all primary mantle-derived magmas must be

250 > 80 mm/yr

50 3 Na8.0 =2.31 E EEE

Figure 12.2. Mg# in basalts from mid-ocean 2 JJ J

514 January 25, 1998

ply rates on slow spreading ZZA

515 January 25, 1998

516 January 25, 1998

nents. Thus intra-oceanic arcs may also contribute to continen-

517 January 25, 1998

magma stagnation and crystallization. Tholeiites occur in ex- 1

518 January 25, 1998

519 January 25, 1998

6 E EE E illustrated in Figure 12.10, using t h e

earth depleted. They gener- C SW Pacific

ally are enriched relative to

MORB in the alkalis and a l -

kaline earths. Significantly, 30 G

these are the elements most

enriched in sediments, and t h e

520 January 25, 1998

calc-alkaline basalt tion zones, which raises t h e

521 January 25, 1998

522 January 25, 1998

523 January 25, 1998

arc crust sediment

524 January 25, 1998

100 J solution, so the excess abundance of these ÒmobileÓ ele-

A long standing and key question is why magmas are

rine basalts from the Marianas. As mentioned

525 January 25, 1998

526 January 25, 1998

527 January 25, 1998

with having to decide

value is about +5.6. Island arc magmas derived δ O x = 1

than mantle values. However, high δ18O can also

528 January 25, 1998

where Cm is the concentration in the magma, C 0m

529 January 25, 1998

Composition of the Continental

9 The continental crust is extremely heterogene-

8 [Cm0 /Ca]Sr = 10 is a difficult one. Furthermore, only the upper

Table 12.3. Abundance of Igneous and Sedimentary Rocks in the Upper

530 January 25, 1998

Upper Continental Crust 'Shaw et al.'

al., 1986; Wedepohl, 1995). Studies of these

531 January 25, 1998

The advantages of us-

532 January 25, 1998

533 January 25, 1998

Granite/Rhyolite 3.9 16.0 3.6 4.1 0.9× 104 2.67 2.5

534 January 25, 1998

Heat Flow (mW/m 2)

535 January 25, 1998

536 January 25, 1998