the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Century Long Reconstruction of Gridded Phosphorus Surplus Across Europe (1850–2019)
Abstract. Phosphorus (P) surplus in soils significantly contributes to the eutrophication and degradation of water quality in surface waters worldwide. Despite extensive European regulations, elevated P levels persist in many water bodies across the continent. Long-term annual data on soil P surplus are essential to understand these levels and guide future management strategies. This study reconstructs and analyzes the annual long-term P surplus for both agricultural and non-agricultural soils across Europe at a 5 arcmin (≈ 10 km at the equator) spatial resolution from 1850 to 2019. The dataset includes 48 P surplus estimates that account for uncertainties arising from different methodological choices and coefficients in major components of the P surplus. Our results indicate substantial changes in P surplus magnitude over the past 100 years, underscoring the importance of understanding a long-term P surplus. Specifically, the total P surplus across the EU-27 has tripled over 170 years, from 1.19 (±0.28) kg ha−1 yr−1 in 1850 to around 2.48 (±0.97) kg ha−1 yr−1 in recent years. We evaluated the plausibility and consistency of our P surplus estimates by comparing them with existing studies and identified potential areas for further improvement. Notably, our dataset supports aggregation at various spatial scales, aiding in the development of targeted strategies to address soil and water quality issues related to P. The P surplus reconstructed dataset is available at: https://doi.org/10.5281/zenodo.11351028.
- Preprint
(12080 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 15 Sep 2024)
-
RC1: 'Comment on essd-2024-294', Anonymous Referee #1, 26 Aug 2024
reply
This study reconstructs and analyzes long-term phosphorus (P) surplus in Europe from 1850 to 2019 at a high spatial resolution, providing a comprehensive dataset that includes 48 P surplus estimates and accounts for various uncertainties. The findings reveal that P surplus across the EU-27 has tripled over 170 years, highlighting the need for targeted management strategies. Additionally, the study identifies big variation in P surplus estimates by comparing its database with major previously published P databases and highlights areas for improvement. This research contributes to the improvement of quantification of historical and spatial P budgets in Europe by considering varying methods, including more budget terms, extending the temporal span of current studies, increasing resolution, and comparing various datasets. However, a few comments regarding the calculation details and discussions are provided for the authors' consideration.
Specific commentsSection 2.1:
- A more detailed definition of each rate term needs to be provided wherever a rate appears. Since the rates are always expressed in kg ha−1 yr−1 in this study, is the area being referred to the grid area or a specific land type area? Only part of the rate terms in this manuscript indicate which area is used as the denominator.
- Since more definitions of land types and other budget terms are provided in the following sections, I suggest that the authors move Section 2.1 to the end of Section 2.
Section 2.3
- Line 215: Do you mean the IFA report published in 2017 instead of 2013 that contains 2014-2015 data?
- Section 2.3.3: Not all terms in each equation are clearly defined in the text after the equation, in equations 15-23. For example, what is P2O5fercrops in Equation 15, what is nu and i in equations 22-23? It may seem redundant since some of the terms are discussed before the equations, but it would be easier for readers who are not familiar with these terms if each term's definition is clearly listed after the equations. Or, since this manuscript includes many parameters, it would be very helpful to include a table at the end listing all parameters' names, definitions, units, and/or calculation methods.
- Line 230: The fertilizer application rates here are defined as total fertilizer use/harvest area. I was wondering if the authors use the same definition for all fertilizer application rates discussed in the study, such as the data from Batool et al. (2022) mentioned later, and in equations 18-20. It seems that the fertilizer application rates are defined as total fertilizer use/land use area in equations 18-20. However, land use areas may not correspond to harvest areas.
- Equations 37-38: If the unit of Pman is kg yr−1, why are the units of Mantreat and Manleft (=Pman*ratio) kg ha−1 yr−1? Were any details omitted from the equations? Is the gridded area or the cropland+pasture area used as the denominator to calculate the application rate? It seems that you used the grid's physical area as the denominator based on Line 365. However, based on your Equations 41-42, it seems that the country total cropland+pasture area is used as the denominator (since you assume that "the entire amount of treated P manure is applied to soil" according to line 355 and “equal distribution rates for cropland and pasture within each grid cell” in Line 380)?
- Lines 430-435: You may want to add explanations for CW, C, and A after the equations or refer to a table listing all parameters. Do the area parameters C and A correspond only to the areas that are both of a specific land use type and have chemical weathering?
Sections 2.4.1-2.4.2:
- It seems that Section 2.4.2 would be better placed before Section 2.4.1.?
Section 3.1:
- Including a series of maps showing accumulated P surplus in the SI could be informative in illustrating how much P has accumulated or been lost in each area.
Figure 4:
- Are there any hypotheses that could explain the peaks in P surplus around 1980 and the subsequent decrease afterward? In addition, are there any reasons that can explain the varying uncertainty range over time based on your modeling methods? For example, could it be due to changes in crop portfolios over time or differences in data availability? Understanding these reasons could provide insights into improving the modeling in the future. The discussion can be included in a discussion section after the results.
Section 3.2:
- Some of the calculation details could be mentioned in the Methods section.
Lines 645-650:
- The discussion starting here can be included in a new Discussion section. Since your database only includes a limited number of scenarios and your parameter values are borrowed from previous studies, it would be informative to include a brief discussion that lists the key limitations of your database and suggests ways to improve them in the future, instead of simply providing some limitation examples. Other limitations could include, for example, the temporal and spatial variation of parameter values that are not/partially included in the current database and the uncertainty in applying country-level data to a high-resolution map.
Technical corrections
- Abstract & Introduction Line 15: It would be helpful to include a short note defining 'P surplus' (similar to the definitions you provided in Lines 30, 35 and 55) upon its first appearance in the abstract and introduction for readers who are unfamiliar with the term or have different definitions for it.
- Abstract line 5-10: Does the area in ha refer to the country's total land area or specifically to agricultural land?
- Line 35: Provide the full name of FAOSTAT the first time it appears.
- Line 55: You may want to move 'phosphorus (P)' to the first time the letter 'P' appears in the manuscript. In addition, there are multiple places in the manuscript where the same abbreviations are explained repeatedly (e.g., 'phosphorus (P)' appears multiple times).
- Figure 1: I can guess what the arrows in different colors represent, but it would improve clarity if the authors also mention their meanings in the caption.
Citation: https://doi.org/10.5194/essd-2024-294-RC1 -
AC1: 'Comment on essd-2024-294', Masooma Batool, 26 Aug 2024
reply
We have recently received comments on our manuscript via an email. We thank the reviewer for the positive comments and suggestions that would improve the clarity of the paper.
Comments: I recently read your preprint in Earth System Science Data titled “Century Long Reconstruction of Gridded Phosphorus Surplus Across Europe (1850 – 2019)”. Nice topic of research 😉.
I wanted to request a few specific corrections regarding the citation of a paper:
- Ensure that the paper is consistently cited as Muntwyler et al. (2024) in line 39 and Table 3, where currently 2023 is referenced.
- Please add "crop residues" to the description of my paper in Table 3 to accurately reflect the removal of P from the budget with crop residues (you might also want to merge indices a and e in the same table to maintain consistency).
- Distinguish more clearly between the studies by Panagos et al. (2022a) and Einarsson et al. (2020) and my study. Our methodologies and the purposes of the P budgets differ significantly; mine is derived from a process-based model, whereas theirs are empirical.
- Note that the scale of Panagos’ phosphorus budget is at the NUTS2 level, only the soil stock resolution is 500m.
Thank you for addressing these corrections. I believe they will enhance the accuracy and clarity of your research.
Citation: https://doi.org/10.5194/essd-2024-294-AC1 -
AC2: 'Reply on AC1', Masooma Batool, 26 Aug 2024
reply
We thank the Reviewer for the valuable suggestions. We provide a detailed response to the Reviewer’s comments in the Supplement
Data sets
Century Long Reconstruction of Gridded Phosphorus Surplus Across Europe (1850-2019) Masooma Batool et al. https://zenodo.org/records/11351028
Model code and software
Century Long Reconstruction of Gridded Phosphorus Surplus Across Europe (1850-2019) Masooma Batool et al. https://zenodo.org/records/11351028
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
186 | 64 | 18 | 268 | 6 | 6 |
- HTML: 186
- PDF: 64
- XML: 18
- Total: 268
- BibTeX: 6
- EndNote: 6
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1