Next Article in Journal
Accurate Measurement of Tower Grounding Resistance for Single-Tower and Multi-Tower Parallel Scenarios Based on the Clamp Meter Method: For the Sustainable Operation of Towers
Previous Article in Journal
Place-Based Perspectives on Understanding the Value of Sustainable Heritage-Inspired Arts and Crafts in Jordan
Previous Article in Special Issue
Carbon Footprint of Main Grain Crop Production in Hubei and Jiangsu Provinces, 2005–2019
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Use of Indices Applied to Remote Sensing for Establishing Winter–Spring Cropping Areas in the Republic of Kazakhstan

1
Department of Remote Sensing, National Center for Space Research and Technology, 15 Shevchenko, Almaty 050010, Kazakhstan
2
Department of Land Resources and Cadastre, Kazakh National Agrarian Research University, Valikhanov St 137, Almaty 050000, Kazakhstan
3
Department of Cartography and Geoinformatics, Al-Farabi Kazakh National University, 71 al-Farabi, Almaty 050010, Kazakhstan
4
School of Information Technology and Engineering (SITE), Kazakh British Technical University, Almaty 050010, Kazakhstan
5
Department of Geological and Environmental Sciences, Western Michigan University, Kalamazoo, MI 49008, USA
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(17), 7548; https://doi.org/10.3390/su16177548
Submission received: 16 June 2024 / Revised: 17 August 2024 / Accepted: 21 August 2024 / Published: 31 August 2024
(This article belongs to the Special Issue Farmers’ Adaptation to Climate Change and Sustainable Development)

Abstract

:
Farmers in Kazakhstan face unreliable water resources. This includes water scarcity in the summer, high fluctuations in precipitation levels, and an increase in extreme weather events such as snow, rain, floods, and droughts. Wheat production is regulated and subsidized by the Kazakh government to strengthen food security. The proper monitoring of crop production is vital to government agencies, as well as insurance and banking structures. These organizations offer subsidies through different levels support. Some farmers already use farmland soil monitoring combined with adaptive combinations of different crops. These include winter–spring plowing crop programs. Winter wheat crops are generally more adaptive and may survive summer droughts. Kazakhstan is a large country with large plots of farmland, which are complicated to monitor. Therefore, it would be reasonable to adapt more efficient technologies and methodologies, such as remote sensing. This research work presents a method for identifying winter wheat crops in the foothills of South Kazakhstan by employing multi-temporal Sentinel-2 data. Here, the researchers adapted and applied a Plowed Land Index, derived from the Brightness Index. The methodology encompasses satellite data processing, the computation of Plowed Land Index values for the swift recognition of plowed fields and the demarcation of winter wheat crop sowing regions, along with a comparative analysis of the acquired data with ground surveys.

1. Introduction

In modern agriculture, remote sensing technologies, including satellite imagery, are becoming increasingly useful for monitoring and managing agricultural crop production. These technologies are especially valuable for large agricultural lands, aiding farmers and researchers in optimizing yields and ensuring sustainable farming practices. The application of remote sensing technologies in various countries has demonstrated significant advantages. For example, in the United States, satellite images are used to monitor corn and soybean production, providing crucial data on crop conditions and facilitating the timely application of fertilizers and pesticides [1,2]. In Europe, remote sensing supports precision farming by offering high-resolution data that help in efficient land management involving the production of diverse crop types [3,4]. Indian researchers have used satellite data for crop area mapping and production estimation, supporting food security and agricultural planning [5,6,7]. China is advanced in satellite data applications for crop monitoring and yield forecasts [8,9,10]. In Brazil and Argentina, remote sensing is utilized to assess pasture conditions and manage rangelands [11,12]. Russia has intensively applied satellite data for analyzing land-use changes and monitoring forest resources [13,14]. In South Africa, remote sensing is used for vegetation monitoring and drought risk assessment [15,16].
Despite the proven benefits of remote sensing in agriculture, its application in Kazakhstan, particularly for monitoring winter–spring plowing crops, remains underexplored. Kazakhstan, with its vast agricultural landscapes, can benefit from integrating satellite technologies into agricultural practices. The Sentinel-2 satellite, with its high resolution and frequent revisit times [17,18,19] represents a valuable resource for agrarian land monitoring and land use change detection. This study focuses on the application of Sentinel-2 imagery for identifying arable lands and monitoring the timing of the sowing and growth of winter–spring plowing crops in the Turkestan region of Kazakhstan. Appropriate cooperation and remote sensing application technologies for farmland are required. Wheat production is one of the main forms of agricultural crop production in Kazakhstan, and it is regulated and subsidized by the Kazakh government. Monoculture crops often become infected with different types of virus. Moreover, Kazakhstan has very unstable water availability: in some periods, large floods cover planning areas, while in other periods, substantial droughts may occur. Appropriate strategies are necessary to adapt to mitigate these issues. Crop monitoring on large plantation sites would help to resolve these issues.
Winter crops constitute the main activity of Kazakh farmers in the foothill zone. The area sown with winter grain crops varies significantly from year to year, depending on the crop rotations carried out on farms and weather conditions. The development of a method for automatically recognizing the sowing mask of winter grain crops based on remote sensing is the basis for solving various problems in space monitoring and directly affects the condition and yield of winter grain crops and the timing and volume of various types of field work on farms.
Using Sentinel-2 images from September to November 2021 and March to May 2022, along with data from the autumn and spring periods of 2017–2018, the Tasseled Cap Brightness Index [20,21] was processed with the Plowed Land Index, specifically designed to identify arable lands in the foothill zone of South Kazakhstan.
Processed satellite data were used for the winter–spring plowing crop classifications, along with field surveys for the calibration and verification of the processed works. This work mainly emphasizes winter wheat tracking, based on the requests of farmers from the Turkestan region.
Remote sensing technology applications allow the monitoring of large farming areas of soil and crop conditions, in connection with subsidy-support systems provided by state agencies. Through the timely implementation of agricultural measures, such as fertilization and plant protection, the control of land productivity supports the sustainable development of crop yields and soil conditions. Digitized data on the conditions of farm fields allows the more effective planning and management of resources, including water, fertilizers and pesticides. The optimization of winter grain production using remote sensing makes it possible to obtain monitoring data more efficiently. The development of modern agricultural technologies and the need for their maintenance create advancements in remote sensing technology applications.
The purpose of this research is to develop and use remote-sensing-based tools, which makes it possible to efficiently monitor the sowing areas of different crop zones. These efforts will be used to rapidly investigate agricultural production and solve subsequent problems in crop production. The practical applications of these efforts are for the benefit of agricultural producers, as well as land administrations, the Kazakh government, insurance organizations, and banking institutions providing subsidies. The results obtained contribute to optimizing the production of the winter grain crops, increasing the productivity of agricultural land, improving the sustainable development of rural areas and ensuring food security.

