Monte Carlo simulation of mixing in Ca3Fe2Ge3O12–Ca4Ge4O12 garnets and implications for the thermodynamic stability of pyrope–majorite solid solution

Static lattice energy calculations, based on empirical pair potentials, were performed for a large set of structures differing in the arrangement of octahedral cations within the garnet 2 × 2 × 2 supercell. The compositions of these structures varied between Ca₃Fe₂Ge₃O₁₂ and Ca₄Ge₄O₁₂. The energies were cluster expanded using pair and quaternary terms. The derived ordering constants were used to constrain Monte Carlo simulations of temperature-dependent mixing properties in the ranges of 1,073–3,673 K and 0–10 GPa. The free energies of mixing were calculated using the method of thermodynamic integration. The calculations predict a wide miscibility gap between Fe-rich (cubic) and Fe-pure (tetragonal) garnets consistent with recent experimental observations of Iezzi et al. (Phys Chem Miner 32:197–207, 2005). It is shown that the miscibility gap arises due to a very strong cation ordering at the Fe-pure composition, driven by the charge difference between Ca²⁺ and Ge⁴⁺ cations. The structural and thermodynamic analogies between Ca–Ge and Mg–Si systems suggest that a similar miscibility gap should exist between pyrope and Mg–Si-majorite.     

Geogenic and anthropogenic zinc in epiphytic and terricolous lichens in Finland

At 31 sites in Finland the Zn content was determined of the epiphytic lichen Hypogymnia physodes (L.) Nyl. growing on pine and birch, the terricolous lichen Cladina spp. in pine bark, and the inner part of birch bark. The anomalous contents of Zn in H. physodes on pine and Cladina spp. in soil appeared to be strongly related to sulphide ore deposits. Anthropogenic Zn correlated closely with S and Fe in H. physodes on pine and with Fe in Cladina spp. The binding of geogenic Zn was related to iron hydroxides, and was attributed to the different uptake mechanism of S in these lichens. Elevated values of Zn in H. physodes on birch were due to both geogenic and anthropogenic sources. Zinc in pine bark increased with the Zn content in H. physodes and reached a limiting value.     

Comparison of the Ordovician-Carboniferous Boundary Between Korea and NE China: Implications for Correlation and Tectonic Evolution

There is a great hiatus between Ordovician and Carboniferous strata in the Northeast China and Korean Peninsula. In order to understand geology and tectonic evolution, and to find out the similarities and differences in both regions, two sections in the Western Hill near Beijing in NE China and several sections in the Korean Peninsula were selected to examine their geologic boundaries between Lower and Upper Paleozoic strata to compare their characteristic features. At four sites in the two sections in the Western Hill near Beijing were examined their contact relations. The Hui Yu section is the same horizon where one site is top of a quarry hill and the other of down hill. Mid-Carboniferous Qingshuijian Formation rests on the Ordovician Majiagou Formation. Limestone beds are more commonly intercalated with shale and sandstone at site 2 of the Hui Yu section, while at site 1, conglomerate beds are dominant. Site 1 of the Se Shu Fen section shows eroded and concealed karst topography and conglomerate beds are intercalated within shale beds. Silurian and Devonian strata are absent in these areas. In the Korean Peninsula, most O-C contacts occur between Ordovician limestone formation and Carboniferous strata, although Silurian strata occur beneath the Carboniferous strata in the Jeongseon area and Pyeongnam Basin. Most contact relations are parallel unconformity and angular unconformity is rarely seen. The O-C relations in both regions are similar to each other, and these indicate that the Korean Peninsula was located near or belonged to the Sino-Korean paraplatform during Paleozoic time.     

A new bounding surface model for thermal cyclic behaviour

