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].
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.