Impact Factor 3.678

The world's most-cited Plant Sciences journal

Original Research ARTICLE

Front. Plant Sci., 16 April 2018 |

Optimization of Nitrogen, Phosphorus, and Potassium Fertilization Rates for Overseeded Perennial Ryegrass Turf on Dormant Bermudagrass in a Transitional Climate

  • 1Key Laboratory of Horticultural Plant Biology, College of Horticulture and Forestry Science, Huazhong Agricultural University, Wuhan, China
  • 2College of Plant Science and Technology, Huazhong Agriculture University, Wuhan, China
  • 3Agricultural Department, The University of Swabi, Khyber Pakhtunkhwa, Pakistan
  • 4Southwest Engineering Technology and Research Center of Landscape Architecture (State Forestry Administration), Southwest Forestry University, Kunming, China

Bermudagrass [Cynodon dactylon (L.) Pers.] turf loss due to severe cold in transitional climates is a major concern. To overcome this problem, warm-season grass is often overseeded with a cool-season turfgrass. In this study, modeling and efficient nutrient management were used to evaluate this problem. A three-factor and five-level central composite rotatable design (CCRD) with a simulation of a regression model was used to optimize fertilization rates. The study investigated the combined effects of fertilization with nitrogen (N), phosphorus (P), and potassium (K) on both the morphological and physiological attributes and on the integrated turf performance (ITP) of overseeded perennial ryegrass (Lolium perenne). Fertilization with N and P significantly increased turf height, density, color, fresh and dry weights, while N, P, and K significantly affected turf cover, quality and winter-kill. The Spring transition was delayed by fertilization with N and P, and accelerated by fertilization with K. Photosynthesis (Pn), transpiration (Tr), and stomatal conductance (Gs) were considerably enhanced by fertilization with N, P, and K. Protein levels and total chlorophyll levels were substantially increased by fertilization with N and P and with N, P, and K, respectively, during a 2-year period. During two separate experiments conducted during 2 consecutive years, the optimal combinations of N, P, and K were N: 30, P: 24, K: 9, and N: 30, P: 27, K: 6 g m−2. The major conclusion of this study is that a balanced nutrient application utilizing N, P, and K is key to enhancing the winter performance of perennial ryegrass.


In transitional zone environments, it is difficult for both cool and warm-season grasses to perform well throughout the year (Jiang and Huang, 2000). In transitional zones, cool-season grasses are favored relative to warm-season grasses, because of the winter dormancy of warm-season grasses, which reduce turfgrass cover for an extended period of time. Bermudagrass is a warm-season turfgrass species used in sports fields, golf courses, parks, lawns, and prairies, because of its recuperative potential and extreme thermotolerance, however, low temperatures induce extended winter dormancy in bermudagrass (Fan et al., 2014). To overcome this problem, bermudagrass is commonly overseeded with a cool-season turfgrass (Shi et al., 2014a,b). The practice of winter overseeding usually involves seeding a cool-season grass into an established warm-season turfgrass for temporarily restoring turfgrass cover, color, quality, playability, aesthetics, functional quality, and overall performance during the dormant cool period (Pompeiano et al., 2011). Perennial ryegrass (Lolium perenne L.) is a major turf and forage grass, and commonly overseeded on bermudagrass (Cynodon dactylon L.) that is applied in the late summer through the mid-fall (Richardson, 2004; Verdenal et al., 2008; Yang et al., 2014). An ideal turfgrass for overseeding should quickly establish and maintain high turf quality. Additionally, an ideal turfgrass should not impede a smooth transition to bermudagrass in the spring (Beard and Menn, 1988).

An appropriate amount of nitrogen (N) fertilizer promotes photosynthesis, the production of dry matter, and stress tolerance (Guo et al., 2009). However, either under or over fertilization negatively affects plant growth (Hons et al., 2004). Deficiencies in phosphorus (P), another most important macronutrient negatively affects plant growth (Vance, 2001). P promotes winter survival, cold hardiness and dry matter production (Tyler et al., 1981). Similarly, potassium (K) fertilization plays a key role in turfgrass quality, root growth, disease resistance and cold hardiness (Snyder and Cisar, 2000).

Climatic changes (extreme low and high temperature) will cause instability and are anticipated to become more frequent. The low temperatures and short photoperiods of winter reduce grass production (Uleberg et al., 2014). These abiotic stresses adversely affect plant growth by inducing morphological, physiological, and biochemical changes in plants (Larcher, 2003). Inadequate fertilization enhances sensitivity to abiotic stress, even in highly managed turf surfaces (Pompeiano et al., 2011). Although the recommended nutrient rate for each crop is impractical for a number of reasons, and sessile nature of plants force agriculturists to develop more efficient fertilization programs for different seasons and localities.

Previous studies on turf grasses were limited to individual turfgrass species only in warm or cool seasons that tested the effects of nutrients in laboratories and nurseries. Unfortunately, studies on the quantitative relationship between fertilizers and turf quality are rarely reported, and limited research has been done on overseeding practices. There are few reports that investigated the combined effects of N, P, or K on cool-season turfgrasses, but to the best of our knowledge, no one has attempted to optimize the combination of N, P, and K on overseeded perennial ryegrass and quantify the interactive effects of these macronutrients on integrated turf performance (ITP).

Thus, the main objectives of this study were to quantify the interactive effects of N, P, and K fertilization on the morphology, physiology, and ITP of overseeded perennial ryegrass, and to develop a model that predicts the optimal combination of N, P, and K for maximum ITP.


Experimental Site Description and Sample Collection

The experiment was conducted from 2014 to 2016 at an ornamental plant garden at Huazhong Agricultural University (HZAU), Wuhan-China (30°28 N, 114°21 E). Meteorological data for both growing years were collected from a weather station located near the experimental site. This region has a humid subtropical monsoon climate, with a mean annual rainfall of 1,200–1,500 mm. From September of 2014 to May of 2015, the average annual temperature ranged from 5 to 22°C, while it was 4.6 to 21.6°C from September 2015 to May 2016 (Figure 1). There was snowfall in January and February during both years. Soil samples were taken randomly at a depth of 20 cm to test its physiochemical properties. The soil has a loess structure, pH 6.5, 19.80 g kg−1 organic matter content, total N of 1.39 g kg−1, available P of 25.56 mg kg−1, and available K of 151.14 mg kg−1.


Figure 1. Monthly total rainfall and monthly maximum, minimum and mean temperature during the growing seasons of the two experiments.

