Crop Yield and Physicochemical Properties of Wheat Grains as Affected by Tillage Systems

: In this study, we aimed to determine the effects of tillage systems (TS) on wheat grain yield, yield components, and physicochemical properties under rain-fed conditions over 2 years (2018–2019) in the Kuhdasht region in southwestern Iran. The tillage treatments were a combination of conventional tillage (CT), reduced tillage (RT), and no tillage (NT) systems based on randomized complete blocks with three replications. Wheat grain yields and yield components, physical properties (geometrical dimensions, sphericity, thousand-grain mass (TGM), water absorption (WA), volume (V), and density), chemical composition (protein content (PC), water content (WC), ﬁber content (FC), ash content (AC)), surface color parameters (L*, a*, and b*), and correlations among measured parameters were assessed. The results showed that wheat grain yield and yield components were signiﬁcantly affected by TS ( p < 0.01). The wheat grain yield was higher under CT (2.72 t ha − 1 ), with a signiﬁcant difference between RT (1.76 t ha − 1 ) and NT (1.20 t ha − 1 ). The highest TGM (49.19 g) was achieved under CT, followed by RT (43.41 g), with the lowest (39.17 g) in the NT system. TS had signiﬁcant effects on certain physical properties of wheat grains ( p < 0.01). Wheat grain, size, shape, and mass values were higher under CT than RT and NT systems in all experimental years. CT resulted in the highest WA, while the lowest resulted from NT. TS had no signiﬁcant inﬂuence on the density of the grains. TS had signiﬁcant effects on PC, AC, and WC values of grains at the 1% probability level and on FC at the 5% level. Higher PC values were obtained under the CT system (13.07%, against 11.90% and 10.67% for RT and NT, respectively) in all growing seasons. Additionally, the AC was signiﬁcantly lower in the sample grains under RT (2.38%) and NT (2.43%) than in those from CT (2.56). The FC was higher under NT (15.73%) than RT (13.73) and CT (13.71%). The grain WC was signiﬁcantly higher under NT (7.57%) than RT (6.79%) and CT (7.0%). The TS signiﬁcantly affected the surface color parameters of grains, while the L* (lightness) and b* (yellowness) values of grains under CT were higher than in RT and NT.


