Computational Model Library

Our mission is to help computational modelers develop, document, and share their computational models in accordance with community standards and good open science and software engineering practices. Model authors can publish their model source code in the Computational Model Library with narrative documentation as well as metadata that supports open science and emerging norms that facilitate software citation, computational reproducibility / frictionless reuse, and interoperability. Model authors can also request private peer review of their computational models. Models that pass peer review receive a DOI once published.

All users of models published in the library must cite model authors when they use and benefit from their code.

Please check out our model publishing tutorial and feel free to contact us if you have any questions or concerns about publishing your model(s) in the Computational Model Library.

Displaying 10 of 247 results NetLogo clear search

Netlogo Aquaponic Simulation

Derrick | Published Tuesday, October 06, 2026 | Last modified Wednesday, October 07, 2026

This agent-based aquaponics simulation models the biological resilience and complex dynamics of a closed-loop agricultural system. Designed for educational use and agricultural research, the model helps aquaponic designers, farmers, and hobbyists understand systemic dynamics and stress-test species configurations prior to a physical build. The simulation dynamically models the nitrogen cycle by tracking aquatic bio-loads, bacterial filtration, and plant nutrient uptake. Utilizing real-world biological metrics from the Loyola Global Aquatic Excrement Dataset, the model calculates species-specific ammonia excretion rates (toxicity) for various aquatic life, including Trout, Tilapia, and Goldfish. Nitrifying bacteria agents convert this toxicity into usable nitrates, which are then absorbed by specific crop agents—such as Tomatoes, Swiss Chard, or Lettuce—at variable rates. Key features include dynamic stochastic harvesting, measurable agent health, vermaculture and crawdad integration, and reactive biological mitigators like emergency duckweed buffering. The system also introduces chaotic environmental variables, allowing users to test the agricultural setup’s resilience against random shocks like bio-filter pump failures and temperature crashes. This model is released under the GNU General Public License v3.0 (GPL-3.0) to encourage collaborative improvement within the academic and farming communities.

Artificial Anasazi is an agent-based model of farming households in the Long House Valley, Arizona, from 800 to 1350 CE. Each household chooses where to farm and settle based on expected harvest, water sources and stored maize, and simulated household counts are compared with the archaeological record.

This release adapts the NetLogo Models Library version (Stonedahl & Wilensky 2010) of Janssen’s (2009) NetLogo replication of the original Ascape model (Dean et al. 2000; Axtell et al. 2002). It contains four configurations in one NetLogo 7.0.4 file, selected with two switches:

  • Both switches off: the baseline, which is the Models Library model with small fixes that make it work as intended (the arable-uplands zone name, the setup base yield, the initial household count, one-tick-lag capacity reporting, and seed control with observation-only instruments).
  • soil-on?: dynamic soil quality that depletes under cultivation and regrows, following Janssen (2010, eq. 4).
    …

Aquaponics & Multitrophic System Simulation

derrickatasu | Published Friday, September 25, 2026 | Last modified Thursday, October 08, 2026

An agent-based model built in NetLogo that simulates a closed-loop recirculating aquaponics and integrated multitrophic aquaculture (IMTA) ecosystem, modeling the complex dynamics between aquatic species, vermaculture, and hydroponic plant beds. This model captures the operational mechanics of a complex, sustainable production network. It tracks nutrient flows, water recirculation pathways, biological growth stages, and environmental feedback penalties (such as nutrient toxicity or starvation) across multiple specialized tanks and grow beds. Model Architecture & HabitatsAquatic Subsystem: Duckweed Tank, Breeding Tank, Fry Tank, Fingerling Tank, Adult Tank, and a Holding Tank for harvested fish. Benthic & Waste Management: Crawdad Tank 1 & 2 for bottom-feeding and waste processing, alongside a Vermaculture node populated by composting worms. Hydroponic Subsystem: Five sequential Grow Beds and a recirculating loop feeding back into the system. Key Agent DynamicsWater Flow Visualization: Active links between tanks spawn animated droplet agents that visually simulate recirculating water dynamics across the network. Biological Growth & Aging: Fish and plants progress through age-based development thresholds governed by species-specific parameters. Feedback Loops:Toxicity Penalty: Excess nutrients (nutrient-per-day > 80) slow fish growth rates by 50%. Starvation Penalty: Insufficient nutrient levels (nutrient-per-day < 30) delay plant maturation by 50%. Interface Controlsfish-type Chooser: Selects the primary aquatic species (Tilapia, Channel Catfish, Trout, Goldfish), altering base growth thresholds. plant-type Chooser: Selects the cultivated crop (Lettuce, Basil, Tomatoes, Mint). nutrient-per-day Slider: Manages daily nutrient input to balance system load and avoid toxicity or starvation thresholds. Dashboard Plot: Real-time tracking of total-fish-harvested and total-plants-harvested.

