As the concentratareas of carbon emissions, cities are meanwhile the most suitable to deeply participate and effectively implement nationwide emission reduction targets and emission reduction policies numerous city administrations that have strong jurisdictions. ased on the data of 108 cities in the Yangtze River Economic Beltestimate the carbon intensity of energy consumption and analye its evaluation of spatial-temporal features and influencing mechanisms kernel density analysis, STIRPAT model, and quantile regression panel model. We fthat: (1) The cities’ average carbon intensity energy in the Yangtze River Economic Belt hamaintained a downward trend from 2003 to 2020, decreasing from 1.87 yuan in 2003 to 0.51 t/ 10 000 yuan in 2020 with an annual average growth rate of -7.34%. (2) For the spatial evaluation, the cities’ average carbon intensity energy showa continuously increasing spatial feature, which shiftfrom downstream regions with scarce fossil fuels to midstream and upstream regions with abundant fossil fuels. Besides, the cities’ carbon intensity energy s highest in small cities, followed by medium-sized cities, and lowest in large cities “small cities>medium-sized cities>large cities”. (3) Economic growth, population size, government intervention, and industrial structure significantly can reduce the carbon intensity of energy consumption, but the human capital level and energy intensity were the main contributing factorshe positive promotion effect of human capital level no significant for low energy intensity cities. e suggested that government take different carbon intensity levels into full consideration before policy making.

Harmonious urban-water relations are the basis and key to promoting the high-quality development of new urbanization and realizing the harmonious coexistence of human beings and nature. Using exploratory spatial data analysis methods and the GTWR model to construct a cross-scale spatial heterogeneity analysis framework of basin-region-province-municipality, an empirical study was conducted in 110 prefecture-level cities in the Yangtze River Economic Belt from 2011 to 2020 to reveal the inner rules and spatial evolution patterns of the coordinated development of urbanization efficiency and water resources efficiency in the Yangtze River Economic Belt, and to analyze the spatial heterogeneity drivers of its coordinated development. The spatial heterogeneity of the factors driving the coordinated development was also analyzed. The results show that: (1) the urbanization efficiency and water resources efficiency of the Yangtze River Economic Belt as a whole show a fluctuating upward trend during the study period, with the urbanization efficiency fluctuating upward in the range of 0.395 0-0.767 6 and the water consumption of 10 000 Yuan GDP fluctuating downward in the range of 167.57 to 69.63. (2) the urbanization efficiency of the Yangtze River Economic Belt as a whole gradually increases from the upstream to the downstream areas, and the water resources efficiency shows a high level at both ends and a high level in the middle. (3) The coupling between urbanization efficiency and water resources efficiency in the Yangtze River Economic Belt is increasing, evolving from a loosely coupled state to a synergistic coupling state, and the overall coordinated development degree is on the verge of disorder, showing a spiral development process of rise-fall-rise again.(4) The intrinsic drivers of the coordinated development of urbanization efficiency and water resources efficiency in the Yangtze River Economic Belt are, in terms of absolute values of multi-year averages, mainly economic development level>agricultural development status>industrialization level>openness to the outside world>administrative level of cities>water resources endowment.

How to effectively play regional low-carbon emission reduction from the industrial policy level is a topic worthy of attention. What is the effect of green industrial policy as an important institutional driver to promote low-carbon transformation and upgrading on regional low-carbon emission reduction effect? What is the impact of the emotional signals released by the policy text on low-carbon emission reduction? Based on the content of green industrial policies introduced by the central government and the Yangtze River Delta region, this paper examinethe affective preferences of green industrial policies from 2010-2019, and systematically examine the low-carbon emission reduction effects of green industrial policies. that green industrial policies mostly designed based on welfare-oriented preferences The PMC model showthat the internal consistency of local green industrial policies s high and the scores he internal consistency of central green industrial policies s and the scores balanced. These findings suggest that the rational design of green industrial policy content s conducive to the promotion of low-carbon emission reduction effects, which in turn promote regional low-carbon emission reduction levels.

