Empirical models relating viscosity and tracer diffusion in magmatic silicate melts

The Adam-Gibbs equations describing relaxation in silicate melts are applied to diffusion of trace components of multicomponent liquids. The Adam-Gibbs theory is used as a starting point to derive an explicit relation between viscosity and diffusion including non-Arrhenian temperature dependence. The general form of the equation is D_iη = A_iexp{Δ(s_cE_i)/TS_c}, where D is diffusivity, η is melt viscosity, T is absolute temperature, Δ(s_cE_i) is the difference between the products of activation energies and local configurational entropies for viscous and diffusive relaxation, A_i is a constant that depends on the characteristics of the diffusing solute particles, and S_c is configurational entropy of the melt. The general equation will be impractical for most predictive purposes due to the paucity of configurational entropy data for silicate melts. Under most magmatic conditions the proposed non-Arrhenian behaviour can be neglected, allowing the general equation to be simplified to a generalized form of the Eyring equation to describe diffusion of solutes that interact weakly with the melt structure: D_iη/T = Q_iexp{ΔE_i/RT}, where Q_i and ΔE_i depend on the characteristics of the solute and the melt structure. If the diffusing solute interacts strongly with the melt structure or is a network-forming cation itself, then ΔE_i = 0, and the relation between viscosity and diffusion has the functional form of the classic Eyring and Stokes-Einstein equations; D_iη/T = Q_i. If the diffusing solute can make diffusive jumps without requiring cooperative rearrangement of the melt structure, the diffusivity is entirely decoupled from melt viscosity and should be Arrhenian, i.e., D_i = Q_iexp{B_i/T}. A dataset of 594 published diffusivities in melts ranging from the system CAS through diopside, basalt, andesite, anhydrous rhyolite, hydrous rhyolite, and peralkaline rhyolite to albite, orthoclase, and jadeite is compared with the model equations. Alkali diffusion is completely decoupled from melt viscosity but is related to melt structure. Network-modifying cations with field strength Z_i²/r between 1 and 10 interact weakly with the melt network and can be modelled with the extended form of the Eyring equation. Diffusivities of cations with high field strength have activation energies essentially equal to that of viscous flow and can be modelled with a simple reciprocal Eyring-type dependence on viscosity. The values of Q_i, ΔE_i and B_i for each cation are different and can be related to the cation charge and radius as well as the composition of the melt through the parameters Z_i²/r, M/O, and Al/(Na + K + H). I present empirical fit parameters to the model equations that permit prediction of cation diffusivities given only charge and radius of the cation and temperature, composition and viscosity of the melt, for the entire range of temperatures accessible to magmas near to or above their liquidus, for magmas ranging in composition from basalt through andesite to hydrous or anhydrous rhyolite. Pressure effects are implicitly accounted for by corrections to melt viscosity. Ninety percent of diffusivities predicted by the models are within 0.6 log units of the measured values.     

Implications of Nb/U, Th/U and Sm/Nd in plume magmas for the relationship between continental and oceanic crust formation and the development of the depleted mantle

The Nb/U and Th/U of the primitive mantle are 34 and 4.04 respectively, which compare with 9.7 and 3.96 for the continental crust. Extraction of continental crust from the mantle therefore has a profound influence on its Nb/U but little influence on its Th/U. Conversely, extraction of midocean ridge-type basalts lowers the Th/U of the mantle residue but has little influence on its Nb/U. As a consequence, variations in Th/U and Nb/U with Sm/Nd can be used to evaluate the relative importance of continental and basaltic crust extraction in the formation of the depleted (Sm/Nd enriched) mantle reservoir.This study evaluates Nb/U, Th/U, and Sm/Nd variations in suites of komatiites, picrites, and their associated basalts, of various ages, to determine whether basalt and/or continental crust have been extracted from their source region. Emphasis is placed on komatiites and picrites because they formed at high degrees of partial melting and are expected to have Nb/U, Th/U, and Sm/Nd that are essentially the same as the mantle that melted to produce them. The results show that all of the studied suites, with the exception of the Barberton, have had both continental crust and basaltic crust extracted from their mantle source region. The high Sm/Nd of the Gorgona and Munro komatiites require the elevated ratios seen in these suites to be due primarily to extraction of basaltic crust from their source regions, whereas basaltic and continental crust extraction are of subequal importance in the source regions of the Yilgarn and Belingwe komatiites. The Sm/Nd of modern midocean ridge basalts lies above the crustal extraction curve on a plot of Sm/Nd against Nb/U, which requires the upper mantle to have had both basaltic and continental crust extracted from it.It is suggested that the extraction of the basaltic reservoir from the mantle occurs at midocean ridges and that the basaltic crust, together with its complementary depleted mantle residue, is subducted to the core-mantle boundary. When the two components reach thermal equilibrium with their surroundings, the lighter depleted component separates from the denser basaltic component. Both are eventually returned to the upper mantle, but the lighter depleted component has a shorter residence time in the lower mantle than the denser basaltic component. If the difference in the recycling times for the basaltic and depleted components is ∼1.0 to 1.5 Ga, a basaltic reservoir is created in the lower mantle, equivalent to the amount of basalt that is subducted in 1.0 to 1.5 Ga, and that reservoir is isolated from the upper mantle. It is this reservoir that is responsible for the Sm/Nd ratio of the upper mantle lying above the trend predicted by extraction of continental crust on the plot of Sm/Nd against Nb/U.     

Geochemical evidence for diversity of dust sources in the southwestern United States

