AncientS-ABM is an agent-based model for simulating and evaluating the potential social organization of an artificial past society, configured by available archaeological data. Unlike most existing agent-based models used in archaeology, our ABM framework includes completely autonomous, utility-based agents. It also incorporates different social organization paradigms, different decision-making processes, and also different cultivation technologies used in ancient societies. Equipped with such paradigms, the model allows us to explore the transition from a simple to a more complex society by focusing on the historical social dynamics; and to assess the influence of social organization on agents’ population growth, agent community numbers, sizes and distribution.

AncientS-ABM also blends ideas from evolutionary game theory with multi-agent systems’ self-organization. We model the evolution of social behaviours in a population of strategically interacting agents in repeated games where they exchange resources (utility) with others. The results of the games contribute to both the continuous re-organization of the social structure, and the progressive adoption of the most successful agent strategies. Agent population is not fixed, but fluctuates over time, while agents in stage games also receive non-static payoffs, in contrast to most games studied in the literature. To tackle this, we defined a novel formulation of the evolutionary dynamics via assessing agents’ rather than strategies’ fitness.

As a case study, we employ AncientS-ABM to evaluate the impact of the implemented social organization paradigms on an artificial Bronze Age “Minoan” society, located at different geographical parts of the island of Crete, Greece. Model parameter choices are based on archaeological evidence and studies, but are not biased towards any specific assumption. Results over a number of different simulation scenarios demonstrate better sustainability for settlements consisting of and adopting a socio-economic organization model based on self-organization, where a “heterarchical” social structure emerges. Results also demonstrate that successful agent societies adopt an evolutionary approach where cooperation is an emergent strategic behaviour. In simulation scenarios where the natural disaster module was enabled, we observe noticeable changes in the settlements’ distribution, relating to significantly higher migration rates immediately after the modeled Theran eruption. In addition, the initially cooperative behaviour is transformed to a non-cooperative one, thus providing support for archaeological theories suggesting that the volcanic eruption led to a clear breakdown of the Minoan socio-economic system.

…

This is the R code of the mathematical model that includes the decision making formulations for artificial agents. This code corresponds to equations 1-70 given in the paper “A Mathematical Model of The Beer Game”.

This is the R code of the mathematical model used for verification. This code corresponds to equations 1-9, 15-53, 58-62, 69-70, and 72-75 given in the paper “A Mathematical Model of The Beer Game”.

This is the R code of the mathematical model that includes the decision making formulations for artificial agents. Plus, the code for graphical output is also added to the original code.

The modeling includes citizens, bounded into families; firms and governments; all of them interacting in markets for goods, labor and real estate. The model is spatial and dynamic.

The Nice Musical Chairs (NMC) model represent the competition for space between groups of stakeholders of farming and herding activities in the arid Afro-Eurasia.

This model is based on the Narragansett Bay, RI recreational fishery. The two types of agents are piscivorous fish and fishers (shore and boat fishers are separate “breeds”). Each time step represents one week. Open season is weeks 1-26, assuming fishing occurs during half the year. At each weekly time step, fish agents grow, reproduce, and die. Fisher agents decide whether or not to fish based on their current satisfaction level, and those that do go fishing attempt to catch a fish. If they are successful, they decide whether to keep or release the fish. In our publication, this model was linked to an Ecopath with Ecosim food web model where the commercial harvest of forage fish affected the biomass of piscivorous fish - which then became the starting number of piscivorous fish for this ABM. The number of fish caught in a season of this ABM was converted to a fishing pressure and input back into the food web model.

This model extends the original Artifical Anasazi (AA) model to include individual agents, who vary in age and sex, and are aggregated into households. This allows more realistic simulations of population dynamics within the Long House Valley of Arizona from AD 800 to 1350 than are possible in the original model. The parts of this model that are directly derived from the AA model are based on Janssen’s 1999 Netlogo implementation of the model; the code for all extensions and adaptations in the model described here (the Artificial Long House Valley (ALHV) model) have been written by the authors. The AA model included only ideal and homogeneous “individuals” who do not participate in the population processes (e.g., birth and death)–these processes were assumed to act on entire households only. The ALHV model incorporates actual individual agents and all demographic processes affect these individuals. Individuals are aggregated into households that participate in annual agricultural and demographic cycles. Thus, the ALHV model is a combination of individual processes (birth and death) and household-level processes (e.g., finding suitable agriculture plots).

As is the case for the AA model, the ALHV model makes use of detailed archaeological and paleoenvironmental data from the Long House Valley and the adjacent areas in Arizona. It also uses the same methods as the original model (from Janssen’s Netlogo implementation) to estimate annual maize productivity of various agricultural zones within the valley. These estimates are used to determine suitable locations for households and farms during each year of the simulation.

This model is an extension of the Artificial Long House Valley (ALHV) model developed by the authors (Swedlund et al. 2016; Warren and Sattenspiel 2020). The ALHV model simulates the population dynamics of individuals within the Long House Valley of Arizona from AD 800 to 1350. Individuals are aggregated into households that participate in annual agricultural and demographic cycles. The present version of the model incorporates features of the ALHV model including realistic age-specific fertility and mortality and, in addition, it adds the Black Mesa environment and population, as well as additional methods to allow migration between the two regions.

As is the case for previous versions of the ALHV model as well as the Artificial Anasazi (AA) model from which the ALHV model was derived (Axtell et al. 2002; Janssen 2009), this version makes use of detailed archaeological and paleoenvironmental data from the Long House Valley and the adjacent areas in Arizona. It also uses the same methods as the original AA model to estimate annual maize productivity of various agricultural zones within the Long House Valley. A new environment and associated methods have been developed for Black Mesa. Productivity estimates from both regions are used to determine suitable locations for households and farms during each year of the simulation.

This model describes the consequences of limited vision of agents in harvesting a common resource. We show the vulnerability of cooperation due to reduced visibility of the resource and other agents.