Oceanic element inventories derived from marine sedimentary rocks place important constraints on oxidativecontinental weathering in deep time, but there remains a scarcity in complementary observations directly fromcontinental sedimentary reservoirs. This study focuses on better defining continental weathering conditions nearthe Archean-Proterozoic boundary through the multi-proxy (major and ultra-trace element, Fe and Cr stableisotopes, μ-XRF elemental mapping, and detrital zircon U-Pb geochronology) investigation of the ca. 2.45 billionyear old (giga annum, Ga) Cooper Lake paleosol (saprolith) developed on a sediment-hosted mafic dike withinthe Huronian Supergroup (Ontario, Canada).Throughout the variably altered Cooper Lake saprolith, ratios of immobile elements (Nb, Ta, Zr, Hf, Th, Al, Ti)are constant, indicating a uniform pre-alteration dike composition, lack of extreme pH weathering conditions,and no major influence from ligand-rich fluids during weathering or burial metasomatism/metamorphism. Theloss of Mg, Fe, Na, Sr, and Li, a signature of albite and ferromagnesian silicate weathering, increases towards thetop of the preserved profile (unconformity) and dike margins. Coupled bulk rock behaviour of Fe-Mg-Mn and colocalizationof Fe-Mn in clay minerals (predominantly chlorite) indicates these elements were solubilized primarilyin their divalent state without Fe/Mn-oxide formation. A lack of a Ce anomaly and immobility of Mo, V, and Cr further support pervasively anoxic weathering conditions. Subtle U enrichment, if primary, is the onlygeochemical evidence that could be consistent with oxidative element mobilization. The leaching of ferromagnesiansilicates was accompanied by variable mobility and depletion of transition metals with a relativedepletion order of Fe≈Mg≈Zn > Ni > Co > Cu (Cu being significantly influenced by secondary sulfideformation). Mild enrichment of heavy Fe isotopes (δ56/54Fe from 0.169 to 0.492‰) correlating with Fe depletionin the saprolith indicates open-system loss of isotopically light aqueous Fe(II). Minor REE+Y fractionation withincreasing alteration intensity, including a decreasing Eu anomaly and Y/Ho ratio, is attributed to albitebreakdown and preferential scavenging of HREE > Y by clay minerals, respectively. Younger metasomatismresulted in the addition of several elements (K, Rb, Cs, Be, Tl, Ba, Sn, In, W), partly or wholly obscuring theirearlier paleo-weathering trends.The behavior of Cr at Cooper Lake can help test previous hypotheses of an enhanced, low pH-driven continentalweathering flux of Cr(III) to marine reservoirs between ca. 2.48–2.32 Ga and the utility of the stable Crisotope proxy of Mn-oxide induced Cr(III) oxidation. Synchrotron μ-XRF maps and invariant Cr/Nb ratios revealcomplete immobility of Cr despite its distribution amongst both clay-rich groundmass and Fe-Ti oxides.Assuming a pH-dependent, continental source of Cr(III) to marine basins, the Cr immobility at Cooper Lakeindicates either that signatures of acidic surface waters were localized to uppermost and typically unpreservedregolith horizons or were geographically restricted to acid-generating point sources. However, given detritalpyrite preservation in overlying fluvial sequences, it is probable that the oxidative sulfide corrosion required todrive surface pH < 4 lagged behind in this region relative to other early Proterozoic sequences. The entiresaprolith exhibits a consistently light stable Cr isotope composition (δ53/52Cr: −0.321 ± 0.038‰, 2sd, n=34)that cannot be linked to Cr(III) oxidation and is instead interpreted to have a magmatic origin.

Shergottite Northwest Africa (NWA) 8159 is a basaltic rock derived from a mantle source with chemical characteristics that are unique in the martian meteorite suite. To further investigate this source reservoir, the Pb isotope compositions of plagioclase/maskelynite, pyroxene, phosphates, and shock melt-glass in NWA 8159 have been measured in situ by Secondary Ion Mass Spectrometry (SIMS). Due to the limited spread in Pb isotope data, these Pb isotope compositions have been used to calculate an imprecise PbPb isochron age of 3.4 ± 2.1 Ga (2σ), which is broadly consistent with the crystallization age of 2.37 ± 0.25 Ga determined previously by 147Sm143Nd. The lack of radiogenic in-growth within individual minerals since 2.4 Ga means that this sample is depleted in U, which is in agreement with NWA 8159's positive initial ε143Nd. An initial Pb composition was calculated using an x-y weighted average of the least radiogenic Pb isotope population measured in the sample. This initial Pb composition is not consistent with the model for Pb growth in the shergottite mantle at 2.4 Ga. This composition is, however, consistent with the model for the Nakhla-Chassigny mantle. Using the latter model, a source μ (238U/204Pb) of 2.6 ± 0.6 has been calculated. This μ-value is in contrast with the other depleted shergottites (1.4-1.5) and falls significantly off the array of source ε143Nd vs. μ defined by the rest of the martian meteorite suite and thus, necessitates a differentiation history distinct from the other martian meteorites. Sequestering Pb in sulphides during differentiation is the only mechanism to fractionate U from Pb and create a low-μ reservoir. Consequently, the relatively high μ-value of the source of NWA 8159 is in contrast with the positive initial ε143Nd and indicates that its mantle source region likely lacked significant sulphur. This is consistent with the lack of sulphides in the sample itself and could have played a role in its complicated oxidation history.

The triple oxygen isotope compositions of phosphate grains in six martian meteorites have been measured by Secondary Ion Mass Spectrometry (SIMS) and combined together with their chlorine isotope and halogen concentrations have been used to constrain hydrosphere-lithosphere interactions on Mars. These samples include three enriched shergottites (Zagami, Roberts Massif 04262 and Larkman Nunatak 12011), one depleted shergottite (Tissint), an orthopyroxenite (Allan Hills 84001), and a regolith breccia (Northwest Africa 7533). The phosphates measured here have a range in δ18O [(18O/16O)sample/(18O/16O)Standard-1] × 103] from +1.0 to +6.8‰ and could be a result of indigenous mantle values, mixing with martian water, or replacement reactions taking place on the surface of Mars. Three samples have a Δ17O [δ17O-1000(1 + δ18O /1000)0.528-1] in equilibrium with the martian mantle (ALH 84001, Tissint, and Zagami), while three samples (LAR 12011, RBT 04262, and NWA 7533) have an elevated positive Δ17O outside of analytical uncertainty of the martian fractionation line (MFL). The phosphates in the latter group also have positive and negative δ37Cl [(37Cl/35Cl)sample/(37Cl/35Cl)standard – 1] × 103] and enrichments in halogens not seen in the rest of the sample suite. Perchlorate formation on Earth fractionates Cl in both positive and negative directions and generates a correlated positive Δ17O. Further, perchlorate has been detected in wt% amounts on the martian surface. Thus, these results strongly suggest the presence of multiple Cl isotope reservoirs on the martian surface that have interacted with the samples studied here over the last ca. 2 Ga of geologic time. The weighted average of Δ17O measurements from phosphate grains (n = 13) in NWA 7533, which are the explicit result of exchange reactions on the martian surface, yields a statistically robust mean value of 1.39 ± 0.19‰ (2σ, MSWD = 1.5, p = 0.13). This value likely represents an accurate estimate for an oxidized surface reservoir on Mars.

Despite the intense cratering history of the Moon, very few traces of meteoritic material have been identified in the more than 380 kg of samples returned to Earth by the Apollo and Luna missions. Here we show that an ~200-µm-sized fragment collected by the Luna 16 mission has extra-lunar origins and probably originates from an LL chondrite with similar properties to near-Earth stony asteroids. The fragment has not experienced temperatures higher than 400 °C since its protolith formed early in the history of the Solar System. It arrived on the Moon, either as a micrometeorite or as the result of the break-up of a bigger impact, no earlier than 3.4 Gyr ago and possibly around 1 Gyr ago, an age that would be consistent with impact ages inferred from basaltic fragments in the Luna 16 sample and of a known dynamic upheaval in the Flora asteroid family, which is thought to be the source of L and LL chondrite meteorites. These results highlight the importance of extra-lunar fragments in constraining the impact history of the Earth–Moon system and suggest that material from LL chondrite asteroids may be an important component.

