Lithosphere thermal structure at the eastern margin of the Bohemian Massif: a case petrological and geophysical study of the Nied?wied? amphibolite massif (SW Poland)

The crustal section beneath amphibolite Nied?wied? Massif (Fore-Sudetic Block in NE Bohemian Massif), modelled on the basis of geological and seismic data, is dominated by gneisses with subordinate granites (upper and middle crust) and melagabbros (lower crust). The geotherm was calculated based on the chemical analyses of the heat-producing elements in the rocks forming the crust and the measurements of their density and heat conductivity. The results were verified by heat flow calculations based on temperature measurements from 1,600?m deep well in the Nied?wied? Massif and by temperature–depth estimates in mantle xenoliths coming from the nearby ca. 4.5?My basanite plug in Lutynia. The paleoclimate-corrected heat flow in the Nied?wied? Massif is 69.5?mW?m^?2, and the mantle heat flow is 28?mW?m^?2. The mantle beneath the Massif was located marginally relative to the areas of intense Cenozoic thermal rejuvenation connected with alkaline volcanism. This results in geotherm which is representative for lithosphere parts located at the margins of zones of continental alkaline volcanism and at its waning stages. The lithosphere–asthenosphere boundary (LAB) beneath Nied?wied? is located between 90 and 100?km depth and supposedly the rheological change at LAB is not related to the appearance of melt.     

A review of the Irish crustal structure and signatures from the Caledonian and Variscan Orogenies

This paper reviews the complex crustal and upper-mantle seismic velocity structure of Ireland and surrounding seas. Data from 11 seismic refraction profiles reveal that onshore Ireland mean crustal velocities range between 6.25 and 6.5 km s⁻¹ with crustal thickness of 28.5–32 km. Superimposed on a three-layer crust, the sedimentary layer has a thickness of approximately 6–8 km at the southern coastline, but only 3–4 km in the vicinity of the Shannon Estuary in western Ireland. The lateral heterogeneity of the upper-crustal layer is pervasive throughout Ireland, with velocities of 5.7–6.2 km s⁻¹ and a layer thickness of 3–10 km. A low-velocity zone is found in the south-east which is interpreted as the buried south-western extension of the Leinster Granite. The mid-crustal layer (6.3–6.7 km s⁻¹) is between 8 and 16 km thick. Significant changes occur in the vicinity of the Shannon Estuary, around the location of the Iapetus Suture Zone. The lower crust is fairly uniform with velocities of 6.8–7.2 km s⁻¹ and a thickness of approximately 8–10 km except towards the south of Ireland where the Moho appears as a transition zone. Offshore Ireland, a two-layer crust with a thickness of 24–26 km beneath the North Celtic Sea Basin and only 14–15 km beneath the Rockall Trough prevails.     

大陆地震构造系统: 以青藏高原及邻区为例

青藏高原及邻区三角形发震构造域是全球大陆最显著的地震多发区.脆性活动断层及其弹性回跳模式无法合理解释该区深度集中分布在10~40 km的点状震源.针对发震构造和地震机理不明确这一重大科学问题, 以大陆动力学和地球系统动力学新思想为指导, 对青藏高原及邻区发震构造系统进行域、层、带、点相关研究, 阐明大陆地震构造系统的结构型式, 认为下地壳固态流变及其韧性剪切带是提供地震能量的孕震构造, 中地壳韧-脆性剪切带是累积地震能量的发震构造, 上地壳脆性断裂是释放地震能量的释震构造.在研究青藏高原及邻区地震构造系统及其形成背景的基础上, 进一步论证了大陆地震热流体撞击的形成机理: 地幔墙导致大洋中脊之下的软流圈热流物质层流到大陆特定部位汇聚加厚并底辟上升, 造成大陆下地壳部分熔融和固态流变, 并改变莫霍面的产状, 固态流变物质侧向非均匀流动, 形成大陆盆山体系, 流动的韧性下地壳与脆性上地壳之间具有韧-脆性剪切滑脱性质的中地壳不断积累由下地壳热能转换而来的应变能, 形成发震层, 震源定位于下地壳热流物质富集带("热河")中的固态-半固态流变物质撞击到强弱层块之间的构造边界, 不同热构造环境和撞击角度产生5种不同类型的地震.从而为大陆地震的科学预测奠定了全新的理论基础.      

Thermal conductivity of dense hcp iron:Direct measurements using laser heated diamond anvil cell

Thermal conductivity(k) of iron is measured up to about 134 GPa.The measurements are carried out using the single sided laser heated diamond anvil cell,where the power absorbed by a Fe metal foil at hotspot is calculated using a novel thermodynamical method.Thermal conductivity of fee(γ)-Fe increases up to a pressure of about46 GPa.We find thermal conductivity values in the range of 70-80 Wm~(-1)K~(-1)(with an uncertainty of 40%),almost constant with pressure,in the hcp(e) phase of Fe.We attribute the pressure independent k above 46 GPa to the strong electronic correlation effects driven by the electronic topological transition(ETT).We predict a value of thermal conductivity of ε-Fe of about 40±16 Wm~(-1)K~(-1) at the outer core of Earth.     

花岗岩浆形成定位机制的思考与研究进展

花岗岩(广义)是陆壳的标志,也是地球岩石圈区别于其它行星岩石圈的标志。文章介绍了行星探测和大洋调查等方面的成果对花岗岩形成的地质约束:行星从岩浆表壳向岩石表壳转换过程以及现代地幔过程,均没有产生有规模意义的花岗岩;花岗岩及其所标志的陆壳,应是星球出现水圈和沉积岩之后的产物;花岗岩在地球岩石圈二维空间上的平均生长速率,大约为485×10~3km~2/Myr;岩浆主要来自地壳岩石的部分熔融(深熔)。在此基础上,文章介绍了深熔作用方面的研究进展,讨论了部分熔融岩石的流变行为与其内熔体比的关系,并比较了岩浆侵入模型与岩浆对流模型在解释花岗岩形成定位机制方面的异同。侵入模型的困难之一来自岩体与源区分离。由于源区位于岩体下方且远离岩体,因而是不可观察的,除非岩体及其与源区之间的岩石因风化或构造被剥蚀殆尽。文章最后介绍了"深熔-对流"模型的研究进展。该模型认为"源区"与"定位区间"是统一的,当"源区"岩石的熔体比例超过流变学的临界熔体比,岩石转变为"脏"岩浆;"脏"岩浆层内的重力分异诱发热对流,后者引起"顶蚀作用",导致重熔界面(MI)或固-液转换界面(SLT)不断向上移动和岩浆层的逐渐增厚。基本认识是:熔区内的热对流是深熔作用能够形成大规模花岗岩浆的必要条件;没有对流,陆壳岩石的部分熔融只能产生混合岩,不能产生岩基规模的花岗岩。     

Microstructural and experimental constraints on the rheology of partially molten gabbro beneath oceanic spreading centers

Flow laws for high-temperature creep of olivine, plagioclase, and diabase are used to place constraints on the rheology of partially molten lower oceanic crust. This analysis is motivated by the observation of olivine lattice preferred orientations and subgrain microstructures in oceanic gabbros that lack evidence for dislocation creep in coexisting plagioclase and pyroxene. Extrapolation of experimental flow laws indicates that at temperatures above 1100°C and stresses less than 10 MPa, olivine may be the weakest phase in rocks with gabbroic composition. By accounting for variations in the melt fraction (0–10%) and grain size of partially molten plagioclase aggregates we can constrain the rheological conditions where olivine deforms by dislocation creep while plagioclase deforms by diffusion creep. Calculated effective viscosities range from 10¹⁵ to 10¹⁹ Pa s; based on observations of the geometry of the partially molten zone beneath the East Pacific Rise and the microstructural and experimental constraints we favor a value of ∼10¹⁸ Pa s. This value approaches estimates for the viscosity of the upper mantle beneath ridge axes, but is significantly higher than previously suggested for the partially molten lower crust. Such high viscosities are inconsistent with ridge evolution models that require large amounts of lower crustal flow to accommodate melt redistribution. However, the results are compatible with recent models that favor local magma replenishment from the mantle at closely spaced intervals along the spreading center axis in a 2D, ‘sheet-like’ fashion.     

Rheological structure and dynamical response of the DSS profile BALTIC in the SE Fennoscandian Shield

Numerical modelling was applied to study the present-day state of stress and deformation under different tectonic loading conditions at the seismic BALTIC–SKJ profile in south-eastern Finland and in Estonia. The finite element method was used to solve the numerical problem. The two-dimensional model was constructed using the results from both seismic and thermal studies along the profile. The model is 700 km long and 200 km deep, and is roughly divided into an inhomogeneous, laterally layered crust and a homogeneous mantle lithosphere. Both the linear elastic and non-linear elasto-plastic rheologies were used. Elasto-plasticity was achieved by calculating a rheological strength as a function of depth along the profile. Different tectonic load cases were analysed with displacement, force and pressure type boundary conditions. Also, the effect of different strain rates was investigated. The results suggest that even with relatively low compressive stress levels the lower crust deforms in a plastic manner for a wet crustal rheology. When applying a dry crustal rheology, plastic yielding is attained only with much higher stress fields.     

Mantle derived crystal-poor rhyolitic ignimbrites: Eruptive mechanism from geochemical and geochronological data of the Piedra Parada caldera,Southern Argentina

Trace elements, isotopic modeling and U-Th-Pb SHRIMP zircon age constraints are used to reconstruct the eruption history and magmatic processes of the Piedra Parada Caldera. In the early Eocene, the crystal-poor Barda Colorada ignimbrite(BCI), having 15% micro-porphyritic crystals with respect to magmatic components, erupted a volume estimated in more than 300 km~3. The Piedra Parada caldera is located in the Patagonian Andes foreland, at the southern end of the calderas field of the Pilcaniyeu Volcanic Belt(PVB). This belt is related to an extensional tectonic setting as a result of the collision of the Farallon-Aluk ridge with South America, which enabled the development of a transform ocean/continental plate margin followed by the detachment of the Aluk plate and the opening of a slab window. The BCI extra-caldera Plateau is a 100 m thick deposit, having a lower unit with high silica(Si O_2 76 wt.%),potassium poor rhyolitic composition(trondhjemitic like magma), and an upper unit with normal to high potassium rhyolitic composition(granitic like magma). A trace elements modeling of the BCI units shows that the BCI lower and upper units did not evolve from fractionation or immiscibility in the shallow magma reservoir. The BCI also have a primitive isotopic signature(initial87 Sr/86 Sr =0.7031-0.7049 and ε_(Nd)= +3.4 to +3.65). Thus, tectonic, compositional and isotopic constraints suggest the fast ascent of high silica magmas to a shallow reservoir, and point to an upper mantle origin for these rhyolitic magmas in a transitional(Orogenic-Anorogenic) tectono-magmatic setting. U-Th-Pb SHRIMP zircon crystallization ages of the Syn-caldera stage BCI units(56 -51.5 Ma) show a protracted life of 5 Ma for this caldera reservoir. The age of 52.9 ± 0.3 Ma is considered the best fit for the possible maximum age for the caldera collapse. The Late-caldera magmatism has trachyandesitic and rhyolitic compositions.The trace element modeling suggests that these rhyolites evolve from the trachyandesites and do not evolve from the BCI residual magma. The trachyandesites have U-Th-Pb SHRIMP zircon crystallization ages of 52 ± 1 Ma, suggesting that the caldera eruption was triggered by the arrival of the trachyandesitic magma.     

Soil volatile mercury,boron and ammonium distribution at Cañadas caldera,Tenerife, Canary Islands,Spain

Spatial distribution of soil Hg, B and NH₄ was investigated in the soils of Cañadas caldera, Canary Islands, in the summer of 1992. Soil Hg, B and NH₄ were also studied over several transects intersecting structural features. Soil Hg concentration ranged from 10.83 to 45,000 μg kg⁻¹ whereas those of B and NH₄ ranged from 8.45 to 4512 μg kg⁻¹ and from 0.31 to 181 mg kg⁻¹, respectively. Probability plot analysis identified 3 geochemical populations for Hg and B whereas only 2 were detected for NH₄. Multivariate analysis (cluster, factor and multiple regression analysis) was used to identify subtle geochemical characteristics of the groups and factors and to determine the relationship between Hg, B and NH₄ and the soil secondary parameters. High soil Hg, B and NH₄ concentrations are interpreted as being related to areas where a convection system has developed, coinciding with the most recent volcanic centers located along the basaltic rift zones and with the main features inside Cañadas caldera: Teide volcano and Roques de Garcı́a. Background levels of soil volatiles are present in locations where no indication of subsurface thermal activity occurs. Generally, the secondary controls on Hg, B and NH₄ are subtle and are overwhelmed in areas characterized by subsurface geothermal activity.     

