Using Tropical Rainfall Measuring Mission (TRMM) observations from storms collected over the oceans surrounding East Asia, during summer, a method of creating physically consistent cloud-radiation databases to support satellite radiometer retrievals is introduced. In this method, vertical profiles of numerical model-simulated cloud and precipitation fields are optimized against TRMM radar and radiometer observations using a hybrid empirical orthogonal function (EOF)–one-dimensional variational (1DVAR) approach.The optimization is based on comparing simulated to observed radar reflectivity profiles and the corresponding passive microwave observations at the frequencies of the TRMM Microwave Imager (TMI) instrument. To minimize the discrepancies between the actual and the synthetic observations, the simulated cloud and precipitation profiles are optimized by adjusting the contents of the hydrometeors. To reduce the dimension of the hydrometeor content profiles in the optimization, multivariate relations among hydrometeor species are used. After applying the optimization method to modify the simulated clouds, the optimized cloud-radiation database has a joint distribution of reflectivity and associated brightness temperatures that is considerably closer to that observed by TRMM PR and TMI, especially at 85 GHz. This implies that the EOF–1DVAR approach can generate profiles with realistic distributions of frozen hydrometeors, such as snow and graupel. This approach may be similarly adapted to operate with the variety and capabilities of the passive microwave radiometers that compose the Global Precipitation Measurement (GPM) constellation. Furthermore, it can be extended to other oceanic regions and seasons.