# Early crustal processes revealed by the ejection site of the oldest martian meteorite

A. Lagain, Sylvain Bouley, B. Zanda, K. Miljković, A. Rajšić, David

Baratoux, V. Payré, L. Doucet, N. Timms, R. Hewins, et al.

► To cite this version:

A. Lagain, Sylvain Bouley, B. Zanda, K. Miljković, A. Rajšić, et al. Early crustal processes revealed by the ejection site of the oldest martian meteorite. Nature Communications, 2022, 13, pp.3782. 10.1038/s41467-022-31444-8. hal-03740195

HAL Id: hal-03740195

https://hal.science/hal-03740195v1

Submitted on 30 Jul 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

ARTICLE

https://doi.org/10.1038/s41467-022-31444-8

OPEN

<figure><img src="image_3.png" /></figure>

Check for updates

# Early crustal processes revealed by the ejection site of the oldest martian meteorite

A. Lagain ID <fnref n="1" />⁻³, S. Bouley²,³, B. Zanda³,⁴, K. Miljković ID <fnref n="1" />, A. Rajšič<fnref n="1" />, D. Baratoux⁵,⁶, V. Payré<fnref n="7" />, L. S. Doucet ID ⁸, N. E. Timms<fnref n="1" />,⁹, R. Hewins ID ⁴,¹⁰, G. K. Benedix ID <fnref n="1" />,¹¹,¹², V. Malarewic²,⁴, K. Servis ID <fnref n="1" />,¹³ & P. A. Bland<fnref n="1" />

The formation and differentiation of the crust of Mars in the first tens of millions of years after its accretion can only be deciphered from incredibly limited records. The martian breccia NWA 7034 and its paired stones is one of them. This meteorite contains the oldest martian igneous material ever dated: ~4.5 Ga old. However, its source and geological context have so far remained unknown. Here, we show that the meteorite was ejected 5-10 Ma ago from the north-east of the Terra Cimmeria—Sirenum province, in the southern hemisphere of Mars. More specifically, the breccia belongs to the ejecta deposits of the Khujirt crater formed 1.5 Ga ago, and it was ejected as a result of the formation of the Karratha crater 5-10 Ma ago. Our findings demonstrate that the Terra Cimmeria—Sirenum province is a relic of the differentiated primordial martian crust, formed shortly after the accretion of the planet, and that it constitutes a unique record of early crustal processes. This province is an ideal landing site for future missions aiming to unravel the first tens of millions of years of the history of Mars and, by extension, of all terrestrial planets, including the Earth.

<footnote n="1">¹Space Science and Technology Centre, School of Earth and Planetary Science, Curtin University, Perth, WA, Australia. ²Université Paris-Saclay, CNRS, GEOPS, 91405 Orsay, France. ³IMCCE, Observatoire de Paris, 77 avenue Denfert-Rochereau, 75005 Paris, France. ⁴Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie (IMPMC), Muséum national d'Histoire naturelle, Sorbonne Université et CNRS, 75005 Paris, France. ⁵Géosciences Environnement Toulouse, University of Toulouse, CNRS and IRD, Toulouse 31400, France. ⁶Université Félix Houphouët-Boigny, Abidjan, Côte d'Ivoire.</footnote>

<footnote n="7">⁷Department of Astronomy and Planetary Science, Northern Arizona University, Flagstaff, AZ, USA. ⁸Earth Dynamics Research Group, TIGeR, School of Earth and Planetary Sciences, Curtin University, Perth, WA, Australia. ⁹The Institute for Geoscience Research (TIGeR), Curtin University, Perth 6845 WA, Australia. ¹⁰EPS, Rutgers University, Piscataway, NJ 08854, USA. ¹¹Department of Earth and Planetary Sciences, Western Australian Museum, Perth, WA, Australia. ¹²Planetary Sciences Institute, Tucson, AZ, USA. ¹³Pawsey Supercomputing Centre, CSIRO, Kensington, WA, Australia. ¹⁴email: anthony.lagain@gmail.com</footnote>
The geological record of the formation and differentiation of our planet has been destroyed by its subsequent evolution, but extremely rare clues may be obtained from other terrestrial planets. Mars provides a unique and accessible example of an early evolutionary path corresponding to that, inaccessible, of our own world. We can investigate it with spacecraft, and samples are available for in-depth analysis on Earth in the form of martian meteorites. So far, the only available martian samples that appear to have recorded the early conditions and the evolution of the planet until the present time are Northwest Africa (NWA) 7034 and its paired stones. They are the most diverse martian meteorites in terms of composition, containing a variety of igneous, sedimentary, and impact melt clasts, including the most evolved and oldest igneous clasts and zircons (4.44–4.48 Ga old¹⁻¹⁰, grey bars in Fig. 1). These evolved clasts are derived from a variety of magmas (monzonitic or mugearitic) and probably formed by re-melting of the primary martian crust either at various depths in the presence of volatiles or by differentiation of large impact melt sheets⁹,¹⁰. The abundance of trace elements reported in the old zircon population⁵ indicates a variability of U/Yb ratios, which suggest different types of source rocks and processes for the genesis of these magmas. These old evolved clasts have most likely been excavated by an impact event during the Early Amazonian period, ~1.5 Ga ago³,⁵,⁸,¹¹⁻¹³ (green bars in Fig. 1), before being lithified¹² and ejected ~ 5 Ma ago¹²,¹³ (Methods). Hence, this regolith breccia provides evidence for the formation of evolved crustal material 4.48 Ga ago on Mars¹⁻¹⁰,¹⁴,¹⁵, and contains clues to the early environment and evolution of Mars.

However, the source region of this unique meteorite and its geological context have so far remained unknown, and with it, a region where the earliest geological records of the planet²,³ are exposed on the surface. Knowing this source region would provide insights into early Mars geological history and crustal extraction²,³. This source region may therefore become a high-priority target for detailed orbital analyses and in-situ exploration¹⁶.

Following a hypervelocity impact, ejecta materials moving faster than the escape velocity (5 km/s)¹⁷ may get through the martian atmosphere and continue their course into interplanetary space to become martian meteorites. Slower debris fall back on the surface in a radial pattern or ray around the primary crater,

