Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BOUT++”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Impact of equilibrium radial electric field on energy loss process after pedestal collapse

An impact of the equilibrium radial electric field on energy loss processes after pedestal collapse is numerically investigated using the BOUT++ framework. Using linear stability analysis, the resistive ballooning mode is shown to be stabilized by the radial shear of the equilibrium radial electric field. On the other hand, the energy loss level after the pedestal collapse increases if the equilibrium radial electric field is taken into account. The spatio-temporal and phase diagram analyses show that the equilibrium radial electric field partially cancels the fluctuation-driven toroidally axisymmetric radial electric field and weakens the E × B shearing rate after pedestal collapse, weakening the turbulence suppression by vortex shearing. Finally, the equilibrium radial electric field therefore increases turbulence intensity in nonlinear cyclic oscillations among pressure gradient, E × B shearing rate, and turbulence intensity, which gives rise to subsequent bursts of turbulent transport and increases the energy loss level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Pedestal origin and extrapolation of high-density small edge-localised-modes peak parallel energy fluence in ITER and SPARC

Experimental analysis and simulations with the BOUT++ code show that small edge-localised modes (ELMs) in reactor-relevant high-density regimes originate in a region close to the separatrix and only marginally perturb the pedestal structure. The measured divertor peak parallel energy fluence (ε ∥,peak ) for a database of small ELM scenarios in DIII-D and ASDEX Upgrade can be reproduced, within 40 % accuracy on average, if an ad hoc modification of the Eich peak parallel ELM energy fluence model is applied to account for the small ELM pedestal birth location. This allows for first-order extrapolation of small-ELM divertor ε ∥,peak to ITER and SPARC, resulting in values that satisfy the nominal melting threshold of tungsten monoblocks of 12 MJ m −2 . The findings reported in this study, both via modelling and direct measurements, constitute a step forward in assessing small ELMs in high edge-collisionality scenarios as a viable plasma regime for the operation of next-generation fusion machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Remarkable 2024 North Atlantic Mid‐Season Hurricane Lull

Abstract The 2024 North Atlantic (hereafter Atlantic) hurricane season started quickly, with the earliest Category 5 on record (Beryl) and three hurricanes forming through 14 August. Following Ernesto's dissipation on 20 August, the Atlantic hurricane season became extremely quiet during the climatological peak of hurricane season, with only one Category 2 hurricane (Francine) and one tropical storm through 23 September. Several environmental factors likely contributed to this unexpected, prolonged lull. During mid‐to‐late August, subseasonal conditions were broadly favorable for Atlantic hurricanes, but a northward shift in African easterly wave emergence latitude yielded fewer tropical cyclone seed disturbances that also traversed unfavorably cool ocean water. During early‐to‐mid September, subseasonal variability driven by the Madden‐Julian oscillation was less conducive to hurricane activity, with several bouts of increased vertical wind shear across the central Atlantic. Throughout most of the lull, the tropical Atlantic was anomalously dry and subsident, suppressing hurricane formation chances. Plain Language Summary The 2024 North Atlantic (hereafter Atlantic) hurricane season was anticipated to be extremely active, with seasonal forecasts calling for a potentially historic season. Following a busy start including the earliest Category 5 Atlantic hurricane on record (Beryl), the season became very quiet between 20 August–23 September, with only one Category 2 hurricane (Francine) and one tropical storm forming. Following the lull, the Atlantic became extremely busy again, with seven hurricanes forming from 25 September through the end of the season–the most on record for this period. This paper focuses on the surprising lull and attributes it to several different factors, including a late‐August northward shift in African easterly wave tracks (organized thunderstorm complexes spanning thousands of kilometers that often serve as Atlantic hurricane seeds). This northward shift brought these systems out over colder water and into a drier airmass, suppressing their hurricane formation chances. During the first 3 weeks of September, a large‐scale tropical phenomenon known as the Madden‐Julian oscillation, resulted in increased sinking motion, drier mid‐levels and increased vertical wind shear–all of which typically reduce Atlantic hurricane activity. Overall, the lull demonstrates how seed disturbances and other environmental factors can strongly modulate aggregate seasonal tropical cyclone activity. Key Points Following an active start, the 2024 Atlantic hurricane was quiet during the climatological peak with only 1 tropical storm and 1 hurricane The lull was not anticipated by seasonal forecasting agencies, which called for an extremely active 2024 Atlantic hurricane season Factors reducing storm count included a northward shift in African easterly waves, broad‐scale subsidence, and reduced mid‐level moisture

