A virtual experiment technique on the elementary behaviour of granular materials with discrete element method

Multi‐scale investigations aided by the discrete element method (DEM) play a vital role for current state‐of‐the‐art research on the elementary behaviour of granular materials. Similar to laboratory tests, there are three important aspects to be considered carefully, which are the proper stress/strain definition and measurement, the application of target loading paths and the designed experiment setup, to be addressed in the present paper. Considering the volume sensitive characteristics of granular materials, in the proposed technique, the deformation of the tested specimen is controlled and measured by deformation gradient tensor involving both the undeformed configuration and the current configuration. Definitions of Biot strain and Cauchy stress are adopted. The expressions of them in terms of contact forces and particle displacements, respectively, are derived. The boundary of the tested specimen consists of rigid massless planar units. It is suggested that the representative element uses a convex polyhedral (polygonal) shape to minimize possible boundary arching effects. General loading paths are described by directly specifying the changes in the stress/strain invariants or directions. Loading can be applied in the strain‐controlled mode by specifying the translations and rotations of the boundary units, or in the stress‐controlled mode by using a servo‐control mechanism, or in the combination of the two methods to realize mixed boundary conditions. Taking the simulation results as the natural consequences originated from a complex system, virtual experiments provide particle‐scale information database to conduct multi‐scale investigations for better understanding in granular material behaviours and possible development of the constitutive theories provided the qualitative similarity between the simulation results from virtual experiments and observations on real material behaviour. Copyright © 2011 John Wiley & Sons, Ltd.     

Influence of spatial variations of microtopography and infiltration on surface runoff and field scale hydrological connectivity

Surface runoff on agricultural fields arises when rainfall exceeds infiltration. Excess water ponding in and flowing through local microtopography increases the hydrological connectivity of fields. In turn, an increased level of hydrological connectivity leads to a higher surface runoff flux at the field boundaries. We investigated the functional hydrological connectivity of synthetical elevation fields with varying statistical properties. For this purpose, we developed an object-oriented ponding and redistribution model to which Philip’s infiltration model was coupled. The connectivity behaviour is determined by the presence of depressions with a large area and spatial organization of microtopography in rills or channels. The presence of microdepressions suppresses the effect of the spatial variation of infiltration properties. Connectivity behaviour of a field with a varying spatial distribution of infiltration properties can be predicted by transforming the unique connectivity function that was defined for a designated microtopography.     

Contrasting komatiite belts,associated Ni–Cu–(PGE) deposit styles and assimilation histories

The Agnew–Wiluna greenstone belt in the Yilgarn Craton of Western Australia is the most nickel-sulfide-endowed komatiite belt in the world. The Agnew–Wiluna greenstone belt contains two mineralised units/horizons that display very different volcanological and geochemical features. The Mt Keith unit comprises >500 m-thick spinifex-free adcumulate-textured lenses, which are flanked by laterally extensive orthocumulate-textured units. Spinifex texture is absent from this unit. Disseminated nickel sulfides, interstitial to former olivine crystals, are concentrated in the lensoidal areas. Massive sulfides are locally present along the base or margins of the lenses or channels. The Cliffs unit is locally >150 m thick and comprises a sequence of differentiated spinifex-textured flow units. The basal unit is the thickest, and contains basal massive nickel-sulfide mineralisation. The Mt Keith and Cliffs units display important common features: (i) MgO contents of 25–30% in inferred parental magmas; and (ii) Al/Ti ratios of ～20 (Munro-type). However, the Mt Keith unit is highly crustally contaminated (e.g. LREE-enriched, high HFSEs), whereas the Cliffs unit does not display evidence of significant crustal assimilation. We argue that the distinct trace-element concentrations and profiles of the two komatiite units reflect their different emplacement style and country rocks: the Mt Keith unit is interpreted to have been emplaced as an intrusive sill into dacitic volcanic units whereas the Cliffs unit was extruded as lava flow onto tholeiitic basalts in a subaqueous environment. The mode of emplacement and nature of country rock is the single biggest factor in controlling mineralisation styles in komatiites. On the other hand, evidence of crustal contamination does not necessarily provide information of the prospectivity of komatiites to host Ni–Cu–(PGE) mineralisation, despite being a good proxy for the style of komatiite emplacement and the nature of country rocks.     

