Data from: An oceanic mechanism for geyser formation on Enceladus
Data files
Jul 28, 2026 version files 225.22 MB
-
extract_power.ipynb
128.04 KB
-
gaussian_1km_dataset.txt
645 B
-
gaussian_5km_dataset.txt
645 B
-
power_dent_H1_U2.data
12.29 MB
-
power_dent_H1_U2.meta
456 B
-
power_dent_H1_U4.data
12.29 MB
-
power_dent_H1_U4.meta
456 B
-
power_dent_H1_U6.data
12.29 MB
-
power_dent_H1_U6.meta
456 B
-
power_dent_H3_U2.data
12.29 MB
-
power_dent_H3_U2.meta
456 B
-
power_dent_H3_U4.data
12.29 MB
-
power_dent_H3_U4.meta
456 B
-
power_dent_H3_U6.data
12.29 MB
-
power_dent_H3_U6.meta
456 B
-
power_dent_H5_U2.data
12.29 MB
-
power_dent_H5_U2.meta
456 B
-
power_dent_H5_U4.data
12.29 MB
-
power_dent_H5_U4.meta
456 B
-
power_dent_H5_U6.data
12.29 MB
-
power_dent_H5_U6.meta
456 B
-
power_ridge_H1_U2.data
12.29 MB
-
power_ridge_H1_U2.meta
456 B
-
power_ridge_H1_U4.data
12.29 MB
-
power_ridge_H1_U4.meta
456 B
-
power_ridge_H1_U6.data
12.29 MB
-
power_ridge_H1_U6.meta
456 B
-
power_ridge_H3_U2.data
12.29 MB
-
power_ridge_H3_U2.meta
456 B
-
power_ridge_H3_U4.data
12.29 MB
-
power_ridge_H3_U4.meta
456 B
-
power_ridge_H3_U6.data
12.29 MB
-
power_ridge_H3_U6.meta
456 B
-
power_ridge_H5_U2.data
12.29 MB
-
power_ridge_H5_U2.meta
456 B
-
power_ridge_H5_U4.data
12.29 MB
-
power_ridge_H5_U4.meta
456 B
-
power_ridge_H5_U6.data
12.29 MB
-
power_ridge_H5_U6.meta
456 B
-
produce_fig2.ipynb
128.01 KB
-
produce_fig3.ipynb
6.89 KB
-
produce_fig4.ipynb
96 KB
-
produce_fig5.ipynb
201.01 KB
-
produce_fig6.ipynb
230.05 KB
-
README.md
7.40 KB
-
velocity_H40.data
1.61 MB
-
velocity_H40.meta
183 B
-
velocity_H60.data
1.61 MB
-
velocity_H60.meta
183 B
Abstract
The tiger stripes of Enceladus are a sequence of four parallel and evenly spaced fissures 130 km long that erupt water from a global subsurface ocean into space. We present an ocean dynamical theory of their formation. Given one initial fissure, topography on the underside of the ice shell rubs periodically against the ocean due to the satellite's libration. This motion excites internal gravito-inertial waves which propagate down, reflecting at the seafloor, before reimpacting the underside of the shell. Wave breaking upon impact generates thermal energy which promotes melting via a feedback mechanism, forming subsequent fissures in a chain of parallel stripes. Dissipation predicted by linear analytical theory matches nonlinear simulations using the MIT General Circulation Model, with stronger agreement when the topography is less steep. Sufficiently large topography and libration rate excite waves with enough energy to initiate the fissure formation process. For south polar ocean thickness between 30 and 60 km, we constrain the stratification (N < 4.6 x 10-5s-1) at which the resulting wave pattern replicates the observed 35 km spacing of the tiger stripes.
Description of the data and file structure
Files and variables
Python notebooks have been provided for reproducing Figures 2-6 in the main text (named produce_fig#.ipynb).: produce_fig2.ipynb, produce_fig3.ipynb, produce_fig4.ipynb, produce_fig5.ipynb, and produce_fig6.ipynb
The corresponding data for reproducing the non-analytical figures (3 and 4) has been provided as well. In order to read these data files directly, one must have the MITgcm model installed (https://mitgcm.readthedocs.io/) and, in particular, the MITgcmutils directory that comes along with it, providing the functions that read the data. The scripts here rely on having the MITgcm installed and provide example code for reading in the data.
For Figure 4 in particular, use extract_power.ipynb to get the estimated local dissipation from the data and to produce Fig 4a. This value is used to populate the P_local array in produce_fig4.ipynb, which generates Fig 4b and Fig 4c.
The files gaussian_1km_dataset.txt and gaussian_5km_dataset.txt were generated by a finite element model of the ice shell and are used in produce_fig5.ipynb.
For additional data or inquiries, contact dabdulah@mit.edu.
File: power_dent_H1_U2.meta
Description: Diagnosed power for an IGW excited by topography of height 1 km with libration amplitude 2.6 mm/s
File: power_dent_H1_U4.meta
Description: Diagnosed power for an IGW excited by topography of height 1 km with libration amplitude 4.6 mm/s
File: power_dent_H1_U6.meta
Description: Diagnosed power for an IGW excited by topography of height 1 km with libration amplitude 6.5 mm/s
File: power_dent_H3_U2.meta
Description: Diagnosed power for an IGW excited by topography of height 3 km with libration amplitude 2.6 mm/s
File: power_dent_H1_U6.data
Description: Diagnosed power for an IGW excited by topography of height 1 km with libration amplitude 6.5 mm/s
File: power_dent_H1_U2.data
Description: Diagnosed power for an IGW excited by topography of height 1 km with libration amplitude 2.6 mm/s
File: power_dent_H1_U4.data
Description: Diagnosed power for an IGW excited by topography of height 1 km with libration amplitude 4.6 mm/s
File: power_dent_H3_U4.meta
Description: Diagnosed power for an IGW excited by topography of height 3 km with libration amplitude 4.6 mm/s
File: power_dent_H3_U6.meta
Description: Diagnosed power for an IGW excited by topography of height 3 km with libration amplitude 6.5 mm/s
File: power_dent_H3_U2.data
Description: Diagnosed power for an IGW excited by topography of height 3 km with libration amplitude 2.6 mm/s
File: power_dent_H3_U4.data
Description: Diagnosed power for an IGW excited by topography of height 3 km with libration amplitude 4.6 mm/s
File: power_dent_H5_U2.meta
Description: Diagnosed power for an IGW excited by topography of height 5 km with libration amplitude 2.6 mm/s
File: power_dent_H5_U4.meta
Description: Diagnosed power for an IGW excited by topography of height 5 km with libration amplitude 4.6 mm/s
File: power_dent_H5_U2.data
Description: Diagnosed power for an IGW excited by topography of height 5 km with libration amplitude 2.6 mm/s
File: power_dent_H5_U6.meta
Description: Diagnosed power for an IGW excited by topography of height 5 km with libration amplitude 6.5 mm/s
File: power_dent_H3_U6.data
Description: Diagnosed power for an IGW excited by topography of height 3 km with libration amplitude 6.5 mm/s
File: power_ridge_H1_U2.meta
Description: Diagnosed power for an IGW excited by topography of height 1 km with libration amplitude 2.6 mm/s
File: power_ridge_H1_U4.meta
Description: Diagnosed power for an IGW excited by topography of height 1 km with libration amplitude 4.6 mm/s
File: power_dent_H5_U6.data
Description: Diagnosed power for an IGW excited by topography of height 5 km with libration amplitude 6.5 mm/s
File: power_dent_H5_U4.data
Description: Diagnosed power for an IGW excited by topography of height 5 km with libration amplitude 4.6 mm/s
File: power_ridge_H1_U6.meta
Description: Diagnosed power for an IGW excited by topography of height 1 km with libration amplitude 6.5 mm/s
File: power_ridge_H3_U2.meta
Description: Diagnosed power for an IGW excited by topography of height 3 km with libration amplitude 2.6 mm/s
File: power_ridge_H1_U4.data
Description: Diagnosed power for an IGW excited by topography of height 1 km with libration amplitude 4.6 mm/s
File: power_ridge_H3_U4.meta
Description: Diagnosed power for an IGW excited by topography of height 3 km with libration amplitude 4.6 mm/s
File: power_ridge_H3_U6.meta
Description: Diagnosed power for an IGW excited by topography of height 3 km with libration amplitude 6.5 mm/s
File: power_ridge_H5_U2.meta
Description: Diagnosed power for an IGW excited by topography of height 5 km with libration amplitude 2.6 mm/s
File: power_ridge_H1_U2.data
Description: Diagnosed power for an IGW excited by topography of height 1 km with libration amplitude 2.6 mm/s
File: power_ridge_H1_U6.data
Description: Diagnosed power for an IGW excited by topography of height 1 km with libration amplitude 6.5 mm/s
File: power_ridge_H3_U6.data
Description: Diagnosed power for an IGW excited by topography of height 3 km with libration amplitude 6.5 mm/s
File: power_ridge_H3_U2.data
Description: Diagnosed power for an IGW excited by topography of height 3 km with libration amplitude 2.6 mm/s
File: power_ridge_H3_U4.data
Description: Diagnosed power for an IGW excited by topography of height 3 km with libration amplitude 4.6 mm/s
File: power_ridge_H5_U4.meta
Description: Diagnosed power for an IGW excited by topography of height 5 km with libration amplitude 4.6 mm/s
File: power_ridge_H5_U6.meta
Description: Diagnosed power for an IGW excited by topography of height 5 km with libration amplitude 6.5 mm/s
File: power_ridge_H5_U2.data
Description: Diagnosed power for an IGW excited by topography of height 5 km with libration amplitude 2.6 mm/s
File: power_ridge_H5_U4.data
Description: Diagnosed power for an IGW excited by topography of height 5 km with libration amplitude 4.6 mm/s
File: power_ridge_H5_U6.data
Description: Diagnosed power for an IGW excited by topography of height 5 km with libration amplitude 6.5 mm/s
File: velocity_H40.meta
Description: Vertical velocity for the excited IGW in a 36 km domain, N2 = 1e-11
File: velocity_H60.meta
Description: Vertical velocity for the excited IGW in a 57 km domain, N2 = 1.9e-9
File: velocity_H60.data
Description: Vertical velocity for the excited IGW in a 57 km domain, N2 = 1.9e-9
File: velocity_H40.data
Description: Vertical velocity for the excited IGW in a 36 km domain, N2 = 1e-11
Code/software
The MITgcmutils folder enables reading the .data and .meta files provided. The plotting files provided are Python notebooks.
