A method for extending a general constitutive model to consider the electro-chemo-mechanical phenomena of mineral crystals in expansive soils

The mechanical behavior of expansive soils varies according to the chemical composition of the pore fluid. It is well known that electrochemical phenomena on the surface of clay mineral crystals considerably affect their macrostructural behavior. In particular, a change in the pore fluid composition causes osmotic consolidation or swelling. In this study, a model is constructed to describe the characteristic behavior of expansive soils by coupling the interlaminar behavior of clay mineral crystals and the soil skeleton behavior. The interlaminar behavior is derived from the electro-chemo-mechanical coupling equilibrium of mineral crystals, and the soil skeleton behavior is given by a general elastoplastic constitutive model for soils. This modeling approach extends a general model to consider the electro-chemo-mechanical phenomena of mineral crystals. Oedometer tests with the replacement of the cell fluid and the swelling pressure and deformation tests on expansive clays are simulated with the proposed model. The simulation results indicate that the proposed method can reasonably represent the typical behavior of expansive soils.     

Sorption of selenite in a multi-component system using the “dialysis membrane” method

⁷⁹Se is a potentially mobile long-lived fission product, which may make a dominant contribution to the long-term radiation exposure resulting from deep geological disposal of radioactive waste. Its mobility is affected by sorption on minerals. Selenium sorption processes have been studied mainly by considering interaction with a single mineral surface. In the case of multi-component systems (e.g. soils), it is difficult to predict the radioelement behaviour only from the mineral constituents. This study contributes to the understanding of multi-component controls of Se concentrations towards predicting Se behaviour in soils after migration from a disposal site. This goal was approached by measuring selenite sorption on mono and multi-phase systems physically separated by dialysis membranes. To the best of the authors’ knowledge, very few studies have used dialysis membranes to study the sorption competition of selenite between several mineral phases. Other workers have used this method to study the sorption of pesticides on montmorillonite in the presence of dissolved organic matter. Indeed, this method allows measurement of individual K_d in a system composed of several mineral phases. Dialysis membranes allowed (i) determination of the competition of two mineral phases for selenite sorption (ii) and determination of the role of humic acids (HAs) on selenite sorption in oxidising conditions. Experimental results at pH 7.0 show an average Se(IV) sorption distribution coefficient (K_d) of approximately 125 and 9410 L kg⁻¹ for bentonite and goethite, respectively. The average K_d for goethite decreases to 613 L kg⁻¹ or 3215 L kg⁻¹ in the presence of bentonite or HA, respectively. For bentonite, the average K_d decreases slightly in the presence of goethite (60 L kg⁻¹) and remains unchanged in the presence of HA. The experimental data were successfully modelled with a surface complexation model using the PHREEQC geochemical code. The drastic decrease in Se(IV) sorption on goethite in a multi-phase system is attributed to competition with dissolved silica released by bentonite. As with Si the HA compete with Se for sorption sites on goethite.     

湿地“淹埋深-历时-频率”阈值研究的理论和方法

针对"淹埋深-历时-频率"(S, D, F)阈值研究存在的问题,根据土壤学、植物生态学、系统论、地理系统学说以及反演理论等的相关原理,系统分析了湿地(S, D, F)阈值确定的理论和方法。结果表明:在正常情况下,水饱和至地表的最大地下水埋深、湿地植物主要根系分布层下限和沉水植物分布下界淹水深分别是湿地土壤上界、湿地植被上界和湿地下界的S阈值;D阈值只能通过湿地土壤或湿地植被边界界定在水位过程线上反演,稳定湿地边界的D阈值应取与50%频率值对应的D值或D值的多年平均值,波动湿地边界的D阈值是常数;"非正常情况下"湿地的D阈值应参照湿地水文地貌分类的同一区域子类中"正常情况下"的湿地的D阈值;检验D阈值科学性的标准是利用D阈值确定的湿地水文边界必须与湿地土壤边界或湿地植被边界耦合。     

大断面海底隧道CRD法绝对位移控制基准建立及应用研究

