Spatial Distribution of Cropping Systems in South Asia Using Time-Series Satellite Data Enriched with Ground Data

: A cropping system practice is the sequential cultivation of crops in diﬀerent crop seasons of a year. Cropping system practices determine the land productivity and sustainability of agriculture in regions and, therefore, information on cropping systems of diﬀerent regions in the form of maps and statistics form critical inputs in crop planning for optimal use of resources. Although sat-ellite-based crop mapping is widely practiced, deriving cropping systems maps using satellites is less reported. Here, we developed moderate-resolution maps of the major cropping systems of South Asia for the year 2014–2015 using multi-temporal satellite data together with a spectral matching technique (SMT) developed with an extensive set of ﬁeld observation data supplemented with expert-identiﬁed crops in high-resolution satellite images. We identiﬁed and mapped 27 major cropping systems of South Asia at 250 m spatial resolution. The rice-wheat cropping system is the dominant system, followed by millet-wheat and soybean-wheat. The map showing the cropping system practices of regions opens up many use cases related to the agriculture performance of the regions. Comparison of such maps of diﬀerent time periods oﬀers insights on sensitive regions and analysis of such maps in conjunction with resources maps such as climate, soil, etc., enables optimization of resources vis-à-vis enhancing land productivity. Thus, the current study oﬀers new opportunities to revisit the cropping system practices and redesign the same to meet the challenges of food security and climate resilient agriculture.


Introduction
A cropping system indicates the crop types and sequences practiced in a region over the crop seasons of a year.Generally, a crop year consists of multiple crop seasons of roughly defined duration.The crops grown during these seasons comprise the cropping system of that region, characterizing the agricultural practices followed in that area as a whole rather than in an individual parcel of land.In India, there are three crop seasons: kharif (June to November), rabi (November to March) and zaid (March to June).Crops are grown in one, two or all three of these seasons, with variations from region to region.For instance, if rice is grown in the initial season followed by wheat in the second season and no crop in the third, then it is called a rice-wheat cropping system.The agroecological and agrometeorological conditions of a region have a bearing on the cropping system practiced there.
Cropping systems play a role in determining the land productivity of a region by facilitating optimal utilization of resources to maximize crop yields [1].However, resource-intensive cropping systems tend to impact soil health, ground and surface-water hydrology, environmental sustainability, climate resilience of agriculture and income security of farmers.Sustainable agriculture is closely linked to cropping system practices [2,3].Therefore, monitoring cropping systems is of critical importance to sustainable agriculture and food production.In particular, cropping system maps can provide important information on crop types, the number of crops grown and the management practices followed in irrigated or rainfed systems during a crop year at a particular location.Given the growing population and rapid climate change, there is increasingly a need to monitor and reduce emissions in the quest for sustainable agriculture development.
Cereals and legumes, especially short-duration crops, are crucial for food and nutrition security in Asia and Africa.They have the potential to improve livelihoods by enhancing incomes and providing stable employment.However, increasing cropping intensity through be er irrigation technology, fertilizer use and mechanization has led to growing pressure on croplands and other environmental concerns.The need to allocate water to other competing sectors (urban and industrial development) and reduced water availability in a highly varying climate inhibit the expansion of irrigated cropland area [4,5].Additionally, the adverse impacts of accelerated irrigation development and fertilizer use have been manifest in increasing salinity in croplands, leading to loss of soil fertility.
Given that context, mapping cropping systems can provide a crucial input into decision-making to improve cropping systems, explore opportunities to increase production and raise farmers' incomes under conditions where there can be no further expansion of agricultural area.Additionally, accurate crop information from the lowest administrative level is necessary to estimate country-level crop-wise acreage [6][7][8], which is a basic input that feeds into agricultural policies.Since dryland cropping systems are highly mutable and are influenced by biophysical, social and economic factors [9][10][11][12], accurate estimation of cropping pa erns is essential for planning.In such a scenario, geospatial products can provide input parameters for food security studies, and also for developing seasonal cropping pa erns of cereals followed by short-duration legumes [13,14].
There is a need to revisit the cropping system practices of different regions in view of the growing importance of sustainable agriculture and food systems.With climate change leading to increased exposure of agriculture to various risks, and fast-changing food habits altering the demand and supply scenario of food, there is increasing recognition of the importance of carbon-neutral economies and also reassessment of the global trade practices.Mapping and inventory taking of croplands following different cropping systems are basic information products that can aid in the development of more suitable cropping systems in different regions.Spatial information on current cropping systems integrated with layers of information on resources such as soil, weather and water will enable a holistic assessment of the sustainability of current agricultural systems.
There have been numerous studies on mapping cropping pa erns using satellite data, but they are limited to the subdistrict and plot levels [15][16][17] and single or major crops at large scale such as rice, wheat and sugarcane [18][19][20][21][22]. Currently, the various data sources available for use include regional and subnational statistical data on cultivated areas, coarse/medium-scale land-use/land-cover (LULC) maps and paddy maps from the 1980s to the 2010s [23][24][25], maps of rice areas using medium spatial resolution data [19,[26][27][28] and a combination of various other sources [29,30].
Modern remote sensing is characterized by frequent imaging of the earth's surface to produce large streams of spectral data in high to moderate resolutions and making available analysis-ready data and biophysical products in the public domain to map agriculture systems.These developments in remote sensing offer numerous opportunities to produce more accurate and customized information on crop cover and its changes.
Information on spatial distribution of different cropping systems is important for spatial analysis to identify areas appropriate for sustainable food production, for a be er understanding of where croplands are available to introduce/grow short-duration legumes across winter fallows, and to overcome current farmers' income uncertainties.Analysis of near-real-time satellite imagery can be a substitute for ground surveys to estimate the acreage of major crop types; moreover, it is relatively quick, cheap and independent [31][32][33][34][35].Some studies have used satellite data to map agricultural areas at various temporal and spatial resolutions [36][37][38][39][40][41]; moderate resolution imaging spectrometer (MODIS) normalized difference vegetation index (MODIS NDVI) time-series data to map croplands, cropping intensity at lower administrative levels and for river basins [42][43][44]; and land surface water index (LSWI) data to monitor natural vegetation-and soil-related wetness at a regional scale [26,[45][46][47][48].
However, to our knowledge, none of these studies have dealt with entire cropping systems, i.e., crops grown in a sequence over the entire year.This study maps major cropping systems across South Asia with the help of MODIS time-series data and spectral matching techniques (SMTs) to provide precise information on cropping sequences with the objective of enabling improvement of crop water productivity and proper natural resource management.

Study Area
South Asia is located between 5°38′ and 36°54′N latitudes and 61°05′ and 97°14′E longitudes, covering nearly 477 million ha.The region comprises Bangladesh, Bhutan, India, Nepal, Pakistan and Sri Lanka (Figure 1), and is divided, based on similar climatic conditions, into six agroecological zones (AEZs) [49] in which agriculture is strongly influenced by rainfall and seasonal winds.South Asia is home to 65% of the underprivileged people who live in rural and remote areas and depend on agriculture for their livelihood [50].People living in 70% of the countryside depend on coastal fisheries and land.South Asia has nine major river basins: Brahmaputra, Godavari, Indus, Ganges, Tapti, Krishna, Kaveri, Narmada and Mahanadi.Irrigation projects in this region serve a command area of nearly 133 million ha [51].Rice and wheat are the staple foods in the region; the former crop is grown throughout the year in some areas.The major cropping systems in South Asia vary by country (Table 1).India has diversified cropping systems, but rice-wheat is dominant in northern India and rice-rice in southern India.In Sri Lanka, rice-rice, rice-fallow and other crops-fallow are the dominant cropping systems.In Bangladesh, rice-rice, rice-pulses-rice, rice-fallow-rice and jutepulses-rice are popular; Bhutan favors rice-fallow and mixed crops.The dominant cropping systems in Nepal are rice-wheat, maize-wheat and maize-rice, whereas in Pakistan they are rice-wheat, rice-rice, rice-pulses and co on.Thus, rice-rice and rice-wheat are the common major cropping systems across the region.

Ground Reference Data
Ground survey information with geographical coordinates (Table 2) was collected broadly throughout South Asia for the crop year 2014-2015 for identification of classes of cropping systems together with an accuracy assessment exercise (Figure 1).Ground data collection was conducted using two distinct approaches.For the training data, homogeneous patches were selected, and comprehensive crop-related information was collected to generate ideal spectra signatures.In contrast, the validation data were gathered from various random locations, focusing primarily on land use, land cover, crop name rather than detailed information.A total of 2203 field samples were collected, of which 1416 points were used for training to develop an ideal spectral data bank (ISDB) and 883 for class validation (Table 2).Land-use land-cover information was collected for each sample (at least 250 m × 250 m) along with data on crop type, cropping intensity and irrigation water source.During collection of ground data, the area around the sample was categorized into one of three classes: small (≤10 ha); medium (10-15 ha); and large (≥15 ha) to enable acreage estimation.In some areas, farmers were interviewed to glean information on sowing dates, irrigation practices and cropping pa erns to aid class identification.For areas that could not be visited due to lack of proper roads, information was obtained from secondary sources, including local agriculture experts and records.Class names were allo ed in the field using a labeling protocol [31,52].Data limitations in some areas led to classes being identified on the basis of prior experience with crop signatures.
Data on India's agriculture were obtained from the official websites of the Department of Agriculture of the Government of India; other countries' data were obtained from the respective national agricultural/statistics departments.
As per Indian national statistics, the dominant rainfed cropping systems are co onfallow (15.3 M ha) followed by rice-fallow (12.6 M ha) and groundnut-fallow (9.7 M ha).The major irrigated cropping systems are rice-wheat (14.8 M ha) followed by millet-wheat (10.1 M ha).The other major cropping systems are pulses-wheat (9.2 M ha) and soybeanwheat (8.3 M ha) (Table 2).

Methodology
The workflow of this study consisted of three major steps: (1) preparation of satellite data including mosaicking and stacking; (2) running the machine learning algorithm, i.e., unsupervised classification (K-means) on Google Earth Engine (GEE) and identification of classes and labeling with the help of SMTs; and (3) accuracy assessment using validation data and national statistics (Figure 2).

Preparation of Satellite Data and Running Machine Learning algorithm
In GEE, the readily available MOD13Q1 V6.1 product provides preprocessed NDVI data for every 16 days on a per-pixel basis, masked for water, clouds, heavy aerosols and cloud shadows.Twelve MODIS tiles covering the study area for the study year June 2014-May 2015 were considered for our analysis.
The monthly maximum NDVI imagery was calculated using Equation (1).
The images were first stacked as monthly maximum-value composites, which means that the maximum NDVI values for every month were stacked for the study year 2014-15, making an analysis-ready data (ARD) cube of 12 months, one for each month, thus providing the NDVI variation in a crop year and overall coverage of the whole study area.MODIS data has temporal resolution of one day.Since this study considered monthly maximum images, there was less chance of missing pixel values; in an exceptional case, linear regression interpolation was done to fill the gap.
The final ARD cube was set into a K-means unsupervised machine learning algorithm with 100 clusters for every agroecological zone.The machine learning algorithm classifies the ARD cube into 100 clusters with 100 class spectral signatures.

Ideal Spectral Signatures (ISS)
Ideal spectral signatures for cropping systems were generated based on independent training data drawn from extensive field survey information using time-series data.The ground samples were grouped based on their distinctive characteristics and later classified into classes of unique cropping systems.Figure 3 shows the ISS of various classes.
The spectral signature curves explain cropping systems behavior over time.Regularly, NDVI values are low during initial crop development, high through peak growth and low again at harvesting stage, due to low to high reflectance value of near infrared region band during the crop development.In Figure 3a,b, double crop (rice-rice, ricewheat systems, etc.) classes exhibit two peaks in a single curve, signifying that the first peak is in the start season and the second peak is in the later season.Figure 3c,d shows only a single crowning, which means that there was only a single crop (e.g., groundnutfallow, sorghum-fallow, etc.) in a crop year.Depending on the NDVI values and ground data, we labeled classes as a combination of crop intensity and crop names.

Matching CSS with ISS to Group Classes Using SMTs
The initial 100 CSS signatures were grouped using quantitative spectral matching techniques (QSMTs), i.e., based on similar spectral signatures [53,54] (Figure 4a).The CSS (Figure 4c) signatures are then matched with the ISS (Figure 4b) and a match is determined by QSMTs with spectral correlation similarity (SCS) R-square value of 0.70 or higher are grouped and labeled as a class (Figure 4d).Classes with good visual matches both in terms of magnitude and shape of the spectral signatures are considered for class labeling.The preliminary labeling of classes was validated, mainly by using field data, Google Earth and other secondary sources.The process was repeated to classify and label all 100 classes, leading to a finalization of classes.Wherever there was limited ground data, highresolution data from Google Earth was used for reference.In the event of any ambiguity in the labeling of classes-e.g., some classes may not correlate with spectral signaturessuch classes were reclassified using the process above [52,54].

Accuracy Assessment
Accuracy assessment was then carried out with 883 independent field points that were not used in the classification process using an error matrix.The overall classification accuracy of cropping systems and users' and producers' accuracies were calculated for each cropping system.

Sub-Pixel Areas
Highly precise area estimates from the coarse resolution imagery are only possible through sub-pixel areas (SPAs) [13,52].This study used MOD13Q1 in which every pixel covers 250 m × 250 m, i.e., an area of about 6.25 ha.Full pixel areas (FPAs) are not an accurate representation of actual areas, as they may include areas falling under other LULC classes or mixed classes.On the other hand, sub-pixel areas (SPAs) involving actual area calculation, include the percentage of crop area.Therefore, areas based on SPA provide a near approximation of the cropped area.
Within the cropland class, there is a high possibility of discrepancies in the form of other LULC results.However, some classes have >50% area cropped and the proportions can vary widely.The percentage of various classes (water, cropland, etc.) is decided based on observation of ground data, photographs from ground data collection and by visualizing the field in Google Earth imagery.For example, in a rice-rice cropping system, a rice area sub pixel area (SPA), i.e., actual area is calculated based on the crop land area fraction (CAF), i.e., percent of the rice area in 250 m × 250 m area of random rice field observation samples in that class.This process was carried out for each cropping system obtained.To get the actual areas, the FPA was multiplied by the CAF [55].Since major cropping systems were mapped, the areas obtained are nearly equivalent to SPAs, with a few exceptions.Each class is evaluated for its actual area as follows: 4. Results

Cropping Systems
We find there are 27 major cropping systems in South Asia (Figure 5 and Table 3).Non-croplands are found in the extreme northern and western parts of South Asia and the Himalayas in the north, and are ~55% of the land cover of South Asia.The noncroplands consist mainly of deserts, mountains, glaciers and hills.The croplands are in the river valleys and vegetated plains.Around 45 percent of South Asia comprises croplands.There are 14 double cropped cropping systems and 11 single crops cropping systems, 1 mixed crop cropping system, and 1 other land use class, i.e., non-croplands in South Asia.Each cropping system has unique crop or crop-combinations.For example, there are several cropping systems with rice as one of the crops such as rice-rice double crop (i.e., kharif season rice followed by rabi season rice), or rice-fallow system in which there is only a single rice crop (i.e., kharif season rice followed by fallow winter season rice).Some other major classes were rice-wheat, rice-pulses, sorghum-fallow, millet-fallow and mixed crops.

Spatial Distribution of Each Cropping System
There is, in general, a dominating major crop of the rainy (kharif) season such as rice together with a secondary crop in the rabi and/or zaid summer season (not common).Below, we describe our findings of fifteen major cropping systems of South Asia in greater details (Figure 6): Rice-wheat: This is the dominant cropping system in the South Asia covering nearly 19.9 Mha stretching from Pakistan (along the Indus River basin) across into northern India (Indo-Gangetic plains) covering the states of Punjab in the west to Bihar in the east.In this system, rice crop is planted at the start of the rainy season (kharif) in June (analysis of spectral signature-Figure 3) and harvested at the end of summer monsoon season starting in August and lasting through September.As harvest progresses, NDVI reduces more gradually till it reaches a minimum and then around November wheat crop starts to emerge starting first with Bihar and ending last with those in Pakistan.Finally, wheat harvests start in March-April and they proceed much faster than rice harvests.
Rice-rice: Rice followed by a second rice crop can be observed in eastern India down to the south of India, Pakistan and in entire Bangladesh covering nearly 5.1 Mha of area.This cropping system is prevalent in irrigated lands and coastal areas.Unfortunately, the practice of mono-cropping in these regions has led to soil degradation and a deficiency in essential micronutrients.The first rice crop is in the kharif season, and sown during the monsoon months, typically between June and July, and is ready for harvest from September to October.After harvesting the first rice crop, the same fields are used for the second rice crop, which is planted during the dry winter months of the rabi season, usually between November and December.The second rice crop harvest starts from February to March as seen in their spectral signature.
Rice-pulses: Rice cultivation in rainy season and followed by pulses cultivation is observed in east to south India covering about 5.2 Mha of area.This cropping system is practiced in parts of Chha isgarh, Odisha, Tamil Nadu and Bihar states in India.It has led to high yields of both rice and pulses.However, the challenge lies in the irregular availability of irrigation in these regions.Rice cultivation takes place during the kharif season, where the sowing of rice seeds begins from June to July and harvesting between the months of September and October.Our spectral signature shows that subsequently, pulses are sown between November and December in this cropping system regions with the matured pulses harvested between February and March.Nearly 0.3 Mha is observed in Nepal and Bangladesh and no significant area of this cropping system is present in Pakistan, Bhutan and Sri Lanka.Rice-pulses-rice: This is the cropping system of three crops in a single crop year in which rice is followed by pulses and then by another rice crop.This cropping system is observed mostly in Bangladesh, but also in a few parts of Sri Lanka and India where irrigation coverage is excellent.The first rice crop is sown between June and July and is typically harvested from September to October.Subsequently, pulses are sown between November and December and then harvested from February to March.After the pulse harvest, preparations are made for the second rice crop to be sown in April to May, followed by the harvesting of the second rice crop from June to July.It is this la zaid season crop that is irrigated.
Soybean-wheat: The cultivation of soybean followed by wheat is seen majorly in Madhya Pradesh state of India covering nearly 6.6 Mha of area.This system is also found in Rajasthan, and certain regions of southern Maharashtra.Soybean cultivation in India primarily occurs during the kharif season, with our spectral signature showing the sowing of soybean seeds between June and July and harvest taking place from September to October.In the subsequent rabi season, after the soybean harvest, typically occurring in Oc-tober to November, fields are prepared for wheat cultivation.Wheat is sown from November to December, and the mature wheat crop is harvested between March and April.Approximately 0.1 Mha of this cropping system is observed in Pakistan, whereas there is no significant area in other countries.
Pulses-wheat: Unlike the rice-pulse-rice system, this system is concentrated in northern India and into neighboring Pakistan, spanning an area of approximately 4.5 million hectares (Table 1).It is prevalent in several states, especially northern Indian regions, including U ar Pradesh, Madhya Pradesh, Rajasthan, Haryana and Punjab states.The usual pa ern involves planting kharif pulses like chickpeas, pigeon peas and lentils during the monsoon season, followed by the rabi wheat crop.Pulses cultivation follows a specific timeline, with pulses being sown from June to August and harvested between October and November.Wheat is sown from November to December, and the mature wheat crop is ready for harvesting between March and April.
Maize-wheat: Maize followed by wheat cultivation is observed in western parts of India and Pakistan, covering nearly 6.1 Mha of area.In this cropping system, maize is the primary crop during the kharif season, while wheat takes center stage in the rabi season.This system is widely practiced in states like U ar Pradesh, Rajasthan, Madhya Pradesh, Bihar and Punjab.A significant portion of the maize-wheat system is rainfed, which means it heavily relies on rainfall, posing a major challenge due to its unpredictability.In the kharif season of India, maize sowing is carried out between June and July, and the maize harvest occurs during September to October.Subsequently, in the rabi season that follows, wheat is sown between November and December, and the mature wheat crop is usually ready for harvesting between March and April.Approximately 0.1 Mha of this cropping system is observed in Pakistan, whereas there is no significant area in other countries.
Millet-wheat: The cultivation of millet followed by wheat is observed in northern parts of India and also in Pakistan, covering 12.5 Mha of area.These are significant cropping systems in the western part of India, particularly in the semi-arid regions of Gujarat, Rajasthan and Haryana.Wheat is cultivated as a rabi crop, while Millet is grown during the kharif season.Millet cultivation is carried out between June and July, with millet ready for harvesting in September to October.Subsequently, in the following rabi season, wheat is sown from November to December, and the mature wheat crop is typically harvested between March and April.
Soybean-chickpea: This type of cropping system is observed in the central parts of India covering nearly 5.4 Mha of area, especially in the areas of Madhya Pradesh and some areas of U ar Pradesh states.During the kharif season in India, soybean sowing starts between June and July, with the soybean harvest typically taking place from September to October.In the subsequent rabi season, chickpeas are sown from November to December and are ready for harvesting between February and March.
Sesame-mustard: The cultivation of sesame followed by mustard is observed in U ar Pradesh state of India covering nearly 3.9 Mha of area.The sesame-mustard cropping system is adopted in various states across India, including Rajasthan, Gujarat, Madhya Pradesh, U ar Pradesh, Bihar and portions of Haryana and Punjab.In the kharif season, sesame cultivation starts with sowing between June and July, and the sesame crop is typically harvested from September to October.In the subsequent rabi season, Mustard is sown from November to December and is usually ready for harvesting between February and March.
Sugarcane: Sugarcane is a whole-year crop, i.e., cultivated throughout the crop year mostly in U ar Pradesh and Maharashtra states.Sugarcane is grown in several states across India, with U ar Pradesh, Maharashtra, Karnataka, Tamil Nadu, Andhra Pradesh and Bihar being the major sugarcane-producing states.Sugarcane is usually planted in the pre-monsoon or early monsoon season.In northern India, it is planted from March to June, while in southern India, planting can even extend into July.
Rice-fallows: Rice fallows are predominantly found in both rainfed and irrigated areas of the south Asia.The significant rainfed rice-fallows area found in India of about 11.4 Mha and mixed irrigated rice-fallows of nearly 13.3 Mha, whereas Bangladesh and Pakistan have nearly 0.5 Mha.
Pigeon pea-fallows: The pigeon pea-fallow cropping system is mainly practiced in the semi-arid and arid regions of India, including areas in Rajasthan, Gujarat, Maharashtra, Karnataka and Madhya Pradesh, covering about 6.7 Mha.Generally, depending upon the rainfall, the sowing of pigeon pea will be done and cultivation period ranges from 6 months to more than one year in some areas.
Co on-fallow: Co on-fallows areas are also significantly found in India, which is about 16.2 Mha of area.This cropping systems is predominantly in states of Karnataka, Maharashtra and some parts of Gujarat.Co on sowing happens in June-July, and harvesting takes place from October to January.
The remote sensing-based data on full-pixel crop areas of major cropping systems (Table 3) was further refined as the actual area statistics which were calculated using SPA by taking random samples of each class and identified the percentage of crop found in respective class pixels.The SPA was calculated given the possibility of non-crop pixels in crop pixels due to coarser resolution.
Sub-pixel areas in each cropping system were calculated.It was observed that on an average, 10 to 20 percent of each class contains other LULC.We calculated the sub-pixel areas for each class to get the actual area by eliminating 10-20 percent accordingly.

Country-Wise Cropping Systems
Country-wise cropping systems and their respective sub-pixel areas shows the dominant cropping sequence followed in each country (Table 4).In Pakistan, the dominant cropping systems are rice-wheat (4.3 M ha) followed by millet-wheat (3.4 M ha) and ricerice (1.4 M ha).In Nepal, the major cropping systems are rice-wheat (0.4 M ha) and ricepulses (0.36 M ha).In Bhutan, the cropping systems are rice-fallow (0.012 M ha) and pulses-rice-rice (0.007 M ha).In Bangladesh, the dominant cropping systems are pulsesrice-rice (3.8 M ha) followed by rice-rice (1.1 M ha), whereas in Sri Lanka, the dominant cropping system is pulses-rice-rice (0.49 M ha).In India, there are highly diversified cropping systems such as rice-wheat (15.1 M ha), rice-fallow (11.4 M ha), co on-fallow (16.2 M ha) and rice-fallow/mixed crops (13.3 M ha), followed by a significant shared area with different cropping systems.The dominance of rice-fallow cropping systems provides opportunities for intensification with short-duration pulse crops.Data on area by cropping system was extracted from a remote sensing-derived cropping systems map and compared with National Agriculture Statistics (NAS) for the year 2014-15 obtained from the Department of Agriculture.The analysis reveals significant discrepancies between the reported areas in NAS and those derived from the remote sensing map for various crop combinations.Crop combinations like pigeon pea-fallow (77.56%), maize-chickpea (54.95%) and pulses-fallow (67.82%) show substantial positive differences.Conversely, soybean-wheat (−13.02%),millets-fallow (−8.91%) and milletswheat (−3.44%) exhibit negative differences.Some combinations, such as rice-fallow (−0.99%) and sugarcane (6.46%), show minimal differences, indicating a varying level of alignment between reporting and mapping methods (Figure 7).The accuracy assessment was calculated for identified cropping systems with 883 independent ground points.Overall classification accuracy was observed to be 76.59 percent, with kappa coefficient of 0.75.User and producer accuracies of major cropping systems ranged from 70 per cent to 90 percent, whereas mixed cropping system showed less accuracy ranging from 50 percent to 70 percent.
The major rice-wheat cropping system showed a user accuracy of 89 percent and producer accuracy of 72 percent, and pulses-wheat systems revealed a user accuracy of 92 percent and producer accuracy of 100 percent.Likewise, each cropping system was validated using validation points, their individual class accuracies and kappa coefficient determined (Table 5).Due to domination of small-holding farms and mixed crops in major areas of south Asia, the accuracy of some of the minor cropping systems was lower.The mapping accuracy of such cropping systems can be increased by utilizing high resolution satellite data along with more extensive ground data collection to build be er models.

Discussion
While the classification of irrigated, homogenous crops such as rice and wheat were relatively easy throughout South Asia, identification of rainfed crops like groundnut and pulses was difficult.Further, cropping pa erns in most regions of south Asia vary widely, with a high possibility of mixed crops.In India and Bangladesh, rice is the major irrigated crop and cultivated in large homogenous extent, whereas in Pakistan, wheat is the major irrigated crop, comprising more than 25 percent of all crops.Rice and wheat cropping system areas could be identified and classified more correctly.The identification of other crops mainly depends upon extensive ground data, their ideal signatures and crop phenology.Western South Asia comprising Afghanistan and Pakistan has large tracts of uncultivated area, whereas eastern South Asia, i.e., northeastern India, contains hill areas and the north part of south Asia contains ice glaciers.Crop identification in hilly regions was affected by shadows and different cropping pa erns.Nearly 26 percent of South Asian cropland consists of double crops (i.e., rainy and winter season crops).Double crops in irrigated areas were identified using spectral curves and ground data.Only 45 percent of the total area in south Asia was under cultivation, with other LULC comprising mainly deserts, glaciers and hilly areas.Knowledge of regional cropping systems is essential to maximize land productivity and improve crop production towards the objective of sustainable agriculture.Of late, high-intensity cropping practices and overexploitation of resources have resulted in unsustainable agriculture.The 27 cropping systems identified in this study cover all major and dominant agricultural practices in the region.These novel information products can serve as baseline data for developing sustainable agriculture solutions.When integrated with thematic layers such as weather, soil, groundwater availability and land degradation in GIS environments, they can unlock opportunities to develop action-oriented solutions to improve the performance of agriculture systems.It is also possible to produce maps with finer granularity using Sentinel data of 10 m resolution once more ground level data are collected.Sustainable agriculture and food system solutions need such innovative information products, which are potentially the low-hanging benefits from modern remote sensing technology.
Policy Implications: Cropping system maps and statistics form basic inputs for formulating sustainable agriculture plans to optimize land productivity.The major challenges of agriculture in South Asia are food security, income security and climate resilience, which can be addressed only through systematic efforts to optimize crop pa erns and resource consumption in producing crops.The maps of cropping systems can now be superimposed on resources and climatology data layers to offer deeper insights to decision-makers on the efficiency of the current agricultural systems and to prioritize critical areas for immediate interventions.The design and implementation of such interventions needs a strong base of data and information on cropping systems and related analytics that our maps will now provide.These interventions will involve both management and structural measures at the grass root level, and may involve certain incentives and investments, but cannot proceed without understanding the entire cropping system of the local area.Specifically, the government of India's flagship projects on water use efficiency in agriculture, the National Food Security Mission (NFSM) and Agroforestry development mission will become more effective when cropping systems information is integrated within their design and implementation plans.The cropping system determines multiple ecosystem services, pest predation and water quality.For example, the development of agroforestry plans in a region should take into account the existing cropping system.Carbon sequestration and carbon emission potentials of croplands are directly related to the cropping systems followed in a region.Therefore, any effort towards achieving carbon neutral agroecosystems will fundamentally get guided by cropping system practices.
Policymakers across South Asia will find our cropping system maps to be quite useful, since they provide a thorough understanding of the geographical distribution of various crops within an area.These maps will help with all governmental planning of resource allocation, making investments in pest management, irrigation and crop-specific subsidies more effective.Furthermore, they are crucial in guaranteeing food security since they enable policymakers to evaluate the yield of vital food crops and create plans for agricultural output diversification.By highlighting regions with significant production potential for particular crops, cropping sequence type maps will also aid in market growth and strengthen local and regional economies.They also aid in the identification of vulnerable regions that require protection, the evaluation of the effects of farming on natural resources and the development of sustainable agricultural policy.All things considered, these maps are invaluable resources that support the development of evidence-based policies and well-informed decisions about a variety of aspects of land use and agriculture.

Conclusions
Knowledge of regional cropping systems is essential to maximize land productivity and improve crop production towards the objective of sustainable agriculture.Of late, high-intensity cropping practices and overexploitation of resources have resulted in unsustainable agriculture.
This study has developed a new geospatial application in the agriculture sector to address crop planning and sustainable agriculture practices.
Remote sensing imagery from MODIS 250 m was used to generate maps and statistics on different cropping systems practiced in South Asia.The 27 cropping systems identified in this study cover all major and dominant agricultural practices in the region.These novel information products can serve as baseline data for developing sustainable agriculture solutions.When integrated with thematic layers such as weather, soil, groundwater availability and land degradation in GIS environments, they can unlock numerous opportunities to develop action-oriented solutions to improve the performance of agriculture systems.It is also possible to produce maps with finer granularity using Sentinel data of 10 m resolution.Sustainable agriculture and food system solutions need such innovative information products, which are low-hanging benefits from modern remote sensing technology.

Figure 1 .
Figure 1.Spatial distribution of ground data collected for the crop year 2014-2015.

Figure 2 .
Figure 2. The steps to map crop dominance classification.

Figure 3 .
Figure 3. (a,b) Ideal Spectra Signatures of double crops (c,d) Ideal Spectra Signatures of single crop 3.2.2.Class Spectra Signatures (CSS) Class spectra signatures (Figure 4a) were generated for the 100 classes (clusters) from unsupervised classification (K-means) of the 16-day NDVI temporal data.The signatures were generated for every individual class by taking the NDVI mean of all classified pixels of each cluster for every month.

Figure 4 .
Figure 4. Matching CSS with ISS and grouping of similar classes.

Figure 5 .
Figure 5.The distribution of 27 major cropping systems across South Asia.

Figure 6 .
Figure 6.Maps of spatial distribution of individual major cropping systems.

Figure 7 .
Figure 7.Comparison of remote sensing-derived crop areas with that from national statistics.

Table 1 .
Country-wise major cropping systems in South Asia.

Table 2 .
Field samples used for training and validation, and national statistics on cropping systems in India.

Table 3 .
The area covered under major cropping systems of South Asia.

Table 4 .
Country-wise net crop area under different cropping systems.

Table 5 .
Accuracy assessment of the croplands product by error matrix method.