Wednesday, January 17, 2018

OTHR on 6950 kHz

TDoA maps 6950 KHz 20180112T2138Z 

Cross-correlations 6950 kHz 20180112T2138Z 
The bandwidth of this signal is larger than the 12kHz seen by the KiwiSDR. The auto-correlation for this signal has peaks at multiples of about 0.0462938428 sec but it is interesting that the pulses look like this:
abs(z) vs. time 6950 kHz 20180112T2138Z
so an interesting form of modulation is being used.

Sunday, January 14, 2018

Cyprus OTHR

Yesterday I found the OTHR from Cyprus operating with 25 Hz pulse repetition rate and μ=500kHz/s sweep rate. The most likely geo-location is just north of the known location:

TDoA plots 17030 kHz 20180113T1126Z

Cross-correlations 173030 kHz 20180113T1126Z
I was thinking about an alternative way of geo-location using an extended Kalman filter, and found this reference: It is written very well and I was able to implement it quickly in octave using dfpdp from the optim toolbox for the Jacobian. More on this later.

Friday, January 5, 2018

De-chirping applied to a HF over the horizon radar signal (2)

This is a follow up to this blog post. Again the signals of the OTHR in Cyprus were analyzed using GPS time-stamped IQ samples from several KiwiSDRs.

In this case the OTHR signal was found between 13965 and 13985 kHz. The slope of the FMCW signal is μ=1MHz/s, and there are 50 sweeps per second. This results in a bandwidth of about 20kHz.

As is well known, the frequencies Δf of the spectrum of the de-chirped signal correspond to time delays τ=-Δf/μ in the case of a stationary target.
  • for this one has to multiply the received signal with the original waveform
  • tuning offsets w.r.t. center frequency of the OTHR signal correspond to a time delay of the original waveform
  • when Δf is positive the delay corresponds to the frequency Δf-fs (fs is the sampling frequency), i.e. the ranges wrap around the x-axis.
Applying a 2nd Fourier transform to the time direction of time vs. Δf data one obtains the Doppler frequency fd vs. Δf where fd = v•λ/2 (for velocity v and wavelength λ). Here one has to be careful to adjust the phases in each bin to correct for the fact that the sampling times are not exactly aligned to multiples of 0.02 seconds.

Plots for 5 KiwiSDR are shown below. The splitting of the ionospheric reflection is most likely due to O- and X-mode propagation in the F layer.

13975 kHz 20180105T1114Z @CS5SEL; the band at fd=0 is caused by an interference

13975 kHz 20180105T1114Z @DF0KL

13975 kHz 20180105T1114Z @G8JNJ

13975 kHz 20180105T1114Z @Izhvesk

13975 kHz 20180105T1114Z @Julisdalen
Thanks to all KiwiSDR owners who make their receivers available and do connect a GPS antenna.

Sunday, December 17, 2017

TDoA measurements using GPS time-stamped IQ samples from KiwiSDRs (8)

Today I picked up a HF radar signal on 13500+-20 kHz. These are the time-differences between four KiwiSDRs:
TDoA 13490 kHz

This is most likely the Calypso CODAR radar system:
TDoA maps for 14390 kHz
Within the TDoA resolution the three sites in Ta'Sopu (Gozo), Ta'Barkat (Malta) and Pozzallo Harbour (Sicily) cannot be distinguished.

It is still puzzling for me that the TDoA maps seem to provide reasonable results, as for now great-circle distances are used, i.e., delays due to skywave propagation are not taken into account.

Thursday, December 14, 2017

De-chirping applied to a HF over the horizon radar signal

Here is another look at the HF OTHR signal from Cyprus already mentioned in this blog post.

This OTHR signal is (at least to 1st order) FMCW modulated (period 0.02 sec, bandwidth 20kHz), as can be seen in the FM demodulated phase which follows a sawtooth pattern with period 0.02 sec.

Although the bandwidth of the KiwiSDR IQ data stream is smaller than the bandwidth of the OTHR signal it turns out that de-chirping can be done: using KiwiSDR GPS time-stamps t, de-chirping can be implemented by multiplying the IQ samples by exp(-2i*pi*f(t)) with f=@(x) 0.5*1e6*(mod(x,0.02)-0.01).**2-50;  (this is octave notation)

The de-chirped signal looks like in the bottom panel below and is flat as expected:
OTHR FM demodulated phase; bottom: de-chirped

The relative range (which includes a contribution from the doppler effect, see, e.g., here) is then obtained by performing a FFT transformation on the de-chirped signal in each 0.02 sec period (~240 samples):
OTHR de-chirped spectrum vs. time

  • There is a main reflection in the ionosphere which has an interesting, time-varying structure to the left of it, i.e., towards smaller ranges
  • In the light blue regions there is a kind of digital pattern visible; this might be due to the signal processing and limited bandwidth of the KiwiSDR but it could also be caused by some additional modulation of the radar signal

Monday, December 11, 2017

TDoA measurements using GPS time-stamped IQ samples from KiwiSDRs (7)

14370 kHz OTHR 20171211T1040

The autocorrelation function has peaks at multiples of 0.02 seconds: (as is well known)
14370 kHz OTHR autocorrelation
The pulse shapes and the FM demodulated phases are shown below w.r.t. GPS seconds mod 0.02 seconds. Of the full bandwidth of this signal is 20kHz the KiwiSDRs sample 12 kHz. Therefore only part of the chirp signal is seen:
14370 kHz OTHR pulse shape and FM demodulated phase
Correcting for the sample rate the chirpyness becomes 12001*83.326 = 1e6, i.e., one sweep covers 20 kHz.

In the time differences different paths of propagation can be seen; it looks like the ionosphere at this time was quite patchy:
TDoA 14370 kHz

The TDoA location is consistent with the known position of the transmitter:
TDoA maps 14370 kHz

TDoA measurements using GPS time-stamped IQ samples from KiwiSDRs (6)

6030 kHz Radio Romania International in DRM

Time differences:
TDoA 6030 kHz

Note that there are two propagation paths with different delays observed in Khimki. Using the delays indicated in red the following maps are obtained:
TDoA maps 6030 kHz
For the other path the maps look like this:
TDoA maps 6030 kHz
Both are compatible with the known position of the transmitter.