A comparison of helicopter‐borne electromagnetic systems for hydrogeologic studies

The increased application of airborne electromagnetic surveys to hydrogeological studies is driving a demand for data that can consistently be inverted for accurate subsurface resistivity structure from the near surface to depths of several hundred metres. We present an evaluation of three commercial airborne electromagnetic systems over two test blocks in western Nebraska, USA. The selected test blocks are representative of shallow and deep alluvial aquifer systems with low groundwater salinity and an electrically conductive base of aquifer. The aquifer units show significant lithologic heterogeneity and include both modern and ancient river systems. We compared the various data sets to one another and inverted resistivity models to borehole lithology and to ground geophysical models. We find distinct differences among the airborne electromagnetic systems as regards the spatial resolution of models, the depth of investigation, and the ability to recover near‐surface resistivity variations. We further identify systematic biases in some data sets, which we attribute to incomplete or inexact calibration or compensation procedures.     

Modeling unusual eruptive behavior of Mt. Etna,Italy, by means of event bush

One of the best-studied volcanoes of the world, Mt. Etna in Sicily, repeatedly exhibits eruptive scenarios that depart from the behavior commonly considered typical for this volcano. Episodes of intense explosive activity, pyroclastic flows, dome growth and cone collapse pose a variety of previously underestimated threats to human lives in the summit area of the volcano. However, retrospective analysis of these events shows that they were likely caused by the same very sets of premises and starting conditions as “normal” eruptions, yet combined in an unexpected, probably unique, way. To cope with such unexpected consequences, we involve an approach of artificial intelligence developed specially for needs of the geosciences, the event bush. Scenarios inferred from the event bush fit the observed ones and allow to foresee other low-probability events that may occur at the volcano. Application of the event bush provides a more impartial vision of volcanic phenomena and may serve as an intermediary between expert knowledge and numerical assessment, e.g., by means of Bayesian Belief Networks.     

Oxygen isotope compositions of phosphate from arvicoline teeth and Quaternary climatic changes, Gigny, French Jura

Oxygen isotope compositions of biogenic phosphates from mammals are widely used as proxies of the isotopic compositions of meteoric waters that are roughly linearly related to the air temperature at high- and mid-latitudes. An oxygen isotope fractionation equation was determined by using present-day European arvicoline (rodents) tooth phosphate: δ¹⁸O_p = 20.98(±0.59) + 0.572(±0.065) δ¹⁸O_w. This fractionation equation was applied to the Late Pleistocene karstic sequence of Gigny, French Jura. Comparison between the oxygen isotope compositions of arvicoline tooth phosphate and Greenland ice core records suggests to reconsider the previously established hypothetical chronology of the sequence. According to the δ¹⁸O value of meteoric water-mean air temperature relationships, the δ¹⁸O value of arvicoline teeth records variations in mean air temperatures that range from 0° to 15°C.     

Kriging in a global neighborhood

The kriging estimator is usually computed in a moving neighborhood; only the data near the point to be estimated are used. This moving neighborhood approach creates discontinuities in mapping applications. An alternative approach is presented here, whereby all points are estimated using all the available data. To solve the resulting large linear system the kriging estimator is expressed in terms of the inverse of the covariance matrix. The covariance matrix has the advantage of being positive definite and the size of system which can be solved without encountering numerical instability is substantially increased. Because the kriging matrix does not change, the estimator can be written in terms of scalar products, thus avoiding the more time-consuming matrix multiplications of the standard approach. In the particular case of a covariance which is zero for distances greater than a fixed value (the range), the resulting banded structure of the covariance matrix is shown to lead to substantial computational savings in both run time and storage space. In this case the calculation time for the kriging variance is also substantially reduced. The present method is extended to the nonstationary case.     

Analysis of collaboration networks in OpenStreetMap through weighted social multigraph mining

