Pre-peer review manuscript

A reevaluation of the proposed ‘Lairg ’ and its potential implications for the deep structure of northern

Michael J. Simms1 and Kord Ernstson2

1. Department of Natural Sciences, National Museums Northern Ireland, Cultra, Holywood, Co. Down BT18 0EU, Northern Ireland. [email protected] 2. Faculty of Philosophy, University of Würzburg, Germany. [email protected]

Abstract

It has been suggested that the Lairg Gravity Low represents a buried ~40 km across from which the 1.2 Ga Stac Fada Member ejecta deposit originated. The structure is too large to represent a simple crater, and there is no evidence of a central peak. Reanalysis of the point Bouguer gravity anomaly data reveals a ring of positive anomalies around the central negative anomaly and we interpret it as the eroded central part of a peak ring crater. Peak ring craters show a consistent 2:1 relationship between peak ring diameter and total crater diameter, implying that the putative Lairg crater may be ~100 km across. This would place the rim of the crater within a few km of the Stac Fada Member , a proximity that is inconsistent with the thickness and clast size of the ejecta deposit. We propose that the impact crater was originally formed further east, at a substantially greater distance from the Stac Fada Member than today. Subsequently it was translocated westwards, to its present location beneath Lairg, during the Caledonian Orogeny. This suggests that a deep-seated , analogous to the Flannan and Outer Isles thrusts, exists beneath the Moine Thrust in north-central Scotland.

Introduction

The Stac Fada Member is a 4-12m thick unit exposed intermittently along the ~50km outcrop of the fluviatile and lacustrine Stoer Group (Mesoproterozoic, ~1.2 Ga) in north-west Scotland. It contains abundant devitrified angular melt clasts in a matrix and, for decades, was thought to be volcanic in origin (Stewart 2002). However, the discovery of shocked mineral grains (Amor et al. 2008, Osinski et al. 2011, Reddy et al. 2015) prove that it is the product of a km-scale impact. The near continuity of the Stac Fada Member outcrop was interpreted as evidence of its relative proximity to the impact crater (Amor et al. 2008) yet Stoer Group strata beneath the impact deposit are undisturbed, indicating that the present outcrop lies significantly beyond the crater rim. 1

Amor et al. (2008) suggested that the crater was located to the west, perhaps buried beneath a thick cover of younger rocks in what is now Basin, but Stewart (2002) had contended previously that the source of the Stac Fada Member, which he considered volcaniclastic in origin, actually lay to the east. Simms (2015) similarly argued for an eastern source leading to the inference that the crater, if it still exists, may lie beneath mainland Scotland. However, Proterozoic and Archean target rocks across much of the region now lie buried beneath a cover of predominantly Moinian (Neoproterozoic, ~1 Ga) metasediments that were thrust westwards across northern Scotland ~430 Ma ago (McClay and Coward 1981). As such it is unlikely that any physical manifestation of the crater will exist at the surface, even assuming that it survived many millions of years of post-impact erosion and the effects of subsequent tectonism associated with the Moine Thrust. However, geophysical methods offer the potential for locating a buried impact crater (Pilkington and Grieve 1992).

Geophysical surveys (BGS map; Rollin et al. 2009; Leslie et al. 2010) have revealed a deep gravity anomaly centred on the town of Lairg, little more than ~50 km east of the Stac Fada Member outcrop at its closest point. There is a remarkable correspondence between the location of the Lairg Gravity Low, as it is known, and the location of the impact crater as predicted from inferred source directions of the Stac Fada Member impact ejecta sheet (Simms 2015). Comparing the Lairg Gravity Low with gravity signatures of other impact structures led to the suggestion that it might represent an impact structure in the Archean basement that now lies buried by overthrust Moinian metasediments (Simms 2015). The Moine Thrust which underlies these metasediments is considered to be near horizontal, and at relatively shallow depth (~1 km) across much of northern Scotland.

Simms (2015) estimated the putative Lairg impact crater to have a diameter of at least 40 km based on the original geophysical analyses of Rollin et al. (2009), although recognising that erosion and/or tectonic effects might significantly have modified an originally larger structure (Simms 2016). This would place it among the fifteen largest of almost 200 impact craters currently known on Earth (Hergarten and Kenkmann 2015).

Erosion of the impact crater

A pronounced angular exists between the Stoer Group (Mesoproterozoic, ~1.2 Ga) and the Group (Neoproterozoic, ~1 Ga) on the west coast of northern Scotland. Further east the Stoer Group is absent and the Torridon Group rests directly on a deeply eroded surface of Lewisian (Archean, ~3 Ga). These observations testify to the scale of erosion that the crater may have experienced in the almost 200 million year interval between the impact and deposition of the Torridon Group. This is comparable with the time elapsed since the Manicouagan impact structure in eastern

2

Canada was formed 214 million years ago. It is estimated that 2 km or more of post-impact denudation has occurred in this region since the late Triassic (Degeai and Peulvast 2006), reducing the crater from an original estimated diameter of 100 km to the 72 km diameter structure currently visible. As such the Manicouagan Crater might be considered broadly analogous with the putative Lairg crater at the onset of Torridon Group deposition. Post-impact erosion and/or the effects of Caledonian tectonics might have removed shallower parts of the Lairg impact crater, similarly reducing its apparent diameter to what we see today in the published geophysical data. Although nothing of the Lairg structure is visible at the surface, reanalysis of the gravity data may help to clarify aspects of the Archean crust beneath the Lairg Gravity Low and ascertain if what we see today is a true reflection of the original structure.

Crater scale vs. structure