The spatial differences between regional carbon emission efficiency and high-quality development are significantt is to dynamically grasp the evolution of the spatial and temporal patterns and the coupling development trend of the two to promote the balanced development of the Yangtze River Delta urban agglomeration. Taking 26 prefecture-level cities as the research objects, this paper construct the index system and coupling degree model to measure carbon emission efficiency and high-quality developmententropy weight method, SBM model and Tobit panel model to reveal the macro pattern of carbon emission efficiency and high-quality developmentfrom 2000 to 2020 he spatial and temporal differences of their coupling development and the driving factors. The results as follows: (1) The inter-annual changes in carbon emission efficiency and high quality development show a synchronized upward trend fluctuation; he former gradually show “North-South” differences in macro patterns, while the latter’s “East-West” differences gradually becme. (2) The coupling degree and coupling coordination degree of carbon emission efficiency and high-quality development in the period of primary coordination and above, with an overall upward trendhe gap between the two in different provinces and cities haexperienced a change of first narrowing and then widening. The macro pattern of the two differsignificantly, but both show the characteristics of contraction of the high-value area to the core cities and the periphery, and the “southeast-northwest” difference is obvious. (3) Economic base, government intervention, urbanization rate and industrial structure contribute to the degree of coupling coordination, while carbon emission intensity haa negative constraining effect on it. Finally, suggestions put forward in terms of promoting synergy between the two, strengthening inter-regional interaction and differentiated development paths.

The traditional urbanization mode driven mainly by industrialization has brought about rapid growth in carbon emissions. The new-type urbanization mode, characterized with people-centered, urban-rural integration, industry and city interaction, energy-saving and intensive, and suitable living conditions, has the potential advantages of low-carbon development. ow to promote the coordinated development between the construction of new-type urbanization mode and low-carbon development is an important issue in the current low-carbon transformation of Chinese cities and towns. This article tthe Yangtze River Delta urban agglomeration as an example, and use the coupling coordination model and geographic detector detection method to explore the process and mechanism of the coupling coordination relationship between new-type urbanization and low-carbon development. The results indicate that: （1) The overall trend of the coupling and coordination between new-type urbanization and low-carbon development in the is improving, forming a spatial pattern with Shanghai as the core; （2) The evolution and types of the coupling and coordinated development of cities significantly different. The  of the eastern cities faster than that of the western cities. Relatively speaking, the lag of new-type urbanization development s the main obstacle to the improvement of the coupling and coordination level of new-type urbanization and low-carbon development; （3) There s significant interaction between new-type urbanization and low-carbon development. During the study period, there s a significant increase the influence of new-type urbanization system on low-carbon development system, while the influence of the subsystem of low development driving force and low-carbon pressure on new-type urbanization system reduced slightly.

Scientific understanding of urban-rural settlement evolution and rules of coal resource-based city is a crucial part in implementing the strategy of coordinated development of urban and rural as well as promoting the sustainable development of resource-based city. This paper takes Huaibei City, Anhui Province, as an example to explore the urban-rural settlement evolution and coordinated development mode of coral resource-based city since 1990s, from the dimensions of settlement space, economic structure and social life. The methods of field study, remote sensing data and statistical data analysis were adopted. The main findings were as follows: （1) Urban construction land scale increased rapidly. There were significant spatial differences between changed areas. However, rural construction land scale increased slowly, and the spatial agglomeration degree was gradually improved. The urban-rural settlement evolution underwent three stages, including low-speed dispersion, mono-centric agglomeration and poly-centric cluster development stages. (2) Urban industry structure shifted from coal leading industry to a diversified and non-coalized orientation. Rural agricultural industrialization was greatly promoted. The proportion of nonfarm income was significantly increased. Besides, the economic ties between urban and rural areas were gradually strengthened. The economic structure experienced dual development, transformation development and diversified development stages. （3) The system of “small society” in mine was gradually collapsed due to the social space integration of “mine-urban”; The rural social network opened gradually, and the residential mode differentiated. Social life evolution of urban and rural areas underwent three stages, including “urban-mine-rural” segregation stage, largely increased population mobility stage and lifestyle integration stage. (4) Based on the analysis of socioeconomic development of typical rural towns, this study reached to a conclusion of three coordinated development modes of rural and urban according to the difference of development characteristics. These stages included mine-town (village) collaboration mode, mine (town) joint development mode and leading effect of center town mode. This study made attempts to provide reference for the “urban-mine-rural” transformation development of similar urban-rural types.

