All methods included in the research, such as the experimental design, measurement of planting uniformity, yield, and sustainable performance indicators, are under the guidelines of the International Rice Research Institute (IRRI) or global standards, which are indicated in the specific sections and parameters below. The manuscript was internally reviewed and approved by IRRI.
Site and crop descriptions, experimental design and water management. The experiment was conducted in Trung-Thanh Village, Co Do District, Can Tho, Vietnam (10.178103 °N latitude; 105.524434 °E longitude), across two consecutive rice-cropping seasons. These were the Winter-Spring season (WS), or dry season, from 8 November 2018 (sowing) to 14 February 2019 (harvest); and the Summer-Autumn season (SA), or early wet season, from 1 March 2019 (sowing) to 28 May 2019 (harvest). Rice varieties used were OM5451 and Dai-Thom-8 for the WS and SA seasons, respectively. The use of these plants complies with the national guidelines of Vietnam32. At the start of the WS season, fields were drained of floodwater, whereas, at the start of the SA season, irrigation water was required during land preparation before the onset of the monsoon rains. The mean farm size in the study area was 2.1 ± 0.1 ha, with an acid sulphate clay soil type8. The predominant crop establishment method was wet direct-seeding with broadcast pregerminated seed and the majority of farmers used four-wheel tractors for land preparation and combine harvesters for harvesting8.
Table 1
Distribution of the treatments and replications in the experiment plots
Farmer 1
|
Farmer 2
|
Farmer 3
|
Farmer 4
|
BroadC-1
|
BlowS-2
|
BlowS-3
|
MecT-4
|
BlowS-1
|
MecT-2
|
MecT-3
|
DrumS-4
|
DrumS-1
|
BroadC-2
|
DrumS-3
|
BroadC-4
|
MecT-1
|
DrumS-2
|
BroadC-3
|
BlowS-4
|
BroadC = manual broadcasting, BlowS = blower-seeding, DrumS = drum-seeding, MecT = mechanized transplanting; the numbers associated with the treatments in the Table (i.e. 1, 2, 3, and 4) represents for the blocks or farmers, correspondingly.
The four crop establishment methods were considered as separate treatments and implemented in a randomized complete block design (Table 1). The four treatments were: 1) BroadC (Fig. 1a), 2) BlowS using the Kasei 3WF-3A-26L machine (Fig. 1b), 3) DrumS using the Hoang-Thang drum seeder (Fig. 1c), and 4) MecT using the Yanmar VP7D25 transplanter (Fig. 1d). The four different farmers’ fields were considered as blocks or replicates and the four crop establishment methods were applied in each field (block). The field area of each treatment ranged from 3,000 to 4,000 m2. Different fields were used each season. Irrigation and drainage were applied similarly across the four treatments, but were different for the WS and SA seasons depending on the weather and flood conditions at the research site (Fig. 2). Growing time of the rice was 90 and 83 days for DSR and MecT, respectively. However, MecT required the seedlings to be prepared 12 days prior to crop establishment.
Land preparation, fertilizer and pesticide application, and harvesting operations, were the same for all treatments. Land preparation included plowing using locally fabricated rotavators and puddling with wet leveling. Fertilizer and pesticides were applied using Knapsack-blowers and -sprayers, respectively, combined in the Kasei 3WF-3A-26L machine. Harvesting in all treatment plots was done using combine harvesters (Kubota DC-70). Straw after harvest was incorporated and burned in the field before WS and SA, respectively, as a common farmer practice at the research site. 1M5R was applied in all treatments following the criteria shown in Table 2.
Table 2
Specifications of best practices for irrigated rice production (1 Must (certified seed) and 5 Reductions (reduced rates of seed, fertilizer, pesticides and water; reduced post-harvest losses) applied in the field trial at Trung-Thanh Village, Co Do District, Can Tho, for both seasons. (Max. = maximum).
Criteria*
|
Requirements
|
|
Seed rate
|
≤120 kg ha−1
|
Certified seed
|
Nitrogen
|
≤100 kg ha−1
|
Applied with at least three splits
|
Insecticides
|
Max. 1 product application
|
No application within 40 days after sowing
|
Fungicides
|
Max. 2 product applications
|
No application after the flowering crop stage
|
Water management
|
Dry fields during the cultivation following AWD technique
|
Harvesting
|
Combine harvester
|
Harvest when 80–85% of the grains per panicle are straw or yellow-colored
|
*Postharvest processes were excluded for analyses of the findings in this study. |
Measuring planting uniformity. To measure the planting uniformity, five 40- x 50-cm quadrants (randomly placed in a cross diagonal transect) were sampled in each treatment plot 7 days after sowing or transplanting. Seedling density was assessed by counting the number of seedlings in each quadrant divided by the quadrant area. The standard deviation (SD) was then used to compare the variation in seedling density from the mean across all replicate plots.
Quantification of grain yield. Grain yield was determined by the crop-cut method from each experimental plot. In the WS season, the samples for the crop cut were taken from two 5-m2 (2.5 x 2.0 m) quadrants, which were located 5 m from the center of each plot along a cross diagonal transect. In the SA season, the same sampling procedure was applied with one more sampling at the center of each plot (total of three samples for each plot). The threshed paddy grains were cleaned (unfilled spikelets removed), weighed and recorded as fresh weight. The moisture content (MC) of the grain samples was determined using a grinding-type moisture meter (Kett®, product code: F511), which was precalibrated using the oven method33. The grain yield was calculated at 14% MC.
Analysis of energy efficiency and indirect GHGEs. Energy efficency (GJ ha−1) was analyzed based on the net differences between the outputs and inputs of rice production—Equation 1 (Eq. 1):
NEV = Eout – Ein (GJ ha-1) (Eq. 1)
Where NEV is the net energy value for energy efficiency; Eout is the output energy value only accounting for the harvested grains but not including rice straw because this residue was incorporated before WS and burned before SA in this research; Ein is the input energy value accounting for mechanized operations including machine production and fuel consumption, labor and agronomic inputs such as seeds, fertilizer and pesticide. The conversion factors reported in Ecoinvent (2019)34 were used to estimate the energy of the related materials and processes (Table 3). In addition, the energy conversion factor for machine production was calculated through fuel consumption at 15 MJ L-1 (35,36). Fertilizer inputs, such as nitrogen (N), phosphorus (P2O5) and potassium (K2O) were calculated based on the chemical content of N, P2O5 and K2O, such as urea (46-0-0) and DAP (18-46-0). Pesticide and herbicide inputs were converted based on the content of active ingredients and the conversion weight of the applied pesticides. Manpower was calculated based on the metabolic equivalent of task (MET), which is the ratio of human metabolic rate when performing an activity to the metabolic rate at rest, and on a labor energy conversion factor37, with the assumption that the mean weight of a Vietnamese is 55 kg.
GHGE (kgCO2-eq ha-1) is calculated based on Eq. 2, that accounts for the production of agronomic inputs including seeds, fertilizer and pesticide (GHGagro-input); mechanized operations (GHGoperation), soil emissions (GHGsoil) and rice straw management (GHGricestraw).
GHGE = GHGagro-inputs + GHGoperation + GHGsoil + GHGricestraw (kgCO2-eq ha-1 season-1) (Eq. 2)
The GHG conversion factors for agronomic inputs and mechanized operation are shown in Table 3. GHGsoil is calculated based on Eq. 338, accounting for CH4 and N2O emissions. The CH4 emission is affected by water management, pre-season soil management and rice straw incorporation; while the N2O emission is affected by N use for rice cultivation38.
GHG soil = Timegrow*28*EFdefault* SFwater* SFpre*SFricestraw + 265*EF1FR*Ffertilizer (kgCO2-eq ha-1 season-1) (Eq. 3)
Where Timegrow is the rice-growing period; 28 and 265 are the Global Warming Potentials of CH4 and N2O, respectively, for conversion to CO2-eq38; EFdefault, SFwater and SFpre, are the CH4 emission and scaling factors of water management and pre-season soil management, respectively; and SFricestraw is the scale factor for rice straw management. EF1FR is the N2O emission factor in flooded rice systems and fertilizer amount of applied N, calculated based on Eq. 438. Water management was considered as single- and multiple-drainage scenarios during the WS and SA, respectively (Fig. 2). Total growing time of the direct-seeded rice was 90 days while that of transplanted rice in the field was 83 days. The seedling preaparation time of 12 days was accounted for in the transplanted rice scenario. However, the land area used for seedling is only 1:100 for growing compared with the common practice in Vietnam, which was observed to be the case in this study. The emission and scaling factors are shown in Table 3.
SFricestraw = \({\left(1+{\text{R}}_{\text{s}\text{t}\text{r}\text{a}\text{w}}*\text{C}{\text{F}}_{\text{s}\text{t}\text{r}\text{a}\text{w}}\right)}^{0.59}\) (Eq. 4)
Where Rstraw is the incorporation rate of rice straw (dry matter, t ha-1) and CFstraw is the conversion factor of rice straw depending on time of incorporation before the crop establishment. Yield of straw only accounted for top parts of rice plant harvested is 50% of rice yield39. This factor is only applied for the straw incorporation scenario of WS but not for the burning scenario of SA. On the other hand, GHG emission from straw burning is taken into account through the last component (GHGricestraw) in Eq. 2 which is reported in Romasanta (2017)40.
Table 3
Energy and GHGE conversion factors used for calculating the relative energy efficiency of the four crop establishment methods from crop establishment to harvest.
Parameters
|
|
Energy
|
|
GHGE
|
|
|
Unit
|
Value
|
Source
|
|
Unit
|
Value
|
Sources
|
Land use
|
|
MJ ha−1
|
0.0024
|
34,41
|
|
See details under “Soil emissions”
|
Seeds
|
|
MJ kg−1
|
30.1
|
34,41
|
|
kgCO2-eq kg−1
|
1.12
|
34,41,46
|
Grain
|
|
MJ kg−1
|
15.2
|
34,42
|
|
|
|
|
Diesel consumption
|
|
MJ L−1
|
44.8
|
34,35,41
|
|
kgCO2-eq MJ−1
|
0.08
|
34,41,46
|
Gasoline consumption
|
|
MJ L−1
|
39.1
|
34,35,41
|
|
kgCO2-eq MJ−1
|
0.08
|
34,41,46
|
Electric power
|
|
MJ kWh−1
|
3.6
|
34,35,41
|
|
kgCO2-eq kWh−1
|
0.564
|
34,41
|
Machine production
|
|
MJ L−1
|
15.6
|
35,36
|
|
|
|
|
N
|
|
MJ kg−1
|
58.7
|
34,41,43
|
|
kgCO2-eq kg−1
|
5.68
|
34,41,46
|
P2O5
|
|
MJ kg−1
|
17.1
|
34,41,43
|
|
kgCO2-eq kg−1
|
1.09
|
34,41,46
|
K2O
|
|
MJ kg−1
|
8.83
|
34,41,43
|
|
kgCO2-eq kg−1
|
0.52
|
34,41,46
|
Herbicide
|
|
MJ kg−1
|
354
|
34,41,44
|
|
kgCO2-eq kg−1
|
23.3
|
34,41,46
|
Pesticide
|
|
MJ kg−1
|
182
|
34,41,43
|
|
kgCO2-eq kg−1
|
10.4
|
34,41,46
|
Driving 4WT and combine harvesters
|
|
MJ h−1
|
0.44
|
37,45
|
|
|
|
|
Manual labor
|
|
MJ h−1
|
0.89
|
37,45
|
|
|
|
|
Soil emission:
|
|
|
|
|
|
|
|
|
• EFdefault of CH4 in WS
|
|
|
|
kg ha−1 day−1
|
1.7
|
47
|
• EFdefault of CH4 in SA
|
|
|
|
kg ha−1 day−1
|
2.8
|
47
|
• SFpre for pre-season soil management
|
|
|
|
1
|
13
|
• SFwater for single drainage
|
|
|
|
|
0.71
|
38
|
• SFwater for multiple drainage
|
|
|
|
|
0.55
|
38
|
• SFN for Nitrogen use
|
|
|
|
% N applied
|
0–1
|
38
|
• CFincorporation
|
|
|
|
|
1
|
13
|
• CH4 from burning straw
|
|
|
|
kg Mg−1
|
4.51
|
40
|
• N2O from burning straw
|
|
|
|
kg Mg−1
|
0.069
|
40
|
Computation of sustainability performance Indicators. The Sustainable Rice Platform (SRP) has developed 12 sustainability performance indicators for rice production (SRP, 2019). We computed the seven agronomic indicators: productivity (grain yield), nitrogen-use efficiency (NUE), phosphorous-use efficiency (PUE), biodiversity (pesticide use), labor productivity, profitability (net profit) and GHGE as defined by SRP version 213. In addition, we included potassium-use efficiency (KUE) due to its importance in rice productivity. Farmers were asked to record input and economic data in diaries, which were checked and collected by project staff every 3–4 weeks. To compute phosphorus (P) and potassium (K) application rates, the amounts of P2O5 and K2O for each fertilizer application were determined and multiplied by a factor of 0.4364 and 0.8302, respectively, to convert them into the elemental form13. To compute for NUE, PUE and KUE, the total grain yield harvested was divided by the elemental N, P or K rate applied and was expressed in terms of kg grain kg−1 elemental N, P or K. To compare pesticide practices among treatments, we reported the total frequency of application of formulated pesticide products. To compute for labor productivity, both hired and owned (family) male and female laborers were considered and the number of labor days per season (for all activities from land preparation until harvest, including regular field visits by farmers) were estimated by dividing the total labor cost per season by the average daily wage rate (VND 200,000 day−1, collected during this research) at the time taken across all activities. The result was then divided by the grain yield as determined from crop cuts.
Net income was calculated by deducting the total production cost from the gross income obtained from grain yield. Production cost consisted of: 1) land use; 2) service costs of mechanized operations such as land preparation, mechanical transplanting, fertilizer and pesticide applications and combine harvesting; 3) agronomic inputs including seeds, fertilizer and pesticide; and, 4) labor. Gross income consisted of the income from the total fresh harvested grain sold at the field. Costs of inputs and price of paddy are shown in Table 4.
Table 4
Cost of inputs and price of paddy.
Inputs
|
|
Unit
|
Value
|
|
Land use
|
|
$US ha−1 year−1
|
2,000
|
|
Seed
|
|
$US kg−1
|
5.2
|
|
Urea 46-0-0
|
|
$US kg−1
|
58.7
|
|
TSP 18-46-0
|
|
$US kg−1
|
0.6
|
|
MOP 0-0-60
|
|
$US kg−1
|
0.4
|
|
NPK 16-16-16
|
|
$US kg−1
|
0.6
|
|
NPK 16-16-8
|
|
$US kg−1
|
0.5
|
|
Herbicide
|
|
$US L−1
|
4.8
|
|
Molluscicide
|
|
$US L−1
|
6.1
|
|
Fungicide/Insecticide
|
|
$US kg−1
|
12.3
|
|
Fungicide/Insecticide
|
|
$US L−1
|
11.0
|
|
Land preparation
|
|
$US ha−1
|
94.1
|
|
Manual broadcast-seeding
|
|
$US ha−1
|
26.0
|
|
Blower seeding
|
|
$US ha−1
|
26.0
|
|
Drum seeding
|
|
$US ha−1
|
26.0
|
|
Mechanized transplanting
|
|
$US ha−1
|
220.0
|
|
Crop care
|
|
$US ha−1
|
56.5
|
|
Harvesting
|
|
$US ha−1
|
90.3
|
|
Statistical analysis and software. SPSS software and Analysis of Variance (ANOVA) were used to evaluate the effects of the contrasting crop establishment-based scenarios on the measured production and environmental parameters using a Least Significant Difference (LSD) at α = 0.05 to compare the mean values. To compare the different treatments against FP, one-way ANOVA was used to analyze the sustainability performance indicators for each treatment. Each treatment was analyzed separately against FP because the treatments were applied in a block design within fields that were interspersed with the FP fields. Seedling density was analyzed using log transformation due to non-normally distributed residuals. Energy balance analysis was based on the Cumulative Energy Demand 1.09 method by SIMAPRO (2019)41 and CO2-equivalent analysis was based on the GWP-100a of IPCC (2013)46.