Next Article in Journal
The Effect of Trichoderma citrinoviride Treatment under Salinity Combined to Rhizoctonia solani Infection in Strawberry (Fragaria x ananassa Duch.)
Next Article in Special Issue
First Steps to Domesticate Hairy Stork’s-Bill (Erodium crassifolium) as a Commercial Pharmaceutical Crop for Arid Regions
Previous Article in Journal
Licorice, Doum, and Banana Peel Extracts Inhibit Aspergillus flavus Growth and Suppress Metabolic Pathway of Aflatoxin B1 Production
Previous Article in Special Issue
Cultivation of Desert Truffles—A Crop Suitable for Arid and Semi-Arid Zones
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tef (Eragrostis tef) Responses to Phosphorus and Potassium Fertigation under Semi-Arid Mediterranean Climate

1
Gilat Research Center, Agricultural Research Organization, Negev 85280, Israel
2
The Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem, Rehovot 76100, Israel
*
Author to whom correspondence should be addressed.
Equal contribution.
Agronomy 2021, 11(8), 1588; https://doi.org/10.3390/agronomy11081588
Submission received: 27 June 2021 / Revised: 26 July 2021 / Accepted: 5 August 2021 / Published: 10 August 2021
(This article belongs to the Special Issue New Crops for Arid Regions)

Abstract

:
Tef (Eragrostis tef (Zucc.) Trotter) is an annual small grain, panicle bearing, C4 cereal crop native to Ethiopia, where it is a major staple food. The objectives of the present study were to characterize the responses of two tef genotypes to escalating phosphorus (P) and potassium (K) levels and to determine an optimum range for P and K at which tef performance is maximized. Two experiments were carried out in the Gilat Research Station, each testing two different genotypes of tef (405B and 406W), one experiment in pots in controlled conditions, and the other in the field. In both experiments, the highest grain yield increased until 6 mg L−1 P, and declined at 12 mg L−1 P. The decline was precipitous and significant in the pot experiment, and gradual and statistically insignificant in the field experiment. In the pots experiment, the grain yield increased until 40 mg L−1 K, with no significant decrease thereafter. The effect of K concentration was only seen in the grain yield and not in the size of the other plant organs. In the field experiment, grain yield was highest at 80 mg L−1 K, but it was not statistically different from 40 mg L−1. The effect of K on growth was only apparent at maturity and not at flowering.

1. Introduction

Tef (also written as “teff”) (Eragrostis tef (Zucc.) Trotter) is an annual small grain, panicle bearing, C4 cereal crop native to Ethiopia, where it is a major staple food, accounting for two thirds of the daily protein intake of the population, and 11% of the caloric intake [1]. It is resistant to many biotic and abiotic stresses, making it an attractive crop for smallholder farmers [1]. It is known for its many health benefits, including its low glycemic index and high concentrations of essential amino acids, and it is also gluten free [2].
The majority of tef cultivation is located in Ethiopia [1], and it is mostly cultivated under low input, rainfed conditions However, in recent years interest has increased in cultivating tef outside of its native environment [1], including drylands where irrigation is required. Interest has also increased in cultivating tef in a more intensive manner [1,3]. One way of intensifying tef cultivation would be to apply water and fertilizer through irrigation (fertigation), which could allow for a high yield even in conditions of low rainfall or low soil fertility. Applying nutrients through fertigation requires a detailed understanding of how the crop responds to different concentration of nutrients, of which nitrogen (N), phosphorus (P) and potassium (K) are the most prominent. The response of tef to N fertilizer has been described in a previous paper [4], and in this paper we will discuss the response of tef to P and K.
Phosphorus fertilization of tef is common in Ethiopia [5]. A blanket recommendation by the Ethiopian Ministry of Agriculture was to apply between 30 and 40 kg P ha−1, depending on the soil type [5], which was later changed to 100 kg DAP ha−1 regardless of the soil type. In recent years, soil tests have been used to improve accuracy of P fertilization [5]. A number of field trials have been published in recent years showing the effect of different levels of P on tef growth parameters. Dubale [6] conducted a field experiment in the North Shewa Highlands of Ethiopia with three different rates of P fertilizer (0, 50, and 100 kg P2O5 ha−1) and found that broadcasting different amounts of P fertilizer at sowing did not significantly affect P concentration in plant tissue or grain yield, perhaps because there was already sufficient P in the soil in which the experiment was conducted or the P that was applied became unavailable to plants. Dereje et al. [7] conducted a full factorial field experiment in the Assosa Zone in Ethiopia in which five levels of N and four levels of P were tested. They found significant main effects and interactive effects of N and P on many different growth parameters including yield. After economic analysis, they concluded that in soils similar to their test plots, the best practice would be to apply 10 Kg P ha−1, along with 46 Kg N ha−1.
Although K is universally recognized as important for crop growth and yield, for many years tef grown in Ethiopia was usually fertilized only with N and P and not with K [8]. This was due to the common belief, stemming from studies conducted in the 1960s and 1970s, that claimed that K was not limiting in Ethiopian agricultural soils [8,9]. Recently, this idea was called into question with largescale mapping of Ethiopian soils in a project initiated by the Ethiopian government, called the Ethiopian Soil Information System (Ethiosis) [10]. A number of field trials have been recently published indicating that K application can increase tef yield [9,11]. Similarly, a greenhouse trial in which different levels of K were applied to tef growing in pots with different Ethiopian soils was conducted in order to determine an optimum K fertilization rate [12].
The purpose of this study is to elucidate the response of tef to different levels of P and K applied in the irrigation water. It was carried out in two phases. In the first phase, tef was grown under controlled conditions in a greenhouse, in perlite, an inert medium that allowed full control of the nutrient concentration in the root zone. This is in contrast to soil, which can absorb, fix, or release ions of K+ and PO4, making it difficult to determine the exact nutrient concentrations to which the roots are exposed. The second phase involved growing tef out in the field in a semi-arid Mediterranean climate in natural soil. While the first phase was designed to clarify the basic response of tef to different P and K concentrations, the second phase was intended to show how the lessons learned from the first phase play out in a realistic agricultural environment. The same experiments were used to study N fertigation in tef, recently published by Gashu et al. [4].

2. Materials and Methods

2.1. Plant Materials and Experimental Design

