Combining shock barometry with numerical modeling: Insights into complex crater formation—The example of the Siljan impact structure (Sweden)

Siljan, central Sweden, is the largest known impact structure in Europe. It was formed at about 380 Ma, in the late Devonian period. The structure has been heavily eroded to a level originally located underneath the crater floor, and to date, important questions about the original size and morphology of Siljan remain unanswered. Here we present the results of a shock barometry study of quartz‐bearing surface and drill core samples combined with numerical modeling using iSALE. The investigated 13 bedrock granitoid samples show that the recorded shock pressure decreases with increasing depth from 15 to 20 GPa near the (present) surface, to 10–15 GPa at 600 m depth. A best‐fit model that is consistent with observational constraints relating to the present size of the structure, the location of the downfaulted sediments, and the observed surface and vertical shock barometry profiles is presented. The best‐fit model results in a final crater (rim‐to‐rim) diameter of ~65 km. According to our simulations, the original Siljan impact structure would have been a peak‐ring crater. Siljan was formed in a mixed target of Paleozoic sedimentary rocks overlaying crystalline basement. Our modeling suggests that, at the time of impact, the sedimentary sequence was approximately 3 km thick. Since then, there has been around 4 km of erosion of the structure.     

Characterization and morphological reconstruction of the Terny impact structure,central Ukraine

The Terny impact structure, located in central Ukraine, displays a variety of diagnostic indicators of shock metamorphism, including shatter cones, planar deformation features in quartz, diaplectic glass, selective melting of minerals, and whole rock melting. The structure has been modified by erosion and subsequently buried by recent sediments. Although there are no natural outcrops of the deformed basement rocks within the area, mining exploration has provided surface and subsurface access to the structure, exposing impact melt rocks, shocked parautochthonous target rocks, and allochthonous impact breccias, including impact melt‐bearing breccias similar to suevites observed at the Ries structure. We have collected and studied samples from surface and subsurface exposures to a depth of approximately 750 m below the surface. This analysis indicates the Terny crater is centered on geographic coordinates 48.13° N, 33.52° E. The center location and the distribution of shock pressures constrain the transient crater diameter to be no less than approximately 8.4 km. Using widely accepted morphometric scaling relations, we estimate the pre‐erosional rim diameter of Terny crater to be approximately 16–19 km, making it close in original size to the well‐preserved El'gygytgyn crater in Siberia. Comparison with El'gygytgyn yields useful insights into the original morphology of the Terny crater and indicates that the amount of erosion Terny experienced prior to burial probably does not exceed 320 m.     

The Keurusselkä impact structure,Finland—Impact origin confirmed by characterization of planar deformation features in quartz grains

Abstract– Although the meteorite impact origin of the Keurusselkä impact structure (central Finland) has been established on the basis of the occurrence of shatter cones, no detailed microscopic examination of the impactites from this structure has so far been made. Previous microscope investigations of in situ rocks did not yield any firm evidence of shock features (Raiskila et al. 2008; Kinnunen and Hietala 2009). We have carried out microscopic observations on petrographic thin sections from seven in situ shatter cone samples and report here the discovery of planar fractures (PFs) and planar deformation features (PDFs) in quartz and feldspar grains. The detection and characterization of microscopic shock metamorphic features in the investigated samples substantiates a meteorite impact origin for the Keurusselkä structure. The crystallographic orientations of 372 PDF sets in 276 quartz grains were measured, using a universal stage (U‐stage) microscope, for five of the seven distinct shatter cone samples. Based on our U‐stage results, we estimate that investigated shatter cone samples from the Keurusselkä structure have experienced peak shock pressures from approximately 2 GPa to slightly less than 20 GPa for the more heavily shocked samples. The decoration of most of the PDFs with fluid inclusions also indicates that these originally amorphous shock features were altered by postimpact processes. Finally, our field observations indicate that the exposed surface corresponds to the crater floor; it is, however, difficult to estimate the exact diameter of the structure and the precise amount of material that has been eroded since its formation.     

Identification of the Bloody Creek structure,a possible impact crater in southwestern Nova Scotia,Canada

