1 km

Scientific Data volume 10, Article number: 363 (2023) Cite this article

1 Altmetric

Metrics details

Crop-specific, high-resolution phosphorus rate information is essential for sustainable agricultural fertilizer management in China. However, substantial uncertainties exist in the current phosphorus fertilizer dataset because of only coarse national statistics used in dataset development and no crop-specific information provided. This study harmonized provincial and county-level phosphorus and component fertilizer statistics and crop distribution data to generate 1 km gridded maps of phosphorus rate for rice, wheat and maize in the years of 2004–2016 (CN-P). CN-P provides a comparable estimate on phosphorus rate for each crop over 2004–2016, and demonstrates an improved spatial heterogeneity. Existing dataset developed using national statistics tends to smooth out the variability within country and significantly underestimates actual phosphorus rate. CN-P shows that, during 2004–2016, wheat received the most phosphorus rate (8.7 g P2O5 m−2), while maize showed the rapidest increasing trend (2.36% yr−1). The CN-P dataset has the potential to be widely applied in modeling studies on sustainable agricultural fertilizer management strategies and phosphorus pollution.

Phosphorus plays a vital role in growing crops1,2 and meeting the food demands3,4,5. China is the largest producing and consuming country of phosphorus fertilizer6. In 2020, there is around 13 million tons phosphorus fertilizers produced in China, and approximately 9 million tons were used in agriculture, which accounts for 20.6% of the global phosphorus fertilizer consumption7. Since 1980, the total phosphorus fertilizer application has grown, and substantial phosphorus fertilizer demands were projected in 2050 due to growing populations in China4.

Overuse of phosphorus fertilizer results in a series of environmental problems8. Currently, there are twenty provinces are subjecting to different level phosphorus pollution9. Soil available phosphorus was found increased from 17.09 mg L−1 in the 1990s to 33.28 mg L−1 in the 2000s in China10; and 48 Tg phosphorus has leached to water bodies in the past 60 years3,11. As a consequence, eutrophication was found in 23.6% major lakes in China, and 4.5% were moderately eutrophic, and 0.9% were severely eutrophic12. Therefore, sustainable phosphorus fertilizer management is crucial to both food security and environmental conservation.

To investigate the efficient phosphorus fertilizer management, it is critical to understand the state of agricultural phosphorus rate. Previous studies have been examining the historical phosphorus footprint6, phosphorus budge1 and phosphorus losses3. These analyses are typically based on the national-scale phosphorus statistics from FAOSTAT7 and IFASTAT13. A gridded agricultural phosphorus dataset was developed by Lu and Tian14, in which they obtained the gridded phosphorus consumption by multiplying the national-scale phosphorus rate statistics with gridded cropland data. Despite widely used, there are three main deficiencies in the existing dataset: (1) The most-widely available phosphorus statistic is national-scale. This is a very coarse data, and tends to lump the spatial heterogeneities within country; (2) Data source of phosphorus rate was from statistics department and cropland data was obtained from remote sensing method in the data product. Direct composite of the two kinds of data may raise the inconsistency issue of statistical caliber; (3) The current gridded phosphorus dataset only provides the total agricultural phosphorus information, rather than crop-specific. The lack of crop-specific data leads to the difficulty to separate phosphorus fertilizer by various crops and growing seasons in modeling analysis.

To overcome these deficiencies, we constructed a new gridded phosphorus rate dataset, which was referred to as CN-P. CN-P is 1 km-resolution gridded dataset of phosphorus rate for three staple food crops (i.e. rice, wheat and maize) in China over 2004–2016. This contains gridded maps of phosphorus rate in each year by each crop. In data construction, we employed the county-level statistic, which is the finest agricultural statistics in China. After validation, we found that CN-P demonstrates a comparable crop-specific phosphorus rate with previous farmer's surveys, and provides an improved spatially-explicit distribution relative to previous state-of-the-art dataset when compared for individual year. We found that only using national statistics in previous gridded phosphorus dataset tends to smooth out the spatial variability and significantly underestimates agricultural phosphorus rate.

