基于SHAW模型的黄土高原半干旱区农田土壤水分动态模拟

黄土高原半干旱区土壤蒸发强烈,准确地掌握土壤水分动态对于旱地农业水分管理至关重要。应用基于物理基础的一维水热耦合SHAW(The Simultaneous Heat and Water)模型,模拟了陕西子洲岔巴沟流域1964~1967年土壤水分和土壤蒸发的动态特征,以及神木六道沟流域2006年坡地和梯田土壤水分变化。结果表明,除表层土壤水分模拟结果偏差较大,其他土层模拟值与实测值基本吻合,模拟期土壤水分模拟的相对平均绝对误差(Relatively Mean Absolutely Error,RMAE)为5.2%~11.4%。1964~1967年土壤累积蒸发量模拟值与实测值平均相对偏差为0.8%~6.1%,土壤蒸发的模拟值与实测值较为一致。因此,SHAW模型可以用于黄土高原半干旱区农田土壤水分动态规律研究。

黄土高原水蚀风蚀交错区土壤剖面水分动态的数值模拟研究

利用SWAP(Soil-Water-Atmosphere-Plant)模型对黄土高原水蚀风蚀交错区坡地土壤-植被-大气系统中的水循环进行数值模拟。结果显示,SWAP模型很好的模拟了不同土地利用方式条件下的土壤水循环过程。根据模拟结果,水蚀风蚀交错区的丰水年份,农地和种植第一年的紫花苜蓿地季末土壤水分稍有盈余,谷子和紫花苜蓿的日蒸散量分别为1.2~2.6 mm和1.2~2.5 mm。

解集GCM输出模拟黄土塬区土壤水分平衡的潜在变化

壤水分平衡对气候变化存在着响应,在全球变化的大背景下,研究土壤水分的可能变化是气候变化影响评估中非常重要的一项内容。目标是利用经验统计方法解集GCM网格逐月的降水和温度数据,并使用建立的气候变化情景作为WEPP的输入文件评估黄土高原王东沟流域2010～2039年土壤水分平衡(土壤水分、蒸发、渗漏和蒸腾)的可能变化。结果表明,3种情景预测2010～2039年王东沟流域年均降水可能增长1.8%～17.5%,年最高温度和最低温度分别可能增长0.5～0.9℃和2.0～2.3℃。作物蒸腾变化主要在4～6月份,土壤蒸发变化主要发生在7～9月份;作物蒸腾年均变化-5%～19%,土壤水分年均变化-4%～4%,土壤蒸发年均变化-7%～7%,均为A2a减少,B2a和GGal增大;A2a的土壤水分渗漏增长最大,GGal次之,B2a基本不变。这些结果表明气候变化及其导致的作物生长变化对土壤水分平衡存在重要的影响。

不同坡度草地含沙水流水力学特性及其拦沙机理

参照黄土区侵蚀降雨和坡面片蚀产沙特征,采用恒流泥沙输送装置模拟坡面上方来水来沙,探讨不同坡度草地含沙水流的水力学特性及其对上方来沙的拦蓄机理。结果表明,草地坡面的水流弗劳德数随坡度增大而增加,而阻力系数与坡度呈反势。按明渠水流的一般标准,不同坡度草地水流均为层状缓流。草地坡面拦沙效应随坡度增大而减小,且径流前期的减沙作用较后期更为显著。不同坡度草地坡面的出流泥沙平均直径和大颗粒(>10μm)泥沙含量均显著小于上方来沙,这说明草地的拦沙效应主要体现在对大粒径泥沙的拦蓄上。

降水量的BP人工神经网络预测模型及其应用

由于影响因素的复杂性 ,预测降水量具有相当的难度。在假设区域长时间内降水量和蒸发量保持平衡的基础上 ,用 BP人工神经网络建立了陕西省汉中市的降水量预测模型 ,根据前 3个月降水量和蒸发量对降水量资料进行了模拟预测 ,结果认为其准确率为 84% ,合格率为 10 0 %。

Simulation study of China's net primary production

