Joint project between
University of Bonn and University of Mainz, Phase 1
University of Bonn: Armin Blanke (PhD student) and Silke Trömel (PI)
University of Mainz: Manuel Moser (PhD student) and Christiane Voigt (PI)
POLICE will exploit radar polarimetry for quantitative process studies in mixed phase clouds and for model evaluation. In-situ measurements of microphysical cloud parameters will be performed with the research aircraft HALO and the Falcon in the dendritic ice growth layer (DGL) and the melting layer below to evaluate targeted hypotheses on the origin of enhanced specific differential phase KDP in the DGL and quantify different indicators for the discrimination between aggregation and riming processes in the clouds. This will enable the set-up of a discrimination algorithm based on polarimetric weather radar measurements. Both the explanation of KDP-bands and the ability to distinguish between aggregation and riming is crucial for data assimilation and model microphysics above the melting layer. The in-situ measurements are used to evaluate the state-of the art retrievals of particle number concentration Nt, mean particle diameter Dm and ice water content IWC from the polarimetric radar signals and to evaluate the representation of particle type and size in the ICON-LAM model. The use of advanced spectral bin microphysical schemes in the Hebrew University Cloud Model (HUCM) in combination with measured cloud profiles will
also give insights in processes responsible for potential deficiencies in ICON-LAM in order to enhance the
representation of the cloud’s microphysical and radar derived properties in ICON-LAM.
Figure 1: Schematic of the methods and goal of the POLICE project: in-situ measurements of cloud microphysical properties will be evaluated along with polarimetric radar retrievals and cloud and weather models to investigate processes in mixed phase clouds.
The distinction between aggregation and riming below the dendritic growth layer (DGL) is important, because the latter signals the presence of super-cooled liquid water (SLW). SLW maybe partially responsible for additional ice generation via the Bergeron-Findeisen process. Riming also favors secondary ice production through the Hallet-Mossop ice multiplication process (rime splintering), which is active between -3°C and -8°C. Identification and distinction of regions with dominant riming and aggregation has implications for data assimilation and model microphysics above the melting layer. Following the demand for an area-wide algorithm to distinguish between aggregation and riming exploiting national weather radar networks, POLICE exploited and analyzed quasi-vertical profile (QVP) data of reflectivity (ZH), differential reflectivity (ZDR) and in particular depolarization ratio (DR). Similar to ZDR, DR tends to decrease in rimed snow relative to aggregated snow, but the corresponding difference in DR is 2-4 dB larger (e.g., Ryzhkov et al., 2017). DR can be estimated by dual-polarization radars operating in SHV mode (simultaneous transmission/reception of orthogonally polarized waves) and serve as a proxy for circular depolarization ratio (Matrosov, 2004). Naturally, DR combines the information content of ZDR and cross-correlation coefficient (ρhv) in a single quantity. To set up a riming detection algorithm, the Doppler spectra, available from the vertically pointing C-band radar measurements of the German weather radar network (to date the only operational network providing such data) and post-processed with the procedure described in Gergely et al. (2022) served as input. It is the newly developed mean isolated spectra profile (MISP) technique, which provides clutter free and noise reduced spectra data in a height vs. time format. The MISPs of mean Doppler velocity (MDV) are used to identify regions with particles falling faster than 1.5 m/s and accordingly associated with riming. In addition, the degree of riming of the column above the radar is retrieved by making use of the MDV MISPs by determining the rime mass fraction (RMF) according to Kneifel and Moisseev (2020). Figure 1 shows an investigated test case, which further exhibits a correlation of 0.7 between the Doppler velocities and DR, emphasizing the strong potential of this variable for riming detection.
Figure 1: Scatter plot and linear regression of ZDR vs. DR observed with the DWD C-band radar in Essen on 02.01.2022 between 0300 to 0930 UTC. The coloring of the individual data point indicates MISPs mean Doppler velocities. The correlation coefficient r is provided for ZDR vs. DR.
In order to develop the algorithm, relationships between different polarimetric variables and Doppler velocities can be learned supervised. Four methodologies were tested and trained using QVP data of ZH, ZDR, and DR as input from four days containing five precipitation events observed with the DWD Essen radar (ESS) so far. The four methodologies investigated are the gradient boosting model (GBM) based on decision trees, the quadratic discriminant analysis (QDA), the logistic regression (LR) and the multilayer perceptron (MLP) artificial neural network trained with a back propagation algorithm. The ground truth and the results of the riming detection algorithm have been transformed in binary fields: riming (1) and no riming (0).
The performances of the resulting algorithms have been evaluated by computing different widely used scores/performance measures such as the confusion matrix, balanced accuracy, Kappa value, F1-score and Matthews correlation coefficient (MCC; Matthews, 1975).
The algorithms have been evaluated using 70% of the data for model training and 30% of data for the performance assessment. Prior to training, the data set was prepared to obtain the same ratio of riming to non-riming cases as in the original data set.
Due to spatial and temporal mismatches between QVPs and MISPs, due to retrieval uncertainties, and the constraint of pixel-by-pixel comparison arising from the binary analyses, the experimental design does not allow an optimal score of the metrics of 1. However, in summary, the GBM-based riming algorithm performed best on the test data set with a balanced accuracy of 71%, an F1 score of 0.61 and a normalized MCC of 0.73.
Figure 2: Binary time-height plots of mean Doppler velocities faster than 1.5 m/s (top panels, a-d) and corresponding GBM retrieval results (bottom panels, e-h) from four selected example cases on 13 May 2021 between 1515 to 2030 UTC (a, e), 03 November 2021 between 1130 to 1500 UTC (b, f), 02 November 2022 between 0300 to 0900 UTC (c, g) and 02 November 2022 between 1700 to 2100 UTC (d, f). The yellow color displays (predicted) riming, while the purple color indicates no (predicted) riming.
In an ensuing step, the retrieval prediction is additionally smoothed in time and height to account for possible mismatches. Overall, the final GBM-based retrieval shows promising results, with a balanced accuracy of 79 %, an F1 score of 0.64 and a normalized MCC of 0.78 when tested with the complete initial input data set. Figure 2 gives an impression of the performance of the new GBM retrieval. In all cases investigated the riming pattern is nicely represented and the retrieval is able to detect riming with a mean RMF of 0.47 (1.5 m/s corresponds to a RMF of approximately 0.46 at C-band).
The extreme stratiform precipitation event observed on 14 July 2021 in western Germany with devastating floods in the Ahr valley in Rhineland-Palatinate was selected to validate the performance of the final GBM retrieval independently. Despite obvious time shifts leading to double penalties affecting the scores, the prediction shows convincing results (see Fig. 3) with a balanced accuracy of 67 %, an F1 score of 0.47 and a normalized MCC of 0.69. These metrics are in the same order of magnitude as those obtained for the test data set.
Figure 3: Binary time-height plots of mean Doppler velocities faster than 1.5 m/s (left) and corresponding GBM retrieval results (right) from the flooding case on 14 July 2021 between 0100 to 1830 UTC observed by ESS. The yellow color displays (predicted) riming, while the purple color indicates no (predicted) riming.
The validation of the retrieval with more independent cases is still ongoing and the extension of the data set will help to make the algorithm more robust. However, the key development for a radar-based riming detection algorithm has been made.
References:
Within the scope of the POLICE project, the University of Mainz has supported aircraft measurements of clouds in the Arctic and their thermodynamic phase. During the three aircraft field campaigns “Aircraft campaign observing FLUXes of energy and momentum in the cloudy boundary layer over polar sea ice and ocean” (AFLUX), “Atmospheric airborne observations in the Central Arctic” (MOSAiC-ACA; Shupe et al., 2022) and HALO-(AC)3 conducted within the framework of the “Arctic Amplification: Climate Relevant Atmospheric and Surface Processes, and Feedback Mechanisms (AC)3” project (Wendisch et al., 2017), a comprehensive data set of microphysical cloud properties was collected. The research aircraft Polar 5 and 6, operated by the Alfred Wegener Institute (AWI), was used as a platform to conduct in total 33 research flight in the northern Fram Strait between Greenland and Svalbard for remote sensing and in-situ measurements of Arctic clouds.
Figure 1: Maps of the flights during AFLUX, MOSAiC-ACA and HALO-(AC)3 in the vicinity of Svalbard, Longyearbyen (LYR). Back- ground shows the sea ice concentration at the halftime of each campaign recorded by the Advanced Microwave Scanning Radiometer 2 (AMSR2) onboard the GCOM-W1 satellite.
The in-situ data obtained during AFLUX and MOSAiC-ACA (Mech et al., 2022) have been used to study the microphysical properties and thermodynamic phases of Arctic low-level clouds during the season of maximum and minimum sea ice extend. We investigated the distributions of particle number concentration N, effective diameter Deff and cloud water content CWC (liquid and ice) and developed a method to quantitatively derive the occurrence probability of their thermodynamic phase from the combination of microphysical cloud probe and Polar Nephelometer data. Finally, we assessed changes in cloud microphysics and cloud phase related to the ambient meteorological conditions in spring and summer and addressed effects of the sea ice and open ocean surface conditions. A southerly flow from the sea ice in cold air outbreaks dominates cloud formation processes at temperatures mostly below -10°C in spring, while northerly warm air intrusions favor the formation of liquid clouds at warmer temperatures in summer. We find large differences in the particle sizes in spring and summer and an impact of the surface conditions, which modify the heat and moisture fluxes in the boundary layer. The in-situ data show that mixed-phase clouds are the dominant thermodynamic cloud phase in spring with a frequency of occurrence of 61% over the sea ice and 66% over the ocean. Pure ice clouds exist almost exclusively over the open ocean in spring, and in summer the cloud particles are most likely in the liquid water state (Moser et al., submitted to ACPD, 2023).
Figure 2: Frequency of occurrence for each particle regime (1a, 1b: Ice particles; 2a, 2b, 2c: Mixed-phase particles; 3: liquid particles; 4: Aerosol particles), separated by season and surface conditions. The values are normalized by the respective environmental conditions (spring-ice, spring-ocean, summer-ice and summer-ocean).
Further analyses with cloud data recently collected during HALO-(AC)3 in March and April 2022 will be used to extend the in-situ measured Arctic microphysical cloud statistics. During HALO-(AC)3 three aircraft were flying in formation collecting closely spatially collocated and almost coincident in-situ and remote sensing including non polarimetric radar observations. The data can help to evaluate remote sensing and satellite retrievals (e.g. Klingebiel et al., 2023), and may elucidate the role of clouds in the region of the world with strongest anthropogenic climate change (Wendisch et al., 2022).
References:
Joint contribution of University of Bonn and University of Mainz
A large dataset of in-situ microphysical measurements from several flight tracks obtained during the Olympic Mountain Experiment (OLYMPEX; Houze et al., 2017) are analysed to assess the accuracy of state-of-the-art polarimetric microphysical retrievals.
During OLYMPEX, the National Science Foundation (NSF) Doppler On Wheels (DOW; Houze et al. 2018) polarimetric X-band radar performed repeatedly sequences of Range-Height Indicator (RHI) scans within an azimuthal sector of 22 degrees. Similar to Quasi-Vertical Profiles (QVPs; Ryzhkov et al. 2016), Columnar-Vertical Profiles (CVPs; Murphy et al. 2020) and RHI scan-based Vertical Profiles (R-QVP; Allabakash et al. 2019), our technique presents noise-reduced polarimetric variables in a height versus time format. The RHIs performed within an azimuthal sector over a certain time period are averaged within a predefined window size to extract vertical profiles of polarimetric variables at the location of interest. The vertical profiles obtained from ensuing sector-averaged RHIs are combined in a time versus height format and used as data base for polarimetric ice microphysical retrievals. Because this technique considers a more limited spatial domain for averaging, it is well suited to also study flight transects in horizontally less homogeneous cloud fields. During OLYMPEX, the UND Citation research aircraft (Poellot et al. 2017), equipped with advanced measurement devices, flew over the DOW 148 times. The in-situ collected microphysical cloud properties including cloud particle number concentration Nt , particle size distribution, mean volume diameter Dm, and two-dimensional images of hydrometeors are reprocessed and analyzed to enable the direct comparison with the microphysical radar retrievals along 3D trajectories of the research aircraft. For the joint analyses University of Bonn focused on the processing of radar data and the microphysical retrievals, while the University of Mainz processed the in-situ measurements. The combined application of IWC(Zh, Kdp) based on horizontal reflectivity Zh and specific differential phase Kdp (Bukovčić et al. 2018; 2020) in regions where linear differential reflectivity Zdr<0.4 dB and IWC (Zdr, Kdp) (Ryzhkov et al. 2019) elsewhere, as suggested by Carlin et al. (2021), shows the best agreement with the in-situ measurements (Pearson correlation coefficient r= 0.96, RMSE=0.19 g/m³, Fig. 1). The hybrid retrieval by Carlin et al. (2021) outperforms all other aforementioned retrievals and also the non-polarimetric retrieval by Hogan et al. (2006). However, we see a tendency for a slight overestimation of IWC at warmer and an underestimation at colder temperatures.
Figure 1: Comparison of polarimetric IWC (g m-3) retrievals (dots) following Carlin et al. (2021) with collocated airborne in-situ measurements (solid line). Colours indicate the temperature in °C. Vertical bars represent in-situ standard deviations.
Polarimetric retrievals for the particle number concentration Nt following Ryzhkov and Zrnić (2019) show more pronounced deviations from in-situ measurements, but still reach a convincing correlation of r=0.88. Overall, the comparison of most recent polarimetric retrievals of mean volume diameter Dm (Ryzhkov et al. 2019), IWC, and Nt with collocated in-situ observations during the OLYMPEX campaign shows agood agreement when certain requirements are fulfilled, e.g. a flight length of at least 30 seconds through the radar sector used for the comparison.
The valuable data set of collocated radar and in-situ measurements extracted from OLYMPEX is further enlarged by aircraft field campaigns including BLUESKY (Voigt et al. 2021) and CIRRUS-HL where in total five different stratiform cloud situations over the DWD polarimetric radar Hohenpeißenberg were sampled with the DLR research airplanes Falcon and HALO. Figure 2 shows a stratiform precipitation event monitored with Hohenpeißenberg radar on 23 May 2020 (right) and collocated (in space and time) with in-situ cloud measurements conducted by the DLR research aircraft Falcon equipped with a Cloud Imaging Probe (left). Measuring microphysical cloud properties at 4 different altitudes, each flight leg is dominated by a different hydrometeor type: mainly aggregates in the highest leg at 4 km above the radar, rimed particles at 3 km, partly melted particles in the melting layer and droplets below.
Figure 2: Quasi-vertical profiles (QVPs) of logarithmic differential reflectivity ZDR (right panel) monitored with the polarimetric radar at Hohenpeißenberg during the BLUESKY field campaign together with collocated in-situ particle measurements (left panel) recorded at the same time by the Cloud Imaging Probe (CIP).
References:
The work in the first year addresses the evaluation of existing in-situ measurements of clouds for radar intercomparison (WP 1). In-situ cloud measurements performed during an overflight of the polarimetric radar station in Bonn (BoXPol) during the ECLIF2/NDMAX flight campaign in February 2018 have been evaluated. The horizontal reflectivities measured by BoXPol along the flight path of the research aircraft NASA DC-8 and the particle number concentrations in the particle size range from 15 to 960 µm measured with the cloud imaging probe CIP onboard the DC-8 are shown in Figure 1 The radar’s blind spot is marked in grey and the red line indicates the time when the radar scan was conducted. In the area close to the radar station a thick cloud layer with particles up to 700 µm in diameter was observed giving rise to the enhanced radar reflectivity. Further away from the radar station the observed cloud particles were smaller and below the detection limit of the radar.
The High Volume Precipitation Spectrometer HVPS has been ordered which will provide a better resolution of microphysical cloud properties detected near polarimetric radar stations in Germany (WP-2). A flight strategy is developed to augment the quality of the in-situ and radar intercomparisons in clouds.
Figure 1: Horizontal reflectivities ZH (top) along the flightpath of the NASA DC-8 research aircraft and the simultaneously measured in-situ particle number concentrations N (bottom).
Contribution of University of Mainz
The contribution from the University of Mainz is focusing on airborne in-situ cloud particle measurements over polarimetric radar stations using so-called Optical Array Probes (OAPs). Hereby, microphysical cloud properties are extracted from two dimensional images recorded with OAPs inside clouds (see Fig. 3). During flight, cloud particles pass through a collimated laser beam and their shadows are projected onto a linear array of 128 photodetectors. The presence of a particle is registered by a change in light intensity on each diode. The registered changes in the photodetectors are stored at a rate consistent with the aircraft’s velocity and the instrument’s size resolution. Particle images are reconstructed from individual “slices”, where a slice is the state of the 128-element linear array at a given moment in time. A slice must be stored at each time interval that the particle advances through the beam a distance equal to the resolution of the probe (see Fig. 4 for the working principle of an OAP and recorded images of hydrometeors).
Major work contributing to POLICE will be done with the High Volume Precipitation Spectrometer (HVPS) (see Fig. 3), which takes images of cloud particles in the size range from 150 µm to 1.92 cm. From the recorded 2D hydrometeor images, the cloud droplet number concentration, particle size distribution, median volume diameter can be derived. In addition, particles are categorized by hydrometeor type to distinguish between aggregation and riming of large ice crystals. Therefore, this measurement technique is promising to explain the origin of the enhanced specific differential phase KDP and to quantify indicators for the discrimination between aggregation and riming processes in polarimetric radar signals.
Figure 2: An Optical Array Probe (here: High Volume Precipitation Spectrometer) mounted below the wing of a research aircraft (1).
Figure 3: Schematic of an Optical Array Probe (top) and recorded images of hydrometeors.
References