Klotzbach, P J↗
DOE OSTI2022

Modeling study of divertor particle and heat flux asymmetries for EAST H-mode discharges

The BOUT++ transport code is run to study the effects of plasma drifts on the divertor out-in asymmetries (DOIAs) of particle and heat fluxes and their decay widths for EAST lower single null H-mode discharges. The diamagnetic drift seems to have no effects on the DOIAs of total particle and heat fluxes due to its divergence-free nature. However, it could significantly increase the DOIAs of peak particle and heat fluxes and the flux decay widths. The E × B drift is found to induce a large plasma flow to the divertor region, enhancing the DOIAs of both total and peak particle and heat fluxes and the flux decay widths. Both the radial and poloidal components of the E × B drift are necessary in increasing the DOIAs, however, the radial E × B drift seems to play a more important role. The effects on the DOIAs caused by both diamagnetic and E × B drifts are reversed with the reverse of toroidal magnetic field. The heat flux decay width λ q and spreading width S q are important physical and engineering parameters for the divertors and could be obtained by fitting the heat flux profiles at divertor targets. The λ q at the outer target from the simulation case with all drifts could well match with the multi-machine scaling proposed by Eich and the DOIA of λ q is in reasonable agreement with the scaling proposed by Goldston.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulations of heat fluxes in an ELMy H-mode discharge on HL-2A

In order to study the distribution and evolution of the transient heat flux on HL-2A during edge-localized-mode (ELM) bursts, the BOUT++ electromagnetic six-field two-fluid model is used to simulate the pedestal collapse under the lower single-null divertor geometry. The equilibrium profiles of HL-2A ELMy H-mode discharge No. 24 953 are adopted as the initial condition in the original case. In this instance, linear analysis shows that the resistive ballooning mode (RBM) and drift-Alfven wave are unstable to this equilibrium, and RBM is the dominant instability. The evolutions of the radial heat fluxes at the outer mid-plane and heat fluxes to the inner and outer targets during the ELM event are presented. Six more equilibria are constructed based on the original case to find out the influence of the pedestal profiles on the peak electron heat flux. The results indicate that the heat flux increases with temperature and/or density, and the theoretical analysis and simulation results consistently show that the heat flux q ∥e is proportional to n e0,SEP $T^{\frac32}_{e0,SEP}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Understanding the negative triangularity ELM trigger and ELM free state on DIII-D with ECE-imaging

Here, the Electron Cyclotron Emission Imaging (ECEI) diagnostic was used to observe a finite-n interchange mode structure in the edge of negative triangularity shaped plasmas on DIII-D. At a small negative triangularity (δ u = -0.2), the plasma is in the H-mode with ELMs that are triggered by a low-n interchange mode. At a larger negative triangularity (δ u = -0.4) and low NBI power (2 MW), a dithering oscillation is observed that is triggered by a low-n interchange mode, whereas at higher NBI power (>2 MW), the edge reverts to L-mode and the low-n interchange mode is present continuously. In all cases, the edge pressure gradient is clamped when the interchange mode is present. It is concluded that the low-n interchange mode prevents the plasma from transitioning to H-mode at a large negative triangularity. This agrees with linear BOUT++ simulations which suggest that the interchange-type MHD can be a resistive ballooning mode whereby resistivity can significantly increase the finite-n ballooning mode growth rate. The absence of H-mode at large negative triangularity can, thus, be explained by the excitation of low-n pressure driven resistive ballooning modes in the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Two-stage crash process in resistive drift ballooning mode driven ELM crash

