Next Article in Journal
Reproduction Modes and Conservation Implications in Three Polyploid Sorbus Stenoendemics in Eastern Slovakia (Central Europe)
Next Article in Special Issue
Plant-Based Solutions for Non-Productive Sites Useful in the Management of Dry Land
Previous Article in Journal
The Arabidopsis LHT1 Amino Acid Transporter Contributes to Pseudomonas simiae-Mediated Plant Growth Promotion by Modulating Bacterial Metabolism in the Rhizosphere
Previous Article in Special Issue
Effect of Exogenous Plant Debris and Microbial Agents on Phytoremediation of Copper-Contaminated Soil in Shanghai
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Treatment of Landfill Leachate by Short-Rotation Willow Coppice Plantations in a Large-Scale Experiment in Eastern Canada

1
Institut de Recherche en Biologie Végétale, Université de Montréal, Montréal, QC H1X 2B2, Canada
2
Groupe Ramo, 457 Rang du Ruisseau des Anges Sud, Saint-Roch-de-l’Achigan, QC J0K 3H0, Canada
3
Department of Civil, Geological and Mining Engineering, Polytechnique Montréal, Montréal, QC H3C 3A7, Canada
*
Author to whom correspondence should be addressed.
Submission received: 9 December 2022 / Revised: 6 January 2023 / Accepted: 7 January 2023 / Published: 13 January 2023

Abstract

:
The treatment of leachate by vegetative filters composed of short-rotation willow coppice (SRWC) has been shown to be a cost-effective alternative to conventional and costly methods. However, few studies have considered the treatment capability of willow filters at a scale large enough to meet the industrial requirements of private landfill owners in North America. We report here on a field trial (0.5 ha) in which a willow plantation was irrigated with groundwater (D0) or aged leachate at two different loadings (D1 and D2, which was twice that of D1). Additionally, half of the D2-irrigated plots were amended with phosphorus (D2P). The system, which operated for 131 days, was highly efficient, causing the chemical oxygen demand concentration to drop significantly with the total removal of ammonia (seasonal average removal by a concentration of 99–100%). D2P efficacy was higher than that of D2, indicating that P increased the performance of the system. It also increased the willow biomass 2.5-fold compared to water irrigation. Leaf tissue analysis revealed significant differences in the concentrations of total nitrogen, boron, and zinc, according to the treatment applied, suggesting that the absorption capacity of willows was modified with leachate irrigation. These results indicate that the willow plantation can be effective for the treatment of landfill leachate in respect of environmental requirements.

1. Introduction

Waste management is a major issue for municipalities and governments around the world [1]. In Canada, more than 25 million tons of municipal solid wastes (MSWs) are generated each year and have to be eliminated. This represents approximately 729 kg cap−1 yr−1 [2]. One of the major challenges of landfills is the generation of polluted leachate resulting from the percolation of rainwater through the different layers of waste [3]. A multitude of contaminants, including ammonia, dissolved organic matter, phenolic compounds, chlorides, sulfates, trace elements (TE), and xenobiotic organic compounds, can be found in leachate, which can lead to the contamination of ground and surface water near the site [4,5,6,7]. Landfill leachate shows large variability in its pollutant concentration, depending in particular on the age of the waste stored in the landfill cells [8]. In the province of Quebec (Canada), leachate must be captured and treated before it is discharged into the surrounding environment [9]. Conventionally, leachates are treated by expensive physico-chemical, filtration, or biological processes. The current leachate treatment practice in Quebec uses biological processes operated in alternating aerobic or anaerobic modes for carbon, ammonia, and total Kjeldahl nitrogen (TKN) removal [10].
To overcome some of the treatment efficiency and cost challenges encountered in leachate treatment, an innovative approach was proposed and tested on operating landfill sites in Europe over the last twenty years. It consists of leachate treatment systems that combine aerated retention ponds and short-rotation coppice (SRC) using willows or poplars established on the roof of former landfill cells or on adjacent arable land irrigated with leachate [11,12,13,14]. This is a particularly appealing approach. Costs are low since the implementation and treatment both occur on-site [15]. Furthermore, the biomass produced in various sectors of chain values (such as the production of bioenergy or other bioproducts) that are associated with the ecological services offered by such structures (such as CO2 capture and storage, visual landscape improvement, and the increase in biodiversity) can be integrated every 3–5 years.
However, to secure these advantages, the phytofiltration systems must be functioning properly, which implies that (i) plant growth is not affected by the nature or quantities of the leachate applied and that (ii) contaminant leaching into the underground water after filtration is reduced or prevented (if on adjacent land). It has been reported that the efficiency of phytofiltration systems was subject to variations depending on plantation design setup and that symptoms of toxicity related to excessive salt concentrations (chloride ions) and nutrient imbalances (nitrogen and phosphorus) have been observed [11,16]. Other factors can also interfere with the effectiveness of leachate treatment, including soil texture, tree species or cultivars, the quality and quantity of the leachate applied, and irrigation methods, which can make it difficult to apply standard operating procedures in various environments [17,18].
Several studies showed that willows could effectively capture nitrogen, phosphorus, organic matter, and certain metals (Cd, Zn, Pb, Ni) in various applications, such as added organic fertilizers [19,20], municipal wastewater [21], or leachate [22]. More recently, Lachapelle-T et al. [23] carried out a pilot project in southern Quebec (Canada) using willows in SRC to treat municipal wastewater and concluded that the system allowed the removal of more than 95% of the nitrogen and phosphorus contained in the effluent. More importantly, the study also demonstrated that such plantations, characterized by very high evapotranspiration rates, could receive up to 25 mL ha−1 yr−1 of wastewater without negative impacts on willow growth or the environment [24]. Since landfill leachate has a high and variable contaminant concentration, little is known about the treatment ability of willows in SRC.
The objectives of this study were (1) to determine the efficacy of an SRC willow plantation when used as a vegetative filter in the pedoclimatic conditions of Quebec (Canada) for treating various loads of leachate from inoperative landfill cells and (2) to verify the impacts on soil, runoff water, and willow yields.

2. Results

2.1. Water Balance

2.1.1. Rainfall, Leachate and Water Irrigation

Rainfall rates measured in 2020 were below the climate normal 1980–2005 (Environment Canada, 2020), particularly from May to early July (Figure 1A). At the same time, the cumulative reference ETo has been estimated at 796 mm: a significant value that could be attributed mainly to the early heat waves observed in June and early July 2020 (Figure 1B). Data calculated from the evapotranspiration model showed D2, D2P > D1 > D0 rates reaching a peak in mid-August (10.7, 9.6, and 5.5 mm/day, respectively; Figure 1C), indicating that leachate irrigated willows evapotranspired about twice as much as water irrigated ones (D0). Leachate and groundwater irrigation was applied for 131 days in 2020 (from 4 June to 12 October). Under our experimental conditions, the average irrigation target rate of 100% (3 m3 day−1 for D0 and D1; 6 m3 day−1 for D2) was achieved on 49 days during the season. In total, the four irrigated experimental plots per treatment (804 m2) were fed with 237–244 m3 (D0–D1) or 439–453 m3 (D2–D2P), corresponding, respectively, to 2994 m3 ha−1 and 5561 m3 ha−1 (Figure 1D).

2.1.2. Root Zone Depletion Model

Using the RZDM results, the quantity of water remaining in the soil, combining irrigation, rainfall, and evapotranspiration, was determined for each treatment (Figure 2). Soils treated with leachate (D1, D2) had lower moisture levels when compared to those treated with water only (D0). In addition, high volumes of irrigation (D2), combined with high rainfall events that occurred in July and August, decreased water depletion to zero on four occasions during the summer.

2.2. Soil Analytical Results

2.2.1. Ammonia and Nitrate Content

Ammonia availability in the soil remained low during the growing season and did not vary significantly among D1, D2, and D2P treatments despite the high ammonia concentration in the leachate used for irrigation (Table 1). The average concentration of ammonium in the soil ranged from 7.5 to 9.5 mg kg−1 dry weight and was similar to the control treatment D0 (8.8 kg−1 dry weight) (Figure 3A), but no statistical difference was found between the treatments. The average nitrate concentration in the soil was low in June prior to irrigation (Figure 3B). After four months of irrigation (end of October), a significant increase in nitrate concentration was observed only in the soil treated with leachate. D0 was statistically different from D1, D2, and D2P. The relative stability in ammonia concentration throughout the growing season with the concomitant increase in nitrate was most likely due to aerobic bacterial nitrification and plant uptake.

2.2.2. Micro- and Macronutrient Soil Content Variation

