Time scales of magma transport and mixing at Kīlauea Volcano, Hawai'i

Modelling of volcanic processes is strongly limited by a poor knowledge of the timescales of storage, mixing and final ascent of magmas into the shallowest portions of volcanic 'plumbing' systems immediately prior to eruption. It is impossible to measure these timescales directly; however, micro-analytical techniques provide indirect estimates based on the extent of diffusion of species through melts and crystals. Here, diffusion in olivine phenocrysts from the 1959 Kīlauea Iki eruption is used to constrain the timing of mixing events in the crustal plumbing system on timescales of months to years before eruption. The timescales derived from zonation of Fe-Mg in olivines, combined with contemporaneous geophysical data suggests mixing occurred on 3 timescales: (1) up to 2 years prior to eruption in the deep storage system, (2), in a shallow reservoir, between incoming hot melts and cooler, resident melt for several weeks to months prior to eruption, and (3), in the conduit and summit reservoir, between the resident magma and cooled surface lava, draining back into the vent on timescales of hours to several days during pauses between episodes. Synchronous inflation of the shallow reservoir with deep earthquake swarms and mixing suggests a fitfully open transcrustal magmatic system prior to and during eruption.We acknowledge NERC studentship funds (I. Sides) and a United States Geological Survey Jack Kleinman grant, which allowed samples for this study to be collected.This is the final version of the article. It first appeared from the Geological Society of America via https://doi.org/10.1130/G37800.

Ancient and recent collisions revealed by phosphate minerals in the Chelyabinsk meteorite

AbstractThe collision history of asteroids is an important archive of inner Solar System evolution. Evidence for these collisions is brought to Earth by meteorites. However, as meteorites often preserve numerous impact-reset mineral ages, interpretation of their collision histories is controversial. Here, we combine analysis of phosphate U-Pb ages and microtextures to interpret the collision history of Chelyabinsk—a highly shocked meteorite. We show that phosphate U-Pb ages correlate with phosphate microtextural state. Pristine phosphate domain U-Pb compositions are generally concordant, whereas fracture-damaged domains universally display discordance. Combining both populations best constrains upper (4473 ± 11 Ma) and lower intercept (−9 ± 55 Ma, i.e., within error of present) U-Pb ages. All phosphate U-Pb ages were completely reset during an ancient high energy collision, whilst fracture-damaged domains experienced further Pb-loss during mild and recent collisional re-heating. Targeting textural sub-populations of phosphate grains permits more robust reconstruction of asteroidal collision histories.</jats:p

Dynamic Split Tensile Strength of Basalt, Granite, Marble and Sandstone: Strain Rate Dependency and Fragmentation

Funder: Albert-Ludwigs-Universität Freiburg im Breisgau (1016)AbstractThe aim of this study is to understand the strength behaviour and fragment size of rocks during indirect, quasi-static and dynamic tensile tests. Four rocks with different lithological characteristics, namely: basalt, granite, sandstone, and marble were selected for this study. Brazilian disc experiments were performed over a range of strain rates from ~ 10–5 /s to 2.7 × 101 /s using a hydraulic loading frame and a split Hopkinson bar. Over the range of strain rates, our measurements of dynamic strength increase are in good agreement with the universal theoretical scaling relationship of (Kimberley et al., Acta Mater 61:3509–3521, 2013). Dynamic fragmentation during split tension mode failure has received little attention, and in the present study, we determine the fragment size distribution based on the experimentally fragmented specimens. The fragments fall into two distinct groups based on the nature of failure: coarser primary fragments, and finer secondary fragments. The degree of fragmentation is assessed in terms of characteristic strain rate and is compared with existing theoretical tensile fragmentation models. The average size of the secondary fragments has a strong strain rate dependency over the entire testing range, while the primary fragment size is less sensitive at lower strain rates. Marble and sandstone are found to generate more pulverised secondary debris when compared to basalt and granite. Furthermore, the mean fragment sizes of primary and secondary fragments are well described by a power-law function of strain rate.</jats:p

Can Archean Impact Structures Be Discovered? A Case Study From Earth's Largest, Most Deeply Eroded Impact Structure