We report a two-stage crash process in edge localized mode (ELM) driven by resistive drift-ballooning modes (RDBMs) numerically simulated in a full annular torus domain with a scale-separated four-field reduced MHD (RMHD) model using the BOUT++ framework. In the early nonlinear phase, the small first crash is triggered by linearly unstable RDBMs, and m/n=2/1 magnetic islands are nonlinearly excited by nonlinear coupling of RDBMs as well as their higher harmonics. Here, m is the poloidal mode number, n is the toroidal mode number, the q = 2 rational surface exists near the pressure gradient peak, and q is the safety factor. Simultaneously, middle-n RDBM turbulence develops but is poloidally localized around X-points of the magnetic islands, leading to the small energy loss. The second large crash occurs in the late nonlinear phase. Higher harmonics of m/n=2/1 magnetic islands well develop around the q = 2 surface via nonlinear coupling and make the magnetic field stochastic by magnetic island overlapping. Turbulence heat transport develops at X-points of higher harmonics of m/n=2/1 magnetic islands, resulting in the turbulence spreading in the poloidal direction. The large second crash is triggered when the turbulence covers the whole poloidal region so that the magnetic island generation and magnetic field stochastization before the large crash can be interpreted as ELM precursors. It is concluded that the ELM trigger is attributed to the turbulent spreading in the poloidal direction in synchronization with the magnetic field stochastization and the crash is driven by E × B convection rather than the conventional Rechester–Rosenbluth anomalous electron heat transport.

Seto, H. (ORCID:0000000259278643)↗
DOE OSTI2024

Theoretical and global simulation analysis of collisional microtearing modes

Microtearing modes (MTMs) are suggested as a candidate for anomalous thermal transport in tokamak H-mode discharges. This study investigates MTMs in tokamak plasmas, employing simulations in the BOUT++ framework. It simplifies and linearizes the governing equations in detailed linear simulations. The study meticulously evaluates various conductivity models under diverse plasma conditions and collision regimes. The research thoroughly assesses different conductivity models across a range of plasma conditions and collision regimes. A unified dispersion relation that includes both MTM and Drift-Alfvén Wave (DAW) instabilities is derived, showing that DAW and MTM instabilities occur at varying distances from the rational surface. Specifically, MTMs become unstable near the rational surface but stabilize farther away, while drift-Alfvén instability appears farther from the rational surface. The study also re-derives MTM dispersion relations using Ohm's law and the vorticity equation, providing a thorough analysis of electromagnetic and electrostatic interactions in tokamaks. Global simulations demonstrate an inverse correlation between MTM growth rates and collisionality, and a direct correlation with temperature gradients. The nonalignment of the rational surface with the peak ω*e stabilizes the MTMs. Nonlinear simulations highlight electron temperature relaxation as the primary saturation mechanism for MTMs, with magnetic flutter identified as the dominant mode of electron thermal transport.

Fan, K. (ORCID:0000000227518809)↗
DOE OSTI2025

Behavioral resilience via dynamic circuit firing homeostasis

Homeostatic regulation ensures stable neural circuit output under changing conditions. We find that in Drosophila larvae, either presynaptic weakening due to perturbation of transmitter release or postsynaptic weakening due to perturbation of glutamate receptors at synapses between motor neuron (MN) and muscle has little impact on locomotion, suggesting a nonsynaptic compensatory mechanism. In vivo imaging shows that five different forms of synaptic weakening increase the duration of activity bouts in type I MNs. Strikingly, this compensation is input selective: occurring only in the tonic type Ib MN, not the phasic type Is MN that innervates the same muscle. Moreover, an inhibitory class of central pre-MNs that innervates the tonic—but not phasic—input decreases in activity. The adjustment in activity occurs remarkably quickly: within minutes of synapse perturbation. We propose that MN firing is dynamically regulated by two coordinated mechanisms: a cell-autonomous adjustment of MN excitability and a circuit adjustment of inhibitory central drive. The input selectivity of this process suggests homeostatic adjustment to maintain tonic drive but hold constant the phasic drive that organizes locomotory wave patterns.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Observation of quasi-coherent density fluctuation in scrape-off layer enhancing boundary transport in high-βN hybrid plasmas on DIII-D

