The complex impact structure Serra da Cangalha,Tocantins State,Brazil

Abstract– Serra da Cangalha is a complex impact structure with a crater diameter of 13,700 m and a central uplift diameter of 5800 m. New findings of shatter cones, planar fractures, feather features, and possible planar deformation features are presented. Several ring‐like features that are visible on remote sensing imagery are caused by selective erosion of tilted strata. The target at Serra da Cangalha is composed of Devonian to Permian sedimentary rocks, mainly sandstones that are interlayered with siltstone and claystones. NNE–SSW and WNW–ESE‐striking joint sets were present prior to the impact and also overprinted the structure after its formation. As preferred zones of weakness, these joint sets partly controlled the shape of the outer perimeter of the structure and, in particular, affected the deformation within the central uplift. Joints in radial orientation to the impact center did not undergo a change in orientation during tilting of strata when the central uplift was formed. These planes were used as major displacement zones. The asymmetry of the central uplift, with preferred overturning of strata in the northern to western sector, may suggest a moderately oblique impact from a southerly direction. Buckle folding of tilted strata, as well as strata overturning, indicates that the central uplift became gravitationally unstable at the end of crater formation.     

An appraisal of the Serra da Cangalha impact structure using the Euler deconvolution method

Abstract— The applicability of the Euler deconvolution method in imaging impact crater structure vis‐à‐vis delineation of source depth of the circular magnetic anomaly and/or basement depth beneath the crater is addressed in this paper. The efficacy of the method has been evaluated using the aeromagnetic data obtained over the Serra da Cangalha impact crater, northeastern Brazil. The analyses of the data have provided characteristic Euler deconvolution signatures and structural indices associated with impact craters. Also, through the interpretation of the computed Euler solutions, our understanding of the structural features present around the impact structure has been enhanced. The Euler solutions obtained indicate shallow magnetic sources that are interpreted as possibly post‐impact faults and a circular structure. The depth of these magnetic sources varies between 0.8 and 2.5 km, while the Precambrian basement depth was found at ?1.5 km. This is in good agreement with the estimates of the Precambrian basement depth of about 1.1 km, calculated using aeromagnetic data. The reliability of the depth solutions obtained through the implementation of the Euler method was confirmed through the use of the existing information available in the area and the result of previous studies. We find that the Euler depth solutions obtained in this study are consistent with the results obtained using other methods.     

Investigation of Shuttle Radar Topography Mission data of the possible impact structure at Serra da Cangalha,Brazil

Abstract— The Serra da Cangalha crater structure in northeast Brazil, ?13 km in diameter, has long been widely considered to be a confirmed impact structure, based on reports of shatter cone findings. Only very limited field work has been carried out at this crater structure. Landsat Thematic Mapper (TM) and Shuttle Radar Topography Mission (SRTM) data sets for the region around this crater structure are compared here with regard to their suitability to determine first‐order structural detail of impact crater structures. The SRTM data provide very detailed information regarding drainage patterns and topography. A pronounced central ring of up to 300 m elevation above the surrounding area, two comparatively subdued intermediate rings of 6 and 10.5 km diameter, respectively, and the broad, complex crater rim of up to >100 m elevation can be distinguished in the Serra da Cangalha data. The maximum cratering‐related regional deformation (radial and concentric features) seems to be limited to a radial distance of 16–18 km from the center of the structure. A first comparison of macrostructural information from several impact structures with that from Serra da Cangalha does not yield firm trends, but the database is still very small at this stage. The varied nature of the target geology strongly influences the development of structural features in any impact event.     

Insights about the formation of a complex impact structure formed in basalt from numerical modeling: The Vista Alegre structure,southern Brazil

We present the outcomes of simulations of the formation of the Vista Alegre impact structure, Paraná Basin, Brazil. The target comprised a thick sequence of volcanic rocks of predominantly basaltic composition of the Serra Geral Formation that had been deposited on top of sedimentary rocks (sandstones) of the Pirambóia/Botucatu formations. The cratering process was modeled using the iSALE shock physics code. Our best‐fit model suggests that (1) the crater was originally ~10 km in size; (2) it was formed in ~115 s by a stony projectile of 1000 m in diameter, for an assumed impact velocity of 12 km s^?1; (3) target rocks underwent a peak pressure of ~20 GPa, in agreement with previous petrographic investigations of shock deformation. Furthermore, the model points out that the sedimentary strata below the layer of volcanic rocks were raised by ~650 meters at the central part of the crater, which resulted in the current partial exposure of the sandstones at the surface. The outcomes of our modeling suggest that parameters like cohesion and strength of the target rocks, after shock compression, determined the final morphology of the crater, especially the absence of a topographically prominent central peak. Finally, the results of the numerical modeling are roughly in agreement with gravity data over the structure, in particular with respect to the presence of the uplifted sedimentary strata, which are responsible for a low gravity signature at the center of the structure.     

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.     

Hydrocode modeling of the Sierra Madera impact structure

Abstract— We present the first hydrocode simulations of the formation of the Sierra Madera structure (west Texas, USA), which was caused by an impact into a thick sedimentary target sequence. We modeled Sierra Madera using the iSALE hydrocode, and here we present two best‐fit models: 1) a crater with a rim (final crater) diameter of ?12 km, in agreement with previous authors' interpretations of the original structure, and 2) a crater ?16 km in diameter with increased postimpact erosion. Both models fit some of the geologic observational data, but discrepancies with estimates of peak shock pressure, extent of deformation, and stratigraphic displacement remain. This study suggests that Sierra Madera may be a larger crater than previously reported and illustrates some of the challenges in simulating impact deformation of sedimentary lithologies. As many terrestrial craters possess some amount of sedimentary rocks in the target sequence, numerical models of impacts into sedimentary targets are essential to our understanding of target rock deformation and the mechanics of crater formation.     

Ceres internal structure from geophysical constraints

Thermal evolution modeling has yielded a variety of interior structures for Ceres, ranging from a modestly differentiated interior to more advanced evolution with a dry silicate core, a hydrated silicate mantle, and a volatile‐rich crust. Here we compute the mass and hydrostatic flattening from more than one hundred billion three‐layer density models for Ceres and describe the characteristics of the population of density structures that are consistent with the Dawn observations. We show that the mass and hydrostatic flattening constraints from Ceres indicate the presence of a high‐density core with greater than a 1σ probability, but provide little constraint on the density, allowing for core compositions that range from hydrous and/or anhydrous silicates to a mixture of metal and silicates. The crustal densities are consistent with surface observations of salts, water ice, carbonates, and ammoniated clays, which indicate hydrothermal alteration, partial fractionation, and the possible settling of heavy sulfide and metallic particles, which provide a potential process for increasing mass with depth.     

Numerical modeling of impact heating and cooling of the Vredefort impact structure

Abstract— Large meteorite impacts, such as the one that created the Vredefort structure in South Africa?2 Ga ago, result in significant heating of the target. The temperatures achieved in these events have important implications for post‐impact metamorphism as well as for the development of hydrothermal systems. To investigate the post‐impact thermal evolution and the size of the Vredefort structure, we have analyzed impact‐induced shock heating in numerical simulations of terrestrial impacts by projectiles of a range of sizes thought to be appropriate for creating the Vredefort structure. When compared with the extent of estimated thermal shock metamorphism observed at different locations around Vredefort, our model results support our earlier estimates that the original crater was 120–160 km in diameter, based on comparison of predicted to observed locations of shock features. The simulations demonstrate that only limited shock heating of the target occurs outside the final crater and that the cooling time was at least 0.3 Myr but no more than 30 Myr.     

Azimuthal impact directions from oblique impact crater morphology

Planetary impact craters have a high degree of radial symmetry. This hampers efforts to identify the azimuthal impact direction for most craters – the radially symmetric component of an impact crater swamps any asymmetries that may be present. We demonstrate how the asymmetric component can be isolated and the direction of the asymmetries quantified using a two-dimensional eigenfunction expansion over a circular domain. The complex coefficients of expansion describe the magnitude and phase (angular alignment) of each term. From the analysis of hypervelocity impact craters formed in the laboratory, with impact angles ranging from 0° to 50° from the surface normal, we show that asymmetries which reveal the impact direction are still present at just 10° from the surface normal, and that the phase of one complex coefficient of expansion, c ₁₁, indicates the impact direction. Analysis of the lunar crater Hadley shows bilateral symmetry in the radially asymmetric component, which may be due to oblique impact. The 31-km lunar ray crater Kepler has morphological features that indicate the azimuthal impact direction. Coefficient c ₁₁ gives an azimuthal impact direction similar to that expected from the morphology, although post-impact gravitational collapse and slumping obscure the result to some degree. Ray craters may provide a means of testing the method for smaller 'simple' craters when data are available.     

