Extended Application of a Novel Phase Calibration Approach of Multiple-Frequency Range Imaging to the Chung-Li and MU VHF Radars

Transcription

1 2488 J O U R N A L O F A T M O S P H E R I C A N D O C E A N I C T E C H N O L O G Y VOLUME 26 Extended Application of a Novel Phase Calibration Approach of Multiple-Frequency Range Imaging to the Chung-Li and MU VHF Radars JENN-SHYONG CHEN Department of Computer and Communication Engineering, Chienkuo Technology University, Changhua, Taiwan CHING-LUN SU AND YEN-HSYANG CHU Institute of Space Science, National Central University, Jhongli, Taiwan GERNOT HASSENPFLUG* Research Institute for Sustainable Humanosphere, Kyoto University, Uji, Japan MARIUS ZECHA Leibniz-Institut für Atmosphärenphysik, Kühlungsborn, Germany (Manuscript received 11 February 2009, in final form 16 June 2009) ABSTRACT Multiple-frequency range imaging (RIM), designed to improve the range resolution of radar echo distribution, is now available for the recently upgraded Chung-Li VHF radar (24.98N, E). To complete the RIM technique of this radar, a novel phase calibration approach, proposed initially for the Ostsee Wind (OSWIN) VHF radar, was employed to examine the effects of phase bias and the range-weighting function on the received radar echoes. The estimated phase bias indicated a time delay of ;1.83 ms for the signal in the radar system. In contrast, such a time delay is more difficult to determine from the phase distribution of twofrequency cross-correlation functions. The same calibration approach was also applied successfully to the middle and upper atmosphere (MU) radar (34.858N, E) and revealed a time delay of ;0.33 ms for the radar parameters employed. These calibration results for various radars demonstrate the general usability of the proposed calibration approach. With the high-resolution performance of RIM, some small-scale Kelvin Helmholtz (KH) billows, double-layer structures, and plumelike structures in the troposphere that cannot be seen in height time intensity plots have been recognized in present observations. The billows and double layers were found to be closely related to strong vertical wind shear and small Richardson number, supporting the hypothesis of a dynamic process of KH instability. On the other hand, the plumelike structures were observed to grow out of a wavy layer and could be attributed to saturation and breaking of gravity waves. These fine structures have shown some remarkable features resolved by the RIM method applied to VHF radars in the lower atmosphere. 1. Introduction * Current affiliation: National Institute of Information and Communications Technology, Tokyo, Japan. Corresponding author address: Jenn-Shyong Chen, Department of Computer and Communication Engineering, Chienkuo Technology University, No. 1, Jieshou N. Rd., Changhua 500, Taiwan. jschen@ctu.edu.tw Range imaging (RIM) with multiple transmitting frequencies has been successfully used with VHF/UHF atmospheric radars, which improves the range resolution of radar echo distribution greatly (the sampling step of RIM can be as small as 1 m). Many studies have been carried out with the help of high-resolution range imaging of the atmospheric structures: for example, Kelvin Helmholtz instability (KHI) (Chilson et al. 2003; Luce et al. 2007), high-resolution wind profiling (Yu and Brown 2004), and vertical velocity bias associated with KHI (Chen et al. 2008). In the literature, RIM is also known as frequency interferometric imaging (FII) (Luce et al. 2001). RIM is based on some inversion algorithms that were initially applied to multiple-receiver coherent radar DOI: /2009JTECHA Ó 2009 American Meteorological Society