Multiple data were used to construct CN-P (Table 1). Statistical yearbooks, "Cost-benefit statistics of agricultural products", provide the provincial historical data on crop-specific phosphorus and component fertilizer rate over 2004–201615,16,17,18,19,20,21,22,23,24,25,26,27. Historical provincial crop-specific growing areas, total agricultural growing areas and total phosphorus and component fertilizer consumptions were downloaded from National Bureau of Statistics of China28. The above provincial data includes 31 provinces in the mainland of China, excluding Hongkong, Macau and Taiwan. The county-level statistics were based on "China County Statistical Yearbook"29,30,31,32,33,34,35,36,37,38,39,40,41, compiled by Chinese Academy of Agricultural Sciences. Note that the county-level phosphorus and component fertilizer consumption statistics only report the total fertilizer consumption, but not crop-specific. The county-level growing areas include both statistics of crop-specific and total agricultural growing areas in 2267 counties. The crop distribution data are based on 1 km gridded crop growing area distribution maps for rice, wheat and maize between 2004 and 2016 obtained from National Science & Technology Infrastructure42.

Figure 1 illustrates the steps of data processing, in which we harmonized above provincial/county-level phosphorus and component fertilizer statistics and crop distribution data to generate 1 km gridded maps of phosphorus rate for rice, wheat and maize in the years of 2004–2016.

Diagram for processing spatially explicit time series of phosphorus rate of rice, wheat and maize in China over 2004–2016. "P" refers to phosphorus fertilizers and "Com" refers to Component fertilizers.

We firstly multiplied the provincial phosphorus and component fertilizer rate of rice, wheat and maize with their growing areas to generate crop-specific phosphorus and component fertilizer consumption at each province in each year. After that, subtracting above three crops from the total phosphorus and component fertilizer consumption at each province can produce the provincial phosphorus and component fertilizer consumption of other crops. We then divided the phosphorus and component fertilizer consumption of other crops by growing areas of other crops to obtain the phosphorus and component fertilizer rate of other crops at each province in each year. Following above steps, we generated the phosphorus and component fertilizer rate of rice, wheat, maize and other crops at each province over 2004–2016.

Using the county statistics, we firstly calculated the growing area of rice, wheat, maize and other crops in each county-year combination using the same method for provincial data in the above step. For the missing data, cubic spine interpolation method was employed to gap-fill when data were unavailable in less than 3 consecutive years. To harmonize the provincial and county-level statistics, we did not simply multiply the two statistics directly, as did in previous dataset14. Instead, we calculated the ratio of crop-specific phosphorus/component fertilizer consumption to county total as Eq. (1). Then, the phosphorus/component fertilizer consumption for crop-specific in each county-year combination can be estimated by multiplying the county-level total phosphorus/component fertilizer consumption statistics with above ratios as Eq. (2). Finally, these fertilizer consumptions divided by the crop-specific county-level growing area statistics can produce the phosphorus /component fertilizer rate of rice, wheat and maize at each county in each year as Eq. (3). This practice enables us avoiding the potential inconsistency problem of statistical calibers compared to the direct computation method in early study14, and maintaining the spatial heterogeneities between counties. Further comparison can be seen in the technical validation section.

where FRi, p, y is phosphorus /component fertilizer rate for crop type i, at province p, in year y; Ai, j, y is growing area statistics for crop type i, at county j located in the province p, in year y; Rati, j, y is the ratio of phosphorus /component fertilizer consumption for crop type i to the summed fertilizer consumption \(\left({\sum }_{i}{{\rm{FR}}}_{i,p,y}\times {{\rm{A}}}_{i,j,y}\right)\), at county j located in province p, in year y. FCj, y is the phosphorus/component fertilizer consumption statistics at county j, in year y; FCi, j, y is the phosphorus/component fertilizer consumption estimates for crop type i, at county j, in year y. FRi, j, y is phosphorus /component fertilizer rate for crop type i, at county j, in year y