通过改变光热条件分析胡杨群落光合作用对土壤呼吸速率的影响

结合干旱环境下塔里木河下游地区胡杨群落长势、土壤呼吸、土壤水分、温度的实时监测资料,通过改变光热条件对胡杨群落进行人工遮阴控制,分析植物群落长势对土壤呼吸速率的影响,同时对不同光环境下胡杨群落土壤呼吸速率变化特征及其影响因素进行了探讨。结果表明,胡杨群落土壤呼吸速率日变化呈现多峰值曲线,但最大值集中出现在12:00—16:00之间,最低值出现在早上08:00和晚上20:00;植物群落长势很大程度上决定着土壤呼吸速率的大小;胡杨群落土壤呼吸的日变化与水热因子之间联系存在很大的不确定性,但在很大程度上受植物新近合成同化产物控制。     

区域资料四维同化试验

根据国家气象中心有限区分析预报系统的基本方案，设计一个以6小时为周期的间歇区域资料四维同化方案，对我国的两次暴雨过程进行两种不同粗细水平分辨率的资料四维同化和短期模式预报的试验。结果表明，由同化方案初值获得的短期预报明显好于LAFS方案的预报结果，高分辨率同化方案对预报的改进更明显。     

Using a H∞ filter assimilation procedure to estimate root zone soil water content

Root zone soil water content impacts plant water availability, land energy and water balances. Because of unknown hydrological model error, observation errors and the statistical characteristics of the errors, the widely used Kalman filter (KF) and its extensions are challenged to retrieve the root zone soil water content using the surface soil water content. If the soil hydraulic parameters are poorly estimated, the KF and its extensions fail to accurately estimate the root zone soil water. The H‐infinity filter (HF) represents a robust version of the KF. The HF is widely used in data assimilation and is superior to the KF, especially when the performance of the model is not well understood. The objective of this study is to study the impact of uncertain soil hydraulic parameters, initial soil moisture content and observation period on the ability of HF assimilation to predict in situ soil water content. In this article, we study seven cases. The results show that the soil hydraulic parameters hold a critical role in the course of assimilation. When the soil hydraulic parameters are poorly estimated, an accurate estimation of root soil water content cannot be retrieved by the HF assimilation approach. When the estimated soil hydraulic parameters are similar to actual values, the soil water content at various depths can be accurately retrieved by the HF assimilation. The HF assimilation is not very sensitive to the initial soil water content, and the impact of the initial soil water content on the assimilation scheme can be eliminated after about 5–7 days. The observation interval is important for soil water profile distribution retrieval with the HF, and the shorter the observation interval, the shorter the time required to achieve actual soil water content. However, the retrieval results are not very accurate at a depth of 100 cm. Also it is complex to determine the weighting coefficient and the error attenuation parameter in the HF assimilation. In this article, the trial‐and‐error method was used to determine the weighting coefficient and the error attenuation parameter. After the first establishment of limited range of the parameters, ‘the best parameter set’ was selected from the range of values. For the soil conditions investigated, the HF assimilation results are better than the open‐loop results. Copyright © 2010 John Wiley & Sons, Ltd.     

Phase Equilibria of the Lyngdal Granodiorite (Norway): Implications for the Origin of Metaluminous Ferroan Granitoids

