Low‐P and high‐T metamorphism of basalts: Insights from the Sudbury impact melt sheet aureole and thermodynamic modelling

Low‐pressure and high‐temperature (LP–HT) metamorphism of basaltic rocks, which occurs globally and throughout geological time, is rarely constrained by forward phase equilibrium modelling, yet such calculations provide valuable supplementary thermometric information and constraints on anatexis that are not possible to obtain from conventional thermometry. Metabasalts along the southern margin of the Sudbury Igneous Complex (SIC) record evidence of high‐grade contact metamorphism involving partial melting and melt segregation. Peak metamorphic temperatures reached at least ~925°C at ~1–3 kbar near the SIC contact. Preservation of the peak mineral assemblage indicates that most of the generated melt escaped from these rocks leaving a residuum characterized by a plagioclase–orthopyroxene–clinopyroxene–ilmenite‐magnetite±melt assemblage. Peak temperatures reached ~875°C up to 500 m from the SIC lower contact, which marks the transition to metabasalts that only experienced incipient partial melting without melt loss. Metabasalts ~500 to 750 m from the SIC contact are characterized by a similar two‐pyroxene mineral assemblage, but typically contain abundant hornblende that overgrew clino‐ and orthopyroxene along an isobaric cooling path. Metabasalts ~750 to 1,000 m from the SIC contact are characterized by a hornblende–plagioclase–quartz–ilmenite assemblage indicating temperatures up to ~680°C. Mass balance and phase equilibria calculations indicate that anatexis resulted in 10–20% melt generation in the inner ~500 m of the aureole, with even higher degrees of melting towards the contact. Comparison of multiple models, experiments, and natural samples indicates that modelling in the Na₂O–CaO–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (NCFMASHTO) system results in the most reliable predictions for the temperature of the solidus. Incorporation of K₂O in the most recent amphibole solution model now successfully predicts dehydration melting by the coexistence of high‐Ca amphibole and silicate melt at relatively low pressures (~1.5 kbar). However, inclusion of K₂O as a system component results in prediction of the solidus at too low a temperature. Although there are discrepancies between modelling predictions and experimental results, this study demonstrates that the pseudosection approach to mafic rocks is an invaluable tool to constrain metamorphic processes at LP–HT conditions.     

Constraints on the timing and conditions of high‐grade metamorphism,charnockite formation and fluid–rock interaction in the Trivandrum Block,southern India

Incipient charnockites have been widely used as evidence for the infiltration of CO₂‐rich fluids driving dehydration of the lower crust. Rocks exposed at Kakkod quarry in the Trivandrum Block of southern India allow for a thorough investigation of the metamorphic evolution by preserving not only orthopyroxene‐bearing charnockite patches in a host garnet–biotite felsic gneiss, but also layers of garnet–sillimanite metapelite gneiss. Thermodynamic phase equilibria modelling of all three bulk compositions indicates consistent peak‐metamorphic conditions of 830–925 °C and 6–9 kbar with retrograde evolution involving suprasolidus decompression at high temperature. These models suggest that orthopyroxene was most likely stabilized close to the metamorphic peak as a result of small compositional heterogeneities in the host garnet–biotite gneiss. There is insufficient evidence to determine whether the heterogeneities were inherited from the protolith or introduced during syn‐metamorphic fluid flow. U–Pb geochronology of monazite and zircon from all three rock types constrains the peak of metamorphism and orthopyroxene growth to have occurred between the onset of high‐grade metamorphism at c. 590 Ma and the onset of melt crystallization at c. 540 Ma. The majority of metamorphic zircon growth occurred during protracted melt crystallization between c. 540 and 510 Ma. Melt crystallization was followed by the influx of aqueous, alkali‐rich fluids likely derived from melts crystallizing at depth. This late fluid flow led to retrogression of orthopyroxene, the observed outcrop pattern and to the textural and isotopic modification of monazite grains at c. 525–490 Ma.     

基于时相变化的晴空条件下Himawari-8 AHI火点检测

Himawari-8静止气象卫星具有高空间分辨率、高观测频次和高时效特点,对于火点检测具有很强优势。对Himawari-8卫星的3.9μm和11.2μm两通道亮温进行了连续时相变化研究,得出两通道的亮温在时间上的变化差值稳定且规律明显。根据两通道的亮温时相特征,考虑白天可见光对3.9μm通道的影响,并结合火点产生时引起的亮温变化特征,提出了适用于晴空条件下改进的火点检测算法。在多处进行了此算法的实验,例如2018-11-27 T 16:40(UTC时)河北张家口市桥东区一化工厂附近发生的严重爆炸起火事件以及2019-02-28澳大利亚西南部发生的火灾事件,均快速有效的检测到了火点。实验表明,改进的火点检测算法能很好的进行火点检测,并能解决晨昏交界、冰雪下垫面、常规火源点、太阳耀光等火点检测的难题。     

P–T–t–d evolution of orogenic middle crust of the Roc de Frausa Massif (Eastern Pyrenees): a result of horizontal crustal flow and Carboniferous doming?

Structural, petrological and textural studies are combined with phase equilibria modelling of metapelites from different structural levels of the Roc de Frausa Massif in the Eastern Pyrenees. The pre‐Variscan lithological succession is divided into the Upper, Intermediate and Lower series by two orthogneiss sheets and intruded by Variscan igneous rocks. Structural analysis reveals two phases of Variscan deformation. D1 is marked by tight to isoclinal small‐scale folds and an associated flat‐lying foliation (S1) that affects the whole crustal section. D2 structures are characterized by tight upright folds facing to the NW with steep NE–SW axial planes. D2 heterogeneously reworks the D1 fabrics, leading to an almost complete transposition into a sub‐vertical foliation (S2) in the high‐grade metamorphic domain. All structures are affected by late open to tight, steeply inclined south‐verging NW–SE folds (F3) compatible with steep greenschist facies dextral shear zones of probable Alpine age. In the micaschists of the Upper series, andalusite and sillimanite grew during the formation of the S1 foliation indicating heating from 580 to 640 °C associated with an increase in pressure. Subsequent static growth of cordierite points to post‐D1 decompression. In the Intermediate series, a sillimanite–biotite–muscovite‐bearing assemblage that is parallel to the S1 fabric is statically overgrown by cordierite and K‐feldspar. This sequence points to ~1 kbar of post‐D1 decompression at 630–650 °C. The Intermediate series is intruded by a gabbro–diorite stock that has an aureole marked by widespread migmatization. In the aureole, the migmatitic S1 foliation is defined by the assemblage biotite–sillimanite–K‐feldspar–garnet. The microstructural relationships and garnet zoning are compatible with the D1 pressure peak at ~7.5 kbar and ~750 °C. Late‐ to post‐S2 cordierite growth implies that F2 folds and the associated S2 axial planar leucosomes developed during nearly isothermal decompression to <5 kbar. The Lower series migmatites form a composite S1–S2 fabric; the garnet‐bearing assemblage suggests peak P–T conditions of >5 kbar at suprasolidus conditions. Almost complete consumption of garnet and late cordierite growth points to post‐D2 equilibration at <4 kbar and <750 °C. The early metamorphic history associated with the S1 fabric is interpreted as a result of horizontal middle crustal flow associated with progressive heating and possible burial. The upright F2 folding and S2 foliation are associated with a pressure decrease coeval with the intrusion of mafic magma at mid‐crustal levels. The D2 tectono‐metamorphic evolution may be explained by a crustal‐scale doming associated with emplacement of mafic magmas into the core of the dome.     

Dynamic characteristics of marine soft clay under variable phase difference and initial static shear stress

Abstract

Under seismic loading, the soil layer is subjected to multidirectional cyclic shear stress with different amplitudes and frequencies because of the coupling of multiple shear waves and the soil element within a slope or behind a retaining wall is subjected to initial static shear stress before subjected to cyclic loading. Due to the complexity of seismic loading propagation, a phase difference exists between the initial static shear stress and cyclic shear stress. To investigate the influence of the phase difference and initial static shear stress on cyclic shear strain, cyclic modulus, and cyclic strength, a series of laboratory tests are performed on Wenzhou marine soft clay by multi-directional simple shear system, which can simulate the actual state better by controlling the horizontal cyclic stress in the x and y directions simultaneously. As the phase difference varies from 0° to 90°, the dynamic shear modulus increases and cyclic strain accumulation decreases with an increasing number of cycles. The shear strain increases with the initial shear stress.     

北太平洋副热带模态水盐度的年代际变化

副热带模态水(Subtropical Mode Water;STMW)在气候变化中起着重要作用。本文利用全球高分辨率数值模拟结果,研究了北太平洋STMW核心层盐度(Core Layer Salinity;CLS)的年代际变化及其物理机制。结果表明,CLS存在显著的年代际变化,其空间分布则与背景流场分布特征有关。侵蚀区CLS滞后生成区CLS约1~2年,这主要是海流平流输运引起的。生成区内,STMW的季节循环一般可分为生成期(12-4月)、隔离期(5-6月)和侵蚀期(7-11月),生成期混合层盐度(Mixed Layer Salinity;MLS)决定着隔离期和侵蚀期的CLS,而MLS年代际变化则主要由同太平洋年代际涛动存在负相关性的海表面淡水通量的变化引起。     

月球交会对接VLBI差分时延研究

