Big data and machine learning in geoscience and geoengineering: Introduction

In recent years,we have entered the so-called Fourth Paradigm with the regular production of huge amount of observational data.Big data is often characterized by the three‘V's:Volume of data,Variety and Velocity.The concept of big data can potentially address some existing issues in areas of geoscience and geoengineering.Large-scale,comprehensive,multidirectional and multifield geotechnical monitoring is becoming a reality in the very near future.     

Contemporary rates of chemical denudation and atmospheric CO2 sequestration in glacier basins: an Arctic perspective

This paper presents new estimates of solute fluxes from five high Arctic glacier basins in Svalbard. These estimates are combined with data from two other glacier basins to assess the effectiveness of chemical denudation on Svalbard and to estimate rates of temporary (or transient) CO₂ drawdown. We use a solute provenance model to partition solutes into marine, aerosol, atmospheric and crustal components and to estimate their annual fluxes. Crustally derived solute fluxes are equivalent to a mean chemical denudation rate of 350 Σmeq⁺ m⁻² a⁻¹ for Svalbard (range: 160–560 Σmeq⁺ m⁻² a⁻¹), which lies within the global range of 94–4200 Σmeq⁺ m⁻² a⁻¹ for 21 glacier basins in the northern hemisphere, and is close to the continental average of 390 Σmeq⁺ m⁻² a⁻¹. Specific annual discharge is the most significant control upon chemical denudation in the glacierized basins, and basin lithology is an important secondary control, with carbonate‐rich and basaltic lithologies currently showing the greatest chemical denudation rates. Estimates of transient CO₂ drawdown are also directly associated with specific annual discharge and rock type. On Svalbard transient CO₂ drawdown lies in the range 110–3000 kg C km⁻² a⁻¹, whilst the range is 110–13000 kg C km⁻² a⁻¹ for the northern hemisphere glacial data set. Transient CO₂ drawdown is therefore usually low in the Svalbard basins unless carbonate or basalt rocks are abundant. The analysis shows that a large area of uncertainty in the transient CO₂ drawdown estimates exists due to the non‐stoichiometric release of solute during silicate hydrolysis. Silicate hydrolysis is particularly non‐stoichiometric in basins where the extent of glacierization is high, which is most probably an artefact of high flushing rates through ice‐marginal and subglacial environments where K‐feldspars are undergoing mechanical comminution. Copyright © 2000 John Wiley & Sons, Ltd.     

Microseismic Full Waveform Modeling in Anisotropic Media with Moment Tensor Implementation

Seismic anisotropy which is common in shale and fractured rocks will cause travel-time and amplitude discrepancy in different propagation directions. For microseismic monitoring which is often implemented in shale or fractured rocks, seismic anisotropy needs to be carefully accounted for in source location and mechanism determination. We have developed an efficient finite-difference full waveform modeling tool with an arbitrary moment tensor source. The modeling tool is suitable for simulating wave propagation in anisotropic media for microseismic monitoring. As both dislocation and non-double-couple source are often observed in microseismic monitoring, an arbitrary moment tensor source is implemented in our forward modeling tool. The increments of shear stress are equally distributed on the staggered grid to implement an accurate and symmetric moment tensor source. Our modeling tool provides an efficient way to obtain the Green’s function in anisotropic media, which is the key of anisotropic moment tensor inversion and source mechanism characterization in microseismic monitoring. In our research, wavefields in anisotropic media have been carefully simulated and analyzed in both surface array and downhole array. The variation characteristics of travel-time and amplitude of direct P- and S-wave in vertical transverse isotropic media and horizontal transverse isotropic media are distinct, thus providing a feasible way to distinguish and identify the anisotropic type of the subsurface. Analyzing the travel-times and amplitudes of the microseismic data is a feasible way to estimate the orientation and density of the induced cracks in hydraulic fracturing. Our anisotropic modeling tool can be used to generate and analyze microseismic full wavefield with full moment tensor source in anisotropic media, which can help promote the anisotropic interpretation and inversion of field data.     