http://www.sciencedirect.com/science/article/pii/S1674987113001072

The late Permian Emeishan large igneous province(ELIP) covers ~0.3 x 106 km2 of the western margin of the Yangtze Block and Tibetan Plateau with displaced,correlative units in northern Vietnam(Song Da zone).The ELIP is of particular interest because it contains numerous world-class base metal deposits and is contemporaneous with the late Capitanian(~260 Ma) mass extinction.The flood basalts are the signature feature of the ELIP but there are also ultramafic and silicic volcanic rocks and layered maficultramafic and silicic plutonic rocks exposed.The ELIP is divided into three nearly concentric zones(i.e.inner,middle and outer) which correspond to progressively thicker crust from the inner to the outer zone.The eruptive age of the ELIP is constrained by geological,paleomagnetic and geochronological evidence to an interval of 3 Ma.The presence of picritic rocks and thick piles of flood basalts testifies to high temperature thermal regime however there is uncertainty as to whether these magmas were derived from the subcontinental lithospheric mantle or sub-lithospheric mantle(i.e.asthenosphere or mantle plume) sources or both.The range of Sr(I_(Sr) = 0.7040-0.7132),Nd(ε_(Nd)(t) ≈-14 to +8),Pb(~(206)Pb/~(204)Pb_1≈ 17.9-20.6) and Os(γ_(Os) =-5 to +11) isotope values of the ultramafic and mafic rocks does not permit a conclusive answer to ultimate source origin of the primitive rocks but it is clear that some rocks were affected by crustal contamination and the presence of near-depleted isotope compositions suggests that there is a sub-lithospheric mantle component in the system.The silicic rocks are derived by basaltic magmas/rocks through fractional crystallization or partial melting,crustal melting or by interactions between mafic and crustal melts.The formation of the Fe-Ti-V oxide-ore deposits is probably due to a combination of fractional crystallization of Ti-rich basalt and fluxing of C02-rich fluids whereas the Ni-Cu-(PGE) deposits are related to crystallization and crustal contamination of mafic or ultramafic magmas with subsequent segregation of a sulphide-rich portion.The ELIP is considered to be a mantle plume-derived LIP however the primary evidence for such a model is less convincing(e.g.uplift and geochemistry) and is far more complicated than previously suggested but is likely to be derived from a relatively short-lived,plume-like upwelling of mantle-derived magmas.The emplacement of the ELIP may have adversely affected the short-term environmental conditions and contributed to the decline in biota during the late Capitanian.     

秦岭造山带根部地壳结构及流变学演化

在地质和地球物理资料基础上，运用物理学观点，研究和分析了秦岭造山带根部－大别前寒武纪变质地体的三维结构及流变学演化历史，通过现代及中－新元古代时期地壳流变学剖面的构筑，强调地壳流变学分层性及变质变形分解作用对中下地壳结构及地球动力学演化的控制作用，线状强应变带与透镜状弱应变域的规律组合，是秦岭造山带及其根部地壳结构的基本样式，并具尺度不变性，将古老中下地壳近３Ｇａ的流变学演化历史划分为７个阶段。它     

Age and metal source of orogenic gold deposits in Southeast Guizhou Province,China: Constraints from Re-Os and He-Ar isotopic evidence

The orogenic gold deposits in Southeast Guizhou are an important component of the Xuefeng polymetallic ore belt and have significant exploration potential, but geochronology research on these gold deposits is scarce. Therefore, the ore genetic models are poorly constrained and remain unclear. In the present study, two important deposits(Pingqiu and Jinjing) are investigated, including combined Re-Os dating and the He-Ar isotope study of auriferous arsenopyrites. It is found that the arsenopyrites from the Pingqiu gold deposit yielded an isochron age of 400 ± 24 Ma,with an initial ~(187)Os/~(188)Os ratio of 1.24 ± 0.57(MSWD = 0.96). An identical isochron age of 400 ± 11 Ma with an initial ~(187)Os/~(188)Os ratio of 1.55 ± 0.14(MSWD = 0.34) was obtained from the Jinjing deposit. These ages correspond to the regional Caledonian orogeny and are interpreted to represent the age of the main stage ore. Both initial ~(187)Os ratios suggest that the Os was derived from crustal rocks. Combined with previous rare earth element(REE), trace elements, Nd-Sr-S-Pb isotope studies on scheelite, inclusion fluids with other residues of gangue quartz, and sulfides from other gold deposits in the region, it is suggested that the ore metals from Pingqiu and Jinjing were sourced from the Xiajiang Group. The He and Ar isotopes of arsenopyrites are characterized by ~3 He/~4 He ratios ranging from 5.3 × 10~(-4) Ra to 2.5 × 10~(-2) Ra(Ra = 1.4 × 10~(-6), the ~3 He/~4 He ratio of air), 40 Ar=/~4 He ratios from 0.64 × 10~(-2) to 15.39×10~(-2), and ~(40)Ar/~(36)Ar ratios from 633.2 to 6582.0. Those noble gas isotopic compositions of fluid inclusions also support a crustal source origin,evidenced by the Os isotope. Meanwhile, recent noble gas studies suggest that the amount of in situ radiogenic ~4 He generated should not be ignored, even when Th and U are present at levels of only a few ppm in host minerals.     

