Lucky-7 decoded

Lucky-7 is a Czech Amateur 1U cubesat that was launched on July 5 together with the meteorological Russian satellite Meteor M2-2 and several other small satellites in a Soyuz-2-1b Fregat-M rocket from Vostochny. The payload it carriers is rather interesting: a low power GPS receiver, a gamma ray spectrometer and dosimeter, and a photo camera.

It transmits telemetry in the 70cm Amateur satellite band using 4k8 FSK with an Si4463 transceiver chip. Additional details about the signal can be found here.

After some work trying to figure out the scrambler used by the Si4463, I have added a decoder for Lucky-7 to gr-satellites. This post shows some of the technical aspects of the decoder.

DSLWP-B activities for the third week of July

During this week, the Amateur payload of DSLWP-B was active during the following slots:

  • 14 Jul 19:00 to 21:00
  • 15 Jul 12:00 to 14:00
  • 17 Jul 04:40 to 06:50
  • 18 Jul 20:50 to 22:50
  • 20 Jul 14:20 to 16:20
  • 21 Jul 05:30 to 07:30

Among these, the Moon was visible from Europe only on July 14, 18 and 21, so Dwingeloo only observed these days, which were mainly devoted to the download of SSDV images of the lunar surface. As usual, the payload took an image automatically at the start of each slot, so some of the slots were used for autonomous lunar imaging, even though no tracking was made from Dwingeloo.

This post is a detailed account of the activities done with DSLWP-B during the third week of July.

Galileo operational again

In my previous post I wrote about the outage of the Galileo EU satellite navigation constellation that started on 2019-07-12 01:50 GST. For several days, all the Galileo satellites kept transmitting the same batch of ephemeris, which corresponded to 2019-07-11 21:50 GST, rather than sending updated batches of ephemeris.

This situation continued until around 2019-07-16 19:00, when several people reported that they were receiving updated ephemeris from some Galileo satellites. These ephemeris could be used for navigation correctly. The NAVSAS group from University of Torino has published a post where they show, for each satellite, when the updated ephemeris were first received by their equipment in Torino, Italy.

The restoration of the system was publicly announced with NAGU 2019027, published on 2019-07-18 08:20. This NAGU stated that starting from 2019-07-17 20:52 the service was restored, but users might experience some instability.

It is interesting to look at the MGEX broadcast ephemeris to see what happened between 2019-07-16 19:00, when users started to see updated ephemeris, and 2019-07-18, when the system was considered to be fully restored. In this post, I do a detailed analysis of the broadcast ephemeris.

Galileo constellation outage

You may have already heard that the Galileo EU satellite navigation constellation is out of service since last Thursday July 11th. Currently the Galileo constellation status page shows that the status of most Galileo satellites is not usable because of a service outage. The satellites not affected by this outage are E20, which was made unavailable on 2014-05-27, and E14 and E18, the Galileo eccentric satellites, which failed to achieve the circular nominal orbit and are only used for testing purposes.

The European GNSS Agency has given very little information regarding the causes of the problem. The available information boils down to NAGU 2019026, posted on 2019-07-13 20:15, which states that starting from 2019-07-12 01:50 the Galileo signals should not be used.

This has originated many rumours and confusion about the problem. It seems that the major cause was a failure in the Precise Timing Facility, which is in charge of the generation of a realization of the GST, the timescale used by Galileo. This has affected the OSPF, which is the service that generates the orbit and clock products (ephemeris) for the satellites. Thus, since Thursday, no new ephemeris are being computed and uploaded to the satellites.

Taking rumours aside, in this post I will look at some facts about the outage which can be learned by analizing the Galileo signal. Other people have published similar studies, such as the NAVSAS group from University of Torino.

DSLWP-B activities for the second week of July

During this week, the Amateur payload of DSLWP-B has been active on the following slots:

  • 2019-07-09 12:30 to 14:30 UTC
  • 2019-07-10 13:50 to 15:50 UTC
  • 2019-07-12 02:00 to 04:00 UTC
  • 2019-07-12 16:30 to 18:30 UTC

In these activations, the last remaining image of the eclipse was downloaded. Also, images of the lunar surface as well as stars have been taken and downloaded. This post is a summary of the activities made with DSLWP-B over the week.

Lunar features in DSLWP-B eclipse images

Now that the DSLWP-B eclipse images have become widespread, appearing even in some newspapers, I have taken the time to identify the features of the lunar surface that can be seen in two of the images. As with any Earth and Moon pictures of DSLWP-B, the part of the Moon that can be seen in these images belongs to the far side, with longitudes of approximately 100º E and 100º W (the division between the near side and the far side happens at 90º E and 90º W, and 0º is the middle of the near side).

I have compared the images with a simulation of the camera view done in GMAT. Using this simulation as a reference and these lunar surface maps as well as Google Moon, I have labelled the features that are visible in the images. The 4K lunar surface map from Celestia Motherlode has been used in GMAT, instead of the default, lower resolution map.

Occultation entry

The figure below shows the GMAT camera view simulation for one of the images taken as DSLWP-B was hiding behind the Moon. The field of view in this figure is much larger than the field of view of the DSLWP-B Inory eye camera. The up direction is the normal to DSLWP-B’s orbit around the Sun (defined as the plane containing the position and velocity vectors with DLSWP-B with respect to the Sun). Therefore, it points approximately towards the north pole. DSWLP-B is moving towards the upper right corner of this image.

GMAT camera view simulation for image 0xE2

The figure belows shows the ground track view. The satellite has crossed the near side and is now starting to orbit over the far side, soon becoming hidden behind the Moon.

GMAT ground track view for image 0xE2

Image 0xE2, taken on 2019-07-02 18:57:20 is shown below. The image has been rotated to match the orientation of the GMAT camera view simulation.

Image 0xE2, rotated

The craters that are visible in this image are labeled in the figure below. These belong to the Mare Australe region, and several well known craters such as Jenner and Lamb can be seen.

Image 0xE2, rotated and labelled

Tammo Jan Dijkema has done a similar exercise with image 0xE3, which was taken a minute after the image shown above. DSLWP-B has moved towards the upper right corner of the image, so that a larger portion of the lunar surface is visible.

Occultation exit

The figure below shows the GMAT camera view simulation corresponding to image 0xE5, which was taken shortly after DSLWP-B exited the occultation, so that the Earth was visible again.

GMAT camera view simulation for image 0xE5

Below we see the ground track corresponding to this image. The satellite has crossed the far side of the Moon in a south to north direction and soon will cross over to the near side.

GMAT ground track view for image OxE5

The figure below is image 0xE5, which was taken on 2019-07-02 19:33:05. It has been rotated to match the orientation of the GMAT simulation.

Image 0xE5, rotated

The craters are labelled in the figure below. Important craters in the Coulomb-Sarton basin such as Stefan and Bragg are visible. An image of this region with the craters labelled can be seen here.

Image 0xE5, rotated and labelled

Results of DSLWP-B eclipse imaging

Probably you’ve already seen the amazing images of the July 2 solar eclipse taken by DSLWP-B. These have have been shared in Twitter by the Harbin Institute of Technology Amateur Radio club BY2HIT, Tammo Jan Dijkema from the Dwingeloo radiotelescope in the Netherlands, and others.

I have devoted the last few posts to the planning of the imaging times for the eclipse test run and the validation of the test run done on June 30. This post is a detailed account of the results of the eclipse imaging. Between July 3 and 5, five of the six eclipse images taken on July 2 were downloaded, as well as some other images taken later. Here I give a summary of the downloads during these days and compare the images to the predictions I made.

Analysis of DSLWP-B eclipse test run (again)

Yesterday I did an analysis of the image taken and downloaded during the test run for the imaging of the solar eclipse shadow on the Earth with DSLWP-B on July 2. Only one of the six images taken to validate the camera exposure and pointing and orbit errors was downloaded yesterday. The window to download the remaining images was during today’s activation of the Amateur payload between 05:30 and 07:30 UTC.

Three images were downloaded today, with a few missing chunks. Moreover, Wei Mingchuan BG2BHC has sent me better ephemeris, which give a more accurate prediction of the orbit both on the June 30 test run and on the eclipse imaging run on July 2. This post is an analysis of the images downloaded today.

Analysis of DSLWP-B eclipse test run

As described in one of my latest posts, today DSLWP-B has made a test imaging run in preparation for the solar eclipse on July 2. A series of images was taken just before the Moon hid the centre of the camera field of view and just after the Moon left the centre of the image, in approximately the same relative positions as for the July 2 eclipse imaging run.

The activation of the Amateur payload started at 05:30 UTC and the payload was commanded to change the configuration of the camera to use 2x zoom. The satellite was occulted by the Moon at 05:40 UTC, preventing the reception of telemetry until it reappeared at 06:16 UTC. The first series of images was taken automatically between 05:51 and 05:54, with the satellite behind the Moon.

After the satellite reappeared from behind the Moon, telemetry confirmed that three images, with IDs 0xD9, 0xDA and 0xDB had been taken. Between 06:29 and and 06:32, the satellite took the second series of images. Telemetry confirmed that these images were taken correctly with IDs 0xDC, 0xDD and 0xDE.

The priority was to use the rest of the activation to download images 0xDA and 0xDD, taken respectively at 05:52:40 and 06:31:10 UTC (when reading the times given in the planning post, note that the times listed there are the moments when the command is sent to the payload by the satellite, but the payload needs about 20 additional seconds to take the image). However, there were difficulties in commanding the payload, so half an hour was lost trying to command the satellite and only image 0xDA could be downloaded before the payload went off at 07:30 UTC.

Image 0xDA was downloaded without errors in a single transmission. It is shown below.

Image 0xDA, taken on 2019-06-30 05:52:40, downloaded between 07:11 ad 07:26

As we see, the exposure of the image is correct, so this image validates that the camera configuration can be used for the eclipse imaging run. Additionally, the image can be used to evaluate camera pointing and ephemeris errors.

As computed in this Jupyter notebook, the separation between the Moon rim and the centre of the field of view in the image shown above is 6.47 degrees. Using the camera calculations Jupyter notebook that I have shown in previous posts, we see that, according to the 20190630 ephemeris from dslwp_dev and my GMAT calculations, the angular distance between the Moon rim and the centre of the image at 05:52:40 UTC should be 3.25 degrees, assuming that the camera points perfectly away from the Sun.

The rate at which the Moon rim moves through the field of view is approximately 0.029 degrees per second. Thus, if the camera was pointing perfectly away from the Sun, this would indicate that DSLWP-B is 110 seconds earlier in its orbit that what predicted by the ephemeris, so that events concerning the relative position of the satellite and the Moon happen 110 seconds later than predicted.

However, one should take these calculations with a grain of salt. In my astrometry post, I showed that the camera was pointing 3.25 degrees off-axis. Therefore, it is convenient to assume an error of +/-3.25 degrees in the angle measurement done with the image. In units of time, this is +/-111 seconds.

So the data seems to suggest that DSLWP-B is one or two minutes earlier in its orbit and that the imaging times should be compensated by making them one or two minutes later, but there is not enough statistical evidence to support this argument. It will be very interesting to see image 0xDD, which will be downloaded tomorrow. The analysis of this image will give additional data.

In any case, so far it seems that orbit and pointing errors are within the tolerance given by the series of three images, which are taken at -1, 0, and +1 minutes offset from the nominal imaging time computed by Wei.

GPS SVN49 broadcasting non-standard codes again

As a GNSS engineer at my day job in GMV, it’s not uncommon to find myself looking at spectrums of the GNSS signal bands, either on a spectrum analyzer or on an IQ recording done by an SDR. A few days ago, I spotted something in the L1 band (1575.42MHz) that quickly caught my eye: a pair of strong carriers at +/- 511.5kHz away from the L1 centre frequency. This was visible on some of the recordings that we had done in the last few days, and also live on a spectrum analyzer.

When monitoring the signal on the spectrum analyzer, it disappeared a few hours later, making me suspect that it was transmitted by a satellite rather than a local interferer. Back at home, I did some recordings with STRF to try to identify the satellite using its Doppler.

The Doppler signature of the signal was a perfect match for GPS SVN49, also as USA-203 and GPS IIRM-7. This satellite, launched in 2009, was the first to demonstrate the L5 signal. During in-orbit testing, an anomaly with its navigation signals caused by the L5 filter was discovered. The satellite was never put into operational use, and has been used for varied tests ever since.

After searching information about this satellite, I learned that some researchers from Torino, Italy, had already observed back in 2017 the same kind of signals that I was seeing.

This post is a detailed study of the L1 signal that is currently being broadcast by GPS SVN49. The data used here has been published in Zenodo as “RF recordings of GPS SVN49 broadcasting non-standard codes“.