Cassava Yield, Garrification Traits and Economic Returns as Influenced by Different Genotypes in an Ultisol

Abstract

Field trials were conducted at the University of Uyo Teaching and Research Farm in 2009/2010 and 2010/2011 to evaluate cassava yield, garrification traits, and economic returns as influenced by genotypes grown in an ultisol. Uyo is located between latitudes 4°30' and 5°27'N and longitudes 7°50'E and 80°20'E. The area, which lies within the humid tropical rainforest zone of southeastern Nigeria, has an annual mean rainfall of 2500 mm with a bimodal rainfall pattern, monthly sunshine of 3.14 hours, and a mean annual temperature of 28°C. The experimental site was previously put into cultivation of some arable crops such as fluted pumpkin, okra, and waterleaf before it was fallowed for two years. A randomized complete block design, replicated three times was used. Treatments were 16 cassava genotypes (AR1-82, NR02/0028, NR03/0174, CR12-45, NR03/0211, TMS98/2132, TMS01/1206, TMS01/1368, TMS01/1371, CR36-5, NR02/0007, NR03/0155, AR37-108, TMS01/1412, TMS30572 and local best variety “Obubit okpo”). Results obtained via combined (pooled) analysis of the two seasons indicated that cassava yield differed significantly among the genotypes. NR02/0018 (42.50 t/ha) out-yielded others by 17% - 52% followed by NR 02/0007 (35.42 t/ha), and NR 03/0155 (32.08 t/ha). TMS 98/2132 had the most preferred garrification attributes followed by NR 02/0018, NR 02/0007, and NR03/01155. NR 02/0018 genotype gave the highest economic returns followed by NR 02/0007 and NR03/0155. NR02/0018 superseded others in yield, garrification traits, and economic returns, which is therefore recommended.

Share and Cite:

Ndaeyo, N. , Okon, J. , Egesi, C. , Eke-Okoro, O. and Ogbe, F. (2022) Cassava Yield, Garrification Traits and Economic Returns as Influenced by Different Genotypes in an Ultisol. American Journal of Plant Sciences, 13, 744-755. doi: 10.4236/ajps.2022.136050.

1. Introduction

Cassava (Manihot esculenta Crantz), a member of the Euphorbiaceae family,is one of the most important food crops in the tropics. It is a staple food for over 50 million Nigerians, and is cultivated on small and large scales, especially in the southern parts of the country. The crop is amendable to agronomic as well as genetic improvement and has a high yield potential under good conditions and performs better than other crops under sub-optimal conditions. Cassava grows over a wide range of ecological conditions and is more tolerant to low soil fertility, drought, pests, and diseases, and can be grown all year round [1] [2]. It tolerates altitudes of up to 2000 m above sea level [1], requires annual rainfall above 750 mm [3], and can survive 4 - 6 months of dry weather. Cassava appears to be generally photoperiodic neutral for commercial cultivation purposes and the suitable temperature requirement range is 24˚C - 30˚C. It derives its importance from its starchy, thickened, storage roots which are valuable sources of cheap calories, especially in developing countries where calories deficiency and malnutrition are widespread [1] [4]. Over two-thirds of the total production of cassava is consumed in various forms by humans. Cassava storage roots can be processed into garri, fufu, chips, starch, tapioca, and flour among others [5] [6].

Over the years, the International Institute of Tropical Agriculture (IITA), Ibadan, Nigeria, and National Root Crop Research Institute, Umudike, Nigeria have developed improved and new high-yielding cassava varieties under the Tropical Manihot Selection (TMS) and National Root (NR) Series [7] [8]. These improved varieties are said to yield between 25 and 40 tonnes per hectare, are resistant to pests and diseases, and have acceptable sensory/ culinary and industrial qualities [9] [10] compared to local varieties with yields between 5 and 10 t/ha and even those earlier released. Efforts are still continuing in the areas of breeding for specific ecologies, low cyanide content, specific industrial requirements, and early bulking. The desired attributes of a good cassava variety are basically early bulking, high dry matter and starch contents, resistance to pests and diseases, low cyanide content, weed-suppressing ability, good in-ground storability and good processing, and cooking qualities [11].

