Genesis of the mafic granophyre of the Vredefort impact structure (South Africa): Implications of new geochemical and Se and Re-Os isotope data

ABSTRACT This contribution is concerned with the debated origin of the impact melt rock in the central uplift of the world’s largest confirmed impact structure—Vredefort (South Africa). New major- and trace-element abundances, including those of selected highly siderophile elements (HSEs), Re-Os isotope data, as well as the first Se isotope and Se-Te elemental systematics are presented for the felsic and mafic varieties of Vredefort impact melt rock known as “Vredefort Granophyre.” In addition to the long-recognized “normal” (i.e., felsic, >66 wt% SiO2) granophyre variety, a more mafic (<66 wt% SiO2) impact melt variety from Vredefort has been discussed for several years. The hypothesis that the mafic granophyre was formed from felsic granophyre through admixture (assimilation) of a mafic country rock component that then was melted and assimilated into the superheated impact melt has been pursued here by analysis of the two granophyre varieties, of the Dominion Group lava (actually metalava), and of epidiorite mafic country rock types. Chemical compositions, including high-precision isotope dilution–derived concentrations of selected highly siderophile elements (Re, Os, Ir, Pt, Se, Te), and Re-Os and Se isotope data support this hypothesis. A first-order estimate, based on these data, suggests that some mafic granophyre may have resulted from a significant admixture (assimilation) of epidiorite to felsic granophyre. This is in accordance with the findings of an earlier investigation using conventional isotope (Sr-Nd-Pb) data. Moreover, these outcomes are in contrast to a two-stage emplacement model for Vredefort Granophyre, whereby a mafic phase of impact melt, derived by differentiation of a crater-filling impact melt sheet, would have been emplaced into earlier-deposited felsic granophyre. Instead, all chemical and isotopic evidence so far favors formation of mafic granophyre by local assimilation of mafic country rock—most likely epidiorite—by a single intrusive impact melt phase, which is represented by the regionally homogeneous felsic granophyre.

Systematic survey of K, Th, and U signatures in airborne radiometric data from Australian meteorite impact structures: Possible causes of circular features and implications

ABSTRACT Airborne radiometric (gamma-ray) data provide estimates of the concentrations of potassium (K), thorium (Th), and uranium (U) in soil, regolith, and bedrock. Radiometric data constitute an important source of geochemical information, commonly used in mineral exploration and for geological mapping of Earth and other planets. Airborne radiometric data have rarely been applied to the exploration and analyses of impact structures, in contrast with other conventional geophysical tools (e.g., gravimetry, magnetism, and seismic reflection/refraction). This work represents the first systematic survey of the K, Th, and U radiometric signatures of Australian impact structures, based on the continent-wide airborne radiometric coverage of Australia. We first formulated several hypotheses regarding the possible causes of formation of circular radiometric patterns associated with impact structures. Then, the radiometric signatures of 17 exposed impact structures in Australia were documented. Our observations confirmed the supposition that impact structures are commonly associated with circular radiometric patterns. We then selected the five structures with the most prominent circular radiometric patterns (Gosses Bluff, Lawn Hill, Acraman, Spider, and Shoemaker), and we discuss the possible origin of these anomalies. Based on these five case studies, we argue that such patterns result from either crustal deformation induced by the impact event and/or from postimpact superficial processes controlled by the crater topography. This work also suggests that airborne radiometric data may be useful, in combination with other geophysical tools, in the search for new possible impact structures.

Anisotropy of magnetic susceptibility (AMS) of impact melt breccia and target rocks from the Dhala impact structure, India