Abstract— An approximately 0.4 km diameter elliptical structure formed in Devonian granite in southwestern Nova Scotia, herein named the Bloody Creek structure (BCS), is identified as a possible impact crater. Evidence for an impact origin is based on integrated geomorphic, geophysical, and petrographic data. A near‐continuous geomorphic rim and a 10 m deep crater that is infilled with lacustrine sediments and peat define the BCS. Ground penetrating radar shows that the crater has a depressed inner floor that is sharply ringed by a 1 m high buried scarp. Heterogeneous material under the floor, interpreted as deposits from collapse of the transient cavity walls, is overlain by stratified and faulted lacustrine and wetland sediments. Alteration features found only in rim rocks include common grain comminution, polymict lithic microbreccias, kink‐banded feldspar and biotite, single and multiple sets of closely spaced planar microstructures (PMs) in quartz and feldspar, and quartz mosaicism, rare reduced mineral birefringence, and chlorite showing plastic deformation and flow microtextures. Based on their form and crystallographic orientations, the quartz PMs consist of planar deformation features that document shock‐metamorphic pressures ≤25 GPa. The age of the BCS is not determined. The low depth to diameter ratio of the crater, coupled with anomalously high shock‐metamorphic pressures recorded at its exposed rim, may be a result of significant post‐impact erosion. Alternatively, impact onto glacier ice during the waning stages of Wisconsinian deglaciation (about 12 ka BP) may have resulted in dissipation of much impact energy into the ice, resulting in the present morphology of the BCS.     

Cloud Creek structure,central Wyoming,USA: Impact origin confirmed

Abstract— The circular Cloud Creek structure in central Wyoming, USA is buried beneath ?1200 m of Mesozoic sedimentary rocks and has a current diameter of ?7 km. The morphology/morphometry of the structure, as defined by borehole, seismic, and gravity data, is similar to that of other buried terrestrial complex impact structures in sedimentary target rocks, e.g., Red Wing Creek in North Dakota, USA. The structure has a fault‐bordered central peak with minimum diameter of ?1.4 km, composed predominantly of Paleozoic carbonates thickened by thrust faulting and brecciation, and is elevated some 520 m above equivalent strata beyond the outer rim of the structure. There is a ?1.6 km wide annular trough sloping away from the central peak (maximum structural relief, 300 m) and terminated by a detached, fault‐bounded, rim anticline. The youngest rocks within the structure are Late Triassic (Norian?) clastics and these are overlain unconformably by post‐impact Middle Jurassic (Bathonian?) sandstones and shales. Thus, the formation of the Cloud Creek structure is dated chronostratigraphicly as ?190 ± 20 Ma. Reported here for the first time are measurements of planar deformation features (PDFs) in shocked quartz grains in thin sections made from drill cuttings recovered in a borehole drilled at the southern perimeter of the central peak. Other, less definitive microstructures consistent with impact occur in samples collected from boreholes drilled into the central peak and rim anticline. The shock‐metamorphic evidence confirms an impact origin for the Cloud Creek structure.     

The Foelsche structure,Northern Territory,Australia: An impact crater of probable Neoproterozoic age

Abstract— The Foelsche structure is situated in the McArthur Basin of northern Australia (16°40′ S, 136°47′ E). It comprises a roughly circular outcrop of flat‐lying Neoproterozoic Bukalara Sandstone, overlying and partly rimmed by tangentially striking, discontinuous outcrops of dipping, fractured and brecciated Mesoproterozoic Limmen Sandstone. The outcrop expression coincides with a prominent circular aeromagnetic anomaly, which can be explained in terms of the local disruption and removal or displacement of a regional mafic igneous layer within a circular area at depth. Samples of red, lithic, pebbly sandstone from the stratigraphically lowest exposed levels of the Bukalara Sandstone within the Foelsche structure contain detrital quartz grains displaying mosaicism, planar fractures (PFs) and planar deformation features (PDFs). PFs and PDFs occur in multiple intersecting sets with orientations consistent with a shock metamorphic origin. The abundance and angular nature of the shocked grains indicates a nearby provenance. Surface expression and geophysical data are consistent with a partly buried complex impact crater of ?6 km in diameter with an obscured central uplift ?2 km in diameter. The deformed outcrops of Limmen Sandstone are interpreted as relics of the original crater rim, but the central region of the crater, from which the shocked grains were likely derived, remains buried. From the best available age constraints the Foelsche structure is most likely of Neoproterozoic age.     