The ~31-km-wide Hiawatha structure, located beneath Hiawatha Glacier in northwestern Greenland, has been proposed as an impact structure that may have formed after the Pleistocene inception of the Greenland Ice Sheet. To date the structure, we conducted 40Ar/39Ar analyses on glaciofluvial sand and U-Pb analyses on zircon separated from glaciofluvial pebbles of impact melt rock, all sampled immediately downstream of Hiawatha Glacier. Unshocked zircon in the impact melt rocks dates to ~1915 million years (Ma), consistent with felsic intrusions found in local bedrock. The 40Ar/39Ar data indicate Late Paleocene resetting and shocked zircon dates to 57.99 ± 0.54 Ma, which we interpret as the impact age. Consequently, the Hiawatha impact structure far predates Pleistocene glaciation and is unrelated to either the Paleocene-Eocene Thermal Maximum or flood basalt volcanism in east Greenland. However, it was contemporaneous with the Paleocene Carbon Isotope Maximum, although the impact’s exact paleoenvironmental and climatic significance awaits further investigation.

We report on the first zircon hafnium-oxygen isotope and trace element study of a transect through one of the largest terrestrial impact melt sheets. The differentiated melt sheet at the 1.85 Ga, originally ca. 200 km in diameter Sudbury impact crater, Ontario, Canada, yields a tight range of uniform zircon Hf isotope compositions (ε_Hf(1850) of ca. −9 to −12). This is consistent with its well-established crustal origin and indicates differentiation from a single melt that was initially efficiently homogenised. We propose that the heterogeneity in other isotopic systems, such as Pb, in early-emplaced impact melt at Sudbury is associated with volatility-related depletion during the impact cratering process. This depletion leaves the isotopic systems of more volatile elements more susceptible to contamination during post-impact assimilation of country rock, whereas the systems of more refractory elements preserve initial homogeneities. Zircon oxygen isotope compositions in the melt sheet are also restricted in range relative to those in the impacted target rocks. However, they display a marked offset approximately one-third up the melt sheet stratigraphy that is interpreted to be a result of post-impact assimilation of ¹⁸O-enirched rocks into the base of the cooling impact melt.

Given that impact cratering was a more dominant process in the early history of the inner Solar System than it is today, and the possibility that impact melt sheets were sources of ex situ Hadean zircon grains, these findings may have significance for the interpretation of the early zircon Hf record. We speculate that apparent ε_Hf-time arrays observed in the oldest terrestrial and lunar zircon datasets may be related to impact melting homogenising previously more diverse crust.

We also show that spatially restricted partial melting of rocks buried beneath the superheated impact melt at Sudbury provided a zircon crystallising environment distinct to the impact melt sheet itself.

Impact cratering was a fundamental geological process in the early Solar System and, thus, constraining the timescales over which large impact structures cool is critical to understanding the thermal evolution and habitability of early planetary crusts. Additionally, impacts can induce mass extinctions and establishing the precise timing of the largest impacts on Earth can shed light on their role in such events. Here we report a high-precision zircon U–Pb geochronology study of the Morokweng impact structure, South Africa, which appears to have a maximum present-day diameter of ∼80 km. Our work provides (i) constraints on the cooling of large impact melt sheets, and (ii) a high-precision age for one of Earth's largest impact events, previously proposed to have overlapped the ca. 145 Ma Jurassic–Cretaceous (J–K) boundary. High-precision U–Pb geochronology was performed on unshocked, melt-grown zircon from five samples from a borehole through approximately 800 m of preserved impact melt rock. Weighted mean 206Pb/238U dates for the upper four samples are indistinguishable, with relative uncertainties (internal errors) of better than 20 ka, whereas the lowermost sample is distinguishably younger than the others. Thermal modeling suggests that the four indistinguishable dates are consistent with in situ conductive cooling of melt at this location within 30 kyr of the impact. The younger date from the lowest sample cannot be explained by in situ conductive cooling in line with the overlying samples, but the date is within the ∼65 kyr timeframe for melt-present conditions in footwall rocks below the impact melt sheet that is indicated by our thermal model. The Morokweng impact event is here constrained to 146.06 ± 0.16 Ma (2σ; full external uncertainty), which precedes current estimates of the age of the J–K boundary by several million years.

