Research Article | Open Access
Population Viability Analysis of the Endangered Roan Antelope in Ruma National Park, Kenya, and Implications for Management
Population viability analysis (PVA) was used to (1) establish causes of roan population decline for the past 30 years in Ruma National Park (RNP), the only park where wild roans remain in Kenya, and (2) predict the probability of roan persistence under existing and alternative management options. PVA was done using long-term data based on population dynamics, life history, climatic conditions, and expert knowledge. Poaching was identified as the main cause of roan decline in RNP. Several antipoaching and prioritized habitat management interventions to promote population recovery and sustainable conservation of roans are described. PVA predictions indicated that, without these interventions, the roan population cannot persist more than 3 decades. Furthermore, ensuring sustainable conservation of roans in RNP will boost tourism in Western Kenyan and thus alleviate poverty in this part of the country. Improved income from tourism will reduce the possible pressures from hunting and give greater incentives for local people to be actively engaged in roan conservation.
One-third of the total population of roan antelope (Hippotragus equinus) of about 76,000 is now thought to occur in only 4 countries: Burkina Faso, Cameroon, Zambia, and Tanzania . The East African subspecies (H. equinus langheldi) has declined rapidly throughout its range, including Kenya, where it now only survives in Ruma National Park (RNP). Recent roan population decline in RNP from 100 individuals in 1979 to less than 50 currently is of serious concern to the Park management, because a population of less than 50 roan antelopes is not considered viable according to population genetic criteria . Without any interventions, the roan population in RNP will eventually go extinct just as it happened in the other Kenyan protected areas in the past. In fact, the question to ask now is not if but when the extinction will occur. Population viability analysis (PVA) can be used to adequately answer this question and to provide evidence on the likely positive or negative effects of alternative management techniques on roan antelope at RNP.
PVA is a set of modelling techniques that estimates the future size and risk of population extinction [2, 3]. PVA utilizes life history or population growth rate parameters such as survival and fecundity rates as input variables to project dynamics and estimate risk of population extinction . PVA can be used to (i) estimate the probability of extinction ; (ii) predict the future population size ; and (iii) assess risks and benefits of alternative interventions for population recovery .
PVA has been used in the past to model the viability of roan antelope populations and assess consequences of alternative management options . However, some of the findings from these studies are site-specific and cannot be applied generally to areas outside the study area. Hence, further research is needed to devise more effective management interventions for sustainable conservation and management of roan antelopes. This is necessitated by the lack of adequate measures to halt roan population declines in the protected areas previously investigated. The two primary aims of this study were to (i) estimate the future population trends and likelihood of extinction of roan antelopes under the current conditions in RNP and (ii) rank the risks and benefits of management alternatives for roan population recovery.
2. Materials and Methods
2.1. The VORTEX Baseline PVA Model
There are several computer programs available for PVA, but VORTEX was used in this study as it is widely and commonly used. VORTEX is an individual-based model (IBM) that creates a representation of each animal in its memory and follows the fate of the animal through each year of its lifetime . VORTEX keeps track of the sex, age, and parentage of each animal.
Demographic parameters and life history attributes used for the baseline PVA model were based on findings of previous studies on roan antelopes [7–12]. Complementary information about home ranges and age-specific mortality rates were derived from analysis of roan population and distribution data collected in RNP from 1979 to 2008 . A summary of the specific parameter input data used for the baseline PVA model is presented in Table 1.
2.2. Scenarios Modelled
The VORTEX PVA model was used to simulate the population of roan antelopes in Ruma National Park (1) to replicate the observed population decline for past 30 years and (2) to project the future population viability over a 100-year period. In each of these two categories several scenarios were modelled as described below.
2.2.1. Replication of the Observed Population Decline
Simulations of the observed roan population decline for the past 30 years were done to investigate the factors responsible for the decline during that period. The baseline scenario used the roan parameters is described in Section 2.1 above.
All other scenarios were a modification of this baseline scenario. To investigate the effect of inbreeding depression, a model was simulated without the inbreeding depression component. Effects of age-specific mortalities were investigated by simulating several PVA simulation runs using varying mortality rates of 0%, 10%, 20%, 30%, and 40%. Catastrophic effects were studied by considering infrequent severe catastrophes and frequent weak catastrophes. For the infrequent and severe catastrophes scenario, the frequency of each catastrophe in the baseline scenario was halved and the severity was doubled. Likewise, for the frequent and weak catastrophes scenario, the frequency and severity of each catastrophe in the baseline scenario were doubled and halved, respectively. Simulation runs were repeated 1000 times to obtain the most likely population trajectory for each scenario.
2.2.2. Projection of the Future Population Viability
The roan population was simulated to project the population viability for the next 100 years using the baseline scenario and several scenarios of alternative management interventions. Eight management options were investigated: (1) reducing the death of adults; (2) reducing the death of calves; (3) reducing death of subadults; (4) reducing effects of fire; (5) reducing effects of drought; (6) combined interventions consisting of reduction of mortalities and severity of fires and drought; (7) restocking with more roans without changing anything else; and (8) establishing an intensively managed protected sanctuary together with restocking. Again, simulation runs were repeated 1000 times to evaluate the population persistence probabilities for each of these options. Table 2 gives a summary of the demographic parameters used by different scenarios in comparison with the baseline model. Other parameters were the same as those of the baseline model except age-specific mortalities and impacts of catastrophes on drought, fire, and floods.
Adult mortality could be reduced by increasing security patrols, establishing more security outposts, providing adequate suitable habitat through provision of adequate water supplies, controlling the number of competing grazers in the park, and initiating community development projects to alleviate poverty with the aim of reducing poaching for bush meat as a livelihood (management option 1). Calf mortality could be reduced by maintaining conducive habitat for secluding calves from predation by controlled burning (management options 2). Subadult mortality could be reduced by providing adequate suitable habitat for subadults through provision of enough water supplies and controlling the number of competing grazers in the park (management option 3). Some of the competing grazer species include the reedbuck and hartebeest, whose populations have been shown to correlate negatively with roan population in the park [15–17]. The severity and spread of fires could be reduced by implementing prescribed burning programs and increasing road fire breaks (management option 4). The effect of drought (management option 5) could be mitigated by construction of more water dams and water troughs and ensuring that water is pumped to these water points during the periods of drought. The combined intervention scenario (management option 6) assumed that all the above five interventions were to work simultaneously to reduce age-specific mortalities and effects of catastrophic fires, droughts, and floods. Management option 7 assumed running the baseline model with different initial roan population groups from restocking with 4, 6, 8, and 10 roan groups. The typical herd size for a breeding roan group is 12 individuals [10, 15].
The protected sanctuary (management option 8) would promote roan population recovery via several intensive management strategies: (i) prescribed burning; (ii) construction of more road fire breaks; (iii) construction of more water points and ensuring that water is pumped to these water points during the dry season; (iv) putting more effort in controlling poaching by increasing security patrols and establishing more security outposts; (v) limiting the number of competing grazers in the sanctuary; (vi) removal of predators from the sanctuary; (vii) improving the relationship with the surrounding local communities so that they support wildlife conservation; (viii) initiating community development projects to alleviate poverty with the aim of reducing poaching for bush meat as a livelihood; (ix) providing mineral supplements; (x) offering improved veterinary services; (xi) construction of ridges and trenches to control flooding; (xii) ensuring that the sanctuary is fenced and the fence is maintained in a functioning state; and (xiii) restocking with 4, 6, 8, and 10 roan groups. These strategies could increase the carrying capacity from 288 to 640 roans as has been proved in South Africa .
2.2.3. Sensitivity Analysis
The projection of future population viability was repeated to test the sensitivity of age-specific mortalities and catastrophes. A shorter period of 30 years was used as this was considered of relevance to park management objectives at the moment. In the sensitivity analysis, the age-specific mortalities and severity of catastrophic fires and droughts were set from 0% to 40% at intervals of 5%. The population at 30 years for each of these scenarios was analyzed using a generalized linear model (GLM) performed with a Poisson error structure to assess which parameters are most influential in the population viability analysis.
3.1. Replication of the Observed Population Decline
The baseline model derived from the observed roan parameters and those obtained from past studies yielded a population of 43, which is the same as the current population of roans remaining in RNP, but the model was different from the observed trend of population decline for the past 30 years (Figure 1(a)). Nevertheless, it can be used as a basis for projecting future population viability.
(a) Catastrophes and inbreeding
(b) Adult mortality
(c) Subadult mortality
(d) Calf mortality
Models that included inbreeding depression, infrequent severe catastrophes, or frequent weak catastrophes did not differ much from the baseline model (Figure 1(a)). This implies that inbreeding depression and catastrophes were not the major causes of population decline during the past 30 years.
Models with varying age-specific mortality rates seemed to correspond to the observed population in a number of years. The roan population decline from 1979 to 1989 could be attributed to high mortality rates of up to 40% in adults, subadults, or calves. Reduction of adult mortality to 10% and below caused the population to increase accordingly (Figure 1(b)).
However, reduction in mortality rates of subadults and calves to levels as low as 5% did not halt the population decline (Figures 1(c) and 1(d)). This suggests that the observed population decline for the past 30 years may have been primarily controlled by variation in adult mortalities caused by poaching of adult roans and effects of catastrophes such as drought, fire, and floods. Poaching of roans in the study area (RNP) especially using snares has been reported by several researchers over the last 35 years [15, 19–25]. This implies that years with low adult mortality could allow the population to increase and vice versa for years with high adult mortalities. For instance, population increases between 1993 and 2000 and between 2001 and 2004 may reflect a corresponding reduction in adult mortality. In fact the slight increase in roan population observed between 2001 and 2004 was attributed to increased antipoaching security efforts [23, 26].
3.2. Projection of Future Population Viability
Under the current situation (as simulated by the baseline scenario), the roan population has a 100% probability of extinction before 100 years. Extinctions are estimated to begin in 32 years (Table 3). Restocking the park with more roan groups can postpone the median time to extinction up to 48 years but it does not lower the probability of extinction (Figure 2(a), Table 3).
NB: and SD () = population growth rate and its standard deviation; PE (%) = mean probability of extinction; and SD () = mean population size and its standard deviation; TE = median time to extinction, in years.
(a) Restocking with more roan groups
(b) Age-specific mortalities and catastrophes
(c) Protected sanctuary with intensive management
Reducing the rates of calf and subadult mortality to 5% did not have any significant change on the probability of extinction or time to extinction as compared to the baseline scenario (Figure 2(b); Table 3). However, reduction of fire and drought severity could reduce the probability of extinction and enable a persistence probability over 100 years of less than 1% and 3%, respectively (Table 3). Reduction of adult mortality to 10% changed the population growth rate from negative to positive; it caused the population to increase at a rate of 0.1% per year. Also, under this management option, the roan population showed a persistence probability of more than 86% over 100 years (Table 3). Combined management interventions could raise the population growth rate to 2.3% and the persistence probability to 99.8% over 100 years.
Intensive management of the current roan population in a protected sanctuary could raise the population growth rate to 4.3% and allow the population to persist for 100 years with a persistence probability of 100% (Table 3). Under this management option, the population will reach the carrying capacity in 65 years.
Restocking with more roan groups combined with the protected sanctuary intervention did not show any improvement in the overall population growth rate apart from reaching the carrying capacity slightly earlier (Figure 2(c)). However, restocking caused the population to reach high numbers within a short period and to reach the carrying capacity faster. For instance, at 10 years, the simulated mean population size was 130 roans for the management option with current initial population of 3 roan groups (36 individuals) and 370 roans for that with an initial population of 10 roan groups (120 individuals). Also, using the current initial population, the simulated mean population size reached carrying capacity in 65 years whereas that with 10 groups reached the same carrying capacity about 20 years earlier.
3.3. Sensitivity Analysis
Sensitivity analysis was performed for the age-specific mortalities, fire, and drought. Adult mortality was the most sensitive parameter and the most important in determining the change of the roan population (Figure 3). Maintaining adult mortality rates lower than 15% could allow the population to increase and recover. Adult mortality rates higher than 15% could cause the roan population to decrease and consequently go extinct. Maintaining the severity of catastrophic droughts at levels lower than 5% could increase the roan population at a rate lower than that caused by reduction of adult mortalities. The variation in calf and subadult mortalities as well as severity of fires had no substantial impact on the roan population dynamics; whether their effect is eliminated or increased the roan population will eventually decline to extinction.
Further analysis of the impact of these five parameters on roan population dynamics using generalized linear models showed that the adult mortality could account for 78% of the population decline (Table 4). All the other four parameters together accounted for only 18% of the decline. Examination of the regression coefficients (β) showed that drought was the second important factor responsible for roan population changes.
4.1. Replication of the Observed Population Decline
Simulation to replicate the observed roan population decline helped to understand the factors responsible for the decline. Our results indicate that the most important factor affecting the roan antelope population is adult mortality, which confirms existing concerns about the increasing amount of poaching within the park. Records show that poaching in the park was heavy as early as 1970s  and has continued until now [22–25]. The highest rate of population decline from 1979 to 1989 coincided with a period of transition from National Reserve to National Park status . This transition may have caused escalation of poaching by the local community in retaliation to eviction from the park. In Kenya, a national reserve allows limited access of its resources to local communities but national parks are strictly managed for wildlife with total exclusion of humans. Park records also show that the slight increase in roan population observed between 2001 and 2004 was attributed to increased antipoaching security efforts [23, 26]. Past studies have shown that high levels of predation, especially of adults, were responsible for roan population decline in Kruger National Park [7, 9, 25].
4.2. Projection of Future Population Viability
Projections of future roan population viability based on the baseline scenario showed a high probability that at least three decades could pass before extinction eventually occurred, even if no intervention is undertaken. However, this should not send the wrong signal that interventions are not urgently needed since extinctions have been projected to happen in the near future. In Masai Mara National Reserve, a remnant roan population of 45 (±17 SE) animals in 1971 became extinct after about two and half decades . Therefore, the predicted extinction of the RNP remnant roan population of 43 animals in about three decades seems realistic. Hence, this period should be used to implement management interventions.
The projections of future population viability indicated that various management interventions can halt the population decline and cause the population to recovery to healthy levels and probably to carrying capacity. These include reduction of adult mortality, combined management interventions, and establishment of an intensively managed protected roan sanctuary. Reduction of adult mortality appears to be the most crucial intervention without which the roan population is destined to go extinct in the future. Measures for reducing adult mortality were included as part of the combined interventions and protected sanctuary management options. Adult mortality needs to be maintained at levels far below 15% for the roan population to recover. Studies in South Africa have found out that adult mortality rates of 15% or more severely restricts the recovery potential of roan antelopes .
Unfortunately, some management interventions such as prescribed burning, construction of more water points, and restocking under the current conditions could be a waste of time, effort, and money. Although these interventions are still important and necessary they can only be effective when combined with measures of reducing adult mortality.
4.3. Management and Conservation Implications
The next crucial stage after identification of the cause of roan population decline in RNP is to implement effective strategies that can halt the decline and propel the population to recovery. The high adult mortality rate on roan antelopes in Kenya is most likely caused by poaching [15, 19–24] and not by predation. There are no lions in RNP and the only potential predators are hyenas, which may prey on calves and not adult roans . Therefore, the best interventions in this case should include methods of reducing adult mortality by controlling poaching in RNP. Such methods can involve increasing security patrols, opening more security outposts, involving the surrounding communities in management of the park, and initiating community development projects to alleviate poverty with the aim of reducing poaching for bush meat as a livelihood. Involvement of local communities has proven to be effective in curbing elephant poaching in Tanzania’s Selous ecosystem , reducing rhino poaching in Kenya’s Lewa Conservancy , enhancing conservation of many wildlife species in Kenya’s Nakuru Wildlife Conservancy , and reducing human-wildlife conflicts in Kenya’s nonprotected areas .
Apart from adult mortality, the decline in roan population was also contributed by other factors such as drought, fire, and calf and subadult mortality [16, 34, 35]. Therefore, for sustainable conservation of roans in RNP to be achieved, a combined management intervention seems to be the way forward. Such intervention could include (i) measures of reducing adult mortality, (ii) establishment of more water points to alleviate effects of drought, (iii) prescribed burning to control the effects of fires, and (iv) more active manipulative management of the roan habitat. For instance, while reducing adult mortality to 10% could increase the current roan population in RNP to about 100 individuals in 30 years, employing a combined management option could propel the population to about 250 individuals over the same period.
With the current roan habitat conditions coupled with the observed habitat decline over the past 30 years, interventions in RNP under the prevailing management regime may not manage to propel the roan population to numbers that can be considered viable in isolation. Furthermore, putting all efforts of population recovery of a critically locally endangered species in one isolated population is too risky. There is need to establish other roan populations in Kenyan protected areas, where they became locally extinct in the past. Establishing an intensively managed protected sanctuary for roan antelopes in RNP can provide a “seed population” that can be used to establish several other populations in the roan’s former known ranges in Kenya. This may offer a more lasting solution to the problem of roan population decline in Kenya. This concept of breeding roan antelopes in an enclosed sanctuary has already proved to be very successful in other countries . In South Africa, 7 roans were placed in a fenced enclosure within the roan range in 1994 and the population grew to 41 roans by 2001 .
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
The author acknowledges the Commonwealth Scholarship Commission for funding this research and appreciates the management of Ruma National Park (RNP) for providing secondary data for this research.
- IUCN, IUCN Red List of Threatened Species, 2017, http://www.iucnredlist.org/.
- M. E. Gilpin and M. E. Soule, “Minimum viable populations: processes of species extinction,” in Conservation Biology: The Science of Scarcity and Diversity, pp. 19–34, Sinauer Associates, Sunderland, Mass, USA, 1987.
- R. Frankham, C. J. A. Bradshaw, and B. W. Brook, “Genetics in conservation management: Revised recommendations for the 50/500 rules, Red List criteria and population viability analyses,” Biological Conservation, vol. 170, pp. 56–63, 2014.
- D. Ludwig, “Is it meaningful to estimate a probability of extinction?” Ecology, vol. 80, no. 1, pp. 298–310, 1999.
- P. S. Miller and R. C. Lacy, VORTEX: A Stochastic Simulation of the Extinction Process. Version 9.50 User's Manual, Conservation Breeding Specialist Group (SSC/IUCN), Apple Valley, Minn, USA, 2005.
- S. R. Beissinger and D. R. McCullough, Population Viability Analysis, University of Chicago Press, Chicago, USA, 2002.
- C. A. McLoughlin and N. Owen-Smith, “Viability of a diminishing roan antelope population: Predation is the threat,” Animal Conservation, vol. 6, no. 3, pp. 231–236, 2003.
- J. K. Kimanzi, Mapping and modelling the population and habitat of the roan antelope (Hippotragus equinus langheldi) in Ruma National Park, Kenya [PhD Dissertation], Newcastle University, UK, 2011, https://theses.ncl.ac.uk/dspace/bitstream/Kimanzi/.
- R. Harrington, N. Owen-Smith, P. C. Viljoen, H. C. Biggs, D. R. Mason, and P. Funston, “Establishing the causes of the roan antelope decline in the Kruger National Park, South Africa,” Biological Conservation, vol. 90, no. 1, pp. 69–78, 1999.
- J. D. Skinner and R. H. N. Smithers, The Mammals of the Southern African Subregion, University of Pretoria, Pretoria, South Africa, 1990.
- S. C. J. Joubert, The population ecology of the roan antelope, Hippotragus equinus equinus (Desmarest, 1804), in the Kruger National Park, University of Pretoria, 1976.
- P. Pinto, P. Beja, N. Ferrand, and R. Godinho, “Hybridization following population collapse in a critically endangered antelope,” Scientific Reports, vol. 6, Article ID 18788, pp. 1–9, 2016.
- R. M. Poche, “Notes on the roan antelope (Hippotragus equinus) in West Africa,” Journal of Applied Ecology, vol. 11, no. 3, pp. 963–968, 1974.
- K. Ralls, J. D. Ballou, and A. Templeton, “Estimates of lethal equivalents and the cost of inbreeding in mammals,” Conservation Biology, vol. 2, no. 2, pp. 185–193, 1988.
- J. K. Kimanzi, Roan Antelope Population and Habitat Evaluation in Ruma National Park: Implications for Management, Lambert Academic Publishing, Saarbrucken, Germany, 2012.
- J. K. Kimanzi, R. A. Sanderson, and S. P. Rushton, “Habitat suitability modelling and implications for management of roan antelopes in Kenya,” African Journal of Ecology, vol. 52, no. 1, pp. 111–121, 2014.
- J. K. Kimanzi, J. N. Wanyingi, H. Ipara, and P. O. Odwori, “Factors influencing population dynamics of five antelope species in Ruma National Park, Kenya,” African Journal of Education, Science and Technology, vol. 3, no. 1, pp. 59–70, 2016.
- W. G. Dörgeloh, W. Van Hoven, and N. F. G. Rethman, “Population growth of roan antelope under different management systems,” African Journal of Wildlife Research, vol. 26, no. 4, pp. 113–116, 1996.
- R. Allsopp, “Roan antelope population in the lambwe valley, Kenya,” The Journal of Applied Ecology, vol. 16, no. 1, pp. 109–115, 1979.
- Youth for Conservation, A Status Report on Snares in Ruma National Park, Youth for Conservation (YFC), Nairobi, Kenya, 2006.
- J. K. Kimanzi and J. N. Wanyingi, “The declining endangered roan antelope population in Kenya: what is the way forward?” Conference Papers in Science, vol. 2014, pp. 1–6, 2014.
- D. Kones, The Status of Kenya’s Roan Antelope (Hippotragus equinus-langheldi) in Ruma National Park 1974–2005, vol. 49, Kenya Widlife Service (KWS), Nairobi, Kenya, 2005.
- J. K. Kimanzi, R. A. Sanderson, S. P. Rushton, and M. J. Mugo, “Spatial distribution of snares in Ruma National Park with implications for management of roan antelope and other wildlife,” Oryx (Fauna and Flora International) Journal, vol. 49, no. 2, pp. 295–302, 2015.
- G. George Ariya, “Wildlife snaring by the local community in Ruma National Park, Kenya: can conservation tourism be an alternative livelihood strategy?” International Journal of Business and Social Science, vol. 6, no. 2, pp. 141–149, 2015.
- C. P. Havemann, T. A. Retief, C. A. Tosh, and P. J. N. de Bruyn, “Roan antelope Hippotragus equinus in Africa: A review of abundance, threats and ecology,” Mammal Review, vol. 46, no. 2, pp. 144–158, 2016.
- J. Asila, Ruma National Park Animal Census Research Report, Kenya Wildlife Service (KWS), Nairobi, Kenya, 2004.
- KWS, A policy framework and development programme 1991–1996, Kenya Wildlife Service (KWS), Nairobi, Kenya, 1990.
- M. D. Broten and M. Y. Said, “Population Trends of Ungulates in and around Kenya's Masai Mara Reserve,” in Serengeti II: Dynamics, Management, and Conservation of an Ecosystem, A. R. E. Sinclair and P. Arcese, Eds., pp. 169–193, The University of Chicago Press, Chicago, USA, 1995.
- F. K. Waweru, C. Alina, and C. M. Musyoki, Rhino Habitat Suitability Assessment in Ruma National Park, Kenya Wildlife Service, Nairobi, Kenya, 1995.
- W. Lotter and K. Clark, “Community involvement and joint operations aid effective anti-poaching in Tanzania,” PARKS, vol. 20, no. 1, pp. 19–29, 2014.
- Daily Nation, “Report shows community involvement helping reduce human-wildlife conflict,” http://www.nation.co.ke/counties/Laikipia/community-wildlife-conflict-Lewa-report/1183290-2694848-jq9msu/index.html.
- J. O. Ogutu, B. Kuloba, H.-P. Piepho, and E. Kanga, “Wildlife population dynamics in human-dominated landscapes under community-based conservation: The example of Nakuru Wildlife Conservancy, Kenya,” PLoS ONE, vol. 12, no. 1, Article ID e0169730, 2017.
- D. Western, J. Waithaka, and J. Kamanga, “Finding space for wildlife beyond national parks and reducing conflict through community based conservation: the Kenya experience,” PARKS, vol. 21, no. 1, pp. 51–62, 2015.
- J. O. Ogutu, H.-P. Piepho, M. Y. Said et al., “Extreme wildlife declines and concurrent increase in livestock numbers in Kenya: What are the causes?” PLoS ONE, vol. 11, no. 9, Article ID e0163249, 2016.
- E. Gandiwa, I. M. A. Heittkönig, P. H. C. Eilers, and H. H. T. Prins, “Rainfall variability and its impact on large mammal populations in a complex of semi-arid African savanna protected areas,” Tropical Ecology, vol. 57, no. 2, pp. 163–180, 2016.
Copyright © 2018 Johnstone K. Kimanzi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.