Project GitHub Repository: https://github.com/kotack1/DukeForest-TackRiosMoyer

List of Figures

  • Figure 1. Deer Sightings by Hour over Three Month Period in Spring 2023
  • Figure 2. Histogram of Deer Group Size
  • Figure 3. Deer Group Size by Time of Day Observed
  • Figure 4. Deer Group Size by Phase of Moon
  • Figure 5. Moon Phases and Time of Day
  • Figure 6. Deer Observations by Month and Time of Day
  • Figure 7. Deer Observations by Month and Moon Phase
  • Figure 8. Heatmap of Deer Observations by Moon Phase and Time of Day
  • Figure 9. Heatmap of Deer Observations by Time of Day and Camera Trap Location

Rationale and Research Questions

Hypotheses

This research focused on the intricate dynamics of white-tailed deer movements within Duke Forest, particularly in relation to the time of day and the position of the moon. The objective was to ascertain the impact of these environmental factors on the spatial behavior of deer populations. Additionally, the study aimed to measure the frequency of deer visits to various areas of the forest, with a specific emphasis on understanding how these patterns correlate with nearby developmental activities. This approach provided valuable insights into the adaptability and movement patterns of white-tailed deer in response to anthropogenic changes in their habitat. The methodology and findings of this study offer a significant contribution to the understanding of wildlife ecology, particularly in areas experiencing urban development (Oleniacz).

  1. Hypothesis 1:
    • H0: The time of day does not have an impact on observed deer
    • Ha: The time of day being dawn/dusk results in an increase of observed deer
  2. Hypothesis 2:
    • H0: Development has no impact on observed deer
    • Ha: An increase in development results in a decrease of observed deer
  3. Hypothesis 3:
    • H0: The phase of the moon has no impact on observed deer
    • Ha: The moon being full/gibbous results in an increase of observed deer
  4. Hypothesis 4:
    • H0: The phase of the moon has no effect on the time of day deer are observed
    • Ha: The moon being full/gibbous results in an increase of observed deer at dawn/dusk

Dataset Information

In this study, data collection was facilitated through the utilization of trail cameras strategically positioned within Duke Forest, a 7,000-acre research and teaching laboratory managed by Duke University (Duke University, n.d.). Under the guidance of Dr. Roberts’ laboratory, known for its expertise in monitoring deer populations and their movements, a total of 50 cameras were deployed along established migration routes within the forest. These cameras were programmed to capture a sequence of ten images over a span of ten seconds whenever motion was detected, continuing this process until no further movement was observed.

Following the retrieval of the cameras, the collected data was uploaded to Wildlife Insights, a platform that leverages artificial intelligence to initially categorize the species captured in the images. To ensure the accuracy of species identification, manual verification and correction by trained personnel were subsequently carried out, amending any miss classifications as necessary. This methodological approach provided a comprehensive and accurate assessment of the wildlife within Duke Forest.

The data set was chosen due to Katie Tack’s position as an assistant to Dr. Sarah Roberts. It is not yet publicly available and was given to our team directly by Dr. Roberts.

Table summary of deer_data structure
Column Name Description
common_name Common known name of species observed
deployment_id Sampling Point of the camera trap which recorded the data
start_date Date animal was first observed
start_time Time animal was first observed
group_size Amount of deer observed in the trap image
month Numerical number of the month
month_name Month name
cam_id Identification number of the camera
division Duke Forest division where the camera is located
TOD Time of day image was recorded - created column
moon_phase Phase of the moon when the image was recorded - created column

Data Wrangling

In this study, a meticulous data wrangling process was employed to refine and optimize the dataset for analysis. Initially, the data underwent a mutation process, wherein superfluous information was systematically filtered out, ensuring that only pertinent data elements were retained. This step was crucial for enhancing the quality and relevance of the dataset, thereby facilitating more accurate and focused analyses. Following this, the streamlined data was strategically merged with an additional dataset containing geographic coordinates. This integration was instrumental in enriching the dataset with spatial context, allowing for more comprehensive and nuanced interpretations of the data, particularly in analyses that required geographic or locational insights. The combination of these data wrangling techniques significantly improved the dataset’s utility for the research objectives, demonstrating the importance of effective data management in the extraction of meaningful insights from complex datasets.

