SOCFARMSOCFARM
AL-6222052520 · 2023–2026

Supporting organic and conventional farming through geospatial methods

SOCFARM is an international applied project of Uzbek–Hungarian cooperation. It combines Sentinel‑2, GIS, field spectrometry, soil and yield measurements to develop a methodology for assessing winter wheat condition, improving fertilizer and water management, and predicting yield.

NDVI Sentinel‑2 vegetation status and fieldwork
1791 m² yield samples
320+field spectra
112soil samples
R² 0.71ExtraTrees model
Project codeAL-6222052520
Duration1 October 2023 — 30 September 2026
Uzbek partnerTIIAME National Research University
Hungarian partnerUniversity of Nyíregyháza
Study areasMiddle Chirchiq and Boka districts, Tashkent region
Main cropWinter wheat
Key technologiesSentinel‑2, GEE, ArcGIS Pro/QGIS, GNSS, ATP9110, Python/ML
Application areasAgriculture and water management, agrochemical services, higher education

The practical problem addressed

Three information gaps in field management were identified:

  1. Crop condition was assessed mostly by visual observation and experience; satellite monitoring was rarely used.
  2. Soil fertility data were outdated or insufficient; fertilizers were often applied at a uniform rate across the field.
  3. A locally validated methodology was lacking to justify irrigation timing, yield forecasting and management zones.

Farmers asked for simple, mobile-friendly maps and decision support. These needs shaped the project workflow: diagnostics → monitoring → prediction → web dissemination.

Needs assessment with farmers, December 2023
Needs assessment with farmers, December 2023
Goal

To support organic and conventional farming with geospatial methods and develop a methodology for assessing crop condition, fertilizer and water demand, and yield from multi-source data.

Key objectives

1

Needs assessment and selection of trial fields

2

Field-journal monitoring of phenology, irrigation and fertilization

3

Measuring and mapping N, P, K, humus and moisture

4

Building a cereal spectral library with ATP9110

5

Computing Sentinel‑2 time series and vegetation indices

6

Collecting 1 m² yield samples and evaluating the ML model

7

Disseminating results via website, trainings and scientific publications

Complementary roles of the partners

TIIAME NRU

TIIAME National Research University

  • GIS and remote sensing
  • Field, soil and yield measurements
  • Sentinel‑2 and machine learning
  • Conventional irrigated farming
University of Nyíregyháza

University of Nyíregyháza

  • Organic farming expertise
  • Organic fertilizer and triticale trials
  • Training and international conferences
  • Joint publications and exchange of experience

Study areas

Middle Chirchiq — university teaching and research farm
Middle Chirchiq — university teaching and research farm
Boka — “Hasan Hamza” and neighbouring farms
Boka — “Hasan Hamza” and neighbouring farms

2025: sample segments of 10.07 ha / 13.8 ha · 2026: 39.8 ha / 8 fields, ≈69 ha

Scientific workflow

Field → laboratory/processing → analysis → result
Field → laboratory/processing → analysis → result
  • All measurements at georeferenced points
  • Field–lab–Sentinel‑2 linkage through a common ID
  • Quality control: duplicates, blanks, calibration, QC
  • Model evaluation with grouped cross-validation
  • Turning results into maps and web products
Timeline 2023–2026: year 1 — design and baseline data → year 2 — denser sampling and spectra → year 3 — soil profiles, neighbouring fields and final integration
Timeline 2023–2026: year 1 — design and baseline data → year 2 — denser sampling and spectra → year 3 — soil profiles, neighbouring fields and final integration
Stakeholder needs assessment
Stakeholder needs assessment
Stakeholder needs assessment
Stakeholder needs assessment
01 · 2023

Stakeholder needs assessment

In December 2023, surveys and face-to-face interviews covered 5 farms and 28 questions: water and soil management, use of satellite monitoring, yield forecasting and economic potential, and interest in organic methods.

Phenological and agrotechnical observations
Phenological and agrotechnical observations
Phenological and agrotechnical observations
Phenological and agrotechnical observations
02 · Monitoring

Phenological and agrotechnical observations

The field journal consistently recorded irrigation, fertilization, crop protection and growth stages from sowing to harvest, season by season.

Soil sampling: 2024–2025
Soil sampling: 2024–2025
Soil sampling: 2024–2025
Soil sampling: 2024–2025
03 · 2024–2025

Soil sampling: 2024–2025

112 samples · WGS84 coordinates · duplicate and blank control · sample–lab–GIS consistency.

2026: soil profiles and area expansion
2026: soil profiles and area expansion
2026: soil profiles and area expansion
2026: soil profiles and area expansion
04 · 2026

2026: soil profiles and area expansion

22 profiles × 2 layers. In Boka, 8 wheat fields of 6 farms were covered; in Middle Chirchiq, 39.8 ha of the university farm.

Laboratory preparation and agrochemical analysis
Laboratory preparation and agrochemical analysis
05 · Laboratory

Laboratory preparation and agrochemical analysis

Sample reception, drying, sieving and analysis. Indicators: humus (Tyurin method), N–NO₃ and P₂O₅ (spectrophotometry), K₂O (flame photometry), oven 105 °C, mechanical composition.

QC: blank, control standard, ≥10% duplicates, calibration R² ≥ 0.995

Field spectroradiometry
Field spectroradiometry
Field spectroradiometry
Field spectroradiometry
06 · Spectrometry

Field spectroradiometry

Spectroradiometer data were collected in the field with the OptoSky ATP9110, using a white reference panel and repeated measurements.

Ground truth for yield
Ground truth for yield
Ground truth for yield
Ground truth for yield
07 · Yield

Ground truth for yield

Yield was harvested within 1 × 1 m quadrats: ear count, grain count and grain mass were measured at each point and linked via GNSS coordinates to the Sentinel‑2 pixel and the yield model.

112topsoil samples
22soil profiles
180+gravimetric moisture readings
320+field spectra
1791 m² yield samples

All data types were collected on the same fields at georeferenced points. This made it possible to link the soil–plant–satellite–yield chain and independently validate the model.

P

Available phosphorus is the main limiting factor

S2

Dynamics of NDVI, EVI, NDRE, GNDVI, SAVI, NDMI

SPEC

Hyperspectral field library for cereals

ML

ExtraTrees: R² 0.71; RMSE 10.8 c/ha

IP

2 intellectual property objects

Soil nutrient supply classes, 2025

Available phosphorus fell into the “very low” or “low” class at many points and was identified as the main limiting factor. Potassium was generally high to very high, humus mostly low to medium, and nitrogen varied with irrigation and agronomic background.

Share of classes for N, P, K and humus
Share of classes for N, P, K and humus

2025 — spatial distribution of nitrogen and phosphorus

Middle Chirchiq — N
Middle Chirchiq — N
Boka — N
Boka — N
Middle Chirchiq — P
Middle Chirchiq — P
Boka — P
Boka — P

Spatial variability of nitrogen depends on irrigation and fertilization history; phosphorus deficiency in both areas calls for zonal, differentiated fertilization.

2025 — potassium and humus

Middle Chirchiq — K
Middle Chirchiq — K
Boka — K
Boka — K
Middle Chirchiq — humus
Middle Chirchiq — humus
Boka — humus
Boka — humus

Potassium reserves are sufficient in many zones; to preserve humus, organic matter, residue return and minimum tillage are recommended.

2026 — Middle Chirchiq soil maps (nitrogen, phosphorus, potassium, humus)

Nitrogen
Nitrogen
Phosphorus
Phosphorus
Potassium
Potassium
Humus
Humus

2026 — Boka area soil maps (nitrogen, phosphorus, potassium, humus)

Five neighbouring farms were covered together with the main field; results enabled management zones to be delineated at a wider scale.

Nitrogen
Nitrogen
Phosphorus
Phosphorus
Potassium
Potassium
Humus
Humus

Soil properties, yield and three-year comparison

Relationship between soil indicators and yield
Relationship between soil indicators and yield
Three-year and cross-area comparison
Three-year and cross-area comparison

Results show that a universal recommendation based on a single parameter is insufficient: fertilizer, moisture, texture, irrigation and phenology must be interpreted together.

Important negative result. Soil moisture under dense canopy could not be reliably estimated from Sentinel‑2 optical indices; in 2026 all cross-validated models gave R²CV < 0. For irrigation decisions, Sentinel‑1 radar, permanent sensors and synchronous field measurements are recommended.
Middle Chirchiq — moisture / NDMI
Middle Chirchiq — moisture / NDMI
Boka — moisture / NDMI
Boka — moisture / NDMI
2025 regression
2025 regression
2026 cross-validation
2026 cross-validation

Time series of NDVI, EVI, GNDVI, NDRE, SAVI and NDMI supported spatio-temporal monitoring of growth, heading, grain filling and ripening.

Boka: indices and phenological stages
Boka: indices and phenological stages
Middle Chirchiq: indices and phenological stages
Middle Chirchiq: indices and phenological stages

Spatial variation of NDVI

