Low-level stratiform clouds have extensive areal coverage over the subtropical oceans, especially near the western continental boundaries (Wood, 2012). These shallow cloud layers moderate the ocean–atmosphere radiative exchange, exerting a strong impact on the sea surface temperature and the cloud–climate feedback (Brient et al., 2016; Zelinka et al., 2016). Variability in the response of these clouds to anthropogenic warming remains a significant source of spread in climate change projections (Zelinka et al., 2020). Due to their extensive coverage and strong interaction with radiation, shallow clouds can significantly affect the weather on seasonal and shorter timescales (Fast et al., 2019). Understanding the thermodynamic behavior and structure of shallow clouds is therefore of major importance to improve the representation of the ocean–atmosphere interaction and generate more robust weather forecast models and climate predictions (Bony and Dufresne, 2005; Zelinka et al., 2020).

Given their characteristic thickness, shallow clouds are challenging to accurately represent in global-scale simulations. Incorporating additional experimental data can be of use to develop more appropriate parameterizations for climate models (Randall, 2013; Medeiros et al., 2023; Schneider et al., 2024). A key property that, in large part, describes the radiative impact of shallow clouds is their liquid water content (LWC). The vertical profile of the LWC provides valuable information about the cloud's internal structure, including the turbulent mixing within the cloud layer. Accurate retrievals of LWC in shallow clouds, however, remain a complex endeavor (Ebell et al., 2010). The usual small amounts of liquid water in shallow clouds (< 1 gm-3) are hard to measure and require very sensitive instruments and an accurate methodology.

Radiometer and single-frequency radar observations have been used in the past to estimate the LWC by relying on assumptions of the vertical LWC profile, the drop size distribution, or accurate absolute instrument calibration (Frisch et al., 1998; Löhnert et al., 2001; Matrosov et al., 2004; Illingworth et al., 2007; Küchler et al., 2018). A more direct technique consists of using a dual-frequency radar system (Atlas, 1954; Eccles and Mueller, 1971; Hogan et al., 2005; Ellis and Vivekanandan, 2011). The estimation of the LWC profile using the differential frequency radar reflectivity (DFR) technique does not depend on the assumptions made in single-frequency radar approaches. If gaseous attenuation is accounted for, the measured differential attenuation signal from a pair of radar frequencies is proportional to the amount of liquid water in the shallow cloud if drop size effects are negligible (Lhermitte, 1990). In addition, the signal attenuation due to liquid droplets becomes stronger as the operating frequency increases. While the differential technique is not affected by systematic errors, random measurement errors are inversely proportional to the signal frequency. Both facts suggest the use of high frequencies and large frequency pair separations to measure greater differential attenuation and facilitate the retrieval of small LWC values with improved accuracy (Hogan et al., 2005; Williams and Vivekanandan, 2007). This has been demonstrated over the last few decades, as the LWC retrieval accuracy has progressively improved when higher-frequency radars became available, with the most recent works using Ka-band (35 GHz) and W-band (94 GHz) frequencies (Hogan et al., 2005; Zhu et al., 2019).

G-band radar instruments represent a recently developed capability in atmospheric remote sensing research, offering new possibilities to probe clouds and precipitation and leading to novel measurement techniques (Battaglia et al., 2014, 2020). They have been recently used to profile water vapor content in clouds (Cooper et al., 2021) and to provide novel insight into the submillimeter drop size distribution through Doppler measurements (Courtier et al., 2022; Socuellamos et al., 2024a) and multifrequency radar observations (Lamer et al., 2021; Socuellamos et al., 2024b). The synergy of a G-band and a Ka-band radar brings several distinct advantages to the LWC retrieval compared to the W-band and Ka-band pair: (1) increased differential attenuation for submillimeter drop sizes, (2) reduction in the reflectivity measurement random error, and, consequently, (3) improved sensitivity in the LWC determination.

The US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF) deployment in the context of the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE; Russell et al., 2021) offered a convenient opportunity to study shallow clouds through the deployment of several instruments at the Scripps Pier in La Jolla, CA, USA. A ground-based prototype radar called CloudCube participated in EPCAPE for a 6-week period in March and April 2023 (Socuellamos et al., 2024b). CloudCube is a multifrequency (Ka-, W-, and G-band) profiling radar developed at the Jet Propulsion Laboratory (JPL) under the National Aeronautics and Space Administration Earth Science Technology Office (NASA-ESTO) Instrument Incubator Program (IIP) that includes a 239 GHz (G-band) Doppler module (Socuellamos et al., 2024a). For the analysis developed in this work, the Ka-band and W-band channels of CloudCube did not have the sensitivity needed to detect the cloud formations of interest. Therefore, we have complemented CloudCube's G-band data with observations from the ARM Ka-band Zenith-pointing Radar (KAZR; Kollias et al., 2016), as well as ceilometer (Morris, 2016), microwave radiometer (Cadeddu, 2024), and radiosonde measurements. Using the Ka-band and G-band DFR observations, more than 100 profiles of LWC, averaged over 60 s, have been retrieved for two cases on 2 different days.

The article starts with a brief description of the radar instruments and the data selected for the analysis (Sect. ). This is followed in Sect.  by an explanation of the methodology developed to retrieve the LWC and a discussion of the results. Section  presents a reflectivity-based approach to retrieve the LWC. Vertically integrated LWC is used in Sect.  to validate the results against radiometer liquid water path (LWP) measurements. Finally, in Sect.  we perform an analysis of the adiabaticity of the retrieved LWC profiles, while Sect.  describes the benefits of using G-band frequencies for the LWC retrieval.

The LWC retrieval described in this article has benefited from the EPCAPE field campaign, which allowed different instruments to gather and operate simultaneously to enhance the scientific outcomes of the proposed research. The KAZR (Kollias et al., 2016) and the CloudCube G-band radar module (Socuellamos et al., 2024a) are the two main instruments used in this study. For simplicity, CloudCube's G-band module will be abbreviated as CloudCubeG in the following. KAZR and CloudCubeG operational parameters during the EPCAPE campaign are described in Table .

Besides both radars' data, which form the basis for the LWC retrieval, observational data from other instruments have also been used. In particular, we have worked with ceilometer data to identify the cloud base and have made use of microwave radiometer (MWR) measurements to validate our retrieval approach by comparing the integrated LWC, i.e., the liquid water path (LWP), with that derived by the radiometer. Additionally, radiosonde data have also been considered to input vertical profiles of temperature, pressure, and relative humidity into our model and to calculate the gaseous and liquid attenuation experienced by the radar signals. Figure  shows a picture of the location of all the instruments while deployed during the EPCAPE field campaign, also including the location from which the radiosondes were released. The distance between them was not larger than a few tens of meters.

We have selected two case studies on 2 different days for the analysis developed in this work. The first case, recorded on 17 April 2023 at 21:19:00 UTC starting time, lasts 100 min, while the second dataset, taken on 25 April 2023 at 13:12:00 UTC starting time, continues for 20 min. The radar data that have been used for the LWC retrieval consist of close to 18 000 and 1800 profiles of observed echo power at G-band and Ka-band, respectively, which have been averaged over 60 s to retrieve 120 LWC profiles. The data correspond to shallow stratocumulus formations with thickness between 120 and 360 m, containing characteristic LWC and LWP amounts smaller than 1 gm-3 and 200 gm-2, respectively. These cases, particularly challenging for the DFR technique, were chosen to demonstrate the advantages of the Ka-band and G-band proposed retrieval.

Numerous sources of uncertainty in the LWC retrieval have been identified related to the use of the dual-frequency radar technique (Williams and Vivekanandan, 2007), including different radar parameters and sampling mismatches between the radars, the presence of ice particles in the cloud volume, potential Mie scattering contamination, and errors in the differential absorption factors used as inputs in the retrieval model.

The KAZR and CloudCubeG operated with different spatiotemporal resolutions as well as different beamwidths, as shown in Table . While these mismatches can contribute to the inaccuracy of the LWC retrieval due to cloud inhomogeneity, we have taken various steps to minimize these errors. First, we have integrated CloudCubeG range resolution and interpolated KAZR time resolution to obtain data with a common spatiotemporal gate spacing, resulting in 30 m × 0.4 s range-time-resolution cells. Second, to mitigate possible errors originating from this artificial matching and the beamwidth mismatch, we have averaged the data along different range and time gates, as detailed in Sect. . Averaging over time and range has been proven to improve the errors due to illumination volume mismatch (Hogan et al., 2005; Williams and Vivekanandan, 2007).

