BACKGROUND
the asian elephant  has been captive in thailand for approximately 4000 years, and is classified as a working animal under the beast of burden act of  <dig> .
although traditionally used in war and religious ceremonies, in the 18th and 19th centuries, the asian elephant was used as a source of transportation for logs and people, primarily in the logging industry .
however, over the past century there has been a substantial decline in the population of working elephants, from approximately  <dig>  <dig> to fewer than  <dig>  <dig> elephants today  <cit> .
this decline was due to an increase in the use of machinery instead of elephants, in conjunction with the logging ban of 1989; thus, the need to capture wild elephants as replacement animals was significantly reduced .
the remaining elephants are now primarily found working in the tourism industry, with  <dig> camps identified across thailand in  <dig> by the national elephant institute .
although tourism offers a limited solution, many researchers are in agreement that the overall welfare of these elephants depends greatly on their economic value, utility, and the financial position of their owners; all of which are considerably enhanced as a result of the tourism industry .
concerns over poor welfare due to factors such as artificially created social groups, inadequate quantity and quality of food, and poor veterinary care have been investigated to some extent .
however, relatively few studies have investigated the physical effects of work on this population of working elephants.
there is only a brief mention in one report of the presence of abrasions caused by the movement of the tourism saddle and girth strap  <cit> .
this is in contrast to the domestic equid industry involving donkeys, horses and mules, where a number of studies have implicated overworking/overloading and improper saddle/harness design as major causal factors for external lesions .
similar welfare concerns are likely relevant to working elephants.
consequently, the objectives of this study were to determine the nature of equipment-related skin injuries in a population of asian elephants working in the tourism industry by determining the following: 1) the prevalence of external injuries in four saddle related areas: neck, girth, back and tail, among working tourism elephants in thailand; and 2) potential risk factors for skin lesions related to working conditions, including demographic factors and saddle design.
methods
data for this cross-sectional study were collected between may and july  <dig>  by the lead researcher  who worked in partnership with the elephant hospital of the national elephant institute, formerly known as the thai elephant conservation centre.
data were collected from  <dig> camps across thailand that received routine visitation from veterinarians from the livestock and wildlife hospital of mahidol university, or the national elephant institute’s elephant hospital.
elephants were eligible for the study if they worked in the tourism industry and wore a saddle .
each elephant underwent a veterinary examination and elephant keepers  answered a series of standardized questions pertaining to each elephant and its saddle-wearing history using a translator .
the outcome of interest for the present study was the presence of an ‘active lesion’ in any one of four anatomical areas: neck, girth, back and tail.
lesions were placed initially into  <dig> categories .
however, to achieve adequate power for statistical analyses, lesions were classified as either ‘active’ or ‘inactive/absent’.
active lesions were those that could be described as rubbed/pink, raw, a full depth ulcer, an abscess , a healing lesion, or depigmented skin.
a body region was categorized as having an inactive/absent lesion if there was no lesion, a healed lesion/scar, or a callus.
the following independent variables were investigated to determine their association with the presence of an active lesion: age of the elephant ; region of the body ; number of hours worked per day ; the provision of a break ; weight of the saddle ; and contact material .
to allow for a greater number of observations per category, shredded bark and bark were combined into a single “bark” category.
descriptive analyses were carried out using microsoft office excel,  <dig>  and stata mp <dig>  .
the prevalence proportions of elephants with at least one lesion , with at least one active lesion, and with no lesions were estimated.
the prevalence of active lesions was also summarized by body region.
among camps, the prevalence of elephants with at least one active lesion was estimated.
explanatory variables of interest were summarized and tabulated by camp.
using stata mp <dig> , univariable and multivariable multilevel logistic regression models were built to examine the association between the presence of an active lesion at any of the four anatomical sites  and the explanatory variables: age, working hours, contact material, break, and saddle weight.
these models contained random intercepts for camp and elephant to account for clustering, as multiple elephants were observed per camp, and multiple body regions were examined per elephant.
the models were used to determine variance estimates at both the elephant and camp level.
decisions to include random intercepts were made based on model fit using akaike’s information criterion  score, in combination with a likelihood-ratio test .
correlation among independent variables was tested using spearman’s rank correlations.
any two variables found to be highly correlated at r ≥  <dig>  were assessed for completeness of data and biological relevance.
this information was used to determine the most reliable variable to be included in the model, to avoid issues associated with collinearity.
the linearity assumption was examined for continuous explanatory variables  using locally weighted regression with the lowess command.
if the lowess smoother indicated a curvilinear relationship, a quadratic term was investigated and remained in the model if α ≤  <dig> .
variables with p ≤  <dig>  on univariable analysis  were considered for inclusion in the multivariable model.
variables included in the final multilevel multivariable model showed either statistical significance , or evidence of acting as a confounding variable .
interaction effects were initially examined, but due to estimation issues associated with relatively small sample size, we limited our analysis to a main effects model.
pearson residuals and best linear unbiased predictors  were examined to assess the fit of the multilevel multivariable model.
RESULTS
descriptive statistics
in total,  <dig> body areas, from  <dig> elephants and  <dig> camps were examined for saddle related lesions.
the prevalences of elephants  with at least one active lesion, at least one lesion , and no lesions were found to be  <dig>  % ,  <dig>  % , and  <dig>  % , respectively.
the distribution of active lesions by body area was as follows:  <dig>  % on the neck ,  <dig>  % on the girth ,  <dig>  % on the back , and  <dig>  % on the tail .
among camps, the average prevalence of elephants with at least one active lesion was  <dig>  % , although this value ranged from  <dig> to 100 % .
the median value of prevalence of elephants with at least one active lesion by camp was  <dig>  %.
the number of elephants per camp ranged from  <dig> to  <dig> .
the age of the elephants in years ranged from  <dig> to  <dig>  with an average age of  <dig>  years.
scheduled working hours  were found to be relatively consistent among elephants within a camp, although this value ranged between  <dig> to 10 hours among camps .
whether or not the elephants were provided a break was also found to be a camp level variable .
saddle weight was found to be equivalent for all elephants among the same camp, with 50 % of the camps using heavy  saddles, and 50 % using light  saddles .
the contact material of the saddle pad was found to vary both among camps, and between elephants within the same camp .table  <dig> prevalence of elephants with at least one lesion , compared to prevalence of elephants with at least one active lesion by tourism camp in thailand
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
*an active lesion is defined as a lesion categorized as rubbed/pink, raw, full depth ulcer, abscess , healing lesion, depigmented skin, or other.
an inactive lesion is defined as a lesion categorized as either a healed lesion, or a callus
^exact 95 % confidence interval
#one-sided, exact  <dig>  % confidence interval
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
nd = no data
*break is defined as a period of time when either the entire saddle was completely removed, or the girth strap was loosened
**saddle weight is dichotomized into light  and heavy
***the individual elephant that works 2 hours/day gets the break
^each elephant had only one contact material, but the same material was not necessarily used for all animals within a camp
risk factor analysis
initially, the multilevel model included random intercepts for both elephant and camp.
however, the random intercept for elephant was removed because the variance component was negligible  and model fit, based on aic and a likelihood-ratio test, was not improved by its inclusion.
it should be noted that the neck region variable was removed from both univariable and multivariable analyses.
there were so few lesions found in this region that it resulted in estimation problems for coefficients in our multilevel multivariable model.
i) univariable analysis
based on univariable analysis, active lesions were more likely to be found in the back region in comparison to the other anatomical areas .
increasing age was also identified by univariable analysis as being a risk factor for having an active lesion .
there appeared to be a significant quadratic relationship between working hours/day and the log odds of having an active lesion .
there was no significant difference in the prevalence of active lesions associated with the use of a heavy  or light  saddle, or with whether or not the elephants were provided a break .
we observed no significant differences between the use of the three contact materials  compared to bark.table  <dig> univariable multilevel* logistic regression models of factors associated with the presence of an active lesion on elephants from tourism camps in thailand
ref = referent category
n = number of elephants with data collected for the indicated variable
m = number of camps with data collected for the indicated variable
*models include a random intercept for camp
**break is defined as a period of time when either the girth strap was loosened, or the whole saddle was completely removed
ii) multivariable analysis
in the final multilevel multivariable model, the following variables were included: body region, contact material, saddle weight, the provision of a break, age, and working hours/day and working hours/day squared.
as observed in the univariable analysis, the back region was found to be at higher odds for having an active lesion compared to the girth and tail regions .
in comparison to each of the bark, carpet, and blanket materials, the use of a rice sack was found to be a significant risk factor for having an active lesion .
further comparisons between the bark, carpet, and blanket materials revealed no significant difference in the odds of having an active lesion .
increasing age was found to be associated with increased odds of having an active lesion.
a significant quadratic term was again observed for working hours/day.
initially the odds of having an active lesion increased with working hours/day; however this relationship peaked at 6– <dig> hours, followed by a decline .
after controlling for the other variables in the model, receiving a break significantly increased the odds of having an active lesion .
no significant difference was found between heavy and light saddles in association with the odds of having an active lesion.
however, this variable was found to act as a confounder for both contact material and age, and was retained in the final model.table  <dig> multivariable multilevel* logistic regression model of factors associated with the presence of an active lesion on  <dig> sites on elephants working in tourism camps in thailand
body area
contact material
saddle weight
break**
age
working hours/day
variance =  <dig>
*model includes a random intercept for camp
**break is defined as a period of time when either the girth strap was loosened, or the whole saddle was completely removed
note: odds ratios are calculated using the second listed variable as the referent category
*model used for comparisons includes a random intercept for camp
fig.
<dig> relationship between the predicted log odds of an active lesion and hours worked per day.
predicted log odds of active lesion versus number of hours worked per day based on final multilevel multivariable logistic regression model for factors associated with active lesions on elephants from tourism camps in thailand
an intercept only analysis revealed that 20 % of the model variance was explained at the camp level.
in the final multilevel multivariable analysis, we found 13 % of the variance in the outcome was accounted for at the camp level.
we observed no outliers based on examination of the pearson residuals.
based on visual examination of the blups, we found that model assumptions concerning normality and homogeneity of variance were met.
discussion
the prevalence and associated risk factors for skin lesions were examined using data from  <dig> working elephants from  <dig> thailand tourism camps.
descriptive analyses revealed a high prevalence of active lesions among these elephants especially in the back region.
risk factors associated with increased odds of having an active lesion were identified as: body region , the use of rice sacks as padding material in contact with the skin, increasing age of the elephant, the provision of a break, and longer working days.
a high prevalence of active lesions associated with working conditions was found, as  <dig>  % of these elephants had at least one active lesion, and  <dig>  % had at least one lesion .
this finding is comparable to the prevalence of pack wounds found in various studies on working equids .
interestingly, virtually none of the model variance was explained at the elephant level, suggesting lesions among sites on an elephant occur relatively independently of each other.
approximately 13 % of the model variance was explained at the camp level, even after accounting for several camp level variables.
these variables included working hours/day, the provision of a break, and saddle weight.
this suggests that camp level management has a moderate impact on elephant welfare associated with external injuries, potentially indicating an area for more detailed investigation in future studies.
the back region appears to be a high risk area for the presence of an active lesion relative to the other anatomical sites, followed by the girth region.
these findings correspond to the results of various studies regarding working donkeys and mules.
sells et al.
identified the withers and shoulder as the regions of highest lesion prevalence, at 40 % and 31 %, respectively  <cit> .
biffa and woldemeskel saw the highest lesion prevalence on the withers at  <dig>  %  <cit> .
pearson et al.
documented the back region to be the area of highest lesion prevalence at 90-93 %, and the neck region as a lower risk area at 4-10 %  <cit> .
pritchard et al.
saw lesions on the head, neck, ribs, flank and tail base in less than 10 % of equids studied, and observed the highest prevalence of lesions in mules on the breast/shoulder, withers and girth regions   <cit> .
the results of our study, along with these comparable equid results, suggest a primary area for intervention.
future studies and current management/veterinary care regarding working elephants should focus efforts primarily on the back region and secondarily on the girth region, as opposed to the neck and tail, where lesions were relatively rare.
padding material in contact with the skin was a significant factor for the prevalence of active lesions.
in particular, rice sacks significantly increased the odds of having an active lesion relative to all other padding materials.
avoiding their use in favour of bark, carpet, or blankets would be beneficial in terms of decreasing the odds of having an active lesion.
bark, carpet and blankets were not found to be significantly different from one another; therefore, at this time there are no recommendations for a potential standard material.
increasing age had an expected association with increasing odds of having an active lesion.
although an odds ratio of  <dig>  is relatively small, if we take into account that these are long lived animals, the odds of an active lesion increases  <dig>  times for every 15 year increase in age.
due to the nature of a cross-sectional study, there is an issue concerning prevalence vs. incidence of active lesions on older animals.
the study design did not allow us to discern whether older animals were more likely to have active lesions due to increased susceptibility, decreased healing times, or some combination of both.
regardless, at this time, we recommend that a higher level of monitoring/attention be given to older elephants by their owners or caretakers.
it is important to note that due to the inclusion of the age variable, which had the largest amount of missing data , the number of elephants included in the final model decreased from  <dig> to  <dig>  however, the direction  of the odds ratios of the other variables remained the same; only the magnitude of these effects was altered.
as we thought age was an important biological variable that showed statistical significance, we decided to include it in the final model despite the loss of some observations for the multivariable analysis.
surprisingly, the provision of a break was found to be a risk factor associated with increased odds of having an active lesion.
it is possible that camps that offer more breaks may be doing so in response to a current or previous problem with lesions.
but again, as a result of the limitations of a cross-sectional study, the incidence of active lesions as a result of breaks could not be investigated directly.
control trials or cohort studies should be conducted to investigate the effect of breaks on the incidence of lesions along with more specific descriptions concerning their timing  and the conditions associated with these breaks.
for instance, the loosening of the girth strap during a break may cause more movement of the saddle equipment and, therefore, more abrasions.
the association between the odds of having an active lesion and working hours/day had a quadratic relationship.
the predicted log odds of an active lesion were observed to initially increase with increasing working hours/day, followed by a peak at 6– <dig> hours and subsequent decline.
this drop was an unexpected finding, but may be related to the “healthy worker bias”.
for instance, more resilient elephants may be selected more frequently to work more hours, or those elephants that sustain too many lesions may have decreased numbers of working hours, or be removed from the camp.
this could result in a selection bias and an impression of a healthier population of elephants working longer hours, which may explain the quadratic relationship identified.
sells et al.
found a similar result in their study analysing working equids in morocco, where increased frequency of work/week was found to be associated with a decreased risk of sustaining a pack wound  <cit> .
they also identified a healthy worker selection bias as a possible explanation for these findings  <cit> .
due to the nature of the present study, the prevalence of skin lesions was observed only in animals currently working in the camp, and data regarding total working hours/day were gathered only for current management.
therefore, to control for this potential bias, future studies analysing the role of working hours should examine all elephants that have been used by the camp, and take information on any changes in working hours among the elephants.
it may be prudent at this time to recommend limiting work to <6 hours/day since the odds of observing an active lesion was highest at this time point.
although the weight of the saddle was found to be a confounding variable for the present study, a significant difference between heavy  and light  saddles was not found to be associated with an increased odds of having an active lesion.
however, multiple studies investigating similar factors among working equids have implicated overloading as a major causal factor for skin lesions  <cit> .
kontogeorgopoulos also cautions that although these elephants can carry up to 300 kg, it is recommended they carry no more than 200 kg  <cit> .
therefore, the total weight that elephants carry  should be taken into consideration.
it is possible that relative to the individuals being carried, the weight of the saddle is not significant.
it is also possible that the lack of variability in the weight variable within camps may have impacted our ability to find an effect of weight on lesion outcome.
it is also important to note that by dichotomizing the weight of the saddles due to limitations in their assessment in the field, a great deal of information was lost that could have been identified if modelled as a continuous variable.
CONCLUSIONS
in conclusion, lesions related to working conditions were found to be highly prevalent among the tourism elephants studied.
these lesions were found to be associated with risk factors including the use of rice sacks as padding material in contact with the skin, longer working days , and the provision of a break.
the back region was identified as being at an increased risk of having an active lesion, followed by the girth region.
while the cross-sectional nature of the study posed several inferential limitations, at this time we recommend the following until longitudinal studies can be conducted: increased monitoring of older elephants, increased monitoring of the back region of elephants of all ages, working hours limited to less than 6 hours per day, and the avoidance of rice sacks as padding material in contact with the skin.
additional file
additional file 1:
elephant observation sheet.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
sm, ta, and das planned and designed the study.
sm led the collection and recording of all data.
os and dlp conducted the statistical analyses and wrote the initial draft manuscript.
all authors were involved in writing and editing the final manuscript.
all authors read and approved the final manuscript.

 the asian elephant has been captive in thailand for approximately 4000 years, and is classified as a working animal under the beast of burden act of. although traditionally used in war and religious ceremonies, in the 18th and 19th centuries, the asian elephant was used as a source of transportation for logs and people, primarily in the logging industry. however, over the past century there has been a substantial decline in the population of working elephants, from approximately to fewer than elephants today. this decline was due to an increase in the use of machinery instead of elephants, in conjunction with the logging ban of 1989; thus, the need to capture wild elephants as replacement animals was significantly reduced. the remaining elephants are now primarily found working in the tourism industry, with camps identified across thailand in by the national elephant institute. although tourism offers a limited solution, many researchers are in agreement that the overall welfare of these elephants depends greatly on their economic value, utility, and the financial position of their owners; all of which are considerably enhanced as a result of the tourism industry. concerns over poor welfare due to factors such as artificially created social groups, inadequate quantity and quality of food, and poor veterinary care have been investigated to some extent. however, relatively few studies have investigated the physical effects of work on this population of working elephants. there is only a brief mention in one report of the presence of abrasions caused by the movement of the tourism saddle and girth strap. this is in contrast to the domestic equid industry involving donkeys, horses and mules, where a number of studies have implicated overworking/overloading and improper saddle/harness design as major causal factors for external lesions. similar welfare concerns are likely relevant to working elephants. consequently, the objectives of this study were to determine the nature of equipment-related skin injuries in a population of asian elephants working in the tourism industry by determining the following: 1) the prevalence of external injuries in four saddle related areas: neck, girth, back and tail, among working tourism elephants in thailand; and 2) potential risk factors for skin lesions related to working conditions, including demographic factors and saddle design. methods data for this cross-sectional study were collected between may and july by the lead researcher who worked in partnership with the elephant hospital of the national elephant institute, formerly known as the thai elephant conservation centre. data were collected from camps across thailand that received routine visitation from veterinarians from the livestock and wildlife hospital of mahidol university, or the national elephant institute’s elephant hospital. elephants were eligible for the study if they worked in the tourism industry and wore a saddle. each elephant underwent a veterinary examination and elephant keepers answered a series of standardized questions pertaining to each elephant and its saddle-wearing history using a translator. the outcome of interest for the present study was the presence of an ‘active lesion’ in any one of four anatomical areas: neck, girth, back and tail. lesions were placed initially into categories. however, to achieve adequate power for statistical analyses, lesions were classified as either ‘active’ or ‘inactive/absent’. active lesions were those that could be described as rubbed/pink, raw, a full depth ulcer, an abscess, a healing lesion, or depigmented skin. a body region was categorized as having an inactive/absent lesion if there was no lesion, a healed lesion/scar, or a callus. the following independent variables were investigated to determine their association with the presence of an active lesion: age of the elephant; region of the body; number of hours worked per day; the provision of a break; weight of the saddle; and contact material. to allow for a greater number of observations per category, shredded bark and bark were combined into a single “bark” category. descriptive analyses were carried out using microsoft office excel, and stata mp. the prevalence proportions of elephants with at least one lesion, with at least one active lesion, and with no lesions were estimated. the prevalence of active lesions was also summarized by body region. among camps, the prevalence of elephants with at least one active lesion was estimated. explanatory variables of interest were summarized and tabulated by camp. using stata mp, univariable and multivariable multilevel logistic regression models were built to examine the association between the presence of an active lesion at any of the four anatomical sites and the explanatory variables: age, working hours, contact material, break, and saddle weight. these models contained random intercepts for camp and elephant to account for clustering, as multiple elephants were observed per camp, and multiple body regions were examined per elephant. the models were used to determine variance estimates at both the elephant and camp level. decisions to include random intercepts were made based on model fit using akaike’s information criterion score, in combination with a likelihood-ratio test. correlation among independent variables was tested using spearman’s rank correlations. any two variables found to be highly correlated at r ≥  were assessed for completeness of data and biological relevance. this information was used to determine the most reliable variable to be included in the model, to avoid issues associated with collinearity. the linearity assumption was examined for continuous explanatory variables using locally weighted regression with the lowess command. if the lowess smoother indicated a curvilinear relationship, a quadratic term was investigated and remained in the model if α ≤ . variables with p ≤  on univariable analysis were considered for inclusion in the multivariable model. variables included in the final multilevel multivariable model showed either statistical significance, or evidence of acting as a confounding variable. interaction effects were initially examined, but due to estimation issues associated with relatively small sample size, we limited our analysis to a main effects model. pearson residuals and best linear unbiased predictors were examined to assess the fit of the multilevel multivariable model.  descriptive statistics in total, body areas, from elephants and camps were examined for saddle related lesions. the prevalences of elephants with at least one active lesion, at least one lesion, and no lesions were found to be  %,  %, and  %, respectively. the distribution of active lesions by body area was as follows:  % on the neck,  % on the girth,  % on the back, and  % on the tail. among camps, the average prevalence of elephants with at least one active lesion was  %, although this value ranged from to 100 %. the median value of prevalence of elephants with at least one active lesion by camp was  %. the number of elephants per camp ranged from to. the age of the elephants in years ranged from to with an average age of  years. scheduled working hours were found to be relatively consistent among elephants within a camp, although this value ranged between to 10 hours among camps. whether or not the elephants were provided a break was also found to be a camp level variable. saddle weight was found to be equivalent for all elephants among the same camp, with 50 % of the camps using heavy saddles, and 50 % using light saddles. the contact material of the saddle pad was found to vary both among camps, and between elephants within the same camp.table prevalence of elephants with at least one lesion, compared to prevalence of elephants with at least one active lesion by tourism camp in thailand 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 *an active lesion is defined as a lesion categorized as rubbed/pink, raw, full depth ulcer, abscess, healing lesion, depigmented skin, or other. an inactive lesion is defined as a lesion categorized as either a healed lesion, or a callus ^exact 95 % confidence interval #one-sided, exact  % confidence interval 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 nd = no data *break is defined as a period of time when either the entire saddle was completely removed, or the girth strap was loosened **saddle weight is dichotomized into light and heavy ***the individual elephant that works 2 hours/day gets the break ^each elephant had only one contact material, but the same material was not necessarily used for all animals within a camp risk factor analysis initially, the multilevel model included random intercepts for both elephant and camp. however, the random intercept for elephant was removed because the variance component was negligible and model fit, based on aic and a likelihood-ratio test, was not improved by its inclusion. it should be noted that the neck region variable was removed from both univariable and multivariable analyses. there were so few lesions found in this region that it resulted in estimation problems for coefficients in our multilevel multivariable model. i) univariable analysis based on univariable analysis, active lesions were more likely to be found in the back region in comparison to the other anatomical areas. increasing age was also identified by univariable analysis as being a risk factor for having an active lesion. there appeared to be a significant quadratic relationship between working hours/day and the log odds of having an active lesion. there was no significant difference in the prevalence of active lesions associated with the use of a heavy or light saddle, or with whether or not the elephants were provided a break. we observed no significant differences between the use of the three contact materials compared to bark.table univariable multilevel* logistic regression models of factors associated with the presence of an active lesion on elephants from tourism camps in thailand ref = referent category n = number of elephants with data collected for the indicated variable m = number of camps with data collected for the indicated variable *models include a random intercept for camp **break is defined as a period of time when either the girth strap was loosened, or the whole saddle was completely removed ii) multivariable analysis in the final multilevel multivariable model, the following variables were included: body region, contact material, saddle weight, the provision of a break, age, and working hours/day and working hours/day squared. as observed in the univariable analysis, the back region was found to be at higher odds for having an active lesion compared to the girth and tail regions. in comparison to each of the bark, carpet, and blanket materials, the use of a rice sack was found to be a significant risk factor for having an active lesion. further comparisons between the bark, carpet, and blanket materials revealed no significant difference in the odds of having an active lesion. increasing age was found to be associated with increased odds of having an active lesion. a significant quadratic term was again observed for working hours/day. initially the odds of having an active lesion increased with working hours/day; however this relationship peaked at 6– hours, followed by a decline. after controlling for the other variables in the model, receiving a break significantly increased the odds of having an active lesion. no significant difference was found between heavy and light saddles in association with the odds of having an active lesion. however, this variable was found to act as a confounder for both contact material and age, and was retained in the final model.table multivariable multilevel* logistic regression model of factors associated with the presence of an active lesion on sites on elephants working in tourism camps in thailand body area contact material saddle weight break** age working hours/day variance =  *model includes a random intercept for camp **break is defined as a period of time when either the girth strap was loosened, or the whole saddle was completely removed note: odds ratios are calculated using the second listed variable as the referent category *model used for comparisons includes a random intercept for camp fig. relationship between the predicted log odds of an active lesion and hours worked per day. predicted log odds of active lesion versus number of hours worked per day based on final multilevel multivariable logistic regression model for factors associated with active lesions on elephants from tourism camps in thailand an intercept only analysis revealed that 20 % of the model variance was explained at the camp level. in the final multilevel multivariable analysis, we found 13 % of the variance in the outcome was accounted for at the camp level. we observed no outliers based on examination of the pearson residuals. based on visual examination of the blups, we found that model assumptions concerning normality and homogeneity of variance were met. discussion the prevalence and associated risk factors for skin lesions were examined using data from working elephants from thailand tourism camps. descriptive analyses revealed a high prevalence of active lesions among these elephants especially in the back region. risk factors associated with increased odds of having an active lesion were identified as: body region, the use of rice sacks as padding material in contact with the skin, increasing age of the elephant, the provision of a break, and longer working days. a high prevalence of active lesions associated with working conditions was found, as  % of these elephants had at least one active lesion, and  % had at least one lesion. this finding is comparable to the prevalence of pack wounds found in various studies on working equids. interestingly, virtually none of the model variance was explained at the elephant level, suggesting lesions among sites on an elephant occur relatively independently of each other. approximately 13 % of the model variance was explained at the camp level, even after accounting for several camp level variables. these variables included working hours/day, the provision of a break, and saddle weight. this suggests that camp level management has a moderate impact on elephant welfare associated with external injuries, potentially indicating an area for more detailed investigation in future studies. the back region appears to be a high risk area for the presence of an active lesion relative to the other anatomical sites, followed by the girth region. these findings correspond to the results of various studies regarding working donkeys and mules. sells et al. identified the withers and shoulder as the regions of highest lesion prevalence, at 40 % and 31 %, respectively. biffa and woldemeskel saw the highest lesion prevalence on the withers at  %. pearson et al. documented the back region to be the area of highest lesion prevalence at 90-93 %, and the neck region as a lower risk area at 4-10 %. pritchard et al. saw lesions on the head, neck, ribs, flank and tail base in less than 10 % of equids studied, and observed the highest prevalence of lesions in mules on the breast/shoulder, withers and girth regions. the results of our study, along with these comparable equid results, suggest a primary area for intervention. future studies and current management/veterinary care regarding working elephants should focus efforts primarily on the back region and secondarily on the girth region, as opposed to the neck and tail, where lesions were relatively rare. padding material in contact with the skin was a significant factor for the prevalence of active lesions. in particular, rice sacks significantly increased the odds of having an active lesion relative to all other padding materials. avoiding their use in favour of bark, carpet, or blankets would be beneficial in terms of decreasing the odds of having an active lesion. bark, carpet and blankets were not found to be significantly different from one another; therefore, at this time there are no recommendations for a potential standard material. increasing age had an expected association with increasing odds of having an active lesion. although an odds ratio of is relatively small, if we take into account that these are long lived animals, the odds of an active lesion increases times for every 15 year increase in age. due to the nature of a cross-sectional study, there is an issue concerning prevalence vs. incidence of active lesions on older animals. the study design did not allow us to discern whether older animals were more likely to have active lesions due to increased susceptibility, decreased healing times, or some combination of both. regardless, at this time, we recommend that a higher level of monitoring/attention be given to older elephants by their owners or caretakers. it is important to note that due to the inclusion of the age variable, which had the largest amount of missing data, the number of elephants included in the final model decreased from to however, the direction of the odds ratios of the other variables remained the same; only the magnitude of these effects was altered. as we thought age was an important biological variable that showed statistical significance, we decided to include it in the final model despite the loss of some observations for the multivariable analysis. surprisingly, the provision of a break was found to be a risk factor associated with increased odds of having an active lesion. it is possible that camps that offer more breaks may be doing so in response to a current or previous problem with lesions. but again, as a result of the limitations of a cross-sectional study, the incidence of active lesions as a result of breaks could not be investigated directly. control trials or cohort studies should be conducted to investigate the effect of breaks on the incidence of lesions along with more specific descriptions concerning their timing and the conditions associated with these breaks. for instance, the loosening of the girth strap during a break may cause more movement of the saddle equipment and, therefore, more abrasions. the association between the odds of having an active lesion and working hours/day had a quadratic relationship. the predicted log odds of an active lesion were observed to initially increase with increasing working hours/day, followed by a peak at 6– hours and subsequent decline. this drop was an unexpected finding, but may be related to the “healthy worker bias”. for instance, more resilient elephants may be selected more frequently to work more hours, or those elephants that sustain too many lesions may have decreased numbers of working hours, or be removed from the camp. this could result in a selection bias and an impression of a healthier population of elephants working longer hours, which may explain the quadratic relationship identified. sells et al. found a similar result in their study analysing working equids in morocco, where increased frequency of work/week was found to be associated with a decreased risk of sustaining a pack wound. they also identified a healthy worker selection bias as a possible explanation for these findings. due to the nature of the present study, the prevalence of skin lesions was observed only in animals currently working in the camp, and data regarding total working hours/day were gathered only for current management. therefore, to control for this potential bias, future studies analysing the role of working hours should examine all elephants that have been used by the camp, and take information on any changes in working hours among the elephants. it may be prudent at this time to recommend limiting work to <6 hours/day since the odds of observing an active lesion was highest at this time point. although the weight of the saddle was found to be a confounding variable for the present study, a significant difference between heavy and light saddles was not found to be associated with an increased odds of having an active lesion. however, multiple studies investigating similar factors among working equids have implicated overloading as a major causal factor for skin lesions. kontogeorgopoulos also cautions that although these elephants can carry up to 300 kg, it is recommended they carry no more than 200 kg. therefore, the total weight that elephants carry should be taken into consideration. it is possible that relative to the individuals being carried, the weight of the saddle is not significant. it is also possible that the lack of variability in the weight variable within camps may have impacted our ability to find an effect of weight on lesion outcome. it is also important to note that by dichotomizing the weight of the saddles due to limitations in their assessment in the field, a great deal of information was lost that could have been identified if modelled as a continuous variable.  in conclusion, lesions related to working conditions were found to be highly prevalent among the tourism elephants studied. these lesions were found to be associated with risk factors including the use of rice sacks as padding material in contact with the skin, longer working days, and the provision of a break. the back region was identified as being at an increased risk of having an active lesion, followed by the girth region. while the cross-sectional nature of the study posed several inferential limitations, at this time we recommend the following until longitudinal studies can be conducted: increased monitoring of older elephants, increased monitoring of the back region of elephants of all ages, working hours limited to less than 6 hours per day, and the avoidance of rice sacks as padding material in contact with the skin. additional file additional file 1: elephant observation sheet. competing interests the authors declare that they have no competing interests. authors’ contributions sm, ta, and das planned and designed the study. sm led the collection and recording of all data. os and dlp conducted the statistical analyses and wrote the initial draft manuscript. all authors were involved in writing and editing the final manuscript. all authors read and approved the final manuscript. 