A pot experiment and a field experiment were carried out in the winter of 2015–2016 and the summer of 2016, respectively, in order to determine the response of tef to N, P and K nutrition. These experiments contained different segments in which different nutrients were varied, but which shared a control treatment. The N segment of these experiments has already been reported on by Gashu et al. [4]. The P and K segments of the experiments were conducted at the same time and under the same experimental conditions as the published N segment.
Two tef genotypes from the Israel Gene Bank (Agricultural Research Organization, Volcani Center, Rishon Lezion, Israel) were used in the current study, RTC-405 and RTC 406 (denoted hereafter as 405B and 406W, with the last letter indicating brown or white seed color). The entire collection was propagated and initially phenotyped in a common garden experiment during spring 2015 [13].
The pot experiment was conducted in a walk-in plastic-covered tunnel (6 m-wide, 2.4 m tall and 30 m-long) at the Gilat Research Center, Negev, Israel during winter 2015–2016. The P segment of the experiment comprised 4 treatments of P concentrations in the irrigation solution (1, 3, 6, and 12 mg L−1 P) and the K segment of the experiment comprised 4 treatments of K concentration in the irrigation solution (10, 20, 40, 80 mg L−1 K). Both segments included two tef genotypes (405B and 406W). A factorial (P or K treatments X genotypes) randomized block design was employed with five replicates. The P and K segments of the experiment shared a control group (N, P and K = 40, 6, and 40 mg L−1, respectively). The concentration of the other minerals in the fertigation solution of all treatments were: 40 mg L−1 N, 20 mg L−1 Ca, 20 mg L−1 Mg, 28 mg L−1 S, 0.3 mg L−1 B, 0.6 mg L−1 Fe, 0.3 mg L−1 Mn, 0.15 mg L−1 Zn, 0.02 mg L−1 Cu, and 0.02 mg L−1 Mo. No 0 P or 0 K treatments were included, since the perlite is an inert medium and the plants simply cannot grow without at least a small dose of P and K. More details regarding the environmental conditions and a timeline of the experiment are described in Gashu et al. [4]. A full list of treatments and replications, including those reported on by Gashu et al. [4], are listed in Appendix A.
The field experiment was conducted at the Gilat Research Center during the summer of 2016. The P segment of the experiment comprised 4 treatments of P concentrations in the irrigation solution (0, 3, 6, and 12 mg L−1 P) and the K segment of the experiment comprised 3 treatments of K concentration in the irrigation solution (0, 40 and 80 mg L−1 K). Both segments included two tef genotypes (405B and 406W). A factorial (P or K treatments × genotypes) split-plot block design with 5 replications was used, with fertigation treatment in the main plots and genotypes in subplots. Each main plot (5 m long × 4.2 m wide) consisted of 28 rows (14 rows per genotype) with 15 cm of space between rows. Each main plot was irrigated by 14 drip lines between each pair of rows. Seeds were directly sown on 13 July 2016 into well-prepared dry soil at a depth of ~1 cm, and seeding rate of 800 mg/m2, using a hand-driven precision garden seeder (1001B, Earthway, Bristol, IN, USA). During the first two weeks, plots were irrigated to drainage in order to avoid salt accumulation. Two weeks after sowing, fertigation treatments were started by injecting 1L of custom-made fertilizer solutions into 100 L of water. Fertigation was applied daily via a drip system with water amounts determined according to Penman Monteith potential evapotranspiration [13] multiplied by the following crop coefficients: 1.0 for the first 26 days after sowing, 1.2 for 26–30 days, 1.0 for 30–40 days, 1.2 for 40–60 days, 1.0 for 60–70 days and 0.8 from 70 days after sowing to the final harvest as described by Yihun et al. [14]. Both P and K segments shared a control group (N, P and K = 60, 6, 40 mg L−1, respectively). A full list of treatments and replications, including those reported on by Gashu et al. [4], are listed in Appendix A. Environmental conditions and a timeline of the experiment are described by Gashu et al. [4].

2.2. Data Collection

In both experiments, days to 50% panicle emergence was recorded for each pot and plot by visual observation. The youngest fully expanded leaves from representative plants were selected and used for indirect chlorophyll measurement using SPAD 502 (Minolta Corporation, Ramsey, NJ, USA) chlorophyll meter (85 and 50 days after emergence for pot and field experiment, respectively). In both experiments, plants were sampled twice: at flowering stage (50 and 40 days after sowing for pot and field experiment, respectively) and at maturity stage (final harvest), at 105 (405B) and 114 (406W) days after sowing in the pot experiment, and at 84 and 99 days after sowing (405B and 406W genotypes, respectively) in the field experiment. In the pot experiment, two repetitions, which included 50 pots (25 pots for each sampling time), were sampled by destructive harvest, whereas in the field experiment, tef plants were sampled from 1 × 1 m measured area in each sub-plot. Sampling was carried out for both genotypes at the same time.
The number of tillers per plant was recorded at maturity. The height and panicle length of the plant was measured at the end of the growing season. In the pot experiment, the roots were separated from shoots and any loose perlite was washed off with tap water and placed into a paper bag. The shoots were then rinsed several times in distilled water to avoid any contaminants and placed in a different paper bag. Both plant parts were oven-dried for 72 h at 70 °C and weighed to determine the dry matter (DM). Shoots were subsequently used for nutrient analysis.
Lodging was evaluated visually at the end of the winter pot experiment using the following scale: <20% lodged = 1, 20–40% lodged = 2, 40–60% lodged = 3, 60–80% lodged = 4, 80–100% lodged = 5. In the summer field experiment, all the plots were fully lodged by the end of the experiment and, therefore, were not evaluated.
Once harvested, the grains were separated from the straw by hand (pot experiment) or threshing machine (field experiment). The dried tef shoots were then ground in a grinding machine. Approximately 0.1 g of the ground samples and grain was digested by 2 mL sulfuric acid under 180 °C and supplemented with hydrogen peroxide [15]. The concentration of P was determined by an autoanalyzer (Lachat Instruments, Milwaukee, WI, USA) and the concentration of K was determined using an atomic absorption flame photometer (Corning 400, Corning, New York, NY, USA). Mineral uptake of shoot and grain was calculated by multiplying each mineral concentration with respective vegetative DM and grain yield.
Before treatments were applied in the field trial, soil samples from three layers (0–30 cm, 30–60 cm, and 60–90 cm-depth) were collected at five locations within the experimental plot. Soil samples were then oven-dried for 96 h at 70 °C and ground to pass through a 2 mm sieve. Soil analysis included: pH, electrical conductivity (EC), moisture content, K, P, NO3, Mg, Na and Cl in a saturated paste. Phosphorus was extracted using the Olsen bicarbonate method [16], and measured using an autoanalyzer (Lachat Instruments, Milwaukee, WI, USA). Potassium was extracted using ammonium acetate and measured using an atomic absorption flame photometer (Corning 400, Corning, New York, NY, USA). Average values of these soil properties at each depth were previously published [4]. The highest concentration of P and K was found in the 0–30 cm soil strata, at 1.1 mg/kg and 38.8 mg/kg, respectively.