Furthermore, to avoid errors due to the presence of ice particles, the analysis has been focused on full-liquid shallow clouds, identified after studying the radiosonde temperature vertical profile and the G-band radar Doppler spectra. To reduce Mie scattering contamination, we have excluded regions below the cloud base measured by the ceilometer.

Finally, the differential gaseous and liquid absorption parameters needed as inputs for the LWC retrieval have been calculated based on in situ measurements of pressure, temperature, and relative humidity from radiosondes released next to the radars. On 25 April 2023, we have used the data from a radiosonde released 15 min prior to the start of the observation period. For 17 April 2023, the observation did not coincide in time with a sonde release. However, we have compared the data recorded from two sondes that were released about 2 h before and after the selected observation period. The changes in temperature, pressure, and relative humidity between these two sondes were less than 0.5 K, 7 hPa, and 3 %, respectively, giving us confidence that these properties remained relatively stable during the observation period. Even if there had been a variation of 1 K, Hogan et al. (2005), as well as Williams and Vivekanandan (2007), claimed that the propagation error into the LWC retrieval would be smaller than 0.03 gm-3 using a Ka-band and W-band pair of radars and that for more widely spaced frequencies the error would decrease further.

The reflectivity data obtained from KAZR and CloudCubeG as well as the mean Doppler velocity data from KAZR are plotted in Fig. . The black line describes the cloud base as measured by the internal ceilometer cloud-base detection algorithm (ceil data product; Morris, 2016) with temporal and range resolutions of 16 s and 30 m, respectively. Pixels below the cloud base have been removed for the analysis. The reflectivity shown in Fig.  has been corrected for gaseous attenuation computed following Rosenkranz (1998) using the radiosonde data. It is evident when comparing KAZR and CloudCubeG reflectivity plots that the reflectivity is larger at Ka-band than at G-band and that those differences appear to increase with height above cloud base, which is indicative of the differential attenuation signal used to derive the LWC profile. There is also some evidence for light drizzle sedimentation as observed by the non-monotonic increase in Ka-band reflectivity with height, especially when the cloud is several hundred meters thick. Finally, the Doppler velocity observation shows overturning cells with timescales on the order of 5 to 10 min, with a tendency for the scale of the overturning circulations to shorten as the cloud thickness decreases.

In single-frequency radar systems, the LWC has traditionally been retrieved using an empirical relationship with the radar reflectivity. These relationships provide a straightforward approach to estimate the LWC and can also serve as a point of comparison for different techniques. Assuming a lognormal distribution of the drop population, this relationship is typically of the form LWC=aZb, where Z is the radar reflectivity in linear units and a and b are parameters that describe the drop size distribution (DSD) in terms of the total drop concentration, No, and the distribution width, σ (Atlas, 1954; Sauvegeout and Omar, 1987; Fox and Illingworth, 1997; Miles et al., 2000; Matrosov et al., 2004). Using a Ka-band radar, the coefficients a and b have been discussed to be a=2.4 (for maritime clouds) or a=4.6 (for continental clouds) and b=0.5 (Miles et al., 2000; Matrosov et al., 2004). The values derived in Atlas (1954), Sauvegeout and Omar (1987), and Fox and Illingworth (1997) provide similar results and have not been included in this analysis for simplicity.

Figure  shows the correspondence between KAZR and CloudCube reflectivity values with the retrieved LWC. A clear trend can be deduced, where the LWC increases for greater values of observed reflectivity. In addition, Fig. a includes the empirical curves for the two cases in Miles et al. (2000), which show a large systematic low bias of the reflectivity-based approach compared with the values retrieved from the DFR.

