# A Crop Simulation Model for Tef (Eragrostis tef (Zucc.) Trotter)

Kirsten Paff, Senthold Asseng

► To cite this version:

Kirsten Paff, Senthold Asseng. A Crop Simulation Model for Tef (Eragrostis tef (Zucc.) Trotter). Agronomy, 2019, 9 (12), 10.3390/agronomy9120817. [https://doi.org/10.3390/agronomy9120817](https://doi.org/10.3390/agronomy9120817)

HAL Id: hal-02620840

https://hal.inrae.fr/hal-02620840

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# A Crop Simulation Model for Tef (Eragrostis tef (Zucc.) Trotter)

Kirsten Paff *,† and Senthold Asseng

Department of Agricultural and Biological Engineering, University of Florida; Gainesville, FL 32603, USA;
sasseng@ufl.edu

* Correspondence: kirsten.paff@inra.fr

† Currently at the French National Institute of Agricultural Research (INRA), UMR System,
34000 Montpellier, France.

Received: 29 October 2019; Accepted: 26 November 2019; Published: 28 November 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Tef is an Ethiopian staple grain that provides both food security and income for smallholders. As tef is nutritious and gluten free, it is also gaining popularity as a health food. A tef model was calibrated based on the Decision Support System for Agrotechnology Transfer's (DSSAT) NWheat model and included parameter changes in phenology, photoperiod response, radiation use efficiency, and transpiration efficiency for both standard and elevated atmospheric CO₂, based on published literature for tef and other C4 species. The new DSSAT-Tef model was compared with tef field experiments. DSSAT-Tef accurately simulated phenology and responded to changes in N supply and irrigation, but overestimated growth and occasionally yields. Simulation-observation comparisons resulted in an RMSE of 2.5 days for anthesis, 4.4 days for maturity, 2624 kg/ha (49.6%) for biomass, and 475 kg/ha (41.0%) for grain yield. Less data were available for N uptake, and the model simulated crop N uptake with an RMSE of 45 kg N/ha (46.2%) and 15 kg N/ha (37.3%) for grain N. While more data from contrasting environments are needed for further model testing, DSSAT-Tef can be used to assess the performance of crop management strategies, the suitability of tef for cultivation across growing environments, and food security.

**Keywords:** crop model; tef; DSSAT; Ethiopia

## 1. Introduction

Tef, or *Eragrostis tef* (Zucc.) Trotter, is an Ethiopian staple grain [1]. A strong preference for tef grain in Ethiopia means that smallholder farmers can demand a higher price for tef than for other grains, making tef an important source of income [2]. Tef's high calcium and iron content [3,4], as well as its lack of gluten [5] have resulted in a growing demand for tef as a health food in industrialized nations. In both Ethiopia and abroad, tef straw is valued as a high quality, low input, warm season fodder [6,7]. Crop models are computer programs that simulate crop growth and yields, under varying environmental conditions and management practices [8]. Crop models can be used to virtually test management practices and crop traits using a fraction of the time and money needed to run a field trial [9]. Crop models are particularly important for assessing the effects of climate change on agricultural production and for evaluating food security [10,11].

Two published crop models currently exist for tef. The Food and Agriculture Organization (FAO)-AEZ crop growth simulation model estimates yields using a three step process [12]. First, the model calculates the radiation limited yield, then the water limited yield, and finally, the yield affected by soil and management limitations. Each yield is calculated based on the yield of the previous step [12]. The FAO-AEZ model has certain limitations. As growing season length is treated as an input, the model cannot simulate the effects of temperature, or day length, on phenology [12]. The
FAO-AEZ model also only produces yields at maturity, and not the dynamic daily changes in yields, which restricts its applications [12]. The second tef model, the FAO's AquaCrop model, simulates biomass production using the daily crop transpiration and the normalized crop water productivity [13]. A crop specific harvest index is then used to calculate the grain yield based on the simulated biomass production [14]. The AquaCrop tef model was designed for water limited growing scenarios and therefore does not directly take soil nutrients and fertilizer management practices into account [13]. This limits its application in Ethiopia, where soil degradation [15] and low levels of fertilizer use [16] cause nutrient limited growing conditions. In addition, the AquaCrop tef model was not developed to take differences in daylength into account [1, 13, 17–19], limiting its application to low latitudes.

The Decision Support System for Agrotechnology Transfer (DSSAT) is a globally used crop modeling platform that works on a daily time step and accounts for the effects of nitrogen, CO₂, and daylength on crop growth and development [20,21]. The goal of this paper was to calibrate an existing, well-tested wheat crop module for tef in the DSSAT platform to create a DSSAT-Tef model using published field data for tef and other C4 crops, such as sorghum.

## 2. Materials and Methods

### 2.1. Model Choice

The tef crop model was created using the DSSAT-NWheat model as a basis. The NWheat model was selected, instead of a C4 crop model like DSSAT-Sorghum, as a starting point, because of its robust model structure, which has been successfully applied to many different growing conditions around the world [22]. The parameters and structures of DSSAT-NWheat were preserved unless the literature suggested a significant difference between tef and wheat. The major differences between tef and wheat are that tef is a short-day [6] C4 crop [23], while wheat is a long-day [24] C3 crop [25]. In order to adapt the wheat model for tef, changes were made to the phenology, photoperiod response, radiation use efficiency, and transpiration efficiency for both standard and elevated atmospheric CO₂ based on published research for tef and other C4 crops, such as sorghum, maize, and millet. The Agricultural Production Systems Simulator (APSIM)-Sorghum model [26], in particular, was referenced when there were gaps in the tef literature, as sorghum is a C4 crop, like tef, and the parameters were often already in the correct format for the DSSAT-NWheat model. This paper uses model specific nomenclature, but the DSSAT parameter names are explained so that a reader with a general understanding on crop modeling can understand them. For more detailed descriptions of the general DSSAT model structure and of the DSSAT-NWheat model, the reader should refer to Jones et al. [21] and Kassie et al. [22].

### 2.2. Field Locations

The DSSAT-Tef simulations were compared to observed field data from 9 Ethiopian locations and one location in the United States (Table 1). The Ude 2 location from Tulema et al. [27] was discarded due to reported high weed pressure, which the model could not simulate as it only accounts for abiotic stresses. Many tef field experiments from the literature could not be used for model calibration, because they lacked sufficient data on the environmental conditions or management practices. Due to limited field data, there were no separate calibration and validation datasets. Table 2 lists the soil characteristics for the experiments used. In cases where the drainage parameters were not given, they were calculated in DSSAT's SBuild program [20] using the soil clay, silt, and sand percentages. All soils were standardized to a depth of 90 cm with the same layer depths for each soil. In cases where the properties of the soil layer were unknown, it was assumed that they were the same as the layer above. Root growth for all locations, except for the Fallon, NV, and Ofla District locations, was limited to the top 30 cm of the soil, by setting the root growth factor to 0 below 30 cm. This was based on Araya et al.'s [1] findings that the root depth for early maturing tef varieties was generally less than 30 cm. A 90 cm maximum root depth was assumed for the Fallon, NV, and Ofla District locations due to the later maturity [28,29], which is associated with deeper roots for tef [30]. The high
yields at these locations could not have been achieved with shallow roots [29,31]. Table 3 specifies the management practices used for the simulations. The DSSAT-Tef model uses post-emergence plant density (plants/m²), rather than the more commonly reported seeding rate (kg seed planted/ha), so the plant density had to be assumed for most sites (Table 3).

A combination of satellite and station weather data was used. The Climate Hazards Group InfraRed Precipitation with Station (CHIRPS) weather dataset [32] was used for the Jamma District location and for the Ude, Kajima, and both Gare Arera locations. A combination of CHIRPS temperature data and observed rainfall data [13] was used for the Mekelle and Ilala locations. The NASA Prediction of Worldwide Energy Resource (POWER) dataset [33] was used for the Ofla District location. The Fallon Experimental, NV, USA. Weather Station was used for the Fallon, NV, location, though missing data were filled in using data from the Fallon NAAS, NV, USA. Weather Station [34].

<table><caption>**Table 1.** Source of experimental data for tef model calibration.</caption><thead><tr><th>Location Name</th><th>Latitude/Longitude</th><th>Elevation (m a.s.l.)</th><th>Soil Type</th><th>Avg. Yearly Rainfall (mm)</th><th>Year</th><th>Source</th></tr></thead><tbody><tr><td>Mekelle 1, Ethiopia</td><td>13°3' N/39°6' E</td><td>2130 <sup>1</sup></td><td>Cambisol</td><td>600</td><td>2008</td><td>[13]</td></tr><tr><td>Mekelle 2, Ethiopia</td><td>13°3' N/39°6' E</td><td>2130 <sup>1</sup></td><td>Cambisol</td><td>600</td><td>2009</td><td>[13]</td></tr><tr><td>Ilala, Ethiopia</td><td>13°4' N/39°4' E</td><td>NA</td><td>Vertisol</td><td>650</td><td>2008</td><td>[13]</td></tr><tr><td>Fallon, NV, USA</td><td>39°29'20" N/118°49'41" W <sup>2</sup></td><td>1215 <sup>2</sup></td><td>Dia Loam Soil</td><td>127 <sup>2</sup></td><td>2009</td><td>[29]</td></tr><tr><td>Jamma District, Ethiopia</td><td>10°27' N/39°15' E</td><td>2630</td><td>Vertisol</td><td>868.2</td><td>2011</td><td>[16]</td></tr><tr><td>Ofla District, Ethiopia</td><td>12°31'58" N/39°30'13" E</td><td>2450</td><td>Vertisol</td><td>980.5</td><td>2011</td><td>[31]</td></tr><tr><td>Ude 1, Ethiopia</td><td>8°48' N/39°38' E</td><td>1800</td><td>Vertisol</td><td>775 <sup>3</sup></td><td>1997</td><td>[27]</td></tr><tr><td>Kajima, Ethiopia</td><td>8°48' N/39°38' E</td><td>1920</td><td>Andisol</td><td>775 <sup>3</sup></td><td>1997</td><td>[27]</td></tr><tr><td>Gare Arera 1, Ethiopia</td><td>09°03' N/38°30' E</td><td>NA</td><td>Nitisol</td><td>1046 <sup>4</sup></td><td>2002</td><td>[35]</td></tr><tr><td>Gare Arera 2, Ethiopia</td><td>09°03' N/38°30' E</td><td>NA</td><td>Vertisol</td><td>1046 <sup>4</sup></td><td>2002</td><td>[35]</td></tr></tbody></table>

## 2.3. *Model Initialization*

Not all locations reported the same types of data, so the observed sample size ranged from 9 to 38. As field data for tef were limited and published field experiments rarely gave full details on the environmental conditions and management practices, assumptions had to be made to fill in some of the gaps. Initial conditions were based on observed measurements. If these were not available, the initial conditions were calibrated by individually changing soil water and N to match the simulated yields with the observed yield of the lowest input (N or irrigation) treatment. The same initial conditions calibrated for the lowest input treatment were then applied to all other treatments within this experiment (Table 4). Simulations started one month before planting to allow surface soil water to reach the equilibrium it would have in the actual field.

<footnote n="1">¹ From Araya et al. [1]. ² From a Dia loam, 0 to 1 percent slopes group from the Fallon-Fernley Area, Nevada, parts of Churchill, Lyon, Storey, and Washoe Counties' soil survey [36]. ³ Annual rainfall in 1997. ⁴ Annual rainfall in 2002.</footnote>
<table><caption>**Table 2.** Soil characteristics for a Cambisol soil at Mekelle in 2008 and 2009, a Vertisol soil at Ilala and a Dia Loam soil at Fallon, Vertisol soils in the Jamma and Ofla districts, a Vertisol soil at Ude 1, an Andisol soil at Kajima, and a Nitisol and a Vertisol soil at Gare Arera. Lower limit for plant available soil water (LL), drained upper limit (DUL), saturation (SAT), and soil organic carbon (OC). Data after Araya et al. [13], Davison and Laca [29], Getu [16], Okubay [31], Tulema et al. [27], and Tulema et al. [35].</caption><thead><tr><th colspan="8">Mekelle 1, Ethiopia</th><th colspan="8">Mekelle 2, Ethiopia</th></tr><tr><th colspan="8">Cambisol</th><th colspan="8">Cambisol</th></tr><tr><th>Depth (m)</th><th>LL <sup>1</sup></th><th>DUL <sup>2</sup></th><th>SAT <sup>3</sup></th><th>OC <sup>4</sup></th><th>SRGF <sup>5</sup></th><th>Depth (m)</th><th>LL <sup>1</sup></th><th>DUL <sup>2</sup></th><th>SAT <sup>3</sup></th><th>OC <sup>4</sup></th><th>SRGF <sup>5</sup></th><th>Depth (m)</th><th>LL <sup>1</sup></th><th>DUL <sup>2</sup></th><th>SAT <sup>3</sup></th><th>OC <sup>4</sup></th><th>SRGF <sup>5</sup></th></tr></thead><tbody><tr><td>0.05</td><td>0.14</td><td>0.27</td><td>0.46</td><td>0.72 <sup>7</sup></td><td>1</td><td>0.05</td><td>0.14</td><td>0.27</td><td>0.46</td><td>0.72 <sup>7</sup></td><td>1</td><td>0.19</td><td>0.38</td><td>0.54</td><td>1.02 <sup>9</sup></td><td>1</td><td>0.05</td><td>0.16 <sup>8</sup></td><td>0.17 <sup>10</sup></td><td>0.47 <sup>8</sup></td><td>1.16 <sup>10</sup></td><td>1</td></tr><tr><td>0.10</td><td>0.14</td><td>0.27</td><td>0.46</td><td>0.72 <sup>7</sup></td><td>1</td><td>0.10</td><td>0.14</td><td>0.27</td><td>0.46</td><td>0.72 <sup>7</sup></td><td>1</td><td>0.19</td><td>0.38</td><td>0.54</td><td>1.02 <sup>9</sup></td><td>1</td><td>0.10</td><td>0.16 <sup>8</sup></td><td>0.17 <sup>10</sup></td><td>0.47 <sup>8</sup></td><td>1.16 <sup>10</sup></td><td>1</td></tr><tr><td>0.20</td><td>0.14</td><td>0.27</td><td>0.50</td><td>0.60 <sup>7</sup></td><td>1</td><td>0.20</td><td>0.14</td><td>0.27</td><td>0.50</td><td>0.60 <sup>7</sup></td><td>1</td><td>0.19</td><td>0.38</td><td>0.54</td><td>1.02 <sup>9</sup></td><td>1</td><td>0.20</td><td>0.14 <sup>8</sup></td><td>0.17 <sup>10</sup></td><td>0.45 <sup>8</sup></td><td>0.66 <sup>10</sup></td><td>1</td></tr><tr><td>0.30</td><td>0.22</td><td>0.42</td><td>0.52</td><td>0.60 <sup>7</sup></td><td>0.61</td><td>0.30</td><td>0.22</td><td>0.42</td><td>0.52</td><td>0.60 <sup>7</sup></td><td>0.61</td><td>0.19</td><td>0.38</td><td>0.54</td><td>1.02 <sup>9</sup></td><td>0.61</td><td>0.30</td><td>0.14 <sup>8</sup></td><td>0.17 <sup>10</sup></td><td>0.45 <sup>8</sup></td><td>0.66 <sup>10</sup></td><td>0.61</td></tr><tr><td>0.40</td><td>0.22 <sup>6</sup></td><td>0.42 <sup>6</sup></td><td>0.52 <sup>6</sup></td><td>0.60 <sup>6</sup></td><td>0</td><td>0.40</td><td>0.22 <sup>6</sup></td><td>0.42 <sup>6</sup></td><td>0.52 <sup>6</sup></td><td>0.60 <sup>6</sup></td><td>0</td><td>0.19 <sup>6</sup></td><td>0.38 <sup>6</sup></td><td>0.54 <sup>6</sup></td><td>0.90 <sup>9</sup></td><td>0</td><td>0.40</td><td>0.14 <sup>8</sup></td><td>0.17 <sup>10</sup></td><td>0.45 <sup>8</sup></td><td>0.66 <sup>10</sup></td><td>0.50</td></tr><tr><td>0.50</td><td>0.22 <sup>6</sup></td><td>0.42 <sup>6</sup></td><td>0.52 <sup>6</sup></td><td>0.60 <sup>6</sup></td><td>0</td><td>0.50</td><td>0.22 <sup>6</sup></td><td>0.42 <sup>6</sup></td><td>0.52 <sup>6</sup></td><td>0.60 <sup>6</sup></td><td>0</td><td>0.19 <sup>6</sup></td><td>0.38 <sup>6</sup></td><td>0.54 <sup>6</sup></td><td>0.90 <sup>9</sup></td><td>0</td><td>0.50</td><td>0.14 <sup>8</sup></td><td>0.17 <sup>10</sup></td><td>0.44 <sup>8</sup></td><td>0.44 <sup>10</sup></td><td>0.41</td></tr><tr><td>0.60</td><td>0.22 <sup>6</sup></td><td>0.42 <sup>6</sup></td><td>0.52 <sup>6</sup></td><td>0.60 <sup>6</sup></td><td>0</td><td>0.60</td><td>0.22 <sup>6</sup></td><td>0.42 <sup>6</sup></td><td>0.52 <sup>6</sup></td><td>0.60 <sup>6</sup></td><td>0</td><td>0.19 <sup>6</sup></td><td>0.38 <sup>6</sup></td><td>0.54 <sup>6</sup></td><td>0.90 <sup>9</sup></td><td>0</td><td>0.60</td><td>0.10 <sup>8</sup></td><td>0.12 <sup>10</sup></td><td>0.44 <sup>8</sup></td><td>0.44 <sup>10</sup></td><td>0.33</td></tr><tr><td>0.70</td><td>0.22 <sup>6</sup></td><td>0.42 <sup>6</sup></td><td>0.52 <sup>6</sup></td><td>0.60 <sup>6</sup></td><td>0</td><td>0.70</td><td>0.22 <sup>6</sup></td><td>0.42 <sup>6</sup></td><td>0.52 <sup>6</sup></td><td>0.60 <sup>6</sup></td><td>0</td><td>0.19 <sup>6</sup></td><td>0.38 <sup>6</sup></td><td>0.54 <sup>6</sup></td><td>0.90 <sup>9</sup></td><td>0</td><td>0.70</td><td>0.10 <sup>8</sup></td><td>0.12 <sup>10</sup></td><td>0.44 <sup>8</sup></td><td>0.44 <sup>10</sup></td><td>0.27</td></tr><tr><td>0.80</td><td>0.22 <sup>6</sup></td><td>0.42 <sup>6</sup></td><td>0.52 <sup>6</sup></td><td>0.60 <sup>6</sup></td><td>0</td><td>0.80</td><td>0.22 <sup>6</sup></td><td>0.42 <sup>6</sup></td><td>0.52 <sup>6</sup></td><td>0.60 <sup>6</sup></td><td>0</td><td>0.19 <sup>6</sup></td><td>0.38 <sup>6</sup></td><td>0.54 <sup>6</sup></td><td>0.87 <sup>9</sup></td><td>0</td><td>0.80</td><td>0.10 <sup>8</sup></td><td>0.12 <sup>10</sup></td><td>0.44 <sup>8</sup></td><td>0.44 <sup>10</sup></td><td>0.22</td></tr><tr><td>0.90</td><td>0.22 <sup>6</sup></td><td>0.42 <sup>6</sup></td><td>0.52 <sup>6</sup></td><td>0.60 <sup>6</sup></td><td>0</td><td>0.90</td><td>0.22 <sup>6</sup></td><td>0.42 <sup>6</sup></td><td>0.52 <sup>6</sup></td><td>0.60 <sup>6</sup></td><td>0</td><td>0.19 <sup>6</sup></td><td>0.38 <sup>6</sup></td><td>0.54 <sup>6</sup></td><td>0.87 <sup>9</sup></td><td>0</td><td>0.90</td><td>0.10 <sup>8</sup></td><td>0.12 <sup>10</sup></td><td>0.44 <sup>8</sup></td><td>0.44 <sup>10</sup></td><td>0.18</td></tr><tr><td>2.10</td><td>0.19 <sup>6</sup></td><td>0.38 <sup>6</sup></td><td>0.54 <sup>6</sup></td><td>0.67 <sup>9</sup></td><td>0</td><td>2.10</td><td>0.09 <sup>8</sup></td><td>0.12 <sup>6</sup></td><td>0.40 <sup>8</sup></td><td>0.39 <sup>6</sup></td><td>0.02</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 2. *Cont.*</caption><thead><tr><th colspan="11">Jamma District, Ethiopia</th><th colspan="11">Ofla District, Ethiopia</th></tr><tr><th colspan="11">Vertisol</th><th colspan="11">Vertisol</th></tr><tr><th>Depth (m)</th><th>LL<sup>1</sup></th><th>DUL<sup>2</sup></th><th>SAT<sup>3</sup></th><th>OC<sup>4</sup></th><th>SRGF<sup>5</sup></th><th>Depth (m)</th><th>LL<sup>1</sup></th><th>DUL<sup>2</sup></th><th>SAT<sup>3</sup></th><th>OC<sup>4</sup></th><th>SRGF<sup>5</sup></th></tr></thead><tbody><tr><td>0.05</td><td>0.28<sup>11</sup></td><td>0.51<sup>11</sup></td><td>0.54<sup>11</sup></td><td>1.02</td><td>1</td><td>0.05</td><td>0.26<sup>8</sup></td><td>0.44<sup>8</sup></td><td>0.51<sup>8</sup></td><td>2.55</td><td>1</td></tr><tr><td>0.10</td><td>0.28<sup>11</sup></td><td>0.51<sup>11</sup></td><td>0.54<sup>11</sup></td><td>1.02</td><td>1</td><td>0.10</td><td>0.26<sup>8</sup></td><td>0.44<sup>8</sup></td><td>0.51<sup>8</sup></td><td>2.55</td><td>1</td></tr><tr><td>0.20</td><td>0.28<sup>11</sup></td><td>0.51<sup>11</sup></td><td>0.54<sup>11</sup></td><td>1.02</td><td>1</td><td>0.20</td><td>0.26<sup>8</sup></td><td>0.44<sup>8</sup></td><td>0.51<sup>8</sup></td><td>2.55</td><td>1</td></tr><tr><td>0.30</td><td>0.28<sup>11</sup></td><td>0.51<sup>11</sup></td><td>0.54<sup>11</sup></td><td>1.02</td><td>0.61</td><td>0.30</td><td>0.29<sup>8</sup></td><td>0.40<sup>8</sup></td><td>0.49<sup>8</sup></td><td>0.06</td><td>0.61</td></tr><tr><td>0.40</td><td>0.26<sup>11</sup></td><td>0.48<sup>11</sup></td><td>0.51<sup>11</sup></td><td>0.90</td><td>0</td><td>0.40</td><td>0.29<sup>8</sup></td><td>0.40<sup>8</sup></td><td>0.49<sup>8</sup></td><td>0.06</td><td>0.50</td></tr><tr><td>0.50</td><td>0.26<sup>11</sup></td><td>0.48<sup>11</sup></td><td>0.51<sup>11</sup></td><td>0.90</td><td>0</td><td>0.50</td><td>0.29<sup>8</sup></td><td>0.40<sup>8</sup></td><td>0.49<sup>8</sup></td><td>0.06</td><td>0.41</td></tr><tr><td>0.60</td><td>0.26<sup>11</sup></td><td>0.48<sup>11</sup></td><td>0.51<sup>11</sup></td><td>0.90</td><td>0</td><td>0.60</td><td>0.29<sup>8</sup></td><td>0.40<sup>8</sup></td><td>0.49<sup>8</sup></td><td>0.06</td><td>0.33</td></tr><tr><td>0.70</td><td>0.26<sup>11</sup></td><td>0.48<sup>11</sup></td><td>0.51<sup>11</sup></td><td>0.90</td><td>0</td><td>0.70</td><td>0.29<sup>8</sup></td><td>0.40<sup>8</sup></td><td>0.49<sup>8</sup></td><td>0.06</td><td>0.28</td></tr><tr><td>0.80</td><td>0.21<sup>11</sup></td><td>0.45<sup>11</sup></td><td>0.48<sup>11</sup></td><td>0.87</td><td>0</td><td>0.80</td><td>0.29<sup>8</sup></td><td>0.40<sup>8</sup></td><td>0.49<sup>8</sup></td><td>0.06</td><td>0.22</td></tr><tr><td>0.90</td><td>0.21<sup>11</sup></td><td>0.45<sup>11</sup></td><td>0.48<sup>11</sup></td><td>0.87</td><td>0</td><td>0.90</td><td>0.29<sup>8</sup></td><td>0.40<sup>8</sup></td><td>0.49<sup>8</sup></td><td>0.06</td><td>0.18</td></tr></tbody></table><table><thead><tr><th colspan="11">Ude 1, Ethiopia</th><th colspan="11">Kajima, Ethiopia</th></tr><tr><th colspan="11">Vertisol</th><th colspan="11">Andisol</th></tr><tr><th>Depth (m)</th><th>LL<sup>1</sup></th><th>DUL<sup>2</sup></th><th>SAT<sup>3</sup></th><th>OC<sup>4</sup></th><th>SRGF<sup>5</sup></th><th>Depth (m)</th><th>LL<sup>1</sup></th><th>DUL<sup>2</sup></th><th>SAT<sup>3</sup></th><th>OC<sup>4</sup></th><th>SRGF<sup>5</sup></th></tr></thead><tbody><tr><td>0.05</td><td>0.37<sup>8,9</sup></td><td>0.50<sup>8,9</sup></td><td>0.52<sup>8,9</sup></td><td>1.25</td><td>1</td><td>0.05</td><td>0.25<sup>7,8</sup></td><td>0.41<sup>7</sup></td><td>0.49<sup>7</sup></td><td>1.91</td><td>1</td></tr><tr><td>0.10</td><td>0.37<sup>8,9</sup></td><td>0.50<sup>8,9</sup></td><td>0.52<sup>8,9</sup></td><td>1.25</td><td>1</td><td>0.10</td><td>0.25<sup>7,8</sup></td><td>0.41<sup>7</sup></td><td>0.49<sup>7</sup></td><td>1.91</td><td>1</td></tr><tr><td>0.20</td><td>0.37<sup>8,9</sup></td><td>0.50<sup>8,9</sup></td><td>0.52<sup>8,9</sup></td><td>1.25</td><td>1</td><td>0.20</td><td>0.25<sup>7,8</sup></td><td>0.41<sup>7</sup></td><td>0.49<sup>7</sup></td><td>1.91</td><td>1</td></tr><tr><td>0.30</td><td>0.37<sup>8,9</sup></td><td>0.50<sup>8,9</sup></td><td>0.52<sup>8,9</sup></td><td>1.25<sup>6</sup></td><td>0.61</td><td>0.30</td><td>0.25<sup>7,8</sup></td><td>0.41<sup>7</sup></td><td>0.49<sup>7</sup></td><td>1.91</td><td>0.61</td></tr><tr><td>0.40</td><td>0.39<sup>8,9</sup></td><td>0.49<sup>8,9</sup></td><td>0.54<sup>8,9</sup></td><td>1.25<sup>6</sup></td><td>0</td><td>0.40</td><td>0.25<sup>6</sup></td><td>0.41<sup>6</sup></td><td>0.49<sup>6</sup></td><td>1.91<sup>6</sup></td><td>0</td></tr><tr><td>0.50</td><td>0.39<sup>8,9</sup></td><td>0.49<sup>8,9</sup></td><td>0.54<sup>8,9</sup></td><td>1.25<sup>6</sup></td><td>0</td><td>0.50</td><td>0.25<sup>6</sup></td><td>0.41<sup>6</sup></td><td>0.49<sup>6</sup></td><td>1.91<sup>6</sup></td><td>0</td></tr><tr><td>0.60</td><td>0.39<sup>8,9</sup></td><td>0.49<sup>8,9</sup></td><td>0.54<sup>8,9</sup></td><td>1.25<sup>6</sup></td><td>0</td><td>0.60</td><td>0.25<sup>6</sup></td><td>0.41<sup>6</sup></td><td>0.49<sup>6</sup></td><td>1.91<sup>6</sup></td><td>0</td></tr><tr><td>0.70</td><td>0.39<sup>8,9</sup></td><td>0.49<sup>8,9</sup></td><td>0.54<sup>8,9</sup></td><td>1.25<sup>6</sup></td><td>0</td><td>0.70</td><td>0.25<sup>6</sup></td><td>0.41<sup>6</sup></td><td>0.49<sup>6</sup></td><td>1.91<sup>6</sup></td><td>0</td></tr><tr><td>0.80</td><td>0.39<sup>8,9</sup></td><td>0.50<sup>8,9</sup></td><td>0.52<sup>8,9</sup></td><td>1.25<sup>6</sup></td><td>0</td><td>0.80</td><td>0.25<sup>6</sup></td><td>0.41<sup>6</sup></td><td>0.49<sup>6</sup></td><td>1.91<sup>6</sup></td><td>0</td></tr><tr><td>0.90</td><td>0.39<sup>8,9</sup></td><td>0.50<sup>8,9</sup></td><td>0.52<sup>8,9</sup></td><td>1.25<sup>6</sup></td><td>0</td><td>0.90</td><td>0.25<sup>6</sup></td><td>0.41<sup>6</sup></td><td>0.49<sup>6</sup></td><td>1.91<sup>6</sup></td><td>0</td></tr></tbody></table><table><thead><tr><th colspan="11">Gare Arera 1, Ethiopia</th><th colspan="11">Gare Arera 2, Ethiopia</th></tr><tr><th colspan="11">Nitisol</th><th colspan="11">Vertisol</th></tr><tr><th>Depth (m)</th><th>LL<sup>1</sup></th><th>DUL<sup>2</sup></th><th>SAT<sup>3</sup></th><th>OC<sup>4</sup></th><th>SRGF<sup>5</sup></th><th>Depth (m)</th><th>LL<sup>1</sup></th><th>DUL<sup>2</sup></th><th>SAT<sup>3</sup></th><th>OC<sup>4</sup></th><th>SRGF<sup>5</sup></th></tr></thead><tbody><tr><td>0.05</td><td>0.37<sup>8,12</sup></td><td>0.45<sup>8,12</sup></td><td>0.56<sup>8,12</sup></td><td>3.77<sup>12</sup></td><td>1</td><td>0.05</td><td>0.38<sup>8,13</sup></td><td>0.53<sup>8,13</sup></td><td>0.55<sup>8,13</sup></td><td>2.59<sup>13</sup></td><td>1</td></tr><tr><td>0.10</td><td>0.37<sup>8,12</sup></td><td>0.43<sup>8,12</sup></td><td>0.56<sup>8,12</sup></td><td>3.77<sup>12</sup></td><td>1</td><td>0.10</td><td>0.38<sup>8,13</sup></td><td>0.53<sup>8,13</sup></td><td>0.55<sup>8,13</sup></td><td>2.59<sup>13</sup></td><td>1</td></tr><tr><td>0.20</td><td>0.37<sup>8,12</sup></td><td>0.45<sup>8,12</sup></td><td>0.56<sup>8,12</sup></td><td>3.77<sup>12</sup></td><td>1</td><td>0.20</td><td>0.38<sup>8,13</sup></td><td>0.53<sup>8,13</sup></td><td>0.55<sup>8,13</sup></td><td>2.59<sup>13</sup></td><td>1</td></tr><tr><td>0.30</td><td>0.37<sup>8,12</sup></td><td>0.45<sup>8,12</sup></td><td>0.54<sup>8,12</sup></td><td>2.88<sup>12</sup></td><td>0.61</td><td>0.30</td><td>0.37<sup>8,13</sup></td><td>0.52<sup>8,13</sup></td><td>0.53<sup>8,13</sup></td><td>2.05<sup>13</sup></td><td>0.61</td></tr><tr><td>0.40</td><td>0.37<sup>8,12</sup></td><td>0.44<sup>8,12</sup></td><td>0.54<sup>8,12</sup></td><td>2.88<sup>12</sup></td><td>0</td><td>0.40</td><td>0.37<sup>8,13</sup></td><td>0.52<sup>8,13</sup></td><td>0.53<sup>8,13</sup></td><td>2.05<sup>13</sup></td><td>0</td></tr><tr><td>0.50</td><td>0.34<sup>8,9</sup></td><td>0.46<sup>8,12</sup></td><td>0.52<sup>8,12</sup></td><td>2.20<sup>12</sup></td><td>0</td><td>0.50</td><td>0.37<sup>8,13</sup></td><td>0.52<sup>8,13</sup></td><td>0.54<sup>8,13</sup></td><td>2.09<sup>13</sup></td><td>0</td></tr><tr><td>0.60</td><td>0.34<sup>8,9</sup></td><td>0.46<sup>8,12</sup></td><td>0.52<sup>8,12</sup></td><td>2.20<sup>12</sup></td><td>0</td><td>0.60</td><td>0.37<sup>8,13</sup></td><td>0.52<sup>8,13</sup></td><td>0.54<sup>8,13</sup></td><td>2.09<sup>13</sup></td><td>0</td></tr><tr><td>0.70</td><td>0.34<sup>6</sup></td><td>0.45<sup>8,12</sup></td><td>0.52<sup>6</sup></td><td>2.20<sup>6</sup></td><td>0</td><td>0.70</td><td>0.37<sup>6</sup></td><td>0.52<sup>6</sup></td><td>0.54<sup>6</sup></td><td>2.09<sup>13</sup></td><td>0</td></tr><tr><td>0.80</td><td>0.34<sup>6</sup></td><td>0.46<sup>8,12</sup></td><td>0.52<sup>6</sup></td><td>2.20<sup>6</sup></td><td>0</td><td>0.80</td><td>0.37<sup>6</sup></td><td>0.52<sup>6</sup></td><td>0.54<sup>6</sup></td><td>2.09<sup>13</sup></td><td>0</td></tr><tr><td>0.90</td><td>0.34<sup>6</sup></td><td>0.46<sup>8,12</sup></td><td>0.52<sup>6</sup></td><td>2.20<sup>6</sup></td><td>0</td><td>0.90</td><td>0.37<sup>6</sup></td><td>0.52<sup>6</sup></td><td>0.54<sup>6</sup></td><td>2.09<sup>13</sup></td><td>0</td></tr></tbody></table>

<footnote n="1">¹ Soil lower limit, also known as wilting point (%volume/100). ² Soil drained upper limited, also known as field capacity (%volume/100). ³ Soil saturation point (%volume/100). ⁴ Soil organic carbon (%). ⁵ Soil root growth factor. ⁶ This layer is a copy of the layer above it. ⁷ Based on Habtegebrial et al. [37]. ⁸ This value was calculated using the DSSAT SBuild program [20]. ⁹ Based on the Vertisol description in Getu [16]. ¹⁰ From the Dia loam, 0 to 1 percent slopes group from the Fallon-Fernley Area, Nevada, parts of Churchill, Lyon, Storey, and Washoe Counties' soil survey [36]. ¹¹ The saturation soil water content (SAT) of a soil layer was based on the soil porosity, which is estimated as 100-Bulk Density/2.65*100. To consider a general but small pore space that never saturates, 3 Vol% was subtracted from the soil porosity in each layer [38]. ¹² Based on a Gare Arera Nitisol described in Tulema et al. [39]. ¹³ Based on a Gare Arera Vertisol described in Tulema et al. [39].</footnote>
<table><caption>Table 3. Crop management data.</caption><thead><tr><th>Location</th><th>Planting Date</th><th>Plant Density</th><th>N Fertilizer Application Dates</th><th>N Fertilizer Application (kg N/ha)</th><th>Irrigation</th><th>Cultivar Name</th></tr></thead><tbody><tr><td>Mekelle 1, Ethiopia</td><td>August 2</td><td>900<sup>1</sup></td><td>50% at planting and 50% 30 DAP<sup>2</sup></td><td>60<sup>3</sup></td><td>0,20,47,71, and 95 mm</td><td>DZ-974 and Keyh (local)</td></tr><tr><td>Mekelle 2, Ethiopia</td><td>July 31</td><td>900<sup>1</sup></td><td>50% at planting and 50% 30 DAP<sup>2</sup></td><td>60<sup>3</sup></td><td>0, 40, and 90 mm</td><td>Keyh (local)</td></tr><tr><td>Ilala, Ethiopia</td><td>August 5</td><td>900<sup>1</sup></td><td>50% at planting and 50% 30 DAP<sup>2</sup></td><td>60<sup>3</sup></td><td>No</td><td>DZ-974 and Keyh (local)</td></tr><tr><td>Fallon, NV, USA</td><td>June 4</td><td>900<sup>1</sup></td><td>No fertilizer applied</td><td>0</td><td>Yes</td><td>Dessie, U.K. Brown, U.K. White, PI 193508, PI 193514, PI 195932, PI 273889, PI 329680, PI 347632, PI 494366, PI 494432, PI 494433, PI 494465, PI 494479, PI 557457</td></tr><tr><td>Jamma District, Ethiopia</td><td>July 24<sup>4</sup></td><td>900<sup>1</sup></td><td>At planting</td><td>0, 23<sup>5</sup>, 46<sup>5</sup>, 69<sup>5</sup></td><td>No</td><td>Wajera (local)</td></tr><tr><td>Ofla District Ethiopia</td><td>July 24<sup>4</sup></td><td>900<sup>1</sup></td><td>At planting</td><td>0, 23<sup>5</sup>, 46<sup>5</sup>, 69<sup>5</sup></td><td>No</td><td>DZ-Cr-387 (Kuncho)</td></tr><tr><td>Ude, Ethiopia</td><td>July 22</td><td>900<sup>1</sup></td><td>On day of planting<sup>6</sup></td><td>64<sup>3</sup></td><td>No</td><td>DZ-01-354</td></tr><tr><td>Kajima, Ethiopia</td><td>July 22</td><td>900<sup>1</sup></td><td>On day of planting<sup>6</sup></td><td>64<sup>3</sup></td><td>No</td><td>DZ-01-354</td></tr><tr><td>Gare Arera 1, Ethiopia</td><td>July 24</td><td>900<sup>1</sup></td><td>Before sowing</td><td>40<sup>7</sup></td><td>No</td><td>DZ-01-354</td></tr><tr><td>Gare Arera 2, Ethiopia</td><td>July 25</td><td>900<sup>1</sup></td><td>Before sowing</td><td>60<sup>7</sup></td><td>No</td><td>DZ-01-354</td></tr></tbody></table>

<table><caption>Table 4. Initial estimated soil water and mineral nitrogen conditions for each location at the start of the simulation. Simulations started one month before planting.</caption><thead><tr><th>Location</th><th>Initial Soil Water (%FC<sup>1</sup>)</th><th>Initial Soil N (kg N/ha)</th></tr></thead><tbody><tr><td>Mekelle 1, Ethiopia</td><td>100<sup>2</sup></td><td>10</td></tr><tr><td>Mekelle 2, Ethiopia</td><td>100<sup>2</sup></td><td>10</td></tr><tr><td>Ilala, Ethiopia</td><td>100<sup>2</sup></td><td>10</td></tr><tr><td>Fallon, NV, USA</td><td>100</td><td>40</td></tr><tr><td>Jamma District, Ethiopia</td><td>10</td><td>30</td></tr><tr><td>Ofla District Ethiopia</td><td>10</td><td>150<sup>3</sup></td></tr><tr><td>Ude, Ethiopia</td><td>10</td><td>20</td></tr><tr><td>Kajima, Ethiopia</td><td>10</td><td>20</td></tr><tr><td>Gare Arera 1, Ethiopia</td><td>40</td><td>10</td></tr><tr><td>Gare Arera 2, Ethiopia</td><td>40</td><td>10</td></tr></tbody></table>

## 2.4. Phenology

General phenology information for tef was taken from the existing literature. The base temperature used to calculate developmental thermal time accumulation was increased from the 0 °C used in the wheat model [20] to 7.8 °C for tef [6]. The vernalization requirement was reduced to 0 in the ecotype file, as tef does not have a vernalization requirement [6].

There was no clear optimum temperature for tef phenology available in the literature. Reported optimal growing temperature ranges included 10–27 °C [2], 15–21 °C [41], and >26 °C [7]. It is unclear if these ranges referred to minimum, average, or maximum daily temperatures. Tsegay [42] wrote that the maximum temperature for tef production was 35 °C. As the literature did not specify for tef if the above temperatures were related to phenology or growth, the DSSAT-Tef model used the

<footnote n="1">¹ Rounded from the 923 plants/m² used by Araya et al. [1]. ² Days after planting. ³ The specific type of fertilizer was not specified. It was assumed to be conventional urea. ⁴ Based on Settie et al. [40]. ⁵ As urea N super granule (USG). ⁶ Assumed, as it was not specified in Tulema et al. [27]. ⁷ As a combination of conventional urea and diammonium phosphate.</footnote>

<footnote n="1">¹ Field capacity. ² After Araya et al. [13]. ³ 91% of soil N was concentrated in the top 30 cm.</footnote>
optimum and maximum temperatures from APSIM-Sorghum [26] (Table 5), which is a C4 crop, like tef [23]. APSIM-Sorghum, unlike DSSAT-Sorghum, also used the same three temperature parameters for calculating phenology as DSSAT-NWheat, on which the DSSAT-Tef model was based (Table 5).

<table><caption>**Table 5.** Temperature parameters for phenology for different crop models.</caption><thead><tr><th>Model</th><th>Reference</th><th>TBASE <sup>1</sup></th><th>TOPT <sup>2</sup></th><th>TTOPT <sup>3</sup></th></tr></thead><tbody><tr><td>DSSAT-NWheat</td><td>[20]</td><td>0</td><td>26</td><td>34</td></tr><tr><td>APSIM-Sorghum</td><td>[26]</td><td>11</td><td>30</td><td>42</td></tr><tr><td>APSIM-Millet</td><td>[43]</td><td>10</td><td>33</td><td>47</td></tr><tr><td>DSSAT-Sorghum</td><td>[20]</td><td>8</td><td>34</td><td>NA</td></tr><tr><td>DSSAT-Millet</td><td>[20]</td><td>10</td><td>34</td><td>NA</td></tr><tr><td>DSSAT-Tef</td><td></td><td>7.8</td><td>30</td><td>42</td></tr></tbody></table>

Phenology differences between the individual cultivars were controlled using the photoperiod sensitivity parameter (PPSEN), the thermal time from seedling emergence to the end of tillering (juvenile phase) (P1) parameter, the phyllochron interval parameter (PHINT), and the thermal time from the beginning of grain fill to maturity (P5) parameter in the cultivar file (Table 6). These parameters were individually changed for each cultivar in order to fit observed phenology data.

<table><caption>**Table 6.** Cultivar parameters in the DSSAT model.</caption><thead><tr><th>Cultivar Name</th><th>PPSEN <sup>1</sup></th><th>P1 <sup>2</sup></th><th>P5 <sup>3</sup></th><th>PHINT <sup>4</sup></th><th>GRNO <sup>5</sup></th><th>MXFIL <sup>6</sup></th><th>STMMX <sup>7</sup></th></tr></thead><tbody><tr><td>Dessie</td><td>16.5</td><td>100</td><td>500</td><td>80</td><td>32</td><td>1.2</td><td>1.0</td></tr><tr><td>Keyh (local)</td><td>18.0</td><td>100</td><td>600</td><td>110</td><td>16</td><td>1.0</td><td>3.0</td></tr><tr><td>U.K. Brown</td><td>16.5</td><td>100</td><td>450</td><td>100</td><td>16</td><td>1.0</td><td>1.0</td></tr><tr><td>U.K. White</td><td>18.0</td><td>100</td><td>450</td><td>85</td><td>16</td><td>1.0</td><td>1.0</td></tr><tr><td>Wajera (local)</td><td>18.0</td><td>100</td><td>200</td><td>80</td><td>36</td><td>1.0</td><td>3.0</td></tr><tr><td>Combined DZ-974 and Keyh (local)</td><td>18.0</td><td>100</td><td>500</td><td>105</td><td>16</td><td>1.0</td><td>1.0</td></tr><tr><td>DZ-01-354</td><td>18.0</td><td>100</td><td>550</td><td>100</td><td>16</td><td>1.0</td><td>3.0</td></tr><tr><td>DZ-Cr-387</td><td>18.0</td><td>100</td><td>100</td><td>120</td><td>32</td><td>1.3</td><td>3.0</td></tr><tr><td>PI 193508</td><td>19.0</td><td>100</td><td>200</td><td>100</td><td>30</td><td>1.1</td><td>1.0</td></tr><tr><td>PI 193514</td><td>19.0</td><td>100</td><td>200</td><td>105</td><td>16</td><td>1.0</td><td>1.0</td></tr><tr><td>PI 195932</td><td>19.0</td><td>100</td><td>200</td><td>105</td><td>16</td><td>1.0</td><td>1.0</td></tr><tr><td>PI 273889</td><td>19.0</td><td>100</td><td>190</td><td>105</td><td>27</td><td>1.0</td><td>1.0</td></tr><tr><td>PI 329680</td><td>19.0</td><td>100</td><td>200</td><td>105</td><td>16</td><td>1.0</td><td>1.0</td></tr><tr><td>PI 347632</td><td>19.0</td><td>100</td><td>190</td><td>105</td><td>32</td><td>1.8</td><td>1.0</td></tr><tr><td>PI 494366</td><td>14.0</td><td>100</td><td>700</td><td>70</td><td>16</td><td>1.0</td><td>1.0</td></tr><tr><td>PI 494432</td><td>15.5</td><td>100</td><td>500</td><td>100</td><td>26</td><td>1.0</td><td>1.0</td></tr><tr><td>PI 494433</td><td>19.0</td><td>100</td><td>140</td><td>120</td><td>16</td><td>1.0</td><td>1.0</td></tr><tr><td>PI 494465</td><td>19.0</td><td>100</td><td>140</td><td>120</td><td>16</td><td>1.0</td><td>1.0</td></tr><tr><td>PI 494479</td><td>19.0</td><td>100</td><td>200</td><td>105</td><td>16</td><td>1.0</td><td>1.0</td></tr><tr><td>PI 557457</td><td>17.0</td><td>100</td><td>600</td><td>75</td><td>16</td><td>1.0</td><td>1.0</td></tr></tbody></table>

## 2.5. Photoperiod Response

DSSAT-NWheat uses the parameter nwheats_ppfac to account for photoperiod response. The NWheat equation was changed to the following for tef:

$$ tef_{ppfac} = 1 - 0.002 * PPSEN * (TWILEN - 11)^2 \quad (1) $$

where tef_ppfac is a unitless photoperiod response factor; PPSEN is a unitless photoperiod sensitivity factor that varies between cultivars; and TWILEN is the daylength. For NWheat, the ppfac parameter increases with increasing daylength. Tef is a short-day plant [6], so the photoperiod response needed

<footnote n="1">¹ Base daily mean temperature below which no crop phenology occurs (°C). ² Optimum daily mean temperature at which the maximum rate of crop phenology occurs (°C). ³ Maximum daily mean temperature at which crop phenology stops (°C).</footnote>

<footnote n="1">¹ Sensitivity to photoperiod factor. ² Thermal time from seedling emergence to the end of tillering (juvenile phase).</footnote>

<footnote n="3">³ Thermal time from beginning of grain fill to maturity. ⁴ Phyllochron interval. ⁵ Coefficient of kernel number per stem weight at the beginning of grain filling (100 kernels/g stem). ⁶ Potential kernel growth rate (mg/kernel/day).</footnote>

<footnote n="7">⁷ Potential final dry weight of a single tiller (excluding grain) (g/stem).</footnote>
to be changed so that tef_ppfac decreased with increasing day length, and as a result, the rate of
development decreased in response to longer days. Van Delden [6] defined the critical daylength
for short day plants as the day length above which day length begins to impact the time to heading.
The critical daylength of 20 h used in NWheat was reduced to 11 h, which was the average critical
daylength across four tef cultivars [6]. The photoperiod sensitivity of each cultivar (PPSEN) was
estimated by using phenology data from field trials (Tables 1 and 6). As tef is extremely photoperiod
sensitive, the value for PPSEN used in the tef model was significantly higher than the typical range of
1–5 used by the NWheat model (Table 6). Figure 1 shows the response of tef_ppfac across multiple
daylengths and PPSEN values.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** The photoperiod response factor in response to daylength using a photoperiod sensitivity parameter (PPSEN) ranging from 1 (blue) to 19 (red).</figcaption></figure>

2.6. Growth

Though there were no radiation use efficiency (RUE) data available for tef, it was assumed
that, as a C4 crop, tef would have a higher RUE than wheat. Models for the C4 crops maize, millet,
and sorghum all have higher aboveground RUE values than wheat models [20,44–46]. In order to
calculate the DSSAT-Tef aboveground RUE, the percent difference between the DSSAT-CERES-Wheat
and DSSAT-Sorghum aboveground RUE was used to estimate the ratio RUE ratio between wheat
and a C4 crop. The percent change in RUE between the DSSAT-CERES-Wheat and DSSAT-Sorghum
model was 18.5%. This percent increase was applied to the RUE of DSSAT-NWheat to estimate the
RUE for tef. The RUE for DSSAT-NWheat was 3.2 g plant dry matter/MJ PAR, so the RUE value for
DSSAT-Tef was set to 4.5 g/MJ. The NWheat model used a weighted mean temperature with emphasis
on day time temperature [47]. To account for a higher base temperature for photosynthesis in tef, the
threshold minimum temperature for photosynthesis was set to 4 degrees below the base temperature
for phenology at 3.8 °C, the same difference, but for a different base temperature used in NWheat.
These new parameter values were used in all tef simulations without any additional calibration.

The DSSAT-NWheat model initialized the grain weight at the beginning of grain filling using
the initial grain weight parameter (INGWT). Growing conditions throughout the grain fill period
and cultivar parameters, such as grain number (GRNO), potential kernel growth rate (MXFIL),
photosynthesis, and available water soluble carbohydrates, all affect the change in kernel grain weight
from this initial value. Tef grains are much smaller than wheat grains [31], so the initial grain weight
needed to be decreased for the tef model. The first step to estimating the initial grain weight for
DSSAT-Tef was calculating the ratio of the initial grain weight parameter used by the DSSAT-NWheat
model [20] and the final mature grain weight (mg/grain) at 0% grain moisture of an Ethiopian wheat field experiment [31]. This ratio was then applied to the final grain weight of tef at 0% moisture grown at the same location and under the same management conditions [31], in order to calculate the initial grain weight parameter for tef. Based on the estimated initial grain weights, INGWT was changed from 3.50 mg/grain for NWheat to 2.47 mg/100 grains for DSSAT-Tef. The tef model used mg/100 grains in its input file, but once the INGWT was read into the code, it was divided by 100 to get mg/grain.

The units for the grain number cultivar parameter in NWheat were the number of kernels/g stem weight at the beginning of grain filling and affected the rate of grain filling [20]. As tef has a much higher grain number per m² than wheat, the GRNO units in the cultivar file were changed to 100 kernels/g stem weight at the beginning of grain filling for the DSSAT-Tef model.

Lodging is a concern for tef production, especially for high input growing scenarios [37]. Due to a lack of lodging data in the available experiments, except for an indication of no or little occurrence of lodging in these datasets, no attempt was made to include lodging in the DSSAT-Tef model.

## 2.7. CO₂ Response

Though there were no published field experiments on the effects of elevated atmospheric CO₂ on tef, it is likely that, as a C4 crop, tef would have a less pronounced response to elevated CO₂ than wheat. It has been found that transpiration efficiency increases with increasing CO₂ levels for C4 crops, while transpiration efficiency and RUE increase for C3 crops [48]. DSSAT considered atmospheric CO₂ concentration as an input [20]. DSSAT-NWheat had two major pathways for dealing with the effects of elevated CO₂, transp_eff Coeff, and rue_factor. The coefficient transp_eeff Coeff is the transpiration efficiency coefficient and is calculated using:

$$transp\_eff\_coeff = 0.009 * (1 + (0.28/350) * (WEATHER\% CO_2 - 350)) \quad (2)$$

where WEATHER% CO₂ is the actual atmospheric CO₂ concentration in ppm. Hammer and Muchow [49] used a transpiration efficiency coefficient of 0.009 kPa for their sorghum simulation. As Hammer and Muchow [49] based their calculations only on aboveground biomass, their transpiration efficiency coefficient was compatible with the DSSAT-NWheat approach for modeling transpiration [20]. Pembleton et al. [50] used the APSIM-Wheat model as a basis for estimating the functions for modifying transpiration efficiency for various forage crops. Their function for forage sorghum was y = 0.0008x + 1. Based on Pembleton et al.'s [50] equation, a 350 ppm increase in CO₂ would result in a 28% increase in transpiration efficiency in sorghum. Teixera et al. [51] arbitrarily assumed a 37% increase in transpiration efficiency for maize and a 69% increase for wheat when atmospheric CO₂ was increased from 350 to 1032 ppm. Assuming a linear response, this would correspond to a 19% and 35% increase in transpiration efficiency for maize and wheat, respectively, if CO₂ were doubled from 350 to 700 ppm. As the calculations of Pembleton et al. [50] were based on field data, rather than an arbitrary assumption, their CO₂ response was used for DSSAT-Tef.

As RUE for C4 crops are not influenced by increasing atmospheric CO₂ [48], the rue_factor was set to 1 for DSSAT-Tef. This was accomplished by setting the parameter RUEFAC in the ecotype file to 0, which set rue_factor to a constant of 1.

## 2.8. Cultivar Parameters

After the changes described above were made to the code, the individual cultivars were parameterized one at a time to account for differences in phenology and productivity. Table 6 gives the parameter values used for each cultivar. Cultivar parameters were calibrated in a multistep process, starting with the time to anthesis and followed by time to maturity, total aboveground biomass, and grain yield. PPSEN, P1, and PHINT affected time to anthesis. P5 affected the time to maturity. STMMX affected biomass yields. GRNO and MXFIL affected grain yields. Only one parameter was changed for each run, in order to avoid unintended interaction effects. Tef has high genetic diversity
across cultivars. Dawit and Hirut [52] reported that days to 50% flowering ranged between 43 and 83 days across 506 tef accessions and that days to 50% maturity ranged from 93 to 130 days.

DSSAT did not output the days until heading, which means that observed anthesis dates were needed for calibration. Many of the publications used for the DSSAT-Tef model reported days until heading, rather than days until anthesis, however. For the Ofla District location [31], 10 days were added to the days to 50% heading to approximate the days until anthesis. This estimation was based on the difference between days to heading and days to flowering reported by Getu [16]. No anthesis information was reported for the Ude, Kajima, and both Gare Arera locations [27,35]. The observed days until anthesis were approximated by adding 10 days to the days until heading that were reported for the same cultivar grown under a 12 h daylength by Tadelle Tadesse [53]. The Fallon, NV, location [29] reported days to full seed head emergence, which was assumed to be equivalent to anthesis, so no days were added to the reported dates.

## 2.9. Statistics

All statistical analysis was done using the basic R package [54]. The simulated results were compared to the published observed data using the $r^2$, the root mean squared error (RMSE), and the relative root mean squared error (rRMSE). The $r^2$ was calculated for the 1:1 line to quantify the departure of the model from the ideal reproduction of each observed value [55]. The RMSE and rRMSE were also calculated based on the simulated versus observed 1:1 line.

## 3. Results

The performance of the model compared with observed data is summarized in Table 7. Observed days from planting to emergence ranged from 5–6 days, but the observed data were only available for three locations. The model captured this small range very well, with a coefficient of determination ($r^2$) of 1.00, an RMSE of 0.8 days, and an rRMSE of 14.4%. Unlike for emergence, there was a large range in observed days from planting to anthesis and maturity, which ranged from 44–98 days and 71–150 days, respectively. The model accurately simulated days to anthesis, with an $r^2$ of 0.98, an RMSE of 2.5 days, and days to maturity with an $r^2$ of 0.96 and an RMSE of 4.4 days. Figure 2 shows the simulated versus observed days from planting to anthesis and days from planting to maturity for each location. There were no reported phenology data for Ude or Kajima, so the observed phenology for the same cultivar (DZ-01-354) from a different study [53] was used for calibration. As both Ude and Kajima shared the same coordinates, and therefore weather data, the sites were lumped together as the Ada Region when graphed.

<table><caption>**Table 7.** Summary of the DSSAT-Tef model performance at Mekelle, Ilala, Jamma District, Ofla District, Ude, Kajima, and Gare Arera in Ethiopia, and Fallon, NV, USA.</caption><thead><tr><th>Model Attribute</th><th>Number of Paired Data Points</th><th>Observed Range</th><th>R<sup>1,2</sup></th><th>M<sup>2</sup></th><th>RMSE<sup>3</sup></th><th>rRMSE<sup>4</sup></th></tr></thead><tbody><tr><td>Emergence (DAP)<sup>5</sup></td><td>9</td><td>5–6</td><td>1.00</td><td>0.876</td><td>0.8</td><td>14.4</td></tr><tr><td>Anthesis (DAP)<sup>5</sup></td><td>34</td><td>44–98</td><td>0.98</td><td>1.005</td><td>2.5</td><td>3.6</td></tr><tr><td>Maturity (DAP)<sup>5</sup></td><td>34</td><td>71–150</td><td>0.96</td><td>1.008</td><td>4.4</td><td>3.9</td></tr><tr><td>Total Aboveground Biomass (kg/ha)</td><td>23</td><td>1400–9652</td><td>0.80</td><td>1.396</td><td>2624</td><td>49.6</td></tr><tr><td>Grain yield (kg/ha)</td><td>38</td><td>352–3013</td><td>0.74</td><td>1.159</td><td>475</td><td>41.0</td></tr><tr><td>Crop N Uptake (kg/ha)</td><td>10</td><td>20–206</td><td>0.81</td><td>0.623</td><td>45</td><td>46.2</td></tr><tr><td>Grain N (kg/ha)</td><td>12</td><td>4.3–92</td><td>0.79</td><td>0.794</td><td>15</td><td>37.3</td></tr></tbody></table>

<footnote n="1">¹ $r^2$ for the 1 to 1 line of simulated versus observed data. ² The slope *m* of a linear regression (forced through origin) is shown as a quantitative indication of over- or under-estimation of a model's performance. ³ Root mean squared error. ⁴ Relative root mean squared error. ⁵ Days after planting.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Model performance for (a) anthesis and (b) maturity, in terms of days after planting. Simulated versus observed days after planting at Mekelle, Ethiopia, in 2008 (red filled square), Mekelle, Ethiopia, in 2009 (green filled triangle), Ilala, Ethiopia, in 2008 (blue filled diamond), Fallon, NV, USA, in 2009 (black filled dot), Jamma District, Ethiopia, in 2011 (open purple square), Ofla District, Ethiopia, in 2011 (open green diamond), and Ada Region (Ude and Kajima), Ethiopia, in 1997 (open upward yellow triangle). The solid black line is the 1:1 line. No error bars are available.</figcaption></figure>

Total aboveground biomass simulations were less accurate than the phenology simulations with an r² of 0.80 and an RMSE of 2625 kg/ha (rRMSE 49.7%). Figure 3a shows simulated versus observed total aboveground biomass for each location. The model tended to overestimate total aboveground biomass. Figure 4 indicates that the model consistently captured the increasing trend in total aboveground biomass in response to increasing applied nitrogen. The model overestimated simulated total aboveground biomass more severely for water stressed treatments, such as the zero irrigation treatments at Mekelle and Iala (leftmost red square and blue diamond in Figure 3a), than for irrigated treatments. Figure 5 shows that the model captured the increase in total aboveground biomass, in response to increasing irrigation amounts, for all treatments, except one.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Model performance for (a) total aboveground biomass and (b) grain yield, all at 0% moisture. Simulated versus observed yields at Mekelle, Ethiopia, in 2008 (red filled square), Mekelle, Ethiopia, in 2009 (green filled triangle), Ilala, Ethiopia, in 2008 (blue filled diamond), Fallon, NV, USA, in 2009 (black filled dot), Jamma District, Ethiopia, in 2011 (open purple square), Ofla District, Ethiopia, in 2011 (open green diamond), Ada Region (Ude and Kajima), Ethiopia, in 1997 (open upward yellow triangle), and both soils in Gare Arera, Ethiopia, in 2002 (open downward blue triangle). No observed biomass data were available for the Fallon, NV, USA, site. The solid black line is the 1:1 line. No error bars are available.</figcaption></figure>
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Effects of N fertilizer on (a,b) total aboveground biomass and (c,d) on grain yield for (a,c) Jamma District in 2011 and (b,d) Ofla District in 2011. Bars represent the observed data, while the lines represent the simulated data. Observed data from Getu [16] and Okubay [31]. No error bars are available.</figcaption></figure>

Observed grain yield, at 0% grain moisture, varied from 352-3013 kg/ha. Grain yield simulations had an r² of 0.74 and an RMSE of 475 kg/ha (rRMSE 41.0%). Figure 3b shows simulated versus observed grain yield for each location. Grain yield was, in general, well simulated, except for some overestimations, such as the Ilala and irrigated Mekelle 2009 treatments (Figure 3b). There seemed to be no clear relationship between the amount of water or nitrogen applied and the accuracy of the grain yield simulation (Figure 6), suggesting that there was no model bias towards high or low input scenarios (Figures 4 and 5). The model consistently captured the increasing trend in grain yield in response to increasing nitrogen and irrigation applications (Figures 4 and 5).
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Effects of irrigation on (a,b) total aboveground biomass and (c,d) on grain yield for (a,c) Mekelle, Ethiopia, in 2008 and (b,d) Mekelle, Ethiopia, in 2009. Bars represent the observed data, while the lines represent the simulated data. Observed data from Araya et al. [13]. No error bars are available.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 6.** (a) Simulated-observed grain yield for tef versus applied irrigation for Mekelle in 2008 and 2009. (b) Simulated-observed grain yield for tef versus applied N for Mekelle 2008, Mekelle 2009, the Jamma District in 2011, the Ofla District in 2011, the Ada Region (Ude and Kajima) in 1997, and Gare Arera in 2002. Observed data from Araya et al. [13], Getu [16], Okubay [31], and Tulema et al. [27,35]. No error bars are available.</figcaption></figure>
Observed total aboveground nitrogen content ranged from 20–206 kg N/ha. There were no in-season measurements available. The model simulated total aboveground nitrogen with an r² of 0.80 and an RMSE of 45 kg N/ha (rRMSE 46.2%). Figure 7a shows simulated and observed total aboveground nitrogen for each location that had observed data. The model underestimated total aboveground nitrogen at higher N application treatments, but treatments with no added nitrogen were close to the observed values. The model captured the increase in total aboveground nitrogen with increasing nitrogen applications.

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Model performance for (a) total aboveground nitrogen and (b) grain nitrogen. Simulated versus observed yields at Jamma District, Ethiopia in 2011 (open purple square), Ofla District, Ethiopia, in 2011 (open green diamond), Ada Region (Ude and Kajima), Ethiopia, in 1997 (open upward yellow triangle), and both soils in Gare Arera, Ethiopia, in 2002 (open downward blue triangle). Only grain nitrogen was available for the Gare Arera sites, and only total nitrogen was available for the Ada Region sites. The solid black line is the 1:1 line. No error bars are available.</figcaption></figure>

Observed grain nitrogen ranged from 4–92 kg N/ha. There were no in-season measurements available. The model simulated grain nitrogen with an r² of 0.78 and an RMSE of 15 kg N/ha (rRMSE 37.3%). Figure 7b shows simulated and observed grain nitrogen for each location that had observed data. The model did not consistently over- or under-estimate the grain nitrogen. While treatments with no applied nitrogen fertilizer gave more accurate simulations for the Jamma District and Gare Arera locations, the zero nitrogen treatment for the Ofla District location was more off than the fertilized treatments at that location. This was possibly due to the high initial N used for the Ofla District location.

**4. Discussion**

A new DSSAT-Tef model was created based on DSSAT-NWheat and information from the literature. Phenology, photoperiod response, radiation use efficiency, transpiration efficiency, and atmospheric CO₂ response were modified to take the specifics of tef into account. The DSSAT-Tef model could reproduce observed phenology data and the observed yield responses to water and nitrogen, without any structural changes to the original model, using key crop parameters (like RUE and transpiration efficiency) from the literature without additional calibration and some calibration of a limited set of cultivar parameters. The range in parameter values was in accordance with the high variation in agronomic characteristics found in tef [2].

The comparison of the DSSAT-Tef model with a range of published field data showed an RMSE for total aboveground biomass of 2624 kg/ha and for grain yield of 475 kg/ha. The rRMSE for total aboveground biomass for the DSSAT-Tef model was 49.6%. DSSAT-Tef's biomass yield rRMSE was higher than the average 15% rRMSE of the DSSAT-NWheat model [22]. The APSIM-Sorghum and
DSSAT-Sorghum models [56] also reported lower rRMSEs for total aboveground biomass (11.5% and 12.8%, respectively) than the DSSAT-Tef model.

The rRMSE for grain yield for the DSSAT-Tef model was 41.0%. DSSAT-Tef's grain yield rRMSE was higher than the average 13% rRMSE of the DSSAT-NWheat model [22] and the 20.3% rRMSE of the APSIM-Sorghum model [56]. DSSAT-Tef's grain yield rRMSE was, however, comparable to the 39.5% rRMSE of the DSSAT-Sorghum model [56].

The DSSAT-Tef model's grain yield RMSE was comparable to the two existing tef models, but the model overestimated biomass production more than the other two models [12, 13, 17, 18]. The FAO-AEZ tef model was applied at 13 locations [12]. The FAO-AEZ grain yield RMSE was 402 kg/ha (rRMSE 29.1%). The FAO-AquaCrop tef model has been tested at multiple locations in multiple papers [13,17,18]. The AquaCrop final biomass RMSE ranged from 466 kg/ha to 720 kg/ha (rRMSE 23.6%) [13,18]. The AquaCrop grain yield RMSE ranged from 116 kg/ha to 540 kg/ha (rRMSE 19.1%) [13,17].

The biomass RMSE of DSSAT-Tef was higher than that of AquaCrop. The higher RMSE might be the result of DSSAT-Tef being tested with data that came from a greater range of conditions than the data used for AquaCrop. All AquaCrop test data came from fertilized experiments in northern Ethiopia [1,13,17,18] and did not include the variability caused by N stress and photoperiod response that was included in the DSSAT-Tef test data. The DSSAT-Tef grain yield RMSE fell within the RMSE range for AquaCrop, but the DSSAT-Tef grain yield rRMSE was higher than that of the AquaCrop model [13,17]. The DSSAT-Tef grain yield RMSE was slightly higher than that of the FAO-AEZ model, and the DSSAT-Tef's grain yield rRMSE was also higher than that of the FAO-AEZ model. A lack of variability in test data for the FAO-AEZ model could be why the grain yield RMSE of the FAO-AEZ model was comparable to that of DSSAT-Tef, despite the fact that the FAO-AEZ model assumed constant management practices across all locations and used seasonal averages and totals for weather data, rather than daily inputs [12].

Unlike the FAO-AEZ and AquaCrop tef models, DSSAT-Tef could model total aboveground N and grain N, which could be used to evaluate the fodder quality of tef straw and the protein content of tef grain. DSSAT-Tef could also respond to more environmental variation as it accounted for the effects of N, photoperiod, and atmospheric CO₂, although the later was not tested with field data due to a lack of experiments. Further research is needed to improve the model's ability to simulate N demand and partitioning. Further research is also needed to confirm that tef has a similar response to elevated CO₂ as sorghum. If this research becomes available, the DSSAT-Tef model could be used to assess the impacts of climate change on tef production and food security in Ethiopia. More detailed field data might also confirm the explanations proposed below for some of the DSSAT-Tef model responses.

The DSSAT-Tef model overestimated the days to maturity at the Ilala location [13]. This in turn contributed to the overestimation of total aboveground biomass and grain yield. The reported phenology and yield values for both the Mekelle location in 2008 and the Ilala location were an average of two different cultivars, the improved DZ-974 and the local Keyh variety. It was not reported what the ratio of improved to local variety was for each location, and it is possible that the differences in observed phenology for each location were the result of a difference in the percentage of each cultivar present. The combination of the two cultivars was treated as a single cultivar when calibrating the cultivar parameters. As the same cultivars were used for each location, it was assumed that the cultivar parameters were the same for each location. Therefore, the discrepancy between observed and simulated data at the Ilala location could not be resolved by altering cultivar parameters.

The observed days to anthesis and maturity varied across fertilizer treatments for the Jamma and Ofla District locations [16,31]. The model did not account for this spread, however, as DSSAT used temperature and daylength to determine phenology and not nitrogen availability [20]. Neither publication could confirm the reason for the range in phenology, though it was proposed that tef reduced the time until heading when grown under lower N conditions as an avoidance response. By developing more rapidly, the plant would be able to reproduce before depleting the limited soil nutrients [31].
The simulated total aboveground biomass and grain yield at the Ofla District location followed a strong linear trend for the first three fertilizer treatments, but spiked for the highest N treatment. A possible explanation for this is that the last treatment suffered from lodging (although there was no indication in the report), for which the model did not account. Tef is particularly prone to lodging at fertilizer application rates above 60 kg N/ha [37], and the highest fertilizer rate at the Ofla District location was 69 kg N/ha [31]. Future research should consider including a lodging routine for the DSSAT-Tef model to address this issue.

Other causes for discrepancies between the model and the observed data included differences in initial conditions or crop management practices and biological stresses. Not all studies thoroughly described the initial conditions or crop management practices, so in some cases, assumptions had to be made, like on rooting depth, soil moisture, and plant density. The model only accounted for abiotic stresses, so if there were yield losses due to unreported abiotic stresses, the model would not capture them.

**5. Conclusions**

The DSSAT-Tef model offered new modeling capabilities beyond the two existing tef models FAO-AEZ and AquaCrop tef, by including photoperiod sensitivity, crop management, crop N dynamics, and the impact of elevated atmospheric CO₂ concentration. While the impact of elevated atmospheric CO₂ concentrations on tef growth was based on a general understanding of crop response to CO₂, this was not evaluated with observations due to a lack of field data. Further research is also needed to improve N simulations and to add a lodging routine to the model. Additional field data are needed to do a complete evaluation of the model. The DSSAT-Tef model can be used to assess management practices and the suitability of regions for growing tef both in Ethiopia and other parts of the world.

**Author Contributions:** Conceptualization, K.P. and S.A.; data curation, K.P.; formal analysis, K.P.; funding acquisition, K.P.; investigation, K.P.; methodology, K.P. and S.A.; project administration, S.A.; resources, K.P. and S.A.; software, K.P.; supervision, S.A.; validation, K.P. and S.A.; visualization, K.P. and S.A.; writing, original draft, K.P.; writing, review and editing, K.P. and S.A.

**Funding:** K.P. received funding from the University of Florida's Graduate Assistantship Program.

**Acknowledgments:** We thank B. Zaitchik for providing the CHIRPS weather data from Ethiopia and M. Eggen and N. Huth for providing access to tef field data.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Araya, A.; Stroosnijder, L.; Girmay, G.; Keesstra, S. Crop coefficient, yield response to water stress and water productivity of teff (Eragrostis tef (Zucc.)). *Agric. Water Manag.* **2011**, *98*, 775-783. [CrossRef]

2. Seyfu, K. Tef, Eragrostis tef (Zucc.) Trotter: Promoting the Conservation and Use of Underutilized and Neglected Crops; Institute of Plant Genetics and Crop Plant Research: Gatersleben, Germany; International Plant Genetic Resources Institute: Rome, Italy, 1997.

3. Vohwinkel, F.; Piepho, H.P.; Heiligtag, B.; Richter, C. Comparative Yield and Nutrient Concentration of Six Cultivars of Tef (Eragrostis Tef (Zucc.) Trotter). *Ethiop. J. Nat. Resour.* **2002**, *4*, 21-36.

4. Costanza, S.H.; Dewet, J.M.J.; Harlan, J.R. Literature-Review and Numerical Taxonomy of Eragrostis-tef (Tef). *Econ. Bot.* **1979**, *33*, 413-424. [CrossRef]

5. Spaenij-Dekking, L.; Kooy-Winkelaar, Y.; Koning, F. The ethiopian cereal tef in celiac disease. *N. Engl. J. Med.* **2005**, *353*, 1748-1749. [CrossRef] [PubMed]

6. Van Delden, S.H. *On Seed Physiology, Biomechanics and Plant Phenology in Eragrostis Tef*; Wageningen University: Wageningen, The Netherlands, 2011.

7. Miller, D. *Teff Grass Crop Overview and Forage Production Guide*, 3rd ed.; Seeds, C., Ed.; CalWest Seeds: Woodland, CA, USA, 2011.

8. Asseng, S. *Wheat Crop Systems: A Simulation Analysis*; CSIRO: Melbourne, Australia, 2004.
9. Chenu, K.; Porter, J.; Martre, P.; Basso, B.; Chapman, S.; Ewert, F.; Bindi, M.; Asseng, S. Contribution of Crop Models to Adaptation in Wheat. *Trends Plant Sci.* **2017**, *22*, 472–490. [CrossRef] [PubMed]

10. Asseng, S.; Ewert, F.; Martre, P.; Rotter, R.; Lobell, D.; Cammarano, D.; Kimball, B.; Ottman, M.; Wall, G.; White, J.; et al. Rising temperatures reduce global wheat production. *Nat. Clim. Chang.* **2015**, *5*, 143–147. [CrossRef]

11. Rosenzweig, C.; Elliott, J.; Deryng, D.; Ruane, A.; Muller, C.; Arneth, A.; Boote, K.; Folberth, C.; Glotter, M.; Khabarov, N.; et al. Assessing agricultural risks of climate change in the 21st century in a global gridded crop model intercomparison. *Proc. Natl. Acad. Sci. USA* **2014**, *111*, 3268–3273. [CrossRef]

12. Yizengaw, T.; Verheye, W. Modeling Production Potentials of Tef (*Eragrostis-Tef*) in the Central Highlands of Ethiopia. *Soil Technol.* **1994**, *7*, 269–277. [CrossRef]

13. Araya, A.; Keesstra, S.; Stroosnijder, L. Simulating yield response to water of Teff (*Eragrostis tef*) with FAO's AquaCrop model. *Field Crop. Res.* **2010**, *116*, 196–204. [CrossRef]

14. Steduto, P.; Hsiao, T.C.; Raes, D.; Fereres, E. AquaCrop-The FAO Crop Model to Simulate Yield Response to Water: I. Concepts and Underlying Principles. *Agron. J.* **2009**, *101*, 426–437. [CrossRef]

15. Hailelassie, A.; Priess, J.; Veldkamp, E.; Teketay, D.; Lesschen, J.P. Assessment of soil nutrient depletion and its spatial variability on smallholders' mixed farming systems in Ethiopia using partial versus full nutrient balances. *Agric. Ecosyst. Environ.* **2005**, *108*, 1–16. [CrossRef]

16. Getu, A. Soil Characterization and Evaluation of Slow Release Urea Fertilizer Rates on Yield Components and Grain Yields of Wheat and Teff on Vertisols of Jamma District of South Wollo Zone, Amhara Region. Master's Thesis, Haramaya University, Haramaya, Ethiopia, 2012.

17. Tsegay, A.; Raes, D.; Geerts, S.; Vanuytrecht, E.; Abraha, B.; Deckers, J.; Bauer, H.; Gebrehiwot, K. Unravelling Crop Water Productivity of Tef (*Eragrostis tef* (Zucc.) Trotter) through AquaCrop in Northern Ethiopia. *Exp. Agric.* **2012**, *48*, 222–237. [CrossRef]

18. Tsegay, A.; Vanuytrecht, E.; Abrha, B.; Deckers, J.; Gebrehiwot, K.; Raes, D. Sowing and irrigation strategies for improving rainfed tef (*Eragrostis tef* (Zucc.) Trotter) production in the water scarce Tigray region, Ethiopia. *Agric. Water Manag.* **2015**, *150*, 81–91. [CrossRef]

19. Haileselassie, H.; Araya, A.; Habtu, S.; Meles, K.; Gebru, G.; Kisekka, I.; Girma, A.; Hadgu, K.; Foster, A. Exploring optimal farm resources management strategy for Quncho-teff (*Eragrostis tef* (Zucc.) Trotter) using AquaCrop model. *Agric. Water Manag.* **2016**, *178*, 148–158. [CrossRef]

20. Hoogenboom, G.; Porter, C.H.; Shelia, V.; Boote, K.J.; Singh, U.; White, J.W.; Hunt, L.A.; Ogoshi, R.; Lizaso, J.I.; Koo, J.; et al. *Decision Support System for Agrotechnology Transfer (DSSAT) Version 4.7*; DSSAT Foundation: Gainesville, FL, USA, 2017.

21. Jones, J.; Hoogenboom, G.; Porter, C.; Boote, K.; Batchelor, W.; Hunt, L.; Wilkens, P.; Singh, U.; Gijsman, A.; Ritchie, J. The DSSAT cropping system model. *Eur. J. Agron.* **2003**, *18*, 235–265. [CrossRef]

22. Kassie, B.T.; Asseng, S.; Porter, C.H.; Royce, F.S. Performance of DSSAT-Nwheat across a wide range of current and future growing conditions. *Eur. J. Agron.* **2016**, *81*, 27–36. [CrossRef]

23. Kebede, H.; Johnson, R.C.; Ferris, D.M. Photosynthetic Response of *Eragrostis-Tef* to Temperature. *Physiol. Plant.* **1989**, *77*, 262–266. [CrossRef]

24. Kumar, S.; Sharma, V.; Chaudhary, S.; Tyagi, A.; Mishra, P.; Priyadarshini, A.; Singh, A. Genetics of flowering time in bread wheat Triticum aestivum: Complementary interaction between vernalization-insensitive and photoperiod-insensitive mutations imparts very early flowering habit to spring wheat. *J. Genet.* **2012**, *91*, 33–47. [CrossRef]

25. Hocking, P.J.; Meyer, C.P. Effects of CO2 Enrichment and Nitrogen Stress on Growth, and Partitioning of Dry-Matter and Nitrogen in Wheat and Maize. *Aust. J. Plant Physiol.* **1991**, *18*, 339–356. [CrossRef]

26. Hammer, G.L.; Muchow, R.C. Quantifying climatic risk to sorghum in Australia's semi-arid tropics and subtropics: Model development and simulation. In *Climatic Risk in Crop Production: Models and Management for the Semi-Arid Tropics and Subtropics*; Muchow, R.C., Bellamy, J.A., Eds.; CAB International: Wallingford, UK, 1991; pp. 205–232.

27. Tulema, B.; Zapata, F.; Aune, J.; Sitaula, B. N. fertilisation, soil type and cultivars effects on N use efficiency in tef (*Eragrostis tef* (Zucc.) Trotter). *Nutr. Cycl. Agroecosyst.* **2005**, *71*, 203–211. [CrossRef]

28. Abraha, M.; Shimelis, H.; Laing, M.; Assefa, K.; Amelework, B. Assessment of the genetic relationship of tef (*Eragrostis tef*) genotypes using SSR markers. *S. Afr. J. Bot.* **2016**, *105*, 106–110. [CrossRef]
29. Davison, J.; Laca, M. *Grain Production of 15 Tef Varieties in Churchill County, Nevada During 2009*; University of Nevada Cooperative Extension: Reno, NV, USA, 2010.

30. Ayele, M.; Blum, A.; Nguyen, H. Diversity for osmotic adjustment and root depth in Tef (*Eragrostis tef* (Zucc) Trotter). *Euphytica* **2001**, *121*, 237–249. [CrossRef]

31. Okubay, G. Soil Fertility Characterization and Response of Wheat (*Triticum aestivum*) and Tef (*Eragrostis tef*) to Different Rates of Slow Release and Conventional Urea Fertilizer in Vertisols of Ofla District, Southern Tigray, Ethiopia. Master's Thesis, Haramaya University, Haramaya, Ethiopia, 2012.

32. Funk, C.; Peterson, P.; Landsfeld, M.; Pedreros, D.; Verdin, J.; Shukla, S.; Husak, G.; Rowland, J.; Harrison, L.; Hoell, A.; et al. The climate hazards infrared precipitation with stations—A new environmental record for monitoring extremes. *Sci. Data* **2015**, *2*, 150066. [CrossRef] [PubMed]

33. NASA (Ed.) *NASA Prediction of World Wide Energy Resource (POWER) Climatology Resource for Agroclimatology*; NASA Langley Research Center (LaRC): Hampton, VA, USA, 2017.

34. NOAA. *Fallon Experimental Station, NV, US.(NOAA)*; National Centers for Environmental Information: Asheville, NC, USA, 2018.

35. Tulema, B.; Aune, J.; Breland, T. Availability of organic nutrient sources and their effects on yield and nutrient recovery of tef (*Eragrostis tef* (Zucc.) Trotter) and on soil properties. *J. Plant Nutr. Soil Sci. Z. Fur Pflanzenernahr. Und Bodenkd.* **2007**, *170*, 543–550. [CrossRef]

36. USDA NRCS. *Web Soil Survey*, 9th September 2016 ed.; USDA NRCS: Lincoln, NE, USA, 2016.

37. Habtegebrial, K.; Sing, B.; Haile, M. Impact of tillage and nitrogen fertilization on yield, nitrogen use efficiency of tef (*Eragrostis tef* (Zucc.) Trotter) and soil properties. *Soil Tillage Res.* **2007**, *94*, 55–63. [CrossRef]

38. Rawls, W.; Brakensiek, D.; Saxton, K. Estimation of Soil-water properties. *Trans. Asae* **1982**, *25*, 1316–1320. [CrossRef]

39. Tulema, B.; Aune, J.; Johnsen, F.; Vanlauwe, B. The prospects of reduced tillage in tef (*Eragrostis tef* Zucca) in Gare Arera, West Shawa Zone of Oromiya, Ethiopia. *Soil Tillage Res.* **2008**, *99*, 58–65. [CrossRef]

40. Settie, A.; Tadesse, D.; Fekadu, F.; Seyfu, K. Seed Yield Stability and Adaptability to Different Soil Types of Some Improved Varieties of Tef. *Sebil* **1996**, *7*, 22–28.

41. Teshome, Y.; Verheye, W. Crop growth requirements and approach to yield prediction for barley and tef in the Central Highlands of Ethiopia. *Pedologie* **1993**, *43*, 357–372.

42. Tsegay, A. *Tef*; FAO: Rome, Italy, 2012; pp. 150–154.

43. Van Oosterom, E.; Carberry, P.; Hargreaves, J.; O'Leary, G. Simulating growth, development, and yield of tillering pearl millet II. Simulation of canopy development. *Field Crop. Res.* **2001**, *72*, 67–91. [CrossRef]

44. Liu, H.L.; Yang, J.Y.; He, P.; Bai, Y.L.; Jin, J.Y.; Drury, C.F.; Zhu, Y.P.; Yang, X.M.; Li, W.J.; Xie, J.G.; et al. Optimizing Parameters of CSM-CERES-Maize Model to Improve Simulation Performance of Maize Growth and Nitrogen Uptake in Northeast China. *J. Integr. Agric.* **2012**, *11*, 1898–1913. [CrossRef]

45. ICRISAT. *Modeling the Growth and Development of Sorghum and Pearl Millet*; ICRISAT: Patancheru, India, 1989.

46. White, J.W.; Alagarswamy, G.; Ottman, M.J.; Porter, C.H.; Singh, U.; Hoogenboom, G. An Overview of CERES-Sorghum as Implemented in the Cropping System Model Version 4.5. *Agron. J.* **2015**, *107*, 1987–2002. [CrossRef]

47. Alderman, P.D.; Quilligan, E.; Asseng, S.; Ewert, F.; Reynolds, M. Proceedings of the Workshop on Modeling Wheat Response to High Temperature. In Proceedings of the Workshop on Modeling Wheat Response to Higher Temperature, El Batan, Texcoco, Mexico, 19–21 June 2013.

48. Urban, D.W.; Sheffield, J.; Lobell, D.B. Historical effects of CO2 and climate trends on global crop water demand. *Nat. Clim. Chang.* **2017**, *7*, 901. [CrossRef]

49. Hammer, G.L.; Muchow, R.C. Assessing Climatic Risk to Sorghum Production in Water-Limited Subtropical Environments. 1. Development and Testing of a Simulation-Model. *Field Crop. Res.* **1994**, *36*, 221–234. [CrossRef]

50. Pembleton, K.G.; Cullen, B.R.; Rawnsley, R.P.; Harrison, M.T.; Ramilan, T. Modelling the resilience of forage crop production to future climate change in the dairy regions of Southeastern Australia using APSIM. *J. Agric. Sci.* **2016**, *154*, 1131–1152. [CrossRef]

51. Teixeira, E.I.; de Ruiter, J.; Ausseil, A.G.; Daigneault, A.; Johnstone, P.; Holmes, A.; Tait, A.; Ewert, F. Adapting crop rotations to climate change in regional impact modelling assessments. *Sci. Total Environ.* **2018**, *616*, 785–795. [CrossRef]
52. Dawit, T.; Hirut, K. *Germplasm Collection, Conservation and Characterization of Tef*; Plant Genetic Resources Centre/Ethiopia: Addis Abeba, Ethiopia, 1995.

53. Tadelle Tadesse, K. *Photoperiod Sensitivity of Tef* [Eragrostis tef (Zucc.) Trotter] *Cultivars*; Wageningen University: Wageningen, The Netherlands, 2005.

54. R Core Team. *R: A Language and Environment for Statistical Computing*; R Foundation for Statistical Computing: Vienna, Austria, 2019.

55. Mitchell, P. Misuse of regression for empirical validation of models. *Agric. Syst.* **1997**, *54*, 313–326. [CrossRef]

56. Akinseye, F.M.; Adam, M.; Agele, S.O.; Hoffmann, M.P.; Traore, P.C.S.; Whitbread, A.M. Assessing crop model improvements through comparison of sorghum (*sorghum bicolor* L. moench) simulation models: A case study of West African varieties. *Field Crop. Res.* **2017**, *201*, 19–31. [CrossRef]

<figure><img src="image_12.png" /></figure>

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).