Shocked quartz grains in the early Cambrian Vakkejokk Breccia,Sweden—Evidence of a marine impact

Here we present a study of the abundance and orientation of planar deformation features (PDFs) in the Vakkejokk Breccia, a proposed lower Cambrian impact ejecta layer in the North‐Swedish Caledonides. The presence of PDFs is widely accepted as evidence for shock metamorphism associated with cosmic impact events and their presence confirms that the Vakkejokk Breccia is indeed the result of an impact. The breccia has previously been divided into four lithological subunits (from bottom to top), viz. lower polymict breccia (LPB), graded polymict breccia (GPB), top sandstone (TS), and top conglomerate (TC). Here we show that the LPB contains no shock metamorphic features, indicating that the material derives from just outside of the crater and represents low‐shock semi‐autochthonous bombarded strata. In the overlying, more fine‐grained GPB and TS, quartz grains with PDFs are relatively abundant (2–5% of the grain population), and with higher shock levels in the upper parts, suggesting that they have formed by reworking of more distal ejecta by resurge of water toward the crater in a marine setting. The absence of shocked quartz grains in the TC indicates that this unit represents later slumps associated with weathering and erosion of the protruding crater rim. Sparse shocked quartz grains (<0.2%) were also found in sandstone beds occurring at the same stratigraphic level as the Vakkejokk Breccia 15–20 km from the inferred crater site. It is currently unresolved whether the sandstone at these distal sites is related to the impact or just contains rare reworked quartz grains with PDFs.     

Shock‐metamorphic microfeatures in chert from the Jebel Waqf as Suwwan impact structure,Jordan

Abstract– The petrographic investigation of a shocked, chalcedony‐, quartzine‐, and quartz‐bearing allochthonous chert nodule (probably Upper Cretaceous) recovered from surficial wadi gravels in the inner parts of the central uplift of the approximately 6 km in diameter Jebel Waqf as Suwwan impact structure, Jordan, reveals new potential shock indicators in microfibrous–spherulitic silica, in addition to well‐established shock‐metamorphic effects in coarser crystalline quartz. The microcrystalline chert groundmass exhibits a macroscopic dendritic and suborthogonal fracture pattern commonly associated with thin “recrystallization bands” that intersect the pre‐existing diagenetic chert fabric. Fibrous aggregates of quartzine spherulites in chalcedony‐quartzine‐quartz veinlets locally have a shattered appearance and show conspicuous “curved fractures” perpendicular to the quartzine fiber direction (and parallel to [0001]) that commonly trend subparallel to planar fractures (PFs) in neighboring shocked quartz. Quartz exhibits PFs, feather features (FFs), and mainly single sets of planar deformation features (PDFs) parallel to the basal plane (0001) (Brazil twins) and, rarely, additional PDFs parallel to {101¯3}. Shock petrography indicates shock pressures of ≥10 GPa and high shock‐induced differential stresses that affected the chert nodule. The internal crosscutting relationships of primary diagenetic and impact‐related deformational features together with shockpressure estimates suggest that the curved fractures across quartzine spherulites might represent specific (low‐ to medium‐pressure) shock‐metamorphic features, possibly in structural analogy to basal plane PFs in quartz. The dendritic–suborthogonal fractures in the microcrystalline chert groundmass and recrystallization bands are likely related to impact‐induced shear deformation and recrystallization, respectively, and cannot be considered as definite shock indicators.     

Shatter cones and planar deformation features confirm Santa Marta in Piauí State,Brazil, as an impact structure

A total of 184 confirmed impact structures are known on Earth to date, as registered by the Earth Impact Database . The discovery of new impact structures has progressed in recent years at a rather low rate of about two structures per year. Here, we introduce the discovery of the approximately 10 km diameter Santa Marta impact structure in Piauí State in northeastern Brazil. Santa Marta is a moderately sized complex crater structure, with a raised rim and an off‐center, approximately 3.2 km wide central elevated area interpreted to coincide with the central uplift of the impact structure. The Santa Marta structure was first recognized in remote sensing imagery and, later, by distinct gravity and magnetic anomalies. Here, we provide results obtained during the first detailed ground survey. The Bouguer anomaly map shows a transition from a positive to a negative anomaly within the structure along a NE–SW trend, which may be associated with the basement signature and in parts with the signature developed after the crater was formed. Macroscopic evidence for impact in the form of shatter cones has been found in situ at the base around the central elevated plateau, and also in the interior of fractured conglomerate boulders occurring on the floor of the surrounding annular basin. Planar deformation features (PDFs) are abundant in sandstones of the central elevated plateau and at scattered locations in the inner part of the ring syncline. Together, shatter cones and PDFs provide definitive shock evidence that confirms the impact origin of Santa Marta. Crystallographic orientations of PDFs occurring in multiple sets in quartz grains are indicative of peak shock pressures of 20–25 GPa in the rocks exposed at present in the interior of the crater. In contrast to recent studies that have used additional, and sometimes highly controversial, alleged shock recognition features, Santa Marta was identified based on well‐understood, traditional shock evidence.     

Shocked quartz in polymict impact breccia from the Upper Cretaceous Yallalie impact structure in Western Australia

Yallalie is a ~12 km diameter circular structure located ~200 km north of Perth, Australia. Previous studies have proposed that the buried structure is a complex impact crater based on geophysical data. Allochthonous breccia exposed near the structure has previously been interpreted as proximal impact ejecta; however, no diagnostic indicators of shock metamorphism have been found. Here we report multiple (27) shocked quartz grains containing planar fractures (PFs) and planar deformation features (PDFs) in the breccia. The PFs occur in up to five sets per grain, while the PDFs occur in up to four sets per grain. Universal stage measurements of all 27 shocked quartz grains confirms that the planar microstructures occur in known crystallographic orientations in quartz corresponding to shock compression from 5 to 20 GPa. Proximity to the buried structure (~4 km) and occurrence of shocked quartz indicates that the breccia represents either primary or reworked ejecta. Ejecta distribution simulated using iSALE hydrocode predicts the same distribution of shock levels at the site as those found in the breccia, which supports a primary ejecta interpretation, although local reworking cannot be excluded. The Yallalie impact event is stratigraphically constrained to have occurred in the interval from 89.8 to 83.6 Ma based on the occurrence of Coniacian clasts in the breccia and undisturbed overlying Santonian to Campanian sedimentary rocks. Yallalie is thus the first confirmed Upper Cretaceous impact structure in Australia.     

Shock metamorphism investigations of quartz grains in clasts from impact breccia of the Eyreville B drill core,Chesapeake Bay impact structure,USA

Abstract– The Chesapeake Bay impact structure, approximately 85 km in diameter, has been drilled in 2005–2006 at Eyreville (Virginia, USA), to a total depth of 1766 m. In the drill cores, the abundance of shock metamorphosed material is very variable with depth. Shocked mineral and lithic clasts, as well as melt particles, are most abundant in suevitic impact breccia section (1397–1451 m depth). Shocked quartz (i.e., quartz grains with planar fractures and/or planar deformation features) and melt particles, although rare, are also dispersed in the Exmore Formation unit (444–867 m depth). Other lithologies in the Eyreville drill cores show no clear evidence of shock metamorphism. Here, we report on the investigations of 40 samples from the impact breccia section. A total of more than 27,000 quartz grains were examined in about 200 clasts. The abundance of highly shocked clasts tends to decrease with increasing depth. Crystalline clasts derived from the crystalline basement are commonly only slightly shocked (contain generally <10 rel% of shocked quartz grains). The clasts of metamorphosed sediments show a low proportion of shocked quartz grains (mostly <10 rel%). Sedimentary clasts show a wide range of proportions of shocked quartz grains, with several of them being highly shocked clasts (most values between 0 and 40 rel%). Conglomerates show the highest proportion of shocked quartz grains of all types of clasts (up to 83 rel%). Polycrystalline quartz clasts are also commonly highly shocked (contain mostly between 10 and 40 rel% of shocked quartz grains). These hard nonporous clasts are possibly more liable to show evidence of shock. The investigations suggest that the intensity of shock metamorphism is the result of several parameters, such as original position in the target (both horizontal and vertical) and the properties of each lithology (e.g., grain size, porosity, and amount of matrix). According to the universal‐stage investigations, the dominant orientations of planar deformation features in quartz are , , and also .     

Original size of the Vredefort Structure: Implications for the geological evolution of the Witwatersrand Basin

Abstract— Historically, there have been a range of diameter estimates for the large, deeply eroded Vredefort impact structure within the Witwatersrand Basin, South Africa. Here, we estimate the diameter of the transient cavity at the present level of erosion as ～124–140 km, based on the spatial distribution of shock metamorphic features in the floor of the structure and downfaulted Transvaal outliers. Taking erosion into account (<6 km) and scaling to original final rim diameter, an estimate of close to 300 km for the rim diameter is obtained. Independent estimates of the final rim diameter, based on an empirical relation of central uplift diameter to rim diameter, spatial distribution of pseudotachylites, and concentric large scale structural patterns, give a similar estimate of close to 300 km for the original final rim diameter. An impact structure of this size is expected to have had an original multi-ring form. At this size, the Vredefort impact structure encompasses the bulk of the Witwatersrand Basin, which appears to owe its preservation to the Vredefort impact. In addition, the Vredefort impact event may have been the thermal driver for some of the widespread hydrothermal activity in the area, which, in recent interpretations, is believed to be a component in the creation of the world-class gold deposits of the Witwatersrand Basin.     

Documentation of shock features in impactites from the Dhala impact structure,India

The fundamental approach for the confirmation of any terrestrial meteorite impact structure is the identification of diagnostic shock metamorphic features, together with the physical and chemical characterization of impactites and target lithologies. However, for many of the approximately 200 confirmed impact structures known on Earth to date, multiple scale‐independent tell‐tale impact signatures have not been recorded. Especially some of the pre‐Paleozoic impact structures reported so far have yielded limited shock diagnostic evidence. The rocks of the Dhala structure in India, a deeply eroded Paleoproterozoic impact structure, exhibit a range of diagnostic shock features, and there is even evidence for traces of the impactor. This study provides a detailed look at shocked samples from the Dhala structure, and the shock metamorphic evidence recorded within them. It also includes a first report of shatter cones that form in the shock pressure range from ~2 to 30 GPa, data on feather features (FFs), crystallographic indexing of planar deformation features, first‐ever electron backscatter diffraction data for ballen quartz, and further analysis of shocked zircon. The discovery of FFs in quartz from a sample of the MCB‐10 drill core (497.50 m depth) provides a comparatively lower estimate of shock pressure (~7–10 GPa), whereas melting of a basement granitoid infers at least 50–60 GPa shock pressure. Thus, the Dhala impactites register a strongly heterogeneous shock pressure distribution between <2 and >60 GPa. The present comprehensive review of impact effects should lay to rest the nonimpact genesis of the Dhala structure proposed by some earlier workers from India.     

Shocked quartz grains from the Målingen structure,Sweden—Evidence for a twin crater of the Lockne impact structure

The Målingen structure in Sweden has for a long time been suspected to be the result of an impact; however, no hard evidence, i.e., shock metamorphic features or traces of the impactor, has so far been presented. Here we show that quartz grains displaying planar deformation features (PDFs) oriented along crystallographic planes typical for shock metamorphism are present in drill core samples from the structure. The shocked material was recovered from basement breccias, below the sediment infill, and the distribution of the orientation of the shock‐produced PDFs indicates that the studied material experienced low shock pressures. Based on our findings, we can exclude that the material is transported from the nearby Lockne impact structure, which means that the Målingen structure is a separate impact structure, the seventh confirmed impact structure in Sweden. Furthermore, sedimentological and biostratigraphic aspects of the deposits that fill the depression at Målingen are very similar to features at the Lockne impact structure. This implies a coeval formation age and thus also the confirmation of the first known marine target doublet impact craters on Earth (i.e., the Lockne–Målingen pair).     

The Aouelloul crater,Mauritania: On the problem of confirming the impact origin of a small crater

Abstract— The impact origin of small craters in sedimentary rocks is often difficult to confirm because of the lack of characteristic shock metamorphic features. A case in point is the 3.1 Ma Aouelloul crater (Mauritania), 390 m in diameter, which is exposed in an area of Ordovician Oujeft and Zli sandstone. We studied several fractured sandstone samples from the crater rim for the possible presence of shock metamorphic effects. In thin section, a large fraction of the quartz grains show abundant subplanar and planar fractures. Many of the fractures are healed and are evident only as fluid inclusion trails. A few grains showed sets of narrow and densely spaced fluid inclusions trails in one (rarely two) orientations per grain, which could be possible remnants of planar deformation features (PDFs), although such an interpretation is not unambiguous. In contrast, an impact origin of the crater is confirmed by Re-Os isotope studies of the target sandstone and glass found around the crater rim, which show the presence of a distinct extraterrestrial component in the glass.     

The geology of the Målingen structure: A probable doublet to the Lockne marine‐target impact crater,central Sweden

The Målingen structure is an approximately 700 m wide, rimmed, sediment‐filled, circular depression in Precambrian crystalline basement approximately 16.2 km from the concentric, marine‐target Lockne crater (inner, basement crater diameter approximately 7.5 km, total diameter in sedimentary strata approximately 13.5 km). We present here results from geologic mapping, a 148.8 m deep core drilling from the center of the structure, detailed biostratigraphic dating of the structure's formation and its age correlation with Lockne, chemostratigraphy of the sedimentary infill, and indication for shock metamorphism in quartz from breccias below the crater infill. The drill core reveals, from bottom to the top, approximately 33 m of basement rocks with increased fracturing upward, approximately 10 m of polymict crystalline breccia with shock features, approximately 97 m of slumped Cambrian mudstone, approximately 4.7 m of a normally graded, polymict sedimentary breccia that in its uppermost part grades into sandstone and siltstone (cf. resurge deposits), and approximately 1.6 m of secular sediments. The combined data set shows that the Målingen structure formed in conjunction with the Lockne crater in the same marine setting. The shape and depth of the basement crater and the cored sequence of crystalline breccias with shocked quartz, slumped sediments, and resurge deposits support an impact origin. The stratigraphic and geographic relationship with Lockne suggests the Lockne and Målingen craters to be the first described doublet impact structure by a binary asteroid into a marine‐target setting.     

Complex crater formation: Insights from combining observations of shock pressure distribution with numerical models at the West Clearwater Lake impact structure

Large impact structures have complex morphologies, with zones of structural uplift that can be expressed topographically as central peaks and/or peak rings internal to the crater rim. The formation of these structures requires transient strength reduction in the target material and one of the proposed mechanisms to explain this behavior is acoustic fluidization. Here, samples of shock‐metamorphosed quartz‐bearing lithologies at the West Clearwater Lake impact structure, Canada, are used to estimate the maximum recorded shock pressures in three dimensions across the crater. These measurements demonstrate that the currently observed distribution of shock metamorphism is strongly controlled by the formation of the structural uplift. The distribution of peak shock pressures, together with apparent crater morphology and geological observations, is compared with numerical impact simulations to constrain parameters used in the block‐model implementation of acoustic fluidization. The numerical simulations produce craters that are consistent with morphological and geological observations. The results show that the regeneration of acoustic energy must be an important feature of acoustic fluidization in crater collapse, and should be included in future implementations. Based on the comparison between observational data and impact simulations, we conclude that the West Clearwater Lake structure had an original rim (final crater) diameter of 35–40 km and has since experienced up to ~2 km of differential erosion.     

LUKANGA SWAM: PROBABLY ASTROBLEME

This article summarizes information relevant to the interpretation of Lukanga Swamp — a depression some 52 km in diameter in central Zambia — as a considerably eroded astrobleme. The results of reconnaissance geological mapping, an airborne magnetic survey and petrological study of a small set of samples of impact rocks are presented. Circular topographic features coincide with a remarkably subdued magnetic field. A circular set of tangential faults indicates a diameter of the outer rim approximately 110 km. Polymict breccia occur, carrying marble megaclasts (klippes) up to 50 m long. Strongly shocked quartz clasts with several sets of decorated planar features are present as well as devitrified, flow-banded silica glass exhibiting decarbonation textures around marble clasts. Provisional correlation of the brecciation event with the local geological record suggests a Paleozoic age of the astrobleme.