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.     

A Modeling Study of the Effects of Anomalous Snow Cover over the Tibetan Plateau upon the South Asian Summer Monsoon

The effect of anomalous snow cover over the Tibetan Plateau upon the South Asian summer monsoon is investigated by numerical simulations using the NCAR regional climate model (RegCM2) into which gravity wave drag has been introduced. The simulations adopt relatively realistic snow mass forcings based on Scanning Multi-channel Microwave Radiometer (SNINIR) pentad snow depth data. The physical mechanism and spatial structure of the sensitivity of the South Asian early summer monsoon to snow cover anomaly over the Tibetan Plateau are revealed. The main results are summarized as follows. The heavier than normal snow cover over the Plateau can obviously reduce the shortwave radiation absorbed by surface through the albedo effect, which is compensated by weaker upward sensible heat flux associated with colder surface temperature, whereas the effects of snow melting and evaporation are relatively smaller.The anomalies of surface heat fluxes can last until June and become unobvious in July. The decrease of the Plateau surface temperature caused by heavier snow cover reaches its maximum value from late April to early May. The atmospheric cooling in the mid-upper troposphere over the Plateau and its surrounding areas is most obvious in May and can keep a fairly strong intensity in June. In contrast, there is warming to the south of the Plateau in the mid-lower troposphere from April to June with a maximum value in May.The heavier snow cover over the Plateau can reduce the intensity of the South Asian summer monsoon and rainfall to some extent, but this influence is only obvious in early summer and almost disappears in later stages.     

人工料场的开挖及安全防护措施

结合宝泉抽水蓄能电站龟山石料场开挖施工，介绍了料场开挖施工方案和道路选择，以及相应的施工安全防护体系。     

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.     

2000-2015年三江平原主要作物需水量征及影响因素分析

全生育期内作物需水量的研究是农业水资源有效利用和进行合理灌溉的重要依据。基于三江平原22个气象站点2000—2015年逐日气象观测资料及中国区域地面气象要素数据集，利用国际粮农组织 （FAO）Penman-Monteith模型和分段单值平均作物系数法，分别对三江平原水稻、玉米和大豆的作物需水量进行计算，分析作物需水量年际变化特征，采用通径分析法研究作物需水量的变化成因。结果表明：（1）三江平原16 a来年均参考作物蒸散量为537.4 mm，日均为 3.5 mm，呈波动减少趋势。（2）生长季内，水稻在分蘖期需水量最大，为177.1 mm，玉米在七叶期需水量最大，为99.7 mm，大豆在结荚期需水量最大，为96.1 mm；水稻、玉米和大豆的净灌溉需水量分别为195.4 mm、130.8 mm和72.2 mm，对灌溉的依赖程度水稻>玉米>大豆。（3）由通径分析结果可知，三江平原作物需水量的主要影响因素为净辐射、气温和日照时数。     

Analysis of seafloor backscatter strength dependence on the survey azimuth using multibeam echosounder data

The sediment backscatter strength measured by multibeam echosounders is a key feature for seafloor mapping either qualitative (image mosaics) or quantitative (extraction of classifying features). An important phenomenon, often underestimated, is the dependence of the backscatter level on the azimuth angle imposed by the survey line directions: strong level differences at varying azimuth can be observed in case of organized roughness of the seabed, usually caused by tide currents over sandy sediments. This paper presents a number of experimental results obtained from shallow-water cruises using a 300-kHz multibeam echosounder and specially dedicated to the study of this azimuthal effect, with a specific configuration of the survey strategy involving a systematic coverage of reference areas following “compass rose” patterns. The results show for some areas a very strong dependence of the backscatter level, up to about 10-dB differences at intermediate oblique angles, although the presence of these ripples cannot be observed directly—neither from the bathymetry data nor from the sonar image, due to the insufficient resolution capability of the sonar. An elementary modeling of backscattering from rippled interfaces explains and comforts these observations. The consequences of this backscatter dependence upon survey azimuth on the current strategies of backscatter data acquisition and exploitation are discussed.