Assessing Stormwater Nutrient and Heavy Metal Plant Uptake in an Experimental Bioretention Pond

With the purpose to study a solution based on Sustainable Urban Drainage Systems (SUDS) to reduce and treat stormwater runoff in urban areas, a bioretention pond (BP) was realized in the Agripolis campus of the University of Padova, Italy. The BP collected overflow water volumes of the rainwater drainage system of a 2270 m2 drainage area consisting almost entirely of impervious surfaces. Sixty-six Tech-IA® floating elements, supporting four plants each, were laid on the water surface. Eleven species of herbaceous perennial helophyte plants, with ornamental features, were used and tested. The early growth results of the BP functioning showed that nearly 50% of the total inflow water volume was stored or evapotranspirated, reducing the peak discharge on the urban drainage system. Among plants, Alisma parviflora, Caltha palustris, Iris ‘Black Gamecock’, Lysimachia punctata ‘Alexander’, Oenanthe javanica ‘Flamingo’, Mentha aquatica, Phalaris arundinacea ‘Picta’, and Typha laxmannii had the best survival and growth performances. A. parviflora and M. aquatica appeared interesting also for pollutant reduction in runoff water.


Introduction
The high rate of urbanization has resulted in a large increase of impervious coverage in the landscape which can reach a very high percentage of the urban surface. Impervious surfaces decrease rainfall infiltration into the soil increasing runoff in terms of both peak flow and volume [1,2]. Rainwater in the urban landscape is therefore mainly directed into the municipal drainage system, creating serious problems in case of heavy rains, such as local floods, river inundations, etc., and reducing water availability and quality [3]. Urban runoff can be and often is a significant source of water pollution, causing a decline of fisheries, swimming areas, and other beneficial attributes of water resources [4]. At the same time, climate changes are causing the intensification and concentration of rainfall events, exacerbating the problem [5,6].
To reduce the problem, some environmentally sustainable approaches to urban development have been proposed as an alternative to the traditional ones to better manage the runoff in urban areas [7][8][9][10]. A stormwater best management practice (BMP) is a technique, measure, or structural control that is used to manage the quantity and improve the quality of stormwater runoff in the most cost-effective manner [11,12]. BMPs include stormwater planting, open channels, porous pavements, etc., in addition to a set of overall site design strategies and highly localized, small-scale, decentralized source control techniques, also known as Low-Impact Development (LID) systems in the USA [13,14] or Water-Sensitive Urban Design (WSUD) in Australia [15,16]. To describe stormwater technologies,

