An integrated geophysical study of Vestbakken Volcanic Province, western Barents Sea continental margin, and adjacent oceanic crust

This paper describes results from a geophysical study in the Vestbakken Volcanic Province, located on the central parts of the western Barents Sea continental margin, and adjacent oceanic crust in the Norwegian-Greenland Sea. The results are derived mainly from interpretation and modeling of multichannel seismic, ocean bottom seismometer and land station data along a regional seismic profile. The resulting model shows oceanic crust in the western parts of the profile. This crust is buried by a thick Cenozoic sedimentary package. Low velocities in the bottom of this package indicate overpressure. The igneous oceanic crust shows an average thickness of 7.2 km with the thinnest crust (5–6 km) in the southwest and the thickest crust (8–9 km) close to the continent-ocean boundary (COB). The thick oceanic crust is probably related to high mantle temperatures formed by brittle weakening and shear heating along a shear system prior to continental breakup. The COB is interpreted in the central parts of the profile where the velocity structure and Bouguer anomalies change significantly. East of the COB Moho depths increase while the vertical velocity gradient decreases. Below the assumed center for Early Eocene volcanic activity the model shows increased velocities in the crust. These increased crustal velocities are interpreted to represent Early Eocene mafic feeder dykes. East of the zone of volcanoes velocities in the crust decrease and sedimentary velocities are observed at depths of more than 10 km. The amount of crustal intrusions is much lower in this area than farther west. East of the Kn?legga Fault crystalline basement velocities are brought close to the seabed. This fault marks the eastern limit of thick Cenozoic and Mesozoic packages on central parts of the western Barents Sea continental margin.     

Exploring the impact of diagenesis on (isotope) geochemical and microstructural alteration features in biogenic aragonite

For the Quaternary and Neogene, aragonitic biogenic and abiogenic carbonates are frequently exploited as archives of their environment. Conversely, pre‐Neogene aragonite is often diagenetically altered and calcite archives are studied instead. Nevertheless, the exact sequence of diagenetic processes and products is difficult to disclose from naturally altered material. Here, experiments were performed to understand biogenic aragonite alteration processes and products. Shell subsamples of the bivalve Arctica islandica were exposed to hydrothermal alteration. Thermal boundary conditions were set at 100°C, 175°C and 200°C. These comparably high temperatures were chosen to shorten experimental durations. Subsamples were exposed to different ¹⁸O‐depleted fluids for durations between two and twenty weeks. Alteration was documented using X‐ray diffraction, cathodoluminescence, fluorescence and scanning electron microscopy, as well as conventional and clumped isotope analyses. Experiments performed at 100°C show redistribution and darkening of organic matter, but lack evidence for diagenetic alteration, except in Δ₄₇ which show the effects of annealing processes. At 175°C, valves undergo significant aragonite to calcite transformation and neomorphism. The δ¹⁸O signature supports transformation via dissolution and reprecipitation, but isotopic exchange is limited by fluid migration through the subsamples. Individual growth increments in these subsamples exhibit bright orange luminescence. At 200°C, valves are fully transformed to calcite and exhibit purple‐blue luminescence with orange bands. The δ¹⁸O and Δ₄₇ signatures reveal exchange with the aqueous fluid, whereas δ¹³C remains unaltered in all experiments, indicating a carbonate‐buffered system. Clumped isotope temperatures in high‐temperature experiments show compositions in broad agreement with the measured temperature. Experimentally induced alteration patterns are comparable with individual features present in Pleistocene shells. This study represents a significant step towards sequential analysis of diagenetic features in biogenic aragonites and sheds light on reaction times and threshold limits. The limitations of a study restricted to a single test organism are acknowledged and call for refined follow‐up experiments.     

An approach for the determination of precipitation input for worst-case flood modelling

ABSTRACT

There is a lack of suitable methods for creating precipitation scenarios that can be used to realistically estimate peak discharges with very low probabilities. On the one hand, existing methods are methodically questionable when it comes to physical system boundaries. On the other hand, the spatio-temporal representativeness of precipitation patterns as system input is limited. In response, this paper proposes a method of deriving spatio-temporal precipitation patterns and presents a step towards making methodically correct estimations of infrequent floods by using a worst-case approach. A Monte Carlo approach allows for the generation of a wide range of different spatio-temporal distributions of an extreme precipitation event that can be tested with a rainfall–runoff model that generates a hydrograph for each of these distributions. Out of these numerous hydrographs and their corresponding peak discharges, the physically plausible spatio-temporal distributions that lead to the highest peak discharges are identified and can eventually be used for further investigations.

Editor A. Castellarin; Associate editor E. Volpi     

Impact of loading displacements on SLR-derived parameters and on the consistency between GNSS and SLR results