common_name deployment_id start_date start_time group_size month month_name cam_id division TOD moon_phase
White-tailed Deer Sampling Point 59 (KO) 2023-04-10 14 4 4 April 59 (KO) day waning gibbous
White-tailed Deer Sampling Point 73 (DR) 2023-04-26 11 1 4 April 73 (DR) morning waxing crescent
White-tailed Deer Sampling Point 30 (KO) 2023-04-21 02 1 4 April 30 (KO) night waxing crescent
White-tailed Deer Sampling Point 59 (KO) 2023-05-04 17 1 5 May 59 (KO) day waxing gibbous
White-tailed Deer Sampling Point 65 (KO) 2023-04-03 07 2 4 April 65 (KO) morning waxing gibbous
White-tailed Deer Sampling Point 11 (BW) 2023-03-20 18 2 3 March 11 (BW) day waning crescent

Exploratory Analysis

Understanding the Data

Deer Sightings by Hour over Three Months in Spring 2023

Figure 1: Deer Sightings by Hour over Three Months in Spring 2023

The graphical representation delineates discernible peaks in deer sightings throughout the day, with particular emphasis on a notable surge in April during both the early morning and evening hours.

Histogram of Deer Sightings

Figure 2: Histogram of Deer Sightings

The observational data pertaining to deer frequency underscores a prevailing trend wherein a majority of the observed deer exhibited solitary movement. Deer may choose to move solo instead of in a herd for various ecological and behavioral reasons, as supported by scientific literature. One key study conducted by Kjellander and Nordstrom (2003) sheds light on some key factors influencing deer movement patterns, including: Resource Acquisition and Competition, Territorial Behavior and Mating Strategies, Avoidance of Predadtion and Social Dynamics and Dispersal.

It is important to acknowledge that the observational data utilized in this study pertains specifically to the spring of 2023. This temporal parameter bears significance as it may impact herd size, attributable to one or more of the factors elucidated earlier.

Herd Size Scatter Plots

The following two plots look at factors that may affect the size of the herd observed in the Duke Forest. These are the time of day and the moon phase.

Herd Size & Time of Day

Figure 3: Herd Size & Time of Day

The plot shows how the size of the observed herds changes depending on the time of day. There is a small change, showing that larger herd sizes are observed closer to dawn and dusk rather than around noon.

Herd Size by Moon Phase

Figure 4: Herd Size by Moon Phase

As seen here, larger herds, and thus more overall deer, as seen more often during the full moon and the phases surrounding the full moon (waxing gibbous and waning gibbous). Smaller herd sizes are observed around the new moon. This was exptected to be the result as deer are generally more likely to be found in fields during the new moon and in forests during the full moon due to light levels. There does appear to be the outlier in waning crescent, showing that other factors are likely to be affecting observed herd size at the same time as moon phase.

Time of Day Scatter Plots

Deer Sightings by Moon Phase and Time of Day

Figure 5: Deer Sightings by Moon Phase and Time of Day

This plot shows what time of day deer were observed in based on the phase of the moon. Although there were far fewer deer observed during the new moon than the full moon, this data does not show much correlation between moon phase and time of day. Waning crescent has sightings at all hours of the day but waxing gibbous had the fewest difference in time sightings out of all of the moon phases. More data would be needed to determine why the results are not as was anticipated. Temperature, weather, and movement of predators are all factors that may have affected this data in an unpredicted way.

Deer Sightings by Month and Time of Day

Figure 6: Deer Sightings by Month and Time of Day

