Abstract

Numerous studies have investigated the geodynamic history and lithological composition of the Proterozoic basement, Caledonian nappes, and Devonian extensional basins and shear zones onshore west Norway. However, the offshore continuation of these structures, into the northern North Sea, where they are suspected to have influenced the structural evolution of the North Sea rift, is largely unknown. Existing interpretations of the offshore continuation of Caledonian and Devonian structures are based on simple map-view correlations between changes in offshore fault patterns and pronounced onshore structures, without providing evidence for the presence, nature, and geometry of offshore, basement-hosted structures.

By integrating three-dimensional (3-D) seismic, borehole, and onshore geological and petrophysical data, as well as two-dimensional (2-D) forward modeling of gravity and magnetic data, we reveal the structural architecture and composition of the crystalline basement on the Måløy Slope, offshore west Norway. Based on 3-D mapping of intrabasement reflection patterns, we identified three basement units that can be correlated with the Caledonian thrust belt, and the major Devonian Nordfjord-Sogn detachment zone, located only 60 km to the east, onshore mainland Norway. Similar to that observed onshore, offshore crystalline basement of the Proterozoic basement (Western Gneiss Region) and allochthons is folded into large-scale antiforms and synforms. These units are separated by the strongly corrugated Nordfjord-Sogn detachment zone. Our analyses show that different types of crystalline basement can be distinguished by their seismic reflection character, and density and magnetic properties. We speculate that the main causes of the observed intrabasement reflectivity are lithological heterogeneities and strain-induced structures such as shear and fracture zones. Our interpretation of the architecture of crystalline basement offshore west Norway has important implications for the location of the suture zone between Baltica and Laurentia.

INTRODUCTION

The discovery of hydrocarbon plays related to fractured and weathered crystalline basement (e.g., Nelson, 2001; Cuong and Warren, 2009; Hartz et al., 2013) has driven interest in the nature and three-dimensional (3-D) geometry of intrabasement structures. Furthermore, an increasing number of studies are demonstrating how the composition and structure of crystalline basement can influence the geometry and growth of normal faults developing within overlying cover strata, in addition to the physiography of the associated rift basin (e.g., Morley et al., 2004; Paton, 2006; Morley, 2009, 2010; Kirkpatrick et al., 2013; Phillips et al., 2016). Despite their importance, the geometry and composition of crystalline basement are largely neglected or simply considered homogeneous in analogue and numerical models (e.g., Gupta et al., 1998; Cowie et al., 2000; Huismans and Beaumont, 2008); as such, these properties are not explicitly represented in tectono-stratigraphic models of rift development (e.g., Gawthorpe and Leeder, 2000). A better understanding of basement heterogeneity will, however, allow us to: (1) construct more realistic physical and numerical models of rift tectono-stratigraphic development; (2) understand how basement properties control structural inheritance in a range of tectonic settings (e.g., rifts, orogenic belts); and (3) better account for heat-flow variations and impacts in numerical models of basin geodynamics.

The northern North Sea rift is considered to have formed above a heterogeneous substrate with, in places, structures inherited from previous Proterozoic and Paleozoic contractional and extensional events (e.g., Bartholomew et al., 1993; Færseth et al., 1995; Roberts et al., 1995; Færseth, 1996; Doré et al., 1997; Fossen et al., 2016; Fazlikhani et al., 2017). These structures include strongly deformed Precambrian basement, allochthonous rocks, and large-scale Devonian extensional shear zones and faults that are exposed onshore west Norway, immediately east of the rift (e.g., Andersen and Jamtveit, 1990; Fossen, 1992; Andersen and Andresen, 1994; Milnes et al., 1997; Eide et al., 1997, 1999; Krabbendam and Dewey, 1998; Fossen, 2000; Osmundsen and Andersen, 2001; Johnston et al., 2007). Despite intensive study of the formation and collapse of the Caledonides onshore west Norway, in addition to northern Scotland, Orcadia and West Orkney Basin, and the Greenland Caledonides (e.g., Brewer and Smythe, 1984; McGeary and Warner, 1985; Cheadle et al., 1987; Coward et al., 1989; Snyder, 1990; Bird et al., 2015; Norton et al., 1987; Coward, 1990; Wilson et al., 2010; Andresen et al., 2007; Gilotti and McClelland, 2008), the onshore geology can be projected offshore for only a short distance, meaning the structural architecture and composition of the crystalline offshore basement remain poorly understood.

Most studies base their interpretation of offshore pre-existing basement structures on widely spaced, relatively low-resolution, two-dimensional (2-D) seismic reflection and potential field data, as well as regional-scale, principally map-view correlations of offshore fault geometries with onshore basement structures, without providing evidence for the nature, location, and geometry of intrabasement structures (e.g., Færseth et al., 1995; Færseth, 1996; Doré et al., 1997; Fossen et al., 2016). Few studies have attempted to decipher the composition and location of intrabasement structures beneath the North Sea rift using geophysical data (e.g., Reeve et al., 2014; Fossen et al., 2016; Phillips et al., 2016; Fazlikhani et al., 2017). Using borehole, gravity, and magnetic data, Smethurst (2000) presented an interpretation of basement lithology and structures offshore west Norway, but the 3-D structural architecture of the crystalline basement remains unclear. Furthermore, due to their inherent ambiguity, gravity and magnetic data alone may lead to false interpretations if not constrained by a sound subsurface model and petrophysical data. Recently, Fazlikhani et al. (2017) showed that the crystalline basement beneath the northern North Sea exhibits a varying seismo-acoustic character that can be divided into three primary seismic facies that can be mapped throughout the basin. Despite this advance, Fazlikhani et al. (2017) only used 2-D seismic reflection data, and they did not constrain their interpretation with other available data (e.g., potential field). Borehole data could support seismic reflection–based interpretations by providing direct constraints on basement composition and, if sufficiently closely spaced, structure. However, many boreholes terminate above crystalline basement, with those that do sample basement being widely spaced. This, combined with poor seismic data resolution, complicates the imaging and interpretation of complex 3-D intrabasement structures.

Intrabasement seismic reflectivity, such as that identified and mapped by Fazlikhani et al. (2017), may be the result of compositional and structural heterogeneities within the crust, and the type and amount of ductile and brittle deformation that an area experienced (e.g., Smythe et al., 1982; Jones and Nur, 1984; Wang et al., 1989; Goff and Holliger, 2003; Smithson and Johnson, 2003; Meissner et al., 2006; Lyngsie et al., 2007; Bongajum et al., 2012; Ahmadi et al., 2016; Hedin et al., 2016). Although intrabasement structures have been described from regional, deep-imaging 2-D seismic reflection profiles (e.g., Wang et al., 1989; Allmendinger et al., 1987; Kern and Wenk, 1990; Reston, 1996; Cook et al., 1999; Abramovitz and Thybo, 2000; Juhojuntti et al., 2001; Simancas et al., 2003; Torvela et al., 2013; Ruiz et al., 2017), few studies use borehole data (e.g., Christensen and Szymanski, 1988; Zhang et al., 2005; Hedin et al., 2012, 2014, 2016; Lorenz et al., 2015), synthetic seismogram modeling (e.g., Fountain et al., 1984; Hurich et al., 1985; Reeve at al., 2014; Phillips et al., 2016; Kennett et al., 2017; Veludo et al., 2017), and/or gravity and magnetic forward modeling (e.g., Ebbing et al., 2006; Lyngsie and Thybo, 2007; Ritzmann and Faleide, 2007; Marello et al., 2010; Gernigon and Brönner, 2012; Lassen and Thybo, 2012; Mazur et al., 2015; Olyphant et al., 2017) to help further constrain a largely seismic reflection–based interpretation.