2.3. Statistical Analyses

For each set of nutrient treatments (P or K levels), a two-way analysis of variance (ANOVA) was carried out for each variable to determine the effects of genotype, nutrient (P or K), and genotype X nutrient interaction using JMP 13.0 (JMP, Cary, NC, USA) software. Mean separations were performed by the Tukey–Kramer honest significant difference (HSD) test at p = 0.05.

3. Results

3.1. Effects of P and K on Different Genotypes

P and K mostly had the same effects on growth and nutrient content, regardless of the genotype. This is reflected by the fact that the interactive effect (nutrient X genotype) effect was rarely significant, with the exception of the K X genotype effect on yield in the pot experiment (p = 0.0046) (Figure 1F). Therefore, in the following paragraphs, only the main effects of genotype and nutrient concentration will be presented, except for that single case.

3.2. Pot Experiment

3.2.1. Effect of Genotypes

Overall, 406B produced greater biomass than 405W, with significantly larger shoot DM (45%, 43%) and root DM (70%, 85%), in the P and K treatments of the experiment, respectively (Figure 1A,B,D,E). Likewise, 406B was also a taller than 405W (36%, 29%) in the P and K treatments, respectively (Table 1 and Table 2). 406B flowered about 8 days later (Table 1 and Table 2) and produced 45% less tillers than 405W (Table 1 and Table 2). Differences in panicle length and SPAD between the genotypes were less than 10% (Table 1 and Table 2), and there was no difference in lodging index. In the K segment of the experiment, 406B had a significantly smaller HI (30% Table 2). Genotype had no significant effect on yield in the P segment of the experiment (Figure 1C). In the K segment of the experiment, at 40 and 80 mg L−1 K, 406B had a larger yield, but at 10 and 20 mg L−1 K, 406W had a larger yield (Figure 1F).

3.2.2. Effect of P

P fertilization had a significant effect on the DM of the shoot, root, and grain (Figure 1A–C) in both genotypes. The shoot DM was lowest at 1 mg L−1 P (17.5 g), rose sharply until 3 mg L−1 P and rose gradually thereafter until 12 mg L−1 (Figure 1A). The root DM was lowest at 1 mg L−1 P and rose sharply until 3 mg L−1 P and did not rise significantly as the P rose over 3 mg L−1 (Figure 1B). The grain DM was low at 1 mg L−1 P, rose sharply until 6 mg L−1 P, and declined drastically between 6 mg L−1 P and 12 mg L−1 P.
At flowering, P fertilization had a significant effect on shoot P concentration (Figure 2A) that rose by 400% between the 1 mg L−1 P treatment and the 12 mg L−1 P treatment. Also at maturity, P fertilization had a significant effect on shoot P concentration (Figure 2B). From 1 mg L−1 P to 6 mg L−1 P, there was no significant change in shoot P concentration, ranging between 0.04% and 0.09% (Figure 2A), but between 6 and 12 mg L−1, the P concentration rose sharply to 0.37% (Figure 2A). The grain showed a steady increase between 0.28% and 0.54% P as the P fertilization was increased from 1 to 12 mg L−1 P (Figure 2C).
Several other growth parameters were affected by P fertilization (Table 1). Tef plants grown under the lowest P treatment (1 mg L−1 P) were significantly shorter than the other treatments. The panicle length was unaffected by P fertilization, as was the SPAD. Lower P treatments took more time to reach 50% flowering than the higher P treatments, and the lodging index was higher in the high P treatments. Interestingly, the number of tillers increased as P fertilization was increased, and at the highest P treatment the HI fell dramatically.

3.2.3. Effect of K

K fertilization had no significant effect on shoot or root DM (Figure 1D,E). Grain DM was lowest at 10 mg L−1 K and rose until 40 mg L−1 K, with no significant changes between 40 mg L−1 K and 80 mg L−1 K (Figure 1F).
K fertilization had a significant effect on shoot K concentration at flowering (Figure 2D) and at maturity. At flowering, K concentration in the shoot rose steadily from 1.8% to 3.6% as the K treatment increased from 10 to 80 mg L−1. At maturity, K concentration in the shoot rose from 0.6% to 1.5% as the K treatments increased from 10 to 80 mg L−1 (Figure 2E). However, the effect of K fertilization on grain K concentration was negligible (albeit statistically significant) between the highest and lowest treatments, with the lowest being 0.48% and the highest being 0.55% (Figure 2F).
Other growth parameters were also affected by K fertilization (Table 2). The plants were 11% taller in the highest K treatment compared to the lowest K treatment, and the plants took 4 days more to reach 50% flowering in the lowest K treatment compared to the highest one. The lower K treatments had more tillers and less lodging, and the higher K treatments had a higher HI.

3.3. Field Experiment

3.3.1. Effect of Genotypes

In the field, both cultivars had relatively low yields and low HI regardless of the P and K treatments (Figure 3B,D, Table 3 and Table 4). It is worth noting that 406B had a significantly smaller grain yield than 405W. Overall, the two cultivars responded to P and K treatments in a similar way (Figure 3 and Figure 4, Table 3 and Table 4).

3.3.2. Effect of P and K

At maturity, the shoot DM was larger in the treatments that received P compared to the 0 P treatment, but there was no significant difference in shoot DM between any of the P treatments above 0 (Figure 3B). At flowering, shoot DM follows the same trend but the there are no statistically significant differences (Figure 3A). The grain yield was not significantly affected by P treatment in the field, but there was a trend towards greater grain yield as the P concentration in the irrigation water rose from 0 to 6 mg L−1, and a moderate, statistically insignificant decline thereafter (Figure 3C).
At maturity, the P concentration in the shoot was larger in the treatments that received P compared to the 0 P treatment, but there was no significant difference in the shoot P concentration between any of the P treatments above 0 (Figure 4B). At flowering, P concentration in the shoot showed no clear trend. The P treatments had no significant effect on grain P concentration (Figure 4C).
At flowering, there was no effect of K fertilization on shoot dry weight (Figure 3D). However, at maturity, the shoot dry weight and the grain yield responded favorably to K fertilization. Shoot DM increased by 27% and grain DM and 35%, respectively, between the 0 and 80 mg L−1 K treatments (Figure 3E,F). At flowering, K concentration in the shoots was not significantly affected by K fertilization (Figure 4D). At maturity, the K concentration in the shoots rose by 12% between the 0 and 80 mg L−1 K treatments (Figure 4E). K concentration in the grain was unaffected by K treatments (Figure 4D).
None of the other growth parameters (height, panicle length, days to flowering, harvest index, SPAD, or number of tillers) were significantly affected by P or K fertilization levels in the field experiment (Table 3 and Table 4).

4. Discussion

Similar to the response of tef to N fertilization detailed by Gashu et al. [4], the response to P in the pot experiment can be divided into three ranges: underfertilization, sufficient fertilization, and overfertilization. As fertilization level rose from 1 mg L−1 P to 6 mg L−1 P, the grain yield increased by over threefold, indicating that section of the graph was underfertilized (Figure 1C). At 12 mg L−1, the yield decreased dramatically, indicating overfertilization. The optimal P concentration is, therefore, between 6 and 12 mg L−1, and future experiments should include treatments within this range for a clearer picture. Interestingly, while overfertilization severely impacted the grain yield, it did not have any negative effect on the shoot or root DM (Figure 1A,B).
At flowering, P concentration in the shoot increased linearly with added P regardless of the concentration (Figure 2A). However, at maturity, between 1 mg L−1 P and 6 mg L−1 P, there was no significant increase in shoot P, while between 6 and 12 mg L−1 P the shoot P concentration increased by 400% (Figure 2B). Hawkesford et al. [17] stated that plants rarely show signs of P toxicity since they are able to downregulate P transporters involved in net root P uptake. This seems to be the case between 1 and 6 mg L−1 P, since at maturity (Figure 2B) the concentration in the shoot was not significantly increased in that range. The sharp increase in shoot P concentration between 6 and 12 mg L−1 seems to indicate that the tef plant was unable to regulate the P concentration effectively. That being said, there was little evidence of actual toxicity, since other than a reduction in grain yield, plant growth was unaffected by P overfertilization (Figure 1A,B). Given the large number of tillers in the tef that received 12 mg L−1 P compared to the other treatments (Table 1), it seems that excess P simply caused an increase in tiller production at the expense of grain production. The detrimental effects of P overfertilization were reported by Girma et al. [18] in a field trial in Oklahoma, who found that even low levels of P fertilization decreased yield.
In the field experiment, there was a trend towards higher yields at 6 mg L−1 P, which is similar to the pot experiment, but the difference in grain yield between the minimum and maximum was 33% but statistically insignificant (Figure 3B). It is noteworthy that the grain yield and HI in the field were exceedingly low. Possible reasons for this have been discussed in detail by Gashu et al. [4], including the heat and long daylight hours in mid-summer months.
It is instructive to compare the P concentration in the shoot in the field to that in the pot experiment at flowering in order to determine, at least on a basic level, what degree of P availability occurred in the field under different treatments. In the field, shoot P concentration at flowering under 0 P application, averaged across cultivars, was 0.22% (Figure 4A), while under the 12 mg L−1 P treatment it was 0.27%. Both these concentrations are close to the P concentration found in the pot experiment in the shoot at flowering time (0.24%) for the 3 mg L−1 P treatment (Figure 2A). This suggests that regardless of the treatment, the P availability in the soil was approximately 3 mg L−1 P in the root zone. Evidently, there was already some available P in the soil, somewhat below the 6 mg L−1, which gave an optimum grain yield in the pot experiment. However, adding P did not increase the availability of P in the field at flowering. Perhaps the added P was immobilized through biotic or abiotic processes.
The P concentration at maturity tells a slightly different story. At maturity, plots receiving 0 mg L−1 in the irrigation water had a significantly smaller P concentration in the shoots than those that received even low amounts of P concentrations (Figure 4B). Furthermore, the total P uptake was significantly higher in those plots that received P in the irrigation water, with no significant difference between those treatments (data not shown). Evidently, as the plant grew larger and P became more limiting, there was a significant advantage to fertigating with P.
This may explain why the grain yield at 6 mg L−1 in the field was 33% higher than the grain yield when 0 P was applied, though not a statistically significant difference. This follows the same trend as the pot experiment, in which grain yield reached a maximum at 6 mg L−1 P.
In the pot experiment, the response to K fertilization can be divided into two ranges: underfertilization and sufficient fertilization. The grain DM responded positively to K fertilization up to 40 mg L−1 K, and then plateaued (Figure 1F). Unlike the P (Figure 1A–C) and N (Gashu et al., 2020) segments of this experiment, there was no significant negative effect of K at the highest level of K fertilization, although there was a small, non-significant decline in grain DM between the 40 mg L−1 and 80 mg L−1 treatment (Figure 1F). This is similar to the response of tef grain yield to K fertilization observed by Misskire et al. [19] in a field experiment in Ethiopia. They reported yield rising and then plateauing with increased K, with a slight, non-significant drop at the highest K level. It is interesting that no decline in grain yield occurred in the higher K treatments, since Mulugeta et al. [9] found a decline in grain yield in plants with much lower shoot K concentrations (0.6–0.9%) than the concentrations we measured in the highest K fertilization treatments (Figure 2C). Furthermore, maximum yield was observed at a much higher vegetative K concentration in our experiment (1%) compared to the 0.63% K reported by Mulugeta et al. [9] as the internal K requirement of tef. It is possible that the decreased grain yield in high K treatments reported by Mulugeta et al. [9] was not caused by the increase in K itself, but rather a secondary effect caused by an imbalance of other nutrients. In our experiment, all other nutrients were provided at high rates and in available forms, so imbalances of other nutrients did not occur. This underscores the importance of conducting experiments in perlite, where minerals can be provided in a highly controlled manner.
Interestingly, in the pot experiment, there were no significant changes in shoot or root biomass across the entire K treatment range (Figure 1D,E). In the field experiment there was an increase in shoot DM as the K concentration in the fertigation water increased from 0 to 80 mg L−1, but this difference was not apparent at flowering and became apparent only at maturity (Figure 2E). The lack of response of the vegetative portions of the plant to K fertilization in the pot experiment and the late-stage response in the field experiment shows that the completion of the lifecycle of the plant and measurement of grain yield are essential for showing the effect of K fertilization on tef. The increase in straw yield due to K fertilization which we observed in the field and not in the pot experiment has been reported previously. Gebrehawariat et al. [11] reported an increase in straw with increased K fertilization in tef. Interestingly, while Gebrehawatiat et al. [11] reported that increased K caused an increase in the number of fertile tillers, we observed a decrease in tillers with increased K in the pot experiment (Table 2), and no clear effect of K on tiller number in the field experiment (Table 4).
In the pot experiment, K concentration in the shoot increased with increasing K fertilization (Figure 2D,E), but the K concentration in the grain was not responsive to K fertilization, remaining at around 0.5% regardless of cultivar or K treatment. In the field experiment, there were significant differences between the cultivars (0.46% in 406W and 0.39% in 405B) but no significant effect of K treatment on grain K concentration. Evidently, the K concentration in the shoot is flexible, but the K concentration in the grain needed to be at a certain range in order to allow grain development. This is in contrast to the findings of Gebrehawariyat et al. [11], who found that on average grain K concentration increased by 33% from the lowest K treatment (0 kg K2O ha−1) to the highest K treatment (120 kg K2O ha−1).
In the field, shoot K concentration at flowering under 0 K application, averaged across cultivars, was 2.6%, while under 80 mg L−1 K treatment it was 2.7%. Both these concentrations are close to the K concentration found in the pot experiment in shoot at flowering (2.8%) under the 40 mg L−1 K treatment (Figure 2D). This suggests that at the beginning of the field experiment, regardless of the treatment, the K availability in the soil was approximately 40 mg L−1 K in the root zone. At maturity, however, there was a 16% increase in shoot K concentration between the 0 and 80 mg L−1 K treatments, as well as a 56% increase in total K uptake between those two treatments (data not shown), indicating that fertigation with K increased K availability later on in the lifecycle.