The Proterozoic (950 Ma) Lyngdal granodiorite of southern Norwaybelongs to a series of hornblende–biotite metaluminousferroan granitoids (HBG suite) coeval with the post-collisionalRogaland Anorthosite–Mangerite–Charnockite (AMC)suite. This granitoid massif shares many geochemical characteristicswith rapakivi granitoids, yet granodiorites dominate over granites.To constrain both crystallization (P, T, fO₂, H₂O in melt) andmagma generation conditions, we performed crystallization experimentson two samples of the Lyngdal granodiorite (with 60 and 65 wt% SiO₂) at 4–2 kbar, mainly at fO₂ of NNO (nickel–nickeloxide) to NNO + 1, and under fluid-saturated conditions withvarious H₂O–CO₂ ratios for each temperature. Comparisonbetween experimental phase equilibria and the mineral assemblagein the Lyngdal granodiorite indicates that it crystallized between4 and 2 kbar, from a magma with 5–6 wt % H₂O at an fO₂of NNO to NNO + 1. These oxidized and wet conditions sharplycontrast with the dry and reduced conditions inferred for thepetrogenesis of the AMC suite and many other rapakivi granitesworldwide. The high liquidus temperature and H₂O content ofthe Lyngdal granodiorite imply that it is not a primary magmaproduced by the partial melting of the crust but is derivedby the fractionation of a mafic magma. Lyngdal-type magmas appearto have volcanic equivalents in the geological record. In particular,our results show that oxidized high-silica rhyolites, such asthe Bishop Tuff, could be derived via fractionation of oxidizedintermediate magmas and do not necessarily represent primarycrustal melts. This study underlines the great variability ofcrystallization conditions (from anhydrous to hydrous and reducedto oxidized) and petrogenetic processes among the metaluminousferroan magmas of intermediate compositions (granodiorites,quartz mangerites, quartz latites), suggesting that there isnot a single model to explain these rocks. KEY WORDS: ferroan granitoids; crystallization conditions; experiments; Norway; Sveconorwegian; Bishop Tuff     

在数值预报中应用卫星反演资料的质量控制研究

资料变分同化作为提高数值预报准确率的有效手段之一，受到了越来越多的关注。如何将分辨率较高的卫星资料应用到数值预报中，弥补常规观测网时空分辨率的不足，成为数值天气预报的重要课题。以1996年6月28-29日发生在江淮流域的暴雨过程为例，依据其中尺度场和能量场特征，分别以低云量、急流、湿静力能为标准设计了三套质量控制方案。将质量控制后的反演湿度场变分同化并应用到中尺度MM5模式中，结果表明：经质量控制后，可以明显地改善降水区域及降水量的预报；对降水预报，特别是强降水天气过程的预报，该质量控制方案具有普遍意义。     

Particularités de la contamination crustale des phonolites : exemple du Velay oriental (Massif central)

Sr and Nd isotopic compositions of one trachyte, eight phonolites and five basalts have been measured. The isotopic characteristics of the trachyte can be explained by a combined assimilation–fractional crystallization process within an upper crustal magmatic chamber. Some phonolites display isotopic signatures identical to basalts, suggesting that they have been protected against any crustal assimilation during their formation. Some others have low Sr contents, whereas they are enriched in radiogenic Sr (0.70451<⁸⁷Sr/⁸⁶Sr_i<0.71192), and display basaltic ¹⁴³Nd/¹⁴⁴Nd ratios. Both observations could be explained by very strong alkali feldspar fractionation and by subsequent very low assimilation of surrounding rocks (between 0.3 and 4%) during intrusion. To cite this article: J.-M. Dautria et al., C. R. Geoscience 336 (2004).     

多普勒雷达径向速度同化在淮河暴雨数值模拟中的应用

针对2007年7月淮河流域的一次强降雨过程,利用WRF中尺度数值模式及其三维变分同化系统(WRF-3DVAR),开展了多普勒雷达径向速度的三维变分同化对暴雨过程模拟效果的影响研究。结果表明:WRF-3DVAR能够有效地同化多普勒雷达径向速度资料,同化后使得模式初始场出现了一定的调整,包含更详尽的中尺度特征信息,进而显著改善模式对大暴雨过程前12h降水的模拟效果。在高分辨率中尺度数值模式中有效地利用多普勒天气雷达资料,能较好地提高中尺度降雨预报。