Depositional processes of impactites from the YAX-1 drill core in the Chicxulub impact structure inferred from vertical profiles of PDF orientations and grain size distributions of shocked quartz

Core samples from the Chicxulub impact structure provide insights into the formation processes of a shallow-marine-target, complex crater. Although previous studies investigated the impactites (generally suevitic and polymict breccias) of the Yaxcopoil-1 (YAX-1) drill core in the Chicxulub impact structure, the interpretation of its deposition remains controversial. Here, we analyze planar deformation features (PDFs), grain size, and abundance of shocked quartz throughout the YAX-1 impactite sequence (794–895 m in depth). PDF orientations of most quartz grains in YAX-1 impactites show a distribution of both low angles ({104}, {103}, {102}) and high angles (orientations higher than 55° to c-axis), while the lower part of the impactite sequence contains quartz showing only PDF orientations of low angles. High-abundance, coarse-grained shocked quartz is found from the lower to middle parts of the impactites, whereas it abruptly changes to low-abundance, fine-grained shocked quartz within the upper part. In the uppermost part of the impactites, repeated oscillations in contents of these two components are observed. PDF orientation pattern suggests most of the shocked quartz grains experienced a range of shock pressure, except two samples in the lower part of impactites, which experienced only a high level of shock. We suggest that the base and lower part of the impactite sequence were formed by ejecta curtain and melt surge deposits, respectively. Our results are also consistent with the interpretation that the middle part of the impactite sequence is fallback ejecta from the impact plume. Additionally, we support the contention that massive seawater resurges into the crater occurred during the deposition of the upper and uppermost part of the impactites.     

Saqqar: A 34 km diameter impact structure in Saudi Arabia

Here we present the first proof of an impact origin for the Saqqar circular structure in northwestern Saudi Arabia (Neville et al. 2014 ), with an apparent diameter of 34 km, centered at 29°35′N, 38°42′E. The structure is formed in Cambrian–Devonian siliciclastics and is unconformably overlain by undeformed Cretaceous and Paleogene sediments. The age of impact is not well constrained and lies somewhere between 410 and 70 Ma. The subsurface structure is constrained by 2‐D reflection seismic profiles and six drilled wells. First‐order structural features are a central uplift that rises approximately 2 km above regional datums, surrounded by a ring syncline. The crater rim is defined by circumferential normal faults. The central uplift and ring syncline correspond to a Bouguer gravity high and an annular ring‐like low, respectively. The wells were drilled within the central uplift, the deepest among them exceed 2 km depth. Sandstone core samples from these wells show abundant indicators of a shock metamorphic overprint. Planar deformation features (PDFs) were measured with orientations along (0001), {103}, and less frequently along {101} and {104}. Planar fractures (PFs) predominantly occur along (0001) and {101}, and are locally associated with feather features (FFs). In addition, some shocked feldspar grains and strongly deformed mica flakes were found. The recorded shock pressure ranges between 5 and 15 GPa. The preserved level of shock and the absence of an allochthonous crater fill suggest that Saqqar was eroded by 1–2 km between the Devonian and Maastrichtian. The documentation of unequivocal shock features proves the formation of the Saqqar structure by a hypervelocity impact event.     

Determination of the equilibrium fO2 in bulk samples of H,L, and LL ordinary chondrites by solid‐state electrochemistry

High‐temperature solid‐state electrochemistry techniques (EMF method) were used to measure the oxygen fugacity (fO₂) of the ordinary chondrites Ochansk (H4), Savtschenskoje (LL4), Elenovka (L5), Vengerovo (H5), and Kharkov (L6). The fO₂ results are presented in the form of the following equations: It was found that fO₂ regularly increases from H chondrites to LL chondrites. Measured fO₂ are ~1.5 higher than those previously calculated from mineral assemblages. Kharkov (L6) is a little more oxidized than Elenovka (L5) in agreement with the progressive oxidation model. At the same time, Ochansk (H4) is more oxidized than Vengerovo (H5) and exhibits a slightly different slope compared to other chondrites and at T > 1200 K, becomes more reduced than Kharkov (L6) or Elenovka (L5). Measured oxygen fugacity values of meteorites fall within (0.1–1.0)·log fO₂ of one another. The possible explanation of discrepancies between measured and calculated values is discussed.     