To accurately predict soil volume changes under thermal cycles is of great importance for analysing the performance of many earth structures such as the energy pile and energy storage system. Most of the existing thermo‐mechanical models focus on soil behaviour under monotonic thermal loading only, and they are not able to capture soil volume changes under thermal cycles. In this study, a constitutive model is proposed to simulate volume changes of saturated soil subjected to cyclic heating and cooling. Two surfaces are defined and used: a bounding surface and a memory surface. The bounding surface and memory surface are mainly controlled by the preconsolidation pressure (a function of plastic volumetric strain) and the maximum stress experienced by the soil, respectively. Under thermal cycles, the distance of the two surfaces and plastic modulus increase with an accumulation of plastic strain. By adopting the double surface concept, a new elastoplastic model is derived from an existing single bounding surface thermo‐mechanical model. Comparisons between model predictions and experimental results reveal that the proposed model is able to capture soil volume changes under thermal cycles well. The plastic strain accumulates under thermal cycles, but at a decreasing rate, until stabilization. Copyright © 2017 John Wiley & Sons, Ltd.     

Spatialising the Melanesian Canarium industry: Understanding economic upgrading in an emerging industry among three Pacific small island states

This paper examines economic upgrading in the Canarium indicum (Canarium) nut industry in Vanuatu, the Solomon Islands and Papua New Guinea. Canarium is a tree that is indigenous to Melanesia, and has been the subject of several commercialisation attempts since 1988. The paper assesses the outcomes to various actors in the Canarium supply chain from attempts to upgrade industry products and processes by: (1) increasing the available resource in suitable locations; (2) improving nut products and processing techniques; (3) increasing actors’ knowledge and supply chain coordination; and (4) establishing product standards. A two-phase data generation process was implemented. Document analysis and participant observations of industry workshops initiated a set of four economic upgrading interventions that were adaptively implemented. A second stage of workshops and 76 interviews enabled outcomes to be assessed at the project’s end. Findings suggest that a small number of urban-based entrepreneurs benefit and subsequently are lead actors in industry development, but at the expense of benefits being distributed to a larger, more spatially disparate group of smallholder and small commercial growers. These economic upgrading outcomes are circumscribed by core-periphery relations in Pacific small island states and the scale of industry in each country. Thus, spatial inequalities are reproduced through the emerging Canarium industries. We argue that different routes to industry development are required in each country. Development initiatives that capitalise on the benefits of micro-enterprise clusters, joint action and regional institutional arrangements are proposed to overcome the impediments imposed by the particular geographies of Pacific island states.     

Recent Advances in Prediction and Processing of Strong Ground Motions

We present an overview of our recent results on utilizing small earthquakes in the earthquake engineering practice. Site-specific ground motion time-histories of large earthquakes can be successfully simulated using recordings of small earthquakes which are often referred to as 'empirical Green's functions' in seismology. Another important practical problem is whether and how these observations can be used in seismic risk studies which are based on empirical attenuation relations for ground motion parameters. We study a possibility of extrapolating attenuation relations for small earthquakes, to larger magnitudes using the data from the Garner Valley downhole array in Southern California. Finally we introduce efficient ground motion processing techniques in frequency- and time-domains and apply them to site response estimation.     

Formulation of a unified constitutive model for clays and sands

This paper presents a new generalized effective stress model, referred to as MIT-S1, which is capable of predicting the rate independent, effective stress–strain–strength behaviour of uncemented soils over a wide range of confining pressures and densities. Freshly deposited sand specimens compressed from different initial formation densities approach a unique condition at high stress levels, referred to as the limiting compression curve (LCC), which is linear in a double logarithmic void ratio, e, mean effective stress space, p′. The model describes irrecoverable, plastic strains which develop throughout first loading using a simple four-parameter elasto-plastic model. The shear stiffness and strength properties of sands in the LCC regime can be normalized by the effective confining pressure and hence can be unified qualitatively, with the well-known behaviour of clays that are normally consolidated from a slurry condition along the virgin consolidation line (VCL). At lower confining pressures, the model characterizes the effects of formation density and fabric on the shear behaviour of sands through a number of key features: (a) void ratio is treated as a separate state variable in the incrementally linearized elasto-plastic formulation: (b) kinematic hardening describing the evolution of anisotropic stress–strain properties: (c) an aperture hardening function controls dilation as a function of ‘formation density’; and (d) the use of a single lemniscate-shaped yield surface with non-associated flow. These features enable the model to describe characteristic transitions from dilative to contractive shear response of sands as the confining pressure increases. This paper summarizes the procedures used to select input parameters for clays and sands, while a companion paper compares model predictions with measured data to illustrate the model capability for describing the shear behaviour of clays and sands. Copyright © 1999 John Wiley & Sons, Ltd.     

Quantification of karst aquifer discharge components during storm events through end-member mixing analysis using natural chemistry and stable isotopes as tracers

Karst aquifer components that contribute to the discharge of a water supply well in the Classical Karst (Kras) region (Italy/Slovenia) were quantitatively estimated during storm events. Results show that water released from storage within the epikarst may comprise as much as two-thirds of conduit flow in a karst aquifer following rainfall. Principal components analysis (PCA) and end-member mixing analysis (EMMA) were performed using major ion chemistry and the stable isotopes of water (δ¹⁸O, δ²H) and of dissolved inorganic carbon (δ¹³C_DIC) to estimate mixing proportions among three sources: (1) allogenic river recharge, (2) autogenic recharge, and (3) an anthropogenic component stored within the epikarst. The sinking river most influences the chemical composition of the water-supply well under low-flow conditions; however, this proportion changes rapidly during recharge events. Autogenic recharge water, released from shallow storage in the epikarst, displaces the river water and is observed at the well within hours after the onset of precipitation. The autogenic recharge end member is the second largest component of the well chemistry, and its contribution increases with higher flow. An anthropogenic component derived from epikarstic storage also impacts the well under conditions of elevated hydraulic head, accounting for the majority of the chemical response at the well during the wettest conditions.     

Calculated Phase Relations in the System NCKFMASH (Na2O-CaO-K2O-FeO-MgO-Al2O3-SiO2-H2O) for High-Pressure Metapelites

Petrogenetic grids in the system NCKFMASH (Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O)and the subsystems NCKMASH and NCKFASH calculated with the softwareTHERMOCALC 3.1 are presented for the P–T range 7–30kbar and 450–680°C, for assemblages involving garnet,chloritoid, biotite, carpholite, talc, chlorite, kyanite, staurolite,paragonite, glaucophane, jadeite, omphacite, diopsidic pyroxene,plagioclase, zoisite and lawsonite, with phengite, quartz/coesiteand H₂O in excess. These grids, together with calculated compatibilitydiagrams and P–T and T–X_Ca and P–X_Ca pseudosectionsfor different bulk-rock compositions, show that incorporationof Ca into the NKFMASH system leads to many of the NKFMASH invariantequilibria moving to lower pressure and/or lower temperature,which results, in most cases, in the stability of jadeite andgarnet being enlarged, but in the reduction of stability ofglaucophane, plagioclase and AFM phases. The effect of Ca onthe stability of paragonite is dependent on mineral assemblageat different P–T conditions. The calculated NCKFMASH diagramsare powerful in delineating the phase equilibria and P–Tconditions of natural pelitic assemblages. Moreover, contoursof the calculated phengite Si isopleths in P–T and P–X_Capseudosections confirm that phengite barometry in NCKFMASH isstrongly dependent on mineral assemblage. KEY WORDS: phase relations; metapelites; NCKFMASH; THERMOCALC; phengite geobarometry     

基于物理特征的平面流场拓扑简化算法

基于网格合并思想实现的平面流场拓扑简化,存在物理特征被忽略和并行实现困难的问题,难以满足远程可视化的要求。文中提出1种基于物理特征的平面流场拓扑简化算法。算法通过计算生成平面流场的物理特征参数图像,利用图像分割提取平面流场的特征区域达到流场简化的目的。实验结果表明,算法能从物理特征和宏观的角度出发快速、准确描述平面流场的拓扑结构,为平面流场可视化提供了1种新方法。