ABSTRACT The ~11-km-wide, Paleoproterozoic Dhala impact structure in north-central India comprises voluminous exposures of impact melt breccia. These outcrops are discontinuously spread over a length of ~6 km in a semicircular pattern along the northern, inner limit of the monomict breccia ring around the central elevated area. This study of the magnetic fabrics of impact breccias and target rocks from the Dhala impact structure identified a weak preferred magnetic orientation for pre-impact crystalline target rocks. The pre- and synimpact rocks from Dhala have magnetite and ilmenite as common magnetic phases. The distributions of magnetic vectors are random for most impact melt breccia samples, but some do indicate a preferred orientation. Our anisotropy of magnetic susceptibility (AMS) data demonstrate that the shape of susceptibility ellipsoids for the target rocks varies from prolate to oblate, and most impact melt breccia samples display both shapes, with a slight bias toward the oblate geometry. The average value for the corrected degree of anisotropy of impact melt rock (P′ = 1.009) is lower than that for the target rocks (P′ = 1.091). The present study also shows that both impact melt breccia and target rock samples of the Dhala structure have undergone minor postimpact alteration, and have similar compositions in terms of magnetic phases and high viscosity. Fine-grained iron oxide or hydroxide is the main alteration phase in impact melt rocks. Impact melt rocks gave a narrow range of mean magnetic susceptibility (Km) and P′ values, in contrast to the target rock samples, which gave Km = 0.05–12.9 × 10−3 standard international units (SI) and P′ = 1.036–1.283. This suggests similar viscosity of the source magma, and limited difference in the degrees of recorded deformation. Between Pagra and Maniar villages, the Km value of impact melt breccias gradually decreases in a clockwise direction, with a maximum value observed near Pagra (Km = 1.67 × 10−3 SI). The poor grouping of magnetic fabrics for most impact melt rock samples implies local turbulence in rapidly cooled impact melt at the front of the melt flow immediately after the impact. The mean K1 for most impact melt samples suggests subhorizontal (<5°) flow in various directions. The average value of Km for the target rocks (4.41 × 10−3 SI) is much higher compared to the value for melt breccias (1.09 × 10−3 SI). The results of this study suggest that the melt breccias were likely part of a sheet-like body of sizeable extent. Our magnetic fabric data are also supported by earlier core drilling information from ~70 locations, with coring depths reaching to −500 m. Our extensive field observations combined with available widespread subsurface data imply that the impact melt sheet could have covered as much as 12 km2 in the Dhala structure, with an estimated minimum melt volume of ~2.4 km3.

Australian impact cratering record: Updates and recent discoveries

ABSTRACT There are currently 31 confirmed structures of impact origin in Australia. More than 49 additional structures have been proposed to have formed due to asteroid impact but await confirmation. Many discoveries have been made in Australia in the time since the last comprehensive review of the Australian impact cratering record was published in a peer-reviewed journal in 2005. These include further expanding the record of confirmed craters, and providing new insights into a variety of impact-related processes, such as shock deformation, phase transitions in accessory minerals, new impact age determinations, studies of oblique impacts, and more. This update is a review that focuses principally on summarizing discoveries made since 2005. Highlights since then include confirmation of five new Australian impact structures, identification of Earth’s oldest recognized impact structure, recognition of shock deformation in accessory minerals, discovery of the high-pressure phase reidite in Australia, determination of the links between impact craters and some ore deposits, and publication of the first generation of numerical hydrocode models for some Australian craters.

Dynamics of collapse of an impact central uplift: Evidence from folds and faults in the collar of the Vredefort Dome, South Africa

ABSTRACT Structural analysis of overturned metasedimentary strata of the lower Witwatersrand Supergroup in the inner collar of the Vredefort Dome reveals the presence of tangential folds and faults associated with the 2.02 Ga impact. The folds are distinct from previously identified subradially oriented, vertical to plunging-inclined, gentle folds that are interpreted as the products of convergent flow (constriction) during the initial stages of central uplift formation. The tangential folds comprise disharmonic, open, asymmetric, horizontal to plunging-inclined anticline-syncline pairs with centripetally dipping axial planes and right-way-up intermediate limbs. They display centripetal-down vergence (anticline radially outward of the syncline) that is consistent with steep inward-directed shear of the overturned strata. We attribute this kinematic pattern to subvertical collapse of the Vredefort central uplift during the latter stages of crater modification. The folds are cut by pseudotachylite-bearing steep to vertical tangential faults that display center-down slip of <10 m up to ~150 m. Both the tangential folds and the faults suggest that the large-scale overturning of strata related to outward collapse of the Vredefort central uplift was accompanied by a component of inward-directed collapse via layer-parallel shearing and folding, followed by faulting. Subradially oriented faults with conjugate strike separations of 1–2 km in the NNE collar of the dome suggest penecontemporaneous tangential extension of the inner collar rocks. This evidence indicates that second-order structures in the metasedimentary collar of the Vredefort Dome preserve a complex, multistage record of evolving strain associated with both initial convergent and upward flow (constriction) related to central uplift rise and later divergent and downward flow (flattening) linked to its collapse, and that centripetally directed collapse features may be important components of the structural inventory of very large central uplifts.

Petrographic and chemical studies of the Cretaceous-Paleogene boundary sequence at El Guayal, Tabasco, Mexico: Implications for ejecta plume evolution from the Chicxulub impact crater