This plot combines the previous plots and looks at month and time of day as factors for observed herd size. Although more deer were observed in April and May than in March, this may be due to the fact that all of the cameras had not been placed in the forest yet, and thus there was less data to work off of. Taking that into account, herds being observed towards the beginning and end of the day rather than in the middle remains consistent throughout all of the recorded months. This shows that, at least during the spring months, time of day remains a steady factor. Herd size also varies throughout the day in every month, with no consistent pattern of behavior.

Deer Sightings by Month and Moon Phase

Figure 7: Deer Sightings by Month and Moon Phase

Looking at the previous data in a different way, however, shows more of a pattern. Most of the deer observed in March were observed in the afternoon regardless of the phase of moon, but were still mostly present in the afternoon during a full (or near full) moon. During new (or near new) moon, they were consistently observed at all hours.

In April, this changed so that very few deer were observed in the afternoon at all, regardless of moon phase. Instead, they were mainly seen in the morning across the entire month.

In May, the observed behavior became more random with no consistent pattern in any moon phase of time of day.

These changes in behavior may be due to the hours of daylight increasing as it got closer to the Summer solstice.

Heatmaps

Deer Sightings by Moon Phase and Time of Day

Figure 8: Deer Sightings by Moon Phase and Time of Day

This heatmap is another way of looking at the Deer Sightings by Moon Phase and Time of Day chart seen above. It includes herd size into the data, and further pushes forth the unpredictablity of herd size using the factors that were considered in this research. As it does not match up with expected and known deer behavior, it is likely that other factors are affecting this behavior. More research must be conducted with new variables.

Deer Sightings by Time of day and Camera Trap Location

Figure 9: Deer Sightings by Time of day and Camera Trap Location

The final heatmap shows how the observed herd size changes based on time of day and camera trap. Some camera traps saw more movement than others. Some camera traps also recorded larger herd sizes more consistently than others. This may be due to their location near roads, water bodies, structures, etc.

Spatial Analysis

In order to get a better idea of the behavior of deer over the area covered by the camera traps, we also looked at number of different deer signtings at each of the camera traps. the total number of deer sighted in an area was considered to be the sum of the number of deer in the group in each sighting accross all the sightings at the camera trap. The density of deer was considered overall and also compared over different times of day to see if any patternes were visible.

This image of the Duke Forest (Duke University, n.d.) provides a good framework for understanding the locations we will be discussing in this project.

Maps

Map of locations of camera traps in the Duke Forest

Map of number of deer seen in the Duke forest

Map of number of deer seen during the morning

Map of number of deer seen during the afternoon

Map of number of deer seen during the evening

In referencing these maps, there does not appear to be a strong relationship between where time of the day and the location of the deer. Nor does there seem to a relationship between the number of deer present and the closeness to development. Some of the locations with the highest number of deer sightings are in smaller sections of forest or close to human infrastructure, such as the camera trap near the Eubanks Road. This data of course can not give us a complete idea of spatial behavior of the deer because camera traps only give information about presence or absence at very specific locations.

When starting this project, we hoped to conduct a spatial analysis of deer observations and levels of development - however, we could not find the proper data needed to perform this analysis.

Statistical Analysis

Methods

Our statistical approach involves the application of analysis of variance (ANOVA) tests to investigate hypotheses one and three. As a reminder, these null hypotheses poist there will be no relation between deer observations and time of day or moon phases, respectively. ANOVA has been selected as the preferred statistical method for this portion of the project due to its efficacy in evaluating relationships between quantitative and categorical variables. In our specific context, deer observations (group_size) is treated as a quantitative variable, while time of day (TOD) and moon phase (moon_phase) are both considered categorical.

To test hypothesis four, we have selected to use a chi-squared means test. To repeat, this null hypothesis poists that there is no relationship between moon phase and the time of day deer are observed. A Chi-Squared Means test was selected here due to the non-binary categorical nature of the variables under consideration. In our data, there are eight possible moon phases and three possible times of day, which we created earlier in our data wrangling section. Due to how our data has been collected, we are able to treat each entry of time of day as an deer observation additionally (data was only logged when an animal set off the camera). This analysis does not consider the amount of deer observed.

