# Silicon Improves Yield Performance by Enhancement in Physiological Responses, Crop Imagery, and Leaf and Culm Sheath Morphology in New Rice Line, PadiU Putra

Authors: Zulkarami Berahim; Mohamad Husni Omar; Nurul-Idayu Zakaria; Mohd Razi Ismail; Rhushalshafira Rosle; Nor Athirah Roslin; Nik Norasma Che’Ya
Journal: BioMed Research International
 (2021)

---

## Abstract

The PadiU Putra rice line is a blast-resistant and high-yield rice line with high potential. The application of topdressing and the foliar applied method of silicon (Si) treatments could strengthen the culm to resist breakage and ultimately increase yield production. Treatments which consisted of a control, a Si topdressing, and a Si foliar applied were arranged in a randomised complete block design. At 55 days after transplanting (DAT), the foliar applied Si treatments had 59% higher dry matter partitioning to the roots. Meanwhile, at 75 DAT, both Si foliar applied and topdressing method showed increased assimilate partitioning into the culm sheath by 29% and 49%, respectively. Dark green and light yellowish colours were obtained in both Si treatments using UAV, indicating similar results to physiological responses. Remarkably, Si foliar applied treatments enhanced the diameter and width of the outer and inner layers of the diameter of vascular bundles at 75 DAT by 58, 181, and 80%, respectively. The yield production of rice increased by 53% in the Si foliar applied, compared to the control, and produced a 1.63 benefit-cost ratio.


---

## Body

## 1. Introduction

Rice,Oryza sativa L., is a major cereal crop. Sixty per cent of the global population (4.5 billion) lives in Asia, a figure which will increase to 9.3 billion in the year 2050. The demand for rice will also increase by up to 50% due to the increasing population [1]. In Malaysia, rice is the third most widely planted crop after oil palm and rubber. In 2013, 674,332 hectares (ha) of land was used for plantations with an average yield of only 4.5 t ha-1 and with a self-sufficiency level (SSL) of 72% [2]. Following that, the Malaysian government has further targeted an increase in the average paddy yield from 4.5 t ha-1 to 6.0 t ha-1 [3]. Recently, a new variety of rice, currently called PadiU Putra, was developed by [4, 5], in which resistance genes from the rice variety Pongsu Seribu 2 were pyramided into MR219. This new variety is claimed to be a blast-resistant rice variety with high yield potential and is to be released for commercial cultivation [5]. As a new variety, it is also exposed to the lodging problem due to being a taller plant. Lodging is defined as the permanent displacement of a stem from its upright position [6]. It mostly occurs just before harvest when the plant, in particular the lower portion of the stem, is unable to withstand the weight of the panicle. To date, three types of lodging have been recorded, namely, bending of stems at the base, breakage of stems at any point along the length, and root lodging [7].The lodging resistance of a rice plant could be achieved through a suitable plant height and shorter basal stems [8]. However, the study conducted by [9] showed that the application of silicon (Si) fertiliser increased the plant height. Leaf and culm morphology can be affected by the application of Si fertiliser; however, few studies have focused on the structural changes to the leaf and culm. The present study and [10–12] provide details on variations in the silica cell structure and Si deposition process by using scanning electron microscopy (SEM). SEM provides Si localisation in various parts of the rice plant such as the leaf-blade, leaf-sheath, stem, root, and husk. At the same time, the composition and structure of plant cell walls are ideally suited to the functions they perform [11]. The plant cell wall provides mechanical support to cells, tissues, and the entire plant body against lodging that might occur [13]. Sclerenchyma cells, which have both primary and thick secondary walls, provide the major mechanical support in the mature regions of the plant body [11].Si is one of the main elements in the Earth’s crust, second to oxygen in abundance [14]. Si is primarily present in the epidermal cells which provide structural rigidity to the plant. A remarkable unique feature of Si is that once it is deposited as silica gel, it is not redistributed to other parts of the plant [15]. Si accumulation sometimes exceeds that of crucial plant nutrients, especially in the grass family [16]. Rice is known as a Si accumulator [17], and the accumulation of Si in leaves and tissues can improve rice growth as a result of increasing the photosynthetic rate [18]. It can also enhance biomass partitioning, improve grain quality [9, 19], and increase yield [20]. Si is likely deposited in the cell walls of shoots, the intercellular regions, and silica cells [21]. A mature silica cell with observable Si deposition is referred to as a silica body. Si is deposited in two types of silica cells in the rice leaf to form two corresponding silica bodies, one dumbbell-shaped and the other bulliform-shaped [21]. The number of silica bodies in shoots correlates positively with the Si content in shoots. Si is also considered an environmentally friendly element in relation to soils, fertilisers, and plant nutrition [22].The method of fertiliser application is essential as it affects nutrient accessibility. Si was commonly applied as a base fertiliser [9] or foliage fertiliser [23]. Plants absorb Si most effectively when it is applied at the tillering to booting stages [24]. Interestingly, these methods are encouraging for nutrients that are taken up by means of vascular transportation. On the other hand, foliar is relatively easy to apply at manageable rates. However, the topdressing method for solid fertiliser could be supplied at higher rates [10]. In addition, these methods are favourable for nutrients that are taken up by means of vascular transportation. Moreover, the ease of application with Si could also make this form of application more viable and economical [10].The use of unmanned aerial vehicles (UAVs) for precision agriculture has been increasing, and these machines have been designed specifically to ease problematic agricultural activities [25]. Initially, UAVs were used to take aerial images of a farm [26]. A key advantage of using UAV is that this allows farmers to monitor and obtain valuable information about the field and crop conditions such as crop health, vigour, and growth stages and ultimately to determine yield [27].Thus, the main objective of this experiment was to examine the topdressing and foliar forms of Si applications on the physiological performance and structural changes of leaf and culm sheaths to enable an evaluation of the yield improvement and a cost analysis of the PadiU Putra rice line.


## 2. Materials and Methods

### 2.1. Field Experiment and Soil Properties

The field experiment was conducted in major rice granary areas from October 2018 until February 2019. The location of the study was at Kampung Maharaja, Tunjang, Kubang Pasu, Jitra, and Kedah under the Muda Agricultural Development Authority (MADA)(6° 16′ N 100° 21′ E, 76 m elevation). The soil has a clay loam texture (22.9% sand, 49.1% silt, and 28% clay) with pH 5.3 and was used as the cultivation medium. The soil nutrient status was 0.38% total N, 43 mg kg-1available P, 167 mg kg-1 available K, 1559 mg kg-1 available Ca, 433 mg kg-1 available Mg, 529 mg kg-1 available Fe, 23 mg kg-1 available Mn, 3.7 mg kg-1 available Zn, and 1.1 mg kg-1 available Cu.


### 2.2. Planting Materials and Plant Establishment

A transplanted blast-resistant and high-yielding rice line with high potential (MR219 X Pongsu Seribu 1 through SSR markers), known as PadiU Putra, was used in the experiment. The field was thoroughly prepared and levelled using a leveller machine before transplanting took place. Twenty-day-old seedlings were transplanted at20×15cm spacing. Nine plots for each experiment with 1,637 m2 (24.4×67.1m) areas were prepared and separated by a 0.5 m barrier. Each plot was fertilised with NPK compound fertiliser at 360 kg ha-1, 175 kg ha-1, and 175 kg ha-1 at 15, 50, and 70 days after transplanting (DAT). Urea was top-dressed at a rate of 100 kg ha-1 at 35 DAT. In each plot, a uniform plant stand was maintained, and standard agronomic practices were followed, as recommended by government agency practices.


### 2.3. Experimental Design and Statistical Analysis

In this experiment, treatments comprised two sources of methods of Si application (control treatments/untreated: without any additional Si source; topdressing treatments: plants supplied with Si in tablet form; and foliar applied treatments: plants supplied with Si through the foliar spray method) were arranged in a randomised complete block design (RCBD) with three replications. Si tablets were supplied by the Behn Meyer AgriCare (M) Company, Malaysia, at a rate of 3 kg ha-1 at 45 and 65 DAT, corresponding to the vegetative and reproductive stages, respectively. Two thousand and five hundred parts per million of Si foliar application with 0.01% dimethyl sulfoxide (DMSO) as an adjuvant was performed by spraying plants uniformly to the point of run-off (approximately 140 L ha-1) using a Stihl type mist blower (Stihl Sdn Bhd, Malaysia, 14 L) with a constant flow. The treatments were applied between 9:00 to 11:00 am on a day with clear sky at 45 and 65 DAT. Water was maintained at 5 cm by irrigation throughout the rice cultivation periods, and a 0.5 m boundary was prepared to avoid assortment. The Statistical Analysis System (SAS 9.2) by least significant different (LSD) at P≤0.05 was performed.


### 2.4. Growth Measurements

#### 2.4.1. Plant Height and Tiller Number

The plant height was measured according to the methods described by [28], where the measurement was made from the plant base to the tip of the highest leaf blade. The tiller number per square metre was counted by a fully expanded tiller. Samples with three replications were taken after ten days of Si treatment had been applied to the plants using all parameters. The purpose of this was to allow enough time for the plants to respond to the treatments. At the same time, two critical times in the rice growth stages were involved in this sampling time, the vegetative and reproductive stages. Four samples per square metre for each treatment were taken randomly and counted at 55 and 75 DAT.


#### 2.4.2. Average Diameter of Internodes and Diameter of Internodes at 20 cm

The average diameter of the internodes was measured at each internode midregion using a digital Vernier caliper from three representative plants from each treatment, with three replications, at 55 and 75 DAT. The diameter of the internodes at 20 cm was measured at the section of the internodes at 18 to 22 cm from the basal.


### 2.5. Photosynthesis, Stomatal Conductance, and Chlorophyll Content

The photosynthesis rate was measured on fully expanded young leaves (the third leaf from the top) at 0900-1100 am on a day with clear sky using a portable photosynthesis system (Li-6400XT, LI-COR, Lincoln, Nebraska, USA). The measurements were taken on the abaxial surface at a CO2 reference rate of 400 μmol m-2 s-1 at 55 and 75 DAT. The photosynthetic photon flux density (PPFD) was 900 mmol m-2 s-1. The stomatal conductance was derived from the same photosynthesis measurements described earlier. The chlorophyll content of the leaves was measured by the indirect method using a Portable Minolta SPAD 502 Plus chlorophyll meter (Delta T, UK). The third fully expanded leaf from the top was chosen for data measurements at 55 and 75 DAT. Three replications were taken for each of the parameters.


