前兆异常的非线性识别与地震危险地点的判定

利用多种描述复杂现象的非线性方法处理了多台多项多年的前兆观测资料，通过“场”、“源”特征差异性的对比求得地震危险地点是一个能量“振荡”型积累的地区，具体表现为它是在前兆异常演化过程中出现的异常从属度最高、自相似性变化最大、前兆分维数下降最多的地区，而且还是在地震活动过程中出现的一个熵值相对较低的地区。研究结果还表明，震中位置的不确定性随着发震时刻的逼近而不断减弱     

Model prediction capacity of ephemeral gully evolution in conservation tillage systems

Ephemeral gully (EG) erosion has an important impact on agricultural soil losses and increases field surface hydrology connectivity and transport of pollutants to nearby water bodies. Watershed models including an EG component are scarce and not yet properly evaluated. The objective of this study is to evaluate the capacity of one such tool, AnnAGNPS, to simulate the evolution of two EG formed in a conservation tillage system. The dataset for model testing included runoff measurements and EG morphological characteristics during 3 years. Model evaluation focused on EG evolution of volume, width, and length model outputs, and included calibration and testing phases and a global sensitivity analysis (GSA). While the model did not fully reproduce width and length, the model efficiency to simulate EG volume was satisfactory for both calibration and testing phases, supporting the watershed management objectives of the model. GSA revealed that the most sensitive factors were EG depth, critical shear stress, headcut detachment exponent coefficient b, and headcut detachment leading coefficient a. For EG outputs the model was additive, showing low sensitivity to interactions between the inputs. Prediction of EG spatial evolution on conservation tillage systems requires improved development of gully erosion components, since many of the processes were developed originally for traditional tillage practices or larger channel systems. Our results identify the need for future research when EG form within conservation tillage systems, in particular to study gully headcut, soil erodibility, and width functions specific to these practices.     

单台地震自动定位网格搜索法及其MATLAB试验

提出了一种利用单个地震台站三分向数字地震波形记录进行单台地震自动定位的网格搜索法。该方法以台站为中心,建立辐射状网格,按可变震中距步长和方位角步长进行搜索。对一个检测到的地震事件记录波形,对每个震中距,计算预测波形特征函数与记录波形特征函数的互相关函数值,选其最大值对应的震中距;对每个方位角,合成沿径向和沿横向的波形分量,并计算P波段的径向绝对积分值与横向绝对积分值之差,选其最大值对应的方位。对于选定的方位,如果合成的径向分量与记录垂直分量乘积在P波段的积分值为负,则取该方位角为震中相对于台站的方位角,否则该角要加上180°。在MATLAB环境中对该方法的基本步骤进行了试验,初步结果表明该方法是可行的。     

Mathematical Analysis of a Model of River Channel Formation

The study of overland flow of water over an erodible sediment leads to a coupled model describing the evolution of the topographic elevation and the depth of the overland water film. The spatially uniform solution of this model is unstable, and this instability corresponds to the formation of rills, which in reality then grow and coalesce to form large-scale river channels. In this paper we consider the deduction and mathematical analysis of a deterministic model describing river channel formation and the evolution of its depth. The model involves a degenerate nonlinear parabolic equation (satisfied on the interior of the support of the solution) with a super-linear source term and a prescribed constant mass. We propose here a global formulation of the problem (formulated in the whole space, beyond the support of the solution) which allows us to show the existence of a solution and leads to a suitable numerical scheme for its approximation. A particular novelty of the model is that the evolving channel self-determines its own width, without the need to pose any extra conditions at the channel margin.     

中国数字地震台网（CDSN）数据文件系统：—CDSN应用系列论文一

中国数字地震台网(CDSN)观测并生成的数字地震资料是开展数字地震学研究必不可少的高质量数据源。为了使科学家们更好地使用这些数据,本文主要从数字地震资料产生出的各个环节、数据结构与格式、数据物理量之间的换算向读者作概括的介绍.希望读者对数字化地震资料的数据文件、数据记录、数据字均能有透彻的理解,从而节省时间和精力,提高资料的使用率。     

基于ANSYS软件benchmark模型耗能减振控制分析

本文系统介绍了当前国际结构振动控制的公共平台———Benchm ark模型,首次给出使用通用有限元软件ANSYS对第3阶段20层非线性钢框架抗震结构的有限元建模,结果表明,模态、无控地震反应分析结果和ASCE给出的MATLAB结果相同。利用粘滞阻尼器对此模型的耗能减振数值仿真分析的结果表明:粘滞阻尼器是一种性能很好的阻尼器,使用本文给出ANSYS的Benchm ark有限元模型来进行耗能减振控制分析是方便的、可行的。     

豫北地区视应力时空变化特征

利用2008年以来豫北地区49次M_L 2.5以上地震资料,通过Brune模型,结合遗传算法,反演计算得到豫北地区地震视应力范围为0.029 9 MPa—0.642 MPa,平均视应力为0.174 8 MPa,分析该区视应力时空变化特征,讨论震级和能量的关系,结果表明,豫北地区目前处于应力积累阶段,应力高值集中在聊兰断裂濮阳—鄄城段范县地区。     

Enlightenments from researches and predictions of 2014–2016 super El Nio event

<正>El Nio is a remarkable climate phenomenon with a basinwide warming of sea surface temperatures(SST) in the easterncentral tropical Pacific. El Nio means The Little Boy, or Christ Child in Spanish, and on the contrary, a basinwide cooling of the tropical Pacific SST is called La Nia that means The Little Girl in Spanish. Always, a large-scale SST change in the tropical     

Dissecting large earthquakes in Japan: Role of arc magma and fluids

We synthesized information from recent high-resolution tomographic studies of large crustal earthquakes which occurred in the Japanese Islands during 1995–2008. Prominent anomalies of low-velocity and high Poisson's ratio are revealed in the crust and uppermost mantle beneath the mainshock hypocenters, which may reflect arc magma and fluids that are produced by a combination of subducting slab dehydration and corner flow in the mantle wedge. Distribution of 164 crustal earthquakes ( M 5.7–8.0) that occurred in Japan during 1885–2008 also shows a correlation with the distribution of low-velocity zones in the crust and uppermost mantle. A qualitative model is proposed to explain the geophysical observations recorded so far in Japan. We consider that the nucleation of a large earthquake is not entirely a mechanical process, but is closely related to the subduction dynamics and physical and chemical properties of materials in the crust and upper mantle; in particular, the arc magma and fluids.     

Accelerating the discontinuous Galerkin method for seismic wave propagation simulations using multiple GPUs with CUDA and MPI

We have successfully ported an arbitrary high-order discontinuous Galerkin method for solving the three-dimensional isotropic elastic wave equation on unstructured tetrahedral meshes to multiple Graphic Processing Units (GPUs) using the Compute Unified Device Architecture (CUDA) of NVIDIA and Message Passing Interface (MPI) and obtained a speedup factor of about 28.3 for the single-precision version of our codes and a speedup factor of about 14.9 for the double-precision version. The GPU used in the comparisons is NVIDIA Tesla C2070 Fermi, and the CPU used is Intel Xeon W5660. To effectively overlap inter-process communication with computation, we separate the elements on each subdomain into inner and outer elements and complete the computation on outer elements and fill the MPI buffer first. While the MPI messages travel across the network, the GPU performs computation on inner elements, and all other calculations that do not use information of outer elements from neighboring subdomains. A significant portion of the speedup also comes from a customized matrix–matrix multiplication kernel, which is used extensively throughout our program. Preliminary performance analysis on our parallel GPU codes shows favorable strong and weak scalabilities.