2 NOVEMBER 2009 N O T E S A N D C O R R E S P O N D E N C E 2489 imaging (CRI) to resolve multiple echo centers in the radar volume (Hysell 1996; Palmer et al. 1998; Yu et al. 2000). The success of applying the RIM technique to the atmosphere is also based on precise phase calibration of the radar echoes. Therefore, some approaches have been employed for the phase calibration in the RIM processing (Kilburn et al. 1995; Fernandez et al. 2005). Recently, we also proposed a calibration approach for the OSWIN VHF radar (54.18N, 11.88E) (Chen and Zecha 2009, hereafter CZ). The proposed approach is not only practical for the phase calibration of RIM, but also improves the correction of the range-weighting function effect on the imaged powers. To validate the general usability of the proposed calibration approach, we applied this approach to the Chung-Li VHF radar (24.98N, E; Taiwan), which has been upgraded recently, and the middle and upper atmosphere (MU) radar (34.858N, E; Japan). The first range imaging of the Chung-Li radar is achieved after appropriate phase and range-weighting corrections. In this note, we give a brief description of the proposed calibration approach in section 2. In section 3, we present some calibration results of the Chung-Li and the MU radars. The first RIM processing carried out with the Chung-Li radar is displayed, and, moreover, some remarkable KH billows, double-layer structures, and plumelike structures observed by the MU radar are exhibited and examined. Conclusions are given in section Instruments and methodologies The Chung-Li and MU radar systems are introduced briefly in the appendix, along with the equations of RIM used in this study. The issues of phase bias and rangeweighting effect are also described briefly. The calibration approach proposed by CZ is based on a simple hypothesis: the difference between the imaged powers around the edge of two adjacent range gates should be very small after optimal (or effective) range-weighting correction and phase compensation. Accordingly, the following estimator has been utilized to calculate the difference of the imaged powers: ERR 5å N i51 (P 1i P 2i ) 2 P 1i P 2i 5å N i51 P 1i 2 1 P 2i, (1) P 2i P 1i where P 1i and P 2i are the two sets of imaged powers estimated at the same heights above ground, around the edge of two adjacent range gates; N is the number of the imaged powers. ERR varies with the parameter pair of phase compensation and range-weighting function given in calculation: ERR will be zero when P 1i 5 P 2i.IfP 1i 6¼ P 2i, however, ERR becomes larger as the difference between P 1i and P 2i increases. Therefore, we can find the optimal correction for the imaged powers according to the minimum ERR. Since the phase imbalances of different frequency pairs are supposed to be mainly owing to the time delay of the signal in the radar system and the time delay corresponds to a range delay, for conveniencewecandelaytherangeofthesampling gate [variable r in (A3)] in calculating the imaged powers. To compare with the phase bias obtained from the FDI phase distribution (FDI phase: the phase of the cross-correlation function between two echoes of a frequency pair), however, the range delay is transformed into a phase angle by regarding the change of phase within a range-gate interval as Therefore, in practical calculation we use the term (phase bias/ 3608) 3 (gate interval) (m) for the range delay or (phase bias/3608) 3 (pulse length) (ms) for the time delay, in which the phase bias is the variable of phase angle. Note that a so-defined phase bias is not limited to the phase interval between 08 and 3608, but varies from zero to any positive number to reveal the optimal time delay in the radar system. The choice of lower and upper bounds of phase bias in a calculation depends on the radar system. Moreover, the range of the s z values (standard deviation of a Gaussian range-weighting function; see appendix) for calibration can be given according to the radar pulse length employed. For example, a range of m ( m) was chosen for a 1-ms (2-ms) pulse length in our present calculation. The computational step size of phase bias was 108, and the step size of s z was5m(10m)for1-ms (2-ms) pulse length. Owing to inevitable noise, such obtained optimal phase biases (and s z ) for different gate pairs at different times are sometimes not exactly the same. Therefore, the likely phase bias (and s z ) is revealed from the histogram of the optimal phase biases (and s z ). Readers can refer to CZ for a detailed description of the practical calibration procedure. 3. Calibration results Table 1 lists the RIM experiments examined. Cases 1 4 were carried out with the Chung-Li radar, and case 5 was observed with the MU radar. a. Chung-Li radar Figure 1 shows the distributions of the optimal phase biases and optimal s z values of the lower 30 gate pairs for case 1; each distribution contains 100 estimates. As seen, the phase biases of the lower gate pairs were more concentrated, where their mean signal-to-noise ratios (SNRs)

