Author index volume 40

Tracer diffusion coefficients for Li in glasses of albite, orthoclase, and obsidian composition have been determined by a method involving deposition of a thin source on polished glass wafers, anneal under controlled temperature conditions (300–900°C), and ion-microprobe determination of the concentration profile. All results conform to an Arrhenius-type relationship,D = D₀exp(?Q/RT), whereQ is 23, 17, and 22 kcal mol^?1;D₀ is 0.2, 0.003, and 0.03 cm²s^?1 for albite and orthoclase glasses, and obsidian respectively. Lithium is thus a fast diffusing ion and behaves similarly to sodium in the same glasses. A mechanism involving jumps of the diffusing ions through oxygen hexagonal rings is suggested by consideration of ionic radii ratio of alkali (H, Li, Na, K, Rb, and Cs) ions to the oxygen anions.     

Compositional gradients surrounding spherulites in obsidian and their relationship to spherulite growth and lava cooling

Spherical masses of crystal fibers (spherulites) crystalize from rhyolitic melt/glass mainly in response to significant undercooling while lava cools. Spherulite growth should induce compositional gradients in the surrounding glass from expulsion of incompatible constituents and diffusion of those constituents away from the spherulite. Finite-difference numerical modeling of one-dimensional diffusion, in which diffusivities are allowed to vary with temperature, is used to investigate how compositional gradients reflect spherulite growth and lava cooling. Overall, three forms of gradients are identified. Elements that diffuse quickly are expelled from the spherulite but then migrate away too quickly to become enriched at the boundary of the spherulite. Elements that diffuse slowly are trapped within the growing spherulite. Between those endmembers are elements that are not trapped, yet diffuse slow enough that they become enriched at the contact. Their slow diffusion away then elevates their concentrations in the surrounding glass. How enriched those elements are at the spherulite-matrix interface and how far their enrichments extend outwards into the glass reflect how spherulites grow and thermal conditions during growth. Concentrations of H₂O, Rb, F, Li, Cl, Na, K, Sr, Cs, Ba, and Be were measured in and around spherulites in obsidian from a 4.7?±?1?km³ rhyolite lava dome erupted from Tequila volcano, Mexico. Measurable concentration gradients are found for H₂O, Rb, and F. Attributes of those gradients and the behaviors of the other elements are in accord with their experimentally constrained diffusivities. Spherulites appear to have grown following radial, rather than volumetric, growth. The observed gradients (and lack of others) are more consistent with growth mainly below the glass transition, which would necessitate the dome cooling at ca. 10^?5 to 10^?7?°C?s^?1. Such slow cooling is consistent with the relatively large volume of the dome.     

Ion microprobe studies of water in silicate melts. Concentration-dependent water diffusion in obsidian

The ion microprobe at Johnson Space Center has been calibrated for in situ water determinations on a 10-μm scale over the range 0.2 wt.% H₂O to 1.8, 6.8, and 3.7 wt.%, for basaltic, albitic, and rhyolitic glasses, respectively. The basalt glass calibration curve differs substantially from those of albite and rhyolite glasses, indicating a need to carefully match composition and/or melt structure between H₂O standards and unknowns.A value for the diffusivity of water as a function of concentration and time has been calculated from water diffusion profiles measured in rhyolite glasses prepared at 850°C and 700 barsP_t(H₂O) [1]. Transient diffusion into a semi-infinite medium is described by the equation:?(φ/2)?¸/?φ=?(D_w?¸/?φ)/?φ #x003B8;=1, φ=0, θ→ 0, θ→∞, wherex =distance from the cylinder edge,t =time,C₀ =initial concentration,C_s =concentration at the edge,C =concentration at x,θ = C ? C₀/C_s ? C₀,φ = x/t^1/2, andD_w =diffusivity of water. An iterative technique has been used to calculate solutions to the diffusion equation as a function ofD_w [2]. Comparison of these solutions with the ion probe data indicate that, for0.2wt.% ≤ C ≤ 3.7wt.%H₂O,D_w can be described by an exponential function of θ, of the formD_w = D₀exp(bθ), withD₀ (i.e.,D_w at 0.2%) = (0.8?2.2) × 10^?8 cm²/s and2 ≤ b ≤ 4.     

Ion microprobe studies of water in silicate melts: Temperature-dependent water diffusion in obsidian

