Search USGSSearch

Geology topics

Matthew J. Gould

Publications and source records attributed to Matthew J. Gould.

5 recordsLinked to original sources

A unified approach to long-term population monitoring of grizzly bears in the Greater Yellowstone Ecosystem

Long-term wildlife research and monitoring programs strive to maintain consistent data collections and analytical methods. Incorporating new techniques is important but can render data sets incongruent and limit their potential to discern trends in demographic parameters. Integrated population models (IPMs) can address these limitations by combining data sources that may span different periods into a unified statistical framework while providing a holistic view of population dynamics. We developed an IPM in a Bayesian framework for grizzly bears ( Ursus arctos ) in the Greater Yellowstone Ecosystem. We coupled demographic data with multiple, independent population count data to link annual changes in abundance with vital rates over 4 decades (1983–2023). Abundance increased threefold from an estimated 270 individuals in 1984 to 1030 individuals in 2023. Parameter estimates indicated survival of bears ≥2 years of age was high, contributing to robust population growth during the 1980s (λ = 1.023 [50 % interquartile range = 0.993–1.082]) and 1990s (λ = 1.064 [1.023–1.103]). A slowing of population growth started around 2000 (2000s: λ = 1.030 [0.989–1.068]) and continued into the 2010s (λ = 1.021 [0.985–1.057]), due primarily to reductions in survival of bears <2 years of age. These findings corroborate previous research that identified density-dependent effects as a likely cause. The IPM framework provided greater certainty and understanding regarding the dynamic demographic characteristics of the population and serves as a powerful monitoring tool for this long-lived species. Implementation of the IPM allows timely dissemination of demographic data to help inform adaptive management strategies and policy decisions necessary for the continued management and conservation of this population. This robust and flexible monitoring system allows scientists to investigate the effects of a changing ecosystem on population dynamics, incorporate new data sources and statistical models, and respond to changes in monitoring needs for the population. We highlight the efficacy of the IPM in estimating and tracking demographic parameters for a long-lived species, while accommodating shifts in monitoring techniques and data collections typical of long-term wildlife conservation programs worldwide.

Idaho, Montana, Wyoming

Pleistocene–Holocene vicariance, not Anthropocene landscape change, explains the genetic structure of American black bear (Ursus americanus) populations in the American Southwest and northern Mexico

The phylogeography of the American black bear ( Ursus americanus ) is characterized by isolation into glacial refugia, followed by population expansion and genetic admixture. Anthropogenic activities, including overharvest, habitat loss, and transportation infrastructure, have also influenced their landscape genetic structure. We describe the genetic structure of the American black bear in the American Southwest and northern Mexico and investigate how prehistoric and contemporary forces shaped genetic structure and influenced gene flow. Using a suite of microsatellites and a sample of 550 bears, we identified 14 subpopulations organized hierarchically following the distribution of ecoregions and mountain ranges containing black bear habitat. The pattern of subdivision we observed is more likely a product of postglacial habitat fragmentation during the Pleistocene and Holocene, rather than a consequence of contemporary anthropogenic barriers to movement during the Anthropocene. We used linear mixed-effects models to quantify the relationship between landscape resistance and genetic distance among individuals, which indicated that both isolation by resistance and geographic distance govern gene flow. Gene flow was highest among subpopulations occupying large tracts of contiguous habitat, was reduced among subpopulations in the Madrean Sky Island Archipelago, where montane habitat exists within a lowland matrix of arid lands, and was essentially nonexistent between two isolated subpopulations. We found significant asymmetric gene flow supporting the hypothesis that bears expanded northward from a Pleistocene refugium located in the American Southwest and northern Mexico and that major highways were not yet affecting gene flow. The potential vulnerability of the species to climate change, transportation infrastructure, and the US–Mexico border wall highlights conservation challenges and opportunities for binational collaboration.

Arizona, New Mexico, Utah, Wyoming

