Assessment of Soybean Evapotranspiration and Controlled Water Stress Using Traditional and Converted Evapotranspirometers

: Evapotranspiration ( ETR ), reference evapotranspiration ( ET 0 ), and seasonal ETR totals were determined for soybean over two growing seasons, at Keszthely, Hungary, using traditionally operated and converted evapotranspirometers. The study aimed to document the plant–water response of two soybean varieties ( Sinara : Sin ; Sigalia : Sig ) which have di ﬀ erent water demands. Three water supply treatments were tested: unlimited ( WW ) watering, 50% of crop water requirement ( RO ), and rainfed ( P ). Reconstructed evapotranspirometers allowed crop water deprivation to be simulated under ﬁeld conditions. ETR sums were higher during the cooler 2017 than in the warmer 2018, calling attention to the importance of being informed about more detailed meteorological variables other than monthly (seasonal) means. In addition to variation in daily mean air temperatures ( T a ), maximum T a played a key role in determining ETR under naturally occurring extreme weather conditions in 2018. Irrespective of the variety, daily mean ETR was on average 65–75% greater than in the water-stress treatment. Unexpectedly, water stress-tolerant Sin used slightly more water than Sig , which was bred for standard weather conditions. Measured mean ETR was as much as 10% higher than derived ET 0 rates, causing crop coe ﬃ cient to exceed 1.0 during ﬂowering. Careful selection of the soybean variety when practicing water-saving management may lead to more e ﬃ cient variety improvement in a breeding program. It may also be important for soybean producers and farmers to adopt the best variety, aiming to decrease the use of irrigation water to increase seed yield.


Introduction
In most soybean (Glycine max (L.) Merr.)-growing areas in temperate climatic zones, including Hungary, there is a need to meet nutritional and food demands locally. The import of soybean seed is very expensive, and due to the tropical origin of the crop, soybean growth requires irrigation in some years. Soybean is the world's most important protein and oil source, and its use and consumption in food is recommended due to its high content of Ω6 and Ω3 fatty acids [1]. Soybean yield is strongly dependent on the genotype × environment interaction [2], but it is a crop that is resilient to saline environments, allowing for grain production under drought stress [3]. (b) denote the experimental field and observation garden, respectively. Three water treatments and measuring towers ( • ) were placed, as shown in a, as follows: X, rainfed plots; 1, 1-12 evapotranspirometer pots in the field. The evapotranspirometers' pots with different water and variety treatments are highlighted as follows: ■, pots of a traditional evapotranspirometer with Sinara (Sin WW); •, pots of water-stressed treatment with Sinara (Sin RO); ■, pots of a traditional evapotranspirometer with unlimited water in Sigalia (Sig WW); •, pots of water-stressed treatment with Sigalia (Sig RO). In (b) the following abbreviations were used: 2, measuring cellar; 3, US class "A" pan; 4, wind tower with a pyranometer; 5, QLC-50 automatic climate station; 6, evaporation pan with a surface area of 20 m 2 . The same experimental setup was used in both observation seasons (2017 and 2018). Several plant indicators (phenology, leaf and leaflet areas, chlorophyll index, yield, water use efficiency) were measured in all water and variety treatments (upper part (a): 1 and X). Daily evapotranspiration was measured in (a) 1, 1-12 evapotranspirometer pots. The crop coefficient was derived from these equipment pots. Evapotranspiration totals were also estimated from the field (a). All meteorological variables were derived from observational garden (b). (b) denote the experimental field and observation garden, respectively. Three water treatments and measuring towers ( • ) were placed, as shown in a, as follows: X, rainfed plots; 1, 1-12 evapotranspirometer pots in the field. The evapotranspirometers' pots with different water and variety treatments are highlighted as follows: ■, pots of a traditional evapotranspirometer with Sinara (Sin WW); •, pots of water-stressed treatment with Sinara (Sin RO); ■, pots of a traditional evapotranspirometer with unlimited water in Sigalia (Sig WW); •, pots of water-stressed treatment with Sigalia (Sig RO). In (b) the following abbreviations were used: 2, measuring cellar; 3, US class "A" pan; 4, wind tower with a pyranometer; 5, QLC-50 automatic climate station; 6, evaporation pan with a surface area of 20 m 2 . The same experimental setup was used in both observation seasons (2017 and 2018). Several plant indicators (phenology, leaf and leaflet areas, chlorophyll index, yield, water use efficiency) were measured in all water and variety treatments (upper part (a): 1 and X). Daily evapotranspiration was measured in (a) 1, 1-12 evapotranspirometer pots. The crop coefficient was derived from these equipment pots. Evapotranspiration totals were also estimated from the field (a). All meteorological variables were derived from observational garden (b).
, pots of water-stressed treatment with Sigalia (Sig RO). In (b) the following abbreviations were used: 2, measuring cellar; 3, US class "A" pan; 4, wind tower with a pyranometer; 5, QLC-50 automatic climate station; 6, evaporation pan with a surface area of 20 m 2 . The same experimental setup was used in both observation seasons (2017 and 2018). Several plant indicators (phenology, leaf and leaflet areas, chlorophyll index, yield, water use efficiency) were measured in all water and variety treatments (upper part (a): 1 and X). Daily evapotranspiration was measured in (a) 1, 1-12 evapotranspirometer pots. The crop coefficient was derived from these equipment pots. Evapotranspiration totals were also estimated from the field (a). All meteorological variables were derived from observational garden (b).

Weather, Crop and Soil Characteristics
A QLC-50 climate station (Vaisala, Helsinki, Finland) equipped with a CM-3 pyranometer (Kipp & Zonen Corp., Delft, the Netherlands) was operated at ARS. The combined sensor was placed at a standard height of 2 m above the ground level. Signals from T a , relative humidity (RH), wind speed (u), and global radiation (R s ) were collected every 2 sec, and 10 min averages were logged. The height of the anemometer was at 10.5 m above the ground level. These data were converted to 2 m when computing ET 0 .
Phenological developmental stages were continuously followed by using the Fehr and Caviness (1977) [27] scale in both growing seasons.
In weekly determination of leaf area, each green trifoliate leaf was placed on flat pieces of evenly colored red cardboard of known area, directly under a vertically mounted digital camera (Canon EOS 7D, Tokyo, Japan). Leaf area was determined photographically for five sample plants in each treatment. A histogram, based on the segmentation method, was applied in an image processing program (SGDIP 0.1, on-site designed) to count the area of region of interest [28]. The threshold point was set between the two peaks representing the background and the leaf area. Leaf-area index (LAI) was calculated as the ratio of leaf area divided by total ground area occupied by each plant. To study the distribution of leaflet size inside trifoliate leaves, LAI max was selected and evaluated when leaflet area was relatively invariable. This time period coincided with the peak of the growing season, from the end of July to the first 10 days of August. A high resolution photo of every leaflet from 10 plants was taken for each treatment. About 3000 leaflets per season were measured. At the same time, a SPAD-502 chlorophyll meter (Spectrum Technologies Inc., Plainfield, IL, USA) was applied to detect the "greenness" of individual leaflets. This device transmits light through a leaf at wavelengths of 650 and 940 nm. Any characteristic that alters the color of a crop (e.g., water shortage) may influence the value of SPAD readings. Ten replications were used for each treatment.
The soil type was a Ramann-type brown forest soil (clay loam; Haplic cambisol World Reference Base for Soil Resources 2015-WRB, 2006 [29]) with a mean bulk density of 1.15 Mg m −3 in the top 1 m of the profile and 273 mm m −1 of available water. The pots of the evapotranspirometer were also filled with the upper 0.9 m soil layer from the surroundings of the installed evapotranspirometers. Soil water was gravimetrically measured at 0.1 m depth intervals, down to 1 m, twice in each growing season, at the seeding stage in May and at full maturity at the beginning of September. Samples were weighed immediately after sampling, and oven-dried at 105 • C for 72 h until constant weight. Soil water was converted into available soil water (ASW), which is the amount of water in mm above the wilting point [30].

ETR and Derived Variables for The Field
In the field, to estimate total ETR, water balance was used. The only time when the crop was irrigated was after sowing on May 10th 2017. The amount of irrigation water was equal to 10 mm to obtain uniform emergence in the field. ETR totals in rainfed plots was derived as follows: where PR is precipitation during the growing season (mm), I is the amount of irrigation water (mm), C is the capillary rise of water from below the root zone (mm), R is the surface runoff (mm), D is deep percolation (mm), and ∆S is the cumulative change in soil water storage between planting and harvest (mm). ETR sum in the field estimation was based on the soil water balance method (∆S based on gravimetrically measured soil water at the beginning and end of the season, and PR). C was assumed to be zero, due to the 1.5 m aquitard sand layer that hinder the capillarity. The water table depth was at 8 m. D was also ignored because there was only one occasion when extreme PR occurred. Due to a flat field and previous monthly PR sums, the only exceptional PR value was in July 2018 after a dry period, Atmosphere 2020, 11, 830 6 of 27 which completed the missing soil water without deep percolation. Following the work of Gajić et al., any PR value below 100 mm were assumed to be zero run-off and run-on [14].
Thornthwaite-Mather type compensation evapotranspirometers were used to detect the daily water use of crops. The fixed pots are metal containers with a volume of 4 m 3 (surface area: 4 m 2 ; 1 m in depth). Three replications provided unlimited water supply to each variety. In water-stressed treatments, the number of replicates was the same. The pots were filled with a 0.1 m layer of gravel at the bottom and the remaining volume of pots was filled with soil from the upper layer of an adjacent field. The water level in the compensational vessel was controlled by a swimming ring magnet and a Reed-switch. The Reed-switch was connected to the water supply tap directly by a magnetic switch. A data logger (HYGACQ V1.3 type) was connected to log the amount of water use. The WHYGACQ computer program (HYGACQ, EWS Limited Partnership, Eger, Hungary) was used to upload data. Measured daily ETR values were expressed as a residual term from the water balance equation. By closing the water supplier tap every second day, 50% water withdrawal was ensured. More details are available in Anda et al. (2019a) [31]. In the RO treatment, the contribution of PR was also excluded by means of mobile rainout shelters (2.5 m long, 4.5 m wide, 2-2.5 m in height) covered with 200 µm thick UV transparent polyethylene film (MOL Petrolkémia Ltd., Budapest, Hungary). This mobile equipment was only installed during PR events to avoid increases in T a (extra heat stress of RO crops) under rain-free time periods.
Three water levels were applied in the current study: WW: use of Thornthwaite-Mather type compensation evapotranspirometer; RO: 50% water deprivation from the reproductive stage; P: rainfed plots.
Daily crop ET 0 rate (mm day −1 ) was computed by the widely-applied FAO-56 Penman-Monteith equation [32,33]: where G is soil heat flux density (MJ m −2 day −1 ), T a is mean daily air temperature at 2 m height ( • C), u is wind speed at 10.5 m height (m s −1 ), e s is saturation vapor pressure (kPa), e a is actual vapor pressure [kPa], ∆ is the slope of the vapor pressure curve (kPa • C −1 ) and 0.408 is a conversion factor from MJ m −2 day −1 to equivalent evaporation in mm day −1 . R n was estimated from R s , daily mean T a , mean daily e a , the site latitude and elevation based on Allen et al. (2005) [34]. A fixed value of 0.23 was applied for albedo. The on-site R n derivation was reported in more detail by Soós and Anda (2014) [35]. Soil heat flux density was assumed to be G = 0 on a daily basis. The value of 900 combines several variables from the original FAO-56 equation [32], while 0.34 is an empirical factor. The ratio of daily measured ETR (ET m ) and derived ET 0 rate resulted in the crop coefficient (K c ) for both soybean varieties: WUE for seed yield (y s ) was computed as the ratio of y s and seasonal total ETR [36] expressed as kg m −3 : At maturity, 2 × 2 m subplots were harvested in the highlighted areas of the field to determine the y s of crops. Five replicates per variety were included ( Figure 1). The harvested pods were oven dried at 65 • C for at least 48 h (until constant weight) to obtain pod dry weight. After threshing pods by hand, seed weight was obtained and adjusted to 13% moisture. The process was the same in each evapotranspirometer pot.

Statistical Analyses
The experimental design was a complete block with three replications due to the fixed nature of the evapotranspirometer.
The effect of the treatment, variety and season on LAI as well as SPAD was analyzed by three-way analysis of variance (ANOVA). The starting model was a full-factorial model with all main effects and all two-and three-way interactions. The non-significant interaction with the highest p-value was removed from the model iteratively. Decisions are detailed in the Results and Discussion section.
Differences between experimental groups were assessed using a Student's t-test or one-way ANOVA followed by Tukey's (honestly significant difference) HSD post hoc test (p < 0.05).
A two-tailed one-sample t-test was applied to compare the proportion of LAI, SPAD, K c and WUE. To facilitate the presentation of results, a 95% confidence interval was calculated.
Boxplots were applied to show differences between the LAIs of soybean varieties, water treatments and seasons.
To illustrate the impact of meteorological elements (T a ; RH; solar radiation (R s ); u; VPD) on ETR, Pearson's correlation analysis was applied. Multiple stepwise regression analysis was carried out to obtain the combined effect of different meteorological variables on ETR.
All statistical computations were implemented using SPSS software, version 17.0 [37] and R-statistics software [38].

Weather Conditions of the Studied Growing Seasons
The location of the study, Keszthely, which is positioned in the Carpathian Basin, is expected to have strongly variable inter-and intra-annual climate conditions (mild, continental climate with long warm and dry summers). The long-term (1971-2000) growing season's average T a is 16.9 • C ( Table 1). The climate of Keszthely is typified by highly variable and irregular PR events with long-term seasonal total of 384.3 mm. The growing season's monthly mean T a followed a similar pattern in both 2017 and 2018: they were unusually warm, with a seasonal average T a of 1.3 (p = 0.075) and 2.4 • C (p = 0.006) that exceeded the long-term mean. The only exception in 2017 was the cool (and dry) April that caused a slight delay in soybean sowing. The 2017 growing season was slightly drier than the long-term mean, receiving 10.1% less (p = 0.947) PR. The amount of PR was more frequent and of larger magnitude in 2018 than during the 2017 growing season. The sum of PR in the 2018 growing season was 21.6% higher (p = 0.093) than the long-term average. The end of each seasons was typified by considerable monthly PR, in September, which had no effect on soybean development. However, after disregarding the extremely high monthly sum of PR in September, and considering the soybean's remaining growing season, 2017 received 44.9% less (p = 0.006) while 2018 received 6.4% (p = 0.046) more PR than the climate norm.

Phenological Stages
Independent of the fact that seeding in 2018 was earlier than in 2017 by two weeks, soybean plant growth, development, and phenological stages were similar in both seasons under different irrigation water treatments (p = 0.211-0.441). The length of soybean development ranged from 114 to 120 days across treatments in 2017 (115-121 days in 2018). The only difference between the length of phenological stages was the one and a half times longer flowering (R1-R2) in 2018, probably due to the wetter weather conditions during June (18.6% lower monthly mean VPD and 40 mm higher monthly PR sum), in comparison to the previous season's weather during June ( Figure 2).
Water deprivation at the reproductive stage slightly shortened the length of the vegetative cycle in both varieties. In both Sin RO and Sig RO, seed maturity occurred in 4 (p < 0.03) and 6 days (p < 0.026) earlier in 2017 and 2018, respectively, compared to their WW treatments. Irrespective of the season studied, the length of the vegetative cycle of rainfed (118 days) crops was intermediate between the WW (121 days) and RO (116 days) values over the two-season study. The vegetative cycle of a rainfed soybean cultivar, Don Mario, was completed seven days earlier than its irrigated treatment in Uruguay [38].   VE-emergence; V3-third node; V5-fifth node; R1-start of bloom; R2-full bloom; R3-start of pod; R4-full pod; R5-start of seed; R6-full seed; R7-start of maturity; R8-full maturity. X-axis: day of year (DAS). WW-evapotranspirometer with unlimited water supply; RO-50% water deprivation in evapotranspirometer; P-rainfed. V-vegetative; R-reproductive. DAS-days after sowing.
Water deprivation at the reproductive stage slightly shortened the length of the vegetative cycle in both varieties. In both Sin RO and Sig RO, seed maturity occurred in 4 (p < 0.03) and 6 days (p < 0.026) earlier in 2017 and 2018, respectively, compared to their WW treatments. Irrespective of the season studied, the length of the vegetative cycle of rainfed (118 days) crops was intermediate between the WW (121 days) and RO (116 days) values over the two-season study. The vegetative cycle of a rainfed soybean cultivar, Don Mario, was completed seven days earlier than its irrigated treatment in Uruguay [39].

Leaf (Leaflet) Area of Soybean
Seasonal LAI means ranged from 3.3 ± 0.64 (Sig RO in 2017) to 5.3 ± 0.51 (Sin WW in 2017) in the two-season study. There was no significant impact of season on mean LAI in WW. Decreases in the seasonal LAI means were 44.4% (p < 0.001) and 44.7% (p < 0.001) in Sin RO and Sig RO, respectively, during 2017, in comparison to average LAI of crops grown in WW. In the next growing season, water deprivation reduced LAI by 24.2% (p = 0.002) and 25.3% (p = 0.003) in Sin RO and Sig RO, respectively, compared to WW. Irrespective of the season, mean LAI of P was intermediate to the WW and RO values ( Figure 3).  Independent of the variety and season, the LAI of P was always lower than that of WW. Furthermore, the LAI of RO was even lower than that of P for each season-variety combination. The difference between the RO treatment and the associated WW as well as P treatment was different for each season-variety combination. As a rule, LAI values of P were intermediate between WW and RO values. In 2018, due to higher PR, the LAI of P treatment was closer to that of WW treatment, while in drier 2017 it was closer to that of RO treatment. The lower LAI in P caused by drier weather conditions in 2017 may have occurred due to greater investment in root structure to make more effective use of soil moisture at deeper soil layers [39].
Higher seasonal variation in different LAI characteristics including median and quartiles Independent of the variety and season, the LAI of P was always lower than that of WW. Furthermore, the LAI of RO was even lower than that of P for each season-variety combination. The difference between the RO treatment and the associated WW as well as P treatment was different for each season-variety combination. As a rule, LAI values of P were intermediate between WW and RO values. In 2018, due to higher PR, the LAI of P treatment was closer to that of WW treatment, while in drier Atmosphere 2020, 11, 830 10 of 27 2017 it was closer to that of RO treatment. The lower LAI in P caused by drier weather conditions in 2017 may have occurred due to greater investment in root structure to make more effective use of soil moisture at deeper soil layers [40].
Higher seasonal variation in different LAI characteristics including median and quartiles occurred in water deprived RO and P crops. The lower LAI of RO confirmed that reduced compensation of water shortage in soybean growth might be expected under water stress during the reproductive stage [8]. The two-season mean LAI of Sin RO ranging from 3. The statistical analysis of LAI was carried out by three-way ANOVA. The three factors were season (S), variety (V) and water supply (W). In the full-factorial model, the three-way interaction was not significant (p = 0.75) so it was excluded from the model. In the model containing all two-way interactions, the V × W interaction had the largest p-value (p = 0.51), so it was removed from the model. Similarly, the V × S interaction had the largest p-value (p = 0.49) so it was also removed from the model. The main effect of V was slightly significant (p = 0.04), so it was also excluded. The final model included the main effect of W and S and their interaction ( Table 2). The interaction is illustrated in Figure 4. Analysis on the distribution of the middle and two outer leaflets in different trifoliate heights was carried out during the two experimental seasons ( Figure 5). A previous one-season result showed a significant difference in the size of the middle and two outer leaflets [44], irrespective of variety. The leaflet areas of the two varieties were pooled.
Probably due to the earlier seeding in 2018, the number of trifoliate layers increased with five more levels than in 2017. In addition, the size and level of trifoliates with the highest leaflet area also varied. Trifoliate level with maximum leaf area in the middle (outer) leaflet was observed in the 6th and 10th levels in 2017 and in the 7th and 10th levels in 2018. In 2017, increases in leaflet size were 35.0% (p < 0.001) and 41.4% (p < 0.001) in the middle and outer leaflets, respectively, relative to the 2018 values. Summarizing the results for leaflets, we conclude that different weather conditions in the two seasons caused variation to both the number and leaf areas of trifoliates. In 2017, fewer leaf layers with higher leaf areas, relative to 2018, were observed. The interaction is illustrated in Figure 4. Analysis on the distribution of the middle and two outer leaflets in different trifoliate heights was carried out during the two experimental seasons ( Figure 5). A previous one-season result showed a significant difference in the size of the middle and two outer leaflets [43], irrespective of variety. The leaflet areas of the two varieties were pooled. Probably due to the earlier seeding in 2018, the number of trifoliate layers increased with five more levels than in 2017. In addition, the size and level of trifoliates with the highest leaflet area also varied. Trifoliate level with maximum leaf area in the middle (outer) leaflet was observed in the 6th and 10th levels in 2017 and in the 7th and 10th levels in 2018. In 2017, increases in leaflet size were 35.0% (p < 0.001) and 41.4% (p < 0.001) in the middle and outer leaflets, respectively, relative to the 2018 values. Summarizing the results for leaflets, we conclude that different weather conditions in the two seasons caused variation to both the number and leaf areas of trifoliates. In 2017, fewer leaf layers with higher leaf areas, relative to 2018, were observed.

Chlorophyll Content of the Leaves (SPAD)
SPAD values mirror the "greenness" of leaves and provide an index of the actual chlorophyll content. Increased SPAD may contribute to impaired photosynthetic activity leading to better crop performance and crop yield. The index was documented by Bajgain et al. (2015) [44] to be a determinant factor of soybean photosynthesis as well. Smaller chlorophyll content per unit leaf surface likely contributes to less absorbed radiation that may potentially be converted to crop dry matter. Photophosphorylation controlled by photosystems is triggered by the density of chlorophyll units in the leaf, and is related to SPAD values [45]. SPAD and water supplies were inversely related in both growing seasons. The highest SPAD values ranging from 42.7 ± 3.43 (Sig P; 2017) to 46.3 ± 2.14 less absorbed radiation that may potentially be converted to crop dry matter. Photophosphorylation controlled by photosystems is triggered by the density of chlorophyll units in the leaf, and is related to SPAD values [46]. SPAD and water supplies were inversely related in both growing seasons. The highest SPAD values ranging from 42.7 ± 3.43 (Sig P; 2017) to 46.3 ± 2.14 (Sin P; 2018) were obtained in rainfed crops. Decreases by 33.4% (p < 0.000), 33.5% (p < 0.001), 34.0% (p < 0.001) and 37.1% (p < 0.001) were observed in Sin WW (2017) The effect of treatment, variety and season on SPAD was analysed by three-way ANOVA. The starting model was the full-factorial model with all main effects and all two-and three-way interactions. The non-significant three-way interaction was removed from the model, then the new model was recalculated. This process was iterated until all parts of the model were significant. In each step, the interaction with the highest p-value was removed. The final model included the main effects of water and season as well as their interaction. The pairwise comparison of the six water-season combinations were compared with Tukey's HSD test ( Figure 6).  There were no significant differences in SPAD values between the two seasons in P (Sin: 4.8% (p = 0.0502); Sig: 6.7% (p = 0.067)). Higher SPAD values in rainfed soybean in the present study agreed well with the results of Ergo et al. (2018) [9] for rainfed soybean, variety SPS4 × RR, at Manfredi in Argentina, between 2012 and 2013.

ETR in Soybean
The seasonal patterns of daily mean ETR were similar to each treatment of the studied seasons, however, the size of water losses varied among different (water/season/variety) treatments. Temporal patterns in ETR were the same as in LAI development. An increase in water demand from emergence to the R5 phenological stage occurred. This was followed by a sharp decrease in crop senescence until the R7-R8 stages. For all treatments, lower ETR occurred in the vegetative stage, ranging from 0.1 to about 2.5 mm day −1 . Statistical characteristics in variation of daily mean ETR are summarized in Figure 7. The top daily water demands between 6.6 (Sin RO) and 10.9 (Sig WW) mm during July  There were no significant differences in SPAD values between the two seasons in P (Sin: 4.8% (p = 0.0502); Sig: 6.7% (p = 0.067)). Higher SPAD values in rainfed soybean in the present study agreed well with the results of Ergo et al. (2018) [9] for rainfed soybean, variety SPS4 × RR, at Manfredi in Argentina, between 2012 and 2013.

ETR in Soybean
The seasonal patterns of daily mean ETR were similar to each treatment of the studied seasons, however, the size of water losses varied among different (water/season/variety) treatments. Temporal patterns in ETR were the same as in LAI development. An increase in water demand from emergence to the R5 phenological stage occurred. This was followed by a sharp decrease in crop senescence until the R7-R8 stages. For all treatments, lower ETR occurred in the vegetative stage, ranging from 0.1 to about 2.5 mm day −1 . Statistical characteristics in variation of daily mean ETR are summarized In 2017, unlimited watering significantly increased daily mean ETR by 77.3% (p < 0.001) and 76.7% (p < 0.001) in Sin and Sig, respectively, in comparison to daily ETR of the RO crop (Figure 7). In the next experimental season, 69.4% (Sin: p < 0.001) and 59.4% (Sig: p < 0.001) increments were recorded in WW treatments relative to water-deprived (RO) treatments. Daily mean ETR of RO were similar to each other: 2.3 ± 2.28, 2.2 ± 2.24, 2.  In 2017, unlimited watering significantly increased daily mean ETR by 77.3% (p < 0.001) and 76.7% (p < 0.001) in Sin and Sig, respectively, in comparison to daily ETR of the RO crop (Figure 7). In the next experimental season, 69.4% (Sin: p < 0.001) and 59.4% (Sig: p < 0.001) increments were recorded in WW treatments relative to water-deprived (RO) treatments. Daily mean ETR of RO were similar to each other: 2.3 ± 2.28, 2.2 ± 2.24, 2.2 ± 1.63 and 2.2 ± 1.68 mm day −1 for Sin RO (2017), Sig RO (2017), Sin RO (2018) and Sig RO (2018), respectively (p = 0.524-0.750). Many outliers in RO rooted after stressed crops were watered every second day, and ETR of RO diverged considerably from R1 when water was withheld. Figure 7 were restricted to RO treatments during 2017 because daily ETR means in 2017 were more concentrated around small values while they were evenly distributed in the observed range in 2018. The lower the box height (interquartile range in 2017), the more easily the outliers evolved. The maximum ETR of RO crops ranged from 8.7 (Sin) to 10.1 mm day −1 (Sig) in 2017, while the more moderate maximum ETR ranged from 6.6 (Sin) to 7.0 mm day −1 (Sig) in 2018. The limited emergence of outliers in 2018 could be explained by analyzing the distribution of daily ETR in the two investigated seasons. The violin plots, which contain the density curves, provide additional information about the data distributions (Figure 8). With the exception of Sig WW in 2018, the shape of WW curves with highest LAI values were quite balanced, showing an evenly distributed daily ETR. This aspect was highlighted in the conclusion of Suyker and Verma (2008) [23], who explained that beyond threshold LAI (LAI > 3-4), after full-canopy closure, greater leaf areas would not significantly affect ETR. The long-necked thin "violins" of RO crops showed that the majority of ETR values were relatively low ("paunchy" violins). Surprisingly, the warmer season of 2018 produced shorter violins than 2017 curves, which could be entirely attributed to the variability in weather conditions of the two investigated seasons.

Outliers in
The seasons in the current study were considered to be warm relative to the climatic norm. However, changing patterns of weather perturbations, such as variations in Tamax and its distribution, might also impact the development of soybean. Despite 1 °C cooler weather in the 2017 season, daily mean ETR were 18.7% (p < 0.001) and 20.2% (p < 0.001) higher in Sin WW and Sig WW, respectively, when compared to their daily mean ETR over 2018. Despite this, there was no significant growth in the seasonal mean LAI of WW crops during 2018. Ta is a known driving factor of a crop's water loss. The growing seasons-mainly summers-of 2017 and 2018 were characterized by either considerable variation in Ta (in relation to 30-year averages), or ranges of Ta including distributions in Tamax ( Figure  9) that varied substantially between the two summers (June-August). With the exception of Sig WW in 2018, the shape of WW curves with highest LAI values were quite balanced, showing an evenly distributed daily ETR. This aspect was highlighted in the conclusion of Suyker and Verma (2008) [23], who explained that beyond threshold LAI (LAI > 3-4), after full-canopy closure, greater leaf areas would not significantly affect ETR. The long-necked thin "violins" of RO crops showed that the majority of ETR values were relatively low ("paunchy" violins). Surprisingly, the warmer season of 2018 produced shorter violins than 2017 curves, which could be entirely attributed to the variability in weather conditions of the two investigated seasons. The seasons in the current study were considered to be warm relative to the climatic norm. However, changing patterns of weather perturbations, such as variations in T amax and its distribution, might also impact the development of soybean. Despite 1 • C cooler weather in the 2017 season, daily mean ETR were 18.7% (p < 0.001) and 20.2% (p < 0.001) higher in Sin WW and Sig WW, respectively, when compared to their daily mean ETR over 2018. Despite this, there was no significant growth in the seasonal mean LAI of WW crops during 2018. T a is a known driving factor of a crop's water loss. The growing seasons-mainly summers-of 2017 and 2018 were characterized by either considerable variation in T a (in relation to 30-year averages), or ranges of T a including distributions in T amax (Figure 9) that varied substantially between the two summers (June-August). Atmosphere 2020, 11, x FOR PEER REVIEW 15 of 27 Extreme Tamax values in 2017 simultaneously impacted heat and drought stresses during soybean development. Dornbos and Mullen (1991) [49] documented that the simultaneous occurrence of heat and water stresses additively reduced yield in soybean that equaled the sum of individual stresses. However, this phenomenon was not confirmed by Ergo et al. (2018) [9], who identified water stress to be the dominant factor beside heat stress, motivating the greatest specific response to soybean variety SPS4 × RR. The experimental design used in our study would not enable us to separate the influences of the two stressors, but it is assumed that both might have occurred, mainly during 2017.
Previous investigations on soybean reported, that PR and flowering mainly are very sensitive to air moisture [50]. Two moisture characteristics were included in the analysis, RH and VPD. Both of them varied considerably in the growing seasons. In addition to an increased amount of seasonal PR in 2018, the number of days (101 days) in which daily mean RH exceeded 70% was about twice the number of similarly "humid" days between May and September in 2017 (55 days). This weather perturbation might also affect ETR due to the tropical origin of soybean, therefore higher RH is likely to induce more favorable growing conditions, mainly during flowering [9]. Daily mean RH under 70% was often associated with low temporary RH values of under 30% measured at around solar noon that could damage soybean development, even with proper soil water (i.e., plants grown in evapotranspirometers). Wei et al. (2015) [51] applied average minimum RH > 45% as the limit in a soybean (variety Zhonghuang No. 13) transpiration model experiment, and this required weather adjustment to model inputs in Daxing, China (39°37′ N), between 2008 and 2010. VPD, the "driving force" of ETR, is also a critical weather variable that controls crop transpiration via stomatal regulation. The second likely reason for the decline in daily average water loss in 2018 might be the  [49] found soybean to be very susceptible to extreme heat around the reproductive stage, when hot weather reduced crop growth and flowering. The number of days with daily T amax over 30 • C was higher in 2017 (38 days) than in 2018 (28 days) even though the mode of T amax distribution was below the threshold of 30 • C in both growing seasons.
Extreme T amax values in 2017 simultaneously impacted heat and drought stresses during soybean development. Dornbos and Mullen (1991) [50] documented that the simultaneous occurrence of heat and water stresses additively reduced yield in soybean that equaled the sum of individual stresses. However, this phenomenon was not confirmed by Ergo et al. (2018) [9], who identified water stress to be the dominant factor beside heat stress, motivating the greatest specific response to soybean variety SPS4 × RR. The experimental design used in our study would not enable us to separate the influences of the two stressors, but it is assumed that both might have occurred, mainly during 2017.
Previous investigations on soybean reported, that PR and flowering mainly are very sensitive to air moisture [31]. Two moisture characteristics were included in the analysis, RH and VPD. Both of them varied considerably in the growing seasons. In addition to an increased amount of seasonal PR in 2018, the number of days (101 days) in which daily mean RH exceeded 70% was about twice the number of similarly "humid" days between May and September in 2017 (55 days). This weather perturbation might also affect ETR due to the tropical origin of soybean, therefore higher RH is likely to induce more favorable growing conditions, mainly during flowering [9]. Daily mean RH under 70% was often associated with low temporary RH values of under 30% measured at around solar noon that could damage soybean development, even with proper soil water (i.e., plants grown in evapotranspirometers). Wei et al. (2015) [51] applied average minimum RH > 45% as the limit in a soybean (variety Zhonghuang No. 13) transpiration model experiment, and this required weather adjustment to model inputs in Daxing, China (39 • 37 N), between 2008 and 2010. VPD, the "driving force" of ETR, is also a critical weather variable that controls crop transpiration via stomatal regulation. The second likely reason for the decline in daily average water loss in 2018 might be the 20.3% (p = 0.023) decline in the VPD of the 2018 summer (June-August), which was a wet growing season. One trend was salient: lower VPD was accompanied by a loss in ETR intensity (also see Figure 7).
Seasonal average daily ET 0 of 5.3 ± 1.51 mm day −1 and 5.0 ± 1.45 mm day −1 were recorded in 2017 and 2018, respectively. A seasonal mean ET 0 of 5.2 (41 • N) and 4.9 mm day −1 (33 • 42 N) reported by Payero and Irmak (2013) [17] and Anapalli et al. (2018) [13], respectively were very similar to the ET 0 results in our study. The insignificant increment (7.7%; p = 0.067) in daily mean ET 0 in 2017 was probably due to the increase in mean VPD (from 1.28 to 2.03 kPa) in summer compared to the summer of 2018. Seasonal variation in measured daily mean ETR and daily average ET 0 rates followed the same trend observed by others, including Payero et al. (2005) [52], for soybean.
In 2017, increases in cumulative ETR by 79.3, 81.3 and 23.8 mm were obtained in Sin WW, Sig WW and Sig RO, respectively, relative to 2018; moreover, Sin RO had 9.7 mm lower crop water than the 2018 season. The average evapotranspiration totals of 655.2 mm in Sin WW and 616.7 mm in Sig WW in our two-season study were close to 633 mm reported by Payero et al. (2005) [52] for North Platte, Nebraska (41 • N) in an irrigated soybean, variety Renze 2600 RR, in 2003. It was also similar to the observation of Doorenbos and Kassam (1979) [53], whose cumulative ETR estimations in soybean ranged from 450 to 700 mm, which were controlled by the weather and cultivation conditions, the varieties used, the length of the vegetation cycle, and the observation time period.
Soybean water demand of 431-585 mm reported by Zhang et al. (2018) [16] for full irrigation conditions in Noxubee, Mississippi (31 • N) agreed well with the WW results in the current study. Luo et al. (2018) [15] reported an ETR total of 325 mm in rainfed soybean, at Iowa (41 • 55 N), between 2013 and 2014, which was consistent with the cumulative ETR of the RO crop ranging from 291 to 317 mm observed in the present study. Using soybean variety Zhonghuang, Wei et al. (2015) [51] reported seasonal ETR total ranging from 283 to 347 mm in Daxing, China (39 •  There were deviations in measured and derived daily ET 0 rates in both growing seasons. Figure 10a,b show that the computed ET 0 rates were slightly higher than the measured rates. Due to late sowing in 2017 (9th May), due to a cool April, measured ETR was considerably lower than ET 0 during the early growth cycle (highlighted by an ellipse in Figure 10a). A similar reduction in ETR as a result of delayed sowing was reported by Suyker and Verma (2008) [23] in soybean variety Asgrow at Mead, Nebraska (41 • N). Small, recently emerged crops with a LAI below 1 could not satisfy the increased atmospheric water demand in May, causing a wider variation in the relationship between seasonal measured ETR and derived ET 0 ( Figure 10). This deviation between measured ETR and

Crop Coefficients, Kc
The Kc values showed day-to-day variations in both experimental periods. The variability in seasonal Kc also showed apparent seasonal dynamics, and peaked together with seasonal ETR and

Crop Coefficients, K c
The K c values showed day-to-day variations in both experimental periods. The variability in seasonal K c also showed apparent seasonal dynamics, and peaked together with seasonal ETR and LAI in the most active period of soybean development (R5). Daily 14 and 0.48 in the initial, mid-season and late season, respectively, reported by Candogan et al. (2013) [55] at Bursa, Turkey (40 • N) in a two-season study (2005)(2006) were somewhat lower than the K c results of our experiment. A brief summary of previously observed K c values is presented in Table 3.   The warm spring in 2018 (end of April to the beginning of May) increased initial Kc values slightly. Irrespective of the season, the measured daily Kc values were characterized by rapid increases from early June to about mid-August, followed by a steady decrease due to reduced LAI, solar radiation and ETR towards September. This pattern is evident for variation in daily Kc, although the coefficients varied in each season (Table 4). The warm spring in 2018 (end of April to the beginning of May) increased initial K c values slightly. Irrespective of the season, the measured daily K c values were characterized by rapid increases from early June to about mid-August, followed by a steady decrease due to reduced LAI, solar radiation and ETR towards September. This pattern is evident for variation in daily K c , although the coefficients varied in each season (Table 4). Table 4. The monthly mean crop coefficient (K c ) for three developmental stages (V: vegetative growth; R1-R6: from the start of bloom to full seed; R7-R8: from the start of maturity to full maturity) based on Fehr and Caviness (1977) [27]. Abbreviations: Sinara: Sin; Sigalia: Sig; WW: no limited water supply. Slightly greater K c values obtained in 2017 were probably due to increased maximum T a that was also reflected in the higher ETR obtained in that season. Small K c values in the vegetative stage and throughout maturity were lower than 1, with one exception (Sin WW) in 2018. In Sig WW, K c values below 1 were characteristic of most of the growing season. There were no significant differences in K c between the two varieties, irrespective of the season.
Assuming that there may be large difference in geographic/weather conditions and agronomic management practices between observation sites, K c values cited in the literature might differ significantly. K c developed with an evapotranspirometer (unlimited watering) for two soybean varieties were higher than the value of 0.7 found by Karam et al. (2005) [8] in a lysimeter study, in Lebanon (33 • N), although their seasonal mean ETR of 5.5 mm day −1 was close to our average ETR of 5.3 and 5.0 mm day −1 in Sin WW and Sig WW, respectively, during 2018. However, soybean daily K c values between 0.27 and 1.47 during the 2007 and 2008 growing seasons noted by Irmak et al. (2013) [61] in south-central Nebraska (41 • 1 N) were consistent with the K c results of our study.
Given that K c values might be significantly affected by local environmental conditions and agronomic practices, our study is primarily appropriate for on-site use in the Carpathian Basin. One of the reasonable factors that might induce deviations in actual K c is the size of LAI, due to the modified fraction of soil cover (canopy opening), by impacting ETR [63]. The K c values were slightly influenced by LAI sizes in different treatments in the present study: higher LAI had higher K c caused by a higher ETR.
3.5. The Impact of Weather Variables on Daily Mean ETR ETR is a multivariate process affected by various meteorological factors [64]. In the current study, the influence of weather characteristics (T a , RH, u, R s , VPD and ET 0 ) on measured daily mean ETR of two soybean varieties (Sin and Sig) were investigated under unlimited (WW) watering at Keszthely ARS, during the 2017 and 2018 growing seasons (Table 5). Due to limiting water supply every second day, ETR of RO treatments and PR were excluded from the analysis. Irrespective of the season, the correlation coefficients (r) hardly fluctuated between the two varieties studied. This may be due to the low number of varieties studied. Moreover, differences in the relationship between ETR and weather variables, including VPD, were not confirmed. Table 5. Correlation coefficient (r) between daily measured evapotranspiration rate (ETR) and weather variables including mean air temperature (T a ), relative humidity (RH), wind speed (u), net radiation (R s ), and vapor pressure deficit (VPD) for the two soybean varieties, Sinara (Sin) and Sigalia (Sig). Crops with unlimited watering (WW) were included in the study. ET 0 represents the Penman-Monteith reference evapotranspiration.

Seasons
ET 0 (mm) Similar to the results of Payero and Irmak (2013) [17], with one exception of Sig WW in 2017, the most dominant variable with the highest r of 0.92 was VPD. There was one exception, Sig (2017), which possessed only a slightly higher r (0.01) for T a (0.74) than for VPD (0.73). In this case, the variation in r between T a and VPD was so small that it was not worth considering separately. VPD integrated the impact of both T a and RH, explaining 76, 73 and 74% of variability in ETR for Sin WW 2017, 2018 and Sig WW (2018), respectively. The only exception was Sig WW, in which T a was the dominant meteorological variable (74% variability) during 2017. Lower r values ranging from 0.58 to 0.74 were observed with T a and measured ETR. This meteorological variable, which is related to available energy, was the second most relevant factor in determining daily ETR. The sizes of r for T a were followed by RH (scattering throughout −0.50), which was the third meteorological variable that most influenced measured ETR. The only inefficient meteorological variable related to measured ETR was u. However, it is important to point out that the factors controlling ETR interact within the soil-crop-atmosphere continuum, including soybean ETR. Results in our analysis, in accordance with the findings of Liu and Feng (2012) [64], revealed that the most important factors in ETR were VPD (atmospheric water demand) and two energy sources (T a and R n ).
It is worth noting that Sadok and Sinclair (2009) [65] found genetic variability in the ETR to VPD relationship when studying 22 genotypes from a recombined inbred line (RIL) population and their parents in Gainesville, FL, USA (29 • N). They reported that VPD at which plants limit ETR might vary in different genotypes. Depending on the genotype, VPD reached a maximum ETR at a VPD ranging from 1.3 to 2 kPa. The above conclusion was strengthened by Gilbert et al. (2011) [66], who investigated 11 soybean genotypes at Harvard University, Cambridge, MA, USA (42 • N) in 2009.
Every meteorological variable was included in the multiple stepwise regression analysis (Table 6). On the basis of the current analysis, VPD impacted measured ETR the most with only one exception, Sig WW during 2017, where T a became the dominant variable (first equations in Table 6). The second equation in every treatment-except for Sig WW in 2018-were completed with the meteorological variable T a . In the case of Sig WW in 2018, VPD was completed with RH. Table 6. Multiple stepwise regression analysis between meteorological elements (air temperature, T a ; relative humidity, RH; net radiation, R s ; wind speed, u; vapor pressure deficit, VPD) and measured soybean evapotranspiration rate (ETR) in varieties Sinara (Sin) and Sigalia (Sig). WW denotes unlimited watering.

Relationship between Seed Yield and Seasonal Water Use
Due to the higher and more evenly distributed PR in 2018, 14.2% (p < 0.001) and 6.7% (p < 0.001) higher seed yields were produced in the rainfed crop than in the previous (2017) season ( Figure 11). Unlimited watering and partial water withdrawal equalized seed production between the two studied seasons. Irrespective of the season, seed yield was significantly lower with declining water supply. Compared to the RO treatments, the WW crops with unlimited watering produced highest seed yield. The two-season increments in yield for WW ranged from 50.3 (Sig 2017 p < 0.001) to 68.7% (Sin 2018) (p < 0.001) relative to values for RO. The seed yield of rainfed crops were intermediate between the WW and RO production results. Lower, but highly significant decreases in seed yield of P varied between 30 and 40% (p < 0.001), when compared to the results of WW. However, with one exception of Sin (in 2018; p < 0.001), variation in seed yield between rainfed and RO crops was not statistically significant (p = 0.085-0.315).
Regression analysis between seed yield and cumulative ETR of soybean showed a quadratic relationship for all water/variety treatments for the two growing seasons (Figure 12). The second-order polynomial function, which had a determination coefficient of 0.85 (p < 0.001; RMSE: 0.05), implied that irrigation seemed to be an effective tool to increase soybean yield on the current study site. However, it should be noted that continuously high soil water content in an evapotranspirometer may also expose the soybean to harm. Among others, Bajgain et al. (2015) [45] observed in soybean variety Enrei, in Kyoto, Japan (34 • N) between 2010 and 2011, that leaf gas exchange was negatively affected by excessive water, demonstrating a decline in physiological activity, LAI and seed yield by about 25-30%. The temporary increase (seed-filling stage) in water tables negatively affected the growth and seed yield of soybean, variety Fukuyutaka in Fukuoka, Japan (33 • N) between 2011 and 2013 [42]. In contrast, Rhine et al. (2010) [67] identified significant variation in excessive water tolerance of five different soybean cultivars placed into the Missouri Soybean Variety Testing Program at Hayward, Missouri (36 • N), in a three-year cultivar screening trial.
Regression analysis between seed yield and cumulative ETR of soybean showed a quadratic relationship for all water/variety treatments for the two growing seasons (Figure 12). The secondorder polynomial function, which had a determination coefficient of 0.85 (p < 0.001; RMSE: 0.05), implied that irrigation seemed to be an effective tool to increase soybean yield on the current study site. However, it should be noted that continuously high soil water content in an evapotranspirometer may also expose the soybean to harm. Among others, Bajgain et al. (2015) [44] observed in soybean variety Enrei, in Kyoto, Japan (34° N) between 2010 and 2011, that leaf gas exchange was negatively affected by excessive water, demonstrating a decline in physiological activity, LAI and seed yield by about 25-30%. The temporary increase (seed-filling stage) in water tables negatively affected the growth and seed yield of soybean, variety Fukuyutaka in Fukuoka, Japan (33° N) between 2011 and 2013 [41]. In contrast, Rhine et al. (2010) [67] identified significant variation in excessive water tolerance of five different soybean cultivars placed into the Missouri Soybean Variety Testing Program at Hayward, Missouri (36° N), in a three-year cultivar screening trial. In the case of rainfed crops, to estimate the amount of total ETR, a simplified soil water balance method [62] is applied to the top 1 m of soil, which is 0.4 m deeper than the bulk root zone of soybean. AWS (available water storage) data before seeding and just after harvesting were used to assess the total amount of lost water (ETR). The difference between these two AWS data sets supplemented with seasonal PR sums accounted for seasonal ETR in rainfed crops. It is worth noting that the AWS of rainfed crops in the present study (data not shown) were almost greater than the threshold level fixed to 50% depletion of total AWS as reported by Gajic et al. In the case of rainfed crops, to estimate the amount of total ETR, a simplified soil water balance method [62] is applied to the top 1 m of soil, which is 0.4 m deeper than the bulk root zone of soybean. AWS (available water storage) data before seeding and just after harvesting were used to assess the total amount of lost water (ETR). The difference between these two AWS data sets supplemented with seasonal PR sums accounted for seasonal ETR in rainfed crops. It is worth noting that the AWS of rainfed crops in the present study (data not shown) were almost greater than the threshold level fixed to 50% depletion of total AWS as reported by Gajic et al.  [68]. Surprisingly, the WUE values in that study were only moderately lower than those (1.0-1.7 and 0.8-1.1 kg m −3 ) for rainfed (var. TMG1180 RR) and irrigated (var. M8210 Ipro) soybean, respectively, as noted by Lathuilliére et al. (2018) [69] at Capuaba farm, Lucas do Rio Verde, Mato Grosso (13 • S). Probably due to less PR in 2017, the WUE of rainfed Sin and Sig were 34.5% (p < 0.001) and 36.5% (p < 0.001) higher than that in Sin WW and Sig WW, respectively. As a consequence of even and abundant PR events during 2018, the variation in WUE between WW and P ceased completely (Sig), or was reduced by only a few percent in Sin. We assume that the seasonal PR sum during 2018 might be high enough to allow almost maximum seed yield, even in rainfed crops. The lowest WUE values were obtained in WW irrespective of the variety. Rainfed conditions increased soybean WUE by as much as 35-40%, consistent with the results of Montoya et al. (2017) [39] and Candogan et al. (2013) [55]. WUE above 0.75 under semiarid conditions (Nebraska, USA; 41 • N), for three consecutive growing seasons resulted in maximum seed yield in soybean (Payero et al. 2005) [52]. Regarding the WUE results in our study, before drawing conclusions based on earlier investigations, a critical overview of data derived from other locations with different climatic conditions is necessary. WUE of 0.52-0.72 and 0.65-0.66 in well-watered and water-stressed soybean (variety: Talon), respectively, in Bari, Italy (41° N), were slightly lower than the values of 0.77 published by Katerji et al. (2008) [68]. Surprisingly, the WUE values in that study were only moderately lower than those (1.0-1.7 and 0.8-1.1 kg m −3 ) for rainfed (var. TMG1180 RR) and irrigated (var. M8210 Ipro) soybean, respectively, as noted by Lathuilliére et al. (2018) [69] at Capuaba farm, Lucas do Rio Verde, Mato Grosso (13° S). Probably due to less PR in 2017, the WUE of rainfed Sin and Sig were 34.5% (p < 0.001) and 36.5% (p < 0.001) higher than that in Sin WW and Sig WW, respectively. As a consequence of even and abundant PR events during 2018, the variation in WUE between WW and P ceased completely (Sig), or was reduced by only a few percent in Sin. We assume that the seasonal PR sum during 2018 might be high enough to allow almost maximum seed yield, even in rainfed crops. The lowest WUE values were obtained in WW irrespective of the variety. Rainfed conditions increased soybean WUE by as much as 35-40%, consistent with the results of Montoya et al. (2017) [38] and Candogan et al. (2013) [55]. WUE above 0.75 under semiarid conditions (Nebraska, USA; 41° N), for three consecutive growing seasons resulted in maximum seed yield in soybean (Payero et al. 2005) [52]. Regarding the WUE results in our study, before drawing conclusions based on earlier investigations, a critical overview of data derived from other locations with different climatic conditions is necessary.
WUE was influenced by V (p < 0.001), S (p = 0.001) and W (p < 0.001) and some of their interactions (p < 0.001), with two exceptions, non-significant V × S (p = 0.325) and V × S × W (p = 0.713) effects on WUE. The non-significant V × S interaction revealed that the response of seasons on WUE did not differ between both varieties.

Conclusions
The aim of the present study was to investigate ETR and a few related determinants of ETR in two soybean varieties differing in water demands for three different water supply levels over the two growing seasons (2017 and 2018). Among two of the controlled water treatments used, WW crops were grown in an evapotranspirometer that was traditionally operated. The RO crop received 50% of its water requirement through specially converted evapotranspirometers from the R1 developmental stage. The third treatment was rainfed soybean. WUE was influenced by V (p < 0.001), S (p = 0.001) and W (p < 0.001) and some of their interactions (p < 0.001), with two exceptions, non-significant V × S (p = 0.325) and V × S × W (p = 0.713) effects on WUE. The non-significant V × S interaction revealed that the response of seasons on WUE did not differ between both varieties.

Conclusions
The aim of the present study was to investigate ETR and a few related determinants of ETR in two soybean varieties differing in water demands for three different water supply levels over the two growing seasons (2017 and 2018). Among two of the controlled water treatments used, WW crops were grown in an evapotranspirometer that was traditionally operated. The RO crop received 50% of its water requirement through specially converted evapotranspirometers from the R1 developmental stage. The third treatment was rainfed soybean.
There were three differences in soybean's features related to water demands. 1. Irrespective of the variety and season, declines in seasonal mean LAI of RO were two-fold higher than the mean LAI of P. A higher SPAD of RO and P, but not of all treatments, suggested that the seasonal amount/distribution of PR is one of the main weather variables limiting the dry land SPAD value, whose size might enhance photosynthesis and related determinants of crop yield at the study site. Daily mean ETR reported in the current study for WW was, on average, about 65-75% greater than in RO for both varieties. Regarding the two soybean varieties, Sin WW used only 6.1% more water than Sig WW, averaged across the two seasons. Unexpectedly, water stress-tolerant Sin used slightly more water than Sig, the variety that was bred for "normal" weather conditions, and which should be more vulnerable to water-deprived environments. It is important to note that these ETR values were determined under an unlimited water supply in an evapotranspirometer. The cumulative ETR curves were far from the PR totals in the studied area over the two-season period. In both growing seasons, the impact of season on mean K c in both soybean varieties was insignificant.
2. ETR and the impact of meteorological variables on ETR: The most consistent meteorological variable related to soybean ETR was VPD, irrespective of season. The higher ETR totals in 2017 than in 2018 call attention to the importance of being informed about more detailed meteorological variables than monthly (seasonal) means. Despite an approximately 1.5 • C higher seasonal mean T a in 2018, the ranges of daily mean ETR were narrower than during the cooler 2017 growing season. A detailed analysis based on recognized daily variation and extreme T a certified that daily T amax values could also sometimes play a significant role in soybean ETR. The lower VPD in summer also contributed to the lower ETR during 2018. Interestingly, even when T amax and ETR/ evapotranspiration total increased, the WUE in different water supplies hardly varied in 2018.
3. The adaptability of a previously applied technical conversion of evapotranspirometers was confirmed, affirming thus that the reconstructed devices allow water deprivation under field conditions.
In Hungary, seasonal PR sums were insufficient to cover soybean's water demand and irrigation was required in both seasons. From this standpoint, PR sums might greatly constrain water availability for soybean growth at Keszthely. Irrigation is one possible option that soybean growers might have to adapt to account for PR variability which is relevant to areas that are vulnerable to the negative impacts of global climate, change as the Carpathian Basin. Funding: A special thanks to Karintia Ltd. for their kindness by supporting us by donating good-quality, disease-free soybean seed.