The color of pumice: case study on a trachytic fall deposit, Meidob volcanic field, Sudan

 This paper examines the cause of color variations of trachytic pumices which are essentially uniform in chemical composition and proposes a geological model for their formation. A pyroclastic sequence of distinct subunits with brown, buff, and black pumices was deposited during the 5000-B.P. eruption of a tuff ring in the central Meidob volcanic field (Sudan). Subunits of buff pumices locally contain minor amounts of streaky pumice with pale-gray and dark-gray domains. The combined evidence of petrographic studies, chemical analyses of whole pumices and groundmass separates, electron microprobe analyses, optical spectroscopy, transmission electron microscopy, and magnetic susceptibility measurements show that color variations of the pumice clasts are related to the size and distribution of Fe³⁺-rich oxide microcrysts. Buff pumice and light-gray domains of streaky pumice have a colorless, transparent groundmass with very few microcrysts. Dark-gray domains of streaky pumice contain abundant hematite and/or magnetite microcrysts visible in thin section within a transparent, colorless glass groundmass. The groundmass of the black pumice clasts is brown in thin section which is most likely caused by submicroscopic magnetite microcrysts. Brown pumice clasts have a mixed groundmass consisting of brown domains and domains with opaque microcrysts in transparent glass. Variations in the eruption dynamics have been inferred from lithological observations. Subunits of black pumices are related to eruption pulses with low magma discharge and high water/magma mass ratio, whereas subunits of buff pumice were deposited during eruption pulses with high magma discharge and low water/magma mass ratio. Brown pumices represent the top part of the magma body, and the initial stage of the eruption probably had a low magma discharge. Streaky pumices are interpreted as the product of syn-eruptive mixing of Fe³⁺-rich oxide microcryst-bearing magma and microcryst-free magma. Received: 3 February 1997 / Accepted: 28 July 1997     

Heat transfer and the onset of chaos in a spherical,axisymmetric, anelastic model of whole mantle convection

Abstract

An isochemical uni-phase model of whole mantle convection has been developed in terms of which factors influencing the onset of time dependent chaotic behavior may be assessed. The model is spherical but restricted in generality to the analysis of axisymmetric solutions. In this paper we have employed it to examine the impact of compressibility and sphericity on the nature and onset of time dependence. Particular attention has been given to an examination of the impact that the onset of time dependence has upon the power law relation that connects the heat transfer (represented by the Nusselt number) to the strength of the thermal forcing (represented by the Rayleigh number). In order to obtain these results very extensive numerical simulations were required and the results themselves should be rather useful in the context of models of the thermal history of the planet.     

Constraints on cooling rates and permeabilities of pumice in an explosive eruption jet from colour and magnetic mineralogy