From back-analysis to run-out prediction: a case study in the Western Italian Alps

Rock avalanches are complex phenomena that occur with a low frequency but which have a high destructive potential. As a consequence, the people who are responsible for the management of a territory are more and more interested in predicting the possible evolutions of well-known potential events. Tackling the above problems from a quantitative point of view, the RASH3D code, based on continuum mechanics concepts, has been here used to predict the evolution of a potential rock avalanche in the Western Italian Alps. A calibration-based approach, in which rheological parameters are constrained by systematic adjustment during trial-and-error back-analysis of past events similar to the landslide under investigation, is proposed to set rheological parameter values to be used for prediction purposes. The back-analysis of a 2?10⁶ m³ rock avalanche located in the Divedro Valley, close to the area of the potential event, has then been analysed using both a frictional and a Voellmy rheology. The characteristics of the slope and the dynamics of the event have made the frictional rheology more suitable to come to the correct simulation of the historical case. The back-analysis results have contributed not only in the selection of the rheological parameter values but also in the choice of the type of rheological law to use in the carried out forward-analyses.     

A new tectonic model for the Laurentia−Avalonia−Baltica sutures in the North Sea: A case study along MONA LISA profile 3

We present a new model for the lithospheric structure of the transitions between Laurentia, Avalonia and Baltica in the North Sea, northwestern Europe based on 2¾D potential field modelling of MONA LISA profile 3 across the Central Graben, with constraints from seismic P-wave velocity models and the crustal normal incidence reflection section along the profile. The model shows evidence for the presence of upper-and lower Palaeozoic sedimentary rocks as well as differences in crustal structure between the palaeo-continents Laurentia, Avalonia and Baltica. Our new model, together with previous results from transformations of the gravity and magnetic fields, demonstrates correlation between crustal magnetic domains along the profile and the terrane affinity of the crust. This integrated interpretation indicates that a 150 km wide zone, characterized by low-grade metamorphosis and oblique thrusting of Avalonia crust over Baltica lower crust, is characteristic for the central North Sea area. The magnetic susceptibility and the density across the Coffee Soil Fault range from almost zero and 2715 kg/m³ in Avalonia crust to 0.05 SI and 2775 kg/m³ in Baltica crust. The model of MONA LISA profile 3 indicates that the transition between Avalonia and Baltica is located beneath the Central Graben with a ramp–flat–ramp geometry. Our results indicate that the initial rifting of the Central Graben and the Viking Graben was controlled by the location of the Caledonian collisional suture, located at the Coffee Soil Fault, and that the deep crustal part of Baltica extends further to the west than hitherto believed.     

Crustal Structure of the Crimean Mountains along the Sevastopol–Kerch Profile from the Results of DSS and Local Seismic Tomography