3 2490 J O U R N A L O F A T M O S P H E R I C A N D O C E A N I C T E C H N O L O G Y VOLUME 26 TABLE 1. Radar parameters used in RIM experiments. Cases 1 4: Chung-Li radar and case 5: MU radar. Cases 2 4 are carried out alternately, and each mode is ;110 s. The sampling time resolution is ;0.256 s for cases 1 4 and ;0.512 s for case 5. Case Start end time/date Pulse shape a Rx filter bandwidth (khz) Pulse code Range resolution (m) Sampling step, z o (m) b Frequency set c UTC S 500 No , 1800 f a 30 Mar UTC G 250 No , 1800 f a 11 Apr UTC S 250 No , 1800 f a 11 Apr UTC G 500 No , 1800 f a 11 Apr UTC 8 Feb S bit Spano code , 3000 f b 0700 UTC 10 Feb 2006 a G: Gaussian, S: square. b z o : the start of sampling height. c f a 5 [51.75, , 52.0, , 52.25] MHz (Df MHz); f b 5 [46.0, 46.25, 46.5, 46.75, 47.0] MHz (Df MHz). were larger, too. Rough inspection can tell that the likely phase bias was ;3208. On the other hand, the phase biases of the upper gate pairs were dispersive and their mean SNRs were small; usually the results at low SNR are less reliable. As for the optimal s z shown in the right panel, the distributions of the lower gate pairs centered on ;170 m, which can be regarded as the likely value of s z that is close to the theoretical one. Nevertheless, we find that lower SNR leads to a larger and more dispersive value of s z, as observed clearly at the upper gate pairs. In theory, s z is not related to SNR; however, the echoes at low SNR are statistically more random, causing the value of s z to be also SNR dependent. To obtain smoother imaged powers at gate boundaries, therefore, a value of s z adaptive to SNR can be employed. The proposal of the SNR-dependent s z is one main innovation of the proposed calibration approach. The distributions of the optimal phase biases (and s z ) of different gate pairs can be merged to form a more reliable histogram, as shown in Fig. 2. Figures 2a d display the results of cases 1 4, respectively. As seen, the peak locations of the four histograms are close (Fig. 2, left); that is, between 3208 and Taking 3308 as the representative phase bias, this indicates a time delay of ;1.83 ms (3308/ ms). As for the histogram of s z, we can find differences in the peak locations of the four cases (right panels): ;170, ;260, ;220, and ;210 m for cases 1 4, respectively, which indicate approximately the theoretical value of s z in the absence of noise. According to the pulse shape and filter bandwidth employed in the experiment, the Gaussian range-weighting function for case 1 is narrowest, and that for case 2 is broadest. In view of this, the likely values of s z indicated in Fig. 2 are reasonable. It should be mentioned that all the characteristics of the s z histogram shown here are similar to those of the OSWIN VHF radar demonstrated by CZ. Moreover, the values of optimal s z (not the likely s z ) are also SNR-dependent, and Eq. (3), from CZ, can be employed again to find the mean relationship between optimal s z and SNR. The above calibration suggests a time delay of ;1.83 ms for the Chung-Li radar system. It is expected that such a time delay can also be determined from the FDI phase distribution. Theoretically, the distribution of FDI phases of each gate is a quasi-gaussian shape with the peak located at the phase angle predicted for the center of the radar volume; such a shape is related to the range-weighting function as well as the intensity distribution of the radar main beam. However, the unbiased peak location occurs only when no phase imbalance exists between echoes at two different carrier frequencies. As presented in CZ, the peak location of the FDI phase distribution can be estimated in advance according to the radar parameters used. For our present experiments, the extents of FDI phases and their peak locations are summarized in Table 2 for different frequency separations Df. Note that the phase extents and peak locations of Df khz (250 khz) and 375 khz (750 khz) are different for the four grouped gates; for Df khz (500 khz) the FDI phases dwell in different phase extents at odd and even gates. As for Df khz (1000 khz) the FDI phases always lie between 08 and 3608 with the peak at Taking case 1 as an example, Fig. 3 exhibits the FDI phases observed, where Fig. 3a d display the respective FDI phase distributions of the grouped gates for four frequency pairs. Other frequency pairs result in similar features to those displayed here as long as their frequency separations are the same; for brevity, they are not shown here. Referring to Table 2, an inspection of

