# Mass balance and long-term soil accumulation of trace elements in arable crop systems amended with urban composts or cattle manure during 17 years

Aurélia Marcelline Michaud, Philippe Cambier, Valerie Sappin-Didier, Valentin Deltreil, Vincent Mercier, Jean-Noël Rampon, Sabine Houot

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Aurélia Marcelline Michaud, Philippe Cambier, Valerie Sappin-Didier, Valentin Deltreil, Vincent Mercier, et al.. Mass balance and long-term soil accumulation of trace elements in arable crop systems amended with urban composts or cattle manure during 17 years. Environmental Science and Pollution Research, 2020, 27, pp.5367-5386. 10.1007/s11356-019-07166-8. hal-02624633

HAL Id: hal-02624633

https://hal.inrae.fr/hal-02624633

Submitted on 26 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Mass balance and long-term soil accumulation of trace elements in arable crop systems amended with urban composts or cattle manure during 17 years

Aurélia Marcelline Michaud<fnref n="1" />,<fnref n="2" /> [ID] · Philippe Cambier<fnref n="1" /> · Valérie Sappin-Didier<fnref n="3" /> · Valentin Deltreil<fnref n="4" /> · Vincent Mercier<fnref n="1" /> · Jean-Noël Rampon<fnref n="1" /> · Sabine Houot<fnref n="1" />

Received: 26 July 2019 / Accepted: 24 November 2019

© The Author(s) 2019

## Abstract

Organic waste products (OWP) application to crop lands makes possible nutrients recycling. However, it can result in long-term accumulation of trace elements (TE) in soils. The present study aimed at (i) assessing the impact of regular applications of urban composts and manure on the TE contents of topsoils and crops in a long-term field experiment, (ii) comparing the TE mass balances with the stock variations of TE in soils, and (iii) proposing a prospective evaluation of this practice, based on estimated soil safe threshold values and simulations of soil TE accumulation for 100 years. In the long-term field experiment, physico-chemical properties and TE contents (Cd, Cr, Cu, Hg, Ni, Pb and Zn) have been measured in OWP, soils, plants and leaching waters for the period 1998–2015, and used for mass balance calculations and long-term simulations of TE accumulations. The composts of green wastes and sludge (GWS) and of municipal solid waste (MSW) were the OWP with the largest TE contents, while the farmyard manure tended to have the lowest. Repeated application of OWP led to significant accumulation of Zn and Cu in the topsoil layer (not for Cr, Cd, Hg, Ni, Pb), especially with GWS, without overpassing calculated protective threshold values. No effect of repeated application of OWP has been observed on TE contents in grains (wheat, maize, barley). The positive mass balance has been dominated by the input flux of TE through OWP and resulted in the observed increases of soil stocks for Cu and Zn. Prospective simulation of soil content evolution until 2100 showed that soil content reached 0.4 mg Cd kg⁻<fnref n="1" /> soil (GWS, MSW), 38 mg Cu kg⁻<fnref n="1" /> soil (GWS) and 109 mg Zn kg⁻<fnref n="1" /> soil (GWS), which remained lower than protective threshold values.

**Keywords** Field experiment · Heavy metals · Organic waste products · PNEC · Soil quality · Soil-plant transfer

Responsible Editor: Zhihong Xu

