Skip to content
BY-NC-ND 3.0 license Open Access Published by De Gruyter Open Access January 3, 2016

Effects of the introduction of pre-treated wastewater in a shallow lake reed stand

  • Mária Dinka EMAIL logo , Anita Kiss , Norbert Magyar and Edit Ágoston-Szabó
From the journal Open Geosciences


Reed stands may be employed in the amelioration of water quality or even in the treatment of wastewater. In this study, the nutrient concentrations of (i) the above- and below-ground Common Reed (Phragmites australis) biomass, and (ii) surface and interstitial water were analyzed in a natural stand used in wastewater treatment. The reed stand was located in Hungarian part of Lake Fertő/Neusiedler See, by the shore near Fertőrákos Bay. The nitrate, phosphate and dissolved organic nitrogen concentrations of surface water were found to be higher on the inlet side of the reed stand compared to the outlet. The N and P concentrations in the above-ground biomass and P concentrations in the below-ground biomass increased after the introduction of pre-treated wastewater. The inter-annual differences in the characteristics of sediment interstitial water and in the nutrient content of reed tissues were assessed using statistical methods. The samples taken before and after the introduction of the pre-treated wastewater in the parcel formed different clusters. The results of the study provide further evidence that the nutrient retention capacity of natural stands of P. australis may be employed in the treatment of wastewater while protecting and preserving the valuable natural assets of the lake.

1 Introduction

Phragmites australis (Cav.) Trin ex Steudel is one of the most prolific and widespread emergent macrophytes in the world. While it is considered an invasive species in some places [1], it has also been shown to provide numerous cultural benefits, as well as being an important species in many wetland habitats [213].

It has functionally adapted to anaerobic conditions in sediments, on account of its ability to translocate oxygen into the rhizosphere via its well-developed aerenchima and to oxygenize and increase the redox potential of the otherwise anaerobic surrounding sediment [14, 15]. In this way, it can create favorable conditions for the aerobic and facultative anaerobic microorganisms in its root zone [16].

P. australis fulfills an important role in water purification because of its high filtering capacity: i.e. by bioremediation bacterial action on the surface of the roots, and by the uptake and incorporation of nutrients into its own biomass (which can then be partially removed by harvesting). Thanks to these characteristics, P. australis can be successfully used in phytoremediation water treatment [1721]. However, high water-borne nutrient levels represent stressful conditions for reed stands and may affect the metabolism and the growth of the reed, leading in the end to the decline of reed-stands [22, 23].

The efficiency of P. australis in water purification and nutrient removal, along with its usefulness in the treatment of domestic and agricultural wastewater, has been thoroughly described, especially regarding the use of the Common Reed in constructed wetlands [11, 19, 2426].

The aim of the study therefore was to determine the potential capacity of a natural stand of P. australis to assimilate nutrients from pre-treated wastewater in the reed stands on the western shore of Lake Fertő/Neusiedler See. The approach involved the measurement of (i) the changes in nutrient concentrations in surface and interstitial water at locations within the stand of Common Reeds and at the inlet and outlet of the pre-treated wastewater (ii) the nutrient content of the above- and below-ground biomass before and after the introduction of the wastewater.

2 Materials and Methods

2.1 Site description

The study site is located in Lake Fertő/Neusiedler See, a large, shallow, reed-dominated lake, situated on the Hungarian-Austrian border with a surface area of 309 km2 and with a regulated outflow. Common Reed covers 54% of the lake, and 85% of the Hungarian part (75 km2). P. australis, is the dominant species in the littoral zone and it has been shown to play an important role in the biogeochemical nutrient cycle of the lake through the uptake, storage and decomposition [2731]. The lake has been declared a biosphere reserve by UNESCO, as well as a Ramsar wetland, National Park, and World Heritage site. A detailed description of the lake is given in Löffler [32].

The landward edge of a wide inhomogeneous reed stand–the area of the present study-is situated near Fertőrákos Bay on the lake. A 7.4 ha area of the natural reed stand was set aside for the subsequent cleaning of the pre-treated wastewater (bottom left: 47.715993, 16.667188; bottom right: 47.716195, 16.668347; top left: 47.721653, 16.665602; top right: 47.722043, 16.666504). The parcel was provided with hydro-mechanical devices regulating the introduction of the water (Fig. 1).

Figure 1 a) The Hungarian part of Lake Fertő/Neusiedler See. (area in blue: open water and inner ponds in the reed stands; area in white: reed stands with thin lines representing the canal system). b) Sketch of the filter bed delineated in the natural reed stand and used for the treatment of pre-treated wastewater, with sampling sites (PR: sites closed to the inlet, PV: sites close to the outlet, FR: reference site; black points: sampling sites for surface and interstitial water; black squares: reed sampling sites).
Figure 1

a) The Hungarian part of Lake Fertő/Neusiedler See. (area in blue: open water and inner ponds in the reed stands; area in white: reed stands with thin lines representing the canal system). b) Sketch of the filter bed delineated in the natural reed stand and used for the treatment of pre-treated wastewater, with sampling sites (PR: sites closed to the inlet, PV: sites close to the outlet, FR: reference site; black points: sampling sites for surface and interstitial water; black squares: reed sampling sites).

From May to October a dosing rate of 300 m3 day−1, while from November to April 250 m3 day−1 of pre-treated effluent was discharged to the reed bed. As a consequence of this dosing, the water depth over the wetland increased to 5 and 25 cm in summer and winter, respectively. Water flow across the reed bed was not uniform, because the soil surface of the bed was uneven. The difference in height across the parcel in the direction of the outlet was about 3 cm in every hundred meters [33].

The sampling sites (PR and PV; Fig. 1) were located 20 m from the inlet and outlet at each site. We established sampling plots at 4, 8, and 16 m from the edge of the reed bed to the interior in order to sample the interstitial water. The surface water samples were taken from the inlet, outlet, and at 4, 16 m at PR and PV sites (Fig. 1). At the same time, a reference site was selected at Fertőrákos Bay (FR) for the comparison of surface water characteristics. The reed samples were taken from the PR and PV sites close to the surface and interstitial sampling sites at 4 and 16 m, respectively (Fig. 1).

2.2 Surface water and interstitial water

