Modeling Non-Rotational Ocean Circulation and Heat Distribution in Icy Moons

Subsurface oceans beneath the ice shells of icy moons like Europa and Enceladus are considered promising environments for extraterrestrial life. Their long-term habitability depends on internal heating and efficient vertical heat transport to maintain liquid water beneath the surface. This study models vertical heat diffusion in a non-rotating ocean column to investigate thermal structure and energy balance in such systems. A one-dimensional numerical simulation was developed using temperature-dependent thermal conductivity and fixed Dirichlet boundary conditions, initialized with a linear temperature gradient from −10 K at the surface to +10 K at the base. Over 1000 time steps, the temperature profile became nonlinear, with a kink indicating the transition from ice to water. Despite fixed boundary temperatures, the interior warmed, and the average temperature rose to 2.84 K. This resulted from asymmetric conductivity: efficient heating from below and slow heat loss through the upper ice. These results illustrate how conductivity structure controls thermal evolution and ice shell stability on ocean worlds.

Name thesis:degree_name

Bachelor

Department dc:contributor.department

Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences

Grantor dc:publisher

Massachusetts Institute of Technology

Year dc:date.issued

2025

Author dc:creator

• Nath, Anika

Advisors dc:contributor.advisor

• Teague, Richard
• Hill, Christopher N.

Statement dc:rights

• In Copyright - Educational Use Permitted
• Copyright retained by author(s)

Licence dc:rights.uri

• https://rightsstatements.org/page/InC-EDU/1.0/

Handle dc:identifier.uri

https://hdl.handle.net/1721.1/159901

OAI identifier oai:identifier

oai:dspace.mit.edu:1721.1/159901