FormingWorlds / PROTEUS

Coupled atmosphere-interior framework to simulate the temporal evolution of rocky planets.
https://fwl-proteus.readthedocs.io
Apache License 2.0
12 stars 1 forks source link

Volatile delivery module #168

Open nichollsh opened 2 months ago

nichollsh commented 2 months ago

We now have the initial components of an escape model integrated into PROTEUS. However, it is possible that accretion is also contemporaneous with magma ocean evolution. It would be interesting to explore what role this plays in the changing volatile inventory/distribution of young planets.

It should be relatively easy to "add" additional volatile mass to the inventory as the planet evolves, in the same manner in which escape removes it. We would need some estimate of the rate at which this occurs.

@timlichtenberg is this a good idea and is there a nice way to estimate this?

timlichtenberg commented 2 months ago

Yes, having time-variable inventories of H-C-N-S together with fO2 is a good and important first step for accretion. Ultimately, we would need a varying (growing) size of the aragog grid to additional add ongoing mass accretion. The latter is something that would be highly valuable, but more a long-term idea and not first priority. The former (variable H-C-N-S + fO2) should be pretty straightforward to do I believe.

timlichtenberg commented 2 months ago

Regarding the estimation of mass accretion: this is pretty much unsolved in my view from a planet formation perspective. What we should do from a PROTEUS point-of-view is enabling arbitrary accretion paths (i.e. addition/loss of H-C-N-S and fO2) by read-in from somewhere else (a module, or data/text file that is produced from a different model).

nichollsh commented 2 months ago

Do we have something like an order of magnitude estimate for the duration and rate of accretion?

Handling this with a module is a nice idea. Maybe this could even go into ZEPHYRUS?

timlichtenberg commented 2 months ago

Obviously I am making it complicated, but a simpler initial approach may be useful. One could, for example, just start with some volatile addition rate and fO2/time change and then add a splitting term for H-C-N-S fractionation.

timlichtenberg commented 2 months ago

The volatile accretion can be exponential or linear for example as a start.

timlichtenberg commented 2 months ago

And last comment: I would prefer a different module for accretion. Like escape, volatile and mass accretion is a near-boundless area of research in itself, and such a module would have substantial potential to grow massively over time.

nichollsh commented 2 months ago

Seems like there's a large number of ways in which to take this, so a separate model makes sense. Adding initial support into PROTEUS with a fixed rate (or choosing linear, exponential as you suggest) doesn't sound too tricky.

Now we just need someone to do this! I might have a look at making a stub module for this in the near future, but my physics knowledge on this is lacking.

nichollsh commented 2 months ago

Created an initial stub repo here: https://github.com/FormingWorlds/ABASSI

timlichtenberg commented 2 months ago

Nice, good name. Eq. 1 in Sturtz et al. (2022, Icarus) is a reasonable starting point for an adaptable continuous accretion law. Figure 6 of Johansen et al. (2023, A&A) showcases the pebble-only accretion hypothesis for the Solar System terrestrial planets. These are mass accretion laws, however. I will think about how volatile addition paths can be parameterized or where to start from.

Re: names: Next time there is a choice on a module name we may want to choose a female deity, per suggestion of @EmmaPostolec. There are many dimensions of diversity, however, so adding an African deity is also a good addition.

nichollsh commented 2 months ago

See also this preprint https://arxiv.org/abs/2409.13820

nichollsh commented 3 weeks ago

https://www.nature.com/articles/s41561-021-00835-9

timlichtenberg commented 3 weeks ago

Well well well... 🙃 https://www.nature.com/articles/nature21045 https://iopscience.iop.org/article/10.3847/2041-8213/ac9521/meta It's a complex and unresolved topic.