Agent-based model in biology
Encyclopedia
Agent-based models have many applications in biology
Biology
Biology is a natural science concerned with the study of life and living organisms, including their structure, function, growth, origin, evolution, distribution, and taxonomy. Biology is a vast subject containing many subdivisions, topics, and disciplines...

, primarily due to the characteristics of the modeling method. Agent-based modeling is a rule-based, computational modeling methodology that focuses on rules and interactions among the individual components or the agents of the system. The goal of this modeling method is to generate populations of the system components of interest and simulate their interactions in a virtual world. Agent-based models start with rules for behavior and seek to reconstruct, through computational instantiation of those behavioral rules, the observed patterns of behavior. Several of the characteristics of agent-based models important to biological studies include:
  1. Modular structure: The behavior of an agent-based model is defined by the rules of its agents. Existing agent rules can be modified or new agents can be added without having to modify the entire model.
  2. Emergent properties: Through the use of the individual agents that interact locally with rules of behavior, agent-based models result in a synergy that leads to a higher level whole with much more intricate behavior than those of each individual agent.
  3. Abstraction: Either by excluding non-essential details or when details are not available, agent-based models can be constructed in the absence of complete knowledge of the system under study. This allows the model to be as simple and verifiable as possible.
  4. Stochasticity: Biological systems exhibit behavior that appears to be random. The probability of a particular behavior can be determined for a system as a whole and then be translated into rules for the individual agents.


Before the agent-based model can be developed, one must choose the appropriate software or modeling toolkit to be used. Madey and Nikolai provide an extensive list of toolkits in their paper “Tools of the Trade: A Survey of Various Agent Based Modeling Platforms”. The paper seeks to provide users with a method of choosing a suitable toolkit by examining five characteristics across the spectrum of toolkits: the programming language
Programming language
A programming language is an artificial language designed to communicate instructions to a machine, particularly a computer. Programming languages can be used to create programs that control the behavior of a machine and/or to express algorithms precisely....

 required to create the model, the required operating system
Operating system
An operating system is a set of programs that manage computer hardware resources and provide common services for application software. The operating system is the most important type of system software in a computer system...

, availability of user support, the software license type, and the intended toolkit domain. Some of the more commonly used toolkits include Swarm, NetLogo
NetLogo
NetLogo is a multi-agent programming language and integrated modeling environment.-About:NetLogo was designed in the spirit of the Logo programming language to be "low threshold and no ceiling," that is to enable easy entry by novices and yet meet the needs of high powered users. The NetLogo...

, RePast, and Mason. Listed below are summaries of several articles describing agent-based models that have been employed in biological studies. The summaries will provide a description of the problem space, an overview of the agent-based model and the agents involved, and a brief discussion of the model results.

Forest insect infestations

In the paper titled “Exploring Forest Management Practices Using an Agent-Based Model of Forest Insect Infestations”, an agent-based model was developed to simulate attack behavior of the Mountain Pine Beetle
Mountain pine beetle
The mountain pine beetle Dendroctonus ponderosae, is a species of bark beetle native to the forests of western North America from Mexico to central British Columbia. It has a hard black exoskeleton, and measures about 5 millimeters, about the size of a grain of rice.Mountain pine beetles inhabit...

, Dendroctonus ponderosae, (MPB) in order to evaluate how different harvesting policies influence spatial characteristics of the forest and spatial propagation of the MPB infestation over time. About two-thirds of the land in British Columbia
British Columbia
British Columbia is the westernmost of Canada's provinces and is known for its natural beauty, as reflected in its Latin motto, Splendor sine occasu . Its name was chosen by Queen Victoria in 1858...

, Canada
Canada
Canada is a North American country consisting of ten provinces and three territories. Located in the northern part of the continent, it extends from the Atlantic Ocean in the east to the Pacific Ocean in the west, and northward into the Arctic Ocean...

 is covered by forests that are constantly being modified by natural disturbances such as fire, disease, and insect infestation. Forest resources make up approximately 15% of the province’s economy, so infestations caused by insects such as the MPB can have significant impacts on the economy. The MPB outbreaks are considered a major natural disturbance that can result in widespread mortality of the Lodgepole pine
Lodgepole Pine
Lodgepole Pine, Pinus contorta, also known as Shore Pine, is a common tree in western North America. Like all pines, it is evergreen.-Subspecies:...

 tree, one of the most abundant commercial tree species in British Columbia. Insect outbreaks have resulted in the death of trees over areas of several thousand square kilometers.

