Experimental studies of equilibrium iron isotope fractionation in ferric aquo-chloro complexes

Here we compare new experimental studies with theoretical predictions of equilibrium iron isotopic fractionation among aqueous ferric chloride complexes (Fe(H₂O)₆³⁺, FeCl(H₂O)₅²⁺, FeCl₂(H₂O)₄⁺, FeCl₃ (H₂O)₃, and FeCl₄^-), using the Fe-Cl-H₂O system as a simple, easily-modeled example of the larger variety of iron-ligand compounds, such as chlorides, sulfides, simple organic acids, and siderophores. Isotopic fractionation (⁵⁶Fe/⁵⁴Fe) among naturally occuring iron-bearing species at Earth surface temperatures (up to ∼3‰) is usually attributed to redox effects in the environment. However, theoretical modeling of reduced isotopic partition functions among iron-bearing species in solution also predicts fractionations of similar magnitude due to non-redox changes in speciation (i.e., ligand bond strength and coordination number). In the present study, fractionations are measured in a series of low pH ([H⁺] = 5 M) solutions of ferric chloride (total Fe = 0.0749 mol/L) at chlorinities ranging from 0.5 to 5.0 mol/L. Advantage is taken of the unique solubility of FeCl₄^- in immiscible diethyl ether to create a separate spectator phase, used to monitor changing fractionation in the aqueous solution. Δ⁵⁶Fe_aq-eth = δ⁵⁶Fe (total Fe remaining in aqueous phase)−δ⁵⁶Fe (FeCl₄^- in ether phase) is determined for each solution via MC-ICPMS analysis.Both experiments and theoretical calculations of Δ⁵⁶Fe_aq-eth show a downward trend with increasing chlorinity: Δ⁵⁶Fe_aq-eth is greatest at low chlorinity, where FeCl₂(H₂O)₄⁺ is the dominant species, and smallest at high chlorinity where FeCl₃(H₂O)₃ is dominant. The experimental Δ⁵⁶Fe_aq-eth ranges from 0.8‰ at [Cl^-] = 0.5 M to 0.0‰ at [Cl^-] = 5.0 M, a decrease in aqueous-ether fractionation of 0.8‰. This is very close to the theoretically predicted decreases in Δ⁵⁶Fe_aq-eth, which range from 1.0 to 0.7‰, depending on the ab initio model.The rate of isotopic exchange and attainment of equilibrium are shown using spiked reversal experiments in conjunction with the two-phase aqueous-ether system. Equilibrium under the experimental conditions is established within 30 min.The general agreement between theoretical predictions and experimental results points to substantial equilibrium isotopic fractionation among aqueous ferric chloride complexes and a decrease in ⁵⁶Fe/⁵⁴Fe as the Cl^-/Fe³⁺ ion ratio increases. The effects on isotopic fractionation shown by the modeling of this simple iron-ligand system imply that ligands present in an aqueous environment are potentially important drivers of fractionation, are indicative of possible fractionation effects due to other speciation effects (such as iron-sulfide systems or iron bonding with organic ligands), and must be considered when interpreting iron isotope fractionation in the geological record.     

Efficient methods for computing observation impact in 4D-Var data assimilation

This paper presents a practical computational approach to quantify the effect of individual observations in estimating the state of a system. Such a methodology can be used for pruning redundant measurements and for designing future sensor networks. The mathematical approach is based on computing the sensitivity of the analyzed model states (unconstrained optimization solution) with respect to the data. The computational cost is dominated by the solution of a linear system, whose matrix is the Hessian of the cost function, and is only available in operator form. The right-hand side is the gradient of a scalar cost function that quantifies the forecast error of the numerical model. The use of adjoint models to obtain the necessary first- and second-order derivatives is discussed. We study various strategies to accelerate the computation, including matrix-free iterative solvers, preconditioners, and an in-house multigrid solver. Experiments are conducted on both a small-size shallow-water equations model and on a large-scale numerical weather prediction model, in order to illustrate the capabilities of the new methodology.     

Undersampling power-law size distributions: effect on the assessment of extreme natural hazards