Materials and Methods
The bioretention pond was set in the Agripolis Campus of the University of Padova (Figure 1), in Legnaro (45 • 35 N; 11 • 96 E). The area has an annual average temperature of 12.3 • C and an average minimum and maximum temperatures of −5.5 and 32.8 • C. The average annual rainfall is 811 mm, mostly distributed during the growing season, from April to November.
During Spring 2011, in proximity of a building and a parking lot of the campus, a soil area of about 70 m 2 was dug up to a depth of about 165 cm. With the excavation, a storage basin of about 44.5 m 3 (8.10 × 5.23 m wide and 1.05 m deep) was obtained, and the entire basin was lined with a 1.5 mm thick polyolefin film.
The BP was designed to collect the overflow volumes of the existing rainwater drainage system of an area of 2270 m 2 consisting almost entirely of impervious surfaces (an asphalted road, sidewalks, and a building roof) ( Figure 2). In fact, this drainage system has a good capacity to store and slowly let infiltrate into the soil all the runoff volumes from this area except during heavy or frequent rainfall events. The overflow volumes, which were previously discharged in a channel of the urban rainwater drainage system, were collected in a concrete sump and pumped into the BP. Water exceeding the storage capacity of the pond flowed into another sump from where it was pumped out in the channel of the local urban surface drainage system.
On July 2011, when the pond was almost full of water (102.5 mm in depth), 66 self-floating elements (Tech-IA ® , Padova, Italy), supporting four plants each, were laid down on the water surface. Tech-IA ® is a rectangular panel (0.50 × 0.90 m) produced in ethylene-vinil acetate (EVA), with eight gridded windows in which plants can be anchored. Its mass is 1732 g, and it may support a load  . The rainwater drainage system with the path of water from storm drains is also reported. The overflow water was collected in the inflow sump (IN) and pumped in the BP; from the outflow sump (OUT), the water was pumped in a channel of the urban surface drainage system.   . The rainwater drainage system with the path of water from storm drains is also reported. The overflow water was collected in the inflow sump (IN) and pumped in the BP; from the outflow sump (OUT), the water was pumped in a channel of the urban surface drainage system.  Data collection considered the capacity of the BP in managing the overflow volumes from the existing rainwater drainage system subtracted to the canal of the urban drainage system. The inflow and outflow volumes were calculated on the basis of the operating time of the two drainage pumps (inflow pump and outflow pump) (Submersible pump MC/50-70, Pedrollo S.p.A., San Bonifacio, Italy), knowing their flow rate (1600 L min −1 ).     Data collection considered the capacity of the BP in managing the overflow volumes from the existing rainwater drainage system subtracted to the canal of the urban drainage system. The inflow and outflow volumes were calculated on the basis of the operating time of the two drainage pumps (inflow pump and outflow pump) (Submersible pump MC/50-70, Pedrollo S.p.A., San Bonifacio, Italy), knowing their flow rate (1600 L min −1 ).
The average daily actual evapotranspiration of the system was estimated in order to evaluate the capacity of the pond to ensure an adequate reservoir of water especially during dry periods. The values were calculated by measuring with a water level sensor (Levelogger Edge, Solinist Ltd., Georgetown, ON, Canada) the lowering of the water level in the pond during dry periods in different seasons. The values were compared with the average daily reference evapotranspiration ET 0 values calculated by the Penman-Monteith formula [43] using the data of the local weather station.
Samples of water were collected in the inflow and in the outflow sumps whenever the corresponding hour-counter revealed that pumps had operated. In these occasions, also three samples of water were collected from the pond at 20 cm of depth. Furthermore, samples of the BP water were collected every three weeks with no rainfall event. Water samples were analyzed for the concentration of nutrients and other ions (i.e., Cl − , NO 3− , PO 4 3− , SO 4 2− , Na + , NH 4 + , K + , Ca + , and Mg + ), salinity, and dissolved heavy metals (i.e., Cu, Cb, Zn, and Pb). Nutrients were evaluated by means of ionic chromatography (ICS-900, Dionex, Sunnyvale, CA, USA); salinity, pH and heavy metal were determined with ICP SPECTRO CirOS Vision EOP (SPECTRO Analytical Instruments GmbH & Co., Kleve, Germany). The plant characteristics at planting (i.e., height, leaf number, root length, dry weight) were determined in a sample of four plants per species. On November 2011, plant growth was evaluated by recording only in vivo parameters such as stem and leaf number, height, plant survival. The root growth was evaluated by means of its length and a visual rating (root visual rating RVR; 1-9 scale) based on root number and overall root growth. Furthermore, a visual rating was also adopted to evaluate the aesthetic values of the plants (aesthetic visual rating AVR; 1-9 scale) based on their potential growth in conventional condition. At the end of the second vegetative season (November 2012), the plants were evaluated as previously and, in addition, on a half of the plants, the dry weight of above-ground plant organs (AGPO) and below-ground plant organs (BGPO), comprehensive of rhizomes and stolons when present, was determined. For plants grown on substrate, the concentrations of Cd, Cu, Pb, and Zn in the dry matter were also determined, adopting [44] procedure for mineralization and ICP procedure for reading. The data were also used to calculate heavy metal content in both AGPO and BGPO, multiplying the concentrations by their respective dry weights.
The data were analyzed by mean of the analysis of variance. Statgraphics Centurion XVI software program (Statpoint Technologies, Inc., Warrenton, VA, USA) was used for data analysis. The data from the analysis of water collected in the BP were averaged before statistical analysis. The data on plant survival were analyzed by mean of the chi-square test. Non-linear regression (SigmaPlot for Windows 11.0; Systat Software, Inc., Chicago, IL, USA) was used to describe changes in nutrient concentration over time.

