Abstract
This paper proposes a hybrid, interpretable machine learning framework for the preliminary screening of soil macronutrients using Sentinel-2 and AgroLens data. This study aims not to replace laboratory analysis, but to test the feasibility of obtaining a useful proxy signal for estimating nitrogen (N), phosphorus (P), and potassium (K) content using a limited set of remote sensing and agricultural features. The developed pipeline includes data auditing, leakage control, feature engineering, train-only normalization, group-aware partitioning, baseline/SOTA model comparison, hybrid regression modeling, SHAP interpretation, and uncertainty assessment. The experiment used 4471 AgroLens observations and 126 features derived from Sentinel-2 spectral aggregates, vegetation indices, temporal characteristics, and crop-related parameters. The evaluation indicated that the proposed approach consistently improves forecasting quality relative to baseline models under reduced-input conditions. Linear relationships between target variables ranged from 0.14 to 0.17, while nonlinear relationships reached 0.23. SHAP analysis revealed significant contributions from vegetation indices, crop-specific interactions, and Sentinel-2 spectral channels. The findings support the applicability of the proposed framework for preliminary monitoring, prioritizing field surveys, and decision support in digital agriculture. Although an additional AgroLens control segment was used to assess the robustness of the study, the study did not include independent external validation of the data collected across different geographic or agro-climatic conditions.