The long wavelength part of the Earth’s gravity field can be determined, with varying accuracy, from satellite laser ranging (SLR). In this study, we investigate the combination of up to ten geodetic SLR satellites using iterative variance component estimation. SLR observations to different satellites are combined in order to identify the impact of each satellite on the estimated Stokes coefficients. The combination of satellite-specific weekly or monthly arcs allows to reduce parameter correlations of the single-satellite solutions and leads to alternative estimates of the second-degree Stokes coefficients. This alternative time series might be helpful for assessing the uncertainty in the impact of the low-degree Stokes coefficients on geophysical investigations. In order to validate the obtained time series of second-degree Stokes coefficients, a comparison with the SLR RL05 time series of the Center of Space Research (CSR) is done. This investigation shows that all time series are comparable to the CSR time series. The precision of the weekly/monthly C₂₁ and S₂₁ coefficients is analyzed by comparing mass-related equatorial excitation functions χ^mass_1,2 with geophysical model results and reduced geodetic excitation functions. In case of χ^mass₁, the annual amplitude and phase of the DGFI solution agrees better with three of four geophysical model combinations than other time series. In case of χ^mass₂, all time series agree very well to each other. The impact of C₂₀ on the ice mass trend estimates for Antarctica are compared based on CSR GRACE RL05 solutions, in which different monthly C₂₀ time series are used for replacing. We found differences in the long-term Antarctic ice loss of 12.3 Gt/year between the GRACE solutions induced by the different C₂₀ SLR time series of CSR and DGFI, which is about 13 % of the total ice loss of Antarctica. This result shows that Antarctic ice mass loss quantifications must be carefully interpreted.