Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Vet. Sci., 20 November 2025

Sec. Animal Behavior and Welfare

Volume 12 - 2025 | https://doi.org/10.3389/fvets.2025.1682443

This article is part of the Research TopicAdvancing Animal Welfare Assessment: From Biomarkers to Smart Monitoring TechnologiesView all articles

Use of thermal imaging to collect fresh faeces for non-invasive evaluation of stress levels in the European hare (Lepus europaeus)

Vlastimil Skotk,Vlastimil Skoták1,2Jan Cukor,
Jan Cukor1,3*Marek Sedl
ek,Marek Sedláček1,4Richard &#x;ev
íkRichard Ševčík1Rostislav LindaRostislav Linda1Kate&#x;ina BrynychovKateřina Brynychová1Matj Kostka,Matěj Kostka1,2Sabine Macho-MaschlerSabine Macho-Maschler5Jan Hu&#x;ek,Jan Hušek1,6Rupert PalmeRupert Palme5
  • 1Forestry and Game Management Research Institute, Jíloviště, Czechia
  • 2Faculty of Forestry and Wood Technology, Mendel University in Brno, Brno, Czechia
  • 3Faculty of Forestry and Wood Sciences, Czech University of Life Sciences in Prague, Prague, Czechia
  • 4Central European Institute of Technology, Masaryk University, Brno, Czechia
  • 5Experimental Endocrinology, Department of Biological Sciences and Pathobiology, University of Veterinary Medicine, Vienna, Austria
  • 6Department of Zoology, National Museum of the Czech Republic, Prague, Czechia

Assessing long-term stress in wild animal populations is extremely complicated and, in some species, practically impossible due to the complexity of sampling. Here, we tested and verified a unique non-invasive method for collecting fresh faecal samples located using thermal imaging cameras from a model species, the European hare. Subsequent analysis of faecal glucocorticoid metabolites (fGCMs) allows for the determination of stress levels without capture of individuals in hare populations. fGCM values ranged from 1.8 to 65.8 ng/g. Whilst the average value across locations in winter was 15.6 ng/g ± 9.9 SD, values in spring (April/March) were significantly higher (18.2 ng/g ± 11.2 SD; p < 0.05). These higher values coincided with the peak in European hare reproduction in spring. Significantly higher values were also confirmed in urban environments, indicating increased stress levels compared to natural environments, despite hares appearing to have adapted to the urban landscape. Higher values were also found in structurally poorer landscapes than Austrian rich agroecosystems. Having proved the usefulness of thermal imaging cameras for collecting large numbers of fresh droppings for the non-invasive evaluation of stress in wild hare sub-populations, we suggest the method could be applied to other species where capture and handling exert stress, injury or mortality.

Introduction

Whilst the population density of wild animal species is usually directly related to the quality of the environment in which they live, it can be difficult to assess the direct impact of environmental quality on population density as it is influenced by many factors (1). One method that has been used to assess environmental suitability and population success, particularly in endangered species, is to compare stress levels between individuals (ind.), populations or sub-populations, particularly as regards the effects of long-term stress caused by declines in environmental quality (2, 3). Endocrinological studies that monitor stress hormone concentrations in blood as an indicator of physiological stress are common (4) and are used for a wide range of wild and domestic animal species and for assessing the welfare of animals kept in captivity. However, the method tends to be far more successful in the case of domestic animals, which are easier to handle.

Most methods today analyse stress hormones such as cortisol and corticosterone, steroid hormones produced by the adrenal cortex and derived from cholesterol (5). These glucocorticoids (GCs) directly affect sugar metabolism and thus play an important role in controlling the metabolism of almost all tissues in the body (6). The dominating glucocorticoid (cortisol or corticosterone) differs between species, where most mammals including the European hare (Lepus europaeus) have cortisol (7).

A wide range of sample types may be used to monitor GC levels, these typically include blood, saliva, excrement, milk, fur, feathers and eggs (8, 9). Unfortunately, the sampling process itself can cause considerable stress to wild animals as it usually involves capturing and handling. Furthermore, blood sampling can cause acute stress, which will then be reflected in the levels of stress hormones in the blood, giving a false indication of environmental stress (10). In such cases, it would be more appropriate to use a non-invasive sampling method to determine GC and their metabolites (11). However, a major limitation of non-invasive sampling is the difficulty of reliably identifying freshly deposited faeces. As shown by Lafferty et al. (12), exposure to environmental conditions can rapidly alter glucocorticoid metabolite concentrations, which highlights the importance of collecting only fresh samples to ensure reliable results. Collection of an animal’s faeces offers a particular advantage in that they are easy to collect. However, faecal (and urine) samples can be difficult to (a) identify to species level, and (b) link to a specific individual. Furthermore, circulating hormone levels in faecal samples are integrated over time and thus are representative of cumulative hormone secretion (9). Consequently, most modern studies now monitor faecal GC metabolite (fGCM) concentrations when using faecal (or urine) samples as a non-invasive method to evaluate stress in animal welfare assessments (1316).