5. Conclusions

The purpose of this work was to show the effects of P and K availability on tef growth variables. While the pot experiment showed these effects over a wide range of P and K availabilities, the field experiment demonstrated the effects over a much smaller range, evidently due to minerals available in the soil and interaction of applied nutrients with the soil. The main lessons that can be drawn from these experiments are as follows:
  • Up to 6 mg L−1 P and 40 mg L−1 K, respectively, there was a clear positive effect of P and K fertilization in the pot experiment. P and K concentration in the shoot was positively affected by P and K fertilization. The concentration of P in the grain was positively affected by P fertilization, whereas the K concentration in the grain was barely affected by K fertilization. These observations were true for both genotypes.
  • A clear negative effect of P overfertilization on grain yield was evident when perlite pots were fertigated with 12 mg L−1 P. This effect could not be seen in the field and to the best of our knowledge has not been reported before. It seems to have been caused by overinvestment in tillers and underinvestment in grain. P fertilization tended to have a positive effect on the number of tillers.
  • We observed that the benefit of K fertilization was only evident at the end of the plants’ lifecycle, which has practical ramifications for experimental design when testing the effect of K on tef.
  • No statistically significant negative effect from K overfertilization was observed in the pot or field experiment, even though the K concentration of the shoots at maturity was well above what was considered optimum in the literature.
  • While K concentration in the shoot was clearly affected by K availability, the K concentration in the grain remained more or less constant, around 0.5% regardless of K availability.
  • The response of tef to different P and K doses in the field was attenuated compared to the response in perlite, evidently because of the native ability of the soil to release and fix P and K. Fertigation recommendations in the field will need to take into account the nutrient availability in the soils.
These lessons, as well as quantitative data presented in this paper is, useful in understanding how tef responds to P and K fertilization. Together with the description of tef response to N presented in Gashu et al. (2020), this paper contributes to our understanding of how tef might be grown intensively both inside and outside Ethiopia.

Author Contributions

Conceptualization, K.G. and U.Y.; Data curation, K.G.; Formal analysis, M.H. and K.G.; Investigation, K.G., I.Z. and Y.S.; Methodology, I.Z., Y.S. and U.Y.; Project administration, U.Y.; Supervision, I.Z., Y.S. and U.Y.; Writing—original draft, M.H.; Writing—review and editing, M.H., K.G., Y.S. and U.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This project was financed by the center of fertilizer and plant nutrition (CFPN).

Acknowledgments

We would like to thank Luda, Ina Finegold, Yonatan Saroya, Jajaw Bimro, Behailu Kassahun and Fanosie Mekonen for assistance. Y.S. is the incumbent of the Haim Gvati Chair in Agriculture.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. List of treatments in winter pot experiment, including N treatments reported by Gashu et al. 2020.
Table A1. List of treatments in winter pot experiment, including N treatments reported by Gashu et al. 2020.
SegmentMineral Concentrations in Irrigation Solution (ppm)# of GenotypesRepetitions
Nitrogen (N)Phosphorus (P)Potassium (K)
Nitrogen Segment1064025
2064025
8064025
12064025
Phosphorus Segment4014025
4034025
40124025
Potassium Segment 4061025
4062025
4068025
Shared control group4064025
Table A2. List of treatments in summer field experiment, including N treatments reported by Gashu et al. 2020.
Table A2. List of treatments in summer field experiment, including N treatments reported by Gashu et al. 2020.
SegmentMineral Concentrations in Irrigation Solution (ppm)# of GenotypesRepetitions
Nitrogen (N)Phosphorus (P)Potassium (K)
Nitrogen Segment064025
3064025
12064025
Phosphorus Segment6004025
6034025
60124025
Potassium Segment606025
6068025
Shared Control group6064025

