Scaling of response spectrum and duration for aftershocks

Aftershocks have the potential to cause collapse of a structure that has been already damaged by the preceding main shock. Seismic safety of a structure should therefore be ascertained through a damage analysis using the anticipated main shock and few larger-aftershock motions. Simulation of aftershock motions needs characterization of the seismic hazard due to aftershocks, and therefore it will be useful to develop a conditional scaling model that can predict the response spectrum of an anticipated aftershock motion consistent with the design spectrum of the main shock motion anticipated at the same station. In this study an attempt is made to develop a conditional scaling model for the pseudo spectral velocity spectrum via linear regression analysis on the aftershock and main shock recordings for the 1999 Chi–Chi earthquake. It is shown that it may be possible to obtain a simpler and approximate version of the conditional model from an unconditional model. Damage-causing potential of a ground motion also depends on its strong motion duration (SMD) and therefore a conditional scaling model is developed for SMD of the aftershock motion in several narrow frequency-bands. The model is developed for the larger-aftershock motions and it is shown that a reasonable replacement of such a model may be obtainable directly from an unconditional model. Finally, a simple weighted averaging scheme is proposed to obtain the composite SMD from the SMDs for different frequency bands by using the pseudo spectral acceleration spectrum of the motion.     

A hybrid discrete–finite element model for numerical simulation of geomaterials

A hybrid discrete–finite element model is introduced for simulation of mechanical behavior of geomaterials. The soil or rock is modeled as a system of discrete balls that interact through normal and shear springs. The balls can be bonded at the contact points to withstand the applied deviatoric stresses. The important feature of this model is that the confining walls that can be imagined for example as the surrounding membrane or the mold in a physical test are modeled by deformable finite elements. This allows simulation of laboratory test features more realistically compared to the situations where the surrounding walls are rigid. The relationships between micro- and macro-properties are investigated in this paper as well. These relationships and the corresponding curves are helpful tools in calibration of the numerical model for the macroscopic elastic properties.     

Crustal Growth by Magmatic Accretion Constrained by Metamorphic P-T Paths and Thermal Models of the Kohistan Arc, NW Himalayas

Magmatic accretion is potentially an important mechanism inthe growth of the continental crust and the formation of granulites.In this study, the thermal evolution of a magmatic arc in responseto magmatic accretion is modeled using numerical solutions ofthe one-dimensional heat conduction equation. The initial andboundary conditions used in the model are constrained by geologicalobservations made in the Kohistan area, NW Himalayas. Takingconsideration of the preferred intrusion locations for basalticmagmas, we consider two plausible modes of magmatic accretion:the first involves the repeated intrusion of basalt at mid-crustaldepths (‘intraplate model’), and the second evaluatesthe simultaneous intrusion of basalt and picrite at mid-crustaldepths and the base of the crust respectively (‘double-platemodel’). The results of the double-plate model accountfor both the inferred metamorphic P–T paths of the Kohistanmafic granulites and the continental geotherm determined frompeak P–T conditions observed for granulite terranes. Thedouble-plate model may be applicable as a key growth processfor the production of thick mafic lower crust in magmatic arcs. KEY WORDS: thermal model; magmatic underplating; P–T path; granulite; lower crust     

基于地质记录用时域组合模型预测气候变化趋势的初步研究

用时域分析组合模型建立了１００万年来６０°Ｎ６月份太阳辐射量时间序列、宝鸡黄土粒度时间序列、渭南夏季风指数时间序列的动态模型．研究结果表明,时域分析组合模型较好地提取了时间序列的信息,得到的几个显著周期Ｔ＝１３３,１００,８９,４１,２３,１９ｋａ,与地球轨道三要素的变化周期接近．用组合模型拟合实测数据,精度是高的;用其预测未来气候替代性指标时间序列的变化情况,发现未来气候有向干冷方向发展的趋势．时域分析组合模型为研究气候变化趋势提供了一种定量分析、预测的方法．     

基于准三维有限元方法建立的地震活动模型

基于有限元方法,参照细胞自动机模型,建立了一个地震活动演化的动力学模型.在这个模型中,根据边界条件和单元介质参数的分布,利用有限元方法确定各单元的应力增长速率.每个单元设定初始应力和摩擦因数,利用破裂准则可得到单元破裂所需时间.当单元破裂后,将破裂单元作为内部边界,利用有限元方法计算由于这个单元破裂对系统其他单元造成的应力调整.根据上述规则,建立一个由30×40个节点组成,在定常位移速率边界条件作用下的地震活动模型,探讨了区域应力场的动态演化、地震活动图像和各种参数变化对地震活动的影响.结果表明该模型在研究地震活动方面有一定优势并具有实用价值.     

雷暴探测研究的进展

基于达州市2015年10月-2016年9月的闪电定位和探空观测资料,以850 hPa与500 hPa温差、大气可降水量、K指数、对流有效位能、对流抑制和抬升指数作为雷暴预警因子,利用投影寻踪动态聚类方法对该时段内的雷暴个例构建了预警模型。结果表明：（1）模型预报结果定量评估的临界成功指数为72.00%,该模型对个例达到了识别和预警效果;（2）该预警模型与常规预警方法相比,具有识别率高,计算简便,客观性强等特点,可为雷暴预警提供了一种新的客观预报方法。

     

