Growth Performance of Photoperiod-Sensitive Rice ( Oryza sativa L.) Varieties in Different Soil Types under Rainfed Condition in Cambodia

: In Cambodia, rice is predominantly produced in areas with rainfed lowland conditions where photoperiod-sensitive varieties are cultivated. A number of varieties have been released for rainfed lowland areas, and the rice grain yield has reportedly increased by approximately 12% over the past ten years. Moreover, great ﬂuctuations in yield performance have been observed across different soil types of the rainfed ecosystems of Cambodia. Therefore, the present study aimed to analyze the grain yield and stability among ten popular varieties that were released for rainfed lowland ecosystems across the four different soil types in Cambodia in two years. The grain yield varied 566 g m − 2 as the highest in clay soil and about 220 g m − 2 as the lowest in sandy soil. A combined ANOVA revealed signiﬁcant differences for the main effect of genotype, environment, and genotype-by-environment interaction (GEI) for all yield-related traits and grain yield per square meter. The principal component test results showed that the heterogeneity of grain yield was mainly attributable to the effect of environment, followed by the effect of genotype. In fertile conditions, a higher percentage of ﬁlled grains was supported by higher leaf N until the late stage with a wider ﬂag leaf. In conditions of moderate fertility, larger numbers of panicles were supported by a higher percentage of productive culms with higher leaf N until the late stage. In conditions of poor fertility, a higher percentage of ﬁlled grains was supported by higher leaf N until the late stage, which is considered to be important for higher grain yield. The variety Phka Rumduol showed these preferable traits and produced higher yields in fertile to poor natural soil fertility conditions with moderate variation. This variety is considered to be more desirable and ideal due to its stability and higher grain yield. The other varieties, namely, Phka Mealdei, Phka Rumdeng, and CAR4, were identiﬁed as above-average yielders. Therefore, those varieties potentially may be recommended for cultivation in rainfed lowland rice ecosystems in Cambodia due to their high yields. CAR4 showed moderate variation at the same level as Phka Rumduol. From the point of stability, Phka Rumduol and CAR4 can be expected to excel.


Introduction
Under changing climate conditions, the improvement of varieties previously released is a major challenge. Understanding the yield-determining factors and yield stability

Trial in 2019
The METs were conducted in eight rainfed lowland environments during the wet season from July to December, 2019 at two different locations in each of four provinces-Battambang (BTB), Pursat (PUR), Kampong Thom (KPT), and Siem Reap (SRP)-where different soil types were distributed as follows: (1) Tuol Samroung (light gray to brown, clayey topsoil with big cracks when dry), (2) Krakor (light gray to brown, loamy topsoil over loamy subsoil), (3) Prateah Lang (sandy topsoil with a bright color over loamy subsoil), and (4) Prey Khmer (very sandy soil with a bright color). The eight experimental sites represented the most important rainfed paddy rice ecosystems surrounding the Great Lake of Cambodia.
Topsoil layers 0-20 cm and 20-40 cm in depth were collected from the eight experimental sites prior to cultivation in 2019 and their soil textures and chemical properties were analyzed; the detailed soil properties of the eight test environments are elaborated in Table 1. Soil analysis was conducted at a laboratory of the General Directorate of Agriculture, Ministry of Agriculture, Forestry and Fisheries. We followed the standard protocols in Dobermann and Fairhurst (2000) [27] for the following parameters: pH (1:1 v/v soil-H 2 O mixture), total soil C, total soil N, total soil P, cation-exchange capacity (CEC), exchangeable K, and soil texture. Daily rainfall and air temperature were recorded at the meteorological station in each province. Soil pH was measured at a soil:1 N HCl ratio of 1:2.5 (w v −1 ) by a pH meter. The total nitrogen and carbon concentrations were measured by the Pregl-Dumas Method [28]. The determination of total phosphorus was analyzed considering the former description [29]. The exchangeable K + , Ca 2+ , Mg 2+ , and Na + were measured by the ammonium acetate saturation method [30]. The cation exchange capacity was calculated from the amounts of K + , Ca 2+ , Mg 2+ , and Na + measured in the samples. The ten rice genotypes were evaluated using a randomized complete block (RCB) design with four replications at each experimental site. Eighteen-day-old seedlings of each genotype were transplanted in a 10 m 2 plot at a density of 250 hills per plot (20 cm × 20 cm: 1 seedling/hill) in July of 2019. The specific fertilization rate for each soil type of the eight environments was used based on the recommended doses proposed by CARDI; the recommended rates of fertilizer application (N:P 2 O 5 :K 2 O kg ha −1 ) were 64:23:20, 118:23:0, 57:23:30, and 65:25:25 for the Tuol Samroung, Krakor, Prateah Lang, and Prey Khmer soil groups, respectively. A compound fertilizer of urea, DAP, and KCl was applied as the source of N, P 2 O 5 , and K 2 O at all experiment sites. Briefly, all DAP and KCl applications were done basally, while urea was split for two times, of which half was applied basally and the other half was evenly split and used as top-dressing fertilizer at the early tillering and panicle initiation stages.
The precipitation and air temperature data during the rice cultivation season from June to December 2019 over the past five years were collected from the meteorological station in each province.

Trial in 2020
The experiment was conducted in four rainfed lowland environments in the wet season from July to December, 2020, on location in each of the same four provinces as in the experiment in the previous year, namely, Battambang (BTB: BTB1 in 2019), Pursat (PUR: PUR1 in 2019), Kampong Thom (KPT: KPT1 in 2019), and Siem Reap (SRP: SRP2 in 2019). The methods of fertilizer application management and cultivation were the same as those in the experiment in 2019. The data of daily mean air temperature and monthly rainfall were recorded at the provincial meteorological station in each province.