forming secondary craters. The presence of 100 meter-size secondaries attests to the freshness of their associated primary craters¹⁸. Using the size and spatial distribution of more than 90 million impact craters >50 m (both primaries and secondaries) detected using a Crater Detection Algorithm (CDA)¹⁸⁻²⁰ on the whole surface of Mars from the global Context Camera (CTX) mosaic²¹, a previous work¹⁸ identified ray systems of secondary craters <150 m associated with 19 large primary craters. For each of them, a formation model age was measured using small craters superposed on their ejecta blanket, and 18 were found younger than 10 Ma old. The analysis of the size frequency distribution of these 18 young crater candidates revealed that those larger than 7 km (i.e., 17 out of 18) align with the predicted number and size of craters accumulated on the whole surface of Mars over the last 8.2 ± 2 Ma²². Hence, those impact craters were found to constitute the complete crater population >7 km in diameter formed on Mars over the last ~10 Ma, potentially responsible for the ejection of martian meteorites¹⁸. One of these craters, Tooting, has already been recognized as the most likely ejection site of the depleted olivine-phyric shergottites launched 1.1 Ma ago, located on the Tharsis volcanic province¹⁸.

Past studies that have proposed parent terrains for the unique NWA 7034 regolith breccia all agree that it must come from the heavily cratered southern Noachian highlands of Mars⁷,¹¹,¹³,²³⁻²⁶ based on their notable geophysical and geochemical affinities, such as the elevated concentrations in potassium (K), thorium (Th), and iron (Fe) in the meteorite⁸,⁹,²³⁻²⁶, the ages of the oldest minerals found in the breccia¹,⁵,⁸⁻¹⁰, its unique magnetic signature⁷, and its visible-infrared reflectance spectra²⁴.

In this work, we search for the most likely site of ejection of the regolith breccia by using four criteria based on its geochemical and geophysical properties as well as its geochronological records (Methods), which we compare to potential sites based on their known properties and geological context: (1) high magnetic field intensity and remanent magnetization at the surface from up-to-date orbital dataset²⁷ (Methods, Fig. 2b, c); (2) high elemental Th and K concentrations²⁸,²⁹ (Methods, Fig. 2d, e) of the areas surrounding each crater candidate; (3) superposition on a Noachian geological unit³⁰ and (4) connection with material from an Early Amazonian impact. We show that only one crater candidate characteristics match with the meteorite properties. The oldest clasts of NWA 7034 and paired stones were excavated ~1.5 Ga

<figure><img src="image_4.png" /><figcaption>**Fig. 1** Summary of NWA 7034 and paired stone radiometric ages, and chronology of major events experienced by the breccia. Dates from each study are reported in Supplementary Data 1. The ejection event is constrained from ²²Ne/²¹Ne cosmic ray exposure ages¹², ¹³. Whole rock (w.r), zircon and or baddeleyite (zr.b), apatite (ap), and augite (aug) on which ages have measured are also mentioned. The different chronometers used in these studies are reported (Sm-Nd, Pb-Pb, U-Pb, K-Ar, U-Th/He, U-Th-Sm/He). Note that green, red and blue boxes correspond to resetting ages of the noted chronometer. Reset ages in green are widely interpreted as the disruption induced by an impact-derived heating event that has excavated the oldest components of the breccia ~1.5 ago⁹, ¹¹, ²³, ²⁶, although its precise age is still unconstrained due to the wide range of isotopic dates reported in the literature.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Distribution of the most likely crater sources for martian meteorites. **a** Global context of the 19 crater candidates¹⁸ (Supplementary Table 1) and location of provinces and rovers (yellow triangles) referred to in the present study. Background: Mars Orbiter Laser Altimeter (MOLA) shaded relief (https://astrogeology.usgs.gov/search/map/Mars/GlobalSurveyor/MOLA/Mars_MGS_MOLA_DEM_mosaic_global_463m). **b, c**: Magnetic field intensity and remanent magnetization at the surface from ref. ²⁷. **d, e**: Potassium and Thorium concentration at the surface from ref. ²⁸,³⁰. Beige area corresponds to discarded provinces due to the weathered basaltic surface contribution (Methods).</figcaption></figure>

ago by an impact that has formed a 40 km crater. The ejecta material of this crater were subsequently ejected by a second impact a few million years ago, which led to the formation of a 10 km crater. The geologic context of the ejection site is consistent with one of the oldest province of Mars, a relic of the differentiated primordial martian crust. This region constitutes a unique record of the first tens of millions of years of the history of Mars.

## Results

Only nine craters of the previously identified large and recent primary craters¹⁸ are located on Noachian highlands terrains³⁰ (Fig. 2a and Supplementary Table 1). A high temperature event is required to account for the metamorphic changes experienced by NWA7034 and paired stones 1.5 Ga ago³,⁸⁻¹⁰,¹²⁻¹⁴ (Methods). None of these nine craters are associated with Amazonian volcanic terrains³⁰, thus making implausible a volcanic origin for the 1.5 Ga event that reset radiochronometers in the apatites, feldspars and zircons within the breccia³¹. On the other hand, two of these nine craters, Karratha and Gasa, respectively labeled 11 and 17 in Fig. 2a (Supplementary Table 1) are superposed on the crater floor or an ejecta blanket of an Amazonian impact crater (younger than ~3.2 Ga old). Such an Amazonian crater can account for the 1.5 Ga event that led to the excavation of the Noachian basement, its brecciation and resetting of the radiochronometers, and the lithification in the ejecta deposits.

While Gasa crater has been noticed in previous studies due to its extended rays visible on thermal imagery³²,³³, Karratha crater

is devoid of such thermally visible rays, and, to our knowledge, has never been reported before as a young primary crater. Gasa crater is located within Cilaos crater (20 km). A model age of 572 ± 110 Ma was estimated from crater counts for the Cilaos impact event (Methods and Supplementary Fig. 4). Cilaos crater is hence too young to be associated with the 1.5 Ga resetting event. Moreover, magnetic signatures and elemental abundances of K and Th reported in the region surrounding Gasa crater are lower compared to those associated with Karratha (Supplementary Fig. 3 and Supplementary Table 1).

Secondary craters from Karratha extend over more than 350 km (Supplementary Fig. 5). It is located within a highly degraded 25 km impact structure, Dampier crater, most likely Noachian in age (Fig. 3). This old crater is filled by the ejecta blanket of a nearby impact crater (Khujirt, D = 40 km), located about 30 km away (rim to rim) in the south-west of Karratha. Karratha crater is superposed on this ejecta material, whose formation occurred during the Early Amazonian period, between 1.25±0.38 and 1.87±0.73 Ga (Methods and Supplementary Fig. 4). An estimate for the thickness of the Khujirt ejecta blanket from scaling laws in ref. ³⁴ gives 60 m where Karratha is located. The maximum depth of debris reaching the martian escape velocity following a 350 m asteroid impact forming a 10 km size crater is ~50 m (ref. ¹⁷). Hence, the large majority of the ejected debris capable of escaping martian gravity following the formation of Karratha are from the Khujirt ejecta blanket, not from the underlying material related to the Dampier crater (Fig. 3).

Moreover, orbital datasets indicate that the Karratha crater is associated with particularly high concentrations of K and Th and
<figure><img src="image_6.png" /><figcaption>**Fig. 3** The NWA 7034 launch site geological context. **a** Perspective view (Context Camera mosaic<sup>21</sup>) and cross-section through Karratha along a SW-NE axis (Supplementary Fig. 5b) as interpreted from numerical modelling simulations. Shades of colors denote impactites (impact melt and in situ breccia; ejecta and fall-back breccia). **b** Schematic of chronological events experienced by the host terrain of the regolith breccia.</figcaption></figure>

elevated magnetic intensities, when compared to all other craters, although these values are lower than those measured in NWA 7034 and paired meteorites (Supplementary Fig. 3). The Karratha crater appears therefore a good candidate for the launch site of NWA 7034, if the formation condition of Khujirt was intense enough to account for the ~1.5 Ga resetting ages measured in augite, apatite, and zircon (Methods). We modeled the temperature of ejecta fragments associated with the formation of the Khujirt crater using the iSALE shock physics code<sup>35-37</sup> (Methods). The simulation indicates that the ejecta fragment experienced a temperature up to 1000 °C where Karratha is located (Supplementary Figs. 6 and 7), compatible with the temperature required to reset the U-Pb in apatites, the Ar-Ar in augites and to disturb the metamict zircons in the breccia<sup>3,5,8,10-15</sup> (500–800 °C, see Methods). With the exception of a shock twin observed in a zircon contained in NWA 7034<sup>38</sup>, the paucity of shock deformation above 10 GPa observed in the breccia<sup>39</sup> is also consistent with an origin as poorly consolidated ejecta material for the regolith breccia.

The four criteria used to locate the crater source of the regolith breccia thus allow nailing down the young crater candidate population to a unique solution. The geological context of Karratha matches the chronology, the lithology, and the magnetic and elemental signatures of the NWA 7034 meteorite group. We conclude that the ejecta blanket on which Karratha is superposed, is associated with the Early Amazonian impact event (Khujirt, in green in Fig. 1) that has excavated the oldest zircons and clasts present in the breccia from the Noachian southern highlands (in grey in Fig. 1). The impact that has formed Karratha crater has subsequently ejected material of the Khujirt crater a few million years ago (in yellow in Fig. 1), making the Karratha impact crater the source of NWA 7034 and paired stones (Fig. 3).

**Implications for the early crust extraction.** The bedrock of the Khujirt crater, located in the north-east of the Terra Cimeria—

Sirenum region (here and after noted TCTS) is believed to be composed of basaltic and more evolved lithologies such as those represented by the monzonitic and noritic clasts, containing the concordant 4.4 Ga zircons<sup>1-5,8-10</sup>. The TCTS province is located between Hesperia Planum and the Tharsis bulge (Fig. 2a), and characterized by the highest concentrations in K (>0.35 wt.%) and Th (>0.35 ppm) measured on Mars from the orbit<sup>28,29</sup> (Fig. 2b, c). TCTS is the only highlands region where the high concentrations of both elements are correlated<sup>29</sup>. Furthermore, it presents the highest magnetic field anomaly (> 5000 nT) and the strongest remanent magnetization (>5 A.m<sup>2</sup>.kg<sup>-1</sup>) on Mars<sup>27</sup> (Fig. 2d, e).

By removing the contribution of the largest impact basins and volcanoes from the gravity field and topography, it has been found that the TCTS province is characterized by the highest crustal thickness of the planet, i.e., >50 km<sup>40</sup>. Since this terrain is overprinted by the Hellas and Argyre basins' ejecta materials, formed earlier than ~4.1 Ga, this region is likely a relic of the most ancient crust<sup>40</sup>, which is confirmed by the location of the NWA 7034 ejection site. The compositional effect of the crust on the magnetic intensity<sup>41</sup> suggests that TCTS remained largely unaffected by demagnetizing processes since the pre-Noachian, and that surface material has not been mixed with the surrounding Noachian regolith. This is consistent with a thick crustal block<sup>40</sup>, whose distinct formation and evolution possibly reflect the first stage of differentiation occurring very early in the history of the planet. The geological stability of this region makes it unlikely that Amazonian and Hesperian hydrothermal processes have contributed to forming magnetite to be the dominant magnetization carrier in this highly magnetized region of the martian crust<sup>42,43</sup>. Instead, analysis of magnetic signatures, K and Th enrichment, and depth of large impact craters in Eridania basins within the TCTS region suggest non-magmatic long-lived hydrothermal systems, heat-driven by radiogenic elements with half-lives of billions of years<sup>44</sup> (e.g., <sup>232</sup>U and <sup>40</sup>K), that might have significantly contributed to the observed
crustal magnetic field throughout the pre-Noachian and Noa-
chian eras42–44. Early hydrothermal circulation in this province
would potentially have sustained life-compatible environment for
a long period of time3,44.

The TCTS region covers about 10% of the planet and has been
interpreted as a crustal block characterized by a geochemically
evolved component. Evolved rocks have been observed and
analyzed on the ground in Gale and Gusev craters by Curiosity
and Spirit, respectively, in the immediate vicinity of this province,
(Fig. 2a). Those igneous rocks show felsic alkaline and sub-
alkaline compositions45–49 that might be explained by fractional
crystallization47,48, and may indicate the presence of a differ-
entiated crust early in the history of the planet46.

The U-Pb and Pb-Pb isotopic compositions in monzonic clasts in the meteorite suggest the existence of an isotopically enriched (relative to the martian mantle) and differentiated crust on Mars¹⁻⁵,¹⁴,⁵⁰ that was extracted before 4.547 Ga¹. According to the Hf isotopic signature of 4.43 Ga zircons recovered from NWA 7034, it has been proposed that a magma ocean crystallized within the first 20 Ma after the accretion of the planet¹. The initial εHf value of these old zircons and models of early magma ocean crystallisation⁵¹⁻⁵³ imply an andesitic composition for the early martian crust¹,⁵,⁵³, although the relationship between U-Pb ages and εHf of the oldest zircons suggest that they cristallized from low ¹⁷⁶Lu/¹⁷⁷Hf magmas potentially of basaltic affinity⁵. Such a crust was then reworked 100 Ma later by impacts⁹,¹⁰, producing the melts from which the old zircons crystallized. The analyses of >3.8 Ga evolved rocks in Gale crater⁴⁶, the inversion of the martian gravity field (constrained by petrological data that support the existence of light evolved crustal components -less dense than basalt- in the southern highlands)⁵⁴ and, finally, seismic data from the *Insight* mission (indicating that the rate of P wave against S wave is compatible with basaltic to andesitic crustal materials⁵⁵ and suggesting that the crust density is <3100 g.cm⁻³, so lower if only composed by basaltic rocks⁵⁶), all point out to the presence of highly ancient evolved crustal components in TCTS.

If andesitic in composition, the crust could either be secondary (reworking of the primordial crust) or primordial. If primordial, basaltic to andesitic melts that originated from a deep mantle source might have crystallized⁵⁷. However, isobaric partial melting experiments⁵⁸ and adiabatic ascent of primitive mantle compositions⁴⁸ argue against the formation of andesitic magmas under such a scenario. If secondary, the primordial crust would have been extracted from the magma ocean extremely early, before 4.547 Ga ago, i.e., <20 Ma after solar system formation¹. Alternatively, if the primordial crust was basaltic in composition, its differentiation and re-melting might have resulted in an evolved crust as observed for the continental crust on Earth. Another possibility is the absence of a global magma ocean, as suggested by ref. ⁸, where a low-degree of partial melting of a fertile mantle could produce an enriched crust with rare-earth element patterns similar to those observed within the regolith breccia. In the later case, scattered magma oceans could have occurred in locations different from the source terrain of the breccia. Differentiated primordial crustal blocks of the planet such as that in TSTC would have been formed by simple partial melting shortly after the accretion of the planet.

In any case, we suggest that clasts contained in the regolith breccia are representative of the TCTS province, making this region a relic of the early crustal processes on Mars, and thus, a region of high interest for future missions. The study of TCTS would help us unravel the conditions of formation and the first evolution stage of Mars, and by extension of all terrestrial planets, considering the fact that, in light of these findings, early crustal processes appear uniquely preserved and accessible on Mars. The

flanks and central peaks of large and preserved craters within this
region might constitute outcrops of high interest, containing the
missing geological clues to the early-stage evolution of Mars.

Methods

**NWA 7034 and pair characteristics.** The diversity of clasts contained in the breccia makes this meteorite one of the martian samples with the most complex history, recording multiple events, from the crystallization of the martian primary crust to the ejection of the rock. The crystal clasts contained in the breccia include low-Ca pyroxene, augite, plagioclase and alkali feldspar, with a near absence of olivine. The fine-grained matrix contains pyroxene, plagioclase, iron oxides, Cl-apatite, chromite and pyrite⁵,⁸,¹¹,²⁶,⁵⁹. A vitrophyre melt clast containing the highest Ni abundance ever measured in any martian meteorite or igneous rock (1020 ppm) has also been reported⁵⁰, suggesting contamination by a chondritic impactor, also seen in all lithic clast types except orthopyroxenite⁹,²⁶. Basaltic clasts and bulk matrix compositions of this meteorite have been found to be analogous to igneous rocks analyzed by the Spirit rover within Gusev crater as well as the average martian crust composition determined by the Gamma-Ray Spectrometer (GRS) aboard the Mars Odyssey spacecraft⁵³. However, some evolved clasts exhibiting trachyandesitic, basaltic, and Fe-, Ti-, and P- rich (FTP) lithologies constitute rock types similar to those analyzed by the Curiosity rover in Gale crater⁴⁵,⁴⁶,⁴⁸,⁴⁹.

This meteorite exhibits unique characteristics, including one of the highest concentrations in potassium and thorium ever measured in a martian meteorite (Supplementary Fig. 1a). Coupled with isotopic analysis (¹⁴⁷Sm/¹⁴⁴Nd and ¹⁷⁶Lu/¹⁷⁷Hf), this suggests that the mantle source of NWA 7034 is isotopically and chemically distinct from the other martian meteorites mantle sources²⁻¹⁰. The unique magnetic mineralogy of the breccia makes NWA 7034 the most magnetized martian meteorite⁷ with remanent magnetization (20–60 A/m) one order of magnitude higher than any other martian meteorites (Supplementary Fig. 1b). Metamict zircons⁸,¹⁰apatites¹⁵, augite¹¹ and alkali-feldspars in leucocratic clasts¹² age measurements suggest the occurrence of one single metamorphic event ~1.5 Ga ago (in green in Fig. 1). This event has been interpreted either as a volcanic one³¹ inducing protracted metamorphism or as an impact that would have excavated the oldest component of the meteorite, with temperatures ranging between 500 °C and 800 °C (refs. 3,9,11,26). Grain shape and size distributions of the breccia clasts indicate they were likely deposited by impact-ejecta materials under base surge conditions¹¹. Moreover, the petrographic similarity of the meteorite with terrestrial suevite⁹,¹¹, accretionary dust rims²³), and the presence of stishovite⁶⁰ (a high-pressure polymorph of SiO₂) all favor an impact origin for the resetting of the U-Pb in augite and alkali-feldspars and the disturbance of the metamict zircons in the breccia⁸,¹⁰,¹¹. Suevite is described as a polymict breccia containing lithic and mineral fragments that exhibit a variety of shock metamorphism stages⁶¹,⁶². It usually composes the upper layer of the crater floor cavity and the proximal ejecta layer⁶¹. Based on the grain size and shape, previous work⁹ suggested that the meteorite is representative of a proximal ejecta blanket deposited in a pyroclastic flow regime, consistent with conclusions from ref. ¹¹, according to which the accretionary dust rims seen in the breccia have been formed under base surge conditions following an impact event, ~1.5 Ga ago.