Determining the sources of crustal reflectivity, and mapping crystalline basement structure and composition require a holistic approach, drawing on a range data types, methods, and scales of observation. On the Måløy Slope, offshore west Norway (Fig. 1A), intrabasement structures are well imaged in seismic reflection data due to relatively shallow burial of the pre-Devonian metamorphic basement beneath a thin (<3.5 km) and stratigraphically simple sedimentary cover (Fig. 2). Due to its complex but well-known tectonic history, which yielded a wide variety of crystalline rock types and associated structures, the study area serves as an ideal natural laboratory in which to constrain the geometry, investigate the origin, and establish the regional importance of intrabasement reflectivity. We used borehole and core data, conventional 2-D and 3-D seismic reflection data, gridded potential field data, and 2-D forward modeling of gravity and magnetic data to constrain the nature, composition, and 3-D architecture of the crystalline basement on the Måløy Slope.

Our results show that intrabasement reflections correlate with major structures and basement units comprising the Proterozoic basement and allochthonous nappes of the Caledonian thrust belt and associated Devonian extensional shear zones. By linking onshore and offshore structures, we can reveal the 3-D geometry of large-scale intrabasement structures along the western part of the Norwegian margin. Using the depth and architecture of these intrabasement structures, we briefly discuss the implications of our interpretation for the highly uncertain location of the suture zone between Baltica and Laurentia. In addition, based on the known geological history, we suggest several plausible origins for intrabasement reflectivity mapped within the study area. We show that the primary causes of basement reflectivity are lithological heterogeneities, and associated interference and scattering effects of the seismic wave field, as well as strain-induced structures such as shear and fracture zones. Based on our results, we consider that the approach followed in this study can be applied in other subsurface studies to help us more confidently interpret the origin of intracrystalline basement reflectivity in other sedimentary basins.

GEOLOGICAL FRAMEWORK

Study Area

The study area is located along the Måløy Slope, offshore west Norway, comprising an area of ∼3600 km2 (Fig. 1A). The western boundary of the study area is defined by the Sogn graben; the southern and northern boundaries are approximately coincident with the coastal outlets of the Sognefjord and Nordfjord, respectively, two major fjords onshore western Norway. Onshore, Proterozoic basement (Western Gneiss Region [WGR]), remnants of the Caledonian thrust belt, a major Devonian extensional shear zone, and several intermontane Devonian sedimentary basins are preserved (e.g., Andersen and Andresen, 1994; Milnes et al., 1997; Johnston et al., 2007, Vetti and Fossen, 2012). These basement units are suspected to form the metamorphic, crystalline basement offshore west Norway (e.g., Smethurst, 2000; Reeve et al., 2014), but their offshore continuation and 3-D geometry remain unclear (Færseth et al., 1995; Smethurst, 2000).

Regional Tectonic History

The crystalline basement in the northeastern North Sea is largely the result of terrane accretion during the Sveconorwegian and Caledonian orogenies (Ziegler, 1975; Fossen, 1992; Bingen et al., 2008; Gee et al., 2008). The Caledonian orogeny culminated in SE-directed thrusting of allochthonous and parautochthonous nappes onto the subducting margin of Baltica, which is commonly termed the Western Gneiss Region (WGR; ca. 403 Ma; Andersen et al., 1991; Fossen, 1992; Fossen and Dunlap, 1998; Root et al., 2004; Gee et al., 2008). Early to Middle Devonian postorogenic extension (ca. 400 Ma) resulted in W- to WNW-directed backsliding of the nappe pile along a reactivated low-angle Caledonian décollement zone (Fossen, 1992, 2000; Fossen and Dunlap, 1998), the formation of large-scale WNW-dipping extensional shear zones (Nordfjord-Sogn detachment zone [NSDZ]), exhumation of the WGR, and the formation of Devonian continental basins (Andersen, 1998; Fossen, 2000; Osmundsen and Andersen, 2001). Simultaneously, N-S shortening, related to the interplay of the NSDZ and left-lateral movement along the Møre-Trøndelag fault complex (MTFC) further north, resulted in transtensional folding of the WGR, the detachment zone, and the Caledonian allochthons and Devonian deposits (Chauvet and Séranne, 1994; Krabbendam and Dewey, 1998; Dewey and Strachan, 2003; Osmundsen and Andersen, 2001; Fossen et al., 2013). Offshore, crystalline basement along the Måløy Slope is displaced by Middle to Late Jurassic normal faults and is unconformably overlain by Mesozoic strata (e.g., Færseth et al., 1995; Sømme et al., 2013; Reeve et al., 2015).

Onshore Geological Framework and Basement Units

The WGR is located in the footwall of the NSDZ (Figs. 1C and 2) and is composed of highly deformed Proterozoic (1750–900 Ma; Engvik and Andersen, 2000) granodioritic and migmatitic gneisses, and strongly foliated amphibolites, augen-gneisses, and anorthosites (Hacker et al., 2010). Meter- to kilometer-scale, E-W–elongated lenses of granulites, gabbros, and eclogite bodies are common. Devonian extension and exhumation of the WGR strongly reworked and obliterated many Caledonian structures, with the associated rock foliation being generally subparallel to the NSDZ, following gently W-plunging, WGR-cored antiforms and synforms (Dewey et al., 1993).

The allochthonous nappe stack is located in the hanging wall of the NSDZ and is further divided into the Lower, Middle, and Upper Allochthons (Fig. 1; Roffeis and Corfu, 2014). The Lower Allochthon consists of gneisses and weakly metamorphosed sandstones and graywackes of the Baltica craton (Hacker and Gans, 2005; Johnston et al., 2007). The Middle Allochthon is characterized by orthogneisses, granulites, gabbros, and anorthosites, which formed the basement and cover sequences of the western Baltica margin (Andersen et al., 1994; Krabbendam et al., 2000). The Upper Allochthon is composed of ophiolitic and island-arc–derived material that is exotic to Baltica (Andersen et al., 1994; Johnston et al., 2007). The total thickness of the allochthons is highly variable, but their minimum present-day thickness is estimated to be ∼4 km (Johnston et al., 2007).

The rocks of the allochthons are unconformably overlain by Devonian metasediments of the Hornelen, Håsteinen, Kvamshesten, and Solund basins (Fig. 1; Osmundsen and Andersen, 2001). The Devonian metasediments generally comprise basin margin–attached, upward-coarsening packages of alluvial-fan sandstones and conglomerates of allochthonous provenance (Cuthbert, 1991; Osmundsen and Andersen, 2001; Templeton, 2015). In the basin axes, upward-coarsening, up to 200-m-thick, alluvial plain and lacustrine successions occur (Steel, 1976; Steel et al., 1977; Steel and Aasheim, 1977). The Devonian deposits experienced anchizone to lower-greenschist-facies metamorphism (Seranne and Seguret, 1987; Torsvik et al., 1986; Sturt and Braathen, 2001), corresponding to sediment burial depths of 5–13 km (Svensen et al., 2001) and maximum temperatures of 300 °C (±50 °C).

The NSDZ is the principal shear zone in the study area. This structure formed to accommodate Devonian extensional tectonics and associated exhumation of the WGR (e.g., Fossen, 2000). The NSDZ juxtaposes high-pressure rocks of the WGR in its footwall against low-grade metamorphic rocks of the allochthons and Devonian basins in its hanging wall (Fig. 1C). The shear zone is folded into a series of kilometer-scale, E-W–trending antiforms and synforms that plunge to the west with varying dips (∼5°–25° to W-WNW; Fig. 1C; Norton, 1987; Labrousse et al., 2002). The 2–6-km-thick (Norton, 1986, 1987; Johnston et al., 2007), W-dipping NSDZ is composed of various mylonitic rock types that indicate varying degrees of deformation under retrograde amphibolite- to greenschist-facies conditions; these rock types imply that the detachment initially formed at crustal depths between 20 km (Wilks and Cuthbert, 1994) and 30–40 km (Johnston et al., 2007).

DATA SETS AND METHODS

Borehole, Core, and Petrophysical Data

