Seismic wave propagation in icy ocean worlds
Simon C. Stähler, Mark P. Panning, Steve D. Vance, Ralph Lorenz, Martin van Driel, Tarje Nissen-Meyer, Sharon Kedar
SSeismic wave propagation in icy ocean worlds
Simon C. St¨ahler ∗† , Mark P. Panning ‡ , Steve D. Vance ‡ ,Ralph D. Lorenz § , Martin van Driel ∗ ,Tarje Nissen-Meyer ¶ , Sharon Kedar ‡ November 20, 2017
Seismology was developed on Earth and shaped our model of the Earth’sinterior over the 20th century. With the exception of the Philae lander, allin situ extraterrestrial seismological effort to date was limited to other ter-restrial planets. All have in common a rigid crust above a solid mantle. Thecoming years may see the installation of seismometers on Europa, Titan andEnceladus, so it is necessary to adapt seismological concepts to the settingof worlds with global oceans covered in ice. Here we use waveform analy-ses to identify and classify wave types, developing a lexicon for icy oceanworld seismology intended to be useful to both seismologists and planetaryscientists. We use results from spectral-element simulations of broadbandseismic wavefields to adapt seismological concepts to icy ocean worlds. Wepresent a concise naming scheme for seismic waves and an overview of thefeatures of the seismic wavefield on Europa, Titan, Ganymede and Ence-ladus. In close connection with geophysical interior models, we analyzesimulated seismic measurements of Europa and Titan that might be usedto constrain geochemical parameters governing the habitability of a sub-iceocean.
1. Introduction
Most of the larger icy moons of the solar system’s gas planets are considered to harborliquid oceans below their surface (Nimmo and Pappalardo, 2016). These oceans have ∗ ETH Z¨urich, Institute for Geophysics † Leibniz-Institute for Baltic Sea Research, Rostock, Germany ‡ Jet Propulsion Laboratory, California Institute of Technology, Pasadena, USA § Johns Hopkins University Applied Physics Laboratory, Laurel, USA ¶ Department of Earth Sciences, University of Oxford, United Kingdom a r X i v : . [ phy s i c s . g e o - ph ] N ov igure 1: Comparison of a measured seismograms on Earth (Earthquake in the LakeTanganyika region on 2016/02/24, 32 km depth; measured at BFO in South-ern Germany in 60 degree epicentral distance) with synthetic seismogramsfor Earth, Enceladus, Europa, Ganymede and Titan with equivalent angularsource-receiver distance. The source model is the real moment tensor, with asource duration of 2 seconds. For the icy moons, the depth has been reducedto 3 km, to ensure that the event is in the ice.2een predicted from the estimated internal energy budget by radioactive decay andtides (Lewis, 1971; Cassen et al., 1979) and are very likely to form in a multi-moonsituation around a giant planet. In specific settings, even small worlds like Enceladuscan have enough energy available to sustain a liquid ocean at least temporarily. Theexistence of oceans at Europa, Titan and Enceladus has been identified using densitystructure inferred from gravitational moment-of-inertia measurements (Schubert et al.,2004), tectonic interpretations of the ice thickness, and electromagnetic effects consis-tent with conductive fluids beneath the ice (Kivelson et al., 2000; B´eghin et al., 2010).Due to the low density contrast between ice and water and ambiguity of electrical con-ductance measurements, these methods alone cannot constrain the thicknesses of iceand ocean layers uniquely, which are crucial to the ocean chemistry and therefore hab-itability. Especially in the large worlds Ganymede, Callisto and Titan, high-pressureice phases at the ocean bottom would could limit water-rock interaction which regulateredox and ocean composition(Vance et al., 2014, 2017b), with a potentially negativeeffect on habitability. Briny fluids (Hogenboom et al., 1995) or ice convection mightmitigate this effect (Choblet et al., 2017), nevertheless the presence, structure anddynamics of these ices is therefore of great scientific interest.The thickness of a surface ice layer less than a few kilometers can be estimated fromradar measurements, as planned for the upcoming ESA JUICE and NASA EuropaClipper missions (Bruzzone et al., 2013; Phillips and Pappalardo, 2014). Attenuationin the water will make it difficult to measure the ocean depth with radar; seismo-logical measurements would complement the existing methods. Where radar is sen-sitive to electrical conductivity, and crystal orientation fabric (COF), seismic wavesare influenced by elastic parameters and indirectly by material chemistry, COF, andtemperature (Diez, 2013). A detailed discussion of geophysical measurements linkedto habitability can be found in the companion paper by Vance et al. (2017b).The seismology of an icy ocean world differs from that of Earth. Many of the classicalconcepts and terminology of seismology become meaningless with the addition of athick global ocean and ice covering. The global ocean, with its low seismic velocityand absence of S-waves, partially decouples the seismic wavefield in the ice and theunderlying solid interior, creating seismic wave types that are not observed at a globalscale on Earth. The coming decades will likely bring the first seismological station onan icy moon [e.g.][](Hand et al., 2017), and mission design begins now. This paper givesan introduction to the general features of icy ocean world seismology, which shouldhelp terrestrial seismologists adapt their vocabulary and thinking to this environment.With the possibility of landers on Europa and Titan in the coming decades, these twoworlds are discussed in more detail. Enceladus is also discussed, even though thereis currently no strong proposal for a lander and its small, irregular size distinguish itfrom the other ocean worlds. Ganymede is also discussed as a target for a possiblefuture lander (Vijendran et al., 2010), and as end-member in size and ocean thickness.To keep the paper concise, we will not discuss Callisto and potential ocean worlds inthe Uranus and Neptune systems in detail, even though the general concepts presentedhere apply to them as well. For the four worlds, potential seismic tests for habitabilityare presented and tested using synthetic seismograms.Previous studies of seismology on Europa, mostly older than a decade, were therefore3imited by the computational methods of the time. The detailed study of (Kovachand Chyba, 2001) was based on ray-theoretical methods and analytic derivations anddid not calculate any full seismograms. Lee et al. (2003) used ray theory to discussbody waves and surface waves in the ice and ocean layer of Europa, Panning et al.(2006) used the normal mode method (Woodhouse, 1988) to calculate seismograms,but were limited to long periods (¿10s), where some of the ice layer effects are notevident. Computational seismology has progressed rapidly in the last decade (for anoverview, see Igel, 2017), enabling full modeling of the broadband global wavefieldfor the first time now. We used the spectral-element solver AxiSEM (Nissen-Meyeret al., 2014), which is extremely efficient for layered models, since it reduces the 3Dwave propagation to a 2D problem. A second companion paper (Panning et al., 2017),couples these results with a seismicity model to estimate the seismic background noisecaused by tidal cracking on Europa.Any realistic lander mission on an icy world will face severe constraints in terms ofseismometer mass: Generally, the noise level of a seismometer is controlled by the sizeof its proof mass: Smaller masses mean higher eigenfrequencies (above 1 Hz) and ahigher noise level at long periods. Classical broadband seismometers, like the InSightvery Broadband (VBB) main instrument (Lognonne et al., 2012; Mimoun et al., 2012;Lognonn´e et al., 2016), have a very good signal-to-noise ratio over frequency rangefrom below 0.01 Hz to 10 Hz, but have total masses of several kg. Also, their sensi-tivity to shock and tilted installation impacts mission planning, which is prohibitiveas a secondary payload on a lander. Geophones, as proposed for the Mars96 mission(Lognonn´e et al., 1998) are very lightweight, passive instruments, but have a low sen-sitivity below one Hz and will only record local events and structure. The SESAMEinstrument on the Philae lander (Seidensticker et al., 2007) used piezoelectric sensors,sensitive in the mHz range to see structure on a meter scale around the landing siteon 67P/Churyumov–Gerasimenko (Knapmeyer et al., 2016). The second instrumenton InSight, SP for short period, is a MEMS-based seismometer with a flat response to120 s period and a noise level slightly higher than that of a Trillium compact boreholeinstrument on Earth (Pike et al., 2016). Therefore, the report of the Europa Lan-der Science Definition Team (Hand et al., 2017) defined the baseline instrument asequivalent to the Short Period (SP) instrument of the InSight Mars lander.Figure 1 shows a seismogram from Earth, compared with synthetic seismograms foran identical source-receiver configuration on Europa, Enceladus, Titan and Ganymede.This comparison highlights the diverse waveforms that are created by the icy-oceanenvironment. The measured seismogram on Earth is different from the PREM ref-erence synthetic (Dziewo´nski and Anderson, 1981), since the surface waves traveledalong a path through thick continental crust, while PREM contains an oceanic, thin-crust model. No simple global reference model can explain the complete shape of theseismic waveform. Also, ocean-generated microseismic noise appears in the data, withperiods of 7 seconds, and another noise signal at around 100s period. This may act asa rough guide for the differences one could expect between simulations and real datafrom icy worlds.Seismograms on icy worlds will contain exotic phases, especially the prominent long-period flexural waves in the ice and a strong, non-dispersive longitudinal phase on the4-component (both of which were identified and explained in floating terrestrial icesheets by Press and Ewing (1951)). On Enceladus, the Rayleigh (on Z and R) andLove (on T) surface waves orbit the moon in a relatively brief 800 seconds, and theyshow little energy decrease due to the small size of the object. The seismogram istherefore dominated by its repeated occurrences, which permit determination of thedistance between an ice-quake and the receiver (for an application to Mars, see Panninget al., 2015). On Europa, the ice layer is so thin relative to the size of the moon thatamplitude and shape of the seismogram depends very strongly on the thickness of theice layer. On Ganymede, with a 116km ice layer, the overall seismogram shares manyfeatures with the Earth: Body waves with low amplitudes arrive first, followed by acomplex surface wave train, but very weak exotic phases.The paper starts with a description of the models used. In section 3, the seismicwavefield is described, a concise naming scheme is introduced and the various surfacewave types in the ice layer are discussed. Section 4 applies the methodology to Europa,Titan, Ganymede and Enceladus, and identifies potential seismic measurements on thisrepresentative set of icy ocean worlds.
2. Models and Modeling
The structural models used here were adopted from the companion paper Vance et al.(2017b). As described in detail there, the models include up-to-date properties for ices,saline oceans, rocky interiors and iron cores, and are thermodynamically self-consistentwithin the limits of available equations of state. Radial structures are computed asper Vance et al. (2014), with self-consistent ice and ocean thermodynamics, usingboundary conditions of surface and ice-ocean interface temperature. Thermodynamicsfor rocks have been added as per Cammarano et al. (2006). All modeling tools arefreely available via GitHub (http://github.com/vancesteven/PlanetProfile).Temperature profiles can be tuned to produce different ice shell thicknesses. It isalso possible to consider a range of internal compositions and temperature profilesfor the interior below the ocean. We consider the range of plausible ice thicknessesinferred from isostatic surface features, gravitational moment of inertia constraints,and possible radiogenic and tidal heat inputs.For basic estimates of attenuation, we followed the approach of Cammarano et al.(2006) to obtain temperature and frequency dependent estimates of shear quality fac-tor, Q µ with the expression Q µ ω γ = B a exp (cid:18) γH ( P ) RT (cid:19) (1) H ( P ) = g a T m , (2)in which B a = 0 .
56 is a normalization factor, ω is the seismic frequency, exponent γ = 0 . R is the ideal gas constant. H , the activation enthalpy, scales with the melting temperature T m and with the5nisotropy coefficient g a , and the values of g a chosen for various ices are described inVance et al. (2017b). The bulk quality factor, Q κ , is neglected. Values of Q µ / ( ω γ ) inthe ice range from 10 in the colder upper regions of the ice, to 10 approaching thefreezing point, consistent with measurements in glaciers on Earth (Peters et al., 2012).Table 1: Properties of ocean worlds considered here. a Schubert et al. (2004) b Thomas(2010); Iess et al. (2014) c Jacobson et al. (2006); Iess et al. (2012) d Vanceet al. (2014) e Vance et al. (2017b)Radius Bulk density Moment of Inertia ice thickness ocean depth high-P ice?km (kg m − ) km kmEuropa a ± ±
46 0.346 ± .
005 5–30 ¿100 noGanymede a ± . ± . ± . d ¿100 d yes d Enceladus b,c ± ± e e noTitan c ± .
09 1879.8 ± .
004 0.3438 ± . e ¿80 maybeIt is important to note that the mineral composition and structure of the rockyinterior may be very different from Earth’s. For simplicity, we will use the term”mantle” for the rocky layer between water and a potential iron core, which maycontain parts compositionally similar to terrestrial crust or even contain unknownminerals created by exotic pressure-temperature regimes. The reader can find relevantliterature in the companion paper (Vance et al., 2017b). This article deals with modeling the full wavefield of a layered, i.e. spherically symmet-ric icy world. To this end, the spectral-element solver AxiSEM was used (Nissen-Meyeret al., 2014). It separates the problem of wave propagation in a cylindrically or spheri-cally symmetric object into an analytical solution of the problem in the azimuthal ( ϕ )direction perpendicular to the source-receiver plane, and a numerical spectral-elementdiscretisation within the in-plane r, θ , which reduces the numerical cost to that of a2D method (Nissen-Meyer et al., 2008, 2007) and includes attenuation (van Driel andNissen-Meyer, 2014a) and anisotropy (van Driel and Nissen-Meyer, 2014b). The soft-ware was modified to handle general 1D velocity profiles. The changes are included inthe recent release 1.4. The properties of the SEM meshes and the resulting numericalcost are shown in appendix B Compared to the typical meshes for global wave simula-tion on terrestrial planets, smaller elements are needed due to the low wave velocitiesin the ocean. As these are fluid elements, the maximum allowed time step is largerthan on earth, so that the computational cost is about a factor 2 larger than for anterrestrial planet of similar radius.The reciprocity of the Green’s function permits switching of the locations of sourceand receiver of a seismic wavefield. This has been used in the Python package Instaseis (van Driel et al., 2015), which allows to reconstruct seismograms for arbitrary sourceand receiver locations from a precalculated wavefield database. Instaseis uses thestored displacement wavefield to calculate strain, which allows to simulate arbitrary6uropa Ganymede EnceladusFigure 2: Velocity Profiles of three icy moons (Vance et al., 2017b, electronic supple-ment): Europa (left), with a 20 km thick ice layer on top of a 116 km deepocean. Ganymede with (center) with a 124 km thick ice layer on top of a 116km deep ocean and 600 km of high-pressure ice phases. Enceladus (right)has a 52 km ice layer on top of an extremely shallow ocean and no core.Compared to Earth, the velocity gradients within each layer are very smalldue to the relatively low pressure gradients. The only considerable velocitygradients exist in the ocean.moment tensor and single force point sources with arbitrary source time functions.Theoretical arrival times were calculated with TauP software package (Crotwell et al.,1999) in the implementation of ObsPy 1.0 (Krischer et al., 2015).To estimate the effects of three-dimensional heterogeneities and scattering, we in-troduced lateral heterogeneities into some runs. Since AxiSEM assumes an axially-symmetric model, these heterogeneities have to be symmetric around the source axis.This creates a limited scattering effect, especially since no off-path scattering is sim-ulated, but it serves as a good first-order approximation to estimate the total effectof heterogeneities. For the scattering, we implement a von K´arm´an random mediumwith a correlation length of 5 km and velocity variations of 10%. This means thatthe correlation length is similar to the shortest wave length of P-waves. In modelswith lateral heterogeneities to simulate scattering, these structures become symmetricaround the receiver axis.Input models for all planets discussed in this article are available as an electronicsupplement. The seismic wavefields for the models discussed here are available as In-staseis databases on http://instaseis.ethz.ch/icy ocean worlds/ and can be used freelyfor further analysis. Waveform databases for Earth models can be accessed from theIRIS
Syngine service (Krischer et al., 2017).7 . The seismic wavefield
Body waves
P P-wave in upper ice layer or mantleS S-wave in upper ice layer or mantleF P-wave in F luid layerK P-wave in the (solid) core Boundaries
Surface reflections are not specifically named e Bottom of the ice ( e is) layero O cean floorg Top of the rocky layer ( g ranitic crust), if high pressure iceis present below oceanm Top of the m antle, if extra layers are present between oceanand mantleInterface interactions x /ˆ x Underside reflection on interface xx Topside reflection on interface x Specific phases
Pe/Se P or S-waves bottoming in the uppermost ice layerPeP/SeS P or S-waves being reflected at the ocean-ice boundaryPn/Sn ” n ormal” P- or S-waves traversing all water layers, turningupwards in the mantlePoP/SoS P or S-waves being reflected at the ocean floorPmP/SmS P- or S-waves being reflected at the top of the m antle Prefixes for multiples.
Append by integer N for multiple reverberation. y is a wildcard for any following phase.Pf Ny / Sf Ny P or S-waves being reflected N times at the ocean floorand at the ocean/ice boundary; reverberation in the f luidoceanPe Ny / Se Ny P or S-waves being reflected N times at the ocean/iceboundary and at the surface; reverberation in the ice ( e is) Surface waves . N is integer and indicates wave packets traveling along theminor arcs (odd numbers) or major arc (even numbers) of the great circleC N C rary wave in the ice layerR R ayleigh wave at the surfaceRo R ayleigh-like wave on the o cean-floor: Scholte waveRe R ayleigh-like wave at the ice ( e is)-ocean interfaceLF L ong-period F lexural wave, in the ice layerLL L ong-period L ongitudinal wave, in the ice layerLQ L ong-period Toroidal mode ( Q uerwelle), in the ice layerG Love wave, in the ice layerGo Love wave, in the ocean floor9 ceoceancrust eom iceoceanmantlehigh-p icemantle c r u s t Earth Europa &polar Earth Large Ocean worlds v P ≈ P ≈ P ≈ P ≈ P ≈ P ≈ g crust v P ≈ Figure 3: Comparison of the structure of the uppermost layers in Earth and an oceanworld. The letters “e”, “o”, “g”, “m” mark the interfaces in seismic phasenames. Note that sea ice in the polar regions of Earth qualifies as an oceanworld in this definition.To facilitate seismic analysis, it is necessary to be able to name phases (i.e. distinctsignals in the seismogram) unambiguously. The terrestrial naming scheme adopted byIASPEI (Storchak et al., 2003) needs to be extended in the context of icy moons toaccount for the various ice phases. A completely non-ambiguous scheme, which marksevery interface crossing, was proposed by Knapmeyer (2003), but since it producesrather cumbersome phase names, we propose a simplified version here. Lee et al. (2003)proposed a scheme specifically for local waves, but since it cannot be generalized toglobal wave propagation, we decided not to use it here. In our scheme, all namesfrom the IASPEI scheme keep their meaning, where applicable, but we add a few newletters for ice-specific phases. Because “i” and “I” have already been used for innercore phases, we propose “e” as the letter generally related to ice phases, using theGerman word
Eis (analogous to core phases being called “K” from
Kern ).The main issue at hand is that there is only the Moho as a strong, global near-surface interface on Earth, but there are at least 3 in ocean worlds (a Moho, plus theice-bottom and the ocean floor and possibly an interface between high-pressure iceand the crust). Standard IASPEI nomenclature would denote these interfaces by theirdepths. Since this depth is poorly constrained on ocean worlds, we propose to givethese interfaces descriptive letters: “e” for the bottom of a surface ice layer, “o” forthe bottom of an ocean, “g” for the top of the rocky layers and “m” for the top of themantle. To use this scheme for glaciers on Earth or worlds like Ceres, the bottom ofthe ice layer should also be called “e” if it is in contact with a solid rocky interior, see3 for an overview.The IASPEI standard for marking reflections on these interfaces (appended “+”or “-”) is not widely used in the community, so we stick to the nomenclature of Akiand Richards (2002), where topside reflections are noted by the interface descriptor (ordepth), as in PeP, and underside reflections are noted by the descriptor as a superscript(e.g. P m P). Where the usage of a typographic feature is problematic, the mnemonicTauP convention of PˆmP should be used instead.Terrestrial seismology calls compressional waves in the ocean either ”H”, for a hy-droacoustic wave from a source in the water, or ”T”, for a hydroacoustic wave froma source in the ground. Both wave types are linked to the SOFAR channel of lowsound speed between 500 and 1500 meter depth in Earth’s ocean (SOund Fixing And10igure 4: Paths of various P-phases, including multiples from the ice layer, the wa-ter layer and surface reflections for a Europa-like icy moon. Note that thethickness of the ice layer is not to scale.11anging; Tolstoy and Ewing, 1950; Ewing et al., 1951). This channel is an effect ofthe specific contributions of depth-dependent salinity and pressure to sound speed inearth’s ocean and most likely does not exist in ocean worlds. We therefore proposeto mark a leg of the path that traverses the ocean generally with “F” ( f luid). Theterrestrial “H”-wave would then be “FP” and the “T”-wave “pFP”.Waves in the uppermost ice layer and in the mantle are both called P or S, dependingon their type. A PFPFP wave is therefore one that started within the ice layer, crossedthe ocean, traveled through crust and potentially mantle, turned there and went backthrough the ocean to the surface. To retain a certain elegance in naming phases, wepropose to use “Pn” for this wave type, analogous to the mantle phases on Earth asproposed by (Mohoroviˇci´c, 1909). To specifically highlight local phases that never leftthe ice, we propose to add an “e”: “Pe”, “Se”. Note the distinction from “PeP”, whichwas reflected at the ice bottom.To keep phase names concise, we introduce a few more abbreviations: “Po” is shortfor a topside reflection from the ocean floor (full: PFoFP), “Pm” is short for a topsidereflection from the mantle top (full: PFPmPFP). “Pf” is a single reverberation withinthe water (PFoF e FoFP). These abbreviations can also be used as prefixes to markmultiples on other phases, as in “PfPn” = “PFoF e FoFPFP”. See fig. 4 for an overviewof ice- and ocean-interacting body wave paths. As on Earth, a number followingone descriptor means a multiple reflection. “Pf2Pn” was reflected twice at the oceanbottom and twice at the bottom of the ice layer before traveling as a P-wave throughthe mantle.Up-going seismic waves from an event in or below the water layer should be distin-guished. This differs from the terrestrial usage, where source depth is explicitly notcoded in the phase name. We therefore propose to write the first letter in lowercasefor events below the ice, e.g. pFP for an event originating below the ocean layer or fPfor an acoustic source in the ocean. This is a slight extension of the terrestrial usage,where only rays with an upward leg are noted as such. This naming scheme combinesthe brevity of the IASPEI scheme with the flexibility of Knapmeyer’s scheme, while atthe same time not being limited to local waves, as in the scheme by Lee et al. (2003).See table 2 for examples of phase names in our scheme, Knapmeyer’s scheme and theinput to the popular TauP software.We show a visualization of the global wavefield in fig. 5 that follows the example ofthe IRIS global stack (Astiz et al., 1996). Seismograms for different distances are colorcoded, where vertical motion is color coded blue, motion in R direction (away from theevent) in green and transverse motion in T. Each trace is normalized separately. Wesee that in the coda of Pn-waves (see also fig. 7), every body wave phase is followedby multiples from the mantle-ocean and ocean-ice boundaries.
Seismic events within the ice layer will create a range of seismic wave types that arenot observed in global seismograms on Earth (see fig. 6 for a spectrogram, which12igure 5: Global stack of seismograms (Astiz et al., 1996) with phase names accordingto the scheme in table 3.1 for a Europa model with 5 km ice thickness (top)and 20 km ice thickness (bottom). Vertical motion is blue, transverse redand radial green. 13 l e x u r a l w a v e F l e x u r a l w a v e L o ng i t ud i n a l w a v eTo r i o d a l m o d e L o v e w a v e R a y l e i gh w a v e B o d y w a v e s B o d y w a v e s L o ng i t ud i n a l w a v e ( m a j o r a r c ) R a y l e i gh w a v e C r a r y w a v e Figure 6: Spectrogram of an event in a Europa model with 20 km ice thickness in 60degree distance, showing the different spectral characteristics of the wavetypes. Even though the flexural wave and the longitudinal wave dominate abroadband seismogram (see fig. 1), other phases can be identified, becausethey are occur at higher frequencies.14ur scheme(short) our scheme (full) Knapmeyer TauP LeePn PFPFP POPoPopOp P20P136P136P20P N/APe P P P PPoP PFoFP POPo+POP P20Pv136P PCPPnPn PFPFPPFPFP POPoPopOpPOPoPopOp P136PP136P N/APFSFP PFSFP POPoSopOp P20P136S136P20P N/APf PFoF e FoFP POPo+pO-Po+pOp P20Pv136pˆ20Pv136p20p N/AC C N/A 2950kms N/ATable 2: Examples of phase names in our proposed naming scheme, the scheme pro-posed by Knapmeyer (2003) and the one proposed by Lee et al. (2003). Thelatter is only applicable to local phases that did not enter the rocky part belowthe ocean. Surface-wave-like phases like the Crary phase cannot be named inany of the other schemes, only TauP can describe it by its phase velocity.separates the spectral characteristic better; see Table 3 for a summary of the differentwave types, including their group and phase velocities.). This is especially the case forice layers thinner than 40 km, as expected for Europa. The three most prominent newphases are flexural waves, longitudinal waves, and the Crary phase. As shown by Pressand Ewing (1951), the response of a floating ice layer of thickness d is similar to thatof an infinite slab. Kovach and Chyba (2001) applied this to the situation on Europa,so that we can summarize the most important results here, assuming a compressionalwave velocity v P , ice = 4 km / s and a shear velocity v S , ice = 2 km / s: Rayleigh waves:
Short period surface waves with vertical/radial polarization do notsee the bottom of the ice layer and can therefore be described as a Rayleigh wavewith the usual properties, if the inverse wavenumber k − is smaller than the icethickness d : ( kd = 2 π dλ > c R = 0 . v S , ice ≈ . / s is valid. Rayleigh wave (ice bottom):
At short periods, an interface wave travels at the ice-water interface, called
Stoneley wave by Press and Ewing (1951) and Kovachand Chyba (2001) (even though that term usually describes solid-solid interfacewaves). An interesting property is that its maximum phase velocity is controlledby the sound velocity of the ocean : c St = 0 . v P , water ≈ . / s. Its signalis only detectable for ice shells of a few kilometers, but there, the arrival timedifference from the surface Rayleigh wave constrains the sound velocity of theocean. Flexural waves:
Long-period vertically polarized surface waves, with kd < c Fl ∝ √ f that is independent15f the velocity gradient in the ice. Due to this frequency dependence, the shortperiods arrive first, and the group velocity of long-periods goes to zero. Thecutoff frequency is where the Rayleigh regime (see above) starts. Panning et al.(2006) showed the effect of different ice thicknesses on this wave train. Longitudinal wave:
Another long-period solution is a wave with motion in the radialdirection, in phase over the whole ice layer, described analytically in Press andEwing (1951). This wave has a phase and group velocity defined by the expression v L = 2 v S (cid:115) − (cid:18) v S v P (cid:19) , (3)with v L ≈ .
46 km/s for the example P and S velocities. It therefore travels asa fast, non-dispersive pulse along the surface. This means that it arrives muchbefore all other ice phases at long distances, and even before mantle body wavesat intermediate distances (see tab. 4). It is the most prominent feature in thebroadband seismogram. Both flexural waves and longitudinal waves are solutionsfor wavelengths long relative to the ice shell thickness, and so only show up atrelatively low frequencies ( kd <
Love waves:
No shear waves exist in the fluid ocean, so the ice layer is a wave guidefor SH waves, which are completely reflected at surface and ice-ocean interfaces.This results in a dispersive group velocity of U ( f ) = v S , ice (cid:115) − (cid:18) nv S , ice f d (cid:19) , (4)where n = 1 , , , ... is the mode of the Love wave, v S , the shear wave speed and d the ice thickness. Similar to the flexural wave, we have an inverse dispersion,where infinite frequencies travel with U = v S , ice , while lower frequencies areslower. The minimum frequency of the n -th mode Love wave is f min = nv S , ice d ≈ nd mHz , (5)for ice with v S , ice = 2 km / s. This means that f min = 0 . Toroidal wave:
For long periods below the minimum frequency of the guided Lovewave, the whole ice layer moves in phase. Analogous to the longitudinal wave,this results in a non-dispersive, low-frequency wave with high amplitude and aphase and group velocity equal to the SH velocity of the ice shell.
Crary phase:
As a result of SV waves being completely reflected at critical angle forSV-P conversion on the boundaries of the ice, another phase arises with inter-esting properties (Kovach and Chyba, 2001, p. 282): As first described by Crary16ave type frequencies polarization symmetry group vel. U phase vel. c Love wave high transversal antisymmetric v S (cid:114) − (cid:16) nv S fd (cid:17) v S / (cid:114) − (cid:16) nv S fd (cid:17) Toroidal mode very low transversal symmetric ≈ v S ≈ v S Crary wave monochromatic radial n/a ≈ v P (cid:112) − v /v v P Longitudinal low radial/vertical symmetric ≈ v S (cid:112) − v /v v S (cid:112) − v /v Flexural wave low vertical antisymmetric c − λ d c d k v S (cid:113) ρ ice ( kd ) ρ water − v /v kdρ ice /ρ water Rayleigh wave high radial/vertical n/a 0 . v S . v S Table 3: Wave types encountered in ice floating on water. v P , v S are the wave velocitiesin ice, which are assumed to be constant over the depth. ρ ice , ρ water are thedensities of ice and water The definition of symmetry is analog to Lamb waves,i.e. symmetry with respect the the center axis of the ice layer.(1954), with application to terrestrial floating ice in the Arctic, it has mainlyradial displacement, a phase velocity of v P (for homogeneous ice layers) and aharmonic frequency spectrum characterized by f Cr = ( n + 1) v S d (cid:114) − (cid:16) v S v P (cid:17) , (6)where n = 1 defines the characteristic frequency, and larger integer values of n representing overtones. Its group velocity is U ≈ . v P , which means itarrives slightly after the longitudinal wave. Measurement of f Cr constrains theice thickness d , if v S and v P are known. Since f Cr ( d = 5km) = 0 .
46 Hz and f Cr ( d = 20km) = 0 .
11 Hz, the thickness of Europa’s ice layer can be measuredwith a seismometer that is sensitive at these frequencies. For ice thicknessesabove 40 km, the maximum frequency becomes too low to be observable forrealistic events with a space-ready seismometer.Especially for thin ice layers, the wavefield has two frequency regimes: Most ofthe phases specifically related to the ice layer (Crary phase, longitudinal and flexuralwaves) have a maximum frequency between 0.1 and 0.5 Hz (for thicknesses of 20 and5 km resp.), depending on ice thickness. Body waves from the mantle are not limited inmaximum frequency and should thus be clearly distinguishable (see the spectrogram infig. 6). Scattering within the ice layer enhances this effect by damping high frequencyice body waves, with long paths in the ice layer, generally ice-body waves fade awayafter a short distance.
Due to the large velocity contrast between compressional wave velocity in the ice andin the water, no phases like PFP exist in shallow ocean worlds like Europa, while17hey would be a dominant signal in the late seismogram on Ganymede, Titan. OnEnceladus, they will be present, if its ice is thin enough. Attenuation in water is verysmall, so phases that are multiply reflected at the top and bottom of the water layer(Pf, Pf2...) are observable almost over the whole distance range. Their polarization ismainly on the radial component.
Body waves in the rocky interior cover the whole frequency range. On a thin-iceworld like Europa, they dominate the seismic signal above 0.1 Hz. Because of this,a broadband seismometer, sensitive enough to measure in the whole frequency range,can estimate the wave type from its frequency content alone. It also means that, whiledisplacement amplitudes are dominated by ice phases after their arrival, body wavesmight still be recovered from the high-frequency content of the seismogram.The Scholte wave is one particular crustal phase occurring at the ocean bottom.This wave shares most characteristics of a Rayleigh wave, but is defined on a solid-fluid interface (Scholte, 1947). Its phase velocity is frequency-dependent, but belowthe shear velocity in the solid and the sound velocity in the fluid, whatever is lower(Rauch, 1980), i.e. below 1.8 km/s. It has been used in marine exploration on Earthto estimate the shear wave velocity of the uppermost layers below the sea floor (Noletand Dorman, 1996; Kugler et al., 2005). If it can be detected, a detailed dispersionanalysis might be used to constrain the existence of hydrated low-velocity layers onthe sea floor. As it can be seen in Figure 5, it shows strong water multiples.
In terrestrial seismology, several distance regimes are defined, because seismogramschange their first order shape considerably from one to the next (Kennett, 2002):
Local (0 to 100 km)
For shallow events, the ground displacement is controlled bythe near field of the source and significant permanent deformation occurs atthe surface. The local regime can be considerably larger for mega-thrust events(Grapenthin and Freymueller, 2011).
Regional (100–1000 km)
Crustal body wave phases arrive first. Crustal S-waves andsurface waves form one wave train (Lg).
Far regional (1000–3000 km)
Mantle body wave phases arrive first and are separatedfrom the surface wave train, but due to the upper mantle discontinuities, multiple(triplicated) phases arrive within a short time window.
Teleseismic (3000–9000 km)
Mantle body wave phases are clearly separated fromeach other and are visible in the seismogram as well-defined sharp pulses. Directbody waves do not sense the core.
Core shadow (beyond 9000 km)
Direct P-waves are shadowed by the core. Multiplecomplex core phases arrive first, followed by surface-reflections (PP).18oon regional ocean shadow teleseismicEuropa 5 km ice 0 − ◦ (0.2%) 5 − ◦ (6.5%) 30 − ◦ (86.6%)Titan 33 km ice 0 − ◦ (3%) 20 − ◦ (18.3%) 55 − ◦ (75.7%)114 km ice 0 − ◦ (17.9%) - (0%) 50 − ◦ (75.4%)Ganymede 104 km 0 − ◦ (25%) - (0%) 60 − ◦ (50%)Enceladus 15 km 5 − ◦ (41.1%) 80 − ◦ (58.7%)Earth 1 − ◦ (0.75%) 10 − ◦ (5.9%) 30 − ◦ (56.4%)Table 4: Distance ranges of seismograms regimes on different icy worlds. Epicentraldistances between earthquake and receiver are noted in degrees. The numberin parentheses notes the fraction of the total planetary surface that is coveredby this distance regime. If we assume a random distribution of events onthe surface, this is also the probability that any given event will fall intothis distance range from a lander. For comparison, the distance ranges onEarth are also noted, where ”ocean shadow” corresponds to ”far regionalearthquakes”.To summarize, the distance regimes define whether the first arriving phases havemainly travelled through the crust, the complex upper mantle, the relatively homo-geneous lower mantle or the core. The situation on icy ocean worlds is similar. Wepropose to separate the following regimes: Local
Near field effects are dominant, including permanent displacement.
Regional
Body waves that traveled through the ice layer only arrive first.
Ocean shadow
Direct body waves are masked by the curvature of the planet’s surface.The first- arriving phases are longitudinal and Crary waves and Love waves.
Teleseismic range
Mantle body waves arrive first
Core shadow
The planetary body’s core masks direct body waves.The distance range of these regimes depends mostly on the thickness of the ice layer.An exception is the local regime, whose size depends mainly on the magnitude ofthe event, but also on the local velocity structure. On thin-ice worlds like Europa,considerable coseismic displacement will occur over hundreds of kilometers, but itslong periods will not be detectable with a space-ready seismometer. Therefore, we donot discuss this region in detail here.The distance regimes for a set of icy ocean worlds can be found in Table 4. Forthe thin-ice worlds like Europa, it is unlikely that an event will randomly fall into theregional area around a lander, where direct ice phases are well-recordable (unless aspecifically active region has been chosen for the lander). Roughly 5% of randomlydistributed events will be in the ocean shadow range, where special ice layer phaseswill be the first arrivals. The remaining events will be beyond that in the range, wheremantle and core body waves arrive first. 19 .3. Single-station seismology on icy moons
The first efforts at seismology on icy moons will probably be conducted from a sin-gle station. Single-station seismology was common on Earth for decades, and manyfundamental discoveries, like the existence of the inner core, were derived from obser-vations of earthquakes on single observatories (Lehmann, 1936). Other fundamentalwork studied one earthquake at multiple receivers, which is conceptually similar (e.g.the discoveries of the core and the Mohorovicic discontinuity by Oldham, 1906; Mo-horoviˇci´c, 1909). It has long been known that the parameters of earthquakes can bedetermined from one station alone (Ekstr¨om et al., 1986; Wu et al., 2006). A largebody of work on modern single-station seismology has been done in the context of theMars seismometer on the InSight lander (Panning et al., 2015; Khan et al., 2016; B¨oseet al., 2016). The location of a marsquake is determined by the polarization of P andsurface waves and the time difference between specific phases.Because the longitudinal and the Crary wave are prominent features in the seismo-gram of Europa, they can be used for fast determination of the backazimuth of theevent (the direction as seen from the seismometer; see appendix A for a worked exam-ple). Their energy loss due to attenuation is relatively low, so they propagate aroundthe planet multiple times. From the time difference between these orbits, the distancecan be determined. On thick-ice worlds, the Rayleigh wave can take this role. Thisconstraint is best achieved using distant events, using the following work flow:1. Estimate the ice thickness from converted waves in the coda of body waves ofdistant events. This analysis needs no event localization and not even a clearidentification of phases. Figure 7 shows the coda of Pn waves in different Europamodels. Its coda contains clearly recognizable multiples, whose time difference isthe ice thickness divided by the wave velocity. The effect is clearest for very largedistances (135 degree), where the incident angle is almost vertical. The horizontalcomponent seismograms have their maximum energy several seconds after thevertical seismogram. The two lower models in fig. 7 have been modified tosimulate the effect of a highly attenuative ice layer (lowQ) or of a low-attenuation,high-scattering model, similar to the lunar crust (scat). Here, the reverberationsbecome much less pronounced, but are still visible, especially on the long-distancehorizontal components. This method will work best for large distances, where Pnis a clear arrival. The 45 degree seismograms contain triplicated phases, whichcreates the very complex waveform.2. With the constraint from step 1, refine the thickness estimate from the character-istic harmonic frequency of the Crary wave (for ice layers ¡ 50 km) and estimateP- and S-wave velocities in the ice. Figure 8 shows spectra of the Crary wavefrom simulated waveforms for two Europa models, Titan and Ganymede. A timewindow after the high amplitude longitudinal wave LL was selected, where theCrary wave can easily be found as a harmonic signal, even without identifyingany other phases. For all realistic ice thicknesses of Europa, the Crary peaks arevery prominent in the spectra of the R-component and can be used to determinethe ice thickness reliably. For ice layers of more than 100 km, the peaks are less20rominent and the implicit assumption of a thin ice layer with constant wavespeeds is not valid anymore, as seen for the Ganymede model.3. With an ice layer model, locate local events, using the time difference betweenPe and Se phases, as well as reflections from the ice-ocean boundary (PeS orPeP), once the ice thickness has been constrained. The distance of a local eventwith P and S arrival times t P , t S is x = t S − t P /v S − /v P ≈ ( t S − t P ) · / s . (7)4. From multiple-orbit Crary, longitudinal or Rayleigh waves, locate distant events.Appendix A shows an example of determining the direction of an Europaquakeusing the polarization of body waves and the energy content of the horizontalchannels. If minor arc ( t ), major arc ( t ), and single orbit ( t ) arrival times ofany surface wave are known, the distance of the event is (Panning et al., 2015)∆ = 180 ◦ (cid:18) − t − t t − t (cid:19) . (8)Note that the group velocity of the wave is not needed. If it is known, only twoof the arrival times are needed.5. Identify ocean multiples in the coda of body waves of distant events and constrainthe ocean depth by estimating the two-way travel time of reflected phases.6. Identify Scholte waves at the ocean bottom to constrain the uppermost crustalstructure below the sea floor from spectral peaks or dispersion curves (see Sect.4.2 for an example on Titan).7. Identify as many body wave arrivals as possible to constrain overall mantle ve-locity structure with the event locations from step 4, analogous to the methodproposed by Khan et al. (2016) for Mars.8. If direct ocean phases (PFP, SFS) have been identified together with other bodywaves, the speed of sound in the ocean can be estimated, constraining the oceanchemistry.The process sketched here will be an iterative one in practice. Steps 1 and 2 can be doneeven in noisy seismograms without clearly identified phases or located events. Steps3 and 4 need some identified phases and benefit from first estimates on ice thicknessand ice velocities. The steps 5-8 need first estimates of event distances to work, buttheir results can be used to better constrain the velocity models in turn. While theprocess is described here only for a single event, observations of multiple events willbe necessary to obtain reliable results. Recent estimates of the seismicity induced bytidal cracking on Europa (Panning et al., 2017) suggest at least one globally observableevent per week. The seismicity on other ocean worlds is not yet estimated.21igure 7: Coda of a Pn wave measured at 45, 90 and 135 degree distances for var-ious Europa models. Black lines correspond to vertical displacement (Z-component), red to horizontal pointing away from the source (R-component).Clearly visible are the reverberations from the ice layer, which will serve toestimate the ice thickness from a single seismogram. P o w e r / d B Europa, 5kmEuropa, 5km Titan, 46kmTitan, 46km period / seconds
Europa, 20kmEuropa, 20km period / seconds
Ganymede, 126kmGanymede, 126km
Crary, RCrary, ZCrary peak & overtones
Figure 8: Spectra of a time window around the arrival of the Crary phase for two Eu-ropa models of ice thicknesses of 5 and 20 km and Titan and Ganymede. Theperiod of the main mode and overtones of the Crary phase are proportionalto the ice thickness (eq. 6). 22 . Icy ocean worlds in the solar system
Europa has been the subject of the most detailed studies to date, and a lander missioncould take place in the coming decades (Hand et al., 2017). Therefore, we discuss wavepropagation and measurability of seismic waves in detail here.
Compared to its radius, Europa has a thinner ice layer than all other icy ocean worlds.This creates pronounced ice-layer-specific seismic phases. The transition from a ice-Rayleigh wave to flexural waves occurs between 0.1 Hz (for 20 km ice thickness) and0.5 Hz (for 5 km ice) (Panning et al., 2006), which would be well observable with awide-band seismometer, as the SP instrument used for InSight (Pike et al., 2016) andproposed as the baseline instrument for the NASA Europa lander (Hand et al., 2017).Apart from the thickness of the ice layer, which can be measured as describedabove, or from the Crary phase, key science questions on Europa are the depth andcomposition of the ocean, the presence of liquid or mushy pockets in the ice andthe potential of material transfer between the rocky interior, the ocean and the icelayer. The depth of the ocean can be determined well from the time delay betweenocean multiples (PfPn, Pf2Pn, ...) in the P-wave coda of teleseismic Europaquakes.The ocean chemistry can be estimated from the sound speed in the ocean, whichis increased by higher salt contents. This increases the curvature of the SFS path,resulting in higher amplitudes for SFS at short periods; a phase that would not beobservable in a pure water ocean. This could be used to constrain the ocean chemistry,once the ocean depth and thickness have been determined. Liquid pockets close tothe lander would result in scattering of seismic energy from certain incidence angles.Motion of liquids in the ice shell would create seismic signals similar to geysers on Earth(Kedar et al., 1998). The sharpness of the ocean-rock interface, including potentialhydrated layers could be estimated from Scholte waves on the seafloor.Due to Europa’s modest size, multiple planetary orbits of seismic ice-phases, espe-cially the longitudinal waves, would be observable with an SP-like instrument, allowingfor a straightforward inversion for the distance of a seismic event. Due to the thin icelayer, no direct SH waves are possible beyond 100 km, so that the transverse compo-nent will be very distinct from the radial component, allowing an estimation of theevent backazimuth.In a companion paper (Panning et al., 2017), the background seismicity from tidalcracking is estimated, showing that it is detectable with a seismometer. This seismicitywill create a significant seismic hum, which exceeds noise floor of industry standardTrillium compact instruments (Ringler and Hutt, 2010) at periods between 1 and 10seconds. It could be be used for ambient noise analysis techniques, to constrain thethickness of the ice layer or mantle discontinuities.Another important seismic source will be surface tectonics. The young surface ofEuropa shows clear signals of subduction Kattenhorn and Prockter (2014), collapses23f cavities (Walker and Schmidt, 2015) and potentially plume activity, which would allhave specific seismic signatures (Vance et al., 2017a). While no atmosphere exists tocreate microseismic noise, the vicious circulation in the ocean (Soderlund et al., 2013)might just be strong enough to create a detectable long-period signal (Panning et al.,2017). d i s p l a c e m en t [ n m ] Body waves Crary phase Flexural and Love waves
Displacement seismogram of a M W P o w e r [ m / s / H z ], [ d B ] Spectrum, Body waves
Spectrum, Crary phase P o w e r [ m / s / H z ], [ d B ] Spectrum, Flexural and Love waves
ZRT90% c.l. ofcracking noise
Figure 9: Measurability of seismic phases on Europa: The upper plot shows an ex-ample seismogram of a M W P [km/s] v S [km/s] ρ [kg/m ]Iceh 46 3.9 2 925ocean 410 1.55 – 2.55 – 1020 – 1170IceVI 5.5 4.4 2.26 1394mantle 7.99 4.54 3526for a magnitude 3.3 event at a distance of 90 degree (fig. 9). Body waves exceed theinstrument self-noise by 20 dB between 0.1 and 2 Hz, and will therefore be observable.The longitudinal phase, which is crucial for event location, exceeds the self-noise by30 dB. The late seismogram on R and Z, including Rayleigh and flexural waves, canbe measured for periods up to 50s. The Love wave on T, which is confined to theice layer and therefore strongly affected by attenuation and scattering in the ice, isnot observable for most focal mechanisms. Based on the seismicity models in Panninget al. (2017), a magnitude 3.3 event can be expected to occur somewhere between onceper week and once per month. Compared to Mars or Earth, an instrument on Europawill not be affected by atmospheric noise, so instrument self-noise is the reference tomeasure signal strength against. Titan is predicted to have water and ice layers with total thickness around 480 kmabove a low-density silicate mantle. The ice covering is interpreted to be 55-80 kmbased on the observed Schumann resonance (B´eghin et al., 2012). Following Vanceet al. (2017b), we consider thicknesses in a slightly wider range between 46 and 118 km.Depending on the ocean’s temperature profile with depth, high-pressure ice phasesmight be present below the ocean, with a total thickness of up to 230 km or notpresent at all. Depending on the ocean temperature, the ice thickness may be so highthat the longitudinal, Crary, and flexural waves are restricted to very long periods,unobservable with a realistic instrument. Because the rocky mantle makes up morethan 80% of the diameter in all models, the teleseismic range dominates on Titan (seefig. 10).An important question for habitability in Titan is the existence of high-pressureice phases at the bottom of the ocean. If they are present, they may form a barrierto water-rock interactions that regulate redox and ocean composition, and which arethereby critical for supporting life. While intra-ice convection might be able to trans-fer ions (Choblet et al., 2017), the absence of a high-pressure ice layer could benefithabitability by providing a direct water-rock interface, potentially with life-supportinghydrothermal systems analogous to those on Earth’s seafloor. Because such a high-pressure ice layer is not detectable with radar or gravitational moment-of-inertia mea-surements, seismological measurements provide an unique probe. As was argued inVance et al. (2017b), the thickness of surface ice above the ocean and high-pressureice below are strongly correlated and both mainly controlled by the ocean ion content25igure 10: Global stack of seismograms for a Titan model of 46 km ice thickness (purewater ocean and 270K in table 7 of Vance et al. (2017b), top, compared toone with 119 km ice (bottom, model with 3%NH ocean and 255K in table7 of Vance et al. (2017b)). The thinner ice model has a seismic wavefieldresembling that of Europa, with a dominant longitudinal phase that is thefirst arriving phase between 20 and 55 degrees. In the thick-ice model,seismic phases from high-pressure ice are negligible, and the wavefield ismore similar to a terrestrial planet, where body waves always arrive first.In contrast to Earth or Mars, every phase is followed by a large number ofmultiples from ice and/or the ocean.26
50 800 850 900 950 1000 1050 1100time after earthquake in seconds30201001020304050 d i s p l a c e m e n t / n a n o m e t e r Pn-wavetrain Sn-wavetrain body wave coda in Titan, distance: 90 degree, M W :3.86 5.5 km iceVIno iceVI A m p li t u d e / db ( m / s ) Pn-wave spectrum
Sn-wave spectrum
Figure 11: Detectability of an high-pressure ice layer (IceVI) below the ocean frombody and Scholte waves. The upper plot shows the vertical displacementof an M3.8 event in 90 degree distance and 5 km depth for two interiormodels with 46 km ice thickness. The red lines correspond to a model witha 5.5 km IceVI layer the bottom of the ocean. The model displayed withblack lines is identical, but the mantle starts directly below the ocean. Thepeaks in the S-spectrum correspond to reverberations in the ice layer andare clearly identifiable. 27nd temperature profile. According to these simulations, every model with a surfaceice thickness of more than 50 km must also contain high-pressure ice layers. Seismicmethods could constrain surface ice and high-pressure ice thicknesses independently,and thereby constrain ocean chemistry and temperature.As one example, Figure 11 shows a comparison of waveforms for two thin surface icemodels (46 km) on top of a 410 km deep ocean with 3.3 weight percent NH (similarto the 3 wt percent NH model with 264 K in table 8 of (Vance et al., 2017b)). Forthis surface ice thickness, high pressure ice can neither be confirmed nor excluded.The figure compares seismograms of this model with one where the iceVI layer hasbeen replaced by mantle (see table 5). Because there is no distance where high-pressure ice phases arrive first, their presence will have to be inferred from the codaof mantle body waves. The strong velocity and density increase at the ocean-mantleinterface (assuming that no iceVI layer is present) causes a strong impedance contrast,which is weakened if there is an intermediating layer. The lower two plots in fig. 11show acceleration spectra of a window around the P-wave arrival (left gray patch inupper subfigure) and the S-wave arrival (right gray patch in upper subfigure). Theamplitude of the Pn wave train in a model with an ice layer is two times higherthan in the model without the layer. Since it will be difficult to reliably determinethe magnitude and focal mechanism of a seismic event with just one station, thiseffect may go unnoticed. The spectrum of the S-waves (lower right) however, showsclear effects of the iceVI layer: While the spectrum of the simple model is almostflat between 1 and 10 seconds, the iceVI layer produces two clear peaks at 2 and 5.1seconds, corresponding to reverberations of P and SV waves. This signal should beeasily distinguishable, even without a detailed source model.Titan is the only icy ocean world with an atmosphere. Due to the atmosphere’s lowdensity in comparison with the solid and liquid portions, the direct effect on the seismicwave propagation in the planet will be minimal; however, it could serve as a sourcefor seismic signals. Wind waves on Titans methane lakes have been predicted to reachsignificant wave heights of 0.2 m at periods of 4 seconds (Lorenz and Hayes, 2012)for winds of 1m/s. The handful of Cassini measurements available through spring didnot show evidence of waves above a few millimeters in height e.g. (Hofgartner et al.,2014; Zebker et al., 2014) but these may simply have been on calm days. If or when0.2m waves form, they create pressure on the sea floor of the order of 100 Pa, whichis comparable to pressure variations on the abyssal sea floor on Earth due to oceanwave interference (Cox et al., 1984; Davy et al., 2014; St¨ahler et al., 2016). Becausethe latter creates a globally observable seismic signal (Kedar et al., 2008), it must beassumed that measurable ocean-generated microseisms are possible on Titan. Surfacewaves from these signals could be used to constrain the uppermost layers below thesurface. Note, that while Titan’s seas occupy 1% of Titan’s surface area total, theycover some 12% of the terrain northward of 55 degree N (Hayes, 2016), so microseismsmight be expected to be a much more prominent feature of the seismic environment inTitan’s high northern latitudes than elsewhere, and more so during the summer whenwinds are expected to be strongest.Pressure changes in the atmosphere might also create seismic signals directly, asobserved on long periods on Earth (Peterson, 1993; Beauduin et al., 1996). Local28igure 12: Global stack of seismograms for a Ganymede model of 104 km ice thickness.The thick ice layer creates nicely separated body wave reverberations in theice. Seismic waves that bottomed in the high pressure ice layer arrive in thecoda between 20 and 90 degree distance. They seem to connect to Sn at 90degree, since the P-velocity in the high-pressure ice is very close to the S-velocity in the mantle. The overall wavefield is very similar to a cold-oceanTitan.phenomena, such as ”dust devils” on Earth and potentially Mars (Lorenz et al., 2015)might also be observable. Seismic waves may also create atmospheric signals that mightbe detectable from orbit, as observed on Earth and proposed for Venus (Lognonn´eet al., 2016). The density structure inferred from moment of inertia measurements suggests thatGanymede has water layers in different forms, with total thickness of about 800 kmabove a rocky mantle and an iron core (Vance et al., 2014). The uppermost ice layeris assumed to have a thickness of more than 50 km, floating on a liquid ocean of morethan 100 km depth. Due to the high gravity and abundance of water, pressure atthe ocean floor creates layers of ice V and VI 10s to 100s of km thick, between theliquid ocean and the rocky mantle. Similar to Titan, the existence and thickness ofthese layers can only be measured by seismology. Since both are controlled by ocean29hemistry and temperature, this creates a direct geophysical observable for habitabilityconditions of the ocean (Vance et al., 2017b).The thick ice layer means that longitudinal and Crary phases are relatively weakand constrained to long periods, which makes them unobservable without a true high-sensitivity broadband seismometer. The ice thickness therefore would have to be con-strained from multiples in the coda of body waves. Compared to Titan, the thicknessof the combined water layers is 30% greater, which clears the coda of body waves.As fig. 12 shows, the wavefield is nevertheless dominated by water multiples. Theanalysis presented in the previous section to detect high-pressure ice phases could bedone in a similar fashion for Ganymede, though the high-pressure ice layers are muchthicker and more heterogeneous than on Titan which will shift the resonances to longerperiods.Ganymede’s interior is likely to be at least as seismically active as that of Earth’smoon, where continuous seismicity was observed by Apollo instruments (Vance et al.,2017b). A major question at Ganymede involves the presence of a liquid iron core,which seems to be required by Galileo observations of an intrinsic dipolar magnetic field(Kivelson et al., 2002). Melt could be present in the overlying rock as well. This studynecessarily focuses on the initial problem of identifying radial boundaries closer to thesurface, but a study focused on the goal of long-term exploration of Ganymede wouldneed to evaluate the ability to probe the deeper interior to understand Ganymede’scomposition and current thermal state.
Enceladus’ global subsurface ocean has been confirmed from its measured shape (McK-innon, 2015) and libration (Thomas et al., 2016), but the ocean thickness probablyvaries globally, exceeding 15 km at the south pole due to a thinner ice shell, andless than 10 km everywhere else. For the simulation, we assumed a constant globalthickness of 10 km.Due to its small size, geometrical spreading and attenuation of seismic energy playsa small role on Enceladus. As fig. 1 and the global wavefield stack in fig. 13 shows,even for magnitude 3 events, more than 5 orbits of longitudinal phases are observable,which would allow for easy determination of ice thickness and event locations. Tidalheating from the resonance with Dione suggests a high energy budget for seismicity,consistent with the active tectonic deformation of the South Pole region (Porco et al.,2006).The plumes at Enceladus’ south pole (Porco et al., 2006; Hansen et al., 2006) arespecifically caused by geysers or cryovolcanism. Since geysers are sources of seismicsignals on Earth (e.g. Kedar et al., 1998), the south pole region of Enceladus canbe expected to have a very high seismic background noise. At the same time, thesubterranean volcanic system means that the seismic velocity structure will be veryheterogeneous, which may mean that the concepts of global seismology derived fromspherically symmetric models presented here are not applicable.In addition to these sources, Enceladus should also regularly generate seismic signalsin the ice as a result of tidal flexing (Vance et al., 2017a). The combination of these30 n P w P n P e P e d i ff P e P P e P e P P e P P e P P e P P e d i ff T o r o i da l = L Q R a y l e i gh = R S e S S e P S e S P e P S e S LQ1LQ2 LQ3LQ4LQ6 LQ5LQ8 LQ7LQ9R1 R2R3 R4R5 R6R7 R8 S e 4 S e 5 P n o P n C r a r y Figure 13: Global stack of seismograms for a Enceladus model of 52 km ice thickness.Since this is 20% of the radius, ice phases (Pe, Pediff are the first to arriveup to 75 degree distance) and surface reflected ice phases are prominentin the coda for all distances. The ice layer is too thick for a measurablelongitudinal or Crary phase. Due to the small diameter and the relativelylow attenuation, surface waves orbit the moon multiple times and dominatethe seismogram (see fig. 1). 31ources merits further investigation building on the study of Panning et al. (2017).
5. Discussion
This study gives an overview of general characteristics of icy moon seismology on aglobal scale, but is in no way exhaustive. An obvious limitation is that only spher-ical models were studied so far. The above-mentioned methods to determine the icethickness are likely sensitive to ice undulations (e.g., Nimmo et al., 2007; Lef`evre et al.,2014). The characteristic frequency of the Crary wave is sensitive to the average thick-ness. The Love wave cutoff is sensitive to the minimum thickness, while body wavecoda measures the thickness at the location of the source and the receiver. Enceladusfor example has an ellipticity of 1/50, and probably a variation in ice and ocean thick-ness of more than 50% between the south pole and the equator (Thomas et al., 2016),which will strongly affect propagation, especially of multiply orbiting surface waves.From simulations using cylindrically symmetric models, we found that the longitudi-nal wave is a very robust feature and multiple orbits are detectable, even with strongheterogeneities in the ice. Undulations of the ice bottom may, however, have a verystrong effect on the amplitude and characteristic frequency of the Crary wave andtherefore on its effectiveness in constraining ice thickness.On Europa and Enceladus, the existence of plumes suggests that liquid channelsexist in the ice that will strongly scatter seismic waves. Therefore, full-3D simulationswill be necessary as a next step.Seismic wave propagation on icy ocean worlds is a numerically challenging problem,since it requires high frequencies of up to one Hz to resolve seismic features like theCrary wave for thin ice layers. Since the numerical cost in full 3D methods typicallyscales with frequency f as O ( f ), this prohibits their application on planet scales. Thisis problematic, given that the ice waves may be very sensitive to three-dimensionalstructure, like varying ice thickness. This is rather different from the situation onEarth, were surface waves are strongest at longer periods (typically above 15 seconds),which are possible to simulate on large HPC systems. Seismic wavefield solvers suitablefor smooth three-dimensional structure at planetary scales at frequencies of up to oneHz are just becoming available (Leng et al., 2016) and will offer more insight intothe question which phases are usable. Local structure around a lander may stronglyaffect the measured seismic waveforms, which may be simulated by including local 3Dwavefield simulations into global 1D simulations.This article is intended to set the stage of global seismic wave propagation in icymoons on which more detailed studies of effects of three-dimensional structure canbuild.Since the presented analysis of the seismic wavefields is in alllikelihood incomplete, we make the Instaseis databases available at http://instaseis.ethz.ch/icy ocean worlds/ . The scripts to reproduce the figuresare available on the seismo-live website .Future research should also explore the parameter space of possible mantle compo-sitions and thermal states. As was shown here, seismology in principle provides tools32o constrain both, but detailed studies are necessary for each planetary body, in closecollaboration with planetary geophysics (see table 6 for a brief summary).
6. Conclusion
In this paper, we discussed the general characteristics of global seismology on icy oceanworlds. The existence of the liquid-filled gap between the mantle and the icy crustcreates a different wavefield than the one that is known on earth. In general, thismeans that a large part of seismic energy remains close to the surface, where it can bemeasured, which makes even relatively small events observable globally. Also, severalseismic measurables are specifically sensitive to parameters related to habitability.This is a promising perspective for the short term installations of seismometers on icymoons.
7. Acknowledgements
The authors acknowledge computational support in the project pr63qo ”3D wave prop-agation and rupture: forward and inverse problem” at
Leibniz-Rechenzentrum
Garch-ing. SCS was supported by grant SI1538/4-1 of Deutsche Forschungsgemeinschaft
DFG . MvD was supported by grants from the Swiss National Science Foundation(SNF-ANR project 157133 ”Seismology on Mars”) and the Swiss National Supercom-puting Center (CSCS) under project ID sm682.This work was partially supported by strategic research and technology funds fromthe Jet Propulsion Laboratory, Caltech, and by the Icy Worlds node of NASA’s As-trobiology Institute (13-13NAI7 2-0024). A part of the research was carried out at theJet Propulsion Laboratory, California Institute of Technology, under a contract withthe National Aeronautics and Space Administration.The Instaseis wavefield databases available at http://instaseis.ethz.ch . The scriptsto reproduce the figures are available on the seismo-live website .All rights reserved prior to publication.33cientific objective Seismic observable signalstrength distance range referenceTectonic activity Location of seismic events strong global appendix AIce thickness Resonant frequency ofCrary phase strong global fig. 8Transition frequency be-tween Rayleigh and flexu-ral surface wave intermediateglobal fig. 6Reverberations in bodywave coda strong teleseismic fig. 7Autocorrelation of ambi-ent noise strong ambient noise (Panninget al., 2017)Ocean depth Reverberations in bodywave coda strong teleseismic fig. 7Autocorrelation of noise intermediateambient noise (Panninget al., 2017)Ocean chemistry(from sound veloc-ity in water) Reverberations in bodywave coda strong teleseismic fig. 7High pressure icephases Coda of Sn-waves strong teleseismicfig. 11 Scholte waves strong teleseismicP-to S ratio intermediateteleseismic, de-pends on focalmechanism fig. 11Core diameter Autocorrelation of noise weak ambient noise (Panninget al., 2017)Autocorrelation of seismo-gram intermediate >
100 deg fig. 7Table 6: Potential scientific objectives of a future icy ocean world seismometer andseismic observables that address it.34 . Determine the backazimuth of an event from aEuropan seismogram
As described in section 3.3, the radial and transverse component of ocean world seismo-grams differ strongly, which can be used to estimate the backazimuth of an earthquake(i.e. the direction of the earthquake as seen from the receiver). We tested a simpleautomated implementation of two methods.1. Maximize the energy of the Longitudinal wave. The Longitudinal wave is purelylongitudinally polarized. Therefore, the horizontal components of the seismo-gram can be rotated such that energy in the time window around the longitu-dinal wave is maximal one one horizontal component (R) and minimal on theone perpendicular to it (T). Note that in the presence of noise, the energy in theT-component will not be zero. Also, this method has a 180 degree ambiguity,since the polarity of the longitudinal wave can be positive or negative.2. Maximize coherence between vertical and radial component for body waves. Thefirst arriving body waves are compressional waves with a high incidence angle,which have motion purely along the direction of propagation (which is a lin-ear combination of the Z-direction and the R-direction pointing away from theevent). The direction away from the event is the one, which has maximal corre-lation to the vertical component. This method has a much lower resolution thanthe first one, but has no 180 degree ambiguity and so can used to remove thatambiguity from the higher reolution method 1.Figure 14 shows the application of both methods to a M3 Europaquake in 83 degreegreat circle arc distance with a true backazimuth of 76 degree for a Europa modelwith 5 km ice thickness. It shows that both methods together correctly identify thebackazimuth within a few degree.
B. Properties of the numerical simulations
The numerical simulations presented in this manuscript were done using the opensource spectral-element wavefield solver AxiSEM. Some modifications were done onthe mesher, to make it more stable around very strong velocity contrasts close tothe surface. These modifications were added to the main development branch of thesoftware and are included in the recent version 1.4. The simulations were run onthe Linux HPC systems
SuperMUC of the Leibniz-Rechenzentrum Garching and
PizDaint of the Swiss National Supercomputing Centre CSCS.Due to the low velocity in the ocean, more (small) elements are needed and simula-tion is about 50% more expensive compared to a terrestrial planet of the same radius.A 1 Hz simulation for Ganymede or Titan is therefore slightly more expensive thanone of the Mars runs done for the InSight blindtest (Ceylan et al., 2017; Clinton et al.,2017). 35ody Radius Period nelement time step nsteps costkm s ms 10 CPUhEnceladus 0.25 1406304 2.73 1854.5 39990.5 355224 5.6 228.4 4921 90630 11.4 28.6 62Europa 1565 1 2139072 8.88 867.4 1871(5km ice) 2 546000 15.2 129.1 278Europa 1 2047690 8.95 824.0 1777(20km ice) 2 521136 15.4 121.6 262Titan 2574 1 7663680 12.6 2189.5 4722(33km ice) 2 1929600 25.2 276.1 595Titan 1 7416318 9.94 2687.3 5795(124km ice) 2 1873808 19.9 339.6 732Ganymede 2631 1 8281306 11.5 2584.6 55742 2090130 23 326.8 705Mars 3396 1 5098080 7.7 2371.5 5114Earth(PREM) 6371 1 9520000 7.7 4436.8 9568Table 7: Properties of the AxiSEM simulations for this paper. Since AxiSEM is a 2Dmethod, the numerical cost scales roughly with the inverse cube of the mini-mum period. NCPU is the optimal number of CPUs for mesh decomposition;nelement the number of elements of the SEM mesh. nsteps is the numberof elements times the number of time steps for a 1 hour simulation. Thisis directly proportional to the computational cost, which is shown here onPiz Daint, a Cray XC40 system at the Swiss National Supercomputing Cen-tre CSCS. For comparison, the cost of a 1 hour simulation of Mars (modelDWAK in Clinton et al. (2017)) and Earth (PREM) are shown. These costsare based on using gfortran 6.1 and may be lower for more optimized compilers36 Z originalrotated, body wave 6004002000200400600 Z originalrotated, Longitudinal wave10050050100 d i s p l a c e m e n t / n m N/R d i s p l a c e m e n t / n m N/R
376 378 380 382 384 386 388 time after event / sec
E/T
500 540 580 620 660 time after event / sec
E/T angle / degrees E n e r g y i n L o n g i t u d i n a l t i m e w i n d o w BAZ from Longitudinal energy: 70 degree
True BAZ: 76 degree B o d y w a v e c o h e r e n c e h o r / v e r t BAZ from body wave: 132 degree
Figure 14: Determining the backazimuth with two methods: 1. The energy of theCrary wave should be maximal on the R-component and minimal of T. 2.The first arriving body waves are mainly P-SV waves, therefore coherencebetween Z and R should be maximized by rotating the seismogram intothe correct backazimuth. The seismogram belongs to a M 3.1 event inteleseismic distance with a true backazimuth of 76 degree. The right figureshows the original waveforms in the selected time windows compared to theones in the optimally rotated system.
C. Retrieving seismic waveforms for ocean worlds
While the forward simulations of the seismic wavefield require thousands of CPU-hours; the
Instaseis -database-approach allows to reuse the stored wavefield to calculateseismograms for arbitrary source-receiver locations on the fly (van Driel et al., 2015).The disk space requirement of the databases depends on the maximum depth of anevent, but is generally in the range of several ten Gigabyte. Databases for the modelsshown in table 7 are stored on a server at ETH Z¨urich. Similarly to the terrestrialdatabases stored at the IRIS service
Syngine (Krischer et al., 2017), they can be openlyaccessed via the Python package
Instaseis M rr event in 40 degreedistance on Ganymede, including storage into a MiniSEED file. import instaseisdb_path = ’http://instaseis.ethz.ch/icy_ocean_worlds/Gan126km-00pMS-hQ_hyd30km_2s’db = instaseis.open_db(db_path)src = instaseis.Source(latitude=0.0, longitude=0.0, m_rr=1e17)rec = instaseis.Receiver(latitude=40, longitude=00.0)st = db.get_seismograms(src, rec)st.write(’Ganymede_event.mseed’, format=’MSEED’) eferences Keiiti Aki and Paul G Richards.
Quantitative Seismology , volume II. University ScienceBooks, 2nd editio edition, 2002. ISBN 0935702962.L. Astiz, P. Earle, and Peter M. Shearer. Global Stacking of Broadband Seismo-grams.
Seismological Research Letters , 67(4):8–18, 1996. ISSN 0895-0695. doi:10.1785/gssrl.67.4.8.R. Beauduin, Philippe Lognonn´e, J. P. Montagner, S. Cacho, J. F. Karczewski, andM. Morand. The effects of the atmospheric pressure changes on seismic signals orhow to improve the quality of a station.
Bulletin of the Seismological Society ofAmerica , 86(6):1760–1769, 1996. ISSN 00371106.Christian B´eghin, Christophe Sotin, and Michel Hamelin. Titan’s native ocean re-vealed beneath some 45 km of ice by a Schumann-like resonance.
Comptes Rendus -Geoscience , 342(6):425–433, 2010. ISSN 16310713. doi: 10.1016/j.crte.2010.03.003.URL http://dx.doi.org/10.1016/j.crte.2010.03.003 .Christian B´eghin, Orlien Randriamboarison, Michel Hamelin, Erich Karkoschka,Christophe Sotin, Robert C Whitten, Jean Jacques Berthelier, Rjean Grard, andFernando Simes. Analytic theory of Titan’s Schumann resonance: Constraints onionospheric conductivity and buried water ocean.
Icarus , 218(2):1028–1042, 2012.ISSN 00191035. doi: 10.1016/j.icarus.2012.02.005.Maren B¨ose, John F. Clinton, Savas Ceylan, F. Euchner, M. van Driel, Amir Khan,D. Giardini, Philippe Lognonn´e, and W.B. Banerdt. A Probabilistic Frame-work for Single-Station Location of Seismicity on Earth and Mars.
Physicsof the Earth and Planetary Interiors , 262:48–65, 2016. ISSN 00319201. doi:10.1016/j.pepi.2016.11.003.L. Bruzzone, J.J. Plaut, G. Alberti, D.D. Blankenship, F. Bovolo, B.A.Campbell, A. Ferro, Y. Gim, W. Kofman, G. Komatsu, W. McKin-non, G. Mitri, R. Orosei, G.W. Patterson, D. Plettemeier, and R. Seu.RIME: Radar for Icy Moon Exploration. In , pages 3907–3910. IEEE, 72013. ISBN 978-1-4799-1114-1. doi: 10.1109/IGARSS.2013.6723686. URL http://ieeexplore.ieee.org/document/6723686/ .Fabio Cammarano, V. Lekic, M. Manga, Mark P Panning, and B. Romanowicz.Long-period seismology on Europa: 1. Physically consistent interior models.
Jour-nal of Geophysical Research , 111(E12):E12009, 12 2006. ISSN 0148-0227. doi:10.1029/2006JE002710. URL http://doi.wiley.com/10.1029/2006JE002710 .P. S. Cassen, J. Peale, and R. T. Reynolds. Is there liquid water on Europa?
Geo-physical Research Letters , 6:731–734, 1979.38avas Ceylan, Martin van Driel, Fabian Euchner, Amir Khan, John F. Clinton, LionKrischer, Maren B¨ose, Simon St¨ahler, and Domenico Giardini. From Initial Modelsof Seismicity, Structure and Noise to Synthetic Seismograms for Mars.
Space ScienceReviews , 211(1-4):595–610, 10 2017. ISSN 0038-6308. doi: 10.1007/s11214-017-0380-6. URL http://link.springer.com/10.1007/s11214-017-0380-6 .G Choblet, G Tobie, C Sotin, K Kalousov´a, and O Grasset. Heat transport in thehigh-pressure ice mantle of large icy moons and the dissipation number.
Icarus ,285:252–262, 2017. ISSN 0019-1035. doi: 10.1016/j.icarus.2016.12.002. URL http://dx.doi.org/10.1016/j.icarus.2016.12.002 .John F. Clinton, D Giardini, P Lognonn´e, B Banerdt, M. van Driel, M Dril-leau, N Murdoch, M. Panning, R Garcia, David Mimoun, M Golombek,J Tromp, R Weber, Maren B¨ose, Savas Ceylan, I Daubar, B Kenda, AmirKhan, L Perrin, and A Spiga. Preparing for InSight: An Invitation to Par-ticipate in a Blind Test for Martian Seismicity.
Seismological Research Letters ,88(5):1290–1302, 9 2017. ISSN 0895-0695. doi: 10.1785/0220170094. URL http://srl.geoscienceworld.org/lookup/doi/10.1785/0220170094 .Charles Cox, Thomas Deaton, and Spahr Webb. A Deep-Sea Differential PressureGauge.
Journal of Atmospheric and Oceanic Technology , 1(3):237–246, 9 1984.ISSN 0739-0572. doi: 10.1175/1520-0426(1984)001¡0237:ADSDPG¿2.0.CO;2..A. P. Crary. Seismic studies on Fletcher’s Ice Island, T-3.
Transactions, American Geo-physical Union , 35(2):293, 1954. ISSN 0002-8606. doi: 10.1029/TR035i002p00293.URL .H Philip Crotwell, Thomas J Owens, and Jeroen Ritsema. The TauP Toolkit: Flexibleseismic travel-time and ray-path utilities.
Seismological Research Letters , 70(April):154–160, 1999. doi: 10.1785/gssrl.70.2.154.Cline Davy, Guilhem Barruol, Fabrice R. Fontaine, Karin Sigloch, and ElonoreStutzmann. Tracking major storms from microseismic and hydroacous-tic observations on the seafloor.
Geophysical Research Letters , 41(24):8825–8831, 2014. ISSN 00948276. doi: 10.1002/2014GL062319. URL http://doi.wiley.com/10.1002/2014GL062319 .Anja Diez.
Effects of cold glacier ice crystal anisotropy on seismic data . PhD thesis,Karlsruher Institut f¨ur Technologie, 2013.Adam M. Dziewo´nski and Don L. Anderson. Preliminary reference Earthmodel.
Physics of the Earth and Planetary Interiors , 25(4):297–356,6 1981. ISSN 00319201. doi: 10.1016/0031-9201(81)90046-7. URL http://linkinghub.elsevier.com/retrieve/pii/0031920181900467 .Gran Ekstr¨om, Adam M. Dziewonski, and Joseph M. Steim. Single station CMT; Ap-plication to the Michoacan, Mexico, Earthquake of September 19, 1985.
GeophysicalResearch Letters , 13(3):173–176, 1986. doi: 10.1029/GL013i003p00173.39aurice Ewing, Ivan Tolstoy, and Frank Press. Proposed use of the T phase in tsunamiwarning systems.
Bulletin of the Seismological Society of America , 40(1):53–58,1951.Ronni Grapenthin and Jeffrey T. Freymueller. The dynamics of a seismicwave field: Animation and analysis of kinematic GPS data recorded dur-ing the 2011 Tohoku-oki earthquake, Japan.
Geophysical Research Letters ,38(18), 9 2011. ISSN 00948276. doi: 10.1029/2011GL048405. URL http://doi.wiley.com/10.1029/2011GL048405 .K. P. Hand, A. E. Murray, J. B. Garvin, and and the Science Definition Team. Reportof the Europa Lander Science Definition Team. Technical report, NASA, Pasadena,California, 2017.Thomas Mejer Hansen, Andre G. Journel, Albert Tarantola, and Klaus Mosegaard.Linear inverse Gaussian theory and geostatistics.
Geophysics , 71(6):R101, 2006.ISSN 00168033. doi: 10.1190/1.2345195.Alexander G. Hayes. The Lakes and Seas of Titan.
Annual Re-view of Earth and Planetary Sciences , 44(1):57–83, 2016. ISSN0084-6597. doi: 10.1146/annurev-earth-060115-012247. URL .J. D. Hofgartner, A. G. Hayes, J. I. Lunine, H. Zebker, B. W. Stiles, C. Sotin,J. W. Barnes, E. P. Turtle, K. H. Baines, R. H. Brown, B. J. Buratti, R. N.Clark, P. Encrenaz, Randolph L. Kirk, A. Le Gall, R. M. Lopes, Ralph D.Lorenz, M. J. Malaska, K. L. Mitchell, P. D. Nicholson, P. Paillou, J. Rade-baugh, S. D. Wall, and C. Wood. Transient features in a Titan sea.
NatureGeoscience , 7(7):493–496, 6 2014. ISSN 1752-0894. doi: 10.1038/ngeo2190. URL .D.L. Hogenboom, J.S. Kargel, J.P. Ganasan, and L. Lee. Magnesium sulfate-water to400 MPa using a novel piezometer: Densities, phase equilibria, and planetologicalimplications.
Icarus , 115(2):258–277, 1995. doi: 10.1006/icar.1995.1096.L Iess, R A Jacobson, M Ducci, D J Stevenson, J I Lunine, J W Armstrong, S WAsmar, P Racioppa, N J Rappaport, and P Tortora. The tides of Titan.
Science ,337(6093):457–459, 2012.L Iess, D J Stevenson, M Parisi, D Hemingway, R A Jacobson, J I Lunine, F Nimmo,J W Armstrong, S W Asmar, M Ducci, and P Tortora. The Gravity Field andInterior Structure of Enceladus.
Science , 344(6179):78–80, 2014.Heiner Igel.
Computational seismology: a practical introduction . Oxford UniversityPress, Oxford, 1st edition, 2017. ISBN 9780198717416.R A Jacobson, P G Antreasian, J J Bordi, K E Criddle, R Ionasescu, J B Jones, R AMackenzie, M C Meek, D Parcher, F J Pelletier, and others. The Gravity Field of40he Saturnian System from Satellite Observations and Spacecraft Tracking Data.
The Astronomical Journal , 132(6):2520–2526, 2006.Simon A. Kattenhorn and Louise M. Prockter. Evidence for subduction in the iceshell of Europa.
Nature Geoscience , 7(10):762–767, 9 2014. ISSN 1752-0894. doi:10.1038/ngeo2245. URL .Sharon Kedar, Hiroo Kanamori, and Bradford Sturtevant. Bubble collapse as thesource of tremor at Old Faithful Geyser.
Journal of Geophysical Research B: SolidEarth , 103(10):24283 – 24229, 1998. doi: 10.1029/98JB01824.Sharon Kedar, M. S. Longuet-Higgins, F. Webb, N. Graham, Robert W. Clayton, andC. Jones. The origin of deep ocean microseisms in the North Atlantic Ocean.
Pro-ceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences ,464(2091):777–793, 2008. ISSN 1364-5021. doi: 10.1098/rspa.2007.0277.B. L. N. Kennett.
The Seismic Wavefield: Interpretation of seismograms on regionaland global scales . Cambridge University Press, 2002. ISBN 0521006651..Amir Khan, Martin van Driel, Maren B¨ose, D. Giardini, Savas Ceylan, J. Yan,John F. Clinton, F. Euchner, Philippe Lognonn´e, N. Murdoch, David Mi-moun, Mark P Panning, Martin Knapmeyer, and W.B. Banerdt. Single-station and single-event marsquake location and inversion for structure us-ing synthetic Martian waveforms.
Physics of the Earth and Planetary In-teriors , 7 2016. ISSN 00319201. doi: 10.1016/j.pepi.2016.05.017. URL http://linkinghub.elsevier.com/retrieve/pii/S0031920116300875 .M G Kivelson, K. K. Khurana, C T Russell, M Volwerk, R J Walker, C Zimmer,R. Greeley, R. T. Pappalardo, J. W. Head, G. V. Hoppa, B. R. Tufts, R. Green-berg, P. E. Geissler, M. G. Kivelson, D. S. Colburn, R. T. Reynolds, F. M.Neubauer, D. J. Southwood, R. J. Walker, J. A. Slavin, A. J. Kliore, D. P. Hin-son, F. M. Flasar, A. F. Nagy, T. E. Cravens, D. T. Hall, P. D. Feldman, M. A.McGrath, D. F. Strobel, and J. Saur. Galileo magnetometer measurements: astronger case for a subsurface ocean at Europa.
Science (New York, N.Y.) , 289(5483):1340–3, 2000. ISSN 0036-8075. doi: 10.1126/science.289.5483.1340. URL .M G Kivelson, K. K. Khurana, and M Volwerk. The Permanentand Inductive Magnetic Moments of Ganymede.
Icarus , 157(2):507–522, 6 2002. ISSN 00191035. doi: 10.1006/icar.2002.6834. URL .M. Knapmeyer, H. H. Fischer, J. Knollenberg, K. J. Seidensticker, K. Thiel,W. Arnold, W. Schmidt, C. Faber, F. Finke, and D. M¨ohlmann. TheSESAME/CASSE instrument listening to the MUPUS PEN insertion phaseon comet 67P/ChuryumovGerasimenko.
Acta Astronautica , 125:234–249,2016. ISSN 00945765. doi: 10.1016/j.actaastro.2016.02.018. URL http://dx.doi.org/10.1016/j.actaastro.2016.02.018 .41artin Knapmeyer. Syntax Diagrams for the Body Wave Nomenclature with ex-tensions for other terrestrial planets. In
Geophysical Research Abstracts, , vol. 5,EGS-AGU-EGU Joint Assembly , Nice, 2003.R. L. Kovach and Christopher F. Chyba. Seismic Detectabil-ity of a Subsurface Ocean on Europa.
Icarus , 150(2):279–287,4 2001. ISSN 00191035. doi: 10.1006/icar.2000.6577. URL http://linkinghub.elsevier.com/retrieve/pii/S0019103500965771 .Lion Krischer, Tobias Megies, Robert Barsch, Moritz Beyreuther, Thomas Lecocq,Corentin Caudron, and Joachim Wassermann. ObsPy: a bridge for seismologyinto the scientific Python ecosystem.
Computational Science & Discovery , 8(1):014003, 2015. ISSN 1749-4699. doi: 10.1088/1749-4699/8/1/014003. URL http://stacks.iop.org/1749-4699/8/i=1/a=014003?key=crossref.2a104192c817bd4a435deb73c5877c16 .Lion Krischer, Alexander R Hutko, Martin van Driel, Simon C. St¨ahler,Manochehr Bahavar, Chad Trabant, and Tarje NissenMeyer. OnDemand Cus-tom Broadband Synthetic Seismograms.
Seismological Research Letters , 88(4):1127–1140, 7 2017. ISSN 0895-0695. doi: 10.1785/0220160210. URL https://pubs.geoscienceworld.org/srl/article/88/4/1127-1140/354100 .S. Kugler, T. Bohlen, S. Bussat, and G. Klein. Variability of Scholte-wave Dispersionin Shallow-water Marine Sediments.
Journal of Environmental & Engineering Geo-physics , 10(2):203–218, 6 2005. ISSN 1083-1363. doi: 10.2113/JEEG10.2.203. URL http://jeeg.geoscienceworld.org/cgi/doi/10.2113/JEEG10.2.203 .Sunwoong Lee, Michele Zanolin, Aaron M Thode, Robert T Pappalardo, andNicholas C Makris. Probing Europas interior with natural sound sources.
Icarus ,165:144–167, 2003. doi: 10.1016/S0019-1035(03)00150-7.Axel Lef`evre, Gabriel Tobie, Gal Choblet, and Ondej ˇCadek. Structure and dynamicsof Titan’s outer icy shell constrained from Cassini data.
Icarus , 237:16–28, 2014.ISSN 10902643. doi: 10.1016/j.icarus.2014.04.006.Inge Lehmann. P’.
Bureau Central S´eismologique International, S´erie A, TravauxScientifiques , 14:87–115, 1936.Kuangdai Leng, Tarje Nissen-Meyer, and Martin van Driel. Efficient global wavepropagation adapted to 3-D structural complexity: a pseudospectral/spectral-element approach.
Geophysical Journal International , 207(3):1700–1721, 12 2016. ISSN 0956-540X. doi: 10.1093/gji/ggw363. URL http://gji.oxfordjournals.org/lookup/doi/10.1093/gji/ggw363 .J. S. Lewis. Satellites of the outer planets: Thermal models.
Science , 172:1171–1172,1971. doi: 10.1126/science.172.3988.1127.P. Lognonn´e, V.N. Zharkov, J.F. Karczewski, B Romanowicz, M Menvielle,G. Poupinet, B Brient, C Cavoit, A Desautez, B Dole, D Franqueville,42 Gagnepain-Beyneix, H Richard, P Schibler, and N Striebig. The seis-mic OPTIMISM experiment.
Planetary and Space Science , 46(6-7):739–747, 6 1998. ISSN 00320633. doi: 10.1016/S0032-0633(98)00009-9. URL http://linkinghub.elsevier.com/retrieve/pii/S0032063398000099 .P. Lognonn´e, F. Karakostas, L. Rolland, and Y. Nishikawa. Modeling of atmospheric-coupled Rayleigh waves on planets with atmosphere: From Earth observation toMars and Venus perspectives.
Journal of the Acoustical Society of America , 140(2),2016. ISSN 00014966. doi: 10.1121/1.4960788.Philippe Lognonne, W B Banerdt, K Hurst, David Mimoun, R Garcia, M Lefeuvre,J Gagnepain-Beyneix, M Wieczorek, A Mocquet, and Mark P Panning. INSIGHTand single-station broadband seismology: From signal and noise to interior struc-ture determination. In
Lunar and Planetary Institute Science Conference Abstracts ,volume 43, page 1983, 2012.Ralph D. Lorenz and Alexander G. Hayes. The growth of wind-waves in Titan’s hydrocarbon seas.
Icarus , 219(1):468–475,2012. ISSN 00191035. doi: 10.1016/j.icarus.2012.03.002. URL http://dx.doi.org/10.1016/j.icarus.2012.03.002 .Ralph D. Lorenz, Sharon Kedar, Naomi Murdoch, Philippe Lognonn??, Taichi Kawa-mura, David Mimoun, and W. Bruce Banerdt. Seismometer detection of dust devilvortices by ground tilt.
Bulletin of the Seismological Society of America , 105(6):3015–3023, 2015. ISSN 19433573. doi: 10.1785/0120150133.William B McKinnon. Effect of Enceladus’s rapid synchronous spin on in-terpretation of Cassini gravity.
Geophysical Research Letters , 42(7):2137–2143, 4 2015. ISSN 00948276. doi: 10.1002/2015GL063384. URL http://doi.wiley.com/10.1002/2015GL063384 .David Mimoun, P Lognonn´e, W B Banerdt, K Hurst, S Deraucourt, J Gagnepain-Beyneix, T Pike, S Calcutt, M Bierwirth, R Roll, P Zweifel, D Mance, O Robert,T N´ebut, S Tillier, P Laudet, L Kerjean, R Perez, D Giardini, U Christenssen,and R Garcia. The InSight SEIS Experiment. In
Lunar and Planetary ScienceConference , volume 43, page 1493, 2012.Andrija Mohoroviˇci´c. Das Beben vom 8. Oktober 1909 (translated title).
Jahrbuchdes meteorologischen Observatoriums in Zagreb , 4(1):1–67, 1909.F. Nimmo and R.T. Pappalardo. Ocean Worlds in the Outer Solar System.
Journal ofGeophysical Research: Planets , 2016. ISSN 21699097. doi: 10.1002/2016JE005081.URL http://doi.wiley.com/10.1002/2016JE005081 .F. Nimmo, P. C. Thomas, R. T. Pappalardo, and W. B. Moore. The global shape ofEuropa: Constraints on lateral shell thickness variations.
Icarus , 191(1):183–192,2007. ISSN 00191035. doi: 10.1016/j.icarus.2007.04.021.43arje Nissen-Meyer, Alexandre Fournier, and F. A. Dahlen. A two-dimensionalspectral-element method for computing spherical-earth seismogramsI. Moment-tensor source.
Geophysical Journal International , 168(3):10671092, 32007. ISSN 1365-246X. doi: 10.1111/j.1365-246X.2006.03121.x. URL http://onlinelibrary.wiley.com/doi/10.1111/j.1365-246X.2006.03121.x/abstracthttp://onlinelibrary.wiley.com/doi/10.1111/j.1365-246X.2006.03121.x/full .Tarje Nissen-Meyer, Alexandre Fournier, and F. A. Dahlen. A 2-D spectral-element method for computing spherical-earth seismograms-II. Waves insolid-fluid media.
Geophysical Journal International , 174(3):873–888, 92008. ISSN 0956540X. doi: 10.1111/j.1365-246X.2008.03813.x. URL http://onlinelibrary.wiley.com/doi/10.1111/j.1365-246X.2008.03813.x/abstracthttp://doi.wiley.com/10.1111/j.1365-246X.2008.03813.x .Tarje Nissen-Meyer, Martin van Driel, Simon C. St¨ahler, Kasra Hosseini, StefanieHempel, Ludwig Auer, Andrea Colombi, and Alexandre Fournier. AxiSEM:broadband 3-D seismic wavefields in axisymmetric media.
Solid Earth , 5(1):425–445, 6 2014. ISSN 1869-9529. doi: 10.5194/se-5-425-2014. URL .Guust Nolet and Leroy M. Dorman. Waveform Analysis of Scholte ModesIn Ocean Sediment Layers.
Geophysical Journal International , 125(2):385–396, 5 1996. ISSN 0956540X. doi: 10.1111/j.1365-246X.1996.tb00006.x. URL https://academic.oup.com/gji/article-lookup/doi/10.1111/j.1365-246X.1996.tb00006.x .Richard Dixon Oldham. The Constitution of the Interior of the Earth, as Re-vealed by Earthquakes.
Quarterly Journal of the Geological Society , 62(1-4):456–475, 2 1906. ISSN 0370-291X. doi: 10.1144/GSL.JGS.1906.062.01-04.21. URL http://jgslegacy.lyellcollection.org/cgi/doi/10.1144/GSL.JGS.1906.062.01-04.21 .Mark P Panning, V. Lekic, M. Manga, Fabio Cammarano, and B. Romanowicz. Long-period seismology on Europa: 2. Predicted seismic response.
Journal of GeophysicalResearch , 111(E12):E12008, 12 2006. ISSN 0148-0227. doi: 10.1029/2006JE002712.URL http://doi.wiley.com/10.1029/2006JE002712 .Mark P Panning, ric Beucler, Mlanie Drilleau, Antoine Mocquet, Philippe Lognonn´e,and W Bruce Banerdt. Verifying single-station seismic approaches using Earth-based data: Preparation for data return from the InSight mission to Mars.
Icarus ,248:230–242, 3 2015. ISSN 00191035. doi: 10.1016/j.icarus.2014.10.035. URL http://linkinghub.elsevier.com/retrieve/pii/S001910351400582X .Mark P Panning, Simon C. St¨ahler, Hsin-Hua Huang, Steve Vance, Sharon Kedar,Victor C. Tsai, and William T. Pike. The seismic noise environment of Europa.
Journal of Geophysical Research E: Planets , 2017.L. E. Peters, S. Anandakrishnan, R. B. Alley, and D. E. Voigt. Seismic attenuationin glacial ice: A proxy for englacial temperature.
Journal of Geophysical Research:Earth Surface , 117(2):1–10, 2012. ISSN 21699011. doi: 10.1029/2011JF002201.44on Peterson. Observations and Modeling of Seismic Background Noise. Technicalreport, USGS, Albuquerque, New Mexico, 1993.Cynthia B. Phillips and Robert T. Pappalardo. Europa Clipper Mission Concept: Ex-ploring Jupiter’s Ocean Moon.
Eos, Transactions American Geophysical Union ,95(20):165–167, 5 2014. ISSN 00963941. doi: 10.1002/2014EO200002. URL http://doi.wiley.com/10.1002/2014EO200002 .W T Pike, S Calcutt, I M Standley, A G Mukherjee, J Temple, T Warren, C Char-alambous, H Liu, A Stott, and J B McClean. A Silicon Seismic Package (SSP) forPlanetary Geophysics. In
Lunar and Planetary Science Conference , volume 47, page2081, 2016.C. C. Porco, P. Helfenstein, P. C. Thomas, A. P. Ingersoll, J. Wisdom, R. West,G. Neukum, T. Denk, R. Wagner, T. Roatsch, S. Kieffer, E. Turtle, A. McEwen,T. V. Johnson, J. Rathbun, J. Veverka, D. Wilson, J. Perry, J. Spitale, A. Brahic,J. A. Burns, A. D. DelGenio, L. Dones, C. D. Murray, and S. Squyres. CassiniObserves the Active South Pole of Enceladus.
Science , 311(5766), 2006.Frank Press and Maurice Ewing. Propagation of elastic waves in afloating ice sheet.
Transactions, American Geophysical Union , 32(5):673, 1951. ISSN 0002-8606. doi: 10.1029/TR032i005p00673. URL .Dieter Rauch. Experimental and Theoretical Studies of Seismic Interface Waves inCoastal Waters. In William A Kuperman and Finn B Jensen, editors,
Bottom-Interacting Ocean Acoustics , pages 307–327. Springer US, Boston, MA, 1980. ISBN978-1-4684-9051-0.A. T. Ringler and C. R. Hutt. Self-Noise Models of Seismic Instruments.
SeismologicalResearch Letters , 81(6):972–983, 2010. ISSN 0895-0695. doi: 10.1785/gssrl.81.6.972.J. G. Scholte. The range of existence of Rayleigh and Stone-ley waves.
Geophysical Journal International , 5(s5):120–126, 5 1947. doi: 10.1111/j.1365-246X.1947.tb00347.x. URL https://academic.oup.com/gsmnras/article-lookup/doi/10.1111/j.1365-246X.1947.tb00347.x .G Schubert, J D Anderson, T Spohn, and W B McKinnon. Interior composition,structure and dynamics of the Galilean satellites.
Jupiter: The Planet, Satellitesand Magnetosphere , pages 281–306, 2004.K. J. Seidensticker, D. M¨ohlmann, I. Apathy, W. Schmidt, K. Thiel, W. Arnold,H. H. Fischer, M. Kretschmer, D. Madlener, A. P´eter, R. Trautner, and S. Schieke.SESAME - An experiment of the Rosetta Lander Philae: Objectives and gen-eral design.
Space Science Reviews , 128(1-4):301–337, 2007. ISSN 00386308. doi:10.1007/s11214-006-9118-6. 45. M. Soderlund, B. E. Schmidt, J. Wicht, and D. D. Blankenship. Ocean-driven heating of Europas icy shell at low latitudes.
Nature Geoscience ,7(1):16–19, 2013. ISSN 1752-0894. doi: 10.1038/ngeo2021. URL .Simon C. St¨ahler, Karin Sigloch, Kasra Hosseini, Wayne C. Crawford, Guil-hem Barruol, Mechita C. Schmidt-Aursch, Maria Tsekhmistrenko, J.-R. Scholz,A Mazzullo, and M Deen. Performance report of the RHUM-RUM ocean bot-tom seismometer network around La R´eunion, western Indian Ocean.
Ad-vances in Geosciences , 41:43–63, 2016. doi: 10.5194/adgeo-41-43-2016. URL .Dmitry A. Storchak, Johannes Schweitzer, and Peter Bormann. The IASPEI stan-dard seismic phase list.
Seismological Research Letters , 74(6):761–772, 2003. URL http://srl.geoscienceworld.org/content/74/6/761.short .P C Thomas. Sizes, shapes, and derived properties of the saturnian satellites after theCassini nominal mission.
Icarus , 208(1):395–401, 2010.P.C. Thomas, R. Tajeddine, M.S. Tiscareno, J.A. Burns, J. Joseph,T.J. Loredo, P. Helfenstein, and C. C. Porco. Enceladuss measuredphysical libration requires a global subsurface ocean.
Icarus , 264:37–47, 2016. ISSN 00191035. doi: 10.1016/j.icarus.2015.08.037. URL http://linkinghub.elsevier.com/retrieve/pii/S0019103515003899 .Ivan Tolstoy and Maurice Ewing. The T phase of shallow-focus earthquakes.
Bulletin of the Seismological Society of America , 40(1):25–51, 1950. URL https://pubs.geoscienceworld.org/bssa/article-lookup/40/1/25 .Martin van Driel and Tarje Nissen-Meyer. Optimized viscoelastic wave propa-gation for weakly dissipative media.
Geophysical Journal International , 199(2):1078–1093, 9 2014a. ISSN 0956-540X. doi: 10.1093/gji/ggu314. URL http://gji.oxfordjournals.org/cgi/doi/10.1093/gji/ggu314 .Martin van Driel and Tarje Nissen-Meyer. Seismic wave propagation in fullyanisotropic axisymmetric media.
Geophysical Journal International , 199(2):880–893, 9 2014b. ISSN 0956-540X. doi: 10.1093/gji/ggu269. URL http://gji.oxfordjournals.org/cgi/doi/10.1093/gji/ggu269 .Martin van Driel, Lion Krischer, Simon C. St¨ahler, Kasra Hosseini, and Tarje Nissen-Meyer. Instaseis: instant global seismograms based on a broadband waveformdatabase.
Solid Earth , 6(2):701–717, 2015. ISSN 1869-9529. doi: 10.5194/se-6-701-2015. URL .Steven D Vance, Mathieu Bouffard, Mathieu Choukroun, and Christophe Sotin.Ganymede’s internal structure including thermodynamics of magnesium sulfateoceans in contact with ice.
Planetary and Space Science , 96:62–70, 2014. ISSN00320633. doi: 10.1016/j.pss.2014.03.011.46teven D Vance, Sharon Kedar, Mark P Panning, Simon C. Staehler, Bruce G Bills,Ralph D. Lorenz, Hsin-hua Huang, William T. Pike, Julie C Castillo, PhilippeLognonn´e, Victor C. Tsai, and Alyssa R. Rhoden. Vital Signs: Seismology of oceanworlds.
Astrobiology , in press, 10 2017a. URL http://arxiv.org/abs/1610.10067 .Steven D Vance, Mark P Panning, S. C. St¨ahler, Bruce G Bills, Fabio Cammarano,Victor C. Tsai, Hsin-Hua Huang, Sharon Kedar, Christophe Sotin, W.B. Banerdt,Sridhar Anandakrishnan, William T. Pike, and Ralph D. Lorenz. Geophysical testsfor habitability in ocean worlds.
Journal of Geophysical Research E: Planets , to besubm, 2017b.Sanjay Vijendran, Jeremy Fielding, Johan K¨ohler, Rob Gowen, Phil Church, and PeterFalkner. A penetrator for the Jupiter Ganymede Orbiter Mission. In
Proceedings ofthe 7th International Planetary Probe Workshop: Barcelona, Spain , 2010.C C Walker and B E Schmidt. Ice collapse over trapped water bodies on Enceladus andEuropa.
Geophysical Research Letters , 42(3):712–719, 2 2015. ISSN 00948276. doi:10.1002/2014GL062405. URL http://doi.wiley.com/10.1002/2014GL062405 .J H Woodhouse. The Calculation of Eigenfrequencies and Eigenfunctions of the FreeOscillations of the Earth and the Sun. In
Seismological algorithms , pages 321–370,1988.Yih-min Wu, Hsin-yi Yen, Li Zhao, Bor-shouh Huang, and Wen-tzong Liang. Magni-tude determination using initial P waves : A single-station approach.
GeophysicalResearch Letters , 33(March):4–7, 2006. doi: 10.1029/2005GL025395.Howard Zebker, Alex Hayes, Mike Janssen, Alice Le Gall, Ralph D. Lorenz, and Lau-ren Wye. Surface of Ligeia Mare, Titan, from Cassini altimeter and radiometeranalysis.
Geophysical Research Letters , 41(2):308–313, 1 2014. ISSN 00948276. doi:10.1002/2013GL058877. URL http://doi.wiley.com/10.1002/2013GL058877http://doi.wiley.com/10.1002/2013GL058877