A Top-Cut Model for Deposits with Heavy-Tailed Grade Distribution

In some ore deposits, the grade distribution is heavy-tailed and high values make the inference of first- and second-order statistics nonrobust. In the case of gold data, high values are usually cut and the block estimation is performed using truncated grades. With this method, a part of the deposit is omitted, resulting in a potential bias on resources figures. Ad-hoc corrections can be applied on the final figures to take into account the tail of the grade distribution, but no theoretical guidelines are available. A geostatistical model has been developed to handle high values based on the assumption that for high grade zones, the only tangible information is the geometry. The grade variable z can be split into three parts: the truncated grade ( $\operatorname{Min} (z, z_{\mathrm{e}})$ where z _e is the top-cut grade), a weighted indicator above top-cut grade (1{z≥z _e}), and a residual. Within this framework, the residual is poorly structured, and in most cases is a pure nugget-effect. Moreover, it has no spatial correlation with the truncated grade and the indicator above top-cut grade. This decomposition makes the variographic study more robust because variables (indicator and truncated grade) do not keep high grade values. The underlying hypotheses can be tested on experimental indicator variograms and the top-cut grade can be justified. Finally, the estimation is based on a truncated grade and indicator cokriging. The model is applied to blast holes from a gold deposit and on a simulated example. Both cases illustrate the benefits of keeping the high values in the estimation process via an appropriate mathematical model.     

New geological model of the Lagoa Real uraniferous albitites from Bahia (Brazil)

New evidence supported by petrography (including mineral chemistry), lithogeochemistry, U-Pb geochronology by Laser Ablation Inductively Coupled Plasma Mass Spectrometry (LA-ICP-MS), and physicochemical study of fluid and melt inclusions by LA-ICP-MS and microthermometry, point to an orogenic setting of Lagoa Real (Bahia-Brazil) involving uraniferous mineralization. Unlike the previous models in which uraniferous albitites represent Na-metasomatised 1.75 Ga anorogenic granitic rocks, it is understood here that they correspond to metamorphosed sodium-rich and quartz-free 1.9 Ga late-orogenic syenitic rocks (Na-metasyenites). These syenitic rocks are rich not only in albite, but also in U-rich titanite (source of uranium). The interpretation of geochemical data points to a petrogenetic connection between alkali-diorite (local amphibolite protolith) and sodic syenite by fractional crystallization through a transalkaline series. This magmatic differentiation occurred either before or during shear processes, which in turn led to albitite and amphibolite formation. The metamorphic reactions, which include intense recrystallization of magmatic minerals, led uraninite to precipitate at 1.87 Ga under Oxidation/Reduction control. A second population of uraninites was also generated by the reactivation of shear zones during the 0.6 Ga Brasiliano Orogeny. The geotectonic implications include the importance of the Orosirian event in the Paramirim Block during paleoproterozoic S?o Francisco Craton edification and the influence of the Brasiliano event in the Paramirim Block during the West-Gondwana assembly processes. The regional microcline-gneiss, whose protolith is a 2.0 Ga syn-collisional potassic granite, represents the albitite host rock. The microcilne-gneiss has no petrogenetic association to the syenite (albitite protolith) in magmatic evolutionary terms.     

Observations of Coronal Mass Ejections with the Coronal Multichannel Polarimeter

The Coronal Multichannel Polarimeter (CoMP) measures not only the polarization of coronal emission, but also the full radiance profiles of coronal emission lines. For the first time, CoMP observations provide high-cadence image sequences of the coronal line intensity, Doppler shift, and line width simultaneously over a large field of view. By studying the Doppler shift and line width we may explore more of the physical processes of the initiation and propagation of coronal mass ejections (CMEs). Here we identify a list of CMEs observed by CoMP and present the first results of these observations. Our preliminary analysis shows that CMEs are usually associated with greatly increased Doppler shift and enhanced line width. These new observations provide not only valuable information to constrain CME models and probe various processes during the initial propagation of CMEs in the low corona, but also offer a possible cost-effective and low-risk means of space-weather monitoring.     

Reionization and the Cosmic Dawn with the Square Kilometre Array

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21 cm line from the earliest phases of star and galaxy formation in the Universe. This 21 cm signal provides a new and unique window both on the time of the formation of the first stars and accreting black holes and the subsequent period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.     

A first deep seismic survey in the Sea of Marmara: Deep basins and whole crust architecture and evolution

Increased source strength, streamer length and dense spatial coverage of seismic reflection profiles of the SEISMARMARA Leg 1 allow to image the deep structure of the marine North Marmara Trough (NMT) on the strike-slip North Anatolian Fault (NAF) west of the destructive Izmit 1999 earthquake. A reflective lower crust and the Moho boundary are detected. They appear upwarped on an E-W profile from the southern Central Basin eastwards, towards more internal parts of the deformed region. Thinning of the upper crust could use a detachment suggested from an imaged dipping intracrustal reflector that would allow upper crustal material to be dragged from beneath it and above the lower crust, accounting for the extensional component but also southwest motion of the southern margin of the NMT. Sections across the eastern half of the NMT, crossing the Cinarcik and Imrali basins, reveal several faults that are active reaching into the basement and have varying strike and proportions of normal and strike-slip displacement. They might be viewed as petals of a large scale negative flower-structure that spreads over a width of 30 km at surface and is rooted deeper in the lithosphere. Under the Central Basin a very thick sediment infill is revealed and its extensional bounding faults are active and imaged as much as 8 km apart down to 6 km depth. We interpret them as two deep-rooted faults encompassing a foundering basement block, rather than being merely pulled-apart from a jog in a strike-slip above a décollement. The deep-basin lengthening would account for only a modest part of the proposed 60 km finite motion since 4 Myr along the same direction oblique to the NMT that sidesteps the shear motion from its two ends. Thus differential motion occurred much beyond the deep basins, like subsidence involving the NMT bounding faults and the intracrustal detachments. The complex partitioned motion localized on active faults with diverse natures and orientations is suggested to represent the overburden deformation induced from horizontal plane simple shear occurring in depth at lithospheric scale, and in front of the North Anatolian Fault when it propagated through the region.     

