A virtual voyage through 3D structures generated by chaotic mixing of magmas and numerical simulations: a new approach for understanding spatial and temporal complexity of magma dynamics

In this contribution, we present a virtual voyage through 3D structures generated by chaotic mixing of magmas and numerical simulations with the aim to highlight the power of 3D representations in the understanding of this geological phenomenon. In particular, samples of mixed juveniles from Salina island (Southern Italy) are reconstructed in 3D by serial lapping and digital montage and numerical simulations are performed by using a 3D chaotic dynamical system. Natural and simulated magma mixing structures are visualized by using several multimedia tools including animations and “virtual reality” models. It is shown that magma interaction processes can generate large spatial and temporal compositional heterogeneities in magmatic systems. The same topological structures are observed in both 3D reconstructed rock samples and chaotic numerical simulations, indicating that the mixing of magmas is governed by chaotic dynamics. The use of 3D multimedia models gives the opportunity to penetrate into magma mixing structures and to understand their significance in the context of magma dynamics. Such an approach is very powerful since multimedia tools can strongly capture the attention of the reader bringing him/her into an interactive and memorable geological experience. Electronic supplementary material  enclosed:     

幔源CO_2释出机理、脱气模式及成藏机制研究进展

针对幔源CO2如何从地幔岩浆中脱出并进入沉积地层中形成CO2气藏聚集这一关键问题,总结了国内外研究进展和前缘方向。研究表明,地幔深部的碱性玄武岩浆和碱性岩浆才是深部流体和CO2等挥发份大量赋存、渗滤和释出的场所。浅成侵入岩、次火山岩和火山通道等是CO2释放和聚集的有利位置,岩浆期后和岩浆衰弱期的热液活动阶段是CO2大量释放和聚集的有利时期。幔源CO2进入沉积盆地中具有3种脱气模式,即沿岩石圈断裂直接脱气模式、热流底辟体脱气模式和壳内岩浆房-基底断裂组合脱气模式。CO2的固有物化性质决定其运移相态多样,具有运移和聚集过程同步的特征。只有在满足大量的化学消耗及地层水或原油的溶解和耗散之后才能形成CO2有效聚集。幔源CO2成藏和分布主要受岩浆气源体和气源断裂体系的控制。今后,在超临界CO2及其对油气运移聚集的作用、CO2与深大断裂及火山岩的关系、CO2脱气运移机制、CO与常规烃类油气的耦合差异成藏机制等方面仍需要进一步的研究和探索。     

(ϖμ/ϖT) P of the lower mantle