While the parameter b is usually assumed to be constant and independent of the DSD characteristics, a can adopt a relatively broad range of values depending on No and σ. Using b=0.5 (as in Atlas, 1954, and Miles et al., 2000), a=(π/6)ρNo0.5exp⁡(-4.5σ2), where ρ is the liquid water density (Matrosov et al., 2004). In Miles et al. (2000), the drop concentration was found to be 75 cm-3 for maritime clouds and 280 cm-3 for continental clouds, while the distribution width was 0.38 for both cases. A plausible explanation for the discrepancy between the retrieved LWC and the empirical curves may come from the fact that the cloud selected in our observation is formed by a high number of droplets of similar small sizes. This characteristic is typical of continental clouds, which can have drop concentrations close to 1000 cm-3 and narrower distribution widths (Squires, 1958; Sassen et al., 1999). Furthermore, recent analysis has shown that observed drop distributions narrow substantially as drop concentration increases (Lebsock and Witte, 2023), which was confirmed in a cloud chamber (Chandrakar et al., 2016). This natural correlation will tend to amplify the dependance of the parameter a on the droplet number concentration.

Figure  shows the possible range of values that the parameter a can cover for different drop concentrations and distribution widths. The parameter a approaches values of 16 for realistic concentrations of 1000 cm-3 with σ < 0.2.

Figure  shows the cloud droplet number concentration derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) aboard the AQUA satellite during two overpasses at approximately 20:30:00 UTC on 17 and 25 April 2023. The cloud optical depth and effective radius derived from the 3.7 µm channel are used in the derivation following the approach described in Grosvenor et al. (2018). As seen in Fig. , drop concentrations were close to or higher than 300 cm-3 near the observation site for both days. Based on this information, even if the overpasses did not coincide in time with the radar measurements, we assume that droplet concentrations could have exceeded 300 cm-3 during our observation windows. For a combination of No > 300 cm-3 and σ < 0.35, a can become greater than 9. Taking this value as a reference, we can provide a new relationship for the reflectivity-based LWC retrieval: 6LWC=9Z0.5. The curve obtained by Eq. () is also included in Fig. a. This relationship provides a better fit with the retrieved LWC, which suggests that the reflectivity-based approach can provide a complementary result to the DFR technique and potentially constrain the drop concentration in addition to the LWC.

An alternative method to account for the radar attenuation is based on the Hitschfeld and Bordan (1954) scheme, where the reflectivity would be corrected at each range gate prior to computing the LWC from a reflectivity-based relationship. For a range resolution of 30 m and LWC smaller than 1 gm-3, the attenuation by cloud liquid water at Ka-band is expected to be on the order of 0.01 dB per range gate as computed using ITU (2023). Even over the depth of a 300 m cloud, the cumulative attenuation at Ka-band would be on the order of 0.1 dB. Therefore, the reflectivity-based retrieval is not significantly affected by attenuation in this particular scene.

During the observation periods, a microwave radiometer (MWR) at the EPCAPE site recorded the liquid water path (LWP) among other environmental parameters (mwr3c data product; Cadeddu, 2024). The LWP retrieved by the radiometer has been used to validate the LWC retrieval technique by integrating the LWC across the cloud height.

Figure a shows the comparison between the radiometer and the radar-retrieved LWPs. In addition, we have included the adiabatic LWP, calculated assuming a linear increase in the LWC with range, as explained in Sect. . We have also added the LWP following the empirical relationships between the LWC and KAZR reflectivity, as discussed in Sect. . All LWP curves have been averaged over 60 s to match the time resolution of the radar LWP. The radar-reflectivity-based retrievals based on Miles et al. (2000) are again observed to be clear outliers here, with a systematic underestimate of the LWP relative to all other methods. On the other hand, the LWP calculated following Eq. () shows excellent correlation with the MWR retrieval. The estimated uncertainty in the MWR-derived LWP is 15 gm-2 (Cadeddu et al., 2020), while the estimated retrieval uncertainty of the DFR LWP is on the order of 10 gm-2. On 17 April 2023, the temporal trend of the radiometer LWP is well represented by the DFR retrieval, matching the progressive temporal decline of the LWP associated with a shallower cloud towards the end of the observation. Also on the same day, for both the DFR- and radiometer-retrieved LWP, the values are generally below the adiabatic curve for approximately the first 65 s of the observation window. Towards the end of the recording period, after 65 s, where the cloud becomes shallower, all three curves tend to match. This is an indication that the thermodynamic behavior of the cloud may be related to its thickness and is consistent with the observation that the thicker a cloud becomes, the greater its subadiabaticity (Albrecht et al., 1990; Rauber et al., 2007). A comprehensive analysis of the LWC and LWP adiabaticity is performed in Sect. .

