Twenty-Years of Hop Irrigation by Flooding the Inter-Row Did Not Cause a Gradient along the Row in Soil Properties, Plant Elemental Composition and Dry Matter Yield

: In hops ( Humulus lupulus L.), irrigation by ﬂooding the inter-row can carry away suspended particles and minerals, causing gradients in soil fertility. The effect of more than 20 years of ﬂooding irrigation on soil and plants was evaluated in two hop ﬁelds by measuring soil and plant variables in multiple points along the rows. In a second experiment 1000 kg ha − 1 of lime was applied and incorporated into the soil to assess whether liming could moderate any gradient created by the irrigation. At different sampling points along the rows, signiﬁcant differences were recorded in soil properties, plant elemental composition and dry matter yield, but this was not found to exist over a continuous gradient. The variations in cone yield were over 50% when different sampling points were compared. However, this difference cannot be attributed to the effect of irrigation, but rather to an erratic spatial variation in some of the soil constituents, such as sand, silt and clay. Flooding irrigation and frequent soil tillage resulted in lower porosity and higher soil bulk density in the 0.0–0.10 m soil layer in comparison to the 0.10–0.20 m layer. In turn, porosity and bulk density were respectively positively and negatively associated with crop productivity. Thus, irrigation and soil tillage may have damaged the soil condition but did not create any gradient along the row. The ridge appeared to provide an important pool of nutrients, probably caused by mass ﬂow due to the evaporation from it and a regular supply of irrigation water to the inter-row. Liming raised the soil pH slightly, but had a relevant effect on neither soil nor plants, perhaps because of the small amounts of lime applied.


Introduction
Hop plants (Humulus lupulus L.) require an adequate supply of water during the growing season to sustain their huge canopy [1]. In most of the hop producing regions of the world, the crop needs to be irrigated, particularly in lower latitudes of reduced precipitation in summer. Although hop fields have started to be drip irrigated all over the world, there is a long tradition of surface watering of this crop, by flooding the space between rows [1,2]. In this kind of surface or furrow irrigation system, water is applied at the top end of each furrow (in hops to the inter-row space) and flows down the field under the influence of gravity [3]. This is still the most commonly used irrigation method for hops in northern Portugal [4]. The water use efficiency with this irrigation technique is highly dependent on the field gradient and water infiltration rate, which can vary considerably, inducing spatial and temporal variability in the main soil properties [5]. In addition, flood irrigation can affect the spatial distribution of soil physicochemical properties which may exacerbate the spatial variability in crop growth and yield [6].
Flood irrigation can have a major impact on soil properties by varying salinity, redox potential, compaction and/or porosity [7][8][9][10]. Furthermore, hop fields which are Flood irrigation can have a major impact on soil properties by varying salinity, redox potential, compaction and/or porosity [7][8][9][10]. Furthermore, hop fields which are floodirrigated need to be frequently tilled to control summer weeds and to reduce soil compaction and superficial crusts in the short term. This allows a better infiltration of water, but that can also have a negative impact on the soil in the long term [11,12]. Soil compaction, increased by furrow irrigation, may also reduce soil drainage and aeration, contributing to the reduction of soil redox potential which influences soil chemistry and plant nutrient availability [10,13]. The degree of compaction of a soil can be assessed by measuring some physical properties, such as bulk density and porosity [12,13]. As the soil becomes more compact, bulk density increases and soil porosity decreases, which reduces water and air diffusion into the soil [11,14]. In some hop fields in northern Portugal it was found that the decrease in soil redox potential, associated with an excess of water and/or poor drainage, was the main cause of the spatial variability found in crop growth and yield [4].
Soil pH is another relevant issue in hop production. The range of pH most suited for growing hops is considered to be between 5.7 and 7.5 [15,16]. The application of lime is recommended for acidic soils, and a positive relationship has been found between the increase in soil pH and hop yield [15,17]. However, the effect of liming on crops can also vary with the irrigation system. Some researchers have studied the influence of liming in rice under flooding conditions, since great interactions between flooding, soil acidity and nutritional disorders are usually found [18][19][20]. In hops, these interactions are less well known, or the response to liming, but it is believed that it may be relevant enough to be studied, since the crop continues to be irrigated by flooding in several parts of the world.
This study evaluated the variation in soil properties and nutritional status and the productivity of hop plants created along the rows by flooding irrigation. As a second line of study, the effect of the application of lime on soil properties and on hop nutritional status, growth and yield was evaluated, to ascertain if the application of lime could compensate for the variability created by the irrigation system. Both lines of study were carried out in commercial hop fields which had been flood-irrigated for over 20 years.

