Free oscillations for modern interior structure models of the Moon

1. T.V. Gudkova S. Raevskiy Schmidt Institute of Physics of the Earth RAS Free oscillations for modern interior structuremodels of the Moon

2. Crust thickness 60 km (Nakamura et al.,1976; Goins et al.,1977) 30-45 km( reprocessing or inversion of the Apollo seismic data) (Chenet et al., 2000; Khan et al., 2000; Vinnik et al., 2001;Gagnepain-Beyneix et al. 2004; Chenet et al., 2004) ≈ 20 – 110 km (Hood and Zuber, 2000) Core radius of the Moon Seismic data 170-360 km(Nakamura et al., 1974) Data on gravitational field and topography 220-450 km( Konopliv et al., 1998) Data obtained when the Moon passed through the Earth’s magnetospheric tail 340±90 km(Hood et al., 1999) Seismic lunar modelsI/MR2=0.3931±0.0002 “Lunar Prospector” (Konopliv et al., 1998) • Solid inner core 240 km • Liquid outer core 330 km • Partial melt up to 480 km Weber et al. (2011)

3. Interior structure models used for calculations The models fit both geodetic (lunar mass, polar moment of inertia, and Love numbers) and seismological (body wave arrivals measured by Apollo network) data. MW model (dashed line) : Weber et al. (2011) reanalyzed Apollo lunar seismograms using array processing mehods to search for the presence of reflected and converted seismic energy from the core. The results suggest the presence of a solid inner (240 km) and fluid outer (330 km, 8 g/cm3) core, overlain by a partially molten boundary layer (150 km thick). MG model (solid line) : Garsia et al. (2011) have constructed the Very Preliminary Reference Moon model and estimated the core radius by detecting core reflected S wave arrivals from waveform stacking methods. It has a core radius of 38040 km and an average core mass density of 5.21.0 g/cm3.

4. Free Oscillations Free oscillations, if they are excited, are particularly attractive to probe beneath the surface of an extraterrestrial body into its deep interior.Interpretation of data on free oscillations does not require knowledge of the time or location of the source; thus, data from a single station are sufficient. Since the planet has finite dimensions and is bounded by a free surface, the study of the free oscillations is based on the theory of vibration of an elastic sphere. The planet reacts to a quake (or an impact) by vibrating as a whole, vibrations being the sum of an infinite number of modes that correspond to a set of frequencies.

5. The free oscillations are divided into two types:(a) Torsional oscillations, whose displacement vector is perpendicular to the radius of a sphere - nTl(b) Spheroidal oscillations, whose displacement vector has components in both the radial and azimuthal direction - nSl The radial functionsnWlof torsional andnUlspheroidal oscillations depend on two indicesnandl, n - number of the overtone (the number of nodes along the radius in functionsUandW), l– the degree of the oscillation.

6. Functions 0Wl proportional to the displacements for torsional oscillations for the fundamental mode,l=2 to 10 vesus relative radius r/R. The important feature of free oscillations is that they concentrate towards the surface with increasing the degree l . Therefore different regions of interiors are sounded by different frequency intervals. The fundamental modes sound to those depth in the interiors where its displacement  0.3 The horizontal line drawn at level =0.3 enables one to judge graphically which modes give information about one or another zone of the planet. 0Wlis normalized to unity at the surface.

7. The relative period difference for models MW and MG is about 1.5% for basic tones and about 2-4% for overtones of torsional oscillations. For spheroidal oscillations it is about 2% for basic tones and up to 10% for overtones. If these modes are detected, we could refine the models. Relative period difference Т/Т(%) as a function of the oscillation number forfundamental mode (solid line) and two first overtones (dashed line- 1 and dot-dashed line - 2) of torsional (а) and spheroidal oscillations (b) for models by Weber et al. (2011)and Garsia et al. (2011).

8. The effect of the inner core rigidity on the structure of oscillationsThe model MW has the ‘core oscillation’ - FC. As the rigidity of the inner core increases, its period (42.96 min at μ=0) decreases up to 6.94 min at μ=4.23x1011 dyne/cm2, and its amplitude covers the mantle. At μ=0.5x1011 dyne/cm2 it looks like a regular oscillatiob R with a period of about 16 min. Besides the regular oscillation and its overtones O1,O2, ... , there are ‘inner core’ oscillations, with the energy localised in the core. The curves are not crossed, they change a slope and a type of oscillation. Spheroidal oscillations, l=2. Period T as a function of shear modulus of the inner core (for the model by Weber et al. (2011)  =4.23х1011dyn/cm2) . R – regular oscillations, FC – basic core oscillation, IC – oscillations of the inner core and О – overtones of regular oscillations.

9. The transition of FC type oscillation to R type and R type to FC type for the model by Weber et al. (2011), spheroidal oscillations, l=2 As the lunar core is rather small, the period difference of oscillations for the regular oscillation and the first overtone, for the models with inner core and without it, is very small (0.1 and 0.7%).It increases with the overtone number, and reachs 5-10% for the second and third overtones. The rigidity of the inner core influences mostly the periods of core oscillations, but their amplitudes are very small at the surface.

10. The amplitude spectrum for fundamental tones of torsional oscillations for l= 2, 3, 5, 7, 10, 20, 30, 40, 50 and their first overtones. M0=1: uN ≈10-29-10-26 cm M0=1022dyn cm: uN ≈10-7-10-4 cm l7-10 (up to 500 km) M0=1023dyn cm: uN ≈10-6-10-3 cm l5 (up to 750 km) X-axis - frequency =2/T and the period T (min), Y-axis - the amplitude of the horizontal component uN (cm). M0=1. Focal mechanism is 45o, 45o, 45o for dip, strike and slip angles with a focal depth of 50 km. The epicenter coordinates are 29oN, 98oW, the seismometer coordinates are 0oN, 45oW, the epicentral distance is 58o.

11. The amplitudes of displacements uN , uE and uR for fundamental modes of spheroidal oscillations with l=2-20. M0=1. Focal mechanism is 45o,45o,45o for dip, strike and slip angles at focal depths of 50 km (solid lines), 150 km (dot-dashed lines), 200 km (dashed lines). The epicenter coordinates are 29oN, 98oW, the seismometer coordinates are 0oN, 45oW, the epicentral distance is 58o. For =0.02s-1: uR ≈ 6x10-30x1.6x1022 aR ≈ 0.4x10-10 two orders lower than the value that can be recorded

12. CONCLUSION Free oscillations, if they are excited, are indeed particularly attractive to probe beneath the surface of the Moon into its deep interiors. The spectrum of torsional modes nTl wouldallow noticeable progress to be made in constructing aglobal model of the Moon’s interior structure (up to about 500 - 700 km depth), it was shown that the torsionalmodes with l > 5-7 can be recorded with current instruments. The seismic events on theMoon detected so far are too weak to excite freeoscillations that could be recorded. The accuracy of seismometers has been improved, and the application of new methods for processing seismic data let us hope to identify harmonics with smaller amplitudes in the future.