Datasets for "Europa's ocean translates interior tidal heating patterns to the ice-ocean boundary" by D. G. Lemasquerier, C. J. Bierson and K. M. Soderlund

heatflux-velocities.zip Heat flux, radial velocity, and azimuthal velocity for the simulations presented in the main text (Figure 3), i.e. - 6 simulations with stress-free boundary conditions (folder "NeumannSF") - 6 simulations with no slip boundary conditions (folder "NeumannNS") Each folder contains the 6 simulations corresponding to 6 \(q^*\) values, \(q^*\)=[0.00,0.46,0.87,1.30,1.72,2.16]. For each simulation, we provide the data for the heat fluxes at the top and bottom boundaries, and the radial and azimuthal velocities in the ocean midplane. - "dsTdr_bottom.txt" and "dsTdr_top.txt" contain heat flux data at the bottom and top boundaries respectively. Each 2D map is an array of size Ntheta*Nphi where Ntheta is the number of grid points in latitude (-90° to 90°) and Nphi=2*Ntheta is the number of grid points in longitude (-180° to 180°). In the present files, this array is stored in column-major order. For instance, the first Ntheta lines correspond to latitudes from -90° to 90°, and a fixed longitude of -180°. The next Ntheta lines correspond to the next longitude. The heat flux is given in units of the absolute value of the average bottom heat flux. - "vr_mid.txt" and "vp_mid.txt" contain radial and azimuthal velocity data at mid-height in the ocean. The ordering is the same as described above. Each 2D map is an array of size Ntheta*Nphi where Ntheta is the number of grid points in latitude (-90° to 90°) and Nphi=2*Ntheta is the number of grid points in longitude (-180° to 180°). In the present files, this array is stored in column-major order. Velocities are provided in dimensionless units as a Reynolds number, i.e. they are normalized by \(\nu/D\) where \(\nu\) is the kinematic viscosity of the fluid and \(D\) the spherical shell thickness. icethickness.zip Matlab files for the ice thickness model input and outputs for four cases. The filenames refer to the simulation from which the relative heat flux variations were extracted: - "NeumannSF_q0.mat": Figure 5(c), and Figure C2(c) (Neumann, q*=0.0, stress-free) - "NeumannSF_q087.mat": Figure 5(f), and Figure C2(i) (Neumann, q*=0.87, stress-free) - "NeumannNS_q0.mat": Figure C2(f) (Neumann, q*=0.0, no-slip) - "NeumannNS_q087.mat": Figure C2(l) (Neumann, q*=0.87, no-slip) Each file contains: 1) The ice thickness model inputs (Figure C1): - Lambda: 100*50 array, latitude in radians. - Phi: 100*50 array, longitude in radians. - epsdot2avg: 100*50 array (longitude*latitude), quantity \(\overline{\dot{\epsilon}_{ij}^2}\) in \(s^{-2}\), where \(\epsilon_{ij}\) is the tidal strain within the ice. Computed from Ojakangas and Stevenson (1989), Appendix B. - Fbmat: 100*50 array, oceanic heat flux in \(W/m^2\) . - Ts: 100*50 array, surface temperature in K, computed from Ojakangas and Stevenson (1989), Appendix A. 2) The ice thickness model outputs: - D: 100*50 array, ice thickness in meters. - etafinal: 100*50*50 array (longitude*latitude*depth), ice viscosity in Pa s. The third dimension is the depth, which can be represented by a linearly spaced vector going from 0 to D(i,j). - Hfinal: 100*50*50 array (longitude*latitude*depth), volumetric heating in the ice in W/m^3. - TT: 100*50*50 array (longitude*latitude*depth), temperature in the ice in K.