Hydrological Behaviour of the BP
During the period April 2011-November 2012, 121 rainfall events were recorded (a total of 944 mm of rain), but only 14 events generated overflow volumes from the rainwater drainage system. The total inflow volume in the BP was 245 m 3 , and the water volume leaving the pond as outflow was 126 m 3 , corresponding to 119 m 3 (nearly 50% of the total inflow volume) collected or evapotranspired by the BP system. This volume was therefore subtracted from the urban stormwater drainage system, reducing the peak flow rates in the canal during rainy periods. However, it is interesting to note that only 10% of the total potential runoff volume (about 2140 m 3 calculated by multiplying the drainage area for the rainfall) gave rise to overflow volume, because during the examined period the events were mostly of medium-low amount, perfectly managed by the existing rainwater drainage system.
In these occasions, the average daily evapotranspiration, calculated in no rainfall periods, was of 1.01 mm d In the same periods, the average daily reference evapotranspiration ET 0 values calculated by the Penman-Monteith formula [43] were 0.95, 3.31, and 5.21 mm d −1 , respectively. The values of the actual evapotranspiration were relatively high if we consider the low transpiration of newly established plants: this was probably offset by the evaporation from the water surface left free by the floating elements (about 30% of the pond surface) and from the gridded windows without plants.
The BP was meant to guarantee a sufficient water depth for plant growth also in high-drought condition. However, even if rainfall events during the analyzed period were not frequent, the water depth in BP was high (over 100 cm deep) for most part of the experimental time. Only during summer 2012, the water level reached the lowest level (0.83 m) on 31 August. Nevertheless, as the average daily reference evapotranspiration ET 0 value during that summer was equal to 5.21 mm d −1 , considering a crop coefficient Kc equal to 1.2, as set for reed swamp [43] with a good plant growth, the actual ET of the system during the driest summer period could be up to 450 mm compared to the actual 285 mm that we observed. The hypothetical higher evapotranspiration would have resulted in lower depth of the storage water (about 0.60 m), which would have allowed the survival of the plants, demonstrating an adequate sizing of the BP.

Nutrient and Heavy Metal Concentration in Stormwater
The concentrations of nutrients and other ions in the BP inflow and outflow water were in general very low. Cl − , Na + , Mg + , and Ca + concentrations in the BP water did not change over time (on average, 0.59, 12.4, 1.30, and 11.2 mg L −1 , respectively), and no difference between inflow and outflow water was found. NO 3 − concentration in the BP water was relatively high before plant establishment (1.53 mg L −1 on average until July 2011) and, after plant establishment, it was significantly lower (on average 0.355 mg L −1 ). From April 2012, the values were sometimes lower than the limit of instrument detection (<0.02 mg L −1 ). Figure 5, reporting the box and whisker diagrams of nitrate concentration in the BP inflow and outflow water throughout the experimental period, highlights that the outflow water had, as expected, the same concentration as in the BP water but much lower than in the inflow water.
Land 2018, 7, x FOR PEER REVIEW 6 of 16 Monteith formula [43] were 0.95, 3.31, and 5.21 mm d −1 , respectively. The values of the actual evapotranspiration were relatively high if we consider the low transpiration of newly established plants: this was probably offset by the evaporation from the water surface left free by the floating elements (about 30% of the pond surface) and from the gridded windows without plants.
The BP was meant to guarantee a sufficient water depth for plant growth also in high-drought condition. However, even if rainfall events during the analyzed period were not frequent, the water depth in BP was high (over 100 cm deep) for most part of the experimental time. Only during summer 2012, the water level reached the lowest level (0.83 m) on 31 August. Nevertheless, as the average daily reference evapotranspiration ET0 value during that summer was equal to 5.21 mm d −1 , considering a crop coefficient Kc equal to 1.2, as set for reed swamp [43] with a good plant growth, the actual ET of the system during the driest summer period could be up to 450 mm compared to the actual 285 mm that we observed. The hypothetical higher evapotranspiration would have resulted in lower depth of the storage water (about 0.60 m), which would have allowed the survival of the plants, demonstrating an adequate sizing of the BP.