In the step, we transferred the above phosphorus /component fertilizer rate to phosphorus rate (i.e. P2O5 g/m2). The raw statistics of phosphorus fertilizer are already the amount of P2O5 from phosphorus fertilizer. For the data of component fertilizer, we converted it to grams of P2O5 by multiplying the weighted P2O5 content of imported and domestic-produced component fertilizer each year (Table 2). These steps allow us estimating the crop-specific phosphorus rate at county-level in the period of 2004–2016 by summing up the gram of P2O5 per unit of crop growing areas from both phosphorus and component fertilizers each year.

As the final step, gridded crop growing area distribution maps of rice, wheat and maize were used as the mask raster, and extracted the pixels with the crop growing to produce the gridded maps of phosphorus rate for each crop in each year.

The CN-P43 dataset is publicly available for download from the Zenodo repository. The dataset is saved in the form of GeoTiff format. It is organized in folders according to crops, and named by the format of "CNP_<crop>_<year>.tif" (Table 3). Each file contains a 1 km gridded phosphorus rate map of rice, wheat and maize of China in certain year, with the unit grams of P2O5 per unit growing area (i.e. g P2O5 m−2). These files were saved with the WGS84 geographic coordinate system. As an example of CN-P, Fig. 2 demonstrates the spatial distribution of phosphorus rate for rice, wheat and maize averaged over 2004–2016. Table 4 gives the national average and time trend for each crop calculated based on CN-P. Wheat received the most phosphorus rate (8.7 g P2O5 m−2, p < 0.05), and maize showed the rapidest increasing trend (2.36% yr−1, p < 0.001) (Table 4).

Spatial distribution of the crop-specific phosphorus rate averaged over 2004–2016. (a) map of rice; (b) map of wheat; (c) map of maize.

To test the validity of CN-P, we compared the data with the estimates of farmer's surveys44,45,46,47,48,49,50. As these estimates were conducted by questionnaire approach in individual years (denoted by points in Fig. 3) or average value over a period (denoted by collinear points in Fig. 3), these surveys could be viewed as an independent estimate. There was a large discrepancy between surveys, such as Du et al.45 estimated a much higher phosphorus rate especially for wheat (Fig. 3b) and maize (Fig. 3c) than others. In generally, the phosphorus rate of the three crops we estimated fell within the range of these studies (Fig. 3). For time trends, the farmer's surveys also presented a more significant increasing trends for maize (Fig. 3c) compared with rice (Fig. 3a) and wheat (Fig. 3b), which is consistent with the trends of CN-P. This part of comparison suggests CN-P is comparable with these independent farmer's surveys on phosphorus rate for each crop.

Comparison of crop-specific phosphorus rate with farmer's surveys over 2004–2016. Blue line is the estimate of CN-P; points are estimates in individual years, and the collinear points are the estimates averaged by multiple years in farmer's surveys.

Because there is no previous crop-specific phosphorus dataset, we calculated the total phosphorus consumption of China and spatial distribution of total phosphorus rate per cropland based on county statistics, which were used to compared with previous phosphorus datasets. In this analysis, as our county statistics include the total phosphorus and cropland data, our results based on county statistics should be viewed as observations.

Firstly, we calculated the total phosphorus consumption based on county statistics and compared with the results of FAOSTAT17, IFASTAT13 and NBS28 in Fig. 4. The phosphorus fertilizer results varied between data sources, which reflects uncertainties in phosphorus fertilizer statistics even for these official data sources (Fig. 4). This is because of different survey method used by these databases. For example, IFASTAT database is based on the survey sent to country correspondents, including fertilizer associations, fertilizer companies, consultants, experts etc. But the FAOSTAT database is based on the FAOSTAT fertilizers questionnaire. In generally, our estimates were within the range of the three data sources before 2013. Beyond 2013, our estimates were higher than the other three databases, but close to the results of NBS (Fig. 4). Based on our results, there was an increase in the phosphorus consumption over time until 2015 and then followed by a slight drop after the year. The time trend pattern of our estimates is consistent with NBS (Fig. 4). The year of the changing point is the same with the year when China's Ministry of Agriculture introduced the Action that seek to achieve zero growth in the use of chemical fertilizer in 201551.

