GPROF-NN: a neural-network-based implementation of the Goddard Profiling Algorithm

The Global Precipitation Measurement (GPM) mission measures global precipitation at a temporal resolution of a few hours to enable close monitoring of the global hydrological cycle. GPM achieves this by combining observations from a spaceborne precipitation radar, a constellation of passive microwave (PMW) sensors, and geostationary satellites. The Goddard Profiling Algorithm (GPROF) is used operationally to retrieve precipitation from all PMW sensors of the GPM constellation. Since the resulting precipitation rates serve as input for many of the level 3 retrieval products, GPROF constitutes an essential component of the GPM processing pipeline. This study investigates ways to improve GPROF using modern machine learning methods. We present two neuralnetwork-based, probabilistic implementations of GPROF: GPROF-NN 1D, which (just like the current GPROF implementation) processes pixels individually, and GPROF-NN 3D, which employs a convolutional neural network to incorporate structural information into the retrieval. The accuracy of the retrievals is evaluated using a test dataset consistent with the data used in the development of the GPROF and GPROF-NN retrievals. This allows for assessing the accuracy of the retrieval method isolated from the representativeness of the training data, which remains a major source of uncertainty in the development of precipitation retrievals. Despite using the same input information as GPROF, the GPROF-NN 1D retrieval improves the accuracy of the retrieved surface precipitation for the GPM Microwave Imager (GMI) from 0.079 to 0.059mmh 1 in terms of mean abso- lute error (MAE), from 76.1% to 69.5% in terms of symmetric mean absolute percentage error (SMAPE) and from 0 :797 to 0 :847 in terms of correlation. The improvements for the Microwave Humidity Sounder (MHS) are from 0.085 to 0.061mmh 1 in terms of MAE, from 81% to 70.1% for SMAPE, and from 0 :724 to 0 :804 in terms of correlation. Comparable improvements are found for the retrieved hydrometeor profiles and their column integrals, as well as the detection of precipitation. Moreover, the ability of the retrievals to resolve small-scale variability is improved by more than 40% for GMI and 29% for MHS. The GPROFNN 3D retrieval further improves the MAE to 0.043mmh 1; the SMAPE to 48.67 %; and the correlation to 0:897 for GMI and 0.043mmh 1, 63.42 %, and 0:83 for MHS. Application of the retrievals to GMI observations of Hurricane Harvey shows moderate improvements when compared to co-located GPM-combined and ground-based radar measurements indicating that the improvements at least partially carry over to assessment against independent measurements. Similar retrievals for MHS do not show equally clear improvements, leaving the validation against independent measurements for future investigation. Both GPROF-NN algorithms make use of the same input and output data as the original GPROF algorithm and thus may replace the current implementation in a future update of the GPM processing pipeline. Despite their superior accuracy, the single-core runtime required for the operational processing of an orbit of observations is lower than that of GPROF. The GPROF-NN algorithms promise to be a simple and cost-efficient way to increase the accuracy of the PMW precipitation retrievals of the GPM constellation and thus improve the monitoring of the global hydrological cycle

Microwave Brightness Temperatures of Tilted Convective Systems

Aircraft and ground-based radar data from the Tropical Ocean and Global Atmosphere Coupled-Ocean Atmosphere Response Experiment (TOGA COARE) show that convective systems are not always vertical. Instead, many are tilted from vertical. Satellite passive microwave radiometers observe the atmosphere at a viewing angle. For example, the Special Sensor Microwave/Imager (SSM/I) on Defense Meteorological Satellite Program (DMSP) satellites and the Tropical Rainfall Measurement Mission (TRMM) Microwave Imager (TMI) on the TRMM satellite have an incident angle of about 50deg. Thus, the brightness temperature measured from one direction of tilt may be different than that viewed from the opposite direction due to the different optical depth. This paper presents the investigation of passive microwave brightness temperatures of tilted convective systems. To account for the effect of tilt, a 3-D backward Monte Carlo radiative transfer model has been applied to a simple tilted cloud model and a dynamically evolving cloud model to derive the brightness temperature. The radiative transfer results indicate that brightness temperature varies when the viewing angle changes because of the different optical depth. The tilt increases the displacements between high 19 GHz brightness temperature (Tb(sub 19)) due to liquid emission from lower level of cloud and the low 85 GHz brightness temperature (Tb(sub 85)) due to ice scattering from upper level of cloud. As the resolution degrades, the difference of brightness temperature due to the change of viewing angle decreases dramatically. The dislocation between Tb(sub 19) and Tb(sub 85), however, remains prominent