Surface and interstitial water samples were collected on four dates in 2004, both before (April 4), and after the introduction of the pre-treated wastewater (June 13, August 1, October 17). Interstitial water was collected from depths of 0–20 and 20–40 cm using triplicate PVC tubes (5 cm diameter, previously purged) the lower 20 cm of which was perforated. Reducing interaction time of the samples with air to a minimum was a priority in the sampling process. Subsamples were filtered in the field for certain laboratory analyses and ZnCl2 was added to the sub-samples taken for S2− determination. The temperature, pH, Eh and oxygen concentration of the surface and interstitial water samples was determined in situ with a Hydrolog 2100 field instrument (Grabner, Vienna).

The samples were transported to the laboratory in a cooler box and stored in a refrigerator until the chemical analyses, which began the following day. The NH4+,NO2, NO3, and SO42 concentrations were determined using a Dionex DX-120 ionchromatograph, after filtering the samples at 0.45 μm in the field and then 0.2 μm in the laboratory (using Chromafil filters). Dissolved organic (DOC), inorganic (DIC), total (DTC) carbon, and dissolved total nitrogen (DN) concentrations were measured using a LiquiTOC analyser, from previously filtered samples (using 0.45 μm Chromafil filters). The PO43 was determined spectrophotometrically, using the molybdenum blue method. The S2− concentrations of surface and interstitial water were measured employing iodometric titration [34].

2.3 Reeds

The reeds were sampled on two dates (in August 2003 and August 2004) to determine shoot length, density, inflorescence, basal diameter, living leaf (number, surface), leaf area index (LAI), shoot dry mass, and biomass using the procedures briefly described in the present study. More information can be found in [3] and [35].

The Reeds were harvested from randomly selected 0.25 m2 quadrats at each site in 4 replicates by cutting shoots at the sediment surface. The current year above- ground biomass and the standing litter (old shoots) were determined separately. Eight samples of rhizomes were collected using 1 m long sampling tubes with a diameter of 19.5 cm [36, 37] at the PR and PV sites separately. The rhizomes were separated into living, senescing, and dead categories based on the colour, branching and consistency of the rhizome system.

The C, N, S concentrations of the above-ground and below-ground biomass were determined using a Fisons NA1500 NCS-analyser, and the P concentrations with the molybdenum blue method after the digestion with sulphuric acid.

2.4 Data analysis

Besides the calculation of the descriptive statistics of the datasets, uni-and multivariate statistical methods were applied to investigate the differences between the sampling sites in space and time. The significant differences between the group’s means (i) in time for all the sampling events at one site, and (ii) in space for one sampling time and all the sites were investigated using One-way ANOVA. The homogeneity of variances was assessed using F-tests. In the case of both tests a p = 0.05 significance level was applied [38, 39].

To answer the question of which sampling sites and times are similar to each other, hierarchical cluster analysis was used with Ward’s method [40] and squared Euclidean distance. This is a widely known multivariate classification method in which each case starts in a separate cluster and joins up to the other clusters as the linkage distance grows, until only one cluster remains [41]. This method has been successfully applied in hydrology [4245], hydrogeology [46], geology [4751], chemistry [52] and anthropology [53] to find similar and homogeneous groups of observations.

The validity of the groupings was verified using linear discriminant analysis (LDA), which separates the observations with linear planes resulting in a percentage of correctly classified cases [54, 55]. During the analysis, a predictive model was built for group membership. The model consists of discriminant functions based on linear combinations of the predictor variables [44].

The groupings of the sampling sites were also assessed based on the measured variables’ stochastic connections using a powerful pattern recognition technique, principal component analysis (PCA) [56]. The principal components are uncorrelated and are obtained as a linear combination of the original variables [57]:


where z is the component loading, a is the component score, x is the measured value of the variable, i is the component number, j is the sample number, and m is the total number of variables [58]. The pattern of the observations’ grouping can be visualized by plotting the scores of the first and second principal components on a scatter-plot [59].

The statistical analyses were performed using R Software. Besides number of basic functions of base and stats packages, lda from MASS and PCA from FactoMineR package were used during the analysis.

3 Results

3.1 Surface and interstitial water

3.1.1 Overview

The chemical characteristics of surface and interstitial water were found to be quite diverse (Table 1). The pH of the surface water fell between 7.0 and 8.3, with the highest values measured at the inlet (8.3) and outlet (8.2). The pH values at the reed bed sites were lower than the pH values measured in surface water at the reference site (FR: 8.8–9.3). The pH of the interstitial water was lower than that of the surface water in each case.

Table 1

Chemical characteristics of surface and interstitial water in the investigated reed parcel (PR: sites close to the inlet, PV: sites close to the outlet, FR: reference site; sites for sampling of the surface and interstitial water: PR 4, 8, 16 m, PV 4, 8, 16 m from the edge of the parcel).

Sampling sitessurface waterinterstitial water
ConductivityµS cm−1min14641103154792290579022271219133616761316.513191086.5
DOmg l−1min0.250.640.80.561.
Oxygen saturation%min2.36.685.911.520.87757.
S2−mg l−1min<LOD<LOD4.071.651.53<LOD<LOD0.54.715.260.631.531.65
PO43−mg l−1min60.344.480.
Clmg l−1minn.a.70.43159.1250.3247.85n.a.261.278.0988.5987.3357.3660.1457.7
NO3mg l−1min0.4040.040.2690.040.0870.1960.1210.1630.1180.1070.1080.0860.129
SO42−mg l−1min196.2105.54179.646.476.0810.02365.518.9945.5717.712.813.253.84
Na+mg l−1min245.06122.23229.459.5255.4173.65104.82119.76154.48135.7125.16105.9479.87
K+mg l−1min26.8512.0125.1412.3111.2711.4414.316.59.614.78.3610.6211.48
Mg2+mg l−1min26.3154.5636.9161.4555.7962.3690.657.6548.0860.58113.65116.27109.26
Ca2+mg l−1min78.94120.8689.5287.5869.5876.7218.7796.71114.44128.79134.39138.07117.44
DICmg l−1min88.3490.95100.3394.37101.7299.41120.71111.68135.79144.23172.5173.69145.33
DOCmg l−1min17.5119.2917.8919.6621.7829.3919.0221.5221.8227.6231.1434.8236.04
DNmg l−1min3.191.992.442.

The electrical conductivity of surface water decreased proceeding from the inlet (1464–1965 μS cm−1) towards the outlet (790–1184 μS cm−1). Most of the time the electrical conductivity of interstitial water increased with sediment depth, and EC values were always higher in surface water at the reference site (Table 1).

