diffusion data for silicate minerals, glasses, and …jbrady/papers/diffusionpaper_a.pdfdiffusion...

22
Diffusion Data for Silicate Minerals, Glasses, and Liquids John B. Brady 1. INTRODUCTION Diffusion is an integral part of many geologic pro- cesses and an increasing portion of the geologic literature is devoted to the measurement, estimation, and application of diffusion data. This compilation is intended to be a guide to recent experimentally- determined diffusion coefficients and, through the papers cited, to important older literature. To provide a context for the tables, a brief summary of the equa- tions required for a phenomenological (macroscopic) description of diffusion follows. Although the equa - tions for well-constrained experiments are relatively straightforward, the application of the resulting diffu- sion coefficients to complex geologic problems may not be straightforward. The reader is urged to read one or more diffusion texts [e.g., 99, 121, 32, 147, 13.51 before attempting to use the data presented here. 2. FORCES AND FLUXES Diffusion is the thermally-activated, relative move- ment (flux) of atoms or molecules that occurs in response to forces such as gradients in chemical potential or temperature. Diffusion is spontaneous and, therefore, must lead to a net decrease in free J. B. Brady, Department of Geology, Smith College, North- ampton, MA 01063 Mineral Physics and Crystallography A Handbook of Physical Constants AGU Reference Shelf 2 Copyright 1995 by the American Geophysical Union. energy. For example, the movement of An i moles of chemical component i from a region (II) of high chemical potential (/.Li”) to a region (I) of lower chemical potential (pi’) will cause the system Gibbs energy (G) to fall since and using the definition of pi [ 139, p. 1281, (2) Thus, a chemical potential gradient provides a therm - odynamic force for atom movement. On a macroscopic scale, linear equations appear to be adequate for relating each diffusive flux to the set of operative forces [137, p.451. The instantaneous, one-dimensional, isothermal diffusive flux JR (moles of i/m2s) of component i in a single-phase, n- component system with respect to reference frame R may be described by (3) where x (m) is distance and the n2 terms L: (moles of i/m.J.s) are “phenomenological diffusion coeffi- cients” [36]. Because each component of the system may move in response to a gradient in the chemical potential of any other component, the complexity of describing diffusion in multicomponent systems rises 269

Upload: duongbao

Post on 24-Mar-2018

217 views

Category:

Documents


2 download

TRANSCRIPT

Diffusion Data for Silicate Minerals, Glasses, and Liquids

John B. Brady

1. INTRODUCTION

Diffusion is an integral part of many geologic pro- cesses and an increasing portion of the geologic literature is devoted to the measurement, estimation, and application of diffusion data. This compilation is intended to be a guide to recent experimentally- determined diffusion coefficients and, through the papers cited, to important older literature. To provide a context for the tables, a brief summary of the equa- tions required for a phenomenological (macroscopic) description of diffusion follows. Although the equa - tions for well-constrained experiments are relatively straightforward, the application of the resulting diffu- sion coefficients to complex geologic problems may not be straightforward. The reader is urged to read one or more diffusion texts [e.g., 99, 121, 32, 147, 13.51 before attempting to use the data presented here.

2. FORCES AND FLUXES

Diffusion is the thermally-activated, relative move- ment (flux) of atoms or molecules that occurs in response to forces such as gradients in chemical potential or temperature. Diffusion is spontaneous and, therefore, must lead to a net decrease in free

J. B. Brady, Department of Geology, Smith College, North- ampton, MA 01063

Mineral Physics and Crystallography A Handbook of Physical Constants AGU Reference Shelf 2

Copyright 1995 by the American Geophysical Union.

energy. For example, the movement of An i moles of chemical component i from a region (II) of high chemical potential (/.Li”) to a region (I) of lower chemical potential (pi’) will cause the system Gibbs energy (G) to fall since

and using the definition of pi [ 139, p. 1281,

(2)

Thus, a chemical potential gradient provides a therm - odynamic force for atom movement.

On a macroscopic scale, linear equations appear to be adequate for relating each diffusive flux to the set of operative forces [137, p.451. The instantaneous, one-dimensional, isothermal diffusive flux JR (moles of i/m2s) of component i in a single-phase, n- component system with respect to reference frame R may be described by

(3)

where x (m) is distance and the n2 terms L: (moles of i/m.J.s) are “phenomenological diffusion coeffi- cients” [36]. Because each component of the system may move in response to a gradient in the chemical potential of any other component, the complexity of describing diffusion in multicomponent systems rises

269

n i

i = 1

n

∑ d µ i = 0

BRADY 271

Table 1. Commonly-used solutions to equation (9)

Boundary Conditions Solution

Thin-film Solution Ci + CO for Ix]> 0 as t + 0 c, 3.00 for x = 0 as t + 0

} (y)=&exp[&) where cx=lL”(Ci -CO)dx (10)

Semi-infinite Pair Solution

1

Ci =CO for x>O at t=O C, =Cl for x<O at t=O 1

(E) = ierfc[--&=j (11)

Finite Pair Solution Ci =Cl for O<x<h at t=O Ci =CO for h<x<L at t=O

(ilC,/~x), =0 for x=0&L

Finite Sheet - Fixed Surface Composition

cases such as (12) and (13), the fractional extent of completion of homogenization by diffusion is propor- tional to a. These fi relations provide impor- tant tests that experimental data must pass if diffusion is asserted as the rate-controlling process. They also provide simple approximations to the limits of diffusion when applying measured diffusion coefficients to specific problems [ 1471.

4. DIFFUSION COEFFICIENTS

In general, one must assume that Di is a function of Ci. Therefore, equations (9)-( 13) may be used with confidence only if Ci does not change appreciably during the experiment. This is accomplished either (a) by using a measurement technique (typically involving radioactive tracers) that can detect very small changes in Ci, or (b) by using diffusional exchange of stable isotopes of the same element that leave the element concentration unchanged. Approach (a) yields a “tracer diffusion coefficient” for the element that is specific to the bulk composition studied. Approach (b) yields a “self-diffusion coeffi - cient” for the isotopically doped element that is also specific to the bulk chemical composition. Both approaches generally ignore the opposite or exchange flux that must occur in dominantly ionic phases such

as silicate minerals, glasses, and liquids. If Ci does change significantly in the experiment

and Di is a function of Ci, observed compositional profiles might not match the shape of those predicted by (9)-(13). In such cases the Di calculated with these equations will be at best a compositional “average” and equation (8) should be considered. Experiments in which composition does change significantly are often termed “interdiffusion” or “chemical diffusion” experiments. A commonly-used analytical solution to (8), for the same boundary conditions as for equation (1 1), was obtained by Matano [122] using the Boltzmann [8] substitution (x=x/& : 1

Equation (14) can be evaluated numerically or graphi - tally from a plot of (Ci - CO)/(Cl - CO) versus x, where the point x = 0 (the Matano interface) is selected such that

(15)

272 DIFFUSION DATA

0.0 1.0

(4 s 0.8 u

s u

3

; 0.6

x

9 3 0.4

-1.5 1.0 g 0.2

-2.0 0.0 0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1,

X2 x/L

Figure 1. Graphical solutions to equation (9). (A) The “thin film” solution of equation (10) is shown with a=1 for various values of fi (labels on lines). Plotting ln(Ci - CO) as a function of x2 at any time yields a straight line of slope -l/(4Dit). (B) The “finite pair” solution of equation (12) is shown for h/L=O.4 and various values of Dt/L2 (labels on curves). Plotting (Ci -CO)/(Cl - CO) as a function of x/L normalizes all cases to a single dimensionless graph. The initial boundary between the two phases is marked by a dashed vertical line.

[32, p.230-2341. For binary, cation exchange between ionic crystals, the Matano interface is the original boundary between the two crystals. Diffu- sion experiments that follow the Boltzmann-Matano approach have the advantage of determining Di as a function of Ci and the disadvantage of risking the complications of multicomponent diffusion, for which neither (8) nor (14) is correct.

Darken [35] and Hartley and Crank [83] showed that for electrically neutral species, the binary (AtiB) interdiffusion coefficients (D,,-, ) obtained in a Boltzmann-Matano type experiment and the tracer dif- fusion coefficients ( Da and D*,) for the interdiffus ing species are related by

D AB-, = (N,Di + NADIR)

where N A is the mole fraction and YA the molar activ- ity coefficient of component A. Darken developed his analysis in response to the experiments of Smigelskas and Kirkendall [ 1481, who studied the interdiffusion of Cu and Zn between Cu metal and Cu7&&~~ brass. Smigelskas and Kirkendall observed that MO wires (inert markers) placed at the boundary between the Cu and brass moved in the direction of the brass during their experiments, indicating that more Zn atoms than Cu atoms crossed the boundary. Darken’s analysis showed that in the presence of inert markers the

independent fluxes and, therefore, the “tracer” diffusion coefficients of Cu and Zn can be determined in addition to the interdiffusion coefficient DCuZnm,.

A similar analysis for binary cation interdiffusion (AZ w BZ where AZ represents Al”Z;’ and BZ represents BzbZ,“) in appreciably ionic materials such as silicate minerals yields [5, 112, 113, lo]:

I[

1 +

(17)

which requires vacancy diffusion if a#b. If a=b, then ( 17) simplifies to

D (Daz)( Diz )

AB-, = (NA,D;Z + NBZD;Z) Ii I +

(18)

[ 12 I]. This expression permits interdiffusion coeffi- cients to be calculated from more-easily-measured tracer diffusion coefficients. For minerals that are not

BRADY 213

ideal solutions, the “thermodynamic factor” in brackets on the right side of (l6)-( 18) can significantly change the magnitude (and even the sign!) of interdiffusion coefficients from those expected for an idea1 solution [ 12, 251.

coefficients Dy to calculate multicomponent interdif- fusion coefficient matrices, and its inverse offer the most hope for diffusion analysis of multicomponent problems [ 1 15, 19, 201. It should be noted in using equations ( I6)-( 19) that the tracer diffusivities themselves may be functions of bulk composition.

5. MULTICOMPONENT DIFFUSION 6. EXPERIMENTAL DESIGN

Rarely is diffusion in geologic materials binary. Multicomponent diffusion presents the possibility that diffusive fluxes of one component may occur in response to factors not included in (7)-( 1 S), such as gradients in the chemical potentials of other compo- nents, coupling of diffusing species, etc. In these cases, equation (6) or (3) must be used. The off- diagonal (i#j) diffusion coefficients, Dg or Li, are unknown for most materials and most‘workers are forced to assume (at significant peril!) that they are zero. Garnet is the one mineral for which off- diagonal diffusivities are available [ 19, 201 and they were found to be relatively small.

One approach used by many is to treat diffusion that is one-dimensional in real space as if it were one- dimensional in composition space. This approach was formalized by Cooper and Varshneya [28, 301 who discuss diffusion in ternary glasses and present criteria to be satisfied to obtain “effective binary dif- fusion coefficients” from multicomponent diffusion experiments. In general, diffusion coefficients obtained with this procedure are functions of both composition and direction in composition space. Another approach, developed by Cullinan [33] and Gupta and ,Cooper [75, 291, is to diagonalize the diffusion coefficient matrix for (6) through an eigenvector analysis. Although some simplification in data presentation is achieved in this way, a matrix of diffusion coefficients must be determined before the analysis can proceed.

Lasaga [105, see also 1611 has generalized the relationship (17) to multicomponent minerals. The full expression is quite long, but simplifies to

(19)

if the solid solution is thermodynamically ideal. In (19) zi is the charge on cation i, 411 = 1 if i=j, and 6il = 0 if i#j. This approach, usmg tracer diffusion

A few additional considerations should be mentioned regarding the collection and application of diffusion data.

l In anisotropic crystals, diffusion may be a function of crystallographic orientation and the diffusion coefficient becomes a second rank tensor. The form of this tensor is constrained by point group symmetry and the tensor can be diagonalized by a proper choice of coordinate system [ 13 I].

l Vacancies, dislocations, and other crystal defects may have profound effects on diffusion rates [e.g., 121, 1641. Therefore, it is essential that the mineral being studied is well-characterized. If the mineral, liquid, or glass contains a multivalent element like Fe, then an equilibrium oxygen fugacity should be controlled or measured because of its effect on the vacancy concentration.

l Diffusion in rocks or other polycrystalline materials may occur rapidly along grain boundaries, crystal interfaces, or surfaces and will not necessarily record intracrystalline diffusivities. Experiments involving grain-boundary diffusion [e.g., 102, 103, I 1, 47, 49, 95, 154, 13, 1581 are beyond the scope of this summary. If polycrystalline materials are used in experiments to measure “intrinsic” or “vol- ume” diffusion coefficients, then the contributions of “extrinsic” grain-boundary diffusion must be shown to be negligible.

l Water can have a major effect on diffusion in many geologic materials, even in small quantities [e.g., 73, 150, 46, 167, 711. If water is present, water fugacities are an essential part of the experimental data set.

l Because many kinetic experiments cannot be “reversed,” every effort should be made to demonstrate that diffusion is the rate-controlling process. Important tests include the m relations noted in equations (lo)-( 13) and the “zero time experiment.” The @t test can be used if data are gathered at the same physical conditions for at least two times, preferably differing by a factor of four or more. The often-overlooked “zero time experiment,”

274 DIFFUSION DATA

which duplicates the sample preparation, heating to the temperature of the experiment, quenching (after “zero time” at temperature), and data analysis of the other diffusion experiments, commonly reveals sources of systematic errors.

These and other important features of diffusion experiments are described in Ryerson [141].

7. DATA TABLE

Due to space limitations, the data included in this compilation have been restricted to comparatively recent experimental measurements of diffusion coeffi- cients for silicate minerals, glasses, and liquids. Many good, older data have been left out, but can be found by following the trail of references in the recent papers that are included. Older data may also be found by consulting the compilations of Freer [56, 571, Hofmann [86], Askill [2], and Harrop [81]. Some good data also exist for important non-silicate minerals such as apatite [ 157, 23, 451 and magnetite- titanomagnetite [65, 136,591, but the silicate minerals offered the clearest boundary for this paper.

Almost all of the data listed are for tracer or self- diffusion. Interdiffusion (chemical diffusion) exper- iments involving minerals are not numerous and interdiffusion data sets do not lend themselves to compact presentation. Interdiffusion data for silicate liquids and glasses are more abundant, but are not in- cluded [see 156, 3, 41. Neither are non-isothermal (Soret) diffusion data [see 1081. Diffusion data listed for silicate glasses and liquids have been further restricted to bulk compositions that may be classified as either basalt or rhyolite.

Because diffusion is thermally activated, coefficients for diffusion by a single mechanism at different temperatures may be described by an Arrhenius equation

D = D, exp (20)

and fit by a straight line on a graph of log D (m2/s) as a function of l/T (K-i) [147, Chap. 21. Log D, (m2/s) is the intercept of the line on the log D axis (l/T = 0). AHl(2.303.R) is the slope of the line where R is the gas constant (8.3143 J/mole.K) and

AH (J/mole) is the “activation energy.” D, and AH have significance in the atomic theory of diffusion [see 32, 1211 and may be related for groups of similar materials [162, 82, 1061.

Diffusion data are listed in terms of AH and log D, and their uncertainties. In many cases data were converted to the units &I, m2/s) and form (log D,) of this compilation. Logarithms are listed to 3 decimal places for accurate conversion, even though the original data may not warrant such precision. No attempt was made to reevaluate the data, the fit, or the uncertainties given (or omitted) in the original papers. Also listed are the conditions of the experiment as appropriate including the temperature range, pressure range, oxygen fugacity, and sample geometry. Extrapolation of the data using (20) to conditions outside of the experimental range is not advisable. However, for ease of comparison log D is listed for a uniform temperature of 800°C (1200°C for the glasses and liquids), even though this temperature may be outside of the experimental range. Use these tabulated log D numbers with caution.

Finally, a sample “closure temperature” (T,) has been calculated for the silicate mineral diffusion data. The closure temperature given is a solution of the Dodson [38] equation

for a sphere of radius a=O.l mm and a cooling rate (dT/dt) of 5 K/Ma. The example closure temperatures were included because of the importance of closure temperatures in the application of diffusion data to petrologic [107] and geochronologic [123, Chap. 51 problems. Note that other closure temperatures would be calculated for different crystal sizes, cooling rates, and boundary conditions.

Acknowledgments. This paper was improved by help- ful comments on earlier versions of the manuscript by D. J. Cherniak, R. A. Cooper, R. Freer, J. Ganguly, B. J. Giletti, S. J. Kozak, E. B. Watson, and an anonymous reviewer. N. Vondell provided considerable help with the library work.

Mineral, Glass, Orienta- Diffusing Temperature P 0, 4 log D, (or D) log D ‘IT,” Experiment/Comments Ref. or Liquid tion Component Range (“C) (MPa) (MPa) (kJ/mole) (m*/s) 8oo”c (“C)

m

E 4

Y WI

u Mineral, Glass, Orienta- Diffusing Temperature P 02 ma log D, (or D) log D ‘IT,” Experiment/Comments Ref. 7

or Liquid tion Component Range (“C) (MPa) (MPa) (kJ/mole) (m2/s) 800°C (“C)

[l@l ? 2

[loll

Mineral, Glass, Orienta- Diffusing Temperature P O2 4 log Do (or D) log D “T,” Experiment/Comments Ref. or Liquid tion Component Range (“C) (MPa) (MPa) (kJ/mole) (m*/s) 800°C (“C)

nepheline powder

biotite powder (see paper for camp.) phlogopite powder (see paper for camp.) phlogopite powder (Ann4) (see paper for camp.) biotite (An& powder (see paper for camp. reference) biotite powder

chlorite (sheridamte) (see paper) mucovite

powder

powder

muscovite (see powder paper for camp.)

‘80

180

4oAr

Sr

“wet” NO 650-800 100,20 uB/

2000 NO “wet”

650-800 100 UB/ “wet” NO

750-900 100 UBo “wet” NO

800 200 uB/ “wet” NO

c

r

I 104.6(+10.5) -12.229 -17.3 199 Exchange with CO*, spherical [27]

model, bulk analysis

-9.041

-7.854

-16.0

-16.4

233 Exchange with ‘*O-enriched [541 water, bulk analysis, cylin- drical model, ion probe too

305 Exchange with 180-enriched [541 water, bulk analysis, cylin-

I ’ i

1’2”‘:’ 88

171.7 -5.21 -13.6

121.3 -7.98 -13.9

drical model, ion probe too 373 Degassing into water, bulk t621

analysis, cylindrical model, 671 (AV,zO m3/mole)

281 Degassing into water, bulk [771 analysis, cylindrical model (AV =1.4~10-~ m3/mole) Exchange with 41KC1 solu- F371

t

163 (+21) -8.114 -16.1

163 (f21) -11.699(f1.204) -19.6

172 (+25) -1l.OOO(f1.322) -19.4

239 (k8) -7.523 (f0.452) -19.1

268.2 (k7.1) -5.620 (f0.506) -18.7

D=lxlO-*’

t

t L

c

t

t

c

L

on, bulk analysis & ion

m

F 4

Y 4

ii

Mineral, Glass, Orienta- Diffusing Temperature P O2 4 log Do (or D) log D ‘IT,” Experiment/Comments Ref. u or Liquid tion Component Range (“C) (MPa) (MPa) (kllmole) (m2/s) 800°C (“C) ii

Sermanite (syn) lc-axis 1x0 02MgSi207) lkermanite (svn) I c-axis ‘80 (CazMgSizb;) 1 Sermanite (syn) ///c-axis / 45Ca UhMgSi207) Sermanite (syn) /c-axis 54Mn 02MgSi207) lkermanite (syn) /c-axis 59Fe 02MgSi207) &ermanite (syn) /c-axis 6oCo GbMgSbO,)

500-800 0.1 -12.2- -17.9 UB

1180 “wet” 0.1 w, 1

=I@‘4

1000-1300 0.1 co2

800-1300 0.1 (CO)/ co,:

1100-1300 0.1 N2

1100-1300 0.1 N2

1100-1300 0.1 N2

1100-1300 N2

I

i

robe, Rutherford backscatter,

m, sectioning, Ion probe, Rutherford backscatter,

I

I

I atmosphere, ion probe profile

Mineral, Glass, Orienta- Diffusing Temperature P 02 4 log D, (or D) log D “T,” Experiment/Comments Ref. or Liquid tion Component Range (“C) (MPa) (MPa) (kJ/mole) (m*/s) 800°C (“C)

Mineral, Glass, Orienta- Diffusing Temperature P O2 =a log Do (or D) log D “T,” Experiment/Comments Ref. or Liquid

g tion Component Range (“C) (MPa) (MPa) (kJ/mole) (m2/s) 800°C (“C)

2

olivine sintered 59Fe 1130 0.1 lo-to- n.d. n.d. log D59i+ = -10.143 + 0.2 [851 Wd%o-Faloo~ powder 10-12 , act + 2.705 [Fe/(Fe+Mg)]

olivine //c-axis Ni 1149-1234 evac. UB 193 (f10) -8.959 (f0.36) -18.4 / 404 Thin Ni film, microprobe G'61 (Fo93.7Fa6.3) tube profile, anisotropy found olivine (Fo,,Fas> //c-axis Fe-Mg 1125-1200 0.1 lo-l3 243 -5.759 -17.6 432 Microprobe profile, D varies [I41 fayalite powder (couple) interdiffusion w/direction, composition, fo2 olivine (Fa,,Te, )- //c-axis Fe-Mg 900-1100 evac. UB 208.5(+18.8) Do= 1.5(*0.3) Microprobe profile, D also ~1241 divine Fo91Fa9) (couple) interdiffusion tube +9.1 x [Mg/ x10-4) - 1.1x varies with direction and P

Wg+Fe)l [Mg/(Mg+Fe)] (AV =5,5~10-~m~/mole) a olivine(MgzSi04) //c-axis Co-Mg 1150-1300 0.1 air 196 -8.690 -18.2 ! 403 Microprobe profile, D’s [l-m - (Co2 SiO.+) (couple) interdiffusion 1300- 1400 526 2.288 -23.3 78 1 extrapolated to pure Fo forsterite - //c-axis Ni-Mg 1200-1450 0.1 air 414 to 444 -1.652 -21.8 702 Microprobe profile, other [1261 liebensbergite (couple) interdiffusion interdiffusion coefficients [128]

grossularite C&9 isotropic ‘80 1050 800 UB n.d. D=2.5~10-~~ Exchange with t80-enriched [581 Fe.1 Ah$hPd 850 200 nd. D=4.8xlO-2’ water, ion probe profiles almandine isotropic ‘80 8001000 100 UB 301 (+46) ,-8.222 (*0.740) -22.9 725 Exchange with ‘*O-enriched I311 (Al67SP2xAns b2) water, ion probe profiles ww isotropic 25Mg 750-900 200 uB/ 239 (*16) -8.009 -19.6 5 14 Thin 25Mg0 film, ion probe [341 E’m%GrdJrI ) MH profiles AlmsoPYPzo - isotropic Fe 1300-1480 2900- UBI 275.43 -7.194 -20.6 588 ! / Calculated from interdiffusion VOI Spess94Almg (couple) 4300 GrPC (rt36.49) experiments using model,

(AV,=5.6~10-~ m3/mole) AlmzoPYPzo - isotropic Mg 13001480 2900- UB/ 284.52 -6.959 -20.8 604 Calculated from interdiffusion WI Spessa4Alm6 (couple) 4300 GrPC (T37.55) experiments, (AV, =

5.6(f2.9)x10e6 m3/mole) AlmsoPYPzo - isotropic Mn 1300-1480 2900- UB/ 253.44 -7.292 -19.6 526 Calculated from interdiffusion m Spessa4Alm6 (couple) 4300 GrPC (k37.19) experiments, (AV, =

5.3(+3.0)x10-6 m3/mole) almandine isotropic 86Sr 800- 1000 100 UB 205 (+17) -12.000(&0.602) -22.0 616 Exchange with %r water [311 b’%SPd% pY2 ) solution, ion probe profiles almandine isotropic 145w 800-1000 100 UB 184 (+29) -12.523(*0.602) -21.5 ! 562 Exchange with 145Nd water 1311 b%SP2dnd’Y2) solution, ion probe profiles wope powder tslSm 1300-1500 3000 UB/ 140 -11.585 -18.4 321 Exchange with silicate melt, WI

GTPC autoradiography, “sphere” almandine isotropic 167 ET 800- 1000 100 UB 230 (f38) -10.301(?0.763) -21.5 605 Exchange with 167Er water [311 @l&Pd@Y2) solution, ion probe profiles AlmsoPYPzo - isotropic Fe-Mn 1300- 1480 4000 UBI 224.3 - 10.086 -21.0 57 1 Microprobe profiles, model [411 Spess94Alm6 (couple) interdiffusion GrPC (k20.5) fits of alm-rich composition

9 j!

Mineral, Glass, Orienta- or Liquid tion

Diffusing Temperature P O2 Component Range (“C) (MPa) (MPa)

=a (kJ/mole)

log D, (or D) log D “Tc” Experiment/Comments Ref. (m2/s) 1200°C (“C)

Zircon dissolution,

Mineral, Glass, Orienta- Diffusing Temperature P 0, aa log Do (or D) log D “Tc” Experiment/Comments Ref. or Liquid tion Component Range (“C) (MPa) (MPa) @J/mole) (m*/s) 1200°C (“C)

rhyolite glass and Ce 875-1100 0.1 air 490.4 2.72 (f0.99) -14.7 Thin film, serial sectioning by [91] (dehydrated) NM melt (f23.9) etching, counting surface rhyolite glass and Eu 700-1050 0.1 air 288.7 (k5.0) -3.11 (0.22) -13.3 Thin film, serial sectioning by [91] (dehydrated) NM melt etching, counting surface rhyolite (obsidian) melt 6% LREE 1000-1400 800 UB/ 251.5(*42.3) -4.638 (k1.436) -13.6 Monazite dissolution, [1381 Lake County, OR water GrPC microprobe profiles rhyolite (obsidian) melt 1% LREE 1000-1400 800 UB/ 510.9(+59.0) 3.362 (f0.629) -14.8 , Monazite dissolution, [1381 Lake County, OR water GrPC microprobe profiles rhyolite (obsidian) glass Hz0 400-850 0.1 Nz 103@5) -14.59 (f1.59) -18.2 Dehydration in Nz, FTIR [1661 Mono Craters profile, equilibrium model rhyolite (obsidian) glass H20 510-980 0.1 air 46.48 -10.90 (+0.56) -12.5 Dehydration in air, bulk 1931 New Mexico (fl 1.40) weight loss, low water rhyolite (obsidian) melt with 14CGz 800-1100 1000 UBI 75(?21) -7.187 -9.9 Thin film of Na2 14C03, B- 11551 Lake County, OR 8% water GrPC track profiles of cross section

basalt melt 6Li 1300-1400 0.1 air 115.5 -5.125 -9.2 Thin film of 6LiC1, ion probe [116] (alkali) profile of cross section basalt melt 0 1160-1360 0.1 IW to 215.9 -2.439 -10.1 Oxidation/reduction of bead, [1601 Goose Island co2 (f13.4) thermo-gravimetric balance basalt (alkali melt 0 1280-1400 400 UB/ 293(+29) -0.790 (k2.5 1) -11.2 Reduction by graphite, bulk ]401 olivine) BC GrPC Fe0 analysis by titration basalt (alkali melt 0 1280-1450 1200 UBI 360(+_25) 1.450 (t0.081) -11.3 Reduction by graphite, bulk ]401 olivine) BC GrPC Fe0 analysis by titration basalt (alkali melt 0 1350-1450 2000 uB/ 297 (+59) -0.770 (k1.87) -11.3 Reduction by graphite, bulk ]401 olivine) BC GrPC Fe0 analysis by titration basalt (tholeiite) melt 0 1300-1450 1200 uB/ 213 (k17) -3.010 (YkO.59) -10.6 Reduction by graphite, bulk [401 1921 Kilauea GrPC Fe0 analysis by titration basalt (FeTi) melt ‘80 1320-1500 0.1 Co2 251 (+29) -2.854 -11.8 Exchange with tsO-selected [I51 Galapagos ‘32 02 gas, bulk analysis of sphere “basalt” (Fe-free) melt Ar 1300-1450 looo- 113.2 (k7.5) -6.140 (kO.068) -10.2 [ Method not described r391’ (synthetic) 3000 basalt (“alkali”) melt 24Na 1300-1400 0.1 air 163 (+13) -4.02 (f0.46) -9.8 Thin film, serial sectioning by [116] Tenerife grinding, counting surface basalt (tholeiite) melt 45Ca 1260-1440 0.1 air 184.1 -4.272 -10.8 Thin film, B-track profiles of [891 1921 Kilauea cross-section on film “basalt” (Fe-free) melt 45Ca ilOO-14cKl 0.1 UB 106.3 -6.301 -10.1 Thin film of 45CaCl, B-track [152] ’ (synthetic) profiles of cross section “basalt” (Fe-free) melt 45Ca 1100-1400 1000 UB 141.0 -5.284 -10.3 Thin film of 45CaCl, B-track [152] 1 (synthetic) --I “basalt” (Fe-free) melt 45Ca 1100-1400 2000 UB 208.4 -3.211 -10.6 Thin film of 45CaCl, B-track t15211 (synthetic) profiles of cross section I

Mineral, Glass, Orienta- or Liquid tion

Diffusing Temperature P O2 Component Range (“C) (MPa) (MPa)

fia (kJ/mole)

log D, (or D) log D “Tc” Experiment/Comments Ref. (m2/s) 1200°C (“C)

T

basalt (“alkali”)

1921 Kilauea 1 basalt (“alkali”) 1

melt

melt

Tenerife basalt (“alkali”) melt Tenerife basalt (“alkali”) melt Tenerife basalt (tholeiite) melt 1921 Kilauea basalt (“alkali”) melt Tenerife basalt (tholeiite) melt

46Sc 1300-1400

54Mn

sSSr

1300-1400

59Fe

1 1300-1400

1300-1400

6OCo 1260-1440

6OCo 1300-1400

*%r 1260-1440

Tenerife basalt (tholeiite) melt 133& 1260-1440 1921 Kilauea basalt (“alkali”) melt 133& 1300-1400 Tenerife basalt (“alkali”) melt ‘Ts 1300-1400 Tenerife basalt (tholeiite) melt ‘=Eu, 1320-1440

Key to oxygen fugacity or atmosphere abbreviations in data table

0.1

0.1

0.1

0.1

0.1

0.1

0.1

0.1

0.1

0.1

0.1

air

air

air

air

air

air

air

air

air

air

air

ul3/ GrPC

CQ- GM - IW - h4H- N2 - NO -

Pure carbon dioxide atmosphere Graphite-methane buffer Iron-wustite buffer Magnetite-hematite buffer Pure nitrogen atmosphere Nickel-nickel oxide Pure oxygen atmosphere Quartz-fayalite-magnetite buffer

UB - Unbuffered oxygen fugacity UBlGrPC Unbuffered, but fo2 limited by the graphite-bearing piston-cylinder assembly UB/NO Unbuffered, but near nickel-nickel oxide due to the cold seal pressure vessel

hin film, serial sectionin

profiles of cross section,

284 DIFFUSION DATA

REFERENCES

I. Ashworth, J. R., Birdi, J. J., and Emmett, T. F., Diffusion in coronas around clinopyroxene: modelling with local equilib- rium and steady state, and c non- steady-state modification to account for zoned actinolite- hornblende, Contrib. Mineral. Petrol., 109, 307-325, 1992.

2. Askill, J., Tracer Dtffusion Data for Metals, Alloys, and Simple Oxides, 97 pp., Plenum Press, New York, 1970.

3. Baker, D. R., Interdiffusion of hydrous dacitic and rhyolitic melts and the efficacy of rhyolite contamination of dacitic enclaves, Contrib. Mineral. Petrol., 106, 462-473, 1991.

4. Baker, D. R., The effect of F and Cl on the interdiffusion of peralkaline intermediate and silicic melts, Am. Mineral., 78, 316-324, 1993.

5. Barrer, R. M., Bartholomew, R. F., and Rees, L. V. C., Ion ex- change in porous crystals. Part II. The relationship between self- and exchange-diffusion coefficients, J. Phys. Chem. Solids, 24, 309-317, 1963.

6. Barrer, Richard M., Diffusion in and through Solids, 464 pp., Cambridge University Press, Cambridge, England, 195 1.

7. Besancon, J. R., Rate of dis- ordering in orthopyroxenes, Am. Mineral., 66, 965-973, 1981.

8. Boltzmann, L., Zur Integration der Diffusionsgleichung bei variabeln Diffusion-coefficienten Ann. Physik Chem., 53, 959- 964, 1894.

9. Brabander, D. J., and Giletti, B. J., Strontium diffusion kinetics in amphiboles (abstract), AGU I992 Full Meeting, supplement to Eos Trans. AGU, 73, 64 1, 1992.

IO.

I I.

12.

13.

14.

15.

16.

17.

18.

19.

Brady, J. B., Reference frames and diffusion coefficients, Am. J. Sci., 275, 954-983, 1975. Brady, J. B., lntergranular diffu- sion in metamorphic rocks, Am. J. Sci., 283-A, I8 l-200, 1983. Brady, J. B., and McCallister, R. H., Diffusion data for clino- pyroxenes from homogenization and self-diffusion experiments, Am. Mineral., 68, 95- 105, 1983. Brenan, J. M., Diffusion of chlorine in fluid-bearing quartzite: Effects of fluid composition and total porosity, Contrib. Mineral. Petrol., 115, 215-224, 1993. Buening, D. K., and Buseck, P. R., Fe-Mg lattice diffusion in olivine, J. Geophys. Res., 78, 6852-6862, 1973. Canil, D., and Muehlenbachs, K., Oxygen diffusion in an Fe- rich basalt melt, Geochim. Cosmochim. Acta, 54, 2947- 2951, 1990. Carlson, W. D., and Johnson, C. D., Coronal reaction textures in garnet amphibolites of the Llano Uplift, Am. Mineral., 76, 756-772, 1991. Carslaw, H. S., Jaeger, J. C., Conduction of Heat in Solids, 510 pp., Oxford, Clarendon Press, 1959. Chakraborty, S., Dingwell, D., and Chaussidon, M., Chemical diffusivity of boron in melts of haplogranitic composition, Geo- chim. Cosmochim. Actu, 57, 1741-1751, 1993. Chakraborty, S., and Ganguly, J., Compositional zoning and cation diffusion in garnets, in Diffusion, Atomic Ordering, and Mass Transport, edited by Ganguly, J., pp. 120-175, Springer-Verlag, New York,

1991. 20. Chakraborty, S., and Ganguly,

J., Cation diffusion in alumino- silicate garnets: experimental determination in spessartine- almandine diffusion couples, evaluation of effective binary diffusion coefficients, and applications, Contrib. Minerul. Petrol., 111, 74-86, 1992.

21. Cherniak, D. J., Diffusion of Pb in feldspar measured by Rutherford backscatteing spectroscopy, AGU 1992 Full Meeting, supplement to Eos Trans. AGU, 73, 64 I , 1992.

22. Cherniak, D., Lead diffusion in titanite and preliminary results on the effects of radiation damage on Pb transport, Chem. Geol., 110, 177-194, 1993.

23. Cherniak, D. J., Lanford, W. A., and Ryerson, R. J., Lead diffusion in apatite and zircon using ion implantation and Rutherford backscattering tech- niques, Geochim. Cosmochim. Acta, 5.5, 1663-1673, 1991.

24. Cherniak, D. J., and Watson, E. B., A study of strontium diffusion in K-feldspar, Na-K feldspar and anorthite using Rutherford backscattering spec- troscopy, Earth Planet. Sci. Lett., 113, 41 l-425, 1992.

25. Christoffersen, R., Yund, R. A., and Tullis, J., Inter- diffusion of K and Na in alkali feldspars, Am. Mineral., 68, 1126-l 133, 1983.

26. Clark, A. M., and Long, J. V. P., The anisotropic diffusion of nickel in olivine, in Thomas Graham Memorial Symposium on Diffusion Processes, edited by Sherwood, J. N., Chadwick, A. V., Muir, W. M., and Swinton, F. L., pp. 5 I l-521, Gordon and Breach, New York, 1971.

BRADY 285

27. Connolly, C., and Muehlenbachs, K., Contrasting oxygen diffusion in nepheline, diopside and other silicates and their relevance to isotopic systematics in meteorites, Geochim. Cosmochim. Acta, 52, 15851591, 1988.

28. Cooper, A. R., The use and limitations of the concept of an effective binary diffusion coefficient for multi-component diffusion, in Mass Transport in Oxides, edited by Wachtman, J. B., Jr., and Franklin, A. D., pp. 79-84, National Bureau of Standards, Washington, 1968.

29. Cooper, A. R., Jr. Vector space treatment of multicomponent diffusion, in Geochemical Transport and Kinetics, edited by Hofmann, A. W., Giletti, B. J., Yoder, H. S., Jr., and Yund, R. A., pp. 15-30, Carnegie Institution of Washington, Washington, 1974.

30. Cooper, A. R., and Varshneya, A. K., Dii’fusion in the system &O-SrO-SiO,, effective binary diffusion coefficients, A m. Ceram. Sot. .I., 51, 103-106, 1968.

3 1. Coughlan, R. A. N., Studies in diffusional transport: grain boundary transport of oxygen in feldspars, strontium and the REE’s in garnet, and thermal histories of granitic intrusions in south-central Maine using oxygen isotopes, Ph.D. thesis, Brown University, Providence, Rhode Island, 1990.

32. Crank, J., The Mathematics of Diffusion, 414 pp., Clarendon Press, Oxford, 1975.

33. Cullinan, H. T. , Jr., Analysis of the flux equations of multi- component diffusion, Ind. Engr. Chem. Fund., 4, 133-139, 1965.

34. Cygan, R. T., and Lasaga, A.

35.

36.

37.

38.

39.

40.

41.

42.

43.

C., Self-diffusion of magnesium in garnet at 750-900°C Am. J. Sci., 285, 328-350, 1985. Darken, L. S., Diffusion, mobility and their interrelation through free energy in binary metallic systems, Am. Inst. Min. Metal. Engrs. Trans., 175, 184-201, 1948. de Groot, S. R., Mazur, P., Non-equilibrium Thermodynam- ics, 510 pp., North Holland, Amsterdam, 1962. Dennis, P. F., Oxygen self- diffusion in quartz under hydro- thermal conditions, J. Geophys. Res., 89, 4047-4057, 1984. Dodson, M. H., Closure temperature in cooling geo- chronological and petrological problems, Contrib. Mineral. Petrol., 40, 259-274, 1973. Draper, D. S., and Carroll, M. R., Diffusivity of Ar in haplobasaltic liquid at 10 to 30 kbar (abstract), AGU 1992 Fall Meeting, supplement to Eos Trans. AGU, 73, 642, 1992. Dunn, T., Oxygen chemical diffusion in three basaltic liquids at elevated temperatures and pressures, Geochim. Cosmo- chim. Acta, 47, 1923- 1930, 1983. Elphick, S. C., Ganguly, J., and Loomis, T. P., Experi- mental determination of cation diffusivities in aluminosilicate garnets, I. Experimental methods and interdiffusion data, Contrib. Mineral. Petrol., 90, 36-44, 1985. Elphick, S. C., Graham, C. M., and Dennis, P. F., An ion microprobe study of anhydrous oxygen diffusion in anorthite: a comparison with hydrothermal data and some geological implications, Contrib. tI4ineral. Petrol., 100, 490-495, 1988. Farver, J. R., Oxygen self-

44.

45.

46.

47.

48.

49.

50.

51.

52.

diffusion in diopside with application to cooling rate determinations, Earth Planet. Sci. Lett.. 92, 386-396, 1989. Farver, J. R., and Giletti, B. J., Oxygen diffusion in amphi- boles, Geochim. Cosmochim. Actu, 49, 1403-1411, 1985. Farver, J. R., and Giletti, B. J., Oxygen and strontium diffusion kinetics in apatite and potential applications to thermal history determinations, Geochim. Cos- mochim. Acta, 53, 162 l- I63 1, 1989. Farver, J. R., and Yund, R. A., The effect of hydrogen, oxygen, and water fugacity on oxygen diffusion in alkali feldspar, Geochim. Cosmochim. Acta, 54, 2953-2964, 1990. Farver, J. R., and Yund, R. A., Measurement of oxygen grain boundary diffusion in natural, fine-grained quartz aggregates, Geochim. Cosmochim. Acta, 55, 1597-1607, 1991a. Farver, J. R., and Yund, R. A., Oxygen diffusion in quartz: Dependence on temperature and water fugacity, Chem. Geol., 90, 55-70, 1991b. Farver, J. R., and Yund, R. A., Oxygen diffusion in a fine- grained quartz aggregate with wetted and nonwetted micro- structures, J. Geophys. Res., 97, 14017-14029, 1992. Fick, A., Uber diffusion, Ann. Physik Chem., 94, 59, 1855. Fisher, G. W., Nonequilibrium thermodynamics as a model for diffusion-controlled metamor- phic processes, Am. J. Sci., 273, 897-924, 1973. Foland, K. A., 40Ar diffusion in homogeneous orthoclase and an interpretation of Ar diffusion in K-feldspars, Geochim. Cosmo- chim. Acta, 38, 15 1- 166, 1974a.

286 DIFFUSION DATA

53.

54.

55.

56.

57.

58.

59.

60.

61 GCrar,d. O., and Jaoul, O.,

Foland, K. A., Alkali diffusion in orthoclase, in Geochemical Trcmsport and Kinetics, edited by Hofmann, A. W., Giletti, B. J., Yoder, H. S., Jr., and Yund, R. A., pp. 77-98, Carnegie Institution of Washington, Washington, 1974b. Fortier, S. M., and Giletti, B. J., Volume self-diffusion of oxygen in biotite, muscovite, and phlogopite micas, Geochim. Cosmochim. Acta, 55, 1319- 1330, 1991. Foster, C. T., A thermodynamic model of mineral segregations in the lower sillimanite zone near Rangeley, Maine, Am. Mineral,, 66, 260-277, I98 I. Freer, R., Self-diffusion and impurity diffusion in oxides, J. Materials Sci., 15, 803-824, 1980. Freer, R., Diffusion in silicate minerals and glasses: a data digest and guide to the literature, Contrib. Mineral. Petrol., 76, 440-454, I98 I. Freer, R., and Dennis, P. F., Oxygen diffusion studies. I. A preliminary ion microprobe investigation of oxygen diffusion in some rock-forming minerals, Mineral. Mug., 4.5, 179-192, 1982. Freer, R., and Hauptman, Z., An experimental study of mag- netite-titanomagnetite interdiffu- sion, Phys. Earth Planet. Inter., 16, 223-23 1, 1978. Ganguly, J., and Tazzoli, V., Fe+*-Mg interdiffusion in ortho- pyroxene: constraints from cation ordering and structural data and implications for cooling rates of meteorites (abstract), in Lunar and Planet. Science XXlV, pp. 5 17-5 18, Lunar and Planet. Institute, Houston, 1992.

62

63.

64.

65.

66.

67.

68.

69.

70.

Oxygen diffusion in San Carlos olivine, J. Geophys. Res., 94, 4119-4128, 1989. Giletti, B. J., Studies in diffusion 1: argon in phlogopite mica, in Geochemical Transport and Kinetics, edited by Hofmann, A. W., Giletti, B. J., Yoder, H. S., Jr., and Yund, R. A., pp. 107-l 15, Carnegie Institution of Washington, Washington, 1974. Giletti, B. J., Rb and Sr diffusion in alkali feldspars, with implications for cooling histories of rocks, Geochim. Cosmochim. Acta, 5.5, I33 I - 1343, 1991. Giletti, B. J., Diffusion kinetics of Sr in plagioclase (abstract), AGU 1992 Spring Meeting, supplement to Eos Trans. AGU, 73, 373, 1992, Giletti, B. J., and Hess, K. C., Oxygen diffusion in magnetite, Earth Planet. Sci. Lett., 89, 115-122, 1988. Giletti, B. J., Semet, M. P., and Yund, R. A., Studies in diffusion -- III. Oxygen in feldspars: an ion microprobe determination, Geochim. Cosmochim. Acta, 42, 45-57, 1978. Giletti, B. J., and Tullis, J., Studies in diffusion, IV. Pressure dependence of Ar Diffusion in Phlogopite mica, Earth Planet. Sci. Lett., 3.5, 180-183, 1977. Giletti, B. J., and Yund, R. A., Oxygen diffusion in quartz, J. Geophys. Res., 89, 4039-4046, 1984. Giletti, B. J., Yund, R. A., and Semet, M., Silicon diffusion in quartz, Geol. Sot. Am. Abstr. Prog., 8, 883-884, 1976. Graham, C. M., Experimental hydrogen isotope studies III: diffusion of hydrogen in hydrous

71.

72.

73.

74.

75.

76.

77.

78.

minerals, and stable isotoy‘e exchange in metamorphic rocks, Contrib. Mineral. Petrol., 76, 216-228, 1981. Graham, C. M., and Elphick, S. C., Some experimental con- straints on the role of hydrogen in oxygen and hydrogen diffusion and AI-Si interdiffu- sion in silicates, in Diffusion, Atomic Ordering, and Muss Transport: Selected Problems in Geochemistry, edited by Ganguly, J., pp. 248-285, Springer-Verlag, New York, 1991. Graham, C. M., Viglino, J. A., and Harmon, R. S., Experimental study of hydrogen- isotope exchange between aluminous chlorite and water and of hydrogen diffusion in chlorite, Am. Mineral., 72, 566-579. 1987. Griggs, D. T., Hydrolytic weakening of quartz and other silicates, Geophys. J., 14, 19- 31, 1967. Grove, T. L., Baker, M. B., and Kinzler, R. J., Coupled CaAI- NaSi diffusion in plagioclase feldspar: experiments and applications to cooling rate speedometry, Geochim. Cosmo- chim. Arto, 48, 2113-2121, 1984. Gupta, P. K., and Cooper, A. R., Jr., The [D] matrix for multicomponent diffusion, Physica, 54, 39-59, 197 I. Harrison, T. M., Diffusion of 40Ar in hornblende, Contrib. Mineral. Petrol., 78, 324-33 1, 1981. Harrison, T. M., Duncan, I., and McDougall, I., Diffusion of 40Ar in biotite: temperature, pressure and compositional effects, Geochim. Cosmochim. Acta, 4Y, 2461-2468, 1985. Harrison, T. M., and Watson,

79.

80.

81.

82.

83.

84.

85.

86.

87.

E. B., Kinetics of zircon dissolution and zirconium diffusion in granitic melts of variable water content, Contrib. Mineral. Petrol., 84, 66-72, 1983. Harrison, T. M., and Watson, E. B., The behavior of apatite during crustal anatexis: equilibrium and kinetic considerations, Geochim. Cosmochim. Acta, 48, 1467- 1477, 1984. Harrison, W. J., and Wood, B. J., An experimental investi- gation of the partitioning of REE between garnet and liquid with reference to the role of defect equilibria, Contrib. Mineral. Petrol., 72, 145- 155, 1980. Harrop, P. J., Self-diffusion in simple oxides (a bibliography), J. Materials Sci., 3, 206-222, 1968. Hart, S. R., Diffusion compen- sation in natural silicates, Geachim. Cosmochim. Acta, 45, 279-29 1, 198 1. Hartley, G. S., and Crank, J., Some fundamental definitions and concepts in diffusion processes. Trans. Faraduy Sot., 45, 801-818, 1949. Hayashi, T., and Muehlenbachs, K., Rapid oxygen diffusion in melilite and its relevance to meteorites, Geochim. Cosmo- chim. Actu, 50, 585-591, 1986. Hermeling, J., and Schmalzreid, H., Tracerdiffusion of the Fe- cations in Olivine (Fe,Mg,.,), SiO, (III), Phys. Chem. Miner., II, 161-166, 1984. Hofmann, A. W., Diffusion in natural silicate melts: a critical review, in Physics of Magmatic Processes, edited by Hargraves, R. B., pp. 385-417, Princeton Univ. Press, 1980. Hofmann, A. W., and Giletti,

88.

89.

90.

91.

92.

93.

94.

95.

B. J., Diffusion of geochron- ologically important nuclides in minerals under hydrothermal conditions, Eclogae Geol. Helv., 63, 141-150., 1970. Hofmann. A. W., Giletti, B. J., Hinthorne, J. R., Andersen, C. A., and Comaford, D., Ion microprobe analysis of a potassium self-diffusion experi- ment in biotite, Earth Planet. Sci. Lett., 24, 48-52, 1974. Hofmann, A. W., and Magaritz, M., Diffusion of Ca, Sr, Ba, and Co in a basalt melt, J. Geophys. RES., 82, 5432-5440, 1977.

Jambon, A., Zhang, Y., and

Houlier, B., Cheraghmakani, M., and Jaoul, O., Silicon

Stolper, E. M., Experimental

diffusion in San Carlos olivine, Phys. Earth Planet. Inter., 62,

dehydration of natural obsidian

329-340, 1990. Jambon, A., Tracer diffusion in

and estimation of DH20 at low

granitic melts, J. Geophys. Res., 87, 10797-10810, 1982.

water contents, Geochim.

Jambon, A., and Semet, M. P., Lithium diffusion in silicate glasses of albite, orthoclase, and

Cosmochim. Acta, 56, 293 l-

obsidian composition: an ion-

2935, 1992.

microprobe determination, Earth Planet. Sci. Lett., 37, 445-450,

Joesten, R., Evolution of

1978.

mineral assemblage zoning in diffusion metasomatism, Geo- chim. Cosmochim. Acta, 41, 649-670, 1977. Joesten, R., Grain-boundary diffusion kinetics in silicates and oxide minerals, in Diffusion, Atomic Ordering, and Mass Transport: Selected Problems in Geochemistry, edited bv Ganpulv. J.. DV. 345-

BRADY 287

395, Springer-Verlag, New York, 1991 a.

96. Joesten, R., Local equilibrium in metasomatic processes revisited. Diffusion-controlled growth of chert nodule reaction rims in dolomite, A m . Mineral., 76, 743-755, 199lb.

97. Joesten, R., and Fisher, G., Kinetics of diffusion-controlled mineral growth in the Christmas Mountains (Texas) contact aureole, Geol. Sot. Am. Bull., 100, 714-732, 1988.

98. Johnson, C. D., and Carlson, W. D., The origin of olivine- plagioclase coronas in meta- gabbros from the Adirondack Mountains, New York, J. Metamorph. Geol., 8, 697-7 17, 1990.

99. Jost, W., Diffusion in Solids, Liquids, and Gases, 558 pp., Academic Press, New York, 1960.

101. Kasper, R. B., Cation and oxygen diffusion in albite, Ph.D. thesis, Brown University, Providence, Rhode Island, 1975.

102. Kingery, W. D., Plausible con- cepts necessary and sufficient for interpretation of grain boundary phenomena: I, grain boundary characteristics, structure, and electrostatic potential, A m. Ceram. Sot. J., 57, l-8, 1974a.

103. Kingery, W. D., Plausible concepts necessary and sufficient for interpretation of grain

100. Jurewicz, A. J. G., and Watson,

boundary phenomena: II, solute

E. B., Cations in olivine, part

segregation, grain boundary

2: diffusion in olivine

diffusion,

xenocrysts, with applications to petrology and mineral physics,

and

Contrib. Mineral. Petrol., 99,

general

186-201, 1988.

discussion, Am. Ceram. Sot. J., 57, 74-83, 1974b.

288 DIFFUSION DATA

104. Kronenberg, A. K., Kirby, S. H., Aines, R. D., and Rossman, G. R., Solubility and diffu- sional uptake of hydrogen in quartz at high water pressures: implications for hydrolytic weakening, J. Geophys. Res., Y/. 1272:3-12734, 1986.

105. Lasaga, A. C., Multicomponent exchange and diffusion in silicates, Geochim. Cosmo- chim. Acta, 43, 455-469, 1979.

106. Lasaga, A. C., The atomistic basis of diffusion: defects in minerals, in Kinetics o,f Ceochemical Processes, edited by Lasaga, A. C., and Kirkpatrick, R. J., pp. 261-320, Mineralogical Society of America, Washington, 198 I.

107. Lasaga, A. C., Geospeedometry: an extension of geother- mometry, in Kinetics and Equilibrium in Mineral Reactions, edited by Saxena, S. K., pp. 81-114, Springer- Verlag, New York, 1983.

108. Lesher, C. E., and Walker, D., Thermal diffusion in petrology, in D$iision, Atomic Ordering, and Muss Transport: Selected Problems in Geochemistry, edited by Ganguly, J., pp. 396- 45 I, Springer-Verlag, New York, 1991.

109. Lin, T. H., and Yund, R. A., Potassium and sodium self- diffusion in alkali feldspar, Contrib. Mineral. Petrol., 34, 177-184, 1972.

I IO. Lindner, R., Studies on solid state reactions with radiotracers, J. Chem. Phys., 23, 4 10-4 I I, 1955.

1 I I. Lindner, R., Silikatbildung durch Reaktion im festen Zustand, Z. Physik. Chem., 6, 129-142, 1956.

112. LindstrGm, R., Chemical diffusion in alkali halides, J. Phys. Chem. Solids, 30, 401-

405, 1969. 113. Lindstriim, R., Chemical inter-

diffusion in binary ionic solid solutions and metal alloys with changes in volume, J. Phys. C. Solid St&e, 7, 3909-3929, 1974.

I 14. Liu, M., and Yund, R. A., NaSi-CaAI interdiffusion in plagioclase, Am. Mineral., 77, 275-283, 1992.

I IS. Loomis, T. P., Ganguly, J., and Elphick, S. C., Experimental determination of cation diffusiv- ities in aluminosilicate garnets II. Multicomponent simulation and tracer diffusion coefficients, Contrib. Mineral. Petrol., 90, 45-5 I, 1985.

I 16. Lowry, R. K., Henderson, P., and Nolan, J., Tracer diffusion of some alkali, alkaline-earth and transition element ions in a basaltic and an andesitic melt, and the implications concerning melt structure, Contrib. Mineral. Petrol., 80, 254-26 I, 1982.

1 17. Lowry, R. K., Reed, S. J. B., Nolan, J., Henderson, P., and Long, J. V. P., Lithium tracer- diffusion in an alkali-basalt melt -- an ion-microprobe determina- tion, Earth Planet. Sci. Lett., 53, 36-41, 1981.

I 18. Mackwell, S. J., and Kohlstedt, D. L., Diffusion of hydrogen in olivine: implications for water in the mantle, J. Geophys. Res., 9.5, 5079-5088, 1990.

I 19. Magaritz, M., and Hofmann, A. W., Diffusion of Sr, Ba, and Na in obsidian, Geochim. Cosmo- chim. Actu, 42, 595-605, 1978a.

120. Magaritz, M., and Hofmann, A. W., Diffusion of Eu and Gd in basalt and obsidian, Geochim. Cosmochim. Actu, 42, 847- 858, 1978b.

121. Manning, J. R., Diffusion

Kinetics for Atoms in Crystals, 257 pp., Van Nostrand, Princeton, 1968.

122. Matano, C., On the relation between the diffusion coeffi- cients and concentrations of solid metals, Jupan J. Phys., 8, 109-l 13. 1933.

123. McDougall, I., Harrison, T. M., Geochronology and Thermo- chronolog~~ by, the 40A r/‘9A r Method, 2 I2 pp., Oxford Univ. Press, New York, 1988.

124. Meibnei.. D. J., Cationic diffusion in olivine to 14OO’C and 3.5 kbar, in Geochemical Transport and Kinetics, edited by Hofmann, A. W., Giletti, B. J., Yoder, H. S., Jr., and Yund, R. A., pp. I 17-129, Carnegie Institution of Washington, Washington, 1974.

125. Morioka, M., Cation diffusion in olivine -- I. Cobalt and magnesium, Geochirn. Cosmo- chim. Acta, 44, 759-762, 1980.

126. Morioka, M., Cation diffusion in olivine -- II. Ni-Mg, Mn-Mg, Mg and Ca, Geochim. Cosmochim. Actu, 4S, I 5 73- 1580, 1981.

127. Morioka, M., and Nagasawa, H., Diffusion in single crystals of melilite: II. Cations, Geochim. Cosmochim. Acta, 55, 751-759, 1991a.

128. Morioka, M., and Nagasawa, H., Ionic diffusion in olivine, in Diffusion, Atomic Ordering, and Mass Transport: Selected Problems in Geochemistry, edited by Ganguly, J., pp. l76- 197, Springer-Verlag, New York, 1991b.

129. Morishita, Y., Giletti, B. J., and Farver, J. R., Strontium and oxygen self-diffusion in titanite (abstract), Eos Trans. AGU, 71, 652, 1990.

130. Muehlenbachs, K., and Connolly, C., Oxygen diffusion

BRADY 289

in leucite, in Stuhle Isotope Geochemistry: A Tribute to Samuel Epstein, edited by Taylor, H. P., Jr., O’Nejl, J. R., and Kaplan, 1. R., pp. 27- 34, The Geochemical Society, San Antonio, Texas, 1991.

13 I Nye, J. F., Physical Properties r)f Crystals, 322 pp., Clarendon Press, Oxford, 1957.

132. Onsager, L., Reciprocal rela- tions in irreversible processes, I, Physical Rev., 37, 405-426, 193la.

133. Onsager, L., Reciprocal rela- tions in irreversible processes II, Physical Rev., 38, 2265- 2279, 193 I b.

134. Onsager, L., Theories and problems of liquid diffusion, Ann. N. Y. Acad. Sci., 46, 24 I- 265, 1945.

135. Philibert, Jean, Atom Move- ments: Diffusion and Mass Transport in Solids, 577 pp., Editions de Physique, Les Ulis, France, 199 I.

136. Price, G. D., Diffusion in the titanomagnetite solid solution, Mineral. Mug., 44, 195-200, 1981.

137. Prigogine, I., Introduction to Thermodynamics of Irreversible Processes, 147 pp., Wiley Inter- science, New York, 1967.

138. Rapp, R. A., and Watson, E. B., Monazite solubility an dissolution kinetics: implica- tions for the thorium and light rare earth chemistry of felsic magmas., Contrib. Mineral. Petrol., 94, 304-3 16, 1986.

139. Reiss, Howard, Methods oj Thermodynamics, 217 pp., Blaisdell Publishing Company, Waltham, Massachusetts, 196.5.

140. Rovetta, M. R., Holloway, J. R., and Blacic, J. D., Solubility of hydroxl in natural quartz annealed in water at 900°C and 1.5 GPa, Geophys. Res. Lett.,

13, 145-148, 1986. 141. Ryerson, F. .I., Diffusion

measurements: experimental methods, Methods Exp. Phys., 24A, 89-130, 1987.

142. Ryerson, F. J., Durham, W. B., Cherniak, D. J., and Lanford, W. A., Oxygen diffusion in olivine: effect of oxygen fugacity and implications for creep, J. Geophps. Res., 94, 4105-4118, 1989.

143. Ryerson, F. J., and McKeegan, K. D., Determination of oxygen self-diffusion in Qkermanite, anorthite, diopside, and spinel: Implications for oxygen isotopic anomalies and the thermal histories of Ca-Al-rich inclusions, Geochim. Cosmo- chim. Acta, 58(17), in press, 1994.

144. Sautter, V., Jaoul, O., and Abel, F., Aluminum diffusion in diopside using the 27Al(p,y)28Si nuclear reaction: preliminary results, Earth Planet. Sci. Lett., 89, 109- 1 14, 1988.

145. Shaffer, E. W., Shi-Lan Sang, J., Cooper, A. R., and Heuer, A. H., Diffusion of tritiated water in b-quartz, in Ceochenz- ical Transport and Kinetics, edited by Hofmann, A. W., Giletti, B. J., Yoder, H. S., Jr., and Yund, R. A., pp. 131-138, Carnegie Institution of Wash., Washington, 1974.

146. Sharp, Z. D., Giletti, B. J., and Yoder, H. S., Jr., Oxygen diffu- sion rates in quartz exchanged with CO,, Earth Planet. Sci. Lett., 107, 339-348, 1991.

147. Shewmon, Paul G., Diffusion in Solids, 246 pp., Minerals, Metals & Materials Society, Warrendale, PA, 1989.

148. Smigelskas, A. D., and Kirkendall, E. O., Zinc diffu- sion in alpha brass, Am. Inst.

Min. Metal. Engrs. Trans., 171, 130-142, 1947.

149. Sneeringer, M., Hart, S. R., and Shimizu, N., Strontium and samarium diffusion in diopside, Geochim. Cosmochim. Acta, 48, 1589-1608, 1984.

150. Tiernan, R. J., Diffusion of thallium chloride into single crystals and bicrystals of potassium chloride, Ph.D. thesis, MIT, Cambridge, Massachusetts, 1969.

151. Trull, T. W., and Kurz, M. D., Experimental measurement of 3He and 4He mobility in olivine and clinopyroxene at magmatic temperatures, Geochim. Cosmo- chim. Acta, 57, 1313-1324, 1993.

152. Watson, E. B., Calcium diffusion in a simple silicate melt to 30 kbar, Geochim. Cosmochim. Acta, 43, 3 I 3- 322, 1979a.

IS?. Watson, E. B., Diffusion of cesium ions in H,O-saturated granitic melt, Science, 205, 1259-1260, 1979b.

154. Watson, E. B., Diffusion in fluid-bearing and slightly-melted rocks: experimental and numer- ical approaches illustrated by iron transport in dunite, Contrib. Mineral. Petrol., 107, 417-434, 1991a.

155. Watson, E. B., Diffusion of dissolved CO, and Cl in hydrous silicic to intermediate magmas, Geochim. Cosmochim. Acta, 55, 1897-1902, 1991b.

156. Watson, E. B., and Baker, D. R., Chemical diffusion in magmas: an overview of experimental results and geochemical applications, in Physical Chemistry of Magmas, Advances in Physical Geo- chemistry, Y, edited by Perchuk, L., and Kushiro, I., pp. 120- 15 1, Springer-Vet-lag, New

290 DIFFUSION DATA

York, I YY I. 157. Watson, E. B., Harrison, T. M.,

and Ryerson, F. J., Diffusion of Sm, Ar, and Pb in tluorapatite, Geochim. Cosmochim. Acta, 4Y, 1813-1823, 1985.

1.58. Watson, E. B., and Lupulescu, A., Aqueous tluid connectivity and chemical transport in clino- pyroxene-rich rocks, Earth Planet. Sci. Lett., 117, 27Y- 294, 1993.

159. Watson, E. B., Sneeringer, M. A., and Ross, A., Diffusion of dissolved carbonate in magmas: experimental results and appli- cations, Earth Planet. Sci. Lett., 61, 346-358, 1982.

160. Wendlandt, R. F., Oxygen diffusion in basalt and andesite melts: experimental results and

discussion of chemical versus tracer diffusion, Contrib. Mineral. Petrol., 108, 463-47 I, 1991.

I6 I. Wendt, R. P., The estimation of diffusion coefficients for ternary systems of strong and weak electrolytes, J. Phys. Chem., 69, 1227-1237, 1965.

162. Winchell, P., The compensation law for diffusion in silicates, High Temp. Sci., I, 200-2 15, 1969.

163. Yund, R. A., Diffusion in feldspars, in Feldspar Mineral- ogy, edited by Ribbe, P. H., pp. 203-222, Mineralogical Society of America, Washington, 1983.

164. Yund, R. A., Quigley, J., and Tullis, J., The effect of dislocations on bulk diffusion in

feldspars during metamorphism, J. Metamorph. Geol., 7, 337- 341, 1989.

165. Yurimoto, H., Morioka, M., and Nagasawa, H., Diffusion in single crystals of melilite: I. oxygen, Geochim. Cosmochim. Actu, 53, 2387-2394, 1989.

166. Zhang, Y., Stolper, E. M., and Wasserburg, G. J., Diffusion in rhyolite glasses, Geochim. Cosmochim. Acta, 55, 44 I - 456, 199la.

167. Zhang, Y., Stolper, E. M., and Wasserburg, G. J., Diffusion of a multi-species component and its role in oxygen and water transport in silicates, Earth flunet. Sci. Lett., 103, 228- 240. 199lb.