References

  1. Lee, H.; Teff, A. Rising Global Crop: Current Status of Teff Production and Value Chain. Open Agric. J. 2018, 12, 185–193. [Google Scholar] [CrossRef]
  2. Gebru, Y.A.; Sbhatu, D.B.; Kim, K.P. Nutritional Composition and Health Benefits of Teff (Eragrostis tef (Zucc.) Trotter). J. Food Qual. 2020, 2020. [Google Scholar] [CrossRef]
  3. Yihun, Y.M. Agricultural Water Productivity Optimization for Irrigated Teff (Eragrostic tef) in Water Scarce Semi-Arid Region of Ethiopia; Springer: New York, NY, USA, 2015; ISBN 9781138027664. [Google Scholar]
  4. Gashu, K.; Halpern, M.; Zipori, I.; Bustan, A.; Saranga, Y.; Yermiyahu, U. Tef (Eragrostis tef Trotter) Responses to Nitrogen Fertigation under Semi-Arid Mediterranean Climate. Agronomy 2020, 10, 1870. [Google Scholar] [CrossRef]
  5. Reda, G.T. Verification and Demonstration of Soil Test Based Phosphorus Fertilizer Recommendation Rate on Yield of Teff (Eragrostis Tef (Zucc) Trotter) in Vertisols of Northern Ethiopia. J. Nat. Sci. Res. 2016, 6, 51–55. [Google Scholar]
  6. Dubale, T.K. Soil Test and Tissue Analysis Based Nitrogen and Phosphorus Fertilizer Applications on Teff [Eragrostis tef (Zucc.) Trotter] Cultivars in North Shewa Central Highlands, Ethiopia. Adv. Crop Sci. Technol. 2019, 7, 1–6. [Google Scholar] [CrossRef]
  7. Dereje, G.; Alemu, D.; Adisu, T.; Anbessa, B. Response of Yield and Yield Components of Tef [Eragrostis Tef (Zucc.) Trotter] to Optimum Rates of Nitrogen and Phosphorus Fertilizer Rate Application in Assosa Zone, Benishangul Gumuz Region. Ethiop. J. Agric. Sci. 2018, 28, 81–94. [Google Scholar]
  8. Mulugeta, D.; Mamo, T.; Sokolowski, E.; Nachmansohn, E. Potash Fertilization of Teff and Wheat in the Highlands of Ethiopia. Electron. Int. Fertil. Corresp. 2017, 48, 12–19. [Google Scholar]
  9. Mulugeta, D.; Mamo, T.; Beyene, S.; Kidanu, S. Effect of potassium levels on teff (Eragrostis tef (zucc.) trotter) growth and yield in central highland vertisols of Ethiopia. Eurasian J. Soil Sci. 2020, 9, 105–118. [Google Scholar] [CrossRef] [Green Version]
  10. Ethiopian Agricultural Transformation Agency. Ministry of Agriculture Soil Fertility Status and Fertilizer Recommendation Atlas for Tigray Regional State; Ethiopian Agricultural Transformation Agency: Addis Ababa, Ethiopia, 2014.
  11. Gebrehawariyat, F.M.; Haile, W.; Mamo, T.; Zipori, I.; Sokolowski, E. Response of teff [Eragrostis tef (Zucc.) Trotter] to potassium fertilizer application in four districts of North Shewa, Ethiopia. Electron. Int. Fertil. Corresp. 2018, 53, 3–15. [Google Scholar]
  12. Mulugeta, D.; Tekalign, M.; Sheleme, B.; Selamyihun, K. Potassium critical level in soil for Teff (Eragrostis tef (Zucc.) Trotter) grown in the central highland soils of Ethiopia. SN Appl. Sci. 2019, 1, 1–9. [Google Scholar] [CrossRef] [Green Version]
  13. Ben-Zeev, S.; Bimro, J.; Barak, V.; Saranga, Y. Phenotypic diversity and heritability in Eragrostis tef under irrigated Mediterranean conditions. Isr. J. Plant Sci. 2019, 65, 222–231. [Google Scholar] [CrossRef]
  14. Yihun, Y.M.; Haile, A.M.; Schultz, B.; Erkossa, T. Crop Water Productivity of Irrigated Teff in a Water Stressed Region. Water Resour. Manag. 2013, 27, 3115–3125. [Google Scholar] [CrossRef]
  15. Snell, F.D.; Snell, C.T. Colorimetric Methods of Analysis, 3rd ed.; D. Van Nostran Co Inc.: New York, NY, USA, 1949. [Google Scholar]
  16. Pierzynski, G.; Sharpley, A.N. Methods of Phosphorus Analysis for Soils, Sediments, Residuals, and Waters: Introduction. In Methods of Phosphorus Analysis for Soils, Sediments, and Waters; North Carolina State University: Raleigh, NC, USA, 2000; p. 110. [Google Scholar]
  17. Hawkesford, M.J.; Horst, W.; Kichey, T.; Lambers, H.; Schjoerring, J.K.; Moller, I.S.; White, P. Functions of macronutrients. In Mineral Nutrition of Higher Plants; Marchner, P., Ed.; Elsevier: Amsterdam, The Netherlands, 2012. [Google Scholar]
  18. Girma, K.; Reinert, M.; Ali, M.S.; Sutradhar, A.; Mosali, J. Nitrogen and Phosphorus Requirements of Teff Grown Under Dryland Production System. Crop Manag. 2012, 11, 1–14. [Google Scholar] [CrossRef]
  19. Misskire, Y.; Mamo, T.; Taddesse, A.M.; Yermiyahu, U. The effect of potassium on yield, nutrient uptake and efficiency of teff (Eragrostis tef Zucc. Trotter) on vertisols of North Western Ethiopian Highlands. J. Plant Nutr. 2019, 42, 307–322. [Google Scholar] [CrossRef]