2. Materials and Methods

2.1. Study Area

The territory of the Turkestan region boasts diverse landscapes, soils and agro-climatic resources. Significant areas (up to half of the region’s territory) are occupied by plains and sands in the north and west of the region. In the middle of the region, from northwest to southeast, the Karatau Range runs, with most peaks around 1800 m, separating the plains territory from the mountainous regions. Low-mountain and high-mountain ranges of the Western Tien Shan (with elevations ranging from 1500–2000 to 2500–3500 m and more) are in the southeast of the region. The diversity of the region’s relief has resulted in significant heterogeneity of soil-climatic conditions and agricultural specialization in different areas [22]. The cultivation of winter grain crops is concentrated in the foothills of the Turkestan region.
For this study, the Tyulkubas district, located in the foothill zone in the eastern part of the Turkestan region, bordering the Zhambyl region, was selected (Figure 1).
The agricultural lands of the Tyulkubas district cover an area of 148.3 thousand hectares. This includes arable land, perennial plantations, hayfields and pastures. According to statistical data for 2022, arable lands account for 64.5 thousand hectares, of which fodder crops (predominantly alfalfa) occupy 27.6 thousand hectares (43%), winter wheat crops 16.9 thousand hectares (26%), oilseed crops 12.1 thousand hectares (19%), spring wheat crops 5.9 thousand hectares (9%) and vegetables and melons 1.6 thousand hectares (3%). Thus, winter wheat and fodder crops dominate the cropping structure of the district. The cropping areas of agricultural crops vary from year to year, depending on weather conditions and crop rotation.
The Tyulkubas district possesses unique geographical and soilclimatic characteristics that significantly influence agricultural activities and field productivity. Depending on the vertical zonation and decreasing altitude of mountain ranges, a change in soil type occurs from mountain meadow and subalpine soils (the most fertile, with a high humus content of up to 10%) to mountain brown soils (moderately fertile) and chestnut soils (less fertile, with low humus contents of 3–4%).
Land use in the Tyulkubas district is in the dry foothill zone and in the mountainous agroclimatic zones.
The dry foothill zone is characterized by precipitation amounts of 100–250 mm, and the highest moisture during the vegetation period (HTC ranging from 0.3 to 0.5). The vegetation period is characterized by a prolonged period of warm days (up to 160–180 days per year). This zone provides the most favorable agro-climatic conditions for cultivating various crops (grains, oilseeds, fodder, vegetables, fruit and berry crops), primarily under rainfed conditions, sometimes under irrigation.
The mountainous agroclimatic zone has fewer thermal resources, a shorter growing season and frost-free period (160 days or less) and increased moisture. The territory is primarily used as pasture and hayfields, and cultivation mainly involves cold-resistant crops (grains, fodder, etc.) under rainfed conditions. Due to limited thermal resources, the region is not suitable for the development of irrigated agriculture.
Precipitation during the autumn period is crucial for seedling emergence and wintering of winter crops. In the foothill and mountainous areas, precipitation during the autumn period can reach 100–200 mm, significantly more than on the plains. By the end of the autumn period, winter wheat crops usually reach the third leaf stage or the beginning of tillering, which, with sufficient snow cover, contributes to good crop wintering.
During the winter period, precipitation in the form of snow falls in the foothill and mountainous areas, reaching thicknesses of 20–30 to 35–40 cm, effectively protecting crops from freezing. Vegetative growth of winter wheat crops typically resumes in the second half of March, with active growth and development observed throughout April and May. Heading and flowering usually occur in the third week of May, followed by maturation in June. Harvesting of winter wheat crops, depending on the ripening dates, usually begins in mid-July and continues until the end of August, sometimes extending into early September. The agro-climatic characteristics of winter wheat cultivation and the timing of field operations were used as a basis for satellite monitoring of cropping areas and the development of a methodology for determining the timing and areas of winter wheat crop sowing in the mountainous region [23].

2.2. Satellite and Ground Data

