In this paper, we propose a spatiotemporal prey-predator model with fear and Allee effects. We first establish the global existence of solution in time and provide some sufficient conditions for the existence of non-negative spatially homogeneous equilibria. Then, we study the stability and bifurcation for the non-negative equilibria and explore the bifurcation diagram, which revealed that the Allee effect and fear factor can induce complex bifurcation scenario. We discuss that large Allee effect-driven Turing instability and pattern transition for the considered system with the Holling-Ⅰ type functional response, and how small Allee effect stabilizes the system in nature. Finally, numerical simulations illustrate the effectiveness of theoretical results. The main contribution of this work is to discover that the Allee effect can induce both codimension-one bifurcations (transcritical, saddle-node, Hopf, Turing) and codimension-two bifurcations (cusp, Bogdanov-Takens and Turing-Hopf) in a spatiotemporal predator-prey model with a fear factor. In addition, we observe that the circular rings pattern loses its stability, and transitions to the coldspot and stripe pattern in Hopf region or the Turing-Hopf region for a special choice of initial condition.
Citation: Yuhong Huo, Gourav Mandal, Lakshmi Narayan Guin, Santabrata Chakravarty, Renji Han. Allee effect-driven complexity in a spatiotemporal predator-prey system with fear factor[J]. Mathematical Biosciences and Engineering, 2023, 20(10): 18820-18860. doi: 10.3934/mbe.2023834
Related Papers:
[1]
Azevedo Joaquim, Mendonça Fábio .
Small scale wind energy harvesting with maximum power tracking. AIMS Energy, 2015, 2(3): 297-315.
doi: 10.3934/energy.2015.3.297
[2]
Nour Khlaifat, Ali Altaee, John Zhou, Yuhan Huang .
A review of the key sensitive parameters on the aerodynamic performance of a horizontal wind turbine using Computational Fluid Dynamics modelling. AIMS Energy, 2020, 8(3): 493-524.
doi: 10.3934/energy.2020.3.493
[3]
Hassam Nasarullah Chaudhry, John Kaiser Calautit, Ben Richard Hughes .
The Influence of Structural Morphology on the Efficiency of Building Integrated Wind Turbines (BIWT). AIMS Energy, 2014, 2(3): 219-236.
doi: 10.3934/energy.2014.3.219
[4]
Sri Kurniati, Sudirman Syam, Arifin Sanusi .
Numerical investigation and improvement of the aerodynamic performance of a modified elliptical-bladed Savonius-style wind turbine. AIMS Energy, 2023, 11(6): 1211-1230.
doi: 10.3934/energy.2023055
[5]
Nadwan Majeed Ali, Handri Ammari .
Design of a hybrid wind-solar street lighting system to power LED lights on highway poles. AIMS Energy, 2022, 10(2): 177-190.
doi: 10.3934/energy.2022010
[6]
Ammar E. Ali, Nicholas C. Libardi, Sohel Anwar, Afshin Izadian .
Design of a compressed air energy storage system for hydrostatic wind turbines. AIMS Energy, 2018, 6(2): 229-244.
doi: 10.3934/energy.2018.2.229
[7]
Rashid Al Badwawi, Mohammad Abusara, Tapas Mallick .
Speed control of synchronous machine by changing duty cycle of DC/DC buck converter. AIMS Energy, 2015, 3(4): 728-739.
doi: 10.3934/energy.2015.4.728
[8]
Ayman B. Attya, T. Hartkopf .
Wind Turbines Support Techniques during Frequency Drops — Energy Utilization Comparison. AIMS Energy, 2014, 2(3): 260-275.
doi: 10.3934/energy.2014.3.260
Sunimerjit Kaur, Yadwinder Singh Brar, Jaspreet Singh Dhillon .
Multi-objective real-time integrated solar-wind-thermal power dispatch by using meta-heuristic technique. AIMS Energy, 2022, 10(4): 943-971.
doi: 10.3934/energy.2022043
Abstract
In this paper, we propose a spatiotemporal prey-predator model with fear and Allee effects. We first establish the global existence of solution in time and provide some sufficient conditions for the existence of non-negative spatially homogeneous equilibria. Then, we study the stability and bifurcation for the non-negative equilibria and explore the bifurcation diagram, which revealed that the Allee effect and fear factor can induce complex bifurcation scenario. We discuss that large Allee effect-driven Turing instability and pattern transition for the considered system with the Holling-Ⅰ type functional response, and how small Allee effect stabilizes the system in nature. Finally, numerical simulations illustrate the effectiveness of theoretical results. The main contribution of this work is to discover that the Allee effect can induce both codimension-one bifurcations (transcritical, saddle-node, Hopf, Turing) and codimension-two bifurcations (cusp, Bogdanov-Takens and Turing-Hopf) in a spatiotemporal predator-prey model with a fear factor. In addition, we observe that the circular rings pattern loses its stability, and transitions to the coldspot and stripe pattern in Hopf region or the Turing-Hopf region for a special choice of initial condition.
1.
Introduction
Can cats or dogs barter? Trade is one of the essential features of human intelligence. Brosnan et al., (2008) reported that chimpanzees, who are an intelligent species like humans, can trade but are reluctant to trade. This subsequently led to major divergence in the fates of the two species. The market is nothing but an expression of human intelligence. As Malone and Bernstein (2015) mentioned, intelligence does not arise only in individual brains; it also arises in groups of individuals. The securities market, including crypto-assets, is the ultimate expression of human intelligence. Fama (1965) stated that the efficiency of the securities market in the strong or semi-strong sense causes instantaneous changes in traders' subjective equilibrium with the incorporation of securities information, and the speed at which they respond to this information determines the winners and losers. In the weak sense, market equilibrium follows changes in traders' subjective equilibria due to noise. The efficiency of the securities market can be summarized as the correct and rapid incorporation of information into prices, the presence of arbitrage opportunities, countless traders always searching for such opportunities, the gradual loss of arbitrage opportunities, and the market becoming more efficient. However, it should be noted that without noise, transactions would not occur, and the price distortion would persist in the presence of strong bias.
In an efficient market where noise and bias have no effect, and information is perfectly symmetrical, security prices should accurately reflect only information. However, if traders had perfectly and simultaneously symmetrical knowledge and information, including the asset valuation model, the transaction would not function because the traders' valuation of the asset would be the same. They must be doing some rational calculations via intelligence in natural and digital computing. In this sense, a rational representative agent in macroeconomics corresponds to the perfect symmetry of information (Lucas (1976), Kirman (1992), Hartley (1996)).
While Black (1986) treated this "symmetry breaking" as noise, the effect of noise on a security's price is expected to be symmetrical based on its nature. However if there exist so many irrational noise traders synchronizing erroneous stochastic beliefs that both affect prices and earn higher expected returns, the unpredictability of noise traders' beliefs creates a risk in the price of the asset that deters rational arbitrageurs from aggressively betting against them. As a result, prices can diverge significantly from fundamental values even in the absence of fundamental risk. Moreover, bearing a disproportionate amount of risk that they themselves create enables noise traders to earn a higher expected return than rational investors (De Long et al., (1990)).
On the other hand, Tversky and Kahneman (1974) found that the effect of bias is asymmetric. They described three heuristics that are employed in making judgments under uncertainty: 1) representativeness, 2) availability of instances or scenarios, and 3) adjustment from an anchor. These heuristics are highly economical and usually effective, but they lead to systematic and predictable errors. However, the effect of these heuristics has not been detected as a global bias in the securities market. Consequently, identifying the specific effects of noise and bias on security prices is challenging.
If attention is paid to any statistical property in any complex system, the log-normal distribution is the most natural and appropriate among the standard or "normal" statistics to overview the whole system (Kobayashi et al., (2011)). Log-normality emerges as a familiar and typical example of statistical aspects in various complex systems. Since every member of any complex system has its own history, each member is in the process of growth (or retrogression). The log-normal distribution is realized because of Gibrat's law, or the Matthew effect. It is applied to cities' size and growth rate, where the proportionate growth process may cause a distribution of city sizes that is log-normal. When considering the entire size distribution, not just the largest cities, the city size distribution is log-normal (Samuels (1965)). However, it has been argued that it is problematic to define cities through their fairly arbitrary legal boundaries. According to Gabaix (1999), Zipf's law is a very tight constraint on the class of admissible models of local growth. It says that for most countries, the size distribution of cities strikingly fits a power law: the number of cities with populations greater than S is proportional to 1/S. Suppose that, at least in the upper tail, all cities follow some proportional growth process (this appears to be verified empirically). This automatically leads their distribution to converge to Zipf's law.
According to Gibrat's law of proportionate effect, the relative change in a firm's size is not affected by its absolute size, which implies that both small and large firms experience similar average rates of growth. However, Samuels (1965) provides evidence that contradicts this law by showing that larger firms grow significantly faster. Aoki and Nirei (2017) developed a neoclassical growth model that produces Pareto's law of income distribution and Zipf's law of firm size distribution based on firm-level productivity shocks. Executives and entrepreneurs invest in both risk-free assets and their own firms' risky stocks, which affect their income and wealth. They used the model to investigate how changes in tax rates can explain the evolution of top incomes in the US. The model accounts for the recent decline in the Pareto exponent of the income distribution and the trend of the top 1 percent income share. In the same research direction, Nirei and Aoki (2016) developed a neoclassical growth model that incorporates heterogeneous households and explains the Pareto distributions of income and wealth in the upper tail. They introduced households' business productivity risks and borrowing constraints in a standard Bewley (1977) model to generate Pareto distributions. Low-productivity households rely on wages and safe asset returns, whereas high-productivity households do not diversify their business risks. The model can accurately explain the observed income distribution in the United States with reasonable parameter values. The authors conducted comparative statics to investigate how changes in parameters affect Pareto distributions and found that the increase in top income dispersion in recent decades may be attributed to changes in the top tax rates in the 1980s. Their analytical finding offers a consistent explanation for numerical comparative statics.
Few studies are available on global bias in the stock market. The fundamental indicators of stocks include information about the effects of noise and bias on stock prices; however, distinguishing between them is generally hard. In this article, I present the fundamentals hypothesis based on rational expectations (Muth (1961)) and, using a log-normal distribution model, detect global bias components from the price-earnings (P/E), price-to-book (P/B) and price-to-cash flow (P/CF) ratios. The traditional theory of the firm assumes that the firm acts in the stockholders' interests and that stockholders are interested in profits such that the object of the firm is to maximize profit. However, there is a range in the profit concept (Bodenhorn (1964)). The analysis results support our fundamentals hypothesis as the detected biases show similar characteristics. Additionally, the results show that the cash flow indicators contain relatively few bias components and are closer to the fundamentals. I further demonstrate and examine why the positive P/IC ratio among the indicators analyzed is a proxy for the fundamentals that do not include bias components. Globally, this implies that the efficient market hypothesis holds. The shape of the curve, and in particular the strength of the bias, is stable throughout time and independent of whether the economy is good or bad. The answer is simple: Cash is a fact, and profit is an opinion. Namely, opinions of management and accountants are added as noise to fundamentals.
2.
Materials and method
2.1. Hypothesis
When Xt denotes the true fundamentals of listed companies at time t, and Rt denotes their growth rates, the following Gibrat process represents the growth of those companies
Xt=RtXt−1.
(1)
It is important to note that I assume that the growth rates Rt are mutually independent random variables that follow the same distribution with finite variance. The initial value of the fundamentals being set as X0 yields
XT=X0T∏t=1Rt
(2)
at time T. Taking the log of both sides of the equation results in
logXT=logX0+logR1+⋯+logRT.
(3)
Therefore,
logXT∼LN(μ,σ2)
(4)
would hold for a sufficiently large T based on the central limit theorem. Essentially, the true fundamental XT of listed companies follows the log-normal distribution.
Furthermore, by assuming rational expectations through the future point in time T as of the present point in time 0 on the premise of a going concern, the following equation becomes true:
X0=E[XT]T∏t=1E[R−1t].
(5)
Therefore, the rational expectations X0 for fundamentals follow the log-normal distribution. In other words, the fundamentals Xt at time t follow the log-normal distribution.
2.2. Data
Among the actual fundamental indicators, we analyze the price-earnings (P/E), forward price-earnings (P/FE), price-to-book (P/B), price-to-cash flows from operating activities (P/OC), price-to-cash flows from investing activities (P/IC) price-to-cash flows from financing activities (P/FC), and the price-to-cash equivalents at the year-end (P/CE) ratios. We divide the cash flow ratios into positive and negative data. We use daily data with absolute values less than 1,000 for all companies listed in Japan for the period spanning 1,817 business days from January 2007 to May 2014 (data source: Nikkei NEEDS https://nkbb.nikkei.co.jp/en/service/nikkei-needs/).
2.3. Test
If our fundamentals hypothesis is true, the actual fundamental indicators should follow the log-normal distribution. To test our hypothesis, we calculate the p-values using the Kolmogorov-Smirnov test, Pearson's χ2 test and Anderson-Darling test, to test the goodness of fit to the log-normal distribution.
Our tests can be strongly affected by the sample size that differs by indicator, so we extract samples at the 300 quartiles from each data set to calculate the p-value for testing the goodness of fit of the overall average of each indicator to the log-normal distribution with variance. For all cases, the null hypothesis is that "the indicator follows the log-normal distribution, " or the otherwise worded "the indicator is a proxy of true fundamentals." The indicator reflects the fundamentals more when the p-value is closer to 1, and when p-value is closer to 0, the bias is stronger. We show the lowest p-value for each test in Table 1 for the data spanning 1,817 business days.
Table 1.
The lowest p-values for the data spanning 1,817 business days. Samples at the 300 quartiles. Significant (bias) levels: 1%(***), 3%(**), 5%(*).
The test results indicate that bias is strong on the P/E, P/FE, and P/B ratios. The forward price-earnings ratio (P/FE) is the most biased. Among all price-to-cash flow ratios, the bias is significant only for the P/OC ratio. It can be easier to manipulate than others. I illustrate the time series of p-values since the p-value is the indicator of the strength of bias. From Table 1, Figure 1, and Figure 2, the positive P/IC ratio is the best proxy for the fundamentals that do not include bias components. Probably, it should be hard to manipulate according to the opinions of management and accounting.
Figure 1.
Time series of the goodness of fit test to the log-normal distribution. Blue: Kolmogorov-Smirnov, Yellow: Pearson's χ2, Blue: Anderson-Darling. The higher the value, the less bias there is, and vice versa.
Figure 4.
Left: GDP-1 to fit P/FE (κ=18.82,α=1,γ=0.385,μ=0.993), Right: GDP-2 to fit P/FE (κ=13.70,α=0.515,γ=0.238,μ=0.993). Note that P/FE has the strongest bias.
The detected biases exhibit similar characteristics. Each indicator is shown in Figure 2 and compared to the log-normal distribution. While the figures only show the data for the first day of all periods, January 4, 2007, the biases remain unchanged. It is worth noting that the shape and strength of the curve are independent of whether the economy is good or bad. The artificial stock market appears to have natural intelligence.
3.
Discussion
The results indicate that the fundamental indicators, including the P/E, P/FE, and P/B ratios, are strongly affected by bias. Bias also has a significant effect on the positive and negative P/OC ratios. Additionally, there is a weak bias on the negative P/IC, positive and negative P/FC, and P/CE ratios.
When we compare the test results, the positive P/IC ratio is the stable proxy of fundamentals among all these indicators. Positive cash flows from investing activities represent the realized gain or loss from past investments such as marketable securities, tangible fixed assets, the sales of investment securities, and income from the collection on loans declared at the end of the accounting period. In other words, positive cash flow is the indicator that most directly reflects past business decisions. We can interpret the year-end cash equivalent ratio and the other cash flow ratios reflecting the fundamentals because they have less bias of the components.
On the other hand, although the P/E ratio and P/B ratio are definite values, investors might not view them as indicators that reflect the fundamentals because there is a high degree of freedom in accounting. The P/OC ratio might also have a lower credibility than other cash flow indicators.
The answer is simple: "Cash is a fact, and profit is an opinion." Namely, opinions of management and accountants are added as noise to true fundamentals. As a result, the Kesten (1973) process
Xt=RtXt−1+ϵt,E[ϵt]>0,
(6)
is realized so that the Pareto distribution is obtained. This result means that a positive bias accompanied their opinions.
In fact, these biases fit the Pareto distribution quite well. The following functions represent generalized Pareto distribution (GPD).
F(x)=1−(1+(x−μκ)1/γ)−α
(7)
f(x)=αγκ−1/γ(x−μ)−1+1γ(1+(x−μκ)1/γ)−α−1
(8)
4.
Conclusions
What are the implications of the results? Does the existence of bias negate the efficient market hypothesis? From the viewpoint of intelligence in natural and digital computing, how should we deal with this problem? The artificial stock market, where countless traders use digital computers to trade, appears to have natural intelligence. Intelligence does not arise only in individual brains but also arises in groups of individuals. The interest in whether markets are efficient is based on the idea that inefficiencies are the reason active returns are possible. If our fundamentals hypothesis is true, the fundamental indicators, including the P/E, P/FE, and P/B ratios, are strongly affected by bias. This result means that investors might not view them as indicators that reflect the fundamentals since they have a high degree of freedom in accounting. Can we conclude that the market is inefficient because of this? This issue has important economic implications. Globally, our result implies that the efficient market hypothesis holds. The shape of the curve, and in particular the strength of the bias, is stable throughout time and independent of whether the economy is good or bad. It is the path of the ants (Kirman (1993); Sano (2022)) to avoid obstacles, and traders move back and forth along it in response to noise (Black (1986); De Long et al., (1990)). This noise, among other things, plays a role in activating the market system, which knows the fundamentals of firms and represents a stable global bias in the stock market. Nevertheless, when the biases revealed by this study disappear, the stock market will become more efficient. It is the responsibility of firms and financial authorities to disclose fair information.
Use of AI tools declaration
The author declares he has not used Artificial Intelligence (AI) tools in the creation of this article.
Acknowledgments
The author is grateful to the Department of Economics of Fukui Prefectural University. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
Conflicts of interest
The author declares no conflicts of interest in this paper.
References
[1]
F. Courchamp, T. Clutton-Brock, B. Grenfell, Inverse density dependence and the Allee effect, Trend. Ecol. Evol., 14 (1999), 405–410. https://doi.org/10.1016/S0169-5347(99)01683-3 doi: 10.1016/S0169-5347(99)01683-3
[2]
P. A. Stephens, W. J. Sutherland, R. P. Freckleton, What is the Allee effect, Oikos, 87 (1999), 185–190. https://doi.org/10.2307/3547011 doi: 10.2307/3547011
[3]
P. A. Stephens, W. J. Sutherland, Consequences of the Allee effect for behaviour, ecology and conservation, Trends Ecol. Evol., 14 (1999), 401–405. https://doi.org/10.1016/s0169-5347(99)01684-5 doi: 10.1016/s0169-5347(99)01684-5
[4]
F. Courchamp, L. Berec, J. Gascoigne, Allee Effects in Ecology and Conservation, Oxford University Press, Oxford, 2008. https://doi.org/10.1093/acprof:oso/9780198570301.001.0001
[5]
D. Scheel, C. Packer, Group hunting behavior of lions: a search for cooperation, Anim. Behav., 41 (1991), 697–709. https://doi.org/10.1016/S0003-3472(05)80907-8 doi: 10.1016/S0003-3472(05)80907-8
[6]
M. T. Alves, F. M. Hilker, Hunting cooperation and Allee effects in predators, J. Theor. Biol., 419 (2017), 13–22. https://doi.org/10.1016/j.jtbi.2017.02.002 doi: 10.1016/j.jtbi.2017.02.002
[7]
R. Han, G. Mandal, L. N. Guin, S. Chakravarty, Dynamical response of a reaction-diffusion predator-prey system with cooperative hunting and prey refuge, J. Stat. Mech., 2022 (2022), 103502. https://doi.org/10.1088/1742-5468/ac946d doi: 10.1088/1742-5468/ac946d
[8]
L. Y. Zanette, A. F. White, M. C. Allen, M. Clinchy, Perceived predation risk reduces the number of offspring songbirds produce per year, Science, 334 (2011), 1398–1401. https://doi.org/10.1126/science.1210908 doi: 10.1126/science.1210908
[9]
X. Y. Wang, L. Y. Zanette, X. F. Zou, Modelling the fear effect in predator-prey interaction, J. Math. Biol., 73 (2016), 1179–1204. https://doi.org/10.1007/s00285-016-0989-1 doi: 10.1007/s00285-016-0989-1
[10]
P. A. Schmidt, L. D. Mech, Wolf pack size and food acquisition, Am. Nat., 150 (1997), 513–517. https://doi.org/10.1086/286079 doi: 10.1086/286079
[11]
F. Courchamp, D. W. Macdonald, Crucial importance of pack size in the African wild dog Lycaon pictus, Anim. Conserv., 4 (2001), 169–174. https://doi.org/10.1017/S1367943001001196 doi: 10.1017/S1367943001001196
[12]
S. Vinoth, R. Sivasamy, K. Sathiyanathan, B. Unyong, G. Rajchakit, R. Vadivel, et al., The dynamics of a Leslie type predator-prey model with fear and Allee effect, Adv. Differ. Equ., 2021 (2021), 1–22. https://doi.org/10.1186/s13662-021-03490-x doi: 10.1186/s13662-021-03490-x
[13]
B. Dennis, Allee effects: population growth, critical density, and the chance of extinction, Nat. Resour. Model., 3 (1989), 481–538. https://doi.org/10.1111/j.1939-7445.1989.tb00119.x doi: 10.1111/j.1939-7445.1989.tb00119.x
[14]
A. D. Bazykin, Nonlinear Dynamics of Interacting Populations, World Scientific Publishing Company Incorporated, Singapore, 1998.
[15]
F. Courchamp, B. T. Grenfell, T. H. Clutton-Brock, Impact of natural enemies on obligately cooperative breeders, Oikos, 91 (2000), 311–322. https://doi.org/10.1034/j.1600-0706.2000.910212.x doi: 10.1034/j.1600-0706.2000.910212.x
[16]
A. Kent, C. P. Doncaster, T. Sluckin, Consequences for predators of rescue and Allee effects on prey, Ecol. Modell., 162 (2003), 233–245. https://doi.org/10.1016/S0304-3800(02)00343-5 doi: 10.1016/S0304-3800(02)00343-5
[17]
S. Zhou, Y. Zhou, G. Wang, The stability of predator-prey systems subject to the Allee effects, Theor. Popul. Biol., 67 (2005), 23–31. https://doi.org/10.1016/j.tpb.2004.06.007 doi: 10.1016/j.tpb.2004.06.007
[18]
K. H. Elliott, G. S. Betini, D. R. Norris, Fear creates an Allee effect: experimental evidence from seasonal populations, Proc. R. Soc. B Biol. Sci., 284 (2017), 20170878. https://doi.org/10.1098/rspb.2017.0878 doi: 10.1098/rspb.2017.0878
[19]
S. K. Sasmal, Population dynamics with multiple Allee effects induced by fear factors-A mathematical on prey-predator interactions, Appl. Math. Modell., 64 (2018), 1–14. https://doi.org/10.1016/j.apm.2018.07.021 doi: 10.1016/j.apm.2018.07.021
[20]
Y. Li, M. He, Z. Li, Dynamics of a ratio-dependent Leslie-Gower predator-prey model with Allee effect and fear effect, Math. Comput. Simul., 201 (2022), 417–439. https://doi.org/10.1016/j.matcom.2022.05.017 doi: 10.1016/j.matcom.2022.05.017
[21]
Y. Shi, J. H. Wu, Q. Cao, Analysis on a diffusive multiple Allee effects predator-prey model induced by fear factors, Nonlinear Anal. Real World Appl., 59 (2021), 103249. https://doi.org/10.1016/j.nonrwa.2020.103249 doi: 10.1016/j.nonrwa.2020.103249
[22]
L. Y. Lai, Z. L. Zhu, F. D. Chen, Stability and bifurcation analysis in a predator-prey model with the additive Allee effect and the fear effect, Mathematics, 8 (2020), 1280. https://doi.org/10.3390/math8081280 doi: 10.3390/math8081280
[23]
J. P. Suraci, M. Clinchy, L. M. Dill, D. Roberts, L. Y. Zanette, Fear of large carnivores causes a trophic cascade, Nat. Commun., 7 (2016), 10698. https://doi.org/10.1038/ncomms10698 doi: 10.1038/ncomms10698
[24]
X. Wang, X. Zou, Modeling the fear effect in predator-prey interactions with adaptive avoidance of predators, Bull. Math. Biol., 79 (2017), 1325–1359. https://doi.org/10.1007/s00285-016-0989-1 doi: 10.1007/s00285-016-0989-1
[25]
A. J. Wirsing, W. J. Ripple, A comparison of shark and wolf research reveals similar behavioral responses by prey, Front. Ecol. Environ., 9 (2011), 335–341. https://doi.org/10.1890/090226 doi: 10.1890/090226
[26]
M. J. Sheriff, C. J. Krebs, R. Boonstra, The sensitive hare: sublethal effects of predator stress on reproduction in snowshoe hares, J. Anim. Ecol., 78 (2009), 1249–1258. https://doi.org/10.1111/j.1365-2656.2009.01552. doi: 10.1111/j.1365-2656.2009.01552
[27]
J. Wang, Y. Cai, S. Fu, W. Wang, The effect of the fear factor on the dynamics of a predator-prey model incorporating the prey refuge, Chaos: Interdiscip. J. Nonlinear Sci., 29 (2019), 083109. https://doi.org/10.1063/1.5111121 doi: 10.1063/1.5111121
[28]
S. Samaddar, M. Dhar, P. Bhattacharya, Effect of fear on prey-predator dynamics: Exploring the role of prey refuge and additional food, Chaos: Interdiscip. J. Nonlinear Sci., 30 (2020), 063129. https://doi.org/10.1063/5.0006968 doi: 10.1063/5.0006968
[29]
B. F. Xie, Z. C. Zhang, Impact of Allee and fear effects in a fractional order prey-predator system incorporating prey refuge, Chaos: Interdiscip. J. Nonlinear Sci., 33 (2023), 013131. https://doi.org/10.1063/5.0130809 doi: 10.1063/5.0130809
[30]
J. Liu, Qualitative analysis of a diffusive predator-prey model with Allee and fear effects, Int. J. Biomath., 14 (2021), 2150037. https://doi.org/10.1142/S1793524521500376 doi: 10.1142/S1793524521500376
[31]
D. Pal, D. Kesh, D. Mukherjee, Qualitative study of cross-diffusion and pattern formation in Leslie-Gower predator-prey model with fear and Allee effects, Chaos, Solitons Fractals, 167 (2023), 113033. https://doi.org/10.1016/j.chaos.2022.113033 doi: 10.1016/j.chaos.2022.113033
[32]
H. Amann, Dynamic theory of quasilinear parabolic equations Ⅱ. reaction-diffusion systems, Differ. Integr. Equations, 3 (1990), 13–75. https://doi.org/10.57262/die/1371586185 doi: 10.57262/die/1371586185
[33]
H. Amann, Dynamics theory of quasilinear parabolic systems Ⅲ. Global existence, Math. Z., 202 (1989), 219–250. https://doi.org/10.1007/BF01215256 doi: 10.1007/BF01215256
[34]
H. Amann, Nonhomogeneous linear and quasiliear elliptic and parabolic boundary value problems, in Function Spaces, Differential Operators and Nonlinear Analysis (eds. H.J. Schmeisser and H. Triebel), Teubner-Texte zur Mathematik, 133 (1993), 9–126. https://doi.org/10.1007/978-3-663-11336-2_1
[35]
S. L. Hollis, R. H. Martin, M. Pierre, Global existence and boundedness in reaction-diffusion systems, SIAM J. Math. Anal., 18 (1987), 744–761. https://doi.org/10.1137/0518057 doi: 10.1137/0518057
[36]
S. Wiggins, Introduction to Applied Nonlinear Dynamical Systems and Chaos, Springer, New York, 2003.
[37]
L. Perko, Differential Equations and Dynamical Systems, Springer-Verlag, New York, 2001.
[38]
Y. A. Kuznetsov, Elements of Applied Bifurcation Theory, Springer-Verlag, New York, 1998.
[39]
W. Govaerts, Y. A. Kuznetsov, H. G. E. Meijer, B. Al-Hdaibat, V. De Witte, A. Dhooge, et al., MATCONT: Continuation toolbox for ODEs in Matlab, Universiteit Gent, Belgium, 2019.
[40]
M. G. Crandall, P. H. Rabinowitz, Bifurcation from simple eigenvalues, J. Funct. Anal., 8 (1971), 321–340. https://doi.org/10.1016/0022-1236(71)90015-2 doi: 10.1016/0022-1236(71)90015-2
[41]
F. Q. Yi, J. J. Wei, J. P. Shi, Bifurcation and spatiotemporal patters in a homogeneous diffusive predator-prey system, J. Differ. Equations, 246 (2009), 1944–1977. https://doi.org/10.1016/j.jde.2008.10.024 doi: 10.1016/j.jde.2008.10.024
[42]
R. Han, S. Dey, M. Banerjee, Spatio-temporal pattern selection in a prey-predator model with hunting cooperation and Allee effect in prey, Chaos, Solitons Fractals, 171 (2023), 113441. https://doi.org/10.1016/j.chaos.2023.113441 doi: 10.1016/j.chaos.2023.113441
[43]
A. Morozov, S. Petrovskii, B. L. Li, Spatiotemporal complexity of patchy invasion in a predator-prey system with the Allee effect, J. Theor. Biol., 238 (2006), 18–35. https://doi.org/10.1016/j.jtbi.2005.05.021 doi: 10.1016/j.jtbi.2005.05.021
[44]
P. R. Chowdhury, S. Petrovskii, M. Banerjee, Oscillation and pattern formation in a slow-fast prey-predator system, Bull. Math. Biol., 83 (2021), 1–41. https://doi.org/10.1007/s11538-021-00941-0 doi: 10.1007/s11538-021-00941-0
[45]
S. Chow, C. Li, D. Wang, Normal Forms and Bifurcation of Plannar Vector Fields, Cambridge University Press, 1994. https://doi.org/10.1017/CBO9780511665639
This article has been cited by:
1.
Z. Kaseb, H. Montazeri,
Data-driven optimization of building-integrated ducted openings for wind energy harvesting: Sensitivity analysis of metamodels,
2022,
258,
03605442,
124814,
10.1016/j.energy.2022.124814
2.
Tesfaye (Shiferaw) Sida,
Will Ethiopia be a springboard or a stonewall for GM crops in Africa?,
2021,
39,
1087-0156,
147,
10.1038/s41587-021-00827-5
Ngwarai Shambira, Golden Makaka, Patrick Mukumba,
Velocity Augmentation Model for an Empty Concentrator-Diffuser-Augmented Wind Turbine and Optimisation of Geometrical Parameters Using Surface Response Methodology,
2024,
16,
2071-1050,
1707,
10.3390/su16041707
5.
Debela Alema Teklemariyem, Eshetu Tadesse Yimer, Venkata Rammaya Ancha, Balewgize Amare Zeru,
A CFD parametric investigation on velocity increment induced by an empty diffuser with circumferential holes for horizontal axis wind turbine applications,
2024,
11,
2331-1916,
10.1080/23311916.2024.2386382
6.
Debela Alema Teklemariyem, Eshetu Tadesse Yimer, Venkata Rammaya Ancha, Balewgize Amare Zeru,
Parametric study of an empty diffuser geometric parameters and shape for a wind turbine using CFD analysis,
2024,
10,
24058440,
e26782,
10.1016/j.heliyon.2024.e26782
7.
Saad A. Mutasher, Husham M. Ahmed,
2023,
CFD Analysis of Brimmed Diffuser Augmented Wind Turbine,
979-8-3503-2709-0,
1,
10.1109/ICETAS59148.2023.10346393
8.
Zahra Adnan Shawket, Suad Hassan Danook,
Overview Improving the Efficiency of a Wind Turbine by Using a Nozzle and Solar Radiation,
2025,
26665190,
100191,
10.1016/j.uncres.2025.100191
Table 1.
The lowest p-values for the data spanning 1,817 business days. Samples at the 300 quartiles. Significant (bias) levels: 1%(***), 3%(**), 5%(*).
Figure 1. Nullclines representation of (3.2) for the system parameter values (i)δ=0.04-no interior equilibrium points; δ=0.06-unique interior equilibrium point; δ=0.05-two interior equilibrium points. (ii)β=0.075-no interior equilibrium points; β=0.01-unique interior equilibrium point; β=0.055-two interior equilibrium points. (iii)γ=0.235-no interior equilibrium points; γ=0.21-unique interior equilibrium point; γ=0.225-two interior equilibrium points. (iv)ρ=0.045-no interior equilibrium points; ρ=0.08-unique interior equilibrium point; ρ=0.055-two interior equilibrium points. The remaining system parameter values for each case are chosen from the following set of parameter values α=0.38,β=0.01,γ=0.21,η=0.3,σ=0.2,ρ=0.08,δ=0.06. Here, the solid blue line is the prey nullcline and the red dashed line is the predator nullcline
Figure 2. Existence of limit cycle for (3.2) at the critical parameter value (i)γ=γ[HB]=0.211421 around the coexistence equilibrium point E∗=(0.352582,0.228873), (ii)δ=δ[HB]=0.058389 around the coexistence equilibrium point E∗=(0.350447,0.233725), (iii)α=α[HB]=0.402793 around E∗=(0.351249,0.226873), (iv)β=β[HB]=0.014360 around E∗=(0.351251,0.226876), (v)ρ=ρ[HB]=0.077660 around E∗=(0.353390,0.230085), (vi)σ=σ[HB]=0.193995 around E∗=(0.351248,0.233894), (vii)η=η[HB]=0.303690 around E∗=(0.350605,0.232377). The remaining parameter values are provided in α=0.38,β=0.01,γ=0.21,ρ=0.08,η=0.3,δ=0.06,σ=0.2 for each case
Figure 3. Identification of saddle-node bifurcation in (3.2) for the specific system parameter value γ[SN]=0.227779 with the coexistence equilibrium point E∗=(0.258868,0.088301). The remaining system parameter values are α=0.38,β=0.01,ρ=0.08,η=0.3,σ=0.2 and δ=0.06
Figure 4. Phase portraits of Bogdanov-Takens bifurcation of (3.2) in the neighborhood of the unique coexistence equilibrium point E∗=(0.403015,0.045993) in the uv-plane. The critical bifurcating parameter values are (γ[BT],σ[BT])=(0.268561,1.324208). For a small perturbation, (i)ϵ1=0,ϵ2=0, (ii)ϵ1=−0.009472,ϵ2=−0.523889, (iii)ϵ1=−0.009472,ϵ2=−0.520889, (iv)ϵ1=−0.01172,ϵ2=−0.523889, (v)ϵ1=−0.007825,ϵ2=−0.523889, (vi)ϵ1=−0.009472,ϵ2=−0.53489. The remaining parameter values are α=0.38,β=0.01,ρ=0.08,η=0.3,δ=0.06
Figure 5. Phase portraits of Bogdanov-Takens bifurcation of (3.2) in the vicinity of the unique coexistence equilibrium point E∗=(0.281675,0.107900) in the uv-plane. The critical bifurcating parameter values are (η[BT],σ[BT])=(0.818409,1.580402). In a small neighborhood (η[BT]+ϵ1,σ[BT]+ϵ2), where (i)ϵ1=0,ϵ2=0, (ii)ϵ1=−0.124587,ϵ2=−0.4438945, (iii)ϵ1=−0.126367,ϵ2=−0.4438945, (iv)ϵ1=−0.146887,ϵ2=−0.4438945, (v)ϵ1=−0.10887,ϵ2=−0.4438945, (vi)ϵ1=−0.124587,ϵ2=−0.4538945. The remaining parameter values are α=0.38,β=0.01,γ=0.21,ρ=0.08,δ=0.06
Figure 6. One parameter bifurcation diagrams corresponding to the system parameter as specified in the x-axis of each figures. The other system parameters used here are provided as (ⅰ) α=0.38,β=0.01,ρ=0.08,η=0.3,δ=0.06,σ=0.2; (ⅱ) α=0.38,β=0.01,γ=0.21,ρ=0.08,η=0.3,σ=0.2; (ⅲ) α=0.38,γ=0.21,ρ=0.08,η=0.3,δ=0.06,σ=0.2. Here, SN represents saddle-node bifurcation point, BP represents transcritical bifurcation point, H represents Hopf bifurcation point
Figure 7. Two parameter bifurcation diagrams corresponding to the system parameter as specified in the axis of each figure. The other system parameter values used here are provided as (ⅰ) α=0.38,β=0.01,ρ=0.08,δ=0.06,σ=0.2; (ⅱ) α=0.38,β=0.01,γ=0.21,ρ=0.08,δ=0.06; (ⅲ) α=0.38,γ=0.21,ρ=0.08,η=0.3,δ=0.06. Here, BT represents Bogdanov-Takens bifurcation point, CP represents cusp bifurcation point
Figure 8. Turing-Hopf bifurcation diagram in d-α plane for different Allee effect threshold values. Here the other system parameters are β=0.048,γ=0.0826,ρ=0.135,η=0.8,δ=0.012,σ=0.128. The Hopf bifurcation threshold values α[HB] are 0.3254 and 0.316 for γ=0.0826 and γ=0.0829, respectively
Figure 9. Pattern selection of the prey species u for different diffusion coefficients in the super-critical bifurcation case 0<g1<g2 with random initial condition. The other system parameters used here are provided as α=0.32,β=0.048,γ=0.0826,ρ=0.135,η=0.8,δ=0.012,σ=0.128 for each case. Here the space stepsizes are chosen as Δx=Δy=1 and the time step size is chosen as Δt=0.001 and dc=37.9638. The parameter set is located in the TR (see, e.g., Figure 8 as a reference)
Figure 10. Turing-Hopf bifurcation diagram and coldspot pattern solutions of the prey species u for different diffusion coefficients in the sub-critical bifurcation case g1<0 and g2<0 with random initial condition. The other system parameters used here are provided as α=0.4,β=0.042,γ=0.1,ρ=0.12,η=0.86,δ=0.06,σ=0.18 for each case. Here the space stepsizes are chosen as Δx=Δy=1 and the time step size is chosen as Δt=0.004 and dc=53.6568. The parameter set is located in the TR
Figure 11. Spatiotemporal patterns of the prey species u at different moments with initial condition (4.21). The other system parameters used here are provided as α=0.5,β=0.042,γ=0.1,ρ=0.12,η=0.86,δ=0.06,σ=0.18 and d=40 for each case. Here, the space step sizes are chosen as Δx=Δy=1 and the time step size is chosen as Δt=0.004 and dc=42.6939. The parameter set is located in the HR
Figure 12. Spatiotemporal patterns of the prey species u at different moments with initial condition (4.21). The other system parameters used here are provided as α=0.5,β=0.042,γ=0.1,ρ=0.12,η=0.86,δ=0.06,σ=0.18 and d=45 for each case. Here the space step sizes are chosen as Δx=Δy=1 and the time stepsize is chosen as Δt=0.004 and dc=42.6939. The parameter set is located in the THR