ALBERT

Description: Remote sensing observations at sub-millimeter wavelengths provide higher sensitivity to small hydrometeors and low water content than observations at millimeter wavelengths, which are traditionally used to observe clouds and precipitation. They are employed increasingly in field campaigns to study cloud microphysics and will be integrated into the global meteorological observing system to measure the global distribution of ice in the atmosphere with the launch of the Ice Cloud Imager (ICI) radiometer on board the second generation of European operational meteorological satellites (Metop-SG). Observations at these novel wavelengths provide valuable information not only on their own but also in combination with complementary observations at other wavelengths. This study investigates the potential of combining passive sub-millimeter radiometer observations with a hypothetical W-band cloud radar for the retrieval of frozen hydrometeors. An idealized cloud model is used to investigate the information content of the combined observations and establish their capacity to constrain the microphysical properties of ice hydrometeors. A synergistic retrieval algorithm for airborne observations is proposed and applied to simulated observations from a cloud-resolving model. Results from the synergistic retrieval are compared to equivalent radar- and passive-only implementations in order to assess the benefits of the synergistic sensor configuration. The impact of the assumed ice particle shape on the retrieval results is assessed for all retrieval implementations. We find that the combined observations better constrain the microphysical properties of ice hydrometeors, which reduces uncertainties in retrieved ice water content and particle number concentrations for suitable choices of the ice particle model. Analysis of the retrieval information content shows that, although the radar contributes the largest part of the information in the combined retrieval, the radiometer observations provide complementary information over a wide range of atmospheric states. Furthermore, the combined observations yield slightly improved retrievals of liquid cloud water in mixed-phase clouds, pointing towards another potential application of combined radar–radiometer observations.

Description: Representation of the drop size distribution (DSD) of rainfall is a key element of characterizing precipitation in models and observations, with a functional form necessary to calculate the precipitation flux and the drops' interaction with radiation. With newly available oceanic disdrometer measurements, this study investigates the validity of commonly used DSDs, potentially useful a priori constraints for retrievals, and the impacts of DSD variability on radiative transfer. These data are also compared with leading satellite-based estimates over ocean, with the disdrometers observing a larger number of small drops and significantly more variability in number concentrations. This indicates that previous appraisals of raindrop variability over ocean may have been underestimates. Forward model errors due to DSD variability are shown to be significant for both active and passive sensors. The modified gamma distribution is found to be generally adequate to describe rain DSDs but may cause systematic errors for high-latitude or stratocumulus rain retrievals. Depending on the application, an exponential or generalized gamma function may be preferable for representing oceanic DSDs. An unsupervised classification algorithm finds a variety of DSD shapes that differ from commonly used DSDs but does not find a singular set that best describes the global variability.

Description: A neural-network-based method, quantile regression neural networks (QRNNs), is proposed as a novel approach to estimating the a posteriori distribution of Bayesian remote sensing retrievals. The advantage of QRNNs over conventional neural network retrievals is that they learn to predict not only a single retrieval value but also the associated, case-specific uncertainties. In this study, the retrieval performance of QRNNs is characterized and compared to that of other state-of-the-art retrieval methods. A synthetic retrieval scenario is presented and used as a validation case for the application of QRNNs to Bayesian retrieval problems. The QRNN retrieval performance is evaluated against Markov chain Monte Carlo simulation and another Bayesian method based on Monte Carlo integration over a retrieval database. The scenario is also used to investigate how different hyperparameter configurations and training set sizes affect the retrieval performance. In the second part of the study, QRNNs are applied to the retrieval of cloud top pressure from observations by the Moderate Resolution Imaging Spectroradiometer (MODIS). It is shown that QRNNs are not only capable of achieving similar accuracy to standard neural network retrievals but also provide statistically consistent uncertainty estimates for non-Gaussian retrieval errors. The results presented in this work show that QRNNs are able to combine the flexibility and computational efficiency of the machine learning approach with the theoretically sound handling of uncertainties of the Bayesian framework. Together with this article, a Python implementation of QRNNs is released through a public repository to make the method available to the scientific community.

Description: A two-dimensional variational retrieval (2D-Var) is presented for a passive microwave imager. The overlapping antenna patterns of all frequencies from the Advanced Microwave Scanning Radiometer 2 (AMSR2) are explicitly simulated to attempt retrieval of near-surface wind speed and surface skin temperature at finer spatial scales than individual antenna beams. This is achieved, with the effective spatial resolution of retrieved parameters judged by analysis of 2D-Var averaging kernels. Sea surface temperature retrievals achieve about 30 km resolution, with wind speed retrievals at about 10 km resolution. It is argued that multi-dimensional optimal estimation permits greater use of total information content from microwave sensors than other methods, with no compromises on target resolution needed; instead, various targets are retrieved at the highest possible spatial resolution, driven by the channels' sensitivities. All AMSR2 channels can be simulated within near their published noise characteristics for observed clear-sky scenes, though calibration and emissivity model errors are key challenges. This experimental retrieval shows the feasibility of 2D-Var for cloud-free retrievals and opens the possibility of stand-alone 3D-Var retrievals of water vapour and hydrometeor fields from microwave imagers in the future. The results have implications for future satellite missions and sensor design, as spatial oversampling can somewhat mitigate the need for larger antennas in the push for higher spatial resolution.