Stress levels will be affected by several factors, including environmental issues, which can have a significant impact on reproduction and population dynamics, especially in endangered species (17). Over the last century landscapes have been profoundly affected by human management worldwide, especially those associated with agroecosystems (18), which require a constant supply of energy in the form of fertilisers, pesticides and other agrotechnical interventions (19). These, often huge, agricultural operations have had a profound impact on wild animal populations, often representing the main limiting factor in present population dynamics (20, 21). Animal populations mainly affected by such operations include ground-nesting field birds, such as the grey partridge (Perdix perdix), the northern lapwing (Vanellus vanellus) and corncrake (Crex crex) (22), and mammals such as the European hare, whose population dynamics show a clear negative reaction to changes in the agricultural landscape (23). Though a relatively small species, the hare has relatively high environmental demands and, consequently, has been classified as a bioindicator species (24, 25). Just 60 years ago, the hare was a very common species in Central Europe (26), as documented by hunting catches of more than 1 million individuals per year in the Czech Republic alone. In recent decades, agricultural changes (e.g., a shift to large ‘prairie’ field systems and crop homogenisation) have had a demonstrable negative impact on hares, with numbers declining throughout Central Europe, even to the point of virtual extinction in some area (2729). Despite this, there is a lack of information regarding the impact of environmental stress on European hare (sub-) populations, most such studies having focused on stress monitoring in snowshoe hares (Lepus americanus) (3033) or mountain hares (Lepus timidus) (34, 35). In all studies, cortisol (or its faecal metabolites) has been monitored as the prime stress indicator in hares (3640).

The aim of this study was to (i) test and verify a new method for identifying fresh faeces (infra-red monitoring), which can then be used to assess European hare stress levels based on analysis of fGCM concentrations (i.e., metabolites of cortisol); (ii) to compare stress levels in sub-populations from different agricultural landscapes with differing hare densities; (iii) to compare stress levels in hares from agricultural environments with those from urban environments, and (iv) to compare stress levels during the period of food scarcity (December/January) and the breeding season (February/March). We hypothesise that fGCM concentrations will decrease as landscape suitability (i.e., heterogeneity) increases.

Materials and methods

Study area

Samples were obtained from two agricultural landscapes and one urbanised area in the Czech Republic and an additional agricultural site in Austria (Figure 1). The first agricultural site, located close to the village of Cerhenice (50.060278 N, 15.066389E; 240 m a.s.l.; Czech Republic), has a relatively homogeneous landscape and is characterised by conventional modern farming practices. In 2023, the average field size in this area was 13.13 ± 13.51 SD ha. The proportion of arable land is relatively high at 87.15%, with other non-agricultural landscape elements, such as wildflower strips (0.73%), almost absent (41). The second agricultural site, situated 1.1 km from the first site near the village of Ratenice (50.098611 N, 15.058750E; 188 m a.s.l.; Czech Republic), is managed using extensive farming practices and is characterised by a more heterogeneous landscape. The average field block size is smaller (4.64 ± 5.03 SD ha), the proportion of arable land reduced (83.01%) and the proportion of features such as wildflower strips higher (4.68%). The third agricultural site, located near Wildendürnbach (48.779167 N, 16.493611E; 193 m a.s.l.; Lower Austria), lies in a highly heterogeneous landscape with significantly smaller strip field blocks (1.78 ± 1.70 SD ha) with higher edge density, a lower proportion of arable land (64.83%) and a higher proportion of wildflower strips (5.93%) and other non-cropped elements (41). The final urbanised site is situated in the city of Prague (50.087500 N, 14.421389E; 192 m a.s.l.; Czech Republic), a large Central European city with a population of approximately 1,350,000 people.

Figure 1
Map of the Czech Republic showing elevation levels with color coding. Three inset images at a 1:10,000 scale display agricultural landscapes corresponding to study areas: Cerhenice (brown dot), Ratenice (red dot), and Wildendürnbach (yellow dot). Prague is marked with a black dot. A legend describes the elevations in meters above sea level, ranging from fifty to sixteen hundred.

Figure 1. Location and landscape structure of the study areas.

Note that, whilst fGCM concentrations were determined for the urban location, they were not compared with those from the other sites due to the clear difference in environment type and the absence of current data on hare population density in such urban environments. Previous monitoring efforts have recorded average adult hare spring population densities in more heterogeneous localities (Ratenice, Wildendürnbach) are significantly higher (mean density ± SD = 117 ± 68 and 114 ± 49 ind./100 ha, respectively) than population density (30 ± 24 ind./100 ha) in more homogeneous locality (Cerhenice) (41).

European hare pellet samples were obtained during night patrols (i.e., the natural activity period of the species, between 21:00 and 02:00 during two consecutive nights at each site) at known activity centres with the use of Merger LRF XP50 thermal imaging binoculars (640 × 480 pixel microbolometer, NETD < 25 mK, F50/1.0 lens, 2.5–20 × magnification; Pulsar, Lithuania), which enabled real-time detection of both individuals and fresh pellets based on their heat signature, which made them clearly distinguishable from older pellets. The pellets were usually obtained from the hare’s last resting place as defecation inevitably occurred immediately after the individual was disturbed by the approach of a researcher. This method allowed for the collection of fresh droppings only, thereby minimising the potential risk of steroid degradation due to weather exposure. If multiple faecal pellets were produced by the same individual, all pellets were collected, thoroughly homogenized, and treated as one composite sample representing that individual (weight of the sample is specified below). Upon collection, the samples were immediately placed in test tubes, frozen at −20 °C in a field freezer and then stored at −20 °C until laboratory analysis.

Sample analysis

For quantification of fGCMs, the samples were first dried at 75 °C for 5 h and then mechanically homogenised. For extraction, 0.20 g (± 0.005 g) of dry matter was mixed with 4.0 mL of 100% methanol and 1.0 mL of distilled water. The samples were then placed on a shaker for 30 min and subsequently centrifuged at 2500 g for 10 min, after which the supernatant was diluted at a 1:10 ratio with test buffer (13). Determination of fGCMs was performed using a group-specific 11-oxoaetiocholanolone enzyme immunoassay, which detects 11,17-dioxoandrostanes, a group of cortisol metaolites (13, 36, 40). For statistical evaluation purposes, samples with concentrations below the detection limit were set to 1.8 ng/g.

