Weak-contrast approximation of the elastic scattering matrix in anisotropic media

The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995.     

Evolution of cool skin and direct air-sea gas transfer coefficient during daytime

Absorption of solar radiation within the thermal molecular sublayer of the ocean can modify the temperature difference across the cool skin as well as the air-sea gas transfer. Our model of renewal type is based on the assumption that the thermal and diffusive molecular sublayers below the ocean surface undergo cyclic growth and destruction, the heat and gas transfer between the successive burst events are performed by molecular diffusion. The model has been upgraded to include heating due to solar radiation. The renewal time is parameterized as a function of the surface Richardson number and the Keulegan number. A Rayleigh number criterion characterizes the convective instability of the cool skin under solar heating. Under low wind speed conditions, the solar heating can damp the convective instability, strongly increasing the renewal time and correspondingly decreasing the interfacial gas exchange. In the ocean, an additional convective instability caused by salinity flux due to evaporation becomes of importance in such cases. The new parameterization is compared with the cool skin data obtained in the western equatorial Pacific during the Tropical Ocean Global Atmosphere Coupled Ocean Atmosphere Response Experiment in February 1993. In combination with a model of the diurnal thermocline it describes main features of the field data both in nighttime and daytime. Under low wind speed conditions (< 5 m s^-1) diurnal variations of the sea surface temperature due to the formation of a diurnal thermocline were substantially larger than those across the cool skin. Under wind speeds > 5 m s^-1, diurnal variations of the surface temperature due to the variations of the thermal molecular sublayer become more important.     

Modelling of the vertical fluxes of nitric acid,ammonia, and ammonium nitrate

Results from numerical investigations regarding the exchange of HNO₃, NH₃, and NH₄NO₃ between the atmosphere and the biosphere are presented. The investigations were performed with a modified inferential method which is based on the generally accepted micrometeorological ideas of the transfer of momentum, sensible heat and matter near the Earth's surface and the chemical reactions among these nitrogen compounds. This modified inferential method calculates the micrometeorological quantities (such as the friction velocity and the fluxes of sensible and latent heat), the height-invariant fluxes of the composed chemically conservative trace species with group concentrationsc ₁=[HNO₃]+[NH₄NO₃] (total nitrate),c ₂=[NH₃]+[NH₄NO₃] (total ammonia), andc ₃=[HNO₃]-[NH₃] as well as the fluxes of the individual nitrogen compounds. The parameterization of the fluxes is based on the flux-gradient relationships in the turbulent region of the atmospheric surface layer. The modified inferential method requires only the data of wind velocity, temperature, humidity and concentrations (HNO₃, NH₃, and NH₄NO₃) measured at a reference height by stations of a monitoring network.     

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

使用价值分析优化方法在功能—技术矩阵求解中的应用

功能—技术矩阵,经众多工作者的努力研究,已渐趋成熟与实用。但如何利用功能—技术矩阵构造在广义现代设计领域中适用的科学模式,解决多学科设计方案的组成与评价问题,还存在薄弱环节。本文利用“使用价值分析优化”的概念,从设计方案技术先进性和经济合理性的协调统一的角度出发,较好地实现了功能—技术矩阵可行优化组合链的组成与评价,弥补了以往的缺陷。     

Energy transfer among Pb, Ce and Mn in fluorescent calcite from Kuerle, Xinjiang, China

Natural calcite from Kuerle, Xinjiang, China, shows orange-red fluorescence when exposed to short-wave ultraviolet (UV) light (Hg 253.7 nm). Photoluminescence (PL) emission and excitation spectra of the calcite are observed at room temperature in detail. The PL emission spectrum under 208 nm excitation consists of three bands: two UV bands at 325 and 355 nm and an orange-red band at 620 nm. The three bands are ascribed to Pb²⁺, Ce³⁺ and Mn²⁺, respectively, as activators. The Pb²⁺ excitation band is observed at 243 nm, and the Ce³⁺ excitation band at 295 nm. The Pb²⁺ excitation band is also observed by monitoring the Ce³⁺ fluorescence, and the Pb²⁺ and Ce³⁺ excitation bands, in addition to six Mn²⁺ excitation bands, are also observed by monitoring the Mn²⁺ fluorescence. These indicate that four types of the energy transfer can occur in calcite through the following processes: (1) Pb²⁺ → Ce³⁺, (2) Pb²⁺ → Mn²⁺, (3) Ce³⁺ → Mn²⁺ and (4) Pb²⁺ → Ce³⁺ → Mn²⁺.     

模糊理论在公路边坡稳定性分析中的应用

边坡工程是一个复杂的系统工程。边坡的稳定与否直接关系到其防护工作的程度。因此,如何有效、合理地判断边坡的稳定性对整个工程的经济与社会效益极为重要。影响边坡稳定的因素众多,判断其稳定与否应考虑各因素的综合效果。论文运用模糊理论中的多因素综合评判模型对某高速公路边坡的稳定性进行研究。通过全面考虑影响因素,建立相应的评判因素集及单因素评判矩阵,运用专家评判法和判断矩阵法确定出各因素的权重。最后通过多层次模糊评判得到该边坡的稳定性评判结果。研究证明,该法在边坡稳定性分析评判中是切实可行的。     

引黄水在北大港水库蓄存期水质咸化的研究

依据引黄水在北大港水库蓄存期水质咸化的事实以及传质理论,定性地分析了引黄水水质咸化的主要原因,包括盐碱土传质作用、蒸发浓缩作用以及风的扰动加速传质作用等.在北大港水库旁修建了一个模拟试验水池,进行了现场模拟试验研究.依据水量平衡和氯离子总量平衡的原理,分别定量地计算了传质作用和蒸发浓缩作用对水质咸化影响的比重.结果表明,在模拟水池蓄水后的30 d内,传质作用对水质咸化影响的比重有随时间增加而增长的趋势;30 d后,传质作用对水质咸化影响的比重稳定在80%,蒸发浓缩作用对水质咸化影响的比重稳定在20%,说明导致水质咸化的最主要因素是盐碱土传质作用.     

粗差发现和定位能力与相关系数的关系

根据粗差判断方程中的判断矩阵和两个统计检验量之间相关系数的函数式,论证了两种不同的研究方法所确定的粗差不能定位的数学模型实际上是相等的。通过算例,不仅说明两种研究方法对观测量不能定位粗差的判断是一致的,而且使用判断矩阵研究观测量的粗差发现和定位能力会更加方便简单。     

不同辐射传输方案对中尺度降水影响的对比分析

在MM5非静力稳定中尺度气象模式中引进了建立在δ-4流近似和相关-k分布基础上的对云水、雨水、冰晶和霰的辐射特性进行详细描述的辐射传输方案。新建立的辐射传输方案和MM5中原有的辐射传输方案在华南暴雨中的模拟结果相互比较，并与天气实况的对比表明:辐射在中尺度暴雨中起着重要的作用；辐射传输方案对云辐射特性描述的准确程度对于地面降水影响是明显的；不同的辐射传输方案对地面降水的影响存在较大的差异，并且这些差异在白天比在夜间明显；辐射传输过程对地面降水影响的差异主要表现在降水中心上，而对降水的地理分布改变很小；相对而言，不同的辐射传输方案之间对短波描述的差异对地面降水的影响明显，而对长波描述差异的影响不大；新辐射传输方案能够在一定程度上改进MM5对中尺度降水的模拟能力。