Measurements 2.3.1. Trial in 2019
Data collection was conducted based on the Standard Evaluation System (SES) for rice of IRRI. Five plants were randomly selected from each replication for measuring the number of stems. Then, five plants were randomly sampled from each replication for measuring plant length, panicle number per hill, spikelet number per panicle, and percentage of filled grains. The panicle length of the longest culm in each hill was measured. Growth duration was recorded when the panicles in each plot turned brown. To measure dry matter weight as biomass, the shoot dry weights (above-ground biomasses) of five random plants were collected at the full-heading stage and the physiological maturity stage and were subjected to an oven-drying process at 70 • C until the constant weight was obtained.
To estimate the grain yield per square meter, 250 plants were harvested in a 10 m 2 plot and determined based on the basis of a standardized 14% moisture content. The grain yield per square meter along with agronomic traits such as growth duration, plant height, panicle number per hill, spikelet number per panicle, percentage of filled grains, 1000-grain weight (14% moisture content basis), and harvest index (grain yield/(grain yield + straw biomass)) were recorded.

Trial in 2020
Five plants were randomly selected at first from each replication and continuously used for measuring the number of stems and plant length every two weeks. Shoot dry matter weight, panicle length, yield components, and grain yield were measured by the same methods as used in the experiment in the previous year. The flag leaf on the main stem was used for measuring the length and width of the leaf blade. The nitrogen concentration in white rice was measured by the Kjeldahl method.

Statistical Analysis
All data were processed using Microsoft Excel 2019, and a statistical analysis was performed with SPSS (IBM SPSS Statistics, Chicago, IL, USA). The correlation analysis between the soil phyco-chemical parameters or between the agronomic parameters was used to clarify the important factors relating to the difference in the natural soil fertility. An analysis of variance (ANOVA) was conducted to evaluate the variation of genotypes, environments, and genotype-by-environment interaction across the eight environments. Tukey's honestly significant difference test was used to clarify the differences in the grain yield and yield components between the varieties or environments. The principal component analysis was employed to summarize explanatory variables such as agronomic parameters and environmental factors into fewer indicators and to identify the contribution of each parameter or factor to the difference in the grain yield. Pearson product-moment correlation analysis for the grain yield was performed between the locations and provinces to clarify the differences in grain yield between the varieties in different environments. Multiple regression analysis was conducted by the step-wise method to estimate the most important yield components.

Soil Properties at the Experimental Sites
This study was conducted at eight sites (two locations in each of four provinces where different types of soil are distributed). Therefore, soil physical and chemical parameters were compared among the experimental sites. Soil characteristics at the experimental sites differed greatly among the provinces (Table 2). The pH, total C, total N, total P, CEC, and exchangeable K were mostly highest in the Toul Samroung soil at the BTB1 or BTB2 site, and tended to be low in SRP, especially SRP2. The variation in soil texture between the environments appeared large in clay. When exploring the relationship between soil chemical properties and soil texture, a significant positive relationship (p < 0.01) was found between the clay content and the total C, total N, total P, CEC, and exchangeable K (Figure 1). These results suggest that soil with higher clay content tended to be more fertile than soil with lower clay content.   [19] reported that the sand content was negatively correlated with the clay content, CEC, and available K in Battambang, Kampong Thom, Pursat, and Siem Reap. They stated that a high sand content decreased the soil's ability to retain nutrients in Cambodia [19]. Their finding that rice yield in unfertilized plots strongly reflected the inherent soil fertility or nutrient-holding capacity agrees with the findings of previous studies [31,32]. According to White et al. (2000) [33], the Toul Samroung, Prateah Lang, and Prey Khmer soils are classified as Endoaqualfs (loamy/clayey topsoil, clayey subsoil), Plinthustalfs (sandy topsoil, loamy/clayey subsoil), and Psamments (sandy topsoil, sandy subsoil), respectively, in the USDA classification Krakor(loamy/clayey topsoil, loamy/clayey subsoil). Toul Samroung soil proved to be the most fertile of the three soil groups used in this study, probably because of its high clay and CEC [33]. In the case of our current study, the clay content was positively correlated with the total C, total N, total P, CEC, and exchangeable K, as described above. These results are considered to have a tendency similar to those of several former reports. Natural soil fertility was comparatively high in Tuol Samroung soil in Battambang and comparatively low in Prey Khmer soil in Siem Reap. The pH, total C, total N, total P, CEC, and exchangeable K were mostly highest in the Toul Samroung soil at the BTB1 or BTB2 site, and tended to be low in SRP, especially SRP2. The variation in soil texture between the environments appeared large in clay. When exploring the relationship between soil chemical properties and soil texture, a significant positive relationship (p<0.01) was found between the clay content and the total C, total N, total P, CEC, and exchangeable K ( Figure 1). These results suggest that soil with higher clay content tended to be more fertile than soil with lower clay content.   [19] reported that the sand content was negatively correlated with the clay content, CEC, and available K in Battambang, Kampong Thom, Pursat, and Siem Reap. They stated that a high sand content decreased the soil's ability to retain nutrients in Cambodia [19]. Their finding that rice yield in unfertilized plots strongly reflected the inherent soil fertility or nutrient-holding capacity agrees with the findings of previous studies [31,32]. According to White et al. (2000) [33], the Toul Samroung, Prateah Lang, and Prey Khmer soils are classified as Endoaqualfs (loamy/clayey topsoil, clayey subsoil), Plinthustalfs (sandy topsoil, loamy/clayey subsoil), and Psamments (sandy topsoil, sandy subsoil), respectively, in the USDA classification Krakor (loamy/clayey topsoil, loamy/clayey subsoil). Toul Samroung soil proved to be the most fertile of the three soil groups used in this study, probably because of its high clay and CEC [33]. In the case of our current study, the clay content was positively correlated with the total C, total N, total P, CEC, and exchangeable K, as described above. These results are considered to have a tendency similar to those of several former reports. Natural soil fertility was comparatively high in Tuol Samroung soil in Battambang and comparatively low in Prey Khmer soil in Siem Reap.

