11 Mar 2021
11 Mar 2021
Physical characteristics of frozen hydrometeors inferred with parameter estimation
- ECMWF, Shinfield Park, Reading, RG2 9AX, UK
- ECMWF, Shinfield Park, Reading, RG2 9AX, UK
Abstract. Frozen hydrometeors are found in a huge range of shapes and sizes, with variability on much smaller scales than those of typical models or satellite observations. Neither models nor in-situ measurements can fully describe this variability, so assumptions have to be made in applications including atmospheric modelling and radiative transfer. Here parameter estimation has been used to optimise six different assumptions relevant to frozen hydrometeors in passive microwave radiative transfer. This covers cloud overlap, convective water content and particle size distribution (PSD), the shapes of large-scale snow and convective snow and an initial exploration of the ice cloud representation (particle shape and PSD combined). These parameters were simultaneously adjusted to find the best fit between simulations from the European Centre for Medium-range Weather Forecasts (ECMWF) assimilation system, and near-global microwave observations covering the frequency range 19 GHz to 190 GHz. The choices for the cloud overlap and the convective particle shape were particularly well constrained (or identifiable) and there was even constraint on the cloud ice PSD. The practical output is a set of improved assumptions to be used in version 13 of the Radiative Transfer for TOVS microwave scattering package (RTTOV-SCATT), taking into account newly available particle shapes such as aggregates and hail, and additional PSD options. The parameter estimation explored the full parameter space using an efficient assumption of linearly additive perturbations. This helped illustrate issues such as multiple minima in the cost function, and non-Gaussian errors, that would make it hard to implement the same approach in a standard data assimilation system for weather forecasting. Nevertheless, as modelling systems grow more complex, parameter estimation is likely to be a necessary part of the development process.
Alan J. Geer
Status: open (until 06 May 2021)
-
RC1: 'Comment on amt-2021-50', Anonymous Referee #1, 08 Apr 2021
reply
This study develops a physically consistent means of combining NWP model simulations with a simplified radiative transfer (RT) model to estimate both the micro- and macrophysical properties of frozen atmospheric particles with the ultimate goal of improving physical constraints in operational forecasts of clouds and precipitation. I have felt for some time that such an approach would be useful in data assimilation where these physical parameters would be part of the state vector solution. However, the author has done a good job of highlighting and explaining some of the significant challenges of using parameter estimation within a DA system.
One of my main concerns with using coarse-resolution model simulations and simplified RT models is that the parameter estimation will simply reflect the NWP model biases (such as caused by sub-grid-scale parameterizations) and uncertainties in the RT model. In addition, the parameters will require continual updating based on changes to the model parameterizations and complexity of the RT model. However, the author has fully addressed these concerns in section 6.
An important conclusion of this study is the importance of utilizing a wide range of microwave frequencies to constrain the strong frequency dependence of single-scattering properties, which are sensitive to the particle size distribution, mass-size relationship, and particle density and shape. However, the author makes clear that “the particle habit itself is not being precisely identified, but rather the full ‘hydrometeor model’ which is controlled by the choice of particle habit in the current framework but includes also an assumed mass-size relation.” This is an important clarification because there is often a misconception that particle habit is solely responsible for reducing biases in forward RT calculations at high microwave frequencies when compared to observations.
Finally, this study not only discusses the limitations of the method but also offers a way for improvement and describes how it builds on earlier work.
My only complaint is that the paper is rather long. The concern is that readers may be less inclined to read the paper in its entirety and may overlook some of the nuances and important results of the study.
In summary, the author has done an excellent job of articulating a complex problem in simple, understandable terms and in demonstrating the general applicability of the method to other NWP models and RT models.
Minor points:
Line 53: Typo: “There is also a need …”
Line 230-231: Actually, the South Dakota School of Mines and Technology T-28 storm-penetrating aircraft has, for decades, flown through convective cores to take in situ measurements of hail and graupel (e.g., Detwiler et al., 2012, Bull. American Meteorol. Soc.; Field et al., 2019, J. Appl. Meteorol. Climatol., 58, 231-245).
Line 309: “approx.” should be spelled out.
Line 377: Typo: “shape of the cost function is very different”
Line 582: Typo: “If they had done that …”
Line 624: Typo: “it is important to …”
Line 635: Typo: “a little bit like …”
Line 723: Typo: “in the future …”
-
RC2: 'Comment on amt-2021-50', Anonymous Referee #2, 09 Apr 2021
reply
Review for AMT of “Physical characteristics of frozen hydrometeors inferred with parameter estimation” by Alan Geer.
Paper Summary:
Our understanding of ice clouds and ice particle characteristics globally, and especially in deep convection, is sorely lacking. The author aims to address aspects of this problem via perturbation of six parameters relevant to frozen hydrometeors in the ECMWF assimilation system, with the goal of simultaneously adjusting parameters in such a way that simulated microwave radiances agree with observed microwave radiances from F17 SSMIS (for 13 channels, extending from 19 – 183 GHz) over a 13 – 22 June 2019 test period. An output of this study is improved frozen hydrometeor assumptions to be used in version 13 of RTTOV-SCATT.
Major Issue(s) and Comments:
One major issue I had was the lack of discussion on mixed phase microphysical processes and presence of super-cooled liquid in the “frozen” tops of convective cores. A neglect of this implies that ice particle parameter settings will likely have unphysical settings or values due to scattering signatures in convective regions having to be matched by varying only ice parameter degrees of freedom. If variations in liquid (at the expense of ice) were allowed, simulated microwave Tb signatures would be substantially different since liquid would cause emission instead of scattering, thus impacting the choice of ice parameter settings (and potentially allowing graupel to occur much more frequently). We can comfortably neglect liquid at cold icy temperatures in the weak ascent regions of stratiform clouds, but since convective ice is a particular focus of this paper (which is certainly welcomed given the lack of convective ice papers), then at the least, some discussion of the role of liquid at cold temperatures in convection must be discussed. Clearly some of the results for perturbing convective-ice related parameters signify a compensation for lack of liquid allowed in convection, and thus, parameter values decidedly returned are also likely not realistic due to this neglect. The number of studies documenting the prevalence of liquid down to minus 37.5 deg C in convection are numerous and increasing (e.g., in general, independent of aerosol effects or discussion: Kumjian et al. (2012; JAS); Dolan et al. (2013, JAMC); Xu and Zipser (2015, JGR); van Lier-Walqui et al. (2016, MWR); Fuchs et al. (2018, JGR), and within the context of aerosols: Rosenfeld and Woodley (2000; Nature); Fan et al. (2018; Science)). In short, importantly, liquid/graupel/hail is a microphysical phenomenon common in the cold tops of convection, further evidenced by the commonality of lighting in convection globally. I would not be terribly surprised if the lack of choices regarding choosing graupel/hail in this study (and others) relate to our underestimation of liquid in the upper reaches of convection. At the very least, this should be given distinct discussion within the context of the parameter settings determined.
Minor Issue(s) and Comments:
- Is noise in observed microwave radiances considered at any point?
- When it is said that mixing ratios are halved/doubled, does this mean for all altitudes equally, from above the melting level to cloud top?
- The 20K discrepancies in Fig. 4 are sort of glossed over mostly due to occurrence frequencies, but I suspect they are the convective regions (convection by frequency of occurrence is of course small relative to ice cloud in general, so this seems plausible). Would a map of the large discrepancies be correlated with deep convection patterns or no?
- Related to the issue of large fractions of liquid in convective cores, if something so simple as an effective cloud fraction (C ) can be used as a parameter, could a simple factor that governs the amount of convective condensate partitioned into ice vs liquid (with the latter describe bed using a M-P DSD for simplicity) be incorporated quite easily as well? Or, even, couldn’t the temperature for freezing in convection simply be modified to vary from 0 deg C to approx.. minus 37 deg C?
- Abstract, first sentence, Line 2: I recommend replacing ‘models or satellite observations’ with ‘model grid boxes or satellite observation fields-of-view’.
-
RC3: 'Comment on amt-2021-50', Anonymous Referee #3, 13 Apr 2021
reply
This article describes a major upgrade and its physical justification of the frozen hydrometeor parameters in the radiative transfer model of RTTOV relevant to microwave and submillimeter wave spectrum. In the macro-scale, this upgrade revised a cloud overlap scheme and introduced a new “convective snow” that is analogous to the cumulus parameterization in non-cloud-resoling models. In the micro-scale, new particle size distribution and ice/snow particle shapes are introduced. In the parameter selection procedure, the cost function is configured to match the PDF between calculated and observed brightness temperatures, which is innovative. I found that the author has done a thorough examination to this parameter selection problem and worked his best in this inter-tangled web of uncertain parameters to find the best configuration for practical use. The effort is commendable. The structure of the paper is logical, and the writing is easy to follow (although a little lengthy). I recommend accept after some minor revisions. The following are my concerns.
- Although the new (Final) configuration seems to be more physically sound, as a reader, I am not clear how the cost decrease from 0.932 to 0.911 translates to real improvement. I believe that Figure 5 is intended to show the impact, but neither the figure nor the explanation make me feel clearer about the impact. Therefore, I’d like to suggest: (1) improve Figure 5, so that it can show the result from the “Best” is closer to “Observations” than “Control”, (2) add more explanation to argue the “Best” is better than “Control”, and (3) add a PDF figure to show the PDF shifted more closer to “Observations” when “Control” is replaced by “Best”. The author may also try some other ideas. The point is: how to translate the cost decrease to easy-to-understand improvements.
- The “snow mixing ratio” dimension is peculiar. As I understood it, it is rooted from the uncertainty in particles’ falling speed when you convert from mass flux to water content. But still, altering mixing ratio doesn’t seem to be a sound approach. Fortunately, in the “Final” configuration it is decided not to change it. I’d like to propose to eliminate this dimension in the first place.
- The situation-dependent results listed in Table 6 are interesting and make us think more about the physics in terms of how to select ice particles. Clearly, the “snow” particles in the tropics and higher latitudes are different – all the results seem to point out that tropical snow particles are denser than those in mid-latitudes. Or, maybe the “convective” and “large-scale” separation of clouds is not an ideal one after all. Microphysically, the anvil associated with tropical convections are distinctively different from the stratiform clouds associated with frontal large-scale uplifting. This may be an area of future improvement of the “best configuration”. Adding some statements along this line in the conclusion section will be helpful.
Alan J. Geer
Alan J. Geer
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
185 | 65 | 8 | 258 | 5 | 2 |
- HTML: 185
- PDF: 65
- XML: 8
- Total: 258
- BibTeX: 5
- EndNote: 2
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1