Trace-element composition and zoning in clinopyroxene- and amphibole-group minerals: Implications for element partitioning and evolution of carbonatites

The present work is a first comprehensive study of the trace-element composition and zoning in clinopyroxene- and amphibole-group minerals from carbonatites, incorporating samples from 14 localities worldwide (Afrikanda, Aley, Alnö, Blue River, Eden Lake, Huayangchuan, Murun, Oka, Ozernaya Varaka, Ozernyi, Paint Lake, Pinghe, Prairie Lake, Turiy Mys). The new electron-microprobe data presented here significantly extend the known compositional range of clinopyroxenes and amphiboles from carbonatites. These data confirm that calcic and sodic clinopyroxenes from carbonatites are not separated by a compositional gap, instead forming an arcuate trend from nearly pure diopside through intermediate aegirine–augite compositions confined to a limited range of CaFeSi₂O₆ contents (15–45 mol%) to aegirine with < 25 mol% of CaMgSi₂O₆ and a negligible proportion of CaFeSi₂O₆. A large set of LA-ICPMS data shows that the clinopyroxenes of different composition are characterized by relatively low levels of Cr, Co and Ni (≤ 40 ppm) and manifold variations in the concentration of trivalent lithophile and some incompatible elements (1–150 ppm Sc, 26–6870 ppm V, 5–550 ppm Sr, 90–2360 ppm Zr, and nil to 150 ppm REE), recorded in some cases within a single crystal. The relative contribution of clinopyroxenes to the whole-rock Rb, Nb, Ta, Th and U budget is negligible. The major-element compositional range of amphiboles spans from alkali- and Al-poor members (tremolite) to Na–Al-rich Mg- or, less commonly, Fe-dominant members (magnesiohastingsite, hastingsite and pargasite), to calcic–sodic, sodic and potassic–sodic compositions intermediate between magnesio-ferrikatophorite, richterite, magnesioriebeckite, ferri-nyböite and (potassic-)magnesio-arfvedsonite. In comparison with the clinopyroxenes, the amphiboles contain similar levels of tetravalent high-field-strength elements (Ti, Zr and Hf) and compatible transition elements (Cr, Co and Ni), but are capable of incorporating much higher concentrations of Sc and incompatible elements (up to 500 ppm Sc, 43 ppm Rb, 1470 ppm Sr, 1230 ppm Ba, 80 ppm Pb, 1070 ppm REE, 140 ppm Y, and 180 ppm Nb). In some carbonatites, amphiboles contribute as much as 25% of the Zr + Hf, 15% of the Sr and 35% of the Rb + Ba whole-rock budget. Both clinopyroxenes and amphiboles may also host a significant share (~ 10%) of the bulk heavy-REE content. Our trace-element data show that the partitioning of REE between clinopyroxene (and, in some samples, amphibole) and the melt is clearly bimodal and requires a revision of the existing models assuming single-site REE partitioning. Clinopyroxenes and amphiboles from carbonatites exhibit a diversity of zoning patterns that cannot be explained exclusively on the basis of crystal chemistry and relative compatibility of different trace-element in these minerals. Paragenetic analysis indicates that in most cases, the observed zoning patterns develop in response to removal of selected trace elements by phases co-precipitating with clinopyroxene and amphibole (especially magnetite, fluorapatite, phlogopite and pyrochlore). With the exception of magnesiohastingsite–richterite sample from Afrikanda, the invariability of trace-element ratios in the majority of zoned clinopyroxene and amphibole crystals implies that fluids are not involved in the development of zoning in these minerals. The implications of the new trace-element data for mineral exploration targeting REE, Nb and other types of carbonatite-hosted rare-metal mineralization are discussed.     

Review of vibroseis data acquisition and processing for better amplitudes: adjusting the sweep and deconvolving for the time-derivative of the true groundforce

