A timing and precise diagnosis of crop nutrient status is essential for optimizing management practices that promote environmentally friendly and enhanced crop yields. Although plant tissue analysis has conventionally been employed to evaluate the nutritional status of crops, this method cannot capture the spatial variability of crop nutrients. In contrast, satellite-based remote sensing can monitor the nutrient status of crops across expansive areas. This study explored the capability of multi-source satellite images (PlanetScope-4: 3 m, 4 bands; PlanetScope-8: 3 m, 8 bands; Sentinel-2: 10–60 m, 13 bands; PRISMA: 30 m, 239 bands) in mapping 12 foliar nutrients in cranberries. Three machine learning approaches, including partial least squared regression (PLSR), support vector regression (SVR), random forest regression (RFR), were used to relate foliar nutrients to different types of satellite-derived features (SR: surface reflectance; VI: vegetation indices; TF: texture features) or their combinations (SR+VI, VI+TF and SR+VI+TF). Model performance was compared across different foliar nutrients, modelling approaches and combinations of model input features using R2 (the coefficient of determination) and RRMSE (relative root mean square error, = root mean square error/nutrient range × 100 %). Input features that were important to foliar nutrient modelling were identified. The model performance difference among nutrients was consistent between Planet-4 and Sentinel-2, as well as between Planet-8 and PRISMA. In the Planet-4 and Sentinel-2 derived models, N was best predicted (average R2 = 0.77, average RRMSE=15 %), followed by macronutrients S (0.60–0.63, 11 %), Mg (0.58–0.65, 10–11 %), Ca (0.49–0.51, 9 %), Na (0.69, 22 %), P (0.49, 9 %) and K (0.20, 8 %), and then by all micronutrients(i.e., Fe, Mn, B, Cu and Zn: R2 = 0.04–0.61; RRMSE=16–28 %). In the Planet-8 and PRISMA derived models, macronutrients (i.e., N, P, K, Mg, Ca, S and Na) had lower R2 and RRMSE (R2 = 0.06–0.59; RRMSE=7–57 %) than micronutrients (i.e., Fe, Mn, B, Cu and Zn: R2 = 0.18–0.60; RRMSE=19–66 %). The successful retrieval of foliar nutrients from satellite imagery was influenced by many factors, including the intercorrelation between nutrients and model input features, the data availability at critical growth stages, and satellite images characteristics (e.g., spatial and spectral resolutions). Except for foliar nitrogen, foliar nutrients typically do not exhibit distinct absorption features associated with C, H, N, or O molecular bonds in the 400–2500 nm range. Our results indicate that their successful retrieval can be primarily attributed to the association between foliar nutrients and other leaf components (e.g., pigments, water, and dry matter) that do display spectral features within this range. Our study demonstrated the potential of integrating multi-source satellite data for precise nutrient monitoring over large scales.