Sector collapse,sedimentation and clast population evolution at an active island-arc volcano: Stromboli,Italy

Sciara del Fuoco is the subaerial part of a partially filled sector-collapse scar that extends to 700 m below sea level on Stromboli volcano. The collapse occurred <5000 years ago, involved 1.81 km³ of rock and is the latest of a series of major collapses on the north-west flank of Stromboli. A north-east trending arc-axial fault system channels magmas into the volcano and has caused tilting and/or downthrow to the north-west. The slope of the partial cone constructed between the lateral walls of the collapse scar acts as a channelway to the sea for most eruptive products. From 700 m below sea level and extending to >2200 m and >10 km from the shore to the NNW, a fan-shaped mounded feature comprises debris avalanche deposits (>4 km³) from two or more sector collapses. Volcaniclastic density currents originating from Sciara del Fuoco follow the topographic margin of the debris avalanche deposits, although overbank currents and other unconfined currents widely cover the mounded feature with turbidites. Historical (recorded) eruptive activity in Sciara del Fuoco is considerably less than that which occurred earlier, and much of the partial fill may have formed from eruptions soon after the sector collapse. It is possible that a mass of eruptive products similar to that in the collapse scar is dispersed as volcanogenic sediment in deep water of the Tyrhennian basin. Evidence that the early post-collapse eruptive discharge was greater than the apparent recent flux (2kg/s) counters suggestions that a substantial part of Stromboli's growth has been endogenous. The partial fill of Sciaria del Fuoco is dominated by lava and spatter layers, rather than by the scoria and ash layers classically regarded as main constituents of Strombolian (cinder) cones. Much of the volcanic slope beneath the vents is steeper than the angle of repose of loose tephra, which is therefore rapidly transported to the sea. Delicate pyroclasts that record the magmatic explosivity are selectively destroyed and diluted during sedimentary transport, mainly in avalanches and by shoreline wave reworking, and thus the submarine deposits do not record well the extent and diversity of explosive activity and associated clast-forming processes. Considerable amounts of dense (non-vesicular) fine sand and silt grains are produced by breakage and rounding of fragments of lava and agglutinate. The submarine extension of the collapse scar, and the continuing topographic depression to >2200 m below sea level, are zones of considerable by-passing of fine sand and silt, which are transported in turbidity currents. Evidently, volcanogenic sediments dispersed around island volcanoes by density currents are unlikely to record well the true spectrum and relative importance of clast-forming processes that occurred during an eruption. Marine sedimentary evidence of magmatic explosivity is particularly susceptible to partial or complete obliteration, unless there is a high rate of discharge of pyroclastic material into the sea.     

午后多重极光弧观测研究

本文利用2001-2003年南极中山站175天全天空摄像机观测,对午后多重极光弧的出现率及其与Kp指数的关系进行了统计分析,结果表明午后多重极光弧出现率呈一单峰分布,最大发生率出现在1445UT(1645MLT),其位置在1500MLT极光热点(1300-1700MLT)近夜侧的部分。与地磁活动指数Kp的相关统计分析表明,Kp值为2-3之间时多重极光弧有较大的出现率,这说明中等地磁活动情形下午后多重极光弧有较高的出现率。事件分析表明多重极光弧的强度变化与地磁Pc5脉动具有较高的相关性,并且有类似的频谱特征,这说明午后多重极光弧可能与同时出现的Pc5地磁脉动有关。     

Early Cretaceous volcanic rocks and Early Cenozoic extrusions of Cape Mary, Schmidt Peninsula, north Sakhalin: Geochemical study

The Early Cretaceous volcanic rocks of the Mariisky sequence and Early Cenozoic extrusive-vent rocks of Cape Mary are exposed at the northwestern extremity of Schmidt Peninsula, north Sakhalin. In chemical composition, all the rocks are subdivided into four groups. Three groups include the volcanic rocks of the Mariisky sequence, which consists, from bottom to top, of calc-alkaline rocks, transitional calc-alkaline-tholeiite rocks, and incompatible element-depleted tholeiites. These rocks show subduction geochemical signatures and are considered as a fragment of the Moneron-Samarga island arc system. Trace-element modeling indicates their derivation through successive melting of garnet-bearing mantle and garnet-free shallower mantle sources containing amphibole; pyroxene; and, possibly, spinel. The mixed subduction and within-plate characteristics of the extrusive vent rocks of Cape Mary attest to their formation in a transform continental margin setting.     

Geochemical Evidence for Slab Melting in the Trans-Mexican Volcanic Belt