We estimate (/T) _P of the lower mantle at seismic frequencies using two distinct approaches by combining ambient laboratory measurements on lower mantle minerals with seismic data. In the first approach, an upper bound is estimated for |(/T) _P | by comparing the shear modulus () profile of PREM with laboratory room-temperature data of extrapolated to high pressures. The second approach employs a seismic tomography constraint ( lnV _S / lnV _P ) _P =1.8–2, which directly relates (/T) _P with (K _S /T) _P . An average (K _S /T) _P can be obtained by comparing the well-established room-temperature compression data for lower mantle minerals with theK _S profile of PREM along several possible adiabats. Both (K _S /T) and (/T) depend on silicon content [or (Mg+Fe)/Sil of the model. For various compositions, the two approaches predict rather distinct (/T) _P vs. (K _S /T) _P curves, which intersect at a composition similar to pyrolite with (/T) _P =–0.02 to –0.035 and (K _S /T) _P =–0.015 to –0.020 GPa/K. The pure perovskite model, on the other hand, yields grossly inconsistent results using the two approaches. We conclude that both vertical and lateral variations in seismic velocities are consistent with variation due to pressure, temperature, and phase transformations of a uniform composition. Additional physical properties of a pyrolite lower mantle are further predicted. Lateral temperature variations are predicted to be about 100–250 K, and the ratio of ( lnp/ lnV _S ) _P around 0.13 and 0.26. All of these parameters increase slightly with depth if the ratio of ( lnV _S / lnV _P ) _P remains constant throughout the lower mantle. These predicted values are in excellent agreement with geodynamic analyses, in which the ratios ( ln / lnV _S ) _P and ( / lnV _S ) _P are free parameters arbitrarily adjusted to fit the tomography and geoid data.     

Dynamics of the mid-ocean ridge plate boundary: Recent observations and theory

The global mid-ocean ridge system is one of the most active plate boundaries on the earth and understanding the dynamic processes at this plate boundary is one of the most important problems in geodynamics. In this paper I present recent results of several aspects of mid-ocean ridge studies concerning the dynamics of oceanic lithosphere at these diverging plate boundaries. I show that the observed rift valley to no-rift valley transition (globally due to the increase of spreading rate or locally due to the crustal thickness variations and/or thermal anomalies) can be explained by the strong temperature dependence of the power law rheology of the oceanic lithosphere, and most importantly, by the difference in the rheological behavior of the oceanic crust from the underlying mantle. The effect of this weaker lower crust on ridge dynamics is mainly influenced by spreading rate and crustal thickness variations. The accumulated strain pattern from a recently developed lens model, based on recent seismic observations, was proposed as an appealing mechanism for the observed gabbro layering sequence in the Oman Ophiolite. It is now known that the mid-ocean ridges at all spreading rates are offset into individual spreading segments by both transform and nontransform discontinuities. The tectonics of ridge segmentation are also spreading-rate dependent: the slow-spreading Mid-Atlantic Ridge is characterized by distinct bulls-eye shaped gravity lows, suggesting large along-axis variations in melt production and crustal thickness, whereas the fast-spreading East-Pacific Rise is associated with much smaller along-axis variations. These spreading-rate dependent changes have been attributed to a fundamental differences in ridge segmentation mechanisms and mantle upwelling at mid-ocean ridges: the mantle upwelling may be intrinsically plume-like (3-D) beneath a slow-spreading ridge but more sheet-like (2-D) beneath a fast-spreading ridge.     

南海海盆区莫霍面分布规律及其对深部钻探的意义

钻遇莫霍面是人类一直以来的梦想。深海海底是地球上离莫霍面最近的地方，目前有研究推测南海是世界上莫霍面深度最浅的海域之一，但缺乏足够的直接证据。深反射地震探测可以直接揭示岩石圈的构造形态，是莫霍面探测的重要手段。本文基于长达15000 km的深反射多道地震剖面的解释、处理、制图和分析，结合前人的研究，形成了南海海盆区莫霍面反射特征和空间分布的初步认识。① 南海东部次海盆南部早期经历了较快速扩张，岩浆供应充足，受扩张停止后岩浆活动影响较小，基底平坦，地质构造相对简单，同时洋壳地震速度结构不存在异常，且有较强的广角莫霍面反射波和可识别的地幔顶部折射波，具备莫霍面钻探的基本条件。② 南海海盆不同区域的莫霍面反射强度存在较大差异。其中东部次海盆莫霍面反射最为强烈且清晰，西北次海盆次之，西南次海盆仅有零星出现的清晰莫霍面反射且可信度不高。③ 识别南海海盆区莫霍面地震反射长度超过3500 km，首次形成了海盆区深度域莫霍面地震反射空间分布图。与重力反演的莫霍面深度相比，利用深反射多道地震计算的莫霍面深度细节更为丰富，并且可以在垂向上清晰刻画莫霍面的结构。整体上，南海海盆区莫霍面地震反射强烈和可信度高的区域中，深度较浅的区域之一是东部次海盆南部，最浅处仅约9. 5 km，其中水深4. 01 km，洋壳厚度仅5. 54 km。综合判断，东部次海盆南部是南海重要的莫霍面钻探备选区，这对南海莫霍面钻探选址具有重要意义。     

施工顺序对软土地坪复合地基的影响探讨

场地回填土与搅拌桩施工顺序将直接影响软土地坪复合地基的处理效果,其中沉降控制是影响地坪安全运行的关键。以近海软土地区工业厂房水泥土搅拌桩复合地基工程为例,探讨了复合地基水泥土搅拌桩、回填土两种不同施工顺序的影响因素及相应处理效果;通过数值模拟试验,计算分析了不同填土厚度下两种工况的地基固结总沉降、施工沉降、工后沉降。结果表明,采用“先土后桩”的施工顺序的施工期沉降大于“先桩后土”,但工后沉降大大减小,提高了软土地坪地基处理效果,有利于处理后地坪的安全运行,验证了前述分析结论,为合理设计方案的确定提供了重要的参考依据。     

短周期仪器记录在远震分析中的应用

地震泼由多种频率成份组成,在其传播过程中。由于地球介质的滤波作用,随着传播距离的增大。高频成份逐渐被衰减。周期越大的成份传播越远。因此,在地震观测上,不同频带的地震仪用于观测不同震中距范围的地震,使仪器较好地响应,就我国台网而言,短周期仪用于地方震、近震的监测,中长周期仪和长周期仪用于远震,极远震的监测,其中长周期仪偏重于记录极远震。但作者在实际分析工作中注意到,短周期仪的记录应用于远震,极远震的分析,显示出独特的作用,充分挖掘和利用短仪资料,将会更加丰富远震、极远震的震相资料。本文以高台地震台的资料为依据,从震中距和震级角度,与中长仪和长仪作比较,对短仪记录远震、极远震的震相及特征作了统计分析与初步探讨,并对新疆地区“影区”地震S波的记录、日本地区地震_PPcP等震相的出现提出了现象的存在,给出了相应的观测结果,以待步一步研究。     

南极长城站地区现代环境地球化学特征

本文通过环境背景值、风化壳地球化学、对流层(大气气溶胶)地球化学和人为地球化学异常,初步探讨了南极长城站地区的现代环境地球化学特征。分析表明:环境要素固有的地球化学性质、区域环境条件和自然环境演变之间具有深刻的内在联系。     

云南大地热流及地热地质问题

本文报道了云南20多个大地热流测试数据。据热流分布特点并参考前人有关资料,结合地球物理测深资料、岩石热物性与放射性生热率测试数据,对各区地壳深度范围的温度、热流配分组成作了粗略分析,划分出各有特征的若干地质-热流区。云南普遍地热偏高(滇东南小范围除外)的主要原因是新构造运动强烈;怒江以西近代岩浆活动亦强烈,具有典型的板块交汇带的基本地热特点     

Na depletion in modern adakites via melt / rock reaction within the sub-arc mantle

Interaction between slab-derived melt and mantle peridotite and the role of slab melt as a metasomatizing agent in the sub-arc mantle is being increasingly recognized. Adakite, the slab melt erupted on the surface, usually exhibits anomalously high MgO, CaO, Cr and Ni contents that indicate interaction with mantle peridotitite. Here we note that Cenozoic adakites have Na₂O contents below 5.8 wt.% with ∼95% samples lower than 5.0 wt.%, and are generally depleted in this component relative to experimental basalt partial melts (mostly beyond 5.0 wt.% and up to 9.0 wt.% Na₂O) produced under 1.5-3.0 GPa conditions that are most relevant to adakite production. We interpret the adakite Na depletion to be also a consequence of the melt / rock reaction that takes place within the hot mantle wedge. During ascent and reaction with mantle peridotite, primary adakite melts gain mantle components MgO, CaO, Cr and Ni but lose Na₂O, SiO₂ and perhaps K₂O to the mantle, leading to Na-rich mantle metasomatism. Selective assimilation of predominately mantle clinopyroxene, some spinel and minor olivine at high T/P has been considered to be an important process in producing high-Mg adakites from primary low-Mg slab melts [Killian, R., Stern, C. R., 2002. Constraints on the interaction between slab melts and the mantle wedge from adakitic glass in peridotite xenoliths. Eur. J. Mineral. 14, 25-36]. In such a process, Na depletion in the assimilated melt is the result of dilution due to the increase in melt mass. Phase relationships in the reaction system siliceous melt + peridotite and quantitative calculation suggest that assimilation of mantle clinopyroxene, olivine and spinel and fractional crystallization of sodic amphibole and orthopyroxene, under conditions of moderate T/P and increasing melt mass, is also an important process that modifies the composition of adakites and causes the Na depletion.