the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Coincident in situ and triple-frequency radar airborne observations in the Arctic
Cuong M. Nguyen
Mengistu Wolde
Alessandro Battaglia
Leonid Nichman
Natalia Bliankinshtein
Samuel Haimov
Kenny Bala
Dirk Schuettemeyer
Download
- Final revised paper (published on 10 Feb 2022)
- Preprint (discussion started on 22 Jul 2021)
Interactive discussion
Status: closed
-
AC1: 'Updated preprint', Cuong Nguyen, 27 Jul 2021
Since the original submission, we’ve revised the manuscript significantly based on comments from the AMT editors and to reflect the fact that the second part of the paper is submitted to AMT (amt-2021-227). The manuscript preprint is updated to the latest version. Major changes include the revision of the abstract and summary sections. Section 5 is now merged into the Summary and Discussion section. Figures 12, 13, 15, 16, 18, and 19 in section 4 are updated. Our apologies for any inconvenience.
Citation: https://doi.org/10.5194/amt-2021-148-AC1 -
RC1: 'Comment on amt-2021-148', Anonymous Referee #1, 03 Aug 2021
The paper presents the measurement results from the airborne RadSnowExp campaign which offers near-simultaneous and coincident triple-frequency radar observation and in-situ (in-cloud) ice particle characterization. The observational setup is rather unique as it basically shifts what is commonly done on the ground (e.g. BAECC campaign in Finland, von Lerber et al 2017) on an airplane allowing for direct in-cloud particle imaging, but posing new challenges regarding the observational constraint given by the airborne platform.
The analysis of the dataset focuses on the connection between triple-frequency radar signatures at the X-Ka-W band and particle properties which is a very relevant subject for snow microphysical studies. The study confirms the existence of such connection as it is predicted by various modeling studies which are presented in the paper introduction. The results of the study are supporting the idea of using multifrequency radars for microphysical retrievals.Given the significance of the dataset presented I think the paper constitutes a valuable contribution to AMT. However, I have a few major comments that I suggest to be addressed before the paper is published.
1) Figure quality. I do not think that the presentation quality is sufficient for a final publication. Many figures are very hard to evaluate due to the fact that they are quite compressed in terms of the range of values. Also, text and labels are often hardly visible. Some significant work must be put on the figure quality.
The following are just suggestions connected to the aforementioned "readability" point, but it is up to the authors to take it or not. The number and size of the figures are significant, perhaps some work can be done also in this direction to rationalize the figure-load and facilitate the reading. As an example, Figures 12, 15, and 18 occupy an area comparable to the one occupied by Figure 13 and not equivalently discussed. Perhaps they could be probably be accommodated as subpanels of their respective "event-dashboards" (figures 13, 16, and 19) making this a complete overview of the measurements.
Figure 6 can be combined with Figure 5 giving a single general overview of the flight path and atmospheric conditions which connects well to the description given in the text.
Finally, I see a little relevance of figures 1 and 6 which are not necessary for the paper and can be moved to supplementary material or even left out. The study focuses on the measurements taken during the flight of 22 November and should only present data from that flight in my opinion. I basically had problems during the reading in following the various hierarchical groups: campaign->flights->segments->sections(A, B, C ...).2) Data availability. I did not find the data availability section. There are occasions where the paper specifically states the importance of the presented dataset which makes the data availability not only highly recommended, but quite essential to deliver the value of the study to the scientific community.
3) Scope and Uncertainties. The abstract (and in part also the summary) states that there is a "close relationship" between triple-frequency and particles' bulk density, level of riming, aggregation, and characteristic size of the PSD.
The degrees of aggregation and riming are not evaluated if not only qualitatively, but I do not understand from this paper how to use DFR to make a quantitative estimation of aggregation and riming degree.
Regarding bulk density, I do not see such close relation. Judging from figures 14, 17, and 20 it seems that bulk density is connected to mean size but can take various values at the same DFR range. Looking at Fig 21b it seems that high-density values are found for small DFR and on both left and right sides of the histogram. The "rotation" feature in the triple-frequency plot is not really evident. The range of bulk density values is very limited and skewed towards low densities which suggests a general problem in estimating this quantity. This also suggests that higher density values are found at the borders of the histograms due to problems of statistical representativeness (rare values are found in small samples). Density values in Fig 21 seem to correlate mostly with MVD rather than DFR.
Regarding MVD I think that a correlation with DFRs is clear. However, Fig 21a only shows the mean MVD for a combination of DFRs and does not show other important quantities such as the variance of MVD which I believe is essential for the retrieval study of Mroz (2021).Detailed Points
Line 115. Figure 2 - it is very difficult to connect the curves to the legend symbols. Perhaps enlarge the legend fonts or group the legend labels in different blocks according to their respective main group (already color-coded). Also, it is not totally clear to me how this is used in the study. If it is only for illustrative purposes or it is actually an attempt to connect with microphysical properties?
For example, It would be nice to connect the triple-frequency characteristic of these modeled particles with microphysical quantities as they are defined in section 3.2 lines 230-240. What are the MVD and bulk density of these modeled particles? How do they compare with the mean values measured for the same DFRs?Line 161-168 If the Ka and W band radars are absolutely calibrated, their return for Rayleigh ice particles should be around 1.2 dB and not 0. This is because of the frequency-dependent difference between the dielectric factor K for ice and water. The radars cannot be simultaneously calibrated in an absolute sense and have DFR=0 for small ice particles. Please clarify the calibration procedure.
This point is also discussed in Dias Neto et al. 2019 and Ori et al 2020Line 185 Figure 4 It is very difficult to evaluate a bias of 0.8 dB on a small scale that spans over 100 dBZ. Considering the objective of the figure I would cut it between -15 and 10 dBZ focusing on the upper part only
Line 236 It would be nice to include a formula also here like it is done for the other quantities. Usually it is define as
int_0^MVD V(D) N(D) dD = int_MVD^inf V(D)N(D)dD
where V(D) is the volume as a function of size.
Since the video disdrometer cannot measure the volume of snowflakes (which is also ill-defined considering that snowflakes' shapes are irregular) it is better to say also how volume is calculated here. Is it still assumed to be a spheroid with a 0.6 aspect ratio?
The given citation seems inappropriate to me. Leroy (2016) describes a methodology to calculate Median Mass Diameter (MMD) and it is not clear how this is connected to MVD.
Finally, the statement "This is the characteristic diameter that contributes most to cloud liquid water or mass" is confusing and incorrect to me. By definition, the size contributing the most to the mass should be the one that maximizes the function m(D)N(D)dD (i.e. the mode of the mass distribution). Even considering the volume equivalent to mass (by assuming constant density) stating that MVD is the size that contributes the most to the total volume would be again incorrect. The mode and the median value of distribution are in general diverse, this is especially true for multimodal distributions such it is the case in the presented case studies.Lines 291-294. I am not really sure if I can understand these sentences. First, a 10-minute running window corresponds to roughly 6 km considering the average ground speed. Is homogeneity important for this thresholding technique? How is the threshold of 0.6 identified? What do the authors mean by "accurate analysis"?
I guess that a good correlation is one easy indicator that the authors can use in order to connect measurements on-board of the aircraft and apart from it, but I wonder if this analysis could be biased by the characteristics of the measurements. As an example: If the cloud field analyzed is very homogeneous both measurements would result in a signal mostly dominated by random noise and thus even if the two signals are connected in reality the correlation coefficient would be close to 0.Line 320 Figure 10 This figure is not readable. I suggest the authors make much better use of the page real estate; increasing the vertical size of the figure, allowing for a better evaluation of the various curves, and significantly enlarging the font sizes.
Line 338 Figure 11 Enlarge axes font size of the legend.
Line 346 Is it possible that the 30um peak is due to the shattering of ice particles at the probes? Shattering is not discussed in the text. The reference list includes Lawson (2011) but that reference is not present in the text (Line 600).
Line 380 Figure 13 (the same applies to figure 16 and 19). I like these overview plots, but the Figures are barely readable at maximum magnification on a screen. Also, the subpanels are not labeled and it is difficult to follow the discussion on them. I suggest significantly increase the size of the figures. An idea to make better use of the page surface could be to put all time-plots on the left column sharing the same x-time axis and the ABCDE-sections classification. The left column could take up to 2/3 or even 3/4 of the figure width. Then, the snow images could be arranged on the right column. Also, I suggest reducing the number of ice images including only a few significant ones.
Line 480 Figure 22. What are the black lines?
Line 493 It is not clear to me where to find the relationships between ice particle properties and triple-frequency signature in this study. The paper presents a qualitative assessment of relations among these quantities
Lines 503 and 505. I guess here it refers to Figure 21 and not 22
Line 505 I actually see a very little sensitivity of estimated bulk density to triple frequency. From Fig 21b I do not see a transition from more reddish colors to blue/grey while "rotating" counterclockwise in the triple-frequency plot. Can the authors illustrate better this point?
Line 506 I saw that the "rotation feature" was much better illustrated in the first version of the manuscript uploaded. And I think that the text got it the other way around, or? A decrease in effective bulk density is expected when DFR X-Ka increases (counterclockwise rotation). For high values of DFR KaW and the low value of DFR X-Ka, we expect denser particles.
Minor PointsLine 17 Please introduce the CPI acronym
Line 18 DFR acronym is introduced later at line 21
Line 22 Double period ..
Line 22 I guess the phrase was intended without the word "that", but I would suggest rephrasing it anyway to make it easier to understand.
Line 56 Mismatched parenthesis ))
Line 95 I think there is a sign problem in the attenuation component of Eq 1. Assuming the attenuation to be semi defined positive such that the measured reflectivity Z=Ze-A then DFR = Z1-Z2 - (A1-A2) = Z1-Z2 + (A2-A1) [Lehrmitte 1990, Tridon 2020]
Line 232 Missing year in Heymsfield et al.
Line 288 misspelled Gans?
Line 350 Figure 12. The caption refers to panels (a) and (b) but the figure panels are not labeled. The same applies to Figures 15 and 18
Line 351 Text refers to left/right panels, but it is better to use panels labels (a) (b) according to AMT guidelines
Line 436 Figure 15. I guess the caption refers to sections selected from figure 16
Line 497 The MEAN particle diameter.
Line 546-550 I think that usually, 2020a comes before 2020b in the reference list.
Citation: https://doi.org/10.5194/amt-2021-148-RC1 - AC2: 'Reply on RC1', Cuong Nguyen, 02 Oct 2021
-
RC2: 'Comment on amt-2021-148', Anonymous Referee #2, 18 Aug 2021
The comment was uploaded in the form of a supplement: https://amt.copernicus.org/preprints/amt-2021-148/amt-2021-148-RC2-supplement.pdf
- AC3: 'Reply on RC2', Cuong Nguyen, 02 Oct 2021
-
RC3: 'Comment on amt-2021-148', Anonymous Referee #3, 26 Aug 2021
Review of the article titled “Coincident In-situ and Triple-Frequency Radar Airborne Observations in the Arctic” by Cuong M. Nguyen et al.
This article shows promising results of triple-frequency radar observations from the Radar Snow Experiment (RadSnowExp). Part of the uniqueness of this article is that both in situ and remotely sensed observations were collected by the same aircraft. Some complications arise when combining these datasets for their analysis, but the authors carefully and thoroughly describe the methodology used for volume matching and range calibration.
The authors studied the relationship between in situ sampled cloud microphysical properties and radar triple-frequency signals particularly, they centered their analysis on the relationships between median volume diameters, effective bulk density and dual frequency ratios. Lastly, the authors suggest a path forward with the possibility for quantitative retrieval of particle size using measured DFR but more in-depth analysis is needed to reach that step.
This article is generally well written, it is interesting, and it will be of interest to the scientific community. Nevertheless, I have a few concerns I suggest be addressed before this article is published. Below are some general and specific comments.
General Comments:
The authors analyzed one flight (22 November) and divide it into 3 different segments (1948-2000 UTC, 2005-2028 UTC and 2121-2135 UTC) giving a total number of 49 minutes of DFR observations. Please consider this to put into context the overall findings that derive from this dataset that are stated in this article.
Figures need significant improvement, both in the actual quality of the figure (suggest improving ppi) and the legibility of axis labels and legends. For figure with more than one panel the different panels should be labeled. This was done for some of the figures but not in all. I suggest a uniform way to address multipanel figures and their caption. Finally, figure captions need to carefully state what is plotted in each panel/figure.
I’m not sure figures 1 and 6 are necessary if the article centers around the Nov 22 flight (in case fig 1 is kept, I made I suggestion below to have all the flights drawn or the domain where all the flights took place). Similarly, I don’t think figure 5 is necessary either, no extra analysis is done of how the overlapping sizes were treated or any extra analysis that would make this plot needed in the article, just stating in the text the size ranges each instrument measures should be sufficient.
Specific Comments:
Lines 14-16: Consider stating the amount of data that was used to reach the conclusions stated in the article and not the overall flight hours of RadSnowExp.
Line 18: I’m not entirely sure this article showed how to accurately derive the level of riming from the DFR plane, consider rewriting this sentence.
Line 26: Add the definition of GPM-DPR
Lines 29-31: Reference needed.
Line 34: Suggest modifying this sentence as “The GPM Core Observatory carries…” or similar.
Lines 39-41: Reference needed.
Lines 63-65: Figure 1 does not support this statement, a figure showing all the flight tracks, or the experiment domain would be more useful.
Line 71: Consider replacing ‘uniquely’ by ‘unique’
Line 71: Figure 3 should be Figure 2, please reorder figures.
Line 110: Should be figure 3 when re-ordered.
Line 110: consider replacing ‘points‘ by ‘retrievals’ or similar.
Lines 149-154: The second bullet point is not clear. What do you mean by ‘three radar data’? How was it mapped into a common range axis? How ‘reasonable’ is the homogeneity assumption? Was any sensitivity analysis done to evaluate it? A schematic of all the smoothing methods could be valuable to assess the actual volumes that are being compared.
Lines 175-179: I consider this to be an important factor in the data analyzed here that could grant the inclusion of a figure showing these values to support this statement.
Line 180: Figures should be numbered in sequential order they are referenced in the text. Consider reorganize the figures or the text.
Line 181: Were there more than 1 flight on 22 Nov? How much data does this represent? i.e.: How many data points were used to reach this conclusion?
Lines 193-194: Consider adding sizes to contextualize ‘small cloud droplets’ and ‘large precipitation hydrometeors’.
Line 197: Consider replacing ‘or’ by ‘and’.
Line 206: Is the difference between all the 9 groups and the subset of ice habits only the inclusion or not of the Drops and Artifacts category? How are small particles treated?
Lines 212-213: It would be good to state the uncertainty values that are within the range presented by Baumgardner et al. (2017).
Lines 219-224: This paragraph is confusing, if the wiring had little effect on the estimated water content, then what was the factor that made the accuracy drop from 0.002 g/m3 to 0.05 g/m3? Was this value taken as a constant value regardless of the size of the hydrometeors sampled? How was the estimation of the accuracy of the Nevzorov probe done?
Line 226: If the minimum of 50 μm was used as lower bound then FCDP was not used in this analysis?
Line 229: Consider adding “(The definition of) several bulk…”
Line 232: Add year to Heymsfield et al.
Lines 242-248: I think this paragraph and figure 6 could be removed from the article. Just adding a sentence that states that 22Nov will be analyzed and why should be sufficient.
Line 258: Consider adding “In (current or past) literature…” and adding several references of this literature.
Lines 266-267: Convair average ground speed is 100m/s, this makes the volume for the in situ sampling of 200-500 m. How is this mismatch between the radar and in situ volumes handled?
Lines 271-272: I suggest analyzing if reflectivity shows a different signal comparing Nadir vs Zenith samples and then compare DFR.
Section 3.3.1: This paragraph and associated plot is confusing, the almost 500 m difference between the observations at nadir and zenith show that the vertical structure of the cloud has a large impact on DFR. This is particularly clear when comparing DFR X/Ka, where the Zenith ratio has an almost constant value regardless of the values measured at Nadir. Is there a reason for this? As mentioned before, joint distributions of reflectivity would be better to analyze this 500 m effect.
Line 288: Add ‘s’ to Gans.
Line 294: Figures should be numbered in sequential order.
Lines 292-294: I’d consider regions where both the Nadir and the Zenith correlation coefficient are good. This will hint at a more homogeneous cloud in the vertical and thus a more reliable comparison between what is sample in situ and at a 245 m difference in height. For example, the 4th region marked in figure 10 shows a big difference between Nadir and Zenith (based on correlation coefficient) this is most probably hinting that the part of the cloud sampled has a notable vertical structure so, I’m not sure it is a good case to analyze DFR because this gives an extra reason for the differences between the DFRs.
Lines 296-305: How long was this flight? How many samples were analyzed?
Line 333: Replace ‘similarity measurements‘ by ‘correlation coefficient filter’ or cc threshold or similar.
Line 339: How where the different sections within segment 1, 2 and 4 defined? Was this breakdown into sections defined by the aircraft sampling pattern? Because figure 13 shows that there are different processes occurring in these different sections, for example, section D shows clearly different behaviors in DFR and CPI particle fraction near the beginning of the section when compared to the end of the section.
Lines 340-346: It’d be best if this description of the first segment has figure 13 as reference, it’d help contextualize the differences in the different segments.
Lines 345-346: The bimodality in the PSD distributions it difficult to see for all sections, especially the referred maxima at 1 mm. Please clarify.
Line 351: IWC should be TWC or is the legend in the figure incorrect?
Line 357: Why/How was the aircraft at 6 km height? I assumed the black line in Figure 10a is the aircraft path, so in sections b-d shouldn’t the aircraft be at around 2.4 km height? It would be extremely beneficial to have the different sections shown in figure 10a.
Line 364: Consider rewriting ‘remarkably mirrors’
Line 365: Similar to the comment before, from Figure 10a after the first few minutes of the first segment (first few minutes of section a that the aircraft descended) the aircraft seems to be flying at a constant altitude of ~2.4 km, is this not the case? Aren’t sections A-D correspond to segment 1 that is shown by the first box in figure 10a?
Lines 372-373: This sentence that the fraction of dendrites and rimed particles drops to its lowest level in section D can be misleading, this is the case near the beginning of the section, but by the end of the section this is clearly not the case. Consider rephrasing to avoid confusion.
Lines 390-391: Would it be possible to add to figure 14 the line representing graupel particles using discrete dipole approximation? It could be helpful to add the relevant curves from figure 2.
Line 405: Consider ordering figures in sequential order that are mentioned in the article (Figure 15 should come before Figure 16). Also, how where these different sections defined?
Lines 409-412: From figure 16 top right panel it seems like the fraction of dendrites is higher in C than in B?
Line 419: MDV does not reach 6 mm in section A, please clarify.
Lines 420-421: ZDR was not mentioned before in the article, was there not any ZDR signature in the previous segment?
Line 421-423: Consider rewriting this sentence, variables do not mimic other variables. Also, MDV is not > 8 mm for all the times that DFR is ~ 10 dB, this occurs just for the second maxima in DFR in section B.
Line 423: consider adding a time series of ZDR to figure 16.
Lines 426-427: What do you mean by ‘fluctuations in the DFRs’? Also, what datapoints correspond to section A, B or C is not clear from the plot, consider making the markers edge a different color linked with each section. Also, I suggest adjusting the limits of the plot to better fit the data plotted this will make the differences in the markers size and colors clearer.
Lines 429-435: These few sentences are confusing consider rewriting. The hook feature is present in data from section C not B. Also, remove parenthesis for the references inside the larger parenthesis, like in the Petty and Huang, 2010 reference.
Line 444: Why was this segment chose to be analyzed in depth? Please consider the previous comment regarding the difference between observations at Nadir and Zenith with respect to the vertical structure of the cloud and how different processes can be playing a role at different heights of the clouds that could give difference in DFRs that are not exclusively related to the factors analyzed here. This variability in height can clearly be seen in figure 10a where different microphysics might be acting between the lowest trusted range sampled by the radar and the in situ observations sampled 245 m below.
Lines 484-485: This sentence is misleading as a summary and discussion part of this article where 1 day was analyzed and 3 segments of that flight that resulted in 49 minutes of DFR observations.
Lines 488-490: This sentence is confusing, please consider rewriting it.
Lines 496-499: Add ‘for the flight we analyzed’ or similar phrase here as the results described could be case dependent.
Lines 503-504: Figure 22 should be figure 21?
Line 503: I’m not sure I understand how figure 22 (or 21) is a first attempt for quantitative retrieval of particle size using measured DFR Ka/W and DFR X/Ka. There are a lot mean MVD and density values that are linked with different DFRs values. This looks more like a qualitative analysis.
Line 505: Figure 21b shows that equivalent density increases with decreasing both DFR
Line 510: Remove extra ‘.’
Lines 515-517: Change ‘demonstrated’ by ‘shown’, otherwise, this is too strong of a statement based in the case analyzed in this article.
Lines 523-524: This sentence is not clear, consider rewriting it for clarity.
Fig. 1: please improve figure 1a, it is very difficult to visualize the location of the flights, what is the color in the track representing? Not sure how convoluted a figure showing the path of all the flights done in RadSnowExp, but it’ll sure be useful to see such figure.
Fig. 4: Consider center this figure around the heights that are referenced in the article. For example, if the ground is not referenced why extend the y-axis all the way 5 km?
Fig. 7: Please add in the figure caption what is shown in the figure (what is the dashed line and the error? bars).
Fig. 8: Consider adding lat and lon to the plots and make the maps larger the synoptic map is not very useful in the context given for people not familiarize with the area the flights took place to easily link all 3 maps.
Fig. 9: What is the F9 legend? What do you mean by ‘the ground to air temperature’ ground temperature is ~ -23C and why is this important information to have?
Fig. 10: Improve figure caption to add what are the red, blue and black lines in panel (b) and also improve the legibility of the legend in panels b-d. Also, please consider not using similarity measurement and just use correlation coefficient.
Fig. 12: Please correct what figure is referencing the five sections. Also, I don’t see a reference in the article that requires panel b to be part of the figure.
Fig. 13: Please label each panel of the figure for an easier read and to improve the connection between the text and each plot in the figure.
Fig 14: Consider modifying the colorscheme of the scatter plot to improve the readability of the figure (pinkish colors represents both low and high values).
Fig. 15b: I don’t see a mention to this panel in the article.
Fig. 16: Similar to Figure 13 please label each panel.
Fig. 21: A joint distribution with number of samples would give reference to the mean MVD and density.
Fig. 22: Add what the black line and vertical black lines are in the figure caption.
Citation: https://doi.org/10.5194/amt-2021-148-RC3 - AC4: 'Reply on RC3', Cuong Nguyen, 02 Oct 2021