Geochemical studies of Plio-Quaternary volcanic rocks from theValle de Bravo–Zitácuaro volcanic field (VBZ) incentral Mexico indicate that slab melting plays a key role inthe petrogenesis of the Trans-Mexican Volcanic Belt. Rocks fromthe VBZ are typical arc-related high-Mg andesites, but two differentrock suites with distinct trace element patterns and isotopiccompositions erupted concurrently in the area, with a traceelement character that is also distinct from that of other Mexicanvolcanoes. The geochemical differences between the VBZ suitescannot be explained by simple crystal fractionation and/or crustalassimilation of a common primitive magma, but can be reconciledby the participation of different proportions of melts derivedfrom the subducted basalt and sediments interacting with themantle wedge. Sr/Y and Sr/Pb ratios of the VBZ rocks correlateinversely with Pb and Sr isotopic compositions, indicating thatthe Sr and Pb budgets are strongly controlled by melt additionsfrom the subducted slab. In contrast, an inverse correlationbetween Pb(Th)/Nd and ¹⁴³Nd/¹⁴⁴Nd ratios, which extend to lowerisotopic values than those for Pacific mid-ocean ridge basalts,indicates the participation of an enriched mantle wedge thatis similar to the source of Mexican intraplate basalts. In addition,a systematic decrease in middle and heavy rare earth concentrationsand Nb/Ta ratios with increasing SiO₂ contents in the VBZ rocksis best explained if these elements are mobilized to some extentin the subduction flux, and suggests that slab partial fusionoccurred under garnet amphibolite-facies conditions. KEY WORDS: arcs; mantle; Mexico; sediment melting; slab melting     

Lead angles and emitting electron energies of Io-controlled decameter radio arcs

The Io-controlled radio arcs are emissions in the decametric radio range which appear arc shaped in the time-frequency plane. Their occurrence is controlled by Io's position, so it has been for long inferred that they are powered by the Io-Jupiter electrodynamic interaction. Their frequency ranges correspond to the electron cyclotron frequencies along the Io Flux tube, so they are expected to be generated by cyclotron maser instability (CMI). The arc shape was proposed to be a consequence of the strong anisotropy of the decametric radio emissions beaming, combined with the topology of the magnetic field in the source and the observation geometry. Recent papers succeeded at reproducing the morphologies of a few typical radio arcs by modeling in three dimensions the observation geometry, using the best available magnetic field model and a beaming angle variation consistent with a loss-cone driven CMI. In the continuation of these studies, we present here the systematic modeling of a larger number of observations of the radio arcs emitted in Jupiter's southern hemisphere (including multiple arcs or arcs exhibiting abrupt changes of shape), which permits to obtain a statistical determination of the emitting field line localization (lead angle) relative to the instantaneous Io field line, and of the emitting particle velocities or energies. Variations of these parameters relative to Io's longitude are also measured and compared to the location of the UV footprints of the Io-Jupiter interaction. It is shown that the data are better organized in a reference frame attached to the UV spot resulting from the main Alfvén wing resulting from the Io-Jupiter interaction. It is proposed that the radio arcs are related to the first reflected Alfvén wing rather than to the main one.     

The general solution of the HENON–HEILES problem

The general solution of the Henon–Heiles system is approximated inside a domain of the (x, C) of initial conditions (C is the energy constant). The method applied is that described by Poincaré as ‘the only “crack” permitting penetration into the non-integrable problems’ and involves calculation of a dense set of families of periodic solutions that covers the solution space of the problem. In the case of the Henon–Heiles potential we calculated the families of periodic solutions that re-enter after 1–108 oscillations. The density of the set of such families is defined by a pre-assigned parameter ε (Poincaré parameter), which ascertains that at least one periodic solution is computed and available within a distance ε from any point of the domain (x, C) for which the approximate general solution computed. The approximate general solution presented here corresponds to ε = 0.07. The same solution is further improved by “zooming” into four square sub-domain of (x, C), i.e. by computing sufficient number of families that reduce the density parameter to ε = 0.003. Further zooming to reduce the density parameter, say to ε = 10⁻⁶, or even smaller, although easily performable in both areas occupied by stable as well as unstable solutions, was found unnecessary. The stability of all members of each and all families computed was calculated and presented in this paper for both the large solution domain and for the sub-domains. The correspondence between areas of the approximate general solution occupied by stable periodic solutions and Poincaré sections with well-aligned section points and also correspondence between areas occupied by unstable solutions and Poincaré sections with randomly scattered section points is shown by calculating such sections. All calculations were performed using the Runge-Kutta (R-K) 8th order direct integration method and the large output received, consisting of many thousands of families is saved as “Atlas of the General Solution of the Henon–Heiles Problem,” including their stability and is available at request. It is concluded that approximation of the general solution of this system is straightforward and that the chaotic character of its Poincaré sections imposes no limitations or difficulties.     