隧道位移控制基准的建立是确保隧道安全施工的一项重要内容。以厦门－翔安海底隧道陆域浅埋暗挖段CRD法施工为工程背景,陆域浅埋段总长约2 525 m,主要穿越全、强风化花岗岩地带,含水率高、自稳能力差。依据全站仪非接触量测获得初期支护三维绝对位移分布规律,确定三维绝对位移基准控制指标,同时对86个围岩位移监测断面的监测数据进行变形正常断面和变形异常或险情断面的位移量统计分析,最终确定CRD法（center cross diaphragm method,交叉中隔壁法）三维绝对位移控制基准。研究结果表明,将CRD1、CRD3部拱顶竖直位移和上中隔墙、下中隔墙水平位移作为CRD法位移控制指标;CRD1部和CRD3部拱顶下沉极限位移分别为200 mm和130 mm;CRD3部拱顶竖直位移约为CRD1部的0.5～1倍;上中隔墙水平极限位移约为100 mm,下中隔墙水平极限位移约为15 mm。通过现场97个断面应用,共发现58次异常和险情,由于发现及时,处置有力,保证了施工安全。     

Stable isotopes in deep-sea corals and a new mechanism for “vital effects”

Offsets from isotopic equilibrium in biogenic carbonates have complicated paleoclimate reconstructions for decades. A new archive of climate, deep-sea corals, is used to evaluate the calcification processes, independent of photosynthesis, that contribute to these offsets. Carbon and oxygen stable isotope data from six modern deep-sea corals show strong linear trends between δ¹³C and δ¹⁸O. Slopes of these trends between samples are similar and range between 1.9 to 2.6 for Δδ¹³C/Δδ¹⁸O. Linear trends intersect isotopic equilibrium for δ¹⁸O and are slightly depleted for δ¹³C. Variations in the isotopic ratios are strongly correlated with the density banding structure. Isotopically depleted aragonite is associated with light, quickly precipitating bands, whereas isotopically enriched points correspond to slowly accumulating, less dense aragonite. The densest white band at the trabecular center is furthest from isotopic equilibrium for both carbon and oxygen. Data from this region fall off the linear trend between δ¹⁸O and δ¹³C. This deviation, where δ¹³C remains constant while the δ¹⁸O continues to decrease, does not support “vital effect” mechanisms that call upon kinetic fractionation to explain offsets from isotopic equilibrium. We propose a new mechanism for vital effects in these deep-sea corals that is based on a thermodynamic response to a biologically induced pH gradient in the calcifying region.     

A method for generating virtual soil profiles with complex,multi-layer stratigraphy

This paper presents a framework for generating multi-layer, unconditional soil profiles with complex stratigraphy, which simulates the effects of natural erosion and sedimentation processes. The stratigraphy can have varying degrees of randomness and can include features such as lenses, as well as sloped and undulating layers. The method generates the soil comprising the layers using local average subdivision (LAS), and a random noise component that is added to the layer boundaries. The layers are created by generating coordinates of key points in the simulated ground profile, which are then interpolated with a customised, 2D, linear interpolation algorithm. The resulting simulations facilitate more accurate probabilistic modelling of geotechnical engineering systems because they provide more realistic geologies, such as those usually encountered in the ground. Fortran code implementing this framework is included as supplementary material.     

Analysis of single rock blocks for general failure modes under conservative and non‐conservative forces