Despite the wide utility of apatite, Ca5(PO4)3(F,Cl,OH), in the geosciences, including tracing volatile abundances on the Moon and Mars, little is known about how the mineral responds to the extreme temperatures and pressures associated with hypervelocity impacts. To address this deficiency, we here present the first microstructural analysis and chemical mapping of shocked apatite from a terrestrial impact crater. Apatite grains from the Paasselkä impact structure, Finland, display intragrain crystal-plastic deformation as well as pervasive recrystallization—the first such report in terrestrial apatite. A partially recrystallized grain offers the opportunity to investigate the effect of shock recrystallization on the chemical composition of apatite. The recrystallized portion of the fluorapatite grain is depleted in Mg and Fe relative to the remnant non-recrystallized domain. Strikingly, the recrystallized region alone hosts inclusions of (Mg,Fe)2(PO4)F, wagnerite or a polymorph thereof. These are interpreted to be a product of phase separation during recrystallization and to be related to the reduced abundances of certain elements in the recrystallized domain. The shock-induced recrystallization of apatite, which we show to be related to changes in the mineral’s chemical composition, is not always readily visible in traditional imaging techniques (such as backscattered electron imaging of polished interior surfaces), thus highlighting the need for correlated microstructural, chemical, and isotopic studies of phosphates. This is particularly relevant for extraterrestrial phosphates that may have been exposed to impacts, and we urge the consideration of microstructural data in the interpretation of the primary or secondary nature of elemental abundances and isotopic compositions.

Uranium-lead (U-Pb) geochronology of individual shocked zircon grains has unique potential for dating bolide impact events. Neoblasts in granular-textured zircon have been recognized as the shock-related feature most effective at recording the impact age. Here we report the discovery of large neoblasts (5–100 µm in dimension) in shocked zircon at the Sudbury impact structure, Canada—the first report of in situ coarsely granular zircon from a terrestrial impact site other than the Vredefort structure, South Africa. The neoblast-bearing sample was taken from a heterogeneous, lithic clast–rich igneous unit associated with the roof rocks of the impact melt sheet, making this the first time a crater has been dated using neoblastic zircon from the upper part of its stratigraphy. Previous in situ discoveries of coarsely granular zircon at Vredefort were all in impact-generated mafic melt emplaced beneath the impact melt sheet. Electron backscatter diffraction analysis of the impact-aged neoblasts indicates that the high-pressure conditions inferred in the formation of many small neoblasts were not necessarily involved in the formation of these large ones. Their large size, internal zonation, and occurrence in a slowly cooling environment collectively suggest that large neoblasts at Sudbury formed by relatively protracted, post-impact growth in shocked zircon incorporated into impact-related melt. Based on insight from large neoblast growth in terrestrial settings, we suggest that the ca. 4.33 Ga neoblasts recently reported in lunar zircon may imply a major basin-forming event on the Moon at that time. New knowledge of the cratering environments in which large neoblasts form also raises the prospect of possibly linking ex situ granular zircon in lunar breccias with specific impact structures—and thus better calibrating the lunar cratering record with radiometric ages.