中国探月3期任务中,月球交会对接技术是任务成功的重要保障.利用嫦娥3号(CE03)绕月飞行的VLBI (Very Long Baseline Interferometry)时延数据,模拟仿真绕月交会对接过程中,同波束VLBI观测模式下,差分群时延的变化情况.仿真结果显示,在远程导引段,轨道器和上升组合体轨道距离保持100 km,持续半小时,差分群时延很好地反映了两者的轨道信息,可以用于定轨定位;自主控制段,上升组合体靠近轨道器,在轨道距离从5 km减小到20 m过程中,上升组合体加速追赶轨道器时,差分群时延快速趋近于0,上升组合体减速远离轨道器时,差分群时延绝对值快速变大.最后,利用嫦娥3号奔月段同时发射两个DOR (Differential One-Ranging)信号的VLBI时延数据,计算差分相时延,初步展示了月球交会对接过程中同波束VLBI差分相时延的误差情况.     

2016与1998年春季北大西洋海表温度异常的差异及成因

利用再分析资料以及混合层海温诊断方程, 研究1997—1998与2015—2016年超级厄尔尼诺次年北大西洋海表温度异常(sea surface temperature anomalies, SSTA)的差异及成因。结果显示, 北大西洋SSTA在1998年春季呈明显正负正三极型式分布, 而在2016年呈弱的负正负型态。诊断热带北大西洋SSTA的影响因素表明, 1998年春季暖SSTA除了之前研究强调的海洋表面向大气的潜热输送异常减少, 以及吸收太阳辐射的增加外, 海洋动力过程即Ekman纬向漂流也起着重要的作用。热力过程与厄尔尼诺峰值后出现的北大西洋涛动(North Atlantic Oscillation, NAO)负位相有关, 其可引起亚速尔高压减弱, 产生西南风异常, 通过风-蒸发-海表温度(sea surface temperature, SST)反馈机制使热带北大西洋蒸发减弱, 海表增暖, 沃克环流下沉支的东移对这一增暖也有贡献。与1997—1998厄尔尼诺事件不同, 2015—2016厄尔尼诺事件没有强迫出负位相NAO, 而是出现弱NAO正位相, 热带北大西洋为弱的东风异常, 使海表发生一定的冷却, 形成2016春季北大西洋SSTA与1998年的明显差异。     

Phase equilibria modelling and LASS monazite petrochronology: P–T–t constraints on the evolution of the Priest River core complex,northern Idaho

Phase equilibria modelling, laser‐ablation split‐stream (LASS)‐ICP‐MS petrochronology and garnet trace‐element geochemistry are integrated to constrain the P–T–t history of the footwall of the Priest River metamorphic core complex, northern Idaho. Metapelitic, migmatitic gneisses of the Hauser Lake Gneiss contain the peak assemblage garnet + sillimanite + biotite ± muscovite + plagioclase + K‐feldspar ± rutile ± ilmenite + quartz. Interpreted P–T paths predict maximum pressures and peak metamorphic temperatures of ~9.6–10.3 kbar and ~785–790 °C. Monazite and xenotime ²⁰⁸Pb/²³²Th dates from porphyroblast inclusions indicate that metamorphism occurred at c. 74–54 Ma. Dates from HREE‐depleted monazite formed during prograde growth constrain peak metamorphism at c. 64 Ma near the centre of the complex, while dates from HREE‐enriched monazite constrain the timing of garnet breakdown during near‐isothermal decompression at c. 60–57 Ma. Near‐isothermal decompression to ~5.0–4.4 kbar was followed by cooling and further decompression. The youngest, HREE‐enriched monazite records leucosome crystallization at mid‐crustal levels c. 54–44 Ma. The northernmost sample records regional metamorphism during the emplacement of the Selkirk igneous complex (c. 94–81 Ma), Cretaceous–Tertiary metamorphism and limited Eocene exhumation. Similarities between the Priest River complex and other complexes of the northern North American Cordillera suggest shared regional metamorphic and exhumation histories; however, in contrast to complexes to the north, the Priest River contains less partial melt and no evidence for diapiric exhumation. Improved constraints on metamorphism, deformation, anatexis and exhumation provide greater insight into the initiation and evolution of metamorphic core complexes in the northern Cordillera, and in similar tectonic settings elsewhere.     

InSAR相位解缠算法的比较

相位解缠作为合成孔径雷达干涉测量(InSAR,Interferometric Synthetic Aperture Radar)中一个关键的步骤,为获取地形高程提供了基础。相位解缠的质量将直接影响着DEM生成的质量。对现有的几类算法进行实验分析,并对解缠结果的质量进行比较,分析其精确性。结果表明:统计耗费网络流算法具有较好的解缠连续性以及较高的精度,可以获得一个较优的全局解。