The values for dissolved oxygen (DO) concentration at the inlet and outlet did not differ remarkably from each other, and in most cases were higher than in the reed stand. The oxygen concentration also decreased from the edge (4 m) towards to interior (16 m) of the parcel. The DO concentration in the interstitial water was lower than in the surface water in most cases. At reference site, the DO concentration of the surface water was 4-5 times higher than in the reed parcel.

Higher PO43 concentrations were measured at the inlet (6.0–15.25mg l−1) than at the outlet (0.17–1.17mg l−1). In the reed stand the PO43 concentrations of the surface and interstitial water were higher at the PR sites in most cases. The phosphate concentration in surface water at the reference site was almost always lower (avg. 0.12 mg l−1) than the values measured in the reed stand.

The SO42 concentration in the surface water was higher at the inlet (range = 196.2–306.7 mg l−1) than at the outlet (range = 10.0–41.7mg l−1; Table 1). The SO42 concentrations of surface and interstitial water were both higher at the PR sites. The SO42 concentrations of surface water were almost always higher than in the interstitial water, and these concentrations decreased with depth. The SO42 concentrations were highest at the reference site (Table 1).

The concentration of S2− was always beneath the limit of detectability in the surface water at the inlet, outlet, or even at the reference site. However, in surface water in the reed stand it varied between 1.5 and 12.7 mg l−1. It was almost always lower than in the interstitial water (Table 1).

Ammonium and NO2 were not detectable in either surface or interstitial water. Nitrate concentrations were higher at the inlet compared to the outlet (Table 1) and decreased from the inlet to the outlet, but there was no consistent pattern across the distances and sites when the four sampling dates were compared. At the reference site, the NO3 concentrations were always lower than the values measured at the outlet (Table 1).

The DN (total dissolved nitrogen) of the surface water at the inlet varied between 3.19 and 6.77 mg l−1 and it was higher than at the outlet (2.3–2.7 mg l−1). DN concentrations in the surface water within the reed bed were highest at the PR sites and decreased from the inlet towards the outlet. The DN in the surface water at the outlet was almost the same as that at the reference site (FR range = 1.50–2.62 mg l−1).

Higher DOC concentrations were measured in the surface water at the outlet compared to the inlet. The DOC concentration of the interstitial water increased with the depth and it was higher than that in surface water.

3.1.2 Multivariate results

The cluster results of the surface water samples indicated that the samples taken before and after the introduction of the pre-treated wastewater in the parcel (PR and PV) formed two different groups, regardless of which site they came from, but in the meanwhile this process did not affect the reference site in the open water of Fertőrákos Bay (Fig. 2a). Here, the samples remained in one group before and after the introduction of pre-treated wastewater. With the grouping of the interstitial water samples, it became clear that before the introduction of the pre-treated wastewater all the samples from all the sites formed one group, while after it, these separated according to which site they were measured at, near the outlet and the inlet, PV and PR, respectively (Fig. 2b). These groupings were justified by the results of the linear discriminant analysis.

Figure 2 Grouping of the sampling sites based on the results of the cluster analysis conducted on the surface (a) and interstitial water samples (b).
Figure 2

Grouping of the sampling sites based on the results of the cluster analysis conducted on the surface (a) and interstitial water samples (b).

Using PCA the principal component (PC) loadings and scores were obtained for the interstitial water samples. The first two PCs explained 62–74% of the total variance. The scores of these two PCs were plotted against each other on scatterplots in the case of all sampling events. These plots confirmed the results of the cluster analysis; the samples originating from the sites near the inlet (PR) and the outlet (PV) separated from each other (A1) in spite of the fact that not the total variance of the dataset but only a part of it was considered during the analysis.

3.2 Characteristics of the Common Reed

3.2.1 Above ground

Regarding the point when the above-ground biomass reached its maximum, shoot height varied between sites and from 2003 to 2004 (Table 2). The average shoot height was smaller in most cases at the PV than at the PR site. The differences in shoot basal diameter between the PR and PV sites were notable in both years (Table 2). The number of living leaves varied between 7 and 10 per shoot at both sites; the only year when they differed significantly was 2004. The morphometry of the shoots (shorter and thinner shoots) indicated a change in 2004, leading to a considerable decrease in dry mass. This also affected biomass as well (Table 2). As a combined result of shoot density increase and the previously discussed decrease in shoot morphometry, a decrease in biomass was recorded at the PV sites (Table 2). Considerable differences in the assimilating surface were recorded between the sites and from year to year (Table 2). The leaf area index (LAI) was higher at the PR than at the PV sites.

Table 2

Biometrical parameters and biomass of the Phragmites australis (PR: sites close to the inlet, PV: sites close to the outlet)

(* represents the significant differences (p = 0.05) between the investigated years at the same sampling sites; the letters indicate the significant differences between the sampling sites in the given year; > marks the significant differences between the investigated years; values in bold indicate that the sampling site differed significantly from another site(s) in the investigated year or from the values from among the same sites in time).

Parameters/Sampling sitesUnit20032004
densityNo m−2847254.766.7921189264
inflorescenceNo m−29.3162818.710161025
rate of inflorescence%10.7a19.7ab50.0b*27.3b11.5a13.8a12.5a*41.3b
basal diametermm10.1b*7.6a*6.5a7.2a*>7.1c*5.9ab*6.1b5.4a*
living leafshoot−19.39.310.110.17.4a8.0ab9.7bc10.6c
LAI (leaf area index)m2m−22.41.511.
shoot dry massg shoot−138.1b*25.1a18.3a22.3a>26.3b*19.6ab17.5a18.0a
biomasskg m−
nutrient standing stock
Pg m−21.9c0.8b0.5a0.7ab1.
Ng m−236.7c15.6b9.6a16.1ab28.822.420.913.2
Cg m−21070.0b696.7b382.8a648.1ab1058.61017700.1470.6
Sg m−25.2c*2.1b1.2a3.5abc>1.2*
living rhizomekg m−
living rootkg m−20.3ab0.3b0.2a0.
senescence rhizomekg m−20.7*1.41.0*1.4b*1.4b0.3a*0.2a
senescence rootkg m−
bud biomasskg m−

3.2.2 Below-ground

The majority of the living rhizomes and roots is situated in the upper sediment layer, above the first horizontal rhizome (20–30 cm). Differences between the investigated years were found in the course of the analysis of variance of the different parameters (Table 2). Living rhizome biomass varied between 2.5–4.3 and 1.1–1.8 kg m−2 at sites PR and PV respectively, while the total living biomass (rhizome, root, buds) varied between 3.1–4.9 kg m−2 at the PR sites, and between 1.4 and 2.1 kg m−2 at the PV sites.