We report a new age for the Sääksjärvi impact structure, Finland, a 6 km diameter feature that formed in crystalline rocks of the Precambrian Baltic Shield. Two previous studies reported 40Ar/39Ar data for Sääksjärvi and suggested conflicting formation ages of ≤330 Ma or c. 560 Ma. The former age represents a possible complication for models which indicate that the region was covered by sediments of the Caledonian foreland basin throughout much of the Phanerozoic. We conducted a study combining imaging, microstructural analysis and U–Pb dating of shocked zircon from Sääksjärvi. The U–Pb dataset indicates a c. 600 Ma impact into predominantly c. 1850 Ma target rocks. A concordia age of 608 ± 8 Ma (2σ) confirms Sääksjärvi as the first known Ediacaran impact structure in the Baltic Shield and only the second worldwide. Our data indicate that the Sääksjärvi impact structure formed in exposed crystalline basement rocks of the Baltic Shield prior to the development of the Caledonian foreland basin. Given that most impact structures on Earth are relatively small features, radiometric dating of small impact structures in crystalline basement may place boundaries on the timing and spatial extent of palaeobasins that might otherwise be difficult to constrain.

Hypervelocity impacts can produce features in zircon that are not normally produced by endogenic processes. However, lightning can also induce extreme pressure–temperature excursions, and itsefect on zircon has not been studied. With the aim to recognise features that form in response to extreme pressure–temperature excursions but are not unique to hypervelocity impacts, we imaged and undertook microstructural characterization of zircon in a fulgurite (a tubular body of glass andfused clasts that formed in response to a lightning strike). We document zircon with granular ZrO2 and rims of vermicular ZrO2, features which vary in abundance with increasing distance from the fulgurite’s central void. This indicates that these features formed in response to the lightning strike. Zircon dissociation to ZrO2 and SiO2 is a high-temperature, relatively low-pressure phenomenon, consistent with previous suggestions that lightning strikes involve extreme temperatures as well as pressures greater than those usually generated in Earth’s crust but rarely > 10 GPa. The rims of monoclinic ZrO2 record crystallographic evidence for precursor cubic ZrO2, demonstrating that cubic ZrO2 is not unique to hypervelocity impacts. Given the likelihood that this fulgurite experienced pressures of, at most, a few GPa, evidence for cubic ZrO2 indicates peak temperatures > 2000 °C.

Constraining the origin and history of very ancient detrital zircons has unique potential for furthering our knowledge of Earth’s very early crust and Hadean geodynamics. Previous applications of the Ti-in-zircon thermometer to >4 Ga zircons have identified a population with relatively low crystallization temperatures () of ~685 °C. This could possibly indicate wet minimum-melting conditions producing granitic melts, implying very different Hadean terrestrial geology from other rocky planets. Here we report the first comprehensive ion microprobe study of zircons from a transect through the differentiated Sudbury impact melt sheet. The new zircon Ti contents and corresponding  fully overlap with those of the Hadean zircon population. Previous studies, which measured Ti in impact melt sheet zircons did not find this wide range because they analyzed samples only from a restricted portion of the melt sheet and because they used laser ablation analyses that can overestimate true Ti content. It is important to note that internal differentiation of the impact melt is likely a prerequisite for the observed low  in zircons from the most evolved rocks. On Earth, melt sheet

Granular neoblastic zircon (ZrSiO4) with systematically oriented granules has been proposed as evidence for extreme shock pressures (>30 GPa) and subsequent high temperatures (>1200 ◦C). It is widely agreed to reflect the solid-state phase transition from zircon to its high-pressure polymorph reidite and subsequent reversion to zircon. This model is based on crystallographic relationships between granules of a single type of granular zircon and does not explain the formation of other types of granular zircon textures, for example, grains with randomly oriented granules or with large, often euhedral granules. Here we report the first nano-scale observations of granular neoblastic zircon and the surrounding environment. We conducted combined microstructural analyses of zircon in the lithic clast from an impact melt dike of the Vredefort impact structure. Zircon granules have either random or systematic orientation with three mutually orthogonal directions of their c-axes coincident with [110] axes. Each1-2 μm zircon granule is a mosaic crystal composed of nanocrystalline subunits. Granules contain round inclusions of baddeleyite (monoclinic ZrO2) and amorphous silica melt. Tetragonal and cubic ZrO2 also occur as sub-μm-sized inclusions (<50 nm). Filament-like aggregates of nanocrystalline zircon are present as “floating” in the surrounding silicate matrix. They are aligned with each other, apparently serving as the building blocks for the mosaic zircon crystals (granules). Our results indicate shock-related complete melting of zircon with the formation of immiscible silicate and oxide melts. The melts reacted and crystallized rapidly as zircon granules, some of which experienced growth alignment/twinning and parallel growth, causing the characteristic systematic orientation of the granules observed for some ofthe aggregates. In contrast to the existing model, in which this type of granular zircon is considered to be a product of reversion from the high-pressure polymorph reidite, our nano-scale observations suggesta formation mechanism that does not require phase transition via reidite but is indicative of instant incongruent decomposition, melting and rapid crystallization from the melt

