The volcanic succession of the Paleogene Faroe Islands Basalt Group (FIBG; Passey & Jolley Reference Passey and Jolley2009) comprises three dominantly basaltic lava flow formations separated by sedimentary rock formations (Fig. 1). The oldest lava flow succession, the Beinisvørð Formation, comprises a thick succession of simple, sheet-like, lava flows (Passey & Bell Reference Passey and Bell2007; Passey & Jolley Reference Passey and Jolley2009). The exposed succession shows an increase in the dominance of columnar jointing towards the top, including thick multi-tiered columnar jointed flows in central east Suðuroy. The Malinstindur Formation is a thick succession of dominantly compound lava flows which is separated from the Beinisvørð Formation by two volcano–sedimentary rock formations (Passey & Bell Reference Passey and Bell2007; Passey & Jolley Reference Passey and Jolley2009). The oldest of these, the Prestfjall Formation, is composed of carbonaceous and dominantly argillaceous sedimentary rocks (Passey & Jolley Reference Passey and Jolley2009) which crop out most extensively on the island of Suðuroy (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970). The formation is generally poorly exposed and has been mainly characterised by observations made from the coal mines on Suðuroy (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970). To the north, the Prestfjall Formation crops out on the islands of Tindhólmur and western Vágar (Fig. 1), but is of a notably more arenaceous character. The lower boundary of the Prestfjall Formation has been interpreted to mark a significant hiatus in volcanic activity following the cessation of Beinisvørð Formation eruptions (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970). Overlying the Prestfjall Formation, the volcaniclastic rocks of the Hvannhagi Formation crop out in eastern Suðuroy, at isolated locations in north-west Suðuroy, west Vágar, and on Tindhólmur. The formation comprises both pyroclastic and volcaniclastic sedimentary rocks that were deposited in a terrestrial environment (Passey & Jolley Reference Passey and Jolley2009), and has been regarded as marking the resumption of large-scale volcanism in the region (Passey Reference Passey2004).
Recent analysis of upper Beinisvørð Formation lava flow geochemistry (Millett et al. Reference Millett, Hole, Jolley, Passey and Rossetti2020) has provided evidence for both an increase in the overall degree of partial melting and a reduction in the initial pressure of melting. An increase in partial melting was identified from data recorded in lava flows from the upper Beinisvørð Formation. This was interpreted as evidence for the early onset of northeast Atlantic rifting. Millett et al. (2020) also proposed a linkage between the reduction in initial pressure of melting, a transient reduction in mantle temperature, and the North Atlantic Igneous Province (NAIP) wide volcanic hiatus. This was correlated with a proposed increase in thickness and duration of interbedded sedimentary rock units prior to cessation of volcanism and deposition of the carbonaceous Prestfjall Formation.
Attributed to the transition from pre-rift to syn-rift (Larsen et al. Reference Larsen, Waagstein, Pedersen and Storey1999; Millett et al. Reference Millett, Hole, Jolley, Passey and Rossetti2020), FIBG strata between the extrusive successions of the Beinisvørð and Malinstindur formations are stratigraphically complex and spatially variable. This study contributes to the understanding of this complexity which will help to answer questions regarding the rate of slowing, cessation, and resumption of eruption at a period of major reorganisation of the NAIP.
1. Methods
1.1 Photogrammetry and field data collection
Primary lithological records were collected from outcrops in the field. Field transects and spot locations were recorded (Vágar and Suðuroy) and borehole core was logged from Suðuroy (KOLBH1-04) and Vágar (Gásadalur 1989-1). Many of these are detailed in Passey & Jolley (Reference Passey and Jolley2009), but are supplemented here with additional data. Lithology, structure, and thickness of rock units were used to construct maps, with photogrammetry interpretation undertaken on models for cliff sections from Vágar and Tindhólmur. The photogrammetry models were prepared following the workflow described in Passey et al. (Reference Passey, Vosgerau, Sørensen and Vautravers2016), Vosgerau et al. (Reference Vosgerau, Passey, Svennevig, Strunck, Jolley, Bowman, Smyth, Good, Passey, Hirst and Jordan2016), and Sørensen & Dueholm (Reference Sørensen and Dueholm2018). Briefly summarised, stereo-images (overlapping images) of the coastal cliffs were acquired using a single handheld digital camera deployed on a helicopter flying with an open door. The camera was pointed approximately perpendicular to the slope of the cliffs and the pilot was instructed to fly parallel to the cliffs. The images were prepared for analysis using the Agisoft software Metashape for tie-point generation. The triangulation or bundle adjustment of the images were then undertaken in the Anchorlab software three-dimensional (3D) Stereo Blend using a combination of the tie-points, Global Positioning System data, sea-level levering points, and control points taken from the Faroes Environment Agency online portal (www.kortal.fo). The orientation of the images (results of the triangulations) was exported from 3D Stereo Blend and used to generate 3D mesh models of the cliff phases using nFrames software SURE. These, which in the following will be referred to as the photogrammetry models, were subsequently interpreted using Lime v2.2.1. (Buckley et al. Reference Buckley, Ringdal, Naumann, Dolva, Kurz, Howell and Dewez2019). Photogrammetry models were interpreted with cross-reference to air, boat, and shore photographs to clarify details. The models were interpreted using colour-coded lines to characterise different volcanic and sedimentary beds. More laterally extensive surfaces (e.g., large flow unit tops and unconformities) were highlighted to allow for tracing depositional packages. These extensive surfaces together with some indicative bed junctions were then exported and are presented here. Surface co-ordinates for significant facies changes and unconformities were exported into a spreadsheet and then to Surfer10 to create surface and isochron maps.
1.2 Palynology
Samples were taken for palynological analysis from shales, coals, and poorly sorted sandstone units in field sections and from cored boreholes. These were processed following standard techniques, including hydrofluoric acid digestion, boiling in 20% hydrochloric acid to remove precipitates, and oxidation for 5 min in dilute 40% nitric acid where necessary. Coal samples were prepared by dissolution using fuming nitric acid. The resultant >7 μm residues were mounted in a permanent petropoxy mounting medium and examined under a transmitted light Olympus BX53 microscope. For each sample, counts of 250 specimens were targeted, but with the highly variable recovery from interlava field sedimentary rocks, frequently not attained. Accordingly, data are normalised as percentages in sections where larger counts were obtained, and as square roots in more variable recovery sections.
Ecological analysis of palynological data was used to underpin depositional environment interpretations of sedimentary rocks. Many terrestrial sedimentary sequences demonstrate progression in plant community seral succession through the duration of their deposition; these stages of community ecology are referred to as early, mid, and late successions (e.g., Vitousek Reference Vitousek2004; Jolley et al. Reference Jolley, Bell, Williamson and Prince2009). Within the NAIP, newly erupted lava flow surfaces and volcaniclastic sediment deposits would, with the passage of time, begin to be colonised by early successional plant communities. These were dominated by disturbance and stress tolerant species (e.g., Equisetum, ‘mares’ tails’) and in humid environments, ferns, and mosses. In the absence of disturbance, these communities would have diversified to form early–mid successional communities (Whittaker et al. Reference Whittaker, Bush and Richards1989; Thornton Reference Thornton and Sigurdsson2000). Palynofloras derived from these communities were dominated by polypodiaceous fern spores (e.g., Laevigatosporites haardtii), but include moss and liverwort spores, Ginkgo pollen (Monocolpopollenites tranquilus), and walnut/hickory pollen (Juglandaceae; Caryapollenites, Momipites and Platycaryapollenites species). In Miocene and younger volcanic terrains, primary succession has involved grasses and herbaceous angiosperms as key early colonists (Whitaker et al. Reference Whittaker, Bush and Richards1989; Thornton Reference Thornton and Sigurdsson2000). Evidence indicates that this niche was dominated by ferns in the Paleogene (Jolley Reference Jolley and Widdowson1997; Jolley et al. Reference Jolley, Bell, Williamson and Prince2009).
Continuing substrate stability over time would have led to the development of mid-successional vegetation. Palynofloras derived from these communities are higher in diversity, and are typified by Fagaceae (chestnut types), Platanaceae (plane types) Juglandaceae (walnut and hickory types) and Myricaceae (myrtles and hazel types). Many of these taxa occur within other communities, but at lower frequencies. Fagaceous pollen attributed to Cupuliferoipollenites, which is one of the dominant taxa in this grouping, occur in a wide variety of assemblages but decline in significance in late successional vegetation.
Late successional vegetation within the NAIP was most often dominated by a Cupressaceae–Nyssaceae association (swamp cypress and black gum), which is of low diversity (e.g., Jolley et al. Reference Jolley, Bell, Williamson and Prince2009, Reference Jolley, Passey, Hole and Millett2012). This community often formed the climax vegetation. A number of species important in early–mid and mid successional vegetation decline in importance and diversity in this community as they were eliminated by effective competition from Cupressaceae–Nyssaceae mire species. Patchworks of disturbed areas within this late-successional community are inferred from the incursion of other species; these include the taxon Caryapollenites, which is associated with mid-successional mires (Ellis et al. Reference Ellis, Jolley, Passey and Bell2009; Jolley et al. Reference Jolley, Bell, Williamson and Prince2009). Other such taxa, including Monocolpopollenites tranquilus (Ginkgo), Retitricolpites retiformis (Platanaceae) and Tricolpites cf. hians (Platanaceae), have previously been noted as part of riparian vegetation groupings (Streigler Reference Striegler1990; Jolley Reference Jolley and Widdowson1997), characteristic of disturbed fluvial margin environments. Accompanying these species was a range of taxa of lesser significance that may represent disturbed vegetation or understory.
Distributed across the range of environments is the betulaceous grain Alnipollenites verus. This taxon has been suggested as a nitrogen-fixing early colonist in nutrient depleted soils (Chapin et al. Reference Chapin, Walker, Fastie and Sharman1994; Hobbie et al. Reference Hobbie, Macko and Shugart1998; Jolley et al. Reference Jolley, Widdowson and Self2008). Its occurrence suggests that the soils were nutrient deficient across a range of profiles.
In addition to pollen and spores from higher plants, the palynofloras recovered from the FIBG include green algae, acritarchs, and dinoflagellate cysts. Chlorophycean algae occur commonly in samples from lava field sedimentary rock interbeds, where they are indicative of bodies of low-energy fresh to weakly brackish water. Small acanthomorph acritarchs are recovered frequently from tidally influenced fluvial deposits (Jolley et al. Reference Jolley, Passey, Hole and Millett2012, Reference Jolley, Millett, Schofield, Broadley and Hole2021), where they probably tolerated weakly saline water masses. Dinoflagellate cysts have been recorded from a range of normal salinity marine and brackish water environments in the Paleogene (e.g., Bujak et al. Reference Bujak, Downie, Eaton and Williams1980; Pross & Brinkhuis Reference Pross and Brinkhuis2005; Vieira & Jolley Reference Vieira and Jolley2020), following similar distribution patterns to those of extant species. However, the small size of these palynomorphs means that they behave as silt size sedimentary particles in depositional systems. This taphonomic factor is reflected in their low frequency occurrence within the tidal limit of estuarine systems.
The palynological taxonomy used in this study is that of dinoflaj3 (dinoflagellate cysts), Tappan (Reference Tappan1980; other algae), and Jolley & Morton (Reference Jolley and Morton1992; pollen and spores).
1.3 Whole rock geochemistry
A limited study of major and trace element geochemistry of the cored sedimentary interbeds was conducted using X-ray fluorescence (XRF) and flame atomic absorption photospectrometry (FAAS). Samples collected for palynology were analysed for major element geochemistry with some trace element analysis. Standard XRF techniques were used to analyse the range of major oxides and trace elements present (see Supplementary data S1 and S2; available at: https://doi.org/10.1017/tre2200005). These data are presented either as element oxide ratios or as major element weathering indices (Sheldon & Tabor, Reference Sheldon and Tabor2009). Ratios and indices are correlated to palynology-derived plant community ecology and changes in the depositional system to highlight links to depositional system dynamics and plant nutrient uptake (Vitousek, Reference Vitousek2004; Jolley et al., Reference Jolley, Widdowson and Self2008, Reference Jolley, Passey, Hole and Millett2012).
2. Pre-rift to syn-rift rock succession on Suðuroy
The Prestfjall and Hvannhagi formations of the island of Suðuroy have received the majority of attention from past scientific studies. This was principally because of the lithological evidence that the carbonaceous sediments presented for a major hiatus in FIBG eruption, but also for the resource implications of the coal beds. The current study initially focusses on the Prestfjall and Hvannhagi formations of Suðuroy within the context of this published research. This research is expanded by the addition of palynological and contextual lithological evidence from previously unstudied exposures and boreholes. Integration of these new data with the published record has allowed understanding of the processes controlling deposition of both the Prestfjall and Hvannhagi formations in Suðuroy. From this data-rich geographical area, the current study extends the understanding of these important interlava units into the little studied northerly exposures of Tindhólmur and Vágar.
2.1. Prestfjall Formation
On Suðuroy, the Prestfjall Formation covers an area of ~23 km2 and where complete, ranges in thickness between 3 m and 15 m with an average thickness of ~9 m (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970). The Beinisvørð–Prestfjall Unconformity (Passey & Varming Reference Passey and Varming2010) exhibits some incision into the uppermost flows of the Beinisvørð Formation and is overlain by the sedimentary rocks of the Prestfjall Formation (Passey & Jolley Reference Passey and Jolley2009). The Prestfjall Formation (Fig. 2d) is usually only accessible and exposed in completeness in mine workings, which except for one on the east side of the Rokhagi valley are now all closed. Fortunately, a complete description of the Prestfjall Formation derived from now closed mine workings and test pits was made by Rasmussen & Noe-Nygaard (Reference Rasmussen and Noe-Nygaard1970), who described 41 separate profiles from the island of Suðuroy (Supplementary Data S3).
These authors also subdivided the Prestfjall Formation into 5–6 informal beds based on the position of two commonly occurring coal seams (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970). The combined coal seam thickness rarely exceeded 1 m, having an average thickness of between 70 cm and 80 cm (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970; Lund Reference Lund1989), although an exceptional thickness of ~1.7 m was observed locally (Lund Reference Lund1989). The base of the Prestfjall Formation is commonly composed of a whitish-grey ‘underclay’ that has a maximum thickness of ~1 m (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970). However, in profiles examined here, KOLBH1-04 borehole and Holið í Helli (Fig. 1), the ‘underclay’ is represented by a sideritic claystone succession with a thickness of 4 m–10 m (Passey Reference Passey2014).
The ‘underclay’ bed is usually overlain by the ‘lower coal’, which has an average thickness of between 38 cm and 45 cm, and a maximum thickness that rarely exceeds 1 m (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970). The coal of the lower coal seam is usually composed of alternating layers of dull (durain) and bright (vitrain) coal. The lower coal seam is commonly overlain by a dark shale that has an average thickness of ~24 cm, and which frequently contains streaks and lenses of bright (vitrain) coal (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970; Passey Reference Passey2014). The upper coal seam overlies the dark shale and has an average thickness of ~38 cm, and a maximum thickness in the region of ~80 cm (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970; Passey Reference Passey2004). The coal seam is commonly composed of bright (vitrain) coal.
The ‘upper coal’ is generally overlain by the ‘roof clay’ that reaches a maximum observed thickness of ~13 m (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970), and is locally termed ‘ranin’. The lower part of the roof clay sometimes contains streaks or lenses of coal (Passey Reference Passey2014). The upper part is locally a reddish claystone unit, informally referred to as ‘takleir’. At several locations in western Suðuroy the ‘ranin’ roof clay is incised by channels containing greenish-brown volcaniclastic sandstones and granule-grade conglomerates with an observed maximum thickness of ~4 m (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970; Passey Reference Passey2004). The volcaniclastic sandstones are composed of palagonitised basaltic glass with a large proportion, particularly in the conglomerates, of finely crystalline basalt clasts of various compositions and textures.
2.2. Hvannhagi Formation
The lower rock units of the Hvannhagi Formation on Suðuroy comprise a sequence of air-fall tuffs and lapillistone fall deposits at least 25 m thick, or a succession of volcaniclastic conglomerates up to 8 m in thickness (Passey & Jolley Reference Passey and Jolley2009). An accurate thickness for these units is difficult to obtain because of thick, irregular dolerite sills, particularly in the northeast of Suðuroy (Fig. 2a, c). The pyroclastic sequence of greyish to red (where weathered) coarse tuffs and lapillistones are overlain by a 3 m–17 m thick sequence of volcaniclastic mudstones, sandstones, and conglomerates. These volcaniclastic sedimentary rocks are immature and are predominantly composed of reworked palagonitised basaltic glass clasts with rare clasts of finely crystalline basalt. Occasionally, reworked clasts of the underlying greyish olivine-phyric coarse tuff are observed within this sequence. Locally, this 3 m–17 m thick volcaniclastic sedimentary section contains a heterolithic unit of thinly-bedded sandstones and mudstones up to 5.5 m thick (Fig. 2b).
This volcaniclastic sedimentary succession is overlain by a 1 m–4 m thick greyish olivine-phyric welded tuff (Passey & Jolley Reference Passey and Jolley2009). The coarse tuff is homogeneous and contains abundant elongated dark brown flattened fiamme, which define a planar foliation or eutaxitic texture (Passey Reference Passey2004). The tuff is, in turn, overlain by a 1 m–9 m thick sequence of poorly sorted volcaniclastic sedimentary conglomerates and sandstones, composed dominantly of reworked palagonitised basaltic glass clasts. Other clasts are composed of coal, volcaniclastic mudstone, finely crystalline basalt, and olivine-phyric tuff.
2.3 New palynological observations on the Prestfjall and Hvannhagi formations, Suðurøy
Records of the Prestfjall Formation succession in exposures and mine workings were made by Rasmussen & Noe-Nygaard (Reference Rasmussen and Noe-Nygaard1970); these focussed principally on their value as a coal resource. These records have been supplemented by additional field sections and a borehole (KOLBH1-04, Fig. 1c).
2.3.1. Holið í Helli profile
The Holið í Helli profile in northeast Suðuroy (Fig. 3) was described in detail by Passey (Reference Passey2014). Overlying the Beinisvørð Formation basaltic lava flows and a post-Prestfjall Formation sill are <4 m of volcaniclastic sandstones and conglomerates which are, in turn, overlain by >4 m of clay ironstone. The uppermost beds of the Prestfjall Formation comprise 5.5 m of nodular mudstone, bright coal, and organic-rich mudstone. A sharp and planar contact separates the Prestfjall Formation from the overlying 3 m thick sequence of devitrified tuff beds attributed to the Hvannhagi Formation (Passey Reference Passey2004, Reference Passey2014).
Palynological analysis of a suite of samples from the Prestfjall Formation at this exposure yielded palynofloras completely dominated by septate fungal hyphae, with rare occurrences of gymnosperm pollen (Fig. 3). An abrupt change in the character of the assemblages was recorded above 10 m from the base of the section, within the clay ironstone bed (Fig. 3). Samples above this level are largely barren of palynomorphs, which given the nature of the organic rich mudstone is noteworthy, possibly reflecting oxidation of organic matter in situ. A single sample at the top of the organic-rich mudstone (Bed 6, Fig. 3) yielded a palynoflora including swamp cypress pollen (Inaperturopollenites hiatus), the juglandaceous angiosperm pollen Caryapollenites circulus, and the early successional pteridophyte spore Laevigatosporites haardtii. This assemblage is indicative of early–mid succession vegetation growing within the catchment, but transported into the lacustrine depositional environment (Jolley et al. Reference Jolley, Bell, Williamson and Prince2009).
Variations in clay ironstone facies were used by Passey (Reference Passey2014) to demonstrate that the Prestfjall Formation at Holið í Helli was deposited in a shallow lacustrine environment. The absence of rootlets in the succession indicates that the shallow lacustrine facies did not become emergent at any time, although Passey, (Reference Passey2014) suggested that oxidation rims in siderite spherules indicated seasonal water depth variations. Deposition in a shallow lava field lake is supported by the palynological data, accumulations of fungal debris being characteristic of lacustrine environments in Cenozoic lava fields (Jolley et al. Reference Jolley, Bell, Williamson and Prince2009). The low diversity, early–mid succession flora of Bed 6 (Fig. 3) is interpreted as having been deposited in a proximal, shallower lake. Subsequent pyroclastic activity forming the oldest units of the Hvannhagi Formation resulted in interaction of tephra with the Prestfjall Formation-lake, leading to devitrification and the layering of the oldest of the three tuffaceous units. From the evidence available at the Holið í Helli location, it is not possible to determine unequivocally if the onset of pyroclastic activity caused a cessation to lacustrine carbonaceous sediment deposition. However, the presence of common allochthonous petrified wood fragments in the oldest tuffaceous bed (Bed 7, Fig. 3) indicates that deposition of the Prestfjall Formation lacustrine sediments was potentially coeval with the onset of Hvannhagi Formation pyroclastic activity. Reworking of these wood fragments into the Hvannhagi Formation from older Prestfjall Formation strata is thought unlikely, as there is no accompanying coal or volcaniclastic siltstone lithoclasts recorded in Bed 7.
2.3.2. Rokhagi mine section
This section was recorded and palynology samples collected in an active mine adit in 1997, where the working exposed a section similar to that recorded by Lund (Reference Lund, Bott, Saxov, Talwani and Thiede1983, Reference Lund1989) and Parra et al. (Reference Parra, Delmont, Dumon, Ferragne and Pons1987), although the basal part of the ‘underclay’ and the ‘roof clay’ were not exposed (Fig. 4). Palynofloras from samples taken from the ‘underclay’ were dominated by the swamp cypress pollen Inaperturopollenites hiatus, with few Pityosporites (Pinaceae). Inaperturopollenites hiatus is abundant in many interbed successions, the parent Cupressaceae, swamp cypress plants being pollen overproducers, resulting in the taxon occurring abundantly in a range of flow top and floodplain depositional environments in lava fields (Jolley Reference Jolley and Widdowson1997; Jolley et al. Reference Jolley, Bell, Williamson and Prince2009). The rest of the ‘underclay’ assemblage is diverse with low taxon dominance. Assemblages include floodplain fern-derived pteridophyte spores (Deltoidospora), juglandacous angiosperm pollen (Momipites species and Caryapollenites species), and Alnipollenites verus (Betulaceae and Alnus types). These were derived from a mid-seral succession transitional swamp plant community (Jolley et al. Reference Jolley, Bell, Williamson and Prince2009, Reference Jolley, Gilmour, Gilmour, Kemp and Kelley2015), indicating longer-term substrate stability within the catchment area. Occurrences of the chlorophycean algae Botryococcus braunii, indicate deposition in a freshwater lacustrine environment (Tappan Reference Tappan1980).
The overlying dull coal unit lacks a leached seat earth and has no preserved root structures. It yielded a diverse palynoflora of similar composition to the ‘underclay’. The common wetland angiosperm pollen, Alnipollenites verus declines in abundance up section, with I. hiatus becoming dominant. Dominance of this swamp cypress pollen overproducer is characteristic of distal lacustrine facies (Jolley et al. Reference Jolley, Bell, Williamson and Prince2009; Daly et al. Reference Daly, Jolley, Spicer and Ahlberg2011a). This lower, fusain rich coal is separated from the upper vitrain rich bright coal by a thin grey shale unit (‘mid clay’), again with no root traces present. The lower ‘mid clay’ is dominated by Laevigatosporites haardtii, spores derived from a primary colonist pteridophyte fern. These ferns were common at riparian fluvial margins and disturbed forested areas (Collinson Reference Collinson2002; Jolley et al. Reference Jolley, Bell, Williamson and Prince2009; Daly et al. Reference Daly, Jolley, Spicer and Ahlberg2011a). In the upper part of the ‘mid-clay’ these fern spores are replaced by a diverse angiosperm pollen flora, which includes Momipites spp., Caryapollenites circulus, and C. veripites (Fig. 5). These taxa were derived from juglandaceous angiosperms, characteristic of early mid-successional vegetation swamp communities (Wing & Hickey Reference Wing and Hickey1984; Jolley et al. Reference Jolley, Bell, Williamson and Prince2009; Daly et al. Reference Daly, Jolley, Spicer and Ahlberg2011b). Samples from the ‘upper coal’ are totally dominated by swamp cypress pollen (I. hiatus) in low diversity, allochthonous assemblages.
Occurrences of the acritarch Chomotriletes and the chlorophycean algae Schizosporis parvus are indicative of waterlogged soils and lacustrine biofacies (Tappan Reference Tappan1980). These algae, the lack of root fossils, the absence of soil profiles, and the dominance of swamp cypress pollen in low diversity assemblages support interpretation as a lacustrine depositional environment (Lund Reference Lund1989; Passey Reference Passey2014). This interpretation is further supported by the finely laminated nature of the Prestfjall Formation units at exposure (Fig. 4a). Within the lacustrine basin, the durain and fusain rich coals were formed from concentrations of drifted plant material, although dispersed woody plant material occurred in the majority of lithological units.
Although >2 m of the ‘roof clay’ was sampled in the Rokhagi mine exposed faces, the palynofloras show a similar, apparently seral successional trend to that recorded in the underlying ‘mid clay’ to ‘upper coal’ interval. Primary successional fern dominated (Laevigatosporites haardtii) assemblages were replaced by diverse, angiosperm pollen (dominantly Momipites and Caryapollenites species) rich palynofloras (Figs 4, 5). These were in turn succeeded by palynofloras wholly dominated by swamp cypress pollen (I. hiatus). These repeated apparent seral succession trends reflect increasing distance from the input point to the lacustrine system and potentially, increased lake water depth.
2.3.3. KOLBH1-04 borehole
The KOLBH1-04 borehole [61.56888N, 6.93527W] was drilled in 2004 on the west side of the Rokhagi valley to test coal reserves away from the then open mines on the eastern side of the valley (Fig. 1). A number of thin coaly intervals occur in the lower ~7 m of the Prestfjall Formation, which overlie an aphyric basalt lava flow of the Beinisvørð Formation. There is no differentiated ‘underclay’ although the named coal units of Rasmussen & Noe-Nygaard (Reference Rasmussen and Noe-Nygaard1970) were identified at the time of drilling (Fig. 6). The ‘roof clay’ unit of the Prestfjall Formation is overlain by the Hvannhagi Formation consisting of 0.6 m of volcaniclastic sandstone succeeded by 0.3 m of peperite beneath the overlying olivine-phyric flows of the Malinstindur Formation lava flow field. The red-brown mudstones of the ‘takleir’ unit of the roof clay are not present in this borehole although it is <0.5 km from the Rokhagi mine adits.
Palynofloras recovered from this succession are similar to those recorded in the Rokhagi mine section (Fig. 4). All were dominated by I. hiatus (Cupressaceae, Taxodium) and Pityosporites spp. (Pinaceae, pines) associated with common occurrences of Caryapollenites circulus (Juglandaceae, Carya) and Laevigatosporites haardtii (a primary colonist polypodiaceous fern). These taxa are associated with frequently occurring Caryapollenites species (C. veripites, C. inelegans and C. triangulus, see Figs 5, 6), and fagaceous pollen (Cupuliferoipollenites cinglum and Cupuliferoidaepollenites liblarensis). In the KOLBHI-04 borehole, these more diverse assemblages were recovered from four separate intervals within the Prestfjall Formation succession (Fig. 6).
Samples taken at 217.4 m, 218.7 m, and 219.0 m all yielded marine dinoflagellate cysts including Areoligera cf. coronata and A. cf. senonensis (Fig. 6) associated with brackish to freshwater acritarchs and chlorophycean algae (Pediastrum bifidites, Botryococcus braunii, Schizosporis parvus and Leiosphaeridia spp.). The dinoflagellate cysts in this assemblage are often associated with shallow marine sediments, particularly shoreface deposits (Pross & Brinkhuis Reference Pross and Brinkhuis2005; Jolley et al. Reference Jolley, Millett, Schofield, Broadley and Hole2021). Because they behave as silt-size sedimentary particles, they are prone to transportation by tidal action within tidally influenced fluvial channel depositional systems to the extent of the tidal limit (Anderson Reference Anderson1998; Debenay et al. Reference Debenay, Carbonel, Morzadec-Kerfourn, Cazaubon, Denèfle and Lézine2003).
At the upper boundary of this microplankton bearing interval are occurrences of spores and pollen reworked from older sedimentary rocks (Fig. 6). These taxa indicate reworking of Carboniferous, Triassic, and Late Cretaceous palynofloras and transportation into the FIBG lava field, probably within lithoclasts.
FAAS analysis was conducted on the same sample set used for palynological analysis. Attention was focussed on relative trends within the FAAS data, particularly ratios which highlight volcaniclastic input and weathering (Retallack Reference Retallack1999; Sheldon, Reference Sheldon2003; Jolley et al. Reference Jolley, Widdowson and Self2008; Reference Jolley, Passey, Hole and Millett2012). Both the barium:strontium (Ba:Sr), and aluminium:phosphate ratios measured mobile element loss through weathering and plant uptake (Vitousek Reference Vitousek2004). These show repeated depletion trends (Fig. 6), indicative of fluctuation in substrate stability and biomass in the catchment. Mobile element depleted samples correspond to more diverse, mid-seral succession palynofloras (Vitousek Reference Vitousek2004; Jolley et al. Reference Jolley, Widdowson and Self2008). These were derived from moderately stable substrates within the catchment at periods of lower lake water levels. Low weathering ratios are positively correlated with low diversity, high dominance palynofloras, and are characteristic of deeper lake states (Fig. 6).
Using an approach similar to that of Vitousek (Reference Vitousek2004), Sheldon (Reference Sheldon2003) and Sheldon & Tabor (Reference Sheldon and Tabor2009), residual magnesium oxide (MgO) concentrations were calculated for the sedimentary rocks sampled. In the uppermost part of the KOLBH1-04 cored succession, concentrations of MgO rise in relation to the average value for Beinisvørð Formation basalts (Fig. 6). This rise is partly due to the occurrence of a volcaniclastic sandstone at 219.4 m, but increasing values from 217.0 m, continuing up section imply increasing volcanic input in the uppermost Prestfjall Formation.
3. Thickness variation of the Prestfjall and Hvannhagi formations on Suðuroy
3.1. Prestfjall Formation
Observations on the thickness and facies of the Prestfjall Formation on Suðuroy are largely dependent on limited exposures (e.g., Fig. 2) and the measurements of Rasmussen & Noe-Nygaard (Reference Rasmussen and Noe-Nygaard1970; Supplementary Data S3). Many of their measured sections (see Fig. 7a) were established in now long-closed coal mine adits and coal prospecting excavations. Additional data from the KOLBH1-04 borehole, field observations, and from the photographic profile of the west Suðuroy cliff line are included here (Passey et al. Reference Passey, Vosgerau, Sørensen and Vautravers2016). With these additional data, the creation of isopach maps and integration of palynological data, a modified depositional environment model can be proposed.
An isopach map for the gross thickness of the Prestfjall Formation on Suðuroy (Fig. 7a), shows the thickest deposits to be centred on the west side of the Rokhagi valley, in the north and in the extreme east at Holið í Helli [61.547N, 6.746W]. Overall, the Prestfjall Formation thins abruptly across the centre and northeast of the island centred around Stapin ([61.583N, 6.829W] Figs 1, 7a). There is no evidence for the Prestfjall Formation having been deposited in situ in the Stapin area, although allochthonous rafts of Prestfjall Formation occur in thick sills intruding the Hvannhagi Formation (Fig. 4a). At Stapin, volcaniclastic Hvannhagi Formation conglomerates and shales overlie an older, Beinisvørð Formation volcanic vent (the ‘Stapin Vent’; Passey Reference Passey2004; Passey & Jolley Reference Passey and Jolley2009). Northwest of Stapin, thick Prestfjall Formation coals and shale beds are inferred from peperitised rocks in the base of the Hvannhagi Formation succession at Reyðibarmur ([61.610N, 6.924W] Figs 1, 7a). These exposures are dissimilar to others of the Prestfjall Formation in being overlain by a thick volcaniclastic Hvannhagi Formation succession (Fig. 7h).
Thin profiles of the Prestfjall Formation have been recorded over central Suðuroy, extending from Tjørnunes [61.598N, 6.890W] on the northeast coast to Fámjin [61.525N, 6.881W] on the central west coast. Although exposures are limited, a profile comprising basaltic lava flows and thin volcaniclastic sediments at Hvannagjógv [61.569N, 6.889W] indicates that a localised flow terrain may have influenced part of the western limit of the Prestfjall Formation lacustrine basin. Difficulties in differentiating intrusive from extrusive basalt units also appear to have resulted in profile 20 near Tjørnunes being underestimated or misinterpreted (Fig. 7a; Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970). These central Suðuroy extrusive volcanics and the intrusion of sills into the Prestfjall Formation are orientated from Reyðibarmur in the north, south southeast to profiles 29, 30 and 31 on Oyrnafjall ([61.547N, 6.849W] Fig. 7a, d, e, f; Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970). Here, basalt intrusions invade the Prestfjall Formation in a manner similar to that at Reyðibarmur. Deposition within the southern area of the lacustrine basin was thickest at Rangabotnur [61.551N, 6.890W] and to the west (Fig. 7a), while east of profile 36, pre-existing topography resulted in a thin sedimentary succession.
Isopach mapping of the gross thickness of coaly sediments of the Prestfjall Formation was undertaken by Ellis et al. (Reference Ellis, Jolley, Passey and Bell2009). In addition to mapping the gross thickness, isopach mapping of the five informal units of the Prestfjall Formation (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970) has highlighted incompatibilities between individual profiles, particularly with reference to the KOLBH1-04 borehole. Isopach mapping of the ‘underclay’ (Fig. 7b) shows that the borehole contains an anomalously thick unit in comparison to those only a few tens of metres away (e.g., Figs 4, 6). In contrast, the ‘lower coal’ (Fig. 6) is significantly thinner in the KOLBH1-04 borehole and in profiles 10–13 (Fig. 7c) immediately to the south, than in the Rokhagi profile (Fig. 4). The regular occurrence of coal lenses and flasers in the ‘underclay’ of KOLBH1-04 and profiles 10–13 indicate lateral facies change from these lacustrine proximal shales to lacustrine distal, durain-rich coals. Comparison of palynological data from the Rokhagi mine profile with that in the KOLBH1-04 borehole (Figs 4, 6) supports this interpretation of the lithological data. Based on palynological data, both the ‘upper’ and ‘lower’ coals in the KOLBH1-04 profile are age equivalent to the ‘upper coal’ in the Rokhagi mine profile. Similarly, the ‘lower’ coal of the Rokhagi mine profile is correlative with the ‘underclay’ of the KOLBH1-04 borehole.
Given that the carbonaceous content of the Prestfjall Formation is allochthonous, lateral variation in the concentration of drifted organic debris would have resulted from water current and wind effect on the lake basin. Prevailing wind directions from the northwest in the Selandian to lower Ypresian are indicated from tephra distribution in the North Sea and Faroe–Shetland basins (Knox & Morton Reference Knox, Morton, Morton and Parson1988). Amplified by fluvial input down the regional slope from the north northwest, this would have resulted in accumulations of woody debris along the distal, eastern margin of the lacustrine basin, against the Reyðibarmur–Hvannagjógv–Oyrnafjall zone (Fig. 7c). Although the highly carbonaceous profiles recorded in mines on the east side of the Rokhagi valley are not the thickest of the Prestfjall Formation successions (Fig. 7a, c), the thicker coal units in these locations would have been subjected to greater compaction (compaction factors of 6–10; Elliott, Reference Elliott1984) than areas of dominantly clastic sedimentary rocks (compaction factors of 1.2–2.0; Mondol et al. Reference Mondol, Bjørlykke, Jahren and Høeg2007; Bjørlykke et al. Reference Bjørlykke, Jahren, Mondol, Marcussen, Croize, Peltonen and Thyberg2009). Decompaction of these profiles suggests that the eastern margin of the Prestfjall Formation basin was at its deepest on the eastern side of Rokhagi.
Conglomerates and normally graded pebbly sandstones recorded at the top of the Prestfjall Formation in the northwest of the basin (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970) mark a period of uplift and erosion. These deposits have also been recorded on the west coast of Suðuroy near Grímsfjall ([61.598N, 6.987W] Fig. 2d). Here, a paleovalley (15–20 m in depth) eroded into the Prestfjall Formation was filled by pillow-based, columnar jointed lava flows of the Malinstindur Formation. The erosive, conglomerate and sandstone filled fluvial channels present in the northwest of the lacustrine basin are not observed to cut the red-brown claystones of the ‘takleir’ unit of the ‘roof clay’. Combination of conglomerate–sandstone filled channel thicknesses with those of the ‘taklier’ claystone (Fig. 7g) indicate that the claystones were overbank facies of the fluvial channel fill sandstones and conglomerates. This combined ‘taklier’ and conglomerate unit is confined to the northwest of the Prestfjall Formation lacustrine basin (Fig. 7g), coincident with the lava filled valley eroded into the top of the Prestfjall Formation at Grímsfjall (Fig. 2d). Together, these indicate incision of a limited fluvial drainage system into the carbonaceous Prestfjall Formation. Uplift and incision of this area could have been a response to localised thermal uplift of the area following the onset of Malinstindur Formation volcanism. However, this erosion occurred in the same period as uplift and incision of the eastern margin of the Faroe–Shetland Basin following relative sea level fall. Here, basin flank uplift resulted in the rejuvenation of shelf margin drainage systems to form the age equivalent Flett Unconformity (Shaw Champion et al. Reference Shaw Champion, White, Jones and Lovell2008; Hartley et al. Reference Hartley, Roberts, White and Richardson2011; Jolley et al. Reference Jolley, Millett, Schofield, Broadley and Hole2021).
3.2. Hvannhagi Formation
The gross thickness isopach for the Hvannhagi Formation (Fig. 7h) clearly identifies a depocentre in the northeast of Suðuroy. Composed of >25 m thickness of tuffs and volcaniclastic sedimentary rocks (Passey & Jolley Reference Passey and Jolley2009), these are at the thickest between Stapin and Hvannhagi proximal to the source. These rocks thin significantly to the south with distance from source, where ~3 m of tuffaceous rocks are recorded at Holið í Helli (Fig. 3). Mapping the western limit of these rocks is problematic, their exposure being obscured by thick dolerite sills. However, they appear to pass laterally into volcaniclastic sedimentary rocks and basaltic lava flows at Hvannagjógv. Overlying the intruded, fluidal peperitised Prestfjall Formation strata exposed at Reyðibarmur, the Hvannhagi Formation includes normally graded volcaniclastic sandstones, lapillistones, and tuffs. These are composed of palagonitised volcanic glass, larger fragments exhibiting cuspate edges. The absence of intrusions in these sedimentary rocks suggests that the peperitisation of the underlying Prestfjall Formation took place during shallow burial, prior to Hvannhagi Formation deposition.
Approximately 2 m above the base of the Hvannhagi Formation at Leirbakki ([61.581N, 6.831W] Fig. 1), clast supported conglomerates made up of 2 cm–8 cm pebbles occur in lensoid bodies up to 10 m across (Fig. 2a, c). Resting unconformably on the surface of the late Beinisvørð Formation Stapin Vent, the conglomerates include large blocks up to 50 cm wide, but are overlain by up to 1 m of finely bedded mudstones which include some thin, bright coals up to 2 cm thick (Fig. 2a). These shales and drifted coals yielded a lacustrine palynoflora dominated by taxa common in transported assemblages (Inaperturopollenites hiatus and Pityosporites spp.). Occurrences of Caryapollenites circulus and Cupuliferoipollenites cingulum subsp. oviformis in these shales indicate coeval deposition of lowermost Hvannhagi Formation volcaniclastic and upper Prestfjall Formation carbonaceous beds.
A second, different flora was recorded from the upper part of the Hvannhagi Formation succession in the Svalbaá stream [61.571N, 6.835W], near Trongisvágur (Fig. 1). Here, C. circulus occurs commonly with Caryapollenites inelegans and Caryapollenites veripites in a coal rip-up clast within the upper conglomerate beds. This clast along with others in the uppermost Hvannhagi Formation are clearly derived from the Prestfjall Formation coals, presumably from an area in the northwest where they have been largely eroded away.
4. Pre-rift to syn-rift rock succession on Vágar and adjacent islands
The northern exposure of the Prestfjall and Hvannhagi formations crops out on the islands of Tindhólmur and Vágar (Fig. 1). Exposure is largely restricted to coastal cliffs, which are mostly inaccessible, other than when viewed from out to sea. Because of this inaccessibility combined with poor exposure, there are no detailed published records of the Beinisvørð–Malinstindur transition in this area. A borehole (Gásadalur 1989-1) drilled near Bøur on Vágar penetrated this succession and this provides the most comprehensive record of Prestfjall and Hvannhagi formations sedimentary facies on Vágar. The small island of Tindhólmur presents the least disrupted exposure of the Beinisvørð–Malinstindur transition and is considered first.
4.1. Tindhólmur
The south facing cliffs of the island of Tindhólmur expose a section through the Prestfjall and Hvannhagi formations and the base of the Malinstindur Formation (Fig. 8). On Mykineshólmur (Fig. 1) and the on western end of Tindhólmur, thick Beinisvørð Formation simple lava flows dip to the southeast. On Tindhólmur, these are overlain by the Prestfjall Formation which is composed of horizontally bedded volcaniclastic coarse-grained to fine-grained sandstones in beds 10 cm–30 cm thick. These sandstones are ~5 m thick in the west, thickening abruptly over a Beinisvørð Formation lava flow termination to ~28 m over the centre and east of the island. There is no evidence of a significant erosive unconformity at the base of the Prestfjall Formation, the sedimentary rocks having been deposited at the eastern side of a series of Beinisvørð Formation simple lava flow terminations (Fig. 8).
Numerous shallow invasive flows and sills are exposed at the boundary between the Prestfjall Formation and the overlying Hvannhagi Formation, obscuring the nature of the contact (Fig. 8b). In the centre of the island, a 260 m wide channel eroded into the underlying Prestfjall Formation was filled with crudely bedded lava flows, which overstep the edge of the channel structure onto the upper surface of the Prestfjall Formation (Figs 8, 9a). Orientated 160o south southest, this channel occurs in the middle of the south face of Tindhólmur, its eastern margin being influenced by a north northeast–south southwest trending fault (Figs 8, 9, 10), interpreted as dominantly strike slip with some minor downthrow component to the west. The upper surface of these valley fill compound lava flows is reddened to a few centimetres (Fig. 9a), resulting from an eruptive hiatus and weathering of the lava flow surface.
A second steep sided valley fill is partly exposed at the eastern tip of the island (Figs 8b and 9b). The valley is >150 m wide, the eastern margin having been eroded away. The valley was incised through Prestfjall Formation sandstones, and was subsequently filled by compound lava flows of the Malinstindur Formation (Figs 8b, d, 9b). Overlying the valley fill of compound flows is a 4.8 m thick succession of mostly clast supported, coarse polymictic basalt clast conglomerates of the Hvannhagi Formation (Fig. 8b). These thicken into the centre of Tindhólmur, reaching 24 m–26 m thickness on either side of the central valley, thinning to 7.2 m in the west (Fig. 8b). In this unit, low angle, trough cross bedding is evident across the island, with evidence of incision and erosion of earlier Hvannhagi Formation sediments. A well-exposed section (Fig. 9c) in the west of the island shows that the upper Hvannhagi Formation transitioned upsection from the clast supported conglomerates to clast poor conglomerates and volcaniclastic sandstones, some with cross bedding.
The Hvannhagi Formation is interpreted to have been deposited as a volcanic lahar (Smith & Lowe Reference Smith, Lowe, Fisher and Smith1991). Clast rich, ungraded beds were deposited under hyperconcentrated flow. Clast poor, weakly horizontally bedded and cross bedded sandstones in the upper beds (Fig. 9c) reflect the transition from deposition by hyperconcentrated flow to streamflow. Shallow incisions of younger, clast rich conglomerate beds into older Hvannhagi Formation conglomerates are associated with wide channels and clinoforms, suggesting that these were deposited by a dilute lahar flow (cf. Pierson et al. Reference Pierson, Pringle and Cameron2011).
A subsequent period of erosion removed some of the previous valley fill flows and created accommodation space along the axis of the central valley (Figs 9a, 10). Rejuvenation of the Tindhólmur central valley and differential uplift of the east flank by 8 m with respect to the west flank, was calibrated from the photogrammetry model (Figs 8, 9a). Differential movement of one valley flank with respect to the other is attributed to further movement on the valley edge fault (Fig. 10). This irregular topography was onlapped and overstepped by the prograding Malinstindur Formation compound flow field (Figs 8, 9, 10). In addition to infilling the central valley structure, these flows covered the western flank of the central valley before the uplifted eastern flank was submerged (Fig. 10).
4.2. Vágar
The Gásadalur 1989-1 borehole (Fig. 1; [62.09471N, 7.40663W]) was drilled on the island of Vágar in 1989, approximately halfway between the villages of Gásadalur and Bøur, penetrating the Malinstindur, Hvannhagi and Prestfjall formations before terminating in the uppermost lava flows of the Beinisvørð Formation. The Prestfjall Formation rests on a basalt lava flow where a sandy paleosol developed (Fig. 11). Overlying this are approximately 11 m of volcaniclastic sandstones and conglomerates attributed to the Prestfjall Formation. These are, in turn, overlain by matrix supported, very poorly sorted, polymictic, granulite to pebble grade volcaniclastic conglomerates of the Hvannhagi Formation. Low-angle cross bedding occurs in the sandstone units of the Prestfjall Formation while the conglomerates of the Hvannhagi Formation appear as unstructured or normally graded beds.
Fifty-two samples from the borehole were analysed for palynology, although recovery was concentrated within the Prestfjall Formation. These palynofloras are similar in composition to those recorded in the Prestfjall Formation of the KOLBH1-04 borehole. Dominated by Pityosporites spp. and Inaperturopollenites hiatus (Figs 5, 11) these floras are typical of fluvially transported assemblages (Daly et al. Reference Daly, Jolley, Spicer and Ahlberg2011a). In comparison to Prestfjall Formation records from Suðuroy, Caryapollenites species, typical of transitional swamp floodplain communities, occur at lower frequency in the Gásadalur borehole. Palynofloras also decline in diversity and frequency in the upper beds of the Prestfjall Formation. Common Laevigatosporites haardtii and Deltoidospora adriennis fern spores were recorded in the upper beds of the Prestfjall Formation (Fig. 11). Laevigatosporites haardtii is derived from ferns which exploited disturbed substrates and newly exposed overbank areas (Collinson Reference Collinson2002; Jolley et al. Reference Jolley, Bell, Williamson and Prince2009), while the parent plant of D. adriennis was associated with colonisation of crevasse splays and inter channel sand and mud banks (Jolley et al. Reference Jolley, Bell, Williamson and Prince2009; Daly et al. Reference Daly, Jolley, Spicer and Ahlberg2011a). These beds also yielded common Triporopollenites coryloides (Myricaceae, Myrtle types) which indicate dryer, better-drained substrates associated with sandy floodplains or sandy channel margins.
In addition to these indications of a dynamic depositional environment, the dinocyst Impletosphaeridium ligospinosum, the prasinophycean algae Leiosphaeridia spp., and reworked Mesozoic palynomorphs occurred in the uppermost siltstones of the Prestfjall Formation (40.29 m–40.20 m; Fig. 11). Characteristic of low salinity marginal marine conditions, their occurrence indicates that the upper Prestfjall Formation at Gásadalur was deposited within tidally influenced fluvial channels. As with the uppermost bed of the Prestfjall Formation in the KOLBH1-04 borehole on Suðuroy, the presence of reworked Mesozoic palynomorphs (Fig. 11), indicates a drainage system integrated outside the lava field. The source area of these reworked palynomorphs is currently unknown, as with other occurrences of Mesozoic and some Carboniferous taxa in the upper Malinstindur Formation and Sneis Formation interbeds (Passey & Jolley Reference Passey and Jolley2009). Given the transition into a marine environment in the Faroe–Shetland Basin to the southeast (Jolley et al. Reference Jolley, Millett, Schofield, Broadley and Hole2021), a clastic source to the northwest is thought probable.
The pollen and spore record from the Prestfjall Formation in the Gásadalur borehole was derived from disturbed fluvial margin and flow top habitats, deposited in a fluvial environment. Because the overlying Hvannhagi Formation was deposited by a high-energy system derived from an active volcanic zone, the palynofloras are poor. Only low frequencies of fluvially transported taxa occur. However, they do indicate that the catchment area of the Hvannhagi Formation conglomerates was vegetated to some extent.
Samples used for palynological analysis were also subjected to XRF analysis of major and trace elements (Fig. 11, Supplementary data S2). The plot of MgO values against an average value for the Malinstindur Formation highlights the first evidence of volcanism in the uppermost few centimetres of the Prestfjall Formation prior to the comprehensive volcanism of the Malinstindur Formation. This is supported by the Ba:Sr and Al2O4:base ratios which drop significantly in this interval.
The architecture of the Prestfjall and Hvannhagi formations deposits around Vágar is obscured by the inaccessible nature of the exposures. Boat surveys have proven to be of limited value because of the low-angle of view, although Rasmussen & Noe-Nygaard (Reference Rasmussen and Noe-Nygaard1970) recorded complex disruption of the Prestfjall Formation from basaltic intrusions. Because of this lack of accessibility, a photographic survey was undertaken from a helicopter which provided the basis for a photogrammetry model (Passey et al. Reference Passey, Vosgerau, Sørensen and Vautravers2016). This model provided coverage of the Prestfjall and Hvannhagi formations over the islands of Tindhólmur and western Vágar.
The succession drilled by the Gásadalur 1989-1 borehole is exposed in inaccessible sea cliffs southeast and northwest of the village of Bøur on Vágar. Creation of a photogrammetry model from the stereo photographic survey of the western coastline of Vágar has allowed these sections to be placed in the wider context of the transition between the Beinisvørð and Malinstindur formations (Fig. 12).
Mapping of the Beinisvørð–Prestfjall Unconformity surface (Passey & Varming Reference Passey and Varming2010) identified an east southeast dip across Vágar. The intra Malinstindur Formation Kvívík Beds (Passey & Jolley Reference Passey and Jolley2009), also showed a dip to the east, indicating little overall structural change during the deposition of the Prestfjall and Hvannhagi formations. Examination of photogrammetry profiles of the west coast of Vágar showed that the interval between the Beinisvørð–Prestfjall Unconformity surface and the base of the Malinstindur Formation thins from Bøur [62.086N, 7.369W] in the south to Barðið [62.143N, 7.457W] in the north. From Barðið to Ritunøva [62.145N, 7.341W] on the north coast, the same interval thins to a few metres before it passes below sea level to the west of Viðvík ([62.124N, 7.356W]; Figs 1, 12).
Thinning of the Prestfjall–Hvannhagi formation interval to the north and east of Vágar is linked to the underlying topography at the cessation of Beinisvørð Formation eruptions. In north-western Vágar, the exposed upper Beinisvørð Formation is composed of a thick sequence of compound basalt flows onlapped by a succession of dipping simple lava flows and volcaniclastic rocks, which thicken to the south (Figs 12, 13). These compound lava flows are interpreted to have originated from a low-angle shield volcano point source (Passey & Bell Reference Passey and Bell2007), here termed the Mykines–Vágar low-angle shield volcano. An isopach map for the upper unit of the compound flow field demonstrates a source to the northwest (Fig. 13a). This is honoured by the onlap pattern of the subsequent Beinisvørð Formation simple lava flows, which overlie the compound flow field (Figs 12, 13b). A gravity anomaly, some 15 km north of Mykines was identified by Schrøder (Reference Schrøder1971); this probably reflects the location of the Mykines–Vágar low-angle shield.
In the cliffs north and south of Gásadalur village, the simple lava flows of the upper Beinisvørð Formation are interbedded with dark red volcaniclastic sandstones and siltstones which thicken to the south (Fig. 12b). Where exposed at Reyðastíggjatangi [62.104N, 7.436W], these are composed of blocky, fissile, red-brown volcaniclastic shales and poorly sorted, fine sandstones ~1.3 m thick. Low angle, trough bedding occurs in exposures of this unit tens of metres to the north. Two of these shale units are exposed, separated by a simple lava flow, the oldest shale resting on south southeast prograding hyaloclastite foresets (Fig. 12a). These hyaloclastites are correlative with a simple lava flow exposed to the north (Fig. 12b). Uppermost Beinisvørð Formation red volcaniclastic sedimentary interbeds persist and thicken to the south southeast towards Bøur before dipping below sea level (Figs 12, 14). Together, the change in volcanic facies and south southeast thickening sandstones reflect an increase in accommodation space south of Gásadalur village. This exploited topography is inherited from the southern margin of the Beinisvørð Formation, Mykines–Vágar low-angle shield volcano. This accommodation space was subsequently filled by the deposition of the Prestfjall and Hvannhagi formations.
Exposure of the Prestfjall Formation on Vágar is limited, cropping out in inaccessible sea cliffs north and south of Gásadalur. A record of 90 cm of shales with coaly material included was made by Rasmussen & Noe-Nygaard (Reference Rasmussen and Noe-Nygaard1970) at Hellisgjógv ([62.144N, 7.449W] Fig. 1) in the northwest. No beds of carbonaceous sediment or disseminated woody debris were observed during fieldwork for this study.
Along the coastal cliffs from Gásadalur to Ritunøva, the Prestfjall and Hvannhagi formations are very poorly exposed and were regarded as one unit in the photogrammetry interpretation of the Gásadalur to Bøur section (Fig. 12). From Gásadalur to Bøur, there are isolated clear exposures of the Prestfjall Formation, but as noted by Rasmussen & Noe-Nygaard (Reference Rasmussen and Noe-Nygaard1970), the volcaniclastic sedimentary rocks are extensively intruded by sills and invasive lava flows. Near Reyðastíggjatangi (Fig. 14a), a ~18 m thick unit of the Prestfjall Formation crops out. Overlying a basal sandstone bed <2 m thick, red to orange, crudely bedded, poorly sorted, and fine-grained sandstones and mudstones occur. These sedimentary units have a sharp upper contact with the overlying clast-rich conglomerates of the Hvannhagi Formation.
The Hvannhagi Formation in the Gásadalur to Bøur section is thicker and more complex than any recorded previously (Figs 12, 14). Dominantly composed of clast supported, ungraded and normally graded, crudely trough bedded volcaniclastic conglomerates, and thin planar bedded sandstones, these beds are ~4 m thick (Figs 12, 14). Clasts are angular to sub-rounded and composed of basalt derived from flow top and inter flow facies from different lava flows. Examination in the field shows that the matrix of these beds is composed of altered volcanic glass and sand grade volcaniclastic material.
Two depositional phases of the Hvannhagi Formation have been interpreted from the photogrammetry model, but exhibit no obvious lithological difference. The oldest unit occurs along the whole of the Gásadalur to Bøur cliff exposure, except where the latter has been penecontemporaneously eroded (Fig 12, Gásadalur to Bøur section, 14b). In the north, this unit attains 53 m thickness, expanding to ~80 m in the south of the exposure. Overlying this unit is a second, younger conglomerate unit. Separated from the older Hvannhagi Formation conglomerates by shallow invasive flows and sills which follow the boundary, the upper volcaniclastic conglomerates reach a maximum thickness of 94 m (Fig. 14b). The edge of this upper Hvannhagi Formation unit appears to have been penetrated by the Gásadalur 1989-1 borehole. Separated from underlying conglomerates by an intrusive sheet (Fig. 11, 75.6–81.28 m), the upper conglomerate unit (81.28 m–92.62 m) is ~11 m in thickness. However, there is a significant angular unconformity at the boundary between the Hvannhagi Formation and the overlying Malinstindur Formation, with a topographical relief of up to 119 m over 1.2 km of outcrop (Fig. 12 Gásadalur to Bøur section). Due to this unconformity, the depositional architecture of the upper Hvannhagi Formation is difficult to reconstruct.
5. Synthesis
Examination of the Beinisvørð–Malinstindur transition exposed on the Faroe Islands has identified a trend in proximal to distal depositional environments from north to south. While all areas were impacted by volcanic activity, those in the Vágar to Tindhólmur area were dominated by lava field processes. Sections exposed on Suðuroy demonstrate dominantly lower energy depositional environments. Accordingly, integration of the depositional environments across these two areas is discussed down the environmental gradient from north to south.
5.1. Mykines–Vágar area
Photogrammetry and lava facies evidence of a large compound flow field in northwest Vágar and north Mykines, have confirmed the existence of the Mykines–Vágar low-angle shield volcano. Thick successions of compound lava flows seen on the easternmost tip of Mykines (Fig. 13), indicate that the structure was to the north of Mykines and northwest of Vágar, corresponding to an earlier gravity anomaly interpretation (Schrøder Reference Schrøder1971). Subsequently, upper Beinisvørð Formation simple lava flows and volcaniclastic sedimentary rock units onlapped the low-angle shield volcano. These crudely stratified sandstones and shales are lacking in any palynoflora, supporting an origin as weathered and remobilised volcanic ash derived from fissures that erupted the late-stage simple lava flows. Ash, scoria, and weathered basalt clasts from the fissure zone would have been remobilised by rainfall and transported by hyperconcentrated flow as a muddy slurry (cf. Smith & Lowe Reference Smith, Lowe, Fisher and Smith1991; Zernack Reference Zernack, Roverato, Dufresne and Procter2021). These muddy lahars progressively infilled a pre-existing low, previously a zone of lava–water interaction south of Gásadalur village. The occurrence and thickness of these volcaniclastic deposits is clearly not related to the duration between eruptions (Millett et al. Reference Millett, Hole, Jolley, Passey and Rossetti2020). Their stratigraphical and paleogeographical location is instead a depositional system response to the presence of a substantial low-angle shield volcano to the northwest.
Despite the topographical annealing effect of the late-stage Beinisvørð Formation simple lava flows and interbedded lahar deposits, the geometry of the succeeding Prestfjall and Hvannhagi formations indicates persistence of a paleovalley in the Gásadalur to Bøur area. In the area of Tindhólmur, the Prestfjall Formation depositional system occupied a topographical low to the east of successive Beinisvørð Formation simple lava flow terminations. The finely bedded volcaniclastic sandstones of the Prestfjall Formation exposed on Tindhólmur (Fig. 9b) indicate deposition of volcaniclastic sediment in a localised lacustrine basin. This contrasts with the Prestfjall Formation fluvial sandstone/shale deposits drilled in the Gásadalur 1989-1 borehole and exposed within the Gásadalur to Bøur paleovalley, indicating a source of remobilised volcanic ash from a primary volcanic accumulation.
No lithological change is apparent in the Prestfjall Formation on Tindhólmur, but the blocky claystones at the top of the formation in the Gásadalur 1989-1 borehole contain evidence of a tidally influenced fluvial channel and a link to a source of extra-basinal clastic reworking. Several rises in relative sea level have been recorded within the FIBG lava pile (Jolley et al. Reference Jolley, Passey, Hole and Millett2012), potentially linked to drowning of the low-lying basalt plain in response to thermal subsidence at the cessation of Beinisvørð Formation eruption (Millett et al. Reference Millett, Hole, Jolley, Passey and Rossetti2020). This would have resulted in the formation of ria-like drowned valleys rather than sheet flooding of the basalt plain (Shaw Champion et al. Reference Shaw Champion, White, Jones and Lovell2008; Jolley et al. Reference Jolley, Passey, Hole and Millett2012, Reference Jolley, Millett, Schofield, Broadley and Hole2021).
The 150 m thick conglomeratic succession of the Hvannhagi Formation on Vágar is interpreted as a sequence of lahar deposits derived from remobilisation of a sector collapse volcanic debris avalanche (VDA) in the northwest of Vágar. Evidence of a VDA source for the Hvannhagi Formation on Vágar is provided by the high volume of very poorly sorted basalt clasts, derived from multiple lava flows. These VDA deposits were reworked as dilute lahars (cf. Zernack Reference Zernack, Roverato, Dufresne and Procter2021) down the pre-existing valley system (Figs 12, 14–16). Repeated beds of very poorly sorted, unstructured to faintly bedded, and normally graded conglomerates indicate deposition under hyperconcentrated flow. Horizontally bedded sandstone units deposited in transition to fluvial flow occur in some units, suggesting dilution of the hyperconcentrated flow front, waning lahar states, or as lateral overbank facies. Further down the drainage slope at Tindhólmur, lahars initially deposited debris flow conglomerates. Here, waning flow resulted in the deposition of cross-bedded sandstones and pebbly sandstones showing transition from hyperconcentrated to streamflow (Fig. 9).
Considered by previous authors to be a transitional contact between the volcaniclastic deposits and compound lava flow field (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970; Waagstein Reference Waagstein, Morton and Parson1988; Passey & Jolley Reference Passey and Jolley2009), new photogrammetry data of the Hvannhagi–Malinstindur Unconformity indicate a more nuanced relationship. This unconformity marked a significant relative uplift event associated with dominantly strike-slip faulting both on Vágar and Tindhólmur. Faulting and erosion occurred along the fault scarp associated with erosion of a significant proportion of the upper Hvannhagi Formation, leaving an irregularly eroded topography of 119 m on Vágar. On Tindhólmur this relationship is seen to be more complex. Erosion of valley features and their subsequent partial infilling by basalt lava flows, was followed by a hiatus reflected in weathering of the central valley fill lavas. Subsequent deposition of the Hvannhagi Formation conglomerates was, in turn, followed by a rejuvenation of the central valley, with Malinstindur Formation compound lavas filling the central valley and onlapping the surrounding topography (Figs 8, 10). This relationship confirms that eruption of lava flows characteristic of the Malinstindur Formation began following Prestfjall Formation deposition, prior to the VDA and associated lahar deposits of the Hvannhagi Formation in the north of the islands. The volcaniclastic deposits of the Hvannhagi Formation therefore represent dynamic, landscape scale events within the early Malinstindur Formation lava field.
Mapping of the basal Malinstindur Formation lava flow package (Fig. 15), illustrates that the flow field originated to the northwest, indicating continued eruptive activity sourced from the Mykines–Vágar shield. Subsequent basalt plains eruptions drowned the complex erosional topography of the Hvannhagi–Malinstindur Unconformity. The eroded surface of the upper Hvannhagi Formation created a series of ‘mini-basins’ which were the last to be overstepped by the compound flow field.
5.2. Suðuroy
A thick succession of eastwards dipping Beinisvørð Formation simple lava flows is exposed on the west coast of Suðuroy. These are regarded as having been sourced from an eruptive fissure system west of Suðuroy (Waagstein Reference Waagstein, Morton and Parson1988; Passey & Varming Reference Passey and Varming2010). Consequently, late Beinisvørð Formation lava flows and sedimentary rock interbeds vary in composition and thickness west to east across the island. Typical of this are the ponded columnar jointed (e.g., Kulagjógv Flow) and multi-tiered (e.g., Hov Flow) lava flows exposed in the interior of the island and north of Hov, respectively (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970). These flows are not present on the west coast, being ponded down the depositional slope to the east. Similarly, interbedded volcanic debris flow units (e.g., at Hov) and volcaniclastic shale units (e.g., at Froðba) are recorded between lava flows in the east of the island (Fig. 1), but not seen in the west.
Vent proximal volcanic and volcaniclastic deposits occur in the northeast of Suðuroy, sourced from the Stapin Vent (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970; Passey & Jolley Reference Passey and Jolley2009). The pyroclastic deposits of this vent form an annular structure, which elevated the area above the surrounding Beinisvørð Formation surface, a topography inherited during deposition of the Prestfjall Formation. Exposures inland are limited, but photogrammetry interpretation of the west coast has shown no evidence of significant erosion at the base of the Prestfjall Formation (Passey et al. Reference Passey, Vosgerau, Sørensen and Vautravers2016). As a consequence of this, deposition of the Suðuroy Prestfjall Formation lacustrine system was controlled by this landscape inheritance. Within the lacustrine Prestfjall Formation, repeated cycles in palynofloral ecology, geochemical ratios, and sedimentology, indicate that longer duration environmental forcing occurred (Fig. 6). There is no evidence that these cycles were related to lava field processes, but rather that they reflect orbital or solar forcing controlling lacustrine water depth via precipitation. Assuming 21 k obliquity or 40 ky eccentricity cycling (Berger Reference Berger1988), this would indicate that the lacustrine record is 105–205 ky in duration.
Lateral facies variation between the uppermost units of the Prestfjall Formation and the basal Hvannhagi Formation/Malinstindur Formation is suggested by palynological data from Stapin, fossil wood in the lowermost Hvannhagi tuffs at Holið í Helli, and raised residual MgO values in the upper Prestfjall Formation of the KOLBHI-04 borehole. Incised channels filled with normally graded conglomerates and sandstones are linked here to overbank reddened clay stones which indicate localised erosion following deposition of the Prestfjall Formation in northwest Suðuroy. Elsewhere, the Prestfjall Formation is overlain by apparently conformably bedded volcanic, or thin volcaniclastic rocks of the Hvannhagi and Malinstindur formations.
The syn-eruptive facies of the Hvannhagi Formation on Suðuroy include primary pyroclastic tuffs, conglomerates, sandstones, and shales. These, primary olivine-phyric and welded tuffs confirm the linkage between the Hvannhagi Formation and the olivine-phyric lavas of the basal Malinstindur Formation (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970; Hald & Waagstein Reference Hald and Waagstein1984; Waagstein Reference Waagstein, Morton and Parson1988; Larsen et al. Reference Larsen, Waagstein, Pedersen and Storey1999; Passey & Jolley Reference Passey and Jolley2009).
Exposures of the base of the Malinstindur Formation are infrequent on Suðuroy, and they are often poor or invaded by numerous dolerite sills. Previous authors (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970; Waagstein Reference Waagstein, Morton and Parson1988; Passey & Jolley Reference Passey and Jolley2009) identified a transitional zone where lower Malinstindur Formation lava flow lobes are interbedded with volcaniclastic rocks attributed to the Hvannhagi Formation. This is supported by the succession recorded in the KOLBH1-04 borehole, where these volcaniclastic interbeds occur in the Malinstindur Formation above the Hvannhagi Formation.
6. Conclusions
The Hvannhagi Formation of the FIBG is a rare deep time example of a sector collapse-generated debris flow deposit associated with a low-angle shield (3–5° slope). Many recorded modern and deep time examples are from more steeply sloped stratovolcanoes in continental settings, dissimilar to the Paleogene rifting margin of the northeast Atlantic (Waresback & Turbeville Reference Waresback and Turbeville1990; Smith Reference Smith, Fisher and Smith1991; Smith & Lowe Reference Smith, Lowe, Fisher and Smith1991; Dufresne et al. Reference Dufresne, Siebert, Bernard, Roverato, Dufresne and Procter2021; Zernack Reference Zernack, Roverato, Dufresne and Procter2021). Thick lahar successions in western Vágar were deposited in two main phases, isopach mapping indicating an origin up to >40 km north of the lahar deposits (Fig. 15; Cronin et al. Reference Cronin, Lecointre, Palmer and Neall2000). These lahars were partly confined to a shallow north northeast–south southwest orientated topographical low, reflecting continued influence of late Beinisvørð Formation topography. A synchronous vent structure (a Malinstindur Formation rejuvenated Stapin Vent) northeast of Suðuroy (Fig. 16) sourced volcaniclastic deposits from pyroclastic eruptions. These eruptions initiated before the termination of Prestfjall Formation lacustrine shale deposition, suggesting that volcanism slowed considerably at the end of Beinisvørð Formation eruption, but did not stop completely.
The depositional geometries of the Hvannhagi Formation across the Faroe Islands exhibited no overriding single controlling factor. Proximity to volcanogenic source and eruptive morphologies interacted with syn-eruptive brittle fracturing and the inherited Beinisvørð–Prestfjall Unconformity surface to create a complex topography. This constrained debris flows, fluvial channels, and the interdigitating, onlapping lava flows of the Malinstindur Formation. The role of low-angle shield volcanoes in constraining the architecture and facies of sedimentary depositional systems is shown clearly. This study also highlights for the first time that shield volcanoes were coeval with and responsible for eruption of part of the Beinisvørð Formation, previously regarded as being composed solely of simple lava flows sourced from fissure erruptions (Passey & Bell. Reference Passey and Bell2007; Passey & Jolley Reference Passey and Jolley2009). Exposures on Vágar and Mykines confirm the presence of a low-angle shield to the north which was previously identified from gravity data. Additional Beinisvørð Formation compound flow fields sourced from low-angle shields may exist beneath the extensive cover of later fissure sourced simple flows. This field observation also implies that the dominantly simple flow characteristic of the exposed Beinisvørð Formation is attributable to it being distal to erruption source.
The oldest volcaniclastic strata of the Hvannhagi Formation exposed in northeast Suðuroy are shown to be partly laterally equivalent to the upper beds of the Prestfjall Formation. Hvannhagi Formation volcaniclastic units subsequently also interdigitated with the prograding Malinstindur Formation compound flow field (Rasmussen & Noe-Nygaard Reference Rasmussen and Noe-Nygaard1970; Passey & Jolley Reference Passey and Jolley2009). To the north, eruption of localised Malinstindur Formation compound lava flows is recorded in topographical lows prior to deposition of the Hvannhagi Formation on Tindhólmur. The lack of widespread erosion at the base of the Prestfjall Formation and the diachronism of the Prestfjall–Hvannhagi and Hvannhagi–Malinstindur formations’ boundaries all indicate a broad continuum of depositon and eruption. The pre-rift to syn-rift transition represented by the Beinisvørð to Malinstindur formations interval is stratigraphically equivalent to the the Flett Unconformity recorded on the eastern margin of the Faroe–Shetland Basin (Shaw Champion et al. Reference Shaw Champion, White, Jones and Lovell2008; Hartley et al. Reference Hartley, Roberts, White and Richardson2011; Jolley et al. Reference Jolley, Millett, Schofield, Broadley and Hole2021). This major unconformity rejuvenated the earlier Upper Thanetian Unconformity surface, resulting in the incision of basin-facing dendritic draininge systems in response to rift flank uplift. Within the stratigraphically equivalent Beinisvørð to Malinstindur formations interval, there are complex facies changes between lava flow field, fluvial, lacustrine to primary volcaniclastic, and lahar. These deposits, while locally sculpted as part of a complex lava field terrain, reperesent the correlative conformity of the major Flett Unconformity.
Cessation of volcanism at the end of Beinisvørð Formation eruption in the FIBG area, was regarded as a response to a transient reduction in mantle temperature (Millett et al. Reference Millett, Hole, Jolley, Passey and Rossetti2020). The resultant volcanic hiatus is however, short or a period of very low eruption frequency less than the duration of Prestfjall Formation deposition. Characterised by sedimentary depositional systems, the eruptive hiatus or slow-down was followed by an initially episodic resumption of volcanism focussed on multiple low-angle shields. This style of volcanism replaced the fissure-fed and distal to source lava plain facies of the Beinisvørð Formation. It is this more varied, local sourcing that gives rise to the complex interactions of VDA's, lahars and pyroclastic deposits forming the Malinstindur Formation, and the interdigitation of these deposits with the uppermost Prestfjall Formation and oldest Malinstindur Formation lava flows.
7. Supplementary material
Supplementary material is available online at https://doi.org/10.1017/S1755691022000056
8. Acknowledgements
The photogrammetry survey for this study was carried out as part of the Cambridge Arctic Shelf Project (CASP)'s Flood Basalt Impact on Hydrocarbon Systems Project 2014–16 and CASP's sponsors' financial support is gratefully acknowledged as well as the Geological Survey of Denmark and Greenland’ contribution. We would also like to thank Atlantic Airways for their help in the acquisition of the oblique photographs used in preparation and interpretation of the models. Davie Brown and an anonymous reviewer are thanked for helpful reviews of the manuscript. John Howell is thanked for discussions on photogrammetry and decompaction of the sedimentary rock record.