Our serpentinization model is based on [1]} which is an improved serpentinization model from heat balance [2]}, [3]} and dynamic models [4]}. [1]} used a kinetic approach to estimate the rate of serpentinization and follow its evolution.
Their model implemented interfacial water but included some simplifications and neglected several effects that did not cause a significant difference in the results for the smaller object (\(\sim \) 15 km radius), low water ice content and low porosity they studied.
In our cases, however, these neglected effects are important as we aim to consider larger planetesimals and icy bodies as well.
In the description of the serpentinization model (in Sect. REF . and REF .), we present which effects were taken into account in the current model as an improvement, compared with the previous model in which these parameters were simplified.