4 NOVEMBER 2009 N O T E S A N D C O R R E S P O N D E N C E 2491 FIG. 1. Distributions of (left) optimal phase biases and (right) optimal s z of different gate pairs. Each distribution is self-normalized. the distributions in Fig. 3 reveals that the differences (phase imbalance) of observed and predicted peak locations are variable: they are about 2458, 1608 (odd gates) and (even gates), 808 and and 808 and 708, and and and 2208 for the frequency pairs Df 5 500, 250, 125, and 375 khz, respectively. Such a set of phase imbalances is peculiar, different from that of the OSWIN radar, and thus needs further clarification. In fact, the negative values of phase imbalances are due to phase ambiguity. Therefore, we conclude that the phase imbalances are about 3158, 1608, 808, and 2308 for Df 5 500, 250, 125, and 375 khz, respectively. Such a quasi-linear variation in phase imbalances corresponds to a time delay of ;1.78 ms, which is close to the ;1.83 ms revealed by the proposed calibration approach. The above examination suggests again that we should be careful sometimes in using the FDI phase distribution for calibration; phase ambiguity may occur when the time delay of the signal in the radar system is large. In comparison and to its advantage, the calibration approach as used does not encounter this problem of phase ambiguity. Moreover, one may wonder if the time delay is so large that it could be ; n ms, where n is zero or a positive integer. In the case of n 5 1 or 2, the time delay is ;2.83 or ;4.83 ms, corresponding to a phase imbalance of ;6908 ( )or;10508 ( ) for Df khz. We examined this and concluded that such large time delays can be discarded, as demonstrated in Fig. 4. In Fig. 4, we cannot find representative peaks of phase biases (the peaks at the leftmost edge of

5 2492 J O U R N A L O F A T M O S P H E R I C A N D O C E A N I C T E C H N O L O G Y VOLUME 26 FIG. 2. Histograms of (left) optimal phase biases and (right) optimal s z. (a) (d) Cases 1 4 listed in Table 1, respectively. The phase bin is 108, and the s z bin is 10 m.

6 NOVEMBER 2009 N O T E S A N D C O R R E S P O N D E N C E 2493 TABLE 2. Chung-Li radar and MU radar (in brackets) for Df range from 500 to 125 khz. 500 khz 250 khz 125 khz 375 khz [1000 khz] [500 khz] [250 khz] [750 khz] Extent Peak Extent Peak Extent Peak Extent Peak Gates (deg) (deg) (deg) (deg) (deg) (deg) (deg) (deg) 1, 5, 9, [0 360] [180] [0 180] [90] [0 90] [45] [0 270] [135] 2, 6, 10, [0 360] [180] [ ] [270] [90 180] [135] [ ] [45] 3, 7, 11, [0 360] [180] [0 180] [90] [ ] [225] [ ] [315] 4, 8, 12, [0 360] [180] [ ] [270] [ ] [315] [90 360] [225] the abscissa are unreliable), and the optimal values of s z are not appropriate. Based on the above calibration results, one of the RIM observations carried out with the Chung-Li radar is shown in Fig. 5. Figure 5a discloses some echo layers at a range resolution of 300 m; these echo layers are range smeared. Figure 5b displays the imaged powers without calibration; as seen, we cannot resolve the echo layers correctly. In comparison, Fig. 5c shows that the range imaging is improved greatly after proper corrections of phase bias and range-weighting effects. As seen, thin layered and turbulent/convective structures can be identified above and below the height of ;4 km, clearly indicating different atmospheric characteristics in the middle troposphere and in the boundary layer. In the literature, similar structures have been observed and addressed in detail with the MU radar (e.g., Luce et al. 2006, 2007). b. MU radar The general usability of the proposed calibration approach was further tested using MU radar data, as FIG. 3. For case 1: distributions of FDI phases of four frequency pairs with different frequency separations Df. The phase bin is 108.

