| Sumario: | 25 pages, 14 pages, supporting information https://doi.org/10.1029/2025JB031510.-- Data Availability Statement: All the model results, data processing and plotting scripts relevant to this study have been deposited in the public repository of Zenodo at Zhang and Li (2025) (https://doi.org/10.5281/zenodo.14969626). For visualizing the 3D model results, we utilized the open source software Paraview (Ahrens et al., 2005), which can be downloaded from https://www.paraview.org. Figures 2, 10 and 11a were crafted using Adobe Illustrator (https://www.adobe.com/products/illustrator.html). Figure 14 was generated with GMT6 (Wessel et al., 2019) obtained from https://www.generic-mapping-tools.org. The remaining figures were created and composed using Matplotlib (Hunter, 2007) from https://matplotlib.org, with scientific colormaps (Fabio Crameri et al., 2020) sourced from https://www.fabiocrameri.ch/colourmaps/. Curve fitting, optimization, special function calculations, and numerical integration were carried out using SciPy (Virtanen et al., 2020) from https://scipy.org. Array computations and data storage were managed by NumPy from https://numpy.org (Harris et al., 2020). The computer algebra software Wolfram Mathematica (https://www.wolfram.com/mathematica/) was employed to conduct the analytical solution of Equations 15a-16. The I2VIS code could be found in https://doi.org/10.17605/OSF.IO/bnvth (Gerya et al., 2021), and the I3VIS code, with all methods published in Gerya (2019), can be directly requested from the main code developer (taras.gerya@erdw.ethz.ch) or the corresponding author (li.zhonghai@ucas.ac.cn)
|