Pb,Hf and Nd isotope compositions of the two Réunion volcanoes (Indian Ocean): A tale of two small-scale mantle “blobs”?

Pb, Hf, Nd and Sr isotopes of basaltic lavas from the two Réunion Island volcanoes are reported in order to examine the origin of the sources feeding these volcanoes and to detect possible changes through time. Samples, chosen to cover the whole lifetime of the two volcanoes (from 2 Ma to present), yield a chemically restricted (compared to OIB lavas) but complex distribution. Réunion plume isotopic characteristics have been defined on the basis of the composition of uncontaminated shield-building lavas from the Piton de la Fournaise volcano. The average ?_Nd, ?_Hf, ⁸⁷Sr/⁸⁶Sr and ²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb, and ²⁰⁸Pb/²⁰⁴Pb isotope ratios calculated for this component are + 4.4, + 9.1, 0.70411, 18.97, 15.59 and 39.03, respectively. In Pb–Pb isotope space, each volcano defines a distinct linear trend but slight variations are also detected within the various volcanic sequences. The Piton des Neiges volcano yields a distinct and significantly more scattered isotopic distribution than Piton de la Fournaise for both Pb, Hf and Nd isotope tracers. A principal component analysis of the Pb isotope data from Piton de la Fournaise reveals a major contribution of the C and EM-1 components (with a clear Dupal flavor) as main components for the modern Réunion plume. The same components have been identified for Piton des Neiges but with a stronger participation of a depleted mantle component and a weaker EM-1 contribution. The compositional change of the lavas erupted by the Piton des Neiges and Piton de la Fournaise volcanoes is attributed to the impingement of two small-scale blobs of plume material at the base of the Réunion lithosphere. Compared to other hot-spots worldwide, in particular Hawaii and Kerguelen, magmas beneath Réunion are generated from a considerably more homogeneous, compositionally more primitive plume higher in ²⁰⁶Pb. Although shallow-level contamination processes have been locally detected they did not alter significantly the composition of the plume magmas. This is tentatively attributed to mantle dynamics producing small, high-velocity blobs that ascend rapidly through the lithosphere, and to the lack of a well-developed magma chamber at depth in the lithosphere.     

Inversion of time-dependent nuclear well-logging data using neural networks

The purpose of this work was to investigate a new and fast inversion methodology for the prediction of subsurface formation properties such as porosity, salinity and oil saturation, using time‐dependent nuclear well logging data. Although the ultimate aim is to apply the technique to real‐field data, an initial investigation as described in this paper, was first required; this has been carried out using simulation results from the time‐dependent radiation transport problem within a borehole. Simulated neutron and γ‐ray fluxes at two sodium iodide (NaI) detectors, one near and one far from a pulsed neutron source emitting at ～14 MeV, were used for the investigation. A total of 67 energy groups from the BUGLE96 cross section library together with 567 property combinations were employed for the original flux response generation, achieved by solving numerically the time‐dependent Boltzmann radiation transport equation in its even parity form. Material property combinations (scenarios) and their correspondent teaching outputs (flux response at detectors) are used to train the Artificial Neural Networks (ANNs) and test data is used to assess the accuracy of the ANNs. The trained networks are then used to produce a surrogate model of the expensive, in terms of computational time and resources, forward model with which a simple inversion method is applied to calculate material properties from the time evolution of flux responses at the two detectors. The inversion technique uses a fast surrogate model comprising 8026 artificial neural networks, which consist of an input layer with three input units (neurons) for porosity, salinity and oil saturation; and two hidden layers and one output neuron representing the scalar photon or neutron flux prediction at the detector. This is the first time this technique has been applied to invert pulsed neutron logging tool information and the results produced are very promising. The next step in the procedure is to apply the methodology to real data.     

Using a groundwater quality negotiation support system to change land-use management near a drinking-water abstraction in the Netherlands

A negotiation support system (NSS) was developed to solve groundwater conflicts that arose during land-use management. It was set up in cooperation with the stakeholders involved to provide information on the impact of land use, e.g., agriculture, nature (forested areas), recreation, and urban areas, on the quality of both infiltrating and abstracted groundwater. This NSS combined simulation programs that calculate (1) the concentrations of nitrate in shallow groundwater for each land-use area and (2) the transport of nitrate in the groundwater-saturated zone. The user interface of the NSS enabled scenario analyses. The NSS was validated at a drinking-water abstraction near Holten (the Netherlands) using a spatial planning process aimed at sustainable land-use and groundwater-resource management. Two land-use scenarios were considered: a base scenario reflecting the autonomous development and an adapted land-use scenario. The calculated results for shallow groundwater provided an explicit spatial overview of the impact of historical land use and N application on the quality of abstracted groundwater as well as insight into the impact of changes in land use and N application. Visualization of the conflicting interests of agriculture and the drinking-water abstraction helped all stakeholders accept the necessary changes in land use identified by the adapted land-use scenario of the NSS. These changes were included in the preferred land-use management option in the regional planning process, which has since been formalized. The NSS provided system insight, scoping analyses, and education, in addition to generating quantitative information on the impact of land-use functions on groundwater quality.