7 2494 J O U R N A L O F A T M O S P H E R I C A N D O C E A N I C T E C H N O L O G Y VOLUME 26 FIG. 4. As in Fig. 2, but the interval of phase bias of the abscissa is changed. The data are from case 1. demonstrated in Fig. 6. Figure 6a indicates that the phase bias is ;1208, corresponding to a time delay of ;0.33 ms (1208/ ms). The likely value of s z is ;90 m, and the values of s z are also SNR dependent (not shown). Figures 6b e show some FDI phase distributions for several frequency pairs; the other distributions, not shown, possess similar features to the panels here, depending on their frequency separations. Referring to the predicted peak locations listed in Table 2, the phase imbalances can be estimated from the mean peak locations of the humps, except for the panel of Df khz: ;608, ;258, and ;908 for Df 5 500, 250, and 750 khz, respectively. A quasi-linear relationship between the phase imbalance and the frequency separation is evident, although the FDI phase distribution of Df khz does not show an observable hump. The cause of the concealed hump for Df khz is currently unknown. Again, the proposed calibration approach can avoid the need to deal with such an indistinct phase distribution. Figure 7 exhibits a special RIM example from our observations. Note that the time starts from the right side of the abscissa, which makes it more convenient to compare with the wind field depicted later. As seen, the diagram of the received power shows nothing exciting (Fig. 7a), and the Capon imaged powers without any calibration are range-biased severely (Fig. 7b). After appropriate calibration, however, some KH billows are revealed between 6.0 and 6.6 km, during h (Fig. 7c). Before and after these KH billows are doublelayer structures. Note that the observations were interrupted for tens of seconds each 4-min interval for data storage, making the imaged power vary abruptly sometimes (e.g., at 1200 UTC). The close relationship between double-layer structures and KH billows has been addressed in detail (Browning and Watkins 1970; Worthington and Thomas 1997) and has also been demonstrated by simulations (Fritts et al. 2003, 2006). Vertical wind shear is the most important factor responsible for KH billows, and this can be recognized from Fig. 7d, where the horizontal wind field, derived from full-correlation analysis (FCA), is shown. As seen, significant vertical wind shear existed between 6.0 and 6.6 km. An estimate of mean Richardson number during the period displayed here is about 0.1 ; 1.0, by taking the values of ; rad 22 s 22 for the square of the assumed Brunt Väisälä frequency (a typical range reported in the literature). This magnitude of Ri is around the threshold (0.25) required for the occurrence of active KH billows. Recently, Luce et al.

8 NOVEMBER 2009 N O T E S A N D C O R R E S P O N D E N C E 2495 FIG. 5. A portion of the observations of case 1. (a) Received radar echoes with a range resolution of 300 m. (b) Range imaging without phase and range-weighting corrections. (c) Range imaging with adaptable s z (varying with SNR) and phase correction (3308). Sampling step of range is 2 m in the imaging process and time resolution is ;30 s. (2008) reported some KH billows above the jet stream and attributed it to the trigger of an inertia gravity wave generated by the jet stream. Different from that case, the KH billows shown in Fig. 7 are far away from the core of the jet stream and persist only for a short period; in view of this, it seems that our case here is mainly related to the local wind shear. According to the authors knowledge, the case shown in Fig. 7 could be the clearest observation so far of the relationship between KH billows, double-layer structures, and vertical wind shear demonstrated by RIM with a VHF radar. Figure 7c also discloses a double-layer structure around the height of 6.8 km. In fact, such a structure appears quite often in our present MU radar data. In the literature, the model of turbulent layers addressed by Woodman and Chu (1989) can be one interpretation of the double-layer structure observed; that is, highly anisotropic turbulence located at both edges of a turbulent layer can reflect/ scatter the radar echoes much stronger than the refractive index irregularities inside the main and more wellmixed turbulent layer. In addition to this model, we suspect that small-scale KH billows may occur very often and intermittently and could be the main mechanism of double-layer structures observed within a small vertical range. Unfortunately, even the RIM used with VHF radars cannot resolve such small KH billows yet; only the consequent double-layer structure is recognizable. One more novel observation with the MU radar is shown in Fig. 8. The received powers with 150-m resolution do not reveal any notable features (Fig. 8a). After