Spatial and temporal distribution of vegetation net primary production (NPP) in China was studied using three light-use efficiency models (CASA, GLOPEM and GEOLUE) and two mechanistic ecological process models (CEVSA, GEOPRO). Based on spatial and temporal analysis (e.g. monthly, seasonally and annually) of simulated results from ecological process mechanism models of CASA, GLOPEM and CEVSA, the following conclusions could be made: (1) during the last 20 years, NPP change in China followed closely the seasonal change of climate affected by monsoon with an overall trend of increasing; (2) simulated average seasonal NPP was: 0.571 +/- 0.2 GtC in spring, 1.573 +/- 0.4 GtC in summer, 0.6 +/- 0.2 GtC in autumn, and 0.12 +/- 0.1 GtC in winter. Average annual NPP in China was 2.864 +/- 1 GtC. All the five models were able to simulate seasonal and spatial features of biomass for different ecological types in China. This paper provides a baseline for China's total biomass production. It also offers a means of estimating the NPP change due to afforestation, reforestation, conservation and other human activities and could aid people in using for-mentioned carbon sinks to fulfill China's commitment of reducing greenhouse gases.

Modelling raindrop impact and splash erosion processes within a spatial cell: a stochastic approach

A new approach is proposed to simulate splash erosion on local soil surfaces. Without the effect of wind and other raindrops, the impact of free-falling raindrops was considered as an independent event from the stochastic viewpoint. The erosivity of a single raindrop depending on its kinetic energy was computed by an empirical relationship in which the kinetic energy was expressed as a power function of the equivalent diameter of the raindrop. An empirical linear function combining the kinetic energy and soil shear strength was used to estimate the impacted amount of soil particles by a single raindrop. Considering an ideal local soil surface with size of I m x I m, the expected number of received free-failing raindrops with different diameters per unit time was described by the combination of the raindrop size distribution function and the terminal velocity of raindrops. The total splash amount was seen as the sum of the impact amount by all raindrops in the rainfall event. The total splash amount per unit time was subdivided into three different components, including net splash amount, single impact amount and re-detachment amount. The re-detachment amount was obtained by a spatial geometric probability derived using the Poisson function in which overlapped impacted areas were considered. The net splash amount was defined as the mass of soil particles collected outside the splash dish. It was estimated by another spatial geometric probability in which the average splashed distance related to the median grain size of soil and effects of other impacted soil particles and other free-falling raindrops were considered. Splash experiments in artificial rainfall were carried out to validate the availability and accuracy of the model. Our simulated results suggested that the net splash amount and re-detachment amount were small parts of the total splash amount. Their proportions were 0.15% and 2.6%, respectively. The comparison of simulated data with measured data showed that this model could be applied to simulate the soil-splash process successfully and needed information of the rainfall intensity and original soil properties including initial bulk intensity, water content, median grain size and some empirical constants related to the soil surface shear strength, the raindrop size distribution function and the average splashed distance. Copyright (c) 2007 John Wiley & Sons, Ltd.

A Video Quality Assessment Metric Based on Human Visual System

It is important for practical application to design an effective and efficient metric for video quality. The most reliable way is by subjective evaluation. Thus, to design an objective metric by simulating human visual system (HVS) is quite reasonable and available. In this paper, the video quality assessment metric based on visual perception is proposed. Three-dimensional wavelet is utilized to decompose video and then extract features to mimic the multichannel structure of HVS. Spatio-temporal contrast sensitivity function (S-T CSF) is employed to weight coefficient obtained by three-dimensional wavelet to simulate nonlinearity feature of the human eyes. Perceptual threshold is exploited to obtain visual sensitive coefficients after S-T CSF filtered. Visual sensitive coefficients are normalized representation and then visual sensitive errors are calculated between reference and distorted video. Finally, temporal perceptual mechanism is applied to count values of video quality for reducing computational cost. Experimental results prove the proposed method outperforms the most existing methods and is comparable to LHS and PVQM.

Numerical simulation of mesoscopic mechanical behaviors of gradual multi-fiber-reinforced polymer matrix composites

To simulate the deformation and the fracture of gradual multi-fiber-reinforced matrix composites, a numerical simulation method for the mesoscopic mechanical behaviors was developed on the basis of the finite element and the Monte Carlo methods. The results indicate that the strength of a composite increases if the variability of statistical fiber strengths is decreased.

Stress transfer and damage evolution simulations of fiber-reinforced polymer-matrix composites

The stress transfer from broken fibers to unbroken fibers in fiber-reinforced thermosetting polymer-matrix composites and thermoplastic polymer-matrix composites was studied using a detailed finite element model. In order to check the validity of this approach, an epoxy-matrix monolayer composite was used as thermosetting polymer-matrix composite and a polypropylene (PP)-matrix monolayer composite was used as thermoplastic polymer-matrix composite, respectively. It is found that the stress concentrations near the broken fiber element cause damage to the neighboring epoxy matrix prior to the breakage of other fibers, whereas in the case of PP-matrix composites the fibers nearest to the broken fiber break prior to the PP matrix damage, because the PP matrix around the broken fiber element yields. In order to simulate composite damage evolution, a Monte Carlo technique based on a finite element method has been developed in the paper. The finite element code coupled with statistical model of fiber strength specifically written for this problem was used to determine the stress redistribution. Five hundred samples of numerical simulation were carried out to obtain statistical deformation and failure process of composites with fixed fiber volume fraction.

