Received 22 February 2016; accepted 22 April 2016; published 25 April 2016
Tea (Camellia sinensis (L). O. Kuntze), a perennial tree belonging to the Theacea family   , is managed as a low bush in continuous phase of vegetative growth. The crop provides income to farmers and employment to rural populations in many parts of the world. In Kenya, 377 million kilograms of made tea was produced from 171,916 Ha, in 2011, contributing to 4% of GDP and 28% of export earnings  . This was the highest export earnings from a single commodity and/or crop. The crop is grown in the high rainfall areas of Kenya, namely, foothills of the Aberdare Ranges and Mt Kenya, in the East of the Great Rift Valley and the Mau Ranges; Nandi, Kisii and Kakamega hills and slopes of Mt. Elgon in the West of the Rift Valley  . These areas range in altitude from 1300 m to 2700 m amsl   . These areas lie almost along the equator, and shoots are harvested at 7 to 14 days intervals throughout the year  . However, yields   and quality  -  are affected by fluctuations in weather factors within and between the years in any one location. Despite production close to the equator, the differences in geographical areas of production have been recorded to influence leaf nutrients levels  ; growth rates  -  , productivity  -  , precursors of tea quality parameters  -  and black tea quality   -  of tea. Despite use of same genotypes and management practices, yields and quality vary across the region. This suggests growth patterns in different genotypes may not be the same in the region. There is need to assess possible variations in growth patterns of different genotypes to determine whether yields can be predicted from the yields components in different locations.
Tea productivity    and quality   varied with altitude and weather factors   . There was yield reduction during dry and cool seasons   . Indeed, the variations in yield performance of genotypes relative to each other and to the environment have been widely documented    with the phenomenon having been formulated in the Genotype by Environment (G × E) interaction concept  . However, the mechanisms causing the variations in yields are not well defined, especially since the level and extents of the variations change with genotypes. The understanding of the variations in responses in tea yield components with genotypes and locations may provide an explanation in the observed yield responses and lead to development of management practices that enhance productivity in different clones in different locations.
2.1. Experimental Treatments and Design
The trial was set up in 2011 in three different tea growing geographic regions (Kangaita, Timbilil and Kipkebe) varying in proximity by between 45 and 245 Km and altitudinal range differences of between 78 an 435 m, namely (Table 1). Slopes at all the sites are gentle to slightly sloping (0% - 15%). The trials were set up on plots of mature tea of the same age, planted in 1991 and comprising twenty clones namely; TRFK 6/8, TRFK 31/8, AHP S15/10, EPK TN14-3, BBK 35, TRFK 54/40, TRFK12/12, TRFK 12/19, TRFK 31/27, TRFK 11/26, TRFK 57/15, TRFK 7/3, TRFK 7/9, TRFK 56/89, STCK 5/3, TRFK 303/259, TRFK 303/577, TRFK 303/999, TRFK 303/1199 and TRFK 2X1/4 set in Randomized Complete Block design replicated three times at each site  . Each plot comprised of 20 bushes. The tea was managed under standard management practices in Kenya  . The yield components recorded were shoot density, dry weight, growth rate and water potential. The experiment was analysed as a 2 × 3 factorial 2 design with the clones and sites as the treatments.
2.2. Data Sampling and Recording
2.2.1. Soil Characteristics
Disturbed soil samples were collected from two sites from each trial location in approximate diagonal line across the experiment. The samples were collected in the dry season between January and March at depths of 0 - 20, 20 - 40, 40 - 60   using a Jarret auger. These were subjected to full chemical (pH and nutrients) analysis and physical (soil texture) analysis for site characterization.
Table 1. Study sites geographic location.
2.2.2. Soil Chemical and Physical Analysis
Fresh (un-dried) sub samples were subjected to analysis of soil pH by reading off a Jenway 3305 pH meter. Nitrogen content was determined using the Kjeldahl method. For the mineral nutrient analysis, the soils were air dried, ground and sieved through a 2 mm. sieve. The ground samples were then extracted using the Mehlich III method then analysed for nutrients K, P, Ca, Mg, Mn, Na, Cu, Fe and Zn using a plasma atomic emission spectrophotometer (ICPE-9000, Shimandzu). Sub samples from the disturbed soil samples were subjected to particle size analysis using the pipette method  , taking 63 μm as the sand/silt boundary.
2.2.3. Site Weather Characteristics
Rainfall and temperature data were recorded from weather stations located at each trial site and accompanying meteorological data derived as follows; Rainfall was recorded daily using a standard rain gauge. Maximum, minimum, wet and dry bulb temperatures were recorded everyday at 09.00 h and 15.00 h local time using mercury in glass thermometers (Cassella (London) Ltd., UK). The relative humidity (RH%) was derived from the wet and dry bulb temperature readings using relevant tables  . The wet and dry thermometer readings recorded as described above were used to derive saturated vapour pressure deficit (SVPD) using the formula. SVPD = ew − e' 
e' = air vapour pressure (mb) at t' (Table 94),
ew = air vapour pressure (mb) at t (Table 94),
t' = dry bulb temperature (˚C),
t = wet bulb temperature (˚C).
Green leaf comprising of mostly two leaves and a bud were plucked every 7 - 10 days and converted to made tea (mt) by multiplying by a factor of 0.225  .
2.2.5. Shoot Density (SD), and Shoot Dry Weight (SDWT)
Shoot density was determined by getting the mean of the number of mature harvestable shoots (two leaves and a bud) captured within a 0.25 m2 grid randomly thrown on to the plucking tables  of five randomly selected bushes at every plucking round. Shoots falling within the grid were plucked, weighed and counted. The shoots were then oven dried at 105˚C for 48 hours and weighed. The SDWT was determined by dividing the dry weight by the number of shoots harvested.
2.2.6. Shoot Growth Rate (SGR)
The rate of shoot growth (millimetres per day (mmd−1)) was determined by tagging five shoots from each of three randomly selected bushes per plot. Growth was monitored by measuring the length from the tip of the auxiliary bud to the base every three days until the new shoot developed into a mature harvestable two leaves and a bud. The total length measured at each interval was divided by the number of days between two successive measurements to determine the growth rate (mmd−1).
2.2.7. Shoot/Xylem Water Potential (SWP)
The xylem water potential of pluckable shoots was measured between 11.00 am and 2.00 pm, using the pressure chamber technique  . Three shoots each from five randomly selected bushes per plot were cut and measured one by one on each plucking day. The shoots were and transferred to the site of the pressure chamber where one centimetre was cut off the stalk and the shoot immediately inserted into the gas chamber. The key of the compressed nitrogen gas turned on until the first gas bubbles were released from the cut shoot stalk surface and the pressure reading taken. These measurements were taken twice a season and averaged to get the season mean.
2.3. Data Analysis
The data collected was subjected to analysis of variance (ANOVA) using MSTAT-C (Version 2.10) statistical package, as a factorial two design, with clone (genotype) as the main factor and location as the second factor. Correlations between yield and yield components and weather parameters were done using SPSS (Version 17.0) statistical software.
3. Results and Discussion
3.1. Soil Variability and Suitability
The soil characteristics of the experimental sites are given in Table 2 and Table 3. The textural properties of the soils showed variations, with Kangaita soils having coarser texture than Timbilil and Kipkebe soils. Kangaita had higher sand and lower clay contents than the other two sites, which were similar. The porosity ranged from 38% to 56%, with little variation between the sites. The soils from all the sites were of volcanic origin  . Tea is grown in soils of varying texture and in soils with clay content as high as 83% in Kenya and as low as 1.7% in Taiwan  . The soils in this study fell within these ranges and were similar to those observed in Kericho  ,
Table 2. Soil physical characteristics of the trial sites.
*Soil description after Jaetzold and Schmidt, 2010.
Table 3. Soil chemical characteristics of trial sites.
and sandy loam in Kangaita in Kangaita and Nyambene  . These results demonstrate the ability of tea to grow in a wide variety of soil textural types. The pH of the soils from the three sites ranged from 4.0 to 3.4 but were predominantly lower than 4.0 (Table 3). There was very little variation in the soil pH between the sites. Tea grows in soils of optimal pH of 4.0 to 6.0  , but can grow in pH below 4.0  . Indeed, optimal growth of tea had been reported at pH ranges between 3.8 and 5.7 in land newly cleared from primary, secondary forests and tree plantations   . The soil mineral contents (Table 3) were within the ranges observed in the major tea growing areas  . The soil nitrogen contents were adequate for tea growth. However, site variations were evident with Kangaita recording higher nitrogen levels than the other locations. These results demonstrate that the soils were suitable for tea growing. Thus, constraints or low tea productivity at any of the sites would be due to other factors. These results demonstrate that though tea can be successfully and commercially grown in a wide variety of soil types, the variations in soil types contribute significantly to locational yield variations.   observed soil differences between sites and even significantly big differences between sites across East Africa tea growing regions. He observed that soil water deficits, compaction and soil water deficits could restrict productivity of tea in different locations. Indeed results of soil water potential which is directly related to plant water status shows significant site variations in this study.
3.2. Weather and Geographical Locations
The weather components from the three trial sites during the duration of the trial are shown in Table 6. There were noticeable differences between the sites in all the weather parameters measured and derived (rainfall, rain days, ambient temperatures, relative humidity and vapour pressure deficit).
Mean monthly temperatures generally rose with declining altitude from Timbilil to Kipkebe (Table 1 and Table 4). However, Kipkebe had the highest mean ambient temperatures but mean Kangaita and Timbilil temperatures were similar despite the difference in altitude (Table 1). This was due the close proximity of Kangaita site to the Mt. Kenya and the Mt Kenya Forest which lowered mean daily temperatures. The mean temperature difference between Kangaita and Timbilil was only 1˚C while between Kipkebe and Timbilil and Kipkebe and Kangaita were 3˚C and 4˚C, respectively. These results were similar to those obtained earlier  where temperature difference between locations at high altitude (over 2000 m) was lower than those between high and lower altitude (below 2000 m) locations. Therefore the two locations, Kangaita and Timbilil, both lying above 2000 m were similar in terms of mean annual temperatures despite their altitudinal difference, but significantly different from the low altitude Kipkebe site.
The rainfall pattern in Kangaita was bimodal with peaks in April-May and October-November, separated by cold and hot dry seasons from June to August and January to March, respectively. Rainfall pattern in Timbilil was unimodal with rains starting in April and continuing up to December. In Kipkebe the rainfall was weakly bimodal with peaks in April to July and in December (Table 4) as had been reported earlier for Kericho  . Total rainfall declined with decreasing altitude from Timbilil to Kipkebe (Table 1 and Table 4). However, these relationships were not significant. Timbilil and Kipkebe had more rain days during the year but were drier in the January-March season, with no rain in the month of January. These seasonal variations are likely to impact on total annual yields variations between the sites.
3.5. Relative Humidity (Rh) and Vapour Pressure Deficit (VPD)
Relative humidity (Rh) and vapour pressure deficit (VPD) also showed variations between sites. Mean monthly relative humidity increased from high to low altitude. The high humidity recorded at Kipkebe was attributed to its closer proximity to the Lake Victoria. Monthly VPD was highest in Timbilil and lowest in Kangaita. Highest VPDs were recorded during the January March period at all sites but was highest at Timbilil. This suggests more severe drought in Timbilil than at Kipkebe and Kangaita. Seasonal variations of the weather parameters within and between sites were also evident. Soil water deficits are determined by ambient temperatures and humidity, similar to shoot water potential (SWP). Atmospheric humidity is inversely related to vapour pressure
Table 4. Monthly weather parameters at all the trial sites, Jan-Dec 2012.
deficits   . An inverse linear relationship between VPD and SWP in tea has been reported    , although the SWP of tea shoots were more closely related to VPD than to soil moisture   reported variations in soil water deficits between sites in Kericho, which affected the difference in yield variations.
3.6. Yield Components
3.7. Shoot Growth Rate
The shoot growth rates (Table 5) in Kangaita and Timbilil were similar, but lower (p ≤ 0.001) than that of Kipkebe. Plant growth rates at the two sites followed the locational temperature patterns (Table 4) and were thus expected to be similar or insignificantly different since temperature is a major determinant of growth as reported in previous studies  -  . The clones x sites interactions were also significant (p ≤ 0.001) demonstrating the shoot growth rates of the various cultivars were not responding in similar patterns in all sites. There was a positive response of shoot growth rates to temperature but an inverse response to vapour pressure deficit across three locations (Table 5 and Table 6). Clonal variations in shoot growth rates varied both within and across locations. Indeed, variations in clonal growth rates were not consistent across all sites as the shoot growth rate rankings showed variation with location (Table 6). Some clones however, were ranked in the same quartile across two or three locations. The differences between the highest and lowest clonal growth rates also varied with location with the highest (0.56 mm/day) in Kipkebe and the lowest in Kangaita (0.15 mm/day). Similar results were reported in earlier findings of clonal variation within one location   across seasons  and across locations    . Shoot growth rate will vary with genotype and the genotypic variation will also vary from one location to the other.  opined that changes in shoot growth rates could be explained by the inherent effects of temperature and relative humidity, which varied across seasons. Temperature and relative humidity indeed, determine the prevailing vapour pressure deficits in any one location. The findings by  therefore explain the relations between shoot growth rates and saturated vapour pressure deficits observed in this study. Vapour pressure deficits are therefore a major weather parameter other than temperature that determines the shoot growth rate of the tea plant.
3.8. Shoot Dry Weight and Shoot Density
Shoot densities showed significant (p ≤ 0.001) variations due to genotype and location. The genotype x location interactions was also significant (p ≤ 0.001). Whereas shoot dry weights varied (p ≤ 0.001) between sites there was no apparent response to temperature or vapour pressure deficit. The shoot densities variation (p ≤ 0.001) due to location however, correlated positively with vapour pressure deficits but not temperature, across the sites. The highest mean shoot density was recorded in Kangaita (119 shoots m−2) and lowest at Timbilil (32 shoots m−2).  reported varying clonal response of shoot density with altitude.
Table 5. Effect of geographical location and genotype on yield components.
Cl = clone; St. = site.
Table 6. Effect of geographical location on ranking of genotype yield and components.
3.9. Shoot Water Potential
Shoot water potential showed no variation between genotypes but showed highly significant (p ≤ 0.001) site variations. The changes followed the pattern of vapour pressure deficits across the sites (Table 4). Indeed,  demonstrated that shoot water potential of tea shoots were more closely related to vapour pressure deficits than to soil moisture. Shoot water potential influences shoot growth by determining the cellular turgidity required for cell expansion  . These findings suggest a threshold of plant water status for tea shoot survival and growth, and this does not vary with genotypes but varies with location. The threshold for plant water status is therefore determined by location weather attributes that dictate the vapour pressure deficit. Selection for optimal yield performance across locations should therefore take cognisance of prevailing vapour pressure deficits as the shoot growth rates determine the potential yields attainable  . Management practices to optimise yields may also yield better results if measures to alleviate vapour pressure deficits are effected.
3.10. Annual Yields
There were significant (p ≤ 0.001) genotypic and site differences in annual yields (Table 7). The clones x site interactions were also significant (p ≤ 0.001). There were significant yield variations due to genotype. The genotypic variations were manifested in all locations and the extents of the variations changed from location to location. Clonal mean yields were lowest (p ≤ 0.05) in Timbilil and highest in Kipkebe. This yield trend mimicked the locational temperature patterns (Table 4).
Table 7. Effect of geographical location and genotype on annual yields (kg∙ha−1) and their rankings.
*Consistent superior performance (1st quartile) across all three locations; #Consistent poor performance (4th quartile) across two locations; +Poor performance (4th quartile) across all three locations.
The clones and sites interaction were illustrated in the variations in clonal yield rankings between sites (Table 9); clone TRFK 303/577 recorded highest yield (5162 kg mt/ha) at Kipkebe while TRFK 7/3 recorded the lowest yield (763 kg mt/ha) at Timbilil. All clones recorded highest yields at Kipkebe. There were also variations in yield rankings among the clones between sites. Yield differences between the highest and lowest yielding clones also varied between sites. These differences were highest in Timbilil at 2410 kg and lowest in Kipkebe at 1473 kg. This phenomenon elucidated the variation in genotype yield response to the different environments. Clonal yield rankings showed clonal yield stability across only two environments in clones TRFK 303/577, TRFK 12/19, AHP S15/10 and TRFK 54/40.
Some clones showed consistent good yield performance across all environments, falling in the 1st quartile in all locations. However, consistent inferior yield performance (4th quartile) across all three locations was exhibited by only two clones (TRFK 7/3 and 6/8). Similarly, consistent poor performance under low temperature conditions (Kangaita and Timbilil) was exhibited by only one clone (STC 5/3). One clone, TRFK 11/26 showed consistent poor performance across two temperature regimes of Timbilil and Kipkebe. Majority of the clones did not show any consistent yield correlations with locations. These results demonstrate the variability in clonal yield response to environment as earlier reported, on yields in Kenya    , Tanzania  , across East Africa region  in Sri Lanka  and on catechin content   . However, the results could not be extrapolated to other regions unless the differences in yield were related to specific environmental variables such as soil water deficit (SWD) and temperature  .
The differences between the three geographical locations can be explained in terms of the weather differences (Table 4), as had been observed in earlier studies  -  . The earlier reported studies however, involved sites around Kericho, within a radius of 20 Km   and two only sites, with different weather patterns (Kericho and Kangaita)  . The findings from this study replicate findings reported from single sites and also from multiple sites above. Clonal yield variations were observed to occur under different weather patterns (Kericho and Kangaita) and different temperature regimes (Kericho/Kangaita and Sotik). This means that despite close similarities in weather parameters between locations, tea genotypes will respond differently in the locations. Hence, different clones with similar performance under the same management in one site will require different management options to optimize yields in a separate location.  made similar observations. Therefore, importing clones selected for high yield in a different location may not optimise yields.
3.11. Shoot Components and Yield Interactions
Yield correlations showed that yield components’ correlation to yields varied with environment (Table 8). The magnitude of the correlation varied with locations. Just as  observed that shoot density and shoot weights may be important in determining yields between clones in one site, this does not apply to different clones in different geographical locations. The correlations showed that different yield components may be important in determining yields in varying locations. The magnitude of the correlations between yields and growth components varied with locations. Significant relationship between shoot growth rate and yield was observed in Timbilil (r = 0.476, p ≤ 0.05) but not at Kipkebe (r = 0.384) and Kangaita (r = 0.001) while shoot density showed significant relationship in Kangaita (r = 0.652, p ≤ 0.05) (Table 8). Observations by  that shoot density and shoot weights are important in determining yields between clones concurred with observation from Timbilil only. Shoot dry weights and shoot water potentials showed no significant correlation to yield. Thus the relationship is not universal. These results imply variation of yield components contribution to yield amongst genotypes. Yield selection using yield components may be genotype specific and therefore weighting of the components may need to be considered during genotype selection.
In Kangaita, all yield components in the model were additive to yield (Equation (1)). In Kipkebe and Timbilil however, the high temperatures resulted in low plant shoot water potentials thereby reducing potential yields (Equations (2) and (3)). The negative shoot water potential was attributed to the prolonged droughts experienced in the two locations during the period of the study.
S = Shoot growth rate,
D = Shoot density,
W = Shoot dry weight,
P = Shoot water potential.
Table 8. Effect of geographical location on relationship (r) between yields and yield components.
N = 20, **Correlation is significant at the 0.01 level; *Correlation is significant at the 0.05 level; Sgr = Shoot growth rate; Sd = Shoot density; Sdwt = Shoot dry weight; Swp = Shoot water potential.
Stepwise regression further revealed the yield components that significantly contributed to yields in each location. In Kangaita, shoot density was significant (r = 0.652; 0.425) as was the combination of shoot density and shoot dry weight (r = 0.792; 0.628). In Timbilil only shoot growth rate was significant in determining yield (r = 0.476; 0.226). For Kipkebe, the regression could not be generated, suggesting that at that location, the individual and combination of yield components do not met the significant threshold to determine yield. This therefore suggests that at Kipkebe all yield components may combine equally to determine potential yields.  reported similar results from 4 clones in Timbilil in Kericho and further reported that effects of individual components were highly varied and did not relate with potential clonal tea yields.
 showed a relationship between shoot sizes and shoot mass. The size of shoots harvested is an important factor in determining yield   . The regression models above indicate the importance of shoot weight in yield determination.  however, reported that shoot size contributed only 11% of the total seasonal yield variation, the remaining 89% being accounted for by the number of shoots. Seasonal yield variations mainly occur due to differences in shoot density  . However, later investigations   hypothesized that shoot growth rate was the major component causing seasonal fluctuation in yield while shoot density was the major factor determining yield difference between varieties. In this study the genotype and environment interactions show that the yield components contribution to yield in any environment varies significantly with the genotype. The findings reported above are results of studies conducted in single site with the varying environmental factor being the climatic seasons.
Yield components’ contribution to yield varied with genotype within the same environment and between environments (Table 8). This could explain the variation in yield stabilities across different environments. Similar results have been reported   illustrating the variations in clonal growth response to temperatures across four sites within Kericho region. Findings from this study indicate that variations in tea response to weather parameters may apply even across locations with different climatic regimes. The results therefore suggest that not all yield components can be used for yield selection in all environments. Indeed, in Kangaita clonal yield selection would best be done using shoot density as a yield indicator while in Timbilil shoot growth rate would be the best selection indicator.
Variation in clonal yield components’ response to environments explains the variability clonal yield response among environments. The interaction and the dominant component response determine the actual clonal yields realised in every location. Yield optimisation should therefore be targeted at the management practices that reduces the yield limiting factors and maximises the dominant responsive component e.g. sprinkler irrigation during periods of moisture stress to reduce soil water deficits and raise the relative humidity hence lowering the SVPD which is the main factor limiting growth in such conditions. Clonal selections will also have to be location specific. Clonal selection for universal planting will not always optimise yields across all locations. Yield modelling even for a single clone across all sites will at best be estimate not accurate. Perhaps this would explain the difficulty in modelling of yield predictions encountered in past works.
3.12. Yield and Yield Components’ Response to Weather
All yield components studied (shoot growth rates, shoot dry weights and shoot densities), showed significant clonal variations in their sensitivity to the different environments. Generally, shoot growth rate increased with temperature across the locations while shoot dry weights varied significantly (p ≤ 0.05) between sites but showed no response to temperature or vapour pressure deficit. All yield components gave significant (p ≤ 0.001) genotypic variations in their responses to the environment (Table 8). Yield components and weather trends presented in Table 4 and Table 5, respectively indicate general linear responses of yield components to weather parameters from low to high altitude (Kipkebe to Kericho). Shoot growth rates increase with increasing ambient temperatures across the sites, while shoot dry weights and shoot densities showed no apparent response to temperatures across the sites. Correlation analysis showed the same relationship but shoot density tended to reduce with rising ambient temperature. Across the three trial sites the temperatures tended to drop with increased rainfall leading to yield depression at Kangaita and Timbilil where ambient temperatures were lower. At Kipkebe however, lower rainfall did not depress yields as the temperatures remained warmer.
Correlation analysis however, show that the yield and individual yield component response to weather parameters vary significantly (p ≤ 0.05) between locations (Table 9). These correlations explain the locational variation in the yield models above (Table 9). The results imply that it is not possible to give an accurate yield model
Table 9. Effects of location on correlations of yield, yield components and weather parameters across three geographical regions of tea production.
*Significant at 0.05; **Significant at 0.01; Yld = Annual yield; Sgr = shoot growth rate; Sd = Shoot density; Sdwt = Shoot dry weight; Swp = Shoot water potential; Temp = Ambient temperature; Rain = Annual rainfall; Rh = Relative humidity; Rdays = Rain days; Svpd = Saturated vapour pressure deficit; Tir = Total incident radiation.
for yield prediction across all environments for all tea genotypes and that the general response of yield and yield component across environments is not universal.
Soil chemical and physical parameters varied between sites, but were suitable for tea growing. Yield components contribution to yield varied with location. The response of tea yield components to weather factors also varied with location. Genotype yields and yield components varied significantly in their response to environment, variations in correlations between yield and yield components indicate that the response of tea yield components to weather factors varied with location. These variations give rise to the yield variability between locations and between genotypes. Shoot density and shoot growth rates were more closely related to yields only in Kangaita and Timbilil, and might therefore be used as yield selection indicators in Kangaita and Timbilil. Yield components may not be used universally for clonal selection. Genotypes selected in any one location may perform optimally in another location. Selected clones should therefore be subjected to testing before adoption for commercial planting in another location.
The authors wish to thank Kenya Tea Development Authority (KTDA) and Kipkebe Tea Company Ltd. For providing trial sites, the Inter-University Council of East Africa: Lake Victoria Research (VicRes) Initiative and National Council of Science and Technology (NACOSTI) for funding of the study and the staff of Crop Environment and Physiology section, and Kangaita station of the Tea research Institute for their assistance in data collection. This paper is published with the approval of the Director of the Tea Research Institute of Kenya.
 Bokuchava, M.A. and Skobelava, N.I. (1969) The Chemistry and Biochemistry of Tea and Tea Manufacture. Advances in Food Research, 17, 215-292.
 Owuor, P.O., Obanda, M., Nyirenda, H.E. and Mandala, W.L. (2008) Influence of Region of Production on Clonal Black Tea Characteristics. Food Chemistry, 108, 263-271.
 Eden, T. (1931) Studies in the Yield of Tea I. The Experimental Errors of Field Experiments with Tea. Journal of Agricultural Science, 21, 547-573.
 Squire, G.R. (1979) Weather, Physiology and Seasonality of Tea (Camellia sinensis) Yield in Malawi. Experimental Agriculture, 15, 321-330.
 Tanton, T.W. (1992) Tea Crop Physiology. In: Willson, K.C. and Clifford, M., Eds., Tea: Cultivation to Consumption, Chapman and Hall, London, 173-199.
 Kwach, B.O., Kamau, D.M., Msomba; S.W., Muhoza, C. and Owuor, P.O. (2014) Effects of Location of Production, Nitrogenous Fertilizer Rates and Plucking Intervals on Clone TRFK 6/8 Tea in East Africa: II. Mature Leaf Nutrients. International Journal of Tea Science, 10, 25-40
 Ng’etich, W.K. and Stephens, W. (2001) Responses of Tea to Environment in Kenya. II. Genotype × Environment Interactions for Total Dry Matter Production and Yield. Experimental Agriculture, 37, 333-342.
 Squire, G.R., Obaga, S.M.O. and Othieno, C.O. (1993) Altitude, Temperatures and Shoot Production of Tea in the Kenyan Highlands. Experimental Agriculture, 29, 107-120.
 Wachira, F., Ng’etich, W., Omolo, J. and Mamati, G. (2002) Genotype & Environment Interactions for Tea Yields. Euphytica, 127, 289-296.
 Ng’etich, W.K., Stephens, W. and Othieno, C.O. (2001) Responses of Tea to Environment in Kenya. 3. Yield and Yield Distribution. Experimental Agriculture, 37, 361-372.
 Msomba, S.W., Kamau, D.M., Uwimana, M.A., Muhoza, C. and Owuor, P.O. (2014) Effects of Location of Production, Nitrogenous Fertilizer Rates and Plucking Intervals on Clone TRFK 6/8 Tea in East Africa: I. Yields. International Journal of Tea Science, 10, 14-24.
 Cherotich, L., Kamunya, S.M., Alakonya, A., Msomba, S.W., Uwimanna, M.A. and Owuor, P.O. (2013) Characterization of 20 Clones of Tea (Camellia Sinensis (L.) O. Kuntze) Using ISSR and SSR Markers. Agricultural Science Research Journal, 3, 292-302.
 Cherotich, L., Kamunya, S.M., Alakonya, A., Msomba, S.W., Uwimanna, M.A. and Owuor, P.O. (2014) Genotypic Stability and Adaptability of Tea Cultivars in Relation to Catechin Levels across Four Environments. Tea, 35, 8-16.
 Kwach, B.O., Owuor, P.O., Kamau, D.M., Msomba, S.W. and Uwimana, M.A. (2016) Variations in the Precursors of Plain Black Tea Quality Parameters Due to Location of Production and Nitrogen Fertilizer Rates in Eastern African Clonal Tea Leaves. Experimental Agriculture, 52, 266-278.
 Okal, A.W., Owuor, P.O., Kamau, D.M. and Mang’uro, L.O.A. (2012) Effects of Production Location and Plucking Interval on the Fatty Acids Concentration of Clonal Tea. Food Science and Technology Research, 18, 351-356.
 Okal, A.W., Owuor, P.O., Kamau, D.M. and Mang’uro, L.O.A. (2012) Variations of Fatty Acids Levels in Young Shoots of Clonal Tea with Location of Production and Nitrogenous Fertilizer Rates in the Kenya Highlands. Journal of Agricultural Science and Technology, 14, 1543-1554.
 Owuor, P.O., Kamau, D.M., Kamunya, S.M., Msomba, S.W., Uwimana, M.A., Okal, A.W. and Kwach, B.O. (2011) Effects of Genotype, Environment and Management on Yields and Quality of Black Tea. In: Lichtfouse, E., Ed., Genetics, Biofuels and Local Farming Systems: Sustainable Agriculture Reviews, Vol. 7, Springer, Heidelberg, 277-307.
 Owuor, P.O., Okal, A.W., Kamau, D.M., Msomba, S.W., Uwimana, M.A. and Kamunya, S.M. (2013) Influence of Nitrogen Fertilizer Rates and Harvesting Intervals on Clonal Tea Green Leaf Fatty Acids Levels in the Lake Victoria Basin of Kenya. Journal of Agriculture, Food and Environment, 11, 667-674.
 Owuor, P.O., Obaga, S.O. and Othieno, C.O. (1990) The Effects of Altitude on the Chemical Composition of Black Tea. Journal of Science Food and Agriculture, 50, 9-17.
 Owuor, P.O., Wachira, F.N. and Ng’etich, W.K. (2010) Influence of Region of Production on Relative Clonal Plain Tea Quality Parameters in Kenya. Food Chemistry, 119, 1168-1174.
 Owuor, P.O., Kamau, D.M., Kamunya, S.M., Msomba, S.W., Jondiko, E.O. and Uwimana, M.A. (2013) The Response of Clone BBK 35 Tea to Nitrogen Fertilizer Rates and Harvesting Intervals in the Lake Victoria Basin of Kenya. Journal of Agriculture, Food and Environment, 11, 757-763.
 Balasuriya, J. (1999) Shoot Population Density and Shoot Weight of Clonal Tea (Camellia Sinensis) at Different Altitudes in Sri Lanka. European Journal of Agronomy, 11, 123-130.
 Anandacoomaraswamy, A., De Costa, W.A.J.M., Shyamalie, H.W. and Campbell, G.S. (2000) Factors Controlling Transpiration of Mature Field-Grown Tea and Its Relationship with Yield. Agricultural and Forest Meteorology, 103, 375-386.
 Mahanta, P.K., Baruah, S., Owuor, P.O. and Murai, T. (1988) Flavour Volatiles of Assam Black Teas Manufactured from Different Plucking Standards and Orthodox Teas Manufactured from Different Altitudes of Darjeeling. Journal of the Science of Food and Agriculture, 45, 317-324.
 Tanton, T.W. (1979) Some Factors Limiting Yields of Tea (Camellia sinensis). Experimental Agriculture, 15, 187-191.
 Odhiambo, H.O., Nyabundi, J.O. and Chweya, J. (1993) Effects of Soil Moisture and Vapour Pressure Deficits on Shoot Growth and Yield of Tea in The Kenya Highlands. Experimental Agriculture, 29, 341-350.
 Wickramaratne, M.R.T. (1981) Genotype-Environmental Interaction in Tea (Camellia sinensis L.) and Their Implications in Tea Breeding Selections. The Journal of Agricultural Science, 96, 471-478.
 Gee, W. and Baader, J.W. (1986) Particle Size Analysis. In: Klute, A., Ed., Methods of Soil Analysis: I: Physical and Mineralogical Methods, Agronomy Services, No. 9. American Society of Agronomy, Inc., Soil Science Society of America, Inc., Madison, 384-411.
 List, J.R. (1949) Hygrometric and Physchrometric Tables. In: List, J.R., Ed., Smithsonian Meteorological Tables, Vol. 114, 6th Edition, 1971 Reprint, Section VII, Smithsonian Institution Press, Washington DC, 347-385.
 Odhiambo, H.O. (1989) Nitrogen Rates and Plucking Frequency on Tea: The Effect of Plucking Frequency and Nitrogenous Fertilizer Rates on Yield and Yield Components of Tea (Camellia sinensis (L.) O. Kuntze) in Kenya. Tea, 10, 90-96.
 Jaetzold, R., Schmidt, H., Hornetz, B. and Shisanya, C. (2010) Farm Management Handbook of Kenya, Vol. II, Natural Conditions and Farm Management Information. 2nd Edition, PART II/B Central Kenya, Subparts B1 and B2 Central and Rift Valley Province Middle and Ministry of Agriculture, Kenya, in Cooperation with the German Agency for Technical Cooperation (GTZ), Nairobi.
 Othieno, C.O. (1992) Soils. In: Willson, K.C. and Clifford, M.N., Eds., Tea: Cultivation to Consumption, Chapman and Hall, London, 137-171.
 Kebeney, S.J., Kamau, D.M., Othieno, C.O., Ng’etich, W.K. and Owuor, P.O. (2010) Changes in Soil Chemical and Leaf Nutrient Content in Tea Due to Nitrogen Fertiliser Rates and Application Intervals. Tea, 31, 22-37.
 Stephens, W., Othieno, C.O. and Carr, M.K.V. (1992) Climate and Weather Variability at Tea Research Foundation of Kenya. Agricultural and Forest Meteorology, 61, 219-235.
 Williams, E.N.D. (1971) Investigations into Certain Aspects of Water Stress in Tea. In: Carr, M.K.V. and Carr, S., Eds., Water and the Tea Plant, Tea Research Institute of East Africa, Kericho, 79-87.
 Tanton, T.W. (1982) Environmental Factors Affecting the Yield of Tea (Camellia sinensis). II. Effects of Soil Temperature, Day Length and Dry Air. Experimental Agriculture, 18, 53-63.
 Burgess, P.J. and Carr, M.K.V. (1997) Response of Young Tea (Camellia Sinensis) Clones to Drought and Temperature. Experimental Agriculture, 33, 367-383.
 Stephens, W. and Carr, M.K.V. (1990) Seasonal and Clonal Differences Shoot Extension Rates and Numbers in Tea (Camellia sinensis). Experimental Agriculture, 26, 83-98.
 Odhiambo, H.O. (1991) The Effects of Seasonal Changes on Soil Moisture, Atmospheric Humidity, Ambient Temperature and Radiation on Shoot Water Status, Growth and Yield of Four Clones of Tea (Camellia sinensis L.). Msc. Thesis, University of Nairobi, Nairobi.
 Kamunya, S.M., Msomba, S., Makola, L., Korir, R., Kamau, P., Wachira, F.N. and Ndunguru, B.J. (2012) Performance and Genetic Stability for Yield and Quality of Improved Tea Clones in Kenya and Tanzania. Tea, 33, 5-17.
 Tanton, T.W. (1981) Growth and Yield of the Tea Bush. Experimental Agriculture, 17, 177-191.