The availability of reliable information on the physicochemical properties of soils is a necessary tool for maintaining and improving fertility and effective optimization of agricultural land management in many countries. However, ground-based research methods require significant financial and time resources. It has been established that methods based on remote sensing data are an efficient, accurate, and less costly solution for studying different types of soil cover parameters. This work aims to determine the predicted indicator of humus content in soils in selected regions of the Kyiv region (Ukraine) with the corresponding soil types. For this, the spectral properties of chernozem soils were investigated based on Landsat 8 OLI satellite images. A mosaic of the mean spectral reflectance values for the study period (2013-2015) was created using the Google Earth Engine. The vegetation indices NDSI, NDWI, NDBI, MSAVI, and NDVI were used to identify bare soils. Using multiple linear regression, an optimal F-Comparing Nested Model was created for predicting humus content in soils, including seven parameters. The model's accuracy was estimated with such indicators R=0.95, R2= 0.90, σy = 0.16 %. The approach based on the proposed model can be used to support the adoption of the necessary management decisions to improve soil fertility and maintain balanced land use.