This NetLogo model simulates the movement and foraging behavior of 20 bird species (classified as urban avoiders or utilizers) across a 20×20m grid landscape representing southern Bogotá, Colombia. Agents follow an optimal foraging strategy (90% of the time) or move randomly (10%), gaining energy from land-cover-dependent food sources and dying if energy falls below their basal metabolic rate. The model compares two landscape scenarios — current land cover (2016) and a proposed scenario incorporating the ‘Media Luna del Sur’ ecological corridor from Bogotá’s 2022–2035 Land Use Plan (POT) — to evaluate whether the proposed green infrastructure improves functional connectivity for avian biodiversity across the urban matrix. Outputs include cumulative patch visitation counts, exported as raster files, used to generate connectivity heatmaps.

HousingABM_Japan is a NetLogo agent-based model of the residential market of Tokyo’s 23 wards. It evaluates whether a single parameter configuration can jointly reproduce key features of prices, rents, yields, and market turnover across distinct market regimes from 2001 to 2025, with particular attention to demand- and supply-side trend-following during the 2021–2025 price surge.

The model builds on the Bank of England housing-market lineage (Baptista et al. 2016; Carro et al. 2023) and introduces four extensions: (1) dynamic linkages between the sale and rental markets through vacancy, rents, and yields; (2) heterogeneous demand-side trend-following; (3) supply-side trend-following through construction-cost trend anchoring and a momentum-dependent dynamic premium; and (4) housing-equity borrowing that converts unrealized equity into additional borrowing capacity.

Twenty parameters are calibrated using 2001–2015 data and held fixed for post-calibration evaluation over 2016–2020 and 2021–2025, while annual exogenous inputs follow observed historical paths. The model reproduces the shift from moderate price growth to the 2021–2025 surge, as well as rent acceleration, surge-period yield compression, and persistently low market turnover, although it understates the intermediate acceleration of 2016–2020.

…

This NetLogo model simulates the impact of mowing frequency, mower type and the percentage of refuge strips on the abundance of arthropods in managed grasslands. It represents five functional arthropod groups: holometabolous plant-dwelling arthropods (HoP, represented by butterflies); holometabolous ground-dwelling arthropods (HoG, represented by ants); nesting pollinators (NP, represented by bees); hemimetabolous plant-dwelling arthropods (HeP, represented by grasshoppers); and ground-resident arthropods (GR, represented by spiders). Three mowing frequencies are simulated: intensive (four mowings per year), intermediate (two mowings per year) and extensive (one mowing per year). Bar and disc mowers differ in their effects on arthropod mortality and vegetation height. Refuge strip coverage can be set to 0%, 10% or 20%. Arthropod abundance changes over time through movement, reproduction, natural mortality, mower-related mortality and vegetation regrowth.

SERDUX-MARCIM simulates the propagation of a cyberattack over the computational network of an organization in the maritime sector, at the strategic level of decision-making. It is the instantiation of ABM-MARCIM, the agent-based model of the MARCIM framework for the modeling and simulation of maritime cyberdefense.

Every computational asset of the target organization – servers, endpoints, routers, gateways, vessel systems, radars – is an agent that occupies one of six states at each time step: Susceptible, Exposed, Resistant, Degraded, Unavailable or Destroyed, the initials of which give the model its name. The states Degraded, Unavailable and Destroyed are associated with the D5 cyberattack effects (disrupt, degrade, deny, destroy, deceive) as a function of the degree and the duration of the attack. Thirteen transitions between states are admissible.

Unlike a conventional agent-based model, the local update function is not an individual behavioral rule. It is a system of six ordinary differential equations with eight time-dependent transition rates – propagation, cyberattack (degraded), cyberattack (unavailable), cyberattack (destroyed), recovery, sanitation, loss of resistance, and unavailability by other causes. The values of those rates derive from the capabilities of the target organization, the capabilities of the attacker, and the degree and duration of the cyberattack, computed through a cyber risk approach aligned with the OWASP Risk Rating Methodology, the ISACA categorization of security controls and the IMO Guidelines on Maritime Cyber Risk Management.

…

Peer reviewed LSCs-CF: An agent-based model of public space conditions and the collective formation of local sports communities

Yijia JIANG | Published Thursday, August 13, 2026 | Last modified Thursday, August 20, 2026

An agent-based model of how public space conditions shape the collective formation of local sports communities in high-mobility cities. A small gathering is already exercising in an abstract, systematically configurable public space; the model asks how strangers join them and whether the loose gathering consolidates into a stable community. Three conditions are manipulated, each acting through exactly one channel: spatial legibility governs whether a passer-by notices the activity, routine compatibility governs whether a member attends in a given week, and social anchor availability governs how much one episode of co-presence builds a tie. Ties carry a continuous strength with two thresholds, so that relationships of differing strength cross their thresholds at different condition levels. Because the model does not represent the origin of an activity, dissolution is irreversible and survival is a clean binary outcome. Across 32,790 runs in eight experiments the model yields two well separated thresholds along routine compatibility — one for survival, one for consolidation — with a wide band between them in which a space is used continuously but no community forms. The tie submodel is adapted from Jin, Girvan and Newman (2001) and verified against it by docking. Robustness is established against the seeded gathering, the pre-existing social network, a joint ±20 per cent perturbation of the behavioural parameters, and the parameters of the tie machinery.

This is a NetLogo 7.0.1 agent-based model of interdependent team productivity under O-ring production logic. The model asks when collective problem-solving capacity (CPS) improves team output and when its effect is constrained by trust, burnout, specialization diversity, weak-link quality, and team-formation rules. Agents are heterogeneous in skill, CPS, trust, burnout, effort, learning rate, openness, specialization, adaptability, and aspiration. Each tick forms temporary teams, computes individual contribution quality, combines contributions through a geometric O-ring production function with an explicit weakest-link term, records output and success, and updates agents through feedback, learning, trust change, burnout, recovery, and turnover. The paper follows the ODD protocol, reports exact implementation equations, and analyzes six BehaviorSpace experiments with 50 stochastic repetitions per condition. Scenario results show that CPS is beneficial but not sufficient: productivity is highest when CPS is combined with trust, low burnout, effective coordination, and a reliable weakest-link floor. Diversity yields only modest gains at low CPS but larger gains when CPS is high. Formation-mode results are especially informative because mixed CPS-skill formation maximizes total output, whereas random formation has the highest binary success rate, showing that output magnitude and threshold success can diverge. The study contributes a transparent computational mechanism linking collective intelligence with O-ring production and identifies boundary conditions under which high-CPS teams can still be limited by weak links and interdependence.

This model tests whether local housing supply elasticity governs crash severity inside a single metropolitan housing market. Saiz (2010) established that across US metros, regions constrained by geography and regulation experience deeper boom-bust cycles than flexible ones. That finding is routinely applied downward to neighborhoods and ZIP codes as though the mechanism scaled without qualification.

The empirical record for the Washington DC and Northern Virginia region says it does not. Across 84 ZIP codes, measured supply elasticity ranges from 0.35 to 4.95 with a median of 1.21. The worst single-year price decline between 2007 and 2012 averaged 9.6 percent in constrained ZIP codes and 9.0 percent in flexible ones, a gap that cannot be distinguished from noise. Wide variation in the proposed cause, no meaningful separation in the proposed effect.

The model embeds households, houses and a metro-wide credit condition in the real ZIP geography of the region using three GIS layers and an empirical price panel. Local elasticity governs construction, exactly as theory predicts. Prices are driven by a shared macro drift schedule and, under the credit-amplified mode, by a leverage cycle with a financial accelerator and a deviation penalty. The design question is whether those shared forces are sufficient to override local supply differences at the sub-metropolitan scale.

…

Displaying 10 of 247 results NetLogo clear search

This website uses cookies and Google Analytics to help us track user engagement and improve our site. If you'd like to know more information about what data we collect and why, please see our data privacy policy. If you continue to use this site, you consent to our use of cookies.
Accept