Climatic Conditions of the Experimental Sites in 2019
The four provinces are located west, south, east, and north of Tonle Sap Lake, where the climate is considered to be varied. Thus, the recorded daily mean air temperature and monthly rainfall are shown in Figure 2 and compared. The daily mean air temperature tended to be lower in Siem Reap (SRP) at the later part of growth (during the reproductive stage after October) than in the other provinces in 2019 as compared with the average for the most recent five years. The air temperature tended to be high in Battambang (BTB) in 2019. It was observed that the highest monthly rainfall was recorded during the vegetative stage in September rather than by the late growth stage. There was more precipitation in 2019 in September in Siem Reap (SRP) and less in all the provinces after October as compared with the average over the most recent five years.
to year and within years [19]. While rice is grown mostly during the wet season, there is frequent intermittent drought [34]. The area of rice cultivation in Cambodia has almost tripled from 1984 (0.978 × 10 6 ha) to 2014 (2.859 × 10 6 ha), and it continued to increase by 2.57% per year until 2020 [35]. Although the grain yield level in average has increased from 1.29 t ha −1 (1984) to 3.39 t ha −1 (2020), the yield of rainfed lowland rice in the wet season remains lower than 3.0 t ha −1 [19]. The variation in the amount of precipitation between provinces in 2019 appears larger compared with the average for five years due to the large value in Siem Reap in 2019 as opposed to the average in August and September and the smaller value in all provinces from October to December in 2019.

Yield-Related Traits in 2019
The grain yield and its components have large variation in the rainfed rice varieties grown at different environments. The grain yield per unit land area and its yield components are initially listed in Table 3 for comparison. The yield components and yield of the ten photoperiod-sensitive varieties across the eight environments were investigated ( Table  3). The ANOVA results detected a highly significant difference for the main effects of genotype and environment on yield components; moreover, the genotype-by-environment interaction effect was significant for all these components. As shown in Table 3, the heaviest mean grain yield of 449.8 g m −2 was recorded in Phka Rumduol (G1), followed by Phka Mealdei (G4), yielding 430.4 g m −2 , while the lowest mean grain yield of 309.1 g m −2 was observed in CAR11. The highest mean panicle number per hill and percentage of filled grains were observed in Phka Rumduol (G1), followed by Phka Mealdei (G4). The mean spikelet number per panicle was largest in Raing Chey (G6), and the mean 1000-grain The climate in this region is tropical monsoonal, with a wet season (from June to November) followed by a prolonged dry season, and rainfall is irregular both from year to year and within years [19]. While rice is grown mostly during the wet season, there is frequent intermittent drought [34]. The area of rice cultivation in Cambodia has almost tripled from 1984 (0.978 × 10 6 ha) to 2014 (2.859 × 10 6 ha), and it continued to increase by 2.57% per year until 2020 [35]. Although the grain yield level in average has increased from 1.29 t ha −1 (1984) to 3.39 t ha −1 (2020), the yield of rainfed lowland rice in the wet season remains lower than 3.0 t ha −1 [19]. The variation in the amount of precipitation between provinces in 2019 appears larger compared with the average for five years due to the large value in Siem Reap in 2019 as opposed to the average in August and September and the smaller value in all provinces from October to December in 2019.

Yield-Related Traits in 2019
The grain yield and its components have large variation in the rainfed rice varieties grown at different environments. The grain yield per unit land area and its yield components are initially listed in Table 3 for comparison. The yield components and yield of the ten photoperiod-sensitive varieties across the eight environments were investigated ( Table 3). The ANOVA results detected a highly significant difference for the main effects of genotype and environment on yield components; moreover, the genotype-by-environment interaction effect was significant for all these components. As shown in Table 3, the heaviest mean grain yield of 449.8 g m −2 was recorded in Phka Rumduol (G1), followed by Phka Mealdei (G4), yielding 430.4 g m −2 , while the lowest mean grain yield of 309.1 g m −2 was observed in CAR11. The highest mean panicle number per hill and percentage of filled grains were observed in Phka Rumduol (G1), followed by Phka Mealdei (G4). The mean spikelet number per panicle was largest in Raing Chey (G6), and the mean 1000-grain weight was heaviest in Phka Mealdei (G4). Environment E2 (Battambang, BTB2) showed the highest mean for all yield components and grain yields, followed by E1 (BTB1). The lowest mean grain yield level was recorded in E7 and E8 in Siem Reap ( Table 3). All of the yield components were high in Battambang and low in Siem Reap.
As the next analysis, the principal component test was conducted, employing various agronomic and environmental factors, such as the maximum stem number per hill; plant length at harvesting; shoot dry matter weight at heading and harvesting; panicle length; panicle number per hill; spikelet number per panicle; percentage of filled grains; 1000-grain weight; harvest index (HI); clay content; soil pH; total C, N, and P contents; cation exchange capacity (CEC); exchangeable K content; average daily mean air temperature; and accumu-Agriculture 2023, 13, 1602 8 of 18 lated precipitation during each variety's growth period to summarize many explanatory variables into fewer indicators and identify the contribution extent of each parameter or factor to the difference in the grain yield. Three interaction principal components were extracted; the first component (PC1) and the second component (PC2) accounted for 58.1% and 15.1%, respectively, of the total variation of grain yield in the METs ( Figure 3A). The loading value for PC1 was large in the clay content, exchangeable K content, CEC, and total P content ( Figure 3B). Then, the loading value of precipitation for PC1 became a large value on the negative side ( Figure 3B), meaning that high precipitation leads to a smaller grain yield. However, this was due to originally high precipitation from August to November in Siem Reap, which has comparatively poor natural fertility compared with the other provinces. Temperature and precipitation are known to be the key determinant factors affecting rice production in the context of climate change [36]. According to Pheakdey et al. (2017) [36], the three climate variables (TCV), that is, maximum temperature, minimum temperature, and rainfall by the Ordinary Least Squares analysis, explained approximately 63% of rice yields in the wet season, which was large compared to the results in the dry season, at 56% for the period of 1993-2012 in Cambodia. The amount of rainfall did not have a direct influence on the grain yield in the cases of these four provinces in the current study. The remaining percentage could be explained by non-climate factors such as soil fertility and other factors. In the current study, the rainfall condition was considered to have a relationship with natural soil fertility partly at the selected experimental sites as described above. The total precipitation in Siem Reap is 1072 mm in the ordinal year, the value of which is larger by approximately 30, 20, and 10%, respectively, than that of BTB, PUR, and KPT ( Figure 2). As the results of the principal component test indicate, PC1 consisted of parameters relating to the environment and PC2 consisted agronomic parameters ( Figure 3A,B). condition was considered to have a relationship with natural soil fertility partly at the selected experimental sites as described above. The total precipitation in Siem Reap is 1072 mm in the ordinal year, the value of which is larger by approximately 30, 20, and 10%, respectively, than that of BTB, PUR, and KPT ( Figure 2). As the results of the principal component test indicate, PC1 consisted of parameters relating to the environment and PC2 consisted agronomic parameters ( Figure 3A,B).  In Figure 3C, the relationship between the coefficient of variance (CV) in grain yield and grain yield is shown. The CV in grain yield among the genotypes in each environment was largest by 15.1% in E2, followed by E1, and was smallest by 7.3% in E8. Additionally, the CV values showed a significantly positive correlation with the grain yield (r = 0.7171, p = 0.046). From this relationship, it can be seen that the variation in grain yield between genotypes is enlarged under fertile conditions rather than in conditions that have comparatively poor natural soil fertility. The CV in grain yield across different environments in each variety was highest by 24.0% in G4 and smallest by 15.5% in G7, values that did not show specific relationships with grain yield. The variation in grain yield depending on genotypes was apparently small in comparison with that depending on different environments. These results suggest that the eight trial environments in the current study were diverse and constituted the major cause of the heterogeneity of grain yield performances of those genotypes across the eight sites tested. This finding agrees with those of several authors who reported that environment is the main cause of the variation in yield performance [37][38][39]. Gauch and Zobel (1997) [40] stated that when the genotype-by- explained approximately 63% of rice yields in the wet season, which was large compared to the results in the dry season, at 56% for the period of 1993-2012 in Cambodia. The amount of rainfall did not have a direct influence on the grain yield in the cases of these four provinces in the current study. The remaining percentage could be explained by nonclimate factors such as soil fertility and other factors. In the current study, the rainfall condition was considered to have a relationship with natural soil fertility partly at the selected experimental sites as described above. The total precipitation in Siem Reap is 1072 mm in the ordinal year, the value of which is larger by approximately 30, 20, and 10%, respectively, than that of BTB, PUR, and KPT ( Figure 2). As the results of the principal component test indicate, PC1 consisted of parameters relating to the environment and PC2 consisted agronomic parameters ( Figure 3A,B). In Figure 3C, the relationship between the coefficient of variance (CV) in grain yield and grain yield is shown. The CV in grain yield among the genotypes in each environment was largest by 15.1% in E2, followed by E1, and was smallest by 7.3% in E8. Additionally, the CV values showed a significantly positive correlation with the grain yield (r = 0.7171, p = 0.046). From this relationship, it can be seen that the variation in grain yield between genotypes is enlarged under fertile conditions rather than in conditions that have comparatively poor natural soil fertility. The CV in grain yield across different environments in each variety was highest by 24.0% in G4 and smallest by 15.5% in G7, values that did not show specific relationships with grain yield. The variation in grain yield depending on genotypes was apparently small in comparison with that depending on different environments. These results suggest that the eight trial environments in the current study were diverse and constituted the major cause of the heterogeneity of grain yield performances of those genotypes across the eight sites tested. This finding agrees with those of several authors who reported that environment is the main cause of the variation in yield performance [37][38][39]. Gauch and Zobel (1997) [40] stated that when the genotype-by-

450, 450 >
to the results in the dry season, at 56% for the period of 1993-2012 in Cambodia. The amount of rainfall did not have a direct influence on the grain yield in the cases of these four provinces in the current study. The remaining percentage could be explained by nonclimate factors such as soil fertility and other factors. In the current study, the rainfall condition was considered to have a relationship with natural soil fertility partly at the selected experimental sites as described above. The total precipitation in Siem Reap is 1072 mm in the ordinal year, the value of which is larger by approximately 30, 20, and 10%, respectively, than that of BTB, PUR, and KPT ( Figure 2). As the results of the principal component test indicate, PC1 consisted of parameters relating to the environment and PC2 consisted agronomic parameters ( Figure 3A,B). In Figure 3C, the relationship between the coefficient of variance (CV) in grain yield and grain yield is shown. The CV in grain yield among the genotypes in each environment was largest by 15.1% in E2, followed by E1, and was smallest by 7.3% in E8. Additionally, the CV values showed a significantly positive correlation with the grain yield (r = 0.7171, p = 0.046). From this relationship, it can be seen that the variation in grain yield between genotypes is enlarged under fertile conditions rather than in conditions that have comparatively poor natural soil fertility. The CV in grain yield across different environments in each variety was highest by 24.0% in G4 and smallest by 15.5% in G7, values that did not show specific relationships with grain yield. The variation in grain yield depending on genotypes was apparently small in comparison with that depending on different environments. These results suggest that the eight trial environments in the current study were diverse and constituted the major cause of the heterogeneity of grain yield performances of those genotypes across the eight sites tested. This finding agrees with those of several authors who reported that environment is the main cause of the variation in yield performance [37][38][39]. Gauch and Zobel (1997) [40] stated that when the genotype-by-400, 400 > 350, 350 > 300, and 300 > to the results in the dry season, at 56% for the period of 1993-2012 in Cambodia. The amount of rainfall did not have a direct influence on the grain yield in the cases of these four provinces in the current study. The remaining percentage could be explained by nonclimate factors such as soil fertility and other factors. In the current study, the rainfall condition was considered to have a relationship with natural soil fertility partly at the selected experimental sites as described above. The total precipitation in Siem Reap is 1072 mm in the ordinal year, the value of which is larger by approximately 30, 20, and 10%, respectively, than that of BTB, PUR, and KPT ( Figure 2). As the results of the principal component test indicate, PC1 consisted of parameters relating to the environment and PC2 consisted agronomic parameters ( Figure 3A,B).  In Figure 3C, the relationship between the coefficient of variance (CV) in grain yield and grain yield is shown. The CV in grain yield among the genotypes in each environment was largest by 15.1% in E2, followed by E1, and was smallest by 7.3% in E8. Additionally, the CV values showed a significantly positive correlation with the grain yield (r = 0.7171, p = 0.046). From this relationship, it can be seen that the variation in grain yield between genotypes is enlarged under fertile conditions rather than in conditions that have comparatively poor natural soil fertility. The CV in grain yield across different environments in each variety was highest by 24.0% in G4 and smallest by 15.5% in G7, values that did not show specific relationships with grain yield. The variation in grain yield depending on genotypes was apparently small in comparison with that depending on different environments. These results suggest that the eight trial environments in the current study were diverse and constituted the major cause of the heterogeneity of grain yield performances of those genotypes across the eight sites tested. This finding agrees with those of several authors who reported that environment is the main cause of the variation in yield performance [37][38][39]. Gauch and Zobel (1997) [40] stated that when the genotype-by-250. SN: maximum stem number per hill; PL: plant length at harvesting; DWh: shoot dry matter weight at heading; DWf: shoot dry matter weight at harvesting; PnL: panicle length; PN: panicle number per hill; SpN: spikelet number per panicle; FG: percentage of filled grains; GW: 1000-grain weight; HI: harvest index; Cl: clay content; pH: pH(HCl); C: total C; N: total N; P: total P; CE: CEC; K: exchangeable K; Temp: average daily mean air temperature; Rain: accumulated precipitation during growth period of each variety.
In Figure 3C, the relationship between the coefficient of variance (CV) in grain yield and grain yield is shown. The CV in grain yield among the genotypes in each environment was largest by 15.1% in E2, followed by E1, and was smallest by 7.3% in E8. Additionally, the CV values showed a significantly positive correlation with the grain yield (r = 0.7171, p = 0.046). From this relationship, it can be seen that the variation in grain yield between genotypes is enlarged under fertile conditions rather than in conditions that have comparatively poor natural soil fertility. The CV in grain yield across different environments in each variety was highest by 24.0% in G4 and smallest by 15.5% in G7, values that did not show specific relationships with grain yield. The variation in grain yield depending on genotypes was apparently small in comparison with that depending on different environments. These results suggest that the eight trial environments in the current study were diverse and constituted the major cause of the heterogeneity of grain yield performances of those genotypes across the eight sites tested. This finding agrees with those of several authors who reported that environment is the main cause of the variation in yield performance [37][38][39]. Gauch and Zobel (1997) [40] stated that when the genotype-by-environment is large, genotypes should be selected for local/specific adaptation. However, the effect of genotype-by-environment interaction (GEI) was relatively small as compared with that of environment (E) and genotype (G) in the current study (Table 3). Table 4 shows the grain yield of each variety in each environment for comparison across the varieties in each environment or across the environment in each variety. Phka Rumduol (G1) in BTB, PUR, SRP, and SRP showed grain yields (g m −2 ) of approximately 546 to 566, 512 to 525, 399 to 419, and 301 to 331, respectively. The grain yield was apparently high in BTB, at the same level or above in PUR, at the same or slightly lower level in KPT, and distinctively higher in SRP as compared with that in rainy rainfed fields in 2016 and 2017, as reported by Kong et al. [21]. All of the varieties used were released for rainfed conditions; Phka Mealdei (G4), CAR3 (G7), and CAR4 (G8) seem to be drought tolerant, and Phka Rumduol (G1), Phka Romeat (G2), Phka Rumdeng (G3), and CAR6 (G9) seem to be submergence tolerant (Kong, K., Director of Rice Department, General Directorate of Agriculture, MAFF of Cambodia, personal communication). As described above, the precipitation was comparatively small at the later stage of growth in 2019, while Phka Rumduol (G1), which might be recognized as submergence tolerant, showed preferable growth. Therefore, it was considered that the grain yield of all varieties in the current study in 2019 would be at ordinary levels or above their averages. Values followed by different small letters are significantly different among varieties (genotypes) in each environment and those followed by different capital letters are significantly different among varieties (genotypes) in each environment at a 0.05 probability level by a Tukey's honest significant difference test.
The influence of environmental differences such as natural soil fertility on grain yield was similar in all varieties used. The value tended to be large in BTB and small in SRP regardless of the variety (Table 4).
Moreover, the Pearson product-moment coefficient of the grain yield is shown in Table 5, which was performed between the locations and provinces in order to clarify the difference in the grain yield between the varieties in different environments. The correlation coefficient of the relationship in the order of grain yield among ten varieties at a single environment and that in a different environment was significantly positive in any combination. This means that high-yielding varieties in a fertile environment show comparatively high yields in poor-fertility environments or intermediate environments.
In contrast, low-yielding varieties show comparatively low yields everywhere. The ideal genotypes must be both high-yielding and stable. In the current study, Phka Rumduol showed the highest grain yield in all environments (E1 to E8) and its variation depending on different environments was 22.5%, which was in the intermediate level among the ten varieties, which averaged 20.8% (Table 4, Figure 3C).

Yield-Determining Factors in 2019
The relationships between the yield components and grain yield used to investigate the yield-determining factors are shown in Figure 4. The panicle number per hill, percentage of filled grains, and 1000-grain weight showed significant positive correlations with the grain yield in all environments and in the case of all of the data pooled. On the contrary, the spikelet number per panicle showed a positive correlation with the grain yield only at E5 (KPT1).   [19] reported with regard to the growth of a popular high-quality rice cultivar "Phka Rumduol" that the effect of location (L: four provinces where are same with the current study) X management (M: kind and amount of fertilizer) on the spikelet number per panicle was not clear during the wet season of 2017. Thus, as the next step, a multiple regression analysis was conducted with the yield components as applied variables for grain yield according to the stepwise method. The percentage of filled grains was extracted as the most important variable for explaining the difference in the grain yield at E1 (BTB1) and E1 (BTB2), where there were comparatively fertile soil conditions, and at E7 (SRP1) and E8 (SRP2), where there were comparatively poor natural soil fertility conditions ( Table 6). In contrast, the panicle number per hill was extracted as the most important variable at E3 (PUR1), E4 (PUR2), E5 (KPT1), and E6 (KPT2), where there are intermediate conditions of natural soil fertility. However, the limiting factors of grain yield in each environment were not clear. Study of these ten popular photoperiod-sensitive varieties should be repeated for another year in different environments in order to confirm their real yield and stability performances. To this end, we conducted an experiment with the same concept in the following year.

Yield-Determining Factors in 2019
The relationships between the yield components and grain yield used to investigate the yield-determining factors are shown in Figure 4. The panicle number per hill, percentage of filled grains, and 1000-grain weight showed significant positive correlations with the grain yield in all environments and in the case of all of the data pooled. On the contrary, the spikelet number per panicle showed a positive correlation with the grain yield only at E5 (KPT1).   [19] reported with regard to the growth of a popular high-quality rice cultivar "Phka Rumduol" that the effect of location (L: four provinces where are same with the current study) X management (M: kind and amount of fertilizer) on the spikelet number per panicle was not clear during the wet season of 2017. Thus, as the next step, a multiple regression analysis was conducted with the yield components as applied variables for grain yield according to the stepwise method. The percentage of filled grains was

Yield and Its Components with Relating Parameters in 2020
In the 2020 study, four photoperiod-sensitive varieties were examined: Phka Rumduol (G1 in 2019), which was the highest yielding variety in the previous experiment in 2019; Phka Mealdei (previous G4), which was high-yielding on a level similar to that of Phka Rumduol; CAR4 (previous G8), which was a high-yield variety after G1 and G4; and Raing Chey (previous G6), which was comparatively low yielding and is a popular variety cultivated in over 70% of provinces in Cambodia [41]. were used. These four varieties were cultivated at one site each in the four provinces and the cultivation management was the same as in 2019. Thus, similar analytical procedures as were used for 2019, such as ANOVA for investigating the influence of genotype or environment on the grain yield, its difference between the varieties or environments, the relationships between the grain yield per unit land area, and the yield components between the important yield components and agronomic parameters or plant nutrient status were employed for the results in 2020. Climate conditions such as the change in daily mean air temperature and monthly rainfall in in July to December 2020 were very similar to the average for the most recent five years, and it can be said that the rice cultivation season in 2020 was ordinal. The results of twoway factorial ANOVA for grain yield are shown in Table 7. The effects of variety (genotype: G), province (environment: E), and G × E on the grain yield were significant; the results are the same as those from the previous year 2019. The order of varieties in grain yield in each environment and that of environment in each variety are shown in Table 8. Their tendencies were similar to those determined in the previous year's experiment ( Table 4). The effect of environment on the variation in grain yield was larger than that of genotype, and this tendency was similar to that in 2019 ( Table 4). As shown in Figure 2, the monthly rainfall in 2019 tended to be smaller compared to the average for the most recent 5 years except in BTB in July, KPT in September, and SRP in August and September. However, the difference in grain yield between the varieties in each province and between the provinces in each variety was confirmed to be similar in 2019 and 2020, which means the effect of the soil environment on the grain yield can be considered large compared to that of the precipitation.  Table 8. Grain yield (g m −2 ) of four varieties in four different environments (2020). Values followed by different small letters are significantly different among varieties (genotypes) in each environment, and those followed by different capital letters are significantly different among varieties (genotypes) in each environment at a 0.05 probability level by a Tukey's honest significant difference test.

Variety
The relationships between the yield components and grain yield are shown in Figure 5. The panicle number per hill had a positive correlation with the grain yield in four environments as well as in the case when all of the data were pooled. The spikelet number per panicle did not show any specific relationship with grain yield in four environments. The percentage of filled grains showed a positive relationship with yield at BTB (fertile condition) and SRP (poor fertility) and in all the data pooled, and the 1000-grain weight showed a positive relationship with the yield at PUR and KPT (comparatively intermediate fertility) and in all of the pooled data. The trend in the relationships between the yield components and grain yield was similar with that in the previous year, 2019. According to Banayo et al. (2018) [42], in the absence of any specific growth constraints, sink size (the number of spikelets m −2 ) is a major yield constraint for rainfed lowland rice. In the current study, the variation in spikelet number per panicle was comparatively small in sandy soil condition (SRP) rather than in the other soil types (BTB, PUR, KPT) ( Figure 5). However, the difference in spikelet number per panicle between the varieties did not reflect in the difference in grain yield in any soil types. Therefore, it was considered that sink size was determined primarily by the panicle number per hill in rainfed condition in Cambodia. Thus, we investigated the parameters that have a unique relationship with the panicle number per hill and the percentage of filled grains that were determined to be the important parameters depending on the genotype in different environments obtained by multiple regression analysis in the previous experiment (Table 6). Then, we found that the days from transplanting to the first flowering showed a negative relationship with grain yield in the four environments ( Figure 6A). Additionally, the stem number per hill 42 days after transplanting (DAT) showed a negative relationship with panicle number per hill at BTB; from this, it can be understood that a smaller stem number around the maximum tillering stage in a fertile soil condition is better for securing a larger number of panicles ( Figure 6B). On other hand, the stem numbers at 42DAT showed a positive correlation with panicle number per hill, which means that a larger number of stems around the maximum tillering stage should be secured for a larger number of panicles ( Figure 6B). At the same time, in the relationship between the stem number per hill at 42 DAT and the percentage of productive culms, an optimal number of stems, specifically around ten stems per hill, was found to produce a preferable higher percentage of productive culms ( Figure  6C). In addition, the percentage of productive culms showed a positive relationship with the panicle number per hill ( Figure 6D). Thus, we investigated the parameters that have a unique relationship with the panicle number per hill and the percentage of filled grains that were determined to be the important parameters depending on the genotype in different environments obtained by multiple regression analysis in the previous experiment (Table 6). Then, we found that the days from transplanting to the first flowering showed a negative relationship with grain yield in the four environments ( Figure 6A). Additionally, the stem number per hill 42 days after transplanting (DAT) showed a negative relationship with panicle number per hill at BTB; from this, it can be understood that a smaller stem number around the maximum tillering stage in a fertile soil condition is better for securing a larger number of panicles ( Figure 6B). On other hand, the stem numbers at 42DAT showed a positive correlation with panicle number per hill, which means that a larger number of stems around the maximum tillering stage should be secured for a larger number of panicles ( Figure 6B). At the same time, in the relationship between the stem number per hill at 42 DAT and the percentage of productive culms, an optimal number of stems, specifically around ten stems per hill, was found to produce a preferable higher percentage of productive culms ( Figure 6C). In addition, the percentage of productive culms showed a positive relationship with the panicle number per hill ( Figure 6D). Furthermore, the percentage of productive culms had a positive relationship with the leaf N concentration of the flag leaf at the late stage (after heading) ( Figure 7A). Ehara et al. (1997) [43] reported that the leaf N concentration indicates the extent of nitrogen absorption capacity of each variety, especially NH 4+ absorption on a root weight basis. The nitrogen concentration in rice seedlings correlates with tiller emergence [44], and the soluble nitrogen in plants is reflected by the amount of total nitrogen.
With regard to tillering [45], Ehara et al. (1992a) [46] stated that securing sufficiently high nitrogen content per unit leaf area of rice seedlings is important for attaining vigorous tillers that will become productive culms. Our current findings in the relationship between the leaf N concentration and the percentage of productive culms are in good agreement with the findings of previous reports [44][45][46]. On the other hand, the leaf N concentration tended to be in a positive relationship with the leaf width of the flag leaf blade at flowering; this correlation coefficient was significant in BTB and PUR as well as when all of the data were pooled ( Figure 7B). At same time, leaf width showed a positive correlation with the percentage of filled grains in BTB ( Figure 7C). In fertile soil conditions, the larger leaf width, leading to a larger leaf area, contributes to an increased percentage of filled grains through maintenance of the photosynthetic rate until a later stage of growth. Moreover, leaf width did not have a positive relationship with leaf length of the flag leaf blade at flowering in any soil environment in the current study ( Figure 7D). Of course, a longer leaf blade will lead to a larger leaf area in general, although in the case of BTB the larger leaf area of the flag leaf supported by the wider leaf blade, not by the longer leaf blade, was considered to be important. Ehara et al. (1992b) [47] reported that a larger leaf area supporting the elongation of the leaf blade with its erect form is preferable for attaining a higher net assimilation rate in rice plants. In the current study, we found a tendency that differed from that found in a former study [47], which showed that due to Furthermore, the percentage of productive culms had a positive relationship with the leaf N concentration of the flag leaf at the late stage (after heading) ( Figure 7A). Ehara et al. (1997) [43] reported that the leaf N concentration indicates the extent of nitrogen absorption capacity of each variety, especially NH 4+ absorption on a root weight basis. The nitrogen concentration in rice seedlings correlates with tiller emergence [44], and the soluble nitrogen in plants is reflected by the amount of total nitrogen.
With regard to tillering [45], Ehara et al. (1992a) [46] stated that securing sufficiently high nitrogen content per unit leaf area of rice seedlings is important for attaining vigorous tillers that will become productive culms. Our current findings in the relationship between the leaf N concentration and the percentage of productive culms are in good agreement with the findings of previous reports [44][45][46]. On the other hand, the leaf N concentration tended to be in a positive relationship with the leaf width of the flag leaf blade at flowering; this correlation coefficient was significant in BTB and PUR as well as when all of the data were pooled (Figure 7B). At same time, leaf width showed a positive correlation with the percentage of filled grains in BTB ( Figure 7C). In fertile soil conditions, the larger leaf width, leading to a larger leaf area, contributes to an increased percentage of filled grains through maintenance of the photosynthetic rate until a later stage of growth. Moreover, leaf width did not have a positive relationship with leaf length of the flag leaf blade at flowering in any soil environment in the current study ( Figure 7D). Of course, a longer leaf blade will lead to a larger leaf area in general, although in the case of BTB the larger leaf area of the flag leaf supported by the wider leaf blade, not by the longer leaf blade, was considered to be important.   [47] reported that a larger leaf area supporting the elongation of the leaf blade with its erect form is preferable for attaining a higher net assimilation rate in rice plants. In the current study, we found a tendency that differed from that found in a former study [47], which showed that due to specific characteristics of Cambodian rice different varieties tend to have extremely long leaves.  Therefore, in the case of Cambodian rice varieties a larger leaf area supported by wider leaf blades will be preferable. Furthermore, the leaf N concentration after heading showed an especially positive relationship with the percentage of filled grains ( Figure 8A). Then, the leaf N concentration at harvesting had a positive relationship with the percentage of filled grains in all of the environments ( Figure 8B). From these results, we are certain that higher leaf nitrogen nutrition is valuable for securing a higher percentage of filled grains. However, as Khema et al. (2022) [41] suggested, a higher nitrogen concentration in rice plants leads to a higher protein concentration in the grains until the late stage of growth, which is a negative factor when discussing the quality of rice grain, that is to say, the taste of the cooked rice. It may be necessary to consider a kind of trade-off between rice grain quantity and quality for rice cropping in rainfed conditions.  Therefore, in the case of Cambodian rice varieties a larger leaf area supported by wider leaf blades will be preferable. Furthermore, the leaf N concentration after heading showed an especially positive relationship with the percentage of filled grains ( Figure 8A). Then, the leaf N concentration at harvesting had a positive relationship with the percentage of filled grains in all of the environments ( Figure 8B). From these results, we are certain that higher leaf nitrogen nutrition is valuable for securing a higher percentage of filled grains. However, as Khema et al. (2022) [41] suggested, a higher nitrogen concentration in rice plants leads to a higher protein concentration in the grains until the late stage of growth, which is a negative factor when discussing the quality of rice grain, that is to say, the taste of the cooked rice. It may be necessary to consider a kind of trade-off between rice grain quantity and quality for rice cropping in rainfed conditions.  Therefore, in the case of Cambodian rice varieties a larger leaf area supported by wider leaf blades will be preferable. Furthermore, the leaf N concentration after heading showed an especially positive relationship with the percentage of filled grains ( Figure 8A). Then, the leaf N concentration at harvesting had a positive relationship with the percentage of filled grains in all of the environments ( Figure 8B). From these results, we are certain that higher leaf nitrogen nutrition is valuable for securing a higher percentage of filled grains. However, as Khema et al. (2022) [41] suggested, a higher nitrogen concentration in rice plants leads to a higher protein concentration in the grains until the late stage of growth, which is a negative factor when discussing the quality of rice grain, that is to say, the taste of the cooked rice. It may be necessary to consider a kind of trade-off between rice grain quantity and quality for rice cropping in rainfed conditions.

Conclusions
Our principal component test results showed that the heterogeneity of grain yield was mainly attributable to the effect of the environment, followed by the effect of genotype. In fertile conditions, a higher percentage of filled grains was supported by the higher leaf N until late stage, with a wider flag leaf; in moderate fertility conditions, more panicles were supported by a higher percentage of productive culms with a higher leaf N until the late stage, while in conditions of poor fertility a higher percentage of filled grains were supported by higher leaf N until the late stage. These results are important for research into higher grain yields. The Phka Rumduol variety has these preferable traits, and can be expected to produce higher yields in fertile to poor natural soil fertility conditions with moderate variations. This variety may be considered more desirable and ideal due to its stability and high grain yield. Other varieties, namely, Phka Mealdei, Phka Rumdeng, and CAR4, were identified as above-average yielders. Therefore, these varieties are potentially recommended for cultivation in the rainfed lowland rice ecosystems of Cambodia. CAR4 showed moderate variation at the same level as Phka Rumduol. Thus, from the point of stability, CAR4 can be expected to excel.