If pre-Torridon erosion removed the outer parts of the Lairg Crater, how might we determine if what remains is a reasonable representation of the original crater or is merely part of a once larger structure? The answer lies in identifying some of the fundamental differences in crater structure that accompany an increase in crater size, from simple bowl-shaped craters to complex multi-ring structures (Morgan et al. 2016), and ascertaining if any of these critical features are evident in the gravity data.

Bowl-shaped Simple Craters on Earth, such as in Arizona, USA, are no more than a few kilometres across and approximate more closely to the shape and dimensions of the original transient crater than do larger impact structures (Melosh 1989). In larger and deeper structures the walls of the transient crater collapse and, coupled with uplift of the central part of the floor, generate Complex Craters that are substantially wider and shallower than the original transient cavity. Failure of the walls along concentric fractures may produce ring-shaped troughs and terraces in the upper part of the crater while accommodation factors associated with inward slumping commonly give rise to radial transpressional ridges (Kenckmann and van Dalwigk 2000). Impact structures become increasingly complex with size and at diameters greater than ~25 km the central peak is replaced by a basin surrounded by a raised ring (Osinski and Pierazzo 2013). The very largest terrestrial structures have multiple rings but on Earth only Sudbury, Chicxulub and Vredefort fall into this category (French 1998). It is these observed changes in crater structure with scale that have implications for reevaluating the putative impact crater beneath Lairg.

The Lairg Gravity Low appears to represent a bowl-like structure, yet it is far too large to be a Simple Crater. The extent and thickness of the Stac Fada Member impact deposit at outcrop also implies a structure at least several tens of kilometres across, which is broadly consistent with the 40 km diameter estimated for the putative Lairg crater (Simms 2015). Such a crater might be anticipated to have either

3

a central peak or a peak ring, yet neither can be discerned in the original analyses of the Lairg Gravity Low by either Rollin (2009) or Leslie et al. (2010). Our reanalysis of the gravity data aims to shed light on the some of the finer details of the structure responsible for the Lairg gravity anomaly. This need not detract from the basic hypothesis that the Lairg Gravity Low represents an entire impact crater since the coarseness of previous analyses may have masked critical diagnostic features or, alternatively, that it actually represents the central basin of a substantially larger peak ring impact structure.

Re-evaluation of the Lairg gravity field

The original gravity dataset analysed by Rollin et al (2009) is freely available from the British Geological Survey. It is this data that has been reanalysed by one of us (KE).

1. Data and new processing

From the Bouguer gravity anomaly database of the British Geological Survey (Rollin 2009) a rectangular field of gravity stations was selected (Fig. 2), located with respect to the negative anomaly centred on Lairg (Fig. 3). The size of the rectangle is somewhat arbitrary and is constrained by the presence of sea in the north-western and south-eastern corners of the region, and by the distribution of the surrounding regional gravity anomalies that may influence the Lairg local anomaly. The window comprises about 10,000 data, corresponding to approximately one gravity station per square kilometre.

For the Bouguer gravity field (Fig. 3) a regional trend field was computed applying strong low-pass filtering of the data by a multifold moving average. We recognise that constructing a gravity regional field is an ambiguous process that may preferentially accentuate long-wavelength or short-wavelength anomalies but, after several attempts with various filters, the regional field shown in Fig. 4 is considered a reasonable compromise. A Bouguer residual anomaly (Fig. 5) was then derived by subtracting the regional field (Fig. 4) from the measured Bouguer field (Fig. 3). This generated a complex picture of short-wavelength anomalies within which a ring of relatively positive anomalies surrounds the conspicuous negative anomaly centred on Lairg.

2. Gravity profiles and model calculations

Three profiles of Bouguer data were constructed from the map of the Lairg residual anomalies, shown in fig. 5, with each crossing the centre of the Lairg negative anomaly (Fig. 6). Their orientation was chosen to evaluate various features of the Lairg anomaly, with the three gravity profiles compared in Fig. 7. All three profiles have a similar shape with a somewhat structured slope towards the central gravity minimum. The peripheral ring of positive anomalies is striking and is broken only by a gap in

4

the north. The NW - SE and SW - NE gravity profiles were selected for a very simple 2.5D model calculation (Fig. 8). In the absence of more specific density data, a straightforward modelling produces a two-layer density distribution that assumes a mass deficit responsible for the Lairg gravity anomaly. Because of this simple assumption the shape of the negative mass follows more or less the shape of the gravity curve and reveals a step-like slope of the central depression with a depth extending to about 3 km below the thrust plane. This is surrounded by a rim wall with a diameter of approximately 50 km, and a peripheral shallow mass deficit.

From potential theory it is known that the integration over a gravity anomaly provides a measure of the total mass responsible for the anomaly, regardless of the density distribution within the mass. Once a model adaptation to a measured gravity anomaly has been performed, its mass is the same for all other density distributions providing they fit the measured gravity. This is the basis for a very coarse evaluation of the central mass deficit producing the Lairg central negative anomaly. Because of the simple mass calculation for a spherical segment the geological mass of the Lairg central structure has been replaced by such a segment with a density of 0.15 g/cm3 (Fig. 9) and a total mass defecit of 7 x 1014 kg.

3. The Lairg Gravity Low as an impact structure

Gravity surveys across terrestrial impact craters have contributed to understanding their specific structure and impact cratering processes in general (Pilkington and Grieve 1992, Sharpton et al. 1993, Pohl 2015, and others).