A new calibration to determine the closure temperatures of Fe‐Mg ordering in augite from nakhlites

Recently it has been shown that the relatively low closure temperature (T_c) of 500 (100)°C calculated for augite from Miller Range nakhlite (MIL 03346,13) using the available geothermometers would correspond to a slow cooling rate inconsistent with the petrologic evidence for an origin from a fast‐cooled lava flow. Moreover, previous annealing experiments combined with HR‐SC‐XRD on an augite crystal from MIL 03346 clearly showed that at 600 °C, the Fe²⁺‐Mg degree of order remained unchanged, thus suggesting that the actual T_c is close to this temperature. In order to clarify this discrepancy, we undertook an ex situ annealing experimental study at 700, 800, and 900 °C, until the equilibrium in the intracrystalline Fe²⁺‐Mg exchange is reached, using an augite crystal from Miller Range nakhlite (MIL 03346,13) with a composition of about En₃₆Fs₂₄Wo₄₀. These data allowed us to calculate the following new geothermometer calibration for Martian nakhlites: where The application of this new equation to other Martian nakhlites (NWA 988 and Nakhla) suggests that for augite with composition close to that of MIL 03346, the T_c is up to 170 °C higher with respect to the one calculated using the previous available geothermometer equation, thus suggesting a significantly faster cooling in agreement with petrologic evidence.     

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.     

Basaltic volcanism on the angrite parent body: Comparison with 4 Vesta

Carbon and nitrogen data from stepped combustion analysis of eight angrites, seven eucrites, and two diogenites, alongside literature data from a further nine eucrites and two diogenites, have been used to assess carbon and nitrogen incorporation and isotope fractionation processes on the angrite parent body (APB), for comparison with volatile behavior on the HED parent body (4 Vesta). A subset of the angrite data has been reported previously (Abernethy et al. 2013 ). Two separate families of volatile components were observed. They were (1) moderately volatile material (MVM), mostly combusting between ~500 and 750 °C and indistinguishable from terrestrial contamination and (2) refractory material (RM), mainly released above 750 °C and thought to be carbon (as ) and nitrogen (as N₂ or ) dissolved within the silicate lattice, fitting with the slightly oxidized (~IW to IW+2) angrite fO₂ conditions. Isotopic fractionation trends for carbon and nitrogen within the plutonic and basaltic (quenched) angrites suggest that the behavior of the two volatile elements is loosely coupled, but that the fractionation process differs between the two angrite subgroups. Comparison with results from eucrites and diogenites implies similarities between speciation of carbon and nitrogen on 4 Vesta and the APB, with the latter being more enriched in volatiles than the former.     

Chondrule formation during low-speed collisions of planetesimals: A hybrid splash–flyby framework

Chondrules probably formed during a small window of time ~1–4 Ma after CAIs, when most solid matter in the asteroid belt was already in the form of km-sized planetesimals. They are unlikely, therefore, to be “building blocks” of planets or abundant on asteroids, but more likely to be a product of energetic events common in the asteroid belt at that epoch. Laboratory experiments indicate that they could have formed when solids of primitive composition were heated to temperatures of ~1600 K and then cooled for minutes to hours. A plausible heat source for this is magma, which is likely to have been abundant in the asteroid belt at that time, and only that time, due to the trapping of ²⁶Al decay energy in planetesimal interiors. Here, we propose that chondrules formed during low-speed ( $\lesssim 1\mathrm{km}{\mathrm{s}}^{-1}$) collisions between large planetesimals when heat from their interiors was released into a stream of primitive debris from their surfaces. Heating would have been essentially instantaneous and cooling would have been on the dynamical time scale, $1/\sqrt{\left(\mathit{G\rho }\right)}$ ~30 min, where $\rho$ is the mean density of a planetesimal. Many of the heated fragments would have remained gravitationally bound to the merged object and could have suffered additional heating events as they orbited and ultimately accreted to its surface. This is a hybrid of the splash and flyby models: We propose that it was the energy released from a body's molten interior, not its mass, that was responsible for chondrule formation by heating primitive debris that emerged from the collision.     