Eight young detrital zircons were recently discovered⁵ in NWA 7533 with ages from 1548.0 ± 8.8 Ma to 299.5 ± 0.6 Ma (in red in Fig. 1). Analysis of their isotopic composition indicates a common mantle source, sampled by deep-seated magmatic activity, possibly representative of Tharsis or Elysium volcanic provinces⁵. These zircons were likely transported by eolian processes to the source region of the breccia⁵. The oldest grain being 1548 Ma old⁵, this zircon population is consistent with an excavation ~1.5 Ga ago, possibly triggered by an impact. The breccia lithification would have occurred subsequently, 225 Ma ago¹² (in blue in Fig. 1), while the ejection took place ~5 Ma ago¹²,¹³ (in yellow in Fig. 1).

**Constraints from orbital dataset.** We compare the abundance of K and Th as well as the magnetic field intensity and the magnetization of the surface of Mars derived from orbital measurements at the immediate vicinity of each crater candidate with those of the breccia. The concentrations in K and Th are from the Gamma-Ray Spectrometer (GRS) aboard the Mars Global Surveyor (MGS)²⁸,²⁹. Although the spatial resolution of the data is low (5°×5°, corresponding to ~296 km at the equator), this dataset offers a consistent method to compare the concentration of both elements qualitatively between several regions where the crater candidates are located⁶³. We report here both the pixel value of the concentration of K and Th at the crater location (Supplementary Table 1) and the bilinear interpolation computed with a radius of 296 km around each crater centroid (Supplementary Table 2).

