# How to calculate the thickness of the ice-layer of a frozen ocean planet/moon?

How do I calculate the thickness of the upper ice layer on cold ocean worlds like Europa, Enceladus, Ganymede, ...? I'm asking this for a programm I'm currently writing.

Given/Known is: mass, radius, average density of metal (assume Fe) core, rock (assume MgSiO3) layer, ice (assume H2O) layer; internal heat from tidal heating, radioactive decay and residual formation heat; external heating from solar radiation (assume equal constant heating distributed over the entire surface or ignore it if it is irrelevant/minuscule); pressure of the atmosphere above (assume pure N2 with variable pressure or no atmosphere)

Wanted is: thickness frozen (crust) and molten layer

Bonus stuff: If easy adaptation of the formula for the calculation of the rocky layers crust thickness is possible I would appreciate it. Should it be possible, an adaptation to find the rocky crust thickness and the ice crust thickness on wolds where both are needed would be helpful.

Thanks for any answer and time spent of solving this.

This question asks for hard science. All answers to this question should be backed up by equations, empirical evidence, scientific papers, other citations, etc. Answers that do not satisfy this requirement might be removed. See the tag description for more information.

• I honestly think this might be more suited to Astronomy's site. You might find someone here who can answer it, but it seems to be a real-world question, not about worldbuilding. If you can show us how it relates to your world that might at least stop the question getting closed. – Measure of despare. Mar 25 at 18:33
• @Agrajag It appears to me that the question is specifically about building planets (worlds). I don't think it could be any more on topic. – Arkenstein XII Mar 25 at 19:49
• Are you looking for the dynamic or the steady state equation? – L.Dutch Mar 25 at 19:52
• @ArkensteinXII Well not about a world at least, or is the "How do I balance 15 variables I'm not going to define, all of which have an undefined error bar?" sort of question on topic, d'y' rekon? – Measure of despare. Mar 25 at 19:53
• I usually do the first one when I move questions, although this time there's no need to delete the original. – SealBoi Mar 25 at 21:51

For determining the ice layer thickness you can use a simplified 1D model, as follows:

• You want the ice to be present above a layer of water. This sets the temperature $$T_b$$ at the ice-water interface. This would be 0 Celsius for pure water a 1 bar.
• You know the thermal flow coming from the core, $$Q$$
• you want the condition to be stationary, this means that the flow is entirely dissipated to the outer environment. Therefore you have that $$Q=\lambda \cdot A \cdot (T_b - T_a)\over d$$, where $$d$$ is the thickness of ice layer, $$\lambda$$ is the thermal conductivity of ice, $$A$$ the surface and $$T_a$$ the surface temperature.

If you have the possibility of setting the other interface temperatures as well, you can use the same equation to determine the thicknesses of all the layers.

Mind that this model can be used for a spherical shell as long as the thickness is small with respect to the radius of the shell. Else a spherical model would be more appropriate.

In that case it would be, assuming only radial flux, $$Q=4 \pi \lambda R_a \cdot R_bT_a - T_b \over R_b - R_a$$

• Not knowing the thermal flow from the core, you cannot evaluate the thickness of ice, because for its counting you need thermal flow through the lower and upper surfaces, (you know only one of them), and their temperatures (about 0 and the measured one). – Gangnus Mar 28 at 8:06

While I'm sure there's hard-set rules, and someone with a degree in astrophysics may answer this question more thoroughly, I think you'll find that theres an incredible amount of factors to something like this.

Take Europa for example.

We estimate anywhere between 19 to 25 kilometers of ice, but mind you. this is an estimate based on a bunch of data we have compiled over time, and even then, the estimate is only down to 6km, which, all-in-all is a big distance (not on an astrophysical scale, but you get what I mean).

We assume there's quite a bit of Liquid water under Europa's surface because of observed ejection of liquid and also how impact craters look/are formed. We also know that if there is liquid water, it is kept liquid simply due to Jupiter's gravitation [citation needed]. We can measure its size and can infer to how much water there is around it, regardless of state. But even in the real world, we come up short answering hard-set questions. We don't know how much not-water the inside of the moon is.

We also have evidence that the ice sheet around Europa might have been different thicknesses throughout its life.

And then we still haven't talked about tectonics or other geological events.

My advice? Let the user of your app, be it you, or anyone else, define themselves, how much ice they want, it's probably better that way, and if someone's worlds have different rules, then your app can accommodate to that by not setting any hard-set rules.

Interested to see what the app looks like/what it will be able to do though!

• Thes goal is to be realistic. I'm not willing to go down the quick-and-drity approximation road just yet. – TheDyingOfLight Mar 25 at 20:45
• That's fair, and I admire you wanting to do this as realistic as possible. I'd echo the other commenters and say to post it in Astronomy instead. You can either delete the question here or ask to close it. – Cristian C. Mar 25 at 20:57

The layer under ice is the layer of possibly dirty water. Really, if you have tidal and residual warmth and temperature on the surface, "all" you need is "merely" the thermal conductivity of the ice and that of the water layer. The problem is, that the thermal conductivity of water is defined by convection cells structure and is unknown. And it cannot be known without inner observations. You cannot guess it theoretically. The contemporary science cannot predict the behaviour of the turbulent flows. Not only counting by formulae is impossible yet, but even the computer modelling. It is simply one of the unsolved tasks.

Edit. The situation is even worse. Here (https://en.wikipedia.org/wiki/Turbulence_modeling) you can see the maximum where the recent science reached. Only statistical evaluation of the evolution of flows. So, maximally you can guess how the structure of flows could evolve. And that structure is dynamic, and because of that the thickness of the ice layer will dynamically change, and again it will influence that turbulence structure. So, you cannot say what the thickness is, because it simply changes.

Using the equations from the reference, you could guess the dynamical range of the changes. But, I am afraid, that is a great work for a doctor of Sc. degree.

This question asks for hard science. All answers to this question should be backed up by equations, empirical evidence, scientific papers, other citations, etc. Answers that do not satisfy this requirement might be removed. See the tag description for more information.

• The question is tagged hard science. All answers to this question should be backed up by equations, empirical evidence, scientific papers, other citations, etc. Answers that do not satisfy this requirement might be removed. – L.Dutch Mar 26 at 6:29
• @L.Dutch I think the answer "there are no known equations to describe this" is a valid one in this case. After all, it's hard to provide hard science for a topic of which there is no science (yet). – Rekesoft Mar 26 at 9:30
• @Rekesoft, that might be true for the liquid layer. However, OP is asking for the thickness of the solid layer, and also the statement "this is not possible" should be backed up by appropriate citation/reference. Mind, I am not discussing the content, but the form, due to the tag involved. – L.Dutch Mar 28 at 7:32
• @L.Dutch I had changed the question, so, formally your requirement is fulfilled. But really, the older variant of the post, equal to the first paragraph, did not need reference, for the information mentioned there belongs to the middle school program. We don't support 2*2=4 by references. It is enough widely known information to be checked by public knowledge. – Gangnus Mar 28 at 8:00