Field Preparation and Crop Husbandry

Both experiments were initiated in May of 2014 and 2015. The field selected for experimentation was sprayed with glyphosate, a pre-emergence herbicide, and well plowed. Plant remnants, debris, and weeds were removed. A thick layer of washed sand was broadcasted and thoroughly mixed with a rotavator and leveled. A total of 24 plots with dimensions of 2 × 1.5 m (length × width) were arranged in 4 rows. Each row contained 6 plots. There were 50 cm distance between each row. This arrangement created 6 columns that contained 4 plots each. There were 60 cm distance between each column. Hybrid bermudagrass [Cynodon dactylon (L.) Pers.] “Tifway419” was plugged in the second week of May from one square inch plugs at a distance of 5–6 inches. Upon new growth, plots were fertilized with 10 g m−2 of N (urea, 46%), P (superphosphate, P2O5 ≥ 15%), and K (potassium sulfate, K2O ≥ 50%). Management practices like weeding, insect pest examination and irrigation were performed to maintain high quality turf for overseeding.

Overseeding experiments were conducted in the first week of October 2014 and 2015. Bermudagrass was vertically mowed in two directions to a height of <1 cm and slightly raked with brushes to ensure that seeds dropped in-to the soil. Perennial ryegrass (Lolium perenne L.) seeds were obtained from the ornamental plant nursery at HZAU and were overseeded at a density of 30 g m−2. Plots were irrigated with an automatic sprinkler, immediately after overseeding. After germination, the plots were fertilized with combinations of N, P, and K (Table 2).

Urea (N ≥ 46%), superphosphate (P2O5 ≥ 15%), and potassium sulfate (K2O ≥ 50%) were used as fertilizers. An additional 20 g N m−2 was also applied in late April of each year to promote the recovery of the bermudagrass.

Data Collection and Measurements

Morphological Data

Percent ryegrass cover was measured three times biweekly using the line-intersect method (Jellicorse et al., 2012) in the second week of November during both experiments. When the percent cover reached >80%, we analyzed turf height, density, color, and quality. Turf height was measured four times on five randomly selected tillers. Turf density based on average leaf counts was measured thrice on a monthly basis by uprooting five randomly selected plants. A biweekly visual rating scale (1–9) for turf color and quality was used 9 times. Winter-kill was measured quantitatively on a rating of 0–9. A 10 cm2 grid was randomly placed on the turf at three different dates in January of each year. The shoots within the grid were cut above the crowns and fresh weights were determined and then normalized to square decimeters (g/dm2). The same procedure was opted for aboveground dry weight except that samples were oven dried at 65°C for 48 h. For both experiments, the spring transition of bermudagrass from overseeded ryegrass was measured three times every week beginning in the first week of May using the line-intersect method.


Photosynthetic parameters, i.e., photosynthesis (Pn), transpiration (Tr), stomatal conductance (gs), inter cellular CO2 concentration (Ci), and water use efficiency (WUE) were measured by LI-6400 portable photosynthesis system (LI-COR Inc., Lincoln, NE, USA). Measurements were taken by holding multiple intact leaves of same culm in the chamber. The constant airflow rate was 500 μmol s−1, saturation PFD was 1,200 μmol m−2 s−1, the CO2 concentration was 350 cm3 m−3, and the temperature was 20°C with 80% relative humidity. Measurements were made during sunny days between 09:00 to 11:30 a.m. (Beijing time) on four different dates at one-week intervals during January 2015 and 2016. WUE is the ratio of Pn to Tr.

Protein and Chlorophyll

Protein content in the fresh leaves of ryegrass was quantified using a spectrophotometer and a commercial assay kit (A045-2) from Nanjing Jiancheng Bioengineering Institute, China. Total chlorophyll content was quantified as previously described (Esmaeili et al., 2015).

Experimental Design and Statistical Analysis

A three-factor (N, P, and K) and five-level central composite rotatable design (CCRD) was used for optimization and evaluation of the combined effects of N, P, and K on the morphology, physiology of the ryegrass, and the interactive effects on ITP. The actual amount of fertilizer varied between a minimum value of 0 and a maximum value of 30 g m−2. Treatment levels and ranges of independent variables in the experiment were coded (Table 1) according to the following Equation (1):

Xi=xi-xi0Δxi    (1)

where Xi is the dimensionless coded value of the ith variable, xi is the real value of the ith variable, xi0 is the real value of the ith variable at the center point, and Δxi is the scaling factor (Montgomery, 2001). The independent variable limits and their values are shown in (Table 2).


Table 1. Independent variables, codes, and both their coded and actual values in the central composite rotatable design.


Table 2. Experimental design and layout of treatment combinations in the central composite rotatable design.

In the CCRD, the number of experimental units was based on a complete 2k factorial. The total number of experimental runs was calculated with Equation (2)

T=2k+2k+m0    (2)

where T is the number of experimental runs, k is the number of factors, and m0 is the number of replicates of the center point.

The center point of the CCRD needs to be replicated more than five times to accurately determine its value. In contrast, the other treatments do not need to be replicated (Montgomery, 2001). The center point treatment was replicated 9 times, with a control (runs = 23 + 2 × 3 + 9 + 1 = 24).

Morphological, physiological, and ITP responses to N, P, and K fertilization ratios were fit to a second-order polynomial Equation (3)

Y=a0+a1x1+a2x2+a3x3+a4x1x2+a5x1x3+a6x2x3       +a7x1x1+a8x2x2+a9x3x3    (3)

where Y is the dependent variable; a0 is a constant; a1, a2, and a3 are linear coefficients of x1, x2, and x3, respectively; a4, a5, and a6 are interaction coefficients; a7, a8, and a9 are the quadratic coefficients of x1, x2, and x3, respectively; and x1, x2, and x3 are the coded values for N, P, and K respectively.

Linear regression of the above equation was performed for each variable. The constant, regression coefficient of linear, quadratic, and interaction terms were determined in SPSS (16.0). Significance of the linear, quadratic, and interaction forms was judged using an F-value at P < 0.05. The adequacy of the regression equations was evaluated using the coefficient of determination, R2. MS Excel 2016 and SPSS (16.0) were used for statistical analysis and modeling. Sigma plot (8.0) was used to generate 3D surface plots to show interactive effects on ITP.


Percent Ryegrass Cover