The Lairg Gravity Low lacks any evidence of a central uplift but our reanalysis of the gravity data does reveal a ring of relatively positive anomalies surrounding the central negative anomaly. We interpret this circle of positive anomalies as an Inner Peak Ring, developed through elastic rebound and slumping of the transient crater walls at the end of the excavation stage (Melosh 1989, Kenkmann et al. 2013). Rock from greater depths may be uplifted in these ring structures which, because of the general increase in density with depth, can account for the positive gravity anomalies. It is not uncommon for the inner ring of a peak ring crater to be incomplete (Morgan 2016) and hence the presence of gaps in the ring-like structure that we have identified does not detract from our interpretation of it as an impact crater. On the contrary, the existence of a ring of negative gravity anomalies beyond the ring of positive anomalies (Fig. 7) supports our model of a peak ring crater. To interpret the ring of positive anomalies as representing the rim of a smaller 50 km-diameter impact structure is inconsistent with density reduction that might be expected in the rim region through deposition of low-density ejecta and slumping in the modification stage.

5

We have compared the gravity profiles across the Lairg anomaly with gravity profiles for several proven impact structures of diameters ranging between 25 km (Ries crater) and ~180 km (Chicxulub impact structure) (Fig. 10). Craters comparable in size to that postulated for Lairg are Manicouagan (~100 km), Popigai (100 km) and Chesapeake (85 km), but there is not a standard gravity profile for these large complex impact structures. This is underlined by the statement of Wünnemann et al. (2005) that pre-impact target properties may exert a considerable influence on the structure of complex craters on Earth. This is reflected in considerable variations in gravity signature even among structures of similar size, such as the Manicouagan and Popigai impact structures (Fig. 11). From this it is evident that the Lairg gravity anomaly has most in common with the 100-km diameter Popigai impact structure but also bears similarities to smaller complex craters, particularly the Rochechouart impact structure in south-west France. The Rochechouart crater has experienced significant erosion such that clear morphological features of a crater structure are no longer evident at the surface. Nonetheless there is clear lithological evidence at the surface for its impact origin, in the form of impact melt rocks, , and dikes. Originally thought to have a diameter of between 20 and 35 km (e.g., Pohl 2015), others have considered a diameter of 40 or 50 km to be more realistic (Sapers et al. 2014). This larger estimate corresponds more closely to the pronounced gravity signature shown in Fig. 10, although not actually referred to by Sapers et al. (2014). In less eroded impact structures in sedimentary or mixed targets, peak ring uplift may be recognized through lithological contrasts, but this has proven difficult for the deeply eroded Rochechouart structure where the target is entirely crystalline. Hence, although the Rochechouart gravity signature suggest a peak-ring uplift there is, as yet, no direct lithological evidence to confirm this (Fig. 10). At the Ries crater, in southern Germany, low-density, post-impact, lake sediments up to 400 metres thick occur within the central part of the crater and contribute significantly to the magnitude of the negative gravity anomaly. The presence of an inner peak ring at Ries is proven by inliers of uplifted crystalline basement rock projecting through the post-impact lake sediments, and also indicated by various geophysical measurements (e.g. Pohl et al. 1977). After subtracting the lake sediment contribution from the overall Ries gravity anomaly the resultant gravity curve closely resembles the Rochechouart and Lairg anomalies. The results from our reanalysis of gravity data across the Lairg region compare favourably with gravity signatures across proven impact structures. Although it does not prove conclusively that an impact crater exists at this location, it does lend support to our hypothesis that the Lairg Gravity Low may represent part of a once much larger impact structure.

4. Extrapolating crater size from peak-ring diameter

6

A consistent relationship has been observed for the relative dimensions of peak ring craters, with total crater diameter approximately twice that of peak ring diameter (Fig 11). Although there are relatively few well-documented examples on Earth, many others are known from rocky bodies elsewhere in the Solar System and this relationship appears to be independent of both crater size and location on other Solar System bodies (Melosh 2015). This ‘Factor 2’ rule for complex impact craters with inner peak rings has significant implications for interpreting the gravity signature of the Lairg structure. If our estimate of 50 km diameter for the Lairg Gravity Low represents the full crater diameter, then we might anticipate a peak ring 25 km across. By comparison with other impact crater gravity signatures, previously discussed, then we would expect to recognize this in the Lairg gravity data as at least a subtle inflexion or change of gravity gradientsuggesting a diameter of roughly 100 km for the complete Lairg impact crater.

Applying this estimated figure to the approximate mass deficit calculated from the gravity anomaly indicates that the Lairg structure is comparable with the mass deficit associated with many other terrestrial impact structures (Fig. 12) and, as such, this would make Lairg the largest impact structure yet discovered in Europe.

Constraints on the size of the Lairg impact structure

The Stac Fada Member impact deposit was described as a 'proximal ' on account of its near continuous outcrop (notwithstanding recent erosion) (Amor et al. 2008). The absence of large scale disruption of strata or significant soft-sediment deformation in the pre-impact Stoer Group succession indicates that the present outcrop, which at its closest point is just 55 km from the centre of the Lairg Gravity Low, lay beyond the area materially affected by seismic shocks associated with the impact and ensuing crater formation. Claims have been made for impact-induced seismites associated with, but hundreds of kilometres from, several large impact craters, among them Chicxulub (Renne et al. 2018), Manicouagan (Clutson et al. 2018) and Sudbury Addison et al. 2005). For many other craters no such claims have been made, although this may reflect an absence of appropriate strata in which they might be preserved. However, the absence of any impact-induced seismite in the Stoer Group, among sediments that demonstrably are prone to soft sediment deformation (Stewart 2002) remains enigmatic and suggests that the impact was a considerable distance from the present outcrop of the impact deposit.

But what other constraints might be placed on the original extent of the Lairg impact structure? In particular, how might a potentially 100 km diameter impact crater be accommodated in the gap between Lairg, at the approximate centre of the structure, and the present outcrop of the Stac Fada

7