Abstract Constraining the duration of magmatic activity on the Moon is essential to understand how the lunar mantle evolved chemically through time. Determining age and initial isotopic compositions of mafic lunar meteorites is a critical step in defining the periods of magmatic activity that occurred during the history of the Moon and to constrain the chemical characteristics of mantle components involved in the sources of the magmas. We have used the in situ Pb-Pb SIMS technique to investigate eight lunar gabbros and basalts, including six meteorites from the Northwest Africa (NWA) 773 clan (NWA 2727, NWA 2700, NWA 3333, NWA 2977, NWA 773, and NWA 3170), NWA 4734, and Dhofar 287A. These samples have been selected as there is no clear agreement on their age and they are all from the dominant low titanium chemical group. We have obtained ages of 2981 ± 12 Ma for NWA 4734 and 3208 ± 22 Ma for Dhofar 287. For the NWA 773 clan, four samples (the fine-grained basalt NWA 2727 and the three gabbros NWA 773, NWA 2977, NWA 3170) out of six yielded isochron-calculated ages that are identical within uncertainties and yielding an average age of 3086 ± 5 Ma. The age obtained for the fine-grained basalt NWA 2700 is not precise enough for comparison with the other samples. The gabbroic sample NWA 3333 yielded an age of 3038 ± 20 Ma suggesting that two distinct magmatic events may be recorded in the meteorites of the NWA 773 clan. The present study aims to identify and assess all potential issues that are associated with different ways to date lunar rocks using U-Pb?based methods. To achieve this, we have compared the new ages with the previously published data set. The entire age data set from lunar mafic meteorites was also screened to identify data showing analytical issues and evidence of resetting and terrestrial contamination. The data set combining the ages of mafic lunar meteorites and Apollo rocks suggests pulses of magmatic activity with two distinct phases between 3950 and 3575 Ma and between 3375 and 3075 Ma with the two phases separated by a gap of approximately 200 Ma. The evolution of the Pb initial ratios of the low-Ti mare basalts between approximately 3400 and 3100 Ma suggests that these rocks were progressively contaminated by a KREEP-like component.

Apatite's ubiquity in crystalline rocks, variable trace element contents (particularly with regard to the REE, actinides and Sr), and amenability to various datingtechniques based on the decay of the radioisotopes U and Th, permit specific provenance determinations. In this study, we first present a comprehensive descriptionof the trace element behaviour of apatite in various kinds of bedrocks (igneous rocks from felsic through to ultramafic compositions, metamorphic rocks from low tohigh grades and of diverse protolith composition, and authigenic apatite) in which we explain why apatite is so highly diverse in terms of its trace elementcomposition. Next, we present a synthesis of bedrock apatite trace-element compositional data from previous work, assembling a library of apatite compositions thatincludes the most abundant apatite-bearing lithologies in the Earth's crust, and many other less abundant rock types. Compositional statistics, classification, and amachine learning classifier are then applied to this dataset to generate biplots that can be used to determine the broad source lithology of detrital apatite, withmisclassification averaging 15%. This methodology is tested in three case studies to demonstrate its utility. In these examples, detrital apatite can be convincinglylinked to different lithology types, and combined apatite trace-element and UePb data can determine the terranes from which individual apatites were likely derived.The addition of apatite trace-element information therefore enables the determination of the source lithology, making the extraction of novel information and morespecific provenance determinations possible, and opening up new avenues in source-to-sink modelling.

