Outcrop-scale manifestations of reactivation during multiple superimposed rifting and basin inversion events: the Devonian Orcadian Basin, northern Scotland

The Devonian Orcadian Basin in Scotland hosts extensional fault systems assumed to be related to the initial formation of the basin, with only limited post-Devonian inversion and reactivation. However, a recent detailed structural study across Caithness, underpinned by published Re–Os geochronology, shows that three phases of deformation are present. North–south- and NW–SE-trending Group 1 faults are related to Devonian ENE–WSW transtension associated with sinistral shear along the Great Glen Fault during the formation of the Orcadian Basin. Metre- to kilometre-scale north–south-trending Group 2 folds and thrusts are developed close to earlier sub-basin-bounding faults and reflect late Carboniferous–early Permian east–west inversion associated with dextral reactivation of the Great Glen Fault. The dominant Group 3 structures are dextral oblique NE–SW-trending and sinistral east–west-trending faults with widespread syndeformational carbonate mineralization (± pyrite and bitumen) and are dated using Re–Os geochronology as Permian (c. 267 Ma). Regional Permian NW–SE extension related to the development of the offshore West Orkney Basin was superimposed over pre-existing fault networks, leading to local oblique reactivation of Group 1 faults in complex localized zones of transtensional folding, faulting and inversion. The structural complexity in surface outcrops onshore therefore reflects both the local reactivation of pre-existing faults and the superimposition of obliquely oriented rifting episodes during basin development in the adjacent offshore areas

Active normal faulting during the 1997 seismic sequence in Colfiorito, Umbria: Did slip propagate to the surface?

In order to determine whether slip during an earthquake on the 26th September 1997 propagated to the surface, structural data have been collected along a bedrock fault scarp in Umbria, Italy. These collected data are used to investigate the relationship between the throw associated with a debated surface rupture (observed as a pale unweathered stripe at the base of the bedrock fault scarp) and the strike, dip and slip-vector. Previous studies have suggested that the surface rupture was produced either by primary surface slip or secondary compaction of hangingwall sediments. Some authors favour the latter because sparse surface fault dip measurements do not match nodal plane dips at depth. It is demonstrated herein that the strike, dip and height of the surface rupture, represented by a pale unweathered stripe at the base of the bedrock scarp, shows a systematic relationship with respect to the geometry and kinematics of faulting in the bedrock. The strike and dip co-vary and the throw is greatest where the strike is oblique to the slip-vector azimuth where the highest dip values are recorded. This implies that the throw values vary to accommodate spatial variation in the strike and dip of the fault across fault plane corrugations, a feature that is predicted by theory describing conservation of strain along faults, but not by compaction. Furthermore, published earthquake locations and reported fault dips are consistent with the analysed surface scarps when natural variation for surface dips and uncertainty for nodal plane dips at depth are taken into account. This implies that the fresh stripe is indeed a primary coseismic surface rupture whose slip is connected to the seismogenic fault at depth. We discuss how this knowledge of the locations and geometry of the active faults can be used as an input for seismic hazard assessment

Igneous sills record far-field and near-field stress interactions during volcano construction: Isle of Mull, Scotland

Sill emplacement is typically associated with horizontally mechanically layered host rocks in a near-hydrostatic far-field stress state, where contrasting mechanical properties across the layers promote transitions from dykes, or inclined sheets, to sills. We used detailed field observations from the Loch Scridain Sill Complex (Isle of Mull, UK), and mechanical models to show that layering is not always the dominant control on sill emplacement. The studied sills have consistently shallow dips (1°–25°) and cut vertically bedded and foliated metamorphic basement rocks, and horizontally bedded cover sedimentary rocks and lavas. Horizontal and shallowly-dipping fractures in the host rock were intruded with vertical opening in all cases, whilst steeply-dipping discontinuities within the sequence (i.e. vertical fractures and foliation in the basement, and vertical polygonal joints in the lavas) were not intruded during sill emplacement. Mechanical models of slip tendency, dilation tendency, and fracture susceptibility for local and overall sill geometry data, support a radial horizontal compression during sill emplacement. Our models show that dykes and sills across Mull were emplaced during NW–SE horizontal shortening, related to a far-field tectonic stress state. The dykes generally accommodated phases of NE–SW horizontal tectonic extension, whereas the sills record the superposition of the far-field stress with a near-field stress state, imposed by emplacement of the Mull Central Volcano. We show that through detailed geometric characterisation coupled with mechanical modelling, sills may be used as an indication of fluctuations in the paleostress state

An evaluation of Mesozoic rift-related magmatism on the margins of the Labrador Sea: Implications for rifting and passive margin asymmetry