### 2.6. Biomass Partitioning

Three plants were harvested from each treatment at 55 and 75 DAT. They were partitioned into roots, culms, and leaves to determine the dry weights of each part. The dry weights of the plant parts were measured using a digital balance (QC35EDE-S Sartorius, Germany) after drying them in an oven at 72°C for three days until the weight became constant. The total biomass was calculated based on the total dry weights of the leaf, culm, root, and total yield. The root to shoot ratio was calculated using the following formula [29]: Root:ShootRatio=TotalRootDryWeight/TotalShootDryWeight.


### 2.7. Rice Crop Monitoring Using a Multirotor Unmanned Aerial Vehicle (UAV) and Red, Green, and (RGB) Digital Camera

For aerial imaging, a multirotor UAV DJI Phantom 4 Pro V2.0 with a gimbal-stabilised 4 K60 and 20 megapixel RGB digital camera attached was used to fly above the experimental field area. The flight plan was designed before the data acquisition by using theDroneDeploy software on a tablet. The altitude for the data acquisition of RGB images was set at 80 m, equivalent to a 2.19 cm spatial resolution. To avoid biased colours and lighting due to cloud shadows, the data collection was conducted in the morning under clear skies and at low wind speed conditions between 08:30 am and 11:00 am local time at 55 and 75 DAT. The camera settings were adjusted according to the light conditions and set to a fixed exposure for each flight.


### 2.8. Image Processing

Agisoft Metashape Professional software (http://www.agisoft.com/) was used to develop and align the imagery mosaic using Structure from Motion (SfM) algorithms. For each set of images, Agisoft PhotoScan software aligned the images and built point cloud models of the surface. Agisoft allows the generation and visualisation of a dense point cloud model, based on the estimated camera positions, to combine into a single dense point cloud [30]. The software provides a user-friendly process for mosaicking the imagery. The imagery was added and aligned using the Align Photo function. Then, the imagery generated and visualised a dense point cloud model based on the estimated camera position using the Build Dense Cloud function. It calculated the depth information for each camera, which could be combined into a single dense point cloud [30]. The geometrics of the map were reconstructed due to the poor texture of some elements of the scene and noisy or poorly focused images (known as outliers among the points) by using the Build Mesh function. The images were used to build the texture exported as a mosaicked orthophoto image [30]. Finally, the mosaicked orthophoto generated a Digital Surface Model (DSM), and the DSM and orthoimage were imported to build the 2.5 digital models.


### 2.9. Scanning Electron Microscopy Analysis of Leaf and Culm

At 55 and 75 DAT, the middle section of the flag leaves (approximately 0.5 cm in width) and the section of internodes at 18 to 22 cm from the basal (approximately 0.3 cm in width) from each treatment were cut with a sharp blade and fixed in fixative (4% glutaraldehyde) for 2 days at 4°C. The samples were washed in a buffer (0.1 M sodium cacodylate) with three changes of 30 minutes each. The samples were postfixed in osmium tetraoxide for 2 hours at 4°C and were washed again in the buffer (0.1 M sodium cacodylate) with three changes of 30 minutes each. The samples were then dehydrated in each graded acetone series of 35, 50, 75, and 95% (30 minutes) and 100% (1 hour each for three changes) followed by critical point drying (Baltec CPD 030). The leaves and culm specimens were mounted on aluminium stubs, covered with double-sided adhesive tape, and sputter-coated with gold with a sputter coater (Baltec SCD 005). They were viewed under a scanning electron microscope (JEOL JSM 6400, Japan) at an accelerating voltage of 15 kV. Three replications were taken for each of the parameters.


### 2.10. Yield Components, NPV, and BCR

After 80% grain maturity, all plants were harvested from each plot of treatments in the field trials experiment. Sampling was carried out by harvesting plants at maturity within a quadrant of 1 m2 to determine the yield components. The panicle number of each collected plant was counted to calculate the total number of panicles from the unit area (1 m2). Panicle per hill, grain number per panicle, and percentage of filled grains per panicle were counted and calculated manually. The thousand-grain weight (g) was also obtained using an electronic balance (QC 35EDE-S Sartorius, Germany). Ten panicle-bearing tillers from each treatment were sampled. Prior to weighing the grains, fully filled grains were manually separated from the unfilled grains. The percentage of filled grains per panicle was derived from the ratio of the number of fully ripened grains (filled grains) to the total number of grains per panicle per average hill (Yoshida 1981).The yield was determined with three replications for each treatment at the end of the experiment at 120 DAT. The yield was separately harvested by a harvester machine (World Star, WS70 Plus, China), loaded into separate lorries, and weighed in the rice mills (BERNAS Sdn Bhd). The yield of rice divided by the area in hectares, noted on the purchasing receipt, was collected from the rice mills. Moreover, the Malaysian government gives an additional incentive for every metric ton of rice yield increase, compared with the usual yield per hectare [31]. To determine the economic feasibility of both methods of Si application, a benefit-cost analysis and a net present revenue were carried out. The benefit-cost ratio (BCR) was calculated using the method of [32], whereas the net present value (NPV) was determined as described by [33].
(1)Thebenefit−costratioBCR=TCTR,NetpresentvalueNPV=ΣTR−TCt1+rtnt=1,where TC is the total cost, TR is the total revenue, and r is the discount rate of rice field for per season (t). If NR or NPV>0, then the total revenue is greater than the total cost; if NR or NPV=0, then the total revenue is equal to the total cost; and if the NR or NPV<0, then the total revenue is less than the total cost. In this study, NR and NPV are measured in Malaysian Ringgit (RM) and are based on one hectare. If BCR>1, then the total revenue is greater than the total cost; if BCR=1, then the gross revenue is equal to the total cost; and if the BCR<1, then the gross revenue is less than the total cost.


## 2.1. Field Experiment and Soil Properties

The field experiment was conducted in major rice granary areas from October 2018 until February 2019. The location of the study was at Kampung Maharaja, Tunjang, Kubang Pasu, Jitra, and Kedah under the Muda Agricultural Development Authority (MADA)(6° 16′ N 100° 21′ E, 76 m elevation). The soil has a clay loam texture (22.9% sand, 49.1% silt, and 28% clay) with pH 5.3 and was used as the cultivation medium. The soil nutrient status was 0.38% total N, 43 mg kg-1available P, 167 mg kg-1 available K, 1559 mg kg-1 available Ca, 433 mg kg-1 available Mg, 529 mg kg-1 available Fe, 23 mg kg-1 available Mn, 3.7 mg kg-1 available Zn, and 1.1 mg kg-1 available Cu.


## 2.2. Planting Materials and Plant Establishment

A transplanted blast-resistant and high-yielding rice line with high potential (MR219 X Pongsu Seribu 1 through SSR markers), known as PadiU Putra, was used in the experiment. The field was thoroughly prepared and levelled using a leveller machine before transplanting took place. Twenty-day-old seedlings were transplanted at20×15cm spacing. Nine plots for each experiment with 1,637 m2 (24.4×67.1m) areas were prepared and separated by a 0.5 m barrier. Each plot was fertilised with NPK compound fertiliser at 360 kg ha-1, 175 kg ha-1, and 175 kg ha-1 at 15, 50, and 70 days after transplanting (DAT). Urea was top-dressed at a rate of 100 kg ha-1 at 35 DAT. In each plot, a uniform plant stand was maintained, and standard agronomic practices were followed, as recommended by government agency practices.


## 2.3. Experimental Design and Statistical Analysis

In this experiment, treatments comprised two sources of methods of Si application (control treatments/untreated: without any additional Si source; topdressing treatments: plants supplied with Si in tablet form; and foliar applied treatments: plants supplied with Si through the foliar spray method) were arranged in a randomised complete block design (RCBD) with three replications. Si tablets were supplied by the Behn Meyer AgriCare (M) Company, Malaysia, at a rate of 3 kg ha-1 at 45 and 65 DAT, corresponding to the vegetative and reproductive stages, respectively. Two thousand and five hundred parts per million of Si foliar application with 0.01% dimethyl sulfoxide (DMSO) as an adjuvant was performed by spraying plants uniformly to the point of run-off (approximately 140 L ha-1) using a Stihl type mist blower (Stihl Sdn Bhd, Malaysia, 14 L) with a constant flow. The treatments were applied between 9:00 to 11:00 am on a day with clear sky at 45 and 65 DAT. Water was maintained at 5 cm by irrigation throughout the rice cultivation periods, and a 0.5 m boundary was prepared to avoid assortment. The Statistical Analysis System (SAS 9.2) by least significant different (LSD) at P≤0.05 was performed.


## 2.4. Growth Measurements

### 2.4.1. Plant Height and Tiller Number

The plant height was measured according to the methods described by [28], where the measurement was made from the plant base to the tip of the highest leaf blade. The tiller number per square metre was counted by a fully expanded tiller. Samples with three replications were taken after ten days of Si treatment had been applied to the plants using all parameters. The purpose of this was to allow enough time for the plants to respond to the treatments. At the same time, two critical times in the rice growth stages were involved in this sampling time, the vegetative and reproductive stages. Four samples per square metre for each treatment were taken randomly and counted at 55 and 75 DAT.


### 2.4.2. Average Diameter of Internodes and Diameter of Internodes at 20 cm

The average diameter of the internodes was measured at each internode midregion using a digital Vernier caliper from three representative plants from each treatment, with three replications, at 55 and 75 DAT. The diameter of the internodes at 20 cm was measured at the section of the internodes at 18 to 22 cm from the basal.


## 2.4.1. Plant Height and Tiller Number

The plant height was measured according to the methods described by [28], where the measurement was made from the plant base to the tip of the highest leaf blade. The tiller number per square metre was counted by a fully expanded tiller. Samples with three replications were taken after ten days of Si treatment had been applied to the plants using all parameters. The purpose of this was to allow enough time for the plants to respond to the treatments. At the same time, two critical times in the rice growth stages were involved in this sampling time, the vegetative and reproductive stages. Four samples per square metre for each treatment were taken randomly and counted at 55 and 75 DAT.


## 2.4.2. Average Diameter of Internodes and Diameter of Internodes at 20 cm

The average diameter of the internodes was measured at each internode midregion using a digital Vernier caliper from three representative plants from each treatment, with three replications, at 55 and 75 DAT. The diameter of the internodes at 20 cm was measured at the section of the internodes at 18 to 22 cm from the basal.


## 2.5. Photosynthesis, Stomatal Conductance, and Chlorophyll Content

The photosynthesis rate was measured on fully expanded young leaves (the third leaf from the top) at 0900-1100 am on a day with clear sky using a portable photosynthesis system (Li-6400XT, LI-COR, Lincoln, Nebraska, USA). The measurements were taken on the abaxial surface at a CO2 reference rate of 400 μmol m-2 s-1 at 55 and 75 DAT. The photosynthetic photon flux density (PPFD) was 900 mmol m-2 s-1. The stomatal conductance was derived from the same photosynthesis measurements described earlier. The chlorophyll content of the leaves was measured by the indirect method using a Portable Minolta SPAD 502 Plus chlorophyll meter (Delta T, UK). The third fully expanded leaf from the top was chosen for data measurements at 55 and 75 DAT. Three replications were taken for each of the parameters.


## 2.6. Biomass Partitioning

Three plants were harvested from each treatment at 55 and 75 DAT. They were partitioned into roots, culms, and leaves to determine the dry weights of each part. The dry weights of the plant parts were measured using a digital balance (QC35EDE-S Sartorius, Germany) after drying them in an oven at 72°C for three days until the weight became constant. The total biomass was calculated based on the total dry weights of the leaf, culm, root, and total yield. The root to shoot ratio was calculated using the following formula [29]: Root:ShootRatio=TotalRootDryWeight/TotalShootDryWeight.


## 2.7. Rice Crop Monitoring Using a Multirotor Unmanned Aerial Vehicle (UAV) and Red, Green, and (RGB) Digital Camera

For aerial imaging, a multirotor UAV DJI Phantom 4 Pro V2.0 with a gimbal-stabilised 4 K60 and 20 megapixel RGB digital camera attached was used to fly above the experimental field area. The flight plan was designed before the data acquisition by using theDroneDeploy software on a tablet. The altitude for the data acquisition of RGB images was set at 80 m, equivalent to a 2.19 cm spatial resolution. To avoid biased colours and lighting due to cloud shadows, the data collection was conducted in the morning under clear skies and at low wind speed conditions between 08:30 am and 11:00 am local time at 55 and 75 DAT. The camera settings were adjusted according to the light conditions and set to a fixed exposure for each flight.


## 2.8. Image Processing

Agisoft Metashape Professional software (http://www.agisoft.com/) was used to develop and align the imagery mosaic using Structure from Motion (SfM) algorithms. For each set of images, Agisoft PhotoScan software aligned the images and built point cloud models of the surface. Agisoft allows the generation and visualisation of a dense point cloud model, based on the estimated camera positions, to combine into a single dense point cloud [30]. The software provides a user-friendly process for mosaicking the imagery. The imagery was added and aligned using the Align Photo function. Then, the imagery generated and visualised a dense point cloud model based on the estimated camera position using the Build Dense Cloud function. It calculated the depth information for each camera, which could be combined into a single dense point cloud [30]. The geometrics of the map were reconstructed due to the poor texture of some elements of the scene and noisy or poorly focused images (known as outliers among the points) by using the Build Mesh function. The images were used to build the texture exported as a mosaicked orthophoto image [30]. Finally, the mosaicked orthophoto generated a Digital Surface Model (DSM), and the DSM and orthoimage were imported to build the 2.5 digital models.


## 2.9. Scanning Electron Microscopy Analysis of Leaf and Culm

At 55 and 75 DAT, the middle section of the flag leaves (approximately 0.5 cm in width) and the section of internodes at 18 to 22 cm from the basal (approximately 0.3 cm in width) from each treatment were cut with a sharp blade and fixed in fixative (4% glutaraldehyde) for 2 days at 4°C. The samples were washed in a buffer (0.1 M sodium cacodylate) with three changes of 30 minutes each. The samples were postfixed in osmium tetraoxide for 2 hours at 4°C and were washed again in the buffer (0.1 M sodium cacodylate) with three changes of 30 minutes each. The samples were then dehydrated in each graded acetone series of 35, 50, 75, and 95% (30 minutes) and 100% (1 hour each for three changes) followed by critical point drying (Baltec CPD 030). The leaves and culm specimens were mounted on aluminium stubs, covered with double-sided adhesive tape, and sputter-coated with gold with a sputter coater (Baltec SCD 005). They were viewed under a scanning electron microscope (JEOL JSM 6400, Japan) at an accelerating voltage of 15 kV. Three replications were taken for each of the parameters.


## 2.10. Yield Components, NPV, and BCR

After 80% grain maturity, all plants were harvested from each plot of treatments in the field trials experiment. Sampling was carried out by harvesting plants at maturity within a quadrant of 1 m2 to determine the yield components. The panicle number of each collected plant was counted to calculate the total number of panicles from the unit area (1 m2). Panicle per hill, grain number per panicle, and percentage of filled grains per panicle were counted and calculated manually. The thousand-grain weight (g) was also obtained using an electronic balance (QC 35EDE-S Sartorius, Germany). Ten panicle-bearing tillers from each treatment were sampled. Prior to weighing the grains, fully filled grains were manually separated from the unfilled grains. The percentage of filled grains per panicle was derived from the ratio of the number of fully ripened grains (filled grains) to the total number of grains per panicle per average hill (Yoshida 1981).The yield was determined with three replications for each treatment at the end of the experiment at 120 DAT. The yield was separately harvested by a harvester machine (World Star, WS70 Plus, China), loaded into separate lorries, and weighed in the rice mills (BERNAS Sdn Bhd). The yield of rice divided by the area in hectares, noted on the purchasing receipt, was collected from the rice mills. Moreover, the Malaysian government gives an additional incentive for every metric ton of rice yield increase, compared with the usual yield per hectare [31]. To determine the economic feasibility of both methods of Si application, a benefit-cost analysis and a net present revenue were carried out. The benefit-cost ratio (BCR) was calculated using the method of [32], whereas the net present value (NPV) was determined as described by [33].
(1)Thebenefit−costratioBCR=TCTR,NetpresentvalueNPV=ΣTR−TCt1+rtnt=1,where TC is the total cost, TR is the total revenue, and r is the discount rate of rice field for per season (t). If NR or NPV>0, then the total revenue is greater than the total cost; if NR or NPV=0, then the total revenue is equal to the total cost; and if the NR or NPV<0, then the total revenue is less than the total cost. In this study, NR and NPV are measured in Malaysian Ringgit (RM) and are based on one hectare. If BCR>1, then the total revenue is greater than the total cost; if BCR=1, then the gross revenue is equal to the total cost; and if the BCR<1, then the gross revenue is less than the total cost.


## 3. Results

### 3.1. Plant Growth Parameters

#### 3.1.1. Plant Height and Tiller Number

The plant height and tiller number per square metre for different Si method applications at 55 and 75 DAT in rice plants are shown in Table1. The results showed that the use of Si foliar applied and topdressing treatments showed significantly reduced plant height by 17 and 13% at 55 DAT, respectively, compared to the control. The results further revealed that both methods of Si treatments imposed caused slight differences in plant height during the reproductive stages of the rice plant (75 DAT). According to the results, the tiller number of both methods of Si treatments was not statistically different from the control at both DAT.Table 1
Plant height and tiller number on different Si method applications at 55 and 75 DAT in rice plants.
TreatmentPlant height (cm)Tiller number (m-2)DAT55755575Control79.7a84.5a438a513aTopdressing66.3b83.0a429a460aFoliar spray69.3b81.3a379a432aLSD(P=0.05)6.844.5517.0517.59CV4.212.424.9912.00Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.


#### 3.1.2. Diameter and Length of Internodes

The diameter of internodes was not statistically significant among treatments at 55 DAT (as shown in Table2). The diameter of internodes was significantly different at 75 DAT in the control, followed by the Si foliar applied and topdressing treatments, with 4.9, 4.0, and 3.5 mm, respectively. On the other hand, the length of internodes for Si topdressing treatments was significantly different, by 22%, at 55 DAT compared to the other treatments. The results revealed a reduction in the internode length at 75 DAT in both methods of Si treatments. Compared to the control, the reductions were by 14% in the topdressing and 9% in the foliar applied treatments.Table 2
Diameter and length of internodes on different Si method applications at 55 and 75 DAT in rice plants.
TreatmentDiameter of internodes at 20 cm (mm)Length of internodes (mm)DAT55755575Control2.72a4.9a35.83b63.74aTopdressing2.66a3.5c43.85a55.99bFoliar spray2.04a4.0b30.51c58.27bLSD(P=0.05)0.750.455.383.52CV13.314.737.732.62Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.


### 3.2. Physiological Responses

The photosynthesis rate (Pn) on different Si method applications at 55 and 75 DAT in rice plants is presented in Table3. The photosynthesis rates of plants treated with the Si topdressed and Si foliar spray methods were significantly different at 55 DAT, with the values of 11.3 and 11.0 μmol CO2 m-2 s-1, respectively. Moreover, both Si treatments increased the Pn by 30 (Si foliar applied methods) and 26% (Si topdressed methods).Table 3
Photosynthesis rate, stomatal conductance, and relative chlorophyll content from different Si method applications at 55 and 75 DAT in rice plants.
TreatmentPhotosynthesis rate (μmol CO2 m-2 s-1)Stomatal conductance (mmol H2O m-2 s-1)Relative chlorophyll content (SPAD value)DAT557555755575Control8.7b17.8a111.0b314.7a29.5b34.6aTopdressing11.3a16.2b183.0a269.0a26.6b31.2aFoliar spray11.0a15.8b102.9b218.5a32.5a35.6aLSD(P=0.05)2.031.3412.5818.165.233.53CV8.693.554.1912.893.614.61Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.Stomatal conductance measurements shared similarities with the photosynthesis rate (as shown in Table3). Plants treated with Si topdressing treatments were significantly different in stomatal conductance at 55 DAT, with a value of 183 mmol H2O m-2 s-1, or by a 65% increase compared to other treatments.The relative chlorophyll content of different Si method applications at 55 and 75 DAT in rice plants is shown in Table3. Si foliar-treated plants increased in chlorophyll content, with the highest mean value of 32.52 SPAD unit at 55 DAT compared to the control. Likewise, the relative chlorophyll content showed no significant difference at 75 DAT.


### 3.3. Biomass Partitioning

Figure1 and Table 4 illustrate the pattern of biomass partitioning in PadiU Putra rice plants with different Si method applications at 55 and 75 DAT. At 55 DAT, a statistical analysis revealed that Si foliar applied and Si topdressing treatments resulted in biomass partitioning to the leaves, culm, and roots being significantly different compared to the control (Table 4). Remarkably, Si foliar applied treatments were significantly different in biomass partitioning to the leaves, culm, and roots, following samples at 75 DAT. As for the total biomass partitioning, plants treated with Si foliar spray showed the highest dry matter production per hill at 55 and 75 DAT, where 43.4 and 63.2 g hill-1 were observed. Similar to the results from the root to shoot ratio, the application of Si foliar applied treatments showed the highest value compared to other treatments, with 1.41. Meanwhile, the application of Si topdressing treatments had the highest root to shoot ratio at 75 DAT, with 0.97.Figure 1
Whole plant percentage dry matter partitioning from different Si method applications at 55 and 75 DAT in rice plants.Table 4
Dry matter production on different Si method applications at 55 and 75 DAT in rice plants.
Biomass partitioningLeaf (g hill-1)Culm (g hill-1)Root (g hill-1)Total (g hill-1)Root : shoot ratioDAT55755575557555755575Control2.9b16.9a4.7b14.6b8.7c22.3b16.3c53.8b0.49c0.72cTopdressing6.4a14.4b5.9b10.5c12.0b24.0b24.2b48.9b0.98b0.97aFoliar spray7.5a16.3a10.5a18.8a25.4a28.1a43.4a63.2a1.41a0.80bLSD(P=0.05)0.572.121.881.990.693.771.370.693.332.12CV7.988.625.677.657.646.0310.097.465.643.49Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.In whole plant percentage dry matter partitioning, the highest dry matter partitioning occurred with the leaves, by 27%, in the Si topdressing; culm, by 29% in the control; and roots, by 59% in the Si foliar applied treatments. These were observed at 55 DAT (Figure1). In contrast, biomass was partitioned to the leaves by 31% in the control; culm by 29% in the Si foliar applied, and roots by 49% in the Si topdressing treatments. These were observed at 75 DAT. The leaves and roots were most affected by the partitioning changes under both Si applied treatments, respectively, compared to other parts of the plant, especially at 55 DAT. In addition, at 75 DAT, Si foliar applied and Si topdressing treatments affected greater culm and root development of the rice plants compared to the control. According to the biomass partitioning, Si topdressing treatments improved dry matter partitioning into the leaves by 9% and 8% at 55 and 75 DAT, respectively, compared to the control. A similar pattern was observed in the Si foliar applied treatments, where 5% and 3% of dry matter partitioning into roots and culm, respectively, were enhanced, compared to the control treatments.


### 3.4. Crop Growth Map Using Unmanned Aerial Vehicle (UAV) and RGB Digital Camera

Figure2 shows the boundary of the study area in this research during the main planting season, made using a multirotor UAV at an altitude of 80 metres from the ground to monitor rice crop management and yield prediction, at 55 and 75 DAT. The red, green, and blue (RGB) image illustrates that the whole of the study area was about 1.42 ha. Based on Figures 2(a)–2(d), the rice plants showed uniformly as a dark green colour in the plot, indicating better plant growth performance at 55 DAT, especially from foliar applied treatments. Similar findings were observed at 75 DAT, where the rice plants consistently demonstrated a light yellowish colour in all plots. However, in the control treatment, darker greenish and yellow colours were spotted around the plot.Figure 2
Boundary of the study area in this research at Tunjang, Jitra, Kedah, Malaysia at (a) 55 and (b) 75 DAT. Boundary of each of the treatment plots (T1: control/untreated plants; T2: Si topdressing plants; and T3: Si foliar applied plants) at (c) 55 and (d) 75 DAT.


### 3.5. Scanning Electron Microscopy (SEM) Observations

#### 3.5.1. SEM Investigation of Rice Leaves at 55 and 75 DAT

In SEM observations, various structures were identified in the rice leaves, including silica bodies, ladder-like structures or dumbbell-shaped trichomes (fine outgrowths or appendages on plants), and stomas (as shown in Figure3). Essentially, identical structures were found on all treatments of the leaves, regardless of the sampling times.Figure 3
SEM image of the epidermal region of a rice plant leaf. The white granular areas of Si were detected on the abaxial leaf surfaces (a–c) at 55 DAT and (d–f) at 75 DAT at magnification 350x. LS: ladder-like structures; SB: silica bodies; TR: trichomes; ST: stomas.
(a)(b)(c)(d)(e)(f)Si foliar applied treatments were significantly different in terms of the number of trichomes per field of view at 75 DAT (as shown in Table5). Trichome distribution in Si foliar applied treatments was the highest by 29%, compared to other treatments. Similar findings were noted with the length of the ladder-like structures. The Si topdressing and Si foliar applied treatments were significantly different to the control at 75 DAT, with 9.95 and 9.97 μm, respectively. The leaf samples of plants treated with both Si treatments showed a 13% extralength ladder-like structure in a row of silica cells and small silica bodies all around the sample.Table 5
Number of trichomes, ladder-like structures, and length of ladder-like structures of rice leaves with different Si method applications at 55 and 75 DAT.
TreatmentNumber of trichomes per field of view (FOV)Number of ladder-like structure per field of view (FOV)Length of ladder-like structure (μm)DAT557555755575Control6a5b26a15a5.75a8.67bTopdressing5a4b26a16a4.09b9.95aFoliar spray8a7a24a14a3.49b9.97aLSD(P=0.05)4.161.153.831.731.582.87CV2.4512.502.796.394.257.39Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation. Number of trichomes and ladder-like structures (FOV) were counted at magnification 350x.The adaxial surfaces of the plants with Si topdressing treatments showed the highest number of silica bodies per view with 248 (18%) and 262 (24%), respectively, at both sampling times (Table6). However, Si foliar applied treatments and control treatments were slightly different in terms of the number of silica bodies per view, observed at 55 and 75 DAT. In another parameter, the distance between the adjacent silica bodies was significantly different in Si topdressing treatments (5.40 μm) at 55 DAT, compared to other treatments. Interestingly, plants applied with Si topdressing and Si foliar spray at 75 DAT showed significantly higher distances between adjacent silica bodies than the control; the former two treatments showed 6.37 and 6.32 μm, respectively.Table 6
Characteristics of rice leaves in terms of number of silica bodies and distances between silica bodies from different Si method applications at 55 and 75 DAT in rice plants.
TreatmentNumber of silica bodies per field of view (FOV)Distance between adjacent silica bodies (μm)DAT55755575Control211b211b4.49b5.22bTopdressing248a262a5.40a6.37aFoliar spray211b216b4.36b6.32aLSD(P=0.05)12.9311.961.121.77CV7.892.906,401.97Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation. Number of silica bodies (FOV) was counted at magnification 1000x.Table7 shows the stomatal characteristics of rice leaves with different Si method applications at 55 and 75 DAT in rice plants. The highest number of stomata per field of view was obtained when Si was applied as foliar spray and for Si topdressed plants at 55 and 75 DAT, respectively. Si foliar applied and Si topdressing treatments showed an increased number of stomata per view, by 50 and 19% at 55 and 75 DAT, respectively. As for the method of application, the stomata diameters of Si foliar applied and Si topdressing were significantly higher than the control at 75 DAT, at 17.40 μm. Similar results were found with stomatal width, where Si foliar applied (2.48 μm) and Si topdressing (2.36 μm) were significantly different than the control at 55 DAT. Remarkably, both Si treatments improved stomata diameter by 117% and stomata width by 39% in Si foliar applied treatments at 75 DAT. Si foliar applied at 75 DAT had the highest stomatal width, followed by Si topdressing treatments and the control, though they were insignificantly different. Stomatal width in rice plant leaves had a markedly increased photosynthesis rate in this study using the Si treatments. There was a positive linear relationship between stomatal width and the photosynthesis rate for all treatments, with correlation y=0.2948x−0.4951. The significant coefficients of 0.85 in the relationship showed that the highest stomatal width was obtained with a 5.63 μm photosynthesis rate in Si foliar applied treatments at 75 DAT.Table 7
Stomatal characteristics of rice leaves with different Si method applications at 55 and 75 DAT.
TreatmentNumber of stomata per field of view (FOV)Stomatal diameter (μm)Stomatal width (μm)DAT557555755575Control4b46b12.00a8.02b2.06b4.06bTopdressing5b57a12.18a17.40a2.36a4.08bFoliar spray8a47b10.80a17.40a2.48a5.63aLSD(P=0.05)2.1111.961.816.061.791.53CV8.632.405.508.325.199.10Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation. Number of stomata (FOV) was counted at magnification 1000x.


#### 3.5.2. SEM Investigation in Rice Culm at 55 and 75 DAT

To further identify the factors responsible for the effects of different Si method applications on the rice plant strength, the anatomical structure of the plants was analysed for each of the treatments with ×30 and ×200 magnification, as shown in Figure4. In the transverse sections, the characteristics of layer thickness and vascular bundles were observed.Figure 4
SEM image of the epidermal regions of rice plant culm sheaths. (a–c) at 55 DAT, (d–f) at 75 DAT at magnification 30x, and (g, h) at 75 DAT at magnification 13x. vb: vascular bundles.
(a)(b)(c)(d)(e)(f)(g)(h)(i)The diameter of the culm was significantly different in both Si method applications (Table8). The diameter of the culm of the topdressing treatments (4171 mm) was significantly higher than the Si foliar applied treatments and control plants at 55 DAT. Interestingly, similar findings relating to the width and middle width of the culm were obtained in Si topdressing treatments at 55 DAT, with 2888 and 2147 mm, respectively. The diameter and width of the culm were higher when Si was applied as topdressing; they increased by 33% and 9%, respectively, compared to the Si foliar applied treatments and the control. Like the middle with of the culm, Si topdressing treatments showed the broadest rice culm at 31%.Table 8
The diameter and width of the culm with different Si method applications at 55 and 75 DAT in rice plants.
TreatmentCulm diameter (mm)Culm width (mm)Culm middle diameter (mm)Culm middle width (mm)DAT5575557555755575Control3130b7927a2639b4888a2254a3678a1782a1942bTopdressing4171a5684b2888a3981b2147a2782b1571a2543aFoliar spray3307b7310a2146c4218a1653b3216a1261a767cLSD(P=0.05)5.647.404.995.703.146.8013.2112.00CV6.509.196.825.678.378.906.925.62Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.The culm wall thickness of the first layer was significantly different depending on the method of Si application (Table9). Si applied as topdressing treatments had the highest culm wall thickness of the first layer, followed by the Si foliar treatments and the control at 55 DAT. Meanwhile, Si foliar applied treatments were significantly higher than other treatments at 75 DAT. Likewise, the culm wall thickness of the second layer showed no significance at 55 and 75 DAT. In the Si topdressing and Si foliar applied treatments, the culm wall thickness of the first layer was the highest, with 9 and 24% at 55 and 75 DAT, respectively. The culm wall thickness of the third layer was significantly different for the Si foliar applied treatments, with 387 μm at 55 DAT (24%). The equation correlation between the culm wall thickness and the height of the plant was y=−1.9362x+461.79 with R2=0.5193. There was a negative linear relationship between them. The negative correlation indicated that the thickness of the culm increased, while the height of the plant decreased.Table 9
Culm wall thickness of rice plants with different Si method applications at 55 and 75 DAT.
TreatmentCulm wall thickness of the 1st layer (μm)Culm wall thickness of the 2nd layer (μm)Culm wall thickness of the 3rd layer (μm)DAT557555755575Control308b266c322a305a312b339aTopdressing337a295b338a293a295b320aFoliar spray294b330a324a286a387a344aLSD(P=0.05)10.718.0110.514.8712.7315.04CV2.456.307.783.457.503.98Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.The images and results showed the presence of vascular bundles in the culm crosssections of rice plants in the outer and inner layers with 200x of magnification, regardless of treatment, are shown in Figure5 and Table 10. However, there were differences in terms of the diameter and width of vascular bundles on both layers. In the control plants, the diameter and width of vascular bundles on the outer layers were widened more at early data collection (55 DAT) compared to other treatments. This was in contrast to Si foliar applied treatments, which showed significant differences in terms of diameter and width of the outer layer of vascular bundles and the inner layer of the diameter of vascular bundles with 161.5 (58%), 134.0 (181%), and 204.0 (80%) μm at 75 DAT, respectively.Figure 5
SEM image of the epidermal regions of rice plants. (a–c) At 55 DAT and (d–f) at 75 DAT at magnification 200x. L1: first layer; L2: second layer; Vb1: vascular bundle at first layer; Vb2: vascular bundle at second layer.
(a)(b)(c)(d)(e)(f)Table 10
Vascular bundle characteristics of rice culm with different Si method applications at 55 and 75 DAT.
TreatmentOuter layer diameter of vascular bundles (μm)Outer layer width of vascular bundles (μm)Inner layer diameter of vascular bundles (μm)Inner layer width of vascular bundles (μm)DAT5575557555755575Control79.8a102.0b58.8a47.7c59.5a113.6b48.8a121.0aTopdressing53.8b96.9b26.7b102.0b77.8a107.1b49.4a100.2bFoliar spray60.4b161.5a27.8b134.0a70.9a204.0a40.3a117.0abLSD(P=0.05)17.8825.9518.984.6119.4825.5332.2817.16CV7.959.556.989.767.805.768.807.91Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.


### 3.6. Yield Component at 120 DAT

The panicle length, number of grains per panicle, and percentage of filled grain were significantly different between treatments at 120 DAT (as shown in Table11). Following that, Si foliar applied and Si topdressing treatments had the highest panicle length, with 26.67 and 26.53 cm, respectively. Meanwhile, the number of grains per panicle was increased in Si foliar applied treatments by 14% compared to the control. Similar results were obtained with the grain filling percentage in Si foliar applied treatments, where an increase of 95.2% was obtained. On the other hand, Si topdressing treatments improved by 8% in terms of thousand-grain weight compared to other treatments. Interestingly, the yield production of rice significantly increased in Si foliar applied, followed by the Si topdressing treatments, by 53 and 39%, respectively, compared to the control. A relatively higher rice yield was associated with the Si foliar applied and Si topdressing treatments, with 6.8 and 6.15 t ha-1, compared to the national Malaysia average yield of 4.8 t ha-1.Table 11
Yield components of rice plants at 120 DAT with different Si method applications at 55 and 75 DAT.
TreatmentPanicle length (cm)Spikelets per panicleFilled grain (%)1,000-grain weight (g)Yield (tha-1)Control24.77b175b89.53b26.03b4.44cTopdressing26.53a175b91.87b28.13a6.15bFoliar spray26.67a200a95.20a26.33b6.80aLSD(P=0.05)3.335.888.522.126.84CV5.648.634.083.494.21Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. CV: coefficient of variation.


### 3.7. Economic Viability of Adopting Si Topdressing and Si Foliar Applied Treatments in Rice Plants

Table12 summarises the benefit-cost ratio (BCR) and net present value (NPV) of the different Si treatments used in rice cultivation with the Padiu Putra rice variety (Figures 6–8). The Si foliar applied and Si topdressing treatments produced higher yield and were more profitable during this planting season compared to the control. From an economic perspective (BCR), for every RM 1 spent in adopting the Si foliar applied and Si topdressing treatments, RM 1.59 and RM 1.47 were gained in return, respectively. In terms of NPV, the Si foliar applied and Si topdressing treatments were more economically viable compared to the control, producing an additional income of RM 2843.27 ha-1 and RM 2101 ha-1, respectively. The income generation of farmers improved by 486 and 359% with the Si foliar applied and Si topdressing treatments.Table 12
Effects of different Si method applications on benefit-cost ratio (BCR) and net present value (NPV) in rice plants at 120 DAT.
TreatmentBenefit-cost ratioNet present value (RM ha-1)Control1.09c558.04cTopdressing1.49b2749.95bFoliar spray1.63a3583.13aLSD(P=0.05)1.881.55CV5.672.45Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. CV: coefficient of variation.Figure 6
Financial analysis of paddy project using transplanting technique for one hectare of area (control).Figure 7
Financial analysis of paddy project using transplanting technique for one hectare of area (Si topdressing).Figure 8
Financial analysis of paddy project using transplanting technique for one hectare of area (Si foliar applied).


## 3.1. Plant Growth Parameters

### 3.1.1. Plant Height and Tiller Number

The plant height and tiller number per square metre for different Si method applications at 55 and 75 DAT in rice plants are shown in Table1. The results showed that the use of Si foliar applied and topdressing treatments showed significantly reduced plant height by 17 and 13% at 55 DAT, respectively, compared to the control. The results further revealed that both methods of Si treatments imposed caused slight differences in plant height during the reproductive stages of the rice plant (75 DAT). According to the results, the tiller number of both methods of Si treatments was not statistically different from the control at both DAT.Table 1
Plant height and tiller number on different Si method applications at 55 and 75 DAT in rice plants.
TreatmentPlant height (cm)Tiller number (m-2)DAT55755575Control79.7a84.5a438a513aTopdressing66.3b83.0a429a460aFoliar spray69.3b81.3a379a432aLSD(P=0.05)6.844.5517.0517.59CV4.212.424.9912.00Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.


### 3.1.2. Diameter and Length of Internodes

The diameter of internodes was not statistically significant among treatments at 55 DAT (as shown in Table2). The diameter of internodes was significantly different at 75 DAT in the control, followed by the Si foliar applied and topdressing treatments, with 4.9, 4.0, and 3.5 mm, respectively. On the other hand, the length of internodes for Si topdressing treatments was significantly different, by 22%, at 55 DAT compared to the other treatments. The results revealed a reduction in the internode length at 75 DAT in both methods of Si treatments. Compared to the control, the reductions were by 14% in the topdressing and 9% in the foliar applied treatments.Table 2
Diameter and length of internodes on different Si method applications at 55 and 75 DAT in rice plants.
TreatmentDiameter of internodes at 20 cm (mm)Length of internodes (mm)DAT55755575Control2.72a4.9a35.83b63.74aTopdressing2.66a3.5c43.85a55.99bFoliar spray2.04a4.0b30.51c58.27bLSD(P=0.05)0.750.455.383.52CV13.314.737.732.62Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.


## 3.1.1. Plant Height and Tiller Number

The plant height and tiller number per square metre for different Si method applications at 55 and 75 DAT in rice plants are shown in Table1. The results showed that the use of Si foliar applied and topdressing treatments showed significantly reduced plant height by 17 and 13% at 55 DAT, respectively, compared to the control. The results further revealed that both methods of Si treatments imposed caused slight differences in plant height during the reproductive stages of the rice plant (75 DAT). According to the results, the tiller number of both methods of Si treatments was not statistically different from the control at both DAT.Table 1
Plant height and tiller number on different Si method applications at 55 and 75 DAT in rice plants.
TreatmentPlant height (cm)Tiller number (m-2)DAT55755575Control79.7a84.5a438a513aTopdressing66.3b83.0a429a460aFoliar spray69.3b81.3a379a432aLSD(P=0.05)6.844.5517.0517.59CV4.212.424.9912.00Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.


## 3.1.2. Diameter and Length of Internodes

The diameter of internodes was not statistically significant among treatments at 55 DAT (as shown in Table2). The diameter of internodes was significantly different at 75 DAT in the control, followed by the Si foliar applied and topdressing treatments, with 4.9, 4.0, and 3.5 mm, respectively. On the other hand, the length of internodes for Si topdressing treatments was significantly different, by 22%, at 55 DAT compared to the other treatments. The results revealed a reduction in the internode length at 75 DAT in both methods of Si treatments. Compared to the control, the reductions were by 14% in the topdressing and 9% in the foliar applied treatments.Table 2
Diameter and length of internodes on different Si method applications at 55 and 75 DAT in rice plants.
TreatmentDiameter of internodes at 20 cm (mm)Length of internodes (mm)DAT55755575Control2.72a4.9a35.83b63.74aTopdressing2.66a3.5c43.85a55.99bFoliar spray2.04a4.0b30.51c58.27bLSD(P=0.05)0.750.455.383.52CV13.314.737.732.62Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.


## 3.2. Physiological Responses

The photosynthesis rate (Pn) on different Si method applications at 55 and 75 DAT in rice plants is presented in Table3. The photosynthesis rates of plants treated with the Si topdressed and Si foliar spray methods were significantly different at 55 DAT, with the values of 11.3 and 11.0 μmol CO2 m-2 s-1, respectively. Moreover, both Si treatments increased the Pn by 30 (Si foliar applied methods) and 26% (Si topdressed methods).Table 3
Photosynthesis rate, stomatal conductance, and relative chlorophyll content from different Si method applications at 55 and 75 DAT in rice plants.
TreatmentPhotosynthesis rate (μmol CO2 m-2 s-1)Stomatal conductance (mmol H2O m-2 s-1)Relative chlorophyll content (SPAD value)DAT557555755575Control8.7b17.8a111.0b314.7a29.5b34.6aTopdressing11.3a16.2b183.0a269.0a26.6b31.2aFoliar spray11.0a15.8b102.9b218.5a32.5a35.6aLSD(P=0.05)2.031.3412.5818.165.233.53CV8.693.554.1912.893.614.61Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.Stomatal conductance measurements shared similarities with the photosynthesis rate (as shown in Table3). Plants treated with Si topdressing treatments were significantly different in stomatal conductance at 55 DAT, with a value of 183 mmol H2O m-2 s-1, or by a 65% increase compared to other treatments.The relative chlorophyll content of different Si method applications at 55 and 75 DAT in rice plants is shown in Table3. Si foliar-treated plants increased in chlorophyll content, with the highest mean value of 32.52 SPAD unit at 55 DAT compared to the control. Likewise, the relative chlorophyll content showed no significant difference at 75 DAT.


## 3.3. Biomass Partitioning

Figure1 and Table 4 illustrate the pattern of biomass partitioning in PadiU Putra rice plants with different Si method applications at 55 and 75 DAT. At 55 DAT, a statistical analysis revealed that Si foliar applied and Si topdressing treatments resulted in biomass partitioning to the leaves, culm, and roots being significantly different compared to the control (Table 4). Remarkably, Si foliar applied treatments were significantly different in biomass partitioning to the leaves, culm, and roots, following samples at 75 DAT. As for the total biomass partitioning, plants treated with Si foliar spray showed the highest dry matter production per hill at 55 and 75 DAT, where 43.4 and 63.2 g hill-1 were observed. Similar to the results from the root to shoot ratio, the application of Si foliar applied treatments showed the highest value compared to other treatments, with 1.41. Meanwhile, the application of Si topdressing treatments had the highest root to shoot ratio at 75 DAT, with 0.97.Figure 1
Whole plant percentage dry matter partitioning from different Si method applications at 55 and 75 DAT in rice plants.Table 4
Dry matter production on different Si method applications at 55 and 75 DAT in rice plants.
Biomass partitioningLeaf (g hill-1)Culm (g hill-1)Root (g hill-1)Total (g hill-1)Root : shoot ratioDAT55755575557555755575Control2.9b16.9a4.7b14.6b8.7c22.3b16.3c53.8b0.49c0.72cTopdressing6.4a14.4b5.9b10.5c12.0b24.0b24.2b48.9b0.98b0.97aFoliar spray7.5a16.3a10.5a18.8a25.4a28.1a43.4a63.2a1.41a0.80bLSD(P=0.05)0.572.121.881.990.693.771.370.693.332.12CV7.988.625.677.657.646.0310.097.465.643.49Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.In whole plant percentage dry matter partitioning, the highest dry matter partitioning occurred with the leaves, by 27%, in the Si topdressing; culm, by 29% in the control; and roots, by 59% in the Si foliar applied treatments. These were observed at 55 DAT (Figure1). In contrast, biomass was partitioned to the leaves by 31% in the control; culm by 29% in the Si foliar applied, and roots by 49% in the Si topdressing treatments. These were observed at 75 DAT. The leaves and roots were most affected by the partitioning changes under both Si applied treatments, respectively, compared to other parts of the plant, especially at 55 DAT. In addition, at 75 DAT, Si foliar applied and Si topdressing treatments affected greater culm and root development of the rice plants compared to the control. According to the biomass partitioning, Si topdressing treatments improved dry matter partitioning into the leaves by 9% and 8% at 55 and 75 DAT, respectively, compared to the control. A similar pattern was observed in the Si foliar applied treatments, where 5% and 3% of dry matter partitioning into roots and culm, respectively, were enhanced, compared to the control treatments.


## 3.4. Crop Growth Map Using Unmanned Aerial Vehicle (UAV) and RGB Digital Camera

Figure2 shows the boundary of the study area in this research during the main planting season, made using a multirotor UAV at an altitude of 80 metres from the ground to monitor rice crop management and yield prediction, at 55 and 75 DAT. The red, green, and blue (RGB) image illustrates that the whole of the study area was about 1.42 ha. Based on Figures 2(a)–2(d), the rice plants showed uniformly as a dark green colour in the plot, indicating better plant growth performance at 55 DAT, especially from foliar applied treatments. Similar findings were observed at 75 DAT, where the rice plants consistently demonstrated a light yellowish colour in all plots. However, in the control treatment, darker greenish and yellow colours were spotted around the plot.Figure 2
Boundary of the study area in this research at Tunjang, Jitra, Kedah, Malaysia at (a) 55 and (b) 75 DAT. Boundary of each of the treatment plots (T1: control/untreated plants; T2: Si topdressing plants; and T3: Si foliar applied plants) at (c) 55 and (d) 75 DAT.


## 3.5. Scanning Electron Microscopy (SEM) Observations

### 3.5.1. SEM Investigation of Rice Leaves at 55 and 75 DAT

In SEM observations, various structures were identified in the rice leaves, including silica bodies, ladder-like structures or dumbbell-shaped trichomes (fine outgrowths or appendages on plants), and stomas (as shown in Figure3). Essentially, identical structures were found on all treatments of the leaves, regardless of the sampling times.Figure 3
SEM image of the epidermal region of a rice plant leaf. The white granular areas of Si were detected on the abaxial leaf surfaces (a–c) at 55 DAT and (d–f) at 75 DAT at magnification 350x. LS: ladder-like structures; SB: silica bodies; TR: trichomes; ST: stomas.
(a)(b)(c)(d)(e)(f)Si foliar applied treatments were significantly different in terms of the number of trichomes per field of view at 75 DAT (as shown in Table5). Trichome distribution in Si foliar applied treatments was the highest by 29%, compared to other treatments. Similar findings were noted with the length of the ladder-like structures. The Si topdressing and Si foliar applied treatments were significantly different to the control at 75 DAT, with 9.95 and 9.97 μm, respectively. The leaf samples of plants treated with both Si treatments showed a 13% extralength ladder-like structure in a row of silica cells and small silica bodies all around the sample.Table 5
Number of trichomes, ladder-like structures, and length of ladder-like structures of rice leaves with different Si method applications at 55 and 75 DAT.
TreatmentNumber of trichomes per field of view (FOV)Number of ladder-like structure per field of view (FOV)Length of ladder-like structure (μm)DAT557555755575Control6a5b26a15a5.75a8.67bTopdressing5a4b26a16a4.09b9.95aFoliar spray8a7a24a14a3.49b9.97aLSD(P=0.05)4.161.153.831.731.582.87CV2.4512.502.796.394.257.39Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation. Number of trichomes and ladder-like structures (FOV) were counted at magnification 350x.The adaxial surfaces of the plants with Si topdressing treatments showed the highest number of silica bodies per view with 248 (18%) and 262 (24%), respectively, at both sampling times (Table6). However, Si foliar applied treatments and control treatments were slightly different in terms of the number of silica bodies per view, observed at 55 and 75 DAT. In another parameter, the distance between the adjacent silica bodies was significantly different in Si topdressing treatments (5.40 μm) at 55 DAT, compared to other treatments. Interestingly, plants applied with Si topdressing and Si foliar spray at 75 DAT showed significantly higher distances between adjacent silica bodies than the control; the former two treatments showed 6.37 and 6.32 μm, respectively.Table 6
Characteristics of rice leaves in terms of number of silica bodies and distances between silica bodies from different Si method applications at 55 and 75 DAT in rice plants.
TreatmentNumber of silica bodies per field of view (FOV)Distance between adjacent silica bodies (μm)DAT55755575Control211b211b4.49b5.22bTopdressing248a262a5.40a6.37aFoliar spray211b216b4.36b6.32aLSD(P=0.05)12.9311.961.121.77CV7.892.906,401.97Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation. Number of silica bodies (FOV) was counted at magnification 1000x.Table7 shows the stomatal characteristics of rice leaves with different Si method applications at 55 and 75 DAT in rice plants. The highest number of stomata per field of view was obtained when Si was applied as foliar spray and for Si topdressed plants at 55 and 75 DAT, respectively. Si foliar applied and Si topdressing treatments showed an increased number of stomata per view, by 50 and 19% at 55 and 75 DAT, respectively. As for the method of application, the stomata diameters of Si foliar applied and Si topdressing were significantly higher than the control at 75 DAT, at 17.40 μm. Similar results were found with stomatal width, where Si foliar applied (2.48 μm) and Si topdressing (2.36 μm) were significantly different than the control at 55 DAT. Remarkably, both Si treatments improved stomata diameter by 117% and stomata width by 39% in Si foliar applied treatments at 75 DAT. Si foliar applied at 75 DAT had the highest stomatal width, followed by Si topdressing treatments and the control, though they were insignificantly different. Stomatal width in rice plant leaves had a markedly increased photosynthesis rate in this study using the Si treatments. There was a positive linear relationship between stomatal width and the photosynthesis rate for all treatments, with correlation y=0.2948x−0.4951. The significant coefficients of 0.85 in the relationship showed that the highest stomatal width was obtained with a 5.63 μm photosynthesis rate in Si foliar applied treatments at 75 DAT.Table 7
Stomatal characteristics of rice leaves with different Si method applications at 55 and 75 DAT.
TreatmentNumber of stomata per field of view (FOV)Stomatal diameter (μm)Stomatal width (μm)DAT557555755575Control4b46b12.00a8.02b2.06b4.06bTopdressing5b57a12.18a17.40a2.36a4.08bFoliar spray8a47b10.80a17.40a2.48a5.63aLSD(P=0.05)2.1111.961.816.061.791.53CV8.632.405.508.325.199.10Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation. Number of stomata (FOV) was counted at magnification 1000x.


### 3.5.2. SEM Investigation in Rice Culm at 55 and 75 DAT

To further identify the factors responsible for the effects of different Si method applications on the rice plant strength, the anatomical structure of the plants was analysed for each of the treatments with ×30 and ×200 magnification, as shown in Figure4. In the transverse sections, the characteristics of layer thickness and vascular bundles were observed.Figure 4
SEM image of the epidermal regions of rice plant culm sheaths. (a–c) at 55 DAT, (d–f) at 75 DAT at magnification 30x, and (g, h) at 75 DAT at magnification 13x. vb: vascular bundles.
(a)(b)(c)(d)(e)(f)(g)(h)(i)The diameter of the culm was significantly different in both Si method applications (Table8). The diameter of the culm of the topdressing treatments (4171 mm) was significantly higher than the Si foliar applied treatments and control plants at 55 DAT. Interestingly, similar findings relating to the width and middle width of the culm were obtained in Si topdressing treatments at 55 DAT, with 2888 and 2147 mm, respectively. The diameter and width of the culm were higher when Si was applied as topdressing; they increased by 33% and 9%, respectively, compared to the Si foliar applied treatments and the control. Like the middle with of the culm, Si topdressing treatments showed the broadest rice culm at 31%.Table 8
The diameter and width of the culm with different Si method applications at 55 and 75 DAT in rice plants.
TreatmentCulm diameter (mm)Culm width (mm)Culm middle diameter (mm)Culm middle width (mm)DAT5575557555755575Control3130b7927a2639b4888a2254a3678a1782a1942bTopdressing4171a5684b2888a3981b2147a2782b1571a2543aFoliar spray3307b7310a2146c4218a1653b3216a1261a767cLSD(P=0.05)5.647.404.995.703.146.8013.2112.00CV6.509.196.825.678.378.906.925.62Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.The culm wall thickness of the first layer was significantly different depending on the method of Si application (Table9). Si applied as topdressing treatments had the highest culm wall thickness of the first layer, followed by the Si foliar treatments and the control at 55 DAT. Meanwhile, Si foliar applied treatments were significantly higher than other treatments at 75 DAT. Likewise, the culm wall thickness of the second layer showed no significance at 55 and 75 DAT. In the Si topdressing and Si foliar applied treatments, the culm wall thickness of the first layer was the highest, with 9 and 24% at 55 and 75 DAT, respectively. The culm wall thickness of the third layer was significantly different for the Si foliar applied treatments, with 387 μm at 55 DAT (24%). The equation correlation between the culm wall thickness and the height of the plant was y=−1.9362x+461.79 with R2=0.5193. There was a negative linear relationship between them. The negative correlation indicated that the thickness of the culm increased, while the height of the plant decreased.Table 9
Culm wall thickness of rice plants with different Si method applications at 55 and 75 DAT.
TreatmentCulm wall thickness of the 1st layer (μm)Culm wall thickness of the 2nd layer (μm)Culm wall thickness of the 3rd layer (μm)DAT557555755575Control308b266c322a305a312b339aTopdressing337a295b338a293a295b320aFoliar spray294b330a324a286a387a344aLSD(P=0.05)10.718.0110.514.8712.7315.04CV2.456.307.783.457.503.98Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.The images and results showed the presence of vascular bundles in the culm crosssections of rice plants in the outer and inner layers with 200x of magnification, regardless of treatment, are shown in Figure5 and Table 10. However, there were differences in terms of the diameter and width of vascular bundles on both layers. In the control plants, the diameter and width of vascular bundles on the outer layers were widened more at early data collection (55 DAT) compared to other treatments. This was in contrast to Si foliar applied treatments, which showed significant differences in terms of diameter and width of the outer layer of vascular bundles and the inner layer of the diameter of vascular bundles with 161.5 (58%), 134.0 (181%), and 204.0 (80%) μm at 75 DAT, respectively.Figure 5
SEM image of the epidermal regions of rice plants. (a–c) At 55 DAT and (d–f) at 75 DAT at magnification 200x. L1: first layer; L2: second layer; Vb1: vascular bundle at first layer; Vb2: vascular bundle at second layer.
(a)(b)(c)(d)(e)(f)Table 10
Vascular bundle characteristics of rice culm with different Si method applications at 55 and 75 DAT.
TreatmentOuter layer diameter of vascular bundles (μm)Outer layer width of vascular bundles (μm)Inner layer diameter of vascular bundles (μm)Inner layer width of vascular bundles (μm)DAT5575557555755575Control79.8a102.0b58.8a47.7c59.5a113.6b48.8a121.0aTopdressing53.8b96.9b26.7b102.0b77.8a107.1b49.4a100.2bFoliar spray60.4b161.5a27.8b134.0a70.9a204.0a40.3a117.0abLSD(P=0.05)17.8825.9518.984.6119.4825.5332.2817.16CV7.959.556.989.767.805.768.807.91Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.


## 3.5.1. SEM Investigation of Rice Leaves at 55 and 75 DAT

In SEM observations, various structures were identified in the rice leaves, including silica bodies, ladder-like structures or dumbbell-shaped trichomes (fine outgrowths or appendages on plants), and stomas (as shown in Figure3). Essentially, identical structures were found on all treatments of the leaves, regardless of the sampling times.Figure 3
SEM image of the epidermal region of a rice plant leaf. The white granular areas of Si were detected on the abaxial leaf surfaces (a–c) at 55 DAT and (d–f) at 75 DAT at magnification 350x. LS: ladder-like structures; SB: silica bodies; TR: trichomes; ST: stomas.
(a)(b)(c)(d)(e)(f)Si foliar applied treatments were significantly different in terms of the number of trichomes per field of view at 75 DAT (as shown in Table5). Trichome distribution in Si foliar applied treatments was the highest by 29%, compared to other treatments. Similar findings were noted with the length of the ladder-like structures. The Si topdressing and Si foliar applied treatments were significantly different to the control at 75 DAT, with 9.95 and 9.97 μm, respectively. The leaf samples of plants treated with both Si treatments showed a 13% extralength ladder-like structure in a row of silica cells and small silica bodies all around the sample.Table 5
Number of trichomes, ladder-like structures, and length of ladder-like structures of rice leaves with different Si method applications at 55 and 75 DAT.
TreatmentNumber of trichomes per field of view (FOV)Number of ladder-like structure per field of view (FOV)Length of ladder-like structure (μm)DAT557555755575Control6a5b26a15a5.75a8.67bTopdressing5a4b26a16a4.09b9.95aFoliar spray8a7a24a14a3.49b9.97aLSD(P=0.05)4.161.153.831.731.582.87CV2.4512.502.796.394.257.39Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation. Number of trichomes and ladder-like structures (FOV) were counted at magnification 350x.The adaxial surfaces of the plants with Si topdressing treatments showed the highest number of silica bodies per view with 248 (18%) and 262 (24%), respectively, at both sampling times (Table6). However, Si foliar applied treatments and control treatments were slightly different in terms of the number of silica bodies per view, observed at 55 and 75 DAT. In another parameter, the distance between the adjacent silica bodies was significantly different in Si topdressing treatments (5.40 μm) at 55 DAT, compared to other treatments. Interestingly, plants applied with Si topdressing and Si foliar spray at 75 DAT showed significantly higher distances between adjacent silica bodies than the control; the former two treatments showed 6.37 and 6.32 μm, respectively.Table 6
Characteristics of rice leaves in terms of number of silica bodies and distances between silica bodies from different Si method applications at 55 and 75 DAT in rice plants.
TreatmentNumber of silica bodies per field of view (FOV)Distance between adjacent silica bodies (μm)DAT55755575Control211b211b4.49b5.22bTopdressing248a262a5.40a6.37aFoliar spray211b216b4.36b6.32aLSD(P=0.05)12.9311.961.121.77CV7.892.906,401.97Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation. Number of silica bodies (FOV) was counted at magnification 1000x.Table7 shows the stomatal characteristics of rice leaves with different Si method applications at 55 and 75 DAT in rice plants. The highest number of stomata per field of view was obtained when Si was applied as foliar spray and for Si topdressed plants at 55 and 75 DAT, respectively. Si foliar applied and Si topdressing treatments showed an increased number of stomata per view, by 50 and 19% at 55 and 75 DAT, respectively. As for the method of application, the stomata diameters of Si foliar applied and Si topdressing were significantly higher than the control at 75 DAT, at 17.40 μm. Similar results were found with stomatal width, where Si foliar applied (2.48 μm) and Si topdressing (2.36 μm) were significantly different than the control at 55 DAT. Remarkably, both Si treatments improved stomata diameter by 117% and stomata width by 39% in Si foliar applied treatments at 75 DAT. Si foliar applied at 75 DAT had the highest stomatal width, followed by Si topdressing treatments and the control, though they were insignificantly different. Stomatal width in rice plant leaves had a markedly increased photosynthesis rate in this study using the Si treatments. There was a positive linear relationship between stomatal width and the photosynthesis rate for all treatments, with correlation y=0.2948x−0.4951. The significant coefficients of 0.85 in the relationship showed that the highest stomatal width was obtained with a 5.63 μm photosynthesis rate in Si foliar applied treatments at 75 DAT.Table 7
Stomatal characteristics of rice leaves with different Si method applications at 55 and 75 DAT.
TreatmentNumber of stomata per field of view (FOV)Stomatal diameter (μm)Stomatal width (μm)DAT557555755575Control4b46b12.00a8.02b2.06b4.06bTopdressing5b57a12.18a17.40a2.36a4.08bFoliar spray8a47b10.80a17.40a2.48a5.63aLSD(P=0.05)2.1111.961.816.061.791.53CV8.632.405.508.325.199.10Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation. Number of stomata (FOV) was counted at magnification 1000x.


## 3.5.2. SEM Investigation in Rice Culm at 55 and 75 DAT

To further identify the factors responsible for the effects of different Si method applications on the rice plant strength, the anatomical structure of the plants was analysed for each of the treatments with ×30 and ×200 magnification, as shown in Figure4. In the transverse sections, the characteristics of layer thickness and vascular bundles were observed.Figure 4
SEM image of the epidermal regions of rice plant culm sheaths. (a–c) at 55 DAT, (d–f) at 75 DAT at magnification 30x, and (g, h) at 75 DAT at magnification 13x. vb: vascular bundles.
(a)(b)(c)(d)(e)(f)(g)(h)(i)The diameter of the culm was significantly different in both Si method applications (Table8). The diameter of the culm of the topdressing treatments (4171 mm) was significantly higher than the Si foliar applied treatments and control plants at 55 DAT. Interestingly, similar findings relating to the width and middle width of the culm were obtained in Si topdressing treatments at 55 DAT, with 2888 and 2147 mm, respectively. The diameter and width of the culm were higher when Si was applied as topdressing; they increased by 33% and 9%, respectively, compared to the Si foliar applied treatments and the control. Like the middle with of the culm, Si topdressing treatments showed the broadest rice culm at 31%.Table 8
The diameter and width of the culm with different Si method applications at 55 and 75 DAT in rice plants.
TreatmentCulm diameter (mm)Culm width (mm)Culm middle diameter (mm)Culm middle width (mm)DAT5575557555755575Control3130b7927a2639b4888a2254a3678a1782a1942bTopdressing4171a5684b2888a3981b2147a2782b1571a2543aFoliar spray3307b7310a2146c4218a1653b3216a1261a767cLSD(P=0.05)5.647.404.995.703.146.8013.2112.00CV6.509.196.825.678.378.906.925.62Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.The culm wall thickness of the first layer was significantly different depending on the method of Si application (Table9). Si applied as topdressing treatments had the highest culm wall thickness of the first layer, followed by the Si foliar treatments and the control at 55 DAT. Meanwhile, Si foliar applied treatments were significantly higher than other treatments at 75 DAT. Likewise, the culm wall thickness of the second layer showed no significance at 55 and 75 DAT. In the Si topdressing and Si foliar applied treatments, the culm wall thickness of the first layer was the highest, with 9 and 24% at 55 and 75 DAT, respectively. The culm wall thickness of the third layer was significantly different for the Si foliar applied treatments, with 387 μm at 55 DAT (24%). The equation correlation between the culm wall thickness and the height of the plant was y=−1.9362x+461.79 with R2=0.5193. There was a negative linear relationship between them. The negative correlation indicated that the thickness of the culm increased, while the height of the plant decreased.Table 9
Culm wall thickness of rice plants with different Si method applications at 55 and 75 DAT.
TreatmentCulm wall thickness of the 1st layer (μm)Culm wall thickness of the 2nd layer (μm)Culm wall thickness of the 3rd layer (μm)DAT557555755575Control308b266c322a305a312b339aTopdressing337a295b338a293a295b320aFoliar spray294b330a324a286a387a344aLSD(P=0.05)10.718.0110.514.8712.7315.04CV2.456.307.783.457.503.98Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.The images and results showed the presence of vascular bundles in the culm crosssections of rice plants in the outer and inner layers with 200x of magnification, regardless of treatment, are shown in Figure5 and Table 10. However, there were differences in terms of the diameter and width of vascular bundles on both layers. In the control plants, the diameter and width of vascular bundles on the outer layers were widened more at early data collection (55 DAT) compared to other treatments. This was in contrast to Si foliar applied treatments, which showed significant differences in terms of diameter and width of the outer layer of vascular bundles and the inner layer of the diameter of vascular bundles with 161.5 (58%), 134.0 (181%), and 204.0 (80%) μm at 75 DAT, respectively.Figure 5
SEM image of the epidermal regions of rice plants. (a–c) At 55 DAT and (d–f) at 75 DAT at magnification 200x. L1: first layer; L2: second layer; Vb1: vascular bundle at first layer; Vb2: vascular bundle at second layer.
(a)(b)(c)(d)(e)(f)Table 10
Vascular bundle characteristics of rice culm with different Si method applications at 55 and 75 DAT.
TreatmentOuter layer diameter of vascular bundles (μm)Outer layer width of vascular bundles (μm)Inner layer diameter of vascular bundles (μm)Inner layer width of vascular bundles (μm)DAT5575557555755575Control79.8a102.0b58.8a47.7c59.5a113.6b48.8a121.0aTopdressing53.8b96.9b26.7b102.0b77.8a107.1b49.4a100.2bFoliar spray60.4b161.5a27.8b134.0a70.9a204.0a40.3a117.0abLSD(P=0.05)17.8825.9518.984.6119.4825.5332.2817.16CV7.959.556.989.767.805.768.807.91Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. DAT: days after transplanting; CV: coefficient of variation.


## 3.6. Yield Component at 120 DAT

The panicle length, number of grains per panicle, and percentage of filled grain were significantly different between treatments at 120 DAT (as shown in Table11). Following that, Si foliar applied and Si topdressing treatments had the highest panicle length, with 26.67 and 26.53 cm, respectively. Meanwhile, the number of grains per panicle was increased in Si foliar applied treatments by 14% compared to the control. Similar results were obtained with the grain filling percentage in Si foliar applied treatments, where an increase of 95.2% was obtained. On the other hand, Si topdressing treatments improved by 8% in terms of thousand-grain weight compared to other treatments. Interestingly, the yield production of rice significantly increased in Si foliar applied, followed by the Si topdressing treatments, by 53 and 39%, respectively, compared to the control. A relatively higher rice yield was associated with the Si foliar applied and Si topdressing treatments, with 6.8 and 6.15 t ha-1, compared to the national Malaysia average yield of 4.8 t ha-1.Table 11
Yield components of rice plants at 120 DAT with different Si method applications at 55 and 75 DAT.
TreatmentPanicle length (cm)Spikelets per panicleFilled grain (%)1,000-grain weight (g)Yield (tha-1)Control24.77b175b89.53b26.03b4.44cTopdressing26.53a175b91.87b28.13a6.15bFoliar spray26.67a200a95.20a26.33b6.80aLSD(P=0.05)3.335.888.522.126.84CV5.648.634.083.494.21Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. CV: coefficient of variation.


## 3.7. Economic Viability of Adopting Si Topdressing and Si Foliar Applied Treatments in Rice Plants

Table12 summarises the benefit-cost ratio (BCR) and net present value (NPV) of the different Si treatments used in rice cultivation with the Padiu Putra rice variety (Figures 6–8). The Si foliar applied and Si topdressing treatments produced higher yield and were more profitable during this planting season compared to the control. From an economic perspective (BCR), for every RM 1 spent in adopting the Si foliar applied and Si topdressing treatments, RM 1.59 and RM 1.47 were gained in return, respectively. In terms of NPV, the Si foliar applied and Si topdressing treatments were more economically viable compared to the control, producing an additional income of RM 2843.27 ha-1 and RM 2101 ha-1, respectively. The income generation of farmers improved by 486 and 359% with the Si foliar applied and Si topdressing treatments.Table 12
Effects of different Si method applications on benefit-cost ratio (BCR) and net present value (NPV) in rice plants at 120 DAT.
TreatmentBenefit-cost ratioNet present value (RM ha-1)Control1.09c558.04cTopdressing1.49b2749.95bFoliar spray1.63a3583.13aLSD(P=0.05)1.881.55CV5.672.45Mean values followed by the same letters within a column are not significantly different atP≤0.05 by the LSD test. CV: coefficient of variation.Figure 6
Financial analysis of paddy project using transplanting technique for one hectare of area (control).Figure 7
Financial analysis of paddy project using transplanting technique for one hectare of area (Si topdressing).Figure 8
Financial analysis of paddy project using transplanting technique for one hectare of area (Si foliar applied).


## 4. Discussion

Rice lodging is associated with height reduction, as stems weaken when supporting the weight of the grain, causing it to fall over. This is becoming a significant factor in rice production and often results in significant yield losses. Si application treatment is one of the approaches to fortify the leaf and culm sheaths and ultimately enhance rice yield. Broadcasting techniques, such as topdressing and foliar spray, are the common fertilisers used by the farmers, as they cost less but are very effective.Reducing the plant height and internode length and the diameter of the internode are the first responsive events of plants to lodging and commonly regulated by Si amendments [34]. There were positive indications with Si application in terms of physiological measurement such as photosynthesis, chlorophyll content [12, 35], and biomass partitioning [36]. These physiological measurements were supported with images using the UAV and can be used for rice yield predictions [30, 37].Si foliar applied treatments that are capable of altering physiological and anatomical structures are preferred, for they could fortify the leaf and culm sheaths and ultimately increase rice yield production. The greater improvement of biomass partitioning to the leaf and root in the early stages using Si topdressing and foliar applied was followed by the culm and root being affected at the reproductive stages by both Si treatments, respectively. It also improved the biomass partitioning, mostly to the culm and root, through improvements in photosynthetic efficiency [38]. The assimilation rates in photosynthetic leaves increased due to enhanced photosynthetic metabolites and enzyme activity, which directly improved grain production. The ratio of dry matter accumulation in roots to shoots with Si foliar applied treatments was greater than in the control plants, which clearly indicated the altered partitioning of resources within the rice plants. This might be due to the higher translocation of nutrients in roots, which influenced the movement of ions from the soil to the plant leaves. Similar results were reported by [39–41].Using scanning electron microscopy (SEM), an extralong ladder-like structure and an intense number of silica bodies on the leaf structures for both Si treatments were observed. Silica content could be well-deposited in the leaf epidermis of rice plants and form a thick layer as a support structure and for the protection of the leaf [42, 43]. It was assumed that the consistently intense lignin accumulation in the leaf was for self-defence mechanism purposes; similar findings were observed by [12]. Interestingly, Si foliar applied treatments might alter the leaves stomatal width, thereby enhancing the stomatal aperture and improving photosynthetic efficiency for better dry matter accumulation [40, 41].The strengthened culm sheath structures in Si-treated plants might be due to the thickness of the cell walls in the sclerenchyma and more vascular bundles in both the peripheral and the inner sections of the outer layers [44]. Similar findings with both Si treatments were found according to [12, 45]. That indicates an efficient assimilated translocation process from root to shoot in Si treatments.Greater yield production with Si treatment showed that rice plants improve their photosynthesis and assimilated distribution [46, 47]. Spikelets per panicle and grain filling percentage were higher when Si was applied as foliar spray, compared to the other treatments. Compared to this parameter, the percentage of filled grain could potentially reveal more about the effects of Si supplementation as it is an indicator of sterility and fertility. Therefore, yield seems to be maximised when Si was supplied as foliar spray at the reproductive stage, which is also referred to as preheading [10]. Interestingly, BCR and NPV showed the better economic viability of adopting Si foliar applied treatments in rice plants. These results suggest the ability of Si foliar applied treatments to enhance farmers’ income generation tremendously [47–53]. Thus, the Si foliar spray is the most suitable method of application and could be incorporated in combination with fungicides and insecticides to aid the supply of nutrients to crops in the form of sprays [10].Further study should be conducted on the mechanism of Si foliar applied treatments on the plant uptake through the leaf and root; the enzymes and hormonal regulation related to yield improvement, such as sucrose synthase enzyme; and the molecular identification of genes associated with hormonal regulation.


## 5. Conclusions

Si foliar applied treatments affect the plant growth, physiological process, leaf and culm morphology, and the yield of rice. The adjustment of the leaves’ stomatal width, the enhancement of the stomatal aperture, the more vascular bundles in the culm, and the efficient photosynthesis rate lead to better dry matter accumulation from root to shoot with Si foliar applied treatments. Therefore, rice yield production increased tremendously with a higher benefit-cost ratio and net present value.



---

*Source: 6679787-2021-06-01.xml*