Validating the performance of occupancy models for estimating habitat use and predicting the distribution of highly-mobile species: A case study using the American black bear

Occupancy models have become a valuable tool for estimating wildlife-habitat relationships and for predicting species distributions. Highly-mobile species often violate the assumption that sampling units are geographically closed shifting the probability of occupancy to be interpreted as the probability of use. We used occupancy models, in conjunction with noninvasive sampling, to estimate habitat use and predict the distribution of a highly-mobile carnivore, the American black bear ( Ursus americanus ) in New Mexico, USA. The top model indicated that black bears use areas with higher primary productivity and fewer roads. The predictive performance of such models is rarely validated with independent data, so we validated our model predictions with 2-independent datasets. We first assessed the correlation between predicted and observed habitat use for 28 telemetry-collared bears in the Jemez Mountains. Predicted habitat use was positively correlated with observed use for all 3 years (2012: ρ = 0.81; 2013: ρ = 0.87; 2014: ρ = 0.90). We then predicted the probability of use within a cell where a bear mortality was documented using 2043 mortality locations from sport harvest, depredation, and vehicle collisions. The probability of habitat use at a mortality location was also positively correlated with observed use by the species (2012: ρ = 0.74; 2013: ρ = 0.89; 2014: ρ = 0.93). Our validation procedure supports the notion that occupancy models can be an effective tool for estimating habitat use and predicting the distribution of highly-mobile species when the assumption of geographic closure has been violated. Our findings may be of interest to studies that are estimating habitat use for highly-mobile species that are secretive or rare, difficult to capture, or expensive to monitor with other more intensive methods.

New Mexico

Density of American black bears in New Mexico

Considering advances in noninvasive genetic sampling and spatially explicit capture–recapture (SECR) models, the New Mexico Department of Game and Fish sought to update their density estimates for American black bear ( Ursus americanus ) populations in New Mexico, USA, to aide in setting sustainable harvest limits. We estimated black bear density in the Sangre de Cristo, Sandia, and Sacramento Mountains, New Mexico, 2012–2014. We collected hair samples from black bears using hair traps and bear rubs and used a sex marker and a suite of microsatellite loci to individually genotype hair samples. We then estimated density in a SECR framework using sex, elevation, land cover type, and time to model heterogeneity in detection probability and the spatial scale over which detection probability declines. We sampled the populations using 554 hair traps and 117 bear rubs and collected 4,083 hair samples. We identified 725 (367 male, 358 female) individuals. Our density estimates varied from 16.5 bears/100 km 2 (95% CI = 11.6–23.5) in the southern Sacramento Mountains to 25.7 bears/100 km 2 (95% CI = 13.2–50.1) in the Sandia Mountains. Overall, detection probability at the activity center (g0) was low across all study areas and ranged from 0.00001 to 0.02. The low values of g0 were primarily a result of half of all hair samples for which genotypes were attempted failing to produce a complete genotype. We speculate that the low success we had genotyping hair samples was due to exceedingly high levels of ultraviolet (UV) radiation that degraded the DNA in the hair. Despite sampling difficulties, we were able to produce density estimates with levels of precision comparable to those estimated for black bears elsewhere in the United States.

New Mexico

Estimating black bear density in New Mexico using noninvasive genetic sampling coupled with spatially explicit capture-recapture methods

During the 2004–2005 to 2015–2016 hunting seasons, the New Mexico Department of Game and Fish (NMDGF) estimated black bear abundance (Ursus americanus) across the state by coupling density estimates with the distribution of primary habitat generated by Costello et al. (2001). These estimates have been used to set harvest limits. For example, a density of 17 bears/100 km2 for the Sangre de Cristo and Sacramento Mountains and 13.2 bears/100 km2 for the Sandia Mountains were used to set harvest levels. The advancement and widespread acceptance of non-invasive sampling and mark-recapture methods, prompted the NMDGF to collaborate with the New Mexico Cooperative Fish and Wildlife Research Unit and New Mexico State University to update their density estimates for black bear populations in select mountain ranges across the state. We established 5 study areas in 3 mountain ranges: the northern (NSC; sampled in 2012) and southern Sangre de Cristo Mountains (SSC; sampled in 2013), the Sandia Mountains (Sandias; sampled in 2014), and the northern (NSacs) and southern Sacramento Mountains (SSacs; both sampled in 2014). We collected hair samples from black bears using two concurrent non-invasive sampling methods, hair traps and bear rubs. We used a gender marker and a suite of microsatellite loci to determine the individual identification of hair samples that were suitable for genetic analysis. We used these data to generate mark-recapture encounter histories for each bear and estimated density in a spatially explicit capture-recapture framework (SECR). We constructed a suite of SECR candidate models using sex, elevation, land cover type, and time to model heterogeneity in detection probability and the spatial scale over which detection probability declines. We used Akaike’s Information Criterion corrected for small sample size (AICc) to rank and select the most supported model from which we estimated density. We set 554 hair traps, 117 bear rubs and collected 4,083 hair samples. We identified 725 (367 M, 358 F) individuals; the sex ratio for each study area was approximately equal. Our density estimates varied within and among mountain ranges with an estimated density of 21.86 bears/100 km2 (95% CI: 17.83 – 26.80) for the NSC, 19.74 bears/100 km2 (95% CI: 13.77 – 28.30) in the SSC, 25.75 bears/100 km2 (95% CI: 13.22 – 50.14) in the Sandias, 21.86 bears/100 km2 (95% CI: 17.83 – 26.80) in the NSacs, and 16.55 bears/100 km2 (95% CI: 11.64 – 23.53) in the SSacs. Overall detection probability for hair traps and bear rubs, combined, was low across all study areas and ranged from 0.00001 to 0.02. We speculate that detection probabilities were affected by failure of some hair samples to produce a complete genotype due to UV degradation of DNA, and our inability to set and check some sampling devices due to wildfires in the SSC. Ultraviolet radiation levels are particularly high in New Mexico compared to other states where NGS methods have been used because New Mexico receives substantial amounts of sunshine, is relatively high in elevation (1,200 m – 4,000 m), and is at a lower latitude. Despite these sampling difficulties, we were able to produce density estimates for New Mexico black bear populations with levels of precision comparable to estimated black bear densities made elsewhere in the U.S. Our ability to generate reliable black bear density estimates for 3 New Mexico mountain ranges is attributable to our use of a statistically robust study design and analytical method. There are multiple factors that need to be considered when developing future SECR-based density estimation projects. First, the spatial extent of the population of interest and the smallest average home range size must be determined; these will dictate size of the trapping array and spacing necessary between hair traps. The number of technicians needed and access to the study areas will also influence configuration of the trapping array. We believe shorter sampling occasions could be implemented to reduce degradation of DNA due to UV radiation; this might help increase amplification rates and thereby increase both the number of unique individuals identified and the number of recaptures, improving the precision of the density estimates. A pilot study may be useful to determine the length of time hair samples can remain in the field prior to collection. In addition, researchers may consider setting hair traps and bear rubs in more shaded areas (e.g., north facing slopes) to help reduce exposure to UV radiation. To reduce the sampling interval it will be necessary to either hire more field personnel or decrease the number of hair traps per sampling session. Both of these will enhance detection of long-range movement events by individual bears, increase initial capture and recapture rates, and improve precision of the parameter estimates. We recognize that all studies are constrained by limited resources, however, increasing field personnel would also allow a larger study area to be sampled or enable higher trap density. In conclusion, we estimated the density of black bears in 5 study areas within 3 mountains ranges of New Mexico. Our estimates will aid the NMDGF in setting sustainable harvest limits. Along with estimates of density, information on additional demographic rates (e.g., survival rates and reproduction) and the potential effects that climate change and future land use may have on the demography of black bears may also help inform management of black bears in New Mexico, and may be considered as future areas for research.

New Mexico