Zuccheretti 647_659 647 ANNALS OF GEOPHYSICS, VOL. 46, N. 4, August 2003 The new AIS-INGV digital ionosonde Enrico Zuccheretti (1), Giuseppe Tutone (1), Umberto Sciacca (1), Cesidio Bianchi (1) and Baskaradas James Arokiasamy (2) (1) Istituto Nazionale di Geofi sica e Vulcanologia, Roma, Italy (2) TRIL fellow, The Abdus Salam International Center for Theoretical Physics, Trieste, Italy Abstract A new digital ionosonde called AIS-INGV (Advanced Ionospheric Sounder) was designed both for research and for routine service of HF radio wave propagation forecast. Nearly the entire system was developed in the Laboratorio di Geofi sica Ambientale at the Istituto Nazionale di Geofi sica e Vulcanologia (INGV), Rome. It exploits advanced techniques for signal analysis, recent technological devices and PC resources. This paper describes design concepts and performance of the new ionosonde. Mailing address: Dr. Enrico Zuccheretti, Istituto Nazionale di Geofi sica e Vulcanologia, Via di Vigna Murata 605, 00143 Roma, Italy; e-mail: zuccheretti@ingv.it Key words ionosonde – pulse compression – phase coherent integration 1. Introduction Since their origin, radio and radar techniques have been used to sound the ionosphere. They were suitable both to study the higher atmosphere and to forecast the propagation conditions useful for applications on radio communications. From 1970 onwards, modern systems have been added to the older ones: making use of analogue techniques they were able only to measure the electron density and the virtual height of the refl ecting layers. The more recent and modern ionosondes are called Advanced Ionospheric Sounders (AIS). As the RF de- signers started using techniques such as pulse compression, coherent integration etc., the use of digital AIS reduces the transmitted power maintaining a favourable S/N ratio. The new sounders perform analysis on the received echo which enables the extraction of more parameters (such as polarisation conditions, electron density and doppler shift) from the echo. A more detailed review of the ionosondes evolution can be found in Davies (1990), Hunsuker (1991) and Bibl (1998). Based on the above mentioned characteristics, a digital low power pulse compressed ionosonde was developed in the Laboratorio di Geofi sica Ambientale (LGA) of the Istituto Nazionale di Geofi sica e Vulcanologia (INGV), Rome. The objective of this ionosonde called AIS-INGV was to reduce the transmitted peak power (less than 250 W) and consequently its weight and size, power consumption and hardware complexity. The AIS-INGV ionosonde is completely programmable. The control, data acquisition, on-line processing and storage of the acquired data are supported by a directly interfaced PC. This ionosonde employs a 16 bit complementary phase code and exploits the most advanced HF- radar techniques, such as pulse compression and phase coherent integration. The process- ing gain due to correlation process and phase coherent integration compensates the reduction of the transmitted power. The fi rst prototype is capable of measuring the virtual height. This 648 Enrico Zuccheretti, Giuseppe Tutone, Umberto Sciacca, Cesidio Bianchi and Baskaradas James Arokiasamy basic system allows future expansions like polarization information, doppler analysis, etc. This paper describes design concepts, digital signal processing and system performances of the new ionosonde AIS-INGV. 2. System design concepts 2.1. Specifi cations The specifi cations of the new system are presented in table I. The new ionosonde has been designed and built on the base of the theory of radar systems, applied to the study of ionosphere (Skolnik, 1980, 1997; Hunsucker, 1991). Figures 2 and 3, together with table II, show the confi guration of the ionosonde with arrangement of its subsystems and boards. The complete design report by Arokiasamy et al. (2002), gives a detailed description of the structure and operation of the various constitutive boards. Since a detailed analysis is beyond the scope of this paper, only a synthetic description of the ionosonde referring to fi g. 1 (which is a simplifi ed functional diagram of the whole sys- tem) will be presented here. 2.2. General description In fi g. 1, the ionosonde has been divided into functional blocks representing the main functions; one block does not necessarily correspond to physical blocks or circuits. The thicker lines refer to digital buses. The blocks are grouped into two sections corresponding approximately to the transmitting and the receiving sections: on the left, there is a power amplifi er (PWA), a frequency synthesiser (SYN) and the code generator (CODE GEN); on the right there is the radio receiver (Rx), the Analog to Digital Converter (ADC) and the Digital Signal Processor (DSP). Every block is controlled by a personal computer (PC control and storage), which can store and display data. The antennas (Tx ant. and Rx ant.) complete the system. In order to increase the signal-to-noise ratio without increasing the transmitted peak power, the pulse compression technique (encoded pulse) is applied. After the echo is received, it is compared to the transmitted code by means of a correlation process. To this purpose a digital bi- phase complementary code is used. It lets the RF carrier assume only two values for its phase: 0° Table I. AIS-INGV ionosonde specifi cations. The symbol (*) indicates the maximum range of input received powers, in all possible conditions; it is not required that the system maintains the linearity over it, provided that it allows the recognition of the echo position. Parameter Requirement Height range 90 ∏ 750 km Distnace resolution 5 km Max. peak Tx power (average power) 250 W (5 ~ 10 W) Receiver sensitivity ~ - 85 dBm for 0 dB S/N Input dynamic range ~ 80 dB (*) Frequency range 1 ∏ 20 MHz Frequency resolution (step) 25, 50, 100 kHz Frequency scan duration (max.) 3 min (with 50 kHz step) Acquisition sampling rate ~ 100 kHz Acquisition quantization 8 bit Storage data rate (max.) 60 kbytes (with 50 kHz step) 649 The new AIS-INGV digital ionosode and 180° with respect to a reference, following a pattern designed to achieve the cancellation of the correlation «sidelobes». The duration of a single phase of the code (called «sub-pulse») was set to 30 ms after the design considerations based on the specifi cations. 2.3. Frequency synthesis and code generation The AIS-INGV adopts a receiving heterodyne system based on multiple conversions. This kind of system allows the reduction of internally generated noise due to the «image» interference. In order to make the system effective, it is im- portant to raise the Intermediate Frequency (IF) to move the images far enough from the input passband. The IFs have been chosen at 35.9, 4.1 and 0.1 MHz. In this way, the image related to the fi rst conversion (from 1-20 MHz to 35.9 MHz) lies between 72.8 and 91.8 MHz, well beyond the RF input band. Similarly the second conversion (from 35.9 to 4.1 MHz) has its image at 44.1 MHz, the third (from 4100 to 100 kHz) at 3.9 MHz; all the images are easily removable by means of a fi lter before the conversion. To accomplish such conversions, three Local Oscillators (LO) at 36.9-55.9, 40 and 4 MHz were devised (note that only the fi rst has variable frequency). Since digital techniques that implement a direct synthesis of a sinusoid with a determined frequency and phase were available, it was decided to use such «DDS» devices (Direct Digital Synthesis). Information from the received signal is ex- tracted by means of the «quadrature» demodu- lation. It demands maximum phase coherence between all signals generated in the system. So, though different DDS devices are used, their reference is unique: an OX 125 MHz quartz oscillator. All the frequencies and timing signals are derived from this unique reference oscillator. So the whole system is phase locked and a drift (few Hertz over some MHz) of the reference oscillator does not affect the coherence of the system. All the DDS devices were mounted on a single board (SYN): two of them generate the fi rst two local oscillator outputs; the third Fig. 1. Simplifi ed functional diagram. 650 Enrico Zuccheretti, Giuseppe Tutone, Umberto Sciacca, Cesidio Bianchi and Baskaradas James Arokiasamy generates the Tx carrier. The third local oscillator output at 4 MHz is derived from the second at 40 MHz by means of frequency division by ten. All the LO outputs are fi ltered to clean the waveforms, to make them as pure as possible. All DDS devices can output a specifi ed frequency value by means of proper programming, coming from the PC bus. A supplementary 400 kHz square wave is generated from the 4 MHz; it acts as a reference for other sub-systems (e.g., the ADC), that must all be phase-locked with respect to the 125 MHz reference. The reference clock (400 kHz in fi g. 1) is used to create the codes (30 ms period corresponding to the sub-pulse length). This function was implemented in a different board, the code generator (CODE GEN). The 400 kHz clock (2.5 ms period) is also used to generate the quadrature sampling clock. The codes are implemented as digital waveforms, with the 1 s and the 0 s corresponding to the carrier phase rotations. While the proper timing is derived inside the board, the PC has only to send the exact set of 1 s and 0 s. In this way, the system is very fl exible; it is possible to easily change the type and length of the code. ) ) ( p , ) Fig. 2. The ionosonde subsystems. Interconnections (antennas not included). 651 The new AIS-INGV digital ionosode T x O sc il la to r ( ) gg ( ) ( ) ( ) ( ) g F ig . 3. F un ct io na l di ag ra m o f th e m ai n un it o f A IS -I N G V i on os on de . 652 Enrico Zuccheretti, Giuseppe Tutone, Umberto Sciacca, Cesidio Bianchi and Baskaradas James Arokiasamy A A A A A A A A A A Once generated, the codes (CODE in fi g. 1) are sent to the synthesiser board to modulate the RF carrier. The modulated or «coded» RF carrier (CRF in fi g. 1) is then sent to the amplifi er. The modulation of the RF carrier by the code could be suffi cient to avoid the amplifi er energy output out of the transmission phase. Anyway, in order to interdict the amplifi er completely, a trigger pulse is generated inside the board (AMP trig in fi g. 1). It is obtained starting from the same reference used for the code, considering that a complete Tx pulse is a multiple of 16 sub-pulses. Table II. Symbols used in the block diagram. 653 The new AIS-INGV digital ionosode 2.4. Power amplifi cation and antenna system In fi g. 2 the transmitting power amplifi er (PWA) is considered as a separate subsystem. To increase the effi ciency and reliability of the system, an industrial product was bought outside INGV. The peak power in linear conditions was set at 250 W. The linearity specifi cations are not critical, because of their low impact on the useful RF band. On the other hand, the harmonics and spurious signals are to be maintained low in order to increase the power effi ciency and to limit EMI. The antenna system employed is a couple of crossed «delta» antennas, to minimise the direct coupling between the Tx and Rx (fi g. 2). Each antenna has a 42 m base and is mounted on a single 24 m mast. Such an antenna type is similar to the rhombic antenna (to which it is equivalent considering the image effect due to the ground) and its impedance remains almost constant over a wide range allowing good matching with a gain (1.5-3 dB). A couple of baluns are used to match the balanced antennas with the unbalanced coaxial transmission lines. 2.5. Receiver and A/D conversion The demodulation is performed by means of a multiple conversion heterodyne system. The three local oscillator outputs are to be used inside the receiver as LO inputs of mixers which actually perform the conversions. Before and after the conversion fi lters have to be added to clean the signal from noise. Specifi cally, a separate board was thought to be put at the Rx antenna input to select a relatively narrow RF band. In detail, a preliminary fi ltering, wide band, is put as the fi rst stage. It has a pass band in the ionosonde operative range (1-20 MHz) and some protection against high voltage bursts. Then six band pass fi lters follow; they are switchable, i.e. only one of them is actually working at a time. This combination reduces the number of RF fi lters maintaining the response relatively fl at. The fi lters have different bandwidths, the lower frequency fi lters have a narrower bandwidth with respect to the higher frequency ones (the bandwidth is around the half of the central frequency). Switches are used in the receiving chain, in order to inter- dict the reception during transmission. The output of the RF band selection fi lter is fed to the fi rst mixer, inside the real receiver (the switchable fi lters and the receiver are two boards reported in fi g. 1 as the Rx block). In the receiver chain, the three mixers are followed by fi lters and amplifi ers. As previously mentioned, these fi lters limit the noise bandwidth and avoid the image interference. The last fi lter is the narrowest, with its pass band limited by the bandwidth of the code: 66 kHz. All hardware fi lters (i.e. excluding fi ltering accomplished by means of digital signal processing) are passive, analog, with lumped elements (capacitors and inductors). Inside the receiver, a variable attenuator is added to increase the dynamic range. In the fi rst realisation of the ionosonde, this feature is used only as a way to manually calibrate the sys- tem and let it work well, without dynamically changing the attenuation (and the overall gain of the receiver) to match the system to the actual ionospheric attenuation and the actual noises. The attenuator is programmable by means of the PC bus, so a future expansion of the system capability is possible. Obviously, many amplifi ers are introduced in the receiving chain, to bring the signal to a good level, suitable to be processed by the following stage. The analog to digital conversion is performed in a different board (ADC). The clock is furnished by the previously seen frequency synthesiser. Some circuitry is necessary to accomplish the «quadrature», or I-Q demodulation. Two con- verters sample the analog receiver output at a 100 kHz sampling frequency (10 ms period), but one is delayed by one fourth of a period with respect to the other (2.5 ms). The sampled voltages are stored in two temporary memories, whose address is generated by a digital circuit, driven by the same 100 kHz clock used by the ADCs. Once the acquisition process is completed (about 5.2 ms) data are sent to the DSP using the same address generator. 2.6. PC control and DSP board The Personal Computer (PC control and storage) supervises the operation of the whole 654 Enrico Zuccheretti, Giuseppe Tutone, Umberto Sciacca, Cesidio Bianchi and Baskaradas James Arokiasamy system. A control program, written in a high level language, provides the interface with the operators, letting them decide the operating parameters of the soundings (starting and ending frequency, repetition rate, etc.) and displaying the ionogram while sounding. It also provides the storage of the acquired data on the local hard disk, using the internal timer to mark the time of the sounding. An interface between the PC bus and the ex- ternal ionosonde bus was introduced. A dedi- cated bus was designed to treat properly all con- trol and data fl owing from the PC. The control and data lines in fi g. 2 are a «buffered» bus that actually brings the control signals and some data packets (e.g., the code or the words to program the Tx frequency) to the Main Unit. A commercial Digital Signal Processing board (DSP) is employed for on-line analysis. It is programmed at a low level assembler lan- guage. The functions to be performed are better explained later. DSP output is a stream of values (an array for each frequency of sounding) re- presenting the time domain signal amplitude. The location of the amplitude peaks is related to the heights of the refl ecting layers. The PC analyses each trace, extracting the position of the peaks and displaying them on the screen. Independently of the display, the complete traces are stored on disk to be used in off-line processing (extraction of relevant parameters as critical frequencies, electron density, and so on). 2.7. Additional blocks In fi g. 1, some additional blocks were not included because of their minor conceptual relevance. Among these, there is a power supply (+ 12 V, - 12 V and + 5 V) to furnish the supply voltages to the various blocks, except the PC and the power amplifi er that are self powered. In some cases, a – 5 V is generated locally in boards. In all cases, the incoming supply lines have to be fi ltered to reduce the electromagnetic susceptibility to duct noise. An additional card was designed to calibrate the system during the assembly phase. It is the «synthetic echo board» (SYN ECHO), referring to its capability to output an echo-like signal. Its internal constitution is very similar to the code generator; the difference lies in the fact that the RF code is produced after a desired delay (user programmable) from the Tx trigger. In this way, the entire system can easily be check without transmitting signals on the air (the use of an external physical delay line would have needed a very long cable, more than 33 km to simulate a 50 km echo!). Referring to fi g. 1, in this operating mode the signals are to be connected as follows: the CODE does not come from «CODE GEN» but from the «SYN ECHO», while the CRF does not go to the power amplifi er but directly to the receiver through a proper attenuator. 3. Digital Signal Processing Much of the on-line processing is done at software level, inside the DSP board. After reading the data stored in the ADC board, a Complex Fast Fourier Transform (CFFT) is performed. From this point on every calculation is carried out in the frequency domain (Texas Instruments, 1994; Oppenheim et al., 1999). A frequency based fi ltering follows the CFFT. This fi ltering reduces the IF bandwidth a little more; a better S/N ratio is reached without affecting the code decoding. A second type of fi ltering is then performed to reduce the high amplitude signals with respect to a reference threshold. In fact, overcoming a specifi ed threshold is a sign that a narrow band interfering signal has passed through the last IF fi lter. The slight degradation of the spectrum due to the digital fi ltering is compensated by the great improvement of the Signal to Noise ratio. At this point, the fl ow branches, for the sys- tem must perform different correlation in accord- ance with the transmitted code. As previously mentioned, the code is «complementary» to elim- inate the sidelobes of the correlation process. Once the type of transmitted pulse is recognised, the correlation is performed. After the code 1 correlation, the system prepares itself for the reception of code 2. The process is repeated N times (where N is the number of integrations) after that the two integrated correlations are summed so that the complementary code becomes effective. 655 The new AIS-INGV digital ionosode Board Parameter Values FSY 125 MHz ( fref) reference oscillator stability 15 ppm (0° ∏ 70°) + 1 ppm per year FSY Rx oscillator output frequency (LO#1) 36.9 ∏ 55.9 MHz FSY Tx oscillator output frequency 1 ∏ 20 MHz FSY DDS oscillators frequency resolution 4 ·10– 9 of fref (30 mHz) FSY CRF output level – 12 dBm FSY TF#1 center frequency / 3 dB bandwidth / min. attenuation 46.4 MHz / 28 MHz / 2 dB FSY TF#2 center frequency / 3 dB bandwidth / min. attenuation 40 MHz / 5 MHz / 5 dB FSY TF#3 center frequency / 3 dB bandwidth / min. attenuation 4 MHz / 500 kHz / 8 dB FSY LO#1∏3 output level + 11, + 8, + 5 dBm CTM Pulse compression 16 bit bi-phase complementary codes CTM Sub-pulse duration 30 ms CTM CODE 1 / 0 / Tx off voltages + 1 / – 1 / 0 V PWA RF Power Amplifi er max. gain / 3 dB band limits 68 dB / 0.5 ∏ 50 MHz PWA RF Power Amplifi er max. output / max. linear output + 54 dBm / + 53 dBm PWA RF Power Amplifi er max. spurious level (linear) – 14 dBc @16 MHz, – 27 dBc @2 MHz PWA RF Tx residual output power out of transmission SWF LBI RFF pass band / insertion loss 1 ∏ 25 MHz / 3 dB SWF Rx on/off switches attenuation during Rx / Tx 3 dB / 25 dB SWF 1P6T insertion loss / min. isolation 1 dB / 38 dB SWF NBS RFF central frequency 1.3 - 2.2 - 3.3 - 6.0 - 10 - 17 MHz SWF NBS RFF 3 dB bandwidth 0.6 - 1.1 - 1.8 - 3.0 - 5.0 - 8.5 MHz SWF NBS RFF typ. insertion loss 2 dB SWF Overall SWF gain / attenuation during Rx / Tx 8 dB max. / 36 dB typ. RCV RF#1 center frequency / 3 dB bandwidth / attenuation 35.9 MHz / 2.5 MHz / 5 dB RCV RF#2 center frequency / 3 dB bandwidth / attenuation 4.1 MHz / 200 kHz / 4 dB RCV RF#3 center frequency / 3 dB bandwidth / attenuation 100 kHz / 66 kHz / 4 dB RCV Variable gain amplifi er dynamics / resolution 30 dB / 2 dB RCV Overall RCV board gain (0 dB attenuation) 52 dB ADC ADC input range 1.55 ∏ 3.26 V ADC ADC resolution 8 bit = 6.7 mV ADC Sampling rate 100 kHz ADC Oversampling factor (samples per sub-pulse) 3 ADC Demodulation system Synchronous with quadrature carriers DSP Input data stream 290 ns per sample (16 bit) DSP Processor speed 40 MIPS-80 MHz Table III. Functional blocks characteristics summary. 656 Enrico Zuccheretti, Giuseppe Tutone, Umberto Sciacca, Cesidio Bianchi and Baskaradas James Arokiasamy Fig. 4a,b. Ionograms recorded at 7 a.m. 2002/5/29: a) AIS-INGV ionosonde; b) other digital ionosonde. a b 657 The new AIS-INGV digital ionosode This phase coherent integration of the re- ceived signal is a powerful technique, because random noises are partially cancelled while the in-phase signals are amplifi ed. The processing gain depends on the number of integrations N and this integration remains effective in accordance with the ionospheric coherence time. Lastly an inverse CFFT is performed, the amplitude of the signal is computed and the controlling PC is called, so that it can store and display the data. The processing gain, correlation and phase coherent integration, is about 25 dB. 4. System performance and conclusions After the manufacturing and assembling of the boards and subsystems, tests were done, determining the performance values reported in table III. Unless differently stated all analog lines are matched to 50 W impedance, all digital levels are TTL compatible. The main unit absorbs less than 30 W of electric power (excluding the power amplifi er and the PC) and the weight of the main unit is 9 kg. Total weight and power consumption depend on the particular PC (with monitor) and power amplifi er used to complete the system. At the present stage the total volume is limited ap- proximately: 50 cm width, 50 cm depth, 80 cm height. The ionosonde system has been tested for months, producing ionograms of good quality, comparable with ionograms recorded by other ionosondes. Some results will be presented, making comparisons with the two ionosonde working at INGV. In all ionograms on the hor- izontal axes there are the frequencies in MHz, while the units on the vertical axes are km; they were scaled graphically and placed in the pages to allow an easy comparison among scales; all soundings were performed in the Rome ion- ospheric observatory. Figure 4a is a ionogram recorded by the AIS- INGV ionosonde; fi gs. 4b and 4c are ionograms recorded by a digital and an analogue ionosonde respectively. The soundings were performed within a few minutes around 7.00 a.m. (local time). Figure 5a,b show another sample of iono- grams recorded by the AIS-INGV ionosonde and another digital ionosonde. The shown ionograms were chosen as re- presentative of the ionograms recorded by the various ionosondes; the deviations from the «normal» cases (the ones shown) are rare. It is possible to note the good quality of the iono- grams recorded by the AIS-INGV ionosonde, accompanied by a good background noise level (compared to analog ionosonde). For some conditions, the AIS-INGV ionosonde is even better than the digital ionosonde used or comparison, e.g., at the lower frequencies (E layer). These results encourage us to improve the performance of the ionosonde; some areas in which such improvements may occur are: the new functions to recognise the polarization of the echo, doppler analysis, improvement in S/N ratio and digital signal processing using a fl oating point processor and the software adaptation for remote access. Fig. 4c. Ionogram recorded at 7 a.m. 2002/5/29: other analog ionosonde. 658 Enrico Zuccheretti, Giuseppe Tutone, Umberto Sciacca, Cesidio Bianchi and Baskaradas James Arokiasamy Fig. 5a,b. Ionograms recorded at 5.45 p.m. 2002/5/29: a) AIS-INGV ionosonde; b) other digital ionosonde. a b 659 The new AIS-INGV digital ionosode REFERENCES AROKIASAMY, B.J., C. BIANCHI, U. SCIACCA, G. TUTONE, E. ZUCCHERETTI (2002): The new AIS-INGV digital ionosonde design report, INGV Internal Technical Note. BIBL, K. (1998): Evolution of the ionosonde, Ann. Geofi s., 41 (5/6), 667-680. DAVIES, K. (1990): Ionospheric Radio (P. Peregrinus Ltd., London), pp. 580. HUNSUCKER, R.D. (1991): Radio technique for probing the terrestrial ionosphere, in Physics and Chemistry in Spa- ce, edited by M.C.E. HUBER, L.J. LANZEROTTI and D. STOFFLER (Springer-Verlag, Berlin Heidelberg), vol. 22. OPPENHEIM V.A,. W.R. SCHAFER, J.R. BUCK (1999): Discrete-Time Signal Processing (Prentice Hall, New Jersey), 2nd edition, pp. 870. SKOLNIK, M.I. (1980): Introduction to Radar Systems (Mc Graw-Hill, New York), pp. 422. SKOLNIK, M.I. (1990): Radar Handbook (Mc Graw-Hill, New York), pp. 1200. TEXAS INSTRUMENTS (1994): Digital Signal Processing Ap- plications with the TMS320 Family (Prentice-Hall). (received January 7, 2003; accepted April 28, 2003)