To address the research objectives, multispectral data from the Sentinel-2 satellite system were utilized, obtained using the Level-1C (L1C) optical multispectral instrument (MSI). These data were downloaded from ESA’s Copernicus Open Access Hub, which provides access to Sentinel-2 satellite data.
The Sentinel-2 satellite mission, comprising two satellites, Sentinel-2A and Sentinel-2B, began on 23 June 2015, with the launch of Sentinel-2A. This mission is characterized by a wide swath, and high resolution, and it is based on optical multispectral visualization. Sentinel-2A and Sentinel-2B orbit the Earth in the same orbit, with a phase difference of 180 degrees, providing high image repeatability, with a 5-day revisit period.
The primary objective of the Sentinel-2 mission is to provide multispectral images for various applications in different domains, such as land use monitoring (e.g., agriculture, vegetation, forests, water bodies), emergency management (e.g., monitoring and managing natural disasters), climate research (e.g., glacier and coastal monitoring) and security provision (e.g., border and maritime surveillance).
The Multispectral Instrument (MSI) on Sentinel-2A satellites comprises 13 spectral bands (VIS, NIR, SWIR) with varying spatial resolutions of 10 m, 20 m and 60 m for different spectral bands. The satellites’ orbit is sun-synchronous, at an average altitude of 786 km, allowing for a clear 5-day revisiting cycle between the two satellites [24,25,26].
For this study, Sentinel-2 images were processed using ArcGIS software, version 8 (provided by ESRI). The list of Sentinel-2 data used, for example, for the years 2021–2022, is provided in Table 1. The coverage scenes of the Sentinel-2 imagery over the Tyulkubas District are illustrated in Figure 2.
The images named in Table 1 are the standard identifiers for the Sentinel-2 satellite datasets. Each identifier consists of several parts that provide information about the location and date of acquisition. For example, in T42TWN_20210913, T42TXN_20210928, T42TXM_20210918 images, the first letter, T, indicates that this is a level 1C image, atmospherically corrected and ready for use.
The next number, 42, is related to the UTM (Universal Transverse Mercator) zone code.
The TWN, TXN and TXM are related to the different areas within the UTM zone, including the letter T, the latitude band strip within the zone 42. The two-letter codes WN, XN and XM are the territory squares within the UTM zone 42, designed by the MGRS (Military Grid Reference System). Each MGRS square (WN, XN, XM) is 100 × 100 km within UTM zone 42. The left part of the combination of letters W, X and X represents the square image’s position on the longitude axis. The right part of the combination of letters N, N and M represents the square image’s position on the latitude axis.
WN: W—indicates a specific longitude position (east–west) within the zone. N—indicates a latitude position (north-south) within the zone.
XN: X—another longitude square adjacent to square “W”. N—the same latitude position as square “WN”.
XM: X—the same longitude square as “XN”. M—a square located further south in latitude compared to “N”.
The square codes WN, XN and XM represent different areas within the same UTM zone 42. These squares are necessary to indicate exactly which part of the Earth’s surface was captured by the Sentinel-2 satellite. Such designations allow precise identification of the territory covered by a specific image.
The final numbers represent the image collected in the format year–month–day, YYYYMMDD, such as 20210913 (year 2021, September 13), 20210928 (year 2021, September 09) and 20210918 (year 2021, September 18).
Within each territory, square scenes (WN, XN, XM) have differences in cloudiness. The scenes without cloudiness were used for the research applications. The WN scenes were without cloudiness. Some XN and XM scenes had high cloudiness and were excluded from the data processing, which is why some columns within the XN and XM scenes are empty.
The use of Sentinel-2 mission data provided an effective tool for analyzing changes on the land surface in the Tyulkubas district, allowing us to obtain valuable information about plowed lands, sown areas and the condition of winter grain crops to achieve our research objectives.
For the analysis, calibration and verification of the processed satellite data results, the ground farmlands survey data sets were used, covering various agricultural plots, including carefully selected test GPS (global positioning system) points. The map of field survey routes and test farms in the Tyulkubas district of the Turkestan region during the period from 2 June to 7 June 2022, is presented in Figure 3. The data were collected from the farms “Bratya”, “Makulbek” and “Kyzylzhar”. The farmers of these lands provided the related field support and access to their lands. Farmers provided important information on land use, cultivation of agricultural crops, sowing and harvesting dates and explanations of their important agricultural activities, as well as problems in their fields.
The route surveys were conducted from 2 June to 7 June 2022 (160 GPS points) and during April–May 2018 (over 450 GPS points, 188 GPS points for winter and spring grain crops). The survey results included determining the crop type in the field, GPS observation coordinates, photography, developmental stage and condition of winter and spring grain crops, with fall and spring plowing activities. Visual differences in the spring conditions of fields of the main cultivated agricultural crops during the April 2018 survey in the foothill zone of the region are illustrated, for example, in Figure 4 and Figure 5. During the spring period (early April), the fields of winter grain crops stand out well due to their possessing the most developed green biomass and plant height among the general array of fields, compared to other spring-planted crops. In contrast to winter grain crops, fields of perennial alfalfa also have sufficiently developed biomass, but their tonality is lighter greenish, and there are no pronounced plant rows in the fields in spring. In fields of spring grain crops, only seedlings emerge, so these fields have a weak greenish tonality. Mass field work, as well as plowing and sowing of safflower are actively underway in safflower fields during this period, and in some early fields, the first seedlings of safflower appear (Figure 5d).
A combination of indices based on satellite data was used to identify winter crops.

2.3. Methods

To analyze changes in brightness on the Earth’s surface and extract key spectral components, the Tasseled Cap transformation was applied to multi-zone Sentinel-2 satellite data in the study. The Tasseled Cap transformation is a mathematical algorithm developed by R. J. Kauth and G. S. Thomas [27], which allows for the extraction of primary spectral components from multi-zone satellite images such as Sentinel-2 data. The Tasseled Cap transformation involves a set of coefficients applied to each spectral channel of the image to derive new components characterizing Earth’s surface features. As a result of the Tasseled Cap transformation, three key indices are obtained: the Brightness Index, the Greenness Index and the Wetness Index [28,29,30,31].
The Brightness Index derived from the Tasseled Cap transformation is a linear combination of values from various spectral bands of the image. The linear combination works as follows: the values of each spectral band are multiplied by their respective weighting coefficient, and then the products are summed together. The formula of the Brightness Index (BI) from the Tasseled Cap transformation for Sentinel-2 data is as follows [21]:
BI = 0.3037 × b2 + 0.2793 × b3 + 0.4743 × b4 + 0.5585 × b8 + 0.5082 × b11 − 0.1863 × b12,
where
  • b2—is the blue band of the Sentinel-2 sensor;
  • b3—is the green band;
  • b4—is the red band;
  • b8—is the near-infrared band;
  • b11—is the shortwave infrared band 1;
  • b12—is the shortwave infrared band 2.
The Plowed Land Index (PLI) is derived from the Brightness Index of the Tasseled Cap transformation. The formula for the Plowed Land Index represents an original linear combination of values from six Sentinel-2 spectral bands: B2, B3, B4, B8, B11 and B12. In the process of forming this index, each spectral band was carefully weighted and transformed using specific weighting coefficients. The signs in the formula were chosen to optimize the detection of plowed lands in the foothill zone of the Turkestan region. As a result, a new linear combination with calculated weighting coefficients for detecting plowed lands in the studied region was obtained, as follows:
PLI = −0.0037 × b2 − 0.1793 × b3 + 0.5403 × b4 − 0.5585 × b8 − 0.1082 × b11 + 0.3013 × b12.
Plowed lands are known to be characterized by low positive values. For instance, NDVI values for plowed fields typically range from 0.00 to 0.10 [32,33]. In the created PLI, its minimum positive values (0.00–0.01) also indicate the presence of plowed fields, while higher index values (0.10–0.50) indicate fires [34,35]. For graphical analysis, the following land surface types were selected: plowed fields, bare soil, dense and sparse vegetation and fires, with 5 fields of each type analyzed using 30 pixel points within each field, resulting in a total sample size of 750 points. The graphical representation of the pixel analysis of the PLI for different types of objects varied from −0.180 to 0.260 (Figure 6). The graph shows the lowest negative values of the Plowed Land Index (from −0.180 to −0.130), indicating dense vegetation, values from −0.080 to −0.040, indicating sparse vegetation, values from −0.020 to −0.006, indicating bare soil, values from 0.004 to 0.030, indicating plowed soil and values from 0.050 to 0.260, indicating fires. Thus, small positive values of the PLI allow for a clear distinction of plowed territories from the general set of analyzed objects.
Comparative results of identification of plowed fields using the Plowed Land Index and NDVI for 11 November 2021 are presented in Figure 7. The analysis of these indices showed that the most effective tool for delineating fields with autumn plowing is the Plowed Land Index (0.00–0.010), which clearly defines the boundaries of plowed lands and accurately determines their quantity and area. Determining plowed fields based on NDVI data (0.00–0.10) does not always coincide with actual plowing, as it less distinctly defines the boundaries of plowed fields. Some fields with NDVI values (0.00–0.10), as shown in Figure 7, were found to be without plowing, either in a post-harvest state or with sparse vegetation of perennial grasses or weeds. The delineation of plowed fields using the Plowed Land Index in autumn and spring most clearly defines the boundaries of plowed fields (Figure 8).
Indeed, the analysis of surface brightness changes using the Sentinel-2 images’ Plowed Land Index enables accurate identification of plowed areas in fields and their boundaries in foothill and mountainous regions.

2.4. Methodology

The methodology for recognizing the final mask of winter wheat crop sowing is based on the utilization of Sentinel-2 imagery, along with analysis of the dynamics of calculated values of the Plowed Land Index and NDVI in the research area section, covered by the collected Satellite datasets (see Figure 9). The approach involves identifying spectral changes during autumn and spring associated with soil cultivation and crop sowing, and considering the timing of field operations in the studied region.
In the foothill regions of southern Kazakhstan, autumn soil cultivation and sowing of winter wheat crops typically commence in the second half of September and continue until November. The pace and quality of autumn fieldwork are largely determined by soil moisture and precipitation [36]. However, on some fields, only soil cultivation for spring sowing of spring crops may occur in autumn. Therefore, not all plowed fields in autumn may be sown with winter wheat crops [37]. This factor has been considered in the development of the methodology for automated recognition of the mask of winter wheat crops based on satellite information.
Figure 10 illustrates the visualization of spectral changes occurring in the array of fields in the Tyulkubas district during the autumn and spring periods using Sentinel-2 imagery. A combination of channels, including SWIR2 (short-wave infrared 2), NIR (near infrared) and the green channel, was utilized for the analysis. The results of autumn plowing of fields in Figure 10a are presented as a vector layer, where plowed fields are depicted in shades of brown. Fields with shades of green characterize the absence of autumn soil cultivation or sparse vegetation.
In the spring period (late March–early April), vegetation resumes for overwintering agricultural crops (winter wheats, perennial grasses), while in other fields, spring plowing is carried out for sowing spring crops. Visualization of the spring changes in the same field array using Sentinel-2 images is presented in Figure 10b. Here, spring plowing of the fields is depicted in shades of pink, while fields with shades of green indicate spring vegetation of winter wheats and perennial grasses. Additionally, fields of winter wheat crops in spring have a darker green tone compared to fields of perennial grasses (alfalfa), which were plowed in autumn. In contrast to winter wheat crops, plowing is not observed in autumn for fields of perennial seeded grasses. These autumn and spring agronomic features of field conditions allow for their effective use in recognizing fields of winter wheat crops.
Sentinel-2 satellite data for the autumn (September–November) and spring (April–May) periods provide multi-zone images of medium resolution with high repeatability. Atmospheric correction was applied to all satellite images to eliminate the influence of atmospheric conditions on the images. Subsequently, the agricultural land field masks were digitized based on the Sentinel-2 images.
According to the methodology, the first step involves calculating the Plowed Land Index, which is an effective tool for assessing and analyzing changes in agricultural land. Subsequently, based on the identified ranges of the Plowed Land Index, the detection of autumn plowing fields was carried out in the territory of the Tyulkubas district (Figure 11).
The next stage of the methodology involves determining the spring plowing of fields for sowing spring crops based on the Plowed Land Index and excluding such fields from the mask with autumn plowing. The remaining satellite mask of fields with only autumn tillage will characterize the mask of sown winter wheats.
The methodology involves additional correction of the satellite mask for winter wheat crops based on the analysis of the increase in NDVI values in the spring period (Figure 12). To confirm the vegetation of winter wheat crops in April-May, the Normalized Difference Vegetation Index (NDVI) was calculated. Fields with winter crops in the spring period were characterized by a noticeable increase in NDVI values, while for other fields with spring sowing, NDVI values remained low. Fields plowed in autumn with very low NDVI values in spring were excluded from the winter wheat crop mask. Therefore, the active growth of NDVI values in the spring period for fields plowed in autumn served as an additional tool confirming that winter wheat crop fields were correctly identified.
The results of determining the mask of winter wheat crops using satellite data involved verifying the satellite data against ground-based survey materials and actual information from individual farms. In Figure 10, all fields identified based on satellite data are confirmed by ground-based survey data. Thus, the developed methodology allows for the successful and efficient recognition of the final mask of winter wheat crop sowing based on satellite information in real-time operation.

3. Results

Based on the developed methodology, maps of the distribution of plowed lands in the Tyulkubas district of the Turkestan region were constructed for the autumn period of 2021 and the spring period of 2022, with separate highlighting of fields where autumn and spring plowing was observed (Figure 13), as well as a map of the location of the final areas of sown winter wheat crops in 2022 (Figure 14).
The advantage of using satellite data for monitoring tasks is the ability to obtain standardized operational reports in the form of maps, tables and graphs, reflecting the dynamics of the main agricultural activities on fields at various administrative levels. As an example, Table 2 illustrates the determination of the areas and progress of field plowing based on calculated values of the Plowed Land Index in the Tyulkubas district for the autumn period of 2021. The graphical data in Figure 15 show that during the period from September to October, the main plowing of fields was carried out on 75% of the area, in November on 23% of the area, and in the first half of December on a negligible area of 2%.
The investigated remote sensing datasets included the Autumn 2021 and Spring 2022 periods. The Autumn period of 2021 was the main investigated period of time, shown in Table 2 and Figure 15, since this was when the main agrotechnical works were provided for the winter crops. Spring 2022 was used to exclude other spring crops, which is why they are now included in Table 2 and Figure 15. The main research goal was dedicated to the Autumn 2021 period’s winter crops. Figure 15 is a bar chart illustrating the percentage change in plowing during the Autumn period of 2021. The chart displays the data by period to highlight the key stages and intensity of plowing activities: from September 5 to September 13, 29% of the total plowing work activities took place, from September 14 to September 26, 14% of the total plowing work occurred, and so on for the following periods. This approach is important in terms of the period of time, days and how long the plowing works lasted. In southern Kazakhstan, it is critical to complete the autumn plowing by the end of October, with some flexibility in terms of days, depending on weather conditions, soil moisture and temperature. This is necessary to ensure that winter crops can germinate before the onset of winter. Successful germination ensures the deep rooting of plants, which contributes to their good overwintering and, subsequently, has a positive effect on crop yields. If plowing and sowing are carried out outside the optimal period, this leads to poor wintering and, as a consequence, a decrease in yield in the spring. Thus, the data presented in Figure 15 allow us to monitor and analyze the periods of the main plowing activity and their compliance with the agronomic requirements for the optimal growth of winter crops in South Kazakhstan.

