1 Introduction

Due to the lack of water in the arid and semiarid areas in the coastal regions, potable and fresh water luxurious uses are increasingly contested. In order to solve this problem, groundwater gained an increasing role in the planning and in the development of additional water supplies (Hamdy 2002). Its use for irrigation purposes has reduced the crop yields (Ben Asher et al. 2002; Ayars et al. 2006), increased the soil salination and deteriorated the soil structure and permeability (Shainberg et al. 2002). Yield reduction occurs when the salts accumulate in the root zone to such an extent that the crop is no longer able to extract sufficient water from the salty soil solution, resulting in a water stress for a significant period of time (Ayers and Westcot 1985). The equation relating relative yield (Yr, dimensionless 0–1) and the electrical conductivity of the irrigation water ECw (dS m−1) is the following (Maas and Hoffman 1977):

$$ \mathrm{Yr}=1-\mathrm{b}\left({\mathrm{EC}}_{\mathrm{w}}-\mathrm{a}\right) $$
(1)

where a is the crop tolerance represented by the threshold salinity value (dS m−1), b is the crop sensitivity represented by the slope (rate of yield reduction per unit of salinity), and ECw is the electrical conductivity of irrigation water. Crops produce at their full potential yield until salinity reaches a threshold value, which is the maximum average salt concentration in the root zone that does not reduce the yield; above the threshold value, the yield decreases as a linear function of salinity until the plants die when it drops sharply to zero.

Beside the fertilization capacity of the different nutrients, there is clear nutrient–salt interaction and plant growth (Silberbush et al. 1988). Hence, the use of enhanced nitrogen fertilization is a potential tool in precluding salt induced stress in many crops (Feigin 1985). In addition, potassium has also a high contribution in alleviating salinity problems (Wakeel 2013).

Experimental studies related to nutrient-salt interactions, and plant growth and quality were completed for several plants, including: horticultural crops (Prazeres et al. 2013; Saykhul et al. 2013), leguminous crops (Correia et al. 2010; Bulut et al. 2011) and grasses (Morshedi and Farahbakhish 2012). It was shown that yield response to salinity tolerance is not always a constant as the model of Maas and Hoffman (1977) suggests, but fertilization increases tolerance to salinity (salinity threshold value); however, at high salinity levels, the osmotic effect is more pronounced and enhanced fertilization may provoke a negative effect on crop yield (Beltrão et al. 2002). The general plant physiological mechanisms for these nutrient-salt interaction phenomena were explained by Grattan and Grieve (1999). Studies related to the combined effects of salinity and fertilizers, were also undertaken for grasses (Sagi et al. 1996, 1997).

The reuse of treated wastewater is considered as an alternative to the use of potable water in Mediterranean agriculture and landscape (Carr et al. 2011), especially in golf courses (Salgot et al. 2012). Hence, according to the development of the construction of golf courses in coastal areas, several studies were undertaken on the response of Bermuda grass (Cynodon dactylon, L. Pers) to the reuse of treated municipal wastewater. This response was generally expressed by the yield. It was observed that the minimal yield for the good visual appearance (GVA) of Bermuda grass was around 9,000 kg ha−1 month−1 (Beltrão et al. 1999). However, Costa (2004) observed that grass quality is not only directly related to grass yield, but also it depends on some other environmental factors involved, such as the climate, water and nutrients, and soil salinity, being the grass quality one of the best parameters to qualify the appearance of the lawns.

This study focused on the response of Bermuda grass irrigated with municipal treated wastewater in two different locations in southern Portugal (Algarve). Taking into consideration the chemical composition of this irrigation water, which includes some essential macronutrients and other ions, it is expected that yield can be ensured avoiding the excessive use of fertilizers.

2 Materials and Methods

2.1 Site and Turfgrass

The experiments were carried out in the Salgados golf course and on the Campus of Penha, University of Algarve, Portugal, from May up to the end of July.

The Bermuda turfgrass Cynodon dactilon, was selected for this study, single or consociated with other species, because it is the most frequently used cultivar in fairways and roughs of golf courses, and it represents their largest area. Cynodon dactilon, ‘Bermuda 419’ and ‘Savanna’ are the most commonly used species on the Mediterranean golf courses. “Savanna” is exclusively applied on the roughs; on the other hand, “Bermuda 419”, which is most used, is grown on tees and fairways (Beltrão et al. 2009). They are warm-season species, adapted to a wide range of soil conditions, even though Bermuda grass is generally not cold tolerant. Its intolerance to shade necessitates the use of alternative warm-season species on sites where trees and other structures restrict sunlight penetration (Beard 1982); these species are perennial, have a long life, quick growth at 18–35 ºC, but become brown when temperature decreases. Its growth is extremely vigorous compared to other species (Mateo Box and Urbano 1985).

2.2 Soil and Water Quality Parameters