General Experimental Conditions
The field experiments were carried out during two growing seasons (2017 and 2018) on a commercial farm located in Pinela (41°40′33.6″ N; 6°44′32.7″ W), Bragança, northeastern Portugal. A detailed location of field experiments is shown in Figure 1. and the two hop fields identified in this study as field 1 and field 2 (right). Images from https://www.google.com/maps/place/Pinela (accessed on 7 July 2021). Figure 1. Map of Portugal indicating Pinela (left) and the two hop fields identified in this study as field 1 and field 2 (right). Images from https://www.google.com/maps/place/Pinela (accessed on 7 July 2021).
The region benefits from a Mediterranean-type climate, with an annual average temperature and accumulated precipitation of 12.7 • C and 772.8 mm, respectively. The average monthly temperatures and precipitation recorded during the experimental period are shown in Figure 2.
Horticulturae 2021, 7, x FOR PEER REVIEW 3 of 16 The region benefits from a Mediterranean-type climate, with an annual average temperature and accumulated precipitation of 12.7 °C and 772.8 mm, respectively. The average monthly temperatures and precipitation recorded during the experimental period are shown in Figure 2. The hop plots where the study was undertaken are ⁓2 ha in size each, with the rows having a length ranging from 150 to 180 m, and established with the cultivar Nugget. The fields are arranged in a 7 m conventional high trellis system, with concrete poles connected with steel cables, in a "V" design system. The farmer has managed the fields by flooding irrigation since the hop crop was installed more than 20 years ago. Several tillage passes (3 to 4) are performed every year to remove the crusts and facilitate water infiltration. The fertilization programme includes the application of a compound nitrogen (N): phosphorus (P): potassium (K) fertilizer (7:14:14, 7% N, 14% P2O5, 14% K2O) early in the spring, followed by two applications of N fertilizer (ammonium nitrate, 27% N) as a side-dressing, totalling ⁓150, 44 and 83 kg ha −1 of N, P and K, respectively. The farmer also follows a phytosanitary programme for crop protection against pests and diseases.

Field Experiments and Soil and Plant Sampling
The first experiment (Experiment 1) was carried out during the growing season of 2017 in two hop fields. It consisted of the evaluation of soil properties, plant nutritional status and crop yield, searching for any gradient along the rows created by the irrigation system. The rows used in this experiment were divided into nine segments of equivalent length, creating nine positions (P1, P2, …, P9) for soil and plant sampling. The soil was sampled between rows and on the ridges to a depth of 0.0 to 0.2 m. Three rows and interrows of hops were used to create three replicates for each position. Each soil sample for analysis resulted from six sampling points (composite samples). The soil was sampled by using an open-face auger.
For the determination of soil bulk density and porosity, a different approach to soil sampling was followed. It was found unnecessary to sample in the ridges since no compaction was expected in this part. Instead, the soil was sampled at two different depths, 0.0 to 0.10 m and 0.10 to 0.20 m. Due to the increased difficulty of sampling, particularly in the 0.10 to 0.20 m layer, only five positions were considered (P1, P3, P5, P7 and P9) and sampled in three replicates. For these analyses, undisturbed soil cores were taken by using appropriate cylinders of 100 cm 3 . Soil samplings were carried out on 10 March 2017.
The plants used in this experiment for the evaluation of their nutritional status and crop productivity were randomly selected and marked when plant height was close to 3 The hop plots where the study was undertaken are~2 ha in size each, with the rows having a length ranging from 150 to 180 m, and established with the cultivar Nugget. The fields are arranged in a 7 m conventional high trellis system, with concrete poles connected with steel cables, in a "V" design system. The farmer has managed the fields by flooding irrigation since the hop crop was installed more than 20 years ago. Several tillage passes (3 to 4) are performed every year to remove the crusts and facilitate water infiltration. The fertilization programme includes the application of a compound nitrogen (N): phosphorus (P): potassium (K) fertilizer (7:14:14, 7% N, 14% P 2 O 5 , 14% K 2 O) early in the spring, followed by two applications of N fertilizer (ammonium nitrate, 27% N) as a side-dressing, totalling~150, 44 and 83 kg ha −1 of N, P and K, respectively. The farmer also follows a phytosanitary programme for crop protection against pests and diseases.

Field Experiments and Soil and Plant Sampling
The first experiment (Experiment 1) was carried out during the growing season of 2017 in two hop fields. It consisted of the evaluation of soil properties, plant nutritional status and crop yield, searching for any gradient along the rows created by the irrigation system. The rows used in this experiment were divided into nine segments of equivalent length, creating nine positions (P1, P2, . . . , P9) for soil and plant sampling. The soil was sampled between rows and on the ridges to a depth of 0.0 to 0.2 m. Three rows and inter-rows of hops were used to create three replicates for each position. Each soil sample for analysis resulted from six sampling points (composite samples). The soil was sampled by using an open-face auger.
For the determination of soil bulk density and porosity, a different approach to soil sampling was followed. It was found unnecessary to sample in the ridges since no compaction was expected in this part. Instead, the soil was sampled at two different depths, 0.0 to 0.10 m and 0.10 to 0.20 m. Due to the increased difficulty of sampling, particularly in the 0.10 to 0.20 m layer, only five positions were considered (P1, P3, P5, P7 and P9) and sampled in three replicates. For these analyses, undisturbed soil cores were taken by using appropriate cylinders of 100 cm 3 . Soil samplings were carried out on 10 March 2017.
The plants used in this experiment for the evaluation of their nutritional status and crop productivity were randomly selected and marked when plant height was close to 3 m (to avoid using very atypical plants) and close to each of the positions used for soil sampling. Leaf sampling for crop nutritional status assessment was done at~2 m in height, on 17 July 2017. At harvest (1 September 2017), plant biomass was cut at ground level. Subsequently, the aboveground biomass was separated into leaves, stems and cones and weighed fresh. Subsamples of each plant part were weighed fresh again and then oven-dried at 70 • C and weighed dry for determination of dry matter yield.
The second experiment (Experiment 2) consisted of the application of 1000 kg ha −1 of lime (55% CaCO 3 , 28% CaO and 20% MgO) in February 2017, to assess the liming effects on soil properties and plants in comparison to the untreated control. This experiment was also carried out in two hop plots. The general methodology for soil and plant sampling was similar to that reported for Experiment 1, consisting of marking nine positions along the rows. The soil was sampled on 4 January 2019, only between the rows, at 0.0-0.20 m soil depth, using an open-face auger. Leaf samples were taken at~2 m in height, on 17 July 2017 and 18 July 2018. At harvest (1 September 2017 and 31 August 2018), plant biomass was cut at ground level and treated as reported for Experiment 1.

Laboratory Analyses
The undisturbed soil samples from Experiment 1 were oven-dried at 105 • C and weighed. Soil bulk density was estimated from the weight of dry soil divided by the volume of the cylinder. Soil porosity was determined as the ratio of nonsolid volume (soil particle density-bulk density) to the total volume of soil (soil particle density) [21]. The other soil samples from Experiments 1 and 2 were oven-dried at 40 • C and sieved in a mesh of 2 mm. The samples were analysed for pH (H 2 O and KCl), electrical conductivity (soil:solution, 1:2.5), exchangeable complex (ammonium acetate, pH 7.0) and organic carbon (C) (Walkley−Black method). Extractable P and K were determined by a combination of ammonium lactate and acetic acid buffered at pH 3.7. Soil boron (B) was extracted by hot water and the extracts analysed by the azomethine-H method. More details of these analytical procedures are given in Van Reeuwijk [22]. Other micronutrients [copper (Cu), iron (Fe), zinc (Zn), and manganese (Mn)] were determined by atomic absorption spectrometry after extraction with ammonium acetate and EDTA, following the methodology reported by Lakanen and Erviö [23].
Tissue samples (leaves, stems and cones) from both experiments were oven-dried at 70 • C and ground. Elemental tissue analyses were performed by Kjeldahl (N), colorimetry (B and P), flame emission spectrometry (K) and atomic absorption spectrophotometry (calcium (Ca), magnesium (Mg), Cu, Fe, Zn and Mn) methods after nitric digestion of the samples [24].

Data Analysis
Data was subjected to analysis of variance, according to the experimental designs, using SPSS program version 25. When significant differences were found between the experimental treatments, the means were separated by the Tukey HSD (sampling position) and Student's-t (field, sampling site, lime treatment) tests (α = 0.05). Linear regression analysis was performed to understand the effects of gradient on soil properties and plant nutritional status and productivity in Experiment 1 and the relationship between soil pH and plant variables in Experiment 2. The relation between the variables was obtained through correlation analysis with the Pearson coefficient, when the assumption of normality and linearity was accomplished; when this was not the case, the Spearman coefficient was used.

Soil Properties
The silt and sand contents varied significantly between the sampling positions ( Table 1). The two fields also differed significantly in clay and sand content. The soil bulk density and soil porosity varied significantly between the sampling positions and fields but in the opposite way. The interaction between sampling position and field was significant for soil porosity, which means that the effect of the irrigation on this variable depended on the field. Table 1. Soil separates and soil bulk density and porosity from samples collected at 0.0-0.20 m depth, in March 2017, as a function of sampling position (1, . . . , 9), and field. Means followed by the same letter are not significantly different by Tukey HSD (sampling position) or Student's t (field) tests (α = 0.05). The soil bulk density was higher in the soil surface (0.0-0.1 m) when compared to the deeper (0.1-0.2 m) layer ( Figure 3). The soil bulk density did not vary significantly along the rows for both soil depths. The soil porosity, in turn, was lower in the surface layer, and the gradient found along the rows was not significant for any of the soil layers. The soil bulk density and porosity varied significantly between the two fields, but the gradients found along the rows were not statistically significant. 1). The two fields also differed significantly in clay and sand content. The soil bulk density and soil porosity varied significantly between the sampling positions and fields but in the opposite way. The interaction between sampling position and field was significant for soil porosity, which means that the effect of the irrigation on this variable depended on the field. Table 1. Soil separates and soil bulk density and porosity from samples collected at 0.0-0.20 m depth, in March 2017, as a function of sampling position (1, …, 9), and field. Means followed by the same letter are not significantly different by Tukey HSD (sampling position) or Student's t (field) tests (α = 0.05). The soil bulk density was higher in the soil surface (0.0-0.1 m) when compared to the deeper (0.1-0.2 m) layer ( Figure 3). The soil bulk density did not vary significantly along the rows for both soil depths. The soil porosity, in turn, was lower in the surface layer, and the gradient found along the rows was not significant for any of the soil layers. The soil bulk density and porosity varied significantly between the two fields, but the gradients found along the rows were not statistically significant. Some other soil properties determined from the samples collected at 0.0-0.20 m depth varied significantly between sampling sites, sampling positions and fields ( Table 2). Extractable P and K, conductivity, organic C, CEC and extractable Zn and B showed significantly higher values in the samples collected in the ridges. However, soil pH (H 2 O and KCl), base saturation and extractable Mn were significantly higher in the samples collected in the inter-rows. Most of the soil properties varied significantly between the sampling positions, the exceptions being soil pH, conductivity and extractable K. Soil properties also differed significantly between fields, except for soil conductivity. Significant interaction between the sampling site and the field was found for extractable P, conductivity, exchangeable Ca and extractable Fe. Significant interaction between the sampling position and the field was found for organic C and extractable Fe, Mn, Zn, Cu and B. No significant interaction was found between the three factors of this experiment.

Hop Dry Mater Yield and Leaf Nutrient Concentration
Aboveground dry biomass (stems, leaves, cones and total) in Field 1 showed a clear tendency for a decrease along the rows (Figure 4). However, the decrease was only statistically significant for stem dry matter yield (DMY). For all plant parts, the coefficients of determination (R 2 ) were not particularly high, which helps to explain the lack of significant correlation between the two variables. In Field 2, no clear tendency was found in aboveground DMY. Aboveground dry biomass (stems, leaves, cones and total) in Field 1 showed a clear tendency for a decrease along the rows (Figure 4). However, the decrease was only statistically significant for stem dry matter yield (DMY). For all plant parts, the coefficients of determination (R 2 ) were not particularly high, which helps to explain the lack of significant correlation between the two variables. In Field 2, no clear tendency was found in aboveground DMY. N concentration in the leaves taken at 2 m height did not vary significantly along the rows in any of the fields ( Figure 5). Leaf P also did not vary significantly along the rows but the values in Field 1 were lower than in Field 2. Leaf K levels did not vary significantly along the rows in Field 1 but increased significantly in Field 2. Leaf Ca and Mg levels showed a slight tendency to increase in both fields but without statistical significance. In general, the micronutrients showed even more erratic tendencies when the values of the two fields were compared and only the values of leaf Cu showed a significant decrease along the rows in Field 1.  N concentration in the leaves taken at 2 m height did not vary significantly along the rows in any of the fields ( Figure 5). Leaf P also did not vary significantly along the rows but the values in Field 1 were lower than in Field 2. Leaf K levels did not vary significantly along the rows in Field 1 but increased significantly in Field 2. Leaf Ca and Mg levels showed a slight tendency to increase in both fields but without statistical significance. In general, the micronutrients showed even more erratic tendencies when the values of the two fields were compared and only the values of leaf Cu showed a significant decrease along the rows in Field 1.

Correlation Analysis between Soil Properties and Plant Dry Matter Yield
Soil bulk density and porosity correlated in a different way with soil pH (H2O and KCl), leaf P and total DMY (Table 3). That is, the correlations of soil pH were positive for soil bulk density and negative for soil porosity at 0.0-0.10 m depth. Leaf P concentration was significantly and negatively correlated with soil bulk density at 0.10-0.20 m depth, in contrast to the positive correlation found with soil porosity. Leaf Fe concentration was found significant and negatively correlated only with soil porosity at 0.10-0.20 m depth. The strongest correlations were found for total DMY with soil bulk density (r = −0.706) and soil porosity (r = 0.714), both at 0.10-0.20 m depth.

Correlation Analysis between Soil Properties and Plant Dry Matter Yield
Soil bulk density and porosity correlated in a different way with soil pH (H 2 O and KCl), leaf P and total DMY (Table 3). That is, the correlations of soil pH were positive for soil bulk density and negative for soil porosity at 0.0-0.10 m depth. Leaf P concentration was significantly and negatively correlated with soil bulk density at 0.10-0.20 m depth, in contrast to the positive correlation found with soil porosity. Leaf Fe concentration was found significant and negatively correlated only with soil porosity at 0.10-0.20 m depth. The strongest correlations were found for total DMY with soil bulk density (r = −0.706) and soil porosity (r = 0.714), both at 0.10-0.20 m depth. Significant correlations were found for soil clay content, positive for soil pH and leaf K, and negative for leaf P, leaf Cu, total DMY and cone DMY. In contrast, soil sand content correlated significantly and negatively with soil pH (H 2 O) and leaf K, and positively with leaf Cu and B. Soil silt content correlated significantly and negatively with leaf B.

Soil Properties
Most soil properties, such as extractable K, P, Mn, Zn, Cu, B, conductivity and pH presented significantly higher values in the limed plot in comparison to the untreated control (Table 4). Exchangeable Ca and CEC showed higher values in the limed plot but not significantly different to those observed in the control. Significant differences between the two fields used in this experiment were also found for most of the soil properties, the values of extractable K, P, Zn, Cu and B, conductivity, pH, exchangeable Ca, CEC and base saturation being significantly higher in Field 1. Only extractable Fe was significantly higher in Field 2. The interaction between liming and field was significant for extractable K, conductivity, and pH.

Plant Response to Liming
The concentration of nutrients in the leaves taken at 2 m height showed significant differences between treatments for leaf P in 2017 and for leaf Fe and B in 2018 ( Table 5). The values reported for P and Fe were significantly higher in the limed plots, and those reported for B were significantly higher in the control. Total and cone DMY were significantly lower in the limed plots with the exception of total DMY in 2017, whose differences between treatments were not statistically significant. When comparing fields, significant differences were found for some nutrients and total and cone DMY. However, only leaf concentrations of K, Cu and B, and total and cone DMY, maintained the same trend in both years and fields. In 2017, significant interaction between the liming treatment and the field was only found for leaf N and Mn and in 2018 for leaf P and total DMY.

Correlation Analysis between Soil pH and Plant Variables
Significant correlations between the soil pH (H 2 O and KCl) and leaf nutrient concentration were found for several nutrients, but a similar trend over the two years was found only for leaf Cu and B, both presenting negative correlations with soil pH (Table 6). Soil pH and leaf P, for instance, showed a negative correlation in 2017 and a positive correlation in 2018. Significant and negative relations between soil pH and total and cone DMY were also found for the first year of plant sampling. Table 4. Soil properties from samples collected at 0.0-0.2 m depth, in January 2019, in the inter-rows, as a function of liming treatment and field. Means followed by the same letter are not statistically different by Student's t test (α = 0.05).

Discussion
The results of Experiment 1 showed significant differences in some soil properties at different positions along the rows, but not over a continuous gradient. Thus, the results cannot be attributed to the flooding irrigation, but they were probably caused by heterogeneity in spatial variability of important soil constituents such as clay, sand and silt, since it is well-known that soil texture determines many other soil physical and chemical properties [25]. Variations in soil properties were also found when comparing different soil layers. The soil bulk density was higher in the soil surface layer (0.0-0.1 m), and porosity was found to be higher in the deeper (0.1-0.2 m), layer. The soil bulk density and porosity in agricultural fields are influenced not only by soil texture but also by external loads which cause soil compaction [13,26,27]. In this particular case, it seems that the effects of frequent irrigation and soil tillage prevailed, which may have prevented a proper soil aggregation, leading to an increase in soil bulk density and a reduction in soil porosity on the surface layer which was directly impacted by the cultivator. The variation in soil properties was also significant when comparing fields. The field higher in clay and lower in sand presented significantly higher soil bulk density. Usually, clayey soils tend to have a lower bulk density and higher porosity than sandy soils [28]. However, these results indicate an opposite trend, probably because of the negative effect on soil aggregation and compaction caused by frequent soil tillage. Other studies have also found spatial variability in bulk density and water infiltration on flooded fields caused mainly by tillage practices, particularly when heavy machinery is used [8,29].
The soil samples collected from the ridges showed significantly higher values of extractable P and K. In the ridges, the conditions for nutrient uptake were poor since they are created every year by soil pushed from the inter-rows, which means that they contain nutrients barely taken up by the plant due to the limited expansion of roots in this position. In addition, in this irrigation system, the water flows from the inter-row to the ridge due to the gradient of water potential caused by the evapotranspiration from the latter and the continuous water supply to the inter-row. This means that nutrients tend to accumulate in the ridge, carried by mass flow, in contrast to what happens in the inter-row, from which nutrients tend to be leached out. Mass flow is the main driving force causing the movement of most nutrients in the soil [30][31][32]. Thus, soil conductivity was higher in the ridge, due to the increased presence of salts as demonstrated by the increase in CEC. Organic C also appeared higher in the ridge, probably because this zone is not tilled so frequently, which reduces the exposure of organic matter to the heterotrophic microorganisms that cause its oxidation [33]. This zone also contains the remaining bines (those that do not climb) and weeds, which are incorporated into the soil when the ridge is created, which usually represent more debris than that incorporated in the inter-row. B also increased in the ridge, perhaps due to higher levels of organic C, which have the ability of retaining B in the soil [34,35].
Soil pH (H 2 O and KCl), base saturation and extractable Mn were significantly higher in the samples collected in the inter-rows. These results are probably related to the decrease in the potential redox, which may have increased the pH of the soil [36]. The increase in soil pH in the inter-rows was probably also related to the increase in the concentration of cation ions, such as Ca and Fe [37]. Base saturation increased in the inter-rows probably due to the presence of the divalent cations, less available to move into the ridge by mass flow. The higher concentration of Mn in the inter-rows might have also been due to the reduction of Mn that occurred at the beginning of the reduction process. This can occur when the redox potential is still positive [38].
A clear gradient along the rows was not observed for total and cone DMY. These results did not corroborate the hypothesis that flood irrigation is creating a spatial variation in plant performance along the rows. The differences detected in the plants seem to be due to spatial variability in the soil constituents, namely the soil separates which, in turn, influence soil bulk density and porosity. The results from the correlation analysis showed significant and negative relations between total DMY and soil bulk density (r = −0.706) and between total DMY and clay content (r = −0.676). In contrast, total DMY and soil porosity at 0.10-0.20 m correlated significantly and positively (r = 0.714). The soil surface layer presented a higher bulk density, which has already been explained by the effect of irrigation and frequent soil tillage, which reduces the stability of soil aggregates, increasing bulk density and decreasing porosity [10,26]. On the other hand, it seems that the higher porosity in the 0.10-0.20 m layer was an important factor affecting DMY, likely because in the surface layer the diffusion of oxygen to ensure the biological processes of the soil is always easier. Soils with a higher clay content tend to retain more water, decreasing soil aeration which negatively affects the function of root and plant metabolism [13,39]. Under the conditions of this experiment, the clay content in the soil seemed to be negatively associated with hop DMY, mainly because clay is a determinant factor of soil bulk density and porosity, which were identified in this study as determinant factors in crop productivity.
Irrigation also did not cause any relevant gradient in tissue nutrient concentration as detected by the analysis of variance. However, correlation analysis provided some data that deserves to be commented on. Leaf P was significantly and positively correlated with soil porosity at 0.10-0.20 m, but was negatively correlated with soil bulk density at 0.10-0.20 m and clay content. Leaf P did not show any consistent gradient along the rows, but was lower in the field presenting a higher soil bulk density and clay content. This reveals that P uptake was enhanced by the increased porosity of the soil at the deeper layer and by the lower clay content. Similarly, on barley (Hordeum vulgare L.) there was reported a reduction in P uptake and yield associated with heavy soil compaction [13]. The higher porosity of soil may have facilitated P root uptake from the deeper layer, which is richer in P, probably due to the increase in the vertical movement of P as the result of fertilization and flooding as reported by [40]. In turn, the higher clay content may have resulted in higher P adsorption and lower P availability. In contrast, leaf Fe was significantly and negatively correlated with soil porosity at 0.10-0.20 m depth. Leaf Fe also presented an opposite tendency between fields, decreasing along the rows in the field with a lower clay content and higher soil porosity. This result is probably related to soil reduction conditions, as the availability of Fe decreases when soil oxygen and redox potential increases [37]. Leaf K showed a significant and positive correlation with soil clay content and a negative one with sand content. The availability of K in the soil is not directly affected by redox potential, but its fixation in 2:1 clay minerals is facilitated by the increase in soil pH [36]. There has also been reported an antagonistic effect between Fe and K in paddy fields [41,42], an aspect that may also have influenced these results.
In Experiment 2, the application of lime increased several variables of soil fertility, including pH, but did not significantly increase exchangeable Ca and CEC. In fact, the rate of lime applied in this experiment was too low to cause important changes to soil properties, as is usually achieved when using high rates of lime [32]. In a previous study, Ceh andČremožnik [17] applied 2.3 t lime ha −1 and reported similar results, that is, a reduced effect on soil properties due to the application of lime.
The main effect on the elemental composition of the leaves resulting from the application of lime would have been the significant increase of leaf P in the first growing season after the lime application. This raised the soil pH contributing to a reduction in P fixation, which in acidic soils is due to reactions with Al and Fe oxides, which precipitate P as AlPO 4 and FePO 4 [43].
Total and cone DMY did not increase with the application of lime, but rather showed a decreasing trend. It is generally considered that the optimal pH for hop growth is between 5.7 and 7.5 [15,44]. In this study, soil pH was below the lowest value of the reported range, which would have favoured a positive effect on the vegetation. However, the lime application influenced some soil properties, but not enough to have a high impact on the elemental composition of the leaves. In general, the nutrient content of the leaves was found to be within the sufficiency ranges established for hops [45], both in the limed and in the control treatments. Regarding total DMY, a significant interaction between lime treatment and field was recorded, which may also have contributed to difficulties in the interpretation of these results.
Correlation analysis, in turn, also did not show coherent trends over the two years of the study. Perhaps the most relevant result was the negative correlation between soil pH and biomass production in the first year, which again refers to diverse interactions which may have occurred between environmental variables (year) and factors under study (field and liming). The effect of environmental variables on the performance of the hop plant is well known [46][47][48], although in this study it was not possible to clarify the isolated effects of any of them.

Conclusions
Irrigation by flooding the space between rows over more than 20 years was not responsible for any gradient in soil properties, plant elemental composition and plant performance, although variations in those variables were found at different positions in the row caused by erratic spatial variability of some constituents of the soil, such as sand, silt and clay. However, irrigation followed by soil tillage on repeated occasions during the growing season seems to have reduced soil porosity and increased soil bulk density in the surface 0.0-0.1 m soil layer. These variables were found to be related to crop productivity in positive and negative ways, respectively.
This study also showed that the ridge is a point of nutrient accumulation, particularly for those that move more easily in the soil by mass flow, thereby showing also higher conductivity and CEC. The reduced water potential in the ridge created by evapotranspiration is the driving force causing the water flow from the inter-row. Organic C was also higher in the ridge in comparison with the inter-row, probably due to the annual incorporation of weeds and weaker hop bines (those that did not climb) when the ridge is created in early spring.
Although the original soil was acidic, and the application of 1000 kg ha −1 of lime caused a small increase in pH, this did not lead to other relevant changes in soil properties, nor in plant nutrition status or total and cone DMY. The liming effect might not have been enough to nullify the effects of the interaction between factors that always occur in field experiments.