Characterisation of the 1997 Vulcanian explosions of Soufrière Hills Volcano,Montserrat, by video analysis

The activity of Convention at Montserrat Soufrière Hills Volcano, Montserrat, during the period 1995–1999 included numerous violent explosions. Two major cycles of Vulcanian explosions occurred in 1997: a first of 13 explosions between 4 and 12 August and a second of 75 between 22 September and 21 October. The explosions were short-lived events lasting a few tens of seconds during which partial fountain collapse generated pyroclastic surges and pyroclastic flows, and buoyant plumes ascended 3–15 km into the atmosphere. Each explosion discharged on average 3×10⁵ m³ (dense-rock equivalent, DRE) of magma, draining the conduit to depths of 1–2 km. The paper focuses on the first few seconds of three explosions of the 75 that occurred in September/October 1997: 6 October 1997 at 17:50, 7 October 1997 at 16:02 and 9 October 1997 at 12:32. Physical parameters such as exit velocities, magmatic water contents and magma pressures at fragmentation are estimated by following and modelling the ascent of individual momentum-dominated finger jets visible on videos during the initial stages of each explosion. The model treats each finger jet as an incompressible flow sustained by a steady flux of gas and particles during the few seconds of ascent, and produces results that compare favourably with those using a multiphase compressible code run using similar eruptive parameters. Each explosion reveals a progressive increase in eruptive intensity with time, jet exit velocities increasing from 40 m s^–1 at the beginning of the explosion up to 140 m s^–1 after a few seconds. Modelling suggests that the first magma to exit was largely degassed, whereas that discharged after a few seconds contained up to 2 wt% water. Magma overpressures up to ~10 MPa are estimated to have existed in the conduit immediately prior to each explosion. Progressive increases in jet exit velocity with time over the first few seconds of each explosion provide direct evidence for strong pre-eruptive gradients in water content and magma pressure in the upper reaches (probably 100–500 m) of the conduit. Fountain collapse occurred during the first 10–20 s of each explosion because the discharging jets had bulk densities up to 100 times that of the atmosphere and were unable to entrain enough air to become buoyant. Such high eruptive densities were due to the presence of partially degassed magma in the conduit.Editorial responsibility: A. Woods     

Study on elastic response of structures to near-fault ground motions through record decomposition

Accelerograms recorded near active faults have some important characteristics that make them different from those recorded in far-fault regions. High-frequency components in acceleration records and long-period velocity pulses are among notable specifications of such ground motions. In this paper, a moving average filtering with appropriate cut-off frequency has been used to decompose the near-fault ground motions into two components having different frequency contents: first, Pulse-Type Record (PTR) that possesses long-period pulses; second, the relatively high-frequency BackGround Record (BGR), which does not include large velocity pulses. Comparing the results with those extracted through wavelet analysis shows that moving average filter is an appropriate and efficient tool for near-fault records decomposition. The method is applied to decompose a suite of 91 selected near-fault records and the elastic response of structures is examined through their response to the decomposed parts. The results emphasizes that in contrast with ordinary far-fault earthquake records, response spectra of near-fault ground motions typically have two distinct local peaks, which are representatives of the high- and low-frequency components, i.e., BGR and PTR, respectively. Moreover, a threshold period is identified below which the response of structures is dominated by BGR while PTR controls the response of structures with periods longer than this period.     

大兴安岭多年冻土泥炭地无机氮动态对秋季冻融的响应北大核心CSCD

大兴安岭多年冻土泥炭地是对全球变暖响应敏感的地区之一。在全球变暖、多年冻土退化背景下,为了探明秋季冻融对多年冻土泥炭地无机氮时空变化的影响,本研究于2019年9—11月以大兴安岭三种多年冻土泥炭地为研究对象进行野外原位实验,分析了秋季冻融前、中和后期多年冻土泥炭地浅层和深层土壤无机氮的时空变化特征以及浅层和深层土壤含水量和温度的变化规律,建立了土壤无机氮含量与土壤温度和含水量间的多元线性回归模型。研究表明:多年冻土小叶章泥炭地(XY)、兴安落叶松-泥炭藓泥炭地(XA)和白毛羊胡子苔草泥炭地(BM)的土壤铵态氮(NH_(4)^(+)-N)含量变化范围:(1.00±0.00)~(20.60±0.20)mg·kg^(-1),硝态氮(NO_(3)^(-)-N)含量的变化范围:(0.02±0.01)~(14.64±1.11)mg·kg^(-1),且无机氮以土壤NH_(4)^(+)-N为主;秋季冻融后期无机氮含量明显高于前期。尽管水热交互作用对该时期无机氮没有显著影响,但是在不同冻融阶段,无机氮对环境因子的响应程度存在差异:在秋季冻融前、中和后期浅层无机氮动态分别与浅层温度和含水量的变化相关,但在整个秋季冻融期间BM浅层无机氮含量仅对10~20 cm含水量存在响应(R^(2)=0.344,P<0.01)。研究表明,秋季冻融期内,多年冻土泥炭地无机氮发生初步累积,且浅层环境因子对无机氮响应程度最大。本研究可补充大兴安岭多年冻土泥炭地秋季冻融对土壤无机氮影响研究的相关数据,并为多年冻土泥炭地响应全球变暖的温室气体释放的研究提供基础数据支撑。     