Statistical analysis

For each location, the relationship between hare abundance (ind./100 ha) and (a) mean size of land patch, and (b) fGCM concentration across both study seasons (December/January and February/March), were assessed using the non-parametric Kruskal-Wallis test (assumptions of ANOVA not met), with subsequent multiple comparisons. Likewise, as assumptions of sample heterogeneity were not met in any sample, fGCM concentrations between study seasons for each location separately and study seasons across all locations were compared using the Wilcoxon rank-sum test, with results presented as boxplots. The correlation between acreage of particular field blocks and hare abundance was assessed using Pearson correlation coefficient and respective statistical test (test for correlation between paired samples).

All statistical procedures were performed using R software (42), with plots created using the R package “ggplot2” (43). An alpha level of 0.05 was used for all statistical tests.

Results

A total of 373 faecal samples were analysed, comprising 73 from Cerhenice, 132 from Prague, 59 from Wildendürnbach and 109 from Ratenice. Concentrations of fGCM ranged from 1.8 to 65.8 ng/g (average = 16.68 ng/g), though high values were relatively rare (Figure 2). Whilst values below the minimum detection limit of 1.8 ng/g were recorded at all locations except Prague, Prague had the maximum fGCM concentration of 65.8 ng/g recorded in the Czech Republic (spring), the maximum recorded in Wildendürnbach was 37.6 ng/g (Figure 3).

Figure 2
Bar graph showing the frequency of fecal glucocorticoid metabolites (fGCMs) in nanograms per gram across two seasons: December/January (black) and February/March (gray). The December/January data peaks at a frequency over one hundred for fGCMs around ten nanograms per gram, while February/March data peaks under fifty for similar levels, with both decreasing steadily as fGCMs increase.

Figure 2. Histogram of fGCM concentrations detected across all locations.

Figure 3
Four bar charts display frequency of fGCMs in ng/g for locations: Cerhenice, Wildendürnbach, Prague, and Ratenice. Bars represent two seasons, December/January and February/March. Each chart shows varying fGCM distribution.

Figure 3. Histograms of fGCM concentrations detected at each location.

Lowest (mean ± SD) fGCM values were recorded in Wildendürnbach (10.8 ± 8.3 ng/g), whilst highest values were recorded in Prague (21.0 ± 10.5 ng/g), i.e., around twice those in Wildendürnbach. Comparable concentrations were recorded at the other Czech localities, with an average of 15.2 ± 11.5 ng/g at Cerhenice and 15.6 ± 8.8 ng/g at Ratenice. Overall (all sites combined; data not shown), fGCM concentrations were significantly higher in spring (winter: 15.6 ± 9.9 ng/g, spring: 18.27 ± 11.2 ng/g; W = 14,662, p < 0.05). However, no consistent seasonal trend was observed between localities, with concentrations comparable between seasons at Cerhenice and Ratenice and higher values in spring at Wildendürnbach and Prague (Figure 4), though the difference was only significant at Prague (p < 0.05).

Figure 4
Box plot showing fGCMs levels (nanograms per gram) across four localities: Cerhenice, Wildendürnbach, Prague, and Ratenice, comparing Dec/Jan and Feb/Mar seasons. Each locality displays two boxes representing different seasons, highlighting variations in stress markers. Outliers are depicted as dots.

Figure 4. Boxplots of fGCM concentrations by location and sampling season. Symbols above the groups indicate statistical significance: * = p < 0.05, ○ = marginally significant (0.05 < p < 0.1), no symbol = no significant difference.

A comparison of average field block size, fGCM concentration (across sampling dates) and hare abundance at the three rural sites (urban location excluded; Figure 5) indicated a high negative correlation between field block size and hare abundance (r = −0.975), though the result was not statistically significant due to the low number of sites assessed (p = 0.14). There was a highly significant difference in fGCM concentrations between sites (Kruskal-Wallis test, chi-squared = 15.499, df = 2, p < 0.001), with multiple comparisons showing values from Wildendürnbach to be significantly lower than those at Czech sites, which showed no significant difference.

Figure 5
Box plot comparing mean land patch size in hectares and hare abundance in individuals per 100 hectares across three localities: Cerhenice, Ratenice, and Wildendürnbach. Cerhenice shows the largest land patch size variance, with red asterisks indicating outliers in hare abundance for Cerhenice and Wildendürnbach.

Figure 5. Average field block size (columns, primary axis on left), hare abundance (red stars, secondary axis on right) and fGCM concentration (boxes, secondary axis on right) at each of the rural localities.

Discussion

Over the past 20 years, the determination of fGCMs has become an increasingly popular method for assessing stress in animals, because faeces collection is non-invasive and avoids additional stress due to capture and handling (8). However, faeces need to be collected fresh, because fGCMs were found to degrade easily (16), which is especially challenging in free ranging, wild animals. In this study, we tested a new method for ensuring the collection of fresh wild hare pellets, and used those to determine environmental suitability for hares by relating fGCM concentrations (as an indicator of stress) (35, 44) to the density of a given sub-population.

