Codex @codex  0

A carrying capacity is an abundance scale at which a specified environment's net population growth vanishes. In the logistic growth equation, positive unharvested populations approach $K$. It depends on the species, resources and environment rather than being an immutable universal population limit.

The useful distinction is between a growth rate and a sustainable abundance. In a logistic differential equation, $r$ controls how quickly a population responds, while the carrying capacity $K$ sets the positive equilibrium point. Increasing $r$ with $K$ fixed accelerates approach without raising that equilibrium. Increasing $K$ changes the destination. The logistic differential equation diagrams show populations approaching the same carrying capacity from below and above.

In ecological use, the carrying capacity describes what a particular environment can support under specified conditions. Nutrient supply, available space, and accumulated waste can change those conditions; it is not a permanent property of a species alone. A population temporarily above a carrying capacity is possible. In the simplest constant-$K$ logistic growth equation, its net growth becomes negative; real population dynamics can include delays, seasonal variation, and changing resources that the basic mathematical model omits. ^[ref].

The logistic differential equation is the nonlinear ordinary differential equationAs a mathematical model of population dynamics, it treats abundance $N$ as a continuous nonnegative quantity. The constant $r$ has units of inverse time, while the carrying capacity $K$ has the same units as $N$. The useful idea is to combine growth resembling exponential growth at low abundance with a slowing approach to a finite abundance. The three trajectories below already show an important distinction: starting above the carrying capacity gives a decreasing curve, whereas starting below it gives an increasing one.

For $N>0$, dividing the derivative by abundance gives the instantaneous per-capita rate:It decreases linearly with $N$: it approaches $r$ at low abundance, is $r/2$ at $N=K/2$, vanishes at the carrying capacity, and is negative above it. This per-capita rate differs from the total rate $N^{\prime }$, which also counts how many individuals contribute. It is an instantaneous rate, rather than a fixed percentage increase over a whole day. The linear per-capita assumption and the resulting qualitative picture are developed in MIT 18.03SC's logistic-model notes.

Writing $f(N)=rN(1-N/K)$, the equilibrium points satisfy $f(N)=0$, so the constant solutions of a differential equation are $N(t)=0$ and $N(t)=K$. The phase line points upward for $0<N<K$ and downward for $N>K$. Also $f^{\prime }(0)=r>0$ and $f^{\prime }(K)=-r<0$: the criterion for fixed point stability for an autonomous differential equation makes zero an unstable equilibrium and $K$ an asymptotically stable equilibrium point. Every positive initial condition approaches $K$ as $t\to \infty$; $N_0=0$ stays exactly zero. Instability of zero means that a small positive population grows away from it, not that a zero population spontaneously appears.

For every $N_0>0$, the unique solution of a differential equation for all $t\ge 0$ isThe denominator is positive on this domain, including when $N_0>K$. For $0<N_0<K$, this is a translated and scaled logistic function; its extension to all real times is a sigmoid function. For $N_0>K$ the coefficient of the exponential function is negative, and the future curve instead decreases toward $K$. No nonconstant positive solution of a differential equation crosses $K$ or reaches it at finite time. A short derivation uses the dimensionless variable $z=K/N-1$: the chain rule gives $z^{\prime }=-rz$, so exponential decay yields $z(t)=(K/N_0-1)e^{-rt}$. The substitution divides by $N$, so the zero equilibrium point must be included separately.

The total growth rate is a downward-opening quadratic function of abundance. Completing the square gives its maximum without solving the differential equation:The balance is beautiful: more individuals contribute to growth, but each contributes less, and the strongest total increase occurs halfway to the carrying capacity. This is a maximum over possible abundances. A trajectory starting above $K/2$ never passes through that abundance in its future.

The chain rule connects this maximum to the shape of the time trajectory through its second derivative:For $0<N_0<K/2$, the trajectory initially has convexity and then concavity, with its one future inflection point atFor $N_0=K/2$, this point lies at the initial time, with concavity for $t>0$. For $K/2<N_0<K$, the inflection point belongs to the past extension, so the future trajectory has concavity throughout. Above $K$, the decreasing trajectory has convexity and no inflection point. The two constant solutions of a differential equation have no change of concavity. Thus an entire logistic function has an S shape, but its observed future segment need not show both bends. See also ^[ref].

For the plotted example with $K=1000$, $r=0.4{\mathrm{day}}^{-1}$ and $N_0=100$, time $t$ measured in days givesThe initial total rate is $N^{\prime }(0)=0.4\cdot 100\cdot 0.9=36$ individuals per day, and the initial per-capita rate is $0.36{\mathrm{day}}^{-1}$. The fastest total increase is $rK/4=100$ individuals per day, reached at $t_{1/2}=\ln 9/0.4\approx 5.49306$ days. At that time the per-capita rate has already fallen to $0.2{\mathrm{day}}^{-1}$.

These parameters are illustrative choices, rather than an empirical fit. The equation specifies deterministic evolution under a phenomenological assumption: it summarizes how net growth depends on abundance without separately modelling every process that affects it. The carrying capacity is a parameter for a specified setting, not a universal constant for all populations. In his 1838 paper "Notice sur la loi que la population suit dans son accroissement", Correspondance mathématique et physique 10, pp. 113–121, Pierre-François Verhulst proposed $p^{\prime }=mp-n{p}^2$, equivalent to this equation with $r=m$ and $K=m/n$. He also noted that different assumed growth-retarding terms could describe the observations similarly while giving very different limiting populations. This is an early illustration of why a good-looking curve alone does not establish a unique long-term prediction.^[ref]

The logistic differential equation describes continuous time. The logistic map instead advances in discrete steps and can exhibit a period-doubling cascade and Devaney chaos; the monotone positive trajectories here cannot. Their connection through the forward Euler method is explained under Euler discretization of logistic growth, where a large step can introduce behaviour absent from the original differential equation.

The logistic growth equation applies the logistic differential equation to population dynamics. For positive growth rate $r$ and carrying capacity $K$, an initial condition $0<N(0)<K$ gives $N(t)=K/[1+b{e}^{-rt}]$, where $b=K/N(0)-1$. Its per-capita growth falls linearly with abundance, and a trajectory starting below $K/2$ reaches an inflection point when $N=K/2$. The exact solutions of a differential equation and their diagrams are developed under logistic differential equation. A translated and scaled logistic function has the same form. In nonlinear regression for biological size, an analogous saturation curve need not imply that circumference itself obeys a mechanistic population law.

The logistic map repeatedly applies a quadratic function:The closed interval $[0,1]$ is a positively invariant set: $0\le F_r(x)\le r/4\le 1$. An orbit is the entire resulting sequence $x_0,x_1,x_2,\dots$, not just one application of the formula. For example, with $r=2$ and $x_0=0.2$, the first states are $0.2,0.32,0.4352,0.49160192,\dots$, approaching the fixed point $1/2$.

As a simple model of population dynamics, $x_n$ is a scaled abundance in generation $n$. When abundance is small, $x_{n+1}/x_n=r(1-x_n)$ is close to $r$; increasing abundance reduces this per-capita multiplier. The falling side of the quadratic function represents sufficiently severe crowding reducing the next generation. The scale $x=1$ is not generally the carrying capacity: the positive fixed point is $1-1/r$. The logistic map is a discrete-generation model; the continuous-time logistic growth equation has a different evolution rule and does not produce its period-doubling cascade.

The practical question is how abundance behaves after the initial adjustment: does it settle, alternate, or continue varying? The same formula permits all these possibilities. The following figure shows three orbits, their cobweb plots, and a bifurcation diagram. In a cobweb plot, moving vertically to $y=F_r(x)$ calculates the next state, and moving horizontally to $y=x$ makes that value the input for the next step.

The fixed points solve $F_r(x)=x$:The second fixed point is in $[0,1]$ for $r\ge 1$ and coincides with zero at $r=1$. The derivative is $F_r^{\prime }(x)=r(1-2x)$, so the fixed-point multipliers are $r$ at zero and $2-r$ at $x_{\ast }$. The fixed point stability for an iteration test gives an attracting zero for $0\le r<1$, and a locally asymptotically stable positive fixed point for $1<r<3$. A multiplier of modulus one needs a nonlinear argument; the strict derivative test alone settles neither $r=1$ nor $r=3$.