working hours had a quadratic relationship with the log odds of an active lesion where the probability of an active lesion initially increased with the number of working hours per day and then declined possibly reflecting a “healthy worker” bias where only animals without lesions continue to be able to work these longer hours. the objectives of the present study were to determine the prevalence of cutaneous lesions in anatomic regions in contact with saddle-related equipment among elephants in thailand working in the tourism industry, and to identify potential risk factors associated with these lesions. there was a high prevalence of active lesions, most often located on the back region. specifically, we recommend the following until longitudinal studies can be conducted: increased monitoring of older elephants and the back region of all elephants, working less than 6 hours per day, and the avoidance of rice sacks as padding material in contact with skin. while we recognize that the cross-sectional nature of the study posed some inferential limitations, our results offer several potential intervention points for the prevention of these lesions. using multilevel multivariable logistic regression modelling containing a random intercept for camp we identified the following risk factors: increasing elephant age, the use of rice sacks as padding material in contact with the skin, and the provision of a break for the elephants. data for this cross-sectional study were collected between may and july on elephants from tourism camps across thailand.
 lesions related to working conditions and improper saddle design are a concern for a variety of working animals including elephants. the objectives of the present study were to determine the prevalence of cutaneous lesions in anatomic regions in contact with saddle-related equipment among elephants in thailand working in the tourism industry, and to identify potential risk factors associated with these lesions. data for this cross-sectional study were collected between may and july on elephants from tourism camps across thailand.  there was a high prevalence of active lesions, most often located on the back region. using multilevel multivariable logistic regression modelling containing a random intercept for camp we identified the following risk factors: increasing elephant age, the use of rice sacks as padding material in contact with the skin, and the provision of a break for the elephants. working hours had a quadratic relationship with the log odds of an active lesion where the probability of an active lesion initially increased with the number of working hours per day and then declined possibly reflecting a “healthy worker” bias where only animals without lesions continue to be able to work these longer hours.  while we recognize that the cross-sectional nature of the study posed some inferential limitations, our results offer several potential intervention points for the prevention of these lesions. specifically, we recommend the following until longitudinal studies can be conducted: increased monitoring of older elephants and the back region of all elephants, working less than 6 hours per day, and the avoidance of rice sacks as padding material in contact with skin. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords elephantssaddle soresthailandtourismmulti-level modellingissue-copyright-statement© the author 2015 