2024 season
2024 season
2025, Middle Chirchiq
2025, Middle Chirchiq
2025, Middle Chirchiq
2025, Middle Chirchiq

NDVI differences over time and space reflected variations in irrigation technology, cultivar, agronomic background and canopy density.

Field vs Sentinel‑2 NDVI: scale mismatch

Leaf-level spectroradiometer NDVI vs 10 m Sentinel‑2 pixel
Leaf-level spectroradiometer NDVI vs 10 m Sentinel‑2 pixel
R² ≤ 0.03three campaigns
1.8–3.2×field NDVI higher
  • Leaf level: centimetre scale
  • Sentinel‑2: 100 m² mixed pixel
  • Soil background, shadow, ditches and weeds
  • Date and GNSS errors as additional factors

Field spectra and satellite values cannot be equated directly: band equivalents and a scale-appropriate protocol are required.

2026 cereal spectral library: 23 narrow-band indices — vegetation, chlorophyll, red-edge, pigments, photosynthetic efficiency and water status.

Smoothed spectra — area 1
Smoothed spectra — area 1
Smoothed spectra — area 2
Smoothed spectra — area 2
Mean spectrum and variability — area 1
Mean spectrum and variability — area 1
Mean spectrum and variability — area 2
Mean spectrum and variability — area 2

Red-edge, PCA and Sentinel‑2 band equivalents

Red-edge spectra
Red-edge spectra
First derivative and O₂‑A artefact
First derivative and O₂‑A artefact
PCA loadings — Boka
PCA loadings — Boka
PCA loadings — Middle Chirchiq
PCA loadings — Middle Chirchiq
Boka — Sentinel‑2 band equivalents
Boka — Sentinel‑2 band equivalents
Middle Chirchiq — Sentinel‑2 band equivalents
Middle Chirchiq — Sentinel‑2 band equivalents

PC1 + PC2 explained 88.6–94.1 % of total variance; the O₂‑A artefact near 760 nm was separately controlled in red-edge calculations.

179samples
16–106range, c/ha
R² 0.71grouped CV
10.8RMSE, c/ha
7.7MAE, c/ha
≈19%rRMSE

Model: seven mid-May Sentinel‑2 bands + five indices → ExtraTrees regression → cross-validation grouped by field and season.

Yield distribution across three seasons and irrigation variants
Yield distribution across three seasons and irrigation variants

Most important spectral features

ExtraTrees permutation importance (ΔR²)
ExtraTrees permutation importance (ΔR²)

SWIR1 reflected canopy water status; green and NDRE reflected chlorophyll/nitrogen status; NIR reflected biomass and leaf area. Highest importance: SWIR1 (0.40), Green (0.34), NIR (0.28), NDRE (0.21), SWIR2 (0.13).

Spatial wheat yield prediction maps

Middle Chirchiq — variant 1
Middle Chirchiq — variant 1
Middle Chirchiq — variant 2
Middle Chirchiq — variant 2
Boka — 2025
Boka — 2025
Middle Chirchiq — 2025
Middle Chirchiq — 2025
Methodological caveat. The predicted range is narrower than the measured yield range. The maps are therefore presented as relative productivity and management zones, not absolute yield values. To become a practical forecasting tool, the model must be tested on independent 2026 samples and retrained with three years of data.
Distribution of map values
Distribution of map values
In-field distribution
In-field distribution
Boka distribution
Boka distribution
Middle Chirchiq distribution
Middle Chirchiq distribution

Practical recommendations and implementation directions

1

Phosphorus

Soil-test-based, contour-wise differentiated P fertilization in phosphorus-deficient zones.

2

Water

Do not use optical indices in place of soil moisture; combine Sentinel‑1, sensors and the irrigation log.

3

Monitoring

Regularly track Sentinel‑2 time series in the heading–grain filling window.

4

Management zones

Use yield and soil maps for control sampling and variable-rate application.

5

Model validation

Update the model with independent 2026 data, spatial CV and calibrated absolute yield.

6

Training and services

Integrate the spectral library, protocols and GIS products into teaching and advisory services.

Final report

PDF/DOCX, 2023–2026

Monograph

Cover, ISBN and PDF

Field protocols

Soil, spectra, yield, QC

Maps

PNG/PDF/GeoTIFF and legends

Datasets

CSV/XLSX, metadata and licence

Newsletters

By year and language

Intellectual property

Certificate and short description

Presentations

Seminars, conferences, trainings

Monograph cover
Monograph cover
Monograph

New experience of geospatial and Geo‑AI technologies in agriculture

Sh.S. Shokirov, M.P. Reymov, I.I. Abdurahmonov, O.E. Xaqberdiyev, Z.J. Mamatkulov, X.X. Hasanov

Tashkent: TIIAME NRU, 2025 — 86 pp. · ISBN 978-9910-8486-7-4

  • Field and data collection methodology
  • Soil analysis and spatial mapping
  • Sentinel‑2 and soil moisture
  • Spectral library
  • Yield prediction with machine learning

Datasets, maps and protocols

Spectra files and metadata — CSV/XLSX: sample ID, coordinates, layer, units and QC
Spectra files and metadata — CSV/XLSX: sample ID, coordinates, layer, units and QC
Sentinel‑2 indices and time series
Sentinel‑2 indices and time series
Soil nutrient maps — GeoTIFF/PNG/PDF: projection, legend, date
Soil nutrient maps — GeoTIFF/PNG/PDF: projection, legend, date

Protocols: sampling, laboratory, spectra, 1 m² yield. A data dictionary and version log are provided; licence and personal/restricted data are separated.

Intellectual property objects

“GeoCalc” — certificate No. 010103, 03.04.2026
“GeoCalc” — certificate No. 010103, 03.04.2026
“Geodata Visualizator” — DGU 64712, 22.05.2026
“Geodata Visualizator” — DGU 64712, 22.05.2026

Textbooks and study guides

  1. Geoinformation cartography
  2. Modelling reservoirs with geoinformation technologies
  3. Geodetic metrology
  4. Cadastre of protected natural areas
  5. Geoecological mapping
  6. Space geodesy
  7. Scientific Foundations of Geoinformation Systems

Scientific capacity and training

On the project topic and related fields, 1 DSc (S.N. Abdurakhmonov, 30.04.2026), 1 PhD (R.I. Makhsudov) and 2 master’s theses were defended.

1Monograph
7Textbooks / guides
14Scopus articles
3Foreign journals
10International conferences
31National journals
10National conferences
2Intellectual property

Annual dynamics: 12 → 20 → 44 scientific outputs (2023 → 2025).

One of the project publications
One of the project publications

DOI-verified publications

  1. Leveraging Vegetation Indices and Random Forest for Soil Nutrient Monitoring in Winter Wheat

    Journal of the Indian Society of Remote Sensing, 2025

    DOI: 10.1007/s12524-025-02198-9 ↗
  2. The effect of fermented poultry manure on the yield and quality parameters of triticale cultivated organic way on acidic sandy soil

    E3S Web of Conferences 590, 01009, 2024

    DOI: 10.1051/e3sconf/202459001009 ↗
  3. Temporal Evolution of Climate-Related Land Changes in Zarafshan Valley: MODIS and Google Earth Engine Perspective

    AIP Conference Proceedings 3256, 040016, 2025

    DOI: 10.1063/5.0266743 ↗
  4. Search for the lost Baghi Chinar garden in the Samarkand region and creation of a Digital Terrain Model

    E3S Web of Conferences 590, 05002, 2024

    DOI: 10.1051/e3sconf/202459005002 ↗
  5. Sustainable development of Turkmenistan in the context of climate change

    E3S Web of Conferences 590, 05007, 2024

    DOI: 10.1051/e3sconf/202459005007 ↗
  6. Land use and land management in Kazakhstan

    Problems of AgriMarket, No. 3, 2024

    DOI: 10.46666/2024-3.2708-9991.14 ↗
  7. Theoretical and Methodological Basis of Mapping the Ecological State of Irrigated Lands in the Conditions of Climate Change: Bibliometric Analysis

    ISPRS Archives, 2025

    DOI: 10.5194/isprs-archives-XLVIII-5-W3-2025-63-2025 ↗

2026 publications

  • Multi-temporal monitoring of vegetation dynamics using Sentinel‑2 and Google Earth Engine — Proceedings of SPIE, 2026.
  • Preparation and integration of spatial data in a GIS environment: the case of Boka and Middle Chirchiq — 2026.
  • Spatial variability of soil agrochemical indicators and Trimble Juno 3B / ArcGIS Survey123 — Nukus, 2026.
  • Soil fertility and fertilization efficiency using remote sensing and field analyses — 2026.
  • Yield components of triticale with granulated poultry manure — Nyíregyháza international conference, 2026.
  • Remote Sensing for Precision Agriculture in Arid and Semi-Arid Regions: A Case Study of Uzbekistan — Springer Nature book chapter, in preparation.
Visit · 2–6 June 2025

Scientific and practical visit to Hungary

Nyíregyháza, Budapest and Debrecen — conference, trainings, meetings with universities and agencies.

Nyíregyháza: conference and training
Nyíregyháza: conference and training
Óbuda University: Institute of Geoinformatics
Óbuda University: Institute of Geoinformatics
University of Debrecen: scientific meeting
University of Debrecen: scientific meeting

The “Sustainable Geospatial Solutions for a Changing World” conference and trainings at the University of Nyíregyháza.

Meeting with the Drought Prevention Institute and HEPA in Budapest. At Óbuda University — talks on remote sensing, drought mapping and a master’s programme.

At the University of Debrecen — cooperation in agriculture, water management, geoinformatics and IT.

Results of international cooperation

3

Scopus publications with Hungarian partners

2026

Joint triticale results at the Nyíregyháza conference

1

Book chapter for Springer Nature (collaboration with Tsukuba University)

4

Meetings with universities and agencies in Hungary

Partner scientists: Judit Csabai, Zoltán Kovács, Edit Kosztyuné Krajnyák, Béla Szabó and others.

Discussion with farmers
Discussion with farmers
Soil sampling
Soil sampling
Field spectrometry
Field spectrometry
Debrecen — scientific meeting
Debrecen — scientific meeting
Winter wheat before harvest
Winter wheat before harvest
International cooperation
International cooperation
Laboratory work
Laboratory work
Working with the spectroradiometer in the field
Working with the spectroradiometer in the field
Conference and training
Conference and training
1 m² yield sample
1 m² yield sample
Óbuda University
Óbuda University
Fieldwork
Fieldwork
Needs assessment
Needs assessment
Nyíregyháza — conference
Nyíregyháza — conference
Wheat field
Wheat field
Field experiments
Field experiments
Soil analysis
Soil analysis

Leadership and leading researchers

MR
Mamanbek ReymovProject leader, leading researcherPhD, Associate Professor
SS
Shuhrat ShokirovChief researcherDSc, Professor
SI
Shavkat ImomovLeading researcherDr.Sc. (Tech.), Professor
OX
Obid XaqberdiyevSenior researcherCand.Sc. (Biol.), Assoc. Prof.

Senior researchers

SA
Sarvar AbduraxmonovSenior researcherDSc, Associate Professor
RO
Rustam OymatovSenior researcherPhD, Associate Professor
OM
O‘zbekxon MuxtorovSenior researcherPhD, Associate Professor
IA
Ilhom AbdurahmanovSenior researcherPhD
ZM
Zoxid MamatqulovSenior researcherPhD
BS
Barno SoipovaSenior researcherPhD
AA
A.F. AshurovSenior researcherPhD
SN
Sh.K. NarbayevSenior researcherPhD
AB
A.R. BabajanovSenior researcherDr.Sc. (Econ.), Assoc. Prof.
MR
M.X. RajapbayevSenior researcherPhD
UQ
U.T. Qo‘ziyevSenior researcherPhD
IM
I.M. MusayevSenior researcherCand.Sc. (Tech.), Assoc. Prof.
OI
O‘.P. IslomovSenior researcherAssociate Professor
AI
A.N. InamovSenior researcherPhD, Associate Professor
AJ
A.N. JumanovSenior researcherPhD, Associate Professor

Junior researchers

ZA
Zoir AbdurahmonovJunior researcher
NT
N.N. TeshayevJunior researcher
XX
X.X. XasanovJunior researcher
SH
Sh.Sh. HamidullayevJunior researcher
SI
Sh.R. IsmonovJunior researcher
MS
M.S. SobirovJunior researcher
AO
A.R. OrazbayevJunior researcher
NE
N.A. EshmatovaJunior researcher
TR
T.P. ReymovJunior researcher
DZ
Doniyor ZarifboyevGIS engineer
BO
B.O. OralbayevJunior researcher
MA
M.S. AkbarovJunior researcher, senior lecturer
ZX
Z.X. XafizovaJunior researcher, senior lecturer
YU
Yu.A. UsmanovJunior researcher
IA
I.M. AslanovJunior researcher

Engineers

BB
B.N. BektashovEngineer
SJ
S.X. JamardovEngineer
AD
A.D. DavlatovEngineer
MF
Mirzohid FayzullayevEngineer, co-executor “Tashgistechnologies”

International partner scientists

JC
Judit CsabaiNyíregyháza University — partner scientist
ZK
Zoltán KovácsNyíregyháza University — partner scientist
EK
Edit Kosztyuné KrajnyákNyíregyháza University — partner scientist
BS
Béla SzabóNyíregyháza University — partner scientist
Project leader

Mamanbek Reymov

PhD, Associate Professor
TIIAME National Research University