地磁场能量在地球内部的分布及其长期变化

用国际参考地磁场模型(IGRF)分析了地磁场能量在地球内部的分布及其长期变化.结果表明,从1900年到2005年,地核以外地磁场总能量由6.818×1018J减少到6.594×1018J,减小了3.3%,地表以外地磁场总能量由8.658×101J减小到.63×101J,减小了11.4%.分析地球内部不同圈层地磁场能量的变化表明,地壳(A层)、上地幔(B层)、转换带(C层)、下地幔D′层的地磁场总能量在减小,但是下地幔"层的地磁场总能量却在快速增加.磁能密度随时间的变化更清楚地显示出磁能增加和减小的分界面在r=3840km处.上述结果表明,地核和地表以外地磁场总能量在趋势性减小的同时,也在进行重新分配.进一步分析表明,下地幔D"层磁能快速增长,主要是由高阶磁多极子的增强引起的.在地磁场倒转前,偶极矩减小而多极性相对增强在能量分布上的表现就是磁能向下地幔底部(特别是D"层)集中.     

非均质地基中平面应变隧道开挖面稳定上限分析

非均质是软黏土地基中比较普遍的现象,而目前隧道开挖面稳定研究中比较成熟的理论主要是针对均质土体。因此,从塑性极限分析上限法的基本原理出发,采用平面应变隧道刚体平动破坏模式(多块体上限法),考虑软黏土地基的非均质性,推导了平面应变隧道极限支护压力关于隧道埋深、土体重度及土体强度的上限公式。通过与其他方法的比较分析,证明了极限分析方法在隧道开挖面稳定性方面的可行性;利用该方法的计算结果详细探讨了隧道开挖面稳定的影响因素;而且由计算结果可知,地基土的非均质性在影响隧道开挖面极限支护压力的同时,也影响着隧道开挖破坏面的位置和形状,为工程实践提供重要的理论依据。     

初值及物理过程对\

利用国防科技大学全球中期数值天气预报模式（YinHe Global Spetral model,YHGS）产品驱动WRF对2018年7月4日华中地区暴雨过程进行模拟,并与ERA-interim资料作初始场模拟结果对比,评估YHGS模式产品在此次暴雨过程预报中的应用能力。结果表明：（1）WRF-YHGS对2018年7月4日华中地区暴雨过程有一定的预报能力,其模拟的大尺度环流形势、水汽收支量变化趋势与WRF-ERA有着很好的一致性,YHGS模式产品驱动中尺度数值预报是可行的。（2）WRF-YHGS模拟效果较WRF-ERA差,但大雨量级WRF-ERA湿偏差较大,两组试验各物理量模拟结果存在一定差距,且随着积分时间的增加差异逐渐增大。（3）WRF-YHGS、WRF-ERA模拟结果的差异主要来自YHGS与ERA初始场中差异较大的次天气尺度运动和YHGS全球模式预报场误差两个方面。

     

地电场对青藏高原东北部季节冻土的响应特性

地电场的变化与台址环境的水文、气象及地质背景等相关,在青藏高原东北部季节冻土区,11个地电场台站处于较高的海拔,据台址下覆场地属性分为A类（黄土型）和B类（高原草场型）台站。通过对青藏高原季节冻土区域的地电场和大地电流场的计算和分析,联系区域构造活动和地质环境得出以下认识：青藏高原东北部季节冻土区地电场变化对水热环境响应明显,冬、夏两季测值可能发生跃变;长周期的地电场变化曲线可能与台址附近气温变化相关;台站大地电流矢量在冻土部分冻融交替过程中发生方向和幅度值的改变。A类和B类台址显示出不同的季节变化规律,地电场曲线上升和下降的时间节点各异,这种现象可应用于监测该区域冻土冻融情况和冻土的时空演变。     

拉萨地磁台地磁相对记录室磁场梯度对比分析

测量拉萨地磁台观测区内自然地磁场梯度和建成地磁记录室后的地磁场梯度,并测量新建记录室内不同高度磁场梯度,对比分析新建地磁记录室内的磁场稳定性,进而找出潜在干扰源,并采取有效措施予以去除,为产出连续可靠的地磁相对记录数据提供环境条件.     

呼和浩特—包头盆地震源应力场状态分析

为了分析呼和浩特—包头盆地现今应力场状态,基于内蒙古地震台网记录的地震波形资料,运用振幅比方法获取2000—2016年该盆地及周边M_L≥2.7地震的震源机制解,计算其一致性参数,并反演区域应力场。结果表明,震源机制解特征主要表现为走滑特征,一致性参数空间分布处于低值状态,力轴方位基本与该区的应力场背景特征相符,即主压应力方位为北东、张应力方位为北西。     

An examination of the validation of a model of the hydro/thermo/mechanical behaviour of engineered clay barriers

This paper focuses attention on the development of a numerical model of the hydro/thermo/mechanical behaviour of unsaturated clay and its consequent verification and validation. The work presented describes on-going collaboration between the Cardiff School of Engineering and Atomic Energy of Canada. The model development, which was carried out at Cardiff, can be described as being based on a mechanistic approach to coupled heat, moisture and air flow. This is then linked to a deformation analysis of the material within a ‘consolidation’ type of model. The whole is solved via the finite element method to yield a computer software code named COMPASS (COde for Modelling PArtly Saturated Soil). Some aspects of verification and validation of the model have been addressed in-house. However, the purpose of current AECL work is to provide an independent, rigorous, structured programme of validation and the paper will also explore the further validation of COMPASS within this context. © 1998 by John Wiley & Sons, Ltd.     

Predicting reach‐specific properties of fluvial terraces to guide future fieldwork. A case study for the Late Quaternary River Allier (France) with the FLUVER2 model

Numerical models have not yet systematically been used to predict properties of fluvial terrace records in order to guide fieldwork and sampling. This paper explores the potential of the longitudinal profile model FLUVER2 to predict testable field properties of the relatively well‐studied, Late Quaternary Allier system in France. For the Allier terraces an overlapping ¹⁴C and U‐series chronology as well as a record of ¹⁰Be erosion rates exist. The FLUVER2 modelling exercise is focused on the last 50 ka of the upper Allier reach because for this location and period the constraints of the available dating techniques are tightest. A systematic calibration based on terrace occurrence and thicknesses was done using three internal parameters related to (1) the sediment erodibility; (2) the sediment transport distance; and (3) the sediment supply derived from the surrounding landscape. As external model inputs, the best available, reconstructed, tectonic, climatic and base‐level data were used. Calibrated model outputs demonstrate a plausible match with the existing fluvial record. Validation of model output was done by comparing the modelled and measured timing of aggradation and incision phases for the three locations. The modelled range of landscape erosion rates showed a reasonably good match with existing erosion rate estimates derived from ¹⁰Be measurements of fluvial sands. The quasi‐validated model simulation was subsequently used to make new testable predictions about the timing and location of aggradation and erosion phases for three locations along the Allier river. The validated simulations predict that along the Allier, reach‐specific dynamics of incision and aggradation, related to the variations in sediment supply by major tributaries, cause relevant differences in the local fluvial terrace stratigraphy. Copyright © 2016 John Wiley & Sons, Ltd.     

Stress Field Estimations for Colima Volcano, Mexico, Based on Seismic Data

 For first time, during 1991, seismic activity was recorded during an eruption at Colima volcano. We analyze these data to obtain a stress pattern using a composite focal mechanism technique. From the analysis of regional seismicity, the Tamazula Fault and the Armeria River appear as active features and the dip of the slab east of the Jalisco Block is approximately 12°. Southwest of Colima volcano a vertical alignment of seismic events was observed. We estimate five different composite focal mechanism solutions from our data set, which indicate a change of the stress field at the volcano after the 1991 eruption. These solutions suggest that the stress field in the volcanic edifice was controlled by stresses related to the emplacement of magma superimposed on the regional stress field. No evidence of active local faults in the volcanic edifice was found. We propose a model for the eruptive process that involves tilting of the volcanic edifice. Received: 15 October 1995 / Accepted: 26 October 1998     

叠加变形有限应变分解及其应用探讨

本文提出了二维不同类型叠加变形的有限应变分解方法。只要测定有限应变及能代表各期主应变面的面理或剪切指向,就可以计算出各期构造应变。