The importance of defining chemical potentials,substitution mechanisms and solubility in trace element diffusion studies: the case of Zr and Hf in olivine

The diffusion, substitution mechanism and solubility limits of Zr and Hf in synthetic forsterite (Mg₂SiO₄) and San Carlos olivine (Mg_0.9Fe_0.1)₂SiO₄ have been investigated between 1,200 and 1,500 °C as a function of the chemical potentials of the components in the system MgO(FeO)–SiO₂–ZrO₂(HfO₂). The effect of oxygen fugacity and crystallographic orientation were also investigated. The solubilities of Zr in forsterite are highest and diffusion fastest when the coexisting three-phase source assemblage includes ZrSiO₄ (zircon) or HfSiO₄ (hafnon), and lower and slower, respectively, when the source assemblage includes MgO (periclase). This indicates that Zr and Hf substitute on the octahedral sites in olivine, charge balanced by magnesium vacancies. Diffusion is anisotropic, with rates along the crystal axes increasing in the order a < b < c. The generalized diffusion relationship as a function of chemical activity (as \(a_{{{\text{SiO}}_{2} }}\)), orientation and temperature is: \(logD_{\text{Zr}} = \frac{1}{4}loga_{{{\text{SiO}}_{2} }} + logD_{0} - \left( {\frac{{368 \pm 17\;{\text{kJ}}\;{\text{mol}}^{ - 1} }}{{2.303\;{\text{RT}}}}} \right)\) where the values of log D ₀ are ?3.8(±0.5), ?3.4(±0.5) and ?3.1(±0.5) along the a, b and c axes, respectively. Most experiments were conducted in air (fO₂ = 10^?0.68 bars), but one at fO₂ = 10^?11.2 bars at 1,400 °C shows no resolvable effect of oxygen fugacity on Zr diffusion. Hf is slightly more soluble in olivine than Zr, but diffuses slightly slower. Diffusivities of Zr in experiments in San Carlos olivine at 1,400 °C, fO₂ = 10^?6.6 bars are similar to those in forsterite at the same conditions, showing that the controls on diffusivities are adequately captured by the simple system (nominally iron-free) experiments. Diffusivities are in good agreement with those measured by Spandler and O’Neill (Contrib Miner Petrol 159:791–818, 2010) in San Carlos olivine using silicate melt as the source at 1,300 °C, and fall within the range of most measurements of Fe–Mg inter-diffusion in olivine at this temperature. Forsterite–melt partitioning experiments in the CaO–MgO–Al₂O₃–SiO₂–ZrO₂/HfO₂ show that the interface concentrations from the diffusion experiments represent true equilibrium solubilities. Another test of internal consistency is that the ratios of the interface concentrations between experiments buffered by Mg₂SiO₄ + Mg₂Si₂O₆ + ZrSiO₄ or Mg₂SiO₄ + ZrSiO₄ + ZrO₂ (high silica activity) to those buffered by Mg₂SiO₄ + MgO + ZrO₂ (low silica activity) agree well with the ratios calculated from thermodynamic data. This study highlights the importance of buffering chemical potentials in diffusion experiments to provide constraints on the interface diffusant concentrations and hence validate the assumption of interface equilibrium.     

Thermodynamic properties of chlorite and berthierine derived from calorimetric measurements

In the context of the deep waste disposal, we have investigated the respective stabilities of two iron-bearing clay minerals: berthierine ISGS from Illinois [USA; (Al_0.975FeIII_0.182FeII_1.422Mg_0.157Li_0.035Mn_0.002)(Si_1.332Al_0.668)O₅(OH)₄] and chlorite CCa-2 from Flagstaff Hill, California [USA; (Si_2.633Al_1.367)(Al_1.116FeIII_0.215Mg_2.952FeII_1.712Mn_0.012Ca_0.011)O₁₀(OH)₈]. For berthierine, the complete thermodynamic dataset was determined at 1 bar and from 2 to 310 K, using calorimetric methods. The standard enthalpies of formation were obtained by solution-reaction calorimetry at 298.15 K, and the heat capacities were measured by heat-pulse calorimetry. For chlorite, the standard enthalpy of formation is measured by solution-reaction calorimetry at 298.15 K. This is completing the entropy and heat capacity obtained previously by Gailhanou et al. (Geochim Cosmochim Acta 73:4738–4749, 2009) between 2 and 520 K, by using low-temperature adiabatic calorimetry and differential scanning calorimetry. For both minerals, the standard entropies and the Gibbs free energies of formation at 298.15 K were then calculated. An assessment of the measured properties could be carried out with respect to literature data. Eventually, the thermodynamic dataset allowed realizing theoretical calculations concerning the berthierine to chlorite transition. The latter showed that, from a thermodynamic viewpoint, the main factor controlling this transition is probably the composition of the berthierine and chlorite minerals and the nature of the secondary minerals rather than temperature.     