The agent-based model developed for this study was designed to simulate the MPB attack behavior in order to evaluate how management
Forest management
200px|thumb|right|[[Sustainable development|Sustainable]] forest management carried out by [[Complejo Forestal y Maderero Panguipulli|Complejo Panguipulli]] has contributed to the preservation of the forested landscape around [[Neltume]], a sawmill town in Chile...

 practices influence the spatial distribution and patterns of insect population and their preferences for attacked and killed trees. Three management strategies were considered by the model: 1) no management, 2) sanitation harvest and 3) salvage harvest. In the model, the Beetle Agent represented the MPB behavior; the Pine Agent represented the forest environment and tree health evolution; the Forest Management Agent represented the different management strategies. The Beetle Agent follows a series of rules to decide where to fly within the forest and to select a healthy tree to attack, feed, and breed. The MPB typically kills host trees in its natural environment in order to successfully reproduce. The beetle larvae
Larvae
In Roman mythology, lemures were shades or spirits of the restless or malignant dead, and are probably cognate with an extended sense of larvae as disturbing or frightening...

 feed on the inner bark of mature host trees, eventually killing them. In order for the beetles to reproduce, the host tree must be sufficiently large and have thick inner bark. The MPB outbreaks end when the food supply decreases to the point that there is not enough to sustain the population or when climatic conditions become unfavorable for the beetle. The Pine Agent simulates the resistance of the host tree, specifically the Lodgepole pine tree, and monitors the state and attributes of each stand of trees. At some point in the MPB attack, the number of beetles per tree reaches the host tree capacity. When this point is reached, the beetles release a chemical to direct beetles to attack other trees. The Pine Agent models this behavior by calculating the beetle population density per stand and passes the information to the Beetle Agents. The Forest Management Agent was used, at the stand level, to simulate two common silviculture
Silviculture
Silviculture is the practice of controlling the establishment, growth, composition, health, and quality of forests to meet diverse needs and values. The name comes from the Latin silvi- + culture...

 practices (sanitation and salvage) as well as the strategy where no management practice was employed. With the sanitation harvest strategy, if a stand has an infestation rate greater than a set threshold, the stand is removed as well as any healthy neighbor stand when the average size of the trees exceeded a set threshold. For the salvage harvest strategy, a stand is removed even it is not under a MPB attack if a predetermined number of neighboring stands are under a MPB attack.

The study considered a forested area in the North-Central Interior of British Columbia of approximately 560 hectare
Hectare
The hectare is a metric unit of area defined as 10,000 square metres , and primarily used in the measurement of land. In 1795, when the metric system was introduced, the are was defined as being 100 square metres and the hectare was thus 100 ares or 1/100 km2...

. The area consisted primarily of Lodgepole pine with smaller proportions of Douglas fir and White spruce
White Spruce
Picea glauca is a species of spruce native to boreal forests in the north of North America, from central Alaska east to Newfoundland, and south to northern Montana, Minnesota, Wisconsin, Michigan, upstate New York, Vermont, New Hampshire, and Maine; there is also an isolated population in the...

. The model was executed for five time steps, each step representing a single year. Thirty simulation runs were conducted for each forest management strategy considered. The results of the simulation showed that when no management strategy was employed, the highest overall MPB infestation occurred. The results also showed that the salvage forest management technique resulted in a 25% reduction in the number of forest strands killed by the MPB, as opposed to a 19% reduction by the salvage forest management strategy. In summary, the results show that the model can be used as a tool to build forest management policies.

Invasive species

Invasive species
Invasive species
"Invasive species", or invasive exotics, is a nomenclature term and categorization phrase used for flora and fauna, and for specific restoration-preservation processes in native habitats, with several definitions....

 refers to “non-native” plants and animals that adversely affect the environments they invade. The introduction of invasive species may have environmental, economic, and ecological implications. In the paper titled “An Agent-Based Model of Border Enforcement for Invasive Species Management”, an agent-based model is presented that was developed to evaluate the impacts of port
Port of entry
In general, a port of entry is a place where one may lawfully enter a country. It typically has a staff of people who check passports and visas and inspect luggage to assure that contraband is not imported. International airports are usually ports of entry, as are road and rail crossings on a...