Soil chemical analyses were conducted before (June) and after (October) the irrigation of willows to investigate the effect of each treatment on the variation of mineral concentrations in the soil at the end of the irrigation period. Although data presented in Table 1 revealed a relative variability in soil mineral concentration patterns, some tendencies emerged when each treatment was analyzed individually (June vs. October). The concentration of As, Cr, Mg, and Pb increased in all treatments, while B, Cd, and K concentrations increased only in D1, D2, and D2P treatments. The concentration of Ca decreased in all the treatments except D0. The concentration of the remaining minerals (Al, Ba, Co, Fe, Mn, Ni, P, Se, Zn) showed slight or no variation during the growing season. The total phosphorus concentration was higher in D2P than in D2 in June. A decrease in phosphorus was observed in October with D1 and D2P treatments but not with the D2 treatment. The application of a 2X leachate irrigation volume (D2/D2P) in the soil did not double the concentration of nutrients in the soil. K and B were the only nutrients for which a significant difference was observed between leachate-irrigated and water-irrigated treatments. In October, their concentration was around two to five times higher in leachate-treated willows than in water-treated ones.

2.3. Porewater Analysis

2.3.1. Chemical Oxygen Demand

The average seasonal influent concentration of chemical oxygen demand (COD) was 288 ± 43 mg L−1 in leachate and 9.8 ± 6.1 mg L−1 in groundwater. COD applied to D1; D2-D2P was 20 or 40 times more than in D0, respectively. The statistical analysis of porewater data revealed significant differences between D1, D2, and D2P treatments. The best removal performance was with D1 (80%), compared to D2 (42%) or D2P (62%), while technically, no removal was observed with DO treatment (Table 2). Interestingly, the percentage of removal for COD remained stable over the season (data not shown).

2.3.2. TKN and Ammonia

The average seasonal influential concentration of TKN was 0.22 mg/L in groundwater and 221 mg/L in leachate: very similar to that measured for NH4-N. This indicated that TKN was mainly added to plots by leachate and that the main nitrogen source present in leachate was NH4-N. Although leachate irrigated plots received between 967 (D1) and 1824 (D2 and D2P) times more TKN than groundwater irrigated ones, the efficiency of NH4-N removal was high for D1, D2, and D2P treatments: nearly 100%, approaching the NH4-N concentration measured in the porewater of D0 treatment (Table 2). For each treatment, the percentage of NH4-N removal remained stable over the growing season (data not shown).

2.3.3. Variation in Micro- and Macro-Element Porewater Content

A reduced number of minerals were analyzed in the porewater collected at two irrigation times (mid-season in August and end of the season in October). The total phosphorus concentration decreased throughout the season. However, the addition of phosphorus fertilizer in treatment D2P did not appear to significantly increase the porewater phosphorus content (Table 3). As a general rule, the concentration of Ca, Fe, Mg, Na, and K increased progressively in the porewater during the growing season for D2 and D2P treatments. Conversely, the concentration of these elements decreased in the D1 treatment except for Ca. Interestingly, the concentration of Fe was lower in leachate treatments than in the water treatment. Except for P, the data clearly distinguished two statistical groups (D1/D2/D2P and D0). Other elements (Al, Cu, Fe, Mn, Ni, and Zn) had no clear trend, and treatments could be discriminated against (data not shown).

2.4. Plant Analysis

2.4.1. Biomass Productivity

Very few willows died after the first growing season following their establishment in 2018. After two years of growth (2020), the average final yields of the aboveground biomass were 44.0, 32.8, 18.2, and 17.1 tons of dry mass per hectare (tons dry weight/ha) for D2P, D2, D1, and D0, respectively (Figure 4). Leachate irrigation had a significant effect on the biomass yield (p = 0.0136), which was, on average, 2.5-fold higher in D2P in comparison to D0.

2.4.2. Leaf Surface Area

The leaf surface area, from the largest to smallest area, ranked as follows: D2, D2P > D1, D0 (Figure 5A). Statistical analyses revealed that willows irrigated with the highest leachate volumes (D2, D2P) had larger leaves. However, phosphorus fertilization did not significantly affect the leaf surface area since similar yields were observed for D2 and D2P.

2.4.3. Leaf Pigment Content

Greater variability in the chlorophyll (a + b) and carotenoid concentrations were observed in the leaves of willows irrigated with water (Figure 5B,C). Pigment concentration increased during the summer, peaking in August for all treatments before declining notably in October. Peaks were significantly higher for leachate-irrigated D1, D2, and D2P treatments.

2.4.4. Leaf Nutrient Analysis

A chemical analysis of the leaves and stems harvested at the end of August 2020 was carried out in order to determine the concentration of a limited range of nutrients. A higher nitrogen percentage (between 56 and 86%) was found in the leaves from leachate-irrigated willows than in the leaves of plants irrigated with water (Table 4). A similar finding was observed for B, distinguishing two statistical groups, D0/D1 and D2/D2P, with an increase of about 40% for B in the second group. On the contrary, Zn was at a lower concentration in the leaves of leachate-irrigated willows (between 53% and 62% for the D1/D2/D2P statistical group). There was no clear statistical difference between treatments for the other elements.

3. Discussion