Tsunami hazard assessment of Canada

We present a preliminary probabilistic tsunami hazard assessment of Canadian coastlines from local and far-field, earthquake, and large submarine landslide sources. Analyses involve published historical, palaeotsunami and palaeoseismic data, modelling, and empirical relations between fault area, earthquake magnitude, and tsunami run-up. The cumulative estimated tsunami hazard for potentially damaging run-up (≥1.5 m) of the outer Pacific coastline is ~40–80 % in 50 years, respectively one and two orders of magnitude greater than the outer Atlantic (~1–15 %) and the Arctic (<1 %). For larger run-up with significant damage potential (≥3 m), Pacific hazard is ~10–30 % in 50 years, again much larger than both the Atlantic (~1–5 %) and Arctic (<1 %). For outer Pacific coastlines, the ≥1.5 m run-up hazard is dominated by far-field subduction zones, but the probability of run-up ≥3 m is highest for local megathrust sources, particularly the Cascadia subduction zone; thrust sources further north are also significant, as illustrated by the 2012 Haida Gwaii event. For Juan de Fuca and Georgia Straits, the Cascadia megathrust dominates the hazard at both levels. Tsunami hazard on the Atlantic coastline is dominated by poorly constrained far-field subduction sources; a lesser hazard is posed by near-field continental slope failures similar to the 1929 Grand Banks event. Tsunami hazard on the Arctic coastline is poorly constrained, but is likely dominated by continental slope failures; a hypothetical earthquake source beneath the Mackenzie delta requires further study. We highlight areas susceptible to locally damaging landslide-generated tsunamis, but do not quantify the hazard.     

High-resolution aquifer analog of fluvial–aeolian sediments of the Guarani aquifer system

The Guarani aquifer system (GAS) represents one of the biggest aquifers in the world and is the most relevant groundwater resource in South America. For the first time, by combining field and laboratory measurements, a high-resolution aquifer analog model of fluvial–aeolian sediments of the GAS in São Paulo State (Brazil) is constructed. Three parallel sections of frontal outcrops, 28 m × 5.8 m, and two parallel sections of lateral outcrops, 7 m × 5.8 m, are recorded during open-pit mining of sandy sediments and describe in detail the three-dimensional distribution of the local lithofacies and hydrofacies. Variations of hydraulic conductivity, K, and porosity, n, are resolved on the centimeter scale, and the most permeable units of the fluvial–aeolian facies association are identified. The constructed aquifer analog model shows moderate hydraulic heterogeneity and a mean K value of 1.36 × 10^?4 m/s, which is greater than the reported range of K values for the entire GAS in São Paulo State. The results suggest that the examined sedimentary unit constitutes a relevant portion of the GAS in São Paulo State in the context of groundwater extraction and pollution. Moreover, the constructed aquifer analog is considered an ideal basis for future numerical model experiments, aiming at in-depth understanding of the groundwater flow and contaminant transport patterns at this GAS portion or at comparable fluvial–aeolian facies associations.     

Web information monitoring and crowdsourcing for promoting and enhancing the Algerian geoheritage

The documentation and management of information concerning geological sites are important activities in all the steps of the geoheritage conservation and its day-to-day management. This documentation is mostly available on the web; it includes varied types of data such as pictures, scientific publications, news, and blog articles. In order to ensure the efficient management of the geoheritage resources, scientists and people involved in the promotion of the geoheritage and geotourism should have access to this documentation and kept updated of the most recent and relevant information. In this study, we propose the development of a web information monitoring system that automates the search and the collection of different types of documents related to the Algerian geosites available on the web. It is designed as a tool that continuously browses the web using the most popular search engines and then disseminates the results to the users through periodic newsletters. This proposed system will also serve to chronologically organize all the pictures of geosites gathered from the web and thus in order to study the evolution of Algerian geosites over time. For this purpose, we use one of the crowdsourcing methods in the field of knowledge management that is games with a purpose (GWAP). The idea here is to involve people in sorting pictures of geosites through a playful interface. The results obtained show that this system is a valuable tool that can assist geoscientists and decision makers in enhancing the geoheritage.     

Selecting time windows of seismic phases and noise for engineering seismology applications: a versatile methodology and algorithm

Seismic signal windowing is the preliminary step for many analysis procedures in engineering seismology (standard spectral ratio, quality factor, general inversion techniques, etc.). Moreover a noise window is often necessary for the data quality control through the signal-to-noise verification. Selecting the noise window can be challenging when large heterogeneous datasets are considered, especially when they include short pre-event noise signals. This study proposes a fully automatic and configurable (i.e., with default parameters that can also be user-defined) algorithm to windowing the noise and the P, S, coda and full signal once the P-wave (T _P ) and S-wave (T _S ) first arrivals are known. An application example is given on a KiK-net dataset. A Matlab language implementation of this algorithm is proposed as an online resource.     