In this work, we report the observation of a quasi-coherent density fluctuation (QCF) by the Doppler backscattering system in the scrape-off layer (SOL) region of the DIII-D tokamak. This QCF is observed in high power, high performance hybrid plasmas with near double-null divertor (DND) shape during the electron cyclotron heating period. This mode is correlated with a steepened SOL density profile, and leads to significantly elevated particle and heat fluxes between ELMs. The SOL QCF is a long-wavelength ion-scale fluctuation (k θ ρ s ≈0.2-0.4 and k r ρ s ≈0.03), and propagates in the ion diamagnetic direction in the plasma frame. Its radial expanse is about 1.5–2 cm, well beyond the typical width of heat flux λ q on DIII-D. Also, the SOL QCF does not show any clear dependence on the effective SOL collisionality, and thus may raise issues on the control of plasma-material interactions in low collisionality plasmas in which the blob-induced transport is reduced. A linear simulation using BOUT++ with a 5-field reduced model is performed and compared with experimental observations. In simulation results, an interchange-like density perturbation can be driven by the SOL density gradient, and its peak location and the radial width of the density perturbation are in agreement with the experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulations of tokamak edge plasma turbulent fluctuations based on a minimal 3D model

In this study, a new simulation model for tokamak boundary plasma, SOLT3D, is implemented in the BOUT++ framework. The simulation model includes a set of dynamic equations describing collisional boundary plasma and neutral gas in the tokamak scrape-off layer and divertor region. The model is verified against standard linear plasma instabilities and available nonlinear results. For L-mode like conditions, SOLT3D reproduces characteristics of boundary plasma turbulent fluctuations that are similar to published experimental data, in terms of the amplitude and spatial dependence of the fluctuations. It also reproduces realistic plasma fluxes on material surfaces and expected Bohm-like effective radial transport. Plasma fluctuations inferred from the simulations imply inevitably a significant level of intrinsic 'noise' for axisymmetric tokamak plasma transport modeling, introducing errors on the order of unity. In particular, the toroidally averaged atomic rates below 5–10 eV are strongly modified by turbulent plasma fluctuations, which should significantly affect the standard axisymmetric modeling of the tokamak edge plasma and divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

One-dimensional simulation and validation of divertor detachment induced through nitrogen seeding on HL-2A

Divertor detachment is a promising method to solve the power exhaust problem in tokamak devices or even in future magnetic fusion reactors. In this work, a detachment experiment (HL-2A shot #38008) with mixed gas seeding (60% nitrogen and 40% deuterium) is simulated using the SD1D module in BOUT++. In the process from attachment to detachment, the target electron temperature and the target ion saturation current in simulations are found to be consistent with the experimental results measured by Langmuir probes on the target plate. In order to understand the underlying detachment mechanism on HL-2A, this work analyses the role of different particle species in the cases with different seeding rates. It shows that the plasma density varies little and the density of neutrals (D and D 2 ) slightly decreases with the increase of nitrogen seeding rate, such that plasma–neutral interactions cannot effectively reduce plasma energy and plasma momentum in the divertor. The case with a high seeding rate predicts that increasing seeding rate cannot reach a target temperature lower than 2.5 eV, which is the required temperature for strong plasma–molecule interactions. Thus, the plasma–molecule interactions may not be important in the divertor during nitrogen seeding. Here this work also studies the parallel forces including (1) force due to the parallel electric field, (2) friction forces, (3) ion- and electron-thermal forces, and (4) collisional reactions (e.g. charge exchange recombination and ionisation). It is found that the friction force between nitrogen ions and other particle species (primarily D + ) is the dominant force towards the target, while ion- and electron-thermal forces are the dominant force pushing nitrogen ions back to upstream. Parallel forces determine the parallel distribution of nitrogen impurities, and therefore decide the region of nitrogen radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Achieving a robust grassy-ELM operation regime in CFETR

We have identified a robust grassy-edge localized mode (ELM) operation regime for future tokamak reactors. The regime exists within a pedestal top electron collisionality (ν*) window at high global poloidal beta (β p ). The existence of an upper ν* limit for grassy-ELMs is consistent with results previously reported in experiments (Oyama et al 2010 Nucl. Fusion 50 064014), while the existence of a lower ν* limit has not been reported previously. Using EPED and BOUT++, a theoretical model that quantitatively explains the physics of the grassy-ELMs within the window, which distinguishes them from the small mixed-ELMs at lower ν*, is presented for the first time. A peeling-ballooning stability boundary is obtained by scanning the operating density space. The change in density corresponds to a change in ν* that affects the pedestal bootstrap current. High β p leads to a strong Shafranov shift, which affects the flux surface averaged pressure drive. The two effects combine to create a peeling-dominated window in intermediate ν* buffered by ballooning-dominated regimes. Only the peeling-dominated regime shows a cyclic behavior in the perturbed pressure during the nonlinear simulation of an ELM crash, reminiscent of grassy-ELM dynamics. Similarly, the energy released across the separatrix is demonstrated to be significantly smaller. The quick recovery of the ELM crash is explainable by the rapid rise of a low n kink-peeling instability when the pedestal current I ped exceeds a threshold at high β p . It minimizes the excursion beyond marginal stability and is absent in the ballooning-dominated regime. Comparison with recent experiments over a range of β p and ν* strongly supports the physical picture proposed by the modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulation study of particle transport by weakly coherent mode in the Alcator C-Mod tokamak

A simulation study has been conducted of the physical mechanisms behind the weakly coherent mode (WCM) and its produced particle transport in the I-mode edge plasmas by using the BOUT++ code. The WCM is identified in our simulations by its poloidal and radial distributions as well as its frequency and wavenumber spectra. Its produced radial particle flux is calculated and compared with the experimental value. The good agreement indicates that the WCM is an important particle transport channel in the I-mode pedestal. It is found that the WCM can transport particles across the strong outer shear layer of the Er well established in the formation of I-mode, based on which a possible explanation is provided why I-mode does not feature a density pedestal. The key point lies in the change of the cross-phase between the electric potential and density fluctuations induced by the E × B Doppler shift. In the strong shear layer, although the electric potential fluctuation is significantly suppressed, the cross-phase is close to π/2, resulting in a strong drive of the density fluctuation and particle transport. To identify the physical nature of the WCM, a linear dispersion relation for drift Alfvén modes is derived in the slab geometry. A drift Alfvén wave instability is found to have similar dependence to the simulated linear instability behind the WCM on the resistivity and the parallel electron pressure gradient and thermal force terms in the parallel Ohm's law.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electromagnetic turbulence simulation of tokamak edge plasma dynamics and divertor heat load during thermal quench

Abstract The edge plasma turbulence and transport dynamics, as well as the divertor power loads during the thermal-quench phase of tokamak disruptions, are numerically investigated with BOUT++’s flux-driven six-field electromagnetic turbulence model. Here, transient yet intense particle and energy sources are applied at the pedestal top to mimic the plasma power drive at the edge induced by a core thermal collapse, which flattens the core temperature profile. Interesting features, such as surging of divertor heat load (up to 50 times) and broadening of heat-flux width (up to four times) on the outer-divertor target plate, are observed in the simulation, in qualitative agreement with experimental observations. The dramatic changes in divertor heat load and width are due to the enhanced plasma turbulence activities inside the separatrix. Two cross-field transport mechanisms, namely, the E × B turbulent convection and the stochastic parallel advection/conduction, are identified to play important roles in this process. First, an elevated edge pressure gradient drives instabilities and subsequent turbulence in the entire pedestal region. The enhanced turbulence not only transports particles and energy radially across the separatrix via the E × B convection, which causes the initial divertor heat-load burst, but it also induces amplified magnetic fluctuation B ˜ . Once themagnetic fluctuation is large enough to break the magnetic flux surface, magnetic flutter effect provides an additional radial transport channel. In the late stage of our simulation, | B ˜ r / B 0 | reaches to 10 −4 level that completely breaks magnetic flux surfaces such that stochastic field lines are directly connecting pedestal top plasma to the divertor target plates or first wall, further contributing to the divertor heat-flux width broadening.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Exploring the transition from continuous turbulence fluctuations to bursting ELMs in high SOL density regimes

BOUT++ turbulence simulations of the DIII-D reveal that the density profile between the separatrix and pedestal plays a crucial role in the dynamics of edge localized modes (ELMs) and edge plasma turbulent transport. Nonlinear simulations demonstrate that small ELMs in the DIII-D hybrid scenario under high SOL density conditions are predominantly driven by local ballooning modes near the separatrix, stabilizing global instabilities while enhancing localized pressure fluctuations. A key control parameters for ELM dynamics is the separatrix-to-pedestal density ratio, n e,sep /n e,ped . A high ratio indicates a shallow gradient, favoring small ELMs, while a lower ratio signals a steep gradient, which increases the likelihood of large ELMs. Comprehensive parameter scans, including n e,sep /n e,ped , density gradient profiles near the separatrix, and resistivity, reveal the critical role of these parameters in shaping transitions between turbulence-driven transport and ELM bursting. The scans demonstrate that in high SOL density regimes, small ELMs can result from either global resistive MHD instabilities or local ballooning modes near the separatrix, depending on the steepness of the separatrix density gradient. These findings also highlight the transition from continuous turbulence to small ELMs. The post-crash peak in pressure fluctuations, δP rms serves as a critical metric for identifying transition from continuous turbulence fluctuations to ELM bursting. Larger δP rms values correlate with ELM bursts driven by local or global instabilities, whereas smaller values indicate turbulence-dominated transport. Drift-Alfvén and resistive ballooning turbulence enhance the entrainment of fluctuations from the pedestal to the SOL, contributing to the complex interplay of dynamics in this regime. These findings emphasize the importance of separatrix density shaping and pedestal gradient control for optimizing ELM behavior in ITER and future fusion devices.

Li, Nami [Lawrence Livermore National Laboratory (↗
DOE OSTI2022

Estimation of China Fusion Engineering Test Reactor performance and burning fraction in different pellet fueling scenarios by a multi-species radial transport model

Tritium self-sufficiency in future deuterium–tritium fusion reactors is a crucial challenge. As an engineering test reactor, the China Fusion Engineering Test Reactor requires a burning fraction of 3% for the goal to test the accessibility to the future fusion plant. To self-consistently simulate burning plasmas with profile changes in pellet injection scenarios and to estimate the corresponding burning fraction, a one-dimensional multi-species radial transport model is developed in the BOUT++ framework. Here several pellet-fueling scenarios are then tested in the model. The results show that the increased fueling depth improves the burning fraction by particle confinement improvement and fusion power increase. Nevertheless, by increasing the depth, the pellet cooling-down may significantly lower the temperature in the core region. Taking the density perturbation into consideration, the reasonable parameters of the fueling scenario in these simulations are estimated as pellet radius ${r}_{{\rm{p}}}=3\,{\rm{mm}},$ injection rate $=\,4\,\mathrm{Hz},$ and pellet injection velocity $=\,1000\mbox{-}2000\,{\rm{m}}\,{{\rm{s}}}^{-1}$ without drift or $450\,{\rm{m}}\,{{\rm{s}}}^{-1}$ with high-field-side drift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

White-tailed deer responses to acoustic predator cues are contingent upon past land use and contemporary fire regime

Prey can assess the immediate risk of predation by detecting cues of predator presence, and it is expected that prey should invest in costly antipredator behaviors when a cue of predator presence is detected. Features of the habitat in which the cue is detected, such as vegetative concealment, serve as indirect cues of risk and can mediate how prey respond to direct cues of predator presence. Past agricultural land use and contemporary fire regimes are common disturbances that may modify prey perceptions of risk and could, therefore, alter prey responses to direct cues of predator presence. Here we examined whether the overlap of these two disturbances affected white-tailed deer (Odocoileus virginianus) responses to cues of predator presence by measuring deer vigilance and foraging bout duration in response to coyote (Canis latrans) vocalizations across 20 woodlands that varied in past land use and contemporary fire regime. Frequent fire regimes consistently increased deer visibility to predators across both land-use history contexts. Deer exhibited no behavioral response to the predator cue in habitats containing infrequent fire regimes or agricultural legacies. Deer responded to the cue in frequently burned woodlands without agricultural legacies through increased vigilance and time spent at a foraging location. These findings reveal that land-use legacies and contemporary fire regimes can mediate how prey respond to direct cues of risk. They also suggest that prey may balance the uncertainty associated with cues of predation risk with the urgency of responding to a potential attack by being vigilant and remaining in place.

60 APPLIED LIFE SCIENCES↗