The paper discusses the velocity structure of the crust beneath the Crimean Mountains from the results of active and passive seismic experiments. Based on a new interpretation of seismic data from the old Sevastopol–Kerch DSS profile by modern full-wave seismic modeling methods, a velocity model of the crust beneath the Crimean Mountains has been constructed for the first time. This model shows the significant differences in the structure of two crustal blocks: (1) one characterized by higher velocities and located in the western and central Crimean Mountains, and (2) the other characterized by lower velocities and located in the east, in the Feodosiya–Kerch zone, which are subdivided by a basement uplift (Starokrymskoe Uplift). The former block is characterized by a more complex structure, with the Moho traced at depths of 43 and 55 km, forming two Moho discontinuities: the upper one corresponds to the platform stage, and the lower one, formed presumably at the Alpine stage of tectogenesis as a result of underthrusting of the East Black Sea microplate beneath the southern margin of the Scythian Plate in Crimea. At depths of 7–11 km, velocity inversion zone has been identified, indicating horizontal layering of the crust. Local seismic tomography using the data on weak earthquakes (m_b ≤ 3) recorded by the Crimean seismological network allowed us to obtain data on the crustal structure beneath the Crimean Mountains at depths of 10–30 km. The crustal structure at these depths is characterized by the presence of several high-velocity crustal bodies in the vicinity of cities Yalta, Alushta, and Sudak, with earthquake hypocenters clustered within these bodies. Comparison of this velocity model of the Crimean Mountains with the seismicity distribution and with the results from reconstruction of paleo- and present-day stress fields from field tectonophysical study and earthquake focal mechanisms allowed the conclusion that the Crimean Mountains were formed as a result of on mature crust at the southern margin of the East European Platform and Scythian Plate, resulting from processes during various phases of Cimmerian and Alpine tectogenesis in the compressional and transpressional geodynamic settings. The collisional process is ongoing at the present-day stage, as supported by high seismicity and uplift of the Crimean Mountains.     

青海共和盆地岩石圈热-流变结构及地热意义

地球深度热状况是深部地球动力学和岩石圈活动性研究的重要内容, 岩石圈热结构和热-流变结构可以很好地揭示岩石圈范围内的热状况。近年来, 在青海共和盆地钻探揭露了深部高温干热岩体, 关于其热源机制尚未有定论。本文以青海共和盆地为研究对象, 分析壳内温度分布和流变强度, 探讨壳内低速体的地质属性。结果表明, 共和盆地的地壳流变结构从上而下分为脆性和韧性两层, 韧性层又包括中地壳和下地壳两层韧性层, 在上地壳尺度均表现为脆性破裂为主, 并逐渐过渡为韧性流变; 恰卜恰地区在脆性破裂的上地壳延伸至中下地壳时, 破裂沿一系列滑脱面发生韧性滑动, 局部地段形成壳内熔融, 为恰卜恰地区提供了额外的热源, 使其大地热流值(109.6 mW/m2)显著高于贵德地区(77.6 mW/m2)。这一认识为共和盆地壳内低速体存在提供了新的佐证, 也为区内干热岩热源分析以及高温地热资源探测开发提供了科学依据。     

Seismic structure and rheology of the crust under mainland China

The crust and upper mantle in mainland China were relatively densely probed with wide-angle seismic profiling since 1958, and the data have provided constraints on the amalgamation and lithosphere deformation of the continent. Based on the collection and digitization of crustal P-wave velocity models along related wide-angle seismic profiles, we construct several crustal transects across major tectonic units in mainland China. In our study, we analyzed the seismic activity, and seismic energy releases during 1970 and 2010 along them. We present seismogenic layer distribution and calculate the yield stress envelopes of the lithosphere along the transects, yielding a better understanding of the lithosphere rheology strength beneath mainland China. Our results demonstrate that the crustal thicknesses of different tectonic provinces are distinctively different in mainland China. The average crustal thickness is greater than 65 km beneath the Tibetan Plateau, about 35 km beneath South China, and about 36–38 km beneath North China and Northeastern China. For the basins, the thickness is ~ 55 km beneath Qaidam, ~ 50 km beneath Tarim, ~ 40 km beneath Sichuan and ~ 35 km beneath Songliao. Our study also shows that the average seismic P-wave velocity is usually slower than the global average, equivalent with a more felsic composition of crust beneath the four tectonic blocks of mainland China resulting from the complex process of lithospheric evolution during Triassic and Cenozoic continent–continent and Mesozoic ocean–continent collisions. We identify characteristically different patterns of seismic activity distribution in different tectonic blocks, with bi-, or even tri-peak distribution of seismic concentration in South Tibet, which may suggest that crustal architecture and composition exert important control role in lithosphere deformation. The calculated yield stress envelopes of lithosphere in mainland China can be divided into three groups. The results indicate that the lithosphere rheology structure can be described by jelly sandwich model in eastern China, and crème brulee models with weak and strong lower crust corresponding to lithosphere beneath the western China and Kunlun orogenic belts, respectively. The spatial distribution of lithospheric rheology structure may provide important constraints on understanding of intra- or inter-plate deformation mechanism, and more studies are needed to further understand the tectonic process(es) accompanying different lithosphere rheology structures.     