冻融循环作用下含水率对粗颗粒填料水分迁移影响的宏细观试验研究北大核心CSCD

冻融循环作用和初始含水率是影响粗颗粒填料水分迁移特征和冻胀融沉变形的两个主要因素。为明确冻融循环作用下,不同初始含水率粗颗粒填料的水分迁移特征及细观机制,采用荧光素为追踪剂,以CT细观机理观测为研究手段,开展了不同初始含水率条件下粗颗粒填料的一系列冻融循环试验,探究了温度场和冻深变化、水分迁移图像、补水量变化、最终含水率分布及CT值变化规律等。试验结果表明,冻深基本随冻融循环的次数增加不断加深,初始含水率越大,不同冻融循环作用下冻深的变化越稳定。外界补水量和液态水迁移高度与初始含水率呈负相关。CT扫描结果分析表明,经历多次冻融循环后,试样中的水分迁移导致土体孔隙结构及颗粒构造发生变化,试样的密度普遍增加,孔隙率以减小为主,进而导致土体发生相应变形。     

考虑含水率影响的结构性黄土临界状态模型

根据不同初始含水率原状黄土结构性演化规律,基于临界状态土力学理论,以应力、初始含水率和应变为基本变量,提出了一个可以反映原状黄土结构性演化规律及软化特性的临界状态本构模型。模型通过对比不同初始含水率下重塑与原状黄土等向压缩曲线,建立了不同含水率下原状黄土的结构性参数及其演化方程。此外,模型采用非相关联流动法则,以剪胀方程的形式来求解塑性偏应变。模型共计9个材料参数,均可由压缩试验和常规三轴剪切试验求得。通过与已有试验数据的对比发现,本模型不仅可以较好地体现初始含水率对原状黄土结构强度、变形特性以及结构破坏规律的影响,而且能较为合理地预测原状黄土的硬化及软化特性。所建立的结构性黄土临界状态模型为深化黄土力学特性研究提供了可能,同时为有效计算黄土地基湿陷变形提供了一定的理论依据。     

Differential Thermal Regimes of the Tarim and Sichuan Basins in China: Implications for Hydrocarbon Generation and Conservation

The uncertainty surrounding the thermal regimes of the ultra-deep strata in the Tarim and Sichuan basins, China, is unfavorable for further hydrocarbon exploration. This study summarizes and contrasts the present-day and paleo heat flow, geothermal gradient and deep formation temperatures of the Tarim and Sichuan basins. The average heat flow of the Tarim and Sichuan basins are 42.5 ± 7.6 mW/m2 and 53.8 ± 7.6 mW/m2, respectively, reflecting the characteristics of ‘cold’ and ‘warm’ basins. The geothermal gradient with unified depths of 0–5,000 m, 0–6,000 m and 0–7,000 m in the Tarim Basin are 21.6 ± 2.9 °C/km, 20.5 ± 2.8 °C/km and 19.6 ± 2.8 °C/km, respectively, while the geothermal gradient with unified depths of 0–5,000 m, 0–6,000m and 0–7,000 m in the Sichuan Basin are 21.9 ± 2.3 °C/km, 22.1 ± 2.5 °C/km and 23.3 ± 2.4 °C/km, respectively. The differential change of the geothermal gradient between the Tarim and Sichuan basins with depth probably results from the rock thermal conductivity and heat production rate. The formation temperatures at depths of 6,000 m, 7,000 m, 8,000 m, 9,000 m and 10,000 m in the Tarim Basin are 80°C–190°C, 90°C–220°C, 100°C–230°C, 110°C–240°C and 120°C–250°C, respectively, while the formation temperatures at depths of 6,000 m, 7,000 m, 8,000 m and 9,000 m in the Sichuan Basin are 120°C–200°C, 140°C–210°C, 160°C–260°C and 180°C–280°C, respectively. The horizontal distribution pattern of the ultra-deep formation temperatures in the Tarim and Sichuan basins is mainly affected by the basement relief, fault activity and hydrothermal upwelling. The thermal modeling revealed that the paleo-heat flow in the interior of the Tarim Basin decreased since the early Cambrian with an early Permian abrupt peak, while that in the Sichuan Basin experienced three stages of steady state from Cambrian to early Permian, rapidly rising at the end of the early Permian and declining since the late Permian. The thermal regime of the Sichuan Basin was always higher than that of the Tarim Basin, which results in differential oil and gas generation and conservation in the ultra-deep ancient strata. This study not only promotes theoretical development in the exploration of ultra-deep geothermal fields, but also plays an important role in determining the maturation phase of the ultra-deep source rocks and the occurrence state of hydrocarbons in the Tarim and Sichuan basins.     

