Modernized Solar Radio Spectrograph in the L Band Based on Software Deﬁned Radio

: The paper presents the concept, implementation, and test operation of a modernized solar radio spectrograph for an investigation of the solar emission and solar bursts in radio frequency bands. Besides having a strong diagnostic signiﬁcance for studying the ﬂare energy release, the solar radio bursts can also cause strong interference for radio communication and navigation systems. The current spectrograph for the Ondrejov observatory (Astronomical Institute of Czech Academy of Sciences) was modernized by using a direct-conversion receiver connected to a ﬁeld-programmable gate array (FPGA) for the fast Fourier transform (FFT) spectrum estimation and put into the test operation. The higher time and frequency resolution and lower noise in comparison with the existing analog instrument were reached by the implementation of the latest optimal signal processing methods. To reduce the costs for such modernization, the operating frequency range was divided into four sub-bands of bandwidth 250 MHz, which brings another beneﬁt of greater scalability. The ﬁrst observations obtained by the new spectrograph and their comparison with the analog device are presented in the paper with future steps to put the spectrograph into the regular operation.


Radio Diagnostics in Solar Flare Research
Solar eruptions and flares represent the most powerful energy release processes in the entire solar system. Disturbances connected with these phenomena (magnetic clouds, beams of accelerated particles, enhanced radiation-namely in X-ray, UV, and radio bands) propagating from the Sun through the interplanetary space, are now known as effects of the space weather. In the case of reaching the Earth and its vicinity, they represent a serious risk for technologies (e.g., enhanced geo-currents-a risk for power distribution grids, ionospheric/TEC disturbances-an impact to the navigation systems and communications [1]) and also for humans (astronauts in the open space, aircraft crews and passengers for near-pole flights, etc.). In order to be able to give deterministic predictions of the space weather, we need to understand its drivers first. As the eruptive flares belong to those most effective, a research in the solar-flare physics represents a key to the future space weather forecasts.
Among all the radiation produced by solar flares the radio emission in decimetric range is the most closely related to the plasma processes of the flare energy release. Indeed, a radio emission connected with the beams of particles accelerated during the flare via a magnetic reconnection represents the first "messenger" of the solar flare ever [2]. Radio spectrographs, which are capable to record the dynamics of spectra on very short time scales, thus represent a key tool for remote radio diagnostics of the processes that started the flare. Modern approaches to the flare energy release via a magnetic reconnection [3] involve the energy cascading and fragmentation of the magnetic field structures in the solar atmosphere towards small scales. The small-scale structures, which are expected (as indicated by our numerical simulations) to be the actual locations of an energy release in flares, have also short dynamic timescales (millisecond and shorter). In order to be able to observe their radio signatures, the radio spectrograph should have a comparable temporal resolution. This represents a very strong driver for the spectrograph upgrade from the point of view of the astrophysical research.

Direct Influence of the Solar Radio Bursts on the Radio Systems
In addition to its important diagnostic role in the solar flare research, the solar radio emission (especially high-energy solar radio bursts) has a big impact on the navigation, communication, mobile, and Internet of things (IoT) systems and in a wider consequence to current and future information and transportation systems as well [4][5][6][7][8][9]. The influence of solar emissions can be observed as an increase of the effective noise temperature (or relevant noise factor) of the receiver as described in [10].The statistics of the solar radio bursts was published for example in [11] using a Poisson distribution for a description of the solar radio burst events together with an estimation of the distribution parameters from the long-term observation data. It can serve as a source for an estimation of the probability of high intensity events and their possible threats to the radio systems.
The influence of a solar radio burst to the popular satellite navigation system GPS (global positioning system) was documented in [4,6,7,12]. The unique radio burst from December 2006 with the peak level of one million SFU (solar flux unit) caused a significant reduction of the signal to noise ratio of the most GPS receivers, some receivers had lost tracking of the satellites and the acquisition of the new satellites was not possible due to the higher noise background. However, the performance degradation can occur for lower levels of radio bursts as well, for example a solar radio burst with the intensity 100,000 SFU causes in a common GPS receiver with a typical noise temperature 150 K and antenna gain 3 dB a drop of signal to noise ratio (SNR) by 15 dB that can prevent the receiver to acquire new signals [13].
The performance degradation can be also expected for the satellite-based augmentation systems (SBAS) [6]. These systems were developed as an improvement of the precision and safety of the current satellite navigation systems [14]. The precision improvement is realized by the dissemination of wide range differential corrections while the safety is ensured by the provision of integrity information. The information about the navigation system failure must be detected within 6s; therefore, the drop of the message reception probability is critical. The probability of a reception of the message depends, inter alia, on a signal to noise ratio and a bit rate, which is 10 times higher than the GPS message bit rate. Thus, the SBAS receiver requires a much higher signal to noise ratio for proper reception of the navigation messages than the GPS [14], so the aforementioned SNR drop can surely endanger the integrity of the satellite navigation system. Even though the observations of a degradation of the SNR of GPS and SBAS caused by the solar bursts were documented in [6], the detailed information about the impact of observed phenomena onto the system performance is however missing and should be an object of a future research.

