Experimental solidification of anhydrous latitic and trachytic melts at different cooling rates: The role of nucleation kinetics

Two sets of cooling experiments were run at atmospheric conditions for two anhydrous starting latitic and trachytic melts: 1) five cooling rates (25, 12.5, 3, 0.5, and 0.125 °C/min) between 1300° and 800 °C, and 2) a 0.5 °C/min cooling rate from 1300 °C with quench temperatures at 1200°, 1100°, 1000° and 900 °C. Trachytic run-products are invariably glassy. Nucleation is also suppressed in the latitic run-products at the three highest cooling rates. Conversely, in the 0.5 and 0.125 °C/min runs, latites have a crystal content of  90 vol.%. The phases are: plagioclase, clinopyroxene, glass and iron-bearing oxide (in order of abundance). The variable quench temperatures, investigated by coupling experiments with Pt wire and Pt capsule sample containers in set 2, again did not produce crystallization of trachyte, whereas latitic samples are characterized by  10 vol.% of oxides, pyroxenes and plagioclase (in order of appearance), at temperature < 1000 °C. Effects of (preferential) heterogeneous nucleation on sample holders, of superheating degree, and chemical species loss during cooling are absent for both melt compositions. The difference of solidification paths between these two silicate melts can be ascribed only to their small chemical differences. In comparison with calculated equilibrium conditions all the experimental latitic and trachytic run-products revealed strong kinetic effects, interpretable in the light of the nucleation theory. The glass-forming ability (GFA) of trachyte is higher, whereas their critical cooling rate (Rc) is lower (< 0.125 °C/min), in comparison to latitic melts (Rc > 0.5 °C/min). The experimental results carried out in this study can be applied to lava flows and domes; trachytic lavas are able to flow for longer period with respect to latitic ones in a metastable condition. Glass-rich terrestrial lavas, i.e. obsidians, can be the result of sluggish nucleation kinetics due to the relative high polymerisation of evolved silicate melts.     

非均匀陆面条件下区域蒸散量计算的遥感模型

非均匀陆面条件下的区域蒸散计算是一个复杂的问题。文中首先在利用遥感资料求取地表特征参数 (如植被覆盖度、地表反照率等 )的基础上 ,建立了裸露地表条件下的裸土蒸发和全植被覆盖条件下植被蒸腾计算模型 ,然后结合植被覆盖度 (植被的垂直投影面积与单位面积之比 )给出非均匀陆面条件下的区域蒸散计算方法。实测资料验算表明该模型具有较高的计算精度。文章最后利用该模型对中国北方地区的蒸散量进行了计算 ,并对该研究区蒸散的特点进行了分析     

Signatures in flowing fluid electric conductivity logs

Flowing fluid electric conductivity logging provides a means to determine hydrologic properties of fractures, fracture zones, or other permeable layers intersecting a borehole in saturated rock. The method involves analyzing the time-evolution of fluid electric conductivity (FEC) logs obtained while the well is being pumped and yields information on the location, hydraulic transmissivity, and salinity of permeable layers. The original analysis method was restricted to the case in which flows from the permeable layers or fractures were directed into the borehole (inflow). Recently, the method was adapted to permit treatment of both inflow and outflow, including analysis of natural regional flow in the permeable layer. A numerical model simulates flow and transport in the wellbore during flowing FEC logging, and fracture properties are determined by optimizing the match between simulation results and observed FEC logs. This can be a laborious trial-and-error procedure, especially when both inflow and outflow points are present. Improved analyses methods are needed. One possible tactic would be to develop an automated inverse method, but this paper takes a more elementary approach and focuses on identifying the signatures that various inflow and outflow features create in flowing FEC logs. The physical insight obtained provides a basis for more efficient analysis of these logs, both for the present trial and error approach and for a potential future automated inverse approach. Inflow points produce distinctive signatures in the FEC logs themselves, enabling the determination of location, inflow rate, and ion concentration. Identifying outflow locations and flow rates typically requires a more complicated integral method, which is also presented in this paper.     

非均匀介质热蠕变流动的数值求解

针对非均匀介质中热蠕变流动问题,给出了有限单元方法与网格－粒子方法联合求解新技术,即有限单元方法求解欧拉网格节点上的未知量,分布于单元内部作为物质成分标记的粒子反映变形过程.有限元法求解动量方程和连续性方程时引入了速度场和压力场等阶插值的压力场稳定的Petrov Galerkin方法,求解能量方程时采用了流线迎风Petrov Galerkin方法,网格－粒子算法中采用双线性插值与有限单元插值函数对应.有限单元计算与网格－粒子计算相对独立,两种方法计算的数据通过有限单元节点传递.同时,实现了三角形单元的算法和程序,解决了复杂结构条件下不规则网格计算的问题.通过经典方腔热对流问题验证了程序,给出了不规则形态块体沉降算例,并分析了数值解的稳定性.     

矿震成核过程的公里尺度研究

成核过程是地震机理研究的重要部分.通过岩石实验、高温高压和数值模拟认为,地震成核的临界尺度为大地震总破裂长度的十分之一,一般为10~20 km左右;地震成核时间从年到秒、甚至毫秒.本文利用山东省宁阳县华丰煤矿立体台网的岩爆观测,进行了公里尺度的破裂成核过程实验研究.通过对华丰煤矿1995年2月~1996年5月期间的地脉动记录和12000次小震记录资料的分析,计算了它们的地脉动、S波和P波振幅比以及S波波峰的到时差,得到ML=2.9矿震的成核时间为50天,成核临界尺度为5 km.     