基于粒子滤波的行车轨迹路网匹配方法

行车轨迹是一种时间序列的地理空间位置采样数据,而传统的轨迹—路网匹配方法主要以全局或局部寻优的方式建立轨迹—路网匹配关系,影响了时空场景中数据的匹配计算过程的相对独立性。针对这个问题,本文基于粒子滤波（Particle Filter,PF）原理建立行车轨迹与道路网络之间的匹配关系。首先,沿轨迹中车辆运动方向在道路网络中搜索邻近道路节点,在与道路节点拓扑邻接的道路弧段上初始化随机生成粒子,根据轨迹中车辆运动模型将粒子沿所在道路弧段移动;然后,基于PF原理计算各时刻粒子运动状态及与行车轨迹采样点之间的距离误差,根据高斯概率密度函数计算粒子权重并利用随机重采样方法进行粒子重采样,迭代更新粒子运动状态;最后,计算与搜索到的道路节点拓扑邻接的每条道路弧段中累计粒子权重,通过各道路弧段累计权重计算轨迹—路网匹配关系。以行车轨迹进行实验表明,利用本文方法可以通过粒子时空变化反映采样点的移动,行车轨迹—路网匹配结果的正确率大于85%,能够实现行车轨迹和路网的准确匹配。     

中国东北地区显生宙岩浆作用和洋-陆格局及其与气候演变的关系SCIEI北大核心CSCD

显生宙期间,地球经历了温室-冰室气候的周期性交替变化。在数百万年的时间尺度,这种古气候的转变被认为是碳源和碳汇过程耦合的结果,但一直以来关于两者贡献程度的认识尚不明确。通过全球统计分析,不同学者提出大陆弧火山脱气模型和热带弧-陆碰撞模型用于解释整个显生宙古气候的演变,分别强调了碳源和碳汇的一级控制作用。为了检验上述模型,更好地理解古气候的转变机制和演化细节,本文系统总结了中国东北地区显生宙岩浆作用-矽卡岩型矿床的时空展布和构造背景,以及弧-陆碰撞的时代、规模和古地理位置,通过数据统计和作图对比,发现东北地区岩浆作用-矽卡岩成矿峰期、弧-陆碰撞缝合带的时空迁移与大气圈CO_(2)浓度和大陆冰川沉积有很好的对应关系,暗示东北显生宙构造-岩浆过程和古气候演变的内在联系。综合东北地区及全球的研究进展,本文提出如下倾向性认识:1)洋-陆俯冲过程中火山-变质脱气的强度决定了CO_(2)排放量,而热带区域弧-陆碰撞缝合带的规模决定了全球硅酸盐风化速率和CO_(2)吸收量,在地质演化过程中两者紧密联动,共同控制了显生宙古气候的演变;2)大陆弧岩浆作用的全球爆发不一定能造成温室气候的出现,如果缺乏充分矽卡岩变质脱碳反应,大陆弧CO_(2)排放通量与岛弧、大洋中脊和板内并无显著区别;3)SO_(2)属于短期效应气体,理论和实例研究均暗示爆发式火山作用难以诱发大冰期的形成,火山作用之于长期气候应该仍是促使地球升温而非变冷。     

Horizontal layered structure with heterogeneity beneath continents and island arcs from particle orbits of long-period P waves

We investigate the particle orbits of long-period (about 20 s) P waves observed with the global seismic network. By analysing 84 three-component seismograms recorded at 25 stations from 60 earthquakes occurring beneath 300 km, we quantitatively evaluate the orbits by three sets of eigenvalues and eigenvectors, using a covariance matrix method. The eigenvalues for P waves recorded at stations located on continents are explained by the standard horizontal layered structure model (iasp91). On the other hand, the orbits observed at stations close to island arcs are affected not only by the horizontal layered structure but also by heterogeneity due to subducting plates, mantle diapirs and so on. On the basis of a single-scattering model for a plane P wave, we quantify the heterogeneities by an isotropic scattering coefficient g0. Fitting the theoretical eigenvalues to the observed ones, we estimate g0 for the crust and upper mantle beneath continents to be less than 0.0005 km^-1, and the mean g₀ for the structure beneath island arcs to be about 0.0015 to 0.003 km^-1.