4. Discussion

The calibration and verification of the processed satellite data results were carried out using the field survey data from individual farms in the Tyulkubas district (“Bratya”, “Makulbek”, and “Kyzylzhar’’ farms). According to the farm data, 68 fields were involved in the verification process, of which 22 were fields of winter grain crops and 46 were fields of other agricultural crops. All the fields of winter grain crops identified by the satellite data were confirmed by agricultural formations (Figure 16). The verification also involved data from route surveys and the results of expert interpretations of Sentinel-2 satellite imagery, as well as official statistical data.
For satellite processed data output verification, individual field arrays from route surveys were used (Figure 17).
For these purposes, 60 fields were selected, out of which fifty-six fields with winter wheat (more than 93%) showed matching results between the satellite and ground data, with only four fields having discrepancies, mostly due to small areas. The verification of the satellite-derived winter wheat crop mask based on the Plowed Land Index included a comparison with results from expert interpretation of the images. For this purpose, the areas of autumn and spring plowing fields and the winter wheat crop mask were calculated using ArcGIS 10.8 software. As a result, the area of autumn plowing fields according to the developed methodology in the Tyulkubas district in 2022 amounted to 22,108.3 hectares. For some fields, re-plowing was conducted in spring, covering an area of 2556.4 hectares, and for a small number of fields, covering 90.0 hectares during the April–May period, there was no increase in NDVI or a very weak increase, which is atypical for winter crops. Consequently, according to the methodology, the total area of fields to be subtracted from the mask of fields with autumn plowing amounted to 2646.4 hectares. Taking this into account, the final area of the winter wheat crop mask, according to the developed technology, was 19,461.9 hectares. Additionally, visual expert interpretation of the satellite images was conducted for verification. The area of winter wheat according to the expert interpretation was 19,020.8 hectares. The difference between the areas according to the developed technology and expert interpretation was 441.1 hectares (2.3%). The official statistical data indicate that the area under winter grain crops in the Tyulkubas district in 2022 amounted to 16,885.5 hectares. The satellite results exceed the statistical data by 2576.4 hectares (15.3%). As is known, actual crop areas often differ slightly from official statistical data, a discrepancy that was confirmed by our results in this case. Farmers from the Tyulkubas district are also reviewing to develop more groundwater use programs using groundwater wells and climate change precipitation prediction analyses, which will require additional research programs, similar to those our research group’s colleagues developed in the neighboring Zhambul [38,39] and Almaty [40] regions.

5. Conclusions

The developed methodology, based on the RS applications, analysis of the Plowed Land Index (PLI) for various fields in the autumn and spring periods and ground information data, can be a supportive tool for the rapid automatic monitoring of crop production, including winter grain crops. The efficient approach involving RS data processing for the farming areas and timing of sowing, as well as the generation of maps and tabular data in the context of regions, farms and individual fields, allows researchers to identify the depression areas in fields and carry out timely agro-technical measures to improve the condition of crops. Wheat production is one of the main forms of agricultural crop production in Kazakhstan, and it is regulated and subsidized by the Kazakh government. The determination of the areas for cultivating fields and sowing crops on farms is based on taking into account the field area according to their land use map, which has not changed on farms for many years. Data from the space monitoring of land cultivation show that there is a change in the area of individual fields due to the development of new fields or changes in the configuration of fields (adding or reducing areas). Space data can be used to record all the changes occurring in fields in real time. The determination of the actual sowing mask for winter grain crops in a particular year can serve as the basis for solving further problems in the space monitoring of grain crops (assessing the condition, yield and volume of agricultural work in fields). Satellite data will allow the real-time assessment of the actual areas and timings of the sowing of winter grain crops, which will later be used to assess the conditions and carry out agricultural work in depressed areas of fields. These efforts will also help farmers in communications and ensure efficiency in the provision of information to government agencies, as well as insurance and banking institutions, in order to receive the proper level of subsidy support. The obtained results of the verification of the Sentinel-2 satellite data confirm the accuracy and reliability of using the arable land index (PLI) in identifying winter wheat planting areas. The developed methodology for using satellite data for crop monitoring can be expanded for larger regions in Kazakhstan, which is a major agricultural crop producer.
Farmland insurance institutions, banking companies and land developers have shown the most interest in expanding this methodology to Kazakhstan’s other regions.

Author Contributions

All authors contributed to the concept and design of the study. Preparation of materials, collection and analysis of data were performed by A.A., N.K. and J.S., research supervision, negotiations with and support for farmers were provided by M.N. and R.B. and English editing support was provided by R.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the local communities and co-authors own expenses.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The field survey farmland data presented in this study were collected with approval from Kazakhstan’s local authorities’ and farmers’ support and are available on request from the corresponding author.

Acknowledgments

Research authors acknowledge and are very thankful to the local farms, “Brothers”, “Makulbek” and “Kyzylzhar” from Tyulkubas district, South Kazakhstan, for the field survey support, their openness to sharing data and their interest in expanding cooperation, to continue these research efforts.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Doraiswamy, P.C.; Hatfield, J.L.; Jackson, T.J.; Akhmedov, B.; Prueger, J.; Stern, A. Crop condition and yield simulations using Landsat and MODIS. Remote Sens. Environ. 2004, 92, 548–559. [Google Scholar] [CrossRef]
  2. Peng, Y.; Nguy-Robertson, A.; Arkebauer, T.; Gitelson, A.A. Assessment of Canopy Chlorophyll Content Retrieval in Maize and Soybean: Implications of Hysteresis on the Development of Generic Algorithms. Remote Sens. 2017, 9, 226. [Google Scholar] [CrossRef]
  3. Atzberger, C. Advances in remote sensing of agriculture: Context description, existing operational monitoring systems and major information needs. Remote. Sens. 2013, 5, 949–981. [Google Scholar] [CrossRef]
  4. Immitzer, M.; Vuolo, F.; Atzberger, C. First Experience with Sentinel-2 Data for Crop and Tree Species Classifications in Central Europe. Remote Sens. 2016, 8, 166. [Google Scholar] [CrossRef]
  5. Panigrahy, S.; Sharma, S. Mapping of crop rotation using multidate Indian Remote Sensing Satellite digital data. ISPRS J. Photogramm. Remote Sens. 1997, 52, 85–91. [Google Scholar] [CrossRef]
  6. Yogita, Y.; Mohite, J.; Jagyasi, B. Crop classification based on multi-temporal satellite remote sensing data for agro-advisory services. In SPIE Asia-Pacific Remote Sensing; SPIE-The International Society for Optical Engineering: Bellingham, DC, USA, 2014; p. 926004. [Google Scholar] [CrossRef]
  7. Agilandeeswari, L.; Prabukumar, M.; Radhesyam, V.; Phaneendra, K.L.N.B.; Farhan, A. Crop Classification for Agricultural Applications in Hyperspectral Remote Sensing Images. Appl. Sci. 2021, 12, 1670. [Google Scholar] [CrossRef]
  8. Liu, H.; Yang, X.; Wang, N. Remote sensing-based estimation system for winter wheat yield in North China Plain. Chin. Geogr. Sci. 1999, 9, 40–48. [Google Scholar] [CrossRef]
  9. Zhang, L.; Zhang, Z.; Luo, Y.; Cao, J.; Tao, F. Combining Optical, Fluorescence, Thermal Satellite, and Environmental Data to Predict County-Level Maize Yield in China Using Machine Learning Approaches. Remote Sens. 2019, 12, 21. [Google Scholar] [CrossRef]
  10. Fu, Y.; Huang, J.; Shen, Y.; Liu, S.; Huang, Y.; Dong, J.; Han, W.; Ye, T.; Zhao, W.; Yuan, W. A Satellite-Based Method for National Winter Wheat Yield Estimating in China. Remote Sens. 2020, 13, 4680. [Google Scholar] [CrossRef]
  11. Numata, I.; Roberts, D.A.; Chadwick, O.A.; Schimel, J.P.; Galvão, L.S.; Soares, J.V. Evaluation of hyperspectral data for pasture estimate in the Brazilian Amazon using field and imaging spectrometers. Remote. Sens. Environ. 2008, 112, 1569–1583. [Google Scholar] [CrossRef]
  12. Kurtz, D.B.; Schellberg, J.; Braun, M. Ground and satellite-based assessment of rangeland management in sub-tropical Argentina. Appl. Geogr. 2010, 30, 210–220. [Google Scholar] [CrossRef]
  13. Potapov, P.; Turubanova, S.; Hansen, M.C. Regional-scale boreal forest cover and change mapping using Landsat data composites for European Russia. Remote. Sens. Environ. 2011, 115, 548–561. [Google Scholar] [CrossRef]
  14. Shabanov, N.V.; Egorov, V.A.; Miklashevich, T.S.; Stytsenko, E.A.; Bartalev, S.A. Utility of Leaf Area Index for Monitoring Phenology of Russian Forests. Remote. Sens. 2022, 15, 5419. [Google Scholar] [CrossRef]
  15. Graw, V.; Ghazaryan, G.; Dall, K.; Gómez, A.D.; Abdel-Hamid, A.; Jordaan, A.; Piroska, R.; Post, J.; Szarzynski, J.; Walz, Y.; et al. Drought Dynamics and Vegetation Productivity in Different LandManagement Systems of Eastern Cape, South Africa—A Remote Sensing Perspective. Sustainability 2017, 9, 1728. [Google Scholar] [CrossRef]
  16. Mukhawana, M.B.; Kanyerere, T.; Kahler, D. Review of In-Situ and Remote Sensing-Based Indices and Their Applicability for Integrated Drought Monitoring in South Africa. Water 2022, 15, 240. [Google Scholar] [CrossRef]
  17. Drusch, M.; Del Bello, U.; Carlier, S.; Colin, O.; Fernandez, V.; Gascon, F.; Hoersch, B.; Isola, C.; Laberinti, P.; Martimort, P.; et al. Sentinel-2: ESA’s Optical High-Resolution Mission for GMES Operational Services. Remote Sens. Environ. 2012, 120, 25–36. [Google Scholar] [CrossRef]
  18. Phiri, D.; Simwanda, M.; Salekin, S.; Nyirenda, V.R.; Murayama, Y.; Ranagalage, M. Sentinel-2 Data for Land Cover/Use Mapping: A Review. Remote Sens. 2019, 12, 2291. [Google Scholar] [CrossRef]
  19. Jiang, Y.; Lu, Z.; Li, S.; Lei, Y.; Chu, Q.; Yin, X.; Chen, F. Large-Scale and High-Resolution Crop Mapping in China Using Sentinel-2 Satellite Imagery. Agriculture 2020, 10, 433. [Google Scholar] [CrossRef]
  20. Lamqadem, A.A.; Saber, H.; Pradhan, B. Quantitative Assessment of Desertification in an Arid Oasis Using Remote Sensing Data and Spectral Index Techniques. Remote Sens. 2018, 10, 1862. [Google Scholar] [CrossRef]
  21. Shi, T.; Xu, H. Derivation of Tasseled Cap Transformation Coefficients for Sentinel-2 MSI At-Sensor Reflectance Data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2019, 12, 4038–4048. [Google Scholar] [CrossRef]
  22. Gidrometeoizdat, Agroclimatic Resources of the Chimkent Region of the Kazakh SSR; USSR State Committee for Hydrometeorology and Natural Environment Control, Administration of the Hydrometeorological Service of the Kazakh SSR; Alma-Ata Hydrometeorological Observatory. Gidrometeoizdat: Leningrad, Russia, 1979. (In Translation)
  23. Gringof, I.; Popova, V.; Strashny, V. Agrometeorology; Gidrometeoizdat: Moscow, Russia, 1987. [Google Scholar]
  24. European Space Agency (ESA). Sentinel-2 User Handbook, Revision 2, ESA Standard Document; European Space Agency (ESA): Paris, France, 2015. [Google Scholar]
  25. Hagolle, O.; Sylvander, S.; Huc, M.; Claverie, M.; Clesse, D.; Dechoz, C.; Lonjou, V.; Poulain, V. SPOT-4: Simulation of Sentinel-2 time series on 45 large sites. Remote. Sens. 2015, 7, 12242–12264. [Google Scholar] [CrossRef]
  26. Segl, K.; Guanter, L.; Gascon, F.; Kuester, T.; Rogass, C.; Mielke, C. S2eteS: An End-to-End Modeling Tool for the Simulation of Sentinel-2 Image Products. IEEE Trans. Geosci. Remote. Sens. 2015, 53, 5560–5571. [Google Scholar] [CrossRef]
  27. Kauth, R.; Thomas, G.S. The tasseled cap-a graphic description of the spectral-temporal development of agricultural crops as seen by Landsat. In Proceedings of the Symposium on Machine Processing of Remotely Sensed Data, West Lafayette, IN, USA, 29 June 29–1 July 1976; pp. 4B41–4B51. [Google Scholar]
  28. Baig, M.H.A.; Zhang, L.; Shuai, T.; Tong, Q. Derivation of a Tasselled Cap Transformation Based on Landsat 8 At-Satellite Reflectance. Remote Sens. Lett. 2014, 5, 423–431. [Google Scholar] [CrossRef]
  29. Zhang, X.; Schaaf, C.; Friedl, M.; Strahler, A.; Gao, F.; Hodges, J. MODIS Tasseled Cap Transformation and Its Utility. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Toronto, ON, Canada, 24–28 June 2002. [Google Scholar]
  30. Qiu, B.; Liu, Z.; Tang, Z.; Chen, C. Developing Indices of Temporal Dispersion and Continuity to Map Natural Vegetation. Ecol. Indic. 2016, 64, 335–342. [Google Scholar] [CrossRef]
  31. Crist, E.P.; Cicone, R.C. A Physically Based Transformation of Thematic Mapper Data-The TM Tasseled Cap. IEEE Trans. Geosci. Remote Sens. 1984, 3, 256–263. [Google Scholar] [CrossRef]
  32. Kalisa, W.; Igbawua, T.; Henchiri, M.; Ali, S.; Zhang, S.; Bai, Y.; Zhang, J. Assessment of climate impact on vegetation dynamics over East Africa from 1982 to 2015. Sci. Rep. 2019, 9, 16865. [Google Scholar] [CrossRef] [PubMed]
  33. Measho, S.; Chen, B.; Pellikka, P.; Guo, L.; Zhang, H.; Cai, D.; Sun, S.; Kayiranga, A.; Sun, X.; Ge, M. Assessment of vegetation dynamics and ecosystem resilience in the context of climate change and drought in the horn of Africa. Remote Sens. 2021, 13, 1668. [Google Scholar] [CrossRef]
  34. Qarallah, B.; Othman, Y.A.; Al-Ajlouni, M.; Alheyari, H.A.; Qoqazeh, B.A. Assessment of Small-Extent Forest Fires in Semi-Arid Environment in Jordan Using Sentinel-2 and Landsat Sensors Data. Forests 2022, 14, 41. [Google Scholar] [CrossRef]
  35. Sawalhah, M.N.; Al-Kofahi, S.D.; Othman, Y.A.; Cibils, A.F. Assessing rangeland cover conversion in Jordan after the Arab spring using a remote sensing approach. J. Arid Environ. 2018, 157, 97–102. [Google Scholar] [CrossRef]
  36. Sydykov, M.A. Direct sowing of winter wheat on rainfed lands of southern Kazakhstan. In Proceedings of the International Congress, Global Climate Change and Biodiversity, Almaty, Kazakhstan, 5–6 November 2015; pp. 177–182. [Google Scholar]
  37. Karabkina, N. Operational Solution of Space Monitoring of Grain Crops in South Kazakhstan Based on Data of Sentinel-2, Landsat-8, PlanetScope. In Proceedings of the 16th All-Russia Open Conference on Current Problems in Remote Sensing of the Earth from Space, Moscow, Russia, 12–16 November 2018; pp. 321–328. [Google Scholar] [CrossRef]
  38. Adenova, D.; Sarsekova, D.; Absametov, M.; Murtazin, Y.; Sagin, J.; Trushel, L.; Miroshnichenko, O. The Study of Groundwater in the Zhambyl Region, Southern Kazakhstan, to Improve Sustainability. Sustainability 2024, 16, 4597. [Google Scholar] [CrossRef]
  39. Adenova, D.; Tazhiyev, S.; Sagin, J.; Absametov, M.; Murtazin, Y.; Trushel, L.; Miroshnichenko, O.; Zaryab, A. Groundwater Quality and Potential Health Risk in Zhambyl Region, Kazakhstan. Water 2023, 15, 482. [Google Scholar] [CrossRef]
  40. Kyrgyzbay, K.; Kakimzhanov, Y.; Sagin, J. Climate data verification for assessing climate change in Almaty region of the Republic of Kazakhstan. Clim. Serv. 2023, 32, 100423. [Google Scholar] [CrossRef]