The Labrador Sea is a small (∼900 km wide) ocean basin separating southwest Greenland from Labrador, Canada. It opened following a series of rifting events that began as early as the Late Triassic or Jurassic, culminating in a brief period of seafloor spreading commencing by polarity chron 27 (C27; Danian) and ending by C13 (Eocene-Oligocene boundary). Rift-related magmatism has been documented on both conjugate margins of the Labrador Sea. In southwest Greenland this magmatism formed a major coast-parallel dike swarm as well as other smaller dikes and intrusions. Evidence for rift-related magmatism on the conjugate Labrador margin is limited to igneous lithologies found in deep offshore exploration wells, mostly belonging to the Alexis Formation, along with a postulated Early Cretaceous nephelinite dike swarm (ca. 142 Ma) that crops out onshore, near Makkovik, Labrador. Our field observations of this Early Cretaceous nephelinite suite lead us to conclude that the early rift-related magmatism exposed around Makkovik is volumetrically and spatially limited compared to the contemporaneous magmatism on the conjugate southwest Greenland margin. This asymmetry in the spatial extent of the exposed onshore magmatism is consistent with other observations of asymmetry between the conjugate margins of the Labrador Sea, including the total sediment thickness in offshore basins, the crustal structure, and the bathymetric profile of the shelf width. We propose that the magmatic and structural asymmetry observed between these two conjugate margins is consistent with an early rifting phase dominated by simple shear rather than pure shear deformation. In such a setting Labrador would be the lower plate margin to the southwest Greenland upper plate

A comparison of terrestrial laser scanning and structure-from-motion photogrammetry as methods for digital outcrop acquisition

Terrestrial laser scanning (TLS) has been used extensively in Earth Science for acquisition of digital outcrop data over the past decade. Structure-from-­motion (SfM) photogrammetry has recently emerged as an alternative and competing technology. The real-world performance of these technologies for ground-based digital outcrop acquisition is assessed using outcrops from North East England and the United Arab Emirates. Both TLS and SfM are via­ble methods, although no single technology is universally best suited to all situations. There are a range of practical considerations and operating conditions where each method has clear advantages. In comparison to TLS, SfM benefits from being lighter, more compact, cheaper, more easily replaced and repaired, with lower power requirements. TLS in comparison to SfM provides intrinsically validated data and more robust data acquisition in a wide range of operating conditions. Data post-processing is also swifter. The SfM data sets were found to contain systematic inaccuracies when compared to their TLS counterparts. These inaccuracies are related to the triangulation approach of the SfM, which is distinct from the time-of-flight principle employed by TLS. An elaborate approach is required for SfM to produce comparable results to TLS under most circumstances

Determining Histories of Slip on Normal Faults With Bedrock Scarps Using Cosmogenic Nuclide Exposure Data

Cosmogenic exposure data can be used to calculate time-varying fault slip rates on normal faults with exposed bedrock scarps. The method relies on assumptions related to how the scarp is preserved, which should be consistent at multiple locations along the same fault. Previous work commonly relied on cosmogenic data from a single sample locality to determine the slip rate of a fault. Here we show that by applying strict sampling criteria and using geologically informed modeling parameters in a Bayesian-inference Markov chain Monte Carlo method, similar patterns of slip rate changes can be modeled at multiple sites on the same fault. Consequently, cosmogenic data can be used to resolve along-strike fault activity. We present cosmogenic 36Cl concentrations from seven sites on two faults in the Italian Apennines. The average slip rate varies between sites on the Campo Felice Fault (0.84 ± 0.23 to 1.61 ± 0.27 mm yr−1), and all sites experienced a period of higher than average slip rate between 0.5 and 2 ka and a period of lower than average slip rate before 3 ka. On the Roccapreturo fault, slip rate in the center of the fault is 0.55 ± 0.11 and 0.35 ± 0.05 mm yr−1 at the fault tip near a relay zone. The estimated time since the last earthquake is the same at each site along the same fault (631 ± 620 years at Campo Felice and 2,603 ± 1,355 years at Roccapreturo). These results highlight the potential for cosmogenic exposure data to reveal the detailed millennial history of earthquake slip on active normal faults

Spatial migration of temporal earthquake clusters driven by the transfer of differential stress between neighbouring fault/shear-zone structures

Uncertainty concerning the processes responsible for slip-rate fluctuations associated with temporal clustering of surface faulting earthquakes is a fundamental, unresolved issue in tectonics, because strain-rates accommodated by fault/shear-zone structures are the key to understanding the viscosity structure of the crust and seismic hazard. We constrain the timing and amplitude of slip-rate fluctuations that occurred on three active normal faults in central Italy over a time period of 20–30 kyrs, using in situ 36Cl cosmogenic dating of fault planes. We identify five periods of rapid slip on individual faults lasting a few millennia, separated time periods of up to 10 millennia with low or zero slip-rate. The rapid slip pulses migrated across the strike between the faults in two waves from SW to NE. We replicate this migration with a model where rapid slip induces changes in differential stress that drive changes in strain-rate on viscous shear zones that drive slip-rate variability on overlying brittle faults. Earthquakes increase the differential stress and strain-rate on underlying shear zones, which in turn accumulate strain, re-loading stress onto the overlying brittle fault. This positive feedback produces high strain-rate episodes containing several large magnitude surface faulting earthquakes (earthquake clusters), but also reduce the differential stress on the viscous portions of neighbouring fault/shear-zones slowing the occurrence of large-magnitude surface faulting earthquakes (earthquake anticlusters). Shear-zones on faults experiencing anticlusters continue to accumulate viscous strain at a lowered rate, and eventually this loads the overlying brittle fault to failure, initiating a period of rapid slip through the positive feedback process described above, and inducing lowered strain-rates onto neighbouring fault/shear-zones. We show that these patterns of differential stress change can replicate the measured earthquake clustering implied by the 36Cl data. The stress changes are related to the fault geometry in terms of distance and azimuth from the slipping structure, implying that (a) strain-rate and viscosity fluctuations for studies of continental rheology, and (b) slip-rates for seismic hazard purposes are to an extent predictable given knowledge of the fault system geometry