Nutrient and Heavy Metal Concentration in Stormwater
The concentrations of nutrients and other ions in the BP inflow and outflow water were in general very low. Cl − , Na + , Mg + , and Ca + concentrations in the BP water did not change over time (on average, 0.59, 12.4, 1.30, and 11.2 mg L −1 , respectively), and no difference between inflow and outflow water was found. NO3 − concentration in the BP water was relatively high before plant establishment (1.53 mg L −1 on average until July 2011) and, after plant establishment, it was significantly lower (on average 0.355 mg L −1 ). From April 2012, the values were sometimes lower than the limit of instrument detection (<0.02 mg L −1 ). Figure 5, reporting the box and whisker diagrams of nitrate concentration in the BP inflow and outflow water throughout the experimental period, highlights that the outflow water had, as expected, the same concentration as in the BP water but much lower than in the inflow water. In contrast to what reported for NO3 − , no differences were noted among inflow, outflow, and BP water samples. SO4 2− , NH4 + , and K + concentrations had the same pattern described for NO3 − , but differences resulted significant only for K + . As reported in Figure 6, K + concentrations in inflow water were higher than those in the BP, while outflow water had intermediate values. These low values of nutrients and other what reported for NO 3 − , no differences were noted among inflow, outflow, and BP water samples. SO 4 2− , NH 4 + , and K + concentrations had the same pattern described for NO 3 − , but differences resulted significant only for K + . As reported in Figure 6, K + concentrations in inflow water were higher than those in the BP, while outflow water had intermediate values.
These low values of nutrients and other ions are more comparable with those of rainwater and runoff from roofs than with those of the runoff from trafficked areas [45]. ions are more comparable with those of rainwater and runoff from roofs than with those of the runoff from trafficked areas [45]. The water pH also did not change over time, and the values (ranging from 7.02 to 7.62) were similar to those of rainwater [45,46].
The concentrations of dissolved heavy metals in the water were in general very low and, in the case of Cd and Pb, the values were below the sensibility of the instrument (<0.001 and <0.005 mg L −1 , respectively). Concentrations of Cu and Zn did not apparently change over time or in inflow or outflow waters and averaged 0.007 mg L −1 (range 0.002-0.012 mg L −1 ) and 0.004 mg L −1 (range 0.001-0.028 mg L −1 ), respectively. As seen for ions, the concentrations of heavy metals were more comparable with those of rainwater and runoff from roofs than with those of the runoff from trafficked areas [45,47] or at the outlet of storm sewers [46].
The FTWs have proven to be efficient in ameliorating stormwater quality both at mesocorm experiment level [48] and at field level [41,42,49,50], leading to improvements as high as 14% for NO3-N, 65-75% for Cu, and 40% for Zn after seven days of treatment. In the present study, the apparently low or no effect was probably due to the low concentrations of both nutrients and heavy metals in the inlet water, little surface coverage, and growth stage of the plants.

Plant Growth and Heavy Metal Accumulation
The characteristics of the plants at planting and their performance at the end of the first and second year are reported in Tables 1, 2, 3 and 4. A. parviflora (ALSSU) plant material used for transplant arrived from a traditional cultivation in soil. Nevertheless, all plants were alive at the end of the first year, and a small growth occurred, without any difference among treatments. At the end of the experiment, differences were found only for the biomass of AGPO and BGPO and, of course, of the whole plant (WP) ( Table 1). The presence of the substrate improved ALSSU growth and, in fact, the related values were almost doubled ( Table 1).
The plant material utilized for the transplant of B. caroliniana (BAOCA) was merely stem cutting with preformed roots. The presence of a substrate facilitated plant establishment. Furthermore, the plants had 90% more developed buds and a better overall appearance ( Table 2). The plants of BAOCA grown without substrate died during winter, while two-thirds of plants cultivated in the substrate The water pH also did not change over time, and the values (ranging from 7.02 to 7.62) were similar to those of rainwater [45,46].
The concentrations of dissolved heavy metals in the water were in general very low and, in the case of Cd and Pb, the values were below the sensibility of the instrument (<0.001 and <0.005 mg L −1 , respectively). Concentrations of Cu and Zn did not apparently change over time or in inflow or outflow waters and averaged 0.007 mg L −1 (range 0.002-0.012 mg L −1 ) and 0.004 mg L −1 (range 0.001-0.028 mg L −1 ), respectively. As seen for ions, the concentrations of heavy metals were more comparable with those of rainwater and runoff from roofs than with those of the runoff from trafficked areas [45,47] or at the outlet of storm sewers [46].
The FTWs have proven to be efficient in ameliorating stormwater quality both at mesocorm experiment level [48] and at field level [41,42,49,50], leading to improvements as high as 14% for NO 3 -N, 65-75% for Cu, and 40% for Zn after seven days of treatment. In the present study, the apparently low or no effect was probably due to the low concentrations of both nutrients and heavy metals in the inlet water, little surface coverage, and growth stage of the plants.

