Lunar mapping with 6-meter wavelength
One of the efforts we are working on is mapping the Moon with 6-meter wavelength. This is done using the Jicamarca Radio Observatory 6-meter wavelength radar. The radar is useful for a number of reasons. The radar is huge. The dimensions of the antenna is 300 x 300 meters. It's got 9216 dipole antennas chained as groups of 12 colinear coaxial antennas. It can transmit up to 4 MW of power towards the Moon. It can be split into smaller sections for interferometry on receive. The antenna is fully polarimetric: you can transmit and receive any combination of two linear orthogonal polarizations.
The most important reason for using Jicamarca is frequency. The low frequency of 49.92 MHz allows the radio wave to penetrate very deep into the subsurface of the Moon. We estimate that in highland regions, the wave can penetrate up to 1 km beneath the surface. The low frequency does come with an associated cost. It is not possible to make a very high resolution map of the Moon. We can only observe 10-15 minute passes, and the Moon only has a Doppler width of about 1 Hz at this frequency. The low frequency is also severely affected by ionospheric radio propagation, which modifies the polarization state of the wave, and contains phase noise due to small fluctuations in ionospheric electron density.
This summer, we published our first paper on the topic in the journal Icarus containing first results. Now we are planning to write a follow up study that goes more in depth. In order to compare our results with other observations of the Moon, we have been working towards projecting our range-Doppler map into lunar coordinates. This is now finally done, and we can now project our observations into any projection we want! The gif animation above shows the Moon projected into orthographic projection using Basemap and Matplotlib, which has a nifty warpimage function.
fig = plt.figure() fig.patch.set_facecolor('black') b = Basemap(projection='ortho', lat_0 = lat_0, lon_0 = lon_0, resolution = 'l', area_thresh = 1000.) # plot image with cylindrical projection onto basemap scene b.warpimage(image=compf) # draw the edge of the map b.drawmapboundary(color="white") # draw lat/lon grid lines every 30 degrees. b.drawmeridians(np.arange(0, 360, 30),color="black") b.drawparallels(np.arange(-90, 90, 30),color="black") plt.title("JRO 6-m depolarized return",color="white") print(ofname) # ensure that black color is used for the background plt.savefig(ofname,facecolor=fig.get_facecolor()) plt.close()