Garri, a fermented gelatinous granular flour is a major source of dietary energy for low-income consumers in many parts of Tropical Africa, including major urban [12] [13]. It is eaten in Nigeria and other parts of West Africa. However, acceptability relies on the conformity with the major quality attributes demanded in garri as affected by the cassava cultivars or varieties [14]. Garri colour, texture, and overall garri qualities acceptability such as swelling index, suitability for Sipping and making into “eba” (a native bolus) were reported to be controlled by the duration of fermentation and methods of processing and not on cultivars [15].

Although Nigeria is currently the world-leading producer of cassava, low yields are still being recorded by cassava farmers, particularly at the peasant level. This has been attributed to several factors such as the inability of farmers to select and/or cultivate improved variety, pest and disease attacks, availability of adaptable and appropriate varieties, as well as low soil fertility status. Besides these factors, there is also the challenge of having cassava varieties that can meet the desired qualities of the various end-users. Against these backgrounds, a study was conducted to evaluate the yield and yield components, garrification attributes, cost of production, and economic returns (profit) of some newly bred cassava genotypes under ultisol in a rainforest ecology.

2. Materials and Methods

2.1. Experimental Site and Cropping History

This study was conducted at the University of Uyo Teaching and Research Farm located at Use Offot, Uyo, during the 2009/2010 and 2010/2011 planting seasons. Uyo is located between latitudes 4˚30' and 5˚27'N and longitudes 7˚50'E and 80˚20'E [16]. The area, which lies within the humid tropical rainforest zone of southeastern Nigeria, has an annual mean rainfall of 2500 mm with a bimodal rainfall pattern. The rain usually starts in March and ends in November with a short period of relative moisture stress in August, traditionally referred to as “August Break” [17]. The average monthly sunshine is 3.14 hours and annual temperature is 28˚C. The temperature is generally high in the months of February through April [18]. Uyo has an annual average relative humidity of 79% and evaporation rate of 2.6 cm2. The soil is an ultisol which is moist, highly weathered and acidic but good enough to support cassava production. The experimental site was previously put into cultivation of some arable crops such as fluted pumpkin (Telfairia occidentalis Hook. F.), okra (Abelmoschus esculentus L. Moench), and waterleaf (Talinum triangulare (L.) Juss) before it was fallowed for two years.

2.2. Experimental Design, Treatment, and Cultural Details

The experiment site was laid out in a randomized complete block design with three replicates. The entire experimental size measured 45 m × 25 m. Each replicate measured 45 m × 4 m and was demarcated from each other by 1 m path. Each plot measured 4 m × 4 m and demarcated from the other by a 1 m path. The treatments were 15 cassava genotypes (AR 1-82, NR 02/0018, NR 03/0174, CR 12-45, NR 03/0211, TMS 98/2132, TMS 01/1368, TMS 01/1371, CR 36-5, NR 03/0155, NR 02/0007, AR 37-108, TMS 01/1206, TMS 01/1412, and TMS 30572) obtained from National Root Crops Research Institute, Umudike, Abia State, Nigeria and 1 local best (Obubit okpo) obtained from Use Offot, Uyo, Nigeria.

The experimental site was mechanically ploughed, harrowed and ridged with tractor mounted implements (plough, harrower, ridger, respectively), before it was marked out using measuring tape, rope and pegs. Planting was done on 12th June in both seasons using cassava stem cuttings measuring 25 cm and inserted on the crest of the ridge in a vertical position at a spacing of 1 m × 1 m (10,000 stands per hectare). Weeding was done three times using native weeding hoe at 3, 6 and 9 months after planting (MAP). NPK (15:15:15) fertilizer was manually applied to all plots using ring method at the rate of 400 kg·ha−1. The cassava roots were harvested manually at 12 MAP in line with local farmers’ practice. Yield and yield components determined were: number of stands per hectare (counted at harvest), number of rotten tubers (counted per plot at harvest), number of fresh tubers per plant (counted at harvest), tuber length (measured in cm from proximal to distal end) and circumference (measured in cm at the middle of each tuber), tuber yield per hectare (determined using a top-load weighing balance) and dry matter accumulation in the tubers (drying with oven and weighing with a top-load weighing balance). Cassava stems were harvested at soil surface level at harvest and 50 good stems, each at least of 1 m length constituted a bundle.

2.3. Cassava Processing and Garrification