Many studies have demonstrated that willow vegetation filters can effectively treat various sources of wastewater on-site under a diversity of climates [25,26,27,28,29]. In this study, up to 7.5 mm d−1 (D2/D2P treatments) of leachate was applied to willows (Salix miyabeana ‘SX64’) grown on the clay layer topping an inoperative landfill cell. Over a cycle of 131 days, from June to October, plants were irrigated for about 110 days, thereby showing that willows are well adapted to water-saturated and compacted soils, with no negative impact on their growth and development [12,30,31]. Similar large-scale studies (thousands of m2) were carried-out with Salix miyabeana ‘SX67’ under Quebec’s humid continental climate to treat primary municipal wastewater [20,22]. Cultivar ‘SX67’, a cultivar of the same willow species utilized in this study, has been shown to have a very high-water demand, allowing a high irrigation rate of up to 2961 mm yr−1 with a former agricultural soil, which enables the treatment of up to 30,000 m3 ha−1 yr−1 of municipal wastewater [23,32]. However, it was shown that high hydraulic loading rates (HLR) were applied to willows at the end of the season, while evapotranspiration was low and growth started reducing, causing deep percolation and nitrogen leaching into the groundwater, indicating that HLR should be modulated according to its willow filtering capability [33].
Evapotranspiration is one of the most critical factors in large-scale industrial plant-based treatments whose main objective is to minimize the wastewater volumes to be treated conventionally [34]. Willows can be classified among the species displaying high evapotranspiration rates [35]. Under the climate of Eastern Canada, many willow cultivars remain photosynthetically active very late in the season, when the majority of deciduous species have lost their foliage [36]. This is another advantage of using willows in the treatment of wastewater or leachate. In the current study, the vegetation filter was able to treat, by evapotranspiration, up to 5560 m3 ha−1 of leachate, exhibiting moderate electroconductivity (3.7–4.2 mS cm−1). However, the RZDM model predicted that, over the growing season, an average of 4.5% of the total water (rainfall and irrigation) received by willows would potentially end up as runoff outside the operational filtration unit (data not shown). This water balance model also revealed root zone depletion during most of the irrigation campaign (between 10 and 70 mm for D1 and D2/D2P; also validated by in situ soil moisture tension results; data not shown), indicating that additional work would be needed to maximize the retention capability of the clay soil used in the cultivation of such landfill leachate filter systems. The relative similarity in the RZDM pattern between D1 and D2 treatments, in contrast to the D0 treatment, could be partially explained by the gain in the aboveground biomass production of leachate-irrigated willows, which promotes the evapotranspiration rate of plants as described above and, consequently, decreases soil moisture.
A significant increase in the biomass yield following leachate addition is reported here. The average biomass production in our experimental conditions after two years reached 18 tons ha−1 for intermediate (D1) loading but exceeded values of 44 tons ha−1 for the D2P treatment. In contrast, the biomass yield was only 17 tons ha−1 in the groundwater treatment (D0). However, these values remain less impressive than those obtained by Jerbi et al. [32]. In that study, carried out in a similar climatic context but on agricultural land, annual yields of 40 tons ha−1 y−1 were obtained when the willows were irrigated with wastewater from municipal primary effluents. It is likely that the soil conditions, as well as the various characteristics of the two effluents, explain this difference. Still, leachate had a positive effect on willow growth and induced no signs of toxicity. Apart from the water supply, fertilization is described as a driving factor that influences biomass yield and, consequently, evapotranspiration [29]. Jerbi et al. [32] have shown that high nitrogen inputs to the soil by wastewater irrigation resulted in a significant increase in stomatal conductance, leaf area, as well as chlorophyll and nitrogen leaf content, while the fine root:aboveground biomass ratio decreased [32,37]. Our data closely match these findings since leaf area and pigment content was significantly higher with the highest leachate supply (providing 1 478 kg ha−1 NH4-N cumulative loading over the growing period) compared to water treatment. Moreover, the leachate-irrigated willows maintained their green pigments and leaves for approximately two to three weeks longer than the control treatment at the end of the growing season, which may be part of the reason for their additional growth. It would also explain why the maximal evapotranspiration rate measured in August was two times higher in D2-treated plants (10.7 mm day−1) than in D0-treated ones (5.5 mm day−1). Finally, the high nitrogen supply to the soil during leachate irrigation significantly increased the growth of the longest willow stems by the end of the growing season, contributing to a higher aboveground biomass yield.
In the current study, different removal patterns were observed, with leachate contributing up to 1907 kg COD ha−1 yr−1 and 1478 kg NH4-N ha−1 yr−1. Removal varied, on average, from 42% (D2) or 62% (D2P) to 80% (D1) for COD and was close to 100% for NH4-N. Furthermore, the system’s filtration capacity remained effective throughout the irrigation period. The high efficiency of willow vegetation filters for removing pollutants from wastewater is well documented. A trial conducted in Sweden in 2010 showed a removal efficiency of about 80% for organic carbon (TOC) and total N by the cultivar ‘Tora’ (S. schwerinii x S. viminalis) irrigated with leachate whose total loadings applied were 378 kg ha−1 yr−1 and 2021 kg ha−1 yr−1 for the total N and TOC, respectively [22]. More recently, Jerbi et al. [32] reported high removal efficiencies (93.4% for COD, 99.3% for NH4-N, and 98.3% for P) with S. miyabeana ‘SX67’ irrigated with municipal wastewater contributing to 8700 kg COD ha−1 yr−1, 1245 kg N ha−1 yr−1 and 121 kg P ha−1 yr−1. While willows can take up NH4-N for nutritional purposes, a bacterial nitrification process could also explain its removal efficiency. The increase in NO3–N in the superficial soil horizon was between 3.75-fold for D1 and 6.0-fold for D2/D2P at the end of the growing season. It cannot be excluded that the nitrification process could be underestimated due to the uptake by willows of this preferred nutritional nitrogen form from the soil [38]. Interestingly, we have noticed that the NO3–N concentration decreased by about 2–3 fold in May 2021 in leachate irrigated soil (data not shown), suggesting that a nitrification process could occur during the winter (from November 2020 to May 2021) before initiating a new irrigation cycle. Jones et al. [18] pointed out that highly concentrated leachates (EC > 4.0 mS cm−1) could affect the filtering capacity of the plant system, thereby limiting the volumes that could be processed (<250 m3 ha−1 yr−1). This could be one of the explanations for the lower COD removal observed in the porewater of the D2 leachate treatment (EC 4.2 mS cm−1). However, the filtering capacity of our willow filtering system was not adversely impacted by this EC range for the removal of NH4-N, regardless of the leachate dose applied to the willows. This finding concurs with data obtained with the leachate-irrigated S. matsudana ‘Levante’, indicating that the NH4-N concentration significantly decreased under the lowest leachate supply (i.e., around EC 4.0 mS/cm; 29). Moreover, concentrations detected in the soil porewater after irrigation with moderately concentrated leachate (i.e., 111 mg COD L−1 and 0.82 mg NH4-N L−1 for D2P treatment) did not exceed the standard limits imposed by Quebec regulations, which are 25 mg NH3-N L−1 and 150 mg BOD5 L−1 for the release of treated leachate in the environment [39].
Phosphorus fertilization improved leachate treatment efficiency, increasing COD, TKN, and NH4-N removal by 44%, 32%, and 64%, respectively. Soil matrix physico-chemical characteristics of the site (i.e., high clay content, low organic matter percentage, low initial phosphorus soil concentration, and soil compaction) could negatively affect the phosphorus availability of the soil and, therefore, its assimilation by plants. Furthermore, leachate added very low amounts of phosphorus to the irrigated willow filters. Consequently, one treatment (D2P) was added to the experimental design, consisting of a phosphorus amendment of 40 kg P2O5 ha−1 yr−1 when applied to willows irrigated with the highest leachate dose (D2 treatment), which corresponds to willows requirements in SRC [40,41,42]. Phosphorus is essential to growth and root proliferation at the initial stages of establishment. When associated with nitrogen, it is highly effective on willow aboveground and with a root-dry biomass yield [43]. In the present study, higher aboveground biomass yield (+25%) and leachate treatment efficiency (+44% for COD, +32% for TKN, and +64% for NH4-N compared to D2 treatment) were achieved with P fertilization. These findings suggest that phosphorus could influence the soil bacterial community, especially species involved in degrading the organic matter and nitrogen present in the leachate. Indeed, previous research has reported phosphorus availability to be one of the major variables regulating the abundance and diversity of the bacterial community, and this shaping occurred only when phosphorus was associated with carbon or nitrogen [44,45,46].
In the present study, higher amounts of nutrients were measured in leaf tissues after leachate irrigation (especially B, Fe, and N). The significant strong increase in the N content in leaf tissues of leachate-irrigated willows (+56% in D1, +86% in D2, and +76% in D2P) was likely correlated with the high amount that nitrogen leachate contributed to the soil, thus suggesting the better availability and recovery of essential nutrients by plants for their growth. The high removal efficiency of nitrogen achieved by willows has also been documented in other studies [21,22,32,47]. Nutrient interactions in crop plants are multiple and complex, acting at several physiological levels, including transportation and signaling [48]. They may induce deficiencies or toxicities, modify growth responses, and/or modify nutrient composition [49]. For instance, a positive interaction between phosphorus and nitrogen was reported in the context of their mutual absorption by plants and synergetic effects on plant stimulation. On the other hand, applied high levels of phosphorus may induce an iron deficiency in plants [49]. In the context of the phytoremediation described here, a well-balanced supply of nutrients could be critical for optimizing both leachate treatment and the production of renewable high-quality aboveground wood biomass (willows are usually coppiced every three years) for commercial use.
Although the leachates used in this pilot study exhibited moderate sodium and chloride concentrations (around 300 mg L−1) and EC, we observed that the EC of porewater increased gradually over the course of the season with leachate treatments, while it remained low and stable with the water control. This increase was faster for D2/D2P than for D1 at the beginning of the irrigation campaign, reflecting the different volumes of leachate applied for irrigation. At the end of the irrigation period, it was 3.1, 4.2, and 3.7 mS cm−1 for D1, D2, and D2P, respectively, thus closely approaching the value measured in the leachate (data not shown). In addition, in the porewater of the D2 treatment, the total sodium reached 402 mg L−1 at the end of the irrigation campaign. As such, there is no evidence that this salinity rate was a factor limiting the growth of willows and their efficiency in treating leachate. Nevertheless, from a long-term perspective, the progressive salinization of the soil in which the willows grow could have an antagonist effect on phytoremediation performance. Although several willow cultivars have shown good tolerance to moderate salinity treatment (7 dS m−1), Canadian indigenous willows were reported to have a higher potential for long-term survival under severe salinity treatment (14 dS m−1) [50]. Moreover, recent works have reported an increased adaptation to soil salinity among tetraploid hybrids, thanks to a different Na+ balance between the leaves and roots, which is also characteristic of the cultivar ‘SX64’ used in this study, compared to diploid ones [51].

4. Materials and Methods

4.1. Site Description

The experiment was conducted at the engineered landfill site operated by waste management (WM) located in the municipality of Sainte-Sophie in Southern Quebec (45°47′02.98″ N, 73°54′12.38″ W) (Figure 6). In early July 2018, a 0.5 ha willow plantation was planted on top of a former (inoperative) landfill cell previously capped with a 1 m high clay layer and grassed in the early 1980s. One-year-old, 20 cm long, unrooted dormant cuttings of Salix miyabeana ‘SX64‘ were planted at a density of 16,000 plants/ha (in double rows spaced 1.5 m apart) in early July and coppiced in mid-November 2018. Willow crops were irrigated with leachate from the former landfill cells or groundwater for the two subsequent growing seasons (2019–2020). The climate of Sainte-Sophie is characterized by humid cold winters and hot, humid summers with a monthly average temperature of −12.7 ± 3.1 °C (January, coldest month) to 19.5 ± 0.9 °C (July, hottest month) over a 20-year period [52]. Weather conditions, including wind speed, humidity, rainfall, temperature, pressure, and solar radiation, were measured at 5 s intervals by a Vantage Pro2 Plus weather station, which was installed on site according to the manufacturer’s instructions (Davis Instruments Corp., 2004). Rain on site for the May-October growing seasons in 2019 (479 mm) and 2020 (526 mm) was below the climate normal for the same period (556 mm).