The goal of vibroseis data acquisition and processing is to produce seismic reflection data with a known spatially-invariant wavelet, preferably zero phase, such that any variations in the data can be attributed to variations in geology. In current practice the vibrator control system is designed to make the estimated groundforce equal to the sweep and the resulting particle velocity data are cross-correlated with the sweep. Since the downgoing far-field particle velocity signal is proportional to the time-derivative of the groundforce, it makes more sense to cross-correlate with the time-derivative of the sweep. It also follows that the ideal amplitude spectrum of the groundforce should be inversely proportional to frequency. Because of non-linearities in the vibrator, bending of the baseplate and variable coupling of the baseplate to the ground, the true groundforce is not equal to the pre-determined sweep and varies not only from vibrator point to vibrator point but also from sweep to sweep at each vibrator point. To achieve the goal of a spatially-invariant wavelet, these variations should be removed by signature deconvolution, converting the wavelet to a much shorter zero-phase wavelet but with the same bandwidth and signal-to-noise ratio as the original data. This can be done only if the true groundforce is known. The principle may be applied to an array of vibrators by employing pulse coding techniques and separating responses to individual vibrators in the frequency domain. Various approaches to improve the estimate of the true groundforce have been proposed or are under development; current methods are at best approximate.     

Single-crystal elasticity of stishovite: New experimental data obtained using high-frequency resonant ultrasound spectroscopy and a Gingham check structure model

Single-crystal elasticity of stishovite was examined using a new experimental technique and an empirical macroscopic model. Employing high-frequency resonant ultrasound spectroscopy, single-crystal elastic constants of stishovite were determined: C₁₁ = 443(3), C₃₃ = 781(4), C₁₂ = 193(2), C₂₃ = 199(2), C₄₄ = 256(2), and C₃₃ = 316(2) GPa. The frequency range of the resonant ultrasound spectroscopy techniques was 6–20 MHz, which is much lower than the ～10 GHz range of the Brillouin scattering technique. Of the elastic constants, the shear elastic constants C₄₄ and C₆₆ are consistent with the average value of the previously mentioned Brillouin scattering. Conversely, the four elastic constants, C₁₁, C₃₃, C₁₂, and C₂₃, slightly deviate outside the range of previous Brillouin scattering results. The present results, except those for C₁₂, are consistent with recent lattice dynamic analysis of inelastic X-ray scattering data. The adiabatic bulk modulus was calculated as 298 GPa, which is smaller and more consistent with the result of compression experiments than any other Brillouin scattering results (301–312 GPa). The present result shows greater P-wave velocity anisotropy (24.7%) than any preceding work. To understand the unique elastic properties of stishovite, the Gingham check model was proposed and examined. The result shows that the octahedron of 6-coordinated Si in stishovite crystal has stiffness comparable to that of diamond.     

A model for the estimation of water filtration and nutrient regeneration by exposed sandy beaches

Exposed sandy beaches filter huge volumes of sea water daily as a result of wave and tide action. A regression model is developed which predicts the volume of water filtered daily by a beach as a function of tide range, beach slope and sand particle size. From these parameters estimates are also made of intertidal distance, the degree of expansion of the intertidal above the still water tide range by wave action and the average path length through the sand interstices percolated by filtered sea water. An equation is derived relating the degree of mineralisation of organic matter in the filtered sea water to the distance filtered through the sand and the sand particle size. Employing all this information, a simple model is used to estimate the volumes of sea water filtered and the amounts of inorganic nitrogen regenerated by intertidal sandy beaches over a range of physical conditions. The implications of this nutrient generation for the inshore zone are discussed.     

Single-crystal EPR study of three radiation-induced defects (Al–O2 3?, Ti3+ and W5+) in stishovite

Single-crystal electron paramagnetic resonance spectra of electron-irradiated stishovite, measured at temperatures from 3.5 to 294?K, reveal three S?=?1/2 radiation-induced defects: an aluminum-associated oxygen hole center and two nd ¹ centers (Ti³⁺ and W⁵⁺). The aluminum-associated oxygen hole center, characterized by an orthorhombic site symmetry, coaxial matrices of the electronic Zeeman g, nuclear hyperfine A(²⁷Al) and nuclear quadrupole P(²⁷Al), and the orientation of the g-minimum axis along an O–O direction and those of the unique A(²⁷Al) and P(²⁷Al) axes perpendicular to the O–O direction, is an Al–O₂ ^3? center, with the unpaired electron equally distributed on two equatorial oxygen atoms of a substitutional Al³⁺ ion at the octahedral Si site. Fully optimized Al-doped structure, theoretical ²⁷Al nuclear hyperfine and quadrupole coupling constants and directions, and 3D spin densities from periodic hybrid density functional theory calculations provide further support for this structural model. Spin Hamiltonian parameters of the Ti³⁺ and W⁵⁺ centers, which are confirmed by their diagnostic ⁴⁷Ti, ⁴⁹Ti and ¹⁸³W hyperfine structures, arise from electron trapping on substitutional Ti⁴⁺ and W⁶⁺ ions at the octahedral Si site.     

Predicting catchment sediment yield in Mediterranean environments: the importance of sediment sources and connectivity in Italian drainage basins

Predicting sediment yield at the catchment scale is one of the main challenges in geomorphologic research. The application of both physics‐based models and regression models has until now not provided very satisfying results for prediction of sediment yield for medium to large sized catchments (c. >50 km²). The explanation for this lies in a combination of the large data requirements of most models and a lack of knowledge to describe all processes and process interactions at the catchment scale. In particular, point sources of sediment (e.g. gullies, mass movements), connectivity and sediment transport remain difficult to describe in most models. From reservoir sedimentation data of 44 Italian catchments, it appeared that there was a (non‐significant) positive relation between catchment area and sediment yield. This is in contrast to what is generally expected from the theory of decreasing sediment delivery rates with increasing catchment area. Furthermore, this positive relation suggests that processes other than upland erosion are responsible for catchment sediment yield. Here we explore the potential of the Factorial Scoring Model (FSM) and the Pacific Southwest Interagency Committee (PSIAC) model to predict sediment yield, and indicate the most important sediment sources. In these models different factors are used to characterize a drainage basin in terms of sensitivity to erosion and connectivity. In both models an index is calculated that is related to sediment yield. The FSM explained between 36 and 61 per cent of the variation in sediment yield, and the PSIAC model between 57 and 62 per cent, depending on the factors used to characterize the catchments. The FSM model performed best based on a factor to describe gullies, lithology, landslides, catchment shape and vegetation. Topography and catchment area did not explain additional variance. In particular, the addition of the landslide factor resulted in a significantly increased model performance. The FSM and PSIAC model both performed better than a spatially distributed model describing water erosion and sediment transport, which was applied to the same catchments but explained only between 20 and 51 per cent of the variation in sediment yield. Model results confirmed the hypothesis that processes other than upland erosion are probably responsible for sediment yield in the Italian catchments. A promising future development of the models is by the use of detailed spatially distributed data to determine the scores, decrease model subjectivity and provide spatially distributed output. Copyright © 2006 John Wiley & Sons, Ltd.     

The spatio-temporal development of Copenhagen’s bicycle infrastructure 1912–2013