Plant Growth and Heavy Metal Accumulation
The characteristics of the plants at planting and their performance at the end of the first and second year are reported in Tables 1-4. A. parviflora (ALSSU) plant material used for transplant arrived from a traditional cultivation in soil. Nevertheless, all plants were alive at the end of the first year, and a small growth occurred, without any difference among treatments. At the end of the experiment, differences were found only for the biomass of AGPO and BGPO and, of course, of the whole plant (WP) ( Table 1). The presence of the substrate improved ALSSU growth and, in fact, the related values were almost doubled ( Table 1).
The plant material utilized for the transplant of B. caroliniana (BAOCA) was merely stem cutting with preformed roots. The presence of a substrate facilitated plant establishment. Furthermore, the plants had 90% more developed buds and a better overall appearance ( Table 2). The plants of BAOCA grown without substrate died during winter, while two-thirds of plants cultivated in the substrate remained alive. However, the biomass data collected at the end of the second year indicate that growth was very poor (Table 1). Table 1. Dry weight of above-ground plant organs (AGPO), below-ground plant organs (BGPO), and whole plant (WP) at the beginning and at the end of the experiment (g plant −1 ± sd). All plants of C. palustris (CTAPA) survived and, at the end of the first year, the plants grown with substrate showed better parameters, with the only exception of root length (Table 2). At the end of the experiment, the improvement shown by CTAPA grown with substrate was less evident (Tables 1  and 2).

Beginning of the Experiment
Rhizome cuttings with shoot and poor root system were used for C. indica (CNNIN). At the end of the first year, the plants grown with substrate differed from those without it only in the root system features: the former had shorter but more numerous roots (Table 2). Despite the American Horticultural Society considers CNNIN quite hardy for our winter temperature (www.ahs.org), the plants did not survive through the winter (minimum temperature registered −7.4 • C). Table 2. Plants characteristics at the beginning of the experiment, at the end of the first year, and at the end of the second year (RVR = root visual score; AVR = aesthetic visual score). The plant material used for Iris (IRISS) and L. punctata (LYSPU) transplant was not adapted to the floating system, but no death was observed for both species (Table 3). Some parameters were significantly higher in plants grown with substrate (e.g., plant height of IRISS and stem number of LYSPU), but no differences were observed in biomass accumulation ( Table 1).

