Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

Estimation of annual actual evapotranspiration from nonsaturated land surfaces with conventional meteorological data

Land surface evapotranspiration is an important component both in earth surface heat and water bal-ance, on whose budgets weather and climate depend, to a great extent, for their changes are responsible for the formation and variation of vegetation features on the globe. Besides, the evapotranspiration is an im-portant topic of short-term flood forecasting and the estimation of runoff from mountainous sides. As a result, the problem as to the evapotranspiration has been one of the concerns in …     

山西省主要河流流域面雨量预报业务流程

以T213、HLAFS模式、MM5中尺度模式输出的格点资料以及日本降水量格点资料为基础，将影响山西降水的天气动力模型归纳为诊断模型，从中引出多个能够全面反映降水模型特征的综合物理因子；根据各种数值模式输出的降水量预报性能和质量优劣特点，依据数值模式的形势场预报优于要素场预报的现实，构造在不同环流形势背景下，启动不同预报方程的面雨量预报业务流程，有效地遏止了在环流形势调整时预报输出不能快速响应的弱点，提高了点和面雨量预报的准确度。     

CAPPS在太原市本地化过程中的改进和应用

该文针对CAPPS在太原市运行的误差问题，在分析其主要原因的基础上，变换地形背景场，依据不同季节、不同天气形势，在模式中增加对综合误差系数和扬尘系数的调节功能，进行多种数值模拟试验。结果表明:在误差分析基础上确定的模拟试验方案，使主要污染物SO₂、NO_X和PM₁₀的预报值与监测值的相关系数分别提高到0.96、0.91和0.89，最终使CAPPS成为预报准确率稳定在90%以上的业务化数值模式系统。     

Tetrahedral Plots of the Phase Relations for Basalts

The phase relations of quaternary systems are generally represented by projections onto ternary compositional planes. Such projections often obscure relationships that would only be evident in a three-dimensional tetrahedral plot. The tetrahedral plot requires that compositions of the minerals and melts be transformed into Cartesian coordinates. It is shown here how this transformation is carried out. The application is demonstrated by tetrahedral plots of experimental melt compositions of partially molten lherzolite. Furthermore, the plot can be used to evaluate whether or not a particular basaltic composition represents a primary melt. The methods are applicable to any four-component system.