The experiments were conducted in fluvisols, containing approximately 30 % of clay, similar to the soil types of the local fairways and roughs of most of the golf courses located in the region of Algarve. According to water sources and to the location of the experimental sites, there were three different water treatments, as follows: Potable water S0 (ECw = 0.3 dS m−1) on Campus da Penha, Wastewater WW1 S1 (ECw = 1.6 dS m−1) in Salgados Golf course, and Wastewater WW2 S2 (ECw = 2.4 dS m−1) on Campus da Penha. Table 1 shows the concentration of typical chemical parameters. The values of the electrical conductivity ECw were higher in wastewater, as expected. Potassium (K) concentration of the wastewater was very high, equivalent to a very high K mineral fertilization; on the other hand, nitrogen (N) concentration was low, probably due to prior biological treatment, where microorganisms consumed high quantities of N needed for their metabolism.

Table 1 Concentration of some parameters in irrigation water

2.3 Climate and Water Regimes

The climate of Algarve, particularly the south shore, can be considered as Mediterranean. According to Köppen (1936), it is classified as Csa, with semi-arid characteristics, identified by mild rainy winters and by warm and dry summers. Snow and frost are not common. The annual average temperature, in Faro, ranges between 12 ºC in January and 23 ºC in August, a difference of only 11 ºC.

Table 2 shows the most important climatic parameters - temperature, Penman potential evapotranspiration and precipitation (30 years average) of the coastal area of Algarve, where the study was carried out. The annual value for potential evapotranspiration was 1,280 mm, and the precipitation was around 500 mm. In addition, the values of potential evapotranspiration and precipitation show opposite distribution during the year. During the hottest semester, when the photosynthesis is more efficient and the water is more essential for plant growth, precipitation is less available; irrigation is the only solution to balance this situation (Costa et al. 2011). The average of water shortage, defined as the difference between the rates of evapotranspiration and precipitation, verified between March and October, may reach 940 mm per year. The daily water balance equation for the root zone is given by:

Table 2 Climatic parameters (1970–2000)
$$ \mathrm{I}+\mathrm{P}+\mathrm{CR}={\mathrm{ET}}_{\mathrm{a}}+{\mathrm{D}}_{\mathrm{r}}+\mathrm{R}+\varDelta \mathrm{S} $$
(2)

where I is the net irrigation, P is the precipitation, CR is the capillary rise from the groundwater table to the root zone, ETa is the actual seasonal average evapotranspiration, Dr is the drainage below the root zone, R is the runoff, and ΔS is the change in water storage within the root zone; the units for all these parameters are mm d−1. Capillary rise CR and runoff R were zero during the experimental period, natural precipitation P was negligible, and ∆S ~ 0 due to the high irrigation frequency (daily). Hence, Eq. (2) is simplified to:

$$ \mathrm{I}={\mathrm{ET}}_{\mathrm{a}}+{\mathrm{D}}_{\mathrm{r}} $$
(3)

During the experimental period, the net irrigation water rate I was applied daily and its average value was computed as follows:

  1. a)

    For I < 6 mm d−1, ETa = ETc.,adj, Dr = 0 (no drainage effects)

$$ \mathrm{I}={{\mathrm{ET}}_{\mathrm{c}.}}_{,\mathrm{adj}}={\mathrm{K}}_{\mathrm{s}}.{\mathrm{K}}_{\mathrm{c}}.{\mathrm{ET}}_0 $$
(4)
  1. b)

    For I = 6 mm d−1, ETa = ETc., Dr = 0 (no drainage effects)

$$ \mathrm{I}={\mathrm{ET}}_{\mathrm{c}.}={\mathrm{K}}_{\mathrm{c}}.{\mathrm{ET}}_0 $$
(5)
  1. c)

    For I > 6 mm, ETa = ETc., Dr > 0 (drainage effects)

$$ \mathrm{I}={\mathrm{ET}}_{\mathrm{c}.}+{\mathrm{D}}_{\mathrm{r}}={\mathrm{K}}_{\mathrm{c}}.{\mathrm{ET}}_0+{\mathrm{D}}_{\mathrm{r}} $$
(6)

where ETc.,adj is the adjusted evapotranspiration under non-standard conditions and its value is used when the crop is under stress conditions (Allen et al. 1998); Kc is the crop coefficient under the standard conditions (Doorenbos and Pruitt 1984), ET0 is the average reference evapotranspiration (Penman-Monteith), and Ks (dimensionless 0–1) describes the effect of crop stress on crop evapotranspiration ETc. The estimation of these crop coefficients and crop factors, useful to the definition of crop production systems and characterized by the link between the production and the water application, have been developed by Allen and Pereira (2009), Allen et al. (2011a, b) and Pereira et al. (2012).

In order to simulate different irrigation regimes, a sprinkler point source experimental design was used (Fig. 1). This design is based on the assumption that a point creates a linear irrigation gradient from the water point source obtained from a triangular sprinkler water diagram, producing a gradual change in water application, (Hartman and Zengerl 1985; Or and Hanks 1992). The wettest zone was near the sprinkler and the treatments near the borders were the driest. The irrigation system was stopped when wind speed was larger than 1 m s−1. A 323/90 Naan sprinkler was used, with 2.4 x 4.8 mm diameter nozzles, using a 250 kPa sprinkler pressure, with a wetting radius of about 10 m. Christiansen uniformity coefficient CUC (Christiansen 1942) was used to determine the uniformity of water distribution at 2, 4, 6, and 8 m from the sprinkle point source; the respective CUC values obtained were 94, 93, 92, and 80 %. The plots were irrigated once a day.

Fig. 1
figure 1

Point source experimental design applied in the Salgados golf course. Irrigation gradient was expressed by the daily net irrigation rates (I, mm d−1) and the distance from the sprinkler in m

The colour is one of the best indicators of the aspect-quality of turfgrass (Turgeon 1980). Hence, a colour visual method was used to define the turfgrass quality (Mantell and Stanhill 1981), associated to the sprinkler point source design. The appearance of the lawn was compared with the observation of the colour values of “Standard Soil Colour Charts” (Munsell 1991). The grass yields, under different irrigation quantities, were sampled monthly. The grass was mowed with scissors and collected at the predetermined distance from the sprinkler.

2.4 Statistical Analysis

The effects of treatments were evaluated using the analysis of variance (ANOVA), when no cause and effect relationship is known. The statistical test Dunnett T3 was selected, in order to identify the statistical difference among multiple mean values at the 95 % significance level. SPSS 11.0 (2002) was used. When problems of lack of randomization were known due to the point source experimental design, a geostatistical approach was undertaken, according to Panagopoulos et al. (2006).

3 Results and Discussion

Figure 2 shows the Bermuda turfgrass response to wastewater and potable water application, from May up to the end of July, for the three trials (potable water PWS0 - ECw = 0.3 dS m−1), wastewater WWS1 - ECw = 1.6 dS m−1), and wastewater WWS2 - ECw = 2.4 dS m−1). For low irrigation water quantities (I < 3 mm d−1), water was the limiting factor; the yield was sharply enhanced with increasing rate of water application. Grass yield response was lower on the WWS2 treatment near the dry level and higher on the PWS0 treatment, due to the higher osmotic pressure of the soil. For intermediate water quantities (3 mm d−1 < I < 6 mm d−1), wastewater application triggered slightly higher yields compared to potable water application, particularly if I is close to 6 mm d−1. This was probably due to the lower concentration of nutrients of the potable water when compared to wastewater (Table 1). For high water quantity application (I > 6 mm d−1), there were increased yields in the three treatments. The sensitivity of the yield to water was higher in the WWS1 treatment (sharp increase), due to the lower salinity level when compared to treatment WWS2, explained by soil leaching. On the other hand, the enhancement was lower for the PWS0 treatment (monotonic increase), due to its lower concentration of nutrients combined with soil leaching.

Fig. 2
figure 2

Comparative turfgrass responses to potable water (PWS0) and wastewater (WWS1 and WWS2). Irrigation is expressed by the net irrigation water rates I (mm d−1)

By comparing visually the same grass yield quantity coming from the different treatments, it is recognized that grass quality was not only related to the grass yield but also on several environmental factors involved on grass growth, such as climate, water, nutrients (particularly N and K) and salinity. Hence, generally, at values of I > 5.4 mm d−1, a good quality of the lawns was verified. But also in some cases, a good visual appearance of the grass was seen for lower values of I, i.e., 4 or 5 mm d−1, according to Costa (2004).

The presence of nutrients in wastewater can decrease or even avoid mineral fertilization. However, at low water quantity application, the grass yield response was lower with wastewater treatments, due to the higher osmotic pressure of the soil. On the other hand, for wetter regimes, when the salinity of wastewater reached higher values (ECw = 2.4 dS m−1), the rate of the yield increase was lower due to the detrimental effect of salinity. This effect was most probably attributed to the presence of Cl in WW2 treatment. In golf courses the costs of lawns maintenance (such as cuts frequency and fertilizers application) are high. Therefore, reducing turfgrass growth but simultaneously maintaining a good visual appearance, may contribute to reduce these costs. This can be achieved by using wastewater application providing that chemical and microbiological properties are not limiting.

4 Conclusions

A high research priority must be given on the response of grass quality to environmental factors (including irrigation and fertilization) in order to maintain the sustainability and the good appearance of golf courses in the Mediterranean area. Experimental results and discussion have confirmed that yield improves as net irrigation increases. Higher concentrations of macronutrients, like K, Ca and Mg, in the wastewater, either WWS1 or WWS2, enhanced turfgrass yield. The high K concentration of the wastewater may also improve the soil structure, with benefits to the soil hydraulic properties, which interfere with the transport of water to the roots, and therefore, with turfgrass development.