Member on the coast to the west? A crater of 50 km diameter centred on Lairg would not impinge on the main Lewisian outcrop west of the Moine Thrust, which is consistent with the absence of any impact related structures there. However, this apparent absence of evidence need not necessarily discount the possibility that the crater was originally twice this size, at ~100 km across, extending significantly beyond the present limit of the Moine Thrust.

The scale of post-impact erosion experienced by the crater can be inferred from the angular unconformity between the Stoer Group, dipping at 23o to the west, and the near horizontal Torridon Group that truncates it. This may be significant for delimiting the original extent of the crater. Hundreds of metres, and perhaps even 2 km or more, of Stoer Group strata were removed from areas to the east of its present outcrop where Torridon Group strata now rest directly on the Lewisian basement. No direct evidence of impact, in the form of shatter cones, pseudotachylite, or large-scale brecciation, has been identified in any of the Lewisian between the Moine Thrust and the Stoer Group outcrop, while published geological maps do not reveal any pattern of large-scale fractures that might reflect the existence of an impact crater. However, if we consider the processes of formation it is possible that the outer edge of the crater, once extending into this foreland region, might have been removed by erosion.

During the initial moments of the impact process a deep transient crater is formed, of the order of 8 km for a final crater width of ~50 km. Gravity-driven collapse processes rapidly transform this transient cavity into a wider, and shallower, complex crater in which shallow structures, such as listric faults and subhorizontal basement detachments, develop towards the margins (Kenkmann and Dalwigk 2000). Inevitably these shallow marginal structures are more vulnerable to erosion than the deeper more central parts of the crater. The potential scale of erosion already described may then have removed the shallow outer edges of the crater across the current Lewisian outcrop before deposition of the Torridon Group commenced. The westward dip of the Stoer Group beneath the basal Torridon Group unconformity (Stewart 2002) has further implications for the potential destruction of any westward extension of the crater. It implies that the Lewisian basement, on which the Stoer Group rests, was similarly tilted for an unknown distance east of the present Stoer Group outcrop. By inference, the sub-horizontal basal Torridonian unconformity, passing as it does onto successively older Stoer Group strata and then onto the underlying Lewisian, reflects an eastward increase in the scale of erosion. Assuming that the westward tilt evident in the Stoer Group did not extend across the entire outcrop/subcrop of the Lewisian basement then some sort of 'hinge' structure, or perhaps a series of rotational blocks tilted to the west, must exist to the east of the present Stoer Group outcrop. Were this not the case then simple geometry would necessitate more than 30km depth of pre-Torridonian erosion towards the eastern side of the subcrop. Cambrian strata rest uncomformably upon the

8

Torridon Group and today dip eastwards (Peach et al. 1907). This implies a still steeper westward tilting of the Stoer Group during the early Palaeozoic and, as such, the depth of erosion of a hypothetical ‘unhinged’ basement would be greater still.

Resolving the paradox

It is impossible to accommodate a 50 km radius impact crater between the Stoer Group outcrop and its centre at Lairg solely by invoking post-impact erosion of the outer parts of the crater. Our upper estimate of crater size, approximately 100 km, reduces the distance between the original crater edge and the current Stoer Group outcrop to a minimum of just a few kilometres. This is incompatible with what is, or rather what is not, observed within the Stac Fada Member and contiguous strata. The absence of any impact-generated seismite within the of the immediately pre-impact Stoer Group, sediments that demonstrably are prone to soft-sediment deformation (Stewart 2002), has already been alluded to. Furthermore, the Earth Impact Effects programme (Collins et al. 2005) predicts that an ejecta deposit in such close proximity would be several hundred metres thick and would incorporate very much coarser debris than is found in the Stac Fada Member ejecta deposit. The inevitable conclusion is that if the putative Lairg impact crater is as large as our interpretation of the gravity data suggests, then the spatial relationship between the crater and the Stac Fada Member outcrop must have changed since they formed.

The original interpretion of the Lairg Gravity Low as an impact crater (Simms 2015) followed the consensus view that the Moine Thrust did not descend significantly into the Archean basement. As such the relationship between the crater and the present Stoer Group outcrop would not have changed significantly. With an estimated crater diameter of just 40 km this did not present a major issue, since the closest approach of the Stoer Group outcrop was more than 30 km beyond the supposed rim of the crater. However, to reconcile a 100 km diameter crater centred on Lairg with an impact ejecta deposit only ~10 metres thick and just 60-70 km to the west of Lairg, it would be necessary to invoke crustal shortening of the Archean basement sufficient to translocate the crater westwards for at least several tens of kilometres from its original position. Is this possible and, if so, is there any evidence that might support such a model?

The Moine Thrust has been the focus of geological research for more than a century (Law et al. 2010). Its outcrop has been mapped in extraordinary detail and is the subject of numerous papers interpreting its structure. The current prevailing view is that the Moine Thrust is part of a thin-skinned thrust belt in which movement is located largely above the Archean basement (in which the putative crater is

9

located) and extends east at relatively shallow depth for tens of kilometres from its present surface trace (Coward 1980, 1983, Elliott and Johnson 1980, Butler and Coward 1984, Butler 2010). Interpretation of the deep crustal structure further east in northern Scotland remains far from resolved (Butler 2007) and is based to a significant extent on extrapolation from observations at outcrop to the west and from direct observations of basement inliers and shallow structures within the Moinian further east. However, Brewer and Smythe (1984) specifically state that "surface mapping alone...is incapable of revealing all the fundamental structures of an orogenic belt".