9 2496 J O U R N A L O F A T M O S P H E R I C A N D O C E A N I C T E C H N O L O G Y VOLUME 26 FIG. 6. For case 5 (MU radar). (a) As in Fig. 2. (b) (e) Distributions of FDI phases of four frequency pairs with different frequency separations Df. The phase bin is 108. range imaging, however, two thin layers are disclosed around the heights of 5.2 and 6.2 km. (Note again that the observations were interrupted for tens of seconds each 4-min interval, causing abrupt change at times in the imaged power between data files.) The lower layer is wave modulated and its shape is slightly tilted due to vertical wind shear, which has been examined for bias of vertical velocity by Chen et al. (2008). Existence of atmospheric waves can be revealed from the vertical wind field shown in Fig. 8c. As shown, vertical velocities below the height of ;5.5 km indicate a wave activity having a period of ;4 min, especially between 20.1 and 20.4 h. However, the wave activity stopped at the height of ;5.5 km, indicating a possible critical level at ;5.5 km for this wave. An absorbing phenomenon is the plumelike structures between the two thin layers, which seem to grow out of the lower wavy layer. We suspect that these plumelike structures arose from saturation and breaking of gravity waves. This has shown a remarkable feature observed only by the RIM of a VHF radar in the lower atmosphere so far. 4. Conclusions An extended employment of a novel calibration approach for range imaging (RIM), which was initially applied to the OSWIN VHF radar recently, was made in

10 NOVEMBER 2009 N O T E S A N D C O R R E S P O N D E N C E 2497 FIG. 7. For case 5 (MU radar). (a) Received radar echoes with a range resolution of 150 m. (b) Ranging imaging without phase and range-weighting corrections. (c) Range imaging with adaptable s z (varying with SNR) and phase correction (1208). Sampling step of range is 1 m in the imaging process and time resolution is ;32 s. (d) Horizontal wind field. this study for both the Chung-Li and the MU VHF radars. The first RIM experiment with the Chung-Li radar was also carried out successfully. Our studies have suggested that the major cause of the phase biases of RIM could be the same for the three radar systems mentioned: time delay of the signal in the system, which leads to shift in FDI phase. Moreover, it is revealed that the calibration using the FDI phase distributions of different frequency pairs could be more difficult when the time delay in the radar system is large,

11 2498 J O U R N A L O F A T M O S P H E R I C A N D O C E A N I C T E C H N O L O G Y VOLUME 26 FIG. 8. For case 5 (MU radar). (a) Received radar echoes with a range resolution of 150 m. (b) Range imaging with adaptable s z (varying with SNR) and phase correction (1208). Sampling step of range is 1 m in the imaging process and time resolution is ;32 s. (c) Vertical winds: gray (black) indicate upward (downward) velocities. as happened in the case of the Chung-Li radar. By comparison, the proposed calibration approach is capable of indicating the value of time delay without the problem of phase ambiguity. With the high-resolution performance of RIM, a double-layer structure accompanied by some KH billows has been observed clearly with the MU radar, and their formation is closely related to significant vertical wind shear and small Richardson number, supporting the hypothesis of a dynamic process of KH instability. Moreover, some plumelike structures associated with a wavy layer are first reported in this note and are suspected to be the consequence of saturation and breaking of gravity waves. Based on these initial demonstrations, the RIM technique used with VHF/UHF pulsed radars indeed aids in the physical interpretation of atmospheric phenomena by displaying high-resolution radar echoes in the range direction. More scientific studies on the atmosphere can thus be expected in the future in connection with microstructures and dynamics of the atmosphere. Acknowledgments. This work was supported by the National Science Council of ROC (Taiwan) under Grants NSC M MY3 and NSC M The Chung-Li radar is maintained by the Institute of Space Science, National Central University, ROC (Taiwan), and the MU radar is operated by the Research Institute for Sustainable Humanosphere, Kyoto University, Japan.