A recent model of the crustal magnetic field of Mars has been computed using the Mars Atmosphere and Volatile EvolutionN (MAVEN) magnetometer²⁷. It is the highest resolution model of the martian magnetic field at the surface (spatial resolution: ~100 km/px, magnetic field resolution: <1 nT), allowing the identification of small-scale features associated with geological signatures (Fig. 2b). From this model, the authors also derived the equivalent magnetization distribution (Fig. 2c). Because magnetic fields and magnetization spatial variation exist at a small scale (smaller than the spatial resolution of the model and thus
lower than the magnetization currently detectable from orbit²⁵), we report for each crater candidate two values of the magnetic field intensity and equivalent magnetization: (1) the pixel value associated with the centroid of the crater (Supplementary Table 1 and Supplementary Fig. 3) and (2) a bilinear interpolation computed from a 100 km buffer (resolution of the model) computed around each crater centroid (Supplementary Table 2).

In order to distinguish crater candidates associated with a host terrain exhibiting relatively high values of the four orbital datasets considered here, we compute the first standard deviation from the kernel density distribution of the data. Values higher than +1σ are considered relatively high (Supplementary Fig. 2) and used to discriminate areas with elemental concentration and magnetic properties that might be consistent with the exceptional characteristics of the breccia qualitatively (Supplementary Fig. 3). Enrichment in K and Th in the meteorite being related to magmatic processes, we chose to discard the northern lowlands, and more generally areas above the dichotomy, to compute the distribution of the K and Th concentration due to the weathered basaltic surface contribution⁶⁴ that may not represent the bedrock chemistry on those provinces.

Supplementary Fig. 3 presents the chemical and magnetic signatures for each crater candidate (in pixel values). The high values range associated with each dataset are represented by the white area and impact craters are color-coded as followed: craters located on Noachian geological unit³⁰ and superposed on the material of an Amazonian impact crater (cavity or ejecta) are in green. Those that are only located on Noachian material are in orange, and craters superposed on the ejecta of an Amazonian impact crater are in red. Finally, if none of these two criteria are respected, craters are shown in grey.

**Model age derivation of impact events.** Karratha crater is superposed on the ejecta blanket of a 40 km crater (Khujirt, in blue in Fig. 3) that filled the cavity of an older crater, most likely of Noachian age (Dampier, in orange in Fig. 3). Using the Context Camera (CTX) mosaic²¹, we mapped impact craters superposed on both the 40 km crater cavity (for D > 100 m) and its ejecta blanket (for D > 500 m) to estimate the age of the material surrounding Karratha, i.e., the ejecta blanket of the 40 km crater (Supplementary Fig. 4c, d). For this, crater mapping is performed by using the CraterTools software⁶⁵ and Crater-Size Frequency Distributions (CSFD) are loaded into CraterStats II⁶⁶ and fitted with an isochron using a standard chronology model²². We derived all model ages using the differential representation⁶⁷,⁶⁸. Compared to a cumulative plot, the differential representation allows easier recognition of any resurfacing event contribution, the presence of potential overprinting craters formed prior to the ejecta blanket in the population of mapped craters, or secondary craters⁶⁷. Each point in the CSFD is independent of the subsequent larger diameter bins. Like other representations, the binning of the data can lead to biases when the CSFD is fitted with an isochron⁶⁹ when a small number of impact craters is used to derive model ages⁷⁰. We solved this statistical disadvantage by using the Poisson probability-density function-fitting technique⁶⁹. This solution allows an exact prediction of the model crater chronology model according to the CSFD considered, whatever the chosen binning technique. The crater count on the ejecta blanket of Khujirt crater leads to a model age of $1.87 \pm 0.73_{-0.65}^{+0.65}$ Ga, consistent with the model age obtained from crater count on the crater cavity $1.25 \pm 0.38_{-0.32}^{+0.38}$ Ga.

Gasa crater is located within the Cilaos impact crater. The ejecta blanket of the latter has been mapped (Supplementary Fig. 4c) and crater counts have been performed using the CTX global mosaic²¹ down to ~250 m to estimate the age of the Cilaos impact crater. A middle Amazonian model age²² ($572 \pm 110$ Ma) is obtained using the same technique and chronology model used for Karratha crater (Supplementary Fig. 4d). Even taking into account larger uncertainties of the crater count method linked to the influence of the terrain rheology on the crater size⁷¹ or potential fluctuation in the impact cratering rate and crater production⁶⁷,⁷²⁻⁷⁴, the model age derived for the Cilaos impact event is inconsistent with the $1.5$ Ga resetting age of some minerals in the meteorite⁸,¹⁰⁻¹². In summary, the impact that has formed the ejecta blanket of Khujirt on which Karratha is superposed occurred most likely during the early Amazonian period²⁰, ~1.5 Ga ago. This impact event is the only one that could match the $1.5$ Ga resetting age observed in the meteorite (in green on Fig. 1) as well as the age of the oldest detrital zircon⁵ (in red on Fig. 1), consistent with the excavation of the parent rock of NWA 7034.

