Beach-ridge formation as a possible indicator for an open Limfjord – North Sea connection

Raised beach ridges are prograded sequences of wave-built deposits that may provide valuable information about past relative sea-level changes, climate change and coastal evolution. In the Limfjord in northern Denmark, the Early and Middle Holocene sea-level changes are well-constrained. However, our understanding of Late Holocene sea-level fluctuations is limited, and the exact period when the coastal barrier between the Limfjord and the North Sea formed remains uncertain. In this study, we use optically stimulated luminescence (OSL) dating to determine the age of raised beach ridges at Gjellerodde in the western part of the Limfjord. The OSL ages presented here indicate that the beach ridges formed during three periods at 3.3–2.7, 1.4–1.0, 0.2–0.1 ka. In addition our data suggest a c. 0.2 mm/yr relative sea-level fall during the Late Holocene. The three distinct periods of beach-ridge formation coincide with periods when the Limfjord was open towards the North Sea as documented in historical records and marine records. This suggests that OSL dating of beach ridges can be used as a potential indicator for determining when the connection between the Limfjord and the North Sea was open in the Late Holocene

Preface

The Palaeogene Skaergaard intrusion, East Greenland, has since its first description by Wager & Deer (1939) been a foremost natural laboratory for the study of low-pressure fractionation of basaltic melt. Ocean floors are composed of basalt and the processes that control compositions of basaltic melts are fundamental to the dynamics of the Earth. This special issue of GEUS Bulletin by Peter Thy, Christian Tegner and Charles E. Lesher is the most recent in a more than eighty year succession of trendsetting works on the evolution of the Skaergaard intrusion and evolution of basaltic melt. The early sample collections, housed in the universities of Oxford and Cambridge, were used for the development of fractionation models, mapping of the distribution as well as the partition of major and trace elements between melts and liquidus phases in basaltic melts. These and many other studies resulted in the monumental “Layered Igneous Rocks”, edited by L.R. Wager and G.M. Brown (Wager & Brown 1968). By the 1970s, the original collection in the Oxford University Museum of Natural History and the Sedgwick Museum of Earth Sciences, University of Cambridge had seen extensive use. The science developed and new and more detailed sampling was now required. Fieldwork and sampling in the later part of the 20th century led to a flurry of new studies by more research groups, including Neil Irvine of Carnegie’s Geophysical Laboratory and the University of Oregon group lead by Alex McBirney. Despite all the efforts, no consensus was reached on the many fundamental chemical and physical processes in basaltic magma chambers. The available collections did not provide sufficiently systematic and detailed information for the modelling of the evolution of the basaltic melt and the genesis of the precious metal deposit that had been discovered in the intrusion. Exploration drill cores offered the possibility for tight stratigraphic sampling through >1000 m of gabbro. In combination with surface sampling, a new standard profile with superior systematics could be established. In the summer of 2000, the Danish Lithosphere Centre (DLC) organised fieldwork and the transport of the drill cores to the Museum of Natural History in Copenhagen. This special issue is part of the legacy of the DLC and a fulfilment of the then envisaged potential of the Skaergaard intrusion for fundamental studies of processes in basaltic magma chambers. Thy et al. (2023, this volume) analyse and compile in unprecedented detail the variation in major and trace element compositions of bulk-rock samples and electron microprobe compositions of rock-forming minerals in the Layered Series on the floor of the intrusion. In addition, they include a wealth of information on all samples from the Skaergaard intrusion that are housed in museum, university and survey collections in Denmark, the UK and US, topographic and geological maps, aerial photographs in archives and images of thin sections used in the study. With all its data and appendices, the monograph is the most comprehensive collection of petrological information for the Skaergaard intrusion and the Layered Series on the floor of the intrusion, ever to have been made available for research. The detailed modelling presented re-evaluates petrogenetic constraints and tests petrogenetic models in the literature. The modelling is based on liquidus proportions established in experimental studies of appropriate melt compositions. Thy and co-authors conclude that evolution observed in the gabbros of the Layered Series is predominantly controlled by crystal fractionation of bulk liquid. Late in the evolution and above the boundary between Upper Zone a and b, the gabbros crystallised from ponded Fe-rich, immiscible melt. Redistribution of granophyric melts and decoupling of included and excluded trace elements, not only in bulk liquid but also in mush melts, complicates the modelling of trace element distributions in the late and evolved stages of crystallisation in the crystal mushes of the intrusion. This special issue is a treasure trove for all who study the intrusion and are searching for a data set to be matched with experimental and numeric modelling of low-pressure fractionation of basaltic melt. With this contribution, the Skaergaard intrusion continues to be a foremost natural laboratory for modelling of melts formed in solidifying basaltic magma chambers. Troels F.D. NielsenEmeritus, GEUS, Denmar