Computer simulation of human iodothyronine deiodinase single-chain antibody

Iodothyronine plays a major role in growth, basic metabolism and organ formation. It has an extremely limited source in the body. In this thesis, we designed iodothyronine(T4) as hapten. Then a single chain antibody displayed on phange was obtained from a human phage displaying a single chain antibody library. The specific genes of E3 was subcloned in P-5E vector. According to its amino acid sequences, we simulate its three dimention structure by computer. It has never been reported in PDB.

Digital simulation in a thin layer spectroelectrochemical cell with minigrid platinum electrode by microregion approximation explicit finite difference method

The microregion approximation explicit finite difference method is used to simulate cyclic voltammetry of an electrochemical reversible system in a three-dimensional thin layer cell with minigrid platinum electrode. The simulated CV curve and potential scan-absorbance curve were in very good accordance with the experimental results, which differed from those at a plate electrode. The influences of sweep rate, thickness of the thin layer, and mesh size on the peak current and peak separation were also studied by numerical analysis, which give some instruction for choosing experimental conditions or designing a thin layer cell. The critical ratio (1.33) of the diffusion path inside the mesh hole and across the thin layer was also obtained. If the ratio is greater than 1.33 by means of reducing the thickness of a thin layer, the electrochemical property will be far away from the thin layer property.

Simulation of characteristic electrochemical behavior of azobenzene derivative self-assembled monolayers

The characteristic electrochemical mechanics of azobenzene derivative self-assembled monolayers is discussed in present paper. It is presented that the structure inhibition is one of the most important factors in the increase of electrochemical reactive energy. A corresponding mathematical model was established based on Levich and Marcus's theory. Moreover, computational program was written to simulate the decrease of apparent rate constant (k(app)) of electron transfer with increasing surface concentration.

Radiation effects on the immiscible polymer blend of nylon1010 and high-impact polystyrene (HIPS) I: Gel/dose curves, mathematical expectation theorem and thermal behaviour

This paper studies the radiation properties of the immiscible blend of nylon1010 and HIPS. The gel fraction increased with increasing radiation dose. The network was found mostly in nylon1010, the networks were also found in both nylon1010 and HIPS when the dose reaches 0.85 MGy or more. We used the Charleby-Pinner equation and the modified Zhang-Sun-Qian equation to simulate the relationship with the dose and the sol fraction. The latter equation fits well with these polymer blends and the relationship used by it showed better linearity than the one by the Charleby-Pinner equation. We also studied the conditions of formation of the network by the mathematical expectation theorem for the binary system. Thermal properties of polymer blend were observed by DSC curves. The crystallization temperature decreases with increasing dose because the cross-linking reaction inhibited the crystallization procession and destroyed the crystals. The melting temperature also reduced with increasing radiation dose. The dual melting peak gradually shifted to single peak and the high melting peak disappeared at high radiation dose. However, the radiation-induced crystallization was observed by the heat of fusion increasing at low radiation dose. On the other hand, the crystal will be damaged by radiation. A similar conclusion may be drawn by the DSC traces when the polymer blends were crystallized. When the radiation dose increases, the heat of fusion reduces dramatically and so does the heat of crystallization. (C) 1999 Elsevier Science Ltd. All rights reserved.

Determination of surface pK(a) of SAM using an electrochemical titration method

A self-assembled monolayer (SAM) of 3-mercaptopropionic acid (HSCH2CH2COOH) was formed on a gold electrode. The effect of the charge of the end group on the electrochemical response of Fe(CN)(6)(3-) at the SAM modified electrode was studied by using cyclic voltammetry. At high pH, when the -COOH groups are dissociated, Fe(CN)(6)(3-) current is suppressed; as the solution pH is lowered, Fe(CN)(6)(3-) current increases. The electrochemical titration curve was obtained by correlating the currents to the different electrolyte pH values, from which the surface pK(a) was obtained to be 5.2+/-0.1. Furthermore, a calculation equation was presented to simulate the electrochemical titration. As comparison, the surface pK(a) was also measured by contact angle titration as 5.6+/-0.1. The surface pK(a) values determined by the two methods in our work are consistent and accurate.