At least one large-body (diameter > 1.1 km) hypervelocity cratering event occurred during ~ 0.8-0.90 Ma (Zhamanshin, Kazakhstan) in the Middle Pleistocene Transition period. Analysis designed to reduce uncertainty in the dimensions of the Zhamanshin structure employing high resolution topography demonstrated that it likely generated a ~ 26.5 km diameter multi-ring crater. This is at least two times larger than the current best estimates. Using a range of accepted impactor sizes, velocities, compositions, and angles of impact, such impacts typically yield kinetic energies of impact over 240,000 Megatons (TNT). Explosive energetic events of this magnitude (e.g., Yellowstone Caldera) at other times (K-Pg) have created global environmental effects. The factor of two discrepancy in the dimensions of Zhamanshin increases the kinetic energy yield by factors of 7-10, with significantly larger environmental consequences. This justifies examination of rapid climate transitions linked to biological consequences, including those related to environmental perturbations, at ~0.9 Ma.
Large language models can exhibit unexpected behavior in the blink of an eye. In a recent computer use demo, a language model switched from coding to Googling pictures of Yellowstone, and these sudden shifts in behavior have also been observed in reasoning patterns and jailbreaks. This phenomenon is not unique to autoregressive models: in diffusion models, key features of the final output are decided in narrow ``critical windows'' of the generation process. In this work we develop a simple, unifying theory to explain this phenomenon using the formalism of stochastic localization samplers. We show that it emerges generically as the generation process localizes to a sub-population of the distribution it models. While critical windows have been studied at length in diffusion models, existing theory heavily relies on strong distributional assumptions and the particulars of Gaussian diffusion. In contrast to existing work our theory (1) applies to autoregressive and diffusion models; (2) makes no distributional assumptions; (3) quantitatively improves previous bounds even when specialized to diffusions; and (4) requires basic tools and no stochastic calculus or statistical-physics-based
In this study, we investigate the prey predator dynamics of the elk wolf system in northern Yellowstone National Park, USA, using a data driven modeling approach. We used yearly population data for elk and wolves from 1995 to 2022 to construct a mathematical model using a sparse regression modeling framework. To the best of our knowledge, no previous work has applied this framework to capture elk wolf interactions over this time period. Our modeling pipeline integrates Gaussian process regression for data smoothing, sparse identification of nonlinear dynamics for model discovery, and model selection techniques to identify the most suitable mathematical representation. The resulting model is analyzed for its nonlinear dynamics with ecologically meaningful parameters. Stability and bifurcation analyzes are then performed to understand the systems qualitative behavior. A saddle node bifurcation identifies parameter ranges where both species can coexist, while regions outside this range may lead to the extinction of one or both populations. Hopf and saddle node bifurcations together delineate zones of stable co existence, periodic oscillations, and extinction scenarios. Furthermore, co
This study aims to evaluate the accuracy of predicted eruption times of popular geysers in the Yellowstone National Park. The Yellowstone National Park was the first national park in the United States and is known for its geothermal features consisting of many highly popular geysers such as the Old Faithful. Geysers are fascinating to national park visitors because their eruptions could range from small bubbles to jets of water that are hundreds of meters high, and their eruptions could last from seconds to hours. To help tourists plan their visits, the US National Park Service and other independent groups publish predicted eruption times of popular geysers. We hypothesized that the models developed by the US National Park Service are very accurate with little discrepancy from independent analysis, as park rangers monitor the geysers constantly and likely adjust their models over time according to changing conditions underground, and patterns observed. In addition, since researchers in the park likely rely on these predictions, the models would need to be fine-tuned to ensure that no unnecessary effort or resources are wasted in probing the geysers for variables such as temperature
This paper models brucellosis transmission between elk, cattle, and bison, of high conservation value, in the Greater Yellowstone Ecosystem (GYE). It aims to show how landscape changes in the GYE concomitantly impact brucellosis prevalence in the three species. The approach allows us to see how landscape changes in one location influence disease prevalence in populations elsewhere. The model uses the fact that the populations are configured in such a way that elk are an intermediary for disease transmission between cattle and bison. Using landscape ecology metrics applied to the habitat overlaps between elk and cattle and between elk and bison, the landscape parameters are varied to determine how disease propagates throughout the ecosystem as land-use change occurs. Result aim to provide insights into how land management used for the control of disease spread between cattle and elk, may impact (and be impacted by) disease prevalence in bison.
Brucellosis is a zoonotic bacterial infectious disease that affects livestock and wildlife. The Greater Yellowstone Ecosystem is the last area in the United States where cattle are regularly infected with brucellosis. Even though livestock are vaccinated, interactions with reservoir species still result in spill-over cases to cattle. The National Academy of Sciences, Engineering, and Medicine has indicated that modeling efforts should focus on the transmission between elk and cattle and should be conducted to better understand the effects of land-use change and landscape configuration on disease risk. This chapter determines how the landscape's configuration, using the shape and amount of habitat overlap between elk and cattle, contributes to cross-species brucellosis transmission, and how land-use change translates to disease prevalence. A mathematical-epidemiological model is combined with landscape ecology metrics to estimate transmission rates between the species and model disease spread. The results of this study can ultimately help stakeholders develop policy for controlling brucellosis transmission between livestock and elk in the Greater Yellowstone Area. In turn, this coul
Geyser eruption is one of the most popular signature attractions at the Yellowstone National Park. The interdependence of geyser eruptions and impacts of covariates are of interest to researchers in geyser studies. In this paper, we propose a parametric covariate-adjusted recurrent event model for estimating the eruption gap time. We describe a general bivariate recurrent event process, where a bivariate lognormal distribution and a Gumbel copula with different marginal distributions are used to model an interdependent dual-type event system. The maximum likelihood approach is used to estimate model parameters. The proposed method is applied to analyzing the Yellowstone geyser eruption data for a bivariate geyser system and offers a deeper understanding of the event occurrence mechanism of individual events as well as the system as a whole. A comprehensive simulation study is conducted to evaluate the performance of the proposed method.
A dynamical theory of geophysical precipitation pattern formation is presented and applied to irreversible calcium carbonate (travertine) deposition. Specific systems studied here are the terraces and domes observed at geothermal hot springs, such as those at Yellowstone National Park, and speleothems, particularly stalactites and stalagmites. The theory couples the precipitation front dynamics with shallow water flow, including corrections for turbulent drag and curvature effects. In the absence of capillarity and with a laminar flow profile, the theory predicts a one-parameter family of steady state solutions to the moving boundary problem describing the precipitation front. These shapes match well the measured shapes near the vent at the top of observed travertine domes. Closer to the base of the dome, the solutions deviate from observations, and circular symmetry is broken by a fluting pattern, which we show is associated with capillary forces causing thin film break-up. We relate our model to that recently proposed for stalactite growth, and calculate the linear stability spectrum of both travertine domes and stalactites. Lastly, we apply the theory to the problem of precipita
The recent history of The On-Line Encyclopedia of Integer Sequences (or OEIS), describing developments since 2009, and discussing recent sequences involving interesting unsolved problems and in many cases spectacular illustrations. These include: Peaceable Queens, circles in the plane, the earliest cube-free binary sequence, the EKG and Yellowstone permutations, other lexicographically earliest sequences, iteration of number-theoretic functions, home primes and power trains, a memorable prime, a missing prime, Post's tag system, and coordination sequences.
Define a sequence of positive integers by the rule that a(n) = n for 1 <= n <= 3, and for n >= 4, a(n) is the smallest number not already in the sequence which has a common factor with a(n-2) and is relatively prime to a(n-1). We show that this is a permutation of the positive integers. The remarkable graph of this sequence consists of runs of alternating even and odd numbers, interrupted by small downward spikes followed by large upward spikes, suggesting the eruption of geysers in Yellowstone National Park. On a larger scale the points appear to lie on infinitely many distinct curves. There are several unanswered questions concerning the locations of these spikes and the equations for these curves.
Social foraging shows unexpected features such as the existence of a group size threshold to accomplish a successful hunt. Above this threshold, additional individuals do not increase the probability of capturing the prey. Recent direct observations of wolves in Yellowstone Park show that the group size threshold when hunting its most formidable prey, bison, is nearly three times greater than when hunting elk, a prey that is considerably less challenging to capture than bison. These observations provide empirical support to a computational particle model of group hunting which was previously shown to be effective in explaining why hunting success peaks at apparently small pack sizes when hunting elk. The model is based on considering two critical distances between wolves and prey: the minimal safe distance at which wolves stand from the prey, and the avoidance distance at which wolves move away from each other when they approach the prey. The minimal safe distance is longer when the prey is more dangerous to hunt. We show that the model explains effectively that the group size threshold is greater when the minimal safe distance is longer. Although both distances are longer when the
Hierarchical modeling of abundance in space or time using closed-population mark-recapture under heterogeneity (model M$_{h}$) presents two challenges: (i) finding a flexible likelihood in which abundance appears as an explicit parameter and (ii) fitting the hierarchical model for abundance. The first challenge arises because abundance not only indexes the population size, it also determines the dimension of the capture probabilities in heterogeneity models. A common approach is to use data augmentation to include these capture probabilities directly into the likelihood and fit the model using Bayesian inference via Markov chain Monte Carlo (MCMC). Two such examples of this approach are (i) explicit trans-dimensional MCMC, and (ii) superpopulation data augmentation. The superpopulation approach has the advantage of simple specification that is easily implemented in BUGS and related software. However, it reparameterizes the model so that abundance is no longer included, except as a derived quantity. This is a drawback when hierarchical models for abundance, or related parameters, are desired. Here, we analytically compare the two approaches and show that they are more closely relate
As the magma ascends from its depth of generation to the surface, it is often stored in a series of chambers along the way. The rheological contrast between the viscous magma in the magmatic chambers and the surrounding rocks disturbed the stress field which can give rise to various modes of rock failure at magmatic pressures less than the lithostatic stress, leading to an eruption. Different modes of mechanical failure of the chamber walls are considered depending on the geometry and the sign the relative pressure. Relaxation of viscous stress around magmatic chambers, which is important on the time scale of weeks to months is considered in the analysis of stability with application to both large and extra-large magmatic chambers such as Yellowstone. The effects of a strong deep CO2 flux in Yellowstone are considered in detail. The analysis shows that variations in the flow rate around the observed mean value of 40 kg/m2/yr in the hydrothermally active areas can change the composition of the magma for several hundred thousand years, and cause periodic uplift and subsidence of the caldera surface with a period of several decades.
We present a novel and powerful Compressible High-ORder Unstructured Spectral-difference (CHORUS) code for simulating thermal convection and related fluid dynamics in the interiors of stars and planets. The computational geometries are treated as rotating spherical shells filled with stratified gas. The hydrodynamic equations are discretized by a robust and efficient high-order Spectral Difference Method (SDM) on unstructured meshes. The computational stencil of the spectral difference method is compact and advantageous for parallel processing. CHORUS demonstrates excellent parallel performance for all test cases reported in this paper, scaling up to 12,000 cores on the Yellowstone High-Performance Computing cluster at NCAR. The code is verified by defining two benchmark cases for global convection in Jupiter and the Sun. CHORUS results are compared with results from the ASH code and good agreement is found. The CHORUS code creates new opportunities for simulating such varied phenomena as multi-scale solar convection, core convection, and convection in rapidly-rotating, oblate stars.
Gray wolves were reintroduced to Yellowstone National Park (YNP) in 1995. The population initially flourished, but since 2003 the population has experience significant reductions due to factors that may include disease-induced mortality, illegal hunting, park control pro- grams, vehicle induced deaths and intra-species aggression. Despite facing similar conditions, and interference competition with the wolves, the coyote population at YNP has persisted. In this paper we introduce an epidemiological framework that incorporates natural, human-caused and disease-induced mortality as well as interference competition between two species of predators. The outcomes generated by this theoretical framework are used to explore the impact of competition and death-induced mechanisms on predators coexistence. It is the hope that these results on the competitive dynamics of carnivores in Yellowstone National Park will provide park management insights that result in policies that keep the reintroduction of wolves successful.