This paper aims to qualify the behaviour of contributors to OpenStreetMap (OSM), a volunteered geographic information (VGI) project, through a multigraph approach. The main purpose is to reproduce contributor’s interactions in a more comprehensive way. First, we define a multigraph that combines existing spatial collaboration networks from the literature with new graphs that illustrate collaboration based on specific aspects of the VGI modes of contribution through semantics, geometry and topology. Indeed, the ways that contributors interact with one another through editing, completion, or even consumption may provide additional information on each user’s operation mode and therefore, on the quality of the contributed data. Social collaborations drawn from indirect criteria – for example, comparisons between contributors’ activity areas – can also be contemplated under another network. Second, the resulting multigraph is analysed using data mining approaches to characterise individuals and identify behavioural groups. The implementation of a multiplex network based on an OSM data sample and an initial analysis make it possible to identify useful behaviours for data qualification. The initial results characterise some contributors as pioneers, moderators and truthful contributors, according to their special roles in the graphs. Mapping elements that include these contributors’ participation are likely to be reliable data     

Composition and origin of holotype Al‐Cu‐Zn minerals in relation to quasicrystals in the Khatyrka meteorite

We investigated the khatyrkite–cupalite holotype sample, 1.2 × 0.5 mm across. It consists of khatyrkite (Cu,Zn)Al₂, cupalite (Cu,Zn)Al, and interstitial material with approximate composition (Zn,Cu)Al₃. All mineral phases of the holotype sample contain Zn and lack Fe that distinguishes them from khatyrkite and cupalite in the Khatyrka meteorite particles (Bindi et al. 2009 , 2011 , 2012 , 2015 ; MacPherson et al. 2013 ; Hollister et al. 2014 ). Neither highly fractionated natural systems nor geo‐ or cosmochemical processes capable of forming the holotype sample are known so far. The bulk chemistry and thermal history of khatyrkite–cupalite assemblage in the holotype sample hint for its possible industrial origin. Likewise, the aluminides in the Khatyrka meteorite particles may also be derived from industrial materials and mixed with extraterrestrial matter during gold prospecting in the Listvenitovy Stream valley.     

A parallel computing thin‐sheet inversion algorithm for airborne time‐domain data utilising a variable overburden

Accurate modelling of the conductivity structure of mineralisations can often be difficult. In order to remedy this, a parametric approach is often used. We have developed a parametric thin‐sheet code, with a variable overburden. The code is capable of performing inversions of time‐domain airborne electromagnetic data, and it has been tested successfully on both synthetic data and field data. The code implements an integral solution containing one or more conductive sheets, buried in a half‐space with a laterally varying conductive overburden. This implementation increases the area of applicability compared to, for example, codes operating in free space, but it comes with a significant increase in computational cost. To minimise the cost, the code is parallelised using OpenMP and heavily optimised, which means that inversions of field data can be performed in hours on multiprocessor desktop computers. The code models the full system transfer function of the electromagnetic system, including variable flight height. The code is demonstrated with a synthetic example imitating a mineralisation buried underneath a conductive meadow. As a field example, the Valen mineral deposit, which is a graphite mineral deposit located in a variable overburden, is successfully inverted. Our results match well with previous models of the deposit; however, our predicted sheet remains inconclusive. These examples collectively demonstrate the effectiveness of our thin‐sheet code.     

Structure and Dynamical Influence of Water Vapor in the Lower Tropical Troposphere

In situ, airborne and satellite measurements are used to characterize the structure of water vapor in the lower tropical troposphere—below the height, \(z_*,\) of the triple-point isotherm, \(T_*.\) The measurements are evaluated in light of understanding of how lower-tropospheric water vapor influences clouds, convection and circulation, through both radiative and thermodynamic effects. Lower-tropospheric water vapor, which concentrates in the first few kilometers above the boundary layer, controls the radiative cooling profile of the boundary layer and lower troposphere. Elevated moist layers originating from a preferred level of convective detrainment induce a profile of radiative cooling that drives circulations which reinforce such features. A theory for this preferred level of cumulus termination is advanced, whereby the difference between \(T_*\) and the temperature at which primary ice forms gives a ‘first-mover advantage’ to glaciating cumulus convection, thereby concentrating the regions of the deepest convection and leading to more clouds and moisture near the triple point. A preferred level of convective detrainment near \(T_*\) implies relative humidity reversals below \(z*\) which are difficult to identify using retrievals from satellite-borne microwave and infrared sounders. Isotopologues retrievals provide a hint of such features and their ability to constrain the structure of the vertical humidity profile merits further study. Nonetheless, it will likely remain challenging to resolve dynamically important aspects of the vertical structure of water vapor from space using only passive sensors.