RUS  ENG
Full version
JOURNALS // Computer Research and Modeling // Archive

Computer Research and Modeling, 2025 Volume 17, Issue 5, Pages 941–961 (Mi crm1306)

ANALYSIS AND MODELING OF COMPLEX LIVING SYSTEMS

A minimal model of density-dependent population dynamics incorporating sex structure: simulation and application

O. L. Revutskayaa, G. P. Neverovab, E. Ya. Frismana

a Institute for Complex Analysis of Regional Problems, Far Eastern Branch of RAS, 4 Sholom-Aleikhem st., Birobidzhan, 679016, Russia
b Institute of Automation and Control Processes, Far Eastern Branch of RAS, 5 Radio st., Vladivostok, 690041, Russia

Abstract: This study proposes and analyzes a discrete-time mathematical model of population dynamics with seasonal reproduction, taking into account the density-dependent regulation and sex structure. In the model, population birth rate depends on the number of females, while density is regulated through juvenile survival, which decreases exponentially with increasing total population size. Analytical and numerical investigations of the model demonstrate that when more than half of both females and males survive, the population exhibits stable dynamics even at relatively high birth rates. Oscillations arise when the limitation of female survival exceeds that of male survival. Increasing the intensity of male survival limitation can stabilize population dynamics, an effect particularly evident when the proportion of female offspring is low. Depending on parameter values, the model exhibits stable, periodic, or irregular dynamics, including multistability, where changes in current population size driven by external factors can shift the system between coexisting dynamic modes. To apply the model to real populations, we propose an approach for estimating demographic parameters based on total abundance data. The key idea is to reduce the two-component discrete model with sex structure to a delay equation dependent only on total population size. In this formulation, the initial sex structure is expressed through total abundance and depends on demographic parameters. The resulting one-dimensional equation was applied to describe and estimate demographic characteristics of ungulate populations in the Jewish Autonomous Region. The delay equation provides a good fit to the observed dynamics of ungulate populations, capturing long-term trends in abundance. Point estimates of parameters fall within biologically meaningful ranges and produce population dynamics consistent with field observations. For moose, roe deer, and musk deer, the model suggests predominantly stable dynamics, while annual fluctuations are primarily driven by external factors and represent deviations from equilibrium. Overall, these estimates enable the analysis of structured population dynamics alongside short-term forecasting based on total abundance data.

Keywords: sex structure, density-dependent factors, discrete-time model, parameter estimation, population dynamics

UDC: 51-76:574.34

Received: 11.08.2025
Revised: 26.09.2025
Accepted: 06.10.2025

DOI: 10.20537/2076-7633-2025-17-5-941-961



© Steklov Math. Inst. of RAS, 2026