Comparison of phosphorus consumption derived from CN-P, IFASTAT, FAOSTAT and NBS over 2004–2016.

Secondly, we further compared the spatial distribution map of total agricultural phosphorus rate per cropland with Lu and Tian14 averaged over 2004–2013 (Note: the latest year reported in Lu and Tian14 is 2013). The definition of cropland here is different from crop growing area used in CN-P. Cropland could be used to grow more than one crop seasons in the same field. Using our county statistics, we calculated the total agricultural phosphorus rate per cropland. For Lu and Tian's14 data, we converted grams of phosphorus in their dataset to grams of P2O5 by multiplying by the ratio of 142/62. We found that the spatial heterogeneity differed between the two datasets (Fig. 5). County statistic is higher than the results estimated by Lu and Tian14. Using county statistics, some regions reached to 18–26 g P2O5 per square meter cropland (Fig. 5b), while most areas were less than 18 g P2O5/m2 in the dataset of Lu and Tian14 (Fig. 5a). Another disagreement is that there was a high phosphorus rate hotspot in the central of China presented by our county statistics (Fig. 5b), while this was not shown in Lu and Tian14 (Fig. 5a). This reflects the data construction approach of Lu and Tian14, which multiplied a national-level phosphorus rate with gridded cropland area and adjusted with the national data of IFASTAT inventory. Such approach will lead to the underestimation of phosphorus rate because the national level phosphorus rate could smooth out the variability in the sub-national scale and underestimate the hotspot within country. Therefore, the presence of spatial heterogeneity of Lu and Tian14 is majorly from the gridded cropland areas, and cannot reflect the actual spatial heterogeneities. This part of analysis suggests that introducing county-level statistics is critically needed to show an improved estimates and spatial distribution in phosphorus rate.

Comparison of average phosphorus rate per cropland over 2004–2013. (a) maps of Lu and Tian14; (b) maps based on county statistics. To compare them, the same mask raster was used.

The uncertainties of CN-P are mainly from the following aspects: (1) We applied the provincial crop-specific phosphorus consumption ratio at county scales. Currently, there is no available statistics on county-level phosphorus consumption for each crop in China, and it is unlikely that such data will ever become available in a near future. (2) For conversing component fertilizer to phosphorus contents, we used a static phosphorus content in each year for the whole China (Table 2). Map accuracy can be further improved if the crop-specific county and phosphorus ratio of component fertilizer information are available. (3) We note that data construction before the year of 2004 is quite challenge as the "Cost-benefit Statistics of Agricultural Products" do not record crop-specific phosphorus fertilizer before 2004 and only document the data of total fertilizer (i.e. the sum of nitrogen, phosphorus and potassium fertilizers). If we could collect longer-term datasets, or more accurate approaches to gap-fill our data, we will provide future improvements. Therefore, a continuous survey of crop-specific fertilizer data and development of dynamic crop-type maps to meet the needs of current study are imminently required. The longer time series of crop-specific phosphorus rate maps will improve the characterization of geospatial and temporal patterns of phosphorus fertilizer management in China.

The code of CN-P is archived at the Zenodo repository: https://doi.org/10.5281/zenodo.7460564.

Zou, T., Zhang, X. & Davidson, E. Global trends of cropland phosphorus use and sustainability challenges. Nature. 611, 81–87 (2022).

Article ADS CAS PubMed Google Scholar

Zhang, J. et al. Supply-side drivers of phosphorus emissions from phosphorus supply chains in China. Ecosyst. Health Sust. 8, 2107578 (2022).

Article Google Scholar

Bai, Z. et al. Changes in phosphorus use and losses in the food chain of China during 1950–2010 and forecasts for 2030. Nutr. Cycling Agroecosyst. 104, 361–372 (2016).

Article Google Scholar