Cycling plays an important role in low-carbon transitions. Around the globe, cities are constructing bicycle infrastructure. The city of Copenhagen has a bicycle-friendly infrastructure celebrated for its fine-meshed network. This study documents the spatio-temporal development of Copenhagen’s bicycle infrastructure and explores how the development corresponds to other processes of urban transformation. The study builds on historical maps of bicycle infrastructure that are digitised into geographical information, which allows for a comprehensive analysis of the formation of the network. In search for identifying drivers, the study analyses the city’s spatial growth pattern, migration pattern, development of road network and changes in the transport culture. Analyses reveal that the bicycle infrastructure expanded at a relatively constant pace during distinct periods of urban transformation, including periods when the city suffered from spatial, economic and demographic decline, and dominance of car traffic. By discussing reasons and demands for constructing bicycle infrastructure, the study identifies four distinct periods in which bicycle infrastructure was constructed to enhance comfort and safety (first cycling city); the flow for cars (car city); urban liveability for soft transport (liveable city); and, finally, to improve the flow for cyclists as part a strategic re-design of urban space (liveable cycling city).     

微波消解法在水泥工业样品分析中的应用

在过去的几年中,辅助燃料(如旧轮胎和塑料垃圾)在水泥行业应用越来越重要,不仅是价格原因,更重要的是这些应用在经济和环境方面对于处理废弃材料起到了重要作用。水泥回转窑的高温是塑料垃圾焚烧的最好条件。旧轮胎作为辅助燃料的应用已经很长时间了,因为它们的材料、反应行为及对最终产品的影响都已经得到了很好的认知,而塑料垃圾的分析相对比较困难。如塑料瓶子及塑料薄膜往往添加各种不同的色素及添加剂。水泥对锌和其他痕量元素非常敏感,尤其是氯元素。因此,在燃烧前,必须先将PVC材料挑出,而食品包装材料必须先经过清洗避免氯化钠残留而带入氯元素。这就是粉碎垃圾需要完全消解及分析的关键原因。因为垃圾的组成是不均一的,完全     

Risk perception and decision-making: do farmers consider risks from climate change?

This paper provides one of the first empirical studies that examine the impact of climate change adaptation practices on technical efficiency (TE) among smallholder farmers in Nepal. An adaptation index is used to explore the impact of farmers’ adaptation on TE using the stochastic frontier analysis framework. Data for six districts of Nepal representing all three agro-ecological regions (terai, hill, and mountain) were collected from a focus group discussion, a stakeholder workshop and a household survey. The survey shows that about 91% of the farming households have adopted at least one practice to minimize the adverse impacts of climate change. Empirical results reveal that adaptation is an important factor explaining efficiency differentials among farming households. Those adopting a greater number of adaptation practices on a larger scale are, on average, found to be 13% more technically efficient than those adopting fewer practices on smaller scale. The empirical results also show that average TE is only 0.72, indicating that there are opportunities for farming households in Nepal to further improve productive efficiency, on average by 28%. Other important factors that explain variations in the productive efficiency across farming households include farmer’s education level, irrigation facilities, market access, and social capital such as farmer’s participations in relevant agricultural organizations and clubs. This study provides empirical evidence to policy makers that small scale adjustments made by farmers in response to climate change impacts are effective in improving farmers’ efficiency in agriculture production. This indicates a need for farmers’ involvement in climate change adaptation planning.     

Cache-efficient parallel eikonal solver for multicore CPUs

Numerical solution of the eikonal equation is frequently used to compute first-arrival travel times for a given velocity model in seismic applications. Computations for large three-dimensional models become expensive requiring the use of efficient parallel solvers. We present new parallel implementations of the fast sweeping and locking sweeping methods optimized for shared memory systems such as multicore CPUs; we call them block fast sweeping method (BFSM) and block locking sweeping method (BLSM). Proposed methods are based on the domain decomposition approach with a special attention paid to high efficiency of the cache utilization and task execution synchronization. Performance tests on realistic models show high parallel efficiency of 85–95% on modern multicore CPUs and require the same number of iterations to converge as do the serial sweeping methods. We also highlight the importance of properly selecting the stopping criterion in the iterative sweeping methods aiming for a balance between computational time and accuracy of the result required by an application. In particular, we show that in seismic applications one can reach reasonable accuracy of computed travel times while dramatically reducing the number of iterations compared to the case of using the full convergence stopping criterion.