A method of alternating characteristics with application to advection-dominated environmental systems

We present a numerical method for solving a class of systems of partial differential equations (PDEs) that arises in modeling environmental processes undergoing advection and biogeochemical reactions. The salient feature of these PDEs is that all partial derivatives appear in linear expressions. As a result, the system can be viewed as a set of ordinary differential equations (ODEs), albeit each one along a different characteristic. The method then consists of alternating between equations and integrating each one step-wise along its own characteristic, thus creating a customized grid on which solutions are computed. Since the solutions of such PDEs are generally smoother along their characteristics, the method offers the potential of using larger time steps while maintaining accuracy and reducing numerical dispersion. The advantages in efficiency and accuracy of the proposed method are demonstrated in two illustrative examples that simulate depth-resolved reactive transport and soil carbon cycling.     

Defining the Potential of Nanoscale Re‐Os Isotope Systematics Using Atom Probe Microscopy

Atom probe microscopy (APM) is a relatively new in situ tool for measuring isotope fractions from nanoscale volumes (< 0.01 μm³). We calculate the theoretical detectable difference of an isotope ratio measurement result from APM using counting statistics of a hypothetical data set to be ± 4δ or 0.4% (2s). However, challenges associated with APM measurements (e.g., peak ranging, hydride formation and isobaric interferences), result in larger uncertainties if not properly accounted for. We evaluate these factors for Re‐Os isotope ratio measurements by comparing APM and negative thermal ionisation mass spectrometry (N‐TIMS) measurement results of pure Os, pure Re, and two synthetic Re‐Os‐bearing alloys from Schwander et al. (2015, Meteoritics and Planetary Science, 50, 893) [the original metal alloy (HSE) and alloys produced by heating HSE within silicate liquid (SYN)]. From this, we propose a current best practice for APM Re‐Os isotope ratio measurements. Using this refined approach, mean APM and N‐TIMS ¹⁸⁷Os/¹⁸⁹Os measurement results agree within 0.05% and 2s (pure Os), 0.6–2% and 2s (SYN) and 5–10% (HSE). The good agreement of N‐TIMS and APM ¹⁸⁷Os/¹⁸⁹Os measurements confirms that APM can extract robust isotope ratios. Therefore, this approach permits nanoscale isotope measurements of Os‐bearing alloys using the Re‐Os geochronometer that could not be measured by conventional measurement principles.     

Use of photo-based 3D photogrammetry in analysing the results of laboratory pressure grouting tests

This paper presents a non-destructive, low-cost, photo-based, 3D reconstruction technique for characterizing geo-materials with irregular shapes of a relatively large size. After being validated against two traditional volume measurement methods, namely the vernier caliper method and the fluid displacement method for regular and irregular shapes, respectively, 3D photogrammetry was used to analyse the grout bulbs formed in laboratory pressure grouting tests. The reconstructed 3D mesh model of the sample provides accurate and detailed 3D vertex data, which allowed the volume, densification efficiency and bleeding behaviour of the grout bulbs to be analysed. Comparing the bulb section views at different grouting pressures also offers an intuitive observation of the grout development and propagation process. Moreover, the 3D vertex data and surface area included in the model are of great importance in validating numerical predictions of the pressure grouting process and analysing the interface shear resistance of grouted soil nails or anchors. Compared to existing approaches, the new 3D photogrammetry method possesses several key advantages: (a) it does not require expensive, specialized equipment; (b) samples are not destroyed or modified during testing; (c) it allows to reconstruct objects of various scales and (d) the software is public domain. Therefore, the adoption of this 3D photogrammetry method will facilitate research in the pressure grouting process and can be extended to other problems in geotechnical engineering.     

Magma evolution beneath Bequia,Lesser Antilles,deduced from petrology of lavas and plutonic xenoliths

Extrusive and intrusive igneous rocks represent different parts of a magmatic system and ultimately provide complementary information about the processes operating beneath volcanoes. To shed light on such processes, we have examined and quantified the textures and mineral compositions of plutonic and cumulate xenoliths and lavas from Bequia, Lesser Antilles arc. Both suites contain assemblages of iddingsitized olivine, plagioclase, clinopyroxene and spinel with rare orthopyroxene and ilmenite. Mineral zoning is widespread, but more protracted in lavas than xenoliths. Plagioclase cores and olivine have high anorthite (An?≤?98) and low forsterite (Fo?≤?84) compositions respectively, implying crystallisation from a hydrous mafic melt that was already fractionated. Xenolith textures range from adcumulate to orthocumulate with variable mineral crystallisation sequences. Textural criteria are used to organize the xenoliths into six groups. Amphibole, notably absent from lavas, is a common feature of xenoliths, together with minor biotite and apatite. Bulk compositions of xenoliths deviate from the liquid line of descent of lavas supporting a cumulate origin with varying degrees of reactive infiltration by evolved hydrous melts, preserved as melt inclusions in xenolith crystals. Volatile saturation pressures in melt inclusions indicate cumulate crystallization over a 162–571 MPa pressure range under conditions of high dissolved water contents (up to 7.8 wt% H₂O), consistent with a variety of other thermobarometric estimates. Phase assemblages of xenoliths are consistent with published experimental data on volatile-saturated low-magnesium and high-alumina basalts and basaltic andesite from the Lesser Antilles at pressures of 200–1000 MPa, temperatures of 950–1050 °C and dissolved H₂O contents of 4–7 wt%. Once extracted from mid-crustal mushes, residual melts ascend to higher levels and undergo H₂O-saturated crystallization in shallow, pre-eruptive reservoirs to form phenocrysts and glomerocrysts. The absence of amphibole from lavas reflects instability at low pressures, whereas its abundance in xenoliths testifies to its importance in mid-crustal differentiation processes. A complex, vertically extensive (6 to at least 21 km depth) magmatic system is inferred beneath Bequia. Xenoliths represent fragments of the mush incorporated into ascending magmas. The widespread occurrence of evolved melts in the mush, but the absence of erupted evolved magmas, in contrast to islands in the northern Lesser Antilles, may reflect the relative immaturity of the Bequia magmatic system.     

Characterisation of primary and secondary carbonates in hypabyssal kimberlites: an integrated compositional and Sr-isotopic approach

Carbonates in fresh hypabyssal kimberlites worldwide have been studied to understand their origin [i.e. primary magmatic (high T) versus deuteric (‘low T’) versus hydrothermal/alteration (‘low T’)] and identify optimal strategies for petrogenetic studies of kimberlitic carbonates. The approach presented here integrates detailed textural characterisation, cathodoluminescence (CL) imaging, in situ major- and trace-element analysis, as well as in situ Sr-isotope analysis. The results reveal a wide textural diversity. Calcite occurs as fine-grained groundmass, larger laths, segregations, veins or as a late crystallising phase, replacing olivine or early carbonates. Different generations of carbonates commonly coexist in the same kimberlite, each one defined by a characteristic texture, CL response and composition (e.g., variable Sr and Ba concentrations). In situ Sr isotope analysis revealed a magmatic signature for most of the carbonates, based on comparable ⁸⁷Sr/⁸⁶Sr values between these carbonates and the coexisting perovskite, a robust magmatic phase. However, this study also shows that in situ Sr isotope analysis not always allow distinction between primary (i.e., magmatic) and texturally secondary carbonates within the same sample. Carbonates with a clear secondary origin (e.g., late-stage veins) occasionally show the same moderately depleted ⁸⁷Sr/⁸⁶Sr ratios of primary carbonates and coexisting perovskite (e.g., calcite laths-shaped crystals with ⁸⁷Sr/⁸⁶Sr values identical within uncertainty to those of vein calcite in the De Beers kimberlite). This complexity emphasises the necessity of integrating detailed petrography, geochemical and in situ Sr isotopic analyses for an accurate interpretation of carbonate petrogenesis in kimberlites. Therefore, the complex petrogenesis of carbonates demonstrated here not only highlights the compositional variability of kimberlites, but also raises concerns about the use of bulk-carbonate C-O isotope studies to characterise the parental melt compositions. Conversely, our integrated textural and in situ study successfully identifies the most appropriate (i.e. primary) carbonates for providing constraints on the isotopic parameters of parental kimberlite magmas.

     

Elucidating hydrological complexity and climate-influenced diagenesis during the early Carboniferous of South Wales

