An experimental study of the grain-scale processes of peridotite melting: implications for major and trace element distribution during equilibrium and disequilibrium melting

The grain-scale processes of peridotite melting were examined at 1,340°C and 1.5 GPa using reaction couples formed by juxtaposing pre-synthesized clinopyroxenite against pre-synthesized orthopyroxenite or harzburgite in graphite and platinum-lined molybdenum capsules. Reaction between the clinopyroxene and orthopyroxene-rich aggregates produces a melt-enriched, orthopyroxene-free, olivine + clinopyroxene reactive boundary layer. Major and trace element abundance in clinopyroxene vary systematically across the reactive boundary layer with compositional trends similar to the published clinopyroxene core-to-rim compositional variations in the bulk lherzolite partial melting studies conducted at similar P–T conditions. The growth of the reactive boundary layer takes place at the expense of the orthopyroxenite or harzburgite and is consistent with grain-scale processes that involve dissolution, precipitation, reprecipitation, and diffusive exchange between the interstitial melt and surrounding crystals. An important consequence of dissolution–reprecipitation during crystal-melt interaction is the dramatic decrease in diffusive reequilibration time between coexisting minerals and melt. This effect is especially important for high charged, slow diffusing cations during peridotite melting and melt-rock reaction. Apparent clinopyroxene-melt partition coefficients for REE, Sr, Y, Ti, and Zr, measured from reprecipitated clinopyroxene and coexisting melt in the reactive boundary layer, approach their equilibrium values reported in the literature. Disequilibrium melting models based on volume diffusion in solid limited mechanism are likely to significantly underestimate the rates at which major and trace elements in residual minerals reequilibrate with their surrounding melt. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Assessing soil erosion in Mediterranean karst landscapes of Lebanon using remote sensing and GIS

Maps showing the potential for soil erosion at 1:100,000 scale are produced in a study area within Lebanon that can be used for evaluating erosion of Mediterranean karstic terrain with two different sets of impact factors built into an erosion model. The first set of factors is: soil erodibility, morphology, land cover/use and rainfall erosivity. The second is obtained by the first adding a fifth factor, rock infiltration. High infiltration can reflect high recharge, therefore decreasing the potential of surface runoff and hence the quantity of transported materials. Infiltration is derived as a function of lithology, lineament density, karstification and drainage density, all of which can be easily extracted from satellite imagery. The influence of these factors is assessed by a weight/rate approach sharing similarities between quantitative and qualitative methods and depending on pair-wise comparison matrix.The main outcome was the production of factorial maps and erosion susceptibility maps (scale 1:100,000). Spatial and attribute comparison of erosion maps indicates that the model that includes a measure of rock infiltration better represents erosion potential. Field investigation of rills and gullies shows 87.5% precision of the model with rock infiltration. This is 17.5% greater than the precision of the model without rock infiltration. These results indicate the necessity and importance of integrating information on infiltration of rock outcrops to assess soil erosion in Mediterranean karst landscapes.     

Exurbia from the bottom-up: Confronting empirical challenges to characterizing a complex system

We describe empirical results from a multi-disciplinary project that support modeling complex processes of land-use and land-cover change in exurban parts of Southeastern Michigan. Based on two different conceptual models, one describing the evolution of urban form as a consequence of residential preferences and the other describing land-cover changes in an exurban township as a consequence of residential preferences, local policies, and a diversity of development types, we describe a variety of empirical data collected to support the mechanisms that we encoded in computational agent-based models. We used multiple methods, including social surveys, remote sensing, and statistical analysis of spatial data, to collect data that could be used to validate the structure of our models, calibrate their specific parameters, and evaluate their output. The data were used to investigate this system in the context of several themes from complexity science, including have (a) macro-level patterns; (b) autonomous decision making entities (i.e., agents); (c) heterogeneity among those entities; (d) social and spatial interactions that operate across multiple scales and (e) nonlinear feedback mechanisms. The results point to the importance of collecting data on agents and their interactions when producing agent-based models, the general validity of our conceptual models, and some changes that we needed to make to these models following data analysis. The calibrated models have been and are being used to evaluate landscape dynamics and the effects of various policy interventions on urban land-cover patterns.     

Magma Ascent and Crustal Accretion at Ultraslow-Spreading Ridges: Constraints from Plagioclase Ultraphyric Basalts from the Arctic Mid-Ocean Ridge

Plagioclase ultraphyric basalts (PUBs) with up to 54% plagioclasephenocrysts were dredged in the rift valley and adjacent flanksof the ultraslow-spreading Mohns and Knipovich ridges. The PUBsshow large variations in crystal morphologies and zoning. Thelarge variations suggest that single basalt samples containa mixture of plagioclase crystals that aggregated at differentlevels in the magma conduits. Resorbed crystals and repeatedreverse zones suggest that the magma reservoirs were replenishedand heated several times. Thin concentric zones with melt inclusions,and sharp reductions in the anorthite content of 3–7%,are common between the reverse zones. These zones, and skeletalcrystals with distinctly lower anorthite contents than massivecrystals, are interpreted to be the result of rapid crystalliztionduring strong undercooling. The changes between short periodsof cooling and longer periods with reheating are explained bymultiple advances of crystal-rich magma into cool regions followedby longer periods of gradual magma inflow and temperature increase.The porphyritic basalts are characterizd by more depleted andmore fractionated compositions than the aphyric basalts, withlower (La/Sm)_N, K₂O and Mg-numbers. This relationship, and theobservation that PUBs are sampled only close to segment centresalong these ridges, suggests that the PUBs formed by higherdegrees of melting and evolved in more long-lived magma reservoirs.We propose that the zoning patterns of plagioclase crystalsand crystal morphologies of these PUBs reflect the developmentand flow of magma through a stacked sill complex-like conduitsystem, whereas the aphyric equivalents represent later flowof magma through the conduit. The formation of voluminous higher-degreemelts may trigger the development of the magma conduits andexplain the generally depleted compositions of PUB magmas. KEY WORDS: basalt; mineral chemistry; MORB; magma mixing; magma chamber; major element     

攀西地区铂族元素地球化学异常分布及其筛选

攀西地区不同时代地层的水系沉积物中铂族元素地球化学背景值具有明显的规律性分布.笔者通过对该地区铂族元素含量在各地层中的分布以及地球化学异常的圈定、归并和分类,结合有关铂族矿床地球化学特征的认识,通过对水系沉积物中Pt、Pd综合异常的筛选,在铂族元素综合衬值异常图的基础上共划分出Ⅰ类异常4个、Ⅱ类异常19个、Ⅲ类异常23个以及Ⅳ类异常14个.初步筛选出的19个Ⅱ类异常子区可为进一步开展铂族矿床找矿工作提供依据.     

金川镍铜矿床数学模型对深部矿化变化趋势的指示

按矿化特征将金川X矿区主矿体分为致密块状特富矿、海绵状富矿和浸染状贫矿.应用MicroMine软件对3类矿化体分别建立了数学模型,并应用距离平方反比法对模型单元块进行了镍、铜品位估值.通过分析和对比矿床数学模型900m至1300m不同水平剖面上的Ni、Cu品位及Ni/Cu比值的变化,得出在X矿区的中部西侧是一个Ni和Cu的矿化中心,并且向深部Cu的矿化强度高于Ni.     

Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Detection of free vertical convection and double-diffusion in groundwater monitoring wells with geophysical borehole measurements

Two algorithms for in-situ detection and identification of vertical free convective and double-diffusive flows in groundwater monitoring wells or boreholes are proposed. With one algorithm the causes (driving forces) and with the other one the effects (convection or double-diffusion) of vertical transport processes can be detected based on geophysical borehole measurements in the water column. Five density-driven flow processes are identified: thermal, solutal, and thermosolutal convection leading to an equalization, as well as saltfingers and diffusive layering leading to an intensification of a vertical density gradient. The occurrence of density-driven transport processes could be proven in many groundwater monitoring wells and boreholes; especially shallow sections of boreholes or groundwater monitoring wells are affected dramatically by such vertical flows. Deep sections are also impaired as the critical threshold for the onset of a density-driven flow is considerably low. In monitoring wells or boreholes, several sections with different types of density-driven vertical flows may exist at the same time. Results from experimental investigations in a medium-scale testing facility with high aspect ratio (height/radius = 19) and from numerical modeling of a water column agree well with paramters of in-situ detected convection cells.