12 NOVEMBER 2009 N O T E S A N D C O R R E S P O N D E N C E 2499 APPENDIX Radar Systems and Range Imaging a. Radar systems Figure A1 shows the configurations of the main antenna arrays of the Chung-Li and MU VHF radars. The Chung- Li radar is operated at a central frequency of 52 MHz, in which the ST array is used mainly for observations of the lower atmosphere (;2 20 km) and the IONO array is for the ionosphere. The ST and IONO arrays consist of three subarrays formed by and four-element linear Yagi antennas, which are connected to three identical but independent transmitter receiver modules. The upgraded oscillator can currently generate carrier frequencies from to MHz for RIM experiments (0.5-MHz imaging bandwidth so far; a larger imaging bandwidth is possible but it needs further testing). Moreover, radar pulse shape can be Gaussian or rectangular, and various filter bandwidths are available. [Currently, the source for more radar characteristics is via at clsu@jupiter.ss.ncu. edu.tw or from ChungLiVHF_home.htm (presented only in Chinese).] The MU radar, operated at a central frequency of 46.5 MHz and 1-MHz imaging bandwidth from 46 to 47 MHz, possesses an array consisting of 475 crossed- Yagi antennas, and each antenna is equipped with an independent but identical type of transmitter receiver module. The radar array is divided into 25 independent groups with 19 crossed-yagi antennas each. Each antenna group has its own digital receiver for the combined signal. Furthermore, software can combine the outputs of these receivers arbitrarily. For our present study, the combined output of the full antenna array was employed. Readers can refer to Hassenpflug et al. (2008) for a detailed description of the MU radar system. b. Range imaging Capon s method (Capon 1969) has been demonstrated to be excellent for RIM (Palmer et al. 1999; Luce et al. 2001; Chilson et al. 2003), and one set of the Capon equations can be simply expressed as in which 1 P(r) 5 e H R 1 e 2 3 R 11 R 12 R 1n R 21 R 22 R 2n R R n1 R n2 R nn (A1) (A2) FIG. A1. Configurations of the Chung-Li and MU VHF radars. and e 5 [e j2k 1 r,e j2k 2 r,...,e j2k n r ] T, (A3) where P(r) is the imaged power at range r; the superscripts H and 21 in (A1) and T in (A3) represent the Hermitian, inverse, and transpose operators, respectively, and k n is the wavenumber of the nth carrier frequency. The term R mn is the nonnormalized cross-correlation function of the signals calculated at zero time lag for a pair of transmitting frequencies; the phase of R mn is termed frequency-domain-interferometry (FDI) phase in this paper. Equations (A1) (A3) present the processing without Doppler frequency sorting of the echoes. More algorithms used with RIM and their comparisons can be found in the literature (e.g., Luce et al. 2006). The phase of R mn may be biased owing to the phase imbalance between the two echoes of a frequency pair. Such phase bias/imbalance creates offsets in the imaged power location and/or smears the image, making true

