The Kabanga Ni sulfide deposit, Tanzania: I. Geology, petrography, silicate rock geochemistry, and sulfur and oxygen isotopes

The Kabanga Ni sulfide deposit represents one of the most significant Ni sulfide discoveries of the last two decades, with current indicated mineral resources of 23.23 Mt at 2.64% Ni and inferred mineral resources of 28.5 Mt at 2.7% Ni (Nov. 2008). The sulfides are hosted by a suite of ∼1.4 Ga ultramafic–mafic, sill-like, and chonolithic intrusions that form part of the approximately 500 km long Kabanga–Musongati–Kapalagulu igneous belt in Tanzania and Burundi. The igneous bodies are up to about 1 km thick and 4 km long. They crystallized from several compositionally distinct magma pulses emplaced into sulfide-bearing pelitic schists. The first magma was a siliceous high-magnesium basalt (approximately 13.3% MgO) that formed a network of fine-grained acicular-textured gabbronoritic and orthopyroxenitic sills (Mg# opx 78–88, An plag 45–88). The magma was highly enriched in incompatible trace elements (LILE, LREE) and had pronounced negative Nb and Ta anomalies and heavy O isotopic signatures (δ¹⁸O +6 to +8). These compositional features are consistent with about 20% contamination of primitive picrite with the sulfidic pelitic schists. Subsequent magma pulses were more magnesian (approximately 14–15% MgO) and less contaminated (e.g., δ¹⁸O +5.1 to +6.6). They injected into the earlier sills, resulting in the formation of medium-grained harzburgites, olivine orthopyroxenites and orthopyroxenites (Fo 83–89, Mg# _opx 86–89), and magmatic breccias consisting of gabbronorite–orthopyroxenite fragments within an olivine-rich matrix. All intrusions in the Kabanga area contain abundant sulfides (pyrrhotite, pentlandite, and minor chalcopyrite and pyrite). In the lower portions and the immediate footwall of two of the intrusions, namely Kabanga North and Kabanga Main, there occur numerous layers, lenses, and veins of massive Ni sulfides reaching a thickness of several meters. The largest amount of high grade, massive sulfide occurs in the smallest intrusion (Kabanga North). The sulfides have heavy S isotopic signatures (δ³⁴S wr = +10 to +24) that broadly overlap with those of the country rock sulfides, consistent with significant assimilation of external sulfur from the Karagwe–Ankolean sedimentary sequence. However, based partly on the relatively homogenous distribution of disseminated sulfides in many of the intrusive rocks, we propose that the Kabanga magmas reached sulfide saturation prior to final emplacement, in staging chambers or feeder conduits, followed by entrainment of the sulfides during continued magma ascent. Oxygen isotope data indicate that the mode of sulfide assimilation changed with time. The heavy δ¹⁸O ratios of the early magmas are consistent with ingestion of the sedimentary country rocks in bulk. The relatively light δ¹⁸O ratios of the later magmas indicate less bulk assimilation of the country rocks, but in addition the magmas selectively assimilated additional S, possibly through devolatization of the country rocks or through cannibalization of magmatic sulfides deposited in the conduits by preceding magma surges. The intrusions were tilted at ca. 1.37 Ga, during the Kibaran orogeny and associated synkinematic granite plutonism. This caused solid-state mobilization of ductile sulfides into shear zones, notably along the base of the intrusions where sulfide-hornfels breccias and lenses and layers of massive sulfides may reach a thickness of >10 m and can extend for several 10 s to >100 m away from the intrusions. These horizons represent an important exploration target for additional nickel sulfide deposits.     

Why GPS makes distances bigger than they are

Global navigation satellite systems such as the Global Positioning System (GPS) is one of the most important sensors for movement analysis. GPS is widely used to record the trajectories of vehicles, animals and human beings. However, all GPS movement data are affected by both measurement and interpolation errors. In this article we show that measurement error causes a systematic bias in distances recorded with a GPS; the distance between two points recorded with a GPS is – on average – bigger than the true distance between these points. This systematic ‘overestimation of distance’ becomes relevant if the influence of interpolation error can be neglected, which in practice is the case for movement sampled at high frequencies. We provide a mathematical explanation of this phenomenon and illustrate that it functionally depends on the autocorrelation of GPS measurement error (C). We argue that C can be interpreted as a quality measure for movement data recorded with a GPS. If there is a strong autocorrelation between any two consecutive position estimates, they have very similar error. This error cancels out when average speed, distance or direction is calculated along the trajectory. Based on our theoretical findings we introduce a novel approach to determine C in real-world GPS movement data sampled at high frequencies. We apply our approach to pedestrian trajectories and car trajectories. We found that the measurement error in the data was strongly spatially and temporally autocorrelated and give a quality estimate of the data. Most importantly, our findings are not limited to GPS alone. The systematic bias and its implications are bound to occur in any movement data collected with absolute positioning if interpolation error can be neglected.     

Turbulent flux observations and modelling over a shallow lake and a wet grassland in the Nam Co basin,Tibetan Plateau

The Tibetan Plateau plays an important role in the global water cycle and is strongly influenced by climate change. While energy and matter fluxes have been more intensely studied over land surfaces, a large proportion of lakes have either been neglected or parameterised with simple bulk approaches. Therefore, turbulent fluxes were measured over wet grassland and a shallow lake with a single eddy-covariance complex at the shoreline in the Nam Co basin in summer 2009. Footprint analysis was used to split observations according to the underlying surface, and two sophisticated surface models were utilised to derive gap-free time series. Results were then compared with observations and simulations from a nearby eddy-covariance station over dry grassland, yielding pronounced differences. Observations and footprint integrated simulations compared well, even for situations with flux contributions including grassland and lake. The accessibility problem for EC measurements on lakes can be overcome by combining standard meteorological measurements at the shoreline with model simulations, only requiring representative estimates of lake surface temperature.     

Preparing a seismic hazard model for Switzerland: the view from PEGASOS Expert Group 3 (EG1c)

The seismic hazard model used in the PEGASOS project for assessing earth-quake hazard at four NPP sites was a composite of four sub-models, each produced by a team of three experts. In this paper, one of these models is described in detail by the authors. A criticism sometimes levelled at probabilistic seismic hazard studies is that the process by which seismic source zones are arrived at is obscure, subjective and inconsistent. Here, we attempt to recount the stages by which the model evolved, and the decisions made along the way. In particular, a macro-to-micro approach was used, in which three main stages can be described. The first was the characterisation of the overall kinematic model, the “big picture” of regional seismogenesis. Secondly, this was refined to a more detailed seismotectonic model. Lastly, this was used as the basis of individual sources, for which parameters can be assessed. Some basic questions had also to be answered about aspects of the approach to modelling to be used: for instance, is spatial smoothing an appropriate tool to apply? Should individual fault sources be modelled in an intraplate environment? Also, the extent to which alternative modelling decisions should be expressed in a logic tree structure has to be considered.     

Formation and alteration of plagiogranites in an ultramafic-hosted detachment fault at the Mid-Atlantic Ridge (ODP Leg 209)

We examined small-scale shear zones in drillcore samples of abyssal peridotites from the Mid-Atlantic Ridge. These shear zones are associated with veins consisting of chlorite + actinolite/tremolite assemblages, with accessory phases zircon and apatite, and they are interpreted as altered plagiogranite melt impregnations, which originate from hydrous partial melting of gabbroic intrusion in an oceanic detachment fault. Ti-in-zircon thermometry yields temperatures around 820°C for the crystallization of the evolved melt. Reaction path modeling indicates that the alteration assemblage includes serpentine of the adjacent altered peridotites. Based on the model results, we propose that formation of chlorite occurred at higher temperatures than serpentinization, thus leading to strain localization around former plagiogranites during alteration. The detachment fault represents a major pathway for fluids through the oceanic crust, as evidenced by extremely low δ¹⁸O of altered plagiogranite veins (+3.0–4.2‰) and adjacent serpentinites (+ 2.6–3.7‰). The uniform oxygen isotope data indicate that fluid flow in the detachment fault system affected veins and adjacent host serpentinites likewise. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

One‐dimensional dynamics of saturated incompressible porous media: analytical solutions and influence of inertia terms

The complexity of formulations for the hydromechanical coupled mechanics of porous media is typically minimised by simplifying assumptions such as neglecting the effect of inertia terms. For example, three formulations commonly employed to model practical problems are classified as fully dynamic, simplified dynamic and quasi‐static. Thus, depending on the porous media conditions, each formulation will have advantages and limitations. This paper presents a comprehensive analysis of these limitations when solving one‐dimensional fully saturated porous media problems in addition to a new solution that considers a more general loading situation. A phase diagram is developed to assist on the selection of which formulation is more appropriate and convenient regarding particular cases of porosity and hydraulic conductivity values. Non‐dimensional formulations are proposed to achieve this goal. Results using the analytical solutions are compared against numerical values obtained with the finite element method, and the effect of porosity is investigated. Copyright © 2016 John Wiley & Sons, Ltd.     

Holocene climate evolution of the Ugii Nuur basin,Mongolia

In order to evaluate the Holocene palaeoenvironmental evolution of the Ugii Nuur basin, central Mongolia, investigations on chemical and mineralogical properties of lacustrine sediments were carried out on a 630 cm sediment core from lake Ugii Nuur. The interpretation of the record is based on a principal component analysis （PCA） of the elemental composition of the samples. The results show that lacustrine deposition started at 10.6 kyr BP. Low lake level conditions were identified during the Early Holocene （10.6-7.9 kyr BP）. The Mid Holocene （7.9-4.2 kyr BP） was characterized by generally higher lake levels and thus higher moisture supply, but it experienced strong climatic fluctuations. Arid conditions prevailed from 4.2-2.8 kyr BP and were followed by a stable, more humid phase until today.     

Late Glacial and Holocene environmental history of the Pirin Mountains (SW Bulgaria): a paleolimnological study of Lake Dalgoto (2310 m)

Diatoms, Cladocera, and chironomids preserved in the sediments of Lake Dalgoto were studied to reconstruct the history of the lake ecosystem in the context of the vegetation history as represented by the pollen stratigraphy. Younger Dryas silty sediments at the base of the core are characterized by low diversity of aquatic organisms. The transition to the Holocene is indicated by a sharp change from silt to clay-gyttja. The migration and expansion of trees at lower elevations between 10200 and 8500 ¹⁴C-yr BP, along with higher diversities and concentrations of aquatic organisms and the decreased proportion of north-alpine diatoms, point to rapidly rising summer temperatures. After 6500 ¹⁴C-yr BP the expansion of Pinus mugo in the catchment coincides with signs of natural eutrophication as recorded by an increase of planktonic diatoms. In the late Holocene (4000–0 ¹⁴C-yr BP) Pinus peuce and Abies are reduced and Picea expands. Cereal grains and disturbance indicators suggest late-Holocene human modification of the vegetation.