ABSTRACT The temporal stability and spatial dependence of pixel-specific red-NIR soil line coefficients were studied using analyses of geostatistics and covariance. We used time series of MODIS 8-day composite reflectance data from 2000 to 2019 of bare soil surfaces near Kabul, Afghanistan. The goal was to explore the feasibility of reducing the soil background signal in red-NIR vegetation indices (VIs). The results were that 1) Red-NIR pixel-specific soil line coefficients can be obtained using multi-temporal remote sensing imagery; 2) pixel-specific Red-NIR soil line coefficients derived from multi-temporal remote sensing imagery exhibit clear spatial dependency related to the soil line coefficients of soils sampled in the field, and 3) Red-NIR pixel-specific soil line coefficients were stable for years, but showed instability for periods longer than 16 years. These results imply that A) some changes in soil surface characteristics can be deduced using long-term observations, B) pixel-specific Red-NIR soil lines at any location can be estimated by kriging interpolation. Consequently, soil line coefficients derived from multi-temporal satellite observations can be used to derive information about vegetation, such as the germination stages pixel-by-pixel or above-ground biomass in sparsely vegetated lands. The derived values should be valid over wide areas, minimizing soil background effects.