ABSTRACT A combined petrographic and chemical study of ejecta particles from the Cretaceous-Paleogene boundary sequence of El Guayal, Tabasco, Mexico (520 km SW of Chicxulub crater), was carried out to assess their formation conditions and genetic relation during the impact process. The reaction of silicate ejecta particles with hot volatiles during atmospheric transport may have induced alteration processes, e.g., silicification and cementation, observed in the ejecta deposits. The various microstructures of calcite ejecta particles are interpreted to reflect different thermal histories at postshock conditions. Spherulitic calcite particles may represent carbonate melts that were quenched during ejection. A recrystallized microstructure may indicate short, intense thermal stress. Various aggregates document particle-particle interactions and intermixing of components from lower silicate and upper sedimentary target lithologies. Aggregates of recrystallized calcite with silicate melt indicate the consolidation of a hot suevitic component with sediments at ≳750 °C. Accretionary lapilli formed in a turbulent, steam-condensing environment at ~100 °C by aggregation of solid, ash-sized particles. Concentric zones with smaller grain sizes of accreted particles indicate a recurring exchange with a hotter environment. Our results suggest that during partial ejecta plume collapse, hot silicate compo nents were mixed with the fine fraction of local surface-derived sediments, the latter of which were displaced by the preceding ejecta curtain. These processes sustained a hot, gas-driven, lateral basal transport that was accompanied by a turbulent plume at a higher level. The exothermic back-reaction of CaO from decomposed carbonates and sulfates with CO2 to form CaCO3 may have been responsible for a prolonged release of thermal energy at a late stage of plume evolution.

Extreme plastic deformation and subsequent Pb loss in shocked xenotime from the Vredefort Dome, South Africa

ABSTRACT Accessory mineral U-Pb geochronometers are crucial tools for constraining the timing of deformation in a wide range of geological settings. Despite the growing recognition that intragrain age variations within deformed minerals can spatially correlate to zones of microstructural damage, the causal mechanisms of Pb loss are not always evident. Here, we report the first U-Pb data for shock-deformed xenotime, from a detrital grain collected at the Vredefort impact structure in South Africa. Orientation mapping revealed multiple shock features, including pervasive planar deformation bands (PDBs) that accommodate up to 40° of lattice misorientation by <100>{010} slip, and also an ~50-µm-wide intragrain shear zone that contains {112} deformation twin lamellae in two orientations. Twenty-nine in situ secondary ion mass spectrometry (SIMS) U-Pb analyses from all microstructural domains yielded a well-defined discordia with upper-intercept age of 2953 ± 15 Ma (mean square of weighted deviates [MSWD] = 0.57, n = 29, 2σ), consistent with derivation from Kaapvaal craton bedrock. However, the 1754 ± 150 Ma lower concordia intercept age falls between the 2020 Ma Vredefort impact and ca. 1100 Ma Kibaran orogenesis and is not well explained by multiple Pb-loss episodes. The pattern and degree of Pb loss (discordance) correlate with increased [U] but do not correlate to microstructure (twin, PDB) or to crystallinity (band contrast) at the scale of SIMS analysis. Numerical modeling of the Pb-loss history using a concordia-discordia-comparison (CDC) test indicated that the lower concordia age is instead best explained by an alteration episode at ca. 1750 Ma, rather than a multiple Pb-loss history. In this example, the U-Pb system in deformed xenotime does not record a clear signature of impact age resetting; rather, the implied high dislocation density recorded by planar deformation bands and the presence of deformation twins facilitated subsequent Pb loss during a younger event that affected the Witwatersrand basin. Microstructural characterization of xenotime targeted for geochronology provides a new tool for recognizing evidence of deformation and can provide insight into complex age data from highly strained grains, and, as is the case in this study, elucidate previously unrecognized alteration events.

U-Pb geochronology of apatite crystallized within a terrestrial impact melt sheet: Manicouagan as a geochronometer test site

