Table S1. Average Physico-Chemical Parameters of Influent and Effluent of Each Chamber

Table S1. Average physico-chemical parameters of influent and effluent of each chamber of the HMAE®. CW with CW with Iris pseudacorus Pre-treatment Stabilization pond Parameter Inflow (mean ± SE) Typha latifolia and Salix atrocinerea (mean ± SE) (mean ± SE) (mean ± SE) (mean ± SE) Temperature (ºC) 13.52 ± 1.04 13.10 ± 1.06 12.97 ± 1.20 12.70 ± 1.11 11.28 ± 0.97 Conductivity (µS cm- 777.10 ± 73.99 668.95 ± 63.24 602.67 ± 48.56 622.29 ± 47.31 498.65 ± 44.45 1) pH 7.17 ± 0.15 7.35 ± 0.13 7.06 ± 0.19 7.10 ± 0.13 7.36 ± 0.10 Dissolved oxygen 28.33 ± 4.36 59.85 ± 7.27 55.85 ± 6.82 41.61 ± 5.00 69.84 ± 5.64 (mg L-1) Redox potential (mV) -81.41 ± 21.19 -43.75 ± 16.96 -35.31 ± 16.51 -76.16 ± 22.32 -4.62 ± 25.54 Table S2. Average organic matter concentrations (mg L-1), in terms of chemical organic demand (COD), of influent and effluent of each chamber of the HMAE®. CW with Iris Stabilization CW with Typha Yea Influent Pre-treatment pseudacorus and N pond (mean ± latifolia (mean ± r (mean ± SE) (mean ± SE) Salix atrocinerea SE) SE) (mean ± SE) 2000 18 201.68 ± 40.88 87.71 ± 24.39 138.99 ± 17.40 118.67 ± 12.93 86.01 ± 14.61 2001 9 164.11 ± 27.33 106.55 ± 23.00 110.65 ± 13.44 70.28 ± 11.81 61.66 ± 16.52 2002 12 349.15 ± 56.60 151.80 ± 22.48 205.88 ± 19.02 160.77 ± 23.63 100.77 ± 8.39 2003 8 209.46 ± 20.84 149.38 ± 8.38 132.60 ± 17.90 101.29 ± 7.94 81.01 ± 12.09 2004 12 301.24 ± 60.10 160.86 ± 24.10 118.07 ± 10.42 105.30 ± 11.37 75.14 ± 9.83 2005 12 290.44 ± 45.42 132.07 ± 16.49 101.97 ± 15.07 82.10 ± 6.67 58.42 ± 4.48 2006 12 129.87 ± 16.85 110.50 ± 14.78 111.83 ± 20.76 103.54 ± 15.67 63.88 ± 11.28 2007 12 169.51 ± 32.35 118.83 ± 27.24 97.14 ± 20.44 68.63 ± 14.00 29.90 ± 5.64 2008 12 176.12 ± 38.17 95.86 ± 10.76 93.95 ± 15.44 61.94 ± 8.25 42.78 ± 8.06 2009 12 145.47 ± 24.36 142.62 ± 25.49 119.46 ± 13.50 105.34 ± 18.91 26.89 ± 8.98 2010 4 158.59 ± 40.10 173.23 ± 15.08 154.41 ± 62.40 112.69 ± 52.08 49.78 ± 16.83 2011 4 219.44 ± 9.83 146.22 ± 27.81 122.52 ± 12.13 124.149 ± 4.35 57.99 ± 11.06 Table S3. Average organic matter concentrations (mg L-1), in terms of biological

® organic demand (BOD5), of influent and effluent of each chamber of the HMAE . CW with Iris Stabilization CW with Typha Yea Influent (mean Pre-treatment pseudacorus and N pond (mean ± latifolia (mean ± r ± SE) (mean ± SE) Salix atrocinerea SE) SE) (mean ± SE) 2000 18 60.44 ± 15.66 20.01 ± 3.63 31.67 ± 3.22 40.61 ± 5.71 15.94 ± 2.86 2001 9 85.00 ± 28.43 51.78 ± 18.94 35.67 ± 5.48 24.67 ± 5.01 4.95 ± 5.09 2002 12 162.67 ± 18.06 101.42 ± 17.55 92.34 ± 12.37 36.09 ± 6.21 23.56 ± 3.20 2003 8 130.00 ± 16.33 81.79 ± 7.09 62.43 ± 12.76 46.79 ± 7.03 18.00 ± 1.97 2004 12 131.11 ± 25.47 72.34 ± 14.12 42.92 ± 7.62 32.75 ± 5.07 20.09 ± 1.31 2005 12 155.00 ± 30.99 79.17 ± 9.63 52.83 ± 8.91 36.83 ± 6.82 20.18 ± 2.79 2006 12 88.58 ± 12.22 66.08 ± 11.47 59.67 ± 10.33 52.33 ± 8.04 25.75 ± 7.79 2007 12 102.92 ± 23.74 64.17 ± 15.12 44.83 ± 7.81 31.50 ± 6.50 12.58 ± 3.39 2008 12 114.08 ± 25.64 53.67 ± 10.28 29.42 ± 5.11 23.27 ± 6.10 12.08 ± 3.59 2009 12 90.50 ± 18.88 58.50 ± 10.99 40.00 ± 8.16 30.01 ± 8.96 8.13 ± 2.39 2010 4 70.02 ± 17.32 60.00 ± 14.15 26.34 ± 5.67 34.75 ± 7.43 9.00 ± 4.36 2011 4 101.67 ± 4.41 78.33 ± 9.28 46.00 ± 6.00 32.00 ± 6.11 17.39 ± 7.68 Table S4. Average suspended solids concentrations (TSS, mg L-1) of influent and effluent of each chamber of the HMAE®.

