Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published April 29, 2024 | in press
Journal Article Open

Jet activity on Enceladus linked to tidally driven strike-slip motion along tiger stripes

  • 1. ROR icon California Institute of Technology

Abstract

At Saturn’s moon Enceladus, jets along four distinct fractures called ‘tiger stripes’ erupt ice crystals into a broad plume above the South Pole. The tiger stripes experience variations in tidally driven shear and normal traction as Enceladus orbits Saturn. Here, we use numerical finite-element modelling of a spherical ice shell subjected to tidal forces to show that this traction may produce quasi-periodic strike-slip motion in the Enceladus crust with two peaks in activity during each orbit. We suggest that friction modulates the response of tiger stripes to driving stresses, such that tidal traction on the faults results in a difference in the magnitudes of peak strike slip and delays the first peak in fault motion following peak tidal stress. The simulated double-peaked and asymmetric strike-slip motion of the tiger stripes is consistent with diurnal variations in jet activity inferred from Cassini spacecraft images of plume brightness. The spatial distribution of strike-slip motion also matches Cassini infrared observations of heat flow. We hypothesize that strike-slip motion can extend transtensional bends (for example, pull-apart structures) along geometric irregularities over the tiger stripes and thus modulate jet activity. Tidally driven fault motion may also influence longer term tectonic evolution near the South Pole of the satellite.

Copyright and License

© The Author(s), under exclusive licence to Springer Nature Limited 2024.

Acknowledgement

This research was supported by the Future Investigators in NASA Earth and Space Science and Technology (FINESST) Program (80NSSC22K1318)(A.B., M.S.). We thank the Keck Institute for Space Studies (KISS) at the California Institute of Technology for organizing two workshops about ‘Next-Generation Planetary Geodesy’ which provided insight, expertise and discussions that inspired this research. We also thank M. Knepley, B. Aagaard and C. Williams for providing valuable advice on how to modify PyLith for our simulations. A portion of this research was supported by a Strategic Research and Technology Development task led by J. T. Keane and R. S. Park at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (80NM0018D0004)(J.T.K., R.S.P.).

Contributions

A.B. conceived and designed this study under the supervision of M.S. A.B. drafted the article and constructed Figs. 14. M.S., J.T.K. and A.B. developed numerical models used in the study. R.S.P. provided the shape model necessary for numerical simulations. J.T.K., with the aid of A.B., constructed Fig. 5. E.J.L. provided expertise regarding the geological evolution of the SPT. All authors discussed the results of the study and commented on the manuscript at each stage of revision.

Data Availability

Data used for Fig. 3 (slab densities derived from Cassini ISS images) is publicly available via https://doi.org/10.1016/j.icarus.2019.06.006 (ref. 7). Data used for Fig. 4 (heat flow along the tiger stripe faults derived from Cassini Composite Infrared Spectrometer measurements) is publicly available via https://doi.org/10.2458/azu_uapress_9780816537075-ch008 (ref. 4).

Extended Data Fig. 1 Similar to upper right and upper left panels of Fig. 2 of the main text except values (driving tractions) are evaluated for a model with slipping faults (μ = 0.4).

Extended Data Fig. 2 Example snapshots of mesh geometry.

 

Extended Data Fig. 3 Spin-up parameter Ξ(t) (see Equation (21)) modelled as a function of time (in units of tidal periods) for several prescribed values for the static coefficient of friction μ.

Extended Data Fig. 4 Simplified model relating driving shear traction τD, normal traction σn, μ, and slip s along the tiger stripe faults.

Extended Data Fig. 5 Correlation of timing of lateral slip and plume brightness for several values of modelled μ7.

Extended Data Fig. 6 Crustal thickness assumed for finite-element models.

Extended Data Fig. 7 Correlation of spatial distribution of lateral slip and radiated power per unit length4 for several values of modelled μ.

 

Supplementary Table 1

Supplementary Video

Code Availability

The results used in this study were generated using the software package PyLith44,50. PyLith is an open-source finite-element code for modelling geodynamic processes and is available via Zenodo at https://doi.org/10.5281/zenodo.3269486 (ref. 50). The specific PyLith version used in this study was v2.2.2. The full code used to simulate deformation in this work (including modifications Pylith, CUBIT and PETSc and a user manual) are publicly available via https://doi.org/10.5281/zenodo.10585162 (ref. 51). The mesh geometries utilized in this study were created using CUBIT (v15.2), a node-locked licensed software, which is available through the developer Sandia National Laboratories via https://cubit.sandia.gov/2021/09/14/cubit-15-7-released-october-20-2020/ (ref. 48).

Conflict of Interest

The authors declare no competing interests.

Files

41561_2024_1418_Fig9_ESM.jpg
Files (18.0 MB)
Name Size Download all
md5:f480b02ee59b5048bfa90458bec0e9b3
89.2 kB Preview Download
md5:3b8b66cafc7cc94ab5903e01d63a6ded
223.7 kB Preview Download
md5:a87cd04087e78624e504347d2a0f339b
117.1 kB Preview Download
md5:6b8fb96d61dcfb13812910696c383ce7
16.0 MB Download
md5:133bb04efb04d055a87e7a1d3a4c4869
260.7 kB Preview Download
md5:6b5e1065306a65509253192111188bc3
191.7 kB Preview Download
md5:80cd3e87a68c08dde8d9b91939f788f0
344.8 kB Preview Download
md5:bd5b039445e02cbbce7c4d3e7f915247
117.8 kB Preview Download
md5:7fd543553146d4f8810240feb169f858
710.1 kB Preview Download

Additional details

Created:
April 30, 2024
Modified:
April 30, 2024