Introduction
Wheat is one of the most important cereal products, which has a significant effect on human nutrition. Wheat has been cultivated for more than 10,000 years as a primary human food and it has had a significant impact on human social, demographic, and cultural development [1]. In addition to rice, bread made from wheat grains is a leading human food. As with other grains, wheat grains are a major source of energy, carbohydrates (such as starch), protein, minerals, vitamins, and fiber. The protein content of wheat grains varies between 8 and 20%, depending on the type and class [1]. Wheat grain quality parameters heavily depend on growth conditions, soil fertility, fertilizer application, water ability, genotype, grain transport, and storage situations [2][3][4][5]. Some of the essential quality factors of grains such as wheat include the grain structure, grain cover, resistance to breakage (mechanical damage), and grain composition, which are affected by the rates of inorganic fertilizers [6,7]. Additionally, the yield and quality parameters of grains are affected by agronomic practices, such as tillage systems [8][9][10][11][12][13][14]. Grain size is one of the most important yield components for cereals such as wheat. The grain protein content is the main factor in terms of wheat quality [15]. Some studies have shown little correlation between wheat grain protein and soil tillage systems [16][17][18], while others have found positive effects and reported that CT systems tend to produce higher protein contents than other systems [19][20][21][22].
Color is an important factor in the food industry. The surface color of wheat grain is an important attribute in the sorting and marketing of the grains [23,24]. The color of wheat grains is the key to the milling and baking industry, which determines the color of the whole wheat products [24]. In the literature, several studies on the color measurement of wheat grains are available, in which the researchers reported that the color and appearance of the final product made from wheat grains, such as bread and tortillas, are influenced by the color of the wheat grains [25][26][27].
Soil tillage is one of the most important processes in the production of agricultural products. More than half of the total energy required for cereal production is consumed during soil tillage [15]. On the other hand, soil tillage is also an important factor affecting soil structure properties and crop yield [28][29][30]. Wheat production in Iran is based on the CT system, which includes the use of a moldboard plow and seed bed peroration using a disk harrow. The CT system is characterized by high energy consumption and has the most destructive effects on soil properties [31]; therefore, there is an ongoing effort to replace the CT system with conservation systems such as RT or NT systems. However, this effort may be limited by certain problems such as environmental conditions and economic factors [32,33]. As stated by many researchers, in the conventional tillage and conservation tillage systems, the yield and quality parameters of crops are functions of the interactions of many parameters, which are difficult to predict [34][35][36]. Generally, from the review of the results of research studies conducted in this paper, it can be stated that when reducing the levels of tillage from the conventional system to conservation tillage systems (RT or NT systems) in cereals, there is a decrease in crop yield [37]. The soil tillage system also affects the quality parameters of grains, such as the physical dimensions and chemical composition [34,38,39]. Lopez-Bellodo et al. [19] stated that higher protein content for wheat grains was obtained for the CT system than in the NT system. Wozniak et al. [40] reported that the herbicide tillage system (use of Roundup 360 SL herbicide to control weeds) reduced the volume weight and increased ash content in the wheat grains as compared to CT and RT systems. In addition to the cases mentioned above, the tillage system influences weed infestations. Many investigations show that the NT system increased the weed infestation of crops [41].
The effects of several tillage systems on the yield, yield components, and physicochemical attributes of wheat grains are required to achieve more stable wheat production. According to the needs mentioned above, the purpose of this study was to investigate the effects of tillage systems, namely conventional tillage (CT), reduced tillage (RT), and no tillage (NT) systems, on wheat grain yield, yield components, physical properties (grains geometrical dimensions, arithmetic mean diameter, geometric mean diameter, equivalent mean diameter, sphericity, 1000-grain mass, volume, water absorption capacity, and density), chemical composition (protein content, fiber content, ash content, and water content), and surface color parameters (whiteness, redness, and yellowness) of wheat grains over 2 experiment years. 16 • C, and with a semi-arid climate. The monthly rainfall during the experimental years is shown in Figure 1. The sum of annual rainfall during the growing season in 2017-2018 was 835 mm, which was much higher (about 2 times) than the long-term rainfall (450 mm). Additionally, the sum of annual rainfall during the 2018-2019 growing season was 526 mm, which was more than the long-term rainfall data (450 mm). Field experiments were conducted during 2017-2018 and 2018-2019 growing seasons in the Kohdasht region in the southwest of Iran (33°34′, 47°36′ E), at an elevation of 1195 m above sea level, with a mean annual precipitin of 450 mm, average yearly temperature of 16 °C, and with a semi-arid climate. The monthly rainfall during the experimental years is shown in Figure 1. The sum of annual rainfall during the growing season in 2017-2018 was 835 mm, which was much higher (about 2 times) than the long-term rainfall (450 mm). Additionally, the sum of annual rainfall during the 2018-2019 growing season was 526 mm, which was more than the long-term rainfall data (450 mm).

Experimental Procedure
The cropping system used in this experiment was dry land continuous cultivation in three tillage systems. Trials were established using a randomized complete block design over 2 years, with three tillage treatments and three replications. The tillage treatments were a combination of three systems: (1) conventional tillage (CT), (2) reduced tillage (RT), and (3) no-tillage (NT) in the autumn. Conventional tillage (CT) included using a moldboard plow (depth close to 25 cm) + disk harrow (depth close to 10 cm) + seed drilling using a conventional seed driller. Reduced tillage (RT) included using field cultivation (depth close to 10 cm) after harvest of the previous crop (chickpea) + seed drilling using a direct drilling machine. No-tillage (NT) included using a direct drilling machine after harvest of the previous crop (chickpea), without any prior tillage.
Wheat seeds of Aftab, a popular winter bread wheat variety in Iran, which is suitable for dry land cultivation, were chosen for use in this research and planted at a rate of 160 kg ha −1 (the rate for rain-fed conditions). The row spacing (11 rows in each plot) was 0.2 m. Individual tillage plots were 15 m in length and 30 m in width (450 m 2 plot size), which were divided into three sub-blocks. The wheat cultivar and planting dates were chosen based on the best recommendations for farmers in the area of the study. Crop management exercises (except tillage methods) were diverse among locations depending on the requirements based on soil test results. Fertilizers were applied at the planting date (100 kg ha −1 P2O5 and 75 kg ha −1 N). Herbicides were used for weed control by application of 20 g ha −1 of glyphosate and 0.8 L ha −1 of clodinafop-propargyl in spring.

Soil Sampling
Soil samples from four locations were collected from each main plot at 0-30 cm soil depth after harvest of the previous crop (chickpea) to determine the properties of the studied field soil [42]. Soil bulk density was determined by taking soil samples at 0-10

Experimental Procedure
The cropping system used in this experiment was dry land continuous cultivation in three tillage systems. Trials were established using a randomized complete block design over 2 years, with three tillage treatments and three replications. The tillage treatments were a combination of three systems: (1) conventional tillage (CT), (2) reduced tillage (RT), and (3) no-tillage (NT) in the autumn. Conventional tillage (CT) included using a moldboard plow (depth close to 25 cm) + disk harrow (depth close to 10 cm) + seed drilling using a conventional seed driller. Reduced tillage (RT) included using field cultivation (depth close to 10 cm) after harvest of the previous crop (chickpea) + seed drilling using a direct drilling machine. No-tillage (NT) included using a direct drilling machine after harvest of the previous crop (chickpea), without any prior tillage.
Wheat seeds of Aftab, a popular winter bread wheat variety in Iran, which is suitable for dry land cultivation, were chosen for use in this research and planted at a rate of 160 kg ha −1 (the rate for rain-fed conditions). The row spacing (11 rows in each plot) was 0.2 m. Individual tillage plots were 15 m in length and 30 m in width (450 m 2 plot size), which were divided into three sub-blocks. The wheat cultivar and planting dates were chosen based on the best recommendations for farmers in the area of the study. Crop management exercises (except tillage methods) were diverse among locations depending on the requirements based on soil test results. Fertilizers were applied at the planting date (100 kg ha −1 P 2 O 5 and 75 kg ha −1 N). Herbicides were used for weed control by application of 20 g ha −1 of glyphosate and 0.8 L ha −1 of clodinafop-propargyl in spring.

Soil Sampling
Soil samples from four locations were collected from each main plot at 0-30 cm soil depth after harvest of the previous crop (chickpea) to determine the properties of the studied field soil [42]. Soil bulk density was determined by taking soil samples at 0-10 and 11-20 cm soil depths using a standard core sampler (4 cm long and 8 cm in diameter) and oven-dried at 105 • C for 48 h, then the bulk density was calculated by dividing the sample dry weight to the sample's total volume [43]. Soil-available P was determined by extracting samples with 0.5 mol L −1 NaHCO 3 and by determining P calorimetrically [44].
Soil-available K was measured using the 1 mol L −1 ammonium acetate extraction-flame photometry method [45]. Table 1 shows the properties of the studied field soil.

Crop Sampling
After physiological maturity, all plants in the central area of plots were cut above soil using a hand sickle. The edges (0.5 m) of the plots were cut and discarded before harvest. The plant samples were air-dried, threshed, and cleaned in an air screen cleaner to obtain grain yield and yield component values. The grain yield per hectare, number of plants per unit area (m 2 ), plant height, spike length, number of grains per spike, spike mass, and thousand-grain mass were measured and examined. The number of plants per unit area of plots was determined using a wooden frame with 50 × 350 cm over three replications. Plant height and spike length values were measured using a ruler. The spike mass was measured using a digital scale with an accuracy of 0.01 g. The thousand-grain mass (TGM) of each grain sample treatment was measured by counting 1000 grains and weighing them in an electronic balance with a precision of 0.001 g.
The quality properties of the grains are affected by the moisture content [46]. Harvested grain samples were stored for six months (at 4 • C and a humidity of about 60%) to standard and equalize their moisture content and determine the effects of the tillage system on the water content (WC) (stored moisture for six months) of grains [47,48], then their quality properties were measured. To investigate the effects of the year of study and tillage systems on the physical properties of wheat grains, including grain size, 100 grains randomly selected from each sample and their three linear dimensions (namely length (L), width (W), and thickness (T)) were measured using a digital micrometer over three replications, with an accuracy of 0.01 mm. The arithmetic mean diameter (D a ), geometric mean diameter (D g ), equivalent mean diameter (D P ), and sphericity (ϕ), respectively, of the grain samples were computed using the following relationships [49]: To measure the density of tested samples, firstly 100 grains of each treatment were weighed, then their volume (V) was measured using the liquid displacement method (toluene). The density of the grains was computed by dividing the mass of the sample by the actual volume. To measure the water absorption (WA) capacity of grain samples, 50 whole grains of each sample were weighed and soaked in 100 mL of distilled water in a cylindrical container covered with foil, which were left at room temperature (22 • C) for 24 h. The grains were then dried using absorbent paper and weighed. The water absorption capacity of each sample was computed as (weight after soaking − weight before soaking)/50. The chemical compositions of whole wheat grains, including the protein content (PC), fiber content (FC) (natural detergent fiber (NDF)), ash content (AC) and water content (WC) (stored moisture for six months), were computed using near-infrared reflectance (NIR) technique using a Perten DA 7250 NIR analyzer instrument (Perten-Instruments, Ab, Huddinge, Sweden) [5].
For surface color parameter measurements, glass dishes containing wheat grain samples were placed under 4 LED lamps (12 W, 230 V, AC, and 50 Hz). The images of whole grain surfaces were captured using a 24.5 megapixel digital color camera (Nikon digital SLR camera D53000). The camera was positioned 20 cm vertically above the grain sample surface before capturing pictures with no flash or zoom. The color properties of captured images were analyzed using Image J software version 1.32 g. Finally, for quantization of the grain surface color, L*, a*, and b* coordinates of the CIE color scale were calculated from RGB values from digital images of the grains using a ColorMine color converter

Data Analysis
Analysis of variance (ANOVA) was performed to determine the effects of different treatments (year of study and tillage system on wheat grain yield and properties. All statistical analyses were carried out using SPSS software version 19 (SPSS Institute Inc., Chicago, IL, USA). The levels of significance were shown as * p < 0.05 and ** p < 0.01 using Duncan's multiple range tests.

Wheat Grain Yield and Yield Components
Wheat grain yield and yield components were significantly affected by tillage systems (p < 0.01). There were no significant effects for the year of research or the interaction between the year and tillage system on the wheat grain yield or yield components (p > 0.05) ( Table 2). This may be due to the similar average rainfall observed in both experimental years, which exceeded the long-term rainfall data (450 mm) ( Figure 1).  Table 2 shows the results of the analysis of variance (F-Values) and the comparison of means for the wheat grain yield and yield components in different tillage systems. Conservation tillage systems (especially the NT system) were shown to reduce the wheat grain yield and yield components ( Table 2). The results of many studies have demonstrated that better wheat grain yield and quality parameters have been achieved in CT system in moderately humid soils, in NT system in dry soils, or in the years with low rainfall [19,21]. The wheat grain yield was higher under CT with a significant difference between RT and NT ( Table 2). The poor grain yield under NT may be related to the higher soil bulk density and lower soil moisture in root zones under NT [55]. It seems that the reasons for the decrease in wheat grain yield due to the decrease in tillage operations are the decrease in the number of plants per m 2 , decrease in the number of grains per spike, and also the increase in weeds. Similar results were reported for the reduction of wheat grain yield due to the reduction of tillage operations [35]. In contrast, some researchers [56][57][58] have reported increases in wheat grain yield under the NT system.
With  (Table 2). Wheat plant height and spike length values were significantly greater under CT than RT or NT. The CT system produced the maximum spike mass, while the minimum spike mass observed under NT and the intermediate amount was observed under RT. The number of grains per spike also decreased with decreasing tillage levels. The highest number of grains per spike was related to the conventional tillage system, with no significant difference with RT, while the lowest number was related to NT. Zamani et al. [59] reported a higher number of grains per spike for wheat grain under the CT system.
Wheat thousand-grain mass (TGM) was significantly influenced by the tillage system (F = 503.65). The sowing of wheat seeds in the CT and RT systems produced grains with higher mass values than wheat grown in the NT system. Wheat thousand-grain mass values reduced by about 20% as the level of tillage reduced from CT to the NT system. The highest mean thousand-grain mass value was achieved in the CT system, while the value was lower in the RT and the lowest in the NT system ( Table 2). The lower TGM values under RT and NT systems may have been related to soil surface compaction, less rainwater penetration into the soil, and water stress at the grain filling stage, as both treatments received lower moisture than in CT. Similar results have been reported regarding reduced TGM values with decreasing tillage level by Wozniak et al. [40] and Wozniak [37] for wheat, by Wozniak and Soroka [35] for triticale, and by Wozniak and Gaweda [33] for pea plants. Wozniak et al. [40] reported a significant difference of TGM for wheat grains in tillage treatments with mean values of 42.3 g, 41.5 g, and 40.4 g for CT, RT, and NT systems, respectively. Tabatabaeefar et al. [60] reported TGM mean values of 38.81 g for wheat grains under the CT system and 38.2 g under NT treatments.

Wheat Grain Physical Properties
The results of the study illustrated that the physical properties of wheat grains hinged on the tillage system. The year of research and the interaction between the year and tillage system had no significant effect on the wheat grains' physical properties (p > 0.05). This may be due to the similarly high annual rainfall levels in both experimental years, which were greater than the long-term rainfall data (450 mm) (Figure 1). Table 3 shows the results of the analysis of variance (F-Values) and mean comparison (average for 2 experimental years) for the physical properties of wheat grains in different tillage systems. The tillage systems appeared to have significant effects (p < 0.01) on all physical properties of wheat grains except grain density. Conservation tillage systems (especially the NT system) were shown to reduce the physical properties of wheat grains (Table 3).  Wheat grains dimensions include length (L), width (W), thickness (T) and following them arithmetic mean diameter (D a ), geometric mean diameter (D g ), and equivalent mean diameter (D P ) of the grains were significantly affected by tillage system (p < 0.01) ( Table 2). Wheat grains in the CT system had the highest L value, while grains in the NT had the lowest. The maximum grain width was related to the RT system, while the lowest value was for the NT system. The highest and lowest grain thickness (T) values were 2.90 and 2.63 mm, which were obtained in the RT and NT systems, respectively. Wheat grains in the CT system had the lowest arithmetic mean diameter (D a ) mean value, while grains in the NT system had the highest. The geometric mean (D g ) of grains was lower with no tillage than with reduced and conventional tillage treatments. Wheat grains in the RT system had the highest equivalent mean diameter (D P ) value and grains in the NT system had the lowest. Additionally, the maximum equivalent mean diameter (D P ) of the grains was related to the RT system and the lowest value was for the NT. Wheat grain sphericity (ϕ) was significantly affected by the tillage system (p < 0.01) ( Table 3). Wheat grains in the RT system had the highest sphericity value (0.60%) and grains in the CT and NT systems had the lowest values (0.58%). According to the above results, which are shown in Table 3, almost all of the maximum values for the physical dimensions of wheat grains were for the CT and RT systems and the lowest values were for the NT, showing that the RT system provides better results in terms of impacts on the physical dimensions of wheat compared to other systems. The results are exciting, however unfortunately no research report on the effects of tillage systems on the physical dimensions of wheat grain can be found in the literature to compare with the results obtained in this research. Taner et al. [15] reported the highest particle size index value for wheat grains with the RT system and the lowest with the NT system. Javadi et al. [61] recommended RT for wheat production in Iran.
Wheat grain volume (V) was significantly affected by the tillage system (p < 0.01) ( Table 3). Wheat grains in the CT system had the highest volume and grains in the NT system had the lowest. Additionally, with the reduction of the tillage level from CT to the NT system, the water absorption capacity of wheat grains was decreased by about 23% (from 0.023 g to 0.018 g). This reduction in the water absorption capacity may be explained by the decrease of the grains' geometric dimensions and volume by reducing the tillage level from the CT to the NT system. The density of wheat grains was approximately the same (about 1.3 g mm −3 ) for all tillage systems used in the present study.

Wheat Grain Chemical Compositions
Different chemical compositions of wheat grains were measured and compared in different tillage systems over 2 experimental years, including the protein content, fiber content, ash content, and water content. Table 4 shows the results of the analysis of variance (p-Values) and comparison means for the chemical compositions of wheat grains in different tillage systems. From the data in Table 4, it is evident that all of the measured chemical compositions of wheat grains were significantly influenced by the tillage systems. As seen from the data in Table 4, the experiment year, tillage system, and the interaction between the year and tillage system had significant influences on the wheat grains' protein content at the 1% probability level (p < 0.01). The F value of the analysis of variance demonstrated that wheat grains' protein content was determined to a great extent by the tillage system (F = 259.267), more so than the experimental year (F = 64.067) ( Table 4). The protein content was significantly higher in the sample grains harvested from the plots under the CT system than from RT and NT systems. According to the data in Table 4, the mean values for the protein content were 13.07, 11.90, and 10.67% for grains from the CT, RT, and NT systems, respectively, which were statically different from each other based on Duncan's multiple range tests ( Table 4). The protein content value for the CT system met the industry standard (PC ≥ 12.5%) [21]. These observations are in agreement with other studies that have reported significant decreases in protein content in grains under conservation tillage systems such as RT and NT [11,21,22,38,62,63]. Lopez-Bellido et al. [19,20] reported higher protein content under CT than NT. Additionally, Pringas and Koch [64] reported that the protein content and sedimentation value decreased significantly without tillage compared to the tillage-based system.
As seen from the data in Table 4, the year and interaction between the year and tillage system had no significant influence on the wheat grains' fiber content (p > 0.05); however, the tillage system had a significant effect on this property of at the 5% probability level (p < 0.05) ( Table 4). From the data in Table 4, it is evident that the fiber content of wheat grains increased significantly from 13.71 to 15.73% by decreasing the tillage level from the CT system to the NT system. The mean values of the fiber content of wheat grains were found to be 13.71, 13.73, and 15.73% for the CT, RT, and NT systems, respectively. Moreover, based on Duncan's multiple range tests, the mean value of the fiber content of grains in the NT system was statically different (higher) from the fiber contents of grains in the CT and RT systems. Wozniak and Makarski [38] reported that tillage systems affected the chemical compositions of wheat grain, such as the fiber content. Additionally, Wozniak [11] reported that the fiber content of spring triticale grains was higher for NT than for RT and CT systems.
The effects of the experimental year, tillage system, and the interaction between the year and tillage system on the ash content of harvested wheat grans were significant at Sustainability 2021, 13, 4781 9 of 15 the 1% probability level (p < 0.01) ( Table 4). The ash content of grains was significantly lower in the sample grains harvested from the plots under RT and NT systems than from the CT system. According to the data in Table 4, the mean values for the ash content were 2.56, 2.38, and 2.43% for grains in the CT, RT, and NT systems, respectively, while the ash contents of the grains in the RT and NT systems were statically different (lower) than that of the CT system based on Duncan's multiply range tests (Table 4). A similar result was reported by Wozniak (2016), who reported that the ash content of spring triticale grains was higher for CT than for RT and NT systems.In contrast, Kraska [65] and Wozniak and Makarski [38] reported higher ash content in wheat grains cultivated in a plowless system rather than in the conventional system.
The effect of the tillage system on the water content (stored water in grains over six months) of wheat grains was significant at the 1% probability level (F = 43.86). However, the effects of the experimental year and the interaction between the year and tillage system on the water content of harvested wheat grins were not significant. The results of the mean comparison in Table 4 indicate that by reducing the level of soil tillage from CT to NT, the average values of the water content of grains significantly increase from 7.0% to 7.57%. The mean values of the water content of grains were 7.0, 6.79, and 7.57% for CT, RT, and NT tillage systems, respectively. Moreover, based on Duncan's multiple range tests, the mean value of the water content of grains in the NT system was statically different (higher) than in the CT and RT systems. No statistical difference was found for the water content of grains in the CT and RT systems (Table 4). No reported research result was found in the published literature for the effects of tillage systems on the water content of wheat grains for comparison with the results of this research.
The protein and ash contents of grains and the interactions between the tillage systems and the years of study are shown in Table 5. The protein and ash contents of grains were higher in the CT and RT systems compared to the NT system in both experimental years (Table 5). Additionally, higher protein and ash contents were achieved in the first year (2018) of the study than in the second year (2019). This may be due to higher rainfall in the 2018 growing season.

Wheat Grain surface Color Properties
The results of the data analysis (ANOVA) indicated that the tillage system significantly influenced all of the surface color parameters of wheat grains, including L* (whiteness), a* (redness), and b* (yellowness), at the 1% portability level. The experimental year and the interaction between the year and tillage system had no significant influence on the wheat grain's surface color parameters (p > 0.05) ( Table 6). The impact of the tillage system on the b* (F = 825.48) color value of wheat grains was more than for the L* (F = 32.49) and a* (F = 25.30) parameters within the range of variables studied (Table 6). From the mean comparison results for surface color parameters of wheat grains under the different tillage systems presented in Table 6, it is evident that the whiteness (L*) of the grains increased when increasing the level of tillage, thus changing the surface color of the grains toward a whiter (lighter) hue. These results show that the tillage system selection, which has considerable effects on the physicochemical properties of grains, also has significant effects on other grain's quality properties, such as the surface color. As the level of tillage increased from the NT to CT system, the mean values of the L* (whiteness) in wheat grains increased significantly from 59.33 to 60.15. No study was found in the literature for comparison of the effects of tillage system selection on the color parameters of wheat grains in this study.
From the data in Table 6, it is evident that the mean a* (redness) values of wheat grains ranged from 4.62 (in CT system) to 5.12 (in NT system). The mean a* values of the tested grains for CT, RT, and NT systems were 4.62, 4.76 and 5.12, respectively, while the mean value of a* in the NT system was statically different (higher) than in CT and RT systems ( Table 6). This result shows that as the level of tillage decreased from the CT to NT system, the redness color parameter of the wheat grain surfaces increased.
For the b* (yellowness) surface color scale of the wheat grains, the result showed that as the level of soil tillage decreased from the CT to NT system, the yellowness color parameter of the wheat grains decreased. From the data in Table 5, it is evident that the mean values of the b* ranged from 39.60 (in the NT system) to 44.65 (in the CT system). The mean values of the b* surface color scale of the tested grains were 44.65, 41.00, and 39.60 for CT, RT, and NT systems, respectively, which were statically different from each other based on Duncan's multiple range tests (Table 6).
No studies in the published literature reported results for the relationships between wheat grain surface color parameters and tillage variants and their effects on the surface color parameters of wheat grain for comparison with the results obtained in this study; however, several researchers have reported the relationship between the color of wheat grains and other agronomic factors, such as the relationships between the pre-harvest sprouting and wheat grain color [66] and between the genotype-environment interaction effect for grain color in hard white spring wheat [67]. McCaig et al. [68] reported a significant wheat kernel color change associated with field weathering. Lukow et al. [69] reported a significant effect of the environment on the wheat grain color. Horvath and Veha [70] investigated the dependence of wheat grain color characteristics on the grain size. Tanska et al. [71] investigated wheat grain color in relation to the cultivation regime.

Correlation Matrix
The Pearson correlation coefficients between the quality parameters of the wheat grains under the different experimental conditions (two experimental years and three tillage systems) are shown in Table 7. From the data in Table 7, it is clear that wheat grain yield and yield components are significantly and positively correlated. Additionally, the grain yield and yield components of wheat have significant and positive correlations with all of the physical dimensions of wheat grains except sphericity and density. Wheat grain yield and yield components have significant and positive correlations with the protein content; however, they have significant and negative correlations with fiber and water contents. Grain yield and thousand-grain mass components have significant and positive correlations with the ash content; however, the correlations of other yield components with ash content were not significant. Wheat grain yield and yield components have significant and negative correlations with all of the grain surface color properties ( Table 7). The physical dimensions of wheat grains except sphericity and density have significant and positive correlations together. Additionally, the physical dimensions of wheat grains have significant and positive correlations with the chemical composition, including protein and ash contents, but have significant and negative correlations with fiber and water contents. The wheat grain surface parameters, including whiteness (L*), redness (a*), and yellowness (b*), can be seen to have significant and negative correlations with the physical dimensions. Additionally, the surface color parameters of grains have significant and negative correlations with the chemical composition, including protein and ash contents; however, they have significant and positive correlations with fiber and water contents.

Conclusions
The results of this study demonstrated that the physicochemical properties of wheat grain depended on the tillage system. The no tillage (NT) system was shown to reduce the physicochemical properties of grains. Lower grain size and shape were shown with NT versus RT and CT systems. Additionally, lower grain protein and ash contents and higher fiber and water contents were obtained from the NT system. Grains from the CT system had higher whiteness and yellowness, while those from the NT system had higher red surface coloring. In conclusion, the RT system produced wheat grains with similar or slightly lower-quality properties to the CT system, so this system could be presented as an even-handed replacement for conventional tillage in areas with similar environmental conditions to the climate and site conditions in this study.