CW with Stabilization CW with Typha Influent (mean Pre-treatment pseudacorus Year N pond (mean ± latifolia (mean ± ± SE) (mean ± SE) Salix atrocinerea SE) SE) (mean ± SE) 2000 18 92.98 ± 19.16 36.69 ± 8.03 79.14 ± 10.20 76.29 ± 4.71 34.37 ± 7.45 2001 9 54.91 ± 17.57 28.98 ± 7.50 54.70 ± 9.22 43.69 ± 13.72 27.24 ± 12.76 2002 12 303.76 ± 43.50 124.05 ± 18.58 95.42 ± 14.53 51.14 ± 8.46 28.30 ± 8.14 2003 8 120.19 ± 46.98 43.62 ± 10.58 33.38 ± 5.37 33.37 ± 4.74 28.78 ± 6.05 2004 12 181.33 ± 47.32 40.48 ± 7.17 35.10 ± 5.16 26.67 ± 3.44 16.38 ± 2.48 2005 12 54.50 ± 8.40 57.38 ± 11.26 53.36 ± 11.77 58.08 ± 9.81 20.28 ± 4.22 2006 12 66.64 ± 20.39 47.25 ± 11.67 38.75 ± 9.38 32.09 ± 6.61 19.42 ± 6.04 2007 12 121.50 ± 29.33 39.50 ± 8.73 43.67 ± 10.54 39.67 ± 6.23 14.67 ± 5.44 2008 12 84.73 ± 14.94 42.17 ± 2.90 55.33 ± 10.62 45.67 ± 8.67 15.75 ± 2.71 2009 12 122.00 ± 67.59 46.50 ± 8.99 74.75 ± 29.58 53.25 ± 26.64 30.75 ± 5.28 2010 4 136.00 ± 48.01 48.00 ± 17.01 31.33 ± 27.22 69.33 ± 23.13 12.00 ± 5.29 Table S5. Average Kjeldahl nitrogen concentrations (mg L-1) of influent and effluent of each chamber of the HMAE®. CW with Iris Stabilization CW with Typha Yea Influent Pre-treatment pseudacorus and Salix N pond (mean ± latifolia (mean ± r (mean ± SE) (mean ± SE) atrocinerea (mean ± SE) SE) SE) 2000 18 42.23 ± 8.61 35.91 ± 7.13 30.69 ± 4.15 24.86 ± 3.60 20.31 ± 3.01 2001 9 27.76 ± 13.40 15.61 ± 2.96 20.93 ± 3.75 18.85 ± 4.47 20.56 ± 4.19 2002 12 52.47 ± 8.32 39.35 ± 7.51 42.72 ± 6.19 33.87 ± 5.21 36.09 ± 5.91 2003 8 40.69 ± 5.65 35.28 ± 6.52 35.62 ± 3.79 32.02 ± 6.02 26.76 ± 3.93 2004 12 49.60 ± 8.38 40.48 ± 6.43 29.85 ± 3.58 28.83 ± 2.36 23.41 ± 2.18 2005 12 67.25 ± 4.83 49.31 ± 2.60 47.22 ± 3.40 33.73 ± 2.55 26.58 ± 2.50 2006 12 70.31 ± 11.31 47.62 ± 7.43 36.55 ± 4.81 36.33 ± 6.17 35.37 ± 7.51 2007 12 37.60 ± 7.16 29.43 ± 4.39 19.01 ± 2.27 16.01 ± 2.38 15.83 ± 2.40 2008 12 24.11 ± 3.57 28.42 ± 4.88 22.31 ± 3.42 17.62 ± 2.28 13.89 ± 2.17 2009 12 31.43 ± 4.70 31.79 ± 5.26 23.66 ± 2.08 21.96 ± 2.32 15.87 ± 1.82 2010 4 30.80 ± 3.56 35.83 ± 2.74 33.06 ± 2.38 23.85 ± 5.68 10.37 ± 2.49 2011 4 50.75 ± 12.42 50.65 ± 12.46 40.56 ± 6.78 40.60 ± 5.77 27.56 ± 4.07 Table S6. Average Ammonium nitrogen concentrations (mg L-1) of influent and effluent of each chamber of the HMAE®.

CW with Stabilization CW with Typha Influent Pre-treatment pseudacorus Year N pond (mean ± latifolia (mean ± (mean ± SE) (mean ± SE) Salix atrocinerea SE) SE) (mean ± SE) 2000 18 36.05 ± 7.32 24.72 ± 6.15 16.08 ± 3.31 13.15 ± 3.39 11.23 ± 2.40 2001 9 17.84 ± 6.51 12.87 ± 3.01 16.79 ± 3.31 15.34 ± 3.72 15.89 ± 3.66 2002 12 37.54 ± 6.26 25.71 ± 5.35 24.41 ± 4.46 18.48 ± 4.29 18.89 ± 3.79 2003 8 24.71 ± 6.50 23.86 ± 7.22 16.82 ± 4.80 12.98 ± 5.61 7.61 ± 3.50 2004 12 32.37 ± 5.81 27.80 ± 4.68 20.70 ± 3.15 19.57 ± 2.92 14.28 ± 2.40 2005 12 52.30 ± 3.58 40.94 ± 4.47 34.60 ± 2.84 25.26 ± 2.043 19.52 ± 1.79 2006 12 50.58 ± 11.56 29.28 ± 6.02 23.15 ± 3.61 24.10 ± 5.81 16.01 ± 4.53 2007 12 23.84 ± 5.31 17.55 ± 3.28 12.66 ± 2.25 11.63 ± 1.79 10.46 ± 1.96 2008 12 18.70 ± 3.35 19.03 ± 4.08 13.82 ± 2.58 14.18 ± 2.27 11.18 ± 1.74 2009 12 21.55 ± 2.68 26.91 ± 3.62 16.70 ± 1.90 15.64 ± 1.32 14.60 ± 1.40 2010 4 15.79 ± 3.80 17.14 ± 5.54 11.43 ± 3.01 8.88 ± 2.95 3.49 ± 1.76 2011 4 40.58 ± 11.00 39.55 ± 10.85 33.09 ± 9.34 32.86 ± 3.86 21.22 ± 3.87 Table S7. Average total phosphorous concentrations (mg L-1) of influent and effluent of each chamber of the HMAE®.

CW with Stabilization CW with Typha Yea Influent Pre-treatment pseudacorus N pond (mean ± latifolia (mean ± r (mean ± SE) (mean ± SE) atrocinerea SE) SE) 2000 18 7.77 ± 0.92 7.28 ± 0.89 7.61 ± 0.67 7.43 ± 0.67 6.80 ± 0.42 2001 9 13.56 ± 2.92 11.97 ± 2.40 13.23 ± 3.60 2.55 ± 0.18 4.09 ± 0.54 2002 12 7.84 ± 0.75 6.50 ± 0.59 6.83 ± 0.45 6.64 ± 0.57 5.95 ± 0.56 2003 8 9.19 ± 0.56 7.56 ± 0.93 6.80 ± 0.84 6.73 ± 0.75 5.70 ± 0.74 2004 12 6.28 ± 0.76 4.85 ± 0.46 4.35 ± 0.42 4.31 ± 0.32 3.67 ± 0.49 2005 12 7.99 ± 0.86 6.93 ± 0.68 5.53 ± 0.44 4.73 ± 0.32 3.65 ± 0.39 2006 12 8.69 ± 0.93 7.05 ± 1.10 7.96 ± 1.57 3.37 ± 0.62 4.18 ± 0.58 2007 12 6.21 ± 1.31 6.17 ± 1.70 4.05 ± 0.60 3.77 ± 0.59 3.29 ± 0.72 2008 12 2.28 ± 0.32 1.95 ± 0.24 2.03 ± 0.23 2.19 ± 0.36 2.00 ± 0.19 2009 12 3.83 ± 1.10 4.26 ± 1.19 4.00 ± 0.85 4.01 ± 0.98 2.38 ± 0.55 2010 4 2.05 ± 0.68 2.25 ± 0.62 2.28 ± 0.39 2.15 ± 0.43 1.06 ± 0.45 2011 4 4.67 ± 0.37 4.43 ± 0.83 4.13 ± 0.66 4.43 ± 0.28 3.53 ± 0.49