A garnet population in Yellowknife schist, Canada

Abstract Data are presented on a garnet population in a specimen of garnet-biotite-plagioclase-quartz schist from the cordierite zone of an Archaean thermal dome in the Southern Slave Province of the Canadian Shield. Garnet crystals are bounded by planar dodecahedral faces and by trapezohedral faces which on the 10-μm scale are corrugated. Crystal distribution, as revealed by dissection of a small cubic volume of rock, is random. The size distribution is normal, with a mean diameter of 0.81 mm and a standard deviation of 0.32 mm. In the largest crystal of the population (mean radius 0.83 mm), [Mn] = 100 Mn/(Fe + Mg + Mn + Ca) decreases from 14.5 at the centre to 7.5 and then increases in the outer margin to 8.5; [Fe] increases continuously from 67 at the centre to 77 at the surface; [Mg] increases from 12.5 to 13.5 and then falls sharply to 11; [Ca] remains unchanged at 4.0 and then drops to 3.3. Progressively smaller crystals have progressively lower [Mn] and higher [Fe] concentrations at their centres, while all crystals have the same margin composition. Growth vectors extending from given concentration contours to crystal surfaces are of equal length regardless of the size of the crystal in which the vector is located. A garnet-forming model is presented in which reaction was initiated by a rise in temperature. Nucleation sites were randomly selected. The nucleation rate increased with time and then declined. Crystal faces advanced at a constant linear rate, which implies an increase in volume proportional to surface area. Initially, the composition of garnet deposited on crystal surfaces was determined by van Laar equations of equilibrium, which demanded the withdrawal of Mn and Fe from within chlorite crystals. This transfer reaction was then accompanied by an ion exchange reaction which moved Mn and Fe to garnet surfaces from biotite, in exchange for Mg. The exchange reaction provides an explanation for the high overall concentration of Mn and Fe in garnet and for the observed Mn and Mg reversals in the margins of crystals. The increase of garnet volume in the garnet population is found to be parabolic, i.e. Vαα⁵.     

The influence of aerosol particle composition on cloud droplet formation

A difference in partitioning between cloud droplets and interstitial air for two chemical species (elemental carbon and sulphur) with different expected behaviour in nucleation scavenging was observed in clouds at Mt. Kleiner Feldberg (825 m asl), near Frankfurt, Germany. The fraction of sulphur incorporated in cloud droplets was always higher than the fraction of elemental carbon. This difference in partitioning has also been observed in fog but under different pollution conditions in the Po Valley, Italy. Both these studies were based on bulk samples. In the present study at Kleiner Feldberg, impactor samples of the particles in the interstitial air and the cloud droplet residuals were taken and a single particle analysis was done on the samples. It was found that, for a given particle size, the majority of particles forming cloud droplets were soluble and that insoluble particles preferentially remained in the interstitial air.     

200hPa非地转风和非均匀层结与梅雨暴雨的关系

本文通过1991年梅雨后期，7月1—13日200hPa非地转场分析得出，梅雨期200hPa非地转风主要由积云对流动量转换形成。其次是风速平流，并通过数值试验得以证明。非地转风形成的高层辐散，可进一步使暴雨增强与维持。另外梅雨期非均匀层结对非地转风及暴雨也有较大的影响。     

Digital Image Based Approach for Three-Dimensional Mechanical Analysis of Heterogeneous Rocks

Summary This paper presents a digital image based approach for three-dimensional (3-D) numerical simulation and failure analysis of rocks by taking into account the actual 3-D heterogeneity. Digital image techniques are adopted to extract two-dimensional (2-D) material heterogeneity from material surface images. The 2-D image mesostructures are further extrapolated to 3-D cuboid mesostructures by assuming the material surface as a representation of the inner material heterogeneity within a very small depth. The iterative milling and scanning system is set up to generate the 3-D rock mesostructures. A Hong Kong granite specimen is used as an example to demonstrate the procedure of 3-D mesostructure establishment. The mechanical responses and failure process under the conventional Brazilian tensile test condition are examined through numerical analyses. The stress distribution, crack propagation process and failure model of heterogeneous material cases are simulated with a finite difference software. The numerical results indicate that material heterogeneity plays an important role in determining the failure behavior of rocks under external loading.     

Density-dependent dispersion in heterogeneous porous media Part I: A numerical study

In this paper, we describe carefully conducted numerical experiments, in which a dense salt solution vertically displaces fresh water in a stable manner. The two-dimensional porous media are weakly heterogeneous at a small scale. The purpose of these simulations, conducted for a range of density differences, is to obtain accurate concentration profiles that can be used to validate nonlinear models for high-concentration-gradient dispersion. In this part we focus on convergence of the computations, in numerical and statistical sense, to ensure that the uncertainty in the results is small enough.Concentration variances are computed, which give estimates of the uncertainty in local concentration values. These local variations decrease with increasing density contrast. For tracer transport, obtained longitudinal dispersivities are in accordance with analytical findings. In the case of high-density contrasts, stabilizing gravity forces counteract the growth of dispersive fingers, decreasing the effective width of the transition zone. For small log-permeability variances, the decrease of the apparent dispersivity that is found is in agreement with laboratory results for homogeneous columns.