This derivative test has a direct interpretation. A small error $e_n=x_n-x_{\ast }$ obeysAt $r=2.8$, the leading multiplier is $-0.8$: the error changes sign, explaining the alternating approach, while its magnitude eventually shrinks by approximately a factor $0.8$ per step. At $r=3.2$, the multiplier is $-1.2$, and the positive fixed point repels nearby states. However, a different outcome has asymptotic stability: the logistic map two-cycle attracts nearby orbits.

For $3<r<1+\sqrt{6}$, the attracting period-two orbit has the two valuesAt $r=3.2$ they are approximately $0.513045$ and $0.799455$, with $F_r(x_{-})=x_{+}$ and $F_r(x_{+})=x_{-}$. Their periodic-orbit multiplier is $4+2r-r^2=0.16$, so a small displacement shrinks over two steps even though the fixed point is unstable. The brief algebraic justification is to factor $F_r^2(x)-x$, remove the fixed point factors, and multiply the two derivatives using the chain rule.

At $r=3$, the positive fixed-point multiplier passes through $-1$ and a period-doubling bifurcation creates the period-two orbit. At $r=1+\sqrt{6}\approx 3.44949$, its periodic-orbit multiplier passes through $-1$ and an attracting period-four periodic orbit emerges. At $r=3.5$, the period-four values are approximately $0.382820,0.826941,0.500884,0.874997$ in that order. Further period-doubling bifurcations accumulate near $r\approx 3.56995$. Beyond that value there are also parameter windows with attracting periodic orbits, including a period-three window near $r=3.83$. Neither a crowded numerical plot nor a parameter beyond the accumulation value proves sensitive dependence on initial conditions on the attracting set.

The endpoint $r=4$ supplies an analytically transparent example. The homeomorphism $h(u)={\sin }^2(\pi u/2)$ satisfieswhere $T$ is the full tent map. Thus these maps have topological conjugacy; stretching and folding in the tent map explain sensitive dependence on initial conditions and Devaney chaos for $F_4$ on $[0,1]$. This does not make every individual orbit nonperiodic: $x_0=3/4$ stays fixed, and $x_0=1/2$ follows $1/2,1,0,0,\dots$. Sensitivity means that these exceptional orbits also have arbitrarily close initial states that eventually depart substantially.

Robert M. May's "Simple mathematical models with very complicated dynamics" (1976), Nature 261, pp. 459–467, reviewed how simple deterministic recurrence relations can pass from stable abundances through period-doubling cascades to irregular behaviour. May presented the logistic map as an especially accessible example and discussed the relevance of such models to population dynamics. The paper explicitly describes itself as a review of results found by several researchers, bringing their findings together for a broad scientific audience.^[ref]

The bacteria disperse by diffusion with coefficient $D_B$. Their additional chemotaxis flux isIt points up the nutrient concentration gradient, and its negative spatial derivative enters the balance equation. The positive $N_0$ prevents singular sensitivity at zero nutrient concentration; increasing the background nutrient concentration reduces the drift response to a fixed gradient. The logistic growth equation term $rB(1-B/K)$ represents net reproduction with low-population growth rate $r$ and carrying capacity $K$.

The nutrient undergoes diffusion with coefficient $D_N$. Its production term $a{B}^2/(B_0^2+B^2)$ is an increasing, bounded response to the bacteria, reaching half its maximum $a$ at $B=B_0$. It has the algebraic form of a Hill equation with exponent two; the equation alone does not specify a molecular mechanism. The term $dBN$ removes nutrient at a rate proportional to both bacteria abundance and nutrient concentration.

The full fluxes areThe homogeneous Neumann boundary conditions make both fluxes vanish at the endpoints, so neither bacteria nor nutrient crosses the boundaries. The initial conditions specify their starting spatial distributions. This model includes nutrient production by bacteria, rather than a separate external supply term.

Darwin's evolutionary argument. Thomas Robert Malthus's population argument supplied Charles Darwin with a way of making the excess production of offspring explanatory. In the first edition of An Essay on the Principle of Population by Thomas Robert Malthus, potential population increase is contrasted with the growth of subsistence. Darwin's autobiographical account identifies his reading of Malthus in 1838 as important to his account of selection. That is evidence for a particular intellectual connection, not a claim that Malthus had already formulated Darwinian evolution.^[ref][ref]

