Lidar depolarization measurements distinguish between spherical and non-spherical aerosol particles based on the change of the polarization state between the emitted and received signal. The particle shape information in combination with other aerosol optical properties allows the characterization of different aerosol types and the retrieval of aerosol particle microphysical properties. Regarding the microphysical inversions, the lidar depolarization technique is becoming a key method since particle shape information can be used by algorithms based on spheres and spheroids, optimizing the retrieval procedure. Thus, the identification of the depolarization error sources and the quantification of their effects are crucial. This work presents a new tool to assess the systematic error of the volume linear depolarization ratio (delta), combining the Stokes-Muller formalism and the complete sampling of the error space using the lidar model presented in Freudenthaler (2016a). This tool is applied to a synthetic lidar system and to several EARLINET lidars with depolarization capabilities at 355 or 532 nm. The lidar systems show relative errors of delta larger than 100% for delta values around molecular linear depolarization ratios (similar to 0.004 and up to similar to 10% for delta = 0.45). However, one system shows only relative errors of 25 and 0.22% for delta = 0.004 and delta = 0.45, respectively, and gives an example of how a proper identification and reduction of the main error sources can drastically reduce the systematic errors of delta. In this regard, we provide some indications of how to reduce the systematic errors.