Funder: Trinity College, University of Cambridge; doi: http://dx.doi.org/10.13039/501100000727AbstractThe record of terrestrial impact events is incomplete with no Archean impact structures discovered, despite the expected abundance of collisions that must have occurred. Because no Archean impact structures have been identified, the necessary conditions to preserve an impact structure longer than 2 Byr are unknown. One significant effect of shock metamorphism is that the physical properties of the target rocks change, resulting in distinctive geophysical signatures of impact structures. To evaluate the preservation potential of impact structures, we evaluate the deeply eroded Proterozoic Vredefort impact structure to examine the changes in physical properties and the remnant of the geophysical signature and compare the results with the well‐preserved Chicxulub impact structure. The major structural features of Vredefort are similar to the expected profile of the Chicxulub structure at a depth of 8–10 km. The Vredefort target rocks, while shocked, do not preserve measurable changes in their physical properties. The gravity signature of the impact structure is minor and is controlled by the remnant of the collapsed transient crater rim and the uplifted Moho surface. We anticipate that erosion of the Vredefort structure by an additional 1 km would remove evidence of impact, and regardless of initial size, erosion by &gt;10 km would result in the removal of most of the evidence for any impact structure from the geological record. This study demonstrates that the identification of geologically old (i.e., Archean) impact structures is limited by a lack of geophysical signatures associated with deeply eroded craters.</jats:p

The Distribution of Impactor Core Material During Large Impacts on Earth-like Planets

Abstract Large impacts onto young rocky planets may transform their compositions, creating highly reducing conditions at their surfaces and reintroducing highly siderophile metals to their mantles. Key to these processes is the availability of an impactor’s chemically reduced core material (metallic iron). It is, therefore, important to constrain how much of an impactor’s core remains accessible to a planet’s mantle/surface, how much is sequestered to its core, and how much escapes. Here, we present 3D simulations of such impact scenarios using the shock physics code iSALE to determine the fate of impactor iron. iSALE’s inclusion of material strength is vital in capturing the behavior of both solid and fluid components of the planet and thus characterizing iron sequestration to the core. We find that the mass fractions of impactor core material that accretes to the planet core (f core) or escapes (f esc) can be readily parameterized as a function of a modified specific impact energy, with f core &gt; f esc for a wide set of impacts. These results differ from previous works that do not incorporate material strength. Our work shows that large impacts can place substantial reducing impactor core material in the mantles of young rocky planets. Impact-generated reducing atmospheres may thus be common for such worlds. However, through escape and sequestration to a planet’s core, large fractions of an impactor’s core can be geochemically hidden from a planet’s mantle. Consequently, geochemical estimates of late bombardments of planets based on mantle siderophile element abundances may be underestimates.</jats:p

Reduced atmospheres of post-impact worlds: The early Earth

Impacts may have had a significant effect on the atmospheric chemistry of the early Earth. Reduced phases in the impactor (e.g., metallic iron) can reduce the planet's H$_2$O inventory to produce massive atmospheres rich in H$_2$. Whilst previous studies have focused on the interactions between the impactor and atmosphere in such scenarios, we investigate two further effects, 1) the distribution of the impactor's iron inventory during impact between the target interior, target atmosphere, and escaping the target, and 2) interactions between the post-impact atmosphere and the impact-generated melt phase. We find that these two effects can potentially counterbalance each other, with the melt-atmosphere interactions acting to restore reducing power to the atmosphere that was initially accreted by the melt phase. For a $\sim10^{22}\,\mathrm{kg}$ impactor, when the iron accreted by the melt phase is fully available to reduce this melt, we find an equilibrium atmosphere with H$_2$ column density $\sim10^4\,\mathrm{moles\,cm^{-2}}$ ($p\mathrm{H2}\sim120\,\mathrm{bars}\mathrm{,}~X_\mathrm{H2}\sim0.77$), consistent with previous estimates. However, when the iron is not available to reduce the melt (e.g., sinking out in large diameter blobs), we find significantly less H$_2$ ($7\times10^2-5\times10^3\,\mathrm{moles\,cm^{-2}}$, $p\mathrm{H2}\lesssim60\,\mathrm{bars}\mathrm{,}~X_\mathrm{H2}\lesssim0.41$). These lower H$_2$ abundances are sufficiently high that species important to prebiotic chemistry can form (e.g., NH3, HCN), but sufficiently low that the greenhouse heating effects associated with highly reducing atmospheres, which are problematic to such chemistry, are suppressed. The manner in which iron is accreted by the impact-generated melt phase is critical in determining the reducing power of the atmosphere and re-solidified melt pool in the aftermath of impact