Crystallization temperature determination of Itokawa particles by plagioclase thermometry with X‐ray diffraction data obtained by a high‐resolution synchrotron Gandolfi camera

The crystallization temperatures of Itokawa surface particles recovered by the space probe Hayabusa were estimated by a plagioclase geothermometer using sodic plagioclase triclinicity. The Δ131‐index required for the thermometer, which is the difference in X‐ray diffraction peak positions between the 131 and 11 reflections of plagioclase, was obtained by a high‐resolution synchrotron Gandolfi camera developed for the third generation synchrotron radiation beamline, BL15XU at SPring‐8. Crystallization temperatures were successfully determined from the Δ131‐indices for four particles. The observed plagioclase crystallization temperatures were in a range from 655 to 660 °C. The temperatures indicate crystallization temperatures of plagioclases in the process of prograde metamorphism before the peak metamorphic stage.     

Hydrogen and major element concentrations on 433 Eros: Evidence for an L‐ or LL‐chondrite‐like surface composition

A reanalysis of NEAR X‐ray/gamma‐ray spectrometer (XGRS) data provides robust evidence that the elemental composition of the near‐Earth asteroid 433 Eros is consistent with the L and LL ordinary chondrites. These results facilitated the use of the gamma‐ray measurements to produce the first in situ measurement of hydrogen concentrations on an asteroid. The measured value,  ppm, is consistent with hydrogen concentrations measured in L and LL chondrite meteorite falls. Gamma‐ray derived abundances of hydrogen and potassium show no evidence for depletion of volatiles relative to ordinary chondrites, suggesting that the sulfur depletion observed in X‐ray data is a surficial effect, consistent with a space‐weathering origin. The newfound agreement between the X‐ray, gamma‐ray, and spectral data suggests that the NEAR landing site, a ponded regolith deposit, has an elemental composition that is indistinguishable from the mean surface. This observation argues against a pond formation process that segregates metals from silicates, and instead suggests that the differences observed in reflectance spectra between the ponds and bulk Eros are due to grain size differences resulting from granular sorting of ponded material.     

Mg and Si isotopic fractionation patterns in types B1 and B2 CAIs: Implications for formation under different nebular conditions

Magnesium and silicon isotopic profiles across melilite grains in two type B1 and two type B2 calcium‐aluminum‐rich inclusions (CAIs) reveal differing but constant enrichments in heavy isotopes everywhere except ≤1000 μm from the CAI margins. There is no close correlation in the B1s or the B2s between isotopic composition and åkermanite content of the melilite, a measure of progressive igneous crystallization, yet such a correlation might be expected in a type B2: without a melilite mantle (as in B1s) to seal the interior off and prevent further evaporation, the melt would have maintained communication with the external gas. These observations indicate a model in which B1s and B2s solidified under differing conditions. The B2s solidified under lower hydrogen pressures ( ≤ 10^?4 – 10^?5 bars) than did B1s ( > 10^?4 bars), so surface volatilization was slower in the B2s and internal chemical and isotopic equilibrium was maintained over the interval of melilite crystallization. The outermost zones of the CAIs (≤1000 μm from the edge) are not consistently enriched in heavy isotopes relative to the interiors, as might be expected from diffusion‐limited surface evaporation of the melt. In all cases, the magnesium in the CAI margins is lighter than in the interiors. In one case, silicon in the margin also is lighter, but locally in some CAIs, it is isotopically heavier near the surface. If melt evaporation played a role in the formation of these outer zones, a later event in many cases caused isotopic re‐equilibration with an external and isotopically near‐normal reservoir.     

Bidirectional reflectance distribution function measurements of the Winchcombe meteorite using the Visible Oxford Space Environment Goniometer

