Estimation of Actual Evapotranspiration and Crop Coefﬁcient of Transplanted Puddled Rice Using a Modiﬁed Non-Weighing Paddy Lysimeter

: Lysimetric and eddy covariance techniques are commonly used to directly estimate actual crop evapotranspiration ( ET a ). However, these technologies are costly, laborious, and require skills which make in situ ET estimation difﬁcult, particularly in developing countries. With this in mind, an attempt was made to determine ET a and stagewise crop coefﬁcient (K c ) values of transplanted puddled rice using a modiﬁed non-weighing paddy lysimeter. The results were compared to indirect methods, viz., FAO Penman–Monteith and pan evaporation. Daily ET a ranged from 1.9 to 8.2 mmday − 1 , with a mean of 4.02 ± 1.35 mmday − 1 , and their comparison showed that the FAO Penman–Monteith equation performed well for the coefﬁcient of determination (R 2 of 0.63), root mean squared error (RMSE = 0.80), and mean absolute percentage error (MAPE = 13.6 %), and was highly correlated with ET a throughout the crop season. However, the pan evaporation approach was underestimated (R 2 of 0.24; RMSE = 0.98; MAPE = 22.13%) due to a consistent pan coefﬁcient value (0.71), vegetation role and measurement errors. In addition, actual K c values were obtained as 1.13 ± 0.13, 1.27 ± 0.2, 1.23 ± 0.16, and 0.93 ± 0.18 for the initial, crop development, mid-season, and end-season stages, respectively. These estimated crop coefﬁcient values were higher than FAO K c values. Statistical analysis results revealed that the overall stagewise-derived average K c values were in line with FAO values, but different from the derived pan K c values, although found insigniﬁcant at a 5% signiﬁcance level. In addition, water productivity and agro-meteorological indices were derived to evaluate the cultivar performance in this experiment. Therefore, such a methodology may be used in the absence of weighing lysimeter-derived K c values. The derived regional K c values can be applied to improve irrigation scheduling under similar agro-climatic conditions.


Introduction
Water is one of the most critical natural resources for almost all living organisms existing on the blue planet, and is needed for socio-economic development activities [1][2][3]. Freshwater resources are dynamically fluctuating globally owing to global warming, climate change and unsustainable anthropogenic actions [4]. Freshwater scarcity is a major challenge in today's world, and it has been projected to worsen by 2050 when the world population may reach between 9.4 and 10.2 billion [1]. Worldwide, agriculture water use currently accounts for 70% of the total water resource availability, and is mostly used for irrigation [5]. India accounts and provides satisfactory results [52,53]. Time-averaged crop coefficient values for low land paddy in the Senegal river valley were 1.01, 1.31, and 1.12 during the crop development, midseason, and late-season stages, respectively, and actual crop evapotranspiration was 841.5 and 855.4 mm in 2014 and 2015, respectively [29]. Doorenbos and Pruitt [54] proposed K c values for various crops cultivated in diverse agro-climatic regions.
Furthermore, previous studies showed that there is no database of crop coefficients for different crops in diverse agro-climatic regions of the MIGP. Therefore, FAO K c values are commonly utilized for irrigation scheduling in different crops. FAO-56, on the other hand, is recommended for developing in situ stagewise K c values for various crops for precise irrigation scheduling at the regional scale [20,30]. Therefore, the present investigation was conducted with two specific objectives, i.e., (1) to estimate ET a and develop stagewise K c values of transplanted puddled rice using modified non-weighing paddy lysimeter, and (2), to compare ET a and K c values with indirect methods such as Penman-Monteith and pan evaporation to improve in situ water productivity. In addition, water productivity and agrometeorological indices were also computed for drought-tolerant rice cultivar, viz., Swarn Shreya. The developed knowledge of rice ET a and K c values may aid in improving water use efficiency in similar agro-climatic regions by reducing the water wastage attributed to deep percolation and runoff produced in the vicinity of over irrigation due to the inadequacy of efficient irrigation scheduling information. Daily weather data were collected from ICAR-RCER, an agro-meteorological observatory, and their variation throughout the cropping season was depicted in Table 1.   Daily weather data were collected from ICAR-RCER, an agro-meteorological observatory, and their variation throughout the cropping season was depicted in Table 1. The region receives an annual mean precipitation of about 1213 mm. Soil samples were also collected using a core sampler, and soil physical properties were determined using the standard methods presented in Table 2. Paddy seedlings were transplanted on 26 July and harvested on 25 October 2021. All the recommended fertilizer doses, irrigation, and agronomical practices were followed across the entire field and within lysimeters to ensure uniform growth and check the oasis effect on crop ET.

Estimation of Actual Crop Evapotranspiration Using Paddy Lysimeter
Actual evapotranspiration was determined using a modified paddy non-weighing lysimeter based on the water budget equation expressed in Equation (1). where ∑ In f low = I + P; ∑ Out f low = ET a + R + DP; ∆Storage = ∆S within the root zone; *I is irrigation; P is rainfall; DP is percolation beyond the root zone; R is runoff; and ∆S is the change in soil moisture storage in the root zone.
Both I and P were directly measured. To measure ET a , a paddy lysimeter with slight modifications (outlet at recommended bund height for this region) was used in which three drums (namely, A, B, and C) with a diameter of 50 cm and height of 128.75 cm were buried in a rice field with roughly a fourth of their height above ground level ( Figure 2). Drum A had a closed bottom, while Drum B and C had open bottoms. The outlet pipe was installed at bund height (25 cm) in container C for precise runoff estimation.

Estimation of Actual Crop Evapotranspiration Using Paddy Lysimeter
Actual evapotranspiration was determined using a modified paddy non-weighing lysimeter based on the water budget equation expressed in Equation (1).
where ∑ = + ; ∑ = + + ; ∆ = ΔS within the root zone; *I is irrigation; P is rainfall; DP is percolation beyond the root zone; R is runoff; and ΔS is the change in soil moisture storage in the root zone.
Both I and P were directly measured. To measure , a paddy lysimeter with slight modifications (outlet at recommended bund height for this region) was used in which three drums (namely, A, B, and C) with a diameter of 50 cm and height of 128.75 cm were buried in a rice field with roughly a fourth of their height above ground level ( Figure 2). Drum A had a closed bottom, while Drum B and C had open bottoms. The outlet pipe was installed at bund height (25 cm) in container C for precise runoff estimation. In this experiment, the modified paddy lysimeter was a low-cost volumetric-type set of three drums used to measure and soil water balance components. The major modifications include setting the outlet at the recommended bund height (25 cm) of these re- In this experiment, the modified paddy lysimeter was a low-cost volumetric-type set of three drums used to measure ET a and soil water balance components. The major modifications include setting the outlet at the recommended bund height (25 cm) of these regions to estimate runoff from the field. Twenty-seven-day-old seedlings of the rice cultivar Swarn Shreya were transplanted using a standard agronomic package of practices inside the drums. The entire area of outside drums was also transplanted with the same rice variety. The water level inside the drums was maintained at the same level as it was in the cropped field. The difference in water level over two successive days in drum A represents actual evapotranspiration. In contrast, percolation loss is defined as the daily difference in water levels between drums A and B. The plots were bunded (25 cm) and there was no overflow from drum C throughout the season, resulting in no runoff from the paddy field.

Estimation of Reference Evapotranspiration
Reference evapotranspiration (ET 0 ) from 2017 to 2021 was computed using daily weather data collected from a local agro-meteorological observatory using FAO Penman-Monteith and US Class A pan evaporimeter techniques.

FAO Penman-Monteith Equation
Daily observations of air temperature (maximum and minimum), relative humidity, daily solar hours, and wind speed were considered for determining ET 0 using FAO's CROPWAT 8.0 software [20,24]. For estimating reference evapotranspiration (ET 0 ), the Penman-Monteith equation was used as follows: where ET 0 = a reference evapotranspiration (mm day -1 ), u 2 = a wind speed at 2 m height (m s -1 ), R n = a net radiation at crop surface (MJ m -2 day -1 ), G = a soil heat flux density (MJ m -2 day -1 ) Tmean = a mean daily air temperature at 2 m height ( • C), γ = a psychrometric constant (kPa • C -1 ), ∆ = a slope of vapor pressure curve (kPa • C -1 ), and e s − e a = a saturation vapor pressure deficit (kPa).

Pan Evaporation Approach
The pan evaporation method is a physical approach that may be applied worldwide without requiring parameter adjustments. In this method, the Class A evaporation pan from the US Weather Bureau, having a circular cross-section with a diameter of 120.7 cm and a depth of 25 cm, is widely used to monitor daily evaporation losses. The relationship between reference evapotranspiration and pan evaporation is expressed in Equation (3) ET 0 = K P × E pan (3) where ET 0 = a reference evapotranspiration (mm day -1 ), K p = a pan coefficient (0.71), and E pan = a daily pan evaporation (mm day −1 ).

Estimation of Stagewise Crop Coefficient
The impacts of actual crop evapotranspiration and soil evaporation are integrated into the single crop coefficient approach. In order to estimate stage-K c , daily reference evapotranspiration was estimated using CROPWAT 8.0 software, and modified paddy lysimeter was used to estimate actual daily evapotranspiration by observing the daily water level change in these three drums. Then, the crop coefficient was determined as the ratio of actual evapotranspiration to reference evapotranspiration, as given by Equation (4) where ET a = an actual crop evapotranspiration mm −1 and ET 0 = a reference crop evapotranspiration (mm −1 ). Pan-crop coefficient was estimated by the ratio of actual evapotranspiration to reference evapotranspiration based on the pan evaporation approach.

Statistical Analysis
Various statistical performance parameters such as mean absolute percentage error (MAPE), root mean square error (RMSE) and coefficient of determination (R 2 ) were calculated for error analysis (Equations (5)-(7)).
where P i , O i and O are the observed, estimated, and mean estimated values, respectively. Then, the significant test between the measured and estimated K c values was performed using the non-parametric Mann-Whitney U test. The flowchart of the whole methodology was presented in Figure 3.
where Pi, Oi and O ̅ are the observed, estimated, and mean estimated values, respectively. Then, the significant test between the measured and estimated Kc values was performed using the non-parametric Mann-Whitney U test. The flowchart of the whole methodology was presented in Figure 3.

Water Productivity and Agro-Meteorological Indices
In this study, two different approaches were used to estimate water productivity: crop water productivity (CWP) and total water productivity (TWP). Crop water productivity was determined by dividing the marketable grain yield by the water used to meet crop evapotranspiration. Then CWP of the rice cultivar (Swarn Shreya) was estimated by the following equation: Irrigation water use efficiency (IWUE), or total water productivity(TWP), originated from drought resistance and tolerance ideas [61]. In crop production, IWUEis defined as the crop yield ratio to the water applied in the field [62].
Water used to produce yield m 3 ha −1 (9) The daily meteorological data were used to calculate agro-meteorological indices, viz., growing degree days (GDD), heliothermal units (HTU), photothermal units (PTU), heat use efficiency (HUE), and heliothermal use efficiency (HTUE) (Equations (10)- (14), respectively). Daily GDD was calculated using a base temperature of 10 • C. The daily mean temperature above the base temperature, expressed in degree days, was added to calculate cumulative growing degree days.
Growing degree days (GDD) is calculated as: where T max -daily maximum temperature ( • C), T min -daily minimum temperature ( • C), and T b -minimum threshold/base temperature ( • C).
Heliothermal units (HTU): the product of the growing degree day and the corresponding bright sunshine hours has been termed heliothermal units (HTU). HTU = ∑ GDD × actual bright sunshine hours (n) (11) Photothermal units (PTU): the product of the growing degree day and the corresponding maximum possible sunshine hours has been termed photothermal units (PTU).
N-maximum possible sunshine hours. Heat use efficiency (HUE): heat use efficiency is also represented by thermal time use efficiency (TTUE), which indicates the amount of dry matter produced per unit of growing degree days or thermal time, expressed in g/m 2 / • day. This was computed by using the following formula: Heliothermal use efficiency (HTUE): heliothermal use efficiency was calculated by dividing the total dry matter recorded on respective days by the accumulated heliothermal units, and is expressed as kg ha −1 • C Day h.

Variation in Reference Evapotranspiration, Water Balance Components and Actual Crop Evapotranspiration
Daily mean ET 0 -PM varied from 1.72 to 5.67 mm day −1 with a mean value of 3.39 (±0.95) mm day −1 , while the corresponding daily ET 0 -Pan ranged between 1.56 and 3.98 mm day −1 with an average value of 3.06 ± 0.52 mm day −1 throughout the crop season. However, stagewise ET 0 -PM revealed 3.42 ± 0.73 mm day −1 , 3.56 ± 0.94 mm day −1 , 3.40 ± 1.03 mm day −1 , and 2.77 ± 0.55 mm day −1 during the initial, crop development, mid-season and late-season stages, respectively. Daily mean ET 0 -Pan findings showed their values to be 2.99 ± 0.66 mm day −1 , 3.16 ± 0.47 mm day −1 , 3.02 ± 0.56 mm day −1 , and 2.97 ± 0.29 mm day −1 during the initial, crop development, mid-season and end/lateseason stages, respectively. In addition, the daily variations in ET 0 -PM and ET 0 -Pan throughout the last five years (2017-2021) of the crop period were shown in Figure 4. The results showed that their values ranged between 2.69 and 4.92 mm day −1 , with an average value of 3.59 ± 0.56 mm day −1 for ET 0 -PM and 2.03 to 3.78 mm day −1 with an average of 2.78 ± 0.33 mm day −1 for ET 0 -Pan.
The overall relationship between ET 0 -PM and ET 0 -Pan presented in Figure 5 showed a poor relationship, with R 2 = 0.27 (RMSE:0.28, MAPE:8.34%) in the MIGP region. Throughout the crop period, lower reference evapotranspiration had been reported compared to ET 0 -PM, suggesting the development of local revised pan coefficient values in these regions.
During the cropping season, ET a using a modified paddy lysimeter varied from 1.9 to 8.2 mm day −1 , with a mean value of 4.02 ± 1.35 mm day −1 , and reported a total of 369.4 mm, accounting for 49.77% of the total water applied after transplanting. On the other hand, the percolation beyond the crop root zone was 262.48 mm, or around 35.36% of the total water applied, and the rest (14.87%) was stored within the root zone of the crop ( Figure 6). A total of 111.09 mm evapotranspiration was estimated during the nursery. In addition to this, the stagewise-measured daily average ET a was 3.86 ± 0.95 mm day −1 , 4.51 ± 1.41 mm day −1 , 4.05 ± 1.28 mm day −1 , and 2.52 ± 0.51 mm day −1 during the initial, crop development, mid-season, and late-season stages, respectively (Table 3). be 2.99 ± 0.66 mm day −1 , 3.16 ± 0.47 mm day −1 , 3.02 ± 0.56 mm day −1 , and 2.97 ± 0.29 mm day −1 during the initial, crop development, mid-season and end/late-season stages, respectively. In addition, the daily variations in ET0-PM and ET0-Pan throughout the last five years (2017-2021) of the crop period were shown in Figure 4. The results showed that their values ranged between 2.69 and 4.92 mm day −1 , with an average value of 3.59 ± 0.56 mm day −1 for ET0-PM and 2.03 to 3.78 mm day −1 with an average of 2.78 ± 0.33 mm day −1 for ET0-Pan. The overall relationship between ET0-PM and ET0-Pan presented in Figure 5 showed a poor relationship, with R 2 = 0.27 (RMSE:0.28, MAPE:8.34%) in the MIGP region. Throughout the crop period, lower reference evapotranspiration had been reported compared to ET0-PM, suggesting the development of local revised pan coefficient values in these regions. During the cropping season, using a modified paddy lysimeter varied from 1.9 to 8.2 mm day −1 , with a mean value of 4.02 ± 1.35 mm day −1 , and reported a total of 369.4 mm, accounting for 49.77% of the total water applied after transplanting. On the other hand, the percolation beyond the crop root zone was 262.48 mm, or around 35.36% of the total water applied, and the rest (14.87%) was stored within the root zone of the crop (Figure 6). A total of 111.09 mm evapotranspiration was estimated during the nursery. In addition to this, the stagewise-measured daily average was 3.86 ± 0.95 mm day −1 , 4.51 ± 1.41 mm day −1 , 4.05 ± 1.28 mm day −1 , and 2.52 ± 0.51 mm day −1 during the initial, crop development, mid-season, and late-season stages, respectively (Table 3). During the cropping season, using a modified paddy lysimeter varied from 1.9 to 8.2 mm day −1 , with a mean value of 4.02 ± 1.35 mm day −1 , and reported a total of 369.4 mm, accounting for 49.77% of the total water applied after transplanting. On the other hand, the percolation beyond the crop root zone was 262.48 mm, or around 35.36% of the total water applied, and the rest (14.87%) was stored within the root zone of the crop (Figure 6). A total of 111.09 mm evapotranspiration was estimated during the nursery. In addition to this, the stagewise-measured daily average was 3.86 ± 0.95 mm day −1 , 4.51 ± 1.41 mm day −1 , 4.05 ± 1.28 mm day −1 , and 2.52 ± 0.51 mm day −1 during the initial, crop development, mid-season, and late-season stages, respectively (Table 3).   Note: All values are given in mean ± standard deviation (SD).

Comparation of Actual Evapotranspiration with FAO-PMand Pan Evaporation-Derived Evapotranpiration during Crop Period
During the crop growing season, FAO-PM derived evapotranspiration (ET c -PM) varied from 1.55 to 6.8 mm/day with a mean value of 3.84 ± 1.21 mm day −1 , while the corresponding daily pan-derived evapotranspiration (ET c -Pan) varied from 1.87 to 4.77 mm day −1 with a mean value of 3.45 ± 0.69 mm day −1 . Similarly, daily ET a loss during the crop growing season varied from 1.9 to 8.2 mm/day, with a mean value of 4.02 ± 1.35 mm day −1 (Table 3). Overall, the variations in stagewise-estimated crop evapotranspiration using Penman-Monteith, the pan evaporation approach, and measured daily mean ET a were shown in Figure 7.

Comparation of Actual Evapotranspiration with FAO-PMand Pan Evaporation-Derived Evapotranpiration during Crop Period
During the crop growing season, FAO-PM derived evapotranspiration (ETc-PM) v ied from 1.55 to 6.8 mm/day with a mean value of 3.84 ± 1.21 mm day −1 , while the co sponding daily pan-derived evapotranspiration (ETc-Pan) varied from 1.87 to 4.77 m day −1 with a mean value of 3.45 ± 0.69 mm day −1 . Similarly, daily loss during the c growing season varied from 1.9 to 8.2 mm/day, with a mean value of 4.02 ± 1.35 mm da (Table 3). Overall, the variations in stagewise-estimated crop evapotranspiration us Penman-Monteith, the pan evaporation approach, and measured daily mean w shown in Figure 7.

Comparison of Stagewise Actual K c vs. Single FAO-56 K c Values and Derived Pan-Crop Coefficient Values
At the regional scale, actual crop coefficient values were estimated to be 1.13 ± 0.13, 1.27 ± 0.2, 1.23 ± 0.16, and 0.93 ± 0.18 for the initial, crop development, mid-season and late-season stages, respectively. The derived crop coefficient values were 7.62%, 5.83%, 2.5%, and 3.33%, higher during the initial, crop development, mid-season and end/late season stages than the FAO K c values for the MIGP region, while the pan-crop coefficient values estimated were 1.33 ± 0.41, 1.42 ± 0.36, 1.38 ± 0.48, and 0.85 ± 0.16, i.e., 26.67%, 18.33%, and 15% higher and 5.5% lower for the initial, crop development, mid-season and late-season stages respectively in comparison to the FAO K c values. Overall, the stagewise variation in measured K c compared to the FAO-56 K c values and Pan-K c values depicted in Figure 10. However, the non-parametric Mann-Whitney U test was performed to compare the association between these stagewise-derived averaged actual K c values with FAO-56 K c values and pan-crop coefficient values. The results showed that, at different crop growth stages, the actual average K c values were in line with FAO values but slightly deviated from the pan-coefficient values, though were non-significant at the 5% level of significance ( Figure 11). Overall, the coefficients of determination (R 2 ) between the derived average stagewise K c vs. single FAO-56 K c values and derived pan-crop coefficient values were shown in Table 4.        However, based on the average ET 0 -PM and ET 0 -Pan estimated from 2017 to 2021 in reference to the actual measured crop evapotranspiration during 2021, the average actual K c values were calculated to be 1.08 ± 0.29, 1.14 ± 0.33, 1.18 ± 0.37, and 0.84 ± 0.09, while the corresponding average K c -Pan values were 1.42 ± 0.35, 1.57 ± 049, 1.52 ± 053, and 0.90 ± 1.67 during the initial, crop development, mid-season and end/late-season stages, respectively.

Influence of Weather Variables on Reference and Actual Evapotranspiration
From the correlation matrix involving weather variables and evapotranspiration, it was confirmed that ET 0 -PM was positively and significantly influenced by solar radiation, sunshine hours, Tmax, and wind speed and negatively affected by RH mean , whereas in the case of ET 0 -Pan, pan evaporation, Tmax and solar radiation were the dominant factors that affected it positively and RH mean had a negative effect. More or less similar trends were also seen in ET a as well, where Tmax, sunshine hours, and solar radiation were the main determinants ( Figure 12).

Water Productivity and Agro-Meteorological Indices
The grain yield was measured to be 0.548 kg/m 2 . Irrigation water use efficiency (IWUE) was estimated to be 0.74 kg/m 3 , whereas crop water productivity (CWP) was 1.14 kg/m 3 . Agrometeorological indices, viz., accumulated GDD ( • C days), HTU ( • C days h), PTU ( • C days day −1 ), and RTD during different phenophases, as well as HUE (kg ha −1 • C Day), HTUE (kg ha −1 • C Day h) and PTUE (kg ha −1 • C Day hrs) for grain and dry matter production under different phenophases, were calculated and presented in Tables 5 and 6. The results showed that the accumulated GDD, HTU, PTU and RTD values of the rice cultivar (Swarn Shreya) were estimated to be 2321 • C days, 8138 • C days hr, 25,995 • C days day −1 , and 26.3 • C, respectively, throughout the crop period (transplanting to harvesting). However, HUE, HTUE and PTUE were estimated to be 2.40 kg ha −1 • C day, 0.62 kg ha −1 • C day h, and 0.21 kg ha −1 • C day h, respectively, which will help in studying the crop response to ambient temperature at different phenological stages.
From the correlation matrix involving weather variables and evapotranspiration, it was confirmed that ET0-PM was positively and significantly influenced by solar radiation, sunshine hours, Tmax, and wind speed and negatively affected by RHmean, whereas in the case of ET0-Pan, pan evaporation, Tmax and solar radiation were the dominant factors that affected it positively and RHmean had a negative effect. More or less similar trends were also seen in as well, where Tmax, sunshine hours, and solar radiation were the main determinants (Figure 12). .

Water Productivity and Agro-Meteorological Indices
The grain yield was measured to be 0.548 kg/m 2 . Irrigation water use efficiency (IWUE) was estimated to be 0.74 kg/m 3 , whereas crop water productivity (CWP) was 1.14 kg/m 3 . Agrometeorological indices, viz., accumulated GDD (°C days), HTU (°C days h), PTU (°C days day −1 ), and RTD during different phenophases, as well as HUE (kg ha −1 °C Day), HTUE (kg ha −1 °C Day h) and PTUE (kg ha −1 °C Day hrs) for grain and dry matter production under different phenophases, were calculated and presented in Tables 5 and  6. The results showed that the accumulated GDD, HTU, PTU and RTD values of the rice cultivar (Swarn Shreya) were estimated to be 2321 °C days, 8138 °C days hr, 25,995 °C  Table 5. Accumulated GDD ( • C days), HTU ( • C days hr), PTU ( • C days day −1 ), and RTD ( • C) during different phenophases of rice cultivar.

Variation in Reference Evapotranspiration, Water Balance Components, and Actual Evapotranspiration
Overall, a poor relationship was found between the FAO-PM and pan evaporation approaches for estimating reference evapotranspiration, with R 2 = 0.27, RMSE = 0.28, and MAPE = 8.34% for the last five years (2017-2021). The reason behind the lower reference evapotranspiration in the case of the pan approach in comparison to ET 0 -PM is mainly due to a constant pan coefficient (0.71) that needs to be revised at the local scale to imitate the evaporation from the free water surface. This resulted in the underestimation of vegetation role as well as human errors. Similarly, Chatterjee et al. [30] conducted research in Cuttack, Odisha, and part of Eastern India, and reported a poor relationship between FAO-modified PM and Epan-derived ET 0 , with R 2 varying from 0.16 to 0.24, which is in close agreement with our findings. The decrease in Epan could be due to a phenomenon known as the "Pan Evaporation Paradox" [63,64]. However, previous studies also revealed that ET a can exceed Epan for short periods of time (e.g., at different stages of plant growth) when there is a lot of advected heat and a drop in the diurnal temperature range [63,[65][66][67][68][69]. Daily ET a used 369.4 mm of water, accounting for 49.77% of the total water applied after transplanting to the harvesting stage. Likewise, Dash et al. [70] reported that around half of the water applied in irrigated puddled rice fields is used in the form of evapotranspiration. However, daily ET a using a modified paddy lysimeter varied from 1.9 to 8.2 mm/day, with a mean value of 4.02 ± 1.35 mm day −1 . This variation is mainly due to various biophysical parameters such as vegetation type and density, growth phase, stomatal conductance, weather variations, local field management practices, etc. [34]. Several researchers have indicated that daily mean ET a rates of rice in the Asian continent ranged from 4 to 7 mm per day [30,71], which is in line with the ET a values obtained in our study. However, overall, daily measured ET a increased erratically with the increase in days after transplanting (DAT), reached peak value at the crop development to mid-season stages, mainly due to a larger leaf area index leading to maximum transpiration losses, and eventually decreased at the end-season stage due to leaf senescence [16,18].
ET a estimation using a modified lysimeter, however, had its own limitations. For example, because the drums were installed about 36 cm above the ground, the microclimate (Tmin slightly increased and Tmax slightly decreased by about 0.1-0.3 • C compared to the open field, while relative humidity varied by 1-2%) within the drum differed from the entire field and affected ET a estimation, but only up to 15-20 days after transplanting. After that, plant height increased and came outside the drum. Therefore, the microclimate did not have too much affected ET a estimation. However, the exposed rim of modified lysimeter also caused local advection, but did not significantly affect the ET a estimation. Similarly, research on the potential impacts of advection on flooded rice has not revealed any substantial advection implications [72]. The FAO-recommended Penman-Monteith method is considered to be the best method for ET 0 estimation due to the fact that it considers all the weather parameters which influence ET 0 [20].

Comparison of Crop Evapotranspiration and Stagewise Crop Coefficients
Based on the results, overall ET a measured using a modified paddy lysimeter showed a close association with FAO PM-derived evapotranspiration (R 2 : 0.63). However, the ET 0-Pan approach underestimated (R 2 : 0.24) in these regions. Maina et al. [73] also reported FAO-PM as a better approach than the pan approach due to its good approximation of trustworthy lysimeter readings. In these areas, there has been a documented underestimation of the pan method. The main cause is the constant pan coefficient, which should be updated for these regions [74][75][76]. Similarly, the bulk of the estimated pan coefficients, according to Sabziparvar et al. [77], were not statistically used in the pan ET 0 conversion technique. As a result, for each climate, a good model for calculating the pan coefficient is required [78]. However, the overall actual evapotranspiration estimated using the FAO-PM and pan approaches increased significantly as the phenological phases of growth progressed, peaked at mid-season, and subsequently declined at maturity during the crop growing season, similar to previous Indian and international studies [42,79,80]. Tyagi et al. [18] also reported similar results to those obtained in our study. These variations are mainly due to changes in ET a because of weather conditions, crop type and variety, crop growth stages, and local management practices in these areas [33,72]. Rice crop coefficient values for the initial, crop development, mid-season and late-season stages were calculated as 1.13 ± 0.13, 1.27 ± 0.2, 1.23 ± 0.16, and 0.93 ± 0.18, respectively, and results showed that the overall average crop coefficient values at different stages showed at par with the FAO K c values, but slightly deviated with pan-crop coefficient values. Statistically, these variations were insignificant (at a 5% level of significance). Higher crop coefficient values had been observed, except for end/late-season stage, mainly due to the stagnation of water during the crop season. Depending on the irrigation method, the K c in the initial stage, with the exception of permanent flooding, declined in comparison to the FAO 56 K c value [81]. Various findings also reported that stagewise crop coefficients for rice varied from 0.61 to 1.15, 0.97 to 1.44, 1.14 to 1.59, and 0.9 to 1.02, in the initial, crop development, mid-season, and late-season stages, respectively, using a lysimetric and other water balance approaches which were consistent with our findings [14,18,24,26,27,50,51]. Many distinct crop varieties and diverse water management practices, ranging from deep water to floating up to aerobic sprinkler settings, have been observed to cause a significant difference in the paddy crop coefficient [82]. In addition, higher pan-crop coefficient values were found owing to the constant and fixed pan coefficient value (0.71), resulting in significant ET 0 estimation errors which need to be revised in these locations, and this approach also underestimated the vegetation contribution [77,81,83]. Undertaking studies to revise the pan coefficient values locally considering the influence of weather parameters on pan evaporation in different climatic conditions is recommended.

Influence of Weather Variables on Reference and Actual Evapotranspiration
From the results, it was concluded that solar radiation, sunshine hours, Tmax, and wind speed positively affected evapotranspiration, whereas RHmean negatively affected evapotranspiration (ET 0 -PM, ET 0 -Pan and ET a ). Similar patterns of influence of these weather variables on evapotranspiration estimation were also reported in semi-arid IGPs [84]. Additionally, a significant factor influencing the weekly fluctuation in ET a was reported as net solar radiation variation (r ≥ 0.87) in IGPs by [7]. Chatterjee et al. [30] also reported a strong positive relationship between ET a and Tmean, as well as solar radiation in Eastern India. These variables mainly affected the sensible heat flux from soil to canopy and substantially affected their estimations during crop season.

Assessment of Water Productivity (WP) and Agro-Meteorological Indices
The results revealed that the IWUEand CWPof rice were 0.74 and 1.14 kg/m 3 , respectively. This figure fits within the range of the internationally recorded WP of rice (0.6-1.6 kg/m 3 ) [85]. Similarly, the crop water productivity ranged between 0.9 and 1.54 kg/m 3 , while total WP was varied from 0.68 to 0.98 kg/m 3 for different rice varieties [85]. Rice water productivity ranged from 0.6 to 1.5 kg/m 3 , which supports our findings according to Mboyerwa et al. [86]. These variations are mainly due to climate, soil type, soil structure, crop type, variety, duration, irrigation methods, fertilizer applications, etc. Agrometeorological indices showed that temperature affects crop growth and yield. These indices also varied due to different transplanting dates and varieties [84]. In the present study, various agrometeorological indices for the rice cultivar Swarna Shreya had been calculated and their variations showed the pattern of biomass accumulation at different phenological stages, which ultimately affect crop productivity. This type of study helps in understanding crop response to temperature variation for future yield forecasting in these regions under climate change scenarios.

Conclusions
Agriculture water management is aided by accurate ET quantification and crop coefficient computation during different crop growth stages. In this study, the stagewise average daily crop evapotranspiration rates (mm day −1 ) were estimated as 3.86 ± 0.95, 4.51 ± 1.41, 4.05 ± 1.28, and 2.52 ± 0.51, and their corresponding K c values were 1.13 ± 0.13, 1.27 ± 0.2, 1.23 ± 0.16, and 0.93 ± 0.18 for the initial, crop development, mid-season and late-season stages, respectively. The derived average crop coefficient values were reported 7.62%, 5.83%, 2.5%, and 3.33%, higher during the initial, crop development, mid-season, and end-season stages compared to FAO K c values. However, the average pan-crop coefficients estimated as 1.33 ± 0.41, 1.42 ± 0.36, 1.38 ± 0.48, and 0.85 ± 0.16 were 26.67%, 18.33%, and 15% higher and 5.5% lower than the initial, crop development, mid-season and late-season stages, respectively, in comparison to FAO K c values. Overall, the results showed good agreement (R 2 = 0.63) between ET a and FAO-PM-calculated ETc. Conversely, the pan evaporation approach underestimated the crop evapotranspiration compared to ET a , which suggests revising the pan coefficient values for the study region will achieve better accuracy. Thus, we recommend revising the pan coefficient values considering the influence of weather parameters. The IWUEand CWP for rice were estimated as 0.74 and 1.14 kg/m 3 , respectively. Similarly, agrometeorological indices were also estimated for Swarna Shreya to help in forecasting the occurrence of different phenophases in similar agro-climatic regions. Therefore, this type of methodology could also be used in other edaphoclimatic and data-scarce regions where lysimeters are unavailable to estimate actual evapotranspiration in paddy areas. In addition, the derived stagewise K c values at the regional scale may help in efficient irrigation planning for similar agro-climatic conditions. Overall, K c and ET a values may be beneficial to improve irrigation scheduling in water-intensive rice crop.