After describing the kinematics of a generic rigid block subjected to large rotations and displacements, the Udwadia's General Principle of Mechanics is applied to the dynamics of a rigid block with frictional constraints to show that the reaction forces and moments are indeterminate. Thus, the paper presents an incremental‐iterative algorithm for analysing general failure modes of rock blocks subject to generic forces, including non‐conservative forces such as water forces. Consistent stiffness matrices have been developed that fully exploit the quadratic convergence of the adopted Newton–Raphson iterative scheme. The algorithm takes into account large block displacements and rotations, which together with non‐conservative forces make the stiffness matrix non‐symmetric. Also included in the algorithm are in situ stress and fracture dilatancy, which introduces non‐symmetric rank‐one modifications to the stiffness matrix. Progressive failure is captured by the algorithm, which has proven capable of detecting numerically challenging failure modes, such as rotations about only one point. Failure modes may originate from a limit point or from dynamic instability (divergence or flutter); equilibrium paths emanating from bifurcation points are followed by the algorithm. The algorithm identifies both static and dynamic failure modes. The calculation of the factor of safety comes with no overhead. Examples show the equilibrium path of a rock block that undergoes slumping failure must first pass through a bifurcation point, unless the block is laterally constrained. Rock blocks subjected to water forces (or other non‐conservative forces) may undergo flutter failure before reaching a limit point. Copyright © 2007 John Wiley & Sons, Ltd.     

橄榄石水溶性与含铁量相关性的实验研究

本文首次报导了含铁量不同的橄榄石单晶在不同温压条件下水溶性的试验研究结果。橄榄石单晶试件在300MPa静水压和1000℃至1300℃的温度条件下进行热压，每隔50℃进行一组试验，氧逸度被控制在Ni/NiO水平上。对热压后试件的抛光无包裹体区域做了显微红外光谱分析，红外光束平行于试件的[010]方向，最小直径为15－25μm。分析结果表明，在实验温压范围内，橄榄石中的水含量随铁含量的增加而增加。对于同样铁含量的橄榄石，则随湿度的增加，橄榄石的水含量也增加。最后，本文就实验结果对地幔动力学过程的启示进行了简单的讨论。     

A general method for simulating reactive dissolution in carbonate rocks with arbitrary geometry

The two-scale continuum model is widely used in simulating the reactive dissolution process and predicting the optimum injection rate for carbonate reservoir acidizing treatment. The numerical methods of this model are currently based on structured grids, which are not applicable for complicated geometries. In this study, a general numerical scheme for simulating a reactive flow problem on both structured and unstructured grids is presented based on the finite volume method (FVM). The convection and diffusion terms involved in the reactive flow model are discretized by using the upwind scheme and two-point flux approximation (TPFA), respectively. The location of the centroid node inside each control volume is moved by using an optimization algorithm to make the connections with the surrounding elements as orthogonal as possible, which systematically improves the accuracy of the TPFA scheme. Additionally, in order to avoid the computational complexity resulting from the discretization of the non-linear term, the mass balance equation is only discretized in the spatial domain to get a set of ordinary differential equations (ODEs). These ODEs are coupled with the reaction equations and then solved using the numerical algorithm on ODEs. The accuracy and efficiency of the proposed method are studied by comparing the results obtained from the proposed numerical method with previous experimental and numerical results. This comparison indicates that, compared with the previous methods, the proposed method predicts the wormhole structure more accurately. Finally, the presented method is used to check the effect of the domain geometry, and it is found that the geometry of the flow domain has no effect on the optimum injection velocity, but the radial domain requires a larger breakthrough volume than the linear domain when other parameters are fixed.     

Phase-equilibrium constraints on titanite and rutile activities in mafic epidote amphibolites and geobarometry using titanite–rutile equilibria

Activities of titanite (Ttn, CaTiSiO₅) and/or rutile (Rt, TiO₂) phase components were calculated for 45 well‐characterized natural titanite‐ or rutile‐undersaturated epidote–amphibolites by using the equilibria: (i) 3 anorthite + 2 zoisite/clinozoisite + rutile + quartz = 3 anorthite + titanite + water (referred to as TZARS) and (ii) anorthite + 2 titanite = grossular + 2 rutile + quartz (referred to as GRATiS). In titanite‐bearing and rutile‐absent samples a_Rt is 0.75 ± 0.26. In titanite‐absent, rutile‐bearing samples a_Ttn is 0.89 ± 0.16. Mean values derived for a_Rt/a_Ttn are 0.92 ± 0.12 for rutile + titanite‐bearing samples and 0.42 ± 0.27 for samples lacking both titanite and rutile. Use of these values with TZARS yields pressure estimates for epidote–amphibolites that differ on average by <0.5 kbar from those recorded by established mineral barometers, even where both titanite and rutile are lacking. Despite rather large uncertainties in the average values obtained for a_Rt, a_Ttn or a_Rt/a_Ttn, application of TZARS yields pressure estimates that agree with independent estimates to within ±0.5 kbar for titanite‐ and/or rutile‐saturated samples, and to within ±0.8 kbar for samples that contain neither Ti‐phase. The accuracy and precision of the TZARS barometer are comparable to that of many well‐calibrated barometers. TZARS offers a much‐needed barometer for mafic rocks metamorphosed at epidote‐bearing amphibolite and blueschist facies conditions. In addition, the results provide a basis for application of other thermobarometers, such as Ti‐in‐zircon, where rutile activity is required as input.     

