Climate variables effect on fruiting pattern of Kinnow mandarin (Citrus nobilis Lour × C. deliciosa Tenora) grown at different agro-climatic regions

Kinnow orchards grown in different agro-ecological regions of Punjab, Pakistan, namely Sargodha, Toba Tek Singh (TTS) and Vehari districts, were selected to assess the effect of climate variables on fruit-bearing patterns. Experiment was laid out in RCBD while selecting identical features Kinnow plants and labeled twigs at analogous canopy positions in all three sites. Temperature was reported higher in TTS and Vehari areas, while relative humidity in Sargodha accounted for different levels of agrometeorological indices by computing more variations in warm districts. Climate variables influenced fruit-bearing habits and vegetative growth trend in all three flushes while recording heavy fruit-bearing plants during on-year and light fruit-bearing in off-year at Vehari. Similarly, three vegetative flushes were recorded unevenly in all three sites due to different fruit-bearing patterns induced by climate variables. Harvesting pattern of orchards began earlier in Sargodha, where maximum orchards were harvested before new flowering to add evenness to fruiting habits during on & off-years. In warm conditions, fruit ripening arrived in the peak of winter and mostly domestic market-driven harvesting resulted in late start of fruit picking with more erratic fruit-bearing habits. Both physiological and pathological fruit drops have been significantly affected by climate variables with a higher degree of physiological drop in warm regions and pathological effects in the humid conditions of Sargodha on heavy fruit-bearing plants. Fruit yield and grading quality were also affected in both seasons by showing more asymmetrical trend in yield and fruit grading in warm areas of TTS and Vehari due to an irregular fruiting pattern compared to Sargodha. From now on, the climate variables of the three sites directly influenced the fruiting patterns, vegetative flushes, fruit drops, yields and grades of Kinnow mandarin.


Selection of orchards.
Kinnow orchards were selected in block form with similar characteristics of plants like age, health, vigor, planting system/ geometry (square), density (250-260 plants/ha) and grafted on Rough lemon(Citrus Jambhiri Lush.) rootstock 10 by tagging branches/ twigs to reflect uniform canopy positions of the individual plant 11 . Basic soil properties were analyzed in three sites with organic matters ranges (0.75-0.80%), available phosphorous (6.0-6.5 mg kg −1 ), available potassium (230-260 mg kg −1 ) and loamy structure soil 11 . In each orchard, uniform dose of fertilizers i.e., Nitrogen (1000 g), phosphorous and potash (500 g each) were applied.
Plant material identification and not deposit in herbarium. The plants of Kinnow mandarin were selected by researcher (R.N) as being used in research. Plants having age 12-15 years old were not deposited in any public herbarium as it is not a new species, with no need to deposits as from orchard not uprooted.
Climatic/weather data. Weather data were collected from the Pakistan Meteorological Department (PMD), Islamabad, of three experimental sites and the office of the Deputy Director Agriculture (Extension), Vehari, used in computation of climate variables. Temperature data are shown in Fig. 1 and rainfall as well as relative humidity in Fig. 2. Total annual precipitation (511, 349 and 144 mm), average annual temperature (23.65, 25.19 and 27.11 °C) and average annual relative humidity (66.8, 63.20 and 55.6 percent) were reported in the districts of Sargodha, TTS and Vehari during 2017 and 2018 respectively. Climate variables computation. Climate variables like relative temperature disparity (RTD), phototemperature (T p ), nyctotemperature (T n ), relative humidity disparity (RHD), photo relative humidity (RH p ), nycto relative humidity (RH n ), vapor pressure deficit (VPD), photo vapor pressure deficit (VPD p ) and nycto vapor pressure deficit (VPD n ) were calculated from weather data using computation of [30][31][32]59 . Relatively humidity recorded at midnight (0000 UTC) and midday (1200 UTC) was used while computing climate variables. Thermal energy use efficiency. Thermal energy use efficiency was computed through Kinnow plant yield (kg ha −1 ) basis. Accumulated thermal indices used in computation of thermal energy use efficiency were calculated from fruit-set till arrival of maturity in fruit 2,5,[9][10][11] .

