The population dynamics of pest insects is a subject of interest to farmers, agricultural economists, ecologists, and those concerned with animal welfare.
A life table shows how and how many insects die as they mature from eggs to adults. It helps with pest control by identifying at what life stage pest insects are most vulnerable and how mortality can be increased. [2] A cohort life table tracks organisms through the stages of life, while a static life table shows the distribution of life stages among the population at a single point in time. [3]
Following is an example of a cohort life table based on field data from Vargas and Nishida (1980). [4] The overall mortality rate was 94.8%, but this is probably an underestimate because the study collected the pupae in cups, and these may have protected them from birds, mice, harsh weather, and so on. [4]
Stage () | Number alive at stage () | Cause of mortality | Number dying | Proportion of original cohort surviving () | Proportion of original cohort dying () | Mortality rate at this step () |
---|---|---|---|---|---|---|
Egg | 1000 | Trichogramma confusum wasp parasitism | 386 | 1 | .386 | .386 |
Egg | 614 | Orius insidiosus predation and unknown | 232 | .614 | .232 | .378 |
Larva 1 | 382 | O. insidiosus, cannibalism, and unknown | 90 | .382 | .090 | .236 |
Larva 2 | 292 | Cannibalism and unknown | 59 | .292 | .059 | .202 |
Larva 3 | 233 | Cannibalism and unknown | 34 | .233 | .034 | .146 |
Larva 4 | 199 | Cannibalism, disease, parasitism, and unknown | 59 | .199 | .059 | .296 |
Larva 5 | 140 | Cannibalism, disease, parasitism, and unknown | 69 | .140 | .069 | .493 |
Larva 6 | 71 | Cannibalism, disease, parasitism, and unknown | 1 | .071 | .001 | .014 |
Pupa | 70 | Disease and unknown | 18 | .070 | .018 | 0.257 |
Adult | 52 (28 male, 24 female) | .052 |
From a life table we can calculate life expectancy as follows. Assume the stages are uniformly spaced. The average proportion of organisms alive at stage between beginning and end is [3]
The total number of future stages to be lived by individuals at age and older is [3]
Then the life expectancy at age is [3]
We could have done the same computation with raw numbers of individuals rather than proportions.
If we further know the number of eggs produced (fecundity) at age , we can calculate the eggs produced per surviving individual as
where is the number of individuals alive at that stage.
The basic reproductive rate, [3] also known as the replacement rate of a population, is the ratio of daughters to mothers. If it's greater than 1, the population is increasing. In a stable population the replacement rate should hover close to 1. [2] We can calculate it from life-table data as [3]
This is because each product computes (first-generation parents at age )/(first-generation eggs) times (second-generation eggs produced by age- parents)/(first-generation parents at age ).
If is the initial population size and is the population size after a generation, then [3]
The cohort generation time is the average duration between when a parent is born and when its child is born. If is measured in years, then [3]
If remains relatively stable over generations, we can use it to approximate the intrinsic rate of increase for the population: [3]
This is because
where the approximation follows from the Mercator series. is a change in time, . Then we have
which is the discrete definition of the intrinsic rate of increase.
In general, population growth roughly follows one of these trends: [1]
Insect pest growth rates are heavily influenced by temperature and rainfall, among other variables. Sometimes pest populations grow rapidly and become outbreaks. [5]
Because insects are ectothermic, "temperature is probably the single most important environmental factor influencing insect behavior, distribution, development, survival, and reproduction." [6] As a result, growing degree-days are commonly used to estimate insect development, often relative to a biofix point, [6] i.e., a biological milestone, such as when the insect comes out of pupation in spring. [7] Degree-days can help with pest control. [8]
Yamamura and Kiritani approximated the development rate as [9]
with being the current temperature, being the base temperature for the species, and being a thermal constant for the species. A generation is defined as the duration required for the time-integral of to equal 1. Using linear approximations, the authors estimate that if the temperature increased by (for instance, maybe = 2 °C for climate change by 2100 relative to 1990), then the increase in number of generations per year would be [9]
where is the current annual mean temperature of a location. In particular, the authors suggest that 2 °C warming might lead to, for example, about one extra generation for Lepidoptera, Hemiptera, two extra generations for Diptera, almost three generations for Hymenoptera, and almost five generations for Aphidoidea. These changes in voltinism might happen through biological dispersal and/or natural selection; the authors point to prior examples of each in Japan. [9]
Sunding and Zivin model population growth of insect pests as a geometric Brownian motion (GBM) process. [10] The model is stochastic in order to account for the variability of growth rates as a function of external conditions like weather. In particular, if is the current insect population, is the intrinsic growth rate, and is a variance coefficient, the authors assume that
where is an increment of time, and is an increment of a Wiener process, with being standard-normal distributed. In this model, short-run population changes are dominated by the stochastic term, , but long-run changes are dominated by the trend term, . [10]
After solving this equation, we find that the population at time , , is log-normally distributed:
where is the initial population. [10]
As a case study, the authors consider mevinphos application on leaf lettuce in Salinas Valley, California, for the purpose of controlling aphids. Previous research by other authors found that daily percentage growth of the green peach aphid could be modeled as an increasing linear function of average daily temperature. Combined with the fact that temperature is normally distributed, this agreed with the GBM equations described above, and the authors derived that and . [10] Since the expected population based on the log-normal distribution grows with , this implies an aphid doubling time of days. Note that other literature has found aphid generation times to lie roughly in the range of 4.7 to 5.8 days. [11]
A 2013 study analyzed population dynamics of the smaller tea tortrix, a moth pest that infests tea plantations, especially in Japan. The data consisted of counts of adult moths captured with light traps every 5–6 days at the Kagoshima tea station in Japan from 1961 to 2012. Peak populations were 100 to 4000 times higher than at their lowest levels. A wavelet decomposition showed a clear, relatively stationary annual cycle in the populations, as well as non-stationary punctuations between late April and early October, representing 4–6 outbreaks per year of this multivoltine species. [12] The cycles result from population overshoot. [13]
These moths have stage-structured development life cycles, and a traditional hypothesis suggests that these cycles should be most synchronized across the population in the spring due to the preceding effects of cold winter months, and as the summer progresses, the life stages become more randomly assorted. [12] This is often what's observed in North America. [13] However, this study observed instead that populations were more correlated as the season progressed, perhaps because temperature fluctuations enforced synchrony. The authors found that when temperatures first increased above about 15 °C (59 °F) in the spring, the population dynamics crossed a Hopf bifurcation from stability to repeated outbreak cycles, until stabilization again in the fall. Above the Hopf threshold, population-cycle amplitude increased roughly linearly with temperature. This study affirmed the classic concept of temperature as a "pacemaker of all vital rates." [12]
Understanding life-cycle dynamics is relevant for pest control because some insecticides only work at one or two life stages of the insect. [13]
B. Chaney, a farm advisor in Monterey County, CA, estimates that mevinphos would kill practically all aphids in a field upon application. [10] Wyatt, citing data from various Arthropod Management Tests, estimates that the percent of lettuce aphids killed is 76.1% for endosulfan and 67.0% for imidacloprid. [14]
Insecticides used on gypsy moths in the 1970s had roughly a 90% kill rate. [15]
Temperature change is argued to be the biggest direct abiotic impact of climate change on herbivorous insects. [16] In temperate regions, global warming will affect overwintering, and warmer temperatures will extend the summer season, allowing for more growth and reproduction. [16]
A 2013 study estimated that on average, crop pests and pathogens have moved to higher latitudes at a rate of about 2.7 km/year since 1960. [17] This is roughly in line with estimates of the rate of climate change in general. [18]
In a chemical reaction, chemical equilibrium is the state in which both the reactants and products are present in concentrations which have no further tendency to change with time, so that there is no observable change in the properties of the system. This state results when the forward reaction proceeds at the same rate as the reverse reaction. The reaction rates of the forward and backward reactions are generally not zero, but they are equal. Thus, there are no net changes in the concentrations of the reactants and products. Such a state is known as dynamic equilibrium.
In mathematical physics and mathematics, the Pauli matrices are a set of three 2 × 2 complex matrices that are Hermitian, involutory and unitary. Usually indicated by the Greek letter sigma, they are occasionally denoted by tau when used in connection with isospin symmetries.
The Pareto distribution, named after the Italian civil engineer, economist, and sociologist Vilfredo Pareto, is a power-law probability distribution that is used in description of social, quality control, scientific, geophysical, actuarial, and many other types of observable phenomena; the principle originally applied to describing the distribution of wealth in a society, fitting the trend that a large portion of wealth is held by a small fraction of the population. The Pareto principle or "80-20 rule" stating that 80% of outcomes are due to 20% of causes was named in honour of Pareto, but the concepts are distinct, and only Pareto distributions with shape value of log45 ≈ 1.16 precisely reflect it. Empirical observation has shown that this 80-20 distribution fits a wide range of cases, including natural phenomena and human activities.
Noether's theorem states that every continuous symmetry of the action of a physical system with conservative forces has a corresponding conservation law. This is the first of two theorems proven by mathematician Emmy Noether in 1915 and published in 1918. The action of a physical system is the integral over time of a Lagrangian function, from which the system's behavior can be determined by the principle of least action. This theorem only applies to continuous and smooth symmetries of physical space.
Flight dynamics is the science of air vehicle orientation and control in three dimensions. The three critical flight dynamics parameters are the angles of rotation in three dimensions about the vehicle's center of gravity (cg), known as pitch, roll and yaw. These are collectively known as aircraft attitude, often principally relative to the atmospheric frame in normal flight, but also relative to terrain during takeoff or landing, or when operating at low elevation. The concept of attitude is not specific to fixed-wing aircraft, but also extends to rotary aircraft such as helicopters, and dirigibles, where the flight dynamics involved in establishing and controlling attitude are entirely different.
The Lotka–Volterra equations, also known as the Lotka–Volterra predator–prey model, are a pair of first-order nonlinear differential equations, frequently used to describe the dynamics of biological systems in which two species interact, one as a predator and the other as prey. The populations change through time according to the pair of equations:
Population dynamics is the type of mathematics used to model and study the size and age composition of populations as dynamical systems.
The equilibrium constant of a chemical reaction is the value of its reaction quotient at chemical equilibrium, a state approached by a dynamic chemical system after sufficient time has elapsed at which its composition has no measurable tendency towards further change. For a given set of reaction conditions, the equilibrium constant is independent of the initial analytical concentrations of the reactant and product species in the mixture. Thus, given the initial composition of a system, known equilibrium constant values can be used to determine the composition of the system at equilibrium. However, reaction parameters like temperature, solvent, and ionic strength may all influence the value of the equilibrium constant.
The Clausius–Clapeyron relation, in chemical thermodynamics, specifies the temperature dependence of pressure, most importantly vapor pressure, at a discontinuous phase transition between two phases of matter of a single constituent. It is named after Rudolf Clausius and Benoît Paul Émile Clapeyron. However, this relation was in fact originally derived by Sadi Carnot in his Reflections on the Motive Power of Fire, which was published in 1824 but largely ignored until it was rediscovered by Clausius, Clapeyron, and Lord Kelvin decades later. Kelvin said of Carnot's argument that "nothing in the whole range of Natural Philosophy is more remarkable than the establishment of general laws by such a process of reasoning."
In general relativity, the Gibbons–Hawking–York boundary term is a term that needs to be added to the Einstein–Hilbert action when the underlying spacetime manifold has a boundary.
In calculus, the Leibniz integral rule for differentiation under the integral sign states that for an integral of the form
The Newman–Penrose (NP) formalism is a set of notation developed by Ezra T. Newman and Roger Penrose for general relativity (GR). Their notation is an effort to treat general relativity in terms of spinor notation, which introduces complex forms of the usual variables used in GR. The NP formalism is itself a special case of the tetrad formalism, where the tensors of the theory are projected onto a complete vector basis at each point in spacetime. Usually this vector basis is chosen to reflect some symmetry of the spacetime, leading to simplified expressions for physical observables. In the case of the NP formalism, the vector basis chosen is a null tetrad: a set of four null vectors—two real, and a complex-conjugate pair. The two real members often asymptotically point radially inward and radially outward, and the formalism is well adapted to treatment of the propagation of radiation in curved spacetime. The Weyl scalars, derived from the Weyl tensor, are often used. In particular, it can be shown that one of these scalars— in the appropriate frame—encodes the outgoing gravitational radiation of an asymptotically flat system.
Expected shortfall (ES) is a risk measure—a concept used in the field of financial risk measurement to evaluate the market risk or credit risk of a portfolio. The "expected shortfall at q% level" is the expected return on the portfolio in the worst of cases. ES is an alternative to value at risk that is more sensitive to the shape of the tail of the loss distribution.
In computational chemistry, a constraint algorithm is a method for satisfying the Newtonian motion of a rigid body which consists of mass points. A restraint algorithm is used to ensure that the distance between mass points is maintained. The general steps involved are: (i) choose novel unconstrained coordinates, (ii) introduce explicit constraint forces, (iii) minimize constraint forces implicitly by the technique of Lagrange multipliers or projection methods.
A ratio distribution is a probability distribution constructed as the distribution of the ratio of random variables having two other known distributions. Given two random variables X and Y, the distribution of the random variable Z that is formed as the ratio Z = X/Y is a ratio distribution.
In financial mathematics, tail value at risk (TVaR), also known as tail conditional expectation (TCE) or conditional tail expectation (CTE), is a risk measure associated with the more general value at risk. It quantifies the expected value of the loss given that an event outside a given probability level has occurred.
The term generalized logistic distribution is used as the name for several different families of probability distributions. For example, Johnson et al. list four forms, which are listed below.
The Larson–Miller relation, also widely known as the Larson–Miller parameter and often abbreviated LMP, is a parametric relation used to extrapolate experimental data on creep and rupture life of engineering materials.
Thermo-mechanical fatigue is the overlay of a cyclical mechanical loading, that leads to fatigue of a material, with a cyclical thermal loading. Thermo-mechanical fatigue is an important point that needs to be considered, when constructing turbine engines or gas turbines.
The modified lognormal power-law (MLP) function is a three parameter function that can be used to model data that have characteristics of a log-normal distribution and a power law behavior. It has been used to model the functional form of the initial mass function (IMF). Unlike the other functional forms of the IMF, the MLP is a single function with no joining conditions.