**Impact crater modeling and ejecta temperature.** We model the shock and post-shock temperature of the ejecta curtain associated with the formation of a ~40 km crater (i.e., Khujirt) using iSALE-2D (refs. ³⁵⁻³⁷). (Supplementary Fig. 6 and Supplementary Table 3). For this, we use a 4.5 km diameter impacting the surface at $9.6$ km/s (ref. ⁷⁴). Impactor and target were modeled using the equation of state for dunite⁷⁵ and basalts⁷⁶, respectively. The temperature gradient in the upper crust is assumed to be $15$ K/km (ref. ⁷⁷), which is on the hotter side of present-day Mars and could be appropriate for Mars at $1.5$ Ga. We find that ejecta material experienced a heating ranging between $0$ and $>1000$ °C where Karratha was formed (~1.5 crater radii from the rim of Khujirt crater, Supplementary Fig. 7). This is compatible with the temperature required to disturb the Ar-Ar and U-Pb radio-chronometers in the breccia³,⁵,⁸,¹⁰⁻¹³,¹⁵ (500–800 °C). Depth of origin for this ejecta material does not exceed 5 km. We note that the heating of the landing ejecta is directly dependent on the impactor speed. Increasing impact speed and lowering the size of the projectile so that the resulting crater diameter is the same, the mean

temperature range in ejecta will elevates but remains in the observed range. Furthermore, any macro voids within the falling ejecta could significantly elevate the temperature in the falling ejecta⁷⁸. It is therefore difficult to estimate the duration of the shock temperatures necessary for chronometer resetting, as that would depend on a number of parameters such as mineral composition, sample/ejecta fragment size, porosity, heterogeneity, etc. Direct comparison between numerical modelling outcomes and laboratory measurements is not directly applicable, and merits further work.

## Data availability

The data that support the findings of this study are available within the paper and the Supplementary Data file.

## Code availability

The numerical impact crater formation were made using the iSALE shock physics hydrocode. At present, iSALE is not fully open source. Application for use of iSALE can be made via https://isale-code.github.io/. Any recent stable release can be used to reproduce the data presented. We used the IDL 5.2 software (L3Harris geospatial https://www.l3harrisgeospatial.com/Software-Technology/IDL) to run the craterStats II software available at https://www.geo.fu-berlin.de/en/geol/fachrichtungen/planet/softwarealgorithm, and the ESRI's ArcGIS 10.8.1 software suite (ESRI https://www.esri.com/en-us/arcgis/about-arcgis/overview) and Matlab (https://au.mathworks.com/products/matlab.html) to produce the maps.

Received: 29 December 2021; Accepted: 8 June 2022;

Published online: 12 July 2022

## References

1. Bouvier, L. C. et al. Evidence for extremely rapid magma ocean crystallization and crust formation on Mars. *Nature* **558**, 586–589 (2018).

2. Udry, A. et al. What martian Meteorites Reveal About the Interior and Surface of Mars. *JGR: Planets* **125**, e2020JE006523 (2020).

3. Goodwin, A., Garwood, R. J., Tartèse, R. A review of the “Black Beauty” martian regolith breccia and its martian habitability record. Astrobiology, ahead of print (2022). https://doi.org/10.1089/ast.2021.0069

4. Armytage, R. M., Debaille, V., Brandon, A. D. & Agee, C. B. A complex history of silicate differentiation of Mars from Nd and Hf isotopes in crustal breccia NWA 7034. *EPSL* **502**, 274–283 (2018).

5. Costa, M. M. et al. The internal structure and geodynamics of Mars inferred from a 4.2-Gyr zircon record. *PNAS* **117**, 30973–30979 (2020).

6. Nyquist, L. E. et al. Rb-Sr and Sm-Nd isotopic and REE studies of igneous components in the bulk matrix domain of martian breccia Northwest Africa 7034. *MaPS* **51**, 483–498 (2016).

7. Gattacceca, J. et al. martian meteorites and martian magnetic anomalies: A new perspective from NWA 7034: martian meteorites & magnetic anomalies. *GRL* **41**, 4859–4864 (2014).

8. Humayun, M. et al. Origin and age of the earliest martian crust from meteorite NWA 7533. *Nature* **503**, 513–516 (2013).

9. McCubbin, F. M. et al. Geologic history of martian regolith breccia Northwest Africa 7034: Evidence for hydrothermal activity and lithologic diversity in the martian crust: Geologic History of NWA 7034. *JGR: Planets* **121**, 2120–2149 (2016).

10. Hu, S. et al. Ancient geologic events on Mars revealed by zircons and apatites from the Martian regolith breccia NWA 7034. *MaPS* **54**, 850–879 (2019).

11. MacArthur, J. L. et al. Mineralogical constraints on the thermal history of martian regolith breccia Northwest Africa 8114. *Geochim. Cosmochim. Acta* **246**, 267–298 (2019).

12. Cassata, W. S. et al. Chronology of martian breccia NWA 7034 and the formation of the martian crustal dichotomy. *Sci. Adv.* **4**, eaap8306 (2018).

13. Cartwright, J. A., Ott, U., Herrmann, S. & Agee, C. B. Modern atmospheric signatures in 4.4 Ga martian meteorite NWA 7034. *EPSL* **400**, 77–87 (2014).

14. Deng, Z. et al. Early oxidation of the martian crust triggered by impacts. *Sci. Adv.* **6**, eacb4941 (2020).

15. Bellucci, J. J. et al. Pb-isotopic evidence for an early, enriched crust on Mars. *EPSL* **410**, 34–41 (2015).

16. Rapin, W. et al. Critical knowledge gaps in the Martian geological record: A rationale for regional-scale in situ exploration by rotorcraft mid-air deployment. Bulletin of the AAS. 53 (2021). https://doi.org/10.3847/25c2cfeb.4986bd82

17. Artemieva, N. & Ivanov, B. Launch of martian meteorites in oblique impacts. *Icarus* **171**, 84–101 (2004).
18. Lagain, A. et al. The Tharsis mantle source of depleted sherгottites revealed by 90 million impact craters. *Nat. Comm.* **12**, 6352 (2021).

19. Benedix, G. K. et al. Deriving surface ages on mars using automated crater counting. *ESS* 7, e2019EA001005 (2020).