A revised age provided by conodonts from the Gilwern Oolite of the Clydach area of South Wales allows a clearer understanding of the palaeohydrology, palaeoclimatic history and diagenesis of previously correlated oolitic units. Earlier uncertainty over the apparent sub-regional differences in climate during the early Visean (Chadian—Arundian) lowstand is resolved. Previously a humid interval evidenced by prominent palaeo-epikarst capping the Gilwern Oolite along the northern outcrops of the South Wales synclinorium, was not recognized in what were regarded as correlative outcrops in the Chadian Gully Oolite in the southern part of the synclinorium. Conodont dating now shows that the Gilwern Oolite is much older (Courceyan) than the Chadian Gully Oolite, and during the prolonged exposure of the former there was an interval of weathering under a humid climate. This also explains the contrast between the diagenesis seen between the Gilwern and Gully oolites, which are no longer seen as correlatives.     

When hazard avoidance is not an option: lessons learned from monitoring the postdisaster Oso landslide,USA

On 22 March 2014, a massive, catastrophic landslide occurred near Oso, Washington, USA, sweeping more than 1 km across the adjacent valley flats and killing 43 people. For the following 5 weeks, hundreds of workers engaged in an exhaustive search, rescue, and recovery effort directly in the landslide runout path. These workers could not avoid the risks posed by additional large-scale slope collapses. In an effort to ensure worker safety, multiple agencies cooperated to swiftly deploy a monitoring and alerting system consisting of sensors, automated data processing and web-based display, along with defined communication protocols and clear calls to action for emergency management and search personnel. Guided by the principle that an accelerating landslide poses a greater threat than a steadily moving or stationary mass, the system was designed to detect ground motion and vibration using complementary monitoring techniques. Near real-time information was provided by continuous GPS, seismometers/geophones, and extensometers. This information was augmented by repeat-assessment techniques such as terrestrial and aerial laser scanning and time-lapse photography. Fortunately, no major additional landsliding occurred. However, we did detect small headscarp failures as well as slow movement of the remaining landslide mass with the monitoring system. This was an exceptional response situation and the lessons learned are applicable to other landslide disaster crises. They underscore the need for cogent landslide expertise and ready-to-deploy monitoring equipment, the value of using redundant monitoring techniques with distinct goals, the benefit of clearly defined communication protocols, and the importance of continued research into forecasting landslide behavior to allow timely warning.

     

Tectono-sedimentary evolution of Jurassic–Cretaceous diapiric structures: Miravete anticline,Maestrat Basin,Spain

Integration of extensive fieldwork, remote sensing mapping and 3D models from high-quality drone photographs relates tectonics and sedimentation to define the Jurassic–early Albian diapiric evolution of the N–S Miravete anticline, the NW-SE Castel de Cabra anticline and the NW-SE Cañada Vellida ridge in the Maestrat Basin (Iberian Ranges, Spain). The pre shortening diapiric structures are defined by well-exposed and unambiguous halokinetic geometries such as hooks and flaps, salt walls and collapse normal faults. These were developed on Triassic salt-bearing deposits, previously misinterpreted because they were hidden and overprinted by the Alpine shortening. The Miravete anticline grew during the Jurassic and Early Cretaceous and was rejuvenated during Cenozoic shortening. Its evolution is separated into four halokinetic stages, including the latest Alpine compression. Regionally, the well-exposed Castel de Cabra salt anticline and Cañada Vellida salt wall confirm the widespread Jurassic and Early Cretaceous diapiric evolution of the Maestrat Basin. The NE flank of the Cañada Vellida salt wall is characterized by hook patterns and by a 500-m-long thin Upper Jurassic carbonates defining an upturned flap, inferred as the roof of the salt wall before NE-directed salt extrusion. A regional E-W cross section through the Ababuj, Miravete and Cañada-Benatanduz anticlines shows typical geometries of salt-related rift basins, partly decoupled from basement faults. These structures could form a broader diapiric region still to be investigated. In this section, the Camarillas and Fortanete minibasins displayed well-developed bowl geometries at the onset of shortening. The most active period of diapiric growth in the Maestrat Basin occurred during the Early Cretaceous, which is also recorded in the Eastern Betics, Asturias and Basque-Cantabrian basins. This period coincides with the peak of eastward drift of the Iberian microplate, with speeds of 20 mm/year. The transtensional regime is interpreted to have played a role in diapiric development.