Simulation results for smooth water surface - cloudy atmosphere system

A computational experiment has been performed in order to evaluate systematic errors of atmospheric Total Water Vapor (TWV) and integral Liquid Water Content (LWC) microwave radiometric retrieval from satellites by means of dual-frequency method (inverse problem).  The errors under consideration may arise due to the non-linearity of brightness temperature level on true liquid water and effective cloud temperature dependencies and due to neglecting the spatial distribution of cumulus clouds in the satellite microwave radiometer antenna field-of-view (FOV). Comprehensive evaluation of these errors requires sequential forward and then inverse problem solving in different variations for numerous of possible cloudy atmosphere states (3D), satellite locations and observation angles, various FOV sizes, used frequency channel pairs, etc. It is clear that bruteforce approach would be purblind and short-sighted here. On the other hand, it turns out that cumuli are not fully randomly distributed in the atmosphere, but their number density and power distribution is subject to quite understandable statistics (see, e.g. Plank, 1969).In this dataset, we constrain ourselfs to considering standard atmosphere 3D-cells (i.e. altitude profiles of thermodynamic temperature, pressure and air humidity are standard), but additionally filled with cumuli distributed according to several represantative cases taken from (Plank, 1969), see table 3. A number of modifications for these cases are also considered. Herewith, the shape of clouds is assumed to be cylindrical. The size of the cell is set to 50x50x10 km (Ox x Oy x Oz). An introduced computational grid has 300x300x100 nodes.  A smooth water surface with temperature of 15 deg, Cels and zero salinity is considered as an underlying one. An observation angle is taken to be 0 deg. from zenith. Satellite FOV sizes are 1x1 and 60x60 nodes. The altitude profile of liquid water inside a cloud is approximated according to (Mazin, 1983). The effective cloud temperature parameter used in the course of TWV and LWC retrieval is chosen constant and equal to 0 deg. Cels. All the other parameters necessary while retrieval process are either known in advance or determined from the results of solving the forward problem in case of clear sky. The retrieval process is carried out for three frequency pairs: 22.2 and 27.2 GHz, 22.2 and 36 GHz, 22.2 and 89 GHz.The dataset is generated by this script and using armrad framework (Python + TensorFlow GPU).