In the rhizome system, the senescence dry mass varied between 0.5 and 1.9 kg m−2 at the investigated sites. The senescence rhizomes made up less than half of the total below-ground biomass.

3.2.3 Nutrient concentrations and standing stock

The highest C concentrations were measured in the living roots and senescence rhizome and the lowest ones in living leaves (Table 3), while higher P and N concentrations were measured in the inflorescence, leaves, the living rhizome, and roots than in the stem and senescence rhizome. The S concentration was relatively variable in the examined reed organs. The living and senescence roots and living rhizomes contain the highest S concentration in contrast to above-ground organs.

Table 3

The element concentration (mg g−1) of the above-ground and below-ground reed organs (Legend see Table 2).

P. australisNutrient20032004
living leafP10.8*1.1*1<1.3ab1.3ab*1.4b*1.2a
living rhizomeP0.9*0.60.7*0.6a*0.6a1.0b*0.7a
senescence rhizomeP0.30.2*0.2<0.30.3*0.50.3
living rootP0.8b0.6a0.6a0.
senescence rootP0.8b0.6a*0.7b<0.81.1*0.90.8

In 2004 the stem contained more P and C, the living leaves and senescence roots more N, P and C, the living roots more N and C, the senescence rhizomes more P andN and the living rhizomes more C than in 2003. In the above- ground organs there were no remarkable differences in the S concentrations, nevertheless in the living leaf and the underground organs the S concentration was higher in 2003.

The C standing stock of living shoots varied between 383 and 1070 gm−2 in 2003 and between 471 and 1059 gm−2 in 2004. Similar to the case of the above-ground biomass, the C standing stock was also lower at PV sites than at the PR sites (Tables 2 and 4). As for the below-ground measurements, differences were found between the C standing stock of living (440–2.210 g m−2) and senescence biomass (220–912 g m−2).

Table 4

The element standing stock (g m−2) of the below-ground reed organs (Legend see Table 2).

P. australisNutrient20032004
living rhizomeP2.
senescence rhizomeP0.20.2*0.2<0.4bc0.5c*0.2ab0.1a
living rootP0.
senescence rootP0.

The N standing stock of the annual above-ground phytomass varied between 9.6 and 36.7 g m−2 in 2003 and between 13.2 and 28.8 g m−2 in 2004. It was higher in both years at the PR sites. Considerable differences were recorded at the individual sites between the N standing stock of the below-ground living and decaying biomass (4–63 g m−2 and 3.9–15.4 g m−2, respectively).

As for the P standing stock of the new shoots, it varied between the PR and PV sites (Table 2). It was similar in 2004 (0.7–1.9 gm−2) and in 2003 (0.5–1.9 gm−2). There was 0.8–3.0 and 0.3–1.2 g m−2 P standing stock in the below-ground organs (living and decaying biomass, separately).

As for the above- and the below-ground organs, their samples were separated by cluster analysis into two groups, before and after the introduction of pre-treated wastewater (Figs. 3a and 3b respectively). These groupings were confirmed by the results of the linear discriminant analysis as in the case of the analysis of the water samples.

Figure 3 Grouping of the sampling sites based on the results of the cluster analysis conducted on the above-ground (a) and below-ground reed samples (b).
Figure 3

Grouping of the sampling sites based on the results of the cluster analysis conducted on the above-ground (a) and below-ground reed samples (b).

4 Discussion

Nutrients in water and in sediment interstitial water

Simultaneous physical, chemical, and biological processes affected the pre-treated wastewater as it flowed through the reed parcel, towards the Virágosmajor-Canal. Nutrient removal by P. australis is achieved via two major processes: the absorption by the plant itself and microbial activity in the rhizosphere [16, 60, 61]. As a result, the concentration of the nutrients noticeably decreased between the inlet and the outlet. It was shown that the nutrient removal efficiency of the reed beds depends on (i) the loading rates, (ii) the distribution/ spread of the water at the inlet, (iii) the nutrient species occurring (organic or inorganic N and P forms), and (iv) the abiotic environmental conditions [62].

At the sites in the vicinity of the inlet (PR), the PO43 concentration of the water was significantly higher than at the sites close to the outlet (PV). This happened presumably via the phosphate uptake of the reeds, bacteria, and algae and through phosphate absorption, complexation, and precipitation with metals and clay particles [60, 63].

The decrease in nitrate concentration in the water from the inlet to the outlet suggested that N removal was taking place in the system. The major N removal mechanisms in the reed bed systems used for wastewater treatment are the combined nitrification and denitrification processes [16, 60]. The internal oxygen transport of P. australis to its rhizosphere makes the simultaneous existence of oxidised, anoxic, and reduced zones possible, as also the co-occurrence of these processes [16, 64, 65]. Aerobic bacteria are responsible for oxidizing ammonia, heavy metals and other toxic compounds. Nitrogen can also be taken up by the reeds and incorporated into its own biomass [60]. Ammonia was not detectable in the surface and interstitial water during the investigation period, which may be explained by the efficient ammonium-N up take of reeds and algae and by its removal through the nitrification-denitrification processes.

The phosphate and nitrate concentrations measured in the surface and interstitial water in this study were similar to those of e.g. [66], which were obtained in different reed stands of Lake Fertő/Neusiedler See; however, both of these were higher than the results of [67]. The reason for this may lie in the great fall in the water level of the whole lake and the partial drying out of the parcel on the shoreline in 2003.

The concentration of DOC in the surface water was higher at the outlet than at the inlet. However, in the reed stand even these values were exceeded. This might have been due to the contribution of the DOC originating from the decaying plant material accumulating in the sediment, as reflected by the high organic matter content of the sediment. [68], studying the nutrient release from integrated constructed wetland sediments, also found that sediments released substantially more organic matter than the incoming organic matter that could be degraded. The living reed roots and rhizomes are also capable of emitting DOC to a greater degree during the active growth period [69]. The sulphate concentration of water in most cases decreased as a function of the sediment depth and displayed an inverse relationship to the vertical profile of sulphide concentrations. This happened because sulphate reduction processes can also precipitate heavy metals, where they occur, in the root-zones [16, 70].