Results and discussion
A significant difference in fruit bearing habits was seen during the on-& off-year period, which is explained as below.
Climate variables of experimental sites. Relative temperature disparity (RTD) is wider in winter season and narrows down in summer with a higher trend in Sargodha and lower in Vehari. Phototemperature recorded from daylight maximum temperature and nyctotemperature to depict night duration lower temperature and both showing an increasing level in warm conditions of Vehari during summer and lower at Sargodha in winter months. However, higher phototemperature was reported in Vehari (40.1 °C) in June 2018 and lower nyctotemperature was observed in Sargodha (8.0 °C) in January 2018. Relative humidity disparity (RHD) is the difference between the relative humidity reported at midnight (0000 UTC) and midday (1200 UTC) with widening differences between April and May in both years in Sargodha, followed by TTS and narrowed down in August-September in all three districts.  10 . Sargodha district is to the north with a higher elevation, while Vehari is to the south with a low elevation. The climatic conditions are therefore different on a monthly as well as an annual basis in all three sites by having a direct effect on plant physiology 7 and indirectly on fruit bearing habits and maturity arrival on Kinnow fruit 3 www.nature.com/scientificreports/ and also depict short and long weather conditions of a particular region 34 . Temperature has a cumulative effect on growth, yield and productivity span of plants 35,36 as well as used in computation of climatic indices which are further utilized in management process 61,62 . Phototemperature is climate variable to depict a daytime temperature which also termed as active or positive temperature just above base temperature of particular crop in most days of the year 62,63 and has been widely used in the measurement of different heat units 64,65 . Nyctotemperature derived from minimum temperature is referred as inactive temperature, which usually fall below the threshold temperature of sub-tropical crops like citrus in all winter, spring and autumn months except for a few summer months in warm regions 9,66 and also decide on the accretion in heat units used in plant phenological studies as well as future strategies for controlling insect-pests and diseases 67 . Temperature directly affects plant tissues and organs with death in extreme conditions 68 while its intermediate level affects physiological processes, including reproductive organs 69 that cause imperfections in flower and fruit development 70 . Temperature-based climate variables influence fruit growth, quality and seasonal yield 71 while temperature fluctuations directly affect citrus physiology 72 and indirectly reduce yield 73 . Temperature and relative humidity (RH) are used to compute vapor pressure deficit (VPD), is a difference in air saturation capacity (es) to actual water content/vapor (e) in air 74,75 , which is used to identify eco-physiology and hydraulic traits of plant growth 76 , since global warming has shown a fluctuating year-round pattern and a rising level 77 . VPD regulates stomatal conductance 78 which triggers its cavity to open at lower and closer at higher levels 79 . Leaf to air VPD is widening in high temperature and low RH regime and vice versa, to affect photosynthesis and transpiration 74,80,81 , is also an influential tool to judge ecological behavior of a plant under varying environmental conditions 76 while its fluctuation or elevation directly affects plant growth 82 . High VPD was recorded in the warm conditions of Vehari and TTS, especially in the starting months of summer season, to cause stress to Kinnow plant in both seasons, resulting in more fruit drops, which further developed unevenness in both fruiting seasons. In elevated VPD regime, more water loss due to evapotranspiration, resulting in plant-water imbalance 83 causes plant physiological process malfunction 84 with low carbon assimilation and high transpiration rates resulting in carbon starvation and hydraulic break-down 85 . VDP not only indicates temperature and RH regime of a particular area, but also dictates leaf stomatal conductance, transpiration rate, carbon assimilation, uptake of nutrients, plant-water hydraulic control and stress conditions 86,87 . It also controls gases exchange and defines stomatal limitation by recording a higher range in summer and a lower in winter 80,[88][89][90] . However both RH and temperature different regimes decide the extent of diseases, pests and fruit quality of citrus 15 by recording more pest damage in dry conditions and diseases infestation in warm-humid areas to affect fruit quality 11 which has justified the present work of showing more unevenness in fruiting patterns in warm conditions due to more fluctuation in climate variables. Climate variables also influence on citrus external quality features based on color development 6 as well as internal changes in the levels of juice contents, acidity, sugars and polyphenols 12 .
Thermal indices and energy use efficiency. Growing  Thermal use efficiency or thermal energy use efficiency showed a higher trend in heavy-fruiting plants during on-year and lower in light-fruiting in off-year due to significant differences in Kinnow fruit yield. However, Table 1. Climate variables of experimental sites. RTD (relative temperature disparity), T p (thermo temperature), T n (nytco temperature), RHD (relative humidity disparity), RH p (photo relative humidity), RH n (nytco relative humidity), VPD (vapor pressure deficit), VPD p (photo vapor pressure deficit), VPD n (nytco vapor pressure deficit).  Tables 2 and 3. Temperature counted directly GDDs, different heat units and indirectly PTI, HYTUs, PTU and HTU of a specific region on the basis of the threshold temperature of certain crops 11,60,91 . Citrus growth ceases below 13 °C 92 by keeping this temperature as threshold when computing growing degree days (GDDs) and other heat units of subtracting base temperature from mean daily temperature 10,93 . Vehari region has higher mean daily temperature, followed by TTS and lower was recorded in Sargodha 15 , accordingly GDDs, CHU, mCHU, p Heat unit and n Heat unit were counted in all three districts with leading counts in Vehari and lesser in Sargodha 6,12 . Table 2. Thermal indices of experimental sites. GDDs (growing degree days), CHU (citrus heat unit), mCHU (modified citrus heat unit), HYTUs (hydrothermal units), PTI (photothermal index), PTU (photo thermal unit), HTU (helio thermal unit), p (photo period), n (nytco period).   9 . More GDDs, CHU, mCHU, pHu and nHu were counted in the warm region of Vehari, thus influencing fruit growth and development from fruit-setting till arrival of maturity and subsequently the fruiting habit of the Kinnow mandarin. Hydrothermal units (HYTUs) were computed directly from average relative humidity (RHa) by multiplying with GDDs 9,11 and recorded more in warm-humid regions like TTS 6,12 . In this work, the GDDs were more computed at Vehari around the year, but with low level of RHa and vice versa in the case of Sargodha region, henceforth the counts of HYTUs in these two districts were lower than those of TTS. However, a fluctuating trend in HYTUs counts has been seen in three districts from fruit-setting to maturity due to climate variables supported by the work of 9 on Kinnow mandarin in India to record different HYTUs counts during fruit development. There was an increasing trend in HYTUs count during the summer season due to high daily means temperature and RHa during the rainy season compared to the rest of the months 12 . A fluctuating degree of HYTUs has been seen in all three districts due to the fact that they are located in different agro-ecological and agro-climatic areas, affecting the growth and development of Kinnow fruit and ultimately affecting fruiting habits as well as harvesting patterns. Similar results have been found 5,9 in India.
Photothermal index (PTI) is directly computed from mean daily temperature divided by time span in particular regions 5,96 and varies in different growth phases of citrus fruit 9 by showing a fluctuating trend across the year and locations 6 . In present work, PTI was more counted in warm regions and also in summer months due to higher mean daily temperature as seen in the Vehari district. PTI indicates daily photothermal index of a particular area 9 and can show a fluctuating trend in fruit different phenophases by recording more in cell division and cell enlargement stages 10 due to elevated temperature when Kinnow fruit these stages are continue, especially in summer months. PTI directly decides the span of citrus fruit different growth phases and indirectly fruit-setting and harvesting trend 10 and eventually bearing habit.
Photo and Helio thermal units reflect day-length, bright sun-light period and temperature levels of specific area and are directly computed by multiplying GDDs with day length and bright sunshine hours 5,9 . PTU and HTU were recorded more in warm regions and also in summer period due to elevated temperature and additional day-length than the cool season in the winter months 6 . Similar findings were seen in present work with additional PTU and HTU counts in Vehari district in both seasons and lesser in Sargodha from fruit-set till maturity. Same results of fluctuating PTU and HTU levels in different plum genotypes were reported by 5 in different climatic conditions of semi-arid regions in India to endorse this work. Climate variables in the three districts have distinct differences in mean daily temperature and length of the day; as a result, major changes have been seen in their monthly and annual counts. Climate variables decides meteorological indicators or thermal indices of a certain crop grown in a particular regions 2,4 to directly influence citrus fruit-setting, growth and development phases 10 , quality and yield features 11 , color-break and color development 6 as well as biochemical properties and maturity index 12 . Similarly, climate variables also decides fruit-setting time 9 , maturity arrival 3 and harvesting and marketing of Kinnow fruit 15 which significantly contributed to the fruiting pattern to justify this work. Thermal/ energy use efficiencies are computed directly on yield basis by dividing thermal indices to depict higher levels in high yielding crop and low thermal indices regions 2,96 by indicating positive yield correlation and inversely with different heat units accretion based on temperature, RHa, bright sun-light and day-length 4,9 . Temperature and RH based thermal use efficiencies were recorded higher in TTS during on-year on heavy-fruiting plants due to added yield but more thermal indices than Sargodha where yields were low. Although, thermal indices were counted more at Vehari except HYTUs but yield on all three fruiting habit plants was less during both on-& off-years than TTS and Sargodha. Similar findings regarding low yield and higher agrometeorological indices to computed lesser energy use efficiency was reported by 9 and 4 to support present work of calculating lower thermal use efficiency on low yielding plants and districts. More yields on heavy-fruiting plants during on-year were recorded in all three sites, henceforth, energy use efficiencies were computed to be higher than light-fruiting plants in off-year period to induce irregular bearings and further fluctuations in fruiting habits were caused by extreme climate variables which adversely affected plant physiology and tree potency in uneven fruit bearings.  Table 4, indicating significant differences in both the three experimental sites and the on-& off-years fruiting patterns. Fruiting patterns of citrus plants are directly influenced by environmental conditions 13 , cultivars with a bearing habit either on a single plant or on any branch or in a whole cluster / block form 28 thus exhibiting a superfluous predisposition in the mandarin group 37,39 . Heavy fruiting plant depletes carbohydrates and leaves low photoassimilates for next season flowering or fruiting 44 which imposes competition for carbohydrates and eventually resulted in alternate bearing habits in citrus 45 , mango 48 and avocado 46 . Harvesting time decides next season crop 97 while timely harvesting of heavy fruited plants induces evenness in fruit-bearing habit 29 . In this work, Kinnow plants were harvested late in warm conditions due to inland market consumption and late arrival of maturity indicators for the choice of the native consumers, which resulted in more disproportion of fruiting patterns in the Vehari and TTS districts. On the other hand, early harvesting and spot-picking began in the Sargodha area, which had reduced on-tree fruit load by controlling symmetry patterns in fruit-setting and fruitbearing habits for both on-and off-year fruiting seasons. In addition, extreme climatic variables also prompt asymmetry in the fruiting pattern 28,29 and its magnitude was more recorded in warm regions 15 to substantiate the present work implying unevenness in Kinnow plant fruiting in Vehari and TTS.

Vegetative flush trend in different fruiting habit plants. Kinnow plants have three vegetative flushes,
like other citrus cultivars, which showed a highly significant difference amongst spring, summer and autumn flushes and a slight difference in orchards grown in three different agro-climatic zones, as shown in Table 5. In the spring, vegetative flush was quantified as maximum on heavy-fruiting plants at Sargodha (60%) during on-year and minimum on light-fruiting plants at TTS (48%) in off-year. Maximum summer flush counts were recorded at TTS (35%) on light-fruiting plant during off-year and minimum at Sargodha (28%) on heavyfruiting plant in on-year. However, the higher autumn flush was counted at TTS (18%) on light-fruiting plants in off-year and the lower on heavy-fruiting plants in Sargodha (12%) during on-year.
Kinnow plants like rest of citrus cultivars have more than half percentage of spring, one-third summer and one-sixth autumn flushes 98 . In this work, spring vegetative flush was recorded higher in on-& off-years in all three fruiting habits plants and less was recorded in the autumn season. Similar findings of less summer flush than spring were attributed to harsh environmental conditions 99 . In all three districts, summer vegetative flush was also low than spring due to harsh external conditions, as well as on-tree fruits competed for carbohydrates during the summer months, when fruit cell division and cell enlargement stages were ongoing, with fruit expanding to maximum size 10 . Prior completion of cell division with the earlier commencement of the cell enlargement stages was observed in warm conditions 100 to justify this work of relatively more quantification of summer flush in warm districts Vehari and TTS. In addition, GDDs, PTI, PTU and HTU were accrued as higher in Vehari and TTS during cell division stage, thus this stage earlier wrap up its process than Sargodha and then begin prior cell enlargement 10 . Reported 43 that large chunks of photoassimilates were consumed during fruit growth and development in heavy-fruiting plants, while minutes share was streamlined towards vegetative growth to endorse this study of quantifying less spring and summer vegetative flushes on heavy-fruiting plants. The source-sink relationship for carbohydrates also defines the reproductive and vegetative growth habits of citrus plants 40,41 , while the fruit-load also restricts vegetative and root growth in mandarin to impute alternative bearings 42 .
In the off-year era, fewer carbohydrates consumed by fruit for growth and development, hence, more vegetative flush was recorded in light-fruiting plant during the off-season period and vice versa in heavy-fruiting plants during the on-year. In the off-season, fewer carbohydrates were consumed; thus, the next season (on-year) spring vegetative flush was recorded as more in all three experimental sites in all fruiting habit plants. Similar findings  101 , which also justify this study. In autumn season, meteorological indices like GDDs, HYTUs, PTI, PTU and HTU were less available to Kinnow plant in all three sites due to low means daily temperature with day-length squeezing by slowing down the net photosynthesis rate and retarding vegetative growth, wherein less photoassimilates synthesized with low carbohydrates accumulation in plant parts, although on-tree hanging fruits were in maturing and ripening phase to compete less for carbohydrates. These findings are Harvesting trend of Kinnow orchards. Harvesting trend data for Kinnow orchards are shown in Table 6, which shows major variations in harvesting times as well as trends among three districts. Harvesting of Kinnow orchards was recorded as maximum during off-year period (46%) in Sargodha and minimum during on-year (20%) in Vehari at the end of December. Harvesting trend increased during off-year season at Sargodha (56 and 70%) and recorded low during on-year in Vehari (35 and 48%) in the midst of January and the end of January, respectively. Same increasing trend of orchard harvesting in on-year season was seen in Sargodha (85 and 92%) in mid and end of February, respectively. However, rapid harvesting began in February in both the TTS and Vehari districts, reaching over 80% during off-year period and above 70% in on-year season. In the mid-March, maximum orchards were harvested in Sargodha during the off-year (98%) compared to other districts in same period, slightly above 90%. In certain fruit crops, including citrus, a cyclical carry-over effects of previous year's fruit present on trees are dominantly competed for carbohydrate reserves 45 and/ or trigger phyto-inhibiting effects on next season's floral bud-break/ initiation 50 can stimulate irregular fruiting patterns, commonly called as biennial bearing 38 while on-tree fruits influence on plant metabolism, like changes in coenzymes, sugars and amino acids, which are being accelerated when fruits are harvested late 49 . Similar findings of late harvesting in avocado induce biennial bearing with reduced fruit yields 106 . In this research endeavor, more alternative fruiting was recorded in TTS and Vehari due to delays in harvesting while heavy-fruiting during on-year season exhausted sugar and carbohydrate reserves in Kinnow plant to simulate low induction of floral buds and hence light fruiting on trees for the coming year (off-season). The fluctuating trend in carbohydrates reserves in citrus cultivars contributes to uneven in fruiting 107 , is not an exclusive phenomenon of citrus 101 but also observed in avocado 108 and olive 47 . Fruits load depletes carbohydrate in all plant components 45 with the leading role of roots carbohydrates of supplying energy for next floral and vegetative buds initiation in citrus 109,110 . In this work, both heavy-fruiting and late harvesting plants were deprived off carbohydrates and energy reserves, especially for Kinnow orchards of TTS and Vehari, in order to impute a more biennial fruiting pattern. In Sargodha, both timely harvesting and reduction of fruit loads by spot picking used in the export of Kinnow fruit have resulted in even fruiting for the next season by maintaining a balance of carbohydrates for both seasonal fruits. Similarly, in Satsuma mandarin (Citrus unshiu Marc.), Nishikawa et al. 49 found more sugar phosphate in light bearing vegetative stems and an additional accumulation of amino acids in heavy-fruiting trees to infer that fruiting habits had a direct impact on coenzymes, sugars and amino acids and had an indirect propensity on flowering and fruit-setting to justify this work. Heavy-fruiting and late harvesting plants have shown more alternative fruiting habits, as observed in this work, equally justified nutritional theory on the accessibility of carbohydrate for next season fruiting 45,[51][52][53]104 as well as on-tree fruit load inhibitory effect on coming season flowering 50,52,94,102 . The present work is in line with the findings of Monselise and Goldschmidt 28 Table 7. In both fruiting years and fruiting habits, the physiological fruit drops were seen higher in all three districts and the maximum was recorded in Vehari (60%) on heavy-fruiting plants during on-year season and minimum in Sargodha (47%) on light-fruiting plants during off-year period. However, pathological fruit drops were more recoded on heavy-fruiting plants at Sargodha (37%) during on-year and less on light-fruiting plants in Vehari (25%) during off-year. Total fruit Table 6. Harvesting trend of Kinnow orchards at experimental sites during on-& off-years. Results are shown in means (± SD) and sharing different letters are significantly differed to each other according to LSD test (P ≤ 0.05).  111 and sweet cherry 112 while among other favorable conditions; the availability of carbohydrates also increases flowering as well as fruit set with reduced drops in many fruits, including citrus fruits 113 . Fruits thinning, either natural or artificial, can increase the supply of carbohydrates and further reduce the drop of developing fruitlets 114,115 . Spiegel-Roy and Goldschmidt, 1996 116 referred to as physiological fruit drop in flowering and developing fruitlets during initial fruit stage, while 43 concluded that carbohydrates, especially soluble sugars, were available during fruit cell division for fruitlets retention on tree. In Kinnow mandarin physiological fruit drop ranges 40 to 63% and pathological fruit drops from 5 to 25% 117 . Physiological fruit drop is dominant during cell division stage 118 which lasted for 70-75 days in Kinnow 10 and has been linked to environment, nutrition and plant-water-balance 119,120 as well as pest pressure 121 . Similarly, the pathological drop in citrus fruit continued throughout the fruit growth and development stages before harvesting 11,118 with a dominant effect of adverse weather to proliferate diseases and mature fruit pests such as fruit fly infestation 15 . Same trend of more physiological and pathological fruit drop were seen on heavy bearing plants during on-year and less on light bearing trees in off-years, is in line with previous researchers based on the availability of carbohydrates and developing fruitlets competition for photoassimilates. In addition, dry warm conditions during fruit initial growth stages caused more drops, termed as physiological drop, whereas extended spell of warm-humid conditions in the autumn season has proliferated diseases and resulted in increased pathological drop. Citrus plant physiology was affected by a fluctuating pattern in agrometeorological indices, with different extents of fruit drops seen under climatic variable conditions at three experimental sites, as well as altering the fruiting habit of Kinnow mandarin. Along with weather vagaries 61 , carbohydrate deficiency reduced the induction of floral buds in the coming season 122 to support the hypothesis that climate variables induce fruit drops as well as alternative fruiting patterns in citrus fruits.
Yield and fruit grade quality of Kinnow plants. Yield  Climate variables have a direct impact on fruit growth and development 10 and consequently decreased yields with a decline in fruit quality attributes 123 . Holland et al. 124 recorded a decrease in fruit quality in the era of global warming, while 73 estimated a 1/4th reduction in citrus yields in the US. In this work, Kinnow fruits of inferior quality were harvested in warm regions during on-and off-years from light, medium and heavy bearing plants. In addition, higher temperatures during the fruit cell division stage led to even more physiological drops 10,11 resulting in alternate bearings 13 and finally yielding low-quality fruit 15 . Chelong and Sdoodee 125 found a direct effect of climate variables on fruit yield and quality in their work on Shogan (Citrus reticulata Blanco) in Thailand. In present work, warm dry spells at TTS and Vehari have resulted in a more premature stage Kinnow fruit drops while extending the summer season period has exacerbated stem-end rot disease, causing additional drops in maturing fruits. In fact, both light and heavy fruiting seasons, the fruit drop pattern was related to the external climate, which decreases yield and grade quality while showing further variations in three experimental sites. www.nature.com/scientificreports/ More declines in yield and grading quality have been seen in light fruiting seasons with high intensity in warm regions. However, during heavy fruiting season, hanging fruits were least affected by the vagaries of weather conditions. Citrus fruit A-grade quality is linked to fruit size, shape and apparent view 15,126 or blemishes free outer peel 14 , which was more harvested in Sargodha due to timely harvesting and early spot-picking for export while observing less fluctuating climate variables. The physiological mechanism of biennial fruiting trends in citrus 129 is directly linked with external conditions, especially changing vapor pressure deficit 130 during fruit growth and development phases 10 to determine fruit-load for succeeding season. In addition, the relationship between the sources (leaves) photoassimilates and the sink (fruits) also has an effect on the yield and quality of the citrus fruit 127,128 which has justified this work of recording more A-grade quality fruit on medium bearing trees during on-year. More pest pressure, particularly citrus mites and thrips, was seen in warm areas that directly affected the external outlook of the fruit 11 , reducing the exportable volume 15 , which was recorded higher in warm TSS and Vehari districts to justify less A-grade quality Kinnow fruits produced in more climate variable regions. Present work is in line with previous works on citrus 13,129 .

Conclusion
Climate variables determine the fruiting habit, yield and quality attributes of Kinnow Mandarin. In warm regions, the fluctuating trend in thermal indices has not only influenced plant phenophases, but also fruiting habits, fruit drops at different stages and, consequently, yield and quality characteristics. In TTS and Vehari, the more unpredictable weather behavior resulting in more variations in thermal indices causes an alternating pattern of fruiting by disrupting the source-sink relationship and deteriorating fruit quality, henceforth affecting plant thermal use energy efficiencies. Biennial fruiting pattern is an inherent character linked to citrus is dominantly induced by external environment has become more prevalent in more climate variable regions as seen in this work. This research endeavor may be fruitful in future to decide particular region regarding citrus fruiting habit, quality as well as yield and also pinpoint major management practices in future where more fluctuations in climate variables arise.