Abstract. Understanding the uncertainties in the retrieval of aerosol and surface properties is very important for an adequate characterization of the processes that occur in the atmosphere. However, the reliable characterization of the error budget of the retrieval products is a very challenging aspect that currently remains not fully resolved in most remote sensing approaches. The level of uncertainties for the majority of the remote sensing products relies mostly on post-processing validations and intercomparisons with other data, while the dynamic errors are rarely provided. Therefore, implementations of fundamental approaches for generating dynamic retrieval errors and the evaluation of their practical efficiency remains of high importance. This study describes and analyses the dynamic estimates of uncertainties in aerosol-retrieved properties by the GRASP (Generalized Retrieval of Atmosphere and Surface Properties) algorithm. The GRASP inversion algorithm, described by Dubovik et al. (2011, 2014, 2021), is designed based on the concept of statistical optimization and provides dynamic error estimates for all retrieved aerosol and surface properties. The approach takes into account the effect of both random and systematic uncertainties propagations. The algorithm provides error estimates both for directly retrieved parameters included in the retrieval state vector and for the characteristics derived from these parameters. For example, in the case of the aerosol properties, GRASP directly retrieves the size distribution and the refractive index that are used afterwards to provide phase function, scattering, extinction, single scattering albedo, etc. Moreover, the GRASP algorithm provides full covariance matrices, i.e. not only variances of the retrieval errors but also correlations coefficients of these errors. The analysis of the correlation matrix structure can be very useful for identifying less than obvious retrieval tendencies. This appears to be a useful approach for optimizing observation schemes and retrieval set-ups. In this study, we analyse the efficiency of the GRASP error estimation approach for applications to ground-based observations by a sun/sky photometer and lidar. Specifically, diverse aspects of the error generations and their evaluations are discussed and illustrated. The studies rely on a series of comprehensive sensitivity tests when simulated sun/sky photometer measurements and lidar data are perturbed by random and systematic errors and inverted. Then, the results of the retrievals and their error estimations are analysed and evaluated. The tests are conducted for different observations of diverse aerosol types, including biomass burning, urban, dust and their mixtures. The study considers observations of AErosol RObotic NETwork (AERONET) sun/sky photometer measurements at 440, 675, 870 and 1020 nm and multiwavelength elastic lidar measurements at 355, 532 and 1064 nm. The sun/sky photometer data are inverted alone or together with lidar data. The analysis shows overall successful retrievals and error estimations for different aerosol characteristics, including aerosol size distribution, complex refractive index, single scattering albedo, lidar ratios, aerosol vertical profiles, etc. Also, the main observed tendencies in the error dynamic agree with known retrieval experience. For example, the main accuracy limitations for retrievals of all aerosol types relate to the situations with low optical depth. Also, in situations with multicomponent aerosol mixtures, the reliable characterization of each component is possible only in limited situations, for example, from radiometric data obtained for low solar zenith angle observations or from a combination of radiometric and lidar data. At the same time, the total optical properties of aerosol mixtures are always retrieved satisfactorily. In addition, the study includes an analysis of the detailed structure of the correlation matrices for the retrieval errors in mono- and multicomponent aerosols. The conducted analysis of error correlation appears to be a useful approach for optimizing observation schemes and retrieval set-ups. The application of the approach to real data is provided.