-specific and importer-specific enforcement regimes for a given agricultural commodity
Crop
Crop may refer to:* Crop, a plant grown and harvested for agricultural use* Crop , part of the alimentary tract of some animals* Crop , a modified whip used in horseback riding or disciplining humans...

 that presents invasive species risk. Ultimately, the goal of the study was to improve the allocation of enforcement resources and to provide a tool to policy makers to answer further questions concerning border enforcement and invasive species risk.

The agent-based model developed for the study considered three types of agents: invasive species, importers, and border enforcement agents. In the model, the invasive species can only react to their surroundings, while the importers and border enforcement agents are able to make their own decisions based on their own goals and objectives. The invasive species has the ability to determine if it has been released in an area containing the target crop, and to spread to adjacent plots of the target crop. The model incorporates spatial probability maps that are used to determine if an invasive species becomes established. The study focused on shipments of broccoli
Broccoli
Broccoli is a plant in the cabbage family, whose large flower head is used as a vegetable.-General:The word broccoli, from the Italian plural of , refers to "the flowering top of a cabbage"....

 from Mexico
Mexico
The United Mexican States , commonly known as Mexico , is a federal constitutional republic in North America. It is bordered on the north by the United States; on the south and west by the Pacific Ocean; on the southeast by Guatemala, Belize, and the Caribbean Sea; and on the east by the Gulf of...

 into California
California
California is a state located on the West Coast of the United States. It is by far the most populous U.S. state, and the third-largest by land area...

 through the ports of entry Calexico, California
Calexico, California
Calexico is a city in Imperial County, California. The population was 38,572 at the 2010 census, up from 27,109 at the 2000 census. Calexico is about east of San Diego and west of Yuma, Arizona...

 and Otay Mesa, California. The selected invasive species of concern was the crucifer flea beetle
Beetle
Coleoptera is an order of insects commonly called beetles. The word "coleoptera" is from the Greek , koleos, "sheath"; and , pteron, "wing", thus "sheathed wing". Coleoptera contains more species than any other order, constituting almost 25% of all known life-forms...

 (Phyllotreta cruciferae). California is by far the largest producer of broccoli in the United States and so the concern and potential impact of an invasive species introduction through the chosen ports of entry is significant. The model also incorporated a spatially explicit damage function that was used to model invasive species damage in a realistic manner. Agent-based modeling provides the ability to analyze the behavior of heterogeneous actors, so three different types of importers were considered that differed in terms of commodity infection rates (high, medium, and low), pretreatment choice, and cost of transportation to the ports. The model gave predictions on inspection rates for each port of entry and importer and determined the success rate of border agent inspection, not only for each port and importer but also for each potential level of pretreatment (no pretreatment, level one, level two, and level three).

The model was implemented and ran in NetLogo
NetLogo
NetLogo is a multi-agent programming language and integrated modeling environment.-About:NetLogo was designed in the spirit of the Logo programming language to be "low threshold and no ceiling," that is to enable easy entry by novices and yet meet the needs of high powered users. The NetLogo...

, version 3.1.5. Spatial information on the location of the ports of entry, major highways, and transportation routes was included in the analysis as well as a map of California broccoli crops layered with invasive species establishment probability maps. BehaviorSpace, a software tool integrated with NetLogo, was used to test the effects of different parameters (e.g. shipment value, pretreatment cost) in the model. On average, 100 iterations were calculated at each level of the parameter being used, where an iteration represented a one-year run.

The results of the model showed that as inspection efforts increase, importers increase due care, or the pretreatment of shipments, and the total monetary loss of California crops decreases. The model showed that importers respond differently to an increase in inspection effort. Some importers responded to increased inspection rate by increasing pretreatment effort, while others chose to avoid shipping to a specific port, or shopped for another port. An important result of the model results is that it can show or provide recommendations to policy makers about the point at which importers may start to shop for ports, such as the inspection rate at which port shopping is introduced and the importers associated with a certain level of pest risk or transportation cost are likely to make these changes. Another interesting outcome of the model is that when inspectors were not able to learn to respond to an importer with previously infested shipments, damage to California broccoli crops was estimated to be $150 million. However, when inspectors were able to increase inspection rates of importers with previous violations, damage to the California broccoli crops was reduced by approximately 12%. The model provides a mechanism to predict the introduction of invasive species from agricultural imports and their likely damage. Equally as important, the model provides policy makers and border control agencies with a tool that can be used to determine the best allocation of inspectional resources.

