2000年苏门答腊M_S8.0地震对缅甸M_S6.5地震动态应力触发作用及影响因素

根据地震动态触发理论,计算了2000年苏门答腊MS8.0,地震在缅甸MS6.5地震破裂面上产生的动态库仑应力:正峰值为1.357 5 MPa,负峰值为-1.113 5 MPa。峰值大于地震触发阈值0.01 MPa,苏门答腊地震的发生对缅甸地震具有一定的触发作用。在计算过程中发现,主震基本参数、传播过程中介质参数和被触发地震断层参数的准确性对计算结果有较大影响。     

Devonian integrative stratigraphy and timescale of China

The Global Boundary Stratotype Sections and Points(GSSPs) for the bases of all seven international Devonian stages have been formally defined and ratified by IUGS till 1996, and nowadays, the main tasks for Devonian stratigraphers include further subdivision of these standard stages, strictly constrained absolute ages for the boundaries, and precise neritic-pelagic and marine-terrestrial correlations using multidisciplinary stratigraphy methods. Establishment of high-resolution Devonian integrative stratigraphy framework and timescale of China would play an important role in improving regional and international correlation, facilitating the recognition of important stratigraphic levels in different paleogeographic settings, and understanding the evolution pattern of biota, paleoclimate and paleoenvironment during this critical interval. Based on well-studied bio-and chronostratigraphy of Devonian in South China and adjacent areas, in combination with recent achievements in carbon isotope stratigraphy, event stratigraphy and radioactive isotope ages, this paper briefly summarize the research history and current status of Devonian chronostratigraphy of China, and for the first time introduce Devonian integrative stratigraphy framework of China.Up to date, few studies have been conducted on the astronomical cyclostratigraphy and high-resolution radioactive isotope dating in Devonian of China, which should be our main focuses in the near future.     

Time-averaged sub-layer concentration profiles in oscillatory sheet flows

Although thin on the order of several to tens of millimeters, sheet flows normally comprise a lower pick-up sub-layer and an upper contact-load sub-layer, separated at the bed level (z = 0). The time-averaged concentration profile in the pick-up sub-layer shows a ‘convex upward’ curvature, but ‘concave upward’ shape characterizes the time-averaged concentration profile in the contact-load sub-layer. The time-dependent concentration in the contact-load sub-layer is approximately in-phase with the free stream flow velocity, whereas it is nearly in anti-phase with the free stream flow velocity in the pick-up sub-layer. Two distinct analytical expressions of the time-averaged concentration profiles for the respective sub-layers are proposed. The expressions are validated with detailed observation datasets collected in the Groβer Wellenkanal (GWK) prototype wave flume in Hannover, Germany. The agreement between the predicted and the measured values is excellent. Interparticle collisions in the pick-up sub-layer and convective lifting processes associated with vortex shedding in the contact-load sub-layer are considered responsible for the opposite curvatures and in-phase/anti-phase concentration variations. Both transitional boundary and reference concentration are also elaborated.     

Integrated Two-Dimensional Modeling of Fluid Flow and Compaction in a Sedimentary Basin

In this paper we employ mixed finite elements and numerically study an integrated two-dimensional model of fluid flow and compaction in a sedimentary basin. This model describes a single phase incompressible flow in a two-dimensional section of a sedimentary basin with vertical compaction. At each time step, an iterative algorithm is used to solve this model. The determination of the grid movement is based on the mass conservation and movement of sediments in the basin, while the mixed method is utilized to solve the fluid flow over the moving grid. Numerical experiments are presented to verify this iterative algorithm and show representative solutions for the model under consideration.     

Microseism Induced by Transient Release of In Situ Stress During Deep Rock Mass Excavation by Blasting

During deep rock mass excavation with the method of drill and blast, accompanying the secession of rock fragments and the formation of a new free surface, in situ stress on this boundary is suddenly released within several milliseconds, which is termed the transient release of in situ stress. In this process, enormous strain energy around the excavation face is instantly released in the form of kinetic energy and it inevitably induces microseismic events in surrounding rock masses. Thus, blasting excavation-induced microseismic vibrations in high-stress rock masses are attributed to the combined action of explosion and the transient release of in situ stress. The intensity of stress release-induced microseisms, which depends mainly on the magnitude of the in situ stress and the dimension of the excavation face, is comparable to that of explosion-induced vibrations. With the methods of time–energy density analysis, amplitude spectrum analysis, and finite impulse response (FIR) digital filter, microseismic vibrations induced by the transient release of in situ stress were identified and separated from recorded microseismic signals during a blast of deep rock masses in the Pubugou Hydropower Station. The results show that the low-frequency component in the microseismic records results mainly from the transient release of in situ stress, while the high-frequency component originates primarily from explosion. In addition, a numerical simulation was conducted to demonstrate the occurrence of microseismic events by the transient release of in situ stress, and the results seem to have confirmed fairly well the separated vibrations from microseismic records.     