4.2. Experimental Setup and Willow Treatment

The 0.4 ha strip-block experimental setup was divided into four blocks (B1–B4) comprising a total of 20 equally sized 9.14 m × 22 m rectangular plots, with four treatments (D0, D1, D2, and D2P) replicated four times (Figure 1). Each plot was assigned one of four treatments corresponding to a maximum irrigation rate of 3 m3 day−1 of leachate for D1, 6 m3 day−1 of leachate for D2, 6 m3 day−1 for leachate, plus 40 kg P2O5 ha−1 yr−1 of phosphorus fertilizer (D2P) or 3 m3 day−1 of water for D0. Plot assignments for treatments D1 and D0 were selected based on the irrigation lines system design, while plot assignments between treatments D2 and D2P were randomized (Figure 6). The experimental area was buffered by outer rows of willows of the same width to prevent edge effects and was surrounded by an unplanted zone 8–10 m wide containing a polymer-lined soil berm to collect rainwater and snowmelt runoff. Leachate was pumped from former (inoperative) landfill cells and diverted into an RT1 storage reservoir (5800 L). The stored leachate was then pumped to the irrigation staging area, which was equipped with a decanting tank and a leachate storage tank (5800 L). Water was pumped from groundwater wells to another storage tank (5800 L) at the irrigation staging area (Figure 7). The irrigation system consisted of four distribution lines (one for D0 and D1; two for D2) which were installed with irrigation sprinklers every 7 m. An irrigation control and recording system comprising flowmeters, pressure gauges, and programmable valves were installed to modulate the volume of leachate/water to be applied to the experimental area according to rainfall levels and the evapotranspiration capacity of willows during the growing season. One suction-cup lysimeter (model 1900, Soil Moisture Corp.) was installed at a depth of 50 cm in the middle of each plot to collect soil porewater. Two irrigation campaigns were conducted during the growing season from 2019 to 2020 (88 and 131 days of irrigation, respectively), but only data for the second year are presented.

4.3. Soil, Leachate, Groundwater and Porewater Sampling and Analysis

An initial characterization of the soil surface (0–30 cm) was performed prior to irrigation in early July 2019, consisting of four samples per block, each of which comprised 12 homogenized subsamples: four per plot (Table 5). In 2020, the surface soil was sampled twice (before and after the growing season, i.e., early June and late October). Soil samples were collected using a soil core sampler, previously washed with alconox detergent, and rinsed with distilled water between each plot sampling. Pre-irrigation and post-irrigation sampling consisted of 16 analytical samples (one per plot), each comprising four homogenized subsamples collected from each plot (n = 4). For some analyses, sampling consisted of four analytical samples (one per treatment), each comprising three homogenized subsamples collected from each plot (n = 1). These samples were analyzed for pH, macronutrients (NH4-N, NO3-N, Ca, K, Mg, Mn, P), and trace elements (Al, As, B, Ba, Cd, Co, Cr, Cu, Fe, Pb, Se, Zn). Phosphorus, K, Ca, Mg, Al, Mn, Cu, Zn, B, and Fe were extracted by the Mehlich III method. A total acid extraction method was used for the other trace elements.
The leachate and groundwater used for irrigation and the soil porewater were sampled every 2–3 weeks from the start of irrigation to the end of the growing season for 10 campaigns. Leachate and groundwater were sampled directly from the storage tanks for total suspended solids (TSS), volatile suspended solids (VSS), biochemical oxygen demand (BOD5), chemical oxygen demand (COD), ammonia nitrogen (NH4-N), nitrate nitrogen (NO3-N), total Kjeldahl nitrogen (TKN), total and orthophosphates, Ca, Cl, Fe, K, Mg, Na, Zn, total phenolic compounds, pH, electroconductivity (EC), oxidation-reduction potential (ORP), and total dissolved solids (TDS). pH, ORP, EC, and TDS were measured in the field using a Hanna combo tester (model HA-98129 from ERE Inc. Saint-Leonard, QC, Canada). The rest of the parameters were analyzed in a certified laboratory. The main physicochemical characteristics are presented in Table 6.
Soil porewater was sampled from lysimeters, which were depressurized to −65 centibars (cb) between 24 and 48 h in advance of the sampling by extraction with a low-flow peristaltic pump. The samples were analyzed for COD, NH4-N, and NO3-N in eight campaigns and for alkalinity, both dissolved and in total metals (Ca, Fe, Mg, K, Na, and Zn).

4.4. Plant Biomass

Biomass productivity was measured at the end of the season by weighing three randomly selected fresh trees per plot and per treatment (48 trees in total, cut at 10 cm from the ground surface) which were weighed on-site to determine the wet biomass using a portable electronic scale (GRAM SAFIR X50, Barcelona, Spain). Dry biomass was calculated based on a dry/wet ratio determined by drying wood samples at 70 °C for at least 72 h.

4.5. Plant Sampling and Analysis

The longest stem length of six randomly preselected willows (A-F) per experimental plot and per treatment was measured twice in 2020 (June and October). Three of these preselected willows (A-C) were used for plant analysis. For chlorophyll content determination, ten leaves were collected from each plant three times (mid-July, mid-August, and early October), starting at the fourth leaf from the tip of the longest stem, which was pooled and maintained in a cooler and then stored in a freezer at −80 °C until analysis using the method of Hiscox and Israelstan [54]. Twenty-five leaves of each plant were harvested in mid-July, five of them for measuring the leaf area index (WinFOLIA software program) and the rest determining the N, P, K, Mg, S, Al, B, Cu, Zn, Fe, Mn, and Zn leaf content based on the Mehlich extraction method. Before nutrient determination, the leaves were pooled by treatment and by block and kept in a cooler.

4.6. Data Analysis

A hydrological balance model was developed for the willow plantation based on a daily root zone depletion model (RZDM) [55]. The model consisted of a unit volume of soil with the inputs and outputs calculated for each day the system was in operation, readily available for the water (RAW), and total available water (TAW). The capillary rise was assumed to be zero. Irrigated volumes of water or leachate were measured by the electronic flowmeters installed downstream of the pumps. Rainfall, solar radiation, temperature, and wind speed were measured by the on-site weather station. Runoff volumes were calculated based on the soil conservation service (SCS) method for clayey soils [56]. Deep percolation was the unknown variable calculated in the water balance. Initial depletion was assumed to be zero (saturated condition) based on snowmelt, cool temperatures, and rainfall in May.
The reference evapotranspiration (ETo) was estimated by the Penman-Monteith equation [57], and the evapotranspiration estimate was calculated by multiplying the reference evapotranspiration by the monthly average crop coefficient for willows. The unfertilized treatment (D0) used the unfertilized willow crop coefficient, while the fertilized treatments (D1, D2, D2P) used the fertilized crop coefficient from Guidi et al. [58].
The ability of the vegetative filter to degrade or remove pollutants (i.e., treatment efficacy by concentration) was determined by comparing their concentration in the last leachate (for D1, D2, and D2P) or groundwater (D0) to the ones in the soil porewater using a relative percent difference (% RPD), calculated after each sampling. The last leachate referred to the one applied to willows just before soil porewater sampling. The treatment efficacy for the whole season was estimated by averaging the % RPD of each treatment.

4.7. Statistical Analysis

Groundwater, soil, and other experimental data were analyzed by a two-way analysis of variance, with no interaction statistical model, by using JMP (version 14). The models had concentration as the dependent variable and treatment and block as the independent variables. Biomass productivity data were analyzed by a two-way analysis of variance, with an interaction statistical model, and also by using JMP 14. The model had biomass productivity as a dependent variable and treatment, block, and treatment*block as independent variables, of which block and treatment*block were considered to have a random effect. In all cases, Tukey’s honest significant difference (HSD) post hoc tests were performed to determine the significance of treatment effects at the traditional threshold of 5%.

5. Conclusions

Short-rotation willow coppice filters are a valuable alternative to other current approaches for treating inoperative landfill leachates under a cold humid climate. Under these experimental climatic conditions, willow cultivars (Salix ssp.) have a relatively long growing season, making the process potentially economically valuable at an industrial scale and, thereby, more attractive. This represents a less energy-consuming and more sustainable alternative for companies that manage landfill sites.
The maximum leachate loading for this pre-industrial pilot reached 5560 m3/ha in 2020, and we showed that willows were effective at decreasing the concentration of NH4-N in the effluent and in the soil. However, the efficacy of the vegetative filter for treating COD tended towards saturation when a 2x leachate loading was used for irrigation. The application of leachate on willow plantations modified the concentration of several nutrients present in the leaves, suggesting altered assimilation patterns. As the willows were established on the roof of a landfill cell, mainly composed of clay, some additional tests should be carried out to improve the texture and nutritional properties of this topsoil, increase the leachate treatment capability of this type of vegetative filtering process, and prevent any off-site leakage. Similarly, the presence of contaminants of emerging concern is not investigated in the current study, as well as the possible land quality degradation by salts and trace elements supplied; these are other factors to be considered in order to ensure the long-term survival of willow shrubs and the functionality of the filtering system when operating at an industrial scale.

Author Contributions

P.B.: conceptualization, formal analysis, funding acquisition, administration, writing the original draft; A.P.: methodology, formal analysis, investigation, data curation, analysis and writing; X.L.-T. and L.-C.B.: conceptualization, formal analysis, methodology, funding acquisition, review of the manuscript; Y.C.: conceptualization, validation, review of the manuscript, administration, funding acquisition; F.E.P.: conceptualization, formal analysis, methodology, funding acquisition, review of the manuscript; M.L.: conceptualization, formal analysis, methodology, funding acquisition, administration; writing and critical review of the final version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Funding was provided through a federal NSERC Collaborative Research and Development Grant (RDCPJ-525307-18) and a provincial CRIBIQ Grant (2017-079-C35).

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank Waste Management Québec Inc. for its financial and technical support of this project. Our thanks to Noémie Legault and her technical team for their helpful assistance in processing the samples. Thanks also to Karen Grislis for the critical review of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Das, S.; Lee, S.H.; Kumar, P.; Kim, K.H.; Lee, S.S.; Bhattacharya, S. Solid waste management: Scope and the challenge of sustainability. J. Clean. Prod. 2019, 228, 658–678. [Google Scholar] [CrossRef]
  2. Giroux, L. State of Waste Management in Canada (No. PN 1528); Conseil Canadien des Ministres de L’environnement: Winnipeg, MB, Canada, 2014. [Google Scholar]
  3. Lema, J.M.; Mendez, R.; Blazquez, R. Characteristics of landfill leachates and alternatives for their treatment: A review. Water Air Soil Pollut. 1988, 40, 223–250. [Google Scholar] [CrossRef]
  4. Christensen, T.H.; Baun, D.L. Speciation of heavy metals in landfill leachate: A Review. Waste Manag. Res. 2004, 22, 3–23. [Google Scholar]
  5. Ray, C.; Chan, P. Heavy metals in landfill leachate. Int. J. Environ. Stud. 2008, 27, 225–237. [Google Scholar] [CrossRef]
  6. Peng, Y. Perspectives on technology for landfill leachate treatment. Arab. J. Chem. 2017, 10, S2567–S2574. [Google Scholar] [CrossRef] [Green Version]
  7. Luo, H.; Zeng, Y.; Cheng, Y.; He, D.; Pan, X. Recent advances in municipal landfill leachate: A review focusing on its characteristics, treatment, and toxicity assessment. Sci. Total Environ. 2020, 703, 135468. [Google Scholar]
  8. Kjeldsen, P.; Barlaz, M.A.; Rooker, A.P.; Baun, A.; Ledin, A.; Christensen, T.H. Present and long-term composition of MSW landfill leachate: A review. Crit. Rev. Environ. Sci. Technol. 2002, 32, 297–336. [Google Scholar] [CrossRef]
  9. MELCC, Règlement sur L’Enfouissement ET L’Incinération des Matières Résiduelles. Available online: https://www.environnement.gouv.qc.ca/matieres/reimr.htm (accessed on 1 April 2022).
  10. Renou, S.; Givaudan, J.G.; Poulain, S.; Dirassouyan, F.; Moulin, P. Landfill leachate treatment: Review and opportunity. J. Haz. Mat. 2008, 150, 468–493. [Google Scholar] [CrossRef]
  11. Dimitriou, I.; Aronsson, P.; Weih, M. Stress tolerance of five willow clones after irrigation with different amounts of landfill leachate. Bioresour. Technol. 2006, 97, 150–157. [Google Scholar] [CrossRef] [PubMed]
  12. Bialowiec, A.; Wojnowska-Baryla, I.; Agopsowicz, M. The efficiency of evapotranspiration of landfill leachate in the soil-plant system with willow Salix amygdalina L. Ecol. Eng. 2007, 30, 356–361. [Google Scholar] [CrossRef]
  13. Zalesny, R.S.; Bauer, E.O. Evaluation of Populus and Salix continuously irrigated with landfill leachate II. Soils and early tree development. Int. J. Phytoremediation 2007, 9, 307–323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Rosenkranz, T. Phytoremediation of Landfill Leachate by Irrigation to Willow Short-Rotation Coppice. Master’s Thesis, Swedish University of Agricultural Sciences, Uppsala, Sweden, 2013. [Google Scholar]
  15. Livingstone, D.; Smyth, B.M.; McDonald, J.; Johnston, C. Willow for landfill leachate remediation and bioenergy, is it an effective solution? In Proceedings of the 15th Conference on Sustainable Development of Energy, Water and Environment Systems, Cologne, Germany, 1–5 September 2020. [Google Scholar]
  16. Stephens, W.; Tyrrel, S.F.; Tiberghien, J. Irrigating short rotation coppice with landfill leachate: Constraints to productivity due to chloride. Bioresour. Technol. 2000, 75, 227–229. [Google Scholar] [CrossRef]
  17. Dimitriou, I. Performance and Sustainability of Short-Rotation Energy Crops Treated with Municipal and Industrial Residues. Ph.D. Thesis, Acta Universitatis Agriculturae Sueciae, Swedish University of Agricultural Sciences, Uppsala, Sweden, 2005. [Google Scholar]
  18. Jones, D.L.; Williamson, K.L.; Owen, A.G. Phytoremediation of landfill leachate. Waste Manag. 2006, 26, 825–837. [Google Scholar] [CrossRef]
  19. Cavanagh, A.; Gasser, M.O.; Labrecque, M. Pig slurry as fertilizer on willow plantation. Biomass Bioenergy 2011, 35, 4165–4173. [Google Scholar] [CrossRef]
  20. Guidi, W.; Pitre, F.; Labrecque, M. Short-rotation coppice of willows for the production of biomass in Eastern Canada. In Biomass Now-Sustainable Growth and Use; Matovic, M.D., Ed.; Tech Open Science: London, UK, 2013; ISBN 978-953-51-1105-4. [Google Scholar]
  21. Nissim, W.G.; Jerbi, A.; Lafleur, B.; Fluet, R.; Labrecque, M. Willows for the treatment of municipal wastewater: Performance under different irrigation rates. Ecol. Eng. 2015, 81, 395–404. [Google Scholar] [CrossRef]
  22. Dimitriou, I.P.; Aronsson, P. Landfill leachate treatment with willows and poplars–Efficiency and plant response. Waste Manag. 2010, 30, 2137–2145. [Google Scholar] [CrossRef]
  23. Lachapelle-T, X.; Labrecque, M.; Comeau, Y. Treatment and valorization of a primary municipal wastewater by a short rotation willow coppice vegetation filter. Ecol. Eng. 2019, 130, 32–44. [Google Scholar] [CrossRef]
  24. Sas, E.; Hennequin, L.M.; Frémont, A.; Jerbi, A.; Legault, N.; Lamontagne, J.; Fagoaga, N.; Sarrazin, M.; Hallett, J.P.; Fennell, P.S.; et al. Biorefinery potential of sustainable municipal wastewater treatment using fast-growing willow. Sci. Total. Environ. 2021, 792, 148146. [Google Scholar] [CrossRef] [PubMed]
  25. Hänel, M.; Istenič, D.; Brix, H.; Arias, C.A. Wastewater-Fertigated Short-Rotation Coppice, a Combined Scheme of Wastewater Treatment and Biomass Production: A State-of-the-Art Review. Forests 2022, 13, 810. [Google Scholar] [CrossRef]
  26. McCracken, A.R.; Chris, R.; Johnston, C.J. Potential for wastewater management using energy crops. Manag. Econ. Eng. Agric. Rural. Dev. 2015, 15, 275–284. [Google Scholar]
  27. Aronsson, P.; Dahlin, T.; Dimitriou, I. Treatment of landfill leachate by irrigation of willow coppice–Plant response and treatment efficiency. Environ. Pollut. 2010, 158, 795–804. [Google Scholar] [CrossRef] [Green Version]
  28. Holm, B.; Heinsoo, K. Municipal wastewater application to Short Rotation Coppice of willows-Treatment efficiency and clone response in Estonian case study. Biomass Bioenerg 2013, 57, 126–135. [Google Scholar] [CrossRef]
  29. Guidi Nissim, W.; Palm, E.; Camilla Pandolfi, C.; Mancuso, S.; Azzarello, E. Willow and poplar for the phyto-treatment of landfill leachate in Mediterranean climate. J. Environ. Manag. 2021, 277, 111454. [Google Scholar] [CrossRef]
  30. Kuzovkina, Y.A.; Knee, M.; Quigley, M.F. Soil compaction and flooding effects on the growth of twelve Salix L. species. J. Environ. Hort. 2004, 22, 155–160. [Google Scholar]
  31. Amlin, N.A.; Rood, S.B. Inundation tolerances of riparian willows and cottonwoods. J. Am. Water Resour. Assoc. 2001, 37, 1709–1720. [Google Scholar] [CrossRef]
  32. Jerbi, A.; Brereton, N.J.B.; Sas, E.; Amiott, S.; Lachapelle-T, X.; Comeau, Y.; Pitre, F.E.; Labrecque, M. High biomass yield increases in a primary effluent wastewater phytofiltration are associated to altered leaf morphology and stomatal size in Salix miyabeana. Sci. Total Environ. 2020, 738, 139728. [Google Scholar] [CrossRef]
  33. Amiot, S.; Jerbi, A.; Lachapelle-T, X.; Frédette, C.; Labrecque, M.; Comeau, Y. Optimization of the wastewater treatment capacity of a short rotation willow coppice vegetation filter. Ecol. Eng. 2020, 158, 106013. [Google Scholar] [CrossRef]
  34. Frédette, C.; Labrecque, M.; Comeau, Y.; Brisson, J. Willows for environmental projects: A literature review of results on evapotranspiration rate and its driving factors across the genus Salix. J. Environ. Manag. 2018, 246, 526–537. [Google Scholar] [CrossRef]
  35. Curneen, S.J.; Gill, L.W. A comparison of the suitability of different willow varieties to treat on-site wastewater effluent in an Irish climate. J. Environ. Manag. 2014, 133, 153–161. [Google Scholar] [CrossRef] [PubMed]
  36. Palm. E.; Guidi Nissim; W.; Gagnon-Fee, D.; Labrecque, M. Photosynthetic patterns during autumn in three different Salix cultivars grown on a brownfield site. Photosynth. Res. 2022, 154, 155–167. [Google Scholar] [CrossRef] [PubMed]
  37. Jerbi, A.; Guidi Nissim, W.; Fluet, R.; Labrecque, M. Willow root development and morphology changes under different irrigation and fertilization regimes in a vegetation filter. Bioenergy Res. 2015, 8, 775–787. [Google Scholar] [CrossRef]
  38. Tischner, R. Nitrate uptake and reduction in higher and lower plants. Plant Cell Environ. 2000, 23, 1005–1024. [Google Scholar] [CrossRef]
  39. MELCC. Règlement sur L’Enfouissement et L’Incinération de Matières Résiduelles, Q2, r. 19. 2022. Available online: https://www.legisquebec.gouv.qc.ca/fr/document/rc/Q-2 (accessed on 1 April 2022).
  40. Caslin, B.; Finnan, J.; Johnston, C.; McCracken, A.; Walsh, L. Short Rotation Coppice Willow Best Practice Guidelines; AFBI, Agri-Food and Bioscience Institute: Newforge Lane, Belfast, Antrim, 2010. [Google Scholar]
  41. Ericsson, T. Nutrient cycling in energy forest plantations. Biomass Bioenerg 1994, 6, 115–121. [Google Scholar] [CrossRef]
  42. Adegbidi, H.; Volk, T.; White, E.; Abrahamson, L.; Briggs, R.; Bickelaupt, D. Biomass and nutrient removal by willow clones in experimental bioenergy plantation in New York State. Biomass Bioenerg 2001, 20, 399–411. [Google Scholar] [CrossRef]
  43. Martin, P.J.; Stephens, W. Willow growth in response to nutrients and moisture on a clay landfill cap soil. I. Growth and biomass production. Bioresour. Technol. 2006, 97, 437–448. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Cheng, H.; Yuan, M.; Duana, Q.; Suna, R.; Shena, Y.; Yub, Q.; Li, S. Influence of phosphorus fertilization patterns on the bacterial community in upland farmland. Ind. Crops Prod. 2020, 155, 112761. [Google Scholar] [CrossRef]
  45. Siegenthaler, M.B.; Ramoneda, J.; Frossard, E.; Meszaros, E. Microbial community responses to phosphorus and nitrogen inputs in the organic soil horizons of two contrasting temperate beech forests. Appl. Soil Ecol. 2022, 172, 104357. [Google Scholar] [CrossRef]
  46. Heuck, C.; Weig, A.; Spohn, M. Soil microbial biomass C:N:P stoichiometry and microbial use of organic phosphorus. Soil Biol. Biochem. 2015, 85, 119–129. [Google Scholar] [CrossRef]
  47. Godley, A.; Alker, G.; Hallett, J.; Marshall, R.; Riddell-Black, D. Landfill leachate nutrient recovery by Willow short rotation coppice. I. Yield, Tissue composition and wood quality. Arboric. J. 2004, 27, 281–295. [Google Scholar] [CrossRef]
  48. Bouain, N.; Shahzad, Z.; Rouached, A.; Abbas Khan, G.; Berthomieu, P.; Abdelly, C.; Poirier, Y.; Rouached, H. Phosphate and zinc transport and signalling in plants: Toward a better understanding of their homeostasis interaction. J. Exp. Bot. 2014, 65, 5725–5741. [Google Scholar]
  49. Fageria, V.D. Nutrient interactions in crop plants. J. Plant Nutr. 2001, 24, 1269–1290. [Google Scholar] [CrossRef]
  50. Huang, X.; Soolanayakanahally, R.J.; Guy, R.D.; Shunmugam, A.S.K.; Mansel, S.D. Differences in growth and physiological and metabolic responses among Canadian native and hybrid willows (Salix spp.) under salinity stress. Tree Physiol. 2020, 40, 652–666. [Google Scholar] [CrossRef] [PubMed]
  51. Cseri, A.; Borbely, P.; Poor, P.; Feher, A.; Sass, L.; Jancso, M.; Penczi, A.; Radi, F.; Gyuricza, C.; Digruber, T.; et al. Increased adaptation of an energy willow cultivar to soil salinity by duplication of its genome size. Biomass Bioenerg 2020, 140, 105655. [Google Scholar] [CrossRef]
  52. Environment Canada. Canadian Climate Normals. 2020. Available online: https://climate.weather.gc.ca/climate_normals/index_e.html (accessed on 1 August 2022).
  53. Saxton, K.E.; Rawls, W.J. Soil water characteristic estimates by texture and organic matter for hydrologic solutions. Soil Sci. Soc. Am. J. 2006, 70, 1569–1578. [Google Scholar] [CrossRef] [Green Version]
  54. Hiscox, J.D.; Israelstam, G.F. A method for the extraction of chlorophyll from leaf tissue without maceration. Can. J. Bot. 1979, 57, 1332–1334. [Google Scholar] [CrossRef]
  55. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements; Food and Agriculture Organization of the United Nations, Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998; ISBN 92-5-104219-5. [Google Scholar]
  56. Mishra, S.K.; Singh, V.P. Soil Conservation Service Curve Number (SCS-CN) Methodology, 2003rd ed.; Water Science & Technology Library Book 42; Springer Science & Business Media: Berlin, Germany, 2003. [Google Scholar]
  57. Zotarelli, L.; Dukes, M.D.; Romero, C.C.; Migliaccio, K.W.; Morgan, K.T. Step by Step Calculation of the Penman-Monteith Evapotranspiration (FAO-56 Method); Doc. AE459, UF/IFAS extension; University of Florida: Gainesville, FL, USA, 2010; pp. 1–10. [Google Scholar]
  58. Guidi, W.; Piccioni, E.; Bonari, E. Evapotranspiration and crop coefficient of poplar and willow short-rotation coppice used as vegetation filter. Bioresour. Technol. 2008, 99, 4832–4840. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Irrigation parameters in 2020. (A) Daily rainfall quantities, cumulative rainfall, and 1980–2005 climate normal; (B) A 10-day moving average reference evapotranspiration (ETo) and cumulative ETo; (C) Modeled 10-day moving average of evapotranspiration; (D) Cumulative daily irrigation volumes (m3) applied by treatment (corresponding to an irrigation area of 804 m2).
Figure 1. Irrigation parameters in 2020. (A) Daily rainfall quantities, cumulative rainfall, and 1980–2005 climate normal; (B) A 10-day moving average reference evapotranspiration (ETo) and cumulative ETo; (C) Modeled 10-day moving average of evapotranspiration; (D) Cumulative daily irrigation volumes (m3) applied by treatment (corresponding to an irrigation area of 804 m2).
Plants 12 00372 g001
Figure 2. Root zone depletion model applied to irrigated treatments with water (D0) or leachate (D1, D2) in 2020.
Figure 2. Root zone depletion model applied to irrigated treatments with water (D0) or leachate (D1, D2) in 2020.
Plants 12 00372 g002
Figure 3. (A,B): Ammonia (NH4-H) and nitrate (NO3-N) concentration (expressed in mg N kg−1 of dry weight, DW) in the soil before (June) and after irrigation (October). Different letters above groupings indicate significant differences at p < 0.05.
Figure 3. (A,B): Ammonia (NH4-H) and nitrate (NO3-N) concentration (expressed in mg N kg−1 of dry weight, DW) in the soil before (June) and after irrigation (October). Different letters above groupings indicate significant differences at p < 0.05.
Plants 12 00372 g003
Figure 4. Comparison of the yield measured at the end of the growing season in 2020. Different letters above box plots indicate significant differences at p < 0.05. The horizontal line inside each box-whisker represents the median, vertical small rectangles represent interquartile ranges, and circles are data-outliers.
Figure 4. Comparison of the yield measured at the end of the growing season in 2020. Different letters above box plots indicate significant differences at p < 0.05. The horizontal line inside each box-whisker represents the median, vertical small rectangles represent interquartile ranges, and circles are data-outliers.
Plants 12 00372 g004
Figure 5. Leaf parameters during the growing season in 2020. (A) Average leaf surface area expressed in cm2. (B,C) Average leaf chlorophyll and carotenoids concentrations (expressed in mg g−1 dry weight). Lowercase letters indicate statistical groupings (p < 0.05). The horizontal line inside each box-whisker represents the median, vertical small rectangles represent interquartile ranges, and circles are data-outliers.
Figure 5. Leaf parameters during the growing season in 2020. (A) Average leaf surface area expressed in cm2. (B,C) Average leaf chlorophyll and carotenoids concentrations (expressed in mg g−1 dry weight). Lowercase letters indicate statistical groupings (p < 0.05). The horizontal line inside each box-whisker represents the median, vertical small rectangles represent interquartile ranges, and circles are data-outliers.
Plants 12 00372 g005
Figure 6. Daily temperature 10-day average for the 2020 season and climate norm 1981–2010, and schematic diagram of the experimental setup.
Figure 6. Daily temperature 10-day average for the 2020 season and climate norm 1981–2010, and schematic diagram of the experimental setup.
Plants 12 00372 g006
Figure 7. Diagram representing the irrigation system.
Figure 7. Diagram representing the irrigation system.
Plants 12 00372 g007
Table 1. Variation in the seasonal concentration of trace elements and macronutrients in the soil before (June) and after (October) irrigation with groundwater (D0) and leachate (D1, D2, and D2P) in 2020.
Table 1. Variation in the seasonal concentration of trace elements and macronutrients in the soil before (June) and after (October) irrigation with groundwater (D0) and leachate (D1, D2, and D2P) in 2020.
Soil Reference D0D1D2D2P
ParametersUnitsJuly 2019nJuneOctoberJuneOctoberJuneOctoberJuneOctober
Almg/kg1181 ± 13341255 ± 2101273 ± 1761094 ± 1411024 ± 1601239 ± 1091298 ± 2071304 ± 2521228 ± 235
Asmg/kg<5.0 10.712.210.871.430.991.561.011.94
Bmg/kg0.5 ± 0.140.6 ± 0.20.4 ± 0.1 a1.5 ± 0.82.0 ± 0.9 b 1.3 ± 0.22.9 ± 0.5 b1.6 ± 0.4 2.3 ± 0.5 b
Bamg/kg198 ± 241140.2160162155.4132.2133.6147.8146.4
Cakg/ha5793 ± 72945064 ± 10135197 ± 9634508 ± 4694484 ± 3374393 ± 13834280 ± 14884978 ± 8794700 ± 444
Cdmg/kg< 0.5010.00250.1550.570.030.00250.090.572.27
Comg/kg16 ± 2110.812.51413.510.211.19.410.7
Crmg/kg81 ± 10 156726571.441.561.349.962.4
Cumg/kg10.3 ± 2.6410.8 ± 3.210.8 ± 2.38.5 ± 1.89.1 ± 1.28.7 ± 3.78.5 ± 4.710.4 ± 3.310.4 ± 2.0
Femg/kg345 ± 174340 ± 28358 ± 30351 ± 30372 ± 31367 ± 63325 ± 53328 ± 29371 ± 26
Kkg/ha531 ± 1044483 ± 88487 ± 122 a590 ± 217898 ± 96 b522 ± 197976 ± 308 b640 ± 1331002 ± 168 b
Mgkg/ha1292 ± 19441041 ± 1371139 ± 1881155 ± 2171362 ± 861050 ± 312 1083 ± 3041097 ± 1041248 ± 115
Mnmg/kg22.1 ± 6.2420.0 ± 11.717.0 ± 2.024.8 ± 3.220.5 ± 4.120.9 ± 11.716.0 ± 7.219.2 ± 5.3 19.5 ± 4.9
Nimg/kg47 ± 6135.940.243.140.4383426.834.2
Pkg/ha312 ± 2084441 ± 239412 ± 183141 ± 8786 ± 71326 ± 170341 ± 238418 ± 291330 ± 214
Pbmg/kg7.0 ± 0.510.48.430.056.580.0257.440.929.2
Note: underlined numbers are significantly different. Different letters in rows indicate significant differences (p < 0.05); ± indicates standard deviation.
Table 2. Efficacy of COD and nitrogen removal by the willow filter. Efficacy was calculated as the seasonal concentration average % RPD between influent (groundwater or leachate) and effluent (porewater).
Table 2. Efficacy of COD and nitrogen removal by the willow filter. Efficacy was calculated as the seasonal concentration average % RPD between influent (groundwater or leachate) and effluent (porewater).
ParametersTreatmentD0D1D2D2P
CODInfluent (mg L−1)9.8 ± 6.1288 ± 43288 ± 43288 ± 43
Effluent (mg L−1)28.0 ± 18 a58 ± 20 b198 ± 69 d111 ± 42 c
Efficiency %−185.079.831.261.6
NTKInfluent (mg L−1)0.22 ± 0.06221 ± 37221 ± 37221 ± 37
Effluent (mg L−1)0.90 ± 0.68 a2.3 ± 1.6 b5.6 ± 4.0 c3.8 ± 2.0 c
Efficiency %−309.199.997.598.3
NH4-NInfluent (mg L−1)0.23 ± 0.05218 ± 35218 ± 35218 ± 35
Effluent (mg L−1)0.06 ± 0.04 a0.07 ± 0.08 a,b2.29 ± 3.53 b,c0.82 ± 1.30 c
Efficiency %741009999.6
Note: Different letters in rows indicate significant differences (p < 0.05); ± indicates standard deviation.
Table 3. Variation in the seasonal concentration of macronutrients in the porewater of plots irrigated with groundwater (D0) and leachate (D1, D2, and D2P) in 2020.
Table 3. Variation in the seasonal concentration of macronutrients in the porewater of plots irrigated with groundwater (D0) and leachate (D1, D2, and D2P) in 2020.
ParameterUnits InfluentEffluent
WaterLeachateD0D1D2D2P
AugustOctoberAugustOctoberAugustOctoberAugustOctober
Camg/kg78.7232139.3 ± 28151 ± 39 a344.0 ± 77368 ± 101 b396 ± 80386 ± 108 b433.0 ± 111456 ± 116 b
Femg/kg1.4534.53.8 ± 4.32.7 ± 1.8 a0.09 ± 0.070.06 ± 0.01 b0.10 ± 0.090.15 ± 0.14 b0.04 ± 0.020.04 ± 0.02 b
Kmg/kg4.5267.78.4 ± 5.710 ± 1.6 a38 ± 2235 ± 11 a87 ± 42119 ± 89 b49 ± 1955 ± 36 a
Mgmg/kg39.2122.462 ± 1963 ± 22 a194 ± 63190 ± 66 b210 ± 15201 ± 27 b212 ± 82220 ± 73 b
Namg/kg21.733312 ± 311 ± 5 a262 ± 180250 ± 99 b396 ± 129467 ± 172 b371 ± 209377 ± 210 b
Pmg/kg0.00080.0060.17 ± 0.110.09 ± 0.020.50 ± 0.550.07 ± 0.050.14 ± 0.050.07 ± 0.040.24 ± 0.150.13 ± 0.09
Note: underlined numbers are significantly different. Different letters in rows indicate significant differences (p < 0.05); ± indicates standard deviation.
Table 4. Leaf micro- and macro-element concentrations in 2020 expressed as % of dry weight or in mg/kg dry weight with comparison to D0.
Table 4. Leaf micro- and macro-element concentrations in 2020 expressed as % of dry weight or in mg/kg dry weight with comparison to D0.
ParametersUnitsTreatment
DOD1% D0D2% D0D2P% D0
N%2.08 a3.25 b563.87 b863.66 b76
P%0.260.23−120.3150.2912
K%1.531.5941.647.21.669
Ca%1.431.32−81.25−131.23−14
Mg%0.290.33140.32100.290
S%0.360.43190.3980.360
Almg/kg22.5272022.3−0.917.8−21
Bmg/kg28.8 a31 a840.8 b4240 b39
Cumg/kg6.810.3511162818
Femg/kg59.576.52985.34368.315
Mnmg/kg45.348.36.753.51842−7.3
Znmg/kg137.5 a60.8 b−5652.8 b−6262.8 b−54
Note: underlined numbers are significantly different. Different letters in rows indicate significant differences (p < 0.05).
Table 5. Soil physical and chemical properties at a depth 0–30 cm (n = 16, July 2019).
Table 5. Soil physical and chemical properties at a depth 0–30 cm (n = 16, July 2019).
ParametersUnits Values ParametersUnitsValues
Density g/cm30.87 ± 0.02Crmg/kg78 ± 8
Porosity% vol 66.0 ± 0.7Cumg/kg44 ± 6
Sand wt%34.6 ± 5.9Femg/kg345 ± 4
Silt wt%13.3 ± 3.0Hgmg/kg0.091 ± 0.005
Claywt%52.1 ± 3.0Kkg/ha524 ± 43
Texture-ClayMgkg/ha1278 ± 188
Field capacity (FC) 1v%40 ± 3Mnmg/kg22 ± 3
Wilting point (WP) 1v%29 ± 2Namg/kg895 ± 113
Total available water (TAW) 1v% 12 ± 0.8Nimg/kg45 ± 5
TAW 1mm69.0 ± 3.5Pbmg/kg7.0 ± 0.5
Organic matter%4 ± 2Pkg/ha312 ± 219
pH 7.8 ± 0.2Smg/kg406 ± 152
Saturated hydraulic conductancecm/h0.12 ± 0.00Semg/kg<1.0
Agmg/kg<2.0 SO4%0.008 ± 0.001
Almg/kg1176 ± 129Znmg/kg85 ± 11
Asmg/kg<5.0 Chloride mg/kg20 ± 5
Bmg/kg0.5 ± 0.1Sulphidemg/kg327 ± 159
Bamg/kg193 ± 24Organic carbon%g/g1.5 ± 0.4
Cakg/ha5739 ± 622Total carbon%g/g1.7 ± 0.6
Cdmg/kg<0.50 NH4-N mg N/kg 8.0 ± 2.0
Comg/kg16 ± 2NO3-N mg N/kg 1.8 ± 0.6
Note: 1: FC, WP, and TAW calculated according to Saxton and Rawls [53]; ± indicates standard deviation.
Table 6. Physicochemical properties of aged leachate and groundwater used for the irrigation of willows in 2019 and 2020.
Table 6. Physicochemical properties of aged leachate and groundwater used for the irrigation of willows in 2019 and 2020.
LEACHATEGROUNDWATER
ParametersUnitsnInitial Characterization 2019nAverage Aged Leachate Values 2019nAverage Aged Leachate Values 2020nAverage Groundwater Values 2019nAVERAGE Groundwater Values 2020
TSS mg/L357 ± 39653 ± 37879 ± 1354.0 ± 2.481.8 ± 0.6
VSS mg/L323 ± 17922 ± 12831.0 ± 8.056.0 ± 2.28<5
BOD5mg/L317. ± 8.2631 ± 2979.4 ± 3.055.2 ± 6.070.68 ± 0.67
CODmg/L3204 ± 256302 ± 179288 ± 43619 ± 2999.8 ± 6.1
NH3-Nmg N/L 3177 ± 326260 ± 279213 ± 3560.17 ± 0.0890.22 ± 0.06
NO3-Nmg N/L 30.09 ± 0.0310.05--10.02--
TKNmg N/L 3173 ± 461240118011.111.1
o-PO4mg P/L 30.270 ± 0.01210.18530.047 ± 0.07310.005530.0050 ± 0.0003
P totalmg P/L 30.03 ± 0.0110.02030.103 ± 0.05810.0230.07 ± 0.04
Ca totalmg/L3202 ± 193178 ± 232232 ± 28372 ± 8279 ± 6
Clmg/L3198 ± 1513372297 ± 85136240.0 ± 0.7
SO4mg/L320 ± 2145--136--
Al totalmg/L40.033 ± 0.02310.211--10.003--
As totalmg/L40.0012 ± 0.0005--------
Ba totalmg/L40.370 ± 0.095--------
B totalmg/L1-11.53--10.005--
Cd totalmg/L20.00012--------
Cu totalmg/L40.001 ± 0.00110.000025--10.000025--
Fe total mg/L422 ± 14424 ± 19234.5 ± 5.030.63 ± 0.2120.95 ± 0.50
Hg totalmg/L4<0.00001--------
K total mg/L3147 ± 303220 ± 362268.0 ± 8.234.8 ± 0.524.5 ± 0.9
Mg totalmg/L--399 ± 102122 ± 16338.0 ± 5.5222 ± 6
Mn totalmg/L--10.78--10.10--
Na totalmg/L4170 ± 153264 ± 322333.0 ± 2.8322 ± 2222 ± 6
Ni totalmg/L40.0075 ± 0.003010.013--10.000025--
Pb totalmg/L4<0.003--------
Se totalmg/L40.003 ± 0.002 --------
Zn total mg/L30.34 ± 0.6040.0015 ± 0.001720.006 ± 0.00730.004 ± 0.00620.0008 ± 0.0004
Phenolsmg/L30.009 ± 0.00520.009 ± 0.00420.033 ± 0.01 020.0012<0.003
pH ---56.8 ± 0.1106.78 ± 0.2157.3 ± 0.4107.0 ± 0.2
ECmS/cm--54.2 ± 0.4103.76 ± 0.4150.66 ± 0.10100.62 ± 0.07
ORPmV--4−12.33 ± 781039.4 ± 28.14−29 ± 4710163 ± 89
TDSmg/L--42077 ± 21681893 ± 1994360 ± 9810307 ± 36
Note: BOD5 (biochemical oxygen demand); COD (chemical oxygen demand); EC (electroconductivity); ORP (oxidation-reduction potential); TDS (total dissolved solids); TSS (total suspended solids; VSS (volatile suspended solids); ± indicates standard deviation.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Benoist, P.; Parrott, A.; Lachapelle-T., X.; Barbeau, L.-C.; Comeau, Y.; Pitre, F.E.; Labrecque, M. Treatment of Landfill Leachate by Short-Rotation Willow Coppice Plantations in a Large-Scale Experiment in Eastern Canada. Plants 2023, 12, 372. https://0-doi-org.brum.beds.ac.uk/10.3390/plants12020372

AMA Style

Benoist P, Parrott A, Lachapelle-T. X, Barbeau L-C, Comeau Y, Pitre FE, Labrecque M. Treatment of Landfill Leachate by Short-Rotation Willow Coppice Plantations in a Large-Scale Experiment in Eastern Canada. Plants. 2023; 12(2):372. https://0-doi-org.brum.beds.ac.uk/10.3390/plants12020372

Chicago/Turabian Style

Benoist, Patrick, Adam Parrott, Xavier Lachapelle-T., Louis-Clément Barbeau, Yves Comeau, Frédéric E. Pitre, and Michel Labrecque. 2023. "Treatment of Landfill Leachate by Short-Rotation Willow Coppice Plantations in a Large-Scale Experiment in Eastern Canada" Plants 12, no. 2: 372. https://0-doi-org.brum.beds.ac.uk/10.3390/plants12020372

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop