Assimilation and Fractional Crystallization Controlled by Transport Process of Crustal Melt: Implications from an Alkali Basalt-Dacite Suite from Rishiri Volcano, Japan

Mechanisms of fractional crystallization with simultaneous crustalassimilation (AFC) are examined for the Kutsugata and Tanetomilavas, an alkali basalt–dacite suite erupted sequentiallyfrom Rishiri Volcano, northern Japan. The major element variationswithin the suite can be explained by boundary layer fractionation;that is, mixing of a magma in the main part of the magma bodywith a fractionated interstitial melt transported from the mushyboundary layer at the floor. Systematic variations in SiO₂ correlatewith variations in the Pb, Sr and Nd isotopic compositions ofthe lavas. The geochemical variations of the lavas are explainedby a constant and relatively low ratio of assimilated mass tocrystallized mass (‘r value’). In the magma chamberin which the Kutsugata and Tanetomi magmas evolved, a strongthermal gradient was present and it is suggested that the marginalpart of the reservoir was completely solidified. The assimilantwas transported by crack flow from the partially fused floorcrust to the partially crystallized floor mush zone throughfractures in the solidified margin, formed mainly by thermalstresses resulting from cooling of the solidified margin andheating of the crust. The crustal melt was then mixed with thefractionated interstitial melt in the mushy zone, and the mixedmelt was further transported by compositional convection tothe main magma, causing its geochemical evolution to be characteristicof AFC. The volume flux of the assimilant from the crust tothe magma chamber is suggested to have decreased progressivelywith time (proportional to t^–1/2), and was about 3 x 10^–2m/year at t = 10 years and 1 x 10^–2 m/year at t = 100years. It has been commonly considered that the heat balancebetween magmas and the surrounding crust controls the couplingof assimilation and fractional crystallization processes (i.e.absolute value of r). However, it is inferred from this studythat the ratio of assimilated mass to crystallized mass canbe controlled by the transport process of the assimilant fromthe crust to magma chambers. KEY WORDS: assimilation and fractional crystallization; mass balance model; magma chamber; melt transport; Pb isotope     

Continuous Gradations among Primary Carbonatitic, Kimberlitic, Melilititic, Basaltic, Picritic, and Komatiitic Melts in Equilibrium with Garnet Lherzolite at 3-8 GPa

Multianvil melting experiments in the system CaO–MgO–Al₂O₃–SiO₂–CO₂(CMAS–CO₂) at 3–8 GPa, 1340–1800°C, involvingthe garnet lherzolite phase assemblage in equilibrium with CO₂-bearingmelts, yield continuous gradations in melt composition betweencarbonatite, kimberlite, melilitite, komatiite, picrite, andbasalt melts. The phase relations encompass a divariant surfacein P–T space. Comparison of the carbonatitic melts producedat the low-temperature side of this surface with naturally occurringcarbonatites indicates that natural magnesiocarbonatites couldbe generated over a wide range of pressures >2·5 GPa.Melts analogous to kimberlites form at higher temperatures alongthe divariant surface, which suggests that kimberlite genesisrequires more elevated geotherms. However, the amount of waterfound in some kimberlites has the potential to lower temperaturesfor the generation of kimberlitic melts by up to 150°C,provided no hydrous phases are present. Compositions resemblinggroup IB and IA kimberlites are produced at pressures around5–6 GPa and 10 GPa, respectively, whereas the compositionsof some other kimberlites suggest generation at higher pressuresstill. At pressures <4 GPa, an elevated geotherm producesmelilitite-like melt in the CMAS–CO₂ system rather thankimberlite. Even when a relatively CO₂-rich mantle compositioncontaining 0·15 wt % CO₂ is assumed, kimberlites andmelilitites are produced by <1% melting and carbonatitesare generated by even smaller degrees of melting of <0·5%. KEY WORDS: carbonatite; CO₂; kimberlite; melilitite; melt generation     

Origin of Exceptionally Abundant Phonolites on Ua Pou Island (Marquesas, French Polynesia): Partial Melting of Basanites Followed by Crustal Contamination

On the basis of the first systematic mapping of Ua Pou, longknown for its exceptionally abundant phonolites, we estimatethat these rocks cover 65% of the surface of the island whereasmafic lavas cover 27% and intermediate ones 8%. The silica-undersaturatedsuite was erupted in a restricted time span (2·9–2·35Myr), following the emplacement of tholeiites derived from ayoung HIMU-type source at c. 4 Ma. Primitive basanites, derivedfrom a heterogeneous mantle source with a dominant EM II + HIMUsignature, represent likely parental magmas. The series is characterizedby a Daly gap defined by a lack of phonotephrites. We considerthat the most likely model for the origin of evolved lavas ispartial melting at depth of primitive basanites, leaving anamphibole-rich residuum and producing tephriphonolitic magmas.These tephriphonolitic magmas may have evolved by closed-systemfractional crystallization towards Group A phonolites. Threeother groups of phonolites could have been derived from tephriphonoliticmagmas by open-system fractional crystallization processes,characterized respectively by seawater contamination (GroupB), assimilation of nepheline syenite-type materials (GroupC) and extreme fractionation coupled with assimilation of theunderlying oceanic crust (Group D). The prominence of evolvedlavas is a consequence of their origin from partial meltingof mafic precursors followed by crustal contamination. KEY WORDS: Marquesas; French Polynesia; phonolite; partial melting; contamination     

Experimental determination of argon solubility in silicate melts: An assessment of the effects of liquid composition and temperature

The argon solubility of 38 liquids in the system Na₂O-CaO-MgO-Al₂O₃-SiO₂ (NCMAS) has been determined at 1873 K and 1 bar, the argon concentration of presaturated glasses being measured using a static mass spectrometer. For compositions in the subsystem diopside (CaMgSi₂O₆), nepheline (NaAlSiO₄), albite (NaAlSi₃O₈), anorthite (CaAl₂Si₂O₈), argon solubility is generally a linear function of the relative proportion of each end member, solubility being lowest in diopside melt (1.53 10⁻⁵ cm³ STP · g⁻¹ · bar⁻¹) and highest in albite melt (2.88 10⁻⁴ cm³ STP · g⁻¹ · bar⁻¹). For the tectosilicate joins studied (SiO₂-Na₂Al₂O₄, SiO₂-CaAl₂O₄, SiO₂-MgAl₂O₄) solubility decreases with decreasing silica content in all cases, being highest for Na-bearing liquids and lowest for Mg-bearing liquids at constant molar silica content. Where comparison is possible our results are in good agreement with data from the literature. When our data are considered in isolation we find that argon solubility shows an excellent correlation with calculated ionic porosity. The covariation of argon solubility and liquid density is also reasonable, that with molar volume less convincing and that with polymerization state (as defined by the ratio of the number of nonbridging oxygens and tetrahedral network forming cations; NBO/T) nonexistent. However, when our data are combined with those from the literature no well constrained correlation between argon solubility and ionic porosity is apparent. Based upon this observation and consideration of the temperature dependence of noble gas solubility it is concluded that ionic porosity is not a universally applicable parameter which may be used to predict noble gas solubility as a function of composition, temperature and pressure. Two new models for calculating argon solubility are proposed, both employing the notion of partial molar argon solubilities. The first uses oxide components, for which partial molar argon solubility is directly proportional to partial molar ionic porosity calculated at 1873 K, irrespective of the temperature of experimental equilibration. The second model, which offers the best fit to the available data, employs tetrahedral units rather than oxides as the proposed melt components. This latter model successfully accounts for reported argon solubilities in simple Al-free systems, in simple Al-bearing systems and in natural liquids. This is interpreted to infer that argon is incorporated in large sites in the liquid structure (such as the space within rings of n-tetrahedra) although further work is required to understand the quantitative links between melt structure and noble gas solubility.     

Fluid generation and evolution during exhumation of deeply subducted UHP continental crust: Petrogenesis of composite granite–quartz veins in the Sulu belt,China

Composite granite–quartz veins occur in retrogressed ultrahigh pressure (UHP) eclogite enclosed in gneiss at General's Hill in the central Sulu belt, eastern China. The granite in the veins has a high‐pressure (HP) mineral assemblage of dominantly quartz+phengite+allanite/epidote+garnet that yields pressures of 2.5–2.1 GPa (Si‐in‐phengite barometry) and temperatures of 850–780°C (Ti‐in‐zircon thermometry) at 2.5 GPa (~20°C lower at 2.1 GPa). Zircon overgrowths on inherited cores and new grains of zircon from both components of the composite veins crystallized at c. 221 Ma. This age overlaps the timing of HP retrograde recrystallization dated at 225–215 Ma from multiple localities in the Sulu belt, consistent with the HP conditions retrieved from the granite. The ε_Hf(t) values of new zircon from both components of the composite veins and the Sr–Nd isotope compositions of the granite consistently lie between values for gneiss and eclogite, whereas δ¹⁸O values of new zircon are similar in the veins and the crustal rocks. These data are consistent with zircon growth from a blended fluid generated internally within the gneiss and the eclogite, without any ingress of fluid from an external source. However, at the peak metamorphic pressure, which could have reached 7 GPa, the rocks were likely fluid absent. During initial exhumation under UHP conditions, exsolution of H₂O from nominally anhydrous minerals generated a grain boundary supercritical fluid in both gneiss and eclogite. As exhumation progressed, the volume of fluid increased allowing it to migrate by diffusing porous flow from grain boundaries into channels and drain from the dominant gneiss through the subordinate eclogite. This produced a blended fluid intermediate in its isotope composition between the two end‐members, as recorded by the composite veins. During exhumation from UHP (coesite) eclogite to HP (quartz) eclogite facies conditions, the supercritical fluid evolved by dissolution of the silicate mineral matrix, becoming increasingly solute‐rich, more ‘granitic’ and more viscous until it became trapped. As crystallization began by diffusive loss of H₂O to the host eclogite concomitant with ongoing exhumation of the crust, the trapped supercritical fluid intersected the solvus for the granite–H₂O system, allowing phase separation and formation of the composite granite–quartz veins. Subsequently, during the transition from HP eclogite to amphibolite facies conditions, minor phengite breakdown melting is recorded in both the granite and the gneiss by K‐feldspar+plagioclase+biotite aggregates located around phengite and by K‐feldspar veinlets along grain boundaries. Phase equilibria modelling of the granite indicates that this late‐stage melting records P–T conditions towards the end of the exhumation, with the subsolidus assemblage yielding 0.7–1.1 GPa at <670°C. Thus, the composite granite–quartz veins represent a rare example of a natural system recording how the fluid phase evolved during exhumation of continental crust. The successive availability of different fluid phases attending retrograde metamorphism from UHP eclogite to amphibolite facies conditions will affect the transport of trace elements through the continental crust and the role of these fluids as metasomatic agents interacting with the mantle wedge in the subduction channel.     

Micromechanical approach to swelling behavior of capillary-porous media with coupled physics

A detailed multiscale analysis is presented of the swelling phenomenon in unsaturated clay-rich materials in the linear regime through homogenization. Herein, the structural complexity of the material is formulated as a three-scale, triple porosity medium within which microstructural information is transmitted across the various scales, leading ultimately to an enriched stress-deformation relation at the macroscopic scale. As a side note, such derived relationship leads to a tensorial stress partitioning that is reminiscent of a Terzaghi-like effective stress measure. Otherwise, a major result that stands out from previous works is the explicit expression of swelling stress and capillary stress in terms of micromechanical interactions at the very fine scale down to the clay platelet level, along with capillary stress emerging due to interactions between fluid phases at the different scales, including surface tension, pore size, and morphology. More importantly, the swelling stress is correlated with the disjoining forces due to electrochemical effects of charged ions on clay minerals and van der Waals forces at the nanoscale. The resulting analytical expressions also elucidate the role of the various physics in the deformational behavior of clayey material. Finally, the capability of the proposed formulation in capturing salient behaviors of unsaturated expansive clays is illustrated through some numerical examples.     

Constitutive modelling of fine-grained gassy soil: A composite approach