Several potential dust sources, including generic sources of sparsely vegetated alluvium, playa deposits, and anthropogenic emissions, as well as the area around Owens Lake, California, affect the composition of modern dust in the southwestern United States. A comparison of geochemical analyses of modern and old (a few thousand years) dust with samples of potential local sources suggests that dusts reflect four primary sources: (1) alluvial sediments (represented by Hf, K, Rb, Zr, and rare-earth elements, (2) playas, most of which produce calcareous dust (Sr, associated with Ca), (3) the area of Owens (dry) Lake, a human-induced playa (As, Ba, Li, Pb, Sb, and Sr), and (4) anthropogenic and/or volcanic emissions (As, Cr, Ni, and Sb). A comparison of dust and source samples with previous analyses shows that Owens (dry) Lake and mining wastes from the adjacent Cerro Gordo mining district are the primary sources of As, Ba, Li, and Pb in dusts from Owens Valley. Decreases in dust contents of As, Ba, and Sb with distance from Owens Valley suggest that dust from southern Owens Valley is being transported at least 400 km to the east. Samples of old dust that accumulated before European settlement are distinctly lower in As, Ba, and Sb abundances relative to modern dust, likely due to modern transport of dust from Owens Valley. Thus, southern Owens Valley appears to be an important, geochemically distinct, point source for regional dust in the southwestern United States.     

Oxygen fugacity and geochemical variations in the martian basalts: implications for martian basalt petrogenesis and the oxidation state of the upper mantle of Mars

The oxygen fugacity of the Dar al Gani 476 martian basalt is determined to be quartz-fayalite-magnetite (QFM) −2.3 ± 0.4 through analysis of olivine, low-Ca pyroxene, and Cr-spinel and is in good agreement with revised results from Fe-Ti oxides that yield QFM −2.5 ± 0.7. This estimate falls within the range of oxygen fugacity for the other martian basalts, QFM −3 to QFM −1. Oxygen fugacity in martian basalts correlates with ⁸⁷Sr/⁸⁶Sr, ¹⁴³Nd/¹⁴⁴Nd, and La/Yb ratios, indicating that the mantle source of the basalts is reduced and that assimilation of crust-like material controls the oxygen fugacity. This allows constraints to be placed on the oxidation state of the martian mantle and on the nature of assimilated crustal material. The assimilated material may be the product of early and extensive hydrothermal alteration of the martian crust, or it may be amphibole- or phlogopite-bearing basaltic rock within the crust. In either case, water may play a significant role in the oxidation of basaltic magmas on Mars, although it may be secondary to assimilation of ferric iron-rich material.     

Effect of adsorbed metal ions on the transport of Zn- and Ni-EDTA complexes in a sand and gravel aquifer

Adsorption, complexation, and dissolution reactions strongly influenced the transport of metal ions complexed with ethylenediaminetetraacetic acid (EDTA) in a predominantly quartz-sand aquifer during two tracer tests conducted under mildly reducing conditions at pH 5.8 to 6.1. In tracer test M89, EDTA complexes of zinc (Zn) and nickel (Ni), along with excess free EDTA, were injected such that the lower portion of the tracer cloud traveled through a region with adsorbed manganese (Mn) and the upper portion of the tracer cloud traveled through a region with adsorbed Zn. In tracer test S89, Ni- and Zn-EDTA complexes, along with excess EDTA complexed with calcium (Ca), were injected into a region with adsorbed Mn. The only discernable chemical reaction between Ni-EDTA and the sediments was a small degree of reversible adsorption leading to minor retardation. In the absence of adsorbed Zn, the injected Zn was displaced from EDTA complexes by iron(III) [Fe(III)] dissolved from the sediments. Displacement of Zn by Fe(III) on EDTA became increasingly thermodynamically favorable with decreasing total EDTA concentration. The reaction was slow compared to the time-scale of transport. Free EDTA rapidly dissolved aluminum (Al) from the sediments, which was subsequently displaced slowly by Fe. In the portion of tracer cloud M89 that traveled through the region contaminated with adsorbed Zn, little displacement of Zn complexed with EDTA was observed, and Al was rapidly displaced from EDTA by Zn desorbed from the sediments, in agreement with equilibrium calculations. In tracer test S89, desorption of Mn dominated over the more thermodynamically favorable dissolution of Al oxyhydroxides. Comparison with results from M89 suggests that dissolution of Al oxyhydroxides in coatings on these sediment grains by Ca-EDTA was rate-limited whereas that by free EDTA reached equilibrium on the time-scale of transport. Rates of desorption are much faster than rates of dissolution of Fe oxyhydroxides from sediment-grain surfaces and, therefore, adsorbed metal ions can strongly influence the speciation of ligands like EDTA in soils and sediments, especially over small temporal and spatial scales.     

Climate change impacts on U.S. Coastal and Marine Ecosystems

Increases in concentrations of greenhouse gases projected for the 21st century are expected to lead to increased mean global air and ocean temperatures. The National Assessment of Potential Consequences of Climate Variability and Change (NAST 2001) was based on a series of regional and sector assessments. This paper is a summary of the coastal and marine resources sector review of potential impacts on shorelines, estuaries, coastal wetlands, coral reefs, and ocean margin ecosystems. The assessment considered the impacts of several key drivers of climate change: sea level change; alterations in precipitation patterns and subsequent delivery of freshwater, nutrients, and sediment; increased ocean temperature; alterations in circulation patterns; changes in frequency and intensity of coastal storms; and increased levels of atmospheric CO₂. Increasing rates of sea-level rise and intensity and frequency of coastal storms and hurricanes over the next decades will increase threats to shorelines, wetlands, and coastal development. Estuarine productivity will change in response to alteration in the timing and amount of freshwater, nutrients, and sediment delivery. Higher water temperatures and changes in freshwater delivery will alter estuarine stratification, residence time, and eutrophication. Increased ocean temperatures are expected to increase coral bleaching and higher CO₂ levels may reduce coral calcification, making it more difficult for corals to recover from other disturbances, and inhibiting poleward shifts. Ocean warming is expected to cause poleward shifts in the ranges of many other organisms, including commercial species, and these shifts may have secondary effects on their predators and prey. Although these potential impacts of climate change and variability will vary from system to system, it is important to recognize that they will be superimposed upon, and in many cases intensify, other ecosystem stresses (pollution, harvesting, habitat destruction, invasive species, land and resource use, extreme natural events), which may lead to more significant consequences.     

Multiple stressors in an estuarine system: Effects of nutrients, trace elements, and trophic complexity on benthic photosynthesis and respiration

The effects of nutrients, trace elements, and trophic complexity on benthic photosynthesis and respriation were studied in the Paxtuxent River estuary near St. Leonard, Maryland. Experiments were conducted over three years (1995–1997) in mesocosms containing riverine sediment and water. The experimental design was 2×2×3 factorial with two levels of nutrients (ambient and + nutrients), two of trace elements (ambient and + trace elements) and three of trophic complexity (plankton, plankton + fish, and plankton + fish + benthos). Trace elements included arsenic (As), copper (Cu), and cadmium (Cd). The experiment was conducted three times in 1995 and 1997 and four times in 1996. In 1995 and 1996, sediments were muddy, while in the final year sediments were sandy. In mesocoms with sandy sediments, nutrient additions increased benthic photosynthesis overall, while trace element additions increased benthic photosynthesis in two of three experimental runs. Benthic photosynthesis in these mesocosms appeared to be related to nitrogen loading. Benthic respiration increased in nutrient and trace element amended mesocosms with sandy sediments, apparently in response to higher benthic photosynthesis. Increasing trophic complexity, particularly the presence of benthic organisms, also increased benthic respiration in mesocosms with sandy sediments. There were no effects of nutrient or trace element additions on benthic photosynthesis and respiration when the sediments were muddy. The lack of consistent responses to nutrient additions was surprising given that benthic respiration rates (and presumably nutrient regeneration) were similar in all three years, regardless of sediment type. Muddy, sediments did not mask, the effects of nutrient addition by supplying more nutrients to benthic microalgae than sandy sediments. In 1996, the presence of filter feeding bivalves increased the relative heterotrophy of sediments, measured as production: respiration. Consistent with increased heterotrophy, effluxes of ammonium and soluble reactive phosphorus from sediments were greater in mesocosms containing benthic organisms. Anthropogenically-induced changes in estuaries, such as loading of nutrients and trace elements or reduced trophic complexity, can have important effects on benthic processes and potentially pelagic processes through feedback mechanisms.     

Method of Data Reduction and Uncertainty Estimation for Platinum-Group Element Data Using Inductively Coupled Plasma-Mass Spectrometry

A data reduction method is described for determining platinum-group element (PGE) abundances by inductively coupled plasma-mass spectrometry (ICP-MS) using external calibration or the method of standard addition. Gravimetric measurement of volumes, the analysis of reference materials and the use of procedural blanks were all used to minimise systematic errors. Internal standards were used to correct for instrument drift. A linear least squares regression model was used to calculate concentrations from drift-corrected counts per second (cps). Furthermore, mathematical manipulations also contribute to the uncertainty estimates of a procedure. Typical uncertainty estimate calculations for ICP-MS data manipulations involve: (1) Carrying standard deviations from the raw cps through the data reduction or (2) calculating a standard deviation from multiple final concentration calculations. It is demonstrated that method 2 may underestimate the uncertainty estimate of the calculated data. Methods 1 and 2 do not typically include an uncertainty estimate component from a regression model. As such models contribute to the uncertainty estimates affecting the calculated data, an uncertainty estimate component from the regression must be included in any final error calculations. Confidence intervals are used to account for uncertainty estimates from the regression model. These confidence intervals are simpler to calculate than uncertainty estimates from method 1, for example. The data reduction and uncertainty estimation method described here addresses problems of reporting PGE data from an article in the literature and addresses both precision and accuracy. The method can be applied to any analytical technique where drift corrections or regression models are used.