Digital monitoring of agro-ecosystems indicators on the basis of space and unmanned technologies

The aim of the research is to develop a digital monitoring methodology for agro-ecosystem indicators based on space and unmanned technologies with the conversion of digital aerial photography results into real parameters of agrophytocenosis indicators as physical units. The paper uses the results of research on digital monitoring of agrophytocenosis of singleand poly-species crops of grain and leguminous crops on the experimental field of the Samara State Agrarian University in 2018 (southern forest-steppe of Trans-Volga region). The soil is a typical medium heavy black loamy chernozem (humus content is 5.9 %, easily hydrolysed nitrogen is 80–120 mg/kg, labile phosphorus is 135–145 mg/kg, exchangeable potassium is 150–195 mg/kg, salt extract pH is 6.9). The research conducted in 2018 revealed that the stated aim is achievable in principle. According to the results of 2018, the greatest correlation between the vegetative biomass index NDVI and grain yield for winter and spring cereal crops was detected at the milky stage of the grains (up to r = 0.82).


Introduction
In 2018, the national project "Digital Economy" was adopted in Russia at the national legislative level. It was one of thirteen Russian national projects for the period from 2019 to 2024 (the passport of the project was approved by the decision of the Presidium of the Presidential Council on Strategic Development and National Projects on December 24, 2018) [1].
The national project "Digital Economy" was created to solve, among others, the following problems: the transformation of priority sectors of the economy and social sphere, including health care, education, industry, agriculture, building construction, urban economy, transport and energy infrastructure, financial services, through the introduction of digital technology and platforms. Thus, the digital transformation of the agroindustrial complex becomes a priority for state authorities, research institutions and agrarian educational institutions [2].
Development of space and unmanned technologies, unmanned aircraft systems (UAS) and special software (including platform solutions) in agriculture and crop production will allow increasing mobility, efficiency and quality of monitoring indicators of agroecosystems and associated natural ecosystems in the framework of adaptive landscape agriculture as a whole [3][4][5][6][7][8].
The following areas are the most promising for the use of unmanned aerial photography: building digital elevation models, analyzing the heterogeneity of soil fertility by spectral characteristics, monitoring work in the fields, assessing crop variability by vegetation indexes, thermal maps of crops, identifying damage to vegetation from external influences [3,4]. Digital monitoring of agroecosystem indicators is carried out on the basis of assessment of crop variability by vegetation indexes, using aerial photographs on a multispectral camera.
Having the data of multispectral photography, we can calculate a variety of vegetation indexes: Normalized Difference Vegetation Index (NDVI), which allows to conduct the qualitative and quantitative assessment of the volume of biomass, to assess the intensity of plant vegetation; Enhanced Normalized Difference Vegetation Index (ENVI), which is similar to NDVI, but also uses part of the visible spectrum for more effective indication of plant health status; Green Normalized Difference Vegetation Index (GNDVI), which allows to assess the content of chlorophyll, the degree of aging, the presence of stress in the plant; Difference Vegetation Index (DVI), which allows you to determine areas of darkening, stagnation of water, soil and vegetation; Green Difference Vegetation Index (GDVI), which shows the amount of nitrogen in the leaves, which can optimize fertilizer application based on actual needs [3,4].
Vegetation indexes give an overall picture of the status of certain indicators of agrophytocenosis, as they are relative values (for example, the values of the NDVI index from 0 to 1). However, for agricultural producers it is important to have specific data in physical terms, for example, the volume of biomass in the field or the potential crop productivity in tons per hectare.
Key words: digital transformation of the agroindustrial complex, space and unmanned technologies, unmanned aircraft systems, digital monitoring of agroecosystem indicators, multispectral imaging, vegetation indexes, yield structure and crop productivity.

Study objectives
Objectives of this work are: 1. Development of a methodology for using software and hardware for space and unmanned technologies and unmanned aircraft systems (UAS) for digital monitoring of agroecosystem indicators. 2. Development of a digital platform for online conversion of vegetation indexes into specific parameters of agrophytocenosis indicators in physical units (for example, t/ha, mg/kg and others), applicable to different crops at different vegetation phases in curtain soil and climatic conditions.