Global Precipitation Measurement

This chapter begins with a brief history and background of microwave precipitation sensors, with a discussion of the sensitivity of both passive and active instruments, to trace the evolution of satellite-based rainfall techniques from an era of inference to an era of physical measurement. Next, the highly successful Tropical Rainfall Measuring Mission will be described, followed by the goals and plans for the Global Precipitation Measurement (GPM) Mission and the status of precipitation retrieval algorithm development. The chapter concludes with a summary of the need for space-based precipitation measurement, current technological capabilities, near-term algorithm advancements and anticipated new sciences and societal benefits in the GPM era

Towards variational retrieval of warm rain from passive microwave observations

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\u27 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

Combined Radar and Radiometer Analysis of Precipitation Profiles for a Parametric Retrieval Algorithm

A methodology to analyze precipitation profiles using the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) and precipitation radar (PR) is proposed. Rainfall profiles are retrieved from PR measurements, defined as the best-fit solution selected from precalculated profiles by cloud-resolving models (CRMs), under explicitly defined assumptions of drop size distribution (DSD) and ice hydrometeor models. The PR path-integrated attenuation (PIA), where available, is further used to adjust DSD in a manner that is similar to the PR operational algorithm. Combined with the TMI-retrieved nonraining geophysical parameters, the three-dimensional structure of the geophysical parameters is obtained across the satellite-observed domains. Microwave brightness temperatures are then computed for a comparison with TMI observations to examine if the radar-retrieved rainfall is consistent in the radiometric measurement space. The inconsistency in microwave brightness temperatures is reduced by iterating the retrieval procedure with updated assumptions of the DSD and ice-density models. The proposed methodology is expected to refine the a priori rain profile database and error models for use by parametric passive microwave algorithms, aimed at the Global Precipitation Measurement (GPM) mission, as well as a future TRMM algorithms

Long-Term Arctic Snow/Ice Interface Temperature from Special Sensor for Microwave Imager Measurements

The Arctic sea ice region is the most visible area experiencing global warming-induced climate change. However, long-term measurements of climate-related variables have been limited to a small number of variables such as the sea ice concentration, extent, and area. In this study, we attempt to produce a long-term temperature record for the Arctic sea ice region using Special Sensor for Microwave Imager (SSM/I) Fundamental Climate Data Record (FCDR) data. For that, we developed an algorithm to retrieve the wintertime snow/ice interface temperature (SIIT) over the Arctic Ocean by counting the effect of the snow/ice volume scattering and ice surface roughness on the apparent emissivity (the total effect is referred to as the correction factor). A regression equation was devised to predict the correction factor from SSM/I brightness temperatures (TBs) only and then applied to SSM/I 19.4 GHz TB to estimate the SIIT. The obtained temperatures were validated against collocated Cold Regions Research and Engineering Laboratory (CRREL) ice mass balance (IMB) drifting buoy-measured temperatures at zero ice depth. It is shown that the SSM/I retrievals are in good agreement with the drifting buoy measurements, with a correlation coefficient of 0.95, bias of 0.1 K, and root-mean-square error of 1.48 K on a daily time scale. By applying the algorithm to 24-year (1988⁻2011) SSM/I FCDR data, we were able to produce the winter-time temperature at the sea ice surface for the 24-year period