The tubers (20 kg from each genotype) were manually peeled after harvesting using kitchen knife. The peeled tubers were washed with water and grated using mechanical grater. The grated cassava mash was packed into a white porous sack bags and pressed using hydraulic press (5 hp diesel powered engine) for two days to dewater the pulp. The dewatered lump of the pulp was manually pulverized/crushed and sieved. The resulting semi-dry fine pulp was fried in a shallow iron cast fryer using farmer’s method. During the frying process, a piece of calabash was used in toasting the particles within the fryer. The dried garri particles were cooled and arranged on genotype basis by spreading them on different clean polythene sheets. They were then packed into different polythene bags, weighed and labeled for easy identification [19]. Sensory evaluation of the garri from the genotypes by consumers was done in Department of Crop Science laboratory. A taste panel comprising ten persons inspected and ranked garri from the sixteen genotypes based on preference for texture, sipping, making of “eba” (a local bolus) and swelling index. Garri (1 kg) from each genotype was poured into sixteen separate plates of the same size for the taste panel members to examine by sighting of the garri followed by feeling the texture of the samples with their hands. Garri samples were ranked on a scale of 1 - 6, where 1 = very fine, 2 = fine, 3 = fairly fine, 4 = fairly rough, 5 = rough and 6 = very rough [11] for the traits. The same scale was also used in ranking preference based on texture, and sipping. Sipping of garri was done after 60 cl of cold water was added to 30 g of garri from each of the genotypes. Preference based on making 40 g of garri into“eba” from each of the genotypes was done by pouring 70 cl of hot water to it and rated using the same scale [11]. Swelling index was obtained by measuring 100 g of garri from each of the genotypes into a well calibrated glass cylinder. Cold water was poured into each cylinder to a level of 1000 ml. The swelling index was observed at 5, 10, 15, 20, 25, 30 and 60 minutes [20] which was later converted and expressed in percent.

2.4. Data Analysis

All yield data collected from the two seasons were combined (pooled) for analysis and subjected to analysis of variance procedure and treatments means that indicated significant difference were compared using Duncans Multiple Range Test (DMRT) at 5 percent level of probability [21].

2.5. Cost-Benefit Analysis

Cost-benefit of producing fresh cassava tubers, stem cuttings and garri as affected by the different genotypes were analysed using partial budgeting to determine economic returns [22].

3. Results

3.1. Yield and Yield Components of Cassava as Influenced by Different Genotypes

Table 1 shows yield and yield components of cassava as influenced by different genotypes. The number of stands per hectare at harvest as influenced by different genotypes was not significantly different among the genotypes. However, the highest number of stands per hectare on the average in both seasons was from TMS98/2132 (9792) and the CR36-5 (9792) genotypes while the least was from AR1-82 (8542). The TMS98/2132 and CR36-5 had 2% - 11% more number of stands per hectare than other genotypes. There were significant differences (P ≤ 0.05) in the number of cassava fresh tubers per plant among the genotypes in both seasons (Table 1). The CR36-5 and NR03/0155 had the highest number of fresh tubers per plant (7.17) whereas AR37-108 recorded the least number (4.17). The CR37-5 and NR03/0155 genotypes had 13% - 33% more fresh tubers than other genotypes in both seasons.

There was no significant difference in the length of fresh tubers (cm) as influenced by different genotypes (Table 1) but CR12-45 genotype had the longest fresh tubers (47.67 cm) on the average while the shortest fresh tubers were from AR1-82 (37.78 cm). There was significant difference (P ≤ 0.05) in circumference of fresh tubers as influenced by the different genotypes (Table 1). The TMS 01/1368 has the widest fresh tubers (18.33 cm) whereas the least fresh tubers circumference was from CR12-45 (12.22 cm). The TMS01/1368 had 2% - 32% wider fresh tuber circumference than other cassava genotypes.

Table 1 indicates that there were significant differences (P ≤ 0.05) in the number of rotten tubers per plot as influenced by different genotypes. TMS98/ 2132, NR02/0007 and AR37-108 genotypes had the least number of rotten tubers (0.33), whereas the highest was from TMS01/1371 (3.00). There were significant differences (P ≤ 0.05) in tuber yield as influenced by different genotypes. The NR02/0018 had the highest fresh tuber yield (42.50 t/ha) whereas the local best variety “Obubit okpo” had the least yield (17.00 t/ha). The NR02/0018 genotype out-yielded other cassava genotypes by 16% - 52%.

Similarly, Table 1 shows that the NR02/0018 had the highest garri yield (15.96 t/ha) while the least garri was from the local variety “Obubit okpo” (3.71 t/ha). The NR02/0018 had 24% - 70% more garri yield than others. Table 1 and Table 2 also show that cassava dry matter accumulation as influenced by different genotypes was significantly different (P ≤ 0.05). Dry matter accumulated most in the local variety “Obubit okpo” (350.67 g/kg) while TMS 01/1412 had the least (210.00 g/kg). The local variety, “Obubit okpo”, had 3% - 38% higher dry matter accumulation than other cassava genotypes.

3.2. Garrification Attributes of Cassava as Influenced by Different Genotypes

Table 2 shows garri swelling characteristics of cassava as influenced by different genotypes after soaking in water for 5, 10, 15, 20, 25, 30 and 60 minutes. The NR02/0018 genotype had the highest swelling index of 190, 340, 360, 390, 400 and 410 ml at, 10, 15, 20, 25, 30 and 60 minutes, respectively (which was later converted and expressed in percent) whereas the least swelling was from TMS01/1368. The NR02/0018 and NR03/0174 genotypes had the highest swelling index (293.3%). The swelling index of NR02/0018 and NR03/0174 was 1% - 35% more than those of other cassava genotypes.

Table 1. Yield and yield components of cassava as influenced by different genotypes (means of two seasons).

a, b, c mean the same superscript in each column is not significantly different (P < 0.05) according to DMRT.

Table 2. Garri swelling characteristics (ml) of cassava as influenced by different genotypes (means of two seasons).

Table 3 shows other sensory evaluation attributes of the genotypes. The TMS98/2132 genotype was the most preferred by sight (1.4 score—very fine) and NR03/0174 was the least preferred. On preference based on texture, TMS98/2132 had the finest texture (1.4 score—very fine) while CR12-45 had the worst texture (6.0—very rough).

The TMS98/2132 genotype was the most preferred in sipping (1.6—very fine) while the local best variety “Obubit okpo” was the least preferred (6.0—very rough). Table 3 further shows that TMS98/2132 was the most preferred when made into “eba” and moulded into garri balls with a score of 1.2 (very fine) while TMS01/1412 was the least preferred (6.0—very rough). On the basis of overall acceptability, TMS98/2132 had the most preferred garrification attributes followed by NR02/0018 and NR03/0155 whereas TMS01/1368 genotype had the least acceptance.

3.3. Cost of Production and Economic Return of Management (N /Ha) as Influenced by Different Genotypes

Table 4 indicates that the total cost of production per hectare for fresh tuber production for the cassava genotypes was NI53,000, while for the production of garri an extra N31,250 was required for transportation of the fresh tubers to processing points and final processing of the fresh tubers to garri, all totaling N184,250. The returns to management or net profit showed that NR02/0018 genotype had the highest economic returns (N4,977,000) while the least was from the local best variety (Obubit okpo) on the basis of fresh cassava tubers and stems. When the cassava genotypes are considered on the basis of fresh tuber yield alone, NR02/0018 produced the highest net profit (N4,947,000) and the least was from the local variety (Nl,887,000).

Table 3. Some garrification attributes of cassava as influenced by different genotypes (means of two seasons).

Scale: 1 = Very fine, 2 = Fine, 3 = Fairly fine, 4 = Fairly rough, 5 = Rough, 6 = Very rough.

Table 4. Cost of production and economic return of management (N/ha) as influenced by different genotypes (means of two season.

A yield × mean unit price of N120000.00 per tonne of fresh cassava tuber, N1100.00 per bundle of cassava stem cuttings sand N360000.00 per tonne of garri based on the prevailing market price at the time of harvest.

On the basis of garri and stem production, NR 02/0018 had the highest returns (N5,591,350) whereas the least was from the local variety (NI,180,150). When the genotypes were considered on the basis of garri alone, NR02/0018 equally had the highest economic returns (N5,561,350) whereas the local variety had the least returns (NI,151,350). On the basis of stem yield alone, there was no profit, though NR02/0018 produced the highest stem yield (300 bundle/ha), while the least bundle per hectare was from CR12-45 (120 bundle/ha).

4. Discussion

The results of the swelling index (percentage) of the garri samples showed that swelling capacity is apparently a function of genetic makeup of the different genotypes which perhaps varied in such traits as the starch content and the degree of gelatinization. This reflected in the swelling properties of the garri and was expressed during sipping and making into eba. Sanni et al. (2005) [23] reported that the swelling index of granules reflect the extent of associative forces within the granules, therefore the higher the swelling index, the lower the associative forces.

It has been reported that a good quality garri should have swelling capacity value above three times its original volume [24] [25] [26]. This does not however suggest that those whose swelling index is below three times its original volume may not be cherished for other attributes.

Cost of production and economic returns to management show promising returns to investment on the basis of tuber and garri yields. This is because all the improved varieties had high tuber yielding capacity than the local variety. This is apparently because the improved cassava genotypes are bred to withstand adverse conditions and yet perform considerably well including resistance and tolerance to pests and diseases [2]. However, in terms of stem cutting yield, returns to management were poor, probably because of the prevailing market price of cassava stem cuttings at the time of harvest, which was a period of low demand. However, processing cassava into garri produced the highest returns to investment, with NR02/0018 giving the highest returns, followed by NR02/0007 and NR03/0155 while the local variety produced the least return to management. Eke-Okoro et al. (2005) [27] had reported a very high-profit margin from cassava production and the findings from this study are in tandem with the findings of Ikeh [28].

5. Conclusion

Based on the findings of this study, it is concluded that planting NR02/0018 genotype would produce higher total tuber yield and economic returns and thus be more beneficial to the farmers. On the hand, if the interest of the farmer is on garrification, the TMS 98/2132 ranked highest in garrification overall attributes and as such is preferred followed by NR02/0018, NR02/0007, and NR03/0155. Therefore, the two genotypes (NR02/0018 and TMS 98/2132) which showed superiority in tuber yield, garrification attributes, and economic returns are recommended to the farmers in Uyo, Southeastern Nigeria.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Opeke, L.K. (2006) Essentials of Crop Farming. Spectrum Books Limited, Spectrum House, Ibadan, 140-145.
[2] Nnodu, E.C., Ezulike, T.O. and Asumugha, G.N. (2006) Cassava. In: Idem, N.U.A. and Showemimo, F.A., Eds., Tuber and Fibre Crops of Nigeria: Principles of Production and Utilization, Ade Commercial Press, Zaria, 22-44.
[3] Udoh, D.J., Ndon, B.A., Asuquo, P.E. and Ndaeyo, N.U. (2005) Crop Production for the Tropics. Concept Publishers, Lagos, 442 p.
[4] International Institute of Tropical Agriculture (IITA) (2006) Cassava Recipes for Household Food in Security Ibadan, Nigeria.
https://www.iita.org/
[5] International Institute of Tropical Agriculture (IITA) (1990) Cassava in Tropical Africa: A Reference Manual. IITA, Ibadan, 1-176.
[6] Agwunobi, L.N. and Agbakoba, A. (1990) The Use of Cassava in the Diets of Monogastric Animals. Proceedings 1990 NAERLS National Home Economics Workshop, NRCI, Umudike, 9-13 July 1990, 43-46.
[7] Bokanga, A. and Otoo, E. (1991) Cassava Based Food, How Safe Are They? Tropical Root Crops in a Developing Economy, Proceedings of the 9th Symposium of the International Society for Tropical Root Crops, Accra, 20-26 October 1991, 225-232.
[8] Okogbenin, E., Egesi, C., Mba, C., Espinosa, E., Gullermo, L.S., Ospina, C., Marn, J., Barrera, E., Gutierrez, J. and Fregene, M. (2009) Bridging the Cassava Productivity Gap between Africa and the World: Identification of Elite CIAT Lines Adapted to Nigeria Agro-Ecologies. NRCRI, Umudike, 3-5.
[9] Breckelbaum, T., Belotti, A. and Lozano, J.C. (1978) Cassava Protection Workshop. CIAT Cali, Colombia, 7-12.
[10] International Institute of Tropical Agriculture (IITA) Annual Reports 1976-1996, Ibadan, Nigeria.
[11] Udoh, A.J. and Ndaeyo, N.U. (2000) Crop Productivity and Land Use Efficiency in Cassava-Maize Systems as Influenced by Cowpea and Melon Population. Tropical Agriculture, 77, 44-48.
[12] Ihekoronye, A.I. and Ngoddy, P.O. (1985) Integrated Food Science and Technology for the Tropics. Macmillan Press, London, 189, 241-379.
[13] Makanjuola, O.M., Ogunmodede, A.S., Makanjuola, J.O. and Awonorin, S.O. (2012) Comparative Study on Quality Attributes of Gari Obtained from Some Processing Centres in South West, Nigeria. Advance Journal of Food Science and Technology, 4, 135-140.
[14] Olakunle, M.M., Akinwale, S.O., John, O.M. and Samuel, O.A. (2012) Comparative Study on Quality Attributes of Garri Obtained from Some Processing Centres in South West, Nigeria. Advance Journal of Food Science and Technology, 4, 135-140.
[15] Irtwange, S.V., and Achimba, O. (2009) Effect of the Duration of Fermentation on the Quality of Garri. Current Research Journal of Biological Sciences, 1, 150-154.
[16] UCCDA (Uyo Capital City Development Authority) 1988.
[17] Peters, S.W., Usoro, E.J., Udo, E.J., Obot, U.W. and Okpon, S.N. (1989) Akwa Ibom State: Physical Background, Soils and Land Use and Ecological Problems. A Technical Report of the Task Force on Soils and Land Use Survey, Akwa Ibom State, 189-191.
[18] Enwezor, W.O., Udo, E.J., Ayotode, K.A., Adepetu, J.A. and Chude, V.O. (1990) A Review of Soil and Fertilizer Use Research in Nigeria. Vol. 4, Middle Belt Zone, Federal Ministry of Agriculture and Natural Resources, Abuja, 217 p.
[19] Ndirika, V.I. (2006) Post-Harvest and Processing Technologies of Root and Tuber Crops. In: Idem, N.U.A. and Showemimo, F.A., Eds., Tuber and Fibre Crops of Nigeria: Principles of Production and Utilisation, Ade Commercial Press, Zaria, 196-207.
[20] Ekpe, E.O. (1998) Production of Cassava (Manihot esculenta Crantz) Morphotypes in Cassava-Maize Mixture in the Humid Tropics. Ph.D. Thesis, Federal University of Technology, Owerri, 184 p.
[21] Wahua, T.A. (1999) Applied Statistics for Scientific Studies. Africa-Link Books, Ibadan, 72-25.
[22] Ndaeyo, N.U. and Aiyelari, E.A. (2010) Yields, Cost of Production and Economic Returns to Management of Maize/Cassava Intercrop as Influenced by Different Tillage Practices. Researcher, 2, 68-74.
[23] Sanni, L., Maziya-Dixon, B., Akanya, J., Okoro,Y., Egwuonwu, C.V., Okechukwu, R., Ezedinma, C., Akoroda, M., Lemchi, J., Okoro, E. and Dixon, A. (2005) Standards for Cassava Products and Guidelines for Export. IITA, Ibadan, 93 p.
[24] Ingram (1975) Standards, Specification and Quality Requirements for Processed Cassava Products. Reports G102. Tropical Product Institute, London.
[25] Ukpabi, U.J. and Ndimele, C. (1990). Evaluation of the Quality of Gari Produced in Imo State. Nigerian Food Journal, 8, 105-110.
[26] Sanni, L.O., Adebowale, A.A., Awoyale, W. and Fetuga, G.O. (2008) Quality of Gari (Roasted Cassava Mash) in Lagos State, Nigeria. Nigerian Food Journal, 26, 125-134.
https://doi.org/10.4314/nifoj.v26i2.47446
[27] Eke-Okoro, O.N., Ekwe, K.C. and Nwosu, K.I. (2005) Cassava Stem and Root Production. A Practical Manual. National Root Crops Research Institute, Umudike, 22-30.
[28] Ikeh, A.O. (2017) Effects of Time of Harvesting and Fertilization on Yield, Garrification Traits and Economic Returns of Cassava (Manihot esculenta Crantz) Genotypes. A Doctoral Thesis, Postgraduate School, University of Uyo, Uyo.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.