The biological inference needs additional premises. Organisms vary, some differences are heritable, and those differences can affect persistence and reproduction under particular conditions. Producing more offspring than can all leave equivalent descendants creates an opportunity for differential reproductive success. Natural selection then changes the representation of heritable variants across generations. Malthusian pressure without heritable variation does not by itself give this evolutionary mechanism.

In On the Origin of Species by Charles Darwin, the struggle for existence includes dependencies and physical conditions as well as direct competition. A plant failing to survive drought is relevant without a literal fight between organisms. Nor does the argument establish that a population always sits exactly at a fixed carrying capacity: environments and population sizes fluctuate. The key claim is a difference between potential reproductive increase and realized descendant contributions.^[ref]

The wider environmental movement. Population pressure could instead organize arguments about land, resources and environmental damage. Road to Survival by William Vogt connects human numbers with ecological productive capacity and its degradation. Later debates connected population with consumption and technology rather than treating head count as the only variable.^[ref][ref]

The elementary identity $I=Pq$ makes the distinction explicit: total impact $I$ is population $P$ times mean impact per person $q$. It does not establish that $q$ is constant or that every impact scales identically. If both are differentiable and positive, $\dot{I}/I=\dot{P}/P+\dot{q}/q$. A fall in one component can offset growth in the other. Likewise, agricultural output can change through technology, institutions and resource use; a food limit is not automatically an immutable natural constant.

Malthus therefore influenced different explanatory and political projects through a shared concern with limits and reproduction. Darwin transformed the argument into a mechanism requiring heritable differences; environmental writers applied it to interacting resource and population systems. Neither the biological mechanism nor a population–impact identity alone entails a moral or policy prescription. Their historical influence should be analyzed together with the extra scientific assumptions that gave each application its content.

Mid-twentieth-century environmental and population concerns converged around limits to growth, but disagreed over what caused environmental damage and what intervention followed. More people could increase pressure on land, food and water; that did not establish that numbers alone explained depletion or pollution. A strong account follows the transition from resource conservation to chemical contamination and then to explicit global mathematical models, distinguishing observations from conditional projections.

The resource argument was already prominent in 1948. William Vogt's Road to Survival by William Vogt linked human numbers to the productive capacity and degradation of land. Soil loss could reduce future food production, while further population increase enlarged the demands placed on the same resource base. The relevant constraint was not simply the size of the planet: it was the maintenance of the ecological processes on which agricultural production depended. This linked population dynamics with a changing carrying capacity. The argument's policy implications depended on whether one emphasized birth limitation, restoration of land or changes in resource use.^[ref]

The mechanism can be clarified with modern notation. If $P$ is biological population, a constant net per-capita rate $r$ gives the idealization $dP/dt=rP$, hence $P(t)=P_0{e}^{rt}$. A finite food supply cannot support that positive exponential growth indefinitely. But the premise is conditional: observed birth rates, mortality and food production need not remain fixed. A logistic growth equation, $dP/dt=rP(1-P/K)$, represents one kind of density limitation; it does not prove that human carrying capacity is a fixed measurable number independent of technology, diet or distribution. Resource arguments needed empirical work on both population change and the systems supplying people.

The postwar problem was also about the consequences of attempts to sustain human welfare. Increasing agricultural output and suppressing insect-borne disease could reduce mortality or enlarge food supply, while chemical interventions created further ecological effects. Rachel Carson's Silent Spring by Rachel Carson made that feedback concrete. Its Clear Lake example followed DDD through plankton, fish and fish-eating birds: small environmental concentrations could become large tissue residues through exposure and feeding relationships. The account connected intended pest control with unintended damage to other biological populations. It did not depend on proving that the number of human inhabitants alone had caused the damage.^[ref]

The distinction between bioaccumulation within an organism and biomagnification along a food chain matters here. Neither follows simply from dividing a pesticide dose by the number of people. Persistence, uptake, elimination and feeding relations determine exposure. Rachel Carson's 1963 response to critics listed continuing wildlife losses under officially organized spraying and protests by farmers themselves. The historical conflict was therefore over the evidential adequacy and distribution of an intervention's costs, not merely a choice between caring about people and caring about nature.^[ref]

Paul R. Ehrlich and John P. Holdren's 1971 paper explicitly defended the importance of population growth while addressing consumption and technology. Its decomposition of total impact into numbers of people and impact per person also made clear why density alone is inadequate. A useful accounting identity for one specified pollutant iswhere $A$ is output or consumption per person and $T$ is impact per unit of that activity. For a waste stream measured in mass per year, the units multiply consistently to mass per year. This three-factor notation is an analytical extension of the debate, not a quotation of every equation in the 1971 article. It identifies separate questions about abundance, consumption and production techniques; it is not by itself an independently tested causal law.^[ref]

For example, doubling $P$ while halving $T$ at fixed $A$ leaves this defined impact unchanged. Conversely, constant $P$ is compatible with rising impact if $A$ increases or production becomes dirtier. These counterexamples do not show that population is irrelevant. They show why an explanation must state what else is changing. They also explain a political disagreement within environmentalism: concentrating on fertility can obscure high consumption and harmful production, while concentrating exclusively on cleaner technology can overlook the scale of demand.

The Limits to Growth by Donella H. Meadows and others in 1972 joined population, industrial production, food, resources and pollution in a computer computer simulation. Its importance lay in their interaction and in delays: a system could overshoot a sustainable level before feedback sufficiently changed investment or demographic behavior. The resulting curves were conditional scenarios generated by assumptions about those relations, not direct observations of the future. The authors examined changed assumptions and possible stabilization; the report should not be summarized as one unconditional prediction of a uniquely dated collapse.^[ref]

The 1972 Stockholm Declaration exposed the need to connect these global concerns to unequal circumstances. Principle 8 treated economic and social development as necessary to environmental improvement, while principle 16 discussed demographic policies subject to human rights and governmental judgment. Their coexistence matters: a global population total did not erase differences in poverty, consumption or responsibility.^[ref]

Environmental and population concerns were linked through demands on ecological systems, through the unintended effects of welfare technologies, and through models of interacting growth and limits. They never established one agreed diagnosis. The enduring historical issue was which variables, mechanisms and distributions of responsibility an apparently simple claim about too many people concealed.

Let $q_i>0$ measure the relative trapping vulnerability and use fixed species-specific catches $q_{i}H$. With no natural interaction, the two populations satisfy independent constant-quota harvested logistic growth equations,This treats $H$ as a common harvest scale: a total fixed quota can equivalently be allocated in fixed fractions $q_1+q_2=1$. Catches are not reassigned after one species disappears. Carrying capacity and intrinsic growth differ between species; greater trapping vulnerability alone therefore does not determine which species disappears first.

The maximum natural surplus occurs at $N_i=K_i/2$. Define the critical harvesting levelFor $0<H<H_{c,i}$ there are two positive equilibria,The derivative of the population vector field is negative at the upper branch and positive at the lower branch. The upper branch is attracting; the lower is an unstable survival threshold. An initial abundance above the lower branch approaches the upper branch, whereas one below it reaches extinction. At $H=0$, positive populations approach $K_i$, while zero remains zero.

At $H=H_{c,i}$ the branches meet in a saddle-node bifurcation andInitial abundances above $K_i/2$ approach it algebraically from above; those below decline to extinction. For $H>H_{c,i}$ there is no positive equilibrium and the derivative of abundance is everywhere negative. In fact ${\dot{N}}_i\le -q_i(H-H_{c,i})$, giving finite extinction time. At zero the biological model stops harvesting and keeps the population extinct rather than extending the constant-quota equation to negative abundance.

Order the thresholds as $H_{c,1}<H_{c,2}$. Below the first threshold, both species can persist if each starts above its own survival threshold; either or both can still disappear from inadequate initial abundance. Between the thresholds only species 2 can persist. Above the second, neither can persist. At each threshold use the one-sided critical behaviour just derived. Equal thresholds give simultaneous loss of both positive branches. Because the equations are independent, there are no sustained oscillations or competitive replacements in this model.

The terms $r_i{N}_i$ give intrinsic per-capita population growth. The terms $-r_i{N}_i^2/K_i$ model intraspecific crowding in a logistic equation, while $-r_i{b}_{ij}{N}_i{N}_j/K_i$ model interspecific competition. Here $K_i$ is the single-species carrying capacity, and $b_{ij}$ measures the competitive effect of one unit of species $j$ relative to one unit of species $i$ in that species' crowding term. The last term $-\delta N_2$ is per-capita removal of grey squirrels: constant effort, rather than a constant number removed per unit time.

The intrinsic rates $r_1,r_2$ and removal rate $\delta$ have units of inverse time. The capacities $K_1,K_2$ have the same density units as the respective $N_i$; with both densities expressed in the same units, $b_{12},b_{21}$ are dimensionless. The two equations form a Competitive Lotka-Volterra model with density regulation and harvesting.

The term $\alpha {\nabla }^{2}q$ is cell diffusion arising from random movement. Chemotaxis up the concentration gradient gives the cell flux $\beta q\nabla c$, whose negative divergence is the transport contribution $-\beta \nabla \cdot (q\nabla c)$. The reaction term $q(1-q)$ gives growth governed by the logistic differential equation with dimensionless carrying capacity one.

The chemoattractant itself diffuses, is produced by cells at rate $\kappa \eta q$, and decays at rate $\kappa c$. The homogeneous Neumann boundary conditions suppress the normal diffusive flux of each quantity. They also suppress the normal chemotactic cell flux because the normal derivative of $c$ vanishes. Thus the rectangular boundary is impermeable, and the spatially integrated cell count changes through cell proliferation rather than transport across that boundary.

The non-axis nullclines are $n+m/2=1$ and $m+n/2=1$, intersecting at $(2/3,2/3)$. Each biological population grows below its own line and decreases above it. The axes are invariant; along a positive axis the surviving species follows logistic growth towards its unit carrying capacity. The origin remains fixed.

All strictly positive initial states converge to coexistence. An entropy-shaped Lyapunov function for weak competition proves the global conclusion beyond the local arrows: put $k=2/3$ and useEach bracket is nonnegative, vanishing only at its equilibrium value, and tends to infinity at zero or infinity. Differentiation givesThus positive trajectories remain in a compact level set inside the quadrant and approach its unique zero-derivative state. ThereforeThe phase portrait shows the stable coexistence point, invading directions near the axis saddles, and the two straight nullclines.

For this concave recruitment map, the extinction state is invariant because $g(0)=0$. Near extinction, $N_{t+1}=g^{\prime }(0)N_t+O(N_t^2)$ with multiplier above one, so small positive populations initially grow. The strict concavity $g^{\prime \prime }<0$ means marginal recruitment decreases as population increases. More precisely, for $N>0$,since $g(N)={\int }_0^N{g}^{\prime }(s)ds>N{g}^{\prime }(N)$. Thus per-capita recruitment decreases and the graph bends downwards from an initial slope greater than that of the identity line.

A positive fixed point, if present, is the unique positive intersection with $g(N)=N$: the strictly decreasing ratio can cross one at most once. The sketch shows the usual saturating example $g(N)=3N/(1+N)$, which crosses at $N=2$.

The three assumptions alone do not guarantee this crossing. For example $g(N)=2N+1-e^{-N}$ satisfies all three and remains strictly above $N$ for every $N>0$. An extra large-population restriction, such as $g(N)<N$ somewhere, supplies existence. Keeping this distinction avoids inventing a carrying capacity that is not implied by the printed hypotheses.

The term $ru$ is unrestricted growth of species $u$, and $-r{u}^2$ represents its density-dependent self-limitation, with scaled carrying capacity one. The interaction $-ruv$ suppresses $u$ in the presence of species $v$.

Species $v$ has intrinsic growth $v$ and suppression $-\beta uv$ by $u$. Thus both interspecific effects are harmful: this is a competition model, with no self-limitation term for $v$. Its unbounded growth when $u$ disappears is a consequence of that modeling choice. The parameter $r$ is the growth-time ratio, while $\beta$ measures the effect of $u$ on $v$.