The empirical results included the following three parts: the analysis of the farmers’ preferences and WTA for participating in the fertilizer reduction program, the heterogeneity analysis of farmers’ preference, and the potential classification of farmers.
5.1. Farmers’ Preference for Chemical Fertilizer Reduction Programs
We used Stata 16.0 to conduct an empirical analysis of the farmers’ choice experiment data with the random parameter logit model. The coefficients of the annual cash subsidy in the fertilizer reduction scheme were set as random parameters with normal distribution, and the coefficients of other attribute variables were set as fixed parameters. Then, 1000 Halton sampling was used for estimation. The estimation results are shown in
Table 4, where Model (1) only includes attribute variables, and Model (2) includes attribute variables and household characteristic variables. The coefficients of the farmers’ characteristic variables in Model (2) are based on the “status quo” scheme. Since these two schemes are designed and randomly set, the coefficients and significance obtained by regression are basically the same. We reported a set of coefficients that shows the influence of the farmers’ characteristics on their participation in the fertilizer reduction scheme.
It can be seen from
Table 4 that the results of Model (1) and Model (2) show that the coefficients of the attribute variables of contract length, land area, and fertilizer use are significantly negative, indicating that the longer the program implementation period, the more the increase in the proportion of participating land area and the increase in the proportion of reducing fertilizer will reduce the willingness of farmers to participate. The coefficient of the attribute variable of the annual cash subsidy is significantly positive, indicating that providing cash compensation can improve farmers’ willingness to participate.
The influence of farmers’ characteristics on their participation in the fertilizer reduction program was represented by the results of Model (2). In terms of individual characteristics, male farmers were more willing to participate in the fertilizer reduction program. Age, education level, and whether they were village leaders did not affect the willingness to participate. In terms of household characteristics, the coefficient of the proportion of non-agricultural income of farmers was significantly negative, indicating that the higher the proportion of non-agricultural income of farmers, the lower the willingness of farmers to participate. In terms of farmers’ farmland management, the higher the amount of fertilizer per ha, the lower the willingness of farmers to participate in the fertilizer reduction program. Farmers may think that a higher amount of fertilizer can maintain a greater crop yield, making them less willing to participate in the fertilizer reduction program. The shorter the land distance from Poyang Lake, the higher the farmers’ attention to Poyang Lake, and the higher their willingness to participate. In terms of cognition, the higher the farmers’ cognition of the harm of chemical fertilizer was, the more willing they were to participate. However, farmers’ cognition of the importance of wetland did not significantly affect their willingness to participate in the fertilizer reduction program.
In order to ensure the accuracy of the benchmark regression results, we switched to the conditional logit model and the multiple logit model to test the robustness of the experimental data regarding the farmers’ selection. Models (3) and (4) in
Table 5 represent the conditional logit regression and multiple logit regression, respectively, and the regression results are shown below. It can be seen that the regression results of the conditional logit and multiple logit are not much different from the benchmark regression results, especially for the annual cash subsidy attribute variable that we focused on.
According to the parameter estimation results of the regressions and Formula (9), farmers’ willingness to accept the attributes of the fertilizer reduction program were calculated, as shown in
Table 6. Taking the results of Model (2) as an example, regarding the implementation years of the fertilizer reduction program, for each additional year, the farmers’ WTA was 63.75 CNY/ha/year. For every 10% increase in the proportion of farmers participating in the land area, the WTA was 73.875 CNY/ha/year. For every 10% reduction in fertilizer application, the WTA was 413.505 CNY/ha/year.
The results of this study were compared with those of the related literature, as shown in
Table 7. It can be seen that the WTA calculated in this paper is close to the calculated results of Li [
18] and Zhu [
22], which indicates that the results of this paper are reliable to some extent.
5.2. Individual Characteristics of Farmers and Preferences
We further analyzed the heterogeneity of farmers’ preferences to participate in fertilizer reduction programs with respect to individual characteristics. On the basis of Model (2), the interaction terms of the attribute variables and individual characteristic variables were added, in which the individual characteristics were household support burden, proportion of non-agricultural income, average amount of chemical fertilizer per ha, farmers’ cognition of chemical fertilizer harm, and farmers’ cognition of wetland importance. The regression results are shown in
Table 8 and
Table 9. In general, only a few of the coefficients of the interaction terms between the individual characteristic variables and the attribute variables were significant; that is, some of the characteristic variables of the farmers showed heterogeneity in the WTA values for the attributes of the program.
In Model (5), which included the interaction terms between the household support burden and attribute variables, only the interaction term between household support burden and land area was significant and negative at the level of 10%. Based on Model (3), the farmers’ WTA values for the attributes of the fertilizer reduction program were calculated. The calculations are as follows: First, the interaction coefficient between the land area and household support burden was multiplied by the sample mean of the household support burden, plus the coefficient of land area (if significant), and then divided by the coefficient of the annual cash subsidy, and finally, the absolute value of the result is the farmers’ WTA for the attributes of the fertilizer reduction program. The results showed that, for every 10% increase in the proportion of farmers’ participated land area, the WTA was 39 CNY/ha/year, while the WTA increased by 14.145 CNY/ha/year for every 10% increase in the support burden of the farmers.
In Model (6), which includes the interaction terms between the proportion of non-agricultural income and the attribute variables, only the interaction term between the proportion of non-agricultural income and the proportion of fertilizer reduction is significant and negative; that is, the higher the proportion of non-agricultural income, the higher the WTA value for fertilizer reduction. Even if farmers can obtain higher non-agricultural income, they still hope to maintain the current agricultural income and are not willing to reduce the amount of chemical fertilizer used. Based on Model (4), we calculated the farmers’ WTA for reducing the usage of chemical fertilizers. It can be seen that for every 10% reduction in fertilizer use, the WTA is 441.15 CNY/ha/year, and for every 10% increase in the ratio of non-agricultural income of the household, the WTA rises by 70.8 CNY/ha/year.
In Model (7), which includes the interaction terms of average fertilizer use per hectare and attribute variables, the coefficients of the three interaction terms were not significant. In other words, the amount of chemical fertilizer per ha did not affect the WTA of farmers in terms of contract length, land area, or fertilizer reduction.
In Model (8), which includes the interaction terms between farmers’ cognition of wetland importance and attribute variables, the regression coefficients of two interaction terms are significant. Among them, the coefficient of the interaction term between the farmers’ cognition of the importance of wetland and the contract length is significantly negative, which indicates that the higher the farmers’ cognition of wetland importance, the higher their WTA value for the implementation period of the program. That is, although the farmers recognize the importance of wetland protection, they still believe that the implementation period of the program will have a negative impact on their agricultural yield. Therefore, they hope to obtain a higher level of economic compensation. The coefficient of the interaction term between the farmers’ cognition of wetland importance and the proportion of participating land area was significantly positive, indicating that the higher the farmers’ cognition of wetland importance, the lower their WTA value for the proportion of participating land area; that is, they are more willing to invest a larger proportion of land in the fertilizer reduction program to strengthen wetland protection. In addition, Model (9) included the interaction terms of farmers’ cognition of fertilizer hazards and the attribute variables, and the regression coefficients of the three interaction terms were not significant, indicating that farmers’ cognition of fertilizer hazards did not significantly affect the WTA value for the fertilizer reduction program. However, the regression coefficients were insignificantly positive, and the WTA value for the program may be reduced to some extent.
5.3. Individual Characteristics and Preferences for Fertilizer Reduction Programs
The latent-class logit model was used to analyze the heterogeneity of the sampled farmers in terms of categories. The latent-class logit model can divide the sampled farmers into several potential categories and estimate the coefficients of the attribute variables for each category. In order to test the suitability of the number of sample categories, the Bayesian information criterion (BIC) and Akaike information criterion (AIC) are used to measure the model’s goodness-of-fit. The Bayesian information criterion can effectively prevent the model from overfitting and is the main basis for the latent classification. We classified the sampled farmers into from two to five categories, with AIC values of 942.8227, 886.0922, 849.9880, and 855.4841, and BIC values of 969.4251, 927.4738, 906.1487, and 926.4239. When the farmers were divided into four categories, the AIC value and BIC value were the lowest, and the goodness-of-fit of the model was the best. Therefore, we divided the farmers into four categories and regressed. The estimated results of the latent-class logit model are shown in
Table 10. The basic characteristics of the four types of farmers were further statistically analyzed, and the results are shown in
Table 11. The latent-class logit regression results and the descriptive statistics of the four types of farmers were used to analyze the main characteristics of each category.
There were 45 households in the first category, accounting for 31.69% of the total sample. The regression results showed that, among the attribute variables of the fertilizer reduction program, only the contract length and the annual cash subsidy coefficients were significant, and they were both positive. According to the farmer characteristic statistics, the non-agricultural income for farmers in the first category is in third place, but is close to the sample mean. These farmers also own a certain amount of agricultural land (second place) and agricultural income (second place), the proportion of non-agricultural income is relatively high, and the proportion of agricultural income is relatively low. Therefore, this category can be referred to as the “farm-oriented hybrid type”. This category of farmers is more dependent on labor income, but also has a certain dependence on agricultural income, so these farmers prefer a higher subsidy amount and a long-term fertilizer reduction program.
There were 14 households in the second category, accounting for 9.86% of the total sample. The regression results showed that none of the four coefficients of attribute variables were significant, that is, the farmers were not interested in the fertilizer reduction program. This type of farmer has the lowest mean non-agricultural income, the highest agricultural income, the most farmland area, and the lowest level of education, so this category can be called “ farm type”. These farmers mainly rely on agricultural income and are not interested in participating in the fertilizer reduction program, which may imply that they are not willing to participate in any form of fertilizer reduction program.
There were 18 households in the third category, accounting for 12.68% of the total sample. The regression results showed that only the coefficient of the proportion of participating land area was significant and positive, that is, the larger the proportion of participating land area, the higher the willingness of the farmers to participate in the fertilizer reduction program. According to the characteristics of the third category of farmers, the non-agricultural income is in second place, the agricultural income is the lowest, and the agricultural land area of the farmers is generally the smallest, implying that they may prefer to put more land into the fertilizer reduction program. Hence, they can be called “off-farm workers”. This category of farmer has the smallest land area, mainly relies on labor income, and prefers a fertilizer reduction program with a larger proportion of land area.
There were 65 households in the fourth category, accounting for 45.77% of the total sample, which was the highest proportion among the four categories. The regression results showed that the coefficients of the proportion of participating land area and the proportion of reducing fertilizer were significantly negative; the coefficient of contract length was negative, but not significant; and the coefficient of annual cash subsidy was positive and not significant. In terms of the characteristics of this category, the level of education is the highest, the household non-agricultural income is the highest, the household agricultural income is relatively low (ranking third), the agricultural land area is relatively small (ranking third), and the per ha fertilizer consumption is in the middle level of the sampled farmers (ranking third). Therefore, this category can be called “off-farm-oriented hybrid type”. This category of farmer mainly earns income from non-agricultural work, but also has certain income from agriculture, so they are not willing to participate in the fertilizer reduction program.
In general, the “farm type” farmers seemed to be uninterested in participating in the fertilizer reduction program; the “off-farm workers” had the lowest agricultural income and preferred to put more land into the fertilizer reduction program. The “farm-oriented hybrid type” and “off-farm-oriented hybrid type” farmers accounted for the highest proportion, and their regression results were the closest in terms of the benchmark regression results, indicating that these two types of farmers are the most typical.