Reidite, the high-pressure zircon (ZrSiO4) polymorph, is a diagnostic indicator of impact events. Natural records of reidite are, however, scarce, occurring mainly as micrometer-sized lamellae, granules, and dendrites. Here, we present a unique sequence of shocked zircon grains found within a clast from the Chassenon suevitic breccia (shock stageIII) from the ~200 Ma, 20–50 km wide Rochechouart impact structure in France. Our study comprises detailed characterization with scanning electron microscopy coupled with electron backscatter diffraction with the goal of investigating the stability and response of ZrSiO4 under extreme P–T conditions. The shocked zircon grains have preserved various amounts of reidite ranging from 4% up to complete conversion. The grains contain various variants of reidite, including the common habits: lamellae and granular reidite. In addition, three novel variants have been identified: blade, wedge, and massive domains. Several of these crosscut and offset each other, revealing that reidite can form at multiple stages during an impact event. Our data provide evidence that reidite can be preserved in impactites to a much greater extent than previously documented. We have further characterized reversion products of reidite in the form of fully recrystallized granular zircon grains and minute domains of granular zircon in reidite-bearing grains that occur in close relationship to reidite. Neoblasts in these grains have a distinct crystallography that is the result of systematic inheritance of reidite. We interpret that the fully granular grains have formed from prolonged exposure of temperatures in excess of 1200 °C. Reidite-bearing grains with granular domains might signify swift quenching from temperatures close to 1200 °C. Grains subjected to these specific conditions therefore underwent partial zircon-to-reidite reversion, instead of full grain recrystallization. Based on our ZrSiO4 microstructural constraints, we decipher the grains evolution at specific P–T conditions related to different impact stages, offering further understanding of the behavior of ZrSiO4 during shock.

Lunar meteorites provide a potential opportunity to expand the study of ancient (>4000 Ma) basaltic volcanism on the Moon, of which there are only a few examples in the Apollo sample collection. Secondary Ion Mass Spectrometry (SIMS) was used to determine the Pb isotopic compositions of multiple mineral phases (Ca-phosphates, baddeleyite K-feldspar, K-rich glass and plagioclase) in two lunar meteorites, Miller Range (MIL) 13317 and Kalahari (Kal) 009. These data were used to calculate crystallisation ages of 4332 ± 2 Ma (95% confidence level) for basaltic clasts in MIL 13317, and 4369 ± 7 Ma (95% confidence level) for the monomict basaltic breccia Kal 009. From the analyses of the MIL 13317 basaltic clasts, it was possible to determine an initial Pb isotopic composition of the protolith from which the clasts originated, and infer a 238 U/204 Pb ratio (μ-value) of 850 ± 130 (2σ uncertainty) for the magmatic source of this basalt. This is lower than μ-values determined previously for KREEP-rich (an acronym for K, Rare Earth Elements and P) basalts, although analyses of other lithological components in the meteorite suggest the presence of a KREEP component in the regolith from which the breccia was formed and, therefore, a more probable origin for the meteorite on the lunar nearside. It was not possible to determine a similar initial Pb isotopic composition from the Kal 009 data, but previous studies of the meteorite have highlighted the very low concentrations of incompatible trace elements and proposed an origin on the farside of the Moon. Taken together, the data from these two meteorites provide more compelling evidence for widespread ancient volcanism on the Moon. Furthermore, the compositional differences between the basaltic materials in the meteorites provide evidence that this volcanism was not an isolated or localised occurrence, but happened in multiple locations on the Moon and at distinct times. In light of previous studies into early lunar magmatic evolution, these data also imply that basaltic volcanism commenced almost immediately after Lunar Magma Ocean (LMO) crystallisation, as defined by Nd, Hf and Pb model ages at about 4370 Ma.