ALBERT

Description: An experimental retrieval of oceanic warm rain is presented, extending a previous variational algorithm to provide a suite of retrieved variables spanning non-raining through predominantly warm raining conditions. The warm rain retrieval is underpinned by hydrometeor covariances and drizzle onset data derived from CloudSat. Radiative transfer modelling and analysis of drop size variability from disdrometer observations permit state-dependent observation error covariances that scale with columnar rainwater during iteration. The state-dependent errors and nuanced treatment of drop distributions in precipitating regions are novel and may be applicable for future retrievals and all-sky data assimilation methods. This retrieval method can effectively increase passive microwave sensors' sensitivity to light rainfall that might otherwise be missed. Comparisons with space-borne and ground radar estimates are provided as a proof of concept, demonstrating that a passive-only variational retrieval can be sufficiently constrained from non-raining through warm rain conditions. Significant deviations from forward model assumptions cause non-convergence, usually a result of scattering hydrometeors above the freezing level. However, for cases with liquid-only precipitation, this retrieval displays greater sensitivity than a benchmark operational retrieval. Analysis against passive and active products from the Global Precipitation Measurement (GPM) satellite shows substantial discrepancies in precipitation frequency, with the experimental retrieval observing more frequent light rain. This approach may be complementary to other precipitation retrievals, and its potential synergy with the operational passive GPM retrieval is briefly explored. There are also implications for data assimilation, as all 13 channels on the GPM Microwave Imager (GMI) are simulated over ocean with fidelity in warm raining conditions.

Description: This study assesses the global distribution of mean atmospheric ice mass from current state-of-the-art estimates and its variability on daily and seasonal timescales. Ice water path (IWP) retrievals from active and passive satellite platforms are analysed and compared with estimates from two reanalysis data sets, ERA5 (European Centre for Medium-range Weather Forecasts Reanalysis 5, ECMWF) and MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications 2). Large discrepancies in IWP exist between the satellite data sets themselves, making validation of the model results problematic and indicating that progress towards a consensus on the distribution of atmospheric ice has been limited. Comparing the data sets, zonal means of IWP exhibit similar shapes but differing magnitudes, with large IWP values causing much of the difference in means. Diurnal analysis centred on A-Train overpasses shows similar structures in some regions, but the degree and sign of the variability varies widely; the reanalyses exhibit noisier and higher-amplitude diurnal variability than borne out by the satellite estimates. Spatial structures governed by the atmospheric general circulation are fairly consistent across the data sets, as principal component analysis shows that the patterns of seasonal variability line up well between the data sets but disagree in severity. These results underscore the limitations of the current Earth observing system with respect to atmospheric ice, as the level of consensus between observations is mixed. The large-scale variability of IWP is relatively consistent, whereas disagreements on diurnal variability and global means point to varying microphysical assumptions in retrievals and models alike that seem to underlie the biggest differences.

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 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 study assesses the global distribution of mean atmospheric ice mass from current state-of-the-art estimates and its variability on daily and seasonal timescales. Ice water path (IWP) retrievals from active and passive satellite platforms are compared and analysed against estimates from two reanalysis datasets, ERA5 (European Centre for Medium-range Weather Forecasts Reanalysis 5) and MERRA-2 (Modern-era Retrospective Reanalysis for Research and Applications 2). Large discrepancies in IWP exist between the satellite datasets themselves, making validation of the model results problematic and indicating that progress towards consensus on the distribution of atmospheric ice has been limited. Comparing the datasets, zonal means of IWP exhibit similar shapes but differing magnitudes. Diurnal analysis centred on A-Train overpasses shows homologous structures in some regions, but the degree and sign of the variability varies widely; the reanalyses exhibit noisier and higher amplitude diurnal variability than borne out by the satellite estimates. Spatial structures governed by the atmospheric general circulation are fairly consistent across the datasets, as principal component analysis shows that the patterns of seasonal variability line up well between the datasets but disagree in severity. These results underscore the limitations of the current Earth observing system with respect to atmospheric ice, as the level of consensus between observations is mixed. The large-scale variability of IWP is relatively consistent, whereas disagreements on diurnal variability and global means point to varying microphysical assumptions in retrievals and models alike that seem to underlie the biggest differences.

Description: An experimental retrieval of oceanic warm rain is presented, extending a previous variational algorithm to provide a suite of retrieved variables spanning non-raining through predominantly warm raining conditions. The warm rain retrieval is underpinned by hydrometeor covariances and drizzle onset data derived from CloudSat. Radiative transfer modelling and analysis of drop size variability from disdrometer observations permit state-dependent observation error covariances that scale with columnar rainwater during iteration. The state-dependent errors and nuanced treatment of drop distributions in precipitating regions are novel and may be applicable for future retrievals and all-sky data assimilation methods. This retrieval method can effectively increase passive microwave sensors' sensitivity to light rainfall that might otherwise be missed. Comparisons with space-borne and ground radar estimates are provided as a proof of concept, demonstrating that a passive-only variational retrieval can be sufficiently constrained from non-raining through warm rain conditions. Significant deviations from forward model assumptions cause non-convergence, usually a result of scattering hydrometeors above the freezing level. However, for cases with liquid-only precipitation, this retrieval displays greater sensitivity than a benchmark operational retrieval. Analysis against passive and active products from the Global Precipitation Measurement (GPM) satellite shows substantial discrepancies in precipitation frequency, with the experimental retrieval observing more frequent light rain. This approach may be complementary to other precipitation retrievals, and its potential synergy with the operational passive GPM retrieval is briefly explored. There are also implications for data assimilation, as all 13 channels on the GPM Microwave Imager (GMI) are simulated over ocean with fidelity in warm raining conditions.

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: A methodology based on quantile regression neural networks (QRNNs) is presented that identifies and corrects the cloud impact on microwave humidity sounder radiances at 183 GHz. This approach estimates the posterior distributions of noise-free clear-sky (NFCS) radiances, providing nearly bias-free estimates of clear-sky radiances with a full posterior error distribution. It is first demonstrated by application to a present sensor, the MicroWave Humidity Sounder 2 (MWHS-2); then the applicability to sub-millimetre (sub-mm) sensors is also analysed. The QRNN results improve upon what operational cloud filtering techniques like a scattering index can achieve but are ultimately imperfect due to limited information content on cirrus impact from traditional microwave channels – the negative departures associated with high cloud impact are successfully corrected, but thin cirrus clouds cannot be fully corrected. In contrast, when sub-mm observations are used, QRNN successfully corrects most cases with cloud impact, with only 2 %–6 % of the cases left partially corrected. The methodology works well even if only one sub-mm channel (325 GHz) is available. When using sub-mm observations, cloud correction usually results in error distributions with a standard deviation less than typical channel noise values. Furthermore, QRNN outputs predicted quantiles for case-specific uncertainty estimates, successfully representing the uncertainty of cloud correction for each observation individually. In comparison to deterministic correction or filtering approaches, the corrected radiances and attendant uncertainty estimates have great potential to be used efficiently in assimilation systems due to being largely unbiased and adding little further uncertainty to the measurements.