## 1. Introduction

At low microwave frequencies, radars respond to spatial variations in the index of refraction of the air often characterized by the atmospheric refractive-index structure parameter,

Microwave backscattering from the clear-air ABL is generally a combination of Bragg scattering from small-scale spatial variations in refractive index caused by turbulence and Rayleigh scattering from insects, birds, dust, or other airborne particles whose size is much smaller than the radar wavelength. At S-band frequencies, radars measure clear-air backscatter from both refractive index turbulence and from insects. Based on the experimental observations (Wilson et al., 1994; Eaton et al., 1995; Ince et al., 2003; Martin and Shapiro, 2007), Rayleigh scattering from insects can be the dominant source of clear-air returns in both the convective boundary layer (CBL) and the nocturnal boundary layer (NBL). Depending on an application, radar return from airborne insects may be considered either as source of clutter (such as contamination of weather radar measurements) or as desirable passive tracers of the clear air (Wilson et al., 1994). Due to its very fine spatial and temporal resolution, S-band FMCW radars are able to resolve individual insects and discriminate between discrete particle echoes and distributed clear-air echo.

Despite FMCW radars being excellent tools for turbulence studies in the convective ABL and in the capping inversion region (Eaton et al., 1995), relatively few atmospheric studies reported in the literature have used FMCW radar data extensively for both quantitative and qualitative analysis of clear-air features observed. While S-band FMCW radars have been designed with the capability to obtain height and time resolutions of 1 m and 1 s, respectively, the degree to which these resolution limits are obtained in practice depends upon the properties of the atmospheric echo itself. Range measurement errors due to non-zero Doppler velocities and finite coherence of the clear-air echo have implications on both spatial resolution and Doppler estimation of FMCW radars. Although the theoretical background for FMCW radars is well established, their measurement limitations for atmospheric targets and performance impact for atmospheric applications has received limited attention in the literature.

The detailed review of the theory of operation of S-band FMCW radar and its application to atmospheric boundary layer profiling can be found in the literature (Eaton et al., 1995; Ince et al., 2003). In this study, performance aspects of atmospheric FMCW radars including the effects of range-Doppler ambiguity, parallax typical of two-antenna radar systems, and near-field operation are discussed. The data collected during recent field experiments by the high-resolution S-band FMCW radar, described in Ince et al. (2003), is used to illustrate system performance. Analysis of FMCW radar signatures of atmospheric targets detectable at S-band frequencies has been performed by extracting both qualitative morphological information and quantitative information on the intensity of backscattered power from clear-air turbulence.

## 2. Principles

Frequency-modulated, continuous-wave (FMCW) radars may be thought of as a limiting case of pulse-compression radar where the duty cycle of the transmitted waveform approaches 100%. They operate by transmitting a long, coded waveform of duration *T* and bandwidth *B*. The improvement factor they gain over pulsed radars of equivalent range resolution is given by the time-bandwidth product of the waveform *BT*, which is often referred to as the compression gain. In FMCW systems, this gain can be very large, exceeding 60 dB. While several types of frequency coding may be used to yield the bandwidth *B*, linear frequency modulation is the simplest and most commonly used method in atmospheric FMCW radars.

Consider an FMCW radar that transmits constant amplitude linear FM signal of the form

where *w* is the radian frequency and *a* is the chirp rate in rad/s^{2}. The echo from an atmospheric target is essentially a delayed, attenuated, and possibly Doppler-shifted replica of the transmitted signal. The echo is demodulated by mixing it with a portion of the transmitted signal and low-pass filtering the result. The resulting beat frequency for a point target at range

where

Signal processing of the echo usually involves matched filtering, which is most commonly implemented for all ranges simultaneously through spectral analysis via an FFT algorithm. In this case, the output of the matched filter, which gives the response of the linear FMCW radar to a moving point target, can be expressed as

where

For FMCW Doppler radar, Doppler information can be retrieved on a sweep-to-sweep basis by analyzing the sequence of echoes from a particular range, as discussed by Strauch et al. (1976). In this case, the sampling frequency is the reciprocal of the sweep period *T*, the Nyquist Doppler frequency is 1/2*T*, and resulting in an unambiguous Doppler velocity interval of

The received mean power of a bistatic weather radar from a distributed atmospheric target of uniform reflectivity

where

This equation assumes that the Bragg wavenumber,

## 3. Performance Analysis

### 3.1. Range and Doppler

FMCW radar performance characteristics, such as range-Doppler coupling, can be seen from analysis of (3), which is also ambiguity function for linear FMCW signal. The presence of both range and Doppler terms in the argument of the sinc function illustrates the effect of target motion on the radar’s ability to locate. Figure 1 shows a typical two-dimensional time delay (range) – Doppler (velocity) ambiguity surface for a linear FMCW signal with time-bandwidth product of 1000. As seen from Figure 1, maximum unambiguous Doppler can be increased (by decreasing sweep period) at the expense of range ambiguity, or alternatively increasing maximum range results in Doppler ambiguity. By rearranging the terms in the argument in (3), the apparent range of the target can be expressed as

Additionally, the coherence of the atmospheric target during the sweep interval will limit range resolution of linear FMCW radar. For complex moving targets and for volume scattering, the coherence time (or the reciprocal of the Doppler spectral width) of the echo will limit resolution. A distribution of Doppler velocities observed over an integration time, *T*, will yield a distribution of apparent ranges. The resulting rms spread in range is determined by the transformation of the Doppler spectrum to the range domain using the relationship

where

