SCIENTIFIC REPORTS

<figure><img src="image_1.png" /></figure>

OPEN

# Wild-captive interactions and economics drive dynamics of Asian elephants in Laos

Received: 27 March 2017

Accepted: 25 September 2017

Published online: 01 November 2017

Gilles Maurer<sup>1,2</sup>, Benjamin S. Rashford<sup>3</sup>, Vatsana Chanthavong<sup>4,5</sup>, Baptiste Mulot<sup>2</sup> & Olivier Gimenez<sup>1</sup>

The interactions between wild and captive populations of Asian elephants (*Elephas maximus*) persist in most countries of the species distribution, notably through the reproduction between captive females and wild males. However, these complex interactions have been poorly studied, despite their relevance for conservation of this endangered species. Laos has a centuries-long tradition of raising Asian elephants. Besides being cultural icons, captive elephants are inextricably linked to economics through their work in forestry. Using an ecological-economic model, we investigated the effect of socio-economic strategies on fecundity of the Lao population whose dynamics is shaped by human practices. We demonstrated that fecundity is impacted by: i) the dynamics of the wild elephant pool through mating of captive females by wild males, and ii) the financial incentive of elephant owners to breed their animals. As a result, we expect fecundity to rise in response to increases in elephant prices. The captive population will tend towards an asymptotic limit determined by the wild pool growth rate. However, the population will tend to extinction if exports continue. Our ecological-economic approach, by accounting for economic incentives, allows us to predict new equilibria that can serve as a baseline for designing sustainable management strategies for the species.

The population of Asian elephants, *Elephas maximus*, is in overall decline and classified as endangered by the International Union for Conservation of Nature. Total population abundance ranges from 30,000 to 50,000 individuals; however, these figures are probably overestimates<fnref n="1" />⁻³. One-quarter to one-third of the population consists of captive elephants⁴,⁵. Captive elephants were traditionally sourced from wild populations though bans on wild capture have been in place since the 1970s⁶⁻⁸. These so-called "domestic" elephants are, however, not a domesticated species and are similar to their wild counterparts in terms of genetics and behaviour⁹.

The interactions between wild and captive populations persist in most countries of the species distribution, mainly through reproduction between captive females and wild males. In Laos, India and Myanmar, captive females are released in the forest and pair with wild males⁹⁻¹². Nearly 80% of calves born in captivity in Laos during the past decade are from wild genitors from the Nam Pouy Protected Area (V. Chanthavong, pers. comm.). In Myanmar, live capture of wild elephants has been used to sustain the captive pool. This practice may, however, pose a serious threat to the survival of the country's wild population⁴. Elephant calves also have great appeal to tourists and, hence, are highly valued¹³. With corruption, a lack of law enforcement, and deficient registration systems of captive individuals in many Asian countries, the incentive is high to illegally capture and smuggle calves to supply the tourism sector¹⁴. As a result, the captive population partially relies on wild males for breeding, while the wild pool is threatened by illegal captures to supplement the captive pool. Both captive and wild pools also share common threats and conservation challenges, such as habitat destruction, poaching, increased conflicts with local human populations, and risks of inbreeding depression¹². Despite their relevance for the conservation of the species, the complex interactions between wild and captive pools have been poorly studied.

The understanding of Asian elephant demography and ecology has recently improved thanks to a longitudinal study of a large semi-captive population in Myanmar. For example, it has been shown how life histories could

<footnote n="1">¹Centre d'Ecologie Fonctionnelle et Evolutive (UMR 5175), CNRS – Université de Montpellier – Université Paul-Valéry Montpellier - EPHE, campus CNRS, Montpellier, France. ²Zooparc de Beauval & Beauval Nature, Saint-Aignan, France. ³Department of Agricultural and Applied Economics, University of Wyoming, Laramie, USA. ⁴Department of Livestock and Fisheries, Ministry of Agriculture and Forestry, Vientiane, Lao PDR. ⁵ElefantAsia non-profit organization, Paris, France. Correspondence and requests for materials should be addressed to G.M. (email: gilles.maurer@gmail.com)</footnote>
impact reproductive success, survival or senescence<sup>15-18</sup>. Climate variations also impact age-specific survival of Asian elephants<sup>19,20</sup>. Besides survival, large variations in average fecundity rates have been observed in different populations across the region and over time, with a fecundity rate of 0.023 (2.3% of breeding females) in Laos<sup>21</sup>, 0.083 in Myanmar<sup>11</sup>, and varying from 0.095 to 0.155 for a captive elephant population in India over the past century<sup>10</sup>. No systematic assessment explains why the captive birth rates are highly variable and lower than birth rates observed in wild populations. It has been hypothesized to be the consequence of management and husbandry practices, which in effect place breeding secondary to work<sup>4</sup>

The alteration of vital rates in the dynamics of wild-captive Asian elephant populations illustrates the intricacies of mixing natural and human-induced mechanisms<sup>22</sup>. Our objective was to assess the effect of these interactions on the viability of Asian elephant populations using an ecological-economic modelling approach with Lao populations as a case study.

Lao captive elephants are traditionally used as draught animals in the logging industry and share the same habitat as their wild conspecifics. In the last twenty years, the intensification of elephants' workload has affected the fecundity of the population, leading to a sharp decline from 1,000 to 540 individuals<sup>21</sup>. The long gestation and lactation period for female elephants is not compatible with work<sup>4</sup>. Owners consider breeding to represent a loss of income for nearly three years and therefore avoid breeding their elephants, preferring the high incomes generated by the logging industry<sup>4</sup>. Meanwhile, the contiguous Nam Pouy wild population has experienced a sharp decrease from 400–500 individuals in 1990 to 60–80 in 2010<sup>23-25</sup>, further hindering the reproductive potential between wild males and captive females. Finally, the Lao captive population has also been affected by the withdrawal of individuals through the export of captive juveniles to third countries. These exports violate the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES) despite its ratification by the Lao government in 2014. While exports are highly regulated, the trade of captive elephants is unregulated on the domestic market<sup>26</sup>.

To develop a demographic model incorporating economic constraints, we considered the processes impacting the dynamics of captive elephants with a focus on the determinants of fecundity. Therefore, we developed a model including three components: (i) a demographic component structured by sex and age, including the exports of living captive elephants, to estimate the captive population abundance at a given time; (ii) a wild elephant component to account for the wild male breeding potential; (iii) and an economic component to model the elephant owner's decision-making process. Finally, we defined several demographic and socioeconomic scenarios to investigate new ecological-economic equilibria that can serve as a basis to design sustainable management strategy for the species.

## Results

**Changes in breeding practices.** Breeding practices radically changed since Laos transitioned to an open market economy in 1989. Traditionally, captive elephants were used for occasional village work while spending most of their time in small herds, grazing freely in the village surroundings. Mating of captive females with either captive or wild males was uncontrolled.

In 2015, we conducted 52 semi-structured interviews among elephant owners, keepers, and officials throughout the distribution of captive elephants in Laos – 41 out of 50 respondents indicated difficulties finding a captive mate. Among owners of females, the main reason given for the difficulty in finding a captive male was the absence of captive breeding males (21 out of 43) because “there is no male nearby”, “males are kept far away”, “they (males) are engaged in logging and therefore not available for mating”, or “the few males living nearby” are “too young” or “sterile”. Two other reasons were given equally (11): stud fees being too expensive, and their elephant being too tired for breeding. Owners of males indicated having difficulties finding a captive female for mating because “females are not available” (8), “their male is too tired for mating” (8), and “females owners do not want to pay stud fees” (6).

Meanwhile 28 out of 42 respondents considered that mating with a wild male is easier than with a captive one, including all respondents (11) living at the edge of Nam Pouy park and having access to the wild pool. Such a mating system was easier because wild males are “available”, “free of their movements”, “in good health”, “not exhausted unlike the working captive males”, and “do not ask for stud fee payment”.

As a consequence, the reproduction of captive females depends mainly on the presence of wild males in Nam Pouy, and to a lesser extent, on the availability of a captive mate. Unlike in wild populations, the reproduction of captive females is density-dependent towards males.

**How much does an elephant cost?** The captive population dynamics are also driven by the willingness of owners to breed their female. According to economic theory, owners maximize profits, and therefore have to make an economic choice between breeding their female or having them work in the logging industry. The owner estimates how much he can earn by selling a calf after weaning and compares these earnings with the foregone alternative income from logging during the three-year gestation and weaning period.

Our model incorporated the price of an elephant following supply and demand theory – that is, a product (i.e., elephant) is a function of its quantity (i.e., the population). Thus, the price of an elephant can be expressed in terms of the population size given by the demographic component of the model. We associated the sale price given by 26 owners who sold their elephants between 1994 and 2012, to the population size produced by the model at the time of sale. Elephant prices in constant USD were exponentially related to total population (all parameters being different from 0 with p-value < 0.035).

The price of an elephant in constant USD increased with the decrease in population (Fig. 1). From 1990 to 2012, the increasing scarcity of elephants was accompanied by a rapid increase in sale prices exceeding 20,000 constant USD<sup>27</sup>. This result illustrates that the relationship between quantity and price is consistent with economic
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Elephant price is a function of total population abundance. Prices are given in constant USD, base year 2000.</figcaption></figure>

theory. Lower populations (quantity available to the market) lead to higher prices, which lead more owners to conclude it is profitable to breed their elephants rather than using them for logging.

**Population projections depending on ecological-economic scenarios.** Using our ecological-economic model, we investigated the impacts of economic and government policies on future trends in the captive population. The government may impact the population dynamics by revising its policy on the controversial export of elephants to third countries and by fighting illegal trade. The government could also reinforce its wild elephant conservation policy by enforcing the law and effectively protecting the Nam Pouy park wild pool. Finally, the government and private stakeholders could promote pro-birth policies for captive elephants (e.g., by compensating owners for foregone logging returns).

Based on the 2012 population census of captive elephants and on estimates of the Nam Pouy wild population, we selected five scenarios that are representative of the different dynamics of the wild pool and levels of exports. All scenarios are based on a price-dependent domestic market. The reference scenario (i) showed the current level of exports of juveniles and a decreasing wild pool at a rate of 8%. The four other scenarios were based on an effective prevention of exports with (ii) a decreasing wild pool, (iii) a stabilized wild pool to investigate the effect of a conservation policy, (iv) a growing wild pool, and (v) a maternity-leave scheme that financially compensates owners for the loss of income during the breeding period. We ran all scenarios for a period of 100 years with stochasticity in the economic growth rate.

Whatever the scenario, the population decreases during the first 35 years until it recovers from the transient and unstable structure observed in 2012 (Fig. 2a). If current export practices continue (i), the population decline slows, but tends to extinction over the long term with an estimated population of 10 individuals after 480 years. If exports are prevented, the population reaches an asymptotic equilibrium that depends on the dynamics of the wild pool. After 100 years, the captive population reaches an abundance of 355, 371 and 408 individuals for a decreasing (ii), stabilized (iii), or growing (iv) wild pool. Maternity leaves (v) are the only scenario allowing the population to exceed 2012 abundance with 554 individuals. In this case, compensating the owners' loss of income counteracts the prevailing effect of economic incentives on elephant fecundity.

The population decrease due to exports (i) tends to sharply pull prices upward (Fig. 2c), reaching prices over 800,000 constant USD. The population increases for the other scenarios equalize prices from USD 206,000 [s.d. 20,000] (scenario v) to USD 428,000 [s.d. 39,000] (scenario ii). Consequently, the profit ratio (Fig. 2d) stabilizes respectively at 1.36 [s.d. 0.02] (iv), 1.56 [s.d. 0.02] (iii) and 1.65 [s.d. 0.03] (ii). The rising prices in the export scenario (i) brings the profit ratio to 3.35 [s.d. 0.05] after 100 years. As maternity leaves cover the loss of income, profit ratio of scenario (v) tends to infinity.

Fecundity rates (Fig. 2b) follow the same growing trend as the profit ratio, but are constrained by density dependence in reproduction (see Methods). Export scenario (i) displays the higher fecundity rate of 0.117 followed by the maternity leave (v) fecundity rate 0.105. The fecundity rates in the three other scenarios (ii to iv) are approximately 0.095.

**Changes in demand for elephants.** Because we study elephant population dynamics, we built an economic model centered on supply. Variations in supply, the total population, explain price fluctuations given that demand for elephants is proportional to supply. We investigated the effect of a change in demand by applying a 25% decrease and increase in elephant prices. If demand fluctuates over time, using fixed rates allowed us to define a maximum interval bounding the total population estimates under supply and demand equilibrium (Fig. 3).

In the export scenario (i), population estimates range from 141 to 172 for a 25% decrease and a 25% increase in demand for elephants. The declining (ii) and growing (iv) wild pool scenarios have respective population
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Model projections depending on ecological-economic scenarios (a) Population estimates; (b) Fecundity; (c) Elephant prices in constant thousand USD; (d) Profit ratio. The different scenarios within each figure are: (i) Current level of exports and declining wild pool; (ii) No export and declining wild pool at a growth rate of 0.92; (iii) No export and stable wild pool; (iv) No export and growing wild pool at a growth rate of 1.03; (v) Maternity-leave compensation scheme with no export and declining wild pool. The median value of the 500 simulations is shown in contrasted color. The gray line in (d) shows one-year working income as a reference. (d) does not show the profit ratio for the maternity-leave scenario (v) as it tends to infinity.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Population size projections under various scenarios of demand variation. (i) Current level of exports and declining wild pool. (ii) No export and declining wild pool at 8% per year. (iv) No export and growing wild pool at 3% per year. (v) Maternity-leave compensation scheme with no export and declining wild pool. Dot-dashed lines give the total population size for demand increase (upper estimate), dashed lines the total population size for demand decrease (lower estimate) and plain lines the total population size with constant demand (as shown in Fig. 2a). Demand variations do not affect maternity-leave scenario (v).</figcaption></figure>

estimates of 280 to 402 and 330 to 470. Once again, population size in the maternity-leave scenario (vi) is not affected by the variation in demand.

**Elasticity analysis.** We projected the population out 20 years in each scenario and calculated the elasticity of the resulting population growth rate to model parameters (see Supp. Note 3). The parameter showing the highest elasticity (0.033) is the wild pool growth rate (*grw*) for the growing wild pool scenario (iv). It decreases to 0.021 for the stable wild pool scenario (iii) and 0.005 for the declining wild pool (ii), export (i) and maternity-leaves
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Map of captive and wild elephant populations in Lao PDR. Map has been designed with QGIS version 2.14.0-Essen (www.qgis.org) and Photoshop version 9.0 (www.adobe.com). Protected areas shape files are downloaded from The World Database on Protected Areas (WDPA), IUCN and UNEP-WCMC (2017), available at: www.protectedplanet.net.</figcaption></figure>

(v). The female export rate shows a medium-range elasticity (0.012) for the export scenario (i) while it is null for all other scenarios. All other parameters show similar elasticities whatever the scenario with maximum fecundity and sex ratio at a rate of approximately 0.03. Elasticity to female survival rates decline continuously from 0.03 to zero. Elasticity to male survival shows much lower values with a maximum of 0.007 for the 35 years age class. Other model parameters have an elasticity rate lower than 0.001, including all economic parameters.

Discussion

No study has investigated the effects of socio-economic strategies on the fecundity of an animal population whose dynamics are shaped by human economic practices. We demonstrate that fecundity of captive elephants is impacted by: i) the dynamics of the wild elephant pool through mating of captive females by wild males, and ii) the economic incentives of elephant owners to breed their animals. As a result, we expect fecundity to rise again gradually in response to the rapid increase in the price of elephants. The captive population will tend towards an asymptotic limit determined by the dynamics of the wild pool. In any case, the population will tend to extinction if the government does not stop export practices at current rates. Elephant exports, mostly to Chinese zoos and tourism enterprises, threaten the survival of the Laotian population.

Our model illustrates the strong interactions between the wild and captive pools. The value of captive ele-
phants for conservation have been recognized, not only because of their numbers and the threat that their capture
could pose to wild pools, but also for their potential roles in reintroduction, patrolling, human-elephant conflict
mitigation, and scientific research<sup>9,28</sup>. Beyond the risks of extinction of wild populations, we showed that wild
pool dynamics directly impact the survival of the captive population through crossbreeding. The elasticity analy-
sis confirms the influential contribution of a stable or growing wild pool to the growth of the captive population.
The Nam Pouy wild pool was modelled using a geometric growth rate, based on published data<sup>23-25</sup>; however,
these data are scarce and approximate. A better estimate and monitoring of the abundance of this population
would reinforce the recognition of how dynamics of the wild pool affect the demography of the captive pool. For
the time being, it seems urgent to strictly protect the wild populations and consider organizing new wild-captive
breeding areas in Laos. Policy-makers should also consider the side-benefits of a strict conservation strategy for
wild elephants in regards to the cultural and economic role played by captive elephants.

Our ecological-economic approach also emphasizes the impacts of human practices on the dynamics of the captive elephant pool. Several models have been developed to search for potential balance between ecological and economic dynamics. Most of these models favored a macroeconomic approach – searching for a global equi- librium between total profit maximization and population survival<sup>29-31</sup>. In our case, we directly linked fecundity with owner profit maximization. In the absence of reliable macroeconomic data, this microeconomic approach allowed us to model the owner's decision-making process and study its impact on species' demography. However, our elephant price dataset was derived from data for a population between 530 and 900 individuals. Therefore, the estimated price of elephants in a population of less than 200 individuals remains questionable, being far from the original data range. This was one reason we limited our projections to 100 years.

As shown by our model, changes in handling practices have long-term effects on population structure and
dynamics. In 2012, the population was already transient and unstable, with a strong bias in favour of adults over
20 years old (see the age pyramid in Methods - Model calibration). As a result, the model indicates it will take 35
years before the population grows again. Long-lived species, such as the Asian elephant, are highly sensitive to
population inertia³², with long periods before the population may recover from any change in a vital rate. Further investigation of such transient status and its consequences on elephant population dynamics are needed. But our model demonstrates that the commodification of captive elephants and its accompanying changes in handling practices have probably dominated the ecological dynamics of this population, as shown by the short-term elasticity analysis (Supp. Note 3).

Our maternity-leave compensation scheme provides additional insights into the conservation questions surrounding captive and semi-captive populations. If owners are compensated for their loss of income during breeding, the steep economic trade-off between using an elephant for breeding versus work is removed, the elasticity of all economic parameters being null under this scenario. By covering the opportunity cost of breeding, the maternity-leave scheme we analysed allows the population to exceed its 2012 abundance. Such a policy could be funded either by a tax on elephant incomes or directly by the tourism industry, for which calves are highly prized. Because the tourism industry is also a buyer of elephants, supporting such population growth would also put downward pressure on elephant prices (Fig. 2c).

Overall, our ecological-economic model does not explicitly predict the size of the elephant population in a century. It rather provides a framework for policy-makers, conservationists and managers to anticipate the impacts of socio-economic determinants and their policies —such as exports, or wild elephant conservation— on the fate of the population. Moreover, the interactions between captive and wild pools should be further investigated as they set in motion poorly understood demographic, genetic, epidemiological, economic and socio-cultural processes that strongly influence the species' population viability.

## Methods

**Study Species and Area.** Laos has a centuries-long tradition of raising Asian elephants²³. Most captive elephants live in remote areas and are used to extract logs in the forestry industry²⁷,³³. After their daily working hours, they are tied to trees where they can forage by themselves and are moved twice a day by their keepers. Keepers also often release their elephants in the forest for a few months while the keepers are engaged in rice farming³⁴. Most of the country's captive elephants live in Sayaboury province (Fig. 4), which also hosts the Nam Pouy National Park (Lat: 18°12′–18°58′N, Long: 101°04′–101°32′E) – home of the second largest wild population in the country²³,²⁷,³⁴.

**Population data.** Since 2006, the Lao Care and Management Program (LECMP) of the Department of Livestock and Fisheries of Lao PDR has implemented a national campaign to register all captive elephants. Captive elephants are microchipped by the LECMP. Their records, including name, microchip number, sex, age, physical description, medical history, ownership details and location, are saved in a database. Data are updated twice-yearly during LECMP veterinary care visits. We used the data on captive populations structured by sex and age extracted from the database with the year 2012 as a reference²¹. Data on the Nam Pouy wild elephant population from 1990 to 2012 were extracted from the literature²³⁻²⁵. More details on the data used throughout are provided in Supplementary Note 4.

**Trade in captive elephants.** Captive elephants are privately owned in Laos. Domestic trade is allowed after the payment of a tax to local authorities²⁶. Elephant prices are driven by demand on the domestic market. Domestic sales or re-sales do not change the size of the population unlike exports. CITES Convention does not allow international trade in elephants. Despite its ratification by the Lao government in 2004, the export of living elephants continues. Juvenile females illegally cross the border into Thailand and China⁹. National authorities also export young elephants in the form of “rentals” or diplomatic gifts²¹. Being either legal or illegal, exports are considered as permanent, individuals are permanently removed from the captive population.

**Breeding practices in Lao PDR.** We studied mahout (elephant handler) practices and their impacts on wild-captive interactions using ethno-ecological methods. The topics covered in this study included the perception of the species and its wild-captive status, the description of husbandry practices, including working, temporary releasing, breeding and interacting with the wild pool, and the respondent perception of the future. In particular, this study described the free-grazing and husbandry practices in Laos that are potentially interfering with the wild elephant pool, and described the mating systems between captive mates. The breeding practices observed during the ethno-ecological study served as a basis to model the drivers of captive population fecundity. Another article is in preparation that will present the global results of the ethno-ecological study.

In 2015, we conducted 52 semi-structured interviews among elephant owners, handlers, and officials throughout the distribution of captive elephants in Laos. Interviews were conducted in Lao language by the first author, recorded on a Zoom H1 handy recorder, and fully translated into French and English. All methods were carried out in accordance with the guidelines for anthropological surveys³⁵, and in accordance with the legal requirements of France and the National Center for Scientific Research (CNRS) institutional protocols for research. Research protocol has been approved by Montpellier University. All participants provided informed verbal consent before participation. Sampling was done using a stratified method based on the district of residence and ownership status of the respondent. We coded answers using MaxQDA Analytics Pro 12 software for qualitative data analysis.

**Demographic model for the captive population.** We developed a population model composed of a two-sex and 55-age-class matrix to account for the sexual dimorphism in vital rates and sex-specific animal handling practices (Fig. 5).
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Life cycle of captive females.</figcaption></figure>

**Survival.** Demographic data from Laos did not allow estimation of survival rates. We used sex- and age-specific survival rates for a semi-captive elephant population in neighboring Myanmar that shares similar working practices³⁶. Survival is given per year from age 0 to 55 with $sf^n$ the female survival at age n and $sm^n$ the male survival (see data in Supp. Note 4).

**Exports.** From a population perspective, exports to foreign countries are considered permanent migration. We apply an export function to survival for individuals aged 5 to 15 years with $sf^n(1-ef)$ for females and $sm^n(1-em)$ for males. We defined 2 scenarios: a ban on exports (ef and em equal 0) and the reference scenario showing the current level of exports. Since females are preferred to males due to their docility and reproductive potential, the respective rates for ef and em are set to 0.04 and 0.01 for the reference scenario.

**Birth.** Yearly recruitment is given by the product ($sf^0 sr d F$) for females aged from 15 to 50, with $sf^0$ the female survival before their first year, sr a balanced sex ratio at birth, d the density dependence function and F the fecundity function, described hereafter. The same formula is used for males with $sm^0$ as the male survival rate.

**Density-dependent reproduction.** As captive males are engaged in logging activities, many owners of females do not find male mates. Working males are chained to a tree during their resting time and are therefore not available for mating. A prior breeding agreement between female and male owners is therefore required. But female owners found the costs of stud fees too expensive and male owners are reluctant to engage their males in mating because mating is time-consuming and believed to reduce a male's working strength. Most births therefore result from mating between a captive female and a wild male from Nam Pouy Park. When the breeding process depends on external factors, such as the availability of captive males or the presence of a wild mating pool, the fecundity function is weighted by a coefficient $d^{37}$. We introduced a density-dependent reproductive function in relation to the proportion of the mating pool³⁸ where the probability for a captive female to reproduce is dependent on the ratio of the number of breeding males from the captive and wild pools to the number of captive breeding females:

$$d = \frac{1}{1 + e^{-b \frac{N_m}{N_f}}} \quad (1)$$

with $N_f$ the number of captive females of breeding age (15 to 50 years), $N_m$ the number of captive breeding males from 30 to 55 years plus the estimated number of wild males living in Nam Pouy Park, and $b$ the stretching exponent of the function. The parameter $b$ allows the growth of the function to be adjusted. In the case of this long-lived polygynous species, we set the stretching exponent to 2. Thus, for a given population of 200 reproductive females and 20, 60 or 100 breeding males respectively, fecundity will be affected by a coefficient $d$ equal to 0.57, 0.71 or 0.92, illustrating a steep decrease in the reproductive potential due to the unavailability of males (Fig. 6).

**Population dynamics of the wild pool.** The reproduction of captive females depends partially on the Nam Pouy wild pool, as explained above. Given the limited information available on this population, we model the change in Nam Pouy wild population using a population growth rate $wm(t+1) = grw \times gm(t)$, with $gm(t)$ the total population at time t ($gm(0) = 68$ the population at time 0) and $grw$ the yearly average growth rate of Nam Pouy population. We select three different growth rates to explore how wild elephant conservation policies could impact the captive dynamics: an increasing wild pool at a rate of 3% per year, a stable wild pool, and a declining wild pool at a rate of 8% per year. This annual decrease of 8% reflects the sharp decline of this population over the last 20 years (see data in Supplementary Note 4).

**Economic decision-making.** According to microeconomic theory, owners tend to maximize profit. Because breeding is not compatible with female elephants working, elephant owners make an economic choice between having their female elephants breed or work. Thus, our economic model assumes owners maximize expected profits by choosing whether or not to breed their females, where the expectation arises because owners do not know future elephant prices or logging income with certainty. This decision framework implies a decision rule where owners will choose to breed up until the expected net present value of breeding exactly equals the
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Density-dependent reproductive coefficient depending on the number of females. The density-dependent reproductive function is plotted for several male breeding stock levels: from 20 males (blue line) up to 400 (red line), at 20-male intervals. For a given male stock (one line in the chart), the density-dependent coefficient decreases when the number of females increases, meaning that fecundity is negatively impacted by the excess of females relative to breeding males.</figcaption></figure>

opportunity cost of breeding (i.e., foregone timber income). This is equivalent to choosing to breed until the profit ratio (2) is equal to or greater than 1.

At any point in time, the owner compares how much they could receive in net present value by selling a calf after weaning (*Price*<sub>t+3</sub>), knowing that the survival at 3 years (*surv*<sub>3</sub>) is 0.75, with the alternative income from logging during the 3 years. The profit ratio *P* is therefore:

$$P = \left[ \text{surv}_3 * \frac{\text{Price}_{t+3}}{(1+r)^4} \right] \left[ \frac{\text{Inc}_t}{1+r} + \frac{\text{Inc}_{t+1}}{(1+r)^2} + \frac{\text{Inc}_{t+2}}{(1+r)^3} \right] \quad (2)$$

with *r* = 10%, the discount rate that takes into account the time value of money or the risk aversion for agricultural projects in Asian developing countries<sup>39–41</sup>. Thus, a profit ratio greater than 1 implies there is an economic incentive to breed females, whereas a profit ratio less than 1 implies that the opportunity cost of foregone logging profits is too large to justify breeding. More information on owners’ price estimation methods are given in Supplementary Note 2. Finally, the costs related to each activity do not appear in the calculation of the profit because they are limited – elephants are foraging freely in Laos, and the cost of guarding is identical with or without a calf. As a result, feeding and labor costs of the two alternatives do not influence the marginal decision to breed or not.

All monetary values, including incomes and profit, have been converted to constant US dollars (USD) in the year 2000 using the methods of Williamson<sup>42</sup> in order to correct prices for the high inflation in Laos.

**Price and Income.** According to supply and demand theory, the price of an elephant can be expressed in terms of the population size given by the demographic components of the model. From the ethno-ecological study conducted among 26 owners who sold elephants between 1994 and 2012, we assigned each reported price to the population estimate given by the model at the time of sale. Elephant price in constant USD follows an exponential relationship:

$$Price = v_1 e^{v_2 \text{ population}} \quad (3)$$

with $v_1$ (Est. = 119200, $t$ = 2.234, $p$ = 0.035) and $v_2$ (Est. = -0.003633, $t$ = -4.697, $p$ ≪ 0.01).

Working incomes are estimated from IC Suter’s 2011 survey<sup>27</sup> among 133 elephant owners and keepers, which indicated that average income was equivalent to 4,070 constant USD (i.e., 6.11 times the constant per capita Gross Domestic Product (GDP)). Economists generally use GDP per capita as an indicator of changes in household incomes<sup>43</sup>.

To estimate prices and incomes after 2012, we applied an average economic growth rate (*eco*<sub>t</sub>) to the owner income and price of an elephant according to the long-term GDP projections from the Asian Development Bank<sup>44</sup>. Moreover, we added stochastic variation to the projected economic growth rate following a normal distribution with mean 1 and a variance of 0.01. Thus, from 2012 and forward, prices at time *t* equal: *Future prices*<sub>*t*</sub> = *Reg(price*<sub>*t*</sub>) * ∏*<sub>*j*=1</sub><sup>*t*</sup> *eco*<sub>*j*</sub> * N(1, 0.01)*, with *Reg(price)* the price function (Equation 3). Income at time *t* + 1 is given by *Income*<sub>*t*+1</sub> = *Income*<sub>*t*</sub> * *eco*<sub>*t*</sub> * N(1, 0.01)*.

**Linking fecundity with elephant owners’ profit.** Fecundity is driven by the financial incentive of owners to breed their elephants instead of having them work. When working income exceeds that of breeding (profit ratio *P* < 1), the incentive to breed is low and fecundity is minimal or opportunistic (*f*<sub>min</sub>). Fecundity increases as
<figure><img src="image_8.png" /><figcaption>**Figure 7.** Fecundity as a function of the owner's profit. Profit is the ratio of the calf's selling price to the income generated by the foregone work. Fecundity is maximum when breeding generates twice as much as working.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 8.** Age pyramid of captive elephants in 2012. Dotted lines show the number of males and females given by the 2012 census. Solid lines indicate the results of the model projection for year 2012. The green line shows the structure of the population in absence of export, confirming the effect of this practice on the number of females aged from 20 to 30 years.</figcaption></figure>

profit increases up to a maximum ($f_{max}$) set by the biological constraints of the species (Fig. 7). Thus, we model fecundity (F) as a logistic function:

$$F = f_{max} \frac{1}{1 + \left( \frac{f_{max}}{f_{min}} - 1 \right) e^{-p\sigma}} \quad (4)$$

where P is the profit ratio (2). The function is bounded respectively by $f_{max}$ the maximal fecundity and $f_{min}$ the opportunistic fecundity. The parameter $\sigma$ is used to fit the steepness of the function⁴⁵ so that fecundity is maximized when the profit ratio equals 2. We retained a profit ratio of 2 for $f_{max}$ to account for the most extreme case in breeding arrangements (i.e., where all owners consider it profitable to breed their females). Such a rate is a strong incentive for any owner able to pair his female with a wild male. However, the owner of a captive sire usually claims a maximum of one third of the calf value or calf ownership as stud fee, leaving 2/3 of the calf price to the female's owner (Netvalue = 2/3). In this case the profit for the female's owner should be positive if: Price * surv₃ * Netvalue ≥ $\sum_{t=1}^{3} Inc_t$ (derived from Equation 2), or profit ratio $P \ge 2$. As a result, parameter values used in our model are $f_{max} = 0.155^{10}$, $f_{min} = 0.02^{21}$, and $\sigma = 3$.
**Sensitivity and elasticity analyses.** The sensitivity analysis allows studying the effect of variation in the model parameters on the population dynamics as measured by their asymptotic growth rate (λ). The asymptotic population growth rate, calculated as the dominant eigenvalue of the matrix, is calculated using the R popbio package⁴⁶. We used elasticities instead of sensitivities because these metrics take into account the unit in which the demographic parameters are measured and thus make comparisons possible between these parameters. The elasticity (EL) of parameter α is given by:

$$EL(\alpha) = \frac{(\lambda' - \lambda)}{\lambda(1 + \Delta\alpha)}$$

where λ' is the new value of λ after applying a variation Δ to the parameter α.

**Model calibration.** To calibrate our model, we projected the population over 20 years starting from a stable structured population of one thousand individuals in 1992 and compared its demographic structure with the observed structure in 2012. Population abundance reached 530 individuals, consistent with 2012 estimate of 527 individuals. Moreover, the population structure is highly transient and unstable. A strong bias is observed in favor of elephants over 20 years old and a deficit of females in the 20–29 years age class resulting from our hypotheses (sharp drop in fecundity and export of juvenile females during the last twenty years). The population structure derived in the model (Fig. 8) cannot be distinguished from the structure inferred from the 2012 census data (Chisq = 13.6, df = 11, p = 0.25). All analyses were performed with R version 3.2.3⁴⁷.

**Data Availability.** All data generated or analysed during this study are included in this published article (and its Supplementary Information files).

## References

1. Hedges, S., Tyson, M. J., Sitompul, A. F., Kinnaird, M. F. & Gunaryadi, D. Distribution, status, and conservation needs of Asian elephants (Elephas maximus) in Lampung Province, Sumatra, Indonesia. *Biol. Conserv.* **124**, 35–48 (2005).

2. Sukumar, R. A brief review of the status, distribution and biology of wild Asian elephants Elephas maximus. *Int. Zoo Yearb.* **40**, 1–8 (2006).

3. Choudhury, A. et al. Elephas maximus. *Iucn redlist* Available at: http://www.iucnredlist.org/details/7140/0 (2008).

4. Leimgruber, P. et al. Modeling population viability of captive elephants in Myanmar (Burma): implications for wild populations. *Anim. Conserv.* **11**, 198–205 (2008).

5. Rees, P. A. The welfare and conservation of Asian elephants – a reply to Sukumar. *Oryx* **37**, 20–22 (2003).

6. Csuti, B. In *Biology, Medicine, and Surgery of Elephants* (eds Fowler, M. E. & Mikota, S. K.) 15–22 (Blackwell publishing, 2006).

7. Pimmanrojnagool, V. & Wanghongsa, S. A study of street wandering elephants in Bangkok and the socio-economic life of their mahouts. In *Giants on our hands: Proceedings of the international workshop on the domesticated Asian elephant* (eds. Baker, I. & Kashio, M.) 34–42 (Food and Agriculture Organization of the United Nations, Regional Office for Asia and the Pacific, 2002).

8. Sukumar, R. *The living elephants: evolutionary ecology, behavior, and conservation*. (Oxford University Press, 2003).

9. Lair, R. *Gone astray: the care and management of the Asian elephant in domesticity*. (Food and Agriculture Organization of the United Nations, 1997).

10. Sukumar, R., Krishnamurthy, V., Wemmer, C. & Rodden, M. Demography of Captive Asian Elephants (Elephas maximus) in Southern India. *Zoo Biol.* **27**, 263–272 (1997).

11. Mar, K. U. The demography and life history strategies of timber elephants in Myanmar. (University College London, 2007).

12. Thitaram, C. et al. Genetic assessment of captive elephant (Elephas maximus) populations in Thailand. *Conserv. Genet.* **11**, 325–330 (2010).

13. Godfrey, A. & Kongmuang, C. Distribution, Demography and Basic Husbandry of the Asian Elephant in the Tourism Industry in Northern Thailand. *Gajah* **30**, 13–18 (2009).

14. Nijman, V. *An assessment of the live elephant trade in Thailand*. (Traffic International, 2014).

15. Hayward, A. D., Mar, K. U., Lahdenperä, M. & Lummaa, V. Early reproductive investment, senescence and lifetime reproductive success in female Asian elephants. *J. Evol. Biol.* **27**, 772–83 (2014).

16. Robinson, M. R., Mar, K. U. & Lummaa, V. Senescence and age-specific trade-offs between reproduction and survival in female Asian elephants. *Ecol. Lett.* **15**, 260–6 (2012).

17. Mumby, H. S., Courtiol, A., Mar, K. U. & Lummaa, V. Birth seasonality and calf mortality in a large population of Asian elephants. *Ecol. Evol.* **3**, 3794–803 (2013).

18. Lahdenperä, M., Mar, K. U. & Lummaa, V. Reproductive cessation and post-reproductive lifespan in Asian elephants and pre-industrial humans. *Front. Zool.* **11**, 54 (2014).

19. Mumby, H. S., Courtiol, A., Mar, K. U. & Lummaa, V. Climatic variation and age-specific survival in Asian elephants from Myanmar. *Ecology* **94**, 1131–41 (2013).

20. Mumby, H. S. et al. Stress and body condition are associated with climate and demography in Asian elephants. *Conserv. ation Physiol.* **3**, 1–14 (2015).

21. Suter, I. C., Maurer, G. & Baxter, G. Population viability of captive Asian elephants in the Lao PDR. *Endanger. Species Res.* **24**, 1–7 (2014).

22. Lebreton, J. Dynamical and statistical models for exploited populations. *Aust. N. Z. J. Stat.* **47**, 49–63 (2005).

23. Khounboline, K. Current Status of Asian Elephants in Lao PDR. *Gajah* **35**, 62–66 (2011).

24. Santiapillai, C. Elephants under threat in Laos. *Gajah* **10**, 10 (1993).

25. Duckworth, J., Salter, R. & Kounboline, K. Wildlife in Lao PDR: 1999 status report. (IUCN Lao PDR, 1999).

26. Norachack, B. In *Giants at our hands: Proceedings of the International workshop on the domesticated Asian Elephant* (eds. Baker, I. & Kashio, M.) 173–180 (FAO Regional Office for Asia and the Pacific, 2002).

27. Suter, I. C., Hockings, M. & Baxter, G. S. Changes in Elephant Ownership and Employment in the Lao PDR: Implications for the Elephant-Based Logging and Tourism Industries. *Hum. Dimens. Wildl.* **18**, 279–291 (2013).

28. Dublin, H. et al. *Elephant Range States Meeting Report*. (Kuala Lumpur, 2006).

29. Horan, R. D. & Bulte, E. H. Resource or nuisance? Managing African elephants as a multi-use species. in *American Agricultural Economics Association Conference* (2001).

30. Rakotoarison, H. A cost benefit analysis of big game: the case of wild boar in the Aquitaine region. In *Colloque SFER Chasse, Territoires et Développement durable: outils d'analyse, enjeux et perspectives* 1–27 (2009).
31. Timmermann, K. et al. Development and Use of a Bioeconomic Model for Management of Mussel Fisheries under Different Nutrient Regimes in the Temperate Estuary of the Limfjord, Denmark. *Ecol. Soc.* **19** (2014).

32. Koons, D. N., Holmes, R. R. & Grand, J. B. Population inertia and its sensitivity to changes in vital rates and population structure. *Ecology* **88**, 2857–2867 (2007).

33. Lair, R. In *Giant on our hands: Proceedings of the International Workshop on the Domesticated Asian Elephant* (eds Baker, I. & Kashio, M.) (FAO Regional Office for Asia and the Pacific, 2002).

34. Lassausaie, J. et al. Tuberculosis in Laos, who is at risk: the mahouts or their elephants? *Epidemiol. Infect.* **143**, 922–931 (2015).

35. Baude, O. et al. Corpus Oraux. Guide des bonnes pratiques 2006. (Presses Universitaires, 2006).

36. Mar, K. U. The studbook of timber elephants of Myanmar with special reference to survivorship analysis. in *Giants on our hands: Proceedings of the international workshop on the domesticated Asian elephant* 195–211 (FAO Regional Office for Asia and the Pacific, 2002).

37. Legendre, S. et al. Demographic Stochasticity and Social Mating System in the Process of Extinction of Small Populations: The Case of Passerines Introduced to New Zealand. *Am. Nat.* **153**, 449–463 (1999).

38. Sundelöf, A. & Åberg, P. Birth functions in stage structured two-sex models. *Ecol. Modell.* **193**, 787–795 (2006).

39. Fisher, B., Edwards, D. P., Giam, X. & Wilcove, D. S. The high costs of conserving Southeast Asia's lowland rainforests. *Front. Ecol. Environ.* **9**, 329–334 (2011).

40. Tomish, T., Noordwijk, M. van & Budidarsono, S. in *Tradeoffs or Synergies?* (eds Lee, D. & Barret, C. B.) 221–245 (Centre for Agricultural Bioscience International, 2001).

41. Markandya, A. & Pearce, D. W. Development, the Environment, and the Social Rate of Discount. *World Bank Res. Obs.* **6**, 137–152 (1991).

42. Williamson, S. H. & Officer, L. H. Seven Ways to Compute the Relative Value of a U.S. Dollar Amount, 1774 to present. *Measuring Worth* Available at: www.measuringworth.com/worthmeasures.php (2010).

43. Boarini, P. R. & Johansson, Å. Les indicateurs alternatifs du bien-être. *Cah. Stat. l'OCDE* (2006).

44. Asian Development Bank. *Long term projections of Asian GDP and trade*. (ADB, 2011).

45. Verhulst, P.-F. Recherches mathématiques sur la loi d'accroissement de la population. *Nouv. Mémoires l'Académie R. des Sci. B.-lett. Bruxelles* 1–42 (1845).

46. Stubben, C. & Milligan, B. Estimating and Analyzing Demographic Models. *J. Stat. Softw.* **22** (2007).

47. R Development Core Team. R: A Language and Environment for Statistical Computing. (2011).

## Acknowledgements

We thank Ingrid C. Suter for sharing data on the mahout survey, Doyle McKey and Sarah Cubaynes for reviewing the manuscript. This work was partly supported by a grant from the French National Research Agency, reference ANR-16-CE02-0007.

## Author Contributions

G.M. and O.G. designed the study. V.C. and G.M. collected the data. G.M. built the model and analysed the data with O.G. and B.R. G.M. wrote the paper with contributions from O.G. and B.R. All authors reviewed the manuscript.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-017-13907-x.

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017