The Motivation for Modernization of Solar Radio Spectrography Instruments
The solar radio emission in L band is in the research focus of the Astronomical Institute of the Czech Academy of Science in Ondrejov [15] for many years. The current scientific instrumentation RT5 ( Figure 1) is more than 25 years old. Its performance and architecture reflect the time of development. The spectrograph based on an analog technology is in the operation since 1993 without significant design changes. Its performance like a resolution in time and frequency as well as its dynamic range and linearity is currently below the current requirements demanded by both the advance in technology and the increasing level of knowledge in the solar plasma physics. The aim of the modernization project is to develop a modern but cost-effective instrument comparable or superior to the equipment of other observatories. The modern telescopes are usually designed as wideband radio receivers with a digital signal processing [16,17] and their technology enables us to estimate the frequency spectrum in parallel for all frequency bins in comparison with the analog technology that is based on a scanning principle.
The analog spectrum estimation is usually based on the application of a super heterodyne receiver that sequentially scans the frequency spectrum channel-by-channel and measures the signal power [18]. The modern digital methods are based on an application of the discrete Fourier transform (DFT) or computationally optimized fast Fourier transform (FFT) [19][20][21]. The common method for a consequent spectrum estimation based on the Fourier transform is a short-time Fourier transform (STFT) [22,23].
There are several possible ways how to implement such spectrum computation algorithms. Therefore, the critical design analysis has to be made to choose an optimal design not only to fulfill the system requirements but also to respect implementation issues and realization costs.
The observations and records from the modernized spectrograph will be used both for a study of solar processes and for an investigation of the impact of a solar radio emission onto the current and future communication and radio navigation systems in the frame of several research projects.

Technical Requirements
The technical requirements on the solar radio spectrograph were defined on the base of the long-term experiments with the current spectrograph with respect to the parameters of the similar instruments [17,18,24,25].
The main attention during the conceptual design was paid to a high sensitivity and a wide dynamic range. The secondary design parameters superseding the parameters of the current analog instrument [15] were:

•
Time resolution ≤1 ms. The aim of the modernization project is to develop a modern but cost-effective instrument comparable or superior to the equipment of other observatories. The modern telescopes are usually designed as wideband radio receivers with a digital signal processing [16,17] and their technology enables us to estimate the frequency spectrum in parallel for all frequency bins in comparison with the analog technology that is based on a scanning principle.
The analog spectrum estimation is usually based on the application of a super heterodyne receiver that sequentially scans the frequency spectrum channel-by-channel and measures the signal power [18]. The modern digital methods are based on an application of the discrete Fourier transform (DFT) or computationally optimized fast Fourier transform (FFT) [19][20][21]. The common method for a consequent spectrum estimation based on the Fourier transform is a short-time Fourier transform (STFT) [22,23].
There are several possible ways how to implement such spectrum computation algorithms. Therefore, the critical design analysis has to be made to choose an optimal design not only to fulfill the system requirements but also to respect implementation issues and realization costs.
The observations and records from the modernized spectrograph will be used both for a study of solar processes and for an investigation of the impact of a solar radio emission onto the current and future communication and radio navigation systems in the frame of several research projects.

Technical Requirements
The technical requirements on the solar radio spectrograph were defined on the base of the long-term experiments with the current spectrograph with respect to the parameters of the similar instruments [17,18,24,25].
The main attention during the conceptual design was paid to a high sensitivity and a wide dynamic range. The secondary design parameters superseding the parameters of the current analog instrument [15] were: • Processing frequency range min. 1/2 GHz. • Frequency bins spacing <1 MHz.
The high sensitivity requirement means that the noise floor of the telescope is below the quiet Sun radiation (QSR) level to be able to recognize its variation. The spectral flux density of the quiet Sun varies from 50 to 300 SFU while the maximal observed solar burst flux is up to 1,000,000 SFU [11].
The classical approach of the radio spectrograph design is based on a super heterodyne receiver that sequentially scans individual frequency bins. The number N of the investigated bins can be expressed as: where f max − f min is a frequency range and ∆ f is a frequency spectra bins spacing. The processing time t p of an individual bin is: where t r is a time resolution and t t is a time for returning to the next frequency bin. Let us note, that t t consists of a time for local oscillator retuning and a time for the removing of the transmission of the receiver filters. The time of measurement (i.e., the averaging in each frequency bin) is very short and thus the sensitivity (noise suppression) is poor. A modern approach oriented on the digital processing is based on the digitalization of the whole signal in early stages of the reception chain and frequency spectrum is then estimated by FFT. The results are further smoothed (a Bartlett's approach [26]) with the aim of reduction of the power spectra estimation variance. The averaged samples of the power spectral density C k are expressed by: where C k,m is a discrete power spectral density calculated from the m-th vector of input samples and M is a number of averaged consecutive vectors of the signal samples. The variance of the frequency spectra estimation is reduced M times.
Since M ≥ N and number of frequency bins is several hundreds or thousands, the typical improvement of SNR in the modern "parallel" approach is about 20 dB or more in comparison to the "serial scanning" approach in the past.

Receiver Hardware Concept
There are several existing projects of analog scanning and digital FFT based solar spectrographs [18,24,25,27,28], where usually prevails an approach of the reception, detection, and sensing of signal in the analog domain with a final conversion of the detector output to digital data for the storage and viewing. On the other hand, a modern digital spectrograph is de facto a wideband receiver equipped with the FFT processor for the signal spectrum calculation from samples produced by an analog-to-digital converter (ADC) working in baseband or intermediate-frequency band output of receiver. As the spectrograph is a unique device, there are no usable application specific integrated circuits (ASICs) for its realization (apart from the extremely high costs for their development), so the digital processor is usually programmed to a field-programmable gate array (FPGA) or a digital signal processor (DSP) [24]. This architecture is often called a software radio or a software defined radio (SDR).
The preliminary design study identifies several critical problems: 1.
Need of the high-speed analog-to-digital converters; 2.
Cost of the high-performance FPGA; 3. Software IP cores for FPGA for interfacing the ADC; 4.
Challenging design of high selectivity wideband radio frequency (RF) filters or intermediate frequency (IF) filters.
The first design concept assumed to use a high-speed ADC with the several gigahertz sampling rate to be able to process a full 1 GHz bandwidth in a single channel (Figure 2a). Such a choice would however lead to a very high cost of such ADC and a necessity to use a commercial proprietary IP (intellectual property) core for the interface between the ADC and the FPGA through JESD204B standard [29] that was also very expensive. The further cost analyses had shown that if the processing bandwidth would be divided into several sub-bands (a multi-channel approach, Figure 2b), the much cheaper ADCs with LVDS (low-voltage differential signaling) interface that does not need an expensive IP core can be used and thus the final cost of the instrument would be then significantly reduced. Let us note, that there is no need to process the working bandwidth in a single channel neither from signal processing nor scientific point of view. The cost saving brought by the multi-channel solution is so high that the total hardware cost of four channels is lower than the price of JESD204B IP core for the FPGA. Therefore, we decided to use the multi-channel concept in our modernization project. The first design concept assumed to use a high-speed ADC with the several gigahertz sampling rate to be able to process a full 1 GHz bandwidth in a single channel (Figure 2a). Such a choice would however lead to a very high cost of such ADC and a necessity to use a commercial proprietary IP (intellectual property) core for the interface between the ADC and the FPGA through JESD204B standard [29] that was also very expensive. The further cost analyses had shown that if the processing bandwidth would be divided into several sub-bands (a multi-channel approach, Figure 2b), the much cheaper ADCs with LVDS (low-voltage differential signaling) interface that does not need an expensive IP core can be used and thus the final cost of the instrument would be then significantly reduced. Let us note, that there is no need to process the working bandwidth in a single channel neither from signal processing nor scientific point of view. The cost saving brought by the multi-channel solution is so high that the total hardware cost of four channels is lower than the price of JESD204B IP core for the FPGA. Therefore, we decided to use the multi-channel concept in our modernization project. The problem with a design and realization of high selectivity RF or IF filters, required in the front end for a suppression of the overlapping spectral components in case of the heterodyne concept, can be solved by a selection of the direct conversion receiver ( Figure 3) that requires high selectivity low-pass filters instead of the band-pass ones. Let us note, that the low-pass filter with a bandwidth of several hundred megahertz can be realized as a lumped elements filter from standard SMD components while the band-pass IF filter of a super heterodyne receiver or the RF filter of a tuned radio receiver front end must be realized as a customized product based on the expensive ceramic resonator filter technology or on other special filter technology. The selectivity of the RF bandpass filter at the input of front end is not critical because it shall only additionally suppress the products outside the processed band 1-2 GHz and its realization is common for all channels. The problem with a design and realization of high selectivity RF or IF filters, required in the front end for a suppression of the overlapping spectral components in case of the heterodyne concept, can be solved by a selection of the direct conversion receiver ( Figure 3) that requires high selectivity low-pass filters instead of the band-pass ones. Let us note, that the low-pass filter with a bandwidth of several hundred megahertz can be realized as a lumped elements filter from standard SMD components while the band-pass IF filter of a super heterodyne receiver or the RF filter of a tuned radio receiver front end must be realized as a customized product based on the expensive ceramic resonator filter technology or on other special filter technology. The selectivity of the RF bandpass filter at the input of front end is not critical because it shall only additionally suppress the products outside the processed band 1-2 GHz and its realization is common for all channels. On the other hand, a direct conversion receiver [30][31][32][33] suffers from the following problems: 1. Problem with the semiconductor noise 1/f. 2. DC offsets, and an interference of the local oscillator that oscillates in the middle of the processing frequency band. 3. Amplitude and phase unbalance that causes a crosstalk between upper and lower sidebands.
Fortunately, the first problem can be solved by the careful design of the level plan. The signal must be sufficiently amplified before the conversion to the baseband to avoid an influence of the 1/f noise. In our case, the power consumption of the receiver front end is not critical so using a wide dynamic range quadrature mixer with higher power consumption enables successful mixing of the high-level signal from the amplified input stage.
The problem with a DC offset can be solved by compensation feedback circuits, however paid by the distortion of the central frequency on which oscillates a local oscillator. The most radio systems count with this fact and use modulation schemas that are resistive to this distortion.
The solution to this problem in our design is to use a higher frequency resolution of the signal. The distorted frequency bin is ignored and replaced by an interpolated value.
The crosstalk between upper and lower sidebands caused by the receiver imperfection or insufficient suppression is a problem that can completely spoil the measurement or can bring unwanted effects in the spectrum. The crosstalk suppression is expressed as: where and are the amplitude and phase unbalance, respectively. A typical amplitude and phase unbalance of an uncompensated quadrature direct conversion mixer is about 0.5 dB and 1°. The crosstalk suppression is then about 30 dB.
The higher suppression can be achieved by the mixer calibration typically with the help of the feedback circuits. The receiver front end must be equipped with a possibility to control the DC offset and gain of the baseband I and Q branches and adjust the phase shift of the local oscillator signal. The properly designed compensation can increase the cross talk suppression up to 60 dB that is enough for the most applications.
The outputs of DC offset, amplitude, and phase imbalance detectors can be expressed: , On the other hand, a direct conversion receiver [30][31][32][33] suffers from the following problems: 1.
Problem with the semiconductor noise 1/f.

2.
DC offsets, and an interference of the local oscillator that oscillates in the middle of the processing frequency band.

3.
Amplitude and phase unbalance that causes a crosstalk between upper and lower sidebands.
Fortunately, the first problem can be solved by the careful design of the level plan. The signal must be sufficiently amplified before the conversion to the baseband to avoid an influence of the 1/f noise. In our case, the power consumption of the receiver front end is not critical so using a wide dynamic range quadrature mixer with higher power consumption enables successful mixing of the high-level signal from the amplified input stage.
The problem with a DC offset can be solved by compensation feedback circuits, however paid by the distortion of the central frequency on which oscillates a local oscillator. The most radio systems count with this fact and use modulation schemas that are resistive to this distortion.
The solution to this problem in our design is to use a higher frequency resolution of the signal. The distorted frequency bin is ignored and replaced by an interpolated value.
The crosstalk between upper and lower sidebands caused by the receiver imperfection or insufficient suppression is a problem that can completely spoil the measurement or can bring unwanted effects in the spectrum. The crosstalk suppression is expressed as: where A U and θ U are the amplitude and phase unbalance, respectively. A typical amplitude and phase unbalance of an uncompensated quadrature direct conversion mixer is about 0.5 dB and 1 • . The crosstalk suppression is then about 30 dB. The higher suppression can be achieved by the mixer calibration typically with the help of the feedback circuits. The receiver front end must be equipped with a possibility to control the DC offset and gain of the baseband I and Q branches and adjust the phase shift of the local oscillator signal. The properly designed compensation can increase the cross talk suppression up to 60 dB that is enough for the most applications.
The outputs of DC offset, amplitude, and phase imbalance detectors can be expressed: d p = n I n Q n n I 2 n + n Q 2 n .
The values, produced by these detectors, are then used for a compensation of negative effects described above by controlling the oscillator block and variable gain amplifiers (VGA) in both I and Q branches.
The outputs of I and Q branches of the front-end receiver are then passed to dual-channel ADC and the samples are processed in a FPGA realizing the FFT computation of the full-bandwidth spectrum. The computed samples of the spectrum are passed through a standard interface to a connected computer or server for further processing, viewing, and storage. The FPGA can be also programmed to realize a detection of the receiver imperfections described above and provide the results to the control and setup processor for the compensation purposes.
To be able to simply combine measurements of several spectrograph channels, the measurements must be synchronized with 1 PPS (pulse per second) pulses. The source of such synchronization signal can be for example a 1 PPS output provided by a GPS receiver or an arbitrary frequency standard. If the pulses are not present, the spectrograph channels are switched to the free running mode. The measurement cycle is prolonged on to 1.1 s. The complete block scheme of one channel is at Figure 4. The values, produced by these detectors, are then used for a compensation of negative effects described above by controlling the oscillator block and variable gain amplifiers (VGA) in both I and Q branches.
The outputs of I and Q branches of the front-end receiver are then passed to dual-channel ADC and the samples are processed in a FPGA realizing the FFT computation of the full-bandwidth spectrum. The computed samples of the spectrum are passed through a standard interface to a connected computer or server for further processing, viewing, and storage. The FPGA can be also programmed to realize a detection of the receiver imperfections described above and provide the results to the control and setup processor for the compensation purposes.
To be able to simply combine measurements of several spectrograph channels, the measurements must be synchronized with 1 PPS (pulse per second) pulses. The source of such synchronization signal can be for example a 1 PPS output provided by a GPS receiver or an arbitrary frequency standard. If the pulses are not present, the spectrograph channels are switched to the free running mode. The measurement cycle is prolonged on to 1.1 s. The complete block scheme of one channel is at Figure 4. The communication between blocks and the output interface can be based on common standards (SPI, I2C, USB) according to parameters of chosen hardware components. The general control and setup are then usually provided by some kind of microcontroller unit (MCU).

Spectrograph Hardware
The described hardware concept was applied for a four-channel spectrograph in a frequency range of 1-2 GHz ( Figure 5) using the same RF band-pass filters and flat gain amplifiers both for the common RF input and for the front-end part of the channels. The RF filter is an of-the-shelf LTCC (low temperature cofired ceramics) filter (BFCG-162W+) with the pass band 950-2200 MHz, pass band insertion loss 1.8 dB, and the stop band insertion loss better than 25 dB. The flat gain amplifier is based on an off-the-shelf 40-4000 MHz RF gain block (TRF37C75) with 18 dB gain. These common blocks are placed together with an RF channel splitter, power supply modules, and a GPS module (for 1 PPS signal provision) to separate the board. Each spectrograph channel is then constructed as an independent block in its own case ( Figure 6). The hardware parameters are summarized in Table 1. The communication between blocks and the output interface can be based on common standards (SPI, I2C, USB) according to parameters of chosen hardware components. The general control and setup are then usually provided by some kind of microcontroller unit (MCU).

Spectrograph Hardware
The described hardware concept was applied for a four-channel spectrograph in a frequency range of 1-2 GHz ( Figure 5) using the same RF band-pass filters and flat gain amplifiers both for the common RF input and for the front-end part of the channels. The RF filter is an of-the-shelf LTCC (low temperature cofired ceramics) filter (BFCG-162W+) with the pass band 950-2200 MHz, pass band insertion loss 1.8 dB, and the stop band insertion loss better than 25 dB. The flat gain amplifier is based on an off-the-shelf 40-4000 MHz RF gain block (TRF37C75) with 18 dB gain. These common blocks are placed together with an RF channel splitter, power supply modules, and a GPS module (for 1 PPS signal provision) to separate the board. Each spectrograph channel is then constructed as an independent block in its own case ( Figure 6). The hardware parameters are summarized in Table 1.

FPGA Signal Processor Content
The receiver is equipped with a medium rate FPGA Artix7 by Xilinx (San Jose, CA, USA). The block diagram of the FPGA configuration in Figure 7 demonstrates a sequence of processing

FPGA Signal Processor Content
The receiver is equipped with a medium rate FPGA Artix7 by Xilinx (San Jose, CA, USA). The block diagram of the FPGA configuration in Figure 7 demonstrates a sequence of processing

FPGA Signal Processor Content
The receiver is equipped with a medium rate FPGA Artix7 by Xilinx (San Jose, CA, USA). The block diagram of the FPGA configuration in Figure 7 demonstrates a sequence of processing operations.
In order to meet the timing constraints, one second of a sampled signal data stream is separated in the splitter block to two consecutive half sized sub-blocks that are put into two parallel branches. The FFT operation with length 1024 is therefore calculated only with a half of processed data block and can work on a half of the sampling frequency, i.e., 155 MHz (details in Table 2). operations. In order to meet the timing constraints, one second of a sampled signal data stream is separated in the splitter block to two consecutive half sized sub-blocks that are put into two parallel branches. The FFT operation with length 1024 is therefore calculated only with a half of processed data block and can work on a half of the sampling frequency, i.e., 155 MHz (details in Table 2).  The next block executes a signal spectrum averaging of 300 FFT measurements with the averaging time of 1 ms. This block is synchronized to the 1 PPS signal to be able to synchronize the spectrograph measurement onto an external time signal.
The measurement is then sent via the USB 2.0 interface to the server for storage and further processing. The data management software uses a standard FITS (flexible image transport system) format for data files [34].
The FPGA also realizes the detectors for a determination of the DC offset, amplitude imbalance, and phase imbalance. The compensation loops are closed through an adjacent ARM processor that controls through a SPI interface DC offsets, amplitude, and phase of the quadrature mixer. The other function of the ARM processor is a receiver general control and a setup of the front end and ADC.

Testing
During the spectrograph development and realization, the device was tested to verify the design parameters fulfillment. The test operation had started in the summer of 2018 at the Ondrejov astronomical observatory (Figure 8) with the aim of verification of the spectrograph performance and a comparison of its measurements with the current analog device. The developed spectrograph was tested with the backup 8 m dish antenna originated from the Second World War German radar while the current analog spectrograph is installed at the modern 9.5 m dish antenna that profits from the higher gain. The reason for such installation was to be able to test the new design parameters and operation setup without interruption of the current astronomical experiments. Due to a preliminary character of the tests, only one channel tuned to the center frequency 1525 MHz was put in the operation. The test frequency band was chosen to cover the common satellite navigation systems' frequencies to be able to investigate the solar burst impacts on these systems as stated in the previous part of the paper.  The next block executes a signal spectrum averaging of 300 FFT measurements with the averaging time of 1 ms. This block is synchronized to the 1 PPS signal to be able to synchronize the spectrograph measurement onto an external time signal.
The measurement is then sent via the USB 2.0 interface to the server for storage and further processing. The data management software uses a standard FITS (flexible image transport system) format for data files [34].
The FPGA also realizes the detectors for a determination of the DC offset, amplitude imbalance, and phase imbalance. The compensation loops are closed through an adjacent ARM processor that controls through a SPI interface DC offsets, amplitude, and phase of the quadrature mixer. The other function of the ARM processor is a receiver general control and a setup of the front end and ADC.

Testing
During the spectrograph development and realization, the device was tested to verify the design parameters fulfillment. The test operation had started in the summer of 2018 at the Ondrejov astronomical observatory (Figure 8) with the aim of verification of the spectrograph performance and a comparison of its measurements with the current analog device. The developed spectrograph was tested with the backup 8 m dish antenna originated from the Second World War German radar while the current analog spectrograph is installed at the modern 9.5 m dish antenna that profits from the higher gain. The reason for such installation was to be able to test the new design parameters and operation setup without interruption of the current astronomical experiments. Due to a preliminary character of the tests, only one channel tuned to the center frequency 1525 MHz was put in the operation. The test frequency band was chosen to cover the common satellite navigation systems' frequencies to be able to investigate the solar burst impacts on these systems as stated in the previous part of the paper.

Test Results
During the early one-year test period only one observable event at 20 March 2019 was registered. The reason is that the Sun activity is at its long-term minimum. The spectrum registered by the current spectrograph is on Figure 9, while the measurements of the new digital spectrograph are on Figures 10-13. All spectrograms were depicted relative to reference time (tr = 0 s) 11:11 UTC. Figures 10 and 11 represent the same spectrum part as a spectrogram to respect the style of Figure 9 and as a 3D spectrum waterfall plot. There was chosen an interval of 40 s that symmetrically covers the reference time for a demonstration of a modernized spectrograph data in full resolution. The datafile of the related interval represents a file of size 200 MB, therefore the long-time records are stored as one-second averages ( Figure 13). They are also used for a coarse detection of solar events by an expert observer. The reference value of quiet Sun for all graphs is 100% (shown as 20 dB of relative power Pr for new spectrograph images and as a ratio 1 for the old analog spectrograph). The drops of the spectrum on the sides of the evaluated band were caused by channel low-pass filter characteristics. This drop would be compensated by software processing during combination of overlapping adjacent channels.

Test Results
During the early one-year test period only one observable event at 20 March 2019 was registered. The reason is that the Sun activity is at its long-term minimum. The spectrum registered by the current spectrograph is on Figure 9, while the measurements of the new digital spectrograph are on Figures 10-13. All spectrograms were depicted relative to reference time (t r = 0 s) 11:11 UTC. Figures 10  and 11 represent the same spectrum part as a spectrogram to respect the style of Figure 9 and as a 3D spectrum waterfall plot. There was chosen an interval of 40 s that symmetrically covers the reference time for a demonstration of a modernized spectrograph data in full resolution. The datafile of the related interval represents a file of size 200 MB, therefore the long-time records are stored as one-second averages ( Figure 13). They are also used for a coarse detection of solar events by an expert observer. The reference value of quiet Sun for all graphs is 100% (shown as 20 dB of relative power P r for new spectrograph images and as a ratio 1 for the old analog spectrograph). The drops of the spectrum on the sides of the evaluated band were caused by channel low-pass filter characteristics. This drop would be compensated by software processing during combination of overlapping adjacent channels.

Test Results
During the early one-year test period only one observable event at 20 March 2019 was registered. The reason is that the Sun activity is at its long-term minimum. The spectrum registered by the current spectrograph is on Figure 9, while the measurements of the new digital spectrograph are on Figures 10-13. All spectrograms were depicted relative to reference time (tr = 0 s) 11:11 UTC. Figures 10 and 11 represent the same spectrum part as a spectrogram to respect the style of Figure 9 and as a 3D spectrum waterfall plot. There was chosen an interval of 40 s that symmetrically covers the reference time for a demonstration of a modernized spectrograph data in full resolution. The datafile of the related interval represents a file of size 200 MB, therefore the long-time records are stored as one-second averages ( Figure 13). They are also used for a coarse detection of solar events by an expert observer. The reference value of quiet Sun for all graphs is 100% (shown as 20 dB of relative power Pr for new spectrograph images and as a ratio 1 for the old analog spectrograph). The drops of the spectrum on the sides of the evaluated band were caused by channel low-pass filter characteristics. This drop would be compensated by software processing during combination of overlapping adjacent channels.         To consider the system initial parameters, the measurements for a quiet Sun were executed as well ( Figure 14) and compared to the values obtained by a system with the antenna pointed to the clear sky and during the night to determine a QSR-to-noise system ratio.   To consider the system initial parameters, the measurements for a quiet Sun were executed as well ( Figure 14) and compared to the values obtained by a system with the antenna pointed to the clear sky and during the night to determine a QSR-to-noise system ratio. To consider the system initial parameters, the measurements for a quiet Sun were executed as well ( Figure 14) and compared to the values obtained by a system with the antenna pointed to the clear sky and during the night to determine a QSR-to-noise system ratio. The observed difference of the signal level of about 3 dB for quiet Sun vs. clear sky does not sound very optimistic, but we have to take into account the non-optimal state of the used old antenna in comparison to the final installation. The most significant improvement can be achieved by the placement of a low noise amplifier directly to the antenna feed output as is in the case of the RT5 antenna. Nevertheless, we were still be able to observe quantitative changes of the spectrum level in comparison with the quiet Sun, therefore we can expect even better results for the full operation state. The expected values in the future are about 10% of the quiet Sun level, closer to the pure QSR to noise ratio obtained from the comparison of a QSR measurement with a night operation, where the ratios were 10 dB for single bin spikes of the spectrum and 15 dB for the average value of noise.
The frequency distribution of the noise level in the realized one channel band is mainly determined by the used low-pass filter characteristics as can be seen in 3D plots. Besides this impact, there were observed two spurs on frequency bins 1.436 GHz and 1.615 GHz caused by internal circuits of the receiver. These single-bin spurs will be excluded from the final data during the signal processing and substituted by an interpolated value from adjacent bins.
The testing measurements were also used to identify main sources of radio interference. The antenna placement is on the south-western slope of the hill in an elevation of about 510 m. The nearby Ondrejov village is about 1 km from the observatory and with an elevation of about 450 m, which helps to slightly reduce the impact of a possible close industrial interference. Main distant sources of the radio interference (considering the final full frequency range 1-2 GHz) can be a cellular phone network transmitter or a nearby mobile phone (850 MHz, 950 MHz, and 1800 MHz band and higher harmonics of 850 MHz band) and a digital terrestrial TV transmitter (higher harmonics of 562-706 and 714-858 MHz bands). The impact of these interference sources can be suppressed by preventing the antenna pointing directly to the transmitter locations. The other source of the interference signal can be an ADS-B (automatic dependent surveillance-broadcast) signal at the frequency 1090 MHz. The signal is a product of an aircraft surveillance system transmitted in regular intervals during the flight and it can pose the notable source of interference because the location of a transmitter can be directly in the antenna beam due to the aircraft movement across the sky. Fortunately, such interference will be active only for a short time interval.
The notable interference was during the test operation observed at frequencies 1.458 GHz with relative power 36 dB and 1.467 GHz with relative power 27 dB. Although the source of the interference was not identified, the observed levels are within designed dynamic range of the system and therefore do not cause a saturation of ADC. They will be digitally filtered in the processed signal The observed difference of the signal level of about 3 dB for quiet Sun vs. clear sky does not sound very optimistic, but we have to take into account the non-optimal state of the used old antenna in comparison to the final installation. The most significant improvement can be achieved by the placement of a low noise amplifier directly to the antenna feed output as is in the case of the RT5 antenna. Nevertheless, we were still be able to observe quantitative changes of the spectrum level in comparison with the quiet Sun, therefore we can expect even better results for the full operation state. The expected values in the future are about 10% of the quiet Sun level, closer to the pure QSR to noise ratio obtained from the comparison of a QSR measurement with a night operation, where the ratios were 10 dB for single bin spikes of the spectrum and 15 dB for the average value of noise.
The frequency distribution of the noise level in the realized one channel band is mainly determined by the used low-pass filter characteristics as can be seen in 3D plots. Besides this impact, there were observed two spurs on frequency bins 1.436 GHz and 1.615 GHz caused by internal circuits of the receiver. These single-bin spurs will be excluded from the final data during the signal processing and substituted by an interpolated value from adjacent bins.
The testing measurements were also used to identify main sources of radio interference. The antenna placement is on the south-western slope of the hill in an elevation of about 510 m. The nearby Ondrejov village is about 1 km from the observatory and with an elevation of about 450 m, which helps to slightly reduce the impact of a possible close industrial interference. Main distant sources of the radio interference (considering the final full frequency range 1-2 GHz) can be a cellular phone network transmitter or a nearby mobile phone (850 MHz, 950 MHz, and 1800 MHz band and higher harmonics of 850 MHz band) and a digital terrestrial TV transmitter (higher harmonics of 562-706 and 714-858 MHz bands). The impact of these interference sources can be suppressed by preventing the antenna pointing directly to the transmitter locations. The other source of the interference signal can be an ADS-B (automatic dependent surveillance-broadcast) signal at the frequency 1090 MHz. The signal is a product of an aircraft surveillance system transmitted in regular intervals during the flight and it can pose the notable source of interference because the location of a transmitter can be directly in the antenna beam due to the aircraft movement across the sky. Fortunately, such interference will be active only for a short time interval.
The notable interference was during the test operation observed at frequencies 1.458 GHz with relative power 36 dB and 1.467 GHz with relative power 27 dB. Although the source of the interference was not identified, the observed levels are within designed dynamic range of the system and therefore do not cause a saturation of ADC. They will be digitally filtered in the processed signal data. Let us note that the RT5 spectrograph is more covered by the hill slope from interference sources and we can expect even better results in the future. Of course, the next interference analysis has to be made for the final installation to the location of current analog system.

Discussion and Conclusions
Even though the developed spectrograph used a lower gain old antenna during the tests, the obtained measurements carry out more details and a comparable background noise with respect to the analog spectrograph. The comparison of the parameters of both spectrographs is presented in the Table 3. The measurement of the current spectrograph is saved in a proprietary format while the measurement of the developed spectrograph is archived in the FITS format that allows further data analyses and easy sharing and manipulation. The indisputable advantage of the new spectrograph is a precise amplitude scale in contrast to the current one, which has its amplitude scale distorted because of the imperfect analog processing. Moreover, the higher time and frequency resolution of the new device enables a better identification and analysis of possible interference sources, such as microwave data links, mobile phone systems signals, air traffic radio navigation signals (e.g., ADS-B), terrestrial broadcast of digital television signals, etc.
The test of the L band solar spectrograph based on the modern digital signal processing at Ondrejov astronomical observatory confirmed the expected higher time and frequency resolution, and a high linearity of the amplitude measurement. The low processing noise that can be achieved by the massive averaging allowed by a parallel design (in contrast to the spectrum scanning approach of the current analog device) was however affected by the lower performance of the used antenna system and therefore the observed noise parameters persevered to stay at the original values. This problem is expected to be suppressed using the final installation antenna with better parameters. The crucial point was a proof of a validity of such design, which was accomplished according to the expectations.
As the test operation was evaluated as successful, the next step, expected in the end of 2019, will be an installation of the complete multi-channel system to the modern 10 m dish antenna.
The reasonable price of the spectrograph (approximately 2 kEuro per 250 MHz channel) was achieved by splitting the operating frequency range to several synchronized 250 MHz channels. This approach allows the use of moderately priced components and manufacture technologies without any negative impact on the developed instrument performance.
Author Contributions: Initiation of the spectrograph development and preliminary parameter specification, M.B.; development of the technical concept, P.K.; application of the spectrograph for investigation of the solar radio burst effect on the navigation and communication systems, P.K.; hardware development, FPGA programming and laboratory testing, P.K.; software for signal processing and the measurement format conversions and storage, P.P.; spectrograph installation P.P., P.K., M.B.; spectrograph operation, M.B.; preliminary measurement analyses, P.P.; writing-original draft, P.K.; writing-review and editing, P.P.; project administration P.K.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.