Effect of Fe on garnet–melt trace element partitioning: experiments in FCMAS and quantification of crystal-chemical controls in natural systems

Garnet–melt trace element partitioning experiments were performed in the system FeO–CaO–MgO–Al₂O₃–SiO₂ (FCMAS) at 3 GPa and 1540°C, aimed specifically at studying the effect of garnet Fe²⁺ content on partition coefficients (D^Grt/Melt). D^Grt/Melt, measured by SIMS, for trivalent elements entering the garnet X-site show a small but significant dependence on garnet almandine content. This dependence is rationalised using the lattice strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal–melt partition coefficients from elastic moduli. Nature 372, 452–454], which describes partitioning of an element i with radius r_i and valency Z in terms of three parameters: the effective radius of the site r₀(Z), the strain-free partition coefficient D₀(Z) for a cation with radius r₀(Z), and the apparent compressibility of the garnet X-site given by its Young's modulus E_X(Z). Combination of these results with data in Fe-free systems [Van Westrenen, W., Blundy, J.D., Wood, B.J., 1999. Crystal-chemical controls on trace element partitioning between garnet and anhydrous silicate melt. Am. Mineral. 84, 838–847] and crystal structure data for spessartine, andradite, and uvarovite, leads to the following equations for r₀(3+) and E_X(3+) as a function of garnet composition (X) and pressure (P):

r₀(3+) [Å]=0.930X_Py+0.993X_Gr+0.916X_Alm+0.946X_Spes+1.05(X_And+X_Uv)−0.005(P [GPa]−3.0)(±0.005 Å)

E_X(3+) [GPa]=3.5×10¹²(1.38+r₀(3+) [Å])^−26.7(±30 GPa)

Accuracy of these equations is shown by application to the existing garnet–melt partitioning database, covering a wide range of P and T conditions (1.8 GPa<P<5.0 GPa; 975°C<T<1640°C). D^Grt/Melt for all 3+ elements entering the X-site (REE, Sc and Y) are predicted to within 10–40% at given P, T, and X, when D^Grt/Melt for just one of these elements is known. In the absence of such knowledge, relative element fractionation (e.g. D_Sm^Grt/Melt/D_Nd^Grt/Melt) can be predicted. As an example, we predict that during partial melting of garnet peridotite, group A eclogite, and garnet pyroxenite, r₀(3+) for garnets ranges from 0.939±0.005 to 0.953±0.009 Å. These values are consistently smaller than the ionic radius of the heaviest REE, Lu. The above equations quantify the crystal-chemical controls on garnet–melt partitioning for the REE, Y and Sc. As such, they represent a major advance en route to predicting D^Grt/Melt for these elements as a function of P, T and X.     

The effect of Ca-Tschermaks component on trace element partitioning between clinopyroxene and silicate melt

We have studied the influence of Ca-Tschermaks (Calcium Tschermaks or CaTs) content of clinopyroxene on the partitioning of trace elements between this phase and silicate melt at fixed temperature and pressure. Ion probe analyses of experiments carried out in the system Na₂O–CaO–MgO–Al₂O₃–SiO₂, at 0.1 MPa and 1218°C, produced crystal-melt partition coefficients (D) of 36 trace elements (Li, Cl, Sc, Ti, V, Cr, Fe, Co, Ge, Sr, Y, Zr, Nb, Mo, Ru, Rh, In, Sn, Sb, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Ta and W), for clinopyroxene compositions between 10 and 32 mol% CaTs. Partition coefficients for 2+ to 5+ cations show, for each charge, a near parabolic dependence of log D on ionic radius of the substituting cation, for partitioning into both the M1 and M2 sites of clinopyroxene. Fitting the results to the elastic strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal-melt partition coefficients from elastic moduli. Nature 372, 452–454] we obtain results for the strain-free partition coefficients of theoretical cations (D₀), with site radius r₀, and for the site's Young's Modulus (E).

In agreement with earlier data our results show that increasing ^ivAl concentration in cpx is matched by increasing D, E_M1, E_M2 and D₀ for tri-, tetra- and pentavalent cations. The degree of fractionation between chemically similar elements (i.e. Ta/Nb, Zr/Hf) also increases. In contrast, D values for mono-, di- and hexavalent cations decrease with increasing ^ivAl in the cpx. The large suite of trace elements used has allowed us to study the effects of cation charge on D₀, r₀ and E. We have found that D₀ and r₀ decrease with increasing cation charge, e.g. r₀=0.66 Å for 4+ cations and 0.59 Å for 5+ cations substituting into M1. Values of E_M1 and E_M2 increase with cation charge as well as with increasing ^ivAl content. The increase in E_M2 is linear and close to the trend set by Hazen and Finger [Hazen, R.M., Finger, L.W., 1979. Bulk modulus-volume relationship for cation–anion polyhedra. J. Geophys. Res. 84 (10) 6723–6728] for oxides. E_M1 values are much higher and do not fit the trend predicted by the Hazen and Finger relationship.     

The influence of Cr on the garnet–spinel transition in the Earth''s mantle: experiments in the system MgO–Cr2O3–SiO2 and thermodynamic modelling

The position of the transition from spinel peridotite to garnet peridotite in a simplified chemical composition has been determined experimentally at high pressures and high temperatures. The univariant reaction MgCr₂O₄+2Mg₂Si₂O₆=Mg₃Cr₂Si₃O₁₂+Mg₂SiO₄, has a negative slope in P–T space between 1200 °C and 1600 °C. The experimental results, combined with assessed thermodynamic data for MgCr₂O₄, MgSiO₃ and Mg₂SiO₄ give the entropy and enthalpy of formation of knorringite garnet (Mg₃Cr₂Si₃O₁₂). Thermodynamic calculations in simplified chemical compositions indicate that Cr shifts the garnet-in reaction to much higher pressures than previously anticipated. Moreover, in Cr-bearing systems a pressure–temperature field exists where garnet and spinel coexist. The width of this divariant field strongly depends on the Cr/(Cr+Al) of the system.     

下垫面对WRF模式模拟黑河流域区域气候精度影响研究

利用高精度的土地覆盖、土壤质地类型和地形高度值替换了天气研究和预报模式Weather Research and Forecasting Model(WRF)中的相关数据,通过数值模式试验检验了下垫面数据对WRF模拟精度的影响。同时,通过与黑河综合遥感联合试验中7个测站观测数据的比较,以平均误差、均方根误差和相关系数为指标,分析了WRF模式下垫面数据改变对近地表气象要素的模拟精度的影响。结果表明:(1)WRF模式本身的地形高度信息在黑河流域上游地区有较大误差,造成了一定的模拟误差。而使用高精度的下垫面数据可以提高WRF模式在黑河流域上游复杂区域的模拟能力;(2)2m气温除了随地形高度递减外,还受土壤质地和土地覆盖小幅度影响,而且进行地形订正后的2m气温与2m湿度的模拟在下垫面为水体的区域对比强烈,因此为模式提供准确的水体分布信息也至关重要;(3)2m气温和湿度等要素的模拟差异值与地形高度资料的差异呈负相关,而降雨量的差异与地形高度差异呈微弱的正相关,与土壤质地差异和土地覆盖差异的相关性也比较弱。     

生态水文研究前沿问题及生态水文观测试验

自1987年Ingram HAP提出生态水文学概念以来,生态水文学得到了快速的发展。2007年,“生态水文学与环境可持续性”已成为UNESCO/IHP第7阶段计划的主题之一。回顾了生态水文学概念的变迁,综述了有关术语及其科学内涵,简要分析了水循环与生物地球化学循环、水文与生态系统相互作用、水文过程与生命过程耦合、绿水及其生态作用等方面关注的科学问题及其观测研究进展;介绍了黑河流域生态水文学研究的初步认识。     

Sea-surface wave growth under extraterrestrial atmospheres: Preliminary wind tunnel experiments with application to Mars and Titan

We describe for the first time the generation and measurement of capillary waves in a water surface in a wind tunnel running with air at pressures of 15-1000 mbar. These experiments suggest a stronger dependence of wave generation on atmospheric density than the simple proportionality that might be expected from energy transfer arguments. Additionally, airflow over a nonaqueous fluid (kerosene) was found to produce waves of higher amplitude than for water under the same conditions. These preliminary results may indicate different efficiencies of wave generation on other planets, for which empirical terrestrial relations therefore do not apply, and thus may have a bearing on the lack of strong shoreline features on Mars and the possibility of specular glints from hydrocarbon lakes on Titan.     

实验流水地貌学研究的回顾与展望

本文简述了地貌实验研究的概况,罗列了我国主要的地貌实验室、站,叙述了我国30多年来,河型造床试验、边界条件及构造运动对河型发育影响实验,应用河流地貌实验、试验理论和技术方面的进展情况,指出实验研究中存在的问题,在回顾基础上,展望了流水地貌实验研究的未来发展趋势。     

Sloshing in a three-dimensional rectangular tank: Numerical simulation and experimental validation

Pressure variations and three-dimensional effects on liquid sloshing loads in a moving partially filled rectangular tank have been carried out numerically and experimentally. A numerical algorithm based on the volume of fluid (VOF) technique is used to study the non-linear behavior and damping characteristics of liquid sloshing. A moving coordinate system is used to include the non-linearity and avoid the complex boundary conditions of moving walls. The numerical model solves the complete Navier–Stokes equations in primitive variables by using of the finite difference approximations. In order to mitigate a series of discrete impacts, the signal computed is averaged over several time steps. In order to assess the accuracy of the method used, computations are compared with the experimental results. Several configurations of both baffled and unbaffled tanks are studied. Comparisons show good agreement for both impact and non- impact type slosh loads in the cases investigated.     

Physical behavior of the SEEDS iron-fertilized patch by sulphur hexafluoride tracer release

The first iron (Fe) – fertilization experiment in the western North Pacific was carried out using SF₆ to trace the Fe-fertilized water mass. A solution in 10,800 liters of seawater of 350 kg of Fe and 0.48 M of SF₆ tracer was released into the mixed layer over a 8 × 10 km area. On the first underway transects through the patch after the Fe release, we observed a significant increase of dissolved Fe (ave. 2.89 nM). The fertilized patch was traced for 14 days by on-board SF₆ analysis. A Lagrangian frame of reference was maintained by the use of a drogued GPS buoy released at the center of the patch. The patch moved westward at a rate of 6.8 km d⁻¹. Mixed layer depth increased from 8.5 to 15 m during the experiment. Horizontal diffusivity was determined by the change of SF₆ concentration in the patch. The horizontal diffusivity increased during the experiment. We evaluate here the fate of Fe in a Fe-fertilized patch using the dilution rate determined from sulphur hexafluoride (SF₆) concentration. Dissolved Fe concentrations subsequently decreased rapidly to 0.15 nM on Day 13. However, the dissolved Fe half-life of 43 h was relatively longer than in previous Fe-enrichment studies, and we observed a larger increase of the centric diatom standing stock and corresponding drawdown of macro-nutrients and carbon dioxide than in the previous studies. The most important reason for the larger response was the phytoplankton species in the western North Pacific. In addition, the smaller diffusivity and shallower mixed layer were effective to sustain the higher dissolved Fe concentration compared to previous experiments. This might be one reason for the larger response of diatoms in SEEDS.     

铰接式浮标性能研究

铰接式浮标作为一种新型助航标志越来越受到海监部门的欢迎。铰接式浮标具有良好的水动力性能，在风流中静倾角及在波浪中摇摆角都较小，而且定位准确。本文对铰接式浮标的结构及其水动力性能计算方法作了介绍，并在国家重点实验室──海洋工程实验室作了系列试验，得出一些有益的结果，为优化设计提供了依据。