Across the world, it is witnessed by economists that agriculture plays a fundamental role in the economic achievements of a country. But in modern days, the share of agriculture is waning that ultimately led to a lower rate of saving and employment. Kirby et al. (2017) have defined that in Pakistan, forty-seven per cent of the population is food insecure, and food production mainly depends on traditional irrigation system from already stressed water resources. This largely is due to the higher elasticity of demand in non-agriculture products and services. Agriculture subsidiaries like crops, dairy, fisheries and animal husbandry are a major source of income in Pakistan (GOP, 2017). Agriculture is a relatively well-developed sector of Pakistan's economy since Pakistan's independence, and it is considered as a driving wheel of the economy. It accounts for nineteen per cent of GDP (GOP, 2017) and directly supports more than half of the total population. Pakistan is heavily dependent on its major crops. The main crops of Pakistan include cotton, sugarcane, rice and wheat. Pakistan has one of the world's largest irrigation systems to support the production of its agricultural crops. There are mainly two seasons for the production of crops, i.e. May to November and November to April. The crops of cotton, sugarcane and rice are sown in May and harvested in November, whereas the crop of wheat is sown in November and harvested in April. There is an extensive production gap between the actual and required output of production due to the low level of knowledge of peasants about modern and efficient agriculture systems (Aslam, 2016). Declining water resources (Faruqui, 2004), un-equal land fertility and less education about pest and plant diseases control are other reasons for the said gap (Amanullah et al., 2009; Yang et al., 2005; Kamal & Moghal, 1968). Unwisely use of fertilizers, and modern machinery in the agricultural sector is yet another reason for this mismatch that underpins the urgent need for efficient resource utilization (Rehman et al., 2016).

Punjab is a major province of Pakistan, largest by population and second-largest by area. It has 36 districts, with Lahore being its capital. Punjab has a population of 99.9 million people, i.e. 57 % of the total population of Pakistan. This province has 57% cultivated and 69% of the total crop sown area in Pakistan. Its contribution to the agriculture of the economy of Pakistan is 83% of cotton, 80% of wheat, 97% of rice, 63% of sugar cane and 51% of maize. Agriculture is the main source of income for almost 77% of people in Punjab. Agriculture resources mainly include land and water to produced crops of wheat, cotton, rice and sugar cane. These crops account for more than 75% of the value of total crop production (Finance Division; Govt. of Pakistan, 2015). The crops under study account for 64 % of the annual crop area and about 90 % of major agricultural crops in Pakistan (Finance Division; Govt. of Pakistan, 2015).

Wheat is a major cereal crop in Pakistan, particularly in Punjab. It occupies 37% of crop area and 45% of other agricultural resources (Amanullah et al., 2019). It involves investment in canal irrigation and market development. Wheat has a total share of 14% in agriculture. Punjab followed the green revolution movement in the production of wheat since 1965, due to which wheat per hectare production has increased, and rural society has transformed (Hussain et al., 2002). Despite all revolutionary measures, Punjab Government still needs to reform wheat production in various varieties (Smale et al., 2002). Rice is the second major crop of Pakistan. It is being exported worldwide from Pakistan, which accounts for approximately 10% of the world's trade annually. Only basmati rice has a share of 25% of Pakistan exports. Rice export is the second major source of income in Pakistan that fulfils at least 60% of food need. Cotton is another major cash crop in Pakistan. Pakistan efficiently plays the role of raw cotton producer (Khan & Khaliq, 2010). Cotton exports account for 46% of total exports, and it engages 35% of labour. Production has increased by multiple times from 1952 to 2004 (Ali et al., 2013). Although Punjab and Sindh are the major producers of cotton in Pakistan, Punjab province is the leading producer with a larger cultivated area. Out of total production, more than 80 % comes from Punjab, whereas only 18 % comes from Sindh. Sugar cane is also one of the most important cash crops of Pakistan and the driver of the sugar industry in Pakistan. It accounts for 0.7% of GDP and 3.7% of agriculture value (Godfray et al., 2010), and in rural areas, sugarcane cultivation is a major sign of socioeconomic development because it provides high income and employment opportunities.

The above representation on agriculture reveals that a lot of research has already been surpassed on evaluation of overall production of agricultural crops and their contribution to GDP. However, little research work has been conducted about performance evaluation of administrative units or constituencies. It has become imperative and important to assess the contribution of administrative units (particularly at the level of districts) that to what extent they are efficient to contribute to GDP. It is a crucial decision to evaluate geo-politics using a multitude of indicators instead of anyone measure (Latruffe et al., 2016). Therefore, this study aims to evaluate the performance of the districts of Punjab. It focuses on providing insight for formulating policies regarding sustainable food security, making the agriculture process cost-effective and wisely use of fertilizers and modern machinery in the agricultural sector (Amanullah et al., 2019). The rest of the paper is arranged into the literature review, methodology, data analysis and concluding remarks.

There is an influx of literature on agricultural performance and sustainability at countries' level and/or provinces' level, but there is scarcity at divisional and district levels, particularly with reference to Pakistan. Therefore, a fairly rigorous literature review has been conducted that provides plausible ground for this comparative study of agricultural performance. The relevant studies are therefore iterated. Anastasia and Wim (2017) conducted a research study in Russia using secondary data from Federal State Statistics Service and evaluated different municipalities. It emphasized that there are only a few studies evaluating the agricultural performance of different territories like municipalities despite the fact that it is important to study them in comparison using composite indicator indices. Gezici and Hewings (2004) carried a research study in Turkey using the secondary data of the Istanbul Chamber of Industry and data of the State Statistics Institute. It compared different regions of Turkey from an economic perspective. Trukhachev et al. (2014) proclaimed that it is crucial for the economies to evaluate their constituents on an objective basis. This study was conducted on the comparative agricultural performance of different federal district regions of Russia. It used archival statistical secondary data. This study used mixed methods for data analyses and accentuated the importance of comparison of districts on the basis of agriculture production. Talukder et al. (2017) carried a research study on indicators of agricultural performance. They asserted that there are more than fifty indicators of agricultural performance that can be categorized into productivity, stability, efficiency, durability, compatibility and equity. There is also a large array of techniques for measuring the integrated performance of agricultural units. It also argued that it is important to determine which indicator to include, what weight to be given to each indicator and what technique of analysis to be employed in order to fairly represent the comparative regional positions. This study used composite normalization, weighting and aggregation techniques on primary and secondary data and recommended rather different and relatively newer techniques of analyses to obtain fairly generalizable results.

Literature is rich in different aspects of agriculture; Lewis (2017) stated that agriculture always played a key role in economic growth and eliminating food insecurity. In various economies, household's income is primarily based on agriculture and their ability to exploit natural resources efficiently (Hussain et al., 2015). Keesstra et al. (2016) evaluated that the world population was growing fast in the last few decades, and food resources were shrinking that exposed the world to a challenge that how to properly utilize available resources, particularly soil, water and natural fertilizers. Araya et al. (2021) found that to counter these emerging challenges, they need to use available resources wisely. Due to probable food insecurity, it is perilously significant and resultantly irresistible for developing countries (Lewis, 2017). Borrelli et al. (2015) inclined that ever-mounting population strain on planet and decreasing water resources leading to filth condition of these resources. Fahad et al. (2017) have employed analysis of covariance to evaluate the genetic diversity in wheat crop production with respect to enhancement in production. Results found positive relation among studied variables. This results in the shape of reduced productive capacity, soil degradation, soil erosion and teeming food insecurity. Tanveer et al. (2017) asserted that soil degradation and depletion of natural resources is one of the most high-flying challenges for the modern agriculture sector. Byerlee et al. (2009) stated that the low share of agriculture in economic development and in GDP has a foreseeable upshot in economic progress. According to the agricultural economies witnessed a fall in food production due to low-income elasticity of demand of agriculture production. Maraşlıoğlu (2017) bolstered that scarcity of water resources is the major factor for limiting agriculture development, and it drastically affects food production. Under changing climate, conserving soil and water resources are critical and important to sustain the agriculture sector. It is important to highlight that scarcity of water will alternatively lead to the overuse of underground water, consuming a lot of energy to pumping it out, whereas Pakistan is already having a problem of an energy shortage, specifically in rural areas. Ladha et al. (2003) buttressed that rice is one of the most important food crops in Pakistan. Jalota et al. (2018) concluded that in central Punjab of India, more groundwater withdrawal to irrigate rice, and wheat crop lead to the reduced water table. By improving the irrigation system, shifting planting date, and laser levelling reduced the water losses. Gangwar and Prasad (2005) claimed that in the last few years, both systems suffered from belated planting of subsequent crops resulting in loss of crop yield. Losses can be reduced with a proper cropping system. Khan and Khaliq (2005) found that wheat crop sown before cotton significantly gives a high yield then it is harvested after cotton. It is attributed to high plant density, more fertile tillers and grain spike. Watoo and Mugera (2014) stressed that tube well users and water buyers could enhance their production by 19% and 28% with the use of quality inputs as seeds or fertilizers. Cotton production in Pakistan is directly dependent on fertilizers price, competing for crop price and fertility of the land. Ali et al. (2017) argued that in addition to all these factors, there is a need to educate farmers, enhance credit facilities and improve transportation. Ali et al. (2013) suggested that there is a severe need to educate growers about plant protection measures and the use of fertilizers to get higher production of cotton.

Copious literature has been surpassed on evaluation of production of (output) wheat, cotton, rice and sugar cane along with inputs, e.g. fertilizers, tube wells, harvesters, mounters, reapers, pumps etc.

However, the researchers could not find any study that has evaluated major crops along with their inputs with reference to administrative units of Punjab. Deng (1982) introduced revolutionary grey system theory with its analytical methodology, grey relational analysis. This approach is specifically used to cross analyze a multitude of variables and to evaluate the best alternatives in multi-criteria decision making. The main goal of the theory is to calculate the grey relational coefficient and obtain a grey relational grade among variables in order to evaluate the vague research problems (Ju-Long, 1982). This theory was first introduced for science discipline only, but with the passage of time, it came out as a bridge between natural and management sciences (Liu & Lin, 2010). It is now known as an inter-disciplinary approach. (Chen & Ting, 2002) argued that a grey system exists between black and white systems where white means information needed is exactly available and black system means information that is needed but not available. Therefore, it is "grey". Under a grey system, a connection is established to discover values that are not clear or poorly stated (Hussain et al., 2002). It has five components in grey system theory, including grey predictions, grey decisions, grey programming, grey control and grey relational analysis. Grey Relational Analysis (GRA) is a popular method used with discrete data set. It is a mathematical technique useful in decision making of multi attribute cases. The advantage of GRA is that it is based on real data, simple mathematical calculation and one of the highest reliable methods of decision making in the business world (Chen & Ting, 2002). It quantitatively compares the variable in a vibrant way and establishes the relationship among variables based on similarity and inconsistency. In the past, GRA was used to evaluate: behaviour of energy consumption (Yu et al., 2011), decision making in the fuzzy system (Alcantud, 2018), supplier selection (Yang & Chen, 2006), portfolio selection (Bijarniya et al., 2020), software project performance evaluation (Shepperd & MacDonell, 2012) , safety assessment of e-commerce system (Radziszewska, 2018), wealth management (Wu et al., 2010) and medical data analysis (Xuerui & Yuguang, 2004). Ertugrul and Prof (2016) stated that grey relational grade is used to assess the degree of relationship among factors. Kao and Hocheng (2003) asserted that grey relational analysis helps to analyze the uncertainties and disparities effectively. In contemporary research, authors used various Multi-Criteria Decision-Making Techniques (MCDMT) to evaluate the financial and economic performance of various firms and countries, such as a technique for order of preference by similarity to an ideal solution, analytical hierarchy process, and data envelopment analysis etc. depending upon the nature of the problem under study it was found appropriate to use GRA (one of the techniques from MCDMT). Grey relational theory is used to employ GRA in order to rank the performance of districts of Punjab with respect to agricultural indicators.

This study follows the positivist research philosophy with a deductive approach. It is a cross-sectional study based on secondary data taken from Punjab Development Statistics 2016. It follows the mono method mathematical approach, i.e. GRA. To evaluate the agricultural performance of districts of Punjab, grey relational system theory is used. Dogan (2013) revealed that this system is frequently employed to use incomplete and impure information for analyzing the relations among a multitude of variables. Other statistical techniques like regression analysis could also be used to measure the performance of the agriculture sector, but these have many limitations as a large amount of data is needed, which generates an unsatisfied level of results (Uckun et al., 2012). The Grey relation system is an appendage to all those typical statistical techniques. GRA progresses stepwise, i.e. seven steps in accordance with Kuo et al. (2008); Wu (2002); Hamzacebi et al. (2011); Niazi et al. (2021), and Tayyar et al. (2014).

In the first step, data of all districts were obtained from Punjab Development statistics 2016, followed by the creation of reference series and comparable sequence in the next step. Reference series was formed by identifying the best alternative from the normalizing matrix. In the next (third) step, a normalized matrix was created, i.e. a single data set for healthy comparing and grey relation generation (Tsai et al., 2003). In this study, all four major crops and total area of sown was taken as larger the better degree because the more crops production we have, the better it is. However, other indicators like electric pumps, diesel pumps, tube wells, tractors, threshers, harvester, mounter were normalized on ideal value criteria as these variables are considered to use an ideal amount, not less or more in quantity. Normalized values of fertilizers were obtained by using "smaller the better" criteria as fertilizers are the cost of production, and in economic production, it is needed to use the least available cost and produces more output from it. Subsequently, absolute values were obtained in step 4 by calculating the deviation sequence. A deviation sequence is created to show how much a value can differ from its desirable value of series. Absolute values were created with the difference in the reference sequence and comparable sequence. Step 5 was completed by establishing a co-efficient matrix of grey relation system. Grey relational grade was calculated in the ensuing step, which is equal to all weighted sum of values, and the value with the highest grade is selected as the best alternative choice. In the final step, grey relational grades of all districts were arranged in ascending order to get the highest and lowest grades, respectively.

Following steps of GRA were used to access the best performer among different districts of Punjab:

Step 1: Created a data set and established a decision matrix of data set using the following formula:

x_i(k)=(x₁(1)x₁(2)⋯x₁(m)⋮⋱⋮x_n(1)x_n(2)⋯x_n(m)) Eq. (1)

Where i=1, 2, 3............n k = 1, 2, 3 ..........m

The results showed that district Bahawalpur and Bahawalnagar have the highest grey relational grades, hence depict the best agricultural performance in Punjab, whereas district Mianwali has the lowest grey relational grade and accordingly least performance Table 6. The justification of results can be explored from inspecting the original data comparison of criterion variable of each district one by one can support the rationale of the results. Since the study uses equal weights for all the variables while during the application of GRA and Mianwali has clearly lesser use of the inputs and produces a lesser quantity of crops, therefore, it seems quite logical to attain a lower GRA relational grade as compared to rival districts that have high inputs and outputs. However, if the unequal weights are assigned to criterion variables, there may vary in results. Other reasons for this difference can be like ignoring the production of pulses, fruits, tobacco and vegetables that may again have an impact on the ranking. There is also an array of techniques of assessments like data envelopment analysis etc., that can also produce different results. The results of the study are comparable with some contemporary studies. Cao et al. (2017) have employed the water footprints (WF) technique to measure the use of water in rice crop and its impact on quality and water volume. Findings concluded that the quality of crop and usage of water volume differs accordingly to the availability of water drainage in districts. Kirby et al. (2017) examined the historical trend of crop production, water availability and food insecurity in Pakistan and projected them forward to 2050. Results showed that sown area, water usage and fodder for crops is increasing gradually. It also predicted that if this trend continues for a long time, it will lead to a doubling of groundwater usage. It further suggested that to secure already diminishing water resources; we need more dams, irrigation infrastructure, increasing crop yields, change in crop mix and import more or export less food. They recommended using more grounded water in the short to medium term. Rehman et al. (2017) examined the relation of rice output with water availability, production area and its effect on GDP. Regression analysis was employed, which revealed that area production and rice output have a positive, while water resources have a negative relation with the GDP of Pakistan. Ali et al. (2019) said that in India and Pakistan, agriculture is a major sector of the economy. They focus on improving the sustainability of the agricultural sector and understand that it is important to secure food resources for the local population as well as for other economies that rely on food import from these countries. Research demonstrated that they need to conserve water resources that are already declining. There are plenty of other studies investigating a multitude of perspectives on agriculture performance, importance and issues. However, contemporary research could not focus the output and input variables compositely like crops, threshers, mounters and reapers' contribution to agriculture produce (specifically at district level). Since it is a unique type of study, therefore, its results could not be compared with previous research. However, based on the statistics given in archival documents of the government of Pakistan (i.e. periodical publications of Punjab development statistics and the federal bureau of statistics), the results of the study are consistent but different in their nature and magnitude. Because it has taken not the only output of the agriculture sector of Punjab but also considered the input variables and evaluated their efficiencies to how much they contributed to agricultural production with given resources.

The agriculture sector is the driving wheel of Pakistan's economy because Pakistan is heavily dependent on major crops. It accounts for approximately 19 % of GDP (GOP, 2017) and directly supports more than half of the total population. The province of Punjab is the most important agricultural constituent of Pakistan that is renowned for cotton, rice and wheat crops. This study was thus aimed to evaluate the districts of Punjab with respect to the use of resources and agricultural performance. For the evaluation, a wide range of multi-criteria decision-making techniques were considered but, since the purpose of the study was to evaluate and ultimately hierarchical districts on the basis of incomplete information of agricultural indicators, therefore, Grey Relational Analysis (GRA) was found to be appropriate. The analysis showed that district Bahawalpur and Bahawalnagar have the highest grey relational grades, hence depict the best agricultural performance in Punjab, whereas district Mianwali has the lowest grey relational grade and accordingly least performance. The results are plausible because they coincide with quantities of inputs and outputs. The districts having high quantities in all variables secured high grey relational grade and vice versa. It apparently seems simple because larger acceptable seemingly dominate; however, it is a bit complex due to some smaller acceptable. Other reasons include equal weights to variables in GRA and ignoring some of the crops as aforementioned. It is a seminal study that uses a unique methodology that has integrated thirteen different indicators of agricultural development in one mathematical model and assigned a distinctive composite grade to every district. It provides insights to policymakers to take corrective measures and/or adjust agriculture development policies. The results of the study can work as ingredients to design policies for sustainable food security. The study encompassed thirteen variables, including four major crops and nine non-crop indicators but, could not include fruits, vegetables, tobacco and pulses, which may be incorporated in future studies for sanctifying the results.

Bahawalpur = γ-1, Bahawalnagar = γ-2, Rahimyarkhan = γ-3, Dera gazi khan = γ-4 , Layyah = γ-5, Muzafarabad = γ-6, Rajunpur = γ-7, Faislabad = γ-8, Chiniot = γ-9, Jhang = γ-10, Toba take singh = γ-11, Gujranwala = γ-12, Gujrat = γ-13, Hafizabad = γ-14, Mandi bahawdin = γ-15, Narowal = γ-16, Sialkot = γ-17, Lahore = γ-18, Kasur = γ-19, Nankanasahib = γ-20, Sheikhupura = γ-21, Multan = γ-22, Khanewal = γ-23, Lodhran = γ-24, Vehari = γ-25, Rawalpindi = γ-26, Attock = γ-27, Chakwal = γ-28, Jehlum = γ-29, Sahiwal = γ-30, Okara = γ-31, Pakptan = γ-32, Sargodha = γ-33, Bhakkar = γ-34, Khushab = γ-35, Mianwali = γ-36

Fertilizers= χ-1, Electric pumps = χ-2, Diesel pumps = χ-3, Tube wells = χ-4, Tractor = χ-5, Threshers χ-6, Harevester = χ-7, Mounter = χ-8, Total area sown = χ-9, Rice = χ-11, Cotton = χ-12, Sugarcane = χ-13, Wheat = χ-14.