Scaling the Danish national water resources model for a pan-European quasi-3D groundwater resources model

In this study, we upscale and simplify hydrostratigraphic information from a detailed model for Denmark to a pan-European scale. This is part of a larger project to develop a harmonised overview of the volume and depth of groundwater resources in a quasi-3D European groundwater resource model. A 10 km grid and a maximum of c. 10 hydrostratigraphic layers were chosen as the common scale for the European database. The Danish information is based on the national water resources model (the DK-model), where the information is significantly more detailed (100 m grid and up to 26 layers). Information was transferred from the DK-model to the quasi-3D model by a method involving computations of mean volumes and expert assessment to reduce layers in each cell. In this process, detailed hydrostratigraphic information is lost, which could otherwise be used for local groundwater flow modelling in Denmark. However, the strength of the quasi-3D model is that it still contains the volumes of all hydrostratigraphic units, both the saturated and unsaturated parts. Hence, the upscaled model can contribute to a relatively precise calculation of European groundwater resources for the quantitative assessment of groundwater status across Europe at a 10 × 10 km scale

Mudstone diagenesis and sandstone provenance in an Upper Jurassic – Lower Cretaceous evolving half-graben system, Wollaston Forland, North-East Greenland

The influence of rifting on the composition of Kimmeridgian to Barremian mudstones from northern Wollaston Forland, North-East Greenland is investigated by petrographic and mineralogical analyses of the Brorson Halvø-1 and Rødryg­gen-1 cores, and provenance of analysis by zircon U-Pb age dating of nearby sandstones. Mudstone composition varies systematically as a function of the timing of rifting progression and position in the half-graben depositional system. Pyrite primarily precipitated in the early rift to rift climax phases. Euhedral pyrite overgrowths on framboids formed only during the rift climax phase (Lindemans Bugt Formation). Dolomite is the dominant carbonate cement, except for the sediments deposited in the early waning rift phase (Palnatokes Bjerg Formation) where calcite is dominant, and in the late waning rift phase (Stratumbjerg Formation) where siderite dominates. The highest-temperature reactions with precipitation of illite, quartz, ankerite and barite signify sediment burial depths of >2 km prior to exhumation. Uplift-induced fracturing occurred mainly in the early rift to rift acceleration succession (Bernbjerg Formation). Mudstones in the proximal part of the half-graben (Rødryggen-1) include more detrital kaolinite than the distal mudstones (Brorson Halvø-1), which contain more mixed-layer illite-smectite and illite. Vermiculite was deposited only in the proximal part of the basin in the rift climax and waning rift successions. Chlorite was deposited proximally and distally during the waning rift phase, though supply began earlier in the distal part. Fine-grained sediment in the distal part of the half-graben was therefore probably supplied by axial transport from Palaeoproterozoic crystalline rocks and Meso- to Neoproterozoic metamorphic rocks located to the north and north-west. This agrees with the zircon provenance signature from outcropping sand-rich facies, where zircon grains with U-Pb ages of 2.0–1.6 Ga are dominant, in addition to common 1.6–0.9 Ga ages, and fewer 2.8–2.6 Ga and 0.47–0.36 Ga ages

The PGE-Au Mineralisation of the Skaergaard intrusion: precious metal minerals, petrography and ore genesis