In the case of European hare, studies on captured hares have shown that the related acute stress poses a particularly serious risk as it can lead to trauma or mortality (45). Such influences can be minimised, if fGCMs are measured in samples collected without disturbing the animal. However, age of faecal samples can be crucial in ensuring accurate fGCM readings as levels may change over time due to degradation by bacterial enzymes (16, 43, 4648). Time after defecation and environmental conditions have been found to exert a strong influence on measured fGCM levels (16, 43, 49). Therefore, the generally accepted advice is to collect fresh faeces (13). This may, however, be a serious problem in wildlife studies. Here we report for the first time the use of thermal imaging camera not only for detecting the animals, but also the freshly deposited faecal pellets. By using mobile freezer boxes, we were also able to freeze samples immediately on site, preventing any subsequent degradation (50).

Likewise, by sampling always at night with the thermal camera, a possible influence of a diurnal rhythm in GC secretion was avoided (8, 51). One limitation of the thermal camera method is that it is not possible to determine the sex of individuals due to a lack of sexual dimorphism in the case of European hare, and this may affect GC concentrations (8, 35, 50). This will not be an issue, however, if the method is applied to other wild species with more obvious sex dimorphism, easily observable with a thermal camera. For example, in some mammalian species such as red deer (Cervus elaphus), fallow deer (Dama dama), or mouflon (Ovis musimon), infrared thermography can reliably distinguish males by the presence of antlers or horns and by differences in body morphology.

In our study, which evaluated 373 faecal samples from wild European hares, fGCM levels ranged from 1.8 to 65.8 ng/g, with an average concentration of 16.7 ng/g. These values correspond with those from a comparable study by Cybulska et al. (36), who monitored fGCM concentrations in European hares in an agricultural landscape in Wildendürnbach during the grain harvest. Using 591 faecal samples, they recorded concentrations ranging between 5 and 65.6 ng/g. In another study on captive European hares, Teskey-Gerstl (40) recorded baseline concentrations ranging from 43 to 274 nmol/kg, equivalent to 13.0 to 82.9 ng/g. By applying a stress stimulus (repeated disturbance of rest), they were able to induce metabolite levels up to five times higher than original baseline values, with maximum concentrations reached within 1–2 days of stress stimulus, after which levels gradually returned to baseline. However, that these values refer to cage-bred hares only, where levels of stress can be induced or influenced and subsequently monitored relatively easily (51). To determine levels of stress associated with artificial breeding. Janicki et al. (52) monitored fGCM levels in pellets of European hares kept as pairs in cages. They recorded average values of around 201 ng/g in May and 180 ng/g in November, with maxima exceeding 300 ng/g. Such concentrations were significantly higher than those from wild individuals. The authors attributing this to constant mild stress caused by the unnatural cohabitation of males and females in a confined space. Consequently, stress hormone concentrations from captive studies should not be compared directly with wild animal studies or those using other types of samples.

In the wild, stress levels are assessed to evaluate intraspecific stress or the welfare status of a given species in relation to environmental quality. In our study, we recorded average fGCM concentrations of 10.8 ng/g at a highly heterogenous site in Wildendürnbach, and a significantly higher value of 21.0 ng/g at the homogenous urban site in Prague, i.e., around double the Austrian average. Four outliers (> 50 ng/g) are likely explained by possible trauma, possibly resulting from attempted predation. In a study by Boudreau et al. (30), however, no increase in fGCM levels (against a control) was recorded in hares exposed to predation risk. Further analysis was based on a comparison of the rural sites only as it was not possible to determine hare population density accurately in Prague due to the highly urbanised nature of the city. Consequently, we base our findings solely on comparisons of the Wildendürnbach locality (avg. 10.8 ng/g), Cerhenice (avg. 15.15 ng/g) and Ratenice (avg. 15.63 ng/g). Whilst the population density of the European hare differed significantly (p < 0.001) between the Czech homogenous landscape (Cerhenice) and two heterogenous landscapes in the Czech Republic (Ratenice) and Austria (Wildendürbach) (41), the differences in metabolite levels were not significant. Thus, our results do not indicate that stress levels in hares are lower at localities with lower population density. Instead, stress level changes are more likely to be caused by unsuitable agricultural landscape structure (23, 44, 53). This is particularly clear when comparing urban environments like Prague with highly structured rural landscapes, which represent ideal environments for hares (26). Surprisingly, previous studies have recorded a relatively high hare population density in urban environments (ca. 8 ind./100 ha, including building areas) (54), with the hare showing considerable adaptability to human presence in urban environments, their flight distance being less than half of that in agricultural landscapes (55). Nevertheless, evaluation of stress hormones shows that, despite their apparent tolerance, individuals living in urban environments are likely to suffer from long-term, low-level stress. Our results also showed a notable effect of field block size on stress levels, with a trend of fGCM concentrations increasing as field block area increased. However, this relationship was only significant for the Wildendürnbach site (p < 0.05). This suggests that large land blocks, often associated with intensive agricultural management and low landscape diversity, increase stress levels in hares, probably due to the relative lack of shelter options and the higher risk of disturbance/predation.

In the present study, slight seasonal differences in fGCM levels were observed, with higher concentrations in February/March compared to December/January, most likely related to hares mating in spring. Whilst male hares are fertile for most of the year (except October and November), females are limited to an extended reproductive season stretching from January to August, peaking in March (56). At this point, females become highly receptive, prompting high levels of aggressive competitive activity, a phenomenon often referred to as “March madness” [e.g., see Lincoln (56)]. During this month, when 100% of adult females usually become pregnant, the high levels of competition are likely to result in both physical exhaustion and associated stress, especially in females. This is the most likely explanation for the higher stress levels observed in sub-populations with higher population densities observed in this study, where competition for females is likely to be more intense. As such, the results partially confirm our hypothesis that higher density populations show higher levels of stress, albeit seasonally. However, we compared only spring data with that from December/November, when hares may already be entering their breeding season. Increased fGCM concentrations in spring may partly reflect elevated physiological stress in pregnant females, potentially linked to adaptive maternal effects as shown in other mammal species (57). Thus, we suggest that further sampling may be appropriate in October/November, allowing a better separation of mating-related stress and environmental stress levels.

Conclusion

The results of this study confirm high quality thermal imaging binoculars as a valuable tool for identifying a model animal species, the European hare, at night in the wild, allowing for the non-invasive collection of fresh faecal pellets for determining fGCM, without the need for trapping or handling. As such, we highly recommend its use for collecting fresh faeces in other sensitive or endangered species where disturbance or handling could adulterate subsequent stress levels.

The fGCM values recorded indicated a relatively balanced level of adrenocortical activity in open agricultural landscapes, but higher stress values in urban environments. Whilst the hare shows a high level of adaptability in urban areas, they are clearly not comfortable in such environments. Further, we recorded lower stress levels in hares from a heterogenous extensively farmed landscape (Wildendürnbach) compared with more conventionally managed agricultural landscapes (Czech Republic), confirming that larger, more intensively managed field blocks represent less favourable wildlife environments. Finally, our study confirmed increased stress levels in hares during the spring reproductive period, though we were unable to separate data by sex and our comparative autumn data did overlap slightly with potential reproduction.

Overall, our study confirms the suitability of thermal imaging as a non-invasive monitoring tool and that fGCM assessment in faecal matter obtained in this way can be a sensitive indicator of physiological stress when monitoring the impact of landscape structure and agricultural management on wildlife species. Longer-term monitoring at more locations and at different times of year could be beneficial in gaining a deeper understanding of the relationships between stress physiology, environment and hare behaviour.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

Ethical approval was not required for the study involving animals in accordance with the local legislation and institutional requirements because The analysis was conducted using faecal samples collected by a non-invasive method.

Author contributions

VS: Data curation, Writing – original draft, Investigation, Validation, Conceptualization, Writing – review & editing, Methodology. JC: Writing – review & editing, Project administration, Methodology, Validation, Writing – original draft, Resources, Investigation, Supervision, Funding acquisition, Conceptualization. MS: Writing – original draft, Writing – review & editing, Data curation, Validation, Methodology. RŠ: Data curation, Writing – review & editing, Writing – original draft, Methodology. RL: Writing – original draft, Writing – review & editing, Formal analysis, Data curation, Visualization, Validation, Methodology. KB: Writing – original draft, Writing – review & editing, Data curation, Methodology. MK: Writing – original draft, Data curation, Writing – review & editing. SM-M: Writing – review & editing, Writing – original draft, Data curation. JH: Writing – original draft, Writing – review & editing. RP: Writing – original draft, Supervision, Methodology, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. The research was funded by the Ministry of Agriculture of the Czech Republic, institucional support MZE-RO0123, by the Technological Agency of the Czech Republic (SS05010238) and by by Czech University of Life Sciences Prague, Faculty of Forestry and Wood Sciences (Excellent Team 2025).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Shultz, MT, and Kitaysky, AS. Spatial and temporal dynamics of corticosterone and corticosterone binding globulin are driven by environmental heterogeneity. Gen Comp Endocrinol. (2008) 155:717–28. doi: 10.1016/j.ygcen.2007.11.002

PubMed Abstract | Crossref Full Text | Google Scholar

2. Wingfield, JC, O’reilly, KM, and Astheimer, LB. Modulation of the adrenocortical responses to acute stress in arctic birds: a possible ecological basis. Integr Comp Biol. (1995) 35:285–94. doi: 10.1093/icb/35.3.285

Crossref Full Text | Google Scholar

3. Wingfield, JC, and Hunt, KE. Arctic spring: hormone-behavior interactions in a severe environment. Comp Biochem Physiol B Biochem Mol Biol. (2002) 132:275–86. doi: 10.1016/s1096-4959(01)00540-1

PubMed Abstract | Crossref Full Text | Google Scholar

4. Romero, LM, and Romero, LM. Physiological stress in ecology: lessons from biomedical research physiological stress in ecology: lessons from biomedical research. Trends Ecol Evol. (2016) 19:249–55. doi: 10.1016/j.tree.2004.03.008

Crossref Full Text | Google Scholar

5. Holst, JP, Soldin, OP, Guo, T, and Soldin, SJ. Steroid hormones: relevance and measurement in the clinical laboratory. Clin Lab Med. (2004) 24:105–18. doi: 10.1016/j.cll.2004.01.004

PubMed Abstract | Crossref Full Text | Google Scholar

6. Kuo, T, McQueen, A, Chen, TC, and Wang, JC. Regulation of glucose homeostasis by glucocorticoids. Adv Exp Med Biol. (2015) 872:99–126. doi: 10.1007/978-1-4939-2895-8_5

PubMed Abstract | Crossref Full Text | Google Scholar

7. Botía, M, Escribano, D, Martínez-Subiela, S, Tvarijonaviciute, A, Tecles, F, López-Arjona, M, et al. Different types of glucocorticoids to evaluate stress and welfare in animals and humans: general concepts and examples of combined use. Meta. (2023) 13:106. doi: 10.3390/metabo13010106

PubMed Abstract | Crossref Full Text | Google Scholar

8. Sheriff, MJ, Dantzer, B, Delehanty, B, Palme, R, and Boonstra, R. Measuring stress in wildlife: techniques for quantifying glucocorticoids. Oecologia. (2011) 166:869–87. doi: 10.1007/s00442-011-1943-y