Positive values for the parameter coefficients N (a1), P (a2), and K (a3) during the first experiment (2014–2015) and the second experiment (2015–2016) demonstrate the enhancing effect of these nutrients on ryegrass cover. The relative magnitude of absolute values was: a1 > a2 > a3 (Table 4). The lowest percent cover was observed in the control plots, 45.33 and 51.33% during both experiments (Table 3). In the first and second experiments, the highest percent cover was 84.33 and 80.33% (Table 3), which was observed with N, P, K = 23.918 and N, P, K = 15, respectively.


Table 3. Morphological performance of overseeded perennial ryegrass on bermudagrass in response to different treatment combinations of nitrogen, phosphorus, and potassium fertilizer in the experiments conducted from 2014 to 2015 and 2015 to 2016.


Table 4. Parameter coefficients of the regression equation (Y = a0 + a1x1 + a2x2 + a3x2+ a4x1x2 + a5x1x3 + a6x2x3 + a7x12 + a8x22 + a9x32) for % ryegrass cover, turf height (cm), turf color, turf quality, turf density, winter-kill, fresh weight (g/dm2), dry weight (g/dm2), and spring transition during the first experiment (2014–2015) and the second experiment (2015–2016). x1, x2, x3 are the coded values of the fertilizer rates for N, P, and K.

Turf Height

For turf height, the coefficients a1 and a2 were positive with a1 > a2. In contrast, the coefficient a3 was negative during both years (Table 4). The lowest turf heights were recorded in the control plots (Table 3) during both experiments and were 4.57 and 4.67 cm. In both experiments, the highest turf heights of 13.68 and 14.10 cm (Table 3) were observed when N, P = 23.918 and K = 6.082.

Turf Density

N and P enhanced the turf density 10-fold in the first experiment and 17-fold in the second experiment and a1 > a2. In contrast, K had no effect (Table 4). The highest turf density of 24.27 and 27.53 were achieved at N, P = 23.918, K = 6.082 treated plots during both experiments (Table 3). In the first experiment, the control treatment and in the second experiment, the N = 0 and P, K = 15 treatment yielded the lowest turf densities of 6.93 and 8.33, respectively (Table 3).

Turf Color

The parameter coefficients for turf color a1 (maximum among all) and a2 were positive during the study period. In contrast, a3 (minimum among all) was positive and negative in consecutive years. (Table 4). Poorest color ratings of 4.23 and 4.27 were observed in the control plots (Table 3). We observed the highest color coefficients of 8.56 and 8.53 (Table 3) in response to the N = 30 and P, K = 15 treatment during the first experiment and in response to the N, P = 23.918 and K = 6.082 treatment during the second experiment.

Turf Quality

Turf quality of the overseeded ryegrass was similar during both experiments. In other words, a1, a2, and a3 were all positive with N having the greatest effect, followed by P and then K (Table 4). The N, P = 23.918 and K = 6.028 treatment yielded the highest turf quality of 8.68 and 8.82, and the control plots yielded the lowest turf quality of 3.49 and 3.44 during both years (Table 3).


Throughout both experiments, the ryegrass survived well during winter before the onset of snow fall. The worst winter-killings were detected in unfertilized or zero N plots (Table 3). In general, all of the three factors ameliorated the worst effects of low temperature on ryegrass (Table 4), with the relative magnitudes of the absolute values rank as follows: a1 > a2 > a3.

Fresh and Dry Weight (g dm−2)

During both experiments, the a1 and a2 coefficients were positive and a1 > a2. In contrast, the a3 coefficient was negative (Table 4). For both the fresh and dry weights, the highest yields from both experiments were obtained from the N, P = 23.918 and K = 6.082 treated plots (Table 3). The lowest yields were obtained from the control plots during both experiments (Table 3).

Spring Transition

For the spring transition, in the first experiment, a1 and a2 were negative, and a3 was positive (Table 4). For the second experiment, a1 was negative and both a2 and a3 were positive. The spring transition to bermudagrass was generally more effective when fertilization was reduced.

Gas-Exchange, Protein, and Total Chlorophyll

We analyzed a variety of physiological parameters for the ryegrass under different fertilization treatments (Table 5) and calculated regression coefficients for each parameter (Table 6). All coefficients for Pn were positive. The highest values we obtained for Pn were 6.25 and 4.27 (mmol m−2 s−1) from the N, P, K = 23.918 and N, P = 23.918, K = 6.082 treatments, respectively. The lowest values obtained for Pn were 0.76 and 0.72 (mmol m−2 s−1) from the control conditions during both years. Absolute values showed that N was more effective than P and K. The regression equation coefficients for Tr were positive. We observed a maximum of 2.32 and 1.29 (mmol m−2 s−1) from the N, K = 23.918, P = 6.082, and N = 23.918, P, K = 6.082 treatments, respectively, and a minimum of 0.16 and 0.17 (mmol m−2 s−1) from the control conditions. We observed a similar pattern for Ci and Gs. The maximum values for Ci 319.36 and 292.05 (mmol mol−1) and Gs 0.37 and 0.13 (mol m−2 s−1) in the N, K = 23.918, P = 6.082 conditions, and the lowest values for Ci and Gs in the control conditions (Table 5). Although a1 (N) and a3 (K) positively influenced Ci in the first experiment, only a1 (N) effectively improved Ci in the second experiment. Gs was positively affected by all the three factors in the first experiment, whereas only a2 (P) negatively affected in the second experiment. In contrast to these data, we observed the highest values for WUE in the treatments with lower quantities of fertilizer and lowest values for WUE in the treatments with higher quantities of fertilizer.


Table 5. Physiological analysis of perennial ryegrass that was overseeded on bermudagrass using different treatment combinations of nitrogen, phosphorus, and potassium fertilizers in two experiments (i.e., 2014–2015 and 2015–2016).


Table 6. Parameter coefficients of the regression equation (Y = a0 + a1x1 + a2x2 + a3x2+ a4x1x2 + a5x1x3 + a6x2x3 + a7x12 + a8x22 + a9x32) for photosynthesis (Pn), transpiration (Tr), inter cellular CO2 concentration (Ci), stomal conductance (Gs), water use efficiency (WUE), protein levels and total chlorophyll during the first experiment (2014–2015) and the second experiment (2015–2016). x1, x2, x3 are the coded values of the fertilizer rates for N, P, and K.

We observed highest concentration of protein 2.87 and 2.40 (g L−1) in extracts prepared from the N, P = 23.918, K = 6.082 treated plants and that of lowest concentrations of protein 0.26 and 0.18 (g L−1) in the control conditions. Regression coefficients indicated that although N and P consistently boosted protein concentration, K fertilization had no significant effects in the first experiment.