A laboratory study was performed using the Visible Oxford Space Environment Goniometer in which the broadband (350–1250 nm) bidirectional reflectance distribution function (BRDF) of the Winchcombe meteorite was measured, across a range of viewing angles—reflectance: 0°–70°, in steps of 5°; incidence: 15°, 30°, 45°, and 60°; and azimuthal: 0°, 90°, and 180°. The BRDF dataset was fitted using the Hapke BRDF model to (1) provide a method of comparison to other meteorites and asteroids, and (2) to produce Hapke parameter values that can be used to extrapolate the BRDF to all angles. The study deduced the following Hapke parameters for Winchcombe: w = 0.152 ± 0.030, b = 0.633 ± 0.064, and h_S = 0.016 ± 0.008, demonstrating that it has a similar w value to Tagish Lake (0.157 ± 0.020) and a similar b value to Orgueil (0.671 ± 0.090). Importantly, the surface profile of the sample was characterized using an Alicona 3D® instrument, allowing two of the free parameters within the Hapke model φ and $\overline{\theta }$, which represent porosity and surface roughness, respectively, to be constrained as φ = 0.649 ± 0.023 and $\overline{\theta }$ = 16.113° (at 500 μm size scale). This work serves as part of the characterization process for Winchcombe and provides a reference photometry dataset for current and future asteroid missions.     

The cosmic‐ray exposure history of the Twannberg iron meteorite (IIG)

The Twannberg iron meteorite is one out of only six members of the group IIG. The combined noble gas and radionuclide data obtained in this new systematic study indicate that Twannberg with its ~570 recently recovered specimens was a large object with a preatmospheric radius in the range of ~2 m, which corresponds to ~250 × 10³ kg. The cosmic‐ray exposure age for Twannberg is 182 ± 45 Ma. The most surprising result is the long terrestrial age of T_terr =  ka, which is unexpected considering the humid conditions in Switzerland. However, this age is in accord with glaciation events, indicating that the less shielded samples from Mt. Sujet were found close to the position of the original strewn field, whereas the samples from Gruebmatt and Twannbach, which are from more shielded positions, were glacially transported to the east–northeast during the second last ice age (185–130 ka ago) from an original position west of Mt. Sujet.     

The Erbisberg drilling 2011: Implications for the structure and postimpact evolution of the inner ring of the Ries impact crater

The 26 km diameter Nördlinger Ries is a complex impact structure with a ring structure that resembles a peak ring. A first research drilling through this “inner crystalline ring” of the Ries was performed at the Erbisberg hill (SW Ries) to better understand the internal structure and lithology of this feature, and possibly reveal impact‐induced hydrothermal alteration. The drill core intersected the slope of a 22 m thick postimpact travertine mound, before entering 42 m of blocks and breccias of crystalline rocks excavated from the Variscan basement at >500 m depth. Weakly shocked gneiss blocks that show that shock pressure did not exceed 5 GPa occur above polymict lithic breccias of shock stage Ia (10–20 GPa), with planar fractures and planar deformation features (PDFs) in quartz. Only a narrow zone at 49.20–50.00 m core depth exhibits strong mosaicism in feldspar and {102} PDFs in quartz, which are indicative of shock stage Ib (20–35 GPa). Finally, 2 m of brecciated Keuper sediments at the base of the section point to an inverse layering of strata. While reverse grading of clast sizes in lithic breccias and gneiss blocks is consistent with lateral transport, the absence of diaplectic glass and melt products argues against dynamic overthrusting of material from a collapsing central peak, as seen in the much larger Chicxulub structure. Indeed, weakly shocked gneiss blocks are rather of local provenance (i.e., the transient crater wall), whereas moderately shocked polymict lithic breccias with geochemical composition and ⁸⁷Sr/⁸⁶Sr signature similar to Ries suevite were derived from a position closer to the impact center. Thus, the inner ring of the Ries is formed by moderately shocked polymict lithic breccias likely injected into the transient crater wall during the excavation stage and weakly shocked gneiss blocks of the collapsing transient crater wall that were emplaced during the modification stage. While the presence of an overturned flap is not evident from the Erbisberg drilling, a survey of all drillings at or near the inner ring point to inverted strata throughout its outer limb. Whether the central ring of the Ries represents remains of a collapsed central peak remains to be shown. Postimpact hydrothermal alteration along the Erbisberg section comprises chloritization, sulfide veinlets, and strong carbonatization. In addition, a narrow zone in the lower parts of the polymict lithic breccia sequence shows a positive Eu anomaly in its carbonate phase. The surface expression of this hydrothermal activity, i.e., the travertine mound, comprises subaerial as well as subaquatic growth phases. Intercalated lake sediments equivalent to the early parts of the evolution of the central crater basin succession confirm a persistent impact‐generated hydrothermal activity, although for less time than previously suggested.     