In contrast to the thin-skinned models, others have used essentially the same field observations, often coupled with geophysical data, to argue for a much steeper thrust extending down to the Moho (Watson and Dunning 1979, Stewart 1982, Soper and Barber 1982), with some of these models not far removed from current thinking on the nature of thick-skinned tectonics (e.g. Soper and Barber 1982). Significantly in this respect, geophysical evidence does suggest that a component of thick-skinned Caledonian thrusting may exist beneath northern Scotland.

Several large-scale geophysical surveys of parts of northern Britain, including LISPB (Lithospheric Seismic Profile in Britain; Barton 1992), MOIST (Moine and Outer Isles Seismic Traverse; Blundell et al. 1985) and BIRPS (British Institutions Reflection Profiling Syndicate; Brewer et al. 1983, Brewer and Smythe, 1984), have identified prominent reflectors broadly parallel to the Moine Thrust and descending eastwards at angles between about 20o and 45o. Some flatten out at 17-20 km depth while others penetrate a conspicuous reflector that has been interpreted as the Moho. At least one of these inclined reflectors can be attributed to a known thrust fault, the Outer Isles Thrust, that is developed entirely in Lewisian basement. The others are assumed to represent similar structures. McBride and England (1994) interpreted them as evidence for thick-skinned thrusting associated with the Moine Thrust, which suggests that additional thrusts might exist within the Lewisian basement beneath the Moine Thrust itself. This view is broadly in accordance with that previously suggested by Butler and Coward (1984) who envisaged that thrust imbrication was developed within the Lewisian along ductile shear zones.

The scale of crustal shortening attributable to the Moine Thrust has been a key aspect of many publications, with estimates of between 25 km and 100 km suggested (Elliott and Johnson 1980, Butler and Coward 1984, Barr et al. 1986, Butler et al. 2006). Many of these figures are based on thin- skinned models, in which movement has occurred predominantly in the cover rocks (Moine metasediments and Proterozoic to Lower Palaeozoic sediments) rather than in the Lewisian basement and, as such, they have little relevance to the question we raise here of how the Lairg crater might have been translocated westwards.

10

However, the Crustal Duplex model of Soper and Barber (1982), based on the LISPB data, envisaged a sigmoidal form in which the Moine Thrust steepened eastward to a maximum of 40o to 45o before flattening out at a depth of perhaps 10 km where it joined a floor thrust in the lower crust. This is more akin to current thinking on thick-skinned tectonics than the steep structure envisaged by Watson and Dunning (1979), with the latter model rendered untenable if the substantial displacements measured on the Moine Thrust Zone are accepted. Significantly, Soper and Barber (1982) cited a value of 66% for Caledonian crustal shortening across the northern Highlands. As such a thick-skinned model such as this might provide a mechanism by which the crater, or at least its eroded remnant, could be moved westwards towards the Stoer Group outcrop.

Discussion

Our suggestion that the Lairg Gravity Low represents the central part of a complex peak ring crater has come from a reanalysis of the gravity data. The most parsimonious interpretation might be that the Lairg Gravity Low represents a 50 km crater with the outer rim represented by the ring of positive anomalies surrounding the central negative anomaly, and the peak ring perhaps represented by subtle changes in the gravity gradients within the central low. However, this is consistent neither with the gravity signatures for known peak ring craters nor with the specific gravity signature of the Lairg Gravity Low where a ring of positive anomalies is itself surrounded by slight negative anomalies. Instead the Lairg Gravity Low compares more closely with the peak ring and central low of a substantially larger complex crater. In this respect it is interesting to note that Schedl (2015) used the thickness of the Stac Fada Member and the size of accretionary lapilli within it to arrive at an estimate of 80 to 160 km for crater diameter and a location 225 to 325 km away.

Meteorite impacts on the scale of that suggested by the extent of the Stac Fada Member will cause extensive and deep fracturing of the lithosphere. These fractures potentially may persist as major crustal weaknesses for hundreds of millions of years (Norman 1984). Substantial erosion of the putative 'Lairg Crater', between its formation in the Mesoproterozoic and its burial in the Neoproterozoic, may have reduced it from 100 km across to the 50 km diameter remnant now evident from the gravity data but deep fracturing associated with the impact may have extended significantly beyond this eroded remnant. These fractures may have acted subsequently as foci for deep-seated thrust movements during the Caledonian Orogeny. As such a 100 km diameter Mesoproterozoic impact crater may have been the architect of its own translocation from a site substantially further east to its current resting place beneath Lairg. This process may have further truncated the outer parts of the crater, perhaps accounting for the markedly straight western flank of the Lairg Gravity Low that runs virtually parallel to the strike of the Moine Thrust. 11

The deep structure of the crust to the east of the remains unclear. Present consensus, based on structures seen at outcrop in the west, favours a thin-skinned model. However, steep thrust planes descend deep into the crust to the north and west of the mainland outcrop of the Moine Thrust and demonstrate that thick-skinned tectonics have played an equally significant role in crustal shortening across this region during the Caledonian Orogeny. As such we suggest that movement on a similar deep structure, or structures, further east could accommodate an impact structure on this scale and might account for the seeming mismatch between the inferred scale of the impact, as deduced from the Lairg Gravity Low, and its proximity to the Stac Fada Member impact deposit.

Other types of geophysical survey across the region have proven inconclusive in understanding the nature of the Lairg Gravity Low. There is a minor magnetic anomaly in the same area, but it covers a substantially smaller area than the gravity low and has been attributed to the Rogart-Grudie granite pluton (Rollin 2009). At sites of proven impact structures magnetic surveys generally have proven less useful than gravity surveys. For some, such as Chicxulub, there is a distinct magnetic signature (Rebolledo-Vieyra and Urrutia-Fucugauchi 2004) but for many others the signal is at best equivocal and varies greatly according to target rock composition, impact-related magnetization, and the effects of subsequent crater fill (Cowan and Cooper 2005).

Ultimately answers to these various questions, concerning the scale and even the very existence of the proposed ‘Lairg impact structure’, and the postulated existence of deep-seated thrust faults passing beneath it, must await a programme of deep drilling. Even if ultimately the ‘Lairg impact hypothesis’ proves to be unfounded, such an investigation could nonetheless throw considerable light on the structure of northern Scotland east of the Moine Thrust Belt

Acknowledgements

We thank Renegade Pictures and Lairg & District Community Initiatives for their support of fieldwork associated with this project.

12

References

Addison, W.D., Brumpton, G.R.,Vallini, D.A., McNaughton, N.J., Davis, D.W., Kissin, S.A., Fralick, P.W. & Hammond, A.L. 2005. Discovery of distal ejecta from the 1850 Ma Sudbury . , 33, 193-196.

Amor,K., Hesselbo, S.P., Porcelli, D., Thackrey, S. and Parnell, J. 2008. A proximal ejecta blanket from Scotland. Geology, 36, 303-306.

Barr, D., Holdsworth, R. E. & Roberts. A. M. 1986. Caledonian ductile thrusting in a Precambrian metamorphic complex: the Moine of NW Scotland. Bulletin of the Geological Society of America, 97, 754-764.

Barton, P. J. 1992. LISPB revisited: a new look under the Caledonides of northern Britain. Geophysical Journal International, 110: 371–391. doi:10.1111/j.1365-246X.1992.tb00881.x

Blundell, D.J., Hurich, C.A. & Smithson, S.B. 1985. A model for the MOIST seismic reflection profile, N Scotland. Journal of the Geological Society of London, 142, 245-258.

Brewer, J.A., Matthews, D.H., Warner M.R., Hall, J., Smythe, D.K. & Whittington, R.J. 1983. BIRPS deep seismic reflection studies of the British Caledonides. Nature 305, 206 - 210.

Brewer, J.A. & Smythe, D. K. 1984. MOIST and the continuity of crustal reflector geometry along the Caledonian-Appalachian orogen. Journal of the Geological Society of London, 141, 105-120.

Butler, R. W. H. 2010. The role of thrust tectonic models in understanding structural evolution in NW Scotland. Geological Society, London, Special Publications, 335, 293-320.

Butler, R. W. H., Holdsworth , R. E. & Matthews, S. 2006. Styles of basement involvement in the Moine Thrust Belt, NW Scotland. In: Mazzoli, S. and Butler, R. W. H. (eds) Styles of continental contraction. Geological Society of America, Special Paper, 414.

Clutson M.J., Brown D.E. & Tanner L.H. 2018. Distal processes and effects of multiple Late Triassic terrestrial bolide impacts: Insights from the Norian Manicouagan Event, northeastern Quebec, Canada. In: Tanner L. (ed.) The Late Triassic World. Topics in Geobiology, 46, 127-187. Springer, Cham.

Collins, G.S., Melosh, H.J. & Marcus, R.A. 2005. Earth Impact Effects Program: A Web-based computer program for calculating the regional environmental consequences of a meteoroid impact on Earth. Meteoritics and Planetary Science, 40, 817–840.

Cowan, D.R. & Cooper, G.R.J. 2005. Enhancement of magnetic signatures of impact structures. In: Kenkmann, T, Hörz, F., Deutsch, A. (eds.) Large Meteorite Impacts III, Geological Society of America, Special Papers, 384, 51-65.

Coward, M.P. 1980. The Caledonian thrust and shear zones of N.W. Scotland. Journal of Structural Geology, 2, 11–17. 13

Coward, M. P. 1983. The thrust and shear zones of the Moine Thrust Zone and the NW Scottish Caledonides. Journal of the Geological Society of London, 140, 795–811.

Degeai, J-P. & Peulvast, J-P., 2006. Calcul de l’erosion à long terme en region de socle autour de grands astroblèmes du Quebec et de France. Géographi physique et Quaternaire, 60, 131-148.

Elliott, D., & Johnson, M.R.W. 1980, Structural evolution in the northern part of the Moine thrust belt, NW Scotland. Transactions of the Royal Society of Edinburgh, 71, 69–96.

Ernstson, K. 1984. A gravity-derived model for the Steinheim impact structure. International Journal of Earth Sciences, 73, 483-498.

Ernstson, K. & Fiebag, J. 1992. The (Spain): New insights from geophysical and geological investigations. Geologische Rundschau, 81, 403–425.

French, B.M. 1998. : A Handbook of Shock-Metamorphic Effects in Terrestrial Meteorite Impact Structures. Lunar and Planetary Institute, Contribution no. 954, Lunar and Planetary Institute, Houston, Texas, USA. 120 pp.

Hergarten, S. & Kenkmann, T., 2015. The number of impact craters on Earth: Any room for further discoveries? Earth and Planetary Science Letters, 425, 187-192. doi: 10.1016/j.epsl.2015.06.009.

Hildebrand, A.R., Pilkington, M., Ortiz-Aleman, C., Chavez, R., Fucugauchi, J.U., Connors, M., Graniel-Castro, E., Camara-Zi, A., Halpenny, J.F. & Niehaus, D. 1998. Mapping structure with gravity and seismic reflection data. Geological Society, London, Special Publications, 140, 155-176.

Kahle, H.-G., 1969. Abschätzung der Störungsmasse im Nördlinger Ries. Zeitschrift für Geophysik, 33, 317-345

Kenkmann, T., Collins, G. S. & Wünnemann, K. 2013. The modification stage of crater formation. In: Osinski, G. R. & Pierazzo, E. (eds.), Impact Cratering. Processes and Products. Wiley & Sons, 60-75.

Kenkmann, T. and von Dalwigk, I. 2000. Radial transpression ridges: A new structural feature of complex impact craters. Meteoritics and Planetary Science, 35, 1189-1201.

Leslie, A.G., Krabbendam, M., Kimbell, G.S. & Strachan, R.A. 2010. Regional-scale lateral variation and linkage in ductile thrust architecture: the Oykel Transverse Zone, and mullions in the Moine Nappe, NW Scotland. In: Law, R.D., Butler, R.W.H., Holdsworth, R.E., Krabbendam, M. and Strachan, R.A. (eds.) Continental tectonics and mountain building: The legacy of Peach and Horne. Geological Society, London, Special Publications, 335, 359-381.

McClay, K.R. & Coward, M.P. 1981.The Moine Thrust Zone: an overview. Geological Society, London, Special Publications, 9, 241-260.

Melosh, H. Jay. 1989. Impact cratering: A geologic process. Research supported by NASA. New York, Oxford University Press (Oxford Monographs on Geology and Geophysics, No. 11, 1989, 253 pp.

Melosh, H.J. 2015. Peak-ring craters and multiring basins. Bridging the Gap III, Abstract 1003.pdf. Morgan, J. V., Gluick, S.P.S. et al. 2016. The formation of peak rings in large impact craters. Science, 354, 878-882.

14

Norman, J.W. 1984. Tectonic effects of old very large meteoritic impacts on Earth showing on satellite imagery: a review and speculations. Journal of Structural Geology, 6, 737-747.

Osinski, G. R., Tornabene, L.L. & Grieve, R.A.F. 2011. Impact ejecta emplacement on terrestrial planets. Earth and Planetary Science Letters, 310, 167-181.

Osinski, G.R. & Pierazzo, E. (eds) 2013. Impact cratering: processes and products. Wiley-Blackwell, Chichester, xiii + 316 pp. Peach, B.N., Horne, J., Gunn, W., Clough, C.T., Hinxman, L.W. & Teall, J.J.H. 1907. The Geological Structure of the NW Highlands of Scotland. Geological Survey of Great Britain, Memoir, 668 pp.

Pilkington, M. and Grieve, R.A.F. 1992. The geophysical signature of terrestrial impact craters. Reviews of Geophysics, 30, 161-181.

Pilkington, M., Pesonen, L.J., Grieve, R.A.F. & Masaitis, V.L., 2002. Geophysics and Petrophysics of the Popigai Impact Structure. In: Plado, J. & Pesonen, L.J. (eds) Impacts in Precambrian Shields, Springer, Heidelberg, 87-107. DOI 10.1007/978-3-662-05010-1_4

Pohl, J. 2015. Modelling the gravity anomaly of the Rochechouart impact structure. Bridging the Gap III, abstract 1052 (https://www.hou.usra.edu/meetings/gap2015/pdf/1052.pdf).

Pohl, J., Stoeffler, D., Gall, H., & Ernstson, K. 1977. The Ries impact crater. In: Roddy, D.J., Pepn, R.O. & Merrill, R.B. (eds) Impact and explosion cratering: Planetary and terrestrial implications. Proceedings of the Symposium on Planetary Cratering Mechanics, Flagstaff, Ariz., September 13-17, 1976. New York, Pergamon Press. 343-404.

Pohl, J., Ernstson, K. & Lambert, P. 1978. Gravity measurements in the Rochechouart impact structure (France). Meteoritics, 13, 601-604.

Rebolledo-Vieyra, M. & Urrutia-Fucugauchi, J. 2004. Magnetostratigraphy of the impact breccias and postimpact carbonates from borehole Yaxcopoil1, Chicxulub impact crater, Yucatán, Mexico. Meteoritics and Planetary Science, 39, 821-829.

Reddy, S.M., Johnson, T.E., Fischer, S., Rickard, W.D.A. & Taylor, R. J. M. 2015. Precambrian reidite discovered in shocked zircon from the Stac Fada , Scotland. Geology, 43, 899-902.

Renne, P.R., Arenillas, I., Arz, J.A., Vajda, V., Gilabert, V. & Bermúdez, H.D. 2018. Multi-proxy record of the Chicxulub impact at the Cretaceous-Paleogene boundary from Gorgonilla Island, Colombia. Geology, 46, 547-550.

Rollin, K., 2009. Regional geophysics of northern Scotland. British Geological Survey, CD-ROM.

Sapers, H.M., Osinski, G.R., Banerjee, N.R., Ferrière, L., Lambert, P. & Izawa, M.R.M. 2014. Revisiting the Rochechouart impact structure, France. Meteoritics and Planetary Science, 49, 2152- 2168.

Sharpton, V.L., Burke, K., Camargo-Zanoguera, A., Hall, S.A., Lee, D.S., Marín, L.E., Suáarez- Reynoso, G., Quezada-Muñeton, J.M., Spudis, P.D. & Urrutia-Fucugauchi, J. 1993. Chicxulub multiring impact basin: size and other characteristics derived from gravity analysis. Science, 261, 1564-1567.

Schedl, A. 2015. Searching for distal ejecta on the craton: The sedimentary effects of meteorite impact. Journal of Geology, 123, 201-232.

15

Simms, M.J. 2015. The Stac Fada impact ejecta deposit and the Lairg Gravity Low: Evidence for a buried Precambrian impact crater in Scotland? Proceedings of the Geologists' Association, 126, 742- 761.

Simms, M.J. 2016. Reply to comment by Trevor Faulkner on “The Stac Fada impact ejecta deposit and the Lairg Gravity Low: evidence for a buried Precambrian impact crater in Scotland? [Proc. Geol. Assoc. 126, 742–761 (2015)] and the consequence for the formation of the caves within the Durness Limestone outcrops at , ”. Proceedings of the Geologists' Association, 127, 109.

Soper, N. J. & Barber, A. J. 1982. A model for the deep structure of the Moine thrust zone. Journal of the Geological Society of London, 139, 127-138.

Stewart, A. D. 1982. Late Proterozoic rifting in NW Scotland: the genesis of the ‘Torridonian’. Journal of the Geological Society of London, 139, 415–420.

Stewart, A.D. 2002. The later Proterozoic Torridonian rocks of Scotland; their sedimentology, geochemistry and origin. Geological Society, London, Memoir, 24, 5-21.

Sweeny, J.F., 1978. Gravity study of great impact. Journal of Geophysical Research, 83, 2809-2815.

Watson, J. & Dunning, F.W. 1979. Basement-cover relations in the British Caledonides. In: Harris, A.L., Holland, C.H. & Leake, B.E. (eds) Caledonides of the British Isles – reviewed. Geological Society, London, Special Publication, 8, 67-92.

Wünnemann, K., Morgan, J.V. & Jödicke, H. 2005. Is Ries crater typical for its size? An analysis based upon old and new geophysical data and numerical modeling. In: Kenkmann, T, Hörz, F. & Deutsch, A. (eds.) Large Meteorite Impacts III, Geological Society of America, Special Papers, 384, 67-84.

16

Figure captions

Fig 1. Regional geology of northern Scotland showing the outcrop of the Stoer Group and its relationship to the residual gravity field (contoured at 2mGal intervals) for the Lairg Gravity Low (from Rollin, 2009). Arrows indicate directional azimuths within the Stac Fada Member.

17

Fig. 2. Gravity stations of the British Geological Survey in northern Scotland. The rectangle centred with respect to the Lairg Gravity Low (cross) frames the station used for the re-evaluation of the Lairg gravity field.

18

Fig. 3. Bouguer anomaly map for the framed gravity stations in Fig. #1 displaying the roughly circular Lairg negative anomaly. Please note that because of the sea the most northwesterly and southeasterly parts of the Bouguer map lack any gravity stations.

19

Fig. 4. A regional trend field computed from the Bouguer anomalies in Fig. #2 by radical moving average low-pass filtering.

20

Fig. 5. Bouguer residual anomalies for the Lairg gravity field, with contours at 1 mgal intervals. The residual field results from subtracting the regional field (Fig. #3) from the measured Bouguer field (Fig. #2). Note that in the residual field the central negative anomaly is enclosed by a roughly circular ring of relatively positive anomalies. Dashed red line is circle of radius 25 km centred on the gravity low.

21

Fig. 6. Three dashed lines selected for diametric gravity profiles across the Lairg gravity minimum (see Fig. #6). From the circle a diameter of nearly 50 km for the ring of positive anomalies can be deduced.

22

Fig. 7. The Bouguer gravity profiles taken from the Bouguer map in Fig. #5 revealing roughly similar shape with regard to the central negative anomaly and the surrounding ring of relatively positive anomalies that are only faintly developed on the NNW - SSE profile. Also note the ring of continued gravity lows beyond the suggested peak ring (arrowed) strongly supporting the peak ring character of a much larger Lairg impact structure.

23

Fig. 8. Results of 2.5D model calculation of a very simple two-layer density model for the NW - SE and the SW - NE gravity profiles. For reasons of simplicity, and because effectively the gravity evaluation lacks any "true" zero level, a constant regional shift of -10 mgal has been applied to the curves in Fig. 7.

Fig. 9. Rough approximation of the mass deficit related to the central anomaly by a spherical segment of -0.15 g/cm3 density defecit. Model without vertical exaggeration.

24

Fig. 10. Gravity residual anomalies of impact structures of various sizes compared with the Lairg residual anomaly. Note the different scales. Sources: Chicxulub (modified from Hildebrand et al. 1998), Manicouagan (modified from Sweeny 1978), Popigai (modified from Pilkington 2002), Lairg (this paper), Chesapeake (digitised from Gravity Map, ), Rochechouart (modified from Schmidt 1984), Ries (digitised from gravity map in Kahle 1969).

25

Fig 11. Some peak-ring impact structures with an approximate double ratio of crater/peak-ring diameter.

26

Fig. 12. Gravity-derived mass deficits of terrestrial impact structures as a function of diameter. Modified and supplemented from Ernstson and Fiebag (1992; with data from Pohl et al. [1978] and references therein; Ernstson [1984]).

27

Fig. 13. Schematic sequence of the main events allowing accommodation of a 100 km diameter impact crater centred upon the Lairg Gravity Low. 1170 Ma ago. A recently formed ~100 km diameter peak ring impact crater, surrounded by a thick and extensive ejecta deposit blanketing the sediments of the Stoer Group, is filled with lake sediments above a primary fill of impact breccia and impact melt. 1000 Ma ago. Prolonged erosion, prior to deposition of the Formation (Neoproterozoic, Torridon Group), removes the outer rim of the crater and much of the ejecta deposit and Stoer Group. Just a remnant survives in a downfaulted block far to the west. 400 Ma ago. During the Caledonian Orogeny thin-skinned thrusting emplaces westwards a thick cover of Moinian metasediments across the region. Thick-skinned thrusts extend deep into the Lewisian basement, perhaps nucleated on fractures associated with the original impact, and translocate the impact crater tens of kilometres westwards to its present position relative to the Stoer Group outcrop. Today. Post-Caledonian erosion has stripped away the Moinian metasediments in the west to expose the Stoer Group. The impact crater remains deeply buried beneath Lairg.

28