ANOVA and Chi Square Tests

Hypothesis One - The time of day does not have an impact on observed deer

Table 1: One-Way ANOVA
Df Sum Sq Mean Sq F value Pr(>F)
TOD 2 3.726723 1.8633613 3.929286 0.0200383
Residuals 800 379.379130 0.4742239 NA NA

The p-value (Pr(>F)) is less than 0.05, indicating that there is evidence to reject the null hypothesis. Therefore, there is a statistically significant difference in means among the groups defined by TOD. However, the specific interpretation of which groups are different will require further post-hoc tests or examination of the group means.

Hypothesis Three - The phase of the moon has no impact on observed deer

Table 2: One-Way ANOVA
Df Sum Sq Mean Sq F value Pr(>F)
moon_phase 7 3.637784 0.5196834 1.088756 0.3683484
Residuals 795 379.468069 0.4773183 NA NA

The F value of 1.089 is associated with a p-value of 0.368. Since the p-value is greater than the conventional significance level (0.05), we fail to reject the null hypothesis. This suggests that there is no significant difference in the means of the groups based on the moon phases. Based on the results, moon phases do not appear to have a statistically significant effect on the amount of deer observed in the Duke Forest.

Hypothesis Four - The phase of the moon has no effect on the time of day deer are observed

## 
##  Pearson's Chi-squared test
## 
## data:  contigency_table
## X-squared = 237.95, df = 161, p-value = 7.687e-05

The Chi-Square statistic (237.95) is substantial, indicating a substantial deviation from the expected frequencies based on the assumption of independence. The extremely small p-value of 7.687e-05 suggests strong evidence against the null hypothesis of independence. With a p-value below the conventional significance level, we reject the null hypothesis and conclude that there is a significant association between the phase of the moon and the time of day deer are observed.

Visualizing the Models

One-Way ANOVA for Time of Day and Deer Observations

(#tab:ANOVA Model)Tukey Multiple Comparison of Means
diff lwr upr p adj
morning-day 0.0909629 -0.0441878 0.2261137 0.2547153
night -day -0.0754023 -0.2225017 0.0716971 0.4513530
night -morning -0.1663652 -0.3073211 -0.0254094 0.0157434

In summary, the Tukey HSD test suggests that the group_size differs significantly between night and morning, while there is no significant difference between morning and day or night and day.

Chi-Squared Model for Moon Phases and Time of Day

The presented bar plot indicates a discernible influence of moon phase, particularly during the waning gibbous phase, on the temporal patterns of deer sightings. Notably, there is a pronounced surge in deer sightings during daylight hours when the moon is in a waning gibbous phase. While morning and night also exhibit an increase in deer sightings, the magnitude of this increase is comparatively more pronounced during the day.

Summary and Conclusions

Interpretations

The outcomes of our statistical analyses yield discernible insights, prompting the rejection of null hypotheses one and four. Null hypothesis one’s rejection signifies a discernible influence of time of day on deer sightings within the Duke Forest. Employing the Tukey Honest Significant Difference (HSD) test revealed a significant disparity in the number of deer observed between nighttime and morning periods. Although the ANOVA tests did not precisely pinpoint the time period with the highest deer sightings, they collectively provide compelling evidence suggesting increased deer activity during the morning and night compared to the day.

The rejection of null hypothesis four implies a notable association between the moon’s phase and the timing of deer observations in the Duke Forest. Specifically, a conspicuous rise in deer sightings during the daytime is evident during the waning gibbous phase. While morning and night also experience heightened deer sightings, the magnitude is comparatively subdued. Additional observations from the plot further highlight noteworthy observations concerning other moon phases — new moon, last quarter, full moon, and first quarter. Each are seemingly associated with some degree of limitation in deer observations. Particularly intriguing is the observation that the new moon phase corresponds to the lowest incidence of deer sightings. This could potentially be attributed to the absence of ambient light, resulting in total darkness and potentially obscuring the vision of the camera traps, unless equipped with night vision capabilities.

The ANOVA conducted on the data pertaining to deer sightings in relation to lunar phases led to the acceptance of null hypothesis three. The results indicate an absence of statistically significant differences in the means of deer sightings across distinct phases of the moon. Despite an exhaustive literature search, no scholarly articles were identified that either supported or contradicted our specific findings. However, anecdotal evidence sourced from conversations within deer hunting communities suggested an emerging interest in exploring a potential correlation between deer activity and lunar phases, warranting further investigation with an expanded dataset.

Regrettably, hypothesis two, concerning the influence of surrounding development on deer movements, could not be empirically tested due to unavailability of adequate datasets defining development uniformly. The intended spatial analysis, incorporating variables such as housing, highways, and airport proximity, was hindered by the absence of comprehensive and standardized datasets. This limitation underscores the need for future researchers to address the challenges associated with acquiring robust data on development factors. The identification and compilation of such datasets would enable a more comprehensive exploration of the potential influence of development on deer movements around Duke Forest, thereby contributing to the existing body of knowledge in this domain.

Limitations

The present study is subject to certain limitations that warrant consideration. Foremost among these limitations is the relatively modest sample size and the narrow time frame within which observations were conducted. The dataset spans a duration of three months, and the restricted number of observations poses challenges in establishing robust relationships. The limited temporal scope of the study may not capture the full spectrum of factors influencing deer behavior during the spring of 2023.

The brevity of the observational period raises concerns regarding the generalizability of findings, particularly in discerning nuanced patterns or relationships. The intricate interplay of various environmental and ecological factors, beyond the parameters of time of day and moon phase, may have influenced deer behavior during the specified period. The absence of an extended temporal context restricts the ability to differentiate between temporal idiosyncrasies specific to the spring of 2023 and broader trends recurring across multiple spring seasons.

Furthermore, the inability to compare the spring of 2023 with analogous seasons over a more protracted time frame limits the broader applicability of the study’s findings. Without a comprehensive assessment across multiple spring seasons, the outcomes may be confounded by seasonal variations, hindering a conclusive determination of causality. In order to comprehensively explain the factors influencing deer behavior, future studies with extended observation periods and larger datasets encompassing diverse temporal contexts are necessary. Addressing these limitations would contribute to a more nuanced and robust understanding of the dynamics underpinning deer behavior in relation to time of day and moon phase.

Future Analysis

The current investigation lays the groundwork for a continuous and collaborative research initiative, closely aligned with the Duke Forest team’s focus on camera trap data. As the dataset expands, we anticipate an enriched foundation for exploring intricate patterns and relationships within the observed phenomena.

While our primary emphasis rested on the behavior of White-Tail Deer, it is noteworthy that several other species were documented in the observations. This diversity prompts an avenue for future exploration, suggesting a valuable opportunity to compare the findings related to White-Tail Deer with those of other species thriving within the Duke Forest ecosystem. This approach could potentially unveil broader ecological insights and contribute to a more holistic understanding of the dynamics governing wildlife interactions in this location. Spatial considerations also must be made, inviting an exploration of whether the behavior exhibited by White-Tail Deer in the Duke Forest is congruent with that of their counterparts in the surrounding Durham area.


References

Duke University. (n.d.). Duke Forest – Teaching and Research Laboratory. Retrieved from https://dukeforest.duke.edu

Duke University. (n.d.). General Locator Map. Duke Forest. https://dukeforest.duke.edu/recreation/maps/

Oleniacz, L. (2022, October 27). Scientists track triangle deer to learn how they deal with development. NC State University. Retrieved from https://news.ncsu.edu/2022/10/scientists-track-triangle-deer-to-learn-how-they-deal-with-development/

Kjellander, P., & Nordstrom, J. (2003). Cyclic voles, prey switching in red fox, and roe deer dynamics—a test of the alternative prey hypothesis. Oikos, 101(2), 338-344. doi:10.1034/j.1600-0706.2003.12118.x