The effect of undersampling on estimating the size of extreme natural hazards from historical data is examined. Tests using synthetic catalogs indicate that the tail of an empirical size distribution sampled from a pure Pareto probability distribution can range from having one-to-several unusually large events to appearing depleted, relative to the parent distribution. Both of these effects are artifacts caused by limited catalog length. It is more difficult to diagnose the artificially depleted empirical distributions, since one expects that a pure Pareto distribution is physically limited in some way. Using maximum-likelihood methods and the method of moments, we estimate the power-law exponent and the corner size parameter of tapered Pareto distributions for several natural hazard examples: tsunamis, floods, and earthquakes. Each of these examples has varying catalog lengths and measurement thresholds, relative to the largest event sizes. In many cases where there are only several orders of magnitude between the measurement threshold and the largest events, joint two-parameter estimation techniques are necessary to account for estimation dependence between the power-law scaling exponent and the corner size parameter. Results indicate that whereas the corner size parameter of a tapered Pareto distribution can be estimated, its upper confidence bound cannot be determined and the estimate itself is often unstable with time. Correspondingly, one cannot statistically reject a pure Pareto null hypothesis using natural hazard catalog data. Although physical limits to the hazard source size and attenuation mechanisms from source to site constrain the maximum hazard size, historical data alone often cannot reliably determine the corner size parameter. Probabilistic assessments incorporating theoretical constraints on source size and propagation effects are preferred over deterministic assessments of extreme natural hazards based on historical data.     

Evaluating social vulnerability indicators: criteria and their application to the Social Vulnerability Index

Natural Hazards - As a concept, social vulnerability describes combinations of social, cultural, economic, political, and institutional processes that shape socioeconomic differentials in the...     

Total individual ion activity coefficients of calcium and carbonate in seawater at 25°C and 35%. salinity,and implications to the agreement between apparent and thermodynamic constants of calcite and aragonite

We have calculated the total individual ion activity coefficients of carbonate and calcium, and , in seawater. Using the ratios of stoichiometric and thermodynamic constants of carbonic acid dissociation and total mean activity coefficient data measured in seawater, we have obtained values which differ significantly from those widely accepted in the literature. In seawater at 25°C and 35%. salinity the (molal) values of and are $\text{0.038 ± 0.002}$ and $\text{0.173 ± 0.010}$, respectively. These values of and are independent of liquid junction errors and internally consistent with the value . By defining and on a common scale (), the product is independent of the assigned value of and may be determined directly from thermodynamic measurements in seawater. Using the value and new thermodynamic equilibrium constants for calcite and aragonite, we show that the apparent constants of calcite and aragonite are consistent with the thermodynamic equilibrium constants at 25°C and 35%. salinity. The demonstrated consistency between thermodynamic and apparent constants of calcite and aragonite does not support a hypothesis of stable Mg-calcite coatings on calcite or aragonite surfaces in seawater, and suggests that the calcite critical carbonate ion curve of Broecker and Takahashi (1978, Deep-Sea Research25, 65–95) defines the calcite equilibrium boundary in the oceans, within the uncertainty of the data.     

Trace metals in sediments from the inner continental shelf of the western Beaufort Sea

Increased development and changing climate have enhanced global interest in the Arctic Ocean and adjacent seas. Using a large, 3-yr data base, we have determined the distribution of trace metals in sediments from the inner shelf of the western Beaufort Sea. Metal concentrations in these sediments reflect pristine conditions, consistent with those for most coastal areas in the Arctic and with predictions based on average continental crust. Geographic variations in metal values are primarily related to sediment grainsize distribution. However, grain-size patterns are a complex function of several variables including source areas, ice and water movement, as well as physical and chemical weathering. Future identification of any metal contamination in Beaufort Sea sediments will be simplified by understanding the predictable natural patterns.     

Instrumental neutron activation in geoanalysis

Theoretical aspects of instrumental neutron activation analysis (INAA) are discussed. Various applications of INAA to exploration geochemistry are described demonstrating its capabilities and its limitations.Cost, turnaround time, high sensitivity for certain elements, the non-destructive nature of analysis and its precision and accuracy have combined to make INAA an indispensible method for multielement determination on virtually all geological matrices. Humus, vegetation, heavy minerals, lake bottom sediments as well as rocks and soils comprise the major sample types analyzed by INAA. As many as 50 elements can be determined routinely and easily by INAA.