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.     

Modelling of paleo-saltwater intrusion in the northern part of the Nubian Aquifer System, Northeast Africa

A numerical groundwater model of the Nubian Aquifer System was established to prove the influence of rising seawater levels on the groundwater salinity in northern Egypt over the last 140,000 years. In addition, the impact of a groundwater recharge scenario for these 140,000 years, involving climatic change, on the saltwater/freshwater interface was investigated. Saltwater intrusion induced by rising water levels of the Mediterranean Sea led to salinisation from the Mediterranean Sea to the Qattara depression. This modeling approach was supported by a density-driven model setup and calculation. The modelled saltwater/freshwater interfaces partially fitted the observed ones, especially in the southern half of the Qattara depression. In other parts of the northern Nubian Aquifer System, the ingression of salt water was modelled adequately, but in the west, small regions of the measured interface were not. The development in the Qattara depression (Egypt) and Sirte basin (Libya) were investigated in more detail. The different behaviour in the Sirte basin may be due to high evapotranspiration rates in some former periods, salt solutions from the pre-Quaternary layers or saltwater infiltration from sabkha-like recent salt-bearing sediments.     

Zircon Hf isotope perspective on the origin of granitic rocks from eastern Bavaria, SW Bohemian Massif

The petrogenetic potential of in situ laser ablation Hf isotope data from melt precipitated zircons was explored through the analyses of about 700 individual crystals derived from about 20 different granitic intrusions covering the Variscan basement segment of eastern Bavaria, SE Germany. In combination with geochemical features, four major suites of granitic rocks can be distinguished: (1) NE Bavarian redwitzites (52–57 wt% SiO₂, intrusion ages around 323 Ma) have chondritic εHf(t) values (+0.8 to –0.4). The redwitzites are hybrid rocks and the Hf data are permissive of mixing of a mantle progenitor and crustal melts. (2) Various intermediate rock types (dioritic dyke, granodiorite, palite, 59–63 wt% SiO₂, 334–320 Ma) from the Bavarian Forest yield negative εHf(t) values between –3.4 and –5.1. These values which apparently contradict a mantle contribution fingerprint an enriched (metasomatized) mantle component that was mixed with crustal material. (3) Voluminous, major crust forming granites sensu stricto (67–75 wt% SiO₂, 328–298 Ma) are characterized by a range in εHf(t) values from –0.5 to –5.6. Different crustal sources and/or modification of crustal melts by various input of juvenile material can explain this variation. (4) Post-plutonic (c. 299 Ma) porphyritic dykes of dacitic composition (64–67 wt% SiO₂) from the southern Bavarian Forest have chondritic εHf(t) values (+0.6 to –1.1) and display large intergrain Hf isotope variation. The dykes form a separate petrogenetic group and the Hf data suggest that the zircons crystallized when a pristine mantle-derived parental melt was modified by infiltration of crustal material. The zircon Hf data form a largely coherent positive array with the whole-rock Nd data and both systems yield similar two-stage depleted mantle model ages (1.1–1.7 Ga).     

Measures of Parameter Uncertainty in Geostatistical Estimation and Geostatistical Optimal Design

Studies of site exploration, data assimilation, or geostatistical inversion measure parameter uncertainty in order to assess the optimality of a suggested scheme. This study reviews and discusses measures for parameter uncertainty in spatial estimation. Most measures originate from alphabetic criteria in optimal design and were transferred to geostatistical estimation. Further rather intuitive measures can be found in the geostatistical literature, and some new measures will be suggested in this study. It is shown how these measures relate to the optimality alphabet and to relative entropy. Issues of physical and statistical significance are addressed whenever they arise. Computational feasibility and efficient ways to evaluate the above measures are discussed in this paper, and an illustrative synthetic case study is provided. A major conclusion is that the mean estimation variance and the averaged conditional integral scale are a powerful duo for characterizing conditional parameter uncertainty, with direct correspondence to the well-understood optimality alphabet. This study is based on cokriging generalized to uncertain mean and trends because it is the most general representative of linear spatial estimation within the Bayesian framework. Generalization to kriging and quasi-linear schemes is straightforward. Options for application to non-Gaussian and non-linear problems are discussed.     

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.     

Comparison of seismic activity for Llaima and Villarrica volcanoes prior to and after the Maule 2010 earthquake

Llaima and Villarrica are two of the most active volcanoes in the Chilean Southern Volcanic Zone and presently show contrasting types of activity. Llaima is a closed vent edifice with fumarolic activity, while Villarrica has an open vent with a lava lake, continuous degassing and tremor activity. This study is focused on characterizing the relationships between volcanic and seismic activity in the months before and after the 2010 M8.8 Maule earthquake, which was located in NNW direction from the volcanoes. Time series for tremors, long-period and volcano-tectonic events were obtained from the catalogue of the Volcanic Observatory of the Southern Andes (OVDAS) and from the SFB 574 temporary volcanic network. An increase in the amount of tremor activity, long-period events and degassing rates was observed at Villarrica weeks before the mainshock and continued at a high level also after it. This increase in activity is interpreted to be caused by enhanced magma influx at depth and may be unrelated to the Maule event. In Llaima, an increase in the volcano-tectonic activity was observed directly after the earthquake. The simultaneous post-earthquake activity at both volcanoes is consistent with a structural adjustment response. Since this enhanced activity lasted for more than a year, we suggest that it is related to a medium-term change in the static stress. Thus, the Maule earthquake may have affected both volcanoes, but did not trigger eruptions, from which we assume that none of the volcanoes were in a critical state.     

Influence of plastic fines content on the liquefaction susceptibility of sands: monotonic loading

Acta Geotechnica - The paper presents an experimental study on the effect of plastic fines content on the undrained behavior and liquefaction susceptibility of sand–fines mixtures under...