20. Lagain, A. et al. Model age derivation of large martian impact craters, using automatic crater counting methods. *ESS* 8, e2020EA001598 (2021).

21. Dickson, J. L., Kerber, L. A., Fassett, C. I. & Ehlmann, B. L. A Global, blended CTX Mosaic of Mars with vectorized seam mapping: a new mosaicking pipeline using principles of non-destructive image editing. 49th LPSC, Abstract #2480 (2018). https://www.hou.usra.edu/meetings/lpsc2018/pdf/2480.pdf

22. Hartmann, W. K. martian cratering 8: Isochron refinement and the chronology of Mars. *Icarus* **174**, 294–320 (2005).

23. Wittmann, A. et al. Petrography and composition of martian regolith breccia meteorite Northwest Africa 7475. *MaPS* **50**, 326–352 (2015).

24. Beck, P. et al. A Noachian source region for the “Black Beauty” meteorite, and a source lithology for Mars surface hydrated dust? *EPSL* **427**, 104–111 (2015).

25. Cannon, K. M., Mustard, J. F. & Agee, C. B. Evidence for a widespread basaltic breccia component in the martian low-albedo regions from the reflectance spectrum of Northwest Africa 7034. *Icarus* **252**, 150–153 (2015).

26. Hewins, R. H. et al. Regolith breccia Northwest Africa 7533: Mineralogy and petrology with implications for early Mars. *MaPS* **52**, 89–124 (2017).

27. Langlais, B., Thébault, E., Houliez, A., Purucker, M. E. & Lillis, R. J. A new model of the crustal magnetic field of Mars using MGS and MAVEN. *JGR: Planets* **124**, 1542–1569 (2019).

28. Boynton, W. V. et al. Concentration of H, Si, Cl, K, Fe, and Th in the low- and mid- latitude regions of Mars. *JGR* **112**, 12–99 (2007).

29. Taylor, G. J. et al. Bulk composition and early differentiation of Mars. *JGR* **112**, 03–10 (2007).

30. Tanaka, K. L. et al. Geologic map of Mars. U.S.G.S. Scientific Investigations Map 3292, scale 1:20,000,000, pamphlet 43, (2014). https://doi.org/10.3133/sim3292

31. Lindsay, F. N. et al. 40Ar/39Ar ages of Northwest Africa 7034 and Northwest Africa 7533. *MaPS* **56**, 515–545 (2021).

32. Schon, S. C., Head, J. W. & Fassett, C. I. Unique chronostratigraphic marker in depositional fan stratigraphy on Mars: Evidence for ca. 1.25 Ma gully activity and surficial meltwater origin. *Geology* **37**, 207–210 (2009).

33. Head, J. W., Mustard, J. F., Kreslavsky, M. A., Milliken, R. E. & Marchant, D. R. Recent ice ages on Mars. *Nature* **426**, 797–802 (2003).

34. Sturm, S., Kenkmann, T. & Hergarten, S. Ejecta thickness and structural rim uplift measurements of martian impact craters: Implications for the rim formation of complex impact craters. *JGR: Planets* **121**, 1026–1053 (2016).

35. Amsden, A. A., Ruppel, H. M. & Hirt, C. W. SALE: a simplified ALE computer program for fluid flow at all speeds. United States: N. p. (1980). https://doi.org/10.2172/5176006

36. Collins, G. S., Melosh, H. J. & Ivanov, B. A. Modeling damage and deformation in impact simulations. *MaPS* **39**, 217–231 (2004).

37. Wünnemann, K., Collins, G. S. & Melosh, H. J. A strain-based porosity model for use in hydrocode simulations of impacts and implications for transient crater growth in porous targets. *Icarus* **180**, 514–527 (2006).

38. Cox, M. A. et al. Impact and habitability scenarios for early Mars revisited based on a 4.45-Ga shocked zircon in regolith breccia. Science Advances 8 https://doi.org/10.1126/sciadv.abl7497.

39. Leroux, H. et al. Exsolution and shock microstructures of igneous pyroxene clasts in the Northwest Africa 7533 martian meteorite. *MaPS* **51**, 932–945 (2016).

40. Bouley, S. et al. A thick crustal block revealed by reconstructions of early Mars highlands. *Nat. Geos.* **13**, 105–109 (2020).

41. AlHantoobi, A., Buz, J., O'Rourke, J. G., Langlais, B. & Edwards, C. S. Compositional enhancement of crustal magnetization on Mars. *GRL* 48 (2021). https://doi.org/10.1029/2020GL090379

42. Ruiz, J. The very early thermal state of Terra Cimmeria: Implications for magnetic carriers in the crust of Mars. *Icarus* **203**, 454–459 (2009).

43. Michalski, J. R. et al. Ancient hydrothermal seafloor deposits in Eridania basin on Mars. *Nat. Com.* **8**, 15978 (2017).

44. Ojha, L. et al. Amagmatic hydrothermal systems on Mars from radiogenic heat. *Nat. Com.* **12**, 1754 (2021).

45. Cousin, A. et al. Geochemistry of the Bagnold dune field as observed by ChemCam and comparison with other aeolian deposits at Gale Crater: ChemCam Results From Bagnold Dunes, Mars. *JGR: Planets* **122**, 2144–2162 (2017).

46. Sautter, V. et al. In situ evidence for continental crust on early Mars. *Nat. Geos* **8**, 605–609 (2015).

47. McSween, H. Y. et al. Alkaline volcanic rocks from the Columbia Hills, Gusev crater, Mars. *JGR* **111**, E09591 (2006).

48. Payré, V. et al. Constraining ancient magmatic evolution on Mars using crystal chemistry of detrital igneous minerals in the sedimentary bradbury group, gale crater, Mars. *JGR: Planets* **125**, e2020JE006467 (2020).

49. Payré, V. et al. Alkali trace elements in Gale crater, Mars, with ChemCam: Calibration update and geological implications. *JGR: Planets* **122**, 650–679 (2017).

50. Udry, A., Lunning, N. G., McSween, H. Y. & Bodnar, R. J. Petrogenesis of a vitrophyre in the martian meteorite breccia NWA 7034. *Geochim. Cosmochim. Acta* **141**, 281–293 (2014).

51. Elkins-Tanton, L. T. Linked magma ocean solidification and atmospheric growth for Earth and Mars. *EPSL* **271**, 181–191 (2008).

52. Debaille, V., Brandon, A. D., O'Neill, C., Yin, Q. Z. & Jacobsen, B. Early martian mantle overturn inferred from isotopic composition of nakhlite meteorites. *Nat. Geos.* **2**, 548–552 (2009).

53. Sautter, V. et al. Magmatic complexity on early Mars as seen through a combination of orbital, in-situ and meteorite data. *Lithos* **254-255**, 36–52 (2016).

54. Baratoux, D. et al. Petrological constraints on the density of the Martian crust. *J. Geophys. Res.: Planets* **119**, 1707–1727 (2014).

55. Deng, S. & Levander, A. Autocorrelation reflectivity of Mars. *Geophys. Res. Lett.* **47**, e2020GL089630 (2020).

56. Knapmeyer-Endrun, B. et al. Thickness and structure of the martian crust from InSight seismic data. *Science* **373**, 438–443 (2021).

57. Elkins-Tanton, L. T., Hess, P. C., & Parmentier, E. M. Possible formation of ancient crust on Mars through magma ocean processes. *J. Geophys. Res. Planets.* **110** (2005). https://doi.org/10.1029/2005JE002480

58. Collinet, M., Médard, E., Charlier, B., Vander Auwera, J. & Grove, T. L. Melting of the primitive Martian mantle at 0.5–2.2 GPa and the origin of basalts and alkaline rocks on Mars. *EPSL* **427**, 83–94 (2015).

59. Santos, A. R. et al. Petrology of igneous clasts in Northwest Africa 7034: implications for the petrologic diversity of the martian crust. *Geochim. Cosmochim. Acta* **157**, 56–85 (2015).

60. Daly, L. et al. High Pressure Excursions in the Matrix of martian Meteorite Northwest Africa (NWA) 11522. 81st Annual Meeting of the Met. Soc., Abstract #6237 (2018). https://www.hou.usra.edu/meetings/metsoc2018/pdf/6237.pdf

61. Fettes, D. & Desmons, J. Metamorphic rocks: A classification and glossary of terms. *Recommendations of the International Union of Geological Sciences Subcommission on the Systematics of Metamorphic Rocks* 256 (2007).

62. Kenkmann, T., Poelchau, M. H. & Wulf, G. Structural geology of impact craters. *J. Struct. Geol.* **62**, 156–182 (2014).

63. Baratoux, D. et al. The impact of measurement scale on the univariate statistics of K, Th, and U in the Earth crust. *ESS* **8**, e2021EA001786 (2021).

64. Wyatt, M. B. & McSween, H. Y. Spectral evidence for weathered basalt as an alternative to andesite in the northern lowlands of Mars. *Nature* **417**, 263–266 (2002).

65. Kneissl, T., van Gasselt, S. & Neukum, G. Map-projection-independent crater size-frequency determination in GIS environments - New software tool for ArcGIS. *PSS* **59**, 1243–1254 (2011).

66. Michael, G. G. & Neukum, G. Planetary surface dating from crater size-frequency distribution measurements: Partial resurfacing events and statistical age uncertainty. *EPSL* **294**, 223–229 (2010).

67. Fassett, C. I. Analysis of impact crater populations and the geochronology of planetary sur- faces in the inner solar system: Crater Populations and Surface Chronology. *JGR: Planets* **121**, 1900–1926 (2016).

68. Crater Analysis Techniques Working Group. Standard techniques for presentation and analysis of crater size- frequency data. *Icarus* **37**, 467–474 (1979).

69. Michael, G. G., Kneissl, T. & Neesemann, A. Planetary surface dating from crater size- frequency distribution measurements: Poisson timing analysis. *Icarus* **277**, 279–285 (2016).

70. Williams, J.-P. et al. Dating very young planetary surfaces from crater statistics: A review of issues and challenges. *MaPS* **53**, 554–582 (2018).

71. Prieur, N. C. et al. The effect of target properties on transient crater scaling for simple craters. *JGR: Planets* **122**, 1704–1726 (2017).

72. Lagain, A., Bouley, S., Baratoux, D., Costard, F. & Wieczorek, M. Impact cratering rate consistency test from ages of layered ejecta on Mars. *Planet. Space Sci.* **180**, 104755 (2020).

73. Lagain, A. et al. Has the impact flux of small and large asteroids varied through time on Mars, the Earth and the Moon? *EPSL* **579**, 117362 (2022).

74. Ivanov, B. A. Mars/Moon Cratering Rate Ratio Estimates. *Sp. Sci. Rev.* **96**, 87–104 (2001).

75. Benz, W., Cameron, A. G. W. & Melosh, H. J. The origin of the Moon and the single-impact hypothesis III. *Icarus* **81**, 113–131 (1989).

76. Benz, W. & Asphaug, E. Catastrophic disruptions revisited. *Icarus* **142**, 5–20 (1999).

77. Plesa, A. et al. The thermal state and interior structure of Mars. *GRL* **45**, 198–12,209 (2018).

78. Güldemeister, N. et al. Propagation of impact-induced shock waves in porous sandstone using mesoscale modeling. *MaPS* **48**, 115–133 (2013).
## Acknowledgements

The authors acknowledge Benoit Langlais for providing guidance and magnetic field data, and Denis Fougerouse for discussions and guidance regarding age measurements used in this study. We gratefully acknowledge the developers of iSALE-2D, including Gareth Collins, Kai Wünnemann, Dirk Elbeshausen, Tom Davison, Boris Ivanov and Jay Melosh. This research was funded by the Australian Research Council grants DP170102972, DP210100336, DP180100661, DE180100584, and FT210100063, Curtin University, the Western Australian Government, and the Australian Government.

## Author contributions

A.L. S.B., and B.Z. conceived the project. A.L. performed the orbital dataset analyses, derived model ages, and investigated the geological context of all crater candidates. K.M. and A.R. performed the impact simulations. A.L., S.B., B.Z., D.B., V.P., L.S.D., and R.H. discussed the results. A.L. drafted the manuscript with contributions and critical feedback from S.B., B.Z., K.M., A.R., D.B., V.P., L.S.D., N.E.T., R.H., G.K.B., V.M., K.S., and P.A.B.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-022-31444-8.

**Correspondence** and requests for materials should be addressed to A. Lagain.

**Peer review information** *Nature Communications* thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022