Description: Representation of the drop size distribution (DSD) of rainfall is a key element of characterizing precipitation in models and retrievals, with a functional form necessary to calculate the precipitation flux and the drops' interaction with radiation. With newly available oceanic disdrometer measurements, this study investigates the validity of commonly used DSDs, potentially useful a priori constraints for retrievals, and the forward model errors caused by DSD variability. These data are also compared to leading satellite-based estimates of oceanic DSDs. Forward model errors due to DSD variability are shown to be significant for both active and passive sensors. The modified gamma distribution is found to be generally adequate to describe rain DSDs, but may cause systematic errors for high latitude or stratocumulus rain retrievals; depending on the application, an exponential or generalized gamma function may be preferable for representing oceanic DSDs. An unsupervised classification algorithm finds a variety of DSD shapes that differ from commonly used DSDs, but does not find a singular set that best describes the global variability. Finally, DSD shapes are found to be not particularly distinctive of regional or large-scale environments, but rather occur at varying frequencies over the global oceans.

Description: This work is concerned with the retrieval of physical quantities from remote sensing measurements. A neural network based method, Quantile Regression Neural Networks (QRNNs), is proposed as a novel approach to estimate the a posteriori distribution of Bayesian remote sensing retrievals. The advantage of QRNNs over conventional neural network retrievals is that they not only learn to predict a single retrieval value but also the associated, case specific uncertainties. In this study, the retrieval performance of QRNNs is characterized and compared to that of other state-of-the-art retrieval methods. A synthetic retrieval scenario is presented and used as a validation case for the application of QRNNs to Bayesian retrieval problems. The QRNN retrieval performance is evaluated against Markov chain Monte Carlo simulation and another Bayesian method based on Monte Carlo integration over a retrieval database. The scenario is also used to investigate how different hyperparameter configurations and training set sizes affect the retrieval performance. In the second part of the study, QRNNs are applied to the retrieval of cloud top pressure from observations by the moderate resolution imaging spectroradiometer (MODIS). It is shown that QRNNs are not only capable of achieving similar accuracy as standard neural network retrievals, but also provide statistically consistent uncertainty estimates for non-Gaussian retrieval errors. The results presented in this work show that QRNNs are able to combine the flexibility and computational efficiency of the machine learning approach with the theoretically sound handling of uncertainties of the Bayesian framework. Together with this article, a Python implementation of QRNNs is released through a public repository to make the method available to the scientific community.

Description: A two-dimensional variational retrieval (2DVAR) is presented for a passive microwave imager. The overlapping antenna patterns of all frequencies from the Advanced Microwave Scanning Radiometer-2 (AMSR2) are explicitly simulated to attempt retrieval of near surface wind speed and surface skin temperature at finer spatial scales than individual antenna beams. This is achieved, with the effective spatial resolution of retrieved parameters shown by analysis of 2DVAR averaging kernels. Sea surface temperature retrievals achieve about 30 km resolution, with wind speed retrievals at about 10 km resolution. It is argued that multi-dimensional optimal estimation permits greater use of total information content from microwave sensors than other methods, with no compromises on target resolution needed; instead, various targets are retrieved at the highest possible spatial resolution, driven by the channels' sensitivities. All AMSR2 channels can be simulated within near their published noise characteristics for observed clear-sky scenes, though calibration and emissivity model errors are key challenges. This experimental retrieval shows the feasibility of 2DVAR for cloud-free retrievals, and opens the possibility of standalone 3DVAR retrievals of water vapour and hydrometeor fields from microwave imagers in the future. The results have implications for future satellite missions and sensor design, as spatial oversampling can somewhat mitigate the need for larger antennas in the push for higher spatial resolution.

Description: Satellite microwave remote sensing is an important tool for determining the distribution of atmospheric ice globally. The upcoming Ice Cloud Imager (ICI) sensor will provide unprecedented measurements at sub-millimetre frequencies, employing channels up to 664 GHz. However, the utilization of such measurements requires detailed data on how individual ice particles scatter and absorb radiation, i.e., single scattering data. Several single scattering databases are currently available, with the one by Eriksson et al. (2018) specifically tailored to ICI. This study attempts to validate and constrain the large set of particle models available in this database, to a smaller and more manageable set. A combined active and passive model framework is developed and employed, which converts CloudSat observations to simulated brightness temperatures (TBs) measured by the GPM Microwave Imager (GMI) and ICI. Simulations covering about one month in the tropic pacific ocean are performed, assuming different microphysical settings realized as combinations of particle model and particle size distribution (PSD). Firstly, it is found that when the CloudSat inversions and passive forward model are considered separately, assumed particle model and PSD have a considerable impact on both radar retrieved ice water content (IWC) and simulated TBs. Conversely, when the combined active and passive framework is employed instead, the uncertainty due to assumed particle model is significantly reduced, essentially due to a compensation effect between bulk extinction at passive frequencies and radar reflectivity. Furthermore, simulated TBs for almost all the tested microphysical combinations, from a statistical point of view, agree well with GMI measurements (186.31 and 190.31 GHz), indicating the robustness of the simulations. However, it is difficult to identify a particle model that outperforms any other. One aggregate particle model, composed of columns, yields marginally better agreement to GMI compared to the other particles, mainly for the most severe cases of deep convection. Of tested PSDs, the one by McFarquhar and Heymsfield (1997) is found to give the best overall agreement to GMI and also yields radar dBZ-IWC relationships closely matching measurements by Protat et al. (2016). Only one particle, modelled as an air-ice mixture spheroid, performs poorly overall. On the other hand, simulations at the higher ICI frequencies (328.65, 334.65, and 668.2 GHz) show significantly higher sensitivity to the assumed particle model. This study thus points to the potential use of combined ICI and 94 GHz radar measurements to constrain ice hydrometeor properties in radiative transfer (RT), using the method demonstrated in this paper.

Description: Satellite microwave remote sensing is an important tool for determining the distribution of atmospheric ice globally. The upcoming Ice Cloud Imager (ICI) will provide unprecedented measurements at sub-millimetre frequencies, employing channels up to 664 GHz. However, the utilization of such measurements requires detailed data on how individual ice particles scatter and absorb radiation, i.e. single scattering data. Several single scattering databases are currently available, with the one by Eriksson et al. (2018) specifically tailored to ICI. This study attempts to validate and constrain the large set of particle models available in this database to a smaller and more manageable set. A combined active and passive model framework is developed and employed, which converts CloudSat observations to simulated brightness temperatures (TBs) measured by the Global Precipitation Measurement (GPM) Microwave Imager (GMI) and ICI. Simulations covering about 1 month in the tropical Pacific Ocean are performed, assuming different microphysical settings realized as combinations of the particle model and particle size distribution (PSD). Firstly, it is found that when the CloudSat inversions and the passive forward model are considered separately, the assumed particle model and PSD have a considerable impact on both radar-retrieved ice water content (IWC) and simulated TBs. Conversely, when the combined active and passive framework is employed instead, the uncertainty due to the assumed particle model is significantly reduced. Furthermore, simulated TBs for almost all the tested microphysical combinations, from a statistical point of view, agree well with GMI measurements (166, 186.31, and 190.31 GHz), indicating the robustness of the simulations. However, it is difficult to identify a particle model that outperforms any other. One aggregate particle model, composed of columns, yields marginally better agreement with GMI compared to the other particles, mainly for the most severe cases of deep convection. Of the tested PSDs, the one by McFarquhar and Heymsfield (1997) is found to give the best overall agreement with GMI and also yields radar dBZ–IWC relationships closely matching measurements by Protat et al. (2016). Only one particle, modelled as an air–ice mixture spheroid, performs poorly overall. On the other hand, simulations at the higher ICI frequencies (328.65, 334.65, and 668.2 GHz) show significantly higher sensitivity to the assumed particle model. This study thus points to the potential use of combined ICI and 94 GHz radar measurements to constrain ice hydrometeor properties in radiative transfer (RT) using the method demonstrated in this paper.

Description: Atmospheric tides are important for vertical coupling in the atmosphere, from the stratosphere down to the troposphere and up to the thermosphere. They are planetary-scale gravity waves with well-known periods that are integer fractions of a day and can be observed in the temperature or wind fields in the atmosphere. Most lidar techniques and satellites measure atmospheric tides only in the temperature field and continuous measurements of the tides in the wind field of the stratosphere and lower mesosphere are rare, even though, with modern lidars, they would be feasible. In this study, we present measurements of the diurnal tide in the wind field in the stratosphere and lower mesosphere by ground-based microwave wind radiometry for two different campaigns in tropical and polar regions. Further, we compare our measurements to MERRA-2 reanalysis data. In the tri-monthly mean, we find a good correspondence in the amplitude and phase of the diurnal tide between measurements and reanalysis with the most important features of the diurnal tides represented in both data sets. When looking at shorter timescales, we find significant differences in the data sets. We make an attempt to examine these differences and discriminate between atmospheric variability and noise, and we present some hints for intermittent diurnal tides. We conclude that continuous ground-based observations of tides in the middle atmospheric wind field are feasible, and they deliver consistent results for the mean amplitude and phase of the diurnal tide in the tri-monthly mean. We further discuss the limitations in regards to short timescale observations of tides and the possibility to provide additional insight into middle atmospheric dynamics that is complementary to temperature observations and reanalysis data.