Total chlorophyll content in leaves was effectively enhanced by the NPK fertilization (first experiment) and NP fertilization (second experiment). N fertilization enhanced chlorophyll content the most (Table 6). We observed the maximum levels of chlorophyll 1.90 and 1.91 (mg Chl g−1 fw) in response to the N = 30, P, K = 15 and N, P = 23.918, K = 6.082 treatments and the lowest levels of chlorophyll 0.82 and 0.18 (mg Chl g−1 fw) in the control conditions (Table 5).

Effects of N, P, and K on Integrated Turf Performance (ITP)

The data in (Table 3) was transformed according to the experimental design using Equation (4) yielding (Tables 7, 8).

X=x-xminxmax- xmin    (4)

Table 7. Transformed values of morphological parameters according to the central composite rotatable design and integrated turf performance (ITP) of overseeded perennial ryegrass on bermudagrass from 2014 to 2015.


Table 8. Transformed values of morphological parameters according to the central composite rotatable design and integrated turf performance (ITP) of overseeded perennial ryegrass on bermudagrass from 2015 to 2016.

In Equation (4), X is the coded value for each treatment, x is the measured value of each treatment, xmin is the minimum, and xmax is the maximum value recorded for each parameter.

The ITP of perennial ryegrass was best (0.787 and 0.803) when treated with medium-high N, P and medium-low K (N, P = 23.918 and K = 6.082). The ITP was the worst (0.316 and 0.320) when the perennial ryegrass was treated with zero N and medium P and K (N = 0, P, K = 15) during both experiments. These data indicate that N-containing fertilizer has a substantial effect on overall turf quality (Tables 7, 8).

Regression Model of ITP Under N, P, and K

We applied regression analysis to the values we obtained for ITP from the first experiment and a quadratic regression model Equation (5) was generated in which Y represents ITP.

Y=0.617+0.128x1+0.024x2-0.005x3+0.016x1x2       -0.01x1x3-0.015x2x3-0.048x12-0.031x22-0.035x32    (5)

In Equation (5), the terms x1, x2, x12, x22, and x32 were highly significant.

Using the same protocol, we generated a regression model Equation (6) for second experiment:

Y=0.622+0.125x1+0.027x2-0.004x3+0.015x1x2       -0.06x1x3-0.022x2x3-0.041x12-0.028x22-0.023x32    (6)

In Equation (6), the terms x1, x2, x2x3, x12, x22, and x32 were highly significant.

The coefficient of determination r2 (0.98 and 0.97) of the regression Equations (5) and (6), respectively, were significant with p ≤ 0.001. This shows that the r2 of the regression equations were highly significant during both years, signifying that the models accurately predicted ITP and that approximately 98 and 97% of the variations in the responses could be explained by the models. The lower value of “P” (probability) indicates a highly significant regression model (Halder et al., 2015).

Interactive Effects on ITP

To investigate the interactive effects of the three factors using CCRD, 3D surface plots were generated by setting one variable at zero in the regression models (Equations 5 and 6). The interaction between the two factors (x and y) with ITP at the z-axis is shown in (Figure 2) for both experiments.


Figure 2. Response surfaces showing the effects of fertilization with N and P, N and K, and P and K on ITP for overseeded perennial ryegrass in (A) the first experiment and (B) the second experiment.

Effect of x1 (N) and x2 (P)

Figure 2 shows that ITP first increased and then decreased as N and P increased. The optimal position for ITP is where its value remained constant. The effect of N fertilizer on ITP was obvious for both experiments. ITP rapidly increased as N fertilizer increased either at lower or high P rates followed by decline with further increase in N fertilizer. On the other hand, at lower N rates the increasing rates of P did not markedly increase ITP, whereas at higher rates of N with increasing P rates, the ITP increased distinctly, followed by gradual decline with consistent increase in N rates.

Effect of x1 (N) and x3 (K)

The response surface plots of N and K fertilization in Figure 2 indicate that at lower amount of K, incremental increase in N led to rapid increases in ITP. When the amount of K was high, ITP first increased and then decreased with a continuous increase in N. In contrast, at low levels of N and increased levels of K, the lowest ITP was recorded. ITP reached a maximum at high levels of N and increased levels of K. Further increases in the levels of N and K induced decreases as indicated by red-colored and shaded areas (Figure 2).

Effect of x2 (P) and x3 (K)

The interactive effects of fertilization with P and K depended on the amounts of each nutrient. ITP slightly increased at low levels of P and increasing levels of K, whereas it gradually decreases for K added above a particular level. In contrast, when P was added at high levels, the maximum values for ITP were obtained with increasing amounts of K up to a particular range. On the other hand, when low levels of K were added and the amounts of P were also increasing, ITP markedly increased. However, when K was added at high levels, the ITP decreased when P was increased consistently.

According to the qualitative analysis of factor interactions, it is clear that N is the most important nutrient, followed by P and K, respectively. To achieve a higher ITP, N fertilizer should be applied in high amounts, P fertilizer in medium to high amounts, and K fertilizer in low amounts.


The optimization of N, P, and K fertilization was performed using regression models (Equations 5 and 6) based on a computer simulation of the formula according to the regression model, within the range of each variable [−1. 682~1. 682], based upon a step length of “0.3346.” A total of 1,331 combinations and their corresponding ITP values were obtained, with the maximum value of 0.72 and 0.75 in first and second experiment, respectively, which corresponded to actual fertilizer values of N: 30, P: 24, K: 9, and N: 30, P: 27, K: 6 g m−2, respectively. The lowest ITP values of 0.02 and 0.07 were obtained in the control conditions during both experiments.


Nutrients, especially N, P, and K play a vital role in plant growth and production, as they are involved in many biological processes (Deng et al., 2014). The vegetative growth of perennial ryegrass is vigorous and hence, perennial ryegrass has high nutritional demands (Turgeon, 2008). The current study confirms that there are substantial effects of N, P, and K fertilization on the morphology, physiology, and ITP of overseeded perennial ryegrass on bermudagrass during the winter. During this study, the effect of N was more obvious followed by P and then K, (i.e., N>P>K).

In this study, we observed that applying N, P, and K fertilizer increased the percent of ryegrass cover in two distinct experiments (Table 4). The fact that N, P, and K are major macronutrients nutrients that are required for growth and development may explain these findings. N and P are required for the biosynthesis of abundant and essential biomolecules, such as proteins and nucleic acids, respectively, while K is the major cation of plant cells. Indeed, N-deficient soils often limit plant growth (Rajasekar et al., 2017). Consistent with our findings, fertilization with N, P, and K was previously reported to increase grass cover (Rodriguez et al., 2002; Guertal and Evans, 2006). We observed that fertilization with N and P significantly increased turf height and density, which defines the growth speed and the winter survival rate. A higher growth rate from tillers is a highly desirable turfgrass trait that improves aesthetics, functionality, and recovery from stress (Turgeon, 2008). We demonstrated that fertilization with N and P promoted higher rates of photosynthesis, which drives higher growth rates. Earlier reports demonstrated that the interactive effects of N and P enhanced that uptake and efficiencies of N and P in rice, wheat and grasses (Aulakh and Malhi, 2005) and documented similar findings in turfgrasses (Ahmad et al., 2003; Schlossberg and Schmidt, 2007).

Fertilization with N was major factor promoting the accumulation of chlorophyll in both experiments. In contrast, fertilization with P and K only promoted the accumulation of chlorophyll in the first experiment. Fertilization with N probably promotes the accumulation of chlorophyll in part because N is a major component of chlorophyll (Munshaw et al., 2007; Çelebi et al., 2010). Indeed, N-deficient plants typically exhibit chlorophyll-deficient phenotypes. K and P are not components of chlorophyll. In the first experiment, K and P probably promoted increases in the levels of chlorophyll by an indirect mechanism that is related to the general growth promoting properties of these important macronutrients. Our finding that fertilization with P and K has a minor effect on the accumulation of chlorophyll in turfgrass is consistent with previous work (Baldi et al., 2013). Consistent with previous work (Snyder and Cisar, 2000; McCarty and Miller, 2002), we found that fertilization with N, P, and K affected turf quality, which is an index that considers density, color, growth, and stress tolerance. Fertilization with N, P, and K protected ryegrass from cold stress (i.e., winter-kill). Indeed, plots that received higher rates of these macronutrients performed well. Our finding that fertilization with N, P, and K induced increases in growth during periods of cold stress might be explained by increases in water potential. The maintenance of water potential under high N conditions that included P and K was previously reported to promote leaf processes during stress conditions (Abid et al., 2016). Additionally, high rates of fertilization with N, P, and K were previously reported to promote winter hardiness (Razmjoo and Kaneko, 1993).

Fresh and dry weight of overseeded ryegrass increased as the amounts of N and P fertilizer increased. The increased shoot biomass was mainly due to enhanced lateral growth (i.e., density and height), as discussed above, that maximized leaf area thereby promoting the conversion of solar energy into chemical energy, which is essential for the production of biomass. Photosynthesis is the key driving force for the accumulation of dry matter, organ formation, and plant production (Zlatev and Lidon, 2012). The increase in the leaf area index, leaf photosynthetic area, and the accumulation of carbohydrates enhanced overall plant biomass (Din et al., 2017). Indeed, fertilization with N and P was reported to increases the accumulation of fresh and dry matter (Snyder and Cisar, 2000; Rodriguez et al., 2002; Stanford et al., 2005; Guertal, 2006).

High levels of fertilization attenuated the transition from ryegrass in the winter to bermudagrass in the spring. Generally, smooth, and rapid transitions occurred in plots that received lower applications of N and higher applications of both P and K. One explanation might be that because of the broad impact of N on plant growth and development, high N enhances the vigor of ryegrass by promoting processes that drives growth, such as photosynthesis, and promotes various types of stress tolerance, such as thermotolerance. Indeed, difficulties in the spring transition were previously attributed to the enhanced thermotolerance of new perennial ryegrass cultivars (Mazur and Wagner, 1987). The reduced availably of P and K, vertical mowing, the increase in temperature, and the additional N application in late April enhanced the recovery of bermudagrass. Consistent with these data, K fertilization was previously reported to promote an early spring transition (Duble, 2001; Webster and Ebdon, 2005).

Cold stress often decreases photosynthesis because of a variety of factors, such as stomatal behavior, photoinhibition, and the inhibition of carbon fixation (Wise et al., 2004). Fertilizers containing N, P, and K are well known for affecting photosynthesis in leaves (Wang et al., 2008). Deficiency and excess quantities of these fertilizers decrease the capacity of photosynthesis (Ashraf et al., 2001; Cechin and de Fátima Fumis, 2004). In general, fertilization with N, P, and K induced increase in the gas-exchange parameters of overseeded ryegrass during the winter. Higher rates of transpiration and the regulation of stomatal opening may induce the rate of photosynthesis. Greater biomass, elevated chlorophyll and higher protein content that were observed in response to NPK fertilization might be other factors that drove increases in Pn, Tr, and Gs (Alderman et al., 2011; Grygierzec, 2012). P is a vital structural part of nucleic acids (Wang et al., 2013) and promotes photosynthesis (Fredeen et al., 1990). Elevated levels of P could promote photosynthesis by contributing to the production of the products of noncyclic electron transport, adenosine triphosphate (ATP) and NADPH (Giersch and Robinson, 1987). K was suggested to induce increases in photosynthesis by regulating transpiration, stomatal behavior and turgor pressure (Römheld and Kirkby, 2010).

Elevated values of Ci were due to increased Pn and Gs. These data indicate that overseeded ryegrass leaves were capable of enhancing CO2 influx and fixing CO2 into carbohydrates. K fertilization regulates stomatal functions and N fertilization increases Pn. We suggest that acting together, these nutrients induce increases in CO2 fixation. Higher levels of N, P, and K reduced the WUE. The coefficient values indicated that P increased the WUE only in the second experiment. In our study, we found that increases in Tr were functions of N and K and that increases in Tr reduced the WUE. The lower WUE in the high-nutrient treated plot could also be attributed to lower increases in Pn relative to Tr (Fahad et al., 2016).

Protein and total chlorophyll were affected by N, P, and K in the first year and by N and P during second year (Table 6). The increase in protein and total chlorophyll in response to N, P, and K were correlated with increases in photosynthesis. This relationship is consistent with chlorophyll serving as the main light-harvesting pigment for photosynthesis and with the observation that proteins associated with photosynthesis accumulate to extremely high levels, such as Rubisco and the light-harvesting chlorophyll a/b-binding proteins. Another reason may be that all of the three factors, especially N decreased stress induced damage by regulating enzymatic and non-enzymatic activities (Hussain et al., 2016). Thus, the plants were required to devote fewer resources to stress tolerance and therefore, they were able to invest more resources in processes, such as photosynthesis, that drive growth and development (Abid et al., 2016). N was reported to induce leaf development and photosynthetic efficiency (Dordas and Sioulas, 2008), P was reported to encourage carbon assimilation and the phosphorylation that enhances the regeneration of Rubisco (Liu et al., 2010), and K was reported to regulate stomata (Tsonev et al., 2011). We suggest that all of these factors contributed to higher protein and chlorophyll concentrations that we observed. Thus, our findings are consistent with previous work that documented increases in protein and chlorophyll levels following fertilization with N and P (Dordas and Sioulas, 2008; Nemadodzi et al., 2017).


The results of this study demonstrated that different combinations of N, P, and K affected the morphology and physiology of overseeded perennial ryegrass on bermudagrass with a number of interactive effects. N, P, and K had positive effects on morphological traits with the exception of spring transition and physiological traits with the exception of WUE, and had interactive effects on ITP. The magnitudes of these nutrients effects on ITP followed the order of: N>P>K. Regarding the morphological attributes, growth parameters were mostly affected by fertilization with N and P, and qualitative parameters were affected by fertilization with N, P, and K. N and P delayed the spring transition while, K accelerated the spring transition. Fertilization with N, P, and K markedly increased Pn, Tr, Gs, levels of protein and total chlorophyll. Regression models between ITP and the three nutrients were established, and N, P, and K were optimized by simulation. The optimized scheme with maximum turf performance (ITP) was N: 30, P: 24, K: 9 g m−2 in first experiment and N: 30, P: 27, K: 6 g m−2 in the second experiment. Based on our results, we conclude that this approach can be effectively used to establish nutrient management strategies to maximize performance and yield. To cope with cold climates, turfgrasses management will depend not only on breeding new varieties that are well adopted to cold conditions, but also on the development of nutrient management practices, that will contribute to sustainable crop production.

Author Contributions

LC and SY: conceived the idea of the experiment; MI: performed the experiment and wrote manuscript; RL: revised and edited the manuscript; SF and TL: contributed to the collection of the data and statistical analysis; RL, SY, and LC: provided the technical guidance and editing support.


This research work was supported by Fundamental Research Funds for the Central Universities, China (Grant No. 2662016QD016).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


Abid, M., Tian, Z., Ata-Ul-Karim, S. T., Cui, Y., Liu, Y., Zahoor, R., et al. (2016). Nitrogen nutrition improves the potential of wheat (Triticum aestivum L.) to alleviate the effects of drought stress during vegetative growth periods. Front. Plant Sci. 7:981. doi: 10.3389/fpls.2016.00981

PubMed Abstract | CrossRef Full Text | Google Scholar

Ahmad, I., Khan, M., and Qasim, A. (2003). Growth and development of different turfgrasses as influenced by nitrogen application and leaf nitrogen contents. Int. J Agric. Biol. 5, 175–178.

Alderman, P., Boote, K., Sollenberger, L., and Coleman, S. (2011). Carbohydrate and nitrogen reserves relative to regrowth dynamics of ‘Tifton 85’bermudagrass as affected by nitrogen fertilization. Crop Sci. 51, 1727–1738. doi: 10.2135/cropsci2010.09.0516

CrossRef Full Text | Google Scholar

Ashraf, M., Ahmad, A., and McNeilly, T. (2001). Growth and photosynthetic characteristics in pearl millet under water stress and different potassium supply. Photosynthetica 39, 389–394. doi: 10.1023/A:1015182310754

CrossRef Full Text | Google Scholar

Aulakh, M. S., and Malhi, S. S. (2005). Interactions of nitrogen with other nutrients and water: effect on crop yield and quality, nutrient use efficiency, carbon sequestration, and environmental pollution. Adv. Agron. 86, 341–409. doi: 10.1016/S0065-2113(05)86007-9

CrossRef Full Text | Google Scholar

Baldi, A., Lenzi, A., Nannicini, M., Pardini, A., and Tesi, R. (2013). Growth and nutrient content of hybrid bermudagrass grown for nursery purposes at different nitrogen, phosphorus, and potassium rates. Horttechnology 23, 347–355.

Google Scholar

Beard, J., and Menn, W. (1988). The Texas System of Winter Overseeding. Grounds Maintenance (USA).

Cechin, I., and de Fátima Fumis, T. (2004). Effect of nitrogen supply on growth and photosynthesis of sunflower plants grown in the greenhouse. Plant Sci. 166, 1379–1385. doi: 10.1016/j.plantsci.2004.01.020

CrossRef Full Text | Google Scholar

Çelebi, Z., Andiç, N., and Yilmaz, H. (2010). Determination of proper species mixtures for established turfgrass field in Van region. YYU. J. Agr. Sci. 20, 16–25.

Deng, G., Liu, L. J., Zhong, X. Y., Lao, C. Y., Wang, H. Y., Wang, B., et al. (2014). Comparative proteome analysis of the response of ramie under N, P and K deficiency. Planta 239, 1175–1186. doi: 10.1007/s00425-014-2040-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Din, M., Zheng, W., Rashid, M., Wang, S., and Shi, Z. (2017). Evaluating hyperspectral vegetation indices for leaf area index estimation of Oryza sativa L. at diverse phenological stages. Front. Plant Sci. 8:820. doi: 10.3389/fpls.2017.00820

PubMed Abstract | CrossRef Full Text | Google Scholar

Dordas, C. A., and Sioulas, C. (2008). Safflower yield, chlorophyll content, photosynthesis, and water use efficiency response to nitrogen fertilization under rainfed conditions. Ind. Crop Prod. 27, 75–85. doi: 10.1016/j.indcrop.2007.07.020

CrossRef Full Text | Google Scholar

Duble, R. L. (2001). Turfgrasses: Their Management and Use in the Southern Zone. College Station, TX: Texas A&M University Press.

Google Scholar

Esmaeili, S., Salehi, H., and Eshghi, S. (2015). Silicon ameliorates the adverse effects of salinity on turfgrass growth and development. J. Plant Nutr. 38, 1885–1901. doi: 10.1080/01904167.2015.1069332

CrossRef Full Text | Google Scholar

Fahad, S., Hussain, S., Saud, S., Hassan, S., Tanveer, M., Ihsan, M. Z., et al. (2016). A combined application of biochar and phosphorus alleviates heat-induced adversities on physiological, agronomical and quality attributes of rice. Plant Physiol. Biochem. 103, 191–198. doi: 10.1016/j.plaphy.2016.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, J., Ren, J., Zhu, W., Amombo, E., Fu, J., and Chen, L. (2014). Antioxidant responses and gene expression in bermudagrass under cold stress. J. AM. Soc. Hortic. Sci. 139, 699–705. doi: 10.1007/s10646-016-1696-9

CrossRef Full Text | Google Scholar

Fredeen, A. L., Raab, T. K., Rao, I. M., and Terry, N. (1990). Effects of phosphorus nutrition on photosynthesis in Glycine max (L.) Merr. Planta 181, 399–405. doi: 10.1007/BF00195894

PubMed Abstract | CrossRef Full Text | Google Scholar

Giersch, C., and Robinson, S. P. (1987). Regulation of photosynthetic carbon metabolism during phosphate limitation of photosynthesis in isolated spinach chloroplasts. Photosynth. Res. 14, 211–227. doi: 10.1007/BF00032706

PubMed Abstract | CrossRef Full Text | Google Scholar

Grygierzec, B. (2012). Productivity of selected grasses in mixtures with Trifolium repens L. at two levels of nitrogen fertilization. Fragm. Agron. 29, 31–36.

Google Scholar

Guertal, E. (2006). Phosphorus movement and uptake in bermudagrass putting greens. USGA Turfgrass Environ. Res. Online 5, 1–7.

Guertal, E., and Evans, D. (2006). Nitrogen rate and mowing height effects on TifEagle bermudagrass establishment. Crop Sci. 46, 1772–1778. doi: 10.2135/cropsci2006.01-0006

CrossRef Full Text | Google Scholar

Guo, X., Sun, W., Wei, W., Yang, J., Kang, Y., Wu, J., et al. (2009). Effects of NPK on physiological and biochemical characteristics of winter rapaseed in northwest cold and drought region. Acta Agric. Boreali-Occident. Sin. 2, 27.

Google Scholar

Halder, G., Dhawane, S., Barai, P. K., and Das, A. (2015). Optimizing chromium (VI) adsorption onto superheated steam activated granular carbon through response surface methodology and artificial neural network. Environ. Prog. Sustain. 34, 638–647. doi: 10.1002/ep.12028

CrossRef Full Text | Google Scholar

Hons, F., McFarland, M. L., Lemon, R. G., Nichols, R. L., Mazac, F. Jr., Boman, R., et al. (2004). Managing Nitrogen Fertilizer in Cotton. College Station, TX: Texas A&M University Press.

Google Scholar

Hussain, S., Khan, F., Cao, W., Wu, L., and Geng, M. (2016). Seed priming alters the production and detoxification of reactive oxygen intermediates in rice seedlings grown under sub-optimal temperature and nutrient supply. Front. Plant Sci. 7:439. doi: 10.3389/fpls.2016.00439

PubMed Abstract | CrossRef Full Text | Google Scholar

Jellicorse, W. R., Richardson, M. D., McCalla, J. H., Karcher, D. E., Patton, A. J., and Boyd, J. W. (2012). Seeded bermudagrass establishment in an overseeded perennial ryegrass stand as affected by transition herbicide and seeding date. Appl. Turfgrass Sci. 9. doi: 10.1094/ATS-2012-0721-01-RS

CrossRef Full Text | Google Scholar

Jiang, Y., and Huang, B. (2000). Effects of drought or heat stress alone and in combination on Kentucky bluegrass. Crop Sci. 40, 1358–1362. doi: 10.2135/cropsci2000.4051358x

CrossRef Full Text | Google Scholar

Larcher, W. (2003). “Plants under stress,” in Physiological Plant Ecology, Ecophysiology and Stress of Functional Groups, 4th Edn, (Berlin; Heidelberg: Springer), 345–450.

Google Scholar

Liu, H., Hu, C., Sun, X., Tan, Q., Nie, Z., and Hu, X. (2010). Interactive effects of molybdenum and phosphorus fertilizers on photosynthetic characteristics of seedlings and grain yield of Brassica napus. Plant Soil 326, 345–353. doi: 10.1007/s11104-009-0014-1

CrossRef Full Text | Google Scholar

Mazur, A., and Wagner, D. (1987). Influence of aeration, topdressing, and vertical mowing on overseeded bermudagrass putting green turf. HortScience 22, 1276–1278.

Google Scholar

McCarty, L. B., and Miller, G. (2002). Managing Bermudagrass Turf: Selection, Construction, Cultural Practices, and Pest Management Strategies. Ann Arbor, MI: John Wiley & Sons.

Google Scholar

Montgomery, D. C. (2001). Design and Analysis of Experiments. New York, NY: John Wiley & Sons.

Google Scholar

Munshaw, G., Ervin, E., Parrish, D., Shang, C., Askew, S., Zhang, X., et al. (2007). Influence of late-season iron, nitrogen, and seaweed extract on fall color retention and cold tolerance of four bermudagrass cultivars. Crop Sci. 47, 463–463. doi: 10.2135/cropsci2007.01.0001er

CrossRef Full Text | Google Scholar

Nemadodzi, L., Araya, H., Nkomo, M., Ngezimana, W., and Mudau, N. F. (2017). Nitrogen, phosphorus and potassium effects on the physiology and response biomass yield of baby spinach (Spinacia oleracea L.). J. Plant Nutr. 40, 2033–2044. doi: 10.1080/01904167.2017.1346121

CrossRef Full Text | Google Scholar

Pompeiano, A., Guglielminetti, L., and Volterrani, M. (2011). Freeze tolerance of Zoysia matrella (L.) Merrill as affected by late-season nitrogen application, and changes in carbohydrates during cold acclimation. Plant Biosyst. 145, 885–892. doi: 10.1080/11263504.2011.623191

CrossRef Full Text | Google Scholar

Rajasekar, M., Nandhini, D. U., Swaminathan, V., and Balakrishnan, K. (2017). A review on role of macro nutrients on production and quality of vegetables. Int. J. Chem. Sci. 5, 304–309.

Google Scholar

Razmjoo, K., and Kaneko, S. (1993). Effect of fertility ratios on growth and turf quality of perennial ryegrass (Lolium prenne L.) in winter. J. Plant Nutr. 16, 1531–1538. doi: 10.1080/01904169309364630

CrossRef Full Text | Google Scholar

Richardson, M. (2004). Morphology, turf quality, and heat tolerance of intermediate ryegrass. HortScience 39, 170–173.

Google Scholar

Rodriguez, I. R., Miller, G. L., and McCarty, L. (2002). Bermudagrass establishment on high sand-content soils using various NPK ratios. HortScience 37, 208–209.

Google Scholar

Römheld, V., and Kirkby, E. A. (2010). Research on potassium in agriculture: needs and prospects. Plant Soil 335, 155–180. doi: 10.1007/s11104-010-0520-1

CrossRef Full Text | Google Scholar

Schlossberg, M. J., and Schmidt, J. P. (2007). Influence of nitrogen rate and form on quality of putting greens cohabited by creeping bentgrass and annual bluegrass. Agron. J. 99, 99–106. doi: 10.2134/agronj2006.0136

CrossRef Full Text | Google Scholar

Shi, H., Jiang, C., Ye, T., Tan, D.-X., Reiter, R. J., Zhang, H., et al. (2014a). Comparative physiological, metabolomic, and transcriptomic analyses reveal mechanisms of improved abiotic stress resistance in bermudagrass [Cynodon dactylon (L). Pers.] by exogenous melatonin. J. Exp. Bot. 66, 681–694. doi: 10.1093/jxb/eru373

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, H., Ye, T., Zhong, B., Liu, X., and Chan, Z. (2014b). Comparative proteomic and metabolomic analyses reveal mechanisms of improved cold stress tolerance in bermudagrass (Cynodon dactylon (L.) Pers.) by exogenous calcium. J. Integr. Plant Biol. 56, 1064–1079. doi: 10.1111/jipb.12167

PubMed Abstract | CrossRef Full Text | Google Scholar

Snyder, G. H., and Cisar, J. L. (2000). Nitrogen/potassium fertilization ratios for bermudagrass turf. Crop Sci. 40, 1719–1723. doi: 10.2135/cropsci2000.4061719x

CrossRef Full Text | Google Scholar

Stanford, R., White, R., Krausz, J., Thomas, J., Colbaugh, P., and Abernathy, S. (2005). Temperature, nitrogen and light effects on hybrid bermudagrass growth and development. Crop Sci. 45, 2491–2496. doi: 10.2135/cropsci2005.0070

CrossRef Full Text | Google Scholar

Tsonev, T., Velikova, V., Yildiz-Aktas, L., Gürel, A., and Edreva, A. (2011). Effect of water deficit and potassium fertilization on photosynthetic activity in cotton plants. Plant Biosyst. 145, 841–847. doi: 10.1080/11263504.2011.560199

CrossRef Full Text | Google Scholar

Turgeon, A. J. (2008). Turfgrass Management, 8th Edn. Upper Saddle River, NJ; Columbus, OH: Pearson Prentice Hall.

Tyler, N., Gusta, L., and Fowler, D. (1981). The influence of nitrogen, phosphorus and potassium on the cold acclimation of winter wheat (Triticum aestivum L.). Can. J. Plant Sci. 61, 879–885. doi: 10.4141/cjps81-131

CrossRef Full Text | Google Scholar

Uleberg, E., Hanssen-Bauer, I., van Oort, B., and Dalmannsdottir, S. (2014). Impact of climate change on agriculture in Northern Norway and potential strategies for adaptation. Clim. Change 122, 27–39. doi: 10.1007/s10584-013-0983-1

CrossRef Full Text | Google Scholar

Vance, C. P. (2001). Symbiotic nitrogen fixation and phosphorus acquisition. Plant nutrition in a world of declining renewable resources. Plant Physiol. 127, 390–397. doi: 10.1104/pp.010331

PubMed Abstract | CrossRef Full Text | Google Scholar

Verdenal, A., Combes, D., and Escobar-Gutierrez, A. J. (2008). A study of ryegrass architecture as a self-regulated system, using functional-structural plant modelling. Funct. Plant Biol. 35, 911–924. doi: 10.1071/FP08050

CrossRef Full Text | Google Scholar

Wang, J., Hui, D., Ren, H., Liu, Z., and Yang, L. (2013). Effects of understory vegetation and litter on plant Nitrogen (N), Phosphorus (P), N: P Ratio and their relationships with growth rate of indigenous seedlings in subtropical plantations. PLoS ONE 8:84130. doi: 10.1371/journal.pone.0084130

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Yang, J. F., Han, X. R., Liu, X. H., Zhan, X. M., and Liu, S. G. (2008). Effect of fertilizer application on photosynthetic traits of spring maize. Soil Fert. Sci. Chin. 6, 23–27.

Google Scholar

Webster, D., and Ebdon, J. (2005). Effects of nitrogen and potassium fertilization on perennial ryegrass cold tolerance during deacclimation in late winter and early spring. HortScience 40, 842–849.

Google Scholar

Wise, R., Olson, A., Schrader, S., and Sharkey, T. (2004). Electron transport is the functional limitation of photosynthesis in field-grown pima cotton plants at high temperature. Plant Cell Environ. 27, 717–724. doi: 10.1111/j.1365-3040.2004.01171.x

CrossRef Full Text | Google Scholar

Yang, Z., Miao, Y., Yu, J., Liu, J., and Huang, B. (2014). Differential growth and physiological responses to heat stress between two annual and two perennial cool-season turfgrasses. Sci. Hortic. Amsterdam 170, 75–81. doi: 10.1016/j.scienta.2014.02.005

CrossRef Full Text | Google Scholar

Zlatev, Z., and Lidon, F. C. (2012). An overview on drought induced changes in plant growth, water relations and photosynthesis. Emir. J. Food Agric. 24, 57–72. doi: 10.9755/ejfa.v24i1.10599

CrossRef Full Text | Google Scholar

Keywords: turfgrass, overseeding, fertilizer optimization, central composite rotatable design, integrated turf performance, modeling

Citation: Ihtisham M, Fahad S, Luo T, Larkin RM, Yin S and Chen L (2018) Optimization of Nitrogen, Phosphorus, and Potassium Fertilization Rates for Overseeded Perennial Ryegrass Turf on Dormant Bermudagrass in a Transitional Climate. Front. Plant Sci. 9:487. doi: 10.3389/fpls.2018.00487

Received: 24 October 2017; Accepted: 29 March 2018;
Published: 16 April 2018.

Edited by:

Abraham J. Escobar-Gutiérrez, Institut National de la Recherche Agronomique (INRA), France

Reviewed by:

Corina Carranca, Instituto Nacional de Investigação Agrária e Veterinária, Portugal
Fu Gui Mi, Inner Mongolia Agricultural University, China

Copyright © 2018 Ihtisham, Fahad, Luo, Larkin, Yin and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Shaohua Yin,
Longqing Chen,