1. Introduction
Clouds play a major role in Earth’s hydrological cycle and radiation budget (Boucher et al. 2013); thus, it is important to characterize them with high accuracy and spatial resolution. Reliable and robust ground-based Doppler radars and radiometers are key instruments for detailed observations of cloud processes, long-term data collection, and model and satellite evaluations. Here, we present a system [Jülich Observatory for Cloud Evolution (JOYCE) Radar-94 GHz (JOYRAD-94)] that combines frequency-modulated continuous wave (FMCW) Doppler radar with passive radiometer observations at millimeter wavelengths (94 and 89 GHz, respectively), providing new opportunities for cloud observations.
Radars operating at millimeter-wavelength can be used—among other applications—to determine cloud microphysical processes (e.g., Kollias et al. 2007a; Luke and Kollias 2013; Kneifel et al. 2011; Acquistapace et al. 2017; Tridon and Battaglia 2015) to estimate rain rates (Maahn et al. 2015), or to obtain information about in-cloud turbulence (e.g., Kollias et al. 2002; Borque et al. 2016).
Moreover, radars are combined with other remote sensors in order to enhance the sensitivity and the accuracy of the derived cloud properties. For example, Frisch et al. (1998) used microwave radiometer (MWR) and radar data to retrieve the liquid water content (LWC) profile of a cloud. The vertically integrated amount of liquid water [also liquid water path (LWP)] was constrained by the MWR, and the range-resolved radar signal was used to reconstruct the liquid water vertical distribution. Illingworth et al. (2007) exploited the synergy of collocated radar, microwave radiometer, and lidar observations to obtain further refined cloud properties, such as the hydrometeor phase, type, and position in the vertical column (Cloudnet algorithm).
Currently, the majority of cloud radars—especially at comprehensive observatories (e.g., Löhnert et al. 2015; Kollias et al. 2007b; Haeffelin et al. 2005)—are pulsed systems. Illingworth et al. (2015) pointed out that FMCW radars can be a more economical solution, which then could further increase the number of observatories for comprehensive cloud studies that are needed for climate monitoring. However, FMCW radar requires two antennas, which could in turn increase system costs, especially when considering low-frequency radars where larger antenna diameters are required (Heijnen et al. 2000).
A further disadvantage of FMCW technology is that it requires a large bandwidth at great distance to sustain vertical and temporal resolution, which increases the background noise. Nevertheless, several publications showed that FMCW W-band radars are suitable for both ground-based (Yamaguchi et al. 2009; Huggard et al. 2008; Bennett et al. 2009; Thies et al. 2010; Delanoë et al. 2016) and airborne fog and cloud studies (Pazmany et al. 1994; Mead et al. 2003).
FMCW radars often allow greater flexibility of varying vertical resolution. For JOYRAD-94 a vertical resolution of 5 m can be used (note that in principal, even finer vertical resolutions are possible), which is beneficial for evaluating high-resolution model simulations such as large-eddy simulations (LES) (e.g., van der Dussen et al. 2013) or for investigating processes with strong vertical gradients (e.g., melting layer).
The increased time on target, because of continuous transmission, provides similar average output powers as pulsed systems and leads to comparable sensitivities. In general, FMCW systems are cheaper to build because of their simpler electronic design and lack of high-voltage components (Ligthart and Nieuwkerk 1980), which further leads to better long-term stability and easier calibration. Also, their smaller instrument dimensions and lower transmitter power facilitate their installation in various environments.
High-frequency radars operating in the W-band provide greater sensitivity to small hydrometeors, which allow, for example, for investigating thin liquid or ice clouds and their microphysical processes (Lhermitte 1987). The disadvantage of using high frequencies is increasing signal attenuation as a result of dry gases, water vapor, and frozen and liquid hydrometeors.
Figure 1 shows the W-band radar radiometer at JOYCE (Löhnert et al. 2015) in Jülich, Germany, and at the Alfred Wegener Institute for Polar and Marine Research and the French Polar Institute Paul Emile Victor (AWIPEV) Arctic Research Base in Ny-Ålesund, Norway. At JOYCE (www.joyce.cloud) it is located next to a pulsed 35-GHz system (JOYRAD-35) for high-accuracy dual-frequency retrievals of liquid water content. The Arctic deployment (2016/2017) took place within the German research initiative Arctic Amplification: Climate Relevant Atmospheric and Surface Processes and Feedback Mechanisms [(AC)3] (http://www.ac3-tr.de/), where the objective is to detect, quantify, and model the feedback mechanism leading to the observed accelerated warming in the Arctic compared to global mean.
JOYRAD-94 at (a) JOYCE and (b) AWIPEV station in Ny-Ålesund.
Citation: Journal of Atmospheric and Oceanic Technology 34, 11; 10.1175/JTECH-D-17-0019.1
JOYRAD-94 is a very compact system with a horizontal extent of 1.0 × 1.1 m2, a height of 1.3 m, and is installed and recoverable with two persons. It operates under all weather conditions as a result of two high-powered blowing systems that nominally keep the radomes free of liquid and ice. A dry radome is important because a thin liquid layer can contaminate the radar and radiometer measurements significantly; for example, Hogan et al. (2003) observed a two-way radome attenuation of 15 dB. Also, a thin liquid layer on the receiver radome can affect the radiometer brightness temperature measurements by several kelvins, biasing the retrieval of LWP.
As will be explained in detail later in the manuscript, the vertical resolution of JOYRAD-94 is adjustable within up to 10 different layers, making it possible to adapt the resolution and sensitivity to different measurement objectives at different heights. JOYRAD-94 receives the passive and active signals over the same antenna—that is, both sensors observe the same scene and have the same beamwidth of 0.53°—which reduces the uncertainty of retrievals when combining radar and radiometer measurements (Frisch et al. 1998). Such a synergetic instrument was design by Skou (1995) for airborne observations; however, it has not yet been used for ground-based studies.
Common radiometers have beamwidths between 3° and 6°, and therefore smooth out horizontally varying cloud structures (Huang et al. 2014). The passive channel of JOYRAD-94 provides new opportunities for observing spatial and temporal cloud variability, for example, questioning Taylor’s hypothesis of “frozen” turbulence, which was found to be valid for radiometers with beamwidths larger than 1° (e.g., Huang et al. 2008). Moreover, combining matched radar–radiometer measurements with small beamwidths and high vertical resolution can add information on processes at cloud edges, such as turbulence and evaporation induced by lateral entrainment (Heus and Jonker 2008). However, the latter and LWP variability will not be part of the following discussion; but it is left as a topic for future research.
In the following, we will introduce the instrument’s design, which enables continuous and unattended observations while receiving both the active and passive signals with the same antenna. In section 3, the signal processing is discussed, including the adjustable measurement setup for different atmospheric layers. The performance of the instrument is shown in section 4 followed by a comparison to a 35-GHz pulsed radar. In section 6, we introduce the novelties of JOYRAD-94 for cloud observations. The last section gives a summary and an outlook considering software and hardware updates that will be or have already been included in the system.
2. Instrument design
The instrument is designed for long-term unattended observations of cloud macro- and microphysical properties by combining passive and active microwave radiometry at center frequencies of 89 and 94 GHz, respectively. FMCW radar is used to provide high flexibility of range and Doppler resolution variation along the ranging path. Moreover, low transmitter power (typically 1.5 W) allows for solid-state technology, leading to a long lifetime and low costs.
The main components of the radar–radiometer are the transmitter with adjustable power to protect the receiver from saturation, the Cassegrain two-antenna system, and the receiver containing both the radar receiver channel at 94 GHz and the passive broadband channel at 89 GHz. Both channels are thermally stabilized within a few millikelvin (mK) to ensure accurate output power measurement and radiometric accuracy. A block diagram of the full system and technical specifications are given in Fig. 2 and Table 1, respectively.
Block diagram of JOYRAD-94.
Citation: Journal of Atmospheric and Oceanic Technology 34, 11; 10.1175/JTECH-D-17-0019.1
Radar specifications.
a. Antenna
The Cassegrain antenna, which was recently measured using the solar sun scan technique (Reimann and Hagen 2016), produces a beam of 95% Gaussian shape. Moreover, the modeled and measured half-power beamwidths (HPWB) are 0.48° and 0.53°, respectively. The modeled and measured antenna gains are 51.8 and 51.5 dB, respectively. JOYRAD-94 was calibrated with the modeled values, which leads to an overestimation of reflectivity of 0.26 dB that has been accounted for in the data analysis. Future versions will use the measured parameters instead. The two antenna systems protect the receiver from saturation, as during continuous transmitter operation it is not possible to separate transmitted and received power. Transmitting and receiving antennas are separated by 568 mm, which leads to a signal loss of 10% at 300 m distance as a result of imperfect beam overlap. Because of the small antenna separation, backscattering theory is still fully applicable. The far field of the antenna begins at about 50 m. The antenna system losses are dominated by subreflector blocking—that is, the subreflector shading the output beam—which was found to cause a power loss of 11%.
b. Transmitter
The maximum output power is continuously monitored by the radar with a retraceable precision Rhode and Schwarz W-band power head (accuracy of 0.1 dB). The feedhorn insertion loss is 0.3 ± 0.05 dB, which was determined by a Rhode and Schwarz network analyzer (accuracy of 0.05 W). To ensure linear detection while observing a high dynamic range of input power signals (from thin cirrus to precipitating liquid clouds), the radar’s transmitter power can be reduced in digital steps by a maximum of 16 dB. Sawtooth chirps are generated by a chirp generator module based on a direct digital synthesizer (DDS) that is clocked by a fundamental oscillator of frequency of 916 MHz. The chirp generator produces an output signal in the range of 7.79–7.88 GHz that is then multiplied by 12 to become the transmitter signal between 94.48 and 94.56 GHz.
c. Receiver calibration



During the “hot–cold” calibration, the radiometer channel at 89 GHz is calibrated too. The 89-GHz channel has an RMS noise smaller than 0.5 K when integrating on the hot calibration target for any integration times between 5 and 100 s. Küchler et al. (2016) found an uncertainty of ±0.5 K for direct detection radiometers of a similar type that were calibrated with liquid nitrogen. Maschwitz et al. (2013) suggested repeating absolute calibrations every 3–6 months to sustain measurement accuracy.
d. Calibration cross-check
A calibration cross-check with a metal sphere reflector showed an agreement within 0.5 dB between the expected power ratio of the transmitted and received signals and the measured power ratio. The sphere had a diameter of 48 mm; was lifted by a helium balloon to 104 m above ground at a horizontal distance of approximately 950 m to the radar, resulting in a radial distance of 960 m; and was strapped to the ground by three strings. The distance between the balloon and the metal sphere was about 10 m. The radar was mounted on a positioner that was manually operated using a telescope for precise pointing. To account for the horizontal swinging of the sphere, the integration time was set to 0.2 s and a peak-hold function in the radar software was run, which samples the maximum value observed within 60 s. The procedure was repeated 100 times, resulting in a standard deviation of 0.15 dB of the return signal. Atmospheric attenuation was found to be 0.6 dB over a path of 1920 m, which was determined from radiative transfer calculations assuming homogeneity on the radiation’s path and using surface data from the radar’s meteorological weather station.
e. Blowers
A high-powered blower system was installed on JOYRAD-94 for both antennas (Fig. 1, white fixtures) to prevent the adherence of hydrometeors on the radome. The volume flux is 2000 m3 h−1, corresponding to approximately 20 m s−1 wind speed in full-power mode. The blowers are triggered by a Vaisala rain sensor that is part of the meteorological weather station attached to the radar (Fig. 1, yellow). Moreover, the radome material has a hydrophobic coating and is tilted so that condensate nominally runs off. Therefore, JOYRAD-94 can operate continuously during all weather conditions.
3. Signal processing
In this section we briefly describe only the main processing steps and refer the reader to the cited literature (e.g., Strauch 1976; Delanoë et al. 2016; Ulaby and Long 2014, 615–621) for details of FMCW signal processing.
a. FMCW principle
















Example of a sawtooth chirp
Citation: Journal of Atmospheric and Oceanic Technology 34, 11; 10.1175/JTECH-D-17-0019.1
b. Acquisition and processing
The IF signal from the mixer output is sampled by a fast ADC board, whereas the sampling clock of the ADC is also generated by the chirp generator to ensure an optimal synchronization of ramp generation and sampling. The ADC’s maximum sampling rate is 28.625 MHz. According to Parsivel’s theorem (Hughes 1965), the mean power of the voltage time series, being the sum of the squared sample voltages, equals the sum of the power spectrum over the sampling interval. Thus, applying a Fourier transform to the digitized mixer output voltage provides the power spectrum

















The radar software conducts “level 1” processing, that is, it determines the noise level of the spectra and calculates spectral moments. Moreover, visualization tools are available for instantaneous data monitoring. The detection threshold for a valid sample is a multiple of the noise density standard deviation above the mean noise density of the Doppler spectrum, which must be exceeded in the spectrum. For more details, we refer the reader to Görsdorf et al. [2015, parameter Q in Eq. (14) is the “noise filter”].
c. Doppler spectra
Example spectra are shown in Fig. 4 at three different heights of a precipitating cloud. Particles moving toward the radar have a negative velocity. The noise level was determined using the Hildebrand–Sekhon scheme (Hildebrand and Sekhon 1974). The top panel illustrates a Gaussian-shaped power spectrum of a sample volume that contained ice particles exhibiting a mean Doppler velocity of −1 m s−1. The second power spectrum was sampled within the melting layer showing a negative skewness, a broadening of the spectrum, and a shift to more negative velocities, which indicate the acceleration of melting ice particles. The third panel exhibits a broad spectrum that is typical for regions below the melting layer, where a mixture of liquid particles with different velocities was present.
Example Doppler spectra of JOYRAD-94 (top) above, (middle) within, and (bottom) below the melting layer. Negative velocities indicate motion toward the radar. The integration time was 0.5 s. The mean noise level is indicated (dashed line).
Citation: Journal of Atmospheric and Oceanic Technology 34, 11; 10.1175/JTECH-D-17-0019.1
The lhs of the third panel indicates that the sample volume contained particles whose fall velocities were close to the Nyquist velocity. At lower range gates, particles were present that exceeded the Nyquist velocity (not shown here). If particles move faster than
d. The chirp table
The current measurement modes that are operated on the radar run with a sampling rate of approximately 8 MHz. The Nyquist limit restricts the maximum IF, which can be sampled distinctly, to 4 MHz. Moreover, the chirp generator’s phase noise contaminates IFs below 300 kHz. Therefore, the usable IF signal band is limited to frequencies between 300 and 4 MHz. Because the IF bandwidth determines the maximum unambiguous range, a range from 0 to 12 km cannot be covered by a single chirp sequence. Therefore, several chirp sequences are executed consecutively, whereas for each chirp sequence the IFs between 300 kHz and 4 MHz span different ranges, for example, 100–400 m in the first chirp sequence, 400–2000 m in the second chirp sequence, etc. As a consequence, the final profile is a composite of several chirp sequences.
On the one hand, this is a disadvantage, since the data in different chirp sequences were recorded at different times. On the other hand, the sequence of different chirps provides the opportunity to define several parameters, such as
Currently, the radar is operated in two modes: the standard mode (SM) and the high-vertical-resolution mode (HRM). The former operates with vertical resolutions between 16 and 34 m, with Doppler velocity resolutions (
Main attributes of two radar operation programs: SM and HRM. Several parameters, such as range resolution,
To avoid second-trip echoes, the maximum unambiguous range must be sufficiently large. However, increasing the maximum unambiguous range decreases the maximum unambiguous velocity (see section 3a). On the one hand, this can be a disadvantage, as aliasing becomes likely in lower layers, though it can be corrected, as will be shown in section 6b. On the other hand, it is likely that the transmitted signal at 94 GHz will be attenuated entirely before reaching regions above 10 km, if a convective precipitating cloud system is present. Nevertheless, there is the opportunity to define a measurement mode with two different chirp tables being executed alternatingly to enable observations of both high clouds and precipitation.
4. Radar–radiometer performance
Note that the data presented in the following were collected over a time frame (October 2015–May 2017) in which the sensitivity of JOYRAD-94 was decreased by two reasons: (i) the transmitter power chip degraded rapidly because of a production failure until it was replaced in April 2016, which caused a decrease in sensitivity; and (ii) the detection threshold (see section 3b) was calculated wrongly in the first software version, by not considering coherent averaging, which led to an additional loss of sensitivity by about 3 dB. The software was updated in February 2017.
a. Case study
The performance of the radar–radiometer is evaluated using a 24-h measurement period of 5 November 2015, which is shown in Fig. 5. Several precipitation events passed JOYCE during this day. These were mainly triggered by the southwesterly inflow of warm and moist air over central Europe during the days before and an overpass of an elongated trough at 500 hPa from 4 to 6 November. The precipitation event in the morning was the tail of a convective system occurring the night before. The intensification of precipitation between 0100 and 0200 UTC is clearly visible by the increasing reflectivities in Fig. 5 (top).
The 24-h time series of JOYRAD-94 measurements at JOYCE on 5 Nov 2015. (top ) Radar reflectivity. Cloud-base height from a collocated ceilometer at 3-m distance (black dots). (bottom) Brightness Temperature
Citation: Journal of Atmospheric and Oceanic Technology 34, 11; 10.1175/JTECH-D-17-0019.1
The overlay of the ceilometer cloud-base height reveals a liquid cloud layer at the top of the cloud that was present between 0300 and 0500 UTC. It also shows the advantage of a radar in comparison to a ceilometer, as the radar captured ice particles below that were not detected by the ceilometer. JOYRAD-94 even detects clouds having a thickness of a few tens of meters; however, very thin liquid layers—for example, between 0800 and 0900 UTC—were beyond the sensitivity of the radar, mainly because of the too-conservative detection threshold in the first software version.
The second precipitating system crossed JOYCE between 2000 and 0200 UTC in night from 5 to 6 November. The transmitter power is automatically regulated to prevent the receiver from saturation when high reflecting particles, such as large raindrops, are present. Such events involve strong attenuation so that complete attenuation is possible at low altitudes, for example, as visible in Fig. 5 at 2230 UTC, where the reflectivity is significantly reduced after 2 km (yellow vertical stripe). Here, the transmitted power was decreased from 1.5 to 0.2 W between 2200 and 2300 UTC when precipitation produced an increased reflectivity signal. On the one hand, this causes a loss of information from higher cloud levels; on the other hand, it keeps the detector in its linear regime, thus, it increases data quality in lower levels.
The brightness temperature
The
In this section, we showed that the passive channel complements the radar measurements while observing the same scene. In section 6, we will further discuss the performance of LWP retrievals using
b. Sensitivity
The sensitivity of JOYRAD-94 is illustrated in Fig. 6a. It shows the occurrence frequencies of observed reflectivities at each range bin when the radar has operated in the SM at JOYCE. Only signals that exceeded the reflectivity threshold were taken into account. The figure reveals several features that correspond to certain cloud types or specific positions within the cloud.
(a) Histogram of
Citation: Journal of Atmospheric and Oceanic Technology 34, 11; 10.1175/JTECH-D-17-0019.1
The area around the maximum at about 8.5 km and −25 dBZ corresponds to high-altitude ice clouds that were often connected to precipitating systems (e.g., as shown in Fig. 5) that occurred frequently during the observation period between November 2015 and March 2016. The occurrence of precipitating systems is reflected in the area around 10 dBZ between 0- and 2-km height when precipitation caused strong echoes. While combining Cloudnet data and synoptic maps, we found that the region of increased occurrences, connecting the maximum at the ground and at about 8.5 km, is mainly caused by cumulonimbus clouds associated with frontal passages over western Germany. Further cloud types that occurred frequently during the observation period were nonprecipitating cumulus and stratocumulus containing hydrometeors in the liquid phase. The fingerprints of these cloud events are evident in the local maximum between 500- and 2000-m altitude, exhibiting typical values from −40 to −20 dBZ (e.g., Frisch et al. 1995; Wang and Geerts 2003). The distinct bimodal signature in this region was also observed by a collocated 35-GHz pulsed radar, which is an interesting meteorological feature; however, a detailed analysis of the data and the synoptic conditions used to explain the bimodality is not part of this study.
A gap is visible at the sensitivity limit of JOYRAD-94 at about 7.8 km. It was found that a signal originating from ADC board electronics contaminates the signal processing. The signal’s frequency is equal to the IF frequency that corresponds to the range gate at about 7.8 km. To remove this contamination, the intrinsic noise underground in each IF bin was determined and stored. During operation, the contamination’s amplitude is subtracted from the monitored signal, which causes the gap in Fig. 6a, indicating that the logged value of the contamination signal was overestimated.
Fitting a curve to the lowest detected reflectivities (Fig. 6a) provides the sensitivity of the radar for the SM at JOYCE after filtering the data with the detection threshold discussed in section 3b. Moreover, the sensitivity was determined for each month of operation individually, showing a shift of the sensitivity curve toward larger reflectivities with advancing time. The loss of sensitivity is due to the degradation of an erroneous transmitter power chip implemented in the prototype. The transmitter power chip was replaced in April 2016 and since then the transmitter power has been stable. Note that the curve of the HRM is also influenced by sensitivity losses caused by a higher vertical resolution.
In June 2016, the radar was installed in Ny-Ålesund (NyA), where the radar operated in the HRM and the SM. A software update in February 2017, with improved signal detection thresholds, further increased the sensitivity by 3 dB. The comparison of the resulting sensitivity profiles before and after the replacement, and the software update in Fig. 6a reveals that the sensitivity in the SM has improved with respect to the early time of installation at JOYCE when the power chip was still working properly. We find that the sensitivity loss in the HRM observed at JOYCE is mainly a result of the chip degradation and a too conservative detection threshold; the differences of the HRM and the SM with the new chip are smaller than 6 dB in the lowest 500 m and smaller than 4 dB above.
As clearly visible in Fig. 6a, jumps in the sensitivity occur at several heights. These are caused by different chirp sequences used for probing different altitude regions. In general, a finer vertical resolution leads to a lower sensitivity and vice versa. Because the vertical resolution often changes from chirp to chirp, the sensitivity changes accordingly. The sensitivity limits for all layers in both operation modes are summarized in Table 2.
We further assess the performance of JOYRAD-94 by comparing it to a pulsed 35-GHz radar (JOYRAD-35), as described by Görsdorf et al. (2015), which was located at 3 m distance to JOYRAD-94 at JOYCE. The sensitivity curve of JOYRAD-35 is also included in Fig. 6 as a reference. Clearly evident is the higher sensitivity of JOYRAD-35 in the lowest 3 km (Fig. 6b), which is due to the much higher average transmitter power of 24 W (Acquistapace et al. 2017, their Table 2) and the short integration times in the HRM and SM in the lowest layer. These differences in sensitivity are on the order of 10 dB (15 dB) in the lowest kilometer for the SM (HRM) but reduce to 3 dB (8 dB) above 3 km. These sensitivity differences also affect the analysis in the following section. Note that increasing the integration time of JOYRAD-94 by 0.5 s in the lowest layer would increase the sensitivity by about 3 dB.
5. Assessment of cloud detection
The following analysis is based on the Cloudnet product of Illingworth et al. (2007), combining several remote sensing instruments (e.g., radar, lidar, and MWR) and reanalysis data to classify cloud properties such as vertical thickness, phase, and hydrometeor type. Cloudnet provides a target classification with a vertical resolution of 30 m, allowing for checking which hydrometeor types are present in the observed atmospheric column. We compared the performances of JOYRAD-94 and JOYRAD-35 for two cloud types: (i) single-layer nonprecipitating liquid clouds with a cloud base higher than 200 m and (ii) single-layer nonprecipitating ice clouds with a cloud base higher than 200 m. These cloud types were selected to ensure that both radars operate in the Rayleigh regime so that a comparison of reflectivities is reasonable.
The minimum cloud base was selected to ensure that the clouds were observed in the far field of both radars. A cloud was considered to be precipitating if either drizzle, rain, or melting ice particles, or a combination of these three, were connected to the cloud or were detected within the cloud. Falling ice particles were not defined as precipitation (Illingworth et al. 2007).
The Cloudnet product is derived with 30-s temporal and 30-m vertical resolution. We considered only data that were recorded by both radars simultaneously and within ±10 s of a Cloudnet data point while either a single-layer liquid cloud or a single-layer ice cloud was present. The vertical resolution of JOYRAD-35 was 29 m. JOYRAD-94 was operated in the SM; that is, the resolution varied between 16 and 34 m (see Table 2). Height matching was performed using nearest neighbor comparison. The integration time of JOYRAD-35 was 2 s. The integration time of JOYRAD-94 varied with height between 0.34 and 1.77 s (see Table 2). No corrections concerning gaseous and liquid or ice attenuation were applied to the data, since neither measurements of humidity and temperature profiles nor measurements of integrated water vapor (IWV) are available for this period at JOYCE.
a. Liquid clouds
After selecting single-layer liquid clouds using the Cloudnet classification scheme (a total of 41 791 cases), we investigated the performance of both radars at the cloud center. The cloud center was chosen to exclude the effects of partial beamfilling at cloud edges. The position of the cloud center was calculated from the Cloudnet cloud base and cloud top. Only in 73% of the cases was a signal simultaneously detected by both radars. This is due to the higher sensitivity of JOYRAD-35 compared to JOYRAD-94, especially because of transmitter chip degradation in JOYRAD-94 during the initial deployment at JOYCE and the too conservative detection threshold. If JOYRAD-94 had observed continuously with its nominal reflectivity (see Fig. 6a, “NyA SM”), liquid clouds would have been detected simultaneously in 97% of the cases.
Figure 7 shows scatter histograms of the radar reflectivity and the mean Doppler velocity for simultaneously detected single-layer liquid clouds. Most cases reveal reflectivities between −40 and −30 dBZ, corresponding to the area of enhanced occurrence between 1 and 2 km in Fig. 6a. Both radars agree with a mean difference of 0.52 dB and a standard deviation of 1.82 dB. The mean Doppler velocities agree well, showing a mean difference of 0.04 m s−1 with a standard deviation of 0.34 m s−1.
Scatter histograms of (a) reflectivities and (b) mean Doppler velocities of JOYRAD-94 and JOYRAD-35, which were observed at the cloud center of single-layer liquid clouds having a cloud base above 200 m.
Citation: Journal of Atmospheric and Oceanic Technology 34, 11; 10.1175/JTECH-D-17-0019.1
b. Ice clouds
We classified 49 556 cases as single-layer ice clouds; 93% of the cases were observed simultaneously by JOYRAD-35 and JOYRAD-94. The higher percentage with respect to liquid clouds is due to higher reflectivities that are associated with these clouds, as can be seen in Fig. 8a; that is, the sensitivity limit of JOYRAD-94 was mostly exceeded. If JOYRAD-94 had measured with its nominal reflectivity, ice clouds would have been detected simultaneously in 99% of the cases.
Scatter histograms of (a) reflectivities and (b) mean Doppler velocities of JOYRAD-94 and JOYRAD-35, which were observed at the cloud center of single-layer ice clouds having a cloud base above 200 m.
Citation: Journal of Atmospheric and Oceanic Technology 34, 11; 10.1175/JTECH-D-17-0019.1
The majority of reflectivities are found between −30 and −20 dBZ, which corresponds to the area of enhanced occurrence in Fig. 6a between 6 and 8 km. The values between −20 and −5 dBZ predominantly originate from the area in Fig. 6a connecting the maximum at 8.5 km and at the ground. The distribution in the scatterplot in Fig. 8a starts to deviate from the 1:1 line for
The mean Doppler velocities are concentrated between −1 and 0 m s−1, which are characteristic for ice particles in high altitudes (e.g., Heymsfield and Iaquinta 2000). Again, a positive offset is evident, but still both radars agree with a mean difference of 0.08 ± 0.15 m s−1.
c. Comparative uncertainty estimate
The velocity offset has the same sign for liquid and ice clouds. This indicates that the velocity offset is due to a relative mispointing, that is, either one of the radars or both were not exactly aligned to zenith. Kneifel et al. (2016) found that a slight relative mispointing of 1° elevation produces a velocity offset of 0.1 m s−1, which is already larger than the offset found here.
Since corrections for gaseous and liquid attenuation could not be applied, a quantitative estimate for calibration accuracies cannot be drawn from the comparison of the radars’ reflectivities. All attenuating effects—that is, liquid and gaseous attenuation effects, as well as differential scattering at large particles—cause positive reflectivity offsets, that is,
When comparing Fig. 8a to Fig. 7a, the differential attenuation becomes larger by about 0.6 dB; that is, the offset turns positive. On the one hand, this offset could be explained by attenuation by liquid water. At JOYCE, nonprecipitating single-layer liquid clouds exhibit averaged LWC values of about 0.1 g m−3 and smaller while having vertical extents of about 500–1000 m (Löhnert et al. 2015). This induces a two-way attenuation from the cloud base to the cloud center of about 0.4–0.7 dB (Hogan et al. 2005). On the other hand, it is possible that the average magnitude of the gaseous attenuation differed between the liquid cloud and ice cloud cases. Because the mean altitudes of the cloud centers were 1.2 and 6.1 km for liquid and ice clouds, respectively, we assume that gaseous attenuation has been weaker during liquid cloud events, which could compensate for the differential attenuation by liquid water. The differential two-way attenuation by water vapor is about 1 dB for the lowest 5 km for values of about 10 kg m−2 IWV (Kneifel et al. 2015), being typical for fall and winter at JOYCE (Löhnert et al. 2015). Further, attenuation by oxygen can be up to 0.5 dB for the entire atmosphere. We estimate the decrease in reflectivity difference from Fig. 8a to Fig. 7a to be smaller than 1 dB, though sufficient to compensate attenuation effects by liquid water; that is, we would expect a similar value for
6. New opportunities for cloud observations
This section presents the novelties for cloud observations that arise from JOYRAD-94 that have optimally matched active and passive components. Moreover, we show that it is possible to overcome the trade-off between Nyquist velocity and range resolution without losing temporal resolution when combining data from two radars with different settings.
a. Liquid water path retrieval at 89 GHz
Optimal beam matching between a radar and a radiometer provides instantaneous signatures as was shown in section 4. To make use of these signatures, retrievals are necessary that add information to the radar measurement. Frisch et al. (1998) demonstrated that the vertical distribution of liquid water can be obtained when combining the LWP from a radiometer and radar reflectivity measurements. The liquid water path can be derived from a microwave radiometer at window frequencies, where clear-sky attenuation is low, for example, at 31 or 89 GHz. However, because of the water vapor absorption continuum in the microwave (Rosenkranz 1998), additional information on IWV is needed to constrain retrievals.






Figure 9 shows the performance of an LWP retrieval using only the 89-GHz channel as the input parameter in Eq. (9) (Fig. 9a) and using the 89-GHz channel with additional information on the integrated water vapor from an external source with a random uncertainty of ±2 kg m−2 (e.g., short-term model forecast; Fig. 9b). We assumed an uncertainty of 0.5 K for
(a) LWP retrieval using only 89 GHz deviates on average by 118 ± 44 g m−2 from the model truth (x axis). (b) LWP retrieval using 89 GHz plus IWV from an external source with ±2 kg m−2 accuracy (“noisy IWV”) deviates on average by 118 ± 15 g m−2 from the model truth. Retrievals were performed with 15 175 radiosonde profiles to which clouds were added using a modified adiabatic approach.
Citation: Journal of Atmospheric and Oceanic Technology 34, 11; 10.1175/JTECH-D-17-0019.1
For comparison we ran the retrieval algorithm using only seven frequencies between 20 and 31 GHz, which are commonly used for LWP (e.g., Löhnert and Crewell 2003). The RMS uncertainty of this retrieval was 25 g m−2, being 60% larger than the RMS uncertainty of the retrieval combining 89 GHz and the IWV from an external source. This is due to the 89-GHz channel being more sensitive to liquid water than the 31-GHz channel. Combining
b. Dual-radar dealiasing of Doppler spectra
On the one hand, operating JOYRAD-94 in the HRM provides vertical highly resolved data; on the other hand, the Nyquist velocity is small; hence, aliasing is likely when observing strongly turbulent clouds or precipitation. Maahn and Kollias (2012) introduced a method to dealias Doppler spectra recorded by a rain FMCW radar operating at 24 GHz. A stand-alone pulsed radar must vary the pulse repetition frequency to obtain the real atmospheric signal, which decreases the effective observation time (e.g., Holleman and Beekhuis 2003; Sosnytskiy 2014). In FMCW radar data, aliasing produces a folding into the upper (lower) range gate for particles having an absolute radial velocity larger than
Using two radars that operate simultaneously, having opposite settings with respect to
Flowchart of the DRDM using
Citation: Journal of Atmospheric and Oceanic Technology 34, 11; 10.1175/JTECH-D-17-0019.1
Figure 11 shows two comparisons of mean Doppler velocities of JOYRAD-94 (
Scatter histograms of mean Doppler velocities of JOYRAD-94 (
Citation: Journal of Atmospheric and Oceanic Technology 34, 11; 10.1175/JTECH-D-17-0019.1
The DRDM is simple and computationally inexpensive. Moreover, a prior evaluation of signal quality can further decrease the computing time (Maahn and Kollias 2012). The disadvantages are clearly that bimodal spectra with fully separated peaks could be cut when the third step of the DRDM is conducted. This can be compensated for when including a further processing step that looks for secondary peaks in the spectra. Moreover, if the spectrum is broader than the
The choice of the maximum unambiguous velocities (
To get a copy of the MATLAB code of the DRDM, please contact the corresponding author of this manuscript.
7. Summary and outlook
We presented a new 94-GHz FMCW cloud radar–radiometer (JOYRAD-94) suited for studying cloud microphysical processes and cloud macrophysical properties. JOYRAD-94 is a compact instrument for unattended long-term measurements in all weather conditions. A strong blowing system minimizes adherence of liquid and ice particles on the radome and therefore attenuation effects are avoided. Additionally, the transmitter power can be regulated to prevent the detector from saturation and therefore provide accurate measurements also when strongly reflecting particles are present in the observed column, for example, during precipitation. A further novelty is that the active and passive components receive over the same antenna; that is, optimal beam matching is accomplished.
The active component of JOYRAD-94 can be accurately calibrated with an uncertainty of ±0.5 dBZ. This was shown by observing a metal sphere in the far field of the radar and by a relative comparison to a collocated 35-GHz radar. The calibration of the receiver should be repeated every 3–6 months. The passive broadband channel at 89 GHz has an uncertainty of about ±0.5 K. JOYRAD-94 was tested in two different measurement configurations: a standard configuration (SM) with vertical resolutions between 16 and 34 m and a high-vertical-resolution mode (HRM) resolving the lowest 3 km with about 6-m range resolution. In both modes JOYRAD-94 is sensitive enough to capture the majority of clouds. The sensitivity ranges from −67 (−61) dBZ at 100 m to −38 (−35) dBZ at 10 km when operating in the SM (HRM). The HRM particularly allows for improving observations of small-scale microphysical processes (e.g., processes affected by turbulence) and processes comprising a strong vertical gradient such as melting layer or riming within layers of supercooled liquid water.
It is also possible to adjust the measurement mode within up to 10 layers in the vertical column to account for different cloud particle properties at different heights. For example, it can be beneficial to run a coarse range resolution with a large unambiguous
It was shown that the passive channel can provide accurate LWP estimates. A model study revealed that the liquid water path can be retrieved with an RMS uncertainty of about 15 g m−2 when measuring brightness temperatures at 89 GHz and knowing the IWV from an external source with an uncertainty of ±2 kg m−2. In addition, systematic offsets—for example, as a result of calibration offsets—have to be taken into account. The information of the 89-GHz channel can be used to identify liquid constituents within the vertical column and to retrieve vertical profiles of liquid water when including radar measurements. The impact of optimally matched beams on the retrieval quality is content of a current study.
We presented a new method to dealias Doppler spectra (DRDM) without losing temporal resolution by using the data of the collocated radar as the initial guess. The DRDM provides vertically highly resolved spectra with an increased
As JOYRAD-94 is a prototype, future versions of this instrument will have to overcome initial teething troubles: the first version of the transmitter power chip was degrading very fast, which caused a significant loss of sensitivity of 10 dB. However, this problem has already been solved and we showed that the new power chip produces a stable output power. Furthermore, the latest version of the manufacturer’s software includes an update of the threshold detection procedure, the optimization of the suppression of the ADC board electronics’ signal that contaminates the signal processing, and a dealiasing routine to provide accurate spectra. Next generations will also have the possibility of polarimetric detection and can be mounted on a scanner to derive three-dimensional cloud fields. Moreover, a frequent switching between the two measurement modes will be enabled so that the DRDM can also be applied to single radar measurements when assuming constant conditions within two data samples.
Acknowledgments
This work was funded by the project “Energy Transitions and Climate Change” (http://et-cc.uni-koeln.de/project.html) of the “Excellence Initiative” of the University of Cologne by the German Research Foundation (DFG) under Grant ZUK 81/1 and by the German research initiative (AC)3 by the DFG under Grant TRR172/1 (http://www.ac3-tr.de/). Contributions by S. Kneifel were carried out within the Emmy-Noether Group OPTIMIce funded by the DFG under Grant KN 1112/2-1. We thank Rainer Haseneder-Lind, Pavel Krobot, Birger Bohn, Alexander Myagkov, and the AWIPEV team in Ny-Ålesund for the technical support; and in particular the three anonymous reviewers for their constructive comments, which helped to improve the manuscript.
REFERENCES
Acquistapace, C., S. Kneifel, U. Löhnert, P. Kollias, M. Maahn, and M. R. Bauer-Pfunstein, 2017: Optimizing observations of drizzle onset with millimeter-wavelength radars. Atmos. Meas. Tech., 10, 1783–1802, doi:10.5194/amt-10-1783-2017.
Bennett, A., J. Nash, C. Gaffard, B. Moyna, M. Oldfield, and P. Huggard, 2009: Observations from the UK Met Office 94 GHz FMCW cloud radar. Proceedings of the Eighth International Symposium on Tropospheric Profiling, A. Apituley, H. W. J. Russchenberg, and W. A. A. Monna, Eds., Organizing Committee of the 8th International Symposium on Tropospheric Profiling, S06-O02, 4 pp., http://www.cesar-observatory.nl/istp8/pages/index.htm.
Borque, P., E. Luke, and P. Kollias, 2016: On the unified estimation of turbulence eddy dissipation rate using Doppler cloud radars and lidars. J. Geophys. Res. Atmos., 120, 5972–5989, doi:10.1002/2015JD024543.
Boucher, O., and Coauthors, 2013: Clouds and aerosols. Climate Change 2013: The Physical Basis, T. F. Stocker et al., Eds., Cambridge University Press, 571–657.
Delanoë, J., and Coauthors, 2016: BASTA: A 95-GHz FMCW Doppler radar for cloud and fog studies. J. Atmos. Oceanic Technol., 33, 1023–1038, doi:10.1175/JTECH-D-15-0104.1.
Frisch, A. S., C. W. Fairall, and J. B. Snider, 1995: Measurement of stratus cloud and drizzle parameters in ASTEX with a Kα-band Doppler radar and a microwave radiometer. J. Atmos. Sci., 52, 2788–2799, doi:10.1175/1520-0469(1995)052<2788:MOSCAD>2.0.CO;2.
Frisch, A. S., G. Feingold, C. W. Fairall, T. Uttal, and J. B. Snider, 1998: On cloud radar and microwave radiometer measurements of stratus cloud liquid water profiles. J. Geophys. Res., 103, 23 195–23 197, doi:10.1029/98JD01827.
Görsdorf, U., V. Lehmann, M. Bauer-Pfundstein, G. Peters, D. Vavriv, V. Vinogradov, and V. Volkov, 2015: A 35-GHz polarimetric doppler radar for long-term observations of cloud parameters—Description of system and data processing. J. Atmos. Oceanic Technol., 32, 675–690, doi:10.1175/JTECH-D-14-00066.1.
Haeffelin, M., L. Barth, and J. Delano, 2005: SIRTA, a ground-based atmospheric observatory for cloud and aerosol research. Ann. Geophys., 23, 253–275, doi:10.5194/angeo-23-253-2005.
Heijnen, S. H., L. P. Ligthart, and H. W. Russchenberg, 2000: First measurements with TARA; An S-band transportable atmospheric radar. Phys. Chem. Earth, 25B, 995–998, doi:10.1016/S1464-1909(00)00140-4.
Heus, T., and H. J. J. Jonker, 2008: Subsiding shells around shallow cumulus clouds. J. Atmos. Sci., 65, 1003–1018, doi:10.1175/2007JAS2322.1.
Heymsfield, A. J., and J. Iaquinta, 2000: Cirrus crystal terminal velocities. J. Atmos. Sci., 57, 916–938, doi:10.1175/1520-0469(2000)057<0916:CCTV>2.0.CO;2.
Hildebrand, P. H., and R. S. Sekhon, 1974: Objective determination of the noise level in Doppler spectra. J. Atmos. Sci., 13, 808–811, doi:10.1175/1520-0450(1974)013<0808:ODOTNL>2.0.CO;2
Hogan, R. J., D. Bouniol, D. N. Ladd, E. J. O’Connor, and A. J. Illingworth, 2003: Absolute calibration of 94/95-GHz radars using rain. J. Atmos. Oceanic Technol., 20, 572–580, doi:10.1175/1520-0426(2003)20<572:ACOGRU>2.0.CO;2.
Hogan, R. J., N. Gaussiat, and A. J. Illingworth, 2005: Stratocumulus liquid water content from dual-wavelength radar. J. Atmos. Oceanic Technol., 22, 1207–1218, doi:10.1175/JTECH1768.1.
Holleman, I., and H. Beekhuis, 2003: Analysis and correction of dual PRF velocity data. J. Atmos. Oceanic Technol., 20, 443–453, doi:10.1175/1520-0426(2003)20<443:AACODP>2.0.CO;2.
Huang, D., Y. Liu, and W. Wiscombe, 2008: Determination of cloud liquid water distribution using 3D cloud tomography. J. Geophys. Res., 113, D13201, doi:10.1029/2007JD009133.
Huang, D., E. Campos, and Y. Liu, 2014: Statistical characteristics of cloud variability. Part 1: Retrieved cloud liquid water path at three ARM sites. J. Geophys. Res. Atmos., 119, 10 813–10 828, doi:10.1002/2014JD022001.
Huggard, P. G., and Coauthors, 2008: 94 GHz FMCW cloud radar. Millimetre Wave and Terahertz Sensors and Technology, K. A. Krapels and N. A. Salmon, Eds., Society of Photo-Optical Instrumentation Engineers (SPIE Proceedings, Vol. 7117), 711704, doi:10.1117/12.800347.
Hughes, H. K., 1965: The physical meaning of Parseval’s theorem. Amer. J. Phys., 33, 99, doi:10.1119/1.1971337.
Illingworth, A. J., and Coauthors, 2007: Cloudnet: Continuous evaluation of cloud profiles in seven operational models using ground-based observations. Bull. Amer. Meteor. Soc., 88, 883–898, doi:10.1175/BAMS-88-6-883.
Illingworth, A. J., D. Cimini, C. Gaffard, M. Haeffelin, V. Lehmann, U. Löhnert, E. J. O’Connor, and D. Ruffieux, 2015: Exploiting existing ground-based remote sensing networks to improve high resolution weather forecasts. Bull. Amer. Meteor. Soc., 96, 2107–2125, doi:10.1175/BAMS-D-13-00283.1.
Kneifel, S., M. S. Kulie, and R. Bennartz, 2011: A triple-frequency approach to retrieve microphysical snowfall parameters. J. Geophys. Res., 116, D11203, doi:10.1029/2010JD015430.
Kneifel, S., A. von Lerber, J. Tiira, D. Moisseev, K. Kollias, and J. Leinonen, 2015: Observed relations between snowfall microphysics and triple-frequency radar measurements. J. Geophys. Res. Atmos., 120, 6034–6055, doi:10.1002/2015JD023156.
Kneifel, S., P. Kollias, A. Battaglia, J. Leinonen, M. Maahn, H. Kalesse, and F. Tridon, 2016: First observations of triple-frequency radar Doppler spectra in snowfall: Interpretation and applications. Geophys. Res. Lett., 43, 2225–2233, doi:10.1002/2015GL067618.
Kollias, P., B. A. Albrecht, and F. Marks, 2002: Why Mie? Bull. Amer. Meteor. Soc., 83, 1471–1483, doi:10.1175/BAMS-83-10-1471.
Kollias, P., E. E. Clothiaux, M. A. Miller, B. A. Albrecht, G. L. Stephens, and T. P. Ackerman, 2007a: Millimeter-wavelength radars: New frontier in atmospheric cloud and precipitation research. Bull. Amer. Meteor. Soc., 88, 1608–1624, doi:10.1175/BAMS-88-10-1608.
Kollias, P., E. E. Clothiaux, M. A. Miller, E. P. Luke, K. L. Johnson, K. P. Moran, K. B. Widener, and B. A. Albrecht, 2007b: The atmospheric radiation measurement program cloud profiling radars: Second-generation sampling strategies, processing, and cloud data products. J. Atmos. Oceanic Technol., 24, 1199–1214, doi:10.1175/JTECH2033.1.
Küchler, N., D. D. Turner, U. Löhnert, and S. Crewell, 2016: Calibrating ground-based microwave radiometers: Uncertainty and drifts. Radio Sci., 51, 311–327, doi:10.1002/2015RS005826.
Lhermitte, R. M., 1987: A 94-GHz Doppler radar for cloud observations. J. Atmos. Oceanic Technol., 4, 36–48, doi:10.1175/1520-0426(1987)004<0036:AGDRFC>2.0.CO;2.
Ligthart, L. P., and L. R. Nieuwkerk, 1980: F.M.-C.W. Delft atmospheric research radar. IEE Proc., 127F, 421–426, doi:10.1049/ip-f-1:19800061.
Löhnert, U., and S. Crewell, 2003: Accuracy of cloud liquid water path from ground-based microwave radiometry 1. Dependency on cloud model statistics. Radio Sci., 38, 8041, doi:10.1029/2002RS002654.
Löhnert, U., and Coauthors, 2015: JOYCE: Jülich Observatory for Cloud Evolution. Bull. Amer. Meteor. Soc., 96, 1157–1174, doi:10.1175/BAMS-D-14-00105.1.
Luke, E. P., and P. Kollias, 2013: Separating cloud and drizzle radar moments during precipitation onset using Doppler spectra. J. Atmos. Oceanic Technol., 30, 1656–1671, doi:10.1175/JTECH-D-11-00195.1.
Maahn, M., and P. Kollias, 2012: Improved Micro Rain Radar snow measurements using Doppler spectra post-processing. Atmos. Meas. Tech., 5, 2661–2673, doi:10.5194/amt-5-2661-2012.
Maahn, M., U. Löhnert, P. Kollias, R. C. Jackson, and G. M. McFarquhar, 2015: Developing and evaluating ice cloud parameterizations for forward modeling of radar moments using in situ aircraft observations. J. Atmos. Oceanic Technol., 32, 880–903, doi:10.1175/JTECH-D-14-00112.1.
Maschwitz, G., U. Löhnert, S. Crewell, T. Rose, and D. D. Turner, 2013: Investigation of ground-based microwave radiometer calibration techniques at 530 hPa. Atmos. Meas. Tech., 6, 2641–2658, doi:10.5194/amt-6-2641-2013.
Mead, J. B., I. PopStefanija, P. Kollias, B. A. Albrecht, and R. Bluth, 2003: Compact airborne solid-state 95 GHz FMCW radar system. 31st Int. Conf. on Radar Meteorology, Seattle, WA, Amer. Meteor. Soc., 4A.3, https://ams.confex.com/ams/32BC31R5C/webprogram/Paper63494.html.
Pazmany, A. L., R. E. McIntosh, R. D. Kelly, and G. Vali, 1994: A 94-GHz Doppler radar for cloud observations. IEEE Trans. Geosci. Remote Sens., 32, 731–739, doi:10.1109/36.298002.
Reimann, J., and M. Hagen, 2016: Antenna pattern measurements of weather radars using the sun and a point source. J. Atmos. Oceanic Technol., 33, 891–898, doi:10.1175/JTECH-D-15-0185.1.
Rosenkranz, P. W., 1998: Water vapor microwave continuum absorption: A comparison of measurements and models. Radio Sci., 33, 919–928, doi:10.1029/98RS01182.
Skou, N., 1995: A Scanning Microwave Radar and Radiometer. Conference Proceedings: Second Topical Symposium on Combined Optical-Microwave Earth and Atmosphere Sensing, IEEE, 215–216, doi:10.1109/COMEAS.1995.472318.
Sosnytskiy, S. V., 2014: Restoration of aliased average Doppler spectra. 2014 20th International Conference on Microwaves, Radar and Wireless Communications; MIKON-2014, A. Bikonis and K. Nyka, Eds., IEEE, 4 pp., doi:10.1109/MIKON.2014.6899960.
Strauch, R. G., 1976: Theory and application of the FM-CW Doppler radar. Ph.D. thesis, University of Colorado Boulder, 24 pp.
Thies, B., K. Müller, F. Maier, and J. Bendix, 2010: Fog monitoring using a new 94 GHz FMCW cloud radar. Proc. Fifth Int. Conf. on Fog, Fog Collection and Dew, Münster, Germany, Copernicus Gesellschaft mbH, FOGDEW2010-103, http://meetingorganizer.copernicus.org/FOGDEW2010/FOGDEW2010-103.pdf.
Tridon, F., and A. Battaglia, 2015: Dual-frequency radar Doppler spectral retrieval of rain drop size distributions and entangled dynamics variables. J. Geophys. Res. Atmos., 120, 5585–5601, doi:10.1002/2014JD023023.
Turner, D. D., M. P. Cadeddu, U. Löhnert, S. Crewell, and A. M. Vogelmann, 2009: Modifications to the water vapor continuum in the microwave suggested by ground-based 150-GHz observations. IEEE Trans. Geosci. Remote Sens., 47, 3326–3337, doi:10.1109/TGRS.2009.2022262.
Ulaby, F. T., and D. G. Long, Eds., 2014: Microwave Radar and Radiometric Remote Sensing. 1st ed. University of Michigan Press, 984 pp.
van der Dussen, J. J., and Coauthors, 2013: The GASS/EUCLIPSE model intercomparison of the stratocumulus transition as observed during ASTEX: LES results. J. Adv. Model. Earth Syst., 5, 483–499, doi:10.1002/jame.20033.
Wang, J., and B. Geerts, 2003: Identifying drizzle within marine stratus with W-band radar reflectivity. Atmos. Res., 69, 1–27, doi:10.1016/j.atmosres.2003.08.001.
Yamaguchi, J., and Coauthors, 2009: Evaluation of radar reflectivity (Z) for FMCW millimeter-wave cloud radar “FALCON-I.” IEEJ Trans. Fundam. Mater., 129, 183–189, 10.1541/ieejfms.129.183.