The water slowly flowing through the outlet, leaves the parcel, eventually reaching the Virágosmajor Canal and subsequently Fertőrákos Bay. As a result of the dilution of the outflow water its Eh, pH and oxygen content is lower, while its nutrient content (PO43,NO3) is higher than that of the open water at the reference site, showing a similar situation to the 6 km-wide reed stand which lies across the Bozi Canal [71].

The inhomogeneity of the reed stand area set aside for the subsequent cleaning of the pre-treated wastewater made it difficult to work out a sampling strategy which would allow the preparation of a mass balance with appropriate accuracy.

At the data evaluation stage, another circumstance urged caution, namely the quantity of wastewater discharged daily into the reed parcel. It did not usually inundate the parcel, partly due to the current water level fluctuation of Lake Fertő/Neusiedler See. Consequently, the pretreated wastewater did not leave the parcel at all times as surface water, i.e. the flow of the wastewater was not continuous. Therefore the determination of the residence time became highly difficult and in most cases impossible.


There are three main processes which interact and simultaneously influence each other: (i) nutrient uptake by Phragmites from both the sediment and the water, (ii) the mineral concentrations of the surrounding water and soil, and (iii) the biomass [12]. According to previous studies, nutrient enrichment of the water raises the N and P concentrations in the reed tissue [5, 79, 61, 72, 73]. In the present study, the N and P concentrations of the leaves were higher and the P concentration of the stem, senescence rhizomes and roots was also significantly higher after the introduction of pre-treated wastewater (in 2004) than before (2003). This indicated that the P. australis took up and stored more P from its nutrient-enriched environment (sediment interstitial water and surface water), which was caused by the introduction of the pre-treated wastewater. This is seems to be supported by the decrease in the N and P concentrations of the water from the inlet to the outlet. This result is in harmony with work of Meuleman et al. [20]. They also found that nutrient concentration and nutrient storage in P. australis vegetation in an infiltration wetland used for wastewater treatment was significantly higher than in the natural wetland. This statement is in harmony with the results of the study [74]. In previous studies, it has been pointed out that the uptake of nutrients by reeds is only of quantitative importance in low loaded systems [11].

The significant influence of the wastewater on the biometrical parameters of the reed was described by [13]. They found that the maximal density of shoots was higher, the biomass was twice as high and the shoot diameter was significantly greater in the treated water than in the natural reed stand. Ahigher shoot density was also recorded in our study, which, on the one hand can be attributed to the effect of wastewater, while on the reeds, and on the other hand to the result of the reeds’ winter harvesting.

All the above findings were reflected in and in harmony with the multivariate results as well, specifically that the reed stand functioned as an effective filtering area for the pre-treated wastewater. It retained it, thus leaving the processes of the open water unaffected. This is why its samples remained in one group while those in the parcel separated in both time (before and after) and space (inlet and outlet) after the introduction of the pre-treated wastewater. This situation is similar to the one witnessed in a constructed wetland, the Kis-Balaton Water Protection System [7577].

5 Conclusions

The nutrient concentrations of the reed organs, surface water, sediment interstitial water and the effect of pretreated wastewater were analysed in a partly separated reed stand used for pre-treated wastewater treatment, near Fertőrákos Bay at the landward edge of a wide reed belt of Lake Fertő/Neusiedler See.

Based on the results it can be concluded that the Phragmites australis stand, is (i) able to fulfil an important role in efficiently and sustainably removing nutrients (NO3,PO43, DN) from pre-treated wastewater, while (ii) conserving and protecting the natural processes of the open water. The investigation underlines the efficient applicability of “close-to-natural” reed parcels for the purification of pre-treated wastewater in a highly unstable hydrological system.


The authors are thankful to János Török (co-worker of Hydrometeorological Station at Fertőrákos) for his help in the fieldwork and to Gábor Horváth for the laboratory work. We would like to also thank Paul Thatcher for his work on our English versions. This work was supported by NKFP 3B/0014/2002 project.


[1] Hazelton E.L.G., Mozdzer T.J., Burdick D.M., Kettenring K.M., Whigham D.F., Phragmites australis management in the United States: 40 year of methods and outcomes. AoB plants 6, 2014, plu001; DOI:10.1093/aobpla/pl00110.1093/aobpla/pl001Search in Google Scholar

[2] Björk S., Ecologic investigations of Phragmites communis, Folia Limn. Scand. 1967, 14, 1–248Search in Google Scholar

[3] Květ J., Growth analysis approach to the production ecology of reed swamp plant communities, Hydrobiologia Bucuresti 1971, 12, 15–40Search in Google Scholar

[4] Haslam S., Some aspects of the life history and autecology of Phragmites communis Trin. A review, Pol. Arch. Hydrobiol. 1973, 20, 79–100Search in Google Scholar

[5] Dykyjová D., Hradecká D., Production ecology of Phragmites communis relations of two ecotypes to the microclimate and nutrient conditions of habitat, Folia Geobot. Phytotax. 1976, 11, 23–6110.1007/BF02853314Search in Google Scholar

[6] Raghi-Atri F., Bornkamm R., Wachstum und chemische Zusammensetzung von Schilf (Phragmites australis) in Abhängigkeit von der Gewässereutrophierung, Arch. Hydrobiol. 1979, 85, 192–228Search in Google Scholar

[7] Ho Y.B., Mineral composition of Phragmites australis in Scottish Lochs as related to eutrophication. I. Seasonal changes in organs, Hydrobiologia 1981, 85, 227–23710.1007/BF00017612Search in Google Scholar

[8] Hocking P.J., Seasonal dynamics of production, and nutrient accumulation and cycling by Phragmites australis (Cav.) Trin. ex Steudel in a nutrient-enriched swamp in inland Australia. I. Whole plants, Austr. J. Mar. Freshwater Res. 1989, 40, 421–44410.1071/MF9890421Search in Google Scholar

[9] Hocking P.J., Seasonal Dynamics of Production, and Nutrient Accumulation and Cycling by Phragmites australis (Cav.) Trin. ex Steudel in a Nutrient-enriched Swamp in Inland Australia. II. Individual Shoot, Aust. J.Mar. Freshwater Res. 1989, 40, 445–46410.1071/MF9890445Search in Google Scholar

[10] Ksenofontova T., Morphology, Production and Mineral Contents in Phragmites australis in Different Waterbodies of the Estonian SSR, Folia Geobot. Phytotax. 1988, 23, 17–4310.1007/BF02853293Search in Google Scholar

[11] Brix H., Functions ofmacrophytes in constructed wetlands, Wat. Sci. Tech. 1994, 29, 471–47810.2166/wst.1994.0160Search in Google Scholar

[12] Dykyjová D., Úlehlová B., Mineral economy and cycling of minerals in wetlands. In: Westlake D.F., Květ J., Szczepańsky A. (Eds.), The production ecology of wetlands. Cambridge University Press, New York, 1998, 78–16810.1017/CBO9780511549687.008Search in Google Scholar

[13] Hardej M., Ozimek T., The effect of sewage flooding on growth and morphometric parameters of P. australis (Cav.) Trin.ex Steudel, Ecol. Eng. 2002, 18, 343–35010.1016/S0925-8574(01)00095-7Search in Google Scholar

[14] Armstrong J., ArmstrongW., Phragmites australis–A pleliminary study of soil oxidizing sites and internal gas transport pathway, New Phytol. 1988, 108, 373–38210.1111/j.1469-8137.1988.tb04177.xSearch in Google Scholar

[15] Tulbure M.G., Ghjioca-Robrecht D.M., Johnston C.A., Whigham D.F., Inventory and ventilation eflciency of nonnative and native Phragmites australis (Common Reed) in tidal wetlands of the Chesapeake Bay, Estuaries and Coasts 2012, 35(5), 1353–135910.1007/s12237-012-9529-4Search in Google Scholar

[16] Brix H., Treatment of wastewater in the rhizosphere of wetland plants–The root zone method, Wat. Sci. Tech. 1987, 19, 107–11810.2166/wst.1987.0193Search in Google Scholar

[17] Schierup H.H., Brix H., Danish experience with emergent hydrophyte treatment systems and prospects in the light of future requirements on outlet water quality, Wat. Sci. Tech. 1990, 22, 65–7210.2166/wst.1990.0184Search in Google Scholar

[18] Brix H., Treatment wetlands: an overview. Constructed wetlands for wastewater treatment, Proc. Sci-tech. Conf. Gdańsk, Poland, 1995, 166–176Search in Google Scholar

[19] Kern J., Idler C., Treatment of domestic and agricultural wastewater by reed systems, Ecol. Eng. 1999, 12, 13–2510.1016/S0925-8574(98)00051-2Search in Google Scholar

[20] Meuleman A.F.M., Beekman J.H.P., Verhoeven J.T.A., Nutrient retention and nutrient use eflciency in Phragmites australis stands after wastewater application, Wetlands 2002, 22/4, 712–72110.1672/0277-5212(2002)022[0712:NRANUE]2.0.CO;2Search in Google Scholar

[21] Vymazal K., Greenway M., Tonderski K., Brix H., Mander Ü., Constructed Wetlands for wastewater treatment. In: Verhoeven J.T.A., Beltman B., Bobbink R., Whigham D.F. (Eds.), Wetlands and Natural Resource Management, Springer Verlage, Berlin, 2006, 69–9610.1007/978-3-540-33187-2_5Search in Google Scholar

[22] Ostendorp W., ’Die-back’ of reeds in Europe-a critical review of literature, Aquat. Bot. 1989, 5, 5–2610.1016/0304-3770(89)90063-6Search in Google Scholar

[23] Van der Putten W.H., Die–back of Phragmites australis in European wetlands: an overview of the European Research Programme on reed die-back and progression (1993–1994), Aquat. Bot. 1997, 59, 263–27510.1016/S0304-3770(97)00060-0Search in Google Scholar

[24] Vymazal J., Brix H., Cooper P.F., Green M.B., Haberl R., Constructed Wetlands for Wastewater Treatment in Europe. Backhuys Publishers, Leiden, 1998Search in Google Scholar

[25] Květ J., Dusek J., Husák S., Vascular plants suitable for wastewater treatment in temperate zones. In: Vymazal J. (Ed.), Nutrient Cycling and Retention in Natural and Constructed Wetlands. Backhuys Publisher, Leiden, 1999, 101–110Search in Google Scholar

[26] Németh N., Németh T., Comparison of natural and constructed reed habitats, Commun. Soil Sci. Plan. 2006, 37(15–20), 2553–256310.1080/00103620600823000Search in Google Scholar

[27] Sieghardt H., Maier R., Produktionsbiologische Untersuchun-gen and Phragmites-Beständen im geschlossenen Sclhilfgürtel des Neusiedler Sees. In: AGN Forschungsbericht 1981–1984. Burgenland, 1984, 191–221Search in Google Scholar

[28] Sieghardt H., Hammer L., Teuschl G., Das Schilfrohr (Phragmites australis (Cav.) Trin. ex Steudel), Wachstum und Production in verschiedenen Zonen des Schilfgürtels am Neusiedler See. 1984, BFB-Bericht 51, 37–47Search in Google Scholar

[29] Hietz P., Decomposition and nutrient dynamics of reed (Phragmites australis (Cav.) Trin ex. Steud) litter in Lake Neusiedl, Austria, Aquat. Bot. 1992, 43, 211–23010.1016/0304-3770(92)90068-TSearch in Google Scholar

[30] Dinka M., Szabó E., Tóth I., Changes in nutrient and fibre content of decomposing Phragmites australis litter, Internat. Rev. Hydrobiol. 2004, 89, 519–53510.1002/iroh.200410772Search in Google Scholar

[31] Ágoston-Szabó E., Dinka M., Némedi L., Horváth G., Decomposition of Phragmites australis rhizome in a shallow lake, Aquat. Bot. 2006, 85, 309–31610.1016/j.aquabot.2006.06.005Search in Google Scholar

[32] Löffler H., Neusiedler See: The limnology of a shallow lake in Central Europe. Dr. W. Junk B.V. Publishers, Hague, 197910.1007/978-94-009-9168-2Search in Google Scholar

[33] Kucsara M., Gribovszki Z., Kalicz P., Kiss K.A., A „Fertőrákosi biológiai szűrőmező” próbaüzemének vizsgálata [Investigation of the biological filter field at Fertőrákos during the trial run operation]. Nyugat-Magyarországi Egyetem, Erdőmérnöki Kar, Erdőfeltárási és Vízgazdálkodási Tanszék, Sopron, 2005 (in Hungarian)Search in Google Scholar

[34] Golterman H.L., Clymo R.S., Ohnstand M.A.M., Methods for physical and chemical analysis of freshwaters–1BP Handbook No. 8. Blackwell Sci. Publ. 1978, 1–213Search in Google Scholar

[35] Dykyjová D., Hejny S., Květ J., Proposal for international comparative investigations of productions by stands of reed (Phragmites communis), Folia Geobot. Phytotax. 1973, 8, 435–44210.1007/BF02852069Search in Google Scholar

[36] Schierup H.H., Biomass and primary production in Phragmites communis Trin. swamp in North Jutland, Denmark, Int. Ver. Theor. Angew. Limnol. Verh. 1978, 20, 94–9910.1080/03680770.1977.11896491Search in Google Scholar

[37] Fiala K., Growth and production of underground organs of Typha angustifolia L., Typha latifolia L. and Phragmites communis Trin. Pol. Arch. Hydrobiol. 1973, 20, 59–66Search in Google Scholar

[38] Krueger K.L., Ungar P.S., Anterior dental microwear texture analysis of the Krapina Neandertals, Cent. Eur. J. Geosci. 2012, 4(4), 651–66210.2478/s13533-012-0111-1Search in Google Scholar

[39] Vince T., Szabó G., Csoma Z., Sándor G., Szabó S., The spatial distribution pattern of heavy metal concentrations in urban soils-a study of anthropogenic effects in Berehove, Ukraine. Cent. Eur. J. Geosci. 2014, 6(3), 330–34310.2478/s13533-012-0179-7Search in Google Scholar

[40] Ward Jr J.H., Hierarchical grouping to optimize an objective function, J. Am. Stat. Assoc. 1963, 58(301), 236–24410.1080/01621459.1963.10500845Search in Google Scholar

[41] Day W.H., Edelsbrunner H., Eflcient algorithms for agglomerative hierarchical clustering methods, J Classif. 1984, 1, 7–2410.1007/BF01890115Search in Google Scholar

[42] Kovács J., Nagy M., Czauner B., Kovács I.S., Borsodi A.K., Hatvani I.G., Delimiting sub-areas in water bodies using multivariate data analysis on the example of Lake Balaton (W Hungary), J. Environ. Manage. 2012, 110, 151–15810.1016/j.jenvman.2012.06.002Search in Google Scholar PubMed

[43] Kovács J., Kovács S., Hatvani I.G.,Magyar N., Tanos P., Korponai J. et al., Spatial optimization of monitoring networks on the examples of a river, a lake-wetland system and a sub-surface water system, Water Resour. Manag., 2015, 29, 5275–529410.1007/s11269-015-1117-5Search in Google Scholar

[44] Magyar N., Hatvani I.G., Kovácsné Sz.I., Herzig A., Dinka M., Kovács J., Application of multivariate statistical methods in determining spatial changes in water quality in the Austrian part of Neusiedler See, Ecol. Eng. 2013, 55, 82–9210.1016/j.ecoleng.2013.02.005Search in Google Scholar

[45] Magyar N., Trásy B., Dinka M., Kutrucz Gy., Delineating water bodies on the Hungarian side of Lake Ferto/Neusiedler See. In: Geiger J., Molnár E.P., Malvić T. (Eds.), Theories and Applications in Geomathematics. Selected Studies of the 2012 Croatian-Hungarian Geomathematical Convent, Opatija. GeoLitera Publishing House, Szeged, 2013, 103–115Search in Google Scholar

[46] Hatvani I.G., Magyar N., Zessner M., Kovács J., Blaschke A.P., The Water Framework Directive: Can more information be extracted from groundwater data? A case study of Seewinkel, Burgenland, eastern Austria, Hydrogeol. J. 2014, 22(4), 779–79410.1007/s10040-013-1093-xSearch in Google Scholar

[47] Bradák B., Thamó-Bozsó E., Kovács J., Márton E., Csillag G., Horváth E., Characteristics of Pleistocene climate cycles identified in Cérna Valley loess–paleosol section (Vértesacsa, Hungary), Quaternary International 2011, 234(1), 86–9710.1016/j.quaint.2010.05.002Search in Google Scholar

[48] Bíró L., Polgári M., Tóth T., Vigh T., Refinement of genetic and structural models of the Úrkút manganese ore deposit (W-Hungary, Europe) using statistical evaluation of archive data, Cent. Eur. J. Geosci. 2012, 4(3), 478–49410.2478/s13533-011-0079-2Search in Google Scholar

[49] Kohlmayer N., Grasemann B., Ambiguous versus non ambiguous characterization of components from pictures of cataclasites, Cent. Eur. J. Geosci. 2013, 5(4), 576–58710.2478/s13533-012-0156-1Search in Google Scholar

[50] Kovács J., Bodnár N., Török Á., The application of multivariate data analysis in the reinterpretation of engineering geological parameters; Miocene sediments, cored along the metro line in Budapest, Cent. Eur. J. Geosci., in press, DOI:10.1515/geo-2016- 000510.1515/geo-2016-0005Search in Google Scholar

[51] Torghabeh A., Rezaee R., Moussavi-Harami R., Pradhan B., Kamali M., Kadkhodaie-Ilkhchi A., Electrofacies in gas shale from well log data via cluster analysis: A case study of the Perth Basin, Western Australia, Cent. Eur. J. Geosci. 2014, 6(3), 393–40210.2478/s13533-012-0177-9Search in Google Scholar

[52] Novák M., Kirchkeszner Cs., Palya D., Bodai Zs., Nyiri Z., Magyar N., et al., Optimization of Solid State Urea Clathrate Formation as a Chemical Separation Method Coupled to Compound Specific Stable Carbon Isotope Analysis, Int. J. Environ. An. Ch., 95(15), 1471–148810.1080/03067319.2015.1101644Search in Google Scholar

[53] Williams F., Holmes N., Dental microwear texture analysis of late Pliocene Procynocephalus subhimalayanus (Primates: Cercopithecidae) of the Upper Siwaliks, India. Cent. Eur. J. Geosci. 2012, 4(3), 425–43810.2478/s13533-011-0076-5Search in Google Scholar

[54] Duda R.O., Hart P.E., Stork D.G., Pattern classification. John Wiley & Sons, 2012Search in Google Scholar

[55] Kovács J., Kovács S., Magyar N., Tanos P., Hatvani I.G., Anda A., Classification into homogeneous groups using combined cluster and discriminant analysis, Environ. Modell. Softw. 2014, 57, 52–5910.1016/j.envsoft.2014.01.010Search in Google Scholar

[56] Simeonov V., Stratis J.A., Samara C., Zachariadis G., Voutsa D., Anthemidis A. et al., Assessment of the surface water quality in Northern Greece, Water Res. 2003, 37(17), 4119–412410.1016/S0043-1354(03)00398-1Search in Google Scholar

[57] Kovács J., Tanos P., Korponai J., Kovácsné Sz.I., Gondár K., Gondár-Sőregi K. et al., Analysis of water quality data for scientists. In: Kostas V., Dimitra V. (Eds.), Water quality and water pollution: evaluation of water quality data. InTech Open Access Publisher, Rijeka, 2012, 65–94Search in Google Scholar

[58] Juahir H., Zain S.M., Yusoff M.K., Hanidza T.T., Armi A.M., Toriman M.E. et al., Spatial water quality assessment of Langat River Basin (Malaysia) using environmetric techniques, Environ. Monit. Assess. 2011, 173, 625–64110.1007/s10661-010-1411-xSearch in Google Scholar

[59] Riebe D.J., Niziolek L.C., Investigating Compositional Variation of Ceramic Materials during the Late Neolithic on the Great Hungarian Plain–Preliminary LA-ICP-MS Results, Cent. Eur. J. Geosci. 2015, 7(1), 423–44510.1515/geo-2015-0044Search in Google Scholar

[60] Brix H., Use of constructed wetland in water pollution control: historical development, present status, and future perspectives, Wat. Sci. Tech. 1994, 30, 209–22310.2166/wst.1994.0413Search in Google Scholar

[61] Ruiz J., Velasco M., Nutrient bioaccumulation in Phragmites australis: Management Tool for Reduction of Pollution in the Mar Menor, Water Air Soil Poll. 2010, 205, 173–18510.1007/s11270-009-0064-2Search in Google Scholar

[62] Meuleman A.F.M., van Logtestijn R., Rijs G.B.J., Verhoeven J.T.A., Water and mass budgets of a vertical-flow constructed wetland used for wastewater treatment, Ecol. Eng. 2003, 20, 31–4410.1016/S0925-8574(03)00002-8Search in Google Scholar

[63] de Bashan L.E., Bashan Y., Recent advances in removing phosphorus from wastewater and its future use as a fertilizer (1997–2003), Wat. Res. 2004, 38, 4222–424610.1016/j.watres.2004.07.014Search in Google Scholar PubMed

[64] Brix H., Do macrophytes play a role in constructed treatment wetlands? Wat. Sci. Tech. 1997, 35/5, 11–1710.2166/wst.1997.0154Search in Google Scholar

[65] Vretare V., Internal oxygen transport to below-ground parts: importance for emergent macrophytes. PhD thesis, Lund University, Sweden, 2001Search in Google Scholar

[66] Ágoston-Szabó E., Dinka M., Changes in sediment and sediment interstitial water characteristics in Lake Fertő/ Neusiedler See, Opuscula Zoologica 2004, 35, 3–17Search in Google Scholar

[67] Ágoston-Szabó E., Dinka M., Chemical Properties of the Sediment Interstitial Water at Lake Fertő/Neusiedler See. In: Vymazal J. (Ed.), Water and Nutrient Management in Natural and Constructed Wetlands, Springer, Dordrecht, 2010, 237–25010.1007/978-90-481-9585-5_17Search in Google Scholar

[68] Dong Y., Kayranli B., Scholz M., Harrington R., Nutrient release from integrated constructed wetlands sediment receiving farmyard run-off and domestic wastewater, Water Environ. J. 2013, 27: 439–45210.1111/j.1747-6593.2012.00361.xSearch in Google Scholar

[69] Pinney L.M., Westerhof P.K., Baker L., Transformations in dissolved organic carbon through constructed wetlands, Wat. Res. 2000, 34, 1897–191110.1016/S0043-1354(99)00330-9Search in Google Scholar

[70] Ággoston-Szabó E., Chemical characteristics of sediment interstitial water at Lake Fertő/Neusiedler See, Ecohydrol. Hydrobiol. 2004, 4, 67–76Search in Google Scholar

[71] Dinka M., Ágoston-Szabó E., Berczik Á., Kutrucz Gy., Influence of water level fluctuation on the spatial dynamic of the water chemistry at Lake Fertő/Neusiedler See, Limnologica 2004, 34, 48–5610.1016/S0075-9511(04)80021-5Search in Google Scholar

[72] Květ J., Mineral nutrients in shoots of reed (Phragmites communis Trin.), Pol. Arch. Hydrobiol. 1973, 20, 137–147Search in Google Scholar

[73] Dinka M., The effect of mineral nutrient enrichment of Lake Balaton on the common reed (Phragmites australis), Folia Geobot. Phytotax. 1986, 21, 65–8410.1007/BF02853399Search in Google Scholar

[74] Dinka M., Ággoston-Szabó E., Szeglet P., Comparison between biomass and C, N, P, S contents of vigorous and die-back reed stands of Lake Fertő/ Neusiedler See, Biologia 2010, 65, 237–24710.2478/s11756-010-0006-xSearch in Google Scholar

[75] Hatvani I.G., Kovács J., Székely Kovács I., Jakusch P., Korponai J., Analysis of long-term water quality changes in the Kis-Balaton Water Protection System with time series-, cluster analysis and Wilks’ lambda distribution, Ecol. Eng. 2011, 37, 629–63510.1016/j.ecoleng.2010.12.028Search in Google Scholar

[76] Hatvani I.G., Clement A., Kovács J., Székely Kovács I., Korponai J., Assessing water-quality data: the relationship between the water quality amelioration of Lake Balaton and the construction of its mitigation wetland, J. Great Lakes Res. 2014, 40, 115–12510.1016/j.jglr.2013.12.010Search in Google Scholar

[77] Kovács J., Hatvani I.G., Korponai J., Székely Kovács I., Morlet wavelet and autocorrelation analysis of long-term data series of the Kis-Balaton water protection system (KBWPS), Ecol. Eng. 2010, 36, 1469–147710.1016/j.ecoleng.2010.06.028Search in Google Scholar


Figure A1 Plot of the first and second principal component scores obtained from the interstitial water samples.
Figure A1

Plot of the first and second principal component scores obtained from the interstitial water samples.

Received: 2015-8-19
Accepted: 2015-10-25
Published Online: 2016-1-3
Published in Print: 2016-1-1

© M. Dinka et al., published by De Gruyter Open.

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.

Downloaded on 8.12.2023 from
Scroll to top button