Methods of the research
We performed digital monitoring using the method of obtaining high-resolution images at different crop vegetation stages using unmanned aircraft systems (UAS) using vegetation indexes, including the biomass index (NDVI). Indicators of agrophytocenosis in the form of field statistics were determined according to generally accepted methods based on national standard (GOST). The results of aerial photography were calibrated by the collected data of the field statistics of the relevant parts of experimental fields.

Results
Creating a digital platform to convert the results of various types of digital aerial photography into real physical indicators of crops depends directly on the degree of fullness of statistical databases. That is why the work on creating a platform is intended to continue for 3-4 years, as far as the collection and accumulation of field statistics on various crops is made.
In 2018, the laboratory "AgroEcology" at the Department of Land Management, Soil Science and Agrochemistry of the Samara State Agrarian University together with the company "GIS-R" Ltd. (Samara) conducted research in experimental fields on the following topic "The digital monitoring of agroecosystems indicators based on space and unmanned technologies as the basis of Digital Organic Farming", in accordance with the signed Contract N810/K/2018 dated 03.08.2018 [9].
The experimental plot is located in the central zone of the Samara oblast, in the southern forest-steppe of the Trans-Volga region. The soil of the experimental plot is a typical medium heavy black loamy chernozem (agrochemical parameters of the layer 0-40 cm: humus content is 5.9 %, easily hydrolyzed nitrogen is 80-120 mg/kg, labile phosphorus is 135-145 mg/kg, exchangeable potassium is 150-195 mg/kg , salt extract pH is 6.9) (Fig. 1). The agrometeorological conditions of the agricultural year 2017-2018 can be characterized as not entirely favourable for the production of most field crops. The winter period was warmer than usual at 9.20 °C and the precipitation level exceeded the average multiyear level by 137.4 %. May was arid with a precipitation of 20.7 mm, which is 1.6 times lower than normal; July and August were warmer than usual at 3.1 and 1.40 °C, respectively. The precipitation during the summer was extremely unequal: June and August were extremely dry with the amount of precipitation in these months 2.1 and 3.6 times lower than the multiyear average.
Multispectral and thermal imaging (early in the morning and in the afternoon) surveys were carried out on the fields that make up the polygon 2 (area is 4 hectares): field 1winter wheat; field 2winter wheat + winter vetch; field 3barley; field 4spring hard wheat + barley + peas (Fig. 2).