The Paleo－Environmental Significance of Sr and Ca Contained in Acid－Soluble Fraction of Plateau Loess Deposits in Xining Basin

The Xining Basin,located on northeastern Qinghai-Xizang Plateau,is adjacent to the Qaidam Bsim and Gobi-desert.Therefore,thick,integrated deposits of loess or paleo-sol sequence were well developed.As is well known,there is much paleo-environmental information recorded in loess and plaeo-sol.So if we deal with paleo-environmental index from different an gles ,especially from loess geochemistry,some important and detailed paleo-environmental evolutional information will be acquired.Based on this,the authors determined the contents of Sr and Ca contained in acid-soluble fraction from the bottom to the top of the Panzishan section.This section was obtained through making a well about 17 m in depth,including S1 paleo-sol that was developed during the last interglacial period.According to the results,it is found that variations in contents of Sr and Ca are well correlated with conversion of loess and paleo-sol.In loess deposits,the contents of Sr and Ca are higher than those in paleo-sol deposits.It is known that loess was deposited under cold-dry climate conditions and paleo-sol was developed during warm-humid climate episodes.It is concluded that due to their mobile behaviors,under cold-dry climate conditions,Sr and Ca are stagnant during weak weathering,but under warm-humid climate conditions,Sr and Ca are easily leached and transported in aqueous solution.Therefore,they can be utilized as indicatiors of climate changes.Moreover,variations in contents of Sr and Ca are more sensitive than susceptibility.The mobile behaviors of Sr and Ca contained in acid-soluble fraction can serve as a new quantitative index to evaluate paleo-environmental changes on the plateau.     

局部区域模型重力异常快速算法

地球重力场位系数模型可以用于计算局部重力扰动场元。然而随着地球重力场模型阶次的提高、局域重力场计算范围的增大,其计算速度往往不能满足工程需求。针对这一问题,在对位系数模型泰勒级数展开的基础上提出了采用向量运算、混合编程的方法,同时对连带勒让德函数Belikov递推方法中与经纬度无关的量进行了预先计算,有效提高了计算速度。提出的方法对于利用超高阶次重力场模型快速解算大范围、高分辨率重力场元数据以及累加求和计算具有一定的参考与借鉴意义。     

A discrete global grid system for earth system modeling

To support Earth system modeling, we propose a discrete global grid system that expresses multi-resolution spatial data. Specifically, a unified coding model that expresses a grid of nodes, edges, and cells is constructed for a triangular discrete global grid system. To fulfill the requirements of practical applications, we design a code-based topological query method for this grid system and an algorithm to transform between grid codes and geographic coordinates. We evaluate the Global Finite Volume Community Ocean Model (Global-FVCOM) on the triangular discrete global grid system in the proposed uniform coding model. The ocean tidal waves simulated by the Global-FVCOM running on the coded grid are then compared with results obtained using a traditional irregular spherical grid system, and the results display comparable accuracy. The uniform coding model proposed in this paper provides a triangular discrete global grid system that can represent multi-resolution spatial data and can be used in Earth system models. This unified coding model can also be applied to the geographic coordinate system made up of latitudes and longitudes, as well as diamond and hexagonal grids.     

Gamma-ray bursts: post-burst evolution of fireballs

The post-burst evolution of fireballs that produce γ-ray bursts (GRBs) is studied, assuming the expansion of fireballs to be adiabatic and relativistic. Numerical results as well as an approximate analytic solution for the evolution are presented. Owing to the adoption of a new relation between t R and γ, our results differ markedly from previous studies. Synchrotron radiation from the shocked interstellar medium is carefully calculated, using a conventional set of equations. The observed X-ray flux of GRB afterglows can be reproduced easily. Although the optical afterglows seem much more complicated, our results can still present a rather satisfactory agreement with observations. We also find that the expansion will no longer be highly relativistic about 4 d after the main GRB. We thus suggest that the marginally relativistic phase of the expansion should be investigated so as to check the afterglows observed a week or more later.     

The tunnel seismic advance prediction method with wide illumination and a high signal-to-noise ratio

Tunnel seismic prediction is widely used in the field of tunnel seismic advance detection. The illumination of the target and the signal-to-noise ratio of the data are two key factors affecting the precision of data interpretation. Current seismic prospecting has shortcomings on sites: (1) The lighting shots are solely towards one side of the tunnel wall, (2) the geophones are placed far away from the tunnel face and (3) the surface waves from the tunnel wall dominate over the reflection waves, lowering the signal-to-noise ratio of the data at the tunnel wall. This paper proposes a tunnel symmetrical geometry to tackle the above challenges. The arrangement is to place 12 sources uniformly on each side of the tunnel wall and six geophones on the tunnel wall and face. Results of simulated data and measured data show that the proposed method enables (1) broad illumination of the target body, (2) the enhancement of illumination energy of the target body, and (3) higher data signal-to-noise ratio. The proposed symmetrical geometry method provides better interpretation in terms of broader coverage, higher quality and greater distance of investigation.