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.     

SIMS determination of trace element partition coefficients between garnet, clinopyroxene and hydrous basaltic liquids at 2–7.5 GPa and 1080–1200°C

Trace element partition coefficients (D's) for up to 13 REE, Nb, Ta, Zr, Hf, Sr and Y have been determined by SIMS analysis of seven garnets, four clinopyroxenes, one orthopyroxene and one phlogopite crystallized from an undoped basanite and a lightly doped (200 ppm Nb, Ta and Hf) quartz tholeiite. Experiments were conducted at 2–7.5 GPa, achieving near-liquidus crystallization at relatively low temperatures of 1080–1200°C under strongly hydrous conditions (5–27 wt.% added water). Garnet and pyroxene D_REE show a parabolic pattern when plotted against ionic radius, and conform closely to 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). Comparison, at constant pressure, between hydrous and anhydrous values of the strain-free partition coefficient (D₀) for the large cation sites of garnet and clinopyroxene reveals the relative importance of temperature and melt water content on partitioning. In the case of garnet, the effect of lower temperature, which serves to increase D₀, and higher water content, which serves to decrease D₀, counteract each other to the extent that water has little effect on garnet–melt D₀ values. In contrast, the effect of water on clinopyroxene–melt D₀ overwhelms the effect of temperature, such that D₀ is significantly lower under hydrous conditions. For both minerals, however, the lower temperature of the hydrous experiments tends to tighten the partitioning parabolas, increasing fractionation of light from heavy REE compared to anhydrous experiments.

Three sets of near-liquidus clinopyroxene–garnet two-mineral D values increase the range of published experimental determinations, but show significant differences from natural two-mineral D's determined for subsolidus mineral pairs. Similar behaviour is observed for the first experimental data for orthopyroxene–clinopyroxene two-mineral D's when compared with natural data. These differences are in large part of a consequence of the subsolidus equilibration temperatures and compositions of natural mineral pairs. Great care should therefore be taken when using natural mineral–mineral partition coefficients to interpret magmatic processes.

The new data for strongly hydrous compositions suggest that fractionation of Zr–Hf–Sm by garnet decreases with increasing depth. Thus, melts leaving a garnet-dominated residuum at depths of about 200 km or greater may preserve source Zr/Hf and Hf/Sm. This contrasts with melting at shallower depths where both garnet and clinopyroxene will cause Zr–Hf–Sm fractionation. Also, at shallower depths, clinopyroxene-dominated fractionation may produce a positive Sr spike in melts from spinel lherzolite, but for garnet lherzolite melting, no Sr spike will result. Conversely, clinopyroxene megacrysts with negative Sr spikes may crystallize from magmas without anomalous Sr contents when plotted on mantle compatibility diagrams. Because the characteristics of strongly hydrous silicate melt and solute-rich aqueous fluid converge at high pressure, the hydrous data presented here are particularly pertinent to modelling processes in subduction zones, where aqueous fluids may have an important metasomatic role.     

Partitioning coefficients between olivine and silicate melts

Variation of Nernst partition coefficients (D) between olivine and silicate melts cannot be neglected when modeling partial melting and fractional crystallization. Published natural and experimental ^{olivine/liquid}D data were examined for covariation with pressure, temperature, olivine forsterite content, and melt SiO₂, H₂O, MgO and MgO/MgO + FeO^total. Values of ^{olivine/liquid}D generally increase with decreasing temperature and melt MgO content, and with increasing melt SiO₂ content, but generally show poor correlations with other variables. Multi-element ^{olivine/liquid}D profiles calculated from regressions of D REE–Sc–Y vs. melt MgO content are compared to results of the Lattice Strain Model to link melt MgO and: D₀ (the strain compensated partition coefficient), E_M³⁺ (Young's Modulus), and r₀ (the size of the M site). Ln D₀ varies linearly with Ln MgO in the melt; E_M³⁺ varies linearly with melt MgO, with a dog-leg at ca. 1.5% MgO; and r₀ remains constant at 0.807 Å. These equations are then used to calculate ^{olivine/liquid}D for these elements using the Lattice Strain Model. These empirical parameterizations of ^{olivine/liquid}D variations yield results comparable to experimental or natural partitioning data, and can easily be integrated into existing trace element modeling algorithms. The ^{olivine/liquid}D data suggest that basaltic melts in equilibrium with pure olivine may acquire small negative Ta–Hf–Zr–Ti anomalies, but that negative Nb anomalies are unlikely to develop. Misfits between results of the Lattice Strain Model and most light rare earth and large ion lithophile partitioning data suggest that kinetic effects may limit the lower value of D for extremely incompatible elements in natural situations characterized by high cooling/crystallization rates.     

Strontium diffusion in sanidine and albite, and general comments on strontium diffusion in alkali feldspars

Strontium chemical diffusion has been measured in albite and sanidine under dry, 1 atm, and QFM buffered conditions. Strontium oxide-aluminosilicate powdered sources were used to introduce the diffusant and Rutherford Backscattering Spectroscopy (RBS) used to measure diffusion profiles. For the 1 atm experiments, the following Arrhenius relations were obtained:

Sanidine (Or₆₁), temperature range 725–1075°C, diffusion normal to (001): D=8.4 exp(−450±13 kJ mol⁻¹/RT) m²s⁻¹. Albite (Or₁), temperature range 675–1025°C, diffusion normal to (001): D=2.9 × exp(−224±11 kJ mol⁻¹/RT) m²s⁻¹.

The alkali feldspars in this and earlier work display a broad range of activation energies for Sr diffusion, which may be a consequence of the thermodynamic non-ideality of the alkali feldspar system and/or the mixed alkali effect.     

A palaeomagnetic study of some basaltoids and ores from the Pontic Ranges, Northern Turkey: Palaeogeographic implications

The palaeomagnetism of basic rocks and sulphide ores has been studied in the Küre area, Pontic Ranges, Turkey. Progressive alternating-field demagnetization revealed a characteristic remanent magnetization in all investigated rock types except a dacite. The following virtual geomagnetic poles were obtained:

Basalt and quartz diabase (oldest): D = 59°, I = +66°, ₉₅ = 4.8, pole 49°N, 93°E. Diabase: D = 210°, I = −15°, ₉₅ = 15.0, pole 47°N, 167°E. Massive sulphide ores: D = 107°, I = +63°, ₉₅ = 8.7, pole 18°N, 80°E. Peridotite: D = 131°, I = +54°, ₉₅ = 10.9, pole 2°S, 72°E. Amphibolitized diabase (youngest): D = 293°, I = +59°, ₉₅ = 12.6, pole 40°S, 145°E.

The longidutinal difference in pole positions between the oldest and the youngest rocks is interpreted as being due to a post-Permian counterclockwise rotation of the studied region in relation to the African continent. In addition, there are indications of local rotational movements within the Küre area.     

基于多重分形理论的多年冻土区高寒草甸退化过程中土壤粒径分析

为阐明青藏高原多年冻土区高寒草甸退化过程中土壤粒径分布(PSD)非均匀性和异质性的变化特征, 在青藏高原长江源区, 根据高寒草甸的退化梯度, 选取了未退化区域、 轻度退化区域、 中度退化区域、 重度退化区域和极重度退化区域, 测定了高寒草甸退化过程中土壤的粒径分布、 饱和导水率、 孔隙度与有机质含量. 运用多重分形理论, 并结合土壤颗粒分布与土壤理化特性等参数的相关性进行分析, 为高寒草甸退化对长江源高寒土壤性质变化的影响的定量研究提供一种精确的分析方法. 结果表明: 随着青藏高原多年冻土区高寒草甸退化程度的增加, 土壤颗粒呈粗粒化趋势, 多重分形参数中容量维数(D₀)随之增大, 表征PSD宽度随之增大; 信息维数(D₁)、 信息维数/容量维数(D₁/D₀)、 关联维数(D₂)、 奇异谱宽(Δα)可从不同角度反映的土壤PSD的非均匀性与局部异质性随着高寒草甸退化有先增大后减小的趋势, 中度退化区域的土壤PSD不均匀性最大. 研究发现, 研究区土壤多重分形参数与细砂含量、 土壤的孔隙度、 有机质含量具有较明显的相关性. 多重分形参数能准确描述高寒草甸退化过程中土壤粒径分布的细微差别, 可作为反映土壤性质的潜在指标.     

Status report on stability of K-rich phases at mantle conditions

Experimental research on K-rich phases and observations from diamond inclusions, UHP metamorphic rocks, and xenoliths provide insights about the hosts for potassium at mantle conditions. K-rich clinopyroxene (Kcpx–KM³⁺Si₂O₆) can be an important component in clinopyroxenes at P>4 GPa, dependent upon coexisting K-bearing phases (solid or liquid) but not, apparently, upon temperature. Maximum Kcpx content can reach 25 mol%, with 17 mol% the highest reported in nature. Partitioning ^(K)D_(cpx/liquid) above 7 GPa=0.1–0.2 require ultrapotassic liquids to form highly potassic cpx or critical solid reactions, e.g., between Kspar and Di. Phlogopite can be stable to about 8 GPa at 1250 °C where either amphibole or liquid forms. When fluorine is present, it generally increases in Phl upon increasing P (and probably T) to about 6 GPa, but reactions forming amphibole and/or KMgF₃ limit F content between 6 and 8 GPa. The perovskite KMgF₃ is stable up to 10 GPa and 1400 °C as subsolidus breakdown products of phlogopite upon increasing P. ^(M4)K-substituted potassic richterite (ideally K(KCa)Mg₅Si₈O₂₂(OH,F)₂) is produced in K-rich peridotites above 6 GPa and in Di+Phl from 6 to 13 GPa. K content of amphibole is positively correlated with P; Al and F content decrease with P. In the system 1Kspar+1H₂O K-cymrite (hydrous hexasanidine–KAlSi₃O₈·nH₂O–Kcym) is stable from 2.5 GPa at 400 to 1200 °C and 9 GPa; Kcym can be a supersolidus phase. Formation of Kcym is sensitive to water content, not forming within experiments with H₂O2O>Kspar. Phase X, a potassium di-magnesium acid disilicate ((K_1−x−n)₂(Mg_1−nM_n³⁺)₂Si₂O₇H_2x), forms in mafic compositions at T=1150–1400 °C and P=9–17 GPa and is a potential host for K and H₂O at mantle conditions with a low-T geotherm or in subducting slabs. The composition of phase-X is not fixed but actually represents a solid solution in the stoichiometries □₂Mg₂Si₂O₇H₂–(K□)Mg₂Si₂O₇H–K₂Mg₂Si₂O₇ (□=vacancy), apparently stable only near the central composition. K-hollandite, KAlSi₃O₈, is possibly the most important K-rich phase at very high pressure, as it appears to be stable to conditions near the core–mantle boundary, 95 GPa and 2300 °C. Other K-rich phases are considered.     

Angstrom公式参数对ET0的影响及FAO建议值适用性评价

作为计算太阳总辐射(R_s)的主要公式,Angstrom公式参数(a、b)的合理取值是计算参照腾发量(ET₀)的重要前提。针对FAO所提出的a、b建议值(a=0.25、b=0.5)在中国无辐射观测资料地区被大量使用,而其合理性尚未得到系统评价的情况,基于中国104个地面站的观测数据,在逐月时间尺度上,讨论了a、b变化对ET₀的影响,分析了a、b的地区分布规律,评价了FAO建议值所导致的ET₀计算误差,进而阐明了该建议值在中国7个区域的适用性。提出了无辐射资料情况下a、b的地区综合取值方法。主要研究结论是:①参数a、b偏差对ET₀的计算有重要影响,在中国无资料地区采用FAO建议值将导致较大的太阳总辐射(R_s)和ET₀计算误差。②大多数站点,a的率定值较FAO建议值明显偏小,而b的率定值明显偏大。新疆地区和华南地区a、b率定值分布比较集中,而在其它区域比较分散。③FAO建议参数值在东北、西北和新疆3个区域计算ET₀的适用性较好,而在西南和华南两个区域的适用性很差,计算的ET₀偏高较大。④提出的地区综合取值方法,能使R_s和ET₀的计算精度较FAO建议值显著提高。     

基于区域农业用水量的干旱重现期计算方法

为研究实际水利条件下农业干旱的发生规律,简化农业干旱事件的评估方法,提出基于区域农业用水量的干旱重现期计算方法。通过构建农业用水量距平百分率干旱指标W_A,在基于降雨量距平百分率干旱指标P_A识别干旱事件的基础上,提取W_A干旱指标下的干旱历时和干旱烈度特征变量,并根据以P_A为干旱指标的干旱烈度频率分布曲线F_S（x）和干旱历时频率分布曲线F_D（x）,运用Copula的简化方法计算基于W_A的干旱事件重现期T,最后结合基于P_A的干旱事件重现期T₀,回归分析出T与T₀间关系的计算公式。选取干旱灾害影响严重的亳州市为实证区域开展应用研究,计算得到1975-2007年各场干旱事件的T₀和T以及T₀与T的经验关系式。结果表明:T比T₀更合理地反映区域农业实际受旱状况,重现期T₀和T间存在高度的相关关系,采用T的回归方程可简化计算考虑区域实际抗旱能力下的干旱事件重现期,在区域防旱减灾实践中具有推广应用价值。     

Application of in situ-produced cosmogenic Be and Al to the study of lateritic soil development in tropical forest: theory and examples from Cameroon and Gabon

Depth profiles of in situ-produced cosmogenic nuclides, including ¹⁰Be (T_1/2=1.5×10⁶ years) and ²⁶Al (T_1/20.73×10⁶ years), in the upper few meters of the Earth's crust may be used to study surficial processes, quantifying denudation and burial rates and elucidating mechanisms involved in landform evolution and soil formations. In this paper, we discuss the fundamentals of the method and apply it to two lateritic sequences located in African tropical forests.     

Mineralogy and petrology of cretaceous subsurface lamproite sills, southeastern Kansas, USA

Cores and cuttings of lamproite sills and host sedimentary country rocks in southeastern Kansas from up to 312 m depth were analyzed for major elements in whole rocks and minerals, certain trace elements in whole rocks (including the REE) and Sr isotopic composition of the whole rocks. The lamproites are ultrapotassic (K₂O/Na₂O = 2.0–19.9), alkalic [molecular (K₂O/Na₂O)/Al₂O₃ = 1.3-2.8], enriched in mantle-incompatible elements (light REE, Ba, Rb, Sr, Th, Hf, Ta) and have nearly homogeneous initial Sr isotopic compositions (0.707764-0.708114).

These lamproites could have formed by variable degrees of partial melting of harzburgite country rock and cross-cutting veins composed of phlogopite, K-Ti richterite, titanite, diopside, K-Ti silicates, or K-Ba-phosphate under high H₂O/CO₂ ratios and reducing conditions. Variability in melting of veins and wall rock and variable composition of the metasomatized veins could explain the significantly different composition of the Kansas lamproites.

Least squares fractionation models preclude the derivation of the Kansas lamproites by fractional crystallization from magmas similar in composition to higher silica phlogopite-sanidine lamproites some believe to be primary lamproite melts found elsewhere. In all but one case, least squares fractionation models also preclude the derivation of magmas similar in composition to any of the Kansas lamproites from one another. A magma similar in composition to the average composition of the higher SiO₂ Ecco Ranch lamproite (237.5–247.5 m depth) could, however, have marginally crystallized about 12% richterite, 12% sanidine, 7% diopside and 6% phlogopite to produce the average composition of the Guess lamproite (305–312 m depth).

Lamproite from the Ecco Ranch core is internally fractionated in K₂O, Al₂O₃, Ba, MgO, Fe₂O₃, Co and Cr most likely by crystal accumulation-removal of ferromagnesian minerals and sanidine. In contrast, the Guess core (305–312 m depth) has little fractionation throughout most of the sill except in several narrow zones. Lamproite in the Guess core has large enrichments in TiO₂, Ba, REE, Th, Ta and Sc and depletions in MgO, Cr, Co and Rb possibly concentrated in these narrow zones during the last dregs of crystallization of this magma.

The Ecco Ranch sill did not show any evidence of loss of volatiles or soluble elements into the country rock. This contrasts to the previously studied, shallow Silver City lamproite which did apparently lose H₂O-rich fluid to the country rock. Perhaps a greater confining pressure and lesser amount of H₂O-rich fluid prevented it from escaping.     

Regional time-predictable modeling in the Hindukush–Pamir–Himalayas region

The seismic characteristic of Hindukush–Pamir–Himalaya (HPH) and its vicinity is very peculiar and has experienced many widely distributed large earthquakes. Recent work on the time-dependent seismicity in the Hindukush–Pamir–Himalayas is mainly based on the so-called “regional time-predictable model”, which is expressed by the relation log T=cM_p+a, where T is the inter-event time between two successive main shocks of a region and M_p is the magnitude of the preceded main shock. Parameter a is a function of the magnitude of the minimum earthquake considered and of the tectonic loading and c is positive (0.3) constant. In 90% of the cases with sufficient data, parameter c was found to be positive, which strongly supports the validity of the model. In the present study, a different approach, which assumes no prior regionalization of the area, is attempted to check the validity of the model. Nine seismic sources were defined within the considered region and the inter-event time of strong shallow main shock were determined and used for each source in an attempt at long-term prediction, which show the clustering and occurrence of at least three earthquakes of magnitude 5.5≤M_s≤7.5 giving two repeat times, satisfying the necessary and sufficient conditions of time-predictable model (TP model). Further, using the global applicability of the regional time- and magnitude-predictable model, the following relations have been obtained: log T_t=0.19 M_min+0.52M_p+0.29 log m₀−10.63 and M_f=1.31M_min−0.60M_p−0.72 log m₀+21.01, where T_t is the inter-event time, measured in years; M_min the surface wave magnitude of the smallest main shock considered; M_p the magnitude of preceding main shock; M_f the magnitude of the following main shock; and m₀ the moment rate in each source per year.

These relations may be used for seismic hazard assessment in the region. Based on these relations and taking into account the time of occurrence and the magnitude of the last main shock in each seismogenic source, time-dependent conditional probabilities for the occurrence of the next large (M_s≥5.5) shallow main shocks during the next 20 years as well as the magnitudes of the expected main shocks are determined.     

Marcasite precipitation from hydrothermal solutions

Pyrite and marcasite were precipitated by both slow addition of aqueous Fe²⁺ and SiO₃²⁻ to an H₂S solution and by mixing aqueous Fe²⁺ and Na₂S₄ solutions at 75°C. H₂S₂ or HS₂⁻ and H₂S₄ or HS₄⁻ were formed in the S₂O₃²⁻ and Na₂S₄ experiments, respectively. Marcasite formed at pH < pK₁ of the polysulfide species present (for H₂S₂, pK₁ = 5.0; for H₂S₄, pK₁ = 3.8 at 25°C). Marcasite forms when the neutral sulfane is the dominant polysulfide, whereas pyrite forms when mono-or divalent polysulfides are dominant. In natural solutions where H₂S₂ and HS₂ are likely to be the dominant polysulfides, marcasite will form only below pH 5 at all temperatures.

The pH-dependent precipitation of pyrite and marcasite may be caused by electrostatic interactions between polysulfide species and pyrite or marcasite growth surfaces: the protonated ends of H₂S₂ and HS₂ are repelled from pyrite growth sites but not from marcasite growth sites. The negative ions HS₂ and S₂²⁻ are strongly attracted to the positive pyrite growth sites. Masking of 1π_g^* electrons in the S₂ group by the protons makes HS₂ and H₂S₂ isoelectronic with AsS²⁻ and As₂²⁻, respectively ( et al., 1981). Thus, the loellingitederivative structure (marcasite) results when both ends of the polysulfide are protonated.

Marcasite occurs abundantly only for conditions below pH 5 and where H₂S₂ was formed near the site of deposition by either partial oxidation of aqueous H₂S by O₂ or by the reaction of higher oxidation state sulfur species that are reactive with H₂S at the conditions of formation e.g., S₂O₃²⁻ but not SO₄²⁻. The temperature of formation of natural marcasite may be as high as 240°C ( and , 1985), but preservation on a multimillion-year scale seems to require post-depositional temperatures of below about 160°C ( , 1973; and , 1985).     

利用数字岩心技术评价含黏土砂岩导电模型

黏土附加导电性使得岩石导电机理复杂化,影响了测井解释中饱和度计算的准确性。随着石油勘探的发展,测井解释研究者针对特定地区和储层条件的饱和度求解提出了多种导电模型。本文基于贝雷砂岩的数字岩心,构建多种含黏土的砂岩数字岩心,利用有限元方法模拟岩心饱和水时的电导率C₀,并将数值模拟C₀与5种模型计算C₀对比,分析各种导电模型的适用情况。结果表明,黏土体积分数越小,黏土的阳离子交换容量越小,岩石孔隙度越大,公式计算C₀越接近数值模拟C₀。Doll公式和Indonesia公式计算C₀与数值模拟C₀接近,基本适用于不同孔隙度、黏土体积分数和阳离子交换容量的情况。     

Weathering effects on the strength and deformational behaviour of crystalline rocks under uniaxial compression state

Uniaxial compression tests were performed on different categories of weathering of three lithological units: Malanjkhand granite; Nagpur basalt; and Delhi quartzite, occurring in central and northern parts of India. The deformational behaviour is studied in terms of variation in tangent modulus (E_t50) and initial modulus (E_i) due to weathering. The power relationship between uniaxial compressive strength (σ_c) and E_t50 shows strong correspondence for weathering sequence of common rock types. This relationship has been established by regression analysis and significant correlation parameter (coefficient of determination, r²=0.87) for crystalline rocks. It is shown that there is a systematic decrease in stiffness ratio, that is, ratio of tangent modulus and uniaxial compressive strength with increased weathering state. Comparison of E_t50 and E_i values has shown that E_t50 decreases more gradually than E_i, and reduction is more drastic for E_i values with an increased degree of weathering in all the three rock types. The mode of failure has been found to be influenced by weathering extent in rocks. A brief account is given of the intrinsic characteristics of fresh and weathered rocks and mineralogical changes produced by weathering investigated quantitatively. Correlation drawn between the petrographical and mechanical indices has shown that mechanical properties are apparently dependent on the intrinsic characteristics of weathered rocks.     

荒漠绿洲区人工梭梭林土壤水分空间异质性的定量研究

利用12×12m²样地中1×1m²、0～100cm剖面的土壤水分调查数据,采用地统计学原理与方法,研究了人工梭梭林(Haloxylon ammodendron)在栽植20a后的土壤水分格局的空间异质性程度、异质性组成、尺度以及与梭梭生长的关系.结果表明:人工梭梭林土壤水分空间异质性的96%～88%是由空间自相关因素引起的,随机因素起的作用较小.除60～80cm土层土壤水分的块金值与基台值比值较高(C₀/(C₀+C)=0.5),其它各层都较小(0.04～0.12),变程为1.57～2.97m.在较小(<2m)和较大(>8m)的尺度上,土壤水分的空间相关性较强.沿垂直剖面土壤含水量差异显著,10～20cm土层含水量最高(2.82%),其它各层较小(1.30%～1.67%).     

钻孔岩心黏土混浊水电导率、黄铁矿、pH相关性分析及其在古沉积环境复原的应用:以渤海湾西岸平原DC01孔为例

为探讨复原海岸平原沉积环境和精确划分海陆相地层的有效方法,对渤海湾西岸海河以南平原深达30 m的DC01孔岩心以约20 cm间距取样,测试黏土混浊水电导率（EC）、黄铁矿质量分数（w（FeS₂））和pH值,并开展相关分析研究。结果显示,EC与w（FeS₂）正相关（相关系数r为0.47）,而EC与pH以及w（FeS₂）与pH均负相关（r分别为-0.43和-0.52）。根据EC值的大小将DC01孔自下至上分成5个带（Ⅰ—Ⅴ）,进一步分别对5个带的EC、w（FeS₂）及pH做相关分析,结果显示,Ⅱ带相关性最为突出：EC与w（FeS₂）的相关系数为0.83,极强正相关;EC与pH的相关系数为-0.77,强负相关;w（FeS₂）与pH的相关系数为-0.45,呈中等程度负相关。根据Ⅱ带的EC和w（FeS₂）极强正相关、且EC和w（FeS₂）均明显偏高以及多发育黑色泥炭和腐殖质黏土等特征,推断Ⅱ带为受海水影响的盐沼环境;结合AMS¹⁴C测年结果,推测8 260~7 470 cal.B.P.期间,海水曾经影响到DC01孔的位置。另外,Ⅳ带的EC值也偏高,但w（FeS₂）却较低,其EC与w（FeS₂）的相关系数为-0.03,不呈正相关关系,其原因是因为接近地表的上部沉积物因淋溶作用,在Ⅳ带形成淀积层导致EC异常偏高;最上层的Ⅴ带为淋溶层,这样就导致钻孔上部（深6.7 m以上）Ⅳ带和Ⅴ带EC与w（FeS₂）的相关性出现异常。因此在讨论沉积物上部沉积环境时,将EC与w（FeS₂）综合分析可以更精确地划分钻孔岩心的沉积环境。     

Significance of Nb/Ta as an indicator of geochemical processes in the crust-mantle system

A mantle value of 17.5 for Nb/Ta appears well established; less well established are crustal values of 11–12, although it appears that Nb/Ta for crustal-derived melts is less than mantle Nb/Ta, demonstrating fractionation of these two elements during crustal evolution, and suggesting that Nb/Ta variation may be indicative of a particular chemical process within the crust-mantle system.

Experimental studies on silicate and carbonatitic liquids at high pressure indicate that, although silicate minerals such as garnet, amphibole and clinopyroxene do fractionate Nb and Ta, the partition coefficients (D's) for both elements are very low. Thus involvement of these minerals may explain relatively small changes in Nb/Ta, but appears inadequate to explain the crust-mantle variation. However, high-quality data for Nb, Ta may be used to provide information on mantle melting or metasomatic processes (e.g., amphibole in the source region decreases Nb/Ta in derived melts, while carbonatitic metasomatism will increase Nb/Ta in affected mantle). Titanate minerals have high D's for Nb and Ta, and do fractionate these elements (e.g., D_Nb/D_Ta rutile/liquid of 0.5–0.8), and their involvement in crystal fractionation would increase Nb/Ta in derivative liquids. In contrast, D_Nb/D_Ta for rutile/fluid is 1.25, so that rocks affected by fluid equilibrated with residual rutile will show a decrease in Nb/Ta

Some Archaean gneisses appear to have high Nb/Ta, and may be a complementary component to that part of the crust which has a relatively low Nb/Ta, such as crustal-derived magmas (e.g., A- ad I-type granites and silicic volcanics). Within the crustal system pegmatites are known to have extremely high and variable Nb, Ta contents, often with low Nb/Ta. A fluid is generally considered to be involved in the generation of these rocks. Thus it is possible that fluid/melt partitioning may be the key to fractionating Nb and Ta, with preference for Ta in the fluid, and enrichment of Ta relative to Nb into the mid-upper crustal system, as the crust evolved, through upward movement of fluid.     

Intermediate high-pressure exhumation of the northern segment of the Sanbagawa belt, Saruta-gawa area, central Shikoku, Japan

The retrograde chemical zonal structure of amphibole in hematite-bearing basic and quartz schists from the higher grade zone in the Saruta-gawa area of the Sanbagawa belt was studied to investigate the relationships between the prograde and retrograde P–T paths of the Sanbagawa metamorphism. This amphibole coexists with chlorite, epidote, muscovite, albite, quartz and hematite, and is composed of Al-rich core and Al-poor mantle. The core is fairly homogeneous and has a barroisitic composition. In the mantle part, ^[B]Na increases with decreasing ^[4]Al towards the margins, which have winchite–magnesioriebeckite compositions. The barroisite–winchite–magnesioriebeckite composite crystal is sometimes rimmed by actinolite and/or winchite with low ^[4]Al and ^[B]Na. The Al-rich core and Al-poor mantle are regarded as prograde and retrograde products, respectively. The retrograde mantle in the Saruta-gawa area: (1) is systematically richer in ^[B]Na [0.40–1.73 per formula unit (pfu; for O=23)] than that from the same grade zone in the Asemi-gawa area (0.19–0.78 pfu), about 8 km south of the studied area; (2) tends to be ^[B]Na-poorer (less than 1.73 pfu) than prograde sodic amphibole (up to 1.93 ^[B]Na pfu) produced in the peak temperature stage from the lower grade zone in the same and other areas; and (3) extends its compositional range towards higher ^[B]Na and lower ^[4]Al than prograde-formed amphibole from the same grade zone in the same area. These zonal characteristics imply that (1) the Saruta-gawa samples experienced retrograde metamorphism under higher P/T conditions than the Asemi-gawa samples, (2) the retrograde P–T path of the Saruta-gawa area passes on the lower pressure side of the metamorphic field gradient, and (3) the Saruta-gawa samples underwent retrograde metamorphism under higher P/T conditions than the prograde metamorphism. The higher P/T conditions of the retrograde metamorphism suggests an increasing dP/dT of the geotherm during exhumation. Retrograde P–T conditions during the formation of magnesioriebeckite can be roughly estimated at 7–8 kbar, 400–450°C based on semi-quantitative phase relations of actinolite–winchite–magnesioriebeckite–barroisite series associated with chlorite, epidote, muscovite, albite, quartz and hematite.     

Paleomagnetism of Cretaceous units of the Mametchinskiy Peninsula, Kuyul Region, Northeastern Russia: implications for development and evolution of the Northwest Pacific Basin

The northernmost part of the Kamchatka Peninsula of northeastern Russia, located along the northwestern margin of the Bering Sea, consists of zones of complexly deformed accreted terranes. Progressing from the northwestern Bering Sea inland are the Olyutorskiy, Ukelayat, and Koryak superterranes, which were accreted to the Okhotsk–Chukotsk volcanic–plutonic belt (OChVB) during the Campanian–Maastrichtian (Koryak) to Middle Eocene (Olyutorskiy), respectively. To constrain the accretion paleolatitude of the Koryak superterrane, we paleomagnetically sampled a sedimentary series on the Mametchinskiy Peninsula. At the Mametchinskiy Peninsula, in the northeastern Penzhinskaya Guba (61.45° N, 163.75° E), a gently deformed, well-bedded section of fine-grained Lower to lower Upper Cretaceous turbidites, the Mametchinskaya and Tylakrylskaya Formations are exposed. These strata, which represent the lower part of the sedimentary cover of the terranes in this region and the forearc of OChVB, were sampled at 39 sites (three to seven samples per site). Within the Ainyn terrane, more than 1000 m of section of Cenomanian–Turonian age was sampled at a basal locality (sample groups I and II, sites 1–18, 19–29) and at an upper locality of Valanginian–Barremian age (sample group III, sites 30–39) along the western shore of the Peninsula. Thermal demagnetization and principal component analysis of the demagnetization data show lower-temperature (A) and higher-temperature (B) magnetic components. Although group III samples did not display a coherent A component, the A component of group I and II samples was observed as a single-polarity lower-unblocking temperature component generally removed by 100–400 °C. This component failed the fold test at the 95% confidence level. With respect to direction, the A component is similar to both the present-day field and axial–geocentric dipole directions expected at this site. The B component was observed during thermal demagnetization steps up to 580 °C and was always of downward-directed inclination. Coherence of bedding corrections within each section do not allow statistically meaningful fold tests within groups I, II or III. Assuming the B component represents a Cretaceous magnetization, two overall models are proposed. In the first model (preferred), with the highest clustering of directions (k-value=36.7, N (sites)=36), indicates significant poleward motion of the Ainyn terrane (observed paleolatitude λ_M1=61.0±6.5°; expected North America plate reference site paleolatitude λ_E=74.0±3.5°). In the second model, no significant poleward displacement is implied (λ_M2=72.0±9.6).