End of the First Year
The plants of M. aquatica (MENAQ) and P. arundinacea (PHAAP) showed a better growth with substrate already at the end of the first year, with significant higher values for all the observed parameters (see Tables 1, 3 and 4). These species responded similarly also in plant survival, with 75 and a 50% of death for the two species, occurred only in absence of substrate.
Plants of O. javanica (OENJA) had also better performance in the presence of substrate (Table 3), and, at the end of the experiment, only the plants grown with the substrate were still alive.
The same behavior was observed for S. sagittifolia (SAGSA) and T. laxmannii (THYLX) plants. During the first year, a 16.7% of SAGSA plants growth without substrate died, and the remaining did not survive winter. However, a very poor growth was observed also for plants grown in the substrate (Tables 1 and 4). Regarding TYHLX, 91.7% of plants grown with the substrate survived, while only 58.3% of those grown without the substrate died ( Table 4). The growth of the remaining plants was greatly improved by the substrate and, with the only exception of root length, all parameters showed values increased by over 100% (Tables 1 and 4). From these results, it is clear that both survival and growth of plants in the Tech_IA ® elements were promoted by the presence of a substrate (Figure 7). This was probably due to a better root environment (e.g., humidity around the collar point) which initially favoured rooting and promoted plant establishment. Apart from that, plants like SAGSA, with clumping habit and weak structures, need a suitable anchorage for their growth. On the contrary, the attitude to produce stolons or rhizomes (e.g., Mentha and Iris) favours anchorage to Tech-IA ® elements.  From these results, it is clear that both survival and growth of plants in the Tech_IA ® elements were promoted by the presence of a substrate (Figure 7). This was probably due to a better root environment (e.g., humidity around the collar point) which initially favoured rooting and promoted plant establishment. Apart from that, plants like SAGSA, with clumping habit and weak structures, need a suitable anchorage for their growth. On the contrary, the attitude to produce stolons or rhizomes (e.g., Mentha and Iris) favours anchorage to Tech-IA ® elements.  When evaluating species survival to winter, it appeared that C. indica is not useful in our environment, as well as Bacopa and Sagittaria, whose poor growth indicates their poor adaptability to the employed BP system.
In order to evaluate the ability of the selected plants to improve water quality, the accumulation of heavy metal was also measured. The evaluation considered only the species with good growth results and high plant survival, cultivated with the substrate. As reported in Table 3, the heavy metal with the highest concentration was Zn followed by Cu, Pb, and Cd. Furthermore, in general, higher concentrations were found in the AGPO than in the BGPO.
Among plants species, A. parviflora had the highest values of all heavy metals. On the contrary, Iris had the lowest values of Cu, Pb, and Zn. Mentha, Phalaris, and Typha had the overall lowest values of Cd. The concentration of Zn found in this research is higher than those found by another study [51] in P. arundinacea and other species grown in normal condition, or comparable with those found in plants grown with nutrient solutions. Furthermore, according to reference [52], the concentrations of heavy metals are within normal levels even if the values found in water are relatively low.
If we consider the heavy metals accumulated on a mass basis, the highest values were found in C. palustris and M. aquatica. It is worth noting that ALSSU, as well as MENAQ, CTAPA, LYSPU, and OENJA had a good accumulation of these elements in the AGPO ( Table 5).
As the management of the BP includes an annual cleaning of Tech-IA elements with the removal of the aerial part of plants, the content of heavy metals in AGPO can be the most important factor to consider if the plant selection is done on the basis of its phytoremediation ability. In fact, the heavy metals are removed with the removal of the aerial part of plants, while they remain in BGPO and, because of the decay of the root system, they can be released in the water of the pond.
A last consideration has to be made. In a sustainable approach to manage storm water in an urban context, no nutritional elements were provided to the plants. As the water arriving at the pond was very poor in nutrients, plant growth of all species was poorer than in normal nutritional conditions, as evidenced by the nutritional deficiency symptoms that were observed in all species (Figure 8). It is possible that, if a controlled-release fertilizer was applied on the substrate during active plant growth (i.e., in springtime), the surface growth and aesthetic appearance of the plants could be improved.
When evaluating species survival to winter, it appeared that C. indica is not useful in our environment, as well as Bacopa and Sagittaria, whose poor growth indicates their poor adaptability to the employed BP system.
In order to evaluate the ability of the selected plants to improve water quality, the accumulation of heavy metal was also measured. The evaluation considered only the species with good growth results and high plant survival, cultivated with the substrate. As reported in Table 3, the heavy metal with the highest concentration was Zn followed by Cu, Pb, and Cd. Furthermore, in general, higher concentrations were found in the AGPO than in the BGPO.
Among plants species, A. parviflora had the highest values of all heavy metals. On the contrary, Iris had the lowest values of Cu, Pb, and Zn. Mentha, Phalaris, and Typha had the overall lowest values of Cd. The concentration of Zn found in this research is higher than those found by another study [51] in P. arundinacea and other species grown in normal condition, or comparable with those found in plants grown with nutrient solutions. Furthermore, according to reference [52], the concentrations of heavy metals are within normal levels even if the values found in water are relatively low.
If we consider the heavy metals accumulated on a mass basis, the highest values were found in C. palustris and M. aquatica. It is worth noting that ALSSU, as well as MENAQ, CTAPA, LYSPU, and OENJA had a good accumulation of these elements in the AGPO ( Table 5).
As the management of the BP includes an annual cleaning of Tech-IA elements with the removal of the aerial part of plants, the content of heavy metals in AGPO can be the most important factor to consider if the plant selection is done on the basis of its phytoremediation ability. In fact, the heavy metals are removed with the removal of the aerial part of plants, while they remain in BGPO and, because of the decay of the root system, they can be released in the water of the pond.
A last consideration has to be made. In a sustainable approach to manage storm water in an urban context, no nutritional elements were provided to the plants. As the water arriving at the pond was very poor in nutrients, plant growth of all species was poorer than in normal nutritional conditions, as evidenced by the nutritional deficiency symptoms that were observed in all species (Figure 8). It is possible that, if a controlled-release fertilizer was applied on the substrate during active plant growth (i.e., in springtime), the surface growth and aesthetic appearance of the plants could be improved.

Conflicts of Interest:
The authors declare no conflict of interest.