Displacements of the Earth’s surface caused by tidal and non-tidal loading forces are relevant in high-precision space geodesy. Some of the corrections are recommended by the international scientific community to be applied at the observation level, e.g., ocean tidal loading (OTL) and atmospheric tidal loading (ATL). Non-tidal displacement corrections are in general recommended not to be applied in the products of the International Earth Rotation and Reference Systems Service, in particular atmospheric non-tidal loading (ANTL), oceanic and hydrological non-tidal corrections. We assess and compare the impact of OTL, ATL and ANTL on SLR-derived parameters by reprocessing 12 years of SLR data considering and ignoring individual corrections. We show that loading displacements have an influence not only on station long-term stability, but also on geocenter coordinates, Earth Rotation Parameters, and satellite orbits. Applying the loading corrections reduces the amplitudes of annual signals in the time series of geocenter and station coordinates. The general improvement of the SLR station 3D coordinate repeatability when applying OTL, ATL and ANTL corrections are 19.5 %, 0.2 % and 3.3 % respectively, w.r.t. the solutions without loading corrections. ANTL corrections play a crucial role in the combination of optical (SLR) and microwave (GNSS, VLBI, DORIS) space geodetic observation techniques, because of the so-called Blue-Sky effect: SLR measurements can be carried out only under cloudless sky conditions—typically during high air pressure conditions, when the Earth’s crust is deformed, whereas microwave observations are weather-independent. Thus, applying the loading corrections at the observation level improves SLR-derived products as well as the consistency with microwave-based results. We assess the Blue-Sky effect on SLR stations and the consistency improvement between GNSS and SLR solutions when ANTL corrections are included. The omission of ANTL corrections may lead to inconsistencies between SLR and GNSS solutions of up to 2.5 mm for inland stations. As a result, the estimated GNSS–SLR coordinate differences correspond better to the local ties at the co-located stations when applying ANTL corrections.     

The Basic Formal Ontology as a Reference Framework for Modeling the Evolution of Administrative Units

In information systems, ontologies promise advantages such as enhanced interoperability, knowledge sharing, and integration of data sources. In this article, we show that the upper ontology Basic Formal Ontology can facilitate the modeling of an evolution of administrative units. This is demonstrated by creating a spatiotemporal ontology for the administrative units of Switzerland. The ontology tackles the problem that the geometric data is typically captured by taking snapshots at regular intervals while the thematic data is continually updated. The ontology presented merges time‐stamped geometries with a formally described history of administrative units, allowing for complex spatiotemporal queries neither standard approach would support. The resulting populated knowledge base was evaluated against a set of spatiotemporal test queries. The evaluation showed that this knowledge base supports a wide range of queries on the evolution of the administrative units of Switzerland between 1960 and 2010.     

Hydro-geochemical processes in an area with saline groundwater in lower Shire River valley, Malawi: An integrated application of hierarchical cluster and principal component analyses

Our ability to adapt to changes in groundwater quality, arising from a changing climate and/or local pressures, is dependent on comprehension of the governing controls of spatial variation in groundwater chemistry. This paper presents results of an assessment of dominant hydro-geochemical processes controlling groundwater chemical composition, using an integrated application of hierarchical cluster analysis (HCA) and principal component analysis (PCA) of a major ion dataset of groundwater from lower Shire River valley, Malawi. The area is in the southernmost part of the western section of the East African Rift System (EARS) and has localised occurrence of saline groundwater. HCA classified samples into three main clusters (C1-C3) according to their dominant chemical composition: C1 (dominant composition: Na-Cl; median TDS: 3436 mg L⁻¹), C2 (dominant composition: Na-HCO₃; median TDS: 966 mg L⁻¹) and C3 (dominant composition: alkali earths-HCO₃; median TDS: 528 mg L⁻¹). These clusters were in turn described by the principal components PC1, PC3 and PC2, respectively, resulting from the PCA. The results of the PCA and geochemical interpretation suggest that the spatial variation of groundwater quality in the area is influenced by the following processes: C3 samples result mainly from H₂CO₃ weathering of aluminosilicate minerals by percolating water supersaturated with CO₂. In addition to aluminosilicate weathering, C2 samples are influenced by the processes of cation exchange of Ca²⁺ and Mg²⁺ in the water for Na⁺ on clay minerals, and carbonate precipitation. The increase in ionic strength of C2 samples is attributed to mixing with high TDS groundwater in proximity with C2 samples. The saline/brackish C1 groundwater results from the processes of evaporation (for samples with high water table close to the Shire marshes) and dissolution of Cl⁻ and SO₄-evaporative salts followed by mineralised seep from sedimentary Karoo and Cretaceous Lupata sandstones.     

Uranium–lead ages of apatite from iron oxide ores of the Bafq District,East-Central Iran