Figure b shows a more detailed comparison of the KAZR–CloudCubeG DFR and the MWR-retrieved LWPs. Overall, the relative bias between the two retrievals shows large variability along both observation periods. The mean bias, represented as a dotted line in Fig. b, is -15.9 gm-2 on 17 April 2023 and -33.5 gm-2 on 25 April 2023, which indicates that the KAZR–CloudCubeG-retrieved LWP is, on average, greater than that retrieved by the MWR. This relatively low mean error suggests that a comparison with more temporal averaging will progressively improve the agreement between the two retrievals. The root mean square (rms) deviation, plotted with a dashed line in Fig. b, is 33.3 gm-2 on 17 April 2023 and 79.9 gm-2 on 25 April 2023, which are larger than what can be explained by the analytic uncertainty estimates in the radar (10 gm-2) and radiometer (15 gm-2) uncertainties. Given the disperse variability in the sign of the relative bias and the relatively large rms error, we speculate that there are random sources of uncertainty in the radar retrieval which have not been accounted for. Two likely causes could be radar–volume mismatches or an unidentified influence of drizzle drops on the DFR. However, following a similar validation approach, Hogan et al. (2005) compared the retrieved LWP from a Ka-band and W-band radar system with that from collocated MWRs. For similar amounts of LWP and the same time average as in this work, they found the rms deviation to be 40 and 70 gm-2 for two different cases and radar systems. While the Ka-band and W-band radars in Hogan et al. (2005) had different beamwidths, as with the KAZR and CloudCubeG here, the radars in that study were operated with the same range resolution. Furthermore, profiles with Mie scattering were identified, corresponding to differential Doppler velocities greater than 0.1 ms-1, and excluded from the analysis, which suggests minimal Mie scattering contamination. This speaks for the greater accuracy of the LWC retrieval using a Ka-band and G-band radar system. Even if radar–volume mismatches and Mie scattering contamination may be happening, the rms deviation shows similar agreement between the KAZR–CloudCubeG DFR and the MWR-retrieved LWPs compared to that reported in Hogan et al. (2005). If in future deployments better care was taken to match the volume sampling of the CloudCubeG and KAZR instruments, we suspect that a significant fraction of the unresolved scatter between the radiometer and the DFR retrieval would be eliminated.

To assess the thermodynamic behavior of the shallow cloud, we can compare the linearity of the LWC retrieval to the adiabatic case where the LWC increases linearly with range. In that scenario, the change in the adiabatic liquid water content, LWCad, is driven by the adiabatic lapse rate, Γad, as LWCad(r)=Γadr. Then the ratio between the retrieved LWC and the adiabatic LWCad describes the adiabatic behavior of the particular cloud layer analyzed, which is defined through the degree of adiabaticity, fad. There are three possibilities: (1) if fad=1, the LWC in the cloud parcel is exactly the same as if lifted adiabatically within the cloud range step, and therefore there is no exchange of energy with the environment; (2) if fad<1, the LWC is subadiabatic, indicating than the cloud layer is losing liquid water due to either mixing of warm dry air or sedimentation of drops; and (3) if fad>1, the LWC is superadiabatic, which can plausibly only be explained by sedimentation of small drizzle drops from above. The degree of adiabaticity is often summarized as a layer mean quantity. In such a case, the LWC is replaced by the LWP, and the degree of adiabaticity is equivalently defined as 7fad=LWPLWPad.

To compute Γad, we followed Albrecht et al. (1990), using the temperature, pressure, and humidity measurements of the released radiosondes. Γad was found to have an average value of 2.85 gm-3km⁻¹ within the cloud volumes on both days.

Figure  shows the LWC and the LWCad for the different cloud thicknesses measured along the observation period, which range from 120 to 360 m in steps of 30 m. For this representation, we took all the retrieved 60 s profiles that had the same cloud thickness and averaged them. We have also included the standard deviation of the averaged LWC profiles as shaded areas in Fig. . For all the cases except the one corresponding to 300 m on 17 April 2023, more than eight profiles were available. The first value at cloud base of the retrieved LWC curves was found, on average, to be biased high by approximately 0.1 gm-3, which could correspond to the LWP bias observed in Sect. . For a better visual comparison, we have shifted the retrieved LWC profiles in Fig.  so that they start exactly at 0 gm-3. This shift has not been applied to compute fad.