Mogollon, J. M., Beusen, A. H. W., Grinsven, H. J. M., Westhoek, H. & Bouwman, A. F. Future agricultural phosphorus demand according to the shared socioeconomic pathways. Global Environ. Change. 50, 149–163 (2018).

Article Google Scholar

Gong, H. et al. Toward the sustainable use of mineral phosphorus fertilizers for crop production in China: From primary resource demand to final agricultural use. Sci. Total Environ. 804, 150183 (2022).

Article ADS CAS PubMed Google Scholar

Jiang, S. et al. Phosphorus footprint in China over the 1961–2050 period: Historical perspective and future prospect. Sci. Total Environ. 650, 687–695 (2019).

Article ADS CAS PubMed Google Scholar

FAOSTAT. Food and Agriculture Organization Statistics. http://www.fao.org/faostat/en/ (2023).

Mekonnen, M. & Hoekstra, A. Y. Global anthropogenic phosphorus loads to freshwater and associated grey water footprints and water pollution levels: A high-resolution global study. Water Resour. Res. 54, 345–358 (2018).

Article ADS CAS Google Scholar

Liu, Q. Spatial-temporal changes of fertilization environmental risks of China. J. Agro-Environ. Sci. 36, 1247–1253 (2017).

Google Scholar

Ma, J. et al. Temporal and spatial changes in soil available phosphorus in China (1990-2012). Field Crop Res. 192, 13–20 (2016).

Article Google Scholar

Huang, J., Xu, C., Ridoutt, B. G., Wang, X. & Ren, P. Nitrogen and phosphorus losses and eutrophication potential associated with fertilizer application to cropland in China. J. Clean. Prod. 159, 171–179 (2017).

Article Google Scholar