Orientations of planar cataclasite zones in the Chicxulub peak ring as a ground truth for peak ring formation models

Hypervelocity impact cratering is an important geologic process but the rarity of large terrestrial impact craters on Earth and the limited technical options to study cratering processes in the laboratory hinders our understanding of large-scale impact processes. Drill core recovered from the peak ring of the Chicxulub impact crater during International Ocean Discovery Program (IODP)/International Continental scientific Drilling Program (ICDP) Expedition 364 provides an opportunity to examine target rock deformation and thus, to assess cratering models in this regard. Using oriented computer tomography (CT) scans and line scan images of the core, we present the orientations of mm-to-cm-scale planar cataclasite and ultracataclasite zones in the deformed granitoid target rock of the peak ring. In the upper 470 m of the target rock, the cataclasite zones dip towards the crater center, whereas the dip directions for the ultracataclasite zones are approximately tangential to the peak ring. These two orientations are consistent with deformation expected from hydrocode-modeled principal stress directions for the outward movement of rocks as the transient crater develops, and the inward movement of rocks associated with collapse of the transient crater. Near the base of the core is a 96 m-thick interval of highly-deformed target rock with impact melt rock and rock fragments not observed elsewhere in the core; this interval has previously been interpreted as an imbricate thrust zone within the peak ring. The cataclasite zones below this thrust zone have different orientations than those in the 470 m-thick block above. This observation implies a differential rotation from the overlying block during the final stages of peak-ring formation. Our results support an acoustic fluidization process, wherein blocks that vibrate or slide relative to each other allow the target rock to flow during transient crater collapse, and that the size of coherent rock blocks increases over the course of crater modification as the target rock regains its cohesive strength and acoustic fluidization decreases

New shock microstructures in titanite (CaTiSiO5) from the peak ring of the Chicxulub impact structure, Mexico

Accessory mineral geochronometers such as apatite, baddeleyite, monazite, xenotime and zircon are increasingly being recognized for their ability to preserve diagnostic microstructural evidence of hypervelocity-impact processes. To date, little is known about the response of titanite to shock metamorphism, even though it is a widespread accessory phase and a U–Pb geochronometer. Here we report two new mechanical twin modes in titanite within shocked granitoid from the Chicxulub impact structure, Mexico. Titanite grains in the newly acquired core from the International Ocean Discovery Program Hole M0077A preserve multiple sets of polysynthetic twins, most commonly with composition planes (K1) = ~  {1¯11} , and shear direction (η1) =  , and less commonly with the mode K1 = {130}, η1 = ~  . In some grains, {130} deformation bands have formed concurrently with the deformation twins, indicating dislocation slip with Burgers vector b =  can be active during impact metamorphism. Titanite twins in the modes described here have not been reported from endogenically deformed rocks; we, therefore, propose this newly identified twin form as a result of shock deformation. Formation conditions of the twins have not been experimentally calibrated, and are here empirically constrained by the presence of planar deformation features in quartz (12 ± 5 and ~ 17 ± 5 GPa) and the absence of shock twins in zircon (< 20 GPa). While the lower threshold of titanite twin formation remains poorly constrained, identification of these twins highlight the utility of titanite as a shock indicator over the pressure range between 12 and 17 GPa. Given the challenges to find diagnostic indicators of shock metamorphism to identify both ancient and recent impact evidence on Earth, microstructural analysis of titanite is here demonstrated to provide a new tool for recognizing impact deformation in rocks where other impact evidence may be erased, altered, or did not manifest due to generally low (< 20 GPa) shock pressure

Rock fluidization during peak-ring formation of large impact structures

Large meteorite impact structures on the terrestrial bodies of the Solar System contain pronounced topographic rings, which emerged from uplifted target (crustal) rocks within minutes of impact. To flow rapidly over large distances, these target rocks must have weakened drastically, but they subsequently regained sufficient strength to build and sustain topographic rings. The mechanisms of rock deformation that accomplish such extreme change in mechanical behaviour during cratering are largely unknown and have been debated for decades. Recent drilling of the approximately 200-km-diameter Chicxulub impact structure in Mexico has produced a record of brittle and viscous deformation within its peak-ring rocks. Here we show how catastrophic rock weakening upon impact is followed by an increase in rock strength that culminated in the formation of the peak ring during cratering. The observations point to quasi-continuous rock flow and hence acoustic fluidization as the dominant physical process controlling initial cratering, followed by increasingly localized faulting