Figure 1. Shoot (A), root (B) and grain (C) dry matter at maturity in the winter pot experiment by P fertigation level, and Shoot (D), root (E) and grain (F) dry matter at maturity in the winter pot experiment by K fertigation level. Dashed line represents 405B genotype, and solid line represents 406W genotype. Different letters represent significant (<0.05) differences for main effects. Lower case letters adjacent to the trend line represent the results of the Tukey test for nutrient concentration averaged across genotypes, and capital letters adjacent to the legend represent the results of the Student-t test for the two different genotypes.
Figure 1. Shoot (A), root (B) and grain (C) dry matter at maturity in the winter pot experiment by P fertigation level, and Shoot (D), root (E) and grain (F) dry matter at maturity in the winter pot experiment by K fertigation level. Dashed line represents 405B genotype, and solid line represents 406W genotype. Different letters represent significant (<0.05) differences for main effects. Lower case letters adjacent to the trend line represent the results of the Tukey test for nutrient concentration averaged across genotypes, and capital letters adjacent to the legend represent the results of the Student-t test for the two different genotypes.
Agronomy 11 01588 g001
Figure 2. Shoot P concentration at flowering (A), shoot P concentration at maturity (B), and grain P concentration (C) in the winter pot experiment by P fertigation level. Shoot K concentration at flowering (D), shoot K concentration at maturity (E) and grain K concentration (F) in the winter pot experiment by K fertigation level. Dashed line represents 405B genotype, and solid line represents 406W genotype. Different letters represent significant (<0.05) differences for main effects. Lower case letters adjacent to the trend line represent the results of the Tukey test for nutrient concentration averaged across both genotypes, and capital letters adjacent to the legend represent the results of the Student-t test for the two different genotypes.
Figure 2. Shoot P concentration at flowering (A), shoot P concentration at maturity (B), and grain P concentration (C) in the winter pot experiment by P fertigation level. Shoot K concentration at flowering (D), shoot K concentration at maturity (E) and grain K concentration (F) in the winter pot experiment by K fertigation level. Dashed line represents 405B genotype, and solid line represents 406W genotype. Different letters represent significant (<0.05) differences for main effects. Lower case letters adjacent to the trend line represent the results of the Tukey test for nutrient concentration averaged across both genotypes, and capital letters adjacent to the legend represent the results of the Student-t test for the two different genotypes.
Agronomy 11 01588 g002
Figure 3. Shoot dry matter at flowering (A), Shoot dry matter at maturity (B) and grain dry matter (C) in the summer field experiment by P fertigation level, and Shoot (D), root (E) and grain (F) dry matter in the summer field experiment by K fertigation level. Dashed line represents 405B genotype, and solid line represents 406W genotype. Different letters represent significant (<0.05) differences for main effects. Lower case letters adjacent to the trend line represent the results of the Tukey test for N concentration averaged across genotypes, and capital letters adjacent to the legend represent the results of the Student-t test for the two different genotypes.
Figure 3. Shoot dry matter at flowering (A), Shoot dry matter at maturity (B) and grain dry matter (C) in the summer field experiment by P fertigation level, and Shoot (D), root (E) and grain (F) dry matter in the summer field experiment by K fertigation level. Dashed line represents 405B genotype, and solid line represents 406W genotype. Different letters represent significant (<0.05) differences for main effects. Lower case letters adjacent to the trend line represent the results of the Tukey test for N concentration averaged across genotypes, and capital letters adjacent to the legend represent the results of the Student-t test for the two different genotypes.
Agronomy 11 01588 g003
Figure 4. Shoot P concentration at flowering (A), shoot P concentration at maturity (B), and grain P concentration (C) in the summer field experiment by P fertigation level. Shoot K concentration at flowering (D), shoot K concentration at maturity (E) and grain K concentration (F) in the summer field experiment by K fertigation level. Dashed line represents 405B genotype, and solid line represents 406W genotype. Different letters represent significant (<0.05) differences for main effects. Lower case letters adjacent to the trend line represent the results of the Tukey test for N concentration averaged across genotypes, and capital letters adjacent to the legend represent the results of the Student-t test for the two different genotypes.
Figure 4. Shoot P concentration at flowering (A), shoot P concentration at maturity (B), and grain P concentration (C) in the summer field experiment by P fertigation level. Shoot K concentration at flowering (D), shoot K concentration at maturity (E) and grain K concentration (F) in the summer field experiment by K fertigation level. Dashed line represents 405B genotype, and solid line represents 406W genotype. Different letters represent significant (<0.05) differences for main effects. Lower case letters adjacent to the trend line represent the results of the Tukey test for N concentration averaged across genotypes, and capital letters adjacent to the legend represent the results of the Student-t test for the two different genotypes.
Agronomy 11 01588 g004
Table 1. Effect of P fertigation levels on growth parameters in the winter pot experiment. Different letters represent significant (<0.05) differences for main effects. Upper case letters next to the values in the 3rd row represent the results of the Tukey test for nutrient concentration averaged across genotypes, and lower case letters next to the values in the 4th row represent the results of the Student-t test for the two different genotypes.
Table 1. Effect of P fertigation levels on growth parameters in the winter pot experiment. Different letters represent significant (<0.05) differences for main effects. Upper case letters next to the values in the 3rd row represent the results of the Tukey test for nutrient concentration averaged across genotypes, and lower case letters next to the values in the 4th row represent the results of the Student-t test for the two different genotypes.
Phosphorus Treatment (mg L−1)Height (cm)Panicle Length (cm)Days to FloweringHarvest IndexSPADNumber of Tillers per PlantLodging Index
Genotype 405B135.532.850.60.3041.65.31.2
341.639.447.60.2440.010.71.1
645.940.647.00.2940.610.21.9
1245.443.245.00.0840.016.61.7
Genotype 406W 1 50.532.857.00.2939.12.21.1
359.136.256.40.1738.86.31.4
660.037.853.80.2838.05.51.9
1258.937.251.60.0337.89.41.8
Tukey for P Treatments143.0 B32.8 A53.8 A0.29 A40.3 A3.8 C1.1 B
350.3 A37.8 A52.0 AB0.21 A39.9 A8.5 B1.2 B
652.9 A39.2 A50.1 BC0.28 A39.3 A7.8 B1.9 A
1252.1 A40.4 A48.3 C0.05 B38.8 A13.0 A1.7 A
Tukey for Genotype405B42.1 b39.0 a47.5 b0.23 a40.8 a10.7 a1.5 a
406W57.1 a36.0 b55.7 a0.19 a38.4 b5.8 b1.5 a
P x GenotypeN.S.N.S.N.S.N.S.N.S.N.S.N.S.N.S.
N.S. = Non-significant.
Table 2. Effect of K fertigation levels on growth parameters in the winter pot experiment. Different letters represent significant (<0.05) differences for main effects. Upper case letters next to the values in the 3rd row represent the results of the Tukey test for nutrient concentration averaged across genotypes, and lower case letters next to the values in the 4th row represent the results of the Student-t test for the two different genotypes.
Table 2. Effect of K fertigation levels on growth parameters in the winter pot experiment. Different letters represent significant (<0.05) differences for main effects. Upper case letters next to the values in the 3rd row represent the results of the Tukey test for nutrient concentration averaged across genotypes, and lower case letters next to the values in the 4th row represent the results of the Student-t test for the two different genotypes.
Potassium Treatment (mg L−1)Height (cm)Panicle Length (cm)Days to FloweringHarvest IndexSPADNumber of Tillers per PlantLodging Index
Genotype 405B 10 42.137.049.20.1841.014.21.8
20 45.640.847.60.2041.513.51.9
40 45.840.647.00.2940.610.21.9
80 46.737.245.60.2440.810.02.5
Genotype 406W 10 54.038.057.80.0639.27.31.6
20 56.437.256.40.1138.98.01.7
40 56.737.853.80.2838.05.51.9
80 61.840.253.80.2138.05.82.7
Tukey for K Treatments1048.0 B37.5 A53.0 A0.12 C40.1 A10.8 A1.7 B
2051.0AB39.0 A52.0 AB0.16 BC40.2 A10.7 A1.8 B
4052.9 A39.2 A50.4 B0.28 A39.3 A7.8 B1.9 AB
8054.2 A40.9 A49.7 B0.22 AB39.4 A7.8 B2.6 A
Tukey for Genotype 405B45.1 b40.0 a47.3 b0.23 a41.0 a12.0 a2.0 a
406W58.0 a38.0 a55.4 a0.16 b38.6 b6.6 b2.0 a
K x GenotypeN.S.N.S.N.S.N.S.N.S.N.S.N.S.N.S.
N.S. = Non-significant.
Table 3. Effect of P fertigation levels on growth parameters in the summer field experiment. Different letters represent significant (<0.05) differences for main effects. Upper case letters next to the values in the 3rd row represent the results of the Tukey test for nutrient concentration averaged across genotypes, and lower case letters next to the values in the 4th row represent the results of the Student-t test for the two different genotypes.
Table 3. Effect of P fertigation levels on growth parameters in the summer field experiment. Different letters represent significant (<0.05) differences for main effects. Upper case letters next to the values in the 3rd row represent the results of the Tukey test for nutrient concentration averaged across genotypes, and lower case letters next to the values in the 4th row represent the results of the Student-t test for the two different genotypes.
Phosphorus Treatment (mg L−1)Height (cm)Panicle Length (cm)Days to Flowering (50%)Harvest IndexSPADNumber of Tillers per Plant
Genotype 405B195.632.240.40.0935.414.6
395.931.537.80.0833.517.2
6100.733.137.40.1037.413.8
12102.332.037.00.0734.919.2
Genotype 406W194.724.854.20.0436.88.8
397.330.754.20.0334.912.8
694.634.756.40.0537.411.2
1297.734.656.40.0337.010.8
Tukey for P Treatments195.6 A28.5 A47.3 A0.06 A36.1 A11.7 A
397.3 A31.1 A46.0 A0.06 A34.2 A15.0 A
694.6 A33.9 A46.9 A0.08 A37.3 A12.5 A
1297.7 A33.3 A46.7 A0.05 A37.015.0 A
Tukey for Genotype 405B98.9 a32.2 a38.2 b0.09 a35.3 a16.2 a
406W96.1 a31.2 a55.3 a0.04 b36.5 a10.9 a
P x GenotypeN.S.N.S.N.S.N.S.N.S.N.S.N.S.
N.S. = Non-significant.
Table 4. Effect of K fertigation levels on growth parameters in the summer field experiment. Different letters represent significant (<0.05) differences for main effects. Upper case letters next to the values in the 3rd row represent the results of the Tukey test for nutrient concentration averaged across genotypes, and lower case letters next to the values in the 4th row represent the results of the Student-t test for the two different genotypes.
Table 4. Effect of K fertigation levels on growth parameters in the summer field experiment. Different letters represent significant (<0.05) differences for main effects. Upper case letters next to the values in the 3rd row represent the results of the Tukey test for nutrient concentration averaged across genotypes, and lower case letters next to the values in the 4th row represent the results of the Student-t test for the two different genotypes.
Potassium Treatment (mg L−1)Height (cm)Panicle Length (cm)Days to Flowering (50%)Harvest IndexSPADNumber of Tillers per Plant
Genotype 405B0101.031.837.80.1135.017.0
40100.733.137.40.1034.413.8
8092.626.638.60.1137.220.8
Genotype 406W 0 108.032.560.00.0537.111.0
4094.634.656.40.0537.311.2
8097.232.755.00.0637.310.2
Tukey for K Treatments0104.5 A32.1 A48.9 A0.08 A36.0 A14.0 A
4097.6 A33.9 A46.9 A0.08 A37.4 A12.5A
8094.9 A29.7 A46.8 A0.08 A37.3 A15.5 A
Tukey for Genotype 405B98.1 a30.5 a37.9 b0.11 a36.5 a17.2 a
406W99.9 a33.3 a57.1 a0.05 b37.2 a10.8 b
K x GenotypeN.S.N.S.N.S.N.S.N.S.N.S.N.S.
N.S. = Non-significant.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Halpern, M.; Gashu, K.; Zipori, I.; Saranga, Y.; Yermiyahu, U. Tef (Eragrostis tef) Responses to Phosphorus and Potassium Fertigation under Semi-Arid Mediterranean Climate. Agronomy 2021, 11, 1588. https://doi.org/10.3390/agronomy11081588

AMA Style

Halpern M, Gashu K, Zipori I, Saranga Y, Yermiyahu U. Tef (Eragrostis tef) Responses to Phosphorus and Potassium Fertigation under Semi-Arid Mediterranean Climate. Agronomy. 2021; 11(8):1588. https://doi.org/10.3390/agronomy11081588

Chicago/Turabian Style

Halpern, Moshe, Kelem Gashu, Isaac Zipori, Yehoshua Saranga, and Uri Yermiyahu. 2021. "Tef (Eragrostis tef) Responses to Phosphorus and Potassium Fertigation under Semi-Arid Mediterranean Climate" Agronomy 11, no. 8: 1588. https://doi.org/10.3390/agronomy11081588

APA Style

Halpern, M., Gashu, K., Zipori, I., Saranga, Y., & Yermiyahu, U. (2021). Tef (Eragrostis tef) Responses to Phosphorus and Potassium Fertigation under Semi-Arid Mediterranean Climate. Agronomy, 11(8), 1588. https://doi.org/10.3390/agronomy11081588

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop