Spatio-Temporal Evolution and Driving Factors of “Non-Grain Production” in Hubei Province Based on a Non-Grain Index

: Non-grain production (NGP) on cultivated land has become a common phenomenon due to the prosperity of the rural economy and the optimisation of the agricultural structure. However, the excessive use of cultivating land for NGP has threatened food production and the sustainable use of cultivated land. To halt this trend and to ensure food security, the authors of this paper applied a novel non-grain index to measure NGP, which could reﬂect multiple NGP activities; designated Hubei Province as its object of research; and revealed NGP’s spatio-temporal patterns of the past 30 years. We then assessed the characteristics of NGP based on spatial autocorrelation analysis, the Theil index, and geographically weighted regression. The results showed that the value of the non-grain index grew from 0.497 to 1.113 as NGP increased signiﬁcantly in Hubei Province. The number of high-NGP counties increased, spatial agglomeration became obvious, and the eastern and western sides of Hubei Province witnessed an observable growth in NGP. As a result, the NGP in the eastern and western regions overtook production in the central region. Despite a series of historical subsidy policies and agricultural modernisation initiatives that promoted the planting of grain crops, the policy of “grain on valuable cultivated land” could be better implemented. We conclude by making some suggestions for reducing NGP and protecting cultivated land.


Introduction
The COVID-19 pandemic and global climate change have exacerbated the pressure on global food security [1]. The United States, India, and Brazil are among the top four grainproducing countries with the highest number of COVID-19 cases and deaths, affecting grain planting, and countries have restricted food exports to cope with the shortage of food during the epidemic [2]. Global warming has also led to a lack of crops in many areas. From 1983 to 2019, 75% of wheat-growing regions, 82% of corn-growing regions, and 62% of rice-growing regions around the world had reduced yields due to drought [3]. Additionally, cultivated land is the principal carrier for grain production [4], but urbanisation and industrialisation have transformed cultivated land [5]. For example, a lot of cultivated land has been diverted to construction areas [6]. Some studies have suggested that grain crop production can only fulfil the demand for less than one-third of the population [7]. The report "2022 World Grain Safety and Nutrition Status" estimates that 828 million people worldwide are hungry [8,9].
As a country with scarce land per capita, food safety and cultivated land protection have always been the focus of the Chinese government's attention. The Chinese government has established a relatively comprehensive policy system for cultivated land protection in recent decades, which has helped reduce large-scale farmland loss due to urbanisation [10,11]. In addition to the dominant loss of cultivated land, the recessive loss of grain production capacity has become increasingly difficult to ignore [12]. Due to low grain crop planting income and the optimisation of agricultural structures, part of the cultivated land used to plant grain crops and cash crops has been turned into plant land-use model is the process of comparing grain-crop income and non-grain income. Moreover, cultivated land is widely distributed, with large differences in socioeconomic and resource endowment; some studies suggest that these factors have had different effects on the NGP in each region [13,21], leading to the spatial heterogeneity of NGP. Related agricultural policy has had the most rapid impact on NGP [30,31]. Factors such as grain subsidies and the optimisation of agricultural structures have driven the change in NGP in the short term. However, the research on factors influencing NGP is mainly based on the measurement of the grain ratio on a large scale, whereas the factors that affect Production NGP are often ignored. Despite a large number of studies, some limitations of the previous research are worth noting.
Owing to differences in the type and degree of the NGP in each region [27,28,36], it is necessary to reveal the internal mechanisms of NGP development and the spatio-temporal evolution of NGP. To measure NGP from multiple activities and provide policy suggestions for controlling NGP according to local conditions, we constructed a non-grain index based on the cultivated land areas versus the sown-grain areas, which is suitable for larger studies. In this study, we chose Hubei Province as our research object and revealed the regulation of the spatio-temporal dynamics of NGP from 1991 to 2019, and explored its tipping point. We applied the Theil index, autocorrelation analysis, and the geographically weighted regression model to explore the potential mechanisms of the spatio-temporal dynamics of NGP. Finally, we put forward suggestions for reducing non-grain activities and protecting the sustainability of cultivated land.

Study Area
Hubei Province is located in the centre of China and has various landforms: the highest altitude in western Shennongjia is 3105 m; the central Jianghan Plain is 35 m above sea level with flat terrain and fertile soil. Hubei Province is located in a subtropical zone, with an average temperature of 15 • C-17 • C and an average annual precipitation of 800-1600 mm (http://www.hubei.gov.cn). It receives plenty of light and precipitation. The water network is vertical and horizontal, with abundant water resources and excellent agricultural production conditions. Hubei Province's cultivated land area is 4.77 × 10 6 hm 2 , accounting for 3.73% of the total area of the country. It is an important producing and exporting grain region [37], and it plays a pivotal role in China's grain production. However, in recent years, NGP has become a salient issue in terms of food security, and there are different types and degrees of NGP in each area. For example, the type of NGP of Yichun County in 2020 was mainly planted non-grain crops. Among them, the planting area of medicinal material exceeded 2 × 10 4 hm 2 , accounting for 29.8% of the cultivated land area of Yichun County. Abandoned land is the main reason for the NGP of Xiangrong Village in Honghu City, with a wasteland area of 46.7 hm 2 , accounting for more than 40% of the cultivated land in the village. In Jianli, Qianjiang, and Honghu, the conversion from double-to single-season rice cropping systems (referred to as the rice "double-to-single" phenomenon) increased the risk of NGP. The rise in NGP in Hubei Province has seriously impacted the production of grain. According to the statistics in the China Statistical Yearbook and the Hubei Rural Statistical Yearbook, from 1991 to 2019, Hubei Province's grain production decreased from 5.16% to 4.10%. The annual output of food in Hubei Province decreased by 1.873 million tons, and the status of major grain production areas has continued to decline.

Data Collection and Processing
The data gathered in this study include socioeconomic statistics, land-use facts, and other online data. The data on the cultivated land area, sown-grain area, sown area of cash crops, irrigated area of cultivated land, population degree, agricultural employment, total power of agricultural machinery, and GDP are obtained from the China Statistical Yearbook (1992-2020) and the China Land and Hubei Rural Statistical Yearbook (1992-2020), which are available on the China National Knowledge Infrastructure Platform (https://www.cnki.net/). Administrative boundary data were sourced from the 3rd National Land Resource Survey on Hubei Province. Hubei Province's slope data and road data are derived from the Resource and Environment Science and Data Center of the Chinese Academy of Sciences (https://www.resdc.cn/), and the average annual precipitation data are derived from the monitoring data of the Hubei Meteorological Service (http://hb.cma.gov.cn/). In terms of data processing, to avoid outliers in the NGP measurements, we consolidated small municipal districts, such as merging Wuhan's seven principal urban areas. At the same time, in order to facilitate the compilation of statistics, we merged counties with changing administrative divisions. For example, the Duodao District of Jingmen City was divided into Dongbao District in 2001, and its statistics before 2001 are missing, so we merged them into Jingmen municipal district in this study. As a result, we obtained 83 subjects in 103 county-level divisions of Hubei Province, which we divided into four zones built on geographic location ( Figure 1). crops, irrigated area of cultivated land, population degree, agricultural employment, tot power of agricultural machinery, and GDP are obtained from the China Statistical Yearbo (1992-2020) and the China Land and Hubei Rural Statistical Yearbook (1992-2020), which a available on the China National Knowledge Infrastructure Platfor (https://www.cnki.net/). Administrative boundary data were sourced from the 3rd N tional Land Resource Survey on Hubei Province. Hubei Province's slope data and road da are derived from the Resource and Environment Science and Data Center of the Chine Academy of Sciences (https://www.resdc.cn/), and the average annual precipitation da are derived from the monitoring data of the Hubei Meteorological Servi (http://hb.cma.gov.cn/). In terms of data processing, to avoid outliers in the NGP mea urements, we consolidated small municipal districts, such as merging Wuhan's seve principal urban areas. At the same time, in order to facilitate the compilation of statistic we merged counties with changing administrative divisions. For example, the Duoda District of Jingmen City was divided into Dongbao District in 2001, and its statistics befo 2001 are missing, so we merged them into Jingmen municipal district in this study. As result, we obtained 83 subjects in 103 county-level divisions of Hubei Province, which w divided into four zones built on geographic location ( Figure 1).

Framework
In order to reveal the characteristics of the spatio-temporal evolution of NGP and explore its mechanisms and influencing factors, we constructed a framework as follow ( Figure 2). Step 1: Build datasets of NGP and construct our non-grain index based on th area of cultivated land and sown-grain area.
Step 2: Analyse the spatio-temporal evolutio characteristics of NGP with spatial autocorrelation and the Theil index.
Step 3: Identi the main influencing factors with an ordinary least squares (OLS) model; moreover, an lyse the spatial heterogeneity of factors using a geographically weighted regressio (GWR) model. Step 4: Put forward policy suggestions to control NGP.

Framework
In order to reveal the characteristics of the spatio-temporal evolution of NGP and to explore its mechanisms and influencing factors, we constructed a framework as follows ( Figure 2). Step 1: Build datasets of NGP and construct our non-grain index based on the area of cultivated land and sown-grain area.
Step 2: Analyse the spatio-temporal evolution characteristics of NGP with spatial autocorrelation and the Theil index.
Step 3: Identify the main influencing factors with an ordinary least squares (OLS) model; moreover, analyse the spatial heterogeneity of factors using a geographically weighted regression (GWR) model.
Step 4: Put forward policy suggestions to control NGP.

Non-Grain Index
NGP refers to changes in the use of cultivated land. The ratio of NGP is the proportion of cultivated land used for NGP to the total cultivated land of a region [12], and the

Non-Grain Index
NGP refers to changes in the use of cultivated land. The ratio of NGP is the proportion of cultivated land used for NGP to the total cultivated land of a region [12], and the formula is: where k is the ratio of NGP; S f is the area of cultivated land; S n is the area of cultivated land for NGP; and S G is the area of cultivated land for grain. Planting NGP (planting cash crops) and Production NGP (e.g., breeding, forestry, fruit farming, land abandonment, and leisure and tourism) results in an increase in the ratio of NGP. However, it is very difficult to obtain statistics for the "area of cultivated land for grain"; some studies use grain ratio (Formulas (2) and (3)) to measure Planting NGP (see Table 1), and some studies replace it with "sown-grain area/multiple-crop index" (Formulas (4) and (5)). Formally, the latter seems to be a good way to measure NGP, but it is not. The multiple-crop index is calculated by the ratio of the sown-grain area to the crop-sown area, ignoring the Production NGP. However, some cultivated land is used for breeding, forestry, and fruit farming instead of planting crops [38][39][40]. Thus, the multiple-crop index cannot express the recurrence of crops. Further, we transform Formula (5) and find that it is exactly the same as the grain ratio (Formula (6)). Meanwhile, previous measurements can only express the ratio of grain crops to cash crops, which are Planting NGP.

Formula
Description Reference Number S g is the sown-grain area, and S c is the crop-sown area. Expresses the ratio of grain to crops. [41] (2) Expresses the ratio of cash crops to crops. [35,42,43] (3) ρ is the multiple-crop index, but it cannot measure the recurrence of crops. [38][39][40] (4) The measurement result is the same as that of formula 3. [12,21,44] For example, the actual level of NGP was unchanged when the cultivated land was converted from Planting NGP to Production NGP, which are both NGP activities. However, due to the decrease in crop-sowing area, the grain ratio falsely indicates that NGP is reduced. Similarly, the actual level of NGP rises when the cultivated land is converted from crops (both grain and cash crops) to Production NGP, but the trend of grain ratio is related to the ratio of the sowing area of two crops, with possible increases or decreases or no changes. That is, when cultivated land is used for Production NGP, such as breeding, forestry, fruit farming, land abandonment, and leisure and tourism, there is a discrepancy between the estimate of grain ratio and the actual situation. In order to cover all types of NGP, we applied a so-called non-grain index (Formula (7)). The recurrence of crops is similar due to the natural environment in Hubei Province [38]. As the "area of cultivated land for grain" increases or decreases, the "sown-grain area" changes accordingly. Therefore, we replace "area of cultivated land for grain" with "sown-grain area" to reveal the evolution of NGP. As the "sown-grain area" may be larger than the cultivated area due to the recurrence of crops, to ensure a positive value of K (≥0), we add the constant term "a". According to the repeated cropping of crops in Hubei Province, "a" is set to 1, a high value for K indicates a high level for NGP, and the formula is: From the perspective of the non-grain index, when cultivated land used to plant grain crops turns into NGP, whether it is Planting NGP or Production NGP, the sown-grain area is reduced, S f − S g increases, and the non-grain index also increases, which indicates the growth of NGP. On the contrary, when S f − S g decreases, the non-grain index also decreases, indicating that NGP is shrinking. At the same time, when the recurrence of crops decreases, such as the rice "double-to-single" phenomenon, the decrease in the sown-grain area leads to an increase in S f − S g , and the index indicates the growth of NGP. The current approaches to NGP are complex, including planting cash crops, poultry breeding, tree planting, and abandoning cultivated land due to lack of a labour force, for example. In addition, the rice "double-to-single" phenomenon also increases the risk of NGP, but previous research paid less attention to it. Therefore, compared with the previous formulas that can only measure NGP for planting, a novel non-grain index utilised to measure NGP, which could reflect multiple NGP activities, has great significance for revealing how NGP develops in Hubei Province.

Spatial Autocorrelation Analysis
Spatial autocorrelations include global autocorrelations and local autocorrelations, which test whether observations are correlated between units. This study used the global Moran's I index to reveal the overall distribution of NGP in Hubei Province. A Moran's I value represents a range between −1 and 1. An increase in spatial correlation is observed when the absolute index value approaches 1. Positive values indicate positive spatial correlation, whereas negative values indicate negative spatial correlation. If Moran's I is equivalent to 0, there is no spatial correlation, and only randomness exists [22]. The formula is: where n is the number of counties in Hubei Province; y i and y j are the non-grain index values of counties i and j, respectively; y is the mean of y i and y j ; and the spatial weight matrix w ij indicates the strength of the relationship between counties i and j.

Theil Index Model
The Theil index is derived from a generalised entropy measure [45]. The Theil index can illustrate the characteristics of additive decomposition. The overall difference relates to intra-and inter-regional differences. It is possible to measure how much both contribute to the overall difference [46,47]. The formula is: where n is the number of counties in Hubei Province, y p is the non-grain index value of county p, and u y is the non-grain index value of the province. The value for a Theil index is between 0 and ln n. If the non-grain index of each county is the same, the value of the Theil index is 0, illustrating that there is no difference between counties. If a county has a small area of cultivated land and a large sown-grain area, the value of the Theil index is ln n. The relationship between the intra-regional differences (T wr ) and inter-regional differences (T br ) is: The formula for inter-regional differences (T br ) is: where n is the number of city-level units in Hubei Province; p i is the number of county-level units of one city; p is the number of counties in Hubei Province; y i is the non-grain index value of city i; and µ is the non-grain index value of the province.

GWR Mode
The GWR model represents a development of the OLS test. GWR can reflect changes in the relationship between independent variables and dependent variables with different geographical locations [48,49]. We examined the spatial heterogeneity of the influence of factors on NGP by embedding the non-grain index and factor data into regression models according to spatial location [17,50]. The formula for this calculation is: where y i is the non-grain index value of the county-level unit; (u i , v i ) is the geographic coordinate of county i; β 0 is the regression constant in county i; x im is the value of the interpretation variable x m in county i; the function β m (u, v) has the value β m (u i , v i ) in county i; m is the number of variables affecting the non-grain index; and ε i is the error term. Before GWR modelling, factor data were normalised to eliminate dimensional differences among independent variables. The data processing formula is: where X imax and X imin are the maximum and minimum values of interpretation variable x i , respectively. If the regression coefficient is positive, it means that the explanatory variable has a positive impact on the non-grain index, and if the regression coefficient is negative, it means that the explanatory variable has a negative impact on the non-grain index.

Indicators Selection
Numerous factors affect NGP. Via cultivated land operators, the use of cultivated land is affected by factors such as the production capacity of cultivated land and the efficiency of crop planting. This study refers to previous research results (see Table 2), and, taking into account the availability and classicality of factor data, we identified the principal selection factors of NGP. These factors can better explain the impetus driving NGP from different viewpoints, including the natural environment, socioeconomic development, and agricultural modernisation perspectives [51]. Resource endowment is an important factor that affects the use of cultivated land. Rice needs to grow in an environment with sufficient water and a suitable temperature [58]. Crops are difficult to grow on sloping cropland because of the lack of water and soil. Dry land is not conducive to the growth of crops, resulting in a decline in production [59]. This study selected the average annual rainfall, the slope of cultivated land, and the proportion of dry land as the natural condition indicators.
Socioeconomic development pertains to the economic benefits of agriculture [1], and the relatively low monetary benefit of grain growing is an important factor that leads to NGP. Normally, there is a huge gap between grain income and average income in economically developed areas. Farmers prefer to grow commercial crops with higher investment income. Complete transportation facilities can reduce the costs of transportation of agricultural products [60], and the area of cultivated land per capita also affects a farmer's income, which influences the land-use decision of cultivated land operators [61]. This study selected per capita cultivated land area, per capita GDP, and road density as the socioeconomic indicators.
Agricultural modernisation is an important consideration affecting the efficiency of agricultural production, and mechanisation of agriculture has led to a shift in production methods, which can effectively replace family labour [42]. Irrigation technology provides crops with water and maintains the production of crops. Agricultural modernisation can increase a farmer's income and enlarge their asset investment, thus changing incomes and affecting the costs and choices of land operators. This study selected machinery power and the proportion of irrigated cultivated land as indicators to measure the level of agricultural modernisation.

Temporal Dynamics Changes in NGP
The value of the non-grain index grew from 0.497 to 1.113 from 1991 to 2019, while NGP showed noticeable growth in Hubei Province ( average yearly increase of 0.118, while NGP showed a trend of first rising, then falling and rising again in Hubei Province. Further observing the changes in NGP from the partitions of Hubei Province, from 1991 to 2005-except for northern Hubei Province-the nongrain index continued to rise, and the NGP remained the highest in southern Hubei Province during this period. From 2005 to 2019, the non-grain index of the four regions decreased first and then increased; NGP decreased the most in the southern and northern regions, whereas it increased the most in the eastern and western regions. The value of the non-grain index in the eastern and western regions increased by 0.716 and 0.985, respectively. Whereas the NGP in the northern and southern regions mainly increased at low rates, the value of the non-grain index of the two regions increased by 0.457 and 0.475, respectively. To reveal the origin of the differences in the NGP of Hubei Province, the non-grain index was decomposed by using the Theil index to describe the differences within (intragroup) and between (intergroup) prefecture-level cities. The results are provided in Figure 4. It is evident that NGP differs across Hubei Province. The Theil index decreased from 0.123 to 0.077, and the total difference in the non-grain index among counties decreased gradually. From the perspective of the composition of the Theil index, the NGP's aggregate function is mainly due to the differences between prefecture-level cities, which have average, maximum, and minimum contribution rates of 73.4%, 88.7%, and 64.9%, respectively. The differences between prefecture-level cities are the main source of the total difference observed. In 2001, the NGP difference between Wuhan and Huangshi was the largest, with the difference between the value of their non-grain indexes being 0.676. The difference within prefecture-level cities for the levels of NGP first increased and then decreased before finally stabilising at a low level. After 2001, it stabilised at about 0.22, while its contribution to the total difference fluctuated between 20.05% and 40.9% due to changes in the differences between prefecture-level cities. This indicated that the variation in the total difference of NGP was mainly due to the variation in the difference between prefecture-level cities. The total difference reached its minimum in 2019, indicating that the NGP of each county was close. To reveal the origin of the differences in the NGP of Hubei Province, the non-grain index was decomposed by using the Theil index to describe the differences within (intragroup) and between (intergroup) prefecture-level cities. The results are provided in Figure 4. It is evident that NGP differs across Hubei Province. The Theil index decreased from 0.123 to 0.077, and the total difference in the non-grain index among counties decreased gradually. From the perspective of the composition of the Theil index, the NGP's aggregate function is mainly due to the differences between prefecture-level cities, which have average, maximum, and minimum contribution rates of 73.4%, 88.7%, and 64.9%, respectively. The differences between prefecture-level cities are the main source of the total difference observed. In 2001, the NGP difference between Wuhan and Huangshi was the largest, with the difference between the value of their non-grain indexes being 0.676. The difference within prefecture-level cities for the levels of NGP first increased and then decreased before finally stabilising at a low level. After 2001, it stabilised at about 0.22, while its contribution to the total difference fluctuated between 20.05% and 40.9% due to changes in the differences between prefecture-level cities. This indicated that the variation in the total difference of NGP was mainly due to the variation in the difference between prefecture-level cities. The total difference reached its minimum in 2019, indicating that the NGP of each county was close.

Spatial Dynamics Changes in NGP
From the perspective of spatial patterns, the number of high-NGP counties (nongrain index > 1) in Hubei Province increased from 4 to 62, demonstrating the expansion technique of infilling (see Figure 5). The high-NGP area moved from the middle to the south and then to the east, finally forming a spatial pattern of concentration in the east and a multi-core distribution pattern in the west. In 1991, there were four counties that could be classified as high-NGP counties, which were relatively concentrated in the inner part of Hubei Province:

Spatial Dynamics Changes in NGP
From the perspective of spatial patterns, the number of high-NGP counties (non-grain index > 1) in Hubei Province increased from 4 to 62, demonstrating the expansion technique of infilling (see Figure 5). The high-NGP area moved from the middle to the south and then to the east, finally forming a spatial pattern of concentration in the east and a multi-core distribution pattern in the west. In 1991, there were four counties that could be classified as high-NGP counties, which were relatively concentrated in the inner part of Hubei Province: From the perspective of spatial patterns, the number of high-NGP counties (nongrain index > 1) in Hubei Province increased from 4 to 62, demonstrating the expansion technique of infilling (see Figure 5). The high-NGP area moved from the middle to the south and then to the east, finally forming a spatial pattern of concentration in the east and a multi-core distribution pattern in the west. In 1991, there were four counties that could be classified as high-NGP counties, which were relatively concentrated in the inner part of Hubei Province: Huangshi Municipal District (1.618), Yichang Municipal District (1.077), Wuhan Municipal District (1.072), and Xiangyang City District (1.039). In 2000, the NGP in central Hubei Province increased significantly; 14 counties could be classified as high-NGP counties, which were relatively concentrated in the flat areas in central Hubei Province: Xiangyang City District, Wuhan City District, Huangshi City, Yichang City, Hanchuan City, Qianjiang City, Laohekou City, Xiantao City, Yicheng, Jiayu County, Dangyang City, Zhijiang City, Shayang County, and Honghu County. In 2010, 32 counties could be considered high-NGP counties, which were relatively concentrated on the Jianghan Plain in southern Hubei Province. The low-rate expansion of NGP units was still observed in the northeast at this time, while NGP was high in the south and low in the north, including areas in Huangpi District, Fang County, Changyang County, Xuan'en County, and other counties. In addition, four counties could be considered to be low-NGP counties: Zaoyang City (0.754), Xiangzhou District (0.745), Honghu City (0.744), and Jiangling County (0.701). In the past 30 years, the eastern and western parts of Hubei Province witnessed an obvious growth of NGP; as a result, the NGP in the eastern and western regions overtook the NGP of the central region.  As shown in Figure 6, from 1991 to 2019, NGP showed a significant spatial correlation, and the Moran's I value was greater than 0 in all periods. In 2010, the Moran's I value was 0.461, which indicated that the spatial agglomeration effect of NGP was the most pronounced. As shown in Figure 6, from 1991 to 2019, NGP showed a significant spatial correlation, and the Moran's I value was greater than 0 in all periods. In 2010, the Moran's I value was 0.461, which indicated that the spatial agglomeration effect of NGP was the most pronounced.   Generally, the spatial agglomeration types of the NGP in Hubei Province are mainly high-high clustering (H-H) and low-low clustering (L-L), indicating that the NGP had a tough radiation effect (Figure 7). The distribution of H-H changed from multi-core to single-core and moved eastward during the past 30 years, moving south and then east. In 1991, H-H was assigned to central Hubei Province, including Shayang County, Honghu City, Jianli City, and six other counties. In 2000, H-H included 19 counties, which were still mainly located in central Hubei Province, indicating an expansion trend. In 2000, H-H was mostly concentrated on the Jianghan Plain in southern Hubei Province, in 10 counties. In 2019, H-H moved to the eastern part of Hubei Province, and the clustering was relatively concentrated in the surrounding areas of Caidian District, Xinzhou District, and Ezhou City. The Jianghan Plain gradually changed from H-H to L-L clustering, with the radiation effect in traditional agricultural areas shifting from NGP to growing grain.

Results of the GWR
In order to test the multi-collinearity of various factors, a large number of factors were selected as preliminaries. After calculating the variance inflation factors (VIFs) for the various driving factors using OLS, the significant collinearity factors with a VIF > 7.5 were eliminated, as were the invalid and dominant factors. The VIF values of the various driving factors were all below 5.0, suggesting that there were no redundant variables in the model and no multi-collinearity among certain factors. Various factors were at a high significance level, and the values of X2, X5, and X8 passed the significance test at the 0.05

Results of the GWR
In order to test the multi-collinearity of various factors, a large number of factors were selected as preliminaries. After calculating the variance inflation factors (VIFs) for the various driving factors using OLS, the significant collinearity factors with a VIF > 7.5 were eliminated, as were the invalid and dominant factors. The VIF values of the various driving factors were all below 5.0, suggesting that there were no redundant variables in the model and no multi-collinearity among certain factors. Various factors were at a high significance level, and the values of X2, X5, and X8 passed the significance test at the 0.05 level, while the value of X7 passed the significance test at a level of 0.01 (see Table 3). Therefore, we added the above four indicators into our GWR model to analyse their spatial effect on the influencing factors. According to Table 4, the R 2 and adjusted R 2 values of the GWR model are higher than the OLS model's finding, and the Akaike information criterion (AIC) value of the GWR model is much smaller than that of the OLS model. These results indicate that the GWR model has a better goodness of fit, as it produced a better spatial coefficient regression than the OLS model did. The factors with a significant correlation in the regression of 83 counties and the related effects on the spatial heterogeneity of NGP were analysed and then visualised with the natural breaking point method (see Figure 7). As shown in Figure 8a, the influence of the natural environment on NGP was spatially heterogeneous. The regression coefficient of the factor of NGP increased gradually from east to west. The slope of cultivated land exhibited a positive correlation with NGP in 37 counties; it showed a negative correlation in 46 counties. The former is mostly located in the counties of western Hubei, where the altitude is high, the temperature difference between day and night is considerable, and the increase in cultivated land slope improved the NGP. In Zhushan County and Gucheng County, with the highest regression coefficient, planting cash crops is the main type of NGP, and walnut planting is the native characteristic agriculture. The latter was mainly located in the flat eastern parts of Hubei, which have higher slopes to grow more grain. NGP increased in the flat cultivated land areas, which was particularly significant in Ezhou City, where the annual output of vegetables exceeded 1 million tons.

Spatio-Temporal Characteristics of NGP and Effects
In the past 30 years, the value of the non-grain index has shown an upward trend in Hubei Province, indicating that the NGP has been continuously increasing in Hubei Province, although it is a major grain-producing area. However, some studies suggested that the NGP in major grain-producing areas declined in recent years [12]. The main reason is that the grain ratio formula used in other studies can only indicate that the proportion of grain sown in crops in Hubei Province has increased [35,[41][42][43], neglecting the increase in NGP caused by breeding on cultivated land. NGP growth accelerated after 2016, and high-NGP counties spread across the province. Other studies reached the same conclusion [23]. This was mainly affected by the adjustment policy of the agricultural structure in 2015, where a transformation of agricultural structures and interventions in industrial sectors resulted in an increase in NGP [42]. From the perspective of the decline in the national proportion of grain, NGP has become a salient issue in Hubei Province.
To be specific, the NGP of the central region was less than that of the eastern and As shown in Figure 8b, the influence of GDP on NGP gradually increased from east to west in Hubei Province. In most areas of central and eastern Hubei Province, the correlation coefficient between the economic development level and NGP was around 0, suggesting that economic progress had no obvious impact on NGP. The regions with a significant positive correlation between GDP and NGP were concentrated in northwest Hubei Province, where agricultural income accounts for a large proportion of GDP. The rural revitalisation strategy has promoted local development as well as NGP.
As shown in Figure 8c,d, agricultural modernisation exhibited a negative correlation with NGP. The impact of two indicators on NGP, including the average amount of agricultural machinery and the proportion of irrigated land, increased from west to east in Hubei Province, while agricultural modernisation has gradually promoted the cultivation of grain. This demonstrated that the operating entity was more inclined to plant grain crops on cultivated land with a high level of irrigation and mechanisation. The eastern part of Hubei Province has a higher level of agricultural modernisation. In addition, Yunmeng County, Dongxihu District, Qianjiang City, and Jianli City are agricultural modernisation demonstration zones where the increase in NGP was effectively controlled.

Spatio-Temporal Characteristics of NGP and Effects
In the past 30 years, the value of the non-grain index has shown an upward trend in Hubei Province, indicating that the NGP has been continuously increasing in Hubei Province, although it is a major grain-producing area. However, some studies suggested that the NGP in major grain-producing areas declined in recent years [12]. The main reason is that the grain ratio formula used in other studies can only indicate that the proportion of grain sown in crops in Hubei Province has increased [35,[41][42][43], neglecting the increase in NGP caused by breeding on cultivated land. NGP growth accelerated after 2016, and high-NGP counties spread across the province. Other studies reached the same conclusion [23]. This was mainly affected by the adjustment policy of the agricultural structure in 2015, where a transformation of agricultural structures and interventions in industrial sectors resulted in an increase in NGP [42]. From the perspective of the decline in the national proportion of grain, NGP has become a salient issue in Hubei Province.
To be specific, the NGP of the central region was less than that of the eastern and western parts of Hubei Province in 1991. The growth of the non-grain index in the eastern and western regions exceeded that in the other regions in the past 30 years. As a result, the NGP of the eastern and western regions eventually surpassed that of the focal region. The observed low-rate expansion of the NGP in northern and southern Hubei Province was primarily due to the flat terrain, fertile soil, and abundant water resources in this region, which are conducive to crop growth. Including the Jianghan Plain, the central and southern parts of Hubei Province are important agricultural production areas. The Chinese government's policy of "grain on valuable cultivated land" promoted farmers' enthusiasm to plant grain [62]. The differences in resource endowment are the main reason for the rise of NGP in western Hubei Province [63], where the climate is dry, and the high altitude is not suitable for growing grain. Counties in Enshi Prefecture grow organic vegetables and Chinese herbs due to the advantageous towering mountains, which have deep valleys and large temperature differences between day and night. Related analyses showed that the NGP of Hubei Province presents obvious spatial agglomeration results. The H-H clustering of NGP occurred in Wuhan, Xiaogan, and Ezhou City in the east of Hubei Province, which might result from the development of the Wuhan Metropolitan Area [64], promoting the surrounding population to gather in that central urban area. A loss of labour may lead to NGP or even the abandonment of cultivated land. At the same time, the Wuhan Metropolitan Area has attracted a large number of enterprises that are involved in cultivated land management. The profit-seeking behaviour of enterprises has led to the rise of NGP in eastern Hubei Province, and, along with promoting rural revitalisation, a large amount of cultivated land formerly used to plant grain was converted to cash crops, breeding, forestry, and fruit farming [21]. The Theil index analysis showed that the NGP differences between counties decreased, and an intergroup difference was the main component of the total difference, suggesting that the geographical differences of natural, economic, and social factors are important reasons for the spatial differences in NGP. The increase in the NGP level in Hubei Province has a direct impact on grain production. From 1991 to 2019, the proportion of grain production in China decreased by 1.06% from 5.16% to 4.1% (see Figure 9), and the growth rate of grain production in Hubei has continued to decline. Since 2012, Hubei's per capita share of grain has reached the tipping point of food security (per capita share of grain reached 400 kg [65]), but it has recently declined by more than 33 kg from the peak in 2015. Excessive NGP may pass the tipping point for food security, but the threat is obscured by the stabilisation of total food production.
tion. From 1991 to 2019, the proportion of grain production in China decreased by 1.06% from 5.16% to 4.1% (see Figure 9), and the growth rate of grain production in Hubei has continued to decline. Since 2012, Hubei's per capita share of grain has reached the tipping point of food security (per capita share of grain reached 400 kg [65]), but it has recently declined by more than 33 kg from the peak in 2015. Excessive NGP may pass the tipping point for food security, but the threat is obscured by the stabilisation of total food production.

Mechanism of Factors Affecting NGP
According to previous research, in addition to the relatively low income from grain planting, the low economic efficiency of grain planting [18], the intervention of enterprise and commercial capital [66], the resource endowment [67][68][69] of cultivated land, and related policies all led to NGP in China [12]. Using a GWR model, this study explored the effects of social, economic, and natural factors on NGP, and it showed the effects of various factors on NGP in different counties.
Overall, our research showed that the influence of resource endowment on NGP correlates with obvious spatial differences [70], similar to a prior conclusion [42]. The positive correlation between cultivated land slope and NGP is mainly located in western Hubei [71], which is mountainous and located at a high altitude. Baokang is a famous

Mechanism of Factors Affecting NGP
According to previous research, in addition to the relatively low income from grain planting, the low economic efficiency of grain planting [18], the intervention of enterprise and commercial capital [66], the resource endowment [67][68][69] of cultivated land, and related policies all led to NGP in China [12]. Using a GWR model, this study explored the effects of social, economic, and natural factors on NGP, and it showed the effects of various factors on NGP in different counties.
Overall, our research showed that the influence of resource endowment on NGP correlates with obvious spatial differences [70], similar to a prior conclusion [42]. The positive correlation between cultivated land slope and NGP is mainly located in western Hubei [71], which is mountainous and located at a high altitude. Baokang is a famous mountainous county with more than 3100 hilltops, so it is more suitable for forestry and fruit farming since the singular topography promotes the development of NGP. The central and eastern parts of Hubei Province are abundant in rainfall and heat, which help the growth of crops [72]. Jianghan Plain, which has a large number of agricultural exports, has fertile soil and an extensive water network. However, depending on the results of our GWR model, cultivated land slope was negatively correlated with NGP. This indicates that cultivated land operators tend to choose NGP with superior returns on cultivated land that has a gentle slope and good natural quality. The policy of "grain on valuable cultivated land" has controlled the increase in NGP on the Jianghan Plain [62], but this policy still has room to grow in terms of its further implementation in the southern and eastern parts of Hubei Province.
The regions with a significant positive correlation between GDP and NGP are concentrated in the northwest of Hubei Province, indicating that an increase in GDP promotes NGP. The GDP of Zhuxi County in 2019 was RMB 8.45 billion, and the income of the primary industry accounted for more than 27% of the total. Sending industrial and commercial capital to the countryside promoted the development of local economies, while profit-seeking enterprises converted to cash crops [73,74]. In other areas of Hubei Province, the correlation coefficient between economic progress and NGP was not significant. There are two different groups of studies on this: some suggest that economically developed regions have higher prices of seeds, pesticides, fertilisers, and higher rents for cultivated land [12]. However, grain purchase prices remain relatively stable, which reduces farmers' enthusiasm to grow grain. At the same time, the secondary and tertiary industries in this region developed rapidly, so agricultural production is no longer part of a feeding option for cultivated land operators [41]. Other studies found that economic development had a completely opposite effect on food cultivation in Zhejiang and Jiangsu Provinces [75]. The rise in labour costs instead promoted the cultivation of rice. Combined with our research results, we suggest that the influence of economic progress on NGP is complicated and does not always produce one result [42].
There was a negative correlation between agricultural modernisation and NGP in Hubei Province. The influence of the proportion of agricultural machinery and irrigated cultivated land in NGP increased gradually from west to east in Hubei Province, showing that agricultural modernisation had an obvious control effect on NGP. The reason is that agricultural modernisation has increased the efficiency of crop production and encouraged farmers to grow grain. Since 2003, local governments, including Hanchuan City, have provided subsidies to farmers to purchase farm machinery, thus reducing their financial burden [62]. Agricultural modernisation has had the most obvious effect on NGP in the central and eastern regions, where the degree of agricultural modernisation is high. Yunmeng County, Dongxihu District, Qianjiang City, and Jianli City were identified in this region as agrarian modernisation demonstration zones. This shows that the higher the level of agricultural modernisation is, the more obvious the control effect on NGP is. However, some studies came to a different conclusion: agricultural modernisation significantly improves NGP [42,43,63]. Utilising agricultural machinery has increased agricultural costs, making farmers more inclined to plant cash crops with higher returns [13]. This result is different because previous studies used the grain ratio to measure NGP and ignored the variation of Production NGP [42,43]. Improvements in crop efficiency have resulted in the conversion of land from ponds and wastelands to crops. The grain ratio declines indicate that the proportion of newly sown food crops has decreased. Nevertheless, from the result of the non-grain index, the area of cultivated land for grain is increasing, and agricultural modernisation promotes the planting of grain crops.
Moreover, multiple cropping of grain crops will be captured by the non-grain index; we conducted further data validation to identify this effect. The sown area of rice tends to be stable before 2016 (see Figure 10a), and the average area of "double-to-single" was 6.9 × 10 3 hm 2 (see Figure 10b), accounting for 2.14% of the average grain sowing area change, which had no obvious influence on the index. In 2016, "double-to-single" covered an area of 1.98 × 10 5 hm 2 , accounting for 70.7% of the average grain sowing area change. As a result, the measurement of the non-grain index is higher than the actual NGP level in 2016. As a result, the measurement of the non-grain index is higher than the actual NGP level in 2016. From the actual situation, the "double-to-single" phenomenon in Hubei Province is caused by the conversion of double-season rice to rice-crayfish farming and is most obvious in Jianli, Qianjiang, and Honghu. Because of the lack of relevant regulations, some farmers are digging pits to expand the breeding space on cultivated land. The destruction of the arable layer leads to a decline in rice yield and affects the sustainable utilisation of cultivated land [76].

Conclusions
In the past 30 years, grain production in Hubei Province has declined, and NGP has threatened food security. This study introduces a non-grain index based on the cultivated land area and sown-grain area with Hubei Province as the research object and reveals the spatio-temporal patterns of the past 30 years. We then assessed the characteristics of NGP based on the Theil index and a spatial autocorrelation analysis and explored the effects of different factors on NGP through the application of a GWR model.
Conclusions: (1) For nearly 30 years, while the value of the non-grain index grew from 0.497 to 1.113, NGP has not reached the tipping point for food security in Hubei Province, but its threat is growing. (2) High-NGP counties are increasing in number, and their spatial agglomeration is obvious. The east and west of Hubei Province witnessed an observable growth of NGP. As a result, the NGP in the eastern and western regions overtook the central region's NGP. (3) The factors with the greatest impact on the spatial differentiation of NGP were the slope of cultivated land, per capita GDP, average agricultural machinery used, and proportion of irrigated land. Generally, although a series of subsidy policies and agricultural modernisation promoted the planting of grain crops, NGP is still increasing in Hubei Province. The policy of "grain on valuable cultivated land" still has From the actual situation, the "double-to-single" phenomenon in Hubei Province is caused by the conversion of double-season rice to rice-crayfish farming and is most obvious in Jianli, Qianjiang, and Honghu. Because of the lack of relevant regulations, some farmers are digging pits to expand the breeding space on cultivated land. The destruction of the arable layer leads to a decline in rice yield and affects the sustainable utilisation of cultivated land [76].

Conclusions
In the past 30 years, grain production in Hubei Province has declined, and NGP has threatened food security. This study introduces a non-grain index based on the cultivated land area and sown-grain area with Hubei Province as the research object and reveals the spatio-temporal patterns of the past 30 years. We then assessed the characteristics of NGP based on the Theil index and a spatial autocorrelation analysis and explored the effects of different factors on NGP through the application of a GWR model.
Conclusions: (1) For nearly 30 years, while the value of the non-grain index grew from 0.497 to 1.113, NGP has not reached the tipping point for food security in Hubei Province, but its threat is growing. (2) High-NGP counties are increasing in number, and their spatial agglomeration is obvious. The east and west of Hubei Province witnessed an observable growth of NGP. As a result, the NGP in the eastern and western regions overtook the central region's NGP. (3) The factors with the greatest impact on the spatial differentiation of NGP were the slope of cultivated land, per capita GDP, average agricultural machinery used, and proportion of irrigated land. Generally, although a series of subsidy policies and agricultural modernisation promoted the planting of grain crops, NGP is still increasing in Hubei Province. The policy of "grain on valuable cultivated land" still has room for implementation, while Hubei Province's status as a major grain-producing region continues to decline. Given the prevalence of NGP activities in China, therefore, the following suggestions are put forward for reducing NGP and achieving the sustainable use of cultivated land.
First, the government should give more subsidies to grain growers and guide farmers and enterprises to plant grain crops. The important factor affecting the decision of a cultivated land operator is economic benefit. Increasing grain-planting subsidies and the purchase price of grain crops could enhance farmers' enthusiasm for growing grain. Moreover, guiding enterprises to plant grain crops by limiting the use of cultivated land and specifying penalties for violations could also be beneficial.
Second, the results of this study suggest that agricultural modernisation can guide breeding and other Production NGP to be converted to grain crops. On the one hand, this can be accomplished by reducing the costs of agricultural modernisation and promoting irrigation facility subsidies and agricultural machinery purchase subsidies. On the other hand, fixing finely fragmented land facilitates the application of heavy agricultural machinery, thereby promoting the appropriate scale and modernisation of grain farming. This conclusion has important implications for Production NGP areas such as Xinjiang, Qinghai, and other places.
Finally, we suggest determining cultivated land use according to local natural environments, for example, by implementing clear regulations for rice-crayfish farming in JianLi, Qianjiang, and Honghu. This can be accomplished by ensuring that grain is grown on valuable cultivated land, which limits the use of cultivated land with high production efficiency according to law, and strictly controlling NGP to protect domestic food security. At the same time, the restrictions on the use of cultivated land with poor irrigation conditions and high altitude can be relaxed, such as in western Hubei Province. This would reduce inefficient labour productivity and increase capital investment.
The limitation of the non-grain index is that its measurements overestimate the NGP when "double-to-single" is evident. At the same time, the specific impact of rice-crayfish farming on food security remains unclear, and it lacks relevant regulations. Therefore, it is recommended to focus future research activities on countermeasures and threats to "double-to-single".