13 2500 J O U R N A L O F A T M O S P H E R I C A N D O C E A N I C T E C H N O L O G Y VOLUME 26 physical phenomena like KH billows hard to see or creating the impression of false phenomena, and thus need to be compensated. To find suitable phase compensation is essential. According to our previous studies (Chen 2004; CZ), the reason for phase biases of different frequency pairs is suspected to be mainly the time delay of the signal in the radar system. The exact cause of such time delay is unknown, but possibly originates from many parts of the radar system: response times of the elements, mixing and filtering of the echoes, propagation delay of the signals in the cable lines/circuits, among other causes. In addition to phase bias, the range-weighting effect on P(r) may also be considerable, especially toward the edges of the gates, and should be removed. This can be accomplished approximately by using the inverse of the range-weighting function W 2 (r) 5 exp(2r 2 /s z 2 ) (Franke 1990; Luce et al. 2001). Theoretically s z can be estimated from the radar pulse shape and filter bandwidth; however, considering the ubiquitous noises embedded in the echoes, the use of a noise-dependent s z has been recommended by CZ to improve the continuity of the imaged powers at gate boundaries. REFERENCES Browning, K. A., and C. D. Watkins, 1970: Observations of clear air turbulence by high power radar. Nature, 227, Capon, J., 1969: High-resolution frequency-wavenumber spectrum analysis. Proc. IEEE, 57, Chen, J.-S., 2004: On the phase bias of multiple-frequency radar returns of mesosphere-stratosphere-troposphere radar. Radio Sci., 39, RS5013, doi: /2003rs , and M. Zecha, 2009: Multiple-frequency range imaging using the OSWIN VHF radar: Phase calibration and first results. Radio Sci., 44, RS1010, doi: /2008rs , G. Hassenpflug, and M. Yamamoto, 2008: Tilted refractiveindex layers possibly caused by Kelvin Helmholtz instability and their effects on the mean vertical wind observed with multiple-receiver and multiple-frequency imaging techniques. Radio Sci., 43, RS4020, doi: /2007rs Chilson, P. B., T.-Y. Yu, R. G. Strauch, A. Muschinski, and R. D. Palmer, 2003: Implementation of range imaging on the Platteville 915-MHz troposphere profiler. J. Atmos. Oceanic Technol., 20, Fernandez, J. R., R. D. Palmer, P. B. Chilson, I. Häggström, and M. T. Rietveld, 2005: Range imaging observations of PMSE using the EISCAT VHF radar: Phase calibration and first results. Ann. Geophys., 23, Franke, S. J., 1990: Pulse compression and frequency domain interferometry with a frequency hopped MST radar. Radio Sci., 25, Fritts, D. C., C. Bizon, J. A. Werne, and C. K. Meyer, 2003: Layering accompanying turbulence generation due to shear instability and gravity-wave breaking. J. Geophys. Res., 108, 8452, doi: /2002jd , P. M. Franke, K. Wan, T. Lund, and J. Werne, 2006: Comparison of radar backscatter moments with real turbulence and large-scale flow parameters throughout the lifecycle of a Kelvin-Helmholtz instability. Technical and Scientific Aspects of MST Radar (MST 11), Macmillan Publishers India, 14. Hassenpflug, G., M. Yamamoto, H. Luce, and S. Fukao, 2008: Description and demonstration of the new Middle and Upper atmosphere Radar imaging system: 1-D, 2-D and 3-D imaging of troposphere and stratosphere. Radio Sci., 43, RS2013, doi: /2006rs Hysell, D. L., 1996: Radar imaging of equatorial F region irregularities with maximum entropy interferometry. Radio Sci., 31, Kilburn, C., S. Fukao, and M. Yamamoto, 1995: Extended period frequency domain interferometry observations at stratospheric and tropospheric heights. Radio Sci., 30, Luce, H., M. Yamamoto, S. Fukao, D. Hélal, and M. Crochet, 2001: A frequency domain radar interferometric imaging (FII) technique based on high-resolution methods. J. Atmos. Sol.- Terr. Phys., 63, , G. Hassenpflug, M. Yamamoto, and S. Fukao, 2006: Highresolution vertical imaging of the troposphere and lower stratosphere using the new MU radar system. Ann. Geophys., 24, ,,, M. Crochet, and S. Fukao, 2007: Range-imaging observations of cumulus convection and Kelvin-Helmholtz instabilities with the MU radar. Radio Sci., 42, RS1005, doi: /2005rs ,,, and S. Fukao, 2008: High-resolution observations with MU radar of a KH instability triggered by an inertia gravity wave in the upper part of a jet stream. J. Atmos. Sci., 65, Palmer, R. D., S. Gopalam, T.-Y. Yu, and S. Fukao, 1998: Coherent radar imaging using Capon s method. Radio Sci., 33, , T.-Y. Yu, and P. B. Chilson, 1999: Range imaging using frequency diversity. Radio Sci., 34, Woodman, R. F., and Y.-H. Chu, 1989: Aspect sensitivity measurements of VHF backscatter made with the Chung-Li radar: Plausible mechanisms. Radio Sci., 24, Worthington, R. M., and L. Thomas, 1997: Long-period unstable gravity-waves and associated VHF radar echoes. Ann. Geophys., 15, Yu, T.-Y., and W. O. J. Brown, 2004: High-resolution atmospheric profiling using combined spaced antenna and range imaging techniques. Radio Sci., 39, RS1011, doi: /2003rs , R. D. Palmer, and D. L. Hysell, 2000: A simulation study of coherent radar imaging. Radio Sci., 35,