Numerical modeling of the formation of the Eltanin submarine impact structure

Two-and three-dimensional numerical simulations of the formation of the Eltanin submarine impact structure are described. Based on the numerical results, the impactor size can be estimated, its destruction and the subsequent motion of fragments can be described, and the initial amplitude of the tsunami wave can be determined.     

Insights into Ceres's evolution from surface composition

Inspired by the recent results of the Dawn mission, thermodynamic models of rock alteration and brine evaporation have been used to help understand the conditions under which water–rock interaction took place within the dwarf planet Ceres. This analysis constrains Ceres's early history and offers a framework within which future observations may be interpreted. A broad range of alteration conditions have been simulated using the Geochemist's Workbench and PHREEQC software, associated with the FREZCHEM model that constrains the consequences of freezing the liquid phase in equilibrium with the observed mineralogical assemblage. Comparison of the modeling results with observed surface mineralogy at Ceres indicates advanced alteration under a relatively high fugacity of hydrogen, a conclusion that is consistent with predictions for, and observations of, large ice‐rich bodies. The simulations suggest production of methane that could help regulate the redox environment and possibly form clathrate hydrates upon freezing of the early ocean. The detection of localized occurrences of natrite (sodium carbonate) at the surface of Ceres provides key constraints on the composition of fluids that are necessarily alkaline. In addition, the combined hydrothermal and freezing simulations suggest that hydrohalite may be abundant in Ceres's subsurface, similar to Earth's polar regions. The global homogeneity of Ceres's surface, made of material formed at depth, suggests a large‐scale formation mechanism, while local heterogeneities associated with impact craters and landslides suggest that some form of sodium carbonate and other salts are accessible in the shallow subsurface.     

Insights into the formation of silica‐rich achondrites from impact melts in Rumuruti‐type chondrites

Ancient, SiO₂‐rich achondrites have previously been proposed to have formed by disequilibrium partial melting of chondrites. Here, we test the alternative hypothesis that these achondrites formed by fractional crystallization of impact melts of Rumuruti (R) chondrites. We identified two new melt clasts in R chondrites, one in Pecora Escarpment (PCA) 91241 and one in LaPaz Icefield (LAP) 031275. We analyzed major, minor, and trace element concentrations, as well as oxygen isotopes, of these two clasts and a third one that had been previously recognized (Bischoff et al. 2011) as an impact melt in Dar al Gani (DaG) 013. The melt clast in PCA 91241 is an R chondrite impact melt closely resembling the one previously recognized in DaG 013. The melt clast in LAP 031275 has an L chondrite provenance. We show that SiO₂‐rich melts could form from the mesostases of R chondrite impact melts. However, their CI‐normalized rare earth element patterns are flat, whereas those of ancient SiO₂‐rich achondrites (Day et al. 2012; Srinivasan et al. 2018) and those of disequilibrium partial melts of chondrites (Feldstein et al. 2001) have positive Eu anomalies from preferential melting of plagioclase. Thus, we conclude that ancient SiO₂‐rich achondrites were probably formed by disequilibrium partial melting (due to an internal heat source on their parent bodies), rather than from impact melts.     

Neptune's satellites: Implications from impact crater morphology

The presence of central peak craters and the absence of central pit craters on Triton implies a surface rigidity similar to the Saturnian and Uranian satellites and stronger than that of the Jupiter satellites Ganymede and Callisto. Tectonically degraded terrain may exist at the antipode of the large impact structure on 1989N1. Dome craters on Triton may represent a form of solid state volcanism.     

Summanen structure: Further geological and geophysical evidence of a meteorite impact event in Central Finland

The Summanen structure is located in Central Finland and is one of Finland's 12 known meteorite impact structures. In 2017, the discovery of Summanen was based on numerous shatter cone boulders with planar deformation features (PDFs) and a circular electromagnetic anomaly, which is 2.6 km in diameter. The site was revisited in 2020 and 2022, and shatter cone-bearing outcrops were discovered. PDFs and feather features were identified in samples from these outcrops. A total of 38 PDF sets in 27 quartz grains resulted in rational crystallographic orientations concentrating on {10$\overline{1}$4}, {10$\overline{1}$3}, {10$\overline{1}$2}, and {11$\overline{2}$2}, implying shock pressures of 2–20 GPa. Gravity measurements were taken, and the electrical conductivity of the structure was studied. The gravimetric results revealed a circular negative anomaly of about 4 km in diameter, with an amplitude of −3.5 mGal. Excluding the gravitational effect of water and Quaternary sediments reduces the anomaly to −1.6 mGal. A bowl-shaped conductive layer, likely containing relict saline water in the impact-fractured bedrock, was identified to a depth of 240 m. Topographic and bathymetric data were combined to determine the impact's effect and interpret the level of erosion. Cobbles of sedimentary sand- and siltstones were found on the coastline of Lake Summanen. Based on their similarity to those found in the Söderfjärden impact crater with a Cambrian age, it is likely that these rocks and post-impact infill are also of a similar age.     

Through the impact glass: Insight into the evolution of melt at the Mistastin Lake impact structure

By analyzing impact glass, the evolution of the impact melt at the Mistastin Lake impact structure was investigated. Impact glass clasts are present in a range of impactites, including polymict breccias and clast‐rich impact melt rock, and from a variety of settings within the crater. From the glass clasts analyzed, three petrographic subtypes of impact glass were identified based on their clast content, prevalence of schlieren, color, texture, and habit. Several alteration phases were also observed replacing glass and infilling vesicles; however, textural observations and quantified compositional data allowed for the identification of pristine impact glass. Although the various types of glasses show significant overlap in their major oxide composition, several subtle variations in the major oxide chemistry of the glass were observed. To investigate this variation, a least‐squares mixing model was implemented utilizing the composition of the glass and the known target rock chemistry to model the initial melt composition. Additionally, image analysis of the glass clasts was used to investigate whether the compositional variations correlated to textural difference in the lithologies. We propose that the textural and compositional dichotomy observed is a product of the evolution, assimilation, and emplacement of the glass. The dichotomy is reflective of the melt either being ballistically emplaced (group 2 glasses: occurring in melt‐poor polymict breccias at lowermost stratigraphic position outside the transient crater) or the result of late‐stage melt flows (group 1 glasses, occurring in melt‐bearing polymict breccias and impact melt rocks at higher stratigraphic positions outside the transient crater).     

Eminescu impact structure: Insight into the transition from complex crater to peak-ring basin on Mercury

Peak-ring basins represent an impact-crater morphology that is transitional between complex craters with central peaks and large multi-ring basins. Therefore, they can provide insight into the scale dependence of the impact process. Here the transition with increasing crater diameter from complex craters to peak-ring basins on Mercury is assessed through a detailed analysis of Eminescu, a geologically recent and well-preserved peak-ring basin. Eminescu has a diameter (∼125 km) close to the minimum for such crater forms and is thus representative of the transition. Impact crater size-frequency distributions and faint rays indicate that Eminescu is Kuiperian in age, geologically younger than most other basins on Mercury. Geologic mapping of basin interior units indicates a distinction between smooth plains and peak-ring units. Our mapping and crater retention ages favor plains formation by impact melt rather than post-impact volcanism, but a volcanic origin for the plains cannot be excluded if the time interval between basin formation and volcanic emplacement was less than the uncertainty in relative ages. The high-albedo peak ring of Eminescu is composed of bright crater-floor deposits (BCFDs, a distinct crustal unit seen elsewhere on Mercury) exposed by the impact. We use our observations to assess predictions of peak-ring formation models. We interpret the characteristics of Eminescu as consistent with basin formation models in which a melt cavity forms during the impact formation of craters at the transition to peak ring morphologies. We suggest that the smooth plains were emplaced via impact melt expulsion from the central melt cavity during uplift of a peak ring composed of BCFD-type material. In this scenario the ringed cluster of peaks resulted from the early development of the melt cavity, which modified the central uplift zone.     

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.     

Constraints on the size of the Vredefort impact crater from numerical modeling

Abstract— The Vredefort structure in South Africa was created by a meteorite impact about two billion years ago. Since that time, the crater has been deeply eroded; so to estimate its original size, researchers have had to rely heavily upon comparison to other terrestrial impact structures. Recent estimates of the original crater diameter range from 160 km to as much as 400 km. In this study, we combined the capabilities of both hydrocode and finite-element modeling, using the former to predict where the pressure of an impact-generated shock wave would have been high enough to form planar deformation features (PDFs) and shatter cones and the latter to follow the subsequent displacement of these shock isobars during the collapse of the crater. We established constraints on the sizes of the projectile and the transient crater (and, thus, on the size of the final crater) by comparing the observed locations of PDFs around Vredefort to the results of our simulations of impacts by projectiles of various sizes. These simulations indicate that a rocky projectile with a diameter of ～10 km, impacting vertically at a velocity of 20 km/s generates shock pressures that are consistent with the distribution of PDFs around Vredefort. These projectile parameters correspond to a transient crater ～80 km in diameter or a final crater ～120–160 km in diameter. Allowing for uncertainties in our modeling procedures, we consider final craters 120 to 200 km in diameter to be consistent with the observed locations of PDFs at Vredefort. The shock pressure contour corresponding to the formation of shatter cones is almost horizontal near the surface, making the locations of these features less useful constraints on the crater size. However, they may provide a constraint on the amount of erosion that has occurred since the impact.     

Shocked rocks and impact glasses from the El'gygytgyn impact structure,Russia

Abstract— The El'gygytgyn impact structure is about 18 km in diameter and is located in the central part of Chukotka, arctic Russia. The crater was formed in volcanic rock strata of Cretaceous age, which include lava and tuffs of rhyolites, dacites, and andesites. A mid‐Pliocene age of the crater was previously determined by fission track (3.45 ± 0.15 Ma) and ⁴⁰Ar/³⁹Ar dating (3.58 ± 0.04 Ma). The ejecta layer around the crater is completely eroded. Shock‐metamorphosed volcanic rocks, impact melt rocks, and bomb‐shaped impact glasses occur in lacustrine terraces but have been redeposited after the impact event. Clasts of volcanic rocks, which range in composition from rhyolite to dacite, represent all stages of shock metamorphism, including selective melting and formation of homogeneous impact melt. Four stages of shocked volcanic rocks were identified: stage I (≤35 GPa; lava and tuff contain weakly to strongly shocked quartz and feldspar clasts with abundant PFs and PDFs; coesite and stishovite occur as well), stage II (35–45 GPa; quartz and feldspar are converted to diaplectic glass; coesite but no stishovite), stage III (45–55 GPa; partly melted volcanic rocks; common diaplectic quartz glass; feldspar is melted), and stage IV (>55 GPa; melt rocks and glasses). Two main types of impact melt rocks occur in the crater: 1) impact melt rocks and impact melt breccias (containing abundant fragments of shocked volcanic rocks) that were probably derived from (now eroded) impact melt flows on the crater walls, and 2) aerodynamically shaped impact melt glass “bombs” composed of homogeneous glass. The composition of the glasses is almost identical to that of rhyolites from the uppermost part of the target. Cobalt, Ni, and Ir abundances in the impact glasses and melt rocks are not or only slightly enriched compared to the volcanic target rocks; only the Cr abundances show a distinct enrichment, which points toward an achondritic projectile. However, the present data do not allow one to unambiguously identify a meteoritic component in the El'gygytgyn impact melt rocks.     

Coesite in suevites from the Chesapeake Bay impact structure

The occurrence of coesite in suevites from the Chesapeake Bay impact structure is confirmed within a variety of textural domains in situ by Raman spectroscopy for the first time and in mechanically separated grains by X‐ray diffraction. Microtextures of coesite identified in situ investigated under transmitted light and by scanning electron microscope reveal coesite as micrometer‐sized grains (1–3 μm) within amorphous silica of impact‐melt clasts and as submicrometer‐sized grains and polycrystalline aggregates within shocked quartz grains. Coesite‐bearing quartz grains are present both idiomorphically with original grain margins intact and as highly strained grains that underwent shock‐produced plastic deformation. Coesite commonly occurs in plastically deformed quartz grains within domains that appear brown (toasted) in transmitted light and rarely within quartz of spheroidal texture. The coesite likely developed by a mechanism of solid‐state transformation from precursor quartz. Raman spectroscopy also showed a series of unidentified peaks associated with shocked quartz grains that likely represent unidentified silica phases, possibly including a moganite‐like phase that has not previously been associated with coesite.