PubMed Abstract | Crossref Full Text | Google Scholar

9. Palme, R. Monitoring stress hormone metabolites as a useful, non-invasive tool for welfare assessment in farm animals. Anim Welf. (2012) 21:331–7. doi: 10.7120/09627286.21.3.331

Crossref Full Text | Google Scholar

10. Kim, S, Foong, D, Cooper, MS, Seibel, MJ, and Zhou, H. Comparison of blood sampling methods for plasma corticosterone measurements in mice associated with minimal stress-related artefacts. Steroids. (2018) 135:69–72. doi: 10.1016/j.steroids.2018.03.004

PubMed Abstract | Crossref Full Text | Google Scholar

11. Peter, ID, Haron, AW, Jesse, FFA, Ajat, M, Han, MHW, Fitri, WN, et al. Opportunities and challenges associated with fecal progesterone metabolite analysis. Vet World. (2018) 11:1466. doi: 10.14202/vetworld.2018.1466-1472

PubMed Abstract | Crossref Full Text | Google Scholar

12. Lafferty, DJ, Zimova, M, Clontz, L, Hackländer, K, and Mills, LS. Noninvasive measures of physiological stress are confounded by exposure. Sci Rep. (2019) 9:19170. doi: 10.1038/s41598-019-55715-5

PubMed Abstract | Crossref Full Text | Google Scholar

13. Palme, R, Touma, C, Arias, N, Dominchin, MF, and Lepschy, M. Steroid extraction: get the best out of faecal samples. Wien Tierarztl Monatsschr. (2013) 100:238–46.

Google Scholar

14. Prola, L, Cornale, P, Renna, M, Macchi, E, Perona, G, and Mimosi, A. Effect of breed, cage type, and reproductive phase on fecal corticosterone levels in doe rabbits. J Appl Anim Welf Sci. (2013) 16:140–9. doi: 10.1080/10888705.2013.768903

PubMed Abstract | Crossref Full Text | Google Scholar

15. Möstl, E. Glucocorticoids, their metabolites and their measurement in various animal species. Med Weter. (2014) 70:524–9.

Google Scholar

16. Palme, R. Non-invasive measurement of glucocorticoids: advances and problems. Physiol Behav. (2019) 199:229–43. doi: 10.1016/j.physbeh.2018.11.021

PubMed Abstract | Crossref Full Text | Google Scholar

17. Milligan, SR, and Lloyd, R. Impacts of climate change and environmental factors on reproduction and development in wildlife. Philos Trans R Soc Lond Ser B Biol Sci. (2009) 364:3313–9. doi: 10.1098/rstb.2009.0175

PubMed Abstract | Crossref Full Text | Google Scholar

18. Li, X, Tian, H, Lu, C, and Pan, S. Four-century history of land transformation by humans in the United States (1630-2020): annual and 1gkm grid data for the HIStory of LAND changes (HISLAND-US). Earth Syst Sci Data. (2023) 15:1005–35. doi: 10.5194/essd-15-1005-2023

Crossref Full Text | Google Scholar

19. Dhuldhaj, UP, Singh, R, and Singh, VK. Pesticide contamination in agro-ecosystems: toxicity, impacts, and bio-based management strategies. Environ Sci Pollut Res. (2023) 30:9243–70. doi: 10.1007/s11356-022-24381-y

PubMed Abstract | Crossref Full Text | Google Scholar

20. Pywell, RF, Heard, MS, Woodcock, BA, Hinsley, S, Ridding, L, Nowakowski, M, et al. Wildlife-friendly farming increases crop yield: evidence for ecological intensification. Proc R Soc Lond B Biol Sci. (2015) 282:20151740. doi: 10.1098/rspb.2015.1740

PubMed Abstract | Crossref Full Text | Google Scholar

21. Emmerson, M, Morales, MB, Oñate, JJ, Batary, P, Berendse, F, Liira, J, et al. How agricultural intensification affects biodiversity and ecosystem services In: Dumbrell AJ, Kordas R, Woodward G, (eds). Advances in ecological research, London: Academic Press. vol. 55 (2016). 43–97.

Google Scholar

22. Evans, KL. The potential for interactions between predation and habitat change to cause population declines of farmland birds. Ibis. (2004) 146:1–13. doi: 10.1111/j.1474-919X.2004.00231.x

Crossref Full Text | Google Scholar

23. Ševčík, R, Krivopalova, A, and Cukor, J. Přímý vliv struktury zemědělské krajiny na výměru domovských okrsků zajíce polního: předběžné výsledky z České republiky. Zpr Lésn Výzk. (2023) 68:159–67. doi: 10.59269/zlv/2023/3/701

Crossref Full Text | Google Scholar

24. Can, SÜ, Demirbaş, Y, Erduran, N, and Beuković, D. The European hare (Lepus europaeus) as a biomonitor of selected heavy metal pollution in the Marmara region, Türkiye. Contemp Agric. (2025) 74:35–44. doi: 10.2478/contagri-2025-0006

Crossref Full Text | Google Scholar

25. Marinković, D, Relić, R, Lakić, N, Aničić, M, Beuković, D, Davidović, V, et al. Influence of age, habitat elevation, and distance to a thermal power plant on pathomorphological findings in the European brown hare (Lepuseuropaeus P.). Vet Q. (2023) 43:1–10. doi: 10.1080/01652176.2023.2273887

PubMed Abstract | Crossref Full Text | Google Scholar

26. Smith, RK, Jennings, NV, and Harris, S. A quantitative analysis of the abundance and demography of European hares Lepus europaeus in relation to habitat type, intensity of agriculture and climate. Mamm Rev. (2005) 35:1–24. doi: 10.1111/j.1365-2907.2005.00057.x

Crossref Full Text | Google Scholar

27. Sliwinski, K, Ronnenberg, K, Jung, K, Strauß, E, and Siebert, U. Habitat requirements of the European brown hare (Lepus europaeus Pallas 1778) in an intensively used agriculture region (Lower Saxony, Germany). BMC Ecol. (2019) 19:31. doi: 10.1186/s12898-019-0247-7

PubMed Abstract | Crossref Full Text | Google Scholar

28. Smith, RK, Jennings, NV, Robinson, A, and Harris, S. Conservation of European hares Lepus europaeus in Britain: is increasing habitat heterogeneity in farmland the answer? J Appl Ecol. (2004) 41:1092–102. doi: 10.1111/j.0021-8901.2004.00976.x

Crossref Full Text | Google Scholar

29. Vizzarri, F, Slamecka, J, Sladecek, T, Jurcik, R, Ondruska, L, and Schultz, P. Long-term monitoring of European brown hare (Lepus europaeus) population in the Slovak Danubian lowland. Diversity. (2024) 16:486. doi: 10.3390/d16080486

Crossref Full Text | Google Scholar

30. Boudreau, MR, Seguin, JL, Boonstra, R, Palme, R, Boutin, S, Krebs, CJ, et al. Experimental increase in predation risk causes a cascading stress response in free-ranging snowshoe hares. Oecologia. (2019) 191:311–23. doi: 10.1007/s00442-019-04500-2

PubMed Abstract | Crossref Full Text | Google Scholar

31. Lavergne, SG, Smith, K, Kenney, A, Krebs, CJ, Palme, R, and Boonstra, R. Physiology and behaviour of juvenile snowshoe hares at the start of the 10-year cycle. Anim Behav. (2019) 157:141–52. doi: 10.1016/j.anbehav.2019.09.003

Crossref Full Text | Google Scholar

32. Lavergne, SG, Krebs, CJ, Kenney, AJ, Boutin, S, Murray, D, Palme, R, et al. The impact of variable predation risk on stress in snowshoe hares over the cycle in North America’s boreal forest: adjusting to change. Oecologia. (2021) 197:71–88. doi: 10.1007/s00442-021-05019-1

PubMed Abstract | Crossref Full Text | Google Scholar

33. Sheriff, MJ, Bosson, CO, Krebs, CJ, and Boonstra, R. A non-invasive technique for analyzing fecal cortisol metabolites in snowshoe hares (Lepus americanus). J Comp Physiol B. (2009) 179:305–13. doi: 10.1007/s00360-008-0314-4

PubMed Abstract | Crossref Full Text | Google Scholar

34. Rehnus, M, Wehrle, M, and Palme, R. Mountain hares Lepus timidus and tourism: stress events and reactions. J Appl Ecol. (2014) 51:6–12. doi: 10.1111/1365-2664.12174

Crossref Full Text | Google Scholar

35. Rehnus, M, and Palme, R. Yearly variations in gcm concentrations in female mountain hares (Lepus timidus) and the effect of pregnancy. Animals. (2021) 11:2710. doi: 10.3390/ani11092710

PubMed Abstract | Crossref Full Text | Google Scholar

36. Cybulska, N, Hackländer, K, Palme, R, Frey-Roos, A, and Schai-Braun, SC. Influence of cereal harvest on adrenocortical activity in European hares (Lepus europaeus). Eur J Wildl Res. (2020) 66:95. doi: 10.1007/s10344-020-01430-w

Crossref Full Text | Google Scholar

37. Esposito, L, Auletta, L, Ciani, F, Pelagalli, A, Pasolini, MP, Lamagna, B, et al. Hair cortisol levels in captive brown hare (Lepus europaeus): potential effect of sex, age, and breeding technology. Eur J Wildl Res. (2017) 63:62. doi: 10.1007/s10344-017-1121-6

Crossref Full Text | Google Scholar

38. Lavergne, SG, Peers, MJL, Mastromonaco, G, Majchrzak, YN, Nair, A, Boutin, S, et al. Hair cortisol as a reliable indicator of stress physiology in the snowshoe hare: influence of body region, sex, season, and predator–prey population dynamics. Gen Comp Endocrinol. (2020) 294:113471. doi: 10.1016/j.ygcen.2020.113471

PubMed Abstract | Crossref Full Text | Google Scholar

39. Pedersen, S, Ciesielski, TM, Stawski, C, Reid, N, Sait, STL, and Pedersen, HC. Hunting method affects cortisol levels in harvested mountain hares (Lepus timidus). Wildlife Biol. (2024):e01366. doi: 10.1002/wlb3.01366

Crossref Full Text | Google Scholar

40. Teskey-Gerstl, A, Bamberg, E, Steineck, T, and Palme, R. Excretion of corticosteroids in urine and faeces of hares (Lepus europaeus). J Comp Physiol B. (2000) 170:163–8. doi: 10.1007/s003600050271

PubMed Abstract | Crossref Full Text | Google Scholar

41. Ševčík, R, Cukor, J, Brynychová, K, Šálek, M, and Hruška, J. Agroecosystem quality as an essential driver of European hare (Lepus europaeus) population density in contrasting farmlands. Oikos. (2025) e11695. doi: 10.1002/oik.11695, (in press)

Crossref Full Text | Google Scholar

42. R Core Team. (2019) R: A language and environment for statistical computing. Vienna, Austria: R foundation for statistical computing.

Google Scholar

43. Wickham, H. ggplot2: Elegant graphics for data analysis. New York: Springer-Verlag (2016).

Google Scholar

44. Washburn, BE, and Millspaugh, JJ. Effects of simulated environmental conditions on glucocorticoid metabolite measurements in white-tailed deer feces. Gen Comp Endocrinol. (2002) 127:217–22. doi: 10.1016/S0016-6480(02)00056-4

PubMed Abstract | Crossref Full Text | Google Scholar

45. Mayer, M, Haugaard, L, and Sunde, P. Scared as a hare: effects of capture and experimental disturbance on survival and movement behavior of European hares. Wildlife Biol. (2021) 2021:wlb–00840. doi: 10.2981/wlb.00840

Crossref Full Text | Google Scholar

46. Osburn, KR, Crossey, B, Majelantle, TL, and Ganswindt, A. Examining alterations in fGCM concentrations post-defaecation across three animal feeding classes (ruminants, hindgut fermenters and carnivores). J Zool. (2025) 326:37–44. doi: 10.1111/jzo.13257

Crossref Full Text | Google Scholar

47. Donini, V, Iacona, E, Pedrotti, L, Macho-Maschler, S, Palme, R, and Corlatti, L. Temporal stability of fecal cortisol metabolites in mountain-dwelling ungulates. Sci Nat. (2022) 109:20. doi: 10.1007/s00114-022-01792-y

PubMed Abstract | Crossref Full Text | Google Scholar

48. Eriksson, H, and Gustafsson, J. Steroids in germfree and conventional rats: steroids in the mono- and Disulphate fractions of Faeces from female rats. Eur J Biochem. (1970) 16:252–2606. doi: 10.1111/j.1432-1033.1970.tb01079.x

Crossref Full Text | Google Scholar

49. Carbillet, J, Palme, R, Maublanc, ML, Cebe, N, Gilot-Fromont, E, Verheyden, H, et al. Instability of fecal glucocorticoid metabolites at 4°C: time to freeze matters. J Exp Zool A Ecol Integr Physiol. (2023) 339:625–32. doi: 10.1002/jez.2704

Crossref Full Text | Google Scholar

50. Touma, C, and Palme, R. Measuring fecal glucocorticoid metabolites in mammals and birds: the importance of validation. Ann N Y Acad Sci. (2005) 1046:54–74. doi: 10.1196/annals.1343.006

PubMed Abstract | Crossref Full Text | Google Scholar

51. Baias, A, Bodnariu, A, Nichita, I, and Cristina, R. Stress in laboratory juvenile rabbits: physiological indicators. Sci Pap Anim Sci Biotechnol. (2012) 45:142.

Google Scholar

52. Janicki, Z, Konjević, D, Pintur, K, Severin, K, Slavica, A, and Mašek, T. Non-invasive monitoring of cortisol metabolites level in farmed brown hare (Lepus europaeus). Vet Arh. (2006) 76:S251–7.

Google Scholar

53. Pavliska, PL, Riegert, J, Grill, S, and Šálek, M. The effect of landscape heterogeneity on population density and habitat preferences of the European hare (Lepus europaeus) in contrasting farmlands. Mamm Biol. (2018) 88:8–15. doi: 10.1016/j.mambio.2017.11.003

Crossref Full Text | Google Scholar

54. Mayer, M, and Sunde, P. Colonization and habitat selection of a declining farmland species in urban areas. Urban Ecosyst. (2020) 23:543–54. doi: 10.1007/s11252-020-00943-1

Crossref Full Text | Google Scholar

55. Krivopalova, A, Mikula, P, Cukor, J, Ševčík, R, Brynychová, K, and Šálek, M. Adaptation of farmland mammalian specialist to urban life: escape behavior of European hare along the urban-rural gradient. Sci Total Environ. (2024) 951:175779. doi: 10.1016/j.scitotenv.2024.175779

PubMed Abstract | Crossref Full Text | Google Scholar

56. Lincoln, GA. Reproduction and “march madness” in the Brown hare, Lepus europaeus. J Zool. (1974) 174:1–14. doi: 10.1111/j.1469-7998.1974.tb03140.x

Crossref Full Text | Google Scholar

57. Dantzer, B, Newman, AEM, Boonstra, R, Palme, R, Boutin, S, Humphries, MM, et al. Density triggers maternal hormones that increase adaptive offspring growth in a wild mammal. Science. (2013) 340:1215–7. doi: 10.1126/science.1235765

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: animal welfare, biodiversity conservation, cortisol, intrapopulation stress, non-invasive monitoring, wildlife monitoring

Citation: Skoták V, Cukor J, Sedláček M, Ševčík R, Linda R, Brynychová K, Kostka M, Macho-Maschler S, Hušek J and Palme R (2025) Use of thermal imaging to collect fresh faeces for non-invasive evaluation of stress levels in the European hare (Lepus europaeus). Front. Vet. Sci. 12:1682443. doi: 10.3389/fvets.2025.1682443

Received: 08 August 2025; Accepted: 04 November 2025;
Published: 20 November 2025.

Edited by:

Livia D'Angelo, University of Naples Federico II, Italy

Reviewed by:

Timm Konold, Animal and Plant Health Agency, United Kingdom
Marco Zaccaroni, University of Florence, Italy

Copyright © 2025 Skoták, Cukor, Sedláček, Ševčík, Linda, Brynychová, Kostka, Macho-Maschler, Hušek and Palme. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jan Cukor, Y3Vrb3JAZmxkLmN6dS5jeg==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.