The Skaergaard PGE-Au Mineralisation, alias the Platinova Reef, is hosted in a series of mineralisation levels within a suite of bowl-shaped macrorhythmic layers in the upper Middle Zone of the Skaergaard intrusion. The intrusion is exposed 68°N in East Greenland. The occurrence defines its own type due to its exceptional structure and mineralogy. A wealth of mineralogical data is available in laboratory reports for individual samples and in peer-reviewed publications, but none of these account for the lateral and stratigraphic distribution of PGE and Au parageneses in the gabbros of the intrusion. In this study, we collate and describe the mineralogical data for the first-formed PGE-rich and last-formed gold-rich mineralisation levels and integrate these with petrogenetic models. Recovery of >4000 grains of precious metal phases allow a detailed study of their distribution and compositions throughout the mineralisation, re-equilibration during cooling, inter-grain relationships and relationships to Cu-Fe sulphides and the gabbroic host rocks. The sulphides are dominated by bornite, chalcocite and minor chalcopyrite. All other sulphides, such as pentlandite, are very rare. Fifty-four different precious metal phases are identified in this study, and include the new IMA approved minerals skaergaardite (PdCu), nielsenite (Pd3Pb) and naldrettite (Pd2Sb). Precious metal phases include (1) intermetallic compounds and alloys of Cu and Pd; (2) intermetallic compounds and alloys of Au and Cu (Ag); (3) sulphides of Pd, Cu (Ag, Cd, Hg, Tl); (4) arsenides of Pd (Pt, Ni) and (5) intermetallic compounds of Pd, Cu with Sn, Pb, Te (Sb, Bi). Skaergaardite (PdCu) is the dominant PGE mineral in the lower and main PGE mineralisation level (Pd5). It is accompanied at the western margin of the intrusions by the sulphides vasilite (Pd16S7) and vysotskite (PdS) but is rare at the eastern margin, which is dominated by plumbide zvyagintsevite (Pd3Pb). Gold phases include a suite of intermetallic compounds and alloys from AuCu3 to native gold and are dominated by tetra-auricupride (AuCu). Gold is concentrated in the tops of individual mineralisation levels and in the uppermost precious metal–bearing mineralisation level, followed by stratiform Cu-rich mineralisation levels. Precious metal parageneses demonstrate formation and re-equilibration from liquidus to subsolidus temperatures and control by local geochemical environments. The mineralisation is syn-magmatic and the result of fractionation and evolution in the remaining bulk-silicate liquid and crystal mushes. Fractionation led to sulphide saturation and formation of immiscible sulphide melt droplets. This was followed by reaction with mush melts and re-equilibration to lower temperatures, first under the roof and subsequently after slumping to the floor in mushes of macrorhythmic layers. Droplets of sulphide melt formed between 1030–1050°C and trapped precious metals. The subsequent reaction between sulphide melt and interstitial Fe-rich immiscible melt at c. 1015°C, and redistribution to coexisting melt and fluid, led to the separation of PGE, Au and Cu and their up- and inward transport. Magmatic fluids as well as volatile-rich residual silicate melts were retained in gabbros at the margins and resulted in precious metal parageneses in equilibrium with hydrous low-temperature silicate parageneses

Preface

In this special issue of GEUS Bulletin, the many riddles regarding the platinum group elements and gold (PGE-Au) mineralisation of the East Greenland Skaergaard intrusion are untangled and discussed. The Skaergard PGE-Au mineralisation, as defined in this study, embodies an enigmatic and rich ore-formation that arguably could have been an economic resource, had it not been for its ice-locked position in central East Greenland. The authors of this study (Rudashevsky et al. 2023, this volume) characterise the systematic variability in the precious metal mineralogy from the contact towards the interior of the intrusion based on the analysis of more than 4000 individual PGE-Au grains. This variability is interpreted in the light of 90 years of research and over 1000 publications pertaining to magma chamber processes in the Skaergaard intrusion.  With such an impressive library of knowledge, on a comparatively simple magmatic system such as the Skaergard intrusion, we should have discovered a few islands of truth in igneous petrology and ore-deposit formation. And indeed, we have. But we are also enriched with an evolving story, where answering one question only serves to raise three new questions.  This study demonstrates the variability of PGE-Au phases throughout the ore-forming zone of the Skaergaard intrusion. As previously observed, PGE-Au mineralisation in the central parts is divided into several layers over 30–40 m of the cumulus stratigraphy with increasing Pd/Pt ratios upwards, an Au-rich upper part and a low sulphide content throughout all layers. Close to the contact the precious metal zonation is less pronounced, and it is significantly more sulphide rich. The PGE mineralogy deviates significantly from the centre to the margin. These complex lateral and vertical variations cannot be explained by one genetic model but require an intricate combination of igneous processes including silicate-melt liquid immiscibility, sulphide-melt immiscibility, sulphide-melt resorptions, precious metal transport by volatile-rich fluids and, finally, the solidification rate of the cumulus mushes.  For other well-preserved PGE-Au deposits throughout the world, we observe a great variation of ore-forming models. Remarkably, most of these models may be applied to various parts of the Skaergaard mineralisation. The authors suggest that the Skaergaard intrusion preserves different steps in PGE-Au ore-genesis which, in many other intrusions are obliterated by later igneous events. Therefore, the legacy of the Skaergaard mineralisation is the preservation of igneous ore-forming events that may also precede the genesis of other PGE-Au deposits in the world. After the last conclusion, I guarantee that you will be confused and perhaps a bit triggered but hopefully also inspired and bursting with new questions on the genesis of PGE-Au deposits in mafic and ultramafic igneous complexes. In light of the recent study, you may even be encouraged to look at your favourite PGE-Au deposit with fresh eyes. This study beautifully demonstrates that turning the next page in the book of magma chamber processes is more important than seeing ‘The End’

Machine learning-based estimation and clustering of statistics within stratigraphic models as exemplified in Denmark

Estimating a covariance model for kriging purposes is traditionally done using semivariogram analyses, where an empirical semivariogram is calculated, and a chosen semivariogram model, usually defined by a sill and a range, is fitted. We demonstrate that a convolutional neural network can estimate such a semivariogram model with comparable accuracy and precision by training it to recognise the relationship between realisations of Gaussian random fields and the sill and range values that define it, for a Gaussian type semivariance model. We do this by training the network with synthetic data consisting of many such realisations with the sill and range as the target variables. Because training takes time, the method is best suited for cases where many models need to be estimated since the actual estimation itself is about 70 times faster with the neural network than with the traditional approach. We demonstrate the viability of the method in three ways: (1) we test the model’s performance on the validation data, (2) we do a test where we compare the model to the traditional approach and (3) we show an example of an actual application of the method using the Danish national hydrostratigraphic model

Organic geochemistry of an Upper Jurassic – Lower Cretaceous mudstone succession in a narrow graben setting, Wollaston Forland Basin, North-East Greenland

The Oxfordian–Ryazanian was a period of widespread deposition of marine organic-rich mudstones in basins formed during the early phases of the rifting that heralded the formation of the present-day North Atlantic. Occasionally, uninterrupted deposition prevailed for 20 million years or more. Today, mudstones of this time interval are found on the shelves bordering the North Atlantic and adjacent areas from Siberia to the Netherlands. Here, we report data on two fully cored boreholes from Wollaston Forland (North-East Greenland, approx. 74° N), which represent an uninterrupted succession from the upper Kimmeridgian to the Hauterivian. The boreholes record basin development at two different positions within an evolving halfgraben, located at the margin of the main rift, and thus partially detached from it. Although the overall depositional environment remained an oxygen-restricted deep-shelf setting, rifting-related changes can be followed through the succession. The Kimmeridgian was a period of eustatic highstand and records the incipient rifting with a transgressive trend straddling the transition to the lower Volgian by a gradual change from deposits with high levels of total organic carbon (TOC) and kerogen rich in allochthonous organic matter to deposits with lower TOC and a higher proportion of autochthonous organic matter. This is followed by a slight regressive trend with lower TOC and increased proportions of allochthonous organic matter until rifting culminated in the middle Volgian–Ryazanian, indicated by increasing autochthonous organic matter and higher TOC, which prevailed until basin ventilation occurred towards the end of the Ryazanian. The properties of the reactive kerogen fraction remained rather stable irrespective of TOC, underlining the effect of terrigenous matter input for TOC. These variations are also captured by biological markers and stable carbon isotopes. The deposits are very similar to equivalent successions elsewhere in the proto-North Atlantic region, albeit the proportion of terrigenous kerogen is greater

The Rødryggen-1 and Brorson Halvø-1 fully cored boreholes (Upper Jurassic – Lower Cretaceous), Wollaston Forland, North-East Greenland – an introduction

Two fully cored boreholes, the Rødryggen-1 and the Brorson Halvø-1, were drilled in Wollaston Forland, North-East Greenland, in 2009 and 2010, respectively. The objective was to test the stratigraphic development of the Upper Jurassic – Lower Cretaceous mud-dominated succession in two different settings within the same fault block of a developing half-graben: centrally (Rødryggen-1 borehole) and near the uplifted crest of the rotating fault block (Brorson Halvø-1 borehole). The drilled deposits are equivalent to the principal petroleum source-rock sequence of the petroliferous basins of North-West Europe, Siberia, and basins off eastern Canada and provide a new record of an important phase of marine deoxygenation in the proto-North Atlantic region.

Paleo sea-level indicators and proxies from Greenland in the GAPSLIP database and comparison with modelled sea level from the PaleoMIST ice-sheet reconstruction

One of the most common ways to assess ice-sheet reconstructions of the past is to evaluate how they impact changes in sea level through glacial isostatic adjustment. PaleoMIST 1.0, a preliminary reconstruction of topography and ice sheets during the past 80 000 years, was created without a rigorous comparison with past sea-level indicators and proxies in Greenland. The basal shear stress values for the Greenland ice sheet were deduced from the present day ice-sheet configuration, which were used for the entire 80 000 years without modification. The margin chronology was based on previous reconstructions and interpolation between them. As a result, it was not known if the Greenland component was representative of its ice-sheet history. In this study, I compile sea–level proxy data into the Global Archive of Paleo Sea Level Indicators and Proxies (GAPSLIP) database and use them to evaluate the PaleoMIST 1.0 reconstruction. The Last Glacial Maximum (c.20 000 years before present) contribution to sea level in PaleoMIST 1.0 is about 3.5 m, intermediate of other reconstructions of the Greenland ice sheet. The results of the data-model comparison show that PaleoMIST requires a larger pre-Holocene ice volume than it currently has to match the sea-level highstands observed around Greenland, especially in southern Greenland. Some of this mismatch is likely because of the crude 2500 year time step used in the margin reconstruction and the limited Last Glacial Maximum extent. Much of the mismatch can also be mitigated if different Earth model structures, particularly a thinner lithosphere, are assumed. Additional ice in Greenland would contribute to increasing the 3–5 m mismatch between the modelled far-field sea level at the Last Glacial Maximum and proxies in PaleoMIST 1.0