ABSTRACT The Manicouagan impact event has been the subject of multiple age determinations over the past ~50 yr, providing an ideal test site for evaluating the viability of different geochronometers. This study highlights the suitability of Manicouagan’s essentially pristine impact melt body as a medium for providing insight into the U-Pb isotope systematics of geochronometers in the absence of shock-related overprinting. We performed in situ laser-ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) U-Pb geochronology on apatite and zircon, both of which crystallized as primary phases. This study is the first application of U-Pb geochronology to apatite crystallized within a terrestrial impact melt sheet. U-Pb analyses were obtained from 200 melt-grown apatite grains (n = 222 spots), with a data subset providing a lower-intercept age of 212.5 ± 8.0 Ma. For melt-grown zircon, a total of 30 analyses from 28 grains were obtained, with a subset of the data yielding a lower-intercept age of ± 1.6 Ma. The lower precision (±8.0 Ma; ±3%) obtained from apatite is a consequence of low U and a high and variable common-Pb composition. This resulted from localized Pb*/PbC heterogeneity within the impact melt sheet that was incorporated into the apatite crystal structure during crystallization (where Pb*/PbC is the ratio of radiogenic Pb to common Pb). While considered a limitation to the precision obtainable from melt-grown apatite, its ability to record local-scale isotopic variations highlights an advantage of U-Pb studies on melt-grown apatite. The best-estimate ages from zircon and apatite overlap within error and correlate with previously determined ages for the Manicouagan impact event. An average formation age from the new determinations, combined with previous age constraints, yields a weighted mean age of 214.96 ± 0.30 Ma for the Manicouagan impact structure.

Terrestrial and extraterrestrial chemical components of early Archean impact spherule layers from Fairview Gold Mine, northern Barberton greenstone belt, South Africa

ABSTRACT Early Archean spherule layers, widely accepted to represent distal ejecta deposits from large-scale impact events onto the early Earth, have been described from several stratigraphic levels of the Barberton greenstone belt in South Africa. Recently, exploration drilling at the Fairview Gold Mine (25°43′53″S, 31°5′59″E) in the northern domain of the belt resulted in the discovery of a new set of spherule layer intersections. The Fairview spherule layers in drill cores BH5901, BH5907, BH5911, and BH5949 were intersected just a few meters apart, at about the same stratigraphic position within the transition from the Onverwacht Group to the Fig Tree Group. The Fairview spherule layers have petrographic and chemical similarities to at least three other well-known Barberton spherule layers (S2–S4), and multiple spherule layer bed intersections in drill cores BARB5 and CT3, all from about the same stratigraphic position. They are not uniform in composition, in particular with respect to abundances of highly siderophile elements. The highest concentrations of moderately (Cr, Co, Ni) and highly siderophile (Ir) elements are within the range of concentrations for chondrites and, thus, reinforce the impact hypothesis for the generation of the Fairview spherule layers. Iridium peak concentrations and Cr/Ir interelement ratios for spherule layer samples from drill cores BH5907, BH5911, and BH5949 suggest admixtures of 50%–60% chondritic material, whereas for the BH5901 spherule layer, only an admixture of 1% chondritic material is indicated. We discuss whether these four Fairview spherule layers represent the same impact event, and whether they can be correlated to any of the S2–S4, CT3, and BARB5 intersections.

Revealing microstructural properties of shocked and tectonically deformed zircon from the Vredefort impact structure: Raman spectroscopy combined with SEM microanalyses

ABSTRACT Finite deformation patterns of accessory phases can indicate the tectonic regime and deformation history of the host rocks and geological units. In this study, tectonically deformed, seismically deformed, and shocked zircon grains from a granite sample from the core of the Vredefort impact structure were analyzed in situ, using a combination of Raman spectroscopy, backscatter electron (BSE) imaging, electron backscattered diffraction (EBSD) mapping, electron probe microanalyses (EPMA), energy-dispersive X-ray spectroscopy (EDS) qualitative chemical mapping, and cathodoluminescence (CL) imaging. We aimed to reveal the effects of marginal grain-size reduction, planar deformation bands (PDBs), and shock microtwins on the crystal structure and microchemistry of zircon. Deformation patterns such as PDBs, microtwins, and subgrains did not show any significant effect on zircon crystallinity/metamictization degree or on the CL signature. However, the ratio of Raman band intensities B1g (1008 cm–1) to Eg (356 cm–1) slightly decreased within domains with low misorientation. The ratio values were affected in shocked grains, particularly in twinned domains with high misorientation. B1g/Eg ratio mapping combined with metamictization degree mapping (full width at half maximum of B1g peak) suggest the presence of shock deformation features in zircon; however, due to the lower spatial resolution of the method, they must be used in combination with the EBSD technique. Additionally, we discovered anatase, quartz, goethite, calcite, and hematite micro-inclusions in the studied zircon grains, with quartz and anatase specifically being associated with strongly deformed domains of shocked zircon crystals.