To maximize resolution and sensitivity, it is desirable to make both bandwidth and sweep time as large as possible. The effective value of sweep time is however constrained by the coherence time of the atmospheric echo. For example, a sampling volume with rms radial velocity of 1 m/s has a spectral width of 20 Hz at 3 GHz implying a coherence time of approximately 50 ms.

### 3.2. Parallax and Near-field Operation

From the weather radar equation (4), the received mean power can be correctly estimated only if radar targets are in the far field of antenna, that is*D* is the diameter of the antenna. Additionally, for bistatic FMCW radar employing two spatially separated antennas, a correction of the received power is necessary at near ranges to account for the reduced beam overlap. For Gaussian shaped beams aligned with their axes parallel, the antenna parallax function, or fractional beam overlap is given by Ince et al. (2003),

where *d* is the separation distance between transmit and receive antennas, and

FMCW atmospheric radars can measure boundary layer targets at minimum distances as close as 50 m. However, this minimum range lies in the Fresnel region in which antenna gain and pattern shape vary significantly with distance, and from (4) radar reflectivity is underestimated for targets in the near field. Based on the universal near-field reflectivity correction by Sekelsky (2002), 1 dB reduction in reflectivity at

## 4. Experimental Measurements

Data collected by the University of Massachusetts’ high-resolution FMCW radar (Figure 3) during recent field experiments is used to illustrate system performance and analyze radar signatures of atmospheric targets detectable at S-band frequencies. Detailed description of this system and data processing can be found in Ince et al. (2003). The S-band FMCW radar system implements an internal calibration loop (by injecting into the receiver an attenuated and delayed sample of the high power amplifier output) to perform calibration of the atmospheric echo and estimate volume reflectivity.

Figure 4 shows one day continuous radar record of diurnal cycle of the ABL beginning at 14:27 local time (CDT, 19:27 UTC). The radar echo is expressed in terms of the logarithm of microwave

From Figure 4, four-hour time period of convective ABL (between 14:30 to 18:30 local time) is extracted to perform quantitative analysis of FMCW radar reflectivity. In this case, the collocated radiosonde (operated by NCAR/ATD) measurements of temperature and humidity show typical characteristic of convective ABL bounded above by dry air. In the three panels of Figure 5 respectively reflectivity, Doppler (vertical) velocity, and the correlation coefficient of successive echoes for one-hour of convective ABL are shown. The latter two products are the result of the pulse-pair processing (Doviak & Zrnic, 1993) averaged over 20 pulses (~1 s averaging).

The reflectivity image shows that initially significant Rayleigh backscatter is observed both above and below the capping inversion which peaks near 1000 m altitude at about 17:00 local time. After this time, both the Rayleigh scatter and the distributed Bragg scatter below the inversion decrease significantly, and the strong echo at the top of the boundary layer disappears. The velocity image which is derived from the phase of the single-lag covariance shows structure for some, but not all of the clear-air echo (due to velocity aliasing). The single-lag covariance panel shows extremely high correlations for the Rayleigh scattering while the correlation coefficient for the Bragg scattering is relatively lower on average (due to larger spectral width of clear-air echo).

To discriminate the distributed clear-air echo from particulate scatter, a signal processing method based on the correlation coefficient of successive echoes can be applied to time-series of backscattered power. According to this method, magnitude of the sweep-to-sweep correlation coefficient is compared to an empirically determined threshold (in this case, a value of 0.8 is used) and high-correlation points (due to Rayleigh scatter) are removed and filled with estimated data by using least-squares interpolation (without affecting the statistics of the atmospheric clear-air echo). This processing allows estimation of atmospheric component of the clear air backscatter. Alternatively, median filtering can be applied to the two-dimensional radar reflectivity image to differentiate morphologically differing Bragg scatter (distributed echoes) and Rayleigh scatter (dot echoes). However, in case of large density of particle echoes, their morphology will be similar to distributed scatter, hence the filter will not be able to separate two sources of scattering.

In Figure 6, the computed vertical profiles of mean

In this case, FMCW radar estimates of*z* is height above ground level (Wyngaard & LeMone, 1980). From Figure 6, radar estimated

In Figure 7, log-variance of estimated

In Figure 8, radar reflectivity time-height image (starting at 0650 UTC in Figure 4) of nocturnal boundary layer turbulence below and insect layer within a low-level jet in the NBL is shown. Collocated radiosonde profile of wind speed taken at 0658 UTC is shown in Figure 9. The FMCW reflectivity shows intensified turbulence and wave activity from surface extended to heights of about 150 m where sounding shows evidence of strong shear. Additionally, radiosonde wind measurements shows a low-level jet formed between heights 150 and 550 m with maximum winds of over 16 m/s at approximately 300 m height. In this case, a layer of insects detected by the radar at heights of the wind speed maximum helps reveal qualitatively this additional boundary layer structure.

## 5. Conclusion

In this study we have discussed and illustrated system performance of FMCW radar detection of atmospheric targets at S-band frequencies. In particular, the effects of non-zero Doppler velocities and finite coherence of the atmospheric echo on both spatial resolution of FMCW radars and on Doppler estimation have been presented. Additionally, we have considered parallax errors in reflectivity and near-field operation. For the clear-air atmosphere measurements, the S-band radar’s sensitivity to Bragg and Rayleigh scattering and the effect of Rayleigh scatterers on vertical profiles of mean radar reflectivity have been illustrated. Experimental data collected by the high-resolution atmospheric FMCW radar has been used to perform analysis of S-band radar signatures of atmospheric targets and to quantitatively test theoretical predictions on the atmospheric boundary layer.