The degree of adiabaticity based on the LWPs is studied separately in Fig. . On 17 April 2023, fad is greater than 1 for clouds shallower than 200 m even within the uncertainty limits, manifesting a superadiabatic behavior of such clouds. While the uncertainty in the LWC retrieval is larger for shallower cloud sections, which may impact this outcome, this feature agrees with the results reported in Zhu et al. (2019) and has been speculated to be associated with lateral detrainment of cloud water from nearby convective elements (Miller et al., 1998). This analysis is supported by Fig. , which shows the mean column Doppler velocity and the predominant direction of the cloud profiles for the different cloud thicknesses on 17 April 2023. With the exception of the shallowest cloud section, which has the greatest uncertainty, the superadiabatic behavior is correlated with vertical updrafts, where the mixing of lifted air with humid air near the cloud top may be inducing the sedimentation of cloud droplets. On the other hand, clouds thicker than 200 m show subadiabatic behavior and negative Doppler velocity, suggesting that they sustain more acute mixing with a drier environment, presumably due to the evaporation of liquid water towards the cloud top or the loss of cloud water due to coalescence and sedimentation processes.

As described throughout the previous sections, the use of CloudCube's 239 GHz radar has allowed us, to the best of our knowledge, to retrieve the LWC with improved accuracy in comparison to any other results reported in the literature. This is a direct product of the smaller random error in the reflectivity measurement and the use of a complementary lower-frequency radar, with more than 200 GHz of separation between the two frequencies of operation, that provides increased differential attenuation. To quantify the benefits of using a higher-frequency radar at G-band, we have theoretically compared the LWC uncertainty to a pair of radars where the lower-frequency one operates at either 35 or 94 GHz. To perform this analysis, we have used Eq. (), keeping the same radar parameters for all the cases and only sweeping the frequency of operation. By taking a closer look at Eq. (), we can infer that only the random error in the differential frequency reflectivity ratio, ΔDFR, and the differential liquid attenuation, Al, are frequency-dependent parameters. These are precisely the two parameters that will eventually change the LWC retrieval accuracy based on the operating frequency of the radars. To understand their behavior, Fig.  shows their relative improvement (meaning a decrease in ΔDFR and an increase in Al) with respect to the Ka-band and W-band pair for the two scenarios previously mentioned.

These two parameters have an opposite effect depending on the lower radar frequency. The relative improvement of the DFR random error is greater if the lower radar frequency is 94 GHz. For the case where the higher-frequency radar operates at 239 GHz, as in this work, pairing with a 35 GHz radar offers an improvement of 8 %, while ΔDFR drops by 37 % if it is combined with a 94 GHz radar. This can be understood by bearing in mind that, for the same radar parameters, the lower frequency possesses a greater reflectivity random error that dominates the ΔDFR term (see Eq. ). On the other hand, the differential attenuation between 35 and 239 GHz offers a substantial increase in comparison with the 35 and 94 GHz or 94 and 239 GHz pairs. In particular, the improvement in Al is about 175 % or 2.75 times with the former and about 60 % or 1.6 times with the latter.

While the use of a higher-frequency radar at 239 GHz clearly offers an improvement in the LWC retrieval accuracy with regard to the 35 and 94 GHz radar system, whether combining it with a lower-frequency radar at 35 or 94 GHz depends on how ΔDFR and Al modify the LWC uncertainty. The results for the LWC after introducing the variation shown in Fig.  are presented in Fig. . For any of the lower-frequency radars, the LWC retrieval accuracy is always substantially improved if the higher-frequency radar is at 239 GHz compared to the 35 and 94 GHz radar pair. Interestingly, both curves tend to meet as the higher frequency is increased. At 239 GHz, the relative improvement for both cases is around 65 %, meaning that the LWC can be theoretically retrieved to values 3 times smaller in comparison to the Ka-band and W-band radars.

The proximity of both curves at 239 GHz also indicates that simultaneously retrieving the LWC from a triple-frequency radar system can be beneficial for the LWC retrieval. The Ka-band and W-band pair could be used at regions where the G-band signal suffers from Mie scattering, while combining Ka-band and G-band as well as W-band and G-band in Rayleigh scattering locations would provide two independent retrievals that could potentially improve the accuracy by a factor of 21/2. Furthermore, the triple-frequency system can be used to identify Mie scattering regions if Doppler is not available and to facilitate the removal of biases and isolate sporadic instrument artifacts.

Shallow warm clouds have extensive areal coverage over the tropical and subtropical oceans, and their strong radiative effects make them an important contributor to Earth's energy balance and hydrological cycle. The vertical profile of the liquid water content (LWC) in these warm clouds is key to better understanding their radiative properties, their ability to precipitate, and entrainment and mixing with free-tropospheric air.

Single- and dual-frequency radar techniques have been previously used to estimate the vertical profile of LWC in these shallow cloud systems. Single-frequency algorithms suffer from assumptions regarding the presence or absence of drizzle particles, while dual-frequency algorithms with radar frequencies up to W-band do not have the sensitivity needed to capture the fine structure of the LWC. The dual-frequency radar retrieval accuracy is intrinsically related to the frequencies of operation of the radar system, where higher frequencies contribute to reducing the measurement error and increasing the liquid water attenuation, then leading to more accurate retrievals. We have demonstrated in Sect.  that the synergy between a 35 and a 239 GHz radar offers distinct advantages to the LWC retrieval accuracy in comparison to a 35 and a 94 GHz radar, eventually allowing measurements of 3 times smaller liquid water amounts to LWC values as low as 0.05 gm-3 in the inner cloud structure and 0.22 gm-3 at the cloud base.

We have accomplished this by deploying the unique 239 GHz radar module of the CloudCube ground-based prototype instrument during the EPCAPE field campaign and taking advantage of simultaneous observations of the 35 GHz KAZR and a variety of additional instruments. In total, the analysis included 2 different days and more than 100 profiles of shallow clouds with typical amounts of LWC and liquid water path (LWP) smaller than 1 gm-3 and 200 gm-2, respectively, that were selected to demonstrate the benefits of the proposed technique in a particularly challenging scenario. After deriving the dual-frequency reflectivity ratio for a 60 s average time to reduce artifacts and errors, we retrieved the LWC across the cloud thickness after averaging a series of sliding and overlapping LWC short profiles. The dual-frequency radar reflectivity (DFR) technique was able to capture the vertical variability of the LWC, as demonstrated by the comparison of the LWC profiles for different cloud thicknesses with their adiabatic behavior. We found that clouds shallower than 200 m tend to be in a superadiabatic state, while clouds deeper than 200 m behave subadiabatically. On the other hand, we showed how a complementary reflectivity-based retrieval can be useful to constrain the cloud droplet concentration number in addition to the LWC. We compared the retrieved LWP to that measured by the collocated microwave radiometer (MWR). While the direct comparison of 60 s data points showed large variability, the consideration of the mean LWPs showed improved agreement, suggesting that longer average times would lead to better agreement although negatively affecting the temporal variability of the retrieval. The rms difference between the MWR- and DFR-derived LWP was 33 and 79 gm-2 for the 2 observation days, which is a few times larger than can be explained by the analytic uncertainty estimates in either retrieval. We speculate that small differences in the sample volume between the G-band and Ka-band radars are the root cause of this residual uncertainty and recommend more attention be paid to matching the radar sample volumes in future studies.

The development of G-band radars that are sensitive to smaller characteristic drop sizes is offering new possibilities within the atmospheric remote sensing field. We have proven that G-band radars are a valuable instrument for the LWC retrieval that, when complemented with a Ka-band or a W-band radar, can potentially determine very small amounts of liquid water content in low-level shallow clouds. These results can be useful to better understand the vertical structure of shallow clouds and how they interact with the surrounding atmosphere.

Future work will be devoted to testing, improving, and validating the results retrieved in this work. KAZR and CloudCubeG will be participating in the Cloud and Precipitation Experiment at Kennaook (CAPE-k), where close to 1 year of data will be available by the end of 2025. For CAPE-k, KAZR and CloudCubeG will be accompanied by the W-band ARM Cloud Radar (WACR), and we will exploit the availability of a triple-frequency radar system to benefit the LWC retrieval as described in Sect. . Finally, we intend to compare the results against a radiometer and single-frequency radar joint retrieval, as combining the radiometric observation of total optical depth with the radar capability of profiling can also be used to constrain the cloud LWC.