Search

The thermal and magnetic histories of planetesimals provide unique insights into the formation and evolution of Earth's building blocks. These histories can be gleaned from meteorites by using numerical models to translate measured properties into planetesimal behaviour. In this paper, we present a new 1D planetesimal thermal evolution and dynamo generation model. This magnetic field generation model is the first of a differentiated, mantled planetesimal that includes both mantle convection and non-eutectic core solidification. We have improved fundamental aspects of mantle heat transport by including a more detailed viscosity model and stagnant lid convection parametrisations consistent with internal heating. We have also added radiogenic heating from $^{60}Fe$ in the metallic Fe-FeS core. Additionally, we implement a combined thermal and compositional buoyancy flux, as well as the latest magnetic field scaling laws to predict magnetic field strengths during the planetesimal's thermal evolution until core solidification is complete. We illustrate the consequences of our model changes with an example run for a 500 km radius planetesimal. These effects include more rapid erosion of core thermal stratification and longer duration of mantle convection compared to previous studies. The additional buoyancy from core solidification has a marginal effect on dynamo strength, but for some initial core sulfur contents it can prevent cessation of the dynamo when mantle convection ends. Our model can be used to investigate the effects of individual parameters on dynamo generation and constrain properties of specific meteorite parent bodies. Combined, these updates mean this model can predict the most reliable and complete magnetic field history for a planetesimal to date, so is a valuable tool for deciphering planetesimal behaviour from meteorite properties., Comment: 38 pages, 12 figures, 3 tables. Accepted for publication in Icarus. Supplementary material is a pdf in the TeX source file

Despite being pivotal to the habitability of our planet, the process by which Earth gained its present-day hydrogen budget is unclear. Due to their isotopic similarity to terrestrial rocks across a range of elements, the meteorite group that is thought to best represent Earth's building blocks is the enstatite chondrites (ECs). Because of ECs' nominally anhydrous mineralogy, these building blocks have long been presumed to have supplied negligible hydrogen to the proto-Earth. However, recent bulk compositional measurements suggest that ECs may unexpectedly contain enough hydrogen to readily explain Earth's present-day water abundance. Together, these contradictory findings mean the contribution of ECs to Earth's water budget is currently unclear. As such, it is uncertain whether appreciable water is a systematic outcome of Earth's formation. Here, we explore the amount of hydrogen in ECs as well as the phase that may carry this element using sulfur X-ray absorption near edge structure (S-XANES) spectroscopy. We find that most (~80 wt%) of the measured bulk hydrogen concentration of these meteorites is bonded to sulfur. Moreover, the concentration of the S-H bond is intimately linked to the abundance of micrometre-scale pyrrhotite (Fe1-xS, 0

Accreting in the first few million years (Ma) of the Solar System, planetesimals record conditions in the protoplanetary disc and are the remnants of planetary formation processes. The meteorite paleomagnetic record carries key insights into the thermal history of planetesimals and their extent of differentiation. The current paradigm splits the meteorite paleomagnetic record into three magnetic field generation epochs: an early nebula field ($\lesssim$5 Ma after CAI formation), followed by thermal dynamos ($\sim$5-34 Ma after CAI formation), then a gap in dynamo generation, before the onset of core solidification and compositional dynamos. These epochs have been defined using current thermal evolution and dynamo generation models of planetesimals. Here, we demonstrate these epochs are not as distinct as previously thought based on refined thermal evolution models that include more realistic parametrisations for mantle convection, non-eutectic core solidification, and radiogenic $^{60}Fe$ in the core. We find thermal dynamos can start earlier and last longer. Inclusion of appreciable $^{60}Fe$ in the core brings forward the onset of dynamo generation to $\sim$1-2 Ma after CAI formation, which overlaps with the existence of the nebula field. The second epoch of dynamo generation begins prior to the onset of core solidification, suggesting this epoch is not purely compositionally driven. Planetesimal radius is the dominant control on the strength and duration of dynamo generation, and the choice of reference viscosity can widen the gap between epochs of dynamo generation from 0-200 Ma. Overall, variations in planetesimal properties lead to more variable timings of different planetesimal magnetic field generation mechanisms than previously thought. This alters the information we can glean from the meteorite paleomagnetic record about the early Solar System., Comment: 29 pages, 8 figures. Accepted for publication in Earth and Planetary Science Letters

Astronomical observations and isotopic measurements of meteorites suggest that substructures are common in protoplanetary disks and may even have existed in the solar nebula. Here, we conduct paleomagnetic measurements of chondrules in CO carbonaceous chondrites to investigate the existence and nature of these disk sub-structures. We show that the paleomagnetism of chondrules in CO carbonaceous chondrites indicates the presence of a 101 $\pm$ 48 $\mu$T field in the solar nebula in the outer solar system ($\sim$3 to 7 AU from the Sun). The high intensity of this field relative to that inferred from inner solar system ($\lesssim$3 AU) meteorites indicates a factor of $\sim$5 to 150 mismatch in nebular accretion between the two reservoirs. This suggests substantial mass loss from the disk associated with a major disk substructure, possibly due to a magnetized disk wind., Comment: 16 pages, 3 figures, published in Science Advances

Metallic microstructures in slowly-cooled iron-rich meteorites reflect the thermal and magnetic histories of their parent planetesimals. Of particular interest is the cloudy zone, a nanoscale intergrowth of Ni-rich islands within a Ni-poor matrix that forms below 350{\deg}C by spinodal decomposition. The sizes of the islands have long been recognized as reflecting the low-temperature cooling rates of meteorite parent bodies. However, a model capable of providing quantitative cooling rate estimates from island sizes has been lacking. Moreover, these islands are also capable of preserving a record of the ambient magnetic field as they grew, but some of the key physical parameters required for recovering reliable paleointensity estimates from magnetic measurements of these islands have been poorly constrained. To address both of these issues, we present a numerical model of the structural and compositional evolution of the cloudy zone as a function of cooling rate and local composition. Our model produces island sizes that are consistent with present-day measured sizes. This model enables a substantial improvement in the calibration of paleointensity estimates and associated uncertainties. In particular, we can now accurately quantify the statistical uncertainty associated with the finite number of islands and the uncertainty on their size at the time of the record. We use this new understanding to revisit paleointensities from previous pioneering paleomagnetic studies of cloudy zones. We show that these could have been overestimated but nevertheless still require substantial magnetic fields to have been present on their parent bodies. Our model also allows us to estimate absolute cooling rates for meteorites that cooled slower than 10000{\deg}C My-1. We demonstrate how these cooling rate estimates can uniquely constrain the low-temperature thermal history of meteorite parent bodies., Comment: Manuscript resubmitted after revisions