Aphid population dynamics

In the article titled “Aphid Population Dynamics in Agricultural Landscapes: An Agent-based Simulation Model”, an agent-based model is presented to study the population dynamics
Population dynamics
Population dynamics is the branch of life sciences that studies short-term and long-term changes in the size and age composition of populations, and the biological and environmental processes influencing those changes...

 of the bird cherry-oat aphid
Rhopalosiphum
Rhopalosiphum is a genus of aphid of the family Aphididae that includes 16 species worldwide. In addition to direct behavior as pests on host plants, some species are vectors for viral pathogens.-Species:Species include:...

, Rhopalosiphum padi (L.). The study was conducted in a five square kilometer region of North Yorkshire
North Yorkshire
North Yorkshire is a non-metropolitan or shire county located in the Yorkshire and the Humber region of England, and a ceremonial county primarily in that region but partly in North East England. Created in 1974 by the Local Government Act 1972 it covers an area of , making it the largest...

, a county located in the Yorkshire
Yorkshire
Yorkshire is a historic county of northern England and the largest in the United Kingdom. Because of its great size in comparison to other English counties, functions have been increasingly undertaken over time by its subdivisions, which have also been subject to periodic reform...

 and the Humber
Humber
The Humber is a large tidal estuary on the east coast of Northern England. It is formed at Trent Falls, Faxfleet, by the confluence of the tidal River Ouse and the tidal River Trent. From here to the North Sea, it forms part of the boundary between the East Riding of Yorkshire on the north bank...

 region of England
England
England is a country that is part of the United Kingdom. It shares land borders with Scotland to the north and Wales to the west; the Irish Sea is to the north west, the Celtic Sea to the south west, with the North Sea to the east and the English Channel to the south separating it from continental...

. The agent-based modeling method was chosen because of its focus on the behavior of the individual agents rather than the population as a whole. The authors propose that traditional models that focus on populations as a whole do not take into account the complexity of the concurrent interactions in ecosystems, such as reproduction and competition for resources which may have significant impacts on population trends. The agent-based modeling approach also allows modelers to create more generic and modular models that are more flexible and easier to maintain than modeling approaches that focus on the population as a whole. Other proposed advantages of agent-based models include realistic representation of a phenomenon of interest due to the interactions of a group of autonomous agents, and the capability to integrate quantitative variables, differential equations, and rule based behavior into the same model.

The model was implemented in the modeling toolkit Repast using the JAVA
Java
Java is an island of Indonesia. With a population of 135 million , it is the world's most populous island, and one of the most densely populated regions in the world. It is home to 60% of Indonesia's population. The Indonesian capital city, Jakarta, is in west Java...

 programming language. The model was run in daily time steps and focused on the autumn and winter seasons. Input data for the model included habitat data, daily minimum, maximum, and mean temperatures, and wind speed and direction. For the Aphid agents, age, position, and morphology
Morphology (biology)
In biology, morphology is a branch of bioscience dealing with the study of the form and structure of organisms and their specific structural features....

 (alate
Alate
An alate is a winged reproductive of a social insect . Alate females are typically those destined to become queens , whereas alate males are occasionally referred to as "drones"...

 or apterous) were considered. Age ranged from 0.00 to 2.00, with 1.00 being the point at which the agent becomes an adult. Reproduction by the Aphid agents is dependent on age, morphology, and daily minimum, maximum, and mean temperatures. Once nymphs
Nymph (biology)
In biology, a nymph is the immature form of some invertebrates, particularly insects, which undergoes gradual metamorphosis before reaching its adult stage. Unlike a typical larva, a nymph's overall form already resembles that of the adult. In addition, while a nymph moults it never enters a...

 hatch, they remain in the same location as their parents. The morphology of the nymphs is related to population density and the nutrient quality of the aphid’s
Aphid
Aphids, also known as plant lice and in Britain and the Commonwealth as greenflies, blackflies or whiteflies, are small sap sucking insects, and members of the superfamily Aphidoidea. Aphids are among the most destructive insect pests on cultivated plants in temperate regions...

 food source. The model also considered mortality among the Aphid agents, which is dependent on age, temperatures, and quality of habitat. The speed at which an Aphid agent ages is determined by the daily minimum, maximum, and mean temperatures. The model considered movement of the Aphid agents to occur in two separate phases, a migratory phase and a foraging phase, both of which affect the overall population distribution.

