-
A Novel Higher Order Artificial Neural Networks
NASA Astrophysics Data System (ADS)
Xu, Shuxiang
2010-05-01
In this paper a new Higher Order Neural Network (HONN) model is introduced and applied in several data mining tasks. Data Mining extracts hidden patterns and valuable information from large databases. A hyperbolic tangent function is used as the neuron activation function for the new HONN model. Experiments are conducted to demonstrate the advantages and disadvantages of the new HONN model, when compared with several conventional Artificial Neural Network (ANN) models: Feedforward ANN with the sigmoid activation function; Feedforward ANN with the hyperbolic tangent activation function; and Radial Basis Function (RBF) ANN with the Gaussian activation function. The experimental results seem to suggest that the new HONN holds higher generalization capability as well as abilities in handling missing data.
-
Spectral approach to homogenization of hyperbolic equations with periodic coefficients
NASA Astrophysics Data System (ADS)
Dorodnyi, M. A.; Suslina, T. A.
2018-06-01
In L2 (Rd ;Cn), we consider selfadjoint strongly elliptic second order differential operators Aε with periodic coefficients depending on x / ε, ε > 0. We study the behavior of the operators cos (Aε1/2 τ) and Aε-1/2 sin (Aε1/2 τ), τ ∈ R, for small ε. Approximations for these operators in the (Hs →L2)-operator norm with a suitable s are obtained. The results are used to study the behavior of the solution vε of the Cauchy problem for the hyperbolic equation ∂τ2 vε = -Aεvε + F. General results are applied to the acoustics equation and the system of elasticity theory.
-
Current from a nano-gap hyperbolic diode using shape-factors: Theory
NASA Astrophysics Data System (ADS)
Jensen, Kevin L.; Shiffler, Donald A.; Peckerar, Martin; Harris, John R.; Petillo, John J.
2017-08-01
Quantum tunneling by field emission from nanoscale features or sharp field emission structures for which the anode-cathode gap is nanometers in scale ("nano diodes") experience strong deviations from the planar image charge lowered tunneling barrier used in the Murphy and Good formulation of the Fowler-Nordheim equation. These deviations alter the prediction of total current from a curved surface. Modifications to the emission barrier are modeled using a hyperbolic (prolate spheroidal) geometry to determine the trajectories along which the Gamow factor in a WKB-like treatment is undertaken; a quadratic equivalent potential is determined, and a method of shape factors is used to evaluate the corrected total current from a protrusion or wedge geometry.
-
Cascaded K-means convolutional feature learner and its application to face recognition
NASA Astrophysics Data System (ADS)
Zhou, Daoxiang; Yang, Dan; Zhang, Xiaohong; Huang, Sheng; Feng, Shu
2017-09-01
Currently, considerable efforts have been devoted to devise image representation. However, handcrafted methods need strong domain knowledge and show low generalization ability, and conventional feature learning methods require enormous training data and rich parameters tuning experience. A lightened feature learner is presented to solve these problems with application to face recognition, which shares similar topology architecture as a convolutional neural network. Our model is divided into three components: cascaded convolution filters bank learning layer, nonlinear processing layer, and feature pooling layer. Specifically, in the filters learning layer, we use K-means to learn convolution filters. Features are extracted via convoluting images with the learned filters. Afterward, in the nonlinear processing layer, hyperbolic tangent is employed to capture the nonlinear feature. In the feature pooling layer, to remove the redundancy information and incorporate the spatial layout, we exploit multilevel spatial pyramid second-order pooling technique to pool the features in subregions and concatenate them together as the final representation. Extensive experiments on four representative datasets demonstrate the effectiveness and robustness of our model to various variations, yielding competitive recognition results on extended Yale B and FERET. In addition, our method achieves the best identification performance on AR and labeled faces in the wild datasets among the comparative methods.
-
Special discontinuities in nonlinearly elastic media
NASA Astrophysics Data System (ADS)
Chugainova, A. P.
2017-06-01
Solutions of a nonlinear hyperbolic system of equations describing weakly nonlinear quasitransverse waves in a weakly anisotropic elastic medium are studied. The influence of small-scale processes of dissipation and dispersion is investigated. The small-scale processes determine the structure of discontinuities (shocks) and a set of discontinuities with a stationary structure. Among the discontinuities with a stationary structure, there are special ones that, in addition to relations following from conservation laws, satisfy additional relations required for the existence of their structure. In the phase plane, the structure of such discontinuities is represented by an integral curve joining two saddles. Special discontinuities lead to nonunique self-similar solutions of the Riemann problem. Asymptotics of non-self-similar problems for equations with dissipation and dispersion are found numerically. These asymptotics correspond to self-similar solutions of the problems.
-
DichotomY IdentitY: Euler-Bernoulli Numbers, Sets-Multisets, FD-BE Quantum-Statistics, 1 /f0 - 1 /f1 Power-Spectra, Ellipse-Hyperbola Conic-Sections, Local-Global Extent: ``Category-Semantics''
NASA Astrophysics Data System (ADS)
Rota, G.-C.; Siegel, Edward Carl-Ludwig
2011-03-01
Seminal Apostol[Math.Mag.81,3,178(08);Am.Math.Month.115,9,795(08)]-Rota[Intro.Prob. Thy.(95)-p.50-55] DichotomY equivalence-class: set-theory: sets V multisets; closed V open; to Abromowitz-Stegun[Hdbk.Math.Fns.(64)]-ch.23,p.803!]: numbers/polynomials generating-functions: Euler V Bernoulli; to Siegel[Schrodinger Cent.Symp.(87); Symp.Fractals, MRS Fall Mtg.,(1989)-5-papers!] power-spectrum: 1/ f {0}-White V 1/ f {1}-Zipf/Pink (Archimedes) HYPERBOLICITY INEVITABILITY; to analytic-geometry Conic-Sections: Ellipse V (via Parabola) V Hyperbola; to Extent/Scale/Radius: Locality V Globality, Root-Causes/Ultimate-Origins: Dimensionality: odd-Z V (via fractal) V even-Z, to Symmetries/(Noether's-theorem connected)/Conservation-Laws Dichotomy: restored/conservation/convergence=0- V broken/non-conservation/divergence=/=0: with asymptotic-limit antipodes morphisms/ crossovers: Eureka!!!; "FUZZYICS"=''CATEGORYICS''!!! Connection to Kummer(1850) Bernoulli-numbers proof of FLT is via Siegel(CCNY;1964) < (1994)[AMS Joint Mtg. (2002)-Abs.973-60-124] short succinct physics proof: FLT = Least-Action Principle!!!
-
Artificial neural network and classical least-squares methods for neurotransmitter mixture analysis.
PubMed
Schulze, H G; Greek, L S; Gorzalka, B B; Bree, A V; Blades, M W; Turner, R F
1995-02-01
Identification of individual components in biological mixtures can be a difficult problem regardless of the analytical method employed. In this work, Raman spectroscopy was chosen as a prototype analytical method due to its inherent versatility and applicability to aqueous media, making it useful for the study of biological samples. Artificial neural networks (ANNs) and the classical least-squares (CLS) method were used to identify and quantify the Raman spectra of the small-molecule neurotransmitters and mixtures of such molecules. The transfer functions used by a network, as well as the architecture of a network, played an important role in the ability of the network to identify the Raman spectra of individual neurotransmitters and the Raman spectra of neurotransmitter mixtures. Specifically, networks using sigmoid and hyperbolic tangent transfer functions generalized better from the mixtures in the training data set to those in the testing data sets than networks using sine functions. Networks with connections that permit the local processing of inputs generally performed better than other networks on all the testing data sets. and better than the CLS method of curve fitting, on novel spectra of some neurotransmitters. The CLS method was found to perform well on noisy, shifted, and difference spectra.
-
Demonstrator of atmospheric reentry system with hyperbolic velocity—DASH
NASA Astrophysics Data System (ADS)
Morita, Yasuhiro; Kawaguchi, Jun'ichiro; Inatani, Yoshifumi; Abe, Takashi
2003-01-01
Among a wide variety of challenging projects planned for the coming decade is the MUSES-C mission designed by the ISAS of Japan. Despite huge amount of data collected by the previous interplanetary spacecraft and probes, the origin and evolution of the solar system still remains unveiled due to their limited information. Thus, our concern has been directed toward a sample return to carry sample from an asteroid back to the earth, which will contribute to better understanding of the system. One of the keys to success is considered the reentry technology with hyperbolic velocity, which has not been demonstrated yet. With this as background, the demonstrator of atmospheric reentry system with hyperbolic velocity, DASH, has been given a commitment to demonstrate the high-speed reentry technology, which will be launched in summer of next year by Japan's H-IIA rocket in a piggyback configuration. The spaceship, composed of a reentry capsule and its carrier, will be injected into a geostationary transfer orbit (GTO) and after several revolutions it will deorbit by burn of a solid propellant deorbit motor. The capsule, identical to that of the sample return mission, can experience the targeted level of thermal environment even from the GTO by tracing a specially designed reentry trajectory.
-
An improved numerical method for the kernel density functional estimation of disperse flow
NASA Astrophysics Data System (ADS)
Smith, Timothy; Ranjan, Reetesh; Pantano, Carlos
2014-11-01
We present an improved numerical method to solve the transport equation for the one-point particle density function (pdf), which can be used to model disperse flows. The transport equation, a hyperbolic partial differential equation (PDE) with a source term, is derived from the Lagrangian equations for a dilute particle system by treating position and velocity as state-space variables. The method approximates the pdf by a discrete mixture of kernel density functions (KDFs) with space and time varying parameters and performs a global Rayleigh-Ritz like least-square minimization on the state-space of velocity. Such an approximation leads to a hyperbolic system of PDEs for the KDF parameters that cannot be written completely in conservation form. This system is solved using a numerical method that is path-consistent, according to the theory of non-conservative hyperbolic equations. The resulting formulation is a Roe-like update that utilizes the local eigensystem information of the linearized system of PDEs. We will present the formulation of the base method, its higher-order extension and further regularization to demonstrate that the method can predict statistics of disperse flows in an accurate, consistent and efficient manner. This project was funded by NSF Project NSF-DMS 1318161.
-
An approach to the development of numerical algorithms for first order linear hyperbolic systems in multiple space dimensions: The constant coefficient case
NASA Technical Reports Server (NTRS)
Goodrich, John W.
1995-01-01
Two methods for developing high order single step explicit algorithms on symmetric stencils with data on only one time level are presented. Examples are given for the convection and linearized Euler equations with up to the eighth order accuracy in both space and time in one space dimension, and up to the sixth in two space dimensions. The method of characteristics is generalized to nondiagonalizable hyperbolic systems by using exact local polynominal solutions of the system, and the resulting exact propagator methods automatically incorporate the correct multidimensional wave propagation dynamics. Multivariate Taylor or Cauchy-Kowaleskaya expansions are also used to develop algorithms. Both of these methods can be applied to obtain algorithms of arbitrarily high order for hyperbolic systems in multiple space dimensions. Cross derivatives are included in the local approximations used to develop the algorithms in this paper in order to obtain high order accuracy, and improved isotropy and stability. Efficiency in meeting global error bounds is an important criterion for evaluating algorithms, and the higher order algorithms are shown to be up to several orders of magnitude more efficient even though they are more complex. Stable high order boundary conditions for the linearized Euler equations are developed in one space dimension, and demonstrated in two space dimensions.
-
Stability of semidiscrete approximations for hyperbolic initial-boundary-value problems: An eigenvalue analysis
NASA Technical Reports Server (NTRS)
Warming, Robert F.; Beam, Richard M.
1986-01-01
A hyperbolic initial-boundary-value problem can be approximated by a system of ordinary differential equations (ODEs) by replacing the spatial derivatives by finite-difference approximations. The resulting system of ODEs is called a semidiscrete approximation. A complication is the fact that more boundary conditions are required for the spatially discrete approximation than are specified for the partial differential equation. Consequently, additional numerical boundary conditions are required and improper treatment of these additional conditions can lead to instability. For a linear initial-boundary-value problem (IBVP) with homogeneous analytical boundary conditions, the semidiscrete approximation results in a system of ODEs of the form du/dt = Au whose solution can be written as u(t) = exp(At)u(O). Lax-Richtmyer stability requires that the matrix norm of exp(At) be uniformly bounded for O less than or = t less than or = T independent of the spatial mesh size. Although the classical Lax-Richtmyer stability definition involves a conventional vector norm, there is no known algebraic test for the uniform boundedness of the matrix norm of exp(At) for hyperbolic IBVPs. An alternative but more complicated stability definition is used in the theory developed by Gustafsson, Kreiss, and Sundstrom (GKS). The two methods are compared.
-
Diffusive instabilities in a hyperbolic activator-inhibitor system with superdiffusion
NASA Astrophysics Data System (ADS)
Mvogo, Alain; Macías-Díaz, Jorge E.; Kofané, Timoléon Crépin
2018-03-01
We investigate analytically and numerically the conditions for wave instabilities in a hyperbolic activator-inhibitor system with species undergoing anomalous superdiffusion. In the present work, anomalous superdiffusion is modeled using the two-dimensional Weyl fractional operator, with derivative orders α ∈ [1,2]. We perform a linear stability analysis and derive the conditions for diffusion-driven wave instabilities. Emphasis is placed on the effect of the superdiffusion exponent α , the diffusion ratio d , and the inertial time τ . As the superdiffusive exponent increases, so does the wave number of the Turing instability. Opposite to the requirement for Turing instability, the activator needs to diffuse sufficiently faster than the inhibitor in order for the wave instability to occur. The critical wave number for wave instability decreases with the superdiffusive exponent and increases with the inertial time. The maximum value of the inertial time for a wave instability to occur in the system is τmax=3.6 . As one of the main results of this work, we conclude that both anomalous diffusion and inertial time influence strongly the conditions for wave instabilities in hyperbolic fractional reaction-diffusion systems. Some numerical simulations are conducted as evidence of the analytical predictions derived in this work.
-
Near-perfect broadband absorption from hyperbolic metamaterial nanoparticles
NASA Astrophysics Data System (ADS)
Riley, Conor T.; Smalley, Joseph S. T.; Brodie, Jeffrey R. J.; Fainman, Yeshaiahu; Sirbuly, Donald J.; Liu, Zhaowei
2017-02-01
Broadband absorbers are essential components of many light detection, energy harvesting, and camouflage schemes. Current designs are either bulky or use planar films that cause problems in cracking and delamination during flexing or heating. In addition, transferring planar materials to flexible, thin, or low-cost substrates poses a significant challenge. On the other hand, particle-based materials are highly flexible and can be transferred and assembled onto a more desirable substrate but have not shown high performance as an absorber in a standalone system. Here, we introduce a class of particle absorbers called transferable hyperbolic metamaterial particles (THMMP) that display selective, omnidirectional, tunable, broadband absorption when closely packed. This is demonstrated with vertically aligned hyperbolic nanotube (HNT) arrays composed of alternating layers of aluminum-doped zinc oxide and zinc oxide. The broadband absorption measures >87% from 1,200 nm to over 2,200 nm with a maximum absorption of 98.1% at 1,550 nm and remains large for high angles. Furthermore, we show the advantages of particle-based absorbers by transferring the HNTs to a polymer substrate that shows excellent mechanical flexibility and visible transparency while maintaining near-perfect absorption in the telecommunications region. In addition, other material systems and geometries are proposed for a wider range of applications.
-
Probing low-energy hyperbolic polaritons in van der Waals crystals with an electron microscope.
PubMed
Govyadinov, Alexander A; Konečná, Andrea; Chuvilin, Andrey; Vélez, Saül; Dolado, Irene; Nikitin, Alexey Y; Lopatin, Sergei; Casanova, Fèlix; Hueso, Luis E; Aizpurua, Javier; Hillenbrand, Rainer
2017-07-21
Van der Waals materials exhibit intriguing structural, electronic, and photonic properties. Electron energy loss spectroscopy within scanning transmission electron microscopy allows for nanoscale mapping of such properties. However, its detection is typically limited to energy losses in the eV range-too large for probing low-energy excitations such as phonons or mid-infrared plasmons. Here, we adapt a conventional instrument to probe energy loss down to 100 meV, and map phononic states in hexagonal boron nitride, a representative van der Waals material. The boron nitride spectra depend on the flake thickness and on the distance of the electron beam to the flake edges. To explain these observations, we developed a classical response theory that describes the interaction of fast electrons with (anisotropic) van der Waals slabs, revealing that the electron energy loss is dominated by excitation of hyperbolic phonon polaritons, and not of bulk phonons as often reported. Thus, our work is of fundamental importance for interpreting future low-energy loss spectra of van der Waals materials.Here the authors adapt a STEM-EELS system to probe energy loss down to 100 meV, and apply it to map phononic states in hexagonal boron nitride, revealing that the electron loss is dominated by hyperbolic phonon polaritons.
-
Sustained fitness gains and variability in fitness trajectories in the long-term evolution experiment with Escherichia coli
PubMed Central
Lenski, Richard E.; Wiser, Michael J.; Ribeck, Noah; Blount, Zachary D.; Nahum, Joshua R.; Morris, J. Jeffrey; Zaman, Luis; Turner, Caroline B.; Wade, Brian D.; Maddamsetti, Rohan; Burmeister, Alita R.; Baird, Elizabeth J.; Bundy, Jay; Grant, Nkrumah A.; Card, Kyle J.; Rowles, Maia; Weatherspoon, Kiyana; Papoulis, Spiridon E.; Sullivan, Rachel; Clark, Colleen; Mulka, Joseph S.; Hajela, Neerja
2015-01-01
Many populations live in environments subject to frequent biotic and abiotic changes. Nonetheless, it is interesting to ask whether an evolving population's mean fitness can increase indefinitely, and potentially without any limit, even in a constant environment. A recent study showed that fitness trajectories of Escherichia coli populations over 50 000 generations were better described by a power-law model than by a hyperbolic model. According to the power-law model, the rate of fitness gain declines over time but fitness has no upper limit, whereas the hyperbolic model implies a hard limit. Here, we examine whether the previously estimated power-law model predicts the fitness trajectory for an additional 10 000 generations. To that end, we conducted more than 1100 new competitive fitness assays. Consistent with the previous study, the power-law model fits the new data better than the hyperbolic model. We also analysed the variability in fitness among populations, finding subtle, but significant, heterogeneity in mean fitness. Some, but not all, of this variation reflects differences in mutation rate that evolved over time. Taken together, our results imply that both adaptation and divergence can continue indefinitely—or at least for a long time—even in a constant environment. PMID:26674951
-
Micro-scale extensional rheometry using hyperbolic converging/diverging channels and jet breakup
PubMed Central
Keshavarz, Bavand
2016-01-01
Understanding the elongational rheology of dilute polymer solutions plays an important role in many biological and industrial applications ranging from microfluidic lab-on-a-chip diagnostics to phenomena such as fuel atomization and combustion. Making quantitative measurements of the extensional viscosity for dilute viscoelastic fluids is a long-standing challenge and it motivates developments in microfluidic fabrication techniques and high speed/strobe imaging of millifluidic capillary phenomena in order to develop new classes of instruments. In this paper, we study the elongational rheology of a family of dilute polymeric solutions in two devices: first, steady pressure-driven flow through a hyperbolic microfluidic contraction/expansion and, second, the capillary driven breakup of a thin filament formed from a small diameter jet (Dj∼O(100 μm)). The small length scale of the device allows very large deformation rates to be achieved. Our results show that in certain limits of low viscosity and elasticity, jet breakup studies offer significant advantages over the hyperbolic channel measurements despite the more complex implementation. Using our results, together with scaling estimates of the competing viscous, elastic, inertial and capillary timescales that control the dynamics, we construct a dimensionless map or nomogram summarizing the operating space for each instrument. PMID:27375824
-
Nonlinear Conservation Laws and Finite Volume Methods
NASA Astrophysics Data System (ADS)
Leveque, Randall J.
Introduction Software Notation Classification of Differential Equations Derivation of Conservation Laws The Euler Equations of Gas Dynamics Dissipative Fluxes Source Terms Radiative Transfer and Isothermal Equations Multi-dimensional Conservation Laws The Shock Tube Problem Mathematical Theory of Hyperbolic Systems Scalar Equations Linear Hyperbolic Systems Nonlinear Systems The Riemann Problem for the Euler Equations Numerical Methods in One Dimension Finite Difference Theory Finite Volume Methods Importance of Conservation Form - Incorrect Shock Speeds Numerical Flux Functions Godunov's Method Approximate Riemann Solvers High-Resolution Methods Other Approaches Boundary Conditions Source Terms and Fractional Steps Unsplit Methods Fractional Step Methods General Formulation of Fractional Step Methods Stiff Source Terms Quasi-stationary Flow and Gravity Multi-dimensional Problems Dimensional Splitting Multi-dimensional Finite Volume Methods Grids and Adaptive Refinement Computational Difficulties Low-Density Flows Discrete Shocks and Viscous Profiles Start-Up Errors Wall Heating Slow-Moving Shocks Grid Orientation Effects Grid-Aligned Shocks Magnetohydrodynamics The MHD Equations One-Dimensional MHD Solving the Riemann Problem Nonstrict Hyperbolicity Stiffness The Divergence of B Riemann Problems in Multi-dimensional MHD Staggered Grids The 8-Wave Riemann Solver Relativistic Hydrodynamics Conservation Laws in Spacetime The Continuity Equation The 4-Momentum of a Particle The Stress-Energy Tensor Finite Volume Methods Multi-dimensional Relativistic Flow Gravitation and General Relativity References
-
On the hyperbolicity and stability of 3+1 formulations of metric f( R) gravity
NASA Astrophysics Data System (ADS)
Mongwane, Bishop
2016-11-01
3+1 formulations of the Einstein field equations have become an invaluable tool in Numerical relativity, having been used successfully in modeling spacetimes of black hole collisions, stellar collapse and other complex systems. It is plausible that similar considerations could prove fruitful for modified gravity theories. In this article, we pursue from a numerical relativistic viewpoint the 3+1 formulation of metric f( R) gravity as it arises from the fourth order equations of motion, without invoking the dynamical equivalence with Brans-Dicke theories. We present the resulting system of evolution and constraint equations for a generic function f( R), subject to the usual viability conditions. We confirm that the time propagation of the f( R) Hamiltonian and Momentum constraints take the same Mathematical form as in general relativity, irrespective of the f( R) model. We further recast the 3+1 system in a form akin to the BSSNOK formulation of numerical relativity. Without assuming any specific model, we show that the ADM version of f( R) is weakly hyperbolic and is plagued by similar zero speed modes as in the general relativity case. On the other hand the BSSNOK version is strongly hyperbolic and hence a promising formulation for numerical simulations in metric f( R) theories.
-
Baseline hearing abilities and variability in wild beluga whales (Delphinapterus leucas).
PubMed
Castellote, Manuel; Mooney, T Aran; Quakenbush, Lori; Hobbs, Roderick; Goertz, Caroline; Gaglione, Eric
2014-05-15
While hearing is the primary sensory modality for odontocetes, there are few data addressing variation within a natural population. This work describes the hearing ranges (4-150 kHz) and sensitivities of seven apparently healthy, wild beluga whales (Delphinapterus leucas) during a population health assessment project that captured and released belugas in Bristol Bay, Alaska. The baseline hearing abilities and subsequent variations were addressed. Hearing was measured using auditory evoked potentials (AEPs). All audiograms showed a typical cetacean U-shape; substantial variation (>30 dB) was found between most and least sensitive thresholds. All animals heard well, up to at least 128 kHz. Two heard up to 150 kHz. Lowest auditory thresholds (35-45 dB) were identified in the range 45-80 kHz. Greatest differences in hearing abilities occurred at both the high end of the auditory range and at frequencies of maximum sensitivity. In general, wild beluga hearing was quite sensitive. Hearing abilities were similar to those of belugas measured in zoological settings, reinforcing the comparative importance of both settings. The relative degree of variability across the wild belugas suggests that audiograms from multiple individuals are needed to properly describe the maximum sensitivity and population variance for odontocetes. Hearing measures were easily incorporated into field-based settings. This detailed examination of hearing abilities in wild Bristol Bay belugas provides a basis for a better understanding of the potential impact of anthropogenic noise on a noise-sensitive species. Such information may help design noise-limiting mitigation measures that could be applied to areas heavily influenced and inhabited by endangered belugas. © 2014. Published by The Company of Biologists Ltd.
-
Role of mass-kill hunting strategies in the extirpation of Persian gazelle (Gazella subgutturosa) in the northern Levant
PubMed Central
Bar-Oz, Guy; Zeder, Melinda; Hole, Frank
2011-01-01
Continuous and intensive exploitation of wildlife resources by early agricultural societies had major ecological consequences in the ancient Near East. In particular, hunting strategies of post-Neolithic societies involving the mass killing of wild ungulates contributed to the eventual extirpation of a number of wild species. A remarkable deposit of bones of Persian gazelle (Gazella subgutarosa) from fourth millennium BCE levels at Tell Kuran in northeastern Syria provides insight into the unsustainable hunting practices that disrupted gazelle migratory patterns and helped set the course for the virtual extinction of this species and possibly other steppe species in the Levant. The social context of mass kills conducted during periods when people relied primarily on domestic livestock for animal resources sets them apart from the more targeted and sustainable practices of earlier periods, when wild animals were the major or sole source of animal protein. PMID:21502520