Risk and trust in the cultural industries

Preliminary claims have been made that working practices within cultural industries such as fashion, music, design and the night time economy may differ from Fordist or modernist arrangements. Cultural firms are often imagined to be more innovative, information-rich, dynamic, flexible, non-hierarchical and dependent on local clusters and networks than their Fordist counterparts (Lash and Urry, 1994). As their impact and significance increase, understanding how creative and cultural industries actually work is of high priority. This paper presents preliminary findings from an on-going ESRC funded study of cultural Micro and Small Enterprises (MSEs) within Manchester, England. Drawing on one element of the project, this paper considers the significance of risk and the importance of social trust for the cultural entrepreneur. Following a discussion of Beck’s development of risk as an analytical concept, and its intersection with Giddens’ notion of ‘active trust’, the paper examines how risk and trust are defined, experienced and negotiated by entrepreneurs in Manchester’s cultural industries. It is suggested that senses of risk are constitutive and often pivotal to the whole economic and social basis of cultural entrepreneurship – risk being central to choices made not only in business but in the lifeworld more generally. The paper then investigates the importance of trust for facilitating as well as countering or offsetting risk. Empirical evidence is presented to show how risk and trust co-relate and interact as constitutive elements within a wider set of shifting relationships between work, leisure and lifestyle in the ‘creative city’.     

High-resolution radio observations of Seyfert galaxies in the extended 12-μm sample – I. The observations

We present 8.4-GHz VLA A-configuration observations of 87 sources from the mid-infrared-selected AGN sample of Rush et al. These 0.25-arcsec-resolution observations allow elongated radio structures tens of pc in size to be resolved, and enable radio components smaller than 3.5 arcsec to be isolated from diffuse galactic disc emission. When combined with previous data, matched radio observations covering 90 per cent of the sample have been collected, and these represent the largest subarcsecond–resolution radio imaging survey of a homogeneously selected sample of Seyfert galaxies to date.
We use our observations to identify five radio-loud AGN in the sample. The nature of the radio emission from Seyfert nuclei will be discussed in subsequent papers.     

Validation of the CME Geomagnetic Forecast Alerts Under the COMESEP Alert System

Under the European Union 7th Framework Programme (EU FP7) project Coronal Mass Ejections and Solar Energetic Particles (COMESEP, http://comesep.aeronomy.be ), an automated space weather alert system has been developed to forecast solar energetic particles (SEP) and coronal mass ejection (CME) risk levels at Earth. The COMESEP alert system uses the automated detection tool called Computer Aided CME Tracking (CACTus) to detect potentially threatening CMEs, a drag-based model (DBM) to predict their arrival, and a CME geoeffectiveness tool (CGFT) to predict their geomagnetic impact. Whenever CACTus detects a halo or partial halo CME and issues an alert, the DBM calculates its arrival time at Earth and the CGFT calculates its geomagnetic risk level. The geomagnetic risk level is calculated based on an estimation of the CME arrival probability and its likely geoeffectiveness, as well as an estimate of the geomagnetic storm duration. We present the evaluation of the CME risk level forecast with the COMESEP alert system based on a study of geoeffective CMEs observed during 2014. The validation of the forecast tool is made by comparing the forecasts with observations. In addition, we test the success rate of the automatic forecasts (without human intervention) against the forecasts with human intervention using advanced versions of the DBM and CGFT (independent tools available at the Hvar Observatory website, http://oh.geof.unizg.hr ). The results indicate that the success rate of the forecast in its current form is unacceptably low for a realistic operation system. Human intervention improves the forecast, but the false-alarm rate remains unacceptably high. We discuss these results and their implications for possible improvement of the COMESEP alert system.     

Germanium isotope measurements of high-temperature geothermal fluids using double-spike hydride generation MC-ICP-MS

We present a double-spike isotope dilution MC-ICP-MS technique for the determination of germanium (Ge) isotope fractionation. Using this technique we determined Ge isotope compositions of geothermal spring fluids, a Columbia River Basalt sample, and an in-house diatom standard. Our technique uses a ⁷³Ge/⁷⁰Ge double spike in combination with hydride generation for Ge extraction from the sample matrix. Fractionation is determined on the ⁷⁴Ge/⁷²Ge mass ratio. The double spike allows us to effectively correct analytical isotope fractionation. Our external standard reproducibility is 0.4‰ (2 SD) over the course of several months. The minimum quantity of Ge needed for isotope analysis is approximately 2 ng. Consistent with previous work on geothermal fluids, Ge in the geothermal spring samples presented here is enriched over Si as compared to low temperature weathering signatures. This observation is typically interpreted as Ge exclusion during silicate mineral precipitation (e.g., quartz). Our isotope results indicate that the analyzed high temperature fluids fractionate Ge isotopes with a range in δ⁷⁴Ge between −0.4‰ and −1.4‰ relative to a Columbia River basalt. We cautiously interpret the observed fractionation as preferential removal of heavy Ge isotopes out of solution during cooling of the hydrothermal fluid and subsequent precipitation of quartz.     

The case for a cognate,polybaric origin for kimberlitic olivines

Kimberlitic olivines typically show a continuous range in size and texture rather than two discrete populations. The cores of small euhedral olivines commonly provide the template for the final crystal shape, which in turn closely matches morphologies produced by crystallization from a moderately under-cooled magma. Cores and edges of the majority of all olivines define a continuous compositional field, which can be interpreted in terms of Raleigh crystallization. Marked chemical gradients at the olivine margins are linked to rapid physico-chemical changes to the magma associated with loss of volatiles during the late stages of emplacement. Thus, rapid crystallization of groundmass olivines would deplete the magma in Ni, but increase Ca activity. The latter would be enhanced by decreasing pressure coupled with loss of CO₂ from the carbonate-bearing kimberlite magma.For mantle olivines and the most refractory olivines in kimberlites (~ Fo₉₄) to be in equilibrium with bulk rock compositions matching those of Mg-rich macrocrystic and aphanitic kimberlites (Mg# ~ 88) requires a mineral-melt Mg–Fe distribution coefficient of 0.47. This is well within the experimentally determined range for this distribution coefficient in carbonate-bearing systems. In southern African post-Gondwana alkaline pipe clusters, the average bulk rock Mg# and composition of the associated most Mg-rich olivine both decrease sympathetically from the interior to the continental margin, which is also consistent with a cognate origin for the olivines.A kimberlite magma following a plausible P-T trajectory relative to the CO₂/H₂O peridotite solidus would initially experience superheating, resulting in partial resorption of early-formed olivines that crystallized on the cool conduit walls. It would become supersaturated as it crossed the carbonated peridotite “ledge”, resulting in tabular and hopper growth forms typical of euhedral olivine cores. With further ascent, the magma would once again become superheated, resulting in partial resorption of these cores. Thus, apparently complex textures and internal zonation patterns of kimberlitic olivines are predicted by a plausible magma P-T trajectory.     

Optimal Field Approaches for Electrokinetic In Situ Oxidation Remediation

Numerical simulations were used to identify and evaluate optimum electrode configurations and approaches for electrokinetic in situ chemical oxidation (EK‐ISCO) remediation of low‐permeability sediments. A newly developed groundwater and EK flow and reactive transport numerical model was used to conduct two‐dimensional scenario simulations of the coverage of an injected oxidant, permanganate, and the oxidation of a typical organic contaminant (tetrachloroethene, PCE). For linear configurations of vertical electrodes, the spacing of same‐polarity electrodes is recommended to be about one‐third to one‐quarter of the anode–cathode spacing. Greater coverage could also be achieved by locating additional oxidant injection wells at the divergence of the electric field in linear electrode configurations. Horizontal electrodes allowed greater contact between the injected permanganate and PCE and resulted in faster degradation of PCE compared to vertical electrodes. Pulsed oxidant injection, closer electrode spacing, and electric field reversal also resulted in faster EK‐ISCO remediation.