**Electronic supplementary material** The online version of this article (https://doi.org/10.1007/s11356-019-07166-8) contains supplementary material, which is available to authorized users.

<figure><img src="image_4.png" /></figure>

Aurélia Marcelline Michaud
aurelia.michaud@inra.fr

## Introduction

The biggest challenge that agriculture faces is to produce food enough in quantity and quality to feed increasing population, in the complex context of climate change. Agriculture sustainability needs to consider soil fertility, facing soil threats such as decreasing soil organic carbon (Richard et al. 2019). After industrial revolution, nutrient recycling has been neglected due to the worldwide increased use of mineral fertilizers (Richard et al. 2019), although manure application on soils has traditionally contributed to the recycling of nutrients they contain. Recycling of organic waste products (OWP) in agriculture contributes to circular economy through the recycling of nutrients and organic matter and represents a valuable alternative to their landfilling or incineration (Houot et al. 2014; Noirot-Cosson et al. 2016; Alvarenga et al. 2017). The large variety of waste sources and treatment processes result in large

<footnote n="1">¹ UMR ECOSYS, INRA, AgroParisTech, Université Paris Saclay, 78850 Thiverval-Grignon, France</footnote>

<footnote n="2">² UMR SAS, INRA, AgroCampusOuest, 65 rue de Saint Brieuc, 35042 Rennes, France</footnote>

<footnote n="3">³ UMR ISPA, INRA, 33883 Villenave d'Omon, France</footnote>

<footnote n="4">⁴ Veolia recherche & innovation, Chemin de la digue, 78603 Maisons-Laffitte, France</footnote>
diversity of OWP with various properties (Houot et al. 2014; Sharma et al. 2017). The use of OWP from urban and industrial origin may replace the traditional use of manure in areas where animal breeding is lacking and may be used as a substitute to chemical fertilizers (Noirot-Cosson et al. 2016; Sharma et al. 2017).

However, the recycling of OWP may represent potential risks to the environment and alimentation, due to their contents in contaminants, such as organic contaminants, biological pathogens and trace elements (Gigliotti et al. 1996; Weber et al. 2007; Alvarenga et al. 2017; Bourdat-Deschamps et al. 2017). The repeated use of OWP as organic fertilizer or soil amendments can thus result in long-term accumulation of contaminants in the soil (Delschen 1999) which can cause problems if certain concentration levels are exceeded (Moolenaar et al. 1997a). Prevention of trace elements (TE) accumulation is needed for sustainable agricultural productions (Moolenaar et al. 1997a). Among TE, the most concern is for non-essentials and toxic ones (e.g. As, Cd, Hg, Pb), while trace elements also considered as micronutrients (e.g. B, Cu, Zn) can be phytotoxic at concentrations that pose little risks to human and animal health (Lopez-Rayo et al. 2016). Nevertheless, Cu and Zn are micronutrients of most concern. Indeed, studies showed that TE accumulation including Cu and Zn was associated with decreases in diversity and activity of soil microorganisms, suggesting that the soil microorganisms are more sensitive to TE than higher plants, which can impair biogeochemical processes such as specific pathways of nutrient cycling and organic matter decomposition (McGrath et al. 1995; Charlton et al. 2016). Even with quality thresholds regulating the application of OWP onto soil, the regular application of OWP may increase TE concentrations in soils (Kirchmann et al. 2017) to levels toxic to soil microorganisms and affect biological processes (Moolenaar and Beltrami 1998). Recently, tools and framework have been developed to estimate thresholds values and provide a risk assessment based on eco-toxicological studies (Heemsbergen et al. 2009; Smolders et al. 2009): the predicted-no-effect concentrations (PNECs) calculator and the “threshold calculator for metals in soil” (Oorts et al. 2018). These methods not only consider total TE concentration in soils which is a poor indicator of TE toxicity risks, but also soil physico-chemical properties known to influence the speciation of TE and alter their toxicity and the specific sensitivity of various species to TE toxicity. Risks associated to relatively low TE concentrations found in soils after repeated applications of OWP have seldom been studied, except for example by Heemsgergen et al. (2009) about the guidelines for amended soils with biosolids.

The balance of nutrients, including elemental major inputs and outputs within the field boundary, is an important approach to evaluate the sustainability of an agrosystem and the effects of long-term activities on the soil (Yu et al. 2011). This balance approach has proven to be useful in soil fertility

studies to demonstrate depletion of nutrients from soils. It has also been used in cases of soil pollution to demonstrate accumulation of TE in agricultural soils (Moolenaar and Lexmond 1998). Further use of the mass balance is to calculate the long-term change of TE content in the upper layer of soil. As presented by Moolenaar and Lexmond (1998), the change in total TE content in the topsoil layer results from the balance between inputs and outputs. Recently this mass balance approach has been used at the regional and national scale to simulate the evolution of TE concentration in the topsoil (Six and Smolders 2014; Sterckeman et al. 2018). Analysis at national level cannot consider relevant processes on smaller scales because of the average values and estimations used at larger scale. Belon et al. (2012) recognized that their approach at the department scale was not sufficient to help farmers in making decisions. They also pointed out that at the plot scale the relative importance of the diverse sources of TE may differ from the national and department ones. Furthermore, a balance assessment at the larger scale may lead to wrong conclusions locally. To that respect, it is important to complete and validate national-scale studies by field-scale studies considering precise agronomical practices and actual data for the agricultural compartments (i.e. fertilizer and amendment inputs, atmospheric depositions, plant offtake, leaching) (Moolenaar et al. 1997b; Moolenaar and Lexmond 1998). Nevertheless, the assessment of TE balance is often impossible in many situations due to the lack of outputs data with crop offtakes and leaching (Moolenaar et al. 1997a; Yu et al. 2011). In addition, models and simulations are not a substitute for real data and data obtained from long-term studies are essential for setting legislative limits for concentrations of TE in agricultural soils (Uprety et al. 2009).

Most field experimentations devoted to the evaluation of OWP application effects have been carried out in short-term studies (Douglas et al. 2003; Warman and Termeer 2005), often focusing on nitrogen fertilization (Bell and Leclerc 2015). As underlined by Lopez-Rayo et al. (2016), long-term experiments can give realistic overview of cumulative effects on soil properties and plant quality when applying OWP. Among short- and long-term field experiments devoted to the study of OWP application effects on TE dynamic, sewage sludge and manures were mainly studied (McGrath et al. 1995; Chaudri et al. 2000, 2007; Uprety et al. 2009; Hamner and Kirchman 2015; Charlton et al. 2016). Few short-term experiments were carried out with municipal solid waste or sewage sludge composts (Giusquiani et al. 1994; Businelli et al. 1996; Gigliotti et al. 1996; Weber et al. 2007). Much less studies considered various urban OWP including composts in long-term experiments: the work of Delschen (1999), Riber et al. (2014) and associated studies carried out in the CRUCIAL field experiment, and Houot et al. (2002) and related studies on the “QualiAgro” field experiment. To our knowledge, none of these studies provided actual TE mass
balance at the field scale with in situ measured data (i.e. inputs, outputs, and total balance) after repeated application of various urban OWP, and tested the mass balance model to simulate long-term accumulation of TE in the topsoil layer.

Thus, the aim of the present study was to (i) assess the impact of regular applications of urban composts and manure on the TE contents of topsoils and crops in the “QualiAgro” field experiment, (ii) compare mass balances with the stock variations of TE in the topsoil layer, and (iii) propose a prospective evaluation based on soil safe threshold values and simulation of TE accumulation in the topsoil layer for 100 years of agricultural repeated spreading of various OWP.

# Materiel and methods

## QualiAgro field experiment

The QualiAgro field experiment was set up in 1998 to investigate the long-term effects of repeated applications of urban composts on soil fertility and potential contamination with inorganic and organic contaminants (Houot et al. 2002; Cambier et al. 2019). It is located within the Plateau des Alluets (Yvelines, France), on a silt loam Luvisol, according to IUSS-FAO classification and USDA texture triangle, developed on aeolian loess. The soil of about more than 1.2 m deep is developed on carbonated loess deposits. The upper horizons include a plough layer (i.e. topsoil) and a plough pan layer or ancient ploughed layer between 28 and 35 cm. Horizons of clay accumulation appear around 50–60 cm depth. All these horizons are decarbonated, the carbonated silt loess appearing below 1.4–1.6 m depth. The initial physico-chemical characteristics of the topsoil layer are presented in Table 1. From 1999 to 2015, the 6-ha field experiment has been cropped with a biennial rotation of maize and winter wheat, except in 2007, when spring barley was inserted because of infestation by chrysomela in the region (Appendix A in Electronic

supplementary material). After grain harvest, the maize residues are incorporated into the topsoil layer but the wheat crop residues are exported. The climate is oceanic, with mean annual precipitation of 572 mm and mean annual temperature of 11 °C. The experiment could be considered as representative of cultivated soils of the central Paris Basin and of extended agricultural areas of NW Europe.

Figure 1 presents the experimental plan of the QualiAgro field experiment, with the distribution of the 20 plots of 45 × 10 m² in 4 blocks of replicates of the following 5 treatments randomly distributed within each block:

(1) Co-compost of sewage sludge with green waste and wood chips (GWS),

(2) Co-compost of the home-sorted fermentable fraction of municipal solid waste and green waste, also called biowaste compost (BIOW),

(3) Compost of residual municipal solid waste after sorting of recyclable packaging and papers (MSW),

(4) Farmyard manure from a dairy farm, used as a reference organic amendment (FYM),

(5) No organic amendment (control or CN).

For the period 1998–2013, additional mineral nitrogen fertilization was applied in the form of “Solution 39” (urea + ammonitrate), at doses between 100 and 170 kg N ha⁻¹ on wheat and between 0 and 80 kg N ha⁻¹ on maize, depending on mineral N stocks measured in plots down to 90 cm, at the end of winter before fertilization (Noirot-Cosson et al. 2016). No mineral fertilization of potassium and phosphorus was applied during the period 1998–2013. All treatments and plots received similar pesticide treatments and cropping practices. In 2014 and 2015, according to the new management strategy of the experiment which has been converted into organic farming practices, the additional fertilization was applied with the organic commercial product “Axe N+P”.

Each OWP application was made on the wheat stubbles every 2 years in September, at doses equivalent to 4 t organic C ha⁻¹. This corresponded to average doses of 18.3 (MSW),

<table><caption>**Table 1** Chemical properties and trace element contents measured in the topsoil layer (0–25 cm, sieved < 2 mm) at the beginning of the experiment in 1998, with the following: texture (g kg⁻¹ DM; DM for dry matter), pH, cationic exchange capacity (CEC; cmol⁺ kg⁻¹ DM),</caption><thead><tr><th>Clay</th><th>Silt</th><th>Sand</th><th>pH</th><th>CEC</th><th>Org. C</th><th>Total N</th><th>Olsen P<sub>2</sub>O<sub>5</sub></th></tr></thead><tbody><tr><td rowspan="2">145.8 (7.5)</td><td>g kg<sup>-1</sup> DM</td><td>69.2 (4.5)</td><td>7.06 (0.19)</td><td>cmol<sup>+</sup> kg<sup>-1</sup> DM</td><td>g kg<sup>-1</sup> DM</td><td>1.12 (0.05)</td><td>0.09 (0.01)</td></tr><tr><td>785.0 (5.4)</td><td></td><td></td><td>9.7 (0.8)</td><td>10.6 (0.6)</td><td></td><td></td></tr><tr><td></td><td>Cd</td><td>Cr</td><td>Cu</td><td>Hg</td><td>Ni</td><td>Pb</td><td>Zn</td></tr><tr><td rowspan="2">Referential<sup>a</sup></td><td>0.24 (0.01)</td><td>45.6 (3.1)</td><td>12.0 (1.0)</td><td>mg kg<sup>-1</sup> DM</td><td></td><td></td><td></td></tr><tr><td>0.25 (0.16)</td><td>41 (12)</td><td>12 (11)</td><td>0.10 (0.02)</td><td>14.8 (0.7)</td><td>25.7 (4.1)</td><td>51.6 (3.5)</td></tr><tr><td></td><td></td><td></td><td></td><td>0.07 (0.16)</td><td>18 (6)</td><td>23 (12)</td><td>52 (20)</td></tr></tbody></table>

content of organic carbon (Org. C), total nitrogen (Total N) and Olsen extractable phosphorus (Olsen P₂O₅) (g kg⁻¹ DM), trace element content for Cd, Cr, Cu, Hg, Ni, Pb and Zn (mg kg⁻¹ DM); mean values (standard deviations for all plots), compared to the mean region contents

<footnote>ª Topsoil layers of agricultural soils of “Ile-de-France” (central Paris Basin): medians (standard deviations), *n* between 2000 and 4000, according to Baize et al. (2007)</footnote>
<figure><img src="image_5.png" /><figcaption>**Fig. 1** Experimental plan of the QualiAgro device with the following: GWS for co-compost of sewage sludge with green waste and wood chips, BIOW for co-compost of biowaste and green waste, MSW for compost of</figcaption></figure>

organic fraction of municipal solid waste, FYM for farmyard manure of
dairy cows on straw, CN for control without amendments

26.2 (GWS), 27.7 (BIOW) and 34.6 (FYM) tons of fresh matter per hectare, i.e. about 1.5- to 3-fold larger than usual practices (average recommended amounts: 20–30 t ha⁻¹ per 3 years for farmyard manure and 10–20 t ha⁻¹ per 3 years for composts). More precise characteristics of each amendment are given in “Characteristics of applied organic waste products” section.

Sampling and analyses

Soils were sampled before each OWP spreading. Each plot was sampled in the topsoil layer between 0 and 25 cm depth (10 spadesful further pooled). After air drying, the fine earth (< 2 mm) was prepared according to the standard ISO 11464 and kept for further analyses. In 1998, 2010 and 2015, soils samples were also taken from the layer located below the topsoil layer (25–35 cm) for physico-chemical analyses.

In 2013–2014 and 2014–2015, soil solutions were sampled for further analyses at the depth of 45 cm during drainage seasons with lysimeters with fibreglass wicks (45 cm deep) (Cambier et al. 2014; Chabauty et al. 2016).

Between 1998 and 2013, a total of 9 OWP spreading have been carried out (Appendix A in Electronic supplementary material). All OWP were sampled during application for analysis.

Just before the mechanical harvest, plants were sampled (i.e. for each plot harvested on 10 subplots of 0.5 m triple row for wheat and 5 subplots of 2.5 m double row for maize, pooled to obtain one composite sample per plot) and grains and aerial crop residues separated to measure crop yields and biomasses of residues.

The soil bulk density was measured in the top soil layers simultaneously to soil sampling before OWP spreading, in each plot in 1998, 2004, 2009 and 2013, by using vertical stainless cylinders (French standard X31-501 or ISO 11272).

Soils, OWP and leachates were analysed by the central soil testing laboratory of INRA (INRA LAS, Arras) according to standardized French (AFNOR) or international procedures (ISO). Soil and OWP pH were measured in water using a 1:5 soil/solution ratio (NF EN 15933). Soil organic carbon was determined by dry combustion (corrected from the carbonates) (NF ISO 10694); effective cation exchange capacity (CEC) was determined by cobalthexamine chloride extraction (NF ISO 23470). In OWP and soil samples, total trace element concentrations of Cd, Cr, Cu, Ni, Pb and Zn were determined on 250 µm crushed samples by HF–HClO₄ extraction (NF X31-147). Trace elements analyses were performed by inductively coupled plasma atomic emission spectroscopy (ICP-AES) for Cr, Cu, Ni and Zn or inductively coupled plasma mass spectrometry (ICP-MS) for Cd and Pb. Total Hg concentration was estimated by dry combustion on 250 µm crushed samples. Other chemical properties were determined as follows: in soils, N total content by NF ISO 13878, P extracted by the Olsen method NF ISO 11263; in OWP, organic N was calculated from total N (dry combustion) and mineral N, organic C by dry combustion, P by HF–HClO₄ extraction NF X31-147, CaCO₃ content by NF ISO 10693. TE concentrations in leachates were analysed according to standard methods by ICP-MS.

All plant analyses were carried out in the central plant testing laboratory of INRA (INRA USRAVE, France) according to standardized French (AFNOR) or international procedures (ISO). Grains and crop residues were oven dried (T° < 40 °C) then ground under 0.5 mm. Plant samples were digested in a microwave oven with a sequential procedure including digestion with concentrated nitric acid, H₂O₂ and HF–HClO₄. Trace elements analyses were performed by ICP-AES for Cu, Zn or ICP-MS for Cd, Cr, Ni and Pb. Total Hg concentration was analysed by dry combustion on < 0.5 mm crushed samples.
# Calculations

## Soil stocks of trace elements in the topsoil layer

Organic waste products have been incorporated into the topsoil layer. The stock $Q_{top,i,n}$ (kg ha⁻¹) of any element contained in the topsoil layer is given by:

$$Q_{top,i,n} = [TE]_{top,i,n} \times W_{top,i,n} \times 0.001 \\ = [TE]_{top,i,n} \times h \times d_{top,i,n} \times 10 \quad (1)$$

Where,

top refers to the topsoil layer, i to the treatment and n to the year.

$W_{top,i,n}$ is the soil mass (t ha⁻¹) in the topsoil layer (Appendix B in Electronic supplementary material).

$[TE]_{top,i,n}$ is the trace element content (mg kg⁻¹ or g t⁻¹) (Appendix C in Electronic supplementary material).

h is the height of the soil layer.

$d_{top,i,n}$ is the soil density.

Since the beginning of the experiment, the soil bulk density was measured in 1998, 2004, 2009 and 2013, and calculated using an adjusted polynomial equation for other sampling dates (Appendix B in Electronic supplementary material). It increased in the topsoil layers, more in the control treatment than in the organic treatments. This led to an increase of the soil mass since 1998, especially in the control treatment. Nevertheless, soil stocks must be calculated on a fixed soil mass basis (Ellert and Bettany 1995; Peltre et al. 2012). As the soil bulk density of the control treatment in 2015 was the largest one of all, its mass was considered as the basis for further calculations of corrected soil stock. The mass $W_{top,CN,2015}$ of 3800 t ha⁻¹ was thus used, corresponding to the bulk soil density of 1.52 and a soil depth of 0.25 cm.

According to Peltre et al. (2012), we calculated corrected stocks in the topsoil layer of year n ($Q_{top,i,n}$) by adding a stock of an element present in the sub-layer progressively incorporated between 1998 and year n ($Q_{sub,i,n}$).

The stock of TE from the sub-layer to be added to the topsoil layer stock was calculated with:

$$Q_{sub,i,n} = (W_{top,CN,2015} - W_{top,i,n}) \times [TE]_{sub,i,n} \quad (2)$$

Where,

sub refers to the sub-layer (25–35 cm), i to the treatment and n to the year.

$W_{top}$ is the soil mass in the topsoil layer.

$[TE]_{sub,i,n}$ is the trace element content (mg kg⁻¹) in the sub-layer (Appendix D in Electronic supplementary material).

The TE contents of sub-layers under the topsoil layer have been measured in 1998 in a few places, and in 2010 and 2015 for each plot (Table D.1. in Appendix D, Electronic supplementary material). Considering the lack of data for the

beginning of the experiment in 1998, $[TE]_{sub,CN,2015}$ (with sub for sub-layer, CN for the control treatment and 2015 for the considered year) was considered as the reference for $[TE]_{sub,i,1998}$ at the beginning of the experiment (i.e. 1998) in all treatments, while the $[TE]_{sub,i,2015}$ measured in the OWP treatments in 2015 was considered as the final contents in 2015. Then, we used linear functions to interpolate $[TE]_{sub,i,n}$ between 1998 and 2015, $[TE]_{sub,CN,n}$ being assumed constant (detailed calculations and equations in Appendix D):

$$[TE]_{sub,i,n} = a \times (n-1998) + [TE]_{sub,CN,2015} \quad (3)$$

The corrected stock of TE in the topsoil layer ($Q_{top}$), corresponding to the constant mass of soil of 3800 t ha⁻¹, was calculated with (detailed calculations in Appendix E):

$$Q_{top,i,n} = Q_{top,i,n} + Q_{sub,i,n} \quad (4)$$

Where, $Q_{top,i,n}$ (kg ha⁻¹) is the stock of TE contained in the topsoil layer and $Q_{sub,i,n}$ is the stock of TE from the sub-layer to be added to the stock of the topsoil layer.

## Mass balance calculations

Input and output fluxes of TE were calculated at the field scale (kg ha⁻¹): inputs through OWP applications and atmospheric deposits; outputs through leaching water and crop uptakes. The input fluxes due to organic fertilizers in 2014 and 2015 represented less than 0.7% of the total input of TE with OWP application so that it has been neglected in the following calculations (data in Appendix G).

A balance of trace elements in soil relates the rates of accumulation input – output (Moolenaar et al. 1997b),

$$\text{Mass balance} = \sum Q_{\text{input}} - \sum Q_{\text{output}} \quad (5)$$

Where,

$$Q_{\text{input}} = \sum Q_{\text{atm},i,n} + \sum Q_{\text{OWP},i,n} \quad (6)$$

$$Q_{\text{output}} = \sum Q_{\text{lea},i,n} + \sum Q_{\text{crop},i,n} \quad (7)$$

$Q_{atm,i,n}$ (kg ha⁻¹) is the input flux related to atmospheric deposits of TE, calculated from Azimi et al. (2004) with mean TE atmospheric deposition rate obtained in the Parisian region (g ha⁻¹ year⁻¹).

$Q_{OWP,i,n}$ (kg ha⁻¹) is the TE input flux related to OWP application, calculated for each application from measured applied amount of OWP ($q_{OWP,i,n}$ in tons DM ha⁻¹) and TE measured in the corresponding OWP ($[TE]_{OWP,i,n}$ in mg kg⁻¹ DM) ($Q_{OWP,i,n} = [TE]_{OWP,i,n} \times q_{OWP,i,n}$) (detailed OWP characteristics used for $Q_{OWP}$ calculations in Appendix F).

$Q_{lea,i,n}$ (kg ha⁻¹) is the output flux of TE associated to leached waters collected at 45 cm depth, based on the mean
annual flux measured during 2 drainage periods (2013–2015) (detailed leaching waters characteristics used for $Q_{lea,i,n}$ calculations in Appendix H).

$Q_{crop,i,n}$ (kg ha⁻¹) is the output flux of TE related to crop uptake, calculated each year from the yields measured for each crop ($q_{crop,i,n}$, ton DM ha⁻¹) for exported grains and wheat crop residues and TE concentration measured in corresponding plant samples (mg kg⁻¹ DM) ($Q_{crop,i,n} = [TE]_{crop,i,n} \times q_{crop,i,n}$ × 0.001) (detailed crop characteristics used for $Q_{crop,i,n}$ calculations in Appendix I). For the $[TE]_{crop,i,n}$ lower than the quantification limits, we used raw values when they were given by the laboratory or quantification limit divided by 2 were considered to estimate $Q_{crop,i,n}$.

## Simulation of soil trace element content

The increase in TE content in the topsoil layer is the result of the net balance between input and output flows divided by the plough layer weight (Moolenaar et al. 1997b).

We adapted the equations of Six and Smolders (2014) and Sterckeman et al. (2018) to calculate annually the TE content in the topsoil layer and simulate long-term evolution of trace element contents in the topsoil layer, i.e. $[TE]_{top,i,n}$:

$$[TE]_{top,i,n} = [TE]_{top,i,n-1} + \frac{Qm \text{ input} - Qm \text{ output}}{W_{top}} \quad (8)$$

Where,

The soil TE content over the year $n$ ($[TE]_{top,i,n}$, in mg TE kg⁻¹ soil) is the result of that in year $n-1$ ($[TE]_{top,i,n-1}$, in mg TE kg⁻¹ soil) and of the balance between the mean input ($Q_{m,i}$ input, g ha⁻¹) and the mean output ($Q_{m,i}$ output, g ha⁻¹) of TE in the considered topsoil layer, whose weight is $W_{top}$ (3 800 tonnes ha⁻¹).

$Q_{m,i}$ input and $Q_{m,i}$ output correspond to average annual quantity of TE calculated from data presented in “Mass balance calculations” section (detailed mean annual fluxes data in Appendix J):

$$Q_{m,i} \text{ input} = \frac{\sum Q_{atm} + \sum Q_{OWP}}{\text{Number of years}} \quad (9)$$

$$Q_{m,i} \text{ output} = \frac{\sum Q_{lea} + \sum Q_{crop}}{\text{Number of years}} \quad (10)$$

## Threshold contents for trace elements in soil

Ecological threshold contents in soils for trophic levels (i.e. plants, invertebrates and microbial processes) were estimated with the threshold calculator for metal in soil V2.0® (Arche Consulting, https://www.arche-consulting.be/) (Oorts et al. 2018). The calculator has been developed using the studies of Smolders et al. (2009) and Checkai et al. (2014) based on

datasets derived from scientific literature, research projects and European REACH dossiers (Registration, Evaluation, Authorisation and Restriction of Chemicals EC No 1907/2006).

The 5th percentile of a chronic toxicity distribution has been chosen under the EU REACH Regulation as a concentration that is protective for most species in a community (namely 95%). Therefore, the PNEC value (i.e. predicted no effect concentration) is based on the 5% probability level (i.e. 5% Hazardous Concentration or $HC_5$) of the reliable EC10 (i.e. Effective Concentration for which 10% of the community is affected by the concentration of the trace element) divided by an assessment factor between 1 and 5 (ECHA 2008; OECD 2016). Finally, the $HC_5$ of EC 10 values, without additional assessment factor, has been chosen as a protective threshold value for Cd, Cu, Ni, Pb and Zn for the trophic levels plants, invertebrates and microorganisms. No sufficient soil eco-toxicity data are available for Hg and Cr to use a species sensitivity distribution approach, nor for implementation of bioavailability corrections for the derivation of ecological threshold concentrations of these elements.

The soil properties required as input parameters for bio-availability corrections were measured by the adequate methods mentioned by Oorts (2018) for organic C, clay content and effective cation exchange capacity. Only pH 0.01 M $CaCl_2$ was estimated from the measured pH$_{H2O}$ by using the proposed equation $pH_{CaCl_2} = -0.54 + 1.00 \times pH_{H2O}$.

## Statistic treatments

Parametric or non-parametric tests were applied after testing normality of residues and equality of variances. Tests were performed as follows using StaBoxAgri® (2009) and XLSTAT® (V 2018.2, Addinsoft): ANOVA parametric and Kruskal-Wallis non-parametric tests were performed to identify the treatment effect per year; Newman-Keuls parametric or Bilateral Dunn non-parametric tests were performed to identify significant differences between treatments per year; Friedman's non-parametric test was performed to evaluate significant effect of years with Bilateral Nemenyi non-parametric method used for multiple comparisons for paired samples.

## Results and discussion

### Characteristics of applied organic waste products

The 4 OWP applied in the QualiAgro field experiment differed in their physico-chemical characteristics (Table 2) with average characteristics in agreement with previous studies (Sen Tran et al., 1996; Houot et al. 2014). As already stressed by Obriot et al. (2016), the largest pH was found in FYM and
BIOW, largest content in carbonates in BIOW, largest organic carbon content in FYM, largest organic nitrogen content in GWS and FYM, and largest content in phosphorus and potassium were found in GWS and FYM, respectively.

Average contents in total TE in the OWP are presented in Table 2 and Fig. 2. Over the period 1998–2015, all GWS composts presented TE contents lower than regulatory thresholds defined in NFU 44-095 in application since 2002 for such composts. In 2000, the BIOW composts were above the thresholds for Pb and Ni defined in NFU 44-051 in application since 2006. Surprisingly, the contents in Cd, Pb, Zn were above the thresholds for FYM in 1998. The MSW composts overpassed the limits more often: for Cr and Pb in 1998, 2002, 2004; for Cu, Hg and Ni in 1998; for Zn in 2000 and 2013; but MSW composts conformed to the regulation in 2006, 2007, 2009, 2011, and Cd was always under the regulatory limit.

For most TE, the contents in GWS and MSW composts often appeared above those in BIOW and FYM, although not always significantly (Table 2; statistics in Appendix F in Electronic supplementary material). On average, total contents in TE were larger in the MSW compost for Cd, Cr, Ni and Pb compared to the other OWP, in the GWS compost for Cu and Zn. The TE contents were in agreement with previous studies (McBride and Spiers 2001; Zhang et al. 2012; Alvarenga et al. 2015, 2017).

An overall decrease of contents of most TE was observed with time, in the different OWP except for GWS compost, suggesting improvement of the OWP quality with time (Fig. 2). Therefore, a more meaningful comparison between OWP is given below through the cumulated TE inputs (Table 3).

The TE input for each spreading and their accumulation were calculated and reported to a 10-year period for comparison with the relevant regulatory French limits, i.e. in application since 2002 for NFU 44-095 and since 2006 for NFU 44-051. Despite a total input of Zn and Cu which exceeded the French regulation for GWS, and of Pb with the MSW and BIOW composts, the OWP studied in QualiAgro field

experiment were globally in accordance with French national regulations. The input of Pb was about two times less with GSW applications than with other products; the input of Ni was low with FYM; the input of Hg was lower with BIOW and FYM than with GWS and MSW applications. The general observation, with a few exceptions, is that the inputs for each TE were comparable among treatments and within the range of the regulatory maximum. It must be noticed that amounts of applied OWP have been about 1.5- to 3-fold larger than usual practices, which led to maximize risks related to input fluxes of TE. Copper and zinc were added in greater amounts (6.6–15.0 kg ha⁻¹ 10 years⁻¹ for Cu; 23.9–35.1 kg ha⁻¹ 10 years⁻¹ for Zn) compared to the other TE. Mercury and cadmium were added in the rates respectively of 0.01–0.07 kg ha⁻¹ 10 years⁻¹ and 0.07–0.09 kg ha⁻¹ 10 years⁻¹. Repeated spreading of GWS led to greater amounts of Cd, Cu, Hg and Zn. Spreading of BIOW led to greater amounts of Ni and Pb, and MSW presented the largest input fluxes of Cd and Cr. Spreading of FYM led globally to the lowest inputs of TE.

## Effect of organic waste product applications on soil physico-chemical properties and trace elements contents

The physico-chemical characteristics of the topsoil layer sampled in 1998 before the first OWP application have been summarized in Table 1. Before the first OWP application, the total contents of TE did not significantly differ between the plots selected for each treatment (P = 0.05, n = 4, test of Newman-Keuls for Cd, Cr, Cu, Hg, Ni and Zn, non-parametric test for Pb). Initial TE contents in the topsoil layer were in the range of those found in similar soils of the Parisian region (Baize et al. 2007), and with median contents found in France (BDETM 2019) (Table 4) with the following sequence: Zn > Cr > Pb > Ni > Cu > Cd > Hg. The initial soil pH was 7.06, and the organic carbon content was in average 10.61 g kg⁻¹ of dry soil.

<caption>**Table 2** Physico-chemical properties of the organic waste products studied in the QualiAgro field experiment, as compared to the French regulation level NFU 44 051/095 for each trace element, with the following: dry matter (DM) expressed as percentage of fresh matter (FM), pH, content of CaCO₃, organic carbon (Org. C) and nitrogen (Org. N), total content of phosphorus (P₂O₅) and potassium (K) (g kg⁻¹ DM), trace element content for Cd, Cr, Cu, Hg, Ni, Pb and Zn (mg kg⁻¹ DM); average values for the period 1998–2013 ± standard deviations (n = 9); co-compost of sewage sludge and green waste (GWS), biowaste compost (BIOW), compost of residual municipal solid waste (MSW), farmyard manure (FYM)</caption>

<table><thead><tr><th>OWP</th><th>DM</th><th>pH</th><th>CaCO<sub>3</sub></th><th>Org. C</th><th>Org. N</th><th>P<sub>2</sub>O<sub>5</sub></th><th>K</th><th>Cd</th><th>Cr</th><th>Cu</th><th>Hg</th><th>Ni</th><th>Pb</th><th>Zn</th></tr><tr><td colspan="2">% FM</td><td colspan="6">g kg<sup>-1</sup> DM</td><td colspan="6">mg kg<sup>-1</sup> DM</td></tr></thead><tbody><tr><td>GWS</td><td>63 ± 9</td><td>7.5 ± 0.6</td><td>28 ± 12</td><td>265 ± 45</td><td>21 ± 2</td><td>29 ± 8</td><td>15 ± 5</td><td>1.1 ± 0.7</td><td>40 ± 7</td><td>174 ± 43</td><td>0.8 ± 0.3</td><td>27 ± 7</td><td>62 ± 11</td><td>406 ± 89</td></tr><tr><td>BIOW</td><td>70 ± 9</td><td>8.1 ± 0.5</td><td>90 ± 62</td><td>211 ± 47</td><td>17 ± 4</td><td>11 ± 4</td><td>21 ± 3</td><td>0.7 ± 0.7</td><td>38 ± 15</td><td>62 ± 26</td><td>0.2 ± 0.1</td><td>25 ± 24</td><td>87 ± 56</td><td>235 ± 100</td></tr><tr><td>MSW</td><td>68 ± 13</td><td>7.5 ± 0.5</td><td>69 ± 23</td><td>310 ± 45</td><td>17 ± 2</td><td>8 ± 2</td><td>10 ± 3</td><td>1.4 ± 0.7</td><td>86 ± 57</td><td>143 ± 87</td><td>0.8 ± 0.8</td><td>30 ± 18</td><td>140 ± 87</td><td>388 ± 206</td></tr><tr><td>FYM</td><td>40 ± 10</td><td>9.1 ± 0.3</td><td>47 ± 17</td><td>324 ± 68</td><td>21 ± 3</td><td>13 ± 3</td><td>35 ± 3</td><td>1.1 ± 1.1</td><td>34 ± 24</td><td>95 ± 70</td><td>0.1 ± 0.0</td><td>12 ± 9</td><td>115 ± 127</td><td>342 ± 189</td></tr><tr><td>NFU 44 051/095</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>3</td><td>120</td><td>300</td><td>2</td><td>60</td><td>180</td><td>600</td></tr></tbody></table>
<figure><img src="image_6.png" /><figcaption>**Fig. 2** Evolution of trace element contents in applied organic waste products (OWP) over the period 1998–2015 for Cd, Cr, Cu, Hg, Ni, Pb and Zn (mg kg⁻¹ DM; DM for dry matter). Diamonds stand for co-compost of sewage sludge and green waste (GWS), squares stand for biowaste compost (BIOW), triangles stand for compost of residual municipal solid waste (MSW) and circles stand for farmyard manure (FYM). The hatched rectangle corresponds to the French regulation level for each trace element. The error bars are standard deviations</figcaption></figure>

Several significant and lasting differences of soil physico-chemical properties appeared among the treatments at dates comprised between 2004 and 2015, depending on the element

considered (Peltre et al. 2012; Obriot et al. 2016) (Appendix C in Electronic supplementary material). A decrease in soil organic C content (SOC) occurred in the control treatment CN and it
<caption>**Table 3** Input fluxes by organic waste products (OWP) studied in the QualiAgro field experiment, with the following: the applied quantity expressed as tons of dry matter per hectare (t DM ha⁻¹), the input flux for organic carbon (Org. C) and the trace elements (Cd, Cr, Cu, Hg, Ni, Pb, Zn) expressed as kg ha⁻¹; average input per spreading, cumulated</caption>

input over the period 1998–2015, cumulated input relative to 10 years of spreading as compared to the French regulation level NFU 44 051/095 (NFU) for each trace element; co-compost of sewage sludge and green waste (GWS), biowaste compost (BIOW), compost of residual municipal solid waste (MSW), farmyard manure (FYM)

<table><thead><tr><th>OWP</th><th>Input type</th><th>Applied quantity</th><th>Org. C</th><th>Cd</th><th>Cr</th><th>Cu</th><th>Hg</th><th>Ni</th><th>Pb</th><th>Zn</th></tr><tr><td></td><td></td><td>t DM ha<sup>-1</sup></td><td>kg ha<sup>-1</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></thead><tbody><tr><td>GWS</td><td>Average input</td><td>16.5</td><td>4357</td><td>0.017</td><td>0.7</td><td>2.8</td><td>0.013</td><td>0.4</td><td>1.0</td><td>6.6</td></tr><tr><td>BIOW</td><td>Average input</td><td>19.1</td><td>3922</td><td>0.014</td><td>0.8</td><td>1.2</td><td>0.004</td><td>0.5</td><td>1.8</td><td>4.7</td></tr><tr><td>MSW</td><td>Average input</td><td>12.2</td><td>3745</td><td>0.017</td><td>1.1</td><td>1.8</td><td>0.010</td><td>0.4</td><td>1.7</td><td>5.0</td></tr><tr><td>FYM</td><td>Average input</td><td>12.8</td><td>4058</td><td>0.015</td><td>0.5</td><td>1.3</td><td>0.001</td><td>0.2</td><td>1.5</td><td>4.5</td></tr><tr><td>GWS</td><td>Cumulated input</td><td>149</td><td>39213</td><td>0.150</td><td>5.9</td><td>25.6</td><td>0.119</td><td>4.0</td><td>9.1</td><td>59.6</td></tr><tr><td>BIOW</td><td>Cumulated input</td><td>172</td><td>35302</td><td>0.126</td><td>6.8</td><td>10.9</td><td>0.032</td><td>4.8</td><td>15.9</td><td>42.3</td></tr><tr><td>MSW</td><td>Cumulated input</td><td>110</td><td>33705</td><td>0.154</td><td>9.5</td><td>16.0</td><td>0.089</td><td>3.4</td><td>15.5</td><td>45.3</td></tr><tr><td>FYM</td><td>Cumulated input</td><td>115</td><td>36524</td><td>0.133</td><td>4.1</td><td>11.3</td><td>0.012</td><td>1.5</td><td>13.9</td><td>40.6</td></tr><tr><td>GWS</td><td>Relative to 10 years</td><td>87</td><td>23066</td><td>0.088</td><td>3.5</td><td>15.0</td><td>0.070</td><td>2.4</td><td>5.4</td><td>35.1</td></tr><tr><td>BIOW</td><td>Relative to 10 years</td><td>101</td><td>20766</td><td>0.074</td><td>4.0</td><td>6.4</td><td>0.019</td><td>2.8</td><td>9.4</td><td>24.9</td></tr><tr><td>MSW</td><td>Relative to 10 years</td><td>65</td><td>19826</td><td>0.091</td><td>5.6</td><td>9.4</td><td>0.052</td><td>2.0</td><td>9.1</td><td>26.7</td></tr><tr><td>FYM</td><td>Relative to 10 years</td><td>68</td><td>21485</td><td>0.078</td><td>2.4</td><td>6.6</td><td>0.007</td><td>0.9</td><td>8.2</td><td>23.9</td></tr><tr><td>NFU</td><td>Relative to 10 years</td><td></td><td></td><td>0.15</td><td>6.0</td><td>10.0</td><td>0.10</td><td>3.0</td><td>9.0</td><td>30.0</td></tr></tbody></table>

has become significantly higher in GWS and BIOW treatments compared to CN since 2004 and since 2011 in the case of MSW and FYM treatments. In 2015, SOC content was significantly different among treatments and increased as follows: CN < MSW < FYM < BIO = GWS (Table 4). In 2015, the amended treatments always presented significant increases compared to the CN for the soil pH (up to 1 unit for BIOW and MSW), the organic carbon content (up to + 69% for GWS and + 65% for BIOW), and the cation exchange capacity (up to + 41% for BIOW). Those soil properties are known to influence the availability and mobility of TE, thus their bioavailability and toxicity (Harmsen 2007; Smolders et al. 2009; Smolders et al. 2012).

Figure 3 summarizes the evolution of TE content in the topsoil layer over the period 1998–2015 (detailed statistics are available in Appendix C):

- For Cd, a slight and generally not significant decrease occurred in the control treatment (– 17% between 1998 and 2015), while significant differences have been observed between OWP treatments and CN since 2011.

- A large increase of all Cr values appeared in 2009 compared to previous years probably coming from a change in the analytical methods in the laboratory. This made difficult to analyse the evolution in Cr content in the different treatments. No lasting significant differences appeared among treatments.

- Cu became significantly higher in the amended treatments compared to CN in 2004, namely after 3 OWP applications, with no significant differences between amended treatments until 2011, when the plots receiving

the GWS compost became significantly above the other amended plots. Significant increase with time was observed for the soils receiving OWP, up to + 69% for GWS.

- For Hg, despite the low contents and relatively high standard errors compared to other elements, tendencies were observed with decrease observed in CN, BIOW and FYM treatments; significant differences have been observed with MSW and GWS > CN in 2004, 2006, 2007, 2011 and 2015 (only for GWS).

- For Ni, no significant differences were ever observed between the treatments and with years. Nevertheless, a slight but not significant increase was observed for the BIOW compost over the period 1998–2015.

- With respect to Pb, relatively large standard errors have been observed, since the beginning of the experiment, which can be due to initial spatial heterogeneity. Nevertheless, significant differences between all organic treatments and CN appeared only in 2015.

- Finally Zn content has become significantly larger in most organic treatments compared to CN in 2006, with the largest contents observed in the GWS treatment since 2013. Zn content increased significantly with time in soils receiving GWS (+ 32%), BIOW (+ 20%) and MSW (+ 15%).

In summary, the most conclusive significant differences of TE contents in the OWP compared to CN treatment have been observed: since 2004 for Cu and since 2006 for Zn due to significant increase with time in soils receiving OWP, since
<table><caption>**Table 4** Chemical properties and total trace element contents measured in 2015 in the topsoil layer (0–25 cm, sieved < 2 mm), compared to the average French and European trace element contents, with the following: mean values ± standard deviation of the four replicates (Abc letters stand for significant difference between treatments); co-compost of sewage sludge and green waste (GWS), biowaste compost (BIOW), compost of residual municipal solid waste (MSW), farmyard manure (FYM), no organic amendment (CN); na for not available</caption><thead><tr><th>Treat.</th><th>pH</th><th>Org. C</th><th>Total N</th><th>Olsen P<sub>2</sub>O<sub>5</sub></th><th>CEC</th><th>Cd</th><th>Cr</th><th>Cu</th><th>Hg</th><th>Ni</th><th>Pb</th><th>Zn</th></tr><tr><td></td><td></td><td>g kg<sup>-1</sup> DM</td><td></td><td></td><td>cmol<sup>+</sup> kg<sup>-1</sup> DM</td><td>mg kg<sup>-1</sup> DM</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></thead><tbody><tr><td>GWS</td><td>7.0 ± 0.0 (c)</td><td>16.9 ± 0.9 (a)</td><td>1.7 ± 0.1 (a)</td><td>0.20 ± 0.02 (b)</td><td>10.7 ± 0.4 (b)</td><td>0.25 ± 0.01 (ab)</td><td>73.4 ± 10.2 (ab)</td><td>19.9 ± 1.1 (a)</td><td>0.12 ± 0.02 (a)</td><td>15.6 ± 0.2 (a)</td><td>27.3 ± 5.5 (a)</td><td>68.2 ± 2.2 (b)</td></tr><tr><td>BIOW</td><td>7.8 ± 0.1 (a)</td><td>16.5 ± 0.3 (a)</td><td>1.6 ± 0.0 (a)</td><td>0.10 ± 0.01 (ab)</td><td>12.0 ± 0.4 (a)</td><td>0.25 ± 0.04 (b)</td><td>75.1 ± 7.5 (b)</td><td>15.0 ± 0.9 (b)</td><td>0.08 ± 0.01 (bc)</td><td>15.9 ± 0.8 (a)</td><td>26.1 ± 2.1 (a)</td><td>60.5 ± 1.1 (ab)</td></tr><tr><td>MSW</td><td>7.8 ± 0.1 (a)</td><td>13.2 ± 0.2 (c)</td><td>1.3 ± 0.0 (a)</td><td>0.07 ± 0.00 (ab)</td><td>10.5 ± 0.5 (b)</td><td>0.23 ± 0.01 (ab)</td><td>74.3 ± 4.7 (ab)</td><td>16.2 ± 0.6 (b)</td><td>0.10 ± 0.02 (bc)</td><td>15.5 ± 0.6 (a)</td><td>27.7 ± 4.1 (a)</td><td>61.3 ± 2.7 (ab)</td></tr><tr><td>FYM</td><td>7.4 ± 0.0 (b)</td><td>14.8 ± 0.4 (b)</td><td>1.4 ± 0.0 (a)</td><td>0.11 ± 0.01 (ab)</td><td>10.5 ± 0.4 (b)</td><td>0.24 ± 0.01 (ab)</td><td>74.1 ± 3.1 (ab)</td><td>15.5 ± 1.1 (b)</td><td>0.08 ± 0.02 (bc)</td><td>15.5 ± 0.8 (a)</td><td>25.3 ± 0.9 (a)</td><td>61.1 ± 5.3 (ab)</td></tr><tr><td>CN</td><td>6.8 ± 0.2 (d)</td><td>10.0 ± 0.3 (d)</td><td>1.0 ± 0.0 (a)</td><td>0.06 ± 0.01 (a)</td><td>8.5 ± 0.6 (c)</td><td>0.20 ± 0.02 (a)</td><td>72.5 ± 8.0 (a)</td><td>12.0 ± 0.8 (c)</td><td>0.07 ± 0.01 (c)</td><td>15.3 ± 0.9 (a)</td><td>21.3 ± 0.9 (b)</td><td>49.3 ± 2.1 (a)</td></tr><tr><td>Average contents in EU<sup>a</sup></td><td></td><td></td><td></td><td></td><td></td><td>0.09</td><td>21.72</td><td>13.01</td><td>0.04</td><td>18.36</td><td>15.3</td><td>na</td></tr><tr><td>Median contents in France<sup>b</sup></td><td></td><td></td><td></td><td></td><td></td><td>0.28</td><td>38.3</td><td>13.3</td><td>0.046</td><td>19.5</td><td>21.7</td><td>56.4</td></tr></tbody></table>

2011 for Cd and 2015 for Pb as related to apparent decrease in the control treatment. Losses of TE applied to soil in OWP have been reported in previous studies, as a relation to a combination of small losses from leaching, erosion plus unexplained experimental errors and movement due to tillage operations (McGrath and Lane 1989; Berti and Jacobs 1998). This could explain decrease for Cd and Pb observed in the present study.

In 2015, the contents of Cu and Zn have increased with respect to 1998. As observed in the present study and as already stressed by Smith (2009), composts and sewage sludges contain larger contents of TE than the background values found in soil (Tables 2 and 4), so that their regular application to land gradually raises the total TE content of soil in the long term. These results confirmed the numerous previous studies that showed significant accumulation of Zn, Cu, and less often Cd and Pb in soil, after repeated applications of biosolids (Moolenaar and Beltramy 1998; Richards et al. 1998; Warman and Termeer 2005; Kirchman et al. 2017), or composts (Gigliotti et al. 1996; Delschen 1999; Weber et al. 2007; Alvarenga et al. 2017). McBride and Spiers (2001) showed that sewage sludge can lead to much higher soil accumulations of TE (e.g. Cu, Zn) compared to FYM and mineral fertilization. Our results also confirmed that TE content increases are slow and slight. Thus, long-term evaluations are required to evaluate precisely the effects of repeated applications of OWP on TE accumulation in soils, especially for some toxic TE such as Pb and Cd. Indeed short-term studies usually do not reveal TE increases in soils (Douglas et al. 2003).

To complete the Fig. 3, Table 4 presents average contents found after 9 applications of OWP in 2015 for the 5 treatments, as compared to the average contents found in the topsoil layer in Europe (Toth et al. 2016) and the median contents found in France (BDETM 2019). When comparing to contents found in Europe and France, TE contents measured in 2015 were globally in the range or slightly larger than usual contents found in Europe and France, as follows:

- Cu, Hg, Pb and Zn, in the range or slightly larger than average contents found in Europe and France, knowing that some agricultural practices are known to cause larger Cu concentrations in soils due to inputs of Cu-based pesticides used in vineyard and orchards (Belon et al. 2012).

- Cd, larger than average contents found in Europe and in the range of median contents found in France.

- Ni, slightly lower than average contents found in Europe and France.

## Effect of organic waste product applications on plant trace element contents

Table 5 summarizes the average contents of TE in wheat, maize and barley grains. Yields are also presented. Statistical

<footnote><sup>a</sup> Cited in Toth et al. (2016)</footnote>

<footnote><sup>b</sup> BDETM (https://www.gissol.fr/le-gis/programmes/base-de-donnees-elements-traces-metaliques-bdetm-65, Accessed in 18 July 2019)</footnote>
<figure><img src="image_7.png" /><figcaption>**Fig. 3** Evolution of trace element contents in the topsoil layer over the period 1998–2015 for Cd, Cr, Cu, Hg, Ni, Pb and Zn (mg kg⁻¹ DM; DM for dry matter). Diamonds stand for co-compost of sewage sludge and green waste (GWS), squares stand for biowaste compost (BIOW), triangles stand for compost of residual municipal solid waste (MSW), circles stand for farmyard manure (FYM) and crosses stand for control (CN). The error bars are standard deviations</figcaption></figure>

analyses are detailed in Appendix I (Electronic supplementary data). The TE contents decreased in the order Zn > Cu > Cr > Ni. In addition, for toxic TE, i.e. Cd, Cr, Hg, Ni and Pb, most values were lower than the quantification limit, for the 3 species. A few significant differences appeared between the organic treatments with respect to CN or among the organic

treatments, but not repeatedly in time. Compost application (GWS, BIOW, MSW) did not induce significant increase of TE content in grains for the 3 species studied. Indeed large variations between years were observed, as shown by large standard deviation. This finding was already noticed by Chaudri et al. (2007) who found annual differences in grain
Cd uptake. Indeed agroclimatic conditions and varietal differ-
ences are important factors in the uptake of Cd by crops. For
Cu and Zn, essential elements for plant nutrition, values were
above the quantification limit. Average grain contents of TE
found in QualiAgro were globally equal or below regulatory
(Codex Alimentarius 2015) and published values for Cd, Hg,
Ni and Pb; while contents of Cr, Cu and Zn were in the range
of published values for grains (Kabata-Pendias and Pendias
2001; Feix and Tremel 2005) (Table 5).

These results confirmed that there was few or no effect of
repeated application of OWP on TE contents in cereal grains
(Businelli et al. 1996; Gigliotti et al. 1996; Smith 2009;
Alvarenga et al. 2017). Plant physiology tends to regulate
TE in the shoots and the grains so that most elements do not
reach levels which could cause even chronic toxicity to plants
or consumers (McBride 1995; Chaney 2012). The exception
is for some elements, e.g. Cd, which are readily absorbed and
translocated to shoots without causing phytotoxicity, and can
be transferred in the food chain (Chaney 2012). Another rea-
son explaining the lack of differences between treatments is
the zone of root uptake which is often deeper than the plough
layer (McBride 1995).

Mass balances

The mass balance approach is useful to evaluate the sustain-
ability of a system, and to anticipate depletion and accumula-
tion of TE (Moolenaar et al. 1997b; Yu et al. 2011). Mass
balance calculations can thus be used to assess the impacts
of cropping practices, including spreading of OWP, on TE
contents in soils (i.e. accumulation, steady-state, or depletion).
It considers the major inputs and outputs. In the case of
QualiAgro, we considered atmospheric depositions and
OWP inputs as inputs, and crop uptake and leaching as out-
puts. In the present calculation, the impact of surface runoff
and erosion was not considered. Inputs of TE via pesticides,
mineral N fertilizers and organic fertilizers were also
neglected. All calculations are presented in Electronic supple-
mentary material.

Field mass balance for the period 1998–2015 is presented
in Fig. 4 (Appendix J in Electronic supplementary material),
with cumulated fluxes for inputs and outputs. Input of TE by
OWP were much larger compared to atmospheric inputs and
outputs (i.e. crop offtake and leaching). It represented up to
17% of the initial soil stock for Cd (MSW and GWS), 56% for
Cu (GWS), 34% for Hg (GWS), 8% for Ni (BIOW), 17% for
Pb (BIOW) and 30% for Zn (GWS) (Table 6). Among OWP,
FYM most often presented the smallest TE mass balances
while GWS and MSW presented the largest ones.

Outputs were generally very low. For Cr, Pb and Zn, the
percentages of leaching were lower than 2% of TE OWP
inputs. Leaching represented respectively for Cd, Cu and Ni:
up to 2.3% of Cd inputs flux by GWS, 2.7% of Cu inputs by

FYM and 7.8% of Ni inputs by FYM. For Cu and Zn, for
which TE content was higher than the quantification limits
in grains and residues, crop uptake represented 1.8% (GWS)
to 4.4% (BIOW) of the Cu inputs by OWP; and 5.1% (GWS)
to 7.3% (BIOW) for Zn. For the other TE, for which [TE]crop
in grains was globally lower than the quantification limit, crop
uptake has been assessed to represent less than 2% of TE
inputs by OWP (except for Hg in BIOW).

For all OWP amendments, the resulting positive mass bal-
ance (inputs – outputs) was similar to the input flux by OWP
for Cd, Hg, Ni and Pb (Table 6). For Cu and Zn, the calculated
outputs resulted in slight differences between the positive
mass balance and the input flux by OWP (respectively – 0.2
to 0.3 kg ha⁻¹ and –0.7 to 0.9 kg ha⁻¹). Globally over the
period 1998–2015, mass balance corresponded to the delta
soil stock for Cu, Ni and Zn (Table 6; Appendix E,
Table E.2). For Cd and Hg, mass balance was larger than the
delta soil stock, suggesting losses such as leaching for Cd or
evaporation for Hg (Cobbett and Van Heyst 2004). For Cd,
content in the topsoil layer remained around 0.24 mg kg⁻¹ in
OWP treatments in 1998–2015, while it decreased to 0.20 mg
kg⁻¹ in the control treatment in 2015 (Tables 1 and 4). Putative
leaching of Cd observed in the control treatment could occur
for each treatment. In OWP treatments this could explain the
stability observed in the topsoil layer Cd content in 1998–
2015 despite slight positive mass balance (Table 6). For Pb,
mass balance was also larger than the delta soil stock.
However, there were large variabilities in soil contents so that
it was difficult to conclude. In contrary to amended plots,
negative delta soil stocks have been measured in the control
treatment. This depletion in control soils was as follows: Pb >
Zn > Cu > Cd > Hg. The calculated mass balance in this
control treatment was less negative than the measured de-
creased in soil stocks. Such differences could be due to over-
estimation of atmospheric inputs or underestimation of
leaching outputs, since these fluxes have not been measured
on our site or not during the whole period in case of leaching.

Prospective evaluation: threshold values and “very”
long-term simulation of total soil TE content

As said by Lofts et al. (2007): “At what point does this con-
tamination become pollution (contamination that results in
adverse biological effects)?” We tried to answer that question
at first by (i) estimating threshold values considered as soil
safe limits, and then by (ii) simulating the long-term accumu-
lation of TE in soils exposed to a precise and reasonable sce-
nario of OWP amendments.

In the present study, as explained in “Threshold contents
for trace elements in soil” section, we used the “Threshold
calculator for metals in soils” (Oorts et al. 2018) to estimate
threshold values considered as PNEC. The “Threshold calcu-
lator” is based on species sensitivity distributions (SSD) found
<table><caption>**Table 5** Average yield (t DM ha⁻¹; DM for dry matter) and average trace element contents in consumed exported grains of wheat, barley and maize over the period 1998–2015 for Cd, Cr, Cu, Hg, Ni, Pb and Zn (mg kg⁻¹ DM), with the following: mean values ± standard deviation which stands for the year variation; co-compost of sewage sludge and green waste (GWS), biowaste compost (BIOW), compost of residual municipal solid waste (MSW), farmyard manure (FYM), no organic amendment (CN)</caption><thead><tr><th>Treatment</th><th>Species</th><th>Yield</th><th>Cd</th><th>Cr</th><th>Cu</th><th>Hg</th><th>Ni</th><th>Pb</th><th>Zn</th></tr><tr><td></td><td></td><td>t DM ha<sup>-1</sup></td><td colspan="6">mg kg<sup>-1</sup> DM</td></tr></thead><tbody><tr><td>GWS</td><td>Wheat 1998–2013</td><td>8.0 ± 0.7</td><td>0.015 ± 0.003 (<sup>«</sup>)</td><td>0.22 ± 0.13 (<sup>«</sup>)</td><td>3.0 ± 0.5</td><td>0.003 ± 0.004 (<sup>«</sup>)</td><td>0.11 ± 0.08 (<sup>«</sup>)</td><td>0.12 ± 0.08 (<sup>«</sup>)</td><td>20.0 ± 3.9</td></tr><tr><td>BIOW</td><td>Wheat 1998–2013</td><td>7.9 ± 0.6</td><td>0.014 ± 0.006 (<sup>«</sup>)</td><td>0.26 ± 0.16 (<sup>«</sup>)</td><td>3.3 ± 0.4</td><td>0.003 ± 0.004 (<sup>«</sup>)</td><td>0.12 ± 0.05 (<sup>«</sup>)</td><td>0.12 ± 0.09 (<sup>«</sup>)</td><td>20.4 ± 3.7</td></tr><tr><td>MSW</td><td>Wheat 1998–2013</td><td>7.9 ± 0.7</td><td>0.016 ± 0.005 (<sup>«</sup>)</td><td>0.23 ± 0.15 (<sup>«</sup>)</td><td>3.5 ± 0.5</td><td>0.003 ± 0.005 (<sup>«</sup>)</td><td>0.14 ± 0.05 (<sup>«</sup>)</td><td>0.16 ± 0.09 (<sup>«</sup>)</td><td>21.5 ± 3.3</td></tr><tr><td>FYM</td><td>Wheat 1998–2013</td><td>8.1 ± 0.7</td><td>0.018 ± 0.005 (<sup>«</sup>)</td><td>0.21 ± 0.10 (<sup>«</sup>)</td><td>3.2 ± 0.5</td><td>0.003 ± 0.004 (<sup>«</sup>)</td><td>0.13 ± 0.08 (<sup>«</sup>)</td><td>0.15 ± 0.10 (<sup>«</sup>)</td><td>21.3 ± 3.4</td></tr><tr><td>CN</td><td>Wheat 1998–2013</td><td>7.2 ± 0.7</td><td>0.020 ± 0.009 (<sup>«</sup>)</td><td>0.22 ± 0.12 (<sup>«</sup>)</td><td>3.3 ± 0.5</td><td>0.002 ± 0.004 (<sup>«</sup>)</td><td>0.21 ± 0.11 (<sup>«</sup>)</td><td>0.14 ± 0.09 (<sup>«</sup>)</td><td>20.0 ± 3.6</td></tr><tr><td>GWS</td><td>Barley 2007 & 2015</td><td>7.2 ± 0.3</td><td>0.003 ± 0.002 (<sup>«</sup>)</td><td>0.19 ± 0.27 (<sup>«</sup>)</td><td>4.7 ± 1.1</td><td>0.001 ± 0.001 (<sup>«</sup>)</td><td>0.10 ± 0.13 (<sup>«</sup>)</td><td>0.06 ± 0.03 (<sup>«</sup>)</td><td>29.9 ± 2.5</td></tr><tr><td>BIOW</td><td>Barley 2007 & 2015</td><td>6.7 ± 0.3</td><td>0.003 ± 0.002 (<sup>«</sup>)</td><td>0.17 ± 0.24 (<sup>«</sup>)</td><td>4.8 ± 1.0</td><td>0.001 ± 0.001 (<sup>«</sup>)</td><td>0.06 ± 0.04 (<sup>«</sup>)</td><td>0.05 ± 0.02 (<sup>«</sup>)</td><td>27.9 ± 0.0</td></tr><tr><td>MSW</td><td>Barley 2007 & 2015</td><td>6.6 ± 1.1</td><td>0.006 ± 0.001 (<sup>«</sup>)</td><td>0.19 ± 0.26 (<sup>«</sup>)</td><td>5.0 ± 1.2</td><td>0.001 ± 0.001 (<sup>«</sup>)</td><td>0.08 ± 0.05 (<sup>«</sup>)</td><td>0.06 ± 0.01 (<sup>«</sup>)</td><td>29.7 ± 0.7</td></tr><tr><td>FYM</td><td>Barley 2007 & 2015</td><td>6.8 ± 1.3</td><td>0.005 ± 0.002 (<sup>«</sup>)</td><td>0.18 ± 0.26 (<sup>«</sup>)</td><td>5.2 ± 1.6</td><td>0.001 ± 0.000 (<sup>«</sup>)</td><td>0.09 ± 0.05 (<sup>«</sup>)</td><td>0.08 ± 0.03 (<sup>«</sup>)</td><td>28.7 ± 2.4</td></tr><tr><td>CN</td><td>Barley 2007 & 2015</td><td>5.7 ± 1.0</td><td>0.007 ± 0.003 (<sup>«</sup>)</td><td>0.18 ± 0.22 (<sup>«</sup>)</td><td>4.4 ± 0.6</td><td>0.001 ± 0.000 (<sup>«</sup>)</td><td>0.08 ± 0.00 (<sup>«</sup>)</td><td>0.25 ± 0.22 (<sup>«</sup>)</td><td>25.3 ± 1.7</td></tr><tr><td>GWS</td><td>Maize 1999–2014</td><td>9.1 ± 0.7</td><td>0.005 ± 0.006 (<sup>«</sup>)</td><td>0.07 ± 0.07 (<sup>«</sup>)</td><td>2.0 ± 0.5</td><td>0.0002 ± 0.0004 (<sup>«</sup>)</td><td>0.28 ± 0.26 (<sup>«</sup>)</td><td>0.05 ± 0.05 (<sup>«</sup>)</td><td>17.8 ± 2.4</td></tr><tr><td>BIOW</td><td>Maize 1999–2014</td><td>9.1 ± 0.7</td><td>0.006 ± 0.007 (<sup>«</sup>)</td><td>0.08 ± 0.06 (<sup>«</sup>)</td><td>2.2 ± 0.4</td><td>0.0002 ± 0.0002 (<sup>«</sup>)</td><td>0.20 ± 0.11 (<sup>«</sup>)</td><td>0.06 ± 0.06 (<sup>«</sup>)</td><td>19.2 ± 2.5</td></tr><tr><td>MSW</td><td>Maize 1999–2014</td><td>8.9 ± 1.0</td><td>0.008 ± 0.008 (<sup>«</sup>)</td><td>0.11 ± 0.08 (<sup>«</sup>)</td><td>2.2 ± 0.5</td><td>0.0003 ± 0.0004 (<sup>«</sup>)</td><td>0.22 ± 0.11 (<sup>«</sup>)</td><td>0.06 ± 0.04 (<sup>«</sup>)</td><td>19.9 ± 3.2</td></tr><tr><td>FYM</td><td>Maize 1999–2014</td><td>9.3 ± 1.0</td><td>0.006 ± 0.005 (<sup>«</sup>)</td><td>0.08 ± 0.06 (<sup>«</sup>)</td><td>2.1 ± 0.5</td><td>0.0002 ± 0.0003 (<sup>«</sup>)</td><td>0.22 ± 0.11 (<sup>«</sup>)</td><td>0.06 ± 0.05 (<sup>«</sup>)</td><td>19.2 ± 3.4</td></tr><tr><td>CN</td><td>Maize 1999–2014</td><td>8.1 ± 1.9</td><td>0.006 ± 0.005 (<sup>«</sup>)</td><td>0.09 ± 0.08 (<sup>«</sup>)</td><td>2.2 ± 0.5</td><td>0.0004 ± 0.0006 (<sup>«</sup>)</td><td>0.26 ± 0.09 (<sup>«</sup>)</td><td>0.10 ± 0.10 (<sup>«</sup>)</td><td>19.1 ± 3.5</td></tr><tr><td>Feix and Tremel 2005</td><td>Barley grain</td><td>0.02</td><td>0.02</td><td>0.02</td><td>5.0</td><td>0.01</td><td>0.20</td><td>0.45</td><td>28</td></tr><tr><td>Kabata-Pendias and Pendias 2001 Codex</td><td>Wheat grain</td><td>0.03</td><td>0.10</td><td>4.5</td><td>4.5</td><td>0.01</td><td>0.30</td><td>0.50</td><td>25</td></tr><tr><td>Alimentarius 2015<sup>a</sup></td><td>Grain contents</td><td>0.02–0.04</td><td>0.2</td><td>1.8–10</td><td></td><td></td><td>0.14–0.90</td><td>0.18–0.40</td><td>23–37</td></tr><tr><td></td><td>Cereal grains</td><td>0.1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Wheat</td><td>0.2</td><td></td><td></td><td></td><td></td><td></td><td>0.2</td><td></td></tr></tbody></table>

(«) at least half of values were inferior to the quantification limit

a Contents expressed as milligram per kilogram of fresh matter

<figure><img src="image_8.png" /></figure>
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>
► **Fig. 4** Mass balance of fluxes between compartments over the period 1998–2015, for Cd, Cr, Cu, Hg, Ni, Pb and Zn (kg ha⁻¹). With co-compost of sewage sludge and green waste (GWS), biowaste compost (BIOW), compost of residual municipal solid waste (MSW), farmyard manure (FYM) and control (CN). Inputs fluxes stand for atmospheric deposits and OWP inputs (respectively bold hatchings at the top and fine hatchings at the bottom of input), output fluxes stand for plant exportations and water drainage (respectively zigzag at the top and little squares at the bottom of output)

in the literature and applied to estimate hazardous concentration (HC), taking into account influencing soil physicochemical properties. In the present study, hazardous concentration for 5% of the species (HC₅) derived from EC₁₀ distributions (Effective concentration for which 10% of the community is affected by the concentration of the trace element) has been chosen as it is proposed as a PNEC and “safe” concentration (Oorts 2018).

The results of these calculations, using soil parameters measured in 2015, are presented in Table 7. After 9 applications of OWP, the TE contents in the topsoil layer in 2015 (Table 4) were below the HC₅ of EC₁₀, or PNEC values, for Cd, Cu, Ni, Pb and Zn for the 5 treatments. It can be noticed that PNEC values for the control treatment were always below those found for OWP amended soils, except for Cd whose PNEC are similar in all treatments. This could be related to the so-called sludge protection mechanisms, also valid for all organic amendments (sewage sludge, manure, compost) as shown by Smolders et al. (2012). Addition of organic matter

can have protective effects by decreasing metal bio/availability, through increasing soil retention factors, including pH, soil organic matter and CEC.

The quotients of total TE contents divided by the soil PNEC values stand in the range 0.07–0.09 for Cd, 0.19–0.26 for Cu, 0.27–0.38 for Ni, 0.10–0.12 for Pb and 0.36–0.48 for Zn. They suggest that TE contents in soils receiving repeated applications of OWP would be safe for a long time for Cd, Pb, even for Cu and Ni. For Zn, calculations showed increasing risks associated with the contents of Zn in soils which tend towards the protective threshold in the following order: GWS > CN > FYM > MSW > BIOW. Indeed Zn is the element in OWP-treated soils identified as the main concern in relation to eco-toxicity and long-term soil fertility (Smith 2009; Yoshida et al. 2018).

Results confirmed that repeated application of OWP, even in large amounts such as those applied in the QualiAgro experiment, did not alter the soil quality in terms of TE contents since they remained far below the calculated protective threshold value for Cd, Cu, Ni, Pb and Zn. However, based on the observation of the correspondence between the balance approach and the evolution of TE stocks in the topsoil layer during 17 years, we aimed at completing our evaluation of such practices at a longer time scale, i.e. about 100 years.

The evolution of the total TE concentrations in topsoil layer has thus been simulated until 2100 using the Eq. (8) (“Simulation of soil trace element content” section) and the average measured inputs and outputs summarized in

<table><caption>**Table 6** Cumulated input through organic waste product spreading (Input OWP flux), soil stock variation (Delta stock) and mass balance (Balance inputs – outputs) of trace elements over the period 1998–2015 for Cd, Cr, Cu, Hg, Ni, Pb and Zn (kg ha⁻¹), with the following: the percentage of the considered flux with respect to the initial soil stock</caption><thead><tr><th></th><th></th><th>Cd</th><th>Cu</th><th>Hg</th><th>Ni</th><th>Pb</th><th>Zn</th></tr><tr><th></th><th></th><th colspan="6">kg ha<sup>-1</sup></th></tr></thead><tbody><tr><td rowspan="3">GWS</td><td>Input OWP flux</td><td>0.15 (17.0%)</td><td>25.6 (56.1%)</td><td>0.12 (34.0%)</td><td>4.0 (7.1%)</td><td>9.1 (8.4%)</td><td>59.6 (30.4%)</td></tr><tr><td>Delta stock</td><td>0.05 (5.2%)</td><td>29.2 (63.9%)</td><td>0.09 (25.8%)</td><td>3.0 (5.3%)</td><td>-5.0 (-4.6%)</td><td>60.8 (31.0%)</td></tr><tr><td>Balance inputs – outputs</td><td>0.15 (16.8%)</td><td>25.4 (55.8%)</td><td>0.12 (33.9%)</td><td>4.0 (7.0%)</td><td>9.1 (8.4%)</td><td>58.7 (30.0%)</td></tr><tr><td rowspan="3">BIOW</td><td>Input OWP flux</td><td>0.13 (15.1%)</td><td>10.9 (23.3%)</td><td>0.03 (8.1%)</td><td>4.8 (8.5%)</td><td>15.9 (17.4%)</td><td>42.3 (22.1%)</td></tr><tr><td>Delta stock</td><td>0.07 (8.6%)</td><td>10.0 (21.4%)</td><td>-0.06 (-15.1%)</td><td>3.7 (6.4%)</td><td>8.0 (8.7%)</td><td>38.1 (19.9%)</td></tr><tr><td>Balance inputs – outputs</td><td>0.13 (15.2%)</td><td>10.7 (22.9%)</td><td>0.03 (7.6%)</td><td>4.8 (8.4%)</td><td>15.9 (17.4%)</td><td>41.6 (21.7%)</td></tr><tr><td rowspan="3">MSW</td><td>Input OWP flux</td><td>0.15 (17.1%)</td><td>16.0 (35.2%)</td><td>0.09 (25.5%)</td><td>3.4 (6%)</td><td>15.5 (16.7%)</td><td>45.3 (22.6%)</td></tr><tr><td>Delta stock</td><td>-0.005 (-0.6%)</td><td>15.4 (33.8%)</td><td>0.02 (4.7%)</td><td>2.5 (4.5%)</td><td>11.4 (12.3%)</td><td>30.5 (15.2%)</td></tr><tr><td>Balance inputs – outputs</td><td>0.15 (17.1%)</td><td>15.7 (34.5%)</td><td>0.1 (25.6%)</td><td>3.3 (5.9%)</td><td>15.5 (16.7%)</td><td>44.4 (22.2%)</td></tr><tr><td rowspan="3">FYM</td><td>Input OWP flux</td><td>0.13 (14.8%)</td><td>11.3 (23.4%)</td><td>0.01 (2.8%)</td><td>1.5 (2.6%)</td><td>13.9 (15.0%)</td><td>40.6 (20.7%)</td></tr><tr><td>Delta stock</td><td>0.022 (2.5%)</td><td>10.3 (21.3%)</td><td>-0.05 (-13.8%)</td><td>2.2 (3.9%)</td><td>3.3 (3.5%)</td><td>35.5 (18.1%)</td></tr><tr><td>Balance inputs – outputs</td><td>0.13 (14.9%)</td><td>11.0 (22.8%)</td><td>0.01 (2.6%)</td><td>1.4 (2.5%)</td><td>13.9 (15.0%)</td><td>39.7 (20.3%)</td></tr><tr><td rowspan="3">CN</td><td>Input OWP flux</td><td>nc</td><td>nc</td><td>nc</td><td>nc</td><td>nc</td><td>nc</td></tr><tr><td>Delta stock</td><td>-0.11 (-12.5%)</td><td>-0.6 (-1.2%)</td><td>-0.09 (-26.7%)</td><td>0.1 (0.2%)</td><td>-7.1 (-8.1%)</td><td>-5.7 (-3%)</td></tr><tr><td>Balance inputs – outputs</td><td>-0.0001 (0.0%)</td><td>-0.05 (-0.1%)</td><td>0.0003 (0.1%)</td><td>-0.01 (0.0%)</td><td>0.002 (0.0%)</td><td>-0.5 (-0.3%)</td></tr></tbody></table>

written into brackets; co-compost of sewage sludge and green waste (GWS), biowaste compost (BIOW), compost of residual municipal solid waste (MSW), farmyard manure (FYM), no organic amendment (CN); nc stands for not concerned
<table><caption>**Table 7** Threshold values carried out by the threshold calculator for the 5% probability level (5% Hazardous concentration or HC₅) of the reliable effective concentration EC₁₀ (i.e. effective concentration for which 10% of the community is affected by the concentration of the trace element), with the following: mean values per treatment ± standard deviation of the four replicates (except for Cd for which the calculator calculations were not available for raw data); co-compost of sewage sludge and green waste (GWS), biowaste compost (BIOW), compost of residual municipal solid waste (MSW), farmyard manure (FYM), no organic amendment (CN)</caption><thead><tr><th>Metal</th><th>Treatment</th><th>HC<sub>5</sub> of EC<sub>10</sub> values</th></tr><tr><td></td><td></td><td>mg kg<sup>-1</sup> DM</td></tr></thead><tbody><tr><td rowspan="5">Cd</td><td>GWS</td><td>2.7</td></tr><tr><td>BIO</td><td>2.7</td></tr><tr><td>MSW</td><td>2.7</td></tr><tr><td>FYM</td><td>2.7</td></tr><tr><td>CN</td><td>2.7</td></tr><tr><td rowspan="5">Cu</td><td>GWS</td><td>77.9 ± 2.0</td></tr><tr><td>BIO</td><td>75.2 ± 1.3</td></tr><tr><td>MSW</td><td>67.9 ± 1.1</td></tr><tr><td>FYM</td><td>72.5 ± 1.5</td></tr><tr><td>CN</td><td>64.2 ± 2.0</td></tr><tr><td rowspan="5">Ni</td><td>GWS</td><td>52.5 ± 2.5</td></tr><tr><td>BIO</td><td>59.4 ± 2.3</td></tr><tr><td>MSW</td><td>51.1 ± 2.6</td></tr><tr><td>FYM</td><td>51.2 ± 2.2</td></tr><tr><td>CN</td><td>40.6 ± 3.3</td></tr><tr><td rowspan="5">Pb</td><td>GWS</td><td>240.7 ± 11.7</td></tr><tr><td>BIO</td><td>273.3 ± 10.4</td></tr><tr><td>MSW</td><td>234.4 ± 12.6</td></tr><tr><td>FYM</td><td>234.9 ± 10.8</td></tr><tr><td>CN</td><td>182.5 ± 16.8</td></tr><tr><td rowspan="5">Zn</td><td>GWS</td><td>142.0 ± 4.6</td></tr><tr><td>BIO</td><td>166.3 ± 4.6</td></tr><tr><td>MSW</td><td>151.7 ± 5.8</td></tr><tr><td>FYM</td><td>146.4 ± 4.5</td></tr><tr><td>CN</td><td>116.8 ± 9.1</td></tr></tbody></table>

Appendix J, considering approaching usual agronomic spreading amounts of OWP (i.e. 1.5 lower than amounts of OWP applied in the QualiAgro device). This scenario would be considered as an extreme local scenario in which OWP would be applied during a period of 100 years at an average rate of 30 t ha⁻¹ per 3 years for farmyard manure and 20–30 tons ha⁻¹ per 3 years for composts. In practice, this extreme local scenario would not occur due to the need of maintaining adequate balance of macronutrients in the soils. According to Moolenaar et al. (1997b), the time scale of one century was considered here, since calculations for longer periods should be avoided. Calculations presented in Appendix K showed that the model gave good predictions of [TE]top,i,2015 for Cr, Cu, Ni and Zn (i.e. deviations < 7% in absolute value), acceptable prediction for Cd and Pb (< 21% of deviation) and no

accuracy for Hg which presented large variation in results obtained in QualiAgro, as presented previously. Differences between the model and measured data in QualiAgro could be partly explained by losses for Cd such as underestimated leaching measurements, variability in data for Pb and Hg, and volatilization of Hg. Other limits could also alter simulated results, such as neglecting erosion and using fixed average values, particularly for leaching TE.

The outputs of the calculation using usual agronomic amounts of OWP are presented in Fig. 5 for the evolution of the total content in the topsoil layer for Cd, Cr, Cu, Ni, Pb and Zn. Simulation results showed increases in concentrations in 2100 until 0.4 mg Cd kg⁻¹ DM (GWS, MSW), 56 mg Cr kg⁻¹ DM (MSW), 38 mg Cu kg⁻¹ DM (GWS), 20 mg Ni kg⁻¹ DM (BIOW), 42 mg Pb kg⁻¹ DM (BIOW, MSW) and 109 mg Zn kg⁻¹ DM (GWS). With simulations from 1998 to 2100 using an extreme local scenario, application of OWP in lands could lead to long-term accumulation of TE in the topsoil layer up to + 64% for Cd (MSW), + 22% for Cr (MSW), + 215% for Cu (GWS), + 33% for Ni (BIOW), + 63–65% for Pb (BIOW, MSW) and + 111% for Zn (GWS).

In the case of Cd, Sterckeman et al. (2018) found 3% of increase by the end of the next century for the average Cd content in French soils under annual crops, while Six and Smolders (2014) found a decrease in European soil Cd content for a 100-year period. In both studies, simulations of the soil Cd content was based on average inputs/outputs estimated at larger scale (i.e. sub-region, region), and on Cd concentration in leaching water predicted from soil pH, Cd and organic C content. Contrary to Sterckeman et al. (2018), fixed values for [Cd]lea and Qlea were used in the simulation model in the Eq. (8). The model used in the present study, validated with measured data, could be considered as an extreme local scenario, nevertheless possible. This could explain the larger increase in Cd contents found in the present study compared to the ones cited above. Belon et al. (2012) underlined the lack of output data for crop uptake and leaching. Any accurate simulation approach should also consider suitable inputs and outputs, as well as its outputs validation according to measured data, especially in long-term field experiments such as QualiAgro. In fact in the work of Sterckeman et al. (2018), long-term experiments (including QualiAgro and PROspective devices of the SOERE PRO network) were taken as reference to adjust leaching predictions.

Finally, trace element contents simulated in 2100 (i.e. [TE]top,i,2100) remained below the thresholds values of HC₅EC₁₀ (i.e. PNEC) presented in Table 7. Nevertheless, these results also showed that applications of OWP in usual agronomic amounts would increase the soil Zn content to levels close to the corresponding safe threshold. With TE contents in OWP conform to existing French regulations (NFU 44-051, 2006 and NFU 44-095, 2002), and if OWP are applied in appropriate amounts, repeated application of OWP
<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /><figcaption>**Fig. 5** Long-term simulation of the topsoil layer content of Cd, Cr, Cu, Ni, Pb and Zn (mg kg⁻¹ DM; DM for dry matter), considering the example of an extreme local scenario. Diamonds stand for co-compost of sewage sludge and green waste (GWS), squares stand for biowaste</figcaption></figure>

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /></figure>

<figure><img src="image_21.png" /><figcaption>compost (BIOW), triangles stand for compost of residual municipal solid waste (MSW), circles stand for farmyard manure (FYM) and crosses stand for control (CN)</figcaption></figure>

should not lead to critical contents in soils in the long term (Delschen 1999). Indeed, estimated increases of TE in the topsoil layer were lower than the TE content increase found by Charlton et al. (2016) to affect the microbial biomass in soil.

Conclusion and perspectives

OWP applications to lands offer valuable advantages. The BIOW compost could be considered as a good soil amendment, for acidic soil pH correction and soil organic matter

increase; the GSW compost can be used also as soil organic amendment but as well for phosphorus and nitrogen fertilization; the FYM could be considered as an organic amendment and fertilizer and the MSW compost could be used as source of organic matter implied in stimulation of biological activities and aggregate stability. However, TE contents in OWP are often larger than those found in agricultural soils. The GWS and MSW composts were the OWP with the largest TE contents, while the farmyard manure tended to have the lowest TE contents. In QualiAgro, repeated application of OWP led to significant accumulation of Zn and Cu in the topsoil layer, especially with the compost of green waste and sludges
(GWS). Nevertheless, TE contents in soils did not overpass calculated protective threshold values for Cd, Cu, Ni, Pb and Zn. Among OWP studied, the mass balance represented up to 17% of the initial soil stock for Cd (MSW and GWS), 56% for Cu (GWS), 34% for Hg (GWS), 8% for Ni (BIOW), 17% for Pb (BIOW) and 30% for Zn (GWS). Positive mass balances corresponded to the observed increases in delta soil stocks, especially for Cu and Zn. Weak outputs (i.e. leaching, crop offtakes) were observed compared to the OWP inputs. The QualiAgro data confirmed previous studies showing little or no effect of repeated application of OWP on TE contents in cereal grains. Values were below or in the range of usual TE contents found in grains. Simulation results using the average inputs/outputs measured in the QualiAgro field experiment showed increases in concentrations in 2100 up to 0.4 mg Cd kg⁻¹ DM (GWS, MSW), 56 mg Cr kg⁻¹ DM (MSW), 38 mg Cu kg⁻¹ DM (GWS), 20 mg Ni kg⁻¹ DM (BIOW), 42 mg Pb kg⁻¹ DM (BIOW, MSW) and 109 mg Zn kg⁻¹ DM (GWS).

The mass balance calculations pointed out the difficulty to collect complete data set from the beginning of the experimentation and the necessity to estimate some data. Complete monitoring and saving samples from the beginning of an experiment at same sampling periods for topsoil layer and the soil sub-layer would permit to avoid missing data.

The mass balance approach could be used to calculate input rates of TE through OWP application to avoid potential adverse effects of TE accumulation in soils. To maintain this state, regulation should consider lower TE contents in OWP and/or lower spreading amounts, based on studies such as the QualiAgro one. In the case of Zn, considering usual amounts of OWP (< 10 t DM ha⁻¹ per 2 years) and simulated evolution of soil contents, limiting input flux of Zn by OWP applied in lands would lead to avoid long-term Zn accumulation in the topsoil layer overpassing protective threshold. For instance, the regulatory limit of total Zn input flux could be fixed at a lower level than 30 kg ha⁻¹ per 10 years (e.g. at the level of 20 kg ha⁻¹ per 10 years).

Such study should be applied in the future in other field experiments integrated in the SOERE PRO (network of long-term experiments dedicated to the study of impacts of organic waste product recycling). This would (i) complete mass balances and references of average inputs/outputs that could be used in regional/national-scale evaluations and (ii) validate the simulation model used in the present study for various agropedo-climatic contexts and different OWP. In addition, for complete examination of the TE status in soils receiving repeated applications of OWP, the present study should be done for wider range of TE and completed by measurements of TE bioavailability.

**Acknowledgments** The QualiAgro field experiment forms part of the SOERE PRO (network of long-term experiments dedicated to the study of impacts of organic waste product recycling) certified by ALLENVI

(Alliance Nationale de Recherche pour l'Environnement) and integrated as a service of the "Investment for future" infrastructure AnaEE-France, overseen by the French National Research Agency (ANR-11-INBS-0001). We thank Koen Oorts for his help in calculations with the threshold calculator and editing the "Threshold contents for trace elements in soil" section.

**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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.

## References

Alvarenga P, Mourinha C, Farto M, Santos T, Palma P, Sengo J, Morais MC, Cunha-Queda C (2015) Sewage sludge, compost and other representative organic wastes as agricultural soil amendments: benefits versus limiting factors. Waste Manage 40:44–52

Alvarenga P, Palma P, Mourinha C, Farto M, Dores J, Patanita M, Cunha-Queda C, Natal-da-Luz T, Renaud M, Sousa JP (2017) Recycling organic wastes to agricultural land as a way to improve its quality: a field study to evaluate benefits and risks. Waste Manage 61:582–592

Azimi S, Cambier P, Lecuyer I, Thevenot D (2004) Heavy metal determination in atmospheric deposition and other fluxes in northern France agrosystems. Water Air Soil Poll 157:295–313

Baize D, Deslais W, Saby N (2007) Content of eight trace metals (Cd, Cr, Cu, Hg, Ni, Pb, Se, Zn) in French arable topsoils. ADEME. https://www.ademe.fr/sites/default/files/assets/documents/56610_etalchapitres_1-3.pdf. Accessed 18 July 2019 (in French)

BDETM (2019) French database of trace element contents in topsoil layers. https://www.gissol.fr/le-gis/programmes/base-de-donnees-elements-traces-metalliques-bdetm-65. Accessed 18 July 2019 (in French)

Bell A, Leclerc B (2015) Inventory of French field experiments devoted to the study of the effects of organic waste products from 1974 to 2012 (in French)

Belon E, Boisson M, Deportes IZ, Eglin TK, Feix I, Bispo AO, Galsomies L, Leblond S, Guellier CR (2012) An inventory of trace elements inputs to French agricultural soils. Sci Total Environ 439: 87–95

Berti WR, Jacobs LW (1998) Distribution of trace elements in soil from repeated sewage sludge applications. J Environ Qual 27:1280–1286

Bourdat-Deschamps M, Ferhi S, Bernet N, Feder F, Crouzet O, Patureau D, Montenach D, Moussard GD, Mercier V, Benoit P, Houot S (2017) Fate and impacts of pharmaceuticals and personal care products after repeated applications of organic waste products in long-term field experiments. Sci Total Environ 607:271–280

Businelli M, Gigliotti G, Giusquiani PL (1996) Trace element fate in soil profile and corn plant after massive applications of urban waste compost: a six-year study. Agrochimica 40(4):145–152

Cambier P, Michaud A, Paradelo R, Germain M, Mercier V, Guerin-Lebourg A, Revallier A, Houot S (2019) Trace metal availability in soil horizons amended with various urban waste composts during
17 years—monitoring and modelling. Sci Total Environ 651:2961–2974

Cambier P, Pot V, Mercier V, Michaud A, Benoit P, Revallier A, Houot S (2014) Impact of long-term organic residue recycling in agriculture on soil solution composition and trace metal leaching in soils. Sci Total Environ 499:560–573

Chabauty F, Pot V, Bourdat-Deschamps M, Bernet N, Labat C, Benoit P (2016) Transport of organic contaminants in subsoil horizons and effects of dissolved organic matter related to organic waste recycling practices. Environ Sci Pollut Res 23:6907–6918

Chaney RL (2012) Food safety issues for mineral and organic fertilizers. Adv Agron 117:51–116

Charlton A, Sakrabani R, Tyrrel S, Casado MR, McGrath SP, Crooks B, Cooper P, Campbell CD (2016) Long-term impact of sewage sludge application on soil microbial biomass: an evaluation using meta-analysis. Environ Pollut 219:1021–1035

Chaudri A, McGrath S, Gibbs P, Chambers B, Carlton-Smith C, Godley A, Bacon J, Campbell C, Aitken M (2007) Cadmium availability to wheat grain in soils treated with sewage sludge or metal salts. Chemosphere 66:1415–1423

Chaudri AM, Allain CMG, Barbosa-Jefferson VL, Nicholson FA, Chambers BJ, McGrath SP (2000) A study of the impacts of Zn and Cu on two rhizobial species in soils of a long-term field experiment. Plant Soil 221:167–179

Checkai R, Van Genderen E, Sousa JP, Stephenson G, Smolders E (2014) Deriving site-specific clean-up criteria to protect ecological receptors (plants and soil invertebrates) exposed to metal or metalloid soil contaminants via the direct contact exposure pathway. Integr Environ Assess Manag 10:346–357

Cobbett FD, Van Heyst BJ (2004) Measurements of GEM fluxes and atmospheric mercury concentrations (GEM, RGM and HgP) from an agricultural field amended with biosolids in Southern Ont., Canada (October 2004–November 2004). Atmos Environ 41: 2270–2282

Codex Alimentarius (2015) International standards for food. Food and Agriculture Organization of the United Nations (in French).

Delschen T (1999) Impacts of long-term application of organic fertilizers on soil quality parameters in reclaimed loess soils of the Rhineland lignite mining area. Plant Soil 213:43–54

Douglas JT, Aitken MN, Smith CA (2003) Effects of five non-agricultural organic wastes on soil composition, and on the yield and nitrogen recovery of Italian ryegrass. Soil Use Manage 19:135–138

ECHA (2008) Guidance on information requirements and chemical safety assessment. Chapter R.10: characterisation of dose [concentration]-response for environment. https://echa.europa.eu/documents/10162/13632/information requirements r10_en.pdf/bb902be7a503-4ab7-9036-d866b8ddce69. Accessed 18 July 2019

Ellert BH, Bettany JR (1995) Calculation of organic matter and nutrients stored in soils under contrasting management regimes. Can J Soil Sci 75:529–538

Feix I, Tremel A (2005) Soil contamination, soil-plant transfer. EDP Sciences Edition, France (in French)

Gigliotti G, Businelli D, Giusquiani PL (1996) Trace metals uptake and distribution in corn plants grown on a 6-year urban waste compost amended soil. Agr Ecosyst Environ 58:199–206

Giusquiani PL, Gigliotti G, Businelli D (1994) Long-term effects of heavy-metals from composted municipal waste on some enzyme activities in a cultivated soil. Biol Fert Soils 17:257–262

Hamner K, Kirchmann H (2015) Trace element concentrations in cereal grain of long-term field trials with organic fertilizer in Sweden. Nutr Cycl Agroecosys 103:347–358

Harmsen J (2007) Measuring bioavailability: from a scientific approach to standard methods. J Environ Qual 36:1420–1428

Heemsbergen DA, Warne MSJ, Broos K, Bell M, Nash D, McLaughlin M, Whatmuff M, Barry G, Pritchard D, Penney N (2009) Application of phytotoxicity data to a new Australian soil quality

guideline framework for biosolids. Sci Total Environ 407:2546–2556

Houot S, Clergeot D, Michelin J, Francou C, Bourgeois S, Caria G, Ciesielski H (2002) Agronomic value and environmental impacts of urban composts used in agriculture. In: Insam H, Riddech N, Klammer S (eds) Microbiology of composting. Springer, Berlin Heidelberg, pp 457–472

Houot S, Pons MN, Pradel M, Tibi A (2014) Organic waste materials valorization in soils used for agriculture and forestry, agronomic, environmental and socio-economic impacts. (in French)

Kabata-Pendians A, Pendias H (2001) Trace elements in soils and plants. CRC Press LLC, Third Edition

Kirchmann H, Borjesson G, Katterer T, Cohen Y (2017) From agricultural use of sewage sludge to nutrient extraction: a soil science outlook. Ambio 46:143–154

Lofts S, Chapman PM, Dwyer R, McLaughlin MJ, Schoeters I, Sheppard SC, Adams WJ, Alloway BJ, Antunes PMC, Campbell PGC et al (2007) Critical loads of metals and other trace elements to terrestrial environments. Environ Sci Technol 41(18):6326–6331

Lopez-Rayo S, Laursen KH, Lekfeldt JDS, Delle Grazie F, Magid J (2016) Long-term amendment of urban and animal wastes equivalent to more than 100 years of application had minimal effect on plant uptake of potentially toxic elements. Agr Ecosyst Environ 231: 44–53

McBride MB (1995) Toxic metal accumulation from agricultural use of sludge—are USEPA regulations protective? J Environ Qual 24:5–18

McBride MB, Spiers G (2001) Trace element content of selected fertilizers and dairy manures as determined by ICP-MS. Commun Soil Sci Plant Anal 32(1-2):139–156

McGrath SP, Chaudri AM, Giller KE (1995) Long-term effects of metals in sewage-sludge on soils, microorganisms and plants. J Ind Microbiol 14:94–104

McGrath SP, Lane PW (1989) An explanation for apparent losses of metals in a long-term field experiment with sewage sludge. Environ Pollut 60:235–256

Moolenaar SW, Beltrami P (1998) Heavy metal balances of an Italian soil as affected by sewage sludge and Bordeaux mixture applications. J Environ Qual 27:828–835

Moolenaar SW, Lexmond TM (1998) Heavy-metal balances of agro-ecosystems in the Netherlands. Neth J Agr Sci 46:171–192

Moolenaar SW, Lexmond TM, VanderZee SEATM (1997b) Calculating heavy metal accumulation in soil: a comparison of methods illustrated by a case-study on compost application. Agr Ecosyst Environ 66:71–82

Moolenaar SW, VanderZee SEATM, Lexmond TM (1997a) Indicators of the sustainability of heavy-metal management in agro-ecosystems. Sci Total Environ 201:155–169

NFU 44-051 (2006) French regulation established by AFNOR for organic amendments. (in French)

NFU 44-095 (2002) French regulation established by AFNOR for organic amendments—composts containing sewage sludges. (in French)

Noirot-Cosson PE, Vaudour E, Gilliot JM, Gabrielle B, Houot S (2016) Modelling the long-term effect of urban waste compost applications on carbon and nitrogen dynamics in temperate cropland. Soil Biol Biochem 94:138–153

Obriet F, Stauffer M, Goubard Y, Cheviron N, Peres G, Eden M, Revallier A, Vieuble-Gonod L, Houot S (2016) Multi-criteria indices to evaluate the effects of repeated organic amendment applications on soil and crop quality. Agr Ecosyst Environ 232:165–178

OECD (2016) Guidance on the incorporation of bioavailability concepts for assessing the chemical ecological risk and/or environmental threshold values of metals and inorganic metal compounds. OECD Environmental, Health and Safety Publications, Series on Testing & Assessment, No. 259
Oorts K (2018) Threshold calculator for metals in soil, a global tool for terrestrial risk assessment. Background document threshold calculator for metals in soil V2.0, Arche Consulting.

Oorts K, Schoeters I, Chowdhury J (2018) Threshold calculator for metals in soil: use of available data for derivation of metal soil quality standards for different scenarios and protection goals. 28th annual Meeting of the Society of Environmental Toxicology and Chemistry (SETAC) Europe.

Peltre C, Christensen BT, Dragon S, Icard C, Katterer T, Houot S (2012) RothC simulation of carbon accumulation in soil after repeated application of widely different organic amendments. Soil Biol Biochem 52:49–60

Riber L, Poulsen PHB, Al-Soud WA, Hansen LBS, Bergmark L, Brejnrod A, Norman A, Hansen LH, Magid J, Sorensen SJ (2014) Exploring the immediate and long-term impact on bacterial communities in soil amended with animal and urban organic waste fertilizers using pyrosequencing and screening for horizontal transfer of antibiotic resistance. Fems Microbiol Ecol 90:206–224

Richard G, Stengel P, Lemaire G, Cellier P, Valceshini E (2019) Agriculture and use of nature resources: feeding the world without degrading. QUAE Edition, France (in French)

Richards BK, Steenhuis TS, Peverly JH, McBride MB (1998) Metal mobility at an old, heavily loaded sludge application site. Environ Pollut 99:365–377

Sen Tran T, Côté D, N'Dayegamiye A (1996) Effects of long-term inputs of farmyard manures on the evolution of soils nutrients. Agrosol 4(1):20–30 (in French)

Sharma B, Sarkar A, Singh P, Singh RP (2017) Agricultural utilization of biosolids: a review on potential effects on soil and plant grown. Waste Manage 64:117–132

Six L, Smolders E (2014) Future trends in soil cadmium concentration under current cadmium fluxes to European agricultural soils. Sci Total Environ 485–486:319–328

Smith SR (2009) A critical review of the bioavailability and impacts of heavy metals in municipal solid waste composts compared to sewage sludge. Environ Int 35:142–156

Smolders E, Oorts K, Lombi E, Schoeters I, Ma YB, Zrna S, McLaughlin MJ (2012) The availability of copper in soils historically amended with sewage sludge, manure, and compost. J Environ Qual 41:1–9

Smolders E, Oorts K, Van Sprang P, Schoeters I, Janssen CR, McGrath SP, McLaughlin MJ (2009) Toxicity of trace metals in soil as affected by soil type and aging after contamination: using calibrated bio-availability models to set ecological soil standards. Environ Toxicol Chem 28(8):1633–1642

Sterckeman T, Gossiaux L, Guimont S, Sirguey C, Lin ZB (2018) Cadmium mass balance in French soils under annual crops: scenarios for the next century. Sci Total Environ 639:1440–1452

Toth G, Hermann T, Szatmari G, Pasztor L (2016) Maps of heavy metals in the soils of the European Union and proposed priority areas for detailed assessment. Sci Total Environ 565:1054–1062

Uprety D, Hejcmans M, Szakova J, Kunzova E, Tlustos P (2009) Concentration of trace elements in arable soil after long-term application of organic and inorganic fertilizers. Nutr Cycl Agroecosys 85: 241–252

Warman PR, Termeer WC (2005) Evaluation of sewage sludge, septic waste and sludge compost applications to corn and forage: Ca, Mg, S, Fe, Mn, Cu, Zn and B content of crops and soils. Bioresource Technol 96:1029–1038

Weber J, Karczewska A, Drozd J, Licznar M, Licznar S, Jamroz E, Kocowicz A (2007) Agricultural and ecological aspects of a sandy soil as affected by the application of municipal solid waste composts. Soil Biol Biochem 39:1294–1302

Yoshida H, Ten Hoeve M, Christensen TH, Bruun S, Jensen LS, Scheutz C (2018) Life cycle assessment of sewage sludge management options including long-term impacts after land application. J Clean Prod 174:538–547

Yu WT, Zhou H, Zhu XJ, Xu YG, Ma QA (2011) Field balances and recycling rates of micronutrients with various fertilization treatments in Northeast China. Nutr Cycl Agroecosys 90:75–86

Zhang FS, Li Y, Yang M, Li W (2012) Content of heavy metals in animal feeds and manures from farms of different scales in northeast China. Int J Env Res Pub He 9:2658–2668

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