With the increasing attention of urban development risk governance research, urban resilience has become a hot topic in this field. selecting 26 prefecture-level cities in the urban agglomeration in the middle reaches of the Yangtze River as the research objecton urban resilience indicators from 2006 to 2021 evolutionary resilience theory employ entropy evaluation method to measure the level of urban resiliencehe spatial-temporal evolution characteristics of urban resilience the main influence factors of urban resilience evolution panel data regression models finally the development trend of urban resilience  through gray prediction models. The research results indicate that: (1) During the research period, the overall level of urban resilience in the urban agglomeration in the middle reaches of the Yangtze River s relatively low, with a slow, sustained and steady growth trend in urban resilience significant regional differences, small dispersion, and a slow convergence trend; Regions with high resilience  while low resilience regions widely distributed, which basically formed the ‘core-periphery’ spatial pattern of regular triangular shape with high resilient zone at the apex, medium and low resilient zones at the edge.（2) The intensity of openness, diversification of industrial structure, and urbanization level the main influencfactors the development and change of urban resilience in the region, with urbanization level having the greatest impact, the intensity of openness and diversification of industrial structure having slightly weaker impact, and administrative and innovation capabilities having no significant impact on urban resilience. Therefore, it is necessary to further improve regional administrative management and innovation capabilities. (3) In the future, the urban resilience of the research area still exhibit a steady growth trend, and the overall level of urban resilience will be significantly improved. Regions with high resilience will significantly increase, and regional differences might be expanding. The core position of regional central cities such as Wuhan, Changsha, Nanchang, and Zhuzhou will be further， and the “core-periphery” urban resilience spatial pattern will be further stabilized.

The dynamic change of tourism eco-efficiency of urban agglomeration in the middle reaches of Yangtze River from 2011 to 2020 was measured and analyzed by combining DEA window analysis method with super-efficiency SBM model. The modified gravity model and social network analysis method were used to explore the evolution of spatial network structure of tourism eco-efficiency in urban agglomeration, and the panel data regression model was constructed to analyze its network effects. The results showed that: (1) The overall mean value of tourism eco-efficiency was 0.719, which showed a downward trend in fluctuation with significant regional differences. (2) The fluctuation of the overall network density and the number of network relations increased to 0.162 7 and 123, and the network level and network efficiency decreased slightly but basically stabilized at 0.2 and 0.8. It that the correlation network structure of tourism eco-efficiency in urban agglomeration  loose, with certain hierarchical characteristics, and the overflow channel was not smooth, which needed to be further optimized. (3) The individual indexes of Wuhan, Changsha, Nanchang, Xiaogan and Yichun were always higher than the mean value, which occupcentral position in the correlation network of tourism eco-efficiency. The individual indicators of Changde, Yiyang, Huanggang, Yichang, Xiangyang and other cities were always lower than the mean value and relatively lag behind, at the edge of the network. (4) The increase of overall network density, the decrease of network level and network efficiency could significantly improve the overall tourism eco-efficiency of urban agglomerations and narrow the regional differences. The improvement of individual network centrality index had a significant promoting effect on urban tourism eco-efficiency.

 The monitoring data in Shishou Section of the Yangtze River. There were 41 species in Cypriniformes, accounting for 66.1% of the total species. There was difference in the number of fish species between years, with 2021 seeing more species than any other year. The diversity index and multivariate statistical analysis showed that the Pielou evenness index and Simpson dominance index fluctuated the Shannon-Wiener index and Margalef richness index in 2021 were larger than those in other years. The ABC curve and W statistics showed that the fish community structure was relatively stable in 2018, and the dominant species were mainly large individual species. From 2020 to 2021, the fish community mainly consists of species with short life cycle and small individuals. In 2022, compared 2020-2021, the fish community was relatively stable. Cluster analysis and NMDS analysis showed that the fish community could be divided into two groups ANOSIM test showed no significant difference between the two groups. SIMPER’s analysis showed that the abundance of benthic omnivorous fish adapted to flowing water environment decreased from 2020 to 2022 compared with 2018，while  abundance of benthic carnivorous fish adapted to still water and semi-migratory fish increased， the fish community in Shishou section of the Yangtze River was stable, and the fish size tended to increase.

In order to study the reproductive status of fishes downstream of the Wudongde dam during its initial reservoir impoundment, and to optimize the implementation of ecological operation of Wudongde Hydropower Station, this study investigated the status of early-stage fish resources in the Wudongdedam section, in the lower reaches of the Jinsha River,during April 21 to June 8, 2021, and March 9 to June 28, 2022. The results showed that at least 3 orders, 7 families and 19 species (genera) of fish spawned and bred in the downstream of Wudongde Dam, including 10 species with drifting eggs and 6 endemic fishes in the upper reaches of the Yangtze River. In 2021 and 2022, the total quantity of fish eggs with drifting eggs was 321.72×104ind and 886.40×104ind, respectively. There were three spawning peaks during the survey period in 2021 and four during the survey period in 2022. Correlation analysis between the reproductive scale of fish with drifting eggs and the hydrological process showed that the reproductive scale of fish with drifting eggs in the Wudongde section was significantly correlated with peak flow (R=0.833, P<0.05), flow growth rate (R=0.810, P<0.05) and the accumulated flow increament (R=0.762, P<0.05).The flood peak can somewhat promote the reproduction of fish with drifting eggs. Compared with the previous survey results, our study indicated that the main species of spawning fish increased and the quantity of fish eggs decreased. Further study is suggested on the reproductive requirements of fish in the lower section of Wudongde dam. It is also encouraged to optimize the operation scheme of the dam to maximize the ecological benefits.

The implementation of industrial carbon peak action is an important task to achieve the carbon peaking and carbon neutrality (double carbon) goals.Improving industrial water efficiency is the key to achieve industrial low-carbon transformation, which lays a solid foundation for the carbon peak of industry in 2030. This paper fully considered the task of carbon peaking in the industrial field at the present stage, and incorporated water conservation and emission reduction into the input-output index system. By using the super-efficiency SBM model, Dagum Gini coefficient decomposition and convergence test, this study investigated industrial water efficiency and identified the gaps among  the upper, middle and lower reaches of provinces along the Yangtze River Economic Belt from 2010 to 2021. The empirical results showed that: (1) the overall improvement of industrial water use efficiency was slow, where  Zhejiang and Jiangsu provinces and Shanghai city ranked the top three. Sichuan and Yunnan Provinces showed relatively stable industrial water use efficiency.The industrial water use efficiency in Hunan Province was greatly improved.There existed a downward trend in industrial water use efficiency in Anhui and Guizhou provinces.(2) Both the overall gap and intra-regional gap of industrial water use efficiency demonstrated an enlarging trend.The upstream and downstream inter-regional gap and the fluctuation range were large. The interaction of inter-regional and intra-regional gap contributed the most to the overall gap. (3) The σ convergence characteristics of industrial water use efficiency were not obvious. The dispersion of different regions did not show a consistent decline over time. There was only a brief convergence trend from 2020 to 2021, but there were typical absolute β convergence and conditional convergence characteristics in the whole sampling period. This implied that the regions with low industrial water use efficiency were struggling to narrow the gap between those with high industrial water use efficiency. Meanwhile, the industrial water use efficiency of each region would eventually converge to its own steady-state level. Therefore, in order to improve industrial water use efficiency in the Yangtze River Economic Belt, it is necessary to further enhance the ability of industrial technology innovation, accelerate the regional integration construction, and promote the industrial green development.

The core of ecological compensation for inter-basin water transfer project is to coordinate the interest relationship of various stakeholders.In order to ensure the reasonableness and implementability of compensation and promote the sustainable performance of project benefits, from the perspective of cooperative game, the middle Route of South-to-North Water Diversion (SNWD) project was taken as an example in this study. A horizontal ecological compensation standard of inter-basin water transfer project was constructed,based on the Shapley value method by sorting out the response relationship of water use and sewage emissions between the water receiving area, water source area and downstream affected area during the operation of the project and sharing the cost and benefit among the subjects . The results showed that: (1) from 2019 to 2021, the annual ecological compensation required to be paid by the water receiving area was 5.526 billion yuan, 5.483 billion yuan and 6.396 billion yuan respectively. The compensation amount required to be paid by Beijing and Tianjin was higher than other cities; (2)Under the joint influence of the optimal allocation of transferred water and the level of economic development, the value of water resources in Beijing, Hebei and Henan provinces tended to increase, while that in Tianjin city tended to decrease;(3)The compensation funds received by the source and downstream areas made up for some of their own water transfer losses.The funds increased with the increasing benefits of water transfer in the receiving area, which helped to further improve the stability of cooperation. Overall, this study revealed the mechanism of optimal allocation of water resources in inter-basin water transfer projects from the perspective of cooperation, and provided a feasible idea for the determination of the balancing mechanism of project benefits and the accounting of compensation standards.

Algal-derived off-flavors (ADOFs) pollution is a global issue that severely impacts water quality and the aesthetic value of water bodies. To investigate the temporal and spatial changes、sources and driving factors of this pollution in Lake Yuehu, a city landscape lake, we carried out one year-long study from March 2021 to February 2022. Monthly sampling was performed, and the concentrations of four major extracellular off-flavors (2-methylisoborneol (2-MIB), geosmin (GSM),β-cyclocitral and β-ionone), as well as water quality indicators and algal composition were analyzed. Using Spearman's correlation analysis, redundancy analysis, and partial least squares structural equation modeling (PLS-SEM), we examined the sources and drivers of off-flavors. Our results indicated that the nearshore area of Lake Yuehu was heavily polluted by 2-MIB and GSM, with annual average concentrations as high as 35.04 ng/L and 8.58 ng/L, respectively. High concentration pollution mainly occurred from June to September. The 2-MIB likely originated mainly from Oscillatoria, Pseudanabaena, and Anabaena, while GSM was potentially derived from Anabaena, Aphanizomenon, Pseudanabaena, and Oscillatoria. β-cyclocitral and β-ionone did not appear to be primarily algal in origin. PLS-SEM results revealed that meteorological conditions (radiation and evaporation) served as the principal drivers of off-flavors in Lake Yuehu, while nutrient levels (total nitrogen and total phosphorus) had relatively weak impacts. Meteorological conditions not only had a significant direct influence on algal growth, but also significantly positively moderated the relationship between nutrients and algae. However, they did not moderate the relationship between algae and off-flavors. These findings implied that meteorological conditions mainly induced off-flavors via influencing algal growth, rather than directly affecting the production of off-flavors.

Since the start of operation of the Three Gorges Reservoir impoundment and power generation, the hydrological conditions in the tributary bays have undergone significant changes.Water quality in the bays is closely related to the hydrological conditions. This study took the Xiangxi Bay as the research area to explore the impact of hydrological conditions on water quality. A two-dimensional hydrodynamic water quality model of Xiangxi Baywas built using MIKE21. The changes in water quality was simulated in the bay caused by changes in water level in the main stream and upstream flow of Xiangxi Bay under various operating conditions. The results showed that an increase in the amplitude of water level variation in the main stream and an increase in the upstream flow rate of the Xiangxi Bay would significantly improve the hydrodynamic conditions of the bay. The amplitude of water level variation might increase by 1m, and the average flow rate might increase by 0.22 cm/s. The upstream flow rate might increase by 100 m3/s, and the average flow rate might increase by 0.5 cm/s; An increase in the fluctuation of the water level in the main stream would cause nutrients to migrate upstream, while an increase in the upstream flow of tributaries would accelerate the migration of nutrients to the estuary. Additionally, the enhancement of hydrological conditions in the above two cases would reduce the concentration of chlorophyll-a in the middle section of the bay; Finally, according to the distribution of chlorophyll-a concentration under comprehensive hydrological conditions, thewater area which is 0-17.1 km away from the estuary was divided into level-I water bloom sensitive area, which implied that the monitoring of this reach should be strengthened during the period of high possibility of frequent water blooms. This research results may provide theoretical support for the enhancement of hydrodynamic power, the improvement of water quality, the rational prevention and control of water bloom, and the ecological regulation in the tributaries of the Three Gorges Reservoir.

Taking the Yangtze River Economic Belt as the research object, this paper discussed the relationship between agricultural carbon emissions and per capita net income of rural residents based on the characteristics of the current situation of agricultural carbon emissions in each province from 2005 to 2020. The driving factors of agricultural carbon emissions were also analyzed, and the agricultural carbon emissions of each province from 2023 to 2030 were projected. The results showed that:（1）The agricultural carbon emissions showed significant changes in either growth or decline. During the sampling period, 5 provinces were in the trend of growth, of which Yunnan province was the largest. Six provinces were in the trend of decline, of which Shanghai was the largest. Agricultural carbon emissions were generally "high in the middle and low in the east and west".（2）The EKC test results showed that the EKC relationship existed in all provinces of the Yangtze River Economic Belt, however the overall trend was not completely consistent. Based on the difference of EKC curves, 11 provinces could be divided into three types: inverted "N" type, inverted "U" type and positive "U" type.（3）The results of STIRPAT showed that the number of people employed in the primary industry and the gross output value of agriculture, forestry, animal husbandry and fishery dominated agricultural carbon emissions, while the level of agricultural technology had a suppressive and less powerful effect.（4）The GM (1,1) projections showed that the agricultural carbon emissions showed obvious fluctuations.The effectiveness in reducing emissions from 2023 to 2030 varied from province to province. Based on the forecast results, the 11 provinces could be divided into three types: "continuously declining", "continuously growing" and "rising then falling".

For the strong agricultural country building, the primary task is to ensure national food security.The production efficiency of peasant households is highly related to the stability of national grain output and grain strategic security. Based on the survey data of rice farmers in Jiangsu province, the production efficiency of farmerswas measuredbased on the translogarithmic stochastic frontier production function modelin this paper. Furthermore, the influence of the development of agricultural industry chain on the production efficiency of farmers was investigated from the perspective of agricultural productive service supply. The influencing mechanism of the scale of production investment, production investment structure and formal credit access of farmers on production efficiency of farmerswere also explored.It was found that the average productivity of the sample farmerswasas highas 0.87,although there was still a productivity loss of 13%. The service of agricultural industry chain had significant promoting effects on the production efficiency,especiallyfor small-scale farmers. By encouraging farmers to increase the scale of land operation and agricultural machinery input per unit of land, agricultural industry chain services promoted farmers to replace labor with machinery input, expanded the scale of land operation instead of purchasing agricultural machinery, and improved farmers' access to formal credit to promote their production efficiency. Policy suggestions were put forward to strengthen the service level of the agricultural industry chain, to adjust the supply of agricultural industry chain services to different conditions, and to innovate the supply of rural financial products.

The purpose of this study was to analyze the impacts of comprehensive land consolidation on farmers’ livelihood resilience and to explore a new path for sustainable development of farmers’ livelihoods. Based on the analysis of the influencing mechanism of comprehensive land consolidation on farmers’ livelihood resilience, 466 household survey data from Jingmen City, Qianjiang City, Jingzhou City and Xiantao City in Hubei Province were taken as samples. The analysis was performed using a difference-in-differences model. The baseline regression results showed that comprehensive land consolidation significantly improved the buffering ability, self-organizing ability, learning ability and livelihood resilience of farmers by 1.8%, 2.8%, 1.1% and 5.7%. After a series of robustness tests, the regression results were still valid. The results of heterogeneity analysis showed that the impact of comprehensive land consolidation on farmers’ livelihood resilience decreased with the increase of percentile of farmers’ livelihood resilience, and also affected by differences in landform types. Comprehensive land consolidation in plain areas had a greater effect on the farmers’ livelihood resilience than in hilly and mountainous areas.In conclusion, the comprehensive land consolidation should take the sustainable development of human beings as the core concept, and enhance the efficiency of land by comprehensive management of all factors in the whole area. In addition, it is necessary to fully consider the differences between farmers. Governments should invest more targeted projects to help the most vulnerable farmers to reduce their livelihood vulnerabilities, and establish more opportunities to enhance the endogenous development power of farmers under the impact of risks. Furthermore, combined with resource endowment and environmental characteristics, comprehensive land consolidation for ecological sustainability should be carried out according to local advantages. All these are conducive to the realization of the resilient development of farmers’ livelihoods and the continuous protection of farmers’ livelihood security and quality.

Grey water footprint (GWF) can be used as an indicator to quantitatively evaluate the impact of agricultural production process on water environment. Traditional methods for calculating GWF of crop production are usually based on rough estimation, which cannot accurately quantify the difference between GWFs caused by fertilizer types. In this study, field experiments were conducted to observe the migration process of water and fertilizer in paddy field under conventional fertilization (CF), biochar-based fertilizer (BF) and organic fertilizer partially replacing chemical fertilizer (OF) treatments. This study aimed to identify the most critical pollutants for each runoff event and to explore the changes in the GWF under different fertilization types. The results showed that: (1) Compared with CF treatment the BF treatment and OF treatment both reduced the pollutant losses of total nitrogen, nitrate nitrogen, ammonium nitrogen and total phosphorus, and increased rice yields by 3.48% and 6.81%, respectively. (2) In all runoff events produced during the study period, the occurring frequency of total nitrogen as the key pollutant was the highest, while total phosphorus was also critical pollutants under certain circumstances, implying that the key pollutant determining the amount of grey water was variant. (3) The impacts of fertilizer typeson the GWF of rice production were in the order of: CF (622.85 m³/t) > BF (367.68 m³/t) > OF (318.18 m³/t). OF reduced the grey water footprint of rice production and increased the rice yield, which help achieve the goal of pollution reduction and yield increase. The method for calculating GWF of crop production based on field measurements can accurately quantify the influence of different fertilizer types on grey water footprint.