The study started the simulation run with an initial population of 10,000 alate aphids distributed across a grid of 25 meter cells. The simulation results showed that there were two major population peaks, the first in early autumn due to an influx of alate immigrants and the second due to lower temperatures later in the year and a lack of immigrants. Ultimately, it is the goal of the researchers to adapt this model to simulate broader ecosystems and animal types.

Aquatic population dynamics

In the article titled “Exploring Multi-Agent Systems In Aquatic Population Dynamics Modeling”, a model is proposed to study the population dynamics of two species of macrophytes. Aquatic plants play a vital role in the ecosystems in which they live as they may provide shelter and food for other aquatic organisms. However, they may also have harmful impacts such as the excessive growth of non-native plants or eutrophication
Eutrophication
Eutrophication or more precisely hypertrophication, is the movement of a body of water′s trophic status in the direction of increasing plant biomass, by the addition of artificial or natural substances, such as nitrates and phosphates, through fertilizers or sewage, to an aquatic system...

 of the lakes in which they live leading to anoxic conditions. Given these possibilities, it is important to understand how the environment and other organisms affect the growth of these aquatic plants to allow mitigation or prevention of these harmful impacts.

Potamogeton pectinatus
Potamogeton
Potamogeton is a genus of aquatic, mostly freshwater, plants of the family Potamogetonaceae. Most are known by the common name pondweed, although many unrelated plants may be called pondweed, such as Canadian pondweed...

 is one of the aquatic plant agents in the model. It is an annual growth plant
Annual plant
An annual plant is a plant that usually germinates, flowers, and dies in a year or season. True annuals will only live longer than a year if they are prevented from setting seed...

 that absorbs nutrients from the soil and reproduces through root tubers and rhizomes. Reproduction of the plant is not impacted by water flow, but can be influenced by animals, other plants, and humans. The plant can grow up to two meters tall, which is a limiting condition because it can only grow in certain water depths, and most of its biomass is found at the top of the plant in order to capture the most sunlight possible. The second plant agent in the model is Chara aspera, also a rooted aquatic plant. One major difference in the two plants is that the latter reproduces through the use of very small seeds called oospores and bulbills which are spread via the flow of water. Chara aspera only grows up to 20 cm and requires very good light conditions as well as good water quality, all of which are limiting factors on the growth of the plant. Chara aspera has a higher growth rate than Potamogeton pectinatus but has a much shorter life span. The model also considered environmental and animal agents. Environmental agents considered included water flow, light penetration, and water depth. Flow conditions, although not of high importance to Potamogeton pectinatus, directly impact the seed dispersal of Chara aspera. Flow conditions affect the direction as well as the distance the seeds will be distributed. Light penetration strongly influences Chara aspera as it requires high water quality. Extinction coefficient
Extinction coefficient
Extinction coefficient refers to several different measures of the absorption of light in a medium:*In chemistry, biochemistry, molecular biology, or microbiology, the mass extinction coefficient and the molar extinction coefficient are parameters defining how strongly a substance absorbs light...

 (EC) is a measure of light penetration in water. As EC increases, the growth rate of Chara aspera decreases. Finally, depth is important to both species of plants. As water depth increases, the light penetration decreases making it difficult for either species to survive beyond certain depths.

The area of interest in the model was a lake in the Netherlands
Netherlands
The Netherlands is a constituent country of the Kingdom of the Netherlands, located mainly in North-West Europe and with several islands in the Caribbean. Mainland Netherlands borders the North Sea to the north and west, Belgium to the south, and Germany to the east, and shares maritime borders...

 named Lake Veluwe. It is a relatively shallow lake with an average depth of 1.55 meters and covers about 30 square kilometers. The lake is under eutrophication stress which means that nutrients are not a limiting factor for either of the plant agents in the model. The initial position of the plant agents in the model was randomly determined. The model was implemented using Repast software package and was executed to simulate the growth and decay of the two different plant agents, taking into account the environmental agents previously discussed as well as interactions with other plant agents. The results of the model execution show that the population distribution of Chara aspera has a spatial pattern very similar to the GIS maps of observed distributions. The authors of the study conclude that the agent rules developed in the study are reasonable to simulate the spatial pattern of macrophyte growth in this particular lake.
The source of this article is wikipedia, the free encyclopedia.  The text of this article is licensed under the GFDL.
 
x
OK