Persistence and size of seasonal populations on a consumer–resource relationship depends on the allocation strategy toward life-history functions

The long-term ecological dynamics of a population inhabiting a seasonal environment is analyzed using a semi-discrete or impulsive system to represent the consumer–resource interaction. The resource corresponds to an incoming energy flow for consumers that is allocated to reproduction as well as to maintenance in each non-reproductive season. The energy invested in these life-history functions is used in reproductive events, determining the size of the offspring in each reproductive season. Two long-term dynamic patterns are found, resulting in either the persistence or the extinction of the population of consumers. In addition, our model indicates that only one energy allocation strategy provides an optimal combination between individual consumption and long-term population size. The current study contributes to the understanding of how the individual-level and the population-level are interrelated, exhibiting the importance of incorporating phenotypic traits in population dynamics.

The quality and quantity of resources available for consumption are key factors influencing the abundance dynamics of populations: they affect the extent of intraspecific competition and impose upper limits on the population abundance 1,2 . In addition, the availability of resources is likely to affect life history traits, such as growth, reproduction, and survival, which will ultimately also affect population growth 3 . Understanding how organisms allocate their limited resources toward vital functions necessitates the determination of the different traits that define their life-history strategies, and which are expected to evolve in such a way as to maximize fitness 4 . The allocation strategies that have been studied to determine the optimal life-history traits 5-7 that allow maximizing the population growth rate are considered in the Lotka-Euler equation or through its approximation, the net reproductive rate R 0 (i.e., the total number of offspring that an individual produces during its lifetime 3,5,7,8 ). However, few studies have considered the interaction between consumers and resources within a broad ecological approximation, incorporating phenotypic traits such as the energy allocation strategy to understand its effects on population dynamics 9,10 .
Consumer-resource interactions are one of the most important intraspecific relationships in population ecology, with an extensive research agenda encompassing the prey-predator, plant-herbivore, and host-parasite systems 11,12 . Furthermore, these interactions are an underlying component of any food web 13 . For instance, DeSiervo et al. (2020) using this approach and an experimental model, investigate whether the population dynamics of the Arctic mosquito species, Aedes nigripes, is controlled by its food (aquatic biofilms), or its predators (diving Scientific Reports | (2020) 10:21401 | https://doi.org/10.1038/s41598-020-77326-1 www.nature.com/scientificreports/ beetles). Using the traditional differential equation approach for modeling the population dynamics, Bideault et al. (2019) studied the effects of temperature on the consumer resource interaction strength 14 , whose results are highly relevant under the current global warming scenario. One of the main objectives of population ecology is to describe and determine the changes in the number of individuals that exploit a particular resource over time is, which is frequently done through a theoretical approach using consumer-resource models 15 . In recent years, new mathematical modeling approaches have been formulated and used to address new ways of representing the consumer-resource relationship. For instance, the Integral Projection Models (IPM) are adjustment models that estimate the survival, growth, and reproduction rate, allowing the description of how a population structured by variable and continuous states at the individual level (e.g., body size) changes in discrete time [16][17][18] . Similarly, Individual-Based Models (IBM), which were developed through computer programming, pursue the determination of the impact of a phenotype on the population dynamics and the strategies that emerge as a product of abiotic factors [19][20][21] . These modeling approaches have allowed a deeper understanding of the interaction between individual traits and population dynamics 22,23 .
The incorporation of physiological and energetic principles at the individual level enables the emergence of dynamic patterns at the population level. Thus, the theory of Dynamic Energy Storage (DEB) 24 and the Metabolic Theory of Ecology (MTE) 25,26 provide a framework focused at the individual level. For their formulation, both these theories consider environmental variables and physiological principles related to the optimal acquisition, ingestion, and resource allocation towards the functions of reproduction and somatic growth, which are essential components for the implementation of both IPM and IBM.
Among the mathematical models representing the consumer-resource relationship based on ordinary differential equations [27][28][29] there stand out both biomass conversion (BC) and individual survival (IS) models. These models differ in the representation of the per capita population growth rate of the consumers 30 .The first assumes a rate dependent on the consumption of resources, therefore relating the decrease in the density of resources with the increase in consumer biomass. The second considers a rate dependent on the density of consumers restricted by intraspecific interactions (e.g., in the logistic model, the per capita population growth rate is inversely proportional to the intraspecific competition factor). It is important to note that for both models, the life-history traits of the organisms are a common underlying component in the formulation of the per capita population growth rate (e.g., predation efficiency, prey search and management time, the average number of progeny, probability of survival until some specific age).
The classical approach to the mathematical modeling of biological systems provided by differential (or difference) equations has considered the resource allocation to different functions of life history (e.g., growth and reproduction) to predict optimal life-history traits (e.g., age and body size of maturity) 5,6,31 , and also to determine the optimal strategies that maximize the lifetime reproductive success of individuals. However, these results, focused on the individual level, have not been scaled at a population level: in particular, they ignore the influence of energy resource allocation on population abundance. Thus, Akhmetzhanov et al. (2011) 10 analyzed the allocation energy strategies for reproduction and resource foraging, using a mathematical approach with a semi-discrete model 32 . The long-term dynamics that emerge from the proposed model include the extinction of the population due to an over-exploitation of resources, coexistence in a season-to-season equilibrium, and oscillations in population size and resources. Furthermore, they concluded that the optimal allocation strategies have a stabilizing role in the consumer-resource interaction. However, no link has been established between allocation strategies towards different life-history functions and long-term population abundance.
Considering that the emergence of dynamic patterns are linked to how individuals allocate their resources to different life-history functions, our objective is to analyze the consumer-resource relationship and its effects on the persistence of populations living in seasonal environments, incorporating a constant strategy of energy allocation towards both reproduction and maintenance. To address our objective, we use semi-discrete models, due to the division of the annual cycle into two seasons (namely, reproductive and non-reproductive) acting on two different time scales (namely, continuous and discrete). These models have a common mathematical formalism in terms of impulsive differential equations 33,34 , which are widely used to address topics of interest in epidemiology and population ecology [35][36][37][38][39] . Our consumer-resource model has a "bottom-up" mechanistic formulation [40][41][42] to obtain a better understanding of the dynamic behaviors of the populations by incorporating the individual allocation of energetic resources towards reproduction.

Consumer-resource and energetic model
Our consumer-resource and energetic model considers a population of individuals with generational overlap, in which consumers feed exclusively on a resource that is completely restored every cycle. In addition, individuals experience intraspecific competition for resources. The annual cycle is divided into two seasons: reproductive and non-reproductive. Throughout each non-reproductive season, individuals consume the resource at a rate dependent on the equal division of the resource (ratio-dependent), which is allocated between reproduction and maintenance.
In each temporal cycle of length τ > 0 , the reproductive season represents a relatively small period 43 , which is symbolized by instants t n = nτ , n ≥ 0 . Thus, the development of non-reproductive seasons are taken to be in the intervals (nτ , (n + 1)τ ].
Dynamics in the non-reproductive season. Let P(t) be a measure of the population abundance at time t ≥ 0 , which decreases by natural deaths in each non-reproductive season, according to the Malthus law with rate > 0 . Individuals are exclusive consumers of a resource whose density is denoted by R(t). The individual consumption (per unit time) C increases with respect to the resource per capita rate R(t) /P(t) 44 www.nature.com/scientificreports/ where r 0 is the half-saturation constant, namely, the per capita resource density for which an individual consumption rate, equal to R max /2 , is obtained. Therefore, the resource consumed by the population changes at a rate that is described by The rate (1) is the ratio-dependent functional response. This rate is a Holling type II functional response with the resource per capita rate R(t) /P(t) as variable 44 . The use of the ratio-dependent functional response in predator-prey models provides the simulation of more realistic scenarios with a greater dynamic richness. In addition, it corrects the assumption that predation rates are unsaturated and invalidates the paradoxes of biological control and enrichment 29,46,47 . The resources consumed by each individual are an energy source that is assigned to reproduction and maintenance in the fraction α ( 0 < α < 1 ) and 1 − α , respectively. The energy invested in reproduction and maintenance are quantified by E r (t) and E m (t) at time t ≥ 0 , respectively.
During the non-reproductive season, the individuals forage, obtaining the energy necessary for their maintenance. In foraging, the associated costs can be measured in terms of time and energy spent, without forgetting that vital processes have intrinsic costs 4 . In our mathematical approach, these costs (per unit time) are described by the fraction of the resource (such as energy) which is invested in maintenance: Thus, the energy reproduction and maintenance rates are determined by where ξ c = 1 − ξ . Importantly, E m corresponds to the energy that will be principally used by the organism to maintain its offspring in each reproductive event, not to maintain itself. Indeed, for the sake of simplicity, the mortality rate is taken to be independent of E m . Figure 1, illustrates the relationships among the rates (1), (2) and (3).

Dynamics in the reproductive season.
Reproduction is an essential process in the life of organisms. This is evidenced by the evolution of the different life-histories leading to the maximization of reproductive success. However, since each organism has a finite number of resources, when more energy is invested in reproduction, then less energy is invested in other functions, implying a compromise between the future reproduction and survival of organisms 48 . According to Stearns (1989), the most important trade-off in the life-history of organisms usually involves reproduction costs, composed of both fecundity and survival costs. The fecundity, defined as the potential maximum of the physiological reproductive performance of an individual throughout its useful life, is a concept widely studied in population ecology 49 . The current reproductive performance of organisms is described by fertility, which can vary spatially and temporally among individuals due to variation of both environmental conditions and demographics 49 . Thus, considering that the consumed resource is constituted of energy for individuals, we assume a fertility B proportional to the reproductive energy, namely B (E r (nτ )) = γ E r (nτ ) , where γ is the number of offspring per unit of reproductive energy. According to Reznick (1992), the cost of reproduction is represented by an inverse relation between high fertility and both survival or future reproduction of organisms, which in energetic terms, could imply a reduction in the energy available for maintenance and reproduction due to payment for the costs of survival and fecundity. To represent the fecundity and survival costs, we will consider both increasing and saturated expressions with respect to the increasing fertility, described respectively by C f = κB with κ = κ(γ ) > 0 such that 0 < κγ ≤ 1 and , (2) www.nature.com/scientificreports/ Note that κγ is a non-dimensional constant, which corresponds to the fraction of reproductive energy used in the reproductive events, denominated reproductive effort. Consequently, the fertility costs C f are a fraction of the reproduction energy E r (nτ ). The per capita growth rate r is defined to be the product of the fertility B and the probability of surviving the reproductive season S 15 , which is dependent on the energy allocated to reproduction and maintenance (after the application of costs), respectively. Then, with S (x) = x/(e 1/2 + x) where e 1/2 is the energy at which the survival probability reaches one-half 50,51 . In addition, we assume that a fraction µ ( 0 < µ < 1 ) of the individuals that reproduce die during the reproductive season. Then, the population abundance post-reproductive season is Finally, at the end of each reproductive season, the resource density is restored to the value R(nτ + ) = K.
The model. The dynamics are modeled by a consumer-resource model which divides the annual cycle into two seasons and posits a constant allocation of the resource, as an energy source, towards life-history functions. The relations between energy states, resource density, and population abundance allow defining the following impulse differential system:

Results
The long-term population dynamics (i.e., extinction or persistence in an equilibrium value) of the system (2) is mainly dependent on the constant allocation of energy between life-history functions. In the case of population persistence, there is a unique long-term allocation strategy towards reproductive functions that maximizes the population abundance and minimizes individual consumption. In addition, this strategy is dependent on the parameters associated with both fertility ( κ and B 0 ) and survival (ξ c and e 1/2 ) costs.
Preliminary results. To investigate long-term dynamic patterns, we will relate the energy states of reproduction and maintenance, resource density and population abundance at the end of reproductive seasons, namely in the time sequence {nτ + } n≥0 . In the system (4) the relationships between the state variables are the same for each τ unit of time. Then, there is a transformation that relates to the vector (R, E r , E m , P)((n + 1)τ + ) with (R, E r , E m , P)(nτ + ) such that (n + 1)τ + − nτ + = τ . This relationship is determined by the following discretization (or stroboscopic map) of the impulsive differential system (4): (4) X η : www.nature.com/scientificreports/ and R(nτ + ) = K , where the function defined by (5) is evaluated at t = (n + 1)τ and extended to p = 0: for any n ≥ 0 . Indeed, in non-reproductive seasons the consumer-resource dynamics are described by the continuous component of the system (4), which can be solved. Directly we have P(t) = P(nτ + )e − (t−nτ ) for any t ∈ (nτ , (n + 1)τ ] . In addition, the per capita consumption rate (1) can be expressed in terms of R ′ (t) . Therefore, reproductive and maintenance energy rates are determined by Integrating these functions in the interval (nτ , t] , we obtain Evaluating Eqs. (8)-(9) at t = (n + 1)τ (end of the non-reproductive season), we have with �(n, P(nτ + )) = �((n + 1)τ , R((n + 1)τ ), P(nτ + )) and P((n + 1)τ ) = P(nτ + )e − τ for any n ≥ 0 . In addition, evaluating the discrete component of system (4) at t = (n + 1)τ , we obtain Finally, substituting the Eqs. (10)-(11) into the Esq. (12)- (14) we obtain the discretization given by Eq. (6). In order to obtain the equilibrium points of system (6), we can solve the following equations: w h e r e e r := lim n→+∞ E r (nτ + ) , e m := lim n→+∞ E m (nτ + ) , ρ := lim n→+∞ P(nτ + ) a n d t h e n , �(ρ) := lim n→+∞ �(n, P(nτ + )) . From the first equation, we have e r = (1 − κγ )α�(P)/κγ and then e m = κB 0 (1 − α)ξ c /α for P ∈ {0, ρ} such that �(0) = R max (e τ − 1)/ and Therefore, assuming = 0 , we have Ŵ(ρ) = K − A ρ (from equations (7) and (15)) where Ŵ := lim n→+∞ R((n + 1)τ ) , ρ is the solution of E m ((n + 1)τ + ) e 1/2 + E m ((n + 1)τ + ) P((n + 1)τ ). Long-term population dynamics. From the discretization (6), there are two dynamic behaviors for the long-term population abundance: extinction (see Fig. 2a) and persistence (see Fig. 2b). The differentiation of these behaviors strongly depends on the individual consumption of resources defined throughout each non-reproductive season by At the end of each non-reproductive season, namely at the time t = (n + 1)τ , the individual consumption is given by where R((n + 1)τ ) is the non-consumed resource density by the population, the amount that is obtained from the implicit solution of the resource density equation in the continuous component of the system (4). Thus, projecting the individual consumption of the resource into the long term, and taking Ŵ(ρ) = (K − A ρ + I )e − τ , the expression (17) assumes the form where I := lim n→+∞ (n+1)τ nτ R(s)e (s−nτ ) ds . On the one hand, the individual consumption C ∞ I is composed of a basis amount corresponding to the term A and an amount resulting from the equal division of a resource not consumed by individuals dying during the non-reproductive season, K(e τ − 1) − I . Furthermore, when the population experiences a reduced mortality during the non-reproductive season (i.e., ≈ 0) , the individual consumption is C ∞ I (ρ) ≈ A . On the other hand, whether the population abundance is low, the per capita resource is high, implying an individual consumption close to R max (e τ − 1)/ (equivalent to taking the limit of C I ((n + 1)τ ) as P(nτ + ) → 0 ). Certainly, this quantity does not represent the effective individual consumption, but rather establishes an upper limit for this and therefore represents a value of non-persistence. Thus, behaviors related to the equilibrium solutions of the discrete system (6) can be differentiated by the threshold value , t ∈ (nτ , (n + 1)τ ].
(17) C I ((n + 1)τ ) = [K − R((n + 1)τ )]e τ P(nτ + ) , . www.nature.com/scientificreports/ In particular, when the mortality of the population is low, the long-term abundance is described by Eq. (16) and the threshold value assumes the following form U = R max τ/A . Thus, we conclude that the persistence of the population is established when U > 1 and extinction when U ≤ 1.
Finally, we can see that the stabilization of population size in the long term is in response to a dense-dependent behaviour where the per capita growth rate in the long term is r ∞ := {αA /κ}S (e m ) equivalent to mortality fraction µ , where e m = κB 0 (1 − α)ξ c /α is the equilibrium value of energy maintenance after the reproductive season. In addition, the derivative of r with respect to P is where d�/dP < 0 , dr/dP < 0 and dr/d� > 0 are obtained, which explains the expected dense-dependency.
Then, we can derive the following conclusion.

Theorem 1
We consider the threshold value U = R max τ/A .

1.
If U ≤ 1 then the long-term population behavior is extinction.

If U > 1 then the long-term population behavior is persistence.
Proof We divide the proof into two cases: κγ = 1 and 0 < κγ < 1 . In the first case we define the functions Therefore, the discrete system defined by equations e k+1 m = G(e k m , p k ) and p k+1 = H(e k m , p k ) with k ≥ 0 has two equilibrium solutions: E (P) = (E m , P) where E m = B 0 (1 − α)ξ c /αγ and P ∈ {0, ρ} . Thus, with � ′ (p) = −r 0 �(p)/[r 0 p + Ŵ(p)] where ( ) ′ means the derivative with respect to p. Then, the matrix of linearization around the equilibrium point E (P) is given by the Jacobian matrix in E (P), Therefore, E (0) is uniformly asymptotically stable and E (ρ) is locally asymptotically stable.
To analyze the stability of the equilibrium point E := E (ρ) , we consider the characteristic polynomial of the linearization matrix around this point, F (x) = x 3 + ax 2 + bx + c = 0 , where From the sign of the derivatives a, c < 0 and b > 0 are obtained. Thus, the necessary conditions of Jury's criteria: Adding these expressions, ac − b < 1 − c 2 is obtained. Whereas for analyzing the inequality c 2 − 1 < ac − b we consider the expression Note that in Eq. (19), κγ is the amount that provides the greatest variability to this expression. Then, we define the function for any x ∈ [0, 1] . Now, we will show that 0 < g(x) < 1 for any x ∈ [0, 1] . The function g is a fourth degree poly- where the values of B 3 and B 1 depend on the following three disjoint regions (see Fig. 3): Some important properties of the function g are: Fig. 4). Therefore, by the method of Descartes, we have the five admissible cases which are presented in Table 1. On the one hand, in the cases C 1 , C 3 , C 4 and C 5 , g ′ has there is only one positive real root x = x(A 1 , A 2 , A 3 ) where ( ) ′ means the derivative with respect to x. Note that x > 1 , because otherwise there would be at least one additional zero of g ′ on the interval (0, 1) for the basic properties are met. Then, the function g is a positive and decreasing function for any x ∈ [0, 1] with maximum value g(0) = B 0 < 1 . On the other hand, in case C 2 , g ′ has one positive real root x ∈ (0, 1) at which the polynomial g has a local maximum. Taking g(x) to define a function and only if (A 1 , A 2 , A 3 ) ∈ R 3 (using the application Mathematica). Thus, in all cases, 0 < g(x) < 1 for any x ∈ [0, 1] , which implies c 2 − ac + b < 1 and hence c 2 − 1 < ac − b . Therefore, E (ρ) is locally asymptotically stable.
Optimal strategy. The validity of the condition U > 1 depending on α establishes a range of allocation strategies (α m , α M ) in the which the persistence of the population is guaranteed (see Fig. 5). Thus, α ∈ (α m , α M ) if and only if U > 1 , where α m and α M are the limits of this range such that U = 1 is obtained.
In this persistence range we find a unique allocation strategy described by Scientific Reports | (2020) 10:21401 | https://doi.org/10.1038/s41598-020-77326-1 www.nature.com/scientificreports/ which maximizes the population abundance and minimizes the individual consumption in the long term (see Fig. 6). This allocation strategy satisfies ρ ′ (α op ) = 0 if, and only if, also satisfies A ′ (α op ) = 0 where ( ) ′ means the derivative with respect to α and The arrow on the axes x = 0 and x = 1 represent the sign of g ′ (0) and g ′ (1) respectively. www.nature.com/scientificreports/ In addition, the strategy α op is dependent on the parameters associated with both the costs of fertility ( κ and B 0 ) and survival ( ξ c and e 1/2 ), so that if the fertility costs � := κB 0 are greater than the maintenance ratio δ := e 1/2 /ξ c , the allocation is to favor reproduction, i.e., α op > 0.5 (see Fig. 6a). Otherwise, the allocation is to favor maintenance, i.e., 1 − α op ≥ 0.5 (see Fig. 6b).

Discussion
By studying the dependence of the population abundance in seasonal environments on strategies of energy allocation to life-history functions, we have proposed a consumer-resource model built from the division of the annual cycle (into a reproductive season and a non-reproductive season) and the investment of the resource, as an incoming energy flow to reproduction and maintenance.
In stable environments that have a constant amount of resources (or that exhibit only minimal variations) between consecutive annual cycles, long-term persistence is a consequence of negative feedback between the increase in population size and available resources 52 . Thus, our model showed that long-term trends of population size (i.e., extinction and persistence) depends on the per capita consumption of the resource. Such behaviors are differentiated through a threshold value that admits of a representation as the per capita consumption of the resources in the long term.
The differentiation of dynamic behaviors that emerge from biological models is often done through an amount that acts as a threshold value (e.g., 53,54 ). Specifically, in ecological populations with age structure, population dynamics are determined by the net reproductive rate R 0 5,49 that determines the population growth and the Table 1. Admissible cases given by the method of Descartes applied to the polynomial g assuming that its coefficients B i with i ∈ {0, 1, 2, 3, 4} take the values according to (A 1 , A 2 , A 3 ) ∈ (0, 1) 3 . Each row corresponds to the case in which the signs of the coefficients of g are as indicated depending on the region R j with j ∈ {1, 2, 3}. Figure 5. The individual consumption in the long term versus allocation strategy towards reproduction. Note that the size of the persistence range reduces with decreasing length of the cycle τ or of the maximum rate of consumption R max , which is indicated by the respective arrows. In addition, A = A (α) is a convex function, which implies that there is a unique allocation strategy, into the persistence range, that minimizes individual consumption in the long term. www.nature.com/scientificreports/ potential self-replacement when R 0 ≥ 1 . On the contrary, when R 0 < 1 , the long-term dynamics results in extinction.
Our approach assumes that the assigned resource is an energy source for individuals, invested in life-history functions (e.g., reproduction and maintenance) that determine the population growth 52,55 . This corresponds with the outcome of interconnections between food, allocation and life-history 56 that through our model allows us to find classical behaviors for the size of the population, which correspond to extinction and persistence, and in response to a density-dependent mechanism. A wide literature has shown there is an increase (generally, with a maximum value) of the per capita population growth rate (r) when there is an increase in food availability at both the population and individual levels 52,55,[57][58][59][60][61][62] . Additionally, it shows that r is negatively related to the population density (D), (e.g., r ∝ D −1 ) which allows concluding that there is a negative relation between the per capita availability of resources (F) and D (e.g., F ∝ D −1 ).
Our results establish that the individual scale is important for population persistence. Individual consumption C ∞ I (ρ) is the quantity that collects the main ecological parameters of the life-history assumed in this model. The parameter that causes the greatest variability in C ∞ I (ρ) = A determines the allocation of resources towards reproduction. For this parameter, we have found the range of values that allow long-term population persistence. There is a unique strategy α op which implies a maximum population size ρ(α op ) and a minimum consumption per individual A (α op ) in the long term, inversely related quantities since ρ(α op ) · A (α op ) ∝ K . This optimal strategy depends on the parameters associated with fertility and survival costs, α op = α op (κ, B 0 , e 1/2 , ξ c ) . Although our modeling approach considers allocation strategies represented by a parameter and not by a time-dependent variable state or an environmental measure, the optimal allocation strategy α op shows variability, being plastic, either when the fertility cost and/or the maintenance ratio varies. Indeed, when the maintenance ratio decreases ( δ = e 1/2 /ξ c → 0 ) the individuals increase their survival due to (i) a decrease in the maintenance cost factor ξ ending the non-reproductive season with higher maintenance energy E m (nτ ) or (ii) a decrease in the amount of e 1/2 which implies that only a little maintenance energy E m (nτ + ) is needed to guarantee the individual's survival. In this case, the optimal allocation strategy favors the reproductive process, as Fischer et al. (2009Fischer et al. ( ,2010 50,51 point out. Conversely, the optimal allocation strategy favors the maintenance process. In this study, constant allocation strategies establish the persistence and stability of the population when the threshold value is U > 1 . This result differs from those found by Akhmetzhanov et al. (2011), who established that the strategies of constant allocation lead to population extinction due to a gradual loss or over-exploitation of the resource. This difference may be related to the different mathematical representations used to describe the dynamics in the reproductive and non-reproductive seasons, where we highlight the functional response chosen to describe the predator-prey relationship (ratio-dependent vs. linear prey-dependent), and the variability in the resource density available for consumption at the beginning of each non-reproductive season (constant resource vs. cumulative resource). We propose that these and other differences should be studied in a general modeling framework that incorporates resource abundance, internal energy, and consumer abundance as state variables, namely, a consumer-resource and energetic mathematical model that incorporates the individual-level in consumer-resource interactions.
We showed that the consumer population's ecological dynamics is limited both by the quality or quantity of food and by the design of the acquisition and energy allocation process. Consequently, the usefulness of our model lies in studying the occurrence of the potential behavior outcomes in the long term for a consumer population as a function of parameters of ecological significance. In the proposed mathematical formulation, the www.nature.com/scientificreports/ reproductive season's ecological dynamics are relevant, allowing extrapolating some of the conclusions obtained. The main reproductive trade-off related to energy allocation is the offspring size and offspring mass. The results of Veloso and Bozinovic (2000) 63 show that with the same reproduction effort, it is possible to obtain numerous small-size offspring or a few large-size offspring, due to the positive relations between reproductive effort, the organism's body mass, and the resting metabolic rate. From our model, an increased reproductive effort leads to a lower reproduction energy available after the reproductive season, which implies a lower long-term population size. Due to the negative relation between individual consumption and population size, it is possible to conjecture that the optimal life-history strategy results in a small number of large-size offspring. The organisms have a finite amount of resources that are invested in multiple tasks. When the availability of resources is compromised, organisms often respond by making physiological adjustments (e.g., the resting metabolic rate), allowing the allocation of resources to other vital processes 63,64 . Peña et al. (2020) has shown that stochastic food deprivation produces a compromise between energy expenditure rates and locomotion activity performed by individuals, describing it as a paradoxical result. However, we can hypothesize that a lower individual consumption together with a higher maintenance cost ( ξ → 1 ) implies that this compensation is mediated by an optimal energy allocation, favoring naturally the maintenance activities, sustaining the locomotor (or exploratory) behavior. The model proposed in the present paper has made several assumptions that make it both generally and mathematically accessible. First, all individuals in the population have common traits, such as the allocation strategy towards reproduction. Our approach is based on the theory of life-history that is shared and incorporated by numerous investigations 24 . Therefore, an interesting extension of our model arises precisely by incorporating the variability of traits through phenotypic distributions. Second, energetic resources are assigned to reproductive and maintenance processes, with the latter concentrating on the survival and/or growth of individuals during the annual cycle. Specifically, body size as a measure of an individual's growth has been related to life-history parameters 28,65,66 and birth rates and mortality 3 . We conjecture that incorporating body size would allow us to analyze the consumer-resource dynamics in animals that exhibit determined or indeterminate growth 9,31 . Third, for mathematical simplicity, a constant population size was assumed over the non-reproductive season for obtaining an autonomous system, i.e., independence of the temporal variable. The reduction in population size during the non-reproductive season, by assuming > 0 , implies an increase in per capita resource consumption, which reduces the set of strategies that establish the persistence of the population. Therefore, in the general case, the size of the persistence range L = α M − α m depends on the population size. However, the long-term dynamic behaviors remain unchanged. Fourth, at the end of each reproductive season, the individuals that enter the population do not differ by age from the surviving individuals (except in the case µ = 1 ). Such differentiation would permit an advance towards a reformulation of our model with an age structure, incorporating physiological differences between individuals of different groups. These points pose new challenges for future research.