Compositional Zoning of the Bishop Tuff

Compositional data for >400 pumice clasts, organized accordingto eruptive sequence, crystal content, and texture, providenew perspectives on eruption and pre-eruptive evolution of the>600 km³ of zoned rhyolitic magma ejected as the Bishop Tuffduring formation of Long Valley caldera. Proportions and compositionsof different pumice types are given for each ignimbrite packageand for the intercalated plinian pumice-fall layers that eruptedsynchronously. Although withdrawal of the zoned magma was lesssystematic than previously realized, the overall sequence displaystrends toward greater proportions of less evolved pumice, morecrystals (0·5–24 wt %), and higher FeTi-oxide temperatures(714–818°C). No significant hiatus took place duringthe 6 day eruption of the Bishop Tuff, nearly all of which issuedfrom an integrated, zoned, unitary reservoir. Shortly beforeeruption, however, the zoned melt-dominant portion of the chamberwas invaded by batches of disparate lower-silica rhyolite magma,poorer in crystals than most of the resident magma but slightlyhotter and richer in Ba, Sr, and Ti. Interaction with residentmagma at the deepest levels tapped promoted growth of Ti-richrims on quartz, Ba-rich rims on sanidine, and entrapment ofnear-rim melt inclusions relatively enriched in Ba and CO₂.Varied amounts of mingling, even in higher parts of the chamber,led to the dark gray and swirly crystal-poor pumices sparselypresent in all ash-flow packages. As shown by FeTi-oxide geothermometry,the zoned rhyolitic chamber was hottest where crystal-richest,rendering any model of solidification fronts at the walls orroof unlikely. The main compositional gradient (75–195ppm Rb; 0·8–2·2 ppm Ta; 71–154 ppmZr; 0·40–1·73% FeO*) existed in the melt,prior to crystallization of the phenocryst suite observed, whichincluded zircon as much as 100 kyr older than the eruption.The compositions of crystals, though themselves largely unzoned,generally reflect magma temperature and the bulk compositionalgradient, implying both that few crystals settled or were transportedfar and that the observed crystals contributed little to establishingthat gradient. Upward increases in aqueous gas and dissolvedwater, combined with the adiabatic gradient (for the 5 km depthrange tapped) and the roofward decline in liquidus temperatureof the zoned melt, prevented significant crystallization againstthe roof, consistent with dominance of crystal-poor magma earlyin the eruption and lack of any roof-rind fragments among theBishop ejecta, before or after onset of caldera collapse. Amodel of secular incremental zoning is advanced wherein numerousbatches of crystal-poor melt were released from a mush zone(many kilometers thick) that floored the accumulating rhyoliticmelt-rich body. Each batch rose to its own appropriate levelin the melt-buoyancy gradient, which was self-sustaining againstwholesale convective re-homogenization, while the thick mushzone below buffered it against disruption by the deeper (non-rhyolitic)recharge that augmented the mush zone and thermally sustainedthe whole magma chamber. Crystal–melt fractionation wasthe dominant zoning process, but it took place not principallyin the shallow melt-rich body but mostly in the pluton-scalemush zone before and during batchwise melt extraction. KEY WORDS: Bishop Tuff; ignimbrite; magma zonation; mush model; rhyolite     

粗差验后方差的无偏估计与最优稳健估计

在正态粗差假设下导出了粗差验后方差的无偏估计,对误差工膨胀模型和误差均值移动模型,两者的无偏估计公式是相同的。这证明了李德仁验后方差的朱建军方差不是无偏的。由于偏方定义的彭方法是正态粗差假设下的最优稳健估计。     

蒙古塔木察格盆地塔南凹陷下白垩统层序结构类型、控制因素与层序发育模式

利用录井及地震反射资料,采用层序地层学原理对蒙古塔木察格盆地塔南凹陷下白垩统地层层序进行了划分和对比,划分出了4个三级层序（SQ1-SQ4）。在此基础上,结合盆地的构造演化特点,研究了层序结构特征,划分了层序类型,其中,SQ1和SQ2为初始裂谷期湖泊层序,SQ3为裂谷深陷期层序,SQ4为裂谷后过渡层序。层序的结构主要受构造沉降方式控制,塔南凹陷不同时期构造沉降方式不同,形成了3种层序发育模式,即简单箕状斜坡模式、差异沉降构造反转模式、坳陷型层序发育模式,不同的层序发育模式,其地层岩性圈闭的发育和分布不同。     

利用ARC/INFO进行三维地形模型的建立、处理和动态模拟

介绍了利用 ARC/ INFO建立三维地形模型的方法 ,以及利用其新版本的部分新功能处理三维地形模型、进行动态模拟的方法 ,与以往处理三维模型必须借助外部三维图形开发软件包的方法相比 ,具有高速、灵活、可视性强、支持强大分析功能等优点。形成的三维模型可随时与 ARC/ INFO的其它数据模型交换信息 ,与外部进程实时通讯 ,为使用者提供了二维分析工具所不能提供的信息 ,为开发视景仿真系统等与表面分析、立体显示有关的应用系统提供了强有力的支持。