Onshore and offshore core, borehole, and petrophysical data were used to constrain the composition and physical properties of crystalline rocks along the Måløy Slope. Eleven offshore boreholes penetrate crystalline basement (Norwegian Petroleum Directorate [NPD] Web site: www.npd.no/en/Topics/Wells/; Knape, 1996; Fig. 1A), with three providing cores extracted from crystalline basement rock; we analyzed these cores for basement mineral composition and rock type (Fig. 3; Table 1). Using publicly available well reports, we compiled detailed information about different basement lithologies present offshore west Norway (Table 1). In addition, we collated petrological and petrophysical data for various rock types within each basement unit onshore west Norway from literature (Skjerlie, 1985; Norton, 1986; Dewey et al., 1993; Andersen et al., 1994; Osmundsen and Andersen, 1994; Osmundsen et al., 1998; Bryhni and Lutro, 2000a, 2000b; Lutro and Bryhni, 2000; Krabbendam et al., 2000; Ebbing et al., 2006; Slagstad et al., 2008; Vetti, 2008; Olesen et al., 2010; Reynisson et al., 2010; Hedin et al., 2014), and we screened ∼1200 published onshore rock samples and seven boreholes for rock type, geological unit, density, and susceptibility (NGU database, http://aps.ngu.no/kart/geofysikk/; Table 2). Mean values and standard deviations of density and susceptibility were calculated for the 20 most common and statistically significant onshore rock types (Table 2). These density and susceptibility data provided important constraints for our assessment of the basement composition, and the 2-D forward modeling of gravity and magnetic data drawing on these data.

Seismic Reflection Data

The study area is covered by 2-D seismic reflection profiles (imaging to 7–9 s two-way traveltime [TWT]), and four partly overlapping, post-stack time-migrated 3-D seismic surveys of moderate to good quality (imaging to 5–7 s TWT; Fig. 1A). These data were processed in zero-phase and are presented in reversed polarity as adopted by the Society of Exploration Geophysicists (SEG) convention, where red colors (trough) indicate a downward increase in acoustic impedance. The dominant frequency within intrabasement reflections, extracted from the 3-D seismic data, ranges from 10 to 30 Hz. To improve the signal-to-noise ratio, 3-D seismic volumes were frequency filtered to reduce coherent and random noise (for details about the noise filtering process, see GSA Data Repository Item DR11). To estimate the true depth and geometry of basement structures, seismic data were depth converted using a 3-D velocity model (see Data Repository Item DR2). The top of the crystalline basement, directly tied to well data, is characterized by a high-amplitude trough-peak reflection event (e.g., Fig. 2). Seismic interpretation within crystalline basement was based on the mapping of individual continuous intrabasement reflections, broader changes in seismic facies, and distinct geometrical relationships (e.g., boundaries) between seismic reflection patterns. To assess whether the observed basement reflectivity is real, and to largely rule out seismic artifacts, we inspected the data for possible multiples, diffractions, and migration artifacts, and we compared the acoustic character of identical structures imaged by different overlapping seismic data sets.

2-D Forward Modeling of Gravity and Magnetic Data

We analyzed gridded Bouguer gravity and aeromagnetic data, and we undertook 2-D forward modeling of gravity and magnetic data to assess the physical properties, and thus likely composition, age, and spatial distribution, of basement rocks along the Måløy Slope. The Bouguer gravity and aeromagnetic data presented in this study (Fig. 6) were extracted from data compiled by the Norwegian Geological Survey (NGU; Skilbrei et al., 2000; Olesen et al., 2010) and have a grid resolution of 500 m. Aeromagnetic data were reduced-to-the-pole using the International Geomagnetic Reference Field (IGRF). The 2-D forward modeling was conducted in the depth-domain in association with the depth-converted interpretation of a seismic line (Fig. 7) located along the model profile; this seismic line provided geometrical constraints on the size and position of modeled basement bodies. We varied density and magnetic properties along the modeled line within the range of known values for each basement rock type onshore (Table 2). Using results from gravity and magnetic data analysis (Olesen et al., 2010; Ebbing et al., 2012), in addition to seismic tomography and refraction profiling studies (Stratford et al., 2009; Stratford and Thybo, 2011; Maupin et al., 2013; Mjelde et al., 2016; Kvarven et al., 2016) from onshore and offshore mid-Norway, we set the depth, density, and susceptibility of the lower crust and mantle to constant values of 23 km, 2.95 g/cm3, and 0.01 SI, and 30 km, 3.3 g/cm3, and 0 SI, respectively. To avoid a level of model complexity that we could not constrain with our observations from seismic and onshore geological data, we deliberately chose a simplified modeling approach that assumed no changes in the topography or physical properties of lower-crustal and mantle rocks along our 135-km-long model profile. As a result, our model represents admittedly one of several plausible solutions that could explain the origin and properties of upper-crustal gravity and magnetic signatures, which we here argue readily relate to our understanding of the crystalline basement geology onshore western Norway.

Average densities for Mesozoic and Cenozoic sediments offshore were derived from well-log data from 12 boreholes; we kept these values constant in each model. Additional details about the modeling procedure and its uncertainties are summarized in Data Repository Items DR3 and DR4.

COMPOSITION OF THE CRYSTALLINE BASEMENT

Publicly available data for basement-penetrating wells along the Måløy Slope (Table 1) suggest that the crystalline basement is mainly composed of pre-Devonian gneisses and greenschists/mica schists. In borehole 36/1–1 (Fig. 3A), we describe the basement as a biotite augen gneiss, a rock type typifying the WGR onshore (Hacker et al., 2010). Cores from boreholes 35/3–4 and 36/4–1 (Figs. 3B and 3C) recovered banded gneiss, which is characterized by dipping, millimeter- to centimeter-scale, alternating bands of lighter, feldspar and quartz-rich bands, and darker, biotite-rich bands. Banded or mylonitic gneisses signify pronounced ductile deformation, which is widespread within both the WGR and allochthons. However, based on the mineral composition, and geographical (i.e., offshore Bremangerlandet; Fig. 1) and stratigraphic position of the wells relative to the onshore geology, we infer that boreholes 35/3–4 and 36/4–1 most likely penetrated basement related to the Lower or Middle Allochthon.

With a penetration depth of 233 m, borehole 35/3–2 is the longest drilled section of crystalline basement within the study area and therefore provides important insights into the basement compositional heterogeneity. The crystalline basement is reported to consist of foliated mica schist and gneiss (Table 1); alternating mica-rich and quartz-feldspar–rich bands likely indicate a metamorphic texture similar to the cores from 35/3–4 and 36/4–1 (Figs. 3B and 3C). Furthermore, published well logs show pronounced fluctuations of the P-wave velocity and resistivity values in the uppermost 30 m of the basement. These fluctuations, together with increased neutron values, may indicate that the upper part of the crystalline basement is likely weathered and has fluid-filled fractures (Bartetzko et al., 2005; NPD Web site: www.npd.no/en/Topics/Wells/). Based on these observations, the rock types in borehole 35/3–2 appear comparable to the banded, mylonitic gneisses drilled in wells 35/3–4 and 36/4–1 (Figs. 3B and 3C) and, hence, are most likely part of the Lower or Middle Allochthon.

SEISMIC REFLECTION CHARACTER OF CRYSTALLINE BASEMENT

Basement Seismic Facies

Based on analysis of reflection continuity, amplitude character, and geometrical relationships between individual reflections and reflection packets, we identified four distinct seismic facies (SF1–SF4) and three common styles of seismic facies relationships (SFRa, SFRb, and SFRc; Figs. 4 and 5). The seismic character, 3-D geometry, and crustal position of each intrabasement seismic facies and seismic facies relationship allowed for a direct correlation with the main basement units exposed onshore. The characteristics of our identified intrabasement seismic facies are similar to the ones described in Fazlikhani et al. (2017), suggesting that the acoustic properties of different parts of the crystalline basement are similar throughout the basin and relatively independent of the properties of seismic reflection data. We here emphasize that amplitude strength may not be directly indicative of rock type, as the seismic data of different vintages processed by different companies will have different automatic gain control applied. Furthermore, an attempt to differentiate seismic facies based on their frequency character was largely inconclusive due to strong frequency variations throughout the interpreted crystalline basement, and unknown data seismic acquisition and processing parameters. Nonetheless, we infer that locally observed higher-frequency seismic facies at depth likely represent rock types characterized by a higher P-wave velocity than the surrounding rocks (e.g., Yilmaz, 1987).

Seismic Facies Analysis

Seismic facies SF1 is characterized by moderately continuous, subparallel, medium- to high-amplitude reflections (Fig. 4A). This reflection pattern is reminiscent of that typically characterizing layered sedimentary successions, which, together with its stratigraphic position in the uppermost part of the basement (Fig. 7), suggests SF1 may represent sedimentary sequences within Devonian metasedimentary rocks. Coarsening-upward cycles within the Devonian basins onshore are locally within the range of seismic resolution (50–200 m; Steel, 1976; Steel et al., 1977; Steel and Aasheim, 1977). Devonian metasediments, however, are not penetrated and thus are not confirmed by any of the offshore wells. An alternative interpretation of the observed reflectivity is that it represents compositional layering within allochthon-related metamorphic rocks.

Seismic facies SF2 consists of both low-amplitude chaotic, and medium- to high-amplitude folded and discontinuous seismic reflections (Fig. 4B), suggesting strongly deformed material. The observed folds exhibit an upright geometry, wavelengths of 2–5 km, and amplitudes of up to hundreds of meters. Based on its crustal position, correlation with well and core data (Fig. 7), and geometric similarity to structures (i.e., folds) observed in the allochthons onshore (e.g., Andersen and Jamtveit, 1990; Osmundsen and Andersen, 1994; Bryhni and Lutro, 2000b), SF2 is suspected to represent rocks of the allochthonous nappes.

Seismic facies SF3 is characterized by a 2–4-km-thick, W-dipping zone of moderately continuous, subparallel, high-amplitude reflections that are defined by high-amplitude trough-peak-trough wave trains (Fig. 4C). The frequency within SF3 seems to be generally higher than in other seismic facies, suggesting the presence of a higher-velocity rock type. Analogous seismic facies are commonly described in deep seismic reflection studies, where they are interpreted as mylonites within ductile shear zones (e.g., Wang et al., 1989; Fountain et al., 1984; Shaocheng et al., 1993; Reeve et al., 2014; Hedin et al., 2014, 2016; Phillips et al., 2016). Mylonites have average P-wave velocities of 6.9 km/s (e.g., Ji et al., 1993; Law and Snyder, 1997), explaining the observed high-frequency character. We therefore interpret SF3 in a similar way, as a thick zone of mylonites.

SF4 is defined by discontinuous, oppositely dipping, high- to very high-amplitude reflections that crosscut or terminate against each other (Fig. 4D). This seismic facies is only observed in the deepest seismically imaged parts of the basement (>12 km deep). Similar seismic facies are described from midcrustal levels within the Paleoproterozoic orogenic belt in southern Finland, where they are interpreted as kilometer-scale, conjugate S-C structures that form a penetrative fabric (Torvela et al., 2013). Andersen and Jamtveit (1990), Dewey et al. (1993), and Osmundsen and Andersen (1994) also described S-C deformation fabrics within the WGR onshore, with these fabrics becoming more common upward toward the NSDZ. We therefore interpret SF4 to represent a large-scale ductile deformation fabric within the WGR.

Seismic Facies Relationships

Abrupt terminations of reflections, and sudden changes from one seismic facies to the other are interpreted as boundaries between different basement units; we used these spatial relationships to reconstruct the 3-D geometry and structural architecture of the crystalline basement (Fig. 5). Unit boundaries are defined by steeply dipping, nonreflective planes that are interpreted as brittle faults (SFRa; Fig. 5A), or by a distinct change in the overall geometry and dip direction between two seismic facies (SFRb; Fig. 5B). The latter boundary type is exclusively associated with SF3 where it is overlain by SF2 above or underlain by SF4. Since we suspect SF3 to represent mylonites, we interpret this contact as a shear contact. A third type of seismic facies relationship was identified across an ∼10-km-wide zone in the central part of the study area (SFRc; Figs. 5C and 7). Here, subhorizontal, medium-amplitude reflections of SF1 are crosscut and in parts slightly offset by steeply S- to SW-dipping, high-amplitude reflections that are similar to SF3. Furthermore, the position of the steeply dipping reflections can be linked to small-scale normal faulting at top basement level. Hence, we infer that the zone where we observe SFRc (Fig. 7) comprises an area of intense deformation, i.e., a shear or fracture zone.

GRAVITY AND MAGNETIC SIGNATURE OF CRYSTALLINE BASEMENT

To support our borehole and seismic reflection–based interpretation of the structural architecture and composition of crystalline basement offshore west Norway, we used density and magnetic susceptibility data for the most common onshore basement rock types (Table 2), as well as offshore Bouguer gravity and aeromagnetic data (Fig. 6).

The density and magnetic susceptibility data compiled in Table 2 indicate that Devonian metasediments and rocks of the allochthons and WGR have similar mean densities. The WGR can, however, be distinguished from allochthons and Devonian strata by its strong magnetic character (Table 2). Furthermore, Devonian sediments exhibit very similar magnetic properties to their allochthonous source rocks; these two units may thus be indistinguishable in gravity and magnetic data (Table 2).

Figure 6 shows gridded Bouguer gravity and aeromagnetic data covering the offshore part of the study area. Two positive gravity anomalies are observed: (1) a pronounced N-trending, positive gravity high (maximum 50 mGal) extending along the coastline of west Norway (anomaly A1), and (2) a WSW-trending, positive anomaly with maximum values of 34 mGal (anomaly A2) located west of Dalsfjord (Fig. 6A). Anomaly A2 is further characterized by a strong, short-wavelength positive magnetic anomaly (155 nT) indicating a dense and magnetic source body at relatively shallow depth (Fig. 6B). According to completion reports, boreholes 36/7–1 and 35/9–1, which are located above the anomaly, drilled through granitic gneiss (typical for WGR; Hacker et al., 2010) and greenschist/gneiss (typical for allochthon; Andersen et al., 1994; Johnston et al., 2007), respectively (Table 1). Based on its position along strike of large-scale, WGR-cored antiforms observed onshore (Osmundsen and Andersen, 2001; Fig. 1B), and the folded reflections observed in seismic data (Fig. 7), we interpret anomaly A2 as a WSW-plunging WGR-cored antiform. The high gravity values of anomaly A1 point to a dense, shallow source body; with basement rocks being exposed on the islands offshore west Norway, it seems logical that anomaly A1 represents a crystalline basement high sub-cropping the seafloor. The steep, W-dipping gradients along the western boundary of anomaly A1 may be caused by a strong lateral density contrast between the crystalline basement-cored high in the east and Mesozoic and Cenozoic marine sediments in the west, the effect of which is enhanced by a N-S–striking, W-dipping, basin-bounding fault (fault Øygarden 2 in Bell et al., 2014; see also Fig. 6A).

A pronounced ENE-trending, ∼4-km-wide, linear positive magnetic anomaly (165 nT; anomaly A3; Fig. 6B) around the islands of the Florø area is located along strike of the broadly ENE-WSW–striking Haukå fault (HF) and Eikefjord fault (EF), which bound a gneissose, WGR-cored horst (Lutro and Bryhni, 2000). Magnetic anomaly A4 is located 20 km north of the Florø horst, SW of Bremangerlandet, and it is suspected to be associated with a gabbronoritic/dioritic intrusion that has been described from the island of Frøya (443 ± 4 Ma old Gåsøy intrusion; Hansen et al., 2002; see also Fig. 6B).

A subtle NNE-trending gravity low of 3 to −5 mGal in the northern Måløy Slope may represent locally deep-lying basement, although this inference is not supported by seismic data interpretation (Fig. 7), magnetic data, and onshore-offshore correlations. Instead, a strong, laterally extensive positive magnetic anomaly (anomaly A5; maximum 285 nT; Figs. 6A and 6B), covering approximately the same area, indicates the presence of either a basement structural high, similar to anomaly A2, or an intrusive body. The locations of ENE-trending, WGR-cored antiforms and synforms documented for the onshore part of the study area (Osmundsen and Andersen, 2001; Figs. 6A and 6B) support our interpretation that large-scale positive magnetic anomalies (i.e., anomaly A2) that have a similar trend and that are located along strike of the onshore structures are likely to be attributed to WGR-cored antiforms. Furthermore, biotite augen gneiss drilled in well 36/1–1 appears to confirm the presence of a WGR structural high in the northernmost part of the Måløy Slope.

STRUCTURAL ARCHITECTURE OF CRYSTALLINE BASEMENT ALONG THE MÅLØY SLOPE

Based on: (1) the spatial distributions of SF1–SF4 and SFRa–SFRc, and (2) the lithological and geometrical constraints derived from analysis of well-log, core, and gravity and magnetic anomaly maps, we mapped the four principal basement units (i.e., WGR, NSDZ, allochthons, and Devonian sediments) in three dimensions along the Måløy Slope. The interpreted gross structural architecture of the crystalline basement is illustrated on a N-S–trending seismic profile that connects the offshore with the onshore geology (Fig. 7; a large-scale version of Fig. 7 is presented in Data Repository Item DR5). Characteristic, moderately continuous, high-amplitude reflections of SF3 are clearest, separating discontinuous reflections of SF2 and SF4 (i.e., SFRb; Figs. 4B and 4D). Mylonitic rocks (i.e., SF3; Fig. 4C) form a 2–4-km-thick, highly reflective band that is folded into large-scale (maximum amplitude: ∼8 km, wavelength: ∼50 km), gently W-plunging, E-trending antiforms and synforms (Fig. 7). This thick mylonitic band is suspected to represent the offshore continuation of a major onshore shear zone, the NSDZ, which deforms rocks of the WGR and allochthons. Average dip angles of the shear zone mapped offshore are in agreement with published estimates from onshore studies (Norton, 1987; Labrousse et al., 2002; Johnston et al., 2007; Vetti, 2008; Templeton, 2015), varying between 10° and 17° on E-trending sections, and reaching up to 35° on the flanks of antiforms (Fig. 7C). The top of the WGR is approximated by the base of the NSDZ and locally by the top basement reflection (Fig. 7). The overall geometry of the WGR is thus characterized by two, gently W-dipping, antiforms in the northern and southern Måløy Slope, with these structures connected by an ∼35-km-wide synform resembling the large, folded basement structure inferred from gravity and magnetic anomaly maps (Fig. 6). Furthermore, boreholes 36/1–1 and 36/7–1 provide evidence that rocks attributed to the WGR indeed sub-crop top basement in areas where WGR-cored antiforms are expected to be present.

Consequently, the up to 5-km-thick portion of the basement contained within gently W-plunging synforms of the WGR and NSDZ likely consists of allochthonous material and possibly Devonian sediments (Fig. 1B). The dominant intrabasement reflection types in the central part of the study area are classified as SF2, which we interpret to represent allochthonous material. Lithologies typical for the Lower and Middle Allochthon are reported from six boreholes (Table 1), confirming the presence of allochthonous material at top basement level. Furthermore, we observed small- and large-scale internal folding and complex geometrical relationships between individual reflections and reflection packets, which likely indicate a high degree of structural complexity. These findings are consistent with onshore observations from exposed allochthonous units, reflecting the complex, protracted deformation history of the area (e.g., Andersen and Jamtveit, 1990; Andersen et al., 1994; Osmundsen and Andersen, 2001; Johnston et al., 2007).

The uppermost part of the crystalline basement locally comprises reflections of SF1, which may reflect either allochthons or metasediments of the Devonian basins. The ∼2-km-thick unit of subparallel, almost stratified reflections locally appears to unconformably overlie reflections of SF2 (Fig. 7), suggesting the presence of a depositional, albeit unconformable, contact between Devonian sediments and the allochthon. However, because none of the basement-penetrating boreholes on the Måløy Slope encountered Devonian sediments, an alternative source of the observed subparallel reflectivity must be considered (see below).

2-D FORWARD MODELING OF GRAVITY AND MAGNETIC DATA

The largest uncertainties associated with our interpretation of the seismic reflection and gridded gravity and magnetic anomaly data are associated with: (1) the presence or absence of Devonian metasediments offshore west Norway; and (2) whether changes in seismic facies and reflection patterns effectively represent different lithological units with particular density and magnetic properties. We therefore used 2-D forward modeling of gravity and magnetic data to mitigate these uncertainties, choosing a simple approach that drew on our observations of the basement geology onshore and offshore western Norway. Step-by-step details about the modeling procedure and the associated uncertainties are provided in Data Repository Items DR3 and DR4.

Figure 8 shows the final density and magnetic forward model we generated along the same N-trending seismic profile as shown in Figure 7. Our interpretation supports the presence of a WGR-cored antiform in the southern and northern Måløy Slope (Fig. 8). Note that the model requires the presence of a northern, WGR-cored structural high, although it is not clearly reflected in the amplitude of the gravity signal or in the gridded Bouguer gravity anomaly map (anomaly A5 in Fig. 6A). The absence of a strong gravity anomaly is explained by the similar density values of adjacent basement blocks in the area, which do not cause lateral density contrasts large enough to create a pronounced anomaly. Furthermore, the required lower-density values for the WGR in the north compared to the south are in agreement with an overall northward decrease in average rock densities within the WGR onshore (cf. Olesen et al., 2010).

To explain the high-amplitude, long-wavelength magnetic anomaly in the north (Fig. 8), we had to introduce two highly magnetic blocks to the model at deeper basement levels. We infer that the smaller of these two blocks represents an intrusive body forming part of the series of at least four positive magnetic anomalies that are observed further north along the Norwegian coastline (Olesen et al., 2010). This interpretation is supported by boreholes 6306/10–1 and 6305/12–2, which are located above one of those magnetic anomalies, and which encountered intrusive, granitic rocks within basement (for details, see Data Repository Item DR6). However, we could not observe a distinct seismic facies (e.g., high-amplitude, crosscutting or discordant seismic facies) that could be readily attributed to an intrusive igneous body. The lack of reflectivity may be explained by minimal density and P-wave velocity contrasts and therefore reflection coefficients that may be too small to generate reflections distinguishable from noise. Moreover, gradational or irregular contacts between the intrusion and the country rock, or a largely homogeneous intrusion, may cause poor imaging of the postulated intrusion (Johnson, 1990; Smithson and Johnson, 2003). We interpret the second highly magnetic block (2.8 g/cm3/0.013 SI) to reflect a part of the WGR that is more strongly magnetized than its surroundings, and that may therefore comprise a more mafic portion of the WGR.

The high degree of lithological heterogeneity within the allochthons suggests that the density distribution within this unit is more complex than in the WGR (Table 2). To capture this heterogeneity in our 2-D model, we assumed that the juxtaposition of different allochthonous rock types is reflected in the observed large-scale changes in the character and geometry of seismic reflection patterns within the interpreted allochthons. We therefore used boundaries between different reflection patterns to divide the allochthons into several blocks, as indicated by the yellow dashed lines in Figure 7C. Despite introducing more uncertainty to the model with this procedure, dividing the allochthons into several blocks with distinct density and susceptibility allowed us to generate a better fit between the observed and modeled gravity and magnetic anomalies (Data Repository Item DR3).

A comparatively dense block (2.87 g/cm3) in the central part of the profile correlates with the ∼10-km-wide zone of steep S- to SW-dipping reflections of SFRc (Fig. 5C), which we interpret as a zone of intense deformation. If our interpretation is correct, the increased density value may reflect high-density rock types such as anorthosites, amphibolites, gabbros, and/or greenschist within a shear zone (Table 2; see NGU onshore petrophysical database, http://aps.ngu.no/kart/geofysikk/). These rock types are observed within the strongly mylonitized NSDZ in the area of Florø and Eikefjord onshore (Fig. 1A; Bryhni and Lutro, 2000a, 2000b; Lutro and Bryhni, 2000). Alternatively, greenschists or gabbros associated with ophiolites of the Upper Allochthon could explain the elevated density, but not the low magnetic susceptibility of the block. Unfortunately, the model could not be used to either confirm or disprove the presence of Devonian sediments on the Måløy Slope, as average density and magnetic properties of Devonian sediments and allochthons appear too similar to have a noticeable effect on the observed gravity and magnetic anomalies (Table 2).

Despite the ambiguity inherent in the interpretation of potential field data, the model we present here provides one geologically plausible solution that is broadly consistent with the seismic facies variability observed within reflection data, and the geometry of basement structures confidently mapped just onshore. Uncertainties remain with respect to the range of density and susceptibility values that can occur within each rock type (Table 2), the exact depth and geometry of source bodies, and the effects of potential lower-crustal heterogeneity.

DISCUSSION

Link Between Intrabasement Seismic Reflectivity and Basement Structural Architecture and Composition

We have shown that different basement units (WGR, allochthons, and NSDZ) can be distinguished by their seismic characteristics and physical properties (Table 2; Fig. 8). We here discuss briefly the factors that may contribute to these different seismic characters.

SF1

The simplest explanation for subparallel reflections is compositional layering, which could be caused by tens-to-hundred-meter-thick coarsening-upward cycles within Devonian metasediments, comparable to those observed onshore (Steel, 1976; Steel et al., 1977; Steel and Aasheim, 1977). However, borehole data, and gravity and magnetic forward modeling do not confirm the presence of Devonian metasediments on the Måløy Slope. Alternatively, compositional layering may occur within the allochthons. Onshore geological data show that the allochthons are composed of a variety of accreted terranes and stacked felsic and mafic rock types of varying metamorphic grade (e.g., Andersen et al., 1994; Gee et al., 2008; Hacker et al., 2010); as a result, compositional layering can at least partly explain the observed basement reflectivity.

Since SF1 is exclusively observed in the uppermost part of the crystalline basement (up to 2 km below top acoustic basement), we also consider different degrees of weathering and fluid-filled fracture networks as potential causes for subparallel reflection patterns. Our borehole analysis confirmed that the uppermost part of the crystalline basement is likely weathered and/or contains fluid-filled fractures. The presence of fluids can decrease the P-wave velocity of a rock unit by up to 20% (e.g., Goodwin and Thompson, 1988; Johnson and Hartman, 1991; Smithson et al., 2000; Healy et al., 2009), which can lead to strong impedance contrasts, and therefore, high-amplitude reflections that may be stronger than reflections from primary lithological changes (Hyndman and Shearer, 1989; Harjes et al., 1997). Fluid-filled fractures can cause scattering and crosscutting horizontal reflections (Ganchin et al., 1998). Fractures and joints of all scales and varying orientations are furthermore a common feature of the rocks exposed onshore (Fossen et al., 2016). These observations suggest that fluid-filled fractures may play a major role in the reflectivity of the upper parts of the crystalline basement, provided that such fracture systems have a strong horizontal component and thus are imaged by seismic reflection data.

SF2

We interpret the allochthonous basement units to be characterized by often folded, medium- to high-amplitude, laterally discontinuous reflections within otherwise acoustically transparent or chaotic seismic facies (SF2, Fig. 4). Folded reflections are assumed to represent either large-scale migmatitic structures (Mehnert, 1968) or folded layers of different lithologies (Figs. 4B and 7). In response to the extensive Devonian deformation, gneisses are expected to have formed migmatites within the basement, a common metamorphic rock type characterized by alternating millimeter- to kilometer-thick, felsic and mafic layers; such pronounced lithological heterogeneity may cause prominent intracrystalline basement reflections (e.g., Kern et al., 2002; Smithson and Johnson, 2003). Based on onshore observations (e.g., Bryhni and Lutro, 2000a, 2000b; Andersen et al., 1994; Osmundsen and Andersen, 2001; Johnston et al., 2007), subseismic, upright and recumbent isoclinal folds are expected to further contribute to the seismic wave field, forming scatterers or packages of discontinuous reflections, respectively (Smithson and Johnson, 2003).

SF3

High-amplitude, relatively continuous reflections and reflection packets of SF3 are interpreted as mylonitic shear bands; these formed within the NSDZ during Devonian extension (Fig. 4C). Shear zones often form pronounced continuous reflections and reflection packets due to their planar geometry and layering, and high densities and P-wave velocities. However, shear zones can locally also have the characteristics of a diffractor because of their complex internal architecture (folding, lensing, boudinage; e.g., Smythe et al., 1982; Fountain et al., 1984; Hurich et al., 1985; Christensen and Szymanski, 1988; Wang et al., 1989; Kern and Wenk, 1990; McDonough and Fountain, 1988; de Wit et al., 2001). Therefore, abrupt terminations of high-amplitude reflections within SF3 are likely to be caused by lateral changes in composition, strain, orientation to the seismic wave front (Brocher and Christensen, 1990; McDonough and Fountain, 1988), and/or structural complexity of the mylonitic rock body itself. The mineralogical composition of the shear zone remains speculative, but it is likely to consist of several, strongly sheared rock types. Onshore west Norway, different rock types within the NSDZ are classified as “mylonites,” with compositions ranging from quartzites to mica schists, anorthosites, amphibolites, greenschists, and gneisses (Bryhni and Lutro, 2000a, 2000b; Lutro and Bryhni, 2000). These rock types are also likely to be present within the NSDZ offshore, adding to the complexity of the observed seismic reflectivity.

SF4

Reflections classified as SF4, which we infer belong to WGR material, are most clearly observed at depths >12 km (Figs. 4D and 7). The commonly short, apparently crosscutting reflections are suspected to represent large-scale conjugate or doubly vergent S-C fabrics similar to reflection patterns described from mid- and lower-crustal seismic reflection sections elsewhere (Reston, 1988; Varsek and Cook, 1994; Torvela et al., 2013). Furthermore, field studies show that such deformation fabrics exhibit a remarkable self-similarity at a wide range of scales (10−8 to 105 m; Hippertt, 1999; Waters et al., 2003), making their seismic imaging plausible. Conjugate S-C fabrics form under coaxial deformation during the late stages of shear zone evolution (Fossen, 2010) and are commonly associated with mylonitic rocks. Therefore, reflections of SF4 can most likely be attributed to the internal mineral and velocity anisotropy and structural complexity of sheared rocks, rather than lithological variations within the WGR. Angles between S and C planes in SF4 vary between 20° and 40°, which are within the range of typically observed values for such fabrics (Fossen, 2010). Another explanation for the reflection geometry observed for SF4 is anastomosing shear zones characterized by individual bands that exhibit different directions of motion. Both conjugate and anastomosing shear zones are defined by higher-strain zones that encase lower-strain (protolithic) lenses (e.g., Hudleston, 1999); these two types of shear zones may become indistinguishable from one another at very high strains as shear bands rotate into parallelism with the shear direction.

Both interpretations indicate that the upper part of the WGR is mylonitic, an interpretation that is consistent with onshore observations. Based on data from the Kvamshesten area, Andersen and Jamtveit (1990) reported a transition in deformation fabrics within the WGR from penetrative simple shear fabrics close to the NSDZ to inhomogeneous pure shear fabrics that are truncated by anastomosing zones of simple shear 2–6 km below the NSDZ. This transition in deformation fabrics, and decrease in strain intensity with depth, can be correlated with the transition from SF3 (high-strain zone; planar mylonitic rocks) to SF4 (medium-strain zone; anastomosing mylonitic shear bands/S-C fabrics). Hence, we conclude that variations of the intrabasement facies and reflection patterns are not only a result of lithological changes, interference, and scattering effects, but they can also be linked to the style and amount of ductile deformation that took place in this and likely other areas (e.g., Klemperer, 1987).

Onshore-Offshore Correlation of Crystalline Basement Units

Through 3-D mapping of the four, laterally correlatable seismic facies (Fig. 4) and three key seismic facies relationships (Fig. 5), and by integrating our seismic interpretation with the onshore geology, borehole data, and 2-D forward modeling of gravity and magnetic data, we constructed a 3-D geometric model of the basement architecture offshore west Norway, relating it to the onshore geology presented by Osmundsen and Andersen (2011; Fig. 9). The geometry of the NSDZ is represented by the physiography of the WGR (Fig. 9C), which correlates to WGR antiforms and synforms onshore (Fig. 9A). WGR synforms only comprise allochthonous basement rocks, since we were unable to find convincing evidence for the presence of Devonian metasediments offshore. The absence of Devonian metasediments on the Måløy Slope suggests either nondeposition or erosion of the sediments during Devonian extension and exhumation. We favor the latter explanation, since very thick (6–7 km) Devonian metasediments are identified within large half grabens further west, on the East Shetland Platform (Platt, 1995; Platt and Cartwright, 1998; Marshall and Hewett, 2003), and in the Orcadian and West Orkney Basins (Norton et al., 1987; Coward, 1990: Wilson et al., 2010), suggesting Devonian sediments were present over a much larger area in the past.

Figure 10 shows a 1:1 scale, E-W cross section illustrating the geometric relationships among the WGR, NSDZ, and Caledonian allochthons on the Måløy Slope as interpreted from seismic data, in comparison to the onshore geology near the Devonian Hornelen Basin. The onshore part of the cross section was derived from literature and geological maps (based on Norton, 1986, 1987; Bryhni and Lutro, 2000a, 2000b; Lutro and Bryhni, 2000; Fossen et al., 2016), although the exact structural configuration below the Hornelen detachment is uncertain (indicated by question marks in Fig. 10). The Hornelen detachment defines the contact between the Hornelen Basin and the underlying mylonitic NSDZ (Fossen et al., 2016). Both the lower part of the allochthons and the upper part of the WGR are incorporated in the NSDZ (Fig. 10). Pink color gradients illustrate the upward and downward decreases in strain in the allochthons and WGR away from the core of the NSDZ, respectively. Geological maps (e.g., Fig. 1A) and published onshore profiles (e.g., Fig. 1B) indicate the presence of the Lower and Middle Allochthons within the NSDZ south, east, and below the Hornelen Basin. The exact thickness of the allochthonous units onshore remains unknown (maximum 5 km; Osmundsen and Andersen, 2001).

Offshore, we interpret a slightly steeper dipping NSDZ, represented by an ∼2–4-km-thick, highly reflective (SF3) mylonitic band, separating allochthonous rocks above from the WGR below (Figs. 2 and 10). Here, the allochthons are up to 5 km thick within WGR synforms. The ∼23-km-wide coastal area between the Måløy Slope and Hornelen Basin (labeled “Based on magnetic data” in Fig. 10) likely has allochthonous material in the upper part of the crystalline basement (Lutro and Bryhni, 2000). Using gravity and magnetic data, Smethurst (2000) interpreted this area as “Caledonian ophiolite” (= Upper Allochthon, density: 2.85–2.95 g/cm3; susceptibility: 0.0003–0.001; Olesen et al., 2010). Rocks associated with the Upper Allochthon (greenschists and amphibolites) have a depositional contact with the Devonian Hornelen Basin (Lutro and Bryhni, 2000; Fossen et al., 2016), with Upper Allochthon rocks also present in the nearby Eikefjord and Bremangerlandet areas (Fig. 1A). It therefore seems plausible that crystalline basement along the coastal area offshore Hornelen consists of ophiolitic rocks of the Upper Allochthon.

Our interpretation of the depth and architecture of the WGR and NSDZ shown in Figure 10 has important implications for the location of the suture zone between Baltica and Laurentia. Due to a W-NW–directed increase in metamorphic grade of the WGR onshore, several authors have suggested that the suture zone between Baltica and Laurentia is located beneath the northern North Sea, and that it may have a similar, broadly northerly trend as the North Viking and Sogn grabens (e.g., Fossen, 2000; Lyngsie and Thybo, 2007; Fossen et al., 2016). For similar reasons, Cocks et al. (2005) interpreted the suture zone to follow a NNE-NE trend in the Norwegian Sea. Based on our mapping and interpretation of seismic facies, the NSDZ and WGR (i.e., the eastern edge of Baltica) appear to continue into the Sogn graben (Fig. 10). Considering that reconstructed Caledonian paleopressures from the eastern margin of the Hornelen Basin and Nordfjord area are already very high (20–25 kbar = 65–90 km depth; Wain, 1997; Krabbendam and Dewey, 1998), and assuming a constant dip of the NSDZ of 15°, the WGR must have been subducted to extreme depths of 105–130 km underneath the present-day Sogn graben, ∼150 km further east of the Hornelen Basin. However, due to the likelihood of a varying dip angle of the NSDZ and a poor understanding of the postkinematic rotation history of the NSDZ (Norton, 1986, 1987; Hacker, 2003), these paleodepth values are likely overestimated. Alternatively, the high-amplitude reflections of SF3 that reach the Sogn graben (Figs. 2 and 10) may not belong to the NSDZ, but to a different, so-far-unidentified Devonian shear zone (cf. Fazlikhani et al., 2017).

CONCLUSIONS

Using conventional 2-D and 3-D seismic reflection, borehole, core, and petrophysical data, and 2-D gravity and magnetic forward modeling, as well as onshore-offshore correlations, this study resolved the 3-D geometry and composition of crystalline basement structures offshore west Norway. We have shown the following:

  • (1) The crystalline basement along the Måløy Slope is characterized by a pronounced compositional and structural heterogeneity. Offshore well-log and core data indicate the presence of different kinds of gneisses that are distinguished by their mineralogy and internal deformation fabrics.

  • (2) The basement is characterized by a high degree of intrabasement reflectivity that can be divided into four distinct, laterally correlatable seismic facies. Based on our onshore-offshore correlations and well-log and seismic analysis, we suggest that compositional layering and scattering effects caused by subseismic heterogeneity are likely the main causes for the observed intrabasement reflectivity. Furthermore, mylonitic shear zones, zones of intense fracturing, and the limbs of isoclinal and recumbent folds are considered to cause significant reflectivity.

  • (3) The results of seismic data interpretation show that the basement reflectivity is created by the superimposition of the main basement units in the area, i.e., WGR, NSDZ, Caledonian allochthons, and possibly Devonian sediments. Each unit is distinguished by its seismic character, physical properties, 3-D geometry, stratigraphic position, and spatial relationship to other units.

  • (4) Gridded gravity and magnetic data helped us to constrain the location and geometry of large-scale intrabasement synforms and antiforms, to ascertain lateral changes in the physical properties of the basement, and to bridge the gap between offshore 3-D seismic data and the onshore geology.

  • (5) The interpreted large-scale structural architecture of the crystalline basement offshore west Norway correlates well with onshore structural models (e.g., Osmundsen and Andersen, 2001) and is further constrained by 2-D gravity and magnetic forward modeling. In contrast to previous studies, use of the seismic interpretation of basement structures as geometric constraints for the 2-D gravity and magnetic model is found to be a more robust approach for modeling the distribution of density and magnetic properties in the subsurface.

  • (6) The basement rock types and structures described in this study are common within deformed continental crust and may thus occur elsewhere. Hence, we consider that the approach employed here can be used in other subsurface studies and that our results are a valuable input for analogue and numerical modeling studies to more realistically model processes such as rift development above heterogeneous substrate.

ACKNOWLEDGMENTS

We would like to thank Equinor ASA for allowing us to show the seismic and borehole data used in this study, along with Schlumberger, Geosoft, and Senergy for providing academic licenses for use of Petrel, Oasis montaj, and Interactive Petrophysics software, respectively. In addition, we would like to thank the Norwegian Geological Survey for the permission to show the gravity and magnetic data. We thank the editor, Tony Doré, and two anonymous reviewers for their thorough and thought-provoking reviews, which substantially improved our initial submission. This study forms part of the MultiRift Project, funded by the Research Council of Norway’s PETROMAKS program (project number 215591), and Equinor ASA to the University of Bergen and partners Imperial College, University of Manchester, and University of Oslo.

1GSA Data Repository Item 2019086, Item DR1: Details about the noise-filtering process applied to seismic data; Item DR2: 3-D velocity model and depth conversion of seismic data; Item DR3: Details about the 2-D forward modeling of gravity and magnetic data and figures illustrating several steps of the modelling process; Item DR4: Uncertainties of 2-D forward modelling process; Item DR5: Large-scale version of seismic profile shown in Figure 7; Item DR6: Additional core photos of crystalline basement rocks and locations of cores relative to magnetic anomalies, is available at http://www.geosociety.org/datarepository/2019, or on request from editing@geosociety.org.
Gold Open Access: This paper is published under the terms of the CC-BY-NC license.