利用气象地温资料反演大地热流

利用气象台记录的深度为0.8m,1.6m和3.2m的月平均地温资料,进行Fourier频谱分析,得到定常波和子谐波（变化波）两部分.同时取地表以下为水平分层的均匀双层介质,建立数学模型,根据线性热传导和最优化理论对定常波以及子谐波的振幅谱和相位谱进行联合反演,计算气象台地表附近的定常地温梯度以及地表土壤的热扩散系数等物性参数,并据此计算该地的大地热流值.初步计算结果表明,上述方法得到的某些地区的大地热流值与传统的通过钻孔测量法得到的大地热流值符合得相当好.     

基于CFD技术的河谷型城市风环境模拟——以兰州市城关区为例

本文应用CFD软件对多种流场的强大模拟功能,模拟小尺度的兰州市城关区风环境。选用SRTM资料表征兰州市城关区的地形,应用ArcGIS生成兰州市地形模型,计算相关参数,将计算得到的地形粗糙度和兰州多年气候资料转化为CFD模型的数字参数及模拟边界条件等,首次模拟得到河谷型地区兰州城关区的风环境。结果表明:（1）河谷型地区中,气流过山时风速达最大值,依山体高度增大,而在背风坡风速值明显降低,达到0～1 m·s-1;（2）气流过境时在沿河谷前进,且在谷底风速值达0～1 m·s-1;（3）在模拟城市热岛效应的对比实验中,地块加热对气流过境有停滞作用。     

荒漠河岸胡杨(Populus euphratica)树干液流的时滞效应

以2012年4—9月树干液流速率与环境因子的实测数据,分析了荒漠河岸胡杨(Populus euphratica)树干液流的时滞效应。结果表明:不同月份胡杨树干液流速率对不同环境因子的时滞不同;4—9月,胡杨树干液流速率峰值时间比太阳辐射峰值时间晚约1h,比空气温度峰值时间早约1h,比水汽压差峰值时间早约2h,比空气相对湿度谷值时间早约2h;对影响胡杨树干液流速率的环境因子进行了主成分分析,发现胡杨树干液流速率与1h以前的第1主成分、第2主成分相关性最高,相关系数分别为0.864、0.875;对错位前后胡杨树干液流速率与环境因子进行多元线性回归分析发现,在胡杨树干液流的数值模拟中,考虑液流相对于环境因子的时滞效应可以提高模型的拟合精度。     

基于遥感影像的哈尔滨市城市热环境变化

以Landsat TM/ETM+为基本数据源,利用遥感和GIS技术定量反演了哈尔滨市1989和2000年的地面亮温和植被指数(NDVI),对标准化地面亮温分级数据进行热力景观指数计算,研究城市热场分布的空间格局、时空变化规律和不同植被指数(NDVI)等级下热环境的变化特征.结果表明:1989和2000年哈尔滨市建成区均存在显著的热岛效应,2000年建成区热岛效应增强,哈尔滨市总体地面亮温升高;热力景观日趋破碎,各景观类型受干扰程度较大,处于很不稳定的状态.NDVI与地面亮温的负相关性随着植被覆盖度的升高而增大;植被覆盖度与地面亮温的空间变异呈负相关.