9.5DSMar 10
The Geometry of Quasi-Cycles: How Stoichiometric Covariance Alters Pre-Bifurcation SignaturesLouis Shuo Wang, Jiguang Yu, Ye Liang et al.
Environmental enrichment can destabilize predator--prey coexistence through a Hopf bifurcation, yet real ecosystems are finite and intrinsically stochastic. We investigate how mechanistically derived demographic noise shapes near-Hopf dynamics in the Rosenzweig--MacArthur model by systematically comparing two diffusion closures that share identical deterministic drift but differ solely in predation-induced covariance structure. Starting from a continuous-time Markov chain description, we derive a full-covariance stochastic differential equation whose diffusion tensor inherits stoichiometric coupling, generating a negative prey--predator cross-covariance. This model is contrasted with a drift-matched diagonal-noise comparator. Using linear noise approximation, Lyapunov analysis, and matrix-valued power spectral density formulations, we propagate local covariance structure through the entire diagnostic chain, including stochastic sensitivity ellipses and a dimensionless noisy-precursor indicator. The results highlight that drift equivalence does not imply covariance equivalence and show how event-level noise geometry influences macroscopic behavior in nonlinear ecological systems. This work integrates bifurcation theory and stochastic analysis to advance multi-scale modeling of complex interacting systems.
86.2OCMar 28
Size-Selective Threshold Harvesting under Nonlocal Crowding and Exogenous RecruitmentJiguang Yu, Louis Shuo Wang, Ye Liang
In this paper, we formulate and analyze an original infinite-horizon bioeconomic optimal control problem for a nonlinear, size-structured fish population. Departing from standard endogenous reproduction frameworks, we model population dynamics using a McKendrick--von Foerster partial differential equation characterized by strictly exogenous lower-boundary recruitment and a nonlocal crowding index. This nonlocal environment variable governs density-dependent individual growth and natural mortality, accurately reflecting the ecological pressures of enhancement fisheries or heavily subsidized stocks. We first establish the existence and uniqueness of the no-harvest stationary profile and introduce a novel intrinsic replacement index tailored to exogenously forced systems, which serves as a vital biological diagnostic rather than a classical persistence threshold. To maximize discounted economic revenue, we derive formal first-order necessary conditions via a Pontryagin-type maximum principle. By introducing a weak-coupling approximation to the adjoint system and applying a single-crossing assumption, we mathematically prove that the optimal size-selective harvesting strategy is a rigorous bang-bang threshold policy. A numerical case study calibrated to an Atlantic cod (\textit{Gadus morhua}) fishery bridges our theoretical framework with applied management. The simulations confirm that the economically optimal minimum harvest size threshold ($66.45$ cm) successfully maintains the intrinsic replacement index above unity, demonstrating that precisely targeted, size-structured harvesting can seamlessly align economic maximization with long-run biological viability.
NANov 19, 2012
Real root refinements for univariate polynomial equationsYe Liang
Real root finding of polynomial equations is a basic problem in computer algebra. This task is usually divided into two parts: isolation and refinement. In this paper, we propose two algorithms LZ1 and LZ2 to refine real roots of univariate polynomial equations. Our algorithms combine Newton's method and the secant method to bound the unique solution in an interval of a monotonic convex isolation (MCI) of a polynomial, and have quadratic and cubic convergence rates, respectively. To avoid the swell of coefficients and speed up the computation, we implement the two algorithms by using the floating-point interval method in Maple15 with the package intpakX. Experiments show that our methods are effective and much faster than the function RefineBox in the software Maple15 on benchmark polynomials.