Fig. 2. Polygon 2, 4 hectares area
The obtained values of vegetation indexes and temperature (after processing of aerial photography materials by a special digital platform AgroGIS by "GIS-R" Ltd.) must be accompanied by field statistics data to accumulate a database. On the basis of this database a digital survey is to be developed, which can convert the results in the form of vegetation indexes in real physical indicators (optimally in the on-line mode) [9].
To accomplish the task of collecting field statistics, a preliminary methodology was developed, which will be adjusted in 2019: 1. In each field, nine typical plots of 1*1 m in size (area is 1 m 2 ) are selected, significantly differing in the main indicators of agrophytocenosis (weed infestation, plant formation structure, pests and diseases). We place a marker in the centre of the plot: height is 1.50 m, the upper part (15 cm) is coloured red, and a number is assigned (Fig. 3). The coordinates of the marker are determined by WGS84 satellite navigation.
2. The following observations were conducted at each plot on the day of the flight session (either for 1-2, or maximum 3 days after the session): number of plants of the certain crop/plant mass of crops/number of weeds/weight of weeds/species composition of weeds/height of plants and ear length/pests and diseases/tillering capacity/color and condition of plants/photo. The data were filled in special tables during the measurement process directly in the field.
3. After full maturity of the crop plants, they were collected completely (from one m 2 ) in sheaves to process the yield structure, which was determined in the laboratory.
4. Collected field statistics were processed, systematized and accumulated in special tables in the form of databases.
5. The final stage of the work was binding of the collected field statistics to the results of imaging surveys (multispectral, thermal and optical).
According to coordinates on the electronic map, we determined the location of a specific marker, and then cut out a site around the marker, which was entered in a table with indicators of agrophytocenosis from this site.
As an example, we present the results of digital monitoring of field 2 (winter wheat + winter vetch) from polygon 2.
Indicators of agrophytocenosis on nine markers were determined: May 19 (phase is tillering of winter wheat), June 14 (bloom) and June 25 (milk stage) and July 24 (full ripeness, 7 days before harvesting).
During each session, the structure of winter wheat and a separate winter vetch was determined. The data were entered into different tables for each culture, and then a single table for two crops of this field was composed according to the results of a specific session ( Table 1).  (Table 2).
Field statistics, given in table 2, will serve as an integral part of the database, on the basis of which a digital platform will be developed to convert on-line vegetation indexes into specific physical parameters (in t/ha, mg/kg, etc.) of agrophytocenosis indicators for different crops in different vegetation phases in these soil climatic conditions.
According to the statistical data of the final table 2, the correlation coefficients between the biomass, total and the biomass of the grain and NDVI were calculated.
Calculations has shown a high correlation of the NDVI with physical indicators of biomass for all dates of aerial photography. However, the maximum correlation was noted in the milky stage of winter wheat, that is, on the date of June 25. So, for NDVI from May 19biomass, total, t/har = 0.71, biomass, grains, t/har = 0.79, for NDVI from June 14biomass, total, t/har = 0.67, biomass, grains, t/har = 0.70, for NDVI from June 25biomass, total, t/har = 0.85, biomass, grains, t/har = 0.82. Thus, the highest correlation coefficient between NDVI and the grain yield of winter wheat was found in the milk stage of winter wheat (r = 0.82).
In 2019, it is necessary to conduct a multispectral survey of winter wheat crops at the beginning of the milky-wax stage. We used the same method for digital monitoring of the remaining fields of polygon 2: field 1 winter wheat; field 3barley; field 4spring hard wheat + barley + peas.
The calculation of the correlation coefficients between the NDVI and the physical indicators of biomass, total and grain biomass, showed that the greatest correlation in field 1 was also noted in the milky stage of winter wheat (with grain biomass, r = 0.36); in the field 3in the phase the beginning of barley earing (with the biomass of grain, r = 0.97).
On the field 4 on the date of June 25, the spring hard wheat bloom and the beginning of the barley earing, a negative correlation was obtained between the NDVI and total biomass (r = -0.63) and grain biomass (r = -0.59). Negative correlation coefficients were obtained because of markers 3, 6, 7, and 8, which showed a significant predominance of weed biomass over plant biomass. Table 3 presents the different vegetation indexes for marker 6 of field 4 (spring hard wheat + barley + peas) on the date of June 25 (bloom of spring wheat, beginning of barley earing, pea beans filling).
Comparison of several vegetation indexes with field statistics of a specific plot of a field on a certain date of digital monitoring provides a unique opportunity to establish fundamental patterns, as well as to reveal a significant high correlation between relative index coefficients and real physical parameters of agrophytocenosis indicators like, for example, potential yield. The solution of this issue depends on the accumulation of a long-term database on soil-climatic zones, cultures and vegetation stages of crops [10].
Thus, the research conducted in 2018 revealed that the stated main goal to develop a digital monitoring methodology for agro-ecosystem indicators based on space and unmanned technologies with the conversion of the results of various types of digital aerial photography into real physical parameters of agrophytocenosis indicators is solvable.
It is necessary to continue accumulating the database as the basis for developing a digital platform for on-line conversion of vegetation indexes into specific physical parameters of the main indicators of single-species and poly-species agrophytocenosis for different crops and vegetation stages in these soil and climatic conditions in 2019 field season.
According to the results of 2018, the greatest correlation between the NDVI and the grain yield for winter and spring cereal crops was revealed in the milk stage (up to r = 0.82).
We also revealed the main optimal phases for the development of grain crops for digital monitoring: tillering; booting; bloom; the beginning of milk-wax stage.
We need to expand the range of learning crops with the following species: tilled crops (sunflower, grain maize), grain legumes (chickpea, lentil, soybeans, mung bean) and cruciferous (mustard, rape, oilseed radish).

Conclusion
The full analysis of the obtained digital (multispectral and thermal imaging, orthophoto, photo and video, as well as nitrogen content indexes and others) and field statistics can provide significant results to identify the fundamental patterns of functioning of single-species and poly-species crops in the framework of global climate changes, including increased climate aridity.