Next Article in Journal
Real-Time Temperature Distribution Monitoring in Chinese Solar Greenhouse Using Virtual LAN
Next Article in Special Issue
Relationship between Temperate Grass Sward Characteristics and the Grazing Behavior of Dairy Heifers
Previous Article in Journal
A Non-Destructive Method for Grape Ripeness Estimation Using Intervals’ Numbers (INs) Techniques
Previous Article in Special Issue
Variation of Fatty Acids in Cool-Season Grasses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hair Cortisol as a Measure of Chronic Stress in Ewes Grazing Either Hardwood Silvopastures or Open Pastures

School of Plant and Environmental Sciences, College of Agriculture and Life Sciences, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(7), 1566; https://doi.org/10.3390/agronomy12071566
Submission received: 19 May 2022 / Revised: 13 June 2022 / Accepted: 27 June 2022 / Published: 29 June 2022

Abstract

:
Hair cortisol is a relatively non-invasive and reliable measure of chronic stress, but it has received limited use, especially in pasture systems. A two-year study was carried out to compare behavioral and physiological (intravaginal temperature, hair, and blood cortisol) responses of ewes (Ovis aries) that grazed black walnut (Juglans nigra) silvopasture (BSP), honeylocust (Gleditsia triacanthos) silvopasture (HSP), or open pastures (OP) treatments. Ewe weights and intravaginal temperatures were recorded once for every 3-week interval. Plasma and hair cortisol concentrations were determined by ELISA. Trail cameras detected animal behavior. Ewe average daily gain was greater in HSP compared with OP (p = 0.0456) but did not differ with BSP (p = 0.4686) across both years. Ewes on OP had higher (p < 0.0001) hair cortisol concentrations than ewes on silvopasture treatments both summers. Ewes on OP had ≥ 0.4 °C higher (p ≤ 0.03) intravaginal temperatures during portions of the afternoon than ewes managed in silvopasture treatments. Ewes on OP spent 500–700% more (p < 0.0001) time standing and 125–150% less (p < 0.0001) time lying down compared with ewes on silvopasture treatments. Hair cortisol measures could be an effective and relatively non-invasive technique for determining long-term chronic stress in grazing animals.

1. Introduction

Animal welfare has been a major area of concern for the livestock industry [1]. Heat stress in extensively managed grazing animals can compromise animal welfare and productive performance, causing substantial economic consequences to producers. Silvopasture, which involves the deliberate integration of trees, livestock, and forages into a single management area [2], can be a sustainable means of reducing livestock heat stress in grazing systems. Silvopastures provide shading opportunities that can alter the micro-climatic conditions thus creating a favorable environment for grazing livestock [3,4,5]. This can help reduce the heat stress level in animals ultimately optimizing their overall health and well-being.
Quantifying physiological benefits for animals in pastoral systems such as silvopasture has been challenging as there are various external forces influencing animals in the system. Although weight gain typically is the default measure of animal productivity, it may not fully assess the effects of acute and chronic exposure to high-temperature environments. Cortisol, a primary stress hormone in ruminants [6], can be an important indicator of stress due to acute and chronic exposure to high temperatures, which can significantly elevate cortisol levels in animals [7]. Blood is a common matrix for assessing cortisol levels in animals, but the sampling procedure requires capturing and restraining animals which itself increases the cortisol level [8,9] potentially confounding the reliability of the assessment. Hair cortisol can be a reliable and relatively less invasive method of accessing long-term chronic stress levels in animals and this measure is not confounded by activities such as handling/restraining [10,11]. However, it has received limited use, especially in grazing systems. Similarly, the rectal temperature has also been commonly used as means of assessing the internal core body temperature of animals. However, there are limitations to its use in extensively managed animals especially for extended periods mainly due to regular fecal matter flow [5]. The use of vaginal temperature sensors could be a reliable and relatively non-invasive technique for assessing the core body temperature of extensively managed animals for an extended period [12].
Heat stress can significantly alter animal behavior and information on animal behavior can help producers understand the responses to heat stress and help them manage them accordingly. Previously, animal behavior and activity were evaluated through direct observation which is time-intensive and can potentially influence the natural behavior and activity pattern of animals [13]. Time-lapse remote cameras can be used to determine the behavior and activity pattern of animals with minimal disturbance [14]. These methods of assessing physiological and behavioral responses of animals in response to heat stress can be advantageous to traditional methods as they are relatively non-invasive to animals and thus may improve the reliability of the assessment. However, only a few studies have used these methods as a means of assessing heat stress in animals in extensive management systems [5,15,16].
Animal welfare is of increasing concern both for producers and consumers, and the adoption of silvopasture systems may help producers address this issue. The integrative, intensive, interactive, and intentional approach of silvopasture also offers a unique opportunity for landowners to manage their land for both short-term (livestock) and long-term (trees/tree products) economic returns along with a suite of different ecosystem services. Regardless of the opportunities and economic and environmental benefits, the adoption of these systems has been slow in the past [17]. Limited knowledge of these systems, site-specific management requirements, and their economic potential have been constraints to the adoption of silvopasture practices [18]. The understanding of the dynamics between shade and animal behavior/productivity is limited. Silvopastures integrating pine trees have been studied to some extent in the southeastern U.S. [3,4,19,20,21,22,23] and more limited work has been conducted with hardwood systems [5,15,16,24], however, none of these studies have measured animal stress response. The objective of this study was to compare the behavioral and physiological responses of ewes grazing open pastures and black walnut- (Juglans nigra) and honeylocust- (Gleditsia triacanthos) based silvopastures using relatively non-invasive techniques such as hair cortisol. We hypothesized that ewes in silvopastures would have greater average daily gain, lower hair cortisol levels, and cooler intravaginal temperatures when compared with ewes managed in open pastures.

2. Materials and Methods

2.1. Study Site and Its Management

Six-week grazing trials were carried out during the summers of 2020 (late July to early September) and 2021 (early July to mid-August). The studies were carried out at the Whitethorne Agroforestry Demonstration Center at Virginia Tech’s Kentland Farm in Blacksburg, VA (37°12′00.6′′ N 80°34′34.8′′ W). Soil series at the site include Berks-Lowell-Rayne complex, Unison and Braddock, and Weaver, with slopes of 25–65%, 15–25%, and 0–5%, respectively. The study site consists of three experimental pasture systems—open pasture (OP), black walnut silvopastures (BSP), and honeylocust silvopastures (HSP). Trees in the silvopastures were established in 1995 by planting into an existing cool-season pasture. Trees were thinned to an approximate 12.3 m × 12.3 m configuration in 2012. Dominant forage species in both pasture systems include tall fescue (Schedonorus phoenix), orchardgrass (Dactylis glomerata), Kentucky bluegrass (Poa pratensis), and red and white clovers (Trifolium pratense and T. repens).
Pastures were mowed to an approximate height of 15 cm with a rotary mower to remove seed heads and were fertilized at a rate of 56-25-50 kg nitrogen-phosphorus-potassium (NPK) per hectare during late spring in 2020 and 56-0-0 kg NPK per hectare during late spring in 2021. For weed management, ProClova (florpyrauxifen + 2,4-D, Corteva Agriscience, Wilmington, DE) was applied at 4.3 L ha−1 in 2020 before summer grazing to control stickweed (Verbesina occidentalis) and other broadleaf weeds, including creeping thistle (Cirsium arvense) and milk thistle (Silybum marianum) primarily to BSP systems which had the greater presence of these weeds. Trees in the silvopastures have been trimmed periodically to remove side branches (from ~2.5 to 5 m height) to improve access for farm equipment, allow greater light to the understory, and maintain clear boles. Tree density and tree basal area averaged 91 stems ha−1 and 6.1 m2 ha−1 in BSP and 104 stems ha−1 and 4.2 m2 ha−1 in HSP.

2.2. Weather Data

Daily mean ambient temperatures (AT), maximum and minimum temperatures, relative humidity (RH), and rainfall data for the research site were downloaded from Virginia Tech WeatherSTEM Data Mining Tool (http://vt-arec.weatherstem.com) for the entire study period for both years (accessed on 21 December 2021). These data were recorded at a weather station located about 500 m from the study site. Average hourly AT and RH for specific dates when intravaginal temperature and behavior data were recorded were downloaded for both years and used to calculate the average hourly Temperature Humidity Index (THI) using the equation developed by Mader et al. [25].
THI = [(0.8 × AT) + (RH/100) × (AT − 14.4)] + 46.4

2.3. Microclimatic Data

One Spectrum WatchDog 1000 Series Microstation (Spectrum Technologies, Inc., Aurora, IL, USA) was installed in each experimental unit within a replication to collect AT, RH, and photosynthetically active radiation (PAR) data for both summers. These micro stations were preprogrammed to collect measures every 10 min throughout the study period. Data from the loggers were downloaded into a computer using SpecWare 9 Pro software (Spectrum Technologies, Inc., Aurora, IL). To estimate the level of heat stress within the systems, THI was calculated from AT and RH [25].

2.4. Experimental Design, Sheep, and Grazing Management

The experimental systems (OP, BSP, and HSP) were replicated three times and arranged in a randomized complete block design. Thirty-six (36) open Katahdin ewes from Virginia Tech’s Southwest Virginia Agricultural Research and Extension Center, Glade Spring, VA were used for the study each year. In 2020, ewes used for the study were 3 to 6 years old with a mean initial body weight of 53 kg ± 0.9 kg. In 2021, ewes used for the study were 4 to 7 years old with a mean initial body weight of 52 kg ± 1.1 kg. These animals are part of the University’s breeding flock and are raised on pasture and supplemented with grain only during the late gestation/lactation (March to early May). Prior to the study, ewes were vaccinated for Clostridium perfringens type C and D and C. tetani and were dewormed with Cydectin for each year. Ewes were checked for anemia levels by using a FAMACHA protocol [26] prior to the study and at 3-week intervals thereafter. Ewes were stratified by body weight, age, and coat color and four ewes then were randomly assigned to one of the nine experimental units (EU). Each 0.27-ha EU was subdivided into four permanent subpaddocks. Each subpaddock was further subdivided using electrified net fencing to create four approximately equal sections within each subpaddock and ewes were rotational stocked when the forage residual height reached about 7 cm.

2.5. Forage Analysis

Pre-and post-grazing forage biomass within each EU was estimated with a double-sample technique using a rising plate meter (Jenquip, Fielding, NZ) and quadrats [27,28]. Two forage grab samples were also collected randomly by walking across each experimental treatment at the beginning of each season and dried in an oven at 60 °C for 72 h. Samples were ground in a Wiley mill (Thomas Wiley, Philadelphia, PA) to pass through a 2 mm screen and then ground to pass through a 1 mm screen in a cyclone mill (Udy Mill, (UDY Corporation, Fort Collins, CO, USA). Ground samples were scanned in NIRS DS2500F using ISIScan Nova v. 8.0.6.2 (Foss North America, Eden Prairie, MN) to estimate crude protein (CP), neutral detergent fiber (NDF), and acid detergent fiber (ADF) content. Equations for the forage nutritive analysis were standardized and checked for accuracy using the grass hay equation developed by the NIRS Forage and Feed Consortium [29].

2.6. Animal Weight Gain

Animal body weight (BW) of all ewes was measured on two consecutive days at the beginning (days −1 and 0) and end (days 41 and 42) of the study and averaged to determine ewes’ beginning and ending BW. The average daily gain (ADG) of ewes was determined by dividing the total BW gain by the number of days animals were stocked on the treatment pastures.

2.7. Intra-Vaginal Temperatures

Ewe intra-vaginal temperatures were measured with a Star Oddi Data Storage Tag (DST) micro-T temperature logger (Star Oddi, Iceland; sensitivity- 0.2 °C). These loggers were attached to a blank controlled internal drug release (CIDR) device (Eazi-Breed, Zoetis, Parsipanny, NJ, USA) by removing about a 2 cm segment from the middle of the CIDRs, replacing it with the logger, and wrapping the assembled pieces with vinyl electrical tape. CIDRs with loggers were inserted into the vagina of two ewes within each EU for two consecutive days at mid- (days 19 and 20) and end (days 40 and 41) of the study using a CIDR sheep applicator for both years. The loggers were set to collect temperature data at every 10 min interval. Data collected were downloaded through a communication box and Mercury software (Star Oddi, Iceland) and was exported to Microsoft Excel for further analysis.

2.8. Hair and Plasma Collection, Cortisol Extraction, and Analysis

On day 0 of the experiment, hair samples were collected from an approximate 15 cm × 15 cm site in the loin region of all ewes by clipping close to the skin with an electric clipper (900cl Cordless Clipper with Eagle 30 Small Clipper Blade Set, Premier 1 Supplies, Washington, IA, USA). This sample served as a baseline measure of hair cortisol. The same site was trimmed again on days 21 and 42. Hair samples were wrapped in aluminum foil and stored at room temperature until analysis. Cortisol in hair samples was extracted using methanol [11]. Briefly, 250 mg of hair was washed 3 times with a 5-mL aliquot of 100% isopropanol and was left to dry at room temperature for 72 h. Dried hair samples were ground with a mini bead beater (BioSpec Products, Inc., Bartlesville, OK, USA) for 5 min at 30 Hz. Cortisol was extracted from ground hair samples by adding 1 mL of 100% methanol in a microcentrifuge tube with 50 mg of ground hair and rotating in an orbital shaker for ~24 h. Samples were centrifuged and 600 μL of liquid supernatant was aliquoted to a new tube. Methanol was dried off and samples were reconstituted with ELISA buffer. Cortisol concentration was quantified with a commercial salivary cortisol ELISA (Enzyme-Linked Immunosorbent Assay) kit (Cayman Chemical, MI, USA) according to the manufacturer’s instructions. Along with the hair samples, blood samples were also collected from all ewes by jugular venipuncture on days 0, 21, and 42. Samples were kept in a cooler with ice until plasma was separated from blood samples by centrifugation at 3400× g at room temperature for 15 min. The plasma samples collected were stored at −70 °C until further analysis. Cortisol was extracted by adding 5× volume of ethyl ether into the plasma sample, drying the ethyl ether overnight, and reconstituting with ELISA buffer [30]. The plasma cortisol level was determined using the same commercial cortisol ELISA test kit (Cayman Chemical, MI, USA).

2.9. Animal Behavior Data

Time-lapse imagery was collected simultaneously with measures of intravaginal temperature (two consecutive days every 3 weeks) using Moultrie D-500 trail cameras (EBSCO Industries, Inc., Birmingham, AL). One camera was set up in each EU in such a way that it could capture images of the entire section of the sub paddock where ewes were stocked. Two (2) ewes in each EU with temperature loggers were sprayed with orange or blue color fluorescent paint (one color per ewe) for easy visual identification. Images were captured from morning (0700 or 0800) to evening (2100 h) at one-minute intervals. Animal behaviors (grazing, lying, standing up, drinking water, and eating salt) were recorded manually for each captured image, and the total time spent in each activity by each ewe was summed.

2.10. Statistical Analysis

A regression equation of plate meter reading against forage mass was calculated with PROC REG in SAS Studio, v. 3.5 (SAS Inst., Cary, NC, USA) and the best fit was determined with a quadratic equation. A mixed-effect analysis of variance test was carried out to determine the difference in ADG, forage nutritive value, pre-and post-grazing forage biomass, cortisol measures, intravaginal temperature, and animal behavior using PROC MIXED in SAS Studio, v. 3.5. The study was conducted as a randomized complete block design with three replications; year was included as a random effect. Repeated measures analysis by period was used with a standard variance-covariance structure for the analysis of hair cortisol and intravaginal temperatures data. A compound symmetry variance-covariance structure was used for the analysis of plasma cortisol and ewes’ behavior data. Repeated measures analysis by sampling data was used with a compound symmetry variance-covariance structure for the analysis of pre-and post-grazing forage biomass and herbage disappearance data. Variance-covariance structure for the analysis was selected based on the lowest AIC value. Microclimatic data for both summers were analyzed using a one-way analysis of variance test using PROC GLM in SAS Studio, v.3.5. LS- means and Tukey’s adjusted differences were calculated. Differences were considered significant when p < 0.05 and was reported as trends when 0.05 < p < 0.10.

3. Results

3.1. Weather Data

The mean daily AT throughout the study period in 2020 and 2021 was 22.6 °C and 25.2 °C, respectively. Mean minimum and maximum AT were 18.3 °C and 29.0 °C in 2020 and 17.4 °C and 29.7 °C, respectively, in 2021 (Table 1). In 2020 and 2021, the study site received 7.7 and 1.8 cm of total precipitation throughout the respective study periods. The temperature-humidity index at the study site was above 72 during the late morning and afternoon hours (1000 to 1900 h) during the intravaginal temperature and behavior data collection dates for both years (Figure 1).

3.2. Microclimatic Conditions

In 2020, average daily AT in OP tended (p = 0.0599) to be greater than in BSP whereas it was similar compared with HSP (p = 0.8768; Table 2). In 2021, the average daily AT did not differ among treatments (p ≥ 0.4943). In 2020, RH did not differ among treatments, whereas in 2021, RH was higher in OP than in silvopasture treatments (p ≤ 0.0375). The PAR was higher in OP compared with silvopasture treatments for both years (p < 0.0001).

3.3. Forage Measures

3.3.1. Forage Biomass

The year x treatment interaction was significant for pre-grazing forage biomass and herbage disappearance (p < 0.0001) and non-significant for post-grazing forage biomass (p = 0.1634). Pre-grazing forage biomass in OP was greater (p < 0.0001) compared with BSP for both years whereas it was greater (p < 0.0001) in OP compared with HSP in 2020, but not in 2021 (p = 0.4373; Table 3). Post-grazing forage biomass was greater (p < 0.0001) in OP compared with silvopasture treatments across both years and it was greater (p = 0.0031) in HSP compared with BSP treatment across both years. In 2020, forage disappearance did not differ (p ≥ 0.3867) among treatments whereas in 2021, HSP had greater (p < 0.0001) herbage disappearance followed by BSP and OP treatments.

3.3.2. Forage Nutritive Value and Digestibility

There was no year x treatment interaction for forage CP (p = 0.5887), ADF (p = 0.7700), and NDF (p = 0.2217) data. Forage CP concentration did not differ (p ≥ 0.7346) between treatments across both years (Table 4). Forage from BSP had greater (p = 0.0342) ADF and tended (p = 0.0979) to have greater NDF than OP across both years.

3.4. Animal Gain

There was no year x treatment interaction for ADG (p = 0.1781). The overall ADG of ewes in OP, BSP, and HSP was 51.3, 73.6, and 85.9 g d−1, respectively. Greater ADG was measured in HSP (p = 0.0457) compared with OP across both years (Table 5). The ADG of ewes in BSP did not differ from that of ewes on HSP (p = 0.4686) and OP (p = 0.1902) across both years.

3.5. Intravaginal Temperature

There was no year x treatment interaction (p ≥ 0.1089) for intravaginal temperature data. At 1000 h, the intravaginal temperature of ewes on OP tended to be higher (p = 0.0730) compared with ewes on silvopasture treatments (Figure 2). Ewes on OP had higher (p ≤ 0.0343) intravaginal temperatures than ewes on both silvopasture treatments between 1100 h–1700 h and ewes on HSP between 1800 h–2000 h. The intravaginal temperature of ewes on OP tended to be higher than ewes on HSP at 2100 h (p = 0.0980) and ewes on BSP at 1800 h (p = 0.0804).

3.6. Plasma and Hair Cortisol Concentration

There was no year x treatment interaction for plasma cortisol data (p = 0.3613) whereas there was a significant year x treatment interaction for hair cortisol data (p = 0.0356). Plasma cortisol level was greater (p = 0.0400) in ewes on OP compared with BSP whereas there was no difference in plasma cortisol level in ewes on HSP and BSP (p = 0.6954). There was a trend (p = 0.0983) toward greater plasma cortisol level of ewes on OP compared with HSP treatment across both years (Table 6). Ewes on OP had a greater hair cortisol level compared with ewes on silvopasture treatments both in 2020 and 2021 (p < 0.0001). In 2021, ewes on BSP had lower (p = 0.0236) hair cortisol levels compared with ewes on HSP.

3.7. Animal Behavior

There was a trend toward significant year x treatment interaction for time ewes spent grazing (p = 0.0611) and a significant year x treatment interaction for time ewes spent eating salt (p < 0.0001). There was no year x treatment interaction (p ≥ 0.2747) for standing, lying, and drinking water behavior categories. In 2020, ewes on BSP spent about 20% more time grazing than ewes on OP (p = 0.0054) and HSP (p = 0.0106) whereas in 2021 ewes on BSP spent about 36% more time (p = 0.0014) grazing than ewes on HSP and ewes on OP spent about 20% more time (p = 0.0192) grazing compared with ewes on HSP (Table 7). Ewes on OP spent 400% more time standing than ewes on BSP (p < 0.0001) and 750% more time standing than ewes on HSP (p < 0.0001). Ewes on OP spent 20% less time lying down compared with ewes on BSP (p < 0.0001) and 33% less time lying down compared with ewes on HSP (p < 0.0001). Ewes on silvopasture treatments showed some trend toward 70% less time (p ≤ 0.0599) drinking water compared with ewes on OP.

4. Discussion

4.1. Microclimatic Condition

Microclimatic measures (ambient temperature, relative humidity, soil moisture, and PAR) indicate environmental conditions in OP were more stressful than in silvopasture treatments. Ambient temperature and PAR tend to be lower under tree canopies compared with open systems [31], and past studies have reported cooler microclimatic conditions in silvopasture versus OP [3,4,32,33]. The BSP and HSP silvopasture systems used in this study provided shade and lowered the PAR and ambient temperature within the system compared with OPs with no canopy cover. Differences in tree species also were observed, as PAR was greater in HSP compared with BSP. This is likely due to the differences in tree and leaf size and structure, and differences in canopy architecture between honeylocust and black walnut trees. Although we did not make direct measures of tree, leaf, and canopy size in BSP and HSP treatments, past measures and visual observation indicated black walnut trees generally were larger, with larger canopies, larger leaves and leaflet and denser canopies. These differences translated to lower light levels within the BSP system.

4.2. Forage Measures

4.2.1. Forage Biomass

The impact of trees on forage productivity in silvopasture is dependent on various factors such as forage and tree species utilized, age of the system, and tree arrangement [34]. Pre-grazing forage biomass was 20% less in BSP compared with OP. This corresponded with reduced levels of PAR in that system but may also reflect changes in forage species that have been observed under BSP at this site [5,15,16,35]. Changes in forage biomass in response to increasing shade levels are not linear [36] but yield generally declines substantially when PAR is reduced to less than 50% of full sun. Similar yield reductions have been observed in pine-walnut silvopastures in Missouri when compared with OP systems [23]. In contrast, the smaller honeylocust trees, with thinner canopy cover, allowed more PAR to penetrate to the forage understory; thus, forage production in HSP was similar to that in OP. There has been no significant reduction in forage productivity of cool-season grasses with up to 50% shade in some studies [37,38]. These results underscore the importance of light and tree canopy management for maintaining forage productivity under silvopasture systems.
Greater forage disappearance in silvopasture treatments compared with OP might be attributed to greater intake by ewes in silvopasture treatments. However, actual forage intake by ewes was not measured in this study, and changes in the rate of disappearance might also be due to the underestimation of residual forage mass in silvopasture treatments. Greater intake by lambs in BSP compared with OP (as estimated from herbage disappearance) has been observed at this research site [39]. Intake estimates of animals in silvopasture systems have received very limited study and require more investigation to better understand the efficiency of these systems.

4.2.2. Forage Nutritive Value and Digestibility

The trees within silvopasture systems can significantly impact understory forage nutritive value through botanical, morphological, and physiological adaptations [37,40]. Past studies on both cool-season [41] and warm-season [42] forages have reported improvement to forage nutritive value and digestibility underneath trees or shade as compared with open systems, although equal or lower digestibility of forages underneath the shade has also been reported [43,44]. Forage CP content did not differ among treatments whereas forage ADF and NDF content were greater in BSP compared with OP across both years in the present study. A study by Fannon et al. [35] also reported no difference in forage CP and greater forage ADF, NDF, and lignin content in BSP compared with OP during summer at this site. Higher fiber concentration in BSP likely reflects the differences in vegetation composition between treatments as the competition of warm-season grass (nimblewill’ Muhlenbergia scherberi) was greater in BSP compared with OP [35].

4.2.3. Animal Gain

Animal productivity in grazing systems is directly influenced by a number of variables. Ambient temperature, relative humidity, and solar radiation are the key climatic variables that affect animal production [45] and these climatic conditions can severely impact the overall health, metabolism, and physiology of animals [46,47]. The impacts of climatic variables may be lessened or eliminated by integrating trees within the system to promote animals’ comfort and help improve animal productivity [48,49]. In some cases, the reduction in stress associated with grazing under trees could compensate for the reduction in available forage. The average daily gain of ewes in OP was lower compared with ewes in HSP but did not differ with BSP treatment. This is likely due to greater intake by ewes in HSP as the herbage disappearance in HSP was higher than OP. Additionally, a cooler ambient condition in HSP due to tree shade may have reduced heat stress in ewes thus improving their overall productivity. The animal gain was low during the first summer, but ADG improved significantly in the second year. Forage productivity in silvopasture treatments during the second year was about 10% more along with greater CP and lower fiber content of forages compared with the first year which may have created a better opportunity for ewes gain in the second year compared with the first year. The ewes used for the study were completely new to the environment and may have been stressed due to herd separation which may have significantly impacted their behavior and diet selection during the first year of the study. The literature suggests that sheep when transferred to a new environment significantly modify their behavior and diet selection [50,51] thus influencing their overall performance. However, the same ewes were used during the second year of the study, and they were more familiar with the environment, which may have also helped improve their overall performance.

4.2.4. Intravaginal Temperatures

The intravaginal temperature loggers used to measure ewe core body temperatures in the present study can collect data for an extended period with little disturbance to the animals [12]. Intravaginal temperatures are well correlated with rectal temperatures and other measures of core body temperatures [12]. Ewes on silvopasture treatments were cooler than ewes on OP, especially during periods of greater ambient temperature in the OP systems. A study by Pent et al. [5] also reported a 0.3–0.5 °C reduction in the intravaginal temperature of lambs in BSP compared with lambs in OP during the afternoon hours (1200–1900 h). This difference in core body temperature of ewes between treatments especially during the afternoon hours indicates a stronger response to the ambient condition by ewes in OP than for those in silvopasture treatments. These findings support the work of others [5,52] who found greater core body temperatures in animals without shade during the afternoon hours, corresponding to the hottest part of the day. Shade within the silvopasture systems helped lower the core body temperature of ewes thus reducing their heat load and overall stress compared with ewes on OP.

4.2.5. Plasma and Hair Cortisol Concentration

Measures of plasma cortisol did not differ between ewes that grazed silvopasture or OP treatments, and the data were highly variable. The sampling procedure for this assay required that animals be moved from treatment pastures to a central pen and subsequently handled and restrained to draw blood. Such procedures can cause acute stress, thus this procedure itself likely increased the blood cortisol levels in our animals [52] and would potentially confound the reliability of the assessment. This may have resulted in high variability of data in an acute measure and thus a nonsignificant difference in plasma cortisol levels among treatments in the current study.
Hair cortisol levels were lower in ewes on silvopasture treatments compared with ewes on OP. Blood cortisol diffuses into the hair and accumulates over time as animals are exposed to stressors such as elevated temperature. Hair cortisol thus reflects long-term chronic stress levels in animals over weeks to months [53] and is not confounded by activities such as animal handling and restraint [10,11], and various studies have reported a positive correlation between cortisol level in hair and stress [54,55]. Thus, hair cortisol could provide a more reliable means of assessing long-term chronic heat stress of animals in extensive grazing systems. This approach also is relatively less invasive, collecting hair samples is simple, and the samples can be stored at room temperature for long periods [53]. Lower hair cortisol in ewes on silvopasture treatments likely reflects milder microclimatic conditions in silvopasture treatments compared with OP and the lower body temperatures maintained through the hottest part of the day. Shade minimizes heat stress in ewes [56] and, in past work, sheep on this site were seen to actively move to access shade during the day [15].

4.2.6. Animal Behavior

Ewes on BSP spent greater time grazing compared with ewes of HSP and OP treatments Likely due to more limited forage availability. Animals spend less time grazing when there is plenty of forage available and is of good quality whereas animals spend more time grazing when the quantity or quality of forages is limited [57]. Ewes in silvopastures experienced less heat stress and spent more time lying down than ewes in OP. In contrast, ewes on OP were more stressed and thus spent more time standing up. Time spent lying down is considered a traditional metric of animal comfort [58,59] whereas standing up or loafing is considered a behavioral response to heat stress [60,61]. Standing improves airflow and increases the effectiveness of convection heat loss from the body in animals under heat stress, thus ewes on OP may have spent more time standing as a strategy to reduce head load [62]. More time spent standing could also be an adaptive mechanism to limit heat gain from the ground surface. Higher solar radiation and ambient temperature may have increased the surface temperature in OP compared with silvopasture treatments, limiting ewes’ ability to lie down and lose heat through conduction. However, in silvopasture treatments, cooler ambient conditions due to shade may also have reduced surface soil temperatures, thus encouraging ewes to lie down and lose heat through conduction to the soil. Ewes on OP also spent more time drinking water compared with ewes on silvopasture treatments. This might also be a response of ewes to a stressful environment, as drinking water can help reduce heat load in animals [63]. The differences observed in ewes’ behavior among treatments likely reflect both direct effects of altered microclimatic conditions (and greater heat stress in OP systems) as well as indirect effects in terms of changes in forage composition and productivity (which affected grazing behavior, particularly in BSP).

5. Conclusions

Animal welfare is a major concern for the livestock industry. Heat stress in grazing systems can compromise animal welfare and productive performance, causing economic losses to livestock producers. Silvopastures provide shading opportunities that can help reduce heat stress and improve the health and well-being of grazing animals. Black walnut and honeylocust trees in these silvopasture systems altered micro-climatic conditions and created a more favorable environment for sheep. Ewes grazing HSP had greater gains compared with those in OP. Ewes in silvopasture systems had cooler intravaginal temperatures and lower hair cortisol level concentrations, indicating lower stress levels. Because blood cortisol gets elevated quickly in response to acute stressors such as handling and restraint, hair cortisol provided a more reliable and relatively non-invasive method of assessing animal stress in these extensive grazing systems. This approach may be valuable where measures of chronic (vs. acute) stress are appropriate for the given research. Although both tree species used in the study moderated ambient conditions and helped reduce stress levels in ewes, tree species had variable effects on animal physiology and overall productivity of the system, which should be a key consideration for producers adopting silvopasture systems.

Author Contributions

Conceptualization, S.P. and J.H.F.; methodology, S.P. and J.H.F.; formal analysis, S.P., J.H.F. and G.J.P.; investigation, S.P. and J.H.F.; resources, J.H.F. and G.J.P.; writing—review and editing, S.P., J.H.F. and G.J.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by USDA-NIFA Southern SARE Graduate Student Grant Number GS19-202.

Institutional Review Board Statement

The animal study protocol was approved by the Institutional Review Board of the Virginia Tech IACUC Office under protocol # 20-017.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank crew members at SWAREC for providing animals for the study; the Kentland farm crew for helping with site preparation and maintenance; Sierra Guynn and her team for assisting with veterinary care; Swarup Podder, Aleks Halili, and Eric McKenzie for assisting with the fieldwork. The funding for the study was provided by USDA-NIFA Southern SARE Graduate Student Grant Number GS19-202. The authors would like to acknowledge Virginia Tech Libraries for providing Open Access Subvention Fund (OASF) to support article processing charge.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alonso, M.E.; González-Montaña, J.R.; Lomillos, J.M. Consumers’ concerns and perceptions of farm animal welfare. Animals 2020, 10, 385. [Google Scholar] [CrossRef]
  2. Pent, G.J.; Fike, J.; Orefice, J.N.; Sharrow, S.H.; Brauer, D.; Clason, T.R. Silvopasture practices. In North American Agroforestry, 3rd ed.; Garrett, H.E., Jose, H., Gold, H.A., Eds.; American Society of Agronomy: Madison, WI, USA, 2021; pp. 127–162. ISBN 9780891183778. [Google Scholar]
  3. Karki, U.; Goodman, M.S. Cattle distribution and behavior in southern-pine silvopasture versus open-pasture. Agrofor. Syst. 2010, 78, 159–168. [Google Scholar] [CrossRef]
  4. Karki, U.; Goodman, M.S. Microclimatic differences between mature loblolly-pine silvopasture and open-pasture. Agrofor. Syst. 2015, 89, 319–325. [Google Scholar] [CrossRef]
  5. Pent, G.J.; Fike, J.H.; Kim, I. Ewe lamb vaginal temperatures in hardwood silvopastures. Agrofor. Syst. 2021, 95, 21–32. [Google Scholar] [CrossRef]
  6. Binsiya, T.; Sejian, V.; Bagath, M.; Krishnan, G.; Hyder, I.; Manimaran, A.; Lees, A.; Gaughan, J.; Bhatta, R. Significance of hypothalamic-pituitary-adrenal axis to adapt to climate change in livestock. Int. Res. J. Agric. Food Sci. 2017, 2, 1–20. [Google Scholar]
  7. Möstl, E.; Palme, R. Hormones as indicators of stress. Domest. Anim. Endocrinol. 2002, 23, 67–74. [Google Scholar] [CrossRef]
  8. Widmaier, E.P.; Kunz, T.H. Basal, diurnal, and stress-induced levels of glucose and glucocorticoids in captive bats. J. Exp. Zool. 1993, 265, 533–540. [Google Scholar] [CrossRef]
  9. Müller, C.; Jenni-Eiermann, S.; Blondel, J.; Perret, P.; Caro, S.P.; Lambrechts, M.; Jenni, L. Effect of human presence and handling on circulating corticosterone levels in breeding blue tits (Parus caeruleus). Gen. Comp. Endocrinol. 2006, 148, 163–171. [Google Scholar] [CrossRef]
  10. Goymann, W.; Möstl, E.; Van’t Hof, T.; East, M.L.; Hofer, H. Noninvasive fecal monitoring of glucocorticoids in spotted hyenas, Crocuta crocuta. Gen. Comp. Endocrinol. 1999, 114, 340–348. [Google Scholar] [CrossRef]
  11. Nejad, J.G.; Lohakare, J.; Son, J.; Kwon, E.; West, J.; Sung, K. Wool cortisol is a better indicator of stress than blood cortisol in ewes exposed to heat stress and water restriction. Animal 2014, 8, 128–132. [Google Scholar] [CrossRef]
  12. Burdick, N.; Carroll, J.; Dailey, J.; Randel, R.; Falkenberg, S.; Schmidt, T. Development of a self-contained, indwelling vaginal temperature probe for use in cattle research. J. Therm. Biol. 2012, 37, 339–343. [Google Scholar] [CrossRef]
  13. Bridges, A.S.; Noss, A.J. Behavior and activity patterns. In Camera Traps in Animal Ecology; Springer: Berlin/Heidelberg, Germany, 2011; pp. 57–69. [Google Scholar]
  14. Herlin, A.; Brunberg, E.; Hultgren, J.; Högberg, N.; Rydberg, A.; Skarin, A. Animal welfare implications of digital tools for monitoring and management of cattle and sheep on pasture. Animals 2021, 11, 829. [Google Scholar] [CrossRef]
  15. Pent, G.J.; Greiner, S.P.; Munsell, J.F.; Tracy, B.F.; Fike, J.H. Lamb performance in hardwood silvopastures, II: Animal behavior in summer. Transl. Anim. Sci. 2020, 4, 363–375. [Google Scholar] [CrossRef]
  16. Pent, G.J.; Greiner, S.P.; Munsell, J.F.; Tracy, B.F.; Fike, J.H. Lamb performance in hardwood silvopastures, I: Animal gains and forage measures in summer. Transl. Anim. Sci. 2020, 4, 385–399. [Google Scholar] [CrossRef]
  17. Frey, G.E.; Fike, J.H. Silvopasture case studies in North Carolina and Virginia. In e-General Technical Report SRS–236; US Department of Agriculture Forest Service, Southern Research Station: Asheville, NC, USA, 2018; Volume 236, pp. 1–23. [Google Scholar]
  18. Smith, M.M.; Bentrup, G.; Kellerman, T.; MacFarland, K.; Straight, R.; Ameyaw, L.; Stein, S. Silvopasture in the USA: A systematic review of natural resource professional and producer-reported benefits, challenges, and management activities. Agric. Ecosyst. Environ. 2022, 326, 107818. [Google Scholar] [CrossRef]
  19. Clason, T.R. Silvopastoral practices sustain timber and forage production in commercial loblolly pine plantations of northwest Louisiana, USA. Agrofor. Syst. 1998, 44, 293–303. [Google Scholar] [CrossRef]
  20. Kallenbach, R.; Kerley, M.; Bishop-Hurley, G. Cumulative forage production, forage quality and livestock performance from an annual ryegrass and cereal rye mixture in a pine walnut silvopasture. Agrofor. Syst. 2006, 66, 43–53. [Google Scholar] [CrossRef]
  21. Karki, U.; Poudel, S.; Bhattrai, S.; Karki, L.; Singh, A. Performance, Carcass Characteristics, and Meat Quality of Kiko Wethers Finished on Pasture and Browse. J. Anim. Sci. 2018, 96, 458–459. [Google Scholar] [CrossRef]
  22. Poudel, S.; Karki, U.; McElhenney, W.; Karki, Y.; Tillman, A. Behavior and distribution patterns of Katahdin rams in southern-pine silvopastures with cool-season forages. Agrofor. Syst. 2019, 93, 1887–1896. [Google Scholar] [CrossRef]
  23. Poudel, S.; Karki, U.; Karki, Y.; Tillman, A. Confinement influenced the diurnal behavior of Katahdin rams and Kiko wethers in southern-pine silvopastures. Agrofor. Syst. 2020, 94, 29–40. [Google Scholar] [CrossRef]
  24. McDaniel, A.; Roark, C. Performance and grazing habits of Hereford and Aberdeen-Angus cows and calves on improved pastures as related to types of shade. J. Anim. Sci. 1956, 15, 59–63. [Google Scholar] [CrossRef]
  25. Mader, T.L.; Davis, M.; Brown-Brandl, T. Environmental factors influencing heat stress in feedlot cattle. J. Anim. Sci. 2006, 84, 712–719. [Google Scholar] [CrossRef]
  26. Kaplan, R.; Burke, J.; Terrill, T.; Miller, J.; Getz, W.; Mobini, S.; Valencia, E.; Williams, M.; Williamson, L.; Larsen, M. Validation of the FAMACHA© eye color chart for detecting clinical anemia in sheep and goats on farms in the southern United States. Vet. Parasitol. 2004, 123, 105–120. [Google Scholar] [CrossRef]
  27. Macoon, B.; Sollenberger, L.; Moore, J.; Staples, C.; Fike, J.; Portier, K. Comparison of three techniques for estimating the forage intake of lactating dairy cows on pasture. J. Anim. Sci. 2003, 81, 2357–2366. [Google Scholar] [CrossRef]
  28. Pent, G.J.; Fike, J.H. Lamb productivity on stockpiled fescue in honeylocust and black walnut silvopastures. Agrofor. Syst. 2019, 93, 113–121. [Google Scholar] [CrossRef]
  29. McIntosh, D.W.; Kern-Lunbery, R.; Goldblatt, P.; Lemus, R.; Griggs, T.; Bauman, L.; Boone, S.; Shewmaker, G.; Teutsch, C. Guidelines for Optimal Use of NIRSC Forage and Feed Calibrations in Membership Laboratories, 2nd ed.; The University of Tennessee Press: Knoxville, TN, USA, 2022. [Google Scholar]
  30. Negrao, J.; Porcionato, M.; De Passille, A.; Rushen, J. Cortisol in saliva and plasma of cattle after ACTH administration and milking. J. Dairy Sci. 2004, 87, 1713–1718. [Google Scholar] [CrossRef]
  31. Morecroft, M.; Taylor, M.E.; Oliver, H.R. Air and soil microclimates of deciduous woodland compared to an open site. Agric. For. Meteorol. 1998, 90, 141–156. [Google Scholar] [CrossRef]
  32. Valigura, R.A.; Messina, M.G. Modification of Texas clear-cut environments with loblolly pine shelterwoods. J. Environ. Manag. 1994, 40, 283–295. [Google Scholar] [CrossRef]
  33. Karki, U.; Goodman, M.S. Microclimatic differences between young longleaf-pine silvopasture and open-pasture. Agrofor. Syst. 2013, 87, 303–310. [Google Scholar] [CrossRef]
  34. Hawke, M. Pasture production and animal performance under pine agroforestry in New Zealand. For. Ecol. Manag. 1991, 45, 109–118. [Google Scholar] [CrossRef]
  35. Fannon, A.; Fike, J.H.; Greiner, S.P.; Feldhake, C.M.; Wahlberg, M. Hair sheep performance in a mid-stage deciduous Appalachian silvopasture. Agrofor. Syst. 2019, 93, 81–93. [Google Scholar] [CrossRef]
  36. Mercier, K.M.; Teutsch, C.D.; Fike, J.H.; Munsell, J.F.; Tracy, B.F.; Strahm, B.D. Impact of increasing shade levels on the dry-matter yield and botanical composition of multispecies forage stands. Grass Forage Sci. 2020, 75, 291–302. [Google Scholar] [CrossRef]
  37. Kephart, K.D.; Buxton, D.R.; Taylor, S.E. Growth of C3 and C4 perennial grasses under reduced irradiance. Crop Sci. 1992, 32, 1033–1038. [Google Scholar] [CrossRef]
  38. Buergler, A.L.; Fike, J.H.; Burger, J.A.; Feldhake, C.R.; McKenna, J.A.; Teutsch, C.D. Botanical composition and forage production in an emulated silvopasture. Agron. J. 2005, 97, 1141–1147. [Google Scholar] [CrossRef]
  39. Pent, G.J. Lamb Performance, Behavior, and Body Temperatures in Hardwood Silvopasture Systems. Ph.D. Dissertation, Virginia Tech, Blacksburg, VA, USA, 2017. Available online: http://hdl.handle.net/10919/76730 (accessed on 21 December 2021).
  40. Kephart, K.D.; Buxton, D.R. Forage quality responses of C3 and C4 perennial grasses to shade. Crop Sci. 1993, 33, 831–837. [Google Scholar] [CrossRef]
  41. Krueger, W.C. How a forest affects a forage crop. Rangelands 1981, 3, 70–71. [Google Scholar]
  42. Eriksen, F.; Whitney, A. Effects of Light Intensity on Growth of Some Tropical Forage Species. I. Interaction of Light Intensity and Nitrogen Fertilization on Six Forage Grasses 1. Agron. J. 1981, 73, 427–433. [Google Scholar] [CrossRef]
  43. Lin, C.; McGraw, M.; George, M.; Garrett, H. Nutritive quality and morphological development under partial shade of some forage species with agroforestry potential. Agrofor. Syst. 2001, 53, 269–281. [Google Scholar] [CrossRef]
  44. Neel, J.; Felton, E.; Singh, S.; Sexstone, A.; Belesky, D. Open pasture, silvopasture and sward herbage maturity effects on nutritive value and fermentation characteristics of cool-season pasture. Grass Forage Sci. 2016, 71, 259–269. [Google Scholar] [CrossRef]
  45. Habeeb, A.A.; Gad, A.E.; Atta, M.A. Temperature-humidity indices as indicators to heat stress of climatic conditions with relation to production and reproduction of farm animals. Int. J. Biotechnol. Recent Adv. 2018, 1, 35–50. [Google Scholar] [CrossRef]
  46. Nardone, A.; Ronchi, B.; Lacetera, N.; Ranieri, M.S.; Bernabucci, U. Effects of climate changes on animal production and sustainability of livestock systems. Livest. Sci. 2010, 130, 57–69. [Google Scholar] [CrossRef]
  47. Nardone, A.; Ronchi, B.; Lacetera, N.; Bernabucci, U. Climatic effects on productive traits in livestock. Vet. Res. Commun. 2006, 30, 75. [Google Scholar] [CrossRef]
  48. Kendall, P.; Nielsen, P.; Webster, J.; Verkerk, G.; Littlejohn, R.; Matthews, L. The effects of providing shade to lactating dairy cows in a temperate climate. Livest. Sci. 2006, 103, 148–157. [Google Scholar] [CrossRef]
  49. Fisher, A.; Roberts, N.; Bluett, S.; Verkerk, G.; Matthews, L. Effects of shade provision on the behaviour, body temperature and milk production of grazing dairy cows during a New Zealand summer. N. Z. J. Agric. Res. 2008, 51, 99–105. [Google Scholar] [CrossRef]
  50. Hunter, R.; Davies, G. The effect of method of rearing on the social behaviour of Scottish Blackface hoggets. Anim. Sci. 1963, 5, 183–194. [Google Scholar] [CrossRef]
  51. Arnold, G. Some principles in the investigation of selective grazing. Proc. Aust. Soc. Anim. Prod. 1964, 5, 258–271. [Google Scholar]
  52. Säkkinen, H.; Tornbeg, J.; Goddard, P.; Eloranta, E.; Ropstad, E.; Saarela, S. The effect of blood sampling method on indicators of physiological stress in reindeer (Rangifer tarandus tarandus). Domest. Anim. Endocrinol. 2004, 26, 87–98. [Google Scholar] [CrossRef]
  53. Gow, R.; Thomson, S.; Rieder, M.; Van Uum, S.; Koren, G. An assessment of cortisol analysis in hair and its clinical applications. Forensic Sci. Int. 2010, 196, 32–37. [Google Scholar] [CrossRef]
  54. Davenport, M.D.; Tiefenbacher, S.; Lutz, C.K.; Novak, M.A.; Meyer, J.S. Analysis of endogenous cortisol concentrations in the hair of rhesus macaques. Gen. Comp. Endocrinol. 2006, 147, 255–261. [Google Scholar] [CrossRef]
  55. Yamada, J.; Stevens, B.; de Silva, N.; Gibbins, S.; Beyene, J.; Taddio, A.; Newman, C.; Koren, G. Hair cortisol as a potential biologic marker of chronic stress in hospitalized neonates. Neonatology 2007, 92, 42–49. [Google Scholar] [CrossRef]
  56. Blackshaw, J.K.; Blackshaw, A. Heat stress in cattle and the effect of shade on production and behaviour: A review. Aust. J. Exp. Agric. 1994, 34, 285–295. [Google Scholar] [CrossRef]
  57. Lyons, R.K.; Machen, R.V. Interpreting Grazing Behavior. Agrilife Communication and Marketing, Texas A&M University System. 2000. Available online: https://hdl.handle.net/1969.1/86955 (accessed on 21 December 2021).
  58. Haley, D.; Rushen, J.; Passillé, A.d. Behavioural indicators of cow comfort: Activity and resting behaviour of dairy cows in two types of housing. Can. J. Anim. Sci. 2000, 80, 257–263. [Google Scholar] [CrossRef]
  59. Winckler, C.; Capdeville, J.; Gebresenbet, G.; Hörning, B.; Roiha, U.; Tosi, M.; Waiblinger, S. Selection of parameters for on-farm welfare-assessment protocols in cattle and buffalo. Anim. Welf. 2003, 12, 619–624. [Google Scholar]
  60. Tucker, C.B.; Rogers, A.R.; Schütz, K.E. Effect of solar radiation on dairy cattle behaviour, use of shade and body temperature in a pasture-based system. Appl. Anim. Behav. Sci. 2008, 109, 141–154. [Google Scholar] [CrossRef]
  61. Cook, N.; Mentink, R.; Bennett, T.; Burgi, K. The effect of heat stress and lameness on time budgets of lactating dairy cows. J. Dairy Sci. 2007, 90, 1674–1682. [Google Scholar] [CrossRef]
  62. Silanikove, N. Effects of heat stress on the welfare of extensively managed domestic ruminants. Livest. Prod. Sci. 2000, 67, 1–18. [Google Scholar] [CrossRef]
  63. De, K.; Kumar, D.; Sharma, S.; Kumawat, P.; Mohapatra, A.; Sahoo, A. Effect of drinking earthen pot water on physiological response and behavior of sheep under heat stress. J. Therm. Biol. 2020, 87, 102476. [Google Scholar] [CrossRef]
Figure 1. Average temperature-humidity index (THI) by hour during the intravaginal temperature and ewes’ behavior data collection dates for the summers of 2020 and 2021.
Figure 1. Average temperature-humidity index (THI) by hour during the intravaginal temperature and ewes’ behavior data collection dates for the summers of 2020 and 2021.
Agronomy 12 01566 g001
Figure 2. Mean vaginal temperatures (SE = 0.1) of Katahdin ewes by hour that grazed either open pasture (OP), black walnut silvopasture (BSP), or honeylocust silvopasture (HSP) treatments. * p < 0.1, ** p < 0.05, *** p < 0.01—Indicates a trend or significant difference between BSP and OP p < 0.1, †† p < 0.05, ††† p < 0.01—Indicates a trend or significant difference between HSP and OP.
Figure 2. Mean vaginal temperatures (SE = 0.1) of Katahdin ewes by hour that grazed either open pasture (OP), black walnut silvopasture (BSP), or honeylocust silvopasture (HSP) treatments. * p < 0.1, ** p < 0.05, *** p < 0.01—Indicates a trend or significant difference between BSP and OP p < 0.1, †† p < 0.05, ††† p < 0.01—Indicates a trend or significant difference between HSP and OP.
Agronomy 12 01566 g002
Table 1. Mean daily temperature (°C), minimum temperature (°C), maximum temperatures (°C), and total precipitation (cm) throughout the study period during the summers of 2020 and 2021 at Whitethorne Agroforestry Demonstration Center, Blacksburg, VA, USA.
Table 1. Mean daily temperature (°C), minimum temperature (°C), maximum temperatures (°C), and total precipitation (cm) throughout the study period during the summers of 2020 and 2021 at Whitethorne Agroforestry Demonstration Center, Blacksburg, VA, USA.
Temperature °CPrecipitation (cm)
YearMax 1MinMean
202029.018.322.67.7
202129.717.425.21.8
1 Max—Maximum; Min—Minimum.
Table 2. Microclimatic conditions in an open pasture (OP), black walnut silvopasture (BSP), and honeylocust silvopasture (HSP) treatments throughout the study period during the summers of 2020 and 2021 at the Whitethorne Agroforestry Demonstration Center, Blacksburg, VA, USA.
Table 2. Microclimatic conditions in an open pasture (OP), black walnut silvopasture (BSP), and honeylocust silvopasture (HSP) treatments throughout the study period during the summers of 2020 and 2021 at the Whitethorne Agroforestry Demonstration Center, Blacksburg, VA, USA.
YearTreatments 1Tukey’s Adjusted p-Value
BSPHSPOPSEBSP vs. HSPBSP vs. OPHSP vs. OP
Temperature, °C
202021.822.622.90.340.21590.05990.8081
202122.723.023.20.270.79450.49430.8768
Relative Humidity, %
202084.486.586.11.260.47580.61280.9733
202177.176.080.71.030.74610.03750.0047
PAR 2, μMolm−2 s−1
202012126536410.6<0.0001<0.0001<0.0001
202118834651711.8<0.0001<0.0001<0.0001
THI 3
202069.670.971.30.540.19710.06870.8653
202171.071.372.00.440.87840.23510.4853
1 Treatments: BSP—Black walnut silvopasture; HSP—Honeylocust silvopasture; OP—Open pasture. 2 PAR—Photosynthetically active radiation; 3 THI—Temperature humidity index.
Table 3. Estimated pre- and post-grazing forage biomass and herbage disappearance (difference) in an open pasture (OP), black walnut silvopasture (BSP), and honeylocust silvopasture (HSP) treatments during the summers of 2020 and 2021 at the Whitethorne Agroforestry Demonstration Center, Blacksburg, VA, USA.
Table 3. Estimated pre- and post-grazing forage biomass and herbage disappearance (difference) in an open pasture (OP), black walnut silvopasture (BSP), and honeylocust silvopasture (HSP) treatments during the summers of 2020 and 2021 at the Whitethorne Agroforestry Demonstration Center, Blacksburg, VA, USA.
YearTreatments 1Tukey’s Adjusted p-Value
BSPHSPOPSEBSP vs. HSPBSP vs. OPHSP vs. OP
Pre-grazing 2, kg ha−1
202030903380397061.30.0017<0.0001<0.0001
202134004020409069.9<0.0001<0.00010.4373
Post-grazing, kg ha−1
202021302430318080.20.0074<0.0001<0.0001
202124002730331066.60.0056<0.0001<0.0001
Average22702570324053.90.0031<0.0001<0.0001
Difference 2, kg ha−1
202092091081086.00.97180.38670.4060
2021990129077043.4<0.00010.0013<0.0001
1 Treatments: BSP—Black walnut silvopasture; HSP—Honeylocust silvopasture; OP—Open pasture; 2 Significant year x treatment interaction (p-value: Pre-grazing < 0.0001; Difference = 0.0022).
Table 4. Forage nutritive value in an open pasture (OP), black walnut silvopasture (BSP), and honeylocust silvopasture (HSP) treatments during the summers of 2020 and 2021 at the Whitethorne Agroforestry Demonstration Center, Blacksburg, VA, USA.
Table 4. Forage nutritive value in an open pasture (OP), black walnut silvopasture (BSP), and honeylocust silvopasture (HSP) treatments during the summers of 2020 and 2021 at the Whitethorne Agroforestry Demonstration Center, Blacksburg, VA, USA.
YearTreatments 1Tukey’s Adjusted p-Value
BSPHSPOPSEBSP vs. HSPBSP vs. OPHSP vs. OP
CP 2, %
202018.418.619.00.850.89560.64510.7388
202120.821.320.80.850.68350.97440.7067
Average19.619.919.90.670.73460.76000.9730
ADF 2, %
202031.931.630.50.490.66750.09240.1725
202128.627.726.50.960.53370.16830.3997
Average30.329.728.50.450.38350.03420.1238
NDF 3, %
202060.259.658.80.690.55150.20960.4682
202156.454.052.51.630.35030.14440.5306
Average58.356.855.60.950.31230.09790.4251
1 Treatments: BSP—Black walnut silvopasture; HSP—Honeylocust silvopasture; OP—Open pasture. 2 CP—Crude protein; ADF—Acid detergent fiber; 3 NDF—Neutral detergent fiber.
Table 5. Average daily gain (g day−1) of Katahdin ewes that grazed either open pasture (OP), black walnut silvopasture (BSP), or honeylocust silvopasture (HSP) treatments during the summers of 2020 and 2021.
Table 5. Average daily gain (g day−1) of Katahdin ewes that grazed either open pasture (OP), black walnut silvopasture (BSP), or honeylocust silvopasture (HSP) treatments during the summers of 2020 and 2021.
YearTreatments 1Tukey’s Adjusted p-Value
BSPHSPOPSEBSP vs. HSPBSP vs. OPHSP vs. OP
202012.630.019.921.500.43600.74340.6495
2021134.6141.982.816.730.76410.03610.0203
Average73.685.951.311.870.46860.19020.0457
1 Treatments: BSP—Black walnut silvopasture; HSP—Honeylocust silvopasture; OP—Open pasture.
Table 6. Plasma cortisol level (ng ml−1) and hair cortisol level (pg mg−1) of Katahdin ewes that grazed either open pasture (OP), black walnut silvopasture (BSP), or honeylocust silvopasture (HSP) treatments during the summers of 2020 and 2021.
Table 6. Plasma cortisol level (ng ml−1) and hair cortisol level (pg mg−1) of Katahdin ewes that grazed either open pasture (OP), black walnut silvopasture (BSP), or honeylocust silvopasture (HSP) treatments during the summers of 2020 and 2021.
YearTreatments 1Tukey’s Adjusted p-Value
BSPHSPOPSEBSP vs. HSPBSP vs. OPHSP vs. OP
Plasma cortisol, ng ml−1
202020.219.321.11.540.65340.68330.3925
202111.213.317.92.340.53380.04890.1883
Average15.616.419.61.400.69540.04000.0983
Hair cortisol 2, pg mg−1
202021.522.038.92.060.8642<0.0001<0.0001
202124.633.051.12.950.0236<0.0001<0.0001
1 Treatments: BSP—Black walnut silvopasture; HSP—Honeylocust silvopasture; OP—Open pasture. 2 Significant year x treatment interaction (p = 0.0356).
Table 7. Percent of time (0700–2100 h) spent by ewes within different behavior categories in either open pasture (OP), black walnut silvopasture (BSP), or honeylocust silvopasture (HSP) treatments during the summers of 2020 and 2021.
Table 7. Percent of time (0700–2100 h) spent by ewes within different behavior categories in either open pasture (OP), black walnut silvopasture (BSP), or honeylocust silvopasture (HSP) treatments during the summers of 2020 and 2021.
YearTreatments 1Tukey’s Adjusted p-Value
BSPHSPOPSEBSP vs. HSPBSP vs. OPHSP vs. OP
Grazing 2, %
202040.534.833.71.460.01060.00540.6264
202134.125.131.41.800.00140.29500.0192
Standing, %
20205.73.322.91.220.1872<0.0001<0.0001
20214.92.822.91.310.2630<0.0001<0.0001
Average5.33.122.51.360.2434<0.0001<0.0001
Lying, %
202053.561.743.21.770.00320.0008<0.0001
202160.972.145.41.760.0001<0.0001<0.0001
Average57.366.945.11.800.0004<0.0001<0.0001
Drinking Water, %
20200.070.050.090.0270.53020.62140.2864
20210.000.010.140.0340.83060.00950.0158
Average0.030.030.140.0390.95060.05990.0542
Eating Salt 3, %
20200.240.170.100.0610.39950.12510.4350
20210.030.090.110.0420.33800.18880.7232
1 Treatments: BSP—Black walnut silvopasture; HSP—Honeylocust silvopasture; OP—Open pasture. 2 A trend towards significant year x treatment interaction (p = 0.0611). 3 Significant year x treatment interaction (p < 0.0001).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Poudel, S.; Fike, J.H.; Pent, G.J. Hair Cortisol as a Measure of Chronic Stress in Ewes Grazing Either Hardwood Silvopastures or Open Pastures. Agronomy 2022, 12, 1566. https://doi.org/10.3390/agronomy12071566

AMA Style

Poudel S, Fike JH, Pent GJ. Hair Cortisol as a Measure of Chronic Stress in Ewes Grazing Either Hardwood Silvopastures or Open Pastures. Agronomy. 2022; 12(7):1566. https://doi.org/10.3390/agronomy12071566

Chicago/Turabian Style

Poudel, Sanjok, John H. Fike, and Gabriel J. Pent. 2022. "Hair Cortisol as a Measure of Chronic Stress in Ewes Grazing Either Hardwood Silvopastures or Open Pastures" Agronomy 12, no. 7: 1566. https://doi.org/10.3390/agronomy12071566

APA Style

Poudel, S., Fike, J. H., & Pent, G. J. (2022). Hair Cortisol as a Measure of Chronic Stress in Ewes Grazing Either Hardwood Silvopastures or Open Pastures. Agronomy, 12(7), 1566. https://doi.org/10.3390/agronomy12071566

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop