Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

Probing mid-mantle heterogeneity using PKP coda waves

We investigate the utility of PKP coda waves for studying weak scattering from small-scale heterogeneity in the mid-mantle. Coda waves are potentially a useful probe of heterogeneity in the mid-mantle because they are not preferentially scattered near the CMB, as PKP precursors are, but are sensitive to scattering at all depths. PKP coda waves have not been used for this purpose historically because of interference with other late-arriving energy due to near-surface resonance and scattering. Any study of deep mantle scattering using coda waves requires the removal of near-surface effects from the data. We have analyzed 3624 recordings of PKP precursors and coda made by stations in the Incorporated Research Institutions for Seismology (IRIS) Global Seismographic Network (GSN). To study the range and time dependence of the scattered waves, we binned and stacked envelopes of the recordings. We have considered precursors that arrive within a 20 s window before PKP and coda waves in a 60 s window after PKP. The PKP scattered waves increase in amplitude rapidly with range as predicted by scattering theory. At ranges below 125°, we predict and observe essentially no scattered energy preceding PKP. Coda amplitudes at these ranges are independent of range and provide an estimate of energy due to near-surface effects that we can expect at all ranges. We use the average coda amplitude at ranges from 120 to 125° to correct coda amplitudes at other ranges. PKP coda waves show a strong dependence on time and range and are clearly influenced by scattering in the lower mantle. PKP coda waves, however, do not provide a tighter constraint on the vertical distribution of mantle heterogeneity than is provided by precursors. This is due, in part, to relatively large scatter in coda amplitudes as revealed by a resampling analysis. Modeling using Rayleigh–Born scattering theory and an exponential autocorrelation function shows that PKP coda amplitudes are not highly sensitive to the vertical distribution of heterogeneity in the mantle. To illustrate this we consider single-scattering in two extreme models of mantle heterogeneity. One allows heterogeneity just at the CMB; the other includes heterogeneity throughout the mantle. The amplitudes of precursors are tightly constrained by our stack and support our earlier conclusion that small-scale heterogeneity is uniformly distributed throughout the lower mantle. The best-fit model includes 8 km scale length heterogeneity with an rms velocity contrast throughout the mantle of 1%.     

沈阳满堂4.3级地震的烈度及地震地质背景

2003年3月30日在沈阳市郊满堂乡发生了ML4．3级地震，本文讨论了这次地震的烈度并对震区的地震地质条件作了初步分析。     

地质速率计—矿物氧扩散的应用

概述了同位素封闭体系内的矿物氧扩散和同位素交换机制及其在地质速率计上的应用。火成岩从高温冷却或变质岩从高峰主质温度冷却过程中，由于冷却速度不同，扩散作用导致的矿物晶体内部及晶粒间氧同位素再平衡也有所不同。通过实测岩石中各组成矿物氧同位素比值，模式含量和颗粒半径，据矿物氧扩散和同位素交换模型，可以估算出岩石的冷却速率。     

依兰-萝北地区地震危险区地震活动跟踪研究

本文选取地震活动图像演化和地震学参数变化两项指标，采用有时间段的时空扫描进行演化研究，逐步将这种演化研究推进到短期阶段进行跟踪研究。认为依兰-萝北地区近期发生5级左右地震可能性较小，黑龙江断裂两侧可能会有4级左右地震发生。     

Study of high temperature targets identification and temperature retrieval experimental model in SWIR remote sensing based Landsat8

For surface features in short-wave infrared (SWIR, 1.3–3.0 μm) in remote sensing imagery, pixel values depict the total energy including reflection and emission. For surface features at normal temperature in SWIR band, emission energy can be ignored. While for surface features at high temperature in SWIR band, emission energy is equal to or even higher than the reflection energy. So remote sensing imagery of SWIR band can be used to separate emission and reflection energy as well as to realize temperature retrieval of high temperature targets. In this study, the seventh band (SWIR band) of Landsat8 OLI remote sensing imagery is used to perform the theoretical model research for temperature retrieval of high temperature targets. In the meantime, it is also used with the corresponding observation experiment of synchronization satellite to check the theoretical model. The result shows that the radiant flux density for mixed pixels with high temperature targets is higher than adjacent pixels without high temperature targets. Thus, the high temperature pixels can be identified in SWIR band. The retrieval results of temperature and fractional area for high temperature targets are consistent with reality. In the study, the result illustrates that it is effective to identify high temperature targets in remote sensing imagery of SWIR band and the model is appropriate for temperature retrieval use.     

A review on advances in seismology in Mexico after 30 years from the 1985 earthquake

The 19 September 1985 (M_w8.1) earthquake, located on the Michoacán coast, Mexico, generated great damage in Mexico City, more than 300 km away from the epicentral area. Other important cities near the coast and in central Mexico also suffered severe damage. Thirty years after this important event, the National Autonomous University of Mexico (UNAM), the Centro de Investigación Científica y de Educación Superior de Ensenada, Baja California (CICESE) and other institutions organized a conference to discuss the scientific advances, particularly in seismology, that had taken place in Mexico since then.     

The North Atlantic Oscillation (NAO) and Europe Blocking and Their Impacts on Extreme Snowstorms: A Review

Large-scale, low frequency modes such as the North Atlantic Oscillation(NAO)and blocking, have an important modulation on the northern hemisphere weather and climate. In this paper, the physical mechanism studies on inter-decadal and decadal variability of NAO and blocking were summarized. The relationship between NAO regime transitions and the interannual variability of NAO in winter during was examined by using a statistical approach. The time-space relationship between NAO and European blocking were discussed. Based on two extreme cold and snowstorm events, the impacts on local weather especially the extreme events within the life cycle (two weeks) of the NAO and blocking were further examined. It was found that the frequently occurrence of the Eurasian extreme snowstorm was closely related to the special combination of NAO and blocking regime. In addition, the development of theoretical modes for NAO and blocking was discussed and issues that remain to be solved were proposed.     

火山地震的理论地震图计算方法

本文首先分析了不同类型火山地震的波形,提出了相应的震源模型.然后发展了计算层状介质理论地震图的部分分离变量—有限差分方法,其最重要的环节是引用了吸收边界条件,使计算工作得以简化.这种方法特别适用于火山地震的理论地震图计算.最后,本文给出了计算实例.