Geothermal structure of the eastern Pontides orogenic belt： Implications oceanic lithosphere Black Sea basin and the eastern for subduction polarity of Tethys

The numerical results of thermal modeling studies indicate that the lithosphere is cold and strong beneath the Black Sea basin.The thermal lithospheric thickness increases southward from the eastern Pontides orogenic belt（49.4 km） to Black Sea basin（152.2 km）.The Moho temperature increases from 367℃in the trench to 978℃in the arc region.The heat flow values for the Moho surface change between 16.4 mW m^-2 in the Black Sea basin and 56.9 mW m^-2 in the eastern Pontides orogenic belt. Along the southern Black Sea coast,the trench region has a relatively low geothermal potential with respect to the arc and back-arc region.The numerical studies support the existence of southward subduction beneath the Pontides during the late Mesozoic-Cenozoic.     

Magmatic underplating beneath the Rajmahal Traps: Gravity signature and derived 3-D configuration

The early Cretaceous thermal perturbation beneath the eastern continental margin of the Indian shield resulted in the eruption of the Rajmahal Traps. To understand the impact of the magmatic process that originated in the deep mantle on the lower crustal level of the eastern Indian shield and adjoining Bengal basin the conspicuous gravity anomalies observed over the region have been modelled integrating with available geophysical information. The 3-D gravity modelling has delineated 10–15 km thick high-density (ρ = 3.02 g/cm³) accreted igneous layer at the base of the crust beneath the Rajmahal Traps. Thickness of this layer varies from 16 km to the west of the Rajmahal towards north to about 12 km near Kharagpur towards south and about 18 km to the east of the Raniganj in the central part of the region. The greater thickness of the magmatic body beneath the central part of the region presents itself as the locus of the potential feeder channel for the Rajmahal Traps. It is suggested that the crustal accretion is the imprint of the mantle thermal perturbation, over which the eastern margin of the eastern Indian shield opened around 117 Ma ago. The nosing of the crustal accretion in the down south suggests the possible imprint of the subsequent magmatic intrusion along the plume path.     

Geochemical evolution of historical lavas from Askja Volcano, Iceland: Implications for mechanisms and timescales of magmatic differentiation

The mechanisms and the timescales of magmatic evolution were investigated for historical lavas from the Askja central volcano in the Dyngjufjöll volcanic massif, Iceland, using major and trace element and Sr, Nd, and Pb isotopic data, as well as ²³⁸U-²³⁰Th-²²⁶Ra systematics. Lavas from the volcano show marked compositional variation from magnesian basalt through ferrobasalt to rhyolite. In the magnesian basalt-ferrobasalt suite (5-10 wt% MgO), consisting of lavas older than 1875 A.D., ⁸⁷Sr/⁸⁶Sr increases systematically with increasing SiO₂ content; this suite is suggested to have evolved in a magma chamber located at ∼600 MPa through assimilation and fractional crystallization. On the other hand, in the ferrobasalt-rhyolite suite (1-5 wt% MgO), including 1875 A.D. basalt and rhyolite and 20th century lavas, ⁸⁷Sr/⁸⁶Sr tends to decrease slightly with increasing SiO₂ content. It is suggested that a relatively large magma chamber occupied by ferrobasalt magma was present at ∼100 MPa beneath the Öskjuvatn caldera, and that icelandite and rhyolite magmas were produced by extraction of the less and more evolved interstitial melt, respectively, from the mushy boundary layer along the margin of the ferrobasalt magma chamber, followed by accumulation of the melt to form separate magma bodies. Ferrobasalt and icelandite lavas in the ferrobasalt-rhyolite suite have a significant radioactive disequilibrium in terms of (²²⁶Ra/²³⁰Th), and its systematic decrease with magmatic evolution is considered to reflect aging, along with assimilation and fractional crystallization processes. Using a mass-balance model in which simultaneous fractional crystallization, crustal assimilation, and radioactive decay are taken into account, the timescale for the generation of icelandite magma from ferrobasalt was constrained to be <∼3 kyr which is largely dependent on Ra crystal-melt partition coefficients we used.