Three-dimensional geometrical ray theory and modelling of transmitted seismic energy of data from the Nevada Test Site

This paper presents a geometrically based algorithm for computing synthetic seismograms for energy transmitted through a 3-D velocity distribution. 3-D ray tracing is performed to compute the traveltimes and geometrical spreading (amplitude). The formulations of both kinematic and dynamic ray-tracing systems are presented. The two-point ray-tracing problem is solved by systematically updating the initial conditions and adjusting the ray direction until the ray intersects the specified endpoint. The amount of adjustment required depends on the derivatives of the position with respect to the given starting angles between consecutive rays. The algorithm uses derivatives to define the steepest-descent direction and to update the initial directions. The convergence rate depends on the complexity of the model.
Test seismograms compare favourably with those from a 2-D asymptotic ray theory algorithm and a 3-D Gaussian-beam algorithm. The algorithm is flexible in modelling arbitrary source and recorder geometries for various smoothly varying 3-D velocity distributions. The algorithm is further tested by simulating surface-to-tunnel vibroseis field data. Shear waves as well as compressional waves may be approximately included. Application of the algorithm to a data set from the Rainier Mesa of the Nevada Test Site produced a good fit to the transmitted (first arrival) traveltimes and amplitudes, with approximately 15 per cent variation in the local 3-D velocity.     

Oxygen isotope fractionation in mantle minerals

The increment method is adopted to calculate oxygen isotope fractionation factors for mantle minerals, particularly for the polymorphic phases of MgSiO₃ and Mg₂SiO₄. The results predict the following sequence of¹⁸O-enrichment:pyroxene (Mg, Fe, Ca)₂Si₂O₆>olivine (Mg, Fe)₂SiO₄ > spinel (Mg, Fe)₂SiO₄> ilmenite (Mg, Fe, Ca) SiO₃>perovskite (Mg, Fe, Ca) SiO₃. The calculated fractionations for the calcite-perovskite (CaTiO₃) System are in excellent agreement with the experimental calibrations. If there would be complete isotopic equilibration in the mantle, the spinel-structured silicates in the transition zone are predicted to be enriched in¹⁸O relative to the perovskite-structured silicates in the lower mantle but depleted in¹⁸O relative to olivines and pyroxenes in the upper mantle. The oxygen isotope layering of the mantle might result from differences in the chemical composition and crystal structure of mineral phases at different mantle depths. Assuming isotopic equilibrium on a whole earth scale, the chemical structure of the Earth’s interior can be described by the following sequence of¹⁸O-enrichment:upper crust>lower crust>upper mantle>transition zone>lower mantle>core. Project supported by the National Natural Science Foundation of China and the Chinese Academy of Sciences.     

北京天文台的WET三通道高速光度计

对北京天文台新近改造成的三通道高速光度计的结构和性能进行介绍,并给出了在兴隆站85公分望远镜的一些实测结果.     

平衡剖面研究中的新方法

随着现代化计算手段的运用的新的地学思维方法的影响,在平衡剖面研究中引入计算机模拟技术,分数维方法,三维模式等新方法,使平衡剖面技术具有了更加广阔的应用前景。另外,本文还对平衡剖面技术的局限性和发展趋势进行了探讨。     

Ground-based FTIR Measurements with High Temporal Resolution

We outline how ground-based Fourier transform infrared (FTIR) measurements of stratospheric trace species, obtained with high temporal resolution, could be used to detect filaments of polar vortex air at mid-latitudes and therefore test high spatial resolution chemical transport models (CTMs). Vertical column abundances of HCl, ClONO2, HNO3, N2O and HF have been obtained from FTIR solar absorption measurements made throughout the day from Aberdeen, UK (57°N, 2°W) on several days during winter/spring 1993/94 and 1994/95. The short-timescale (3366l7966464/xxlarge8764.gif" alt="sim" align="MIDDLE" BORDER="0"> 2 hours) variability observed in the columns is attributed to real atmospheric variations and is often associated with the passage of high latitude air over Aberdeen. This is confirmed by 3D modelling studies which qualitatively reproduce and rationalise the observed changes in the column data on January 19 1994, January 20 1995 and February 26 1995. We describe the viewing geometry of ground-based FTIR measurements and we suggest a measurement strategy which should maximise the information retrieved on horizontal gradients in stratospheric trace species columns from FTIR measurements.     

利用不完全Cholesky共轭梯度法求解点源三维地电场

点源三维地电场的求解是一大型数值计算问题．本文用有限差分方法求解,最后形成一个线性方程组Ax＝b,这里A是大型稀疏的带状对称矩阵．解大型稀疏方程组的完全Cholesky分解直接算法,一般要求巨大的机器内存来存储系数矩阵A,而且计算速度极慢．因此引入不完全Cholesky共轭梯度（ICCG）算法及按行索引的稀疏存储模式,充分利用A的稀疏性,使得计算速度大大提高,而内存要求则大大减少,因此ICCG算法是地电三维正演的强有力工具．     

Hemispheric simulation of the Asian summer monsoon

A three-level, 3w128l5nlp72r7/xxlarge946.gif" alt="beta" align="MIDDLE" BORDER="0">-plane, filtered model is used to simulate the Northern Hemisphere summer monsoon. A time-averaged initial state, devoid of sub-planetary scale waves, is integrated through 30 days on a 5° latitude-longitude grid. Day 25 through day 30 integrations are then repeated on a 2.5° grid. The planetary-scale waves are forced by time-independent, spatially varying diabatic heating. Energy is extracted via internal and surface frictional processes. Orography is excluded to simplify synoptic-scale energy sources.During integration the model energy first increases, but stabilizes near day 10. Subsequent flow patterns closely resemble the hemisphere summer monsoon. Climatological features remain quasi-stationary. At 200 mb high pressure dominates the land area, large-scale troughs are found over the Atlantic and Pacific Oceans, the easterly jet forms south of Asia, and subtropical jets develop in the westerlies. At 800 mb subtropical highs dominate the oceans and the monsoon trough develops over the Asian land mass. The planetary scales at all levels develop a realistic cellular structure from the passage of transient synoptic-scale features, e.g., a baroclinic cyclone track develops near 55°N and westward propagating waves form in the easterlies.Barotropic redistribution of kinetic energy is examined over a low-latitude zonal strip using a Fourier wave-space. In contrast to higher latitudes where the zonal flow and both longer and shorter waves are fed by barotropic energy redistribution from the baroclinically unstable wavelengths, the low-latitude waves have a planetary-scale kinetic energy source. Wave numbers 1 and 2 maintain both the zonal flow and all shorter scales via barotropic transfers. Transient and standing wave processes are examined individually and in combination.Wave energy accumulates at wave numbers 7 and 8 at 200 mb and at wave number 11 in the lower troposphere. The 800-mb waves are thermally indirect and in the mean they give energy to the zonal flow. These characteristics agree with atmospheric observation. The energy source for these waves is the three wave barotropic transfer. The implications of examining barotropic processes in a Fourier wave-space, vice the more common approach of separating the flow into a mean plus a deviation are discussed.     

关于黄土碳酸盐与气候环境关系探讨

黄土中碳酸盐与气候环境密切相关。从碳酸盐的成因看，干冷的气候利于原生碳酸盐的保存，而温湿的气候使大量的次生碳酸盐产生。碳酸盐的存在形式与降水量多少有关，降水量大可形成结核层和钙板等，而降水量小则形成各种胶膜。胶膜中各种不同形态的方解石胶膜又可反映不同的气候条件。黄土一古土壤中形成的碳酸盐实际上在古土壤期形成的碳酸盐要高于黄土期形成的，只是由于古土壤期强烈的风化淋溶作用，大量的碳酸盐从古土壤中淀积到下伏的黄土层中。黄土中CaCO3的淀积可分为连续淀积和不连续淀积两类。依据淀积深度和降水量的关系，可建立二者定量的关系。从黄土碳酸盐的不同方面，得到其与气候环境关系较为全面的认识。     

Chandra observations of neutron stars: an overview

We present an overview of Chandra X-ray Observatory observations of neutron stars. The outstanding spatial and spectral resolution of this great observatory have allowed for observations of unprecedented clarity and accuracy. Many of these observations have provided new insights into neutron star physics. We present an admittedly biased and overly brief review of these observations, highlighting some new discoveries made possible by the Observatory’s unique capabilities. This includes our analysis of recent multiwavelength observations of the putative pulsar and its pulsar-wind nebula in the IC443 SNR.