A method for calculating fractional s-character for bonds of tetrahedral oxyanions in crystals

A method for calculating fractional s-character, f _s , for TO bonds has been devised to apply to TO₄ tetrahedral oxyanions in crystals. These f _s -values rank bond lengths with the better correlations obtained for T atoms associated with larger bond strengths and larger electronegativities. As a simple formula, it is found that 2cot²〈?〉₃ does a good job of estimating f _s where 〈?〉₃ is the triple angle average of the three angles common to a given bond.     

Solute carbon and carbon segregation in magnesium oxide single crystals — a secondary ion mass spectrometry study

If carbon is to be analyzed by secondary ion mass spectroscopy (SIMS) in an oxide such as MgO, one has to know how the carbon is incorporated in the oxide host structure, before a successful experiment can be planned. If the carbon impurities derive from dissolved CO₂ component which form a solid solution while the crystal grew from a melt in equilibrium with CO₂, upon cooling, the solid solution becomes supersaturated with respect to the volatile CO₂ component. This creates a thermodynamic driving force for exsolution leading to carbon segregation towards the surface. At the surface rapid degassing occurs in vacuum, enhanced by ion bombardment and electron irradiation. Using freshly cleaved synthetic MgO single crystals it can be shown by SIMS (i) that contamination during short exposure to air and during evacuation remains slight, (ii) that rapid surface/subsurface segregation of solute carbon seems to compete with rapid degassing so that, while no extended segregation profile builds up, the carbon concentration in the bulk beneath the surface decreases to a constant level, (iii) that electron irradiation speeds up degassing, (iv) that heating speeds up carbon diffusion, hence its segregation from the bulk, and (v) that Ar⁺ ion sputtering for the purpose of removing possible contaminants reduces the driving force for carbon surface segregation to the point that no segregation profile can be observed. By placing freshly cleaved MgO crystals under isotopically 99 percent pure ¹³CO₂ for various periods of time subsequent SIMS analysis reveals extended ¹²C profiles, probably about 1 μm wide, which can only have formed by ¹²C segregation from the bulk. These results confirm earlier reports that solute carbon exists as mobile impurity in synthetic MgO and natural olivine, probably due to dissolved CO₂ component.     

Identifying “free” and “bound” lipid fractions in stalagmite samples: An example from Heshang Cave,Southern China

Stalagmites are good archives for paleoecological change, as they are easy to date, and contain multiple environmental proxy records, including climatic records from oxygen isotopes. Lipid biomarkers preserved within stalagmites have recently been used to investigate changes in the overlying soil and vegetation. However, the understanding of lipid records from stalagmites is still at an early stage, and is hindered by the low abundances of lipids preserved and the complexity of the organic matter signal. Here the first results of a sequential extraction procedure are presented, that enables examination of the distribution patterns of “free” (solvent extraction) and “bound” (including physically bound within the calcite matrix and chemically bound to macromolecules) lipids in a stalagmite from southern China. In both groups the dominant compounds are saturated fatty acids, which are an order of magnitude more abundant in the “bound” phase. n-alkanes and n-alcohols chiefly appear in the “free” lipids. In contrast, 3-hydroxy acids are predominantly released under strong acid reflux conditions, suggesting a principal input from bacterial membrane compounds. A direct comparison between the present results and the published data from an Ethiopian stalagmite shows significant differences in the lipid signals from separate sites, with a stronger microbial signal in the Chinese sample. This preliminary investigation of lipid distributions in different modes highlights the importance of microbial geochemical processes in karst systems and supports the use of stalagmites in paleoecological reconstruction.     

Oxygen,carbon and strontium isotope systematics in two profiles across the Glarus thrust: implications for fluid flow

. The Glarus thrust is a prominent tectonic feature in the eastern Helvetic Alps. It has been recognized as a potential major pathway for syntectonic crustal scale fluid flow. The oxygen, carbon and strontium isotope patterns obtained from two vertical profiles across the thrust indicate fundamentally different flow regimes in the southern section of the thrust, where the footwall is represented by Mesozoic limestones, and in the northern section, where the footwall is represented by Tertiary flysch. At the Grauberg locality in the south, the observed isotope patterns give evidence of a net mass transport component from the hanging wall Verrucano to the footwall limestone with a maximum time-integrated volumetric fluid flux of 6.1 m³/m² In the south, the hydration of the lowermost 10 to 20 m of the hanging wall Verrucano requires introduction of an aqueous fluid by subhorizontal flow along the thrust with a minimum time integrated flux of 240 m³/m². At the Lochseite locality in the north, the isotope patterns indicate a vertical mass transport component from the footwall flysch to the hanging wall Verrucano with a time-integrated fluid flux of 2.6 m³/m². In the north, the fluids were probably derived from compaction and dehydration of the footwall flysch during thrusting. The ascending fluids were ponded below the Verrucano and 'lubricated' the thrust. Short-term pressure drops associated with seismic motion along the thrust led to the precipitation of calcite in veins at the thrust surface contributing material to the Lochseiten calc-mylonite, a thin calc-mylonite layer at the thrust contact. Although cross thrust fluid flow may have been two to three orders of magnitude smaller than flow along the thrust, it had a major impact on the isotopic composition of the Lochseiten calc-mylonite. In particular, it buffered the oxygen isotope composition of the calc-mylonite towards the relatively ¹⁸O-depleted composition of the hanging wall Verrucano in the south and towards the relatively ¹⁸O-enriched compositions of the footwall flysch in the north. By this mechanism a regional south to north ¹⁸O-enrichment trend was simulated within the Lochseiten calc-mylonite.     

大位移条件下水平受荷单桩的简明弹塑性计算方法

软土地层中当桩顶水平荷载较大时,采用传统m法计算容易低估桩身弯矩与挠曲变形,有必要针对该问题提出相关计算方法。将地基土体简化为理想弹塑性体,假定桩身某一深度处存在土体的弹塑性变形临界点,临界点以上的土体进入塑性变形状态,而临界点以下的土体仍处于弹性变形状态,分段建立桩身挠曲微分控制方程,得到水平受荷单桩简明弹塑性计算方法。现场单桩实测和参数敏感性分析结果表明:采用简明弹塑性计算方法得到的桩身最大弯矩较传统m法计算精度提高38.1%;桩身最大水平位移计算精度提高22.3%;桩顶边界条件对桩身水平位移与弯矩沿桩身的分布规律影响显著;桩身最大弯矩和水平位移对土体的极限抗力系数及其形状参数较敏感,设计中宜按下限值选取。     

基于附加层法的非规则剖面的瑞利波频散特性研究

对于某些非规则水平层状地质剖面,采用传递矩阵法获得的瑞利波频散曲线可能会出现截止或截断现象,但采用附加层的方法可避免这些现象的发生。通过采用附加层法前后相速度与激发强度大小的对比,来论证附加层法的正确性。对3个典型的非规则剖面进行模拟计算,利用传递矩阵法算得添加附加层前后瑞利波频散曲线,以及某些特征点处位移分量随深度的变化情况。通过对计算结果的对比分析可知,附加层法不仅可计算出地质剖面的实数导波,还可算出其泄漏模式波;在固定的频率点上,附加层的埋深越大,则添加附加层后算得的激发强度大小与原模型的误差越小;当附加层的埋深固定时,添加附加层后计算出的位移曲线与原位移曲线之间的误差随着对比点频率的增大而减小。可见,今后根据实际工作的精度要求,可用附加层法计算某些非规则剖面的瑞利波频散曲线。     

Rectangular parallellepiped resonance method for piezoelectric crystals and elastic constants of alpha-quartz

A theory is presented on the free osculations of a rectangular parallelepiped of piezoelectric crystal, by extending the theory of the rectangular parallelepiped resonance (RPR) method to determine elastic constants of crystals, as exemplified by an alpha-quartz specimen. The piezoelectric contribution to resonance frequencies was examined numerically on the specimen, and it was revealed that piezoelectricity causes approximately 5 kHz increase around 1 MHz. The resonance frequencies of the specimen were measured and inverted to elastic constants by least squares inversion. The inversion was by both the previous non-piezoelectric or elastic theory and by the present piezoelectric theory. The use of the non-piezoelectric theory resulted in an overestimate of 2σ or 0.6% in c ₁₁ and underestimate of σ or 6% in c ₁₂. These are the constants expected to be most affected by piezoelectricity. Errors are less than σ in the other constants. During measurement, it was found that the force applied to hold the specimen caused deviations from free oscillation and experimental errors of 5 kHz. The correction for this force is of some importance in RPR studies of piezoelectric crystals.     

“This is a Montreal Issue”: Negotiating responsibility in global production and investment networks

While geographers have increasingly focused on how global commodity and production networks create new ‘geographies of responsibility’ there has been little empirical work considering how responsibility is worked into management systems and social activism in such networks. Drawing on literature from global production networks, geographies of responsibility and other literatures, this paper explores the dynamic and contested ways in which concepts of responsibility can play a role in network regulation. Both foreign direct investment and commodity networks (here referred to as ‘global production and investment networks’) are subject to complex negotiations and compromises involving corporate social responsibility and sustainability initiatives as well as shareholder activist, human rights, labor, and environmental activism. This is illustrated by reference to conflicts in Canada over Alcan, Inc.’s investments from 1993 to 2007 in the Utkal Alumina Project in Orissa, India. The project involved significant socio-environmental conflict. In Canada, Alcan’s investment was met by civil society campaigns that tested the company’s commitments to sustainability and corporate social responsibility. The case study suggests revising theories of geographies of responsibility. While foreign direct investment can create new relationships between distant others, these are fluid and contingent and not necessarily desirable. Rather than see networks as a source of responsibility we should work to ensure that the relationships that networks foster be structured to ensure our deeper values are respected.     

基于应力增量的单桩p-y曲线分析方法

p-y曲线研究进展一直缓慢,究其原因是缺乏一套成熟的理论体系作支撑。通过考虑桩侧土的实际应力状态,引入Vesic圆孔扩张理论对水平受荷桩p-y曲线进行研究,建立了新的基于应力增量的p-y曲线分析方法。该方法综合考虑了桩侧土体径向应力增量、纵向应力增量效应及桩体深层转动挠曲对水平土抗力的贡献,是一种考虑多种因素的p-y曲线解,可适用于不同桩型条件下的单桩受力分析。算例结果表明:基于应力增量的p-y曲线分析方法在水平受荷桩应用中赋有很好的正确性和准确性,较Reese法、API法和Kim法对大变形问题有更高计算精度。以上研究成果对p-y曲线的发展及深入研究具有重要意义。     

传递矩阵法分析层状地基中桩的扭转变形

研究了扭矩作用下单桩的扭转变形.采用积分变换和传递矩阵方法,求解了成层土在内部环形荷载作用下的基本解;利用此基本解并考虑桩土位移协调条件,提出了层状地基中单桩扭转变形分析的解析方法;并按此理论方法对匀质地基模型进行了数值计算,其结果与已有经典解答相当吻合.