An upgrade of the software developed at the LMU of Munich and at the ICMMG of Novosibirsk allows for the efficient simulation of CCD lidar returns from highly structured clouds of nonspherical scatterers. The multiple scattering contributions to the lidar return contain additional valuable information about the scattering particles. To make use of it, multi-channel lidars with various fields of view and polarization and CCD lidars have been constructed. The progress in opto-electronics made possible the construction of lidars with charged coupled devices (CCD) which allow for taking time resolved two-dimensional pictures of the diffusion of the lidar beam in a cloud, promising additional gain of information from such CCD lidar returns. Shortly, we discuss the single scattering of polarized light by various types of particles and simulate lidar returns and CCD lidar returns from various types of broken clouds and plumes, compare the two kinds of lidar returns and show the superior capability of CCD lidar returns for the distinction of various types of aerosol clouds and plumes (chimney roses). Here we show simulations of monostatic CCD lidar returns for the distinction of types of clouds and bistatic CCD lidar returns to demonstrate the effect of pulse stretching in different clouds. Our theory is based on the stochastic model of corpuscular multiple scattering of polarized light which is a partially deterministic Markovian jump process and which is equivalent to an appropriate version of the radiative transfer equation for polarized light. For the Monte Carlo simulation we make use of the comfortable and highly effective method of iterative fictitious collisions.