In explosive volcanic eruptions, vesicular magma droplets, produced by fragmentation, are propelled into the atmosphere where they are chilled to form pumices. The thermal history of droplets and the permeability of their internal bubble networks determine how much they are deformed in the eruption jet, and hence what information pumices record about the state of the magma at fragmentation. We study these aspects of the `Minoan' plinian eruption of Santorini Volcano by quantifying the rate of oxidation reactions that took place when air entered the hot magma fragments. In our experiments white Minoan pumices were heated for minutes to hours between 600 and 850°C, either in air, or in an atmosphere with an oxygen fugacity at the Ni–NiO buffer. Pumices were unchanged by heating at Ni–NiO. Those heated in air often became pink to dark pink, depending on heating time, and their Curie temperatures, as determined by magnetic susceptibility measurements, increased. We use oxidation rates deduced from these experiments, in conjunction with calculations of the rate of conductive cooling and of the rate at which air can enter a pumice, to constrain the conditions experienced by pumices during the eruption. Natural Minoan pumices less than about 5 cm in radius are white, whereas larger ones often have white rims and pink interiors with Curie temperatures higher than those of white material. We infer that small pumices were cooled before being oxidized, and that oxidation of the interiors of large clasts mostly took place during flight, at temperatures within a few tens of degrees of magmatic values. White rims of large pumices, despite being permeable, were cooled before oxidation could occur. Permeability developed in the liquid state, but did not develop early enough, with respect to cooling, or was not large enough to allow extreme oxidation. We give measurements of pumice permeabilities that should be close to magmatic values.     

Heat transfer between a thermochemical plume channel and the surrounding mantle in the presence of horizontal mantle flow

Heat and mass transfer processes in the conduit of a thermochemical plume located beneath an oceanic plate far from a mid-ocean ridge (MOR) proceed under conditions of horizontal convective flows penetrating the plume conduit. In the region of a mantle flow approaching the plume conduit (in the frontal part of the conduit), the mantle material heats and melts. The melt moves through the plume conduit at the average velocity of flow v and is crystallized on the opposite side of the conduit (in the frontal part of the conduit). The heat and the chemical dope transferred by the conduit to the mantle flow are carried away by crystallized mantle material at the velocity v.The local coefficients of heat transfer at the boundary of the plume conduit are theoretically determined and the balance of heat fluxes through the side of the plume conduit per linear meter of the conduit height. The total heat generation rate, transmitted by the Hawaiian plume into the upper and lower mantle, is evaluated. With the use of regular patterns of heat transfer in the lower mantle, which is modeled on the horizontal layer, heated from below and cooled from above, the diameter of the plume source, the kinematic viscosity of the melt in the plume conduit, and the velocity of horizontal lower-mantle flows are evaluated and the dependences of the temperature drop, viscosity and Rayleigh number for the lower mantle on the diameter of the plume source are presented.     

Rock magnetic analyses as a tool to investigate diversity of drift pumice clasts: An example from Japan's 2021 Fukutoku-Oka-no-Ba (FOB) eruption

The crystallization of nanolites within magma chambers has recently raised a strong interest due to their impact on increasing melt viscosity and triggering magmatic eruptions. In 2021, the Fukutoku-Oka-no-Ba (FOB) underwater eruption produced large quantities of pumices that eventually formed rafts drifting at the surface of the ocean to the East coasts of Japan. Pumices collected along the shore shortly after grounding show various colors, microscopic and Raman analyses made by Yoshida et al. (Island Arc, 31, 1, 2022) revealed the presence of magnetite nanolites in some of them. In this study, we explore the magnetic properties of a batch of pumices of different colors from the FOB eruption, aiming to refine characterization of iron oxide nanolites. We used various analytical techniques such as SEM and FEG-SEM observations, EDS-X analyses, and rock magnetic experiments, including thermomagnetic analyses, hysteresis curves, coercivity analyses and FORC measurements. Our findings reveal that the iron oxides present in the FOB samples are Ti-magnetite, with minor amounts of Mg and Al. The magnetic crystals show a wide range of sizes, from extra small iron oxide nanolites (ESION) in the pumices with the lighter colors, to more bulky grains reaching the micrometer size in some of the dark color samples, significant diffusion is inferred in that case. Consequently, the magnetic characterization of iron oxide crystals within the Fukutoku-Oka-no-Ba pumices reveals varying stages of nucleation, dissolution, growth, and diffusion processes, providing evidences for the heterogeneous state of the magma during the eruption.     

Permeability measurements of natural and experimental volcanic materials with a simple permeameter: Toward an understanding of magmatic degassing processes

Permeability measurement of quenched volcanic porous materials is an important approach to understand permeability development and degassing of vesicular silicic magmas. In this study, we developed a gas permeameter to measure permeability of natural samples and experimental products. The permeameter has broad measurement ranges of pressure difference (10¹–10⁵ Pa) and gas-flow rate (10^− 9–10^− 5 m³/s). These ranges enable us to measure viscous permeability in the range of 10^− 17–10^− 9 m² for 1 centimeter-scale samples, using the Forchheimer equation, which includes the inertial effect of gas flow permeating through the samples. In addition, we improved the procedure for performing permeability measurements of mm-sized products of decompression experiments. Although a previous study reported the first permeability data for vesicular silicic glass products of decompression experiments, we found an overestimation in their permeability data due to problems in sample preparation, especially for very low permeability samples. Our improved measurements give lower permeability values than those of Takeuchi et al. (2005)(Takeuchi, S., Nakashima, S., Tomiya, A., Shinohara, H., 2005. Experimental constraints on the low gas permeability of vesicular magma during decompression. Geophys. Res. Lett., 32, L10312 doi:10.1029/2005GL022491).     

Isotopic feature and uranium dating of the volcanic rocks in the Okinawa Trough

1 Introduction The Okinawa Trough(OT) extends for ~1200 km between Taiwan and Kyushu, and is an incipient back-arc basin bounded by the Ryukyu Arc to the south and east, and by Diaoyudao Uplift Belt[1,2] (Tai- wan-Sinzi Belt[3]) to the north and west. Drastic faults and volcanic activities are the typical features in the OT. Generally, the OT is divided into three segments (north, middle and south) by the Tokara Strait and Kerama Gap[4]. The geological features, such as tec- tonics, m…     

Recycling of magmatic clasts during explosive eruptions: estimating the true juvenile content of phreatomagmatic volcanic deposits

The juvenile content of phreatomagmatic deposits contains both first-cycle juvenile clasts derived from magma at the instant of eruption, and recycled juvenile clasts, which were fragmented and first ejected by earlier explosions during the eruption, but fell back or collapsed into the vent. Recycled juvenile clasts are similar to accessory and accidental lithics in that they contribute no heat to further magma: water interaction, but previously no effective criteria have been defined to separate them from first-cycle juvenile clasts. We have investigated componentry parameters (vesicularity, clast morphology and extent of mud-coating) which, in specific circumstances, can distinguish between first-cycle juvenile clasts, involved in only one explosion, and such recycled juvenile clasts. Phreatomagmatic fall deposits commonly show gross grainsize and sorting characteristics identical to deposits of purely dry or magmatic eruptions. However the abundance of non-juvenile clasts in pyroclastic deposits is a sensitive indicator of the involvement of external water. If this component is calculated including recycled juvenile clasts with accidental and accessory clasts the contrast is even more striking. Data from a Holocene maar deposit in Taupo Volcanic Zone, New Zealand, suggest that the first-cycle juvenile component of the deposits is less than one-third of that determined by simple juvenile:lithic:crystal componentry.     

气煤弹性各向异性系数实验测试

气煤裂隙丰富,具有弹性各向异性,确定气煤弹性各向异性类型与大小,对煤田各向异性研究具有重要意义. 本文在对气煤煤样孔隙率测试的基础上,测量了顺煤层面、垂直于煤层面、与煤层面成45°三个方向上纵波速度（VP）和横波速度（VSH、VSV）,计算弹性常数. 对于具有垂向对称轴横向各向同性模型,各向异性系数平均值|ε|≤0.2、|δ|≤0.2和|γ|≤0.2,表明气煤是弱各向异性介质;对于横向对称轴横向各向异性模型,各向异性系数平均值ε(V)≈δ(V),按照Thomsen观点气煤各向异性是由椭圆裂隙引起弹性各向异性.本文还讨论了各向异性系数与孔隙率的关系,表明各向异性的大小与孔隙率之间有一定的关系. 因此当煤层厚度与其顶底板岩性不变或变化很小,可以认为气煤各向异性主要由裂隙引起的,气煤各向异性大小和方向能够反映煤的裂缝密度和方位.     

圆钢管混凝土截面延性系数

首先介绍了圆钢管混凝土中核心混凝土极限压应变的计算方法，通过参数分析，给出了核心混凝土极限压应变的限值。在此基础上，建立了截面极限曲率的计算公式，给出圆钢管混凝土截面延性系数的计算方法，并对有关影响因素进行分析，最后给出延性系数的建议值。研究进一步表明，圆钢管混凝土具有良好的抗震性能。     

Quantifying the effect of in-stream rock clasts on the retardation of heat along a stream

With the introduction of Distributed Temperature Sensing (DTS) into the field of hydrology, temperature has become a powerful tracer in both space and time. However, the interpretation of the observed temperature signal is often not straightforward due to its non-conservative behavior. The objective of this research is to explore and quantify the retardation of heat along a small first order stream, with the long-term objective of identifying different runoff mechanisms by using heat as a tracer.     

Site coefficients suitable to China site category

Introduction Site effect on strong earthquake ground motion is an important problem for seismologists and earthquake engineers. The direct approach to study site effect is to use data from a vertical array of seismometers, which is called borehole station…     

水平单轴应力与横波反射系数

地球物理勘探中,地下应力的估测是非常重要的.本文利用应力诱导的弱各向异性的各向异性参数, 借助已有的具有水平对称轴横向各向同性介质顾及高阶小量的横波反射系数近似公式,针对下介质为各向同性受单轴应力作用情况随应力的增加求取应力精度降低的问题,从横波反射系数近似公式的推导出发,联合SH波与SV 波反射系数给出新的改进的单轴应力与横波反射系数的直观表达式;计算并对比了受单轴应力作用的各向同性介质的横波近似和严格的反射系数,并利用理论得到的横波反射系数对应力进行了估值.结果表明,利用改进后横波反射系数与应力的直观表达式, 当应力增加时,应力估值精度有了明显的改善,这为利用横波反射数据直观估测应力提供理论依据.     

地震波反射系数偏导矩阵的精确解及对近似计算的精度分析

目前叠前反演方法大多是基于Zoeppritz方程近似式实现的,它仅适应于弱反射介质界面、中小角度(或小偏移距)的地震数据反演,不能满足勘探开发的地质需求.本文建立了基于zoeppritz方程精确求解反射系数的梯度矩阵,分析了矩阵特点和精度,为研究利用反射系数梯度精确解反演地震参数奠定了基础.     

大型浅水湖泊与大气之间的动量和水热交换系数——以太湖为例

湖泊水面与大气之间垂直方向的动量、水汽和热量通量与风速、湿度和温度梯度之间存在比例关系,因此在湖泊水-气相互作用研究中,这比例系数(交换系数)是关键因子.在以往的研究中,交换系数通常直接采用水面梯度观测法或海洋大气近地层的参数化方案进行计算.本文采用涡度相关系统和小气候系统仪器在太湖平台上直接观测的通量和气象要素,对上述交换系数(最小均方差原则)进行优化,结果为:动量交换系数C_D10N=1.52×10^-3、水汽交换系数C_E10N=0.82×10^-3、热量交换系数CH10N=1.02×10^-3,与其他内陆湖泊涡度相关观测数据的推导结果一致.本文的研究结果表明:与海洋参数化方案相比,在相同的风速条件下,湖面的空气动力学粗糙度比海洋高,这可能是由于受到水深的影响;如果采用海洋参数化方案,会导致湖泊年蒸发量的估算值偏大40%.太湖的动量、水汽和热量交换系数可以视为常数,可以不考虑稳定度和风速的影响.这是因为本文中83%的数据为近中性条件.敏感性分析表明:如果考虑稳定度的影响,LE模拟值的平均误差降低了0.5 W/m²,H的平均误差降低了0.4 W/m²,u*的计算值没有变化;如果考虑风速的影响,u*模拟值的平均误差降低了0.004 m/s,LE的平均误差升高了1.3 W/m²,H的模拟结果几乎不受影响.这一结果能为湖气相互作用研究提供参考.     

Heat flow in the Basin and Range province and thermal effects of tectonic extension

In regions of tectonic extension, vertical convective transport of heat in the lithosphere is inevitable. The resulting departure of lithosphere temperature and thickness from conduction-model estimates depends upon the mechanical mode of extension and upon how rapidly extension is (and has been) taking place. Present knowledge of these processes is insufficient to provide adequate constraints on thermal models. The high and variable regional heat flow and the intense local heat discharge at volcanic centers in the Basin and Range province of the United States could be accounted for by regional and local variations in extensional strain rate without invoking anomalous conductive heat flow from the asthenosphere. Anomalous surface heat flow typical of the province could be generated by distributed extension at average rates of about 1/2 to 1%/m.y., similar to rates estimated from structural evidence. To account for higher heat flow in subregions like the Battle mountain High, these rates would be increased by a factor of about 3, and locally at active bimodal volcanic centers, by an order of magnitude more.     

Heat transport in the Red Lake Bog,Glacial Lake Agassiz Peatlands

We report the results of an investigation on the processes controlling heat transport in peat under a large bog in the Glacial Lake Agassiz Peatlands. For 2 years, starting in July 1998, we recorded temperature at 12 depth intervals from 0 to 400 cm within a vertical peat profile at the crest of the bog at sub‐daily intervals. We also recorded air temperature 1 m above the peat surface. We calculate a peat thermal conductivity of 0·5 W m^?1 °C^?1 and model vertical heat transport through the peat using the SUTRA model. The model was calibrated to the first year of data, and then evaluated against the second year of collected heat data. The model results suggest that advective pore‐water flow is not necessary to transport heat within the peat profile and most of the heat is transferred by thermal conduction alone in these waterlogged soils. In the spring season, a zero‐curtain effect controls the transport of heat through shallow depths of the peat. Changes in local climate and the resulting changes in thermal transport still may cause non‐linear feedbacks in methane emissions related to the generation of methane deeper within the peat profile as regional temperatures increase. Copyright © 2006 John Wiley & Sons, Ltd.     

Ambiguous reflection coefficients for anelastic media

The calculation of reflection and transmission coefficients of plane waves at a plane interface between two homogeneous anelastic media may become ambiguous because it is not always obvious how to determine the sign of the vertical component of the slowness vector of the scattered waves. For elastic media, the sign is determined by applying so-called radiation condition when the slowness vector is complex-valued, but it has long been known that this approach does not work satisfactorily for anelastic media. Other approaches have been suggested, e.g., by requiring that the reflection and transmission coefficients should vary continuously with increasing incident angles, or by relating the sign to the direction of the energy flux. In the present paper, it is shown that these approaches may give different results, and that the results can be inconsistent with the elastic case even for weak attenuation. Instead, it is demonstrated that the ambiguity in the reflection coefficient can be resolved by expressing the seismic response of a point source over an interface as a superposition of plane waves and their reflection coefficients, and solving the resulting integral by the saddle point approximation. Although the saddle point itself (point of stationary phase) does not provide new insight, the ambiguity is removed by considering the steepest descent path through the point. Ray synthetic seismograms computed by this method compare well with synthetics computed by the reflectivity method, which does not suffer from the above-mentioned ambiguity since the integration path is taken along the real axis. This paper concentrates on the isotropic case, but it is discussed how the result may be extended to layered transversely isotropic media. The suggested approach, derived for a point source and plane layers, does not directly apply to 2-D or 3-D laterally inhomogeneous media, or to media of general anisotropy. A generalization of the result found is that the sign of the vertical slowness components should be chosen according to the energy flux direction for subcritical incidence and according to the radiation condition for supercritical incidence, even if this creates a discontinuity in the coefficients at the critical incidence angle. Such a discontinuity is sometimes necessary to get results which are consistent with the elastic case. It is discussed how the generalized result can be obtained by applying certain continuity criteria for the sub-and supercritical angle intervals, but the validity of this approach for general models remains to be proved.