The temperature dependence of water diffusivity in rhyolite melts over the range 650–950°C and [P_T(H₂O] = 700 bars is evaluated from water concentration-distance profiles measured in glass with an ion microprobe. Diffusivities are exponentially dependent on concentration over this temperature range and vary from about 10^?8 cm²/s at 650°C to about 10^?7 cm²/s at 950°C at 2 wt.% water. Water solubility also varies with temperature at a rate of ?0.14 wt. per 100°C increase. The avtivation energy (E_a) appears to be constant at 19 ± 1kal/mole for 1, 2,and 3 wt.% H₂O. Comparison of these data with results for cation diffusion indicates that this value is a minimum E_a for diffusion of any species in a rhyolite melt.Compensation plots of log₁₀D₀ (the frequency factor) versus E_a indicate that hydrous rhyolite melts follow the same trend as anhydrous basalts. D₀ increases for H₂O and Ca²⁺ [1] as E_a decreases. This suggests that these molecules may diffuse by different mechanisms than do monovalent cations, and that hydration of the melt affects diffusion of Ca²⁺ and H₂O differently than it does monovalent cation diffusion. The results imply that dramatic increases in cation diffusivities by hydration [1] may occur with additions of less than 1 wt.% H₂O.     

Color-change processes of a plinian pumice and experimental constraints of color-change kinetics in air of an obsidian

Colors of plinian pumices were measured by spectrocolorimetry, and their quantitative color parameters in the L*a*b* color space were determined. A series of heating experiments of obsidian was conducted to simulate the color-change processes of rhyolitic glasses. In these experiments, following three stages of color-change processes were observed. Stage I showed a rapid b* (yellowishness) increase associated with fast dehydration controlled by water diffusivity (D _water). In stage II, a* (reddishness) increase was accompanied by Fe²⁺ decrease. Both a* increase and Fe²⁺ decrease can be simulated by a diffusion model. Obtained diffusivity D _oxidation were about two orders of magnitude smaller than D _water . The a*-value increase after the oxidation in stage III appeared to be quasi-linear with time, indicating the zeroth order reaction corresponding to the formation of hematite-like structures in rhyolitic glasses. The diffusion-limited a* increase model in stage II was applied to a natural plinian pumice fall unit to evaluate time periods of color-change processes through oxidation by air of fragmented rhyolitic materials.     

Pan evaporation modelling and changing attribution analysis on the Tibetan Plateau (1970–2012)

Pan evaporation (E_p) is an important indicator of water and energy and the decline of E_p has been reported in many regions over the last decades. The climate and E_p are dependent on each other. In this study, the temporal trends of E_p and main E_p drivers, namely mean air temperature (T_a), wind speed (u), global solar radiation (R_s), net long‐wave radiation(R_nl) and vapour pressure deficit (D) from 1970 to 2012, were calculated on the basis of 26 meteorological stations on the Tibetan Plateau. The arithmetic average of E_p from 26 stations decreased with the rate of ?11.91 mm a^?2; the trends of R_s, R_nl, T_a, u and D were ?1.434 w m^?2 decade^?1, 0.2511 w m^?2 decade^?1, 0.3590°C decade^?1, ?0.2376 m s^?1 decade^?1 and 9.523 Pa decade^?1, respectively. The diffuse irradiance is an essential parameter to model E_p and quantify the contribution of climatic factors to changing E_p. 60 724 observations of R_s and diffuse solar irradiance (R_d) from seven of the 26 stations were used to develop the correlation between the diffuse fraction (R_d/R_s), and the clearness index (R_s/R_o). On the basis of the estimation of the diffuse component of R_s and climatic data, we modified the PenPan model to estimate Chinese micro‐pan evaporation (E_p) and assess the attribution of E_p dynamics using partial derivatives. The results showed that there was a good agreement between the observed and calculated daily E_p values. The observed decrease in E_p was mostly due to declining wind speed (?13.7 mm a^?2) with some contributions from decreasing solar irradiance (?3.1 mm a^?2); and the increase of temperature had a large positive effect (4.55 mm a^?2) in total whilst the increase of R_nl had insignificant effect (0.35 mm a^?2) on E_p rates. The change of E_p is the net result of all the climatic variables. Copyright © 2014 John Wiley & Sons, Ltd.     

Absolute18O content of standard mean ocean water

The absolute values of the¹⁸O/¹⁶O ratio (R_s) and the relative¹⁸O content (X_s) in SMOW have been determined by comparing SMOW mass spectrometrically with well-defined synthetic mixtures of pure D₂¹⁸O and H₂¹⁶O. The results are:R_S = (2005.20 ± 0.45) × 10^?6, X_S = (2000.45 ± 0.45) × 10^?6     

Measurement of the adiabats of quartz,forsterite, and magnesium oxide at high pressures and high temperatures and adiabatic gradient in the mantle

The geophysically-important adiabat (?T/?P)_s has been measured at pressures up to 50 kbar and temperatures up to 1000°C. A simple power law describes the relationship between (?T/?P)_s and the compression of the material. The power is independent of the material and of the temperature within the uncertainty. This consistency in the power allows the extrapolation of the adiabat to pressure and temperature conditions of the mantle of the earth. The adiabatic gradient is shown to be significantly smaller than the melting gradient.     

Thermal expansion of a few Indian granitic rocks

The study of thermal expansion by a dilatometer technique on a few granitic rocks from the Peninsular shield and Himalayan regions of India confirms that the linear coefficient of thermal expansion (α) is a function of heating rate, crack porosity, thermal cycling, mineral composition and grain orientation. Permanent set in the samples occurs at the limiting temperature (T_p) and restricts the validity of the apparent thermal-expansion coefficient with rise in temperature. Values of α are determined for a heating rate of ?2°C min^?1 in order to calculate the volume coefficient of expansion (α_v) and the temperature dependence of density (ρ_T).     

Degassing of Cl, F, Li, and Be during extrusion and crystallization of the rhyolite dome at Volcán Chaitén, Chile during 2008 and 2009

We investigated the distribution of Cl, F, Li, and Be in pumices, obsidians, and crystallized dome rocks at Chaitén volcano in 2008?C2009 in order to explore the behavior of these elements during explosive and effusive volcanic activity. Electron and ion microprobe analyses of matrix and inclusion glasses from pumice, obsidian, and microlite-rich dome rock indicate that Cl and other elements were lost primarily during crystallization of the rhyolitic dome after it had approached the surface. Glass in pumice and microlite-free obsidian has 888?±?121?ppm Cl, whereas residual glass in evolved microlite-rich dome rock generally retains less Cl (as low as <100?ppm). Estimated Cl losses were likely >0.7?Mt Cl, with a potential maximum of 1.8?Mt for the entire 0.8-km³ dome. Elemental variations reflect an integrated bulk distribution ratio for Cl?>?1.7 (1.7 times more Cl was degassed or incorporated into crystals than remained in the melt). Because Cl is lost dominantly as the very last H₂O is degassed, and Cl is minimally (if at all) partitioned into microlites, the integrated vapor/melt distribution ratio for Cl exceeds 200 (200 times more Cl in the evolved vapor than in the melt). Cl is likely lost as HCl, which is readily partitioned into magmatic vapor at low pressure. Cl loss is accelerated by the change in the composition of the residual melt due to microlite growth. Cl loss also may be affected by open-system gas fluxing. Integrated vapor-melt distribution ratios for Li, F, and Be all exceed 1,000. On degassing, an unknown fraction of these volatiles could be immediately dissolved in rainwater.     

Quantification of the Mass Flux of H2O Gas (Steam) at Three Active Volcanoes Using Thermal Infrared Imagery

We apply a measurement technique that utilizes thermal video of vapor-dominated volcanic plumes to estimate the H₂O gas flux at three degassing volcanoes. Results are compared with H₂O flux measurements obtained using other methods to verify the thermal camera-derived values. Our estimation of the H₂O emission rate is based on the mass and energy conservation equations. H₂O flux is quantified by extracting the temperature and width of the gas plume from the thermal images, calculating the transit velocity of the gas plume from the thermal video, and combining these results with atmospheric parameters measured on-site. These data are then input into the equations for conservation of mass and energy. Selected volcanoes for this study were Villarrica in Chile, Stromboli in Italy, and Santa Ana in El Salvador. H₂O fluxes estimated from the thermal imagery were 38–250?kg?s^?1 at Villarrica, 4.5–14?kg?s^?1 for Stromboli’s Central Crater, and 168–219?kg?s^?1 at Santa Ana. These compare with H₂O flux values estimated by other methods of 73–220, 3–70 and 266?kg?s^?1, at the three volcanoes, respectively. The good agreement between thermal image-derived results and those estimated by other methods seem to validate this method.     

Influence of advection on the characteristics of turbulence over uneven surface in the oasis and the Gobi Desert

Utilizing experimental data of the atmospheric surface layer in the Gobi Oasis of Jinta in a comparative study, we demonstrate that under the condition of unstable stratification, the normalization variances of temperature in the oasis and Gobi Desert meet 1/3(z/Λ)(z/Λ)while normalization variances of both humidity and CO2 in the oasis meet(z)sz\Λ-1/3;s s z Λ the normalization variance of temperature in the oasis is large due to disturbance by advection, whereas variance of CO2 in the Gobi Desert has certain degree of deviation relative to Monin-Obukhov(M-O) scaling, and humidity variance completely deviates from variance M-O scaling. The above result indicates that under the condition of advection, humidity variance meets the relationsm sA sB D and it is determined by relative magnitude of scalar variance of advection transport. Our study reveals that, if the scalar variance of humidity or CO2 transported by advection is much larger than local scalar variance, observation value of scalar variance will deviate from M-O scaling; when scalar variance of advection transport is close to or less than local scalar variance, the observation value of scalar variance approximately meets M-O scaling.     

The238U spontaneous fission decay constant re-determined by fission tracks

The decay constant(λ_f²³⁸) for the spontaneous fission of²³⁸U was re-determined by means of a man-made uranium glass of known age (126 yr). The spontaneous U fission tracks that had accumulated since the date of manufacture were counted on internal faces of the glass with an error of less than 1.7%. No thermal annealing of the spontaneous tracks was observed. The U content was determined by induced fission tracks. The value obtained forλ_f²³⁸ is(8.57 ± 0.42) × 10^?17yr^?1. Main sources of error are the date of glass melting and the determination of the thermal neutron dose.     

Erodibility of low‐compaction steep‐sloped reclaimed surface mine lands in the southern Appalachian region,USA

The use of loose spoils on steep slopes for surface coal mining reclamation sites has been promoted by the US Department of Interior, Office of Surface Mining for the establishment of native forest, as prescribed by the Forest Reclamation Approach (FRA). Although low‐compaction spoils improve tree survival and growth, erodibility on steep slopes was suspected to increase. This study quantified a combined K_C factor (combining the effects of the soil erodibility K factor and cover management C) for low compaction, steep‐sloped (>20°) reclaimed mine lands in the Appalachian region, USA. The combined K_C factor was used because standard Unit Plot conditions required to separate these factors, per Revised Universal Soil Loss Equation (RUSLE) experimental protocols, were not followed explicitly. Three active coal mining sites in the Appalachian region of East Tennessee, each containing four replicate field plots, were monitored for rainfall and sediment yields during a 14‐month period beginning June 2009. Average cumulative erosivity for the study sites during the monitoring period was measured as 5248.9 MJ·mm·ha^?1·h^?1. The K_C ranged between 0.001 and 0.05 t·ha·h·ha^?1·MJ^?1·mm^?1, with the highest values occurring immediately following reclamation site construction as rills developed (June – August 2009). The K_C for two study sites with about an 18–20 mm spoil D₈₄ were above 0.01 t·ha·h·ha^?1·MJ^?1·mm^?1 during rill development, and below 0.003 t·ha·h·ha^?1·MJ^?1·mm^?1 after August 2009 for the post‐rill development period. The K_C values for one site with a 40 mm spoil D₈₄ were never above 0.008 t·ha·h·ha^?1·MJ^?1·mm^?1 and also on average were lower, being more similar to the other two sites after the rill development period. Based on an initial K_C factor (K_e) measured during the first few storm events, the average C factor (C_e) was estimated as 0.58 for the rill development period and 0.13 for the post‐rill development period. It appears that larger size fractions of spoils influence K_C and C_e factors on low‐compaction steep slopes reclamation sites. Copyright © 2013 John Wiley & Sons, Ltd.     

Apatite/liquid partition coefficients for the rare earth elements and strontium

Sixteen sets of apatite/liquid partition coefficients (D^ap/liq) for the rare earth elements (REE; La, Sm, Dy, Lu) and six values for Sr were experimentally determined in natural systems ranging from basanite to granite. The apatite + melt (glass) assemblages were obtained from starting glasses artificially enriched in REE, Sr and fluorapatite components; these were run under dry and hydrous conditions of 7.5–20 kbar and 950–1120°C in a solid-media, piston-cylinder apparatus. An SEM-equipped electron microprobe was used for subsequent measurement of REE and Sr concentrations in coexisting apatites and quenched glasses. The resulting partition coefficient patterns resemble previously determined apatite phenocryst/groundmass concentration ratios in the following respects: (1) the rare earth patterns are uniformly concave downward (i.e., the middle REE are more compatible in apatite than the light and heavy REE); (2) D_REE^ap/liq is much higher for silicic melts than for basic ones; and (3) strontium (and therefore Eu²⁺) is less concentrated by apatite than are the trivalent REE. The effects of both temperature and melt composition on D_REE^ap/liq are systematic and pronounced. At 950°C, for example, a change in melt SiO₂ content from 50 to 68 wt.% causes the average REE partition coefficient to increase from ～7 to ～30. A 130°C increase in temperature, on the other hand, results in a two-fold decrease in D_REE^ap/liq. Partitioning of Sr is insenstitive to changes in melt composition and temperature, and neither the Sr nor the REE partition coefficients appear to be affected by variations in pressure or H₂O content of the melt.The experimentally determined partition coefficients can be used not only in trace element modelling, but also to distinguish apatite phenocrysts from xenocrysts in rocks. Reported apatite megacryst/host basalt REE concentration ratios [12], for example, are considerably higher than the equilibrium partition coefficients, which suggest that in this particular case the apatite is actually xenocrystic.A reversal experiment incorporated in our study yielded diffusion profiles of REE in apatite, from which we extracted a REEαCa interdiffusion coefficient of 2–4×10^?14 cm²/s at 1120°C. Extrapolated downward to crustal temperatures, this low value suggests that complete REE equilibrium between felsic partial melts and residual apatite is rarely established.     

The thermal boundary-layer interpretation of D″ and its role as a plume source

The “anomalous” layer in the lowermost mantle, identified as D″ in the notation of K.E. Bullen, appears in the PREM Earth model as a 150 km-thick zone in which the gradient of incompressibility with pressure, $\text{d}\text{K}\text{d}\text{P}$, is almost 1.6, instead of 3.2 as in the overlying mantle. Since PREM shows no accompanying change in density or density gradient, we identify D″ as a thermal boundary layer and not as a chemically distinct zone. The anomaly in $\text{d}\text{K}\text{d}\text{P}$ is related to the temperature gradient by the temperature dependence of K_s, for which we present a thermodynamic identity in terms of accessible quantities. This gives the numerical result (?K_s/?T)_P=?1.6×10⁷ Pa K^?1 for D″ material. The corresponding temperature increment over the D″ range is 840 K. Such a layer cannot be a static feature, but must be maintained by a downward motion of the lower mantle toward the core-mantle boundary with a strong horizontal flow near the base of D″. Assuming a core heat flux of 1.6 × 10¹² W, the downward speed is 0.07 mm y^?1 and the temperature profile in D″, scaled to match PREM data, is approximately exponential with a scale height of 73 km. The inferred thermal conductivity is 1.2 W m^?1 K^?1. Using these values we develop a new analytical model of D″ which is dynamically and thermally consistent. In this model, the lower-mantle material is heated and softened as it moves down into D″ where the strong temperature dependence of viscosity concentrates the horizontal flow in a layer ～ 12 km thick and similarly ensures its removal via narrow plumes.     

Magnetic properties of synthetic titanomagnetites

Magnetic properties and crystal structure parameters of synthetic solid solutions Fe₃O₄Fe₃TiO₄, Fe₂O₄MgFe₂O₄ and Fe₃O₄Mg₂TiO₄ have been studied. Basic regularities in the behaviour of saturation magnetisation (I_s), Curie temperature (T_C) and cubic lattice parameter a during the substitution of Ti and Mg ions for Fe ions have been found. As the concentration of Ti ions increases, I_s reduces from 70 Gs·cm³ g^?1 to 0, T_C changes from 580 to 130°C and a from 8.391 to 8.520 Å. Growth of the Mg concentration leads to changes in I_s to 19.8 Gs·cm³, g^?1, T_C, to 440°C and a, to 8.360 Å. The full Fe ions substitution gives $\text{“a”=8.440}\text{A}\text{?}$.Chemical compositions of the samples, in which the valency variation of Fe ions at oxidation leads to an increase in susceptibility and T_C, have been determined.     

CO<Subscript>2</Subscript> flux evaluation over the evergreen coniferous and broad-leaved mixed forest in Dinghushan,China

The Dinghushan flux observation site, as one of the four forest sites of ChinaFLUX, aims to acquire long-term measurements of CO₂ flux over a typical southern subtropical evergreen coniferous and broad-leaved mixed forest ecosystem using the open path eddy covariance method. Based on two years of data from 2003 to 2004, the characteristics of temporal variation in CO₂ flux and its response to environmental factors in the forest ecosystem are analyzed. Provided two-dimensional coordinate rotation, WPL correction and quality control, poor energy-balance and underestimation of ecosystem respiration during nighttime implied that there could be a CO₂ leak during the nighttime at the site. Using daytime (PAR > 1.0 μmol^?1·m^?2·s^?1) flux data during windy conditions (u* > 0.2 m·s^?1), monthly ecosystem respiration (Reco) was derived through the Michaelis-Menten equation modeling the relationship between net ecosystem C0₂ exchange (NEE) and photosynthetically active radiation (PAR). Exponential function was employed to describe the relationship between Reco and soil temperature at 5 cm depth (Ts05), then Reco of both daytime and nighttime was calculated respectively by the function. The major results are: (i) Derived from the Michaelis-Menten equation, the apparent quantum yield (α) was 0.0027±0.0011 mgCO₂·μmol^?1 photons, and the maximum photosynthetic assimilation rate (Amax) was 1.102±0.288 mgCO₂·m^?2·s^?1. Indistinctive seasonal variation of α or Amax was consistent with weak seasonal dynamics of leaf area index (LAf) in such a lower subtropical evergreen mixed forest, (ii) Monthly accumulated Reco was estimated as 95.3±21.1 gC·m^?2mon^?1, accounting for about 68% of the gross primary product (GPP). Monthly accumulated WEE was estimated as ?43.2±29.6 gC·m^?2·mon^?1. The forest ecosystem acted as carbon sink all year round without any seasonal carbon efflux period. Annual NEE of 2003 and 2004 was estimated as ?563.0 and ?441.2 gC·m^?2·a^?1 respectively, accounting for about 32% of GPP.     

The relationship between vertical eddy diffusion and buoyancy gradient in the deep sea

Vertical eddy diffusivities (K_v's) have been estimated at fourteen widely separated locations from fourteen²²²Rn profiles and two²²⁸Ra profiles measured near the ocean floor as part of the Atlantic and Pacific GEOSECS programs. They show an inverse proportionality to the local buoyancy gradient [(g/?)(??_pot/?z)] calculated from hydrographic measurements. The negative of the constant of proportionality is the buoyancy flux [?K_v(g/?)(??_pot/?z)] which has a mean of ?4 × 10^?6 cm²/sec³. Our results suggest that the buoyancy flux varies very little near the ocean floor. K_v's for the interior of the deep Pacific calculated from the relationship K_v = (4 × 10^?6cm²/sec³)/[(g/?)(??_pot/?z)] agree well with published estimates. K_v's calculated for the pycnocline are one to two orders of magnitude smaller than upper limits estimated from tritium and⁷Be distributions.Heat fluxes calculated with the model K_v's obtained from the²²²Rn profiles average 31 μcal cm^?2 sec^?1 in the Atlantic Ocean and 8 μcal cm^?2 sec^?1 in the Pacific Ocean.     

Rare earth element systematics of hydrous liquids from partial melting of basaltic eclogite a re-evaluation

The origin of orogenic andesitic magmas is tested by calculations of REE fractionation in hydrous melts derived from partial melting of subducted ocean basalt in eclogite facies. New data on the subsolidus phase proportions of basaltic eclogite, the enrichment of LREE in altered ocean basalts, and experimentally determined REE partition coefficients (K_D's) between garnet and melt have been included in trace element fractionation equations. Non-modal melting of phases combined with variation inK_D's during melting is a unique feature of these calculations.Variation ofK_D, melting proportions, initial proportion of subsolidus phases, degree of melting, and initial REE concentrations yield a wide range of input parameters that produce REE profiles in partial melts of basaltic eclogite matching REE profiles of some orogenic andesites. The positive correlation of REE concentration with silica content for many andesitic suites can be accounted for by non-modal melting if quartz (or a similar phase with low REEK_D values) melts at a high melting proportion and garnet melts at a low melting proportion during the first stages of fusion. However, no mineralogic fractionation scheme can account for REE/silica systematics if REEK_D values are linearly decreasing with increasing melting. Earlier workers who have used similar calculations to discredit the eclogite fractionation model have set overly strict, and sometimes incorrect, constraints concerning the range in REEK_D values for garnet, the subsolidus proportions of phases in basaltic eclogite, and the relative concentrations of REE in subducted ocean crust undergoing partial melting.