Report on the State of the Ecology and Environment in China. (Ministry of Ecology and Environment, the People's Repoblic of China, 2022).

IFASTAT. International Fertilizer Association Statistics. https://www.ifastat.org/ (2023).

Lu, C. & Tian, H. Global nitrogen and phosphorus fertilizer use for agricultural production in the past half century: shifted hot spot and nutrient imbalance. Earth Syst. Sci. Data 9, 181–192 (2017).

Article ADS Google Scholar

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2005 (China Statistics Press, 2005).

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2006 (China Statistics Press, 2006).

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2007 (China Statistics Press, 2007).

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2008 (China Statistics Press, 2008).

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2009 (China Statistics Press, 2009).

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2010 (China Statistics Press, 2010).

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2011 (China Statistics Press, 2011).

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2012 (China Statistics Press, 2012).

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2013 (China Statistics Press, 2013).

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2014 (China Statistics Press, 2014).

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2015 (China Statistics Press, 2015).

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2016 (China Statistics Press, 2016).

National Development and Reform Commission of China (Eds.): Cost-benefit Statistics of Agricultural Products-2017 (China Statistics Press, 2017).

National Bureau of Statistics of China.: https://data.stats.gov.cn/easyquery.htm?cn=E0103 (2023).

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2005 (China Statistics Press, 2005)

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2006 (China Statistics Press, 2006)

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2007 (China Statistics Press, 2007)

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2008 (China Statistics Press, 2008)

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2009 (China Statistics Press, 2009)

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2010 (China Statistics Press, 2010)

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2011 (China Statistics Press, 2011)

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2012 (China Statistics Press, 2012)

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2013 (China Statistics Press, 2013)

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2014 (China Statistics Press, 2014)

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2015 (China Statistics Press, 2015)

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2016 (China Statistics Press, 2016)

Rural Socioeconomic Survey Team (Eds): China County Statistical Yearbook-2017 (China Statistics Press, 2017).

Luo, Y. & Zhang, Z. Crop planting distribution of three staple food crops in China over 2000–2009. National Science & Technology Infrastructure https://doi.org/10.12199/nesdc.ecodb.rs.2022.016 (2022).

Article Google Scholar

Zhang, T., Zhang, W. & Yang, X. CN-P (version 1): a crop-specific, 1 km-resolution phosphorus rate data product in China over 2004–2016 (Version 1). Zenodo https://doi.org/10.5281/zenodo.7460564 (2023).

Zhang, W. Yi, J., & Zhang, F. Reports of Fertilizer Development in China-2016. (China Agricultural University Press, 2017).

Du, J. et al. Spatial-temporal evolution characteristics and driving factors of partial phosphorus productivity in major grain crops in China. J. Plant Nutr. Fert. 28, 191–204 (2022).

Google Scholar

Li, H., Zhang, W., Zhang, F., Du, F. & Li, L. Chemical fertilizer use and efficiency change of main grain crops in China. Plant Nutr. Fert. Sci. 16, 1136–1143 (2010).

Google Scholar

Meng, Y., Xu, F., Yang, F., Xu, Y. & Gong, X. Situation of fertilizer application and analysis of fertilizer saving potential in crop farming in China. Phosphate Compd. Fert. 30, 1–4 (2015).

CAS Google Scholar

Yang, F. et al. Chemical fertilizer application and supply in crop farming in China in 2013. J. Plant Nutr. Fert. 21, 217–225 (2015).

Google Scholar

Yan, X., Jin, J. & Liang, M. Fertilizer use efficiencies and yield-increasing rates of grain crops in China. Soils. 49, 1067–1077 (2017).

Google Scholar

Xu, Y., Yang, F., Zhang, W., Meng, Y. & Jiang, Y. Status and problems of chemical fertilizer application in crop plantations of China from 2014 to 2016. J. Plant Nutr. Fert. 25, 11–21 (2019).

Google Scholar

Jin, S. & Zhou, F. Zero growth of chemical fertilizer and pesticide use: China's objectives, progress and challenges. J. Resour. Ecol. 9, 50–58 (2018).

Google Scholar

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2005 (China Statistics Press, 2005).

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2006 (China Statistics Press, 2006).

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2007 (China Statistics Press, 2007).

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2008 (China Statistics Press, 2008).

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2009 (China Statistics Press, 2009).

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2010 (China Statistics Press, 2010).

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2011 (China Statistics Press, 2011).

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2012 (China Statistics Press, 2012).

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2013 (China Statistics Press, 2013).

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2014 (China Statistics Press, 2014).

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2015 (China Statistics Press, 2015).

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2016 (China Statistics Press, 2016).

Department of Rural Social and Economic Investigation (Eds.): China Rural Statistical Yearbook-2017 (China Statistics Press, 2017).

Lin, B. & Li, J. Application of phosphate fertilizer and N:P ratio in China. Sciences in Agricultural Resources 3, 13–16 (2002).

ADS Google Scholar

Download references

This work is supported by the National Key Research and Development Project of China (2019YFA0607402), and the National Key Scientific and Technological Infrastructure Project Earth System Science Numerical Simulator Facility (EarthLab).

These authors contributed equally: Wenmeng Zhang, Tianyi Zhang.

College of Resources and Environmental Sciences, China Agricultural University, Beijing, China

Wenmeng Zhang & Xiaoguang Yang

State Key Laboratory of Atmospheric Boundary Layer Physics and Atmospheric Chemistry, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing, China

Tianyi Zhang

Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science & Technology, Nanjing, China

Tianyi Zhang

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Tianyi Zhang conceived the research idea; Tianyi Zhang and Xiaoguang Yang provide the raw data; Wenmeng Zhang calculated and verified the dataset; Tianyi Zhang wrote the paper with contributions from all authors.

Correspondence to Tianyi Zhang.

The authors declare no competing interests.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Zhang, W., Zhang, T. & Yang, X. 1 km-resolution gridded dataset of phosphorus rate for rice wheat and maize in China over 2004–2016. Sci Data 10, 363 (2023). https://doi.org/10.1038/s41597-023-02283-z

Download citation

Received: 07 February 2023

Accepted: 31 May 2023

Published: 07 June 2023

DOI: https://doi.org/10.1038/s41597-023-02283-z

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative