1. Introduction
With the escalating water scarcity in agriculture, the need to enhance water-use efficiency in rice cultivation (Oryza sativa L.) has become paramount. Agriculture accounts for 70% of the world's freshwater consumption, a figure that has been exacerbated by factors such as population growth, rapid increases in urbanization, and industrial expansion [
1,
2]. The dwindling availability of water resources due to these challenges accentuates the need for innovative solutions [
3]. Addressing the pressing challenge of increasing rice production to meet the demand of a growing population amidst the backdrop of diminishing water resources, a novel water-saving technique has emerged: drip irrigation with plastic-film mulch (DI) [
4,
5]. DI is a novel approach to rice cultivation wherein an opaque plastic covering is used to mulch the soil surface without standing water throughout the entire rice growth cycle. This water-efficient method is designed to optimize the utilization of water resources for crop production by inhibiting the transfer of water vapor between the soil and the atmosphere, ultimately minimizing the reduction in crop yields associated with conventional rice farming practices.
As an integral component of a plant’s organs, the root plays a pivotal role by providing anchorage for the plant, enhancing nutrient and water uptake, and synthesizing plant hormones, organic acids, and amino acids [
6,
7]. However, the ability of the root to effectively carry out its role depends on its type, structure, and modification. This significance extends to water-saving rice production systems, where root growth and structure can be instrumental to the development of the rice plant and grain yield production [
8]. Therefore, it is necessary to determine the variation in the root responses of different rice varieties under DI. Enhancements in root morpho-physiological traits are closely intertwined with heightened shoot biomass accumulation, improved efficiency in the use of nitrogen, increased water productivity, and augmented resistance to plant lodging in various crops [
9,
10,
11]. The alteration of root structure stands out as a strategy to foster the development of a deep root system in rice. This can be accomplished by stimulating deeper elongation, particularly in conditions of low soil moisture, and by encouraging increased root branching to optimize soil water extraction [
12,
13]. The ability of roots to penetrate the soil to greater depths holds paramount importance in establishing a robust root system, which holds particular significance in gramineous crops with their robust roots [
14,
15]. While differences among cultivars and the application of breeding practices can influence the formation of deep root systems, the potential to induce such systems through effective cultivation management practices remains a compelling avenue for exploration [
16].
In recent years, up-to-date knowledge has been provided by many researchers on water productivity, plant growth, physiological characteristics, panicle structure, greenhouse gas emissions (GWP), and grain yield potential of DI [
4,
5,
17,
18,
19]. However, only a few studies evaluate the varying responses of rice root systems to DI. Therefore, the objectives of this study were (1) to compare the root characteristics, water productivity, and grain yield performance of rice varieties under CF and DI, (2) to compare varietal differences in root patterns in response to DI, and (3) to examine which of the root factors contributes to the improvement of rice grain yield and water productivity in DI.
2. Materials and Methods
2.1. Information on experimental design, soil, and weather
Rice crops were cultivated at the Faculty of Environmental Horticulture research farm, Chiba University, Japan (35.78N, 139.9E), in the summer seasons of 2020 and 2021. The field experimental design of crop management has been previously reported [
19]. The soil at the experimental site was a low-humic Andosol with 39.9% sand, 29.9% silt, and 30.2% clay at 0-20 cm depth. The pH of the soil was 5.2, while the amount of organic matter and total N at 0-20 cm depth were 33.3 g kg
-1 and 3.6 g kg
-1, respectively. The weather parameters, including the mean temperature, sunshine hours, and rainfall during rice-growing periods, are shown in
Figure 1.
The field experiment was arranged in a split-plot design (main plot: irrigation system; subplot: variety) with three replicates. The two rice irrigation systems were as follows. (1) CF: continuous flooded system. In this rice production system, each plot had concrete boundaries with an area of 7.5m2. The field was flood-irrigated so that a 4-5 cm water layer above the soil surface was kept until two weeks before harvest. (2) DI: water-saving drip irrigation with plastic-film mulch. In this system, raised beds were covered with black polyethylene film. The plots were buried within the rows at a depth of 40 cm to prevent lateral percolation between neighboring plots. Irrigation was applied via the tube when volumetric soil moisture content fell below 75% of the field capacity in both years. The irrigation schedule was monitored using an FDR soil moisture meter (DIK321A, Daiki Rika Kogyo Co., Ltd). The two irrigation systems were located adjacent to each other, and both were adopted for direct seeding.
Four rice varieties were examined. (1) Norin24 (upland rice, abbreviated as Norin); (2) Mitsuhikari 2003 (F1 lowland rice, abbreviated as Mitsu); (3) Hokuriku 193 (animal feed lowland rice, abbreviated as Hoku); (4) Koshihikari (lowland rice, abbreviated as Koshi). Viable seeds were directly sown in the field on June 8 and June 1 in 2020 and 2021, respectively. The plant density was 69.6 hills m-2. Before planting, chemical fertilizers were administered at a rate of 40 kg ha-1 (N:P: K = 20:20:20) on each field, and then top-dressing fertilization was carried out by applying 10 kg ha-1 of nitrogen through foliar spraying at the panicle initiation stage in both years. Other key cultivation management tasks, such as weed control and pesticide spraying, were conducted following local recommendations.
2.2. Sampling and Measurements
Aboveground and root biomass were estimated 4 or 5 times at the panicle-initiation, heading, grain-filling, and maturity stages. During each sampling instance, we gathered samples from ten rice hills, which collectively occupied an area of 0.5 m
2 within each plot. Roots were sampled from each replicate at the panicle-initiation stage (DAS50) and heading stage (DAS85), using the dimensions of the monolith soil sampler (150 mm diameter, 30 cm depth) with the aid of five soil-sampling cores. Soil samples were divided into 0-15 cm and 15-30 cm segments and washed carefully on a 0.5 mm stainless steel sieve under running water to separate the roots. The data ignored the roots below 30 cm for both CF and DI. The portion of each root sample was used for the measurement of root length. Roots were arranged on a glass tray, scanned using a scanner (Epson GT-X970, Seiko Epson Corp., Nagano, Japan), and analyzed using ImageJ as described by Tajima et al., 2013 [
20]. Finally, the remaining roots were oven-dried at 80℃ for 72 h and weighed.
To determine the differences in root bleeding rate among varieties grown under CF and DI, 5 hills were marked on each plot at the panicle-initiation stage on July 20 and July 24 in 2020 and 2021. The rice stems were trimmed to a height of 10 centimeters above the soil level at 6:00 PM local time. We then placed 5 grams of cotton wool on the upper portion of each cut stem and enclosed it with a polyethylene bag. After 12 hours, we retrieved the cotton wool and determined the root bleeding rate by measuring the increase in weight of the cotton wool.
A portable photosynthetic apparatus system of Li-6400 (LI-COR, Lincoln, NE, USA) was used to measure the CO2 assimilation rate (An), stomatal conductance (gs), and transpiration rate in the flag leaf of each rice variety under the two water regimes at the early grain-filling stage in 2021. The measurements were carried out when photosynthesis was most active (09:00 to 12:00 local time). Throughout the measurement period, the photosynthetic photon flux density level, average leaf temperature, and average vapor pressure deficit were maintained at 1000 μmolm⁻²s⁻¹, 25°C, and 1.5 kPa, respectively.
Rice plants from a 1 m2 area in each plot were collected to determine grain yield at physiological maturity each year. We recorded the count of panicles for each treatment. To determine the rough grain yield, the panicles were threshed and weighed for each rice variety. The grain weights were recorded for each variety under CF and DI after the grains were dried to a moisture content of 14%. The percentage of grains that had ripened was calculated by dividing the count of filled grains by the total number of grains. Harvest index was defined as grain yield per biomass at maturity. Water productivity was calculated as the ratio of the total quantity of water supplied (irrigation and rainfall) to each variety per grain yield.
2.3. Data analysis
Analysis of variance (ANOVA) was performed using SPSS 20.0 (IBM, New York, USA). The statistical model used included sources of variation due to replication, year (Y), variety (V), irrigation regime (IR), and the interaction of Y x V, Y x IR, V x IR, and Y x V x IR. Data from each sampling date were analyzed separately. The means were separated by Duncan’s multiple range test when the effects were significant in each year.
4. Discussion
Prior to this study, little information was available describing varietal differences in rice root
morphology and physiology under CF and DI systems. Our results indicated that in DI, root growth at the soil surface was comparable with that under CF during the vegetative phase (
Figure 4). Araki et al., 2005 [
21] and Luo et al., 2023 [
22] revealed that plants first draw water from the surface layer, and subsequently, the area of water extraction gradually shifts downward through the soil profile. When the water content in the shallow layer (0-25 cm depth) was abundant, over 90% of the water uptake occurred in this shallower zone [
23,
24]. Maintaining paddy soils in a wet but not continuously saturated state during the vegetative phase has the effect of preserving predominantly aerobic soil conditions [
25]. These particular soil management practices have significantly facilitated better root growth at the surface layer with the DI approach. The results of our study suggest that current levels of root growth at the surface in CF can be achieved under DI. Our results are consistent with those of Bouman et al., 2006 [
26], who showed no significant differences in root characteristics of rice grown under CF compared with an alternate wetting and drying (AWD) irrigation system. We observed that there were large gen
otypic differences in root distribution patterns
. Regarding the deep-root ratio, the lowland varieties produced comparable biomass to subsurface root growth in DI; however, upland varieties significantly allocated more root growth at the lower soil horizons in DI compared to CF (
Figure 5). This is in line with the conclusion of Kato et al., 2010 [
27], who used aerobic cultivation technology. It is known that the soil moisture content, soil nutrient status, and soil physical properties can profoundly influence plant root growth [
28].
The deep root system in upland rice varieties enables them to harness water from the deeper soil layers, a resource that remains out of reach for cultivars with shallow root systems. This capability has proven effective in preserving yields, especially under reduced soil moisture conditions [
29,
30,
31,
32]. For a more precise assessment of rice's deep-root characteristics under water management, soil samples were taken and measured up to a depth of 60 cm by Kato et al., 2007 [
33]. However, due to the presence of concrete boundaries under CF, our study was limited to measurements up to 30 cm below the surface. This discrepancy revealed that the rice varieties under DI had greater deep root growth potential than under CF. The root bleeding rates conducted during the panicle-initiation stage over two consecutive years identified substantial genotypic differences among the rice varieties (
Figure 6). Interestingly, lowland varieties, such as Mitsu and Koshi, exhibited significantly higher root activity under CF than under DI. This observation aligns with the findings reported by Matsuo et al., 2009 [
34], where a decrease in root activity for lowland variety Koshi under water-saving cultivation systems was similarly noted. Furthermore, significant differences in root activity among upland (Norin) and animal feed varieties (Hoku) under CF and DI were not observed (
Figure 6). This could be attributed to the influence of root- and stem-xylem resistance mechanisms. The presence and impact of root- and stem-xylem resistance mechanisms represent essential factors that warrant consideration when assessing disparities in root activity between two irrigation systems [
35,
36]. These resistance mechanisms can substantially influence the transport of water and sap within the plant, potentially masking or reducing variations in root bleeding rates [
37].
The results of this study showed no significant difference among the photosynthetic-associated parameters under CF and DI. However, it was observed that substantial genotypic variation existed among the tested rice varieties (
Table 1). The transpiration rate of Norin showed a significant increase of 11% under DI due to greater deep-root development (
Figure 7a). Deep roots are instrumental in enhancing a plant's resilience in drought conditions, ensuring that transpiration can persist even when surface soil dries out. This strategy promotes drought tolerance and improves transpiration efficiency [
38,
39]. The regulation of deep-root functioning and distribution, adjustments in canopy size, and the synthesis of the hormone abscisic acid (ABA) in both roots and leaves collectively play a central role in controlling leaf water potential (LWP) and enhancing leaf transpiration efficiency [
40,
41,
42]. Under severe drought-stress conditions, abscisic acid (ABA) typically assumes a more critical and pronounced role in the regulation of plant responses to water stress. However, the significance of ABA in plant responses may not be as pronounced in less severe drought-stress conditions. Shi et al., 2015 [
43] reported that both upland cultivars (Gaoshan1) and lowland cultivars (Nipponbare) demonstrated the ability to open their stomata for transpiration even under slight-to-moderate drought conditions, which allowed the effective loading of root-generated ABA into the xylem through apoplastic by-pass flow, subsequently reducing ABA exudation. Furthermore, root structural distributions, including deep-root ratio and root branching, are likely to be promoted under aerobic conditions [
44]. A study by Kato et al., 2013 [
45] in northeastern Thailand utilized an irrigation system for aerobic cultivation and obtained similar results, supporting the notion that aerobic conditions can enhance root structure distributions. Deep-rooted traits often exhibit better transpiration efficiency, which means they can achieve higher photosynthesis per unit of water transpired [
46]. Our results indicated that the deep roots of the upland variety (Norin) could potentially contribute to a high transpiration rate, thus maintaining a high assimilation rate under DI in comparison to CF (
Figure 7a, Fig 8). Additionally, while deep roots are one of the factors influencing transpiration, other environmental factors (i.e., temperature, humidity, and light intensity) also contribute to the regulation of transpiration rates [
47]. In future research, it will be essential to continue investigating the complex interplay among these factors to gain a more comprehensive understanding of how plants manage water resources and optimize transpiration rates to adapt to various environmental conditions, especially in aerobic and mild water-stress conditions. On the other hand, the lowland variety (Koshi) significantly decreased gas-exchange performance (An, g
s, and E), likely due to the lower plasticity index for lower root bleeding rate (
Table 1,
Figure 6). Morita et al., 2000; He et al., 2019; Ansari et al., 2004 and Liu et al., 2021 [
48,
49,
50,
51] reported that root bleeding is closely linked to the active absorption of water and nutrients, photosynthate transport, and redistribution efficiency in the root system. In a study by Li et al. conducted on winter maize in 2022 [
52], they similarly observed a significant effect of water management on various factors, such as root oxidation activity, cytokinin concentrations in both roots and shoots, the rate of photosynthesis in leaves, and the activities of essential enzymes responsible for converting sucrose to starch in grains. The reduction in yield observed in the lowland varieties under the DI system could be attributed to alterations in root activity, function, and morphology, which may disrupt the effective transport and redistribution of photosynthates within the plants.
In general, the primary goals of water-saving technologies are to achieve higher water-use efficiency while maintaining or even increasing crop yields [
53]. Our research findings indicate that the average grain yield under DI was equivalent to that under the CF system. Meanwhile, DI demonstrated significant water-saving potential, saving approximately 35% of the total water supplied compared to CF, resulting in water productivity that was on average 53% higher (
Table 2). Among the varieties, DI significantly boosted the grain yield of the upland variety (Norin) and the comparable yield of animal-feed lowland variety (Hoku). The increase in yield can be primarily attributed to the higher number of spikelets per square meter (m²) in these cultivars, despite a trade-off in ripening ratio and individual grain weight. However, lowland rice varieties (Mitsu and Koshi) produced significantly higher yields under CF. This higher yield can be attributed to their superior ripening ratio and grain weight under CF conditions, offsetting potential trade-offs in other components. On the other hand, a sizeable genotypic variation was observed between root traits and grain yield. The root systems of the upland variety (Norin) under DI were characterized by well-developed root growth on the surface and produced more deep roots compared to the CF sample. Therefore, it can be inferred that those characteristics contributed to a high transpiration rate, thus maintaining a high assimilation rate resulting in the high-yield performance and water productivity of the upland variety in production. Furthermore, there was a significant positive relationship between deep-root ratio, grain yield (0.745*), and water productivity (0.861
*) under DI (
Figure 7b, c). The results indicated that vertical root distribution (i.e., the deep-root ratio) during the vegetative phase influenced the yield and water productivity. The fact that animal-feed lowland rice (Hoku) maintained a comparable yield under DI compared to CF is likely due to the high root-to-shoot ratio during the early growth stage (
Figure 3). This high root-to-shoot ratio of Hoku during the early growth stage contributed to the prioritization of root development, which can increase a plant's capacity to take up essential nutrients, promoting healthy growth and development [
54]. These adaptations in the root system facilitated a greater allocation of photosynthetic products to the shoots in particular leaves, ultimately leading to a reduced root-to-shoot ratio and an enhanced yield while conserving water [
55,
56]. The low yield of lowland varieties in DI was likely due to the low root vigor and water uptake. The relationship between root bleeding rate and grain yield in DI differed between the lowland varieties (Mitsu and Koshi) and the other two. There was a significant positive relationship between root bleeding rate and grain yield under CF (
Figure 4). This highlights the importance of the root activity of lowland rice in influencing grain yield in CF conditions, where root activity plays a critical role in crop performance.
Development and identification of suitable genotypes and crop management options are underway worldwide to make agriculture more resource-use efficient and productive under drip irrigation with plastic-film mulch [
57]. As a first for future studies, exploring varietal differences under drip irrigation is an essential avenue for future research in agriculture. To further investigate and select adaptable crop varieties, developing guidelines for selecting rice varieties optimized for drip irrigation with plastic-film mulch is crucial, considering local climate, soil conditions, plant density, and fertilizer management. Secondly, the alterations in the root structures of deep root systems result from increased root branches [
14]. Further research is needed to improve the development of the rice root system in drip irrigation systems, including the study of adaptive responses, such as adventitious root emergence, lateral root branching, and water uptake. Thirdly, because of the limited extent of the research on root morphology and physiological characteristics in DI systems, additional studies are essential to elucidate the specific traits that contribute to optimal crop performance under DI conditions.
Figure 1.
Average daily air temperature and sunshine hour during the rice growing seasons of 2020(A) and 2021(B).
Figure 1.
Average daily air temperature and sunshine hour during the rice growing seasons of 2020(A) and 2021(B).
Figure 2.
Rainfall, irrigation, and volumetric soil moisture content under drip irrigation system with plastic-film mulch during the rice growing seasons of 2020(A) and 2021(B).
Figure 2.
Rainfall, irrigation, and volumetric soil moisture content under drip irrigation system with plastic-film mulch during the rice growing seasons of 2020(A) and 2021(B).
Figure 3.
Shoot dry matter (a and b), root dry matter (c and d), and root-to-shoot ratio (e and f) of four rice varieties under CF and DI in 2021. Data indicate the mean of three replications. Open symbols and solid symbols indicate CF and DI, respectively.
Figure 3.
Shoot dry matter (a and b), root dry matter (c and d), and root-to-shoot ratio (e and f) of four rice varieties under CF and DI in 2021. Data indicate the mean of three replications. Open symbols and solid symbols indicate CF and DI, respectively.
Figure 4.
Root length, and root weight at panicle-initiation stage(a)(b) and heading stage(c)(d) as affected by varieties and irrigation regimes at the soil depths of 0-15, and 15-30cm, respectively. Different capital letters beside bars show significant differences (P<0.05) between CF and DI, and different lowercase letters beside bars show significant differences (P<0.05) among varieties for each irrigation system by Duncan’s multiple range test.
Figure 4.
Root length, and root weight at panicle-initiation stage(a)(b) and heading stage(c)(d) as affected by varieties and irrigation regimes at the soil depths of 0-15, and 15-30cm, respectively. Different capital letters beside bars show significant differences (P<0.05) between CF and DI, and different lowercase letters beside bars show significant differences (P<0.05) among varieties for each irrigation system by Duncan’s multiple range test.
Figure 5.
Deep root ratio at panicle-initiation stage(a) and heading stage(b) as affected by irrigation system and rice variety. Different capital letters beside bars show significant differences(P<0.05) between CF and DI, and different lowercase letters beside bars show significant differences(P<0.05) among varieties for each irrigation system by Duncan’s multiple range test. Vertical bars represent the standard error of the mean.
Figure 5.
Deep root ratio at panicle-initiation stage(a) and heading stage(b) as affected by irrigation system and rice variety. Different capital letters beside bars show significant differences(P<0.05) between CF and DI, and different lowercase letters beside bars show significant differences(P<0.05) among varieties for each irrigation system by Duncan’s multiple range test. Vertical bars represent the standard error of the mean.
Figure 6.
Root bleeding rate at the panicle-initiation stage as affected by the irrigation system and rice variety in 2020(a) and 2021(b), respectively. Different capital letters beside bars show significant differences(P<0.05) between CF and DI, and different lowercase letters beside bars show significant differences(P<0.05) among varieties for each irrigation system by Duncan’s multiple range test. Vertical bars represent the standard error of the mean.
Figure 6.
Root bleeding rate at the panicle-initiation stage as affected by the irrigation system and rice variety in 2020(a) and 2021(b), respectively. Different capital letters beside bars show significant differences(P<0.05) between CF and DI, and different lowercase letters beside bars show significant differences(P<0.05) among varieties for each irrigation system by Duncan’s multiple range test. Vertical bars represent the standard error of the mean.
Figure 7.
Relationship between deep root ratio and transpiration rate(a), grain yield(b), and water productivity(c) in 2021. Open and solid symbols indicate CF and DI, respectively. * means significant differences at P<0.05. n=24.
Figure 7.
Relationship between deep root ratio and transpiration rate(a), grain yield(b), and water productivity(c) in 2021. Open and solid symbols indicate CF and DI, respectively. * means significant differences at P<0.05. n=24.
Figure 8.
Pearson’s correlation coefficient value of yield, water productivity (Wp), CO2 assimilation rate (An), stomatal conductance(gs), deep root ratio (DR), root length (RL), root weight (RW), and root bleeding rate (BR) under CF (a), DI (b) and CF&DI (c), respectively.
Figure 8.
Pearson’s correlation coefficient value of yield, water productivity (Wp), CO2 assimilation rate (An), stomatal conductance(gs), deep root ratio (DR), root length (RL), root weight (RW), and root bleeding rate (BR) under CF (a), DI (b) and CF&DI (c), respectively.
Table 1.
Photosynthetic-associated parameters of four rice cultivars at the early grain-filling stage under CF and DI in 2021.
Table 1.
Photosynthetic-associated parameters of four rice cultivars at the early grain-filling stage under CF and DI in 2021.
Table 2.
Water productivity, harvest index, yield, and yield component of four rice varieties under CF and DI in 2020 and 2021.
Table 2.
Water productivity, harvest index, yield, and yield component of four rice varieties under CF and DI in 2020 and 2021.