Behavioral Phenomena and Population Estimation of White-Tailed Deer Based on Camera Trap Data

When shaping management actions and conservation programs, it is essential to understand the population dynamics of a species. One species that sees extensive management and research effort is the white-tailed deer (Odocoileus virginianus), as it is both ecologically and recreationally important in many places. Although there are a multitude of ways to study and estimate the population dynamics of white-tailed deer, one method in particular comes with significant advantages: camera trapping. For this study, nine camera traps were set up at deer scrapes across a private ranch in south Texas. After collecting images during the fall/winter rut of 2015, the cameras’ photos were amassed, and photographed bucks were identified based on unique antler formations. The photo data were then analyzed for a variety of purposes: (a) to determine whether bucks visit scrapes non- randomly; (b) to examine patterns in the ranges of dominant individuals; (c) to assess whether diversity indexes identify fundamental differences in the collections of bucks visiting scrapes. Finally, a novel population estimation technique was applied by adapting and applying the program EstimateS. There were a total of 129 bucks identified over the study region. I determined that bucks did indeed visit camera sites non-randomly and that they visited certain scrapes preferentially. This might have to do with the resources available in specific areas. Additionally, I observed that particular dominant individuals have home ranges that overlap, while others are much more solitary. This suggests that the sites themselves are differentially defensible, and that whether a location is dominated has to do with both the resources it offers and its ability to be monopolized. Results of the population estimation technique suggested that there were a number of bucks that remained undetected by the cameras. Future research could work to ascertain the accuracy of the population estimations made by EstimateS.Integrative Biolog

Deriving Biomass Models for Small-Diameter Loblolly Pine on the Crossett Experimental Forest

Foresters and landowners have a growing interest in carbon sequestration and cellulosic biofuels in southern pine forests, and hence need to be able to accurately predict them. To this end, we derived a set of aboveground biomass models using data from 62 small-diameter loblolly pines (Pinus taeda) sampled on the Crossett Experimental Forest in southeastern Arkansas. Of the 25 equations initially evaluated, we chose 17 that best fit our dataset and compared them using a suite of conventional test statistics, including pseudo-R2 , root mean squared error (RMSE), and bias. Because most of the 17 models varied little in pseudoR 2 (ranging between 0.96 and 0.99), bias (all were within ± 0.01), and RMSE, an additional comparison was done using Akaike’s Information Criterion corrected for small sample size (AICc). This test statistic produced considerably more discrimination between the biomass models. Of the 17 models evaluated, six produced ΔAICc scores that met or exceeded the threshold for substantial support. To recommend a single preferred model, we then extrapolated beyond our actual data and qualitatively compared model predictions with those from the National Biomass Estimator. Our “best” model did not have the minimum AICc score, but rather predicted logically consistent aboveground biomass values at both the upper and lower ends of our extrapolation

Critical exponents of the two-layer Ising model

The symmetric two-layer Ising model (TLIM) is studied by the corner transfer matrix renormalisation group method. The critical points and critical exponents are calculated. It is found that the TLIM belongs to the same universality class as the Ising model. The shift exponent is calculated to be 1.773, which is consistent with the theoretical prediction 1.75 with 1.3% deviation.Comment: 7 pages, with 10 figures include

Efficient simulation of stochastic chemical kinetics with the Stochastic Bulirsch-Stoer extrapolation method

BackgroundBiochemical systems with relatively low numbers of components must be simulated stochastically in order to capture their inherent noise. Although there has recently been considerable work on discrete stochastic solvers, there is still a need for numerical methods that are both fast and accurate. The Bulirsch-Stoer method is an established method for solving ordinary differential equations that possesses both of these qualities.ResultsIn this paper, we present the Stochastic Bulirsch-Stoer method, a new numerical method for simulating discrete chemical reaction systems, inspired by its deterministic counterpart. It is able to achieve an excellent efficiency due to the fact that it is based on an approach with high deterministic order, allowing for larger stepsizes and leading to fast simulations. We compare it to the Euler ?-leap, as well as two more recent ?-leap methods, on a number of example problems, and find that as well as being very accurate, our method is the most robust, in terms of efficiency, of all the methods considered in this paper. The problems it is most suited for are those with increased populations that would be too slow to simulate using Gillespie’s stochastic simulation algorithm. For such problems, it is likely to achieve higher weak order in the moments.ConclusionsThe Stochastic Bulirsch-Stoer method is a novel stochastic solver that can be used for fast and accurate simulations. Crucially, compared to other similar methods, it better retains its high accuracy when the timesteps are increased. Thus the Stochastic Bulirsch-Stoer method is both computationally efficient and robust. These are key properties for any stochastic numerical method, as they must typically run many thousands of simulations

Tracking Aqueous Proton Transfer by Two-Dimensional Infrared Spectroscopy and ab Initio Molecular Dynamics Simulations.

Proton transfer in water is ubiquitous and a critical elementary event that, via proton hopping between water molecules, enables protons to diffuse much faster than other ions. The problem of the anomalous nature of proton transport in water was first identified by Grotthuss over 200 years ago. In spite of a vast amount of modern research effort, there are still many unanswered questions about proton transport in water. An experimental determination of the proton hopping time has remained elusive due to its ultrafast nature and the lack of direct experimental observables. Here, we use two-dimensional infrared spectroscopy to extract the chemical exchange rates between hydronium and water in acid solutions using a vibrational probe, methyl thiocyanate. Ab initio molecular dynamics (AIMD) simulations demonstrate that the chemical exchange is dominated by proton hopping. The observed experimental and simulated acid concentration dependence then allow us to extrapolate the measured single step proton hopping time to the dilute limit, which, within error, gives the same value as inferred from measurements of the proton mobility and NMR line width analysis. In addition to obtaining the proton hopping time in the dilute limit from direct measurements and AIMD simulations, the results indicate that proton hopping in dilute acid solutions is induced by the concerted multi-water molecule hydrogen bond rearrangement that occurs in pure water. This proposition on the dynamics that drive proton hopping is confirmed by a combination of experimental results from the literature

New Detections of Arsenic, Selenium, and Other Heavy Elements in Two Metal-Poor Stars

We use the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope to obtain new high-quality spectra covering the 1900 to 2360 Angstrom wavelength range for two metal-poor stars, HD 108317 and HD 128279. We derive abundances of Cu II, Zn II, As I, Se I, Mo II, and Cd II, which have not been detected previously in either star. Abundances derived for Ge I, Te I, Os II, and Pt I confirm those derived from lines at longer wavelengths. We also derive upper limits from the non-detection of W II, Hg II, Pb II, and Bi I. The mean [As/Fe] ratio derived from these two stars and five others in the literature is unchanged over the metallicity range -2.8 = +0.28 +/- 0.14 (std. dev. = 0.36 dex). The mean [Se/Fe] ratio derived from these two stars and six others in the literature is also constant, = +0.16 +/- 0.09 (std. dev. = 0.26 dex). The As and Se abundances are enhanced relative to a simple extrapolation of the iron-peak abundances to higher masses, suggesting that this mass region (75 < A < 82) may be the point at which a different nucleosynthetic mechanism begins to dominate the quasi-equilibrium alpha-rich freezeout of the iron peak. = +0.56 +/- 0.23 in HD 108317 and HD 128279, and we infer that lines of Cu I may not be formed in local thermodynamic equilibrium in these stars. The [Zn/Fe], [Mo/Fe], [Cd/Fe], and [Os/Fe] ratios are also derived from neutral and ionized species, and each ratio pair agrees within the mutual uncertainties, which range from 0.15 to 0.52 dex.Comment: Accepted for publication in the Astrophysical Journal. 13 pages, 10 figure

QED3_3 with Dynamical Fermions in an External Magnetic Field

In this paper, we present results of numerical lattice simulations of two-flavor QED in three space-time dimensions. First, we provide evidence that chiral symmetry is spontaneously broken in the chiral and continuum limit. Next we discuss the role of an external magnetic field $B$ on the dynamically generated fermion mass. We investigate the $B$-dependence of the condensate through calculations with dynamical fermions using the non-compact formulation of the gauge field, and compare the results with those of a comparable study using the quenched approximation.Comment: 21 pages, 12 figure

Lattice QCD meets experiment in hadron physics

We review recent results in lattice QCD from numerical simulations that allow for a much more realistic QCD vacuum than has been possible before. Comparison with experiment for a variety of hadronic quantities gives agreement to within statistical and systematic errors of 3%. We discuss the implications of this for future calculations in lattice QCD, particularly those which will provide input for B factory experiments.Comment: Review talk at HADRON2003, Aschaffenberg, Germany, September 200