Planar deformation features in quartz from impact‐produced polymict breccia of the Xiuyan crater,China

Abstract– The 1.8 km‐diameter Xiuyan crater is an impact structure in northeastern China, exposed in a Proterozoic metamorphic rock complex. The major rocks of the crater are composed of granulite, hornblendite, gneiss, tremolite marble, and marble. The bottom at the center of the crater covers about 100 m thick lacustrine sediments underlain by 188 m thick crater‐fill breccia. A layer of polymict breccia composed of clasts of granulite, gneiss, hornblendite, and fragments of glass as well as clastic matrix, occurs near the base, in the depth interval from 260 to 295 m. An investigation in quartz from the polymict breccia in the crater‐fill units reveals abundant planar deformation features (PDFs). Quartz with multiple sets of PDFs is found in clasts of granulite that consist of mainly quartz and feldspar, and in fine‐grained matrix of the impact‐produced polymict breccia. A universal stage was used to measure the orientation of PDFs in 70 grains of quartz from five thin sections made from the clasts of granulite of polymict breccia recovered at the depth of 290 m. Forty‐four percent of the quartz grains contain three sets of PDFs, and another 40% contain two sets of PDFs. The most abundant PDFs are rhombohedron forms of , , and with frequency of 33.5, 22.3, and 9.6%, respectively. A predominant PDF form of in quartz suggests a shock pressure >20 GPa. The occurrence of PDFs in quartz from the polymict breccia provides crucial evidence for shock metamorphism of target rocks and confirms the impact origin of this crater, which thus appears to be the first confirmed impact crater in China.     

Primordial,thermal, and shock features of ordinary chondrites: Emulating bulk X‐ray diffraction using in‐plane rotation of polished thin sections

Using an X‐ray diffractometer, powder‐like diffraction patterns were acquired from in‐plane rotation of polished thin sections (PTSs) of 60 ordinary chondrites (23 H, 21 L, and 16 LL), in order to explore the thermal and shock metamorphism and its modifications of primordial features. The olivine (Ol) 130 peak position shown as Bragg indices clearly correlates with the chemical group for equilibrated ordinary chondrites (EOCs), while the peak is split or broad for unequilibrated ordinary chondrites (UOCs). The intensity ratio of kamacite may be useful for distinguishing the chemical group between H and L‐LL, but it is not definite because of heterogeneous terrestrial weathering of kamacite, especially in H chondrites. The summed intensities of the orthoenstatite (Oen) 511 and 421 peaks positively correlates with the metamorphic sequence from 3 to 6, while that of clinoenstatite (Cen) 22 is inversely correlated. The shock stage positively correlates with the summed full width of half maximum values of the Oen 511 and 421 peaks and the FWHM of Ol 130 peak for each class. Significant amount of Oen (Pbca) transformed through Cen (C2/c) finally to Cen (P2₁/c) is stable at high pressure for shock stage S6 (Tenham and NWA 4719). The shock melted LL chondrite is characterized by the occurrence of Cen and abundant homogeneous olivine. The effects of both thermal and shock metamorphism are thus incorporated into the bulk X‐ray diffraction (XRD) data. The bulk XRD method is useful for determining the bulk mineralogy, resulting in the classification of ordinary chondrites. The method is also applicable to samples other than PTS.     

The sulfur species in hot rocky exoplanet atmospheres

The first JWST observations of hot Jupiters showed an unexpected detection of SO ${}_2$ in their hydrogen-rich atmospheres. We investigate how much sulfur can be expected in the atmospheres of rocky exoplanets and which sulfur molecules can be expected to be most abundant and detectable by transmission spectroscopy. We run thermochemical equilibrium models at the crust–atmosphere interface, considering surface temperatures 500–5000 K, surface pressures 1–100 bar, and various sets of element abundances based on common rock compositions. Between 1000 and 2000 K, we find gaseous sulfur concentrations of up to 25% above the rock in our models. SO ${}_2$, SO, H ${}_2$S, and S ${}_2$ are by far the most abundant sulfur molecules. SO ${}_2$ shows potentially detectable features in transmission spectra at about 4 $\mu$m, between 7 and 8 $\mu$m, and beyond 15 $\mu$m. In contrast, the sometimes abundant H ${}_2$S molecule is difficult to detect in these spectra, which are mostly dominated by H ${}_2$O and CO ${}_2$. Although the molecule PS only occurs with concentrations ppm, it can cause a strong absorption feature between 0.3 and 0.65 $\mu$m in some of our models for high surface pressures. The detection of sulfur molecules would enable a better characterization of the planetary surface.     

Shocked quartz grains in the polymict breccia of the Granby structure,Sweden–‐Verification of an impact

Abstract— The Middle Ordovician Granby structure in Sweden is generally considered the result of an asteroidal or cometary collision with Earth, although no hard evidence, i.e., shock metamorphic features or traces of the impactor, have been presented to date. In this study, drill core samples of a sedimentary breccia from the Granby structure have been investigated for microscopic shock metamorphic evidence in an attempt to verify the impact genesis of the structure. The finding of multiple sets of decorated planar deformation features (PDFs) in quartz grains in these samples provides unambiguous evidence that the structure is impact derived. Furthermore, the orientation of the PDFs, e.g., ω {101 }, π {101 } and r, z {101 }, is characteristic for impact deformation. The fact that a majority of the PDFs are decorated implies a water‐bearing target. The shocked quartz grains can be divided into two groups; rounded grains found in the breccia matrix likely originated from mature sandstone, and angular grains in fragments from crystalline target rocks. The absence of melt particles provides an estimated maximum shock pressure for the sedimentary derived quartz of 15–20 GPa and the frequency distribution of PDF orientations in the bedrock quartz implies pressures of the order of 10 GPa.     

Planar deformation features in quartz grains from the resurge deposit of the Lockne structure,Sweden

Abstract— Microscopic planar deformation features (PDFs) in quartz grains are diagnostic of shock meta-morphism during hypervelocity impact cratering. Measurements of the poles of sets of PDFs and the optic axis of 25 quartz grains were carried out for a sample of the Loftarsten deposit from the Lockne area, Sweden. The most abundant PDFs observed in the sample from the Lockne area correspond to those found at known impact craters (i.e., ω (1013} and π (1012). This study confirms the previous suggestion that the Lockne structure is an impact crater. The Loftarsten is, therefore, interpreted as the final stage of resurge deposition after a marine impact at Lockne in the Middle Ordovician.     

Investigating the partonic entanglement entropy at low Bjorken-x within Gluon saturation approach

In this work, the entanglement entropy is examined within the context of deep inelastic scattering in $ep$ collisions. The calculation is based on a formalism where the partonic state at small- $x$ is maximally entangled, consisting of a large number of micro-states occurring with equal probabilities. Analytical expressions for the number of gluons, $N_{\text{gluon}}$, are considered, derived from gluon saturation models for dipole-target amplitudes within the framework of the Quantum Chromodynamics (QCD) color dipole picture. A comparison of the entanglement entropy with thermodynamic entropy measured in $pp$ and $ep$ collisions at high energies is done.     

ANIE: A mathematical algorithm for automated indexing of planar deformation features in quartz grains

Abstract– Planar deformation features (PDFs) in quartz, one of the most commonly used diagnostic indicators of shock metamorphism, are planes of amorphous material that follow crystallographic orientations, and can thus be distinguished from non‐shock‐induced fractures in quartz. The process of indexing data for PDFs from universal‐stage measurements has traditionally been performed using a manual graphical method, a time‐consuming process in which errors can easily be introduced. A mathematical method and computer algorithm, which we call the Automated Numerical Index Executor (ANIE) program for indexing PDFs, was produced, and is presented here. The ANIE program is more accurate and faster than the manual graphical determination of Miller–Bravais indices, as it allows control of the exact error used in the calculation and removal of human error from the process.