中国西部现今地壳形变与地震活动特征

以1999-2004年间中国西部地区近650个GPS站点的观测资料为基础,采用二维“高张力样条”函数内插算法获得了连续地壳形变场。结果表明,相对于稳定欧亚参考框架,中国西部现今地壳运动西强东弱,北向运动从西向东逐渐减弱,东向运动逐渐增强,青藏高原东缘及附近地区是东向运动的消减区带;主应变从西向东由压缩应变转变为拉伸应变;最大剪应变主要位于喜马拉雅地块及其东部地区;中国西部地区地震活动在空间分布上似乎集中发生在印度板块北东向挤压欧亚板块的两条共扼带内,与地壳形变场有一定对应关系,地壳形变的高应变率区为中强地震多发地带。     

基于非充分掺混模式的流域来水组成模型

针对如何减小数值求解对流输运方程耗散误差的问题,引进断面计算浓度的概念来计算河段平均浓度,提出了非充分掺混模式的有限控制体积法离散对流输运方程的新算法,据此构建了模拟流域内的水源组成以及不同水源在流域内的时空变化情况的流域来水组成模型。通过数值试验与具体实例验证,结果表明所提出的非充分掺混模式的新算法可有效地提高流域来水模型的计算精度。     

The Creation of a National Agricultural Land Use Dataset: Combining Pycnophylactic Interpolation with Dasymetric Mapping Techniques

Agricultural Census data is summarised over spatially coarse reporting units for reasons of farm confidentiality. This is problematic for research at a local level. This article describes an approach combining dasymetric and volume preserving techniques to create a national land use dataset at 1 km² resolution. The results for an English county are compared with contemporaneous aggregated habitat data. The results show that the accurate estimates of local agricultural land use (Arable and Grass) patterns can be estimated when individual 1 km squares are combined into blocks of > 9 squares, thereby providing local estimates of agricultural land use. This in turn allows more detailed modelling of land uses related to specific livestock and cropping activities. The dataset created by this work has been subject to extensive external validation through its incorporation into a number of other national models: nitrate leaching (e.g. MAGPIE, NEAP‐N), waste, and pathogen modelling related to agricultural activity.     

Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula

FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data. Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material. The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation strongly improves the simulation results and better represents steep fronts of final deposits.     

基于贡献权重迭加法的滑坡风险区划

贡献权重迭加法是对滑坡本底因子和承灾体因子在滑坡发育中的贡献率进行统计后,通过贡献率均值化、规一化处理,利用权重转换模型计算出每一个因子内部的权重——自权重W和因子相互之间的权重——互权重W＇。将滑坡因子、承灾体因子和两者的自权重、互权重分别相乘迭加,得到滑坡危险度和易损度区划结果,再将2项结果相乘得出滑坡风险区划。     

青藏高原西北缘祁漫塔格山中新世快速抬升 的磷灰石裂变径迹证据

对祁漫塔格山体不同海拔高度所取的9个磷灰石样品的裂变径迹分析结果表明,东昆仑西段中新世早中期为主要的隆升期且隆升速率较高,早期隆升速率为111m／Ma,晚期隆升速率为98m／Ma,总体隆升速率为100m／Ma。样品显示出磷灰石裂变径迹长度大致分2类,一类磷灰石裂变径迹长度为（12．21±10．20）-（13．75±0．30）μm,径迹长度分布图基本上为窄而对称的正态分布,反映具有快的剥露冷却速率,未受到后期热事件的干扰。另一类磷灰石裂变径迹长度为（11．88±0．33）～（13．32±0．27）μm,较前一类具有稍慢的剥露冷却速率,并且受到了后期热事件的干扰。     

努力提高三维地震勘探精度，为西部煤炭工业做出新贡献

以国家重大产业技术开发专项“西部煤炭资源高精度三维地震勘探技术”项目的由来、意义和总体研究目标为引,概括的介绍了项目依托工程中各个专项技术研究完成情况,并对非均匀介质成像技术、高精度三维地震静校正技术、高密度采集技术、特观技术、岩性反演技术、属性体解释技术等六项重大关键技术取得的突破性进展进行了重点说明。指出随着我国煤炭生产重点的逐步西移,应加强诸如叠前、叠后深度偏移技术的研究,以解决复杂山区三维地震面元内地震反射波散射问题,提高其三维地震勘探精度,为西部煤炭工业做出新贡献!     

三维地震小面元采集技术在晋城矿区的应用

依托“西部煤炭资源高精度三维地震勘探技术”工程,对晋城矿区进行了旨在提高小断层,小陷落柱探测能力的高密度三维地震勘探。根据面元选择因素及该区地质任务,采用5m×5m网格进行野外数据采集;考虑炮检距、方位角、覆盖次数、排列片横纵比及煤层埋深（350～500m）等因素,采用中点放炮、60道接收,24次覆盖（横向4次,纵向6次）的8线16炮束状观测系统,基岩中激发。原始资料经同一处理流程后,获得5m×5m×1ms、5m×10m×1ms、10m×10m×1ms及2.5m×2.5m×1ms不同单元的三维数据体多个,通过对比可以发现小断层,小陷落柱在其小面元叠加时间剖面、顺层切片及相干切片都有清晰的反映。实例说明,小面元采集技术可以提高对小构造的纵、横向分辨能力,满足山区对三维地震精确勘探的要求。     

Actinic flux and photolysis frequency comparison computations using the model PHOTOGT

The accurate radiative transfer model GOMETRAN, initially designed to yield radiances at TOA in the wavelength range 240–790 nm, has been extended to allow for the computation of actinic fluxes down to 175 nm and for the calculation of photolysis frequencies in the atmosphere. The capability of the extended model PHOTOGT (PHOTOGOMETRAN) is demonstrated in a number of successful comparison studies both with recent experiments (ground-based, balloonborne, airborne) and model calculations of radiances, actinic fluxes and photolysis frequencies in the stratosphere and troposphere. In an atmospheric case study, the impact of new quantum yield data for the O₃ » O₂+O(¹ d) photodissociation channel on the photolytic production of O(¹ d) atoms in the lower atmosphere has been quantified.