Wave effects in global ocean modeling: parametrizations vs. forcing from a wave model

One of the main challenges of the Copernicus Marine Service is the implementation of coupled ocean/waves systems that accurately estimate the momentum and energy fluxes provided by the atmosphere to the ocean. This study aims to investigate the impact of forcing the Nucleus for European Modelling of the Ocean (NEMO) ocean model with forecasts from the wave model of Météo-France (MFWAM) to improve classical air-sea flux parametrizations, these latter being mostly driven by the 10-m wind. Three wave-related processes, namely, wave-state-dependent stress, Stokes drift-related effects (Stokes-Coriolis force, Stokes drift advection on tracers and on mass), and wave-state-dependent surface turbulence, are examined at a global scale with a horizontal resolution of 0.25°. Three years of sensitivity simulations (2014–2016) show positive feedback on sea surface temperature (SST) and currents when the wave model is used. A significant reduction in SST bias is observed in the tropical Atlantic Ocean. This is mainly due to the more realistic momentum flux provided by the wave model. In mid-latitudes, the most interesting impact occurs during the summer stratification, when the wind is low and the wave model produces a reduction in the turbulence linked with wave breaking. Magnitudes of the large-scale currents in the equatorial region are also improved by 10% compared to observations. In general, it is shown that using the wave model reduces on average the momentum and energy fluxes to the ocean in tropical regions, but increases them in mid-latitudes. These differences are in the order of 10 to 20% compared with the classical parametrizations found in stand-alone ocean models.     

Glacial lake depth and volume estimation based on a large bathymetric dataset from the Cordillera Blanca,Peru

Glacial lakes are most often located in remote places making it difficult to carry out detailed bathymetric surveys. Consequently, lake depths and volumes for unmeasured lakes are often estimated using empirical relationships developed mainly from small bathymetric datasets. In this study, we use the bathymetry dataset of the Cordillera Blanca, Peru comprising 121 detailed lake bathymetries, the most extensive dataset in the world. We assess the performance of the most commonly applied empirical relationships for lake mean depth and volume estimation, but also investigate relationships between different geometric lake variables. We find that lake volume estimation performs better when derived from lake mean depth, which in turn is estimated from lake width. The findings also reveal the extreme variability of lake geometry, which depends on glacio-geomorphological processes that empirical–statistical relationships cannot adequately represent. Such relationships involve characteristic uncertainty ranges of roughly ±50%. We also estimate potential peak discharges of outburst floods from these lakes by applying empirical relationships from the literature, which results in discharges varying by up to one-order of magnitude. Finally, the results are applied to the 860 lakes without bathymetric measurements from the inventory dataset of the Cordillera Blanca to estimate lake mean depth, volume and possible peak discharge for all unmeasured lakes. Estimations show that ca. 70% (610) of the lakes have a mean depth lower than 10 m and very few longer than 40 m. Lake volume of unmeasured lakes represent ca. 32% (5.18 × 10⁸ m³) of the total lake volume (1.15 × 10⁹ m³) in the Cordillera Blanca. Approximately, 50% of the lakes have potential peak discharges > 1000 m³/s in case of lake outburst floods, implying a need for additional studies for risk assessment. © 2020 John Wiley & Sons, Ltd.     

Properties of dust source material and volcanic ash in Iceland

The volcanic origin, primarily basaltic, of most of the surface material in Iceland influences its physical properties and appearance. Size distributions, shape analyses and melting experiments were made for surface material collected in high-erosion dust source areas and fresh volcanic ash deposits to determine whether they differ from one another and from dust from other major dust sources. The major differences found between Icelandic dust and dust from other major dust sources in the world, such as the Sahara, are in the particle shapes, lower density and darker colour. Icelandic dust particles greater than 20 μm retain volcanic morphological properties that are also found in fresh volcanic ash. Dust and fresh volcanic ash particles less than 20 μm are crystalline and blocky in nature, similar to the dust from other global source regions. The finer grained (<20 μm) Icelandic particles will have similar suspension and transport behaviours and be similarly hazardous to health and infrastructure as non-Icelandic dust. The coarser particles (>20 μm) will have different suspension and transport behaviours than other dusts due to the volcanic morphology. Icelandic surface material has between 5% and 30% glassy particles compared to fresh volcanic ash which has more than 50% glassy particles. Glassy particles were observed to melt at a lower temperature than the mineral grains; and, as a result, volcanic ash is found to be more threatening to aircraft engines than the typical dust from Iceland. Icelandic dust was observed to be blocky, or plate-like in the respirable size fraction, suggesting similar health hazards as dust from other regions.