Fine-grained marine sediments containing large undissolved gas bubbles are widely distributed around the world. Presence of the bubbles could degrade the undrained shear strength (s_u ) of the soil, when the gas pressure u_g is relatively high as compared with the effective stress in the saturated soil matrix. Meanwhile, the addition of bubbles may also increase s_u when the difference between u_g and pore water pressure u_w becomes smaller than the water entry value, causing partial water drainage from the saturated matrix into the bubbles (bubble flooding) during globally undrained shearing. A new constitutive model for describing the two competing effects on the stress-strain relationship of fine-grained gassy soil is proposed within the framework of critical state soil mechanics. The gassy soil is considered as a three-phase composite material with compressible cavities, which allows water entry from the saturated matrix. Bubble flooding is modelled by introducing an additional positive volumetric strain increment of the saturated clay matrix, which is dependent on the difference between pore gas and pore water pressure based on experimental observations. A modified hardening law based on that of the modified Cam clay model is employed, which in conjunction with the expression for bubble flooding, can describe both the detrimental and beneficial effects of gas bubbles on soil strength and plastic hardening in shear. Only two extra parameters in addition to those in the modified Cam clay model are used. It is shown that the key features of the stress-strain relationship of three fine-grained gassy soils can be reproduced satisfactorily.     

青藏高原拉萨地块西部念青唐古拉岩群的地球化学特征及构造意义

拉萨地块西部呈断块状沿狮泉河-申扎-嘉黎蛇绿混杂岩带附近分布的念青唐古拉岩群被认为是前寒武纪变质基底。本文对念青唐古拉岩群进行了系统的岩石学、地球化学、同位素年代学及构造地质学研究。研究结果表明片岩-片麻岩-变粒岩含十字石、石榴子石等特征变质矿物,遵循粒度分异规律,其原岩可能为来自冈瓦纳古陆核北缘中新元古代弧盆体系的活动大陆边缘浊积岩。斜长角闪岩具低硅、高铁镁、富钙的基性岩特征,其原岩为岛弧型基性火山岩。念青唐古拉岩群中的花岗伟晶岩锆石LA-ICPMS U-Pb年龄为1150±13Ma,具过铝质S型花岗岩地球化学特征,可能为中元古代(1150±13 Ma)以前就开始沉积的念青唐古拉岩群基底岩石通过部分熔融形成。与花岗伟晶岩渐变过渡接触的二云斜长片麻岩第一组变质重结晶锆石U-Pb年龄为701±15 Ma,结合十字石特征变质矿物,暗示了该地区中温高压变质作用的峰期变质,变质程度达角闪岩相;第二组热液流体锆石UPb年龄为301±8.4 Ma,可能与冈瓦纳大陆北缘古特提斯洋演化过程中的岩浆热液作用有关。     

部分变量误差模型的整体抗差最小二乘估计

部分变量误差模型(partial EIV model)的加权整体最小二乘(weighted total least-squares,WTLS)估计不具备抵御粗差的能力。鉴于粗差可能同时出现在观测值和系数矩阵中,本文在提出部分变量误差模型WTLS估计的两步迭代解法的基础上,运用抗差M估计的等价权方法,发展了一种整体抗差最小二乘(TRLS)估计方法,并采用一致最大功效统计量确定降权因子。针对WTLS估计两步迭代解法的特点,设计了两个不同的降权方案:第1个方案是在估计系数矩阵元素时,不对观测值降权,仅对系数矩阵降权;第2个方案是在估计系数矩阵元素时,既对系数矩阵降权,同时也对观测值降权。通过对模拟2D仿射变换和线性拟合实例进行计算和分析,结果表明第1方案优于第2方案,并且优于基于残差和验后单位权方差的抗差估计和现有的变量误差模型抗差估计。     

BDS／GPS 组合 RTK 定位性能分析

针对单系统RTK存在可见卫星数少等问题,文中研究BDS/GPS站间单差的RTK算法模型,该模型采用二次型函数部分最小化及LAMBDA方法联合搜索模糊度。利用该模型分析BDS/GPS组合RTK的定位性能,通过短基线实测数据分析表明:站间单差RTK模型与双差模型是等价的;BDS/GPS组合系统相比于单一系统,明显提高定位的稳健性和精度,改善模糊度固定的成功率。