Abstract. The eXtensible Bremen Aerosol/cloud and surfacE parameters Retrieval (XBAER) algorithm has been applied on the Top-Of-Atmosphere reflectance measured by the Sea and Land Surface Temperature Radiometer (SLSTR) instrument onboard Sentinel-3 to derive snow properties: Snow Grain Size (SGS), Snow Particle Shape (SPS) and Specific Surface Area (SSA) under cloud-free conditions. This is the first part of the paper, to describe the retrieval method and the sensitivity study. Nine pre-defined ice crystal particle shapes (aggregate of 8 columns, Drontal, hollow bullet rosettes, hollow column, plate, aggregate of 5 plates, aggregate of 10 plates, solid bullet rosettes, column) are used to describe the snow optical properties. The optimal SGS and SPS are estimated iteratively utilizing a Look-Up-Table (LUT) approach. The SSA is then calculated using another pre-calculated LUT for the retrieved SGS and SPS. The optical properties (e.g., phase function) of the ice crystals can reproduce the wavelength-dependent/angular-dependent snow reflectance features, compared to laboratory measurements. A comprehensive study to understand the impact of aerosol, ice crystal shape, ice crystal surface roughness, and cloud contamination on the retrieval accuracy of snow properties has been performed based on SCIATRAN radiative transfer simulations. The main findings are (1) Snow angular and spectral reflectance feature can be described by the predefined ice crystal properties only when both SGS and SPS can be optimally and iteratively obtained; (2) The impact of ice crystal surface roughness plays minor effects on the retrieval results; (3) SGS and SSA show an inverse linear relationship; (4) The retrieval of SSA assuming non-convex particle shape, compared to convex particle (e.g. sphere) shows larger results; (5) Aerosol/cloud contamination due to unperfected atmospheric correction and cloud screening introduces underestimation of SGS, inaccurate SPS and overestimation of SSA.