Iron oxide–apatite (IOA) deposits, often referred to as Kiruna-type iron ore deposits, are known to have formed from the Proterozoic to the Tertiary. They are commonly associated with calc–alkaline volcanic rocks and regional- to deposit-scale metasomatic alteration. In the Bafq District in east Central Iran, economic iron oxide–apatite deposits occur within felsic volcanic tuffs and volcanosedimentary sequences of Early Cambrian age. In order to constrain the age of formation of these ores and their relationship with the Early Cambrian magmatic event, we have determined the U–Pb apatite age for five occurrences in the Bafq District. In a ²⁰⁶Pb/²³⁸U vs. ²⁰⁷Pb/²³⁵U diagram, apatite free of or poor in inclusions of other minerals plots along the Concordia between 539 and 527 Ma with four out of five samples from one deposit clustering at the upper end of this range. For this deposit, we interpret this cluster to represent the age of apatite formation, whereas the spread towards younger ages may reflect either minor Pb loss or several events of IOA formation. Apatite with inclusions of monazite (±xenotime) yields disturbed systems with inclusions having developed after formation of the iron ore–apatite deposits, possibly as late as 130–140 Ma ago. Obtained apatite ages confirms that (IOA) and the apatite-rich rocks (apatites) of the Bafq district formed coevally with the Early Cambrian magmatic (-metasomatic) events.     

High matrix interference affecting detection of PAH metabolites in bile of Atlantic hagfish (Myxine glutinosa) used for biomonitoring of deep-water oil production

The characteristic biology and wide distribution of hagfish species makes them relevant for use in pollution biomonitoring at great water depths, particularly in regions where deep-water oil production may take place. The exposure of fish to petrogenic contaminants can normally be detected from the level of polycyclic aromatic hydrocarbon (PAH) metabolites in bile fluid. Some of these metabolites are strong fluorophores, allowing analytical detection by means of simple fluorometric techniques such as fixed wavelength fluorescence (FF) and synchronous fluorescence scanning (SFS). In the present study bile from Atlantic hagfish (Myxine glutinosa) collected in pristine areas (Barents Sea and southwestern Norway) displayed strong bile fluorescence levels, suggesting the presence of PAH contaminants. However, gas-chromatography-mass spectrometry (GC-MS) analyses ruled out PAHs as the origin for this fluorescence signal. Rather, the bile of Myxine contains components resulting in unusually strong background fluorescence interfering at the wavelength pairs used for detection of PAH metabolites. Possible background for the observed matrix interference and implications for detection of PAH metabolites in hagfish is discussed.     

Modelling thermal transients from magmatic underplating—an example from the Vøring margin (NE Atlantic)

The thermal impact of several kilometre-thick magmatic underplating in the lower continental crust is studied with analytical and numerical methods. Simple analytical solutions are derived for the thermal transient in the case of an infinite depth below the underplate and also for the case of a finite depth (down to the asthenosphere). It is shown that these solutions lead to simple approximations for when the transient surface heat flow is at its maximum, what the maximum is, and for how long the transient lasts. Even though these solutions assume that the underplate is emplaced instantaneously, they are useful in the interpretation of underplating over finite time spans. A numerical scheme is suggested for the modelling of underplating that handles both short time intervals as well as long intervals. The scheme treats magmatic underplating in a mass and energy conservative manner, and it is compared against the analytical solutions. Finally, the analytical and numerical results for thermal transients are applied to a transect from the Vøring margin (NE Atlantic), with respect to various degrees of early Cenozoic magmatic intrusion. It appears that more than half of the lower crustal body (LCB) in the Vøring margin must be magmatic underplating for the vitrinite reflectance to be substantially higher than for the non-magmatic case, where the LCB is assumed to comprise Caledonian crust.     

Protoliths and phase petrology of whiteschists

Whiteschists appear in numerous high- and ultrahigh-pressure rock suites and are characterized by the mineral assemblage kyanite + talc (+-quartz or coesite). We demonstrate that whiteschist mineral assemblages are well stable up to pressures of more than 4 GPa but may already form at pressures of 0.5 GPa. The formation of whiteschists largely depends on the composition of the protolith, which requires elevated contents of Al and Mg as well as low Fe, Ca, and Na contents, as otherwise chloritoid, amphibole, feldspar, or omphacite are formed instead of kyanite or talc. Furthermore, the stability field of the whiteschist mineral assemblage strongly depends on XCO₂ and fO₂: already at low values of XCO₂, CO₂ binds Mg to carbonates strongly reducing the whiteschist stability field, whereas high fO₂ enlarges the stability field and stabilizes yoderite. Thus, the scarcity of whiteschist is not necessarily due to unusual P–T conditions, but to the restricted range of suitable protolith compositions and the spatial distribution of these protoliths: (1) continental sedimentary rocks and (2) hydrothermally and metasomatically altered felsic to mafic rocks. The continental sedimentary rocks that may produce whiteschist mineral assemblages typically have been deposited under arid climatic conditions in closed evaporitic basins and may be restricted to relatively low latitudes. These rocks often contain large amounts of the clay minerals palygorskite and sepiolite. Marine sediments generally do not yield whiteschist mineral assemblages as marine shales commonly have too high iron contents. Sabkha deposits may have too high CO₂ contents. Protoliths of appropriate geochemical composition occur in and on continental crust. Therefore, whiteschist assemblages typically are only found in settings of continental collision or where continental fragments were involved in subduction. Our calculations demonstrate that whiteschists can form by closed-system metamorphism, which implies that the chemical and isotopic composition of these rocks provide constraints on the development of the protoliths.