Figure 1. Boundary of the Tyulkubas district, Turkestan region.
Figure 1. Boundary of the Tyulkubas district, Turkestan region.
Sustainability 16 07548 g001
Figure 2. Sentinel-2 coverage scenes for the territory in the Tyulkubas district.
Figure 2. Sentinel-2 coverage scenes for the territory in the Tyulkubas district.
Sustainability 16 07548 g002
Figure 3. Field survey routes and test farms in the Tyulkubas District of the Turkestan Region during the period of 2–7 June 2022.
Figure 3. Field survey routes and test farms in the Tyulkubas District of the Turkestan Region during the period of 2–7 June 2022.
Sustainability 16 07548 g003
Figure 4. Visual differences between agricultural crops in the foothill zone of the Turkestan region during a route survey in April 2018.
Figure 4. Visual differences between agricultural crops in the foothill zone of the Turkestan region during a route survey in April 2018.
Sustainability 16 07548 g004
Figure 5. Photographic survey of surveyed fields in the foothill zone of the Turkestan region in April 2018: (a) winter wheat; (b) perennial alfalfa; (c) spring barley; (d) winter wheat (beginning of tillering); (e) safflower seedlings; (f) safflower cultivation and sowing.
Figure 5. Photographic survey of surveyed fields in the foothill zone of the Turkestan region in April 2018: (a) winter wheat; (b) perennial alfalfa; (c) spring barley; (d) winter wheat (beginning of tillering); (e) safflower seedlings; (f) safflower cultivation and sowing.
Sustainability 16 07548 g005
Figure 6. Graphical representation of pixel values of the Plowed Land Index for various types of objects: plowed fields 0–5, bare soil 5–10, fires 10–15, dense vegetation 15–20 and sparse vegetation 20–25.
Figure 6. Graphical representation of pixel values of the Plowed Land Index for various types of objects: plowed fields 0–5, bare soil 5–10, fires 10–15, dense vegetation 15–20 and sparse vegetation 20–25.
Sustainability 16 07548 g006
Figure 7. Comparison of results for identification of autumn plowing in the fields of the Tyulkubas district using the Plowed Land Index and NDVI data for 11 November 2021.
Figure 7. Comparison of results for identification of autumn plowing in the fields of the Tyulkubas district using the Plowed Land Index and NDVI data for 11 November 2021.
Sustainability 16 07548 g007
Figure 8. Boundaries of plowed fields based on the Plowed Land Index: (a) for 2 November 2021; (b) for 1 May 2022.
Figure 8. Boundaries of plowed fields based on the Plowed Land Index: (a) for 2 November 2021; (b) for 1 May 2022.
Sustainability 16 07548 g008
Figure 9. Methodology for recognizing the final mask of winter wheat crop sowing in South Kazakhstan.
Figure 9. Methodology for recognizing the final mask of winter wheat crop sowing in South Kazakhstan.
Sustainability 16 07548 g009
Figure 10. Visualization of changes in the field array of the Tyulkubas district on Sentinel-2 images: (a) field processing for the autumn period (September–November 2021); (b) field processing for the spring period (April–May 2022).
Figure 10. Visualization of changes in the field array of the Tyulkubas district on Sentinel-2 images: (a) field processing for the autumn period (September–November 2021); (b) field processing for the spring period (April–May 2022).
Sustainability 16 07548 g010
Figure 11. Autumn plowing mask of fields in 2021 over the territory of the Tyulkubas district, Turkistan region, based on the Plowed Land Index, Sentinel-2 Images.
Figure 11. Autumn plowing mask of fields in 2021 over the territory of the Tyulkubas district, Turkistan region, based on the Plowed Land Index, Sentinel-2 Images.
Sustainability 16 07548 g011
Figure 12. Confirmation of winter wheat crop vegetation based on the NDVI index and ground-based data on the Sentinel-2 satellite image for 24 May 2022.
Figure 12. Confirmation of winter wheat crop vegetation based on the NDVI index and ground-based data on the Sentinel-2 satellite image for 24 May 2022.
Sustainability 16 07548 g012
Figure 13. Map of autumn and spring plowing fields in the Tyulkubas district of the Turkestan region based on the Plowed Land Index in 2021–2022.
Figure 13. Map of autumn and spring plowing fields in the Tyulkubas district of the Turkestan region based on the Plowed Land Index in 2021–2022.
Sustainability 16 07548 g013
Figure 14. Satellite map depicting the final locations of winter grain crop sowing areas in the Tyulkubas district of the Turkestan region for the years 2021–2022.
Figure 14. Satellite map depicting the final locations of winter grain crop sowing areas in the Tyulkubas district of the Turkestan region for the years 2021–2022.
Sustainability 16 07548 g014
Figure 15. Plowing rates (%) in the Tyulkubas district of the Turkestan region in autumn 2021 based on satellite imagery.
Figure 15. Plowing rates (%) in the Tyulkubas district of the Turkestan region in autumn 2021 based on satellite imagery.
Sustainability 16 07548 g015
Figure 16. Fragment of satellite field with field verification of winter grain crops with ground truth data for “Kyzylzhar” farm in 2022.
Figure 16. Fragment of satellite field with field verification of winter grain crops with ground truth data for “Kyzylzhar” farm in 2022.
Sustainability 16 07548 g016
Figure 17. Fragment of satellite processed data output verification with the results of ground field surveys in the Tyulkubas district in 2022.
Figure 17. Fragment of satellite processed data output verification with the results of ground field surveys in the Tyulkubas district in 2022.
Sustainability 16 07548 g017
Table 1. Sentinel-2 satellite data used for Tyulkubas District, Turkestan region, for 2021–2022.
Table 1. Sentinel-2 satellite data used for Tyulkubas District, Turkestan region, for 2021–2022.
DateThe Sentinel-2 Satellite Data Used in 2021–2022
WN SceneXN SceneXM Scene
13.9–18.9T42TWN_20210913T42TXN_20210913T42TXM_20210918
26.9–28.9.T42TWN_20210926T42TXN_20210928
11–13.10T42TWN_20211011T42TXN_20211013
16–18.10T42TWN_20211016T42TXN_20211018T42TXM_20211018
31.10T42TWN_20211031
02.11T42TWN_20211102T42TXN_20211102
05.11–07.11T42TWN_20211105T42TXN_20211107T42TXM_20211107
27.11T42TWN_20211127T42TXN_20211217
04.04–06.04T42TWN_20220404T42TXN_20220406
09.04T42TWN_20220409
14.04T42TWN_20220414
21.04T42TWN_20220421
26.04T42TWN_20220426T42TXN_20220426T42TXM_20220426
01.05T42TWN_20220501T42TXN_20220501T42TXM_20220501
16.05T42TWN_20220516T42TXN_20220516
24.05–26.05T42TWN_20220524T42TXN_20220526
31.05–08.06T42TWN_20220608T42TXN_20220531
Table 2. Area of autumn plowing of fields in the Tyulkubas district in 2021 based on satellite imagery.
Table 2. Area of autumn plowing of fields in the Tyulkubas district in 2021 based on satellite imagery.
Observation PeriodThe Area of Autumn Tillage of Fields in 2021 (ha):The Total Area across 3 Scenes (ha):The Total Area across 3 Scenes (%):
TWN SceneTXN SceneTXM Scene
05.09–13.09.20216159.52272.5 6432.0229
14.09–26.09.20212768298.7 3066.714
27.09–13.10.20212374260.1 2634.112
14.10–18.10.20211821.295.2 1916.49
19.10–31.10.20212475.6 2475.611
01.11–10.11.20211804.5775.31571.434151.2319
11.11–27.11.2021923.9 923.94
01.12–17.12.2021 508.3 508.32
Total for the season (ha)18,326.7272.51571.4322,108.3100
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Arystanov, A.; Karabkina, N.; Sagin, J.; Nurguzhin, M.; King, R.; Bekseitova, R. Use of Indices Applied to Remote Sensing for Establishing Winter–Spring Cropping Areas in the Republic of Kazakhstan. Sustainability 2024, 16, 7548. https://doi.org/10.3390/su16177548

AMA Style

Arystanov A, Karabkina N, Sagin J, Nurguzhin M, King R, Bekseitova R. Use of Indices Applied to Remote Sensing for Establishing Winter–Spring Cropping Areas in the Republic of Kazakhstan. Sustainability. 2024; 16(17):7548. https://doi.org/10.3390/su16177548

Chicago/Turabian Style

Arystanov, Asset, Natalya Karabkina, Janay Sagin, Marat Nurguzhin, Rebecca King, and Roza Bekseitova. 2024. "Use of Indices Applied to Remote Sensing for Establishing Winter–Spring Cropping Areas in the Republic of Kazakhstan" Sustainability 16, no. 17: 7548. https://doi.org/10.3390/su16177548

APA Style

Arystanov, A., Karabkina, N., Sagin, J., Nurguzhin, M., King, R., & Bekseitova, R. (2024). Use of Indices Applied to Remote Sensing for Establishing Winter–Spring Cropping Areas in the Republic of Kazakhstan. Sustainability, 16(17), 7548. https://doi.org/10.3390/su16177548

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop