Mechanical Pruning of ‘Clemenules’ Mandarins in Spain: Yield Effects and Economic Analysis

Pruning is one of the most expensive tasks in citrus production, and its mechanization could increase the productivity and competitiveness of citrus farms. The effect of mechanical pruning on yield depends on the variety, crop condition, and location; among other factors. The ‘Clemenules’ mandarin variety is one of the most important ones; therefore, the aim of this work was to study the effects of twelve pruning strategies on ‘Clemenules’ yield over four years, and to conduct an economic analysis. These strategies included fully manual pruning, annual alternation of mechanical/manual pruning, mechanical pruning with manual follow-up, and fully mechanical pruning with different types of cuttings. The results showed that pruning strategies affected the yield and fruit size in some years of the study but not when the four years overall are considered. Strategies that annually alternated mechanical pruning with manual pruning had higher yields than the fully mechanical approaches, and similar yields to the manual and control (no pruning) strategies. Furthermore, although the mechanical pruning costs per hectare were as minimum two times less than those with manual pruning, higher net value was achieved with the manual pruning approach followed by alternate annual mechanical/manual pruning. In conclusion, after four years of study, the pruning strategy with the best global results for ‘Clemenules’ mandarin was that which involved annually alternating mechanical (topping and one-sided hedging) and manual pruning. This strategy allows maintaining regular tree dimensions without dry branches inside.


Introduction
Spain is the sixth largest world producer of citrus and the biggest exporter of citrus for fresh consumption. The Spanish citrus area takes up around 300,000 ha, with a production of 6.2 × 10 6 t in 2019, and with 54% of this production being of oranges, 31% of mandarins, and 18% of lemons. With a production of 1.9 × 10 6 t of mandarins (tangerines, clementines, and satsumas) in 2019, Spain is the world's second highest producer, after China, and is followed by Turkey and Morocco [1].
In the last decades, the crop experienced important economic problems because fruit prices stalled or even decreased, while costs rose continuously [2,3]. A reduction in production costs is necessary to guarantee the profitability of farmers as well as to compete with other countries, where the production costs, mainly manpower costs, are lower because of the lower salaries [4,5].
A reduction in costs must be achieved to increase the competitiveness of citrus producers, and the mechanization of farming practices is a key factor [6]. The current level of mechanization on Spanish citrus orchards is lower than for other 3D crops such as olives, apples, grapes, etc. [7], mainly due to the focus being on marketing high-quality fresh In California, Kallsen [33] observed that severe manual pruning in 'Frost Nucellar' navel orange significantly reduced total yield and fancy fruit with respect to moderate pruning treatments or no pruning. Severe topping reduced total yield, although no significant differences were found between the intensity levels of topping treatments with respect to yield, number of fruits, and fruit size. Over the four years, the mechanically pruned trees produced as much or more fruits of the most valuable sizes than the severe or moderate manual treatments, with lower pruning costs.
Mendonça et al. [34] observed that after an initial year of severe topping pruning in 'Ponkan' tangerines adult trees in which topping reduced yields, in the following years, the plants recovered. In addition, the plants that were skirted presented higher yields than those that were not trimmed.
Yildirim et al. [35] obtained the highest production levels of grapefruits when using mechanical pruning treatments that combined topping and hedging, and the worst production was observed in those approaches that consisted only of topping.
Velázquez and Fernández [36] found that mechanical pruning with manual follow-up increased yields by 30 kg·tree −1 as compared to manual pruning, but with similar costs, which meant an increase in profits. The productivity of the mechanical pruning systems ranged between 1.3 and 5.4 h·ha −1 , manual follow-up pruning required 66.6-80.5 h·ha −1 , and manual pruning needed 107.2 h·ha −1 ; consequently, all systems involving mechanical pruning showed the ability to reduce the time needed to complete the task.
Martin-Gorriz et al. [22] found that the only treatment that significantly reduced yield compared to manual pruning was the one that applied continuous mechanical pruning during the three years and justified this drop with the fact that the pruning was carried out with the fruits already set, because it is a late harvesting variety. However, when mechanical and manual pruning were alternated in consecutive years, in the manual pruning year, the trees were able to recover from the loss of yield that corresponded to the year of mechanical pruning. Similar trials of mechanical pruning combined with manual pruning in 'Fino' lemons showed that the manual pruning was surpassed in yield by the treatment consisting in mechanical pruning all years (topping, skirting, and one-sided hedging) and by the treatment that annually alternated mechanical pruning (topping, skirting, and bilateral hedging) with manual pruning [10]. Both treatments, in addition to improving yield, substantially reduced costs and pruning time compared to the manual system. The treatments with mechanical pruning combined with manual follow-up pruning did not improve yields, nor did they substantially reduce execution costs.
Chueca et al. [11] carried out a similar trial with orange trees var. 'Navel Foyos', in which no significant differences in yield between the manual pruning and the four mechanical pruning treatments, alone and combined with manual follow-up, were found. However, all the treatments that included mechanical pruning were associated with substantially reduced costs, which decreased from 502.8 €·ha −1 for the manual treatment to 70.4 €·ha −1 in the mechanical pruning treatment that included bilateral skirting and topping and one-sided hedging.
In general, mechanical pruning reduces pruning costs and execution times, but its effects on production can be highly variable, depending on the species and varieties, as well as the type (hedging, topping, or skirting) and intensity of the pruning applied; for this reason, it is not feasible to obtain general formulas and it is instead necessary to experiment for each variety and cultivation system, since factors such as the vigour of the plant and the pruning calendar are also very important.
For these reasons, the objective of this work was to study the effects of different strategies of pruning, involving different combinations of mechanical and/or manual pruning, as well as annual alternations between manual and mechanical pruning, on the yield and the fruit size of 'Clemenules' mandarins (Citrus clementina Hort. ex Tan.) in the Mediterranean area, because there is a lack of experience in mechanical pruning for this variety, which is one of the most widespread in this area. Additionally, the working capacities and the economic costs of the different pruning strategies were studied.

Experimental Site and Vegetal Material
Trials were conducted in a commercial orchard of 'Clemenules' clementine mandarins (on Citrange Carrizo rootstock) (Revacitrus S.L.). It was located in Chiva (Valencia, eastern Spain) (geographic coordinates: 39 • 26 30" N, 0 • 32 60" W; altitude: 93 m above sea level). The trees were planted in 1998, in rows directed North-South with trapezoidal-shaped ridges of 0.57 m height. The orchard had a planting frame with 6 m row spacing and 3 m tree spacing (with a tree density of 556 trees·ha −1 ).

Experimental Design
A randomized block design with five repetitions was followed for four consecutive seasons, from 2017 to 2020. The study factor was the pruning strategy, with twelve levels that are defined in Table 1 and shown in Figure 1. The response variables were the citrus production (yield and fruit size), the pruning working capacity, the economic costs of the pruning operations, and the economic profit (net value).
Manual pruning was carried out by specialized workers who decided the branches to cut and prune from the internal to the external part of the tree, removing dried branches and some main branches to open the trees. Mechanical pruning consisted of cutting the external branches and, depending on the strategy, included hedging and topping. Hedging consisted of a vertical cut of the corresponding side of the tree, depending on the strategy and year, with a vertical angle of 5 • towards the upper part of the canopy. Topping was conducted with two oblique cuttings, with a horizontal angle between 5 and 15 • for each side of the canopy (East and West).
The same strategy was applied to the same trees during the four years ( Figure 1). All the trees in the experiment were subjected to the same fertilization, irrigation, and phytosanitary treatments.
In the second year of the experiment, strategies 4, 7, 9, and 11 were divided into strategies 5, 8, 10, and 12, respectively, of annual alternation between mechanical pruning and manual pruning (Figure 1). A minimum of 10 consecutive trees were used for each strategy and repetition, except for the control, in which each repetition consisted of five consecutive trees. All the evaluations were carried out in a representative tree that was randomly selected from the experimental unit (except the extreme trees) in each strategy and repetition, except in the control, in which the central tree was used for this purpose.  Legend: 1 C = control (no pruning); 2 M = manual pruning; 3 TFF = topping and manual follow-up pruning of a complete tree; 4 TH = topping and one-sided hedging, alternating annually between sides (East/West); 5 TH/M = alternating annually between TH and M; 6 THF = TH plus manual follow-up pruning of the non-hedged side; 7 THH = topping and bilateral hedging; 8 THH/M = alternating annually between THH and M; 9 HH = bilateral hedging; 10 HH/M = alternating annually between HH and M; 11 H = one-sided hedging, alternating annually between sides (East/West); 12 H/M = alternating annually between H and M.
Dates of the pruning treatments and the BBCH (Biologische Bundesanstalt, Bundessortenamt and Chemical industry) citrus stage are shown in Table 2.

Pruning Machines
Manual pruning was performed with handsaws and manual or electric pruning shears (Electrocoup F3015, Infaco s.a.s., Cahuzac-sur-Vère, France), depending on the diameter of the branches.
Mechanical pruning (hedging and topping) was performed using a pruner of a single arm provided with cutting discs. In the 2017-2019 trials, the pruner used was the model PF-605 S-R with five discs, and in 2020, the model PF-606 S-R pruner, with six discs, was used (  In the second year of the experiment, strategies 4, 7, 9, and 11 were divided into strategies 5, 8, 10, and 12, respectively, of annual alternation between mechanical pruning and manual pruning (Figure 1). A minimum of 10 consecutive trees were used for each strategy and repetition, except for the control, in which each repetition consisted of five consecutive trees. All the evaluations were carried out in a representative tree that was randomly selected from the experimental unit (except the extreme trees) in each strategy and repetition, except in the control, in which the central tree was used for this purpose.
Dates of the pruning treatments and the BBCH (Biologische Bundesanstalt, Bundessortenamt and Chemical industry) citrus stage are shown in Table 2.

Pruning Machines
Manual pruning was performed with handsaws and manual or electric pruning shears (Electrocoup F3015, Infaco s.a.s., Cahuzac-sur-Vère, France), depending on the diameter of the branches.
Mechanical pruning (hedging and topping) was performed using a pruner of a single arm provided with cutting discs. In the 2017-2019 trials, the pruner used was the model PF-605 S-R with five discs, and in 2020, the model PF-606 S-R pruner, with six discs, was used ( Figure 2

Characterization of Pruned Biomass and Trees
The sizes of the trees were measured before pruning ( diameter along the row (Ø L ), and the canopy diameter across the row (Ø C ) ( Figure 3). The canopy height was calculated by subtracting the height from the ground to the skirts from the total tree height. All measurements were taken without considering unusual extreme shoots.

Characterization of Pruned Biomass and Trees
The sizes of the trees were measured before pruning ( Table 2) to calculate the canopy volume. The geometrical parameters measured were the canopy height (H), the canopy diameter along the row (ØL), and the canopy diameter across the row (ØC) (Figure 3). The canopy height was calculated by subtracting the height from the ground to the skirts from the total tree height. All measurements were taken without considering unusual extreme shoots.
where H is the canopy height (m), ØC is the canopy across-row diameter (m), and ØL is the canopy along-row diameter (m). Furthermore, fresh pruned biomass was characterized after pruning, to determine the intensity of each pruning strategy for each year, through the following parameters: - The length (cm) and the diameter (mm) of the branch in the cutting area were measured in a minimum of 100 branches for the topping, hedging, and manual follow-up cuttings, and in a minimum of 50 branches for the manual cutting. For this, the branches were chosen from pruning strategies in which there was no annual alternation between manual and mechanical pruning. - The quantity of fresh biomass removed with each pruning strategy was measured (kg•tree −1 ). For this, all the pruned branches in each type of cutting of the evaluation The volume of vegetation was calculated considering the citrus canopy as an ellipsoid, as shown in Equation (1): where H is the canopy height (m), Ø C is the canopy across-row diameter (m), and Ø L is the canopy along-row diameter (m). Furthermore, fresh pruned biomass was characterized after pruning, to determine the intensity of each pruning strategy for each year, through the following parameters: - The length (cm) and the diameter (mm) of the branch in the cutting area were measured in a minimum of 100 branches for the topping, hedging, and manual follow-up cuttings, and in a minimum of 50 branches for the manual cutting. For this, the branches were chosen from pruning strategies in which there was no annual alternation between manual and mechanical pruning. - The quantity of fresh biomass removed with each pruning strategy was measured (kg·tree −1 ). For this, all the pruned branches in each type of cutting of the evaluation tree were weighed using a digital dynamometer (Advanced Force Gauge 500 N, Mecmesin Ltd., England, UK). Afterwards, the total fresh biomass of each sample tree was calculated by adding the weights obtained in the different types of cuts included in each strategy. - The moisture content (%) of the fresh biomass was measured in the laboratory, based on the thermogravimetric analysis method, for the first year of the trial. For this, five random samples of fresh pruned biomass of each type of cutting (topping, hedging, manual follow-up, and manual pruning) were selected. Each sample was divided into leaves and wood, and dried in an oven (Dry-Big 2003741, J.P. Selecta S.A., Abrera, Barcelona, Spain) at 65 • C for the time necessary to reach a constant weight. In this way, the moisture content (%) in the biomass of leaves and wood, respectively, was calculated with Equation (2): where M is the moisture content (%), W is the weight of the wet sample (g), and D is the weight of the dry sample (g).

Pruning Effect on Citrus Production: Yield and Fruit Size
To evaluate the effect of the pruning strategy on citrus production each year, the yield and fruit size were determined at harvest time. To obtain the yield (kg·tree −1 ), all of the fruits of the trees under study (one of each strategy and repetition) were collected and weighed with an electronic scale (model PCE-EP 150P1, 150 kg capacity, PCE Ibérica S.L., Tobarra, Albacete, Spain).
Fruit size was determined by measuring the outer equatorial diameter (mm) of 50 fruits that were randomly selected from each tree under study using an automatic computer vision inspection system developed by IVIA [38,39]. This vision system consists of one industrial colour camera (Manta G-125C, Allied Vision Technologies GmbH, Germany) located on a line of fruit transport rollers, and a lighting system based on LED strips. The system acquired 4 images per fruit, with a resolution of 1292 × 964 pixels per image, and the images were processed on an industrial computer (NUVO 5006E-POE, Neousys Technology Inc., Taiwan) that ran the image analysis algorithms and obtain the fruit size results. With these data, the average fruit size for each pruning strategy and year was obtained.

Pruning Working Capacity and Costs
The theoretical working capacity (TWC; trees·h −1 ·operator −1 ) of each pruning strategy was determined based on the productive time spent to prune a tree (min·tree −1 ·operator −1 ). This time was measured in a different way depending on whether mechanical or manual pruning was used. In the case of mechanical pruning, the tractor forward speed was determined in several rows, for each year, for the topping and hedging cuts conducted in one side. For this, the time taken by the tractor to prune a row of known distance and number of trees was measured with a chronometer. From these values, the average productive time per tree was calculated. In the case of manual pruning, both for the fully manual strategy and the follow-up approaches included in the TFF (whole tree) and THF (half tree) strategies, the productive pruning time was determined by timing with a chronometer the pruning of a single tree with a handsaw and manual pruning shears. A minimum of 13 randomly chosen trees were measured each year. The productive time spent per whole tree was multiplied by the number of operators involved in the operation (usually one or two specialised workers) to calculate the time needed by a single worker (min·tree −1 ·operator −1 ).
Once the four-year average productive time (min·tree −1 ·operator −1 ) spent per type of cutting (one-sided topping, one-sided hedging, and whole-tree manual) was obtained, the productive time of each strategy was estimated by adding all of the productive times used in the different types of cuttings included in each strategy. Based on these values, the average number of trees pruned per hour for each operator (TWC), in each strategy, was calculated.
In addition to the TWC value, the work performance (h·operator·ha −1 ) for each pruning strategy was calculated using Equation (3): where N is the tree density (trees·ha −1 ; in this orchard: 556 trees·ha −1 ) and TWC is the theoretical working capacity (trees·h −1 ·operator −1 ). Finally, the average cost of the pruning operation (€·ha −1 ), without considering dead time, for each strategy was determined with Equation (4): where C is the hourly pruning cost (€·h −1 ). Data of hourly pruning service prices were obtained directly from the orchard technicians who contracted them, and were 7.9 €·h −1 for manual pruning and 39 €·h −1 for mechanical pruning. The cost structure of C includes the labour cost and the equipment cost (manual pruning = labour + manual tools; mechanical pruning = tractor driver + tractor + pruner).

Economic Profit
The economic profit was evaluated by calculating the confidence interval of the net value + k (€·ha −1 ) for each strategy and year, and for the four-years average, where k is the sum of all the production costs except the pruning ones and was considered the same for all pruning strategies. For that, the net value + k per hectare was calculated with Equation (5): where Gross value (€·ha −1 ) was calculated according to Equation (6): where Yield (kg·tree −1 ) is the kilograms of mandarins produced per tree each year; Mandarins price refers to the average citrus prices at the harvest dates of each year (Table 3) obtained from the official citrus price table [40]; N is the tree density of the orchard (556 trees·ha −1 ). Afterwards, the confidence interval of the net value + k per strategy and year was calculated using Equation (7): where A net value + k is the average of the net value (€·ha −1 ), z is the confidence level value (95%), s is the standard deviation, and n is the number of elements in the sample.

Statistical Methods
For each year and for the 4-year overall, the effects of the pruning strategies on the canopy volume, the pruned fresh biomass, the yield, the fruit size, and the net value + k, and the effects of the type of cutting used on the dimensions (diameter and length) of the cut branches, the percentages of fresh leaf biomass and wood biomass, and the contents of leaf and wood biomass moisture, were studied through a one-way analysis of variance (ANOVA). The assumption of homoscedasticity was verified with the Levene test, and normality was evaluated with the normal probability plot of the residuals for each ANOVA. When these assumptions were not fulfilled, data transformation was applied.
If homoscedasticity or normality assumptions were still not fulfilled, a nonparametric Kruskal-Wallis test was applied.
Fisher's least significant difference (LSD) test was used for comparisons of means in the ANOVA analysis. In the Kruskal-Wallis tests, the significance of differences in the medians was determined through the non-overlapping of the median notches in the box plots. A 95% confidence level was considered in all tests.
A descriptive analysis was performed to study the productive pruning time and the pruning costs for each pruning strategy. For the net value, the confidence interval was calculated as explained in Equation (7).

Characterization of Pruned Biomass and Canopy Dimensions
The dimensions of the canopy before pruning (dates in Table 2) for each pruning strategy, and for each year, are shown in Table 4. In 2017, the canopy volume of the trees before starting the assay ranged between 15.90 and 24.23 m 3 ·tree −1 ; in 2018, the canopy volume of the pruned trees (strategies 2 to 12) was between 21.76 and 25.18 m 3 ·tree −1 ; in 2019, the canopy volume of the pruned trees was between 17.63 and 23.75 m 3 ·tree −1 ; and in 2020, this value was between 20.46 and 25.14 m 3 ·tree −1 . The canopy volumes of the unpruned trees (control), from 2017 to 2020, were 18.49, 24.49, 24.19, and 24.99 m 3 ·tree −1 , respectively (Table 4).  The average diameter and length of the cut branches by type of cutting and year are shown in Table 5. The branches that were pruned mechanically (topping and hedging) had a smaller size than the branches that were subjected to manual pruning or manual followup pruning. For the four-years overall, in the hedging cut, the branches were 6.78 mm in diameter and 43.72 cm in length; in the topping cut, they were 7.24 mm in diameter and 64.15 cm in length. The branches that were pruned manually had a diameter of 20.06 mm. These branches had a length of 133.84 cm in exclusively manual pruning (strategy 2), and a length of 112.94 cm in manual follow-up pruning (Table 5). These values evidenced that mechanical pruning removed the younger external branches, while manual pruning removed older and larger branches, and aimed to open windows in the canopy in order to allow aeration, incidence of solar radiation, penetration of pesticides, etc.
The pruned fresh biomass for each pruning strategy and year is shown in Table 6. On average for the four years, the amount of biomass pruned with exclusively mechanical pruning strategies (strategies 4, 7, 9, and 11; between 3.18 and 11.48 kg·tree −1 ) was lower than with pruning strategies that included any manual cutting (between 16.25 and 26.96 kg·tree −1 ). Furthermore, within the mechanical pruning strategies, more than twice the biomass was removed in the strategies that included topping (4 and 7) than in the strategies without topping (9 and 11). Moreover, it was observed that with the manual pruning strategies (strategy 2 (all years), and strategies 5, 8, 10, and 12 (2018 and 2020)), a higher amount of biomass was removed than with the mechanical pruning strategy (topping) with manual follow-up in the complete tree (strategy 3) every year. While with strategy 2 (M), 26.96 kg·tree −1 was pruned on average in the four years, with strategy 3 (TFF), this quantity dropped to 21.72 kg·tree −1 . For the strategies that alternated between manual and mechanical pruning, it was observed that in the years when manual pruning was used (2018 and 2020), as much or more biomass was removed than was the case for wholly manual strategy (strategy 2), which could tentatively be explained as a compensatory response to the lighter mechanical pruning of the previous year. Table 5. Diameter of the cutting area (Ø branches , mm) and length (L branches , cm) of cut branches by type of cutting, for each year and for the four years overall. Branches were chosen from pruning strategies in which there was no annual alternation between manual and mechanical pruning. Data are expressed as averages, with standard errors in parentheses. Different letters between types of cutting for Ø branches and L branches indicate significant differences (p < 0.05) according to the median notched box plot for each year 1 and according to Fisher's LSD test for the 4-year overall 2 .

Year
Type of Cutting Ø branches (mm) L branches (cm) The percentages of fresh leaves and wood in the pruned biomass samples for the different types of cut (hedging, topping, manual follow-up, and manual), as well as their moisture contents, are shown in Table 7. Globally, without differentiating by type of cut, the leaf biomass represented 56.47% of the total fresh biomass, and its moisture content was 63.04%. The percentage of wood was 43.53%, with a moisture content of 45.12%. However, it was observed that the moisture content in the manually pruned wood was lower in comparison with that of the mechanically pruned wood. This could be explained in terms of type of branches; the branches that were cut by mechanical pruning were younger and non-lignified, while the branches that were cut manually were mature and lignified. It could also have been due to the fact that during manual pruning, the operators cut, among others branches, the dried ones that they found in the internal part of the tree. Table 7. Percentage of fresh weight of biomass of leaves and wood of the cut branches and their percentage of moisture, by type of cutting and globally. Data are expressed as averages, with standard errors in parentheses. Different letters in a column indicate significant differences between types of cutting for each parameter 1 according to Fisher's LSD test (p < 0.05).
In the third year (2019), significant yield reductions were observed regardless of strategy, and no significant differences were found between strategies (F = 1.47; df = 11, 58; p = 0.175).   In the fourth year of the trial (2020), statistically significant differences were again observed between strategies (F = 23.13; df = 11, 59; p < 0.001). A significantly greater yield was achieved with pruning strategy 5 (TH/M) (127.67 kg·tree −1 ) than with the others. Trees that were pruned manually (strategy 2 and strategies 5, 8, 10, and 12), with manual follow-up on both sides (strategy 3: TFF), as well as those that were not pruned (control) (between 80.32 and 127.67 kg·tree −1 ), had significantly higher yields than trees that were pruned exclusively mechanically (strategies 4, 7, 9, and 11) or mechanically with manual follow-up only on one side (strategy 6: THF) (between 33.51 and 56.77 kg·tree −1 ). In addition, it was observed that the mechanical pruning strategies that included hedging on both sides of the tree had significantly lower yields than the strategies involving hedging only on one side (20.14 kg·tree −1 less in HH than in H, and 21.70 kg·tree −1 less in THH than in TH) and the strategy that included manual follow-up on both sides (TFF) had significantly higher yield (26,72 kg·tree −1 more) than the strategy of manual follow-up on one side and mechanical hedging cutting on the other side (THF). On the other hand, it was observed that in mechanical pruning strategies, topping did not affect yield (there were no significant differences between THH and HH, nor between TH and H). However, when comparing full manual pruning (M) strategy with that of manual pruning plus topping (TFF), a significant reduction in yield was observed (28.70 kg·tree −1 more in M than in TFF).
Finally, when the average for the four years was considered, non-significant differences between strategies were found (F = 0.45; df = 11, 47; p = 0.9188). However, on average, the no pruning, manual pruning, or mechanical pruning alternated annually with manual pruning strategies had higher yield than that of exclusively mechanical pruning. Furthermore, the mechanical pruning strategy with the highest number of cuts (THH) had the lowest yield (62.93 kg·tree −1 on average) (Figure 4).
With regard to the fruit size, there were also significant differences between years (F = 79.49; df = 3, 43; p < 0.001), with the fruits produced in 2018, when the yield was higher, being significantly smaller (50.89 mm) than the fruits of 2020 (57.23 mm), and in turn, significantly smaller than those of 2017 and 2019 (60.20 and 60.92 mm, respectively) ( Figure 5).
Comparing the strategies by year, in the first year of the experiment (2017), there were significant differences in fruit size between strategies (F = 2.82; df = 7, 39; p = 0.021). It was observed that with mechanical pruning strategies that involved manual follow-up pruning, on one (THF) or on both sides (TFF) of the tree, the fruits had a significantly larger diameter than with mechanical pruning only on one side of the tree (H) or with the trees not being pruned (control).
In 2019, with the hedging cutting strategies (HH and H), the fruits had a significantly smaller diameter than with manual pruning (M), mechanical pruning with manual followup (TFF and THF), or mechanical pruning with topping and hedging on both sides (THH) (F = 3.28; df = 11, 58; p = 0.0021).
In 2020, significant differences were also found between pruning strategies (F = 9.94; df = 11, 59; p < 0.001). During this year, significantly larger fruits were obtained with the TH and TH/M strategies than with the HH, HH/M, H, and H/M strategies and with no pruning (control).
When considering the average fruit size of the four years, no significant differences were found between strategies (F = 0.33; df = 11, 47; p = 0.9727). However, it was observed that the HH and H strategies, as well as the control (without pruning), were associated with smaller fruit sizes (between 55.09 and 55.54 mm) compared to the rest of the strategies, for which similar fruit sizes were observed (between 56.80 and 58.73 mm) ( Figure 5).
With regard to the fruit size, there were also significant differences between years (F = 79.49; df = 3, 43; p < 0.001), with the fruits produced in 2018, when the yield was higher, being significantly smaller (50.89 mm) than the fruits of 2020 (57.23 mm), and in turn, significantly smaller than those of 2017 and 2019 (60.20 and 60.92 mm, respectively) ( Figure 5).

Figure 5.
Fruit size (equatorial diameter, mm) (average value-indicated numerically within the column-with standard error bar) by pruning strategy (1-12) and the average for all strategies, for each year and for the four years overall. Different lowercase letters above the bars for each year and four-year overall indicate significant differences (LSD test, p < 0.05). Different capital letters above the bars of the averages of all strategies for each year indicate significant differences (LSD test, p < 0.05).

Figure 5.
Fruit size (equatorial diameter, mm) (average value-indicated numerically within the column-with standard error bar) by pruning strategy (1-12) and the average for all strategies, for each year and for the four years overall. Different lowercase letters above the bars for each year and four-year overall indicate significant differences (LSD test, p < 0.05). Different capital letters above the bars of the averages of all strategies for each year indicate significant differences (LSD test, p < 0.05).

Pruning Working Capacity and Costs
The productive time to perform the topping cut for one side of the tree was 0.10 min·tree −1 ·operator −1 with a tractor ground speed of 1.85 km·h −1 . For the onesided hedging cut, this was 0.08 min·tree −1 ·operator −1 with a tractor ground speed of 2.30 km·h −1 . Furthermore, the time to manually prune a whole tree was 5.48 min·tree −1 · operator −1 , on average, for the four years. The estimated productive time for each strategy is shown in Table 8. Table 8. Four-year average estimated productive pruning time (min·tree −1 ·operator −1 ), pruning theoretical working capacity (TWC, tree·h −1 ·operator −1 ), work performance (h·operator·ha −1 ), and pruning cost (€·ha −1 ) for each pruning strategy. Data are expressed as averages.
The pruning cost of the exclusively mechanical strategies ranged between 28.75 and 130.01 €·ha −1 depending on the number of pruner passes (one pass in strategy 11, two passes in strategy 9, three passes in strategy 4, and four passes in strategy 7). This represented a reduction of between 67.62% and 92.84% with respect to the costs of the pure manual pruning approach (strategy 2). In the strategy involving mechanical pruning and manual follow-up pruning on one side only (strategy 6), the cost was reduced by around 24.78% with respect to the pure manual pruning strategy. The strategy comprising mechanical topping and manual follow-up of the whole tree (strategy 3) had a cost that was 18.06% higher than the pure manual pruning strategy. Pruning strategies that alternated manual with mechanical pruning annually had costs that were between 33.81% and 46.42% lower than the strategy of manual pruning for four consecutive years (Table 8).

Economic Profit
The confidence interval of the net value + k (€·ha −1 ) per pruning strategy and year is shown in Table 9. In the first year, pruning strategies involving topping (TFF, TH, THF, and THH) had the lowest net values + k, between 7838 and 13,061 €·ha −1 . In the second and third years, in general the intervals of all strategies overlapped and, therefore, no differences were found in the net values + k. In the fourth year, the mechanical pruning strategies on both sides, with and without topping (THH and HH), had the lowest net values + k (between 3575 and 6819 €·ha −1 ). With the pure manual pruning strategy (M; between 15,463 and 17,921 €·ha −1 ), higher net values + k were achieved than with the exclusively mechanical or mechanical pruning strategies involving manual follow-up (TFF, TH, THF, THH, HH, and H; between 3575 and 13,416 €·ha −1 ). Furthermore, with the alternate mechanical/manual pruning strategies, the net values + k were much higher than with the corresponding mechanical pruning strategies. When considering the averages of the four years, as observed in the graph in Table 9, there were no significant differences in the net values + k between the different pruning strategies (F = 0.86; df = 11, 47; p = 0.5861). However, on average, with exclusively manual pruning and mechanical/manual pruning strategies and with no pruning, higher net values + k were achieved than with exclusively mechanical or mechanical pruning strategies involving manual follow-up. Table 9. Net value + k (€·ha −1 ) for each pruning strategy, for each year and for the four years overall. Data are expressed as 95% confidence interval. The graphic represents the data for the 4-year overall (average value with 95% confidence interval bar). Different letters indicate significant differences between pruning strategies (LSD test, p < 0.05).

Discussion
In this study, the pruning strategies can be grouped into two large groups according to the amount of biomass removed, namely the strategies that included some type of manual pruning (2, 3, and 6; and 5, 8, 10, and 12), with the largest amount of biomass pruned (20.99 kg•tree −1 on average) due to bigger branches removed both in length and diameter, and the ones with only mechanical pruning including topping (strategies 4 and 7) or without topping (strategies 9 and 11), with a biomass removed of 10.80 kg•tree −1 and 3.62 kg•tree −1 on average, respectively. These values are in the same order of magnitude as the ones obtained in pruned citrus by Velázquez and Fernández [36] and Velázquez-Martí et al. [41] and the moisture content was also similar, at around 30-45%, but without distinguishing between manual and mechanical pruning; furthermore, in this work, it was observed that the wood moisture content was higher (57.59% on average) with mechanical pruning than with manual pruning (32.64%), which can be explained because the removed branches were external, younger, and greener.
A poor relationship between biomass removed and yield was observed for the four years overall, with a correlation coefficient of only R 2 = 0.3208 (linear regression) ( Figure  A1 of Appendix A). Furthermore, this tendency was similar when analysing the data for each year, except in 2020, with a correlation coefficient of 0.7491. In 2017, with large significant differences in yield between pruning strategies, no pruned-biomass-yield relationship was found (R 2 = 0.0123), and in 2018 and 2019, with no significant differences in yield between pruning strategies, the correlation coefficient was also low (R 2 = 0.0229 in 2018 and R 2 = 0.4604 in 2019).
However, the pruning strategy can affect the yield. In fact, in the first year of assay (2017), the strategies involving topping led to significantly reduced yields in comparison with the rest of the strategies, in a similar way to that observed by Wheaton et al. [32] and Kallsen [33]. However, in this work, it was also observed that after the yield reduction for the first year in the strategies involving topping, in the second and third years, the yield was recovered and this behaviour was also found in 'Ponkan' tangerines, according to Mendonça et al. [34]. Anyway, the effects of pruning strategies seem to be noticed after some years of experience. In the fourth year (2020), trees with mechanical pruning for four consecutive years had lower yield although the differences were not significant. Meanwhile, strategies that annually alternated mechanical pruning with manual pruning had a higher yield, similar to those the manual and control strategies. As observed by Martin-Gorriz et al. [22] in 'Fortune' mandarins, the annual or biannual alternation of manual pruning with mechanical pruning obtained an accumulated yield close to that of continued manual pruning and a higher yield than that obtained with mechanical pruning only. In 'Salustiana' oranges, there was a reduction of 17% in the yield of the trees that were pruned mechanically compared to the hand-pruned trees considering the average yield of the four years [29,42]. Continued mechanical pruning produces zig-zag branches

Discussion
In this study, the pruning strategies can be grouped into two large groups according to the amount of biomass removed, namely the strategies that included some type of manual pruning (2, 3, and 6; and 5, 8, 10, and 12), with the largest amount of biomass pruned (20.99 kg·tree −1 on average) due to bigger branches removed both in length and diameter, and the ones with only mechanical pruning including topping (strategies 4 and 7) or without topping (strategies 9 and 11), with a biomass removed of 10.80 kg·tree −1 and 3.62 kg·tree −1 on average, respectively. These values are in the same order of magnitude as the ones obtained in pruned citrus by Velázquez and Fernández [36] and Velázquez-Martí et al. [41] and the moisture content was also similar, at around 30-45%, but without distinguishing between manual and mechanical pruning; furthermore, in this work, it was observed that the wood moisture content was higher (57.59% on average) with mechanical pruning than with manual pruning (32.64%), which can be explained because the removed branches were external, younger, and greener.
A poor relationship between biomass removed and yield was observed for the four years overall, with a correlation coefficient of only R 2 = 0.3208 (linear regression) ( Figure A1 of Appendix A). Furthermore, this tendency was similar when analysing the data for each year, except in 2020, with a correlation coefficient of 0.7491. In 2017, with large significant differences in yield between pruning strategies, no pruned-biomass-yield relationship was found (R 2 = 0.0123), and in 2018 and 2019, with no significant differences in yield between pruning strategies, the correlation coefficient was also low (R 2 = 0.0229 in 2018 and R 2 = 0.4604 in 2019).
However, the pruning strategy can affect the yield. In fact, in the first year of assay (2017), the strategies involving topping led to significantly reduced yields in comparison with the rest of the strategies, in a similar way to that observed by Wheaton et al. [32] and Kallsen [33]. However, in this work, it was also observed that after the yield reduction for the first year in the strategies involving topping, in the second and third years, the yield was recovered and this behaviour was also found in 'Ponkan' tangerines, according to Mendonça et al. [34]. Anyway, the effects of pruning strategies seem to be noticed after some years of experience. In the fourth year (2020), trees with mechanical pruning for four consecutive years had lower yield although the differences were not significant. Meanwhile, strategies that annually alternated mechanical pruning with manual pruning had a higher yield, similar to those the manual and control strategies. As observed by Martin-Gorriz et al. [22] in 'Fortune' mandarins, the annual or biannual alternation of manual pruning with mechanical pruning obtained an accumulated yield close to that of continued manual pruning and a higher yield than that obtained with mechanical pruning only. In 'Salustiana' oranges, there was a reduction of 17% in the yield of the trees that were pruned mechanically compared to the hand-pruned trees considering the average yield of the four years [29,42]. Continued mechanical pruning produces zig-zag branches that become unproductive over time, and more severe rejuvenation pruning is needed after a few years [16]. Alternating mechanical pruning with selective manual pruning avoids the accumulation of unproductive branches in the canopy, and avoids having to undertake a more severe rejuvenation pruning, with the consequent yield reduction in production that it would entail for that year.
In contrast to these results, in 'Navel Foyos' oranges, no significant differences in yield and fruit size between manual and mechanical pruning strategies were found for the three years [11] and this was also the case with those obtained for 'Washington Navel' oranges, where the average yields for the four years in the manual pruning strategy, in the mechanical pruning strategy, and in the strategy involving a mix of mechanical and manual pruning were similar [29,42]. Moreover, in 'Nadorcott', a mandarin variety with a marked biennial bearing, which means that trees will have great fruiting one year and very little flowering the following year, Mesejo et al. [16] demonstrated that cutting the flowering shoots annually with a pruning strategy based on topping and hedging on both sides and with different severity levels of the topping (greater or lesser depth), between 7.2% and 24.5% greater yields, respectively, were obtained, compared to non-pruned trees, after four years of experimentation.
Agustí et al. [43] indicated that the size of the fruit was inversely related to the number of fruits per tree, although this relationship was only clear under conditions of low productivity. If the yield was high, the correlation between both variables would be low. In this work, it can be observed in Figure A2 of Appendix A that the sizes of the fruits (equatorial diameter) decreased as the yield increased, although the coefficient of determination was low (R 2 = 0.2812, linear regression), in a similar way to that which was observed by Martin-Gorriz et al. [22]. In addition, when considering the different years, it was observed that 2018, the year with the highest yields, had the lowest fruit sizes, and 2019, the year with the lowest yields, had the highest fruit sizes. However, when considering the pruning strategies, this trend could not be observed. In some cases, strategies with lower production also had smaller fruits and vice versa. Ahmad et al. [44] found no significant differences in fruit size between pruned and non-pruned trees in the 'Kinnow' mandarin hybrid. Bevington and Bacon [13] did not find a consistent relationship between hedging and fruit size, although on average, the unpruned trees had smaller fruits. However, Morales et al. [19] observed that pruning (topping, with or without skirting) increased the percentage of large fruits and reduced the number of small fruits, as well as decreasing yields, of 'Orlando' tangelo.
Regarding the productive pruning time, it was observed that in 'Clemenules', the pruning workers took the same time in the full manual follow-up pruning as in the manual pruning, since the operators cut similar numbers of internal branches, although trees were previously pruned mechanically in the external part. In fact, the diameter in the cutting area of the branches pruned in manual and manual follow-up pruning strategies were the same, on average, for the four years (20.06 mm). Meanwhile, in 'Fortune' mandarins, Martin-Gorriz et al. [22] found that manual follow-up pruning after mechanical pruning required between 13% and 15% less time than manual pruning. Anyway, because this time depends directly on the subjective vision of the pruner operator regarding which internal branches should be cut, a high reduction in time was not expected. On the contrary, mechanical pruning reduced the time required for pruning, allowing a larger orchard area to be pruned in less time. Furthermore, although the prices per hour for mechanical pruning were higher than for manual pruning (39 €·h −1 versus 7.9 €·h −1 ), the mechanical pruning costs per hectare were between 3 (strategy THH) and 14 (H) times lower, depending on the number of pruner passes, than those of manual pruning (M).
The net value + k takes into consideration the other variables studied: yield, working capacity, and costs, with yield being the one with the highest weight. For this reason, the higher the yield, the greater the net value + k, and vice versa. For the four years overall, the highest net value + k was achieved with the control (trees not pruned), ranging between 10,948 and 20,305 €·ha −1 . Nonetheless, this should not be an option to consider as a normal practice since pruning is necessary to control the overgrowth of the trees, and to guarantee their correct development, as well as their health and future yield [12,[14][15][16]. In the long term, it is expected that approaches involving no pruning may reduce production and, consequently, net value would decrease. The strategies with higher net values + k were manual pruning (M) and alternate annual mechanical/manual pruning (TH/M; THH/M; HH/M; H/M), and the ones with lower net values + k were mechanical pruning (TH; THH; HH; H) and mechanical pruning with manual follow-up (TFF; THF).

Conclusions
After four years of pruning experimentation in 'Clemenules', it may be concluded that alternating mechanical pruning in one year with manual pruning in the next year may be a better option than using exclusively mechanical pruning or mechanical pruning with manual follow-up every year because the yield obtained with the former approach was higher and similar to those of the manual and control strategies. Furthermore, with this system, the dry branches would be removed from the interior of the tree, internal aeration would be improved, and the problems of non-selective mechanical pruning continued over time would be avoided. At the same time, the dimensions of the trees would be more regular, facilitating the movement of tractors between the tree rows. In addition, this strategy would mean a significant cost reduction compared to traditional manual pruning. For all these reasons, this system should be boosted in farm practices through the diffusion of these results in agricultural extension channels.
Among the strategies that alternated mechanical and manual pruning yearly, the strategy that combined mechanical pruning in topping and one-sided hedging on with completely manual pruning in the following year (TH/M) was the one with the highest yield (91 kg·tree −1 on average for the four years). This strategy, by including the topping cut, makes it possible to control the growth of the trees, since, if the growth in height of citrus trees is not limited, the harvesting tasks, which in the Valencian region are usually carried out manually, and the spraying of plant protection products will become difficult.
The results obtained in this research were focussed on the 'Clemenules' variety and, as has been observed in the literature, were different to the results obtained for other citrus varieties/species; this is a limitation of this work and more research is required on other citrus varieties to observe their behaviour after mechanical pruning. On the other hand, the effect of pruning may affect other aspects of production such as the presence of pests/diseases and its control. This, together the influence of climate change in all aspects of citrus production, are future challenges for this research area. Figure A1. Pruned fresh biomass (kg•tree −1 ) versus total yield (kg•tree −1 ) for each pruning strategy (strategies 2-12), in each year under study (2017-2020), and for the 4-year overall. The model equations that describe the relationships between pruned fresh biomass and yield are indicated. Figure A2. Average fruit size (equatorial diameter, mm) versus total yield (kg•tree −1 ) for each tree under study during the four years of trials (N = 219). The model equation that describes the relationship between fruit size and yield is indicated.